跳转至

第四章 有限弹性 (Finite Elasticity)

§1 作者

本章是把Ch 3抽象的连续介质框架具体化为可计算本构的关键章。Humphrey的策略是"以橡胶 (elastomer) 为模板"——橡胶和软组织都是长链高分子，前者的理论是后者的模板。作者特别强调：4.1节是Ch 7、10的"理论字典"，本构的物理意义、实验验证、与微观机制的桥接都在这一章展开。

§2 内容概述

分六节： - 4.1 不可压缩各向同性弹性：DEICE法在橡胶上的完整示例；从共性观察 (非线性、滞回、预调) → 理论框架 (\(W = W(I_C, II_C)\) + Lagrange 乘子) → 4个具体本构 (neo-Hookean, Mooney-Rivlin, Rivlin-Sawyers, Hart-Smith/Crisp 指数) → 参数拟合 → 预测能力。 - 4.2 三维不可压缩弹性解：单轴/双轴拉伸、圆杆扭转、厚壁圆管 (Ch 7血管) 三个通用解 (universal solution)。 - 4.3 可压缩各向同性弹性：可压缩假设下的额外不变量依赖 + Blatz-Ko 泡沫橡胶型。 - 4.4 膜超弹性：当 \(h \ll\) 其他尺寸时 \(\sigma_{33} = 0\) + 平面应力；Ch 7血管壁/心外膜用此。 - 4.5 习题：15题。 - 4.6 参考文献 (53条)。

§3 核心方程与概念

3.1 不可压缩各向同性本构 (Eq. 4.1—4.10)

(1) 不可压缩各向同性W (Eq. 4.1) $\(W = W(I_C, II_C) - p(III_C - 1)\)$ 代入Ch 3.66, 用Cayley-Hamilton化简 (Ch 2.89) 得： $\(\boldsymbol{\sigma} = -p\mathbf{I} + 2W_1 \mathbf{B} - 2W_2 \mathbf{B}^{-1}, \quad W_1 = \frac{\partial W}{\partial I_C}, \quad W_2 = \frac{\partial W}{\partial II_C}\)$ 意义：本构只剩 2个标量响应函数 \(W_1, W_2\)。Humphrey明确说"两者数值上可由实验提取"——这是Ch 7、10的数据接口。

(2) Neo-Hookean (Eq. 4.4—4.5) $\(W = c_1(I_C - 3), \quad c_1 = nkT\)$ 源于纯熵统计力学 (Treloar 1975)：\(n\) 链/单位体积, \(k\) Boltzmann常数, \(T\) 温度。预测：\(\sigma\) 与 \(\mathbf{B}\) 线性。适用范围：拉伸 ≤ 30%。

(3) Mooney (Eq. 4.6) $\(W = c_1(I_C - 3) + c_2(II_C - 3)\)$ 1940年Mooney现象学假设。目标：拉伸下非线性、剪切下线性（观察到的特性）。比 neo-Hookean 描述范围更大。常被引用为"Mooney-Rivlin"（虽然Rivlin更晚证明它的形式）。

(4) Rivlin-Sawyers (Eq. 4.7) $\(W = c_1(I_C - 3) + f(II_C), \quad f = c_2 + c_3(II_C - 3), \quad c_1, c_2 > 0, c_3 < 0\)$ Rivlin在1951年用不变量的"universal solution"提取 \(W_1, W_2\)，到1976年 (Rivlin-Sawyers)才提出具体函数形式。E-不等式 (Ch 3.4.4) 给出 \(c_1, c_2 > 0, c_3 < 0\)。

(5) Hart-Smith-Crisp 指数型 (Eq. 4.8) $\(W_1 = c\exp[c_1(I_C - 3)^2], \quad W_2 = II_C^{-1}\)$ Fung学派指数型本构的"远亲"。Fung 1967 软组织型的 \(W = c(e^Q - 1)\) 形式是这个思路的拓展，Ch 7、10血管/心脏用。

