搜索ANSYS熱分析指南與案例分析（三）
2017-05-08 by:CAE仿真在線 來源:互聯(lián)網(wǎng)
第四章 瞬態(tài)傳熱分析
一、瞬態(tài)傳熱分析的定義
瞬態(tài)熱分析用于計算一個系統(tǒng)的隨時間變化的溫度場及其它熱參數(shù)。在工程上一般用瞬態(tài)熱分析計算溫度場,并將之作為熱載荷進行應(yīng)力分析。
瞬態(tài)熱分析的基本步驟與穩(wěn)態(tài)熱分析類似。主要的區(qū)別是瞬態(tài)熱分析中的載荷是隨時間變化的。為了表達隨時間變化的載荷,首先必須將載荷~時間曲線分為載荷步。載荷~時間曲線中的每一個拐點為一個載荷步,如下圖所示。
對于每一個載荷步,必須定義載荷值及時間值,同時必須選擇載荷步為漸變或階越。
二、瞬態(tài)熱分析中的單元及命令
瞬態(tài)熱分析中使用的單元與穩(wěn)態(tài)熱分析相同。要了解每個單元的詳細說明,請參閱《ANSYS Element Reference Guide》。要了解每個命令的詳細功能,請參閱《ANSYS Command Reference Guide》。
三、ANSYS 瞬態(tài)熱分析的主要步驟
· 建模
· 加載求解
· 后處理
四、建模
· 確定jobname、title、units, 進入PREP7;
· 定義單元類型并設(shè)置選項;
· 如果需要,定義單元實常數(shù);
· 定義材料熱性能:一般瞬態(tài)熱分析要定義導(dǎo)熱系數(shù)、密度及比熱;
· 建立幾何模型;
· 對幾何模型劃分網(wǎng)格。
關(guān)于建模及劃分網(wǎng)格,請參閱《ANSYS Modeling and Meshing Guide》。
五、加載求解
1、定義分析類型
· 如果第一次進行分析,或重新進行分析
GUI: Main Menu>Solution>Analysis Type>New Analysis>Transient
Command: ANTYPE,TRANSIENT,NEW
· 如果接著上次的分析繼續(xù)進行(例如增加其它載荷)
GUI: Main Menu>Solution>Analysis Type>Restart
Command: ANTYPE,TRANSIENT,REST
2、獲得瞬態(tài)熱分析的初始條件
①、定義均勻溫度場
如果已知模型的起始溫度是均勻的,可設(shè)定所有節(jié)點初始溫度
Command: TUNIF
GUI: Main Menu> Solution>-Loads->Settings>Uniform Temp
如果不在對話框中輸入數(shù)據(jù),則默認為參考溫度,參考溫度的值默認為零,但可通過如下方法設(shè)定參考溫度:
Command: TREF
GUI: Main Menu> Solution>-Loads->Settings>Reference Temp
注意:設(shè)定均勻的初始溫度,與如下的設(shè)定節(jié)點的溫度(自由度)不同
Command: D
GUI: Main Menu>Solution>-Loads->Apply>-Thermal->Temperature>On Nodes
初始均勻溫度僅對分析的第一個子步有效;而設(shè)定節(jié)點溫度將保持貫穿整個瞬態(tài)分析過程,除非通過下列方法刪除此約束:
Command: DDELE
GUI: Main Menu> Solution>-Loads->Delete>-Thermal-Temperature>On Nodes
②、設(shè)定非均勻的初始溫度
在瞬態(tài)熱分析中,節(jié)點溫度可以設(shè)定為不同的值:
Command: IC
GUI: Main Menu> Solution>Loads>Apply>-Initial Condit'n>Define
如果初始溫度場是不均勻的且又是未知的,就必須首先作穩(wěn)態(tài)熱分析確定初始條件:
· 設(shè)定載荷(如已知的溫度、熱對流等)
· 將時間積分設(shè)置為OFF:
Command: TIMINT, OFF
GUI: Main Menu> Preprocessor>Loads>-Load Step Opts-Time/Frequenc>Time Integration
· 設(shè)定一個只有一個子步的,時間很小的載荷步(例如0.001):
Command: TIME
GUI: Main Menu> Preprocessor>Loads>-Load Step Opts-Time/Frequenc>Time and Substps
· 寫入載荷步文件:
Command: LSWRITE
GUI: Main Menu> Preprocessor>Loads>Write LS File
或先求解:
Command: SOLVE
GUI: Main Menu> Solution>Solve>Current LS
注意:在第二載荷步中,要刪去所有設(shè)定的溫度,除非這些節(jié)點的溫度在瞬態(tài)分析與穩(wěn)態(tài)分析相同。
3、設(shè)定載荷步選項
①、普通選項
l 時間:本選項設(shè)定每一載荷步結(jié)束時的時間:
Command: TIME
GUI: Main Menu> Solution>-Load Step Opts-Time/Frequenc>Time and Substps
l 每個載荷步的載荷子步數(shù),或時間增量
對于非線性分析,每個載荷步需要多個載荷子步。時間步長的大小關(guān)系到計算的精度。步長越小,計算精度越高,同時計算的時間越長。根據(jù)線性傳導(dǎo)熱傳遞,可以按如下公式估計初始時間步長:
其中
為沿熱流方向熱梯度最大處的單元的長度,
為導(dǎo)溫系數(shù),它等于導(dǎo)熱系數(shù)除以密度與比熱的乘積(
)。
Command: NSUBST or DELTIM
GUI: Main Menu> Solution>-Load Step Opts->Time/Frequenc>Time and Substps
如果載荷在這個載荷步是恒定的,需要設(shè)為階越選項;如果載荷值隨時間線性變化,則要設(shè)定為漸變選項:
Command: KBC
GUI: Main Menu> Solution>-Load Step Opts->Time/Frequenc>Time and Substps
②、非線性選項
l 迭代次數(shù):每個子步默認的次數(shù)為25,這對大多數(shù)非線性熱分析已經(jīng)足夠。
Command: NEQIT
GUI: Main Menu> Solution>-Load step opts>Nonlinear>Equilibrium Iter
l 自動時間步長:本選項為ON時,在求解過程中將自動調(diào)整時間步長。
Command: AUTOTS
GUI: Main Menu> Solution>-Load Step Opts->Time/Frequenc>Time and Substps
l 時間積分效果:如果將此選項設(shè)定為OFF,將進行穩(wěn)態(tài)熱分析。
Command: TIMINT
GUI: Main Menu> Solution>-Load Step Opts->Time/Frequenc>Time Integration
③、輸出選項
· 控制打印輸出:本選項可將任何結(jié)果數(shù)據(jù)輸出到*.out 文件中
Command: OUTPR
GUI: Main Menu>Solution>-Load Step Opts->Output Ctrls>Solu Printout
· 控制結(jié)果文件:控制*.rth的內(nèi)容
Command: OUTRES
GUI: Main Menu>Solution>-Load Step Opts->Output Ctrls>DB/Results File
?、存盤求解
六、后處理
ANSYS提供兩種后處理方式:
· POST1,可以對整個模型在某一載荷步(時間點)的結(jié)果進行后處理;
Command: POST1
GUI: Main Menu>General Postproc.
· POST26,可以對模型中特定點在所有載荷步(整個瞬態(tài)過程)的結(jié)果進行后處理。
Command: POST26
GUI: Main Menu>TimeHist Postproc
1、用POST1進行后處理
· 進入POST1后,可以讀出某一時間點的結(jié)果:
Command: SET
GUI: Main Menu>General Postproc>Read Results>By Time/Freq
如果設(shè)定的時間點不在任何一個子步的時間點上,ANSYS會進行線性插值。
· 此外還可以讀出某一載荷步的結(jié)果:
GUI: Main Menu>General Postproc>Read Results>By Load Step
然后就可以采用與穩(wěn)態(tài)熱分析類似的方法,對結(jié)果進行彩色云圖顯示、矢量圖顯示、打印列表等后處理。
2、用POST26進行后處理
· 首先要定義變量:
Command: NSOL or ESOL or RFORCE
GUI: Main Menu>TimeHist Postproc>Define Variables
· 然后就可以繪制這些變量隨時間變化的曲線:
Command: PLVAR
GUI: Main Menu>TimeHist Postproc>Graph Variables
或列表輸出:
Command: PRVAR
GUI: Main Menu>TimeHist Postproc>List Variables
此外,POST26還提供許多其它功能,如對變量進行數(shù)學操作等,請參閱《ANSYS Basic Analysis Procedures Guide》
七、相變問題
ANSYS熱分析最強大的功能之一就是可以分析相變問題,例如凝固或熔化等。含有相變問題的熱分析是一個非線性的瞬態(tài)的問題:
· 相變問題需要考慮熔融潛熱,即在相變過程吸收或釋放的熱量。ANSYS通過定義材料的焓隨溫度變化來考慮熔融潛熱(如圖所示)。
焓的單位是J/m3,是密度與比熱的乘積對溫度的積分:
求解相變問題,應(yīng)當設(shè)定足夠小的時間步長,并將自動時間步長設(shè)置為ON;
· 選用低階的熱單元,例如PLANE55或SOLID70。如果必須選用高階單元,請將單元選項KEYOPT(1)設(shè)置為1:
Command: keyopt(1)=1
GUI: Main Menu>Prepocessor>Element Type>Add/Edit/Delete>Options>-Specific heat matrix->Diagonalized
· 在設(shè)定瞬態(tài)積分參數(shù)時,請將THETA值設(shè)置為1(默認為0.5):
Command: TINTP
GUI: Main Menu> Solution>-Load and Step Opts->Time/Frequence>Time intergration>THETA
· 線性搜索將有助于加速相變問題的求解。
Command: LNSRCH
GUI: Main Menu>Solution>-Load and Step Opts->Nonlinear>Line Search
實例1:
一個30公斤重、溫度為70℃的銅塊,以及一個20公斤重、溫度為80℃的鐵塊,突然放入溫度為20℃、盛滿了300升水的、完全絕熱的水箱中,如圖所示。過了一個小時,求銅塊與鐵塊的最高溫度(假設(shè)忽略水的流動)。
材料熱物理性能如下:
|
熱性能 |
單位制 |
銅 |
鐵 |
水 |
|
導(dǎo)熱系數(shù) |
W/m℃ |
383 |
70 |
.61 |
|
密度 |
Kg/m3 |
8889 |
7833 |
996 |
|
比熱 |
J/kg℃ |
390 |
448 |
4185 |
以下列出log文件及菜單操作說明
/filename,transient1
/title, Thermal Transient Exercise 1
!進入前處理
/prep7
et,1,plane77 ! 定義單元類型
mp,kxx,1,383 ! 定義材料熱性能參數(shù)
mp,dens,1,8889 !1~銅,2~鐵,3~水
mp,c,1,390
mp,kxx,2,70
mp,dens,2,7837
mp,c,2,448
mp,kxx,3,0.61
mp,dens,3,996
mp,c,3,4185
rectnag,0,0.6,0,0.5 !創(chuàng)建幾何實體
rectang,0.15,0.225,0.225,0.27
rectang,0.6-0.2-0.058,0.6-0.2,0.225,0.225+0.044
aovlap,all !布爾操作
/pnum,area,1
aplot
aatt,1,1,1 !劃分網(wǎng)格
eshape,2
esize,0.02
amesh,2
aatt,2,1,1
amesh,3
aatt,3,1,1
eshape,3
esize,0.05
amesh,4
/pnum,mat,1
eplot
finish
!加載求解
/solu
antype,trans
timint,off !先作穩(wěn)態(tài)分析,確定初始條件
time,0.01 !設(shè)定只有一個子步的時間很小的載荷步
deltim,0.01
esel,s,mat,,3
nsle,s
d,all,temp,20
esel,s,mat,,2
nsle,s
d,all,temp,80
esel,s,mat,,1
nsle,s
d,all,temp,70
allsel
solve !得到初始溫度分布
time,3600 !進行瞬態(tài)分析
timint,on !打開時間積分
deltim,26,2,200 !設(shè)置時間步長,最大及最小時間步長
autots,on !打開自動時間步長
ddelet,all,temp !刪除穩(wěn)態(tài)分析中定義的節(jié)點溫度
outres,all,1 !將每個子步的值寫入數(shù)據(jù)庫文件
solve
finish
save
!進入POST26后處理
/post26
solu,2,dtime,,dtime !2~每一子步采用的時間步長
nsol,3,node(0.1875,0.2475,0),temp,,T_Copper !3~銅塊的中心點
nsol,4,node(0.371,0.247,0),temp,,T_Iron !4~鐵塊的中心點
nsol,5,node(30,0,0),temp,,T_H2O_Bot !5~水箱的底部
nsol,6,node(30,50,0),temp,,T_H2O_Top !6~水箱的頂部
nsol,7,node(0,25,0),temp,,T_H2O_Left !7~水箱的左部
nsol,8,node(60,25,0),temp,,T_H2O_Right !8~水箱的右部
Plvar,2
plvar,3,4,5,6,7,8
finish
!進入POST1后處理
/post1 !設(shè)置為最后一個載荷子步
set,last
esel,s,mat,,1
nsle,s
plnsol,temp
esel,s,mat,,2
nsle,s
plnsol,temp
finish
菜單操作:
1、Utility Menu: File>Change Jobname, 輸入文件名Transient1;
2、Utility Menu: File>Change Title,輸入Thermal Transient Exercise 1;
3、Main Menu>Preprocessor,進入前處理;
4、Main Menu>Preprocessor>Element Type>Add/Edit/Delete, 選擇單元plane77;
5、Main Menu>Preprocessor>Material Props>Isotropic, 定義材料1(銅)的KXX等于 383、DENS等于8889、C等于1390;
6、Main Menu>Preprocessor>Material Props>Isotropic, 定義材料2(鐵)的KXX等于 70、DENS等于7837、C等于448;
7、Main Menu>Preprocessor>Material Props>Isotropic, 定義材料3(水)的KXX等于0.61、DENS等于996、C等于4185;
8、Main Menu>Preprocessor>-Modeling->Create>-Areas->Retangle>By Dimensions, 輸入X1=0, Y1=0, X2=0.6, Y2=0.5, 選擇Apply, 輸入X1=0.15, Y1=0.225, X2= 0.225, Y2=0.27, 選擇Apply, 輸入X1=0.6-0.2-0.058, Y1=0.225, X2=0.6-0.2, Y2=0.225+0.044, 選擇OK;
9、Main Menu>Preprocessor>-Modeling->Operate>Booleans>Overlap, 選擇Pick All;
10、Utility Menu: Plotctrls>Numbering>Areas, on;
11、Utility Menu:Plot>Areas;
12、Main Menu>Preprocessor>-Attributes->Define->All Areas, 選擇材料1;
13、Main Menu>Preprocessor>Meshing->Size Cntrls->-Manualsize->-Global->Size, 輸入單元大小0.02;
14、Main Menu>Preprocessor>Meshing->Mesh->-Areas->Mapped>3 or 4 sided, 選擇銅塊;
15、Main Menu>Preprocessor>-Attributes->Define->All Areas, 選擇材料2;
16、Main Menu>Preprocessor>Meshing->Mesh->-Areas->Mapped>3 or 4 sided, 選擇鐵塊;
17、Main Menu>Preprocessor>-Attributes->Define->All Areas, 選擇材料3;
18、Main Menu>Preprocessor>Meshing->Size Cntrls->-Manualsize->-Global->Size, 輸入單元大小0.05;
19、Main Menu>Preprocessor>Meshing->Mesh->-Areas->Free, 選擇水箱;
20、Utility Menu>Plot>Area;
21、Main Menu>Solution,進入加載求解;
22、Main Menu>Solution>-Analysis Type->New Analysis, 選擇Transient,定義為瞬態(tài)分析;
23、Main Menu>Solution>-Load Step Opts>Time/Frenquenc>Time Integration, TIMINT, off,首先進行穩(wěn)態(tài)分析;
24、Main Menu>Solution>-Load Step Opts>Time/Frenquenc>Time-Time Step,設(shè)定TIME為0.01、DELTIM也為0.01;
25、Utility Menu: Select>Element>mat, 輸入3, 選擇Apply, 選擇Nodes>Attached to, 選擇Element;
26、Main Menu>Solution>-Loads->Apply>-Thermal->Temperature>On Nodes, 選擇Pick All, 輸入20;
27、Utility Menu: Select>Element>mat, 輸入2, 選擇Apply, 選擇Nodes>Attached to, 選擇Element;
28、Main Menu>Solution>-Loads->Apply>-Thermal->Temperature>On Nodes, 選擇Pick All, 輸入80;
29、Utility Menu: Select>Element>mat, 輸入1, 選擇Apply, 選擇Nodes>Attached to, 選擇Element;
30、Main Menu>Solution>-Loads->Apply>-Thermal->Temperature>On Nodes, 選擇Pick All, 輸入70;
31、Utility Menu: Select Everything
32、Main Menu>Solution>-Solve->Current LS
33、Main Menu>Solution>-Load Step Opts>Time/Frenquenc>Time-Time Step,設(shè)定TIME=3600, DELTIM=26, 最小、最大時間步長分別為2, 200, 將Autots設(shè)置為ON;
34、Main Menu>Solution>-Load Step Opts>Time/Frenquenc>Time Integration, 將TIMINT設(shè)置為ON;
35、Main Menu>Solution>-Loads->Delete>-Thermal->Temperature>On Nodes, 選擇Pick All,刪除穩(wěn)態(tài)分析定義的節(jié)點溫度;
36、Main Menu>Solution>-Load Step Opts>Output Ctrls->DB/Results, 選擇Every Substeps;
37、Main Menu>Solution>-Solve->Current LS;
38、Main Menu>TimeHist PostPro,進入POST26;
39、Main Menu>TimeHist PostPro>Define Variables>Add>Solution summary,在User specified label框中輸入dtime,選擇Solution Items>Step Time,選擇OK,Add>Nodal result>,在Node number框中輸入node(0.1875,0.2475,0),在User specified label框中輸入T_Copper,同樣輸入其它節(jié)點;
40、Main Menu>TimeHist PostPro>Graph Virables, 輸入變量代號,顯示各變量隨時間變化的曲線;
41、Main Menu>General Postproc,進入POST1;
42、Main Menu>General Postproc>-Read Results->Last set
43、Utility Menu>Select>Element>mat,輸入1,選擇Apply,Nodes>Attach to,選擇Element, OK;
44、Main Menu>General Postproc>Plot result>Nodal Solution,選擇temperature;
45、Utility Menu>Select>Element>mat,輸入2,選擇Apply,Nodes>Attach to,選擇Element, OK;
46、Main Menu>General Postproc>Plot result>Nodal Solution,選擇temperature;
實例2:
一鋼鑄件及其砂模的橫截面尺寸如圖所示:
砂模的熱物理性能如下表所示:
|
|
單位制 |
|
|
導(dǎo)熱系數(shù)(KXX) |
Btu/hr.in.oF |
0.025 |
|
密度(DENS) |
lbm/in3 |
0.254 |
|
比熱(C) |
Btu/lbm.oF |
0.28 |
鑄鋼的熱物理性能如下表所示:
|
|
單位制 |
0oF |
2643oF |
2750oF |
2875oF |
|
導(dǎo)熱系數(shù) |
Btu/hr.in.oF |
1.44 |
1.54 |
1.22 |
1.22 |
|
焓 |
Btu/in3 |
0 |
128.1 |
163.8 |
174.2 |
初始條件:鑄鋼的溫度為2875oF,砂模的溫度為80oF;
砂模外邊界的對流邊界條件:對流系數(shù)0.014Btu/hr.in2.oF,空氣溫度80oF;
求3個小時后鑄鋼及砂模的溫度分布。
/Title, Casting Solidification
!進入前處理
/prep7
et,1,plane55 !定義單元
mp,dens,1,0.054 !定義砂模熱性能
mp,kxx,1,0.025
mp,c,1,0.28
mptemp,1,0,2643,2750,2875 !定義鑄鋼的熱性能
mpdata,kxx,2,1.44,1.54,1.22,1.22
mpdata,enth,2,0,128.1,163.8,174.2
mpplot,kxx,2
mpplot,enth,2
save
!創(chuàng)建幾何模型
k,1,0,0,0
k,2,22,0,0
k,3,10,12,0
k,4,0,12,0
/pnum,kp,1
/pnum,line,1
/pnum,area,1
/Triad,ltop
kplot
a,1,2,3,4
save
rectng,4,22,4,8
aplot
aovlap,all
adele,3
aplot
save
!劃分網(wǎng)格
esize,1
amesh,5
mat,2
aplot
amesh,4
eplot
/pnum,elem
/number,1
save
!進入加載求解
/SOLU
antype,trans !設(shè)定為瞬態(tài)分析
esel,s,mat,,2 !設(shè)定鑄鋼的初始溫度
nsle,s
/replot
ic,all,temp,2875
esel,inve !設(shè)定砂模的初始溫度
nsle,s
/replot
ic,all,temp,80
allsel
save
lplot
sfl,1,CONV,0.014,,80 !設(shè)定砂模外邊界對流
sfl,3,CONV,0.014,,80
sfl,4,CONV,0.014,,80
/psf,conv,2
time,3 !設(shè)定瞬態(tài)分析時間
kbc,1 !設(shè)定為階越的載荷
autots,on !打開自動時間步長
deltim,0.01,0.001,0.25 !設(shè)定時間步長
timint,on !打開時間積分
tintp,,,,1 !將THETA設(shè)定為1
outres,all,all !輸入每個子步的結(jié)果
solve
!進入后處理
/post26
/pnum,node,1
/number,0
eplot
nsol,2,204,temp,center !設(shè)定鑄鋼中心點溫度隨時間的變量
plvar,2 !繪制溫度~時間曲線
save
finish
菜單操作:
1.Utility Menu>File>Change Title, 輸入Casting Solidification;
2.定義單元類型:Main Menu>Preprocessor>Element Type>Add/Edit/Delete, Add, Quad 4node 55;
3.定義砂模熱性能:Main Menu>Preprocessor>Material Props>Isotropic, 默認材料編號1, 在Density(DENS)框中輸入0.054, 在Thermal conductivity (KXX)框中輸入0.025, 在Specific heat(C)框中輸入0.28;
4.定義鑄鋼熱性能溫度表:Main Menu>Preprocessor>Material Props>-Temp Dependent->Temp Table, 輸入T1=0,T2=2643, T3=2750, T4=2875;
5.定義鑄鋼熱性能:Main Menu>Preprocessor>Material Props>-Temp Dependent ->Prop Table, 選擇Th Conductivity,選擇KXX, 輸入材料編號2,輸入C1=1.44, C2=1.54, C3=1.22, C4=1.22,選擇Apply, 選擇Enthalpy,輸入C1=0, C2=128.1, C3=163.8, C4=174.2;
6.創(chuàng)建關(guān)鍵點:Main Menu>Preprocessor>-Modeling->Create>Keypoints>In Active CS,輸入關(guān)鍵點編號1,輸入坐標0,0,0, 輸入關(guān)鍵點編號2, 輸入坐標22,0,0, 輸入關(guān)鍵點編號3, 輸入坐標10,12,0, 輸入關(guān)鍵點編號4, 輸入坐標0,12,0;
7.創(chuàng)建幾何模型:Main Menu>Preprocessor>-Modeling->Create>-Areas->Arbitrary>Through KPs,順序選取關(guān)鍵點1,2,3,4;
8.Main Menu>Preprocessor>-Modeling->Create>-Areas->Rectangle>By Dimension,輸入X1=4,X2=22,Y1=4,Y2=8;
9.進行布爾操作:
10.Main Menu>Preprocessor>-Modeling->Operate>-Booleans->Overlap>Area,Pick all;
11.刪除多余面:Main Menu>Preprocessor>-Modeling->Delete>Area and Below,3
12.保存數(shù)據(jù)庫:在Ansys Toolbar中選取SAVE_DB;
13.定義單元大小:Main Menu>Preprocessor>-Meshing->Size Cntrls>-Global->Size, 在Element edge length框中輸入1;
14.對砂模劃分網(wǎng)格:Main Menu>Preprocessor>-Meshing->Mesh>-Areas->Free,選擇砂模;
15.對鑄鋼劃分網(wǎng)格:Main Menu>Preprocessor>-Attributes->Define>Default Attribs, 在Material number菜單中選擇2;
16.Main Menu>Preprocessor>-Meshing->Mesh>-Areas->Free,選擇鑄鋼;
17.定義分析類型:Main Menu>Solution>-Analysis Type->New Analysis, 選擇Transient;
18.選擇鑄鋼上的節(jié)點:Utility Menu>Select>Entities, 選擇element,mat,輸入2,選擇Apply,選擇node, attached to element,選擇OK;
19.定義鑄鋼的初始溫度:Main Menu>Solution>-Loads->Apply>Initial Condit’n>Define, 選擇Pick all,選擇temp, 輸入2875, OK;
20.選擇砂模上的節(jié)點:Utility Menu>Select>Entities,Nodes, inverse
21.定義砂模的初始溫度:Main Menu>Solution>-Loads->Apply>Initial Condit’n>Define, 選擇Pick all, 選擇temp, 輸入80, OK;
22.Utility Menu>Select>Everything;
23.Utility Menu>Plot>Lines;
24.定義對流邊界條件: Main Menu>Solution>-Loads->Apply>-Thermal->Converction>On Lines,選擇砂模的三個邊界1,3,4, 在file coefficent框中輸入80, 在Bulk temperature框中輸入, 80;
25.設(shè)定瞬態(tài)分析時間選項:
Main Menu>Solution>Load Step Opts>Time/Frequenc>Time-Time Step,
Time at end of load step 3
Time Step size 0.01
Stepped or ramped b.c. Stepped
Automatic time stepping on
Minimun time Step size 0.001
Maximum time step size 0.25
26.設(shè)置輸出:Main Menu>Solution>Load Step Opts>Output Ctrls>DB/Results File, 在File write frequency框中選擇Every substep;
27.求解:Main Menu>Solution>-Solve->Current LS;
28.進入后處理: Main Menu>Timehist Postproc;
29.定義鑄鋼中心節(jié)點的溫度變量: Main Menu>Timehist Postproc>Define Variables, Add, Nodal DOF result,2,204;
30.繪制節(jié)點溫度隨時間變化曲線:Main Menu>Timehist Postproc>Graph Variable,2。
《ANSYS Verification Manual》中關(guān)于瞬態(tài)熱分析的實例:
VM104 Liquid-solid phase change
VM109 Temperature gradient across a solid cylinder
VM110 Transient temperature distribution in a slab
VM111 Cooling of a spherical body
VM112 Cooling of a spherical body
VM113 Transient temperature distribution in an orthotropic metal bar
VM114 Temperature response to a linearly rising surface temperature
VM115 Thermal response of heat generating slab
VM116 Heat conducting plate with sudden cooling
第五章 熱輻射
一、什么是熱輻射
輻射是一種通過電磁波傳遞能量的方式。電磁波以光速傳播且無需任何介質(zhì)。熱輻射僅為電磁波譜
中的一小段。因為由于熱輻射引起的熱流與物體表面絕對溫度的四次方成正比,因此熱輻射分析是
高度非線性的。
二、分析熱輻射問題
ANSYS提供了三 種方法分析熱輻射問題:
· 用LINK31,輻射線單元,分析兩個點或多對點之間的熱輻射;
· 用表面效應(yīng)單元SURF19或SURF22,分析點對面的熱輻射;
· 用AUX12,熱輻射矩陣生成器,分析面與面之間的熱輻射
以上三種方法既可用于穩(wěn)態(tài)熱分析,也可用于瞬態(tài)熱分析。
熱輻射分析要注意溫度的單位制,因為計算熱輻射使用的溫度單位是絕對溫度。如果在加載時使用
的是華氏溫度,就要設(shè)置460的差值;如果為攝氏溫度,差值為273。
Command:TOFFST
GUI: Main Menu>Preprocessor>Loads>Analysis Options
GUI: Main Menu>Solution>Analysis Options
三、使用LINK31—輻射線單元
LINK31是一個兩節(jié)點非線性線單元,用于計算由輻射引起的兩點之間的熱傳遞。此單元要求輸入如下的實常數(shù):
· 有效的熱輻射面積;
· 形狀系數(shù)
· 輻射率
· Stefan-Boltzmann 常數(shù)
四、使用表面效應(yīng)單元
表面效應(yīng)單元可以方便地分析點與面之間的輻射傳熱。SURF19用于兩維模型,SURF22用于三維模
型。單元應(yīng)設(shè)置為包含輻射KEYOPT(9)。
五、使用AUX12—輻射矩陣生成器
此方法用于計算多個輻射面之間的輻射傳熱。這種方法生成輻射面之間形狀系數(shù)矩陣,并將此矩陣
作為超單元用于熱分析。
AUX12方法由三個步驟組成:
· 定義輻射面
· 生成輻射矩陣
· 在熱分析中使用輻射矩陣
1、定義輻射面
￡¨1￡©在PREP7中建模、劃分網(wǎng)格。輻射面往往是3D模型中的面或2D模型中的邊,如下圖所示:
￡¨2￡©在輻射表面用SHELL57(3D)或LINK32(2D)劃分網(wǎng)格。最好的方法是先選擇輻射表面的
節(jié)點,然后用如下方法創(chuàng)建SHELL57或LINK32單元:
Command: ESURF
GUI: Main Menu>Preprocessor>Create>Elements>On Free Surf
注意:輻射面上的SHELL57或LINK32單元與節(jié)點必須與實體單元相吻合,否則計算的結(jié)果是不正確的。
生成的SHELL57或LINK32單元的取向也很重要。AUX12假定輻射的方向是SHELL57的+Z向或LINK32的+Y向。因此在生成SHELL57或LJNK32單元時要注意節(jié)點的排列順序。如圖所示:
(3)、如果所分析的系統(tǒng)是開放的,即一個面所輻射的熱能未被模型中其它的面吸收,則必須定義一個空間節(jié)點,用于吸收損失的輻射熱量。這個節(jié)點的位置是任意的。對于封閉的系統(tǒng),不應(yīng)定義
空間節(jié)點。
2、生成輻射矩陣
(1)、進入AUX12:
Command: /AUX12
GUI: Main Menu>Radiation Matrix
(2)、選擇組成輻射面的節(jié)點和單元。比較方便的方法是根據(jù)單元類型選擇單元,并選擇單元上的節(jié)點。有關(guān)選擇目標的技巧,請參閱《ANSYS Basic Analysis Guide》
(3)、確定模型是3D還是2D:
Command: GEOM
GUI: Main Menu>Radiation Matrix>Other Setting
AUX12用不同的算法計算2D或3D模型的形狀系數(shù)。AUX12默認為3D。2D分為純平面或軸對稱,默認為純平面。
(4)定義每個輻射面的輻射率(默認為1):
Command: EMIS
GUI: Main Menu>Radiation Matrix>Emissivities
(5)、定義Stefan-Boltzmann常數(shù)(默認為英制單位0.199E-10Btu/hr-in2-R4)
Command: STEF
GUI: Main Menu>Radiation Matrix>Other Settings
(6)、確定用什么方式計算形狀系數(shù):
Command: VTYPE
GUI: Main Menu>Radiation Matrix>Write Matrix
選擇是隱藏還是非隱藏方法:
· 非隱藏方法計算每兩個單元之間的形狀系數(shù),無論它們之間有無障礙;
· 隱藏方法(默認)用一種隱藏線算法判斷兩輻射面之間是否“可見”,如果可見則計算形狀系數(shù)。
(7)如為開放系統(tǒng),定義空間節(jié)點:
Command: SPACE
GUI: Main Menu>Radiation Matrix>Other Settings
(8)計算輻射矩陣并寫入jobename.sub文件:
Command: WRITE
GUI: Main Menu>Radiation Matrix>Write Matrix
如果要打印此矩陣,在執(zhí)行上述命令之前輸入MPRINT,1。
(9)選擇所有的節(jié)點和單元
3、在熱分析中使用輻射矩陣
(1)、重新進入PREP7,定義一個新的單元類型MATRIX50(超單元);
(2)、將單元類型指向超單元:
Command: TYPE
GUI: Main Menu>Preprocessor>-Modeling-Create>Elements>Elem Attributes
(3)、讀入超單元矩陣:
Command: SE
GUI: Main Menu>Preprocessor>-Modeling-Create>Elements>-Superelements-From .SUB File
(4)、不選擇或刪除用于生成輻射矩陣的SHELL57或LINK32單元:
Command: EDELE
GUI: Main Menu>Preprocessor>-Modeling-Delete>Elements
(5)、進入/SOLUTION施加其它熱載荷并求解。
(綜合自網(wǎng)絡(luò))
相關(guān)標簽搜索：ANSYS熱分析指南與案例分析（三） Ansys有限元培訓(xùn) Ansys workbench培訓(xùn) ansys視頻教程 ansys workbench教程 ansys APDL經(jīng)典教程 ansys資料下載 ansys技術(shù)咨詢 ansys基礎(chǔ)知識 ansys代做 Fluent、CFX流體分析 HFSS電磁分析 Abaqus培訓(xùn)
