跳转至

第4章：有限元方法（The Finite Element Method）

1. 作者

T. Christian Gasser (KTH Royal Institute of Technology) 独立撰写本章, Ch 4 是Ch 3 连续介质力学的数值实现。Springer / KTH 流派 FEM 教材, 与 [30, 34, 53, 457, 617-619] 等经典教材对齐。

章节范围: 4.1 Introduction → 4.14 Summary (PDF p188-p252, 65 页)

核心定位: 提供有限元方法 (FEM) 在血管生物力学中的标准工作流: 空间离散化 → 弱形式推导 → 离散方程 → 约束处理 → 装配 → 稳定化 → 求解 → 案例研究。是 Ch 5-7 中所有具体仿真问题的数值工具基础。

2. 内容概述

本章建立 FEM 的完整工作流: - 4.1-4.2 介绍 PDE 三类 (双曲/椭圆/抛物), 给出空间离散化与形函数设计 - 4.3 弱形式推导 (Calculus of Variations), 处理扩散、对流-扩散、线性/非线性固体、不可压流体 - 4.4 离散化后的代数方程 (Galerkin), 给出单元刚度矩阵的标准推导 - 4.5 约束问题: Penalty / Lagrange / Augmented-Lagrange (与 Ch 3 不可压材料呼应) - 4.6 全局装配: 单元贡献整合到全局系统 - 4.7 数值稳定化: hourglass 模式、AD 问题的 Petrov-Galerkin / SUPG / GLS - 4.8 求解器: 稀疏线性系统、显式/隐式时间积分、Newton-Raphson、load incrementation、arc-length - 4.9-4.13 五个血管生物力学案例 (平面双轴、轴对称血管、球形膨胀、血管网络、3D 圆柱管)

Gasser 的核心立场: FEM 是"求解器", 物理问题才是核心; 工程读者必须理解算法在背景中的运行方式, 才能有效使用商业/开源软件 (Abaqus, FEniCS, COMSOL, ANSYS, FEBio) 并避免常见陷阱 (体积锁定, mesh dependency, 不可压病态)。

3. 核心方程与概念

3.1 空间离散化 (Sect. 4.2)

形函数 \(N_i(\xi)\): 在自然坐标 \(\xi \in [-1, 1]\) (1D) / \([-1, 1]^2\) (2D) / \([-1, 1]^3\) (3D) 上定义, 满足单位性 (\(\sum N_i = 1\)), Delta 性 (\(N_i(\xi_j) = \delta_{ij}\)), 线性独立性。

1D 形函数 (Eq. 4.3, 4.4): - 两节点线性: \(N_1 = (1-\xi)/2\), \(N_2 = (1+\xi)/2\) - 三节点二次: \(N_1 = (\xi^2 - \xi)/2\), \(N_2 = (\xi^2 + \xi)/2\), \(N_3 = 1 - \xi^2\)

2D 形函数 (Eq. 4.5): - 四边形双线性: \(N_i = \frac{1}{4}(1 \pm \xi_1)(1 \pm \xi_2)\) - 三角形 (面积坐标): \(N_i = A_i/A\)

3D 形函数: 六面体三线性, 四面体线性, NURBS, VEM (Virtual Element Method, 不显式引入形函数)。

梯度插值矩阵 \(\mathbf{B}\) (Eq. 4.6-4.7): 形函数梯度 \(\nabla N_i\), 通过雅可比变换 \(J = \partial x / \partial \xi\) 映射到物理坐标。应变 \(\boldsymbol{\varepsilon} = \mathbf{B}\mathbf{h}\)。

混合与杂交单元 (Sect. 4.2.3): 不同物理量用不同插值阶 (例如不可压 Navier-Stokes: 速度 P2, 压力 P1)。

