跳转至

第7章：主动血管壁（The Vascular Wall, an Active Entity）

1. 作者

T. Christian Gasser (KTH Royal Institute of Technology) 独立撰写本章。Ch 7 是 Gasser 教科书的最核心一章, 也是 Gasser 长期科研的主战场 (Gasser 在 2000-2020 年代发表了大量 G&R 律和多尺度血管壁建模的论文)。

章节范围: 7.1 Introduction → 7.6 Modeling Arteriogenesis (PDF p367-p421, 55 页)

核心定位: Ch 7 是Ch 5 (被动血管壁) 的延伸。Gasser 在 7.1 开篇就直言: "血管是动态结构, 通过主动收缩 (vasoreactivity) 和组织重塑 (arteriogenesis) 适应环境变化"。

两套时间尺度: - Vasoreactivity (血管反应性): 秒级, SMC 主动收缩/舒张 - Arteriogenesis (动脉生成): 天-周级, 胶原/弹性蛋白合成与降解

Gasser 的核心立场: 完整的血管生物力学必须"被动 (Ch 5) + 主动 (Ch 7) + 流体 (Ch 6)"三者耦合, 任何一项缺失都是简化。本人课题的 G&R 模拟直接站在 Ch 7 的肩膀上。

2. 内容概述

本章系统建立主动血管壁的生物力学框架: - 7.1 介绍 4 类机械传感器 (核膜, 细胞膜, 微域, 细胞支架) + 3 种主动响应 (contractility, remodeling, growth) - 7.2 Vasoreactivity: SMC 结构, Ca²⁺-MLCK 通路, NO-cGMP 通路, 肌源性反应, latch bridges - 7.3 Arteriogenesis: EC-SMC 通信, WSS 谱与 NF-κB, 胶原/弹性蛋白合成通路 (intracellular + extracellular) - 7.4 Angiogenesis (新血管生成) — Gasser 仅简述 - 7.5 Vasoreactivity 建模: Hill 三参数, 现象学, 钙离子浓度基础 - 7.6 Arteriogenesis 建模: 开放系统热力学, 运动学 G&R, 连续 turnover-based G&R, constrained mixture, 例 (薄壁管 G&R)

Gasser 的核心方法学: 开放系统热力学 描述组织 mass turnover; 运动学 multiplicative decomposition (\(\mathbf{F}_{tot} = \mathbf{F} \cdot \mathbf{G}\)) 描述生长; constrained mixture 描述多组分演化。这是 Gasser 在 G&R 律上的标志性工作 (与 Humphrey, Taber, Rodriguez, Cyron 并列)。

3. 核心方程与概念

3.1 主动血管壁的 3 类响应 (Sect. 7.1)

  1. Contractility (收缩性): SMC 肌张力改变 (秒-分钟)
  2. Remodeling (重塑): 宏观力学性质改变 (周-月)
  3. Growth/Atrophy (生长/萎缩): 组织质量改变 (月-年)

4 类机械传感器: - 核膜分子 (nuclear pores) - 细胞膜分子 (ion channels, receptors, glycocalyx) - 微域 (primary cilia, caveolae) - 细胞支架 (cytoskeleton, lipid bilayer)

3 类稳态目标 (homeostatic targets): - WSS: 引用 [74, 236] - 周向壁应力: 引用 [363, 598] — uniform stress hypothesis (Ch 5 引用 [183]) - 轴向预拉伸: 引用 [215, 283]

3.2 Vasoreactivity (Sect. 7.2)

3.2.1 SMC 基础

  • 形态: 纺锤形, 2-10 μm 直径 × 50-400 μm 长; 沿 MLU 周向排列, 径向倾斜 ~20° (引用 [180, 400])
  • 功能: 收缩 + ECM 合成 (与表型相关)
  • 收缩装置: 肌动蛋白 (actin) + 肌球蛋白 (myosin) + 中间丝 (vimentin, desmin) + 致密体 (dense bodies) + 致密斑 (focal adhesions)
  • 收缩机制: antiparallel cross-bridge (肌动-肌球蛋白反向桥接), 步距 10-12 nm

3.2.2 Ca²⁺-MLCK 通路 (Fig. 7.3)

  1. 静止: 自由 Ca²⁺ ~80-140 nmol/L
  2. 激活 (收缩): 500-700 nmol/L
  3. Ca²⁺ + 钙调蛋白 (calmodulin) → Ca²⁺-CaM 复合物
  4. 激活 MLCK (Myosin Light Chain Kinase)
  5. 磷酸化 MLC (Myosin Light Chain)
  6. 允许 actin-myosin cross-bridge 形成 → 收缩
  7. MLCP (Myosin Light Chain Phosphatase) 去磷酸化 → 舒张

