1、中国组织工程研究 第 18 卷 第 2 期 2014 01 08 出版 Chinese Journal of Tissue Engineering Research January 8, 2014 Vol.18, No.2 P.O. Box 10002, Shenyang 110180 www.CRTER.org 218 www.CRTER.org 吕绍茂，男， 1981 年生，湖南省永州市人，汉族，2007 年昆明医科大学毕业，硕士，主治医师，助理教授，主要从事影像解剖与应用研究。 通讯作者：段少银，博士，教授，主任医师 ， 厦门大学附属中山医院，福建省厦门市 361004 doi:10.3
2、969/j.issn.2095-4344. 2014.02.010 http:/www.crter.org 中图分类号 :R318 文献标识码 :B 文章编号 :2095-4344 (2014)02-00218-07 稿件接受 ： 2013-12-19 Lv Shao-mao, Master, Attending physician, Assistant professor, Zhongshan Hospital of Xiamen University, Xiamen 361004, Fujian Province, China Corresponding author: Duan Shao
3、-yin, M.D., Professor, Chief physician, Zhongshan Hospital of Xiamen University, Xiamen 361004, Fujian Province, China Accepted: 2013-12-19 构建颈内动脉瘤双向流固耦合模型的血流模拟 吕绍茂，钟 华，陈丽君，段少银 (厦门大学附属中山医院，福建省厦门市 361004) 文章亮点： 1 虚拟动脉瘤模型相对于动物模型及体外模型，具有相对简单、经济及较好的一致性 等特点 。得到的数据直观可视化，显示的相关参数很容易进行相关比较与分析，发现参数相互间的异同。 2 实
4、验采用最新开发的双向流固耦合方式 (ansys+fluent)的软件组合，功能上比 ansys+cfx 软件组合更强大，具有更强的处理大形变的网格能力。以往研究中采用刚性壁而且是稳态的研究为多， 忽略了血管壁的弹性作用以及作为瞬态的血流动力学变化。 3 以往的双向流固耦合研究采用动脉 壁 及动脉 瘤 统一厚度，而本实验采用动脉壁及动脉瘤分开定义厚度，得出的计算结果更加精确。 4 研究结果与临床观察或事实相吻合，瘤颈部具有很高的剪切力与压力，同时确认瘤顶部 有一 高压点，可能与动脉瘤的破裂有关，为动脉瘤基础与应用研究提供了依据。 关键词： 组织构建； 组织工程 ； 颈内动脉；后交通动脉瘤；有限元
5、建模；血流模拟 ； 国家自然科学基金 主题词： 颈内动脉；动脉瘤；有限元分析；血流动力学 基金资助： 国家自然科学基金面上项目 (81071214)资助 摘要 背景： 动脉瘤血液动力学变化为动脉瘤的治疗与其破裂的预防研究提供了前提与条件，有限元分析是一种很好的技术手段。 目的： 构建颈内动脉瘤的双向流固耦合模型，并进行相关血流模拟。 方法： 利用 GE Lightspeed 64 排螺旋 CT扫描获得颈内动脉瘤 CT 血管造影数据，分别在 mimics10.01 软件上进行三维模型的实体构建， ansys+fluent 软件进行流固耦合及血流模拟与仿真。 结果与结论： 构建了颈内动脉瘤瘤壁和流
6、场的有限元模型，模型具有良好的解剖形态，并与实体模型一致。模拟展 示整个心动周期的血流变化：瘤体血流呈漩涡样，瘤体部流速较瘤颈部慢；血流引起动脉瘤壁的形变、壁面剪切力、压力及应力以瘤颈部最大，瘤顶部最小。模拟值以 0.16 s(快速射血期 )时刻为最高值， 0.74 s(舒张期 )为最低值。结果显示基于 CT 血管造影扫描数据进行颈内动脉瘤的双向流固耦合建模，其方法简单、实用。模拟动脉瘤血流接近人体的生理状态，结果为研究动脉瘤发生、发展提供新的理论依据。 吕绍茂，钟华，陈丽君，段少银 . 构建颈内动脉瘤双向流固耦合模型的血流模拟 J.中国组织工程研究， 2014，18(2):218-224.
7、Blood flow simulation of internal carotid artery aneurysm using two-way flow-solid coupling method Lv Shao-mao, Zhong Hua, Chen Li-jun, Duan Shao-yin (Zhongshan Hospital of Xiamen University, Xiamen 361004, Fujian Province, China) Abstract BACKGROUND: The hemodynamics changes of aneurysm provide the
8、 basis for aneurysm treatment and ruptured aneurysm prevention, while finite element analysis is a good technical means. OBJECTIVE: To construct the two-way flow-solid coupling model of internal carotid artery aneurysm, and thus simulate its flows. METHODS: CT angiography data of internal carotid ar
9、tery aneurysm were recorded with GE Lightspeed 64 spiral CT scanning, and the corresponding three-dimension model was constructed with mimics10.01 software. Then the fluid-solid coupled flow simulations were done depending on Ansys+Fluent software. RESULTS AND CONCLUSION: The two-way flow-solid coup
10、ling model of internal carotid artery aneurysm was built, with the same morphology and anatomy as the three-dimensional CT imaging. In the whole cardiac cycle, the blood flow of aneurysm body was swirl typed and its velocity was slower than that of the aneurysm neck; the deformation, wall shear stre
11、ss, pressure and von Mises stress of aneurysm wall were the maximum at the neck, and the minimum at the top. Their highest value was at the 0.16 s (rapid ejection period), and their lowest value at the 0.74 s (relaxation period). On the basis of CT angiography data, to construct the two-way flow-sol
12、id coupling model of internal carotid artery aneurysm is a simple and practical method. The blood flow simulation of 吕绍茂，等 . 构建颈内动脉瘤双向流固耦合模型的血流模拟 ISSN 2095-4344 CN 21-1581/R CODEN: ZLKHAH 219 www.CRTER.org aneurysm is close to the human physiological conditions, and the results will provide new theo
13、retical basis for the study on the occurrence and development of aneurysm. Subject headings: carotid artery, internal; aneurysm; finite element analysis; hemodynamics Funding: the National Natural Science Foundation of China, No. 81071214 Lv SM, Zhong H, Chen LJ, Duan SY. Blood flow simulation of in
14、ternal carotid artery aneurysm using two-way flow-solid coupling method. Zhongguo Zuzhi Gongcheng Yanjiu. 2014;18(2):218-224. 0 引言 Introduction 动脉瘤是动脉壁上的异常膨出，大多发生在颅底区动脉，发病率为 1%-14%，是导致自发性蛛网膜下腔出血的最常见原因，具有很高的死亡率和致残率 1-4。目前的研究仍未阐明动脉瘤的发病机制，但动脉瘤多好发于血流动力有异常的血管分叉、弯曲等部位 4，所以血流动力 学被学者们公认为在动脉瘤的发生、发展、破裂过程中起着非常
15、重要的作用，血液动力学因素也逐渐成为影响动脉瘤发生发展以及破裂的研究热点 5-12。本文拟构建个性化颈内动脉后交通段动脉瘤模型，然后在此基础上实施瞬态双向流固耦合血流模拟分析，明确动脉瘤血液流动学特点，探讨其破裂的影响因素。 1 材料和方法 Materials and methods 设计： 探索性实验。 时间及地点： 实验于 2012年 6月至 2013年 6月在厦门大学附属中山医院影像科影像基础研究实验室完成。 材料： 利用医院的 PACS系统中 CT血 管造影扫描数据，选择颈内动脉后交通动脉瘤 3例，当时无破裂致蛛网膜下腔出血。扫描数据采用螺旋 CT(Light Speed VCT， G
16、E Inc，USA)获得，自主动脉弓至头顶采用增强连续扫描，扫描层厚 0.625 mm，图像矩阵 512512象 素，管电压 120 kV，管电流 250 mA 。 注 射对 比剂 使 用双 筒高 压 注射 器(STELLANT， MEDRAD Inc， USA)，非离子对比剂 (欧乃派克， 320 g/L， GEInc， USA)。利用三维建模软件mimics10.01 (MaterialiseInc， USA)，分别 建立颈内动脉交通段动脉瘤三维流场几何模型和管壁几何模型，通过geomagic studio 2012对几何模型进行修补，细分，平滑处理， 得 到表面较光滑的几何模型。流场几
17、何模型在软件ansys14.5 ICEM中划分体网格。 方法： 有限元模型的建立步骤 ： 图像数据编辑：颈内动脉交通段 CT血管造影扫描数据导入 Mimics图像三维重建软件，利用软件阈值工具(threshold)选择血管灰度值范围，形成原始图像。然后用动脉区域增长工具 (Dynamic region growing)在以上基础上对原始图像初步分割， 先使用图像编辑工具 (edit及 edit 3D mask)，手动对原始图像进行修补，使之尽量与周围结构区分出来，并使图像轮廓尽量平滑，填充图像 中的空隙 。最后利 用软件 的三维计 算功能 (3D Calculation)采用高质量模式对编辑好
18、的 mask进行三维重建，得到交通段的动脉瘤几何模型。用软件自带的remesh功能对几何模型适当处理，使模型尽量光滑，得到动脉瘤流场的几何模型，并确保两端轮廓线完整，再次使用 remesh功能去除几何模型的断面封口，得到动脉瘤壁的几何模型。 表面网格的优化：使用 geomagic studio 2012软件分别对动脉瘤流场的表面网格和动脉瘤壁的表面网格进行优化，修复质量较差的网格，重新划分网格使网格较均匀，并使用平滑处理和磨砂功能 得 到表面平滑、网格质量较高的动脉瘤流场和动脉瘤瘤壁的表面网格。 动脉瘤壁和动脉瘤场的体网格划分：动脉瘤流场的网格划分在 ICEM软件中进行，网格划分前要进行边界条
19、件的设置，如入口、出口、耦合壁等等，然后使用混合网格划分，主要为中央流体四面体网格和边界层的三棱柱体网格，具体参数要经过多次试验，找到最佳的网格大小。确保所有网格质量大于 0.3，网格质量越好，越有利于计算收 敛和保证计算结果准确性，最后输入体网格采用 fluent V6求解器。 计算方法 ： 流体部分属性及边界条件： 现有的研究， Aenis等 13计算指出，只要动脉直径大于 0.5 mm，则用牛顿流体代替非牛顿流体所引起的误差不超过 2%。而本文的研究动脉直径的入口为 4.8 mm，出口直径也有 3.8 mm，均大于 0.5 mm，所以在计算中采用牛顿流体代替非牛顿流体，以简化模型。设置血
20、液密度为 1 055 kg/m3，血液黏滞度为0.003 5 Pas，收集到的颈内动脉入口血流最大速度约 0.79 m/s，通过公式 Re=VD/， 为流体密度， V为圆管内直径， 为流体的黏度系数，计算得出的雷诺数约 948，一般管道雷诺数 Re4 000为紊乱状态， Re=2 000-4 000为过渡状态，所以血液流动的方式是层流，故本文采用层流模型，研究中将血液流动模拟为弹性管的绝热、不可压缩和非定常的牛顿流体。血液和血管壁面无滑动，此患者颈内动脉起始处入口速度随时间变化曲线由超声多普勒测得 ， 心动周期为 0.74 s(图 1)。 血管壁属性及边界条件： 血管壁的形变问题是极其复杂的非
21、线性过程，为 了简化计 算量，假设动脉壁与动脉瘤壁为线弹性材料，密度为 =1 150 kg/m3，文献报道的血管组织弹性模量在几十 kPa到几千 kPa之间 14，并没有确切的值，本文弹性模量采用 E=50 kPa，泊松比 =0.45(近似于吕绍茂，等 . 构建颈内动脉瘤双向流固耦合模型的血流模拟 P.O. Box 10002, Shenyang 110180 www.CRTER.org 220 www.CRTER.org 不可压缩材料 )， 文献 15中健康的颈内动脉壁的厚度不超过 0.5 mm，动脉瘤壁也肯定比动脉壁要薄，所 以本文设置动脉壁的厚度为 0.3 mm，文献报道动脉瘤壁厚度 为
22、 0.08-0.2 mm16，故动脉瘤壁厚度为 0.15 mm。血管入口至颈内动脉床突段都位于硬脑膜外，周围结构较多，故采用固定的 方式，出口固定位移为 0，入口采用速度为边界条件，出口使用压力为 0。 流固耦合求解设置： 采用 ansys公司最新版的 14.5版本，利用 workbench和 fluent进行双向流固耦合，采用压力 -速度耦合方式，流场计算时的微分方程组离散格式为二阶精度的迎风格式，使用 10-4均方根残差设定标准。流场和固定计算所取时间步长为 0.002 s。使用处理器为 intel e3 1 230 v2，内存 8 G，经超过 2个周期的迭代运算，耗时 20 h，得到稳定
23、的收敛结果。 2 结果 Results 2.1 模型比较 构 建了颈内后动脉交通段动脉瘤的壁和流场的有限元模型，模型具有良好的解剖形态，并与实体模型一致 (图 2)。 2.2 血液流动情况 截取中间第 2个周期几个典型时刻的血管腔内的流线图及速度体积图，分别是 0.16 s(入口瞬时速度为 94.7 cm/s，心动周期内的最大速度 )； 0.36 s(入口瞬时速度为 53.2 cm/s，为峰谷时刻 )， 0.44 s(入口速度为58.19 cm/s，为第二峰速时刻 )， 0.74 s(舒张期末，心动周期内最小入口速度，为 37.2 cm/s)。模拟图像可以发现 (图3)，在整个心动周期内，动
24、脉瘤腔内均存在较为明显漩涡流动。当血流经动脉瘤远颈部时分成两股血流，一股血流沿着载瘤动脉直接下行至大脑前、大脑中动脉，为主干血流，另一股较少血流流入动脉瘤体内，瘤体内形成逆时针的漩 涡，最后加入主干血流；漩涡流动的类型在整个心动周期内比较稳定， 各个时刻流线图及速度体积图形显示远心端瘤颈部血流速度最快，靠近近心端瘤颈部区域也有小部分相对较快的血流，瘤腔内的其他部分血流相对较慢 (图3)。 2.3 血管壁的形变情况 显示的对整个心动周期的血液流动观察后，动脉瘤壁的变形情况，随着心脏的搏动，收缩期整个动脉瘤扩张， 舒张期恢复，形变的规律随着流场的变化而变化，整个动脉瘤在 0.16 s时刻的变形最大
25、，以动脉瘤远心端颈部血管壁变形量最明显，变形约为 0.27 mm，因为此区域直接受到血流的冲击作用，而动脉瘤两侧壁的形变相对较小 (图 4)。 2.4 血管壁面剪切力 (WSS) 选取第 2个心动周期，作为观察动脉瘤的壁面剪切力的变化情况，取 4个标志性点： a点 (位于远心端瘤颈部 )， b、 c点 位于 远心端瘤至瘤顶中间，d点位于瘤顶。整个心动周期内 a点的壁面剪切力变化范围58.2-246.5 Pa， b点为 20.2-115.0 Pa， c点为 4.8-50.1 Pa，d点的壁面剪切力的变化范围为 0.6-15.3 Pa。观察发现因为受到血流的直接冲击， a点区域的壁面剪切力最大，并
26、且图中可以看出 a点区域的壁面剪切力梯度最大，随着漩涡流动到瘤顶 d点区域时，剪切力最小 (图 5)。 2.5 血液对管壁的压力 显示 4个点 a、 b、 c、 d位置的压力及不同时刻对管壁的压力曲线图 (测定点的位置与壁面剪切力模拟相同 )，心动周期内， a点的压力变化范围为 90.4-365.1 Pa ， b 点的压力变化范围为 161.1- 679.2 Pa， c点的压力变化范围 为 263.5-1 534 Pa， d点压力变化范围为 363.2-3 143 Pa；同时还发现远心端瘤颈 a点往后 e点区域压力较大， e点的压力变化范围为538.8-3 174 Pa。图中和数据显示，远心端
27、瘤颈部是血流直接冲击区域，在整个心动周期内压力变化最剧烈，并且在速度最大时刻附近达到最大值，压力梯度也处于较高水平，随着位置向近心端瘤颈部及两侧壁移动，血流速度趋缓，压力变化范围幅度变小，峰值也减小，压力梯度走低，处于低水平 (图 5)。 图 1 上 图为超声多普勒测 得颈内动脉 (入口处 )的瞬时平均速度， 下 图为 GetData软件获得 上 图第 4个心动周期的速度时间曲线 Figure 1 The instantaneous average speed was measured by Doppler ultrasound at the entrance of internal caro
28、tid artery (upper), the curve of velocity-time from the fourth cardiac cycle was obtained by the software of GetData (below) 吕绍茂，等 . 构建颈内动脉瘤双向流固耦合模型的血流模拟 ISSN 2095-4344 CN 21-1581/R CODEN: ZLKHAH 221 www.CRTER.org 图 2 显示动脉瘤形态与位置 Figure 2 Morphology and location of aneurysm 图注： 图中 A 为 CTA 重建后颈内后交 通段动
29、脉瘤三维图像 ， B 为 几何 模型 ， C 为 有限元模型 。 图 3 0.16 s(A), 0.36 s(B),0.44 s(C),0.74 s (D)时刻流线图 Figure 3 The time-line images of 0.16 s (A), 0.36 s (B), 0.44s (C) and 0.74 (D) 图注： 0.16 s 时刻显示血流速度最快，动脉瘤内涡流明显 。 A B C A B C D 图 4 动脉瘤 4 个典型时刻变形 Figure 4 Four typical time-deformation images of aneurysm 图注： 图中显示 A 为
30、0.16 s 最大形变约 0.27 mm， B 为 0.36 s 最大变形约 0.1 mm， C 为 0.44 s 最大变形约 0.13 mm， D 为 0.74 s 最大变形约 0.05 mm A B C D 图 5 血管壁面剪切力与血液对管壁的压力 Figure 5 The wall shear and wall pressure of blood vessels 图注： 图中 A 为 心动周期内，不同区域 4 个点的壁面剪切力分布图， B 为时刻剪切力曲线图。显示瘤颈远心侧 (a 点 )剪切力最大。 C 为 心动周期内，上述 4 个点加 e 点 (瘤 颈远心侧偏后 )的壁面压力分布图，
31、D 为时刻压力曲线图。显示瘤顶 (d点 )及瘤颈远心侧偏后 (e 点 )压力较大。 a b c d a b c d a b c d e a b c d e A B C D 吕绍茂，等 . 构建颈内动脉瘤双向流固耦合模型的血流模拟 P.O. Box 10002, Shenyang 110180 www.CRTER.org 222 www.CRTER.org 2.6 管壁面应力分布情况 根据动脉瘤壁的特点及参考相关文献 17, 本文采用 von Mises stress作为壁的研究应力。 Von Mises stress是基于剪切应变能的一种等效应力，其大概的含义是当单元的形状改变比 能达到一定程
32、度，材料开始屈服。模拟结果 Von Mises stress 值大小在 10-2-10-4 Pa之间。不同时刻的入口速度不同，血管壁内的应力大小也不相同，分析各个时刻的不同位置 von Mises stress，心动周期内远心端瘤颈部具有较大的应力，在 0.16 s的时候 von Mises stress最大，而两侧瘤体中央至顶部及近心端瘤颈部在各个时候都拥有比较小的 von Mises stress(图 6)。 3 讨论 Discussion 3.1 研究意义 颅内动脉瘤破裂出血具有较高的致残 率和死亡率，因而研究动脉瘤的发生、发展及破裂的机制甚是重要。越来越 多 的研究资料说明，血液动力学
33、在动脉瘤的形成、发展及破裂中占据重要的位置。由于种种的限制 ,实际测量动脉瘤的壁面剪切力存在诸多困难 ， 然而计算流体力学提供了很好的方法 ， 但是前期的研究多基于简化的理想模型，也得出了一些有益的结果 ， 但与真实的个人化模型下的结果存在一定的差距。很多模型的分析都是基于刚性壁的研究 18-22， Perktold等 23采用有限元法分析了弹性颈动脉分支血管的脉动血流，研究发现，不同于刚性管壁，分支旁的区域流动分离有所增大； 壁面剪切力大小降低了 25%；顾媛等 24对弹性管壁及刚性壁两种条件下对比研究，弹性血管内的壁面剪切力，要比刚性壁内的减小明显。随着科学技术的发展 ， 特别在计算机不断
34、进步的情况，采用双向耦合分析动脉瘤血液流动和壁面应力的分布将更接近真实的人体生理状态，有利于观察动脉瘤的血流动力学。同时，为研究动脉瘤血流与瘤壁的相互作用机制提供另一种可靠手段 25-29。 3.2 动脉瘤血流特点 颅内动脉瘤的发生与血流的各种生理因素密切相关，动脉瘤常发生在血管发叉或明显弯曲的血管处，这说明与血流动力学的分布特征有关。血流 的冲击会形成两个方向的力，一个是跨比压力，垂直作用于血管壁。另一个即血管壁面剪切力，是血流对血管壁的黏性摩擦力，平行作用于血管壁。最近的研究表明，壁面剪切力在动脉瘤的发生、形成过程中起着主要因素的作用，血管壁的最内层及血管内皮细胞对血管壁面剪切力最敏 感
35、30-31，而且近年来研究认为动脉瘤的生长和扩大与瘤内低血管壁面剪切力相关 32-34。 本文的研究结果显示，瘤颈部壁面剪切力最高，流入后局部管壁逐渐降低，瘤顶壁面剪切力达到最低水平。这与低血流理论相吻合，即低血管壁面剪切力可导致动脉壁的退行性改变。 有作者认为小于 1.5 Pa的血管壁面剪切力会通过细胞凋亡过程导致内皮细胞退变，进而，局部的血管壁面剪切力不足导致动脉瘤破裂 8, 35。临床观察大多数动脉瘤破裂的位置都多位于瘤顶区 36-38，可能正是由于长期处于低剪切力的状态，使得血管内皮损伤和凋亡，最终导致动脉瘤破裂的原因。另外，瘤顶区发现有一高压力范围，在流体内静压力不变的情况下，局部高
36、压力可能引发动脉瘤的破裂。 图 6 显示典型时刻 von Mises stress 分布图 Figure 6 von Mises stress distribution at typical time 图注：图中 A， B 为 0.16 s(前、后面观 )， C， D 为 0.36 s(前、后面观 )， E， F 为 0.44 s(前、后面观 ) ， G， H 为 0.74 s(前、后面观 )， 其中 0.16 s 时刻图显示 von Mises stress 最大位于瘤颈后部。 A B C D E F G H 吕绍茂，等 . 构建颈内动脉瘤双向流固耦合模型的血流模拟 ISSN 2095-43
37、44 CN 21-1581/R CODEN: ZLKHAH 223 www.CRTER.org 3.3 应力与动脉瘤的关系 作者考虑了 血管 壁对流场的作用，对建立的血流及血管壁进行双向流固耦合分析，可以对管壁的应力场和形态变化进 行研究，通过迭代计算出了各个位置处壁面应力随着入口速度随时间变化的关系，结果显示了 Von mises stress沿着瘤顶到瘤颈部方向，呈现逐渐增大趋势。瘤壁的形变程度与 Von mises应力变化基本一致，瘤颈部最大，而瘤顶部形变相对较小。从理论上分析，应力高的地方、材料抗损坏的能力强，这就能解释因为动脉瘤颈部的压力较大，导致瘤颈部及载瘤动脉的壁增厚 39；而瘤
38、顶部的压力较小，所以壁的增厚没有瘤颈部明显，加上瘤顶部低血管壁面剪切力介导血管内皮退变、凋亡，所以瘤顶部更容易破裂 11。 3.4 研究特色与 展望 采用个体化的入口血流速度，分别对动脉瘤壁和载瘤动脉瘤壁分别设定厚度，使用workbench和 fluent双向流固耦合，使得结果更加接近人体真实生理的血液动力学 40-42。当然本研究也做了一些假设条件具有一定的局限性，血管壁和动脉瘤壁设定为同一种线弹性材料，管壁厚度均匀，动脉瘤厚度也 均匀，使得计算结果和真实情况有所差别，其中比较明显的是：动脉瘤壁相对于载瘤动脉的管壁而言，其刚性应该更大。血液动力学还受到动脉瘤周围环境等诸多因素的影响，如脑组织
39、，周围的结缔组织， 脑脊液等情况。 在今后的研究工作中应该 进一步研究动脉瘤和载瘤动脉 的材质特点，再增加更多的研究方法与技术并进行相关参数的验证，将会产生更大的应用价值 43-48。随着力学理论的完善和材料试验的相关研究，相信通过血液动力学的模拟仿真能够为动脉瘤的基础与临床研究提供更多的方法与技术。 结论 ： 基于 CT血管造影 (CT angiography CTA)扫描数据，能够很好的获得 颈内动脉瘤的双向流固耦合有限元模型。血流模拟具有更接近人体动脉血流状态，获得的相关数据具有很好的应用价值。相信进一步的模型完善与功能模拟将为研究动脉瘤的发生、发展提供新的手段。 致 谢 ： 感谢厦门大
40、学附属中山医院超声科进行颈内动脉血流的测定与验证。 作者贡献 ： 实验设计为段少银，实验实施为吕绍茂、钟华，实验评估为段少银、吕绍茂、钟华，资料收集为吕绍茂、钟华、陈丽君。吕绍茂成文，段少银审校，段少银、吕绍茂、钟华对文章负责。 利益冲 突 ： 文章及内容不涉及相关利益冲突。 伦理要求 ： 没有与相关伦理道德冲突的内容。 学术术语 ： 流固耦合 -是流体与固体之间的交互作用。流固耦合可简单分为单向及双向耦合，相关的计算软件有 ADINA、ANSYS、 ESI-CFD 以及 COMSOL 等。单向耦合常见的有热 -结构耦合分 析，而双向耦合常见的则有流 -结构耦合分析。 作者声明 ： 文章为原创
41、作品，无抄袭剽窃，无泄密及署名和专利争议，内容及数据真实，文责自负。 4 参考文献 References 1 梁长虹 ,赵振军 .多层螺旋 CT血管成像 M.北京 :人民军医出版社 ,2008:121-122. 2 薛庆澄 .神经外科学 M.天津 :天津科学技术出版社 ,1990: 350-370. 3 Pierot L, Spelle L, Vitry F,et al.Similar safety in centers with low and high volumes of endovascular treatments for unruptured intracranial aneury
42、sms: evaluation of the analysis of treatment by endovascular approach of nonruptured aneurysms study.AJNR Am J Neuroradiol. 2010;31(6): 1010-1014. 4 Alnaes MS, Isaksen J, Mardal KA,et al.Computation of hemodynamics in the circle of Willis.Stroke. 2007;38(9): 2500-2505. 5 Cebral JR, Sheridan M, Putma
43、n CM.Hemodynamics and bleb formation in intracranial aneurysms.AJNR Am J Neuroradiol. 2010;31(2):304-310. 6 Karmonik C, Yen C, Diaz O,et al.Temporal variations of wall shear stress parameters in intracranial aneurysms-importance of patient-specific inflow waveforms for CFD calculations.Acta Neurochi
44、r (Wien). 2010;152(8): 1391-1398. 7 于红玉 ,李海云 ,张莹 ,等 .颅内动脉瘤 3D几何模型建立方法的研究 J.计算机工程与应用 ,2008,44(17):175-177. 8 Meng H, Wang Z, Hoi Y, et al.Complex hemodynamics at the apex of an arterial bifurcation induces vascular remodeling resembling cerebral aneurysm initiation.Stroke. 2007;38(6): 1924-1931. 9 陈珍
45、 ,袁奇 ,朱光宇 ,等 .颈内动脉狭窄时 Willis环的血液动力学数值模拟 J.西安交通大学学报 ,2009,43(10):119-123. 10 Castro MA, Putman CM, Sheridan MJ,et al.Hemodynamic patterns of anterior communicating artery aneurysms: a possible association with rupture.AJNR Am J Neuroradiol. 2009;30(2):297-302. 11 Cebral JR, Mut F, Weir J,et al.Quanti
46、tative characterization of the hemodynamic environment in ruptured and unruptured brain aneurysms.AJNR Am J Neuroradiol. 2011;32(1): 145-151. 12 Goubergrits L, Schaller J, Kertzscher U,et al.Statistical wall shear stress maps of ruptured and unruptured middle cerebral artery aneurysms.J R Soc Interf
47、ace. 2012;9(69):677-688. 13 Aenis M, Stancampiano AP, Wakhloo AK,et al.Modeling of flow in a straight stented and nonstented side wall aneurysm model.J Biomech Eng. 1997;119(2):206-212. 14 Tsuzuki K, Hasegawa H, Ichiki M,et al.Optimal region-of- interest settings for tissue characterization based on
48、 ultrasonic elasticity imaging.Ultrasound Med Biol. 2008;34(4): 573-585. 15 Leach JR, Rayz VL, Soares B,et al.Carotid atheroma rupture observed in vivo and FSI-predicted stress distribution based on pre-rupture imaging.Ann Biomed Eng. 2010;38(8):2748-2765. 16 Frsen J, Piippo A, Paetau A,et al.Remodeling of saccular cerebral artery aneurysm wall is associated with rupture: histological analysis of 24 unruptured and 42 ruptured cases.Stroke. 2004;35(10):2287-2293. 17 Vorp DA, Van
