《Journal of Biomechanics》：A cubic Hunt-Crossley contact law and subject-specific calibration improve ground reaction force prediction in musculoskeletal walking simulations
编辑推荐：
预测性肌骨（MSK）步行仿真日益被用于为临床决策和设备设计提供依据，然而其在重现地面反作用力（GRF）方面的准确性仍然有限。一个尚未得到充分审视的 contributory 因素是足-地接触模型及其参数。本研究从三个方面弥补这一空白。第一，研究人员依据在体足跟
预测性肌骨（MSK）步行仿真日益被用于为临床决策和设备设计提供依据，然而其在重现地面反作用力（GRF）方面的准确性仍然有限。一个尚未得到充分审视的 contributory 因素是足-地接触模型及其参数。本研究从三个方面弥补这一空白。第一，研究人员依据在体足跟垫冲击数据重新审视Hunt-Crossley（HC）接触公式，并证明Simbody标准Hertzian弹性律无法刻画足跟垫的双相力学行为：将3/2压痕指数替换为立方指数后，模型对实验力-变形特性的拟合显著改善。第二，研究人员将六个此类改进接触球嵌入OpenSim MSK模型，并通过对六次步态试验的实验接触力旋量进行多目标优化来辨识其力学参数，这一受试者个体化标定步骤在近期基于OpenSim的研究中受到的关注有限。第三，研究人员证明，允许预设的下肢运动学在同时符合运动捕捉测量误差与逆运动学误差的范围内变化，对于恢复准确的GRF至关重要：接触参数与运动学的联合辨识使接触力旋量全部六个分量的归一化均方根误差低于7.1%。在相同流程下使用标准Hunt-Crossley模型重复辨识所产生的误差约为前者的两倍，证实是立方接触律而非仅辨识流程本身带来了改善。留一交叉验证表明所辨识参数在各试验间保持稳定，支持该流程作为受试者个体化的概念验证。
研究背景方面，准确可靠的预测性人体运动仿真有望成为临床医生设计个体化治疗方案、预测治疗效果以及开发骨科与辅助器具的有价值工具，且计算负担正稳步下降。此类工具依赖于具有足够代表性的肌骨（MSK）模型以及生理性神经运动控制策略。尽管足-地接触模型及其参数已被认定为推动预测性步行仿真发展的关键因素，但其作用受到的关注有限，导致对步态关键运动学-动力学特征（如地面反作用力，GRF）的预测不准确。广泛采用的Hunt-Crossley（HC）接触模型参数往往未考虑 underlying 生理结构的力学特性。Pàmies-Vilà等提出的标定流程与本工作思路最为接近，但在与参考接触力的一致性方面仍有改进空间。为此，研究人员提出一种由在体足跟垫冲击数据驱动的新型足-地接触模型，并在先进的3D OpenSim MSK模型中对其参数进行受试者个体化辨识，同时联合进行小幅运动学修正以计入测量与逆运动学不准确性。研究证明，通常采用的Hertzian弹性项无法刻画足跟垫的双相力学行为，并提出立方替代方案；所辨识参数在六次运动试验中对地面反作用力旋量所有分量产生的平均归一化均方根误差（NRMSE）低于7.1%；使用标准HC模型的相同辨识产生的NRMSE约为前者的两倍，留一交叉验证证实基于五次试验标定的力学参数能良好迁移至留出试验。
研究人员开展的研究为：首先依据在体足跟垫冲击数据标定单个接触球，然后通过针对六次实验步态试验的多目标优化辨识OpenSim MSK模型中六个接触球的参数。初步实验揭示所得力对预设运动学高度敏感，促使在符合已报道逆运动学不确定性的范围内对关节轨迹进行联合优化。标准HC模型与留一交叉验证分别作为比较与稳定性检验。
在技术方法上，研究人员采用立方弹性项的HC接触球模型（弹性法向力与压痕立方成正比），摩擦采用连续Stribeck类模型；利用Kinoshita等（1993）提供的在体足跟垫冲击数据（减速度与位移轨迹，初始速度0.93 m/s）进行最小二乘拟合；采用含31个自由度、92块下肢肌肉的3D OpenSim MSK模型，每足布置六个半径32 mm的接触球；实验数据来自一名健康女性受试者（35岁，62 kg，1.70 m）的公开步态数据（六次步态试验）；采用Chebyshev标量化处理六目标优化，通过CasADi与IPOPT求解，并使用支持自动微分的OpenSim版本；运动学修正采用二阶6 Hz零相位Butterworth滤波与中心有限差分；采用留一交叉验证评估参数迁移性。
研究结果方面，第一，仅优化力学参数而直接施加实验运动学轨迹时，GRF、压力中心（COP）与自由矩My预测较差；保留冲击试验参数并额外优化运动学变量后预测显著改善；同时辨识运动学与力学接触参数后获得进一步中等程度改善。第二，对各接触力旋量分量采用NRMSE量化预测精度，最佳配置下所有分量NRMSE低于7.1%，且实验垂直GRF的两个峰值均被良好捕捉，优化后的下肢运动学保持在预设界限内。第三，使用标准HC模型在相同优化配置下预测明显更差，各分量NRMSE约为14.65%–15.68%，证实足跟垫拟合不佳会转化为完整MSK仿真中GRF预测的退化。第四，留一交叉验证显示力学参数在各折间稳定，K的跨折变异性低于11%，表明辨识收敛于稳定解而非单次试验伪影；留出试验的NRMSE在所有分量上均低于标准HC基线。
讨论部分指出，仅优化力学参数时与实验数据一致性差，印证了实验导出运动学的不准确性是预测误差来源的假设；仅运动学精化即带来显著改善，证实运动学准确性是接触力旋量估计质量的关键因素；联合辨识带来进一步改善。本工作总体NRMSE低于7.1%优于Pàmies-Vilà等的报告值，但该比较仅为指示性而非受控比较。标准HC预测明显更差，证实与 underlying 生理力学一致的接触律至关重要。局限性包括：仅使用一个足跟垫冲击数据集推导接触力函数形式的结構假设，需新实验在更多数据与不同受试者上验证立方指数假设；参数辨识按构造为受试者个体化，故在单受试者上展示恰当，而立方函数形式为意图跨受试者推广的结构假设，本工作仅提供概念验证证据；未来工作应研究球特异性参数以捕捉不同解剖足区的非均匀刚度与阻尼，这需要足跟垫以外足区的冲击试验以及区域载荷测量（如压力板或仪器化鞋垫数据），而所用数据集不包含此类数据；采用专用于足部的运动捕捉装置而非全身装置将带来更高运动学精度。
研究结论译文：研究人员提供了（i）一种用于裸足步行的生理学知情接触力模型，以及（ii）一种基于接触球的MSK模型相应参数集辨识方法，同时强调了运动学准确性在估计GRF、COP与自由矩中的重要性。将相同辨识流程应用于标准Hunt-Crossley模型，每个接触力旋量分量的NRMSE约为两倍；因此差距可归因于接触律的选择而非流程本身。该结果同时提升了MSK建模的准确性与GRF预测的可靠性，对预测性肌骨工具以及足部承载矫形器、外骨骼与假体的设计具有潜在意义。留一分析在留出试验上重新辨识运动学修正，与研究人员关于实验运动学携带逆运动学不确定性的立场一致；因此泛化性主张针对力学参数集，其跨折稳定性支持其可迁移性。总之，所提出的接触律与辨识流程将每个力旋量分量的GRF预测误差降至7.1%以内，同时使该结果背后的假设明确且可检验。
