Pulse waves in prestressed arteries

In order to better understand the effect of initial stress in blood flow in arteries, a theoretical analysis of wave propagation in an initially inflated and axially stretched cylindrical thick shell is investigated. For simplicity in the mathematical analysis, the blood is assumed to be an incompressible inviscid fluid while the arterial wall is taken to be an isotropic, homogeneous and incompressible elastic material. Employing the theory of small deformations superimposed on a large initial field the governing differential equations of perturbed solid motions are obtained in cylindrical polar coordinates. Considering the difficulty in obtaining a closed form solution for the field equations, an approximate power series method is utilized. The dispersion relations for the most general case of this approximation and for the thin tube case are thoroughly discussed. The speeds of waves propagating in an unstressed tube are obtained as a special case of our general treatment. It is observed that the speeds of both waves increase with increasing inner pressure and axial stretch.     

A description of arterial wall mechanics using limiting chain extensibility constitutive models

Certain aspects of the mechanical response of arterial walls can be described using nonlinear elasticity theory. Uniaxial tests on vascular walls reveal nonlinear stress-strain behavior, with higher extensibility in the low stretch range and progressively lower extensibility with increasing stretch. This phenomenon is well known in the framework of rubber-like materials where it is called a strain-hardening or strain-stiffening effect. Constitutive models of incompressible hyperelasticity that take this into account include power-law models and limiting chain extensibility models. Our purpose in this paper is to bring to the attention of the biomechanics community some essential features of one such model of the latter type due to Gent. This model is compared with isotropic versions of biomechanical constitutive models by Takamizawa-Hayashi and Fung; the latter is a limiting version of a power-law material. Two particular problems are considered for which experimental data on arterial wall deformations are available. The first concerns small oscillations superposed on a large static stretch of a vertical string of arterial tissue. It is shown that the exponential model of Fung and the Gent model match well with the experimental data. The second problem is the extension of an internally pressurized circular cylindrical tube. It is shown that an inversion phenomenon observed experimentally for the human iliac artery can be described within a membrane theory by the Gent model whereas this cannot be described using the exponential model. The foregoing considerations are carried out for isotropic elastic materials in the absence of residual stress. Extensions to include anisotropy are also indicated.     

A new constitutive model for multi-layered collagenous tissues

Collagenous tissues such as the aneurysmal wall or the aorta are multi-layered structures with the mean fibre alignments distinguishing one layer from another. A constitutive representation of the multiple collagen layers is not yet developed, and hence the aim of the present study. The proposed model is based on the constitutive theory of finite elasticity and is characterized by an anisotropic strain-energy function which takes the material structure into account. The passive tissue behaviour is modelled and the related mechanical response is assumed to be dominated by elastin and collagen. While elastin is modelled by the neo-Hookean material the constitutive response of collagen is assumed to be transversely isotropic for each individual layer and based on an exponential function. The proposed constitutive function is polyconvex which ensures material stability. The model has five independent material parameters, each of which has a clear physical interpretation: the initial stiffnesses of the collagen fabric in the two principal directions, the shear modulus pertaining to the non-collagenous matrix material, a parameter describing the level of nonlinearity of the collagen fabric, and the angle between the principal directions of the collagen fabric and the reference coordinate system. An extension-inflation test of the adventitia of a human femoral artery is simulated by means of the finite element method and an error function is minimized by adjusting the material parameters yielding a good agreement between the model and the experimental data.     

A goal function approach to remodeling of arteries uncovers mechanisms for growth instability

A novel, goal function-based formulation for the growth dynamics of arteries is introduced and used for investigating the development of growth instability in blood vessels. Such instabilities would lead to abnormal growth of the vessel, reminiscent of an aneurysm. The blood vessel is modeled as a thin-walled cylindrical tube, and the constituents that form the vessel wall are assumed to deform together as a constrained mixture. The growth dynamics of the composite material of the vessel wall are described by an evolution equation, where the effective area of each constituent changes in the direction of steepest descent of a goal function. This goal function is formulated in such way that the constituents grow toward a target potential energy and a target composition. The convergence of the simulated response of the evolution equation toward a target homeostatic state is investigated for a range of isotropic and orthotropic material models. These simulations suggest that elastin-deficient vessels are more prone to growth instability. Increased stiffness of the vessel wall, on the other hand, gives a more stable growth process. Another important finding is that an increased rate of degradation of materials impairs growth stability.     

Role of elastin anisotropy in structural strain energy functions of arterial tissue

