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.     

Method for characterizing viscoelasticity of human gluteal tissue

Characterizing compressive transient large deformation properties of biological tissue is becoming increasingly important in impact biomechanics and rehabilitation engineering, which includes devices interfacing with the human body and virtual surgical guidance simulation. Individual mechanical in vivo behaviour, specifically of human gluteal adipose and passive skeletal muscle tissue compressed with finite strain, has, however, been sparsely characterised. Employing a combined experimental and numerical approach, a method is presented to investigate the time-dependent properties of in vivo gluteal adipose and passive skeletal muscle tissue. Specifically, displacement-controlled ramp-and-hold indentation relaxation tests were performed and documented with magnetic resonance imaging. A time domain quasi-linear viscoelasticity (QLV) formulation with Prony series valid for finite strains was used in conjunction with a hyperelastic model formulation for soft tissue constitutive model parameter identification and calibration of the relaxation test data. A finite element model of the indentation region was employed. Strong non-linear elastic but linear viscoelastic tissue material behaviour at finite strains was apparent for both adipose and passive skeletal muscle mechanical properties with orthogonal skin and transversal muscle fibre loading. Using a force-equilibrium assumption, the employed material model was well suited to fit the experimental data and derive viscoelastic model parameters by inverse finite element parameter estimation. An individual characterisation of in vivo gluteal adipose and muscle tissue could thus be established. Initial shear moduli were calculated from the long-term parameters for human gluteal skin/fat: G(∞,S/F)=1850 Pa and for cross-fibre gluteal muscle tissue: G(∞,M)=881 Pa. Instantaneous shear moduli were found at the employed ramp speed: G(0,S/F)=1920 Pa and G(0,M)=1032 Pa.     

A transversely isotropic viscoelastic constitutive equation for brainstem undergoing finite deformation

The objective of this study was to define the constitutive response of brainstem undergoing finite shear deformation. Brainstem was characterized as a transversely isotropic viscoelastic material and the material model was formulated for numerical implementation. Model parameters were fit to shear data obtained in porcine brainstem specimens undergoing finite shear deformation in three directions: parallel, perpendicular, and cross sectional to axonal fiber orientation and determined using a combined approach of finite element analysis (FEA) and a genetic algorithm (GA) optimizing method. The average initial shear modulus of brainstem matrix of 4-week old pigs was 12.7 Pa, and therefore the brainstem offers little resistance to large shear deformations in the parallel or perpendicular directions, due to the dominant contribution of the matrix in these directions. The fiber reinforcement stiffness was 121.2 Pa, indicating that brainstem is anisotropic and that axonal fibers have an important role in the cross-sectional direction. The first two leading relative shear relaxation moduli were 0.8973 and 0.0741, respectively, with corresponding characteristic times of 0.0047 and 1.4538 s, respectively, implying rapid relaxation of shear stresses. The developed material model and parameter estimation technique are likely to find broad applications in neural and orthopaedic tissues.     

Dynamic finite element implementation of nonlinear,anisotropic hyperelastic biological membranes

We present a novel method for the implementation of hyperelastic finite strain, non-linear strain-energy functions for biological membranes in an explicit finite element environment. The technique is implemented in LS-DYNA but may also be implemented in any suitable non-linear explicit code. The constitutive equations are implemented on the foundation of a co-rotational uniformly reduced Hughes-Liu shell. This shell is based on an updated-Lagrangian formulation suitable for relating Cauchy stress to the rate-of-deformation, i.e. hypo-elasticity. To accommodate finite deformation hyper-elastic formulations, a co-rotational deformation gradient is assembled over time, resulting in a formulation suitable for pseudo-hyperelastic constitutive equations that are standard assumptions in biomechanics. Our method was validated by comparison with (1) an analytic solution to a spherically-symmetric dynamic membrane inflation problem, incorporating a Mooney-Rivlin hyperelastic equation and (2) with previously published finite element solutions to a non-linear transversely isotropic inflation problem. Finally, we implemented a transversely isotropic strain-energy function for mitral valve tissue. The method is simple and accurate and is believed to be generally useful for anyone who wishes to model biologic membranes with an experimentally driven strain-energy function.     

