数值计算分两步进行。首先给定各点压力初值，计算各单元的初始阻抗、顺应性。然后计算每一时步的阻抗、顺应性和节点压力，直到整个心动周期结束。最终结果的收敛准则按生理平衡条件考虑：在一定误差范围内(ε=0.02)，一个心动周期内动脉输入的血流容积等于静脉输出的血流容积。
3 仿真结果
对正常心肌所作的仿真实验，表明其冠状动脉、静脉血流和三层心肌中各血管段的血流分布模式符合现有的实验观察和定性分析结论。为了简化心肌灌注的力学因素，仅考虑缺血部分为全部梗塞、梗塞以下的血管处于最大舒张状态。
缺血状态下，心肌时相血流分布如图2所示(篇幅所限，只给出侧支循环发展程度为50%的情况)。心肌各层静脉血流基本上为正向血流(从主动脉到右心房方向为正)，在收缩期出现峰值，其幅值随着侧支循环发展程度的增加而增大。心外膜动脉血流在收缩期出现正向峰值。心内膜动脉血流在收缩期开始段出现反向峰值，然后由于血管的塌陷而出现零流现象。舒张期塌陷血管复原后，又有逆流现象出现，其幅值小于压塌前逆流的幅值。侧支循环发展程度越低，逆流峰值越小，零流持续时间越长。各层血管网的阻抗都会受侧支循环发展程度的影响，其中心内膜血管阻抗变动最大。侧支循环发展程度越低，心内膜血管阻抗值越大，说明心肌内膜最易缺血。
图2 缺血心肌灌注血流分布
Fig 2 The blood flow distribution of the different myocardial and vessel segments in ischemic myocardium
a.large and middle arteries; b.small arteries; c.small veins; d.large and middle veins
4 提高灌注压对血流时相分布的影响
