收藏本站
收藏 | 手机打开
二维码
手机客户端打开本文

数字化人体足踝部三维有限元模型的建立及分析

刘清华  
【摘要】: 研究背景: 随着现代化的进程,交通业和建筑业空前发展,以及人们健康观念的提高,足踝部疾病发病率逐年上升,病种也日趋复杂,使骨科临床面临新的挑战,迫切需求对其进行深入研究。足是人体承重的根本,直接与地面接触,除了行走功能外,还能够吸收运动振荡、缓冲地面反作用力,保护躯体重要器官和组织免受伤害,并使人体的重心稳定。踝关节是人体与地面接触的枢纽,行走、跳跃、跑步和登高都需要踝关节的参与,即便骑自行车或驾驶汽车亦离不开踝关节的协调动作,可以说日常生活中的每一个动作都有踝关节的参与,因此踝关节也是最容易受到损伤的关节之一。踝关节力学机制复杂,各种损伤后都可能打破其周围结构的力学平衡而导致不稳定,诱发创伤性关节炎。而且,足踝部结构异常除导致疼痛、畸形及活动障碍外,还可进一步影响下肢、骨盆、脊柱等力学功能的正常发挥,故对其研究具有重要临床意义。 足踝部生物力学的研究是人体生物力学研究中的重要组成部分。许多研究者指出生物力学因素在探求足踝部疾病的病因机制、治疗和预防等方面起到重要的作用。以往对足踝部研究多采用传统实验生物力学测试的手段,但是用这种方法的缺点在于实验手段复杂,进行多种载荷工况下的实验往往耗资大、周期长、效率低,而且由于无法在人体上直接实验,因此很难准确地反映出各种不同受力情况下的应力分布规律。为了解决这一问题,近年来,随着数字化技术的提高,越来越多的研究者采用数值计算的方法进行分析,即以传统的力学分析理论为基础,应用数值分析手段例如有限元方法(finite element method, FEM)等,进行线性和非线性的应力和变形分析。利用FEM可以仿真模拟各种足踝部疾患,从而使得对足踝部复杂的骨胳几何结构、边界条件和材料的不均匀性等问题的生物力学研究有了可能的解决途径。 FEM具有强大的建模功能,在动静状态下能够对具有复杂的几何形状、材料参数和不同受力条件下的物体进行模拟仿真研究,其已经越来越多的被应用到人体生物力学中。因足踝部骨、关节、韧带、肌腱等结构复杂,在进行有限元建模时工作量较大,且缺乏完善的各种组织详细参数,有限元模拟相对于其他大关节的研究难度更大。现有建立的足踝部有限元模型,大部分对踝关节的解剖结构建立不够完善,主要用于与足底压力布相关的研究。然而踝关节损伤对人体功能的影响更加受到临床的重视,亟待对其更深入地研究。 本研究是在总结以前研究的基础上,对足踝部有限元模型的有效构建方法进行了长时间的探索,然后构建出了解剖结构比较精细的正常足踝部有限元模型,分别在模型中模拟各种踝关节损伤及内固定术式,探讨其对足踝部生物力学的影响。 目的: 1.依据1名正常男性志愿者踝关节中立位下螺旋CT扫描的右侧踝关节面上20 cm以下的足踝部影像学资料,探讨利用Mimics、SolidWorks、ANSYS等软件构建出解剖结构精细的正常踝关节三维有限元模型的方法,并对其有效性进行验证,使其能反映正常足踝部的力学特性。 2.根据所建立的正常足踝部的三维有限元模型,对其施加不同载荷,模拟人体静态中立位单足站立负重状态、踝关节受到内旋力、踝关节受到外旋力、踝关节受到内翻力及踝关节受到外翻力,探讨正常踝关节在不同活动方式中,主要关节周围各组织的应力、位移分布,以及关节活动度(range of motion, ROM)的变化等情况,同时对其结果的有效性进行验证。 3.根据所建立的正常足踝部的三维有限元模型,模拟建立下胫腓联合损伤的模型,并在下胫腓联合损伤模型基础上分别建立踝关节平面上方的2.5 cm、5cm螺钉固定模型,探讨各模型在模拟人体中立位单足站立、踝关节内旋、及踝关节外旋加载中主要关节周围组织的应力、位移分布,以及关节活动度的变化等情况,同时对其结果的有效性进行验证;并通过与正常模型在这3种加载中的计算结果对照,探讨下胫腓联合损伤及踝关节平面上方2.5 cm、5 cm螺钉固定后对踝关节应力分布及稳定性的影响。 4.根据所建立的正常足踝部的三维有限元模型,分别模拟建立后踝骨折累及1/5、1/4、1/3及1/2的踝关节面的损伤模型,探讨各模型在模拟人体中立位单足站立、踝关节内旋、及踝关节外旋加载中主要关节周围组织的应力、位移分布,以及关节活动度的变化等情况,同时对其结果的有效性进行验证;并通过与正常模型在这3种加载中的计算结果对照,探讨后踝骨折分别累及1/5、1/4、1/3及1/2的踝关节面后对踝关节应力分布及稳定性的影响。 方法: 1.数字化足踝部三维有限元模型的构建与验证:采用南方医科大学南方医院影像中心的LightSpeed 16排螺旋CT,扫描参数:电压120~140 KV,电流强度240~300 mA,螺距1.375~1.75,层厚7.5 mm,矩阵512×512,重建层厚0.625mm。对志愿者右足自踝关节上20 cm胫腓骨远端向下扫描至足底,扫描时右足保持中立位,以DICOM格式输入个人计算机的Mimics 10.01软件,经自动或手动阈值分割后三维重建出完整足踝部的28块骨骼及外围软组织的三维结构,再以点云输出并导入SolidWorks 2009,利用网格处理向导及曲面生成向导生成几何模型并重建装配体,然后分别导入两种有限元分析软件:(1)方法一为导入有限元分析软件ANSYS 12.0的Workbench模块建立完整的足踝部有限元模型。在踝关节的接触面两侧根据关节间隙分别建立了0.5 mm关节软骨,其他关节采用仅受压的三维杆单元模拟软骨,并建立了125条弹簧模拟韧带和小腿骨间膜,5条梁单元模拟跖筋膜,材料属性分别参照文献确定。然后参照Anderson等方法,对模型仅取胫骨及距骨所构成的简单踝关节模型加以测试。通过模拟人体单足站立状态下的力学传递,对模型胫骨下端的上截面施加600 N垂直载荷,对距骨予以约束。测量踝关节胫骨下关节面的接触压力、接触面积,并将结果与前者进行对照。(2)方法二为导入Simulation建立简化的踝关节有限元模型。根据研究需要,分别建立了包含胫骨、腓骨、距骨、跟骨及舟骨的5骨装配体,以及还包括骰骨及3块楔骨的9骨装配体的有限元模型。在模型中用“仅延伸”的弹簧模拟韧带连接,其中在踝关节5骨装配体中共建立了31条弹簧,而在9骨装配体中建立了42条弹簧模拟踝关节周围韧带及小腿骨间膜等连接装置。定义各组织的材料属性,并通过自动或手动方法在各关节间生成接触对,设置相应的边界条件分别在踝关节5骨装配体有限元模型上模拟人体单足站立及踝关节内、外旋运动加载,而在踝关节9骨装配体有限元模型上模拟踝关节内、外翻运动加载,然后设置好网格划分密度生成网格并设置算例属性进行运算。 2.下胫腓联合损伤及固定术式的有限元分析:在前述方法二所建立的5骨装配体模型上分别将代表下胫腓联合前韧带、骨间韧带、后韧带及小腿远端8 cm的小腿骨间膜弹簧压缩(suppress),模拟建立下胫腓联合损伤的模型。然后在下胫腓联合损伤模型基础上分别在踝关节平面上方的2.5 cm、5 cm处平行于踝关节面并穿过腓骨、胫骨骨干的中央生成两个¢3.5mm孔,分别采用材料为316不锈钢的销钉连接两孔的圆柱面模拟建立两种固定方式的螺钉固定模型。然后设置好网格划分密度生成网格,对各种模型分别模拟人体中立位单足站立及踝关节内旋、外旋运动三种踝关节受力方式进行加载运算,并将计算结果与正常模型在这3种加载中的计算结果相对照。 3.后踝骨折的有限元分析:在前述方法二所建立的5骨装配体模型上,在模型的胫骨下端踝关节面上按照文献报道的标准,分别切除相应大小的骨块,模拟后踝骨折累及1/5、1/4、1/3及1/2的踝关节面,同时删除累及损伤的韧带,建立4种后踝骨折损伤模型。然后设置好网格划分密度生成网格,对各种模型分别模拟人体中立位单足站立及踝关节内旋、外旋运动三种踝关节受力方式进行加载运算,并将计算结果与正常模型在这3种加载中的计算结果相对照。 结果: 1.在按方法一所建的完整的足踝部有限元模型中,将胫、距骨组成的踝关节测试结果与Anderson等实验结果对照后发现,二者在中立位踝关节垂直加载下应力主要分布于胫骨下关节面中部及前外侧,最大应力处于2.7 MPa~4.0 MPa之间,在分布区域和分布趋势、数值大小上基本一致,初步证明本模型是有效的。按方法二所建立的简化的两个踝关节三维有限元模型经测试可以模拟踝关节的不同受力情况,并可进行正确的加载运算,而且单元数、节点数、运算时间适中。各种加载后所产生的应力、位移的结果合理,比较贴近真实状况。将结果与Liacouras和Wayne在相同边界条件下所作的踝关节外旋模拟实验结果进行了对比,发现本模型在踝关节外旋加载中胫骨的旋转角度为3.85°,与其有限元模型计算的4.28°比较接近,另外本模型的主要关节的接触力大小也与其结果对应性较好,证明了本模型的有效性。 2.对下胫腓联合损伤模型的研究发现,下胫腓联合损伤后会导致腓距关节的接触力降低,其他关节接触力增高,踝关节部分韧带抗阻力承受的载荷增大,胫骨、腓骨下端位移增大,胫骨的旋转角度(模型中代表踝关节旋转角度)有了明显增大,分别从正常的踝关节内旋时的7.110增加到了12.430,外旋时的3.85°增加到了5.46°,其中踝关节外旋加载结果显示与Liacouras和Wayne关于下胫腓联合损伤模型的研究结果,——胫骨的旋转角度在踝关节外旋模拟加载中从正常4.280增加到了损伤后为5.60°,对比非常接近,说明本损伤模型结果具有有效性。而后对两种螺钉固定模型模拟研究发现,各关节的接触力都比正常有所减少,胫骨、腓骨下端位移减少,胫骨的旋转角度明显缩小,腓骨旋转角度趋向于同胫骨接近,但这两种螺钉固定之间的差异不太明显。 3.对累及损伤不同大小关节面的后踝骨折模型的研究发现,在中立位单足站立时,后踝骨折从损伤1/4踝关节面开始踝关节的最大应力明显增高,而在踝关节内旋加载中踝关节的最大应力是从损伤1/3踝关节面以上明显增高;各骨折模型在中立位单足站立及踝关节内旋加载中均表现为腓距关节的接触力降低,胫距(踝)关节接触力增高;另外,胫骨的旋转角度在踝关节内旋加载时随着后踝骨折损伤关节面的增大角度呈递减趋势,而在踝关节外旋加载中,其角度又随着损伤程度加重而表现为递增现象。 结论: 1.本研究利用Mimics、SolidWorks、ANSYS等软件依据螺旋CT扫描的二维图像构建出了解剖结构精细的正常人体足踝部三维有限元模型,相对客观地反映了足踝部的基本解剖结构和力学特性,模型经验证有效,并可在模型上模拟足踝部不同受力方式进行力学分析。 2.通过对下胫腓联合损伤及两种螺钉内固定的模型结果分析显示,下胫腓联合损伤会导致胫、腓骨下端位移增加,使踝关节在内旋与外旋活动中关节活动度增加,从而引起关节不稳定;而在使用两种位置的下胫腓联合螺钉固定后,均表现为关节活动中各骨位移减少,关节活动受限,关节接触力减小,但两种位置的螺钉之间对比差异性不大,因此,其并非属于符合生理的固定,只适宜短期使用。 3.通过对累及损伤不同大小踝关节面的后踝骨折模型的结果分析表明,尽管小的后踝骨折(损伤1/4以下关节面)在本研究中显示没有明显增加踝关节的最大应力,临床中可考虑保守治疗,但各种后踝骨折都会引起关节接触力及关节活动度发生明显改变。


