Numerical Simulation of the Influence ofWidth of a Prefabricated Crack on the Dimensionless Stress Intensity Factor of Notched Semi-Circular Bend Specimens

To analyze the effect of the width of a prefabricated crack on the dimensionless stress intensity factor of notched semi-circular bend (NSCB) specimens, ABAQUS software was employed to perform numerical calibration of the crack tip stress intensity factor for the width of prefabricated cracks in the range of 0.0∼2.0mm.)e relative errors of the dimensionless stress intensity factor for different widths of prefabricated cracks were analyzed. )e results indicate that the dimensionless stress intensity factor shows an approximate linear increase as the width of the prefabricated crack increases. )e longer is the length of the prefabricated crack, the “faster” is the increase in speed. )e effect of the dimensionless support spacing on the increase in the speed of the dimensionless stress intensity factor due to the increase in crack width is minimal. When the prefabricated crack width is 2.0mm, the maximum relative error of the dimensionless stress intensity factor is 4.325%. )e new formula for the dimensionless stress intensity factor that eliminates the influence of the width of a prefabricated crack is given, which provides a theoretical basis for the more accurate fracture toughness value measured using an NSCB specimen.


Introduction
Rock mass is a typical noncontinuous natural medium, whose interior presents many defects, such as joints and cracks [1,2].Because rock fractures cannot be simply determined by strength parameter measurement, the existence of internal cracks should be considered.Fracture toughness can characterize the ability of rock fissures to resist crack initiation and propagation [3][4][5].Regarding the testing method, the International Society for Rock Mechanics (ISRM) has successively presented chevron notched short rod (SR) specimen, chevron notched three-point bend round bar (CB) specimen, cracked chevron notched Brazilian disc (CCNBD) specimen, and notched semi-circular bend (NSCB) specimen-a total of four types of specimens to test the mode I fracture toughness of rock [6][7][8].Besides the four test configurations recommended by the ISRM, some scholars have proposed some specimens of other configurations to test the pure mode I, II, and III or mixed mode fracture toughness of rock materials.Such as diametrally compressed ring (DCR) specimen [9], cracked straight-through flattened Brazilian disk (CSTFBD) [10], edge cracked triangular (ECT) specimen [11], and edgenotched disc bend (ENDB) specimen [12][13][14].e NSCB specimen is a new type of specimen employed in static fracture toughness tests, which was recommended by ISRM in 2014. is specimen was proposed by Chong and Kuruppu [15] and has been extensively employed in fracture toughness tests due to its simple structure, ease of specimen preparation, and loading [16][17][18][19][20][21][22][23][24].However, in the calculation of fracture toughness by most researchers, the dimensionless stress intensity factor was obtained by substitution in the relevant formula recommended by the ISRM.In the process of deriving this formula, the recommended method applied ABAQUS software to establish a two-dimensional NSCB specimen model that contains line cracks (i.e., zero-width cracks) to calibrate the dimensionless stress intensity factor and then obtain a relevant calculation formula by fitting.During the preparation of a fracture toughness test specimen, which is affected by factors such as preparation methods, an artificial prefabricated crack has a certain width.Its crack tip is not ideally sharp but approximates the arc shape.e arc radius may substantially vary due to different processors, which causes a non-strict correspondence relation between the dimensionless stress intensity factor given by the recommended formula substitution and fracture load obtained by experiments.e obtained fracture toughness exhibits a certain deviation from the true fracture toughness value.Currently, the study which regards the influence of width of the prefabricated crack on the stress intensity factor is primarily focused on a disc-type specimen or three-point bending beams.Zhu and Wang [25] performed detailed experimental research and theoretical analysis of the relationship between the curvature radius of a rock specimen for three-point bending tests and the test value of fracture toughness.Dai and Wang [26] analyzed the influence of crack width on the dimensionless stress intensity factor of CCNBD specimens.Dong et al. [27] investigated the effect of crack form on the stress intensity factor of a cracked straight-through Brazilian disc (CSTBD).Zhang et al. [28] and Zhang and Liang [29] successively investigated the influence of the width of a prefabricated crack on the fracture toughness test of a holed-cracked flattened Brazilian disc (HCFBD) specimen.Zhou et al. [30] compared the fracture behavior and toughness of straight-notch disc specimens with those of central sharpnotch disc specimens.
ey discovered that the fracture toughness test values of straight-notch specimens were excessive, the crack usually initiated at the corner of the notch, and the test results were more scattered.Cui et al. [31] used a Hopkinson pressure bar to conduct dynamic impact tests of non-ideal crack disc specimens with center incisions.
e test results showed that the ideal disc specimen can be replaced with a non-ideal crack disc specimen with a center incision when the width of the prefabricated crack is less than 1 mm.Wang and Luo et al. [32][33][34] proposed a method that simultaneously determines the fracture toughness and tensile strength of rock using a series of U-shaped notched beams with the same notch depth and different root curvatures.e feasibility of this method was verified by both theoretical analysis and experimentation.Kolhe et al. [35] explored the influence of the crack tip radius on the fracture toughness of CB specimens.Barati and Alizadeh [36] performed research on the ratio of the notch tip curvature to the notch depth and its influence on the plane strain fracture toughness.Cicero et al. [37] examined the influence of the three-point bending beam notch radius on the dimensionless stress intensity factor and fracture toughness by the theory of critical distances (TCD) and explicated the influencing mechanism of the notch tip radius by scanning electron microscopy.Dehghany et al. [38] investigated the effects of the first non-singular stress terms on the fracture assessment of sharp V-notches under mixed mode loading.Hussain and Murthy [39] have proposed a simple, robust, and efficient point substitution type displacement based technique for finite element estimation of the notch stress intensity factors (NSIFs) of sharp V-notched configurations.
e technique can acquire accurate NSIFs even in course meshes made of quadratic elements without the use of any special singular elements at the notch tip.
Currently, the influence of the width of a prefabricated crack on the dimensionless stress intensity factor of the NSCB specimen has not been reported.Analysis of the extent to which the width of the prefabricated crack influences the dimensionless stress intensity factor is necessary to improve the method that measures the fracture toughness of NSCB specimens.
For this reason, the dimensionless stress intensity factor of the crack tip was calibrated by ABAQUS software for different widths of the prefabricated cracks of the NSCB specimens in this paper.
e influence of the width of a prefabricated crack on the dimensionless stress intensity factor is examined; the relative error of the dimensionless stress intensity factor for different widths of the prefabricated cracks is analyzed; and the corrected formula of the NSCB dimensionless stress intensity factor that eliminates the effect of crack width is obtained.It is worth pointing out that it is not easy to prefabricate a crack on rock, and the crack tip of the manually prepared specimen is not an ideal sharp type, but a U-shaped crack with a certain radian.erefore, the next numerical research work in this paper is to treat the prefabricated crack into a U-shaped crack, which is consistent with the actual crack form in the specimen used by researchers to test rock fracture toughness.

