Mixed Mode Fracture Investigation of Rock Specimens Containing Sharp V-Notches

This work aims to assess both experimentally and analytically the fracture behavior of rock specimens containing sharp V-notches (SV-notches) subjected to mixed mode I/II loading. To this end, firstly, several mixed mode fracture tests were conducted on Brazilian disk specimens weakened by an SV-notch (SVNBD sample), performed in their corresponding center and with various notch opening angles. Secondly, the fracture resistance of the tested samples was predicted using a criterion named MTS-FEM. This approach is based on the maximum tangential stress (MTS) criterion, in which the tangential stress is determined from the finite element method (FEM). Additionally, in the present research, the required critical distance is calculated directly from finite element analyses performed on cracked samples. Comparing the experimental results and the analytical predictions, it is shown that the fracture curves obtained from the MTS-FEM criterion are in agreement with the experimental results. These results are achieved without the need for the calculation of stress series expansion coefficients, as an additional advantage of the proposed approach.


Introduction
The existence of different natural defects such as holes, grain boundaries, pores, cracks, and notches define the mechanical behavior of rock masses. Moreover, diverse artificial notch-like defects with different shapes may be created in rock masses when dams or tunnels are built, or when gas and oil wells and excavations are performed. All these defects act as stress risers and play a key role in fracture initiation processes. Therefore, the fracture of notched and cracked samples made of rock has been extensively investigated in the literature. For this purpose, fracture mechanics, which assess the fracture resistance of cracked or notched parts, have been frequently utilized. Since the stress concentration in the sharp V-notch (SV-notch) is more pronounced than that caused by other types of notches, this shape is more dangerous for rock structures, making the corresponding prediction of the load-bearing capacity in rock samples containing SV-notches even more imperative.
In the context of fracture mechanics, there are two main failure modes for notched specimens: (1) pure mode I, or opening mode, in which the notch flanks open without any sliding; and (2) pure mode II, or in-plane shearing mode, in which the notch flanks slide relative to each other without any opening or closing. Generally, due to the arbitrariness of the applied load direction relative to the notch orientation, the combinations of opening and in-plane shearing modes can appear in notch samples. Although several brittle fracture criteria can be found in the fracture mechanics literature for predicting the fracture load of SV-notched samples, three criteria have received more attention: the strain energy density that the results predicted by the MTS-FEM criterion are in agreement with those obtained from testing SVNBD specimens.

Material, Tests and Experimental Results
As mentioned above, there were no experimental data for investigating the mixed mode fracture resistance of rock SV-notched samples. Hence, in this work several Brazilian disks containing a central SV-notch (SVNBD samples) were prepared and tested under mixed mode I/II loading.
As shown schematically in Figure 1, the SVNBD specimen is a circular disk with dimeter and thickness D and t, respectively. There is a central rhombus with the large diagonal being d and the corner angle being 2α. The upper and lower corners of the rhombic hole in the SVNBD samples can be considered sharp V-notches. A compressive load (P) is applied along the diameter of the disk whose direction relative to the notch bisector line, i.e., the loading angle β, defines the mode mixity ratio. The alignment of the applied load with the notch bisector line, i.e., β = 0, gives the pure mode I loading. By rotating the notch bisector line relative to the applied load, the loading condition changes from pure mode I to mixed mode I/II, until a specific angle β II , for pure mode II loading, is achieved. Torabi and Taherkhani [23] calculated the values of β II as a function of the notch length ratio (d/D) and the notch opening angle 2α using FEM. The dimensions and loading conditions of the SVNBD samples tested in the present study are listed in Table 1. It is noteworthy that the loading angle β for pure mode II (i.e., β II ) for each notch opening angle was extracted from the study carried out by Torabi and Taherkhani [23]. Additionally, the pure mode II loading angle β II depends on the notch angle 2α and, therefore, the values of β II for notch opening angles 2α = 30 • , 60 • , and 90 • are different. The loading angles β in mixed mode loading were defined accordingly to the value of β II for each opening angle 2α. criterion, is a particular form of the MTS criterion in which the tangential stress around the SV-notch border is determined using finite element analysis. In addition, the critical distance, which is an important parameter in MTS-FEM criterion, is calculated directly from FEM. It is revealed that the results predicted by the MTS-FEM criterion are in agreement with those obtained from testing SVNBD specimens.

