跳转至

第2章：循环系统（The Circulatory System）

1. 作者

T. Christian Gasser (KTH Royal Institute of Technology 生物力学教授; 兼任 VASCOPS GmbH 联合创始人, A4clinics 软件开发者) 独立撰写本章; 章节 2.1-2.5 涵盖心血管系统解剖、生理、力学系统性质、宏观/微观循环集总参数建模。

章节范围: 2.1 Introduction → 2.5 Summary (PDF p57-p112, 56 页)

核心定位: 全书的"系统生理学"基础——为后续 Ch 3 连续介质力学、Ch 5 管道血管、Ch 6 血流动力学、Ch 7 血管壁力学提供集总参数描述 (lumped parameter), 是 white-box 建模的"零阶"近似。

2. 内容概述

本章分两大部分: 2.1-2.2 是循环系统的解剖、生理与机械系统性质; 2.3-2.4 是宏观/微观循环的集总参数模型 (WindKessel, 集总血管网络, Starling 过滤)。

Gasser 的核心立场是: 生物工程血管建模需要 0D/1D/2D/3D 的多尺度工具集——0D (WindKessel) 描述整体系统, 1D (血管段) 描述波传播, 2D/3D (轴对称/全三维, 留给 Ch 3-6) 描述局部机制。本章定位在 0D-1D 之间, 强调集总参数模型的解析可解性与多血管段网络的可装配性, 为血管介入 (如 EVAR 治疗, Example 2.3) 提供可工程化的工作流。

3. 核心方程与概念

3.1 循环系统的解剖生理 (Sect. 2.1)

  • 三主功能: Supply (输氧/营养) + Cleaning (代谢废物清除) + Immune response (免疫细胞输送)
  • 大循环 vs 肺循环: 大循环由左心室→主动脉→器官→腔静脉→右心房; 肺循环由右心室→肺动脉→肺微循环→肺静脉→左心房
  • 管径层级: 主动脉 (cm) → 动脉 → 小动脉 (10-100 μm) → 毛细血管 (5-8 μm) → 小静脉 (10-200 μm) → 静脉 (cm)
  • 血流速度变化 4 个数量级: 大动脉 cm/s 级别, 毛细血管 μm/s 级别 (Stokes 流)

关键血流动力学事实: 主动脉几乎贡献整个心血管系统的体积顺应性 (volume compliance), 其中胸主动脉占 85%。静脉系统储存约 60% 血液, 容量约为动脉的 50 倍。

3.2 血管壁的三层结构

  • Intima (内膜): 糖萼 + 单层内皮细胞 (EC) + 亚内膜层; 提供抗凝、低阻力界面
  • Media (中膜): 含 30-60% 平滑肌细胞 (SMC) + 5-25% 弹性蛋白 + 15-40% 胶原; 由 Medial Lamellar Units (MLU) 构成
  • Adventitia (外膜): 含 60-80% 胶原纤维 + 10-25% 其他结缔组织; 含 vasa vasorum (滋养血管)

力学意义: 弹性动脉富含弹性蛋白, 位于血管树近端; 肌性动脉含大量 SMC, 位于远端。中膜主要承受周向应力 (circumferential stress), 外膜锚定血管于周围组织。

3.3 内皮细胞 (EC) 与壁剪切应力 (WSS)

  • EC 关键作用:
  • 感知 WSS → 释放 NO → 弛张 SMC
  • 表达 adhesion molecules → 捕获白细胞
  • 维持 anti-thrombogenic 表面
  • WSS → NO 通路: L-arginine → eNOS → NO → 扩散到中膜 → SMC 弛张
  • Gasser 引用文献 [398, 74, 236]: 血管直径由 WSS 自稳态调控, 任何扰动都会触发重塑以恢复 WSS 至正常值

3.4 血液动力学调控 (Hemodynamic Regulation)

四大机制:

  1. Passive response — 纯被动变形, 毫秒级
  2. Vasoreactivity — SMC 收缩/舒张, 秒级
  3. Arteriogenesis — 血管直径/壁厚变化, 周-月级
  4. Angiogenesis — 新血管生成, 月-年级

自调控机制 (Autoregulation): - Myogenic (Bayliss 效应): 压力升高 → 血管收缩 → 阻力增加; 压力下降 → 血管舒张 → 阻力降低。8-20 kPa 压力范围内工作 - Flow-dependent: WSS 升高 → NO 释放 → 血管舒张 - 代谢性: 局部代谢物 (O₂↓, CO₂↑, 腺苷) 触发血管舒张