3.2 PDE 分类 (Sect. 4.1)

  • 双曲 (Hyperbolic): 波传播; 不连续初值传播; 给出激波
  • 椭圆 (Elliptic): 弹性体变形; 解光滑; 边界数据影响整个域; 锐角处奇异性
  • 抛物 (Parabolic): 热传导; 双曲与椭圆之间

3.3 弱形式推导 (Sect. 4.3)

核心工具: Calculus of Variations — 乘以试函数 \(\delta u\) (在 Dirichlet 边界上为零), 积分, 分部积分, 应用散度定理。

3.3.1 扩散问题 (Eq. 4.8-4.9)

强形式: \(\text{div}(\nabla c) + \alpha = 0\), \(c = \bar{c}\) on \(\partial\Omega_c\), \(\nabla c \cdot \mathbf{n} = \bar{q}\) on \(\partial\Omega_q\)

弱形式:

\[ \int_\Omega \nabla\delta c \cdot \nabla c\, dv = \int_{\partial\Omega_q} \delta c \bar{q}\, ds + \int_\Omega \delta c \alpha\, dv \]

3.3.2 对流-扩散 (Eq. 4.10-4.11)

弱形式: \(\int_\Omega \delta c (\mathbf{v} \cdot \nabla c)\, dv + \int_\Omega \nu \nabla\delta c \cdot \nabla c\, dv = \int_{\partial\Omega_q} \delta c \bar{q}\, ds + \int_\Omega \delta c \alpha\, dv\)

3.3.3 线性固体 (Eq. 4.15-4.17)

弱形式 (= PVW):

\[ \int_{\partial\Omega_t} \mathbf{t} \cdot \delta\mathbf{u}\, ds + \int_\Omega \delta\mathbf{u} \cdot (\mathbf{b}_f - \rho\ddot{\mathbf{u}} - \eta\dot{\mathbf{u}})\, dv = \int_\Omega \boldsymbol{\sigma} : \text{grad}_s \delta\mathbf{u}\, dv \]

3.3.4 非线性固体 (Eq. 4.18-4.21)

内部虚功: \(\delta\mathcal{W}_{\text{int}} = \int_{\Omega_0} \delta\mathbf{E} : \mathbf{S}\, dV\)

方向导数 (Eq. 4.19) = 几何贡献 (Eq. 4.20) + 材料贡献 (Eq. 4.21)

\[ D_u \delta\mathcal{W}_{\text{int}} = \int_\Omega [\text{grad}_s \delta\mathbf{u} : \text{grad}_s \Delta\mathbf{u}\, \boldsymbol{\sigma} + \text{grad}_s \delta\mathbf{u} : \mathbf{C} : \text{grad}_s \Delta\mathbf{u}]\, dv \]

3.3.5 不可压流体 (Eq. 4.24-4.25)

两个独立场 \(\mathbf{v}, p\): - 连续性弱形式: \(\int_\Omega \delta p\, \text{div}\mathbf{v}\, dv = 0\) - 动量弱形式: \(\int_\Omega \delta\mathbf{v} \cdot [\rho(\partial_t \mathbf{v} + \mathbf{v} \cdot \nabla\mathbf{v})]\, dv + \int_\Omega \delta\mathbf{d} : \boldsymbol{\sigma}\, dv - \int_\Omega \text{div}\delta\mathbf{v}\, p\, dv = \int_{\partial\Omega_t} \delta\mathbf{v} \cdot \mathbf{t}\, ds + \int_\Omega \delta\mathbf{v} \cdot \mathbf{b}_f\, dv\)

3.4 离散有限元方程 (Sect. 4.4)

Galerkin 方法: 试函数与解函数用相同形函数 (vs Petrov-Galerkin 不同)。

扩散离散 (Eq. 4.26-4.27): \(\mathbf{D}\mathbf{h} - \mathbf{f} = 0\), 其中 \(D_{ij} = \int_{\Omega^e} B_{ai} B_{aj}\, dv\)

