第15講計(jì)算流體動(dòng)力學(xué)和熱應(yīng)力場(chǎng)分析

1、1第第10章章 計(jì)算流體動(dòng)力學(xué)分析計(jì)算流體動(dòng)力學(xué)分析第15講2一、流體流動(dòng)有限單元法一、流體流動(dòng)有限單元法1 不可壓縮流體流動(dòng)的有限單元法不可壓縮流體流動(dòng)的有限單元法1) 無(wú)升力物體繞流無(wú)升力物體繞流（1）流動(dòng)方程和邊界條件）流動(dòng)方程和邊界條件流函數(shù)流函數(shù) 和勢(shì)函數(shù)和勢(shì)函數(shù) 均滿足拉均滿足拉普拉斯（普拉斯（Laplace）方程）方程0022222222yxyx第一類邊界條件第一類邊界條件滿足迪里西來(lái)滿足迪里西來(lái)(Dirichlet)條件條件第二類邊界條件第二類邊界條件滿足諾伊曼滿足諾伊曼(Neumann)條件條件nsnn322221)(dsdxdyyxs拉普拉斯型流函數(shù)方程的泛函表達(dá)式拉普拉斯

2、型流函數(shù)方程的泛函表達(dá)式dsdxdyyyxxs)(2其極值的必要條件是泛函的一階變分其極值的必要條件是泛函的一階變分=0，由此得，由此得邊界上已知切線方向速度邊界上已知切線方向速度40222dsndxdyyxs迦遼金加權(quán)余量法方程迦遼金加權(quán)余量法方程dsdxdyyyxxs)(2應(yīng)用格林（應(yīng)用格林（Green）公式，上式變換為）公式，上式變換為5（2）單元及整體分析）單元及整體分析圓柱繞流為例圓柱繞流為例6選取三角形單元為子域，則其線性插值函數(shù)可取為選取三角形單元為子域，則其線性插值函數(shù)可取為321321 1aaayxyaxaaei根據(jù)泛函分析，在全區(qū)域內(nèi)滿足的迦遼金表達(dá)式。根據(jù)泛函分析，在全區(qū)

3、域內(nèi)滿足的迦遼金表達(dá)式?？梢詫懗稍谌舾蓚€(gè)單元上滿足的迦遼金表達(dá)式可以寫成在若干個(gè)單元上滿足的迦遼金表達(dá)式Eeeedxdyyyxxdxdyyyxxe1只要求解各個(gè)單元的迦遼金式，然后只要求解各個(gè)單元的迦遼金式，然后求和，即可求得全區(qū)域的流場(chǎng)問(wèn)題求和，即可求得全區(qū)域的流場(chǎng)問(wèn)題7三角形單元結(jié)點(diǎn)三角形單元結(jié)點(diǎn)i，j，m的坐標(biāo)值的坐標(biāo)值mjiiieN,dsNyNyNxNxNisiijij)(迦遼金表達(dá)式迦遼金表達(dá)式系數(shù)矩陣系數(shù)矩陣待求的流函數(shù)變量待求的流函數(shù)變量右端列向量右端列向量82有升力物體繞流有升力物體繞流假設(shè)流函數(shù)表達(dá)式為假設(shè)流函數(shù)表達(dá)式為邊界條件邊界條件v01s9邊界條件邊界條件10在在S處滿

4、足庫(kù)達(dá)處滿足庫(kù)達(dá)(Kutta)條件，則：條件，則： (vx)s=0, (vy)s=0機(jī)翼后緣點(diǎn)機(jī)翼后緣點(diǎn)S處的速度處的速度11解這二元一次代數(shù)方程組解這二元一次代數(shù)方程組已求出流函數(shù)的已求出流函數(shù)的兩個(gè)待定系數(shù)兩個(gè)待定系數(shù)全流場(chǎng)的流函全流場(chǎng)的流函數(shù)即為已知數(shù)即為已知12雅可比行列式雅可比行列式133 不可壓縮粘性流體流動(dòng)平面不可壓縮粘性流體流動(dòng)的連續(xù)性方程和運(yùn)動(dòng)方程平面不可壓縮粘性流體流動(dòng)的連續(xù)性方程和運(yùn)動(dòng)方程14一、一、FLOTRAN計(jì)算流體動(dòng)力學(xué)（計(jì)算流體動(dòng)力學(xué)（CFD）分析概述）分析概述二、二、FLOTRAN分析基礎(chǔ)分析基礎(chǔ)三、三、FLOTRAN不可壓層流和湍流分析不可壓層流和湍流分析四

5、、四、FLOTRAN熱分析熱分析五、五、FLOTRAN瞬態(tài)分析瞬態(tài)分析六、六、FLOTRAN可壓縮流分析可壓縮流分析15一、一、FLOTRAN計(jì)算流體動(dòng)力學(xué)（計(jì)算流體動(dòng)力學(xué)（CFD）分析概述）分析概述二、二、FLOTRAN分析基礎(chǔ)分析基礎(chǔ)三、三、FLOTRAN不可壓層流和湍流分析不可壓層流和湍流分析16一、一、FLOTRAN計(jì)算流體動(dòng)力學(xué)計(jì)算流體動(dòng)力學(xué)(CFD)分析概述分析概述(一一)、 FLOTRAN CFD分析二維及三維流體流動(dòng)場(chǎng)分析二維及三維流體流動(dòng)場(chǎng)可解決的問(wèn)題可解決的問(wèn)題作用于氣動(dòng)翼（葉）型上的升力和阻力；作用于氣動(dòng)翼（葉）型上的升力和阻力；超音速噴管中的流場(chǎng)；超音速噴管中的流場(chǎng)；彎

6、管中流體的復(fù)雜的三維流動(dòng)；彎管中流體的復(fù)雜的三維流動(dòng)；計(jì)算發(fā)動(dòng)機(jī)排氣系統(tǒng)中氣體的壓力及溫度分布；計(jì)算發(fā)動(dòng)機(jī)排氣系統(tǒng)中氣體的壓力及溫度分布；研究管路系統(tǒng)中熱的層化及分離；研究管路系統(tǒng)中熱的層化及分離；使用混合流研究來(lái)估計(jì)熱沖擊的可能性；使用混合流研究來(lái)估計(jì)熱沖擊的可能性；用自然對(duì)流分析來(lái)估計(jì)電子封裝芯片的熱性能；用自然對(duì)流分析來(lái)估計(jì)電子封裝芯片的熱性能；對(duì)含有多種流體的對(duì)含有多種流體的(由固體隔開(kāi)由固體隔開(kāi))熱交換器進(jìn)行研究熱交換器進(jìn)行研究17層流或湍流；層流或湍流；傳熱或絕熱；傳熱或絕熱；可壓縮或不可壓縮；可壓縮或不可壓縮；牛頓流或非牛頓流；牛頓流或非牛頓流；多組份輸運(yùn)。多組份輸運(yùn)。(二二)、

7、FLOTRAN分析的種類分析的種類一個(gè)層流分析可以是傳熱的或者是絕熱的，一個(gè)一個(gè)層流分析可以是傳熱的或者是絕熱的，一個(gè)湍流分析可以是可壓縮的或者是不可壓縮的。湍流分析可以是可壓縮的或者是不可壓縮的。18二、二、FLOTRAN分析基礎(chǔ)分析基礎(chǔ)(一一)、FLOTRAN單元的特點(diǎn)單元的特點(diǎn)FLUID141FLUID142用于解算單相粘性流體的二維用于解算單相粘性流體的二維和三維流動(dòng)、壓力和溫度分布和三維流動(dòng)、壓力和溫度分布單單元元名名稱稱維維數(shù)數(shù)形形狀狀自自由由度度F FL LU UI ID D1 14 41 1二二維維4 4節(jié)節(jié)點(diǎn)點(diǎn)四四邊邊形形或或3 3節(jié)節(jié)點(diǎn)點(diǎn)三三角角形形F FL LU UI I

8、D D1 14 42 2三三維維4 4節(jié)節(jié)點(diǎn)點(diǎn)四四面面體體或或8 8節(jié)節(jié)點(diǎn)點(diǎn)六六面面體體速速度度、壓壓力力、溫溫度度、湍湍流流動(dòng)動(dòng)能能、湍湍流流動(dòng)動(dòng)能能耗耗散散率率，多多達(dá)達(dá)六六種種流流體體的的各各自自質(zhì)質(zhì)量量所所占占的的份份額額FLOTRAN單元的特點(diǎn)單元的特點(diǎn)191. 用于模擬湍流的二方程湍流模式；用于模擬湍流的二方程湍流模式；2. 有很多推導(dǎo)結(jié)果，諸如：流場(chǎng)分析中的馬赫數(shù)、壓力系數(shù)、有很多推導(dǎo)結(jié)果，諸如：流場(chǎng)分析中的馬赫數(shù)、壓力系數(shù)、總壓、剪應(yīng)力、壁面處的總壓、剪應(yīng)力、壁面處的y-plus、以及流函數(shù)、熱分析中的、以及流函數(shù)、熱分析中的熱流、對(duì)流換熱（膜）系數(shù)等；熱流、對(duì)流換熱（膜）系數(shù)

9、等；3. 流體邊界條件，包括：速度、壓力、湍流動(dòng)能以及湍流動(dòng)流體邊界條件，包括：速度、壓力、湍流動(dòng)能以及湍流動(dòng)能耗散率。能耗散率。4. 熱邊界條件，包括：溫度、熱流、體積熱源、熱交換（膜）熱邊界條件，包括：溫度、熱流、體積熱源、熱交換（膜）系數(shù)。系數(shù)。FLUID單元的其它特征單元的其它特征20(二二)、FLOTRAN分析的主要步驟分析的主要步驟1確定問(wèn)題的區(qū)域確定問(wèn)題的區(qū)域2確定流體的狀態(tài)確定流體的狀態(tài)3. 生成有限元網(wǎng)格生成有限元網(wǎng)格4施加邊界條件施加邊界條件5設(shè)置設(shè)置FLOTRAN分析參數(shù)分析參數(shù)6求解求解7檢查結(jié)果檢查結(jié)果估計(jì)流體的特征，流體的特征估計(jì)流體的特征，流體的特征是流體性質(zhì)、幾

10、何邊界以及流是流體性質(zhì)、幾何邊界以及流場(chǎng)的速度幅值的函數(shù)場(chǎng)的速度幅值的函數(shù)應(yīng)使用映射網(wǎng)格劃分應(yīng)使用映射網(wǎng)格劃分湍流模型或求解溫度方程等湍流模型或求解溫度方程等選項(xiàng)，用戶必須激活它們選項(xiàng)，用戶必須激活它們這些變量包括速度、壓力、溫度、動(dòng)能（ENKE自由度）和動(dòng)能耗散率（ENDS自由度）等湍流量以及有效粘性（EVIS）。一個(gè)分析通常需要多次重啟動(dòng)。通過(guò)在觀察求解過(guò)程中相關(guān)變通過(guò)在觀察求解過(guò)程中相關(guān)變量的改變率，可以監(jiān)視求解的量的改變率，可以監(jiān)視求解的收斂性及穩(wěn)定性。收斂性及穩(wěn)定性。21(三三)、FLOTRAN分析中產(chǎn)生的一些文件分析中產(chǎn)生的一些文件結(jié)果文件，結(jié)果文件，Jobname.RFL，包含節(jié)

11、點(diǎn)結(jié)果；，包含節(jié)點(diǎn)結(jié)果；打印文件，打印文件，Jobname.PFL，各量的收斂記錄及進(jìn)出口狀態(tài)，各量的收斂記錄及進(jìn)出口狀態(tài)壁面文件，壁面文件，Jobname.RSW，壁面剪切應(yīng)力以及，壁面剪切應(yīng)力以及Y一一Plus信息信息殘差文件，殘差文件，Jobname.RDF，包含節(jié)點(diǎn)殘差，包含節(jié)點(diǎn)殘差調(diào)試文件，調(diào)試文件，Jobname.DBG，包含數(shù)學(xué)求解器的有關(guān)信息，包含數(shù)學(xué)求解器的有關(guān)信息結(jié)果備份文件，結(jié)果備份文件，Jobname.RFO，結(jié)果文件數(shù)據(jù)的一個(gè)拷貝，結(jié)果文件數(shù)據(jù)的一個(gè)拷貝重啟動(dòng)文件，重啟動(dòng)文件，Jobname.CFD，包含，包含F(xiàn)LOTRAN的數(shù)據(jù)結(jié)構(gòu)的數(shù)據(jù)結(jié)構(gòu)22(四四)、提高收斂性

12、和穩(wěn)定性的常用的工具、提高收斂性和穩(wěn)定性的常用的工具1松弛因子松弛因子2慣性松弛慣性松弛3人工粘性人工粘性4速度限值速度限值5積分階次積分階次慣性松弛系數(shù)是以所加項(xiàng)的分母的形式出現(xiàn)的，故其值越小，所起作用越大，其典型值介于1.0（作用中等）到1.010-7（作用很大）之間。松弛因子是一個(gè)其值介于松弛因子是一個(gè)其值介于0和和1之間的小數(shù)，之間的小數(shù)，它表示舊結(jié)果與附加在舊結(jié)果上以形成新它表示舊結(jié)果與附加在舊結(jié)果上以形成新結(jié)果的最近一次計(jì)算量之間的變化量結(jié)果的最近一次計(jì)算量之間的變化量對(duì)某個(gè)自由度的方程組的慣性松弛就是使對(duì)某個(gè)自由度的方程組的慣性松弛就是使其矩陣的主對(duì)角占優(yōu)以保持求解的穩(wěn)定性。其矩

13、陣的主對(duì)角占優(yōu)以保持求解的穩(wěn)定性。人工粘性用于在梯度較大的區(qū)域平抑速度解人工粘性用于在梯度較大的區(qū)域平抑速度解使所求解量不能超出用戶所定義的值。它可使所求解量不能超出用戶所定義的值。它可對(duì)速度、壓力和溫度自由度進(jìn)行限制對(duì)速度、壓力和溫度自由度進(jìn)行限制(VX，VY，V2，PRES，TEMP）23(五五)、 FLOTRAN分析過(guò)程中應(yīng)考查的項(xiàng)目分析過(guò)程中應(yīng)考查的項(xiàng)目1確定總體迭代的數(shù)目確定總體迭代的數(shù)目2.收斂監(jiān)測(cè)收斂監(jiān)測(cè)3中斷中斷FLOTRAN求解求解FLOTRAN分析是一個(gè)非線性的分析是一個(gè)非線性的序列求解過(guò)程，故每次分析首先序列求解過(guò)程，故每次分析首先要確定讓程序執(zhí)行多少次迭代要確定讓程序執(zhí)

14、行多少次迭代程序在每一個(gè)總體迭代里對(duì)速度、壓力、溫度、程序在每一個(gè)總體迭代里對(duì)速度、壓力、溫度、湍流動(dòng)能、動(dòng)能耗散率以及激活的多組份傳輸湍流動(dòng)能、動(dòng)能耗散率以及激活的多組份傳輸方程計(jì)算出一個(gè)收斂監(jiān)測(cè)量。方程計(jì)算出一個(gè)收斂監(jiān)測(cè)量。ANSYS命令有收斂的缺省值，程序?qū)⒉粰z命令有收斂的缺省值，程序?qū)⒉粰z查未激活的自由度或設(shè)為負(fù)值的自由度。查未激活的自由度或設(shè)為負(fù)值的自由度。24(六六)、 對(duì)對(duì)FLOTRAN分析進(jìn)行評(píng)價(jià)分析進(jìn)行評(píng)價(jià)初始參數(shù)初始參數(shù)邊界條件邊界條件正確正確所有變量的收斂監(jiān)測(cè)量都停止增長(zhǎng)所有變量的收斂監(jiān)測(cè)量都停止增長(zhǎng)所有求解量的平均、最大、最小值都不再升降所有求解量的平均、最大、最小值都

