第3章：连续介质力学（Continuum Mechanics）

1. 作者

T. Christian Gasser (KTH Royal Institute of Technology) 独立撰写本章; Ch 3 是 Gasser 这本教科书中数学密度最高的一章, 涵盖运动学 (kinematics)、应力 (stress)、本构建模 (constitutive modeling)、守恒律 (governing laws)、一般原理 (general principles) 和损伤失效 (damage & failure) 六大模块。

章节范围: 3.1 Introduction → 3.9 Summary (PDF p113-p187, 75 页)

核心定位: Ch 5 (管道血管) 与 Ch 7 (血管壁, 含 G&R) 的数学基础。Ch 4 (有限元) 是 Ch 3 偏微分方程的数值解工具, 与 Ch 3 配套使用。

2. 内容概述

本章系统建立线性与非线性连续介质力学的数学框架: 从 Representative Volume Element (RVE) 的基本假设出发, 通过运动学 (kinematics) 描述物体的形变, 通过应力 (stress) 描述载荷, 通过本构方程 (constitutive law) 闭合守恒律方程, 最终得到 BVP (Boundary Value Problem) 的强形式与弱形式。

Gasser 的核心立场是: 血管生物力学的所有模型, 从最简单的线性弹性到最复杂的 G&R, 都必须建立在这套张量分析框架之上。掌握本章的张量语言 (F, C, σ, S, P, B, 等) 是阅读后续章节的前提。Ch 3 的 21 个 Example (占章节 ~30% 篇幅) 提供了"用张量解决具体问题"的标准工作流。

3. 核心方程与概念

3.1 RVE 与运动学 (Sect. 3.1-3.2)

RVE (Representative Volume Element): 描述材料连续介质性质的基本单元, 必须足够大以均质化材料性质, 足够小以保证宏观问题可解析。血管组织的多层级结构挑战这种尺度分离。

运动 \(\chi(X, t)\): 把材料粒子从参考构型 \(\Omega_0\) 的位置 \(X\) 映射到空间构型 \(\Omega\) 的位置 \(x(t) = \chi(X, t)\)。

3.1.1 形变梯度 (Eq. 3.1)

\[ \mathbf{F}(X,t) = \text{Grad}\chi(X,t) = \frac{\partial \chi(X,t)}{\partial X}, \quad F_{iJ} = \frac{\partial \chi_i}{\partial X_J} \]

\(\det \mathbf{F} > 0\) 防止自穿透
特征值表示: \(\mathbf{F} = \lambda_i \hat{n}_i \otimes \hat{N}_i\), \(\lambda_i\) = 主拉伸
\(\mathbf{F}\) 是两点张量 (two-point tensor): 下标 \(i\) 关联 \(\Omega\), 上标 \(J\) 关联 \(\Omega_0\)

3.1.2 极分解 (Polar Decomposition)

任意运动可分解为刚体旋转 + 形变: - 右极分解: \(\mathbf{F} = \mathbf{R}\mathbf{U}\), \(\mathbf{U}\) = 右拉伸张量 - 左极分解: \(\mathbf{F} = \mathbf{v}\mathbf{Q}\), \(\mathbf{v}\) = 左拉伸张量

\(\mathbf{R}, \mathbf{Q}\) 是旋转张量 (\(\mathbf{R}^{-1} = \mathbf{R}^T\), \(\det \mathbf{R} = +1\)); \(\mathbf{U}, \mathbf{v}\) 是对称正定拉伸张量。

3.1.3 线元/体元/面元变换 (Eq. 3.3-3.5)

线元: \(\lambda_a = dl/dL = \sqrt{\mathbf{A}:\mathbf{C}} = \sqrt{I_4}\), 其中 \(\mathbf{C} = \mathbf{F}^T \mathbf{F}\) 是右 Cauchy-Green 应变, \(\mathbf{A} = a_0 \otimes a_0\) 是结构张量, \(I_4\) 是"第四不变量"
体元: \(dv = J dV\), \(J = \det \mathbf{F}\) 是体积比
面元 (Nanson 公式): \(d\mathbf{s} = J \mathbf{F}^{-T} d\mathbf{S}\)

3.2 应变度量 (Sect. 3.2.7)

应变	公式	坐标系	类型	特征值
工程应变 \(\boldsymbol{\varepsilon}\) (Eq. 3.9)	\(\frac{1}{2}(\nabla \mathbf{u} + \nabla^T \mathbf{u})\)	空间	小变形	\(\varepsilon_i\)
右 Cauchy-Green \(\mathbf{C}\) (Eq. 3.10)	\(\mathbf{F}^T \mathbf{F}\)	参考	相对	\(\lambda_i^2\)
左 Cauchy-Green \(\mathbf{b}\) (Eq. 3.11)	\(\mathbf{F} \mathbf{F}^T\)	空间	相对	\(\lambda_i^2\)
Green-Lagrange \(\mathbf{E}\) (Eq. 3.12)	\(\frac{1}{2}(\mathbf{C} - \mathbf{I})\)	参考	绝对	\((\lambda_i^2-1)/2\)
Euler-Almansi \(\mathbf{e}\) (Eq. 3.13)	\(\frac{1}{2}(\mathbf{I} - \mathbf{b}^{-1})\)	空间	绝对	\((1-\lambda_i^{-2})/2\)

