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.     

Nonlinear elastic analysis of blood vessels

Based on the theory of Green and Adkins [9], a strain energy function is proposed to describe the nonlinear mechanical behavior of arteries. The arterial tissue is assumed to be a nonlinear elastic, incompressible material with local triclinicity and transverse isotropy. Although the arterial tissue shows viscous phenomena, experimental results have indicated that viscosity is only a second-order effect as compared to the nonlinear elasticity of the tissue. The advantage of the formulation presented herein is that it is relatively simple and results in an accurate stress-strain relation that can be used readily for the study of wave propagations in the blood vessels. For nonlinear finite strain elasticity of the order two, ten elastic constants are needed to describe the material nonlinearity of the arterial tissue. Based on the orthogonal, transverse isotropies and the incompressibility conditions, ten constraint equations may be established and the elastic constants can be uniquely determined by correlating with the experimental results. An example of calculating these elastic constants is made by using the experimental data of Patel, et al. [14-17] for the intercoastal arteries in living dogs. The predicted mechanical behavior of canine arteries is quite satisfactory as compared to the experimental data except when the longitudinal and the circumferential stretches exceed 1.60. However, such a strain magnitude may not be expected in in-vivo arteries because of the constraints of peripheral connecting tissues. Otherwise, the strain energy function including higher order strain terms should be used.(ABSTRACT TRUNCATED AT 250 WORDS)     

Myocardial material parameter estimation: a comparison of invariant based orthotropic constitutive equations

This study investigated a number of invariant based orthotropic and transversely isotropic constitutive equations for their suitability to fit three-dimensional simple shear mechanics data of passive myocardial tissue.

A number of orthotropic laws based on Green strain components and one microstructurally based law have previously been investigated to fit experimental measurements of stress-strain behaviour. Here we extend this investigation to include several recently proposed functional forms, i.e. invariant based orthotropic and transversely isotropic constitutive relations.

These laws were compared on the basis of (i) ‘goodness of fit’: how well they fit a set of six shear deformation tests, (ii) ‘variability’: how well determined the material parameters are over the range of experiments. These criteria were utilised to discuss the advantages and disadvantages of the constitutive laws.

It was found that a specific form of the polyconvex type as well as the exponential Fung-type law from the previous study were most suitable for modelling the orthotropic behaviour of myocardium under simple shear.     

Biaxial elastic material properties of porcine coronary media and adventitia

The importance of mechanical stresses and strains has become well recognized in vascular physiology and pathology. To compute the stress and strain on the various components of the vessel wall, we must know the constitutive equations for the different layers of the vessel wall. The objective of the present study is to determine the constitutive equation of the coronary artery treated as a two-layer composite: intima-media and adventitial layers. Twelve hearts were obtained from a local slaughterhouse, and the right coronary artery and left anterior descending artery were dissected free from the myocardium. The vessel wall was initially mechanically tested biaxially (inflation and axial extension) as a whole (intact wall) and subsequently as intima-media or adventitial layer. A Fung-type exponential strain energy function was used to curve fit the experimental data for the intact wall and individual layers for the right coronary artery and left anterior descending artery. Two methods were used for the determination of material constants, including the Marquardt-Levenberg nonlinear least squares method and the genetic algorithm method. Our results show that there were no statistically significant differences in the material constants obtained from the two methods and that either set of elastic constants results in good fit of the data. Furthermore, at an in vivo value of axial stretch ratio, we find that the stiffness is as follows: intima-media > intact > adventitia. These results underscore the composite nature of coronary arteries with different material properties in each layer. The present results are necessary for analysis of coronary artery mechanics and to provide a fundamental understanding of vessel physiology.     

Estimation of the effective transversely isotropic elastic constants of a material from known values of the material's orthotropic elastic constants

A method is illustrated for determining the effective transversely isotropic (or isotropic) elastic constants from measured orthotropic elastic constants. This method consists of constructing upper and lower bounds on the effective transversely isotropic (or isotropic) elastic constants using the known orthotropic values. This method is illustrated using three sets of elastic constants for bone. Fortunately, the upper and lower bounds are very close. Thus very good approximations for the effective transversely isotropic (or isotropic) elastic constants for cortical and cancellous bone are obtained from previously published data on the orthotropic elastic constants for those tissue types. This work is undertaken to build a greater database for the transversely isotropic elastic constants of bone with the intention of employing them in a transversely isotropic model of bone poroelasticity. An interesting aspect of the present result is that the Voigt and Reuss bounds are very tight for these anisotropic materials. This is not always the case for these bounds. Received: 14 November 2001 / Accepted: 25 February 2002     