(6) Valanis-Landel可分型 (Eq. 4.9) $\(W(\lambda_1, \lambda_2, \lambda_3) = f(\lambda_1) + f(\lambda_2) + f(\lambda_3)\)$ 把 \(W\) 写成主伸长 \(\lambda_I\) 的可分函数——比 \(W(I_C, II_C)\) 更易实现。条件：材料各向同性 + 主方向与主应力方向一致。Carmichael-Holdaway 1961 早期提出，Valanis-Landel 1967 完善。

(7) Ogden (Eq. 4.10) $\(W = \sum_{p=1}^{N} \frac{\mu_p}{\alpha_p}\left(\lambda_1^{\alpha_p} + \lambda_2^{\alpha_p} + \lambda_3^{\alpha_p} - 3\right)\)$ \(\mu_p\) 应力量纲，\(\alpha_p\) 无量纲，不必为整数。可恢复 Mooney-Rivlin (\(N=2, \alpha=2,-2\))。Treloar 1944 单轴数据用 \(N=3\) 拟合极佳 (Fig. 4.1, \(N=4\) 略好但参数多了 2 个)。Ch 10心肌 (Bogen-McMahon 1979) 用 Ogden 形式。

3.2 通用解 (Eq. 4.11—4.50)

(8) 双轴拉伸均匀变形 (Eq. 4.11—4.17) $\(x_1 = \lambda_1 X_1, \quad x_2 = \lambda_2 X_2, \quad x_3 = \lambda_3 X_3, \quad \lambda_1 \lambda_2 \lambda_3 = 1\)$ 应力 (Eq. 4.15—4.17)： $\(t_{11} = 2(\lambda_1^2 - \lambda_3^2)(W_1 + \lambda_1^2 W_2), \quad t_{22} = 2(\lambda_2^2 - \lambda_3^2)(W_1 + \lambda_2^2 W_2)\)$ 关键 (Eq. 4.20)：联立解出 \(W_1, W_2\)： $\(W_1 = \frac{\lambda_2^2 t_{11} - \lambda_1^2 t_{22}}{2(\lambda_1^2 - \lambda_2^2)}, \quad W_2 = \frac{t_{11} - t_{22}}{2(\lambda_1^2 - \lambda_2^2)}\)$ 意义：双轴实验 = 实验提取 \(W_1, W_2\) 的唯一方法；单轴只有一个方程 (Eq. 4.26)，无法唯一确定两个响应函数。

(9) 单轴拉伸 (Eq. 4.26) $\(\lambda_1 = \lambda_2 = \lambda^{-1/2}, \quad t_{11} = t_{22} = 0, \quad t_{33} = 2(\lambda^2 - \lambda^{-1})(W_1 + \lambda^{-1} W_2)\)$ 单轴 = 1个方程 + 2个未知响应函数 → 不能同时确定 \(W_1, W_2\)。这是Ch 5实验设计选择"双轴 vs 单轴"的核心理由。

(10) 圆杆扭转+拉伸 (Eq. 4.28—4.50) 用复合变形 \(\mathbf{F} = \mathbf{F}_2 \cdot \mathbf{F}_1\) (先拉后扭) 得到 \(t_{rr}, t_{\theta\theta}, t_{zz}, t_{\theta z}\)，并用径向平衡 (Eq. 4.41) 数值积分求 \(p(r)\)。Eq. 4.54—4.55 给出 \(W_1, W_2\) 的实验提取公式 (类似双轴)。扭转实验 + \(M(\pm y)\) 对称性可验证各向同性假设。

(11) 厚壁圆管 (Eq. 4.56—4.58) $\(r = r(R), \quad \theta = \Theta + \gamma Z, \quad z = \Lambda Z\)$ 不可压缩 \(\Rightarrow\) \(r^2 = R^2/\Lambda + C\) (Eq. 4.58)。意义：Ch 7血管残余应力分析 (开环/闭环) 用的就是这套，只是 \(\gamma=0, \Lambda \neq 1\)。圆管 + 扭转同时还提供 \(\gamma\) 与血流剪切相关——Ch 7用。

3.3 可压缩 + 膜 (Eq. 4.59—4.66)

