Acta Biomaterialia 6 (2010) 1505–1514
Contents lists available at ScienceDirect
Acta Biomaterialia
journal homepage: www.elsevier.com/locate/actabiomat
Mechanistic aspects of the fracture toughness of elk antler bone
M.E. Launey a, P.-Y. Chen b, J. McKittrick b, R.O. Ritchie a,c,* a Materials Sciences Division, Lawrence Berkeley National Laboratory, Berkeley, CA 94720, USA b Materials Science and Engineering Program, Department of Mechanical and Aerospace Engineering, University of California, San Diego, CA 92093, USA c Department of Materials Science and Engineering, University of California, Berkeley, CA 94720, USA article info abstract
Article history: Bone is an adaptive material that is designed for different functional requirements; indeed, bones have a Received 12 September 2009 variety of properties depending on their role in the body. To understand the mechanical response of bone Received in revised form 5 November 2009 requires the elucidation of its structure–function relationships. Here, we examine the fracture toughness Accepted 17 November 2009 of compact bone of elk antler, which is an extremely fast-growing primary bone designed for a totally Available online 24 November 2009 different function than human (secondary) bone. We find that antler in the transverse (breaking) orien- tation is one of the toughest biological materials known. Its resistance to fracture is achieved during crack Keywords: growth (extrinsically) by a combination of gross crack deflection/twisting and crack bridging via Biomaterials uncracked ‘‘ligaments” in the crack wake, both mechanisms activated by microcracking primarily at Elk antler Cortical bone lamellar boundaries. We present an assessment of the toughening mechanisms acting in antler as com- Fracture toughness pared to human cortical bone, and identify an enhanced role of inelastic deformation in antler which Resistance-curves further contributes to its (intrinsic) toughness. Published by Elsevier Ltd. on behalf of Acta Materialia Inc.
1. Introduction mechanical properties [4,5] depending on the growth, structure and adaptation, all of which are interconnected to serve a specific Biological materials are mostly complex systems in which function. large numbers of functionally diverse, and frequently multifunc- We focus here on the fracture resistance of non-structural1 tional, sets of elements interact selectively and nonlinearly to bone, namely the compact bone of elk antler. With the exception produce coherent behavior. One of the most intriguing of these of reindeers, antlers are found only in males, and are grown in the materials is bone, which is a highly hierarchical composite of spring and summer, used in the rut in the fall, and are shed in the assemblies of collagenous protein molecules, water, and mineral winter. Unlike human bone, they provide neither structural support carbonated hydroxyapatite nanoparticles that form a tough, light- nor protection of organs. The functions of antlers are display and weight, adaptive and multi-functional material. Bone is often ste- fighting, with no load-bearing role and low stiffness compared to reotyped as a protective and supportive framework for the body; skeletal bone; however, they are designed to undergo high impact though it performs these functions, it is a dynamic organ that is loading and large bending moments without fracture. constantly remodeling and changing shape to adapt to the forces There have been previous evaluations of the toughness of ant- placed upon it. Like all natural materials, its mechanical proper- ler [4–7], although many of these have been inaccurate due to ties are determined by its structure [1–3], which in turn is problems of inappropriate measurement technique (e.g., mea- defined by its (primarily mechanical) function [4,5]. The adapta- surements based on the area under a compression stress–strain tion of compact bone to its mechanical environment includes curve). In particular, single-value linear-elastic fracture parame- 2 both alteration of its shape and adaptation of its internal struc- ters based on crack initiation, such as KIc, have been used but ture and hence properties. This dual optimization of form and structure is well known in engineering materials; however, in 1 natural materials both are intimately related due to their com- ‘‘Non-structural” refers here to the non-supportive role of antler bones, e.g., antlers do not bear any loads nor support organs. mon origin, the growth of the organ. Different bones grow at dif- 2 For materials that display linear-elastic constitutive behavior, the fracture ferent rates, and the kind of primary bone laid down depends on toughness, Kic (where i = I, II or III), is the critical value of the stress intensity K for this rate of accretion. Accordingly, different bones have different unstable fracture in the presence of a pre-existing crack; under tensile opening ½ conditions (i.e., in mode I) K = Yrapp(pa) = KIc, where rapp is the applied stress, a is the crack length, and Y is a function (of order unity) of crack size and geometry. Alternatively, the toughness can be expressed as a critical value of the strain-energy
* Corresponding author. Tel.: +1 510 486 5798; fax: +1 510 643 5792. release rate, Gc, defined as the rate of change in potential energy per unit increase in E-mail address: [email protected] (R.O. Ritchie). crack area.
1742-7061/$ - see front matter Published by Elsevier Ltd. on behalf of Acta Materialia Inc. doi:10.1016/j.actbio.2009.11.026 1506 M.E. Launey et al. / Acta Biomaterialia 6 (2010) 1505–1514 such measurements cannot capture, or even represent, the multi- terize the toughness associated with both crack initiation and ple length-scale toughening acting in cortical bone that leads to growth.6 We confirm that antler bone is the toughest hard mineral- its characteristic resistance-curve (R-curve3) behavior [8–10], ized tissue reported to date, and provide a description of the where the fracture resistance actually increases with crack exten- toughening mechanisms underlying its exceptional resistance to sion. Antler bone is no exception. Vashishth et al. [11,12] have re- fracture. ported rising R-curve (KR, crack-extension resistance) behavior in antlers of red deer, and demonstrated that the superior toughness 2. Structure and properties of elk Cervus elaphus canadensis of antler bone is due to its enhanced ability to form microcracks antler bone during deformation and fracture. Although such stress-intensity- based R-curves do provide a means to characterize crack propaga- The microstructure of the compact bone of antler is compared tion, the underlying assumptions for such KR calculations are in Fig. 1 with that of human humerus. Elk antler is a young bone based on linear-elastic fracture mechanics (LEFM), which cannot predominantly composed of primary osteons [20] that contain vas- 4 account for the energy associated with plastic deformation dur- cular channels (15–25 lm diameter) surrounded by concentric ing bone fracture (an especially important phenomenon in antler). bone lamellae (Fig. 1a and c). The entire primary osteons are Specifically, for such LEFM measurements, the prevailing mode of 100–200 lm in diameter. In comparison, human bone is a second- deformation is assumed to be linear elasticity; accordingly, any re- ary (replacement) bone that is the product of resorption of previ- gion of ‘‘plasticity” that may form in the vicinity of the crack tip ously existing bone tissue and the deposition of new bone in its (i.e., the plastic zone) must be small enough to ignore. This places place. This process results in the formation of secondary osteons7 restrictions on how large a test specimen has to be for ‘‘valid” that have central vascular channels 50–90 lm in diameter, known toughness measurements; specifically, that the in-plane specimen as Haversian canals (Fig. 1b and d); these are surrounded by a series dimensions of crack size and uncracked ligament width must be of concentric lamellae containing osteocytes arranged in circular at least an order of magnitude larger than the plastic-zone size fashion. The entire secondary osteons (or secondary Haversian sys- (termed ‘‘small-scale yielding”); additionally, for geometry- and tems) are about 200–300 lm in diameter. In antler bone (Fig. 1c), thickness-independent toughness values, the out-of-plane thickness a prominent hyper-mineralized region surrounding the primary ost- dimension must be equally larger than the plastic zone (termed eons [20] is present, whereas in the human bone (Fig. 1d), a thin ‘‘plane strain” conditions). For example for antler bone, a LEFM p mineralized region, the cement line [21], surrounds the secondary KIc value of 10 MPa m [6] would require test specimen dimensions osteons. As prominent sites for microcracking, both have strong (in terms of crack size, ligament depth and thickness) in excess of implications for fracture behavior [22,23]. 50 mm for a valid linear-elastic KIc based on current ASTM validity Although antler has a composition very similar to other mam- criteria [13]. However, because the thickness of the cortical shell in malian long bones, it is the only primary mammalian bone that antler bone is typically 5–10 mm, appropriate section sizes for is capable of regeneration (shedding and re-growth each year). LEFM KIc measurements are not feasible. This means that as the During antlerogenesis (antler growth), the porous core and the vas- test samples used in previous studies in the most part were too cular channels are filled with blood, which provides hydration of small for any form of linear-elastic K measurement, the distribution the bone, along with the naturally occurring water in living bone. of local stresses and displacements near a crack tip (i.e., near the Fully grown antlers had been thought to be dead tissue with all fracture origin) would not be well represented by the K-fields fluid removed once the velvet was shed. More recently, blood-filled [14] and the resulting K-based toughness values would be highly fallow deer antlers, with living osteocytes and active osteoblasts, questionable. Consequently, for materials such as antler that dis- have been found 1 h after casting [24]. play significant plastic deformation prior to fracture [15], LEFM is Antlers exhibit the fastest growth among all natural calcified simply not an appropriate methodology to measure the fracture tissues, growing as much as 14 kg in 6 months, with a peak growth toughness. rate of up to 20–40 mm per day [25]. Such physiological effort of For these reasons, a preferred, indeed essential, strategy to eval- growth necessitate a large import of minerals in a short period of uate the fracture toughness of cortical antler bone is to use nonlin- time, which in turn results in antlers having the lowest mineral ear elastic fracture mechanics. This approach can provide a more content in the bone family (55–60 wt.%) [5], high collagen content, realistic description of the crack-tip stress and displacement fields and consequently low stiffness and yield strength, as compared and furthermore is able to additionally capture the contribution to with human cortical bone (Fig. 2). The organic content, especially the toughness from the energy consumed in ‘‘plastic” deformation the type-I collagen, is also distinctly higher in antler compact bone, prior to and during fracture [16,17]. resulting in stiffness 2–3 times lower than that of human cortical 5 Accordingly, in this work we utilize J-integral measurements to bone [3]; with an ultimate strength of 145–160 MPa in the trans- determine the toughness of elk antler cortical bone using R-curves, verse direction [4,6], this confers more extensibility and a higher in the presence of realistically sized small (<1 mm) cracks, to charac- work to fracture [6] such that it exhibits extensive deformation prior to fracture (a functional adaptation). Indeed, the yield strength of antler compact bone in the transverse direction is as
3 The crack resistance-curve or R-curve provides an assessment of the fracture 6 We note here that although there has been some controversy of late [18,19] of the toughness in the presence of subcritical crack growth. It involves measurements of efficacy of using J-integral methods to characterize the fracture toughness of bone, the the crack-driving force, e.g., K or J, as a function of crack extension (Da). The value of calculation of the value of J at fracture in a bend sample is absolutely identical to that the driving force at Da ? 0 provides a measure of the crack-initiation toughness, of the well-known traditional measure of the toughness in bone, that of the ‘‘work of whereas the slope or the maximum value of the R-curve can be used to characterize fracture”, i.e., the energy involved in the fracture process (area under the load– the crack-growth toughness. displacement curve) divided by twice the area of the fracture surface. The only 4 Plastic deformation here is used as a general term to indicate any of the inelastic, difference is that specimens used for work of fracture measurements may not contain non-recoverable deformation mechanisms, such as local collagen fibrillar shearing, a pre-crack or notch. viscoplasticity, and microcracking, that are active at various length-scales in bone. 7 Antlers undergo limited secondary osteon remodeling and consist mainly of 5 J is the nonlinear strain-energy release rate, i.e., the rate of change in potential primary osteons [20]. Secondary osteon remodeling arises in response to mechanical energy for a unit increase in crack area in a nonlinear elastic solid. It is the nonlinear stress and takes about 2–4 months in human bone [3]. In antlers, secondary osteon elastic equivalent of the strain-energy release rate G. It characterizes the stress and remodeling is unlikely to occur since they do not sustain mechanical loads during the displacement fields at a crack tip in such a solid, and as such can be used to define the growth process and are only used in sporadic combat for 1–2 months before onset of fracture there. shedding. M.E. Launey et al. / Acta Biomaterialia 6 (2010) 1505–1514 1507
Fig. 1. Microstructure of primary and secondary bone. Differential interference contrast (Nomarski) micrograph of transverse section of compact (cortical) bone of (a) elk antler, and (b) human humerus. (c) and (d) Back-scattered SEM images. Morphologically, the main distinction between primary (c) and secondary bone (d) is that primary osteons do not have cement lines because they are not the product of bone remodeling; however, the interfaces of the primary osteons in antler are thick hyper-mineralized regions. Primary osteons have smaller vascular channels and fewer lamellae than secondary osteons.
of 12 mm. Six samples of each orientations were taken from loca- tions longitudinal or transverse to the bone long axis (Fig. 3a). An ini- tial notch was applied with a low-speed diamond saw and was subsequently sharpened by repeatedly sliding a razor blade over the saw-cut notch, while continually irrigating with 1 lm diamond slurry. The final micro-notches had a root radius of 3–5 lm. As a result, sharp cracks with initial crack length, a 1.5 mm (a/ W 0.5), were generated in general accordance with ASTM stan- dards [13]. The orientation of the notch was such that the nominal crack-growth direction was either perpendicular to the long axis of the antler (transverse orientation), along the long axis of the antler (in-plane longitudinal), and parallel to the long axis of the antler but perpendicular to the (nominal) crack-propagation direction (anti- plane longitudinal). Each set of samples was further divided into two groups, three to be tested ex situ and three to be tested in situ inside the scanning electron microscope. Prior to testing, all samples were wet polished with an increasingly higher finish to a final polish Fig. 2. Stress–strain curves from three-point bending tests for hydrated human with a 0.05 lm diamond suspension before being immersed in cortical bone and elk antler in the transverse orientation. Data for elk antler are ambient Hanks’ balanced saline solution (HBSS) for 24 h. taken from Ref. [6].
3.2. Fracture toughness J–R curve measurements low as 71 MPa [6], while the corresponding yield strength of hu- man cortical bone is 110–120 MPa (Fig. 2). J–R curves for compact bone of elk antler were performed under rehydrated conditions in mode I (tensile opening) using single- 3. Experimental procedures edge notched bend SE(B) specimens with a crack-growth direction (Fig. 3a): (i) transverse to the long axis of the osteons (transverse), 3.1. Materials (ii) along the direction parallel to the long axis of the osteons (in- plane longitudinal), and (iii) along the direction parallel to the long Test samples from the compact region of North American elk axis of the osteons but perpendicular to the (nominal) crack-prop- (Cervus elaphus canadensis) were sectioned using a low-speed saw agation direction (anti-plane longitudinal). and machined into eighteen bend samples (N = 18). The antler, from R-curves were measured ex situ in 25 °C HBSS to evaluate the a large, mature bull, was shed approximately one year before testing fracture resistance in terms of the J-integral as a function of crack and stored indoors under air-dry condition. Rectangular samples extension, Da, under a monotonically increasing driving force. had a thickness B of 2.0–2.2 mm, a width W of 3 mm, and a length Tests were conducted in three-point bending with a span 1508 M.E. Launey et al. / Acta Biomaterialia 6 (2010) 1505–1514
Fig. 3. The schematic (a) shows the structure of antler bone, as well as the three anatomical orientations that the specimens were taken from the compact bone. (b) Typical load vs. load-line displacement curve obtained during R-curve testing for the transverse, in-plane longitudinal, and anti-plane longitudinal orientation of antler compact bone. Each partial unloading event during the test corresponds to a data point in (c). (c) Full JR(Da) resistance curves for stable ex situ crack extension in hydrated antler compact bone tested in transverse, in-plane longitudinal, and anti-plane longitudinal orientations. (d) Expansion of (c) showing JR(Da) resistance curves at larger crack extension (Da > 0.4 mm). The R-curves for short crack lengths (Da 0.6 mm) are compared with data taken from human cortical bone in both transverse and in-plane longitudinal orientations of the humerus. Data for human cortical bone were adapted from Ref. [17]. In the transverse orientation, fracture toughnesses, J, of antler compact 2 2 bone were recorded at nearly 60 kJ m , almost two times higher than the critical toughness value, Jc, of human cortical bone (Jc 30 kJ m ). The validity of these data points is defined by the measurement capacity of each specimen in accordance with the ASTM standard [13]. The circles in (c) and (d) correspond to data points for antler bone.
(S = 10 mm) to width (W = 3 mm) ratio of 3, in accordance with Here C is the sample compliance, and F is a calibration factor, ASTM E1820-08 [13]. The specimens were loaded in displacement taken to be that which gives the best agreement between the ini- control in a MTS 810 servo-hydraulic testing machine with a loading tial compliance and crack length at the beginning of the test. rate of 0.015 mm s 1 until the onset of cracking, which was deter- In addition, R-curves were measured on HBSS-saturated speci- mined by non-linearity in the load–displacement curve (Fig. 3b). To mens in situ in a Hitachi S-4300SE/N environmental scanning monitor subsequent subcritical crack growth, after this point during electron microscope (ESEM) using a Gatan Microtest three-point loading, the sample was periodically unloaded (by 10–20% of the bending stage. Crack extension was monitored directly in back- peak load) to record the elastic load-line compliance using a LVDT scattered electron mode at a pressure of 35 Pa and a 30 kV excita- mounted on the load frame. After each increment, the specimens tion voltage. were held for 30 s to allow for crack extension to stabilize, followed R-curve determination was limited to small-scale bridging by unloading compliance measurement. This process was repeated conditions, where the size of the zone of crack bridges behind at regular intervals until the end of the test, at which point the com- the crack tip remained small compared to the in-plane test spec- pliance and loading data were analyzed to determine J-integral as a imen dimensions. As noted above, the use of the J-integral as the function of Da. Crack lengths, a, were calculated from the compli- driving force for crack initiation and growth was employed to ance data obtained during the test using compliance expression of capture the contribution from inelastic deformation in the evalu- a three-point bend specimen at load line [26]: ation of toughness. The stress intensity at each measured crack length was calculated by measuring the nonlinear strain-energy a=W 0:997 3:58U 1:51U2 110U3 1232U4 4400U5; ¼ þ release rate, J. The value of J was calculated from the applied load ð1Þ and instantaneous crack length according to ASTM standards where U is a fitting function, written as: [13], and was decomposed into its elastic and plastic contribu- tions: ffiffiffiffiffiffi1 U ¼ p ð2Þ : FC þ 1 J ¼ Jel þ Jpl ð3Þ M.E. Launey et al. / Acta Biomaterialia 6 (2010) 1505–1514 1509