Research Article  Open Access
Numerical Simulation of the Influence of Width of a Prefabricated Crack on the Dimensionless Stress Intensity Factor of Notched SemiCircular Bend Specimens
Abstract
To analyze the effect of the width of a prefabricated crack on the dimensionless stress intensity factor of notched semicircular 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.0 mm. The relative errors of the dimensionless stress intensity factor for different widths of prefabricated cracks were analyzed. The results indicate that the dimensionless stress intensity factor shows an approximate linear increase as the width of the prefabricated crack increases. The longer is the length of the prefabricated crack, the “faster” is the increase in speed. The 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.0 mm, the maximum relative error of the dimensionless stress intensity factor is 4.325%. The 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.
1. 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–5]. Regarding the testing method, the International Society for Rock Mechanics (ISRM) has successively presented chevron notched short rod (SR) specimen, chevron notched threepoint bend round bar (CB) specimen, cracked chevron notched Brazilian disc (CCNBD) specimen, and notched semicircular bend (NSCB) specimen—a total of four types of specimens to test the mode I fracture toughness of rock [6–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 straightthrough flattened Brazilian disk (CSTFBD) [10], edge cracked triangular (ECT) specimen [11], and edgenotched disc bend (ENDB) specimen [12–14]. The NSCB specimen is a new type of specimen employed in static fracture toughness tests, which was recommended by ISRM in 2014. This 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–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 twodimensional NSCB specimen model that contains line cracks (i.e., zerowidth 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. The arc radius may substantially vary due to different processors, which causes a nonstrict correspondence relation between the dimensionless stress intensity factor given by the recommended formula substitution and fracture load obtained by experiments. The 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 disctype specimen or threepoint 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 threepoint 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 straightthrough 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 holedcracked flattened Brazilian disc (HCFBD) specimen. Zhou et al. [30] compared the fracture behavior and toughness of straightnotch disc specimens with those of central sharpnotch disc specimens. They discovered that the fracture toughness test values of straightnotch 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 nonideal crack disc specimens with center incisions. The test results showed that the ideal disc specimen can be replaced with a nonideal 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–34] proposed a method that simultaneously determines the fracture toughness and tensile strength of rock using a series of Ushaped notched beams with the same notch depth and different root curvatures. The 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 threepoint 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 nonsingular stress terms on the fracture assessment of sharp Vnotches 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 Vnotched configurations. The 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. The 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 Ushaped crack with a certain radian. Therefore, the next numerical research work in this paper is to treat the prefabricated crack into a Ushaped crack, which is consistent with the actual crack form in the specimen used by researchers to test rock fracture toughness.
2. Calibration of Dimensionless Stress Intensity Factor
2.1. Fracture Toughness Principle of NSCB Specimen
Figure 1 shows the loading of an NSCB specimen. The 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 is the concentrated load, is the dimensionless length of the prefabricated crack, is the dimensionless width of the prefabricated crack, is the dimensionless specimen thickness, and is the dimensionless support spacing. The dimensionless stress intensity factor value 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:
(a)
(b)
2.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, the geometric parameters that were employed in the method suggested by the ISRM were selected. In the case of the dimensionless support spacing , the stress intensity factors of the specimens were individually calibrated for different dimensionless crack lengths = 0.1, 0.2, 0.3, 0.4, 0.5, 0.6, 0.7, and 0.8. The 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. Therefore, the dimensionless crack length studied in this paper ranges from 0.4 to 0.6.
3. Finite Element Analysis
3.1. NSCBRelated 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. The geometric and mechanical parameters of the specimens are shown in Table 1. The 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. The geometrical parameter range covers the geometrical dimensions of NSCB specimens used in most cases for testing fracture toughness.

3.2. NSCB Specimen Finite Element Model
The NSCB specimens with prefabricated crack widths of 0.0∼2.0 mm, 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. The 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 6node triangular plane stress element (CPS6) is selected for the element in the crack tip integral region, and the 8node quadrilateral element (CPS8) is selected for the grid element outside the crack tip integral region. The 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.
3.3. Calibration Results and Pattern Analysis of the Dimensionless Stress Intensity Factor
Based on the equivalent region integration method, the Jintegral 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. The 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 ; 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 ; it is almost independent of the support spacing.
(a)
(b)
(c)
(d)
3.4. Analysis of Error Induced by Ignoring the Width of the Prefabricated Crack
The 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. The maximum relative error of the dimensionless stress intensity factor (4.325%) is obtained at the support spacing of and crack length of . 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.
3.5. Dimensionless Stress Intensity Factor Formula Eliminates the Influence of Width of the Prefabricated Crack
To obtain a more general situation, the crack width is nondimensional to obtain a more general relationship between the dimensionless stress intensity factor and the dimensionless crack width. The threedimensional scatter plot and the fitted surface are shown in Figure 7.
(a)
(b)
(c)
(d)
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 in the range of the support spacing of 0.5∼0.8. Therefore, a curved surface formula is used to characterize their relationship. The 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 and represent the dimensionless length and dimensionless width of the crack, respectively. The 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. This 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.
4. 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. The 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. Therefore, 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 rubberlike 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] successfully applied average strain energy density (ASED) criterion to composite fracture testing of rocks. On the other hand, an assumption of twodimensional plane stress or plane strain delivers useful twodimensional 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. They found the stress intensity factor and Tstress results obtained from twodimensional conditions, present lower estimate data for the real threedimensional cracked semicircular bend specimen. Also, the influence of specimen thickness and Poisson’s ratio is more pronounced for dominantly mode I condition in the semicircular bend specimen. Berto et al. [44] made a review of corner point singularities under inplane and outofplane 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 threedimensional Ushaped crack under mixed mode loading can be studied.
5. Conclusions
The influence of the width of a prefabricated crack on the dimensionless stress intensity factor of NSCB specimens was analyzed using ABAQUS numerical software. The 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. The longer is the length of the prefabricated crack, the “faster” is the increase in speed. The 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 and the length of the prefabricated crack is , the maximum relative error is 4.325%. Therefore, the influence of the artificial width of the prefabricated crack on the dimensionless stress intensity factor cannot be disregarded.(3)The 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. The 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:  Thickness of specimen 
E:  Elastic modulus 
J:  Jintegral 
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 mmwidth crack 
Y_{2b}:  Dimensionless stress intensity factor of the 2b mmwidth crack 
α_{0}:  Dimensionless crack length 
α_{2b}:  Dimensionless crack width 
α_{t}:  Dimensionless specimen thickness 
α_{s}:  Dimensionless support spacing 
δ:  Relative error 
CB:  Chevron notched threepoint bend round bar specimen 
CCNBD:  Cracked chevron notched Brazilian disc specimen 
CPS6:  6node triangular plane stress element 
CPS8:  8node quadrilateral element 
CSTFBD:  Cracked straightthrough flattened Brazilian disk 
DCR:  Diametrally compressed ring 
ECT:  Edge cracked triangular 
ENDB:  Edgenotched disc bend 
HCFBD:  Holedcracked flattened Brazilian disc specimen 
ISRM:  International Society for Rock Mechanics 
SR:  Chevron notched short rod specimen 
TCD:  Theory of critical distances. 
Data Availability
The data supporting this research article are available from the corresponding author via email.
Conflicts of Interest
The authors declare that they have no conflicts of interest.
Acknowledgments
This study was financially supported by the National Key Research and Development Program of China (Grant no. 2016YFC0600701), the National Natural Science Foundation of China (Grant nos. 51674101, 51674170, and 51822403), and the Key Laboratory Open Project Fund in Henan Province (Grant no. S201605). The authors express their gratitude for the supportive help and instructions of the managers.
References
 W. S. Kreher, “Influence of domain switching zones on the fracture toughness of ferroelectrics,” Journal of the Mechanics and Physics of Solids, vol. 50, no. 5, pp. 1029–1050, 2002. View at: Publisher Site  Google Scholar
 M. Enderlein, A. Ricoeur, and M. Kuna, “Finite element techniques for dynamic crack analysis in piezoelectrics,” International Journal of Fracture, vol. 134, no. 34, pp. 191–208, 2005. View at: Publisher Site  Google Scholar
 M. R. M. Aliha, M. R. Ayatollahi, D. J. Smith, and M. J. Pavier, “Geometry and size effects on fracture trajectory in a limestone rock under mixed mode loading,” Engineering Fracture Mechanics, vol. 77, no. 11, pp. 2200–2212, 2010. View at: Publisher Site  Google Scholar
 M. R. M. Aliha and M. R. Ayatollahi, “Rock fracture toughness study using cracked chevron notched Brazilian disc specimen under pure modes I and II loading—a statistical approach,” Theoretical and Applied Fracture Mechanics, vol. 69, pp. 17–25, 2014. View at: Publisher Site  Google Scholar
 M. R. M. Aliha, E. Mahdavi, and M. R. Ayatollahi, “The influence of specimen type on tensile fracture toughness of rock materials,” Pure and Applied Geophysics, vol. 174, no. 3, pp. 1237–1253, 2017. View at: Publisher Site  Google Scholar
 F. Ouchterlony, “Suggested methods for determining the fracture toughness of rock,” International Journal of Rock Mechanics and Mining Sciences & Geomechanics Abstracts, vol. 25, no. 2, pp. 71–96, 1988. View at: Publisher Site  Google Scholar
 R. J. Fowell, “Suggested method for determining mode I fracture toughness using Cracked Chevron Notched Brazilian Disc (CCNBD) specimens,” International Journal of Rock Mechanics and Mining Sciences & Geomechanics Abstracts, vol. 32, no. 1, pp. 57–64, 1995. View at: Publisher Site  Google Scholar
 M. D. Kuruppu, Y. Obara, M. R. Ayatollahi, K. P. Chong, and T. Funatsu, “ISRMSuggested method for determining the mode I static fracture toughness using semicircular bend specimen,” Rock Mechanics and Rock Engineering, vol. 47, no. 1, pp. 267–274, 2014. View at: Publisher Site  Google Scholar
 M. R. M. Aliha, M. R. Ayatollahi, and R. Pakzad, “Brittle fracture analysis using a ringshape specimen containing two angled cracks,” International Journal of Fracture, vol. 153, no. 1, pp. 63–68, 2008. View at: Publisher Site  Google Scholar
 M. R. M. Aliha, R. Pakzad, and M. R. Ayatollahi, “Numerical analyses of a cracked straightthrough flattened Brazilian disk specimen under mixedmode loading,” Journal of Engineering Mechanics, vol. 140, no. 1, pp. 219–224, 2014. View at: Publisher Site  Google Scholar
 M. R. M. Aliha, G. R. Hosseinpour, and M. R. Ayatollahi, “Application of cracked triangular specimen subjected to threepoint bending for investigating fracture behavior of rock materials,” Rock Mechanics and Rock Engineering, vol. 46, no. 5, pp. 1023–1034, 2013. View at: Publisher Site  Google Scholar
 M. R. M. Aliha, A. Bahmani, and S. Akhondi, “Determination of mode III fracture toughness for different materials using a new designed test configuration,” Materials & Design, vol. 86, no. 5, pp. 863–871, 2015. View at: Publisher Site  Google Scholar
 M. R. M. Aliha, A. Bahmani, and S. Akhondi, “Numerical analysis of a new mixed mode I/III fracture test specimen,” Engineering Fracture Mechanics, vol. 134, pp. 95–110, 2015. View at: Publisher Site  Google Scholar
 M. R. M. Aliha and A. Bahmani, “Rock fracture toughness study under mixed mode I/III loading,” Rock Mechanics and Rock Engineering, vol. 50, no. 7, pp. 1739–1751, 2017. View at: Publisher Site  Google Scholar
 K. P. Chong and M. D. Kuruppu, “New specimen for fracture toughness determination for rock and other materials,” International Journal of Fracture, vol. 26, no. 2, pp. 59–62, 1984. View at: Publisher Site  Google Scholar
 M. R. Ayatollahi and M. R. M. Aliha, “Fracture toughness study for a brittle rock subjected to mixed mode I/II loading,” International Journal of Rock Mechanics and Mining Sciences, vol. 44, no. 4, pp. 617–624, 2007. View at: Publisher Site  Google Scholar
 M. R. M. Aliha, M. R. Ayatollahi, and J. Akbardoost, “Typical upper boundlower bound mixed mode fracture resistance envelopes for rock material,” Rock mechanics and Rock engineering, vol. 45, no. 1, pp. 65–74, 2012. View at: Publisher Site  Google Scholar
 M. R. M. Aliha, M. HeidariRarani, M. M. Shokrieh, M. R. Ayatollahi, and M. R. Ayatollahi, “Experimental determination of tensile strength and K(IC) of polymer concretes using semicircular bend (SCB) specimens,” Structural Engineering and Mechanics, vol. 43, no. 6, pp. 823–833, 2012. View at: Publisher Site  Google Scholar
 M. R. M. Aliha, M. Sistaninia, D. J. Smith, M. J. Pavier, and M. R. Ayatollahi, “Geometry effects and statistical analysis of mode I fracture in guiting limestone,” International Journal of Rock Mechanics and Mining Sciences, vol. 51, pp. 128–135, 2012. View at: Publisher Site  Google Scholar
 M. R. M. Aliha and M. R. Ayatollahi, “Twoparameter fracture analysis of SCB rock specimen under mixed mode loading,” Engineering Fracture Mechanics, vol. 103, pp. 115–123, 2013. View at: Publisher Site  Google Scholar
 J. Akbardoost, M. R. Ayatollahi, M. R. M. Aliha, M. J. Pavier, and D. J. Smith, “Sizedependent fracture behavior of Guiting limestone under mixed mode loading,” International Journal of Rock Mechanics and Mining Sciences, vol. 71, pp. 369–380, 2014. View at: Publisher Site  Google Scholar
 M. Fakhri, E. Amoosoltani, and M. R. M. Aliha, “Crack behavior analysis of roller compacted concrete mixtures containing reclaimed asphalt pavement and crumb rubber,” Engineering Fracture Mechanics, vol. 180, pp. 43–59, 2017. View at: Publisher Site  Google Scholar
 M. M. Mirsayar, X. Shi, and D. G. Zollinger, “Evaluation of interfacial bond strength between Portland cement concrete and asphalt concrete layers using bimaterial SCB test specimen,” Engineering Solid Mechanics, vol. 5, no. 4, pp. 293–306, 2017. View at: Publisher Site  Google Scholar
 A. S. Fayed, “Numerical evaluation of mode I/II SIF of quasibrittle materials using cracked semicircular bend specimen,” Engineering Solid Mechanics, vol. 6, no. 2, pp. 175–186, 2018. View at: Publisher Site  Google Scholar
 Z. M. Zhu and H. S. Wang, “Effect of notch root radius on rock fracture toughness K_{Ic},” Journal of Experimental Mechanics, vol. 8, no. 1, pp. 92–95, 1993. View at: Google Scholar
 F. Dai and Q. Z. Wang, “Effects of finite notch width on stress intensity factor for CCNBD specimen,” Rock and Soil Mechanics, vol. 3, no. 5, pp. 427–431, 2004. View at: Google Scholar
 S. M. Dong, M. Y. Fang, and Y. M. Xia, “Effect of crack shape and size on stress intensity factors of central cracked circular disk,” Journal of University of Science and Technology of China, vol. 34, no. 4, pp. 449–455, 2004. View at: Google Scholar
 G. Q. Zhang, M. Chen, and X. Y. Yang, “Influence of fracture width on rock toughness measurement,” Journal of the University of Petroleum China (Edition of Natural Science), vol. 26, no. 6, pp. 42–45, 2002. View at: Google Scholar
 S. Zhang and Y. L. Liang, “Influence of prefabricated crack width on determining rock fracture toughness for holecracked flattened Brazilian disc,” Journal of Experimental Mechanics, vol. 28, no. 4, pp. 517–523, 2013. View at: Google Scholar
 J. Zhou, Y. Wang, S. M. Dong, and Y. M. Xia, “Effect of specimen geometry on mode I fracture toughness of PMMA,” Polymer Material Science and Engineering, vol. 21, no. 6, pp. 11–14, 2005. View at: Google Scholar
 Z. L. Cui, N. P. Gong, and L. W. Jing, “Experiment and finite element analysis of rock dynamic fracture toughness test on non ideal crack disc specimens,” Rock and Soil Mechanics, vol. 36, no. 3, pp. 694–702, 2015. View at: Google Scholar
 Q. Z. Wang, D. B. Yan, and L. Luo, “Stress concentration factors for beams with single side Ushaped notched subjected to pure bending and a new method for mechanical test of quasibrittle materials,” Science Technology and Engineering, vol. 8, no. 12, pp. 3096–3102, 2008. View at: Google Scholar
 L. Luo and Q. Z. Wang, “Concurrent measurement of tensile strength and fracture toughness of quasibrittle materials using Unotched beams: theoretical analysis,” Engineering Mechanics, vol. 26, no. 9, pp. 244–250, 2009. View at: Google Scholar
 L. Luo, J. R. Yang, and Q. Z. Wang, “Concurrent measurement of tensile strength and fracture toughness of quasibrittle materials using Unotched beams: experimental analysis,” Journal of Basic Science and Engineering, vol. 25, no. 1, pp. 89–101, 2017. View at: Google Scholar
 R. Kolhe, C.Y. Hui, and A. T. Zehnder, “Effects of finite notch width on the fracture of chevron—notched specimens”,” International Journal of Fracture, vol. 94, no. 2, pp. 189–198, 1998. View at: Publisher Site  Google Scholar
 E. Barati and Y. Alizadeh, “A notch root radius to attain minimum fracture loads in plates weakened by Unotches under mode I loading,” Scientia Iranica, vol. 19, no. 3, pp. 491–502, 2012. View at: Publisher Site  Google Scholar
 S. Cicero, V. Madrazo, and I. A. Carrascal, “Analysis of notch effect in PMMA using the theory of critical distances,” Engineering Fracture Mechanics, vol. 86, pp. 56–72, 2012. View at: Publisher Site  Google Scholar
 M. Dehghany, H. Saeidi Googarchin, and M. R. M. Aliha, “The role of first nonsingular stress terms in mixed mode brittle fracture of Vnotched components: an experimental study,” Fatigue & Fracture of Engineering Materials & Structures, vol. 40, no. 4, pp. 623–641, 2017. View at: Publisher Site  Google Scholar
 M. K. Hussain and K. S. R. K. Murthy, “A point substitution displacement technique for estimation of elastic notch stress intensities of sharp Vnotched bodies,” Theoretical and Applied Fracture Mechanics, vol. 97, pp. 87–97, 2018. View at: Publisher Site  Google Scholar
 M. R. Ayatollahi, M. Heydarimeybodi, M. Dehghany, and F. Berto, “A new criterion for rupture assessment of rubberlike materials under modeI crack loading: the effective stretch criterion,” Advanced Engineering Materials, vol. 18, no. 8, pp. 1364–1370, 2016. View at: Publisher Site  Google Scholar
 F. Berto and M. R. Ayatollahi, “Fracture assessment of Brazilian disc specimens weakened by blunt Vnotches under mixed mode loading by means of local energy,” Materials & Design, vol. 32, no. 5, pp. 2858–2869, 2011. View at: Publisher Site  Google Scholar
 M. R. M. Aliha, F. Berto, A. Mousavi, and S. M. J. Razavi, “On the applicability of ASED criterion for predicting mixed mode I + II fracture toughness results of a rock material,” Theoretical and Applied Fracture Mechanics, vol. 92, pp. 198–204, 2017. View at: Publisher Site  Google Scholar
 M. R. M. Aliha and H. Saghafi, “The effects of thickness and Poisson’s ratio on 3D mixedmode fracture,” Engineering Fracture Mechanics, vol. 98, pp. 15–28, 2013. View at: Publisher Site  Google Scholar
 F. Berto, L. P. Pook, and A. Campagnolo, “Corner point singularities under inplane and outofplane loading: a review of recent results,” Engineering Solid Mechanics, vol. 5, no. 3, pp. 167–176, 2017. View at: Publisher Site  Google Scholar
Copyright
Copyright © 2019 Sheng Zhang et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.