2.1   有限元法在医学中的应用   有限元法是数值计算中的一种离散方法，包括矩阵方法在结构力学和弹性力学领域的发展和应用，通过建立力学模型，将其结构荷载、几何形状、材料性能、边界条件等均用模拟结构来求解。目前有限元分析在骨科生物力学研究中大致有以下几个方面(图2)，具体应用时可进行通用模型构建、运动模拟等指导患者避免高危动作延缓病情，为植入物的设计、选择及术前规划提供参考，预测治疗效果、病情转归及预后，以及护具研发及优化，预测手法治疗的变化等等。"

2.2   有限元模型的建立   目前研究者利用电子计算机断层扫描(CT)、磁共振成像(MRI)所得的图像进行有限元模型的建立。首先将扫描图像数据导入Mimics等软件，选择适宜的灰度值，根据骨骼与软组织成像的差异性，提取出骨骼轮廓，并对图像进行阈值分割、蒙版编辑、区域增长、填充、光滑等等操作后，进行三维重建从而得到最初的髋关节骨骼三维模型。将所建立的模型导入3-matic、 Geomagic studio等软件自动生成面网格，通过优化、过滤网格、填补空洞、拟合曲线等操作，从而创建体网格。在Mimics中赋予材料属性后使用有限元分析软件如 ANSYS、ABAQUS等进行骨骼、关节、软组织、内固定植入物和假体等的力学分析。见图3。"