3.5 机械系统性质 (Sect. 2.2)

3.5.1 波速 (Eq. 2.1) — Moens-Korteweg 方程的变体

\[ c = \frac{1}{\sqrt{\rho D}}, \quad D = \frac{1}{A}\frac{dA}{dp} \]

其中 \(\rho\) = 血液密度, \(D\) = 血管面积可扩张度 (distensibility), \(A\) = 血管截面积。实际波速 = \(v \pm c\) (下/上游方向)。

Example 2.1 (水锤方程): 通过控制体 + 牛顿第二定律推导出

\[ \Delta p = \rho c \Delta v \quad \text{(water hammer equation)} \]

3.5.2 平均血压 (Eq. 2.3)

\[ p_{\text{mean}}(x) = \frac{1}{T}\int_0^T p(x,t)dt \approx \frac{1}{3}p_{\text{syst}} + \frac{2}{3}p_{\text{diast}} \]

关键观察 (Fig. 2.10): 压力降主要发生在最小直径的血管 (小动脉、毛细血管、小静脉); 毛细血管以下无脉动, 静脉系统无脉动。

3.5.3 脉搏压 (pp) 与反射波

  • 脉搏压 \(p_p = p_{\text{syst}} - p_{\text{diast}}\)
  • 前向波 + 后向波 (来自下游分支点, 主要是主动脉分叉) = 压力波
  • 老年人主动脉硬化 → 波速加快 → 后向波提前到达 → 收缩压增大 → 脉压增大

3.5.4 血管容量 (Eq. 2.4)

\[ C = \frac{\Delta V}{\Delta p} \]
  • 动脉总容量 \(C_{\text{art}} \approx 2\) ml/mmHg
  • 静脉总容量 \(C_{\text{ven}} \approx 100\) ml/mmHg (约 50 倍于动脉)
  • 主动脉占几乎全部容量, 其中胸主动脉占 85%
  • 主动脉硬化 → 弹性蛋白断裂 → 载荷转移至胶原 (硬 100-1000 倍) → 容量降低 → 左心室负荷增加

3.5.5 血管阻力 (Eq. 2.31) — Hagen-Poiseuille 定律

\[ R = \frac{\Delta p}{q} = \frac{128 \eta l}{\pi d^4} \]
  • \(d^4\) 依赖 → 直径微小变化引起阻力剧烈变化
  • 动脉阻力 \(R_{\text{art}} \approx 1\) mmHg·s/ml
  • 静脉阻力 \(R_{\text{ven}} \approx 0.06\) mmHg·s/ml

3.5.6 毛细血管传输 (Darcy 定律, Eq. 2.5)

\[ L_p = \frac{k}{\eta L} \]

其中 \(k\) = 毛细血管壁本征渗透率, \(L\) = 壁厚, \(\eta\) = 水的黏度。

3.6 WindKessel (WK) 模型家族 (Sect. 2.3)

核心思想: 把血管网络用 R (电阻) + C (电容) + L (电感) 的电路表示, 用电学类比建立血管系统的 0D 模型。

3.6.1 二元 WK (Eq. 2.6)

\[ q(t) = \frac{p(t)}{R} + C\frac{dp(t)}{dt} \]
  • \(R\) = 血管床总外周阻力, \(C\) = 总动脉容量
  • 给定压力 \(p(t)\) → 代数求 \(q(t)\); 给定 \(q(t)\) → 一阶 ODE 求 \(p(t)\)
  • 问题: 高频下 \(Z \to 0\), 与实验不符 (缺少惯性)

3.6.2 三元 WK (Eq. 2.20) — 工程标配

引入主动脉阻抗 \(Z_a = v_{pw} \rho / A\) (基于脉搏波速度), 解决高频行为:

\[ \frac{dp(t)}{dt} + \frac{p(t)}{RC} = Z_a \frac{dq(t)}{dt} + \frac{R + Z_a}{RC} q(t) \]
  • 特点: 高频 \(Z \to Z_a\) (实测可拟合), 是医疗器械测试回路的标配
  • 参数识别: 通过压力衰减法 (Eq. 2.15) 估 \(C\)、最小二乘法 (Eq. 2.17, 2.18) 估 \(R, C, Z_a\)

3.6.3 四元 WK (Eq. 2.30) — 含总动脉惯量 \(L\)

