Passive deformation analysis of human leukocytes

The following analysis presents an experimental and theoretical study of the passive viscoelastic behavior of human leukocytes. Individual neutrophils in EDTA were observed both during their partial aspiration into a small micropipette and after expulsion from a large micropipette where the cell had been totally aspirated and deformed into a sausage shape. To analyze the data, a passive model of leukocyte rheology has been developed consisting of a cortical shell containing a Maxwell fluid which describes the average properties of the cell cytoplasm. The cortical shell represents a crosslinked actin layer near the surface of the cell and is assumed to be under pre-stressed tension. This model can reproduce the results of experiments using micropipette for both short-time small deformation and slow recovery data after large deformation. In addition, a finite element scheme has been established for the same model which shows close agreement with the analytical solution.     

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.     

An orthotropic viscoelastic model for the passive myocardium: continuum basis and numerical treatment

This study deals with the viscoelastic constitutive modeling and the respective computational analysis of the human passive myocardium. We start by recapitulating the locally orthotropic inner structure of the human myocardial tissue and model the mechanical response through invariants and structure tensors associated with three orthonormal basis vectors. In accordance with recent experimental findings the ventricular myocardial tissue is assumed to be incompressible, thick-walled, orthotropic and viscoelastic. In particular, one spring element coupled with Maxwell elements in parallel endows the model with viscoelastic features such that four dashpots describe the viscous response due to matrix, fiber, sheet and fiber-sheet fragments. In order to alleviate the numerical obstacles, the strictly incompressible model is altered by decomposing the free-energy function into volumetric-isochoric elastic and isochoric-viscoelastic parts along with the multiplicative split of the deformation gradient which enables the three-field mixed finite element method. The crucial aspect of the viscoelastic formulation is linked to the rate equations of the viscous overstresses resulting from a 3-D analogy of a generalized 1-D Maxwell model. We provide algorithmic updates for second Piola–Kirchhoff stress and elasticity tensors. In the sequel, we address some numerical aspects of the constitutive model by applying it to elastic, cyclic and relaxation test data obtained from biaxial extension and triaxial shear tests whereby we assess the fitting capacity of the model. With the tissue parameters identified, we conduct (elastic and viscoelastic) finite element simulations for an ellipsoidal geometry retrieved from a human specimen.     

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.     

Dynamic finite element implementation of nonlinear,anisotropic hyperelastic biological membranes

We present a novel method for the implementation of hyperelastic finite strain, non-linear strain-energy functions for biological membranes in an explicit finite element environment. The technique is implemented in LS-DYNA but may also be implemented in any suitable non-linear explicit code. The constitutive equations are implemented on the foundation of a co-rotational uniformly reduced Hughes-Liu shell. This shell is based on an updated-Lagrangian formulation suitable for relating Cauchy stress to the rate-of-deformation, i.e. hypo-elasticity. To accommodate finite deformation hyper-elastic formulations, a co-rotational deformation gradient is assembled over time, resulting in a formulation suitable for pseudo-hyperelastic constitutive equations that are standard assumptions in biomechanics. Our method was validated by comparison with (1) an analytic solution to a spherically-symmetric dynamic membrane inflation problem, incorporating a Mooney-Rivlin hyperelastic equation and (2) with previously published finite element solutions to a non-linear transversely isotropic inflation problem. Finally, we implemented a transversely isotropic strain-energy function for mitral valve tissue. The method is simple and accurate and is believed to be generally useful for anyone who wishes to model biologic membranes with an experimentally driven strain-energy function.     

Toward a biomechanical tool to evaluate rupture potential of abdominal aortic aneurysm: identification of a finite strain constitutive model and evaluation of its applicability

