计及半月板的人体膝关节动力学仿真及三维有限元研究-dynamic simulation and three-dimensional finite element study of human knee joint considering meniscal.docx

第1章绪论1.1课题选题背景及意义人体关节结构的力学性能一直都是生物力学重要的研究课题,膝关节生物力学的研究，目前主要应用在以下几个方面：可以模拟膝关节不同运动、对线不齐以及疲劳损伤等多种情况下的关节面间的接触应力和应变；对关节假体进行疲劳损伤实验模拟,预测假体的使用寿命；模拟置换后骨的重建过程。目前主要用在研究与应力有关的骨结构、人工关节的设计和固定假体微动松动以及骨改建过程等[1]。近年来随着计算机技术的提高和计算速度的不断提高，有限元分析和应用在工程（包括生物工程）设计和分析中得到了越来越广泛的重视，可模拟活体组织的生物力学实验并多重加载[2]；通过模拟复杂的几何形状结构，提供模型任意部位的应力和位移情况[3]；能得到那些无法通过实际测量和客观实体方法难以得到的研究结果；并且研究可重复性好，模拟的精度高，已经成为复杂生物力学问题研究的重要手段之一。目前国内外的一些研究者，或者忽略了半月板，或忽略了韧带、肌肉的影响，或假定股骨远端其固定，且大多模型都是静力和准静力模型，对膝关节的完整结构进行了假设，使得分析出来的计算结果与实际的差距较大[4-8]。因此，建立完整的人体下肢膝关节的动态生物力学模型，从而确定膝关节内的动态应力应变规律就成为本文研究的主要内容。本文通过建立一个完整步态周期的人体膝关节动力学模型，对人体下肢步态行走进行了数值分析，模拟得到人体关节结构的动态力学参数，用于膝关节的三维有限元模型，通过计算分析，从而确定了膝关节内部的动态应力、应变规律，为研究复杂形状、复杂载荷和复杂材料的人体关节提供一种研究思路。1.2人体膝关节生物力学计算模型的国内外研究概况随着图像学和计算机科学的发展应用,利用计算机重建的虚拟人体膝关节三维模型和有限元模型被越来越广泛地应用于计算生物力学的研究中[9-10],依靠医学影像数据建立符合人体解剖结构的膝关节模型,是目前精度较高的建模方法之一。国外自19世纪50年代起就已经开始医学图像三维重建方面的研究，对CT、MRI逆向以及人体骨骼个性化设计制造的研究起步很早，发展迅速，也取得了丰硕的成果[11]。国外自1965年起，就开始了对膝关节的生物力学模型进行了研究。经过了几个发展阶段：1.通过建立数学模型对其进行分析。利用基础的数学关系式，通过计算得到。Paul[12]早在1965年就做出了关于膝关节受力的研究，他得到的关节受力曲线的结果，目前仍是许多学者普遍采用的一种标准载荷曲线。Taylor[13]等利用下肢的骨骼--肌肉模型,通过变换载荷的大小,估算出了在走路和爬楼梯时胫骨--股骨关节的接触力。2.直接利用有限元平台，在其上利用关联点建立模型，提高了分析的可靠度，但模型的精确度和与实际的人体结构有些误差。Costigan[14]等研究了爬楼梯时髋和膝关节的运动。Bendjaballah[15]等模拟出了四条韧带保留完整的情况下膝关节内翻、外翻时的接触应力大小。关于胫骨平台的强度测试，Goldstein和Volz[16]的研究表明，随着距离胫骨软骨层面的距离大小的变化，关节面的强度也在发生变化。Johnson[17]等描绘了近端胫骨、松质骨的压缩强度示意图，发现了接近近端胫骨前缘的松质骨的强度，由近向远是迅速下降的。3.通过坐标扫描仪得到的模型，进一步提高了模型的精确性。但在扫描进入有限元软件中时，会造成边缘数据的丢失，影响分析的结果。Watanabe[18]等对已做过活动轴关节置换手术的膝关节深度弯曲时的运动进行了分析研究。Martelli[19]等建立了全膝关节置换术后的双关节面的髌骨关节模型,但该模型只考虑了髌骨与股骨滑槽的接触状态。事实上,当屈膝超过某种程度时,髌骨就会离开股骨滑槽而与股骨内外髁接触。Matsuda[20]等使用MRI测量的30个正常膝关节和30个内翻畸形膝关节的胫骨平台（PSA），发现正常膝关节的内侧平台（PSA）为5°-15.5°，平均为10.7°，正常膝关节的外侧平台（PSA）为0°-14.5°，平均为7.2°。4.利用逆向工程建立的三维有限元模型，把建立的模型以文件或代码的形式书写，以节点坐标和节点的关联性即单元信息在有限元软件中进行识别，是现在大多采用的研究方法，其精确性高、可靠度好。如：Tomaso[21]等通过建立膝关节三维有限元假体模型，模拟出了在载荷、弯曲角度下的接触应力、接触面积以及疲劳磨损情况,并与使用富士压敏片测出的结果进行了比拟对照。Winson[22]等利用建立的有限元模型，求出了膝关节各部位的接触应力的分布数据和剪应力的峰值。通过建立的有限元模型，Jia[23]等分析了正常步态下膝关节弯曲各角度时接触面积的大小和位置。Yang[24]等利用ABAQUS的有限元模型分析了在接触力为2800N时，膝关节分别旋转0°，4°，8°时的接触应力