A strain-hardening bi-power law for the nonlinear behaviour of biological soft tissues

Biological soft tissues exhibit a strongly nonlinear viscoelastic behaviour. Among parenchymous tissues, kidney and liver remain less studied than brain, and a first goal of this study is to report additional material properties of kidney and liver tissues in oscillatory shear and constant shear rate tests. Results show that the liver tissue is more compliant but more strain hardening than kidney. A wealth of multi-parameter mathematical models has been proposed for describing the mechanical behaviour of soft tissues. A second purpose of this work is to develop a new constitutive law capable of predicting our experimental data in the both linear and nonlinear viscoelastic regime with as few parameters as possible. We propose a nonlinear strain-hardening fractional derivative model in which six parameters allow fitting the viscoelastic behaviour of kidney and liver tissues for strains ranging from 0.01 to 1 and strain rates from 0.0151 s^?1 to 0.7 s^?1.     

Viscoelastic properties of passive skeletal muscle in compression—Cyclic behaviour

Skeletal muscle relaxation behaviour in compression has been previously reported, but the anisotropic behaviour at higher loading rates remains poorly understood. In this paper, uniaxial unconfined cyclic compression tests were performed on fresh porcine muscle samples at various fibre orientations to determine muscle viscoelastic behaviour. Mean compression level of 25% was applied and cycles of 2% and 10% amplitude were performed at 0.2–80 Hz. Under cycles of low frequency and amplitude, linear viscoelastic cyclic relaxation was observed. Fibre/cross-fibre results were qualitatively similar, but the cross-fibre direction was stiffer (ratio of 1.2). In higher amplitude tests nonlinear viscoelastic behaviour with a frequency dependent increase in the stress cycles amplitude was found (factor of 4.1 from 0.2 to 80 Hz).The predictive capability of an anisotropic quasi-linear viscoelastic model previously fitted to stress-relaxation data from similar tissue samples was investigated. Good qualitative results were obtained for low amplitude cycles but differences were observed in the stress cycle amplitudes (errors of 7.5% and 31.8%, respectively, in the fibre/cross-fibre directions). At higher amplitudes significant qualitative and quantitative differences were evident. A nonlinear model formulation was therefore developed which provided a good fit and predictions to high amplitude low frequency cyclic tests performed in the fibre/cross-fibre directions. However, this model gave a poorer fit to high frequency cyclic tests and to relaxation tests. Neither model adequately predicts the stiffness increase observed at frequencies above 5 Hz.Together with data previously presented, the experimental data presented here provide a unique dataset for validation of future constitutive models for skeletal muscle in compression.     

A comparison of cartilage stress-relaxation models in unconfined compression: QLV and stretched exponential in combination with fluid flow

Cartilage exhibits nonlinear viscoelastic behaviour. Various models have been proposed to explain cartilage stress relaxation, but it is unclear whether explicit modelling of fluid flow in unconfined compression is needed. This study compared Fung's quasi-linear viscoelastic (QLV) model with a stretched-exponential model of cartilage stress relaxation and examined each of these models both alone and in combination with a fluid-flow model in unconfined compression. Cartilage explants were harvested from bovine calf patellofemoral joints and equilibrated in tissue culture for 5 days before stress-relaxation testing in unconfined compression at 5% nominal strain. The stretched exponential models fit as well as the QLV models. Furthermore, the average stretched exponential relaxation time determined by this model lies within the range of experimentally measured relaxation times for extracted proteoglycan aggregates, consistent with the hypothesis that the stretched exponential model represents polymeric mechanisms of cartilage viscoelasticity.     

An extended biphasic model for charged hydrated tissues with application to the intervertebral disc

Finite element models for hydrated soft biological tissue are numerous but often exhibit certain essential deficiencies concerning the reproduction of relevant mechanical and electro-chemical responses. As a matter of fact, singlephasic models can never predict the interstitial fluid flow or related effects like osmosis. Quite a few models have more than one constituent, but are often restricted to the small-strain domain, are not capable of capturing the intrinsic viscoelasticity of the solid skeleton, or do not account for a collagen fibre reinforcement. It is the goal of this contribution to overcome these drawbacks and to present a thermodynamically consistent model, which is formulated in a very general way in order to reproduce the behaviour of almost any charged hydrated tissue. Herein, the Theory of Porous Media (TPM) is applied in combination with polyconvex Ogden-type material laws describing the anisotropic and intrinsically viscoelastic behaviour of the solid matrix on the basis of a generalised Maxwell model. Moreover, other features like the deformation-dependent permeability, the possibility to include inhomogeneities like varying fibre alignment and behaviour, or osmotic effects based on the simplifying assumption of Lanir are also included. Finally, the human intervertebral disc is chosen as a representative for complex soft biological tissue behaviour. In this regard, two numerical examples will be presented with focus on the viscoelastic and osmotic capacity of the model.     

