On the importance of geometric nonlinearity in finite-element simulations of trabecular bone failure.

The finite element method, which has been successfully applied to studies of the elastic properties of trabecular bone, is now being used to simulate its failure. These simulations have used a geometrically linear (linear kinematic) approximation to the total stiffness matrix; nonlinear terms in the total stiffness matrix have been excluded from the computation in order to achieve efficiency. Because trabecular bone appears to be a slender (i.e., geometrically nonlinear) structure, we studied the validity of the linear kinematic approximation for simulating its failure. Two cases, designed to bracket the extremes of stability behavior, were explored: a single representative spicule of trabecular bone (case 1) and a volume of trabecular bone consisting of relatively low aspect ratio members (case 2). For case 1, geometrically linear (GL) and nonlinear (GNL) analyses were performed with two different materials models: a plastic damage model and a brittle damage model. When GNL terms were included in the total stiffness matrix, we found that load-path bifurcation preceded tissue failure regardless of the form of the damage model. This bifurcation was the result of a complex coupling between material yield and structural instability. The nature of this coupling was highly sensitive to the form of the damage model. None of these behaviors was observed in the linear analyses, where failure was insensitive to the form of the damage model and where structural instabilities were prevented from occurring. For case 2, compressive loading of a volume of trabecular bone, geometric nonlinear effects were pronounced. There was a bifurcation in load response that resulted in large apparent strain to failure. The GL simulations, on the other hand, precluded this bifurcation. We hypothesize that trabecular bone is a geometric nonlinear structure; nonlinear terms must be included in the total stiffness matrix to accurately simulate its failure.

[1]  T. Hughes,et al.  An element-by-element solution algorithm for problems of structural and solid mechanics , 1983 .

[2]  Sally J. Marshall,et al.  The X-ray tomographic microscope: Three-dimensional perspectives of evolving microstructures , 1994 .

[3]  G. Niebur,et al.  High-resolution finite element models with tissue strength asymmetry accurately predict failure of trabecular bone. , 2000, Journal of biomechanics.

[4]  M. Swain,et al.  Fracture toughness of bovine bone: influence of orientation and storage media. , 2001, Biomaterials.

[5]  M. Ashby,et al.  Cellular solids: Structure & properties , 1988 .

[6]  L. Gibson,et al.  Modeling the mechanical behavior of vertebral trabecular bone: effects of age-related changes in microstructure. , 1997, Bone.

[7]  F. Eckstein,et al.  Estimation of distal radius failure load with micro-finite element analysis models based on three-dimensional peripheral quantitative computed tomography images. , 2002, Bone.

[8]  Z. Bažant,et al.  Stability Of Structures , 1991 .

[9]  Z. Bažant,et al.  Stability of Structures: Elastic, Inelastic, Fracture, and Damage Theories , 1993 .

[10]  Michael A. Puso,et al.  A highly efficient enhanced assumed strain physically stabilized hexahedral element , 2000 .

[11]  S. Goldstein,et al.  Finite‐element modeling of trabecular bone: Comparison with mechanical testing and determination of tissue modulus , 1998, Journal of orthopaedic research : official publication of the Orthopaedic Research Society.

[12]  John H. Kinney,et al.  Elastic constants of cellular structures , 1997 .

[13]  Lien-Wen Chen,et al.  Vibration and Buckling of Initially Stressed Curved Beams , 1998 .

[14]  van B Bert Rietbergen,et al.  Micro-FE analyses of bone: state of the art. , 2001 .

[15]  P Zioupos,et al.  Changes in the stiffness, strength, and toughness of human cortical bone with age. , 1998, Bone.

[16]  D. Fyhrie,et al.  Finite element calculated uniaxial apparent stiffness is a consistent predictor of uniaxial apparent strength in human vertebral cancellous bone tested with different boundary conditions. , 2001, Journal of biomechanics.

[17]  G. Niebur,et al.  Biomechanics of trabecular bone. , 2001, Annual review of biomedical engineering.

[18]  Whirley.  DYNA3D: A nonlinear, explicit, three-dimensional finite element code for solid and structural mechanics , 1993 .

[19]  P. Levitz,et al.  A new method for three-dimensional skeleton graph analysis of porous media: application to trabecular bone microarchitecture. , 2000, Journal of microscopy.

[20]  Hwj Rik Huiskes,et al.  Prediction of trabecular bone failure parameters using a tissue failure criterion and muFE analysis , 2000 .

[21]  Z. P. Bažant,et al.  Stability of Elastic, Anelastic, and Disintegrating Structures: A Conspectus of Main Results , 2000 .

[22]  J. Z. Zhu,et al.  The finite element method , 1977 .

[23]  D B Kimmel,et al.  Early Estrogen Replacement Therapy Reverses the Rapid Loss of Trabecular Bone Volume and Prevents Further Deterioration of Connectivity in the Rat , 1999, Journal of bone and mineral research : the official journal of the American Society for Bone and Mineral Research.

[24]  J. C. Simo,et al.  Strain- and stress-based continuum damage models—I. Formulation , 1987 .

[25]  C. M. Agrawal,et al.  Microstructural heterogeneity and the fracture toughness of bone. , 2000, Journal of biomedical materials research.

[26]  Stress-strain relations for elastomeric foams in uni-, bi- and tri-axial compression modes , 1996 .

[27]  D P Fyhrie,et al.  Failure mechanisms in human vertebral cancellous bone. , 1994, Bone.

[28]  G. Niebur,et al.  Sensitivity of damage predictions to tissue level yield properties and apparent loading conditions. , 2001, Journal of biomechanics.

[29]  F. Eckstein,et al.  Prediction of distal radius failure load with µFE models based on 3D-pQCT scans , 2001 .

[30]  S. Majumdar,et al.  High-resolution MRI and micro-FE for the evaluation of changes in bone mechanical properties during longitudinal clinical trials: application to calcaneal bone in postmenopausal women after one year of idoxifene treatment. , 2002, Clinical biomechanics.

[31]  X. Guo,et al.  Mechanical consequence of trabecular bone loss and its treatment: a three-dimensional model simulation. , 2002, Bone.

[32]  D Vashishth,et al.  Bone stiffness predicts strength similarly for human vertebral cancellous bone in compression and for cortical bone in tension. , 2000, Bone.

[33]  Sanjay Govindjee,et al.  Anisotropic modelling and numerical simulation of brittle damage in concrete , 1995 .