对流-扩散离散 (Eq. 4.28-4.29): \((\mathbf{K} + \mathbf{D})\mathbf{h} - \mathbf{f} = 0\), \(\mathbf{K}\) 是非对称

线性固体离散 (Eq. 4.30-4.31): \(\mathbf{M}\ddot{\mathbf{h}} + \mathbf{C}\dot{\mathbf{h}} + \mathbf{K}\mathbf{h} - \mathbf{f} = 0\)

3.4.1 杆单元 (Eq. 4.32-4.35)

质量矩阵 (Eq. 4.34): \(\mathbf{M} = \frac{m}{6}\begin{bmatrix} 2 & 1 \\ 1 & 2 \end{bmatrix}\) (一致质量); 集中质量 \(\mathbf{M}_{\text{lumped}} = \frac{m}{2}\mathbf{I}\)

刚度矩阵 (Eq. 4.35): \(\mathbf{K} = \frac{EA}{l}\begin{bmatrix} 1 & -1 \\ -1 & 1 \end{bmatrix}\)

3.4.2 非线性固体 (Eq. 4.36-4.38)

节点力: \(f_i = \int_{\Omega^e} B_{ai} \sigma_a\, dv\) (Eq. 4.36)

材料刚度 (Eq. 4.37): \(K^{\text{mat}}_{ij} = \int_{\Omega^e} B_{ai} C_{ab} B_{bj}\, dv\)

几何刚度 (Eq. 4.38): \(K^{\text{geo}}_{ij} = \int_{\Omega^e} \frac{\partial N_{ai}}{\partial x_c} \sigma_{ab} \frac{\partial N_{cj}}{\partial x_b}\, dv\)

压力 BC 刚度 (Eq. 4.39-4.42): follower load, \(K^p_{ij}\) 不对称。

3.4.3 不可压流 (Eq. 4.43-4.45)

\[ \mathbf{M}\dot{\mathbf{h}} + \mathbf{A}\mathbf{h}^2 + \mathbf{K}\mathbf{h} - \mathbf{D}^T\mathbf{q} - \mathbf{f} = 0 \]
约束: \(\mathbf{D}\mathbf{h} = 0\)

3.5 数值积分 (Sect. 4.4.6, Eq. 4.46-4.48)

Gauss-Legendre 积分: \(\int F(x)\, dx \approx \sum_l F(\xi_l) w_l\), \(n\) 节点精确积分 \(2n-1\) 次多项式。

Hourglass 项: 1 节点 (单点) 积分下, 4 节点四边形有 1 个零能模式 (hourglass mode, 形函数展开的 \(\xi_1 \xi_2\) 项)。需要正则化 (Ch 4.7)。

3.6 约束问题 (Sect. 4.5)

3.6.1 Penalty (Eq. 4.49, Sect. 4.5.1)

\(\mathbf{C} = \bar{\mathbf{C}} + \mathbf{C}^{\text{vol}}\), 令 \(C^{\text{vol}} \gg C\) 近似不可压。 - 优点: 简单, 不增加自由度 - 缺点: 病态 (大 \(C^{\text{vol}}\) → 大条件数)

3.6.2 Lagrange (Eq. 4.50-4.52, Sect. 4.5.2)

势能: \(\Pi_L(\mathbf{u}, p) = \Pi(\mathbf{u}) + \int_\Omega p(J(\mathbf{u}) - 1)\, dv\)

两个独立变量 \((\mathbf{u}, p)\) 的双场变分; \(p\) 是Lagrange 乘子 (= 静水压)。

优点: 约束严格满足 缺点: 刚度矩阵有零对角元 (pivot 零), 不能用快速 Gaussian 消元

3.6.3 Augmented-Lagrange (Eq. 4.53-4.54, Sect. 4.5.3)

\(\Pi_{AL}(\mathbf{u}, p) = \Pi_L(\mathbf{u}, p) - \frac{1}{2\kappa}\int p^2\, dv\)

