問題說明
樣塊上下兩端固定，在室溫20℃下進行夾緊，分析其升溫到150℃時的熱應力。
采用間接法進行分析，溫度場單元選擇278，應力場單元為185
首先進行穩態熱分析，然后讀取溫度場的數據，進行結構分析，具體的編程代碼如下
溫度場分析APDL
/CLEAR
/FILENAME TEST
/PREP7
!定義溫度場實體單元
ET,1,SOLID278
!定義材料1的導熱系數
MPTEMP,,,,,,,,
MPTEMP,1,0
MPDATA,KXX,1,,20
!建立樣件模型
BLC5,0,0,0.2,0.1,0.05
!設置體1的網格尺寸
CM,_Y,VOLU
VSEL, , , , 1
CM,_Y1,VOLU
CMSEL,S,_Y
CMSEL,S,_Y1
VATT, 1, , 1, 0
CMSEL,S,_Y
CMDELE,_Y
CMDELE,_Y1
FLST,5,1,4,ORDE,1
FITEM,5,5
CM,_Y,LINE
LSEL, , , ,P51X
CM,_Y1,LINE
CMSEL,,_Y
LESIZE,_Y1, , ,20, , , , ,1
FLST,5,1,4,ORDE,1
FITEM,5,6
CM,_Y,LINE
LSEL, , , ,P51X
CM,_Y1,LINE
CMSEL,,_Y
LESIZE,_Y1, , ,10, , , , ,1
FLST,5,1,4,ORDE,1
FITEM,5,10
CM,_Y,LINE
LSEL, , , ,P51X
CM,_Y1,LINE
CMSEL,,_Y
LESIZE,_Y1, , ,5, , , , ,1
CM,_Y,VOLU
VSEL, , , , 1
CM,_Y1,VOLU
CHKMSH,'VOLU'
CMSEL,S,_Y
VCLEAR,_Y1
VMESH,_Y1
CMDELE,_Y
CMDELE,_Y1
CMDELE,_Y2
!在一側施加溫度150
FLST,2,1,5,ORDE,1
FITEM,2,5
/GO
DA,P51X,ALL,150
!設置其他表面的對流換熱系數
FLST,2,5,5,ORDE,3
FITEM,2,1
FITEM,2,-4
FITEM,2,6
/GO
SFA,P51X,1,CONV,50,
!設置分析類型為穩態求解
FINISH
/SOL
ANTYPE,0
!設置溫度偏移
TOFFST,273,
!寫溫度場物理分析文件
FINISH
/PREP7
PHYSICS,WRITE,Thermal1, , ,
!清空物理環境數據
PHYSICS,CLEAR
!刪除邊界和載荷
LSCLEAR,ALL
!轉換單元類型，thermal to struc
ETCHG,TTS
!設置材料參數
MPTEMP,,,,,,,,
MPTEMP,1,0
MPDATA,EX,1,,2.1E+011 !彈性模量
MPDATA,PRXY,1,,0.3 !泊松比
UIMP,1,REFT,,, !參考溫度
MPDATA,ALPX,1,,2.35E-5 !熱膨脹系數
!定義結構邊界條件--左右固定
FLST,2,2,5,ORDE,2
FITEM,2,5
FITEM,2,-6
/GO
DA,P51X,ALL,0
!設置參考溫度
TREF,20,
!求解設置
FINISH
/SOL
ANTYPE,0
NLGEOM,1
NSUBST,8,10,6
TIME,1
!寫結構物理分析文件
FINISH
/PREP7
PHYSICS,WRITE,struct1, , ,
!讀溫度場物理分析文件
PHYSICS,READ,THERMAL1
!求解
FINISH
/SOL
/STATUS,SOLU
SOLVE
!設置圖像顯示
/PLOPTS,INFO,1 !inf0為legend on
!顯示溫度雲圖
FINISH
/POST1
/EFACET,1
PLNSOL, TEMP,, 0
!讀取結構場物理分析文件
FINISH
/PREP7
PHYSICS,READ,STRUCT1
!定義溫度邊界條件
LDREAD,TEMP,,, , ,'test_temp1','rth',' '
!求解
FINISH
/SOL
/STATUS,SOLU
SOLVE
!查看結果
FINISH
/POST1
SET,LIST,999
SET,FIRST
/EFACET,1
PLNSOL, S,EQV, 0,1.0
結果