A dual optimization method for the material parameter identification of a biphasic poroviscoelastic hydrogel: Potential application to hypercompliant soft tissues

A dual-indentation creep and stress relaxation methodology was developed and validated for the material characterization of very soft biological tissue within the framework of the biphasic poroviscoelastic (BPVE) constitutive model. Agarose hydrogel, a generic porous medium with mobile fluid, served as a mechanical tissue analogue for validation of the experimental procedure. Indentation creep and stress relaxation tests with a solid plane-ended cylindrical indenter were performed at identical sites on a gel sample with dimensions large enough with respect to indenter size in order to satisfy an infinite layer assumption. A finite element (FE) formulation coupled to a global optimization algorithm was utilized to simultaneously curve-fit the creep and stress relaxation data and extract the BPVE model parameters for the agarose gel. A numerical analysis with artificial data was conducted to validate the uniqueness of the computational procedure. The BPVE model was able to successfully cross-predict both creep and stress relaxation behavior for each pair of experiments with a single unique set of material parameters. Optimized elastic moduli were consistent with those reported in the literature for agarose gel. With the incorporation of appropriately-sized indenters to satisfy more stringent geometric constraints, this simple yet powerful indentation methodology can provide a straightforward means by which to obtain the BPVE model parameters of biological soft tissues that are difficult to manipulate (such as brain and adipose) while maintaining a realistic in situ loading environment.     

Efficient isogeometric thin shell formulations for soft biological materials

This paper presents three different constitutive approaches to model thin rotation-free shells based on the Kirchhoff–Love hypothesis. One approach is based on numerical integration through the shell thickness while the other two approaches do not need any numerical integration and so they are computationally more efficient. The formulation is designed for large deformations and allows for geometrical and material nonlinearities, which makes it very suitable for the modeling of soft tissues. Furthermore, six different isotropic and anisotropic material models, which are commonly used to model soft biological materials, are examined for the three proposed constitutive approaches. Following an isogeometric approach, NURBS-based finite elements are used for the discretization of the shell surface. Several numerical examples are investigated to demonstrate the capabilities of the formulation. Those include the contact simulation during balloon angioplasty.     

Finite element implementation of a generalized Fung-elastic constitutive model for planar soft tissues

Numerical simulations of the anisotropic mechanical properties of soft tissues and tissue-derived biomaterials using accurate constitutive models remain an important and challenging research area in biomechanics. While most constitutive modeling efforts have focused on the characterization of experimental data, only limited studies are available on the feasibility of utilizing those models in complex computational applications. An example is the widely utilized exponential constitutive model proposed by Fung. Although present in the biomechanics literature for several decades, implementation of this model into finite element (FE) simulations has been limited. A major reason for limited numerical implementations are problems associated with inherent numerical instability and convergence. To address this issue, we developed and applied two restrictions for a generalized Fung-elastic constitutive model necessary to achieve numerical stability. These are (1) convexity of the strain energy function, and (2) the condition number of material stiffness matrix set lower than a prescribed value. These constraints were implemented in the nonlinear regression used for constitutive model parameter estimation to the experimental biaxial mechanical data. We then implemented the generalized Fung-elastic model into a commercial FE code (ABAQUS, Pawtucket, RI, USA). Single element and multi-element planar biaxial test simulations were conducted to verify the accuracy and robustness of the implementation. Results indicated that numerical convergence and accurate FE implementation were consistently obtained. The present study thus presents an integrated framework for accurate and robust implementation of pseudo-elastic constitutive models for planar soft tissues. Moreover, since our approach is formulated within a general FE code, it can be straightforwardly adopted across multiple software platforms.     

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.     

Static indentation of anisotropic biomaterials using axially asymmetric indenters--a computational study

