An anisotropic elastic-viscoplastic damage model for bone tissue

• Published in: Biomechanics and modeling in mechanobiology Vol. 12; no. 2; pp. 201 - 213
• Main Authors: Schwiedrzik, J. J., Zysset, P. K.
• Format: Journal Article
• Language: English
• Published: Berlin/Heidelberg Springer-Verlag 01-04-2013; Springer Nature B.V
• Subjects: Algorithms; Anisotropy; Biological and Medical Physics; Biomechanics; Biomedical Engineering and Bioengineering; Biophysics; Bone and Bones - pathology; Elastic properties; Elasticity; Engineering; Finite Element Analysis; Models, Biological; Original Paper; Reproducibility of Results; Rheology; Space life sciences; Strain; Stress, Mechanical; Theoretical and Applied Mechanics; Tissue engineering; Viscosity; Constitutive model; Bone; Damage; Viscoplasticity
• ISSN: 1617-7959; 1617-7940
• DOI: 10.1007/s10237-012-0392-9

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
, 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.