15、不再升降求解過(guò)程就算是完成求解過(guò)程就算是完成不能保證所求解的結(jié)果是惟一正確的，不能保證所求解的結(jié)果是惟一正確的，因?yàn)樽匀唤绫旧聿⒉槐ＷC存在惟一解。因?yàn)樽匀唤绫旧聿⒉槐ＷC存在惟一解。振蕩問(wèn)題振蕩問(wèn)題(例如：柱體繞流的旋渦脫落問(wèn)題例如：柱體繞流的旋渦脫落問(wèn)題)用穩(wěn)態(tài)或瞬態(tài)求用穩(wěn)態(tài)或瞬態(tài)求解技術(shù)都不能得到一個(gè)穩(wěn)定的解。要驗(yàn)證一個(gè)分析是穩(wěn)定的解技術(shù)都不能得到一個(gè)穩(wěn)定的解。要驗(yàn)證一個(gè)分析是穩(wěn)定的或是振蕩的，可以通過(guò)對(duì)它執(zhí)行大量的迭代求解來(lái)實(shí)現(xiàn)。或是振蕩的，可以通過(guò)對(duì)它執(zhí)行大量的迭代求解來(lái)實(shí)現(xiàn)。25(七七)、驗(yàn)證結(jié)果、驗(yàn)證結(jié)果如果一個(gè)FLOTRAN分析得到了非預(yù)期的結(jié)果，則應(yīng)進(jìn)行下列所示的操作，這些操作的

