生物力学论文

研究生课程（论文类）试题2014/2015学年第1学期课程名称：生物力学课程代码：19000011论文题目：气囊减重步行训练有限元分析学生姓名：专业﹑学号：生物医学工程学院：医疗器械与食品学院课程（论文）成绩：课程（论文）评分依据（必填）：任课教师签字：日期：年月日课程（论文）题目：气囊减重步行训练有限元分析内容：1.引言随着社会的发展，人类文明的不断进步，人们的生活质量也越来越好，但是，人们也不得不面临越来越多的疾病困扰。最为显著的是下肢功能障。众所周知，步行是人们日常生活中最基本的功能活动之一，但是下肢功能障碍使得患者丧失这一基本功能，给生活带来了极大的不便，在加重了患者本身的不便的同时也给患者的家庭加重了负担，据统计，患有致残的神经系统疾病的患者如脑卒中、脊髓损伤、外伤、格林-巴利等病后，3年能恢复独立行走50m的患者仅有3%—10%。对许多因各种原因导致暂时不能行走的患者来说，恢复步行这一基本功能是他们的第一愿望，这也是康复医学研究的一个主要内容。减重步行训练是针对下肢功能障碍，改善步行能力的一种新的康复治疗技术。躯干和下肢承重能力下降是导致步行不能的重要原因。传统康复治疗已经采用减重的方式进行早期步行训练，例如利用水的浮力进行水中步行，利用各类拐杖或者步行器减少下肢负重等。但都是存在一些不足，例如在水中运动需要特殊环境，拐杖或者助行器需要患者增加上肢用力，造成步行时身体姿态异常。减重训练是以传统实践为依据，利用悬吊装置不同程度地减少上身体重对下肢的负荷，在理论上有利于支撑能力不足的患者早期进行各种步行训练。目前减重训练已经开始临床应用研究。2.减重步行训练的发展及应用现状现代减重步行训练起源于上世纪80年代脊髓损伤的猫动物模型的步行训练。目前较为认可的理论是“中枢模式发生器”（Centralpatterngenerator，CPG）学说。1985年Finch等.在电动跑台上对脊髓横断猫后腿减重后被动迈步的试验中，发现经过一段时间的训练．可使猫重新获得运动能力：并在进一步研究中发现．其运动能力的恢复与脊髓非完全损伤尚存在特定节段感觉输入机制有关。这种踏板运动有赖于特定的节段感觉传人，类似于人类腰骶运动区的中枢模式发生器.CPG是由脊髓中神经元相互联结组成的网络，能发放节律冲动，肌肉的节律性收缩就由低级中枢控制。也有人提出CPG存在于脊髓的腹侧和中部的两侧，它们相互之间存在信号通讯，以脊髓颈和腰膨大处最多，通过神经环路与其他神经关联。当特定感觉传人后，CPG产生步行中屈肌和肌交替转换的神经冲动，首先屈肌兴奋性冲动通过中间神经元抑制伸肌活动，屈肌兴奋完成后伸肌神经兴奋释放，引起伸肌活动，从而在步行动作启动之后，产生自发性屈肌一伸肌交替兴奋，产生迈步动作。传统步行运动治疗强调诱发下肢关节分离运动，单独训练迈步、平衡、重心转移等分解动作。但真实步行环境与此不同，结果常导致患者分离运动好，步行能力差的矛盾现象。减重步行训练（BWSTT）强调在真实步行环境中进行综合性训练，注重实用性，增加安全性，减少能量消耗，扩大活动范围，提高生活质量。自上世纪80年代始，BWSTT开始临床应用于脑卒中及脊髓损伤患者。患者下肢无充分负重能力时即可开始BWSTT，而无需等训练到能步行后才进行传统步态训练，这显然符合脑卒中早期康复理念。因此开始在脑卒中临床方面开始应用。减重训练的临床应用可以追溯到1958年，Margaret等出版了专著“康复治疗中的悬吊疗法”。但是由于方法的局限和认识不足，没有得到发展。将减重训练用于神经瘫痪患者的新热潮始于加拿大学者Visintin等1989年的报道，他们发现痉挛性瘫痪者进行40%减重训练活动平板训练6周后，平衡功能、步行速度和步行耐力均显著高于常规训练组，随访3个月时训练组的步行速度和运动恢复得分进一步提高。Pillar等报道24例研究对象，包括偏瘫、脊髓损伤和膝下截肢，与6名正常人相对照。减重系统为固定在天花板的滑轨和选吊带，减重范围0-150kg。采用录像分析系统采集支撑相和摆动相的时相，分析步态的对称性和时速，发现受累肢体的支撑相时间减少，对称性改善，步速增加，过去不能步行者现在能够步行；治疗师可以集中精力关注受累肢体，促进步态改善。减重步行练系统由两部分组成，即减重装置(PBWS)和电动活动平板(treadmill)。减重装置(电动)主要包括固定支撑架、减重控制台、电动升降杆、减重吊带。减重控制台控制电动升降杆的升降，随着升降杆的升高，患者被逐渐向上吊起．下肢负重减少．减少的重量可以在减重控制台上显示出来。治疗师可以按需要从下肢O％(完全负重)一100％(完全不负重)调整下肢减重量。目前，初始的减重量大多采用减掉患者身体质量的30％一加％，也有减重60％的报道。训练前需满足两个条件：①患者负重达到可能支撑的最大体重；(函患者髋关节能够完全伸展。一旦进步，尽快减少减重量．直至达到全负重．但必须保证正确的步态模式及安全性。减重吊带类似于降落伞固定带，使用时需紧紧缚于患者的腰臀部．固定带的两端对称固定在悬吊支撑架上。活动平板(电动)用于减重患者的步行训练．平板运行时间、速度和坡度可以根据需要进行调节．一般初始速度设定为0.1-0.5m/s。每次步行训练30--40min．采取间歇训练法，间歇频率及时间因人而异，主要根据患者的耐受程度和疲劳恢复情况。由治疗师具体掌握，以后根据患者恢复情况逐步减少间歇次数、间歇时间，增加平板运动速度。近年来，减重步行训练在国内得到广泛的应用。但是减重步行训练在国内所涉及的病种尚不够深入．使用时需多名治疗师参与治疗。人力投入相对较大．限制了它在临床工作中的使用。尽管有关减重步行训练的研究报道较多，但缺乏大样本。多中心对照研究，对开始减重的时机、减重量、平板运动速度、持续间歇时间、疗程等方面。各家报道不一，还有待进一步循证优化。以确定最佳训练方案。3.减重步行训练的有限元分析3.1有限元的基本理论有限元法是结构分析的一种数值计算方法。它在20世纪50年代初期随着计算机的发展应运而生。最早应用于结构力学，后来随着计算机的发展慢慢用于流体力学的数值模拟。有限元分析的基本概念是用较简单的问题代替复杂问题后再求解。它将求解域看成是由许多称为有限元的小的互连子域组成,对每一单元假定一个合适的近似解,然后推导求解这个域的满足条件,从而得到问题的解。这个解不是准确解,而是近似解,因为实际问题被较简单的问题所代替。由于大多数实际问题难以得到准确解,而有限元不仅计算精度高,而且能适应各种复杂形状,因而成为行之有效的工程分析手段。有限元是那些集合在一起能够表示实际连续域的离散单元。有限元的概念早在几个世纪前就已产生并得到了应用,例如用多边形逼近圆来求得圆的周长,但作为一种方法而被提出,则是最近的事。有限元法最初被称为矩阵近似方法,应用于航空器的结构强度计算,并由于其方便性、实用性和有效性而引起从事力学研究的科学家的浓厚兴趣。经过短短数十年的努力,随着计算机技术的快速发展和普及,有限元方法迅速从结构工程强度分析计算扩展到几乎所有的科学技术领域,成为一种丰富多彩、应用广泛并且实用高效的数值分析方法。有限元方法与其他求解边值问题近似方法的根本区别在于它的近似性仅限于相对小的子域中。20世纪60年代初首次提出结构力学计算有限元概念的克拉夫(Clough)教授形象地将其描绘为:“有限元法=RayleighRitz法+分片函数”。将函数定义在简单几何形状的单元域上,且不考虑整个定义域的复杂边界条件,这是有限元法优于其他近似方法的原因之一。在有限元方法中，把计算域离散剖分为有限个互不重叠且相互连接的单元，在每个单元内选择基函数，用单元基函数的线形组合来逼近单元中的真解，整个计算域上总体的基函数可以看为由每个单元基函数组成的，则整个计算域内的解可以看作是由所有单元上的近似解构成。在河道数值模拟中，常见的有限元计算方法是由变分法和加权余量法发展而来的里兹法和伽辽金法、最小二乘法等。根据所采用的权函数和插值函数的不同，有限元方法也分为多种计算格式。从权函数的选择来说，有配置法、矩量法、最小二乘法和伽辽金法，从计算单元网格的形状来划分，有三角形网格、四边形网格和多边形网格，从插值函数的精度来划分，又分为线性插值函数和高次插值函数等。不同的组合同样构成不同的有限元计算格式。对于权函数，伽辽金(Galerkin)法是将权函数取为逼近函数中的基函数；最小二乘法是令权函数等于余量本身，而内积的极小值则为对代求系数的平方误差最小；在配置法中，先在计算域内选取N个配置点。令近似解在选定的N个配置点上严格满足微分方程，即在配置点上令方程余量为0。插值函数一般由不同次幂的多项式组成，但也有采用三角函数或指数函数组成的乘积表示，但最常用的多项式插值函数。有限元插值函数分为两大类，一类只要求插值多项式本身在插值点取已知值，称为拉格朗日(Lagrange)多项式插值；另一种不仅要求插值多项式本身，还要求它的导数值在插值点取已知值，称为哈密特(Hermite)多项式插值。单元坐标有笛卡尔直角坐标系和无因次自然坐标，有对称和不对称等。常采用的无因次坐标是一种局部坐标系，它的定义取决于单元的几何形状，一维看作长度比，二维看作面积比，三维看作体积比。在二维有限元中，三角形单元应用的最早，近来四边形等参元的应用也越来越广。对于二维三角形和四边形电源单元，常采用的插值函数为有Lagrange插值直角坐标系中的线性插值函数及二阶或更高阶插值函数、面积坐标系中的线性插值函数、二阶或更高阶插值函数等。有限元法的理论基础牢靠，物理概念清晰，解题效率高，目前已成为机械产品动、静、热特性分析的重要手段。3.2有限元分析的基本步骤第一步、创建有限元模型根据我所选的毕业设计的题目以及参考实物，我运用SolidWorks软件画出所要分析的底板的立体几何图（如图3.3所示）。第二步、打开ANSYSWorkbench主界面创建一个新的分析项目打开workbench主界面双击Geometry选项，并用鼠标单击StaticStructural选项，并拖拽至A2项目下生成一个关联分析项目B（如下图3.1所示）。图3.1创建新的分析项目第三步、导入实体模型在A2栏的Geometry上单击鼠标右键，在弹出的快捷菜单中选择ImportGeometry→Browse命令，选择所绘制的底板的立体几何图。此时，A2栏Geometry后的变为√，表示实体模型已经存在。第四步、添加材料属性在Workbench主界面空白区域单击鼠标右键选择EngineeringData选项，进入材料库，在这里可以选择添加一些材料的属性。第五步、划分网格根据实际情况和有限元分析的条件及要求选择适当的网格划分方式，在本文中我选取的网格划分目标的物理环境为结构分析（Mechanical）。第六步、施加载荷与约束根据课题的要求和实际情况的需要选择适当的载荷和约束力。第七步、结构后处理在此步骤中可以添加等效应力等效应变以及总变形选项。第八步、保存与退出单击Mechanical界面右上角关闭按钮，退出Mechanical返回到Workbench主界面，此时主界面中的项目管理区中显示的分析项目均已完成。在Workbench主界面中单击工具栏中的保存按钮，保存含有分析结果的文件。3.3底板的结构设计及介绍在做结构设计的时候我尝试着做了两种不同的结构设计图：一种是一体化整体结构（如下图3.2）图3.2地板结构还有一种是由五个不同的单一零件装配成的装配体的结构（如下图所示）。下面是组成装配体的不同的零件的结构设计。图3.2.1零件一图3.2.2零件二图3.2.3零件三图3.2.4零件四图3.2.5零件五图3.3地板结构比较一下这两种设计的不同以及在之后的工作中遇到的一些问题。首先，在立体效果方面，两幅图大概差不多，都能表达减重步行训练装置的底板的效果。其次，在用SolidWorks画图上面的不同。图一是先确定一基准面好后。画一个长方形，直接进行凸台拉伸，形成一个整体的六面体，然后再针对各个部分进行切除拉伸，最后形成一个如图