关键事实: \(\mathbf{C} = \mathbf{b} = \mathbf{I}\) 在无应变时; \(\mathbf{E} = \mathbf{e} = 0\) 在无应变时; 小变形极限下 \(\mathbf{E} \to \boldsymbol{\varepsilon}\), \(\mathbf{e} \to \boldsymbol{\varepsilon}\)。

纯剪切 vs 简单剪切 (Example 3.5): 在主拉伸坐标系下, 纯剪切是主轴方向的拉伸; 简单剪切是绕某轴旋转 \(\pi/4\) 后的剪切。两者小变形下给出相同 \(\mathbf{C}\), 但只有纯剪切没有刚体旋转。

3.3 应力度量 (Sect. 3.3)

3.3.1 Cauchy 应力定理 (Eq. 3.20)

\[ \mathbf{t}(x) = \boldsymbol{\sigma}(x)\mathbf{n}(x), \quad t_i = \sigma_{ia} n_a \]

Cauchy 应力 \(\boldsymbol{\sigma}\) 是真应力 (true stress), 作用于当前构型。角动量守恒要求 \(\boldsymbol{\sigma} = \boldsymbol{\sigma}^T\) (对称)。

3.3.2 应力变换

第一 Piola-Kirchhoff 应力 (Eq. 3.31): \(\mathbf{P} = J \boldsymbol{\sigma} \mathbf{F}^{-T}\), 是两点张量, 描述当前力作用于参考面积
第二 Piola-Kirchhoff 应力 (Eq. 3.32): \(\mathbf{S} = J \mathbf{F}^{-1} \boldsymbol{\sigma} \mathbf{F}^{-T}\), 是参考构型中对称二阶张量, 不受刚体旋转影响

关系: \(\mathbf{S} = \mathbf{F}^{-1} \mathbf{P}\) (Eq. 3.33), \(\mathbf{P} = \mathbf{F}\mathbf{S}\)。

3.3.3 主应力与 Mohr 圆 (Sect. 3.3.2-3.3.3)

通过特征值问题 \((\boldsymbol{\sigma} - \lambda \mathbf{I})\hat{n} = 0\) (Eq. 3.21) 求解主应力 \(\sigma_i\) 与主方向 \(\hat{n}_i\)。

应力不变量 (Eq. 3.22): - \(I_1 = \text{tr}\boldsymbol{\sigma} = \sigma_{11} + \sigma_{22} + \sigma_{33}\) - \(I_2 = \frac{1}{2}[(\text{tr}\boldsymbol{\sigma})^2 - \text{tr}\boldsymbol{\sigma}^2]\) - \(I_3 = \det\boldsymbol{\sigma}\)

Mohr 圆: 物理空间旋转 \(\alpha\) 对应应力空间旋转 \(2\alpha\); 用于 2D 应力状态分析。

3.3.4 八面体应力与 von Mises 应力 (Eq. 3.28)

八面体剪切应力 \(\tau_{\text{oct}}\) 是评估延性材料屈服风险的关键量:

\[ \sigma_M = \frac{3}{\sqrt{2}} \tau_{\text{oct}} = \frac{1}{\sqrt{2}}\sqrt{(\sigma_1-\sigma_2)^2 + (\sigma_2-\sigma_3)^2 + (\sigma_1-\sigma_3)^2} \]

在主应力空间中, \(\sigma_M\) 描述静水轴 \(\sigma_1 = \sigma_2 = \sigma_3\) 周围半径 \(r = \sqrt{2/3}\sigma_M\) 的圆柱面 (Fig. 3.16b)。

3.4 物质时间导数与目标导数 (Sect. 3.4)

速度梯度 (Eq. 3.34): \(\mathbf{l} = \nabla \mathbf{v} = \dot{\mathbf{F}} \mathbf{F}^{-1}\) - 分解: \(\mathbf{l} = \mathbf{w} + \mathbf{d}\) - 变形率张量 (Eq. 3.35): \(\mathbf{d} = \frac{1}{2}(\mathbf{l} + \mathbf{l}^T)\) - 自旋张量 (Eq. 3.36): \(\mathbf{w} = \frac{1}{2}(\mathbf{l} - \mathbf{l}^T)\)

