Characterization and assessment of hyperelastic and elastic properties of decellularized human adipose tissues

Decellularized adipose tissue (DAT) has shown potential as a regenerative scaffold for plastic and reconstructive surgery to augment or replace damaged or missing adipose tissue (e.g. following lumpectomy or mastectomy). The mechanical properties of soft tissue substitutes are of paramount importance in restoring the natural shape and appearance of the affected tissues, and mechanical mismatching can lead to unpredictable scar tissue formation and poor implant integration. The goal of this work was to assess the linear elastic and hyperelastic properties of decellularized human adipose tissue and compare them to those of normal breast adipose tissue. To assess the influence of the adipose depot source on the mechanical properties of the resultant decellularized scaffolds, we performed indentation tests on DAT samples sourced from adipose tissue isolated from the breast, subcutaneous abdominal region, omentum, pericardial depot and thymic remnant, and their corresponding force–displacement data were acquired. Elastic and hyperelastic parameters were estimated using inverse finite element algorithms. Subsequently, a simulation was conducted in which the estimated hyperelastic parameters were tested in a real human breast model under gravity loading in order to assess the suitability of the scaffolds for implantation. Results of these tests showed that in the human breast, the DAT would show similar deformability to that of native normal tissue. Using the measured hyperelastic parameters, we were able to assess whether DAT derived from different depots exhibited different intrinsic nonlinearities. Results showed that DAT sourced from varying regions of the body exhibited little intrinsic nonlinearity, with no statistically significant differences between the groups.     

Method for characterizing viscoelasticity of human gluteal tissue

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

A rheological network model for the continuum anisotropic and viscoelastic behavior of soft tissue

The mechanical behavior of soft tissue demonstrates a number of complex features including nonlinearity, anisotropy, viscoelasticity, and growth. Characteristic features of the time-dependent and anisotropic behavior are related to the properties of various components of the tissue such as fibrous collagen and elastin networks, large proteins and sugars attached to these networks, and interstitial fluid. Attempts to model the elastic behavior of these tissues based on assumptions about the behavior of the underlying constituents have been reasonably successful, but the essential addition of viscoelasticity to these models has been met with varying success. Here, a new rheological network model is proposed using, as its basis, an orthotropic hyperelastic constitutive model for fibrous tissue and a viscoelastic reptation model for soft materials. The resulting model has been incorporated into numerical and computational models, and is shown to capture the mechanical behavior of soft tissue in various modes of deformation including uniaxial and biaxial tension and simple shear.     

Muscle-driven finite element simulation of human foot movements

This paper describes a finite element scheme for realistic muscle-driven simulation of human foot movements. The scheme is used to simulate human ankle plantar flexion. A three-dimensional anatomically detailed finite element model of human foot and lower leg is developed and the idea of generating natural foot movement based entirely on the contraction of the plantar flexor muscles is used. The bones, ligaments, articular cartilage, muscles, tendons, as well as the rest soft tissues of human foot and lower leg are included in the model. A realistic three-dimensional continuum constitutive model that describes the biomechanical behaviour of muscles and tendons is used. Both the active and passive properties of muscle tissue are accounted for. The materials for bones and ligaments are considered as homogeneous, isotropic and linearly elastic, whereas the articular cartilage and the rest soft tissues (mainly fat) are defined as hyperelastic materials. The model is used to estimate muscle tissue deformations as well as stresses and strains that develop in the lower leg muscles during plantar flexion of the ankle. Stresses and strains that develop in Achilles tendon during such a movement are also investigated.     

Muscle-driven finite element simulation of human foot movements

