Research Papers

Mechanical Evaluation of Hydroxyapatite Nanocomposites Using Finite Element Modeling

[+] Author and Article Information
Ani Ural

e-mail: ani.ural@villanova.edu
Department of Mechanical Engineering
Villanova University
800 Lancaster Avenue
Villanova, PA 19085

A Please provide complete page range for Ref. [32].

1Corresponding author.

Contributed by the Materials Division of ASME for publication in the JOURNAL OF ENGINEERING MATERIALS AND TECHNOLOGY. Manuscript received January 13, 2012; final manuscript received October 24, 2012; published online January 23, 2013. Editor: Hussein Zbib.

J. Eng. Mater. Technol 135(1), 011007 (Jan 23, 2013) (9 pages) Paper No: MATS-12-1010; doi: 10.1115/1.4023187 History: Received January 13, 2012; Revised October 24, 2012

Hydroxyapatite (HA) has been proposed as a candidate material for bone implants because of its similarity to the inorganic phase in bone. However, due to its lower mechanical properties compared to bone, it has not been used in load bearing bone implants. Inclusion of second phase reinforcements in HA such as carbon nanotubes (CNT) and graphene nanosheets is expected to significantly improve its mechanical properties. In this study, a computational framework that will improve the understanding of the mechanical behavior of graphene nanosheet and CNT-reinforced HA-nanocomposites is proposed. The variation of elastic modulus of HA-nanocomposites is assessed based on the nanofiller type, volume fraction, alignment, area, thickness, and aspect ratio using the finite element modeling. The results of the simulations show that graphene nanosheets are more effective in improving the elastic modulus of nanocomposites than CNTs at similar volume fractions. HA-nanocomposites reinforced by graphene nanosheets exhibit transversely isotropic material properties and provide the highest elastic modulus when aligned along a direction or randomly distributed in a plane, whereas CNTs provide the best reinforcement when aligned along an axis. Variation in graphene nanosheet area, thickness, aspect ratio, and carbon nanotube length have negligible effect on elastic modulus of the HA-nanocomposite. In addition, comparison between the finite element simulations and theoretical calculations show that clustering of nanoinclusions reduces the effectiveness of the reinforcement they provide. The simulation results and the computational framework presented in this study are expected to help in determining the best design and manufacturing parameters that can be adapted for developing HA-nanocomposite bone implant materials.

© 2013 by ASME
Your Session has timed out. Please sign back in to continue.


