A Methodology to Determine the Effective Plastic Zone Size Around Blunt V-Notches under Mixed Mode I/II Loading and Plane-Stress Conditions

: The determination of the ductile failure behavior in engineering components weakened by cracks and notches is greatly dependent on the estimation of the plastic zone size (PZS) and, particularly, the effective plastic zone size (EPZS). Usually, time-consuming complex elastic–plastic analyses are required for the determination of the EPZS. Such demanding procedures can be avoided by employing analytical methods and by taking advantage of linear elastic analyses. In this sense, this work proposed a methodology for determining the PZS around the tip of blunt V-notches subjected to mixed mode I/II loading and plane-stress conditions. With this aim, ﬁrstly, existing approximate mathematical expressions for the elastic stress ﬁeld near round-tip V-notches reported in the literature are presented. Next, Irwin’s approach (fundamentally proposed for sharp cracks) and a yield criterion (von Mises or Tresca) were applied and are presented. With the aim of verifying the proposed methodology, elastic–plastic ﬁnite element analyses were performed on virtual AISI 304 steel V-notched specimens. It was shown that the analytical formulations presented cannot estimate the complete shape of the plastic zone. However, the EPZS, which is crucial for predicting the type of ductile failure in notched members, can be successfully estimated.


Introduction
The resistance of structures against failure is defined by their load-carrying capacity (LCC), which may be affected by stress risers such as notches of various shapes (e.g., blunt V-notches). When dealing with non-linear engineering materials (e.g., most metallic alloys), ductile failure criteria depend on the yielding behavior of the notched structure, which can be analyzed through the plastic zone size (PZS) and, alternatively, the effective plastic zone size (EPZS) at the tip of the notch. This EPZS is defined as the size of the plastic zone in a direction perpendicular to the maximum principal stress (MPS). Numerous works have been developed in previous decades analyzing the PZS for cracks and notches, some of which are presented below.
Larsson and Carlsson [1] determined the PZS at the crack tip under plane-strain conditions and improved their solution by considering the effects of T-stress on their calculations. Subsequently, the effect of T-stress on the crack tip plastic zone was studied in [2]. Afterwards, Yeh et al. [3] investigated the PZS in isotropic and anisotropic materials by using various failure criteria, and their results were compared with the results obtained in [1,2]. In 2012, Sousa et al. [4] estimated the PZS and showed that T-stress improvements are limited to medium nominal stress to yield strength ratios. Some researchers used the Dugdale model for the estimation of PZS in cracks [5][6][7][8].
For instance, Petroski [9] used weight functions and proposed a new technique requiring only a single elastic stress intensity factor (SIF) and the elastic stress field in the un-cracked body. Another technique was used in [10], in which, with proper crack-dislocation density functions, the singularity was cancelled and the Dugdale PZS was determined. Additionally, the PZS of a Dugdale crack interacting with a circular inclusion was studied in [11].
Since analytical methods are relatively easier and faster than numerical elastic-plastic finite element (FE) methods for studying and determining the PZS, they have been used extensively in the literature for various crack types and conditions, from estimating PZS in a penny-shaped crack embedded in a transversely isotropic medium [12] to investigating the effects of the T-stress on PZS [13,14]. Kiselev and Rivkin [15] used the von Mises yield criterion to analytically study the dependence of the size of the plastic zone on the type of stress state. Chang and Ohr [16] formulated an integral equation to inspect the effects of thickness on PZS. Rogowski [17] proposed a formula considering plane strain conditions to predict the size of the plastic zone in a cracked thin elastic-plastic layer sandwiched between two ideal elastic adherents. This researcher ultimately verified the suggested equations by FEM calculations. Zhang et al. [18] analytically studied the plastic zone size of a stiffened panel with several collinear cracks. They used the weight function method, displacement compatibility, and superposition principle to derive the suggested equations, and, finally, the results were compared to those achieved from FE analysis and the literature.
Sometimes, the mathematical expressions in analytical studies of the PZS become rather long. In these cases, the best way to solve these equations is by writing a code. Bernahou et al. [19] used von Mises and Tresca yield criteria and wrote a code to analytically estimate the PZS at a crack tip under pure mode I, pure mode II, and mixed mode I/II loading conditions. Shi [20] proposed a highly accurate solution procedure to describe the interaction effect among the doubly periodic rectangular-shaped arrays of cracks and studied the interaction effects among the periodic cracks on the PZS and the crack tip opening displacement (CTOD). Many studies have also been undertaken attempting to improve Irwin's model [21,22] or, based on Irwin's approach, to gain a better estimation of PZS [23]. In 2016, Teimoori and Faal [24] used Irwin's approach for the assessment of the PZS for interacting cracks in a circular plane subjected to anti-plane deformation. Regarding 3D analysis of cracks, Caputo et al. [25], proposed for the first time an analytical solution for determining the PZS in short cracks. Some other investigations worth mentioning herein concerning the 3D analysis of the PZS at the crack neighborhood have been reported in [26][27][28].
As already mentioned, FE analysis has also been used for the calculation of the PZS around the crack tip (e.g., [29,30]). Since experimental studies generally require substantial funding, FE analysis has been used for purposes like the validation of the analytical methods or the investigation about the effects of certain factors on the shape and size of the plastic zone. Kudari et al. [31] used elastic-plastic FE analyses in various specimens to assess the shape and size of the plastic zone ahead of the crack tip and compared the results with those derived from the J-integral. Bouiadjra et al. [32] studied, by utilizing non-linear FE analyses, the effect of the presence of micro-cavities ahead of the crack tip on the shape and size of the plastic zone. In 2010, Yi et al. [33] used the maximum crack opening displacement (MCOD) to suggest a new method for estimating the PZS at the crack tip. To determine the relationship between PZS and MCOD, they used FE simulation to obtain a series of PZS and MCOD values. Marques et al. [34] employed 3D elastic and elastic-plastic FE analyses to numerically predict the size, shape, and volume of the plastic zone in modified cracked specimens under pure mode I and mixed mode I/II loading conditions. Many factors can influence the shape and size of the plastic zone. For example, the influence of deformation anisotropy, the elastic mismatch, and the tension along a crack on the plastic zone in an anisotropic body have been reported in [35][36][37], respectively. Alongside analytical and numerical methods for determination of the size and shape of the plastic zone, experimental studies have also been performed. A number of PZS measurement techniques have been reviewed or employed in [38][39][40][41].
Regarding the plastic zone and its effects on fractures under cyclic loading conditions, many studies have been conducted by various researchers (e.g., [42,43]). For instance, Besel and Breitbarth [44] coupled 3D finite element simulation with experimental results and studied different types of plastic zones generated at the tip of a crack under cyclic loading in a commercial aluminum alloy AA2024-T3. They concluded that their suggested approach of combining experimental and numerical results can offer a practical basis for the quantitative energy-based analysis of the damage accumulation process ahead of a crack tip in elastic-plastic materials under cyclic loading conditions.
Concerning the PZS in the vicinity of notches, Dugdale [45] investigated yielding at the end of a slit in a sheet and obtained a relation between the extent of plastic yielding and the applied external load. Pratap and Pandey [46] proposed a modified Smith's model that can be used for the assessment of the PZS in elliptical and circular notches under plane-strain conditions. Youshi et al. [47] evaluated the PZS for cracks originating from a circular or an elliptical notch by using the Dugdale strip yield concept. Chatterjee [48,49] estimated the PZS in notched specimens under pure mode II loading conditions, in order to choose the specimen dimensions. Moreover, Shi and Puls [50] used the modified Neuber's solution together with the slip-line theory to evaluate the maximum tensile stress and the PZS at a shallow notch under plane-strain and mode I loading conditions. Similar to [3], Yeh et al. [51] conducted research on the effect of the T-stress on the damage zone size prediction of notched laminated composites, reaching similar results. Torabi [52] presented some useful equations for obtaining the first and the second estimation (i.e., the Irwin estimation) of the size of the plastic zone along the bisector line of a blunt V-notch under mode I loading and plane-stress conditions, based on the Tresca yield criterion. In addition, the PZS at the neighborhood of V-and U-shaped notches under various loading conditions has been assessed in [53][54][55][56] by utilizing elastic-plastic FE analysis. Additionally, Gomez and Torabi [57] employed the equivalent material concept (EMC) together with a failure criterion to predict the failure of U-notched components in elastic-plastic materials. Torabi and Shahbazian [58] have recently proposed some explicit formulas for estimating the PZS in V-and U-shaped notches under pure mode I loading and plane-stress/strain conditions, based on the von Mises and Tresca yield criteria. They have utilized Irwin's approach and achieved the second estimation of the PZS in such notches. They have also suggested [59] a semi-analytical approach from which the EPZS in U-shaped notches under mixed mode I/II loading can be determined.
Even though PZS can be determined by elastic-plastic FE analysis, this analysis may be complex and lengthy. This is why in this investigation the effective plastic zone size (EPZS) was analytically estimated and applied for the comprehension of the ductile failure regime in the vicinity of round-tip V-notches subjected to mixed mode I/II loading and plane-stress conditions. With this, firstly, the equations of the elastic stress field around the tip of blunt Vnotches, together with the notch stress intensity factors (NSIFs) in mixed mode I/II loading are reported. Then, the elastic stress components were substituted into the von Mises and Tresca yield criteria to obtain the boundary of the plastic zone around the notch border. Next, the second estimation of the PZS was calculated by extending Irwin's approach (initially proposed for sharp cracks under pure mode I loading) to blunt V-notches under mixed mode I/II loading. Finally, the validation of the proposed formulas was performed by means of elastic-plastic FE analyses. This investigation revealed that, for most of the cases being studied, the EPZSs can be obtained from the proposed closed-form solutions, with a good accuracy while avoiding time-consuming elastic-plastic analyses.

