跳转至

第五章 实验方法 (Experimental Methods)

§1 作者

本章是Ch 4本构理论的"实验对应面"。Humphrey的策略是"先哲学后技术"——5.1节花了3页讲"理论-观察-实验"的循环 (DEICE法)，5.2—5.4节才是具体测量方法 (VDA, marker tracking, 双轴测试机, 多轴测试机, 计算机辅助)。作者特别强调"理论必须走在实验前面，但实验不必等待理论"——这是Bernard 1865 + Duhem 1894的现代版重述。

§2 内容概述

分七节： - 5.1 总论 (3小节)：科学方法哲学 (Bernard-Duhem) → 测量系统组成 (传感器+信号调节+记录) → 量纲分析 (Buckingham Pi)。 - 5.2 应变测量 (2小节)：视频测距仪 (VDA) + marker tracking (MRC, 互相关) → 三角形/四边形标记法求 \(\mathbf{F}\) (Eq. 5.6)。 - 5.3 载荷测量：应变计、力传感器 (load cell)、压力传感器、热电偶的原理。 - 5.4 计算机辅助实验：生物表面双平面相机 (biplane)、马蹄形 (horse-shoe) 试件、伺服液压测试系统。 - 5.5 参数估计与统计：Marquardt-Levenberg非线性最小二乘 + 多变量回归 + E-不等式约束。 - 5.6 习题 (15题)：从Pi定理应用到MRC算法实现。 - 5.7 参考文献 (50条)。

§3 核心方程与概念

3.1 测量系统与误差 (Eq. 5.1—5.5)

(1) 测量系统三件套 (Fig. 5.1) 传感器 (transducer) + 信号调节 (conditioner, 放大+滤波) + 记录器 (recorder, A/D+数字存储)。4个关键特性：线性 (calibration简单)、稳定 (calibration不变)、滞回小、信噪比高。

(2) 准确度 vs 精确度 (5.1.2末段) - 准确度 (accuracy)：测量值逼近真实值的程度 → \(|\bar{x} - R|\) 或单点最差 \(|x - R|\)。 - 精确度 (precision)：重复性 → 标准差或最大偏差。 经典示例：5个电压测量 104.2, 103.4, 105.2, 103.7, 104.5V，真值 100V → \(\bar{x} = 104.2V\)，准确度 = 4.2% (按第一派) 或 5.2% (按最差点)，精确度 = 1.0V/1%。

(3) 三类误差 - Gross (粗大)：违反协议 → 靠设计避免。 - Systematic (系统)：标定错误、仪器偏置 → 例中均值偏离真值就是系统误差。 - Random (随机)：电源波动、电磁干扰 → 滤波或屏蔽消除。

(4) 量纲分析 (Buckingham Pi, Eq. 5.1—5.4) 五步法：列变量 → 写维度 → 选尺度 → 写Pi群 → 重写方程。血管力学尺度选择示例（5.1.3末段）：\(L_* = r_i\) (变形内半径) 或 \(\sqrt{A}\)，\(T_*\) = 周期或松弛时间，\(M_* = \rho V\)。意义：8.2.4节 颅内动脉瘤动力稳定性问题会用此法。

3.2 应变测量 (Eq. 5.5—5.8)

(5) VDA (Video Dimension Analyzer, 5.2.1首段) 最早的非接触有限应变测距设备 (Yin et al. 1972)。线扫描，输出电压 ∝ 两边缘距离。缺点：只能测单方向平均应变，对刚体旋转敏感 → 后被 marker tracking 取代。

(6) MRC (Method of Rows and Columns) 最简标记质心算法。黑白标记 + 灰度8-bit → 每行/列的"最暗"位置即质心 \(x, y\)。条件：对比度必须高。