目标导数问题: \(\dot{\boldsymbol{\sigma}}\) (Cauchy 应力的物质时间导数) 不是目标客观量 (Eq. 3.40); 替代: - Truesdell 率 (Eq. 3.41 下): \(\overset{\circ}{\boldsymbol{\sigma}} = \dot{\boldsymbol{\sigma}} - \mathbf{l}\boldsymbol{\sigma} - \boldsymbol{\sigma}\mathbf{l}^T + (\text{tr}\mathbf{l})\boldsymbol{\sigma}\) - Jaumann 率: \(\overset{\triangle}{\boldsymbol{\sigma}} = \dot{\boldsymbol{\sigma}} + \boldsymbol{\sigma}\mathbf{w} - \mathbf{w}\boldsymbol{\sigma}\)

3.5 功率共轭对 (Sect. 3.4.3)

应力	应变率
Cauchy \(\boldsymbol{\sigma}\)	变形率 \(\mathbf{d}\)
第二 P-K 应力 \(\mathbf{S}\)	\(\dot{\mathbf{C}}/2\) 或 Green-Lagrange 率 \(\dot{\mathbf{E}}\)
第一 P-K 应力 \(\mathbf{P}\)	\(\dot{\mathbf{F}}\)

变形功率: \(\dot{\psi} = \boldsymbol{\sigma} : \mathbf{d}\) (每单位变形体积), \(\dot{\mathcal{W}} = \int_{\Omega} \boldsymbol{\sigma} : \mathbf{d}\, dv\)。

3.6 本构建模 (Sect. 3.5)

3.6.1 线性弹性 (Hooke, Eq. 3.51)

\[ c_{ijkl} = \frac{E\nu}{(1+\nu)(1-2\nu)}\delta_{ij}\delta_{kl} + \frac{E}{2(1+\nu)}(\delta_{ik}\delta_{jl} + \delta_{il}\delta_{jk}) \]

各向同性线弹性由 2 个独立参数描述: \((E, \nu)\) 或 \((G, K)\) 或 \((E, G)\)。

关系: \(E = 2(1+\nu)G\) (Eq. 3.46), \(G = E/[2(1+\nu)]\), \(K = E/[3(1-2\nu)]\)。

刚度张量性质: minor 对称 \(c_{ijkl} = c_{jikl} = c_{ijlk}\); major 对称 \(c_{ijkl} = c_{klij}\) (弹性体); 正定 \(\Delta\boldsymbol{\varepsilon} : \mathbf{C} : \Delta\boldsymbol{\varepsilon} > 0\)。

3.6.2 超弹性 (Hyperelasticity, Sect. 3.5.3)

Helmholtz 自由能 \(\Psi(\mathbf{C}, \mathbf{A}_i) = \Psi(I_1, I_2, I_3, \ldots, I_m)\), 其中 \(I_i\) 由 \(\mathbf{C}\) 和结构张量 \(\mathbf{A}_i\) 形成。

Coleman-Noll 程序 (Eq. 3.129):

\[ \mathbf{S} = 2\frac{\partial \Psi(\mathbf{C}, G_1, \ldots, G_n)}{\partial \mathbf{C}} \]

体积-等容分解 (Eq. 3.55, 3.56): - 运动分解: \(\chi(X) = \bar{\chi}(\chi_{\text{vol}}(X))\) - 形变梯度: \(\bar{\mathbf{F}} = J^{-1/3} \mathbf{F}\), \(\mathbf{F}_{\text{vol}} = J^{1/3}\mathbf{I}\) - 自由能: \(\Psi(J, \bar{\mathbf{C}}) = \Psi^{\text{iso}}(\bar{\mathbf{C}}) + \Psi^{\text{vol}}(J)\) - 应力: \(\mathbf{S} = 2J^{-2/3}\text{Dev}\left[\frac{\partial \Psi^{\text{iso}}(\bar{\mathbf{C}})}{\partial \bar{\mathbf{C}}}\right] - Jp\mathbf{C}^{-1}\), 其中 \(\text{Dev}(\bullet) = (\bullet) - [\mathbf{C}:(\bullet)]\mathbf{C}^{-1}/3\)

不可压材料 (Eq. 3.131): \(\Psi^{\text{vol}}\) 退化为 Lagrange 乘子 \(\kappa(J-1)\):

\[ \mathbf{S} = 2\text{Dev}\left[\frac{\partial \Psi^{\text{iso}}(\mathbf{C})}{\partial \mathbf{C}}\right] - p\mathbf{C}^{-1} \]

重要区分: \(\kappa\) (Lagrange 参数) ≠ \(p\) (静水压) — 严格不可压用 \(p\) 替换 \(\kappa\), 但两者物理意义不同。

3.6.3 黏弹性 (Viscoelasticity, Sect. 3.5.4)

Boltzmann 叠加原理 (Eq. 3.62):

\[ \varepsilon(t) = \int_{-\infty}^{t} J(t-\xi)\dot{\sigma}(\xi)\, d\xi \]

其中 \(J\) 是蠕变柔量, \(G\) 是松弛模量。

Maxwell 单元 (Eq. 3.64): 弹簧 + 阻尼器串联