(12) 可压缩各向同性 (Ch 4.3) $\(W = W(I_C, II_C, III_C), \quad III_C = J^2 \neq 1\)$ 典型 Blatz-Ko (1962) 泡沫橡胶型 (Ch 4.3节)。

(13) 膜超弹性 (Ch 4.4) 薄壁 (\(h \ll R\)) 时 \(t_{33} = 0\) 自动满足 + 平面应力，但 \(W\) 还是依赖全部 \(I_C, II_C\) (3D不变量)。Eq. 4.59附近展示 \(C\) 矩阵 (在物理分量下) 的 \(3\times 3\) 表达式。Ch 7血管壁 (动脉壁 \(h/R \sim 0.05-0.15\)) 是这一节的直接应用。

§4 关键算法或建模方法

  1. 不变量的"通用解" → 响应函数提取 (Eq. 4.20, 4.54)：双轴/扭转实验是唯一能同时给出 \(W_1, W_2\) 的实验。Ch 5 详细给出实验设计。
  2. Lagrange 乘子 \(p\) 的边界条件求解 (Eq. 4.15—4.17, 4.41—4.43)：先由 \(t_{33}=0\) 或径向平衡把 \(p\) 解出来，再求其他应力分量。
  3. 不可压缩验证 (Eq. 4.58, 4.30)：扭转时测外半径 \(r\) 是否 = 拉伸后半径 (Eq. 4.30)；圆管时只需测一个 \(R\) 即可推所有 \(r\)。实验上验证 \(J=1\) 的简便方法。
  4. 非齐次变形 + 复合变形 (Eq. 4.28—4.33)：把"先拉后扭"拆为 \(\mathbf{F}_1\) (拉) × \(\mathbf{F}_2\) (扭)，让两个简单问题分别有解。Ch 7血管"轴向预拉伸 + 环向拉伸"也用这个套路。
  5. Mooney-Rivlin vs Ogden vs Fung 指数型比较 (4.1.5, 4.2.1末段)：Humphrey警告"更多参数 → 描述能力更好，但预测能力未必好"。奥卡姆剃刀应基于"参数数量 vs 数据点"权衡。

§5 关键结论

  1. 橡胶和软组织行为高度相似——4.1.1节明确：非线性、有限应变、近似不可压缩、各向同性 (or 横观各向同性) → 同一套本构理论。
  2. 本构的核心 = 1个标量 \(W(\mathbf{C})\) (Ch 3结论) + 不可压缩 → 1个 Lagrange 乘子 \(p\)。Ch 4把 \(W\) 的具体形式列出 6 大类，每类都有适用场景。
  3. 双轴实验 = 唯一能唯一确定 \(W_1, W_2\) 的实验。单轴只能确定 \(W_1 + \lambda^{-1} W_2\) 这种组合。这是 Ch 5 实验设计的根本依据。
  4. 不可压缩 → 1个未知半径解所有 (Eq. 4.58, 4.30)。实验上验证 \(J=1\) 的简便测试。
  5. 复合变形 (F = F₂·F₁) 是处理"非均匀 + 复杂加载"的标准范式。Ch 7血管的轴向+环向+剪切都用这个思路。
  6. 本构不存在"唯一正确"形式：neo-Hookean、Mooney-Rivlin、Ogden、Fung指数、Valanis-Landel、Rivlin-Sawyers 都是"候选"。DEICE法强调"在适用条件下验证预测能力"，不是"哪条形式更真"。
  7. 复合材料和不可压缩"几乎完全"假设的合理性：橡胶 \(J - 1 \sim 10^{-4}\)，血管 \(J - 1 \sim 10^{-3}\)——\(J=1\) 是合理简化；非 \(J=1\) 时引入 \(W(I_C, II_C, III_C)\)。

