A new one-dimensional constitutive model for human cortical bone is proposed to simulate the viscoelastic–viscoplastic behaviour occurring during creep-recovery tests. The material parameters are determined by fitting experimental results of creep-recovery tests reported in the published literature. An efficient computational algorithm for the integration of the proposed constitutive model at the material point level is derived. The derived algorithm in conjunction with the Jacobian 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.