The validation of a generalized Hooke's law for coronary arteries

The exponential form of constitutive model is widely used in biomechanical studies of blood vessels. There are two main issues, however, with this model: 1) the curve fits of experimental data are not always satisfactory, and 2) the material parameters may be oversensitive. A new type of strain measure in a generalized Hooke's law for blood vessels was recently proposed by our group to address these issues. The new model has one nonlinear parameter and six linear parameters. In this study, the stress-strain equation is validated by fitting the model to experimental data of porcine coronary arteries. Material constants of left anterior descending artery and right coronary artery for the Hooke's law were computed with a separable nonlinear least-squares method with an excellent goodness of fit. A parameter sensitivity analysis shows that the stability of material constants is improved compared with the exponential model and a biphasic model. A boundary value problem was solved to demonstrate that the model prediction can match the measured arterial deformation under experimental loading conditions. The validated constitutive relation will serve as a basis for the solution of various boundary value problems of cardiovascular biomechanics.     

A finite element study of invariant-based orthotropic constitutive equations in the context of myocardial material parameter estimation

A previous study investigated a number of invariant-based orthotropic and transversely isotropic constitutive equations for their suitability to fit three-dimensional simple shear mechanics data of passive myocardial tissue. The study was based on the assumption of a homogeneous deformation. Here, we extend the previous study by performing an inverse finite element material parameter estimation. This ensures a more realistic deformation state and material parameter estimates. The constitutive relations were compared on the basis of (i) ‘goodness of fit’: how well they fit a set of six shear deformation tests and (ii) ‘variability’: how well determined the material parameters are over the range of experiments. These criteria were utilised to discuss the advantages and disadvantages of the constitutive relations. It was found that a specific form of the polyconvex type as well as the exponential Fung-type equations were most suitable for modelling the orthotropic behaviour of myocardium under simple shear.     

A Conewise Linear Elasticity mixture model for the analysis of tension-compression nonlinearity in articular cartilage

A biphasic mixture model is developed that can account for the observed tension-compression nonlinearity of cartilage by employing the continuum-based Conewise Linear Elasticity (CLE) model of Curnier et al. (J. Elasticity, 37, 1-38, 1995) to describe the solid phase of the mixture. In this first investigation, the orthotropic octantwise linear elasticity model was reduced to the more specialized case of cubic symmetry, to reduce the number of elastic constants from twelve to four. Confined and unconfined compression stress-relaxation, and torsional shear testing were performed on each of nine bovine humeral head articular cartilage cylindrical plugs from 6 month old calves. Using the CLE model with cubic symmetry, the aggregate modulus in compression and axial permeability were obtained from confined compression (H-A = 0.64 +/- 0.22 MPa, k2 = 3.62 +/- 0.97 x 10(-16) m4/N.s, r2 = 0.95 +/- 0.03), the tensile modulus, compressive Poisson ratio, and radial permeability were obtained from unconfined compression (E+Y = 12.75 +/- 1.56 MPa, v- = 0.03 +/- 0.01, kr = 6.06 +/- 2.10 x 10(-16) m4/N.s, r2 = 0.99 +/- 0.00), and the shear modulus was obtained from torsional shear (mu = 0.17 +/- 0.06 MPa). The model was also employed to predict the interstitial fluid pressure successfully at the center of the cartilage plug in unconfined compression (r2 = 0.98 +/- 0.01). The results of this study demonstrate that the integration of the CLE model with the biphasic mixture theory can provide a model of cartilage that can successfully curve-fit three distinct testing configurations while producing material parameters consistent with previous reports in the literature.     

Effects of damage on the orthotropic material symmetry of bovine tibial trabecular bone

The macroscopic mechanical properties of trabecular bone can be predicted by its architecture using theoretical relationships between the elastic and architectural properties. Microdamage caused by overloading or fatigue decreases the apparent elastic moduli of trabecular bone requiring these relationships to be modified to predict the damaged elastic properties. In the case of isotropic damage, the apparent level elastic properties could be determined by multiplying all of the elastic constants by a single scalar factor. If the damage is anisotropic, the elastic constants may change by differing factors and the material coordinate system could become misaligned with the fabric coordinate system. High-resolution finite element models were used to simulate damage overloading on seven trabecular bone specimens subjected to pure shear strain in two planes. Comparison of the apparent elastic moduli of the specimens before and after damage showed that the reduction of the elastic moduli was anisotropic. This suggests that the microdamage within the specimens was inhomogeneous. However, after damage the specimens exhibited nearly orthotropic material symmetry as they did before damage. Changes in the orientation of the orthotropic material coordinate system were also small and occurred primarily in the transverse plane. Thus, while damage in trabecular bone is anisotropic, the material coordinate system remains aligned with the fabric tensor.     

