This paper describes the numerical implementation and validation of a constitutive model for simulating the mechanical behaviour of human cortical bone tissue. The model incorporates linear viscoelasticity coupled with damage to predict the creep and creep-recovery responses, respectively. The material parameters are determined by fitting the experimental results reported in the published literature. A computational algorithm for the integration of the proposed constitutive model at the material point level is derived. The derived algorithm in conjunction with the tangent stiffness matrix is implemented in the finite element code ABAQUS. The model predictions are found to be in good agreement with the experimental data presented in literature.
展开▼