iomega版主你好，我按照你给我提的建议改写了程序，请你看看是不是应该这样改

lchylinchy

版主，请问一下，我如果要在上表面同时加载热流密度和对流，是不是要用到表面热效应单元呢?
如果要用到，请问应该怎样用呢
主要是不知道怎样去编写命令流
谢谢版主

lchylinchy

版主，你好！
我现在考虑是把上平面同时加载对流载荷和热流密度；在另外三个面上加载对流载荷，用到了表面热效应单元，请你看看这样对吗？
!假设硅对激光的吸收率为0.9，激光照射直径为150μm
fini
/clear
/prep7
et,1,solid90
et,2,surf152             !定义表面效应单元
blc4,0,0,1.5e-3,0.04e-3,0.2e-3 !建立硅片模型长1.5mm，宽200μm，厚度40μm
mp,kxx,1,140              !定义硅片的导热系数
!mp,dens,1,2230              !定义硅片的密度

KEYOPT,2,3,0
KEYOPT,2,4,1
KEYOPT,2,5,0
KEYOPT,2,6,0
KEYOPT,2,7,0
KEYOPT,2,8,3          !根据表面温度计算对流系数
KEYOPT,2,9,0

!R,1
!R,2



lsel,s,,,2,4,2
lsel,a,,,5,7,2
lesize,all,,,100

lsel,s,,,9,12,3
lsel,a,,,10,11,1
lesize,all,,,20

lsel,s,,,3,6,3
lsel,a,,,1,8,7
lesize,all,,,4

type,1
vmesh,all

asel,s,,,4
nsla,s,1
type,2
esurf,all
alls
fini

/solu

!toffst,300
!ic,all,temp,298

asel,s,,,1,3,1
nsla,s,0
sf,all,conv,10,25          !其它三个面定义对流系数，alls
esel,s,type,,2
nsle,s
sf,all,conv,12.5,25     !上平面施加对流载荷


!alls

!施加一个直径为150μm的激光均匀热流密度
esel,s,type,,1
nsle,s
nsel,r,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
prnsol,TEMP
fini

iomega

R_conv  =  1/ (h*A), with unit of K/W
h = 10 W/m^2-K;
A: area for convection =   1.5mm (length)* [ 0.2mm (width) + 2* 0.04mm(thickness)]

iomega

you have to use surface elements if trying to apply both heat flux and convection on the same area.

In your new case, on top surface,only the center 150um square nodes will have heat flux, while the rest of the area nodes will have h - since applying heat flux operation overwrites the previous applying h process for the center 150um square nodes.

lchylinchy

谢谢楼主
前两天有事没有来论坛
谢谢你的帮助

lchylinchy

楼主
在我的分析中，要不要考虑到热辐射的情况呢？

lchylinchy

两位版主
请问一下，为什么我做的这个热分析中没有.rth文件呢？
在求解过程中，有一个警告是No constraints have been defined using the D command.
这个有什么问题吗

wq1627

我也是做激光材料去除，其中包含相变，为什么我加了焓值之后，似乎没有发生相变，温度到达熔点后没有下降或者保持的趋势。焓值好像并没有发生作用。
希望斑主能帮忙解答一下我的问题