Ab initio predictions of structural and elastic properties of struvite: contribution to urinary stone research

In the present work, we carried out density functional calculations of struvite--the main component of the so-called infectious urinary stones--to study its structural and elastic properties. Using a local density approximation and a generalised gradient approximation, we calculated the equilibrium structural parameters and elastic constants C(ijkl). At present, there is no experimental data for these elastic constants C (ijkl) for comparison. Besides the elastic constants, we also present the calculated macroscopic mechanical parameters, namely the bulk modulus (K), the shear modulus (G) and Young's modulus (E). The values of these moduli are found to be in good agreement with available experimental data. Our results imply that the mechanical stability of struvite is limited by the shear modulus, G. The study also explores the energy-band structure to understand the obtained values of the elastic constants.     

A Bimodular Theory for Finite Deformations: Comparison of Orthotropic Second-order and Exponential Stress Constitutive Equations for Articular Cartilage

Cartilaginous tissues, such as articular cartilage and the annulus fibrosus, exhibit orthotropic behavior with highly asymmetric tensile–compressive responses. Due to this complex behavior, it is difficult to develop accurate stress constitutive equations that are valid for finite deformations. Therefore, we have developed a bimodular theory for finite deformations of elastic materials that allows the mechanical properties of the tissue to differ in tension and compression. In this paper, we derive an orthotropic stress constitutive equation that is second-order in terms of the Biot strain tensor as an alternative to traditional exponential type equations. Several reduced forms of the bimodular second-order equation, with six to nine parameters, and a bimodular exponential equation, with seven parameters, were fit to an experimental dataset that captures the highly asymmetric and orthotropic mechanical response of cartilage. The results suggest that the bimodular second-order models may be appealing for some applications with cartilaginous tissues.     

Analysis of indentation: implications for measuring mechanical properties with atomic force microscopy.

Indentation using the atomic force microscope (AFM) has potential to measure detailed micromechanical properties of soft biological samples. However, interpretation of the results is complicated by the tapered shape of the AFM probe tip, and its small size relative to the depth of indentation. Finite element models (FEMs) were used to examine effects of indentation depth, tip geometry, and material nonlinearity and heterogeneity on the finite indentation response. Widely applied infinitesimal strain models agreed with FEM results for linear elastic materials, but yielded substantial errors in the estimated properties for nonlinear elastic materials. By accounting for the indenter geometry to compute an apparent elastic modulus as a function of indentation depth, nonlinearity and heterogeneity of material properties may be identified. Furthermore, combined finite indentation and biaxial stretch may reveal the specific functional form of the constitutive law--a requirement for quantitative estimates of material constants to be extracted from AFM indentation data.     

Sensitivities of medial meniscal motion and deformation to material properties of articular cartilage, meniscus and meniscal attachments using design of experiments methods

This study investigated the role of the material properties assumed for articular cartilage, meniscus and meniscal attachments on the fit of a finite element model (FEM) to experimental data for meniscal motion and deformation due to an anterior tibial loading of 45 N in the anterior cruciate ligament-deficient knee. Taguchi style L18 orthogonal arrays were used to identify the most significant factors for further examination. A central composite design was then employed to develop a mathematical model for predicting the fit of the FEM to the experimental data as a function of the material properties and to identify the material property selections that optimize the fit. The cartilage was modeled as isotropic elastic material, the meniscus was modeled as transversely isotropic elastic material, and meniscal horn and the peripheral attachments were modeled as noncompressive and nonlinear in tension spring elements. The ability of the FEM to reproduce the experimentally measured meniscal motion and deformation was most strongly dependent on the initial strain of the meniscal horn attachments (epsilon(1H)), the linear modulus of the meniscal peripheral attachments (E(P)) and the ratio of meniscal moduli in the circumferential and transverse directions (E(theta)E(R)). Our study also successfully identified values for these critical material properties (epsilon(1H) = -5%, E(P) = 5.6 MPa, E(theta)E(R) = 20) to minimize the error in the FEM analysis of experimental results. This study illustrates the most important material properties for future experimental studies, and suggests that modeling work of meniscus, while retaining transverse isotropy, should also focus on the potential influence of nonlinear properties and inhomogeneity.     