(7) Downs et al. 1990 二步相关搜索 (Eq. 5.5) $\(\rho = \sum_{m,n=1}^{7} R_{mn} T_{mn}\)$ - 粗搜索：\(35 \times 35\) 像素 → \(5 \times 5\) 网格 \(7 \times 7\) tile → 25个相关值。 - 精搜索：粗定位点 \(13 \times 13\) → 49个相关值。 总计算量 = 74 个相关值（直接 \(35 \times 35\) 是 841 个）→ 80386 CPU 31ms 完成 → 实时 (30Hz)。这是"硬件约束驱动算法设计" 的经典例子。

(8) 三角形标记法求 \(\mathbf{F}\) (Eq. 5.6—5.8) 在 \(\kappa_0\) 三点 \(A, B, C\)、在 \(\kappa_t\) 三点 \(a, b, c\)。设 \(dX^{(1)} = B - A\), \(dX^{(2)} = C - A\), \(dx^{(1)} = b - a\), \(dx^{(2)} = c - a\)。则： $\(\begin{pmatrix} dx_1^{(1)} & dx_1^{(2)} \\ dx_2^{(1)} & dx_2^{(2)} \end{pmatrix} = \begin{pmatrix} F_{11} & F_{12} \\ F_{21} & F_{22} \end{pmatrix} \begin{pmatrix} dX_1^{(1)} & dX_1^{(2)} \\ dX_2^{(1)} & dX_2^{(2)} \end{pmatrix}\)$ \(\mathbf{F} = dx \cdot dX^{-1}\) (2x2 矩阵求逆) → 然后极分解求 \(\mathbf{U, R, \Phi}\)。Eq. 5.8 给 2D 旋转角 \(\phi\) 的直接计算公式 (不用 \(\mathbf{U}\) 平方根)： $\(\phi = \tan^{-1}\left[\frac{F_{12} - F_{21}}{F_{11} + F_{22}}\right]\)$ 意义：Ch 5、Ch 7 心肌 marker 跟踪 (McCulloch et al. 1987) 都用此公式。

3.3 加载测量 (5.3节)

(9) 应变计 (Strain gauge) 金属箔电阻随应变变化 → 电桥读数。关键：温度补偿 (温度系数与待测材料接近)；泊松横向补偿 (双臂/全桥)。

(10) 测压元件 (Load cell) - Pancake (扁平)：测轴向力。 - S-type：S形测拉+压。 - 多轴：3分量力矩。 心血管测试关键：测压元件量程选 1-100 N 时分辨率 0.01-0.1 N；血管 biaxial 试验机多用 Instron 8800 + custom clamps。

(11) 压力传感器 (Catheter-tip, strain-gauge, optical) - Millar catheter：3F 微型压力传感器，频响 > 10 kHz。 - 血管内血压测量的"金标准"。

(12) 热电偶 (Thermocouple) 两种金属接合 → 温差电动势 (Seebeck effect)。心血管：激光/射频消融的局部温度测量。

3.4 计算机辅助实验 (5.4节)

(13) 双平面相机系统 (Biplane) 2台X射线相机同步 → 三维标记位置 (Eq. 5.38 校准公式)。应用：In-vivo 心肌 marker 跟踪 (Waldman, McCulloch 1987)；缺点：辐射剂量。

(14) 马蹄形 (Horse-shoe) 试件 环形血管/输尿管样试件 → 一端固定 + 另一端加力 → 测内压 + 长度 + 扭角。优点：可测"无接触"的天然形状 (zero-stress state)。

(15) 伺服液压 (Servo-hydraulic) 测试机 - 闭环控制：位移/力/应变率反馈控制 (Instron 8800, MTS, EnduraTec)。 - 软件：LabView, TestStar, custom C。 - 关键：试件的对中 (alignment) → 偏心会产生弯矩 → 影响数据。

3.5 参数估计 (5.5节)

