/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，