考虑相变（直接汽化），手动添加了大概的焓值，表面基质相汽化温度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         !选择活单元上的节点