有限元熱分析基本概念課件

有限元熱分析基本概念符號下列符號在全文中的意義如下：fANSYS中標準單位(SI)溫度熱流量熱傳導率密度比熱對流換熱系數(shù)熱流溫度梯度內部熱生成DegreesC(orK)WattsWatts/(meter-degreeC)kilogram/(meter3)(Watt-sec)/(kilogram-degreeC)Watt/(meter2-degreeC)Watt/(meter2)degreeC/meterWatt/(meter3)傳導傳導的熱流由傳導的傅立葉定律決定：負號表示熱沿梯度的反向流動(i.e.,熱從熱的部分流向冷的).Tnq*dTdn對流對流的熱流由冷卻的牛頓準則得出:對流一般作為面邊界條件施加TBTs輻射從平面

i到平面

j的輻射熱流由施蒂芬-玻斯曼定律得出:在ANSYS中將輻射按平面現(xiàn)象處理(i.e.,體都假設為不透明的)。ij熱力學第一定律能量守恒要求系統(tǒng)的能量改變與系統(tǒng)邊界處傳遞的熱和功數(shù)值相等。能量守恒在一個短的時間增量下可以表示為方程形式將其應用到一個微元體上，就可以得到熱傳導的控制微分方程。EEEEstoredinthrutheboundaryoutthrutheboundarygenerated+++=0有限元方法將控制微分方程轉化為等小的積分形式(參閱ANSYS理論手冊第6.1節(jié))。有限元方法(續(xù))假設單元內溫度變化可以用多項式表示。一般情況下，根據(jù)單元類型的不同，應當包含不同的一次項,平方和混合的立方項。多項式假設保證了溫度在單元內部和單元邊界上都是連續(xù)的。

寫出以單元結點溫度為未知數(shù)的多項式:有限元方法(續(xù))將假設的溫度變化代入積分方程，注意到每項都乘上了實際的溫度數(shù)值，將兩邊約去得到有限元方法(續(xù))方程可以重新寫為簡化形式:有限元方法(續(xù))其中,有限元方法(續(xù))系統(tǒng)方程是將單元的貢獻組裝而成有限元方法(續(xù))

尺寸分析由前面的方程可以很快得出我們是否需要使用與幾何尺寸有關的單位:例子:3結點三角形單元(續(xù))有限元單元模型: 2三角形單元 4結點推導單元1矩陣:1243xy121231單元形函數(shù)例子:3結點三角形單元(續(xù))推出梯度-溫度矩陣定義各向同性材料特性矩陣單元傳導矩陣例子:3結點三角形單元(續(xù))同樣得到單元2的矩陣并組合成為總體矩陣矩陣可以分塊如圖，因為T3=T4=0同時求解得到未知的溫度求解單元1結點3的響應熱流Q3

和Q4

