A viscoelastic anisotropic hyperelastic constitutive model of the human cornea

A constitutive model based on the continuum mechanics theory has been developed which represents interlamellar cohesion, regional variation of collagen fibril density, 3D anisotropy and both age-related viscoelastic and hyperelastic stiffening behaviour of the human cornea. Experimental data gathered from a number of previous studies on 48 ex vivo human cornea (inflation and shear tests) enabled calibration of the constitutive model by numerical analysis. Wide-angle X-ray scattering and electron microscopy provided measured data which quantify microstructural arrangements associated with stiffness. The present study measures stiffness parallel to the lamellae of the cornea which approximately doubles with an increase in strain rate from 0.5 to 5%/min, while the underlying stromal matrix provides a stiffness 2–3 orders of magnitude lower than the lamellae. The model has been simultaneously calibrated to within 3% error across three age groups ranging from 50 to 95 years and three strain rates across the two loading scenarios. Age and strain-rate-dependent material coefficients allow numerical simulation under varying loading scenarios for an individual patient with material stiffness approximated by their age. This present study addresses a significant gap in numerical representation of the cornea and has great potential in daily clinical practice for the planning and optimisation of corrective procedures and in preclinical optimisation of diagnostic procedures.     

The anisotropic material constitutive models for the human cornea

This paper presents an anisotropic analysis model for the human cornea. The model is based on the assumption that the fibrils in the cornea are organised into lamellae, which may have preferential orientation along the superior-inferior and nasal-temporal directions, while the alignment of lamellae with different orientations is assumed to be random. Hence, the cornea can be regarded as a laminated composite shell. The constitutive equation describing the relationships between membrane forces, bending moments, and membrane strains, bending curvatures are derived. The influences of lamella orientations and the random alignment of lamellae on the stiffness coefficients of the constitutive equation are discussed.     

Inverse identification of region-specific hyperelastic material parameters for human brain tissue

The identification of material parameters accurately describing the region-dependent mechanical behavior of human brain tissue is crucial for computational models used to assist, e.g., the development of safety equipment like helmets or the planning and execution of brain surgery. While the division of the human brain into different anatomical regions is well established, knowledge about regions with distinct mechanical properties remains limited. Here, we establish an inverse parameter identification scheme using a hyperelastic Ogden model and experimental data from multi-modal testing of tissue from 19 anatomical human brain regions to identify mechanically distinct regions and provide the corresponding material parameters. We assign the 19 anatomical regions to nine governing regions based on similar parameters and microstructures. Statistical analyses confirm differences between the regions and indicate that at least the corpus callosum and the corona radiata should be assigned different material parameters in computational models of the human brain. We provide a total of four parameter sets based on the two initial Poisson’s ratios of 0.45 and 0.49 as well as the pre- and unconditioned experimental responses, respectively. Our results highlight the close interrelation between the Poisson’s ratio and the remaining model parameters. The identified parameters will contribute to more precise computational models enabling spatially resolved predictions of the stress and strain states in human brains under complex mechanical loading conditions.

     

A general framework for the numerical implementation of anisotropic hyperelastic material models including non-local damage

The highly nonlinear mechanical behaviour of soft tissues solicited within the physiological range usually involves degradation of the material properties. Mechanically, having these biostructures undergoing such stretch patterns may bring about pathological conditions related to the steady deterioration of both collagen fibres and material’s ground substance. Tissue and subject variability observed in the phenomenological mechanical characterisation of soft tissues often hinder the choice of the computational constitutive model. Therefore, this contribution brings forth a detailed overview of the constitutive implementation in a computational framework of anisotropic hyperelastic materials with damage. Surmounting the challenge posed by the mesh dependency pathology requires the incorporation of an integral-type non-local averaging, which seeks to include the effects of the microstructure in order to limit the localisation phenomena of the damage variables. By adopting this approach, one can make use of multiple developed material models available in the literature, a combination of those, or even propose new models within the same numerical framework. The numerical examples of three-dimensional displacement and force-driven boundary value problems highlight the possibility of using multiple material models within the same numerical framework. Particularities concerning the considered material models and the damage effect implications to represent the Mullins effect, induced anisotropy, hysteresis, and mesh dependency are discussed.     

