ANSYS apdl命令流耦合场分析案例4--------三维感应加热

前言

感应加热的激励源为365000HZ的交流电，线圈电流密度为2.04e8A/m^2，线圈和管子的几何模型如下图所示，另外再在外面在加上一个空气域。
在进行求解时建议读者分两次求解，也就是分两次复制，先求解谐波磁场，这样就可以得到线圈与管子电流密度的分布矢量图，谐波磁场求解完成后输入:

/post1
set,1,1
plvect,jt,,,,vect,elem,on

就可以得到电流密度矢量图（如下图）：

完整感应加热命令

如果直接完整求解，电磁场的载荷和边界条件会被覆盖或者删除，最后只会得到瞬态热分析的结果。温度单位为开尔文，假定管道的初始温度为常温20摄氏度左右。

/prep7
et,1,solid236,1		!管
et,2,solid236  		!线圈
et,3,solid236		!空气
et,4,surf152,,,,1,1
r,4,0
mp,murx,1,1
mp,rsvx,1,1/2.326e7
mp,murx,2,1
mp,rsvx,2,2e-8
mp,murx,3,1
!===几何模型=============
cylind,15e-3,15.62e-3,0,60e-3
cylind,18e-3,20e-3,30e-3,55e-3
block,-0.04,0.04,-0.04,0.04,0.09,-0.03
vovlap,all
numcmp,all
!=====划分网格===============
type,1
mat,1
esize,0.002
mshape,1,3d
mshkey,0
vmesh,1
type,4
mat,1
real,4
!===创建热辐射区域=======
n
*get,nmax,node,,num,max
esurf,nmax
et,4,0
!======================
type,2
mat,2
esize,0.003
mshape,1,3d
vmesh,2
!=======================
type,3
mat,3
esize,0.015
mshape,1,3d
vmesh,3
!=====边界条件和载荷===========
allsel,all
asel,s,,,13,18
da,all,az,0
cswpla,14,1
vsel,s,,,2
eslv,s
emodif,all,esys,14
vsel,s,,,2
bfv,all,js,0,2.04e7,0
allsel,all
!=============================
/solu
antype,harm
harfrq,365e3
solve
finish
!======瞬态热分析================
/prep7
lsclear,all
et,1,solid90
et,2,0
et,3,0
et,4,surf152,,,,1,1
keyopt,4,9,1
r,4,1,5.67e-8
mp,kxx,1,155
mp,c,1,893
mp,dens,1,2730
mp,emis,1,0.3
d,nmax,temp,293
asel,s,,,3,6
sfa,all,,conv,30,293
esel,s,mat,,1
ldread,hgen,,,,2,,rst
!========================
/solu
antype,trans
toffst,0
tunif,293
cnvtol,heat,1
kbc,1
trnopt,full
time,3
deltim,1e-2
outres,basic,all
solve
finish
!====后处理============
/post1
set,1,last
plnsol,temp

结果温度云图如图所示：

总结

只要把空气域画成圆柱形，那么这个案例就是轴对称的，当然可以采用二维简化模型来建模，因为在ansys帮助文档里已经提供了 一个非线性材料的二维感应加热模型，它只是介绍了二维感应加热，单元的选取，边界条件的施加，所以这个案例只是为读者对三维感应加热单元及边界条件的应用提供参考。