16、大部分都可以在開(kāi)始分析前完成。即使只進(jìn)行了零次迭代，ANSYS也會(huì)生成一個(gè)Jobname.RFL文件并檢查所有的輸入數(shù)據(jù)。1. 檢查作為結(jié)果總結(jié)的一部分而打印出來(lái)的質(zhì)量檢查作為結(jié)果總結(jié)的一部分而打印出來(lái)的質(zhì)量平衡情況，內(nèi)部檢查將確定是否有通過(guò)模型的質(zhì)平衡情況，內(nèi)部檢查將確定是否有通過(guò)模型的質(zhì)量流，允許質(zhì)量流的邊界條件是：量流，允許質(zhì)量流的邊界條件是：確定的速度邊界條件；確定的速度邊界條件；確定的壓力邊界條件；確定的壓力邊界條件；未確定的邊界。未確定的邊界。262. 在在ANSYS里檢查邊界條件，以保證其正確性。里檢查邊界條件，以保證其正確性。3. 檢查所定義的流體性質(zhì)及其隨溫度的變化正確與否，

17、檢查所定義的流體性質(zhì)及其隨溫度的變化正確與否， 這可在這可在.RFL文件中方便地檢查。文件中方便地檢查。4. 檢查用以建立模型的單位制與用以定義流體性質(zhì)的單檢查用以建立模型的單位制與用以定義流體性質(zhì)的單 位制是否一致。位制是否一致。5. 確認(rèn)與所選選項(xiàng)相聯(lián)系的方程的求解是否正確。確認(rèn)與所選選項(xiàng)相聯(lián)系的方程的求解是否正確。6. 如果求解發(fā)散，可能的原因還有：有限元網(wǎng)格不夠精如果求解發(fā)散，可能的原因還有：有限元網(wǎng)格不夠精 細(xì)、或者鄰近出口處流場(chǎng)梯度大大。細(xì)、或者鄰近出口處流場(chǎng)梯度大大。7. 如果僅僅只有某個(gè)特定的量產(chǎn)生發(fā)散，則可將該量重如果僅僅只有某個(gè)特定的量產(chǎn)生發(fā)散，則可將該量重 新初始化到一個(gè)

18、單值，并作重啟動(dòng)分析。新初始化到一個(gè)單值，并作重啟動(dòng)分析。27三、三、FLOTRAN不可壓層流和湍流分析不可壓層流和湍流分析除了進(jìn)口和出口邊界外，流除了進(jìn)口和出口邊界外，流動(dòng)被壁面或?qū)ΨQ面所約束動(dòng)被壁面或?qū)ΨQ面所約束外流邊界通常是遠(yuǎn)場(chǎng)外流邊界通常是遠(yuǎn)場(chǎng)速度或壓力邊界條件速度或壓力邊界條件(一一)、流體流動(dòng)分析的特點(diǎn)、流體流動(dòng)分析的特點(diǎn)流動(dòng)流動(dòng)內(nèi)流內(nèi)流外流外流層流或湍流分析能計(jì)算層流或湍流分析能計(jì)算出在二維或三維幾何中出在二維或三維幾何中的流動(dòng)及壓力分布。的流動(dòng)及壓力分布。需指定密度和粘性需指定密度和粘性28不可壓縮流不可壓縮流可壓縮流可壓縮流區(qū)別區(qū)別在于狀態(tài)方程在于狀態(tài)方程和求解方法和求解方法

