涡流仿真，有程序

jianglang

哪位大神指点一下，为什么我求得的实部和虚部电流会这么小？
/PREP7
!材料类型
ET,1,PLANE53,0,,1     !空气
ET,2,PLANE53,2,,0     !线圈
ET,3,PLANE53,4,,0     !铝板
EMUNIT,MKS
!线圈实常数

N=400
PI=4*ATAN(1)
R,1,1.35E-5,N,,1,1   !
R,2,1.04E-4,,,1

MP,MURX,1,1           !空气材料特性

MP,MURX,2,1             !线圈材料特性
MP,RSVX,2,3.04878E-8

MP,MURX,3,1            !铝板材料特性
MP,RSVX,3,2.7E-7

/PNUM,AREA,1
/PNUM,AREA,1

!建模

l=0.0005
Rin=0.0001
Rout=0.0016
RECTNG,0,0.012,0,0.009                          !铝板
RECTNG,Rin+12*l,Rout+12*l,0.0095,0.0185                   !线圈
RECTNG,0.00605,0.00606,0.005,0.009              !裂纹
PCIRC,0,0.13,0,90                               !空气

AOVLAP,ALL
APLOT,ALL

!属性赋予
ASEL,S,AREA,,2     !线圈
AATT,2,1,2
ASUM
!*GET,A,AREA,,AREA            ! AREA OF COIL CROSS-SECTION
!*GET,COIL_XC,AREA,1,CENT,X   !获取线圈的质心
!CF=2*PI*(3.04878E-8)*COIL_XC*N**2/A*2
ASEL,S,AREA,,6     !铝板
AATT,3,2,3

ASEL,S,AREA,,7     !空气
AATT,1,,1
ASEL,S,AREA,,5     !裂纹
AATT,1,,1

/PNUM,MAT,1
ALLSEL,ALL

!CSYS,1
!网格划分

MSHKEY,1                         !线圈线划分
ESIZE,,18
LESIZE,6,,,18
MSHKEY,1
ESIZE,,2
LESIZE,7,,,2
AMESH,2                 


MSHKEY,1                       !裂缝线划分
ESIZE,,1
LESIZE,9,,,1
MSHKEY,1
ESIZE,,8
LESIZE,10,,,8
AMESH,5

ASEL,S,AREA,,6                !铝板自由网格划分
LCCAT,12,17
LCCAT,9,12
LCCAT,9,10
LCCAT,10,16
MSHKEY,0
MSHAPE,1,2D
AMESH,6


ASEL,S,AREA,,7
!SMRTSIZE,0.5
MSHKEY,0
ESIZE,,121
LESIZE,19,,,61
MSHKEY,0
ESIZE,,118
LESIZE,18,,,59
MSHKEY,0
ESIZE,,150
LESIZE,13,,,75
MSHAPE,1,2D
AMESH,7                         !空气区网格划分


!创建电路部分
N,100001,0.2,0,0
N,100002,0.2,0.1,0
N,100003,0.2,0.05,0
N,100004,0.14,0,0
N,100005,0.14,0.1,0



AV=12
ET,4,CIRCU124,4,0          !独立电压源
TYPE,4
REAL,4
MAT,2
R,4,AV
E,100002,100001,100003
  
ET,5,CIRCU124,0,0           !电阻
TYPE,5
REAL,5
MAT,2
R,5,2E-8
E,100002,100005

ET,6,CIRCU124,5,0           !线圈
TYPE,6
REAL,1
MAT,2
E,100005,100004,88
/SOLU
ANTYP,HARM
HARFRQ,1000
!施加0电位边界条件
D,100001,VOLT,0
D,100004,VOLT,0

!耦合边界条件
ALLSEL,ALL                  !耦合线圈
ASEL,S,AREA,,2
NSLA,S,1
CP,1,CURR,ALL
!*GET,ELM,ELEM,,NUM,MIN      ! GET AN ELEMENT NUMBER IN THE COIL REGION
!CP,2,EMF,ALL

ALLSEL,ALL                  !耦合铝板
ASEL,S,AREA,,6
ESEL,S,MAT,,3
NSLA,S,1
CP,3,CURR,ALL
CP,4,EMF,ALL
ALLSEL,ALL
EPLOT,ALL



DL,13,7,AZ,0
DL,18,7,AZ,0
DL,19,7,AZ,0
DL,4,6,AZ,0
DL,1,6,AZ,0



SOLVE


!/CLEAR
!*DO,l,0,0.012
!l=l+0.001
!*enddo

/POST1
SET,1                        !读入第一步结果实部
PLF2D                        !绘制二维实部磁力线图
*GET,CURREAL,NODE,88,CURR
SET,1,,,1                   !读入第一步结果
PLF2D                         !绘制二维虚部磁力线图
*GET,CURIMAG,NODE,88,CURR
*SET,Vreal,12*SQRT(2)
Imag=CURREAL**2+CURIMAG**2
Zimag=Vreal*CURIMAG/Imag            !阻抗虚部
Zreal=Vreal*CURREAL/Imag            !阻抗实部


!ESEL,S,MAT,,3               !被测对象铝板
!ETABLE,JE,NMISC,6           !存储涡流密度
!PLETAB,JE,1                 !显示被测对象中的涡流密度（云图）
!PRETAB,JE                   !列表显示被测对象中的涡流密度

ESEL,S,MAT,,2                !选择探头
ETABLE,RES,NMISC,8          !存储探头线圈各单元直流电阻
ETABLE,IND,NMISC,9          !存储探头线圈各单元电感
SSUM
*GET,CRES,SSUM,,ITEM,RES     ! GET COIL RESISTANCE
*GET,CIND,SSUM,,ITEM,IND     ! GET COIL INDUCTANCE
CRES=CRES                  ! COIL RESISTANCE
CIND=CIND                  ! COIL INDUCTANCE
!*GET,NCUR,NODE,N2,CURR   ! GET SOLUTION CURRENT
!ESEL,ALL
!FINISH
*enddo
*GET,Zimag,ESEL,2，
*GET,Zreal,ESEL,2，

201010202012

创建的电路部分起什么作用呢 还有没搞懂 电压是怎么加上去的

dzy

同问，又没有大神可以指导一下