1. Introduction
Over the last decade, because of its appealing characteristics, such as high strength-to-weight ratio and excellent corrosion resistance, composite materials such as Glass Fiber Reinforced Polymer (GFRP) have noticed considerable growth in different applications. It is necessary to understand the fracture behavior of composite materials, such as matrix cracking, delamination, translaminar fracture, and fiber breakage, especially in aircraft structures. A fracture in a structure will tend to propagate and eventually lead to the structure collapse. This problem is resolved by fracture mechanics testing, which delivers information in critical stress intensity factors near the crack tip at fracture. Therefore, it needs to use accurate standard tests to measure the fracture properties of composite materials [
1,
2,
3,
4,
5]. Compact tension specimen was used to determine the translaminar fracture toughness of carbon/epoxy composite laminates [
2,
3]. It was found that the major characteristics describing the fracture surfaces in the laminates comprised a mixture of the failure mechanisms, ply splitting, fibers bridging, and fiber pull-out. Haldar et al. [
6] showed that the fiber bundle pull-out was the process that dissipates the most energy. Furthermore, Souza et al. [
7] modified the ASTM E399 test method to account for the orthotropy of the composite materials and specimen geometry effects using a correction function based on a numerical evaluation of the strain energy release rate. They increased the initial notch depth from ≈16 mm to ≈23 mm to avoid compression and shear failures instead of mode I translaminar fracture under cyclic loading. Vantadori et al. [
8] modified Jenq-Shah’s two-parameter model (MTPM) [
9] to predict the fracture toughness of fiber-reinforced polymers (FRP) using three-point bend specimen according to RILEM [
10]. It was found [
8] that the predicted values of the fracture toughness and modulus of elasticity by MTPM are almost constant; consequently, such parameters are proved to be size-independent. The intrinsic and extrinsic mechanisms of mode I crack growth in FRP have been reviewed by Siddique et al. [
11]. They reviewed the main parameters controlling the fracture toughness of such materials. Furthermore, Jia et al. [
12] introduced a valuable strategy based on biomimicking to improve the fracture toughness of brittle materials through an intrinsic to extrinsic (ITE) transition. In the ITE transition, toughness started as an intrinsic parameter at the basic material level, but by designing a protein-like effective stress-strain behavior, the toughness at the system level became an extrinsic parameter that increases with the system size with no limit.
The study of translaminar fracture was carried out under the American Society for Testing Materials (ASTM) standard through a single-edge tension specimen. El-Hajjar and Haj-Ali [
13] evaluated the applicability of the ASTM E1922 [
14] standard for usage with pultruded thick-section composites, which expanded the standard’s scope to encompass these materials. Laffan [
15] proposed a study for translaminar tensile failure utilizing a single-edge tension specimen with different thicknesses (by changing the numbers of the plies). He observed that the fracture toughness for specimens increased in thicker 0° plies, and he attributed that to the increase of fiber pull-out. Furthermore, Laffan et al. [
16] expressed that using a single-edge tension specimen based on ASTM E1922 was appropriate for quasi-isotropic laminates. However, it rendered less accurate results for extremely orthotropic lay-ups. In the unidirectional fiber polymer composite under pure mode I, according to ASTM E1922, Saadati et al. [
17] found that crack propagation started from the notch tip and followed the fiber direction. He et al. [
18] demonstrated that brittle matrix could increase fracture work which reduced notch sensitivity. Zappalorto [
19] also showed that the shift from notch insensitivity to notch sensitivity was unaffected by the notch root radius, the notch depth, and critical notch size. Moreover, Jamali et al. [
20] found crack propagation is perpendicular to the intended development direction under tension.
On the other hand, ASTM D3039 [
21] was used to measure the tensile behavior of fiber-reinforced polymers (FRP). Elbadry et al. [
22] and El-Wazery et al. [
23] studied the tensile behavior of FRP with hand lay-up with different fiber volume fractions, and concluded that increasing the fiber weight contents increased the tensile properties.
The Numerical method is used to simulate processes that cannot be examined or seen in the lab, besides the cost and time of experiments that it saves. Hashin damage model [
24,
25] provides a good prediction for four modes of failure for matrix and fibers. As Pham et al. [
26] demonstrated for the FRP, a discrete crack under a 3D Hashin failure criterion was first used to predict the damage initiation with continuum shell elements. The matrix crack was parallel to the fiber path, and its orientation was determined by the transverse plane’s maximum principal stress, which was validated experimentally. Furthermore, Koloor et al. [
27] used Hashin’s modulation by finite element and energy dissipation to regularly get the multidirectional composite material yield value. Duarte et al. [
28] simulated the FRP using the ABAQUS program under unidirectional tension based on Hashin’s criteria. They showed that plates with a higher plies number with fibers orientation in the direction of applied load plies had the highest stiffness and strength. They [
28] compared the Hashin damage criterion and the eXtended Finite Element Method (XFEM) to predict FRP failure stages. They concluded that the Hashin damage criterion and XFEM predicted the same strength and stiffness of FRP for load levels up to the failure of plies due to matrix cracking. Besides, Contour Integral Method (CIM) is a valid model to estimate J-integrals according to stress intensity factors related to contours above the crack region. El-Sagheer et al. [
29] and Abd-Elhady et al. [
30,
31] used the CIM to determine the stress intensity factor and J-integrals in many applications.
On the other hand, Carpenter et al. [
32] confirmed that the matrix under compressive loads follows the relationships established by linear elastic fracture mechanics (LEFM). They found that a log-log plot of the failure load and initial notch length of the experimental data exhibited a linear trend with a slope of −0.54, while the numerical predictions had a linear trend of slope −0.58. As is already known, this plot was expected from LEFM to be linear with a negative slope of 0.5. Liu et al. [
33] invoked 3-D finite element analysis (FEA) to study the effect of cohesive zone model (CZM) parameters on the post-buckling and delamination behaviors of FRP under compression. They found that the cohesive strengths mainly affected the unstable delamination stage for the laminates under compression and had little effect on local and global buckling loads. Rozylo [
34] concluded that the results obtained based on CZM (numerical) and acoustic emission signals (experimental) showed high agreement. Panettieri et al. [
35] used CZM to simulate delaminations growth in compression after impact. Zhou et al. [
36] used the shear damage initiation criterion, available in ABAQUS/explicit, to model the shear failure due to fracture within shear bands in metal-ceramic functionally graded bolted joint. It is worth noting that the bolt was made of porous ZrO
2/(ZrO
2 + Ni) FGMs. They used Tsai–Wu tensor theory as the failure criteria of the C/SiC plates. They concluded that ZrO
2 + 15 vol% Ni of two mm thickness is the optimal shear band to balance such bolted joint’s shearing strength and heat insulation performance.
The main objective of the present work is to experimentally study the ability of the ASTM E1922 standard test method to measure the real fracture toughness of unidirectional glass fiber reinforced epoxy (GFRE). It is worth noting that two types of cracked ASTM E1922 specimens have been adopted in the present work. The crack types are single edge through-thickness cracked (TTC) specimen (traditional specimen) and single edge matrix cracked (MC) specimen (suggested in the present work as an unprecedented specimen). The pre-crack exists in the matrix without cutting the fibers behind the crack tip in the MC specimen. In the first stage of the experimental work, MC specimens have been used to show the failure at the loading point. While in the second stage, three types of tensile specimens, namely, smooth, TTC, and MC specimens, have been manufactured with dimensions according to ASTM D3039 to measure the exact value of the fracture toughness of GFRE through the MC ASTM D3039 specimen. No such MC specimen in the fibrous composite is suggested before. Therefore, three-dimensional finite element analysis (3D FEA) simulates the composite laminates based on ASTM E1922 and ASTM D3039 standard tests. Moreover, Hashin criteria and CIM are used in the present simulation to predict the progressive damage and the geometry correction factor (Y) of the GFRE cracked specimen, respectively.