关于电磁场线圈1/4模型的边界条件问题

i我飞翔

大家好，本人水货加小鲜肉，初学ansys 这个高大上的软件，想建立一个关于管道用线圈磁化的三维局部1/4模型，但是不知道怎么加载边界条件，加载电流部分已经解决，但是边界条件设置缺遇到问题，用的是静态分析和97号单元，请大神和朋友们帮助，先谢谢大家了，以下是个人仿真的程序，希望大家能够指点一二，以求共同进步，
/clear
/filname,xxx
r1=0.5                        !线圈内径
th=1.85                        !线圈厚度
hc=5.1                        !线圈高度
rd=9            !导体半径
gao=12            !导体高
d=0.5                !线圈距导体距离
turn=400        !线圈匝数
fill=0.95473        !线圈填充因子
/prep7
et,1,solid97,1        !导体
et,2,solid97,2        !线圈
et,3,solid97        !空气

pi=3.1415926

mp,murx,1,1
mp,rsvx,1,2.7e-8

mp,murx,2,1
mp,rsvx,2,1.7241e-8

mp,murx,3,1

area=th*hc*1e-6        !线圈截面积

a_coil=pi*(r1+th)*(r1+th)-pi*r1*r1        
v_coil=(a_coil*hc*1e-9)/4        !线圈体积

R,1,area,turn,v_coil,0,1,0,  
RMORE,4,fill                 


cyl4,0,0,0,0,rd,90,gao                  
wpoffs,0,0,gao+d

cyl4,0,0,r1,0,r1+th,90,hc                  !线圈

wpoffs,0,0,-12.5

SPHERE,25,0,0,90                         !空气场
SPHERE,50,0,0,90

/PNUM,VOLU,1  !显示实体编号
/REP,FAST
NUMCMP,VOLU
/REP,FAST

CSWPLA,11,1,1,1

CSWPLA,12,2,1,1

csys,0

alls
vovlap,all
numcmp,all


vsel,s,,,1      !划分被测体        
vatt,1,,1,0                 !vatt-定义所选单元属性  mat，real，type，esys
esize,0.75                  !mat（material number）材料号，real（real constant set number）实常数，
mshkey,1                    !type（element type number）单元类型编号，esys-单元坐标类型
mshape,0,3d
vmesh,all


vsel,s,,,2              
vatt,2,1,2,11
esize,0.75
mshkey,1
mshape,0,3d
vmesh,all


vsel,s,,,4               !划分近场空气
vatt,3,,3,12
esize,1.5
mshkey,0
mshape,1,3d
vmesh,all


vsel,s,,,3               !划分远场空气
vatt,3,,3,12
esize,4.5
mshkey,0
mshape,1,3d
vmesh,all


alls
VLSCAL,all, , ,0.001,0.001,0.001, ,0,1    !模型缩小，米制
!******************************************扫掠划分远场空气测试
!vsel,s,,,3               !划分远场空气
!vatt,3,,3,12
!vatt,4,,4,12
!esize,4.5
!esize,40                      !infin111作为远场单元，只能划分一层单元，不能为多层单元
!vsel,s,,,3
!vsweep,3
!***********************************************************************开始加载求解
/solu
ALLSEL,ALL
CSYS,0
!!!!!!!!!!=========================AZ,AY,AZ=0    1.7，1.5
!Asel,s,ext                        !!!!!!!
!da,all,ax,0                       !!!!!!!
!da,all,ay,0                       !!!!!!!
!da,all,az,0                       !!!!!!!!
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!

!!!!!!!!!=====================  0.029，0.15
!ASEL,S,AREA,,13
!da,all,ax,0
!da,all,ay,0
!da,all,az,0

!ASEL,S,LOC,X,0
!DA,ALL,AZ,0
!ASEL,S,LOC,Y,
!DA,ALL,AZ,0
!!!!!!!!!=====================  1.09,1.44
!ASEL,S,AREA,,13
!da,all,ax,0
!da,all,ay,0
!da,all,az,0

!ASEL,S,LOC,X,0
!DA,ALL,AX,0
!ASEL,S,LOC,Y,
!DA,ALL,AY,0
!!!!!!!!!========================= 0.016，0.1
!ASEL,S,AREA,,13
!da,all,ax,0
!da,all,ay,0
!da,all,az,0
!lsel,s,loc,,x,0
!lsel,r,loc,,y,0
!nsll,s,
!d,all,az,0
!!!!!!!!!=====================
!ASEL,S,LOC,X,0
!da,all,aX,0
!ASEL,S,LOC,Y,0
!da,all,aY,0

esel,s,mat,,1
nsle,s
cp,1,volt,all

!vsel,s,,,1
esel,s,mat,,2
!eslv,s                          !选中依附于体积的单元
nsle,s                          !选中附着在单元上的节点
cp,2,curr,all                   !耦合线圈电流自由度
BFE,all,VLTG,1,24*1.414         !加电压，峰值

alls
antype,harm
harfrq,10000
solve

/post1
esel,s,mat,,2
SET,1,LAST,1,1, , ,
PRNSOL,CURR                !电流虚部值，保存
SET,1,LAST,1,0, , ,
PRNSOL,CURR                 !电流实部值，保存


自己设置了几个边界条件，但都好像不是正确结果，请大家指点，谢谢，附上几个图，大家看下，
file:///C:\Users\iwofeixiang\AppData\Roaming\Tencent\Users\991204861\QQ\WinTemp\RichOle\OUGMA8`$5)AMQOA5UDX_(GE.png

完整模型仿真结果

局部模型仿真结果