3.2.3 NO-cGMP 通路 (Fig. 7.3d)

  1. EC 释放 NO (eNOS 催化)
  2. NO 扩散进入 SMC, 激活 soluble guanylate cyclase
  3. 增 cGMP → 激活 PKG
  4. PKG 激活 MLCP + 增 SERCA (Ca²⁺ 回收到肌浆网) + 增 NCX/PMCA (Ca²⁺ 外流) → 舒张
  5. Rho kinase 通路 (Fig. 7.3c) 抑制 MLCP → calcium sensitization (收缩敏感性增)

3.2.4 Latch Bridges

SMC 特色: 即使 Ca²⁺ 降低, 仍能维持长时间收缩 (与骨骼肌不同) - 机制: latch bridges (去磷酸化但仍附着的桥) - 意义: 最小 ATP 水解下保持力, 长时间"低能耗收缩"

3.2.5 SMC 对血管壁力学的影响

  • SMC 收缩可减小血管直径 20-50% (引用 [25, 104, 130, 385])
  • 收缩在体内直径 (homeostatic length) 处产生最大主动应力 (引用 [387])
  • muscular 动脉 收缩显著改变 p-d 性质; elastic 动脉 收缩改变较小 (Fig. 7.5, 引用 [105])
  • 冠脉独特: muscular 动脉但收缩效应很小 (Fig. 7.5)
  • Isometric 收缩 (等长) 产生最大力; Isotonic 收缩 (等张) 做外功

3.3 Arteriogenesis (Sect. 7.3)

3.3.1 WSS 谱与 EC 反应 (Fig. 7.6, 引用 [150])

3 类 WSS 模式: - Static WSS: 常量 - Pulsatile WSS: 方向不变 + 时变 (保护性) - Oscillatory WSS: 方向变化 (前/后流) — 致动脉粥样硬化

NF-κB 表达 反映 WSS 致炎性: 振荡 WSS → 高 NF-κB → 内皮激活 → 炎症 → 动脉粥样硬化

3.3.2 胶原合成通路 (Fig. 7.7, Sect. 7.3.3)

Intracellular (细胞内): 1. 核内: DNA → mRNA (转录) 2. 胞质: 核糖体翻译 pre-procollagen 3. 内质网: post-translational modification (去除 N-端信号肽, 羟基化, 糖基化) → procollagen α 链 4. 三条 α 链 → procollagen 分子 5. 高尔基体: 装入分泌囊泡 → 胞吐到胞外

Extracellular (细胞外): 6. 端肽酶 (collagen peptidases) 切去两端 → tropocollagen 7. 赖氨酰氧化酶 (lysyl oxidase, Cu 依赖) → 共价交联 → 胶原原纤维 (fibril) 8. 许多 fibrils → 胶原纤维 (fiber, 5-30 nm) → 血管壁主承力单元

3.3.3 弹性蛋白合成 (Sect. 7.3.4)

  • 仅在生长/伤口愈合/妊娠/重塑中显著
  • 出生后弹性蛋白合成基本停止
  • 同样需要 lysyl oxidase + desmosine/isodesmosine 交联
  • 半衰期数十年 (引用 [7]) — 一旦降解, 难修复

3.4 Vasoreactivity 建模 (Sect. 7.5)

3.4.1 Hill 三参数 (Eq. 7.1, 骨骼肌)

\[ P = \frac{P_0(1 - v/v_0)}{1 - (P_0/a)(v/v_0)} \]
  • \(P_0\): 最大等长张力, \(v_0\): 最大收缩速度, \(a\): 形状参数
  • 应用: 骨骼肌广泛, SMC 中很少用 (Gasser 引用 1: rare)

3.4.2 现象学 (Sect. 7.5.2, Table 7.1)

总应力 = 被动 (ECM) + 主动 (SMC):

\[ P_{total} = P_{ECM}(\lambda_\theta) + P_{SMC}(\lambda_\theta) \]

  • 被动: \(P_{ECM} = c_{ECM}(\lambda_\theta - 1)\exp[(\lambda_\theta - 1)^2]\)
  • 主动: \(P_{SMC} = c_{SMC}\exp[-(\lambda_\theta - \lambda_m)^2/(2\lambda_d^2)]\) — 高斯激活曲线, \(\lambda_m\) 峰值位置, \(\lambda_d\) 分散度
  • 肌张力 \(0 \leq \alpha \leq 1\): \(P_{SMC} = \alpha P_{SMC}^\star\), \(\alpha\) 是激活水平