\[ \frac{d\varepsilon}{dt} = \frac{1}{E}\frac{d\sigma}{dt} + \frac{\sigma}{\eta} \]

- 松弛时间: \(\tau = \eta/E\) - 应力响应: \(\sigma(t) = \Delta\sigma \exp(-t/\tau)\) (Eq. 3.66) - 极限行为: \(\tau \ll t\) → 黏性流体; \(\tau \gg t\) → 弹性固体

Kelvin-Voigt 单元 (Eq. 3.73): 弹簧 + 阻尼器并联

\[ \sigma = E\varepsilon + \eta\frac{d\varepsilon}{dt} \]

- 延迟时间 \(\tau = \eta/E\) - 蠕变响应: \(\varepsilon(t) = (\Delta\sigma/E)(1 - e^{-t/\tau})\)

Standard Solid 单元 (Eq. 3.74): 主弹簧 + Maxwell 单元并联

\[ \frac{d\varepsilon}{dt} = \frac{1}{E_M + E_e}\left[\frac{d\sigma}{dt} + \frac{E_M}{\eta}(\sigma - E_e \varepsilon)\right] \]

黏超弹 (Visco-hyperelastic) Standard Solid (Eq. 3.76-3.77): 隐藏内部变量 \(\mathbf{C}_M\) (Maxwell 弹簧的右 Cauchy-Green 应变)

\[ \mathbf{S} = 2\frac{\partial \Psi^{\text{iso}}_E(\mathbf{C})}{\partial \mathbf{C}} + 2\frac{\partial \Psi^{\text{iso}}_M(\mathbf{C}_M)}{\partial \mathbf{C}_M} : \frac{\partial \mathbf{C}_M}{\partial \mathbf{C}} - \kappa \mathbf{C}^{-1} \]

\[ \dot{\mathbf{E}}_M + \frac{1}{\tau}\mathbf{E}_M = \dot{\mathbf{E}} \]

Example 3.15 (不可压 neo-Hookean): \(\Psi^{\text{iso}}(\mathbf{C}, \mathbf{C}_M) = \frac{G}{2}[(I_1 - 3) + \beta(I_{1M} - 3)]\)

Example 3.16 (应力分解法): \(\mathbf{S} = \mathbf{S}_E + \mathbf{S}_M - p\mathbf{C}^{-1}\), 速率方程 \(\dot{\mathbf{S}}_M + \mathbf{S}_M/\tau = \beta \dot{\mathbf{S}}_E\) (Eq. 3.89) — 数值实现比能量分解法更高效 (引用 [194, 222, 291, 506])。

3.6.4 多相连续介质 (Sect. 3.5.5)

混合理论 (Mixture Theory): n 相同时占据 RVE 同一位置; 用 Fick 扩散定律描述相间传输; 应变能 \(\Psi = \sum_i \xi_i \Psi_i\), \(\sum \xi_i = 1\) (Eq. 3.98)
多孔弹性 (Poroelasticity, Biot 1941): 固相 (骨架) + 液相占据 RVE 分离区域; Darcy 定律 \(q = -L_p \nabla p - \mathbf{b}_f\) 描述流体过滤; 应力分解 \(\boldsymbol{\sigma} = \boldsymbol{\sigma}^e - \alpha p \mathbf{I}\) (Eq. 3.99)

区别: 混合理论更一般 (允许相间耦合), 多孔弹性更适合结构-流体耦合 (如心肌)。

3.7 守恒律 (Sect. 3.6)

3.7.1 质量守恒 (Eq. 3.102, 3.104)

Lagrangian: \(\frac{\partial \rho}{\partial t} + \rho \text{div}\mathbf{v} = 0\) Eulerian 不可压: \(\text{div}\mathbf{v} = 0\)

3.7.2 动量守恒 (Cauchy 动量方程, Eq. 3.107)

\[ \rho \frac{D\mathbf{v}}{Dt} = \text{div}\boldsymbol{\sigma} + \mathbf{b}_f \]

Navier-Stokes 方程 (Eq. 3.115, 代入 Newtonian 黏度):

\[ \rho\left(\frac{\partial \mathbf{v}}{\partial t} + \mathbf{v} \cdot \nabla\mathbf{v}\right) = \eta \text{div}(\nabla\mathbf{v}) - \nabla p + \mathbf{b}_f \]

无量纲形式 (Eq. 3.116): 引入 Strouhal \(\text{Sr} = L/(TV)\), Reynolds \(\text{Re} = VL\rho/\eta\), Froude \(\text{Fr} = V/\sqrt{gL}\) 数, 评估各项相对重要性。

3.7.3 热力学定律 (Sect. 3.6.3)

第一定律 (Eq. 3.121): \(\dot{u} + \text{div}\mathbf{q}_h - r_h - \boldsymbol{\sigma}:\mathbf{d} = 0\)

Clausius-Duhem 不等式 (Eq. 3.124):