Large deformation shear properties of liver tissue

Characterization of the mechanical properties of soft biological tissues is important for establishing the mechanical tolerances of the tissues, and for input to computational models. In this work, the viscoelastic properties of bovine liver tissue in shear loading have been measured using relaxation and constant shear rate loading. The tissue is nonlinearly viscoelastic for strains greater than 0.2%, has a yield strain of approximately 10, and shows moderate strain-rate sensitivity. The response can be modelled using a nonlinear viscoelastic differential model previously developed for brain tissue.     

Nonlinear viscoelastic characterization of bovine trabecular bone

The time-independent elastic properties of trabecular bone have been extensively investigated, and several stiffness–density relations have been proposed. Although it is recognized that trabecular bone exhibits time-dependent mechanical behaviour, a property of viscoelastic materials, the characterization of this behaviour has received limited attention. The objective of the present study was to investigate the time-dependent behaviour of bovine trabecular bone through a series of compressive creep–recovery experiments and to identify its nonlinear constitutive viscoelastic material parameters. Uniaxial compressive creep and recovery experiments at multiple loads were performed on cylindrical bovine trabecular bone samples (\(n = 19\)). Creep response was found to be significant and always comprised of recoverable and irrecoverable strains, even at low stress/strain levels. This response was also found to vary nonlinearly with applied stress. A systematic methodology was developed to separate recoverable (nonlinear viscoelastic) and irrecoverable (permanent) strains from the total experimental strain response. We found that Schapery’s nonlinear viscoelastic constitutive model describes the viscoelastic response of the trabecular bone, and parameters associated with this model were estimated from the multiple load creep–recovery (MLCR) experiments. Nonlinear viscoelastic recovery compliance was found to have a decreasing and then increasing trend with increasing stress level, indicating possible stiffening and softening behaviour of trabecular bone due to creep. The obtained parameters from MLCR tests, expressed as second-order polynomial functions of stress, showed a similar trend for all the samples, and also demonstrate stiffening–softening behaviour with increasing stress.     

On the viscoelastic character of liver tissue: experiments and modelling of the linear behaviour

This paper presents the results of recent investigations into the viscoelastic character of liver tissue, especially the linear response. The linear domain of liver tissue has been determined by a strain sweep oscillation test. The fluid-like viscoelastic behaviour of liver tissue was confirmed by stress relaxation experiments. To demonstrate the consistency and quality of the experimental data, the relaxation spectrum was used to reconstruct the relaxation modulus, as well as the storage and loss moduli. The generalised Maxwell model was successfully employed to model the mechanical behaviour in the linear viscoelastic domain.     

A finite nonlinear hyper-viscoelastic model for soft biological tissues

Soft tissues exhibit highly nonlinear rate and time-dependent stress-strain behaviour. Strain and strain rate dependencies are often modelled using a hyperelastic model and a discrete (standard linear solid) or continuous spectrum (quasi-linear) viscoelastic model, respectively. However, these models are unable to properly capture the materials characteristics because hyperelastic models are unsuited for time-dependent events, whereas the common viscoelastic models are insufficient for the nonlinear and finite strain viscoelastic tissue responses. The convolution integral based models can demonstrate a finite viscoelastic response; however, their derivations are not consistent with the laws of thermodynamics. The aim of this work was to develop a three-dimensional finite hyper-viscoelastic model for soft tissues using a thermodynamically consistent approach. In addition, a nonlinear function, dependent on strain and strain rate, was adopted to capture the nonlinear variation of viscosity during a loading process. To demonstrate the efficacy and versatility of this approach, the model was used to recreate the experimental results performed on different types of soft tissues. In all the cases, the simulation results were well matched ($R^2⩾0.99$) with the experimental data.     

The influence of humidity on the viscoelastic behaviour of human hair

The possibility of curling hair is attributed to the fact that the mechanical properties of hair fibres depend on time, temperature and humidity. In this study the dependence of the viscoelastic behaviour of human hair fibres on humidity is characterised in a bending deformation. An experimental set-up was used to perform bending relaxation measurements at different humidities. The relaxation data were fitted using a so-called "stretched exponential equation". The humidity dependence was incorporated by applying time-humidity superposition. Also, the influence of humidity on the initial E-modulus was found. The influence of humidity on the behaviour of human hair fibres could be modelled by using a general characterisation of the behaviour of human hair fibres.The general characterisation is used to predict the recovery in a human hair fibre after curling. The predicted recovery appeared to agree reasonably well with experimentally determined values.     

Material properties of the human calcaneal fat pad in compression: experiment and theory