\(\kappa \to \infty\) 时趋近 Lagrange; 但避免病态与零 pivot; Gasser 引用 [53] (general displacement implementation)。

Hu-Washizu 三场变分 (Example 4.8): \((\mathbf{u}, p, \theta)\), \(\theta\) 是单元体积变化; 用于消除体积锁定 (volume locking)。

3.7 装配 (Sect. 4.6)

全局刚度矩阵: \(\mathbf{K} = \sum_{e=1}^{n_e} \mathbf{K}_e\), 装配到全局节点。

带宽最小化: 自动节点重编号; 稀疏存储: skyline, CSR, COO。 边界条件施加: - Dirichlet: 修改全局刚度 (删除行/列或惩罚) - Neumann: 修改右端向量 (例 Example 4.11)

3.8 稳定化 (Sect. 4.7)

3.8.1 刚度正定性 (Sect. 4.7.1)

单元刚度 \(\mathbf{K}\) 必须正定, 否则小扰动 \(\Delta\mathbf{f}\) 引起大 \(\Delta\mathbf{h}\)。Lyapunov 稳定要求 \(\Delta\mathbf{h} = \mathbf{K}^{-1}\Delta\mathbf{f}\) 有界。

Hourglass 模式: 单点积分下零能模式; 需要阻尼或附加项。

3.8.2 AD 稳定化 (Sect. 4.7.2)

Péclet 数 \(Pe = vh/(2\nu)\) 大时, Galerkin 出现节点间振荡。

Full Upwind (Eq. 4.61): 加人工扩散 \(\nu^* = vh/2\); 简单但过度扩散。

Petrov-Galerkin (Eq. 4.62): 试函数添加 bubble 模式; 1D 形式 \(S_i = (1\pm\xi)/2 \mp 3\beta(1-\xi^2)/4\)

SUPG (Streamline Upwind Petrov-Galerkin, Eq. 4.64-4.65):

\[ P(\delta c) = \mathbf{v} \cdot \nabla \delta c, \quad \tau = \nu^* / |\mathbf{v}|^2 \]

GLS (Galerkin Least-Squares): 类似 SUPG 但加 \(\tau R(c)\) 项。

各向异性扩散 (Eq. 4.63): \(\nu^* = \nu^* \mathbf{v} \otimes \mathbf{v} / |\mathbf{v}|^2\), 沿流方向加扩散, 垂直于流方向不加。

3.9 求解器 (Sect. 4.8)

3.9.1 稀疏线性系统 (Sect. 4.8.1)

  • 直接法 (LU, LDU): 适合 ~100k 自由度; 包含 pivoting (partial, complete, rook)
  • 迭代法: 大问题 (>100k); 需 preconditioner; 求解器: GMRES, BiCGSTAB, multigrid

3.9.2 时间积分 (Sect. 4.8.2)

Forward/Backward Euler: 一阶精度; 全局误差 \(O(\Delta t)\), 局部误差 \(O(\Delta t^2)\) - Forward Euler: 显式, 条件稳定 (CFL) - Backward Euler: 隐式, 无条件稳定

Runge-Kutta: 更高阶 Newmark 方法: \(\mathbf{K}^* = k_1 \mathbf{K} + k_2 \mathbf{D} + k_3 \mathbf{M} - \mathbf{K}_f\)

3.9.3 Newton-Raphson (Sect. 4.8.6.1)

\(\mathbf{h} \leftarrow \mathbf{h} - \delta\mathbf{h}\), \(\delta\mathbf{h}\) from \(\mathbf{K}^*(\mathbf{h})\delta\mathbf{h} = 0\)。二次收敛 (Eq. 4.83): \(\log|r_{n+1}| / \log|r_n| \approx 2\)。

3.9.4 加载增量 (Sect. 4.8.6.2)

把总载荷分为 \(n_l\) 步施加, \(\Delta\lambda = 1/n_l\)。两层嵌套 loop: 外层载荷增量, 内层 Newton-Raphson。