19、密度的變化只有通過(guò)重力加速度才密度的變化只有通過(guò)重力加速度才能驅(qū)動(dòng)流動(dòng)，溫度變化導(dǎo)致密度變能驅(qū)動(dòng)流動(dòng)，溫度變化導(dǎo)致密度變化?；?。由壓力變化引起的密度變化由壓力變化引起的密度變化明顯地影響動(dòng)量和能量方程明顯地影響動(dòng)量和能量方程29(二二)、激活湍流模型、激活湍流模型層流與層流與湍流的區(qū)別湍流的區(qū)別慣性輸運(yùn)與慣性輸運(yùn)與粘性輸運(yùn)之比粘性輸運(yùn)之比比值的比值的增加增加不穩(wěn)定性增大且開(kāi)不穩(wěn)定性增大且開(kāi)始出現(xiàn)速度脈動(dòng)始出現(xiàn)速度脈動(dòng)湍流模型在控制方程中，使用增大的粘性湍流模型在控制方程中，使用增大的粘性(有效粘性有效粘性)來(lái)考慮這些脈動(dòng)對(duì)平均流動(dòng)的影響來(lái)考慮這些脈動(dòng)對(duì)平均流動(dòng)的影響有效粘性層流粘性有效粘性層流

20、粘性 + 湍流粘性湍流粘性1雷諾數(shù)雷諾數(shù)無(wú)量綱雷諾數(shù)用于測(cè)無(wú)量綱雷諾數(shù)用于測(cè)量慣性力與粘性之比量慣性力與粘性之比Re= VLc 當(dāng)雷諾數(shù)超過(guò)當(dāng)雷諾數(shù)超過(guò)2300時(shí)，時(shí)，通常應(yīng)激活湍流模型通常應(yīng)激活湍流模型302確定分析是不是湍流確定分析是不是湍流3湍流比和進(jìn)口參數(shù)湍流比和進(jìn)口參數(shù)當(dāng)分析開(kāi)始時(shí)，有效粘性初始化為層流粘性的若干倍。不管湍流模型是否被激活，初始化都會(huì)發(fā)生。隱含倍數(shù)（也叫湍流比）對(duì)于大多數(shù)湍流分析的合理值是1000。4. 湍流模型湍流模型FLOTRAN提供六種湍流模型提供六種湍流模型（1）標(biāo)準(zhǔn)）標(biāo)準(zhǔn)k一一 模型模型（2）零方程湍流模型）零方程湍流模型（3）Re一一Normalized

