Microsoft Word - numero_68_art_01_4781.docx Ch. F. Markides et alii, Frattura ed Integrità Strutturale, 68 (2024) 1-18; DOI: 10.3221/IGF-ESIS.68.01 1 Revisiting classical concepts of Linear Elastic Fracture Mechanics - Part II: Stretching finite strips weakened by single edge parabolically- shaped notches Christos F. Markides National Technical University of Athens, School of Applied Mathematical and Physical Sciences, Department of Mechanics, Zografou Campus, 5 Heroes of Polytechneion Avenue, 157 73, Attiki, Greece markidih@maill.ntua.gr, http://orcid.org/0000-0001-6547-3616 Stavros K. Kourkoulis National Technical University of Athens, School of Applied Mathematical and Physical Sciences, Department of Mechanics, Laboratory of Strength and Materials, Zografou Campus, 5 Heroes of Polytechneion Avenue, 157 73, Attiki, Greece stakkour@central.ntua.gr, http://orcid.org/0000-0003-3246-9308 KEYWORDS. Linear Elastic Fracture Mechanics, Parabolically-shaped notches, Single edge notched strip, Mode-I loading, Stress Concentration Factor, Stress Intensity Factor, Complex potentials. Citation: Kourkoulis, S.K., Markides, Ch.F, Revisiting classical concepts of Linear Elastic Fracture Mechanics-Part II: Mode I stress concentrations in single and double notched strips, Frattura ed Integrità Strutturale, 68 (2024) 1-18. Received: 14.12.2023 Accepted: 05.01.2024 Published: 11.01.2024 Issue: 04.2024 Copyright: © 2024 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 his is the second part of a short three-paper series, the aim of which is to revisit some classical concepts of Linear Elastic Fracture Mechanics (LEFM). In the first paper of the series [1] some controversial issues were discussed, related to: (i) the unnatural overlapping of the lips of ‘mathematical’ cracks in ‘infinite plates’ (a phenomenon that is predicted by the solution of the respective first fundamental problem of LEFM, in case it is formalistically applied), and, (ii) the closely associated issue of negative mode-I Stress Intensity Factors (SIFs). It could be counterargumented of course T https://youtu.be/gMy-3uo-x2A Ch. F. Markides et alii, Frattura ed Integrità Strutturale, 68 (2024) 1-18; DOI: 10.3221/IGF-ESIS.68.01 2 that concepts like ‘infinite plate’ and ‘mathematical’ cracks are of limited applicability in practical engineering problems (where neither the size of actual structural members is infinite nor the discontinuities are ‘mathematical’ cuts with singular point tips). In spite of this skepticism, it was highlighted in ref. [1] that the solutions of LEFM, if properly adjusted, can provide interesting results, also, for actual structures. In this context, an attempt is described here to relief the suffocating assumptions of ‘mathematical’ discontinuities and ‘infinite’ media, by considering a uniaxially stretched plane strip of finite dimensions weakened by a notch of parabolic shape. The distinction between cracks and notches (either sharp or blunt) concerns the engineering community long ago. In fact, it can be stated that the magnificent history of Fractures Mechanics started with the pioneering study of Inglis, related to the stress concentration at the apex of elliptical holes [2]. The importance of flaws in the form of ‘mathematical’ cracks (rather than of notches) was highlighted a few years later [3], when Griffith performed his well-known experiments with specimens made of glass [4]. Based on the data gathered, he concluded that the “... weakness (of the specimens tested) is due almost entirely to minute cracks in the surface ...” [5]. From this instant on, the main challenge was the description of the stress fields that are developed in the immediate vicinity of the tips of ‘mathematical’ cracks or sharp notches (V-shaped notches) and around the ‘crowns’ of blunted notches (either U- or hyperbolically-shaped). The main difference between the two cases is that around the tip of sharp discontinuities the stress field components attain infinite values while around the crown of blunt notches the stresses remain bounded. Although it may be considered a paradox, the problem of the stress field around sharp discontinuities was proven to be easier (at least under the assumption of linear elasticity), and its solution was essentially facilitated by the concept of the Stress Intensity Factor, introduced by Irwin [6]. Using the SIF concept and taking advantage of Westergaard’s work [7], Irwin provided, already at the late fifties, compact and relatively flexible equations for the stress components around the tips of ‘mathematical’ cracks as the first terms of a series expansion. Around the same period, Williams [8, 9] published his seminal papers with the familiar equations for the stress field around sharp V-notches, in terms of eigenfunction series expansions. Almost simultaneously, Neuber [10] dealt with the stress concentration factors for notched bodies, of various geometries under various loading schemes, adopting Airy’s biharmonic potential functions. For the next decades Neuber’s books became reference points for engineers dealing with problems of practical interest. Equations similar to the ones introduced by Williams [8, 9] were presented almost thirty years later by Carpenter [11], who adopted the technique of complex potentials, developed by Kolosov [12] and Muskhelishvili [13]. On the other hand, the respective problem for blunt discontinuities was proven rather tougher. Even today, contributions providing reliable, full-filed, closed form solutions of the problem are mostly welcome. Creager and Paris are, perhaps, the first ones to tackle the problem, in their effort to confront stress corrosion cracking, i.e. “...growth of cracks due to the combined and interrelated action of stress and environment ...” [14]. According to their approach, corrosion cracking is responsible for the generation of discontinuities in the form of “... an elliptical or hyperbolic cylinder ... in which the radius of curvature at the tip is small in comparison to the major dimensions of the void ...” and not in the form of the “... usual plane ending with zero radius of curvature ...” [14]. They concluded that the respective stress field in the immediate vicinity of the crown of the blunt notch can be still described by means of a ‘generalized’ SIF, assuming that the Cartesian reference system is translated from the tip of the ellipse to its focal point. Some twenty years later, simplified expressions for the respective stress field components were obtained by Glinka [15], in terms of the stress concentration factor, the radius of the crown of the notch and the distance from it (keeping the exponents of the distance from the tip of the notch unaltered). Although it was reported that these simplified formulae provided results in good agreement with the respective ones obtained numerically by means of the Finite Element technique, Lazzarin and Tovo [16] indicated that Glinka’s statement about the ‘universality’ of the notch stress fields are incompatible to Williams’ [8, 9] conclusions concerning the dependence of the singularity degree on the opening angle of the notch. In other words, Lazzarin and Tovo [16] suggested that Glinka’s formulae are valid exclusively for notches with zero opening angles. Some years later, the problem was tackled, also, by Nui et al. [17], who presented a solution for the stress field in a finite plane weakened by a single edge V-shaped notch with rounded tip. The solution was achieved using a modified Schwartz-Cristoffel transformation in parallel with Muskhelishvili’s technique [13]. The last decade of the 20th century is signaled by the contributions of Lazzarin’s scientific team, which, are, perhaps, the most influential ones on the issue of the stress fields around either blunt or sharp notches. They dealt not only with the description of the stress components but they contributed, also, in the direction of understanding the conditions causing fracture of notched structural members either under static or fatigue loading schemes. Having as starting point the study of welded lap joints under fatigue conditions [18], they presented flexible approximate solutions for the components of the stress field in the vicinity of open notches, adopting the complex potential’s technique [16]. Later, their approximate solutions [16] were improved, taking into account, also, rounded V-shaped notches with large opening angles [19]. In general, they considered a broad variety of geometrical configurations [20] (including even finite domains [21]), discussing their potential influence on the fatigue limit predictions. In addition, they contributed significantly to the development of Ch. F. Markides et alii, Frattura ed Integrità Strutturale, 68 (2024) 1-18; DOI: 10.3221/IGF-ESIS.68.01 3 the ‘generalized stress intensity factor’ concept, both for rounded V-shaped notches [22] and, also, for U-shaped ones [23]. Later on, their attention was focused on the fracture of notched members, both experimentally and analytically, by means of proper fracture criteria, either energy- or stress-based [24-28]. Concluding this short review, it becomes evident that the interest on the issue of the stress field around various types of notches is continuous and uninterrupted since it was discussed for the first time, due to its paramount importance for the solution of practical engineering problems. The topic is even today under intensive study. In most cases the stress field used is that introduced by the scientific team of late Professor Lazzarin [19, 22]. Indicatively only, one could mention the contributions by Chen and Fan [29], who attempted estimation of fracture loads for blunt U-shaped notches by means of two fracture criteria, adopting the equations for the stress field provided by Lazzarin and Filippi [22]. Recently, Ghadirian et al. [30] determined the mode-I fracture toughness of rock like materials using the notched Brazilian disc configuration and the stress field introduced by Filippi et al. [19]. The same field was used by Sangsefidi et al. [31], who determined the mixed-mode fracture toughness of rocks, again by means of the Brazilian disc test with specimens weakened by U-shaped notches. Nowadays, hybrid schemes are extensively used, combining successfully analytical solutions with experimental protocols and numerical tools [32-34]. In the light of the above discussion, it can be definitely stated that the question concerning the stress field in the vicinity of the crown of notched members is by no means closed. In this context, an alternative approach is presented in this study for the analytical confrontation of the problem, based on a proper conformal mapping and the complex potentials technique [13]. The novelty of the present approach is its capability to deal with finite domains and parabolically shaped notches, independently of the specimen-notch relative dimensions. Moreover, the formulae provided for the components of the stresses are full-field and of closed form. Results of the present solution were comparatively considered against the ones obtained by the respective solution by Filippi et al. [19], for similar geometrical configurations and loading schemes. The agreement is proven quite satisfactory, at least from the qualitative point of view. The same is true for the comparison of the results of the present solution against those of a numerical project that is in progress [35]. THEORETICAL CONSIDERATIONS The problem he first fundamental problem of plane elasticity of a finite strip (length: 2b, width: 2h), weakened by a stress-free edge notch, of parabolic shape L, is analytically explored in this study. The strip is uniaxially stretched by means of a uniform stress distribution σxx=σo (σο>0) that is applied all along its width 2h, as it is shown in Fig.1. The axis of symmetry of the parabola is, also, axis of symmetry of the strip and it coincides with the y-axis of the Cartesian reference system xOy, while the load is applied along the x-axis of the system. Figure 1: The configuration of the problem. Assuming further that the material of the strip is linearly elastic, isotropic and homogeneous, Muskhelishvili’s complex potentials technique is adopted. In this context, use is made of the following analytic function ω(ζ) [13]:             2 2 2 2z ω(ζ) i(ζ iα) ; z x iy i(ξ iη iα) ; x 2(α η)ξ; y ξ (α η) (1) This function maps conformally the region of the plane z=x+iy=reiθ that lies outside of the parabola L with equation: x y   c 2h 2b L T Ch. F. Markides et alii, Frattura ed Integrità Strutturale, 68 (2024) 1-18; DOI: 10.3221/IGF-ESIS.68.01 4  2x 2α α y (2) to the mathematical lower half plane of the complex variable ζ=ξ+iη. In other words, orthogonal parabolas in the z plane (with their focus at the origin O of the (x, y) Cartesian reference system), defining a curvilinear reference system (ξ, η) in the z plane (Fig.2a), correspond to orthogonal lines in the mathematical lower half plane (Fig.2b). In this way, the vertex (or base-point or tip), z(0, –α2) of the parabola L, corresponds to the point ζ(0, 0) in the mathematical plane. A finite part, of dimensions 2bx2h, of the z plane (bounded by red lines in Fig.2a), containing a finite part of the parabola L, represents the region of the edge notched strip under study. (a) (b) Figure 2: The conformal mapping of the actual plane z with the notched strip (a) on to the mathematical lower half plane ζ (b). For the sake of generality, the x-axis is not considered as symmetry axis of the strip, but rather it is located at a distance c from its upper side. In this context, the two end points of the parabolically shaped notch are (–xo, c), (xo, c) (Fig.2a), while the notch itself is reflected to the linear segment (marked red in Fig.2b) defined by the two end points (–ξο, 0), (ξο, 0) in the mathematical plane ζ=ξ+iη. Outline of the method To obtain the solution for the stress field developed in the notched stretched strip, an intact strip of the same dimensions 2bx2h stretched by the uniformly distributed tensile stresses σο>0, is considered first. Obviously, at any point of this intact strip (and, therefore, at any point of an internal parabolic locus L), the stress field is described by a uniaxial tensile stress, σxx=σο>0, while it holds that σyy=σxy=0 (Fig.3a). The solution of this auxiliary problem in terms of the variable ζ in the mathematical plane will be determined in terms of the functions Φ1(ζ), Ψ1(ζ). As a next step, consider the opposite stress field, applied at all the points of the internal locus L, i.e., σxx= –σο, analyzed into components σηη, σξη in the (ξ, η) curvilinear system, mentioned previously (Fig.3b). Let us now set these σηη, σξη as the boundary stresses on the parabolic notch L of a notched strip 2bx2h with stress free sides (Fig.3c), and denote by Φ2(ζ), Ψ2(ζ) the functions solving the specific problem. Then, the solution for the edge notched strip into question, i.e., the finite strip subjected to uniform tension σo on its sides and weakened by a stress free parabolically-shaped notch L, is obtained by simply superposing the previous two solutions (Φ1(ζ), Ψ1(ζ)) and (Φ2(ζ), Ψ2(ζ)), as it is shown schematically in Fig.3d).    1 2 1 2Φ(ζ ) Φ (ζ ) Φ (ζ ), Ψ(ζ ) Ψ (ζ ) Ψ (ζ ) (3) The above procedure and the particular solutions Φ1(ζ), Ψ1(ζ) and Φ2(ζ), Ψ2(ζ), are explicitly described in next sections. In order for the above described method to be better understood it is deemed necessary, to highlight some crucial points of the superposition procedure proposed. The underlying idea is to somehow reach the configuration and solution of a parabolically notched strip by ‘introducing’ the notch to the respective intact stretched strip (Fig.3a). This is here achieved by imposing an opposite stress field –σo at the points of the internal locus L of the intact stretched strip. In fact, such an imposition of –σo on the locus L of the intact stretched strip (Fig.3b), cancels the tensile stress σo along L, rendering the -4 -3 -2 -1 0 1 -4 -3 -2 -1 0 1 2 3 4    oo        a d f  g  m m  g f d a   ζ -20 -15 -10 -5 0 5 10 -20 -15 -10 -5 0 5 10 15 20 x      a d f g y  L ox ox   a d f  g  m m 2 c 2b 2h  η ξ [cm] [cm] z   r z =ω(ζ) z =ω(ζ) Ch. F. Markides et alii, Frattura ed Integrità Strutturale, 68 (2024) 1-18; DOI: 10.3221/IGF-ESIS.68.01 5 area enclosed by L and the upper boundary of the strip stress free (or, equivalently, this area is ‘neutralized’). Actually, this can be conceived as if the specific area had been removed from the strip, thus, transforming it into the notched one. Figure 3: (a) The problem of the stretched intact strip under direct uniaxial tension and the stress field along the internal parabolic locus L; (b) The opposite stress field applied at all the points of the internal locus L; (c) The problem of the notched strip under the opposite stress field on the boundary of the notch L and stress free sides; (d) The stretched, parabolically notched strip obtained by the superposition of the problems (a) and (c). The question arising now is how could this ‘neutralization’ be implemented? A direct approach is by using the principle of superposition, holding for a linearly elastic behavior of the strip, i.e., by superposing to the initial problem of stretching of the intact strip a problem that induces a –σo stress along L. But the –σo stress field acting along the internal locus L could only be achieved by considering the problem of the intact strip acted by σxx=–σo on its loaded sides, which upon being superimposed to the stretched strip (by σxx=σo) would lead to the unstressed strip (i.e., there would be no problem). For this reason, it was necessary to consider the notched strip with stress free sides, loaded along its notch L by the boundary stresses σηη, σξη (Fig.3c), which correspond to the –σo stress field (Fig.3b). Then, the superposition of the problems shown in Figs.3a and 3c, suffices for both the zeroing of stresses along L in the intact strip (indirectly transforming it into the notched one), and, also, for the achievement of the required loading scheme (i.e., that of a tensile uniform stress σo at the loaded boundaries of the strip). Obviously, the above-described procedure is essentially a superposition of complex potentials, Φ1(z)+Φ2(z), Ψ1(z)+Ψ2(z), namely, the perturbation to the stress field due to Φ1(z), Ψ1(z) that is caused by Φ2(z), Ψ2(z), which are responsible for the stress vanishing along L. Furthermore, this procedure is well conceived as an extension of the respective one introduced in Muskhelishvili’s [13] milestone book. Indeed, according to Muskhelishvili, the configuration of a perforated infinite strip is achieved by considering the intact strip, an internal boundary L in the location of the perforation, and, then by mapping conformally the area outside of L onto the infinite plane with the unit hole, and, finally, by demanding L to be free from stresses (thus, rendering the intact infinite strip a perforated one). The above assumptions and the respective analytical solutions introduced in next sections are validated by comparing their outcomes against those provided by well-established solutions, for similar geometries and loading schemes, as well as against a numerical solution in progress. Stretching of the intact strip It is easily seen that the complex potentials solving the problem of stretching an intact strip 2bx2h by a stress σo, which is uniformly applied along its vertical edges (Fig.3a), read as: x y   c 2h 2b L   1 1( ), ( )    x y  c L         x y  L   2 2( ), ( )    2b 2h c x y  1 2( ) ( ) ( )        2b 2h c  1 2( ) ( ) ( )        L (a) (b) (c) (d) Ch. F. Markides et alii, Frattura ed Integrità Strutturale, 68 (2024) 1-18; DOI: 10.3221/IGF-ESIS.68.01 6       ο ο 1 1 1 1 1 1 σ σ Φ (z) Φ (ω(ζ)) Φ (ζ) , Ψ (z) Ψ (ω(ζ)) Ψ (ζ) 4 2 (4) Indeed, substituting Eqns.(4) into the familiar formulae providing stresses in terms of Φ1(z), Ψ1(z) in the actual domain, i.e.:        xx yy 1 yy xy 1 1 1σ σ 4 Φ (z), σ iσ 2 Φ (z) zΦ (z) Ψ (z) (5) (where  denotes the real part; prime denotes the first order derivative and over-bar denotes the conjugate complex value [13]), it follows that at every point of the intact strip it holds that σxx=σο>0, σyy=σxy=0, as it was expected, and, obviously, the same stresses appear all along the internal parabola L of the strip (Fig.3a). Consider now the opposite stress field, i.e.:    xx o yy xyσ σ , σ σ 0 (6) Substituting from Eqns.(6) in the transformation formulae [13]:         2 iβ ηη ξξ yy xx ηη ξξ ξη yy xx xyσ σ σ σ , σ σ 2iσ (σ σ 2iσ )e (7) one obtains the respective stress components in the (η, ξ) curvilinear system, rotated by an angle β with respect to the x-axis (Fig3.b), as:   ηη ξξ oσ σ σ (8)    2 iβ ηη ξξ ξη oσ σ 2iσ σ e (9) Taking into account that [13]: 2iβ Eqns.(1)ω (ζ) ζ iα e ζ iαω (ζ)        (10) Eqn.(9) becomes:      ηη ξξ ξη o ζ iα σ σ 2iσ σ ζ iα (11) Adding Eqns.(8) and (11), and recalling that ζ=ξ+iη, yield:        2 ηη ξη o 2 2 ξ iξ(η α) σ iσ σ ξ (η α) (12) In particular, for η=0, Eqn.(12) provides (after separating the real from the imaginary part) the σηη,L and σξη,L components of the stress field at the points of the internal parabola L, due to the opposite stress field of Eqns.(6), as:      2 ηη,L o ξη,L o2 2 2 2 ξ αξ σ σ , σ σ ξ α ξ α (13) The notched strip with stress free sides loaded by certain stresses along the parabolic notch L Consider now the first fundamental problem of a strip 2bx2h with stress free linear edges, baring a parabolically shaped edge notch L acted by the previous stresses σηη,L and σξη,L of Eqns.(13) (Fig.3c). The solution of this problem, in terms of Ch. F. Markides et alii, Frattura ed Integrità Strutturale, 68 (2024) 1-18; DOI: 10.3221/IGF-ESIS.68.01 7 the variable ζ=ξ+iη (via the mapping of Eqns.(1)), is expressed in terms of Φ2(ζ), Ψ2(ζ). Taking into account Eqns.(1) for η=0, and Eqns.(13), the boundary condition for the stresses on the parabolic notch L reads, in terms of the two functions Φ2 and Ψ2, as follows:          2 2 2 2 2 ο (ξ iα) (ξ iα)Φ (ξ ) (ξ iα)Φ (ξ ) Φ (ξ ) (ξ iα)Ψ (ξ ) σ ξ 2 (14) Following Muskhelishvili [13], the solution to the functional equation of Eqn.(14) is:                 o o o o ξ ξ 2 o o 2 2 2 2 ξ ξ 1 σ ξdξ 1 σ ξdξ ζ iα (ζ iα) Φ (ζ ) , Ψ (ζ ) Φ (ζ ) Φ (ζ ) 2πi(ζ iα) ξ ζ 2πi(ζ iα) ξ ζ ζ iα 2(ζ iα) (15) After some relatively lengthy algebra it is obtained that:        o ο 2 ο ο σ ζ ξ Φ (ζ) ζ log 2ξ 2πi(ζ iα) ζ ξ (16)                       2 2 3 2 2 2 o ο 2 ο3 3 2 ο ο ο σ 3iαζ 6α ζ iα ζ ξ ζ 2iαζ 5α ζ(ζ iα) 1 1 Ψ (ζ) log ξ 2πi 2(ζ iα) ζ ξ (ζ iα) 2(ζ iα) ζ ξ ζ ξ (17) Stretching of the strip with the stress free edge notch L Substituting from Eqns.(4), (16) and (17) in Eqn.(3), the solution of the problem in question (i.e., that of stretching a finite strip 2bx2h with the stress free edge notch L, (Fig.3d or Fig.1)), is obtained as:         o o ο ο ο σ σ ζ ξ Φ(ζ) ζ log 2ξ 4 2πi(ζ iα) ζ ξ (18)                         2 2 3 2 2 2 o o ο ο3 3 2 ο ο ο σ σ 3iαζ 6α ζ iα ζ ξ ζ 2iαζ 5α ζ(ζ iα) 1 1 Ψ(ζ) log ξ 2 2πi 2(ζ iα) ζ ξ (ζ iα) 2(ζ iα) ζ ξ ζ ξ (19) Inversing the transformation of Eqns.(1) yields:   ζ iα iz (20) Substituting for ζ from Eqn.(20) in Eqns.(18) and (19), the solution in terms of the variable z=x+iy=reiθ is obtained as:                o o ο ο ο σ σ i iα iz ξ Φ(z) iα iz log 2ξ 4 2πi z iα iz ξ (21)                          3 2 o o ο ο3 2 3 2 ο 2 3 ο 2 2 ο σ σ 3αz 4iα iα iz ξ iz 4α Ψ(z) i log 2 i ξ 2 4π z ziα iz ξ 5αz iz iz 8α iz 4iα 2ξ ziα iz ξ (22) The stress field in the stretched notched strip At one’s convenience, the stress components in the stretched, notched strip (Fig.1 or Fig.3d), may be expressed either in the Cartesian (x, y) or the curvilinear (ξ, η) reference sytem, via the well-known formulae [13]: Ch. F. Markides et alii, Frattura ed Integrità Strutturale, 68 (2024) 1-18; DOI: 10.3221/IGF-ESIS.68.01 8        yy xy xx yyσ iσ 2 Φ(z) zΦ (z) Ψ(z), σ 4 Φ(z) σ (23)           ηη ξη ξξ ηη ω(ζ) ω (ζ) σ iσ 2 Φ(ζ) Φ (ζ) Ψ(ζ), σ 4 Φ(ζ) σ ω (ζ) ω (ζ) (24) In this context, substituting from Eqns.(21, 22) in Eqns.(23), the explicit expressions for the Cartesian stress components σxx, σyy, σxy, in terms of z=x+iy=reiθ, at any point (r, θ) in the notched strip, are obtained. The respective (rather lengthy) full-field and closed-form expressions are given in Appendix. In particular, regarding the stress field along the notch L, the curvilinear system (ξ, η) appears to be more convenient than the Cartesian one. In this context, substituting from Eqns.(18, 19) into Eqns.(24), setting ζ=ξ (i.e., η=0), and taking into account that by the fourth of Eqns.(1) it holds that ξ=(α2+y)1/2 and ξo=(α2+c)1/2, it follows that:                              2 2 2 2 2 2o ξξ o 2 2 2 α c α y2σ σ σ α α y log 2 α c π α y , α y c π 2α y α c α y (25)  ξη ηησ σ 0 (26) Eqns.(26) were expected, expressing the fulfillment of the notch stress free boundary conditions. Regarding σξξ, as it is seen from Eqn.(25), becomes singular for y=c, i.e., at the end points of the notch L (i.e., at the points of intersection of L with the upper edge of the strip). However, as it will be proven in next sections, σξξ remains bounded even infinitesimally close to the notch end points, e.g. considering a value y≤0.99c, the one that will be adopted in the next applications. EXPLORING THE STRESS FIELD ALONG SOME STRATEGIC LOCI – VALIDATION OF THE SOLUTION n the following examples, the variation of the stress field components is plotted along some strategic loci, i.e., along the boundaries of the strip, along the bisector of the notch and, also, along the flanks of the notch, in an attempt to highlight the potentials of the solution and, also, to explore critical features of the respective distributions, regarding the fulfillment of the boundary conditions imposed. Moreover, the variation of the polar stress components is plotted in the immediate vicinity of the notch base (or notch tip), in order for the present solution to be compared against the well- established one provided analytically and numerically by Filippi et al. [19]. Due to symmetry, all plots are realized only for the half strip. Along the strip sides and the notch bisector (i.e., y-axis), use was made of the expressions of the Appendix. As already mentioned, in order to avoid the singularity at the end points of the notch (±xo, c), y was set equal to 0.99c, instead of c, all along the upper side of the strip. For the plots along the half notch, use was made of Eqn.(25) for the only non-zero stress component σξξ. A strip of dimensions 2bx2h=(30x20) cm was considered. The x-axis was constantly located at a fixed distance c=5 cm from the upper strip side. Three characteristic geometries were considered for the notch, depending on the value of α of Eqn.(2); namely, for α=0.25, 0.5 and 1.0 cm1/2, corresponding to notches becoming gradually deeper (longer) and wider (as it can be seen in the figures of the following sections). In all cases, the strip was stretched by a uniformly distributed tensile stress σxx=σo=10 MPa, applied along its vertical edges (normally to the direction of the bisector of the notch). Stresses along the strip edges, the notch bisector and the notch flanks In Fig.4, the variation of the σxx stress (red color) along the sides of the strip and along the bisector of the notch is plotted, together with the variation of the σξξ stress component (green color) along the flank of the notch, for three characteristic α- values. The stresses are properly adjusted to the length scale considered in the plots, providing a better overview of their variations along the 2bx2h strip. Some characteristic numerical values for the stresses are, also, shown in the diagrams. It is seen that even in the case of a strip of dimensions comparable to those of the notch, the applicability of the present solution is well justified. Clearly, for a bigger strip the solution will perfectly fulfill the boundary conditions imposed. It is worth mentioning, also, that the sign change of the σξξ stress component along the notch flanks, resembles the respective one observed for the distribution along the elliptic or circular holes in the case of the infinite, uniaxially stretched plate. I Ch. F. Markides et alii, Frattura ed Integrità Strutturale, 68 (2024) 1-18; DOI: 10.3221/IGF-ESIS.68.01 9 Figure 4: The variation of σxx along the bisector of the notch (y-axis), along the edges sides of the strip, and the variation of σξξ along half of the perimeter of the notch, for a notch with (a) α=0.25, (b) α=0.5 and (c) α=1.0 cm1/2.  10 [MPa]  41.19  x y 16.42 2h 2b 1.81 8.27 8.27 10.56 10.5610.48 10.48 41.19    xx  xx xx xx   -20 -15 -10 -5 0 5 10 -20 -15 -10 -5 0 5 10 15 20 Notched strip [cm], for α=1.0 cm1/2  10 [MPa]  68.35  x y 14.44 2h 2b 0.55 8.83 8.83 10.39 10.3910.37 10.37 68.35    xx  xx xx xx   -20 -15 -10 -5 0 5 10 -20 -15 -10 -5 0 5 10 15 20 Notched strip [cm], for α=0.5 cm1/2 10 [MPa]  124.59  x y 12.55 2h 2b 0.14 9.00 9.00 10.34 10.3410.35 10.35 124.59     xx  xx xx xx -25 -20 -15 -10 -5 0 5 10 -25 -20 -15 -10 -5 0 5 10 15 20 Notched strip [cm], for α=0.25 cm1/2   (a) (b) (c) Ch. F. Markides et alii, Frattura ed Integrità Strutturale, 68 (2024) 1-18; DOI: 10.3221/IGF-ESIS.68.01 10 The variation of the remaining two Cartesian stress components, σyy and σxy, along the same as previously loci are shown in Fig.5. It is seen that the distributions along the strip edges fulfill (according to a quite satisfactory manner) the boundary conditions imposed. The non-zero values of σyy and σxy along the notch flanks (in the interior) are to suffice the zeroing of the boundary stress components σηη and σξη along the notch flank (approaching the strip from its exterior). Figure 5: The variation of σyy and σxy along the bisector of the notch (y-axis), along the linear edges of the strip and along half of the perimeter of the notch, for a notch with (a) α=0.25, (b) α=0.5 and (c) α=1.0 cm1/2. 10 [MPa]  x y 2h 2b 4.71 0.03 0.03 0.03 0.03 0.00 xy xy xy xy xy 0.00 0.00 11.96 y 20.75   4.71 -20 -15 -10 -5 0 5 10 -20 -15 -10 -5 0 5 10 15 20 Notched strip [cm], for α=1.0 cm1/2 10 [MPa]  x y 12.63 2h 2b 12.63 0.17 0.17 0.25 0.25 0.90 yy yy yy yy yy 5.85 0.90 7.50y 22.09  y 20.43     -20 -15 -10 -5 0 5 10 -20 -15 -10 -5 0 5 10 15 20 Notched strip [cm], for α=1.0 cm1/2 10 [MPa]  xy x y 2h 2b 2.61 0.11 0.11 0.04 0.04 0.00 xy xy xy xy xy 0.00 0.00 20.92 y 20.71   -20 -15 -10 -5 0 5 10 -20 -15 -10 -5 0 5 10 15 20 Notched strip [cm], for α=0.5 cm1/2 10 [MPa]  yy x y 12.59 2h 2b 12.59 0.19 0.19 0.16 0.16 0.54    yy yy yy yy yy 11.02 0.54 14.39y 22.36  y 20.27  -20 -15 -10 -5 0 5 10 -20 -15 -10 -5 0 5 10 15 20 Notched strip [cm], for α=0.5 cm1/2 10 [MPa]  yy x y 11.71 2h 2b 11.71 0.19 0.19 0.13 0.13 0.43    yy yy yy yy yy 21.78 0.43 28.53 y 22.60  y 20.15  -20 -15 -10 -5 0 5 10 -20 -15 -10 -5 0 5 10 15 20 Notched strip [cm], for α=0.25 cm1/2 10 [MPa]  xy x y 2h 2b 1.27 0.14 0.14 0.05 0.05 0.00  xy xy xy xy xy 0.00 0.00 39.27 y 20.69  -20 -15 -10 -5 0 5 10 -20 -15 -10 -5 0 5 10 15 20 Notched strip [cm], for α=0.25 cm1/2(a) (b) (c) Ch. F. Markides et alii, Frattura ed Integrità Strutturale, 68 (2024) 1-18; DOI: 10.3221/IGF-ESIS.68.01 11 The variation of the polar stress in the vicinity of the base (tip) of the notch In Fig.6, the variation of the polar stress components σrr, σθθ, σrθ, is plotted along a semi-circular locus (red color) with its center at the origin of the reference system, enclosing a part of the strip in the vicinity of the base (tip) of the notch. Three different stretched, notched strip configurations are studied in accordance with the parameters set in previous section. The radii r of the circular loci around each one of the three notch bases were arbitrarily chosen as: r=10α2 for α=0.25 cm1/2, resulting to an angle range 53.13ο≤θ≤–90ο (Fig.6a), r=4α2 for α=0.5 cm1/2 (resulting to 30ο≤θ≤–90ο) (Fig. 6b), and r=3α2 for α=1.0 cm1/2 (resulting to 19.47ο≤θ≤–90ο) (Fig.6c). For the calculations of stresses, the general expressions given in the Appendix were used, in conjunction with the transformation formulae for the rotation of the reference system at an angle θ. The distributions of σrr, σθθ, σrθ, resemble qualitatively the respective well established ones given in ref.[19], based on completely different analytical and numerical approaches, thus, justifying further the validity of the present solution. Figure 6: The variation of the σrr, σθθ and σrθ, stress components along a circular locus around the base (tip) of the notch, for notches with (a) α=0.25, (b) α=0.5 and (c) α=1.0 cm1/2. -5 0 5 10 -10 -5 0 5 10 x y 1/ 21.0cm  o90   o19.47 r x2 O y 2r 3   0 6 12 18 -90 -60 -30 0 30 Po la r st re ss c om po ne nt s [M Pa ] Angle θ [deg] rr  r o19.47 -5 0 5 10 15 20 25 -90 -60 -30 0 30 60 P ol ar s tr es s co m po ne nt s [M P a] Angle θ [deg] o53.13 rr r  -5 0 5 10 -10 -5 0 5 10 x y 1/20.25cm  r  2r 10  2 x y o53.13 o90    0 6 12 18 24 -90 -60 -30 0 30 P ol ar s tr es s co m po ne nt s [M Pa ] Angle θ [deg] rr  r-5 0 5 10 -10 -5 0 5 10 x y 1/20.5cm  o90   o30 r x2 O y 2r 4   (a) (b) (c) Ch. F. Markides et alii, Frattura ed Integrità Strutturale, 68 (2024) 1-18; DOI: 10.3221/IGF-ESIS.68.01 12 Stress variation along the notch bisector (y-axis) In the present case, where the strip is stretched normally to the notch bisector, the stress field along the y-axis (aligned along the bisector of the notch) is of particular interest, since it provides among others the Stress Concentration Factor. In this direction, the reduced formulae for the Cartesian stress components are next provided along the y-axis. By setting θ=–π/2 in the full-field expressions of the Appendix, it is obtained that:                                   22 2 1o ο ο xx o 2 2 ο o 2σ α r α r α π α ξ α ξ σ σ tan 1 α 1 π α r ξ 2 rr r r ξ r α (27)                                   22 2 1o ο ο yy 2 2 ο o 2σ α r α r α π α ξ α ξ σ 2 tan 1 α 1 π α r ξ 2 rr r r ξ r α (28) xyσ 0 (29) where ξο=(α2+c)1/2 (see the fourth of Eqns.(1) for η=0 and y=c), and r ranges in the interval α2≤r≤(2h–c). Using these formulae, the distributions of the normal stress components along the y-axis are plotted in Fig.7. As previously, in order to draw these diagrams, a strip of dimensions 2bx2h=(30x20) cm was considered and the x-axis was constantly located at a distance c=5 cm from the upper edge of the strip. Four different parabolic notch geometries were considered, with α= 0.25, 0.5, 1.0, 1.4 cm1/2. The strip was stretched by a uniform distribution of tensile stress σo=10 MPa on its lateral edges. The stress variations were plotted again in an appropriate scale to adjust to the strip dimensions for a better overview (obviously, the results in Fig.7 concerning the strips with for α=0.25, 0.5 and 1.0 cm1/2, coincide with those in Figs.4 and 5, as they are essentially based on the same full-field formulae, given in the Appendix, reduced now to Eqns.(27, 28)). Figure 7: The variation of σxx and σyy stress components along the notch bisector (y-axis), for four notch dimensions. The stress field at the base (tip) of the notch and the stress concentration k In the edge-notched strip, uniaxially stretched normally to the axis of symmetry of the notch (i.e., for the “mode I” notch), the main interest is focused on the point at the base (tip) of the notch, where a significant stress concentration is expected -20 -15 -10 -5 0 5 -20 -15 -10 -5 0 5 10 15 Edge notched strip [cm] 21.78 11.02 5.85 yy,max 4.49   10 [MPa]  124.59 68.35 41.19 xx,max 33.99  xx,min o   10.35 10.37 10.48 10.68 yy,min 0  0.43 0.54 0.90 1.35 2 c 5cm x y 0.25 0.5 1.0 2h 2b 2.25 4.58 9.80 14.77 Notch for 1 .4 xx yy Notches Ch. F. Markides et alii, Frattura ed Integrità Strutturale, 68 (2024) 1-18; DOI: 10.3221/IGF-ESIS.68.01 13 to appear, governing, thus, the failure of the strip. In this context, setting r=α2 in Eqn.(28) yields σyy=0, while Eqn.(27) provides the critical tensile stress σxx,cr at the tip (0, –α2) of the notch, as:       ο xx,cr o 4ξ σ σ 1 πα (30) Taking into account that by the third of Eqns.(1), ο oξ x / 2α , and that by Eqn.(2)  2 οx 2α α c , one obtains:        2 xx ,cr o 4 σ σ 1 c α πα (31) In turn, the stress concentration factor, k, at the tip of the blunt (parabolic) notch is given by the simple expression:    xx,cr 2 o σ 4 k 1 c α σ πα (32) From Eqn.(32) the obvious conclusion is drawn that k is proportional to the 1/2 power of the depth (c+α2) of the notch. The strongest dependence of k is on α, which, apart from the depth, dictates the sharpness of the notch base. The smaller the value of α, the higher the sharpness of the notch tip, reaching for α→0 the limiting case of the ‘mathematical’ edge crack. In the latter case, the Stress Concentration concept shall give its place to the Stress Intensity one and the related Mode-I stress intensity factor KI concept at the crack tip, which clearly constitutes a special limiting case of the present solution. For a fixed value c=5 cm, the dependence of k on α, which rules the depth (length) (c+α2) and span 2xo of the notch, is shown in Fig.8. Clearly, k increases rapidly after a value α of about 0.60. The parameter c is also of crucial importance for the numerical values of k, however, it should be dealt with caution since values of c higher than that considered here, could lead to undesirable bending effects which must be taken into account in the analytical solution. In any case, such effects can be pretty well confronted by assuming long strips, in which case the present solution may be, also, applied without further modifications. Figure 8: The stress concentration factor k versus the parameter α (dictating the dimensions of the notch). DISCUSSION AND CONCLUSIONS n alternative full-field and closed form, analytical solution was presented, for the plane problem of a finite strip that is weakened by a single edge notch of parabolic shape. The strip is assumed to be under uniform, uniaxial stretching. The analysis is based on an extension of the pioneering technique introduced by Mushkelishvili [13] for the solution of the problem of a perforated infinite strip, by considering firstly an intact strip and then ‘neutralizing’ (rendering it stress 0 6 12 18 24 30 0.0 0.2 0.4 0.6 0.8 1.0 1.2 1.4 S tr es s co nc en tr at io n k α [cm1/2] -2.5 0.0 2.5 5.0 -8 -4 0 4 8 Notch [cm] x y 1/20.10cm 0.25 0.40 0.56 0.75 0.97 1.10 1.40 k c 2 A Ch. F. Markides et alii, Frattura ed Integrità Strutturale, 68 (2024) 1-18; DOI: 10.3221/IGF-ESIS.68.01 14 free) a locus identical to that of the empty space of the perforation (by demanding the stresses to be zero along the perforation), transforming, thus, the intact infinite strip into a perforated one. Following this approach, an intact finite strip was firstly considered here under uniaxial tension. Then a part of the strip’s area, resembling an edge parabolically-shaped notch, was ‘neutralized’ by superposing on its boundary an auxiliary stress field, opposite to the one developed due to the action of the tensile stress along the edges of the intact strip. The specific procedure permitted determination of flexible (and relatively compact) formulae for the components of the stress field all over the strip. The main advantage of these formulae (apart from the fact that they are of closed form and full-field) is that they are relieved from the assumption of an ‘infinite’ medium. In other words, there is no need to assume that the dimensions of the notch are ‘small’ compared to those of the strip itself. The specific characteristic of the present solution permits relatively easy parametric analyses of various geometric configurations of the notch and the finite strip. It is mentioned, characteristically that by adjusting properly the parameters governing the shape of the parabolic notch, these formulae degenerate to the respective ones for a strip with a ‘mathematical’ crack. The stress concentration factor, k, was proven to depend mainly on the parameter α, of the equation x=2α(α2+y)1/2 (the equation of the parabolic notch), namely the parameter that governs the sharpness of the notch. It was concluded that for values of α higher than 0.6 the stress concentration factor is relatively small, ranging (almost linearly) in a narrow interval 4> /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