This paper describes a finite element scheme for realistic muscle-driven simulation of human foot movements. The scheme is used to simulate human ankle plantar flexion. A three-dimensional anatomically detailed finite element model of human foot and lower leg is developed and the idea of generating natural foot movement based entirely on the contraction of the plantar flexor muscles is used. The bones, ligaments, articular cartilage, muscles, tendons, as well as the rest soft tissues of human foot and lower leg are included in the model. A realistic three-dimensional continuum constitutive model that describes the biomechanical behaviour of muscles and tendons is used. Both the active and passive properties of muscle tissue are accounted for. The materials for bones and ligaments are considered as homogeneous, isotropic and linearly elastic, whereas the articular cartilage and the rest soft tissues (mainly fat) are defined as hyperelastic materials. The model is used to estimate muscle tissue deformations as well as stresses and strains that develop in the lower leg muscles during plantar flexion of the ankle. Stresses and strains that develop in Achilles tendon during such a movement are also investigated.     

Calibration of hyperelastic constitutive models: the role of boundary conditions,search algorithms,and experimental variability

Biomechanics and Modeling in Mechanobiology - The calibration of hyperelastic constitutive models of soft tissue and tissue surrogates is often treated as an exercise in curve-fitting to the...     

Spherical indentation method for determining the constitutive parameters of hyperelastic soft materials

A comprehensive study on the spherical indentation of hyperelastic soft materials is carried out through combined theoretical, computational, and experimental efforts. Four widely used hyperelastic constitutive models are studied, including neo-Hookean, Mooney–Rivlin, Fung, and Arruda–Boyce models. Through dimensional analysis and finite element simulations, we establish the explicit relations between the indentation loads at given indentation depths and the constitutive parameters of materials. Based on the obtained results, the applicability of Hertzian solution to the measurement of the initial shear modulus of hyperelastic materials is examined. Furthermore, from the viewpoint of inverse problems, the possibility to measure some other properties of a hyperelastic material using spherical indentation tests, e.g., locking stretch, is addressed by considering the existence, uniqueness, and stability of the solution. Experiments have been performed on polydimethylsiloxane to validate the conclusions drawn from our theoretical analysis. The results reported in this study should help identify the extent to which the mechanical properties of hyperelastic materials could be measured from spherical indentation tests.     

Three-dimensional finite element analysis of the foot during standing--a material sensitivity study

Information on the internal stresses/strains in the human foot and the pressure distribution at the plantar support interface under loading is useful in enhancing knowledge on the biomechanics of the ankle-foot complex. While techniques for plantar pressure measurements are well established, direct measurement of the internal stresses/strains is difficult. A three-dimensional (3D) finite element model of the human foot and ankle was developed using the actual geometry of the foot skeleton and soft tissues, which were obtained from 3D reconstruction of MR images. Except the phalanges that were fused, the interaction among the metatarsals, cuneiforms, cuboid, navicular, talus, calcaneus, tibia and fibula were defined as contact surfaces, which allow relative articulating movement. The plantar fascia and 72 major ligaments were simulated using tension-only truss elements by connecting the corresponding attachment points on the bone surfaces. The bony and ligamentous structures were embedded in a volume of soft tissues. The encapsulated soft tissue was defined as hyperelastic, while the bony and ligamentous structures were assumed to be linearly elastic. The effects of soft tissue stiffening on the stress distribution of the plantar surface and bony structures during balanced standing were investigated. Increases of soft tissue stiffness from 2 and up to 5 times the normal values were used to approximate the pathologically stiffened tissue behaviour with increasing stages of diabetic neuropathy. The results showed that a five-fold increase in soft tissue stiffness led to about 35% and 33% increase in the peak plantar pressure at the forefoot and rearfoot regions, respectively. This corresponded to about 47% decrease in the total contact area between the plantar foot and the horizontal support surface. Peak bone stress was found at the third metatarsal in all calculated cases with a minimal increase of about 7% with soft tissue stiffening.     

A model of growth and rupture of abdominal aortic aneurysm

