The effect of groove and notch tip angles on testing fracture toughness by SEVNB method: models and experimental validation

ABSTRACT The single-edged V-notch beam (SEVNB) method is considered as an effective method for evaluating the fracture toughness values of brittle materials. In this method, it is assumed that the V-notch is a natural crack. However, this assumption may cause an overestimation of the fracture toughness due to the “notch passivation effect”. To investigate the effects of the V-notch and groove tip angles on the fracture toughness testing of ceramic materials, three typical models were established in this work. The stress intensity factors of these models were calculated using a J-integral based on the linear finite element method (LFEM). The results indicated that the measured fracture toughness values could be overestimated by 0.5%- 13.7% when the angle of the V-notch tip increased from 10° to 60°. Increasing the angle formed by the V-notch and groove from 10° to 60°, fracture toughness was overevaluated by about 0% – 2.0%. When the angle formed by the V-notch and groove increased to 120°, the fracture toughness was overevaluated by about 31%. Finally, two equations were fitted to assess the angles effects on fracture toughness, and the results have been validated by experiments. An important reference for the SEVNB method can be found in this work.


Introduction
As an index of the mechanical properties of brittle ceramic materials, fracture toughness reflects the maximum resistance during crack growth, and it plays a vital role in the damage and destruction of ceramics [1][2][3]. Accurate testing of the fracture toughness has been continuously studied by many researchers in the past few decades. Single-edged V-notch beam (SEVNB) method shows great application potential in this field [4][5][6]. The SEVNB method was developed from the single-edged notch beam (SENB) method to improve the overestimation of the SENB method due to the "notch passivation effect" [5]. In this method, a sharp V-notch at the root of the groove is prefabricated to eliminate the effect. Nishida et al [7] studied the effect of notch root radius on the fracture toughness of Al 2 O 3 ceramic with fine grain size. It is concluded that fracture toughness can be measured when the radius of the notch tip is less than 10 μm. Unfortunately, the traditional SEVNB method still overestimates the fracture toughness values of some fine-grained ceramic materials due to notch passivation effect [6][7][8][9][10][11]. One modified method is to use an ultra-short pulsed laser to ablate a sharper V-notch at the groove root. According to the previous reports [12][13][14], the tip radius of the V-notch ablated by a laser is less than 0.5 μm, which is less than size of the microstructure feature (average grain size) of most ceramics. However, the V-notch produced by a laser is still not a natural crack, the notch tip angle may have a significant effect on the fracture toughness test, for larger angle will shield the induction and propagation of the crack. In addition, the formula used in the SEVNB method by assuming the notch is a crack may be inappropriate. Thus, the applicability of the SEVNB method should be determined.
Fett [15] first studied the SEVNB method with a model of a groove whose root was inserted in a real crack, and he recommended the following function: where K exp is the experimental stress intensity factor of the limited crack length (l) compared to the root radius of the groove root (ρ), K th is the true stress intensity factor with a crack that is sufficiently long compared to the root radius of the groove root, and Y is the appropriate geometric correction factor. Y = 1.12 for a through-thickness edge crack. Fett [10] further studied the relationship of real crack length and the groove root radius. It showed that if the groove root radius ρ > ρ c (groove critical radius), spontaneous crack extension is possible. However, the ideal crack inserted at the root of groove may be not the same as a V-notch. Recently, Wang et al. [16] studied the effect of the angle of the V-notch on fracture toughness testing with the linear finite element method (LFEM). They concluded that the fracture toughness would be overestimated when the V-notch tip angle was higher than 60°. However, they did not provide the accurate relationship between the angle and the stress intensity factor. Wang et al. [17] also studied the angle of the V-notch during fracture toughness testing by assuming the V-notch was a natural crack. In addition, the angle of the notch tip is difficult to be processed and measured accurately in three-dimensional samples. Thus, the results may be not accurate enough for fracture toughness testing. In order to analyze the impact of the angle on the fracture toughness test in the SEVNB method in more detail, the width of the groove, the tip angle and length of the V-notch should be analyzed separately to better guide the fracture toughness test.
In this work, we aimed to study the relationship between the angle of the groove tip with a V-notch and the fracture toughness by the J-integral based on the LFEM. Based on this relationship, it is anticipated that a more accurate fracture toughness can be obtained when using the SEVNB method.

