An orthotropic viscoelastic material model for passive myocardium: theory and algorithmic treatment

This contribution presents a novel constitutive model in order to simulate an orthotropic rate-dependent behaviour of the passive myocardium at finite strains. The motivation for the consideration of orthotropic viscous effects in a constitutive level lies in the disagreement between theoretical predictions and experimentally observed results. In view of experimental observations, the material is deemed as nearly incompressible, hyperelastic, orthotropic and viscous. The viscoelastic response is formulated by means of a rheological model consisting of a spring coupled with a Maxwell element in parallel. In this context, the isochoric free energy function is decomposed into elastic equilibrium and viscous non-equilibrium parts. The baseline elastic response is modelled by the orthotropic model of Holzapfel and Ogden [Holzapfel GA, Ogden RW. 2009. Constitutive modelling of passive myocardium: a structurally based framework for material characterization. Philos Trans Roy Soc A Math Phys Eng Sci. 367:3445–3475]. The essential aspect of the proposed model is the account of distinct relaxation mechanisms for each orientation direction. To this end, the non-equilibrium response of the free energy function is constructed in the logarithmic strain space and additively decomposed into three anisotropic parts, denoting fibre, sheet and normal directions each accompanied by a distinct dissipation potential governing the evolution of viscous strains associated with each orientation direction. The evolution equations governing the viscous flow have an energy-activated nonlinear form. The energy storage in the Maxwell branches has a quadratic form leading to a linear stress–strain response in the logarithmic strain space. On the numerical side, the algorithmic aspects suitable for the implicit finite element method are discussed in a Lagrangian setting. The model shows excellent agreement compared to experimental data obtained from the literature. Furthermore, the finite element simulations of a heart cycle carried out with the proposed model show significant deviations in the strain field relative to the elastic solution.     

Leukocyte deformability: finite element modeling of large viscoelastic deformation.

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

Effect of colchicine on viscoelastic properties of neutrophils

The effect of colchicine (15-60 micrograms/ml) on the viscoelastic properties of human neutrophils was studied by the micropipette technique. The small deformation of the neutrophil in response to a step aspiration pressure was analyzed by using a three-element model in which an elastic element, K1, is in parallel with a Maxwell element composed of another elastic element, K2, in series with a viscous element, mu. Colchicine treatment of neutrophils caused decreases in K2 and mu without affecting K1. The results indicate that the integrity of the microtubules plays a significant role in providing the viscoelastic resistance (as represented by the Maxwell element in the model) of neutrophils to deforming stress.     

An experimental and theoretical approach of elasticity and viscoelasticity of compact and spongy bone with periodic homogenization

An experimental study on samples from a fresh human femur highlighted the viscoelastic behaviour of bone. The memory aspect of the behaviour is weak compared to the elastic aspect. A theoretical and numerical approach based on the finite element method illustrates the existing bond between the structure and the micro behaviour on the one hand, and the macro behaviour on the other. The orthotropic elastic behaviour observed by various authors is shown numerically. The memory functions of compact and spongy bone are given. The numerical results illustrate the behaviour observed in the experiments.     

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.     

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.     

Influence of physicochemical factors on rheology of human neutrophils.

The effects of variations in temperature, pH, and osmolality on the rheological properties of human neutrophils were determined by studying the cell deformation in response to aspirational pressure applied via a micropipette. The time history of the deformation was analyzed by the use of a standard solid viscoelastic model consisting of an elastic element K1 in parallel with a Maxwell element (an elastic element K2 in series with a viscous element mu). With changes in temperature over a range of 9-40 degrees C, only mu varied inversely with temperature, while K1 and K2 did not show significant alterations. Variations in pH over the range of 5.4-7.8 did not significantly affect the viscoelastic coefficients, but K1 and mu rose at pH 8.4. An increase in osmolality caused all three coefficients to rise, but a decrease in osmolality had relatively little effect on the coefficients. These changes in response to physicochemical variations serve to provide insights into the viscoelastic properties of neutrophils and their possible roles in health and disease.     

A novel porous mechanical framework for modelling the interaction between coronary perfusion and myocardial mechanics

The strong coupling between the flow in coronary vessels and the mechanical deformation of the myocardial tissue is a central feature of cardiac physiology and must therefore be accounted for by models of coronary perfusion. Currently available geometrically explicit vascular models fail to capture this interaction satisfactorily, are numerically intractable for whole organ simulations, and are difficult to parameterise in human contexts. To address these issues, in this study, a finite element formulation of an incompressible, poroelastic model of myocardial perfusion is presented. Using high-resolution ex vivo imaging data of the coronary tree, the permeability tensors of the porous medium were mapped onto a mesh of the corresponding left ventricular geometry. The resultant tensor field characterises not only the distinct perfusion regions that are observed in experimental data, but also the wide range of vascular length scales present in the coronary tree, through a multi-compartment porous model. Finite deformation mechanics are solved using a macroscopic constitutive law that defines the coupling between the fluid and solid phases of the porous medium. Results are presented for the perfusion of the left ventricle under passive inflation that show wall-stiffening associated with perfusion, and that show the significance of a non-hierarchical multi-compartment model within a particular perfusion territory.     

