生理系统建模与仿真2011-10.pptVIP

9 血液循环系统的建模与仿真 9.1 循环系统的生理知识 9.2 循环系统建模与仿真的基本方法 9.2.1 物理模型 9.2.2 数学模型 9.3 心脏模型 9.3.1 液压模型 9.3.2 泵模型 9.3.3 几何形体模型 9.4 血管网模型 9.4.1 动静脉血液动力学模型 （1）动脉血管模型 （2）静脉血管模型 9.4.2 循环系统整体模型 9.4.3 冠脉循环模型 思考题19 血压与血流的关系 Q为血流，l为血管长度，L为流体感应系数，R为流体阻尼 顺应性与跨壁压的关系 当存在塌陷时，实验数据还表明其血管顺应性将大幅提高。 血管内有瓣膜的等效电路模型 下肢的静脉血管为了防止血液倒流存在一些不完全闭合的单向膜瓣，在等效电路模型中常常采用切换流阻的方法模拟这些静脉血管中的血流状况。R2表示正向流阻，R1表示反向流阻。根据生理实验数据，一般取反向流阻为正向流阻的10倍。 左心室 在心脏的四个房室结构中，最重要的部分是左心室，这是由于心脏作为全身血液循环的泵的作用主要来自于左心室的周期性收缩和舒张活动。 心室收缩特性的描述 将心室视为用输入输出关系描述的液压装置 将心室视为一个泵，用心室的血压和体积之比来描述 把心室视为心肌纤维的集合，用心肌的机械特性，如张力-速度关系来描述 模型主要基于Frank和Starling原理 Frank实验 采用两个压力计记录心房和心室压力，发现当增加心室容积时，等容收缩压将增加，但当增加到一定限度时，此压力又呈下降趋势。 Starling实验 用水银计测定心室、主动脉、外周的压力，用可塌陷的管子代表外周阻力。 心脏定律（Frank和Starling原理） 在生理极限范围内，当静脉回流增加时，心脏容积逐渐扩大，同时心缩力增强，心输出量增加，直至每次收缩时的输出量与每次舒张时的回流量达到新的平衡为止；当增加心输出阻力，即升高动脉压时，起初每搏输出量减少，心中血液积聚量增多，随后，由于心舒期容积扩大，心缩力增加，能够克服增加的阻力，则每搏输出量恢复，重新取得平衡。这种输入与输出关系类似于一个液压装置，故称之为液压模型。 血压-容积关系（Suga等人的实验） 左心室和右心室的血压-容积关系可用一个时变函数E(t)来表示： 式中E(t)称为时变倒电容，它的生理意义为心肌的弹性系数。在相应的等效电路中，它与电容的倒数等价。plV为左心室血压，VlV为左心室容积，V0为在血压容积平面上E(t)曲线在收缩末期与溶剂轴的截距。 压力-容积曲线（为了求得V0） 时变弹性曲线 心室的压力变化 不同容积时，心室的压力变化： 心室压力变化的等效电路 心室的压力变化的等效电路可以用一个时变电容来表示 改进的心室等效电路等式及模型 同时考虑顺应性和阻抗特性 心室的几何形体模拟 心室可用一个空心的几何形体来模拟，常用的几何形状有： 球形 圆柱形 双焦点椭球形 心室的几何形体模型的子系统 心肌系统 心室几何形状系统 主动脉瓣系统 循环负荷系统 基于几何形体和机械特性的模型 T表示心室壁面张力，V表示心室容积 Avanzolini关系式 Avanzolini在1985年提出的左心室血压与心室几何形状和心肌张力之间的关系： pV为左心室血压，f(V)为一个无量纲的心室形状因子，可表示为心室形状的函数，因为容积是与形状有关的。 心室形状因子f(V) 球形 圆柱形 血管模型 p表示血压，Q表示血流，l表示血管长度，r表示血管半径，h为血管壁厚度。 Navier-Stokes方程 假设血液为不可压缩的牛顿液体，且血管横截面为圆形，则血液在血管中的流动过程可以用流体力学中的Navier-Stokes方程来描述： ?是血液的重力密度，v是血流速度，t是时间，?是血液的粘滞系数，g是重力加速度。 柱坐标系中的近似解 切向 轴向 径向 同心圆分割法（差分方法） 血管分割为N个同心圆柱，其中取中心圆柱的半径以及除最外层的其它各圆筒的厚度为?r，而最外层的厚度取为?r/2，并有 ?r = 2r0/(2N-1)，rn = (2n-1) r0/(2N-1) 血流速度梯形分布 实验表明，血液流动速度在血管外壁处为零，而在血管的中心最大，所以作为近似，假定血流速度沿半径方向为梯形分布，即在每一圆筒内为常数，在最外层筒内为零，然后由外至内逐级上升，在中心筒内为最大。 血压与血流的关系 Q为血流，l为血管长度，L为流体感应系数，R为流体阻尼，Y为杨氏模量，h为管壁厚度。 血管中血压和血流的关系类似于电路中电压与电流之间的关系，可以用一个等效电路来模拟血流在血管中的流动状态。 等效电路 电阻-等效流阻，电感-等效流感，电容-顺应性，电压-血压，电流-血流 血管塌陷过程 实验结果表明，薄壁管在内外压力差变化时将变形。若内外压差为零时的容