Linear Elastic Stress Field for Round-Tip V-Notches under Mixed Mode I/II Loading
Filippi et al. [60] suggested an approximate closed-form expression for the linearelastic stress distributions around a blunt V-notch under mixed mode I/II conditions (see Figure 1). As depicted in Figure 1, the origin of both the Cartesian and polar reference frames is located at the distance r 0 behind the notch tip on the notch bisector line. The resulting equations in the polar coordinate system are the following: where σ θθ , σ rr , σ rθ , K I V,ρ , and K II V,ρ are the tangential stress, the radial stress, the inplane shear stress, and the mode I and the mode II notch stress intensity factors (NSIFs), respectively. The superscript V means that the notch is of a V-shape, and the superscript ρ indicates that the tip of the V-notch is round and not sharp. The functions f ij (θ) and g ij (θ) alongside the eigenvalues λ i and µ i , which depend on the notch opening angle (2α), are reported in Appendix A. The parameters q and r 0 can be obtained from the following equations: where ρ is the notch tip radius. The NSIFs can be calculated by employing the following equations: where ω 1 and ω 2 are two auxiliary parameters (see Appendix A). If their values are known, the NSIFs can be obtained from Equations (8) and (9) [61]: Metals 2021, 11, 1042

of 18
Here, it should be noted that when r is equal to r 0 Equation (9) becomes singular. Here, it should be noted that when r is equal to r0 Equation (9) becomes singular.

