摘 要：

随机挑选200名实验对象，记录他们的基本信息、关节角度和滑行速度，对其中181个有效数据进行分析.根据1 名志愿者的 MR，CT 数据建立了膝关节140°软骨与骨的三维模型.根据运动学分析结果和膝关节解剖学结构定义载荷及约束，在 Abaqus 6. 14中得出转弯膝软骨应力集中出现的位置.研究结果表明体重过大、速度过快时应适当增大转弯半径以减少向心力对膝关节软骨施加的应力；犁式转弯的过弯速度应控制在3 m /s 之内.

关 键 词 ：膝关节软骨;有限元分析;犁式滑雪;转弯侧膝关节;生物力学

双板滑雪初学者多采用犁式滑行.由于身体过度紧张、肌肉僵直，经常出现不能有效减速、转弯等情况,从而导致碰撞或摔倒.因此，对于初学者而言，膝关节存在一定的损伤风险.犁式滑雪普遍速度较慢,小腿呈“内八”姿态，在转弯时大多数体重转移至转弯腿膝关节内侧.在错误动作、速度 、体重、关节角度等多种外界条件的影响下，内侧半月板存在极大的损伤风险.在犁式转弯阶段，对转弯侧膝关节软骨进行有限元分析，揭示犁式转弯过程中膝关节软骨损伤机理.

1 实验对象和方法

1 . 1 实验对象

本研究实验对象来自大众滑雪人群.在沈阳体育学院白清寨滑雪场随机选出男女各100名初学者，并从中挑选1 名无膝关节损伤史的志愿者.

1.2 实验方法

通过影像分析得出犁式滑雪转弯时转弯腿膝关节角度和转弯速度.对1 名志愿者进行膝关节C T 和 M R 检 测 ，将 得 到 的 DICOM数据导人Mimics Research 17. 0 x 6 4 软件中进行图像分割与重建.分别将CT和 M R数据进行阈值分割，从某些像素点出发，以组织边缘为界限进行逐层修复，最终实现目标提取.对滑雪模型进行3 维度的受力分析，模拟实际运动过程并结合解剖学相关内容对各模型施加约束，推导出有限元分析中所需的边界条件，对膝关节模型中的半月板进行有限元分析.

2 结果与讨论

2 . 1 采集结果

对实验对象信息、膝关节角度、转弯速度和问卷等数据进行分析，共整理出有效数据181个.其中男子87个 ，女 子 94个 ，统计数据见表1.成功获 取 1 位无膝关节损伤史的27岁男性志愿者的膝关节C T 和 M R数据，得 到 C T 平扫图片377张， MR矢状面224张，均为DICOM格式.在犁式转弯过程中，人体重心向转弯膝侧偏移，约占体重的6 0 % 〜8 0 % ，取中间值70% .根据采集的关节角度数据可知:在犁式滑行整个过程中膝关节始终在112°〜169°之间往复运动，平均在 140°左右持续时间最长、频率最高.

2. 2 建模结果

以选取的骨组织区域为例，将股骨用单独颜色标记.经过手工填补缺失部分，应 用 Calculate3D 将股骨三维重建[4]，如图1 所示.分别对胫骨、腓骨、髌骨和股骨等关键部位进行三维重建.在实际情况下各骨组织之间通过韧带相互连接.由于本研究只针对膝关节软骨展开，因此骨模型的建立只是为了确定软骨模型的准确性.基于以上模型的建立并利用M R数据对膝关节软骨组织进行提取，最终建立膝关节140°时软骨组织模型，如图2 所示.

2. 3 力学模型

将犁式滑雪模型简化为棍图后进行XKZ三个方向的受力分析，如图 3 所示.其中V 是犁式转弯时的最大速度;G'和 G"分别为转弯侧足底和转弯膝所受的重力;F'为转弯膝所受的向心力;，N'与N" ，R分别为足底和转弯膝所受的支持力为极限时的转弯半径;α和θ分别为转弯侧膝关节角度和场地坡度;f' 为速度方向上的滑动摩擦力.随后分别在XY'坐标轴（由于雪板与雪面之间存在夹角γ，故胫骨平台面与雪面也存在一定夹角γ'，这里λ≈ γ' ，γ 的大小与转弯半径负相关，如图4所示）和 XZ坐标轴下进行分析，如图5,6所示.

