0
Research Papers

The Quartic Piecewise-Linear Criterion for the Multiaxial Yield Behavior of Human Trabecular Bone

[+] Author and Article Information
Arnav Sanyal

Department of Mechanical Engineering,
Orthopaedic Biomechanics Laboratory,
University of California,
Berkeley, CA 94720
e-mail: arnavsanyal@gmail.com

Joanna Scheffelin

Department of Mechanical Engineering,
Orthopaedic Biomechanics Laboratory,
University of California,
Berkeley, CA 94720

Tony M. Keaveny

Department of Mechanical Engineering,
Orthopaedic Biomechanics Laboratory,
University of California,
Berkeley, CA 94720
Department of Bioengineering,
University of California,
Berkeley, CA 94720

Manuscript received August 2, 2014; final manuscript received November 10, 2014; accepted manuscript posted November 17, 2014; published online December 10, 2014. Editor: Beth Winkelstein.

J Biomech Eng 137(1), 011009 (Jan 01, 2015) (10 pages) Paper No: BIO-14-1360; doi: 10.1115/1.4029109 History: Received August 02, 2014; Revised November 10, 2014; Accepted November 17, 2014; Online December 10, 2014

Prior multiaxial strength studies on trabecular bone have either not addressed large variations in bone volume fraction and microarchitecture, or have not addressed the full range of multiaxial stress states. Addressing these limitations, we utilized micro-computed tomography (μCT) based nonlinear finite element analysis to investigate the complete 3D multiaxial failure behavior of ten specimens (5 mm cube) of human trabecular bone, taken from three anatomic sites and spanning a wide range of bone volume fraction (0.09–0.36), mechanical anisotropy (range of E3/E1 = 3.0–12.0), and microarchitecture. We found that most of the observed variation in multiaxial strength behavior could be accounted for by normalizing the multiaxial strength by specimen-specific values of uniaxial strength (tension, compression in the longitudinal and transverse directions). Scatter between specimens was reduced further when the normalized multiaxial strength was described in strain space. The resulting multiaxial failure envelope in this normalized-strain space had a rectangular boxlike shape for normal–normal loading and either a rhomboidal boxlike shape or a triangular shape for normal-shear loading, depending on the loading direction. The finite element data were well described by a single quartic yield criterion in the 6D normalized-strain space combined with a piecewise linear yield criterion in two planes for normal-shear loading (mean error ± SD: 4.6 ± 0.8% for the finite element data versus the criterion). This multiaxial yield criterion in normalized-strain space can be used to describe the complete 3D multiaxial failure behavior of human trabecular bone across a wide range of bone volume fraction, mechanical anisotropy, and microarchitecture.

FIGURES IN THIS ARTICLE
<>
Copyright © 2015 by ASME
Your Session has timed out. Please sign back in to continue.

References