Experimental procedure
Zirconia powder with 2.3 mol% Y 2 O 3 (2.3Y-TZP, Fan Meiya commercial materials Co. Ltd, Jiangxi, China), zirconia powder with 3 mol% Y 2 O 3 (3Y-TZP Fan Meiya commercial materials Co. Ltd, Jiangxi, China), and Al 2 O 3 (99.5% pure, Showa Denko, Japan) powder were used as the starting materials. Granulated raw materials were pressed into disks at 20 MPa in a steel die, and coldpressed isostatically at 250 MPa for 3 min to obtain the green bars. The green bars of the 2.3Y-TZP and 3Y-TZP were sintered in air at 1450°C for 2 h in a MoSi 2 resistance furnace (LHT 08/17, Nabertherm, Germany) at a heating rate of 5°C /min. The green bars of the Al 2 O 3 were sintered in air at 1300°C for 2 h with a heating rate of 5°C /min. The densities of the sintered disks were measured in distilled water according to Archimedes' principle, and the theoretical densities were assumed to be 6.08 g/cm 3 for the TZP ceramic and 3.98 g/cm 3 for the Al 2 O 3 ceramic. All the measured relative densities of the investigated samples were larger than 98%. The sintered samples were ground into test bars with dimensions of 3 × 4 × 35 mm. A groove was machined by a diamond wheel with thickness of 200 μm at the center of 3 × 35 mm surfaces, and sharp V-notches with different lengths were fabricated at the bottom of the groove by a femtosecond laser (Chameleon ultra; Coherent, USA) with a power and scanning speed of 70 mW and 100 μm/s, respectively. Scanning electron microscopy (SEM; ZEISSEVO 18, Oberkochen, Germany) was used to measure the V-notch length (l) on the fractured surface, while optical microscopy (SZX 10, Olympus, Japan) was used to measure the depth of the groove (u) and the groove shape. All the test bars were tested in a mechanical test machine (Model 5567, Instron, USA) using a three-point flexure test configuration with a span of 30 mm and a constant loading rate of 0.05 mm/min. The fracture toughness values were calculated according to the ASTM C-1421 standard [18], which assume the notch as a natural crack: where a is the notch length which is equal to the V-notch length (v n ) plus the depth of the groove (u); F is the fracture load; S is the span width; B is the sample width (3 mm); h is the sample height (4 mm); and Y is the dimensionless correction factor which can be calculated by the following formula: The microstructures of the sintered samples were observed by scanning electron microscopy (SEM) on the polished surfaces. The phase compositions of the samples were analyzed by X-ray diffraction (XRD, X'pert PRO, Panalytical, Netherlands). The monoclinic phase volume fractions (V m ) of the 2.3Y-TZP and 3Y-TZP were calculated using the formula proposed by Toraya et al. [19]. Figure 1 shows the SEM image of test bar with a groove and a V-notch for the SEVNB method test [12]. The tip angle of the V-notch and the width of the groove have significant effects on the stress concentration at the tip of the V-notch. When the tip angle of the V-notch is sharper and the groove is narrower, the stress at the V-notch tip becomes more concentrated. Therefore, the tip angle of the V-notch, θ 1 , was used to study the effect of the V-notch morphology on the stress intensity factor testing, and θ 2 was used to study the effect of the groove morphology and V-notch length on the stress intensity factor testing. The interaction of the laser and the material around the V-notch tip was studied by previous researchers [20][21][22]. The results showed that although a thin molten region with a width of about 2 μm was produced on the V-notch surface by a laser, the effect of the molten region on the fracture toughness testing could be ignored.

Geometric models
In order to analyze the influences of θ 1 and θ 2 on the accuracy of the fracture toughness testing by the SEVNB method in the three-point flexure test, three two-dimensional (2D) models were constructed. Figure 2a shows a schematic of the model with an ideal crack, that was used to calculate the theoretical stress intensity factor (K I/th ) of an ideal crack tip. This model is referred to as the ideal crack model. Figure 2b shows the schematic diagram of the model of a groove with a V-notch root, which was marked as a non-crack model. The groove length and V-notch length are denoted by u and v n , respectively. This model was used to calculate the stress intensity factors (K I/θ1 ) corresponding to different θ 1 values. The stress intensity factor of this model is only related to θ 1 and does not depend on θ 2 , and it can be used to guide the processing of the V-notch. Figure 2c shows the schematic of the model of a groove with a semi-circle root, and a crack was implanted in the middle of the semi-circle root. This model is marked as groove model. The groove and crack lengths are denoted as u and l. This model is used to calculate the stress intensity factors (K I/θ2 ) corresponding to different θ 2 values. The stress intensity factor of this model is only related to θ 2 and does not depend on θ 1 , and it can be used to guide the processing of the groove and V-notch length.

