弯曲动脉内脉动流的有限元分析.PDF

第：?卷第2期 北京工业太学学报 vol二7No2 2001年6月 J0uRNALOFBElJlNGPOLYT阳HNIc0NIVERSITY Jun200】 弯曲动脉内脉动流的有限元分析 乔爱科，伍时桂，刘有军 (北京工业大学生物医学工程中心，北京 100022) 摘要： 给出了主动脉弓内的血液流动的理论模型．边界条件；根据生理脉动流条件，对主动脉弓内血液流动 进行有限元数值模拟．分析了计算结果及动脉粥拌硬化的流体动力学机理． 关键词： 动脉；生理脉动流；有限元分析；动脉粥样硬化 中圈分类号： Q66 文献标识码：A 文章编号：0254一0037(2001)02—024(】一昕 弯管动脉内的血液流动问题是生物流体力学的一个重要研究课题．对弯曲管内的流动问题的研究， n、00d Komai和KazuoTafIishi诅¨。等人的研究【I作，部 等¨1、chang等¨1、Berger等¨H1、Rindt，Y眦等¨1．Yutaka 以弯曲180。刚性圆管内的振荡流或脉动流为基础，来研究主动脉弓内的血液流动问题，而没有考虑升主 动脉对主动脉弓内血液流动的影响． 本文旨在对升主动脉和主动脉弓内发展中的生理性血液脉动流问题进行探讨．在主动脉弓的弯曲血 管部分，动脉粥样硬化、血管瘤等病变的发生和发展，将严重影响外周血管的血液流动特性…．因此．冉必 要对主动脉弓内的血液脉动流进行深入的研究，这对防治心血管疾病具有重要的现实意义．此外，研究主 动脉弓内的血液流动问题将为研究其他弯曲血管(如冠状动脉等)内血液流动问题打下良好的基础．且对 发展动脉中的脉搏波传播理论具有很重要的理论指导意义． 1 几何模型 主动脉阱左心室的主动脉瓣作为起点，开始是一段较为平直的 升主动脉，然后以一种复杂的三维空间方式作大约180。的弯曲(称 为主动脉弓)，并向心脏、头部和上肢伸展出分支；向下相当平直地穿 过隔膜直达腹部”1．为了便于分析，必须将升主动脉和主动脉弓的 几何形状进行简化． 图l是狗的升主动脉及主动脉弓示意图．作为数值计算的几何模 型，本文将升主动脉简化为长的直管，将主动脉弓简化为尼弯曲180。 的圆管，血管各截面的内半径为尺． 图1 升主动脉及主动脉弓，‘唐剖 2 数学模型 对于动脉血管内脉动流的动力学研究，存在很多的困难．主要表现在以下几个方面”…1：血液流动的 非线性、血管壁的非线性弹性(或非线性粘弹性)及有限变形、动脉血管的复杂几何形状，以及血管和血液 的脉动耦合等．由于问题的高度复杂性，作为研究的第1步，暂时仍将问题简化，假设血管为等圆截面刚 收稿日期：2000．04．】4． 基金项日：国家自然科学基金资助项目；国家教育部“高等学校骨干教师资助计划”基金资助项目 万方数据 第!期 乔爱科等：弯曲动脉内脉动流的有限元分析 24 怔甘’一皿敞刀小口J也湘十嘲批佧． 在升主动脉即直管段，采用圆柱坐标系下轴对称流动的守恒型血流控制方程…’． 在主动脉弓即弯曲血管段，采用如图2所示的圆环坐标系，建立守恒型的血流控制方程为如下． 。向动量：p{鲁+南[岳(r胁，¨+品c卅，+未(r～～)+胁，％一V；等】}= 一击考+F{V2～+寺(嚣c刚一嚣sin。一≥)} c，， 连续方程：岳(，m，)+嘉(饥)+丢(r～)=。 (4) 其中，H=凡。+rcos日；v。、v。、v，分别为血流轴向速度、周向速度和径向速度；p为压力；p为血液密度；一 嘶壮唐 边界条件： ①壁面无滑流条件