Microsoft Word - numero_42_art_21.docx M. Peron et alii, Frattura ed Integrità Strutturale, 42 (2017) 196-204; DOI: 10.3221/IGF-ESIS.42.21 196 Notch stress intensity factors under mixed mode loadings: an overview of recent advanced methods for rapid calculation M. Peron, S.M.J. Razavi, F. Berto, J. Torgersen Department of Mechanical and Industrial Engineering, Norwegian University of Science and Technology (NTNU), Richard Birkelands vei 2b, 7491, Trondheim, Norway. mirco.peron@ntnu.no, javad.razavi@ntnu.no, filippo.berto@ntnu.no, jan.torgersen@ntnu.no ABSTRACT. Recently some methods for the rapid calculation of notch stress intensity factors (NSIFs) have been developed and three of them are compared in this work. First, the criteria proposed by Lazzarin et al. and Treifi et al. have been reviewed. The former is based on the calculation of the mean value of SED on two different control volume (characterized by two different radius values) centred at the stress singularity point, whereas the latter takes advantage of the strain energy density averaged within two control volumes (semi-circular sector) centred at the notch tip. Then, a new method based on the evaluation of the total and deviatoric SED averaged in a single control volume has been proposed. Finally, plate specimens weakened by different notch geometries have been subjected to the application of the above mentioned methods and the obtained values of the NSIFs have been compared with those derived according to Gross and Mendelson. KEYWORDS. NSIFs; SED; Control volume; Mixed mode; FE analysis. Citation: Peron, M., Razavi, S.M.J., Berto, F., Torgersen, J., Notch stress intensity factors under mixed mode loadings: an overview of recent advanced methods for rapid calculation, Frattura ed Integrità Strutturale, 42 (2017) 196-204. Received: 02.07.2017 Accepted: 29.07.2017 Published: 01.10.2017 Copyright: © 2017 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 ecause of the detrimental effects played by notches on mechanical behavior [1-3] , several criteria have been proposed in literature for the assessment of material strength of components [4,5]. Concerning brittle and quasi- brittle materials, Notch stress intensity factor (NSIF) has been reported to be an effective criterion for the evaluation of static strength [6], and also for high cycle fatigue strength of structural materials weakened by sharp V- notches [7] and of welded joints [8,9]. Gross and Mendelson [10] expressed mode I and mode II NSIFs, which define the intensity of the asymptotic stress field ahead of the sharp V-notch tip. (1) (2)    111 002 limrK r               212 002 lim rrK r            B M. Peron et alii, Frattura ed Integrità Strutturale, 42 (2017) 196-204; DOI: 10.3221/IGF-ESIS.42.21 197 where r and θ is the radial and angular coordinate of a polar coordinate system centred at the notch tip (Fig. 1), σθθ and τrθ are the stress components according to the coordinate system and the mode I and mode II first eigenvalues in William’s equations [11] are denoted respectively as λ1 and λ2. The main drawback exploiting NSIF criterion for the calculation of (1) and (2) is that very refined meshes are needed using finite element (FE) analysis, and this requirement holds also for cracked component [12-17]. Figure 1: Polar coordinate system centred at the notch tip. However, an approach based on the determination of the total elastic strain energy density (SED) averaged over a control volume of radius R0 surrounding the points of stress singularity has been developed [18] and it has broadly gained interest in the last years, since it allows to exploit course meshes overcoming the main drawback of NSIF criterion both for static static [19-24] and fatigue [1, 25-27] strength evaluation. Moreover, some approximate methods for the rapid calculation of the NSIFs, based on the averaged strain energy density, are available in literature and, referring to mixed mode loading, they required the solution of a system of two equations in two unknowns (K1 and K2). In this work, two of these methods, one proposed by Lazzarin et al. [28] and the other by Treifi and Oyadiji [29], have been reviewed, and, afterwards, a new method based on the evaluation of the total and deviatoric SED averaged in a single control volume has been proposed, allowing also in this case the determination of the SIFs, KI and KII, of cracks under mixed mode loading by means of two independent equations, one linked to the total SED and the other to the deviatoric one. APPROXIMATE METHODS Lazzarin et al. approach he first method has been proposed by Lazzarin et al. [28] and it is based on the evaluation of the averaged SED on two different control volumes (circular sectors) centred at the notch tip and characterized by the radii Ra and Rb (Fig. 2a). Known the SED values ( aW and bW ), by means of a FE analysis, and defined the control radii (Ra and Rb), it is possible to obtain a system of two equations in two unknowns (K1 and K2):         1 2 1 2 2 2 2 21 1 2 2 , 1 22 1 2 1 1 2 2 2 2 21 1 2 2 , 1 22 1 2 1 1 2 1 2 2 2 1 2 2 2 a FE a a a a b FE b b b b I K I K W C K D K E R R I K I K W C K D K E R R                                          (3) where is the Young’s modulus of the material while I1 and I2 are the integrals of the angular stress functions [18], which depend on the notch opening angle, 2α = 2π - 2γ, and the Poisson's ratio ν. This method cannot be applied to a crack subjected to mixed mode loading, since an indeterminate system of equations would be obtained. Solving the system of equations, the values of the NSIFs can be determined: T M. Peron et alii, Frattura ed Integrità Strutturale, 42 (2017) 196-204; DOI: 10.3221/IGF-ESIS.42.21 198 (4) (5) Figure 2: Control volumes in the Lazzarin et al. approach (a) and in the Treifi et al. approach (b). Treifi et al. approach The second method has been proposed by Treifi et al. [29] and it is based on the evaluation of the averaged SED on two different control volumes (semi-circular sectors with a central angle equal to γ) centred at the notch tip and characterized by a radius R (Fig. 2b). Known the SED values ( aW and bW ) by means of a FE analysis, and defined the control radius (R), it is possible to obtain a system of two equations in two unknowns (K1 and K2):         1 21 2 1 21 2 2 2 1, 2, 12, 2 21 2 1 2 , 1 2 1 222 1 2 1 1 2 1 2 2 2 1, 2, 12, 2 21 2 1 2 , 1 2 1 222 1 2 1 1 2 1 2 21 2 ( ) 21 2 ( ) s s s a FE s s s b FE I I IK K K K W MK NK QK K E RR R I I IK K K K W MK NK QK K E RR R                                                      (6) where I1,s, I2,s and I12,s are the integrals of the angular stress functions [29], which depend on the notch opening angle, 2α, the angle defined by the semi-circular sector, γ, and the Poisson's ratio ν. In this case also the contribution of the mixed term (K1.K2 is present because the integration for the strain energy evaluation is not performed on a control volume symmetrical with respect to the notch bisector line (Fig. 2b). Due to the presence of the mixed term it is possible to decouple the contributions of the loading modes, obtaining a solution of the system even in the crack case. Solving the system of equations, the values of the NSIFs can be determined:       2 2 , , , , , , 2 1 4 4 a FE b FE a FE b FE a FE b FE MN W W W W W W Q K M       (7) , , 2 12 a FE b FEW WK QK   (8) New approach based on the deviatoric SED In the present contribution a new method is proposed. It is based on the evaluation of the total and deviatoric SED averaged in a single control volume (circular sector) centred at the notch tip and characterized by a radius R (Fig. 3a). Two independent equations can be obtained: one linked to the total SED and the other to the deviatoric one. In this way it is possible to obtain a solution of the system even in the case of a crack. , , 1 a b FE b a FE a b b a D W D W K D C D C    2 , 1 2 a FE a a W C K K D   M. Peron et alii, Frattura ed Integrità Strutturale, 42 (2017) 196-204; DOI: 10.3221/IGF-ESIS.42.21 199 As previously, knowing the SED values ( and ), by means of a FE analysis, and defining the control radius (R), it is possible to obtain a system of two equations in two unknowns (K1 and K2): (9) where I1,dev and I2,dev are the integrals of the angular stress functions related to the deviatoric strain energy density [18], which depend on the notch opening angle, 2α, and the Poisson's ratio ν. Solving the system of equations, the values of the NSIFs can be determined: (10) (11) Figure 3: Control volumes in the new approach (a) and in the modified version of the new approach (b). As discussed earlier, the total SED can be derived directly from nodal displacements, so that also coarse meshes are able to give sufficiently accurate values for it. On the other hand, the calculation of the deviatoric SED by means of a FE code is based on the von Mises equivalent stress averaged within the element [18]. This quantity is more sensitive to the refinement level of the adopted mesh, so the new proposed method could not be mesh-insensitive. With the aim to improve the results obtained from the application of the new method (based on the deviatoric SED) in the case of coarse meshes, a modified version is proposed. The approach is similar to the previous but it is applied to a control volume consisting of a circular ring (Fig. 3b). Being the calculation of the deviatoric SED by means of a FE code based on the von Mises equivalent stress averaged within the element, that is a parameter sensitive to the refinement level of the adopted mesh, it could be useful to exclude from the calculation the area characterized by the highest stress gradient (i.e. the region close to the notch tip) in the case of coarse meshes. The control volume results to be constituted by a circular ring characterized by an outer radius Ra and by an inner radius Rb (Fig. 3b). As before, knowing the SED values (W and devW ), by means of a FE analysis, and defining the control radii (Ra and Rb), it is possible to obtain a system of two equations in two unknowns (K1 and K2):     1 1 2 2 1 1 2 2 2 2 1 1 2 2 2 2 2 22 2 1 2 2 2 1, 1 2, 2 , 2 2 2 22 2 1 2 1 1 1 1 1 2 22 1 1 1 1 1 2 23 FE a b a ba b dev dev dev FE a b a ba b I K I K W R R R RE R R I K I K W R R R RE R R                                                                   (12) Solving this system of equations, as already shown in the previous cases, the values of the NSIFs can be determined. W devW         1 2 1 2 2 2 2 21 1 2 2 1 22 1 2 1 1 2 2 2 1, 2, 2 21 2 , 1 22 1 2 1 1 2 1 2 2 2 1 3 2 2 FE dev dev dev FE dev dev I K I K W SK TK E R R I IK K W S K T K E R R                                           , 1 dev FE dev FE dev dev TW T W K TS T S    2 , 1 2 dev FE dev dev W S K K T   M. Peron et alii, Frattura ed Integrità Strutturale, 42 (2017) 196-204; DOI: 10.3221/IGF-ESIS.42.21 200 RESULTS he methods described above have been applied to five different geometries of notched plates subjected to mixed mode I+II loading. For each geometry K1 and K2 have been first calculated according to Gross and Mendelson (Eqs. 1 and 2), by means of FE analyses adopting very refined meshes in the close neighbourhood of the notch tip (size of the smallest element of the order of 10-5 mm). Afterwards the approximate methods have been applied, taking into consideration three different values of the control radius R0 (0.1, 0.01 and 0.001 mm) and by using coarse and refined FE meshes. The mesh has been performed by means of 8-nodes elements (PLANE 183), using the FE code ANSYS® 14.5. In the FE analyses a Poisson’s ratio ν equal to 0.3 and a Young’s modulus E equal to 206 GPa have been adopted. In the summary tables the normalized NSIFs are reported, according to the definition: , 1 i i i normalized K K a     (13) Case studies: geometries and results The geometries taken into consideration (Fig. 4) consist in notched finite or infinite plates subjected to mixed mode I+II loading. The geometry of the finite plates is characterized by equal width and height, 2W = H = 10 mm, while plates of infinite extension are characterized by 2W = H = 10 mm. Diamond shape notches (Fig. 4a) are characterized by a projected notch depth 2a = 2 mm and a notch inclination angle φ = 45o. Two different notch opening angle 2 have been analysed: 45° and 30°. The obtained results and the comparison between the different approaches are reported in Table 1 and Table 2. The square hole (Fig. 4b) is characterized by a projected notch depth 2a = 2 mm and a notch opening angle 2α = 90o. The obtained results and the comparison between the different approaches are reported in Table 3. Figure 4: Case studied for the application of the method: diamond shape notch (a), square hole (b), central tilted crack in a finite plate (c) and central tilted crack in a plate of infinite extension (d). DISCUSSION abs. 1-5 show the results obtained from different geometries of notched plates subjected to mixed mode I+II loading, by adopting coarse and refined meshes for the application of the approximate methods. All the considered approaches allow to obtain very good approximations when refined meshes are adopted in the FE analyses, being the deviations, with respect to the values calculated according to Gross and Mendelson [10], lower than 1% in most of the cases. The obtained results become more interesting when coarse meshes, which require shorter calculation time, are adopted in the FE analyses. It can be observed that the approach of Lazzarin et al. [28] enables to obtain very good approximations, being the deviations lower than 1% in most of the cases analyzed in the present contribution. It should be noted that this method has not been applied to cracks subjected to mixed mode loading (Tables 4, 5), since an indeterminate system of equations would be obtained. The percentage error increases to about 3-6% in the case of Treifi et al. approach [29], which reaches a maximum percentage deviation of 12-18% in the case of tilted cracks (Table 5). T T M. Peron et alii, Frattura ed Integrità Strutturale, 42 (2017) 196-204; DOI: 10.3221/IGF-ESIS.42.21 201 The new proposed method, based on the evaluation of the total and deviatoric SED, allows to obtain a percentage error close to that observed in the case of Treifi et al. [29]. The deviation still remains greater than that observed in the case of Lazzarin et al. [28], because of the dependence of the deviatoric SED on the mesh size. This problem is overcome by the modified version of the new method that, through a control volume consisting of a circular ring, enables to exclude the region characterized by the highest stress gradient making the method less sensitive to the refinement level of the adopted mesh. The new method, particularly the modified version, provides very good approximations and a greater applicability than the approach of Lazzarin et al. [28] and [30], so it could be useful for rapid calculation of the NSIFs. Coarse mesh (64 finite elements) Refined mesh (3128 finite elements) R0 [mm] Method K1 K2 ΔK1 (%) ΔK2 (%) K1 K2 ΔK1 (%) ΔK2 (%) Gross and Mendelson 0.656 0.911 0.1 and 0.075 Lazzarin et al. 0.650 0.919 -0.91 0.88 0.650 0.919 -0.91 0.88 0.1 Treifi et al. 0.694 0.878 5.79 -3.62 0.693 0.878 5.64 -3.62 0.1 New method 0.681 0.891 3.81 -2.20 0.666 0.905 1.52 -0.66 0.1 New modified method 0.664 0.908 1.22 -0.33 0.01 and 0.0075 Lazzarin et al. 0.657 0.909 0.15 -0.22 0.655 0.912 -0.15 0.11 0.01 Treifi et al. 0.665 0.895 1.37 -1.76 0.663 0.897 1.07 -1.54 0.01 New method 0.671 0.883 2.29 -3.07 0.657 0.909 0.15 -0.22 0.01 New modified method 0.656 0.911 0.00 0.00 0.001 and 0.00075 Lazzarin et al. 0.659 0.901 0.46 -1.10 0.656 0.909 0.00 -0.22 0.001 Treifi et al. 0.658 0.907 0.31 -0.44 0.657 0.908 0.15 -0.33 0.001 New method 0.671 0.856 2.29 -6.04 0.658 0.901 0.30 -1.10 0.001 New modified method 0.658 0.905 0.30 -0.66 Table 1: Comparison between approximate methods for NSIFs evaluation of diamond-shaped notch (2α = 45°). Coarse mesh (64 finite elements) Refined mesh (3063 finite elements) R0 [mm] Method K1 K2 ΔK1 (%) ΔK2 (%) K1 K2 ΔK1 (%) ΔK2 (%) Gross and Mendelson 0.654 0.813 0.1 and 0.075 Lazzarin et al. 0.648 0.818 -0.92 0.62 0.650 0.817 -0.61 0.49 0.1 Treifi et al. 0.670 0.803 2.45 -1.23 0.681 0.796 4.13 -2.09 0.1 New method 0.686 0.792 4.89 -2.58 0.663 0.808 1.38 -0.62 0.1 New modified method 0.660 0.810 0.92 -0.37 0.01 and 0.0075 Lazzarin et al. 0.657 0.811 0.46 -0.25 0.655 0.813 0.15 0.00 0.01 Treifi et al. 0.621 0.847 -5.05 4.18 0.638 0.829 -2.45 1.97 0.01 New method 0.677 0.789 3.52 -2.95 0.657 0.811 0.46 -0.25 0.01 New modified method 0.656 0.812 0.31 -0.12 0.001 and 0.00075 Lazzarin et al. 0.660 0.807 0.92 -0.74 0.656 0.811 0.31 -0.25 0.001 Treifi et al. 0.673 0.784 2.91 -3.57 0.660 0.805 0.92 -0.98 0.001 New method 0.675 0.779 3.21 -4.18 0.659 0.807 0.76 -0.74 0.001 New modified method 0.658 0.809 0.31 -0.49 Table 2: Comparison between approximate methods for NSIFs evaluation of diamond-shaped notch (2α = 30°). M. Peron et alii, Frattura ed Integrità Strutturale, 42 (2017) 196-204; DOI: 10.3221/IGF-ESIS.42.21 202 Coarse mesh (48 finite elements) Refined mesh (2206 finite elements) R0 [mm] Method K1 K2 ΔK1 (%) ΔK2 (%) K1 K2 ΔK1 (%) ΔK2 (%) Gross and Mendelson 0.618 1.209 0.1 and 0.075 Lazzarin et al. 0.613 1.229 -0.81 1.65 0.613 1.229 -0.81 1.65 0.1 Treifi et al. 0.604 1.249 -2.27 3.31 0.632 1.184 2.27 -2.07 0.1 New method 0.635 1.175 2.75 -2.81 0.625 1.200 1.13 -0.74 0.1 New modified method 0.629 1.196 1.78 -1.08 0.01 and 0.0075 Lazzarin et al. 0.617 1.205 -0.16 -0.33 0.617 1.210 -0.16 0.08 0.01 Treifi et al. 0.617 1.211 -0.16 0.17 0.618 1.200 0.00 -0.74 0.01 New method 0.601 1.394 -2.75 15.30 0.617 1.212 -0.16 0.25 0.01 New modified method 0.618 1.192 0.00 -1.41 0.001 and 0.00075 Lazzarin et al. 0.618 1.166 0.00 -3.56 0.618 1.182 0.00 -2.23 0.001 Treifi et al. 0.618 1.213 0.00 0.33 0.617 1.202 -0.16 -0.58 0.001 New method 0.627 1.167 1.46 -3.54 0.618 1.178 0.00 -2.56 0.001 New modified method 0.619 1.130 0.16 -6.53 Table 3: Comparison between approximate methods for NSIFs evaluation of diamond-shaped notch (2α = 90°). Coarse mesh (64 finite elements) Refined mesh (3395 finite elements) R0 [mm] Method K1 K2 ΔK1 (%) ΔK2 (%) K1 K2 ΔK1 (%) ΔK2 (%) Gross and Mendelson 0.655 0.638 0.1 Treifi et al. 0.636 0.642 -2.90 0.63 0.660 0.637 0.76 -0.16 0.1 New method 0.697 0.620 6.41 -2.82 0.639 0.645 -2.44 1.10 0.1 New modified method 0.639 0.645 -2.44 1.10 0.01 Treifi et al. 0.613 0.654 -6.41 2.51 0.635 0.647 -3.05 1.41 0.01 New method 0.708 0.616 8.09 -3.45 0.653 0.640 -0.31 0.31 0.01 New modified method 0.653 0.640 -0.31 0.31 0.001 Treifi et al. 0.624 0.651 -4.73 2.04 0.644 0.644 -1.68 0.94 0.001 New method 0.712 0.615 8.70 -3.61 0.662 0.636 1.07 -0.31 0.001 New modified method 0.657 0.639 0.31 0.16 Table 4: Comparison between approximate methods for NSIFs evaluation of central tilted crack (2α = 0°) in a plate of finite extension. CONCLUSIONS n the present contribution three methods for the rapid calculation of the NSIFs, based on the averaged strain energy density, are compared. The first method, proposed by Lazzarin et al., is based on the calculation of the SED averaged in two different control volumes centred at the notch tip. This approach cannot be applied to notch with zero opening angle (cracks) subjected to mixed mode loading. The second method, presented by Treifi et al., overcomes this problem taking advantage of the strain energy density averaged within two control volumes (semi-circular sector) centred at the notch tip. Then a new method based on the evaluation of the total and deviatoric strain energy density has been proposed. I M. Peron et alii, Frattura ed Integrità Strutturale, 42 (2017) 196-204; DOI: 10.3221/IGF-ESIS.42.21 203 Coarse mesh (64 finite elements) Refined mesh (3395 finite elements) R0 [mm] Method K1 K2 ΔK1 (%) ΔK2 (%) K1 K2 ΔK1 (%) ΔK2 (%) Gross and Mendelson 0.595 0.595 0.1 Treifi et al. 0.667 0.564 12.10 -5.21 0.703 0.547 18.15 -8.07 0.1 New method 0.649 0.572 9.08 -3.87 0.596 0.595 0.17 0.00 0.1 New modified method 0.598 0.594 0.50 -0.17 0.01 Treifi et al. 0.582 0.599 -2.18 0.67 0.603 0.592 1.34 -0.50 0.01 New method 0.649 0.571 9.08 -4.03 0.597 0.594 0.34 -0.17 0.01 New modified method 0.598 0.594 0.50 -0.17 0.001 Treifi et al. 0.575 0.602 -3.36 1.18 0.593 0.595 -0.34 0.00 0.001 New method 0.649 0.571 9.08 -4.03 0.612 0.588 2.86 -1.18 0.001 New modified method 0.598 0.594 0.50 -0.17 Table 5: Comparison between approximate methods for NSIFs evaluation of central tilted crack (2α = 0°) in a plate of infinite extension. The described methods have been applied to plates subjected to mixed mode I+II loading and weakened by different V- notch geometries. The values of the NSIFs derived according to Gross and Mendelson have been compared with those obtained by means of the approximate methods taking into consideration three different values of the control radius R0 (0.1, 0.01 and 0.001 mm) and by using coarse and refined FE meshes. The comparison of the results shown that the new proposed method provides the best combination between the degree of approximation and the level of applicability, so it could be useful for rapid calculation of the NSIFs. Central tilted cracks in a plate of finite (Fig. 4c) or infinite (Fig. 4d) extension are characterized by a projected crack length 2a = 2 mm and a crack inclination angle φ = 45o. The obtained results and the comparison between the different approaches are reported in Tables 4, 5. It is worth noting that the case of the infinite plate has been modelled as a finite plate characterized by width and height two orders of magnitude greater than the crack length. REFERENCES [1] Lazzarin, P., Comportamento a fatica dei giunti saldati in funzione della densità di energia di deformazione locale: influenza dei campi di tensione singolari e non singolari, Frattura ed Integrita Strutturale, 9 (2009) 13-26. [2] Maragoni, L., Carraro, P. A., Peron, M., Quaresimin, M., Fatigue behaviour of glass/epoxy laminates in the presence of voids, Int. J. Fatigue, 95 (2017) 18–28. [3] Brotzu, A., Felli, F., Pilone, D., Effects of the manufacturing process on fracture behaviour of cast TiAl intermetallic alloys, Frattura ed Integrita Strutturale, 27 (2013) 66-73. [4] Askes, H., Susmel, L., Gradient enriched linear-elastic crack tip stresses to estimate the static strength of cracked engineering ceramics, Frattura ed Integrita Strutturale, 25 (2013) 87-93. [5] Susmel, L., On the overall accuracy of the Modified Wöhler Curve Method in estimating high-cycle multiaxial fatigue strength, Frattura ed Integrita Strutturale, 16 (2011) 5-17, [6] Seweryn, A., Brittle fracture criterion for structures with sharp notches, Eng. Fract. Mech., 47 (1994) 673–681. [7] Boukharouba, T., Tamine, T., Niu, L., Chehimi, C., Pluvinage, G., The use of notch stress intensity factor as a fatigue crack initiation parameter, Eng. Fract. Mech., 52 (1995) 503–512. [8] Lazzarin, P., Tovo, R., A notch intensity factor approach to the stress analysis of welds, Fatigue Fract. Eng. Mater. Struct., 21 (1998) 1089–1103. [9] Atzori, B., Meneghetti, G., Fatigue strength of fillet welded structural steels: Finite elements, strain gauges and reality, Int. J. Fatigue, 23 (2001) 713–721. M. Peron et alii, Frattura ed Integrità Strutturale, 42 (2017) 196-204; DOI: 10.3221/IGF-ESIS.42.21 204 [10] Gross, B., Mendelson, A., Plane elastostatic analysis of V-notched plates, Int. J. Fract. Mech., 8 (1972) 267–276. [11] Williams, M. L., Stress singularities resulting from various boundary conditions in angular corners on plates in extension, J. Appl. Mech., 19 (1952) 526–528. [12] 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. [13] Ayatollahi, M.R., Razavi, S.M.J., Rashidi Moghaddam, M., Berto, F., Mode I fracture analysis of Polymethylmetacrylate using modified energy—based models, Phys. Mesomec., 18 (2015) 53-62. [14] Ayatollahi, M.R., Rashidi Moghaddam, M., Razavi, S.M.J., Berto, F., Geometry effects on fracture trajectory of PMMA samples under pure mode-I loading. Eng. Fract. Mech., 163 (2016) 449–461. [15] Ayatollahi, M.R., Razavi, S.M.J., Sommitsch, C., Moser, C., Fatigue life extension by crack repair using double stop- hole technique, Mater. Sci. Forum, 879 (2017) 3-8. [16] Razavi, S.M.J., Ayatollahi, M.R., Sommitsch, C., Moser, C., Retardation of fatigue crack growth in high strength steel S690 using a modified stop-hole technique, Eng. Fract. Mech., 169 (2017) 226–237. [17] Rashidi Moghaddam, M., Ayatollahi, M.R., Razavi, S.M.J., Berto, F., Mode II Brittle Fracture Assessment Using an Energy Based Criterion, Phys. Mesomec. (in press). [18] 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. [19] 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. [20] 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. [21] Torabi, A. R., Campagnolo, A., Berto, F., Local strain energy density to predict mode II brittle fracture in Brazilian disk specimens weakened by V-notches with end holes, Mater. Des., 69 (2015) 22–29. [22] 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. [23] Berto, F., Ayatollahi, M. R., Fracture assessment of Brazilian disc specimens weakened by blunt V-notches under mixed mode loading by means of local energy, Mater. Des., 32 (2011) 2858–2869. [24] Berto, F., Barati, E., Fracture assessment of U-notches under three point bending by means of local energy density, Mater. Des., 32 (2011) 822–830. [25] Livieri, P., Lazzarin, P., Fatigue strength of steel and aluminium welded joints based on generalised stress intensity factors and local strain energy values, Int. J. Fract., 133 (2005) 247–276. [26] Berto, F., Campagnolo, A., Lazzarin, P., Fatigue strength of severely notched specimens made of Ti-6Al-4V under multiaxial loading, Fatigue Fract. Eng. Mater. Struct., 38 (2015) 503–517. [27] Berto, F., Croccolo, D., Cuppini, R., Fatigue strength of a fork-pin equivalent coupling in terms of the local strain energy density, Mater. Des., 29 (2008) 1780–1792. [28] Lazzarin, P., Berto, F., Zappalorto, M., Rapid calculations of notch stress intensity factors based on averaged strain energy density from coarse meshes: Theoretical basis and applications, Int. J. Fatigue, 32 (2010) 1559–1567. [29] Treifi, M., Oyadiji, S. O., Strain energy approach to compute stress intensity factors for isotropic homogeneous and bi-material V-notches, Int. J. Solids Struct., 50 (2013) 2196–2212. [30] Gallo, P., Berto, F., Glinka, G Generalized approach to estimation of strains and stresses at blunt V-notches under non-localized creep, Fatigue Fract. Eng. Mater. Struct, 39 (2016) 292-306. << /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