一个产品的3d静磁场计算讨论（带apdl命令流）

mwei

该产品设计有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

mwei

模型见图，未显示空气

kamakot

你好像没有做两侧边界的自由度耦合，命令是cpcyc。我做过类似的1/20模型，是用到了这个命令的。