Next Article in Journal
Characterization of the Static, Creep, and Fatigue Tensile Behavior of Basalt Fiber/Polypropylene Composite Rods for Passive Concrete Reinforcement
Next Article in Special Issue
Phytic Acid Doped Polypyrrole as a Mediating Layer Promoting Growth of Prussian Blue on Cotton Fibers for Solar-Driven Interfacial Water Evaporation
Previous Article in Journal
Use of Dicyclopentadiene and Methyl Dicyclopentadiene for the Synthesis of Unsaturated Polyester Resins
Previous Article in Special Issue
Mechanical Performance of Granite Fine Fly Dust-Filled Basalt/Glass Polyurethane Polymer Hybrid Composites
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

An Assessment of ASTM E1922 for Measuring the Translaminar Fracture Toughness of Laminated Polymer Matrix Composite Materials

by
Islam El-Sagheer
1,
Amr A. Abd-Elhady
1,
Hossam El-Din M. Sallam
2,* and
Soheir A. R. Naga
1
1
Mechanical Design Department, Faculty of Engineering, Helwan University, Cairo 11718, Egypt
2
Materials Engineering Department, Zagazig University, Zagazig 44519, Egypt
*
Author to whom correspondence should be addressed.
Polymers 2021, 13(18), 3129; https://doi.org/10.3390/polym13183129
Submission received: 24 August 2021 / Revised: 10 September 2021 / Accepted: 11 September 2021 / Published: 16 September 2021
(This article belongs to the Collection Reinforced Polymer Composites)

Abstract

:
The main objective of this work is to predict the exact value of the fracture toughness (KQ) of fiber-reinforced polymer (FRP). The drawback of the American Society for Testing Materials (ASTM) E1922 specimen is the lack of intact fibers behind the crack-tip as in the real case, i.e., through-thickness cracked (TTC) specimen. The novelty of this research is to overcome this deficiency by suggesting unprecedented cracked specimens, i.e., matrix cracked (MC) specimens. This MC exists in the matrix (epoxy) without cutting the glass fibers behind the crack-tip in the unidirectional laminated composite. Two different cracked specimen geometries according to ASTM E1922 and ASTM D3039 were tested. 3-D FEA was adopted to predict the damage failure and geometry correction factor of cracked specimens. The results of the TTC ASTM E1922 specimen showed that the crack initiated perpendicular to the fiber direction up to 1 mm. Failure then occurred due to crack propagation parallel to the fiber direction, i.e., notch insensitivity. As expected, the KQ of the MC ASTM D3039 specimen is higher than that of the TTC ASTM D3039 specimen. The KQ of the MC specimen with two layers is about 1.3 times that of the MC specimen with one layer.

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 ZrO2/(ZrO2 + Ni) FGMs. They used Tsai–Wu tensor theory as the failure criteria of the C/SiC plates. They concluded that ZrO2 + 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.

2. Materials and Methods

2.1. Experimental Work

2.1.1. Materials Preparation and Properties

The materials used in the conducted experiment are unidirectional GFRE composites. GFREs are composed of several layers that are bonded together to produce a multilayer composite. The hand lay-up technique is used in the present work. Firstly, glass fibers in the form of roving are impregnated with the resin matrix. The matrix is a high resistance against mechanical stresses and chemical effects and is ready to use after mixing the two components (Resin and Hardener), were supplied by CMB (Giza, Egypt). Continuous glass fibers were supplied by Jushi (Zhejiang, China), 13 μm diameter, acquire the reinforcement.
Consequently, an epoxy-impregnated layer containing continuous unidirectional glass fibers is formed. Then, numerous impregnated layers are stacked on top of one another, and all fibers are aligned in the 0 ° direction (unidirectional). Moreover, glass fibers volume fraction is 37.68%, achieved by ignition loss testing based on ASTM D3171 [37]. Table 1 and Table 2 describe the physical and mechanical properties of fibers and epoxy, respectively, used in the present work.

2.1.2. Specimens Geometry and Experimental Setup

The tested specimens are created by cutting them into pieces from the unidirectional GFRE composite. The Specimens are divided into three types based on the crack type as follows:
  • Smooth Specimen: crack-free specimen as shown in Figure 1a.
  • TTC Specimen: It has an edge crack created using a low-speed saw perpendicular to the fiber direction that left an air gap between two surfaces, as shown in Figure 1b.
  • MC Specimen: the matrix has through crack without cutting the fiber within the crack surfaces, as shown in Figure 1c. The following steps were adopted to create the MC in the unidirectional fiber laminated composites: After placing the dry fibers in each lamina in the hand lay-up technique, the region of MC was defined to put wax around the fibers and drop the Dimethylformamide (DMF), was supplied by Gama Labs (Cairo, Egypt), by syringe. The role of DMF is to stop epoxy curing at the region of MC, and that is producing a debonding region between fibers and epoxy, as shown in Figure 2.
The specimen’s geometric sizes are classified based on ASTM E1922, as shown in Figure 3, and ASTM D3039 Figure 4. The number of layers and the specimen thickness of ASTM E1922 specimens are three and 2.5 mm, respectively. However, the number of layers and specimen thickness of ASTM D3039 specimens are (two and 1.88 mm) and (one and 1.146 mm), respectively. For each GFRE composite type, five specimens are tested for each test to achieve precise results. A universal testing machine (WDW-300), was supplied by Jinan Hensgrand (Jinan, China), conducts the tension test, as shown in Figure 5. To fulfill the load rate requirement for ASTM E1922, i.e., the time from zero to peak-load ranges between 30 and 100 s, a 5 mm/sec loading rate for the ASTM E1922 specimens was used. For ASTM D3039 specimens, a standard head displacement rate of 2 mm/min was used.

2.2. Numerical Study

2.2.1. Geometry and Mechanical Properties