We present here a coupled mathematical model of growth and failure of the abdominal aortic aneurysm (AAA). The failure portion of the model is based on the constitutive theory of softening hyperelasticity where the classical hyperelastic law is enhanced with a new constant indicating the maximum energy that an infinitesimal material volume can accumulate without failure. The new constant controls material failure and it can be interpreted as the average energy of molecular bonds from the microstructural standpoint. The constitutive model is compared to the data from uniaxial tension tests providing an excellent fit to the experiment. The AAA failure model is coupled with a phenomenological theory of soft tissue growth. The unified theory includes both momentum and mass balance laws coupled with the help of the constitutive equations. The microstructural alterations in the production of elastin and remodeling of collagen are reflected in the changing macroscopic parameters characterizing tissue stiffness, strength and density. The coupled theory is used to simulate growth and rupture of an idealized spherical AAA. The results of the simulation showing possible AAA ruptures in growth are reasonable qualitatively while the quantitative calibration of the model will require further clinical observations and in vitro tests. The presented model is the first where growth and rupture are coupled.     

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.     

Constitutive formulation and numerical analysis of the heel pad region

The aim of this work is to provide a numerical approach for the investigation of the mechanical behaviour of the heel pad region. A visco-hyperelastic model is formulated with regard to fat pad tissue, while a fibre-reinforced hyperelastic model is considered for the heel skin tissue. Bone components are defined by means of an orthotropic linear elastic model. Particular attention is paid to the evaluation of constitutive parameters within different models adopted in consideration of experimental tests data. Preliminarily, indentation tests on a skinless cadaveric foot are considered with regard to fat pad tissue. Indentation tests on an intact heel pad of a cadaveric foot are subsequently adopted for the final identification of constitutive parameters of fat pad and skin tissues. A numerical model of the rear foot is defined and different loading conditions are assumed according to experimental data. A comparison between experimental and numerical data leads to the evaluation of the real capability of the procedure to interpret the actual response of the rear foot.     

Constitutive formulation and numerical analysis of the heel pad region

The aim of this work is to provide a numerical approach for the investigation of the mechanical behaviour of the heel pad region. A visco-hyperelastic model is formulated with regard to fat pad tissue, while a fibre-reinforced hyperelastic model is considered for the heel skin tissue. Bone components are defined by means of an orthotropic linear elastic model. Particular attention is paid to the evaluation of constitutive parameters within different models adopted in consideration of experimental tests data. Preliminarily, indentation tests on a skinless cadaveric foot are considered with regard to fat pad tissue. Indentation tests on an intact heel pad of a cadaveric foot are subsequently adopted for the final identification of constitutive parameters of fat pad and skin tissues. A numerical model of the rear foot is defined and different loading conditions are assumed according to experimental data. A comparison between experimental and numerical data leads to the evaluation of the real capability of the procedure to interpret the actual response of the rear foot.     

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.     

An approach to the mechanical constitutive modelling of arterial tissue based on homogenization and optimization

This paper is concerned with characterizing the quasistatic mechanical behaviour of arterial tissue undergoing finite deformation through hyperelastic constitutive functions. Commonly the parameters of constitutive functions are established by a process of optimization based on experimental data. Instead we construct a finite element model of a representative volume element of the material and compute its homogenized response to a range of deformations. These data are then used to provide objective functions for optimizing the parameters of two analytical models from the literature.     

In vitro formation of fibrous septa by liver connective tissue cells

Summary Active fibrous septa are a common feature in liver fibrosis and cirrhosis. Their etiology and formation were studied using cultures of tissue fragments or cells included in collagen gels. Liver fragments obtained from patients with cirrhosis or severe schistosomal fibrosis were able to reorganize the gel and to form discrete, interconnecting fibrous septa composed of parallel arrays of collagen, subsequently colonized by migrating connective tissue cells. The same was obtained in cultures of fibrogranulomatous lesions isolated from schistosome-infected mice livers. However, fragments of normal human and murine liver tissue did not show the capacity to form fibrous septa. Septa formation was also obtained in cultures of cell spheroids formed by liver connective tissue cells isolated from human fibrotic or cirrhotic liver tissues, but not with spheroids of normal skin fibroblasts or smooth muscle cells. This experimental model may represent the fibrous septa formation in vivo, depending on the activity of liver connective tissue cells. The ability of tissue fragments or cell spheroids to form septa in collagen gels might reflect the degree of fibrosis present in the liver tissue in vivo. This research was supported by FINEP and CNPq (Brazil) and CNRS (France).     