Knowledge of the wall stresses in an abdominal aortic aneurysm (AAA) may be helpful in evaluating the need for surgical intervention to avoid rupture. This must be preceded by the development of a more suitable finite strain constitutive model for AAA, as none currently exists. Additionally, reliable stress analysis of in vivo AAA for the purposes of clinical diagnostics requires patient-specific values of the material parameters, which are difficult to determine noninvasively. The purpose of this work, therefore, was three-fold: (1) to develop a finite strain constitutive model for AAA; (2) to estimate the variation of model parameters within a sample population; and (3) to evaluate the sensitivity of computed stress distribution in AAA due to this biologic variation. We propose here a two parameter, hyperelastic, isotropic, incompressible material model and utilize experimental data from 69 freshly excised AAA specimens to both develop the functional form of the model and estimate its material parameters. Parametric analyses were performed via repeated finite element computations to determine the effect of varying each of the two model parameters on the stress distribution in a three-dimensional AAA model. The agreement between experimental data and the proposed functional form of the constitutive law was very good (R2 > 0.9). Our finite element simulations showed that the computed AAA wall stresses changed by only 4% or less when both the parameters were varied within the 95% confidence intervals for the patient population studied. This observation indicates that in lieu of the patient-specific material parameters, which are difficult to determine the use of population mean values is sufficiently accurate for the model to be reasonably employed in a clinical setting. We believe that this is an important advancement toward the development of a computational tool for the estimation of rupture potential for individual AAA, for which there is great clinical need.     

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相似文献   

Mechanical properties of brain tissue in tension

This paper contains experimental results of in vitro, uniaxial tension of swine brain tissue in finite deformation as well as proposes a new hyper-viscoelastic constitutive model for the brain tissue. The experimental results obtained for two loading velocities, corresponding to strain rates of 0.64 and 0.64 x 10(-2)s(-1), are presented. We believe that these are the first ever experiments of this kind. The applied strain rates were similar to those applied in our previous study, focused on explaining brain tissue properties in compression. The stress-strain curves are convex downward for all extension rates. The tissue response stiffened as the loading speed increased, indicating a strong stress-strain rate dependence. Swine brain tissue was found to be considerably softer in extension than in compression. Previously proposed in the literature brain tissue constitutive models, developed based on experimental data collected in compression are shown to be inadequate to explain tissue behaviour in tension. A new, non-linear, viscoelastic model based on the generalisation of the Ogden strain energy hyper-elastic constitutive equation is proposed. The new model accounts well for brain tissue deformation behaviour in both tension and compression (natural strain in <-0.3,0.2>) for strain rates ranging over five orders of magnitude.     

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.     

Theory of the Sphering of Red Blood Cells

A rigorous mathematical solution of the sphering of a red blood cell is obtained under the assumptions that the red cells is a fluid-filled shell and that it can swell into a perfect sphere in an appropriate hypotonic medium. The solution is valid for finite strain of the cell membrane provided that the membrane is isotropic, elastic and incompressible. The most general nonlinear elastic stress-strain law for the membrane in a state of generalized plane stress is used. A necessary condition for a red cell to be able to sphere is that its extensional stiffness follow a specific distribution over the membrane. This distribution is strongly influenced by the surface tension in the cell membrane. A unique relation exists between the extensional stiffness, pressure differential, surface tension, and the ratio of the radius of the sphere to that of the undeformed red cell. The functional dependence of this stiffness distribution on various physical parameters is presented. A critique of some current literature on red cell mechanics is presented.     

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.     

The influence of the material properties on the biomechanical behavior of the pelvic floor muscles during vaginal delivery

In this work, a finite element model intends to represent the effects that the passage of a fetal head can induce on the muscles of the pelvic floor, from a mechanical point of view.The finite element method is a valuable tool, that is contributing to the clarification of the mechanisms behind pelvic floor disorders related to vaginal deliveries, although some care is necessary in order to obtain correct results. The present work shows how the variation of the material parameters, used in the constitutive model, can affect the obtained results from a finite element simulation. The constitutive equation adopted in this work for the pelvic floor muscles is a modified form of the incompressible transversely isotropic hyperelastic model proposed earlier by Humphrey and Yin.Results for the pelvic floor strain and stresses obtained during the passage of the fetus head are presented. The results show the importance of the material parameters and the need for a correct constitutive model.     

