Microsoft Word - numero_37_art_28 S. Vantadori et alii, Frattura ed Integrità Strutturale, 37 (2016) 215-220; DOI: 10.3221/IGF-ESIS.37.28 215 Focussed on Multiaxial Fatigue and Fracture Two-parameter fracture model for cortical bone Sabrina Vantadori, Andrea Carpinteri, Giovanni Fortese, Camilla Ronchei, Daniela Scorza University of Parma, Dept. of Civil-Environmental Engineering and Architecture, Parco Area delle Scienze 181/a, 43124 Parma – Italy sabrina.vantadori@unipr.it, http://orcid.org/0000-0002-1904-9301 Filippo Berto University of Padua, Dept. of Management and Engineering, Stradella San Nicola 3, 36100, Vicenza - Italy NTNU, Department of Engineering Design and Materials, Richard Birkelands vei 2b, 7491, Trondheim – Norway ABSTRACT. The analysis of the bone fracture behaviour is fundamental for prevention, diagnosis and treatment of traumas. In the present paper, an experimental analysis of the fracture behaviour of a bovine femoral cortical bone is carried out, where specimens are extracted from a diaphysis. Fracture toughness is computed by employing a two-parameter fracture model originally proposed for concrete. In order to take into account the possibility of crack deflection (kinked crack) due to osteons orientation, a modified version of such a model is here discussed. The fracture toughness results are then compared with those reported in the literature, related to a femur of an 18-month-old bovine. KEYWORDS. Cortical Bone; Fracture Toughness; Quasi-brittle Material; Two-parameter Model. INTRODUCTION one is a specialised tissue which has both metabolic and mechanical functions [1-3]. The load-bearing capacity of bone is limited up to a certain extent, beyond which it fails [4]. The analysis of the fracture behaviour of bone is fundamental for prevention, diagnosis and treatment of traumas. Basic parameters which represent the structure and functions of bone have to be measured, such as its fracture toughness. In the present paper, an experimental analysis of the fracture behaviour of a bovine femoral cortical bone is carried out, where specimens are extracted from a diaphysis. The experimental campaign is conducted to determine the fracture toughness, by employing a two-parameter fracture model originally proposed for concrete [5], that is, for a quasi-brittle material showing a nonlinear slow crack growth before the peak load is reached. In bones, such a behaviour is produced by the mechanism of extrinsic toughening categorised in four classes [6]: (i) constrained microcracking; (ii) crack deflection and twist; (iii) uncracked-ligament bridging; (iv) collagen-fibril bridging. Such a two-parameter fracture method is based on experimental data obtained from three-point bending tests on single edge-notched specimens, and employs linear elastic fracture mechanics Mode I expressions. However, for the bone material, it cannot be applied in its original formulation since the crack starting from notch may deflect. In order to understand the cause of such a deflection under Mode I loading (three-point bending), the bone microstructure level has to be briefly examined [7,8]. In cortical bone, it is represented by osteons (Fig. 1). The osteons B S. Vantadori et alii, Frattura ed Integrità Strutturale, 37 (2016) 215-220; DOI: 10.3221/IGF-ESIS.37.28 216 are oriented parallel to the bone axis, which consists of a vascular canal (or Haversian canal) surrounded by concentric lamellae. The interface between osteons and interstitial lamellae is called cement line (Fig. 1). When the osteons alignment is perpendicular to the loading direction (transversal specimens), the stress state under three- point bending loading is biaxial due to normal stresses produced by bending and shear stresses at the cement line interface between osteons and interstitial lamellae [9]. In such a case, the crack starting from the notch is subjected to Mixed Mode loading (Mode I and Mode II). On the other hand, when the osteons alignment is parallel to the loading direction (longitudinal specimens), the stress state is uniaxial. In such a case, the crack starting from the notch is subjected to Mode I loading. The Two-Parameter Model can be directly applied for the latter case, whereas a modified formulation is hereafter proposed for the former case in order to take into account that the crack is also subjected to Mode II loading. More precisely, the unloading compliance expression related to a single kinked crack is determined by employing the Castigliano theorem in the manner suggested by Paris [10], and the effective crack length is computed. Then, the critical stress- intensity factor at the crack tip is evaluated as a function of that related to a straight crack having length equal to the projected length of the effective kinked crack [11,12]. vascular canal interstitial lamellae cement line osteon Figure 1: Cortical bone microstructure level. Finally, the fracture toughness results here obtained are compared with those related to a femur of an 18-month-old bovine [13], experimentally determined according to the standard ASTM E399-1 [14]. TWO-PARAMETER MODEL ccording to the Two-Parameter Model, the specimens present a notch in the lower part of the middle cross section (Fig. 2). The three-point bending tests are performed under crack mouth opening displacement control (average speed equal to 0.1 mm h-1). W L S B a a0 (a) W L S B a0 (b) a1 a2  Figure 2: Crack propagates under: (a) pure Mode I; (b) Mixed Mode. Each specimen is monotonically loaded: after the peak load is achieved, the post-peak stage follows and, when the force is equal to about 95% of the peak load, the specimen is fully unloaded. Then, the specimen is reloaded up to failure. A S. Vantadori et alii, Frattura ed Integrità Strutturale, 37 (2016) 215-220; DOI: 10.3221/IGF-ESIS.37.28 217 The initial compliance, iC , is used to determine the elastic modulus, E [15]:   BWC VaS E i 2 006  (1) where S, W and B are the loading span, depth and thickness of the specimen, respectively, 0a is the notch length, and iC is the linear elastic compliance. Further, the parameter V can be expressed as follows [15]:   W a V 002 0 3 0 2 000 with 1 66.0 04.287.328.276.0       (2) Therefore, if the crack propagates under pure Mode I, the effective critical crack length, a , is determined from the following equation by employing an iterative procedure:   BWC VaS E u 2 6   (3) where uC is the unloading compliance, and  V is obtained from Eq. 2 by replacing 0a with a . Finally, the Mode I critical stress-intensity factor, SICK , is computed by employing the measured value of the peak load, maxP [15]: )( 2 3 2 max  fa BW SP K SIC  (4) where:    W a f          with 121 )70.293.315.2()1(99.11 )( 2/3 2 (5) MODIFIED TWO-PARAMETER MODEL modified procedure is hereafter proposed when crack propagates under Mixed Mode loading (Mode I and Mode II). Specimens geometry and experimental test procedure are equal to those presented in the previous Section (see Fig. 2(b)). The elastic modulus is determined through Eq. 1. Under Mixed Mode loading, the effective critical crack length, 210 aaaa  , is obtained from the following equation by employing an iterative procedure: )]}()( )][([ ])[(][ { cos cos coscos coscoscossincos cos cos 2 cos 2 sin 2 cos )( 6 10 10 210 210 23 00 10 10 426 002 )( )()( W aa Vaa W aaa Vaaa Va W aa Vaa Va BWC S E u                   (6) A S. Vantadori et alii, Frattura ed Integrità Strutturale, 37 (2016) 215-220; DOI: 10.3221/IGF-ESIS.37.28 218 Eq. (6) is deduced by employing the Castigliano theorem in the manner suggested by Paris [10], being  the crack kinking angle (Fig. 2(b)) and 01 3.0 aa  . Note that, as is shown in Fig. 2(b), the kinked crack path consists of the two segment, named 1a and 2a . If the value of 2a obtained from Eq. 6 is negative, it means that the effective crack length is 10 aaa  with 01 3.0 aa  . Such a length is obtained from the following equation by employing an iterative procedure: ]})[(][{ )()( 001010426002 cos cos 2 cos 2 sin 2 cos)( 6      Va W aa VaaVa BWC S E u      (7) Finally, the critical stress-intensity factor, SICK , is computed through Eqs 4 and 5, by considering a straight crack having length equal to the projected length of the effective kinked crack:        01 210 2102 max 3.0when cos with cos 2 3 aa W aaa faaa BW P K SIC        (8a) or     01 10 102 max 3.0when cos with cos 2 3 aa W aa faa BW P K SIC        (8b) RESULTS AND DISCUSSION alues of fracture toughness, SICK , of the analysed cadaveric femur diaphysis of a 24-month-old bovine are computed with respect to the osteons alignment: perpendicular or parallel to the loading direction. The former specimens are extracted from anterior (FA-1 and FA-2 in Tab. 1), posterior (FP in Tab. 1), medial (FM-1 and FM-2 in Tab. 1) and lateral (FL in Tab. 1) cortical bone, whereas the latter specimens from posterior cortical bone (FP- long1 and FP-long2 in Tab. 1). All specimens exhibit a non-linear slow crack growth before the peak load is reached. It can be observed that the fracture toughness values for specimens characterised by different osteons alignment with respect to the loading direction are significantly different. As a matter of fact, when the osteons alignment is perpendicular to the loading direction, crack grows under Mixed Mode (Fig. 3(a)), and a higher resistance to fracture is observed. When the osteons alignment is parallel to the loading direction, crack grows under Mode I (Fig. 3(b)), and a lower resistance to fracture is observed. For transversal specimens, the average value of SICK ( mMPa15.087.3  ) is in the range of the cortical bone fracture values [16]. Such a value is then compared with those determined by Libonati et al. [13] according to ASTM standards [14], following the LEFM and considering cracks under pure Mode I loading. They found mMPaK SIC 1.06.5  from SE(B) type and mMPaK SIC 6.08.5  from C(T) type. The difference of such results with respect to those here obtained is due to the fact that the reduction of fracture toughness for a kinked crack as compared with the straight counterpart has not been taken into account by Libonati et al. in Ref. [13]. The present study highlights that the value of the near-tip stress-intensity factor of a kinked crack can be considerably lower than that for a straight crack of the same length, and that has to be taken into account for prevention, diagnosis and treatment of bone traumas. V S. Vantadori et alii, Frattura ed Integrità Strutturale, 37 (2016) 215-220; DOI: 10.3221/IGF-ESIS.37.28 219 Specimens No. Elastic modulus E (MPa) Fracture toughness SICK (MPa m1/2) FA-1 13692.54 4.09 FA-2 14393.65 3.68 FP 14038.29 3.74 FM-1 14253.00 3.92 FM-2 13488.94 3.85 FL 14042.47 3.91 FP-long1 12371.99 2.14 FP-long2 12370.47 1.99 Table 1: Elastic modulus and fracture toughness.           (a) FA-1 FL (b) FP-long1 FP-long2 Figure 3: Fracture: (a) Mixed Mode for transversal specimens FA-1 and FL; (b) pure Mode I for longitudinal specimens. CONCLUSIONS n the present paper, the behaviour of a compact bone in terms of fracture toughness has been analysed. Fracture toughness has experimentally been evaluated through specimens obtained from the femur diaphysis of a bovine. The influence of local biaxial stress state on fracture behaviour has been analysed by employing two specimen types. As a matter of fact, when the osteons alignment is perpendicular to the loading direction (transversal specimens), the stress state is biaxial due to normal stresses produced by bending and shear stresses, at the cement line interface between osteons and interstitial lamellae. On the other hand, when the osteons alignment is parallel to the loading direction (longitudinal specimens), the stress state is uniaxial. Then fracture toughness values have been computed by a modified version of the Two-Parameter Model originally formulated for crack propagating under Mode I. Such a modified version is here proposed for Mixed Mode (Mode I and Mode II). The theoretical results obtained are compared with some data available in the literature, by highlighting that the value of the near-tip stress-intensity factor of a kinked crack can be considerably lower than that for a straight crack of the same length. REFERENCES [1] An, H.Y., Draughn, R.A., Mechanical testing of bone and the bone-implant interface, CRC Press, Boca Raton, (2000). [2] Marks, S.C., Popoff, S.N., Bone cell biology: the regulation of development, structure, and function in the skeleton, Am. J. Anat., 183 (1988) 1-44. [3] Keaveny, T.M., Hayes, W.C., Mechanical properties of cortical and trabecular bone, CRC Press, Boca Raton, (1993). [4] Li, S., Abdel-Wahab, A., Silberschmidt, V.V., Analysis of fracture processes in cortical bone tissue, Eng. Frac. Mech., 110 (2013) 448-458. DOI: 10.1016/j.engfracmech.2012.11.020 [5] Jenq, Y., Shah, S., J., Two parameter fracture model for concrete, Eng. Mech., 111 (1985) 1227-1241. I S. Vantadori et alii, Frattura ed Integrità Strutturale, 37 (2016) 215-220; DOI: 10.3221/IGF-ESIS.37.28 220 [6] Nalla, R.K., Kinney, J.H., Ritchie, R.O., Mechanistic fracture criteria for the failure of human cortical bone, Nat. Mater., 2 (2003) 164-168. DOI: 10.1038/nmat832 [7] Katz, J.L., The structure and biomechanics of bone, Symp. Soc. Exp. Biol., 34 (1980) 137-168. [8] Rho, J.Y., Kuhn Spearing, L., Zioupos, P., Mechanical properties and the hierarchical structure of bone, Med. Eng. & Phys., 20 (1998) 92-102. DOI: 10.1016/S1350-4533(98)00007-1 [9] Martin, R.B., Burr, D.B., Sharkey, N.A., Skeletal tissue mechanics, Springer-Verlag, New York, (1998). [10] Paris, P.C., The mechanics of fracture propagation and solutions to fracture arrester problems, Document D2-2195, The Boeing Company, (1957). [11] Kitagawa, H., Yuuki, R., Ohira, T., Crack-morphological aspects in fracture mechanics, Eng. Frac. Mech., 7 (1975) 515-529. DOI: 10.1016/0013-7944(75)90052-1 [12] Cotterell, B., Rice, J.R., Slightly curved or kinked cracks, Int. J. Frac., 16 (1980) 155-169. DOI: 10.1007/BF00012619 [13] Libonati, F., Vergani, L., Bone toughness and crack propagation: An experimental study, Proc. Engin., 74 (2014) 464- 467. DOI: 10.1016/j.proeng.2014.06.298 [14] ASTM E399-1, Standard Test Method for Linear-Elastic Plane-Strain Fracture Toughness KIc of Metallic Materials. [15] Tada, H., Paris, P.C., Irwin, G.R., The stress analysis of cracks handbook, ASME Press, New York, (2000). [16] Khanal, S.P., Mahfuz, H., Rondinone, A.J., Leventouri, T., Improvement of the fracture toughness of hydroxyapatite (HAp) by incorporation of carboxyl functionalized single walled carbon nanotubes (CfSWCNTs) and nylon, Mat. Sci. Eng. C, 60 (2016) 204-210. DOI: 10.1016/j.msec.2015.11.030 << /ASCII85EncodePages false /AllowTransparency false /AutoPositionEPSFiles true /AutoRotatePages /None /Binding /Left /CalGrayProfile (Dot Gain 20%) /CalRGBProfile (sRGB IEC61966-2.1) /CalCMYKProfile (U.S. Web Coated \050SWOP\051 v2) /sRGBProfile (sRGB IEC61966-2.1) /CannotEmbedFontPolicy /Error /CompatibilityLevel 1.4 /CompressObjects /Tags /CompressPages true /ConvertImagesToIndexed true /PassThroughJPEGImages true /CreateJobTicket false /DefaultRenderingIntent /Default /DetectBlends true /DetectCurves 0.0000 /ColorConversionStrategy /CMYK /DoThumbnails false /EmbedAllFonts true /EmbedOpenType false /ParseICCProfilesInComments true /EmbedJobOptions true /DSCReportingLevel 0 /EmitDSCWarnings false /EndPage -1 /ImageMemory 1048576 /LockDistillerParams false /MaxSubsetPct 100 /Optimize true /OPM 1 /ParseDSCComments true /ParseDSCCommentsForDocInfo true /PreserveCopyPage true /PreserveDICMYKValues true /PreserveEPSInfo true /PreserveFlatness true /PreserveHalftoneInfo false /PreserveOPIComments true /PreserveOverprintSettings true /StartPage 1 /SubsetFonts true /TransferFunctionInfo /Apply /UCRandBGInfo /Preserve /UsePrologue false /ColorSettingsFile () /AlwaysEmbed [ true ] /NeverEmbed [ true ] /AntiAliasColorImages false /CropColorImages true /ColorImageMinResolution 300 /ColorImageMinResolutionPolicy /OK /DownsampleColorImages true /ColorImageDownsampleType /Bicubic /ColorImageResolution 300 /ColorImageDepth -1 /ColorImageMinDownsampleDepth 1 /ColorImageDownsampleThreshold 1.50000 /EncodeColorImages true /ColorImageFilter /DCTEncode /AutoFilterColorImages true /ColorImageAutoFilterStrategy /JPEG /ColorACSImageDict << /QFactor 0.15 /HSamples [1 1 1 1] /VSamples [1 1 1 1] >> /ColorImageDict << /QFactor 0.15 /HSamples [1 1 1 1] /VSamples [1 1 1 1] >> /JPEG2000ColorACSImageDict << /TileWidth 256 /TileHeight 256 /Quality 30 >> /JPEG2000ColorImageDict << /TileWidth 256 /TileHeight 256 /Quality 30 >> /AntiAliasGrayImages false /CropGrayImages true /GrayImageMinResolution 300 /GrayImageMinResolutionPolicy /OK /DownsampleGrayImages true /GrayImageDownsampleType /Bicubic /GrayImageResolution 300 /GrayImageDepth -1 /GrayImageMinDownsampleDepth 2 /GrayImageDownsampleThreshold 1.50000 /EncodeGrayImages true /GrayImageFilter /DCTEncode /AutoFilterGrayImages true /GrayImageAutoFilterStrategy /JPEG /GrayACSImageDict << /QFactor 0.15 /HSamples [1 1 1 1] /VSamples [1 1 1 1] >> /GrayImageDict << /QFactor 0.15 /HSamples [1 1 1 1] /VSamples [1 1 1 1] >> /JPEG2000GrayACSImageDict << /TileWidth 256 /TileHeight 256 /Quality 30 >> /JPEG2000GrayImageDict << /TileWidth 256 /TileHeight 256 /Quality 30 >> /AntiAliasMonoImages false /CropMonoImages true /MonoImageMinResolution 1200 /MonoImageMinResolutionPolicy /OK /DownsampleMonoImages true /MonoImageDownsampleType /Bicubic /MonoImageResolution 1200 /MonoImageDepth -1 /MonoImageDownsampleThreshold 1.50000 /EncodeMonoImages true /MonoImageFilter /CCITTFaxEncode /MonoImageDict << /K -1 >> /AllowPSXObjects false /CheckCompliance [ /None ] /PDFX1aCheck false /PDFX3Check false /PDFXCompliantPDFOnly false /PDFXNoTrimBoxError true /PDFXTrimBoxToMediaBoxOffset [ 0.00000 0.00000 0.00000 0.00000 ] /PDFXSetBleedBoxToMediaBox true /PDFXBleedBoxToTrimBoxOffset [ 0.00000 0.00000 0.00000 0.00000 ] /PDFXOutputIntentProfile () /PDFXOutputConditionIdentifier () /PDFXOutputCondition () /PDFXRegistryName () /PDFXTrapped /False /CreateJDFFile false /Description << /ARA /BGR /CHS /CHT /CZE /DAN /DEU /ESP /ETI /FRA /GRE /HEB /HRV (Za stvaranje Adobe PDF dokumenata najpogodnijih za visokokvalitetni ispis prije tiskanja koristite ove postavke. Stvoreni PDF dokumenti mogu se otvoriti Acrobat i Adobe Reader 5.0 i kasnijim verzijama.) /HUN /ITA /JPN /KOR /LTH /LVI /NLD (Gebruik deze instellingen om Adobe PDF-documenten te maken die zijn geoptimaliseerd voor prepress-afdrukken van hoge kwaliteit. De gemaakte PDF-documenten kunnen worden geopend met Acrobat en Adobe Reader 5.0 en hoger.) /NOR /POL /PTB /RUM /RUS /SKY /SLV /SUO /SVE /TUR /UKR /ENU (Use these settings to create Adobe PDF documents best suited for high-quality prepress printing. Created PDF documents can be opened with Acrobat and Adobe Reader 5.0 and later.) >> /Namespace [ (Adobe) (Common) (1.0) ] /OtherNamespaces [ << /AsReaderSpreads false /CropImagesToFrames true /ErrorControl /WarnAndContinue /FlattenerIgnoreSpreadOverrides false /IncludeGuidesGrids false /IncludeNonPrinting false /IncludeSlug false /Namespace [ (Adobe) (InDesign) (4.0) ] /OmitPlacedBitmaps false /OmitPlacedEPS false /OmitPlacedPDF false /SimulateOverprint /Legacy >> << /AddBleedMarks false /AddColorBars false /AddCropMarks false /AddPageInfo false /AddRegMarks false /ConvertColors /ConvertToCMYK /DestinationProfileName () /DestinationProfileSelector /DocumentCMYK /Downsample16BitImages true /FlattenerPreset << /PresetSelector /MediumResolution >> /FormElements false /GenerateStructure false /IncludeBookmarks false /IncludeHyperlinks false /IncludeInteractive false /IncludeLayers false /IncludeProfiles false /MultimediaHandling /UseObjectSettings /Namespace [ (Adobe) (CreativeSuite) (2.0) ] /PDFXOutputIntentProfileSelector /DocumentCMYK /PreserveEditing true /UntaggedCMYKHandling /LeaveUntagged /UntaggedRGBHandling /UseDocumentProfile /UseDocumentBleed false >> ] >> setdistillerparams << /HWResolution [2400 2400] /PageSize [612.000 792.000] >> setpagedevice