zs19811011 发表于 2006-11-19 22:25:34

帮忙纠正我做的三维电磁场瞬态分析单相变压器模型（命令流）

我做的是单相变压器模型，并且只有原边，在原边加载交流电流，但是计算结果有错误，请大家帮忙指出命令流中的错误，谢谢了！
以下是我的命令流（由于对命令流不是很了解，所以可能贴了好多多余的命令，希望大家多多指教）：

/PREP7
ET,1,SOLID97
ET,2,SOLID97
KEYOPT,2,1,2
KEYOPT,2,2,0
KEYOPT,2,5,0
R,1,0.9,100,9.3258,0,1,0,   
RMORE,1,0.9,
UIMP,1,MURX, , ,1,
UIMP,1,MGXX, , , ,
UIMP,1,RSVX, , , ,
UIMP,1,PERX, , , ,
UIMP,2,MURX, , ,2000,   
UIMP,2,MGXX, , , ,
UIMP,2,RSVX, , , ,
UIMP,2,PERX, , , ,
UIMP,3,MURX, , ,1,
UIMP,3,MGXX, , , ,
UIMP,3,RSVX, , ,2e-8,   
UIMP,3,PERX, , , ,
BLOCK,-1,12,-1,1,-5,5,
BLOCK,1,10,-1,1,-3,3,   
/VIEW, 1 ,,,1   
/ANG, 1
/REP,FAST   
/VIEW, 1 ,,1
/ANG, 1
/REP,FAST   
FLST,2,2,6,ORDE,2   
FITEM,2,1   
FITEM,2,-2
VOVLAP,P51X
VDELE,       2, , ,1
BLOCK,-1,1,-1,1,-3,3,   
BLOCK,10,12,-1,1,-3,3,
FLST,2,3,6,ORDE,2   
FITEM,2,1   
FITEM,2,-3
VOVLAP,P51X
CYLIND,1.8,1.5,-1.5,1.5,0,90,   
CYLIND,1.8,1.5,-1.5,1.5,90,180,
CYLIND,1.8,1.5,-1.5,1.5,180,270,
CYLIND,1.8,1.5,-1.5,1.5,270,360,
/USER, 1
/VIEW, 1,0.241807359834    ,0.946033792634    , -0.215752784279
/ANG,1,   131.799871147   
/REPLO
FLST,2,4,6,ORDE,3   
FITEM,2,3   
FITEM,2,6   
FITEM,2,-8
VGLUE,P51X
CYL4, , ,30, , , ,30
/VIEW, 1 ,,,1   
/ANG, 1
/REP,FAST   
/VIEW, 1 ,,1
/ANG, 1
/REP,FAST   
/DIST, 1 ,1.371742,1
/REP,FAST   
/DIST, 1 ,1.371742,1
/REP,FAST   
/DIST, 1 ,1.371742,1
/REP,FAST   
/DIST, 1 ,1.371742,1
/REP,FAST   
/DIST, 1 ,1.371742,1
/REP,FAST   
/DIST, 1 ,1.371742,1
/REP,FAST   
/FOC, 1 ,,-0.300000,,1
/REP,FAST   
FLST,3,1,6,ORDE,1   
FITEM,3,6   
VGEN, ,P51X, , , ,15, , , ,1
/VIEW, 1,0.662558747463E-01,0.997798598758    ,0.284875033368E-02
/ANG,1,   80.9266354507   
/LIG,1,1,1.000,0.277588698958    ,0.930372259688    , -0.239482718819    ,   0.00000000000
/REPLO
FLST,3,1,6,ORDE,1   
FITEM,3,6   
VGEN, ,P51X, , , ,-15,-15, , ,1
/VIEW, 1,0.354070399814E-01,0.529300765168    , -0.847695134770
/ANG,1,   165.896287638   
/LIG,1,1,1.000,0.221205580014    ,0.273209859556    , -0.936175445102    ,   0.00000000000
/REPLO
/VIEW, 1, -0.242800348496E-02,0.345442848401E-01, -0.999400218723
/ANG,1,   166.726655353   
/LIG,1,1,1.000,0.183491678896    , -0.232668831754    , -0.955089534288    ,   0.00000000000
/REPLO
LPLOT   
/VIEW, 1 ,,,1   
/ANG, 1
/REP,FAST   
/VIEW, 1 ,,1
/ANG, 1
/REP,FAST   
/VIEW, 1,0.133874482163    ,0.943559513762    , -0.302940698841
/ANG,1,   156.277825985   
/LIG,1,1,1.000,0.132698591093    , -0.533586981185    , -0.835270026657    ,   0.00000000000
/REPLO
FLST,3,1,6,ORDE,1   
FITEM,3,6   
VGEN, ,P51X, , ,5, , , , ,1
LPLOT   
/VIEW, 1,0.154361698880    ,0.406014640281    , -0.900735575958
/ANG,1,   166.186608393   
/LIG,1,1,1.000, -0.184412486294E-01, -0.982084960857    , -0.187534130248    ,   0.00000000000
/REPLO
/VIEW, 1, -0.894956635730    , -0.245754185573E-01, -0.445475778202
/ANG,1,   156.037821842   
/LIG,1,1,1.000,0.130938757934    , -0.941292032609    ,0.311166114829    ,   0.00000000000
/REPLO
/VIEW, 1, -0.751660531290    , -0.637571083067    ,0.168847741287
/ANG,1,   164.116403760   
/LIG,1,1,1.000,0.796184476405    , -0.598045251320    , -0.918267766223E-01,   0.00000000000
/REPLO
/VIEW, 1, -0.253798776386    , -0.808733032040    ,0.530600663393
/ANG,1,   123.591326251   
/LIG,1,1,1.000,0.922925982766    , -0.212538168418    , -0.320990899095    ,   0.00000000000
/REPLO
/VIEW, 1 ,1,1,1
/ANG, 1
/REP,FAST   
/FOC, 1 ,0.300000,,,1   
/REP,FAST   
VPLOT   
FLST,2,9,6,ORDE,4   
FITEM,2,1   
FITEM,2,-6
FITEM,2,9   
FITEM,2,-11
VOVLAP,P51X
LOCAL,11,1,0,0,0, , , ,1,1,
CSYS,11,
FLST,5,4,6,ORDE,3   
FITEM,5,3   
FITEM,5,9   
FITEM,5,-11
CM,_Y,VOLU
VSEL, , , ,P51X
CM,_Y1,VOLU
CMSEL,S,_Y
CMSEL,S,_Y1
VATT,       3,       1,   2,      11
CMSEL,S,_Y
CMDELE,_Y   
CMDELE,_Y1
FLST,5,4,6,ORDE,4   
FITEM,5,1   
FITEM,5,-2
FITEM,5,4   
FITEM,5,-5
CM,_Y,VOLU
VSEL, , , ,P51X
CM,_Y1,VOLU
CMSEL,S,_Y
CMSEL,S,_Y1
VATT,       2,       1,   1,       0
CMSEL,S,_Y
CMDELE,_Y   
CMDELE,_Y1
CM,_Y,VOLU
VSEL, , , ,       7
CM,_Y1,VOLU
CMSEL,S,_Y
CMSEL,S,_Y1
VATT,       1,       1,   1,       0
CMSEL,S,_Y
CMDELE,_Y   
CMDELE,_Y1
SMRT,6
MSHAPE,1,3D
MSHKEY,0
FLST,5,9,6,ORDE,5   
FITEM,5,1   
FITEM,5,-5
FITEM,5,7   
FITEM,5,9   
FITEM,5,-11
CM,_Y,VOLU
VSEL, , , ,P51X
CM,_Y1,VOLU
CHKMSH,'VOLU'   
CMSEL,S,_Y
VMESH,_Y1   
CMDELE,_Y   
CMDELE,_Y1
CMDELE,_Y2
/UI,MESH,OFF
FLST,5,4,6,ORDE,3   
FITEM,5,3   
FITEM,5,9   
FITEM,5,-11
VSEL,S, , ,P51X
ALLSEL,BELOW,VOLU   
VPLOT   
FLST,4,470,1,ORDE,22
FITEM,4,9   
FITEM,4,-16
FITEM,4,20
FITEM,4,-21
FITEM,4,39
FITEM,4,-40
FITEM,4,45
FITEM,4,-48
FITEM,4,71
FITEM,4,-118
FITEM,4,235
FITEM,4,-240
FITEM,4,244
FITEM,4,-249
FITEM,4,255
FITEM,4,-338
FITEM,4,817
FITEM,4,-893
FITEM,4,978
FITEM,4,-1206   
FITEM,4,2180
FITEM,4,-2183   
CP,1,CURR,P51X
ALLSEL,BELOW,VOLU   
VPLOT   
ALLSEL,ALL
VPLOT   
FINISH
/SOLU   
ANTYPE,4
TRNOPT,FULL
FLST,2,4,5,ORDE,4   
FITEM,2,13
FITEM,2,-14
FITEM,2,33
FITEM,2,-34
DA,P51X,AX,0,   
FLST,2,4,5,ORDE,4   
FITEM,2,13
FITEM,2,-14
FITEM,2,33
FITEM,2,-34
DA,P51X,AY,0,   
FLST,2,4,5,ORDE,4   
FITEM,2,13
FITEM,2,-14
FITEM,2,33
FITEM,2,-34
DA,P51X,AZ,0,   
OUTRES,ALL,ALL,

