Constitutive model for brain tissue under finite compression

While advances in computational models of mechanical phenomena have made it possible to simulate dynamically complex problems in biomechanics, accurate material models for soft tissues, particularly brain tissue, have proven to be very challenging. Most studies in the literature on material properties of brain tissue are performed in shear loading and very few tackle the behavior of brain in compression. In this study, a viscoelastic constitutive model of bovine brain tissue under finite step-and-hold uniaxial compression with 10 s(-1) ramp rate and 20 s hold time has been developed. The assumption of quasi-linear viscoelasticity (QLV) was validated for strain levels of up to 35%. A generalized Rivlin model was used for the isochoric part of the deformation and it was shown that at least three terms (C(10), C(01) and C(11)) are needed to accurately capture the material behavior. Furthermore, for the volumetric deformation, a two parameter Ogden model was used and the extent of material incompressibility was studied. The hyperelastic material parameters were determined through extracting and fitting to two isochronous curves (0.06 s and 14 s) approximating the instantaneous and steady-state elastic responses. Viscoelastic relaxation was characterized at five decay rates (100, 10, 1, 0.1, 0 s(-1)) and the results in compression and their extrapolation to tension were compared against previous models.     

Constitutive modelling of brain tissue: Experiment and theory

Recent developments in computer-integrated and robot-aided surgery—in particular, the emergence of automatic surgical tools and robots—as well as advances in virtual reality techniques, call for closer examination of the mechanical properties of very soft tissues (such as brain, liver, kidney, etc.). The ultimate goal of our research into the biomechanics of these tissues is the development of corresponding, realistic mathematical models. This paper contains experimental results of in vitro, uniaxial, unconfined compression of swine brain tissue and discusses a single-phase, non-linear, viscoelastic tissue model. The experimental results obtained for three loading velocities, ranging over five orders of magnitude, are presented. The applied strain rates have been much lower than those applied in previous studies, focused on injury modelling. The stress-strain curves are concave upward for all compression rates containing no linear portion from which a meaningful elastic modulus might be determined. The tissue response stiffened as the loading speed increased, indicating a strong stress-strain rate dependence. The use of the single-phase model is recommended for applications in registration, surgical operation planning and training systems as well as a control system of an image-guided surgical robot. The material constants for the brain tissue are evaluated. Agreement between the proposed theoretical model and experiment is good for compression levels reaching 30% and for loading velocities varying over five orders of magnitude.     

Inertia effects on characterization of dynamic response of brain tissue

Modeling and simulation of traumatic brain injury (TBI) resulted from collision or blast loading requires characterization of mechanical response over a wide range of loading rates under valid testing conditions. In this study, mechanical response of fresh bovine brain tissue was studied using the two modified Kolsky bar techniques. Radial deformation behavior of annular specimens, which are typically used to characterize the dynamic uniaxial compressive response of biological tissues, was examined using a modified Kolsky bar and a high speed camera to collect images while the specimen deforms at an axial strain rate of 2000s(-1). The high-speed images revealed inhomogeneous specimen deformation possibly brought about by radial inertia and causing a multi-axial stress state. To acquire valid stress-strain results that can be used to produce constitutive behavior of the soft materials, a novel torsion technique was developed to obtain pure shear response at dynamic loading rates. Experimental results show clear differences in the material response using the two methods. These results indicate that the previously demonstrated annular specimen geometry aimed at reducing inertia induced stress components for high rate soft materials uniaxial-compressive testing may still possess a significant component of radial inertia induced radial stress which consequently caused the observed inhomogeneous deformation in brain tissue test samples.     

The mechanical behaviour of brain tissue: large strain response and constitutive modelling

The non-linear mechanical behaviour of porcine brain tissue in large shear deformations is determined. An improved method for rotational shear experiments is used, producing an approximately homogeneous strain field and leading to an enhanced accuracy. Results from oscillatory shear experiments with a strain amplitude of 0.01 and frequencies ranging from 0.04 to 16 Hz are given. The immediate loss of structural integrity, due to large deformations, influencing the mechanical behaviour of brain tissue, at the time scale of loading, is investigated. No significant immediate mechanical damage is observed for these shear deformations up to strains of 0.45. Moreover, the material behaviour during complex loading histories (loading-unloading) is investigated. Stress relaxation experiments for strains up to 0.2 and constant strain rate experiments for shear rates from 0.01 to 1 s(-1) and strains up to 0.15 are presented. A new differential viscoelastic model is used to describe the mechanical response of brain tissue. The model is formulated in terms of a large strain viscoelastic framework and considers non-linear viscous deformations in combination with non-linear elastic behaviour. This constitutive model is readily applicable in three-dimensional head models in order to predict the mechanical response of the intra-cranial contents due to an impact.     

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.     

