基于非牛顿流体的体外反搏下冠状动脉血流动力学研究 (1).pdf
《基于非牛顿流体的体外反搏下冠状动脉血流动力学研究 (1).pdf》由会员分享,可在线阅读,更多相关《基于非牛顿流体的体外反搏下冠状动脉血流动力学研究 (1).pdf(5页珍藏版)》请在咨信网上搜索。
1、基于非牛顿流体的体外反搏下冠状动脉血流动力学研究赵 艺1,陈 熹2 1.哈尔滨工业大学(深圳)教务部,广东深圳 518055;2.首都医科大学生物医学工程学院,北京100069摘 要 目的:基于心血管患者的冠状动脉血管模型与血液的非牛顿黏性特征,通过计算和分析血管内血液流动的各个动力学参数,探讨增强型体外反搏对冠状动脉狭窄处血流动力学的影响。方法:基于心血管患者冠状动脉造影的 CT 图像,用 Mimics 软件对左冠状动脉及其主要分支进行三维重建;利用有限元方法,对冠状动脉模型进行流体力学计算。假设血管壁不发生变形,血液为不可压缩流体,血液黏性符合非牛顿流体规律,分别计算在体外反搏前和处于体外
2、反搏时,在非定常速度进口的条件下的血流动力学参数的变化规律。结果:相对于体外反搏前,在体外反搏时,流线中的回流和二次流基本消失,而且在狭窄处下游,壁面剪切应力明显变大,振荡剪切指数明显变小。结论:体外反搏特征对冠脉血流动力学特征有重大改善,可以起到辅助治疗心血管疾病的作用。关键词 动脉粥样硬化;冠状动脉狭窄;非牛顿流体;体外反搏;血流动力学基金项目:国家自然科学基金青年项目项目编号:11802187;广东省基础与应用基础研究基金项目项目编号:2019A1515111178通讯作者:陈 熹The improvement of hemodynamics index of coronary arte
3、ries during Enhanced External Counterpulsation in patients with coronary heart disease based on non-Newtonian fluid ZHAO Yi1,CHEN Xi2 1.Department of ac-ademic affairs,Harbin Institute of Technology(Shenzhen),Shenzhen 518055,China;2.School of Biomedical Engineering,Capital Medical University,Beijing
4、 100069,ChinaAbstract:Objective Based on the coronary artery model of cardiovascular patients and the non-Newtonian fluid of blood,the effects of Enhanced External Counterpulsation(EECP)on the hemodynamics of coronary artery stenosis were studied by calculating and analyzing various dynamic paramete
5、rs of blood flow in the blood vessels.Method Based on CT images of coronary angiography in car-diovascular patients,3D reconstruction of the left coronary artery and its main branches was performed using Mimics software;Fluid mechanics calculation of coronary artery model was carried out by using fi
6、nite element method.Assuming that the vessel wall does not deform,the blood is incompressible fluid,and the blood viscosity conforms to the law of non-Newtonian fluid,calculate the change of hemodynamic parameters under the condition of unsteady velocity inlet before and during EECP,respectively.Res
7、ults Compared with that before EECP,during EECP,the vortex and secondary flow in the streamline disappeared,and the time-averaged wall shear stress(TAWSS)significantly increased and the oscillatory shear index(OSI)significantly decreased downstream of the stenosis.ConclusionEECP have a significant i
8、mprovement in coronary hemodynamic characteristics,these findings have important implications on EECP as adjuvant therapy before or after percutaneous coronary intervention(PCI)in patients with diffuse atherosclerosis.Key Words:Atherosclerosis;Coronary artery stenosis;Non-Newtonian fluid;Enhanced Ex
9、ternal Counterpulsation;Hemody-namics 弥漫性冠状动脉粥样硬化是心血管疾病中最常见的一种,经常伴随其他病症,如脂蛋白等,严重危害着人类的生命健康1-2。冠状动脉狭窄(尤其是面积狭窄在 75%以上的狭窄)会导致心肌供血不足,心肌湍流杂音,进而导致动脉粥样硬化等各种心血管疾病,对临床诊断和治疗有重要意义3-4。血液黏度也可能会导致心血管疾病死亡风险增加,在临床检验、病兆诊断上都具有一定的意义5-6。而在心血管血流动力学模拟的研究中,前人往往都忽略了血液的非牛顿黏性特征。增强型体外反搏(EECP)在心血管疾病的辅助治疗方面的应用越来越广泛7-9,可以辅助血小板的聚
10、集10,辅助降低顽固性心绞痛患者的血压11,尤其是对老年人的心血管疾病辅助治疗有重大帮助12,甚至对心力衰竭患者有一定的辅助治疗作用13。然而,目前对于体外反搏的血流动力学研究往往将血液简化为牛顿流体,忽略了血液的黏性作用14-16。本文基于心血管患者的 CT 造影,基于血液非牛顿流体的黏性特征,利用三维图像重建和有限元分析的方法,探讨体外反搏对冠状动脉狭窄条件下各流体力学参数的影响,说明体外反搏有助于心血管疾病的精确诊断和辅助治疗。8633吉林医学 2023 年 12 月第 44 卷第 12 期1 资料与方法1.1 资料:本文所用 CT 胸部造影图片来自北京市北京大学第三医院某诊断为动脉粥样
11、硬化的患者,是在获得了患者的知情同意下使用的。1.2 方法:将 CT 图片导入到 Mimics 软件(Materi-alise,10.01)中,对左冠状动脉 包含左前降支(LAD)和左回旋支(LCx)及其主要分支进行三维重建。从 Mimics 软件中,用 STL 格式文件输出,导出到 ICEM 软件中(ANSYS,14.5)进行有限元网格划分(非结构化网格)。将网格导入 Fluent 软件(AN-SYS,14.5)计算,假设血液为不可压缩流体,密度为1 060 kg/m3,牛 顿 流 体 假 设 下 的 黏 度 为0.004 5 kg/ms。在非牛顿 Carreau 流体假设下,流体的表观黏度
12、可写为17:=+(0-)1+()(n-1)/2 (1)其中 为血液的表观黏度,和 0分别为血液在剪切率为无穷大和 0 时的黏度,=0.25 kg/(ms),0=0.002 5 kg/(ms),为剪切率,为黏度随剪切率变化的时间常数(=25s),n=0.25为非牛顿流体的幂律指数。血液的雷诺数为 Re=vD/,其中 为流体密度,v 为流体的速度,D 为圆管内径,为流体黏度或表观黏度。基于本研究的假设,血液的雷诺数约为 600,因此血液流动可假设为层流。血液流动计算基于三维、非定常、不可压缩纳维斯托克斯方程:vt+vv=-p+v (2)以及连续性方程:v=0 (3)其中,v 为速度矢量,t 为时间
13、,p 为压强。为了获得真实的肱动脉压力波,我们用软件进行了脉搏波分析(SphygmoCor Version 9,AtCor Medi-cal Pty.Ltd,Australia)。在 EECP 治疗之前和期间对患者进行了测量,肱动脉搏动产生的压力振荡会传递到肱动脉血压袖带,由传感器测量,然后输入微处理器,计算机软件记录肱动脉的脉搏波,如图 1(A)和图 1(E)所示。入口压力波边界条件采用了获得的真实压力波,其中图 1(A)是 EECP 前患者的主动脉压力波,图 1(E)是 EECP 时患者的主动脉压力波。用 Fluent 软件中的 User defined function 功能编写和读取入
14、口压力波形。出口使用阻力边界条件,与文献18一致。血管壁设为无滑移刚性壁面。计算时间为 4 个心动周期,一共 2.52 s,取最后一个心动周期的结果作为我们的计算结果。将随时间变化 的 结 果 导 入 tecplot 软 件 中(Tecplot 360 2011),用编写的 Equation 计算出各个位置的平均壁 面 剪 切 力(TAWSS)19和 振 荡 剪 切 指 数(OSI)20,并画出流线。其中,OSI 是衡量壁面剪切力脉动特性的重要参数,其表达式为:OSI=121-|T0wdt|T0|w|dt (4)其中,w为壁面剪切力,T 为总时间。同时,本研究也计算了壁面剪切应力危险区域比率(
15、SAR-TAW-SS)和振荡剪切指数危险区域比率(SAR-OSI),分别是指壁面剪切应力危险的区域面积(TAWSS1.5)占整个血管壁面的百分比。2 结果2.1 EECP 前的冠脉血流动力学分析:患者的左冠状动脉在 LAD 支的主干第一级下游和第二级上游分别有一处狭窄,根据测量,其面积狭窄分别为75%和 65%,下面本文研究这两处狭窄对血流动力学参数的影响。图1(B)为患者左冠状动脉的平均壁面剪切力,可以看出,在 LAD 主干的两处狭窄处,TAWSS 有一定的升高,而在狭窄下游,TAWSS 要明显降低,这种剪切力的急剧变化,有可能会形成血管斑块,严重时会造成斑块破裂。而在 LCx 的第1 级分
16、支处有轻微狭窄,下游处的 TAWSS 相对较高,说明 LCx 分支处的病变并不是很严重,流线也相对较平稳,不容易产生涡流和二次流。相对应的,图1(C)为患者左冠状动脉的振荡剪切指数,可以看出,在 LAD 主干的两处狭窄处,OSI 都有一定的升高区域,说明在该区域中,血管有明显的病变发生,壁面剪切应力的变化特别剧烈,较容易发生冠状动脉粥样硬化。而在 LCx分支狭窄的地方,OSI 都有一定的高数值区域,同样对应着病变的发生,虽然是轻微狭窄,但是也应该引起足够的重视。图 1(D)为患者 LAD 主干第 2 级分支一处狭窄附近的流线。可以看出,在主干第二个狭窄下游,出现了回流和涡流,在分支狭窄的下游,
17、出现了二次流。在血管中出现的回流、涡流、二次流往往会导致局部流动紊乱,并对动脉壁产生极大的影响,可能会导致心肌供血不足以及动脉粥样硬化的发生。9633吉林医学 2023 年 12 月第 44 卷第 12 期 注:图1 A:患者 EECP 前的肱动脉压力波波形;B:EECP 前的平均壁面剪切力(单位:dyn/cm2);C:EECP 前的震荡剪切指数(无量纲);D:EECP 前的 LAD 主干 2 级分支处的流线;E:患者 EECP 时的肱动脉压力波波形;F:EECP 时的平均壁面剪切力(单位:dyn/cm2);G:EECP 时的震荡剪切指数(无量纲);H:EECP 时的 LAD 主干 2 级分支
18、处的流线 图 1 EECP 前冠脉血流动力学分析2.2EECP 时的冠脉血流动力学分析:当患者做EECP 时,血管内的 TAWSS、OSI 和流线分别如图 1(F)、图 1(G)、图 1(H)所示。此时患者的肱动脉压力波在舒张压时有一个非常明显的高压出现,如图1(E)所示,此高压提升了主动脉的入口压力,有效增加了患者主动脉的心输出量,对血流动力学参数有一个较为明显的提升。对比图 1(F)和图 1(B),本研究发现在 EECP 时,LAD 主干的两处狭窄处的低剪切力区域虽然仍然存在,但范围缩小,而且剪切力值变化梯度不是特别大,这样就降低了剪切力在局部区域急剧变化的可能性,患者动脉壁的损害在血液黏
19、性自主变化的条件下得到了一定的控制。对比图 1(G)和图 1(C),LAD 主干处的高 OSI 分布区域虽然依然存在,但数值明显降低,说明患者动脉壁上受到的剪切力脉动刺激已经降到了比较低的值,壁面剪切应力的变化,在 EECP 的影响下已经得到了明显控制。而图 1(H)中的主干处流线也变得相对规律,回流和二次流基本消失了,涡流也相对比较平稳,不会再引起心肌供血不足。这就是 EECP 时患者冠脉的血流动力学参数的变化,从各个角度来看,肱动脉的压力波型的变化,提升了冠脉血流动力学的壁面剪切应力,减小了壁面剪切应力的剧烈变化,使得冠脉狭窄上下游处的流线变得平稳,起到了一定的辅助治疗作用。3 讨论本文结
20、果表明,在非牛顿流体黏性血液假设下,冠状动脉主干处的连续狭窄对血液动力学参数的影响十分明显。狭窄会导致壁面剪切力的升高以及局部剪切力急剧变化,导致局部振荡剪切指数的急剧升高,进而导致狭窄下游处的回流、涡流以及二次流的产生。血管狭窄处会自然地产生壁面剪切力升高的现象,而血液的黏弹性记忆效应,会在剪切力急剧变化时,通过改变自身的黏度,对这种变化起到一定的阻碍作用。图 1(D)中 LAD 主干下游处出现的涡流和二次流,是因为分支处血管尺度过小,主干处黏度的变化过大,当血液流到分支狭窄处时,就会产生一定的紊乱。在血液黏性的影响下,患者 LAD 主干0733吉林医学 2023 年 12 月第 44 卷第
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 基于非牛顿流体的体外反搏下冠状动脉血流动力学研究 1 基于 牛顿 流体 体外 反搏下 冠状动脉 血流 动力学 研究
1、咨信平台为文档C2C交易模式,即用户上传的文档直接被用户下载,收益归上传人(含作者)所有;本站仅是提供信息存储空间和展示预览,仅对用户上传内容的表现方式做保护处理,对上载内容不做任何修改或编辑。所展示的作品文档包括内容和图片全部来源于网络用户和作者上传投稿,我们不确定上传用户享有完全著作权,根据《信息网络传播权保护条例》,如果侵犯了您的版权、权益或隐私,请联系我们,核实后会尽快下架及时删除,并可随时和客服了解处理情况,尊重保护知识产权我们共同努力。
2、文档的总页数、文档格式和文档大小以系统显示为准(内容中显示的页数不一定正确),网站客服只以系统显示的页数、文件格式、文档大小作为仲裁依据,平台无法对文档的真实性、完整性、权威性、准确性、专业性及其观点立场做任何保证或承诺,下载前须认真查看,确认无误后再购买,务必慎重购买;若有违法违纪将进行移交司法处理,若涉侵权平台将进行基本处罚并下架。
3、本站所有内容均由用户上传,付费前请自行鉴别,如您付费,意味着您已接受本站规则且自行承担风险,本站不进行额外附加服务,虚拟产品一经售出概不退款(未进行购买下载可退充值款),文档一经付费(服务费)、不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。
4、如你看到网页展示的文档有www.zixin.com.cn水印,是因预览和防盗链等技术需要对页面进行转换压缩成图而已,我们并不对上传的文档进行任何编辑或修改,文档下载后都不会有水印标识(原文档上传前个别存留的除外),下载后原文更清晰;试题试卷类文档,如果标题没有明确说明有答案则都视为没有答案,请知晓;PPT和DOC文档可被视为“模板”,允许上传人保留章节、目录结构的情况下删减部份的内容;PDF文档不管是原文档转换或图片扫描而得,本站不作要求视为允许,下载前自行私信或留言给上传者【自信****多点】。
5、本文档所展示的图片、画像、字体、音乐的版权可能需版权方额外授权,请谨慎使用;网站提供的党政主题相关内容(国旗、国徽、党徽--等)目的在于配合国家政策宣传,仅限个人学习分享使用,禁止用于任何广告和商用目的。
6、文档遇到问题,请及时私信或留言给本站上传会员【自信****多点】,需本站解决可联系【 微信客服】、【 QQ客服】,若有其他问题请点击或扫码反馈【 服务填表】;文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“【 版权申诉】”(推荐),意见反馈和侵权处理邮箱:1219186828@qq.com;也可以拔打客服电话:4008-655-100;投诉/维权电话:4009-655-100。