Indentation tests can be used to characterize the mechanical properties of bone at small load/length scales offering the possibility of utilizing very small test specimens, which can be excised using minimally-invasive procedures. In addition, the need for mechanical property data from bone may be a requirement for fundamental multi-scale experiments, changes in nano- and micro-mechanical properties (e.g., as affected by changes in bone mineral density) due to drug therapies, and/or the development of computational models. Load vs. indentation depth data, however, is more complex than those obtained from typical macro-scale experiments, primarily due to the mixed state of stress, and thus interpretation of the data and extraction of mechanical properties is more challenging. Previous studies have shown that cortical bone exhibits a visco-elastic response combined with permanent deformation during indentation tests, and that the load vs. indentation depth response can be simulated using a visco-elastic/plastic material model. The model successfully captures the loading and creep displacement behavior, however, it does not adequately reproduce the unloading response near the end of the unloading cycle, where a pronounced decrease in contact stiffness is observed. It is proposed that the stiffness reduction observed in bone results from an increase in damage;​ therefore, a plastic-damage model was investigated and shown capable of simulating a typical bone indentation response through an axisymmetric finite element simulation. The plastic-damage model was able to reproduce the full indentation response, especially the reduced stiffness behavior exhibited during the latter stages of unloading. The results suggest that the plastic-damage model is suitable for describing the complex indentation response of bone and may provide further insight into the relationship between model parameters and mechanical/physical properties.