Dynamic response of immature bovine articular cartilage in tension and compression, and nonlinear viscoelastic modeling of the tensile response

Very limited information is currently available on the constitutive modeling of the tensile response of articular cartilage and its dynamic modulus at various loading frequencies. The objectives of this study were to (1) formulate and experimentally validate a constitutive model for the intrinsic viscoelasticity of cartilage in tension, (2) confirm the hypothesis that energy dissipation in tension is less than in compression at various loading frequencies, and (3) test the hypothesis that the dynamic modulus of cartilage in unconfined compression is dependent upon the dynamic tensile modulus. Experiment 1: Immature bovine articular cartilage samples were tested in tensile stress relaxation and cyclical loading. A proposed reduced relaxation function was fitted to the stress-relaxation response and the resulting material coefficients were used to predict the response to cyclical loading. Adjoining tissue samples were tested in unconfined compression stress relaxation and cyclical loading. Experiment 2: Tensile stress relaxation experiments were performed at varying strains to explore the strain-dependence of the viscoelastic response. The proposed relaxation function successfully fit the experimental tensile stress-relaxation response, with R2 = 0.970+/-0.019 at 1% strain and R2 = 0.992+/-0.007 at 2% strain. The predicted cyclical response agreed well with experimental measurements, particularly for the dynamic modulus at various frequencies. The relaxation function, measured from 2% to 10% strain, was found to be strain dependent, indicating that cartilage is nonlinearly viscoelastic in tension. Under dynamic loading, the tensile modulus at 10 Hz was approximately 2.3 times the value of the equilibrium modulus. In contrast, the dynamic stiffening ratio in unconfined compression was approximately 24. The energy dissipation in tension was found to be significantly smaller than in compression (dynamic phase angle of 16.7+/-7.4 deg versus 53.5+/-12.8 deg at 10(-3) Hz). A very strong linear correlation was observed between the dynamic tensile and dynamic compressive moduli at various frequencies (R2 = 0.908+/-0.100). The tensile response of cartilage is nonlinearly viscoelastic, with the relaxation response varying with strain. A proposed constitutive relation for the tensile response was successfully validated. The frequency response of the tensile modulus of cartilage was reported for the first time. Results emphasize that fluid-flow dependent viscoelasticity dominates the compressive response of cartilage, whereas intrinsic solid matrix viscoelasticity dominates the tensile response. Yet the dynamic compressive modulus of cartilage is critically dependent upon elevated values of the dynamic tensile modulus.     

Nonlinear mechanical property of tracheal cartilage: a theoretical and experimental study

BACKGROUND: Despite being the stiffest airway of the bronchial tree, the trachea undergoes significant deformation due to intrathoracic pressure during breathing. The mechanical properties of the trachea affect the flow in the airway and may contribute to the biological function of the lung. METHOD: A Fung-type strain energy density function was used to investigate the nonlinear mechanical behavior of tracheal cartilage. A bending test on pig tracheal cartilage was performed and a mathematical model for analyzing the deformation of tracheal cartilage was developed. The constants included in the strain energy density function were determined by fitting the experimental data. RESULT: The experimental data show that tracheal cartilage is a nonlinear material displaying higher strength in compression than in tension. When the compression forces varied from -0.02 to -0.03N and from -0.03 to -0.04N, the deformation ratios were 11.03+/-2.18% and 7.27+/-1.59%, respectively. Both were much smaller than the deformation ratios (20.01+/-4.49%) under tension forces of 0.02 to 0.01N. The Fung-type strain energy density function can capture this nonlinear behavior very well, whilst the linear stress-strain relation cannot. It underestimates the stability of trachea by exaggerating the displacement in compression. This study may improve our understanding of the nonlinear behavior of tracheal cartilage and it may be useful for the future study on tracheal collapse behavior under physiological and pathological conditions.     

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.     

Material properties of the placenta under dynamic loading conditions

