第4章：骨骼肌的生物力学（Biomechanics of Skeletal Muscle）

章节概述

本章由 Carol A. Oatis 撰写，是肌肉骨骼系统的核心章节之一。本章先回顾骨骼肌的微结构（肌节、肌原纤维、粗细肌丝、肌管系统）和宏观结构（肌束、肌外膜、肌腱），再讨论肌肉的力学性质：长度-张力关系（length-tension relationship）、力-速度关系（force-velocity relationship）、肌腱的串联弹性（series elasticity）、肌肉的募集（recruitment）与发放频率（firing frequency）调节；然后展开肌肉的杠杆效应（moment arm）、肌肉的形态学参数（生理横截面积 PCSA、纤维长度、羽状角 pennation angle）；最后讨论肌肉的训练适应（training adaptation）、去适应（detraining）、废用性萎缩（disuse atrophy）、年龄效应。

关键问题与研究动机

本章要回答的核心问题有七个：(1) 肌力大小由哪些因素决定？同样的肌肉在不同收缩长度、不同收缩速度下能产生多大力？(2) 肌肉的力-长度关系（force-length）和力-速度关系（force-velocity）的机制是什么？这些曲线如何预测肌肉在不同关节角度、不同运动速度下的力输出？(3) 肌肉的力矩（moment）和关节力矩的关系是什么？moment arm 如何随关节角度变化？同一块肌肉在不同角度下可能成为主动肌（agonist）、拮抗肌（antagonist）或中性（neutral）。(4) 肌肉如何通过"募集"和"频率调节"两个机制共同决定总力输出？(5) 肌肉的哪些形态学参数（PCSA、纤维长度、pennation angle）决定了它的功能特性？快肌纤维（type II）和慢肌纤维（type I）的力学差异是什么？(6) 肌肉如何适应训练？增加肌力、增加耐力、增加爆发力的训练对肌纤维类型和形态的影响有何不同？(7) 肌肉力学性质在不同年龄、不同病理状态（痉挛、瘫痪、萎缩）下如何变化？

研究动机：理解肌肉力学是理解所有主动运动的钥匙——从坐到站、走路跑步举重、到精细的手指运动都涉及肌肉收缩的力学。临床上肌肉拉伤、肌腱炎、肌少症（sarcopenia）、脑瘫（CP）的肌张力异常、外科术后肌肉功能恢复，都需要这套肌肉力学知识。

主要公式与推导

肌力与肌节长度的关系由 Gordon 等的 length-tension curve 描述：当肌节长度在 2.0-2.2 μm（粗细肌丝重叠最适区）时主动力最大；超出该范围时（细肌丝在中线重叠或粗细肌丝脱离）主动力下降至零。粗肌丝长度约 1.6 μm、细肌丝长度约 1.05 μm。整条肌肉的长度-张力关系是单肌节曲线在肌纤维上沿纤维长度的积分。

Huxley 1957 的肌丝滑动模型导出 force-velocity 关系

\[ (F + a)(v + b) = b(F_0 + a) \]

其中 \(F_0\) 是最大等长力，\(v\) 是肌节缩短速度，\(a\)、\(b\) 是常数（Hill 特征常数）。离心收缩时 \(v < 0\)、\(F\) 可以超过 \(F_0\)，最大离心力约 \(1.3-1.8 F_0\)。

肌肉功率（power）

\[ P = F \cdot v \]

最大功率约在 \(v = 0.3 v_{\max}\) 处达到。

肌肉的力矩（moment）

\[ M = F_m \cdot d \]

其中 \(d\) 是 moment arm。肌腱的串联弹性遵循 Hooke 定律

\[ F_t = k_t \cdot \Delta L_t \]

其中 \(k_t\) 是肌腱刚度。肌腱与肌纤维串联后总顺应性为两者之和

\[ \frac{1}{k_{\text{total}}} = \frac{1}{k_f} + \frac{1}{k_t} \]

肌肉生理横截面积（physiological cross-sectional area, PCSA）

\[ \text{PCSA} = \frac{m \cos\theta}{\rho L_f} \]

其中 \(m\) 是肌肉质量，\(\rho\) 是密度（约 1060 kg/m³），\(L_f\) 是纤维长度，\(\theta\) 是羽状角（pennation angle）。最大肌力

\[ F_{\max} = \sigma_{\max} \cdot \text{PCSA} \]

其中 \(\sigma_{\max}\) 是最大应力，约 30-40 N/cm²。

关键算法与建模方法

本章介绍的肌肉力学模型主要是 Hill-type 模型（contractile element + series elastic element + parallel elastic element）。contractile element（CE）用 force-velocity 和 force-length 关系描述主动收缩，series elastic element（SEE）模拟肌腱和横桥的弹性，parallel elastic element（PEE）模拟静息时的被动张力。三者并联-串联组合，构成完整的肌腱-肌肉单元模型。

更复杂的肌骨模型（musculoskeletal model）将每块肌肉简化为一条力发生线（line of action）经过起止点，包裹在关节力矩计算中。典型工作流：(1) 用 MRI 或尸体解剖数据定义每块肌肉的起止点、体积、纤维走向、羽状角；(2) 计算 PCSA 和 \(F_{\max}\)；(3) 根据关节角度和肌肉长度反算力臂 \(d\)；(4) 根据关节角速度反算肌纤维缩短速度；(5) 用 force-velocity 计算实际肌力 \(F\)；(6) 用 \(M = F \cdot d\) 计算肌肉对关节的力矩贡献。