有限元法可基于影像学图像和解剖特征建立高精度人体仿真模型，利用有限的、相互关联的单元对人体复杂的关节进行简化，通过计算得到直观准确的分析结果，能降低分析成本；利用有限元法可以获得在体内无法测量的信息，由于其在虚拟环境中操作，不受时间和伦理的限制，因此在医学领域，尤其是骨科、生物力学等方面有重要应用价值。但有限元建模过程中将人体复杂结构简化再进行分析是固有的局限性；操作过程中影像学资料在不同的软件之间导入、切换，可能导致重要信息丢失，应将软件完善、统一，以便在相似研究之间横向比较；其中材料参数赋值、载荷施加等步骤直接影响结果，如何确保模拟实验前输入参数的准确是一大挑战。 2.3   正常与DDH患者的髋关节力学特点   髋关节是由股骨头和髋臼共同构成球窝关节，连接人体躯干和下肢，是最大的承重关节。与髋关节生物力学有关的力学因素主要是外展肌力和关节反应力，其中关节反应力是外展肌力与身体重力的合力，并作用于髋关节[18-19]。目前相关的有限元研究也证实，髋关节正常负重时，应力均匀分布于股骨头上，而DDH患者髋关节头臼不匹配，股骨头半脱位或全脱位，髋关节力学关系异常，导致外展肌无力、关节软骨磨损加重、髋部疼痛等，进而继发骨关节炎[20-21]，CHU等[22]的研究证实DDH患者软骨下骨小梁骨微结构恶化，可能与软骨损伤有关。髋臼盂唇承担了部分髋关节的接触压力，当盂唇撕裂或切除接触压力增加，而修复或重建唇唇后，接触压力降低[23-24]。DDH的发生是由髋关节的先天性发育缺陷和后天的生物力学异常共同作用的[25]，这是一个动态变化的过程，髋关节结构在生长发育过程中持续不断进行的不良改变[26-28]；在儿童DDH有限元模型中，髋臼与股骨头接触面积减小，应力增大并集中，且股骨颈局部也存在应力集中的情况[29]；生长中的骨骼对负荷条件敏感，负荷的变化能影响骨骼的形状和功能，MRI显示儿童DDH股骨头生长发育障碍，通常呈非球形且较健侧小，在复位成功后，可以发生显著的重塑改变，同时治疗的时机与时间长短也非常重要[30-32]。 2.4   DDH各种治疗方法的有限元法研究   0-6个月的DDH患儿使用Pavlik吊带治疗有> 90%的成功率[33-34]。由于婴儿髋关节未完全骨化，MRI相对于X射线片及CT可提供软骨性髋臼指数。VAFAEIAN等[35-36]的研究以MRI图像为基础，为正常婴儿的髋关节建立有限元模型并模拟其在佩戴Pavlik吊带时软骨性髋关节的软骨接触压力，显示髋臼的前部、后部和髋臼上内侧存在马蹄形的应力区域，而上外侧处于无负荷状态，在这种接触压力的分布下髋臼倾斜度增加、髋臼加深，髋臼顶仅承受轻微负荷(正常髋关节)或卸载(发育不良髋关节)；研究模拟髋关节在屈曲90°、外展40°、60°和80°时的接触压力，发现随着髋关节外展角度增加(特别是> 75°)则接触压力增大。MáRQUEZ-FLóREZ等[37]的研究也表明，与站立位相比，婴儿处于仰卧髋关节屈曲外展时髋臼面接触力的分布更加均匀；而在极端外展(80°)时髋臼后壁和邻近的外侧股骨头的接触压力有增大和集中的趋势，外侧股骨头上的接触压力可能会减少股骨头血供，从而导致股骨头坏死[38]。HUAYAMAVE等[39]研究Pavlik吊带治疗严重脱位时前倾角在30°-70°之间的生物力学特性，模型分析结果表明减少股骨前倾角可能会增加治疗高度髋关节发育不良的成功率。这些研究从生物力学角度证明了Pavlik吊带治疗成功和DDH并发症发生的潜在生物力学因素。 闭合复位和髋人字石膏固定是6-18月龄DDH患儿的推荐治疗方法。蔡振存等[40]对1例DDH患者“人”字型石膏固定后的髋关节进行有限元分析，结果表明患侧髋关节头与臼间的平均接触压力、最大接触压力、股骨头所受Von Mises应力峰值均高于健侧，患侧头臼接触面积小于健侧。LIU等[41]的研究采用有限元法比较闭合复位治疗前后髋臼和股骨头上Von Mises应力的大小和分布特征，结果显示治疗后髋臼顶部外侧缘和股骨头外侧柱的高应力点消失，应力呈放射状均匀分布并向外周逐渐减弱，髋关节的应力大小和分布与健康侧相似，从而以力学的角度佐证了闭合复位治疗的有效性。ZHANG等[42]对1例右侧DDH的16月龄女性患儿的髋关节进行了三维有限元分析，发现随着髋关节外展角度的增加，髋臼后部的应力明显增加，股骨头外侧的应力也随之增加，闭合复位时髋臼上区软骨接触压力变化不明显，股骨头受压无明显差异，从力学角度说明过度外展时股骨头外侧应力的增加可能与缺血坏死有关，而屈曲角度的增加与之无关。 对于年龄在12-18个月或更大、以及闭合复位失败的患儿，建议采用骨盆截骨术的方式以实现髋关节同心圆复位[43]。NAKAMURA等[44]比较了26例接受Salter骨盆截骨术的女性单侧DDH患者术前术后的CT、MRI三维形态，发现患侧髋臼软骨经手术治疗后改善至正常形态。以往认为Salter截骨以耻骨联合软骨为铰链，Pemberton截骨以Y型软骨为铰链，Dega截骨则以坐骨大切迹后内侧皮质层为铰链；而SUVOROV等[45]使用2岁和6岁患者的健康骨盆CT数据建立模型并模拟这3种骨盆截骨术进行相关生物力学分析，发现在骨盆这个整体单元中，Salter、Dega和Pemberton截骨术没有单一的铰链；Salter截骨术中患侧Y型软骨后支的负载在2岁模型上较大，对侧半骨盆的负载在2岁和6岁的模型中均较大；Pemberton截骨术中，前、后Y型软骨为主要受力部分，对侧耻骨的负载在2岁模型上出现；Dega截骨术以坐骨大切迹区髂后皮质层为主，对侧半骨盆未见负荷。研究同时发现骨盆韧带在骨盆截骨术中发挥着重要的限制作用；年龄的不同对生物力学有影响，2岁患者的骨盆软骨铰链较多，其模型比6岁患者的更“软”，这表明髋臼畸形更易于矫正。TOYOHARA等[46]对髋关节发育不良手术前后骶髂关节应力分布进行分析，结果表明髋臼形态对骶髂关节有影响，髋臼后覆盖的改善可能有助于恢复骨盆的负荷传递，从而改善髋臼发育不良患者骶髂关节的应力环境。 对于早期未得到正规治疗或者漏诊的DDH患者，其髋臼软骨生长中心闭合，X射线片表明髋臼较浅导致股骨头未覆盖，髋关节一致、X射线片上没有关节炎变性的迹象，采用髋臼周围截骨术治疗以改善髋关节疼痛等症状，早期的有限元研究从关节接触面、关节接触压力、股骨的应力等的角度证明髋臼周围截骨术的有效性[47]。PARK等[48]的研究表明髋臼周围截骨术中冠状面旋转是获得最佳临床结果的主要手术参数之一；而近期的一项有限元研究发现仅冠状面矫正可能无法使髋关节生物力学正常化[49]；KITAMURA等[50]通过虚拟髋臼周围截骨术操作，将髋臼骨块在冠状面上外旋，使得外侧中心边缘角达到30°，在矢状面上对髋臼骨块分别前旋转0°，5°，10°和15°，对这4个角度的髋臼进行有限元分析，结果表明在不同的前旋角度下，外旋髋臼可使关节接触压力正常化，且进一步的前旋髋臼可以进一步降低关节接触压力。KITAMURA等[51]的另一项研究以51例髋关节发育不良患者的三维髋关节模型分析髋臼周围截骨术前和术后单腿站立期间髋臼软骨的关节接触压力，研究结论是髋臼周围截骨术可以使64.7%的患者关节接触压力恢复正常，但对于具有非球形股骨头的患者，关节接触压力的恢复较困难。随着对该病的认识加深以及相关研究的进一步进展，AITKEN等[52]建立了20例髋臼周围截骨患者术前、术后的髋关节模型，在机械最佳、临床最佳和手术达到的定位之间，比较了X射线影像覆盖范围、接触面积、峰值/平均接触应力等的指标，得出X射线片表现不能说明髋臼周围截骨术后的关节接触力学达到了最佳的结论，对此则需要充分的术前规划、计算机模拟最佳截骨位置、3D打印截骨模板等措施以达到临床最佳[11]。 上述DDH髋关节以及各阶段治疗相关的有限元分析(表1)，验证了髋关节发育不良时的应力关系，从力学角度阐释了DDH髋臼及股骨产生异常的原因；对于保守治疗的病例，用有限元的方法验证其治疗的有效性，以及髋关节所处位置的角度，证明了“安全角”的重要意义，“安全角”过小则髋关节复位不稳定，过大则导致髋臼后壁及邻近外侧股骨头的软骨接触压力增大，最终影响股骨头血供。利用有限元方法从力学的角度验证手术治疗有效的原因，同时分析经过手术后骨盆力学关系的变化，以推测预后以及指导康复，并基于模拟实验结果指出存在的不足和需要改进的方向；除此之外，对于需要手术治疗的病例，可以借助有限元法进行术前规划，确定截骨的最佳位置以指导手术操作更加精确，同时可以适当降低手术创伤以及减少术中X射线照射的次数，这对患者是非常有利的。然而这些有限元分析中，较少提及髋臼横韧带、关节囊、肌肉等结构在髋关节中的力学作用，这与有限元建模、参数设置等的难易程度有关，如果在有限元模型中添加韧带、软骨、肌肉等元素，那么必须获得这些组织的3D信息，以及各自准确的材料参数，这是目前需要克服的一大困难。"