Characterisation of Trabecular Bone Behaviour Under Impact

Increasing ageing population around the world implies deteriorated bone quality due to diseases such as osteoporosis. Osteoporosis results in reduced bone density and strength and increases the risk of bone fractures due to impact loads such as those caused by a fall. While numerous studies have been undertaken to evaluate the mechanical properties of bone using experimental or computational approaches most of these have focused on static or cyclic loads. Little is known about the behaviour of bone under dynamic loading, such as that caused by a fall. This study focuses on characterising the impact behaviour of spongy bone which is also known as trabecular or cancellous bone using drop tower tests.

A closed-form solution of drop test on an elastic material is developed in detail to serve as a benchmark. Sequential compressive stress wave equations are derived then coded in programming language to obtain desired number of equations, leading to the resulting force history at the end struck. The number of required intervals to obtain the whole pulse depends on the mass ratio; the developed code can be readily expanded to include a wide range of these ratios. A parametric analysis is undertaken to evaluate the influence of parameters such as elastic modulus and drop height. This evaluation defines the key outputs of the force response, namely peak load, pulse width, interval frequency, number of stress intervals and loading rate; hence it establishes the output variable-parameter relationships, e.g. peak load-elastic modulus relation.

The force response from the closed-form solution is used to verify the finite element (FE) models. The Poisson’s ratio is found to have insignificant effect on the resulting peak load and pulse width. Stiffer response is indicated by an increase in peak load or a decrease in pulse width for the same impact. Fast Fourier transform (FFT) is applied to the pulses obtained from the closed-form solution and FE analysis; the returned frequency values reveal that the discrepancy observed in FE results can be attributed to its solver rather than the FFT technique. The time-independent elastic FE model is extended to include the time-dependent viscoelastic behaviour; the latter provides a stiffer response. The output variable-parameter relations from FE analysis are verified by similar expressions derived from the closed-form solution.

Drop tests are conducted in the lab on trabecular bone samples with varying bone volume fraction (bone volume to total volume ratio or BV/TV) to investigate the resulting force response. These responses are categorised on the basis of BV/TV of the tested samples. Samples with large BV/TV are shown to behave in an apparently elastic manner. Low BV/TV samples respond inelastically which the study attempts to simulate using plasticity. Behaviour of samples with growth plates shows inelastic behaviour even though they have large BV/TV.

Regression analysis is carried out on the outputs of pulses obtained from higher BV/TV samples to find the output variable-parameter relations, e.g. inverse relation of pulse width-BV/TV. Inverse modelling using FE analysis is also performed to estimate a representable elastic modulus by matching the peak load and pulse width of the resulting force response. The elastic moduli are found to be in the range of 650 to 1400 MPa for samples with BV/TV ranging between 35 to 53%. Elastic modulus is also evaluated using the initial load-response curve for both lower and higher BV/TV samples to establish linear and power-law loading rate versus BV/TV relationships.

The inelastic behaviour of low BV/TV samples consistently shows initial peak, followed by a drop to a finite non-zero value in the post-elastic regime which is maintained for a considerable prolonged time before the load returns to zero. This pattern has not been previously demonstrated. The duration of plasticity is found to be larger for lower BV/TV samples, while shorter drop time is observed exhibiting brittle breakage after the peak load is attained. The study also considers simulation of the post-yield inelastic response using strain-softening plasticity. It is shown that strain-softening is capable of approximately replicating the post-yield impact response pattern by assuming appropriate value of tangent modulus.

Four dynamic elastic modulus expressions as functions of BV/TV are derived by correlating the relationships regressed experimentally with their correspondences from the closed-form solution. The performance of these linear and power law equations are compared to several relations from the literature and the response from tests conducted at an apparent strain rate of 0.01 /sec on similar samples. It is found that the dynamic elastic moduli found in this study are higher than their quasi-static and monotonic loading counterparts; a conclusion showing the effect of strain rate magnitude on the stiffness of the trabecular bone. The study evaluates the strain rate experienced by bone due to impact loading considered in this study to be a maximum of 44.3 /sec at the initial stage with secant value of 30 /sec at the peak strain. Assuming the high BV/TV samples behave fully elastic, the apparent strain due to impact loads can be as high as 3% in compression.

