馮曉燕，賀 熙，游小龍，周林凱，單保平

（新余學院建筑工程學院，江西 新余 338004）

熱傳導具體表現(xiàn)為能量從物質(zhì)中的溫度高的一邊轉(zhuǎn)移到溫度低的一邊，從溫度高的物質(zhì)轉(zhuǎn)向溫度低的物質(zhì)，這一現(xiàn)象的產(chǎn)生是因為大量分子之間的相互撞擊引起的。目前熱傳導問題是科技研究領域中的熱點話題，常用的數(shù)值方法包括有限元法、比例邊界有限元法、有限差分法、擴展有限元法、有限體積法及無網(wǎng)格法等。

1 有限元法

當分析工程問題時經(jīng)常需要對各種微分方程描述的物理場進行求解，但由于許多微分方程都很難得到其解析解，這時可以將微分方程的問題區(qū)域進行離散化，即把復雜的問題區(qū)域分割為比較容易分析的有限元來表示，有限元之間通過節(jié)點相互連接，然后使用變分法或加權(quán)余數(shù)法，將偏微分方程轉(zhuǎn)化為代數(shù)方程組進行計算，整個問題區(qū)域的求解由對所有子區(qū)域進行獨立的運算結(jié)果組成。最終的結(jié)果并不是一個準確值，而是一個求極限之后的近似值。在實際問題中很難得到精準值的情況下，采用有限元法的可以得到一個精度較高的近似值，而且它能適用于各種復雜形狀，是一種較為成熟的工程分析手段。現(xiàn)如今互聯(lián)網(wǎng)技術(shù)發(fā)展迅速，涌現(xiàn)出了許多應用有限元法軟件，常見的有限元軟件有ANSYS、MSC、ABAQUS等公司的產(chǎn)品。隨著計算機技術(shù)的發(fā)展，有限元軟件與CAD軟件被結(jié)合使用，可以在CAD軟件上直接完成零部件的繪制后，軟件直接運用有限元法原理進行計算，如果不符合要求工程師可以直接在CAD軟件上對零部件進行修改，從而可輕松快捷地完成以前難以應對的復雜工程分析問題。有限元法在實踐中的應用十分普遍，例如在工程結(jié)構(gòu)、傳熱、流體力學、電磁力學等領域中，有限元分析越來越多地被運用于工程問題的仿真模擬，來求解真實的工程問題。熱傳導是一種在自然界常見的現(xiàn)象，如果存在溫度差，熱量就會由高溫的向低溫物體傳輸，就會存在熱量的增加和減少，從而引起溫度的變化。在分析熱傳導問題時，運用拉普拉斯方程描述穩(wěn)定溫度場的求解是工程中常遇到的問題，但很難求出其解析解，這時可以運用有限元法求解偏微分方程求得一個較為精準的近似解，首先利用變分原理將偏微分方程轉(zhuǎn)化為等價的泛函，假設單元上的場變量變化形式及插值函數(shù)或嘗試函數(shù)，尋找嘗試函數(shù)的系數(shù)-節(jié)點場變量，以使泛函取極值，從而得到有限元的計算格式計算出穩(wěn)定溫度場。魯治誠、李錄賢在復合材料熱傳導問題中提出了多尺度有限元法的應用，對于傳熱屬性突變的復合材料，主要關(guān)心復合材料宏觀尺度，忽略細觀尺度的求解精度時，取得單元形狀函數(shù)滿足問題的微分控制方程，將材料的傳熱突變特性反映在形狀函數(shù)上，可在較小的求解規(guī)模上，取得較高的求解精度。這種采用多尺度有限元法既提高了計算精度，又節(jié)省了計算時間。

2 比例邊界有限元法