Mechanical identification of layer-specific properties of mouse carotid arteries using 3D-DIC and a hyperelastic anisotropic constitutive model

The role of mechanics is known to be of primary order in many arterial diseases; however, determining mechanical properties of arteries remains a challenge. This paper discusses the identifiability of the passive mechanical properties of a mouse carotid artery, taking into account the orientation of collagen fibres in the medial and adventitial layers. On the basis of 3D digital image correlation measurements of the surface strain during an inflation/extension test, an inverse identification method is set up. It involves a 3D finite element mechanical model of the mechanical test and an optimisation algorithm. A two-layer constitutive model derived from the Holzapfel model is used, with five and then seven parameters. The five-parameter model is successfully identified providing layer-specific fibre angles. The seven-parameter model is over parameterised, yet it is shown that additional data from a simple tension test make the identification of refined layer-specific data reliable.     

An inhomogeneous and anisotropic constitutive model of human dentin

Dentin constitutes the major part of human tooth. It is composed of a large number of tubules with both variational radii and radially parallel pattern. In addition, peritubular dentin surrounds each tubule lumen and has a higher elastic modulus than the matrix of dentin, i.e. intertubular dentin. Considering the above microstructural characteristics, a micromechanics model is used in this paper to evaluate the overall elastic properties of dentin. Five independent effective elastic parameters in transverse isotropic elasticity matrix can be expressed analytically by the material parameters of peri- and intertubular dentin and the volume fraction of tubules. To determine the effectivity of this theoretical model, a finite element (FE) model simulating a longitudinal tooth slice in moire fringe testing of Wang and Weiner (J. Biomech. 31 (1998) 135) was performed. Furthermore, the FE model was developed incorporate modeling of variation of tubule's diameter and softer characteristic of intertubular dentin near the dentin-enamel junction and around the pulp chamber. It turned out that the isoline figure of longitudinal displacement by FE calculation has very similar patterns to the moire fringe results. However, the FE results of displacement by traditional stress-strain models which regard dentin as a homogeneous and isotropic material show an obviously different strain distributions as compared to published moire fringes results. Thus the inhomogeneous and anisotropic model developed in this paper more accurately reflects the true physical nature of human dentin.     

Relationship between structural modeling and hyperelastic material behavior: application to CNS white matter

 Recent measurements of the material properties of brain tissue allow an examination of the underlying microstructural basis in both physiological and pathophysiological conditions. The purpose of this study is to develop a mathematical relationship between microstructurally based models of the central nervous system (CNS) white matter and equivalent hyperelastic material models. For simplicity, time dependent material behavior is not included in this formulation. The microstructural representation is used to formulate structural property relationships for highly oriented white matter, and is mathematically compared to one isotropic and two anisotropic hyperelastic formulations. For the anisotropic characterizations, the population of axons in the white matter is assumed to align along one preferred direction of the material, yielding a transversely isotropic formulation. Relatively simple strain–energy functions incorporating material anisotropy provide sufficient flexibility to model the nonlinear behavior predicted from structurally based models, although the tangential stiffness of the hyperelastic approaches does not follow completely the behavior of the structurally based formulations. This analysis is an initial step towards linking microstructural aspects of the tissue to material models commonly used for large deformations, and may be an important step in relating predicted tissue deformation to the deformation and stress of cellular and subcellular structures. Received: 15 October 2001 / Accepted: 30 September 2002 Funds for this work were provided by CDC grant R49/CCR312712 and NIH grants P50 NS08803, NICHD RO1 41699, and NINDS RO1 35712.     

Modeling and experimental investigation of electromechanical properties of scleral tissue; a CEM model using an anisotropic hyperelastic constitutive relation

Biomechanics and Modeling in Mechanobiology - The sclera is a soft tissue primarily consisting of collagen fibers, elastin, and proteoglycans. The proteoglycans are composed of a core protein and...     

Spherical indentation method for determining the constitutive parameters of hyperelastic soft materials