Formulation Based on the von Mises Yield Criterion and Under Plane-Stress Conditions
The von Mises yield criterion (also known as the maximum distortion energy criterion) can be written as: where σY is the yield strength of material and σ1, σ2, and σ3 are the three principal stresses. Under plane-stress conditions, where σ3 is zero, Equation (10) is simplified to: The two principal stresses σ1 and σ2 in Equation (11) can be acquired from the twodimensional Mohr's circle as: Now, note that the stress components σθθ, σrr, and σrθ in Equation (12) are given in Equations (1)-(3), respectively, for the case of blunt V-notches. To calculate the plastic zone size (PZS) based on the von Mises yield criterion under mixed mode I/II loading and plane-stress conditions, first, Equations (1)-(3) should be substituted into Equation (12), obtaining the two principal stresses (i.e., σ1 and σ2). Then, these two stresses should be entered in Equation (11). Here, it should be noted that in the final form of Equation (11) r should be replaced by rY, which is the first estimation of the PZS in the vicinity of a roundtip V-notch, at any angular coordinate θ, and under the mentioned loading and stress conditions. Moreover, it is worth mentioning that by a simple coding in a proper software (e.g., MATLAB), all of these calculations can be performed conveniently and rapidly. In addition, it is crucial to understand that the corresponding values of the NSIFs, KI V,ρ and KII V,ρ , can be computed by a linear elastic FE analysis on the V-notched component being analyzed.

Formulation Based on the Tresca Yield Criterion and under Plane-Stress Conditions
The Tresca yield criterion (also known as the maximum shear stress criterion) can be expressed as:

Formulation Based on the von Mises Yield Criterion and under Plane-Stress Conditions
The von Mises yield criterion (also known as the maximum distortion energy criterion) can be written as: where σ Y is the yield strength of material and σ 1 , σ 2 , and σ 3 are the three principal stresses. Under plane-stress conditions, where σ 3 is zero, Equation (10) is simplified to: The two principal stresses σ 1 and σ 2 in Equation (11) can be acquired from the twodimensional Mohr's circle as: Now, note that the stress components σ θθ , σ rr , and σ rθ in Equation (12) are given in Equations (1)-(3), respectively, for the case of blunt V-notches. To calculate the plastic zone size (PZS) based on the von Mises yield criterion under mixed mode I/II loading and plane-stress conditions, first, Equations (1)-(3) should be substituted into Equation (12), obtaining the two principal stresses (i.e., σ 1 and σ 2 ). Then, these two stresses should be entered in Equation (11). Here, it should be noted that in the final form of Equation (11) r should be replaced by r Y , which is the first estimation of the PZS in the vicinity of a round-tip V-notch, at any angular coordinate θ, and under the mentioned loading and stress conditions. Moreover, it is worth mentioning that by a simple coding in a proper software (e.g., MATLAB), all of these calculations can be performed conveniently and rapidly. In addition, it is crucial to understand that the corresponding values of the NSIFs, K I V,ρ and K II V,ρ , can be computed by a linear elastic FE analysis on the V-notched component being analyzed.