\[ -\dot{\psi} - s\dot{\theta} + \boldsymbol{\sigma}:\mathbf{d} - \frac{\mathbf{q}_h \cdot \nabla\theta}{\theta} \geq 0 \]

材料形式 (Eq. 3.125):

\[ -\dot{\Psi} - S\dot{\theta} + \mathbf{P}:\dot{\mathbf{F}} - \frac{\mathbf{Q}_h \cdot \text{Grad}\theta}{\theta} \geq 0 \]

3.7.4 Coleman-Noll 程序 (Eq. 3.127-3.129)

从自由能 \(\Psi(\mathbf{F}, H_1, \ldots, H_n)\) 导出本构关系:

\[ \mathbf{P} = \frac{\partial \Psi}{\partial \mathbf{F}}, \quad \mathbf{S} = 2\frac{\partial \Psi}{\partial \mathbf{C}} \]

并保证耗散不等式 \(-L_i : \dot{H}_i \geq 0\)。

3.8 一般原理 (Sect. 3.7)

3.8.1 自由体图 (FBD, Sect. 3.7.2)

薄壁球面 (Eq. 3.144): \(\sigma = \frac{rp_i}{2h}\) (静定问题, 不依赖材料性质) 薄壁管 (Laplace 定律): \(\sigma_\theta = \frac{rp_i}{h}\) (Example 3.23)

3.8.2 边界值问题 (BVP, Eq. 3.146, 3.147)

固体力学 (Lagrangian) iBVP:

\[ \text{div}\boldsymbol{\sigma} + \mathbf{b}_f = \rho \frac{D^2\mathbf{u}}{Dt^2} \text{ in } \Omega \]

\[ \mathbf{u} = \bar{\mathbf{u}} \text{ at } \partial\Omega_u, \quad \mathbf{t} = \bar{\mathbf{t}} \text{ at } \partial\Omega_t \]

\[ \frac{D\mathbf{u}}{Dt} = \mathbf{v}_0, \mathbf{u} = \mathbf{u}_0 \text{ in } \Omega \text{ at } t = t_0 \]

流体力学 (Eulerian) iBVP: 主变量是 \(\mathbf{v}\) 而非 \(\mathbf{u}\)。

3.8.3 虚功原理 (PVW, Eq. 3.148-3.150)

小变形 (Eq. 3.148): \(\int_\Omega \boldsymbol{\sigma}:\delta\boldsymbol{\varepsilon}\, dv = \int_\Omega \mathbf{b}_f \cdot \delta\mathbf{u}\, dv + \int_{\partial\Omega} \mathbf{t} \cdot \delta\mathbf{u}\, ds\)

有限变形空间形式 (Eq. 3.149): \(\int_\Omega \boldsymbol{\sigma}:\text{grad}_s \delta\mathbf{u}\, dv = \ldots\)

有限变形材料形式 (Eq. 3.150): \(\int_{\Omega_0} \mathbf{S}:\delta\mathbf{E}\, dV = \ldots\) — FEM 的出发点 (Ch 4)

3.9 损伤与失效 (Sect. 3.8)

3.9.1 连续损伤力学 (CDM, Eq. 3.151)

Kachanov 型: \(\boldsymbol{\sigma} = (1-D)\boldsymbol{\sigma}^{\text{eff}}\), \(D\) = 标量损伤参数 (各向同性) 或二阶张量 (各向异性)。

关键物理效应: \(D\) 增长 → 应变软化 → 刚度张量不再正定 → 强椭圆性丧失 → 失去波传播能力。

3.9.2 应变局部化 (Sect. 3.8.2)

非极性连续介质 (non-polar continuum) 在应变软化下, 应变局部化会收敛到无穷小体积且不耗散能量 — 物理上错误。需要正则化: - 梯度增强损伤模型 [425]: \(D\) 演化受邻域变形影响 (非局部) - 黏性正则化: 用 Hreg = H/n 替代 H, 引入网格依赖

关键观察 (Example 3.25): 双线性应力-应变材料在拉伸时, 平均应变 \(\bar{\varepsilon}\) 依赖于截面数量 \(n\) (网格依赖), 这是经典的"数值伪依赖", 需要正则化。

3.9.3 线弹性断裂力学 (LFM, Sect. 3.8.3)

应力强度因子 \(K_I, K_{II}, K_{III}\) [Pa m^(1/2)]: 描述裂尖应力场 能量释放率 \(D = K_I^2(1-\nu^2)/E + K_{II}^2(1-\nu^2)/E + K_{III}^2(1+\nu)/E\)

限制: 需要小变形 + 锐裂尖 → 不适用于血管组织 (非线弹性 + 钝裂尖)。

3.9.4 J-积分 (Eq. 3.155)

\[ J = \int_\Gamma \left[\psi\, dx_2 - \mathbf{t} \cdot \frac{\partial \mathbf{u}}{\partial x_1}\, ds\right] = D \]