The structural behaviour of the human heel pad has been studied extensively due to its ability to absorb shock, protect against excessive local stress, and reduce plantar pressures. However, the material properties of the tissue have not been adequately measured. These must be known in order to perform a finite element analysis of the effect of factors such as foot geometry and shoe/surface construction on heel pad function. Therefore, the purposes of this study were to (a) measure the viscoelastic behaviour of the fat pad in compression, and (b) to determine an appropriate constitutive equation to model the tissue. A series of unconfined compression tests were performed on 8 mm diameter cylinders of fat pad tissue, consisting of quasi-static, 175, 350 mm/s and stress-relaxation tests to 50% deformation. The tissue exhibited nonlinear, viscoelastic behaviour. No significant difference was found in the quasi-static behaviour between samples from different locations and orientations in the heel. The stress-relaxation tests were used to determine the time constant (τ₁=0.5 s), the 175 mm/s test to determine the relaxation coefficient (g₁=28), and the 350 mm/s compression test to determine the material constants (C₁₀₀=C₀₁₀=0.01, C₂₀₀=C₀₂₀=0.1 Pa) of a single-phase, hyperelastic, linear viscoelastic strain energy function (r²=0.98).     

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.     

Analysis of several biomechanical models for the simulation of lamb liver behaviour using similarity coefficients from medical image

In this study, six biomechanical models for simulating lamb liver behaviour are presented. They are validated using similarity coefficients from Medical Image on reconstructed volumes from computerised tomography images. In particular, the Jaccard and Hausdorff coefficients are used. Loads of 20 and 40 g are applied to the livers and their deformation is simulated by means of the finite element method. The models used are a linear elastic model, a neo-Hookean model, a Mooney–Rivlin model, an Ogden model, a linear viscoelastic model and a viscohyperelastic model. The model that provided a behaviour that is closest to reality was the viscohyperelastic model, where the hyperelastic part was modelled with an Ogden model.     

Viscoelastic properties of aqueous solutions of an exocellular polysaccharide from cyanobacteria

The viscoelastic properties of aqueous solutions of the exocellular polysaccharide of Cyanospira capsulata have been studied, over a wide range of polymer concentrations, using small deformation oscillatory, steady and transient shear methods. The viscoelastic spectra generally resemble those of an entangled network, although notable deviations can be observed in the low frequency dependence of G′ and G″. At higher polymer concentrations, the viscoelastic spectrum shows solid-like behaviour over a wide range of frequencies. The superposition of η^*(ω) and η( ) curves occurs only at low frequencies, at higher frequencies the slope of η^*(ω) is lower than that of η( ). By studying the time evolution of shear stress after the inception of a steady shear rate (stress overshoot), the recovery of non-linear properties after steady shearing flow is seen to occur after times of c. 10³ s (in the case of 1·1% w/v solutions).

The overall viscoelastic properties appear original in comparison with those of the two structurally limiting types of polysaccharide, the ‘ordered’ chain xanthan and the ‘random coil’ guar. A rationale for this ‘anomalous’ viscoelastic behaviour can be tentatively proposed in terms of flickering intermolecular cross-interactions between semi-flexible segments, which occur in addition to the usual topological constraints.     

Measurement of the dynamic shear modulus of mouse brain tissue in vivo by magnetic resonance elastography

In this study, the magnetic resonance (MR) elastography technique was used to estimate the dynamic shear modulus of mouse brain tissue in vivo. The technique allows visualization and measurement of mechanical shear waves excited by lateral vibration of the skull. Quantitative measurements of displacement in three dimensions during vibration at 1200 Hz were obtained by applying oscillatory magnetic field gradients at the same frequency during a MR imaging sequence. Contrast in the resulting phase images of the mouse brain is proportional to displacement. To obtain estimates of shear modulus, measured displacement fields were fitted to the shear wave equation. Validation of the procedure was performed on gel characterized by independent rheometry tests and on data from finite element simulations. Brain tissue is, in reality, viscoelastic and nonlinear. The current estimates of dynamic shear modulus are strictly relevant only to small oscillations at a specific frequency, but these estimates may be obtained at high frequencies (and thus high deformation rates), noninvasively throughout the brain. These data complement measurements of nonlinear viscoelastic properties obtained by others at slower rates, either ex vivo or invasively.     

Inverse analysis and robustness evaluation for biological structure behaviour in FE simulation: application to the liver

To prevent traumas to abdominal organs, the selection of efficient safety devices should be based on a detailed knowledge of injury mechanisms and related injury criteria. In this sense, finite element (FE) simulation coupled with experiment could be a valuable tool to provide a better understanding of the behaviour of internal organs under crash conditions. This work proposes a methodology based on inverse analysis which combines exploration process optimisation and robustness study to obtain mechanical behaviour of the complex structure of the liver through FE simulation. The liver characterisation was based on Mooney–Rivlin hyperelastic behaviour law considering whole liver structure under uniform quasi-static compression. With the global method used, the model fits experimental data. The variability induced by modelling parameters is quantified within a reasonable time.     

