Growth,Anisotropy, and Residual Stresses in Arteries

A simple phenomenological theory of tissue growth is used in order to demonstrate that volumetric growth combined with material anisotropy can lead to accumulation of residual stresses in arteries. The theory is applied to growth of a cylindrical blood vessel with the anisotropy moduli derived from experiments. It is shown that bending resultants are developed in the ring cross-section of the artery. These resultants may cause the ring opening or closing after cutting the artery \textit {in vitro} as it is observed in experiments. It is emphasized that the mode of the arterial ring opening is affected by the parameters of anisotropy.     

Alternate method for the analysis of residual strain in the arterial wall

If an artery is cut transversely into rings, and the rings are then cut radially, they spring open into sectors. This phenomenon implies the existence of residual stresses and strains in the arterial wall in the non-loaded state. In the present paper, we propose a new method to calculate the residual strain from the measured wall dimensions and a polar angle of a specimen in the stress-free state, assuming that the wall is homogeneous and incompressible, and that a radially cut, stress-free specimen forms a circular sector. For this analysis, edge angles were measured at the edges of the opened-up specimen. Residual strains were obtained for the descending thoracic aorta, the common carotid artery, and the femoral artery in the rabbit. The results obtained indicated that the magnitude of residual strain was largest in the femoral artery and smallest in the aorta among the three arteries. The opening angle did not depend upon the length of a ring specimen if the ratio of the length to the diameter was ≤ 3.     

Residual strains in conduit arteries

Residual strains and stresses are those that exist in a body when all external loads are removed. Residual strains in arteries can be characterized by the opening angle of the sector-like cross-section which arises when an unloaded ring segment is radially cut. A review of experimental methods for measuring residual strains and the main results about the variation of the opening angle with arterial localization, age, smooth muscle activity, mechanical environment and certain vascular pathologies are presented and discussed. It is shown that, in addition to their well-established ability to homogenize the stress field in the arterial wall, residual strains make arteries more compliant and thereby improve their performance as elastic reservoirs and ensure more effective local control of the arterial lumen by smooth muscle cells. Finally, evidence that, in some cases, residual strains remain in arteries even after they have been cut radially is discussed.     

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%.     

A three-dimensional finite element model for arterial clamping

Clamp induced injuries of the arterial wall may determine the outcome of surgical procedures. Thus, it is important to investigate the underlying mechanical effects. We present a three-dimensional finite element model, which allows the study of the mechanical response of an artery-treated as a two-layer tube-during arterial clamping. The important residual stresses, which are associated with the load-free configuration of the artery, are also considered. In particular, the finite element analysis of the deformation process of a clamped artery and the associated stress distribution is presented. Within the clamping area a zone of axial tensile peak-stresses was identified, which (may) cause intimal and medial injury. This is an additional injury mechanism, which clearly differs from the commonly assumed wall damage occurring due to compression between the jaws of the clamp. The proposed numerical model provides essential insights into the mechanics of the clamping procedure and the associated injury mechanisms. It allows detailed parameter studies on a virtual clamped artery, which can not be performed with other methodologies. This approach has the potential to identify the most appropriate clamps for certain types of arteries and to guide optimal clamp design.     

Computational comparison of aortic root stresses in presence of stentless and stented aortic valve bio-prostheses

We provide a computational comparison of the performance of stentless and stented aortic prostheses, in terms of aortic root displacements and internal stresses. To this aim, we consider three real patients; for each of them, we draw the two prostheses configurations, which are characterized by different mechanical properties and we also consider the native configuration. For each of these scenarios, we solve the fluid–structure interaction problem arising between blood and aortic root, through Finite Elements. In particular, the Arbitrary Lagrangian–Eulerian formulation is used for the numerical solution of the fluid-dynamic equations and a hyperelastic material model is adopted to predict the mechanical response of the aortic wall and the two prostheses. The computational results are analyzed in terms of aortic flow, internal wall stresses and aortic wall/prosthesis displacements; a quantitative comparison of the mechanical behavior of the three scenarios is reported. The numerical results highlight a good agreement between stentless and native displacements and internal wall stresses, whereas higher/non-physiological stresses are found for the stented case.     

