A Phenomenological Model for Predicting Fatigue Life in Bovine Trabecular Bone

[+] Author and Article Information
P. Ganguly, T. L. A. Moore, L. J. Gibson

Department of Materials Science and Engineering, Massachusetts Institute of Technology, Cambridge, MA 02139Exponent Failure Analysis Associates, Inc. Philadelphia, PA 19104

J Biomech Eng 126(3), 330-339 (Jun 24, 2004) (10 pages) doi:10.1115/1.1762893 History: Received August 01, 2003; Revised October 27, 2003; Online June 24, 2004
Copyright © 2004 by ASME
Your Session has timed out. Please sign back in to continue.


Freeman,  M. A. R., Todd,  R. C., and Pirie,  C. J., 1974, “Role of Fatigue in Pathogenesis of Senile Femoral Neck Fractures,” J. Bone Jt. Surg., 56B, pp. 698–702.
Riggs,  B. L., and Melton,  L. J., 1995, “The Worldwide Problem of Osteoporosis: Insights Offered by Epidemiology,” Bone (N.Y.), 17, pp. S505–S511.
Daffner,  R. H., and Pavlov,  H., 1992, “Stress Fractures: Current Concepts,” Am. J. Roentgenol., 159, pp. 245–252.
Burr,  D. B., Forwood,  M. R., Fyhrie,  D. P., Martin,  R. B., Schaffler,  M. B., and Turner,  C. H., 1997, “Bone Microdamage and Skeletal Fragility in Osteoporotic and Stress Fractures,” J. Bone Miner. Res., 12, pp. 6–15.
Giladi,  M., Milgrom,  C., Kashtan,  H., Stein,  M., Chisin,  R., and Dizian,  R., 1986, “Recurrent Stress Fractures in Military Recruits. One Year Follow-up of 66 Recruits,” J. Bone Jt. Surg., 68B, pp. 439–441.
Matheson,  G. O., Clement,  D. B., McKenzie,  D. C., Taunton,  J. E., Llyod-Smith,  D. R., and MacIntyre,  J. G., 1987, “Stress Fractures in Athletes: A Study of 320 Cases,” Am. J. Sports Med., 15, pp. 46–58.
Lane,  J. M., Russell,  L., and Khan,  S. N., 2000, “Osteoporosis,” Clin. Orthop., 372, pp. 139–150.
Silva,  M. J., Keaveny,  T. M., and Hayes,  W. C., 1997, “Load Sharing Between the Shell and Centrum in the Lumbar Vertebral Body,” Spine, 22, pp. 140–150.
Michel,  M. C., Guo,  X.-D. E., Gibson,  L. J., McMahon,  T. A., and Hayes,  W. C., 1993, “Compressive Fatigue Behavior of Bovine Trabecular Bone,” J. Biomech., 26, pp. 453–463.
Bowman,  S. M., Guo,  X. E., Cheng,  D. W., Keaveny,  T. M., Gibson,  L. J., Hayes,  W. C., and McMahon,  T. A., 1998, “Creep Contributes to the Fatigue Behavior of Bovine Trabecular Bone,” J. Biomech. Eng., 120, pp. 647–654.
Moore,  T. L. A., and Gibson,  L. J., 2003, “Fatigue of Bovine Trabecular Bone,” J. Biomech. Eng. 125, pp. 761–768.
Haddock, S. M., Yeh, O. C., Mummaneni, P. V., Rosenberg, W. S., and Keaveny, T. M., 2000, “Fatigue Behavior of Human Vertebral Trabecular Bone,” 46th Annual Meeting, Orthopaedic Research Society, pg. 733.
Moore,  T. L. A., and Gibson,  L. J., 2003, “Fatigue Microdamage in Bovine Trabecular Bone,” ASME J. Biomech. Eng., 125, pp. 769–776.
Moore,  T. L. A., and Gibson,  L. J., 2002, “Microdamage Accumulation in Bovine Trabecular Bone in Uniaxial Compression,” ASME J. Biomech. Eng., 124, pp. 63–71.
Schaffler,  M. B., Radin,  E. L., and Burr,  D. B., 1989, “Mechanical and Morphological Effects of Strain Rate on Fatigue of Compact Bone,” Bone (N.Y.), 10, 207–214.
Burr,  D. B., Turner,  C. H., Naick,  P., Forwood,  M. R., Ambrosius,  W., Hasan,  M. S., and Pidaprati,  R., 1998, “Does Microdamage Accumulation Affect the Mechanical Properties of Bone?,” J. Biomech., 31, pp. 337–345.
Taylor,  D., O’Brien,  F., and Lee,  T. C., 2002, “A Theoretical Model for the Simulation of Microdamage Accumulation and Repair in Compact Bone,” Meccanica, 37, pp. 397–406.
Moore,  T. L. A., O’Brien,  F. J., and Gibson,  L. J., 2004, “Creep Does Not Contribute to Fatigue in Bovine Trabecular Bone,” ASME J. Biomech. Eng., 126, pp. 321–329.
Guo,  X.-D. 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, pp. 145–155.
Schaffner,  G., Guo,  X. E., Silva,  M. J., and Gibson,  L. J., 2000, “Modelling Fatigue Damage Accumulation in Two-dimensional Voronoi Honeycombs,” Int. J. Mech. Sci., 42, pp. 645–656.
Makiyama,  A. M., Vajjhala,  S., and Gibson,  L. J., 2002, “Analysis of Crack Growth in a 3D Voronoi Structure: A Model for Fatigue in Low Density Trabecular Bone,” J. Biomech. Eng., 124, pp. 512–520.
Paris,  P. C., and Erdogan,  F., 1963, “A Critical Analysis of Crack Propagation Laws,” ASME J. Basic Eng., 85, pp. 528–534.
Carter,  D. R., and Caler,  W. E., 1983, “Cycle-dependent and Time-dependent Bone Fracture With Repeated Loading,” J. Biomech. Eng., 105, pp. 166–170.
Carter,  D. R., and Caler,  W. E., 1985, “A Cumulative Damage Model for Bone Fracture,” J. Orthop. Res., 3, pp. 84–90.
Taylor,  D., 2002, “Modelling of Fatigue Crack Growth at the Microstructural Level,” Comput. Mater. Sci., 25, pp. 228–236.
Keaveny,  T. M., Guo,  X.-D. E., Wachtel,  E. F., McMahon,  T. A., and Hayes,  W. C., 1994, “Trabecular Bone Exhibits Fully Linear Elastic Behavior and Yields at Low Strains,” J. Biomech., 27, pp. 1127–1136.
Keaveny,  T. M., Wachtel,  E. F., Guo,  X.-D. E., and Hayes,  W. C., 1994, “Mechanical Behavior of Damaged Trabecular Bone,” J. Biomech., 27, pp. 1309–1318.
Keaveny,  T. M., Pinilla,  T. P., Crawford,  R. P., Kopperdahl,  D. L., and Lou,  A., 1997, “Systematic and Random Errors in Compression Testing of Trabecular Bone,” J. Orthop. Res., 15, pp. 101–110.
Gibson,  L. J., 1985, “The Mechanical Behavior of Cancellous Bone,” J. Biomech., 18, pp. 317–328.
Keaveny,  T. M., Morgan,  E. F., Niebur,  G. L., and Yeh,  O. C., 2001, “Biomechanics of Trabecular Bone,” Annu. Rev. Biomed. Eng., 3, pp. 307–333.
Morgan,  E. F., Bayraktar,  H. H., and Keaveny,  T. M., 2003, “Trabecular Bone Modulus-Density Relationships Depend on Anatomic Site,” J. Biomech., 36, pp. 897–904.
Keaveny,  T. M., Wachtel,  E. F., and Kopperdahl,  D. L., 1999, “Mechanical Behavior of Human Trabecular Bone after Overloading,” J. Orthop. Res., 17, pp. 346–353.
Suresh, S., 1991, Fatigue of materials, Cambridge University Press, Cambridge, U.K., p. 127.
Lemaitre, J., 1992, A Course on Damage Mechanics, Springer-Verlag, p. 110.


