心血管系统的仿真与建模研究.docx

心血管系统的仿真与建模 前言 心血管循环系统是最复杂的生命系统之一，具有高度的动态特性。当心脏周期性地收缩和舒张时，心室射入主动脉的血流将以波的形式自主动脉根部出发沿动脉管系传播，这种波就是脉搏波。 人体心血管系统由心脏、动脉、毛细血管和静脉组成，构成体循环和肺循环两大回路，其中体循环在人体的血液循环系统中占据极其重要的地位。 体循环开始于左心室。血液从左心室搏出后，流经主动脉及其派生的若干分支动脉，进入各组织器官。动脉血经组织器官内的毛细血管完成氧气和营养物质的交换后变为静脉血，再由各级静脉汇集到上腔静脉和下腔静脉，回到左心房，从而完成了整个体循环过程。 描述心血管系统功能状态的主要生理参数有[1]：血压(收缩压、舒张压、脉压、平均动脉压)、心率、每搏输出量、血管顺应性、血流阻力、心肌收缩能力等，它们与人体状态有着密切的联系：生理参数反映了人体的生理病理状态，而当人体生理病理状态发生改变时，相应的生理参数也会随之改变。 本文的主要内容如下： 单弹性腔仿真模型 在Windkessel理论的基础上，建立心血管系统一阶Simulink仿真模型，结合临床数据计算模型参数进行仿真，分析实验结果与实测生理量的符合程度，从而验证模型对于心血管系统的表征能力。 双弹性腔仿真模型 在一阶模型的基础上，引入频率元件L，细化血管顺应性，建立三阶Simlink仿真模型。 参数性质分析 研究模型参数改变时脉搏波的变化趋势，分析各参数对脉搏波的作用。 2. 单弹性腔基本理论与仿真模型 2.1 Windkessel模型简介 对于图2.1所示的往复泵供水系统，柱塞P在马达M的驱动下做往复运动：当P向前挤压时，吸水阀门2关闭、供水阀门1开启，被挤压的流体经传输管路A，并通过终端阻力R最后流入贮水槽V；当P向后抽吸时，供水阀门1关闭、吸水阀门2开启，贮水槽中的水被吸入缸内，下一个周期又如此重复。 图2.1 水在往复泵及其传输管路中的流动[2] 实际应用中，人们往往在传输管路A中接入一个空气腔K以保证供水阀门关闭期间管路内液体流动的连续性。 流???在往复泵与其传输管路中的流动特性给血液在循环系统中的流动以启示。英格兰生理学家Stephen Hales(1677-1761)指出，由于主动脉的弹性扩张使心脏的周期性射血变成血管中血液的平稳流动。他把主动脉比拟为空气腔，并引入血液流动外周阻力的概念，认为外周阻力主要来自于人体组织中的小血管，由此发展为后来的Windkessel模型。 Windkessel模型在心血管系统与往复泵供水系统间建立起如下的比拟[3]： 心血管系统 往复泵供水系统 左心室 往复泵H 主动脉与大动脉 空气泵K 大动脉与毛细血管 终端阻力R 静脉与心房 注水槽V 主动脉瓣 供水阀门1 房室瓣（二尖瓣） 供水阀门2 当心室收缩、房室瓣关闭、主动脉瓣开启时，心室向主动脉射血：从心室射出的血液一部分经过动脉与毛细血管直接进入静脉腔，另一部分则贮留在主动脉与大动脉之中，使其产生弹性扩张(见图2.2a)；在心室舒张、房室瓣开启、主动脉瓣关闭期间，心脏停止向主动脉射血，但在主动脉与大动脉管壁弹性恢复力的作用下，贮留在其中的血液继续向前流动(见图2.2b)。由此可见，主动脉与大动脉管壁的弹性作用使心脏间断的射血转变为动脉管道中连续的血流。 a 收缩期 b 舒张期 图2.2 主动脉与大动脉的弹性腔作用 Windkessel模型将主动脉与大动脉比拟为一个弹性腔，对于某一确定时刻，认为腔内的各种流动参数(如血液压力P和流量g等)仅是时间，的函数，而与距离心脏的距离x无关。Windkessel模型只是分析循环系统中血液流动的一种相当粗糙的模型，但是它能反映动脉血流的一些重要特性，至今仍广为采用。 2.2 单弹性腔模型的原理 图2.3 动脉弹性腔的血流量图 如图2．3所示的单弹性腔模型[4]，qin表示单位时间内从心脏流入动脉弹性腔的血液体积，qout表示单位时间内从动脉弹性腔经微动脉与毛细血管流入静脉腔的血液体积，V为这段动脉血管的容积。 生理压力范围内，忽略血液体积因压力而严生的变化，并认为血液在血管中的流动是连续的，得到心室收缩期和舒张期内血流流动的连续性方程[4] qout+dVdt=qin (2.1) qout+dVdt=0 (2.2) 在弹性腔模型假设下，动脉管中的压力处处相等记为p；静脉管中的压力也处处相等，记为pv，令△p=p-pv表示血液由动脉管留至静脉产生压力降，R为外周阻力，由△pq=R有 p-pvR=qout (2.3) C=p-pvR (2.