Mechanical analysis of the ring opening test applied to human ascending aortas

This work presents experiments, modelling and numerical simulation aimed at describing the mechanical response of human ascending aortas in the ring opening test. The objective is to quantify, from the opening angles measured in the test, the residual stress distribution along the artery wall and, afterwards, how this stress pattern changes when the artery is subjected to standard physiological pressures. The cases studied correspond to four groups including both healthy and pathological arteries. The tissues are characterized via tensile test measurements that enable to derive the material parameters of two constitutive models adopted in the present analysis. Overall, the numerical results obtained for all groups were found to be a useful data that allow to estimate the residual stress and their influence on the vessels under normal and hypertension physiological conditions.     

Passive material properties of intact ventricular myocardium determined from a cylindrical model

The equatorial region of the canine left ventricle was modeled as a thick-walled cylinder consisting of an incompressible hyperelastic material with homogeneous exponential properties. The anisotropic properties of the passive myocardium were assumed to be locally transversely isotropic with respect to a fiber axis whose orientation varied linearly across the wall. Simultaneous inflation, extension, and torsion were applied to the cylinder to produce epicardial strains that were measured previously in the potassium-arrested dog heart. Residual stress in the unloaded state was included by considering the stress-free configuration to be a warped cylindrical arc. In the special case of isotropic material properties, torsion and residual stress both significantly reduced the high circumferential stress peaks predicted at the endocardium by previous models. However, a resultant axial force and moment were necessary to cause the observed epicardial deformations. Therefore, the anisotropic material parameters were found that minimized these resultants and allowed the prescribed displacements to occur subject to the known ventricular pressure loads. The global minimum solution of this parameter optimization problem indicated that the stiffness of passive myocardium (defined for a 20 percent equibiaxial extension) would be 2.4 to 6.6 times greater in the fiber direction than in the transverse plane for a broad range of assumed fiber angle distributions and residual stresses. This agrees with the results of biaxial tissue testing. The predicted transmural distributions of fiber stress were relatively flat with slight peaks in the subepicardium, and the fiber strain profiles agreed closely with experimentally observed sarcomere length distributions. The results indicate that torsion, residual stress and material anisotropy associated with the fiber architecture all can act to reduce endocardial stress gradients in the passive left ventricle.     

Constitutive relation for red cell membrane. Correction.

The intention of this note is to correct a subtle and somewhat esoteric error that the author discovered in his previous publications on membrane elastic behavior. The consitutive relation between membrane force resultants and large, elastic deformations of a membrane surface involves a strain tensor, characterizing the finite deformations. The original strain tensor that appeared in the equations was the Lagrangian strain tensor; however, the proper strain representation (also Lagrangian in nature because it is "measured" relative to the undeformed material state) is transformed by rotations of coordinates in the deformed material state (whereas the Lagrangian strain tensor is transformed by rotations of coordinates in the undeformed state). The principal membrane tensions are unchanged by this correction; the material elastic constants remain the same; and therefore, the material behavior in shear and isotropic tension is the same. However, the tensor, constitutive relation can be properly applied to coordinate systems other than the principal axis system.     

A model for the nonlinear elastic response of large arteries

The nonlinear elastic response of large arteries subjected to finite deformations due to action of biaxial principal stresses, is described by simple constitutive equations. Generalized measures of strain and stress are introduced to account for material nonlinearity. This also ensures the existence of a strain energy density function. The orthotropic elastic response is described via quasi-linear relations between strains and stresses. One nonlinear parameter which defines the measures of strain and stress, and three elastic moduli are assumed to be constants. The lateral strain parameters (equivalent to Poisson's ratios in infinitesimal deformations) are deformation dependent. This dependence is defined by empirical relations developed via the incompressibility condition, and by the introduction of a fifth material parameter. The resulting constitutive model compares well with biaxial experimental data of canine carotid arteries.     

A modified Holzapfel-Ogden law for a residually stressed finite strain model of the human left ventricle in diastole

In this work, we introduce a modified Holzapfel-Ogden hyperelastic constitutive model for ventricular myocardium that accounts for residual stresses, and we investigate the effects of residual stresses in diastole using a magnetic resonance imaging–derived model of the human left ventricle (LV). We adopt an invariant-based constitutive modelling approach and treat the left ventricular myocardium as a non-homogeneous, fibre-reinforced, incompressible material. Because in vivo images provide the configuration of the LV in a loaded state even in diastole, an inverse analysis is used to determine the corresponding unloaded reference configuration. The residual stress in this unloaded state is estimated by two different methods. One is based on three-dimensional strain measurements in a local region of the canine LV, and the other uses the opening angle method for a cylindrical tube. We find that including residual stress in the model changes the stress distributions across the myocardium and that whereas both methods yield qualitatively similar changes, there are quantitative differences between the two approaches. Although the effects of residual stresses are relatively small in diastole, the model can be extended to explore the full impact of residual stress on LV mechanical behaviour for the whole cardiac cycle as more experimental data become available. In addition, although not considered here, residual stresses may also play a larger role in models that account for tissue growth and remodelling.     

Origin of axial prestretch and residual stress in arteries

The structural protein elastin endows large arteries with unique biological functionality and mechanical integrity, hence its disorganization, fragmentation, or degradation can have important consequences on the progression and treatment of vascular diseases. There is, therefore, a need in arterial mechanics to move from materially uniform, phenomenological, constitutive relations for the wall to those that account for separate contributions of the primary structural constituents: elastin, fibrillar collagens, smooth muscle, and amorphous matrix. In this paper, we employ a recently proposed constrained mixture model of the arterial wall and show that prestretched elastin contributes significantly to both the retraction of arteries that is observed upon transection and the opening angle that follows the introduction of a radial cut in an unloaded segment. We also show that the transmural distributions of elastin and collagen, compressive stiffness of collagen, and smooth muscle tone play complementary roles. Axial prestresses and residual stresses in arteries contribute to the homeostatic state of stress in vivo as well as adaptations to perturbed loads, disease, or injury. Understanding better the development of and changes in wall stress due to individual extracellular matrix constituents thus promises to provide considerable clinically important insight into arterial health and disease.     

Stretch and stress distributions in the human artery based on two-layer model considering residual stresses

Biomechanics and Modeling in Mechanobiology - The objective is to know the stress distributions in the arterial walls under residual stresses based on two-layer model. Human common carotid arteries...     

Stress-modulated growth, residual stress, and vascular heterogeneity.

A simple phenomenological model is used to study interrelations between material properties, growth-induced residual stresses, and opening angles in arteries. The artery is assumed to be a thick-walled tube composed of an orthotropic pseudoelastic material. In addition, the normal mature vessel is assumed to have uniform circumferential wall stress, which is achieved here via a mechanical growth law. Residual stresses are computed for three configurations: the unloaded intact artery, the artery after a single transmural cut, and the inner and outer rings of the artery created by combined radial and circumferential cuts. The results show that the magnitudes of the opening angles depend strongly on the heterogeneity of the material properties of the vessel wall and that multiple radial and circumferential cuts may be needed to relieve all residual stress. In addition, comparing computed opening angles with published experimental data for the bovine carotid artery suggests that the material properties change continuously across the vessel wall and that stress, not strain, correlates well with growth in arteries.     

Effect of viscoelastic relaxation on moisture transport in foods. Part I: Solution of general transport equation

Within the framework of continuum mechanics, Singh et al. [1] developed an integro-differential equation, which applies to both Darcian (Fickian) and non-Darcian (non-Fickian) modes of fluid transport in swelling biological systems. A dimensionless form of the equation was obtained and transformed from moving Eulerian to the stationary Lagrangian coordinates. Here a solution scheme for the transport equation is developed to predict moisture transport and viscoelastic stresses in spheroidal biopolymeric materials. The equation was solved numerically and results used for predicting drying and sorption curves, moisture profiles, and viscoelastic stresses in soybeans. The Lagrangian solution was obtained by assembling together several schemes: the finite element method was used to discretize the equation in space; non-linearity was addressed using the Newton-Raphson method; the Volterra term was handled via a time integration scheme of Patlashenko et al. [2] and the Galerkin rule was used to solve the time-differential term. The solution obtained in Lagrangian coordinates was transformed back to the Eulerian coordinates. In part II of this sequence we present the numerical results.Revised version: 5 October 2003     

Effect of initial stresses on the wave propagation in arteries

A theoretical analysis for the problem of wave propagation in arteries is presented. Blood is treated as a Newtonian, viscous incompressible fluid. On the basis of information derived from experimental investigations on the mechanical properties of wall tissues, the arterial wall is considered to be nonlinearly viscoelastic and orthotropic. The analysis is carried out for a cylindrical artery, under the purview of membrane theory, by taking account the effect of initial stresses. The motion of the wall and that of the fluid are assumed to be axisymmetric. Particular emphasis has been paid to the propagation of small amplitude harmonic waves whose wavelength is large compared to the radius of the vessel. By employing the equations of motion of the fluid and those for the wall, together with the equations of continuity, a frequency equation is derived by exploiting the conditions of continuity of the velocity of the arterial wall and that of blood on the endosteal surface of the wall. In order to illustrate the validity of the derived analytical expressions a quantitative analysis is made for the variations of the phase velocities as well as the transmission coefficient with frequency corresponding to different transmural pressures which relate to different initial stresses. Computational results indicate that phase velocities increase with the increase of transmural pressures.     

Compressibility of arterial wall in ring-cutting experiments

It is common practice in the arterial wall modeling to assume material incompressibility. This assumption is driven by the observation of the global volume preservation of the artery specimens in some mechanical loading experiments. The global volume preservation, however, does not necessarily imply the local volume preservation - incompressibility. In this work, we suggest to use the arterial ring- cutting experiments for the assessment of the local incompressibility assumption. The idea is to track the local stretches of the marked segments of the arterial ring after the stress-relieving cut. In the particular case of the rabbit thoracic artery, considered in this work, the following criteria for radial stretches come from preliminary analysis. If after the radial cut the marked segments shorten at the inner surface of the wall and lengthen at the outer surface while remaining unchanged in the middle of the wall then material is locally incompressible. If, however, the marked segments remain unchanged at the surfaces while lengthening in the middle of the wall then the material is locally compressible. Any other scenario would be an indication of the improper modeling assumptions, i.e. residual stresses are not relieved or material constants are inaccurate etc. It is believed that the proposed approach can be successfully implemented in experiments shedding new light on the arterial incompressibility issue.     

A critical comparison of different residence time measures in aneurysms

Flow stagnation and residence time (RT) are important features of diseased arterial flows that influence biochemical transport processes and thrombosis. RT calculation methods are classified into Eulerian and Lagrangian approaches where several measures have been proposed to quantify RT. Each of these methods has a different definition of RT, and it is not clear how they are related. In this study, image-based computational models of blood flow in an abdominal aortic aneurysm and a cerebral aneurysm were considered and RT was calculated using different methods. In the Lagrangian methods, discrete particle tracking of massless tracers was used to calculate particle residence time and mean exposure time. In the Eulerian methods, continuum transport models were used to quantify RT using Eulerian RT and virtual ink approaches. Point-wise RT and Eulerian indicator RT were also computed based on measures derived from velocity. A comparison of these methods is presented and the implications of each method are discussed. Our results highlight that most RT methods have a conceptually distinct definition of RT and therefore should be utilized depending on the specific application of interest.