知网文化
【相似文献】
中国期刊全文数据库 前20条
1 覃立勤;;腓骨中下骨折并下胫腓联合损伤34例的治疗体会[J];广西医学;2011年07期
2 吴轶群;叶晨;张志愿;彭林法;铁瑛;张志勇;王冬梅;张陈平;王成焘;;双侧上颌骨缺损颧种植体修复的有限元分析[J];中国口腔颌面外科杂志;2011年04期
3 贾玉龙;王旭霞;刘舒扬;张君;;牵引上颌埋伏尖牙移动过程的动态模拟[J];临床口腔医学杂志;2011年08期
4 张韶民;李世清;张志伟;王庆生;闫国钦;苗卫东;;下胫腓螺钉在踝关节骨折脱位中的应用46例分析[J];中国误诊学杂志;2011年21期
5 蔚一博;朱强;曹志中;;Micro-CT技术结合逆向工程软件建立上颌第一前磨牙的三维有限元模型[J];第二军医大学学报;2011年07期
6 邹冬华;李正东;黄平;刘宁国;陈忆九;;股骨有限元模型的建立及损伤生物力学验证[J];法医学杂志;2011年04期
7 刘明武;周少怀;金捷;方红育;周雪松;范明宇;;三踝骨折合并下胫腓联合损伤的临床治疗体会[J];中国中医骨伤科杂志;2011年07期
8 刘永青;李杰峰;苟中坤;林序文;余海滨;何启新;钟文辉;赖西铭;;踝关节镜辅助治疗踝关节骨折并下胫腓分离的疗效分析[J];临床医学工程;2011年07期
9 何剑颖;董谢平;;正常髋部与佩戴保护支具髋部的有限元对比分析(英文)[J];中国组织工程研究与临床康复;2011年35期
10 张银光;王岩;王春生;;四桩型胫骨假体基座的优化设计[J];中国骨与关节外科;2011年03期
11 汪光晔;张春才;许硕贵;;正常步态下髋臼底接触面积与压力分布的三维有限元分析[J];中国组织工程研究与临床康复;2011年22期
12 王亚斌;周小建;任亚军;张云庆;;肱骨远端三维有限元模型的建立及生物力学分析[J];中国组织工程研究与临床康复;2011年22期
13 ;中文文摘[J];口腔生物医学;2011年02期
14 魏威;廖胜辉;赖震;;青年颈椎生理曲度变直全节段有限元建模与分析[J];生物医学工程学进展;2011年02期
15 高应超;郭征;付军;田维军;;坐位骨盆的三维有限元分析[J];中国组织工程研究与临床康复;2011年22期
16 刘大军;王燕一;;有限元法在口腔固定修复研究中的应用[J];中华老年口腔医学杂志;2011年02期
17 胡海量;刘利民;沈惠良;;Maisonneuve骨折的诊断与治疗[J];中国骨与关节损伤杂志;2011年06期
18 卢懿明;邹腊春;;关于旋前/后外旋型三踝骨折内固定治疗失败原因分析[J];内蒙古中医药;2011年16期
19 李笑梅;龚璐璐;王少海;汪大林;朱强;;上颌无牙颌修复中种植体部位和数量对应力分布的影响[J];口腔颌面修复学杂志;2011年02期
20 王江华;余雄武;周兴;孙宏斌;;微创与开放结合治疗三踝骨折[J];云南医药;2011年03期
中国重要会议论文全文数据库 前10条
1 苏佳灿;张春才;禹宝庆;许硕贵;王家林;纪方;;骨盆及髋臼三维有限元模型单元选择及构建[A];第十六届全国中西医结合骨伤科学术研讨会暨中西医结合手法治疗骨伤科疾病新进展学习班论文汇编[C];2008年
2 郑杰;杨永宏;张冬生;楼肃亮;廖胜辉;;退变性脊柱侧弯三维有限元模型的建立及有效性验证[A];浙江省医学会骨科学分会30年庆典暨2011年浙江省骨科学学术年会论文汇编[C];2011年
3 姚其卫;牟月照;;下颌双侧游离端可摘义齿三维有限元模型的建立[A];浙江省口腔医学会第二届会员代表大会暨浙江省第二十三届口腔医学学术会议论文汇编[C];2004年
4 肖进;张美超;原林;邱洁;;腰椎运动节段的三维有限元模型[A];解剖学杂志——中国解剖学会2002年年会文摘汇编[C];2002年
5 牛文鑫;杨云峰;俞光荣;丁祖泉;;应用三维有限元模型于跖筋膜切断术的生物力学分析[A];第八届全国生物力学学术会议论文集[C];2006年
6 李晓娜;李振春;赵伟;陈秉智;;下颌单侧游离端太极扣附着体义齿三维有限元模型的建立[A];第六次全国口腔修复学学术会议论文摘要汇编[C];2009年
7 岑远坤;何佳凝;;全下颌牙种植覆盖义齿及其支持组织的三维有限元模型建模方法的改进[A];中华口腔医学会第三次全国口腔修复学术会议论文集[C];1997年
8 孙健;张富强;张秀娟;王冬梅;;个性化下颌体部钛植入支架三维有限元模型的建立[A];2004年上海市口腔医学学术年会论文汇编[C];2004年
9 张健;樊瑜波;;组合式颈椎段整体三维有限元模型的建立[A];中国力学学会学术大会'2005论文摘要集(上)[C];2005年
10 李吉国;倪龙兴;;上颌第一磨牙根管预备后的有限元模型的建立和受力分析[A];2007年第七次全国牙体牙髓病学学术会议论文集[C];2007年
中国博士学位论文全文数据库 前10条
1 刘清华;数字化人体足踝部三维有限元模型的建立及分析[D];南方医科大学;2010年
2 张劲;下颌角截骨整形手术三维有限元模型的建立及生物力学研究[D];南方医科大学;2010年
3 展磊;股骨头坏死三维有限元分析及其与治疗方法选择的研究[D];广州中医药大学;2008年
4 雷勇华;唇腭裂上颌骨三维有限元模型的构建及上颌前牵引生物力学分析[D];中南大学;2008年
5 杨滨;髌股关节整体解剖—力学效能分型初步研究[D];第三军医大学;2009年
6 刘素梅;丹江口水库续建工程诱发地震研究[D];武汉理工大学;2005年
7 张彤;颅面复合体三维有限元模型的建立及上颌前牵引的生物力学研究[D];中国人民解放军军医进修学院;2007年
8 林冬;一个退变颈椎三维有限元模型的建立和应用[D];四川大学;2007年
9 邓锋;口腔正畸微植体支抗技术的生物力学研究[D];四川大学;2006年
10 郭超峰;Lenke1BN型特发性脊柱侧凸有限元模型的建立及后路三维矫形生物力学研究[D];中南大学;2009年
中国硕士学位论文全文数据库 前10条
1 曾季芳;海洋复合材料柔性管的力学性能和优化设计[D];大连理工大学;2009年
2 于书娟;种植体和附着体修复上颌无牙颌的有限元分析[D];第二军医大学;2005年
3 刘梅;不同材料的套筒冠义齿修复单侧游离端缺失基牙位移的三维有限元研究[D];四川大学;2004年
4 张一卓;桥梁伸缩装置疲劳破坏分析及寿命估算[D];长安大学;2005年
5 陈淑萍;Perfect种植体与Replace种植体颈周骨组织应力的比较分析[D];吉林大学;2008年
6 王震东;与微种植体支抗临床应用相关的多因素分析[D];南京医科大学;2005年
7 韩晓辉;不同牙尖斜度全瓷冠的三维有限元分析[D];山东大学;2006年
8 阿兰;骨质、骨量对即刻种植即刻负重影响的生物力学研究[D];吉林大学;2009年
9 赵梦明;Woven玻璃纤维光固化树脂嵌体对松动前牙的固定及三维有限元应力分析[D];天津医科大学;2006年
10 张晓慧;上颌骨Le FortⅠ型截骨牵张成骨的三维有限元分析研究[D];第四军医大学;2005年
中国重要报纸全文数据库 前2条
1 杜精锐;救治小腿毁损伤有了新术式[N];大众卫生报;2001年
2 杜精锐;救治小腿毁损伤有了新术式[N];中国医药报;2001年
中国知网广告投放
 快捷付款方式  订购知网充值卡  订购热线  帮助中心
  • 400-819-9993
  • 010-62982499
  • 010-62783978