谐波场中用ETABLE求线圈电感电阻的问题（附命令流）

飞，走了

用plane53单元，通过ETABLE可以直接求出电感。加载电压电流均可。但是我求得的结果为什么都是零？在论坛上找了好多相关的帖子，好多朋友都有相同的疑惑，但都没人解决。希望版主高手们指点一二，不胜感激！以下是命令流：

finish
/clear,nostart
/FILNAME,2D_COIL,1
/TITLE,Current-fed thick massive conductor in free space
KEYW,MAGNOD,1
/PREP7
nui=50
mu0=3.1415926*4.0e-7
muri=1/nui/mu0
ET,1,53,,,1            !空气
ET,2,53,2,,1          !线圈
KEYOPT,2,1,1
ET,3,110,,,1          !远场
ET,4,53,,,,1          !铁芯
EMUNIT,MKS        !指定单位制
MP,MURX,1,1        !空气特性
MP,MURX,2,1         !线圈特性
MP,RSVX,2,3e-8       !线圈特性   
MP,MURX,3,muri     !铁芯特性  
Rr=662
DX=10.55
DY=14.42
DYY=5   
M=0
N=0
RECTNG,0,600,0,DYY/2+3*(DY+DYY)+DY
*DO,N,0,3
*DO,M,0,3
X1=Rr+N*DX
Y1=DYY/2+M*(DY+DYY)
X2=X1+DX
Y2=Y1+DY
RECTNG,X1,X2,Y1,Y2
*ENDDO
*ENDDO
PCIRC,0,800,0,90
PCIRC,0,1200,0,90
AOVLAP,ALL
NUMCMP,AREA
/PNUM,AREA,ON
ASEL,S,AREA,,1          !给铁芯分配属性
AATT,3,1,4
ASEL,S,AREA,,2,17     !给线圈分配属性
AATT,2,1,2
ASUM
*GET,A,AREA,,AREA     !得到一个线圈的面积 赋值给A
ASEL,S,AREA,,18          !给远场分配属性
AATT,1,1,3
ASEL,S,AREA,,19          !给空气分配属性
AATT,1,1,1
ASEL,ALL
CSYS,1
LSEL,S,LOC,X,1000
LESIZE,ALL,,,1
ESIZE,,10
MSHAPE,0,2D
MSHKEY,1
ASEL,S,AREA,,18
AMESH,ALL          !划分远场
ALLSEL,ALL
SMRTSIZ,6
MSHAPE,1,2D
MSHKEY,0
ASEL,S,AREA,,1,17
ASEL,A,AREA,,19
AMESH,ALL               !划分空气和线圈
AREFINE,1,17,1,4      !线圈网格优化
ALLSEL,ALL
ASEL,R,,,2
NSLA,R,,2
CP,1,VOLT,ALL
*get,n1,node,,num,min  ! get a node in area 1
F,n1,AMPS,5e8
ALLSEL,ALL
*DO,M,3,17
ASEL,R,,,M
NSLA,R,,M
CP,M,VOLT,ALL
ALLSEL,ALL
*ENDDO
CSYS,1
NSEL,S,LOC,X,1200
SF,ALL,INF      !选择远场上的所有节点，并施加远场标志
CSYS,0
NSEL,S,LOC,X,0
D,ALL,AZ,0              !选择Y轴上的所有节点，并施加边界条件为磁通量与Y轴平行
ALLSEL
FINISH
/SOLU
ANTYPE,3
HARFRQ,1e5
SOLVE
SAVE

FINISH
/POST1
SET,LAST
asel,s,,,2
!esla,s,1
etable,ind,nmisc,9
etable,res,nmisc,8
SSUM
*GET,CRES,SSUM,,ITEM,RES ! GET COIL RESISTANCE
*GET,CIND,SSUM,,ITEM,IND ! GET COIL INDUCTANCE

天涯剑客

没有运行你的程序，但发现几个问题：1、ET,2,53,2,,1          !线圈，这是载压螺线圈，要加载电压，耦合电流自由度，和你后面的加载不符合。2、ET,4,53,,,,1          !铁芯，53后面多了一个逗号。