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