Equivalence between short-time biphasic and incompressible elastic material responses

Porous-permeable tissues have often been modeled using porous media theories such as the biphasic theory. This study examines the equivalence of the short-time biphasic and incompressible elastic responses for arbitrary deformations and constitutive relations from first principles. This equivalence is illustrated in problems of unconfined compression of a disk, and of articular contact under finite deformation, using two different constitutive relations for the solid matrix of cartilage, one of which accounts for the large disparity observed between the tensile and compressive moduli in this tissue. Demonstrating this equivalence under general conditions provides a rationale for using available finite element codes for incompressible elastic materials as a practical substitute for biphasic analyses, so long as only the short-time biphasic response is sought. In practice, an incompressible elastic analysis is representative of a biphasic analysis over the short-term response deltat相似文献   

Flexural and Creep Properties of Human Jaw Compact Bone for FEA Studies

The aim of this work was to improve the constitutive model of the human mandible and dentition system by taking into account the non-linear material properties of the structural boney matrix that forms the human jaw bone or mandible. Due to the specific structure of the jaw bone the time dependence of the mechanical properties also forms an important stage of the quantification process. The lack of specific experimental data of this type of material prevents the implementation of these properties into finite element simulations which results in poor quality modelling. Here an attempt was made to determine elastic and viscoelastic mechanical characteristics of the compact bone tissue forming the mandible. The elastic properties of compact bone were determined experimentally from 3 point bending tests and the viscoelastic properties were evaluated from creep tests in compression. A particular human jaw from this complex study was used to reconstruct a geometric model for further numerical experiments.     

Flexural and creep properties of human jaw compact bone for FEA studies

The aim of this work was to improve the constitutive model of the human mandible and dentition system by taking into account the non-linear material properties of the structural boney matrix that forms the human jaw bone or mandible. Due to the specific structure of the jaw bone the time dependence of the mechanical properties also forms an important stage of the quantification process. The lack of specific experimental data of this type of material prevents the implementation of these properties into finite element simulations which results in poor quality modelling. Here an attempt was made to determine elastic and viscoelastic mechanical characteristics of the compact bone tissue forming the mandible. The elastic properties of compact bone were determined experimentally from 3 point bending tests and the viscoelastic properties were evaluated from creep tests in compression. A particular human jaw from this complex study was used to reconstruct a geometric model for further numerical experiments.     

Constitutive properties of hypertrophied myocardium: cellular contribution to changes in myocardial stiffness

Recent studies have suggested that pressure overload hypertrophy (POH) alters the viscoelastic properties of individual cardiocytes when studied in isolation. However, whether these changes in cardiocyte properties contribute causally to changes in the material properties of the cardiac muscle as a whole is unknown. Accordingly, a selective, isolated, acute change in cardiocyte constitutive properties was imposed in an in vitro system capable of measuring the resultant effect on the material properties of the composite cardiac muscle. POH caused an increase in both myocardial elastic stiffness, from 20.5 +/- 1.3 to 28.4 +/- 1.8, and viscous damping, from 15.2 +/- 1.1 to 19.8 +/- 1.5 s (normal vs. POH, P < 0.05), respectively. Recent studies have shown that cardiocyte constitutive properties could be acutely altered by depolymerizing the microtubules with colchicine. Colchicine caused a significant decrease in the viscous damping in POH muscles (19.8 +/- 1.5 s at baseline vs. 14.7 +/- 1.3 s after colchicine, P < 0.05). Therefore, myocardial material properties can be altered by selectively changing the constitutive properties of one element within this muscle tissue, the cardiocyte. Changes in the constitutive properties of the cardiocytes themselves contribute to the abnormalities in myocardial stiffness and viscosity that develop during POH.     

In vivo liver tissue mechanical properties by Transient Elastography: comparison with Dynamic Mechanical Analysis

Understanding the mechanical properties of human liver is one of the most critical aspects of its numerical modeling for medical applications or impact biomechanics. Generally, model constitutive laws come from in vitro data. However, the elastic properties of liver may change significantly after death and with time. Furthermore, in vitro liver elastic properties reported in the literature have often not been compared quantitatively with in vivo liver mechanical properties on the same organ. In this study, both steps are investigated on porcine liver. The elastic property of the porcine liver, given by the shear modulus G, was measured by both Transient Elastography (TE) and Dynamic Mechanical Analysis (DMA). Shear modulus measurements were realized on in vivo and in vitro liver to compare the TE and DMA methods and to study the influence of testing conditions on the liver viscoelastic properties. In vitro results show that elastic properties obtained by TE and DMA are in agreement. Liver tissue in the frequency range from 0.1 to 4 Hz can be modeled by a two-mode relaxation model. Furthermore, results show that the liver is homogeneous, isotropic and more elastic than viscous. Finally, it is shown in this study that viscoelastic properties obtained by TE and DMA change significantly with post mortem time and with the boundary conditions.     