Vasu, R., Carter, D. R., and Harris, W. H., 1981, “Stress Distributions in the Acetabular Region. I. Before and After Total Joint Replacement,” J. Biomech., 15(3), pp. 155–164. [CrossRef]
Carter, D. R., Orr, T. E., and Fyhrie, D. P., 1989, “Relationships Between Loading History and Femoral Cancellous Bone Architecture,” J. Biomech., 22(3), pp. 231–244. [CrossRef] [PubMed]
Cheal, E. J., Spector, M., and Hayes, W. C., 1992, “Role of Loads and Prosthesis Material Properties on the Mechanics of the Proximal Femur After Total Hip Arthroplasty,” J. Orthop. Res., 10(3), pp. 405–422. [CrossRef] [PubMed]
Stone, J. L., Beaupre, G. S., and Hayes, W. C., 1983, “Multiaxial Strength Characteristics of Trabecular Bone,” J. Biomech., 16(9), pp. 743–752. [CrossRef] [PubMed]
Keaveny, T. M., Wachtel, E. F., Zadesky, S. P., and Arramon, Y. P., 1999, “Application of the Tsai-Wu Quadratic Multiaxial Failure Criterion to Bovine Trabecular Bone,” ASME J. Biomech. Eng., 121(1), pp. 99–107. [CrossRef]
Fenech, C. M., and Keaveny, T. M., 1999, “A Cellular Solid Criterion for Predicting the Axial-Shear Failure Properties of Trabecular Bone,” ASME J. Biomech. Eng., 121(4), pp. 414–422. [CrossRef]
Rincon-Kohli, L., and Zysset, P. K., 2009, “Multi-Axial Mechanical Properties of Human Trabecular Bone,” Biomech. Model. Mechanobiol., 8(3), pp. 195–208. [CrossRef] [PubMed]
Niebur, G. L., Feldstein, M. J., and Keaveny, T. M., 2002, “Biaxial Failure Behavior of Bovine Tibial Trabecular Bone,” ASME J. Biomech. Eng., 124(6), pp. 699–705. [CrossRef]
Bayraktar, H. H., Gupta, A., Kwon, R. Y., Papadopoulos, P., and Keaveny, T. M., 2004, “The Modified Super-Ellipsoid Yield Criterion for Human Trabecular Bone,” ASME J. Biomech. Eng., 126(6), pp. 677–684. [CrossRef]
Wolfram, U., Gross, T., Pahr, D. H., Schwiedrzik, J., Wilke, H. J., and Zysset, P. K., 2012, “Fabric-Based Tsai-Wu Yield Criteria for Vertebral Trabecular Bone in Stress and Strain Space,” J. Mech. Behav. Biomed. Mater., 15, pp. 218–228. [CrossRef] [PubMed]
Cowin, S. C., 1986, “Fabric Dependence of an Anisotropic Strength Criterion,” Mech. Mater., 5(3), pp. 251–260. [CrossRef]
Zysset, P. K., and Rincon-Kohli, L., 2006, “An Alternative Fabric-Based Yield and Failure Criterion for Trabecular Bone,” Mechanics of Biological Tissue, G. A.Holzapfel and R. W.Ogden, eds., Springer, Berlin, Germany, pp. 457–470.
Garcia, D., Zysset, P. K., Charlebois, M., and Curnier, A., 2009, “A Three-Dimensional Elastic Plastic Damage Constitutive Law for Bone Tissue,” Biomech. Model. Mechanobiol., 8(2), pp. 149–165. [CrossRef] [PubMed]
Dall'Ara, E., Pahr, D., Varga, P., Kainberger, F., and Zysset, P., 2012, “QCT-Based Finite Element Models Predict Human Vertebral Strength In Vitro Significantly Better Than Simulated DXA,” Osteoporosis Int., 23(2), pp. 563–572. [CrossRef]
Chevalier, Y., Pahr, D., and Zysset, P. K., 2009, “The Role of Cortical Shell and Trabecular Fabric in Finite Element Analysis of the Human Vertebral Body,” ASME J. Biomech. Eng., 131(11), p. 111003. [CrossRef]
Sanyal, A., and Keaveny, T. M., 2013, “Biaxial Normal Strength Behavior in the Axial-Transverse Plane for Human Trabecular Bone–Effects of Bone Volume Fraction, Microarchitecture, and Anisotropy,” ASME J. Biomech. Eng., 135(12), p. 121010. [CrossRef]
Deshpande, V. S., and Fleck, N. A., 2001, “Multi-Axial Yield Behavior of Polymer Foams,” Acta Mater., 49(10), pp. 1859–1866. [CrossRef]
Wang, D. A., and Pan, J., 2006, “A Non-Quadratic Yield Function for Polymeric Foams,” Int. J. Plast., 22(3), pp. 434–458. [CrossRef]
Kopperdahl, D. L., and Keaveny, T. M., 1998, “Yield Strain Behavior of Trabecular Bone,” J. Biomech., 31(7), pp. 601–608. [CrossRef] [PubMed]
Morgan, E. F., and Keaveny, T. M., 2001, “Dependence of Yield Strain of Human Trabecular Bone on Anatomic Site,” J. Biomech., 34(5), pp. 569–577. [CrossRef] [PubMed]
Niebur, G. L., Yuen, J. C., Hsia, A. C., and Keaveny, T. M., 1999, “Convergence Behavior of High-Resolution Finite Element Models of Trabecular Bone,” ASME J. Biomech. Eng., 121(6), pp. 629–635. [CrossRef]
Charras, G. T., and Guldberg, R. E., 2000, “Improving the Local Solution Accuracy of Large-Scale Digital Image-Based Finite Element Analyses,” J. Biomech., 33(2), pp. 255–259. [CrossRef] [PubMed]
van Rietbergen, B., Odgaard, A., Kabel, J., and Huiskes, R., 1996, “Direct Mechanics Assessment of Elastic Symmetries and Properties of Trabecular Bone Architecture,” J. Biomech., 29(12), pp. 1653–1657. [CrossRef] [PubMed]
Turner, C. H., and Cowin, S. C., 1988, “Errors Introduced by Off-Axis Measurements of the Elastic Properties of Bone,” J. Biomech., 110(3), pp. 213–214. [CrossRef]
Papadopoulos, P., and Lu, J., 2001, “On the Formulation and Numerical Solution of Problems in Anisotropic Finite Plasticity,” Comput. Methods Appl. Mech. Eng, 190(37–38), pp. 4889–4910. [CrossRef]
Bevill, G., Eswaran, S. K., Gupta, A., Papadopoulos, P., and Keaveny, T. M., 2006, “Influence of Bone Volume Fraction and Architecture on Computed Large-Deformation Failure Mechanisms in Human Trabecular Bone,” Bone, 39(6), pp. 1218–1225. [CrossRef] [PubMed]
Sanyal, A., Gupta, A., Bayraktar, H. H., Kwon, R. Y., and Keaveny, T. M., 2012, “Shear Strength Behavior of Human Trabecular Bone,” J. Biomech., 45(15), pp. 2513–2519. [CrossRef] [PubMed]
Fields, A. J., Nawathe, S., Eswaran, S. K., Jekir, M. G., Adams, M. F., Papadopoulos, P., and Keaveny, T. M., 2012, “Vertebral Fragility and Structural Redundancy,” J. Bone Miner. Res., 27(10), pp. 2152–2158. [CrossRef] [PubMed]
Nawathe, S., Akhlaghpour, H., Bouxsein, M. L., and Keaveny, T. M., 2013, “Microstructural Failure Mechanisms in the Human Proximal Femur for Sideways Fall Loading,” J. Bone Miner. Res., 29(2), pp. 507–515. [CrossRef]
Tsai, S., and Wu, E., 1971, “A General Theory for Strength of Anisotropic Materials,” J. Compos. Mater., 5(1), pp. 58–80. [CrossRef]
Chakrabarty, J., 1987, Theory of Plasticity, McGraw-Hill, Inc. New York.
Gibson, L. J., Ashby, M. F., Zhang, J., and Triantafillou, T. C., 1989, “Failure Surfaces for Cellular Materials Under Multiaxial Loads—I. Modelling,” Int. J. Mech. Sci., 31(9), pp. 635–663. [CrossRef]
Walker, J. D., and Thacker, B. H., 1999, Yield Surfaces for Anisotropic Plates, American Institute of Physics, Snowbird, UT, pp. 567–570.
Barlat, F., Lege, D. J., and Brem, J. C., 1991, “A Six-Component Yield Function for Anisotropic Materials,” Int. J. Plast., 7(7), pp. 693–712. [CrossRef]
Karafillis, A. P., and Boyce, M. C., 1993, “A General Anisotropic Yield Criterion Using Bounds and a Transformation Weighting Tensor,” J. Mech. Phys. Solids, 41(12), pp. 1859–1886. [CrossRef]
Gol'denblat, I. I., and Kopnov, V. A., 1965, “Strength of Glass-Reinforced Plastics in the Complex Stress State,” Polym. Mech., 1(2), pp. 54–59. [CrossRef]
von Mises, R., 1928, “Mechanik der Plastischen Formäenderung von Kristallen,” Z. Angew. Math. Mech., 8(3), pp. 161–185. [CrossRef]
Hill, R., 1983, The Mathematical Theory of Plasticity, Clarendon, Oxford University, Oxford, New York.
Hoffman, O., 1967, “The Brittle Strength of Orthotropic Materials,” J. Compos. Mater., 1(2), pp. 200–206. [CrossRef]
Schwiedrzik, J. J., Wolfram, U., and Zysset, P. K., 2013, “A Generalized Anisotropic Quadric Yield Criterion and its Application to Bone Tissue at Multiple Length Scales,” Biomech. Model. Mechanobiol., 12(6), pp. 1155–1168. [CrossRef] [PubMed]
Bower, M. V., and Koedam, D. H., 1997, “Tensor Polynomial Failure Criterion: Coefficient Limits Based on Convexity Requirements,” J. Reinf. Plast. Compos., 16(5), pp. 435–477.
Cowin, S. C., and He, Q.-C., 2005, “Tensile and Compressive Stress Yield Criteria for Cancellous Bone,” J. Biomech., 38(1), pp. 141–144. [CrossRef] [PubMed]
Arramon, Y. P., Mehrabadi, M. M., Martin, D. W., and Cowin, S. C., 2000, “A Multidimensional Anisotropic Strength Criterion Based on Kelvin Modes,” Int. J. Solids Struct., 37(21), pp. 2915–2935. [CrossRef]
Kelly, N., and McGarry, J. P., 2012, “Experimental and Numerical Characterisation of the Elasto-Plastic Properties of Bovine Trabecular Bone and Trabecular Bone Analogue,” J. Mech. Behav. Biomed. Mater., 9, pp. 184–197. [CrossRef] [PubMed]
Troy, K. L., and Grabiner, M. D., 2007, “Off-Axis Loads Cause Failure of the Distal Radius at Lower Magnitudes Than Axial Loads: A Finite Element Analysis,” J. Biomech., 40(8), pp. 1670–1675. [CrossRef] [PubMed]
Bevill, G., Farhamand, F., and Keaveny, T. M., 2009, “Heterogeneity of Yield Strain in Low-Density Versus High-Density Human Trabecular Bone,” J. Biomech., 42(13), pp. 2165–2170. [CrossRef] [PubMed]
Chang, W. C., Christensen, T. M., Pinilla, T. P., and Keaveny, T. M., 1999, “Uniaxial Yield Strains for Bovine Trabecular Bone are Isotropic and Asymmetric,” J. Orthop. Res., 17(4), pp. 582–585. [CrossRef] [PubMed]
Nawathe, S., Juillard, F., and Keaveny, T. M., 2013, “Theoretical Bounds for the Influence of Tissue-Level Ductility on the Apparent-Level Strength of Human Trabecular Bone,” J. Biomech., 46(7), pp. 1293–1299. [CrossRef] [PubMed]
Yeh, O. C., and Keaveny, T. M., 2001,“Relative Roles of Microdamage and Microfracture in the Mechanical Behavior of Trabecular Bone,” J. Orthop. Res., 19(6), pp. 1001–1007. [CrossRef] [PubMed]
Guo, X. E., McMahon, T. A., Keaveny, T. M., Hayes, W. C., and Gibson, L. J., 1994, “Finite Element Modeling of Damage Accumulation in Trabecular Bone Under Cyclic Loading,” J. Biomech., 27(2), pp. 145–155. [CrossRef] [PubMed]
Silva, M. J., and Gibson, L. J., 1997, “The Effects of Non-Periodic Microstructure and Defects on the Compressive Strength of Two-Dimensional Cellular Solids,” Intern. J. Mech. Sci., 39(5), pp. 549–563. [CrossRef]

