感应加热问题中的表面效应单元求助（命令流）！

无悔060961

以下是帮助文档中的感应加热例题
/batch,list
/filenam, induc
/prep7
shpp,off
/title, induction heating of a solid cylinder billet
/com,

!第一步： 选择磁场单元
et,1,13,,,1             ! PLANE13, axisymmetric, AZ dof
et,2,13,,,1
et,3,151,,,1,1,1        ! SURF151, thermal, radiation
r,3,0                   ! Real constant set for SURF151

!第二步： 输入参数
row=.015                ! 部件外径 outer radius of workpiece
ric=.0175               ! 先圈内径 inner radius of coil
roc=.0200               ! 线圈外径 outer radius of coil
ro=.05                  ! 模型外径 outer radius of model
t=.001                  ! 模型厚度 model thickness

freq=150000             ! 频率 frequency (Hz.)
pi=4*atan(1)            ! π pi
cond=.392e7             ! 最大电导率 maximum conductivity
muzero=4e-7*pi          ! 真空磁导率 free-space permeability
mur=200                 ! 最大相对磁导率 maximum relative permeability
skind=sqrt(1/(pi*freq*cond*muzero*mur)) ! 集肤深度 skin depth

ftime=3                   ! 结束时间 final time
tinc=.05                ! 谐波分析时间增量 time increment for harmonic analysis
time=0                  ! 起始时间 initialize time
delt=.01                ! 最大时间步 maximum delta time step

!第三步：定义材料特性
emunit,mks                  ! 设置磁场单位 set magnetic units
mp,murx,1,1             ! 空气相对磁导率 空气为一号材料 air relative permeability
mp,murx,3,1             ! 线圈相对磁导率 线圈为三号材料 coil relative permeability
mptemp,1,25.5,160,291.5,477.6,635,698   ! 钢的相对磁导率与温度的关系 temps for relative permeability
mptemp,7,709,720.3,742,761,1000
mpdata,murx,2,1,200,190,182,161,135,104 ! steel relative permeability
mpdata,murx,2,7,84,35,17,1,1
mptemp
mptemp,1,0,125,250,375,500,625          ! 钢的电阻率与温度的关系 temps for resistivity
mptemp,7,750,875,1000
mpdata,rsvx,2,1,.184e-6,.272e-6,.384e-6,.512e-6,.656e-6,.824e-6
mpdata,rsvx,2,7,1.032e-6,1.152e-6,1.2e-6 ! steel resistivity

!第四步：建立几何模型
rectng,0,row,0,t        ! 建立钢条billet
rectng,row,ric,0,t      ! 空气隙 air-gap
rectng,ric,roc,0,t      ! 线圈 coil
rectng,roc,ro,0,t       ! 远场空气 outer air
aglue,all               ! 布尔运算
numcmp,area             ! 压缩面编号

!第五步： 附特性，并剖分
ksel,s,loc,x,row        ! 选择一新的关于关键点或硬点的子集 select keypoints at outer radius of workpiece
kesize,all,skind/2      ! set meshing size to 1/2 skin depth
ksel,s,loc,x,0          ! 选择一新的关于关键点或硬点的子集 select keypoints at center
kesize,all,40*skind     ! 设置沿着关键点的单元尺寸set meshing size
lsel,s,loc,y,t/2        ! 选择一个新的关于线的子集 select vertical lines
lesize,all,,,1          ! set 1 division through thickness
lsel,all

asel,s,area,,1          !选择一新的关于面的子集
aatt,2,1,1              ! 赋予面积1为2号材料，1号单元  set attributes for billet region
asel,s,area,,3
aatt,3,1,2              ! 赋予面积3为3号材料2号单元set attributes for coil region
asel,s,area,,2,4,2
aatt,1,1,2              ! 赋予面积2,4为1号材料2号单元 set attributes for air region
asel,all
mshape,0,2d             !选择单元形状
mshk,1                   !确定划分网格所用的方式及自由网格还是映射 0表示自由划分
amesh,1                 ! 生成节点和在选定区域划分网格 mesh billet area
lsel,s,loc,y,0
lsel,a,loc,y,t
lsel,u,loc,x,row/2
lesize,all,.001
lsel,all
amesh,all               ! mesh remaining areas