An orthotropic viscoelastic material model for passive myocardium: theory and algorithmic treatment

This contribution presents a novel constitutive model in order to simulate an orthotropic rate-dependent behaviour of the passive myocardium at finite strains. The motivation for the consideration of orthotropic viscous effects in a constitutive level lies in the disagreement between theoretical predictions and experimentally observed results. In view of experimental observations, the material is deemed as nearly incompressible, hyperelastic, orthotropic and viscous. The viscoelastic response is formulated by means of a rheological model consisting of a spring coupled with a Maxwell element in parallel. In this context, the isochoric free energy function is decomposed into elastic equilibrium and viscous non-equilibrium parts. The baseline elastic response is modelled by the orthotropic model of Holzapfel and Ogden [Holzapfel GA, Ogden RW. 2009. Constitutive modelling of passive myocardium: a structurally based framework for material characterization. Philos Trans Roy Soc A Math Phys Eng Sci. 367:3445–3475]. The essential aspect of the proposed model is the account of distinct relaxation mechanisms for each orientation direction. To this end, the non-equilibrium response of the free energy function is constructed in the logarithmic strain space and additively decomposed into three anisotropic parts, denoting fibre, sheet and normal directions each accompanied by a distinct dissipation potential governing the evolution of viscous strains associated with each orientation direction. The evolution equations governing the viscous flow have an energy-activated nonlinear form. The energy storage in the Maxwell branches has a quadratic form leading to a linear stress–strain response in the logarithmic strain space. On the numerical side, the algorithmic aspects suitable for the implicit finite element method are discussed in a Lagrangian setting. The model shows excellent agreement compared to experimental data obtained from the literature. Furthermore, the finite element simulations of a heart cycle carried out with the proposed model show significant deviations in the strain field relative to the elastic solution.     

Orthotropic properties of cancellous bone modelled as parameterized cellular material

Analysis of stresses and strains in bone tissues and simulation of their adaptive remodelling require exhaustive information about distribution of constitutive properties of cancellous bone and their relationships to microstructural parameters. Homogenization of "equivalent" trabecular microstructures appears to be an advantageous tool for this task. In this study, parameterized orthotropic constitutive models of cancellous bone are derived from finite element analysis of repeatable microstructure cells. The models, based on a space-filling dodecahedron, are fully three-dimensional and are parameterized with four shape parameters. Variation of the parameters allows to imitate most of typical microstructure patterns observed in real bones, along with a variety of intermediate geometries. Finite element models of cells are generated by a special-purpose structured mesh generator for any arbitrary set of shape parameter values. Static numerical tests are performed for an exhaustive number of parameter value sets (microstructure instances). Coefficients of elastic orthotropic stiffness matrix are determined as tabularized functions of elastic constants versus the shape parameters. Additionally, they are correlated to apparent density and principal fabric tensor values. Comparison of the results with micro-FE data obtained for a large set of cancellous bone specimens proves a good agreement.     

Orthotropic properties of cancellous bone modelled as parameterized cellular material

Analysis of stresses and strains in bone tissues and simulation of their adaptive remodelling require exhaustive information about distribution of constitutive properties of cancellous bone and their relationships to microstructural parameters. Homogenization of “equivalent” trabecular microstructures appears to be an advantageous tool for this task. In this study, parameterized orthotropic constitutive models of cancellous bone are derived from finite element analysis of repeatable microstructure cells. The models, based on a space-filling dodecahedron, are fully three-dimensional and are parameterized with four shape parameters. Variation of the parameters allows to imitate most of typical microstructure patterns observed in real bones, along with a variety of intermediate geometries. Finite element models of cells are generated by a special-purpose structured mesh generator for any arbitrary set of shape parameter values. Static numerical tests are performed for an exhaustive number of parameter value sets (microstructure instances). Coefficients of elastic orthotropic stiffness matrix are determined as tabularized functions of elastic constants versus the shape parameters. Additionally, they are correlated to apparent density and principal fabric tensor values. Comparison of the results with micro-FE data obtained for a large set of cancellous bone specimens proves a good agreement.     

The importance of nonlinear tissue modelling in finite element simulations of infant head impacts

