Research Article  Open Access
Zhitao Zhang, Haijun Wang, Shuyang Yu, "Numerical Simulation of Fracture Behavior of Brittle Solids under I/III Loading", Advances in Civil Engineering, vol. 2020, Article ID 8612784, 14 pages, 2020. https://doi.org/10.1155/2020/8612784
Numerical Simulation of Fracture Behavior of Brittle Solids under I/III Loading
Abstract
This study used numerical analysis to carry out a large number of numerical model calculations based on a new semicircular bend (SCB) model. Instead of existing numerical computation methods (Jintegral), the Mintegral method was used to calculate the mixed mode fracture parameters (K_{I}, K_{II}, and K_{III}), investigate the influence of the geometry and material parameters on the fracture behavior under mixed mode I/III loading, and predict the fracture path. The results revealed that the limitations in the scope of the mixity parameter M^{e} in the previous studies can be overcome to a certain extent. The range of M^{e} was established under different Poisson ratios and can be used as a reference for actual material testing. The simulation path is in good agreement with the experimentally obtained fracture path, and the proposed method can be used to simulate the fracture path under mixed mode I/III loading.
1. Introduction
Fracture mechanics focuses on the mechanical properties of materials and structures with cracks. The cracks can be inherent in the material or generated during manufacturing, and the existence and growth of these cracks decrease the bearing capacity of the structure or even cause its failure.
The semicircular bend (SCB) is a classic fracture mechanics test that is widely used to investigate solid materials owing to its inherent favorable properties, such as its simplicity, minimal machining requirement, testing convenience, and ease of reaching tensile failure as reported by Mahinda and Kuruppu [1]. The SCB specimen was proposed by Chong [2–4] and has been gradually improved. Additionally, it has been widely used for determining the fracture parameter K_{IC}, and its potential application in the mixed mode fracture test has been investigated [5–9]. The investigation of the SCB specimen mainly includes fracture toughness testing and determining the fracture parameters [10–28], discussing and revising the limitations of the extension criteria [11, 13, 16, 25, 29, 30], and numerically predicting the fracture path [14, 20, 23, 27, 31–33]. Owing to the limitation of the specimen shape, the abovementioned studies focused on the fracture parameters of mode I, mode II, and mixed mode I/II. The outofplane shearing mode (mode III) has not been as extensively investigated, although it exists in actual situations. Additionally, existing testing methods are highly dependent on the machine and test fixture and mainly focus on metal materials, ceramics, and other artificial materials [34–38].
Because engineering accidents, such as rock fracture and pavement rupture, occur frequently, the investigation of mode III has attracted a great amount of interest. Moreover, observations have shown that flat cracks tend to reorient themselves to oblique planes during propagation and can grow under mixedmode I/III conditions [39, 40]. Therefore, test methods and specimens under the abovementioned conditions have been proposed for rocks, concrete, asphalt, and other materials. Berto et al. applied the strain energy density (SED) criterion to pure modes and mixed mode I/III loading conditions [41–45]. Aliha et al. [46–49] proposed a new test configuration called the ENDB specimen for investigating the fracture behavior under the mixed mode I/III. Additionally, this configuration can achieve an arbitrary composition of mode I and mode III. Notably, however, more variables need to be controlled in the experiment. Linul et al. [50, 52] conducted experiments with foam materials and obtained fracture toughness data for the mixed modes I/II and I/III. Subsequently, Pirmohammad [53] proposed a new SCB specimen type, whose mixed mode I/III fracture behavior is easier to investigate using the proposed specimens with a tilted crack and further investigated the fracture behavior of asphalt concrete at low temperature.
In previous studies [49, 53], the change of the geometric parameters and material parameters had a significant effect on the mixed mode I/II loading of SCB specimens. However, only few studies have investigated whether the abovementioned conditions affect the fracture behavior of mode I/III. In this study, a relatively accurate numerical analysis method was used to investigate the feasibility of realizing mixed mode I/III loading and the effects of the specimen’s geometric parameters and Poisson’s ratio on these loading conditions. Finally, the fracture path was simulated and analyzed based on the fracture parameters.
2. 3D Numerical Analyses
In the previous numerical studies [41, 53, 54], the fracture parameters were directly calculated using the Jintegral method in the Abaqus software. However, according to Li, the Jintegral is more suitable for pure mode I fracture problems [55]. For the mixed mode fracture, the equivalent domain equation of the interaction integral (Mintegral) [56, 57] is the most accurate method for calculating the fracture parameters. The SCB macrostructural model was established in Abaqus; then, the mesh was remeshed using adaptive grid technology in franc3D, and the bottom inclined crack was inserted. Finally, different K_{I}, K_{II}, and K_{III} values were obtained by the Mintegral.
The schematic diagram of the semicircular geometry is shown in Figure 1; the radius of the SCB specimen (R), specimen thickness (t), and half of the loading span (S) in the SCB test were determined by modeling in Abaqus. Figure 2(a) shows the Abaqus grid model with C3D10 cells. Figure 2(b) shows the mesh model after defining the crack length (a) and the angle of crack (α). Figure 3 shows the element distribution and crack tip types. A singular element was inserted at the crack tip, and conservation integral evaluation was carried out around the two element rings at the crack front. The inner ring of a 15node singular wedgeshaped element and the outer ring of a 20node hexahedral element with symmetric meshes were used to reduce the local discrete errors. Notably, the fracture parameters (K_{I}, K_{II}, and K_{III}) calculated at the surface contact zone are not reliable, because the local plane strain condition is no longer maintained. The element at the crack tip may be severely distorted, as shown in Figure 2(c), and the grid symmetry may be lost. Moreover, Figure 2(c) shows that the transition region at the crack tip has pyramidal elements, while the global region is a tetrahedral element type.
(a)
(b)
(c)
In the created models, radius of SCB specimen (R), specimen thickness (t), and applied load (P) are invariants, respectively, 75 mm, 32 mm, and 40.3437 kN. Geometrical and material parameters such as the half of loading point span ratio (S/R), crack length ratio (a/R), crack angle (α), and Poisson’s ratio were varied in the following ranges: S/R = {0.4, 0.6, 0.8}, a/R = {0.2, 0.26, 0.3}, α = {0, 5, 10, 15, 20, 25, 30, 35, 40, 45, 48, 50, 52, 55, 60}, and = {0.1, 0.2, 0.3, 0.4}. When α is 0°, it is pure mode I fracture. With the increase of α, mixed mode I/III fracture can be obtained.
The fracture parameters (K_{I}, K_{II}, and K_{III}) can be directly extracted using the following expressions. For linear analysis, we can add two valid solutions, and the result is a valid solution, as follows:
Let us consider the corner mark (1) solutions as the Abaqus results and the corner mark (2) solutions as the solutions that we can select. These can be substituted into the expression for the Jintegral as follows:
According to Betti’s reciprocal theorem, the following relationship holds:
By collecting terms, the following relationship can be obtained:with
The cracktip energy release rates can be determined from Irwin’s crack closure integral for small scale yielding assuming plane strain conditions, as follows:
By substituting into the expression for the energy release rate, we obtain the following relationship:
By equating the two definitions for the MIntegral, we can obtain the following relationship:where
We used the Abaqus results for solution (1) and selected the three simple auxiliary solutions (2a), (2b), and (2c). Table 1 presents the values of the asymptotic solutions for the pure mode cracks.

