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.     

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.     

The Mechanical Buckling of Curved Arteries^*

Though tortuosity and kinking are often observed in various arteries and arterioles, little is known about the underlying mechanisms. This paper presents a biomechanical analysis of bent buckling in long arterial segments with a small initial curvature using a thick-walled elastic cylindrical arterial model. The critical buckling pressure was established as a function of wall stiffness, wall dimensions, and the axial tension (or axial stretch ratio). The effects of both wall dimensions and axial stretch ratio on the critical pressure, as well as the thin-walled approximation were discussed. The buckling equation sheds light on the biomechanical mechanism of artery tortuosity and provides guidance for the development of new techniques to treat and prevent artery tortuosity and kinking.     

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.     

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.     

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.     

Propagation of dissection in a residually-stressed artery model

This paper studies dissection propagation subject to internal pressure in a residually-stressed two-layer arterial model. The artery is assumed to be infinitely long, and the resultant plane strain problem is solved using the extended finite element method. The arterial layers are modelled using the anisotropic hyperelastic Holzapfel–Gasser–Ogden model, and the tissue damage due to tear propagation is described using a linear cohesive traction–separation law. Residual stress in the arterial wall is determined by an opening angle \(\alpha \) in a stress-free configuration. An initial tear is introduced within the artery which is subject to internal pressure. Quasi-static solutions are computed to determine the critical value of the pressure, at which the dissection starts to propagate. Our model shows that the dissection tends to propagate radially outwards. Interestingly, the critical pressure is higher for both very short and very long tears. The simulations also reveal that the inner wall buckles for longer tears, which is supported by clinical CT scans. In all simulated cases, the critical pressure is found to increase with the opening angle. In other words, residual stress acts to protect the artery against tear propagation. The effect of residual stress is more prominent when a tear is of intermediate length (\(\simeq \)90\(^\circ \) arc length). There is an intricate balance between tear length, wall buckling, fibre orientation, and residual stress that determines the tear propagation.     

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.     

Arterial wall remodeling under sustained axial twisting in rats

Blood vessels often experience torsion along their axes and it is essential to understand their biological responses and wall remodeling under torsion. To this end, a rat model was developed to investigate the arterial wall remodeling under sustained axial twisting in vivo. Rat carotid arteries were twisted at 180° along the longitudinal axis through a surgical procedure and maintained for different durations up to 4 weeks. The wall remodeling in these twisted arteries was examined using histology, immunohistochemistry and fluorescent microscopy. Our data showed that arteries remodeled under twisting in a time-dependent manner during the 4 weeks post-surgery. Cell proliferation, MMP-2 and MMP-9 expressions, medial wall thickness and lumen diameter increased while collagen to elastin ratio decreased. The size and number of internal elastic lamina fenestrae increased with elongated shapes, while the endothelial cells elongated and aligned towards the blood flow direction gradually. These results demonstrated that sustained axial twisting results in artery remodeling in vivo. The rat carotid artery twisting model is an effective in vivo model for studying arterial wall remodeling under long-term torsion. These results enrich our understanding of vascular biology and arterial wall remodeling under mechanical stresses.     

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.     

Coupled computational analysis of arterial LDL transport -- effects of hypertension

Hypertension, a risk factor for atherosclerosis, increases the uptake of low density lipoproteins (LDL) by the arterial wall. Our objective in this work was to use computational modeling to identify physical factors that could be partially responsible for this effect. Fluid flow and mass transfer patterns in the lumen and wall of an arterial model were computed in a coupled manner, replicating as closely as possible previous experimental studies in which LDL uptake into the artery wall was measured in straight, excised arterial segments. Under conditions of both flow and no-flow, simulations predicted an increase in concentration polarization of LDL at the artery wall when arterial pressure was increased from 120 to 160 mmHg. However, this led to only a slight increase in mean LDL concentration within the arterial wall. However, if the permeability of the endothelium to LDL was allowed to vary with intra-arterial pressure, then the simulations predicted that the uptake of LDL would be enhanced 1.9-2.6 fold at higher pressure. The magnitude of this increase was consistent with experimental data. We conclude that the concentration polarization effects, enhanced by elevated intra-arterial pressure, cannot explain the increase in LDL uptake seen under hypertensive conditions. Instead, the data are most consistent with a pressure-linked increase in endothelial permeability to LDL.     

Torsion of the central pair microtubules in eukaryotic flagella due to bending-driven lateral buckling

Inspired by recent interest in torsion of the central pair microtubules in eukaryotic flagella, a novel thin-walled elastic beam model is suggested to study critical condition under which uniform bending of a flagellum will cause lateral/torsional buckling of the central pair. The model is directed to the central pair itself and the role of all surrounding cross-linkings inside the flagellum is modeled as an equivalent surrounding elastic medium. The model predicts that bending-driven torsion of the central pair does occur when the radius of curvature of the bent flagellum reduces to a moderate critical value typically of tens of microns. In particular, this critical value is almost independent of the flagellum length, and more sensitive to the parameters defining the surrounding elastic medium than the shear modulus of microtubules. The predicted wavelengths of the torsional buckling mode are insensitive to the flagellum length and comparable to some known related experimental data. These results indicate that torsion of the central pair microtubules in flagella is inevitable as a result of bending-driven lateral buckling. This offers an entirely new insight into the ongoing research on the mechanism of the central pair torsion.     

Pulmonary hemodynamics and gas exchange properties during progressive edema

In this investigation we have studied the effect of increments of pulmonary edema on pulmonary hemodynamics, and physiological and hemodynamic shunt in an isolated lung preparation. Hemodynamic shunt was defined by the slope of the relationship between pulmonary arterial and airway pressures; when the slope decreases, there is a greater degree of shunt. Cardiovascular changes were analyzed using a Starling resistor model of the pulmonary circulation where the effective downstream pressure to flow as seen from the pulmonary artery exceeds the pulmonary venous outflow pressure. This effective downstream pressure is referred to as the critical pressure (Pc), and at low lung inflation the locus of this critical pressure is in extra-alveolar vessels. With 3-4 h of progressive edema to an average of 185% initial lobe weight we found a progressive rise in pulmonary arterial pressure (Ppa) from 12.1 to 21.5 cmH2O. About one-third of this increase in Ppa resulted from an increased Pc and the remainder resulted from an increased resistance upstream from the locus of Pc. These results are consistent with the hypothesis that the interstitial accumulation of fluid creates enough of an increase in interstitial pressure to compress extra-alveolar vessels. There was no significant correlation between the amount of edema and the measured physiologic shunt, but the hemodynamic shunt showed a highly significant correlation. The hemodynamic shunt theoretically measures the extent of obstructed airways and may be a useful index of the degree of pulmonary edema.     

Evaluation of blood flow velocity waveform in common carotid artery using multi-branched arterial segment model of human arteries

Arteriosclerosis is considered to be a major cause of cardiovascular diseases, which account for approximately 30% of the causes of death in the world. We have recently demonstrated a strong correlation between arteriosclerosis (arterial elasticity) and two characteristics: maximum systolic velocity (S1) and systolic second peak velocity (S2) of the common carotid artery flow velocity waveform (CCFVW). The CCFVW can be measured by using a small portable measuring device. However, there is currently no theoretical evidence supporting the causes of the relation between CCFVW and arterial elasticity, or the origin of the CCFVW characteristics. In this study, the arterial blood flow was simulated using a one-dimensional systemic arterial segments model of human artery in order to conduct a qualitative evaluation of the relationship between arterial elasticity and the characteristics of CCFVW. The simulation was carried out based on the discretized segments with the physical properties of a viscoelastic tube (the cross-sectional area at the proximal and terminal ends, the length, and the compliance per unit area of the tube (C_S)). The findings obtained through this study revealed that the simulated CCFVW had shape similar characteristics to that of the measured CCFVW. Moreover, when the compliance C_S of the model was decreased, the first peak of the simulated-CCFVW decreased and the second peak increased. Further, by separating the anterograde pulse wave and the reflected pulse wave, which form the CCFVW, we found that the decrease in the first peak of the simulated CCFVW was due to the arrival of a reflected pulse wave from the head after the common carotid artery toward the arrival of a anterograde pulse wave ejected directly from the heart and that the increase in the second peak resulted from the arrival of the peak of the reflected pulse wave from the thoracic aorta. These results establish that the CCFVW characteristics contribute to the assessment of arterial elasticity.     

A model for saccular cerebral aneurysm growth by collagen fibre remodelling

The first structural model for saccular cerebral aneurysm growth is proposed. It is assumed that the development of the aneurysm is accompanied by a loss of the media, and that only collagen fibres provide load-bearing capacity to the aneurysm wall. The aneurysm is modelled as an axisymmetric multi-layered membrane, exposed to an inflation pressure. Each layer is characterized by an orientation angle, which changes between different layers. The collagen fibres and fibroblasts within a specific layer are perfectly aligned. The growth and the morphological changes of the aneurysm are accomplished by the turnover of collagen. Fibroblasts are responsible for collagen production, and the related deformations are assumed to govern the collagen production rate. There are four key parameters in the model: a normalized pressure, the number of layers in the wall, an exponent in the collagen mass production rate law, and the pre-stretch under which the collagen is deposited. The influence of the model parameters on the aneurysmal response is investigated, and a stability analysis is performed. The model is able to predict clinical observations and mechanical test results, for example, in terms of predicted aneurysm size, shape, wall stress and wall thickness.     

Influence of the distensibility of large arteries on the longitudinal impedance: application for the development of non-invasive techniques to the diagnosis of arterial diseases

Background

This study shows that the arterial longitudinal impedance constitutes a hemodynamic parameter of interest for performance characterization of large arteries in normal condition as well as in pathological situations. For this purpose, we solved the Navier?CStokes equations for an incompressible flow using the finite element analysis method and the Arbitrary Lagrangian Eulerian (ALE) formulation. The mathematical model assumes a two-dimensional flow and takes into account the nonlinear terms in the equations of fluid motion that express the convective acceleration, as well as the nonlinear deformation of the arterial wall. Several numerical simulations of the blood flow in large vessels have been performed to study the propagation along an arterial vessel of a pressure gradient pulse and a rate flow pulse. These simulations include various deformations of the wall artery leading to parietal displacements ranging from 0 (rigid wall) to 15% (very elastic wall) in order to consider physiological and pathological cases.

Results

The results show significant changes of the rate flow and the pressure gradient wave as a function of aosc, the relative variation in the radius of the artery over a cardiac cycle. These changes are notable beyond a critical value of aosc equal to 0.05. This critical value is also found in the evolution of the longitudinal impedance. So, above a variation of radius of 5%, the convective acceleration, created by the fluid-wall interactions, have an influence on the flow detectable on the longitudinal impedance.

Conclusions

The interpretation of the evolution of the longitudinal impedance shows that it could be a mean to test the performance of large arteries and can contribute to the diagnosis of parietal lesions of large arteries. For a blood vessel with a wall displacement higher than 5% similar to those of large arteries like the aorta, the longitudinal impedance is substantially greater than that obtained in the absence of wall displacement. This study also explains the effects of convective acceleration, on the shape of the decline of the pressure gradient wave and shows that they should not be neglected when the variation in radius is greater than 5%.     

Structural changes in rat aortic intima due to transmural pressure.

Huang et al. (1997) propose a new hypothesis and develop a mathematical model to explain rationally the in vitro and in situ measured changes (Tedgui and Lever, 1984; Baldwin and Wilson, 1993) in the hydraulic conductivity of the artery wall of rabbit aorta with transmural pressure. The model leads to the intriguing prediction that this hydraulic conductivity would decrease by one half if the thin intimal layer between the endothelium and the internal elastic lamina volume-compresses approximately fivefold. This paper presents the first measurements of the effect of transmural pressure on intimal layer thickness and shows that the intimal matrix is, indeed, surprisingly compressible. We perfusion-fixed rat thoracic aortas in situ with 2 percent glutaraldehyde solution at 0, 50, 100, or 150 mm Hg lumen pressure and sectioned for light and electron microscopic observations. Electron micrographs show a dramatic, nonlinear decrease in average intimal thickness, i.e., 0.62 +/- 0.26, 0.27 +/- 0.14, 0.15 +/- 0.10, and 0.12 +/- 0.07 (SD) micron for 0, 50, 100, and 150 mm Hg lumen pressure, respectively. The volume strain of the intima is more than 20 times greater than the radial strain of the artery wall due to hoop tension and two orders of magnitude greater than the consolidation of the artery wall as a whole assuming constant medial density (Chuong and Fung, 1984). Moreover, in both light and electron microscopic observations, it is easy to find numerous sites where the endothelium puckers into the fenestral pores at high lumen pressure, as predicted by the theory in Huang et al. (1997). In contrast, the average diameter of a fenestral pore increases only 10 percent as the lumen pressure is increased from 0 to 150 mm Hg. These results indicate that the thin intimal layer comprising less than 1 percent of the wall thickness can have a profound effect on the filtration properties of the wall due to the large change in Darcy permeability of the layer and the large reduction in the entrance area of the flow entering the fenestral pores, though the pores themselves experience only a minor enlargement due to hoop tension.     

Evolving biaxial mechanical properties of mouse carotid arteries in hypertension

Quantifying the time course of load-induced changes in arterial wall geometry, microstructure, and properties is fundamental to developing mathematical models of growth and remodeling. Arteries adapt to altered pressure and flow by modifying wall thickness, inner diameter, and axial length via marked cell and matrix turnover. To estimate particular biomaterial implications of such adaptations, we used a 4-fiber family constitutive relation to quantify passive biaxial mechanical behaviors of mouse carotid arteries 0 (control), 7-10, 10-14, or 35-56 days after an aortic arch banding surgery that increased pulse pressure and pulsatile flow in the right carotid artery. In vivo circumferential and axial stretches at mean arterial pressure were, for example, 11% and 26% lower, respectively, in hypertensive carotids 35-56 days after banding than in normotensive controls; this finding is consistent with observations that hypertension decreases distensibility. Interestingly, the strain energy W stored in the carotids at individual in vivo conditions was also less in hypertensive compared with normotensive carotids. For example, at 35-56 days after banding, W was 24%, 39%, and 47% of normal values at diastolic, mean, and systolic pressures, respectively. The energy stored during the cardiac cycle, W(sys)-W(dias), also tended to be less, but this reduction did not reach significance. When computed at normal in vivo values of biaxial stretch, however, W was well above normal for the hypertensive carotids. This net increase resulted from an overall increase in the collagen-related anisotropic contribution to W despite a decrease in the elastin-related isotropic contribution. The latter was consistent with observed decreases in the mass fraction of elastin.     

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.