zs19811011 发表于 2006-11-19 22:27:23

这是我所加的交流电流的命令流

eqslv,jcg,1e-9
*SET,t,3
*SET,c,0
*SET,n,10   
*SET,pi,2*asin(1)   
*SET,con,2*pi/n
neqit,1
*create,load
time,t
i=4*cos(t)
f,21,csgz,i
t=t+con
c=c+1   
outres,all,1
*if,c,eq,(N*.75),then   
outpr, ,1   
*elseif,c,eq,N,then
outpr, ,1   
*else   
outpr, ,0   
*endif
solve   
*end
*do,i,1,10
*use,load   
*enddo
finish

zs19811011 发表于 2006-11-20 10:14:42

还是没有人回答我啊？？？
我觉得好象是我加载电流的命令有问题，可就是不知道具体错误在哪？

zs19811011 发表于 2006-11-21 11:53:08

加载电流的使用AMPS还是csgz呢？
别的地方有错误吗

inkey 发表于 2006-11-22 20:28:21

我运行过了没有错误，但仅仅只有一个柱体
还有网格不是很好，需进一步进行细化吧！

[ 本帖最后由 inkey 于 2006-11-22 20:29 编辑 ]

zs19811011 发表于 2006-11-22 22:24:11

还是有问题，铁心内部的磁场分布有错误...