Despite recent efforts on the development of finite element (FE) head models of infants, a model capable of capturing head responses under various impact scenarios has not been reported. This is hypothesized partially attributed to the use of simplified linear elastic models for soft tissues of suture, scalp and dura. Orthotropic elastic constants are yet to be determined to incorporate the direction-specific material properties of infant cranial bone due to grain fibres radiating from the ossification centres. We report here on our efforts in advancing the above-mentioned aspects in material modelling in infant head and further incorporate them into subject-specific FE head models of a newborn, 5- and 9-month-old infant. Each model is subjected to five impact tests (forehead, occiput, vertex, right and left parietal impacts) and two compression tests. The predicted global head impact responses of the acceleration–time impact curves and the force–deflection compression curves for different age groups agree well with the experimental data reported in the literature. In particular, the newly developed Ogden hyperelastic model for suture, together with the nonlinear modelling of scalp and dura mater, enables the models to achieve more realistic impact performance compared with linear elastic models. The proposed approach for obtaining age-dependent skull bone orthotropic material constants counts both an increase in stiffness and decrease in anisotropy in the skull bone—two essential biological growth parameters during early infancy. The profound deformation of infant head causes a large stretch at the interfaces between the skull bones and the suture, suggesting that infant skull fractures are likely to initiate from the interfaces; the impact angle has a profound influence on global head impact responses and the skull injury metrics for certain impact locations, especially true for a parietal impact.     

Thermodynamic restrictions on the elastic constants of bone

The thermodynamic restrictions on the elastic coefficients of linear orthotropic elasticity and linear transversely isotropy elasticity are recorded and it is shown that previously reported data for the elastic orthotropic constants of bone satisfy these thermodynamic restrictions.     

Sensitivity analysis and parametric study of elastic properties of an unidirectional mineralized bone fibril-array using mean field methods

The key parameters determining the elastic properties of an unidirectional mineralized bone fibril-array decomposed in two further hierarchical levels are investigated using mean field methods. Modeling of the elastic properties of mineralized micro- and nanostructures requires accurate information about the underlying topology and the constituents’ material properties. These input data are still afflicted by great uncertainties and their influence on computed elastic constants of a bone fibril-array remains unclear. In this work, mean field methods are applied to model mineralized fibrils, the extra-fibrillar matrix and the resulting fibril-array. The isotropic or transverse isotropic elastic constants of these constituents are computed as a function of degree of mineralization, mineral distribution between fibrils and extra-fibrillar matrix, collagen stiffness and fibril volume fraction. The linear sensitivity of the elastic constants was assessed at a default set of the above parameters. The strain ratios between the constituents as well as the axial and transverse indentation moduli of the fibril-array were calculated for comparison with experiments. Results indicate that the degree of mineralization and the collagen stiffness dominate fibril-array elasticity. Interestingly, the stiffness of the extra-fibrillar matrix has a strong influence on transverse and shear moduli of the fibril-array. The axial strain of the intra-fibrillar mineral platelets is 30–90% of the applied fibril strain, depending on mineralization and collagen stiffness. The fibril-to-fibril-array strain ratio is essentially ~1. This study provides an improved insight in the parameters, which govern the fibril-array stiffness of mineralized tissues such as bone.     

Elastic properties of cancellous bone derived from finite element models of parameterized microstructure cells

Evaluation of constitutive properties of cancellous bone and their relationships to microstructural parameters is a crucial issue in analysis of stresses and strains in bone tissues and simulation of their remodelling. Known limitations of experimental methods as well as of the micro-FE techniques make the analysis and homogenization of 'equivalent' trabecular microstructures an advantageous tool for this task. In this study, parameterized orthotropic constitutive models of cancellous bone are derived from finite element analysis of repeatable microstructure cells. Two cell types are analysed: cube- and prism-based. The models are fully three-dimensional, have realistic curvilinear shapes and are parameterized with three shape parameters. Variation of the parameters allows to imitate most of the typical microstructure patterns observed in real bones, along with variety of intermediate geometries. Finite element models of cells are generated by a special-purpose structured mesh generator for any arbitrary set of shape parameter values. Six static numerical tests are performed for an exhaustive number of parameter value sets (microstructure instances). Multi-point boundary conditions imposed on the models ensure mutual fitting of deformed neighbouring cells. Values of computed stresses allow to determine all coefficients of elastic orthotropic stiffness matrix. Results have a form of tabularized functions of elastic constants versus the shape parameters. Comparison of the results with micro-FE data obtained for a large set of cancellous bone specimens proves a good agreement, though evidently better in the case of the prism-based cell model.