Material, Tests and Experimental Results
As mentioned above, there were no experimental data for investigating the mixed mode fracture resistance of rock SV-notched samples. Hence, in this work several Brazilian disks containing a central SV-notch (SVNBD samples) were prepared and tested under mixed mode I/II loading.
As shown schematically in Figure 1, the SVNBD specimen is a circular disk with dimeter and thickness D and t, respectively. There is a central rhombus with the large diagonal being d and the corner angle being 2α. The upper and lower corners of the rhombic hole in the SVNBD samples can be considered sharp V-notches. A compressive load (P) is applied along the diameter of the disk whose direction relative to the notch bisector line, i.e., the loading angle β, defines the mode mixity ratio. The alignment of the applied load with the notch bisector line, i.e., β = 0, gives the pure mode I loading. By rotating the notch bisector line relative to the applied load, the loading condition changes from pure mode I to mixed mode I/II, until a specific angle βII, for pure mode II loading, is achieved. Torabi and Taherkhani [23] calculated the values of βII as a function of the notch length ratio (d/D) and the notch opening angle 2α using FEM. The dimensions and loading conditions of the SVNBD samples tested in the present study are listed in Table 1. It is noteworthy that the loading angle β for pure mode II (i.e., βII) for each notch opening angle was extracted from the study carried out by Torabi and Taherkhani [23]. Additionally, the pure mode II loading angle βII depends on the notch angle 2α and, therefore, the values of βII for notch opening angles 2α = 30°, 60°, and 90° are different. The loading angles β in mixed mode loading were defined accordingly to the value of βII for each opening angle 2α. The specimens were manufactured from white marble sheets with a mean thickness of 30 mm and by means of a water jet machine. The central rhombic hole was also created by water jet machining, with the notch corners sharpened by a razor blade. Figure 2 shows the SVNBD samples with different notch angles prepared for the fracture tests. In order to obtain reliable data, four repetitions were considered for each test category. Accordingly, 48 SVNBD samples were finally tested. After manufacturing the SVNBD samples, corresponding to the dimensions given in Table 1, they were tested in a SANTAM universal test machine (SANTAM company, Tehran, Iran) with a load cell capacity of 150 kN, under displacement control conditions, and with a cross head speed of 0.5 mm/min. Figure 3a shows, as an example, a SVNBD specimen with notch angle 30 • subjected to pure mode II loading inside the test machine, while Figure 3b shows the corresponding broken halves. The fracture loads for each test sample were recorded from the test machine. The average fracture loads for each test category, as well as the standard deviations are listed in Table 1. The specimens were manufactured from white marble sheets with a mean thickness of 30 mm and by means of a water jet machine. The central rhombic hole was also created by water jet machining, with the notch corners sharpened by a razor blade. Figure 2 shows the SVNBD samples with different notch angles prepared for the fracture tests. In order to obtain reliable data, four repetitions were considered for each test category. Accordingly, 48 SVNBD samples were finally tested. After manufacturing the SVNBD samples, corresponding to the dimensions given in Table 1, they were tested in a SANTAM universal test machine (SANTAM company, Tehran, Iran) with a load cell capacity of 150 kN, under displacement control conditions, and with a cross head speed of 0.5 mm/min. Figure 3a shows, as an example, a SVNBD specimen with notch angle 30° subjected to pure mode II loading inside the test machine, while Figure 3b shows the corresponding broken halves. The fracture loads for each test sample were recorded from the test machine. The average fracture loads for each test category, as well as the standard deviations are listed in Table 1.  In addition to the SVNBD samples, Brazilian disk (BD) samples without any type of crack/notch, and cracked Brazilian disk (CBD) samples were tested in order to determine, respectively, the tensile strength and fracture toughness of white marble. Both mechanical properties were utilized in the theoretical estimations of the fracture loads. The diameter and thickness of both types of samples were, respectively, 150 mm and 30 mm. Furthermore, the crack length for all CBD samples was 2a = 75 mm. Figure 4  In addition to the SVNBD samples, Brazilian disk (BD) samples without any type of crack/notch, and cracked Brazilian disk (CBD) samples were tested in order to determine, respectively, the tensile strength and fracture toughness of white marble. Both mechanical properties were utilized in the theoretical estimations of the fracture loads. The diameter and thickness of both types of samples were, respectively, 150 mm and 30 mm. Furthermore, the crack length for all CBD samples was 2a = 75 mm. Figure 4 displays examples of BD and CBD samples before and after the tests. The average of the fracture loads for the BD samples, obtained from four test repetitions, was 28,061 N, while the average value for CBD samples (also four repetitions) was approximately 17,317 N. These two averages were the only mechanical inputs required to calculate the tensile strength and the fracture toughness, respectively. In addition to the SVNBD samples, Brazilian disk (BD) samples without any type of crack/notch, and cracked Brazilian disk (CBD) samples were tested in order to determine, respectively, the tensile strength and fracture toughness of white marble. Both mechanical properties were utilized in the theoretical estimations of the fracture loads. The diameter and thickness of both types of samples were, respectively, 150 mm and 30 mm. Furthermore, the crack length for all CBD samples was 2a = 75 mm. Figure 4 displays examples of BD and CBD samples before and after the tests. The average of the fracture loads for the BD samples, obtained from four test repetitions, was 28,061 N, while the average value for CBD samples (also four repetitions) was approximately 17,317 N. These two averages were the only mechanical inputs required to calculate the tensile strength and the fracture toughness, respectively.  The main results determined from the experiments are tensile strength (BD samples), fracture toughness (CBD samples), and critical notch stress intensity factors (SVNBD samples). The tensile strength f t can be written in terms of fracture load P f , diameter D, and thickness t of BD sample as [24]:

displays examples
By substituting the values of D = 150 mm, t = 30 mm, and P f = 28,061 N in Equation (1), the value of the tensile strength for the tested marble is determined as f t = 3.34 MPa.
To calculate the fracture toughness of the analyzed marble, Equation (2) was utilized in this study [25]: where the dimensionless parameter K * I is usually obtained from finite element analyses and depends on the crack length ratio 2a/D. The value of this parameter for CBD samples with 2a/D = 0.5 is K * I = 0.221, was extracted from the study performed by Akbardoost and Ayatollahi [25]. Substituting the dimensions of the CBD samples, the fracture load, and the dimensionless parameter K * I into Equation (2), the fracture toughness is K Ic = 0.977 MPa √ m. Finally, K SV If and K SV IIf are the critical notch stress intensity factors (NSIFs) for mode I and mode II loading conditions, extensively utilized in the literature to describe the fracture resistance of SV-notched specimens against brittle fracture. The values of NSIFs for SVNBD samples, as well as for other V-notched (rounded or sharp tip) samples, can be obtained from the stress, strain or displacements fields determined for some specific points around the notch border, which can be derived numerically (e.g., finite element method [15,16]) or experimentally (e.g., digital image correlation [13], use of strain gauge [14], photoelasticity [12], etc.). After that, a set of linear equations in which the NSIFs are unknown are obtained according to the formulations proposed for the aforementioned fields (for example, the relations derived by Williams [17] for stress components around the SV-notch tip). Solving these linear equations simultaneously allows the values of NSIFs to be obtained. Torabi and Taherkhani [23] employed the finite element method for determining the stress field around the rounded and sharp V-notched Brazilian disks (RVNBD and SVNBD samples), and calculated the corresponding mode I and mode II NSIFs for a wide range of geometrical ratios and loading conditions. For the sake of more usability, they proposed Equations (3) and (4) for SVNBD samples: Y SV I and Y SV II are dimensionless parameters depending on the relative notch length d/D, the loading angle β, and the notch angle 2α. Their values for the SVNBD samples tested in the present work were extracted from the study performed by Torabi and Taherkhani [23], and are shown in Figure 5, and listed in Table 2. λ 1 and λ 2 are eigenvalues of the stress field near the SV-notch border, and only depend on the notch angle 2α. Their values for notch angles of 30 • , 60 • , and 90 • were obtained from [26], as shown also in Table 2. Now, the values of K SV If and K SV IIf for the tested SVNBD samples were calculated by replacing the corresponding values of the average fracture load, specimen dimensions, eigenvalues λ 1 and λ 2, and dimensionless parameters Y SV I and Y SV II , as gathered in Table 2. It should be noted that the units of K SV If and K SV IIf are calculated in terms of eigenvalues λ 1 and λ 2 , which depend on the notch opening angle 2α (as easily observed from Equations (3) and (4)). Therefore, unlike the crack stress intensity factor, whose dimension is the same for both mode I and mode II SIFs, the units of NSIFs for mode I and mode II loading conditions in SV-notched samples are meaningfully different.     In Section 3, the values of K SV If and K SV IIf will be predicted by using the MTS-FEM criterion.