3.9.5 Arc-Length (Sect. 4.8.6.4, Eq. 4.89-4.90)

引入 \(\lambda\) 作为额外自由度, 约束 \(\|\Delta\mathbf{u}\|^2 + \eta^2 \Delta\lambda^2 \|\mathbf{f}\|^2 = \Delta h^2\)。

  • \(\eta = 0\): 位移控制
  • \(\eta = \infty\): 力控制
  • 通用 \(\eta\): 通过极限点 (snap-through) 和回跳 (snap-back)

伪位移控制 (Example 4.13): 用单个节点的位移作为控制量, 比完整 arc-length 简单。

3.9.6 显式 (Sect. 4.8.5)

\(\mathbf{M}^n \Delta\ddot{\mathbf{h}}^{n+1} = \mathbf{K}_f^n \Delta\mathbf{h}^n - \mathbf{D}^n \Delta\dot{\mathbf{h}}^n - \mathbf{K}^n \Delta\mathbf{h}^n\) (Eq. 4.77)

用集中质量 (对角), 不需解方程组。 CFL 条件 (Eq. 4.78): \(\Delta t < \alpha \min(h/c)\)。

3.9.7 不可压流的投影法 (Sect. 4.8.6.5, Eq. 4.91-4.95)

Chorin 投影法 (Eq. 4.92-4.93): 1. 中间速度 \(\mathbf{h}^*\): \(M\dot{\mathbf{h}}^* + A\mathbf{h}^{*2} + K\mathbf{h}^* - \mathbf{f} = 0\) 2. 压力: \(Dq - f_q = 0\) (解 Poisson) 3. 速度修正: \(M(\mathbf{h} - \mathbf{h}^*)/\Delta t - \mathbf{D}^T \mathbf{q} = 0\)

3.10 案例研究 (Sect. 4.9-4.13)

案例 模型 网格 关键结论
4.9 平面双轴 4 节点四边形 (5×5), 二次三角 69k dof, 不可压 Yeoh FEM σ11=131.85 kPa vs 解析 138.0 kPa (中心区)
4.10 圆柱血管膨胀 轴对称二次三角 2.5k dof, 不可压 Yeoh 壁厚弯曲效应使 FEM 比膜模型更硬
4.11 球形膨胀 轴对称 506 个二次三角 neoHookean, G=422.5 kPa 极限点 → 用伪位移控制
4.12 动脉网络 锥形管 + Stokes 流 468 dof, 主动脉 q=30 ml/s, 出口 p=13 kPa 给压力/流率/速度分布
4.13 圆柱管 3D 流 四面体 250k dof, 不可压 Newtonian Re=2880, 入口 v=40 cm/s 入口长度 lh=345.6 cm, Poiseuille 未建立

4. 关键结论

  • [K1] 弱形式是 FEM 的基础: 不是直接从强形式开始, 而是用 Calculus of Variations 推导
  • [K2] Galerkin = Petrov-Galerkin (物理场): Galerkin 用同形函数, Petrov-Galerkin 不同; Petrov-Galerkin 用于稳定化
  • [K3] 材料 vs 几何刚度: 大变形下两者均需计算; 几何刚度在压力载荷 (follower load) 时尤其重要
  • [K4] 不可压三方法: Penalty 简单但病态; Lagrange 严格但 pivot 零; Augmented-Lagrange 折中
  • [K5] 体积锁定的根源: 等参单元对纯体积变形不敏感 (单点积分 hourglass); 解决: 减积分 + 稳定化, 或用 Q1P0 (杂交) 或 Hu-Washizu 三场
  • [K6] CFL 条件限制显式时间步: \(\Delta t < \alpha h/c\); 隐式无此限制但需解方程组
  • [K7] Newton-Raphson 二次收敛: 3-6 步即收敛, 但需良好初值; 加载增量是标准补救
  • [K8] Arc-length 通过极限点: 位移/力控制都不能过极限点, arc-length 引入 \(\lambda\) 作 DOF
  • [K9] Péclet 数控制 AD 稳定化: \(Pe > 1\) 需 SUPG/GLS; \(\beta\) 太小无效果, 太大过度扩散
  • [K10] 三维血管仿真 = 不可压 + 大变形 + 接触 + 流体-结构耦合: Ch 4 的 5 个案例是基础范式
  • [K11] 模态分析 = 简化的稳态动力学: \(K\phi = \omega^2 M\phi\); 实际血管生物力学中常用于预应力下的模态
  • [K12] 稀疏求解器的选择: 100k dof 以下用直接法 (LU, LDU), 以上用迭代法 (GMRES + AMG preconditioner); 大规模 3D 仿真必须用并行求解器
  • [K13] 多物理耦合 (Ch 4.13.3): 三元 WindKessel 出口 BC + 不可压 Navier-Stokes 入口 = FSI 雏形; 真正的 FSI 留到 Ch 6 之后

