物理場的求解方法概述【轉(zhuǎn)發(fā)】

從哲學(xué)的角度講,結(jié)構(gòu)及其相互作用的矛盾構(gòu)成了世界發(fā)展的根本原因。從物理學(xué)的角度來說,結(jié)構(gòu)和相互作用被譯為物質(zhì)和場。所以,場就是物質(zhì)之間的相互作用。物理學(xué)上所談的場包括引力場、電磁場、核場和弱場。場的相互作用是通過交換量子而實現(xiàn)的,現(xiàn)在所發(fā)現(xiàn)的量子有:光子和8種核場量子。工程應(yīng)用所談的場的與上述概念略有不同,它主要從實用而非量子角度來談場,如結(jié)構(gòu)場、流場、聲場、濃度場、溫度場、靜電場、穩(wěn)恒磁場和電磁場等。

要研究場,必須首先建立其數(shù)學(xué)模型。雖然對各個場的研究劃歸到不同的學(xué)科,但是所有學(xué)科基本上采取了相似的建模方法:微元分析法。該法從研究對象中選取一個微元(通常是無限小的立方體)為對象,列出微元6個面上和體內(nèi)的物理量,根據(jù)能量守恒定律、質(zhì)量守恒定律、動量定理等學(xué)科內(nèi)的特定規(guī)律列方程,得到一個或一個以上的偏微分方程。由于時間和空間是場的存在形式,因此該方程一般是以時間和空間為自變量,以所研究場變量為因變量的偏微分方程。一般而言,在固體力學(xué)領(lǐng)域,得到的是平衡方程(或動力方程)、物理方程和幾何方程;在流體力學(xué)領(lǐng)域,是質(zhì)量守恒方程、動量方程和能量方程;在傳熱學(xué)領(lǐng)域,是導(dǎo)熱方程;在電磁場領(lǐng)域,是麥克斯韋方程組。

對于這類問題主要存在四種解法:解析法、半解析法、數(shù)值法和實驗法。解析法十分精確可靠,但是只能求解比較簡單的問題;數(shù)值法充分利用計算機(jī)的特點(diǎn),對一個問題從空間和時間上進(jìn)行離散,理論上可以求解任意復(fù)雜的問題,但是計算量很大;半解析法一般對控制方程進(jìn)行近似,再采用數(shù)值分析方法得到控制方程中的待定系數(shù),比數(shù)值法求解速度更快,不過只針對某些簡單的問題有效;試驗法主要用于對象復(fù)雜的灰箱或黑箱系統(tǒng),是理論分析和數(shù)值計算的基礎(chǔ),但是花費(fèi)昂貴,周期長而且受到模型尺寸、測量精度等的制約。對于邊界條件或形體復(fù)雜的實際工程領(lǐng)域,用得比較多是數(shù)值法。

數(shù)值法又分為兩類:有網(wǎng)格法和無網(wǎng)格法。根據(jù)網(wǎng)格劃分的方式,有網(wǎng)格法又進(jìn)一步分為加權(quán)余量法、有限元法、邊界元法、有限差分法和有限體積法?；诰W(wǎng)格的方法其數(shù)學(xué)原理各不相同,但是求解過程都分為兩個階段:離散化階段和代數(shù)方程組的求解階段。在離散化階段,用網(wǎng)格線將連續(xù)的定解域化分為有限個結(jié)點(diǎn)集,選取適當(dāng)?shù)耐緩綄⑽⒎址匠碳捌涠ń鈼l件轉(zhuǎn)化為網(wǎng)格結(jié)點(diǎn)上相應(yīng)的代數(shù)方程組,即建立代數(shù)方程組。然后使用計算機(jī)求解代數(shù)方程組,得到結(jié)點(diǎn)上的離散近似解。結(jié)點(diǎn)之間的近似解,一般認(rèn)為光滑變化,原則上可以應(yīng)用插值方法確定,從而得到定解問題在整個定解區(qū)域上的近似解。

