ansys定义矩形线圈，磁场很奇怪，求指导交流~~~
代码如下：
finish
/cle
/prep7
Npts=21
i=1
Pi=acos(-1)
mu0=4*Pi*1e-7
! permeability of free space,H/m
mu0r=1
! relative permeability offree space
Br=1.465
! residual induction,Tesla
AO=15 $ AI=6.5 $ L=10
! lower magnet dimensions
d_x=(i-1)*2*ao/(npts-1)
d_Z=5
N=800
CURNT=1
J=n*i*1e6/((ao-ai)*l)
et,1,236
keyopt,1,7,1
MP,MURX,1,1
MP,MURX,2,1.3
MP,MGZZ,2,Br/(mu0*1.3)
AP=15 $ LP=6
! upper magnet dimensions
xo=D_X $ yo=0 $ zo=D_Z+LP/2+L/2
!upper magnet center offset
block,ai,ao,-ai,ai,-l/2,l/2
block,-ai,-ao,-ai,ai,-l/2,l/2
block,-ai,ai,ai,ao,-l/2,l/2
block,-ai,ai,-ai,-ao,-l/2,l/2
wpoff,ai,ai
cylind,0,ao-ai,-l/2,l/2,0,90
wpoff,-2*ai
cylind,0,ao-ai,-l/2,l/2,90,180
wpoff,0,-2*ai
cylind,0,ao-ai,-l/2,l/2,180,270
wpoff,2*ai
cylind,0,ao-ai,-l/2,l/2,270,360
vglue,all
block,-2*AO,XO+2*AO,-2*AO,2*AO,-3*L,2.5*L+0.5*LP+ZO
vovlap,all
mshk,0
msha,1,3d
esize,1
type,1
mat,1
vmesh,1,2
vmesh,9,14
esize,10
vmesh,4
vsel,,,,1
eslv,r
bfe,all,js,1,,j
allsel
vsel,,,,2
eslv,r
bfe,all,js,1,,-j
allsel
vsel,,,,13
eslv,r
bfe,all,js,1,-j,0
allsel
vsel,,,,14
eslv,r
bfe,all,js,1,j
allsel
wpoff,ai,ai
CSWPLA,51,1,1,1
VSEL,S,,,9
ESLV,R
EMODIF,ALL,ESYS,51
BFE,ALL,JS,1,,J
ALLSEL
wpoff,-2*ai
CSWPLA,61,1,1,1
VSEL,S,,,10
ESLV,R
EMODIF,ALL,ESYS,61
BFE,ALL,JS,1,,J
ALLSEL
wpoff,0,-2*ai
CSWPLA,71,1,1,1
VSEL,S,,,11
ESLV,R
EMODIF,ALL,ESYS,71
BFE,ALL,JS,1,,J
ALLSEL
wpoff,2*ai
CSWPLA,81,1,1,1
VSEL,S,,,12
ESLV,R
EMODIF,ALL,ESYS,81
BFE,ALL,JS,1,,J
ALLSEL
nsel,s,ext
d,all,az,0
! flux parallel boundary conditions
alls
vlscale,all,,,1e-3,1e-3,1e-3,,,1
! scale the volumes to mm
alls
fini
/SOLU
solve
fini
/post1
PLVECT,B, , ,,VECT,ELEM,ON,0
gcjyzdd 