5. 挑战和开放性问题

  • [C1] 不可压处理仍是工程难点: Penalty 病态 + Lagrange 零 pivot + Augmented-Lagrange 经验 \(\kappa\) 选择; 实际工程代码 (FEniCS, Abaqus) 各有不同实现, 标准难定
  • [C2] 稳定化的"过度正则化": SUPG/GLS 在低 Péclet 数时反而降低精度; 选择 \(\beta\) 没有解析最优; 自适应方法存在但不成熟
  • [C3] Hourglass 模式的人工阻尼: 物理上没有对应, 仅是数值技巧; 不同阻尼方案给出不同结果
  • [C4] 大变形下的网格质量: Ch 4.11 (球形膨胀) 提到严重网格畸变; ALE (Arbitrary Lagrangian-Eulerian) 或 r-adaptivity 是标准补救, 但实现复杂
  • [C5] Newton-Raphson 的初值敏感性: 极限点附近不收敛; 弧长法是补救, 但需要额外的代码复杂性
  • [C6] FSI (流体-结构耦合) 的实现: Ch 4.13.3 用 WK 边界是简化; 真正的 ALE FSI 在 Ch 4 没有涵盖, 但 Gasser 引用 [546] 留给 Ch 6
  • [C7] 计算效率: 3D 圆柱管 250k dof 已"高分辨率", 但真实血管 (e.g., 整个冠脉循环) 需要 million+ dof; GPU 加速, model order reduction, ML 替代 FEM 都在发展中
  • [C8] 机器学习替代本构: 4.9 案例用 Yeoh (两参数), 但血管组织需要 8+ 参数的 Holzapfel-Gasser-Ogden 模型; ML 神经网络本构 [256] 已成熟, Gasser 应当提及
  • [C9] 数据驱动 FEM: 4.13.3 假设 WindKessel 参数已知; 实际临床中, 这些参数需要从测量数据反演; 贝叶斯 FEM (Bayesian FEM) 正在发展
  • [C10] 2022 年的现状: NVIDIA Modulus, JuliaFEM, FEniCSx 已广泛使用, 但 Gasser 教材仅提到 FEniCS 一次; 现代工程读者可能用不上 Abaqus
  • [C11] 神经网络-PDE 求解器: PINN (Physics-Informed Neural Networks) 在血管生物力学已开始应用; Gasser 应当明确对比传统 FEM 与 PINN 的精度-效率权衡

6. 个人反思与批判性分析

6.1 与本人血管生物力学研究的直接对接