\[ \frac{dp(t)}{dt} + \frac{p(t)}{RC} = Z_a \frac{dq_Z(t)}{dt} + \frac{Z_a}{RC} q_Z(t) + \frac{q(t)}{C} \]
\[ \frac{dq_Z(t)}{dt} = \frac{dq(t)}{dt} - \frac{Z_a}{L} q_Z(t) \]

  • 作用: 改进低频行为; 用于准确预测压力波形
  • 应用: 临床压力测量校准

3.6.4 阻抗模/相 (Eq. 2.23, 三元 WK)

\[ |Z| = \sqrt{\frac{\alpha\beta}{1 + C^2 R^2 \omega^2}}, \quad \varphi = \arccos\left(\frac{R + Z_a + C^2 R^2 \omega^2 Z_a}{\sqrt{\alpha\beta}}\right) \]

其中 \(\alpha = 2 R Z_a + Z_a^2 + R^2(1 + C^2\omega^2 Z_a^2)\), \(\beta = 1 + C^2 R^2\omega^2\)。

3.7 血管段网络 (Sect. 2.3.2)

单段模型 (Fig. 2.23): 三元件 (R, C, L) 描述血管段的流体动力学

  • 阻力 (Eq. 2.31): \(R = 128\eta l / (\pi d^4)\) — Hagen-Poiseuille
  • 容量 (Eq. 2.32): \(C = 3 d^3 \pi l / (16 h E)\) — 线性弹性假设
  • 惯量 (Eq. 2.33): \(L = 4 \rho l / (d^2 \pi)\) — 牛顿第二定律

Example 2.8 (肾动脉适应): 部分肾切除 → 阻力增 \(\alpha\) 倍 → 流量减 → WSS ↓ + 周向应力 ↑ → 血管重塑直至两者恢复稳态值。导出 \(r/r_0\) 与 \(h/h_0\) 作为 \(\alpha\) 的函数。

装配规则 (Sect. 2.3.2.5): - 单管-单管: \(q_1 = q_2\), \(p_1 = p_2\) - 单管-双管分叉: \(q_1 = \xi_q q_2 + (1-\xi_q) q_3\), \(p_1 = p_2 = p_3\) - 压力降可调: \(p_1 = \xi_p p_2\), \(0 \le \xi_p \le 1\)

3.8 微循环: Starling 过滤与 glycocalyx-cleft 模型 (Sect. 2.4)

3.8.1 跨毛细血管浓度差 (Eq. 2.50)

\[ \Delta c = c_v - c_i = \frac{\sigma c_v}{\exp[Pe(L/d)] - 1} \cdot [\exp[Pe(L/d)] - \sigma] \]

其中 \(\sigma\) = Staverman 渗透反射系数, \(Pe = vd/D\) = Péclet 数。

3.8.2 经典 Starling 过滤 (Eq. 2.51)

\[ q_f = L_p (\Delta p - \sigma \Delta \pi) \]

其中 \(\Delta p = p_v - p_i\) = 跨壁流体静压差, \(\Delta \pi = \pi_v - \pi_i\) = 跨壁胶体渗透压差 (COP)。

四大 Starling 力: \(p_v, p_i, \sigma \pi_v, \sigma \pi_i\)。

3.8.3 当前修正: 几乎无 reabsorption (Sect. 2.4.2.3)

关键发现 (引用 [371]): 当稳态过滤被测量时, Starling 过滤模型失败——没有观察到线性关系, 没有静脉端 reabsorption。修正观点: 大部分间质液体由淋巴系统引流, 不是被血管重吸收。

3.8.4 非线性过滤 (Eq. 2.55)

\[ q_f = L_p \left[ \Delta p - \sigma^2 \pi_v^{\text{alb}} \frac{1 - \xi}{1 - \sigma \xi} \right], \quad \xi = \exp[Pe(L/d_{\text{alb}})] \]
  • 高 \(Pe\) 时, 扩散项消失, 关系退化为 \(q_f = L_p \Delta p\) (线性)
  • 低 \(Pe\) 时, \(\sigma^2 \pi_v^{\text{alb}}\) 项主导 (非线性)

3.8.5 Glycocalyx-cleft 双孔模型 (Sect. 2.4.2.6)

  • 糖萼层作为超细孔 (system of ultra-fine pores), 与 EC 间隙 (larger pores) 串联
  • 糖萼下方白蛋白浓度 \(c_{gc}\) 极低 → 有效 COP 高 → 净过滤压符合实测淋巴引流率
  • 解决了 Starling 经典模型的矛盾: 经典模型 overpredict 过滤通量; glycocalyx-cleft 模型拟合实测