(16) Marquardt-Levenberg 算法 混合 Gauss-Newton (靠近最小值时快) + 梯度下降 (远离时稳定)： $\(\Delta \beta = \frac{1}{\lambda + \text{diag}(J^T J)} J^T (y_{exp} - y_{theory})\)$ \(\lambda \to 0\) 是 Gauss-Newton，\(\lambda \to \infty\) 是梯度下降。优点：稳健、不需二阶导。

(17) 目标函数 (Eq. 4.18, 5.42) $\(\varepsilon = \sum_{\alpha=1}^{m} \omega_\alpha \left[ \left( t_{11}^{exp} - t_{11}^{theory} \right)^2 + \left( t_{22}^{exp} - t_{22}^{theory} \right)^2 \right]\)$ \(\omega_\alpha\) 权重 → 平衡不同尺度数据。

(18) E-不等式约束 (Ch 3.4.4末段) 在 Marquardt-Levenberg 中用先验不等式约束参数： - \(c_1 > 0\) (neo-Hookean 必须正) - \(\lambda_1 > \lambda_2 \Rightarrow t_{11}(\lambda_1, \lambda_2) > t_{22}(\lambda_1, \lambda_2)\) 意义：保证拟合出的参数在所有应变状态下都物理合理。

§4 关键算法或建模方法

  1. MRC + Downs 二步搜索 (Eq. 5.5)：标记质心提取，实时性保证 30Hz。Ch 7 血管实验、Ch 10 心肌 marker tracking 都直接用此思路。
  2. \(\mathbf{F}\) 计算 (Eq. 5.6—5.8)：2D三角形标记的代数公式；Ch 6 有限元结果验证 (与体元 \(\mathbf{F}\) 对比)。
  3. Buckingham Pi 五步法 (Eq. 5.1—5.4)：8.2.4节动脉瘤动力稳定性分析会用此法找出关键无量纲数 (e.g., Reynolds, Womersley, \(\beta\))。
  4. Marquardt-Levenberg (5.5节)：所有"双轴→ \(W\)"数据拟合的工具；Ch 7 实际血管本构拟合也用。
  5. E-不等式 + Marquardt (Ch 3.4.4 + 5.5)：组合起来给"物理合理 + 数学最优"的参数。
  6. 2D 旋转角直接公式 (Eq. 5.8)：避免 \(\mathbf{U}\) 平方根 → 比 3D 算法快 ~10×。Ch 5 的"省时trick"。

§5 关键结论

  1. 理论先于实验，观察-实验-理论循环 (DEICE法)：5.1节用 Bernard/Duhem 的19世纪引文反复强调。
  2. 不可压缩血管/组织的应变测量 = 非接触 → VDA / marker tracking / MRI tags / 散斑干涉。接触式应变计对小应变材料 (金属) 适用，对软组织不适用。
  3. 2D 三角形 + 4个marker = 测 \(\mathbf{F}_{2\times 2}\) (Eq. 5.6) — 最少标记数 = 3，4个更稳。Ch 6 有限元结果验证 (5.6末段给出).
  4. Buckingham Pi 把"十几个变量"压成"3-5个Pi群" → 5.1.3末段给血管尺度选择。
  5. 生物力学实验 = "仪器设计" + "数据采集" + "参数估计" 三阶段。Humphrey 5.1.1首段提醒："如果仪器不满足实验需求，应先设计仪器，而不是迁就仪器"。
  6. 精度-准确度-误差三类必须明确报告 (5.1.2末段)。Humphrey 给出"两派准确度定义"对比——提醒读者明确使用哪种。
  7. Duhem (1894) 警告"理论不应在实验进行时发言"——但 Duhem 区分"应用实验" (演示理论) vs "测试实验" (验证理论)。这与 5.1.1 末段呼应"DEICE 5步是循环"。

