Microsoft Word - numero_63_art_07_3908.docx P. Livieri et alii, Frattura ed Integrità Strutturale, 63 (2023) 71-79; DOI: 10.3221/IGF-ESIS.63.07 72 The Stress Intensity Factor of convex embedded polygonal cracks Paolo Livieri, Fausto Segala University of Ferrara, Department of Engineering, Italy paolo.livieri@unife.it, fausto.segala@unife.it ABSTRACT. In the present work, a simple formula for the evaluation of the stress intensity factor (SIF) of convex embedded polygonal cracks has been proposed. This formula is structured as a correction factor of the Oore-Burns’ equation and is based on accurate three-dimensional FE analysis. Furthermore, a precise formula for a regular polygonal crack has been given. KEYWORDS. Weight function, Stress intensity factor, Three-dimensional crack. Citation: Livieri, P., Segala, F., The Stress Intensity Factor of convex embedded polygonal cracks, Frattura ed Integrità Strutturale, 63 (2023) 72-79. Received: 11.09.2022 Accepted: 15.10.2022 Online first: 25.10.2022 Published: 01.01.2023 Copyright: © 2023 This is an open access article under the terms of the CC-BY 4.0, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited. INTRODUCTION he weight function technique was performed by Bueckner [1] and Rice [2, 3] to solve the stress intensity factor (SIF) of a crack as well as to find out the nominal stress over the geometric discontinuity. It is well known that the O- integral [4] is a good approximation to evaluate the SIF along the contour of three-dimensional planar cracks. For example, British Standard 7910 [5] suggests the use of the O-integral provided that the results are documented. Usually, in order to simplify SIF assessments for embedded flaws, the SIF is determined at the ends of both the minor and major axes of the elliptical idealization of the flaw [6, 7, 8]. The SIF of a semi elliptical crack can be accurately calculated by means of the general procedure of Shen and Glinka [9, 10]. They use four terms of approximation and evaluate the unknown parameters based on reference stress intensity factor expressions taken from the literature. In this way a general weight function can be formulated to overcome some problems [11, 12] in the Petrosky and Achenbach method [13] where an approximate displacement field from a known reference was calculated. Also in Fitness-for-Service procedure crack shapes, idealisation becomes necessary when real cracks have been detected during an inspection. Typically, elliptical cracks, semi-elliptical cracks and through wall cracks or edge cracks with a rectilinear flank are considered [14]. For convex three-dimensional cracks few examples are presented in classical textbooks [15, 16, 17]. Furthermore, in order to obtain an acceptable approximation of the maximum SIF KImax, Murakami [18, 19] took into account many types of convex embedded cracks subjected to nominal tensile stress σn. The SIF was investigated numerically, and a final equation for SIF was given in the form:  ,maxI nK Y A where Y is a coefficient that was evaluated as best fitting of the numerical results; and A is the area of the flaw. T https://youtu.be/DRqmAxPDU70 P. Livieri et alii, Frattura ed Integrità Strutturale, 63 (2023) 71-79; DOI: 10.3221/IGF-ESIS.63.07 73 In a previous paper [20] the authors investigated the failure of the O-integral in presence of high curvature for regular (  of class C2) cracks. In terms of regularity, a corner Q’ means a singularity. This requires a new technique in order to correct the O-integral. In particular, in this paper we are interested in adjusting the O-integral for convex polygons. The analytical results will be compared with numerical ones obtained from an accurate three-dimensional FE analysis. WEIGHT FUNCTION FOR A THREE-DIMENSIONAL CRACK: ANALYTICAL BACKGROUND ore and Burns proposed a general equation for the evaluation of the SIF of a two-dimensional crack inside a three-dimensional body subjected to a nominal tensile stress σn(Q). The nominal stress σn(Q) is evaluated without the presence of the crack. Q is the inner point of the crack. The crack can be considered as an open bounded simply connected subset Ω of the plane as reported in Fig. 1. We define:     2( ) ( ) ds f Q Q P s (1) where   ( , )Q Q x y , s is the arch-length parameter and point P(s) runs over the boundary  . In reference [4] Oore- Burns proposed an empirical formula for the evaluation of the mode I stress intensity factor at each point of the border crack of boundary  :       , 2 ( )2 ( ') , ' ( ) ' n I OB Q K Q d Q f Q Q Q (2) Under reasonable hypotheses on the function σn(Q), the integral (2) is convergent and the proof is based on the asymptotic behaviour of f(Q) [21]. In different papers, the authors analysed the properties of the Oore-Burns integral where the accuracy of the equation was tested in the particular case of an elliptical crack [22, 23]. Figure 1: Inner crack. CONVEX POLYGON n a previous paper [20] we investigated the proprieties of the O-integral in the presence of high curvature for regular ( of class C2) cracks. Henceforth, we denote by =(Q’) the correction factor of Eqn. (2) at Q’  ,( ') (Q') ( ')I I OBK Q K Q (3) In order to draw out the factor , we need to take into account some hints: O I P. Livieri et alii, Frattura ed Integrità Strutturale, 63 (2023) 71-79; DOI: 10.3221/IGF-ESIS.63.07 74 A) a regular polygon with very large N sides, is “indistinguishable” from a disk, of which the SIF is given by   2 n times the square root of the disk radius; B)  is close to one away from the corners; C) known FE results for cracks with a different shape: square, equilateral triangles and rectangular (aspect ratio 1/3). D) the crack is subjected to a uniform tensile stress σn. The key to construct (Q’) is a smart use of the hyperbolic tangent function xtanh(x). At an early stage, the factor  will contain some unknown parameters that will be carefully calibrated on the basis of requirements A, B and C. Let us assume p is the perimeter of the convex polygon  and c is the perimeter of the smallest disk containing . Near a fixed edge P with opening angle α (see Fig. 2, where     0   ,   x Q P ), (Q’) can be given as follows:                                    0Q’     1   1 tanh     2  xc c p p p (4) with , ,  and  chosen in order to satisfy A), B) and C) conditions. On the basis of accurate FE analysis on a three-dimensional model as proposed in reference [20], the best agreement is given by the choice   1   10 , =6.95, =0.8 and =0.4. This means that near the corner, the coefficient (Q’) will be close to the value:                                0.8 0.40.8 01Q’     1   1 tanh 6.95  10 2  xc c p p p (5) Now, it is possible to extend Eqn. (5) to entire contour  in a natural way, by taking into account all distances (on the geometry of ) of Q’ from each corner of the polygon:                                                           0.40.80.8 1 0.8 3/2 1 1 1   1   tanh 6.95  10 2  Q’        1 1   1 tanh 6.95  10 2    N k k N k k Q Pc c p p p c Q Pc p p (6) where Pk, k=1, 2, ..., N are the corners with opening angles αk, k=π- αk and  Q P is the distance between Q’ and Pk on the boundary . Figure 2: Polygonal crack. P. Livieri et alii, Frattura ed Integrità Strutturale, 63 (2023) 71-79; DOI: 10.3221/IGF-ESIS.63.07 75 For example, we explain Eqn. (6) on an equivalent triangle with side L, where       0   ,  2 3 L Q y and  0    0 0.5 y L as shown in Fig. (3):             0.4 01Q’    1.046 tanh 4.14      2 y L ·                         0.4 0.4 0 03 1tanh 4.14      ·tanh 4.14      2 2 y y L L (7) Clearly, in the r.h.s. of the Eqn. (7), the only significant contribution is given by             0.4 01Q’    1.046 tanh 4.14      2 y L (8) Figure 3: Equilateral triangular crack. Figure 4: Square crack. As a second example, Fig. 4 shows the case of a square crack of side L, where      0 1      ,   2 Q L y and  0    0 0.5 y L . Then                                                    0.4 0.4 0 0 0.4 0.4 0 0 1 3 1.033·tanh 4.34      · tanh 4.34      · 2 2 Q’     3 1              tanh 4.34      ·tanh 4.34      2 2 y y L L y y L L (9) In the r.h.s. of the Eqn. (8), the only significant contribution is given by             0.4 01Q’    1.033 tanh 4.34      2 y L (10) P. Livieri et alii, Frattura ed Integrità Strutturale, 63 (2023) 71-79; DOI: 10.3221/IGF-ESIS.63.07 76 Eqns. (8) and (10) were recently announced in a very similar preliminary form in paper [20]. We remark that Eqn. (6) can be reasonably extended to every convex set with corners. In the case of a square-like flaw, the Oore-Burns integral can be analytically expressed in simplified form in the middle of the side and in the middle of the rounded corner [24]. The Oore-Burns integral will be approximated by means of Riemann sums plus a suitable asymptotic correction in terms of mesh size [25, 26]. ASYMPTOTIC BEHAVIOUR OF THE SIF ON REGULAR POLYGONS et  be a regular polygon with N sides, inscribed in the disk of radius a. This means that the length of the side is       2· ·sinL a N , k=2π/N. Then (see Fig. 5) from Eqn. (6) after some simple steps:  Q’  tanh 1.74 1 𝜆 𝑁 . (11) where          0    ,  Q cos a y N , 0     · 2 L y ,   0 1 . Taking into account that on the unitary disk, the O-Integral takes the value   2 n , it follows the asymptotic behavior of the SIF, for N: K  √ √ tanh 1.74 1 𝜆 𝑁 . (12) Figure 5: Polygonal crack. Fig. 6 shows the graph of the SIF in dimensionless form for N equal to 20. In the middle of the side the SIF tends to have a constant value near to the one of a circular crack. At the corner the value of the SIF is null and in the proximity of the corner the trend shows a cup as calculated in reference [20]. In order to check Eq. (12), an accurate numerical FE analysis has been proposed in a case of a regular hexagonal crack. Figs. 7 and 8 show a three-dimensional model and the FE model, respectively. The mesh is refined only near the point of interest where the SIF is evaluated as proposed in previous works [22–25]. The dimensions of smaller elements at the tip of the crack were in the order of 10-5 mm. Finally, Fig. 9 proposes the comparison of the SIF in dimensionless form evaluated along the border in the range [0, L/2]. Eq. (12) appears precise and modifies the value of the SIF only near the corner. The value of KI,OB was evaluated by means of the procedure valid for a star domain proposed in reference [26] where a mesh with a step of 0.0157 radiant in the tangential direction and 80 terms in the Fourier series have been considered. L P. Livieri et alii, Frattura ed Integrità Strutturale, 63 (2023) 71-79; DOI: 10.3221/IGF-ESIS.63.07 77 Figure 6: Stress intensity factor in dimensionless form for a regular polygonal crack with N=20 (L/a=0.313) under uniform tensile load. Figure 7: Three-dimensional model for a hexagonal crack (L=1 mm, h/L= 200, d/L=20). Figure 8: Mesh and boundary conditions for a hexagonal crack of Fig. 7. P. Livieri et alii, Frattura ed Integrità Strutturale, 63 (2023) 71-79; DOI: 10.3221/IGF-ESIS.63.07 78 Figure 9: Stress intensity factor in dimensionless form for a hexagonal crack under uniform tensile load. CONCLUSIONS n this paper an accurate correction factor that multiplies the Oore-Burns stress intensity factor (SIF) is given for embedded polygonal cracks. This correction factor, ranging from zero up to one, essentially depends on the distance from the nearest corner, and it is of fundamental importance near the corner where the value of the Oore-Burns SIF does not reach a null value. Far from the corner, the correction factor quickly increases up to one. Furthermore, an accurate formula for regular polygonal cracks is given. A comparison with the FE results for hexagonal cracks shows satisfactory results. NOMENCLATURE α opening angle a crack size, disk radius Y shape factor  crack shape  crack border Q point of  'Q point of crack border KI mode I stress intensity factor KI,OB mode I stress intensity factor calculated with the Oore-Burns’ equation N number of corners  complementary angle s arc length σn nominal tensile stress in the ,x y Cartesian coordinate system p perimeter of  c perimeter of the smallest disk containing   correction factor ,x y actual Cartesian coordinate system I P. Livieri et alii, Frattura ed Integrità Strutturale, 63 (2023) 71-79; DOI: 10.3221/IGF-ESIS.63.07 79 REFERENCES [1] Bueckner, H.F., (1970). A novel principle for the computation of stress intensity factors, ZAMM 50, pp. 529–546. [2] Rice, J.R. (1972). Some remarks on elastic crack tip stress fields. Int J Solids Struct, 8, pp. 751–8. [3] Rice, J.R., (1989). Weight function theory for three-dimensional elastic crack analysis. ASTM STP1020, Wei R.P. and Gangloff R.P., Eds. Philadelphia, American Society for Testing and Materials, pp. 29–57. [4] Oore, M., Burns, D.J., (1980). Estimation of stress intensity factors for embedded irregular cracks subjected to arbitrary normal stress fields. Journal of Pressure Vessel Technology ASME, 102, pp. 202–211. [5] BS 7910:2019 Guide to methods for assessing the acceptability of flaws in metallic structures. [6] Carpinteri, A., Brighenti, R., Spagnoli, A. (2000). Fatigue growth simulation of part-through flaws in thick-walled pipes under rotary bending. Int J Fatigue 22(1), pp. 1–9. [7] Wang, X., Lambert, S.B. (1995). Stress intensity factors for low aspect ratio semi-elliptical surface cracks in finite- thickness plates subjected to nonuniform stresses. Engng Fract Mech, 51, pp. 517–32. [8] Chai, G., Zhang, K., Wu, D. (1996). Analyses on interactions of two identical semielliptical surface cracks in the internal surface of a cylindrical pressure vessel, Int. J. Press. Vessels Pip. 67 (2), pp. 203–210. [9] Glinka, G., Shen, G. (1991). Universal features of weight functions for cracks in mode I. Engng Fract Mech. 40, pp. 1135–46. [10] Shen, G., Glinka, G. (1991). Determination of weight function from reference stress intensity factor. Theor Appl Fract Mech; 15, pp. 237–45. [11] Gorner, F., Mattheck, C., Morawietz, P., Munz, D. (1985). Limitation of the Petrosky–Achenbach crack opening displacement approximation for the calculation of weight function. Engng Fract Mech. 22, pp. 269–77. [12] Fett, T. (1988). Limitation of the Petrosky–Achenbach procedure demonstrated for a simple load case. Engng Fract Mech; 29, pp. 713–6. [13] Petrosky, H.J., Achenbach, J.D. (1971). Computation of the weight function from a stress intensity factor. Engng Fract Mech. 10, pp. 257–66. [14] Zerbs, U. T., Schödel, M., Webster, S., Ainsworth, R. (2007). Fitness-for-Service Fracture Assessment of Structures Containing Cracks: A Workbook based on the European SINTAP/FITNET procedure, Elsevier, 1st ed. Oxford, Amsterdam, The Netherlands. [15] Murakami, Y. (Chief ed) (2001), Stress Intensity Factors Handbook, 4, 5, Pergamon. Press, Oxford, UK. [16] Fett, T., Munz, D., (1997). Stress intensity factors and weight functions, Computational Mechanics Publications. [17] Tada, H., Paris, C.P., Irwin, G.R., (2000). The stress analysis of cracks handbook. Third edition, ASME press. [18] Murakami, Y, (2002). Metal Fatigue: Effects of small defects and non-metallic inclusions, Elsevier. [19] Murakami, Y. and Endo, M. (1983) Quantitative evaluation of fatigue strength of metals containing various small defects or cracks. Eng. Fract. Mech., 17, pp. 1–15. [20] Livieri, P., Segala, S. (2021) Asymptotic behaviour of the Oore-Burns integral for cracks with a corner and correction formulae for embedded convex defects, Engineering Fracture Mechanics, 252, 107663, DOI: 10.1016/j.engfracmech.2021.107663. [21] Ascenzi, O., Pareschi, L., Segala, F. (2002). A precise computation of stress intensity factor on the front of a convex planar crack. International Journal for numerical methods in Engineering 54, pp. 241–261. [22] Livieri, P., Segala F. (2015). New weight functions and second order approximation of the Oore-Burns integral for elliptical cracks subject to arbitrary normal stress field, Eng. Fract. Mech. 138, pp. 100–117. [23] Livieri, P., Segala, S. (2016). Stress intensity factors for embedded elliptical cracks in cylindrical and spherical vessels Theoretical and Applied Fracture Mechanics 86(1), pp. 260–266. [24] Livieri P., Segala S. (2020). A closed form for the Stress Intensity Factor of a small embedded square-like flaw, Frattura ed Integrità Strutturale, 14 (54), pp. 182-191, DOI: 10.3221/IGF-ESIS.54.13. [25] Livieri, P., Segala, F., (2014). Sharp evaluation of the Oore-Burns integral for cracks subjected to arbitrary normal stress field, Fatigue & Fracture of Engineering Materials & Structures 37, pp. 95–106. [26] Livieri, P., Segala, F., (2018). An approximation in closed form for the integral of Oore–Burns for cracks similar to a star domain, Fatigue & Fracture of Engineering Materials & Structures 41, pp. 3–19. << /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