/PREP7
ET,1,PLANE53,,,1 !空气
ET,2,PLANE53,2,,1 !线圈
ET,3,PLANE53,1,,1 !被测体
EMUNIT,MKS
mp,murx,1,1 !空气
mp,murx,2,1 !线圈
mp,rsvx,2,1.71e-8
mp,murx,3,1000000 !镍
mp,rsvx,3,6.844e-8
*SET,coil_D_OUT,18/1000 !线圈外径
*SET,COIL_D_IN,16/1000 !线圈内径
*SET,COIL_H,5/1000 !线圈高
*SET,COIL_AREA,COIL_H*(coil_D_OUT-coil_D_IN)/2 !线圈截面积
*set,coil_d,0.2/1000
*SET,COIL_TURNS,coil_area/(coil_d**2) !线圈匝数取整
*SET,COIN_D,25/1000 !硬币直径
*set,coin_h,2/1000 !硬币厚度
*SET,GAP,3/1000 !检测距离
*SET,coil_VOLT,5 !线圈加载电压
*SET,coil_FREQ,350000 !频率
/pnum,area,1
rectng,coil_d_in/2,coil_d_out/2,-coil_h/2,coil_h/2 !线圈
rectng,0,coin_d/2,-coin_h-coil_h/2-gap,-coil_h/2-gap !被测体
rectng,0,25/1000,-25/1000,25/1000 !空气
aovlap,all
asel,s,area,,4 !空气
aatt,1,,1
asel,s,area,,1 !线圈
aatt,2,,2
asel,s,area,,2 !被测体
aatt,3,,3
asel,all
csys,1
esize,coil_d
mshape,0,2d
mshkey,1
amesh,1,2
smrtsize,2
esize,1/1000
mshape,0,2d
mshkey,0
amesh,4
r,1,coil_area,coil_turns,,1,1
*set,n1,node(coil_d_in/2,0,0)
esel,s,mat,,2 !选中线圈
nsle,s
cp,1,curr,all !线圈耦合电流
esel,s,mat,,3 !选中被测体
nsle,s
cp,2,volt,all !被测体耦合电压
allsel,all
csys,0
nsel,s,loc,x,0 !选中Y轴
d,all,az,0 !设定边界条件为磁通量与Y轴平行
allsel,all
finish
/solu
antype,harm
esel,s,mat,,2 !选中线圈
bfe,all,vltg,,coil_volt !施加载荷
esel,all
harfrq,coil_freq !指定谐性频率
solve
finish
/post1
set,1
*get,curreal,node,n1,curr
set,1,1,,1
*get,curimag,node,n1,curr
*stabtus
finish
|