From the analytical expressions for the crackfront fields, we can obtain the following:
By substituting the value of K^{(2)}, reported in Table 1, in equation (9), we can obtain three equations for the unknown K^{(1)}’s, as follows:where is Poisson’s ratio, is the integral curve away from the tip of the crack, is the Kronecker delta, q is a function that is equal to one at the crack tip and zero on the boundary of the integration domain and can be interpreted as a virtual crack extension, is the value of the q function along the crack front, and L is the length of the cylindrical domain along the crack front. is the solution provided by the Abaqus results.
K_{I}, K_{II}, and K_{III} can be expressed by the dimensionless parameters Y_{I}, Y_{II}, and Y_{III}, as follows:
To verify that the numerical simulation method can accurately realize mixed I/III mode fracturing, the K_{I}, K_{II}, and K_{III} of the crack tip along the specimen thickness were adopted after normalization processing as shown in Figure 4. Here, K_{n} (n = I, II, III) is the fracture parameter and K_{Im} is the maximum value of mode I. The crack length ratio (a/R) of the model is 0.26, one half of the loading point span ratio (S/R) is 0.8, and the crack angle (α) is 0°, 20°, 40°, and 60°, respectively.
(a)
(b)
(c)
(d)
Figure 4(a) shows that, when α was 0°, fracture occurred under the pure mode I condition; K_{I} along the direction of the specimen thickness remained approximately unchanged within a certain range (0.1 < z/t < 0.9) and declined in the free surface, which is consistent with the numerical simulation results reported by Aliha and Saghafi [54]. Moreover, z/t was considered as 0.5 for the K_{I} value of the specimen.
Figures 4(b)–4(d) shows that, as α increased, the effects of modes I and III dominated within a certain range (0.2 < z/t < 0.8) along the thickness direction of the specimen. The proportion of mode II was small and could be ignored when z/t was 0.5, which is consistent with previously reported numerical analysis results [46–49]. Considering that the surface contact zone calculation of K_{I}, K_{II}, and K_{III} is not reliable, and according to the traditional definition of the singular stress field close to the free surface, it is widely accepted that the crack tip stress field is different, as discussed by Bažant and Estenssoro [58]. Hence, K_{I}, K_{II}, and K_{III} close to the free surface should not be considered, and this numerical method is valid and feasible for testing the mixed mode I/III fracture. Notably, in mixed mode I/III loading, the crack tip along the thickness of the specimen completely corresponds to the point where z/t is 0.5. Therefore, subsequent calculations considered this point as the value of K_{I}, K_{II}, and K_{III} of the specimen.
The mixity parameter M^{e} can be used to describe the relative contributions of mode I and mode III, as expressed by equation (15). Under the condition of pure mode I, the value of M^{e} is equal to one, and decreases as the contribution of mode III increases:
To quantify the role of mode I and mode III in the numerical simulation, we propose to substitute equation (16) for equation (15), and the parameter M^{e} can vary from 0 to 1, corresponding to the pure mode III and mode I loading conditions, respectively:
3. Results and Discussion
3.1. Numerical Calculation Results and Analysis
The crack length ratio (a/R) remained unchanged at 0.267 and a = 20 mm. The influences of the crack angle (α), half of the loading point span ratio (S/R), and Poisson’s ratio for Y_{I}, Y_{III}, and M^{e} were investigated. Figures 5 and 6 show the variations of Y_{I}, Y_{III}, and M^{e}, with a crack inclination angle (α) for different values of S/R and in the SCB specimen.
(a)
(b)
(c)
(d)
(e)
(f)
(a)
(b)
(c)
The following conclusions were drawn from image analysis:(a)Y_{III} remained at zero, only when α was 0°. For the same α, the increase of S/R and promoted the Y_{I}. But for the Y_{III}, the increase of S/R and led to its decrease and increase, respectively.(b)For the pure mode I (α = 0°), the increase of promoted the contribution of mode I, and the growth rate tended to increase.(c)With the increase of α, after entering the mixed mode I/III loading, the values of Y_{I} gradually decreased, and Y_{III} first increased and then decreased. Additionally, α corresponding to Y_{III}_{max} slowly increased with the increase of S/R. The sensitivity of Y_{I} to the changes of decreased as α increased, and the sensitivity of Y_{III} was maximum at Y_{III}_{max}.(d)The increase of S/R only had a minor effect on the M^{e} values when was 0.1 of the SCB specimens. When the range of varied from 0.2 to 0.4, the increase of S/R obviously led to the decrease of M^{e}. Additionally, as α increased, the downtrend became more obvious.
Half of the loading point span ratio (S/R) remained unchanged at 0.8, S = 60 mm, and the influences of the crack angle (α), crack length ratio (a/R), and Poisson’s ratio () on Y_{I}, Y_{III}, and M^{e} were investigated. Figures 7 and 8 show the variation of Y_{I}, Y_{III}, and M^{e} with the crack inclination angle (α) for different values of a/R and in the SCB specimen. Figure 9 shows the M^{e}_{min} obtained with the change of a/R and .
(a)
(b)
(c)
(d)
(e)
(f)
(a)
(b)
(c)
The following conclusions were drawn from image analysis:(a)For the same α, the increase of a/R and promoted the Y_{I} values. But for the Y_{III}, the increase of a/R and led to its decrease and increase, respectively.(b)The value of α corresponding to Y_{III}_{max} stabilized at approximately 49°.(c)For the values of M^{e}, the increase of a/R led to the promotion of M^{e}; when the value of a/R was relatively small, the influence of on M^{e} was smaller under the same α.(d)The value of α was 60°, corresponding to the M^{e}_{min} of all SCB specimens. As and a/R became smaller, a smaller M^{e}_{min} value was obtained and the minimum value was 0.39. In this case, the contribution of model III was obviously higher compared with that of model I. Additionally, for smaller a/R, only had a minor influence on M^{e}_{min}; however, when a/R was larger, M^{e}_{min} was more sensitive to the variations of .
Notably, the rule discussed in the previous section is not described.
3.2. Fracture Path Simulation
Using the K_{I}, K_{II}, and K_{III} values of the crack tip and the maximum tensile stress criterion (MTS), three model groups (S/R = 0.8; a/R = 0.267; α = 0°, 30°, 45°, and 60°) were obtained for the crack failure and growth path, as shown in Figures 10 and 11. To facilitate the comparison, the failure mode after the image was flipped horizontally is in good agreement with the previous results [53].
(a)
(b)
(c)
(d)
Under mixed mode I/III loading conditions, the crack surface growth morphology was different from the inplane fracture. Figure 12(a) shows that, when an inplane fracture (mode I, mode II, and mixed mode I/II) occurred, the crack surface remained in the same plane despite a certain deflection. However, Figure 12(b) shows that, when the mixed mode I/III fracture occurred, the crack surface was twisted and deviated from the direction of the crack surface by forming an angle called the kink angle.
Considering the direction of the crack kink angle (), the normal stress and shear stress components on the new plane can be obtained by carrying out coordinate transformation as follows:where x, y, and z are the original axes and x′, y′, and z′ are the axes of the crack plane after deflection. The transformed stress intensity factor can be expressed as follows:
The relationship between the release rate of the mechanical energy and the angle is expressed as follows:
By superposing the contributions of mode III, the corresponding transformed stress intensity factor can be expressed as follows.
By substituting equation (18) into equation (19), the expression of the release rate of the mechanical energy changing with the angle can be obtained for the mixed mode I/III. Figure 13 shows the normalization of the mechanical energy release rate along with the change of the crack kink angle ().
(a)
(b)
As can be seen in the graph, in the case of pure mode I fracturing (K_{III}/K_{I} = 0), became maximum at the kink angle () of 0°. Therefore, the pure mode I fracture extended in the direction of the original when the kink angle was 0°. When the superposition of the mode III contribution was added (K_{III}/K_{I} = 0.25), as shown in Figure 13, the crack kink angle deviated from the original crack plane by an angle that can be considered as a correction by the contribution of mode III. Finally, the crack moved up in the direction perpendicular to the maximum principal tension and changed from mixed mode I/III to pure mode I. The numerical simulation results presented in Figure 10 reveal that the extension surface gradually returned to a plane from a curved surface close to the loading point and finally fractured close to the loading point. This phenomenon is consistent with the fact that the fracture surface is divided into three regions [53]. Table 2 presents the kink angle between 0.2 and 0.8 at the crack front with different α, as extracted from Franc3D. The SCB specimen parameters were S/R = 0.8, a/R = 0.267, and = 0.2. Additionally, it was observed that the kink angle of the crack is related to α, and a larger α resulted in a greater kink angle.