OpenSim 是当前最常用的肌骨建模软件平台（Delp 等 2007），它实现了上述流程并有现成的全身 200+ 肌肉模型。AnyBody Modeling System、LifeModeler、 Anybody Technology 等是商业替代。

募集（recruitment）和频率（rate coding）：Henneman 的"尺寸原则"（size principle）指出小运动神经元（连接慢肌纤维）先募集，大运动神经元（连接快肌纤维）后募集。增加肌力有两种策略：增加募集数量（recruitment）和增加已募集运动神经元的发放频率（rate coding）。低强度收缩主要靠募集，高强度收缩还靠频率。

主要结论

本章的主要结论是：(1) 肌力大小由三个层面决定：单肌节的力-长度-速度性质、整块肌肉的 PCSA、神经系统的募集与频率调节。(2) 同一块肌肉在不同关节角度下产生不同的力矩（moment），因为 moment arm 随角度变化；这意味着"肌肉在某个角度最强"这种说法只有在明确角度和运动任务下才有意义。(3) 离心收缩（eccentric / lengthening）的最大肌力显著高于向心收缩（concentric / shortening）——这是肌肉能在"刹车"运动中产生大力的机制，也是肌腱炎和延迟性肌肉酸痛（DOMS）常发生于离心运动的原因。(4) 肌肉的形态学参数（PCSA、纤维长度、pennation angle）决定了肌肉的"专长"——大 PCSA、短纤维适合产生大力量（如股四头肌）；小 PCSA、长纤维适合产生大速度/位移（如缝匠肌）。(5) 训练会改变肌肉形态——力量训练增加 PCSA、可能使 pennation angle 增大；耐力训练增加线粒体密度、毛细血管数；爆發力训练两者兼具。(6) 废用、年龄、疾病会导致肌肉萎缩、肌力下降、肌纤维类型转化——这是康复医学的重要依据。

挑战与开放问题

本章的挑战集中在肌力的体内测量和肌肉-肌腱单元动力学的精细建模上。体内肌力不能直接测量（只能通过肌电图 EMG 间接推断或通过肌骨模型反算），这使得肌肉力学的验证依赖循环可逆的力矩推断。肌腱的弹性储存-释放（stretch-shortening cycle, SSC）在快跑、跳跃中起重要作用，但其精确贡献量仍未完全定量。肌肉的"协同"（synergy）问题——多个肌肉如何在功能上分组以简化控制——是当前神经肌肉控制研究的活跃领域。肌少症（sarcopenia）的机制、肌肉干细胞（satellite cell）的活化、年龄相关的肌纤维类型转化，都是临床-力学交叉的开放问题。

个人反思与批判性分析

本章是 Ch4——全书承上启下的核心章。Oatis 在这里把肌肉力学讲得既经典又清晰：Hill-type 模型、PCSA、size principle 都是经典。Ch4 与 Ch3 骨力学一起，构成了"组织-材料"层面生物力学的两个基础。

但有两个批评：(1) Oatis 的"肌肉力 = PCSA × 最大应力"是 1970-80 年代的简化，实际上肌肉的力-长度-速度-激活-疲劳四维空间比这复杂得多。(2) Oatis 没有充分讨论"肌肉协同（muscle synergy）"和"中枢控制下的冗余度（redundancy）"问题——人体有 600+ 肌肉、200+ 自由度，但运动任务只有几十个——为什么有这么多冗余？这是现代运动控制研究的核心问题（Bernstein 的"自由度问题"）。

实际研究含义：Hill-type 模型配合 OpenSim 仍是最实用的肌骨建模工具，但用它预测"个体在特定任务下的精确肌力"误差仍大（30-50%）。要更精确，需要：(1) 个体特异性 MRI 测量每块肌肉的体积和走向；(2) 动态 EMG 验证；(3) 力矩实测校验。这是当前运动医学和临床步态分析的前沿。

重要参考文献

[1] Huxley AF. Muscle structure and theories of contraction. Prog Biophys Biophys Chem. 1957;7:255-318.（肌丝滑动的奠基论文） [2] Hill AV. The heat of shortening and the dynamic constants of muscle. Proc R Soc Lond B. 1938;126(843):136-195.（Hill 力-速度方程） [3] Gordon AM, Huxley AF, Julian FJ. The variation in isometric tension with sarcomere length in vertebrate muscle fibres. J Physiol. 1966;184(1):170-192.（肌节长度-张力关系） [4] Delp SL, Anderson FC, Arnold AS, Loan P, Habib A, John CT, Guendelman E, Thelen DG. OpenSim: open-source software to create and analyze dynamic simulations of movement. IEEE Trans Biomed Eng. 2007;54(11):1940-1950.（OpenSim 软件平台） [5] Henneman E, Somjen G, Carpenter DO. Functional significance of cell size in spinal motoneurons. J Neurophysiol. 1965;28(3):560-580.（size principle） [6] Lieber RL. Skeletal Muscle Structure, Function, and Plasticity. 3rd ed. Lippincott.（肌肉结构-功能-适应的标准教材） [7] Fukunaga T, Roy RR, Shellock FG, Hodgson JA, Day MK, Lee PL, Kwong-Fu H, Edgerton VR. Physiological cross-sectional area of human leg muscles based on magnetic resonance imaging. J Orthop Res. 1992;10(6):926-934.（MRI 测 PCSA 的标准方法） [8] Zajac FE. Muscle and tendon: properties, models, scaling, and application to biomechanics and motor control. Crit Rev Biomed Eng. 1989;17(4):359-411.（Hill-type 模型的工程化）