21、Group湍流模型湍流模型（4）Shih新新k一一 模型（模型（NKE） (5）Girimajd非線性模型非線性模型（6）ShihZhu，Lumley模型模型315標(biāo)準(zhǔn)標(biāo)準(zhǔn)k一一 模型模型6零方程湍流模型（零方程湍流模型（ZeroEq）7Re一一Normalized Group湍流模型（湍流模型（RNG）8. Shih新新k一一 模型（模型（NKE）9Girimaij非線性模型（非線性模型（GIR）10Shih，Zhu，Lumley模型（模型（SZL）32(三三)、 網(wǎng)格要求網(wǎng)格要求對(duì)于湍流網(wǎng)格要求比層流嚴(yán)格。顯然，最重要的區(qū)域就是對(duì)于湍流網(wǎng)格要求比層流嚴(yán)格。顯然，最重要的區(qū)域就是有較大梯度的

22、地方，尤其在壁面附近。與自由網(wǎng)格相反，有較大梯度的地方，尤其在壁面附近。與自由網(wǎng)格相反，在壁面結(jié)構(gòu)網(wǎng)格能提供更相容的模擬。在壁面結(jié)構(gòu)網(wǎng)格能提供更相容的模擬。非結(jié)構(gòu)網(wǎng)格非結(jié)構(gòu)網(wǎng)格結(jié)構(gòu)網(wǎng)格結(jié)構(gòu)網(wǎng)格33(四四)、 流動(dòng)邊界條件流動(dòng)邊界條件1.指定流量指定流量2.指定壓力指定壓力3.靜止壁面靜止壁面4.運(yùn)動(dòng)壁面運(yùn)動(dòng)壁面5.未指定邊界未指定邊界6.周期性邊界周期性邊界在邊界上指定所有速度分量在邊界上指定所有速度分量,需知進(jìn)口處密度需知進(jìn)口處密度一般在出口邊界施加相對(duì)壓力（通常為零）一般在出口邊界施加相對(duì)壓力（通常為零）無(wú)滑移條件，所有速度分量均設(shè)置為零無(wú)滑移條件，所有速度分量均設(shè)置為零指定壁面相切的速度

23、分量，指定壁面相切的速度分量，而所有其它的速度分量為零而所有其它的速度分量為零既不知道相對(duì)壓力，又不知道速度既不知道相對(duì)壓力，又不知道速度在兩個(gè)邊界上，邊界條件未知但又相等在兩個(gè)邊界上，邊界條件未知但又相等34四、實(shí)例應(yīng)用分析四、實(shí)例應(yīng)用分析p260-26735尺寸和特性尺寸和特性 入口長(zhǎng)度入口長(zhǎng)度4 in4 in入口高度入口高度1 in1 in過(guò)渡區(qū)長(zhǎng)過(guò)渡區(qū)長(zhǎng)度度2 in2 in出口高度出口高度2.5 in2.5 in初始出口初始出口長(zhǎng)度長(zhǎng)度4 in4 in擴(kuò)展后出擴(kuò)展后出口長(zhǎng)度口長(zhǎng)度30 in30 in空氣密度空氣密度1.21x101.21x10-7 -7 lblbf f-s-s2 2/i

24、n/in4 4空氣粘度空氣粘度2.642x102.642x10-9 -9 lblbf f-s/in-s/in2 2入口速度入口速度1 in/sec1 in/sec* *出口壓力出口壓力0 psi0 psi* *計(jì)算計(jì)算5050步步后重起動(dòng)后重起動(dòng) 4 42 21 12.52.54 4層流模擬：層流模擬：1.1.2D2D管道，模擬雷諾數(shù)為管道，模擬雷諾數(shù)為9090的層流問(wèn)題的層流問(wèn)題, ,速度為速度為1 in/s1 in/s 。2.2.獲得解后，增加速度到獲得解后，增加速度到50 in/s50 in/s，研究其對(duì)流場(chǎng)的影響獲得，研究其對(duì)流場(chǎng)的影響獲得新的解。新的解。3.3.增加管道長(zhǎng)度增加管道長(zhǎng)

25、度30in30in，研究流場(chǎng)變化。，研究流場(chǎng)變化。紊流模擬：紊流模擬：4.4.計(jì)算雷諾數(shù)大于計(jì)算雷諾數(shù)大于40004000時(shí)，用紊流模型重起動(dòng)求解。時(shí)，用紊流模型重起動(dòng)求解。 管道層流和紊流管道層流和紊流CFDCFD分析實(shí)例分析實(shí)例見(jiàn)見(jiàn)ANSYSHelp中中Step By Step36APDLAPDL程序：程序：2_Fluid.txt2_Fluid.txtfini/clear!* /filname, Fluid_CFD_Exe/PREP7 !* ET,1,FLUID141 !* RECTNG,0,4,0,1, RECTNG,6,10,0,2.5, ! LPLOT !* L2TAN,-3,-7

26、!A,3,2,5,8 ! /REPLOT !* LESIZE,1, , ,15,-2, , , ,1LESIZE,3, , ,15,-2, , , ,1!* LESIZE,9, , ,12,1, , , ,1 LESIZE,10, , ,12,1, , , ,1 !* !* LESIZE,5, , ,15,3, , , ,1 LESIZE,7, , ,15,1/3, , , ,1 !* LESIZE,4, , ,10,-2, , , ,1 LESIZE,2, , ,10,-2, , , ,1 LESIZE,8, , ,10,-2, , , ,1 LESIZE,6, , ,10,-2, , ,

27、,1 !* MSHAPE,0,2D MSHKEY,1amesh,all!DL,4, ,VX,1.0,1 DL,4, ,VY,0.0,1 !* !* DL,1, ,VX,0,1 DL,1, ,VY,0,1 DL,3, ,VX,0,1 DL,3, ,VY,0,1 DL,9, ,VX,0,1 DL,9, ,VY,0,1 DL,10, ,VX,0,1 DL,10, ,VY,0,1 DL,7, ,VX,0,1 DL,7, ,VY,0,1 DL,5, ,VX,0,1 DL,5, ,VY,0,1 !* DL,6, ,PRES,0,1 !* FINISH /SOLFLDATA12,PROP,DENS,8FLD

28、ATA12,PROP,VISC,8!* FLDATA7,PROT,DENS,AIR-INFLDATA8,NOMI,DENS,-1FLDATA7,PROT,VISC,AIR-INFLDATA8,NOMI,VISC,-1FLDATA7,PROT,COND,CONSTANT FLDATA8,NOMI,COND,-1, FLDATA7,PROT,SPHT,CONSTANT FLDATA8,NOMI,SPHT,-1, ! /COM,Steady State Analysis,0 FLDATA2,ITER,EXEC,40, FLDATA3,TERM,VX,0.01, FLDATA3,TERM,VY,0.0

29、1, FLDATA3,TERM,VZ,0.01, FLDATA3,TERM,PRES,1e-008, FLDATA3,TERM,TEMP,1e-008, FLDATA3,TERM,ENKE,0.01, FLDATA3,TERM,ENDS,0.01, FLDATA5,OUTP,SUMF,10, !* FLDATA15,PRES,REFE,14.7,FLDATA16,BULK,BETA,1e+015, FLDATA17,GAMM,COMP,1.4, FLDATA14,TEMP,NOMI,70, FLDATA14,TEMP,TTOT,70, FLDATA14,TEMP,BULK,70, TOFFSE

30、T,460,!* solve提示提示37第第9章章 溫度場(chǎng)和熱應(yīng)力場(chǎng)分析溫度場(chǎng)和熱應(yīng)力場(chǎng)分析 38一、一、 穩(wěn)態(tài)熱傳導(dǎo)問(wèn)題的有限元法穩(wěn)態(tài)熱傳導(dǎo)問(wèn)題的有限元法 1 熱傳導(dǎo)方程與換熱邊界熱傳導(dǎo)方程與換熱邊界 QzTzyTyxTxtTczyx在分析工程問(wèn)題時(shí)，經(jīng)常要了解工件內(nèi)部的溫度分布情況，例如發(fā)動(dòng)機(jī)的工作溫度、金屬工件在熱處理過(guò)程中的溫度變化、流體溫度分布等。物體內(nèi)部的溫度分布取決于物體內(nèi)部的熱量交換，以及物體與外部介質(zhì)之間的熱量交換，一般認(rèn)為是與時(shí)間相關(guān)的。 QzTyTxTtTc222222對(duì)于各向同性材料對(duì)于各向同性材料 對(duì)于各向異性材料對(duì)于各向異性材料 zy,x,TT00t初始條件和邊界條

31、件初始條件和邊界條件 391). 1). 給定物體邊界上的溫度，稱為第一類邊界條件給定物體邊界上的溫度，稱為第一類邊界條件 ssTT),(tzyxTTss或 2). 2). 給定物體邊界上的熱量輸入或輸出，稱為第二類邊界條件給定物體邊界上的熱量輸入或輸出，稱為第二類邊界條件 已知物體表面上熱流密度已知物體表面上熱流密度 ),()(tzyxqnzTnyTnxTsszzyyxx403). 3). 給定對(duì)流換熱條件，稱為第三類邊界條件給定對(duì)流換熱條件，稱為第三類邊界條件 物體與其相接觸的流體介質(zhì)之間的對(duì)流換熱系數(shù)和介質(zhì)的溫度為已知物體與其相接觸的流體介質(zhì)之間的對(duì)流換熱系數(shù)和介質(zhì)的溫度為已知 )(sf

32、zzyyxxTThnzTnyTnxT0zTyTxT222222Q各向同性的材料各向同性的材料 0zTyTxT222222物體不包含內(nèi)熱源物體不包含內(nèi)熱源 溫度場(chǎng)滿足溫度場(chǎng)滿足LaplaceLaplace方程方程 412 穩(wěn)態(tài)溫度場(chǎng)分析的一般有限元列式穩(wěn)態(tài)溫度場(chǎng)分析的一般有限元列式 穩(wěn)態(tài)溫度場(chǎng)計(jì)算是一個(gè)典型的穩(wěn)態(tài)溫度場(chǎng)計(jì)算是一個(gè)典型的場(chǎng)問(wèn)題場(chǎng)問(wèn)題 以二維問(wèn)題為例，說(shuō)明用以二維問(wèn)題為例，說(shuō)明用GalerkinGalerkin法建立穩(wěn)態(tài)溫法建立穩(wěn)態(tài)溫度場(chǎng)的一般有限元格式的過(guò)程度場(chǎng)的一般有限元格式的過(guò)程 0QyTyxTxyx穩(wěn)態(tài)熱傳導(dǎo)方程穩(wěn)態(tài)熱傳導(dǎo)方程 42第一類換熱邊界第一類換熱邊界 ssTT第二類

33、換熱邊界條件第二類換熱邊界條件 syyxxqnyTnxT第三類邊界條件第三類邊界條件 )(sfyyxxTThnyTnxT430)()(1dQyTyxTxwyxe在一個(gè)單元內(nèi)的加權(quán)積分公式在一個(gè)單元內(nèi)的加權(quán)積分公式 )()()(111xTxwxTxwxTwxxxx)()()(111yTywyTywyTwyyyy分分部部積積分分 44應(yīng)用應(yīng)用GreenGreen定理，一個(gè)單元內(nèi)的加權(quán)積分公式寫為定理，一個(gè)單元內(nèi)的加權(quán)積分公式寫為 0)()()(1111dnyTnxTwdQwyTywxTxwyyxxeyxe采用采用GalerkinGalerkin方法，選擇權(quán)函數(shù)為方法，選擇權(quán)函數(shù)為 iNw 1單元單

34、元的加的加權(quán)積權(quán)積分公分公式為式為 0)()(332dhTNdTNhNdqNdQNdTyNyNxNxNfieeieesiieeyixie第二類換熱邊界項(xiàng)第二類換熱邊界項(xiàng) 第三類換熱邊界項(xiàng)第三類換熱邊界項(xiàng) 450)()()()(332dhTNdTNNhdqNdQNdTyNyNxNxNfTeeTeesTTeeyTxTe矩陣形式矩陣形式 有限元格式有限元格式 eeePTKn n個(gè)個(gè)聯(lián)立聯(lián)立的線的線性方性方程組程組 KKe e為單元的導(dǎo)熱矩陣或稱為溫度剛度矩陣，為單元的導(dǎo)熱矩陣或稱為溫度剛度矩陣，TTe e為單元的結(jié)為單元的結(jié)點(diǎn)溫度向量，點(diǎn)溫度向量，PPe e稱為單元的溫度載荷向量或熱載荷向量稱為單元

35、的溫度載荷向量或熱載荷向量 460)()()()(332dhTNdTNNhdqNdQNdTyNyNxNxNfTeeeTeeesTeTeeeyTxTee整個(gè)物體上的加權(quán)積分方程是單元積分方程的和整個(gè)物體上的加權(quán)積分方程是單元積分方程的和 PTK整體方程組為整體方程組為 根據(jù)單元結(jié)點(diǎn)的局部編號(hào)與整體編號(hào)的關(guān)系，直接求和得到整體剛度矩陣根據(jù)單元結(jié)點(diǎn)的局部編號(hào)與整體編號(hào)的關(guān)系，直接求和得到整體剛度矩陣 473 三角形單元的有限元列式三角形單元的有限元列式 與計(jì)算彈性力學(xué)平面問(wèn)題時(shí)所采用的方法一樣，二維溫度場(chǎng)問(wèn)題計(jì)與計(jì)算彈性力學(xué)平面問(wèn)題時(shí)所采用的方法一樣，二維溫度場(chǎng)問(wèn)題計(jì)算中所采用的三角形單元算中所采用

36、的三角形單元( (如圖如圖9-19-1所示所示) )可以使用相同的形函數(shù)可以使用相同的形函數(shù) 圖圖9-1 9-1 三角形單元三角形單元 )(21ycxbaANiiii)(21ycxbaANjjjj)(21ycxbaANmmmmjmimjijmmjixxcyybyxyxamijimjmiimjxxcyybyxyxaijmjimijjimxxcyybyxyxa48 T111mmjjiiyxyxyxA2T mjimjiTTTNNNT在三角形單元上，采用在三角形單元上，采用GalerkinGalerkin法可得法可得 0)()(dAQyTyxTxNyxTA49假定單元內(nèi)的導(dǎo)熱系數(shù)為常數(shù)假定單元內(nèi)的導(dǎo)熱

37、系數(shù)為常數(shù) mjimmjmimjjjimijiixmjimjimjiAxTxATTTbbbbbbbbbbbbbbbAdATTTbbbbbbAdAxTxN2222441mjimmjmimjjjimijiixmjimjimjiAxTyATTTcccccccccccccccAdATTTccccccAdAyTyN222244150單元的剛度矩陣為單元的剛度矩陣為 22222244mmjmimjjjimijiiymmmimjjjijijiixecccccccccccccccbbbbbbbbbbbbbbbKj如果單元的內(nèi)部熱源為常數(shù)，由內(nèi)部熱源產(chǎn)生的溫度載荷項(xiàng)為如果單元的內(nèi)部熱源為常數(shù)，由內(nèi)部熱源產(chǎn)生的溫

38、度載荷項(xiàng)為 1113AQdANNNQdAQNmjiATA51由由GreenGreen公式可得公式可得 dSTNNhdSTNhdSTTNhdAyTNyxTNxeTsfTssfTsyTxTA)()(如果在單元邊如果在單元邊上存在熱交換，上存在熱交換，各條邊上的邊各條邊上的邊界換熱條件在界換熱條件在單元?jiǎng)偠染仃噯卧獎(jiǎng)偠染仃囍猩傻母郊又猩傻母郊禹?xiàng)為式項(xiàng)為式 0000210126ijehlK2101200006jmehlK2010001026miehlK圖圖9-1 9-1 三角形單元三角形單元 52由邊界換熱條件生成的溫由邊界換熱條件生成的溫度載荷向量為式度載荷向量為式 0112ijfelhTP11