比例邊界有限元法是新研究發(fā)展趨勢下的必然產(chǎn)物。該方法在求解熱傳導問題上，擺脫了陳舊的解析方式，結(jié)合多元解法的優(yōu)點，達到高效精準的目的。比例邊界有限元法只需要在求解區(qū)域的邊界使用有限元網(wǎng)格進行離散，無需基本解，便可以精確、快速求解熱傳導問題。該方法最早應用在固體力學領域，并且集合了有限元法和邊界元法的特點于一身。在使用比例邊界有限元法進行數(shù)值計算時，該法通過在計算域引入相似中心，建立含有徑向和環(huán)向的比例邊界坐標系，再經(jīng)過相似變換，物理變量轉(zhuǎn)變成域邊界點和相似中心相連的徑向坐標，偏微分控制方程轉(zhuǎn)變?yōu)橐浑A常微分方程。其在單元離散時同邊界元法一樣，無須在整個結(jié)構(gòu)內(nèi)部劃分單元，只需要在計算區(qū)域的邊界環(huán)向坐標方向上進行離散， 從而降低了一維的空間維度，不需要設立相應的基本解，使得計算量明顯減少，而且比例邊界有限元法的計算效率和模擬精度都非常高，這是由于其在輻射坐標上的解析性。該方法還有自動滿足邊界條件達到無限遠的一大優(yōu)勢。比例邊界有限元方法最初是為了求解無限域問題而提出的，但隨著其理論發(fā)展越來越完善，它的應用范圍也變得越來越廣，在各個方面都取得了不錯的成就。在彈性力學研究方向上，胡志強、林皋、王毅等基于Hamilton體系的彈性力學問題的比例邊界有限元方法研究內(nèi)容中，運用比例邊界有限元的離散方法，結(jié)合彈性力學的Hamilton對偶體系，獲得狀態(tài)空間的控制方程，并由此得到彈性靜力學的比例邊界有限元的控制方程，導出靜力邊界剛度矩陣并求解。在靜電場問題上，劉俊、林皋、王復明等學者在基于靜電場分析的比例邊界有限元法內(nèi)容中，由于該方法輻射坐標上具有解析性，故運用這種方法不僅有較高的模擬精度，并且無限遠的邊界條件可以自動滿足。通過以描述靜電場特性的控制方程，也就是拉普拉斯方程為起點，首先建立起笛卡爾坐標系統(tǒng)和比例邊界有限元坐標系統(tǒng)之間的轉(zhuǎn)換關(guān)系，再通過加權(quán)余量法推導出靜電場分析的比例邊界有限元方程，引入特征值問題對該方程進行求解得到電位計算公式以及電場計算公式。最后通過具體的算例對該數(shù)值方法與其他數(shù)值方法得到的結(jié)果進行對比，發(fā)現(xiàn)該方法的精度和計算效率都遠遠優(yōu)于其他方法。

3 有限差分法

