求解微分方程,特別是橢圓型邊值問題的一種離散化方法,其基礎是變分原理和剖分逼近。有限元方法是傳統的裏茨-加廖金方法的發展,並融會瞭差分法的優點,處理上統一,適應能力強,已廣泛應用於科學與工程中龐大複雜的計算問題。

  作為有限元方法出發點的變分原理,是表達物理基本定律的一種普遍形式。其表述可概括如下:給出一個依賴物理狀態v的變數J(v)(v是函數,J(v)在數學上稱為泛函),同時給出J(v)的容許函數集V,即一切可能的物理狀態,則真實的狀態是V中使J(v)達到極小值的函數。剖分逼近是有限元離散化的手段,把問題的整體(即求解域)剖分為有限個基本塊,稱為“單元”,然後通過單元上的插值逼近,得到一個結構簡單的函數集,稱為“有限元空間”,它一般是容許函數集V的子集或有某種聯系。有限元方法就是在這個有限元空間中尋找J(v)的極小解作為近似解。

  典型問題 為具體說明有限元方法,討論二維有界域Ω上的橢圓型方程

,   (1)

變系數 β表示介質不均勻。物理學中許多平衡態或定常態問題都可歸結為這個典型方程。與方程(1)相配的有如下三類邊界條件:

  第一類:

  第二類:

  第三類:

。這裡的 φg及α均為定義在邊界д Ω上的已知函數, 表示外法向導數,第二類邊界條件是第三類當 α=0時的特例。

  為說明有限元方法能統一處理復雜的情況,假定討論的問題是混合邊值,並且介質有間斷,即дΩ分成Г0和Г1兩部分,分別有邊界條件

,   (2)

,   (3)

βxy)有間斷線,把 Ω分為 Ω -Ω +兩部分,在間斷線上微分方程(1)無定義,而代之以接觸條件

,   (4)

表示間斷線上分別指向 Ω +Ω -的法向導數。

  變分原理 與微分方程(1)及附加條件(2)、(3)、(4)的邊值問題相對應的是物理學中的極小能量原理。構造“能量積分”

並取 J( v)的容許函數集 V為一切滿足邊界條件(2)且一階偏導數平方可積的函數,則使 J( v)達到極小值的 u,即

,   (6)

也必滿足方程(1)及(2)、(3)、(4)。事實上,極小能量原理之類的變分原理是物理問題的原始形式,微分方程是數學推導的結果。在變分問題中,隻有邊界條件(2)是強加到容許函數集上的,邊界條件(3)及間斷介質的接觸條件(4)都是極小解 u自然滿足的,這種情況有利於離散化的統一處理。

  剖分逼近 幾何剖分的基本單元可取為三角形、矩形、四邊形、曲邊形等等,其中三角形最基本常用。

  假定問題的求解區域為多邊形,介質間斷線為折線,作三角剖分如圖

所示。在剖分中需註意介質間斷線與某些三角形的邊重合,不同類邊界條件的交點與某些三角形的頂點重合。單元的頂點稱為網格結點,在д Ω上稱邊界結點,在 Ω內稱內結點。

  幾何剖分之後考慮插值逼近。對三角形單元最簡單的是線性插值,即利用每個單元Δk三頂點的函數值確定線性函數αkx+bkyk的三個系數。把所有單元{Δk}確定的{αkx+bkyk}合在一起,就得到Ω上的一個分片線性插值函數。Г0上的邊界結點取值為零的分片線性插值函數都屬於問題(5)、(6)的容許函數集V,全體這樣的函數構成一個有限維線性空間

,稱為有限元空間。假定內結點和Г 1上的邊界結點共有 N個,以 p j( j=1,…, N)表示,則 的維數就是 N oφ i表示 中滿足條件

   (7)

的成員,則{ φ i}構成線性空間 的一組基。 中任意函數 v,都可表為

,   (8)

V j是結點 p j上的函數值 v( p j)。

  單元上的插值方式除瞭用一次函數外,還可以用二次、三次或更高次的多項式,也可用非多項式函數。插值數據除瞭用函數值的拉格朗日型外,還可以是包括導數的埃爾米特型插值。種種的幾何剖分加上種種的插值方式,就產生眾多形式的有限元空間,使有限元方法可有眾多的選擇。

  有限元的離散化 有限元離散化的出發點是與微分方程等價的變分問題。對於典型問題來說,就是從(5)、(6)出發,用剖分逼近的方法構造有限元空間

(也稱試探函數空間),然後求泛函 J( v)在 中的極小解ũ 作為近似解,即ũ滿足

,   (9)

把(8)的表達公式 代入(5)中的 J( v),得

,   (10)

式中

,   (11)

,   (12)

把(9)的極小解表為 ,則( U 1U 2,…, U N)使二次函數(10)達到極小,由微分學知滿足線性方程組

。   (13)