本章是本人 FEniCS 仿真代码的理论基础:

  • 弱形式推导 (Eq. 4.17, 4.18): 本人 G&R 代码中, PVW 是 \(\int \mathbf{S}:\delta\mathbf{E}\, dV\), 直接对应
  • Newton-Raphson 求解 (Eq. 4.82): 本人代码中标准的非平衡迭代
  • Augmented-Lagrange 不可压处理 (Eq. 4.53): 本人用 FEniCS 的 mixed function space Taylor-Hood (P2/P1), 实际就是 Augmented-Lagrange
  • Arc-length 法 (Eq. 4.89): 本人 G&R 模拟中经常遇到极限点 (e.g., 颈动脉狭窄后的负刚度), arc-length 是关键
  • Case 4.10 (圆柱血管膨胀): 直接对应本人 FEniCS 中的 AAA 仿真模板; Yeoh 本构 \(c_1, c_2\) 与 Gasser 引用 [435] 的标定方法一致
  • Case 4.13 (3D 圆柱管流): 对应本人研究中的流固耦合 (FSI) 简化情形

6.2 Gasser 框架的局限性

  • 缺少 ALE FSI: Ch 4 是"标准 FEM", 但流体域变形 (vessel wall 移动导致 lumen 变化) 没有涵盖; Gasser 引用 [546] 留给 Ch 6
  • 接触/摩擦完全缺失: 血管-周围组织, 支架-血管壁, 斑块-管腔, 都需要接触算法; Ch 4 没讲
  • 多尺度均质化: RVE (Ch 3 引入) 在 FEM 中如何实现? 实际血管组织有微观结构 (胶原纤维, SMC), 但 Ch 4 仅给出 Yeoh (唯象) 案例
  • 本构标定: 4.9 案例用 \(c_1, c_2\) 直接给定, 实际工程需要从双轴实验数据反演; Gasser 引用 [482] 但反演算法缺位
  • 网格质量与收敛性: 实际工程中, "网格是否足够细" 是关键问题; Gasser 提到 "1 节点积分导致 hourglass", 但自适应性 (h/p-adaptivity) 没有
  • 并行化: Ch 4 完全没有提到 MPI/OpenMP/GPU; 现代血管仿真 3D 案例 (250k dof) 在笔记本上需要数小时, 大规模仿真需要并行求解器 (PETSc, Trilinos, Ginkgo)
  • 机器学习替代: 4.14 Summary 提到 "machine learning", 但具体如何用 ML 改进 FEM 没有; 现代综述 (e.g., [256]) 已成熟
  • 错误估计与验证: 4.14 Summary 提到"先简单模型, 再细化", 但量化误差估计 (a posteriori error estimation) 缺位; ASME V&V 40 标准在 Ch 1 提到但 Ch 4 没展开

6.3 一个观察: "FEM 的边界" — Ch 4 结束的位置

Ch 4 在 4.14 Summary 给出"建模哲学": "从简单开始, 逐步细化"。这与 Ch 1 IMA 原则呼应: 模型复杂度不是越高越好, 而是与 IMA 匹配。Ch 5-7 将给出具体应用, Ch 4 提供工具。

但 Ch 4 的最后一个 case study (4.13.3) 实际是fluid-structure interaction 的雏形 (WK 出口 BC + 不可压 Navier-Stokes 入口)。这暗示 Ch 4 之后, Gasser 应当展开FSI——但 Ch 5-7 都是单物理 (结构或流体, 不同时)。这是 Gasser 教材最大盲点: 真正的血管生物力学是 FSI + G&R + 多物理耦合, 但教材没有完整 FSI 章节。

本人课题的研究方向, 正是填补这个空缺: FEniCS 中实现 ALE FSI + 黏弹性 + G&R + 多物理耦合, 这是 2022 年以后血管生物力学数值方法的前沿。

7. 重要参考文献 (chapter-local)