路径无关: 任意闭合路径 \(\Gamma\) 给出相同 \(J\) 值, \(J\) 等于能量释放率。

3.9.5 内聚区模型 (Cohesive Zone, Eq. 3.156)

Traction Separation Law (TSL): 通过内聚势 \(\Psi^c(u_d, \zeta)\) 描述裂尖阻力:

\[ \mathbf{T} = \frac{\partial \Psi^c}{\partial u_d}, \quad D = -\frac{\partial \Psi^c}{\partial \zeta}\dot{\zeta} \geq 0 \]

应用: 金属 [26], 混凝土 [254], 生物组织 [200, 233]。

4. 关键结论

[K1] RVE 是连续介质力学的基础: 尺度分离假设 (宏观-微观) 是所有本构模型的合法性前提; 血管组织的多层级结构挑战此假设
[K2] 形变梯度 \(\mathbf{F}\) 包含刚体运动 + 形变: 极分解 \(\mathbf{F} = \mathbf{R}\mathbf{U}\) 把两者分离, 物理观测只依赖 \(\mathbf{U}\) (右拉伸张量)
[K3] Green-Lagrange \(\mathbf{E}\) 在大变形下保持目标客观: 这是大变形本构建模的首选应变度量; 工程应变 \(\boldsymbol{\varepsilon}\) 仅适用于小变形
[K4] 第二 Piola-Kirchhoff 应力 \(\mathbf{S}\) 是大变形本构的"配对": \(\mathbf{S} = 2\partial\Psi/\partial\mathbf{C}\), 与 \(\mathbf{C}\) 配对形成功率共轭
[K5] von Mises 应力 \(\sigma_M\) 是延性材料屈服判据: 静水轴周围的圆柱面 (Fig. 3.16b) 是屈服的几何图像
[K6] Coleman-Noll 程序从自由能 \(\Psi\) 自动导出本构: \(\mathbf{S} = 2\partial\Psi/\partial\mathbf{C}\) + 耗散不等式 \(-L_i:\dot{H}_i \geq 0\)
[K7] 不可压材料需要 Lagrange 乘子 \(\kappa\) (或 \(p\)): 但不等于静水压; 在 decoupled 形式中 \(p\) 才是静水压
[K8] Maxwell 单元的 \(\tau\) 决定材料行为: \(\tau \ll t\) → 流体; \(\tau \gg t\) → 弹性固体; 真实材料是多个 Maxwell 单元的叠加
[K9] Standard Solid 是血管壁黏弹性的最小模型: 主弹簧 + Maxwell 单元并联, 既有瞬时弹性又有长期松弛; \(\beta\) 调节 Maxwell 与主弹簧的刚度比
[K10] Navier-Stokes 方程的无量纲化揭示流动的物理本质: Reynolds 数控制层流-湍流, Strouhal 数控制瞬态, Froude 数控制重力
[K11] 虚功原理是 FEM 的基础: 材料形式 (Eq. 3.150) 直接对应 Ch 4 的有限元离散
[K12] 应变软化在非极性连续介质中是非物理的: 需要正则化 (梯度增强, 黏性) 才能得到物理上合理的结果
[K13] LFM 不适用于血管组织: 钝裂尖 + 非线弹性 → 必须用 J-积分或 cohesive zone
[K14] 多孔弹性是血管壁内流体-结构耦合的标准工具: Biot 1941 理论在心肌建模 [77] 中仍活跃

5. 挑战和开放性问题

