finish
/filname,eddy current resisdent 20141204
chang=40 !待测件长
kuan=40 !待测件宽
gao=1.25 !待测件搞
r1=0.6 !线圈内径
th=1.0 !线圈厚度
hc=0.8 !线圈高度
d=0.5 !线圈距待测件距离
turn=140 !线圈匝数
fill=0.95 !线圈填充因子
/prep7
et,1,solid97,1 !待测件
et,2,solid97,2 !线圈
et,3,solid97 !空气
mp,murx,1,1
mp,rsvx,1,1e-6
mp,murx,2,1
mp,rsvx,2,1.73e-8
mp,murx,3,1
area=th*hc*1e-6 !线圈截面积
pi=3.1415926
a_coil=pi*(r1+th)*(r1+th)-pi*r1*r1
v_coil=a_coil*hc*1e-9 !线圈体积
R,1,area,turn,v_coil,0,1,0,
RMORE,1,fill
block,-chang/2,chang/2,-kuan/2,kuan/2,0,gao
wpoffs,0,0,gao+d
cyl4,0,0,r1,0,r1+th,360,hc
CSWPLA,11,1,1,1
block,-chang,chang,-kuan,kuan,-10*hc,10*hc
vovlap,all
numcmp,all
/TRLCY,VOLU,0.8,3
/pnum,volu,on
vsel,s,,,1
vatt,1,,1,0 !1号材料,1号单元类型
aslv,s
lsla,s
csys,0
lsel,r,loc,z,gao/2
lesize,all,,,8 !把线分成8分
esize,0.4 !单元边长0.4
vsweep,all !扫掠
vsel,s,,,2
vatt,2,1,2,11
esize,0.2
mshkey,0 !自由划分
mshape,1,3d !四面体
vmesh,all
vsel,s,,,3
vatt,3,,3
esize,2 !空气单元尺寸
!smrtsize,2,,,1.2,15,,,on,on
vmesh,3
alls
VLSCAL,all, , ,0.001,0.001,0.001, ,0,1 !模型所有尺寸缩小为原来的0.001,即单位变成m
!***************
/solu
alls
asel,s,ext
da,all,ax,0
da,all,ay,0
da,all,az,0 !空气最外层认为磁场衰减为零
vsel,s,,,2
eslv,s
nsle,s
cp,1,curr,all !耦合线圈电流自由度
BFE,all,VLTG,1,12*1.414 !加电压,峰值
alls
antype,harm !谐波分析
harfrq,150000 !频率150K
solve
!*******************此处开始后处理
/post1
alls
esel,s,mat,,2
SET,1,LAST,1,1, , ,
PRNSOL,CURR !电流虚部值,保存
SET,1,LAST,1,0, , ,
PRNSOL,CURR !电流实部值,保存