39、02jmfelhTP1012mifelhTP圖圖9-1 9-1 三角形單元三角形單元 53二、熱彈性應(yīng)力問(wèn)題的有限元分析二、熱彈性應(yīng)力問(wèn)題的有限元分析 設(shè)在溫度設(shè)在溫度T0T0時(shí)物體處于無(wú)應(yīng)力狀態(tài)，當(dāng)物體內(nèi)發(fā)生溫度變時(shí)物體處于無(wú)應(yīng)力狀態(tài)，當(dāng)物體內(nèi)發(fā)生溫度變化化 T TT T1 1T T0 0時(shí)，物體中的微元體就要產(chǎn)生熱膨脹，對(duì)各時(shí)，物體中的微元體就要產(chǎn)生熱膨脹，對(duì)各向同性體，自由膨脹情況下的應(yīng)變分量為向同性體，自由膨脹情況下的應(yīng)變分量為 0zxyzxyTzyxC/熱膨脹系數(shù)熱膨脹系數(shù) 54如果自由膨脹受到某種約束，微元體就要產(chǎn)生熱應(yīng)力如果自由膨脹受到某種約束，微元體就要產(chǎn)生熱應(yīng)力 TEzyxx

40、/)(TEzxyy/TExzzz/ 0111eeDTID用矩陣的形式表示為用矩陣的形式表示為 彈性本構(gòu)矩陣彈性本構(gòu)矩陣 變溫而產(chǎn)生的應(yīng)變變溫而產(chǎn)生的應(yīng)變 55 0eD考慮變溫影響的彈性應(yīng)力，通常稱為熱應(yīng)力考慮變溫影響的彈性應(yīng)力，通常稱為熱應(yīng)力 TFFuK有限元的平衡方程有限元的平衡方程 56項(xiàng)目項(xiàng)目國(guó)際單位國(guó)際單位英制單位英制單位代號(hào)代號(hào)長(zhǎng)度長(zhǎng)度mft時(shí)間時(shí)間ss質(zhì)量質(zhì)量Kglbm溫度溫度oF力力Nlbf能量（熱量）能量（熱量）JBTU功率（熱流率）功率（熱流率）WBTU/sec熱流密度熱流密度W/m2BTU/sec-ft2生熱速率生熱速率W/m3BTU/sec-ft3導(dǎo)熱系數(shù)導(dǎo)熱系數(shù)W/m-

41、BTU/sec-ft-oFKXX對(duì)流系數(shù)對(duì)流系數(shù)W/m2-BTU/sec-ft2-oFHF密度密度Kg/m3lbm/ft3DENS比熱比熱J/Kg-BTU/lbm-oFC焓焓J/m3BTU/ft3ENTH三、三、ANSYS熱分析的應(yīng)用熱分析的應(yīng)用571 簡(jiǎn)介簡(jiǎn)介( (一一) )、熱分析的目的、熱分析的目的( (二二) )、 ANSYSANSYS熱分析的特點(diǎn)熱分析的特點(diǎn)( (三三) )、 ANSYSANSYS熱分析分類熱分析分類( (四四) )、 邊界條件、初始條件邊界條件、初始條件( (五五) )、 ANSYSANSYS熱分析誤差估計(jì)熱分析誤差估計(jì)58 熱分析用于計(jì)算一熱分析用于計(jì)算一個(gè)系統(tǒng)或

42、部件的溫度分個(gè)系統(tǒng)或部件的溫度分布及其它熱物理參數(shù)布及其它熱物理參數(shù)(一一)、熱分析的目的、熱分析的目的(1) 基于能量守恒原理的熱平衡方程，用有限元法基于能量守恒原理的熱平衡方程，用有限元法 計(jì)算各節(jié)點(diǎn)的溫度，并導(dǎo)出其它熱物理參數(shù)。計(jì)算各節(jié)點(diǎn)的溫度，并導(dǎo)出其它熱物理參數(shù)。(2) 包括熱傳導(dǎo)、熱對(duì)流及熱輻射三種熱傳遞方式。包括熱傳導(dǎo)、熱對(duì)流及熱輻射三種熱傳遞方式。(二二)、 ANSYS熱分析的特點(diǎn)熱分析的特點(diǎn)如熱量獲取或損失、如熱量獲取或損失、熱梯度、熱流密度熱梯度、熱流密度(熱通量熱通量)等等內(nèi)燃機(jī)、渦輪機(jī)、內(nèi)燃機(jī)、渦輪機(jī)、換熱器、管路系統(tǒng)、換熱器、管路系統(tǒng)、電子元件等電子元件等還可以分析相

43、變、有內(nèi)熱源、接觸熱阻等問(wèn)題。還可以分析相變、有內(nèi)熱源、接觸熱阻等問(wèn)題。59(三三)、 ANSYS熱分析分類熱分析分類(四四)、 邊界條件初始條件邊界條件初始條件熱結(jié)構(gòu)耦合熱結(jié)構(gòu)耦合熱流體耦合熱流體耦合熱電耦合熱電耦合熱磁耦合熱磁耦合熱電磁結(jié)構(gòu)耦合等。熱電磁結(jié)構(gòu)耦合等。熱耦合分析熱耦合分析穩(wěn)態(tài)傳熱穩(wěn)態(tài)傳熱瞬態(tài)傳熱瞬態(tài)傳熱溫度溫度熱流率熱流率熱流密度熱流密度對(duì)流對(duì)流輻射輻射絕熱絕熱生熱。生熱。60(1)僅用于評(píng)估由于網(wǎng)格密度不夠帶來(lái)的誤差。僅用于評(píng)估由于網(wǎng)格密度不夠帶來(lái)的誤差。(2)僅適用于僅適用于SOLID或或SHELL的熱單元的熱單元(1DOF)。(3)基于單元邊界的熱流密度的不連續(xù)?；趩?/p>

44、元邊界的熱流密度的不連續(xù)。(4)僅對(duì)線性、穩(wěn)態(tài)熱分析有效。僅對(duì)線性、穩(wěn)態(tài)熱分析有效。(5)使用自適應(yīng)網(wǎng)格劃分可以對(duì)誤差進(jìn)行控制。使用自適應(yīng)網(wǎng)格劃分可以對(duì)誤差進(jìn)行控制。(五五)、 ANSYS熱分析誤差估計(jì)熱分析誤差估計(jì)612 穩(wěn)態(tài)熱分析穩(wěn)態(tài)熱分析(一一)、定義、定義通常在進(jìn)行瞬態(tài)熱分析以前，進(jìn)行穩(wěn)通常在進(jìn)行瞬態(tài)熱分析以前，進(jìn)行穩(wěn)態(tài)熱分析用于確定初始溫度分布。態(tài)熱分析用于確定初始溫度分布。穩(wěn)態(tài)熱分析用于研究穩(wěn)定的熱穩(wěn)態(tài)熱分析用于研究穩(wěn)定的熱載荷對(duì)系統(tǒng)或部件的影響。載荷對(duì)系統(tǒng)或部件的影響。確定由于穩(wěn)定的熱載荷引起的溫度、確定由于穩(wěn)定的熱載荷引起的溫度、熱梯度、熱流率、熱流密度等參數(shù)。熱梯度、熱流率、

45、熱流密度等參數(shù)。62(二二)、 熱分析的單元熱分析的單元 LINK32 LINK32 二維二維2 2節(jié)點(diǎn)熱傳導(dǎo)單元節(jié)點(diǎn)熱傳導(dǎo)單元 LINK33 LINK33 三維三維2 2節(jié)點(diǎn)熱傳導(dǎo)單元節(jié)點(diǎn)熱傳導(dǎo)單元 LINK34 2LINK34 2節(jié)點(diǎn)熱對(duì)流單元節(jié)點(diǎn)熱對(duì)流單元 LINK31 2LINK31 2節(jié)點(diǎn)熱輻射單元節(jié)點(diǎn)熱輻射單元 PLANE55 4PLANE55 4節(jié)點(diǎn)四邊形單元節(jié)點(diǎn)四邊形單元 PLANE77 8PLANE77 8節(jié)點(diǎn)四邊形單元節(jié)點(diǎn)四邊形單元二維實(shí)體二維實(shí)體 PLANE35 6PLANE35 6節(jié)點(diǎn)三角形單元節(jié)點(diǎn)三角形單元 PLANE75 4PLANE75 4節(jié)點(diǎn)軸對(duì)稱單元節(jié)點(diǎn)軸對(duì)稱

46、單元 PLANE78 8PLANE78 8節(jié)點(diǎn)軸對(duì)稱單元節(jié)點(diǎn)軸對(duì)稱單元 SOLID87 10SOLID87 10節(jié)點(diǎn)四面體單元節(jié)點(diǎn)四面體單元三維實(shí)體三維實(shí)體 SOLID70 8SOLID70 8節(jié)點(diǎn)六面體單元節(jié)點(diǎn)六面體單元 SOLID90 20SOLID90 20節(jié)點(diǎn)六面體單元節(jié)點(diǎn)六面體單元 殼殼 SHELL57 4SHELL57 4節(jié)點(diǎn)節(jié)點(diǎn) 點(diǎn)點(diǎn) MASS71MASS71線性線性1414種專用單元種專用單元63(三三)、 ANSYS穩(wěn)態(tài)熱分析的基本過(guò)程穩(wěn)態(tài)熱分析的基本過(guò)程前處理，建模；前處理，建模；求解，施加載荷計(jì)算；求解，施加載荷計(jì)算；后處理，查看結(jié)果。后處理，查看結(jié)果。641)建模建模(

47、1)(1)確定確定jobnamejobname，titletitle，unitunit。(2)(2)進(jìn)入進(jìn)入PREP7PREP7前處理，定義單元類型，設(shè)定單元選項(xiàng)。前處理，定義單元類型，設(shè)定單元選項(xiàng)。(3)(3)定義單元實(shí)常數(shù)。定義單元實(shí)常數(shù)。(4)(4)定義材料熱性能參數(shù)，對(duì)于穩(wěn)態(tài)傳熱，一般只需定義定義材料熱性能參數(shù)，對(duì)于穩(wěn)態(tài)傳熱，一般只需定義 導(dǎo)熱系數(shù)，它可以是恒定的，也可以隨溫度變化。導(dǎo)熱系數(shù)，它可以是恒定的，也可以隨溫度變化。(5)(5)創(chuàng)建幾何模型并劃分網(wǎng)格。創(chuàng)建幾何模型并劃分網(wǎng)格。2)施加載荷計(jì)算施加載荷計(jì)算直接在實(shí)體模型或單元直接在實(shí)體模型或單元模型上施加五種載荷模型上施加五種載