Example 7.1 (双轴血管壁片): 计算"主动 + 被动"工作, 等长 vs 等张收缩

3.4.3 Calcium-based (Sect. 7.5.6, Ch 3 + Ch 7 整合)

  • Free Ca²⁺ 浓度 \(\beta\) 是主动应力的源头
  • MLC 磷酸化水平决定 cross-bridge 状态
  • Hai-Murphy 4-state 模型 (Eq. 7.32-7.34):
  • \(\alpha_1\): cycling cross-bridges
  • \(\alpha_2\): latch bridges (未磷酸化但附着)
  • \(\alpha_3\): detached cross-bridges
  • \(\alpha_4\): phosphorylated detached myosin
  • 主动应力:
    \[ P_{SMC} = k_{cb}\rho(\lambda_a)(\alpha_1 + \alpha_2)(\lambda_a - \lambda_{a,0}) + k_{lb}(\alpha_2/\lambda_a^{lb}) \]
  • \(\rho(\lambda_a)\): actin-myosin 重叠的 PDF (高斯, Eq. 7.35)

3.5 Arteriogenesis 建模 (Sect. 7.6)

3.5.1 开放系统热力学 (Sect. 7.6.1)

Mass balance (Eq. 7.41): \(\partial\rho/\partial t + \rho \cdot \text{div}\,\mathbf{v} = \rho\varsigma_v\), \(\varsigma_v\) = 体积源

Slow growth assumption: \(t\) (秒, 弹性) \(\ll \tau\) (天, 生长), 两者解耦

Momentum (Eq. 7.42): 与 Ch 3 相同, 生长不影响

First Law (Eq. 7.43): \(\partial u/\partial t + \text{div}\,\mathbf{q}_h - r_h - \boldsymbol{\sigma}:\mathbf{d} = \varsigma_h\), 包含非兼容热源 \(\varsigma_h\)

Second Law (Eq. 7.44, Clausius-Duhem):

\[ \gamma\theta = -\frac{\partial\psi}{\partial t} - s\frac{\partial\theta}{\partial t} + \boldsymbol{\sigma}:\mathbf{d} - \frac{\mathbf{q}_h\cdot\nabla\theta}{\theta} + \varsigma_s \geq 0 \]
- \(\varsigma_s\): 非兼容熵源, 来自组织生长 (合成新纤维, undulation 不同)

3.5.2 Kinematics-based Growth (Sect. 7.6.2)

Multiplicative decomposition (Eq. 7.45): \(\mathbf{F}_{tot} = \mathbf{F} \cdot \mathbf{G}\) - \(\mathbf{G}\): 生长变形 (从 initial reference \(\widetilde{\mathcal{B}}_0\) 到 intermediate reference \(\mathcal{B}_0\)) - \(\mathbf{F}\): 弹性变形 (从 \(\mathcal{B}_0\) 到 spatial \(\mathcal{B}\)) - 关键: intermediate \(\mathcal{B}_0\) 是非兼容的 (不同材料点有自己的 \(\mathcal{B}_0\)), 即使在 \(\mathbf{t} = \mathbf{0}, \mathbf{b}_f = \mathbf{0}\) 也有残余应力

Spatial distribution of volume growth (Sect. 7.6.3): - Constant-volume growth: \(\dot{\mathbf{G}} = \mathbf{0}\), 只改密度 - Constant-density growth (Eq. 7.50): \(\text{tr}(\dot{\mathbf{G}}\mathbf{G}^{-1}) = \varsigma_v\) - 各向同性生长: \(\mathbf{G} = \alpha\mathbf{I}\), \(\dot{\alpha} = \alpha\varsigma_v/3\) - 各向异性生长 (Eq. 7.50 推导): 沿解剖学方向, \(\dot{G}_i/G_i\) 沿主方向

Homeostatic target (Sect. 7.6.3.3): 常用 \(\varsigma_v = \zeta(\lambda_\theta - \overline{\lambda_\theta})\) (Eq. 7.54), 与 homeostatic stretch 的偏离驱动生长

3.5.3 Continuous Turnover-based Growth (Sect. 7.6.5-7.6.6)