§6 挑战和开放性问题

  1. 微观机制 (Treloar 纯熵型) 与现象学 (Fung, Mooney) 的鸿沟：Eq. 4.4 neo-Hookean 在 \(\lambda \leq 1.3\) 表现良好，但 \(c_1 = nkT\) 在 \(T\) 变化时如何直接测量 \(n\)？这是"理论参数"和"实验参数"之间的待解问题。
  2. Valanis-Landel可分性 (Eq. 4.9) 是否普遍成立？ Jones-Treloar 1975 验证橡胶 \(\lambda \in [0.189, 2.62]\) 可分，但血管是否可分至今有争议 (Humphrey 1995 评论过)。
  3. \(I_C, II_C\) 在小应变下数值相近 (exercise 4.7) → 双轴实验提取的 \(W_1, W_2\) 对小应变数据噪声极其敏感。Valanis-Landel (Eq. 4.9) 正是为避开此问题提出。
  4. Criscione et al. 2000 指出 \(W(I_C, II_C)\) 形式有概念问题——建议用"体积变化 + 偏量"分离 (Criscione-Humphrey 2000 框架)。Humphrey 2002 写书时此工作刚出，未在Ch 4 深入展开——但这是后续Ch 7、10 实际选择的依据。
  5. 厚壁管的"管壁 + 残余应力"耦合：4.2.3节给的"管 + 拉 + 扭"是"无残余应力"理想解。Ch 7会用 Chuong-Fung 1983 的"开环"模型扩展，加 \(\sigma_r^{res}(R)\)。Ch 4 是必要但非充分准备。
  6. 复合变形顺序依赖：F = F₂·F₁ ≠ F₁·F₂ (矩阵乘法不交换)，对粘弹性或塑性材料，应变历史比当前应变状态更重要。本章的弹性假设在这一步失效。

§7 个人反思与批判性分析

Humphrey的橡胶-软组织类比是本书最聪明的教学trick：读者先学橡胶的6种经典本构 + 5种通用解，再看到Ch 7的血管"和橡胶几乎一模一样"时就不会再被新物理所吓到。这种"先固化、再泛化"的方法适合工程背景读者。

值得自己复现的推导： - 4.20 (双轴 → \(W_1, W_2\))：解一个 \(2\times 2\) 线性方程，是本构建模的"母公式"。 - 4.43 (圆杆扭转 \(p(r)\))：径向平衡 + \(\sigma_{rr}(r_2) = 0\) 边界条件 → 一阶ODE积分。Ch 7血管 \(p(R)\) 的完全类似推导。 - 4.54—4.55 (扭转 → \(W_1, W_2\))：用Leibnitz法则从 \(M(y), L(y)\) 反解 \(W_1, W_2\)。这是 Ch 5 实验设计"扭转实验"的核心公式。 - 4.58 (厚壁管 \(\Lambda\) 测一个 \(r\) 知所有 \(r\))：实验验证不可压缩的"最小信息需求"。

与作者的潜在对话： - 4.1.4节末段您说"Ogden \(N=4\) 拟合 Treloar 数据比 \(N=3\) 略好，但可能只是因为参数多了两个"——是否建议读者用AIC/BIC等模型选择准则？DEICE的"Evaluate"步骤是否应显式包含奥卡姆信息准则？ - 4.1.3节 列出6种本构，但没说"在 \(\lambda\) 极大 (rupture) 区域哪种最稳定"——对动脉瘤/血管破裂预测至关重要。 - 4.2.2节 \(\pm(M/y)\) 对称性验证各向同性——但"对称"是数学约束，"各向同性"是物理假设。中间缺一个"微观结构 → 各向同性 → \(W(I_C, II_C)\)"的论证。这个论证您会留到 Ch 7 吗？

