基干混合有限元心肌纤维动力学仿真方法.docVIP

基于混合有限元的心肌纤维动力学仿真方法 　　摘 要：心脏的节律性收缩是其泵血功能的动力学基础，研究其收缩行为可以为心脏循环系统运转提供科学指导。心脏的收缩是心脏被动力和主动力共同作用的结果。被动力行为与心肌组织的材料性质有关，主动力行为心脏受到电生理调控。本文基于混合有限元方法建立心肌纤维弹性力学模型，仿真实验表明该模型可以有效仿真心脏组织在体力、牵引力和主动力作用下的形态学变化。综上本文提出的方法可以为今后建立电力耦合模型和心脏动力学研究提供理论指导 关键词：心肌纤维力学建模；有限弹性力学；混合有限元 中图分类号：TP391.41 文献标识码：A 文章编号：2095-2163（2015）05- Abstract： The power of heart pumping origins from rhythmic contraction of the myocardial tissue. Research of heart contraction is benefit to understand the circulation system. Passive property of heart， which is determined by materials features of myocardial issue， in combination with the active force instigated by electrical activity that leads to its contraction. In this paper， it provides that based on mixed finite element methods myocardial fiber mechanic model is developed to simulate elastic features of heart contraction. The simulated results suggest that this model can effectively simulate mechanic actions of heart under body force， traction and active force. The method provide a theoretical guide to set up electro-mechanic coupling model of heart and cardiac mechanic behaviors. Keywords： Myocardial Fiber Mechanic Model； Finite Elasticity； Mixed Finite Element Methods 0引 言 心脏的计算模型可以分为电生理模型和力学模型两大类。心脏电生理模型的研究始于Noble[1]，其中是将Hodgkin[2]关于电兴奋在枪乌贼神经上传导的模型运用到心脏的浦氏纤维中。此后，心脏电生理建模取得了长足的发展。然而，心脏的力学模型研究滞后于电生理模型[3-4]，在心脏力学建模领域存在着诸多尚未解决的难题。目前，国外的文献虽然提出了心脏力学的建模方法，但很多文献存在不一致和不准确的地方，比如在处理不可压缩约束时，某些文献在应力表达式加入了项 ，有些文献加入了项 ；有些文献将弱形式和基于最小能量泛函的变分形式混用。另外，在模型中考虑心脏的不可压缩性[5]会给模型的数值求解带来难度，研究人员主要通过罚函数法[6-8]和多域有限元法[9，10]来解决该问题。罚函数法主要用于基于位移的有限元框架中，其形式简单，只需要计算位移变量，然而在材料趋向于不可压缩时，该方法会出现锁定现象[11]。多域有限元方法能够很好地避免锁定现象，对心脏不可压缩性质进行良好建模。基于以上考虑，本文首先严格按照有限弹性力学的方法建立了心脏的力学模型，然后给出了基于混合有限元方法的模型求解流程，最后通过实验对方法进行了测试 1基于有限弹性力学的心脏力学模型 心脏由左右心房、左右心室四个腔室构成，左心室是四个腔室中体积最大的，左心室收缩时产生的高压将血液泵送至左心房[12]。为了能承受高压，左心室的心室壁比较厚，分别由心内层，心肌层和心外层构成。本文的力学模型主要针对左心室，并在左心室的一个六面体形状的切块上进行仿真分析 图1为心脏切块形变示意图，用 表示未发生形变的心脏切块，其上的点记为 ，称该坐标为材料坐标[13]；用 表示发生形变后的心脏切块，其上的点记为 ，称该坐标为空间坐标[13]。将初始构形中的质点 映射到当前构形中的质点 的函数称为形变函数： 接着