Stress intensity factor
In this work, the stress intensity factor was calculated using the J-integral. The J-integral is an independent path line integral along a freely selected counterclockwise contour, Г, surrounding the crack tip. The J-integral is defined as follows: where W is the strain energy density, which can be evaluated as: T is the traction vector, defined as where σ ij represents the stress components, ε ij represents the strain components, and n i represents the normal vector components. For a plane strain case and linear elastic materials, the stress intensity factor can be calculated based on the J-integral, as follows: where E is Young's modulus, and ν is the Poisson's ratio.

Convergence study
The J-integral values of the models were determined by 2D finite element (FE) analysis which was performed by Figure 1. SEM image of the groove and V-notch tilted by 30° [12]. Two angles, θ 1 , and θ 2 , were used to control the morphology of the groove and V-notch, respectively. using ABAQUS (Dassault Systemes, USA). The ideal crack model is shown in Figure 3a and its crack tip is enlarged in Figure 3b. A circumference C with the crack tip as the center and radius R is proposed to facilitate the subsequent refinement of the grid near the crack tip. The typical meshed ideal crack model is shown in Figure 3c. The refined mesh within the circumference C is shown in Figure 3d. The type of the mesh within the circumference C was a four-node bilinear plane-strain quadrilateral with reduced integration, while that outside the circumference C was four-node bilinear plane strain quadrilateral without reduced integration. As Figure 3d shows, the number of nodes (red point) on the radius R was controlled to adjust the length L between the crack tip and the first point while the number of nodes (cyan point) on the circumference C was controlled to determine the average degree (η) of the circumference C. Global elements were used to limit the element size of the entire model. The convergence study results for the mesh refinement are is shown in Table 1. As shown by item B to H, the J-integral values decreased with the mesh refinement within the circumference C. The details are shown in Figure 4. The J-integral decreased with the decrease in the angle η, and the J-integral did significantly decrease until item E. Although the angle of item E was reduced by half, the J integral was only reduced by 0.05%. Further reducing the angle by half, the J integral only increased by 0.002% (item H). Thus, the error of the J-integral calculated by the ideal crack model with the element parameters of item F was less than 0.1%. In order to study the effect of the global element size on the J-integral calculation, the smaller global element sizes (items I and J) were used to determine the J-integral values, which were less than that of item F by 0.014% and 0.026%, respectively. Therefore, in  this work, the number of nodes on the circle C, the number of nodes on the radius R, and the global element size were 100, 50, and 0.05 mm, respectively. The corresponding length L and angle η were 1.603 μm and 3.616°, respectively. The uncertainties in the calculated J-integral values were assumed to be less than 0.1%. Convergence studies of the non-crack model and the groove model were also performed, but the details are not shown here. Results show that the uncertainties in the calculated J-integral value are less than 0.1%. The number of nodes on the circumference C, the number of nodes on the radius R, and global size were also 100, 50, and 0.05 mm, respectively, for the non-crack model and the groove model. The corresponding length L an angle η were 1.603 μm and 3.616°, respectively. The model material property parameters were those reported by Quinn et al. [23], and they are listed in Table 2. Figure 5 shows the typical stress cloud (S11, along the length of the beam) images of the ideal crack model. As shown in Figure 6, the fitted LFEM results of the ideal crack model with different crack lengths were almost consistent with the ASTM C-1421 standard based on a natural crack (Equation (2)). Figure 7 shows the stress distribution along the length of the beam (S11) in the crack tip of the ideal crack model and the V-notch tip of the non-crack model. As expected, the tensile stress distribution at the crack tip was highly concentrated in the ideal crack model, and its stress contour resembled a dumbbell. Near the crack surface, an equilateral triangle area with a stress close to zero was found, which is consistent with the result reported by Wang et al. [16]. However, when the crack tip was enlarged, the two sides of the triangle near the crack tip became cambered, and closer to the crack tip, the angle (θ 1 ) was smaller. This observation was not found in the report of Wang et al. [16]. Figure 7b shows the non-crack model with a V-notch   tip angle θ 1 of 90°.The stress contour still resembled a dumbbell, although a part of the volume was removed compared to the ideal crack model. In addition, at the V-notch tip, no compressive stress area was found. Compared with the fracture load of sample with real crack, the large-angle V-notch will lead to the loss of part of the region with high elastic strain energy, so it will be unfavorable for crack induction and early propagation, and the fracture toughness of the material will be overestimated.