A comprehensive study on the spherical indentation of hyperelastic soft materials is carried out through combined theoretical, computational, and experimental efforts. Four widely used hyperelastic constitutive models are studied, including neo-Hookean, Mooney–Rivlin, Fung, and Arruda–Boyce models. Through dimensional analysis and finite element simulations, we establish the explicit relations between the indentation loads at given indentation depths and the constitutive parameters of materials. Based on the obtained results, the applicability of Hertzian solution to the measurement of the initial shear modulus of hyperelastic materials is examined. Furthermore, from the viewpoint of inverse problems, the possibility to measure some other properties of a hyperelastic material using spherical indentation tests, e.g., locking stretch, is addressed by considering the existence, uniqueness, and stability of the solution. Experiments have been performed on polydimethylsiloxane to validate the conclusions drawn from our theoretical analysis. The results reported in this study should help identify the extent to which the mechanical properties of hyperelastic materials could be measured from spherical indentation tests.     

Hyperelastic anisotropic microplane constitutive model for annulus fibrosus

In a recent paper, Peng et al. (2006, "An Anisotropic Hyperelastic Constitutive Model With Fiber-Matrix Interaction for the Human Annulus Fibrosis," ASME J. Appl. Mech., 73(5), pp. 815-824) developed an anisotropic hyperelastic constitutive model for the human annulus fibrosus in which fiber-matrix interaction plays a crucial role in simulating experimental observations reported in the literature. Later, Guo et al. (2006, "A Composites-Based Hyperelastic Constitutive Model for Soft Tissue With Application to the Human Fibrosis," J. Mech. Phys. Solids, 54(9), pp. 1952-1971) used fiber reinforced continuum mechanics theory to formulate a model in which the fiber-matrix interaction was simulated using only composite effect. It was shown in these studies that the classical anisotropic hyperelastic constitutive models for soft tissue, which do not account for this shear interaction, cannot accurately simulate the test data on human annulus fibrosus. In this study, we show that the microplane model for soft tissue developed by Caner and Carol (2006, "Microplane Constitutive Model and Computational Framework for Blood Vessel Tissue," ASME J. Biomech. Eng., 128(3), pp. 419-427) can be adjusted for human annulus fibrosus and the resulting model can accurately simulate the experimental observations without explicit fiber-matrix interaction because, in microplane model, the shear interaction between the individual fibers distributed in the tissue provides the required additional rigidity to explain these experimental facts. The intensity of the shear interaction between the fibers can be adjusted by adjusting the spread in the distribution while keeping the total amount of the fiber constant. A comparison of results obtained from (i) a fiber-matrix parallel coupling model, which does not account for the fiber-matrix interaction, (ii) the same model but enriched with fiber-matrix interaction, and (iii) microplane model for soft tissue adapted to annulus fibrosus with two families of fiber distributions is presented. The conclusions are (i) that varying degrees of fiber-fiber and fiber-matrix shear interaction must be taking place in the human annulus fibrosus, (ii) that this shear interaction is essential to be able to explain the mechanical behavior of human annulus fibrosus, and (iii) that microplane model can be fortified with fiber-matrix interaction in a straightforward manner provided that there are new experimental data on distribution of fibers, which indicate a spread so small that it requires an explicit fiber-matrix interaction to be able to simulate the experimental data.     

An anisotropic constitutive model for immersogeometric fluid–structure interaction analysis of bioprosthetic heart valves

This paper considers an anisotropic hyperelastic soft tissue model, originally proposed for native valve tissue and referred to herein as the Lee–Sacks model, in an isogeometric thin shell analysis framework that can be readily combined with immersogeometric fluid–structure interaction (FSI) analysis for high-fidelity simulations of bioprosthetic heart valves (BHVs) interacting with blood flow. We find that the Lee–Sacks model is well-suited to reproduce the anisotropic stress–strain behavior of the cross-linked bovine pericardial tissues that are commonly used in BHVs. An automated procedure for parameter selection leads to an instance of the Lee–Sacks model that matches biaxial stress–strain data from the literature more closely, over a wider range of strains, than other soft tissue models. The relative simplicity of the Lee–Sacks model is attractive for computationally-demanding applications such as FSI analysis and we use the model to demonstrate how the presence and direction of material anisotropy affect the FSI dynamics of BHV leaflets.     

Structural three-dimensional constitutive law for the passive myocardium