A simplified model of scar contraction

The healing of wounds is a complex process and the contraction of the resulting scar can have a negative impact on the neighbouring skin. A finite element model of skin simulating the contraction of a scar and deformation of the surrounding skin is presented. The skin is represented by an orthotropic–viscoelastic constitutive law, which is validated against experimental data in the literature. A simplified experimental model of a contracting scar in real skin is also developed. The pattern and size of the wrinkles formed around the contracting scar in the finite element model compare favourably with those formed in the experimental model. The orthotropic nature of skin plays a significant role in the behaviour of skin around scars—the wrinkles have a preferential orientation that corresponds to a direction perpendicular to the Langer's lines in the skin. The pre-stress in skin (a property that is ignored in many skin models) is shown to be an important factor in wrinkle formation around scars. The proposed model can be used to analyse the suturing and closure of wounds of various shapes.     

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.     

Influence of power-law rheology on cell injury during microbubble flows

The reopening of fluid-occluded pulmonary airways generates microbubble flows which impart complex hydrodynamic stresses to the epithelial cells lining airway walls. In this study we used boundary element solutions and finite element techniques to investigate how cell rheology influences the deformation and injury of cells during microbubble flows. An optimized Prony–Dirichlet series was used to model the cells’ power-law rheology (PLR) and results were compared with a Maxwell fluid model. Results indicate that membrane strain and the risk for cell injury decreases with increasing channel height and bubble speed. In addition, the Maxwell and PLR models both indicate that increased viscous damping results in less cellular deformation/injury. However, only the PLR model was consistent with the experimental observation that cell injury is not a function of stress exposure duration. Correlation of our models with experimental observations therefore highlights the importance of using PLR in computational models of cell mechanics/deformation. These computational models also indicate that altering the cell’s viscoelastic properties may be a clinically relevant way to mitigate microbubble-induced cell injury.     

A model for the nonlinear elastic response of large arteries