Effect of V-notch tip angle (θ 1 ) on fracture toughness measurement
The non-crack model can be used to assess the effect of the angle of the V-notch tip θ 1 on the fracture toughness testing by the SEVNB method.
Before the assessment, the effect of the u/a value on the J-integral of the V-notch tip must be evaluated. Figure 8 shows the J-integral of the non-crack model with different u/a values and a constant V-notch tip angle of 60°. The stress cloud images (S11, along the beam length) of the non-crack model with u/a = 0, 0.5, 0.95 are shown in Figure 8a-c, respectively. The stress distributions at the V-notch tip were similar. Moreover, the J-integral value almost remained constant as the u/a value increased. Angles of the V-notch tip θ 1 of 30°, 90°, and 120° were also studied, and the conclusions were the same as those with 60°. Thus, they are not shown in Figure 8. Therefore, the effect of the u/a value on the J-integral of the V-notch tip in the non-crack model can be ignored. When u was almost equal to a, the width of the groove became very thin similar to the results shown in Figure 8b. Thus, the J-integral of the noncrack model only depended on the angle of the V-notch tip, and it is feasible to study the effect of the θ 1 on the fracture toughness testing using the SEVNB method.
The influence of θ 1 on the fracture toughness testing based on the non-crack model with u/ a = 0 is shown in Figure 9. K Ic/th is the fracture toughness of ideal crack model. K Ic/θ1 is the fracture toughness of the non-crack model with u/a = 0. K Ic/θ1 also kept increasing as θ 1 increased. To study the effect of the notch length (a) on the fracture toughness testing, three non-crack models with a/h values of 0.25, 0.5, and 0.75 were constructed. The results are shown in Figure 9a. The three non-crack models followed the same trend with the increase in θ 1 . This implied that θ 1 had the same effect on the fracture toughness testing when a/h was in the range of 0.25-0.75. To further analyze the impact of θ 1 on the fracture toughness testing using the SEVNB method, the results obtained by the noncrack model were fitted, as shown in Figure 9b. The relationship between θ 1 and the relative fracture toughness, K Ic/θ1 /K Ic/th can be described by the following equation: K Ic/th can be calculated using Equation (5), and the test error of the fracture toughness corresponding to different θ 1 values can also be evaluated. For example, K I c/θ1 was higher than the K Ic/th by 13.71% when θ 1 was 60°. K Ic/θ1 was higher than the K Ic/th by 3.10% when θ 1 was 30°. K Ic/θ1 was just higher than K Ic/th by 0.51% whenθ 1 was 10°. This result was consistent with that of Palacios et al. [24]. In their work, the fracture toughness values of tungsten alloy samples with notches produced by a laser were studied. The results showed that the actual fracture toughness values of the    [16]. In their work, K Ic/θ1 was equal to K Ic/th when the θ 1 was less than 60°. The difference may have resulted from the assumption of perfect dumbbell-shaped stress contours. However, this assumption was not precise. As shown by the stress cloud of an ideal crack tip in Figure 10a, when θ 1 was about 30°, the stress included in the V-notch region was close to zero. However, when θ 1 reached  60°, although the stress was still close to zero in the region away from the V-notch tip, the stress near the V-notch tip reached higher levels. The stress cloud images of the non-crack model with θ 1 value of 10°, 30°, 60°, 90°, and 120° are shown in Figures 10b-f, respectively. For θ 1 value of 10° and 30°, an area with high stress concentration were reserved, and the stress in the removed area was close to zero when comparing with the stress cloud of natural crack (Figure 10a). However, the area with high stress concentration near the crack tip was eliminated when θ 1 was greater than 60°. The fracture toughness of the V-notch with the tip at 60° may be overestimated due to the removal of the high stress area. This observation was basically consistent with the results of the simulation (Equation (5)).

