机械制造实习论文范文-功率键合图方法在多分支血液循环系统建模

1、机械制造实习论文范文-功率键合图方法在多分支血液循环系统建模和仿真中的应用    摘要：应用功率键合图方法，建立了一种多分支血液循环系统的计算机仿真模型，即描述血液循环系统内血流动力学变量变化规律的状态方程。该仿真模型较为细致地刻画了血液循环系统的生理特性，形成了较完整的人体血液循环系统的计算机模型，此模型可模拟血液循环系统的生理和病理特性，得出相应的心血管动力学仿真数据和波形，为进行血液循环系统生理和病理的医学研究提供了新的研究手段。关键词：血液循环系统计算机仿真功率键合图法0 引言 功率键合图法是一种系统动力学建模方法，它以图形方法来表示、描述系统动

2、态结构，是对流体系统进行动态数字仿真时有效的建模工具。通过已有的研究工作表明，功率键合图方法可以较好地应用于生物流体系统仿真，特别是人体循环系统的建模和数字仿真10。 我们在以前的工作当中，建立了一个简化的血液循环系统模型10，验证了功率键合图法的可行性和有效性。键合图建模方法的优点是直观形象，便于获得状态空间方程，有利于数值化计算，避免了电模拟方法中推导状态方程困难的弱点 。本文对血液循环系统进行了较细致和全面的划分，建立了一个包括动脉系统、静脉系统、心脏（左、右心室和心房）以及冠脉循环、外周循环的多分支血液循环系统仿真模型。 应用功率键合图方法对血液循环系统进行建模和仿真的基本规则是，(1

3、)把血液循环系统的结构及各主要动态影响因素以图示模型形式，即功率键合图加以表示，(2)从功率键合图出发，建立系统的动态数学模型状态空间方程，(3)在数字计算机上对状态方程进行求解。   1 多分支血液循环系统模型的建立   1.1 系统描述 血液循环系统模型如图所示4。在心血管循环系统中，血液在心脏“泵”的作用下所进行的循环流动，可以看作是一种功率流的流动、传输、分配和转换的过程。血液在左右心室有节律地收缩作用下，被泵向人体的各个部分，其中包括：体循环区（血液由左心室经主动脉、大动脉、外周循环区和腔静脉，回到右心房），肺循环区（血液由右心室流经肺动脉和肺静脉到左心房。），腹

4、部内循环，颈部和头部循环，以及冠脉循环等。在心房和心室、心室和主动脉之间存在着防止血液倒流的膜瓣，如二尖瓣、三尖瓣、主动脉瓣等。     图1血液循环系统模型  1.2 功率键合图模型 应用功率键合图建模方法的第一步是将原系统表达为功率键合图的图示模型。功率键合图由功率键、结点和作用元等主要元素构成，多分支血液循环系统的功率键合图如图2所示。   Rnv Chv Rhh Cha Rna图多分支血液循环系统功率键合图模型(此图有省略)   参考图2，绘制多分支血液循环系统功率键合图的步骤可简述如下： (1)根据对多分支循环系统各个功率流程分支的分

5、析，依次确定各0结点和1结点。 结点表示集总的流容容腔，如心室腔、主动脉弹性腔，在0结点处血液压力为等值，而该结点输入的血流量等于输出的血流量。结点表示集总的流阻管路或流感管路，如大动脉血管，在1结点处血流量为等值，而该结点的压力降等于上流压力值减去下流压力值。在图2 的循环系统模型中共有15个0结点和21个1结点。 (2)画上各结点周围的功率键，并标注功率流向。 功率键是带有箭头和因果线表示功率的线段。本模型中构成功率的两个变量是血压和血流。箭头表示系统作用元中的功率流向，即循环血液的流动方向。 (3)在功率键的一端标注上相应的C、R、L作用元。 为了能够全面、细致地刻画系统特性，本模型中应

