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.     

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.     

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.     

Mechanical characterisation of in vivo human skin using a 3D force-sensitive micro-robot and finite element analysis

The complex mechanical properties of skin have been the subject of much study in recent years. Several experimental methods developed to measure the mechanical properties of skin in vivo, such as suction or torsion, are unable to measure skin’s anisotropic characteristics. An experiment characterising the mechanical properties of in vivo human skin using a novel force-sensitive micro-robot is presented. The micro-robot applied in-plane deformations to the anterior forearm and the posterior upper arm. The behaviour of the skin in each area is highly nonlinear, anisotropic, and viscoelastic. The response of the upper arm skin is very dependent on the orientation of the arm. A finite element model consisting of an Ogden strain energy function and quasi-linear viscoelasticity was developed to simulate the experiments. An orthogonal initial stress field, representing the in vivo skin tension, was used as an additional model parameter. The model simulated the experiments accurately with an error-of-fit of 17.5% for the anterior lower forearm area, 6.5% for the anterior upper forearm and 9.3% for the posterior upper arm. The maximum in vivo tension in each area determined by the model was 6.2 Nm⁻¹ in the anterior lower forearm, 11.4 Nm⁻¹ in anterior upper forearm and 5.6 Nm⁻¹ in the posterior upper arm. The results also show that a finite element model with a neo-Hookean strain energy function cannot simulate the experiments with the same accuracy.     

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.     

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.     

A Bimodular Theory for Finite Deformations: Comparison of Orthotropic Second-order and Exponential Stress Constitutive Equations for Articular Cartilage

Cartilaginous tissues, such as articular cartilage and the annulus fibrosus, exhibit orthotropic behavior with highly asymmetric tensile–compressive responses. Due to this complex behavior, it is difficult to develop accurate stress constitutive equations that are valid for finite deformations. Therefore, we have developed a bimodular theory for finite deformations of elastic materials that allows the mechanical properties of the tissue to differ in tension and compression. In this paper, we derive an orthotropic stress constitutive equation that is second-order in terms of the Biot strain tensor as an alternative to traditional exponential type equations. Several reduced forms of the bimodular second-order equation, with six to nine parameters, and a bimodular exponential equation, with seven parameters, were fit to an experimental dataset that captures the highly asymmetric and orthotropic mechanical response of cartilage. The results suggest that the bimodular second-order models may be appealing for some applications with cartilaginous tissues.     

Carotid artery mechanical properties and stresses quantified using in vivo data from normotensive and hypertensive humans

The goal of this study was to model the in vivo non-linear mechanical behavior of human common carotid arteries (CCAs) and then to compare wall stresses and associated contributions of micro-constituents in normotensive (NT) and treated hypertensive (HT) subjects. We used an established theoretical model of 3D arterial mechanics that assumes a hyperelastic, anisotropic, active–passive, and residually stressed wall. In vivo data were obtained non-invasively from CCAs in 16 NT (21–64 years old) and 25 treated HT (44–69 years old) subjects. The associated quasi-static boundary value problem was solved semi-analytically over a cardiac cycle while accounting for surrounding perivascular tissue. Best-fit values of model parameters, including those describing contributions by intramural elastin, fibrillar collagen, and vascular smooth muscle, were estimated by a non-linear least-squares method. The model (1) captured temporal changes in intraluminal pressure, (2) estimated wall stress fields that appeared to reflect the presence or absence of age and disease, and (3) suggested changes in mechanical characteristics of wall micro-constituents despite medical treatment of hypertension. For example, age was positively correlated with residual stresses and altered fibrillar collagen in NT subjects, which indirectly validated the modeling, and HT subjects had higher levels of stresses, increased smooth muscle tone, and a stiffer elastin-dominated matrix despite treatment. These results are consistent with prior reports on effects of age and hypertension, but provide increased insight into evolving contributions of cell and matrix mechanics to arterial behavior in vivo.     

Estimating material parameters of a structurally based constitutive relation for skin mechanics