生存函数 \(S(\tau - \tau_0) \in [0, 1]\): 在 \(\tau_0\) 沉积的组织在 \(\tau\) 仍存留的比例 - \(S(0) = 1\), \(dS/dx < 0\) - 典型: \(S(\tau) = \exp(-\tau/\zeta)\) (指数衰减, \(\zeta\) 是半衰期)

参考质量密度 (Eq. 7.58):

\[ \rho_0(\tau) = \rho_0(0)S(\tau) + \int_0^\tau \varsigma_M^+(\xi) S(\tau - \xi) d\xi \]
- 与 Ch 3 黏弹性卷积形式相同 (Gasser 引用 3.5.4.2)

应变速率 (Eq. 7.59): 旧组织 + 新沉积组织的应变能加权和

主 Cauchy 应力 (Eq. 7.62):

\[ \sigma_i(\tau) = \frac{\rho_0(0)S(\tau)}{\rho_0(\tau)} \lambda_i \frac{\partial\Psi(\mathbf{F}(\tau, 0))}{\partial\lambda_i} + \int_0^\tau \frac{\varsigma_M^+ S(\tau - \xi)}{\rho_0(\xi)} \lambda_i(\tau, \xi) \frac{\partial\Psi(\mathbf{F}(\tau, \xi))}{\partial\lambda_i(\tau, \xi)} d\xi - \kappa \]

单轴拉伸 (Sect. 7.6.6, Eq. 7.61-7.64): \(\lambda_{tot} = 1 + k\tau\), 给出解析 \(\rho_0(\tau)\) (Eq. 7.61) - 5 周后达到稳态应力 - 物理图像: 旧组织降解 + 新组织沉积 → 应力平衡

Constrained Mixture Model (Sect. 7.6.5): Gasser 引用 [21, 107, 214, 272, 273, 314, 329, 330, 331, 351, 383] - 多组分 (胶原, 弹性蛋白, SMC) 各自独立跟踪 - 沿 mixture 框架, 每组分的 \(\mathbf{F}\) 可能不同 (constrained)

3.5.4 厚壁管 G&R (Sect. 7.6.4, Table 7.6, Fig. 7.23)

  • 初始 \(\widetilde{R}_i = 3.0\) mm, \(\widetilde{R}_o = 4.0\) mm
  • \(p_i = 13.33\) kPa, neoHookean \(c = 120\) kPa
  • \(\zeta = 0.25\) weeks⁻¹ (生长时间常数)
  • m = 30 层空间离散
  • 结果 (Fig. 7.23a): 40 周后稳态, 内侧 \(\alpha = 1.03046\) (生长), 外侧 \(\alpha = 0.968924\) (萎缩), 稳态 homeostatic stretch \(\overline{\lambda_\theta} = 1.10072\)
  • Fig. 7.23b: 预设 \(\overline{\lambda_\theta} = const\) 给非物理行为 (持续膨胀或收缩) → 启示: homeostatic target 应该是非局部的 (整个壁平均), 不是局部 const