假设在犁式转弯过程中，动能没有损耗，为求出有限元分析中的F 和 F'的大小，需要将最大速度转换为角速度，即ω = V/R 双膝关节承受膝以上体重约为总体重的85. 6 % ，因此转弯膝此时受到的重力约为G'=70% x G x 85.6% =402N，转弯膝所受向心力约为F'= m.ω²R= 195 N，其中m' = G x 85. 6 % 为膝关节以上环节的体重，R=0. 6 m.此时，若只考虑转弯膝受到向心力作用，则膝关节所受F 和 Z 向 分 力 : F'γ=m'ω²R·cosγ'，F'z= m'ω²R*sinγ'.当人体保持直立，膝关节内外侧所受的重力分布均匀.犁式转弯时，膝关节屈曲，胫骨内旋，此时膝关节内侧压力增大，约为总载荷的70% .在图5 中，f"≈f'，为速度方向上滑动摩擦力的近似值;F"为膝关节所受向心力的反作用力，其大小为内侧和外侧所受向心力的合力.本研究中雪的摩擦系数μ =0.1，因此f"≈f' = N'μ.

2 . 4 有 限 元 分 析 结 果

2. 4. 1 有限元模型参数设定与有效性验证

使用Abaqus 6. 14对膝关节软骨进行生物力学特性分析.由于相关文献中鲜有提及人种之间半月板和关节软骨的弹性模量、泊松比和软骨密度差异问题，故设置半月板和关节软骨的弹性模量分别为59,20 MPa，泊松比分别为0. 49,0. 46,密度分别为 1 x 10'1. 98 x 1 0 kg/rnm3[5-9].模型中共计50 105个节点，26 839个单元（C3D10).选取内侧半月板，对其进行破坏性试验.结果内侧半月板最大压力约为789 N，最大抗压强度约为35. 8 MPa.此测试结果与尸体生物力学实验的测试结果较接近，说明本研究所建模型与实际膝关节半月板的材料力学特性具有相似性.

2 . 4 . 2 膝关节软骨有限元分析结果

在本研究中，股骨与股骨软骨、胫骨与胫骨软骨之间定义为面与面接触且无滑动约束，内外侧半月板上下表面与胫骨软骨和股骨软骨之间、股骨软骨与胫骨软骨之间也为面与面接触，其中只有内外侧半月板前脚与后脚处由于有韧带连接故存在绑定约束.利用所建立的有限元模型对犁式滑雪转弯膝下滑转弯过程进行了分析，如图7 所示.可知：在 膝 关 节 角 度 为 140°，下滑速度为3 rn/s 时，在自身重力的作用下，膝关节存在应力集中现象;股骨软骨不存在应力集中现象;内侧半月板出现明显的应力集中现象，位置出现在中部和后脚处;外侧半月板出现较为明显的应力集中现象，位置出现在中部和中后部；内侧胫骨软骨出现较为明显的应力集中现象，位置出现在中后部及内侧边缘中部;外侧胫骨软骨出现明显的应力集中现象，位置出现在外侧边缘中部.

3 结 论

1 ) 在下滑过程中应力集中点出现在犁式转弯时转弯膝内侧半月板中部和后脚处、外侧半月板中部和中后部、内侧胫骨软骨中后部和内侧边缘中部及外侧胫骨软骨外侧边缘中部.

2 ) 从内侧半月板受力并结合解剖结构来看，内侧半月板的应力会使其中部挤压到内侧副韧带 ，使得内侧副韧带更加收紧，此种情况是导致犁式滑雪过程中膝关节内侧酸胀的成因之一.

3 ) 在规范的犁式滑雪技术动作前提下，合理地调整运动强度、滑行速度和转弯半径，可有效降低犁式滑雪时对膝关节软骨的损伤.