A three-dimensional constitutive law is proposed for the myocardium. Its formulation is based on a structural approach in which the total strain energy of the tissue is the sum of the strain energies of its constituents: the muscle fibers, the collagen fibers and the fluid matrix which embeds them. The ensuing material law expresses the specific structural and mechanical properties of the tissue, namely, the spatial orientation of the comprising fibers, their waviness in the unstressed state and their stress-strain behavior when stretched. Having assumed specific functional forms for the distribution of the fibers spatial orientation and waviness, the results of biaxial mechanical tests serve for the estimation of the material constants appearing in the constitutive equations. A very good fit is obtained between the measured and the calculated stresses, indicating the suitability of the proposed model for describing the mechanical behavior of the passive myocardium. Moreover, the results provide general conclusions concerning the structural basis for the tissue overall mechanical properties, the main of which is that the collagen matrix, though comprising a relatively small fraction of the whole tissue volume, is the dominant component accounting for its stiffness.     

Human annulus fibrosus material properties from biaxial testing and constitutive modeling are altered with degeneration

The annulus fibrosus (AF) of the intervertebral disk undergoes large and multidirectional stresses and strains. Uniaxial tensile tests are limited for measuring AF material properties, because freely contracting edges can prevent fiber stretch and are not representative of in situ boundary conditions. The objectives of this study were to measure human AF biaxial tensile mechanics and to apply and validate a constitutive model to determine material properties. Biaxial tensile tests were performed on samples oriented along the circumferential–axial and the radial–axial directions. Data were fit to a structurally motivated anisotropic hyperelastic model composed of isotropic extra-fibrillar matrix, nonlinear fibers, and fiber–matrix interactions (FMI) normal to the fibers. The validated model was used to simulate shear and uniaxial tensile behavior, to investigate AF structure–function, and to quantify the effect of degeneration. The biaxial stress–strain response was described well by the model (R ²?>?0.9). The model showed that the parameters for fiber nonlinearity and the normal FMI correlated with degeneration, resulting in an elongated toe-region and lower stiffness with degeneration. The model simulations in shear and uniaxial tension successfully matched previously published circumferential direction Young’s modulus, provided an explanation for the low values in previously published axial direction Young’s modulus, and was able to simulate shear mechanics. The normal FMI were important contributors to stress and changed with degeneration, therefore, their microstructural and compositional source should be investigated. Finally, the biaxial mechanical data and constitutive model can be incorporated into a disk finite element model to provide improved quantification of disk mechanics.     

A meshless microscale bone tissue trabecular remodelling analysis considering a new anisotropic bone tissue material law

In this work, a novel anisotropic material law for the mechanical behaviour of the bone tissue is proposed. This new law, based on experimental data, permits to correlate the bone apparent density with the obtained level of stress. Combined with the proposed material law, a biomechanical model for predicting bone density distribution was developed, based on the assumption that the bone structure is a gradually self-optimising anisotropic biological material that maximises its own structural stiffness. The strain and the stress field required in the iterative remodelling process are obtained by means of an accurate meshless method, the Natural Neighbour Radial Point Interpolation Method (NNRPIM). Comparing with other numerical approaches, the inclusion of the NNRPIM presents numerous advantages such as the high accuracy and the smoother stress and strain field distribution. The natural neighbour concept permits to impose organically the nodal connectivity and facilitates the analysis of convex boundaries and extremely irregular meshes. The viability and efficiency of the model were tested on several trabecular benchmark patch examples. The results show that the pattern of the local bone apparent density distribution and the anisotropic bone behaviour predicted by the model for the microscale analysis are in good agreement with the expected structural architecture and bone apparent density distribution.     

A three-dimensional elastic plastic damage constitutive law for bone tissue

Motivated by mechanical analysis of bones and bone-implant systems, a 3D constitutive law describing the macroscopic mechanical behaviour of both cortical and trabecular bone in cyclic (not fatigue) overloads is developed. The proposed model which mathematical formulation is established within the framework of generalized standard materials accounts for three distinct material evolution modes where elastic, plastic and damage aspects are closely related. The anisotropic elasticity of bone is described by a morphology-based model and distinct damage behaviour in tension and compression by a halfspacewise generalized Hill criterion. The plastic criterion is based on the intact elastic compliance tensor. The algorithm applies three distinct projections based on the relationship between the internal variables and criteria. Their respective consistent tangent operators are presented. Numerical resolutions of several boundary value problems and a biomechanical application are presented to illustrate the potential of the constitutive model and demonstrate the expected quadratic convergence of the algorithm.     

