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

zs19811011

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

/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

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

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

zs19811011

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

inkey

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

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

zs19811011

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

inkey

加载电流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

对这方面不了解