Microsoft Word - numero_34_art_11 G Meneghetti et alii, Frattura ed Integrità Strutturale, 34 (2015) 109-115; DOI: 10.3221/IGF-ESIS.34.11 109 Focussed on Crack Paths A link between the peak stresses and the averaged strain energy density for cracks under mixed-mode (I+II) loading G. Meneghetti, B. Atzori University of Padova, Department of Industrial Engineering, Via Venezia 1, 35131, Padova (Italy) giovanni.meneghetti@unipd.it, bruno.atzori@unipd.it A. Campagnolo, F. Berto University of Padova, Department of Management and Engineering, Stradella S. Nicola 3, 36100, Vicenza (Italy) campagnolo@gest.unipd.it, berto@gest.unipd.it ABSTRACT. In this work, a link between the averaged strain energy density (SED) approach and the peak stress method in the case of cracks subjected to mixed mode (I+II) loading has been investigated. Some closed-form expressions of the strain energy density, averaged in a volume of radius R0, as function of the stress intensity factors are provided for plane strain conditions under mixed mode I+II loadings. On the basis of the peak stress method (PSM) some expressions useful to estimate the mode I and mode II stress intensity factors (SIFs) have been recently derived. These relationships take advantage of the elastic peak stresses from FE analyses carried out by using a given mesh pattern where the element size and type are kept constants. The evaluation of the SIFs from a numerical analysis of the local stress field usually requires very refined meshes and then large computational effort. The usefulness of the PSM-based expressions is that (i) only the elastic peak stresses numerically evaluated at the crack tip are needed and not a set of stress–distance data; (ii) the employed meshes are rather coarse if compared to those necessary for the evaluation of the whole local stress field. By substituting the PSM-based relationships in the closed-form expressions of the averaged SED it appears that the latter can be directly estimated by means of the elastic peak stresses evaluated at the crack tip. Several FE analyses have been carried out on cracked plates subjected to tension loading considering different geometrical combinations, varying the length 2a and the inclination ϕ of the crack (i.e. the mode mixity) as well as the size d of the adopted finite elements, with the aim to evaluate the local SED and the elastic peak stress components σpeak and τpeak. In all cases the numerical values of the SED derived from the FE analyses have been compared with those analytically obtained by using the expressions for the SED based on the elastic peak stresses, in order to verify the range of applicability of the proposed relationships. KEYWORDS. Local approaches; Stress intensity factors; Strain energy density; Finite elements; Coarse mesh. INTRODUCTION otch stress intensity factors (NSIFs) play an important role in static strength assessments of components made of brittle or quasi-brittle materials and weakened by sharp V-shaped notches [1]. This holds true also for components made of structural materials undergoing high cycle fatigue loading [2] as well as for welded joints N G Meneghetti et alii, Frattura ed Integrità Strutturale, 34 (2015) 109-115; DOI: 10.3221/IGF-ESIS.34.11 110 [3,4]. In plane problems, the mode I and mode II NSIFs for sharp V-notches, which quantify the intensity of the asymptotic stress distributions in the close neighbourhood of the notch tip, can be expressed by means of the Gross and Mendelson’s definitions [5]:   111 002 limrK r             (1)   212 002 lim rrK r             (2)  rr r    notch bisector r A=R0 2 R0 Figure 1: (a) Polar coordinate system centred at the notch tip. (b) Control volume (area) of radius R0 surrounding the V-notch tip. where (r,θ) is a polar coordinate system centred at the notch tip (Fig. 1a), σθθ and τrθ are the stress components according to the coordinate system and λ1 and λ2 are the mode I and mode II first eigenvalues in William’s equations [6], respectively. The condition θ = 0 characterizes all points of the notch bisector line. When the V-notch angle 2α is equal to zero, λ1 and λ2 equal 0.5 and K1 and K2 match the conventional stress intensity factors of a crack, KI and KII, according to the Linear Elastic Fracture Mechanics (LEFM). The main practical disadvantage in the application of the NSIF-based approach is that very refined meshes are needed to calculate the NSIFs by means of definitions (1) and (2). The modelling procedure becomes particularly time-consuming for components that cannot be analysed by means of two-dimensional models. Recently, Nisitani and Teranishi [7,8] presented a new numerical procedure suitable for estimating KI for a crack emanating from an ellipsoidal cavity. Such a procedure is based on the usefulness of the linear elastic stress σpeak calculated at the crack tip by means of FE analyses characterized by a mesh pattern having a constant element size. In particular Nisitani and Teranishi [7,8] were able to show that the ratio KI/σpeak depends only on the element size, so that the σpeak value can be used to rapidly estimate KI, provided that the adopted mesh pattern has been previously calibrated on geometries for which the exact value of KI is known. This approach has been theoretically justified and extended also to sharp V-shaped notches subject to mode I loading [9] giving rise to the so-called Peak Stress Method (PSM), which can be regarded as an approximate FE-based method to estimate the NSIFs. Later on, the PSM has been extended to cracks subjected to mode I as well as mode II stresses [10]. The element size required to evaluate K1 and K2 from σpeak and τpeak, respectively, is several orders of magnitude greater than that required to directly evaluate the local stress field. The second advantage of using σpeak and τpeak is that only a single stress value is sufficient to estimate K1 and K2, respectively, instead of a number of stress-distance FE data, as usually made by applying definitions (1) and (2). Since the units of the mode I and mode II NSIFs, K1 and K2, depend on the notch opening angle, a direct comparison of the NSIF values cannot be performed. This problem was overcome by Lazzarin and Zambardi [11], who proposed to use the total elastic strain energy density (SED) averaged over a sector of radius R0 (Fig. 1b) for static [11-14] and fatigue [11,15,16] strength assessments. With reference to plane strain conditions, the SED value can be evaluated as follows: 1 2 2 2 1 1 2 2 1 1 0 0 e K e K W E R E R                (3) where e1 and e2 [11] are two parameters which depend on the notch opening angle 2α and the Poisson’s ratio ν. In principle, Eq. (3) is valid when the influence of higher order, non-singular terms can be neglected inside the control volume. In the case of short cracks or thin welded lap joints, for example, the T-stress must be included in the local SED evaluation [17]. G Meneghetti et alii, Frattura ed Integrità Strutturale, 34 (2015) 109-115; DOI: 10.3221/IGF-ESIS.34.11 111 Aims of the present contribution are as follows:  to recall the fundamental concepts of the PSM for pure modes of loading;  to present the extension of the PSM to the case of mixed mode (I+II) loading;  to investigate a link between the SED approach and the PSM in the case of mixed mode (I+II) loading. THE PEAK STRESS METHOD FOR PURE MODES OF LOADING he Peak Stress Method (PSM) is a simplified numerical method to estimate the NSIFs parameters. Originally it was formulated for cases where only mode I singular stresses exists (i.e. K2 = 0 or mode II stresses are negligible). It has been based on a link between the exact value of mode I NSIF K1, see Eq. (1), and the linear elastic opening peak stress σpeak calculated at the V-notch tip according to the following expression [9]: 1 * 1 1 1.38FE peak K K d      (4) The PSM according to Eq. (4) was applied to correlate the fatigue strength of fillet-, full penetration and butt welded joints subjected to mode I loading [18-20]. Recently the Peak Stress Method has been extended also to mode II crack problems, linking the exact value of mode II NSIF K2, see Eq. (2) with 2 = 0° and 2 = 0.5, and the linear elastic sliding peak stress τpeak calculated at the crack tip according to the following expression [10]: 2 ** 2 1 3.38FE peak K K d      (5) In previous expressions d is the mean finite element size adopted when using the free mesh generation algorithm available in Ansys numerical code, while “exact NSIF values” must be meant as the values obtained using very refined FE mesh patterns in the numerical analyses and applying definitions (1) and (2) to the numerical results. Eqs. (4) and (5) are useful in practical applications because if the mean element size d is kept constant, then also K1/σpeak and K2/τpeak ratios are constant. Eqs. (4) and (5) are valid under the following conditions:  use of 4-node linear quadrilateral elements, as implemented in ANSYS® numerical code (PLANE 42 of Ansys element library or alternatively PLANE 182 with K-option 1 set to 3);  the pattern of finite elements around the V-notch tip must be that shown in Fig. 2b (see also [9, 10]); in particular, four elements share the node located at the crack tip;  concerning Eq. (4), V-notches characterised by an opening angle 2 ranging from 0° to 135°;  the ratio a/d must be greater than 3 in order to obtain * 1.38 3%FEK   , a being the semi-crack length (or the notch depth when dealing with open V-notches). When mode II (sliding) stresses are of interest, meshes are more refined such that the ratio a/d must be greater than 14 in order to obtain ** 3.38 3%FEK   . THE PEAK STRESS METHOD FOR MIXED MODE (I+II) LOADING n the present paragraph the Peak Stress Method is extended to mixed mode (I+II) crack problems. Consider a crack (2α = 0°) centred in a plate having the geometry reported in Fig. 2a and subjected to tensile loading. By varying the inclination angle ϕ of the crack it is possible to obtain different mode mixities, from pure mode I (ϕ = 0°) to mixed mode I+II (ϕ > 0°) loading. Different geometrical combinations have been considered, varying the projected crack length 2h (from 5 to 80 mm) and the inclination ϕ (from 0° to 60°) as well as the size of the element d, with the aim to investigate to what extent the PSM holds true. Finite element analyses have been performed by using the commercial code Ansys® and 4-node quadrilateral element (PLANE 42). The free mesh algorithm has been used in all numerical analyses and the sole control parameter set to generate the mesh has been the so-called ‘global element size’, i.e. the mean element size of the finite elements, which ranged from 0.5 mm to 10 mm. With the purpose of obtaining the pattern of finite T I G Meneghetti et alii, Frattura ed Integrità Strutturale, 34 (2015) 109-115; DOI: 10.3221/IGF-ESIS.34.11 112 elements oriented along the crack bisector line (see Fig. 2b), the geometry of the plate has been divided into six areas, such that each crack tip is shared by four areas, as shown in Fig 2a. By so doing four elements (each one belonging to a different area) share the node located at the crack tip. For the considered case, K1 = KI, K2 = KII, λ1 = λ2 = 0.5, while σpeak and τpeak represent the maximum elastic normal and tangential stress referred to the bisector line and evaluated at the crack tip according to Fig. 2a. The exact values of the mode I and mode II SIF, KI and KII, have been evaluated by means of further finite element analyses performed on the same geometries, but adopting very refined meshes (size of the smallest element of the order of 10-5 mm) in the close neighbourhood of the crack tip. Fig. 3 plots the results of the numerical analyses in terms of the non-dimensional parameters K*FE and K**FE defined in Eqs. (4) and (5). The results for the cases ϕ = 10°, 45°, 30° and 60° have been reported. From Fig. 3, K*FE and K**FE are seen to converge to the previously calibrated values, that is 1.38 [9] and 3.38 [10], respectively, within a scatter band of the numerical results of ±3% also in the case of mixed mode (I+II) loading. This occurs for a ratio a/d greater than a value between 3 and 6, for mode I loading, and between 12 and 21, for mode II loading. It can be observed that the minimum a/d ratios to assure the validity of PSM under mixed mode (I+II) loading are only slightly greater than the results obtained in [9], in the case of pure mode I (a/d  3), and the results reported in [10] with reference to pure mode II (a/d  14). Furthermore, as highlighted in [10], it should be noted that the mode II loading is more critical to analyse with the PSM than the mode I loading because it requires more refined finite element mesh patterns. Figure 2: (a) Geometry and loading condition of the analysed mixed mode crack problem. 2W = 200 mm. (b) Pattern of finite elements around the singularity point: four elements share the node located at the crack tip. A LINK BETWEEN THE PEAK STRESSES AND THE AVERAGED VALUE OF THE LOCAL STRAIN ENERGY DENSITY n the present paragraph, a link between the averaged SED [11,12] and the peak stresses [9,10] in the case of cracks subjected to mixed mode (I+II) loading is investigated. By substituting the PSM-based relationships, Eqs. (4) and (5), in the closed-form expression of the averaged SED, Eq. (3), it appears that the latter can be directly estimated by means of the elastic peak stresses evaluated at the crack tip, σpeak and τpeak: - -1 2 2 21 1 * **1 2 0 0 PSM FE peak FE peak e ed d W K K E R E R                                  (6) I G Meneghetti et alii, Frattura ed Integrità Strutturale, 34 (2015) 109-115; DOI: 10.3221/IGF-ESIS.34.11 113 Several refined FE analyses have been carried out on the same cracked plates taken into consideration in the previous paragraph, with the aim to evaluate the local SED averaged over a control volume centred at the crack tip. Different geometrical combinations have been considered, varying the length 2a and the inclination ϕ of the crack (i.e. the mode mixity), while the radius of the control volume R0 has been kept constant and equal to 0.1 mm. 0.00 0.50 1.00 1.50 2.00 2.50 3.00 3.50 4.00 1 10 100 a/d K * F E , K * * F E K*FE K**FE +3% -3% 1.38  6 2 = 0° = 10° 3.38 +3% -3%  12 (a) K*FE K**FE 0.00 0.50 1.00 1.50 2.00 2.50 3.00 3.50 4.00 1 10 100 a/d K * F E , K * * F E K*FE K**FE +3% -3% 1.38  4 2 = 0° = 30° 3.38 +3% -3%  14 K*FE K**FE (b) 0.00 0.50 1.00 1.50 2.00 2.50 3.00 3.50 4.00 1 10 100 a/d K * F E , K * * F E K*FE K**FE +3% -3% 1.38  5 2 = 0° = 45° 3.38 +3% -3%  21 (c) K*FE K**FE 0.00 0.50 1.00 1.50 2.00 2.50 3.00 3.50 4.00 1 10 100 a/d K * F E , K * * F E K*FE K**FE +3% -3% 1.38  3 2 = 0° = 60° 3.38 +3% -3%  16 K*FE K**FE (d) Figure 3: Calibration of the PSM approach for a crack (2α = 0°) under mixed mode (I+II) loading: (a) ϕ = 10°; (b) ϕ = 30°; (c) ϕ = 45°; (d) ϕ = 60°. Non-dimensional SIFs related to mode I and mode II. The mode mixity ratio (MM) has been evaluated according to the following definition: II I II K MM K K   (7) Eq. (7) provides as master cases MM = 0 for pure mode I with ϕ = 0°, MM = 0.5 for mixed mode with ϕ = 45° and MM = 1 for pure mode II loading. In all cases the numerical values of the SED calculated from the FE analyses have been compared with those analytically obtained by using the expressions for the SED based on the elastic peak stresses, Eq. (6), in order to verify the range of applicability of the proposed method. Being the exact values of the SIFs available, the mean value of the SED has been evaluated also according to Eq. (3). In particular the maximum difference between the SED parameter evaluated analytically (Eq. (3)) and numerically (by FEM) resulted about 5%, that means that the influence of higher order terms, as the T-stress, can be neglected in these cases, at least from an engineering point of view. The ratio between the SED based on the elastic peak stresses (Eq. 6, PSMW ) and the SED calculated from the FE analyses ( FEMW ) has been reported in Fig. 4, with reference to an inclination ϕ of the crack equal to 0°, 30°, 45° and 60°. From Fig. 4, it can be observed that the ratio /PSM FEMW W converges to unity, within a scatter band of ±10% for all different mode mixities taken into consideration. This occurs for a ratio a/d greater than 3 for the case MM = 0 (ϕ = 0°), 8.50 for G Meneghetti et alii, Frattura ed Integrità Strutturale, 34 (2015) 109-115; DOI: 10.3221/IGF-ESIS.34.11 114 MM = 0.37 (ϕ = 30°), 11 for MM = 0.50 (ϕ = 45°) and 16 for MM = 0.63 (ϕ = 60°). In particular the minimum a/d ratio to assure the validity of the proposed method increases as the mode II loading becomes dominant, that is as the mode mixity ratio (MM) defined by Eq. (7) increases. A similar behavior has been observed in the previous paragraph and in [10], with reference to K*FE and K**FE. 0.00 0.50 1.00 1.50 2.00 2.50 3.00 3.50 1 10 100 1000 a/d W P S M /W F E M   Rmm    a 0.00 0.50 1.00 1.50 2.00 2.50 3.00 3.50 1 10 100 1000 a/d W P S M /W F E M   Rmm    b 0.00 0.50 1.00 1.50 2.00 2.50 3.00 3.50 1 10 100 1000 a/d W P S M /W F E M 2 = 0° = 45° R0 = 0.1 mm   c  0.00 0.50 1.00 1.50 2.00 2.50 3.00 3.50 1 10 100 1000 a/d W P S M /W F E M   Rmm    d Figure 4: Ratio between the SED based on the elastic peak stresses ( PSMW ) and the SED derived from the FE analyses ( FEMW ) for a crack (2α = 0°) under mixed mode (I+II) loading: (a) ϕ = 0° (MM = 0), (b) ϕ = 30° (MM = 0.37), (c) ϕ = 45° (MM = 0.50) and (d) ϕ = 60° (MM = 0.63). CONCLUSIONS n the present contribution (*), a link between the averaged SED approach and the peak stress method in the case of cracks subjected to mixed mode (I+II) loading has been investigated:  On the basis of the peak stress method, some expressions useful to estimate the mode I and mode II SIFs, recently derived for pure modes of loading, have been verified also in the case of mixed mode (I+II) crack problems.  Since the normal and tangential peak stresses are proportional to the mode I and mode II SIFs, a link can immediately be established with the SED parameter by means of Eqs. (3), (4) and (5). By substituting the PSM-based relationships in the closed-form expressions of the averaged SED it appears that the latter can be directly estimated by means of the elastic peak stresses evaluated at the crack tip.  The ratio between the SED based on the elastic peak stresses and the SED derived from the FE analyses converge to a unit value for a ratio a/d greater than a value between 3 (MM = 0) and 16 (MM = 0.63). The minimum a/d ratio to assure the validity of the proposed method increases with increasing the mode mixity ratio (MM).  The usefulness of the SED expression based on the elastic peak stresses is that (i) only the elastic peak stresses numerically evaluated at the crack tip are needed and the definition of a control volume is no longer required; I G Meneghetti et alii, Frattura ed Integrità Strutturale, 34 (2015) 109-115; DOI: 10.3221/IGF-ESIS.34.11 115 (ii) the employed meshes are rather coarse if compared to those necessary for SED calculation, because the mean element size can be significantly greater than the radius of the control volume adopted for SED evaluation. (*) A version of the present contribution has already been presented at the Third IJFatigue & FFEMS Joint Workshop on "Characterisation of Crack Tip Fields" organized by the Italian Group of Fracture (IGF). REFERENCES [1] Seweryn, A., Brittle fracture criterion for structures with sharp notches, Eng. Fract. Mech., 47 (1994) 673–681. [2] Boukarouba, T., Tamine, T., Nui, L., Chemini, C. and Pluvinage, G., The use of notch stress intensity factor as a fatigue crack initiation parameter. Eng. Fract. Mech., 52 (1995) 503–512. [3] Lazzarin, P., Tovo, R., A notch intensity factor approach to the stress analysis of welds, Fatigue Fract. Eng. Mater. Struct., 21 (1998) 1089–1104. [4] Atzori, B., Meneghetti, G., Fatigue strength of fillet welded structural steels: finite elements, strain gauges and reality, Int. J. Fatigue, 23 (2001) 713–721. [5] Gross, R., Mendelson, A., Plane elastostatic analysis of V-notched plates, Int. J. Fract. Mech., 8 (1972) 267–272. [6] Williams, M. L., Stress singularities resulting from various boundary conditions in angular corners on plates in tension, J. Appl. Mech., 19 (1952) 526–528. [7] Nisitani, H., Teranishi, T., KI value of a circumferential crack emanating from an ellipsoidal cavity obtained by the crack tip stress method in FEM, In: Guagliano, M., Aliabadi, M. H., editors. Proceedings of the 2nd international conference on fracture and damage mechanics, (2001) 141–146. [8] Nisitani, H., Teranishi, T., KI value of a circumferential crack emanating from an ellipsoidal cavity obtained by the crack tip stress method in FEM, Engng. Fract. Mech., 71 (2004) 579–585. [9] Meneghetti, G., Lazzarin, P., Significance of the elastic peak stress evaluated by FE analyses at the point of singularity of sharp V-notched components, Fatigue Fract. Engng. Mater. Struct., 30 (2007) 95-106. [10] Meneghetti, G., The use of peak stresses for fatigue strength assessments of welded lap joints and cover plates with toe and root failures, Eng. Fract. Mech., 89 (2012) 40-51. [11] Lazzarin, P., Zambardi, R., A finite-volume-energy based approach to predict the static and fatigue behavior of components with sharp V-shaped notches, Int. J. Fract., 112 (2001) 275-298. [12] Berto, F., Lazzarin, P., Recent developments in brittle and quasi-brittle failure assessment of engineering materials by means of local approaches, Mater. Sci. Eng. R, 75 (2014) 1-48. [13] Lazzarin, P., Campagnolo, A., Berto, F., A comparison among some recent energy- and stress-based criteria for the fracture assessment of sharp V-notched components under Mode I loading, Theor. Appl. Fract. Mech., 71 (2014) 21- 30. [14] Torabi, A. R., Campagnolo, A., Berto, F., Experimental and theoretical investigation of brittle fracture in key-hole notches under mixed mode I/II loading, Acta Mech., 226 (2015) 2313-2322. [15] Livieri, P., Lazzarin, P., Fatigue strength of steel and aluminium welded joints based on generalized stress intensity factors and local strain energy values, Int. J. Fract., 133 (2005) 247-276. [16] Berto, F., Campagnolo, A., Lazzarin, P., Fatigue strength of severely notched specimens made of Ti-6Al-4Vunder multiaxial loading, Fatigue Fract. Engng. Mater. Struct., 38 (2015) 503-517. [17] Berto, F., Lazzarin, P., Multiparametric full-field representations of the in-plane stress fields ahead of cracked components under mixed mode loading, Int. J. Fatigue, 46 (2013) 16-26. [18] Meneghetti, G., The peak stress method applied to fatigue assessments of steel and aluminium fillet welded joints subjected to mode-I loading, Fatigue Fract. Engng. Mater. Struct. 31 (2008) 346–369. [19] Meneghetti, G., Atzori, B., Manara, G., The Peak Stress Method applied to fatigue assessments of steel tubular welded joints subject to mode-I loading, Engng. Fract. Mech., 77 (2010) 2100–2114. [20] Meneghetti, G., Campagnolo, A., Berto, F., Fatigue strength assessment of partial and full‐penetration steel and aluminium butt‐welded joints according to the peak stress method, Fatigue Fract. Engng. Mater. Struct. (in press) doi: 10.1111/ffe.12342.