Year	Entry
2014	Gilchrist S, Nishiyama KK, de Bakker P, Guy P, Boyd SK, Oxland T, Cripton PA. Proximal femur elastic behaviour is the same in impact and constant displacement rate fall simulation. J Biomech. November 28, 2014;47(15):3744-3749.
1985	Gibson LJ. The mechanical behaviour of cancellous bone. J Biomech. 1985;18(5):317-328.
2001	Keyak JH. Improved prediction of proximal femoral fracture load using nonlinear finite element models. Med Eng Phys. April 2001;23(3):165-173.
2004	Bayraktar HH, Keaveny TM. Mechanisms of uniformity of yield strains for trabecular bone. J Biomech. November 2004;37(11):1671-1678.
1997	Li B, Aspden RM. Composition and mechanical properties of cancellous bone from the femoral head of patients with osteoporosis or osteoarthritis. J Bone Miner Res. 1997;12(4):641-651.
2011	Trabelsi N, Yosibash Z, Wutte C, Augat P, Eberle S. Patient-specific finite element analysis of the human femur: a double-blinded biomechanical validation. J Biomech. June 3, 2011;44(9):1666-1672.
2000	Ciarelli TE, Fyhrie DP, Schaffler MB, Goldstein SA. Variations in three‐dimensional cancellous bone architecture of the proximal femur in female hip fractures and in controls. J Bone Miner Res. 2000;15(1):32-40.
1994	Keyak JH, Lee IY, Skinner HB. Correlations between orthogonal mechanical properties and density of trabecular bone: use of different densitometric measures. J Biomed Mater Res. November 1994;A28(11):1329-1336.
2009	de Bakker PM, Manske SL, Ebacher V, Oxland TR, Cripton PA, Guy P. During sideways falls proximal femur fractures initiate in the superolateral cortex: evidence from high-speed video of simulated fractures. J Biomech. August 25, 2009;42(12):1917-1925.
2001	Turner CH, Burr DB. Experimental techniques for bone mechanics. In: Cowin SC, ed. Bone Mechanics Handbook. 2nd ed. Boca Raton, FL: CRC Press; 2001:chap 7.
1994	Goulet RW, Goldstein SA, Ciarelli MJ, Kuhn JL, Brown MB, Feldkamp LA. The relationship between the structural and orthogonal compressive properties of trabecular bone. J Biomech. April 1994;27(4):375-389.
1998	Keyak JH, Rossi SA, Jones KA, Skinner HB. Prediction of femoral fracture load using automated finite element modeling. J Biomech. February 1998;31(2):125-133.
2008	Schileo E, Taddei F, Cristofolini L, Viceconti M. Subject-specific finite element models implementing a maximum principal strain criterion are able to estimate failure risk and fracture location on human femurs tested in vitro. J Biomech. 2008;41(2):356-367.
2008	Charlebois M. A Constitutive Law for Softening and Densification of Trabecular Bone in Compression [PhD thesis]. Vienna University of Technology; 2008.
2013	Gilchrist S, Guy P, Cripton PA. Development of an inertia-driven model of sideways fall for detailed study of femur fracture mechanics. J Biomech Eng. December 2013;135(12):121001.
1979	Lakes RS, Katz JL, Sternstein SS. Viscoelastic properties of wet cortical bone, I: torsional and biaxial studies. J Biomech. 1979;12(9):657-678.
2005	Shim VPW, Yang LM, Liu JF, Lee VS. Characterisation of the dynamic compressive mechanical properties of cancellous bone from the human cervical spine. Int J Impact Eng. December 2005;32(1-4):525-540.
1998	Kopperdahl DL, Keaveny TM. Yield strain behavior of trabecular bone. J Biomech. July 1998;31(7):601-608.
1997	Ouyang J, Yang GT, Wu WZ, Zhu QA, Zhong SZ. Biomechanical characteristics of human trabecular bone. Clin Biomech (Bristol, Avon). October 1997;12(7-8):522-524.
2018	Xie S. Characterisation of Time-Dependent Mechanical Behaviour of Trabecular Bone and Its Constituents [PhD thesis]. Edinburgh, Scotland: University of Edinburgh; 2018.
1977	Carter DR, Hayes WC. The compressive behavior of bone as a two-phase porous structure. J Bone Joint Surg. 1977;59A(7):954-962.
1994	Keller TS. Predicting the compressive mechanical behavior of bone. J Biomech. September 1994;27(9):1159-1168.
2004	Bayraktar HH, Morgan EF, Niebur GL, Morris GE, Wong EK, Keaveny TM. Comparison of the elastic and yield properties of human femoral trabecular and cortical bone tissue. J Biomech. January 2004;37(1):27-35.
1993	Hayes WC, Myers ER, Morris JN, Gerhart TN, Yett HS, Lipsitz LA. Impact near the hip dominates fracture risk in elderly nursing home residents who fall. Calcif Tiss Int. March 1993;52(3):192-198.
2006	Taddei F, Cristofolini L, Martelli S, Gill HS, Viceconti M. Subject-specific finite element models of long bones: an in vitro evaluation of the overall accuracy. J Biomech. 2006;39(13):2457-2467.
2009	Bevill G, Keaveny TM. Trabecular bone strength predictions using finite element analysis of micro-scale images at limited spatial resolution. Bone. April 2009;44(4):579-584.
2018	Enns-Bray WS, Ferguson SJ, Helgason B. Strain rate dependency of bovine trabecular bone under impact loading at sideways fall velocity. J Biomech. June 2018;75:46-52.
2001	Heller MO, Bergmann G, Deuretzbacher G, Dürselen L, Pohl M, Claes L, Haas NP, Duda GN. Musculo-skeletal loading conditions at the hip during walking and stair climbing. J Biomech. July 2001;34(7):883-893.
2007	Schileo E, Taddei F, Malandrino A, Cristofolini L, Viceconti M. Subject-specific finite element models can accurately predict strain levels in long bones. J Biomech. 2007;40(13):2982-2989.
1997	Gibson LJ, Ashby MF. Cellular Solids: Structure and Properties. 2nd ed. Cambridge, UK: Cambridge University Press; 1997. Cambridge Solid State Science Series.
2013	Nishiyama KK, Gilchrist S, Guy P, Cripton P, Boyd SK. Proximal femur bone strength estimated by a computationally fast finite element analysis in a sideways fall configuration. J Biomech. April 26, 2013;46(7):1231-1236.
1986	Park HC, Lakes RS. Cosserat micromechanics of human bone: strain redistribution by a hydration sensitive constituent. J Biomech. 1986;19(5):385-397.
1992	Linde F, Hvid I, Madsen F. The effect of specimen geometry on the mechanical behaviour of trabecular bone specimens. J Biomech. 1992;25(4):359-368.
1998	Zysset PK, Goulet RW, Hollister SJ. A global relationship between trabecular bone morphology and homogenized elastic properties. J Biomech Eng. October 1998;120(5):640-646.
1998	Bowman SM, Guo XE, Cheng DW, Keaveny TM, Gibson LJ, Hayes WC, McMahon TA. Creep contributes to the fatigue behavior of bovine trabecular bone. J Biomech Eng. October 1998;120(5):647-654.
1990	Lotz JC, Gerhart TN, Hayes WC. Mechanical properties of trabecular bone from the proximal femur: a quantitative CT study. J Comput Assist Tomogr. January–February 1990;14(1):107-114.
2006	Sambrook P, Cooper C. Osteoporosis. Lancet. June 17, 2006;367(9527):2010-2018.
2001	Heiner AD, Brown TD. Structural properties of a new design of composite replicate femurs and tibias. J Biomech. January 2001;34(6):773-781.
2015	Zani L, Erani P, Grassi L, Taddei F, Cristofolini L. Strain distribution in the proximal human femur during in vitro simulated sideways fall. J Biomech. July 16, 2015;48(10):2130-2143.
2008	Hansen U, Zioupos P, Simpson R, Currey JD, Hynd D. The effect of strain rate on the mechanical properties of human cortical bone. J Biomech Eng. February 2008;130(1):011011.
2007	Cristofolini L, Juszczyk M, Martelli S, Taddei F, Viceconti M. In vitro replication of spontaneous fractures of the proximal human femur. J Biomech. 2007;40(13):2837-2845.
2004	Kaneko TS, Bell JS, Pejcic MR, Tehranzadeh J, Keyak JH. Mechanical properties, density and quantitative CT scan data of trabecular bone with and without metastases. J Biomech. April 2004;37(4):523-530.
2018	Xie S, Wallace RJ, Callanan A, Pankaj P. From tension to compression: asymmetric mechanical behaviour of trabecular bone’s organic phase. Ann Biomed Eng. June 2018;46(6):801-809.
1995	van Rietbergen B, Weinans H, Huiskes R, Odgaard A. A new method to determine trabecular bone elastic properties and loading using micromechanical finite-element models. J Biomech. January 1995;28(1):69-81.
2008	Helgason B, Perilli E, Schileo E, Taddei F, Brynjólfsson S, Viceconti M. Mathematical relationships between bone density and mechanical properties: a literature review. Clin Biomech (Bristol, Avon). 2008;23(2):135-146.
1998	Ulrich D, van Rietbergen B, Weinans H, Rüegsegger P. Finite element analysis of trabecular bone structure: a comparison of image-based meshing techniques. J Biomech. December 1998;31(12):1187-1192.
1999	Cody DD, Gross GJ, Hou FJ, Spencer HJ, Goldstein SA, Fyhrie DP. Femoral strength is better predicted by finite element models than QCT and DXA. J Biomech. October 1999;32(10):1013-1020.
2005	Mayhew PM, Thomas CD, Clement JG, Loveridge N, Beck TJ, Bonfield W, Burgoyne CJ, Reeve J. Relation between age, femoral neck cortical stability, and hip fracture risk. Lancet. July 9, 2005;366(9480):129-135.
2001	Morgan EF, Keaveny TM. Dependence of yield strain of human trabecular bone on anatomic site. J Biomech. 2001;34(5):569-577.
2003	Morgan EF, Bayraktar HH, Keaveny TM. Trabecular bone modulus–density relationships depend on anatomic site. J Biomech. July 2003;36(7):897-904.
1996	Cristofolini L, Viceconti M, Cappello A, Toni A. Mechanical validation of whole bone composite femur models. J Biomech. April 1996;29(4):525-535.
2007	Bessho M, Ohnishi I, Matsuyama J, Matsumoto T, Imai K, Nakamura K. Prediction of strength and strain of the proximal femur by a CT-based finite element method. J Biomech. 2007;40(8):1745-1753.
2018	Levrero-Florencio F, Pankaj P. Using non-linear homogenization to improve the performance of macroscopic damage models of trabecular bone. Front Physiol. May 2018;9:545.
1991	Linde F, Nørgaard P, Hvid I, Odgaard A, Søballe K. Mechanical properties of trabecular bone: dependency on strain rate. J Biomech. 1991;24(9):803-809.
2001	Lakes R. Viscoelastic properties of cortical bone. In: Cowin SC, ed. Bone Mechanics Handbook. 2nd ed. Boca Raton, FL: CRC Press; 2001:chap 11.
1999	Turner CH, Rho J, Takano Y, Tsui TY, Pharr GM. The elastic properties of trabecular and cortical bone tissues are similar: results from two microscopic measurement techniques. J Biomech. April 1999;32(4):437-441.
2012	Grassi L, Schileo E, Taddei F, Zani L, Juszczyk M, Cristofolini L, Viceconti M. Accuracy of finite element predictions in sideways load configurations for the proximal human femur. J Biomech. January 10, 2012;45(2):394-399.
1993	Keaveny TM, Borchers RE, Gibson L, Hayes WC. Trabecular bone modulus and strength can depend on specimen geometry. J Biomech. August 1993;26(8):991-995.

Year

Entry

2014

Gilchrist S, Nishiyama KK, de Bakker P, Guy P, Boyd SK, Oxland T, Cripton PA. Proximal femur elastic behaviour is the same in impact and constant displacement rate fall simulation. J Biomech. November 28, 2014;47(15):3744-3749.

1985

Gibson LJ. The mechanical behaviour of cancellous bone. J Biomech. 1985;18(5):317-328.

2001

Keyak JH. Improved prediction of proximal femoral fracture load using nonlinear finite element models. Med Eng Phys. April 2001;23(3):165-173.

2004

Bayraktar HH, Keaveny TM. Mechanisms of uniformity of yield strains for trabecular bone. J Biomech. November 2004;37(11):1671-1678.

1997

Li B, Aspden RM. Composition and mechanical properties of cancellous bone from the femoral head of patients with osteoporosis or osteoarthritis. J Bone Miner Res. 1997;12(4):641-651.