Indentation has historically been used by biomechanicians to extract the small strain elastic or viscoelastic properties of biological tissues. Because of the axisymmetry of indenters used in these studies however, analysis of the results requires the assumption of material isotropy and often yields an "effective" elastic modulus. Since most biological tissues such as bone and myocardium are known to be anisotropic, the use of conventional indentation techniques for estimating material properties is therefore limited. The feasibility of using an axially asymmetric indenter to determine material directions and in-plane material properties for anisotropic tissue is explored here using finite element analysis. The load versus displacement curves as would be measured by an indenter depend on the orientation of the indenter cross section relative to the in-plane material axes, thus suggesting a method for determining the underlying material directions. Additionally, the stiffness of the tissue response to indentation is sensitive to the values of the in-plane anisotropic material properties and prestretches, and thus test results can be used to back out relevant constitutive parameters.     

Modeling Soft Tissue Damage and Failure Using a Combined Particle/Continuum Approach

Biological soft tissues experience damage and failure as a result of injury, disease, or simply age; examples include torn ligaments and arterial dissections. Given the complexity of tissue geometry and material behavior, computational models are often essential for studying both damage and failure. Yet, because of the need to account for discontinuous phenomena such as crazing, tearing, and rupturing, continuum methods are limited. Therefore, we model soft tissue damage and failure using a particle/continuum approach. Specifically, we combine continuum damage theory with Smoothed Particle Hydrodynamics (SPH). Because SPH is a meshless particle method, and particle connectivity is determined solely through a neighbor list, discontinuities can be readily modeled by modifying this list. We show, for the first time, that an anisotropic hyperelastic constitutive model commonly employed for modeling soft tissue can be conveniently implemented within a SPH framework and that SPH results show excellent agreement with analytical solutions for uniaxial and biaxial extension as well as finite element solutions for clamped uniaxial extension in 2D and 3D. We further develop a simple algorithm that automatically detects damaged particles and disconnects the spatial domain along rupture lines in 2D and rupture surfaces in 3D. We demonstrate the utility of this approach by simulating damage and failure under clamped uniaxial extension and in a peeling experiment of virtual soft tissue samples. In conclusion, SPH in combination with continuum damage theory may provide an accurate and efficient framework for modeling damage and failure in soft tissues.     

Prediction of fibre architecture and adaptation in diseased carotid bifurcations

Many studies have used patient-specific finite element models to estimate the stress environment in atherosclerotic plaques, attempting to correlate the magnitude of stress to plaque vulnerability. In complex geometries, few studies have incorporated the anisotropic material response of arterial tissue. This paper presents a fibre remodelling algorithm to predict the fibre architecture, and thus anisotropic material response in four patient-specific models of the carotid bifurcation. The change in fibre architecture during disease progression and its affect on the stress environment in the plaque were predicted. The mean fibre directions were assumed to lie at an angle between the two positive principal strain directions. The angle and the degree of dispersion were assumed to depend on the ratio of principal strain values. Results were compared with experimental observations and other numerical studies. In non-branching regions of each model, the typical double helix arterial fibre pattern was predicted while at the bifurcation and in regions of plaque burden, more complex fibre architectures were found. The predicted change in fibre architecture in the arterial tissue during plaque progression was found to alter the stress environment in the plaque. This suggests that the specimen-specific anisotropic response of the tissue should be taken into account to accurately predict stresses in the plaque. Since determination of the fibre architecture in vivo is a difficult task, the system presented here provides a useful method of estimating the fibre architecture in complex arterial geometries.     

Computational modeling of the arterial wall based on layer-specific histological data

Arterial walls typically have a heterogeneous structure with three different layers (intima, media, and adventitia). Each layer can be modeled as a fiber-reinforced material with two families of relatively stiff collagenous fibers symmetrically arranged within an isotropic soft ground matrix. In this paper, we present two different modeling approaches, the embedded fiber (EF) approach and the angular integration (AI) approach, to simulate the anisotropic behavior of individual arterial wall layers involving layer-specific data. The EF approach directly incorporates the microscopic arrangement of fibers that are synthetically generated from a random walk algorithm and captures material anisotropy at the element level of the finite element formulation. The AI approach smears fibers in the ground matrix and treats the material as homogeneous, with material anisotropy introduced at the constitutive level by enhancing the isotropic strain energy with two anisotropic terms. Both approaches include the influence of fiber dispersion introduced by fiber angular distribution (departure of individual fibers from the mean orientation) and take into consideration the dispersion caused by fiber waviness, which has not been previously considered. By comparing the numerical results with the published experimental data of different layers of a human aorta, we show that by using histological data both approaches can successfully capture the anisotropic behavior of individual arterial wall layers. Furthermore, through a comprehensive parametric study, we establish the connections between the AI phenomenological material parameters and the EF parameters having straightforward physical or geometrical interpretations. This study provides valuable insight for the calibration of phenomenological parameters used in the homogenized modeling based on the fiber microscopic arrangement. Moreover, it facilitates a better understanding of individual arterial wall layers, which will eventually advance the study of the structure–function relationship of arterial walls as a whole.     

