Large-Scale 3-D Geometric Reconstruction of the Porcine Coronary Arterial Vasculature Based on Detailed Anatomical Data

The temporal and spatial distribution of coronary blood flow, pressure, and volume are determined by the branching pattern and three-dimensional (3-D) geometry of the coronary vasculature, and by the mechanics of heart wall and vascular tone. Consequently, a realistic simulation of coronary blood flow requires, as a first step, an accurate representation of the coronary vasculature in a 3-D model of the beating heart. In the present study, a large-scale stochastic reconstruction of the asymmetric coronary arterial trees (right coronary artery, RCA; left anterior descending, LAD; and left circumflex, LCx) of the porcine heart has been carried out to set the stage for future hemodynamic analysis. The model spans the entire coronary arterial tree down to the capillary vessels. The 3-D tree structure was reconstructed initially in rectangular slab geometry by means of global geometrical optimization using parallel simulated annealing (SA) algorithm. The SA optimization was subject to constraints prescribed by previously measured morphometric features of the coronary arterial trees. Subsequently, the reconstructed trees were mapped onto a prolate spheroid geometry of the heart. The transformed geometry was determined through least squares minimization of the related changes in both segments lengths and their angular characteristics. Vessel diameters were assigned based on a novel representation of diameter asymmetry along bifurcations. The reconstructed RCA, LAD and LCx arterial trees show qualitative resemblance to native coronary networks, and their morphological statistics are consistent with the measured data. The present model constitutes the first most extensive reconstruction of the entire coronary arterial system which will serve as a geometric foundation for future studies of flow in an anatomically accurate 3-D coronary vascular model.     

锁骨下动脉、颈总动脉和椎动脉分叉处血流动力学数值模拟

目的:探讨锁骨下动脉、颈总动脉和椎动脉分叉处的血流动力学特性,分析该处发生血管狭窄引起大脑供血不足的 血流动力学原因。方法:采用内蒙古民族大学附属医院神经内科提供的CT数据,应用医学建模软件MIMICS20.0将患者 二维CT数据进行三维血管重建,经过网格划分及边界条件设置后导入计算流体力学软件FLUENT14.5中。计算和分析 不同血液入口速度的锁骨下动脉、颈总动脉和椎动脉分叉处的血流动力学特性。结果:在血液入口速度不同的情况下,锁 骨下动脉、颈总动脉和椎动脉分叉处的血液流场分布、血液压力分布和血管壁面切应力分布有显著变化。在血液入口速 度增大时,锁骨下动脉分叉处和颈总动脉分叉处的血液流速快、血管壁压力大,颈总动脉内侧血管壁面切应力大,但锁骨 下动脉分叉处和颈总动脉分叉处血管壁面切应力数值和变化幅度小,属于低切应力区。结论:通过血流动力学数值模拟 研究,分析锁骨下动脉、颈总动脉和椎动脉分叉处易发生粥样斑块病变导致大脑供血不足的血流动力学原因。     

Methods for Three-Dimensional Geometric Characterization of the Arterial Vasculature

Complex vascular anatomy often affects endovascular procedural outcome. Accurate quantitative assessment of three-dimensional (3D) in-vivo arterial morphology is therefore vital for endovascular device design, and preoperative planning of percutaneous interventions. The aim of this work was to establish geometric parameters describing arterial branch origin, trajectory, and vessel curvature in 3D space that eliminate the errors implicit in planar measurements. 3D branching parameters at visceral and aortic bifurcation sites, as well as arterial tortuosity were determined from vessel centerlines derived from magnetic resonance angiography data for three subjects. Errors in coronal measurements of 3D branching angles for the right and left renal arteries were 3.1 ± 3.4° and 7.5 ± 3.7°, respectively. Distortion of the anterior visceral branching angles from sagittal measurements was less pronounced. Asymmetry in branching and planarity of the common iliac arteries was observed at aortic bifurcations. The renal arteries possessed considerably greater 3D curvature than the abdominal aorta and common iliac vessels with mean average values of 0.114 ± 0.015 and 0.070 ± 0.019 mm⁻¹ for the left and right, respectively. In conclusion, planar projections misrepresented branch trajectory, vessel length, and tortuosity proving the importance of 3D geometric characterization for possible applications in planning of endovascular interventional procedures and providing parameters for endovascular device design.     

冠脉树三维重建中几何变换矩阵的优化