inkey 发表于 2006-11-25 12:41:41

加载电流csgz
送给你个例子看看应有所用的!
FINISH
/CLEAR

/PREP7
et,1,96
mp,murx,1,1                  !空气
mp,murx,2,1000           !衔铁
mp,murx,3,10000         !永磁体
mp,mgxx,3,10000

rd1=4
x0=3
x1=6.5
x2=8.5
x3=10

y0=0
y1=1.5
y2=7.5

dt=0.001

z0=0
z1=2
z2=z1+3
z3=z2+0.5
z4=-2
z5=15
z6=17
dz=z3+2

CYL4,0,0,rd1
rect,-x0,-x1,-y1,y1
rect,x0,x1,-y1,y1   
APTN,ALL

ASEL,S,LOC,X,-rd1,rd1
AADD,ALL
ALLS

lsel,s,loc,x,x1
lsel,a,loc,x,-x1
lesi,all,,,8,1

lsel,s,loc,y,dt-y1,-dt+y1
lesi,all,,,4

lsel,s,loc,y,y1+dt,rd1-dt
lsel,a,loc,y,-dt-y1,dt-rd1

lesi,all,,,8

lsel,s,loc,x,dt-x1,dt-rd1
lsel,a,loc,x,-dt+x1,dt+rd1
lesi,all,,,10


VEXT,ALL,,,0,0,2

ASEL,S,LOC,X,-4,4
ASEL,R,LOC,z,2
vext,all,,,0,0,3

ASEL,S,LOC,X,-4,4
ASEL,R,LOC,z,5
vext,all,,,0,0,0.5

lsel,s,loc,z,z0+dt,z1-dt
lesi,all,,,5

lsel,s,loc,z,z1+dt,z2-dt
lesi,all,,,6

lsel,s,loc,z,z2+dt,z3-dt
lesi,all,,,3

wpoff,,,dz
csys,4
alls

VSYMM,Z,all, , , ,0,0

CSYS,0
WPCSYS,-1,0

BLOCK,x1,x2,-y1,y1,z0,z5,
BLOCK,-x1,-x2,-y1,y1,z0,z5,
block,-x3,x3,-y2,y2,z4,z6

vptn,all

vsel,s,loc,z,z0,z1
vsel,a,loc,z,z5-z1,z5
vsel,a,loc,x,x1,x2
vsel,a,loc,x,-x1,-x2
vsel,a,loc,z,z5-z3,z5-z2
vsel,a,loc,z,z2,z3
cm,fe,volume                        !定义衔铁组件
vatt,2,1,1

vsel,s,loc,z,z1,z2
vsel,a,loc,z,z5-z1,z5-z2
cm,forever,volume                !定义永磁体
vatt,3,1,1

alls

cmsel,s,fe
cmsel,a,forever
vsweep,all

vsel,inve
mshkey,0
mshape,1,3d
vmesh,all
alls

WPCSYS,-1,0
CSYS,0
wpro,,,90.000000
CSWPLA,11,0,1,1,
CMSEL,S,FOREVER
ALLSEL,BELOW,VOLU   
EMODIF,all,ESYS,11,                !定义永磁体磁极方向

alls

/SOLU   
ANTYPE,0
EQSLV,ICCG, ,0,                     !定义求解器

MAGSOLV,2, , ,0.001,25,0         !进行静磁求解

zlm0622 发表于 2011-12-31 17:15:41

对这方面不了解
页: [1]
查看完整版本: 帮忙纠正我做的三维电磁场瞬态分析单相变压器模型（命令流）