1. 引言
膝关节的损伤是最常见的运动损伤类型 [1] [2],而前交叉韧带(Anterior Cruciate Ligament, ACL)则是最常发生损伤的膝关节韧带 [3]。在这些损伤中,80%属于非接触损伤 [4]。不恰当的跳跃着陆方式是导致非接触ACL损伤的最主要原因之一 [5],而其典型特征包括过小的膝关节屈曲角度,髋关节屈曲角度,躯干屈曲角度,过大的膝关节过度外展角和力矩以及垂直地面反作用力 [6] [7]。
尽管手术技术和康复项目在过去20年里有了显著的改善,但是长期的跟踪调查结果显示,受到过ACL损伤的患者,无论是否接受过手术,都会呈现出膝关节炎和膝关节退行性病变 [8] [9]。为了使人们能够安全地参加和享受体育锻炼对健康的好处,预防ACL损伤可能是最有效的干预措施 [10]。
借助垂直跳跃下落实验可以预测膝盖受伤的风险 [11]。此方法是一种标准化的测试,测试者从箱子跳下,同时采集运动参数,并根据这些参数计算出膝关节的负荷。目前,该实验已经成为前交叉韧带损伤风险评估的筛查方法之一 [12]。
在过去的二十年间,有大量研究试图确定在生物力学和神经肌肉两个方面风险因素,作为发展预防训练的出发点 [3] [13] - [18]。这些研究共同表明,膝关节运动与其他下肢关节运动相关,膝关节负荷不仅受到自身运动的影响,还受到运动链中其他关节和节段的影响,尤其是髋关节和躯干的运动。后地面反作用力是影响ACL载荷的重要因素,髋关节和躯干运动对它影响显著 [19]。
考虑到非接触性ACL损伤的性别差异,有研究试图通过比较男性和女性运动员的运动模式来确定容易诱发非接触性ACL损伤的危险因素 [13],为了了解性别差异和ACL的损伤机制是如何关联的,有研究者记录了一系列特殊的运动过程中女性受试者个体间的差异,并进行了长时间的跟踪调查,通过比较受伤对象和未受伤对象在其中的动作差异,确定了可能导致ACL损伤的因素。
后来的一项研究建立了一种随机的生物力学模型 [19],这个模型模型模拟了停跳(stop-jump)任务中前交叉韧带负荷峰值时的运动特征,并准确预测了前交叉韧带负荷峰值时膝关节屈曲角的男女损伤率比和变化。Julia Kar等人提出了一种基于计算机仿真的方法 [20],这种方法通过逆向动力学的方法,由运动过程反推受力情况。他们建立了骨肌模型并模拟女性运动员的停跳运动,探讨了膝关节运动学与动力学,肌肉兴奋与激活,ACL张力与内力及对称性问题,用于探究ACL负荷与这些因素之间的关联机制。
这种基于计算机仿真的方法提供了一种前交叉韧带损伤相关变量的数值估计的新思路,即通过采集运动数据,输入计算机进行逆向动力学仿真分析,最终得到受力分析结果。然而,目前最常见的一种运动采集方法是基于多摄像头,受试者穿戴标志物的光学式运动采集。虽然目前对于运动数据的精度没有一个严格的定义,但是这种方法是目前公认的金标准。然而,这种运动采集的方法对操作者有着较高的技术要求,因为对于每一位不同的受试者,标志物的安装位置都是不同的,而标志物的摆放位置会直接影响所采集数据的准确性。同时,由于这种采集方法要求同时使用大量摄像头,就限制了其应用场景的范围和灵活性;此外,这种方法要求受试者在运动全程佩戴大量的标值物,这无疑会对其运动过程产生明显的干扰。
本研究提出一种基于双Kinect采集运动数据并通过AnyBody Modeling System™ (以下简称AnyBody)软件进行逆向动力学仿真的生物力学分析方法,运动捕捉过程无需标志物,具有成本低、耗时短的优势。通过对着陆时不同屈曲角度结果的分析,可以便捷有效地分析出跳跃着陆过程中膝关节的生物力学情况,对预防膝关节以及ACL的损伤具有很高的指导意义。
2. 实验材料和实验方法
本文所设计实验流程如图1所示。实验采用了基于双Kinect v2的无标记物运动捕捉系统,用于采集受试者运动过程。使用AnyBody软件进行仿真分析。Kinect v2是一款由Microsoft公司开发的一款同时集成了彩色和深度摄像头的体感设备,AnyBody是丹麦奥尔堡大学(Aalborg university)开发的一款用于运动过程生物力学分析的人体建模仿真软件。
![](//html.hanspub.org/file/9-2570250x10_hanspub.png)
Figure 1. Flow chart of knee load analysis during knee drop landing
图1. 膝关节跳跃着陆过程膝关节负荷分析流程图
测试者为5名健康成年男性,年龄(21 ± 1)岁,身高(174 ± 9) cm,体重(52.5 ± 11.5) kg,均无下肢骨骼肌肉方面疾病,且未曾接受膝关节外科手术。所有测试者均为自愿加入实验并签署了知情同意书。在实验开始前,测试者均接受了详细的动作指导,并进行充分热身。实验根据着陆后膝关节屈曲所达到的最大角度的不同,分为(1)较大最大屈曲角度组和(2)较小最大屈曲角度组两组。每组实验均采集实验数据20人次,最终分别选取不同组别下5组有效数据用于结果分析。
测试者站在高度为45 cm的平台上,平举双臂,并且在整个跳跃至着陆过程中始终保持双臂的平举姿势,身体正直,面向前方,自高台跳落,落点控制在高台前方地面标记位置5 cm范围内,保证着陆时身体重心维持稳定。小屈曲角度的实验中,双脚与地面接触后主动限制膝关节屈曲角度在80度至100度之间,如图2(a)所示。大屈曲角度实验中,双脚着地后不主动控制膝关节弯曲,放松状态下屈膝至最大屈曲角度约160度,如图2(b)所示。
![](//html.hanspub.org/file/9-2570250x11_hanspub.png)
Figure 2. Landing at different buckling angles
图2. 不同屈曲角度的着陆动作
建立AnyBody下肢模型,将采集到的运动数据导入AnyBody软件,并对肌肉骨骼模型进行个性化处理,肌肉部分采用AnyBody中的三元素肌肉模型 [21] [22] [23] (AnyMuscleModel3E),其原理图如图3所示。这是一个完整的希尔模型,考虑了肌肉的平行被动弹性、肌腱的连续弹性、纤维的夹角以及许多其他特性。CE代表肌肉纤维的活性,T为串联弹性元,代表肌腱的弹性,PE为并联的弹性元,代表肌肉纤维的被动刚性 [24]。
![](//html.hanspub.org/file/9-2570250x12_hanspub.png)
Figure 3. Schematic diagram of muscle model
图3. 肌肉模型原理示意图
为了在韧带松弛区产生非线性的弹性特性 [25],AnyBody模型中的韧带受力与型变的关系由以下公式定义:
式中,f表示韧带上的张力,ε代表韧带应变,ε1 = 0.03为参考应变,k是刚度系数。刚度值和参考应变值均采用文献中的数值 [26] [27],其中ACL的刚度系数为10,000。包裹骨骼表面,以防韧带进入骨骼内部 [28]。
驱动人体模型进行逆向动力学分析,得到运动过程中膝关节各个部分的生物力学参数。使用SPSS 22.0软件对着陆时不同最大屈曲角度引起的差异进行方差分析并对各参数进行相关性分析。所有的数据结果都以均值±方差的方式表示,前交叉韧带拉力的值为前内侧束与后外侧束之和。
3. 实验结果及分析
AnyBody软件提供了一种预测地面反作用力的算法,使得模型在没有实际地面反作用力数据的情况下,根据运动过程预测所产生的地面反作用力。有研究表明,使用预测结果作为输入,所得到的分析结果,其水平地面反作用力,关节屈曲力矩和关节反作用力的准确度足以代替传统受力板采集的数据 [29]。
3.1. 垂直地面反作用力
由图4可知,采用不同的屈曲角度着陆时,垂直地面反作用力峰值出现的时刻大致相同,而小屈曲角度组中,垂直地面反作用力在峰值出现前后的变化更加剧烈,且在着陆后前期的垂直地面反作用力值始终较组一偏高。根据方差分析结果,不同屈曲角度下垂直地面反作用力峰值的差异无统计学意义(p > 0.05)。通过相关性分析,得到其与股直肌收缩力呈相关性系数0.992的正相关性。
![](//html.hanspub.org/file/9-2570250x14_hanspub.png)
Figure 4. Vertical ground reaction force changes with time
图4. 垂直地面反作用力随时间的变化
3.2. 股直肌收缩力
由图5可知,在以小屈曲角度着陆时,股直肌在着陆初期的收缩力较大屈曲角度着陆时更大,且变化更剧烈,且在接触地面后,股直肌立即收缩并迅速达到较大的峰值,之后收缩力的值迅速减小。大屈曲角度组股直肌出现收缩后,收缩力值没有出现剧烈的增大,股直肌收缩力变化趋势较缓,其峰值较小屈曲角度组相对较小。不同屈曲角度下,股直肌收缩力峰值存在统计学差异(p < 0.05)。
3.3. 腱肌收缩力
由图6可知,在着陆初期,腱肌的收缩力没有明显变化,证明腱肌肌群没有或只发生微弱的收缩,这一阶段股四头肌的作用占据主导地位。当屈曲角度达到最大值,膝关节开始伸展,此时腱肌肌群开始收缩,当膝关节逐渐由屈曲变为伸展状态,腱肌的收缩力也逐渐达到最大值。不同屈曲角度下,腱肌肌群收缩力峰值存在统计学差异(p < 0.05)。经相关性分析,其与垂直地面反作用力呈相关性系数0.865的负相关性,与股直肌收缩力呈相关性系数0.886的负相关性。
![](//html.hanspub.org/file/9-2570250x15_hanspub.png)
Figure 5. The contractile force of the rectus femoris varies with time
图5. 股直肌收缩力随时间的变化
![](//html.hanspub.org/file/9-2570250x16_hanspub.png)
Figure 6. The contractile force of tendon muscles varies with time
图6. 腱肌肌群收缩力随时间的变化
3.4. 膝关节前后剪切力
由图7可知,在着陆刻,膝关节所受前后剪切力最大,且着陆后膝关节屈曲角度较小时,受到的膝关节前后剪切力远大于大屈曲角度较大的情况,小屈曲角度组的剪切力峰值几乎是打屈曲角度组的两倍。不同屈曲角度下,膝关节前后剪切力峰值存在显著的统计学差异(p < 0.001)。着陆后膝关节的前后剪切力迅速变小,由于大屈曲角度组的膝关节屈曲程度大,胫骨前移的趋势更剧烈,因此在着陆后的过程中,其膝关节前后剪切力逐渐高于小屈曲角度组。
3.5. 膝关节屈曲力矩
由图8可知,打屈曲角度组中由于达到最大膝关节屈曲所需的时间更长,膝关节屈曲力矩峰值出现的时间也更晚。与小屈曲角度组相比,在着陆后的整个过程中,大屈曲角度组的膝关节屈曲力矩均较大。不同屈曲角度下,膝关节屈曲力矩峰值存在统计学差异(p < 0.05)。
![](//html.hanspub.org/file/9-2570250x17_hanspub.png)
Figure 7. The change of shear force with time before and after knee joint
图7. 膝关节前后剪切力随时间的变化
![](//html.hanspub.org/file/9-2570250x18_hanspub.png)
Figure 8. Change of knee flexion moment with time
图8. 膝关节屈曲力矩随时间的变化
3.6. 前交叉韧带拉力
由图9可知,两组前交叉韧带拉力峰值出现的时间大致相同,而小屈曲角度组的峰值较大屈曲角度组更大,由表1可知,小屈曲角度韧带拉力峰值约为大屈曲角度峰值的2倍。不同屈曲角度下,前交叉韧带拉力峰值存在统计学差异(p < 0.05)。在整个着陆的过程中,小屈曲角度组ACL拉力均明显高于大屈曲角度组。经相关性分析,得到其与垂直地面反作用力呈相关性系数0.798的正相关性,与股直肌收缩力呈0.783的正相关性。
![](//html.hanspub.org/file/9-2570250x19_hanspub.png)
Figure 9. Change of ACL tension with time
图9. 前交叉韧带拉力随时间的变化
![](Images/Table_Tmp.jpg)
Table 1. The peak values of various properties of the knee under different degrees of flexion
表1. 不同屈曲程度下膝关节各项属性峰值
4. 讨论
从力学角度讲,ACL损伤是由于其受到了超过极限负荷能力的拉力。根据Berns,Hul的研究小组和Markolf等人的体外研究发现 [30] [31],胫骨近端施加的前剪切力是前交叉韧带加载的主要因素,而在胫骨近端施加前剪切力时,膝关节外翻、内翻和内力矩对前交叉韧带加载也有显著影响。一系列结果表明,在非接触ACL损伤中,股四头肌力可能是导致ACL负荷的主要因素 [3],尤其是当膝关节屈曲角度较小时,股四头肌更有可能造成更大的剪切力。本实验结果发现,股直肌是膝关节伸肌肌群股四头肌中的一条主要肌肉,在着陆后,人体重心会继续下降,此时膝关节做屈曲运动,股四头肌收缩牵拉髌韧带,通过髌骨的杠杆作用,使膝关节屈曲减缓。在股四头肌收缩时,会使髌韧带紧张产生指向前上方的拉力,由于髌韧带连接在胫骨前侧结节处,髌韧带的拉力会在髌骨上端产生一个前向的剪切力,使胫骨产生前移的运动趋势。以较小的屈曲角度着陆时,股四头肌收缩力更大,且肌肉激活时间更早,胫骨前移的趋势更明显,而由于胫骨前移,前交叉韧带产生了更大的应变,从而加剧了前交叉韧带的负荷。
由相关性分析结果可知,腱肌收缩力与垂直地面反作用力以及股直肌收缩力均呈负相关。在维持膝关节稳定的过程中,股四头肌和腱肌分别会对胫骨近端产生向前和向后的拉力。作为股四头肌的拮抗肌,腱肌收缩时,可以减少股四头肌在胫骨近端造成的前向位移,从而减小前交叉韧带的应变。在采用较大膝关节屈曲角度着陆时,腱肌收缩力的峰值较大且作用时间更长,有效地限制了膝关节在小的屈曲角度下的胫骨前移,减小了前交叉韧带的负荷。
通过实验结果,可以看出垂直地面反作用力与股直肌收缩力呈现正相关的关系。垂直地面反作用力使身体获得向上的加速度,使身体重心减速下落,直至下落至最低点。双脚受到的垂直地面反作用力经过下肢传导至膝关节造成膝关节关节面之间的压力。采用更大屈曲角度着陆时,由于着陆后下落至最低点的时间更长,所需的减速加速度较小,使得两脚受到地面的冲击力更小,从而减小了膝关节的压力。并且,由于股直肌收缩力与垂直地面反作用力呈现正相关的线性关系,更小的垂直地面反作用力也有利于减小股直肌收缩力,从而减小膝关节前向剪切力,从而减轻前交叉韧带负荷。因此,采用更大的屈曲角度着陆可以减轻膝关节压力,减小膝关节前向剪切力和前交叉韧带负荷。
5. 结论
通过本次实验的结果和讨论,可知在跳跃着陆过程中,当着陆后膝关节屈曲角度较大时,垂直地面反作用力较小,因此传导至膝关节的冲击力较小,有利于减小膝关节负荷,保护膝关节。而当膝关节的最大屈曲角度较小时,身体重心下落的距离更短,减速时所需要的加速度更大,因此,股四头肌会产生更大的收缩力,使膝关节停止进一步屈曲。而更大的股四头肌收缩力也就在胫骨近端的髌韧带附着点处产生了更大的拉力,它在膝关节造成的前向剪切力是导致前交叉韧带负荷的主要因素。腱肌是位于大腿后侧的一组肌肉,它们与股直肌(股四头肌)之间存在拮抗作用,在膝关节屈曲角度较大时,腱肌的收缩力相较组二有明显的增大,且出现作用时间早,这有利于减小膝关节前向剪切力。综上所述,跳跃着陆过程中,在着陆后应该使膝关节屈曲至较大角度以减小地面冲击力和股直肌的收缩力,从而起到减小膝关节负荷和前交叉韧带应力的作用。
NOTES
*通讯作者。