The nonlinear elastic response of large arteries subjected to finite deformations due to action of biaxial principal stresses, is described by simple constitutive equations. Generalized measures of strain and stress are introduced to account for material nonlinearity. This also ensures the existence of a strain energy density function. The orthotropic elastic response is described via quasi-linear relations between strains and stresses. One nonlinear parameter which defines the measures of strain and stress, and three elastic moduli are assumed to be constants. The lateral strain parameters (equivalent to Poisson's ratios in infinitesimal deformations) are deformation dependent. This dependence is defined by empirical relations developed via the incompressibility condition, and by the introduction of a fifth material parameter. The resulting constitutive model compares well with biaxial experimental data of canine carotid arteries.     

A 3D active-passive numerical skeletal muscle model incorporating initial tissue strains. Validation with experimental results on rat tibialis anterior muscle

This paper presents a three-dimensional finite element model of skeletal muscle and its validation incorporating inital tissue strains. A constitutive relation was determined by using a convex free strain energy function (SEF) where active and passive response contributions were obtained fitting experimental data from the rat tibialis anterior (TA) muscle. The passive and active finite strains response was modelled within the framework of continuum mechanics by a quasi-incompressible transversely isotropic material formulation. Magnetic resonance images (MRI) were obtained to reconstruct the external geometry of the TA. This geometry includes initial strains also taken into account in the numerical model. The numerical results show excellent agreement with the experimental results when comparing reaction force-extension curves both in passive and active tests. The proposed constitutive model for the muscle is implemented in a subroutine in the commercial finite element software package ABAQUS.     

Wall shear stress in backward-facing step flow of a red blood cell suspension.

An experimental investigation of the wall shear stress distribution downstream of a backward-facing step is carried out. The wall shear stress distribution was determined by measuring the deformation of a gel layer, attached to the wall downstream of the step. Speckle pattern interferometry was applied to measure the deformation of the gel layer. The measured deformation, combined with the properties of the gel layer, served as an input for a finite element solid mechanics computation to determine the stress distribution in the gel layer. The wall shear stress, required to generate the measured deformation of the gel layer, was determined from these computations. A Newtonian buffer solution and a non-Newtonian red blood cell suspension were used as measuring fluids. The deformation of the gel layer was determined for a Newtonian buffer solution to evaluate the method and to obtain the properties of the gel layer. Subsequently, the wall shear stress distribution for the non-Newtonian red blood cell suspension was determined for three different flow rates. The inelastic non-Newtonian Carreau-Yasuda model served as constitutive model for the red blood cell suspension. Using this model, the velocity and wall shear stress distribution were computed by means of a finite element fluid mechanics computation. From the comparison between the numerical and the experimental results, it can be concluded that wall shear stresses, induced by the red blood cell suspension, can be modeled accurately by employing a Carreau-Yasuda model.     

Non-homogeneous analysis of three-dimensional transmural finite deformation in canine ventricular myocardium

A new method has been developed for analyzing transmural distributions of finite deformation in canine ventricular myocardium without the need to assume that the strain in a finite volume of the wall is homogeneous. The three-dimensional nodal geometric parameters of bilinear-cubic or bilinear-quadratic finite elements are fitted by least squares to the measured coordinates of 12-18 radiopaque markers implanted in the left ventricular free wall. For six dog hearts, root-mean-squared errors in the fitted in-plane coordinates ranged from 0.079-0.556 mm in the end-diastolic reference state and 0.142-0.622 mm at end-systole. The corresponding error ranges in the radial coordinate were 0.042-0.264 mm at end-diastole and 0.106-0.279 mm at end-systole. Smoothly continuous transmural profiles of wall strain computed as the element deformed during the cardiac cycle from end-diastole to end-systole showed good agreement with the discrete results of conventional homogeneous analysis. Using the kinematics of a thick-walled incompressible cylinder, overall absolute errors due to the non-homogeneity of myocardial deformation were found to be reduced in the new analysis by 30-35% for typical experimental parameters. Overall relative errors were also reduced (from 23 to 20%). Since measurement errors in the reconstructed marker coordinates were spatially smoothed by the fitting procedure, noise in the computed deformations was also substantially attenuated, and transmural gradients of three-dimensional strain components could be obtained with improved accuracy. Hence physiological factors affected by transmural stress and strain distributions, such as myocardial blood flow, ischemia and hypertrophy, may be better understood.     

An anisotropic elastic-viscoplastic damage model for bone tissue

A new anisotropic elastic-viscoplastic damage constitutive model for bone is proposed using an eccentric elliptical yield criterion and nonlinear isotropic hardening. A micromechanics-based multiscale homogenization scheme proposed by Reisinger et al. is used to obtain the effective elastic properties of lamellar bone. The dissipative process in bone is modeled as viscoplastic deformation coupled to damage. The model is based on an orthotropic ecuntric elliptical criterion in stress space. In order to simplify material identification, an eccentric elliptical isotropic yield surface was defined in strain space, which is transformed to a stress-based criterion by means of the damaged compliance tensor. Viscoplasticity is implemented by means of the continuous Perzyna formulation. Damage is modeled by a scalar function of the accumulated plastic strain ${D(\kappa)}$ , reducing all element s of the stiffness matrix. A polynomial flow rule is proposed in order to capture the rate-dependent post-yield behavior of lamellar bone. A numerical algorithm to perform the back projection on the rate-dependent yield surface has been developed and implemented in the commercial finite element solver Abaqus/Standard as a user subroutine UMAT. A consistent tangent operator has been derived and implemented in order to ensure quadratic convergence. Correct implementation of the algorithm, convergence, and accuracy of the tangent operator was tested by means of strain- and stress-based single element tests. A finite element simulation of nano- indentation in lamellar bone was finally performed in order to show the abilities of the newly developed constitutive model.