De With, G., Vandijk, H. J. A., Hattu, N., and Prijs, K., 1981, “Preparation, Microstructure and Mechanical Properties of Polycrystalline Hydroxyapatite,” J. Mater. Sci., 16, pp. 1592–1598. [CrossRef]
Hench, L. L., 1991, “Bioceramics: From Concept to Clinic,” J. Am. Ceram. Soc., 74(7), pp. 1487–1510. [CrossRef]
White, A. A., Best, S. M., and Kinloch, I. A., 2007, “Hydroxyapatite-Carbon Nanotube Composites for Biomedical Applications: A Review,” Int. J. Appl. Ceram. Technol., 4(1), pp. 1–13. [CrossRef]
Asmus, S. M., Sakakura, S., and PezzottiG., 2003, “Hydroxyapatite Toughened by Silver Inclusions,” J. Compos. Mater., 37(23), pp. 2117–2129. [CrossRef]
Chen, Y., Zhang, Y. Q., Zhang, T. H., Gan, C. H., Zheng, C. Y., and Yu, G., 2006, “Carbon Nanotube Reinforced Hydroxyapatite Composite Coatings Produced Through Laser Surface Alloying,” Carbon, 44, pp. 37–45. [CrossRef]
Guo, H., Khor, K. A., Boey, Y. C., and Miao, X., 2003, “Laminated and Functionally Graded Hydroxyapatite/Yttria Stabilized Tetragonal Zirconia Composites Fabricated by Spark Plasma Sintering,” Biomaterials, 24(4), pp. 667–675. [CrossRef] [PubMed]
Kobayashi, S., and Kawai, W., 2007, “Development of Carbon Nanofiber Reinforced Hydroxyapatite With Enhanced Mechanical Properties,” Composites: Part A, 38, pp. 114–123. [CrossRef]
Lahiri, D., Singh, V., Keshri, A., Seal, S., and Agarwal, A., 2010, “Carbon Nanotube Toughened Hydroxyapatite by Spark Plasma Sintering: Microstructural Evolution and Multi-Scale Tribological Properties,” Carbon, 48(11), pp. 3103–3120. [CrossRef]
Li, H., Zhao, N., Liu, Y., Liang, C., Shi, C., Du, X., and LiJ., 2008, “Fabrication and Properties of Carbon Nanotubes Reinforced Fe/Hydroxyapatite Composites by In Situ Chemical Vapor Deposition,” Composites, Part A, 39, pp. 1128–1132. [CrossRef]
Meng, Y. H., Tang, C. Y., Tsui, C. P., and Chen da, Z., 2008, “Fabrication and Characterization of Needle-Like Nano-HA and HA/MWNT Composites,” J. Mater. Sci.: Mater. Med., 19(1), pp. 75–81. [CrossRef] [PubMed]
Pezzotti, G., Asmus, S. M., Ferroni, L. P., and Miki, S., 2002, “In Situ Polymerization Into Porous Ceramics: A Novel Route to Tough Biomimetic Materials,” J. Mater. Sci.: Mater. Med., 13(8), pp. 783–787. [CrossRef] [PubMed]
Pezzotti, G., and Sakakura, S., 2003, “Study of the Toughening Mechanisms in Bone and Biomimetic Hydroxyapatite Materials Using Raman Microprobe Spectroscopy,” J. Biomed. Mater. Res. Part A, 65(2), pp. 229–236. [CrossRef]
Zhao, L., and Gao, L., 2004, “Novel in Site Synthesis of MWNTS-Hyadoxyapatite Composites,” Carbon, 42, pp. 423–460. [CrossRef]
Ramanathan, T., Abdala, A. A., Stankovich, S., Dikin, D. A., Herrera-Alonso, M., Piner, R. D., Adamson, D. H., Schniepp, H. C., Chen, X., Ruoff, R. S., Nguyen, S. T., Aksay, I. A., Prud'homme, R. K., and Brinson, L. C., 2008, “Functionalized Graphene Sheets for Polymer Nanocomposites,” Nat. Nanotechnol., 3(6), pp. 327–331. [CrossRef] [PubMed]
Rafiee, M., Rafiee, J., Srivastava, I., Wang, Z., Song, H., Yu, Z., and Koratkar, N., 2010, “Fracture and Fatigue in Graphene Nanocomposites,” Small, 6(2), pp. 179–183. [CrossRef] [PubMed]
Rafiee, M., Rafiee, J., Wang, Z., Song, H., Yu, Z., and Koratkar, N., 2009, “Enhanced Mechanical Properties of Nanocomposites at Low Graphene Content,” ACS Nano, 3(12), pp. 3884–3890. [CrossRef] [PubMed]
Cannillo, V., Bondioli, F., Lusvarghi, L., Montorsi, M., Avella, M., Errico, M., and Malinconico, M., 2006, “Modeling of Ceramic Particles Filled Polymer-Matrix Nanocomposites,” Compos. Sci. Technol., 66(7–8), pp. 1030–1037. [CrossRef]
Fisher, F. T., Bradshaw, R. D., and Brinson, L. C., 2002, “Effects of Nanotube Waviness on the Modulus of Nanotube-Reinforced Polymers,” Appl. Phys. Lett., 80(24), pp. 4647–4649. [CrossRef]
Hbaieb, K., Wang, Q., Chia, Y., and Cotterell, B., 2007, “Modelling Stiffness of Polymer/Clay Nanocomposites,” Polymer, 48(3), pp. 901–909. [CrossRef]
Hernández-Pérez, A., and Avilés, F., 2009, “Modeling the Influence of Interphase on the Elastic Properties of Carbon Nanotube Composites,” Comput. Mater. Sci., 47(4), pp. 926–933. [CrossRef]
Li, C., and Chou, T.-W., 2003, “Multiscale Modeling of Carbon Nanotube Reinforced Polymers,” J. Nanosci. Nanotechnol., 3(5), pp. 423–430. [CrossRef] [PubMed]
Liu, Y. J., and Chen, X. L., 2003, “Evaluations of the Effective Material Properties of Carbon Nanotube-Based Composites Using a Nanoscale Representative Volume,” Mech. Mater., 35, pp. 69–81. [CrossRef]
Luo, D., Wang, W., and Takao, Y., 2007, “Effects of the Distribution and Geometry of Carbon Nanotubes on the Macroscopic Stiffness and Microscopic Stresses of Nanocomposites,” Compos. Sci. Technol., 67(14), pp. 2947–2958. [CrossRef]
Wang, H., Zhou, H., Peng, R., and Mishnaevsky, L., 2011, “Nanoreinforced Polymer Composites: 3D FEM Modeling With Effective Interface Concept,” Compos. Sci. Technol., 71(7), pp. 980–988. [CrossRef]
Fertig, R. S., and Garnich, M. R., 2004, “Influence of Constituent Properties and Microstructural Parameters on the Tensile Modulus of a Polymer/Clay Nanocomposite,” Compos. Sci. Technol., 64(16), pp. 2577–2588. [CrossRef]
Sheng, N., Boyce, M. C., Parks, D. M., Rutledge, G. C., Abes, J. I., and Cohen, R. E., 2004, “Multiscale Micromechanical Modeling of Polymer/Clay Nanocomposites and the Effective Clay Particle,” Polymer, 45(2), pp. 487–506. [CrossRef]
Termonia, Y., 2007, “Structure-Property Relationships in Nanocomposites,” Polymer, 48(23), pp. 6948–6954. [CrossRef]
Ji, X., Cao, Y., and Feng, X., 2010, “Micromechanics Prediction of the Effective Elastic Moduli of Graphene Sheet-Reinforced Polymer Nanocomposites,” Model. Simul. Mater. Sci. Eng., 18, p. 045005. [CrossRef]
Liu, H., and Brinson, L., 2008, “Reinforcing Efficiency of Nanoparticles: A Simple Comparison for Polymer Nanocomposites,” Compos. Sci. Technol., 68(6), pp. 1502–1512. [CrossRef]
Shi, D. L., Feng, X. Q., Huang, Y. Y., Hwang, K. C., and Gao, H., 2004, “The Effect of Nanotube Waviness and Agglomeration on the Elastic Property of Carbon Nanotube-Reinforced Composites,” J. Eng. Mater. Technol., 126(3), pp. 250–257. [CrossRef]
Wang, J., and Pyrz, R., 2004, “Prediction of the Overall Moduli of Layered Silicate-Reinforced Nanocomposites—Part I: Basic Theory and Formulas,” Compos. Sci. Technol., 64(7–8), pp. 925–934. [CrossRef]
Lee, C., Wei, X., Kysar, J., and Hone, J., 2008, “Measurement of the Elastic Properties and Intrinsic Strength of Monolayer Graphene,” Science, 321(5887), pp. 385–388. [CrossRef] [PubMed]