Figures

Grahic Jump Location
Fig. 1

Yield strain (top row), normalized-yield-strain (middle row), and normalized-yield-stress (bottom row) in the three normal biaxial planes, showing individual responses for all ten specimens. The same tissue-level material properties were assumed for all specimens. For each biaxial plane, the virtually applied normal strain (top and middle row) and normal stress (bottom row) in the third normal direction is zero. Scatter was much less when the normalized-yield-strain was used to describe the multiaxial failure behavior.

Grahic Jump Location
Fig. 2

Normalized-yield-strain in the nine normal-shear biaxial planes, showing individual responses for all ten specimens. For each biaxial plane, the virtually applied normal and shear strains in the other two normal and shear directions, respectively, are zero.

Grahic Jump Location
Fig. 3

Mean absolute error (±SD, n = 10 specimens) of the specimen-specific mathematical fit versus the finite element computed normalized-yield-strain data for (a) the normal (only) 3D space using the quartic yield criterion, (b) the seven normal-shear planes (ɛ∧11-γ∧12, ɛ∧11-γ∧13, ɛ∧11-γ∧23, ɛ∧22-γ∧12, ɛ∧22-γ∧13, ɛ∧22-γ∧23, ɛ∧33-γ∧12) using the quartic yield criterion, (c) the two normal-shear planes (ɛ∧33-γ∧13, ɛ∧33-γ∧23) using the quartic yield criterion, and (d) the two normal-shear planes (ɛ∧33-γ∧13, ɛ∧33-γ∧23) planes using a combination of the quartic and piecewise linear yield criteria.

