第3章：骨的生物力学（Biomechanics of Bone）

章节概述

本章由 Carol A. Oatis 撰写，进入具体生物组织的力学分析。本章先回顾骨作为结缔组织的生物学（结构、组成、细胞、形成与吸收），再把骨视作工程材料讨论其各向异性（anisotropic）、非均匀性（inhomogeneous）、粘弹性（viscoelastic）的力学性质以及骨强度（strength）、断裂韧性（fracture toughness）等核心参数；接着讨论骨的载荷模式（拉伸、压缩、剪切、弯曲、扭转、组合载荷）下的力学响应，并比较皮质骨（cortical bone）与松质骨（trabecular / cancellous bone）的差异；最后讨论骨的年龄、活动、激素水平对骨力学性质的影响以及骨折愈合（fracture healing）的生物力学背景。本章建立了后续 Ch10、Ch13、Ch16、Ch19、Ch25、Ch28、Ch31、Ch34、Ch37、Ch40、Ch43、Ch46 各章"力分析"章中骨反作用力、骨应力计算的力学基础。

关键问题与研究动机

本章要回答的核心问题有五个：(1) 骨是由什么组成的？哪些成分（羟基磷灰石、胶原、水）决定了骨的硬度和韧性？(2) 骨为什么是各向异性的？哈弗斯系统（osteon）的方向性如何导致骨的纵向模量高于横向模量？(3) 皮质骨和松质骨的力学性质有何不同？为什么长骨骨干是皮质骨、椎体和长骨末端是松质骨？这种分布与载荷传递（load transfer）的关系是什么？(4) 骨的强度和断裂韧性受哪些因素影响？哪些因素是生理可调节的（如活动、营养、激素）？哪些是不可调节的（如年龄、性别、基因）？(5) 骨在体内承受什么类型的载荷？弯曲和扭转的强度为何比纯拉伸/压缩低得多？

研究动机：骨科临床中遇到的骨折、骨质疏松、骨折愈合延迟、应力性骨折（stress fracture）、骨重建（remodeling）等，都需要把骨的力学性质和生理载荷结合来判断。临床医生只有在知道"骨有多硬、能承受多大应变"时，才能判断"这次摔倒会不会骨折"、"这个骨折需要多久固定"、"这个应力集中是否合理"。

主要公式与推导

皮质骨的纵向（longitudinal）杨氏模量约 \(E_L \approx 17\text{ GPa}\)，横向模量约 \(E_T \approx 12\text{ GPa}\)；剪切模量约 \(G \approx 3.3\text{ GPa}\)。松质骨的杨氏模量在 \(0.05\text{ GPa}\)（低密度）到 \(3\text{ GPa}\)（高密度）之间，强烈依赖于表观密度（apparent density）\(\rho\)。经验公式

\[ E = C \rho^{n} \]

其中 \(C\) 是常数，\(n\) 约在 1.5-2.0 之间。极限强度也服从类似标度律

\[ \sigma_u = D \rho^{m} \]

骨在应变率提高时表现得"更刚"、强度更高、断裂应变更小，反映其粘弹性。标准线性固体模型（standard linear solid）

\[ \sigma + \tau_\varepsilon \dot{\sigma} = E_0(\varepsilon + \tau_\sigma \dot{\varepsilon}) \]

其中 \(\tau_\varepsilon\) 和 \(\tau_\sigma\) 是两个特征时间常数。

骨在弯曲下的最大正应力

\[ \sigma = \frac{M c}{I} \]

其中 \(M\) 是弯矩，\(c\) 是距中性轴最远距离，\(I\) 是横截面二阶面积矩（second moment of area）。在扭转下

\[ \tau_{\max} = \frac{Tr}{J} \]

其中 \(T\) 是扭矩，\(r\) 是外半径，\(J\) 是极惯性矩。弯曲和扭转在长骨中往往是组合载荷。

骨的应变能密度（韧性指标）

\[ U = \int_0^{\varepsilon_f} \sigma\, d\varepsilon \]

皮质骨的 \(U\) 约 \(6\text{ MJ/m}^3\)（比木材低、比陶瓷高）。断裂韧性（fracture toughness）\(K_{IC}\) 约 \(2-12\text{ MPa}\sqrt{m}\)，是描述骨抵抗裂纹扩展能力的参数。Griffith 准则

\[ \sigma_f = \sqrt{\frac{2 E \gamma_s}{\pi a}} \]

其中 \(\gamma_s\) 是单位面积表面能，\(a\) 是裂纹长度。

关键算法与建模方法

本章介绍的方法主要是：(1) 单轴拉伸/压缩试验测皮质骨的 \(E\) 和 \(\sigma_u\)，试件方向相对骨长轴的取向决定了测得的是纵向还是横向模量；(2) 三点/四点弯曲试验测弯曲模量和弯曲强度；(3) 扭转试验测剪切模量；(4) 压痕试验（indentation test）测局部硬度；(5) μCT 配合有限元（microFE）测松质骨的三维结构和表观力学性质。

对于松质骨，由于其多孔结构，常用 Ashby-Gibson 模型描述其模量与相对密度的关系

\[ \frac{E^*}{E_s} = C \left(\frac{\rho^*}{\rho_s}\right)^n \]