Analytical Approach: Mixed-Mode Fracture Criterion
Like other quasi-brittle materials, rock materials have low resistance under tensile loading and, hence, rock components often fail perpendicularly to the maximum tensile principal stress. Thus, fracture criteria based on this observation provide good estimates for the onset of fracture in rock structures. The maximum tangential stress (MTS) criterion is among these fracture criteria. In the MTS criterion, brittle fracture in cracked or notched specimens takes place when the value of the tangential stress σ θθ , at a specified critical distance from the crack or the notch tip, attains the critical value σ θθc . Also, the crack growth progresses from the defect tip perpendicularly in the direction of the maximum tangential stress. In order to apply the MTS criterion, it is necessary to determine the tangential stress component around the notch tip. Therefore, the relations for determining the stress components around the SV-notch tip are taken into consideration. A sharp V-shaped notch with a coordinate system originating at the notch tip is shown in Figure 6. Williams [17] made use of the Airy stress function method and derived the elastic stresses around a sharp V-notch. By simplifying the parameters and functions described by Ayatollahi et al. [8], the tangential stress component around the SV-notch border can be rewritten as: where K SV I and K SV I I are the notch stress intensity factors for mode I and mode II, and r and θ are the polar coordinates, respectively. Additionally: Materials 2022, 15, x FOR PEER REVIEW 9 of 17 Figure 6. A sharp V-notch with its corresponding coordinates. Table 2 above summarizes the values of the parameters λ1 and λ2 for the SVNBD specimens with different notch opening angles, as provided by Ayatollahi et al. [8]. The same authors [8] employed only the first or singular terms of Equation (5) where is the notch fracture toughness obtained from the fracture tests on SV-notched samples under pure mode I loading. Besides, the fracture initiation angle θ0 in Equation (9) is determined from the following formula derived by Ayatollahi et al. [8], derived from differentiation of the tangential stress relative to angle θ and combinations with Equations The critical distance rc,V in both Equations (9) and (10) can be determined from the material fracture toughness KIc and tensile strength ft, applying the formulation proposed by Carpinteri et al. as [27]: The previous formulation is applicable to linear-elastic brittle materials. In quasi-brittle materials such as rocks, the critical distance rc,V is large when compared to the specimen size and, therefore, on many occasions, classical fracture criteria like SV-MTS, which takes  Table 2 above summarizes the values of the parameters λ 1 and λ 2 for the SVNBD specimens with different notch opening angles, as provided by Ayatollahi et al. [8]. The same authors [8] employed only the first or singular terms of Equation (5), and proposed the SV-MTS criterion for predicting K SV I f and K SV I I f in SV-notched samples as follows: where K SV Ic is the notch fracture toughness obtained from the fracture tests on SV-notched samples under pure mode I loading. Besides, the fracture initiation angle θ 0 in Equation (9) is determined from the following formula derived by Ayatollahi et al. [8], derived from differentiation of the tangential stress relative to angle θ and combinations with Equations (3) and (4): The critical distance r c,V in both Equations (9) and (10) can be determined from the material fracture toughness K Ic and tensile strength f t , applying the formulation proposed by Carpinteri et al. as [27]: The previous formulation is applicable to linear-elastic brittle materials. In quasibrittle materials such as rocks, the critical distance r c,V is large when compared to the specimen size and, therefore, on many occasions, classical fracture criteria like SV-MTS, which takes into account only the singular terms of stress field, cannot predict accurately the fracture load of rock notched samples. In such cases, it has been suggested in previous studies that higher order terms in the stress field could be taken into consideration [28][29][30][31][32][33][34] or, alternatively, the stress or strain components could be determined directly from FEM analyses. The MTS-FEM approach, utilized recently for predicting the fracture behavior of rock notched samples [22,[35][36][37], corresponds to the second approach. It is based on the MTS criterion, in which the tangential stress at the critical distance is determined from finite element method. Assuming linear elastic behavior in the finite element analysis, the stress field ahead of the SV-notch tip is proportional to the magnitude of the applied external load. Thus, a simple proportion between the maximum value of tangential stress at the critical distance r c,V obtained from the FEM, the corresponding arbitrary load applied in the FEM, and the tensile strength f t provides the fracture load as follows: Here, it should be noted that both r c,V and θ 0 are determined from FEM. This is mainly because the objective of the approach is the prediction of the fracture resistance of SV-notched samples without computing any stress coefficients.
In this way, the cracked Brazilian disk tested for determining the fracture toughness is simulated by FEM at the corresponding fracture load conditions. Then the distance from the crack tip along the crack line to the point where the tangential stress is equal to the tensile strength is determined, and set as the critical distance r c . Additionally, the value of critical distance for the cracked and SV-notched samples is assumed to be equal in this study (i.e., it is a material property) as already reported in previous research works (e.g., [37,38]). Indeed, the critical distances for cracked and notched samples may be considered conceptually different because of the influence of both the notch root radius and the notch opening angle (e.g., [39][40][41]. However, there are studies in the literature in which the critical distance is considered the same for both notched and cracked specimens [37,38]. This simplification allows extra experimental work to be avoided and, as justified in the literature (e.g., [37,38]) and below, provides good agreement with the experimental data.
The fracture initiation angle θ 0 is determined by measuring the tangential stress at a ring with radius of r c,V around the notch tip. The direction perpendicular to the maximum tangential stress determines θ 0 . In other words, the angle between this direction and the notch bisector line is known as θ 0 . Additionally, the value of the tangential stress at θ 0 is substituted into Equation (12) as σ θθ (r c,V ,θ 0 )| FEM . In the next section, details for calculation of r c,V and σ θθ (r c,V ,θ 0 )| FEM are explained, and the corresponding results are gathered.

Results
In both SV-MTS and MTS-FEM criteria, the critical distance r c,V is a key material parameter. When applying the SV-MTS criterion, it is obtained from Equation (11), whereas when applying the MTS-FEM criterion, r c,V is derived from FEM analysis. By substituting the values of K Ic = 0.977 MPa.
√ m and f t = 3.34 MPa into Equation (11), the value of r c,V is calculated as 14 mm. Now, in order to determine the critical distance r c,V in the MTS-FEM approach, the CBD sample was modeled and meshed in Abaqus code with eight-node iso-parametric elements. Then, the singularity at the crack tip was considered by using singular collapsed elements at the first ring around the crack tip. In order to achieve reliable results, three different mesh patterns and sizes were used, and the mesh pattern shown in Figure 7 was finally generated. Furthermore, 10,000 elements were completely generated for CBD samples. The compressive load 17,317 N obtained from fracture tests was applied to the FE model. Figure 7 shows the loading conditions for a CBD sample as well as the generated mesh pattern. After simulation, the tangential stress was measured from the nodes along the crack line. The distance from the crack tip where the tangential stress is equal to 3.34 MPa (the tensile strength of tested rock) is 5 mm, which is considered as r c,V . When applying the SV-MTS criterion, after determining the value of critical distance rc,V, the fracture initiation angle θ0 is calculated by solving Equation (10). Next, the critical notch stress intensity factors (NSIFs), and , are calculated from Equation (9). The results are shown in Figure 8, where they are also compared to the experimental results. A significant discrepancy between the fracture curves obtained from the SV-MTS criterion and the experiments can be observed in the figure, except for the results of pure mode I loading.
Since the value of obtained from fracture tests is used in Equation (9), there are no differences between the theoretical prediction and the experimental result. When applying the SV-MTS criterion, after determining the value of critical distance r c,V , the fracture initiation angle θ 0 is calculated by solving Equation (10). Next, the critical notch stress intensity factors (NSIFs), K SV If and K SV IIf , are calculated from Equation (9). The results are shown in Figure 8, where they are also compared to the experimental results. A significant discrepancy between the fracture curves obtained from the SV-MTS criterion and the experiments can be observed in the figure, except for the results of pure mode I loading. Since the value of K SV Ic obtained from fracture tests is used in Equation (9), there are no differences between the theoretical prediction and the experimental result.
As mentioned earlier, calculating the tangential stress from FEM, as proposed in the MTS-FEM approach, can improve the prediction of the conventional MTS criterion for rocks. For this purpose, all SVNBD samples were modeled in Abaqus code and then discretized by nearly 15,000 eight-node iso-parametric elements. Very fine meshes were utilized around the notch tip for considering the high stress gradient in the vicinity of the notch border. It is noteworthy that, similar to the CBD sample, three mesh patterns were utilized first for attaining reliable results, and finally the mesh pattern displayed in Figure 9 was selected. An arbitrary compressive load with magnitude of 10 kN was applied along the diameter of the disk related to the notch bisector line and loading angle β. Figure 9 also shows the boundary conditions and the applied load employed in the present modeling of the SVNBD samples, in addition to the mesh pattern. The finite element analysis was performed for each model and then the values of the tangential stress component at a ring with radius of r c,V = 5 mm around the SV-notch tip were measured. The value and orientation of the maximum tangential stress were considered as σ θθ (r c,V ,θ 0 )| FEM and θ 0 , respectively. Figure 10 shows a sample with the resulting tangential stresses around the notch tip and the location of the maximum tangential stress. It is necessary to consider the coordinate transformation rules in second-order tensors due to the difference between the coordinate system used in the finite element models and the conventional coordinate system for SV-notched samples shown in Figure 6. After that, the fracture load for each SVNBD sample was computed by substituting the value of σ θθ (r c,V ,θ 0 )| FEM , tensile strength of f t = 3.34 MPa, and applied load of 10,000 N into Equation (12). By replacing the fracture loads predicted by the MTS-FEM criterion in Equations (3) and (4) the fracture curves for each tested rock sample configuration are obtained, with the results being shown in Figure 8. It can be observed that the fracture curves predicted by the MTS-FEM criterion are not only more accurate than those obtained by the SV-MTS criterion, but they are also in agreement with the experimental results.      Moreover, the fracture loads predicted by both the SV-MTS and the MTS-FEM fracture criteria are compared with those obtained in the experimental program, with the results being gathered in Table 3. The value of the fracture load in the SV-MTS criterion is determined by substituting the predicted value of K SV If or K SV IIf into Equation (3) or (4), respectively. According to the discrepancy between the theoretical and the experimental fracture loads shown in Table 3, it can be stated again that the MTS-FEM criterion is significantly more accurate than the SV-MTS criterion. In addition, the discrepancy between the fracture loads predicted by the MTS-FEM criterion and the corresponding experiments is nearly 10% for rock SVNBD samples with different notch angles. Therefore, it can be concluded that the MTS-FEM approach can be employed for predicting the onset of fracture in rock-type SV-notched parts. As an additional advantage of the MTS-FEM approach presented in this study, there is no need to calculate the stress coefficients such as NSIFs for predicting the fracture load of rock notched samples. Table 3. Comparison between the fracture load predictions in rock (marble) SVBD samples obtained from SV-MTS and MTS-FEM criteria and experimental results.  Moreover, it can be mentioned that the approach derived from the MTS-FEM criterion is independent of the loading angle β. In other words, for each loading angle (mixed mode or pure mode II), the tangential stress around the SV-notch tip is determined from FEM, and then, the fracture load is predicted.

Conclusions
The fracture resistance of rock-type sharp V-notched samples was investigated experimentally and theoretically. Several sharp V-notched Brazilian disks (SVNBD samples) were prepared from white marble sheets and then tested under different mixed mode I/II loading conditions. In order to predict the fracture resistance of the SVNBD samples, the MTS-FEM criterion was utilized. This criterion is based on the maximum tangential stress criterion in which the tangential stress around the SV-notch tip is determined from finite element analysis. Moreover, the critical distance in this approach is obtained from the finite element analysis of the cracked Brazilian disk sample under pure mode I loading. The fracture loads predicted by the present MTS-FEM criterion, as well as those predicted by the previously proposed SV-MTS criterion, were compared to those obtained in the experimental program. The results showed that the MTS-FEM criterion is not only able to predict the fracture resistance of rock-type SV-notched samples more accurately than the conventional SV-MTS criterion, but also that its predictions are in agreement with the experimental results. This improvement is due to the tangential stresses in the MTS-FEM criterion being determined from FEM, providing more accurate stresses than those derived analytically, which are used in the conventional SV-MTS criterion. Furthermore, it is not necessary to calculate the stress coefficients for predicting the fracture loads of rock notched samples, as is the case with the notch stress intensity factors in the MTS-FEM approach.