The vascular wall exhibits nonlinear anisotropic mechanical properties. The identification of a strain energy function (SEF) is the preferred method to describe its complex nonlinear elastic properties. Earlier constituent-based SEF models, where elastin is modeled as an isotropic material, failed in describing accurately the tissue response to inflation–extension loading. We hypothesized that these shortcomings are partly due to unaccounted anisotropic properties of elastin. We performed inflation–extension tests on common carotid of rabbits before and after enzymatic degradation of elastin and applied constituent-based SEFs, with both an isotropic and an anisotropic elastin part, on the experimental data. We used transmission electron microscopy (TEM) and serial block-face scanning electron microscopy (SBFSEM) to provide direct structural evidence of the assumed anisotropy. In intact arteries, the SEF including anisotropic elastin with one family of fibers in the circumferential direction fitted better the inflation–extension data than the isotropic SEF. This was supported by TEM and SBFSEM imaging, which showed interlamellar elastin fibers in the circumferential direction. In elastin-degraded arteries, both SEFs succeeded equally well in predicting anisotropic wall behavior. In elastase-treated arteries fitted with the anisotropic SEF for elastin, collagen engaged later than in intact arteries. We conclude that constituent-based models with an anisotropic elastin part characterize more accurately the mechanical properties of the arterial wall when compared to models with simply an isotropic elastin. Microstructural imaging based on electron microscopy techniques provided evidence for elastin anisotropy. Finally, the model suggests a later and less abrupt collagen engagement after elastase treatment.     

The pollen tube: a soft shell with a hard core

Plant cell expansion is controlled by a fine‐tuned balance between intracellular turgor pressure, cell wall loosening and cell wall biosynthesis. To understand these processes, it is important to gain in‐depth knowledge of cell wall mechanics. Pollen tubes are tip‐growing cells that provide an ideal system to study mechanical properties at the single cell level. With the available approaches it was not easy to measure important mechanical parameters of pollen tubes, such as the elasticity of the cell wall. We used a cellular force microscope (CFM) to measure the apparent stiffness of lily pollen tubes. In combination with a mechanical model based on the finite element method (FEM), this allowed us to calculate turgor pressure and cell wall elasticity, which we found to be around 0.3 MPa and 20–90 MPa, respectively. Furthermore, and in contrast to previous reports, we showed that the difference in stiffness between the pollen tube tip and the shank can be explained solely by the geometry of the pollen tube. CFM, in combination with an FEM‐based model, provides a powerful method to evaluate important mechanical parameters of single, growing cells. Our findings indicate that the cell wall of growing pollen tubes has mechanical properties similar to rubber. This suggests that a fully turgid pollen tube is a relatively stiff, yet flexible cell that can react very quickly to obstacles or attractants by adjusting the direction of growth on its way through the female transmitting tissue.     

Biorheology and fluid flux in swelling tissues, II. Analysis of unconfined compressive response of transversely isotropic cartilage disc

The response of a cartilage disc to unconfined compressive loading under small deformations is analyzed. The cartilage is considered as a transversely isotropic bicomponent (solid-fluid) tissue. Concentration effects (commonly termed osmotic pressure) are accounted for. The tissue's permeability is taken to be isotropic. Its concentration force is assumed to vary linearly with volume. The analysis shows that if the tissue's fibrous structure is taken into consideration, then the instantaneous response to a step loading depends on the tissue's elasticity and on its concentration force. The subsequent creep response, under commonly used experimental conditions, has a time constant which depends on the concentration force and permeability, but independent of its elastic response. The equilibrium volume is predicted to depend only on the concentration force. Where data is available it confirms the model's predictions. It is concluded from the present analysis that inclusion of concentration effects and the tissue's fibrous structure has significant consequences in terms of the relative roles of the collagen fibers (solid) vs. the ground substance (fluid) in the response of the cartilage to compressive loading.     

Comparative finite element modelling of aneurysm formation and physiologic inflation in the descending aorta

Abstract

Despite the general interest in aneurysm rupture prediction, the aneurysm formation has received limited attention. The goal of this study is to assess whether an aneurysm may be instigated in a healthy model of an aorta inflated by a supra-physiological pressure. The effect of two main aspects on numerical predictions has been explored: i) the geometric design and ii) the constitutive law adopted to represent the material properties. Firstly, higher values of wall stress and displacement magnitude were generated in the physiologic model compared to the cylindrical one when assigning the same material properties. Secondly, greater deformations are observed in the anisotropic model compared to the isotropic one.     