Formulation Based on the Tresca Yield Criterion and under Plane-Stress Conditions
The Tresca yield criterion (also known as the maximum shear stress criterion) can be expressed as: Likewise in the von Mises yield criterion, the principal stresses σ 1 and σ 2 can be obtained from Equation (12). By presuming that σ 3 = 0 (i.e., plane-stress conditions), and that σ 1 ≥ σ 2 , three different conditions may occur: First, σ 1 ≥ σ 2 ≥ 0. In this case, the Tresca yield criterion becomes: Second, σ 1 ≥ 0 ≥ σ 2 . Then, Equation (13) can be written as: Third, 0 ≥ σ 1 ≥ σ 2 . On this occasion, Equation (13) can be written as: Similarly to Section 2.2.1, herein, the first estimation of the PZS based on the Tresca yield criterion under mixed mode I/II loading and plane-stress conditions can be obtained by substituting the principal stresses σ 1 and σ 2 (derived from Equation (12)) into Equations (14)- (16), where r must be replaced by r Y (the first estimate of the PZS). Although the final resulting equations may be formally complex, a simple coding in well-known software may easily solve them.

Second Estimation of the Plastic Zone Size (PZS)
By extending Irwin's approach from sharp cracks subjected to pure mode I loading [62] to round-tip V-notches subjected to mixed mode I/II loading, the second estimation of the plastic zone size for blunt V-notches was obtained and presented here. Irwin's approach ignores the effects of the material strain-hardening in the formation of the plastic zone. Hence, the tensile behavior of a ductile material is simplified by assuming an elastic perfectly plastic model (EPPM). Figure 2a shows a typical true tensile stress-strain curve of an elastic-plastic (ductile) material. In this figure, the parameters σ Y , σ f , and ε f are the yield strength, the true fracture stress, and the strain-to-failure, respectively. To simplify this curve into an elastic perfectly plastic model, a curve like that shown in Figure 2b needs to be drawn in such a way that the area under the curve is equal to the area under the curve of Figure 2a. In addition, the elastic parts of both curves and the value of ε f should be identical. This means that the value of σ 0 (the average flow stress) in Figure 2b will always be greater than the value of σ Y in Figure 2a. Note that, generally, in Irwin's model σ 0 is used instead of σ Y .
Likewise in the von Mises yield criterion, the principal stresses σ1 and σ2 can be obtained from Equation (12). By presuming that σ3 = 0 (i.e., plane-stress conditions), and that σ1 ≥ σ2, three different conditions may occur: First, σ1 ≥ σ2 ≥ 0. In this case, the Tresca yield criterion becomes: Second, σ1 ≥ 0 ≥ σ2. Then, Equation (13) can be written as: Third, 0 ≥ σ1 ≥ σ2. On this occasion, Equation (13) can be written as: Similarly to Section 2.2.1, herein, the first estimation of the PZS based on the Tresca yield criterion under mixed mode I/II loading and plane-stress conditions can be obtained by substituting the principal stresses σ1 and σ2 (derived from Equation (12)) into Equations (14)-(16), where r must be replaced by rY (the first estimate of the PZS). Although the final resulting equations may be formally complex, a simple coding in well-known software may easily solve them.

Second Estimation of the Plastic Zone Size (PZS)
By extending Irwin's approach from sharp cracks subjected to pure mode I loading [62] to round-tip V-notches subjected to mixed mode I/II loading, the second estimation of the plastic zone size for blunt V-notches was obtained and presented here. Irwin's approach ignores the effects of the material strain-hardening in the formation of the plastic zone. Hence, the tensile behavior of a ductile material is simplified by assuming an elastic perfectly plastic model (EPPM). Figure 2a shows a typical true tensile stress-strain curve of an elastic-plastic (ductile) material. In this figure, the parameters σY, σf, and εf are the yield strength, the true fracture stress, and the strain-to-failure, respectively. To simplify this curve into an elastic perfectly plastic model, a curve like that shown in Figure 2b needs to be drawn in such a way that the area under the curve is equal to the area under the curve of Figure 2a. In addition, the elastic parts of both curves and the value of εf should be identical. This means that the value of σ0 (the average flow stress) in Figure 2b will always be greater than the value of σY in Figure 2a. Note that, generally, in Irwin's model σ0 is used instead of σY. The elastic tangential stress (σθθ) distribution on a blunt V-notch bisector line versus the distance from the notch tip is schematically shown in Figure 3. As stated before, the strain-hardening is neglected. For this reason, the stress does not rise to any values greater The elastic tangential stress (σ θθ ) distribution on a blunt V-notch bisector line versus the distance from the notch tip is schematically shown in Figure 3. As stated before, the strain-hardening is neglected. For this reason, the stress does not rise to any values greater than σ Y , meaning that the effect of zone A in Figure 3 should be taken into consideration by redistributing it into region B. than σY, meaning that the effect of zone A in Figure 3 should be taken into consideration by redistributing it into region B. This can be done by means of Equation (17), which is suitable for round-tip V-notches and is derived analogically from Irwin's equation for sharp cracks.
where σθθ(r,θ) is the tangential stress at any arbitrary θ angle, and rp is the second estimation of the PZS. In Equation (17), rY and r0 are known, whereas rp and θ are unknown. rY can be obtained from solving Equation (11) or Equations (14)- (16) according to the von Mises and Tresca yield criteria, respectively; r0 can be calculated from Equation (5). Moreover, by considering any arbitrary value of θ, the second estimation of the PZS inside the material can be determined in that specified direction. This makes it much easier to compare the present analytical results and the numerical results obtained from the FE analyses. Evidently, obtaining the shape and size of the plastic zone around the notch round border relies on the determination of the PZS at various θ angles. Another important issue regarding Equation (17) that needs to be discussed is the proper stress term to be used. Note that in Equation (17), when using any yield criterion the stress term that reaches the material yield strength should be used instead of σθθ(r,θ). Here, when analyzing blunt V-notches under mixed mode I/II loading and plane-stress conditions, and applying the von Mises and Tresca yield criteria, the left hand side of Equations (11) and (14) to (16) should be employed instead of σθθ(r,θ), respectively.