A visco-hyperelastic constitutive approach for modeling polyvinyl alcohol sponge

This study proposes the quasi-linear viscoelastic (QLV) model to characterize the time dependent mechanical behavior of poly(vinyl alcohol) (PVA) sponges. The PVA sponges have implications in many viscoelastic soft tissues, including cartilage, liver, and kidney as an implant. However, a critical barrier to the use of the PVA sponge as tissue replacement material is a lack of sufficient study on its viscoelastic mechanical properties. In this study, the nonlinear mechanical behavior of a fabricated PVA sponge is investigated experimentally and computationally using relaxation and stress failure tests as well as finite element (FE) modeling. Hyperelastic strain energy density functions, such as Yeoh and Neo-Hookean, are used to capture the mechanical behavior of PVA sponge at ramp part, and viscoelastic model is used to describe the viscose behavior at hold part. Hyperelastic material constants are obtained and their general prediction ability is verified using FE simulations of PVA tensile experiments. The results of relaxation and stress failure tests revealed that Yeoh material model can define the mechanical behavior of PVA sponge properly compared with Neo-Hookean one. FE modeling results are also affirmed the appropriateness of Yeoh model to characterize the mechanical behavior of PVA sponge. Thus, the Yeoh model can be used in future biomechanical simulations of the spongy biomaterials. These results can be utilized to understand the viscoelastic behavior of PVA sponges and has implications for tissue engineering as scaffold.     

Determining constitutive behavior of the brain tissue using digital image correlation and finite element modeling

Biomechanics and Modeling in Mechanobiology - Detailed knowledge about the mechanical properties of brain can improve numerical modeling of the brain under various loading conditions. The success...     

Analysis of prolapse in cardiovascular stents: a constitutive equation for vascular tissue and finite-element modelling

The effectiveness of a cardiovascular stent depends on many factors, such as its ability to sustain the compression applied by the vessel wall, minimal longitudinal contraction when it is expanded, and its ability to flex when navigating tortuous blood vessels. The long-term reaction of the tissue to the stent is also device dependant; in particular some designs provoke in-stent restenosis (i.e., regrowth of the occlusion around the stent). The mechanism of restenosis is thought to involve injury or damage to the vessel wall due to the high stresses generated around the stent when it expands. Because of this, the deflection of the tissue between the struts of the stent (called prolapse or "draping") has been used as a measure of the potential of a stent to cause restenosis. In this paper, uniaxial and biaxial experiments on human femoral artery and porcine aortic vascular tissue are used to develop a hyperelastic constitive model of vascular tissue suitable for implementation in finite-element analysis. To analyze prolapse, four stent designs (BeStent 2, Medtronic AVE; NIROYAL, Boston Scientific; VELOCITY, Cordis; TETRA, Guidant) were expanded in vitro to determine their repeating-unit dimensions. This geometric data was used to generate a finite element model of the vascular tissue supported within a repeating-unit of the stent. Under a pressure of 450 mm Hg (representing the radial compression of the vessel wall), maximum radial deflection of 0.253 mm, 0.279 mm, 0.348 mm and 0.48 mm were calculated for each of the four stents. Stresses in the vascular wall were highest for the VELOCITY stent. The method is proposed as a way to compare stents relative to their potential for restenosis and as a basis for a biomechanical design of a stent repeating-unit that would minimize restenosis.     

Experimentally tractable,pseudo-elastic constitutive law for biomembranes: I. Theory

Although visco-elastic in general, the stress-strain relation of biomembranes is one-to-one or pseudo-elastic when being loaded after preconditioning. This pseudo-elastic relation is hypoelastic (i.e., it is not hyperelastic), yet much of the stress response can be characterized by a scalar function omega that represents the work done (per unit reference volume) on the specimen during loading. (Since a pseudo-strain-energy function W is optimized to fit the test data and not the work done, omega is not equal to W in general.) The remaining part tR of the stress response does no work during loading. With biaxial testing, omega can be definitively determined from data. Moreover, for tests with the stretch directions coaxial to the axes of anisotropy, tR can be accurately characterized by a scalar function omega that depends on the strain. This paper is part 1 of 2 with "I. Theory" and "II. Application."