求助：关于热分析中没有RTH文件

lchylinchy

请问一下
我在同一个面上加载了热流密度，然后通过表面效应单元加载了对流
在分析过程中为什么没有RTH文件呢
为什么啊？

而且有一个警告提示： No constraints have been defined using the D command.

请问应该怎样去解决呢？
我的程序如下：
!假设硅对激光的吸收率为0.9，激光照射直径为150μm
fini
/clear
/filename,danliang
/title,thermal analysis of beam
/prep7
et,1,solid90
et,2,surf152
et,3,PLANE42

mp,kxx,1,140           !定义硅片的导热系数
!mp,dens,1,2230


KEYOPT,2,3,0           ! OMEG 用于全局笛卡儿坐标系X轴
KEYOPT,2,4,0           !有中间节点
KEYOPT,2,5,0           !关于额外节点的辐射或对流计算: 无额外节点
!KEYOPT,2,6,0          !使用体积温度,额外节点温度作为体积温度
KEYOPT,2,7,0           !对于膜层散热系数不增加经验项
KEYOPT,2,8,3           !根据表面温度计算对流系数
KEYOPT,2,9,0           !不包括热辐射



k,1,0,0
k,2,1.5e-3,0
k,3,1.5e-3,0.04e-3
k,4,0,0.04e-3              !定义硅片的密度

l,1,2
l,2,3
l,3,4
l,4,1
al,1,2,3,4

lsel,s,,,1,3,2
lesize,all,,,100

lsel,s,,,2,4,2
lesize,all,,,4

type,3
amesh,all

!通过拉伸生成体单元
type,1
extopt,ESIZE,20     !指定拉伸方向划分数为3
extopt,ACLEAR,1    !拉伸完毕后，清除面单元
vext,all,,,,,0.2e-3   !z方向拉伸0.3e-3

!生成表面效应单元
asel,s,,,5
nsla,s,0
type,2
esurf,all
alls
fini

/solu

!antype,static

toffst,273
allsel,all
ic,all,temp,25

!施加对流载荷
esel,s,type,,1
asel,s,,,1,6,1
asel,u,,,5
nsla,s,1
sf,all,conv,10,25          !定义对流系数，其中10表示对流系数，25表示
alls

esel,s,type,,2
asel,s,,,5
nsla,s,1
sf,all,conv,12.5,25          !定义对流系数，其中12.5表示对流系数，25表示



!施加一个直径为150μm的激光均匀热流密度
esel,s,type,,1
nsle,s
nsel,s,loc,x,0.6e-3,0.75e-3
nsel,r,loc,y,0.04e-3
nsel,r,loc,z,0.025e-3,0.175e-3
sf,all,hflux, 88108               !施加均匀热流密度
alls

solv
fini

/post1
plnsol,TEMP              !画出温度分布
/image,save,temp_dian,bmp
PLNSOL,TF,SUM            !画出热流量分布
/image,save,tf_dian,bmp
PLNSOL,TG,SUM            !画出热梯度分布
/image,save,tg_dian,bmp
prnsol,TEMP
fini

bqc

et,3,PLANE42 这里的单元属性定义有问题
‘PLANE42 is used for 2-D modeling of solid structures.’

更换et3为热分析使用单元比如55或77之后，就可以得到.rth结果文件了。

lchylinchy

哦
谢谢
非常感谢
如果我要讲我做的热分析结果用到热应力分析中去
我用顺序耦合的方式来做，请问应该怎样去定义耦合面呢？

woomii

如果是同一模型，用etchg切换单元，用ldread加载温度就可以了