/title,keti
/filname,keti,1
keyw,magnod,1
/prep7
ET,1,PLANE53 !空气
KEYOPT,1,3,1
ET,2,PLANE53 !线圈
KEYOPT,2,3,1
ET,3,INFIN110
KEYOPT,3,3,1 ! 远场
MP,MURX,1,1 !空气材料
MP,MURX,2,1 !线圈材料
MP,RSVX,2,1/(5.8*10E7)
!!!!!!!!!!!!!!!!!!!!
S=0.2
N1=500
N2=450
N3=400
R=1.5*S
*dim,wind,char,3
RECTNG,S,2*S,0,S
RECTNG,2.5*S,3.5*S,0,S
RECTNG,4*S,5*S,0,S
CYL4,0,0,0,0,18*S,90
CYL4,0,0,0,0,12*S,90
ASEL,ALL
AOVLAP,ALL
!!!!!!!!!!!!!!!!!!!!!!
/pnum,area,1 !!!!进行剖分
aplot,all
numcmp,area
aplot,all
r,1,2*s*s,500,,1,1
r,2,2*s*s,450,,1,1
r,3,2*s*s,400,,1,1
asel,s,area,,1
aatt,2,1,2
smrtsize,2
mshape,0,2d
mshkey,0
amesh,1
asel,s,area,,2
aatt,2,2,2
smrtsize,2
mshape,0,2d
mshkey,0
amesh,2
asel,s,area,,3
aatt,2,3,2
smrtsize,2
mshape,0,2d
mshkey,0
amesh,3
allsel,all
aplot,all
SMRT,OFF
ASEL,S,AREA,,4
AATT,1,,3
LOCAL,11,1
WPCSYS,,11
LSEL,S,LOC,Y,45
LESIZE,ALL,,,10
ESIZE,,1
AMESH,4
ALLSEL,ALL
SAVE
APLOT,ALL
ASEL,S,AREA,,5
smrt,on
SMRTSIZE,5
MSHKEY,0
MSHAPE,1,2D
AMESH,5
allsel,all
finish
!!!!!!!!!!!!!!!!!!!!!!!!!!
!!!!!!!!!!!!!!!!!进行求解
/SOLU
cnvtol,csg,,1.0e-4!!!
eqslv,front
CSYS,0
NSEL,S,LOC,X,0
D,ALL,AZ,0
CSYS,1
NSEL,S,LOC,X,18*S
SF,ALL,INF
CSYS,0
ALLSEL,ALL
SAVE
asel,s,area,,1
esla,s
cm,wind1,elem
bfe,all,js,,,,1.2*500/(s*s)
asel,s,area,,2
esla,s
cm,wind2,elem
bfe,all,js,,,,1.5*450/(s*s)
!bfe,all,volt,,12
asel,s,area,,3
esla,s
cm,wind3,elem
bfe,all,js,,,,1.7*400/(s*s)
allsel,all
antype,static
magsolv
finish
!!!!!!!!!!!!!!!!!!!!!!!!!
/post1
allsel,all
plf2d
finish
/solu
*dim,cur,array,3
cur(1)=1.2
cur(2)=1.5
cur(3)=1.7
lmatrix,1,'wind','cur','ind'
*stat,ind |