§6 挑战和开放性问题

  1. 软组织 marker 跟踪的"标记丢失"问题：Downs 算法假设下一帧在 \(35 \times 35\) 区域内 — 高速应变率测试 (e.g., 颅内动脉瘤破裂 1000/s) 时容易丢失。5.2.1末段提到用高速相机+离线分析，但实时控制仍是难题。
  2. MRI tagging 的分辨率限制：3D 心肌 marker 在 MRI tagging 中 ~1mm 分辨率，但心肌细胞 ~10μm → 测量的是"宏观平均"，不是"细胞级应变"。Ch 10 心肌本构中"宏观 vs 细胞"的桥接是开放问题。
  3. 2D vs 3D 测量的限制：Eq. 5.6是2D → 平面外应变 (out-of-plane) 不能直接测。Ch 7 血管壁 3层 (intima, media, adventitia) 的 3D 应变仍需 CT/MRI。
  4. Buckingham Pi 的"经验"问题：尺度选择 (5.1.3 step 3) 没有系统方法 — "靠经验和创造" (Humphrey原话)。初学者往往选错尺度导致 Pi 群没有物理意义。
  5. 仪器商业化的"市场小"问题：5.1.1首段提到"专门为生物力学设计的测试系统往往不商业化" → 实验室必须自己造 (e.g., Humphrey 自己实验室的 biaxial 试件) — 这门槛很高，研究难推广。
  6. 参数估计中的"过拟合"风险：Mooney-Rivlin (\(c_1, c_2\) 2个参数) vs Ogden (\(N\) 个) vs Fung 指数型 (5+ 参数) — 拟合"描述能力"随参数增加而增加，但"预测能力"未必。这是 Ch 4.1.5 末段"Evaluate"步骤的反复警告。
  7. E-不等式约束的实现细节：5.5节说"在 Marquardt-Levenberg 中用约束"——但没有给出具体代码。读者面临"如何在非线性优化中实现 \(\lambda_1 > \lambda_2 \Rightarrow t_{11} > t_{22}\)"的实际困难。

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

Humphrey的实验方法章有一个被低估的trick：把哲学和工程紧密结合。5.1节引Bernard和Duhem 19世纪的话，看似离题，但实际上为后面的仪器设计选择提供方法论支持——读者在5.4节看到"先设计仪器再实验"时，不会觉得突兀。

值得自己复现的算法： - MRC marker 质心提取 (5.2.1中段)。 - Downs 二步相关搜索 (Eq. 5.5)。 - 2D 三角形 \(\mathbf{F}\) 计算 (Eq. 5.6—5.8)。 - Marquardt-Levenberg (用 SciPy curve_fit 或 MATLAB nlinfit 实现)。

与作者的潜在对话： - 5.2节末段提到"高速相机+离线分析"是高速 marker tracking 的解决方案，但未给具体相机型号和帧率——Humphrey实验室用什么相机？50Hz? 1000Hz? - 5.4节 biplane 系统的辐射剂量问题：对 in-vivo 心脏研究，每次实验 ≤ 5 mGy？5.4节没给具体数字。这对临床应用是关键约束。 - 5.5节 Marquardt-Levenberg 的"\(\lambda\) 选择策略"没有给具体规则——是固定减小？是用 armijo 条件？Humphrey 可能假设读者会查 Press et al. Numerical Recipes，但这章的工程导向应该给一两个具体建议。 - 5.1.1末段"动物实验 → 人类测试"的伦理框架 完全没提 — 5.5节是 homo sapiens 临床研究的入口，应简略提 IRB / Helsinki Declaration。