Passive elastic properties of the rat aorta

The passive anisotropic elastic properties of rat's aorta were studied in vitro by subjecting cylindrical segments of thoracic and abdominal aorta to a wide range of deformations. Using data on pressure, axial stretch, outer diameter, axial force and wall thickness, incremental moduli of elasticity in the circumferential, axial and radial directions were computed. Results indicate that while the elastic behavior of the aortic wall is globally anisotropic, there exists a state of deformation at which the vessel displays incremental isotropy. This state of deformation corresponds approximately to the loading conditions to which the aorta is exposed in situ. Values of the moduli, analyzed as a function of transmural pressure, show that the stiffness of the aortic wall is fairly constant at low pressures but raises steeply for pressures higher than physiological. For axial stretches as occurring in situ, the magnitudes of the circumferential and radial moduli do not differ significantly for the thoracic aorta; hence this vessel can be regarded as transversely isotropic over a wide range of pressures. The same observation is valid also for the abdominal aorta when pressures equal or smaller than physiological are considered. For both the thoracic and abdominal segments of the aorta, the circumferential and radial moduli are smaller than the axial modulus at low pressures, while the reverse is true for large pressures.     

Site specific inelasticity of arterial tissue

Understanding the mechanical behaviour of arterial tissue is vital to the development and analysis of medical devices targeting diseased vessels. During angioplasty and stenting, stress softening and permanent deformation of the vessel wall occur during implantation of the device, however little data exists on the inelastic behaviour of cardiovascular tissue and how this varies through the arterial tree. The aim of this study was to characterise the magnitude of stress softening and inelastic deformations due to loading throughout the arterial tree and to investigate the anisotropic inelastic behaviour of the tissue. Cyclic compression tests were used to investigate the differences in inelastic behaviour for carotid, aorta, femoral and coronary arteries harvested from 3-4 month old female pigs, while the anisotropic behaviour of aortic and carotid tissue was determined using cyclic tensile tests in the longitudinal and circumferential directions. The differences in inelastic behaviour were correlated to the ratio of collagen to elastin content of the arteries. It was found that larger inelastic deformations occurred in muscular arteries (coronary), which had a higher collagen to elastin ratio than elastic arteries (aorta), where the smallest inelastic deformations were observed. Lower magnitude inelastic deformations were observed in the circumferential tensile direction than in the longitudinal tensile direction or due to radial compression. This may be as a result of non-collagenous components in the artery becoming more easily damaged than the collagen fibres during loading. Stress softening was also found to be dependent on artery type. In the future, computational models should consider such site dependant, anisotropic inelastic behaviour in order to better predict the outcomes of interventional procedures such as angioplasty and stenting.     

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.     

Estimation of aneurysm wall stresses created by treatment with a shape memory polymer foam device

In this study, compliant latex thin-walled aneurysm models are fabricated to investigate the effects of expansion of shape memory polymer foam. A simplified cylindrical model is selected for the in-vitro aneurysm, which is a simplification of a real, saccular aneurysm. The studies are performed by crimping shape memory polymer foams, originally 6 and 8 mm in diameter, and monitoring the resulting deformation when deployed into 4-mm-diameter thin-walled latex tubes. The deformations of the latex tubes are used as inputs to physical, analytical, and computational models to estimate the circumferential stresses. Using the results of the stress analysis in the latex aneurysm model, a computational model of the human aneurysm is developed by changing the geometry and material properties. The model is then used to predict the stresses that would develop in a human aneurysm. The experimental, simulation, and analytical results suggest that shape memory polymer foams have potential of being a safe treatment for intracranial saccular aneurysms. In particular, this work suggests oversized shape memory foams may be used to better fill the entire aneurysm cavity while generating stresses below the aneurysm wall breaking stresses.     

Wave Propagation in a Viscous Fluid Contained in an Orthotropic Elastic Tube

The problem of pressure wave propagation through a viscous fluid contained in an orthotropic elastic tube is considered in connection with arterial blood flow. Solutions to the fluid flow and elasticity equations are obtained for the presence of a reflected wave. Numerical results are presented for both isotropic and orthotropic elastic tubes. In particular, the pressure pulse, flow rate, axial fluid velocity, and wall displacements are plotted vs. time at various stations along the ascending aorta of man. The results indicate an increase in the peak value of the pressure pulse and a decrease in the flow rate as the pulse propagates away from the heart. Finally, the velocity of wave propagation depends mainly on the tangential modulus of elasticity of the arterial wall, and anisotropy of the wall accounts in part for the reduction of longitudinal movements and an increase in the hydraulic resistance.     