Effect of groove tip angle (θ 2 ) on fracture toughness measurement
To study the influence of θ 2 on the fracture toughness test, a groove model that was not affected by θ 1 was established. In most cases, the V-notch is assumed to be a crack when testing the fracture toughness by the SEVNB method. Fett, et al [15]. determined that the actual fracture toughness can be obtained when the length of the V-notch is longer than the sum of threegrain size scales of the sample materials. However, the fracture load increases due to the geometric shielding effect at the V-notch tip, leading to an overestimation of fracture toughness. Thus, different widths of the groove may need different lengths of the V-notch [5,15]. In this analysis, grooves with different groove widths were built, as shown in Figure 2c. As shown in Figure 11, although the groove widths of 100 μm, 200 μm, and 500 μm were different, θ 2 had the same effect on the fracture toughness testing results. In addition, different groove lengths were also studied. As shown in Figure 11a, the relative length a/h (0.2-0.6) had the same effect on the fracture toughness testing. To  study the effect of θ 2 on the fracture toughness testing, the data in Figure 11a were fitted, and the result is shown in Figure 9b and can be expressed as follows: (6) Equationn (6) shows that K Ic /θ2 is higher than K Ic/th by 0.5% when θ 2 is 40°, and K Ic/θ2 is higher than K Ic/th by 2.0% when θ 2 is 60°, while K Ic/θ2 is higher than K Ic/th by 31.1% when θ 2 is 120°. Equation (1) predicts that K Ic/θ2 is higher than the K Ic/th by 2.20% when θ 2 is 60° (l = ρ), while K Ic/θ2 is higher than K Ic/th by 0.40% whenθ 2 is 40° [15]. These results were almost consistent with the results obtained by LFEM in this work. In addition, these results indicated that θ 1 had more influence on the fracture toughness testing than θ 2 . The stress cloud images of the groove with θ 2 values of 0, 19°, 60°, 91°, 113°, and 131° are shown in Figure 12. Regardless of the value of θ 2 of the crack tip, the contour of the stress cloud images at the crack tip remained almost the same. Meanwhile, the stress distribution almost remained constant due to the large distance between the crack tip and the groove tip. This is the reason that θ 1 had a greater impact on the fracture toughness test than θ 2 . The stress did not accumulate significantly for the groove tip far from the crack tip. Thus, θ 2 has a weaker effect on the fracture toughness testing when using the SEVNB method. It seems that θ 2  has a weaker effect than θ 1 on the fracture toughness test. Although the width of the groove is usually less than 1 mm, the length of V-notch is generally 100 μm (the process of longer V-notch will lead to a decrease in the sharpness of the V-notch). Thus, a relatively narrow groove may lead to a larger θ 2 due to the too small length of the V-notch. Even the V-notch is sharp enough, the large θ 2 may shield the induction and propagation of the crack at the tip of the notch. Therefore, Both the length of V-notch and the width of the groove affect the stress concentration at the tip. Figure 13 shows the typical thermally etched microstructures of the investigated samples. The grains in all the samples had fairly equiaxed shapes and the average grain sizes were found to be about 500 nm, 550 nm and 1200 nm for 2.3Y-TZP, 3Y-TZP, and Al 2 O 3 , respectively. No residual pores were detected. These three ceramic materials have a dense and uniform microstructure, which is suitable for the research on linear elastic fracture of materials.   [29], and (c) 3Y-TZP [28] samples with grooves and V-notches.