Grahic Jump Location
A characteristic stress-strain curve for bovine trabecular bone under monotonic compression. The initial modulus (E0m) was calculated by the best-fit line in the strain range 0.001–0.004. The secant modulus (Esecm) at higher strains were calculated by the ratio of the stress and strain.
Grahic Jump Location
Calculation of the secant modulus for a given cycle for a sample tested at a maximum compressive normalized stress (σ/E0f) of 0.0066.
Grahic Jump Location
Variation of the normalized secant modulus with increasing strain for representative specimens under monotonic compression. The secant modulus variation for each sample has been represented by the best-fit linear plot.
Grahic Jump Location
Variation of secant moduli with increasing strain for the trabecular bone samples at normalized stress (Δσ/E0f) ranges of (a) (0.005, 0.006], (b) (0.006, 0.007] (c) (0.007, 0.008] and (d) (0.008, 0.009]. The experimental points have been identified by the sample numbers. The upper and lower bounds from the monotonic tests have been included for comparison, and shown in the plots by the solid lines (upper and lower).
Grahic Jump Location
Plots showing the accumulation of total residual strain with increasing number of cycles, for typical samples in the (0.005, 0.006], (0.006, 0.007], (0.007, 0.008] and (0.008, 0.009] normalized stress ranges. (a) The initial part of the residual strain plots, showing the primary and the secondary phases. (b) The secondary residual strain evolution at large number of cycles. The plot for Δσ/E0=(0.008,0.009] has been omitted from (b), the fatigue life at that stress level being comparable to the number of cycles shown in (a). Note that the residual strain analysis for the same normalized stress range shown in (a) and (b) corresponds to different specimens.
Grahic Jump Location
Experimentally observed variation in the secondary residual strain accumulation rate with changes in normalized stress. The best-fit line was obtained by fitting a power-law curve on the entire experimental data (open, solid and gray squares). The points that appear to be bounding the data were identified (upper bound: solid squares, lower bound: gray squares) and separate power-law curves were fit to obtain the upper and lower bound variation equations.
Grahic Jump Location
Upper and lower bounds of the decrease of the normalized secant modulus predicted by the numerical model, for normalized stress ranges (a) (0.005, 0.006], (b) (0.006, 0.007], (c) (0.007, 0.008] and (d) (0.008, 0.009]. Note that upper and lower bound responses for both f1 and f2 have been considered in determining the model responses. The experimental points, identified by their sample numbers, have been included for comparison.
Grahic Jump Location
Numerical model prediction of the upper and lower bounds of the normalized secant modulus reduction, at normalized stress ranges (a) (0.005, 0.006] and (b) (0.008, 0.009]. Note that for the model bound constructions, upper and lower bound responses for f1 and a median response for f2 (see Table 1) was assumed. The experimental points, identified by their sample numbers, have been included for comparison.
Grahic Jump Location
The predicted S-N curve for bovine trabecular bone. The dotted line signifies the number of cycles corresponding to the endurance limit. The lower and upper bounds for the endurance limit are shown in the figure. The experimental data reflect cyclic compressive loadings of waisted samples of bovine trabecular bone 10.



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