Adipose tissue engineering with cells in engineered matrices

Tissue engineering has shown promise for the development of constructs to facilitate large volume soft tissue augmentation in reconstructive and cosmetic plastic surgery. This article reviews the key progress to date in the field of adipose tissue engineering. In order to effectively design a soft tissue substitute, it is critical to understand the native tissue environment and function. As such, the basic physiology of adipose tissue is described and the process of adipogenesis is discussed. In this article, we have focused on tissue engineering using a cell-seeded scaffold approach, where engineered extracellular matrix substitutes are seeded with exogenous cells that may contribute to the regenerative response. The strengths and limitations of each of the possible cell sources for adipose tissue engineering, including adipose-derived stem cells, are detailed. We briefly highlight some of the results from the major studies to date, involving a range of synthetic and naturally derived scaffolds. While these studies have shown that adipose tissue regeneration is possible, more research is required to develop optimized constructs that will facilitate safe, predictable, and long-term augmentation in clinical applications.     

Digital image correlation and finite element modelling as a method to determine mechanical properties of human soft tissue in vivo

The mechanical properties of human soft tissue are crucial for impact biomechanics, rehabilitation engineering, and surgical simulation. Validation of these constitutive models using human data remains challenging and often requires the use of non-invasive imaging and inverse finite element (FE) analysis. Post-processing data from imaging methods such as tagged magnetic resonance imaging (MRI) can be challenging. Digital image correlation (DIC), however, is a relatively straightforward imaging method. DIC has been used in the past to study the planar and superficial properties of soft tissue and excised soft tissue layers. However, DIC has not been used to non-invasive study of the bulk properties of human soft tissue in vivo. Thus, the goal of this study was to assess the use of DIC in combination with FE modelling to determine the bulk material properties of human soft tissue. Indentation experiments were performed on a silicone gel soft tissue phantom. A two camera DIC setup was then used to record the 3D surface deformation. The experiment was then simulated using a FE model. The gel was modelled as Neo-Hookean hyperelastic, and the material parameters were determined by minimising the error between the experimental and FE data. The iterative FE analysis determined material parameters (μ=1.80 kPa, K=2999 kPa) that were in close agreement with parameters derived independently from regression to uniaxial compression tests (μ=1.71 kPa, K=2857 kPa). Furthermore the FE model was capable of reproducing the experimental indentor force as well as the surface deformation found (R²=0.81). It was therefore concluded that a two camera DIC configuration combined with FE modelling can be used to determine the bulk mechanical properties of materials that can be represented using hyperelastic Neo-Hookean constitutive laws.     

Mathematical Model for Skeletal Muscle to Simulate the Concentric and Eccentric Contraction

Skeletal muscles are responsible for the relative motion of the bones at the joints and provide the required strength. They exhibit highly nonlinear mechanical behaviour and are described by nonlinear hyperelastic constitutive relations. It is distinct from other biological soft tissue. Its hyperelastic or viscoelastic behaviour is modelled by using CE, SEE, and PEE. Contractile element simulates the behaviour of skeletal muscle when it is subjected to eccentric and concentric contraction. This research aims to estimate the stress induced in skeletal muscle in eccentric and concentric contraction with respect to the predefined strain. With the use of mathematical model for contraction of skeletal muscle for eccentric and concentric contraction, the stress induced in the skeletal muscle is estimated in this research. Mathematical model is developed for the muscle using EMG signals and Force-velocity relationship calculated. With the use of force-velocity of contraction of muscle, mathematical model is developed. This can be useful to understand the mechanical behaviour of skeletal muscles in eccentric and concentric contraction with clinical relevance. Authors are further working to develop the mathematical model with torsion force with proper activation function of muscle and experimentation for extraction of the anisotropic mechanical properties of skeletal muscle.