Experimental results and Discussion
The relationship between the fracture toughness test values of the three ceramic materials and the angle is shown in Figure 14. This indicates that the fracture toughness tested by the SEVNB method increases as θ 2 increases. In addition, the results predicted using Equation (6) are plotted in Figure 14 to compare with the experimental results. It is clear that the predicted curve fit well with the experimental values for the Al 2 O 3 sample, while the fracture toughness values were slightly different for 2.3Y-TZP and 3Y-TZP in the highangle area. The difference can be explained by the different R-curves of the three samples. As shown in Figure 15, the volume fractions of the m-ZrO 2 on the polished surface were zero, and the phase transformations on the fractured surface were 10 vol%, and 30 vol % for 2.3Y-TZP, and 3Y-TZP, respectively. Therefore, R-curve was steep to prevent the crack growth in the TZP samples, in contrast to the R-curves of Al 2 O 3 samples reported in the literature [25][26][27]. The cracks in the 2.3Y-TZP, and 3Y-TZP samples grew stably before fracture. Thus, the actual θ 2 before fracture of the samples became larger. As shown in Figure 14 the predicted curves of 2.3Y-TZP and 3Y-TZP had light differences from the experimental data, but the predictions were the same as the experimental data in the Al 2 O 3 samples for larger θ 2 value. In addition, the fracture toughness values of the 3Y-TZP and Al 2 O 3 samples were measured by the SEVNB method in our previous work [28][29][30]. The morphologies of the V-notch and groove for 3Y-TZP and Al 2 O 3 samples were also studied, θ 1 was in the range of 20°-30°, while θ 2 was about 105° (Figure 16). The fracture toughness was compared with the predicted result of Equation (6) in Figure 14. Although the fracture toughness reported in our previous work was overestimated because of the larger θ 2 , the relationship between the fracture toughness and θ 2 was consistent with the results of this work, as shown in Figures 14b  and 14c.
As listed in Table 3, the fracture toughness values of the 2.3Y-TZP, 3Y-TZP, and Al 2 O 3 predicted by Equation (6)    Wang et al. [17] studied the effect of the notch tip angle on the fracture toughness testing through experiments with typical brittle ZrB 2 and ZrB 2 -SiC ceramics. When the notch tip angle was small, they used samples with only θ 1 (typical samples are shown in Figure 17b,c). When the notch tip angle was large, they used samples with θ 1 and θ 2 (typical samples are shown in Figure 17d,e), and ignored the effect of θ 1 on the fracture toughness testing. In addition, they used a non-crack model to compare with the experimental data. The results are shown in Figure 16a. The match between the experimental data and the prediction curve was poor when the notch tip angle was large. This deviation mainly arose from ignoring the influence of θ 1 when the angle was large. As shown in Figure 17d, when θ 2 was 50°, θ 1 was only 9° which leads a limited effect for the fracture toughness values. However, as shown in Figure 16h, when θ 2 was 140°, θ 1 reached 21° which should not be ignored. The fracture toughness of the ZrB 2 sample with θ 1 and θ 2 values of 21° and 140° predicted by its non-crack model was 4.5 MPa•m 1/2 , but those predicted by Equation 5 (noncrack model) and Equation 6 (groove model) were both 4.7 MPa•m 1/2 . The result is marked in Figure 17a. The fracture toughness values predicted by Equations 5 and 6 were closer to the average of the experiment values. The experiment data and the predictions by Equations 5 and 6 are listed in Table 4. The fracture toughness (K IC/ref ) can be obtained from the SEPB method, which can be considered to be the ideal crack model. For the non-crack model in the ZrB 2 and ZrB 2 -SiC ceramics, the rate of the fracture toughness (K IC/SEVNB /K IC/Ref ) obtained from the experiment were basically consistent with the prediction by Equationn 5 (K IC/SEVNB /K IC/th ). In addition, in tungsten alloys, Equation 5 was also verified to yield good predictions [24]. When θ 1 was 50°, the fracture toughness of the test was overevaluated by 20%, and that predicted by Equation 5 was overevaluated by 21%. When θ 1 was 3°, the experimental test could obtain the real fracture toughness, which was also consistent with that predicted by Equation 5. For the ZrB 2 material with the mixed model notch (non-crack model and groove model), when θ 1 and θ 2 were 9° and 50°, respectively, the test predicted that the fracture toughness would be overestimated by about 2%, and the fracture toughness values predicted by Equations 5 and 6 would be overestimated by 1.5%. When θ 1 and θ 2 were 21° and 140°, the overestimation was 77.2%, as predicted by Equations 5 and 6, respectively, which was almost consistent with the average value of the test data as shown in Figure 17a. In general, different notch shapes may be used to assess the fracture toughness by the SEVNB method. Both Equations 5 and 6 can be used to analyze the excessive fracture toughness values of samples with only V-notches and samples with mixed notches (V-notch + groove). This can be served as a reference for testing the real fracture toughness values of brittle materials.

Conclusion
In this work, two important angles in the fracture toughness testing of ceramic materials by the SEVNB method were proposed to determine the shape of the prefabricated notch and groove. The effects of the two angles on the measured fracture toughness were analyzed, and quantitative relationships (Equation (5) and (6)) were obtained by LFEM and validated by experiments. The V-notch tip angle (θ 1 ) had a significant influence on the tested fracture toughness, which was overevaluated by 13.7%, when the angle of the V-notch was 60°. When the angle formed by the V-notch and groove (θ 2 ) was 60°, the fracture toughness was only overevaluated by 2.0%. The fracture toughness was overevaluated by about 31% when θ 2 increased to 120°. As a result of this research, we can more accurately and reliably test the fracture toughness values of ceramic materials using the SEVNB method.