4. 关键结论

  • [K1] 血管是动态主动结构: 主动行为分 3 层 (contractility, remodeling, growth), 时间尺度跨 10⁸ (秒到年)
  • [K2] Vasoreactivity 与 Arteriogenesis 时间尺度差 6+ 数量级: 慢生长假设 (slow growth) 使两者解耦, 可分别建模再耦合
  • [K3] SMC 收缩是 NO-cGMP 通路的直接输出: 任何破坏 NO 平衡 (高血压, 糖尿病, 老化) → SMC 失活 → 血管壁 G&R 失衡
  • [K4] Latch bridges 是 SMC 长时间低能耗收缩的机制; 骨骼肌没有 → 长期 vasoconstriction 不消耗 ATP
  • [K5] collagen 半衰期 60-70 天 vs elastin 数十年的巨大差异 → 弹性蛋白一旦破坏, 难修复; 老化血管以胶原代偿, 失去弹性
  • [K6] vasoreactivity 建模有 3 套体系: Hill (骨骼肌专用, SMC 少用), Phenomenological (Eq. 7.3, 简单), Calcium-based (Hai-Murphy 4-state, 物理完整)
  • [K7] arteriogenesis 建模有 2 大体系: Kinematics-based (F·G 分解, 显式 homeostatic target), Turnover-based (survival function + convolution, 与 Ch 3 黏弹性同构)
  • [K8] Multiplicative decomposition \(\mathbf{F}_{tot} = \mathbf{F} \cdot \mathbf{G}\) 源自 elasto-plasticity (Gasser 引用 [313]), Gasser 把它扩展到 G&R
  • [K9] Open-system 热力学 是 G&R 律的热力学基础: Clausius-Duhem 含非兼容熵源 \(\varsigma_s\) (Eq. 7.44), 反映新合成的纤维与旧纤维有不同结构
  • [K10] Intermediate reference \(\mathcal{B}_0\) 是非兼容的 (Sect. 7.6.2): 每个材料点有自己的 reference, 即使无外载也有残余应力 → 解释 Ch 5 的残余应力
  • [K11] Homeostatic target 应该是非局部 (整个壁平均) 而不是局部 const (Fig. 7.23 vs 7.23b) — Gasser 给的关键工程经验
  • [K12] Constrained Mixture Model 是 Gasser 引用 [107, 329, 330] 提出的多组分 G&R 框架, 适合胶原/弹性蛋白/SMC 各自有不同时间常数
  • [K13] Survival function \(S(\tau - \tau_0) = \exp(-(\tau-\tau_0)/\zeta)\) 是 G&R 的核心, 给出解析的 \(\rho_0(\tau)\) (Eq. 7.61)
  • [K14] Turnover-based 与 Kinematics-based G&R 的对偶: 前者显式描述"新组织沉积 + 旧组织降解", 后者显式描述"参考构型演化"; 两者物理等价, 计算效率不同
  • [K15] WSS 振荡性 → NF-κB 表达 → 动脉粥样硬化: 这是机械-化学-病理链条的经典例 (引用 [150])
  • [K16] WSS-driven 血管重塑 (引用 [295, 296, 549]): 流量变化 → 血管直径比例变化 → 恢复 WSS → 稳态
  • [K17] Angiogenesis 主要是 hypoxia 驱动, 而 arteriogenesis 是 WSS/应力驱动 — 两个不同的血管生物力学机制
  • [K18] SMC 表型转换 (Sect. 7.2.1): 收缩型 (静息) ↔ 合成型 (增殖/迁移) ↔ 降解型 (类巨噬细胞) — 可逆, 力学 + 化学触发
  • [K19] muscular vs elastic 动脉的 G&R 律可能不同: Gasser 没明确给不同本构, 但实际工程应该区分
  • [K20] SMC 收缩对刚度的影响是双相的 (Fig. 7.4b, 引用 [622]): 小直径时增刚度, 大直径时减刚度 → 简单的"增刚度"假设错误
  • [K21] Angiogenesis (新血管生成) 是 Gasser 没展开的主题, 仅在 7.4 简述; 完整展开需要 Ch 8+ 但本书无

