0
TECHNICAL PAPERS

Analysis of Crack Growth in a 3D Voronoi Structure: A Model for Fatigue in Low Density Trabecular Bone

[+] Author and Article Information
A. M. Makiyama, S. Vajjhala, L. J. Gibson

Department of Materials Science and Engineering, Massachusetts Institute of Technology, Cambridge, MA 02139

J Biomech Eng 124(5), 512-520 (Sep 30, 2002) (9 pages) doi:10.1115/1.1503792 History: Received December 01, 2001; Revised May 01, 2002; Online September 30, 2002
Copyright © 2002 by ASME
Your Session has timed out. Please sign back in to continue.

References

Riggs,  B. L., and Melton,  L. J., 1995, “The Worldwide Problem of Osteoporosis: Insights Afforded by Epidemiology,” Bone (N.Y.), 17, pp. 505S–511S.
Silva,  M. J., Wang,  C., Keaveny,  T. M., and Hayes,  W. C., 1994, “Direct and Computed-Tomography Thickness Measurements of the Human Lumbar Vertebral Shell and End-Plate,” Bone (N.Y.), 15, pp. 409–414.
McBroom,  R. J., Hayes,  W. C., Edwards,  W. T., Goldberg,  R. P., White,  A. A., 1985, “Prediction of Vertebral Body Compressive Fracture Using Quantitative Computed Tomography,” J. Bone Jt. Surg., 67A, pp. 1206–1214.
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.
Freeman,  M. A., Todd,  R. C., and Pirie,  C. J., 1974, “The Role of Fatigue in the Pathogenesis of Senile Femoral Neck Fractures,” J. Bone Jt. Surg., 56B, pp. 698–702.
Burrows,  H. J., 1956, “Fatigue Infraction of the Middle of Tibia in Ballet Dancers,” J. Bone Jt. Surg., 38B, pp. 83–94.
Giladi,  M., Milgrom,  C., Kashtan,  H., Stein,  M., Chisin,  R., and Dizian,  R., 1986, “Recurrent Stress Fracture 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., Lloyd-Smith,  D. R., and MacIntyre,  J. G., 1987, “Stress Fracture in Athletes: A Study of 320 Cases,” American Journal of Sports Medicine , 15, pp. 46–58.
Guo,  XDE, 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.
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.
Arthur Moore, T. L., 2001, “Microdamage Accumulation in Bovine Trabecular Bone,” PhD Thesis, Health Sciences and Technology, MIT.
Arthur Moore T. L., and Gibson L. J., 2001, “Microdamage Accumulation During Compressive Fatigue of Bovine Trabecular Bone,” ASME Bioengineering Conference, Snowbird, Utah, June 2001, pp. 291–292.
Chang,  W. C. W., 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, pp. 582–585.
Morgan,  E. F., and Keaveny,  T. M., 2001, “Dependence of Yield Strain of Human Trabecular Bone on Anatomic Site,” J. Biomech., 34, pp. 569–577.
Caler,  W. E., and Carter,  D. R., 1989, “Bone Creep-Fatigue Damage Accumulation,” J. Biomech., 22, pp. 625–635.
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.
Wright,  T. M., and Hayes,  W. C., 1976, “The Fracture Mechanics of Fatigue Crack Propagation in Compact Bone,” J. Biomed. Mater. Res., 7, pp. 637–648.
Fondrk,  M., Bahniuk,  E., Davy,  D. T., and Michaels,  C., 1988, “Some Viscoplastic Characteristics of Bovine and Human Cortical Bone,” J. Biomech., 21, pp. 623–630.
Huang, J.-S., and Gibson, L. J., 2001, “Creep of Open-Cell Voronoi Foams,” Mater. Sci. Eng. A, in press.
Voronoi,  G. F., 1908, “Novelles applications des parameters continus a la theorie des formes quadratique. Deuxieme Memoire: recherches sur les parallelloedres primitifs,” J. Reine Angew. Math., 134, pp. 198–287.
Silva,  M. J., Hayes,  W. C., and Gibson,  L. J., 1995, “The Effects of Non-Periodic Microstructure on the Elastic Properties of Two-Dimensional Cellular Solids,” Int. J. Mech. Sci., 37, pp. 1161–1177.
Sullivan, J. M., 1988, “VCS,” software available from University of Minnesota Geometry Center, University of Minnesota, Minneapolis, MN.
Sullivan, J. M., 1988, “Torus munge,” software available from University of Minnesota Geometry Center, University of Minnesota, Minneapolis, MN.
Sullivan, J. M., 1997, “Cells,” software available from University of Minnesota Geometry Center, University of Minnesota, Minneapolis, MN.
Brakke,  K. A., 1992, “The Surface Evolver,” Exp. Math., 1, pp. 141–165.
Kraynik, A. M., Neisen, M. K., Reinelt, D. A., and Warren, W. E., 1999, “Foam Micromechanics,” in Foams and Emulsions, Sadoc, J. F., and Rivier, N., eds., NATO ASI Series, Kluwer Academic, pp. 259–286.
Vajjhala,  S., Kraynik,  A. M., and Gibson,  L. J., 2000, “A Cellular Solid Model for Modulus Reduction Due to Resorption of Trabeculae in Bone,” J. Biomech. Eng., 122, pp. 511–515.
Morita,  M., Ebihara,  A., Itoman,  M., and Sasada,  T., 1994, “Progression of Osteoporosis in Cancellous Bone Depending on Trabecular Structure,” Ann. Biomed. Eng., 22, pp. 532–539.
Muller,  R., Van Campenhout,  H., Van Damme,  B., Van Der Perre,  G., Dequeker,  J., Hildebrand,  T., and Ruegsegger,  P., 1998, “Morphometric Analysis of Human Bone Biopsies: A Quantitative Structural Comparison of Histological Sections and Micro-Computed Tomography,” Bone (N.Y.), 23, pp. 59–66.
Rho,  J. Y., Roy,  M. E., Tsui,  T. Y., and Pharr,  G. M., 1999, “Elastic Properties of Microstructural Components of Human Bone Tissue Measured by Nanoindentation,” J. Biomed. Mater. Res., 45, pp. 48–54.
Fazzalari,  N. L., Forwood,  M. R., Smith,  K., Manthey,  B. A., and Herreen,  P., 1998, “Assessment of Cancellous Bone Quality in Severe Osteoarthrosis: Bone Mineral Density, Mechanics and Microdamage,” Bone (N.Y.), 22, pp. 381–388.
Arthur Moore,  T. L., and Gibson,  L. J., 2002, “Microdamage Accumulation in Bovine Trabecular Bone in Uniaxial Compression,” J. Biomech. Eng., 124, pp. 63–71.
Yeni,  Y. N., Brown,  C. U., and Norman,  T. L., 1998, “Influence of Bone Composition and Apparent Density on Fracture Toughness of the Human Femur and Tibia,” Bone (N.Y.), 22, pp. 79–84.
Norman,  T. L., Yeni,  Y. N., Brown,  C. U., and Wang,  Z., 1998, “Influence of Microdamage on Fracture Toughness of the Human Femur and Tibia,” Bone (N.Y.), 23, pp. 303–306.
Levan,  A., and Royer,  J., 1993, “Part-Circular Surface Cracks in Round Bars under Tension, Bending and Twisting,” Int. J. Fract., 61, pp. 71–99.
Taylor,  D., 1998, “Microcrack Growth Parameters for Compact Bone Deduced from Stiffness Variations,” J. Biomech., 31, pp. 587–592.
Taylor, D., and Prendergast, P., 1997, “A Model for Fatigue Crack Propagation and Remodeling in Compact Bone,” Proceedings of the Institution of Mechanical Engineers, Part H 211, pp. 369–375.
Akkus,  O., and Rimnac,  C. M., 2001, “Cortical Bone Tissue Resists Fatigue Fracture by Deceleration of Microcrack Growth,” J. Biomech., 34, pp. 757–764.
O’Brien,  F. J., Taylor,  D., Dickson,  G. R., and Lee,  T. C., 2000, “Visualization of Three-Dimensional Microcracks in Compact Bone,” J. Anat., 197, pp. 413–420.

