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 influence of the material properties on the biomechanical behavior of the pelvic floor muscles during vaginal delivery

In this work, a finite element model intends to represent the effects that the passage of a fetal head can induce on the muscles of the pelvic floor, from a mechanical point of view.The finite element method is a valuable tool, that is contributing to the clarification of the mechanisms behind pelvic floor disorders related to vaginal deliveries, although some care is necessary in order to obtain correct results. The present work shows how the variation of the material parameters, used in the constitutive model, can affect the obtained results from a finite element simulation. The constitutive equation adopted in this work for the pelvic floor muscles is a modified form of the incompressible transversely isotropic hyperelastic model proposed earlier by Humphrey and Yin.Results for the pelvic floor strain and stresses obtained during the passage of the fetus head are presented. The results show the importance of the material parameters and the need for a correct constitutive model.     

Toward a biomechanical tool to evaluate rupture potential of abdominal aortic aneurysm: identification of a finite strain constitutive model and evaluation of its applicability

Knowledge of the wall stresses in an abdominal aortic aneurysm (AAA) may be helpful in evaluating the need for surgical intervention to avoid rupture. This must be preceded by the development of a more suitable finite strain constitutive model for AAA, as none currently exists. Additionally, reliable stress analysis of in vivo AAA for the purposes of clinical diagnostics requires patient-specific values of the material parameters, which are difficult to determine noninvasively. The purpose of this work, therefore, was three-fold: (1) to develop a finite strain constitutive model for AAA; (2) to estimate the variation of model parameters within a sample population; and (3) to evaluate the sensitivity of computed stress distribution in AAA due to this biologic variation. We propose here a two parameter, hyperelastic, isotropic, incompressible material model and utilize experimental data from 69 freshly excised AAA specimens to both develop the functional form of the model and estimate its material parameters. Parametric analyses were performed via repeated finite element computations to determine the effect of varying each of the two model parameters on the stress distribution in a three-dimensional AAA model. The agreement between experimental data and the proposed functional form of the constitutive law was very good (R2 > 0.9). Our finite element simulations showed that the computed AAA wall stresses changed by only 4% or less when both the parameters were varied within the 95% confidence intervals for the patient population studied. This observation indicates that in lieu of the patient-specific material parameters, which are difficult to determine the use of population mean values is sufficiently accurate for the model to be reasonably employed in a clinical setting. We believe that this is an important advancement toward the development of a computational tool for the estimation of rupture potential for individual AAA, for which there is great clinical need.     

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.     

The non-linear response of a muscle in transverse compression: assessment of geometry influence using a finite element model

Most recent finite element models that represent muscles are generic or subject-specific models that use complex, constitutive laws. Identification of the parameters of such complex, constitutive laws could be an important limit for subject-specific approaches. The aim of this study was to assess the possibility of modelling muscle behaviour in compression with a parametric model and a simple, constitutive law. A quasi-static compression test was performed on the muscles of dogs. A parametric finite element model was designed using a linear, elastic, constitutive law. A multi-variate analysis was performed to assess the effects of geometry on muscle response. An inverse method was used to define Young's modulus. The non-linear response of the muscles was obtained using a subject-specific geometry and a linear elastic law. Thus, a simple muscle model can be used to have a bio-faithful, biomechanical response.     

Fitted hyperelastic parameters for Human brain tissue from reported tension,compression, and shear tests

The mechanical properties of human brain tissue are the subject of interest because of their use in understanding brain trauma and in developing therapeutic treatments and procedures. To represent the behavior of the tissue, we have developed hyperelastic mechanical models whose parameters are fitted in accordance with experimental test results. However, most studies available in the literature have fitted parameters with data of a single type of loading, such as tension, compression, or shear. Recently, Jin et al. (Journal of Biomechanics 46:2795−2801, 2013) reported data from ex vivo tests of human brain tissue under tension, compression, and shear loading using four strain rates and four different brain regions. However, they do not report parameters of energy functions that can be readily used in finite element simulations. To represent the tissue behavior for the quasi-static loading conditions, we aimed to determine the best fit of the hyperelastic parameters of the hyperfoam, Ogden, and polynomial strain energy functions available in ABAQUS for the low strain rate data, while simultaneously considering all three loading modes. We used an optimization process conducted in MATLAB, calling iteratively three finite element models developed in ABAQUS that represent the three loadings. Results showed a relatively good fit to experimental data in all loading modes using two terms in the energy functions. Values for the shear modulus obtained in this analysis (897−1653 Pa) are in the range of those presented in other studies. These energy-function parameters can be used in brain tissue simulations using finite element models.     

Characterisation of Anisotropic and Non-linear Behaviour of Human Skin In Vivo

The aim of this work is to characterise the in-plane mechanical behaviour of human skin in vivo. For this purpose the structural skin model proposed by Lanir [1] is employed and a mixed numerical-experimental method is developed. The numerical-experimental method is based on the confrontation of measured data from an experiment, with calculated data from a finite element model, eventually leading to the determination of some of the parameters of a constitutive model, in the present case Lanir's Skin Model. Since collagen, the main constituent of skin, dominates the anisotropic and non-linear behaviour of skin, the parameters of Lanir's Skin Model concerning the mechanical behaviour of the collagen fibres are estimated. In vivo experiments were carried out on the volar forearm. During the experiments, reaction forces and the displacement field at different states of deformation are measured. Both data sets are used for the determination of the parameters.     

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.     

On modelling large deformations of heterogeneous biological tissues using a mixed finite element formulation