5. 挑战和开放性问题

  • [C1] Vasoreactivity 建模的 IMA 决定: Hill 简单但物理粗糙, Hai-Murphy 物理完整但 4-state + 9 参数; 没有标准
  • [C2] Calcium-based 模型的 4-state 多参数标定: \(k_{cb}, k_{lb}, \lambda_a, \lambda_{a,0}, \sigma\) 等至少 5+ 参数, 实验标定困难 (Fig. 7.6 数据是 EC 数据, 不是 SMC)
  • [C3] SMC 表型转换的 G&R: Gasser 给合成型 (产生 ECM) 与收缩型, 但表型转换本身的数学模型缺位 — 应该是 Ch 8 主题, 但本书无
  • [C4] Arteriogenesis 的 homeostatic target 物理意义: \(\overline{\lambda_\theta} = const\) 给非物理 (Fig. 7.23b) — 正确的 homeostatic target 是什么? Gasser 没明说, 实战多用应力 (Taber-Humphrey) 或 WSS
  • [C5] Turnover-based G&R 的"高计算代价": 每时间步要积分全历史 (Eq. 7.58 卷积), 长期模拟 (年) 代价大; Gasser 引用 [107, 329, 330] 的homogenized constrained mixture 是出路
  • [C6] Constrained Mixture 的"组分"定义: Gasser 引用 [21, 273] 给出胶原, 弹性蛋白, SMC 三组分, 但每组分的"代表性 stretch" 怎么选? 文献中常用 affine transformation, 但实验验证少
  • [C7] Anisotropic growth 的方向选择: 7.6.3.2 给"沿解剖学方向", 但真实血管的解剖学方向本身就是 G&R 的结果, 形成循环依赖
  • [C8] Multiplicative decomposition 的非唯一性: \(\mathbf{F}_{tot} = \mathbf{F}\mathbf{G}\) 在 3D 有 9 个标量未知数, 但 Eq. 7.50 只有一个标量方程 — 必须加 6 个运动学假设 (Gasser 给 isotropic + anisotropic 两种, 但选择标准缺位)
  • [C9] Vasoreactivity + Arteriogenesis 耦合: Gasser 把两者分章, 但实际血管两者强耦合 (SMC 收缩改变 WSS → arteriogenesis → 又影响 SMC) — 统一本构 是开放问题
  • [C10] 3 类传感器 (nuclear, membrane, microdomain, cytoskeleton) 的统一数学描述: Gasser 在 7.1 列出 4 类, 但没给出统一数学模型; 现代 mechanobiology 文献 (Discher, Engler) 给具体实现, 但 Gasser 教材未涉及
  • [C11] NO 信号的扩散-反应动力学: NO 半衰期 ~秒, 扩散距离 ~100 μm; 真实血管壁厚 ~1 mm — NO 不会均匀分布, 但 Gasser 假设均匀 (Eq. 7.43 标量化)
  • [C12] 血管壁的炎症-力学耦合: NF-κB 是炎症主调节因子 (引用 [150]), 与力学反馈强耦合, 但 Gasser 把它归到"血管内皮感知", 没展开力学-炎症双向耦合
  • [C13] Angiogenesis 的机械因素: 7.4 仅一句话 "beyond the scope", 实际 angiogenesis 与 arteriogenesis 关系密切 (Gasser 引用 [505] 指向综述)
  • [C14] G&R 律的临床转化: 临床上仍用"直径"判断 AAA, Gasser 引用 RRED 框架 (Ch 1) 还没 RCT 验证
  • [C15] 多场耦合的多尺度: 分子 (NO, Ca²⁺) → 细胞 (SMC) → 组织 (MLU) → 器官 (主动脉), 跨 6 个数量级; G&R 模型集中在组织级, 细胞级用 Hai-Murphy, 分子级没给
  • [C16] 机器学习 G&R 律: 2022+ 文献中 Neural ODE + PINN 已开始替代经典 G&R, Gasser 教材未反映

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

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

本章是本人研究核心: - Hai-Murphy 4-state 模型 (Sect. 7.5.6): 本人 SMC 收缩本构直接使用 - Survival function + convolution (Eq. 7.58): 本人 turnover-based G&R 的核心 - Constrained Mixture Model (Sect. 7.6.5): 本人多组分 (胶原/弹性蛋白/SMC) G&R 的理论框架 - Homeostatic target (Sect. 7.6.3.3): 本人 G&R 模拟的核心驱动 - Open-system 热力学 (Sect. 7.6.1): 本人热力学一致的 G&R 律 基础 - Multiplicative decomposition (Eq. 7.45): 本人FEniCS 实现的 G&R 框架 - Uniform stress hypothesis (引用 [183]): 本人研究"为什么 homeostatic target 是应力" 的核心

本人研究的 IMA = "长期 G&R (10+ 年)", 本人直接站在 Ch 7 肩膀上。

6.2 Gasser 框架的优势

  • 多尺度统一: Ch 5 (被动) + Ch 6 (流体) + Ch 7 (主动) 在同一种数学语言 (连续介质力学 + 热力学) 下统一
  • 工程友好: 7.6.4 的厚壁管 G&R 算法 (Table 7.6) 是FEniCS 实现的范本
  • 物理完整: 开放系统热力学 + 卷积 survival function + 4-state Hai-Murphy → 几乎所有主流 G&R 律都涵盖
  • 临床对应: Ch 7 与 Ch 1 BRRA / Ch 5 AAA 案例直接对应, 工程落地路径清晰

