A nonlocal constitutive model for trabecular bone softening in compression

Using the three-dimensional morphological data provided by computed tomography, finite element (FE) models can be generated and used to compute the stiffness and strength of whole bones. Three-dimensional constitutive laws capturing the main features of bone mechanical behavior can be developed and implemented into FE software to enable simulations on complex bone structures. For this purpose, a constitutive law is proposed, which captures the compressive behavior of trabecular bone as a porous material with accumulation of irreversible strain and loss of stiffness beyond its yield point and softening beyond its ultimate point. To account for these features, a constitutive law based on damage coupled with hardening anisotropic elastoplasticity is formulated using density and fabric-based tensors. To prevent mesh dependence of the solution, a nonlocal averaging technique is adopted. The law has been implemented into a FE software and some simple simulations are first presented to illustrate its behavior. Finally, examples dealing with compression of vertebral bodies clearly show the impact of softening on the localization of the inelastic process.

[1]  Steven K Boyd,et al.  Bone strength at the distal radius can be estimated from high-resolution peripheral quantitative computed tomography and the finite element method. , 2008, Bone.

[2]  A. Curnier,et al.  An objective anisotropic elastic plastic model and algorithm applicable to bone mechanics , 1991 .

[3]  Howard L. Schreyer,et al.  A thermodynamically consistent framework for theories of elastoplasticity coupled with damage , 1994 .

[4]  A. Curnier,et al.  A three-dimensional elastic plastic damage constitutive law for bone tissue , 2009, Biomechanics and modeling in mechanobiology.

[5]  P J Prendergast,et al.  Finite element models in tissue mechanics and orthopaedic implant design. , 1997, Clinical biomechanics.

[6]  Bert Van Rietbergen,et al.  Finite Element Analysis Based on In Vivo HR‐pQCT Images of the Distal Radius Is Associated With Wrist Fracture in Postmenopausal Women , 2007, Journal of bone and mineral research : the official journal of the American Society for Bone and Mineral Research.

[7]  Milan Jirásek,et al.  Consistent tangent stiffness for nonlocal damage models , 2002 .

[8]  T. Keaveny,et al.  Cortical and Trabecular Load Sharing in the Human Vertebral Body , 2005, Journal of bone and mineral research : the official journal of the American Society for Bone and Mineral Research.

[9]  R. Ogden,et al.  Mechanics of biological tissue , 2006 .

[10]  Philippe K. Zysset,et al.  Multi-axial mechanical properties of human trabecular bone , 2009, Biomechanics and modeling in mechanobiology.

[11]  D. Davy,et al.  Inelastic strain accumulation in cortical bone during rapid transient tensile loading. , 1999, Journal of biomechanical engineering.

[12]  M. Ashby,et al.  Cellular solids: Structure & properties , 1988 .

[13]  Dieter H. Pahr,et al.  ANATOMY- AND MORPHOLOGY-BASED SMOOTH FINITE ELEMENT MODELS OF THE VERTEBRAL BODY , 2008 .

[14]  T. Keaveny,et al.  Effects of Teriparatide and Alendronate on Vertebral Strength as Assessed by Finite Element Modeling of QCT Scans in Women With Osteoporosis , 2006, Journal of bone and mineral research : the official journal of the American Society for Bone and Mineral Research.

[15]  Ralph Müller,et al.  Contribution of In Vivo Structural Measurements and Load/Strength Ratios to the Determination of Forearm Fracture Risk in Postmenopausal Women , 2007, Journal of bone and mineral research : the official journal of the American Society for Bone and Mineral Research.

[16]  D T Davy,et al.  A damage model for nonlinear tensile behavior of cortical bone. , 1999, Journal of biomechanical engineering.

[17]  M. Jirásek,et al.  Plastic model with non‐local damage applied to concrete , 2006 .

[18]  Peter Varga,et al.  A patient-specific finite element methodology to predict damage accumulation in vertebral bodies under axial compression, sagittal flexion and combined loads , 2008, Computer methods in biomechanics and biomedical engineering.

