A biomechanical model of artery buckling

The stability of arteries under blood pressure load is essential to the maintenance of normal arterial function and the loss of stability can lead to tortuosity and kinking that are associated with significant clinical complications. However, mechanical analysis of arterial bent buckling is lacking. To address this issue, this paper presents a biomechanical model of arterial buckling. Using an elastic cylindrical arterial model, the mechanical equations for arterial buckling were developed and the critical buckling pressure was found to be a function of the wall stiffness (Young's modulus), arterial radius, length, wall thickness, and the axial strain. Both the model equations and experimental results demonstrated that the critical pressure is related to the axial strain. Arteries may buckle and become tortuous due to reduced (subphysiological) axial strain, hypertensive pressure, and a weakened wall. These results are in accordance with, and provide a possible explanation to the clinical observations that hypertension and aging are the risk factors for arterial tortuosity and kinking. The current model is also applicable to veins and ureters.     

The Effect of Collagenase on the Critical Buckling Pressure of Arteries^*

The stability of arteries is essential to normal arterial functions and loss of stability can lead to arterial tortuosity and kinking. Collagen is a main extracellular matrix component that modulates the mechanical properties of arteries and collagen degradation at pathological conditions weakens the mechanical strength of arteries. However, the effects of collagen degradation on the mechanical stability of arteries are unclear. The objective of this study was to investigate the effects of collagen degradation on the critical buckling pressure of arteries. Arterial specimens were subjected to pressurized inflation testing and fitted with nonlinear thick-walled cylindrical model equations to determine their stress strain relationships. The arteries were then tested for the critical buckling pressure at a set of axial stretch ratios. Then, arteries were divided into three groups and treated with Type III collagenase at three different concentrations (64, 128, and 400U/ml). Mechanical properties and buckling pressures of the arteries were determined after collagenase treatment. Additionally, the theoretical buckling pressures were also determined using a buckling equation. Our results demonstrated that the buckling pressure of arteries was lower after collagenase treatment. The difference between pre- and post- treatment was statistically significant for the highest concentration of 400U/ml but not at the lower concentrations. The buckling equation was found to yield a fair estimation to the experimental critical pressure measurements. These results shed light on the role of matrix remodeling on the mechanical stability of arteries and developments of tortuous arteries.     

Nonlinear buckling of blood vessels: a theoretical study

Tortuosity and kinking often occur in arteries and veins but the underlying mechanisms are poorly understood. It has been suggested recently that long arteries may buckle and become tortuosity due to reduced axial tension or hypertensive pressure, but very few studies have been done to establish the biomechanical basis for artery buckling. Here we developed the arterial buckling equation using a nonlinear elastic thick-walled cylindrical model with residual stress. Our results demonstrated that arteries may buckle due to high blood pressure or low axial tension and that residual stress in the arteries increases the buckling pressure. These results are in general agreement with the previous linear elastic model. The buckling equation provides a useful tool for studying artery tortuosity and kinking.     

Artery buckling analysis using a four-fiber wall model

Artery bent buckling has been suggested as a possible mechanism that leads to artery tortuosity, which is associated with aging, hypertension, atherosclerosis, and other pathological conditions. It is necessary to understand the relationship between microscopic wall structural changes and macroscopic artery buckling behavior. To this end, the objectives of this study were to develop arterial buckling equations using a microstructure-based 4-fiber reinforced wall model, and to simulate the effects of vessel wall microstructural changes on artery buckling. Our results showed that the critical pressure increased nonlinearly with the axial stretch ratio, and the 4-fiber model predicted higher critical buckling pressures than what the Fung model predicted. The buckling equation using the 4-fiber model captured the experimentally observed reduction of critical pressure induced by elastin degradation and collagen fiber orientation changes in the arterial wall. These results improve our understanding of arterial stability and its relationship to microscopic wall remodeling, and the model provides a useful tool for further studies.     

The effect of collagenase on the critical buckling pressure of arteries

The stability of arteries is essential to normal arterial functions and loss of stability can lead to arterial tortuosity and kinking. Collagen is a main extracellular matrix component that modulates the mechanical properties of arteries and collagen degradation at pathological conditions weakens the mechanical strength of arteries. However, the effects of collagen degradation on the mechanical stability of arteries are unclear. The objective of this study was to investigate the effects of collagen degradation on the critical buckling pressure of arteries. Arterial specimens were subjected to pressurized inflation testing and fitted with nonlinear thick-walled cylindrical model equations to determine their stress strain relationships. The arteries were then tested for the critical buckling pressure at a set of axial stretch ratios. Then, arteries were divided into three groups and treated with Type III collagenase at three different concentrations (64, 128, and 400 U/ml). Mechanical properties and buckling pressures of the arteries were determined after collagenase treatment. Additionally, the theoretical buckling pressures were also determined using a buckling equation. Our results demonstrated that the buckling pressure of arteries was lower after collagenase treatment. The difference between pre- and post- treatment was statistically significant for the highest concentration of 400U/ml but not at the lower concentrations. The buckling equation was found to yield a fair estimation to the experimental critical pressure measurements. These results shed light on the role of matrix remodeling on the mechanical stability of arteries and developments of tortuous arteries.     

Effects of elastin degradation and surrounding matrix support on artery stability

Tortuous arteries are often associated with aging, hypertension, atherosclerosis, and degenerative vascular diseases, but the mechanisms are poorly understood. Our recent theoretical analysis suggested that mechanical instability (buckling) may lead to tortuous blood vessels. The objectives of this study were to determine the critical pressure of artery buckling and the effects of elastin degradation and surrounding matrix support on the mechanical stability of arteries. The mechanical properties and critical buckling pressures, at which arteries become unstable and deform into tortuous shapes, were determined for a group of five normal arteries using pressurized inflation and buckling tests. Another group of nine porcine arteries were treated with elastase (8 U/ml), and the mechanical stiffness and critical pressure were obtained before and after treatment. The effect of surrounding tissue support was simulated using a gelatin gel. The critical pressures of the five normal arteries were 9.52 kPa (SD 1.53) and 17.10 kPa (SD 5.11) at axial stretch ratios of 1.3 and 1.5, respectively, while model predicted critical pressures were 10.11 kPa (SD 3.12) and 17.86 kPa (SD 5.21), respectively. Elastase treatment significantly reduced the critical buckling pressure (P < 0.01). Arteries with surrounding matrix support buckled into multiple waves at a higher critical pressure. We concluded that artery buckling under luminal pressure can be predicted by a buckling equation. Elastin degradation weakens the arterial wall and reduces the critical pressure, which thus leads to tortuous vessels. These results shed light on the mechanisms of the development of tortuous vessels due to elastin deficiency.     

Fluid-structure interaction modeling of aneurysmal arteries under steady-state and pulsatile blood flow: a stability analysis

Tortuous aneurysmal arteries are often associated with a higher risk of rupture but the mechanism remains unclear. The goal of this study was to analyze the buckling and post-buckling behaviors of aneurysmal arteries under pulsatile flow. To accomplish this goal, we analyzed the buckling behavior of model carotid and abdominal aorta with aneurysms by utilizing fluid-structure interaction (FSI) method with realistic waveforms boundary conditions. FSI simulations were done under steady-state and pulsatile flow for normal (1.5) and reduced (1.3) axial stretch ratios to investigate the influence of aneurysm, pulsatile lumen pressure and axial tension on stability. Our results indicated that aneurysmal artery buckled at the critical buckling pressure and its deflection nonlinearly increased with increasing lumen pressure. Buckling elevates the peak stress (up to 118%). The maximum aneurysm wall stress at pulsatile FSI flow was (29%) higher than under static pressure at the peak lumen pressure of 130 mmHg. Buckling results show an increase in lumen shear stress at the inner side of the maximum deflection. Vortex flow was dramatically enlarged with increasing lumen pressure and artery diameter. Aneurysmal arteries are more susceptible than normal arteries to mechanical instability which causes high stresses in the aneurysm wall that could lead to aneurysm rupture.     

On the Onset of Cracks in Arteries¹

We present a theoretical approach to study the onset of failure localization into cracks in arterial wall. The arterial wall is a soft composite comprising hydrated ground matrix of proteoglycans reinforced by spatially dispersed elastin and collagen fibers. As any material, the arterial tissue cannot accumulate and dissipate strain energy beyond a critical value. This critical value is enforced in the constitutive theory via energy limiters. The limiters automatically bound reachable stresses and allow examining the mathematical condition of strong ellipticity. Loss of the strong ellipticity physically means inability of material to propagate superimposed waves. The waves cannot propagate because material failure localizes into cracks perpendicular to a possible wave direction. Thus, not only the onset of a crack can be analyzed but also its direction. We use the recently developed constitutive theories of the arterial wall including 8 and 16 structure tensors to account for the fiber dispersion. We enhance these theories with energy limiters. We examine the loss of strong ellipticity in uniaxial tension and pure shear in circumferential and axial directions of the arterial wall. We find that the vanishing longitudinal wave speed predicts the appearance of cracks in the direction perpendicular to tension. We also find that the vanishing transverse wave speed predicts the appearance of cracks in the the direction inclined (non-perpendicular) to tension. The latter result is counter-intuitive yet it is supported by recent experimental observations.     

Mechanical instability of normal and aneurysmal arteries

Tortuous arteries associated with aneurysms have been observed in aged patients with atherosclerosis and hypertension. However, the underlying mechanism is poorly understood. The objective of this study was to determine the effect of aneurysms on arterial buckling instability and the effect of buckling on aneurysm wall stress. We investigated the mechanical buckling and post-buckling behavior of normal and aneurysmal carotid arteries and aorta’s using computational simulations and experimental measurements to elucidate the interrelationship between artery buckling and aneurysms. Buckling tests were done in porcine carotid arteries with small aneurysms created using elastase treatment. Parametric studies were done for model aneurysms with orthotropic nonlinear elastic walls using finite element simulations. Our results demonstrated that arteries buckled at a critical buckling pressure and the post-buckling deflection increased nonlinearly with increasing pressure. The presence of an aneurysm can reduce the critical buckling pressure of arteries, although the effect depends on the aneurysm’s dimensions. Buckled aneurysms demonstrated a higher peak wall stress compared to unbuckled aneurysms under the same lumen pressure. We conclude that aneurysmal arteries are vulnerable to mechanical buckling and mechanical buckling could lead to high stresses in the aneurysm wall. Buckling could be a possible mechanism for the development of tortuous aneurysmal arteries such as in the Loeys–Dietz syndrome.     

Mechanical buckling of artery under pulsatile pressure

Tortuosity that often occurs in carotid and other arteries has been shown to be associated with high blood pressure, atherosclerosis, and other diseases. However the mechanisms of tortuosity development are not clear. Our previous studies have suggested that arteries buckling could be a possible mechanism for the initiation of tortuous shape but artery buckling under pulsatile flow condition has not been fully studied. The objectives of this study were to determine the artery critical buckling pressure under pulsatile pressure both experimentally and theoretically, and to elucidate the relationship of critical pressures under pulsatile flow, steady flow, and static pressure. We first tested the buckling pressures of porcine carotid arteries under these loading conditions, and then proposed a nonlinear elastic artery model to examine the buckling pressures under pulsatile pressure conditions. Experimental results showed that under pulsatile pressure arteries buckled when the peak pressures were approximately equal to the critical buckling pressures under static pressure. This was also confirmed by model simulations at low pulse frequencies. Our results provide an effective tool to predict artery buckling pressure under pulsatile pressure.     

Twist buckling of veins under torsional loading

Veins are often subjected to torsion and twisted veins can hinder and disrupt normal blood flow but their mechanical behavior under torsion is poorly understood. The objective of this study was to investigate the twist deformation and buckling behavior of veins under torsion. Twist buckling tests were performed on porcine internal jugular veins (IJVs) and human great saphenous veins (GSVs) at various axial stretch ratio and lumen pressure conditions to determine their critical buckling torques and critical buckling twist angles. The mechanical behavior under torsion was characterized using a two-fiber strain energy density function and the buckling behavior was then simulated using finite element analysis. Our results demonstrated that twist buckling occurred in all veins under excessive torque characterized by a sudden kink formation. The critical buckling torque increased significantly with increasing lumen pressure for both porcine IJV and human GSV. But lumen pressure and axial stretch had little effect on the critical twist angle. The human GSVs are stiffer than the porcine IJVs. Finite element simulations captured the buckling behavior for individual veins under simultaneous extension, inflation, and torsion with strong correlation between predicted critical buckling torques and experimental data (R² = 0.96). We conclude that veins can buckle under torsion loading and the lumen pressure significantly affects the critical buckling torque. These results improve our understanding of vein twist behavior and help identify key factors associated in the formation of twisted veins.     

Twist buckling behavior of arteries

Arteries are often subjected to torsion due to body movement and surgical procedures. While it is essential that arteries remain stable and patent under twisting loads, the stability of arteries under torsion is poorly understood. The goal of this work was to experimentally investigate the buckling behavior of arteries under torsion and to determine the critical buckling torque, the critical buckling twist angle, and the buckling shape. Porcine common carotid arteries were slowly twisted in vitro until buckling occurred while subjected to a constant axial stretch ratio (1.1, 1.3, 1.5 (in vivo level) and 1.7) and lumen pressure (20, 40, 70 and 100 mmHg). Upon buckling, the arteries snapped to form a kink. For a group of six arteries, the axial stretch ratio significantly affected the critical buckling torque (\(p<0.002\)) and the critical buckling twist angle (\(p<0.001\)). Lumen pressure also significantly affected the critical buckling torque (\(p<0.001\)) but had no significant effect on the critical twist angle (\(p=0.067\)). Convex material constants for a Fung strain energy function were determined and fit well with the axial force, lumen pressure, and torque data measured pre-buckling. The material constants are valid for axial stretch ratios, lumen pressures, and rotation angles of 1.3–1.5, 20–100 mmHg, and 0–270\(^\circ \), respectively. The current study elucidates the buckling behavior of arteries under torsion and provides new insight into mechanical instability of blood vessels.     

Digital Subtraction Angiography Imaging Characteristics of Patients with Extra–Intracranial Atherosclerosis and Its Relationship to Stroke

To investigate the angiographic characteristics and clinical features in patients with suspected extra–intracranial atherosclerosis in a large cohort of Chinese population. On the basis of digital subtraction angiography characteristics, pathological morphology of extra–intracranial atherosclerosis was divided into tortuosity, kinking, coiling, and stenosis in 2,218 individuals aged 45–89 years. The degree of stenosis was further divided into low-grade (<30 %), intermediate-grade (30–69 %), and high-grade stenosis (≥70 %). Clinical manifestations were divided into transient ischemic attack, cerebral infarction and cerebral hemorrhage. The prevalence of tortuosity and stenosis were significantly higher in the extracranial arterial system than that of intracranial arterial system. The prevalence of tortuosity and kinking were significantly higher on the left side than the right side. The prevalence of mild and moderate stenosis in the internal carotid artery was significantly higher in the left side than the right side. The incidence of cerebral infarction was significantly higher in the internal carotid arterial (ICA) system than the vertebrobasilar arterial (VBA) system. Tortuosity is a common carotid abnormality in the Chinese population. The prevalence of ICA tortuosity is higher than that of VBA. The incidence of cerebral infarction in each atherosclerosis group was significantly higher in ICA than that of VBA. The prevalence of stroke is higher in the ICA system than the VBA system. Kinkings and coilings may not have a clinical significance if these lesions are not associated with atheromatous plaques or carotid stenosis.     

Diameter-dependent axial prestretch of porcine coronary arteries and veins

The pressure-diameter relation (PDR) and the wall strain of coronary blood vessels have important implications for coronary blood flow and arthrosclerosis, respectively. Previous studies have shown that these mechanical quantities are significantly affected by the axial stretch of the vessels. The objective of this study was to measure the physiological axial stretch in the coronary vasculature; i.e., from left anterior descending (LAD) artery tree to coronary sinus vein and to determine its effect on the PDR and hence wall stiffness. Silicone elastomer was perfused through the LAD artery and coronary sinus trees to cast the vessels at the physiologic pressure. The results show that the physiological axial stretch exists for orders 4 to 11 (> 24 μm in diameter) arteries and orders -4 to -12 (>38 μm in diameter) veins but vanishes for the smaller vessels. Statistically, the axial stretch is higher for larger vessels and is higher for arteries than veins. The axial stretch λ(z) shows a linear variation with the order number (n) as: λ(z) = 0.062n + 0.75 (R(2) = 0.99) for artery and λ(z) = -0.029n + 0.89 (R(2) = 0.99) for vein. The mechanical analysis shows that the axial stretch significantly affects the PDR of the larger vessels. The circumferential stretch/strain was found to be significantly higher for the epicardial arteries (orders 9-11), which are free of myocardium constraint, than the intramyocardial arteries (orders 4-8). These findings have fundamental implications for coronary blood vessel mechanics.     

A stress-strain relation for a rat abdominal aorta

Assuming the arterial wall is homogeneous, incompressible, isotropic and elastic, a stress-strain relation has been presented for a rat's abdominal aorta. As an illustrating example, the problem of simultaneous inflation and the axial stretch of a cylindrical artery under physiological loading has been solved and then the material coefficients are determined by comparing theoretical results with the existing experiments. The result indicates that the maximum deviation between the theory and experiment for various pressure levels is 3.7% which seems to be a good approximation of theory to the experiments. The variation of circumferential stress and the incremental pressure modulus with inner pressure are also depicted in the work.     

Mechanical anisotropy of inflated elastic tissue from the pig aorta

Uniaxial and biaxial mechanical properties of purified elastic tissue from the proximal thoracic aorta were studied to understand physiological load distributions within the arterial wall. Stress–strain behaviour was non-linear in uniaxial and inflation tests. Elastic tissue was 40% stiffer in the circumferential direction compared to axial in uniaxial tests and～100% stiffer in vessels at an axial stretch ratio of 1.2 or 1.3 and inflated to physiological pressure. Poisson’s ratio v_θz averaged 0.2 and v_zθ increased with circumferential stretch from ～0.2 to ～0.4. Axial stretch had little impact on circumferential behaviour. In intact (unpurified) vessels at constant length, axial forces decreased with pressure at low axial stretches but remained constant at higher stretches. Such a constant axial force is characteristic of incrementally isotropic arteries at their in vivo dimensions. In purified elastic tissue, force decreased with pressure at all axial strains, showing no trend towards isotropy. Analysis of the force–length–pressure data indicated a vessel with v_θz≈0.2 would stretch axially 2–4% with the cardiac pulse yet maintain constant axial force. We compared the ability of 4 mathematical models to predict the pressure-circumferential stretch behaviour of tethered, purified elastic tissue. Models that assumed isotropy could not predict the stretch at zero pressure. The neo-Hookean model overestimated the non-linearity of the response and two non-linear models underestimated it. A model incorporating contributions from orthogonal fibres captured the non-linearity but not the zero-pressure response. Models incorporating anisotropy and non-linearity should better predict the mechanical behaviour of elastic tissue of the proximal thoracic aorta.     

Blood vessel buckling within soft surrounding tissue generates tortuosity

The stability of blood vessel under lumen pressure load is essential to the maintenance of normal arterial function. Previous mechanical models showed that blood vessels may buckle into a half sine wave but arteries and veins in vivo often demonstrate tortuous paths with multiple waves. The objective of this study was to analyze the buckling of blood vessels under lumen pressure with surrounding tissue support. Blood vessels were modeled as elastic cylindrical vessels within an elastic substrate. Buckling equations were established to determine the critical pressure and the wavelength. These equations and simulation results demonstrated that blood vessels do take higher order mode shapes when buckling inside an elastic substrate while they take the basal mode shape without the substrate. The wave number increases i.e. blood vessels take a higher mode shape, as the stiffness of the substrate increases. These results suggest that mechanical buckling is a possible mechanism for the development of tortuous blood vessels. The current model provides a powerful tool for further studying the tortuosity of arteries and veins.     

Patient-Specific Artery Shrinkage and 3D Zero-Stress State in Multi-Component 3D FSI Models for Carotid Atherosclerotic Plaques Based on <i>In Vivo</i> MRI Data

Image-based computational models for atherosclerotic plaques have been developed to perform mechanical analysis to quantify critical flow and stress/strain conditions related to plaque rupture which often leads directly to heart attack or stroke. An important modeling issue is how to determine zero stress state from in vivo plaque geometries. This paper presents a method to quantify human carotid artery axial and inner circumferential shrinkages by using patient-specific ex vivo and in vivo MRI images. A shrink-stretch process based on patient-specific in vivo plaque morphology and shrinkage data was introduced to shrink the in vivo geometry first to find the zero-stress state (opening angle was ignored to reduce the complexity), and then stretch and pressurize to recover the in vivo plaque geometry with computed initial stress, strain, flow pressure and velocity conditions. Effects of the shrink-stretch process on plaque stress/strain distributions were demonstrated based on patient-specific data using 3D models with fluid-structure interactions (FSI). The average artery axial and inner circumferential shrinkages were 25% and 7.9%, respectively, based on a data set obtained from 10 patients. Maximum values of maximum principal stress and strain increased 349.8% and 249% respectively with 33% axial stretch. Influence of inner circumferential shrinkage (7.9%) was not very noticeable under 33% axial stretch, but became more noticeable under smaller axial stretch. Our results indicated that accurate knowledge of artery shrinkages and the shrink-stretch process will considerably improve the accuracy of computational predictions made based on results from those in vivo MRI-based FSI models.     

Pulsatile arterial wall-blood flow interaction with wall pre-stress computed using an inverse algorithm

Background

The computation of arterial wall deformation and stresses under physiologic conditions requires a coupled compliant arterial wall-blood flow interaction model. The in-vivo arterial wall motion is constrained by tethering from the surrounding tissues. This tethering, together with the average in-vivo pressure, results in wall pre-stress. For an accurate simulation of the physiologic conditions, it is important to incorporate the wall pre-stress in the computational model. The computation of wall pre-stress is complex, as the un-loaded and un-tethered arterial shape with residual stress is unknown. In this study, the arterial wall deformation and stresses in a canine femoral artery under pulsatile pressure was computed after incorporating the wall pre-stresses. A nonlinear least square optimization based inverse algorithm was developed to compute the in-vivo wall pre-stress.

Methods

First, the proposed inverse algorithm was used to obtain the un-loaded and un-tethered arterial geometry from the unstressed in-vivo geometry. Then, the un-loaded, and un-tethered arterial geometry was pre-stressed by applying a mean in-vivo pressure of 104.5 mmHg and an axial stretch of 48% from the un-tethered length. Finally, the physiologic pressure pulse was applied at the inlet and the outlet of the pre-stressed configuration to calculate the in-vivo deformation and stresses. The wall material properties were modeled with an incompressible, Mooney-Rivlin model derived from previously published experimental stress-strain data (Attinger et al., 1968).

Results

The un-loaded and un-tethered artery geometry computed by the inverse algorithm had a length, inner diameter and thickness of 35.14 mm, 3.10 mm and 0.435 mm, respectively. The pre-stressed arterial wall geometry was obtained by applying the in-vivo axial-stretch and average in-vivo pressure to the un-loaded and un-tethered geometry. The length of the pre-stressed artery, 51.99 mm, was within 0.01 mm (0.019%) of the in-vivo length of 52.0 mm; the inner diameter of 3.603 mm was within 0.003 mm (0.08%) of the corresponding in-vivo diameter of 3.6 mm, and the thickness of 0.269 mm was within 0.0015 mm (0.55%) of the in-vivo thickness of 0.27 mm. Under physiologic pulsatile pressure applied to the pre-stressed artery, the time averaged longitudinal stress was found to be 42.5% higher than the circumferential stresses. The results of this study are similar to the results reported by Zhang et al., (2005) for the left anterior descending coronary artery.

Conclusions

An inverse method was adopted to compute physiologic pre-stress in the arterial wall before conducting pulsatile hemodynamic calculations. The wall stresses were higher in magnitude in the longitudinal direction, under physiologic pressure after incorporating the effect of in-vivo axial stretch and pressure loading.