6.3 Gasser 框架的局限性

  • 多组分 constrained mixture 的"组分定义"模糊: 胶原/弹性蛋白/SMC 是粗粒化, 真实 ECM 包含多种亚型 (Type I/III/IV/V/VI), 每种亚型不同动力学 → Gasser 模型不能区分
  • 3 类传感器 (4 类实为) 的数学描述缺位: Gasser 列出 4 类传感器, 但没给出统一 ODE/PDE 模型; 现代文献 (Humphrey, 2020+) 给具体实现
  • 力学-化学-炎症-免疫-基因表达的多场耦合: Gasser 给出力学 + 部分化学 (NO, Ca²⁺), 但炎症 (NF-κB), 免疫 (macrophages), 基因表达 (miRNA) 都没纳入
  • 细胞迁移: G&R 律通常假设无迁移, 但 SMC 表型转换 + 迁移是动脉粥样硬化的关键 (Gasser 引用 [52] 提到迁移但没建模)
  • Angiogenesis 没展开: 7.4 仅 1 段, 完整 angiogenesis 涉及分叉 (tip cell, stalk cell), Gasser 引用 [505] 但没给方程
  • 空间异质性: 大动脉 G&R 是局部不同步 (粥样硬化灶, ILT 区域), Gasser 给出单层均匀模型
  • 个体差异: Gasser 给"群体平均参数", 个体化 G&R 参数难以获得 → 临床转化困难
  • 2022+ 进展: PINN G&R, Neural ODE G&R, Multi-fidelity G&R 已在 2022+ 涌现; Gasser 教材未涉及

6.4 一个观察: "Homeostatic target" 的物理意义

Gasser 在 7.6.3.3 给"homeostatic target" 概念, 但没明说应该是什么物理量。Fig. 7.23 给出关键工程经验:

  • a) 整个壁平均的 \(\overline{\lambda_\theta}\) → 稳态, 物理合理
  • b) 局部 const \(\overline{\lambda_\theta}\) → 非物理 (持续膨胀或收缩)

这是关键工程经验: homeostatic target 应该是全局稳态 (如壁平均应力), 不是局部预设常数。在本人研究中, 本人用"壁平均 Cauchy 应力" 作为 homeostatic target, 与 Gasser 经验一致。

6.5 一个批判: 开放系统热力学的"非兼容"概念模糊

Gasser 在 7.6.1.3-7.6.1.4 引入"非兼容热源 \(\varsigma_h\)"和"非兼容熵源 \(\varsigma_s\)" (Eq. 7.43-7.44)。这个"非兼容"是 Gasser 教科书自创概念 (与"质量交换是否改变系统能量"有关), 但:

  • 没有明确定义: 什么算"非兼容" vs "兼容"?
  • 工程实现困难: \(\varsigma_s\) 怎么从实验测量?
  • 与已有 G&R 律的对应: Humphrey 的"stress-mediated growth"和 Taber 的"kinematics-based growth"都可以用"非兼容"统一, 但统一的代价是引入了额外的标量参数

批判: Gasser 没告诉读者"非兼容是统一工具, 不是物理本质", 容易让初学者过度形式化。

6.6 一个意外发现: Turnover-based G&R 与 Ch 3 黏弹性同构

Gasser 在 7.6.5 写出 Eq. 7.58 的卷积形式, 然后明说 "this involves the solution of a time convolution similar to linear viscoelasticity"。这是 Gasser 教科书一个"暗线":

  • Ch 3.5.4 黏弹性: Prony 级数 + 卷积
  • Ch 5.5.10 黏超弹: Standard Solid + 不可压
  • Ch 7.6.5 连续 turnover: 生存函数 \(S(\tau) = \exp(-\tau/\zeta)\) + 卷积

本质: 这三个"不同物理" (黏弹, 黏超弹, G&R) 数学同构, 都是"历史依赖"问题。本人研究中复用 Ch 3 的黏弹性算法实现 Ch 7 的 G&R 律, 大幅降低代码量。Gasser 没明说"这是同构", 但 Eq. 7.58 与 3.121 的形式几乎相同, 给读者留下自学线索。

6.7 一个工程实践: G&R 律标定

Gasser 给出 G&R 律的理论框架 (Sect. 7.6), 但参数标定没详细给。实战中, 标定 G&R 参数比标定本构参数更困难:

  • 本构参数: 双轴实验可标定 (Ch 5.6)
  • G&R 参数: 需要长时程 (月-年) 动物实验 + 体内多次成像 → 临床难, 动物稀缺
  • Gasser 没展开: G&R 参数的不确定性传播 (Monte Carlo, Bayesian), 标定策略 (个体化 vs 群体)