Trauma during pregnancy especially occurring during car crashes leads to many foetal losses. Numerical modelling is widely used in car occupant safety issue and injury mechanisms analysis and is particularly adapted to the pregnant woman. Material modelling of the gravid uterus tissues is crucial for injury risk evaluation especially for the abruption placentae which is widely assumed as the leading cause of foetal loss. Experimental studies on placenta behaviour in tension are reported in the literature, but none in compression to the authors' knowledge. This lack of data is addressed in this study. To complement the already available experimental literature data on the placenta mechanical behaviour and characterise it in a compression loading condition, 80 indentation tests on fresh placentae are presented. Hyperelastic like mean experimental stress versus strain and corridors are exposed. The results of the experimental placenta indentations compared with the tensile literature results tend to show a quasi-symmetrical behaviour of the tissue. An inverse analysis using simple finite element models has permitted to propose parameters for an Ogden material model for the placenta which exhibits a realistic behaviour in both tension and compression.     

Failure modelling of trabecular bone using a non-linear combined damage and fracture voxel finite element approach

Trabecular bone tissue failure can be considered as consisting of two stages: damage and fracture; however, most failure analyses of 3D high-resolution trabecular bone samples are confined to damage mechanisms only, that is, without fracture. This study aims to develop a computational model of trabecular bone consisting of an explicit representation of complete failure, incorporating damage criteria, fracture criteria, cohesive forces, asymmetry and large deformation capabilities. Following parameter studies on a test specimen, and experimental testing of bone sample to complete failure, the asymmetric critical tissue damage and fracture strains of ovine vertebral trabecular bone were calibrated and validated to be compression damage ?1.16 %, tension damage 0.69 %, compression fracture ?2.91 % and tension fracture 1.98 %. Ultimate strength and post–ultimate strength softening were captured by the computational model, and the failure of individual struts in bending and shear was also predicted. This modelling approach incorporated a cohesive parameter that provided a facility to calibrate ductile–brittle behaviour of bone tissue in this non-linear geometric and non-linear constitutive property analyses tool. Finally, the full accumulation of tissue damage and tissue fracture has been monitored from range of small magnitude (normal daily loading) through to specimen yielding, ultimate strength and post–ultimate strength softening.     

Constitutive behavior of ocular tissues over a range of strain rates

The constitutive behavior of bovine scleral and corneal tissues is measured in tension and compression, at quasi-static and moderate strain rates. Experiments are conducted at strain rates up to about 50 strain per second by a pneumatic testing system developed to overcome noise and measurement difficulties associated with the time dependent test of low impedance materials. Results for the tissues at room and the natural bovine body temperatures are similar and indicate that ocular tissue exhibits nonlinear stiffening for increasing strain rates, a phenomena termed rate hardening. For example, at a tensile strain rate of 29/s, corneal tissue is found to develop 10 times the stress that it does quasi-statically at the same strain. Thus, conventional constitutive models will grossly underpredict stresses occurring in the corneo-scleral shell due to moderate dynamic events. This has implication to the accuracy of ocular injury models, the study of the stress field in the corneo-scleral shell for glaucoma research and tonometry measurements. The measured data at various strain rates is represented using the general framework of a constitutive model that has been used to represent biological tissue mechanical data. Here it is extended to represent the measured data of the ocular tissues over the range of tested strain rates. Its form allows for straightforward incorporation in various numerical codes. The experimental and analytical methods developed here are felt to be applicable to the test of human ocular tissue.     

Large strain behaviour of brain tissue in shear: some experimental data and differential constitutive model.

In this paper, some experimental measurements of the behaviour of bovine brain tissue under large shear strains in vitro are reported, and a constitutive model which is consistent with the data is developed. It was determined that brain tissue is not strain-time separable, showing slower relaxation at higher strains, and that the stresses in shear are not linear with increasing shear strain. The new constitutive model is a differential model, including both an "elastic" term, of the Mooney type and a nonlinear viscoelastic term. The latter allows for the change in relaxation behaviour with strain, by modifying an upper convected multimode Maxwell model with a damping function. The model shows good agreement with the experimental shear results and could be used to describe other types of data.     

Constitutive model of brain tissue suitable for finite element analysis of surgical procedures.

Realistic finite element modelling and simulation of neurosurgical procedures present a formidable challenge. Appropriate, finite deformation, constitutive model of brain tissue is a prerequisite for such development. In this paper, a large deformation, linear, viscoelastic model, suitable for direct use with commercially available finite element software packages such as ABAQUS is constructed. The proposed constitutive equation is of polynomial form with time-dependent coefficients. The model requires four material constants to be identified. The material constants were evaluated based on unconfined compression experiment results. The analytical as well as numerical solutions to the unconfined compression problem are presented. The agreement between the proposed theoretical model and the experiment is good for compression levels reaching 30% and for loading velocities varying over five orders of magnitude. The numerical solution using the finite element method matched the analytical solution very closely.     