Prediction of brain deformations and risk of traumatic brain injury due to closed-head impact: quantitative analysis of the effects of boundary conditions and brain tissue constitutive model

In this study, we investigate the effects of modelling choices for the brain–skull interface (layers of tissues between the brain and skull that determine boundary conditions for the brain) and the constitutive model of brain parenchyma on the brain responses under violent impact as predicted using computational biomechanics model. We used the head/brain model from Total HUman Model for Safety (THUMS)—extensively validated finite element model of the human body that has been applied in numerous injury biomechanics studies. The computations were conducted using a well-established nonlinear explicit dynamics finite element code LS-DYNA. We employed four approaches for modelling the brain–skull interface and four constitutive models for the brain tissue in the numerical simulations of the experiments on post-mortem human subjects exposed to violent impacts reported in the literature. The brain–skull interface models included direct representation of the brain meninges and cerebrospinal fluid, outer brain surface rigidly attached to the skull, frictionless sliding contact between the brain and skull, and a layer of spring-type cohesive elements between the brain and skull. We considered Ogden hyperviscoelastic, Mooney–Rivlin hyperviscoelastic, neo–Hookean hyperviscoelastic and linear viscoelastic constitutive models of the brain tissue. Our study indicates that the predicted deformations within the brain and related brain injury criteria are strongly affected by both the approach of modelling the brain–skull interface and the constitutive model of the brain parenchyma tissues. The results suggest that accurate prediction of deformations within the brain and risk of brain injury due to violent impact using computational biomechanics models may require representation of the meninges and subarachnoidal space with cerebrospinal fluid in the model and application of hyperviscoelastic (preferably Ogden-type) constitutive model for the brain tissue.     

A frame-invariant formulation of Fung elasticity

Fung elasticity refers to the hyperelasticity constitutive relation proposed by Fung and co-workers for describing the pseudo-elastic behavior of biological soft tissues undergoing finite deformation. A frame-invariant formulation of Fung elasticity is provided for material symmetries ranging from orthotropy to isotropy, which uses Lamé-like material constants. In the orthotropic case, three orthonormal vectors are used to define mutually orthogonal planes of symmetry and associated texture tensors. The strain energy density is then formulated as an isotropic function of the Lagrangian strain and texture tensors. The cases of isotropy and transverse isotropy are derived from the orthotropic case. Formulations are provided for both material and spatial frames. These formulations are suitable for implementation into finite element codes. It is also shown that the strain energy function can be naturally uncoupled into a dilatational and a distortional part, to facilitate the computational implementation of incompressibility.     

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.     

A model of fluid-biofilm interaction using a Burger material law

A two-dimensional finite element model of the biofilm response to flow was developed. The numerical code sequentially coupled the fluid dynamics of turbulent, incompressible flow with the mechanical response of a single hemispherical biofilm cluster (approximately 100 microm) attached to the flow boundary. A non-linear Burger material law was used to represent the viscoelastic response of a representative microbial biofilm. This constitutive law was incorporated into the numerical model as a Prony series representation of the biofilm's relaxation modulus. Model simulations illuminated interesting details of this fluid-structure interaction. Simulations revealed that softer biofilms (characterized by lower elastic moduli) were highly susceptible to lift forces and consequently were subject to even greater drag forces found higher in the velocity field. A bimodal deformation path due to the two Burger relaxation times was also observed in several simulations. This suggested that interfacial biofilm may be most susceptible to hydrodynamically induced detachment during the initial relaxation time. This result may prove useful in developing removal strategies. Additionally, plots of lift versus drag suggested that the deformation paths taken by viscoelastic biofilms are largely insensitive to specific material coefficients. Softer biofilms merely seem to follow the same path (as a stiffer biofilm) at a faster rate. These relationships may be useful in estimating the hydrodynamic forces acting on an attached biofilm based on changes in scale and cataloged material properties.     

A finite element study of invariant-based orthotropic constitutive equations in the context of myocardial material parameter estimation

A previous study investigated a number of invariant-based orthotropic and transversely isotropic constitutive equations for their suitability to fit three-dimensional simple shear mechanics data of passive myocardial tissue. The study was based on the assumption of a homogeneous deformation. Here, we extend the previous study by performing an inverse finite element material parameter estimation. This ensures a more realistic deformation state and material parameter estimates. The constitutive relations were compared on the basis of (i) ‘goodness of fit’: how well they fit a set of six shear deformation tests and (ii) ‘variability’: how well determined the material parameters are over the range of experiments. These criteria were utilised to discuss the advantages and disadvantages of the constitutive relations. It was found that a specific form of the polyconvex type as well as the exponential Fung-type equations were most suitable for modelling the orthotropic behaviour of myocardium under simple shear.