精英与休闲长跑者的多速度生物力学签名：基于可穿戴惯性传感器的1:1匹配配对分析

《Scandinavian Journal of Medicine & Science in Sports》：The Biomechanical Signature of Elite Distance Running: A Matched Pairs Analysis Across Multiple Speeds

　　

本研究使用可穿戴惯性传感器（inertial measurement unit，IMU），在10–15 km/h跑台速度下，通过五个测量域比较精英与休闲长跑者的生物力学特征。共104名跑者（52名精英；52名休闲）在递增跑台方案中使用RunScribe传感器评估。每个速度记录23个变量，涵盖时空参数、刚度和力、冲击负荷、足部运动学和垂直动力学。精英跑者按性别、年龄、身高和体重与休闲跑者1:1匹配。对每个变量，拟合线性混合效应模型（linear mixed-effects model，LMM）于配对内差异，包含配对水平随机截距和速度斜率，并将速度和配对内体重与年龄不平衡作为固定效应。对23个组效应和23个交互检验分别应用错误发现率（false discovery rate，FDR）校正（Benjamini–Hochberg）。FDR校正后，23个变量中有12个（52.2%）显示显著组效应，6个（26.1%）显示显著组×速度交互。精英跑者表现出更短触地时间（Cohen's d_z = ?0.73至?0.53）、更高腾空比（d_z = +0.52至+0.66）以及初始触地时更靠前的足着地模式（d_z = +0.58至+0.67）。RunScribe估算的垂直地面反作用力（ground reaction force，GRF）加载率和腿刚度也高于休闲跑者（分别为d_z = +0.48至+0.75和+0.51至+0.62）。步频和步长未见总体组间差异。显著的组×速度交互显示，对于大多数受影响变量，休闲跑者的数值随速度变化更陡峭。这些发现描述了精英跑者在六个评估的次最大速度下具有独特的力学特征。

（一）研究背景与意义
理解长跑表现的决定因素是运动科学的核心问题。历史上，长跑表现主要通过生理学视角解释，东非跑者的优势使研究关注高原适应、人体测量学特征、饮食习惯和族群特征等。近年来，碳纤维板跑鞋的普及使生物力学边际增益受到关注，跑步经济性和生物力学效率被视为同等最大摄氧量（VO₂max）跑者之间的关键区分因素。但现有精英—休闲跑者生物力学比较存在三方面局限：多数研究样本量小，常为每组n < 10–22；比较变量范围窄，主要集中于时空参数，足部运动学和冲击负荷等域未被充分探索；未充分控制体重、身高和性别等人体测量学混杂因素。因此，真正区分精英与休闲跑者的生物力学签名仍未完全明确。研究人员开展匹配配对、多域、多速度比较，意义在于更严格控制人体测量学混杂，系统识别跨测量域的生物力学特征，并为可穿戴设备监测提供依据。

（二）研究设计与主要结论
研究人员对104名跑者（52对）在五个测量域、23个变量和10–15 km/h六个速度下进行比较，采用性别、年龄、身高和体重1:1匹配的配对设计，并以LMM分析配对内差异。结果显示，FDR校正后12/23个变量（52.2%）存在显著组效应，6/23个（26.1%）存在显著组×速度交互。精英跑者表现为更短触地时间、更高腾空比、更靠前的足着地模式，以及更高的RunScribe估算垂直GRF加载率和腿刚度；步频和步长无总体组间差异。研究人员认为，这些结果描述了精英跑者跨六个次最大速度的独特力学特征，并提示其与更硬、更具弹性的地面相互作用模式一致，但该解释在原文中被定位为假设而非直接证明的机制。

（三）关键技术与方法（≤250字）
样本来自172名跑者（55名精英、117名休闲），经1:1匹配与质控后保留52对（41男、11女；精英10 km个人最好成绩男<33 min、女<36 min）。受试者双侧鞋带佩戴RunScribe（500 Hz），在0%坡度跑台完成6 min热身后，以10–15 km/h固定递增方案运动，每级1 min，取每级后30 s数据，左右肢平均。提取5个域23个变量。对每个变量的配对内差异（精英?休闲）拟合LMM，含配对随机截距与速度随机斜率，速度及配对内体重、年龄不平衡为固定效应；分别对23个组效应和23个交互检验作Benjamini–Hochberg FDR校正，并按速度计算Cohen's d_z。研究经西班牙阿拉贡CEICA伦理批准。