[C1] 不可压与近不可压的数值困难: \(J = 1\) 严格约束导致 BVP 病态, 需要混合有限元 (u/p) 或罚方法; Gasser 在 Example 3.21 简单演示, 但 Ch 4 数值实现细节未充分展开
[C2] 大变形各向异性本构的参数化困难: 血管组织是各向异性的 (纤维方向), 需要结构张量 \(\mathbf{A}_i\) 描述; Gasser 引用 [497] 的"a priori polyconvex"框架, 但 21 个独立分量的全各向异性弹性张量在工程实践中几乎不可标定
[C3] 内变量(hidden variables)的选择: 黏超弹中 \(\mathbf{C}_M\) 是 Maxwell 弹簧的右 Cauchy-Green 应变, 但选择有多种 (Simo 1987, Holzapfel 2002), 不同选择给出不同本构; Gasser 引用 [194, 222, 291, 506] 但未给出选择标准
[C4] 应变软化与网格依赖: Example 3.25 显示平均应变 \(\bar{\varepsilon}\) 依赖 \(n\) (截面数), 这是经典的"病态问题"; Gasser 提到 Hreg = H/n 黏性正则化, 但梯度增强损伤 (gradient-enhanced damage) 是更严格的方案, Gasser 引用 [415, 425] 但实现细节缺位
[C5] 多孔弹性的有限变形推广: 3.5.5.2 提到 Biot 理论可推广到有限变形, 但只给出 "Jaumann rate 替代 Cauchy rate" 的指示; 实际实现涉及耦合 BVP 的非线性迭代, Gasser 应当给出 reference
[C6] LFM 在血管生物力学的应用边界: 3.8.3 坦承 LFM 要求小变形+锐裂尖, 但AAA 实际是钝裂尖+大变形; Gasser 引用 [200, 233] 的 cohesive zone 工作, 但具体的 TSL 标定方法 (单轴拉伸 vs 双轴 vs 原位) 缺位
[C7] J-积分路径无关性的数学基础: 3.8.4.1 说"J 积分在闭合路径下路径无关", 但证明需要单连通域+无奇点; 血管组织的各向异性+残余应力使单连通假设复杂化
[C8] 混合理论 vs 多孔弹性的选择标准: Gasser 说"混合理论更一般", 但工程实践中几乎全部用多孔弹性; 何时用混合理论? Gasser 未给指引
[C9] 有限变形 PVW 的空间形式 vs 材料形式: 3.7.4.2 给出空间形式 (Eq. 3.149) 和材料形式 (Eq. 3.150), 但实际 FEM 实现用材料形式; 为什么空间形式仍然重要? (例如, contact 问题中空间形式更直观)
[C10] Cauchy 应力 \(\boldsymbol{\sigma}\) 在大变形下"不是客观量": Gasser 解释 \(\dot{\boldsymbol{\sigma}}\) 非客观, 但 \(\boldsymbol{\sigma}\) 本身是张量, 总是客观的; 这是"应力率" vs "应力"的区别, Gasser 的措辞有歧义
[C11] Maxwell 单元的零耗散问题: Gasser 说"非常高与非常低应变率下不耗散", 但实际材料永远有耗散; 修正需要 fractional Zener 模型 (黏壶用分数阶导数), Gasser 未提
[C12] 2022 年的盲点: Physics-informed neural networks (PINN) 求解 Navier-Stokes 与 BVP, 数据驱动的本构建模 (如 auto-differentiable FEM) 正在颠覆传统有限元; Gasser 应当提及

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

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

本章是本人 G&R 课题的数学支柱:

\(\mathbf{F}, \mathbf{C}, \mathbf{S}, \mathbf{E}\) 符号 在本人 FEniCS 代码中直接使用 (例如 F = Identity(3) + grad(u), C = F.T*F, S = 2*diff(Psi, C))
Coleman-Noll 程序 (Eq. 3.127-3.129) 是本人 SMC G&R 中隐藏内变量演化 (如 \(\mathbf{C}_M\)) 的标准推导工具
不可压处理 (Eq. 3.131-3.134): 本人用 FEniCS 求解血管壁时, 必须用混合空间 Taylor-Hood (P2/P1), 即 \(u\) 用 P2 (二次), \(p\) 用 P1 (一次)
Standard Solid 模型 (Sect. 3.5.4.5): 本人研究血管壁黏弹性的最简模型; 参数 \(\tau, \beta\) 与实验拉伸-松弛数据直接拟合
J-积分 (Eq. 3.155): 本人课题"AAA 破裂"判据的候选; 与 cohesive zone 模型互补

6.2 Gasser 框架的局限性

张量密度过高: Ch 3 包含 21 个 Example + 大量公式, 对初学者非常陡峭; Gasser 引用 [53, 185, 271, 354, 401] 给出了 5 本经典连续介质教材, 但没给"最小路径"。对血管生物力学背景的学生, 更需要的是"哪些公式是必背的" 而不是"全部罗列"
3D 与 2D 简化的不清晰: 实际血管生物力学 80% 是轴对称 (cylindrical) 或广义平面应变 (generalized plane strain); Gasser 在 Example 3.19 (厚壁圆柱) 给出轴对称公式, 但广义平面应变未提及 (这在 G&R 模拟中常用)
有限旋转与有限应变的区别未明: 3.2.7 提到"小形变+大旋转"用"有限旋转分析", 但具体什么公式? Gasser 引用 [354] 留给读者
体积-等容分解的"工程实践": Gasser 在 3.5.3.2 给出标准分解, 但未给出不可压近似的稳定性讨论; 在大变形下, decoupled formulation 的准确度依赖于参考构型的选择 (Gasser 引用 [519] 但未深入)
多孔弹性的可压缩性: 3.5.5.2 说"可推广到有限变形"但没有给出大变形下的流体本构; 实际血管壁内流体 (水 + 白蛋白 + 脂质) 是非牛顿 + 黏弹性, 远复杂于 Darcy 定律
2022 年的方法学转向: Physics-Informed Neural Networks (PINN) 在求解 Navier-Stokes 与 BVP 上已取得突破, 自动微分 (JAX/PyTorch) 改变了本构建模的工作流; Gasser 在 3.9 Summary 仅提及 "Machine Learning (ML)" 一句, 深度不足