Figures

Grahic Jump Location
(a) Construction of a two-dimensional Voronoi structure. (b) A Voronoi honeycomb with about 500 cells (Reproduced from Silva et al. (1995), with the kind permission of Elsevier Science.)
Grahic Jump Location
(a) A three-dimensional open-cell Voronoi structure with 27 cells and roughly 300 struts. In our models, the strut diameter was 150μm and the average length of the struts was 588 μm. (b) A micro-computed tomography image of trabecular bone with a relative density of 9.4%. (Reproduced with the kind permission of Dr. Ralph Muller, ETH, Zurich.)
Grahic Jump Location
Flow chart giving an overview of the fatigue analysis.
Grahic Jump Location
The three crack shapes used in our analysis. The crack shape is characterized by the parameter α.
Grahic Jump Location
Normalized stress intensity factors, ΔK/σπa plotted against relative crack length, a/R, for different crack shapes, for (a) tension (b) bending and (c) torsion. (Reproduced from Levan and Royer, 1993 (Reproduced from Int. J. Fracture with the kind permission of Kluwer Academic Publishers.)
Grahic Jump Location
(a) Young’s modulus, normalized by the initial modulus, E/Eo, plotted against number of cycles of fatigue loading, N on a linear scale. The data shown are for Δσ/Eo=0.006 and α=0.5, for structure A. (b) Young’s modulus, normalized by the initial modulus, E/Eo, plotted against number of cycles of fatigue loading, N on a semi-log scale. Open symbols represent experimental results 11, closed symbols represent model results.
Grahic Jump Location
Average of the normalized Young’s modulus, E/Eo, for the three Voronoi structures plotted against number of cycles of loading, N. The error bars represent one standard deviation. All results for α=0.5.
Grahic Jump Location
Normalized stress range, Δσ/Eo, plotted against the number of cycles to failure, Nf. Failure is defined as the number of cycles to give E/Eo=0.90. (a) For relative densities, ρ*s=1, 5 and 10%, for ao=33 μ,α=0 for Voronoi structure A. (b) For initial crack length, ao=16, 33, 66 μ, for ρ*s=5% and α=0 for Voronoi structure A. (c) For crack shapes, α=0, 0.5, 1.0, for ρ*s=5% and ao=33 μ for Voronoi structure A. (d) For three different Voronoi structures, for ρ*s=5%,ao=33 μ and α=0.
Grahic Jump Location
Stress range plotted against number of cycles to failure, for ρ*s=1, 5 and 10%, for ao=33 μ and α=0 for Voronoi structure A
Grahic Jump Location
(a) Number of cycles to failure, Nf plotted against relative density, ρ*s for Δσ/Eo=0.4%, 0.6%, 1.0%, for ao=33 μ,α=0 for Voronoi structure A. (b) Number of cycles to failure, Nf plotted against normalized stress range, Δσ/Eo, for α=0, 0.5, 1.0, for ao=16 μ, and ρ*s=5% for Voronoi structure A. (c) Number of cycles to failure, Nf, plotted against α for Δσ/Eo=0.4%, 0.6%, 1.0%, for ρ*s=5%,ao=33 μ, for Voronoi structure A. (d) Number of cycles to failure, Nf, plotted against crack size, ao, for Δσ/Eo=0.4%, 0.6%, 1.0%, for ρ*s=5%, and α=0 for Voronoi structure A.
Grahic Jump Location
Normalized stress range, Δσ/Eo plotted against number of cycles to failure, Nf comparing numerical model for open cell Voronoi structure and experimental data for bovine trabecular bone 1011. Finite element model with ao=33 μ and α=0 for Voronoi structure A.

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