冠脉造影图像间的几何变换矩阵是实现冠脉树三维重建的关键。本研究分析了优化几何变换矩阵的必要性，提出利用分支点坐标、分支血管方向矢量和分支夹角三类数据优化几何变换矩阵的方法，并给出了kvenberg-Marquardt算法的优化步骤。最后利用两幅未标定的单面冠脉造影图像，实现冠脉树骨架的三维重建。重建的误差统计结果表明，优化方法有效地提高了冠状动脉树三维重建的精度。     

Physical Factors Effecting Cerebral Aneurysm Pathophysiology

Many factors that are either blood-, wall-, or hemodynamics-borne have been associated with the initiation, growth, and rupture of intracranial aneurysms. The distribution of cerebral aneurysms around the bifurcations of the circle of Willis has provided the impetus for numerous studies trying to link hemodynamic factors (flow impingement, pressure, and/or wall shear stress) to aneurysm pathophysiology. The focus of this review is to provide a broad overview of such hemodynamic associations as well as the subsumed aspects of vascular anatomy and wall structure. Hemodynamic factors seem to be correlated to the distribution of aneurysms on the intracranial arterial tree and complex, slow flow patterns seem to be associated with aneurysm growth and rupture. However, both the prevalence of aneurysms in the general population and the incidence of ruptures in the aneurysm population are extremely low. This suggests that hemodynamic factors and purely mechanical explanations by themselves may serve as necessary, but never as necessary and sufficient conditions of this disease’s causation. The ultimate cause is not yet known, but it is likely an additive or multiplicative effect of a handful of biochemical and biomechanical factors.     

Structural and biophysical simulation of angiogenesis and vascular remodeling.

The purpose of this report is to introduce a new computer model for the simulation of microvascular growth and remodeling into arteries and veins that imitates angiogenesis and blood flow in real vascular plexuses. A C++ computer program was developed based on geometric and biophysical initial and boundary conditions. Geometry was defined on a two-dimensional isometric grid by using defined sources and drains and elementary bifurcations that were able to proliferate or to regress under the influence of random and deterministic processes. Biophysics was defined by pressure, flow, and velocity distributions in the network by using the nodal-admittance-matrix-method, and accounting for hemodynamic peculiarities like Fahraeus-Lindqvist effect and exchange with extravascular tissue. The proposed model is the first to simulate interdigitation between the terminal branches of arterial and venous trees. This was achieved by inclusion of vessel regression and anastomosis in the capillary plexus and by remodeling in dependence from hemodynamics. The choice of regulatory properties influences the resulting vascular patterns. The model predicts interdigitating arteriovenous patterning if shear stress-dependent but not pressure-dependent remodeling was applied. By approximating the variability of natural vascular patterns, we hope to better understand homogeneity of transport, spatial distribution of hemodynamic properties and biomass allocation to the vascular wall or blood during development, or during evolution of circulatory systems.     

Three-dimensional arrangement of the vasa vasorum in explanted segments of the aged human great saphenous vein: scanning electron microscopy and three-dimensional morphometry of vascular corrosion casts

The vasa vasorum of skeletonized and nonskeletonized segments of five human great saphenous veins (GSVs), harvested during coronary bypass grafting, were cannulated, rinsed, and injected (casted) with the polymerizing resin Mercox-Cl-2B. After removal of the dry vascular tissue, the casts were examined using scanning electron microscopy. Stereopaired images (tilt angle, 6 degrees ) were taken, imported into a 3D morphometry system, and the 3D architecture of the vasa vasorum (arterial and venous vasa as well as capillaries) was studied qualitatively and quantitatively in terms of vasa diameters, intervascular and interbranching distances, and branching angles. Diameters of parent (d(0)) and large (d(1)) and small (d(2)) daughter vessels of arterial and venous bifurcations served to calculate asymmetry ratios (alpha) and area ratios (beta). Additionally, deviations of bifurcations and branching angles from optimal branches were calculated for selected arterial vasa. The arrangement of the vasa vasorum closely followed the longitudinally oriented connective tissue fibers in the adventitia and the circularly arranged smooth muscle cell layers within the outer layers of the media. Venous vasa by far outnumbered arterial vasa. Vasa vasorum changed their course several times in acute angles and revealed numerous circular constrictions, kinks, and outpouchings. Due to their spatial arrangement, the vasa vasorum are prone to tolerate vessel wall distension generated by acute increases in blood pressure or stretching of the vessel without severe impact on vessel functions. Preliminary comparisons of data from the bifurcations of cast arterial vasa vasorum, with calculated optimal bifurcations, do not yet give clear insights into the optimality principle(s) governing the design of arterial vasa vasorum bifurcations of the human GSVs.     