在眾多的熱傳導問題數(shù)值分析的方法中，有限差分法是很經(jīng)典、發(fā)展很成熟的方法之一，也是解決此類問題的有效工具之一，直到現(xiàn)在仍然有廣泛的應用。有限差分法的原理是用差商去近似的代替微商，也就是用一個具體的函數(shù)值近似的去替換微分方程中的導數(shù)。其本質(zhì)上是將微分問題直接轉(zhuǎn)換為代數(shù)問題的近似數(shù)值。在數(shù)值分析時，把要分析的問題的解域劃分為差分網(wǎng)格，此網(wǎng)格則是由有限數(shù)量的離散點來構(gòu)成，把連續(xù)的解域用有限數(shù)量的網(wǎng)格節(jié)點所代替，網(wǎng)格節(jié)點也就是構(gòu)成網(wǎng)格的那些有限數(shù)量的離散點，通過一系列數(shù)學方法，用網(wǎng)格上的離散變量函數(shù)近似地替換掉連續(xù)變量函數(shù)，用差商近似地替換掉原方程和條件中的微商，即把網(wǎng)格節(jié)點上的具體函數(shù)值近似的替換掉控制方程中的導數(shù)，這樣就可以得到以網(wǎng)格節(jié)點上的函數(shù)值為未知量的代數(shù)方程或者方程組，稱為有限差分方程組，用它來近似地代替原微分方程和定解條件。通過解有限差分方程組，可以得到原問題在網(wǎng)格節(jié)點上的近似解，進一步地可利用插值等方法求得定解問題在整個區(qū)域上的近似解。張建霞、邵創(chuàng)新、田坤等基于有限差分法對低溫防護服進行優(yōu)化設計與仿真，學者們根據(jù)熱傳導微分方程進行了熱傳導建模后，發(fā)現(xiàn)該問題設置的邊界條件較為繁復，使用解析方法分析較為困難，于是決定使用有限差分法來求解這個熱傳導方程，最終在Matlab平臺上進行了計算與模擬驗證，得到了數(shù)值結(jié)果。同樣的，在林文生關(guān)于相變熱傳導問題的數(shù)值解及應用中，因為該問題難以得到解析解，所以選擇使用有限差分法來進行數(shù)據(jù)模擬分析，最后證實了有限差分法在熱傳導問題上的可行性?？偟膩碚f，有限差分法在解決有關(guān)熱傳導的數(shù)據(jù)分析等問題上運用是可行的，是解決此類問題的得力工具。然而有限差分法也有其自身的局限性，在王秦、雷鈞、谷巖等的基于時空聯(lián)合廣義有限差分法求解三維非穩(wěn)態(tài)熱傳導方程中提到，傳統(tǒng)的有限差分法過于依賴網(wǎng)格的劃分，離散差分的過程步驟量大且復雜，從而導致編程和計算難度比較大。為此其采用了有限差分法延伸出的新方法，即廣義有限差分法，最后通過對一些問題進行模擬和分析，驗證了此方法是可行而且準確的。綜上，在熱傳導問題數(shù)值分析中，有限差分法的表現(xiàn)是十分出色的，其優(yōu)點具體表現(xiàn)為原理簡單，應用成熟，且在各類問題上表現(xiàn)得很通用，甚至在無解析解等困難條件下仍然可以使用。其不足之處表現(xiàn)為依賴網(wǎng)格的建立和劃分，需要預估和評價誤差等。

4 有限體積法

目前，有限體積法是一種比較新的研究熱傳導問題的數(shù)值計算方法，其特點是既具有有限分差法的靈活性，又具有較高的精度和較好的便捷性。有限體積法在傳熱計算中起著重要作用，據(jù)統(tǒng)計，全球每年發(fā)表的相關(guān)計算傳熱問題的文章中約有60%使用有限體積法。有限體積法基于熱積分方程，可以使用儲能常數(shù)法和傅里葉定律直接從溫度場的任何有限區(qū)域分析中獲得。取溫度場中任何一個閉合面F，將周圍區(qū)域記錄為V，并根據(jù)能量保存方法對該區(qū)域的熱平衡進行分析，單位時間內(nèi)輸入和輸出區(qū)域之間的熱差與內(nèi)部熱源產(chǎn)生的熱量必須等于該區(qū)域中物體內(nèi)部能量的增加。在有限體積法中，由于使用了單元中溫差的函數(shù)，結(jié)果為近似值。由于體積循環(huán)的控制方式不同，得到了不同的線性方程和分析結(jié)果，控制體積選擇的效率直接取決于有限體積法計算的準確性。目前，控制矩形元素體積的最常用方法是圓形繪制方法。在有限元法中，溫度場是離散的，每個方程被分解成每個單元的貢獻，最終合成一系列線性方程，元素搜索法計算每個元素對單個方程的貢獻，并且有限體積法，即能量方程，可以圍繞每個節(jié)點進行求解和創(chuàng)建。受有限元法的啟發(fā)，可以使用相同的單位搜索方法，首先執(zhí)行單位分析，計算每個元素對每個方程的貢獻，然后執(zhí)行積分綜合。為了得到節(jié)點溫度未知的線性方程，必須為每個單元建立一個單位插值函數(shù)，單元內(nèi)的溫度分布必須表示為節(jié)點的溫度，有限體積法與有限元法完全相同。對于矩形單元，單元中的溫度分布是任意四個節(jié)點的溫度的線性函數(shù)，熱流的密度是恒定的。有限體積具有很多優(yōu)點，但有限體積方法在某些情況下還存在一些計算誤差，所以對有限體積方法的理解有待提高和細化。