This paper presents a structurally based modeling framework to characterize the structure–function relation in skin tissues, based upon biaxial tensile experiments performed in vitro on porcine skin. Equi-axial deformations were imposed by stretching circular skin specimens uniformly along twelve directions, and the resultant loads at the membrane attachment points were measured. Displacement fields at each deformation step were tracked using an image 2D cross-correlation technique. A modeling framework was developed to simulate the experiments, whereby measured forces were applied to finite element models that were created to represent the geometry and structure of the tissue samples. Parameters of a structurally based constitutive relation were then identified using nonlinear optimization. Results showed that the ground matrix stiffness ranged from 5 to 32 kPa, fiber orientation mean from 2 to 13^◦ from the torso midline, fiber undulation mean from 1.04 to 1.34 and collagen fiber stiffness from 48 to 366 MPa. It was concluded that the objective function was highly sensitive to the mean orientation and that a priori information about fiber orientation mean was important for the reliable identification of constitutive parameters.     

Modeling interlamellar interactions in angle-ply biologic laminates for annulus fibrosus tissue engineering

Mechanical function of the annulus fibrosus of the intervertebral disc is dictated by the composition and microstructure of its highly ordered extracellular matrix. Recent work on engineered angle-ply laminates formed from mesenchymal stem cell (MSC)-seeded nanofibrous scaffolds indicates that the organization of collagen fibers into planes of alternating alignment may play an important role in annulus fibrosus tissue function. Specifically, these engineered tissues can resist tensile deformation through shearing of the interlamellar matrix as layers of collagen differentially reorient under load. In the present work, a hyperelastic constitutive model was developed to describe the role of interlamellar shearing in reinforcing the tensile response of biologic laminates, and was applied to experimental results from engineered annulus constructs formed from MSC-seeded nanofibrous scaffolds. By applying the constitutive model to uniaxial tensile stress–strain data for bilayers with three different fiber orientations, material parameters were generated that characterize the contributions of extrafibrillar matrix, fibers, and interlamellar shearing interactions. By 10 weeks of in vitro culture, interlamellar shearing accounted for nearly 50% of the total stress associated with uniaxial extension in the anatomic range of ply angle. The model successfully captured changes in function with extracellular matrix deposition through variations in the magnitude of model parameters with culture duration. This work illustrates the value of engineered tissues as tools to further our understanding of structure–function relations in native tissues and as a test-bed for the development of constitutive models to describe them.     

Determination of material parameters of the two-dimensional Holzapfel–Weizsäcker type model based on uniaxial extension data of arterial walls

Soft tissues are anisotropic materials yet a majority of mechanical property tests have been uniaxial, which often failed to recapitulate the tensile response in other directions. This paper aims to study the feasibility of determining material parameters of anisotropic tissues by uniaxial extension with a minimal loss of anisotropic information. We assumed that by preselecting a certain constitutive model, we could give the constitutive parameters based on uniaxial extension data from orthogonal strip samples. In our study, the Holzapfel–Weizsäcker type strain energy density function (H–W model) was used to determine the material parameters of arterial walls from two fresh donation bodies. The key points we applied were the relationships between strain components in uniaxial tensile tests and the methods of stochastic optimisation. Further numerical experiments were taken. The estimate–effect ratio, defined by the number of data with the precision of estimation less than 0.5% over whole size of data, was calculated to demonstrate the feasibility of our method. The material parameters for Chinese aorta and pulmonary artery were given with the maximum root mean square (RMS) errors 0.042, and the minimal estimate–effect ratio in numerical experiments was 90.79%. Our results suggest that the constitutive parameters of arterial walls can be determined from uniaxial extension data, given the passive mechanical behaviour governed by H–W model. This method may apply to other tissues using different constitutive models.     

Elastic and viscoelastic properties of porcine subdermal fat using MRI and inverse FEA

There is a scarcity of investigation into the mechanical properties of subdermal fat. Recently, progress has been made in the determination of subdermal stress and strain distributions. This requires accurate constitutive modelling and consideration of the subdermal tissues. This paper reports the results of a study to estimate non-linear elastic and viscoelastic properties of porcine subdermal fat using a simple constitutive model. High-resolution magnetic resonance imaging (MRI) was used to acquire a time series of coincident images during a confined indentation experiment. Inverse finite element analysis was used to estimate the material parameters. The Neo Hookean model was used to represent the elastic behaviour (μ = 0.53 ± 0.31 kPa), while a single-element Prony series was used to model the viscoelastic response (α = 0.39 ± 0.03, τ = 700 ± 255 s).     

In vivo characterization of the mechanical properties of human skin derived from MRI and indentation techniques

