考慮相變（直接汽化），手動(dòng)添加了大概的焓值，表面基質(zhì)相汽化溫度3460°C，纖維相汽化溫度4000°C，模型如圖
加載激光熱流密度，生死單元求解，網(wǎng)格最細(xì)劃到了50等分（100等分不求解，不知道為啥）
熱流密度大的話(huà)，基質(zhì)和纖維燒了一部分，但溫度最高達(dá)到了10000°C
熱流密度小的話(huà)，只燒了基質(zhì)，但表面溫度最高還是有5000°C，按理來(lái)說(shuō)高于汽化溫度的不是應(yīng)該都?xì)⑺懒藛幔?br>是我兩種不同材料的生死單元設(shè)置的有問(wèn)題嗎？生死單元部分apdl代碼如下

etable,temperature,temp   !將單元溫度存儲(chǔ)在temperature表中
esel,s,MAT,,1                                                                                                                !選擇材料1的單元
esel,r,etab,temperature,3460                                                !從材料1單元中選擇溫度大于汽化溫度的單元                                                                                             
cm,vapor1,ELEM                                                                                                        !創(chuàng)建基質(zhì)燒蝕單元集1

esel,s,MAT,,2
esel,r,etab,temperature,4000
cm,vapor2,ELEM
etable,temperature,ERAS                 !刪除temperature表
cmsel,all           
finish              !退出當(dāng)前模塊
/solu               !進(jìn)入求解模塊
*get,elemnum,elem,0,count      !讀取單元總數(shù)，賦給參數(shù)elemnum
                *if,elemnum,gt,0,then    !如果單元數(shù)量大于0
           ekill,all              !殺死選擇的單元
           esel,s,live          !選擇活單元
     *else                       !如果單元數(shù)量小于0(沒(méi)有發(fā)生燒蝕)
          esel,s,live           !選擇活單元
     *endif
nsle,s              !選擇活單元上的節(jié)點(diǎn)
nsel,inve         !選擇死單元上的節(jié)點(diǎn) D,all,all,0        !約束死節(jié)點(diǎn)自由度
nsel,inve         !選擇活單元上的節(jié)點(diǎn)