Verification in the Particular Case of Sharp Cracks
To the authors' best knowledge, the proposed approach has not yet been suggested in the literature, and thus, it is new. Hence, its validation and verification seem necessary. There are no available solutions for the plastic zone geometry in V-notches. However, the shape and the size of the plastic zone around sharp cracks for both the von Mises and Tresca yield criteria (under mixed mode I/II loading and plane-stress conditions) have been well studied in the literature (e.g., see [63]). Therefore, the proposed approach will be applied to this particular case, something that implies considering that both the (2α) and the notch tip radius (ρ) are equal to zero. In [63], the ratio of the mode II stress intensity factor (SIF) to the mode I SIF (Kr1) was defined as: This can be done by means of Equation (17), which is suitable for round-tip V-notches and is derived analogically from Irwin's equation for sharp cracks.
where σ θθ (r,θ) is the tangential stress at any arbitrary θ angle, and r p is the second estimation of the PZS. In Equation (17), r Y and r 0 are known, whereas r p and θ are unknown. r Y can be obtained from solving Equation (11) or Equations (14)- (16) according to the von Mises and Tresca yield criteria, respectively; r 0 can be calculated from Equation (5). Moreover, by considering any arbitrary value of θ, the second estimation of the PZS inside the material can be determined in that specified direction. This makes it much easier to compare the present analytical results and the numerical results obtained from the FE analyses. Evidently, obtaining the shape and size of the plastic zone around the notch round border relies on the determination of the PZS at various θ angles. Another important issue regarding Equation (17) that needs to be discussed is the proper stress term to be used. Note that in Equation (17), when using any yield criterion the stress term that reaches the material yield strength should be used instead of σ θθ (r,θ). Here, when analyzing blunt V-notches under mixed mode I/II loading and plane-stress conditions, and applying the von Mises and Tresca yield criteria, the left hand side of Equations (11) and (14) to (16) should be employed instead of σ θθ (r,θ), respectively.

Verification in the Particular Case of Sharp Cracks
To the authors' best knowledge, the proposed approach has not yet been suggested in the literature, and thus, it is new. Hence, its validation and verification seem necessary. There are no available solutions for the plastic zone geometry in V-notches. However, the shape and the size of the plastic zone around sharp cracks for both the von Mises and Tresca yield criteria (under mixed mode I/II loading and plane-stress conditions) have been well studied in the literature (e.g., see [63]). Therefore, the proposed approach will be applied to this particular case, something that implies considering that both the (2α) and the notch tip radius (ρ) are equal to zero. In [63], the ratio of the mode II stress intensity factor (SIF) to the mode I SIF (Kr 1 ) was defined as:  Figures 4b and 5b are taken from [63]. It can be observed how the shape of the plastic zones obtained from the proposed approach, and for both the von Mises and Tresca criteria, are very similar to those shown in literature. .

Verification of the Analytical Estimations by Finite Element Analysis
It is unquestionable that comparing the results obtained from the proposed analytical approach with the actual shape and size of the plastic zone observed in experimental tests is more genuine and precise. Nevertheless, these tests are usually costly and laborious. Because of the availability of reliable commercial FE software, numerical validation where KI and KII are the applied stress intensity factors in mode I and mode II loading conditions, respectively. Figures 4 and 5 represent the shapes of the plastic zone for a sharp crack based on the von Mises and Tresca yield criteria, respectively, and under plane-stress conditions for different values of Kr1. Figures 4a and 5a represent the results provided by the methodology proposed in this work, particularized to sharp cracks, whereas Figures 4b and 5b are taken from [63]. It can be observed how the shape of the plastic zones obtained from the proposed approach, and for both the von Mises and Tresca criteria, are very similar to those shown in literature. .

Verification of the Analytical Estimations by Finite Element Analysis
It is unquestionable that comparing the results obtained from the proposed analytical approach with the actual shape and size of the plastic zone observed in experimental tests is more genuine and precise. Nevertheless, these tests are usually costly and laborious. Because of the availability of reliable commercial FE software, numerical validation