Quantifying nonlinear anisotropic elastic material properties of biological tissue by use of membrane inflation

Determination of material parameters for soft tissue frequently involves regression of material parameters for nonlinear, anisotropic constitutive models against experimental data from heterogeneous tests. Here, parameter estimation based on membrane inflation is considered. A four parameter nonlinear, anisotropic hyperelastic strain energy function was used to model the material, in which the parameters are cast in terms of key response features. The experiment was simulated using finite element (FE) analysis in order to predict the experimental measurements of pressure versus profile strain. Material parameter regression was automated using inverse FE analysis; parameter values were updated by use of both local and global techniques, and the ability of these techniques to efficiently converge to a best case was examined. This approach provides a framework in which additional experimental data, including surface strain measurements or local structural information, may be incorporated in order to quantify heterogeneous nonlinear material properties.     

A large-strain finite element formulation for biological tissues with application to mitral valve leaflet tissue mechanics

This paper presents a finite element formulation suitable for large-strain modeling of biological tissues and uses this formulation to implement an accurate finite element model for mitral valve leaflet tissue. First, an experimentally derived strain energy function is obtained from literature. This function is implemented in finite elements using the mixed pressure-displacement formulation. A modification is made to aid in maintaining positive definiteness of the stiffness matrix at low strains. The numerical implementation is shown to be accurate in representing the analytical model of material behavior. The mixed formulation is useful for modeling of soft biological tissues in general, and the model presented here is applicable to finite element simulation of mitral valve mechanics.     

A discrete-time approach to the formulation of constitutive models for viscoelastic soft tissues

This paper presents a novel approach to constitutive modeling of viscoelastic soft tissues. This formulation combines an anisotropic strain energy function, accounting for preferred material directions, to define the elastic stress–strain relationship, and a discrete time black-box dynamic model, borrowed from the theory of system identification, to describe the time-dependent behavior. This discrete time formulation is straightforwardly oriented to the development of a recursive time integration scheme that calculates the current stress state by using strain and stress values stored at a limited number of previous time instants. The viscoelastic model and the numerical procedure are assessed by implementing two numerical examples, the simulation of a uniaxial tensile test and the inflation of a thin tube. Both simulations are performed using parameter values based on previous experiments on preserved bovine pericardium. Parameters are then adjusted to investigate the sensitivity of the model. The hypotheses the model relies upon are discussed and the main limitations are stated.     

Non-linear rotation-free shell finite-element models for aortic heart valves

Hyperelastic material models have been incorporated in the rotation-free, large deformation, shell finite element (FE) formulation of (Stolarski et al., 2013) and applied to dynamic simulations of aortic heart valve. Two models used in the past in analysis of such problem i.e. the Saint-Venant and May-Newmann–Yin (MNY) material models have been considered and compared. Uniaxial tests for those constitutive equations were performed to verify the formulation and implementation of the models. The issue of leaflets interactions during the closing of the heart valve at the end of systole is considered. The critical role of using non-linear anisotropic model for proper dynamic response of the heart valve especially during the closing phase is demonstrated quantitatively. This work contributes an efficient FE framework for simulating biological tissues and paves the way for high-fidelity flow structure interaction simulations of native and bioprosthetic aortic heart valves.     

A visco-hyperelastic-damage constitutive model for the analysis of the biomechanical response of the periodontal ligament