Grahic Jump Location
Fig. 4

Yield envelopes in the three biaxial normal normalized-strain planes and the nine normal-shear normalized-strain planes for one specimen from the vertebral body (BV/TV = 0.11) using the quartic yield criterion (solid line, Eq. (4)) and the piecewise linear-yield criterion (dashed line, Eq. (5)). For all plots, the solid circles fail in the mode denoted by the horizontal axis and the hollow circles fail in the mode denoted by the vertical axis. For the quartic criterion alone, the mean error of the fit was 4.2% for the combined ɛ∧33-ɛ∧11, ɛ∧33-ɛ∧22, and ɛ∧22-ɛ∧11 planes, 3.1% for the combined ɛ∧11-γ∧12, ɛ∧11-γ∧13, ɛ∧11-γ∧23, ɛ∧22-γ∧12, ɛ∧22-γ∧13, ɛ∧22-γ∧23, and ɛ∧33-γ∧12 planes, and 7.5% for the combined ɛ∧33-γ∧13 and ɛ∧33-γ∧23 planes. For the latter two planes, when the quartic criterion was combined with the piecewise linear criterion (QPL criterion), the mean error reduced to 2.9%. In this case, the assumed failure envelope is the inner surface of the two individual criteria.

Tables

Errata

Discussions

Some tools below are only available to our subscribers or users with an online account.

Related Content

Customize your page view by dragging and repositioning the boxes below.

Related Journal Articles
Related eBook Content
Topic Collections

Sorry! You do not have access to this content. For assistance or to subscribe, please contact us:

  • TELEPHONE: 1-800-843-2763 (Toll-free in the USA)
  • EMAIL: asmedigitalcollection@asme.org
Sign In