有限元脚部模型的简化.docVIP

有限元脚部模型的简化 　　摘 要：为了缩短患足底腱膜或腱膜炎的病人的康复时间，医生们的制鞋过程，使研发人员更快速地得知脚部受力分析，用有限元分析软件建立脚部模型非常必要。现有建模方法需要插入八万个节点，严格参阅解剖学知识找到八万个节点中对应的两节点进行相连，用单元模拟软骨和韧带。而通过线性化趾骨和跖骨、在胫骨和腓骨上端加入近似刚体圆盘来施加受力和限制自由度、结合MATLAB或无前端显示只需输入文件的有限元软件，可以更快得到精确的受力分布图。 　　关键词：有限元；足部；跖骨；趾骨；胫骨；腓骨 　　一直以来，国内外都有大量饱受足部病痛的患者。在这之中，超过七成的病痛都是足底筋膜炎和腱膜炎导致的。[1] 国外现有的治疗方法，是由医生用患者踩踏软泥得到的脚型，浇注成模，再制造出完全贴合脚型的鞋子，来顶起足弓处，分担受力，使筋膜、腱膜缓慢康复。但此制鞋过程需耗费一个月时间、且价格昂贵。[2] 于是，全新的制鞋理念诞生：在鞋内加入刚度和高度调节装置，配合贴合脚底特征的鞋垫即可。为了更好地设计与落实此理念，需要对脚部进行有限元分析，建立脚部模型，施加受力和位置限制，了解足部各处的受力分布，从而确定两调节器的放设位置。同时，也可为医生治疗提供更科学的受力分析情况。 　　一、现有模型的问题 　　为了得到脚部模型，需要依次实行如下步骤：用3D扫描仪得到脚部特征，将其修复改进成点云集合，导成有限元软件可以导入的模型格式，并且利用线性或非线性数学原理增加节点数，以填满所有点云间的空隙。[3] 　　在过程中的最后一步，建模人员总共需要设置八万五千个节点来建模，工作量十分巨大。而且为了确保精确性，更好的方法是在inp文件中手动修改。[4] 纵然节点的坐标可以在其他软件中由程序得到，一来一回的操作也是十分复杂。这其中稍有运算上的漏洞，就要推倒重来。更为复杂的是，建模的下一步，是要模拟软骨和韧带。两者分布在脚上，可以用设置单元，将两节点连接起来的方式来得到。但为了确保精确性，研究人员需要仔细研读解剖学的书，找出八万个节点中精确的两点，显然非常费时费力。以上的步骤，对每一个患者，都要再来一次，工作量可见一斑。虽然不需医生花一个月时间进行手工制作，但作为此创意鞋的先行步骤，有限元分析的建模就要耗费如此精力，对医生而言，不切实际也无法凭自己能力完成，不利于鞋子的下一步推广。 　　二、解决方案 　　（一）基本概念和简化方法 　　最主要的方法是将五根趾骨和五根跖骨进行线性化替代。显然的是，软骨和韧带都是由代表肉质的两个节点连接而成，与骨头所代表的节点并无关联，故进行线性化，减少的节点不会对受力分析的结果有直接影响。其二，五根趾骨和五根跖骨在自身的整个长度中除了截面面积的变化外，别的性质，例如骨头密度、杨氏模量都未发生变化；惯量、转动惯量则都是随截面面积这一自变量进行变化。故只要在线性化的过程中体现出截面面积的变化，即可大致进行等量转换。[5] 　　值得注意的是，在跖骨和趾骨交界处的籽骨，必须予以保留。借此，两骨与脚部整体的相对性得以保留。当胫骨和腓骨所受的力向下传递至脚跟及脚后部的骰骨、舟骨、跟骨、距骨时，五根跖骨因为都与四骨之一及籽骨之间连接，受力情况可保持与简化前模型相同。[6] 另外，跖骨和趾骨分别用两个方向的弹簧和籽骨相连，模拟出平动和转动两种相对位置的变化。 　　第二个可以简化的是最上端的胫骨和腓骨。许多同领域的前期研究，都基本一致得认为分析受力，需要保留整个脚部直至胫骨和腓骨--脚踝处的距骨与腓骨、胫骨、舟骨都通过软骨有接触，此处受力最为复杂，故必须集体保留腓骨和胫骨才能得到准确结果。[7] 更进一步来说，受力必是施加在胫骨和腓骨的上端。但如何建模更符合实际是需要思考的。鉴于此模型要模拟静止和走路，两种步态的受力，而模型中却无一点是限制了自由度的，在分析软件中一定会得出错误的答案。所以，首次尝试是将腓骨和胫骨各自连接一根杆，向上延伸至同一圆盘，此圆盘的面积大抵为胫骨和腓骨的面积之和，而它的杨氏模量则设为一个极大的值。这样，他可以被近似为一个刚体，很难有移动。其面积的大小在此情况下也不再重要，设为一个合理值即可。另外，胫骨和腓骨各自连接的杆，其杨氏模量与胫骨或腓骨相同。但这样一来，施加人自身重力对脚部的压力的一半，四百五十牛，于每只脚，得到的结果依然呈离散性。经过分析，现有模型，是将膝盖处设为在水平面上两方向无自由度，在竖直方向给予限制的自由度，约为0.01mm。但由于膝盖至胫骨及腓骨的距离较长；即两杆距离过长，而两杆本身又无自由度的限制，导入对整个脚部，依然是没有边界条件的限制。[8] 因此，最后在胫骨和腓骨自身最上端，沿着同样的特征形状向上拓高5mm，建立一个高为5mm的柱体，将三自由度限制和力加在此柱体上。此柱体的杨氏模量设一极大