总结性批评： - 5.2节marker tracking侧重 VDA 和 MRC，但完全没讲现代 deep-learning marker tracking (e.g., DeepLabCut, HRNet)。2002年写书时这些不存在，但2026年回头看是一大缺漏——读者读到这章会问"现代有什么新方法"。 - 5.3节载荷测量只给 strain gauge / load cell / pressure / thermocouple 四种 — 没有 IMU (惯性测量单元) (现代可穿戴设备) 或 光纤布拉格光栅 (FBG) (现代可植入应变传感器)。 - 5.4节 Computer-Aided Experimentation侧重 biplane + horse-shoe + servo-hydraulic — 没有提到"数字孪生" (digital twin) 概念，虽然这是心血管固体力学的未来方向。 - 5.5节 Parameter Estimation and Statistics没给"协方差矩阵"和"参数不确定性"——Marquardt-Levenberg 的输出应包括参数的 95% 置信区间，但 Humphrey 只说"找到 best-fit values"。 - 5.6节习题 15题覆盖偏浅——多数是理论复现 (5.1.1节的Pi定理应用)，缺少"实际数据拟合"或"marker tracking 模拟"。

§8 重要参考文献

  • [X1] Yin FCP, Tompkins WR, Peterson KL, Intaglietta M (1972) A video dimension analyzer for dynamic measurement of the cross-sectional shape of blood vessels in situ. Comput Biol Med 2:283-292. — VDA 设备的最早文献。
  • [X2] Waldman LK, Fung YC, Covell JW (1985) Transmural myocardial deformation in the canine left ventricle. Circ Res 57:152-163. — 心脏 in-vivo marker tracking 的早期工作。
  • [X3] McCulloch AD, Smaill BH, Hunter PJ (1987) Left ventricular epicardial deformation in isolated arrested dog heart. Am J Physiol 252:H233-H241. — Eq. 5.6三角形法的具体应用。
  • [X4] Ingels NB, Daughters GT, Stinson EB, Alderman EL (1989) Measurement of midwall myocardial dynamics in intact man by radiography of surgically implanted markers. Circulation 59:849-855. — in-vivo 心脏 marker 跟踪。
  • [X5] Downs J, Halperin HR, Humphrey JD, Yin FCP (1990) An improved video-based computer tracking system for soft biomaterials testing. IEEE Trans Biomed Eng 37:903-907. — Eq. 5.5 二步相关算法的来源。
  • [X6] Humphrey JD, Vawter DL, Vito RP (1986, 1987) Quantification of strains in biaxial soft tissues. J Biomech 等 — MRC + threshold 算法的应用。
  • [X7] Lanir Y, Fung YC (1974) Two-dimensional mechanical properties of rabbit skin. J Biomech 7:29-34. — 双轴测试的最早应用。
  • [X8] Pirolo JS, Yoganandan N, Pintar FA, Reinartz JM (1992) An algorithm for marker tracking in non-rigid bodies. Comput Methods Biomech Biomed Eng — 后处理 marker 跟踪。
  • [X9] Mandel J (1964) The Statistical Analysis of Experimental Data. Wiley. — 准确度/精确度两派定义的来源 (5.1.2末段)。
  • [X10] Holman JP (1989) Experimental Methods for Engineers. McGraw-Hill. — 误差/不确定度标准的工科教材。
  • [X11] Bell JF (1973) The Experimental Foundations of Solid Mechanics. Handbuch der Physik VIa/1, Springer. — 5.1.1首段引文 "distortions of an overemphasis upon technique"。
  • [X12] Bernard C (1865) Introduction to the Study of Experimental Medicine. Macmillan. — 5.1.1节观察者/实验者区分的来源。
  • [X13] Duhem P (1894) The Aim and Structure of Physical Theory. Princeton. — 5.1.1节末段"理论应在实验门外等待" 的来源。
  • [X14] Smutz M, Downs J, Nicosia MA, Phatak N, Carr R, Huang K, Hwang K, Yin FCP (1996) An algorithm for automated tracking of multiple markers in cine MRI. J Magn Reson Imaging 6:105-114. — 5.2.2末段提到。
  • [X15] Butler DL, Grood ES, Noyes FR, Zernicke RF, Brackett K (1984) Effects of structure and strain measurement technique on the material properties of young human tendons and fascia. J Biomech 17:579-596. — 高速相机+离线分析的应用。