批判: Ch 7 给理论深度, 但工程化 (参数标定 + 不确定性) 弱于 Ch 1 (BRRA 框架) 和 Ch 5 (本构标定)。

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

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

  • [X1] Fung Y.C. (1993) Biomechanics: Mechanical Properties of Living Tissues. Springer — 引用 183, uniform stress hypothesis 的奠基论文
  • [X2] Taber L.A., Humphrey J.D. (2001) Stress-modulated growth, residual stress, and vascular heterogeneity. J Biomech Eng 123:528-535 — 引用 558, stress-mediated G&R 经典
  • [X3] Rodriguez E.K., Hoger A., McCulloch A.D. (1994) Stress-dependent finite growth in soft elastic tissues. J Biomech 27:455-467 — 引用 473, kinematics-based G&R 经典
  • [X4] Himpel G., Kuhl E., Menzel A., Steinmann P. (2005) Computational modelling of isotropic multiplicative growth. Comput Model Eng Sci 8:119-134 — 引用 220, multiplicative growth 综述
  • [X5] Humphrey J.D. (2002) Cardiovascular Solid Mechanics: Cells, Tissues, and Organs. Springer — 引用 275, 心血管 G&R 教材
  • [X6] Holzapfel G.A., Gasser T.C., Ogden R.W. (2004) Comparison of a Gasser-Ogden-Holzapfel-based constitutive model with the classical model for arterial wall mechanics. In: Proceedings of the 2003 Summer Bioengineering Conference — 引用 261, GOH 在 G&R 中的应用
  • [X7] Baek S., Rajagopal K.R., Humphrey J.D. (2006) A theoretical model of enlarging intracranial fusiform aneurysms. J Biomech Eng 128:142-149 — 引用 21, constrained mixture model 在动脉瘤的应用
  • [X8] Cyron C.J., Humphrey J.D. (2017) Growth and remodeling of load-bearing biological soft tissues. Mech Res Commun 83:38-45 — 引用 107, homogenized constrained mixture
  • [X9] Latorre M., Humphrey J.D. (2018) A mechanobiologically equilibrated constrained mixture model for growth and remodeling of soft tissues. Z Angew Math Mech 98:2048-2068 — 引用 329-330, 平衡 constrained mixture
  • [X10] Valentin A., Humphrey J.D., Holzapfel G.A. (2013) A multi-layered computational model of coupled elastin degradation, vasoactive dysfunction, and collagenous stiffening in aortic aging. Ann Biomed Eng 41:957-970 — 引用 580, 老化血管 G&R
  • [X11] Wilson J.S., Baek S., Humphrey J.D. (2014) Importance of initial aortic properties on the evolving regional anisotropy, stresses and residual strains in the thoracic aorta. J Mech Behav Biomed Mater 29:521-532 — 引用 613, 主动脉 G&R
  • [X12] Ramachandra A.B., Humphrey J.D. (2019) Biomechanical characterization of murine aortic aneurysms. J Mech Behav Biomed Mater 99:1-10 — 引用 314, 老鼠动脉瘤 G&R
  • [X13] Rachev A., Stergiopulos N., Meister J.J. (1996) Theoretical study of dynamics of arterial wall remodeling in response to changes in blood pressure. J Biomech 29:635-642 — 引用 443, 应力反馈 G&R
  • [X14] Gleason R.L., Humphrey J.D. (2004) A 2D constrained mixture model for arterial adaption to sustained changes in blood flow. Biomech Model Mechanobiol 3:73-86 — 引用 214, 流量反馈 constrained mixture
  • [X15] Ambrosi D., Ateshian G.A., Arruda E.M., Cowin S.C., Dumais J., Goriely A., Holzapfel G.A., Humphrey J.D., Kemkemer R., Kuhl E., Olberding J.E., Taber L.A., Garikipati K. (2011) Perspectives on biological growth and remodeling. J Mech Phys Solids 59:863-883 — 引用 10, G&R 综述
  • [X16] Menzel A., Kuhl E. (2012) Frontiers in growth and remodeling. Mech Res Commun 42:1-14 — 引用 369, G&R 前沿综述
  • [X17] Hill A.V. (1938) The heat of shortening and the dynamic constants of muscle. Proc R Soc Lond B 126:136-195 — Hill 三参数模型的原始论文
  • [X18] Hai C.M., Murphy R.A. (1988) Cross-bridge phosphorylation and regulation of latch state in smooth muscle. Am J Physiol Cell Physiol 254:C99-C106 — 引用 196, Hai-Murphy 4-state 模型
  • [X19] Murphy R.A., Rembold C.M. (2005) The latch-bridge hypothesis of smooth muscle contraction. Can J Physiol Pharmacol 83:857-864 — 引用 384, latch bridges 综述
  • [X20] Rhodin J.A.G. (1980) Architecture of the vessel wall. In: Bohr D.F., Somlyo A.D., Sparks H.V. (eds) Handbook of Physiology, Section 2: The Cardiovascular System, Vol. II — 引用 400, SMC 几何 经典