发一个凸极电机谐波磁场仿真的命令流

zhaoxiang

开始学习电机的电磁场分析，用了一上午作了一个例子，又用了半下午写了命令流文件！
fini
/clear
/filname,the calculation of electromagnetic field in salient machine
!*  
/NOPR   
/PMETH,OFF,1
KEYW,PR_SET,1   
KEYW,PR_ELMAG,1
KEYW,MAGNOD,1     
/GO
!*
*afun,deg
*SET,d1,77.5/1000   
*SET,d2,62/1000
*SET,d3,47.5/1000   
*SET,d4,47.2/1000   
*SET,d5,35/1000
*SET,d6,19/1000
*SET,turn,50
*SET,an1,30*0.4463  
*SET,an2,45*0.3273  
*SET,a1,d3/2
*SET,b1,sin(an1)*d3/2   
*SET,a2,d4/2
*SET,b2,sin(an2)*d4/2  
/PREP7
CYL4,0,0,d1/2, ,d2/2
CYL4,0,0,d2/2   
wpoff,0,b1,0
wprot,0,90,0
ASBW,       2
wpoff,0,0,b1*2  
ASBW,       4  
wpoff,0,0,-b1   
wprot,0,-90,0
adel,2,3,,1
CYL4,0,0,d3/2
ASBA,       5,       2
CSYS,1
agen,3,3,4,1,,60
NUMCMP,ALL  
CSYS,0  
aadd,all
SAVE
NUMCMP,ALL
CSYS,1   
AGEN, ,1, , , ,45, , , ,1
CSYS,0
RECTNG,-a2,a2,-b2,b2,   
CYL4,0,0,d4/2
AINA,2,3
CSYS,1
AGEN,2,4, , , ,90, , ,0  
aadd,2,4,2
csys,0
CYL4,0,0,d5/2
aadd,2,3
numcmp,all
k,0,0,0
circle,49,d6/2
csys,1
lsel,s,loc,x,d6/2
csys,0
asbl,4,all
CYL4,0,0,d3/2, ,d4/2
CYL4,0,0,d1/2
asel,all
aovlop,all
SAVE
wpro,,90
wpro,,,15
asel,s,,,12,13
ASBW,all
wpro,,,60
asel,s,,,11,14,3
asbw,all
wpro,,,60
asel,s,,,10,15,5
asbw,all
WPCSYS,-1
numcmp,all  
allsel,all
MPTEMP,1,0  
MPDATA,MURX,1,,1        !air
MPDATA,MURX,2,,0        !iron       
MPDATA,MURX,3,,1        !winding for +Z direction of current
MPDATA,MURX,4,,1        !winding for -Z direction of current
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
asel,s,area,,2
asel,a,area,,6,9
asel,a,area,,14
aatt,1
asel,s,area,,3
asel,a,area,,15
aatt,2
asel,s,area,,13,19,6
asel,a,area,,4,16,12
asel,a,area,,11,18,7
aatt,3
asel,s,area,,10,20,10
asel,a,area,,1,5,4
asel,a,area,,12,17,5
aatt,4
/pnum,mat,1
asel,all
/replot
ET,1,PLANE53
SMRT,3  
MSHAPE,1,2D
MSHKEY,0
amesh,all
/save
*GET,sarea,AREA,4,AREA
*SET,jf,2   
*SET,jf,jf*turn/sarea   
antype,3
lsel,s,line,,85,92
DL,all, ,AZ,0,0
asel,s,area,,13
asel,a,area,,19
bfa,all,js,,,jf,0,
asel,s,area,,10
asel,a,area,,20
bfa,all,js,,,-jf,0,

asel,s,area,,4
asel,a,area,,16
bfa,all,js,,,jf,120
asel,s,area,,1
asel,a,area,,5
bfa,all,js,,,-jf,120

asel,s,area,,11
asel,a,area,,18
bfa,all,js,,,jf,240
asel,s,area,,12
asel,a,area,,17
bfa,all,js,,,-jf,240
allsel,all
HARFRQ,50,50,   
NSUBST, ,   
KBC,0
HMAGSOLV,50,3,0.001,0.001,0.001,0.001,50,   
FINISH  
/POST1  
SET,LAST
PLF2D,27,0,10,1

zhaoxiang

斑竹给加分吗？

zhaoxiang