是響應熱流例子:3結點三角形單元(續(xù))計算單元1的溫度梯度向量*計算單元1的熱流向量*注意在單元內部梯度和熱流是均勻的*-向量只有一項，因為梯度/熱流假設在單元中是均勻的。有限元熱分析中的基本符號（續(xù)）一般來說，穩(wěn)態(tài)分析中網(wǎng)格上結點溫度比實際溫度要低。也就是說，如果加密網(wǎng)格，溫度將增加，但加密到一定程度，結果將不顯著增加(i.e.,結果收斂)。T網(wǎng)格密度有限元熱分析中的基本符號（續(xù)）引起奇異性的原因整體求解的奇異性在穩(wěn)態(tài)分析中當有熱量輸入(e.g.,施加結點熱流，熱流，內部熱源)而無熱流流出(指定的結點溫度，對流載荷等)，穩(wěn)態(tài)的溫度將是無限大的。等同于結構分析中的剛體位移。溫度梯度/熱流奇異性如果對點熱源處的網(wǎng)格細分下去的話，梯度/熱流將無限增加。凹角和網(wǎng)格中的“裂縫”。形狀不好的單元。網(wǎng)格劃分誤差實際上任何產(chǎn)生不連續(xù)熱流區(qū)域的有限元模型都是有誤差的。在單元內部邊界上熱流不連續(xù)的大小將作為ANSYS進行誤差估計的基礎。網(wǎng)格劃分誤差估計一般用于實體和殼單元，而且單元所在區(qū)域的單元類型是均一的(e.g.,具有共同的特性)，熱流在該區(qū)域中也就是連續(xù)的。誤差計算的細節(jié)在ANSYS理論手冊,19.7.2部分中有敘述。ANSYS誤差度量(續(xù))網(wǎng)格劃分誤差度量(續(xù))誤差限SMNB和SMXB-當用云圖繪制不連續(xù)數(shù)值(溫度梯度和熱流)時(誤差估計功能處于打開狀態(tài)),SMNB和SMXB將出現(xiàn)在圖例區(qū)域，表示出該數(shù)值不連續(xù)的范圍。ANSYS誤差估計(續(xù))網(wǎng)格劃分誤差度量(續(xù))例:假如云圖顯示的是X方向的平均結點熱流(PLNS,TF,X),SMNB和SMXB將顯示在圖例中，其計算方法如下:如何使熱傳遞分析包括非線性?當比熱矩陣，熱傳導率矩陣和/或等效結點熱流向量是溫度的函數(shù)時，分析就是非線性的，需要迭代求解平衡方程。如果所有三項都是與溫度有關的，那么控制方程可以寫為如下形式:下面幾項都可以使得分析包括非線性:與溫度有關的材料特性與溫度有關的對流換熱系數(shù)使用輻射單元與溫度有關的熱源(熱流或熱流矢量)使用耦合場單元(假設載荷向量耦合)何時需要定義比熱和密度?瞬態(tài)問題,這些數(shù)值用于形成比熱矩陣(該矩陣表示瞬態(tài)分析中需要的熱能存儲效果).穩(wěn)態(tài)分析中包括有熱質量傳遞效果(i.e.,模型中有流動導體介質).與結構分析的比較結構位移力均布載荷應變應力溫度分布

內部載荷塑性基礎無接觸 熱溫度熱流率熱流(施加的)溫度梯度熱流(計算的)內部熱生成(heat/volume)無對流輻射恒溫器對于熟悉結構分析的人來說，下面的表格將是非常有幫助的:單元庫1-Dthermalnetworkelements(可以使用在1-D,2-D和3-D單元中)。LINK31RadiationLinkLINK322-DConductionBarLINK333-DConductionBar單元庫1-Dthermalnetworkelements(續(xù))LINK34Node-NodeConvectionLinkMASS71LumpedThermalMass可以用于定義與溫度有關的熱源單元庫控制單元-允許用戶在有限元模型中加入反饋。最簡單的方法-恒溫器!根據(jù)控制結點K或L的溫度或溫度差，一階或二階導數(shù)，溫度積分，或時間，程序可以打開或關閉結點I和J之間的熱流。COMBIN37Node-NodeControlElement單元庫2-Dsolids-傳導平面或軸對稱(幾何,載荷,材料特性)對于軸對稱性質，全局笛卡兒坐標x為徑向且所有x-坐標必須0PLANE55PLANE35單元庫2-Dsolids-傳導(續(xù))PLANE77單元庫2-Dsolids-傳導(續(xù))諧波單元(幾何和材料特性軸對稱但邊界條件非軸對稱)。全局笛卡兒坐標x為徑向且所有x-坐標必須0。使用傅立葉級數(shù)載荷迭加技術,因此單元限于線性分析。PLANE75PLANE78單元庫3-Dsolids

-conductionSOLID70SOLID87SOLID90單元庫橋單元包括2-D平面內的傳導和平面外的對流SHELL57單元庫1-D熱-流單元同時求解帶泵效應的1-D伯努利方程和1-D帶質量傳遞效果的熱傳遞(耦合場)可以在對流中連接平面效果單元用戶可對邊界