6、用了三种作用元：流容、流阻和流感。 流容反映血管的顺应性，画在结点上，用C来表示，简称C元。例如，图2 中的Cta、Car、Cvn、Cpa、Cpv是分别表示与图1相对应部分的胸主动脉、大动脉、腔静脉、肺动脉和肺静脉顺应性的流容。 流感反映血流的惯性特性，画在1结点上，用L来表示，简称L元。如图2中的Lta、Lar、Lvn、Lpa、Lpv、Lco是分别表示相对应的胸主动脉、大动脉、腔静脉、肺动脉、肺静脉及冠状动脉血流惯性的流感。 流阻反映血流粘滞阻力的特性，简称R元，画在结点上。例如图2中Rta、Rar、Rvn、Rpa、Rpv和Rco是分别表示胸主动脉、大动脉、腔静脉、肺动脉、肺静脉及冠状动脉血

7、流粘滞阻力的阻性作用元。 (4)在各功率键上标注因果线，以便于建立系统的数学模型。 功率键上的因果线表示各作用元上流量与压力两变量之间的因果关系，确定了自变量和因变量，便于建立系统的状态方程。对于C元，其功率键上两个变量间，自变量是流量，因变量是压力；对于L元和R元，其功率键上两个变量间压力是自变量，流量是因变量。   经过以上步骤，就完成了循环系统的功率键合图模型。可以看出，键合图模型就是通过结点、功率键和作用元这些元素对心血管循环系统直观而形象的描述和反映。在将循环系统翻译成键合图模型后，就可以方便、有条不紊地推导系统数学模型。2 系统数学模型 功率键合图建模方法的第二步是推导系

8、统的数学模型。在推导系统动态过程的数学模型状态方程时，首先要确定状态变量。应用键合图方法建模的方便之处就在于对状态变量的确定有一定之规，可遵循固定的法则。 由于系统的状态方程是一阶微分方程组，在其变量间有导数关系，而在键合图中，只有流容C和流感L作用元中的两个变量间才有导数或积分关系，所以应当从C元和L元各自的变量间取一个变量作为状态变量。 对于C元，自变量为流量，因变量为压力，其关系为： (1) 对于L元，自变量为压力，因变量为流量，其关系为： (2) 对于R元，流量和压力之间的关系有： (3) 根据规则，取C元功率键上的压力变量p和L元功率键上的流量变量Q为状态变量,状态变量的一阶导数即为

9、状态方程。 因此，对于0结点，由(1)式两边取导数可得： (4) 其中， 是第i个0结点处的压力， 为输入血流量， 为输出血流量， 是第i个0结点处的流容。 对于1结点，由(2)式和(3)式可得： (5) 其中， 是第i个结点处的血流量， 为上流压力， 为下流压力， 和 分别是第i个结点处的流阻和流感。 对每个0节点和1结点都建立类似(4)和(5)的关系式，则可以得到系统的数学模型。本模型的数学模型是36阶的状态空间方程，即模型由36个一阶微分方程组成。下面列出了主动脉循环部分的状态方程： (6) (7) (8) (9) (10) (11)   其中，Cta、Caa、Car分别是胸主

10、动脉、腹主动脉、外周动脉的流容；Lta、Laa、Lar、Lvn分别是胸主动脉、腹主动脉、外周动脉和腔静脉的流感；Rta、Raa、Rsa、Rpc和Rsv是分别表示胸主动脉、腹主动脉、外周动脉、外周循环和腔静脉的流阻。ptao、paao、psar和Qtao、Qaao、Qsar分别是动脉循环中的胸主动脉、腹主动脉、外周动脉部分的压力和流量。 血液循环是由心脏的舒张收缩动作推动的，本文采用了心室时变流容 来表示这种舒张收缩动作， 是时间的周期函数。 对于循环系统中的膜瓣作用，可以作为模型的约束条件加入到系统数学模型当中：当血液正向流动时，膜瓣阻力为一较小的数值；当血液反向流动时，膜瓣阻力为无穷大，即阻止血液倒流。 本模型中的流容、流阻和流感参数参照文献4。      &#