方程組(13)來自正定二次函數的極小解問題,故系數矩陣一定對稱正定。由於基函數 φ i隻在以 p i為頂點的單元上不為零,故系數 α ij= 隻當結點 p ip j連成三角形一邊時才不為零。系數矩陣這種稀疏性質,加上對稱正定,對方程的求解很有利。

  系數

以及自由項 的實際計算,通常按所謂單元分析與總體合成的方式進行。即逐個分析 Ω內的單元和Г 1上的單元邊對有關的 α izf i的貢獻,然後往上迭加。當 Ω內所有單元及Г 1上所有單元邊都分析之後,方程組(13)的系數矩陣及自由項也就合成出來。間斷介質的影響反映在單元分析中被積函數的 βΩ +Ω -取不同的表達式。單元分析通常都采用某種數值積分公式計算。

  從虛功原理出發的離散化 微分方程邊值問題(1)、(2)、(3)、(4)的解u還同時滿足:對容許函數集V中任一函數v,成立

,   (14)

這裡 α( uv)及 F( v)即表達式(11)、(12)。在物理學中,方程(14)是另一變分原理的數學形式,稱為虛功原理或虛位移原理。有限元方法更一般的形式是從虛功方程(14)出發用剖分插值的方式構造一個試探函數空間 ,並同時構造一個檢驗函數空間 ψ;在 中尋找近似解 ũ,使之對 ψ中的任一函數 ψ,成立

,   (15)

當選取 ψ 相同時,(15)中的 ψ可選為基函數 φ i,同時用 代入,就得到方程組(13)。

  對於非自共軛橢圓算子L,微分方程邊值問題Lu=f不存在等價的極小值問題,但這時仍可建立虛功方程(14),其中α(uv)=(Luv)F)=(fv),(·,·)表示L2(Ω)的內積。因此,有限元方法仍然有效。

  從極小能量原理出發進行離散化又常稱為裡茨法,從虛功原理出發稱為加廖金法。後者是前者的推廣。

  評價 傳統的裡茨-加廖金方法,采取解析函數作為試探函數,不能滿足任意多邊形區域的邊界條件,也不適應間斷介質的要求,對現在的典型例子無能為力。差分方法雖然能夠對付,但由於它對方程(1)及條件(2)、(3)、(4)在處理上不統一,在計算效果及理論分析兩方面都帶來不利。有限元方法正好對這兩者揚長避短,一方面保持瞭裡茨-加廖金方法從變分原理出發的優點,在提法上有極大的概括性,給離散化帶來統一處理的方便;另一方面又吸收瞭差分法剖分逼近的優點,能靈活適應各種幾何形狀和間斷介質等復雜情況。有限元方法除瞭解題效能高強外,還有牢靠的理論基礎,是計算數學理論一大成就。

  回顧與展望 有限元方法在中國與西方從不同的實踐背景,沿著不同的學術道路、各自獨立平行地發展起來。在西方,有限元思想在R.庫朗1943年的一篇論文中明確地提出過,但一直沒有受到重視。20世紀50年代中期,歐美工程界J.H.阿吉裡斯、R.W.克拉夫等以航空工程為背景,在結構分析和矩陣方法基礎上提出瞭結構有限元的雛形。60年代初期,引進連續體的單元剖分;60年代中期,逐漸明確有限元法是變分原理加剖分逼近的思想。1968年,西方數學傢對有限元法進行數學的理論分析,開始瞭有限元法在計算數學中的黃金時代。

  在中國,60年代初期,馮康、黃鴻慈等結合解決一系列大型水壩建設的應力分析問題,開展瞭橢圓型邊值問題數值解的系統研究,為克服問題傳統提法中的幾何復雜性和材料復雜性,把能量法與差分法結合在一起,於1964年建立瞭求解橢圓型邊值問題一套普遍有效的方法,命名為基於變分原理的差分方法,即通稱的有限元方法。與此同時,建立瞭方法的數學理論基礎。而後20年中,周天孝、唐立民對混合元擬協調元的發展,應隆安等對無限元的發展,馮康等對邊界有限元的發展,石鐘慈對非協調元的發展,林群對有限元外推理論的發展,都作瞭重要貢獻。

  有限元方法對於定常態問題的計算已經獲得公認的巨大成功,對不定常態問題也有良好開展。有限元方法是一個發展著的體系,在前述的基本原則下可有種種變化和發展,特別是可和其他方法結合起來,進一步解決更困難更復雜的數學問題。

  

參考書目

 馮康、石鐘慈著:《彈性結構的數學理論》,科學出版社,北京,1981。

 G.斯特朗、G.J.菲克斯同著,崔俊芝、宮著銘譯:《有限元分析》,科學出版社,北京,1983。(G.Strang and G.J.Fix,An Analysis of the Finite Element Method,Prentice-Hall,Englewood Cliffs,New Jersey,1973.)

 P.G.Ciarlet,The Finite Element Method for Elliptic Problems,North-Holland,Amsterdam,1978.

 O.C.Zienkiewicz,The Finite Element Method,3rded.,McGraw-Hill,London,1977.