3D reconstruction of the cerebral arterial network from stereotactic DSA.

The authors present an automatic algorithm for 3D reconstruction of cerebral blood vessels by digital subtracted angiography. The patient is localized by a stereotactic method. The reconstruction algorithm includes two steps: first vessel extraction then 2D matching and reconstruction. Accurate vessel skeletons are generated by a combination of mathematical morphological algorithms and adaptive filters. The 3D reconstruction algorithm is based on the reconstruction of vessels center lines. For that purpose, three different projections of the vascular network are used. Reconstruction is computed segment by segment (a curved line between two nodes). For each segment point, the algorithm defines all epipolar solutions on the other views. These epipolar solutions are sorted and pooled by 2D continuity and 3D proximity criteria resulting in a 3D graph. Optimal 3D segment is defined by a recursive algorithm that looks up the better path in the 3D graph. The algorithms have been implemented on a Compatible-PC computer in C language. More than 95% of static copper phantom was reconstructed in 5 min and with 1 mm 3D accuracy. 70% of arteries (from carotid to the seventh node) of a true patient arterial network were reconstructed is less than 30 min.     

Elevated cerebral blood flow and vascular density in the amygdala after status epilepticus in rats

Cerebrovascular changes following status epilepticus (SE) are not well understood, yet they may contribute to epileptogenesis. We studied hemodynamic changes in the cerebral cortex and amygdala by arterial spin labeling (ASL) and dynamic susceptibility contrast (DSC) MRI at 2 days and 14 days after pilocarpine-induced SE in rats. There were no cortical hemodynamic changes, yet in the amygdala we found prolonged elevation in cerebral blood flow (CBF, 129% of control mean, day 14, p < 0.01). There was a trend towards increased cerebral blood volume (CBV) during the same imaging sessions. Through immunohistochemistry, we observed increased vessel density in the amygdala (127% of control mean, day 14, p < 0.05). In conclusion, epileptogenesis may involve hemodynamic changes that are associated with vascular reorganization during post-SE remodeling in the amygdaloid complex.     

Optimality in the developing vascular system: branching remodeling by means of intussusception as an efficient adaptation mechanism.

The theory of bifurcating vascular systems predicts vessel diameters that are related to optimality criteria like minimization of pumping energy or of building material. However, mechanisms for producing the postulated optimality have not been described so far, and quantitative data on bifurcation diameters during development are scarce. We used an embryonic vascular bed that rapidly grows and adapts to changing hemodynamic conditions, the chicken chorioallantoic membrane (CAM), and correlated vascular cast and tissue section morphology with in vivo time-lapse video monitoring. The bifurcation exponent delta and associated parameters were quantitatively assessed in arterial and venous microvessels ranging in diameter from 30 to 100 microm. We observed emergence of optimality by means of intussusception, i.e., formation of transvascular tissue pillars. In addition to intussusceptive microvascular growth (IMG = expansion of capillary networks) and intussusceptive arborization (IAR = formation of feeding vessels from capillaries) the observed intussusception at bifurcations represents a third variant of nonsprouting angiogenesis. We call it intussusceptive branching remodeling (IBR). IBR occurred in vessels of considerable diameter by means of two alternative mechanisms: either through pillars arising close to a bifurcation, which increased in girth until they merged with the connective tissue in the bifurcation angle; or through pillars arising at some distance from the bifurcation point, which then expanded by formation of ingrowing tissue folds until they became connected to the tissue of the bifurcation angle. Morphologic evidence suggests that IBR is a wide-spread phenomenon, taking place also in lung, intestinal, kidney, eye, etc., vasculature. Irrespective of the mode followed, IBR led to a branching pattern close to the predicted optimum, delta = 3.0. Significant differences were observed between delta at arterial bifurcations (2.70 to 2.90) and delta at venous bifurcations (2.93 to 3.75). IBR, by means of eccentric pillar formation and fusion, was also involved in vascular pruning. Experimental changes in CAM hemodynamics (by locally increasing blood flow) induced onset of IBR within less than 1 hr. Our study provides morphologic and quantitative evidence that a similar cellular machinery is used for all three variants of vascular intussusception, IMG, IAR, and IBR. It thus provides a mechanism of efficiently generating complex blood transport systems from limited genetic information. Differential quantitative outcome of IBR in arteries and veins, and the experimental induction of IBR strongly suggest that hemodynamic factors can instruct embryonic vascular remodeling toward optimality.     

Branching Pattern of the Cerebral Arterial Tree

Quantitative data on branching patterns of the human cerebral arterial tree are lacking in the 1.0–0.1 mm radius range. We aimed to collect quantitative data in this range, and to study if the cerebral artery tree complies with the principle of minimal work (Law of Murray). To enable easy quantification of branching patterns a semi-automatic method was employed to measure 1,294 bifurcations and 2,031 segments on 7 T-MRI scans of two corrosion casts embedded in a gel. Additionally, to measure segments with a radius smaller than 0.1 mm, 9.4 T-MRI was used on a small cast section to characterize 1,147 bifurcations and 1,150 segments. Besides MRI, traditional methods were employed. Seven hundred thirty-three bifurcations were manually measured on a corrosion cast and 1,808 bifurcations and 1,799 segment lengths were manually measured on a fresh dissected cerebral arterial tree. Data showed a large variation in branching pattern parameters (asymmetry-ratio, area-ratio, length-radius-ratio, tapering). Part of the variation may be explained by the variation in measurement techniques, number of measurements and location of measurement in the vascular tree. This study confirms that the cerebral arterial tree complies with the principle of minimum work. These data are essential in the future development of more accurate mathematical blood flow models. Anat Rec, 302:1434–1446, 2019. © 2018 The Authors. The Anatomical Record published by Wiley Periodicals, Inc. on behalf of American Association of Anatomists.     

A Novel Simulation Strategy for Stent Insertion and Deployment in Curved Coronary Bifurcations: Comparison of Three Drug-Eluting Stents

The introduction of drug-eluting stents (DES) has reduced the occurrence of restenosis in coronary arteries. However, restenosis remains a problem in stented coronary bifurcations. This study investigates and compares three different second generation DESs when being implanted in the curved main branch of a coronary bifurcation with the aim of providing better insights into the related changes of the mechanical environment. The 3D bifurcation model is based on patient-specific angiographic data that accurately reproduce the in vivo curvatures of the vessel segments. The layered structure of the arterial wall and its anisotropic mechanical behavior are taken into account by applying a novel algorithm to define the fiber orientations. An innovative simulation strategy considering the insertion of a folded balloon catheter over a guide wire is proposed in order to position the stents within the curved vessel. Straightening occurs after implantation of all stents investigated. The resulting distributions of the wall stresses are strongly dependent on the stent design. Using a parametric modeling approach, two design modifications, which reduce the predicted maximum values of the wall stress, are proposed and analyzed.     

Simulation of abdominal aortic aneurysm growth with updating hemodynamic loads using a realistic geometry

Advances in modeling vascular tissue growth and remodeling (G&R) as well as medical imaging usher in a great potential for integrative computational mechanics to revolutionize the clinical treatment of cardiovascular diseases. A computational model of abdominal aortic aneurysm (AAA) enlargement has been previously developed based on realistic geometric models. In this work, we couple the computational simulation of AAA growth with the hemodynamics simulation in a stepwise, iterative manner and study the interrelation between the changes in wall shear stress (WSS) and arterial wall evolution. The G&R simulation computes a long-term vascular adaptation with constant hemodynamic loads, derived from the previous hemodynamics simulation, while the subsequent hemodynamics simulation computes hemodynamic loads on the vessel wall during the cardiac cycle using the evolved geometry. We hypothesize that low WSS promotes degradation of elastin during the progression of an AAA. It is shown that shear stress-induced degradation of elastin elevates wall stress and accelerates AAA enlargement. Regions of higher expansion correlate with regions of low WSS. Our results show that despite the crucial role of stress-mediated collagen turnover in compensating the loss of elastin, AAA enlargement can be accelerated through the effect of WSS. The present study is able to account for computational models of image-based AAA growth as well as important hemodynamic parameters with relatively low computational expense. We suggest that the present computational framework, in spite of its limitations, provides a useful foundation for future studies which may yield new insight into how aneurysms grow and rupture.     

Method for modelling cerebral blood vessels and their bifurcations using circular,homogeneous, generalised cylinders

A method for automatic modelling of blood vessels and their bifurcations from 3D scans of the brain is presented. The method is three-step procedure. First, a skeleton of the cerebral blood vessels is developed, and then the surfaces of the blood vessels are located using an active contour approach. The active contour approach uses circular homogeneous generalised cylinders (CHGCs) to model the thin, elongated blood vessels. Finally, a novel method for modelling the surfaces of the bifurcations in a vessel tree is presented. The method was tested on simulated data: a computed tomography angiography (CTA) and four magnetic resonance angiography (MRA) volumes. Furthermore, the method was tested on ten magnetic resonance images (MRIs) to demonstrate its robustness. The test on the simulated data indicated that the approach for the surface modelling of vessels had a mean radius error of less than 0.1 mm and a mean localisation error of 0.1 mm. Surface models evaluated by an expert in vascular neurosurgery were found to have a smooth appearance and generally agreed with the image data. The test on the MRI scans indicated that the method performed well in noisy environments.     

Scholarly review of geometry and compliance: biomechanical perspectives on vascular injury and healing

Mechanical stress and strain upon cardiovascular tissue are important factors that influence the ultimate configuration of clinically observed disease entities. Although mechanical forces can stimulate cellular changes and response, structural or geometric alterations introduced by disease processes can, in turn, influence local hemodynamic conditions. Dynamic interactions of structural parameters, such as arterial compliance and geometry, can further contribute to the final determination of the mechanical conditions and outcome of the vessel. Manipulation of vascular compliance and geometry may, therefore, have desirable effects. In this article, fundamental vascular biomechanical forces are defined and their association with cellular response and clinical disease processes are introduced. The interplay between vascular geometry and compliance is emphasized, and the potential for mechanical solutions to vascular diseases are explored.     

Morphometric Studies of Human Coronary Artery Trees in Healthy and Disease

Objective According to the report from American Heart Association(AHA),cardiovascular diseases(CVDs)are the leading causes of death globally,and coronary artery disease(CAD),known as coronary atherosclerotic plaques,accounts for over 30%of cardiovascular diseases.Therefore,it is of great clinical significance to study the relationship between coronary bifurcations morphometrical feature change and coronary artery disease.Although coronary atherosclerosis has been extensively investigated,there is a lack of in-deep study on the differences in morphometric features between optimal and realistic geometry of coronary arterial trees.The purpose of the present paper is to determine the morphological changes in patients with CAD lesion compared with non-coronary artery disease(non-CAD)subjects.Methods Due to the difficulty of studying the coronary bifurcations in vivo,image-based in vitro anatomical 3D models have been widely used as a noninvasive method for morphometric measurement and clinical diagnosis of the coronary bifurcations.With the development of coronary computed tomography angiography(CTA)hardware and software technologies,the CTA imaging technique has been shown a promising application in the characterization,visualization,and identification of coronary artery disease in recent decades.The CTA images used to reconstruct three-dimensional(3D)coronary arterial trees are from Asia populations(Southern Chinese populations),including five cadavers without CAD lesion and 102 patients with CAD lesion.The best fit artery diameter was calculated as twice the average radius between the points in the centerlines and the points on the coronary arterial inner wall.The bifurcation angles between larger daughter artery and smaller daughter artery were determined by the intersection angle of their centerlines.Murray's law was introduced to assess the deviation of the realistic vascular networks from its optimal state.Results Based on the morphometric analysis of coronary artery bifurcations in non-CAD subjects and patients with CAD lesion subjects,the most important finding is that morphological feature parameters of non-CAD subjects are closer to the optimal values than those of patients with CAD lesion.Moreover,by comparing the morphometric data between the left and right coronary arteries,the right coronary artery exhibits a structure closer to the optimal one in morphological feature than the left coronary artery.In addition,coronary arterial trees with CAD lesion have higher asymmetry and larger area expansion ratio(AER)than those of the coronary arterial trees without CAD lesion.Conclusions We morphologically found that the coronary arterial trees with CAD lesion and left are more likely to deviate from the optimal structure predicted by Murray's law than those without CAD lesion and right.The degree to which coronary arterial system deviating from their optimal state may directly affect the incidence of coronary artery disease.This computer morpho-logical analysis strategy is illustrated to be effective in the distinguishing of the geometric differences between the healthy and diseased coronary arteries,and the analysis method may have a large potential in cardiovascular disease earlier diagnosis.     

颅内动脉瘤径颈比对瘤体与分支血管血流动力学影响的数值模拟研究

动脉瘤的血流动力学是影响其生长与破裂的重要因素,尤其是形态学参数径颈比(aspect ratio,AR,瘤体长径/瘤颈宽度)对其血流动力学影响较大。本研究使用基于计算流体力学(computational fluid dynamics,CFD)技术的ANSYS 16.0软件包,数值仿真分析了不同径颈比对颅内动脉瘤瘤体与分支血管血流动力学的影响,为临床上制定合理的形态学与血流动力学指标来筛选高危的动脉瘤患者,并进行积极的干预治疗提供一定的理论依据。通过使用空间直角坐标系建立径颈比为3.33、2.5、2、1.67、1.43、1.25的理想颅内动脉瘤几何模型,分析和比较了包括血液流场与涡量分布、流速与流量、壁面压力、壁面切应力(wall shear stress,WSS)、瘤颈近远侧端与分支血管剪切应变率(shear strain rate,SSR)在内的血流动力学参数。数值模拟结果给出了动脉瘤与分支血管内的流线图、涡核图、压力分布云图、WSS分布云图以及随X轴变化的流速与压力峰值分布曲线。分析得出,径颈比决定瘤内血流模式,径颈比减小,瘤顶的流速与SSR增大,瘤壁上的压力与WSS增大,分支血管壁上的压力增大,且WSS/SSR瘤颈远侧端>WSS/SSR瘤颈近侧端>WSS/SSR分支血管中心,涡核区域由瘤体远侧壁增大至覆盖整个动脉瘤,但对分支血管内血流的阻碍作用减小。     

管道铸型标本断层摄影三维重建盆腔动脉

背景：以往的影像诊断中,主要观察二维图像,随着医学科技的发展及诊疗设备的不断更新,在三维空间里研究人体内小血管的分布,已成为今后解剖学研究的必然。 目的：利用聚乙烯醇-氧化铈血管造影对盆腔动脉进行三维重建,为血管介入提供可靠的解剖学依据。 方法：纳入非妇科疾病死亡成人新鲜标本2例,预处理标本,经腹主动脉注入聚乙烯醇-氧化铈对比剂,用美国QE Lights Speed VCT 64 层CT 连续扫描,利用Mimics软件进行三维重建。观察盆腔动脉的各级分支的出现以及管壁清晰度、饱满度。 结果与结论：横断面原始图像及其三维图像重建均达到满意的效果,可清晰显示4级以分支,重建后的模型可以任意三维旋转。说明聚乙烯醇-氧化铈血管造影术是盆腔动脉进行三维重建的一种有效技术方法,对妇产科血管介入技术起指导性作用。     

Hepatic tumor enhancement in computed tomography: combined models of liver perfusion and dynamic imaging

The objective of this study is to show how computational modeling can be used to increase our understanding of liver enhancement in dynamic computer tomography. It relies on two models: (1). a vascular model, based on physiological rules, is used to generate the 3D hepatic vascular network; (2). the physical process of CT acquisition allows to synthesize timed-stamped series of images, aimed at tracking the propagation of a contrast material through the vessel network and the parenchyma. The coupled models are used to simulate the enhancement of a hyper-vascular tumor at different acquisition times, showing a maximum conspicuity during the arterial phase.     

基于流固耦合计算机流体力学模拟分析人体主动脉弓内血液流动

目的比较分析应用弹性血管壁的流固耦合计算流体力学（CFD）方法和刚性血管壁的CFD方法模拟获得的正常主动脉弓内血流动力学参数,同时比较两种方法的优劣,为深入研究血液流动状态与动脉疾病的关系提供帮助。方法取46岁男性,胸主动脉正常CT图像,格式为Dicom,层间距为0.5mm,每片图像的平面分辨率为512×512,像素大小为0.5mm。应用医学图像后处理软件,对通过临床获得正常人体主动脉CT二维医学图像数据进行重构,得到主动脉血流及血管壁的三维立体模型并应用于模拟计算。结果在设定边界条件和初始条件的基础上,经多次迭代耦合计算,获得血管壁形变、等效应力、血流速度、壁面振荡切应力等相关血流动力学参数。结论在心动周期内弹性血管壁的主动脉内血流情况较刚性血管壁主动脉内血流情况更为复杂,管壁等效压力变化较大,血管壁的振荡切应力更高,表明弹性血管壁的流固耦合的CFD模拟更能体现真实主动脉内复杂血流情况,为深入研究血流动力学与心脑血管疾病的关系提供了一定的技术支持。