!the calculation of single-phase static electromagnetic field in 

!salient machine
!there two cycles,one is applied current density,the other is degree of 

!rator 
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
cyl4,0,0,d1/2,,d4/2 
CYL4,0,0,d3/2, ,d4/2
aovlap,all

!**********************************
RECTNG,-a2,a2,-b2,b2,   
CYL4,0,0,d4/2 
AINA,1,2 
CSYS,1 
AGEN,2,3, , , ,90, , ,0  
aadd,1,3
csys,0
CYL4,0,0,d5/2 
aadd,1,2
k,0,0,0
circle,61,d6/2
csys,1
lsel,s,loc,x,d6/2
csys,0
asbl,3,all
asel,all 
cyl4,0,0,d4/2
asel,s,area,,1,3
aovlop,all
numcmp,all
/pnum,area,1
SAVE
wpro,,90
wpro,,,15
asel,s,,,5,6
ASBW,all
wpro,,,60
asel,s,,,4,7,3
asbw,all
wpro,,,60
asel,s,,,3,8,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,,9,12
asel,a,area,,1
asel,a,area,,7
aatt,1
asel,s,area,,2
asel,a,area,,8
aatt,2
asel,s,area,,3,6
asel,a,area,,15,16
aatt,3
asel,s,area,,13,14
asel,a,area,,17,20
aatt,4
/pnum,mat,1
asel,all
/replot
asel,s,area,,7,8
asel,a,area,,13,20
asel,a,area,,3,6
cm,dingzi,area
ET,1,PLANE53
SMRT,3  
MSHAPE,1,2D 
MSHKEY,0
amesh,all
csys,1
nsel,s,loc,x,d4/2-0.00001,d4/2+0.00001
csys,0
cm,nnn,node
asel,s,area,,9,12
asel,a,area,,1
asel,a,area,,2
cm,rato,area
smrt,off
amesh,all
esla,s
cm,rator,elem
allsel,all
save
!******************************************
/solu
antype,0
lsel,s,line,,57,64
DL,all, ,AZ,0,0
ncurr=4
*dim,jf,,ncurr
jf(1)=10,20,35,50
*GET,sarea,AREA,4,AREA
*dim,cur,,1
symfac=1 
asel,s,area,,15,16
asel,a,area,,19,20
esla,s
cm,wind1,elem
allsel,all
*do,jjj,1,6
   cedele,all,,,any
   cmsel,s,rato
   csys,1
   AGEN, ,all, , , ,15*(jjj-1), , , ,1
   csys,0
   cmsel,s,nnn
   cmsel,a,rator
   CEINTF,0.25,az, , , , , ,0,
   *do,iii,1,ncurr
      cur(1)=jf(iii)   
      jff=jf(iii)*turn/sarea   
      asel,s,area,,15
      asel,a,area,,20
      bfa,all,js,,,jff
      asel,s,area,,16
      asel,a,area,,19
      bfa,all,js,,,-jff
      allsel,all
      EQSLV,FRONT,,0
      PIVCHECK,1
      MAGSOLV,0,3,0.001,,25
      lmatrix,symfac,'wind','cur','ind'  
   *enddo
*enddo
finish