Wide-range dynamic magnetic resonance elastography

Tissue mechanical parameters have been shown to be highly sensitive to disease by elastography. Magnetic resonance elastography (MRE) in the human body relies on the low-dynamic range of tissue mechanics <100 Hz. In contrast, MRE suited for investigations of mice or small tissue samples requires vibration frequencies 10–20 times higher than those used in human MRE. The dispersion of the complex shear modulus (G^?) prevents direct comparison of elastography data at different frequency bands and, consequently, frequency-independent viscoelastic models that fit to G^* over a wide dynamic range have to be employed. This study presents data of G^* of samples of agarose gel, liver, brain, and muscle measured by high-resolution MRE in a 7T-animal scanner at 200–800 Hz vibration frequency. Material constants μ and α according to the springpot model and related to shear elasticity and slope of the G^*-dispersion were determined. Both μ and α of calf brain and bovine liver were found to be similar, while a sample of fibrotic human liver (METAVIR score of 3) displayed about fifteen times higher shear elasticity, similar to μ of bovine muscle measured in muscle fiber direction. α was the highest in fibrotic liver, followed by normal brain and liver, while muscle had the lowest α-values of all biological samples investigated in this study. As expected, the least G^*-dispersion was seen in soft gel. The proposed technique of wide-range dynamic MRE can provide baseline data for both human MRE and high-dynamic MRE for better understanding tissue mechanics of different tissue structures.     

Magnetic resonance elastography in nonlinear viscoelastic materials under load

Characterisation of soft tissue mechanical properties is a topic of increasing interest in translational and clinical research. Magnetic resonance elastography (MRE) has been used in this context to assess the mechanical properties of tissues in vivo noninvasively. Typically, these analyses rely on linear viscoelastic wave equations to assess material properties from measured wave dynamics. However, deformations that occur in some tissues (e.g. liver during respiration, heart during the cardiac cycle, or external compression during a breast exam) can yield loading bias, complicating the interpretation of tissue stiffness from MRE measurements. In this paper, it is shown how combined knowledge of a material’s rheology and loading state can be used to eliminate loading bias and enable interpretation of intrinsic (unloaded) stiffness properties. Equations are derived utilising perturbation theory and Cauchy’s equations of motion to demonstrate the impact of loading state on periodic steady-state wave behaviour in nonlinear viscoelastic materials. These equations demonstrate how loading bias yields apparent material stiffening, softening and anisotropy. MRE sensitivity to deformation is demonstrated in an experimental phantom, showing a loading bias of up to twofold. From an unbiased stiffness of \(4910.4 \pm 635.8\) Pa in unloaded state, the biased stiffness increases to 9767.5 \(\pm \,\)1949.9 Pa under a load of \(\approx \) 34% uniaxial compression. Integrating knowledge of phantom loading and rheology into a novel MRE reconstruction, it is shown that it is possible to characterise intrinsic material characteristics, eliminating the loading bias from MRE data. The framework introduced and demonstrated in phantoms illustrates a pathway that can be translated and applied to MRE in complex deforming tissues. This would contribute to a better assessment of material properties in soft tissues employing elastography.

     

Experimental investigations of the human oesophagus: anisotropic properties of the embalmed mucosa–submucosa layer under large deformation

Mechanical characterisation of the layer-specific, viscoelastic properties of the human oesophagus is crucial in furthering the development of devices emerging in the field, such as robotic endoscopic biopsy devices, as well as in enhancing the realism, and therefore effectiveness, of surgical simulations. In this study, the viscoelastic and stress-softening behaviour of the passive human oesophagus was investigated through ex vivo cyclic mechanical tests. Due to restrictions placed on the laboratory as a result of COVID-19, only oesophagi from cadavers fixed in formalin were allowed for testing. Three oesophagi in total were separated into their two main layers and the mucosa–submucosa layer was investigated. A series of uniaxial tensile tests were conducted in the form of increasing stretch level cyclic tests at two different strain rates: 1% s\(^{-1}\) and 10% s\(^{-1}\). Rectangular samples in both the longitudinal and circumferential directions were tested to observe any anisotropy. Histological analysis was also performed through a variety of staining methods. Overall, the longitudinal direction was found to be much stiffer than the circumferential direction. Stress-softening was observed in both directions, as well as permanent set and hysteresis. Strain rate-dependent behaviour was also apparent in the two directions, with an increase in strain rate resulting in an increase in stiffness. This strain rate dependency was more pronounced in the longitudinal direction than the circumferential direction. Finally, the results were discussed in regard to the histological content of the layer, and the behaviour was modelled and validated using a visco-hyperelastic matrix-fibre model.