3.9 淋巴系统 (Sect. 2.1.5)

  • 8 L/day 淋巴液; 一半被淋巴结吸收
  • 起始淋巴管 (10-60 μm) 由"橡叶状"EC + 弹性纤维锚定组成
  • 弹性纤维是机械感受器, 检测组织压力升高 → 开放 primary valves
  • 收集淋巴管 (1-2 mm) 含有瓣膜, 通过节律性收缩推动淋巴

3.10 压力-速度环 (Sect. 2.2.5)

  • 在 Argand 图上, 压力-速度环的切线斜率 = \(\rho c\) (加权波速)
  • 由水锤方程 \(\Delta p / \Delta v = \rho c\) (Eq. 2.2) 推出
  • 应用: 一种比 foot-to-foot 测量更精确的波速测量方法

4. 关键结论

  • [K1] 血管系统的总容量由主动脉主导, 总阻力由小血管主导: 这种分工决定了 WindKessel 二元模型的物理意义——\(C\) 来自主动脉, \(R\) 来自毛细血管前小动脉
  • [K2] Hagen-Poiseuille 的 \(d^4\) 依赖: 解释了为什么小血管阻力占主导, 主动脉直径即使变化 10% 也不会显著改变外周阻力
  • [K3] WSS 自稳态机制: EC 持续感知 WSS, 任何扰动触发血管重塑恢复 WSS 目标值; 这是 G&R (Ch 7) 的核心机制
  • [K4] 三元 WK 是临床工程标准: 既有解析可解性, 又能拟合实际阻抗-频率关系; 医疗器械测试回路的标配
  • [K5] 脉搏波反射决定压力波形: 老年高血压的脉压增大主要来自主动脉硬化 → 后向反射波提前到达, 而非心输出量增加
  • [K6] 血管段模型 (R, C, L) 的可装配性: 血管树可以像电路一样组装, 处理分叉、汇合、压力降; 这是 EVAR 治疗等介入规划的工作流
  • [K7] Starling 经典过滤模型需要修正: 稳态下没有静脉端 reabsorption, 大部分间质液由淋巴引流, glycocalyx-cleft 双孔模型更符合实测
  • [K8] 糖萼是有效渗透屏障: 而不是整个毛细血管壁; 糖萼损伤会大幅改变过滤特性, 与许多疾病的微循环功能障碍相关
  • [K9] 白蛋白贡献 80% 跨壁 COP: 决定了 Gasser 引用的"\(\pi_v \approx 22\) mmHg 来自白蛋白"的临床数值
  • [K10] 压力-速度环提供无创波速测量: 是连接临床血流动力学 (压力/速度波形) 与血管壁力学 (波速 → distensibility → 弹性模量) 的桥梁
  • [K11] 集总参数模型的能力边界: 不能描述波传播 (因此 Example 2.11 之后 Gasser 引出 1D 模型), 也不能解释局部机制 (留给 Ch 3-6)