The human skin is an exceedingly complex and multi-layered material. This paper aims to introduce the application of the finite element analysis (FEA) to the in vivo characterization of the non-linear mechanical behaviour of three human skin layers. Indentation tests combined with magnetic resonance imaging (MRI) technique have been performed on the left dorsal forearm of a young man in order to reveal the mechanical behaviour of all skin layers. Using MRI images processing and a pre and post processor allows to make numerically individualized 2D model which consists of three skin layers and the muscles. FEA has been applied to simulate indentation tests. Neo-Hookean slightly compressible material model of two material constants (C(10), K) has been used to model the mechanical behaviour of the three skin layers and the muscles. The identification of material model parameters was done by applying Levenberg-Marquardt algorithm (LMA). Our methodology of identification provides a range of values for each constant. Range of values of different material properties of epidermis, dermis, hypodermis are respectively, C10(E)=0.12+/-0.06 MPa, C10(D)=1.11+/-0.09 MPa, C10(H)=0.42+/-0.05 KPa, K(E)=5.45+/-1.7 MPa, K(D)=29.6+/-1,28 MPa, K(H)=36.0+/-0.9 KPa.     

A multiscale model for red blood cell mechanics

The objective of this article is the derivation of a continuum model for mechanics of red blood cells via multiscale analysis. On the microscopic level, we consider realistic discrete models in terms of energy functionals defined on networks/lattices. Using concepts of Γ-convergence, convergence results as well as explicit homogenisation formulae are derived. Based on a characterisation via energy functionals, appropriate macroscopic stress–strain relationships (constitutive equations) can be determined. Further, mechanical moduli of the derived macroscopic continuum model are directly related to microscopic moduli. As a test case we consider optical tweezers experiments, one of the most common experiments to study mechanical properties of cells. Our simulations of the derived continuum model are based on finite element methods and account explicitly for membrane mechanics and its coupling with bulk mechanics. Since the discretisation of the continuum model can be chosen freely, rather than it is given by the topology of the microscopic cytoskeletal network, the approach allows a significant reduction of computational efforts. Our approach is highly flexible and can be generalised to many other cell models, also including biochemical control.     

A novel two-layer, coupled finite element approach for modeling the nonlinear elastic and viscoelastic behavior of human erythrocytes

A novel finite element approach is presented to simulate the mechanical behavior of human red blood cells (RBC, erythrocytes). As the RBC membrane comprises a phospholipid bilayer with an intervening protein network, we propose to model the membrane with two distinct layers. The fairly complex characteristics of the very thin lipid bilayer are represented by special incompressible solid shell elements and an anisotropic viscoelastic constitutive model. Properties of the protein network are modeled with an isotropic hyperelastic third-order material. The elastic behavior of the model is validated with existing optical tweezers studies with quasi-static deformations. Employing material parameters consistent with literature, simulation results are in excellent agreement with experimental data. Available models in literature neglect either the surface area conservation of the RBC membrane or realistic loading conditions of the optical tweezers experiments. The importance of these modeling assumptions, that are both included in this study, are discussed and their influence quantified. For the simulation of the dynamic motion of RBC, the model is extended to incorporate the cytoplasm. This is realized with a monolithic fully coupled fluid-structure interaction simulation, where the fluid is described by the incompressible Navier–Stokes equations in an arbitrary Lagrangian Eulerian framework. It is shown that both membrane viscosity and cytoplasm viscosity have significant influence on simulation results. Characteristic recovery times and energy dissipation for varying strain rates in dynamic laser trap experiments are calculated for the first time and are found to be comparable with experimental data.     

In vivo measurement of surface skin strain during human gait to improve the design of rehabilitation devices

Abstract

When designing any rehabilitation, sportswear or exoskeleton device the mechanical behaviour of the body segment must be known, specifically the skin, because an excessive tissue strain may lead to ulceration and bedsores. To date, it is not known if the kinematic variability between subjects have an effect on the skin strain field, and therefore, in the design and manufacturing of rehabilitation products, such as orthoses. Several studies have analysed the skin deformation during human motion, nevertheless, the comparison between the skin strain field in different subjects during normal or pathological gait has not been reported yet. This work presents a comparison of skin strain analysis for different gait patterns to study the differences between people and, specifically, if it is possible to standardize the orthotic design between subjects with the same gait disorder. Moreover, the areas with relatively minimum strain during the ankle-foot motion are compared to improve the design of structural parts of rehabilitation devices. In this case, a validated 3D digital image correlation system has been used for this purpose combined with strain ellipse theory. The results demonstrate variations in the skin strain field between subjects with the same pathology and similarities between subjects with normal gait. However, more studies and experiments are necessaries to validate this hypothesis and also to test it between different gait pathologies.