The periodontal ligament (PDL), as other soft biological tissues, shows a strongly non-linear and time-dependent mechanical response and can undergo large strains under physiological loads. Therefore, the characterization of the mechanical behavior of soft tissues entails the definition of constitutive models capable of accounting for geometric and material non-linearity. The microstructural arrangement determines specific anisotropic properties. A hyperelastic anisotropic formulation is adopted as the basis for the development of constitutive models for the PDL and properly arranged for investigating the viscous and damage phenomena as well to interpret significant aspects pertaining to ordinary and degenerative conditions. Visco-hyperelastic models are used to analyze the time-dependent mechanical response, while elasto-damage models account for the stiffness and strength decrease that can develop under significant loading or degenerative conditions. Experimental testing points out that damage response is affected by the strain rate associated with loading, showing a decrease in the damage limits as the strain rate increases. These phenomena can be investigated by means of a model capable of accounting for damage phenomena in relation to viscous effects. The visco-hyperelastic-damage model developed is defined on the basis of a Helmholtz free energy function depending on the strain-damage history. In particular, a specific damage criterion is formulated in order to evaluate the influence of the strain rate on damage. The model can be implemented in a general purpose finite element code. The accuracy of the formulation is evaluated by using results of experimental tests performed on animal model, accounting for different strain rates and for strain states capable of inducing damage phenomena. The comparison shows a good agreement between numerical results and experimental data.     

Simulation of soft tissue failure using the material point method

Understanding the factors that control the extent of tissue damage as a result of material failure in soft tissues may provide means to improve diagnosis and treatment of soft tissue injuries. The objective of this research was to develop and test a computational framework for the study of the failure of anisotropic soft tissues subjected to finite deformation. An anisotropic constitutive model incorporating strain-based failure criteria was implemented in an existing computational solid mechanics software based on the material point method (MPM), a quasi-meshless particle method for simulations in computational mechanics. The constitutive model and the strain-based failure formulations were tested using simulations of simple shear and tensile mechanical tests. The model was then applied to investigate a scenario of a penetrating injury: a low-speed projectile was released through a myocardial material slab. Sensitivity studies were performed to establish the necessary grid resolution and time-step size. Results of the simple shear and tensile test simulations demonstrated the correct implementation of the constitutive model and the influence of both fiber family and matrix failure on predictions of overall tissue failure. The slab penetration simulations produced physically realistic wound tracts, exhibiting diameter increase from entrance to exit. Simulations examining the effect of bullet initial velocity showed that the anisotropy influenced the shape and size of the exit wound more at lower velocities. Furthermore, the size and taper of the wound cavity was smaller for the higher bullet velocity. It was concluded that these effects were due to the amount of momentum transfer. The results demonstrate the feasibility of using MPM and the associated failure model for large-scale numerical simulations of soft tissue failure.     

An anisotropic elastic-viscoplastic damage model for bone tissue

A new anisotropic elastic-viscoplastic damage constitutive model for bone is proposed using an eccentric elliptical yield criterion and nonlinear isotropic hardening. A micromechanics-based multiscale homogenization scheme proposed by Reisinger et al. is used to obtain the effective elastic properties of lamellar bone. The dissipative process in bone is modeled as viscoplastic deformation coupled to damage. The model is based on an orthotropic ecuntric elliptical criterion in stress space. In order to simplify material identification, an eccentric elliptical isotropic yield surface was defined in strain space, which is transformed to a stress-based criterion by means of the damaged compliance tensor. Viscoplasticity is implemented by means of the continuous Perzyna formulation. Damage is modeled by a scalar function of the accumulated plastic strain ${D(\kappa)}$ , reducing all element s of the stiffness matrix. A polynomial flow rule is proposed in order to capture the rate-dependent post-yield behavior of lamellar bone. A numerical algorithm to perform the back projection on the rate-dependent yield surface has been developed and implemented in the commercial finite element solver Abaqus/Standard as a user subroutine UMAT. A consistent tangent operator has been derived and implemented in order to ensure quadratic convergence. Correct implementation of the algorithm, convergence, and accuracy of the tangent operator was tested by means of strain- and stress-based single element tests. A finite element simulation of nano- indentation in lamellar bone was finally performed in order to show the abilities of the newly developed constitutive model.