[19]  J. C. Simo,et al.  Consistent tangent operators for rate-independent elastoplasticity☆ , 1985 .

[20]  J. Chaboche,et al.  Mechanics of Solid Materials , 1990 .

[21]  Milan Jirásek,et al.  Nonlocal integral formulations of plasticity and damage : Survey of progress , 2002 .

[22]  Z. Bažant,et al.  Stability Of Structures , 1991 .

[23]  D. Krajcinovic,et al.  Some fundamental issues in rate theory of damage-elastoplasticity , 1995 .

[24]  Felix Eckstein,et al.  The role of fabric in the quasi-static compressive mechanical properties of human trabecular bone from various anatomical locations , 2008, Biomechanics and modeling in mechanobiology.

[25]  Bořek Patzák,et al.  Design of object oriented finite element code , 2001 .

[26]  G. Maugin.  The Thermomechanics of Plasticity and Fracture , 1992 .

[27]  Philippe K Zysset,et al.  A review of morphology-elasticity relationships in human trabecular bone: theories and experiments. , 2003, Journal of biomechanics.

[28]  Marco Viceconti,et al.  In vitro replication of spontaneous fractures of the proximal human femur. , 2007, Journal of biomechanics.

[29]  P. M. Naghdi,et al.  A general theory of an elastic-plastic continuum , 1965 .

[30]  Dieter H. Pahr,et al.  From high-resolution CT data to finite element models: development of an integrated modular framework , 2009 .

[31]  F. Eckstein,et al.  Estimation of distal radius failure load with micro-finite element analysis models based on three-dimensional peripheral quantitative computed tomography images. , 2002, Bone.

[32]  P. Zysset,et al.  An Alternative Fabric-based Yield and Failure Criterion for Trabecular Bone , 2006 .

[33]  Philippe Zysset.  A constitutive law for trabecular bone , 1994 .

[34]  Arturo N Natali,et al.  Constitutive modelling of inelastic behaviour of cortical bone. , 2008, Medical engineering & physics.

[35]  T. Keaveny,et al.  Finite element models predict in vitro vertebral body compressive strength better than quantitative computed tomography. , 2003, Bone.

[36]  Qi-Chang He,et al.  Conewise linear elastic materials , 1994 .

[37]  D. Krajcinovic,et al.  Introduction to continuum damage mechanics , 1986 .

[38]  Yan Chevalier,et al.  Cement Distribution, Volume, and Compliance in Vertebroplasty: Some Answers From an Anatomy-Based Nonlinear Finite Element Study , 2008, Spine.

[39]  Tony M Keaveny,et al.  Locations of bone tissue at high risk of initial failure during compressive loading of the human vertebral body. , 2007, Bone.

[40]  Kozo Nakamura,et al.  Assessment of vertebral fracture risk and therapeutic effects of alendronate in postmenopausal women using a quantitative computed tomography-based nonlinear finite element method , 2009, Osteoporosis International.

[41]  Philippe K. Zysset,et al.  An alternative model for anisotropic elasticity based on fabric tensors , 1995 .

[42]  G. Niebur,et al.  Biomechanics of trabecular bone. , 2001, Annual review of biomedical engineering.

[43]  R Huiskes,et al.  Indirect determination of trabecular bone effective tissue failure properties using micro-finite element simulations. , 2008, Journal of biomechanics.

[44]  W C Hayes,et al.  Postyield behavior of subchondral trabecular bone. , 1976, Journal of biomedical materials research.

[45]  Milan Jirásek,et al.  Damage-plastic model for concrete failure , 2006 .

[46]  A. Curnier,et al.  A 3D damage model for trabecular bone based on fabric tensors. , 1996, Journal of biomechanics.

[47]  Arturo N Natali,et al.  Investigation of bone inelastic response in interaction phenomena with dental implants. , 2008, Dental materials : official publication of the Academy of Dental Materials.