4. Conclusions
The following conclusions were drawn from this study:(a)Numerical modeling combined with Abaqus and Franc3D analysis was used to investigate the influence of half of the loading point span ratio (S/R), crack length ratio (a/R), crack angle (α), and Poisson’s ratio () on the behavior of mixed mode I/III loading.(b)The increase of S/R, a/R, and promoted the dimensionless parameters Y_{I}. But for the Y_{III}, only the increase of a/R and S/R led to its decrease. As α increased, the sensitivity of Y_{I} decreased with the changes of , while the sensitivity of Y_{III} became maximum at Y_{III}_{max}.(c)The obvious approach toward reducing the mixity parameter M^{e} is to decrease α. Additionally, it is possible to reduce M^{e} by changing S/R or a/R. Moreover, the SCB specimen materials with small exhibited the best reduction effect.(d)The M^{e}_{min} values of the SCB specimens were obtained under different a/R and , which can be useful as a references for SCB specimens of different materials subjected to mixed mode I/III loading tests.(e)The numerical parameters obtained by the Mintegral were used to simulate the crack growth and failure path, which are in good agreement with previous experimental results. The relationship between α and the kink angle () was quantitatively analyzed, and their positive correlation was confirmed. The proposed method can be used to simulate the fracture path under mixed mode I/III loading.
Data Availability
All data used to support the findings of this study are included within the article.
Conflicts of Interest
There is no conflicts of interest regarding the publication of this paper.
Acknowledgments
The authors gratefully acknowledge Senior Engineer Haijun Wang and Dr. Shuyang Yu for their supports in this research program, and they also extend their thanks to the National Key R&D Program of China (2017YFC0404902), the National Natural Science Foundation of China (51739008), the Natural Science Foundation of Jiangsu Province of China (BK20171130), and the basic scientific research operating expenses of the Public Welfare Scientific Research Institutes at the Central Level of China (Y419005).
References
 D. Mahinda and K. P. C. Kuruppu, “Fracture toughness testing of brittle materials using semicircular bend (SCB) specimen,” Engineering Fracture Mechanics, vol. 91, no. 1, pp. 133–150, 2012. View at: Publisher Site  Google Scholar
 K. P. Chong and M. D. Kuruppu, “New specimen for fracture toughness determination of rock and other materials,” International Journal of Fracture, vol. 26, no. 2, pp. R59–R62, 1984. View at: Publisher Site  Google Scholar
 K. P. Chong, M. D. Kuruppu, and J. S. Kuszmaul, “Fracture toughness determination of layered materials,” Engineering Fracture Mechanics, vol. 28, no. 1, pp. 43–54, 1987. View at: Publisher Site  Google Scholar
 K. P. Chong, K. D. Basham, and W. D. Qing, “Fracture parameters derived from tensionsoftening measurements using semicircular specimens,” International Journal of Rock Mechanics & Mining Sciences & Geomechanics Abstracts, vol. 28, no. 6, p. A398, 1991. View at: Google Scholar
 K. P. Chong and M. D. Kuruppu, “Mixed mode fracture analysis using new semicircular specimens,” Computers & Structures, vol. 30, no. 4, pp. 905–908, 1988. View at: Publisher Site  Google Scholar
 K. P. Chong and M. D. Kuruppu, “New specimens for mixed mode fracture investigations of geomaterials,” Engineering Fracture Mechanics, vol. 30, no. 5, pp. 701–712, 1988. View at: Publisher Site  Google Scholar
 B. N. Whittaker, R. N. Singh, and Sun Gexin, “Rock fracture mechanics principles, design and applications,” Rock Fracture Mechanics:principles, Design, and Applications, vol. 71, p. 275, 1992. View at: Google Scholar
 I. L. Lim, I. W. Johnston, and S. K. Choi, “Assessment of mixedmode fracture toughness testing methods for rock,” International Journal of Rock Mechanics and Mining Sciences & Geomechanics Abstracts, vol. 31, no. 3, pp. 265–272, 1994. View at: Publisher Site  Google Scholar
 I. L. Lim, I. W. Johnston, S. K. Choi, and J. N. Boland, “Fracture testing of a soft rock with semicircular specimens under threepoint bending. Part 1mode I,” International Journal of Rock Mechanics and Mining Sciences & Geomechanics Abstracts, vol. 31, no. 3, pp. 185–197, 1994. View at: Publisher Site  Google Scholar
 I. L. Boland, I. W. Johnston, and S. K. Choi, “Stress intensity factors for semicircular specimens under threepoint bending,” Engineering Fracture Mechanics, vol. 44, no. 3, pp. 363–382, 1993. View at: Publisher Site  Google Scholar
 M. R. Ayatollahi, M. R. M. Aliha, and M. M. Hassani, “Mixed mode brittle fracture in PMMA—an experimental study using SCB specimens,” Materials Science & Engineering A (Structural Materials:, Properties, Microstructure and Processing), vol. 417, no. 12, pp. 348–356, 2006. View at: Publisher Site  Google Scholar
 M. D. Kuruppu, Y. Obara, 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. Funatsu and M. R. Ayatollahi, “Mixed mode I/II brittle fracture evaluation of marble using SCB specimen,” Procedia Engineering, vol. 10, no. 1, pp. 311–318, 2011. View at: Publisher Site  Google Scholar
 W. Wang, J. E. Olson, M. Prodanović, and R. A. Schultz, “Interaction between cemented natural fractures and hydraulic fractures assessed by experiments and numerical simulations,” Journal of Petroleum Science and Engineering, vol. 167, no. 3, pp. 506–516, 2018. View at: Publisher Site  Google Scholar
 D. Singh, P. K. Ashish, and S. F. Chitragar, “Laboratory performance of recycled asphalt mixes containing wax and chemical based Warm mix additives using semi circular bending and tensile strength ratio tests,” Construction and Building Materials, vol. 158, no. 10, pp. 1003–1014, 2018. 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 & Mining Sciences, vol. 71, pp. 369–380, 2018. View at: Publisher Site  Google Scholar
 M. R. M. Aliha, H. Behbahani, and M. H. Rezaifar, “Study of characteristic specification on mixed mode fracture toughness of asphalt mixtures,” Construction and Building Materials, vol. 54, no. 12, pp. 623–635, 2014. View at: Publisher Site  Google Scholar
 B. Huang, G. X. Shu, and G. Zuo, “Using notched semi circular bending fatigue test to characterize fracture resistance of asphalt mixtures,” Engineering Fracture Mechanics, vol. 109, no. 12, pp. 78–88, 2013. View at: Publisher Site  Google Scholar
 C. Jiang, G.F. Zhao, J. Zhu, Y.X. Zhao, and L. Shen, “Investigation of dynamic crack coalescence using a gypsumlike 3D printing material,” Rock Mechanics and Rock Engineering, vol. 49, no. 10, pp. 3983–3998, 2016. View at: Publisher Site  Google Scholar
 B. Birgisson, A. Montepara, E. Romeo et al., “Determination and prediction of crack patterns in hot mix asphalt (HMA) mixtures,” Engineering Fracture Mechanics, vol. 75, no. 34, pp. 664–673, 2008. View at: Publisher Site  Google Scholar
 G. Feng, Y. Kang, T. Meng, Y.Q. Hu, and X.H. Li, “The influence of temperature on mode I fracture toughness and fracture characteristics of sandstone,” Rock Mechanics and Rock Engineering, vol. 50, no. 8, pp. 2007–2019, 2017. View at: Publisher Site  Google Scholar
 M. Ameri, M. S. Nowbakht, and M. R. M. Aliha, “Investigation of fatigue and fracture properties of asphalt mixtures modified with carbon nanotubes,” Fatigue & Fracture of Engineering Materials & Structures, vol. 39, no. 7, pp. 896–906, 2016. View at: Publisher Site  Google Scholar
 M. R. M. Aliha, A. Bahmani, and S. Akhondi, “Mixed mode fracture toughness testing of PMMA with different threepoint bend type specimens,” European Journal of Mechanics—A/Solids, vol. 58, no. 1, pp. 148–162, 2016. 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
 H. Saghafi, M. R. Ayatollahi, and M. Sistaninia, “A modified MTS criterion (MMTS) for mixedmode fracture toughness assessment of brittle materials,” Materials Science & Engineering: A (Structural Materials: Properties, Microstructure and Processing), vol. 527, no. 2122, pp. 5624–5630, 2010. View at: Publisher Site  Google Scholar
 X.J. Li and M. O. Marasteanu, “Using semi circular bending test to evaluate low temperature fracture resistance for asphalt concrete,” Experimental Mechanics, vol. 50, no. 7, pp. 867–876, 2010. 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
 I. Pavier and H. A. Khalid, “A comparison between beam and semicircular bending fracture tests for asphalt,” Road Materials and Pavement Design, vol. 7, no. 1, pp. 163–180, 2006. View at: Publisher Site  Google Scholar
 M. R. Ayatollahi and M. R. M. Aliha, “On determination of mode II fracture toughness using semicircular bend specimen,” International Journal of Solids and Structures, vol. 43, no. 17, pp. 5217–5227, 2006. View at: Publisher Site  Google Scholar
 C. Li, L. Xie, J. Li Ren et al., “Evaluating the applicability of fracture criteria to predict the crack evolution path of dolomite based on SCB experiments and FEM,” Advances in Civil Engineering, vol. 2013, Article ID 959806, 13 pages, 2013. View at: Publisher Site  Google Scholar
 M. Eftekhari, H. A. Baghbanan, and H. Hashemolhosseini, “Fracture propagation in a cracked semicircular bend specimen under mixed mode loading using extended finite element method,” Arabian Journal of Geosciences, vol. 8, no. 11, pp. 9635–9646, 2015. 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
 T. Xiaoge, R. Zhang, Z. Yang, Y. Chu, S. Zhen, and Y. Xv, “Simulation of bending fracture process of asphalt mixture semicircular specimen with extended finite element method,” Advances in Civil Engineering, vol. 2018, Article ID 4081264, 8 pages, 2018. View at: Publisher Site  Google Scholar
 W. S. Johnson, J. E. Masters, T. K. O’brien, and S. K. Lee, “An edge crack torsion method for mode III delamination fracture testing,” Journal of Composites Technology and Research, vol. 15, no. 3, pp. 193–201, 1993. View at: Publisher Site  Google Scholar
 J. G. Lee, “Characterization of the edge crack torsion (ECT) test for mode III fracture toughness measurement of laminated composites,” Journal of Endocrinology, vol. 82, no. 3, pp. 403–408, 2004. View at: Google Scholar
 D. Pennas, W. J. Cantwell, and P. Compston, “The influence of strain rate on the mode III interlaminar fracture of composite materials,” Journal of Composite Materials, vol. 41, no. 21, pp. 2595–2614, 2007. View at: Publisher Site  Google Scholar
 M. Farshad and P. Flueler, “Investigation of mode III fracture toughness using an anticlastic plate bending method,” Engineering Fracture Mechanics, vol. 60, no. 56, pp. 597–603, 1998. View at: Publisher Site  Google Scholar
 M. F. S. F. de Moura, M. V. C. Fernande, A. B. Morais, and R. D. S. G. Campilho, “Numerical analysis of the edge crack torsion test for mode III interlaminar fracture of composite laminates,” Engineering Fracture Mechanics, vol. 76, no. 4, pp. 469–478, 2009. View at: Publisher Site  Google Scholar
 X. Gao and C. Fong Shih, “A parametric study of mixedmode I/III ductile fracture in tough materials under small scale yielding,” Engineering Fracture Mechanics, vol. 60, no. 4, pp. 407–420, 1998. View at: Publisher Site  Google Scholar
 B. Lin, M. E. Mear, and K. RaviChandar, “Criterion for initiation of cracks under mixedmode I + III loading,” International Journal of Fracture, vol. 165, no. 2, pp. 175–188, 2010. View at: Publisher Site  Google Scholar
 F. Berto and P. Lazzarin, “Recent developments in brittle and quasibrittle failure assessment of engineering materials by means of local approaches,” Materials Science and Engineering: R: Reports, vol. 75, pp. 1–48, 2014. View at: Publisher Site  Google Scholar
 F. Ayatollahi, M. Elices, P. Lazzarin, and M. Zappalorto, “Fracture behaviour of notched round bars made of PMMA subjected to torsion at room temperature,” Engineering Fracture Mechanics, vol. 90, no. 5, pp. 143–160, 2012. View at: Publisher Site  Google Scholar
 F. Zappalorto, P. Lazzarin, and M. R. Ayatollahi, “Brittle fracture of sharp and blunt Vnotches in isostatic graphite under torsion loading,” Carbon, vol. 50, no. 5, pp. 1942–1952, 2012. View at: Publisher Site  Google Scholar
 F. Berto, D. A. Cendon, P. Lazzarin, and M. Elices, “Fracture behaviour of notched round bars made of PMMA subjected to torsion at −60°C,” Engineering Fracture Mechanics, vol. 102, no. 2, pp. 271–287, 2013. View at: Publisher Site  Google Scholar
 F. Elices, M. Ayatollahi, and A. Campagnolo, “Fracture tests under mixed mode I + III loading: an assessment based on the local energy,” International Journal of Damage Mechanics, vol. 26, no. 6, pp. 881–894, 2017. 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, no. 12, pp. 95–110, 2015. 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. 8, pp. 863–871, 2015. View at: Publisher Site  Google Scholar
 M. R. M. Aliha, A. Bahmani, and S. Akhondi, “A novel test specimen for investigating the mixed mode I + III fracture toughness of hot mix asphalt composites—experimental and theoretical study,” International Journal of Solids and Structures, vol. 90, no. 3, pp. 167–177, 2016. 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
 E. Linul and L. Marsavina, “Experimental determination of mixedmode fracture toughness for rigid polyurethane foams,” Fracture at all Scales, Springer, Berlin, Germany, 2016. View at: Publisher Site  Google Scholar
 M. R. M. Aliha, E. Linul, A. Bahmani, and L. Marsavina, “Experimental and theoretical fracture toughness investigation of PUR foams under mixed mode I + III loading,” Polymer Testing, vol. 67, no. 2, pp. 75–83, 2018. View at: Publisher Site  Google Scholar
 S. Marsavina and A. Kiani, “Study on fracture behavior of HMA mixtures under mixed mode I/III loading,” Engineering Fracture Mechanics, vol. 153, no. 4, pp. 80–90, 2016. View at: Publisher Site  Google Scholar
 M. R. Ayatollahi and M. R. M. Aliha, “Wide range data for crack tip parameters in two disctype specimens under mixed mode loading,” Computational Materials Science, vol. 38, no. 4, pp. 660–670, 2007. 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, no. 1, pp. 15–28, 2013. View at: Publisher Site  Google Scholar
 F. Z. Li, C. F. Shih, and A. Needleman, “A comparison of methods for calculating energy release rates,” Engineering Fracture Mechanics, vol. 21, no. 2, pp. 405–421, 1985. View at: Publisher Site  Google Scholar
 L. BanksSills, P. A. Wawrzynek, B. Carter, A. R. Ingraffea, and I. Hershkovitz, “Methods for calculating stress intensity factors in anisotropic materials: part IIarbitrary geometry,” Engineering Fracture Mechanics, vol. 74, no. 8, pp. 1293–1307, 2007. View at: Publisher Site  Google Scholar
 J. H. Ingraffea and A. J. Chien, “Evaluation of Mintegral for anisotropic elastic media with multiple defects,” International Journal of Fracture, vol. 114, no. 3, pp. 267–289, 2002. View at: Publisher Site  Google Scholar
 Z. P. Bažant and L. F. Estenssoro, “Surface singularity and crack propagation,” International Journal of Solids & Structures, vol. 15, no. 5, pp. 405–426, 1979. View at: Google Scholar
Copyright
Copyright © 2020 Zhitao 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.