（四）研究结果
3.1 匹配质量：匹配后性别完全平衡，身高残余不平衡可忽略；体重、年龄和体质指数（body mass index，BMI）仍有有意义的配对内残余差异，研究人员将其作为LMM协变量部分校正，但统计调整不能使两组人体测量学完全等价。
3.2 各速度样本可用性：可用匹配对数由10 km/h的52对降至15 km/h的45对；点水平质控后，单一变量最少可用41对（中位数48对）；各速度下人口学特征保持稳定，说明样本量下降但匹配平衡性得以维持。
3.3 生物力学差异概览：23个LMM组效应检验中12个（52.2%）显著，6个（26.1%）存在显著组×速度交互；显著组效应分布于时空参数（5/7）、刚度与力（4/6）、足部运动学（2/4）、垂直动力学（1/2），冲击/冲击域为0/4。
3.4 时空参数：精英跑者触地时间更短（d_z = ?0.73至?0.53）、接触比更低（d_z = ?0.67至?0.52）、腾空比更高（d_z = +0.52至+0.66），步幅角更大（d_z = +0.42至+0.63）；步频和步长无总体组间差异，但步频存在组×速度交互，休闲跑者随速度增加步频上升更陡。
3.5 刚度与用力：精英跑者RunScribe估算垂直GRF加载率更高（d_z = +0.48至+0.75）、峰值垂直GRF更高（d_z = +0.49至+0.66）、腿刚度更高（d_z = +0.51至+0.62）、垂直刚度更高（d_z = +0.26至+0.43）；水平GRF率和总用力率无显著组主效应，但总用力率存在显著组×速度交互。
3.6 冲击负荷与冲击：无冲击变量达到FDR显著的总体组效应；冲击（Shock）和Impact Gs存在显著组×速度交互，提示两组轴向冲击加速度差异随测试速度范围逐渐增大，而非在各速度均匀存在。
3.7 足部运动学：精英跑者初始触地时足着地模式更靠前（d_z = +0.58至+0.67），摆动幅度更低（d_z = ?0.35至?0.26），但后者效应量较小。
3.8 垂直动力学：垂直振荡无总体组效应但存在显著组×速度交互；精英跑者触地时垂直速度更负（d_z = ?0.43至?0.19），且该组效应在BMI校正敏感性分析中未保留。

（五）讨论总结与研究结论翻译
讨论将结果概括为多域“弹簧—质量”签名：更短触地、更长腾空、更高腿刚度、更快垂直GRF加载率和更靠前足着地共同构成协调的生物力学策略，与弹簧—质量模型关于更高效跑者具有更高“弹簧—质量相似性”的预测一致。讨论将足着地差异解释为较固定的技术适应而非速度后果，并指出足部运动学签名不在于“过度旋前”控制，而在于整个支撑期效率。步频和步长总体无差异，提示精英跑者通过重分配步态周期而非增加步数或步长来提速；步频的速度交互支持其在更高速度下较少依赖步频增量。冲击变量总体不显著但存在速度交互，垂直速度更负，解释需谨慎。实践上，接触比和步幅角被视为可穿戴设备可获取的监测指标；讨论还提及增强式训练可提高肌肉—肌腱刚度并改善长跑表现和跑步经济性，但本研究未直接验证干预效果。局限包括：RunScribe足装IMU非金标准测力台或光学动作捕捉、横断面设计不能推断因果、速度范围限于10–15 km/h、绝对速度匹配而非相对生理强度匹配、未收集训练负荷、受试者穿习惯跑鞋，以及女性配对样本（n = 11）未能在FDR校正后保留显著变量。

研究结论翻译（Discussion首段）：本研究提供了首个大规模、经人体测量学匹配、多域的生物力学比较，比较对象为精英与休闲长跑者。通过分析六个速度下五个测量域的23个变量，研究人员识别出一个连贯的精英跑步生物力学签名，其范围远超传统时空参数。出现五项主要发现：第一，精英跑者表现出更短的地面接触时间和更长的腾空时间，无论绝对值还是相对值均如此。第二，精英跑者表现出更靠前的足着地模式。第三，精英跑者表现出独特的用力特征，表现为RunScribe估算的腿刚度和垂直GRF加载率更高。第四，步频和步长未显示显著的组主效应。第五，23个变量中有6个存在显著的组×速度交互，表明在有限的一组评估生物力学参数中，精英?休闲差异的幅度具有速度依赖性。研究人员将精英跑者中这一趋同模式解释为与地面更硬且更具弹性的相互作用一致，但本研究未直接测量代谢成本、弹性回馈和肌肉—肌腱行为，因此该解释应视为假设而非直接证明的机制。