有限元方法三大分支的深度对比 —— 结构有限元 · 固体力学有限元 · 多物理场有限元
===================================================================================

.. raw:: html

   <style>
   body { font-family: "Noto Sans SC", "Microsoft YaHei", sans-serif; line-height: 1.8; }
   .section-title { color: #1a1a2e; border-bottom: 3px solid #e94560; padding-bottom: 8px; }
   .highlight-box { background: #f0f4ff; border-left: 4px solid #4361ee; padding: 16px 20px; margin: 16px 0; border-radius: 0 8px 8px 0; }
   .warning-box { background: #fff3e0; border-left: 4px solid #e94560; padding: 16px 20px; margin: 16px 0; border-radius: 0 8px 8px 0; }
   .success-box { background: #e8f5e9; border-left: 4px solid #2e7d32; padding: 16px 20px; margin: 16px 0; border-radius: 0 8px 8px 0; }
   .info-box { background: #e3f2fd; border-left: 4px solid #1565c0; padding: 16px 20px; margin: 16px 0; border-radius: 0 8px 8px 0; }
   table { width: 100%; border-collapse: collapse; margin: 16px 0; font-size: 0.9em; }
   th { background: #1a1a2e; color: white; padding: 12px; text-align: center; }
   td { padding: 10px 12px; border: 1px solid #ddd; }
   tr:nth-child(even) { background: #f8f9fa; }
   tr:hover { background: #e8eaf6; }
   .tag { display: inline-block; padding: 2px 8px; border-radius: 4px; font-size: 0.85em; margin: 0 2px; }
   .tag-green { background: #e8f5e9; color: #2e7d32; }
   .tag-orange { background: #fff3e0; color: #e65100; }
   .tag-red { background: #ffebee; color: #c62828; }
   .tag-blue { background: #e3f2fd; color: #1565c0; }
   .tag-purple { background: #f3e5f5; color: #6a1b9a; }
   kbd { background: #f0f0f0; padding: 2px 6px; border: 1px solid #ccc; border-radius: 3px; font-family: monospace; }
   </style>

   <meta name="description" content="有限元方法深度对比研究：结构有限元(Abaqus/CalculiX) vs 固体力学有限元 vs 多物理场有限元(MFEM/Kratos)。从型函数、单元构造、自由度、连续性等角度深入对比。">

.. contents:: 目录
   :depth: 3
   :local:

写在前面
--------

本文针对一个困扰作者多年的问题：同样是有限元，为什么"结构有限元"（如 Abaqus、CalculiX）、"固体力学有限元"（连续介质有限元）、和"多物理场有限元"（如 MFEM、Kratos、deal.II）感觉像是三套不同的东西？它们的型函数（shape functions）是否相同？单元构造是否相同？如果不同，本质区别在哪里？

作者已有固体力学有限元、结构有限元、物质点法（MPM）、有限体积法（FVM）的开发经验。本文试图站在一个"写过代码"的工程师视角，用尽可能精确的语言回答上述问题。

本文参考了 Bathe (2014)、Zienkiewicz (2013)、Hughes (1987/2000)、Monk (2003)、Brenner & Scott (2008) 等经典著作，以及 MFEM、Kratos、CalculiX 的官方文档。

---

.. raw:: html

   <h2 class="section-title">一、三个FEM分支的快速定义</h2>

在深入对比之前，先用一句话定义这三个东西：

**① 结构有限元（Structural FEM）**

以 Abaqus、CalculiX、ANSYS Mechanical 为代表。核心特征是**维度缩减**：用梁（1D）、板壳（2D）模型近似三维连续体。引入了**转动自由度**（rotational DOFs）。单元类型包括：杆（truss）、梁（beam）、板（plate）、壳（shell）、以及各种连接/接触单元。理论基于 Bernoulli-Euler 梁、Kirchhoff/Mindlin 板壳理论。

**② 固体力学有限元（Solid Continuum FEM）**

处理完整的三维连续体。单元类型为：六面体（hexahedron）、四面体（tetrahedron）、五面体等。每个节点仅三个平动自由度。理论基于连续介质力学（continuum mechanics），无维度缩减。这是最"纯粹"的有限元——你写第一行 FEM 代码时通常就是这种。

**③ 多物理场有限元（Multiphysics FEM）**

以 MFEM、deal.II、Kratos、MOOSE 为代表。这是一个**框架**而非单一物理场工具。核心特征包括：高阶单元（p-自适应）、矩阵无关求解（matrix-free）、混合有限元（mixed FEM）、H(div) 和 H(curl) 空间。对比前两者，它的"通用性"是量级的差异——它要处理的不只是力学方程，还有 Maxwell 方程组、Navier-Stokes、Darcy 流、甚至耦合问题。

---

.. raw:: html

   <h2 class="section-title">二、型函数（Shape Functions）的深度对比</h2>

**这是全文最核心的问题。答案是：并不相同，但有交叉。**

.. raw:: html

   <h3>2.1 一张表看清全部差异</h3>

.. raw:: html

   <table>
  <tr>
    <th>型函数类型</th>
    <th>多项式次数</th>
    <th>连续性</th>
    <th>用在哪个分支</th>
    <th>典型单元</th>
  </tr>
  <tr>
    <td><strong>Lagrange 型</strong></td>
    <td>p=1,2,...（节点基）</td>
    <td>C⁰</td>
    <td><span class="tag-blue">三者通用</span></td>
    <td>杆单元、六面体、四面体、板壳</td>
  </tr>
  <tr>
    <td><strong>Hermite 型</strong></td>
    <td>p=3（三次）</td>
    <td>C¹</td>
    <td><span class="tag-green">仅结构有限元</span></td>
    <td>欧拉梁、克希霍夫板</td>
  </tr>
  <tr>
    <td><strong>阶谱型（Hierarchical）</strong></td>
    <td>p=1~15+（模态基）</td>
    <td>C⁰（谨慎可得 C¹）</td>
    <td><span class="tag-purple">多物理场有限元</span></td>
    <td>p-FEM 单元</td>
  </tr>
  <tr>
    <td><strong>Nédélec 型</strong></td>
    <td>p=0,1,...（棱边基）</td>
    <td>H(curl) 切向连续</td>
    <td><span class="tag-purple">仅多物理场</span></td>
    <td>电磁场棱边元</td>
  </tr>
  <tr>
    <td><strong>Raviart-Thomas 型</strong></td>
    <td>p=0,1,...（面基）</td>
    <td>H(div) 法向连续</td>
    <td><span class="tag-purple">仅多物理场</span></td>
    <td>Darcy 流，电位移</td>
  </tr>
  <tr>
    <td><strong>NURBS（等几何分析）</strong></td>
    <td>p=2,3,...</td>
    <td>C¹ 或更高</td>
    <td><span class="tag-orange">IGA 专用</span></td>
    <td>光滑壳、CAD互联</td>
  </tr>
</table>

.. raw:: html

   <h3>2.2 Lagrange 型函数——唯一三者共享的型函数</h3>

Lagrange 型函数是 **三者唯一共享** 的型函数族。但即使都是用 Lagrange 型函数，在不同分支中的角色也不同：

- **固体力学**：完全的三维 Lagrange 插值。每个节点一个基函数，插值位移场。
- **结构有限元（板壳）**：面内用 Lagrange 双线性/二次插值，但法向的**转动自由度**不是通过导数型函数实现的（不像 Hermite），而是通过引入额外的转动 DOF，在单元刚度矩阵中与平动耦合。
- **多物理场有限元**：Lagrange 仅用于 H¹ 空间（标量场和矢量场）。对于电磁场和流体，需要 Nédélec 和 RT 空间——这些在结构有限元和固体力学中根本不存在。

.. raw:: html

   <div class="highlight-box">
   <strong>关键理解：</strong><br><br>
   <p>Lagrange 型函数只是"搭积木"中最基础的那一块。三个分支都用了这一块，但在其之上架构了完全不同的上层结构。</p>
   </div>

.. raw:: html

   <h3>2.3 Hermite 型函数——结构有限元的独门武器</h3>

梁单元是理解结构有限元独特性的最佳窗口。

对于 Euler-Bernoulli 梁，控制方程是四阶微分方程：

.. raw:: html

   <div style="background: #f8f9fa; padding: 12px; border-radius: 4px; margin: 8px 0; text-align: center; font-family: serif; font-size: 1.1em;">
   EI · d⁴v/dx⁴ = q(x)
   </div>

这要求挠度场 **v(x)** 在单元交界处同时保持 **v** 和 **dv/dx** 的连续——即 **C¹ 连续**。Lagrange 型函数只能保证 C⁰ 连续，无法做到。

解决方案是 **Hermite 插值**：在每个节点上同时定义 **v**（挠度）和 **θ = dv/dx**（转角）两个自由度。三次 Hermite 型函数可以保证 C¹ 连续性。具体形式（归一化到 [0,1]）：

.. raw:: html

   <div style="background: #f8f9fa; padding: 8px 12px; border-radius: 4px; margin: 8px 0; font-family: monospace; font-size: 0.95em;">
   H₁(ξ) = 1 − 3ξ² + 2ξ³    （对应节点1的挠度）<br>
   H₂(ξ) = ξ(1 − ξ)²         （对应节点1的转角）<br>
   H₃(ξ) = 3ξ² − 2ξ³          （对应节点2的挠度）<br>
   H₄(ξ) = ξ²(ξ − 1)          （对应节点2的转角）
   </div>

**要点**：这里转角自由度是型函数自然定义的——型函数构造本身决定了"每个节点有两个自由度（v 和 θ）"。这与 Mindlin 板壳完全不同。

.. raw:: html

   <div class="warning-box">
   <strong>⚠️ 一个常见的误解</strong><br><br>
   <p>有人以为 Mindlin 板壳中的转动自由度和 Euler 梁中的转角是一回事。<strong>错。</strong></p>
   <ul>
     <li><strong>Euler 梁（Hermite）</strong>：转角 θ = dv/dx，是挠度的导数。C¹ 连续自然实现。</li>
     <li><strong>Mindlin 板壳</strong>：转角是独立的自由度 β（截面法线转角），与挠度无关。使用 C⁰ Lagrange 型函数即可。转角和平动通过剪切刚度耦合。</li>
   </ul>
   <p>后果：Mindlin 板壳在薄板极限下出现剪切锁死，需要减缩积分或假定自然应变法（ANS）。Euler 梁没有这个问题（但只适用于细长梁）。</p>
   </div>

.. raw:: html

   <h3>2.4 Hierarchical 型函数——多物理场独有的p-自适应能力</h3>

结构有限元和固体力学有限元通常使用**固定阶次**的型函数（线性或二次）。如果需要提高精度，靠的是**h-加密**（细化网格）。

多物理场有限元（MFEM、deal.II）支持 **p-自适应**——在不改变网格拓扑的前提下，提高单元内的多项式阶次。这需要**Hierarchical 型函数**：

- Lagendre 多项式作为基函数
- 低阶基是高阶基的"子集"——提高阶次时只需添加新的基函数，无需修改已有的
- 典型阶次可达 p=15 甚至更高

.. raw:: html

   <div class="info-box">
   <strong>p-自适应的工程意义</strong><br><br>
   <p>对于光滑解（如弹性力学中应力集中区域以外），p-自适应的收敛速度是指数级的，远快于 h-加密的代數收敛。这也是为什么多物理场框架在处理高精度问题时更有优势。</p>
   <p>但在处理强不连续（如裂纹扩展）时，p-自适应效果有限——此时需要 XFEM/GFEM 或重划分。</p>
   </div>

---

.. raw:: html

   <h2 class="section-title">三、单元构造的深度对比（等参 · 亚参 · 超参）</h2>

.. raw:: html

   <h3>3.1 三种映射方式</h3>

.. raw:: html

   <table>
  <tr>
    <th>映射方式</th>
    <th>几何插值阶次</th>
    <th>场插值阶次</th>
    <th>用在哪个分支</th>
    <th>典型例子</th>
  </tr>
  <tr>
    <td><strong>等参（Isoparametric）</strong></td>
    <td>= 场阶次</td>
    <td>= 几何阶次</td>
    <td><span class="tag-blue">固体力学、结构有限元</span></td>
    <td>8节点六面体、4节点四边形壳</td>
  </tr>
  <tr>
    <td><strong>亚参（Subparametric）</strong></td>
    <td>高于场阶次</td>
    <td>低于几何阶次</td>
    <td><span class="tag-purple">多物理场p-FEM、部分壳单元</span></td>
    <td>线性几何网格+高阶场插值</td>
  </tr>
  <tr>
    <td><strong>超参（Superparametric）</strong></td>
    <td>低于场阶次</td>
    <td>高于几何阶次</td>
    <td><span class="tag-green">结构有限元（梁）</span></td>
    <td>直梁+三次 Hermite 场插值</td>
  </tr>
</table>

.. raw:: html

   <h3>3.2 为什么结构有限元的梁单元是超参的？</h3>

这是很多人忽略但极其关键的一点。

一个典型的 2D Euler 梁单元：
- **几何**：两个节点定义一条直线段（线性 Lagrange 插值） → 几何阶次 p=1
- **场**：三次 Hermite 插值（含挠度和转角） → 场阶次 p=3

因此：**几何阶次(1) < 场阶次(3)** → 超参单元。

这带来的后果是：单元几何不能弯曲（只能由多个直梁分段近似曲线）。如果需要精确的曲线梁，必须引入第三个几何节点（二次几何）+ 三次 Hermite 场，但这已经不属于标准 2 节点梁单元了。

.. raw:: html

   <h3>3.3 为什么多物理场 p-FEM 是亚参的？</h3>

多物理场框架的典型做法：
- **几何**：用一个较粗的网格描述几何，每个单元使用二次甚至三次 Lagrange 映射处理边界弯曲
- **场**：单元内部使用 p=5~15 的高阶 Hierarchical 型函数

因此：**几何阶次(2~3) > 场阶次(用 p 表示，但几何的高阶已足够精确边界)**

这避免了"每个几何节点都携带高自由度场"的问题——几何节点的 DOF 数量是固定的（通常 3 个平动），而场自由度在单元内部独立管理。

.. raw:: html

   <div class="highlight-box">
   <strong>核心洞察：映射方式的不同本质上是"谁控制精度"的设计哲学不同</strong><br><br>
   <ul>
     <li><strong>等参（固体力学）</strong>：网格和精度绑定——加密网格同时提升几何和场精度。简单但效率有限。</li>
     <li><strong>超参（结构梁）</strong>：因为物理需要高次场（C¹），但几何不需要。两个节点定义直线就够了。</li>
     <li><strong>亚参（p-FEM）</strong>：几何用一个适中的网格描述，精度由场阶次控制。最优解耦策略。</li>
   </ul>
   </div>

---

.. raw:: html

   <h2 class="section-title">四、自由度与连续性要求</h2>

.. raw:: html

   <h3>4.1 每个节点的 DOF 数量</h3>

.. raw:: html

   <table>
  <tr>
    <th>分支</th>
    <th>单元类型</th>
    <th>每个节点的DOF</th>
    <th>DOF内容</th>
  </tr>
  <tr>
    <td rowspan="2"><strong>固体力学</strong></td>
    <td>六面体 / 四面体</td>
    <td><strong>3</strong></td>
    <td>uₓ, uᵧ, u_z</td>
  </tr>
  <tr>
    <td>平面应力/应变</td>
    <td><strong>2</strong></td>
    <td>uₓ, uᵧ</td>
  </tr>
  <tr>
    <td rowspan="2"><strong>结构有限元</strong></td>
    <td>Euler 梁</td>
    <td><strong>6</strong></td>
    <td>uₓ, uᵧ, u_z, θₓ, θᵧ, θ_z</td>
  </tr>
  <tr>
    <td>Mindlin 壳</td>
    <td><strong>5 或 6</strong></td>
    <td>3平移 + 2或3转动（drilling 可选）</td>
  </tr>
  <tr>
    <td rowspan="3"><strong>多物理场</strong></td>
    <td>标量场（温度/压力）</td>
    <td><strong>1</strong></td>
    <td>T 或 p</td>
  </tr>
  <tr>
    <td>矢量场 H(curl)</td>
    <td><strong>棱边DOF</strong></td>
    <td>非节点自由度，每个棱边1个</td>
  </tr>
  <tr>
    <td>混合格式</td>
    <td><strong>多种</strong></td>
    <td>如位移+压力 (uₓ, uᵧ, u_z, p)</td>
  </tr>
</table>

.. raw:: html

   <h3>4.2 连续性要求的层级</h3>

连续性要求是区分这三个分支最重要的数学维度。

.. raw:: html

   <div style="background: #f8f9fa; padding: 16px; border-radius: 8px; margin: 12px 0;">
   <strong>De Rham 复形（De Rham Complex）</strong><br><br>
   <p style="font-family: monospace; font-size: 0.95em;">
   H¹ —∇→ H(curl) —∇×→ H(div) —∇·→ L²<br><br>
   Lagrange    Nédélec    Raviart-Thomas    Discontinuous<br>
   (C⁰ 连续)  (切向连续)  (法向连续)       (无连续要求)
   </p>
   <p>数值方法只有填满这个复形才能保证无杂散模式的电磁场/流体求解。</p>
   <p><strong>结构有限元和固体力学有限元只覆盖了最左边的 H¹ 空间。</strong>多物理场有限元的本质特征就是覆盖了整个 De Rham 复形。</p>
   </div>

.. raw:: html

   <table>
  <tr>
    <th>连续性</th>
    <th>固体力学</th>
    <th>结构有限元</th>
    <th>多物理场</th>
  </tr>
  <tr>
    <td>H¹（C⁰ 连续）</td>
    <td><span class="tag-green">✅ 全部</span></td>
    <td><span class="tag-green">✅ 板壳+部分梁</span></td>
    <td><span class="tag-green">✅ 标量场</span></td>
  </tr>
  <tr>
    <td>C¹ 连续</td>
    <td><span class="tag-red">❌ 不需要</span></td>
    <td><span class="tag-green">✅ Euler梁/Kirchhoff板</span></td>
    <td><span class="tag-orange">⚠️ 特殊情况下</span></td>
  </tr>
  <tr>
    <td>H(curl) 切向连续</td>
    <td><span class="tag-red">❌ 不存在</span></td>
    <td><span class="tag-red">❌ 不存在</span></td>
    <td><span class="tag-green">✅ Nédélec棱边元</span></td>
  </tr>
  <tr>
    <td>H(div) 法向连续</td>
    <td><span class="tag-red">❌ 不存在</span></td>
    <td><span class="tag-red">❌ 不存在</span></td>
    <td><span class="tag-green">✅ Raviart-Thomas面元</span></td>
  </tr>
  <tr>
    <td>L² 无连续</td>
    <td><span class="tag-orange">⚠️ 不常用</span></td>
    <td><span class="tag-orange">⚠️ DG方法</span></td>
    <td><span class="tag-green">✅ DG-FEM</span></td>
  </tr>
</table>

.. raw:: html

   <div class="warning-box">
   <strong>核心结论：</strong><br><br>
   <p>固体力学有限元和结构有限元本质上只处理 H¹（C⁰）空间的问题。结构有限元中出现的 C¹ 连续（Euler 梁、Kirchhoff 板）是特殊情况的处理。</p>
   <p>多物理场有限元的独特之处在于它覆盖了整个 De Rham 复形——这不仅是"量"的增加，更是"质"的不同。如果你不需要 Maxwell 方程、Darcy 流、或混合格式，那么大可不必上多物理场框架。</p>
   </div>

---

.. raw:: html

   <h2 class="section-title">五、积分策略与数值问题</h2>

.. raw:: html

   <h3>5.1 高斯积分阶次的选择</h3>

.. raw:: html

   <table>
  <tr>
    <th>单元</th>
    <th>型函数</th>
    <th>推荐积分阶次</th>
    <th>特殊处理</th>
  </tr>
  <tr>
    <td>线性六面体（8节点）</td>
    <td>三线性 Lagrange</td>
    <td>2×2×2</td>
    <td>近不可压需B-bar或u-p格式</td>
  </tr>
  <tr>
    <td>二次六面体（20节点）</td>
    <td>Serendipity Lagrange</td>
    <td>3×3×3</td>
    <td>减缩积分2×2×2可导致沙漏</td>
  </tr>
  <tr>
    <td>三次 Hermite 梁</td>
    <td>Hermite (p=3)</td>
    <td>2点Gauss（1D）</td>
    <td>非线性问题可能需要3点</td>
  </tr>
  <tr>
    <td>线性四边形壳（Mindlin）</td>
    <td>双线性 Lagrange</td>
    <td>膜/弯2×2, 剪切1点</td>
    <td><strong>减缩积分</strong>解剪切锁死</td>
  </tr>
  <tr>
    <td>高阶六面体（p=5）</td>
    <td>Hierarchical</td>
    <td>6×6×6</td>
    <td>可用 Gauss-Lobatto 适配矩阵无关</td>
  </tr>
  <tr>
    <td>Nédélec 棱边元（最低阶）</td>
    <td>棱边基</td>
    <td>1点/方向</td>
    <td>需调校积分规则避免角点奇点</td>
  </tr>
</table>

.. raw:: html

   <h3>5.2 锁死问题——结构有限元独有的烦恼</h3>

如果只写过固体力学有限元代码，你可能对"锁死"只有一个抽象的认识。但在结构有限元中，锁死是每天都要面对的实际问题。

**剪切锁死（Shear Locking）**：
- Mindlin 板壳在厚度趋近于零时，剪切刚度被高估
- 原因是剪切应变采用与弯曲应变相同的积分阶次
- 解决：减缩积分（reduced integration）、假定自然应变法（ANS）、增强假定应变法（EAS）

**薄膜锁死（Membrane Locking）**：
- 曲壳单元在纯弯曲时，薄膜应变被高估
- 表现：壳的弯曲刚度偏大
- 解决：同样需要特殊处理

**体积锁死（Volumetric Locking）**：
- 固体力学和结构有限元都会遇到
- 当材料接近不可压缩（ν → 0.5）时，体积应变被锁定
- 解决：B-bar 方法、u-p 混合格式、减缩积分

.. raw:: html

   <div class="highlight-box">
   <strong>作者的开发经验：</strong><br><br>
   <p>如果你已经开发了固体力学有限元代码，再开发结构有限元（板壳单元）时，最大的变化不是型函数——而是需要处理各种锁死问题。80%的调试时间花在了"为什么这个单元在薄板极限下刚度偏大"上。</p>
   </div>

---

.. raw:: html

   <h2 class="section-title">六、结构有限元 vs 固体力学有限元：到底差在哪？</h2>

这是用户最直接的问题。答案是：

.. raw:: html

   <div class="success-box">
   <strong>从连续介质力学的角度看，二者本质上是一样的——都是求解平衡方程 + 本构关系 + 几何方程。</strong><br><br>
   <p><strong>区别在于"假设"的介入位置不同：</strong></p>
   </div>

.. raw:: html

   <h3>6.1 固体力学有限元：无假设</h3>

直接从三维连续介质力学出发：
- 平衡方程：σ_{ij,j} + f_i = 0
- 几何方程：ε_{ij} = (u_{i,j} + u_{j,i})/2
- 本构：σ = D:ε
- 在三维空间离散，每个节点 3 个平动 DOF

**优点**：精确，无额外假设。适用于任何几何和加载条件。
**代价**：计算量大。对于薄壁结构，需要非常细的网格才能捕捉弯曲行为。

.. raw:: html

   <h3>6.2 结构有限元：引入运动学假设</h3>

在从"三维问题"降为"一维/二维问题"时，引入了**运动学假设**：

- **Euler-Bernoulli 梁假设**：平截面假设 + 截面法线保持法线（忽略剪切变形）
- **Timoshenko 梁假设**：平截面假设 + 截面法线转为直线（考虑剪切变形）
- **Kirchhoff 板假设**：中面法线保持法线（忽略剪切 → 四阶方程 → C¹ 要求）
- **Mindlin 板假设**：中面法线转为直线（考虑剪切 → C⁰ 足够）

这些假设的后果：

.. raw:: html

   <table>
   <tr><th>假设</th><th>带来的好处</th><th>引入的问题</th></tr>
   <tr><td>维度缩减（3D→1D/2D）</td><td>DOF减少1~2个数量级</td><td>丢失了厚度方向的三维应力细节</td></tr>
   <tr><td>平截面假设</td><td>弯曲问题闭式解</td><td>不能处理截面翘曲</td></tr>
   <tr><td>忽略剪切（Kirchhoff）</td><td>DOF更少</td><td>不适用于厚梁/板</td></tr>
   <tr><td>保留剪切（Mindlin）</td><td>适用于中厚板</td><td>剪切锁死</td></tr>
   </table>

.. raw:: html

   <div class="info-box">
   <strong>一句话总结：</strong><br><br>
   <p>结构有限元 = 固体力学有限元 + 运动学假设 + 维度缩减 + 转动自由度 + 锁死处理。</p>
   <p>从型函数的角度看，结构有限元比固体力学多了一类（Hermite），但少了许多（不需要处理 H(curl) 和 H(div) 空间）。</p>
   <p>从单元构造的角度看，结构有限元引入了超参映射（梁），而固体力学几乎只用等参。</p>
   </div>

---

.. raw:: html

   <h2 class="section-title">七、多物理场有限元的独特之处</h2>

.. raw:: html

   <h3>7.1 不只是"多个物理场"</h3>

多物理场有限元框架（MFEM、deal.II、Kratos）和结构有限元软件（Abaqus）的区别**不是"一个是多物理场、一个只做力学"这么简单**。

本质区别在于**框架设计的抽象层级**：

.. raw:: html

   <table>
   <tr><th>维度</th><th>Abaqus / CalculiX</th><th>MFEM / deal.II / Kratos</th></tr>
   <tr><td>单元库</td><td>封闭的、预定义的</td><td>开放的、用户可定义任意型函数空间</td></tr>
   <tr><td>型函数</td><td>Lagrange + Hermite</td><td>Lagrange + Hierarchical + Nédélec + RT + ...</td></tr>
   <tr><td>多项式阶次</td><td>固定（1或2次）</td><td>可变（1~15+次）</td></tr>
   <tr><td>线性代数</td><td>全矩阵组装</td><td>矩阵无关（matrix-free），支持 GPU</td></tr>
   <tr><td>网格</td><td>纯 h-自适应</td><td>h- + p- + hp-自适应</td></tr>
   <tr><td>混合格式</td><td>有限支持（u-p ONLY）</td><td>一等公民</td></tr>
   <tr><td>非连续Galerkin</td><td>不支持</td><td>内置支持</td></tr>
   </table>

.. raw:: html

   <h3>7.2 矩阵无关（Matrix-Free）方法的革命性</h3>

传统有限元（Abaqus、CalculiX 和你写的代码）：
1. 遍历所有单元，组装全局刚度矩阵 K
2. 求解 K·u = f
3. K 的规模是 DOF × DOF，对于 100 万 DOF 的问题，仅存储 K 就需要 8GB（如果稀疏）

矩阵无关方法（MFEM、deal.II）：
1. 不组装 K
2. 每次需要 K·v 时，逐单元计算并累加（通过 sum-factorization 技术）
3. 存储需求仅为 O(DOF) 而非 O(DOF²)
4. 特别适合高阶单元（p≥3），因为高阶单元的全矩阵 K 填充率更高，不组装的收益更大

.. raw:: html

   <div class="warning-box">
   <strong>实践建议</strong><br><br>
   <p>如果你只做线弹性力学问题（DOF < 10⁶），传统全矩阵组装+直接稀疏求解器 的效率完全够用。不需要上矩阵无关。</p>
   <p>但如果你做大规模（DOF > 10⁷）或高阶（p ≥ 3）问题，矩阵无关方法可以节省一个数量级的内存和计算时间。</p>
   </div>

---

.. raw:: html

   <h2 class="section-title">八、一张表看懂全部区别</h2>

.. raw:: html

   <table>
   <tr>
     <th width="25%">对比维度</th>
     <th width="25%">固体力学有限元</th>
     <th width="25%">结构有限元</th>
     <th width="25%">多物理场有限元</th>
   </tr>
   <tr>
     <td><strong>代表软件</strong></td>
     <td>自研代码</td>
     <td>Abaqus / CalculiX / ANSYS</td>
     <td>MFEM / deal.II / Kratos / MOOSE</td>
   </tr>
   <tr>
     <td><strong>控制方程</strong></td>
     <td>三维连续介质力学</td>
     <td>维度缩减后（板/壳/梁理论）</td>
     <td>任意PDE（力学/电磁/流体/耦合）</td>
   </tr>
   <tr>
     <td><strong>型函数类型</strong></td>
     <td>Lagrange（C⁰）</td>
     <td>Lagrange（板壳）+ Hermite（梁）</td>
     <td>Lagrange + Hierarchical + Nédélec + RT</td>
   </tr>
   <tr>
     <td><strong>单元构造</strong></td>
     <td>等参</td>
     <td>等参（板壳）+ 超参（梁）</td>
     <td>等参/亚参 + Piola 映射</td>
   </tr>
   <tr>
     <td><strong>每个节点 DOF</strong></td>
     <td>3（纯平动）</td>
     <td>3~6（平动+转动）</td>
     <td>1~多种（视场而定）</td>
   </tr>
   <tr>
     <td><strong>连续性要求</strong></td>
     <td>C⁰（H¹）</td>
     <td>C⁰ + C¹（梁）</td>
     <td>H¹ + H(curl) + H(div) + L² 全部</td>
   </tr>
   <tr>
     <td><strong>锁死问题</strong></td>
     <td>体积锁死</td>
     <td>剪切锁死 + 薄膜锁死 + 体积锁死</td>
     <td>取决于具体格式</td>
   </tr>
   <tr>
     <td><strong>多项式阶次</strong></td>
     <td>1~2（固定）</td>
     <td>1~2（固定）</td>
     <td>1~15+（可变）</td>
   </tr>
   <tr>
     <td><strong>自适应</strong></td>
     <td>h-加密</td>
     <td>h-加密</td>
     <td>h- + p- + hp-</td>
   </tr>
   <tr>
     <td><strong>矩阵组装</strong></td>
     <td>全矩阵</td>
     <td>全矩阵</td>
     <td>全矩阵 / 矩阵无关</td>
   </tr>
   <tr>
     <td><strong>开发门槛</strong></td>
     <td><span class="tag-green">低</span></td>
     <td><span class="tag-orange">中</span></td>
     <td><span class="tag-red">高</span></td>
   </tr>
   <tr>
     <td><strong>适合的问题规模</strong></td>
     <td>中小规模</td>
     <td>各种规模</td>
     <td>大规模 + 高精度 + 多物理场</td>
   </tr>
   </table>

---

.. raw:: html

   <h2 class="section-title">九、总结：给"写过代码的工程师"的结论</h2>

.. raw:: html

   <div style="background: linear-gradient(135deg, #1a1a2e, #16213e); color: white; padding: 28px; border-radius: 12px; margin: 24px 0; font-size: 1.05em; line-height: 2;">
   <p style="font-size: 1.15em;"><strong>🎯 核心结论</strong></p>

   <ol style="margin-top: 8px;">
     <li><strong>型函数并不相同，但有交叉。</strong> Lagrange 是三者唯一的交集。结构有限元独有 Hermite。多物理场覆盖整个 De Rham 复形（Lagrange + Nédélec + Raviart-Thomas + Hierarchical）。</li>
     <br>
     <li><strong>单元构造不同。</strong> 固体力学几乎只用等参。结构有限元引入了超参（梁单元）。多物理场大量使用亚参（p-FEM）和 Piola 映射（H(div)/H(curl) 空间）。</li>
     <br>
     <li><strong>结构有限元 = 固体力学 + 运动学假设 + 维度缩减 + 转动DOF + 锁死处理。</strong> 如果你已经会写固体力学 FEM 代码，结构有限元的核心新知识是：维度缩减理论（板壳/梁假设）、转动 DOF 的处理、和锁死的识别与修复。</li>
     <br>
     <li><strong>多物理场不是"更高级的FEM"，而是"不同抽象的框架"。</strong> 它不是为了替代 Abaqus 而存在的——它服务于高精度、大规模、多物理场耦合的场景。如果只做力学，Abaqus/CalculiX 可能更高效。</li>
     <br>
     <li><strong>三个分支的底层数学是统一的（Galerkin 加权残值法），但上层的工程假设、型函数空间、和实现策略完全不同。</strong> 理解这些区别，比"学会另一个软件"有价值得多。</li>
   </ol>
   </div>

---

**参考文献**

1. Bathe, K.-J. *Finite Element Procedures*, 2nd ed., 2014.
2. Zienkiewicz, O.C., Taylor, R.L., Zhu, J.Z. *The Finite Element Method: Its Basis and Fundamentals*, 7th ed., 2013.
3. Hughes, T.J.R. *The Finite Element Method*, 1987/2000.
4. Monk, P. *Finite Element Methods for Maxwell's Equations*, 2003.
5. Brenner, S.C., Scott, L.R. *The Mathematical Theory of Finite Element Methods*, 3rd ed., 2008.
6. Bangerth, W., Rannacher, R. *Adaptive Finite Element Methods for Differential Equations*, 2003.
7. 王勖成《有限单元法》, 2003.
8. MFEM 官方文档: https://mfem.org
9. Kratos Multiphysics 文档: https://kratosmultiphysics.github.io
10. CalculiX 官方文档: http://www.calculix.de

*本文生成日期：2026年5月30日*