5 無網(wǎng)格法

無網(wǎng)格法是近幾十年分析工程問題的重要數(shù)值分析方法，其在數(shù)值計算時不需要生成網(wǎng)格，只需尋找一些隨意分布的坐標點，然后將這些坐標點通過插值函數(shù)列出離散控制方程，從而可清晰地表達出各種復雜狀態(tài)的熱場。這種方法采用點的近似值，可以從根本上或一部分地消除網(wǎng)格，具有較強的適用性和靈活性。目前已經(jīng)提出了十幾種無網(wǎng)格法，如無單元伽遼金法、重構(gòu)核粒子法、無單元伽遼金比例邊界法等。但無網(wǎng)格法沒有近似的解析表達式，而且大部分都是以伽遼金原理為基礎，導致求解時計算量大、計算較為困難。最近新興的一種無網(wǎng)格方法—插值型無單元伽遼金比例邊界法是一種在改進的插值型移動最小二乘法框架下，結(jié)合了比例邊界法和無單元伽遼金法長處的半解析數(shù)值方法。這種方法通過引入比例邊界坐標系，只需在求解域的環(huán)向上進行數(shù)值離散，在徑向上采用解析的方法進行計算，處理工程結(jié)構(gòu)問題時擁有可觀的精度與效率，且具有一定的靈活性。這主要歸因于以下兩點：（1）引入改進的插值型移動最小二乘插值法，使得形成的形函數(shù)具有高階連續(xù)性，很大程度上提高了計算精度與收斂速度；（2）運用矩陣函數(shù)與Schur分解并行的求解技術(shù)，避免出現(xiàn)平行特征向量問題，保證了數(shù)值解的準確性與穩(wěn)定性。目前這種方法已被成功應用于穩(wěn)態(tài)熱傳導問題分析中，并通過具體的數(shù)值算例驗證了其高精度的特點。王峰、陳佳莉與陳燈紅等建立了以滑動Kriging插值為基礎的無單元伽遼金比例邊界法，并采用這種方法準確地分析了有限域及半無限域熱傳導問題。王曉飛、李興國、任紅萍提出了求解瞬態(tài)熱傳導問題的插值型復變量伽遼金法，為驗證此方法在二維熱傳導問題中的有效性，對兩個瞬態(tài)熱傳導的第一類邊值問題運用CVEFG方法進行求解。首先運用插值型復變量伽遼金法推導二維瞬態(tài)熱傳導問題的控制方程，然后對該控制方程進行伽遼金積分弱離散化，經(jīng)過一系列的流程實施，最后通過數(shù)值算例，對二維熱傳導問題進行求解，最終得出結(jié)論，對求解二維瞬態(tài)熱傳導的第一類邊值問題，運用插值型復變量伽遼金法具有較好的擬合效果?？偟膩碚f，在用無網(wǎng)格法在求解二維熱傳導問題時，具有一定的優(yōu)勢。

6 結(jié)語

有限元法、比例邊界有限元法、有限差分法、有限體積法及無網(wǎng)格法等數(shù)值方法都存在自身的優(yōu)勢與不足，在未來研究發(fā)展方向上，有必要不斷改進方法，完善解析路線減弱缺點帶來的不利影響甚至是消除缺點。其中插值型無單元伽遼金比例邊界法在解決穩(wěn)態(tài)熱傳導問題上相比其他方法具有一定的優(yōu)勢，說明了該方法在熱傳導問題分析上的可行性，故采用該方法進一步處理瞬態(tài)熱傳導問題將具有一定的理論意義和實際價值。