Verification of the Analytical Estimations by Finite Element Analysis
It is unquestionable that comparing the results obtained from the proposed analytical approach with the actual shape and size of the plastic zone observed in experimental tests is more genuine and precise. Nevertheless, these tests are usually costly and laborious. Because of the availability of reliable commercial FE software, numerical validation emerges as an alternative to experimental testing for the verification or validation of analytical models. Here, several two-dimensional (2D) finite element (FE) models were created and analyzed by utilizing the commercial code ABAQUS (Dassault Systèmes Simulia Corp., Johnston, RI, USA). Note that FE analyses are performed for two purposes: firstly, to compute the NSIFs; secondly, to determine (numerically) the PZS. The former can be obtained from linear elastic analyses, whereas the latter may be derived from elastic-plastic analyses by using the true material stress-strain curve.
Several tensile rectangular specimens weakened by inclined blunt V-notches, as shown in Figure 6a, were analyzed under plane-stress conditions. Here, it is important to mention that in reality the plane-stress state does not occur, and this idealization is only considered when the simplification from 3D to 2D is in mind. Generally, it is presumed that when two dimensions are far greater than the other dimension, the plane-stress analysis is sufficiently accurate. As shown in Figure 6a, the length and the width of the specimens were far greater than their thickness. Thus, the plane-stress assumption seems reasonable. According to Figure 6a, it is obvious that, since the notch rotation angle (β) was greater than zero and smaller than 90 • , the specimen was under mixed mode I/II loading conditions. In addition, note that for all cases the notch length was equal to 30 mm. Altogether, eight 2D specimens with two different mode mixity ratios and various notch tip radii were simulated. Four specimens with the notch rotation angle equal to 25 • were studied under tensile loads (P) of 35, 40, 44, and 48 kN and notch tip radii of 0.5, 1, 2, and 4 mm, respectively. In the other four specimens, the notch rotation angle was 50 • , with the applied loads being equal to 40, 43, 50, and 55 kN for the notch tip radii of 0.5, 1, 2, and 4 mm, respectively.
In order to achieve a reliable result from numerical analyses, several simulations must be repeated using different numbers of elements, examining the mesh-sensitivity. As is well known, the stress gradient is high in the vicinity of the notch tip. For this reason, refined eight-node quadratic elements with reduced integration are used in this region. In this work, the mesh-sensitivity in the notched model was checked based on the size of the element at the notch neighborhood rather than the total number of elements. One of the advantages of such an approach is that the reader can easily understand how the mesh is refined in the vicinity of the notch, whereas by presenting the total number of elements, the conditions of the meshed region around the notch would not be clear. It was found that when the element size at the notch neighborhood was smaller than 0.1 mm, optimal results can be achieved. However, for the sake of achieving the highest possible accuracy, elements with a size of 0.05 mm were utilized in the FE analyses hereunder. The bottom line of the model was meticulously fixed and a concentrated tensile force (P) according to Figure 6a was exerted on top of the specimen. For all cases studied, the Poisson's ratio (ν) and the V-notch opening angle (2α) were considered to be equal to 0.3 and 45 • , respectively. As mentioned before, the true stress-strain curve of the material was employed as the material property. The material chosen for all the specimens simulated was AISI 304 steel. Figure 6b illustrates the corresponding tensile stress-strain curve. The elastic modulus, the yield strength, and the ultimate strength of this steel were equal to 200 GPa, 342 MPa, and 702 MPa, respectively [64], with σ 0 (easily derived by following Section 2.3) equal to 611.5 MPa.

Results and Discussion
A comparison between the plastic zones in the vicinity of the blunt V-shaped notches subjected to mixed mode I/II loading and plane-stress conditions obtained from the suggested analytical approach (rp) with those derived from the elastic-plastic FE analyses is