This study addresses the issue of modelling material heterogeneity of incompressible bodies. It is seen that when using a mixed (displacement–pressure) finite element formulation, the basis functions used for pressure field may not be able to capture the nonlinearity of material parameters, resulting in pseudo-residual stresses. This problem can be resolved by modifying the constitutive relation using Flory's decomposition of the deformation gradient. A two-parameter Mooney–Rivlin constitutive relation is used to demonstrate the methodology. It is shown that for incompressible materials, the modification does not alter the mechanical behaviour described by the original constitutive model. In fact, the modified constitutive equation shows a better predictability when compared against analytical solutions. Two strategies of describing the material variation (i.e. linear and step change) are explained, and their solutions are evaluated for an ideal two-material interfacing problem. When compared with the standard tied coupling approach, the step change method exhibited a much better agreement because of its ability to capture abrupt changes of the material properties. The modified equation in conjunction with integration point-based material heterogeneity is then used to simulate the deformations of heterogeneous biological structures to illustrate its applications.     

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

Modelling and simulation of porcine liver tissue indentation using finite element method and uniaxial stress–strain data

We hypothesize that both compression and elongation stress–strain data should be considered for modeling and simulation of soft tissue indentation. Uniaxial stress–strain data were obtained from in vitro loading experiments of porcine liver tissue. An axisymmetric finite element model was used to simulate liver tissue indentation with tissue material represented by hyperelastic models. The material parameters were derived from uniaxial stress–strain data of compressions, elongations, and combined compression and elongation of porcine liver samples. in vitro indentation tests were used to validate the finite element simulation. Stress–strain data from the simulation with material parameters derived from the combined compression and elongation data match the experimental data best. This is due to its better ability in modeling 3D deformation since the behavior of biological soft tissue under indentation is affected by both its compressive and tensile characteristics. The combined logarithmic and polynomial model is somewhat better than the 5-constant Mooney–Rivlin model as the constitutive model for this indentation simulation.     

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 transversely isotropic constitutive model of excised guinea pig spinal cord white matter

Narrowing of the spinal canal generates an amalgamation of stresses within the spinal cord parenchyma. The tissue’s stress state cannot be quantified experimentally; it must be described using computational methods, such as finite element analysis. The objective of this research was to propose a compressible, transversely isotropic constitutive model, an augmentation of the isotropic Mooney–Rivlin hyperelastic strain energy function, to describe the guinea pig spinal cord white matter. Model parameters were derived from a combination of inverse finite element analysis on transverse compression experiments and least squared error analysis applied to quasi-static longitudinal tensile tests. A comparison of the residual errors between the predicted response and the experimental measurements indicated that the transversely isotropic constitutive law that incorporates an offset stretch reduced the error by a factor of four when compared to other commonly used models.     

Predicting the failure response of cement-bone constructs using a non-linear fracture mechanics approach

A non-linear fracture mechanics approach was used to predict the failure response of complex cement-bone constructs. A series of eight mechanical tests with a combination of tensile and shear loading along the cement-bone interface was performed. Each experiment was modeled using the finite element method with non-linear constitutive models at the cement-bone interface. Interface constitutive parameters were assigned based on the quantity of bone interdigitated with the cement. There was a strong correlation (r2 = 0.80) between experimentally measured and finite element predicted ultimate loads. The average error in predicted ultimate load was 23.9 percent. In comparison to the ultimate load predictions, correlations and errors for total energy to failure (r2 = 0.24, avg. error = 38.2 percent) and displacement at 50 percent of the ultimate load (r2 = 0.27, avg. error = 52.2 percent) were poor The results indicate that the non-linear constitutive laws could be useful in predicting the initiation and progression of interface failure of cemented bone-implant systems. However improvements in the estimation of post-yield interface properties from the quantity of bone interdigitated with cement are needed to enhance predictions of the overall failure response.     

Subject-specific non-linear biomechanical model of needle insertion into brain

The previous models for predicting the forces acting on a needle during insertion into very soft organs (such as, e.g. brain) relied on oversimplifying assumptions of linear elasticity and specific experimentally derived functions for determining needle-tissue interactions. In this contribution, we propose a more general approach in which the needle forces are determined directly from the equations of continuum mechanics using fully non-linear finite element procedures that account for large deformations (geometric non-linearity) and non-linear stress-strain relationship (material non-linearity) of soft tissues. We applied these procedures to model needle insertion into a swine brain using the constitutive properties determined from the experiments on tissue samples obtained from the same brain (i.e. the subject-specific constitutive properties were used). We focused on the insertion phase preceding puncture of the brain meninges and obtained a very accurate prediction of the needle force. This demonstrates the utility of non-linear finite element procedures in patient-specific modelling of needle insertion into soft organs such as, e.g. brain.     

Leukocyte deformability: finite element modeling of large viscoelastic deformation.

An axisymmetric deformation of a viscoelastic sphere bounded by a prestressed elastic thin shell in response to external pressure is studied by a finite element method. The research is motivated by the need for understanding the passive behavior of human leukocytes (white blood cells) and interpreting extensive experimental data in terms of the mechanical properties. The cell at rest is modeled as a sphere consisting of a cortical prestressed shell with incompressible Maxwell fluid interior. A large-strain deformation theory is developed based on the proposed model. General non-linear, large strain constitutive relations for the cortical shell are derived by neglecting the bending stiffness. A representation of the constitutive equations in the form of an integral of strain history for the incompressible Maxwell interior is used in the formulation of numerical scheme. A finite element program is developed, in which a sliding boundary condition is imposed on all contact surfaces. The mathematical model developed is applied to evaluate experimental data of pipette tests and observations of blood flow.