注: 数字编号仅在本章内引用, 跨章不连续

  • [X1] Bathe K.J. (2006) Finite Element Procedures. Cambridge University Press — 引用 30, FEM 经典教材
  • [X2] Belytschko T., Hughes T.J.R. (eds) (1983) Computational Methods for Transient Analysis. North-Holland — 引用 34, 显式/隐式时间积分的原始论文集
  • [X3] Bonet J., Wood R.D. (2008) Nonlinear Continuum Mechanics for Finite Element Analysis. Cambridge University Press — 引用 53, 非线性 FEM 标准参考
  • [X4] Hughes T.J.R. (2000) The Finite Element Method: Linear Static and Dynamic Finite Element Analysis. Dover — 引用 269, FEM 标准教材
  • [X5] Zienkiewicz O.C., Taylor R.L., Zhu J.Z. (2013) The Finite Element Method: Its Basis and Fundamentals. Butterworth-Heinemann — 引用 617-619, FEM 三卷本
  • [X6] Donea J., Huerta A. (2003) Finite Element Methods for Flow Problems. Wiley — 引用 113, 流体力学 FEM
  • [X7] Brooks A.N., Hughes T.J.R. (1982) Streamline upwind/Petrov-Galerkin formulations for convection dominated flows with particular emphasis on the incompressible Navier-Stokes equations. Comput Methods Appl Mech Eng 32:199-259 — 引用 59, SUPG 原始论文
  • [X8] Hughes T.J.R., Mallet M. (1986) A new finite element formulation for computational fluid dynamics: III. The generalized streamline operator for multidimensional advective-diffusive systems. Comput Methods Appl Mech Eng 58:305-328 — 引用 268, GLS 原始论文
  • [X9] Riks E. (1979) An incremental approach to the solution of snapping and buckling problems. Int J Solids Struct 15:529-551 — 引用 460, Arc-length 法的原始论文
  • [X10] Crisfield M.A. (1983) An arc-length method including line searches and accelerations. Int J Numer Methods Eng 19:1269-1289 — 引用 106, 现代 arc-length 法的标准参考
  • [X11] Courant R., Friedrichs K., Lewy H. (1928) Über die partiellen Differenzengleichungen der mathematischen Physik. Math Ann 100:32-74 — CFL 条件的原始论文, 引用 100
  • [X12] Newmark N.M. (1959) A method of computation for structural dynamics. J Eng Mech Div ASCE 85:67-94 — Newmark 时间积分法
  • [X13] Chorin A.J. (1967) A numerical method for solving incompressible viscous flow problems. J Comput Phys 2:12-26 — 引用 87, Chorin 投影法
  • [X14] Simo J.C., Taylor R.L. (1991) Quasi-incompressible finite elasticity in principal stretches. Comput Methods Appl Mech Eng 85:273-310 — 引用 506, 不可压有限元
  • [X15] Holzapfel G.A., Gasser T.C., Ogden R.W. (2000) A new constitutive framework for arterial wall mechanics and a comparative study of material models. J Elast 61:1-48 — 引用 351, HGO 血管本构 (Ch 5-7 详述)
  • [X16] Gasser T.C., Auer M., Labruto F., Swedenborg J., Roy J. (2010) Biomechanical rupture risk assessment of abdominal aortic aneurysms: model complexity versus predictability. Eur J Vasc Endovasc Surg 40:176-185 — 引用 193, BRRA 临床验证
  • [X17] Quarteroni A., Valli A. (1994) Numerical Approximation of Partial Differential Equations. Springer — 引用 457, 数值 PDE 经典
  • [X18] Saad Y. (2003) Iterative Methods for Sparse Linear Systems. SIAM — 引用 140, 迭代求解器标准
  • [X19] Hughes T.J.R., Cottrell J.A., Bazilevs Y. (2005) Isogeometric analysis: CAD, finite elements, NURBS, exact geometry and mesh refinement. Comput Methods Appl Mech Eng 194:4135-4195 — 引用 269, 等几何分析
  • [X20] Beirão da Veiga L., Brezzi F., Cangiani A., Manzini G., Marini L.D., Russo A. (2013) Basic principles of virtual element methods. Math Models Methods Appl Sci 23:199-214 — 引用 12, VEM 原始论文