!第六步：创造辐射单元
n                       ! 确定一节点 create space node for SURF151   
*get,nmax,node,,num,max !得到最大节点编号
lsel,s,loc,x,row
type,3
real,3
mat,2
lmesh,all            ! 生成节点和单元在选定的线上 mesh billet outer radius with SURF151      
*get,emax,elem,,num,max !得到最大单元编号
emodif,emax,3,nmax   ! 把单元修改成用于辐射的空间节点 modify element to add space node for radiation
et,3,0               ! 把3号单元转换成0号单元 reset type 3 to null element

!第七步：施加边界条件和载荷
nsel,s,loc,x
d,all,az,0              !施加磁力线平行边界条件 apply flux-normal b.c.
nsel,all
esel,s,mat,,3           !选择一个新的单元的子集
bfe,all,js,,,,15e6      ! 对线圈施加电流 密度 apply current density to coil
esel,all
finish

!第八步：定义求解选项，写入物理环境
/solu
antyp,harm
harfrq,150000
physics,write,emag      ! 写入 emag物理环境 write emag physics file
finish

!第九步：转入热分析，单元转换
/prep7
lsclear,all             ! 清除所有边界条件和选项 clear all b.c.'s and options
et,1,55,,,1             ! PLANE55 thermal element, axisymmetric
et,2,0                  ! null element type for coil and air region
et,3,151,,,1,1,1        ! SURF151 element for radiation
keyopt,3,9,1
r,3,1,5.67e-8           ! form factor, Stefan-Boltzmann constant

mptemp
mptemp,1,0,730,930,1000         ! temps for conductivity
mpdata,kxx,2,1,60.64,29.5,28,28
mptemp                          ! temps for enthalpy
mptemp,1,0,27,127,327,527,727
mptemp,7,765,765.001,927
mpdata,enth,2,1,0,91609056,453285756,1.2748e9,2.2519e9,3.3396e9
mpdata,enth,2,7,3.548547e9,3.548556e9,4.3520e9
mp,emis,2,.68                   ! 发射率 emissivity
finish

!第十步： 定义热分析求解选项，写入物理环境
/solu
antype,trans
toffst,273
tunif,100                       ! initial uniform temperature
d,nmax,temp,25                  ! ambient temperature
cnvtol,heat,1                   ! convergence tolerance
kbc,1                           ! step loads
trnopt,full
autos,on                        ! auto time-stepping
deltim,1e-5,1e-6,delt,on        ! time step control
outres,basic,all                ! save all load step information
physics,write,thermal           ! write thermal physics file
finish

!第十一步： 用命令流实现循环求解
*do,i,1,ftime/tinc   ! solution *do loop
time=time+tinc                             ! increment time

physics,read,emag    ! read emag physics file
/solu
*if,i,eq,1,then
tunif,100                    ! initial temperature
*else
ldread,temp,last,,,,,rth     ! read thermal analysis temperatures
*endif
solve                         ! solve harmonic analysis
finish
physics,read,thermal          ! read thermal physics file

/assign,esav,therm,esav       ! redirect files for use in thermal restart
/assign,emat,therm,emat
/solu
*if,i,gt,1,then
antype,trans,rest            ! thermal restart
*endif
time,time                     ! time at end of thermal run
esel,s,mat,,2                 ! select billet region
ldread,hgen,,,,2,,rmg         ! apply coupled joule heating load from emag
esel,all
solve
finish
/assign,esav                  ! reassign files to default
/assign,emat

*enddo                        ! end of solution looping
finish
save                          ! save database

!第十二步：后处理
/post26                       ! time-history postprocessor
/show
nsol,2,1,temp,,tempcl         ! store temperature at billet centerline
nsol,3,2,temp,,tempsurf       ! store temperature at billet outer diameter
plvar,2,3                     ! plot temperature rise over time
prvar,2,3                     ! print temperature rise over time
finish

但第6步创造辐射单元部分不是很懂！表面效应单元不是可以直接用于热辐射吗，为什么还要创造？既然创造了热辐射单元，为什么又要把3号单元转换成0号单元（et,3,0）！
如果是3维模型，选用surf152单元，在这个地方该怎样创造辐射单元？

zxc0228

楼主问题解决了吗？我也遇到了相同的问题