总结性批评： - 4.1.5节(预测能力)只有3句——是全书"DEICE法"五步中最弱的一节。如果读者做Ch 8 血管瘤风险评估，会发现"拟合好"和"预测好"差距可以很大。需要更多具体例子。 - 4.2.3节厚壁管的"扭转"部分给得很简略 (Eq. 4.56—4.58之后只一页)。Ch 7血管会重做，但Ch 4 这一节的"管 + 拉 + 扭"是 Ch 9 重塑建模的几何基础——应有更详尽的应力场。 - 4.3节可压缩比 4.1 节短得多——但心肌的"可压缩性"是Ch 10 的核心议题 (心肌在收缩时体积变化 ~ 5-8%)。建议在4.3节给至少一个"心肌可压缩"型本构示例。 - 4.4节膜超弹性只给"Ch 7会用"，但Ch 7的血管壁是 3层结构 (intima, media, adventitia)，每层有不同本构。4.4节应至少给出"层合板 + 不同本构叠加"的方程框架。 - 未讨论的"残余应力与本构的耦合"：Ch 4 假设参考构形 \(\kappa_0\) 是"无应力"——但血管的 \(\kappa_0\) 通常有残余应力 (Chuong-Fung 1983)。这是 Ch 4 与 Ch 7 之间的"隐含不连续"，4.4节末段应明确指出。

§8 重要参考文献

  • [X1] Ogden RW (1984) Non-Linear Elastic Deformations. Wiley. — 不可压缩有限弹性的"圣典"，Ch 4 主要追随。
  • [X2] Treloar LRG (1975) The Physics of Rubber Elasticity. Clarendon. — 纯熵统计力学的源头，neo-Hookean (Eq. 4.4) 的来源。
  • [X3] Rivlin RS, Saunders DW (1951) Large elastic deformations of isotropic materials. VII. Philos Trans R Soc Lond A243:251-288. — 双轴"通用解"提取 \(W_1, W_2\) 的来源 (Eq. 4.20)。
  • [X4] Rivlin RS, Sawyers KN (1976) Some constitutive equations for non-linear fluids. J Non-Newtonian Fluid Mech 1:179-202. — 25年后给出 \(f(II_C)\) 的具体形式 (Eq. 4.7)。
  • [X5] Valanis KC, Landel RF (1967) The strain-energy function of a hyperelastic material in terms of the extension ratios. J Appl Phys 38:2997-3002. — 可分型 (Eq. 4.9)。
  • [X6] Carmichael AJ, Holdaway GW (1961) The influence of specimen geometry on the biaxial extension of rubber. J Appl Phys 12:309-314. — Valanis-Landel 的早期先驱。
  • [X7] Kearsley EA, Zappas LM (1980) Some experiments on the torsion of rubber. Trans Soc Rheol 24:483-500. — 扭转实验提取 \(W\)。
  • [X8] Penn RW, Kearsley EA (1976) The scaling law for finite torsion of slightly compressible rubber. Int J Solids Struct 12:21-30. — 扭转反演 \(W_1, W_2\) 的具体公式。
  • [X9] Blatz PJ, Ko WL (1962) Application of finite elasticity to the deformation of rubbery materials. Trans Soc Rheol 6:223-251. — 泡沫橡胶可压缩本构 (Ch 4.3)。
  • [X10] Mullins L, Tobin NR (1957) Stress softening in rubber vulcanizates. J Appl Polym Sci 9:2993-3000. — Mullins 效应 (Ch 4.1.1) 的来源。
  • [X11] Fung YC (1967) Elasticity of soft tissues in simple elongation. Am J Physiol 28:1532-1544. — 软组织 Fung 指数型本构的来源 (Eq. 4.8)。
  • [X12] Hart-Smith LJ, Crisp JDC (1967) Large elastic deformations of thin rubber membranes. Int J Eng Sci 5:1-10. — 指数型本构先于Fung的软组织版本。
  • [X13] Humphrey JD, Yin FCP (1986) Fiber-induced material behavior in composites. Mech Res Comm 13:277-283. — Ch 4.1 准备 Ch 7 纤维方向本构的"过渡"。
  • [X14] Criscione JC, Humphrey JD, Douglas AS, Hunter WC (2000) An invariant basis for natural strain which yields orthogonal stress response tensors in isotropic hyperelastic materials. J Mech Phys Solids 48:2445-2465. — 4.1.5节末段提及的"Criscione 不变量"工作。
  • [X15] Jones DF, Treloar LRG (1975) The properties of rubber in pure homogeneous strain. J Phys D 8:1285-1304. — Valanis-Landel 可分性实验验证。