弱形式到底是什么?為什么有限元離不開它?

更多精彩內容,請關注“鋰電芯動”公眾號

弱形式到底是什么?為什么有限元離不開它?

在學習有限元方法(FEM)時,我們總會遇到一個看似抽象又略顯神秘的概念——微分方程的弱形式(Weak Form)

但教材往往“一上來就講推導”,很多人沒弄明白它到底是什么,就已經陷入了公式的旋渦中。

本篇文章,我們從頭梳理弱形式的本質,力圖回答三個關鍵問題:

  • 什么是“弱形式”?它和我們熟悉的微分方程形式有何不同?
  • 為什么有限元方法必須依賴弱形式才能成立?
  • 弱形式到底怎么推導出來?它的結構、意義和數值實現方式是什么?

什么是“弱形式”?

所謂弱形式,是對微分方程的一種積分等價重寫

在“強形式”中,微分方程需要在每一個點都被嚴格滿足,同時要求解函數具有較高的光滑性(例如二階連續可導)。但在實際工程問題中,物理條件復雜,解往往不夠光滑,這就給解析方法和數值求解帶來了巨大挑戰。

弱形式的出現,就是降低對解的光滑性(可微性)要求,它是許多強大的數值方法(尤其是有限元法 Finite Element Method, FEM)的理論基礎。有限元方法的本質,就是將一個偏微分方程問題轉化為一個代數系統進行求解。這個轉化過程的核心步驟,正是將強形式變成弱形式。

強形式 vs 弱形式

為了說明“強形式”和“弱形式”的區別,我們從一個最簡單的一維邊值問題出發,這樣更容易理解其背后的邏輯結構:

例子:一維泊松方程

考慮以下邊值問題: 弱形式到底是什么?為什么有限元離不開它?的圖1

這是一個非常基礎的微分方程,描述了例如一根受力桿的穩態變形問題,其中 u(x) 是位移, f(x) 是分布載荷。

什么是“強形式”?

強形式就是上面直接寫出來的這個微分方程,它要求:

  1. 函數 u(x) 至少是二階可導的(即 u ∈ C^2(0,1));
  2. 在定義域內的每一個點都嚴格滿足該微分方程;
  3. 邊界條件 u(0)=u(1)=0 也必須嚴格滿足。

這樣的要求非常“嚴格”,因此稱為“強形式”。

如何得到“弱形式”?——四個步驟

我們現在來將這個強形式轉換為弱形式,過程如下:

第一步:引入試函數(test function)

選擇一個試函數 v(x) ∈ V,其中 V = { v ∈ H^1(0,1) ∣ v(0)=v(1)=0 },即滿足同樣邊界條件的“光滑但不一定二階可導”的函數。

第二步:乘以試函數并積分(加權殘差法)

我們將方程兩邊都乘以 v(x),然后對整個定義域積分: 弱形式到底是什么?為什么有限元離不開它?的圖2

第三步:使用分部積分(積分轉移導數)

我們不希望對 u'' 求值,因為太“強”。于是使用積分分部公式:

弱形式到底是什么?為什么有限元離不開它?的圖3因為試函數 v(0) = v(1) = 0 ,邊界項為零。

第四步:得到弱形式表達

最終我們得到弱形式:

弱形式到底是什么?為什么有限元離不開它?的圖4?? 強 vs 弱:真正的差別在哪里?

項目 強形式 弱形式 函數要求

u ∈ C^2

u ∈ H^1(只需一階導數)

成立方式

每一點都要滿足

對所有試函數 v 滿足積分等式

數學處理

偏微分方程

積分方程

是否利于數值解

是(可直接構造有限維子空間)

弱形式更“寬容”,允許函數更粗糙;同時也更適合數值方法處理,如有限元、Galerkin方法等。

小結

這個例子展示了一個重要思想:

我們不是“簡化”微分方程,而是“換個角度理解”它。

從點的約束 → 積分意義的約束,正是有限元能夠成立的關鍵。

下一篇文章,我們將基于Comsol的一個官方案例進一步解釋如何將弱形式應用到仿真建模中。

?? 記得【關注+點贊+轉發】,不錯過任何一篇高能干貨!

登錄后免費查看全文
立即登錄
App下載
技術鄰APP
工程師必備
  • 項目客服
  • 培訓客服
  • 平臺客服

TOP

1