该产品设计有5个线圈 ,1个内线圈位于中心,4个外线圈均布在外圆上,通过调节这5个线圈的电流大小得到希望的磁场。
完整的apdl如下,感觉结果不对,哪位帮我看看,问题出在哪里?
个人觉得取的1/4模型有问题,边界条件设置也有些问题。
fini
/clear
/CWD,'D:\ansys'
/prep7
/title, 3-D Static Magnetic Problem
/com,
/units,si
!
! define analysis parameters
!
N1=61 ! inner coil turns
N2=38 ! outer coil turns
Id=4.72 ! current per turn unit: A
L=0.04 !channel length
W=0.022 !channel width,0.25d
d=0.088 !channel mean diameter
Lth=d !thruster length,0.88
dth=2*d !thruster diameter,0.176
douc=d+W !channel outer diameter,0.11
dinc=d-W !channel inner diameter,0.066
deltaw=0.009 !channel thick
dou=douc+deltaw*2 !0.128
din=dinc-deltaw*2 !0.048
deltabp=0.0014 !back pole thick
deltafpin=0.003 !inner front pole thick
deltafpou=0.003 !outer front pole thick
deltasin=0.001 !inner magnetic screen thick
lsin=0.005 !space between inner magnetic screen and inner front pole
Lscreenin=Lth-lsin-deltabp-deltafpin !inner magnetic screen length,0.0786
deltasou=0.001 !outer magnetic screen thick
lsou=0.005 !space between outer magnetic screen and outer front pole
Lscreenou=Lth-lsou-deltabp-deltafpou !outer magnetic screen length,0.0786
dcoi=0.0092 !inner coil diameter
Lcoi=Lth-deltabp-deltafpin !inner coil length
dcoo=0.0062 !outer coil diameter
Lcoo=Lth-deltabp-deltafpou !outer coil length
js1=N1*Id/(dcoi/2*Lcoi)
js2=N2*Id/(dcoo/2*Lcoo)
/prep7
/pnum,volu,1
csys,0
!back pole plate
cylind,dth/2,0,0,deltabp,0,90
!inner screen
cylind,din/2,din/2-deltasin,deltabp,deltabp+Lscreenin,0,90
!outer screen
cylind,dou/2,dou/2+deltasou,deltabp,deltabp+Lscreenou,0,90
!inner channel
cylind,din/2,dinc/2,Lth-L,Lth,0,90
!outer channel
cylind,douc/2,dou/2,Lth-L,Lth,0,90
!inner iron
cylind,dcoi/2,0,deltabp,deltabp+Lcoi,0,90
!inner coil
cylind,dcoi/2,dcoi/2*2,deltabp,deltabp+Lcoi,0,90
!outer iron
wpoffs,0,dth/2-dcoo/2,deltabp
csys,wp
cyl4,0,0,dcoo/2,-90,dcoo/2,90,Lcoo
cyl4,0,0,dcoo/2,-90,dcoo/2*2,90,Lcoo
!CYL4, XCENTER, YCENTER, RAD1, THETA1, RAD2, THETA2, DEPTH
wpave,0,0,0
WPCSYS,1,0
wpoffs,dth/2-dcoo/2,0,deltabp
csys,wp
cyl4,0,0,dcoo/2,0,dcoo/2,180,Lcoo
cyl4,0,0,dcoo/2,0,dcoo/2*2,180,Lcoo
!CYL4, XCENTER, YCENTER, RAD1, THETA1, RAD2, THETA2, DEPTH
csys,0
wpave,0,0,0
cylind,din/2,0,Lth-deltafpin,Lth,0,90
cylind,dth/2,dou/2,Lth-deltafpou,Lth,0,90
!enclosure
cylind,dth/2*1.2,0,0-Lth/5,Lth*1.5,0,90
vovlap,all
numcmp,volu
!define element type
!
et,1,236
!
!define material properties for air and steel
mp,murx,1,1 !air
tb,bh,2,,40 !steel
tbpt,,355,.7
,,405,.8
,,470,.9
,,555,1.0
,,673,1.1
,,836,1.2
,,1065,1.3
,,1220,1.35
,,1420,1.4
,,1720,1.45
,,2130,1.5
,,2670,1.55
,,3480,1.6
,,4500,1.65
,,5950,1.7
,,7650,1.75
,,10100,1.8
,,13000,1.85
,,15900,1.9
,,21100,1.95
,,26300,2.0
,,32900,2.05
,,42700,2.1
,,61700,2.15
,,84300,2.2
,,110000,2.25
,,135000,2.3
,,200000,2.41
,,400000,2.69
,,800000,3.22
tbcopy,bh,2,3
!
!
!set volume attributes
!
vsel,s,volu,,1,3
vatt,2,1,1 ! pole
vsel,s,volu,,5,7
vatt,2,1,1 ! magnetic material
vsel,s,volu,,13,14
vatt,2,1,1 ! magnetic material
!
vsel,s,volu,,4
vatt,1,1,1 ! innercoil
vsel,s,volu,,8,9
vatt,1,1,1 ! outercoil
vsel,s,volu,,10,12
vatt,1,1,1 ! air
vsel,s,volu,,15
vatt,1,1,1 ! air
!mesh the model
!
allsel,all
smrt,6 ! for production work a finer mesh should be used, (LVL 6)
mshape,1,3d
mshkey,0
vmesh,all
/pnum,mat,1
/number,1
eplot
!create the inner coil
!
!
local,11,1,0,0,deltabp+Lcoi/2 !顺时针
vsel,s,,,4
eslv,s
emodif,all,esys,11 ! EMODIF, IEL, STLOC, I1,...
/PSYMB,ESYS,1
eplot
BFE,ALL,JS,,0,js1,0 ! apply current density on all selected elements
/PBF,JS, ,1
/PBC,ALL, ,1
/rep
!create the outer coil1
!
local,21,1,0,dth/2-dcoo/2,deltabp+Lcoo/2,,,
vsel,s,,,8
eslv,s
emodif,all,esys,21 ! EMODIF, IEL, STLOC, I1,...
/PSYMB,ESYS,1
BFE,ALL,JS,,0,-js2,0 ! apply current density on all selected elements
/PBF,JS, ,1
/PBC,ALL, ,1
/rep
!create the outer coil2
local,22,1,dth/2-dcoo/2,0,deltabp+Lcoo/2,,,
vsel,s,,,9
eslv,s
emodif,all,esys,22 ! EMODIF, IEL, STLOC, I1,...
/PSYMB,ESYS,1
BFE,ALL,JS,,0,-js2,0 ! apply current density on all selected elements
/PBF,JS, ,1
/PBC,ALL, ,1
eplot
alls
save
finish
!
!apply boundary conditions
!
/solu
csys,1
! Boundary Conditions
!nsel,s,loc,x,dth/2*1.2
!nsel,a,loc,z,-Lth/5
!nsel,a,loc,z,Lth*1.5
nsel,s,loc,y,0
nsel,a,loc,y,90
d,all,az,0
alls
fini
! Set potentials to zero (flux-parallel)
/PBC,ALL, ,1
eplot
!
!solve
!
csys,0
allsel,all
magsolv,0,,,,,1
finish
save
!
!summarize armature force results
!
/post1
SET,LAST
wpave,0,0,0
wprota,,90
SUCR,111,CPLANE,,
SUSEL,R,111
SUMAP,222,B,SUM
supl,111,222,0
!suvect,111,
rsys,1
PATH,channel,2,30,20,
PPATH,1,0,d/2,0,0,0,
PPATH,2,0,d/2,0,Lth*1.5,0,
PDEF, ,B,sum,AVG
PLPATH,BSUM
PDEF, ,B,X,AVG
PLPATH,BX
PDEF, ,B,y,AVG
PLPATH,By
PDEF, ,B,z,AVG
PLPATH,Bz
|