The present geometrical model is created as dimensions in experimental work, as shown in Figure 3 and Figure 4. The model has formed from three unidirectional layers of 0° fiber direction perpendicular to the crack surface for ASTM E1922 specimens and two layers for ASTM D3039 specimens. Equations calculate the mechanical properties of glass fiber-epoxy lamina from elasticity approach models [38] as follows:
E 11 = V F E F + V M E M + ( ( 4 V F V M ( ν F ν M ) 2 ) ( V F K M + V M K F + 1 G M ) ) ,
ν 12 = V F ν F + V M ν M + ( ( V F V M ( ν F ν M ) ( 1 K M 1 K F ) ) ( 1 G M ) ) ,
G 12 = G M ( G F ( 1 + V F ) + G M V M G F V M + G M ( 1 + V F ) ) ,
For   G 23   : A ( G 23 G M ) 2 + 2 B ( G 23 G M ) + C = 0 ,
ν 23 = K R G 23 K R G 23 ,
E 22 = 2 G 23 ( 1 + ν 23 ) ,
K F = E F 2 ( 1 2 ν F ) ( 1 + ν F ) ,
K M = E M 2 ( 1 2 ν M ) ( 1 + ν M ) ,
K = K M V M ( K F + G M ) + K F V F ( K M + G M ) V M ( K F + G M ) + V F ( K M + G M ) ,
R = 1 + 4 K ( ( ν 12 ) 2 E 11 )
V is the volume fraction, E is the modulus of elasticity, G is the shear modulus, and ν is the Poisson’s ratio. F is superscript for fiber, and M is superscript for matrix. KF, KM, and K are bulk moduli for fiber, matrix, and composite, respectively, under longitudinal strain. A, B, and C can be obtained from Ref. [38]. Finally, R is a constant.
Table 3 shows the mechanical properties of glass fiber–epoxy lamina, where E11 is the longitudinal modulus of elasticity (parallel to fibers direction). E22 and E33 are the in-plane transverse and out-of-plane moduli of elasticity, respectively. G12 and G13 are the longitudinal shear moduli in two orthogonal planes containing the fibers, and G23 is the out-of-plane shear modulus as well as ν12, ν13, and ν23 are the Poisson’s ratios.

2.2.2. Contact and Boundary Conditions

Axial displacement control was acting on the ASTM D3039 specimens and ASTM E1922 specimens. The ASTM E1922 specimen was loaded from the upper hole in the y-direction and was fixed from the lower hole. The ASTM D3039 specimen was loaded from the upper side in the y-direction and was fixed from the opposite side, as shown in Figure 6. The composite layers are considered perfectly bonded by using the tie contact option in ABAQUS/Standard [39].

2.2.3. Discrete Model

Two element types are used in the present investigation according to the Hashin damage criterion and CIM. Each element in the current model was meshed using continuum shell elements type (SC8R), which confirmed with Hashin damage criterion, and an eight-node linear brick (C3D8R) confirmed with contour integral. The mesh refinement procedure was carried out to ensure that the element’s size does not affect the findings. The composite specimen model’s mesh is shown in Figure 6.

2.2.4. Hashin Damage Model