General tube law for collapsible thin and thick-wall tubes

Modeling the complex deformations of cylindrical tubes under external pressure is of interest in engineering and physiological applications. The highly non-linear post-buckling behavior of cross-section of the tube during collapse attracted researchers for years. Major efforts were concentrated on studying the behavior of thin-wall tubes. Unfortunately, the knowledge on post-buckling of thick-wall tubes is still incomplete, although many experimental and several theoretical studies have been performed. In this study we systematically studied the effect of the wall thickness on post-buckling behavior of the tube. For this purpose, we utilized a computational model for evaluation of the real geometry of the deformed cross-sectional area due to negative transmural (internal minus external) pressure. We also developed an experimental method to validate the computational results. Based on the computed cross-sections of tubes with different wall thicknesses, we developed a general tube law that accounts for thin or thick wall tubes and fits the numerical data of computed cross-sectional areas versus transmural pressures.     

On a nonlinear theory for muscle shells: Part II--Application to the beating left ventricle

This paper specializes the nonlinear laminated-muscle-shell theory developed in Part I to cylindrical geometry and computes stresses in arteries and the beating left ventricle. The theory accounts for large strain, material nonlinearity, thick-shell effects, torsion, muscle activation, and residual strain. First, comparison with elasticity solutions for pressurized arteries shows that the accuracy of the shell theory increases as transmural stress gradients and the shell thickness decrease. Residual strain reduces the stress gradients, lowering the error in the predicted peak stress in thick-walled arteries (R/t = 2.8) from about 30 to 10 percent. Second, the canine left ventricle is modeled as a thick-walled laminated cylinder with an internal pressure. Each layer is composed of transversely isotropic muscle with a fiber orientation based on anatomical data. Using a single pseudostrain-energy density function (with time-varying coefficients) for passive and active myocardium, the model predicts strain distributions that agree fairly well with published experimental measurements. The results also show that the peak fiber stress occurs subendocardially near the beginning of ejection and that residual strains significantly alter stress gradients within each lamina, but the magnitude of the peak fiber stress changes by less than 20 percent.     

Transversely isotropic elasticity imaging of cancellous bone

To measure spatial variations in mechanical properties of biological materials, prior studies have typically performed mechanical tests on excised specimens of tissue. Less invasive measurements, however, are preferable in many applications, such as patient-specific modeling, disease diagnosis, and tracking of age- or damage-related degradation of mechanical properties. Elasticity imaging (elastography) is a nondestructive imaging method in which the distribution of elastic properties throughout a specimen can be reconstructed from measured strain or displacement fields. To date, most work in elasticity imaging has concerned incompressible, isotropic materials. This study presents an extension of elasticity imaging to three-dimensional, compressible, transversely isotropic materials. The formulation and solution of an inverse problem for an anisotropic tissue subjected to a combination of quasi-static loads is described, and an optimization and regularization strategy that indirectly obtains the solution to the inverse problem is presented. Several applications of transversely isotropic elasticity imaging to cancellous bone from the human vertebra are then considered. The feasibility of using isotropic elasticity imaging to obtain meaningful reconstructions of the distribution of material properties for vertebral cancellous bone from experiment is established. However, using simulation, it is shown that an isotropic reconstruction is not appropriate for anisotropic materials. It is further shown that the transversely isotropic method identifies a solution that predicts the measured displacements, reveals regions of low stiffness, and recovers all five elastic parameters with approximately 10% error. The recovery of a given elastic parameter is found to require the presence of its corresponding strain (e.g., a deformation that generates ??? is necessary to reconstruct C????), and the application of regularization is shown to improve accuracy. Finally, the effects of noise on reconstruction quality is demonstrated and a signal-to-noise ratio (SNR) of 40 dB is identified as a reasonable threshold for obtaining accurate reconstructions from experimental data. This study demonstrates that given an appropriate set of displacement fields, level of regularization, and signal strength, the transversely isotropic method can recover the relative magnitudes of all five elastic parameters without an independent measurement of stress. The quality of the reconstructions improves with increasing contrast, magnitude of deformation, and asymmetry in the distributions of material properties, indicating that elasticity imaging of cancellous bone could be a useful tool in laboratory studies to monitor the progression of damage and disease in this tissue.     

