坡度对跑步时下下肢双侧不对称性与动态稳定性的影响

《Frontiers in Physiology》：Effects of slope on lower-limb bilateral asymmetry and dynamic stability during running

　　

目的：在固定速度下检验不同跑步坡度对下肢双侧不对称性及动态稳定性的差异。方法：对公开多坡度跑步生物力学数据集进行重复测量二次分析。19名健康跑者以2.78 m·s-1在-6°至+6°五种坡度跑步。采用百分比差值法计算5个峰值负荷变量的不对称性。稳定裕度（Margin of Stability, MoS）分析中，前后向骨盆速度通过加上2.78 m·s-1跑带速度表示为相对跑带基准；内外侧向计算无需跑带速度校正。坡度效应采用线性混合效应模型分析，并使用Benjamini–Hochberg方法控制错误发现率。结果：跟腱应变、膝关节伸展力矩、髌股关节应力的不对称性以及前后向和内外侧向MoS均存在显著坡度效应（q < 0.05），但垂直地面反作用力与踝跖屈力矩不对称性无显著效应。前后向MoS从-6°的-0.573 ± 0.021 m单调增加至+6°的-0.523 ± 0.033 m（F4,72 = 77.74，q < 0.001，偏η2 = 0.81）；内外侧向MoS从0.062 ± 0.022 m降至0.052 ± 0.031 m（F4,72 = 3.50，q = 0.016，偏η2 = 0.16）。跟腱应变不对称性从下坡到上坡减小，而膝伸展力矩和髌股应力不对称性在轻度上坡最大。结论：下肢双侧负荷分布与方向特异性机械裕度随坡度变化；负荷大小与不对称性不同步变化，说明不对称性可提供绝对负荷之外的补充信息。跑步中强负性前后向MoS在力学上合理，不应单独解释为整体不稳定。

研究背景方面，跑步是全球普及的运动，越野跑等地形运动参与增加，但下肢损伤负担集中于膝、踝等部位。坡度会系统性改变下肢负荷分布：上坡更需要正向做功，踝与跟腱负荷增加；下坡以耗能为主，膝伸展肌离心负荷增加，与过度使用损伤相关。已有研究多关注水平跑速度对双侧不对称性（Bilateral Asymmetry）的影响，或关注坡度对整体关节力学和外部负荷的改变，较少同时考察左右肢负荷分配与动态稳定性（Dynamic Stability）。动态稳定性常用外推质心（XCoM）与稳定裕度（Margin of Stability, MoS）框架评价，但既有坡度研究多针对步行或采用最大李雅普诺夫指数评价局部动态稳定性，未量化XCoM与支撑基底面之间的空间稳定裕度。因此，研究人员在固定速度下利用公开数据集，分析五类峰值负荷变量的双侧不对称性和前后向、内外侧向MoS如何随坡度变化。

该研究发表于《Frontiers in Physiology》。研究人员对Van Hooren and Meijer, 2024等发布的公开跑步生物力学数据集做二次分析，样本为19名健康跑者，男女各约半数，水平从非训练者到奥运级，在2.78 m·s^-1下完成-6°、-3°、0°、+3°、+6°五种坡度跑步。结论上，坡度以部位特异性方式改变下肢双侧负荷分配与方向特异性机械裕度；负荷幅值和不对称性不一定同向变化；前后向MoS虽强负但在跑步中力学合理，不能单独作为整体不稳定证据。

关键技术方法上，研究人员使用公开数据集，样本来自19名健康跑者，在仪器化双跑带跑步机上以固定速度跑步。动作捕捉与测力数据经OpenSim和改良Catelli全身肌骨模型处理，通过逆向运动学、逆向动力学和动态优化获得关节角、净关节力矩、肌肉力及跟腱应变、髌股关节应力等组织负荷。支撑相由原始事件文件确定，并按峰值时序做质量控制。双侧不对称性采用百分比差值法，动态稳定性采用XCoM-MoS框架，其中前后向速度需加跑带速度2.78 m·s^-1。统计采用含受试者随机截距的线性混合效应模型，并对7个主要结局用Benjamini–Hochberg控制错误发现率，显著后做Holm校正两两比较。