修正！上述命令流两处错误，进行了改正！并对改变了绕组的接法！
fini
/clear
/filname,the calculation of electromagnetic field in salient machine
!*  
/NOPR   
/PMETH,OFF,1
KEYW,PR_SET,1   
KEYW,PR_ELMAG,1
KEYW,MAGNOD,1     
/GO
!*
*afun,deg
*SET,d1,77.5/1000   
*SET,d2,62/1000
*SET,d3,47.5/1000   
*SET,d4,47.2/1000   
*SET,d5,35/1000
*SET,d6,19/1000
*SET,turn,50
*SET,an1,30*0.4463  
*SET,an2,45*0.3273  
*SET,a1,d3/2
*SET,b1,sin(an1)*d3/2   
*SET,a2,d4/2
*SET,b2,sin(an2)*d4/2  
/PREP7
CYL4,0,0,d1/2, ,d2/2
CYL4,0,0,d2/2   
wpoff,0,b1,0
wprot,0,90,0
ASBW,       2
wpoff,0,0,b1*2  
ASBW,       4  
wpoff,0,0,-b1   
wprot,0,-90,0
adel,2,3,,1
CYL4,0,0,d3/2
ASBA,       5,       2
CSYS,1
agen,3,3,4,1,,60
NUMCMP,ALL  
CSYS,0  
aadd,all
SAVE
NUMCMP,ALL
CSYS,1   
AGEN, ,1, , , ,45, , , ,1
CSYS,0
RECTNG,-a2,a2,-b2,b2,   
CYL4,0,0,d4/2
AINA,2,3
CSYS,1
AGEN,2,4, , , ,90, , ,0  
aadd,2,4,2
csys,0
CYL4,0,0,d5/2
aadd,2,3
k,0,0,0
circle,49,d6/2
csys,1
lsel,s,loc,x,d6/2
csys,0
asbl,4,all
CYL4,0,0,d3/2, ,d4/2
CYL4,0,0,d1/2
asel,all
aovlop,all
SAVE
wpro,,90
wpro,,,15
asel,s,,,12,13
ASBW,all
wpro,,,60
asel,s,,,11,14,3
asbw,all
wpro,,,60
asel,s,,,10,15,5
asbw,all
WPCSYS,-1
numcmp,all  
allsel,all
MPTEMP,1,0  
MPDATA,MURX,1,,1        !air
MPDATA,MURX,2,,0        !iron       
MPDATA,MURX,3,,1        !winding for +Z direction of current
MPDATA,MURX,4,,1        !winding for -Z direction of current
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
asel,s,area,,2
asel,a,area,,6,9
asel,a,area,,14
aatt,1
asel,s,area,,3
asel,a,area,,15
aatt,2
asel,s,area,,13,19,6
asel,a,area,,4,16,12
asel,a,area,,11,18,7
aatt,3
asel,s,area,,10,20,10
asel,a,area,,1,5,4
asel,a,area,,12,17,5
aatt,4
/pnum,mat,1
asel,all
/replot
ET,1,PLANE53
SMRT,3  
MSHAPE,1,2D
MSHKEY,0
amesh,all
save
*GET,sarea,AREA,4,AREA
*SET,jf,2   
*SET,jf,jf*turn/sarea   
antype,3
lsel,s,line,,85,92
DL,all, ,AZ,0,0
!another connection method of winding
asel,s,area,,10
asel,a,area,,13
bfa,all,js,,,jf,0,
asel,s,area,,1
asel,a,area,,4
bfa,all,js,,,-jf,0,

asel,s,area,,5
asel,a,area,,16
bfa,all,js,,,jf,120
asel,s,area,,17
asel,a,area,,18
bfa,all,js,,,-jf,120

asel,s,area,,11
asel,a,area,,12
bfa,all,js,,,jf,240
asel,s,area,,19
asel,a,area,,20
bfa,all,js,,,-jf,240
allsel,all
HARFRQ,50,50,   
NSUBST, ,   
KBC,0
HMAGSOLV,50,3,0.001,0.001,0.001,0.001,50,   
FINISH  
/POST1  
SET,LAST
PLF2D,27,0,10,1

zhaoxiang

自己顶一下啦！

lanny124

好好看看先,谢谢了

lllwww2020

谢谢楼主分享，还有别的电机的命令流可以分享不？

wxch

谢谢！！看看先