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. U radu je opisana numerička implementacija i validacija konstitutivnog modela za simuliranje mehaničkog ponašanja ljudskog kortikalnog koštanog tkiva. Model spreže viskoelastičnost s oštećenjem za simulaciju puzanja i puzanja s naknadnim rasterećenjem. Parametri materijala određeni su aproksimacijom eksperimentalnih rezultata iz literature. Izveden je računalni algoritam za integriranje predloženog konstitutivnog modela na razini materijalne točke. Izvedeni algoritam zajedno s tangentnom matricom krutosti implementiran je u program ABAQUS koji se temelji na metodi konačnih elemenata. Dobiveno je dobro poklapanje rezultata računalne simulacije i eksperimentalnih podataka iz literature.