In the present work, Hashin’s failure criteria are used to predict the translaminar progressive fracture and get the failure load of the present composite specimen according to ASTM D3039 and ASTM E1922. In general, Hashin’s failure criteria are used to predict anisotropic damage in elastic-brittle materials. It is primarily used with FRP and considers four failure modes, namely: fiber breakage in fiber tension, fiber buckling in fiber compression, matrix cracking in matrix tension, and matrix crushing in matrix compression. The initiation criteria have the following general forms:
Fiber   tension ( σ 11   0 ) : F f t = ( σ 11 X T ) 2 + α ( τ 12 S L ) 2 ,  
Fiber   compression ( σ 11 < 0 ) : F f c = ( σ 11 X C ) 2 ,
Matrix   tension   ( σ 22 0 ) : F m t = ( σ 22 Y T ) 2 + ( τ 12 S L ) 2 ,
Matrix   compression   ( σ 22 < 0 ) : F m c = ( σ 22 2 S T ) 2 + [ ( Y C 2 S T ) 2 1 ] σ 22 Y C + ( τ 12 S T ) 2
In the above equations, X T denotes the longitudinal tensile strength,   X C denotes the longitudinal compressive strength,   Y T denotes the transverse tensile strength,   Y C denotes the transverse compressive strength. Furthermore,   S L denotes the longitudinal shear strength, and S T denotes the transverse shear strength. Those can be obtained from Refs. [40,41] and are tabulated in Table 4. α is a coefficient that determines the contribution of the shear stress to the fiber tensile initiation criterion. The σ 11 ,   σ 22 ,   and   τ 12 components of the effective stress tensor σ. If one of these variables ( F F T , F F C , F M T and F M C ) reaches the unity, one of the previous failure modes is achieved and the damage process in the elasticity matrix is given by the following equations:
[ σ 11 σ 22 τ 12 ] = C d [ ε 11 ε 22 ε 12 ] = 1 D [ ( 1 d f ) E 11 ( 1 d f ) ( 1 d m ) ν 21 E 11 0 ( 1 d f ) ( 1 d m ) ν 12 E 22 ( 1 d m ) E 22 0 0 0 ( 1 d s ) G 12 D ] [ ε 11 ε 22 ε 12 ] ,
where Cd is the damaged elasticity matrix, and
D = 1 ( 1 d f ) ( 1 d m ) ν 12 ν 21 ,
d f = { d f t d f c , , i f i f ( σ 11 0 ) ( σ 11 < 0 ) ,
d m = { d m t d m c , , i f i f ( σ 22 0 ) ( σ 22 < 0 ) ,
d s = 1 ( 1 d f t ) ( 1 d f c ) ( 1 d m t ) ( 1 d m c )
where df, dm, and ds are the damage variables for fibers, matrix, and shear, and they can be obtained from Ref. [28]. The initial values of df, dm, and ds are assumed to be zero and gradually increased until they reach the unity. According to equation 15, as the load increases, the effective strength of the elements decreases. Elements with a damage variable equal to the unity are excluded because they can no longer withstand such stress.

2.2.5. Contour Integral Method

Contour Integral provides a standard scheme of finite elements, usually helping to match the mesh to the cracked frame, specifically describing the crack front Assessment can be viewed as a simulated displacement of a material block across the crack tip in 3D, represented by contours where each contour is a circle of fully surrendered components. Each contour gives an integral assessment of the contours. In the present work, CIM was used to obtain the stress intensity factor (SIF) in mode I (KI) for GFREs for TTC and MC specimens.
In the TTC specimen, the crack cuts the matrix and fiber as one bulk material without slippage between the matrix and the fiber. And the CIM can be used to get the SIF directly. Otherwise, for MC-specimen, the crack at matrix only; therefore, the matrix is slipping on the fiber during the tensile load. A unique and simple method is used to simulate the MC specimen to get the SIF. As is already known, the bond strength between fiber and matrix can be classified into three cases. Namely, perfect bond, i.e., no fiber bridging, and the composite failed as a monolithic material (unpreferable case), strong bond (preferable case), i.e., fiber bridging achieved with different efficiency depending on the debonding length of the fiber (Notch sensitivity achieved), as shown in Figure 7.
The last case is a weak bond, i.e., the crack grows parallel to the applied load and fiber direction (notch sensitivity does not achieve). MC model is utilized to calculate KI under debonding zones between fibers and epoxy. To simulate the fiber bridging simply, the debonding zones translated into degrees as θ is debonding angle in degree between fibers and matrix in MC model. Figure 7 shows a simulation of the deboning zone between fibers and matrix. The fibers in the debonding zone act as a spring with efficiency depending on its stiffness (Modulus of elasticity of fibers) and length (debonding length of fibers). Different values of θ are adopted; namely, 0, 10, 20, 30, 40, 50, 60, 70, 80, 90, 100, 110, and 120, to study the effect of debonding length on the crack closure, i.e., the reduction in SIF.

3. Results and Discussions

3.1. Progressive Failure

3.1.1. Progressive Failure for ASTM E1922 Specimen

Figure 8 shows the photographs of the progressive failure for ASTM E1922 specimens under the tensile load through the two holes. For the TTC specimen, the crack started its propagation perpendicular to the direction of the fibers up to about 1 mm, as shown in Figure 8a. By definition, the fracture toughness can be calculated at the crack initiation load. After that, this propagating crack changed its direction to be parallel to the applied load and fibers direction until final failure, as shown in Figure 8a. This observation is confirmed with results obtained by Ref. [17]. On the other hand, smooth specimen and MC specimen are failed through the loaded holes instead of the middle of cross-section, as illustrated in Figure 8b. It can be concluded that ASTM E1922 only works with TTC specimens and does not work with MC specimens. To avoid failure at the point of loadings in the ASTM E1922 specimens, MC specimens should be subjected to uniform tensile loads as followed in the tensile test in ASTM D3039. Therefore, the fracture toughness of fibrous composite materials has been measured through TTC and MC ASTM D3039 specimens.

3.1.2. Progressive Failure for ASTM D3039 Specimen

Based on ASTM D3039, a flat specimen setup is conducted experimentally for smooth and cracked specimens to examine the tensile and fracture behaviors of GFRE material. The dimensions of TTC and MC specimens followed ASTM D3039 with a/W = 0.5. The progressive damage of the smooth specimen under the tensile load can be depicted in Figure 9, divided into major steps in both fibers and matrix damage. In the first step at displacement = 5.531 mm, a longitudinal crack propagates in one side of the specimen, the matrix has started to be damaged, and fibers are sustaining the load. In the second step, at load 6.914 mm, two longitudinal cracks propagate on opposite sides, the matrix has partly been damaged, and fibers are standstill sustaining the load. Finally, at load-displacement 7.202 mm, the matrix has been damaged and splitting from the fibers, while fibers are partly damaged.
Figure 10 shows several major steps in both fibers and matrix damages in the case of TTC specimens. Firstly, at load-displacement at 1.63 mm, a longitudinal crack is propagating besides the crack-tip side, the matrix has started to be damaged, and fibers are sustaining the load. At applied displacement at 4.61 mm, a longitudinal crack is still propagating toward grip, and the matrix has partly been damaged. The fibers partly sustain the load at applied displacement 5.33 mm, and the matrix has been damaged, splitting from the fibers, while fibers are partly damaged.
Hashin’s failure criteria were executed in the present work to predict the fiber and matrix damages due to tensile or compression stress, as shown in Figure 9 and Figure 10. From Figure 9, it was observed that the matrix is beginning damage in the specimen’s middle after that; the debonding between the matrix and fiber has happened. This debonding and matrix damage grows; finally, the fibers are damaged in the upper and lower part of the specimen. From Figure 10, it can be seen that the damage of matrix and fiber are beginning from the pre-crack tip and is growing in the fiber direction parallel to the applied load that confirms with [17,42]. From Figure 9 and Figure 10, the Hanshin model can simulate the progressive damage of smooth and TTC ASTM D3039 specimens.
From Figure 9 and Figure 10 (experimental and numerical results), it can be concluded that the matrix in the ASTM D3039 smooth and TTC specimens are completely damaged and split from the fibers into all layers. After that, the fibers sustain the load, and the main damage type concerns the matrix in tension load. That is why the FRP matrix carries the major load.
The main conclusion from the final failure shape of the TTC specimen of ASTM E1922 and ASTM D3039 is that the ASTM E1922 is completely fractured into two pieces, while in the ASTM D3039, the specimen is completely damage but not separate.
As shown in Figure 11, the debonding between the epoxy-matrix and glass-fiber in the final failure shape of MC specimen of ASTM D3039 is much more observed than smooth and TTC specimens. Furthermore, the pre-crack in the TTC specimen becomes blunting, and the debonding between the fiber and matrix begins from the crack tip then grows perpendicular to the pre-crack.

3.2. Load Displacement Curve

Figure 12 shows the load-displacement curve of ASTM E1922 and ASTM D3039 specimens under tensile load. It can be observed that the applied load increase to reach a peak value which the matrix is completely damaged (as shown in Figure 8, Figure 9 and Figure 10), then the applied load is suddenly dropped. As expected, the applied load peak value for the smooth specimen is more than the MC and TTC specimens of ASTM D3039. Moreover, the value of the peak load of the TTC specimens is lower than that of MC specimens. For the same a/W, the peak value of the applied load of the TTC specimen of ASTM D3039 (7.11 kN) is much higher than ASTM E1922 (2.24 kN), although the thickness of the specimen of ASTM E1922 (2.5 mm) is higher than that of ASTM D3039 (1.88 mm). This difference may be attributed to the applied load being eccentric in ASTM E1922, while the applied load is uniform distributed in ASTM D3039.
The numerical results of the Hashin model can be verified by comparing it with the experimental result for smooth and TTC ASTM D3039 specimens in Table 5. Table 5 shows a good agreement between the peak applied load values of smooth and TTC specimens, computed from the numerical and experimental methods. Therefore, Hashin numerical method can be used to determine the damage values of composite laminates. It is worth noting that the load at which the crack initiated from the notch root of the TTC ASTM E1922 specimen is indicated in Figure 12.
The thickness effect was also studied by changing the number of layers on the tensile and fracture behavior of smooth, MC, and TTC specimens of ASTM D3039. Two specimen thicknesses are 1.146 mm and 1.88 mm that formed using one and two layers, respectively. The thickness or the number of layers does not affect the mode of failure of the smooth, MC, and TTC specimens of ASTM D3039. The pop-in events can be expressed as the damage growth in the matrix, showing a temporary drop in load, and then the fibers sustain the load again before the matrix is completely damaged. Figure 13a. shows that the modulus of elasticity (E) and the ultimate tensile strength (UTS) of unidirectional GFRE of two layers (1.88 mm) are higher than those of GFRE of on layer (1.46 mm). In the case of two layers, the initial tangent modulus Ein is about 33 GPa, and this slope decreased at stress equals about 150 MPa to be 20 GPa, while the UTS equals 466 MPa at strain equals 0.02. In one layer, the Ein is about 19 GPa, and this slope decreased at stress equals about 270 MPa to be 10 GPa, while the UTS equals 409 MPa at strain equals 0.027. Although GFRE with two layers has higher mechanical properties, the deviation of E in GFRE with one layer occurs at stress higher than that of GFRE with two layers. Furthermore, the maximum stress of GFRE with two layers is higher than that of GFRE with one layer either for MC specimen or for TTC specimen, as shown in Figure 13b,c, respectively. It is worth noting that the load at which the crack initiated from the notch root is indicated in Figure 13b,c.

3.3. The stress Concentrator Factor (SCF)

As mentioned before, the unidirectional GFRE experienced initial notch sensitivity behavior followed by notch insensitivity. Therefore, the cracked specimens of ASTM D3039 were analyzed based on the first principle of mechanics of materials, i.e., SCF, as listed in Table 6. The experimental SCF (KtExperimental) is calculated as the ratio of the peak load of the smooth specimen divided by the peak load of the notched (cracked) specimen, while the analytical SCF (KtAnalytical) is calculated based on the following Equation:
K t Analytical = ( 1 + 2 D ρ )
D (7.5 mm) and ρ (0.5) are the depth and radius of the notch, respectively.
It is clear from Table 6 that the coefficient of variation (CV) of smooth specimens is lower than that of TTC specimens. This difference may be attributed to the precision of notch radius manufacturing. KtExperimental for the one-layer specimen is higher than that of two-layer specimens, while the ratios of both values to those calculated analytically, KtAnalytical, are lower than the quarter. Based on these ratios, it can be concluded that the study of the notch sensitivity by the SCF is acceptable.
Furthermore, KtExperimental is around two (2.12 and 1.91 for both cases), and the ratio of notch depth to specimen width is 0.5, i.e., the net stress of notched specimen is twice its gross stress. Comparing the value of this ratio with the value of KtExperimental reveals that both values are more or less the same. The net stress analysis is evidence to prove the notch insensitivity of the present unidirectional laminated GFRE.

3.4. Fracture Toughness

3.4.1. The Stress Intensity Factor for TTC Specimen

The present work uses the CIM to calculate the SIF for TTC of ASTM E1922 and ASTM D3039. In the beginning, it must be verifying the numerical results of SIF values. To verify the present numerical results, a comparison between the SIF values extracted from the CIM and the SIF and geometry correction factor which extracted from the following analytical Equation:
For TTC ASTM E1922 specimens:
K 1 = 1 ( 1 ( a w ) ) 3 / 2 ( ( P ( B W 1 / 2 ) ) ( 1.4 + ( a w ) ) ( a w ) 1 / 2 ( 3.97 10.88 ( a w ) + 26.25 ( a w ) 2 38.9 ( a w ) 3 + 30.15 ( a w ) 4 9.27 ( a w ) 5 ) ) ,
where:
Y = K I σ   π a = 1 π   [ 1 . 4 + α ] [ 3 . 97     10 . 88 α + 26 . 25 α 2     38 . 9 α 3 + 30 . 15 α 4     9 . 27 α 5 ] / [ 1     α ] 1 . 5
where KI is SIF in (MPa·√m), P is Load in (MN), B is the thickness of the specimen in (m), W is the width of the specimen in (m), a is the crack length in (m), Y is geometry correction factor in (dimensionless), and α is (a/W).
Figure 14 illustrates the comparison between the theoretical results of SIF and Y, which are computed from Equations (21) and (22) of TTC ASTM E1922 for different crack length ratios to confirm that the model works well (a/W = 0, 0.2, 0.3, 0.4. 0.5, and 0.6) and the corresponding numerical results. From Figure 14, Table 7 and Table 8, it can be seen that there is a good agreement between the numerical and analytical results. Table 7 shows the absolute error percentage for each crack length, and the absolute error values are less than 4%.
For TTC ASTM D3039 flat specimens, the following Equation from Ref. [43] can be used:
K I = Y σ π a ,
where:
Y = 1 π ( 1.99 0.41 ( a w ) + 18.7 ( a w ) 2 38.48 ( a w ) 3 + 53.85 ( a w ) 4 )
The numerical model using the CIM is shown in Table 9. As illustrated in Table 9, the error is less than 0.9%, indicating a good agreement between the theoretical and numerical results of the SIF of the TTC ASTM D3039 specimen.

3.4.2. The Stress Intensity Factor for MC Specimen

To the best of the authors’ knowledge, no analytical formula in the literature can predict the SIF or Y for ASTM E1922 and ASTM D3039 MC-specimen. The contour integral model modifications are made to calculate Y for MC specimens for a different debonding length of fibers due to the difficulty of predicting this length. Figure 15 shows the relation between the Y at a/W = 0.5 against the debonding angle (θ) for ASTM E1922 and ASTM D3039. In general, Y and subsequently SIF increased by increasing θ, i.e., increasing the debonding length of the fibers. As it is already known, the efficiency of fiber bridging decreased by increasing this length, see Figure 7. The results show that the fibers bridging phenomenon works very well to exercise crack closure at the small values θ, i.e., the debonding length in the intact fibers behind the crack tip is short. At the large values θ, the closing of the fibers behind the crack tip decreased due to the increase in debonding length, and subsequently, the Y increased rapidly.
The lower bound of Y (<<1) at θ = 0° does not take place (perfect bond), and the composite failed as monolithic materials. This hypothesis is supported by the experimental results of Mode-I fracture toughness of Carbon-Carbon composites [44]. Kumar et al. [44] concluded that their fractured surface clearly showed brittle fracture behavior, and the fracture may be explained by energy-consuming mechanisms like fiber pull-out and fiber debonding. However, this brittle fracture behavior is evidence for the absence of fiber bridging. The fiber bridging is the main reason for converting the brittle failure of brittle composites to ductile failure. As mentioned before, the good bond between the fiber and the matrix is favorable to ensure the ductile failure, not perfect bond, not weak bond; both cause brittle failure. Figure 15 shows that Y becomes higher than the unity at a certain value θ, depending on the specimen type. The values of Y below the unity mean that the SIFs of MC cracked specimens are lower than that of the corresponding TTC specimens with infinite width. The upper bound Y (the effect of fiber bridging is diminished) must not exceed the corresponding value of Y for the TTC specimen, as stated in Figure 15. For a deeper understanding, supposing MC and TTC specimens with the same geometry and a/W are subjected to the same stress. According to Equation (23) and the values of Y in Table 10, the SIF of the MC specimen is much lower than that of the TTC specimen due to the crack closure achieved by fiber bridging. The closing SIF due to fiber bridging can be calculated by subtracting the SIF of the MC specimen from the SIF of the TTC specimen (SIFFB = SIFTTC –SIFMC). Furthermore, the closing SIF decreased by increasing θ, and it is equal to zero when the fiber bridging diminished.

3.4.3. The Apparent Fracture Toughness

The apparent fracture toughness (KQ) of TTC and MC ASTM E1922 and ASTM D3039 specimens are calculated based on the values of crack initiation loads indicated in Figure 12 and Figure 13 and listed in Table 10. Equations (21)–(24) are used to predict the values of apparent fracture toughness of TTC ASTM E1922 and ASTM D3039 specimens, while Figure 15 is used to predict the values of Y for MC ASTM D3039 specimens and substituted these values in Equation (23). It is worth noting that the apparent fracture toughness of MC ASTM D3039 specimens must be calculated using the appropriate value θ (θ = 80° in the present work to be Y of MC specimen is equal to 0.62 times that of TTC specimen, see Table 10. As already known, the crack initiation load decreases with increasing θ; thus, it is wrong to study the change in the θ values with the same crack initiation load because it would give a misleading trend that the apparent fracture toughness increases with increasing θ.
As expected, the KQ of TTC ASTM D3039 is lower than that of MC ASTM D3039, although Y of the MC specimen is much lower than that of the TTC specimen. It is clear from Table 10 that the KQ measured based on TTC ASTM D3039 is about 1.45 times that of TTC ASTM E1922. Furthermore, the KQ of TTC ASTM D3039 is not affected by the number of layers, while the KQ of MC ASTM D3039 is affected by the number of layers. The KQ of the MC ASTM D3039 specimen with two layers is about 1.3 times that of the MC ASTM D3039 specimen with one layer so that the ratio of the MC ASTM D3039 specimen with two layers to the corresponding TTC ASTM D3039 specimen (83.33/51.13 = 1.62) is higher than the ratio for one-layer specimens (62.72/50.7 = 1.24).

3.5. The Crack Mouth Opening Displacement, CMOD

Figure 16 shows the crack mouth opening displacement, CMOD, versus the applied load for TTC specimen of ASTM E1922 and ASTM D3039. The values of CMOD increase by increasing the applied load to reach the peak value, then it drops; this trend is similar for the two specimens. From Figure 16, it can be seen that the shapes of the curves of applied load versus CMOD and the applied load versus displacement are similar.

4. Conclusions

Based on the present experimental and numerical results, it can be concluded that the TTC ASTM E1922 specimen is not an appropriate candidate to get the fracture toughness of laminated fibrous composite materials. The suggested MC specimen is a good candidate to predict the actual fracture toughness of such materials. The MC ASTM E1922 specimen failed to initiate the crack from the root of the pre-notch, and the failure occurred at the point of loading due to the presence of stress concentration. Single edge cracked specimen with a/W = 0.5 under uniform loading is suggested to predict the fracture toughness. The dimensions of this specimen followed ASTM D3039 specimen; both TTC and MC ASTM D3039 specimens show notch sensitivity, i.e., matrix crack growth perpendicular to the fiber direction up to crack length equals about 1 mm. After that, the crack kinks to grow parallel to the fiber direction, i.e., notch insensitivity. It is suggested that the fracture toughness might be calculated before the occurrence of kinked crack.
The present experimental and numerical results reveal several points of interest: (1) The SCF is a good candidate to predict notch sensitivity in such composites. (2) The KQ measured based on TTC ASTM D3039 is about 1.45 times that of TTC ASTM E1922. (3) In the case of MC specimens, the numerical results show that Y and subsequently SIF increased by increasing the debonding length of the fibers. (4) As expected, the ratio of the fracture toughness measured from the MC specimen to that measured from the TTC specimen is greater than the unity, ranging from 1.24 in the one-layer specimen to 1.62 in the two-layer specimen. (5) The KQ measured from MC ASTM D3039 increased with increasing the number of layers.

Author Contributions

Conceptualization, H.E.-D.M.S. and S.A.R.N.; methodology, I.E.-S., H.S and A.A.A.-E.; software, I.E.-S.; validation, I.E.-S. and A.A.A.-E.; formal analysis, I.E.-S. and A.A.A.-E.; investigation, I.E.-S. and A.A.A.-E.; resources, I.E.-S. and H.E.-D.M.S.; writing—original draft preparation, A.A.A.-E. and I.E.-S.; writing—review and editing, H.E.-D.M.S. and S.A.R.N.; visualization, A.A.A.-E. and S.A.R.N.; supervision, S.A.R.N., A.A and H.E.-D.M.S. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

The data presented in this study are available on request from the Corresponding author.

Conflicts of Interest

The authors declare no conflict of interest.

Nomenclature

3D FEAThree-Dimensional finite element analysis
ASTMthe American Society for Testing Materials
CCDCharged Coupled Camera Device
CIMContour Integral Method
CVcoefficient of variation
CZMCohesive Zone Model
DMFDimethylformamide
FEAFinite Element Analysis
FEMFinite Element method
FRPFiber Reinforced Polymer
GFREGlass Fiber Reinforced Epoxy
GFRPGlass Fiber Reinforced Polymer
HDMHashin Damage Model
HSNFTCRTHashin Fiber Tension Criterion
HSNMTCRTHashin Matrix Tension Criterion
ITEIntrinsic To Extrinsic
LEFMLinear Elastic Fracture Mechanics
MCMatrix Crack
MTPMModified Jenq-Shah’s two-parameter model
SCFStress Concentration Factor
SIFStress intensity Factor
TTCThrough Thickness Crack
UTSUltimate Tensile Strength
XFEMeXtended Finite Element Method
aCrack Length
BSpecimen Thickness
DNotch Depth
EModulus of Elasticity
EinInitial Tangent Modulus of Elasticity
Eσ= -- MPaTangent Modulus of Elasticity at certain stress at which the modulus will have deviated
FSuperscript for Fiber
Ff tFiber Tension Criterion
Ff cFiber Compression Criterion
Fm tMatrix Tension Criterion
FmcMatrix Compression Criterion
GIn-plane Shear Modulus
KIMode I SIF
KQThe apparent fracture toughness
KtSCF
KtAnalyticalAnalytical SCF
KtExperimentalExperimental SCF
MSuperscript for Matrix
PApplied Load
VFiber Volume Fraction
WSpecimen Width
SLThe Longitudinal Shear Strength
STThe Transverse Shear Strength
XCThe Longitudinal Compressive Strength
XTThe Longitudinal Tensile Strength
YGeometry Correction Factor
YCThe Transverse Compressive Strength
YTThe Transverse Tensile Strength
αThe Contribution Coefficient of The Shear Stress to The Fiber Tensile Initiation Criterion
νPoisson’s Ratio
ρNotch radius

References

  1. Teixeira, R.; Pinho, S.; Robinson, P. Translaminar ply fracture toughness of advanced composites. In Proceedings of the ICCM—18th International Conference on Composite Materials, ICC, Jeju, Korea, 21–26 August 2011. [Google Scholar]
  2. Teixeira, R.D.F.; Pinho, S.; Robinson, P. Translaminar fracture toughness of CFRP: From the toughness of individual plies to the toughness of the laminate. In Proceedings of the ECCM-15th European Conference on Composite Materials, Venice, Italy, 24–28 June 2012. [Google Scholar]
  3. Mohsin, M.; Iannucci, L.; Greenhalgh, E. Translaminar Fracture Toughness Characterisation of Carbon Fibre Reinforced Thermoplastic Composites. In Proceedings of the ECCM-17th European Conference on Composite Materials, Munich, Germany, 26–30 June 2016. [Google Scholar]
  4. Abdullah, S.S.; Iannucci, L.; Greenhalgh, E. On the translaminar fracture toughness of Vectran/epoxy composite material. Compos. Struct. 2018, 202, 566–577. [Google Scholar] [CrossRef]
  5. Gigliotti, L.; Pinho, S. Translaminar fracture toughness of NCF composites with multiaxial blankets. Mater. Des. 2016, 94, 410–416. [Google Scholar] [CrossRef] [Green Version]
  6. Haldar, S.; Herraez, M.; Naya, F.; Gonzalez, C.; Lopes, C.S. Relations between intralaminar micromechanisms and translaminar fracture behavior of unidirectional FRP supported by experimental micromechanics. Compos. Part B Eng. 2019, 174, 107000. [Google Scholar] [CrossRef]
  7. Souza, R.R.; Nascimento Junior, S.L.; Silveira, N.N.; Arbelo, M.A.; Donadon, M.V. Translaminar fracture toughness and fatigue crack growth characterization of carbon-epoxy plain weave laminates. Polym. Compos. 2019, 40, 3791–3804. [Google Scholar] [CrossRef]
  8. Vantadori, S.; Carpinteri, A.; Głowacka, K.; Greco, F.; Osiecki, T.; Ronchei, C.; Zanichelli, A. Fracture toughness characterisation of a glass fibre-reinforced plastic composite. Fatigue Fract. Eng. Mater. Struct. 2021, 44, 3–13. [Google Scholar] [CrossRef]
  9. Jenq, Y.; Shah, S. Two-parameter fracture model for concrete. J. Eng. Mech. 1985, 111, 1227–1241. [Google Scholar] [CrossRef]
  10. RILEM Technical Committee, 89-FMT. Determination of fracture parameters (KsIC and CTODc) of plain concrete using three-point bend tests, proposed RILEM draft recommendations. Mater. Struct. 1990, 23, 457–460. [Google Scholar] [CrossRef]
  11. Siddique, A.; Abid, S.; Shafiq, F.; Nawab, Y.; Wang, H.; Shi, B.; Saleemi, S.; Sun, B. Mode I fracture toughness of fiber-reinforced polymer composites: A review. J. Indust. Text. 2021, 50, 1165–1192. [Google Scholar] [CrossRef]
  12. Jia, Y.; Wang, H.; Liu, B.; Huang, Y.; Gao, H. Intrinsic-to-extrinsic transition in fracture toughness through structural design: A lesson from nature. Extrem. Mech. Lett. 2020, 37, 100685. [Google Scholar] [CrossRef] [Green Version]
  13. El-Hajjar, R.; Haj-Ali, R. Mode-I fracture toughness testing of thick section FRP composites using the ESE (T) specimen. Eng. Fract. Mech. 2005, 72, 631–643. [Google Scholar] [CrossRef]
  14. ASTM. ASTM E1922-04. Standard Test Method for Translaminar Fracture Toughness of Laminated and Pultruded Polymer Matrix Composite Materials; ASTM: West Conshohocken, PA, USA, 2015. [Google Scholar]
  15. Laffan, M.J. Development of Translaminar Fracture Toughness Testing Methods for Composite Materials. Ph.D. Thesis, Department of Aeronautics Imperial College, London, UK, 2011. [Google Scholar]
  16. Laffan, M.; Pinho, S.; Robinson, P.; McMillan, A. Translaminar fracture toughness testing of composites: A review. Polym. Test. 2012, 31, 481–489. [Google Scholar] [CrossRef]
  17. Saadati, Y.; Lebrun, G.; Bouvet, C.; Chatelain, J.-F.; Beauchamp, Y. Study of translaminar fracture toughness of unidirectional flax/epoxy composite. Compos. Part C Open Access 2020, 1, 100008. [Google Scholar] [CrossRef]
  18. He, M.-Y.; Wu, B.; Suo, Z. Notch-sensitivity and shear bands in brittle matrix composites. Acta Metal. Mater. 1994, 42, 3065–3070. [Google Scholar] [CrossRef]
  19. Zappalorto, M. Static notch sensitivity in orthotropic materials and composites. Eur. J. Mech.-Solids 2021, 85, 104094. [Google Scholar] [CrossRef]
  20. Jamali, J.; Mourad, A.I.; Fan, Y.; Wood, J. Through-thickness fracture behavior of unidirectional glass fibers/epoxy composites under various in-plane loading using the CTS test. Eng. Fract. Mech. 2016, 156, 83–95. [Google Scholar] [CrossRef]
  21. ASTM. ASTM D3039/D3039M-17. Standard Test Method for Tensile Properties of Polymer Matrix Composite Materials; ASTM: West Conshohocken, PA, USA, 2015. [Google Scholar]
  22. Elbadry, E.A.; Aly-Hassan, M.S.; Hamada, H. In Tensile and bending properties of jute fiber mat reinforced unsaturated polyester matrix composites. ASME Int. Mech. Eng. Cong. Exposit. 2011, 3, 751–759. [Google Scholar]
  23. El-Wazery, M.; El-Elamy, M.; Zoalfakar, S. Mechanical properties of glass fiber reinforced polyester composites. Int. J. Appl. Sci. Eng. 2017, 14, 121–131. [Google Scholar]
  24. Hashin, Z. Failure criteria for unidirectional fiber composites. J. Appl. Mech. 1980, 47, 329. [Google Scholar] [CrossRef]
  25. Hashin, Z.; Rotem, A. A fatigue failure criterion for fiber reinforced materials. J. Compos. Mater. 1973, 7, 448–464. [Google Scholar] [CrossRef] [Green Version]
  26. Pham, D.C.; Cui, X.; Ren, X.; Lua, J. A discrete crack informed 3D continuum damage model and its application for delamination migration in composite laminates. Compos. Part B Eng. 2019, 165, 554–562. [Google Scholar] [CrossRef]
  27. Koloor, S.S.; Karimzadeh, A.; Yidris, N.; Petrů, M.; Ayatollahi, M.R.; Tamin, M.N. An energy-based concept for yielding of multidirectional FRP composite structures using a mesoscale lamina damage model. Polymers 2020, 12, 157. [Google Scholar] [CrossRef] [Green Version]
  28. Duarte, A.; Sáez, A.D.; Silvestre, N. Comparative study between XFEM and Hashin damage criterion applied to failure of composites. Thin-Walled Struct. 2017, 115, 277–288. [Google Scholar] [CrossRef]
  29. El-Sagheer, I.; Taimour, M.; Mobtasem, M.; Abd-Elhady, A.; Sallam, H.E.-D.M. Finite Element analysis of the behavior of bonded composite patches repair in aircraft structures. Frat. Integ. Strut. 2020, 14, 128–138. [Google Scholar]
  30. Abd-Elhady, A.A. Effect of location and dimensions of welded cover plate on stress intensity factors of cracked plates. Ain Shams Eng. J. 2013, 4, 863–867. [Google Scholar] [CrossRef] [Green Version]
  31. Abd-Elhady, A.A.; Sallam, H.E.-D.M.; Alarifi, I.M.; Malik, R.A.; El-Bagory, T.M. Investigation of fatigue crack propagation in steel pipeline repaired by glass fiber reinforced polymer. Compos. Struct. 2020, 242, 112189. [Google Scholar] [CrossRef]
  32. Carpenter, K.; Lei, Y.; Asadi, A.; Parmigiani, J.; Tabei, A. Applicability of linear elastic fracture mechanics to compressive damage in carbon fiber reinforced epoxy matrix composites. Mech. Adv. Mater. Struct. 2021. [Google Scholar] [CrossRef]
  33. Liu, P.; Gu, Z.; Peng, X.; Zheng, J. Finite element analysis of the influence of cohesive law parameters on the multiple delamination behaviors of composites under compression. Compos. Struct. 2015, 131, 975–986. [Google Scholar] [CrossRef]
  34. Rozylo, P. Experimental-numerical study into the stability and failure of compressed thin-walled composite profiles using progressive failure analysis and cohesive zone model. Compos. Struct. 2021, 257, 113303. [Google Scholar] [CrossRef]
  35. Panettieri, E.; Fanteria, D.; Danzi, F. Delaminations growth in compression after impact test simulations: Influence of cohesive elements parameters on numerical results. Compos. Struct. 2016, 137, 140–147. [Google Scholar] [CrossRef]
  36. Zhou, W.; Zhang, R.; Fang, D. Design and analysis of the porous ZrO2/(ZrO2+ Ni) ceramic joint with load bearing–heat insulation integration. Ceram. Int. 2016, 42, 1416–1424. [Google Scholar] [CrossRef] [Green Version]
  37. ASTM. ASTM D3171-15. Standard Test Methods for Constituent Content of Composite Materials; ASTM: West Conshohocken, PA, USA, 2015. [Google Scholar]
  38. Systèmes Dassault. ABAQUS 6.14 analysis User’s Manual; Dassault Systems: Providence, RI, USA, 2014. [Google Scholar]
  39. Younes, R.; Hallal, A.; Fardoun, F.; Chehade, F.H. Comparative review study on elastic properties modeling for unidirectional composite materials. Compos. Propert. 2012, 17, 391–408. [Google Scholar]
  40. Lupasteanu, V.; Taranu, N.; Popoaei, S. Theoretical strength properties of unidirectional reinforced fiber reinforced polymer composites. Buletinul Institutului Politehnic din lasi. Sectia Constr. Arhitect. 2013, 59, 83. [Google Scholar]
  41. Atta, M.; Abu-Sinna, A.; Mousa, S.; Sallam, H.; Abd-Elhady, A. Flexural behavior of functionally graded polymeric composite beams. J. Indust. Text. 2021. [Google Scholar] [CrossRef]
  42. Abd-Elhady, A.; Mubaraki, M.A.; Sallam, H. Progressive failure prediction of pinned joint in quasi-isotropic laminates used in pipelines. Latin Am. J. Solids Struct. 2018, 15, e96. [Google Scholar] [CrossRef]
  43. David, B. Elementary Engineering Fracture Mechanics, 4th ed.; Martinus Nijhoff: Boston, MA, USA, 1982; pp. 70–77. [Google Scholar]
  44. Kumar, B.V.S.; Londe, V.N.; Lokesha, M.; Anilas, M.; Surendranathan, A.O. Experimental investigation on Mode-I fracture toughness of Carbon-Carbon composites fabricated by preformed yarn method. Mater. Today Proc. 2021, 46, 2431–2435. [Google Scholar] [CrossRef]
Figure 1. Types of cracks, (a) smooth (crack free), (b) through thickness crack, and (c) matrix crack.
Figure 1. Types of cracks, (a) smooth (crack free), (b) through thickness crack, and (c) matrix crack.
Polymers 13 03129 g001
Figure 2. Preparation of matrix crack by wax and DMF.
Figure 2. Preparation of matrix crack by wax and DMF.
Polymers 13 03129 g002
Figure 3. The geometry of ASTM E1922 specimens, (a) smooth specimen, and (b) TTC Specimen. All dimensions in mm.
Figure 3. The geometry of ASTM E1922 specimens, (a) smooth specimen, and (b) TTC Specimen. All dimensions in mm.
Polymers 13 03129 g003
Figure 4. The geometry of ASTM D3039 specimens, (a) smooth specimen, (b) TTC specimen, and (c) MC specimen. All dimensions in mm.
Figure 4. The geometry of ASTM D3039 specimens, (a) smooth specimen, (b) TTC specimen, and (c) MC specimen. All dimensions in mm.
Polymers 13 03129 g004
Figure 5. Experimental Setup.
Figure 5. Experimental Setup.
Polymers 13 03129 g005
Figure 6. Mesh and boundary conditions of the present model, (a) ASTM E1922, and (b) ASTM D3039.
Figure 6. Mesh and boundary conditions of the present model, (a) ASTM E1922, and (b) ASTM D3039.
Polymers 13 03129 g006
Figure 7. A schematic diagram showing the simulation of debonding zone between fibers and matrix in MC Model.
Figure 7. A schematic diagram showing the simulation of debonding zone between fibers and matrix in MC Model.
Polymers 13 03129 g007
Figure 8. Photographs showing the progressive damage for (a) TTC ASTM E1922 specimen and (b) MC ASTM E1922 specimen.
Figure 8. Photographs showing the progressive damage for (a) TTC ASTM E1922 specimen and (b) MC ASTM E1922 specimen.
Polymers 13 03129 g008
Figure 9. Numerical and experimental progressive damage failure of ASTM D3039 specimen.
Figure 9. Numerical and experimental progressive damage failure of ASTM D3039 specimen.
Polymers 13 03129 g009
Figure 10. Numerical and experimental progressive damage failure of TTC ASTM D3039 specimen.
Figure 10. Numerical and experimental progressive damage failure of TTC ASTM D3039 specimen.
Polymers 13 03129 g010
Figure 11. Close-up view of experimental final failure shape of ASTM D3039 for (a) Smooth specimen, (b) MC specimen, and (c) TTC specimen.
Figure 11. Close-up view of experimental final failure shape of ASTM D3039 for (a) Smooth specimen, (b) MC specimen, and (c) TTC specimen.
Polymers 13 03129 g011
Figure 12. Load-Displacement diagram for ASTM E1922 and ASTM D3039 specimens.
Figure 12. Load-Displacement diagram for ASTM E1922 and ASTM D3039 specimens.
Polymers 13 03129 g012
Figure 13. Specimens of ASTM D3039 made of GFRE with one layer and two layers: (a) Stress-Strain curve of the smooth specimen, and Stress-Displacement diagram of (b) MC specimen, and (c) TTC specimen.
Figure 13. Specimens of ASTM D3039 made of GFRE with one layer and two layers: (a) Stress-Strain curve of the smooth specimen, and Stress-Displacement diagram of (b) MC specimen, and (c) TTC specimen.
Polymers 13 03129 g013
Figure 14. Comparison between the theoretical and numerical results of SIF for TTC ASTM E1922 specimen.
Figure 14. Comparison between the theoretical and numerical results of SIF for TTC ASTM E1922 specimen.
Polymers 13 03129 g014
Figure 15. The geometry correction factor (Y) of MC ASTM E1922 and ASTM D3039 specimens with a/W = 0.5 and different debonding angles (θ).
Figure 15. The geometry correction factor (Y) of MC ASTM E1922 and ASTM D3039 specimens with a/W = 0.5 and different debonding angles (θ).
Polymers 13 03129 g015
Figure 16. Crack mouth opening displacement versus the applied load of TTC ASTM E1922 and ASTM D3039 specimens.
Figure 16. Crack mouth opening displacement versus the applied load of TTC ASTM E1922 and ASTM D3039 specimens.
Polymers 13 03129 g016
Table 1. The physical properties of the used materials.
Table 1. The physical properties of the used materials.
MaterialPhysical PropertiesValue
EpoxyDensity1.11 ± 0.02 Kg/L
Abrasion resistance1–6 cm3/50 cm2
Temperature resistancehumid 90 °C/dry 140 °C
E-glass FiberFilament diameter13 µm
Linear density2400 tex
Table 2. The mechanical properties of E-glass fiber and epoxy resin.
Table 2. The mechanical properties of E-glass fiber and epoxy resin.
Mechanical PropertiesE-Glass FiberEpoxy
Modulus of elasticity, E, (GPa) 823.24
In-plane shear modulus, G, (GPa)30.131.1912
Poisson’s ratio, ν,0.220.36
Table 3. The mechanical properties of E-glass fiber/epoxy composite.
Table 3. The mechanical properties of E-glass fiber/epoxy composite.
E11 (GPa)E22 (GPa)E33 (GPa)ν12ν13ν23G12 (GPa)G13 (GPa)G23 (GPa)
Glass fiber/epoxy32.96.76.70.290.290.492.462.462.25
Table 4. Hashin Model Constants.
Table 4. Hashin Model Constants.
X T   ( MPa ) X C   ( MPa ) Y T   ( MPa ) Y C   ( MPa ) S L   ( MPa ) S T   ( MPa )
Glass fiber/epoxy991.9502.7548.3536526.5
Table 5. The comparison between the experimental and numerical results of peak load.
Table 5. The comparison between the experimental and numerical results of peak load.
Smooth Specimen
Peak Load in (kN)
Through Thickness Crack Specimen
Peak Load in (kN)
Experimental12.867.11
Hashin Damage Model by ABAQUS12.786.93
Table 6. The experimental and numerical stress concentrator factors (KtExperimental) & (KtAnalytical), respectively, for TTC specimens of ASTM D3039.
Table 6. The experimental and numerical stress concentrator factors (KtExperimental) & (KtAnalytical), respectively, for TTC specimens of ASTM D3039.
1-Layer-Specimen2-Layers-Specimen
SmoothTTCSmoothTTC
Avg. Peak Load (kN)7.363.4713.126.87
Standard Deviation (kN)1.20.81.791.6
coefficient of variation, CV0.160.230.140.23
KtExperimental 2.12 1.91
KtAnalytical 8.75 8.75
Table 7. Comparison between the theoretical and numerical results of SIF and Y.
Table 7. Comparison between the theoretical and numerical results of SIF and Y.
a/wKI (MPa. √mm)—ASTM E1922 Specimen% Error
TheoreticalNumerical
0.2565.16563.010.38
0.3742.29737.290.67
0.4953.35924.523.02
0.51078.581046.352.98
0.61479.891470.310.64
Table 8. Comparison between the theoretical and numerical results of Y for TTC ASTM E1922 specimen.
Table 8. Comparison between the theoretical and numerical results of Y for TTC ASTM E1922 specimen.
a/wGeomertry Correction Factor (Y) —ASTM E1922 Specimen
TheoreticalNumerical
0.23.2532354033.094618148
0.33.669133513.433058442
0.44.3829190674.726436044
0.55.5337848165.348836997
0.67.4782893467.640451598
Table 9. The stress intensity factor (K1) and the geometry correction factor (Y) for TTC specimens as (a/W = 0.5).
Table 9. The stress intensity factor (K1) and the geometry correction factor (Y) for TTC specimens as (a/W = 0.5).
KI (MPa·√mm)—ASTM D3039 SpecimenY
Theoretical1701.352.83
Numerical1715.3122.85
Table 10. The apparent fracture toughness (KQ) and the geometry correction factor (Y) for TTC specimens with (a/W = 0.5).
Table 10. The apparent fracture toughness (KQ) and the geometry correction factor (Y) for TTC specimens with (a/W = 0.5).
a/W = 0.5ASTM E1922—TTC SpecimenASTM D3039—TTC Specimen
3-Layers1-Layer2-Layers
Load at crack initiation (kN)2.22.1013.32
Y5.532.832.83
KQ (MPa·√m)35.2350.751.13
a/W = 0.5ASTM D3039—MC SpecimenASTM D3039—MC Specimen
1-Layer2-Layer
θ = 80°θ = 80°
Load at crack initiation (kN)4.188.7
Y1.761.76
KQ (MPa·√m)62.7283.33
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

El-Sagheer, I.; Abd-Elhady, A.A.; Sallam, H.E.-D.M.; Naga, S.A.R. An Assessment of ASTM E1922 for Measuring the Translaminar Fracture Toughness of Laminated Polymer Matrix Composite Materials. Polymers 2021, 13, 3129. https://doi.org/10.3390/polym13183129

AMA Style

El-Sagheer I, Abd-Elhady AA, Sallam HE-DM, Naga SAR. An Assessment of ASTM E1922 for Measuring the Translaminar Fracture Toughness of Laminated Polymer Matrix Composite Materials. Polymers. 2021; 13(18):3129. https://doi.org/10.3390/polym13183129

Chicago/Turabian Style

El-Sagheer, Islam, Amr A. Abd-Elhady, Hossam El-Din M. Sallam, and Soheir A. R. Naga. 2021. "An Assessment of ASTM E1922 for Measuring the Translaminar Fracture Toughness of Laminated Polymer Matrix Composite Materials" Polymers 13, no. 18: 3129. https://doi.org/10.3390/polym13183129

APA Style

El-Sagheer, I., Abd-Elhady, A. A., Sallam, H. E. -D. M., & Naga, S. A. R. (2021). An Assessment of ASTM E1922 for Measuring the Translaminar Fracture Toughness of Laminated Polymer Matrix Composite Materials. Polymers, 13(18), 3129. https://doi.org/10.3390/polym13183129

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop