1、熱分析,熱分析的目的 熱分析用于計算一個系統(tǒng)或部件的溫度分布及其它熱物理參數(shù)，如熱量的獲取或損失、熱梯度、熱流密度(熱通量)等 熱分析在許多工程應用中扮演重要角色，如內(nèi)燃機、渦輪機、換熱器、管路系統(tǒng)、電子元件等,ANSYS的熱分析,在ANSYS/Multiphysics、ANSYS/Mechanical、ANSYS/Thermal、ANSYS/FLOTRAN、ANSYS/ED五種產(chǎn)品中包含熱分析功能 ANSYS熱分析基于能量守恒原理的熱平衡方程，用有限元法計算各節(jié)點的溫度，并導出其它熱物理參數(shù) ANSYS熱分析包括熱傳導、熱對流及熱輻射三種熱傳遞方式。此外，還可以分析相變、有內(nèi)熱源、接觸熱阻等

2、問題,ANSYS的熱分析分類,ANSYS的熱分析分類 穩(wěn)態(tài)傳熱：系統(tǒng)的溫度場不隨時間變化 瞬態(tài)傳熱：系統(tǒng)的溫度場隨時間明顯變化 與熱有關的耦合分析 熱結構耦合 熱流體耦合 熱電耦合 熱磁耦合 熱電磁結構耦合等,熱分析的符號與單位,表征物體吸收的熱量，為一個體系的內(nèi)能與體系的體積和外界施加于體系的壓強的乘積之和,傳熱學經(jīng)典理論回顧,熱傳遞的方式,1、熱傳導 熱傳導可以定義為完全接觸的兩個物體之間或一個物體的不同部分之間由于溫度梯度而引起的內(nèi)能的交換。熱傳導遵循付里葉定律：qn=-k*(dT/dx)，式中qn為熱流密度（W/m2），k為導熱系數(shù)(W/m-)，“-”表示熱量流向溫度降低的方向。 2、

3、熱對流 熱對流是指固體的表面與它周圍接觸的流體之間，由于溫差的存在引起的熱量的交換。熱對流可以分為兩類：自然對流和強制對流。熱對流用牛頓冷卻方程來描述： qn= h*(TS-TB)，式中h為對流換熱系數(shù)(或稱膜傳熱系數(shù)、給熱系數(shù)、膜系數(shù)等)，TS為固體表面的溫度， TB為周圍流體的溫度。,熱傳遞的方式(續(xù)),3、熱輻射 熱輻射指物體發(fā)射電磁能，并被其它物體吸收轉(zhuǎn)變?yōu)闊岬臒崃拷粨Q過程。物體溫度越高，單位時間輻射的熱量越多。熱傳導和熱對流都需要有傳熱介質(zhì)，而熱輻射無須任何介質(zhì)。實質(zhì)上，在真空中的熱輻射效率最高。 在工程中通常考慮兩個或兩個以上物體之間的輻射，系統(tǒng)中每個物體同時輻射并吸收熱量。它們之

4、間的凈熱量傳遞可以用斯蒂芬波爾茲曼方程來計算：q=A1F12(T14-T24)，式中q為熱流率， 為輻射率（黑度）， 為斯蒂芬波爾茲曼常數(shù)，約為5.6710-8W/m2.K4，A1為輻射面1的面積，F(xiàn)12為由輻射面1到輻射面2的形狀系數(shù)，T1為輻射面1的絕對溫度，T2為輻射面2的絕對溫度。由上式可以看出，包含熱輻射的熱分析是高度非線性的。,穩(wěn)態(tài)傳熱,如果系統(tǒng)的凈熱流率為，即流入系統(tǒng)的熱量加上系統(tǒng)自身產(chǎn)生的熱量等于流出系統(tǒng)的熱量：q流入+q生成-q流出=0，則系統(tǒng)處于熱穩(wěn)態(tài)。在穩(wěn)態(tài)熱分析中任一節(jié)點的溫度不隨時間變化。穩(wěn)態(tài)熱分析的能量平衡方程為(以矩陣形式表示):KT=Q 式中： K為傳導矩陣，包

5、含導熱系數(shù)、對流系數(shù)及輻射率和形狀系數(shù)；T為節(jié)點溫度向量；Q為節(jié)點熱流率向量，包含熱生成； ANSYS利用模型幾何參數(shù)、材料熱性能參數(shù)以及所施加的邊界條件，生成K 、 T以及Q 。,瞬態(tài)傳熱,瞬態(tài)傳熱過程是指一個系統(tǒng)的加熱或冷卻過程。在這個過程中系統(tǒng)的溫度、熱流率、熱邊界條件以及系統(tǒng)內(nèi)能隨時間都有明顯變化。根據(jù)能量守恒原理，瞬態(tài)熱平衡可以表達為(以矩陣形式表示)：C +KT=Q 式中:K為傳導矩陣，包含導熱系數(shù)、對流系數(shù)及輻射率和形狀系數(shù)； C為比熱矩陣,考慮系統(tǒng)內(nèi)能的增加; T為節(jié)點溫度向量； 為溫度對時間的導數(shù); Q為節(jié)點熱流率向量，包含熱生成。,線性與非線性,如果有下列情況產(chǎn)生，則為非線

6、性熱分析： 材料熱性能隨溫度變化，如K(T),C(T)等； 邊界條件隨溫度變化，如h(T)等； 含有非線性單元； 考慮輻射傳熱 非線性熱分析的熱平衡矩陣方程為： C(T) +K(T)T=Q (T),邊界條件、初始條件,ANSYS熱分析的邊界條件或初始條件可分為七種： 溫度：模型區(qū)溫度已知 熱流率：熱流率已知的點 對流：表面的熱傳遞給周圍的流體通過對流。輸入對流換熱系數(shù)h和環(huán)境流體的 平均溫度Tb 熱輻射：通過輻射產(chǎn)生熱傳遞的面. 輸入輻射系數(shù)，Stefan-Boltzmann常數(shù)，“空間節(jié)點”的溫度作為可選項輸入 絕熱面： “完全絕熱”面，該面上不發(fā)生熱傳遞 熱通量：單位面積上的熱流率已知的面

7、 熱生成率：體的生熱率已知的區(qū)域,熱分析誤差估計,僅用于評估由于網(wǎng)格密度不夠帶來的誤差； 僅適用于SOLID或SHELL的熱單元(只有溫度一個自由度)； 基于單元邊界的熱流密度的不連續(xù)； 僅對一種材料、線性、穩(wěn)態(tài)熱分析有效； 使用自適應網(wǎng)格劃分可以對誤差進行控制。,穩(wěn)態(tài)傳熱分析,穩(wěn)態(tài)傳熱用于分析穩(wěn)定的熱載荷對系統(tǒng)或部件的影響。通常在進行瞬態(tài)熱分析以前，進行穩(wěn)態(tài)熱分析用于確定初始溫度分布。 穩(wěn)態(tài)熱分析可以通過有限元計算確定由于穩(wěn)定的熱載荷引起的溫度、熱梯度、熱流率、熱流密度等參數(shù),熱分析的單元,熱分析涉及到的單元有大約40種，其中純粹用于熱分析的有14種： 線性：LINK32兩維二節(jié)點熱傳導單元

8、 LINK33三維二節(jié)點熱傳導單元 LINK34二節(jié)點熱對流單元 LINK31二節(jié)點熱輻射單元 二維實體：PLANE55四節(jié)點四邊形單元 PLANE77八節(jié)點四邊形單元 PLANE35三節(jié)點三角形單元 PLANE75四節(jié)點軸對稱單元 PLANE78八節(jié)點軸對稱單元,熱分析的單元(續(xù)),三維實體：SOLID87六節(jié)點四面體單元 SOLID70八節(jié)點六面體單元 SOLID90二十節(jié)點六面體單元 殼：SHELL57四節(jié)點 點：MASS71,穩(wěn)態(tài)熱分析的基本過程,ANSYS熱分析可分為三個步驟： 前處理：建模 求解：施加載荷計算 后處理：查看結果,建模,確定jobname、title、unit； 進入

9、PREP7前處理，定義單元類型，設定單元選項； 定義單元實常數(shù)； 定義材料熱性能參數(shù)，對于穩(wěn)態(tài)傳熱，一般只需定義導熱系數(shù)，它可以是恒定的，也可以隨溫度變化； 創(chuàng)建幾何模型并劃分網(wǎng)格。,幾何尺寸(模型),既可用ANSYS建立模型，也可用其它方法建好模型后導入 模型建好后，以上兩種建模方法的具體過程將不再顯示,劃分網(wǎng)格,首先定義單元屬性: 單元類型, 實常數(shù), 材料屬性. 單元類型 下表給出了常用的熱單元類型 每個結點只有一個自由度: 溫度,劃分網(wǎng)格(續(xù)),材料屬性 必須輸入導熱系數(shù), KXX 如果施加了內(nèi)部熱生成率，則需指定比熱 (C) ANSYS提供的材料庫 (/ansys57/matlib)

10、包括幾種常用材料的結構屬性 和熱屬性, 但是建議用戶創(chuàng)建、使用自己的材料庫 把優(yōu)先設置為 “熱分析” ，使材料模型圖形用戶界面只顯示材料的熱屬性 實常數(shù) 主要應用于殼單元和線單元,劃分網(wǎng)格(續(xù)),劃分網(wǎng)格 存儲數(shù)據(jù)文件 使用 MeshTool 劃分網(wǎng)格，使用缺省的智能網(wǎng)格劃分級別6可以生成很好的初始網(wǎng)格 至此完成前處理，下面開始求解,施加載荷計算, 定義分析類型 如果進行新的熱分析： Command: ANTYPE, STATIC, NEW GUI: Main menuSolution-Analysis Type-New AnalysisSteady-state 如果繼續(xù)上一次分析，比如增加邊

11、界條件等： Command: ANTYPE, STATIC, REST GUI: Main menuSolutionAnalysis Type-Restart,施加載荷計算(續(xù)), 施加載荷 可以直接在實體模型或單元模型上施加五種載荷(邊界條件) ： a、恒定的溫度 通常作為自由度約束施加于溫度已知的邊界上。 Command Family:D GUI：Main MenuSolution-Loads-Apply-Thermal-Temperature,施加載荷計算(續(xù)),b、熱流率 熱流率作為節(jié)點集中載荷，主要用于線單元模型中(通常線單元模型不能施加對流或熱流密度載荷)，如果輸入的值為正，代表熱

12、流流入節(jié)點，即單元獲取熱量。如果溫度與熱流率同時施加在一節(jié)點上則ANSYS讀取溫度值進行計算。 注意：如果在實體單元的某一節(jié)點上施加熱流率，則此節(jié)點周圍的單元要密一些，在兩種導熱系數(shù)差別很大的兩個單元的公共節(jié)點上施加熱流率時，尤其要注意。此外，盡可能使用熱生成或熱流密度邊界條件，這樣結果會更精確些。 Command Family: F GUI：Main MenuSolution-Loads-Apply-Thermal-Heat Flow,施加載荷計算(續(xù)),c、對流 對流邊界條件作為面載施加于實體的外表面，計算與流體的熱交換，它僅可施加于實體和殼模型上，對于線模型，可以通過對流線單元LINK3

13、4考慮對流。 Command Family: SF GUI：Main MenuSolution-Loads-Apply-Thermal-Convection,施加載荷計算(續(xù)),d、熱流密度 熱流密度也是一種面載。當通過單位面積的熱流率已知或通過FLOTRAN CFD計算得到時，可以在模型相應的外表面施加熱流密度。如果輸入的值為正，代表熱流流入單元。熱流密度也僅適用于實體和殼單元。熱流密度與對流可以施加在同一外表面，但ANSYS僅讀取最后施加的面載進行計算。 Command Family: F GUI：Main MenuSolution-Loads-Apply-Thermal-Heat Flu

14、x,施加載荷計算(續(xù)),e、生熱率 生熱率作為體載施加于單元上，可以模擬化學反應生熱或電流生熱。它的單位是單位體積的熱流率。 Command Family: BF GUI：Main MenuSolution-Loads-Apply-Thermal-Heat Generat,施加載荷計算(續(xù)),定載荷步選項 對于一個熱分析，可以確定普通選項、非線性選項以及輸出控制。 a. 普通選項 時間選項：雖然對于穩(wěn)態(tài)熱分析，時間選項并沒有實際的物理意義，但它提供了一個方便的設置載荷步和載荷子步的方法。 Command: TIME GUI: Main MenuSolution -Load Step Opts-

15、 Time/Frequenc Time-Time Step/ Time and Substps,施加載荷計算(續(xù)),每載荷步中子步的數(shù)量或時間步大?。簩τ诜蔷€性分析，每一載荷步需要多個子步。 Command: NSUBST GUI: Main MenuSolution-Load Step Opts- Time/Frequenc Time and Substps Command: DELTIM GUI: Main MenuSolution-Load Step Opts- Time/Frequenc Time-Time Step,施加載荷計算(續(xù)),遞進或階越選項：如果定義階越(stepped)選

16、項，載荷值在這個載荷步內(nèi)保持不變；如果為遞進(ramped)選項，則載荷值由上一載荷步值到本載荷步值隨每一子步線性變化。 Command: KBC GUI: Main MenuSolution-Load Step Opts-Time/FrequencTime-Time Step/Time and Substps,施加載荷計算(續(xù)),b. 非線性選項 迭代次數(shù)：本選項設置每一子步允許的最多的迭代次數(shù)。默認值為25，對大數(shù)熱分析問題足夠。 Command: NEQIT GUI: Main MenuSolution-Load Step Opts-NolinearEquilibrium Iter,施加

17、載荷計算(續(xù)),自動時間步長: 對于非線性問題，可以自動設定子步間載荷的增長，保證求解的穩(wěn)定性和準確性。 Command: AUTOTS GUI: Main MenuSolution-Load Step Opts-Time/FrequencTime-Time Step/Time and Substps,施加載荷計算(續(xù)),收斂誤差：可根據(jù)溫度、熱流率等檢驗熱分析的收斂性。 Command: CNVTOL GUI: Main MenuSolution-Load Step Opts-NolinearConvergence Crit,施加載荷計算(續(xù)),求解結束選項：如果在規(guī)定的迭代次數(shù)內(nèi)，達不到收

18、斂，ANSYS可以停止求解或到下一載荷步繼續(xù)求解。 Command: NCNV GUI: Main MenuSolution-Load Step Opts-NolinearCriteria to Stop,施加載荷計算(續(xù)),線性搜索：設置本選項可使ANSYS用Newton-Raphson方法進行線性搜索。 Command: LNSRCH GUI: Main MenuSolution-Load Step Opts-NolinearLine Search,施加載荷計算(續(xù)),預測矯正：本選項可激活每一子步第一次迭代對自由度求解的預測矯正。 Command: PRED GUI: Main Menu

19、Solution-Load Step Opts-NolinearPredictor,施加載荷計算(續(xù)),c. 輸出控制 控制打印輸出：本選項可將任何結果數(shù)據(jù)輸出到*.out 文件中。 Command: OUTPR GUI: Main MenuSolution-Load Step Opts-Output CtrlsSolu Printout,施加載荷計算(續(xù)),控制結果文件：控制*.rth的內(nèi)容。 Command: OUTRES GUI: Main MenuSolution-Load Step Opts-Output CtrlsDB/Results File,施加載荷計算(續(xù)),確定分析選項 a

20、. Newton-Raphson選項(僅對非線性分析有用) Command: NROPT GUI: Main Menu Solution Analysis Options,施加載荷計算(續(xù)),b. 選擇求解器：可選擇如下求解器中一個進行求解： Frontal solver(默認） Jacobi Conjugate Gradient(JCG) solver JCG out-of-memory solver Incomplete Cholesky Conjugate Gradient(ICCG) solver Pre-Conditioned Conjugate Gradient Solver(PC

21、G) Iterative(automatic solver selection option),施加載荷計算(續(xù)),Command: EQSLV GUI: Main MenuSolutionAnalysis Options 注意：熱分析可選用Iterative選項進行快速求解，但如下情況除外： 熱分析包含SURF19或SURF22或超單元； 熱輻射分析； 相變分析 需要restart an analysis,施加載荷計算(續(xù)),c. 確定絕對零度：在進行熱輻射分析時，要將目前的溫度值換算為絕對溫度。如果使用的溫度單位是攝氏度，此值應設定為273；如果使用的是華氏度，則為460。 Command

22、: TOFFST GUI: Main Menu Solution Analysis Options,施加載荷計算(續(xù)),保存模型: 點擊ANSYS工具條SAVE_DB。 求解 Command: SOLVE GUI: Main MenuSolutionCurrent LS,后處理,ANSYS將熱分析的結果寫入*.rth文件中，它包含如下數(shù)據(jù)： 基本數(shù)據(jù)： 節(jié)點溫度 導出數(shù)據(jù)： 節(jié)點及單元的熱流密度 節(jié)點及單元的熱梯度 單元熱流率 節(jié)點的反作用熱流率 其它,后處理(續(xù)),對于穩(wěn)態(tài)熱分析，可以使用POST1進行后處理，進入POST1后，讀入載荷步和子步： Command: SET GUI: Main

23、 MenuGeneral Postproc-Read Results-By Load Step,后處理(續(xù)),可以通過如下三種方式查看結果： 彩色云圖顯示 Command: PLNSOL, PLESOL, PLETAB等 GUI: Main Menu General Postproc Plot ResultsNodal Solu, Element Solu, Elem Table,后處理(續(xù)),對3-D 實體模型繪制云圖時，選項isosurfaces (等值面)是非常有用的. 用 /CTYPE 命令或 Utility Menu PlotCtrls Style Contours Contour

24、Style.,后處理(續(xù)), 矢量圖顯示 Command: PLVECT GUI: Main Menu General Postproc Plot ResultsPre-defined or Userdefined,后處理(續(xù)), 列表顯示 Command: PRNSOL, PRESOL, PRRSOL等 GUI: Main Menu General Postproc List ResultsNodal Solu, Element Solu, Reaction Solu,后處理(續(xù)),檢查結果是否正確 溫度是否在預期的范圍內(nèi)? 在指定溫度和熱流邊界的基礎上，估計預期的范圍 網(wǎng)格大小是否滿足精度

25、? 和受力分析一樣，可以畫出非均勻分布的溫度梯度 (單元解) 并找出高梯度的單元. 這些區(qū)域可作為重新定義網(wǎng)格時的參考 若節(jié)點溫度梯度（平均的）和單元溫度梯度（非平均的）之間的差別很大，則可能是網(wǎng)格劃分太粗糙,練習-穩(wěn)態(tài)傳熱分析,某一潛水艇可以簡化為一圓筒，它由三層組成，最外面一層為不銹鋼，中間為玻纖隔熱層，最里面為鋁層，筒內(nèi)為空氣，筒外為海水，求內(nèi)外壁面溫度及溫度分布。 幾何參數(shù)：筒外徑30 feet 總壁厚2 inch 不銹鋼層壁厚0.75 inch 玻纖層壁厚1 inch 鋁層壁厚0.25 inch 筒長200 feet 導熱系數(shù)：不銹鋼8.27BTU/hr.ft.oF 玻纖0.028B

26、TU/hr.ft.oF 鋁117.4BTU/hr.ft.oF,邊界條件,邊界條件 ： 空氣溫度 70oF 海水溫度 44.5oF 空氣對流系數(shù) 2.5 BTU/hr.ft2.oF 海水對流系數(shù) 80 BTU/hr.ft2.oF 沿垂直于圓筒軸線作橫截面，得到一圓環(huán)，取其中1度進行分析，如圖示。,log文件,/filename, Steady1 /title, Steady-state thermal analysis of submarine /units, BFT Ro=15!外徑(ft) Rss=15-(0.75/12)!不銹鋼層內(nèi)徑ft) Rins=15-(1.75/12)!玻璃纖維層內(nèi)

27、徑(ft) Ral=15-(2/12)!鋁層內(nèi)徑 (ft) Tair=70!潛水艇內(nèi)空氣溫度 Tsea=44.5!海水溫度 Kss=8.27!不銹鋼的導熱系數(shù) (BTU/hr.ft.oF) Kins=0.028 !玻璃纖維的導熱系數(shù) (BTU/hr.ft.oF),log文件(續(xù)),Kal=117.4!鋁的導熱系數(shù)(BTU/hr.ft.oF) Hair=2.5!空氣的對流系數(shù)(BTU/hr.ft2.oF) Hsea=80!海水的對流系數(shù)(BTU/hr.ft2.oF) /prep7 et,1,plane55!定義二維熱單元 mp,kxx,1,Kss!設定不銹鋼的導熱系數(shù) mp,kxx,2,Kins

28、!設定玻璃纖維的導熱系數(shù) mp,kxx,3,Kal!設定鋁的導熱系數(shù) pcirc,Ro,Rss,-0.5,0.5!創(chuàng)建幾何模型 pcirc,Rss,Rins,-0.5,0.5 pcirc,Rins,Ral,-0.5,0.5 aglue,all,log文件(續(xù)),numcmp,area lesize,1,16!設定劃分網(wǎng)格密度 lesize,4,4 lesize,14,5 lesize,16,2 eshape,2!設定為映射網(wǎng)格劃分 mat,1 amesh,1 mat,2 amesh,2 mat,3 amesh,3,log文件(續(xù)),/SOLU SFL,11,CONV,HAIR,TAIR!施加空

29、氣對流邊界 SFL,1,CONV,HSEA,TSEA!施加海水對流邊界 SOLVE /POST1 PLNSOL！輸出溫度彩色云圖 finish,菜單操作,菜單操作: 1.Utility MenuFilechange jobname, 輸入Steady1； 2.Utility MenuFilechange title,輸入Steady-state thermal analysis of submarine； 3.在命令行輸入:/units, BFT； 4.Main Menu: Preprocessor；,菜單操作(續(xù)),5.Main Menu: Preprocessor Element Type

30、 Add/Edit/Delete ,選擇PLANE55； 6.Main Menu: Preprocessor Material Prop -Constant - Isotropic,默認材料編號為1,在KXX框中輸入8.27,選擇APPLY,輸入材料編號為2,在KXX框中輸入0.028,選擇APPLY,輸入材料編號為3,在KXX框中輸入117.4；,菜單操作(續(xù)),7.Main Menu: Preprocessor-Modeling-Create-Areas-CircleBy Dimensions 在RAD1中輸入15,在RAD2中輸入15-(.75/12),在THERA1中輸入-0.5,在T

31、HERA2中輸入0.5,選擇APPLY; 在RAD1中輸入15-(.75/12),在RAD2中輸入15-(1.75/12),選擇APPLY; 在RAD1中輸入15-(1.75/12),在RAD2中輸入15-2/12,選擇OK；,菜單操作(續(xù)),8.Main Menu: Preprocessor-Modeling-Operate-Booleane-GlueArea,選擇PICK ALL；,菜單操作(續(xù)),9.Main Menu: Preprocessor-Meshing-Size Contrls-Lines-Picked Lines 選擇不銹鋼層短邊,在NDIV框中輸入4,選擇APPLY;,菜單

32、操作(續(xù)),選擇玻璃纖維層的短邊,在NDIV框中輸入5,選擇APPLY; 選擇鋁層的短邊,在NDIV框中輸入2,選擇APPLY; 選擇四個長邊,在NDIV中輸入16；,菜單操作(續(xù)),10.Main Menu: Preprocessor-Attributes-Define Picked Area 選擇不銹鋼層,在MAT框中輸入1,選擇APPLY; 選擇玻璃纖維層,在MAT框中輸入2,選擇APPLY; 選擇鋁層,在MAT框中輸入3,選擇OK；,菜單操作(續(xù)),11.Main Menu: Preprocessor-Meshing-Mesh-Areas-Mapped3 or 4 sided,選擇PI

33、CK ALL；,菜單操作(續(xù)),12.Main Menu: Solution-Loads-Apply-Thermal-Convection On lines 選擇不銹鋼外壁,在VALI框中輸入80,在VAL2I框中輸入44.5,選擇APPLY;,菜單操作(續(xù)),選擇鋁層內(nèi)壁,在VALI框中輸入2.5,在VAL2I框中輸入70,選擇OK；,菜單操作(續(xù)),13.Main Menu: Solution-Solve-Current LS； 14.Main Menu: General PostprocPlot Results-Contour Plot-Nodal Solu,選擇Temperature。

34、,ANSYS Verification Manual中關于穩(wěn)態(tài)熱分析的實例,VM58Centerline temperature of a heat generating wire VM92Insulted wall temperature VM93Temperature dependent conductivity VM94Heat generating plate VM95Heat transfer from a cooling spine VM96Temperature distribution in a short solid cylinder VM97Temperature dist

35、ribution along a straight fin VM98Temperature distribution along a tapered fin VM99Temperature distribution in a trapezoidal fin VM100Heat conductivity across a chimney section,ANSYS Verification Manual中關于穩(wěn)態(tài)熱分析的實例(續(xù)),VM101Temperature distribution in a short solid cylinder VM102Cylinder with temperat

36、ure dependent conductivity VM103Thin plate with a central heat source VM105Heat generation coil with temperature dependent conductivity VM108Temperature gradient across a solid cylinder VM118Centerline temperature of a heat generating wire VM160Solid cylinder with harmonic temperature load VM161Heat

37、 flow from a insulated pipe VM162Cooling of a circular fin of rectangular profile VM193Adaptive analysis of two-dimensional heat transfer with convection,練習-瞬態(tài)傳熱分析,一個30公斤重、溫度為70的銅塊,以及一個20公斤重、溫度為80的鐵塊，突然放入溫度為20、盛滿了300升水的、完全絕熱的水箱中,如圖所示。過了一個小時，求銅塊與鐵塊的最高溫度(假設忽略水的流動)。 材料熱物理性能如下:,練習-瞬態(tài)傳熱分析(續(xù)),log文件,/filen

38、ame,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,log文件(續(xù)),rectnag,0,0.6,0,0.5!創(chuàng)建幾何實體 rectang,0.15,0.225,0.225,0.27 rectang,0.

39、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,log文件(續(xù)),/pnum,mat,1 eplot finish !加載求解 /solu antype,trans timint,off!先作穩(wěn)態(tài)分析,確定初始條件 time,0.01!設定只有一個子步的時間很小的載荷步 deltim,0.01 e

40、sel,s,mat,3 nsle,s d,all,temp,20 esel,s,mat,2 nsle,s d,all,temp,80,log文件(續(xù)),esel,s,mat,1 nsle,s d,all,temp,70 allsel solve !得到初始溫度分布 time,3600 !進行瞬態(tài)分析 timint,on !打開時間積分 deltim,26,2,200 !設置時間步長，最大及最小時間步長 autots,on !打開自動時間步長 ddelet,all,temp !刪除穩(wěn)態(tài)分析中定義的節(jié)點溫度 outres,all,1 !將每個子步的值寫入數(shù)據(jù)庫文件 solve finish sav

41、e,log文件(續(xù)),!進入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,

42、8,node(60,25,0),temp,T_H2O_Right !8水箱的右部 Plvar,2 plvar,3,4,5,6,7,8 finish,log文件(續(xù)),!進入POST1后處理 /post1!設置為最后一個載荷子步 set,last esel,s,mat,1 nsle,s plnsol,temp esel,s,mat,2 nsle,s plnsol,temp finish,菜單操作,1、Utility Menu: FileChange Jobname, 輸入文件名Transient1； 2、Utility Menu: FileChange Title,輸入Thermal Trans

43、ient Exercise 1；,19.16.2 菜單操作(續(xù)),3、Main MenuPreprocessor，進入前處理； 4、Main MenuPreprocessorElement Type Add/Edit/Delete, 選擇單元plane77；,菜單操作(續(xù)),5、Main MenuPreprocessorMaterial Props Isotropic, 定義材料1（銅）的KXX等于 383、DENS等于8889、C等于1390；,菜單操作(續(xù)),6、Main MenuPreprocessorMaterial Props Isotropic, 定義材料2（鐵）的KXX等于 70

44、、DENS等于7837、C等于448； 7、Main MenuPreprocessorMaterial Props Isotropic, 定義材料3（水）的KXX等于0.61、DENS等于996、C等于4185；,菜單操作(續(xù)),8、Main MenuPreprocessor-Modeling- Create-Areas- Retangle By Dimensions 輸入X1=0, Y1=0, X2=0.6, Y2=0.5, 選擇Apply;,菜單操作(續(xù)),輸入X1=0.15, Y1=0.225, X2= 0.225, Y2=0.27, 選擇Apply; 輸入X1=0.6-0.2-0.05

45、8, Y1=0.225, X2=0.6-0.2, Y2=0.225+0.044, 選擇OK.,菜單操作(續(xù)),9、Main MenuPreprocessor-Modeling- Operate Booleans Overlap , 選擇Pick All； 10、Utility Menu: PlotctrlsNumberingAreas, on； 11、Utility Menu:PlotAreas；,菜單操作(續(xù)),12、Main MenuPreprocessor-Attributes-Define-All Areas, 選擇材料1； 13、Main MenuPreprocessorMeshin

46、g-Size Cntrls-Manualsize-Global-Size, 輸入單元大小0.02； 14、Main MenuPreprocessorMeshing-Mesh-Areas- Mapped 3 or 4 sided , 選擇銅塊；,菜單操作(續(xù)),15、Main MenuPreprocessor-Attributes-Define-All Areas, 選擇材料2； 16、Main MenuPreprocessorMeshing-Mesh-Areas- Mapped 3 or 4 sided, 選擇鐵塊； 17、Main MenuPreprocessor-Attributes-De

47、fine-All Areas, 選擇材料3； 18、Main MenuPreprocessorMeshing-Size Cntrls-Manualsize-Global - Size, 輸入單元大小0.05； 19、Main MenuPreprocessorMeshing-Mesh-Areas-Free, 選擇水箱；,菜單操作(續(xù)),菜單操作(續(xù)),20、Utility MenuPlotArea； 21、Main MenuSolution，進入加載求解； 22、Main MenuSolution-Analysis Type-New Analysis, 選擇Transient，定義為瞬態(tài)分析；,

48、菜單操作(續(xù)),23、Main MenuSolution-Load Step Opts Time/FrenquencTime Integration, TIMINT, off，首先進行穩(wěn)態(tài)分析；,菜單操作(續(xù)),24、Main MenuSolution-Load Step Opts Time/FrenquencTime-Time Step，設定TIME為0.01、DELTIM也為0.01；,菜單操作(續(xù)),25、Utility Menu: SelectElementmat, 輸入3, 選擇Apply, 選擇NodesAttached to, 選擇Element；,菜單操作(續(xù)),26、Main MenuSolution-Loads-Apply-Thermal-Temperat