Calibration of Dimensionless Stress Intensity Factor
2.1.Fracture Toughness Principle of NSCB Specimen.
Figure 1 shows the loading of an NSCB specimen.e specimen in Figure 1 carries the concentrated load P; its radius is R; its disc thickness is B; the length and width of the artificially prefabricated crack are a and 2b, respectively; and the support spacing is S. Considering the influence of the width of a prefabricated crack, the calculation of the mode I stress intensity factor is shown as follows: where P is the concentrated load, α 0 � a/R is the dimensionless length of the prefabricated crack, α 2b � 2b/R is the dimensionless width of the prefabricated crack, α t � B/R is the dimensionless specimen thickness, and α S � S/2R is the dimensionless support spacing.e dimensionless stress intensity factor value Y(α 0 , α 2b , α t , α S ) that is associated with length, width, and thickness of the prefabricated crack and the support spacing of the specimen is calculated by the following formula: (2)

Verification of Calibration Accuracy.
To verify the accuracy of the equivalent regional integration method with respect to the stress intensity factor calibration in this paper, 2 Shock and Vibration the geometric parameters that were employed in the method suggested by the ISRM were selected.In the case of the dimensionless support spacing α S � 0.5, the stress intensity factors of the specimens were individually calibrated for different dimensionless crack lengths α 0 � 0.1, 0.2, 0.3, 0.4, 0.5, 0.6, 0.7, and 0.8.e calibrated result is compared with that derived from the internationally recommended formula (formula (3)); the result is shown in Figure 2. ( Figure 2 reveals that the calibrated result obtained by the equivalent region integration method of this paper is similar to that derived from the internationally recommended method formula (formula (3)).When the dimensionless support spacing is 0.5 and the crack length is 0.8, the maximum error is 4.4%.However, when the dimensionless crack length ranges from 0.4 to 0.6, the maximum relative error of dimensionless stress intensity factor is 0.62%.Within this range, the numerical calibration method adopted in this paper has good accuracy.erefore, the dimensionless crack length studied in this paper ranges from 0.4 to 0.6.

NSCB-Related Parameters.
According to the range of support spacing and crack length given by the internationally recommended method, the dimensionless stress intensity factor with the common width of a prefabricated crack (0.0∼2.0 mm) is calibrated for the NSCB specimens.e geometric and mechanical parameters of the specimens are shown in Table 1.e material parameters used in the numerical model in this paper are based on the real material parameters of the limestone from Jiaozuo City, Henan Province, China.e geometrical parameter range covers the geometrical dimensions of NSCB specimens used in most cases for testing fracture toughness.

NSCB Specimen Finite Element Model.
e NSCB specimens with prefabricated crack widths of 0.0∼2.0mm, dimensionless lengths of the prefabricated cracks of 0.4∼0.6, and a dimensionless support spacing of 0.5∼0.8 were analyzed using ABAQUS finite element software.
e NSCB specimen loading is shown in Figure 1; the finite element model is shown in Figure 3.
Since the degree of the density of the mesh near the tip of the prefabricated crack will have a substantial influence on the calibration result of the dimensionless stress intensity factor, the prefabricated crack tip for different widths in Figure 3 is arranged using a local seeding technique to place seeds with the same meshing density, as shown in Figure 4(b).Considering the singularity of the crack tip, the 6-node triangular plane stress element (CPS6) is selected for the element in the crack tip integral region, and the 8-node quadrilateral element (CPS8) is selected for the grid element outside the crack tip integral region.e meshing technique uses a combination of free and sweeping meshing techniques.In terms of the boundary conditions, a load is applied to the top of the semidisc specimen, and a vertical upward constraint is applied to the left and right sides of the bottom of the semidisc with a centrally symmetric position.

Calibration Results and Pattern Analysis of the Dimensionless Stress Intensity
Factor.Based on the equivalent region integration method, the J-integral at the crack tip is obtained by ABAQUS software, and the stress intensity factor of the mode I crack is calculated by the following formula: where E is the elastic modulus.Substituting formula (4) into formula (2), the dimensionless stress intensity factor of the crack tip of the NSCB specimen with different crack widths is obtained for different support spacings and lengths of prefabricated cracks.e following calibration results are shown in Figure 5.
(1) In the case of a fixed support spacing and crack length, when the influence of crack width is considered (i.e., 2b > 0 mm), the calibration value of the dimensionless stress intensity factor is larger than that of the ideal crack (i.e., 2b � 0 mm).With an increase in the width of the prefabricated crack, the dimensionless stress intensity factor continuously increases.
(2) When the support spacing is fixed, the length of the prefabricated crack is within the range α 0 � 0.4∼0.6; the dimensionless stress intensity factor of each width of the prefabricated crack increases with the length of the prefabricated crack; and the longer is the length, the "faster" is the increase in speed.(3) When the length of the prefabricated crack is fixed, the dimensionless stress intensity factor of the NSCB specimen increases with an increase in the width of the prefabricated crack within the range of the support spacing α S � 0.5∼0.8; it is almost independent of the support spacing.

Analysis of Error Induced by Ignoring the Width of the Prefabricated Crack.
e relative error δ between a prefabricated crack of a certain width and the ideal crack is calculated by formula (5).Table 2 shows the relative error between the dimensionless stress intensity factor and the ideal crack for different support spacings and lengths of prefabricated cracks: Figure 6 shows the maximum relative error bar chart of the dimensionless stress intensity factor of the NSCB specimen for the condition that the width of the  prefabricated crack is fixed at 2.0 mm, the dimensionless length of the prefabricated crack ranges from 0.4 to 0.6, and the dimensionless support spacing ranges from 0.5 to 0.8.As shown in Figure 6, the increase in the width of the prefabricated crack significantly affects the dimensionless stress intensity factor of the NSCB specimen, which is simultaneously affected by the support spacing and crack length.e maximum relative error of the dimensionless stress intensity factor (4.325%) is obtained at the support spacing of α S � 0.7 and crack length of α 0 � 0.6.When the dimensionless support spacing is 0.5 and 0.8, the relative error of the dimensionless stress intensity factor continuously increases with an increase in the length of the prefabricated crack.When the dimensionless support spacing is 0.6 and 0.7, the change pattern does not monotonously increase.

Dimensionless Stress Intensity Factor Formula Eliminates the Influence of Width of the Prefabricated Crack.
To obtain a more general situation, the crack width is non-dimensional to obtain a more general relationship between the dimensionless stress intensity factor and the dimensionless crack width.
e three-dimensional scatter plot and the fitted surface are shown in Figure 7.
As shown in Figure 7, the scatter plot between the dimensionless stress intensity factor and the crack length and crack width is distributed in the vicinity of a curved surface

Shock and Vibration
in the range of the support spacing of 0.5∼0.8.erefore, a curved surface formula is used to characterize their relationship.e fitted surface formula is shown in the following formula: where k 0 ∼k 5 are parameters that represent a curved surface, which is related to the support spacing, Y represents the dimensionless stress intensity factor, and α 2b and α 0 represent the dimensionless length and dimensionless width of the crack, respectively.e values of the parameters in formula (6) and the correlation of the surface formulas for different support spacing are shown in Table 3.
When the support spacing ranges from 0.5 to 0.8, the correlation between the fitted surface by formula (6) and the data exceeds 0.99, which indicates that Formula (6) adequately reflects the relationship between the dimensionless stress intensity factor and the dimensionless length and width of the crack.After the preparation of the NSCB specimen, only the length and width of the prefabricated crack should be measured and appropriate parameter values for the formula should be selected according to the support spacing selected during the loading of the specimen to consider the influence of crack width and calculate the corresponding dimensionless stress intensity factor.is method eliminates the influence of the width of a prefabricated crack on the dimensionless stress intensity factor.
In theory, accurate rock fracture toughness test values can be obtained by the correction calculation of Formula ( 6) within the investigated range of the prefabricated crack size.However, the NSCB specimens of different sizes and widths of the prefabricated cracks have dimension effects and specimen heterogeneity effects during testing.

Discussion
Fracture toughness is one of the most important parameters in fracture mechanics.Accurate measurement of fracture toughness is of great significance for safety protection of engineering structures.e test results of fracture toughness are affected not only by the precision of test equipment, sample preparation, and human factors but also by the dimensionless stress intensity factor.erefore, accurate calibration of dimensionless stress intensity factors of corresponding samples is a prerequisite for accurate measurement of fracture toughness.
Generally speaking, the research on dimensionless stress intensity factors mainly focuses on the geometry of the specimen, the length of the crack, and the angle between the crack plane and the load.However, when processing the specimen to test the fracture toughness, the prefabricated crack will inevitably have a certain width, which will make the stress distribution of blunt notched crack different from that of ideal zero width.Ayatollahi et al. [40] proposed an "effective Stretch Criterion" based on the existence of a damage zone around the crack tip and the dominant uniaxial nature of the stress field near the crack for fracture assessment of rubber-like materials weakened by a mode I crack.Berto and Ayatollahi [41] based on the local energy theory to assess the fracture properties of Brazil disk with blunt crack.Aliha et al. [42]  6 Shock and Vibration successfully applied average strain energy density (ASED) criterion to composite fracture testing of rocks.On the other hand, an assumption of two-dimensional plane stress or plane strain delivers useful two-dimensional results with reasonable accuracy.However, threedimensional analysis is needed to fully understand the fracture process.Aliha and Saghafi [43] investigated the effects of thickness and Poisson's ratio on 3D mixed mode fracture.ey found the stress intensity factor and T-stress results obtained from two-dimensional conditions, present lower estimate data for the real threedimensional cracked semi-circular bend specimen.Also, the influence of specimen thickness and Poisson's ratio is more pronounced for dominantly mode I condition in the semi-circular bend specimen.Berto et al. [44] made a review of corner point singularities under in-plane and out-of-plane loading.What is more, in the literature [13], Aliha et al. calibrated the dimensionless stress intensity factor (SIF) of the ENDB specimen comprehensively.In the next research step, the stress intensity factor of the three-dimensional U-shaped crack under mixed mode loading can be studied.

Conclusions
e influence of the width of a prefabricated crack on the dimensionless stress intensity factor of NSCB specimens was analyzed using ABAQUS numerical software.e following main conclusions are obtained: (1) As the width of the prefabricated crack increases, the dimensionless stress intensity factor of the NSCB specimen continuously increases.e longer is the length of the prefabricated crack, the "faster" is the increase in speed.e dimensionless support spacing has a minimal effect on the speed of change in the dimensionless stress intensity factor as the width of the prefabricated crack increases.(2) Based on the dimensionless stress intensity factor obtained by calibration using the ideal crack, the errors in the dimensionless stress intensity factor caused by different widths of the prefabricated cracks are obtained when the widths are less than or equal to 2.0 mm.When the support spacing is α S � 0.7 and the length of the prefabricated crack is α 0 � 0.6, the maximum relative error is 4.325%.erefore, the influence of the artificial width of the prefabricated crack on the dimensionless stress intensity factor cannot be disregarded.
(3) e corrected formula that considers the influence of the width of a prefabricated crack on the dimensionless strength factor is given for any dimensionless support spacing in the range of 0.5∼0.8. e degree of fit using the formula exceeds 0.99.Using this formula and the specimen fracture load and geometric parameters, the fracture toughness test value of an NSCB specimen considering the influence of the width of a prefabricated crack can be obtained.

Nomenclature a:
Crack length b: Half of crack width B: ickness of specimen E: Elastic modulus J: J-integral K I : Mode I stress intensity factor k 0 ∼k 5 : Fitting constants related to formula (6) P: Concentrated load P max : Crack initiation load R: Radius of specimen S: Support spacing μ: Poisson's ratio Y: Dimensionless stress intensity factor Y 0 : Dimensionless stress intensity factor of the 0 mm-width crack Y 2b : Dimensionless stress intensity factor of the 2b mm-width crack α 0 : Dimensionless crack length α 2b : Dimensionless crack width α t : Dimensionless specimen thickness α s : Dimensionless support spacing δ: Relative error CB: Chevron notched three-point bend round bar specimen CCNBD: Cracked chevron notched Brazilian disc specimen CPS6: 6-node triangular plane stress element CPS8: 8-node quadrilateral element CSTFBD: Cracked straight-through flattened Brazilian disk DCR: Diametrally compressed ring ECT: Edge cracked triangular ENDB: Edge-notched disc bend HCFBD: Holed-cracked flattened Brazilian disc specimen ISRM: International Society for Rock Mechanics SR: Chevron notched short rod specimen TCD: eory of critical distances.

Figure 3 :Figure 4 :
Figure 3: NSCB models with different widths of the prefabricated cracks.

Table 2 :
Errors of dimensionless stress intensity factors of the NSCB specimens due to different widths of prefabricated cracks.