研究结果部分，摘要结果显示，跟腱应变、膝关节伸展力矩、髌股关节应力不对称性及前后向与内外侧向MoS有显著坡度效应，垂直地面反作用力（vGRF）和踝跖屈力矩不对称性无显著效应。跟腱应变不对称性从下坡到上坡下降，膝伸展力矩和髌股应力不对称性在轻度上坡较大。

2 Methods部分说明，研究人员用公开数据集做重复测量二次分析，固定速度仅比较坡度；参与者为19名健康跑者；原始数据由计算机辅助康复环境采集，标记点经Woltring五次样条平滑后导入OpenSim 3.3与改良Catelli模型；组织负荷由数据库作者用已发布模型计算；支撑相取自有事件文件，质量控制保留完整且数据流一致的支撑相；不对称性用Higher减Lower除以Higher再乘100；MoS以骨盆标志几何中心近似全身体质心，XCoM=p+v/ω₀，ω₀=√(g/h)，前后向需加跑带速度；统计用lme4拟合线性混合效应模型，lmerTest做Satterthwaite自由度，七个主要检验用Benjamini–Hochberg校正，事后比较用emmeans与Holm法。

2.5至2.7结果显示，经质控后13152个支撑相进入vGRF、关节力矩和动态稳定性计算，12792个进入组织负荷计算。原始左右肢峰值显示，踝跖屈力矩和跟腱应变多右侧更大，膝伸展力矩和髌股应力多左侧更大，vGRF方向较均匀。

Asymmetry index结果显示，跟腱应变不对称性从-6°的10.25%降至+6°的5.91%，F_4,72=6.09，q<0.001；-6°显著大于0°、+3°、+6°。膝伸展力矩不对称性在+3°最大，F_4,72=4.21，q=0.009；0°和+3°均显著大于-6°。髌股应力不对称性F_4,72=3.62，q=0.016，但两两比较未全显著。vGRF不对称性q=0.121，踝跖屈力矩不对称性q=0.064，均无主要显著性。

Dynamic stability结果显示，前后向MoS从-6°的-0.573±0.021 m增至+6°的-0.523±0.033 m，F_4,72=77.74，q<0.001，偏η²=0.81，所有两两比较显著；内外侧向MoS从0.062±0.022 m降至0.052±0.031 m，F_4,72=3.50，q=0.016，仅-6°与+6°差异显著。负性前后向MoS在跑步中力学合理。

讨论部分总结，坡度造成下肢负荷在踝、膝、髋间重新分配，但本研究进一步表明这种再分配具有部位特异性和个体方向异性。跟腱应变随上坡增大而双侧不对称性减小，说明高负荷可能压缩左右差异；膝伸展力矩则负荷最大在下坡，而不对称性最大在轻度上坡，说明负荷幅值与不对称性不平行。双侧不对称性反映个体化和结果特异的神经肌肉控制，不能只用绝对负荷判断单侧暴露风险。MoS方面，前后向强负值来自跑带速度校正后的XCoM框架，不代表整体失稳；前后向MoS变小负、内外侧向MoS略降，与同一数据集局部动态稳定性结果互补但不等价。研究人员也指出局限：单速度、样本异质、无解剖不对称信息、骨盆中心近似COM、跑带标称速度假设、方向性仅描述性分析。

结论部分翻译如下：在固定跑步速度下，下肢双侧负荷分布与方向特异性机械裕度随坡度变化。跟腱应变、膝关节伸展力矩、髌股关节应力的不对称性以及前后向和内外侧向MoS均表现出跨坡度差异，而垂直地面反作用力与踝跖屈力矩不对称性未检测到统计学显著差别。将前后向速度表示为相对跑带运动后，前后向MoS呈强负值并从下坡到上坡逐渐变得不那么负，而内外侧向MoS下降较温和。不对称性幅值并未与绝对负荷同步变化，支持在绝对负荷之外同时考虑双侧负荷分布。但校正后的MoS值代表方向与事件特异的瞬时机械裕度，并非整体运动稳定性的直接指标。由于未评估训练适应、扰动反应、表现和损伤结局，其实用意义仍属初步。