Grahic Jump Location
Fig. 1

Sample nanosheet networks used in the finite element simulations showing the effect of the variation of nanosheet network parameters on the models (a) volume fraction, (b) aspect ratio, and (c) nanosheet area

Grahic Jump Location
Fig. 2

Sample nanofiller networks used in the finite element simulations showing different nanofiller orientations. Graphene nanosheets (a) aligned with z-axis, (b) aligned with xy-plane, (c) random distribution. Carbon nanotubes (d) aligned with z-axis, (e) aligned with xy-plane, and (f) random distribution. The orientations that ranged between completely aligned and random cases were investigated in the finite element simulations.

Grahic Jump Location
Fig. 3

(a) Sample finite element mesh showing the different components of the model. The finite element models were generated by combining the graphene nanosheets/CNTs with HA matrix to create the nanocomposite. (b) Nanosheet volume fraction versus elastic modulus for different mesh sizes. Note that Mesh 1 corresponds to the element size used in all the simulations reported in this study. Mesh 2 is a finer mesh with element size that is half of the element size in Mesh 1. Mesh 3 is an even finer mesh that has an element size that is one third of the element size in Mesh 1.

Grahic Jump Location
Fig. 4

Variation of elastic modulus with nanofiller volume fraction for graphene nanosheets and carbon nanotubes. Note that all other variables except volume fraction is kept constant (Graphene nanosheets: area: 1 μm2, aspect ratio: 2, thickness: 1 nm, random orientation; carbon nanotubes: diameter: 1.5 nm, wall thickness: 0.335 nm, length: 100 nm, random orientation)

