双凸极电机的命令流

haohx

[在论坛里面下载过双凸极电机的ppt，跟着里面操作，但是没有计算结果，也没有提示错误，不知道问题出在哪里，请高手赐教。
这是我由操作得到的命令流
/COM,   
/COM,Preferences for GUI filtering have been set to display:
KEYW,MAGNOD,1
/PREP7
! define material prop
MP,MURX,1,,1
MP,MURX,3,,1
MP,MURX,4,,1
tb,bh,2
tbpt,,290,1.2   
tbpt,,300,1.25  
tbpt,,354,1.3   
tbpt,,450,1.35  
tbpt,,620,1.4   
tbpt,,940,1.45  
tbpt,,1420,1.5  
tbpt,,2220,1.55
tbpt,,3300,1.6  
tbpt,,4830,1.65
tbpt,,6600,1.7  
tbpt,,8800,1.75
tbpt,,11500,1.8
tbpt,,15100,1.85
tbpt,,19300,1.9
tbplot
! moding
*afun, deg
d1=77.5/1000
d2=62/1000
d3=47.5/1000
d4=47.2/1000
d5=35/1000
d6=19/1000
turn=50
an1=30*0.4463
an2=45*0.3273
a1=d3/2
b1=sin(an1)*d3/2
a2=d4/2
b2=sin(an2)*d4/2
! define  stator
CYL4,0,0,d1/2, ,d2/2
CYL4,0,0,d2/2
wpoffs,0,b1,0
wprota,0,90,0
ASBW, 2
WPCSYS,-1,0
/REPLOT
wpoffs,0,-b1,0
wprota,0,90,0
ASBW,  4
ADELE,2,3 , ,1
WPCSYS,-1,0
CYL4,0,0,d3/2
ASBA, 5, 2
CSYS,1
AGEN,3,3,4 ,1 ,0,60,0, ,0
CSYS,0
aglue,all
numcmp,AREAR
CSYS,1
AGEN, ,all, , ,0,45,0, , ,1
/REPLOT

!define rotor
CSYS,0
RECTNG,-a2,a2,-b2,b2,
CYL4,0,0,d4/2           
AINA,8,9
CSYS,1
AGEN,2,10, , ,0,90,0, ,0
AADD,10,8
/REPLOT
CSYS,0
CYL4,0,0,d5/2
AADD,9,8
/REPLOT
K,100,0,0,0,0
CIRCLE,100,d6/2,
FLST,3,4,4,ORDE,4   
FITEM,3,5   
FITEM,3,8   
FITEM,3,17  
FITEM,3,-18
ASBL,      10,all
!define coil
CYL4,0,0,d3/2, ,d4/2
CYL4,0,0,d1/2
AOVLAP,ALL
wprot,0,90,0
CSYS,1
AGEN, ,all, , ,0,-15,0, , ,1
ASBW,18
ASBW,19
AGEN, ,all, , ,0,60,0, , ,1
ASBW,16
ASBW,21
AGEN, ,all, , ,0,60,0, , ,1
ASBW,17
ASBW,20
CSYS,0

!  予 牧???APLOT
asel,s,,,19
asel,a,,,21
asel,a,,,16
asel,a,,,17
asel,a,,,18
asel,a,,,10
AATT, 3, , , 0,
ALLSEL,ALL
asel,s,,,11
asel,a,,,7
asel,a,,,24,27
AATT, 4, , , 0,
ALLSEL,ALL
asel,s,,,1,6
asel,a,,,9
asel,a,,,22
AATT, 2, , , 0,
ALLSEL,ALL  
/PNUM,MAT,1
/REPLOT
! mesh
smrtsize,4
mshape,1,2d
MSHKEY,0
amesh,all
!load
*GET,Sarea,AREA,3,AREA
if=2
jf=if*turn/sarea
/solu   
esel,s,mat,,3   
bfe,all,js,1,,,jf   
esel,s,mat,,4   
bfe,all,js,1,,,-jf  
esel,all
! boundary
CSYS,1  
NSEL,S,LOC,x,D1/2
D,all,AZ,0,
nsel,all
CSYS,0
! solution
MAGSOLV,0,-1,0.001, ,25,