有限差分法是最早采用的數(shù)值方法。它將求解區(qū)域分為矩形或正交曲線網(wǎng)格,在網(wǎng)格線交點(diǎn)(即結(jié)點(diǎn)上),將控制方程中的每一個微商用差商來代替,從而將連續(xù)函數(shù)的微分方程離散為網(wǎng)格結(jié)點(diǎn)上定義的差分方程,每個方程中包含了本結(jié)點(diǎn)及其附近一些結(jié)點(diǎn)上的待求函數(shù)值,通過求解這些代數(shù)方程就可獲得所需的數(shù)值解。

有限元法是目前發(fā)展最為成熟的一種數(shù)值方法。該方法將求解域剖分成相連結(jié)又互不重疊的具有一定規(guī)則幾何形狀的有限個子區(qū)域,這些子區(qū)域稱為單元,單元之間以結(jié)點(diǎn)相聯(lián)結(jié)。函數(shù)值被定義在結(jié)點(diǎn)上,在單元中選擇插值函數(shù),以結(jié)點(diǎn)函數(shù)值與基函數(shù)的乘積的線性組合成單元的近似解來逼近單元中的真解。利用某種方法建立單元方程,然后組裝成為總體的有限元方程,接著考慮邊界條件進(jìn)行求解。

邊界元法是在經(jīng)典積分方程和有限元法基礎(chǔ)上發(fā)展起來的求解微分方程的數(shù)值方法。其基本思想是:將微分方程相應(yīng)的基本解作為權(quán)函數(shù),應(yīng)用加權(quán)余量法并應(yīng)用格林函數(shù)導(dǎo)數(shù)聯(lián)系解域中待求函數(shù)值與邊界上的函數(shù)值與法向?qū)?shù)值之間關(guān)系的積分方程,令積分方程在邊界上成立獲得邊界積分方程,然后把邊界離散化,建立邊界元代數(shù)方程組,求解后可獲得邊界上全部結(jié)點(diǎn)的函數(shù)值和法向?qū)?shù)值,將全部邊界值代入積分方程中,即可獲得內(nèi)點(diǎn)函數(shù)值。

有限體積法又稱為控制體積法,其基本思路是:將計算區(qū)域化分為網(wǎng)格,并使每個網(wǎng)格點(diǎn)周圍有一個互不重疊的控制體積,將待解微分方程對每一個控制體積積分,從而得出一組離散方程。然后求解此離散方程,得到網(wǎng)格點(diǎn)上的因變量的值。

隨著計算機(jī)技術(shù)的發(fā)展和計算機(jī)的普及,上述基于網(wǎng)格的方法在過去幾十年內(nèi)得到了長足的發(fā)展,但是這些基于網(wǎng)格的方法時在使用時也遇到了一些問題,如結(jié)構(gòu)大變形問題,對結(jié)構(gòu)破壞過程的仿真,對新材料性能的分析和模擬,應(yīng)力高梯度和瞬態(tài)動力問題等。而且,基于網(wǎng)格的方法如果在計算過程中出現(xiàn)網(wǎng)格畸變,則導(dǎo)致計算失效。為了處理大變形或隨時間變化的不連續(xù)性,需要重新進(jìn)行網(wǎng)格重構(gòu)。為了解決這種問題,出現(xiàn)了無網(wǎng)格數(shù)值計算方法(Meshless Method),該方法拋棄了網(wǎng)格,直接基于結(jié)點(diǎn)形成近似,在特定領(lǐng)域內(nèi)具有重要的研究價值和應(yīng)用前景。

目前已經(jīng)提出十余種無網(wǎng)格方法,各種無網(wǎng)格方法之間的區(qū)別主要在于所使用的試函數(shù)(如移動最小二乘近似,重構(gòu)核函數(shù)近似,單位分解法,徑向基函數(shù)法,點(diǎn)插值法等)和微分方程的等效形式(如迦遼金法,配點(diǎn)法,最小二乘法等)。建立近似函數(shù)時不借助于網(wǎng)格,基于函數(shù)逼近近似而非插值是無網(wǎng)格法有限元法的主要區(qū)別。采用定義在離散結(jié)點(diǎn)上具有緊支特性的函數(shù)來構(gòu)造近似函數(shù),而不用定義在全域上的級數(shù)展開形式是無網(wǎng)格法與經(jīng)典加權(quán)殘量法的主要區(qū)別。

轉(zhuǎn)自公眾號:宋博士 ANSYS學(xué)習(xí)與應(yīng)用