心血管系统的仿真与建模Word文档下载推荐.docx
- 文档编号:18899048
- 上传时间:2023-01-02
- 格式:DOCX
- 页数:15
- 大小:347.88KB
心血管系统的仿真与建模Word文档下载推荐.docx
《心血管系统的仿真与建模Word文档下载推荐.docx》由会员分享,可在线阅读,更多相关《心血管系统的仿真与建模Word文档下载推荐.docx(15页珍藏版)》请在冰豆网上搜索。
左心室
往复泵供水系统
往复泵H
主动脉与人动脉
空气泵K
大动脉与毛细血管
终端阻力R
静脉与心房
注水槽V
主动脉瓣
供水阀门2
房室瓣(二尖瓣)
当心室收缩、房室瓣关闭、主动脉瓣开启时,心室向主动脉射血:
从心室射出的血液一部分经过动脉与毛细血管直接进入静脉腔,另一部分则贮留在主动脉与大动脉之中,使其产生弹性扩张(见图2.2a);
在心室舒张、房室瓣开启、主动脉瓣关闭期间,心脏停止向主动脉射血,但在主动脉与大动脉管壁弹性恢复力的作用下,贮留在其中的血液继续向前流动(见图2.2b)o由此可见,主动脉与人动脉管壁的弹性作用使心脏间断的射血转变为动脉管道中连续的血流。
图2.2主动脉与人动脉的弼性腔作用
Wiiidkessel模型将主动脉与人动脉比拟为一个弹性腔,对于某一确定时刻,认为腔内的各种流动参数(如血液压力P和流量g等)仅是时间,的函数,而与距离心脏的距离x无关。
Wiiidkessel模型只是分析循环系统中血液流动的一种相当粗糙的模型,但是它能反映动脉血流的一些重要特性,至今仍广为采用。
2.2单弹性腔模型的原理
图23动脉弹性腔的血流量图
如图2・3所示的单弹性腔模型⑷,qm表示单位时间内从心脏流入动脉弹性腔的血液体积,爲讥表示单位时间内从动脉弹性腔经微动脉与毛细血管流入静脉腔的血液体积,V为这段动脉血管的容积。
生理压力范围内,忽略血液体积因压力而严生的变化,并认为血液在血管中的流动是连续的,得到心室收缩期和舒张期内血流流动的连续性方稈⑷
dVs八
Qout+乔—Qin(il)
Qout+鴛=°
(2・2)
在弹性腔模型假设下,动脉管中的压力处处相等记为P;
静脉管中的压力也处处相等,记为令△p=p-p“表示血液由动脉管留至
静脉产生压力降,R为外周阻力,由汁R有
(23)
(2.4)
其中C为血管顺应性。
根据(2.1)・(2.4),可得
畔+屮二伽少之丸)(2.5)
C务+护=0(7;
StST)(2.6)
其中,Ts表示心室收缩期的持续时间,T表示心动周期。
(2.5)、(2.6)便是弹性腔模型基础上,动脉压力p在收缩期与舒张期中分别必须满足的两个方程,各参数取值的情况如下:
(1)表示静脉腔屮的血液压力,通常取恒定值=9.3X102-
17.3x102Pa(7-13nmiHg),方便起见,有时也取为零;
(2)R表示外周阻力,看座与压力p无关的常量;
(3)仏表示心脏收缩射出的血液流量,随时间t而变;
(4)C表示动脉顺应性,可以使与血压无关的常量(线性弹性腔理
论),也可以是一个随血压而变化的变量(非线性弹性腔理论)。
下面我们将基于非线性弹性腔理论建立心血管系统的Siimihiik仿真模型。
2.3基于单弹性腔的Simulink仿真模型
2.3.1人体心脏搏动间歇模型
在仿真试验屮我们用正弦曲线仿真心脏的搏动间歇流qin(t):
心脏
收缩时,心室射出的血液流量按照正弦规律变化:
而当心室舒张时,
血流量为零。
这样能近似反映出心脏间断的射血过程⑴
•nt
Gn(r)=
(0<
t<
7;
)
0(7;
<
T)
(2.9)
Sinniliiik仿真模型如图2.4所示。
2.3.2人体动脉管的p-V关系
对于一阶模型而言,我们近似认为动脉p、v满足如下的指数函数关系⑹
V=aebp+c(2.7)
其中:
a,b,c均为常数,需要通过拟合法求值。
此时动脉顺应性C将表示为
C=J=abebp(2.8)
dp
对于各类动脉,常数3的取值变化范围较人,b值则是比较接近
的。
我们通过查阅资料,得到b的平均值为b=-0.0131o
Sinniliiik的关系图为图2.5所示。
图24心脏搏动间歇模型
eu
<
z>
P
2.3.3系统整体框图
根据以上的讨论,我们得出心血管系统一阶Shmihiik的仿真模型
如图2.6所示。
图26心血管系统一阶仿真模型
2.3.4参数分析
我们通过查阅资料⑺,取b=-0.0131,Ts=0.15s,T=0.8s,其他参数取典型值。
仿真时间设置为20s,系统达到稳定状态后的基本波形如图2.7所不。
图27基本仿真波形图
对于图2.7的基本仿真波形图。
上排曲线表示心脏每搏输出波形
(单位:
inl/beats),下排曲线表示心血管系统一腔室模型内的血压变化(单位:
limiHg),横坐标为时间轴(单位:
s)o经过若干个心脏循环周期,系统逐渐达到稳定的工作状态。
对于不同的人体样本,模型中各参数的取值略有不同,而相应的波形也会有变化。
2.4本章小结
基于Windkessel理论建立起的心血管系统一阶仿真模型,较好地模拟了动脉血压变化,缺点是无法描述脉搏波的波形细节,下一章我们将具体介绍一种双弹性腔模型的建立方法。
3.双弹性腔仿真模型的研究
由于单弹性腔模型过于简化,无法解释舒张期所出现的一些诸如潮波、重搏波Z类的脉搏波波纹,而正是这些波纹在临床上有重要的童义。
为了提高分析的精度,根据对脉搏压力曲线和流量曲线进行傅里叶分析后,人们认为应该在模型中引入一些与频率有关的元件,必须考虑血液流动的黏滞惯性。
在前人的研究基础上,1967年由美国的戈特温和瓦特共同提出了双弹性舲模型⑻,该模型采用两个串联的弹性腔,以表现血管系统的不同压力,同时还在两个弹性腔体之间加入一个表示血液惯性的环节,使模型输出的脉搏波曲线能很好地反映岀舒张期波纹。
双弹性腔模型将主动脉及其主要分支视为两个弹性腔,顺应性分別为C]和C2。
心室收缩时,从心室泵出的血液进入第一个弹性腔(用来表征主动脉弓及其主要分支的集总顺应性效应)与血液惯性元件三(用来表征主动脉中的集总血液惯性效应),而后进入第二个弹性腔(用来表征腹主动脉及其主要分支的集总顺应性效应),最后流经集总的外周阻力R(用来表征外周血管床的总阻力)进入静脉腔⑼。
下面我们将利用这种双弹性腔理论,建立对应的心血管系统三阶模型。
3.1双弹性腔模型的基本方程
心皇弹性腔】弾性腔n
图31循环系统双弹性腔的血流量图
如图3.1所示的双弹性腔模型,假设弹性腔I的顺应性为G,舲内血液压力为內,弹性腔体积为弹性腔II的顺应性为C2,腔内血液压力为P2,弹性腔体积为冬,两个弹性腔Z间连接着一个长为1,截面积为A的血液柱,并用一个集总的流动惯性L来体现血液在主动脉中的惯性效应。
假设心室将血液流量q讯函射进弹性腔I和惯性兀件L,冉进入弹
性腔II,最后流经外周阻力R进入静脉腔。
动脉屮血液的流动是脉动的,其压力和流量都是一些随时间作周期性变化的量,可以看成是由“直流”成分(即平均压力和平均流量)和“交变”成分组成,这样便町以用构成电路相类似的方法引进血流阻力元件、顺应性元件与惯性元件,用这些元件建立起相应的血液流动等效回路〔⑼。
对于血流阻力元件R与惯性元件L,只要将它们依次画在(即串联在)相应的两个压力点內、卩2之间即可。
对于反映血液积聚能力的血管顺应性元件C,虽然它是与电容相比拟的,但是却不能像电容一样可以随意串联在回路的任意地方。
在建立血液流动的等效回路时,顺应性元件必须与其后继元件并联,而不能串联。
这是因为顺应性代表弹性腔室,其压力差为腔室压力与基准压力之差,而基准压力一般都选用人体周围外界的大气压力(可看作零压力),所以顺应性C两端总是一端连接腔室压力另一端连基准压力(零压力)。
如果不考虑连接两个弹性腔的血液柱的粘性阻力,
取静脉作为等效回路中的地(零压力),不难看到对应的等效回路为
图32双弹件腔模型的等效电路
由这个等效回路,我们可以建立起相应的方程【⑴。
根据血流量守恒的
条件,知
同时,根据定义可知
(2.10)
式中:
Pv——静脉腔内血液压力;
qi——通过血液惯性原件1的血流量;
尙和绻2一一贮存在弹性腔I和II中的血液流量;
由(2.7)-(2.12),不难得到描述双弹性舲模型的基木方程
在这组基本方程中,静脉压內和每搏心输出量q汛通常认为是已知的,因而通过这三个方程,可以确定出待求的三个未知量7,山与乃。
3.2Simunnk仿真模型的建立
321双弹性腔仿真模型的构建
对于基本方程中的所有除法运算都用增益表示,所有的微分都用逆向积分运算表示,心输出量仍釆用一阶仿真中的搏动间歇流模型,并假设从心室搏出的血流中96%流入主动脉网,据此可做出心血
...q0sin—(0<
Ts),.s._,,m丄十
其中%:
=Ts'
"
的计算过程如下:
10(Ts<
t<
T)
假设心脏每搏输出量70mL(即一次心搏屮由一侧心室输出的血液量为70mL),有
(TsTit
gosin—dt=70JoLs
fTsntr03ntnt
I%sin耳dt=q°
Jsin-d(-)
70ti
"
芮"
6・519
3.2.2基本仿真波形
动态仿真20秒,记录基本波形如图3.4所示
ml/beat•s
图34双弹性腔模型基木仿真波形
比较以上两图可以看出,由双弹性腔模型仿真得到的脉搏波波形特征与实际测量的波形更为相似,并且相对于一阶模型而言,双弹性腔模型较好地反映了波形细节,尤其是脉搏波卜•降过程屮的重搏波现象。
3.3参数选择
孙冬梅阳等人的实验结果表明,运动过程中脉搏波的幅值、频率和波形特征都会发生相应的变化。
对于25名受试着,运动负荷用自行车功量计模拟,负荷从50w开始,lmin为一个时间间隔,负荷增量为25W,逐级增加到预定负荷,并计算出对应的R、L、Ci和C2的变化趋势。
R反映血管外周阻力的人小。
运动过程中由于人量血管扩张和运动肌肉分支血管增加,外周阻力总体呈减小趋势,是心血管系统对运动刺激的一种适应性反应"
J
R的变化范围:
0.40—1.13
L反映血液流动惯性的大小。
运动中由于大量血管的扩张引起流动惯性的减小,使血液流动变得容易,有利于运动过程中对血流量增加的需要,L呈下降趋势。
/•的变化范围:
0.009—0.019
C]反映主动脉弓及其主要分支的集总顺应性。
运动中C丄的增加表明心脏射血时因弹性扩张而增加对血液的积聚作用,从而缓冲动脉压力运动中的过度升高;
但随着运动负荷量增加,心室的有效泵力增强,动脉血压升高,血流速度加快,主动脉及人动脉对血液的积聚能力降低大量血液直接流向外周血管,从而导致C]的下降。
C1的变化范围:
0.162一1.402
C2反映腹主动脉及外周血管总的顺应性。
由于运动负荷作用,肌肉细胞的新陈代谢加快,使得人量的外周血管舒张,血流通畅,同时供给肌肉血流量的血管分支增多,因而C2呈升高趋势。
的变化范围:
0.082—0.237
对于3.2节中建立的双弹性腔模型,我们选择以下模型参数
Cx=l.l
C2=0.15
/?
二0.85
L=0.014
K=0.96表示主动脉的分流系数
将上述参数代入图3.3所示的双弹性舲仿真模型,输出端波形如图3.4所示,脉搏波幅度范用60一94mmHgo考虑到血液经过主动脉与腹主动脉的阻力衰减作用到达末端动脉时,压力有所降低;
另一方面外周动脉的顺应性比主动脉低,脉搏压较人,即收缩压较高,而舒张压较低。
因此以上仿真波形基本符合生理实际。
3.4本章小结
建立出人体心血管系统双弹性腔模型(即三阶模型),并讨论了参数取值情况。
模型较准确地描述了脉搏波特征,尤其是下降支的重搏波,仿真波形与临床检测的正常脉搏波基本一致。
4模型参数与波形特征的关系
脉搏波由心脏舒缩运动产生,并经主动脉根部传播到全身动脉系统中。
挠动脉脉搏波与心血管系统中许多生理参数有关。
由于脉搏波波形复杂,以及在体实验不可取,直接分析其波形有一定困难,本章我们利用考察心血管参数对脉搏波特征值的影响来代替对脉搏波的分析。
只考虑单一参数对脉搏波的影响,并将仿真结果与正常波形进行分析比较。
在考察参数值变化对脉搏波的影响时,只考察波形特征的变化趋势,即收缩压是随参数值的增加而增加,还是随参数值的增加而减小。
可以证明,这些变化是有理论根据的。
4.1双弹性腔模型等效回路的拉普拉斯分析
对于图3.2所示双弹性腔模型的等效回路,在复频域进行拉普拉斯变换刖
恥)普+警(4.1)
/SC2K
q,(s)=^g^(4.2)
qm=0对应于心室舒张期,与脉搏波的降支部分相联系,包含着心
血管系统丰富的生理信息,如主波、重搏波等振荡波纹,极具有研究价值。
心舒期内心脏停止射血,此时的脉搏波相当于心血管系统的零输入响应,特征方程为
S34-—S24-f—+—)S4-—=0(4.4)
RC2\LCiLC2JRLC\C2
方程(4.4)解的一般形式为
p2(0=心+a2e~a3+a4e~a5tcos(a6t+a7)(4.5)其中:
a^a7为系数。
它是由直流分量指数衰减分量a2e~a3f以指数函数为包络线的衰减分量a4e_Qr5fcos(a6t+a?
)三部分叠加而成。
当。
1~。
7取不同的组合,即C】、C2、L.R取不同值时,可形成具有各种振荡波纹的曲线,较好地描绘出脉搏波的主要特征点如主波、重搏波等,反映不同参数对脉搏波形特征的作用,下面分别介绍G和
4.25对脉搏波波形的影响
图41q=0.162.LI.1.402时的仿真波形
如3.3节所述,C]反映主动脉弓及其主要分支的集总顺应性,C]增加表明主动脉弓及其主要分支的集总顺应性增大,而
dV
c=d^
表示当压强增大一个单位时体积的增犬量增加。
对于动脉管段而言,随着C丄的增加,血管对于血液快速流动对管壁产生冲击力的适应性反应-…管壁收缩或者舒张也越容易,即心室射血期内管壁的扩张效应更加明显,而心室舒张时管壁的紧缩效应也更明显。
这就是说,当Cl增加时,收缩压变小,舒张压变大,脉搏压也随Z减小。
由以上分析可以看出,q主要影响脉搏波幅度。
4.3C2对脉搏波波形的影响
mmHg
0-0OOH.
ft)
10、
知
TO
a】fIQ
IVJ
w
■•
曲
■Mt.
4・
图4.2G=0.082、0.15、0.237时的仿真波形
由3.3节知,$反映腹主动脉及外周血管总的顺应性。
由于分析对象为脉搏波,可将视为上肢动脉及外周血管总的顺应性。
当血液从主动脉流向上肢动脉时,管壁厚度、血管半径、血压等均逐渐减小,血流速度减慢,此时尽管上肢动脉及外周血管总的顺应性增人,但相对于主动脉而言,其功效几乎可以忽略,故收缩压和舒张压的变化范围很小。
理论上来说,当C2逐渐增大时,收缩压应有所减小或大致保持不变,但仿真结果显示收缩压缓慢增加,出现这种差异可能有两个原因:
①仿真中用的是集总参数的电网络模型,也许不能正确模拟动脉管中血液流动情况,而得到与以往研究不同的结果;
②临床上促使Q增大的因素也会对血管造成其他影响,这些影响相互叠加的结果是不定的。
由以上分析可以看出,C2主要影响重搏波波幅。
4.6本章小结
本章讨论了模型参数与波形Z间的关系,并且对双弹性腔模型的等效回路进行了拉普拉斯分析。
在参数分析中主要讨论了C]和C2的不同取值对波形的影响。
C]反映主动脉弓及其主要分支的集总顺应性,它主要影响脉搏波幅度。
C2反映腹主动脉及外周血管总的顺应性,它主要影响重搏波波幅。
5参考文献
1柳兆荣.弹性腔理论及其在心血管系统分析中的应用.北京:
科学出版社,1987
2Frank,ODieGnandfoiTndesartenellenPulsesZfBio:
1899,37483-526
3孙磊.心血管系统仿貞-建模与脉搏波分析研究[D].浙江大学,2008
4宁钢民,代开勇,李英奇等心血管系统键合图模型研究[J].浙江大学学报:
工学版,2007.(5):
864-870
5郝卫亚,昊兴裕,张立藩等.循环系统的数学模型及仿貞•实验[J].医学争鸣,2000,⑴.
6陈丽琳,昊效明,杨艳等.心血管循环系统的建模仿戌[』.北京生物医学工程,2006,(3).
7.宁钢民,代开勇,李英奇等.心血管系统键合图模型研究[J].浙江大学学报:
工学版,2007,(5):
864-870.
8孙磊.心血管系统仿貞•建模与脉搏波分析研尢[D].浙江大学,2008.
9CopeFW.AnElasticreservoirtheoryofthehumansystemicarterialsystemusingcurrentdataonaorticelasticity』.Thebulletinofmathematicalbiophysics,1960,22:
19-40.D0l:
10.1007/BF02477969.
10•周希平,王国卿,邹容.基于Simulink仿克的心脏工作模型[J]中国血液流变学杂志,2004,14(3).DOI:
10.3969/j.issn.1009-881X.2004.03.009.
11.徐礼胜,刘佳,何殿宁等.基于双弹性腔模型与整心动周期脉搏波的血流动力孑参数估计[JJ.北京工业大学学报ISTICElPKU,2014,40(10).
12刘静纨,陈建慧,罗志昌.人体微循环弹性腔模型与模型参数佔计[J].北京生物医学工程,2001,
(2).
13載娟.左心循环系统建模及研Jl[D].哈尔滨工程大2008.DOM0.7666/d.y1437855.
14白净,不详.血液循环系统的数字仿真[J].中国科学基^,1994,
(2):
112-116.
15孙冬梅,罗志昌.人体心血管双弹性腔模型参数佔计及其在运动医学中的应用[J].北京匸业大学学报,1996,
(1):
62-70.
16MAS,刁颖敏.心脏•肺循环•体循环系统建模初探[J].同济大学学报:
自然科学版,2002,30
(1).DOI:
10.3321/j.issn:
0253-374X.2002.01.013.
17徐林,昊效明,陈册琳.心血管循环系统建模[J].中国组织工程研究与临床康复,2009,13(20).DOI:
10.3969/j.issn.1673-8225.2009.20.039.
18郭珞源,苑康上孚刚.人体血液循环动力学系统设计仿貞•研尢[』•中国医学影像技术,2003.
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 心血管 系统 仿真 建模
![提示](https://static.bdocx.com/images/bang_tan.gif)