Experimental verification of the roles of intrinsic matrix viscoelasticity and tension-compression nonlinearity in the biphasic response of cartilage

A biphasic-CLE-QLV model proposed in our recent study [2001, J. Biomech. Eng., 123, pp. 410-417] extended the biphasic theory of Mow et al. [1980, J. Biomech. Eng., 102, pp. 73-84] to include both tension-compression nonlinearity and intrinsic viscoelasticity of the cartilage solid matrix by incorporating it with the conewise linear elasticity (CLE) model [1995, J. Elasticity, 37, pp. 1-38] and the quasi-linear viscoelasticity (QLV) model [Biomechanics: Its foundations and objectives, Prentice Hall, Englewood Cliffs, 1972]. This model demonstrates that a simultaneous prediction of compression and tension experiments of articular cartilage, under stress-relaxation and dynamic loading, can be achieved when properly taking into account both flow-dependent and flow-independent viscoelastic effects, as well as tension-compression nonlinearity. The objective of this study is to directly test this biphasic-CLE-QLV model against experimental data from unconfined compression stress-relaxation tests at slow and fast strain rates as well as dynamic loading. Twelve full-thickness cartilage cylindrical plugs were harvested from six bovine glenohumeral joints and multiple confined and unconfined compression stress-relaxation tests were performed on each specimen. The material properties of specimens were determined by curve-fitting the experimental results from the confined and unconfined compression stress relaxation tests. The findings of this study demonstrate that the biphasic-CLE-QLV model is able to describe the strain-rate-dependent mechanical behaviors of articular cartilage in unconfined compression as attested by good agreements between experimental and theoretical curvefits (r2 = 0.966 +/- 0.032 for testing at slow strain rate; r2 = 0.998 +/- 0.002 for testing at fast strain rate) and predictions of the dynamic response (r2 = 0.91 +/- 0.06). This experimental study also provides supporting evidence for the hypothesis that both tension-compression nonlinearity and intrinsic viscoelasticity of the solid matrix of cartilage are necessary for modeling the transient and equilibrium responses of this tissue in tension and compression. Furthermore, the biphasic-CLE-QLV model can produce better predictions of the dynamic modulus of cartilage in unconfined dynamic compression than the biphasic-CLE and biphasic poroviscoelastic models, indicating that intrinsic viscoelasticity and tension-compression nonlinearity of articular cartilage may play important roles in the load-support mechanism of cartilage under physiologic loading.     

High-resolution finite element models with tissue strength asymmetry accurately predict failure of trabecular bone

The ability to predict trabecular failure using microstructure-based computational models would greatly facilitate study of trabecular structure–function relations, multiaxial strength, and tissue remodeling. We hypothesized that high-resolution finite element models of trabecular bone that include cortical-like strength asymmetry at the tissue level, could predict apparent level failure of trabecular bone for multiple loading modes. A bilinear constitutive model with asymmetric tissue yield strains in tension and compression was applied to simulate failure in high-resolution finite element models of seven bovine tibial specimens. Tissue modulus was reduced by 95% when tissue principal strains exceeded the tissue yield strains. Linear models were first calibrated for effective tissue modulus against specimen-specific experimental measures of apparent modulus, producing effective tissue moduli of (mean±S.D.) 18.7±3.4 GPa. Next, a parameter study was performed on a single specimen to estimate the tissue level tensile and compressive yield strains. These values, 0.60% strain in tension and 1.01% strain in compression, were then used in non-linear analyses of all seven specimens to predict failure for apparent tensile, compressive, and shear loading. When compared to apparent yield properties previously measured for the same type of bone, the model predictions of both the stresses and strains at failure were not statistically different for any loading case (p>0.15). Use of symmetric tissue strengths could not match the experimental data. These findings establish that, once effective tissue modulus is calibrated and uniform but asymmetric tissue failure strains are used, the resulting models can capture the apparent strength behavior to an outstanding level of accuracy. As such, these computational models have reached a level of fidelity that qualifies them as surrogates for destructive mechanical testing of real specimens.     

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

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

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

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