A transversely isotropic hyperelastic constitutive model of the PDL. Analytical and computational aspects

This study describes the development of a constitutive law for the modelling of the periodontal ligament (PDL) and its practical implementation into a commercial finite element code. The constitutive equations encompass the essential mechanical features of this biological soft tissue: non-linear behaviour, large deformations, anisotropy, distinct behaviour in tension and compression and the fibrous characteristics. The approach is based on the theory of continuum fibre-reinforced composites at finite strain where a compressible transversely isotropic hyperelastic strain energy function is defined. This strain energy density function is further split into volumetric and deviatoric contributions separating the bulk and shear responses of the material. Explicit expressions of the stress tensors in the material and spatial configurations are first established followed by original expressions of the elasticity tensors in the material and spatial configurations. As a simple application of the constitutive model, two finite element analyses simulating the mechanical behaviour of the PDL are performed. The results highlight the significance of integrating the fibrous architecture of the PDL as this feature is shown to be responsible for the complex strain distribution observed.     

A Transversely Isotropic Hyperelastic Constitutive Model of the PDL. Analytical and Computational Aspects

This study describes the development of a constitutive law for the modelling of the periodontal ligament (PDL) and its practical implementation into a commercial finite element code. The constitutive equations encompass the essential mechanical features of this biological soft tissue: non-linear behaviour, large deformations, anisotropy, distinct behaviour in tension and compression and the fibrous characteristics. The approach is based on the theory of continuum fibre-reinforced composites at finite strain where a compressible transversely isotropic hyperelastic strain energy function is defined. This strain energy density function is further split into volumetric and deviatoric contributions separating the bulk and shear responses of the material. Explicit expressions of the stress tensors in the material and spatial configurations are first established followed by original expressions of the elasticity tensors in the material and spatial configurations. As a simple application of the constitutive model, two finite element analyses simulating the mechanical behaviour of the PDL are performed. The results highlight the significance of integrating the fibrous architecture of the PDL as this feature is shown to be responsible for the complex strain distribution observed.     

Microplane constitutive model and computational framework for blood vessel tissue

This paper presents a nonlinearly elastic anisotropic microplane formulation in 3D for computational constitutive modeling of arterial soft tissue in the passive regime. The constitutive modeling of arterial (and other biological) soft tissue is crucial for accurate finite element calculations, which in turn are essential for design of implants, surgical procedures, bioartificial tissue, as well as determination of effect of progressive diseases on tissues and implants. The model presented is defined at a lower scale (mesoscale) than the conventional macroscale and it incorporates the effect of all the (collagen) fibers which are anisotropic structural components distributed in all directions within the tissue material in addition to that of isotropic bulk tissue. It is shown that the proposed model not only reproduces Holzapfel's recent model but also improves on it by accounting for the actual three-dimensional distribution of fiber orientation in the arterial wall, which endows the model with advanced capabilities in simulation of remodeling of soft tissue. The formulation is flexible so that its parameters could be adjusted to represent the arterial wall either as a single material or a material composed of several layers in finite element analyses of arteries. Explicit algorithms for both the material subroutine and the explicit integration with dynamic relaxation of equations of motion using finite element method are given. To circumvent the slow convergence of the standard dynamic relaxation and small time steps dictated by the stability of the explicit integrator, an adaptive dynamic relaxation technique that ensures stability and fastest possible convergence rates is developed. Incompressibility is enforced using penalty method with an updated penalty parameter. The model is used to simulate experimental data from the literature demonstrating that the model response is in excellent agreement with the data. An experimental procedure to determine the distribution of fiber directions in 3D for biological soft tissue is suggested in accordance with the microplane concept. It is also argued that this microplane formulation could be modified or extended to model many other phenomena of interest in biomechanics.     

Persistent Symmetry Frustration in Pollen Tubes

Pollen tubes are extremely rapidly growing plant cells whose morphogenesis is determined by spatial gradients in the biochemical composition of the cell wall. We investigate the hypothesis (MP) that the distribution of the local mechanical properties of the wall, corresponding to the change of the radial symmetry along the axial direction, may lead to growth oscillations in pollen tubes. We claim that the experimentally observed oscillations originate from the symmetry change at the transition zone, where both intervening symmetries (cylindrical and spherical) meet. The characteristic oscillations between resonating symmetries at a given (constant) turgor pressure and a gradient of wall material constants may be identified with the observed growth-cycles in pollen tubes.