人体胸廓胸外按压的生物力学测试与有限元分析

论文摘要

一、研究背景及目的猝死在我国的发生率每年大约有60万人。猝死的抢救成功率极低,即使在发达国家,猝死的病人只有5%—15%可以到达医院,而其中只有1%—20%的病人经抢救能够存活。抢救心脏骤停的患者,每延迟一分钟做心肺复苏（Cardiopulmonary Resuscitation,CPR）,其存活率即下降7%-10%。心肺复苏在世界医学领域内所受关注的程度正在逐年上升,现代心肺复苏技术经过近50年的发展,形成了如今国际通用的三级九步法。三级即是心肺复苏的三个基本程序:①基础生命支持（basic life support,BLS）;②高级生命支持（advanced life support,ALS）;③持续生命支持（prolonged life support,PLS）。BLS是整个复苏过程的基础和关键,而胸外按压（external chest compression,ECC）又是BLS中的重要环节,2005国际最新心肺复苏指南特别强调了有效不间断胸外按压的意义。有关胸外按压的频率问题、胸外按压的效果以及胸外按压确切机理等问题一直是CPR学术界关注的焦点:（1）胸外按压频率由早年推荐的60次/分钟,增加到最新要求的100次/分钟,近年来有学者提出了高频按压（120次/min以上）的观点,是否频率越高越好呢?（2）胸外按压的有效性问题:最新心肺复苏指南对胸外按压的部位、力度、下压的深度有了明确的规定,并强调按压应有力而快速:按压部位要求在两乳头之间,按压幅度为4～5cm,婴儿为胸廓的1/3～1/2;每次按压后胸廓完全弹回,保证松开的时间与按下基本相同。然而有临床研究显示:临床胸外按压过程中,近50%的胸外按压不合格:一半的胸外按压幅度太浅,而且CPR过程中24%～49%的时间里未进行胸外按压。按压时,允许胸壁弹性回缩能够增加回心血量,而胸壁复位不充分时,回心血量减少导致心排血量降低;过多中断按压,冠脉和脑血流中断,复苏成功率明显降低。（3）胸外按压的机制问题:最早提出的是心泵机制学说,即认为有节律地按压胸骨可使胸骨与脊柱间的心脏被挤压,关闭房室瓣使心室内压增高,推动血流向前;按压解除时,心室恢复舒张状态产生吸引作用,使血流充盈心脏,反复按压推动血液流动而建立人工循环。20世纪80年代,学者们开始从血流机制方面对CPR机制进行探讨,提出了胸泵机制学说,即心肺复苏时胸外按压使心脏复跳的机制是通过胸腔压力的变化而使血流流动。近年来,国内有作者提出胸腔-大血管-心脏之间的“共振频率”效应学说,究竟哪种学说是胸外按压时血流向前的主要机制?以上三个与临床紧密结合的问题的研究,目前在临床上存在观察例数少、缺乏对比实验、病例选择受限、影响因素较多等限制,较难得出科学的结论。心肺复苏传统研究方法有:动物模型研究、临床研究和人体模型研究。能否找到一种研究方法,既可以对胸外按压的效果作更深入的探讨,还可以进一步研究传统认可的“心泵”、“胸泵”机制和推论的“共振”机制,甚至发现新的胸外按压机制?本研究通过人体标本的生物力学试验,测试出人体胸外按压时作用力与胸廓位移及应变的关系;利用生物力学三维有限元建模的方法,在计算机上重建正常人体胸廓三维有限元模型,并对模型进行有效性检验后,应用有限元分析方法模拟分析外力作用下胸廓的应力及应变分布。本方法避免了以往心肺复苏研究的局限性,首次应用生物力学的研究方法分析心肺复苏胸外按压时按压力度、频率与形变关系等胸外按压的关键技术问题,为CPR的基础研究开辟了一个新的探讨领域。二、材料与方法1.胸外按压的解剖学基础:了解人体胸廓、胸壁,以及心脏的解剖结构。2.人体胸廓胸外按压的生物力学测试试验:选择1具男性新鲜尸体标本,实验前摄片以排除胸廓原有疾病（包括骨折,畸形等）的存在。应用解放军第458医院157临床部双排螺旋CT对其进行胸部扫描。采用已知重力0-200牛顿（Newt,N）,分别模拟人体胸外按压时胸廓所承受的负重工况,使用MTS材料试验机和引伸仪,测试垂直加压情况下人体胸廓的位移和应变。本实验分为A、B两组。A组是在具有完整胸部的标本上完成试验;B组是在A组试验完成后,剔除标本胸壁的肌肉和软组织以及胸腔脏器后进行试验。测试方法设计为静态垂直加压和动态加压;静态垂直加压方法指设定加载范围为0-200 N,实施分级加载,每级为50 N。试验机的加载速度为5.0mm/S。动态加压方法指模拟临床胸外按压条件,设置频率为100次/分钟、施加载荷的时间与胸廓放松时间相等（0.3s）的情况下,测试胸廓按压点在不同外力作用下的位移。A、B两组试验位移的测量点相同:按压点（第7测量点）—胸骨体前正中线第4、5肋间;第6测量点—左侧第5肋距前正中线38.0mm;第5测量点—左侧第4肋距前正中线42.0mm;第4测量点—左侧第3肋距前正中线37.6mm;第3测量点—左侧第2肋距前正中线28.7mm;第2测量点—前正中线第3肋水平面的交界点;第1测量点—前正中线第2肋水平面的交界点。3.构建人体胸廓三维模型:选择与前试验同1例的男性新鲜尸体进行螺旋CT扫描,扫描范围:第一胸椎上缘至第十二肋骨下缘进行连续螺旋扫描,选择骨组织窗的CT影像窗位为396,窗宽为1536。扫描层厚1.0mm,扫描时间为36.53 s,共获取人体胸廓二维图像440层,图像以DICOM格式存贮。将二维DICOM格式的图像导入Mimics10.1软件,确定图像空间方位后,对图像进行分割,建立了包括脊柱、胸骨和肋骨的人体骨性胸廓三维模型。4.人体胸廓三维有限元模型的建立:在Mimics10.1软件中,将人体胸廓进行表面网格划分并进行优化后将模型以Ansys Element Files（.lis）的格式导入ANSYS9.0软件中进行体网格划分,并保存胸廓体网格模型文件返回导入Mimics10.1软件。在Mimics10.1的有限元分析模块中,利用软件自带的“CT图像灰度-组织力学材料性质”关系公式,对三维胸廓体网格模型的各单元进行材质分配。材质分配完成后的体网格模型再次返回导入ANSYS9.0软件中进一步优化,如此多次反复、优化后得出完整人体胸廓有限元网格模型。5.人体胸廓三维有限元模型的有效性检验:将前期生物力学测试的结果与已经建立的胸廓有限元模型进行对照,比较相同作用力下两组的位移、应变关系,通过两组数据的趋势曲线来检验模型的有效性。6.人体胸廓胸外按压的有限元分析:将证实有效的人体胸廓三维有限元模型进行垂直方向模拟按压加载,约束边界条件后进行求解运算,分析压力为0-400N时人体胸廓有限元模型各部位的位移、应力和应变分布,结合临床进行结果分析。三、结果1.人体胸廓具有一定的弹性和活动性,胸廓各骨通过肋椎关节、肋横突关节、胸肋关节相互连接,当胸外按压时,垂直向下的力量作用于胸骨,再通过胸骨与肋骨的连接、肋骨和肋椎关节等使胸廓的矢状径和横径发生变化。肋软骨富于弹性,在抢救呼吸骤停的病人时,可进行体外心脏按压或人工呼吸。人体胸外按压时按压的部位即心脏与胸廓最为接近的位置。2.人体胸廓标本的生物力学测试显示:（1）静态垂直加压时,由A组和B组的力-位移关系曲线可知:①压力由0～200 N变化时,位移呈现压力依赖关系,即随着压力的增加,胸廓各测量点的位移也随之增加。②压力相同时,不同测量点的位移不同,按压点（7点）的位移最大,依次是7点＞6点＞5点＞2点＞4点＞1点＞3点。③从A组的压力-位移关系曲线上看到,胸廓下移达到40mm距离的只有按压点在受到200N的载荷力时才能实现。④A、B两组比较:同力作用下B组位移更大,当外力为200N时,B组按压点位移达到60mm,观察发现标本的胸骨和肋骨未发生骨折。（2）动态加压时胸廓的压力-位移关系曲线结果:当模拟临床胸外按压频率为100次/分钟、按压与放松时间相等时测量到的外力与位移的关系从趋势线上得知,位移随着外力的增加而增大。（3）从动态和静态加压时位移的柱形图可知:相同测量点的动态位移明显小于静态位移。根据力-应变关系曲线,当外力作用于胸廓时,应变呈现压力依赖关系,即随着外力的增大,各测量点的应变也增大。3.建立了包含1158 085个节点、736 022个单元的计算机人体胸廓三维有限元模型,此模型完整地重现了人体胸廓的复杂形态,得到了包括锁骨、肋骨、脊柱等结构鲜明、直观、整体的可计算模型图像。4.由于完整胸廓的单元和网格数据偏大,而且考虑到第11、12浮肋与前胸廓不相连,对胸廓运动的影响不大可忽略等因素,计算前减去第11肋以下部分,得到包含356 562个节点、215 808个单元的胸廓模型。通过比较生物力学测试试验与所建立模型在相同力作用下的位移和应变,我们发现,首次建立的模型在外力作用下的位移与力学测试的位移差别较大,重新调整模型的材料性质后将模型与力学测试再次进行比较,最后得出较为合理的胸廓可计算模型。5.胸廓有限元模型分析结果:正常人体胸廓施加压力后的应力分布均匀,呈对称性分布。按压胸廓时,胸廓的应变集中部位与应力集中部位不同,前者多位于胸骨、肋骨交界处,后者多位于肋骨约束自由度处;即变形程度最大的部位在按压部位的胸骨、肋骨交界处,而受力的集中部位在背部。四、结论1.人体胸腔内心脏与胸廓最接近的位置与胸外按压的位置几乎相同,此结构特征是“心泵机制”理论的解剖基础。胸廓是一个由多种骨和骨连接组成的复杂结构,胸外按压时,垂直向下的力量作用于胸骨,再通过胸骨与肋骨的连接、肋骨、肋椎关节等使胸廓的矢状径和横径发生变化,胸腔内压增大,这是“胸泵机制”理论中胸腔压力变化的解剖学基础。2.胸外按压时,在相同力作用下,按压点的位移最大;距离按压点位置越近,位移越大;在同一肋骨水平,位于胸骨的测量点较肋骨测量点的位移大。按压胸廓时,按压力度越大,胸廓的变形也越大,由此推论胸腔压力的改变也会随之增加。3.在同力按压时,完整胸腔的下压位移小于骨性胸廓的下压位移,由于下压点位置正好是心脏与胸骨接触面最大、最密切的位置,心脏的相反作用力可使按压时的动能减少导致位移减小,从而推论:胸外按压时心脏的形变是存在的。4.模拟人体胸外按压频率的生物力学动态测试试验得知作用力与位移的关系曲线并得出结论:随着按压力的增加,位移逐渐加大;从曲线的趋势公式:Y=8.0173X+38.698（Y为作用力,X为位移）得知:当位移达到40mm时,施加的外力是359.39N（相当于37公斤）,其胸廓的位移与静态测试人体标本以及有限元模型计算出的位移比较都明显减少。根据实验结果可以推论,使胸廓下压40mm需要达到至少300N（相当于31公斤）的按压力量,频率越高,需要的下压力量也越大。按压力量的大致明确对于我们今后的胸外按压实际操作具有重大的意义。5.本研究成功建立了人体胸廓三维有限元模型,此模型的建立在生物力学有限元研究领域和急救医学领域均属创新。6.本研究将建立的人体胸廓有限元模型进行了有效性检验,将修正后的模型计算结果与生物力学测试数据进行对照,在相同外力作用下,两组的位移曲线趋势一致,应变曲线趋势也一致而且极为接近,由此我们认为本模型有效,并且能为我们今后的临床研究提供有意义的有限元分析数据。7.从胸廓有限元模型应力、应变的分析结果得知,胸外按压时应变的集中部位在胸骨、肋骨交界处（即肋软骨处）,而各个肋软骨的应力和应变最大部位在第4肋软骨和第5肋软骨处;从人体标本生物力学测试的位移变化来看,各个肋软骨位移变化由大到小依次为:第5肋软骨＞第4肋软骨＞第3肋软骨＞第2肋软骨;结合人体胸廓解剖结构可以推论:第5肋软骨是临床按压时最容易发生骨折的部位,胸外按压的部位确定在第4肋间隙与胸骨的交界处较为合理。8.本研究不足之处和未来展望。本研究在建立人体胸廓有限元模型时所涉及的生物材料力学特性均假定为均质、连续和各向同性,然而实物材料本身并不是均质、连续的,也不是各向同性,而是呈各向异性的特征。对于这种状况,绝大多数有关骨骼的有限元分析模型都是存在的,这也是该方法本身固有的缺陷。人体胸廓具有胸肋关节、肋椎关节、和肋横突关节,然而本研究建立的胸廓模型不具有以上关节,特别是胸肋关节的缺失,在一定程度上影响了胸廓位移的变化。以往有限元有关骨骼的研究多为静态情况下的局部分析,然而胸外按压是一个动态进行的过程,由于人体胸廓有限元模型数据偏大和有限元分析软件功能的限制,本实验目前还未能在胸廓有限元模型上动态分析有关作用力与位移、频率与胸腔容积改变的关系,以及利用有限元模型进行胸廓固有频率的计算,预计在不远的将来,随着有限元分析技术的不断完善,我们的研究会不断地完善,此模型也能更好地为临床服务。本研究建立的人体胸廓有限元模型的分析和应用,不仅局限在急救复苏的研究领域,还可以用于其他人体胸廓生物力学相关问题的研究,其应用范围很广,可重复性强。此模型的分析可以设计为相关分析软件作进一步的推广和应用,具有一定的社会效益和经济效益。通过研究胸外按压时作用力与位移的关系,为我们今后人体胸外按压的临床研究打下良好的基础。

论文目录

• 摘要
• ABSTRACT
• 前言
• 第一章 人体胸廓的生物力学测试研究
• 1.1 胸外按压的解剖学基础
• 1.2 人体胸廓胸外按压的生物力学测试试验
• 第二章 人体胸廓三维有限元模型的建立
• 2.1 人体胸廓CT图像的提取及构建三维模型
• 2.2 人体胸廓三维有限元模型的建立
• 第三章 人体胸廓三维有限元模型的有效性检验
• 1 材料
• 2 方法与结果
• 3 讨论
• 第四章 人体胸廓胸外按压的有限元分析
• 1 材料
• 2 方法
• 3 结果
• 4 讨论
• 全文总结
• 本研究不足之处和未来展望
• 参考文献
• 综述
• 附录
• 攻读学位期间成果
• 致谢