其中 \(E^*\) 是表观模量、\(E_s\) 是固体材料模量、\(\rho^*/\rho_s\) 是相对密度。

应变率（strain rate）效应：骨的 \(\sigma_u\) 在低应变率（\(10^{-3}\text{ s}^{-1}\)）下约 100-150 MPa，在高应变率（\(10^{3}\text{ s}^{-1}\)）下提高约 50%。这种率敏感性有重要的临床意义：摔倒（低应变率）vs. 车祸（高应变率）下同一根骨可能表现出完全不同的强度。

年龄与活动效应：骨的力学性质随年龄下降（特别是绝经后女性的骨密度下降），规律负重活动（如步行、跑步）可以增加骨密度和骨强度，但极端负荷（过度训练）可能导致应力性骨折。

骨折愈合的生物力学：骨折愈合经历了血肿机化期、纤维骨痂期、骨性骨痂期、改建期四个阶段，每个阶段的力学环境（稳定性、载荷传递）都影响愈合结果。绝对稳定（rigid fixation）下发生直接骨愈合（primary bone healing），相对稳定下发生间接骨愈合（secondary bone healing via callus）。

主要结论

本章的主要结论是：(1) 骨是各向异性的复合材料，纵向刚度和强度显著高于横向——这种各向异性源于哈弗斯系统沿骨长轴排列。(2) 皮质骨与松质骨的力学性质差异巨大；皮质骨致密坚硬负责承载高应力，松质骨多孔轻质负责能量吸收和载荷分布。(3) 骨的强度和模量是"密度依赖"的——松质骨的 \(E \propto \rho^{1.5-2.0}\) 关系意味着密度下降 30% 会导致模量下降约 60%，这是骨质疏松易于骨折的核心机制。(4) 弯曲和扭转是骨折的常见原因——长骨骨干的中段在侧向力下产生弯曲，而扭转时剪应力在中性轴上最大；这两点决定了长骨骨折线的典型位置和方向。(5) 骨的力学性质是动态的——年龄、激素、活动、营养都影响骨的强度；机械刺激（Wolff 定律）下骨会重塑以适应载荷。

挑战与开放问题

本章的挑战集中在骨力学性质的体内测量和骨重建的预测上。体内骨的真实应力状态难以直接测量，目前只能通过应变计（strain gauge）或有限元模型（FEM）估算。骨重建的数学模型（如 Stanford 的 μCT 配 FEM）虽已能预测骨密度变化，但要准确预测患者特异性的骨折风险仍需大量临床验证。另外，应力性骨折（特别是运动员和军人）的预测、骨-假体界面的应力集中、骨质疏松治疗的力学阈值都是未完全解决的问题。

个人反思与批判性分析

本章是 Ch1 刚体力学和 Ch2 材料力学在生物组织的第一个应用。从教学角度，Oatis 把"骨是工程材料"这个观点讲得清楚——把骨视作多孔复合材料、用密度标度律、断裂力学、应变率效应来描述，符合 1990 年代以来的生物力学共识。

但有两点值得批评：(1) Oatis 把骨重建基本局限于"载荷适应性"（Wolff 定律）的论述，没有充分展开"细胞-力学耦合"的多尺度机制——这是 2009 年后生物力学研究的活跃领域（van der Meulen、Huiskes 等的工作）。(2) 骨的微损伤（microdamage）和靶向重建（targeted remodeling）概念在 Oatis 中提及较少，而这对解释疲劳骨折和老年脆性骨折至关重要。

实际研究含义：做骨折风险评估（如髋部骨折），传统 DXA 测骨密度（BMD）只解释了 60-70% 的骨折方差——剩下的变异度由骨的几何结构（髋轴长度 HH）、骨微结构（皮质厚度）、微损伤积累、跌倒倾向等共同决定。Oatis 提到的"几何因素"是后续研究中"股骨近端几何参数"研究的前身。

重要参考文献

[1] Currey JD. Bones: Structure and Mechanics. 2nd ed. Princeton University Press.（骨力学结构的经典参考） [2] Cowin SC. Bone Mechanics Handbook. 2nd ed. CRC Press.（骨力学的工具书，覆盖微观-宏观多尺度） [3] Wolff J. The Law of Bone Remodeling. Springer.（Wolff 定律的原始论述） [4] Ashby MF, Gibson LJ. Cellular solids: structure and properties. Pergamon.（多孔固体材料力学，松质骨的 Ashby-Gibson 模型） [5] Huiskes R, Ruimerman R, van Lenthe GH, Janssen JD. Effects of mechanical forces on maintenance and adaptation of form in trabecular bone. Nature. 2000;405(6787):704-706.（骨重建的细胞-力学耦合） [6] van der Meulen MCH, Huiskes R. Why mechanobiology? A survey article. J Biomech. 2002;35(4):401-414.（力学生物学综述） [7] Hernandez CJ, van der Meulen MCH. Understanding bone strength: size matters, but so does architecture. BoneKEy. 2010;10:17-22.（骨几何与骨强度） [8] Kanis JA, McCloskey EV, Johansson H, Oden A, Ström O, Borgström F. Development and use of FRAX in osteoporosis. Osteoporos Int. 2010;21(2):407-413.（骨折风险评估工具）