A Viscoelastic Modeling for Failure Analysis of Human Vertebral Bone Undergoing Multi-Rate Compression

The calibrated model reproduced the overall trend of specimen-to-specimen mechanical variation observed experimentally. Predicted stiffness values were in reasonable agreement with measured data. Fracture force predictions showed moderate agreement for dynamically tested specimens (R² = 0.60), which improved to R² = 0.88 after exclusion of one statistically identified outlier. Compared with a purely linear elastic formulation, the proposed viscoelastic model demonstrated modest improvement in stiffness prediction and more substantial improvement in fracture force prediction. These findings indicate that incorporating density-dependent viscoelastic effects improves representation of vertebral mechanical behavior, particularly at higher loading rates. Owing to its simplicity and computational efficiency, the proposed model requires only limited imaging input and may be useful for future biomechanical investigations, rapid screening, and injury risk prediction.