6.3 一个观察: "本构的物理 vs 唯象" 张力

Gasser 在 3.5 节给出唯象本构 (Yeoh, neo-Hookean 等) 与基于微观机制的本构 (collagen fiber 方向, SMC 收缩) 的对比。Ch 5-7 大量使用唯象本构 (Gasser 引用 [195, 198, 351]), 但血管生物力学的真正难题是微观机制如何涌现为宏观本构。

本人研究的 G&R 框架, 正是试图在唯象本构与微观机制之间建立桥梁: G&R 律 (mass 平衡) 给出本构的时间演化, 而本构的时间演化又改变应力场。这是一个双向耦合问题, Gasser 在 Ch 3 仅给出静态本构, 没有 G&R 的动态方程。这正是 Ch 7 要补的。

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

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

[X1] Holzapfel G.A. (2000) Nonlinear Solid Mechanics: A Continuum Approach for Engineering. John Wiley & Sons — 引用 271, 整本 Ch 3 的基础教材
[X2] Bonet J., Wood R.D. (2008) Nonlinear Continuum Mechanics for Finite Element Analysis. Cambridge University Press — 引用 53, 经典的有限元导向连续介质
[X3] Simo J.C., Hughes T.J.R. (1998) Computational Inelasticity. Springer — 引用 506, 不可压黏超弹的奠基工作
[X4] Simo J.C., Taylor R.L. (1985) Consistent tangent operators for rate-independent elastoplasticity. Comput Methods Appl Mech Eng 48:101-118 — 引用 506, consistent tangent 算子的奠基
[X5] Gasser T.C., Holzapfel G.A. (2002) A rate-independent elastoplastic constitutive model for (fibrous) biological tissues. Comput Mech 29:340-360 — 引用 194, 黏超弹应力分解法的原始论文之一
[X6] Holzapfel G.A., Gasser T.C. (2001) A viscoelastic model for fiber-reinforced composites at finite strains: Continuum basis, computational aspects and applications. Comput Methods Appl Mech Eng 190:4379-4403 — 引用 222, 黏超弹的有限元实现
[X7] Miehe C. (1996) Numerical computation of algorithmic (consistent) tangent moduli in large-strain computational inelasticity. Comput Methods Appl Mech Eng 134:223-240 — 引用 291
[X8] Coleman B.D., Noll W. (1963) The thermodynamics of elastic materials with heat conduction and viscosity. Arch Ration Mech Anal 13:167-178 — 引用 3.127, Coleman-Noll 程序的原始论文
[X9] Biot M.A. (1941) General theory of three-dimensional consolidation. J Appl Phys 12:155-164 — 多孔弹性理论的原始论文
[X10] Rice J.R. (1968) A path independent integral and the approximate analysis of strain concentration by notches and cracks. J Appl Mech 35:379-386 — J-积分的原始论文, 引用 460
[X11] Kachanov L.M. (1958) Time of the rupture process under creep conditions. Izv Akad Nauk SSSR Otd Tech Nauk 8:26-31 — 连续损伤力学的原始论文, 引用 290
[X12] Peerlings R.H.J., de Borst R., Brekelmans W.A.M., de Vree J.H.P. (1996) Gradient enhanced damage for quasi-brittle materials. Int J Numer Methods Eng 39:3391-3403 — 梯度增强损伤, 引用 425
[X13] Hillerborg A., Modéer M., Petersson P.E. (1976) Analysis of crack formation and crack growth in concrete by means of fracture mechanics and finite elements. Cement Concrete Res 6:773-782 — Cohesive zone 模型的原始论文, 引用 254
[X14] Tvergaard V., Hutchinson J.W. (1992) The relation between crack growth resistance and fracture process parameters in elastic-plastic solids. J Mech Phys Solids 40:1377-1397 — Cohesive zone 在金属中的应用, 引用 26
[X15] Horgan C.O., Saccomandi G. (2003) A molecular-statistical basis for some Gent strain energy functions. Int J Non-Linear Mech 38:1029-1036 — 引用 3.5.3.1 的橡胶分子统计基础
[X16] Belytschko T., Liu W.K., Moran B., Elkhodary K. (2013) Nonlinear Finite Elements for Continua and Structures. John Wiley & Sons — 引用 45, 经典的非线性有限元教材
[X17] Curnier A. (1999) Computational Methods in Solid Mechanics. Kluwer — 引用 102, 混合理论参考
[X18] Ogden R.W. (1997) Non-linear Elastic Deformations. Dover — 引用 354, 非线性弹性经典教材
[X19] Marsden J.E., Hughes T.J.R. (1983) Mathematical Foundations of Elasticity. Prentice-Hall — 引用 185, 严谨的连续介质数学基础
[X20] Fung Y.C. (1993) Biomechanics: Mechanical Properties of Living Tissues. Springer — 引用 256, 生物力学经典 (与血管生物力学直接对接)