48、荷( (邊邊界條件界條件) )1)1)恒定的溫度恒定的溫度( (自由度約束自由度約束) )2)2)熱流率熱流率( (節(jié)點(diǎn)集中載荷節(jié)點(diǎn)集中載荷) )3)3)對(duì)流對(duì)流( (面載荷面載荷) )4)4)熱流密度熱流密度( (一種面載荷一種面載荷) )5)5)生熱率生熱率( (體載荷體載荷) )653)后處理后處理基本數(shù)據(jù)：節(jié)點(diǎn)溫度基本數(shù)據(jù)：節(jié)點(diǎn)溫度導(dǎo)出數(shù)據(jù)：導(dǎo)出數(shù)據(jù)：節(jié)點(diǎn)及單元的熱流密度；節(jié)點(diǎn)及單元的熱流密度；節(jié)點(diǎn)及單元的熱梯度；節(jié)點(diǎn)及單元的熱梯度；單元熱流率；單元熱流率；節(jié)點(diǎn)的反作用熱流率；節(jié)點(diǎn)的反作用熱流率；其它。其它。ANSYSANSYS將熱分析的結(jié)果寫入將熱分析的結(jié)果寫入* *.rth.rth

49、文件中文件中66四、穩(wěn)態(tài)熱分析實(shí)例四、穩(wěn)態(tài)熱分析實(shí)例 某一潛水艇可以簡(jiǎn)化為一圓筒，它由三層組成，最外面一層為不銹鋼，中某一潛水艇可以簡(jiǎn)化為一圓筒，它由三層組成，最外面一層為不銹鋼，中間為玻纖隔熱層，最里面為鋁層，筒內(nèi)為空氣，筒外為海水，求內(nèi)外壁面間為玻纖隔熱層，最里面為鋁層，筒內(nèi)為空氣，筒外為海水，求內(nèi)外壁面溫度及溫度分布。溫度及溫度分布。幾何參數(shù)筒外徑幾何參數(shù)筒外徑 30 ft30 ft 總壁厚總壁厚 2 in2 in 不銹鋼層壁厚不銹鋼層壁厚 0.75 in0.75 in 玻纖層壁厚玻纖層壁厚 1 in1 in 鋁層壁厚鋁層壁厚 0.25 in0.25 in 筒長(zhǎng)筒長(zhǎng) 200 ft200

