M22 · 关节功率与能量流（Joint Power and Energy Flow） 阶段五 · 动力学与能量学 本章回答一个核心问题：肌肉产生的机械能去了哪里？ Winter 在章首直言：关节功率 Pm = Mj·ωj 是所有运动学/动力学变量中对病理步态鉴别力最强的一个。学完本章，应能为任意环节列出完整功率平衡方程，并解释”能量如何在环节之间被动传递、主动转移、被生成与被吸收”。
0. 大白话定义
把人体想象成一串用铰链（关节）连起来的钟摆，肌肉是装在铰链上的小马达。
- **能量（Energy）**是”存款”：某个瞬间身体存了多少机械能（位置高 = 势能存款多，跑得快 = 动能存款多）。
- **功（Work）**是”转账总额”：一段时间内马达往身体里充了多少钱（正功）或从身体里抽走多少钱（负功）。
- **功率（Power）**是”转账速度”：此刻马达充钱/抽钱有多快。
最妙的一条规则：钱可以在环节之间”过路”而不花掉。走路摆动末期，小腿的能量会”借道”大腿上传到躯干再利用——这不是马达发力，只是关节力把能量推过去的（被动传递，全身代数和为零）。而双关节肌还能扮演”搬运工”：主动把能量从一个环节搬到相邻环节。
病理步态之所以在功率曲线上”现形”，是因为异常往往不在”存了多少钱”，而在”马达何时充钱、何时抽钱、充抽是否错位”。
1. Definition 术语表（中英对照）
| 英文术语 | 中文 | 一句话定义 |
|---|---|---|
| Mechanical Energy | 机械能 | 系统某一时刻的状态量，单位焦耳（J），分势能与两种动能 |
| Work | 功 | 能量随时间的流量，是过程量而非状态量 |
| Muscle Power | 肌肉（关节）功率 | 净关节力矩与关节角速度之积，W |
| Positive / Negative Work | 正功 / 负功 | 力矩与角速度同向（向心收缩）做正功；反向（离心收缩）做负功 |
| Internal / External Work | 内功 / 外功 | 内功移动自身体节；外功作用于外部负荷（含举高自身体重） |
| Passive Energy Transfer | 被动能量传递 | 关节反作用力 × 关节点速度产生的跨关节能量流，全身代数和为零 |
| Active Energy Transfer | 主动能量传递 | 肌肉在相邻环节同向转动时把能量从一环节搬到另一环节 |
| Energy Conservation (within segment) | 环节内能量守恒 | 势能与动能在环节内部相互转化、总能量不变（理想钟摆为 100%） |
| Segment Power Balance | 环节功率平衡 | 环节能量变化率 = 近远端关节被动流 + 近远端肌肉主动流 之和 |
| Mechanical / Work Efficiency | 机械效率 / 功效率 | 机械功与代谢消耗之比的不同口径定义 |
| Cocontraction | 共收缩 | 主动肌与拮抗肌同时发力，能量一生一吸的隐形损耗 |
2. 核心框架与公式
2.1 三种能量储存形式（式 6.15–6.18）
- 势能：PE = mgh（h 为质心相对参考基准的高度，基准应取该运动中身体到达的最低点）
- 平动动能：KEₜ = ½mv²（v 为质心速度）
- 转动动能：KEᵣ = ½Iω²（I 为转动惯量）
- 环节总能量：Es = mgh + ½mv² + ½Iω²
注意动能与速度平方成正比，方向无关，静止时为零。
2.2 肌肉功率与功（式 6.1–6.7）
- Pm = Mj·ωj——本章最重要的一个乘积。Mj 与 ωj 同号为正功（生成），异号为负功（吸收）。
- 实际运动中 Mj 与 ωj 常错开约 90° 相位 → 功率曲线正负交替（如前臂屈伸一周两循环）。
- 功 = 功率对时间的积分；要分别得到总正功与总负功，必须精确定位功率曲线过零的时刻，分段积分。
- 一般力做功：P = F·V cosθ = Fx·Vx + Fy·Vy。
2.3 环节间被动能量传递
- 关节点处的能量流 = 关节反作用力与关节点速度的点积：Pj = Fx·Vx + Fy·Vy。
- 对每个环节各算一次，全身代数和恒为零——被动传递不创造也不消灭能量，只是”搬运”。
- 例证：步行摆动末期，足/小腿的能量经大腿被动上传至躯干再利用（Winter & Robertson, 1978）。
2.4 环节内能量交换与守恒（式 6.19–6.23）
- 交换标志：势能与动能反相变化（图 6.12：步行中 HAT 的高度与前向速度几乎完全反相，类钟摆）。
- 近似式：Eex = ΔEp + ΔEkt + ΔEkr − ΔEs。完全无交换时三分量之和 = ΔEs；100% 交换时 ΔEs = 0。
- 精确式（多峰多谷时，按绝对变化累加）：
- Ws = Σ|ΔEs|（实际对环节做的功）
- Ws′ = Σ(|ΔEp| + |ΔEkt| + |ΔEkr|)（假设无交换时需要的功）
- 守恒量 Wc = Ws′ − Ws；守恒百分比 Cs = Wc / Ws′ × 100%
- 算例 6.4：步行中小腿 ΔEs = 16.16 J，三分量之和 = 16.65 J → 仅 0.49 J 被交换，小腿是高度非守恒系统（与 HAT 形成鲜明对比）。
2.5 全身总能量与内功计算的四代方法（式 6.24–6.26，Table 6.1）
全身总能量：Eb = Σ Esi（逐时刻对全部环节求和，Bresler & Berry 1951；Ralston & Lukin 1969；Winter et al. 1976）。全身内功：Wb = Σ|ΔEb|（式 6.25）。
| 方法 | 遗漏的成分 |
|---|---|
| ① 环节能量增量法（Fenn 1929；Saunders 1953） | 环节内交换、环节间传递 → 高估（算出短跑运动员 3 马力！） |
| ② 质心法（Cavagna & Margaria） | 把质心当全身代表；但能量是标量，反向运动环节的同增同减被”对消”而漏算 → 系统性低估（M20 留下的”COM 能量误用”正解在此） |
| ③ 环节能量求和法（Winter 1979） | 承认环节内守恒与环节间传递，但漏算不同关节同时生成/吸收的对消，仍低估 |
| ④ 关节功率积分法 ∫Mj·ωj dt（Winter 1983） | 只漏共收缩与抗重力等长收缩；若有外功会自动计入 |
| ⑤ 肌肉功率法 ∫Fm·Vm dt（Morrison 1970；Yack 1986） | 理论上最全，但需要把净力矩拆分到每块肌肉（难点见第 5 章） |
2.6 关节功率曲线解读：慢跑膝关节五相位（图 6.16，Winter 1983）
- K1：支撑早期，膝伸肌吸收 −53 J（离心缓冲）
- K2：支撑中后期，膝伸肌生成 +31 J（向心蹬伸）
- K3：支撑末/摆动初，小腿后旋减速 −11 J
- K4：摆动末期，膝屈肌吸收 −24 J（减速前摆小腿）
- K5：触地前小幅正向爆发 +5 J
对该跑者：吸收（53 + 24 J）远大于生成（31 J）——膝关节在慢跑中主要是”刹车”而非”发动机”。
2.7 共收缩的量化（COCON，式 6.12–6.14）
- 力矩版、EMG 共同面积版（= 2 × 共同面积 / (A + B) × 100%；步态胫前肌-比目鱼肌算例 24%）、Rudolph (2000) 版三套公式。
- ⚠️ Kellis (2003) 警告：公式选择会剧烈改变数值，跨研究比较前必须核对定义。
- 共收缩造成的隐形损耗：净力矩 40 N·m × 3 rad/s 算得 120 W；若拮抗肌暗中抵抗 10 N·m，则主动肌实际以 50 × 3 = 150 W 生成、拮抗肌以 10 × 3 = 30 W 吸收——净功率法两边都低估。
2.8 肌肉介导的主动能量传递（式 6.28）
两环节同向转动时，跨关节肌肉可同时”一头吸、一头放”，即搬运能量：
Pm = Mj·(ω₁ − ω₂)
- ω₁ = ω₂（等长）：纯传递
- ω₁ > ω₂（肌肉拉长）：吸收 + 传递
- ω₁ < ω₂（肌肉缩短）：生成 + 传递
图 6.18（Robertson & Winter, 1980）总结了关节处全部可能的功率功能组合（生成/吸收/传递）。
2.9 环节功率平衡方程（式 6.29）
对任意环节：
dEs/dt = Pjp + Pmp + Pjd + Pmd
近端关节被动流 + 近端肌肉主动流 + 远端关节被动流 + 远端肌肉主动流 = 环节能量变化率。
算例 6.5（用附录 A 数据，frame 5，Δt = 28.6 ms）：
- 小腿：Σ功率 = 44.81 + 7.19 − 33.3 − 1.36 = 17.34 W；ΔElg/Δt = (20.5 − 20.0)/0.0286 = 17.5 W → 平衡残差仅 0.16 W ✓
- 大腿：Σ功率 = 14.97 + 33.83 − 44.81 − 23.08 = −19.09 W；ΔEth/Δt = (47.4 − 47.9)/0.0286 = −17.5 W → 残差 1.59 W ✓
- 结论：23.08 W 从大腿流入膝伸肌，7.19 W 从膝伸肌流入小腿 → 膝伸肌正以 7.19 W 把能量从大腿主动搬向小腿，同时吸收 15.89 W；另有 44.81 W 经膝关节被动流过。
2.10 机械效率（式 6.8–6.11）
- 代谢（肌肉）效率：不可直接算（无法测单块肌肉代谢）。
- 机械效率 = （内功 + 外功）/（代谢 − 静息代谢）
- 功效率 = 外功 /（代谢 − 零功负荷时代谢）
- 正功与负功的代谢成本不同（Abbot 1952 双人联动自行车实验：正功骑行者额外承担两人内功 150 W + 对方负功），需分别引入 η⁺、η⁻。
- 神经效率陷阱：痉挛型脑瘫患者做更大机械功，按公式反而算出 40% > 健康人 33%——效率分数只反映代谢转换，不反映神经控制优劣。
- 四大低效机制（Winter 1978）：① 共收缩；② 抗重力等长收缩（蹲伏步态股四头肌持续负荷，无标准量化法）；③ 一关节产能另一关节同时耗能（注意：步行双支撑期的这类对消是必要功能，非病态）；④ 急促不连贯运动（能量反复充放）。
3. Physical Meaning（物理意义）
- 能量守恒定律在人体运动分析中的完整落地：能量不生不灭，只在”生成（肌肉正功）—储存（PE/KE）—传递（被动/主动）—吸收（肌肉负功）“之间流动。
- 功率平衡方程 dEs/dt = ΣP 本质上是对单个环节写的能量连续性方程，残差大小直接检验运动学/动力学数据的内部一致性（算例 6.5 残差 < 2 W，说明数据自洽）。
- 动能与 v² 成正比：速度翻倍，动能四倍——这解释了为什么”跑快一点”代谢成本急剧上升。
4. Human Movement Meaning（人体运动意义）
- 正常步行是精心设计的节能系统：HAT 内部势-动能转换（钟摆机制）、小腿能量经大腿上传躯干再利用、摆动期大部分靠被动动力学。
- 但下肢本身高度非守恒（小腿交换率仅 ~3%）：每一步都需要肌肉实际充能/抽能，主要”发动机”是支撑末期踝跖屈肌，主要”刹车”是支撑早期膝伸肌与摆动末期膝屈肌。
- 匀速周期运动中，每步净能量变化为零 → 每步内正功总量 = 负功总量。这条恒等式是检验功率分析正确性的金标准。
5. How It Is Measured（如何测量）
- 运动捕捉 → 环节运动学（位置、速度、角速度）；
- 结合体节参数（质量、质心、转动惯量）→ 按式 6.18 算每个环节每个时刻的 Es；
- 测力台 + 逆向动力学（M21）→ 净关节力矩 Mj；
- Pm = Mj·ωj → 功率曲线过零点分段积分 → 正/负功；
- 关节点速度 × 关节反作用力 → 被动传递流；按式 6.29 做逐环节功率平衡并检查残差；
- 全身功 Wb = Σ|ΔEb|（式 6.25）。
- 共收缩需同步 EMG（M24 展开）。
- 外功只能在人-负荷界面测量（式 6.27：We = ∫Fr·Vc dt，如自行车踏板上的力传感器）。
6. Typical Applications（典型应用）
- 病理步态鉴别：Winter 章首明言 Pm = Mj·ωj 是所有运动学/动力学变量中鉴别力最强的变量——因为力矩、角速度各自的异常会相互补偿掩盖，而功率直接暴露”马达何时充抽钱”。
- 假肢与矫形器评估：储能脚能否替代踝跖屈肌的正功爆发，直接看 A2 相位功率曲线的恢复程度。
- 跑步经济性：比较不同跑姿/鞋/速度下的正负功结构与环节守恒百分比 Cs。
- 能量流审计：定位”能量在哪一刻、哪个关节、被谁生成/吸收/搬运”（算例 6.5 是标准范式）。
- 效率研究：运动处方中比较不同工种的机械效率，但必须警惕神经效率陷阱。
7. Clinical Relevance（康复临床相关性）
- 脑瘫蹲伏步态：股四头肌持续等长抗重力收缩——功率曲线看不到（无关节运动则无功率），但代谢成本巨大。解释：等长收缩耗能不做功，属四大低效机制中”无标准量化法”的一项，需结合 EMG 与耗氧量评估。
- 共收缩（如卒中后、ACL 损伤后）：COCON 指数量化；注意 Rudolph (2000) 用其研究 ACL 损伤后的动态稳定性——共收缩在此可能是代偿性保护策略而非纯损耗，干预目标不一定是”消除”。
- 能量传递中断：股直肌（双关节肌）本应把髋关节能量搬向膝关节；选择性运动控制障碍患者搬运功能受损 → 步态启动摆动期费力。[补充说明：这是”直腿抬高迟滞”等临床现象的能量学解释框架。]
- 效率解读红线：绝不把”机械效率更高”直接等同于”运动控制更好”（脑瘫 40% vs 健康人 33% 的反例）。
8. Performance Relevance（运动表现相关性）
- 短跑运动员的机械功率量级远超日常（Fenn 曾高估出 3 马力，修正方法后仍远超步行），爆发力项目的核心就是单位时间充能速率。
- 投掷项目中鞭打效应 = 能量沿近端→远端环节逐级被动传递 + 主动叠加（式 6.28 机制）。
- 耐力项目追求：提高环节内守恒百分比 Cs、最大化双关节肌的搬运功能、最小化共收缩。
- 跳高/跳远起跳：单关节功率峰值与正功总量是选材与训练监控的核心指标。
9. Possible Interpretations（表现型 → 机制推理）
| 观察（表现型信号） | B 类：测量/技术因素 | C 类：功能性因素 | D 类：结构性/病理性因素 |
|---|---|---|---|
| 踝 A2 正功爆发减小 | 测力台与镜头不同步；关节中心定位偏差；滤波过度削峰 | 步速减慢（功率对速度高度敏感）；疲劳 | 小腿三头肌无力（卒中后足下垂）；跟腱病变；储能脚性能不足 |
| 膝 K1 吸收增大 | 体节参数错误致力矩偏大 | 步速加快；下坡走；落地策略僵硬 | 膝伸肌离心控制差；前交叉韧带术后代偿 |
| 全身正功 ≠ 负功（周期运动） | 功率过零点定位不准；分段步长太大 | 加减速阶段（非稳态）；坡道 | ——（稳态周期运动中此恒等式必须成立，不成立先查数据） |
| 双关节肌传递流缺失 | 模型未做功率平衡分解（只算净功率） | 步速过慢使 ω₁≈ω₂ 窗口消失 | 双关节肌选择性控制障碍；挛缩致跨关节长度异常 |
| COCON 指数升高 | EMG 串扰；归一化方法不同（Kellis 警告） | 新任务学习期；不稳定支撑面；疼痛性保护 | 痉挛；锥体外系疾病；关节不稳的神经代偿 |
10. What It Cannot Tell Us（局限性）
- 净力矩掩盖共收缩：Mj·ωj 看不到主动肌与拮抗肌的同时充抽（6.3.1.5 节算例：120 W 净功率背后是 150 W 生成 + 30 W 吸收），须加 EMG。
- 等长收缩隐形：无关节角速度 → 无机械功率，但代谢成本真实存在（蹲伏步态的股四头肌）。
- 单关节功率无法区分”生成”与”传递”：必须做式 6.29 的功率平衡分解。
- 无法拆分到单块肌肉：∫Fm·Vm dt 需要肌肉力分配模型，目前手段有限。
- 质心法不可作为替代：能量是标量，反向运动环节的同增同减会被对消漏算（M20 悬案的正解）。
- 效率分数不含神经控制信息：高机械效率 ≠ 好的运动模式。
11. Related Metrics（相关概念双链）
- 上游输入：M21_力台与链段模型（Mj 的来源）、M21b_COM与COP（全身层面参照）、M19b_平滑与微分问题（滤波与微分误差对功率曲线的影响）
- 下游应用：M24_肌电信号与肌肉激活（共收缩量化、肌肉层面分解）、M25_肌骨建模与平衡控制（支撑力矩、倒立摆能量）
- 综合案例：M26_步态综合案例
- 论文延伸：M27_Bahr2005_损伤机制（落地冲击的能量吸收机制）
12. Evidence（证据与出处）
- 环节能量公式与四代内功计算方法比较：Winter 第 6 章 6.2–6.3 节，Table 6.1（书页 128–137）。
- 慢跑膝关节功率五相位：Winter (1983), J. Biomech. 16:91–97，图 6.16（书页 134–135）。
- 被动传递与主动搬运机制：Winter & Robertson (1978), Biol. Cybern. 29:137–142；Robertson & Winter (1980), J. Biomech. 13:845–854。
- 功率平衡算例（算例 6.5，图 6.19–6.20）：书页 138–141，数据取自教材附录 A。
- 负功代谢成本：Abbot, Bigland & Ritchie (1952), J. Physiol. 117:380–390。
- 共收缩量化与公式敏感性：Falconer & Winter (1985)；Rudolph et al. (2000)；Kellis et al. (2003), J. Electromyogr. Kinesiol. 13:229–238。
- 病理步态能量评估与脑瘫效率反例：Winter (1978), Physiother. Canada 30:183–191。
13. My Understanding（我的理解，5 点）
- 功率 = 力矩 × 角速度，一乘定乾坤。力矩大但关节不动（等长）= 零功率；关节动但力矩小 = 小功率。只有两个都”在线”且方向一致，肌肉才真正在充能。这是为什么功率曲线比单独的力矩曲线、角度曲线更有鉴别力。
- 能量的”三权分立”：生成/吸收归肌肉（主动），搬运归关节反作用力（被动，全身代数和为零）和双关节肌（主动传递，式 6.28）。看步态曲线时脑子里要有这三本账。
- COM 能量法的死刑判决书：能量是标量不是向量——这是质心法低估全身功的根本原因，也是 M20 遗留问题的正解。全身功必须逐环节算 Es 再求和。
- 功率平衡残差是数据质量的试金石：算例 6.5 残差 0.16 W / 17.5 W，说明运动学+动力学+体节参数自洽；残差大说明哪里测错了。这套审计思路可直接移植到我的系统做数据质控模块。
- “低效”不一定该治：双支撑期一腿产能一腿耗能是步行的必要功能；ACL 术后共收缩可能是保护性代偿。机制推理系统必须区分”病态低效”与”功能性代价”——这是表现型→机制推理的关键分界。
模块自测（7 题）
Q1 某棒球质量 1 kg，在离地 2 m 处以 18 m/s 竖直上抛。求出手瞬间的总机械能；到达最高点时势能为多少？（g = 9.8 m/s²，忽略空气阻力）
参考答案：出手瞬间 PE = mgh = 1 × 9.8 × 2 = 19.6 J；KE = ½mv² = ½ × 1 × 18² = 162 J（注意 162 J 恰等于投掷过程对球做的功，见算例 6.1）；总能量 = 19.6 + 162 = 181.6 J。最高点动能为零，全部转化为势能，故 PE(max) = 181.6 J，对应高度 h = 181.6/9.8 ≈ 18.5 m。
Q2 净关节力矩 Mj = 40 N·m，关节角速度 ωj = 3 rad/s，两者同向。求关节功率。若此时拮抗肌暗中产生 10 N·m 反向力矩，主动肌实际的生成功率与拮抗肌的吸收功率各为多少？净功率法低估了多少？
参考答案：净功率 = 40 × 3 = 120 W（正功，生成）。存在共收缩时：主动肌力矩 = 40 + 10 = 50 N·m，生成功率 = 50 × 3 = 150 W；拮抗肌吸收功率 = 10 × 3 = 30 W。净功率法把正功低估了 30 W（120 vs 150），负功则完全漏算（0 vs 30 W）——这正是共收缩使净功率法双侧低估的机制。
Q3 步行中测得小腿各环节能量极值：Es,max = 29.30 J，Es,min = 13.14 J；ΔEp = 2.16 J，ΔEkt = 13.54 J，ΔEkr = 0.95 J。用近似式计算该环节的交换能量，并判断它是否接近理想钟摆。
参考答案：ΔEs = 29.30 − 13.14 = 16.16 J；三分量之和 = 2.16 + 13.54 + 0.95 = 16.65 J；交换量 Eex = 16.65 − 16.16 = 0.49 J，仅占 16.65 J 的约 3%。理想钟摆应 100% 交换（ΔEs = 0），故小腿是高度非守恒系统——步行中下肢的能量基本靠肌肉实际充抽，而非内部转化。
Q4 为什么用”身体质心的势能与动能变化”计算全身内功会系统性低估？请说出两条原因。
参考答案：① 能量是标量而质心是矢量和：步行中大量反向运动的环节（如左右腿）其加速度在质心处对消，但各自的动能变化是正的标量、不能对消——质心看不到”环节各自忙”，同增同减被漏算；② 质心法不含转动动能，也漏掉不同关节同时生成与吸收的对消损耗。此外它绑定测力台数据，跑步腾空期完全无法计算（Winter, 1979 已系统论证）。
Q5 稳态匀速跑步中，每步的净机械能变化约等于多少？由此可推出每步总正功与总负功之间必须满足什么关系？若你的分析结果不满足该关系，第一步该检查什么？
参考答案：净能量变化 ≈ 0（忽略空气与鞋的微小摩擦损耗），故每步内总正功 = 总负功。若不满足，先查技术因素：功率曲线过零点定位是否精确（积分分段错误）、滤波是否过度削峰、运动学与力数据是否同步、体节参数是否正确——可用单环节功率平衡残差（式 6.29）定位是哪一侧数据出问题。
Q6 某时刻大腿角速度 ω₁ = 3.98 rad/s，小腿角速度 ω₂ = 1.24 rad/s（同向），膝伸肌净力矩 M = 5.8 N·m。分别计算肌肉对大腿、对小腿的功率，并说明此刻膝伸肌的三种功率功能。
参考答案：对大腿：P₁ = M·ω₁ = 5.8 × 3.98 ≈ 23.1 W，方向为能量流出大腿（吸收侧）；对小腿：P₂ = M·ω₂ = 5.8 × 1.24 ≈ 7.2 W，能量流入小腿。因 ω₁ > ω₂，肌肉处于拉长状态：它一边从大腿吸收 23.1 W，一边向小腿释放 7.2 W——其中 7.2 W 属于主动传递（搬运），差值 15.9 W 为净吸收（与算例 6.5 的 23.08 / 7.19 / 15.89 W 一致）。
Q7 一位痉挛型脑瘫患者的步态分析显示：其每公斤体重机械功高于同龄健康人，机械效率算出 40%（健康人 33%）。能否据此认为该患者”步态更优”？请说明理由，并指出该患者步态中功率分析可能漏掉的两类真实代谢成本。
参考答案：不能。机械效率只反映代谢能向机械功的转换比例，不反映神经控制质量；患者效率”更高”恰恰是因为其运动被限制在代谢成本相对低的模式里（神经效率陷阱）。功率分析漏掉的两类成本：① 抗重力等长收缩——蹲伏步态中股四头肌持续发力但关节角速度≈0，Mj·ωj ≈ 0，功率曲线上看不到；② 共收缩——主动肌与拮抗肌同时充抽，净力矩法对两者都低估，须用 EMG + COCON 指数补充量化。
上一模块：M21b_COM与COP ｜ 下一模块：M23_3D变量与事件检测