ansys线圈阻抗的计算

没什麽是绝对的

高频线圈靠近一个被测体时，线圈自身阻抗改变，用的ansys，我对这个没学过，下面是我写的一段程序，求大神们帮我改改，谢谢咯
FINISH
/CLEAR
/PREP7
ET,1,PLANE53,,,1            ！定义空气元素
ET,2,PLANE53,2,,1            ！定义线圈元素
ET,3,PLANE53,4,,1            ！定义被测体元素

EMUNIT,MKS
MP,MURX,1,1                  ！相对磁导率1
MP,MURX,2,1
MP,RSVX,2,3.00E-8           ！定义电阻率
MP,MURX,3,1.01
MP,RSVX,3,7.50E-7

RECTNG,0.3/1000,1/1000,10.5/1000,11.3/1000     ！线圈尺寸
RECTNG,0,20/1000,0,10/1000                  ！被测体尺寸
PCIRC,0,40/1000,0,90                         ！空气域尺寸
AOVLAP,ALL                             ！重叠实体面
NUMCMP,AREA
/PNUM,AREA,1

ASEL,S,AREA,,1                  ！选线圈面对其赋值
AATT,2,1,2
ASUM
*GET,A,AREA,,AREA            ！计算线圈横截面积
ASEL,S,AREA,,3
AATT,1,1,1
ASEL,S,AREA,,2
AATT,3,1,3
ASEL,ALL           
CSYS,1            ！改为柱坐标
LSEL,S,LOC,X,30/1000      
LESIZE,ALL,,,1
ESIZE,,10
MSHAPE,0,2d
MSHKEY,1
AMESH,1,2       ！为面1,2划分网格，每条线等分为10等份
SMRTSIZE,2
MSHKEY,0
AMESH,3           ！为面3划分网格

R,1,a,45,,1,1            ！给线圈赋实数
N1=NODE(0.3,10.5,0)     ！选一个线圈面上的节点
ESEL,S,MAT,,2         
NSLE,S                 ！选择线圈面上的所有节点
CP,1,CURR,ALL          ！电流自由度耦合
N2=NODE(20,10,0)       ！选一个被测体面上的节点
ESEL,S,MAT,,3           
NSEL,S                 ！选择被测量面上所以的节点
CP,NEXT,CURR,ALL          ！电流自由度耦合
CP,NEXT,EMF,ALL           ！压力自由度耦合

CSYS,1
NSEL,S,LOC,X,40/1000        ！设置无线表面标志
SF,ALL,INF
CSYS,0
LSEL,S,LOC,X,0            
D,ALL,AZ,0            ！加磁力线平等边界条件     
NSEL,ALL
FINISH

/SOLU
ANTYPE,HARM
ESEL,S,MAT,,2          ！选择材料号为2的元素
BFE,ALL,VLTG,,24      ！加24V电压
ESEL,ALL
HARFPQ,1000000         ！频率1MK
SOLV
FINISH
/POST1
SET,1
PLF2D         ！画出实部磁力线