Grahic Jump Location
Fig. 5

Variation of elastic modulus along all axes with the nanofiller orientation for (a) carbon nanotubes and (b) graphene nanosheets. Note that Orientation 1 and Orientation 2 refer to the in-plane alignment and axis alignment, respectively, whereas xy and z indicate the direction for which the elastic modulus was evaluated. All other variables except alignment is kept constant (graphene nanosheets: volume fraction: 1%, area: 1 μm2, aspect ratio: 2, thickness: 1 nm; carbon nanotubes: volume fraction: 4%, diameter: 1.5 nm, wall thickness: 0.335 nm, length: 40 nm). For HA-graphene nanosheets, the polynomial fits in Orientation 1 are Ex,y = 10−5x3 − 0.0019x2 − 0.0072 x + 110.08 (R2 = 0.995), Ez = 0.0013x2 − 0.0777 x + 102.35 (R2 = 0.971) and in Orientation 2, Ex,y = −10−5x3 + 0.0015x2 − 0.0441 x + 103.87 (R2 = 0.930) and Ez = (2 × 10−5)x3 − 0.0022x2 − 0.0009 x + 109.96 (R2 = 0.980). For HA-carbon nanotubes, the polynomial fits in Orientation 1 are Ex,y = 1 × 10−5x3 − 0.001x2 − 0.0361 x + 108.42 (R2 = 0.992), Ez = 0.001x2 − 0.0614 x + 102.18 (R2 = 0.931) and in Orientation 2, Ex,y = 0.0009x2 − 0.0619 x + 102.12 (R2 = 0.968) and Ez = (6 × 10−5)x3 − 0.0065x2 − 0.1214 x + 126.4 (R2 = 0.998). Note that the variable, x, in the equations denotes the angle of rotation in degrees.

Grahic Jump Location
Fig. 6

Variation of elastic modulus with (a) aspect ratio of graphene nanosheets (volume fraction: 0.67%, area: 1 μm2; thickness: 1 nm; random orientation) (b) thickness of graphene nanosheets (volume fraction: 1%, area: 1 μm2, aspect ratio: 2, random orientation) (c) area of graphene nanosheets (volume fraction: 1%, aspect ratio: 2, thickness: 1 nm, random orientation) (d) length of carbon nanotubes (volume fraction: 4%, diameter: 1.5 nm, wall thickness: 0.335 nm, aligned along z-axis).

Grahic Jump Location
Fig. 7

Comparison of theoretical predictions with finite element simulations for randomly distributed carbon nanotubes and graphene nanosheets. The solid lines for the finite element show the results obtained directly from simulations. Dotted lines extending from the solid lines were included to show the trendline for volume fractions for which simulations were not performed for better comparison with theoretical calculations.




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