50、ft導(dǎo)熱系數(shù)不銹鋼導(dǎo)熱系數(shù)不銹鋼 8.27 Btu8.27 Btuh hftft0 0F F 玻纖玻纖 0.028 Btu0.028 Btuh hftft0 0F F 鋁鋁 117.4 Btu117.4 Btuh.fth.ft0 0F F邊界條件空氣溫度邊界條件空氣溫度 70 F70 F 海水溫度海水溫度 44.5 F44.5 F 空氣對(duì)流系數(shù)空氣對(duì)流系數(shù) 2.5 Btu/h2.5 Btu/hftft2 2. .0 0F F 海水對(duì)流系數(shù)海水對(duì)流系數(shù) 80 Btu/h80 Btu/hftft2 2. .0 0F F鋁玻璃纖維不銹鋼0.50.5空氣空氣海水海水13/41/4R15英尺英尺67鋁玻

51、璃纖維不銹鋼0.50.5空氣空氣海水海水13/41/4R15英尺英尺68pcirc,Ro, Rss,-0.5,0.5 pcirc,Rss, Rins,-0.5,0.5pcirc,Rins, Ral, -0.5,0.5aglue, allnumcmp, arealesize,1,16lesize,4,4lesize,14,5lesise,16,2eshape,2mat,1amesh,1mat,2amesh,2mat,3amesh,31. 創(chuàng)建幾創(chuàng)建幾何模型何模型69pcirc,Ro, Rss,-0.5,0.5pcirc,Rss, Rins,-0.5,0.5pcirc,Rins, Ral, -0.5,0.5aglue, allnumcmp, arealesize,1,16lesize,4,4lesize,14,5lesise,16,2eshape,2mat,1amesh,1mat,2amesh,2mat,3amesh,32. 設(shè)定劃分設(shè)定劃分網(wǎng)格密度網(wǎng)格密度11416470pcirc,Ro, Rss,-0.5,0.5pcirc,Rss, Rins,-0.5,0.5pcirc,Rins, Ral, -0.5,0.5aglue, allnumcmp, arealesize,1,16lesize,4,4lesize,14,5lesise,16,2eshape,2mat,1a