計算力學的俠義定義通常是指結(jié)構(gòu)力學，在這里我們拓展一下概念把流體力學、結(jié)構(gòu)力學、電磁/力學統(tǒng)稱為計算力學。這三者之間的區(qū)別是控制方程：流體力學對應Navier-Stokes方程組，結(jié)構(gòu)力學對應三大方程組，電磁/力學對應Maxwell方程組。計算力學的本質(zhì)是使用離散的方法把連續(xù)的偏微分方程組轉(zhuǎn)變成離散點上的代數(shù)方程組，進行求解，從而獲得目標結(jié)果的空間分布。

我們常見的仿真軟件也主要集中在上述三組方程的領域，本文討論三大領域的計算技術。

有限三兄弟：FD/FV/FEM

有限差分FD、有限體積FV和有限元素FEM是三種最常見的離散方法。

FD/FV多用在流體方程的求解，主要差別在于求解的是微分形式還是積分形式的方程，格式本身的表現(xiàn)形式是一樣的。

有限元素方法在結(jié)構(gòu)/電磁力學中應用較多。有限元素方法的本質(zhì)是選擇特定的基函數(shù)使其組合在離散點上滿足離散方程（類似于多項式逼近的思想）。

格式三姊妹：離散/精度/格式

國家統(tǒng)計局每個月都會公布國名經(jīng)濟的若干指標包括CPI，PPI，以供人們判斷當前經(jīng)濟形勢。毫無疑問統(tǒng)計局不會統(tǒng)計每一個行業(yè)、每一家企業(yè)的數(shù)據(jù)，他們會根據(jù)經(jīng)驗選用少量有代表意義的商品、企業(yè)、行業(yè)的數(shù)據(jù)進行統(tǒng)計平均。這個過程叫采樣。采樣的含義就是當不需要了解全部樣本空間信息的時候，通過采樣了解部分樣本的信息，并假設這種數(shù)據(jù)代表了全部樣本空間的特征。

離散也是一種采樣。它在物理域中選用部分有特征的點，在這些點上對微分方程對應的離散方程進行解算，獲得這些點上的物理數(shù)據(jù)。

離散是由兩步完成的：空間網(wǎng)格劃分和差分近似。離散的本質(zhì)是近似，再進一步說就是把所有的一階導數(shù)/二階導數(shù)/三階導數(shù)，統(tǒng)統(tǒng)用1階/2階/高階精度的近似來表示。導數(shù)的階和近似的階不是一個概念。一階導數(shù)也可以有2階精度。這個所謂的幾階精度就是把導數(shù)用泰勒展開后獲得的截斷誤差的階數(shù)。

為了達到不同的精度，就用不同的方法來做近似，這種方法就是所謂的格式?？梢择v點求解的格式叫顯式格式。必須全場統(tǒng)一求解的叫隱格式。歷史上數(shù)學家們發(fā)明了Range-Kutta，Jameson，MacCormack，Harten等一些列重要的離散格式。

初邊二將：邊界/初始條件

微分方程在邊界上需要給定恰當?shù)倪吔鐥l件，以滿足方程解的適定性。邊界條件可以直接給定物理量，或者物理量的梯度。

初始條件通常針對非定常分析，需要給定初始時刻的全場物理量分布，據(jù)此進行參數(shù)的物理演化，獲得其他時刻的分布。

九九歸一：收斂準則

計算力學各個領域都是采用迭代求解的方法進行的，即便是定常物理場，也是逐步迭代獲得最終結(jié)果的。這就需要對何時結(jié)束計算進行判斷。收斂準則通常有兩類判別方法：守恒量判斷和變化量判斷。

守恒量判斷比如對于一段封閉的管路，流動定常，那么計算收斂后進出口截面的流量應該是相同或者只存在有極小差值，對于可以容忍的極小差值的定義就是所謂收斂準則。

變化量的判斷是用迭代前后兩部變化量的絕對差值與變化量之比來進行，當比值的絕對值小于某個限值即可判定為收斂。

到現(xiàn)在為止共有5篇文章講解了仿真知識架構(gòu)分別是：方程、直接經(jīng)驗（2篇）、物理想象和計算力學。

為了驗證我們的知識框架結(jié)構(gòu)，讓我們來考慮一艘飛船的模擬。

端午節(jié)那天的天宮一號飛船的熱環(huán)境模擬

物理想象：和大多數(shù)人一樣，我們都沒有飛船上的直接生活經(jīng)歷，所以需要依靠物理想象重構(gòu)飛船的熱仿真計算。由于飛船體積較小，返回艙需要經(jīng)歷燒蝕，軌道艙需要經(jīng)歷太陽輻射，飛船結(jié)構(gòu)中一定存在絕熱層，減少艙內(nèi)外的熱量交換。

方程：熱環(huán)境模擬需要考慮到飛船受到的輻射和導熱，也需要考慮飛船內(nèi)部的乘員艙的空氣溫度分布，所以單獨使用結(jié)構(gòu)力學的導熱算法肯定不行，我們選用基于流體力學的求解器，比如Fluent。

直接經(jīng)驗：飛船主要的熱源來源于太陽，飛船繞地球的周期為1.5小時，所以其相對于太陽的角度是快速變化的，熱量也是快速變化，這意味著需要進行非定常計算。

計算規(guī)模：單艙飛船的流固耦合導熱需要數(shù)千萬網(wǎng)格。

計算設定：構(gòu)成飛船的各種物質(zhì)的比熱、導熱系數(shù)、彈性模量、泊松比?？臻g站空氣的構(gòu)成。需要考慮太陽輻射的吸收/反射。需要考慮人體和設備的發(fā)熱。需要考慮空調(diào)系統(tǒng)的進風參數(shù)（無重力情況下風機的風量/流量曲線）。計算艙室和其他艙室連接部分的溫度條件。

審核編輯：劉清