Ansys apdl 结合生死单元仿真激光烧蚀纤维复合材料,表面温度过高? 200
考虑相变(直接汽化),手动添加了大概的焓值,表面基质相汽化温度3460°C,纤维相汽化温度4000°C,模型如图
加载激光热流密度,生死单元求解,网格最细划到了50等分(100等分不求解,不知道为啥)
热流密度大的话,基质和纤维烧了一部分,但温度最高达到了10000°C
热流密度小的话,只烧了基质,但表面温度最高还是有5000°C,按理来说高于汽化温度的不是应该都杀死了吗?
是我两种不同材料的生死单元设置的有问题吗?生死单元部分apdl代码如下
etable,temperature,temp !将单元温度存储在temperature表中
esel,s,MAT,,1 !选择材料1的单元
esel,r,etab,temperature,3460 !从材料1单元中选择温度大于汽化温度的单元
cm,vapor1,ELEM !创建基质烧蚀单元集1
esel,s,MAT,,2
esel,r,etab,temperature,4000
cm,vapor2,ELEM
etable,temperature,ERAS !删除temperature表
cmsel,all
finish !退出当前模块
/solu !进入求解模块
*get,elemnum,elem,0,count !读取单元总数,赋给参数elemnum
*if,elemnum,gt,0,then !如果单元数量大于0
ekill,all !杀死选择的单元
esel,s,live !选择活单元
*else !如果单元数量小于0(没有发生烧蚀)
esel,s,live !选择活单元
*endif
nsle,s !选择活单元上的节点
nsel,inve !选择死单元上的节点 D,all,all,0 !约束死节点自由度
nsel,inve !选择活单元上的节点