5. 挑战和开放性问题

  • [C1] 集总参数模型的"阻抗"物理意义争议: Gasser 引用 [527] 指出三元 WK 中的 \(Z_a\) (主动脉阻抗) 物理意义有争议——它是基于脉搏波速度的"集总"概念, 不完全对应血管壁的局部阻抗。这影响其在精准血压分析中的应用。
  • [C2] Example 2.8 的 G&R 简化: 假设 WSS 和周向应力恢复至"稳态值", 但稳态值本身是经验值, 没有严格的物理定义。Gasser 引用 [271, 272] (Taber's law) 的工作可能给出更明确的稳态目标。
  • [C3] 糖萼损伤与微循环疾病: Ch 2 提到糖萼损伤会改变微循环过滤特性, 但未深入探讨其在 sepsis, diabetes, ischemia-reperfusion 等疾病中的作用机制。糖萼是当前微循环研究的热点。
  • [C4] Lymphatic 系统的力学建模: Gasser 在 2.1.5 坦承 "the biomechanics of the lymphatic system are not yet very well explored"。2.5 Summary 完全没有涉及淋巴系统的力学模型, 这是全书的盲点。
  • [C5] 零维模型与一维模型的边界: Gasser 引入了 WindKessel (0D) 和 vessel segment (1D, 离散), 但二者之间的波传播描述未提及。1D 连续波传播模型 (如 Langewouters 体系) 是 0D 与 1D 之间的桥梁。
  • [C6] 2.4.2.3 的修正观点 vs 2.4.2.5 的非线性模型 之间的微妙矛盾: Gasser 说"稳态下没有 reabsorption" (引 [371]), 但 2.4.2.5 又给出了基于 Starling 框架的非线性模型, 该模型预测有 reabsorption (Fig. 2.31 的左侧负压区)。这种张力没有显式调和。
  • [C7] Example 2.3 (EVAR 治疗) 的简化: 二元 WK 模型 + 恒定 \(q(t) = Q \sin^2(\pi t)\) 过于理想化。真实 EVAR 后, 支架改变几何, 局部 WSS 分布剧烈变化, G&R 启动, 这些都没有在 Example 2.3 体现。
  • [C8] 网络装配的 \(\xi_q\) 物理意义未明: 集总血管网络分叉处, 流分比 \(\xi_q\) 是状态变量, 由系统求解得到; 但实际生理中 \(\xi_q\) 由局部几何 (Murray 定律) 决定, 不是自由参数。
  • [C9] Example 2.1 的水锤方程 (Eq. 2.2) 假设 \(c \gg v\): 在大动脉中 \(c \sim 5-10\) m/s, \(v \sim 0.5-1\) m/s, 假设成立; 但在主动脉瓣狭窄等高流速病理下, \(v\) 可接近 \(c\), 假设破坏。
  • [C10] 2022 年后的更新: fractional-order WindKessel 模型 (考虑血管壁黏弹性)、深度学习参数识别 (替代最小二乘) 已在文献中出现, Gasser 应当提及。

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

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

本章是本人 SMC G&R 课题的集总参数背景:

  • 2.3.2 节血管段模型 (R, C, L) 与 1D 波传播的过渡 是本人研究的物理基础。Gasser 用了离散 0D, 但 G&R 课题需要在 1D 连续框架下处理波传播-重塑耦合。
  • Example 2.8 (肾动脉适应) 与本人研究的"G&R 重塑律"直接对应: Gasser 用了最简单的"恢复 WSS + 周向应力到稳态值"模型, 但实际 G&R 律 (Taber 1998, Humphrey 2002) 包含 mass 平衡、演化方程等更复杂形式。
  • Starling 过滤与糖萼-cleft 模型 (2.4.2) 与本人研究间接相关: 长期重塑需要考虑间质液压力的变化, 间质液压力受 Starling 过滤控制。
  • 2.3.1.4 压力衰减法 (Eq. 2.15) 是本人课题的"参数识别标准方法", 临床数据处理的参考。

6.2 Gasser 框架的局限性

  • 缺乏脉搏波分离方法: Gasser 用 forward + backward wave 描述脉搏波 (Fig. 2.11), 但没有给出 wave separation 的具体算法 (Westerhof 1972 的波分离是工程标配, Gasser 应当给出)。这是工程读者会卡住的地方。
  • 血管段模型缺少局部几何信息: R, C, L 都是平均几何参数, 不能预测局部狭窄、斑块、分叉处的应力集中。Gasser 在 2.5 提到了"up-scaling or homogenization of local vessel properties determines then the surrogate system parameters", 但没有给具体的多尺度方法。
  • 淋巴系统的力学建模完全缺失: 2.1.5 是生理描述, 2.3-2.4 完全跳过淋巴, 2.5 仍未提及。Gasser 把淋巴当"配角", 但实际上淋巴在 microcirculation 中扮演核心角色 (2.4.2.3 的修正就由淋巴承担)。
  • 0D 模型的物理边界: 0D 假设"血管段是 lumped", 但在大血管与毛细血管之间存在多个长度尺度 (弹性动脉、肌性动脉、小动脉、毛细血管), 每个尺度的物理机制不同, 用一个 0D 参数 (R) 描述会丢失尺度信息。
  • Example 2.1 (水锤方程) 的数学推导简短: Gasser 用了控制体+牛顿第二定律, 但跳过了控制体的非定常质量守恒, 对 0D 经验较少的读者会感到突兀。
  • 2022 年的盲点: 集总参数模型在精准医学时代的角色——Gasser 没有讨论个体化 WindKessel 参数估计 (用机器学习从 wearable 设备数据反演) 的可能性。

6.3 一个观察: "集总参数 vs 1D-3D" 的工具选择哲学

Gasser 在 2.5 Summary 的最后一段坦承: 集总参数不能解释局部机制, 多尺度时空描述在 Ch 3+ 展开。这种"先用 0D 简化, 再 1D, 再 2D-3D" 的工具选择哲学, 是工程生物力学的标准方法论, 但有一个隐含假设: 所有血管段都能用同一种本构/同一种网格密度 描述。

实际上: 主动脉 (cm 级, 弹性动脉, 厚壁) 与毛细血管 (μm 级, 单细胞壁) 的力学机制天差地别。Gasser 应当更明确地说明: Ch 3-7 适用于大-中等血管, 微血管 (毛细血管) 已被本章的 Starling 模型处理, 不能直接套用 Ch 5 (管道血管) 的本构。

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

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

  • [X1] Gasser T.C. (2022) Vascular Biomechanics. Springer, Chapter 2. — 本章作者本人所著
  • [X2] Milnor W.R. (1989) Hemodynamics. Williams & Wilkins — 血流动力学经典教材, 引用 394
  • [X3] Hall J.E., Hall M.E. (2020) Guyton and Hall Textbook of Medical Physiology. Elsevier — 生理学经典, 引用 548
  • [X4] Westerhof N., Lankhaar J.W., Westerhof B.E. (2006) The arterial Windkessel. Med Biol Eng Comput 44:35-42 — 三元 WindKessel 的标准参考, 引用 596
  • [X5] Frank O. (1899) Die Grundform des arteriellen Pulses. Z Biol 37:483-526 — 二元 WindKessel 模型的原始论文, 引用 173
  • [X6] Moens A.I. (1878) Die Pulskurve. Leiden — Moens-Korteweg 波速公式的原始论文
  • [X7] Korteweg D.J. (1878) Über die Fortpflanzungsgeschwindigkeit des Schalles in elastischen Röhren. Ann Phys 241:525-542
  • [X8] Starling E.H. (1896) On the absorption of fluids from the connective tissue spaces. J Physiol 19:312-326 — Starling 过滤的原始论文, 引用 524
  • [X9] Staverman A.J. (1951) The theory of measurement of osmotic pressure. Rec Trav Chim Pays-Bas 70:344-352 — 渗透反射系数的原始论文, 引用 526
  • [X10] Levick J.R., Michel C.C. (2010) Microvascular fluid exchange and the revised Starling principle. Cardiovasc Res 87:198-210 — 现代 Starling 修正, 引用 338
  • [X11] Weinbaum S., Tarbell J.M., Damiano E.R. (2007) The structure and function of the endothelial glycocalyx layer. Annu Rev Biomed Eng 9:121-167 — 糖萼作为有效渗透屏障, 引用 543
  • [X12] Adamson R.H., Lenz J.F., Zhang X., Adamson G.N., Weinbaum S., Curry F.E. (2004) Oncotic pressures opposing filtration across non-fenestrated rat microvessels. J Physiol 557:889-907 — 实验证据, 引用 371
  • [X13] Pries A.R., Secomb T.W. (2014) Making microvascular networks work: angiogenesis, remodeling, and pruning. Physiology 29:446-455 — 微血管重塑的现代综述
  • [X14] Hales S. (1733) Statistical Essays: Containing Haemastaticks. Innys & Manby — 第一次血压测量
  • [X15] Bayliss W.M. (1902) On the local reactions of the arterial wall to changes of internal pressure. J Physiol 28:220-231 — Bayliss 效应 (myogenic regulation), 引用本节
  • [X16] Murray C.D. (1926) The physiological principle of minimum work. I. The vascular system and the cost of blood volume. Proc Natl Acad Sci USA 12:207-214 — Murray 定律 (血管树分叉), 引用本节隐含
  • [X17] Caro C.G., Pedley T.J., Schroter R.C., Seed W.A. (1978) The Mechanics of the Circulation. Oxford University Press — 循环系统力学经典, 引用 91
  • [X18] Laughlin M.H., Newcomer S.C., Bender S.B. (2008) Importance of hemodynamic forces as signals for exercise-induced vascular remodeling. Am J Physiol Heart Circ Physiol 295:H1209-H1210 — WSS 调节血管重塑, 引用 74
  • [X19] Ku D.N. (1997) Blood flow in arteries. Annu Rev Fluid Mech 29:399-434 — 动脉血流的经典综述, 引用 609
  • [X20] Taylor C.A., Figueroa C.A. (2009) Patient-specific modeling of cardiovascular mechanics. Annu Rev Biomed Eng 11:109-134 — 1D-3D 多尺度血管建模综述