Results and Discussion
A comparison between the plastic zones in the vicinity of the blunt V-shaped notches subjected to mixed mode I/II loading and plane-stress conditions obtained from the suggested analytical approach (r p ) with those derived from the elastic-plastic FE analyses is made in this section. With the aim of checking the accuracy of the suggested formulations, a discrepancy index (∆) was defined as follows: where (PZS) FE and (PZS) Analyt. are the PZSs achieved from the FE analysis and the analytical approach, respectively. Here, it should be noted that even though in this section the complete shapes of the plastic zones around the blunt V-notches were obtained from the proposed methodology (also referred to as the closed-form solution) and from FE analyses, only the effective plastic zone size was taken into consideration for comparing the results using Equation (19). The EPZS is important because it can be used to determine the ductile failure regime in notched components. Note that crack initiation on the round tip of blunt V-notches (or U-notches) usually starts from the MPS point [53][54][55][56]. Moreover, it is important to remember that the points of MPS and maximal tangential stress (MTS) are very close. Geometrically, since the blunt tip of a V-notch is a part of a circumference, any perpendicular line to the notch tip will pass through the center of curvature. This is why in the figures of this section both the MPS point and the center of curvature are pointed out with a plus sign inside of a circle. Figures 7-10 show the plastic zones around the V-notches of the specimens described above. These plastic zones were obtained from both the proposed methodology (closed-form solution based on Irwin's approach) and the FE analyses under the plane-stress assumption. Figure 7 illustrates the plastic zones based on the von Mises yield criterion when the notch rotation angle (β) was equal to 25 • and the notch tip radii were 0.5, 1, 2, and 4 mm; Figure 8 shows the results for the same geometries but considering the Tresca yield criterion. Figures 9 and 10 represent the obtained plastic zones in the specimens when β = 50 • . The plastic zones in Figure 9 were based on the von Mises yield criterion, whereas the plastic zones in Figure 10 were derived assuming the Tresca yield criterion. As mentioned in Section 2.3, in Irwin's model, σ 0 is normally used instead of σ Y . However, herein, σ Y was utilized for the analytical calculation of the PZS. The reason behind this is explained below. For both the von Mises and Tresca yield criteria, the discrepancies between the numerical and the analytical results (∆, Equation (19)) are reported in Table 1.
According to Table 1, when β was equal to 25 • , even though the results for both the von Mises and Tresca yield criteria were generally acceptable, the results obtained based on the von Mises yield criterion were more accurate for almost all the notch tip radii. On the contrary, when the notch rotation angle was equal to 50 • , the results for the Tresca yield criterion were more precise.
At this point, some discussion on the reasons that cause the discrepancies between the two results (i.e., FE and analytical) seems important. Basically, they arise from the fact that it is extremely difficult to guarantee the accuracy of the stress solutions out of the direction normal to the maximum tangential stress, given that such accuracy depends on several factors. For example, because of the nature of Equations (1)-(3), which are accurate only in regions with significant stress concentrations, the precision of the solution decreases when going further away from the notch border [60]. It is for this reason that, for most parts in this study, the value of the applied load was not too great, so the generated PZS was small or moderate. Moreover, there is a considerable difference between the actual boundary of the round-tip V-notch and the boundary defined for the stress solution. Based on [60], the proposed stress solutions are satisfied on a limited number of points on the notch edge, and not on the entire edge. In addition, the precision declines when the ratio of the notch length to the notch tip radius decreases [60]. creases when going further away from the notch border [60]. It is for this reason that, for most parts in this study, the value of the applied load was not too great, so the generated PZS was small or moderate. Moreover, there is a considerable difference between the actual boundary of the round-tip V-notch and the boundary defined for the stress solution. Based on [60], the proposed stress solutions are satisfied on a limited number of points on the notch edge, and not on the entire edge. In addition, the precision declines when the ratio of the notch length to the notch tip radius decreases [60].  Earlier in this article, it was stated that to solve the proposed formulas (i.e., obtaining rp), five parameters need to be known. These parameters are σY, θ, σ0, and the NSIFs (KI V,ρ and KII V,ρ ). Here, the yield strength of the material (σY) is known, the average flow stress (σ0) can be calculated as explained in Section 2.3, and θ is an arbitrary angle. If it is in- Earlier in this article, it was stated that to solve the proposed formulas (i.e., obtaining r p ), five parameters need to be known. These parameters are σ Y , θ, σ 0 , and the NSIFs (K I V,ρ and K II V,ρ ). Here, the yield strength of the material (σ Y ) is known, the average flow stress (σ 0 ) can be calculated as explained in Section 2.3, and θ is an arbitrary angle. If it is intended to determine the PZS in a given direction, the definition of θ is straightforward. The NSIFs can be calculated from linear elastic FE analyses. In other words, only the elastic part of the stress-strain curve of the analyzed material is important when dealing with the calculation of the NSIFs. Now, suppose that two materials (material A and material B) with identical elastic parts, and also that one of them (material A) develops large ductility and hardening (like the AISI 304 steel used in this study), whereas the other one (material B) has little hardening and low ductility (it behaves as an EPPM with low strain-to-failure; see Section 2.3). In this case, by considering the same applied load and the same geometry in the notched specimens, the values of the NSIFs would be equal, but the value of the average flow stress would be different, given that the areas under the stress-strain curves are significantly different (obviously, the PZSs would be different too). This also means that the difference between σ Y and σ 0 would be significant for material A and small for material B. In other words, the larger the area under the curve, the greater σ 0 . At the end, when material A is used, the proposed methodology provides estimations of PZS that are notably smaller than those obtained with material B. This shows that the accuracy of the proposed formulation depends on the hardening and the strain-to-failure of the material being analyzed, and it decreases when both parameters increase (and vice versa). Similar results have been reported in [58]. In this sense, it should be noted here that, according to Figure 6b, the hardening and the strain-to-failure of AISI 304 steel was high. Therefore, it was expected that when σ0 is used instead of σY, the PZS calculated from the suggested closed-form solution will become remarkably small in comparison to the size of the plastic zone obtained from the FE analyses. Figure 11 depicts the PZSs obtained from the numerical and the proposed analytical solutions for a blunt V-notched rectangular specimen based on the von Mises yield criterion under plane-stress conditions and using σ0 instead of σY. It is evident that when In this sense, it should be noted here that, according to Figure 6b, the hardening and the strain-to-failure of AISI 304 steel was high. Therefore, it was expected that when σ 0 is used instead of σ Y , the PZS calculated from the suggested closed-form solution will become remarkably small in comparison to the size of the plastic zone obtained from the FE analyses. Figure 11 depicts the PZSs obtained from the numerical and the proposed analytical solutions for a blunt V-notched rectangular specimen based on the von Mises yield criterion under plane-stress conditions and using σ 0 instead of σ Y . It is evident that when σ 0 is used instead of σ Y , the proposed analytical formulation provides a significantly smaller PZS than that derived from the FE analyses. Thus, the authors recommend considering σ Y instead of σ 0 when the material being analyzed develops significant hardening and strain-to-failure.     Figure 11. EPZSs obtained from FE analyses and the proposed methodology (closed-form solution). Notch radius and notch rotation angle were 0.5 mm and 25°, respectively. σ0 was used instead of σY in the analytical solution.

Conclusions
The main concluding remarks of this study are presented in the four bullet points below: • The stress solution around round-tip V-notches under mixed mode I/II loading is an approximation to the actual stress state. This research revealed that it can successfully be used to formulate the plastic zone size (PZS) in the neighborhood of blunt Vshaped notches under plane-stress conditions and based on the von Mises and Tresca yield criteria.

•
The accuracy of the stress solutions out of the direction normal to the maximum tangential stress (MTS) depends on many factors; namely and most importantly, the stress field only satisfies the boundary conditions on the limited area on the round edge of the notch. This makes the shape of the plastic zones acquired from the analytical solution different from that obtained from the FE analyses. However, the effective plastic zone size (EPZS), which is important for the determination of the ductile failure regime in notched components, can be accurately and conveniently obtained by using the proposed analytical formulas.

•
This research showed that when the notch rotation angle (β) was equal to 25°, the accuracy of the suggested analytical solutions for calculating the EPZS was better when the von Mises yield criterion is used. However, when β was equal to 50°, the Tresca yield criterion provided finer results.

•
The proposed methodology was derived from Irwin's approach, which neglects the strain-hardening of the material being analyzed and assumes elastic perfectly plastic behavior. When the plasticity of the material is high, in terms of strain hardening and strain-to-failure, the difference between the average flow stress (σ0) and the yield stress (σY) increases, and the methodology provides worse estimations of the PZS. To avoid this, it is suggested to use σY when the material being analyzed develops significant hardening and strain-to-failure. Otherwise, the average flow stress should be used.

Conclusions
The main concluding remarks of this study are presented in the four bullet points below:

•
The stress solution around round-tip V-notches under mixed mode I/II loading is an approximation to the actual stress state. This research revealed that it can successfully be used to formulate the plastic zone size (PZS) in the neighborhood of blunt Vshaped notches under plane-stress conditions and based on the von Mises and Tresca yield criteria.

•
The accuracy of the stress solutions out of the direction normal to the maximum tangential stress (MTS) depends on many factors; namely and most importantly, the stress field only satisfies the boundary conditions on the limited area on the round edge of the notch. This makes the shape of the plastic zones acquired from the analytical solution different from that obtained from the FE analyses. However, the effective plastic zone size (EPZS), which is important for the determination of the ductile failure regime in notched components, can be accurately and conveniently obtained by using the proposed analytical formulas.

•
This research showed that when the notch rotation angle (β) was equal to 25 • , the accuracy of the suggested analytical solutions for calculating the EPZS was better when the von Mises yield criterion is used. However, when β was equal to 50 • , the Tresca yield criterion provided finer results.

•
The proposed methodology was derived from Irwin's approach, which neglects the strain-hardening of the material being analyzed and assumes elastic perfectly plastic behavior. When the plasticity of the material is high, in terms of strain hardening and strain-to-failure, the difference between the average flow stress (σ 0 ) and the yield stress (σ Y ) increases, and the methodology provides worse estimations of the PZS. To avoid this, it is suggested to use σ Y when the material being analyzed develops significant hardening and strain-to-failure. Otherwise, the average flow stress should be used. Funding: This research received no external funding.

Data Availability Statement:
The data presented in this study are available in the main text of this document.

Conflicts of Interest:
The authors declare no conflict of interest.

Appendix A
Functions used in the stress field for round-tip V-shaped notches (modes I and II) [60]: Expressions for parameters ω 1 and ω 2 [60,61]: