Next Article in Journal
Interface Behavior and Interface Tensile Strength of a Hardened Concrete Mixture with a Coarse Aggregate for Additive Manufacturing
Next Article in Special Issue
Numerical and Experimental Extraction of Dynamic Parameters for Pyramidal Truss Core Sandwich Beams with Laminated Face Sheets
Previous Article in Journal
Recent Progress in Hybrid Biocomposites: Mechanical Properties, Water Absorption, and Flame Retardancy
Previous Article in Special Issue
The Effect of Glue Cohesive Stiffness on the Elastic Performance of Bent Wood–CFRP Beams
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Experimental Investigation of Delamination in Composite Continuous Fiber-Reinforced Plastic Laminates with Elastic Couplings

by
Jakub Rzeczkowski
1,*,
Sylwester Samborski
2 and
Marcelo de Moura
3
1
Department of Fundamentals of Technology, Lublin University of Technology, Nadbystrzycka 38, 20-618 Lublin, Poland
2
Department of Mechanical Engineering, Lublin University of Technology, Nadbystrzycka 36, 20-618 Lublin, Poland
3
Department of Mechanical Engineering and Industrial Management, University of Porto, 4200-465 Porto, Portugal
*
Author to whom correspondence should be addressed.
Materials 2020, 13(22), 5146; https://doi.org/10.3390/ma13225146
Submission received: 21 October 2020 / Revised: 10 November 2020 / Accepted: 12 November 2020 / Published: 15 November 2020
(This article belongs to the Special Issue Advances in Fiber-Reinforced Polymer Composites)

Abstract

:
This paper presents an experimental evaluation of influence of the elastic couplings on the fracture toughness as well as on delamination initiation and propagation in carbon/epoxy composite laminates. For this purpose the mode I double cantilever beam (DCB) tests according to the American Society for Testing and Materials (ASTM) D5528 Standard were performed on specimens with different delamination interfaces and specific lay-ups composition exhibiting the bending-twisting (BT) and the bending-extension (BE) couplings. The critical strain energy release rates (mode I c-SERR, GIC) were calculated by using the classical methods, namely: the modified beam theory (MBT), the compliance calibration (CCM) and the modified compliance calibration (MCC). In order to evaluate an accuracy of the different methods, the values of c-SERR obtained by using standardized data reduction schemes were compared with values calculated by using the compliance based beam method (CBBM). All the methods give rise to comparable values of the GIC, which makes the CBBM an appealing choice, since it does not depend on crack length monitoring during the test. Initiation and propagation of delamination were investigated by using the acoustic emission (AE) technique. Moreover, the scanning electron microscope (SEM) analysis were performed after the experimental tests in order to investigate a fracture surface at delamination plane.

1. Introduction

Carbon fiber reinforced plastic (CFRP) composites made from long and continuous fibers demonstrate excellent properties such as low weight and also lifetime maintenance costs owing to their fatigue and corrosion resistance. Those features make composite laminates widely used in various industry sectors. On the other hand, CFRP laminates are often subjected to different loading conditions, which may cause damage of those materials. The most dangerous form of failure in composites is delamination. The fracture toughness represented by the strain energy release rate (SERR) can be determined in accordance with the American Society for Testing and Materials (ASTM) Standards, namely: for the opening mode I, the double cantilever beam (DCB) [1] test is commonly used, whereas the shearing fracture toughness can be obtained by employing the end-notched flexure (ENF) test [2]. The scope of the DCB test is measurement of the mode I fracture toughness (GIC) of unidirectional [0°]n laminates according to available standards. However, a multidirectional (MD) laminates with different fibers orientation at the delamination plane are commonly used in most composite structures applications. Many experiments were prepared on the mode I fracture of multidirectional laminates with different delamination interfaces [3,4]. Apart from experimental investigation several authors conducted numerical analysis on delamination problem in the continuous fiber-reinforced laminates [5,6,7,8,9,10,11]. When delamination propagates between different oriented layers, various undesirable phenomena such as mode-mixity and crack jumping to neighbor interface may occur [12]. Other parameters that may influence delamination resistance are: specimen geometry, stacking sequence and intralaminar damage [13]. Although the initiation toughness is often considered the most important in design point of view, there is significant interest in the measurements of propagation toughness values which are represented by the R-curve effects [14]. Schön et al. [15] prepared an experimental and numerical investigations of fracture toughness prepared on DCB specimens with different delamination interfaces as well as for different materials. Experimental results revealed, that non-zero interfaces had a plateau value approximately four times higher than the initial value, which was almost the same as for the zero interface. Higher toughness values could be probably caused by initial crack deviation from the original symmetric plane. In reality, anticlastic curvature in DCB specimens caused the highest values of the mode I SERR at the center and the lowest at the specimen edges. This non-uniformity was correlated with a non-dimensional ratios Dc and Bt introduced by Davidson et al [16]. Hence, in order to correlate the difference in strain energy release rate at a specimen’s center versus that at its edges a ratio Dc = D12/(D11D12) was used, whereas a stiffness ratio Bt = D16/D11, was found to correlate with the asymmetry in the SERR distribution about the center of the specimen’s width. In order to neglect the non-uniformity of GI in determination of GIC it was proposed that the Dc parameter should be less than 0.25. Nevertheless, Olsson [17] showed that the specimen was in a state of plane strain for small values of a0/b and globally was in a state of a plane stress for large aspect ratios. It was also revealed, that when specimens transferred from plane stress to plane strain the average value of the mode I fracture toughness decreased. Hence, at parameters of DCB specimen such as thickness (h), initial delamination length (a0) or width (b) significantly affected the GIC value. Therefore, Shokrieh [18] introduced a non-uniform ratio β = (GImaxGIavg)/GIavg%, which take into account the influence of abovementioned factors on the SERR distribution along width of the DCB specimen. In general β was a function of non-dimensional ratio Dc and geometrical ratios a0/b and a0/h, β = f(Dc, a0/b, a0/h). Results showed that the mode I fracture toughness of MD double cantilever beam specimens with 0°//0° fiber angle at delamination plane and β < 20% could be predictably measuring the GIC of the unidirectional UD plies with an error lest than 10%. Moreover, it was proved that β ratio had significant influence on the initiation delamination resistance and it was necessary to be considered in calculation of GIC by analytical relations which based on the beam theories. In practice, determination of the fracture toughness in composites may be difficult when elastic couplings occur in laminate [19]. For example, this is the case for laminates where the bending-twisting or the bending-extension couplings take place. Previous research prepared by the author of this paper [20] conducted on the DCB specimens exhibited elastic couplings under the mode I DCB test revealed, that non-zero components in the bending stiffness matrix [B] caused couplings between the bending moment and the twisting curvature as well as between the twisting moment and the mid-plane normal strains. Presence of shearing deformation and twisting curvature caused non-uniform delamination front which in results could produce an inaccuracies in the fracture toughness determination. In order to investigate the influence of irregular crack front on delamination initiation and the mode I strain energy release rate distributions Jiang et al. [21] prepared a modified DCB specimen. They concluded, that when the effects of visual deviation λ and the non-uniformity coefficient β were not taken into account the fracture toughness values obtained by using classical methods based on the beam theory were underestimated. Moreover, when the effect of curved delamination front during experimental measurement was neglected the errors in the mode I SERR values obtained for multidirectional DCB specimens could reached level up to 47.6% [22]. Another problem that can generate errors in measurements of the initiation and propagation values of SERR is visual observation of crack length through the test. Actually, the fracture process zone (FPZ) which develops ahead of the crack tip in result of micro-cracking and plasticity can cause unreliable monitoring of crack position during the mode I tests [23]. Abovementioned factors are not included in classical calculation schemes of the c-SERR recommended by the ASTM Standard. Moreover, visual detection of crack length are incredibly difficult to perform experimentally. To overcome this difficulties, an author of this paper proposed the compliance-based beam method (CBBM) [24]. This methodology is based on the crack equivalent concept and depends only on the specimens compliance during the test. In addition, this method do not require measurements of actual crack length during propagation which can produces non-negligible errors in the fracture toughness measurements [25,26]. Moreover, this methodology can be a great tool to create an artificial R-curve, which was proved in [27] where an excellent agreement between the experimental and numerical R-curves was achieved demonstrating adequacy of the proposed method.
Determination of crack propagation in structures plays an important role in the fracture mechanics of composites and can be investigated by using the acoustic emission (AE) technique [28]. The non-destructive AE method can be utilized in detecting the real delamination onset in damage initiation analysis. Nikbakht et al. [29] focused on damage analysis of the glass/epoxy multidirectional double cantilever beam specimens. Obtained AE data revealed that delamination initiated before any visual method can even see micro cracks. Saeedifar et al. [30] also used AE method during the mode I DCB tests prepared on glass/epoxy composite specimens to determine of crack tip position and to evaluate the interlaminar delamination resistance. They concluded, that obtained outcomes were similar to results calculated by using the nonlinearity (NL) method recommended in ASTM Standard. Generally, the acoustic emission methodology is a powerful technique to examine the initiation and propagation of delamination in composite laminates.
The main goal of this paper is to evaluate the influence of elastic couplings on behavior of the carbon/epoxy composite laminates subjected to the DCB tests. Specimens with different delamination interfaces and stacking sequences exhibited the bending-twisting and the bending-extension couplings were examined. Values of the critical strain energy release rate obtained by using standardized methods were compared with results obtained by using the CBBM method. In addition, experimental investigations of initiation and propagation of delamination were executed by using the acoustic emission technique and the scanning electron microscopy.

2. Elastic Coupling Phenomenon in Composites

The coupling phenomenon is an intrinsic attribute of the layered CFRP structures with different orientations of reinforcing fibers in subsequent plies when the lay-ups are non-symmetric. The constitutive equations describing all possible couplings in a laminate can be written in a block-matrix form [31]:
{ N M } = [ A B B D ] { ε 0 κ 0 }
where the [N] and the [M] matrices stand for the laminate’s forces and moments, respectively. The last column matrix depicts deformation of the neutral plane of the composite (ε0—strains, κ0—curvatures). The key element in Equation (1) is the block matrix in the square brackets; its constituents are the extensional stiffness matrix [A], the coupling stiffness matrix [B] and the bending stiffness matrix [D]. The influence of elastic coupling phenomena on behavior of CFRP laminates subjected to the mode I DCB test were analytically analyzed in previous paper [20].

3. Experimental Procedures

3.1. DCB Specimens

Experimental DCB specimens were performed on uniform thickness multidirectional (MD) carbon/epoxy composite laminates with different delamination interfaces: [0°19//0°19], [0°18/θ//-θ/0°18], [0°18/θ//θ/0°18] and [0°18/0//θ/0°18]. The fiber orientation angle θ was {30°, 45°, 60°, 90°}. Additionally, the bending-twisting (BT) and the bending-extension (BE) coupled laminates with specific stacking sequences which were presented in Table 1 were examined. In order to initiate delamination a 12 µm thick teflon PTFE foil was placed at mid-thickness of the laminate. In addition, both edges of specimens were coated white paint to increase accuracy of visual detection of crack tip. Specimen with the geometry: width B =25 mm, total thickness 2h = 4.98 mm and total length L =175 mm was presented in Figure 1. The loading piano hinges according to the Standard ASTM D5528 were adhesively bonded to the specimen allowing to apply crack opening force. Material properties of tested laminates were obtained during tensile tests and were summarized in Table 2.

3.2. Test Procedure

The double cantilever beam (DCB) tests were performed on Shimadzu ASG-X (Lublin, Poland) tensile testing machine (equipped with 1kN load cell) in accordance with the ASTM D5528 Standard procedure with a constant crosshead speed equal 1 mm/min. During all tests the load-displacement data were recorded by the Trapezium-X software. In addition, delamination onset as well as all propagation of delamination values were visually observed and marked on both specimen edges which were coated with a white paint. In order to ensure high precision of detecting the very beginning initiation of delamination - and in consequence to obtain high accuracy of the initial GIC determination, all experiments were supported by using acoustic emission (AE) technique. The piezoelectric (Fujicera 1045S) AE sensor was adhesively bonded at the top of laminates in 25 mm distance from specimen’s end edge. Elastic waves of the AE signal were recorded continuously by using AMSY-5 Acoustic Emission acquisition system at 10 MHz sampling rate. During all experiments the energy of elastic waves as well as the number of events were registered. The first increase of cumulative counts of the acoustic emission signal was taken into account as initiation of delamination and was called the AE point; note that this point was taken into account in calculations of the mode I critical strain energy release rates. In addition, the real delamination surfaces of the DCB specimens were investigated by using the scanning electron microscope (SEM) (Quanta FEG 450, Pisa, Italy).

3.3. Data Reduction Schemes

3.3.1. Classical Methods

Three classical data reduction schemes of calculation of the mode I critical strain energy release rate (c-SERR) is being proposed by the ASTM D5528 Standard, namely: the modified beam theory (MBT), the compliance calibration method (CCM) and the modified compliance calibration (MCC).
The MBT method uses an experimentally determined correction parameter Δ which is generating as a least squares plot of the cube root of compliance versus delamination length a. Subsequent parameters are: B specimen width, P load and δ displacement. The mode I interlaminar fracture toughness can be determined with the following equation:
G IC = 3 P δ 2 B ( a + | Δ | )
The CCM method uses an additional experimentally obtained parameter n which is created as the slope of the ln(C) versus ln(a). Then, the mode I critical strain energy release rate is calculated as follows:
G IC = n P δ 2 a B
In the MCC method the critical SERR can be calculated by using the following equation in which the A1 is the slope of a least squares plot of delamination length normalized by the specimen thickness a/h as a function of the cubic root of compliance:
G IC = 3 P 2 C 2 / 3 2 A 1 B h

3.3.2. Compliance-Based Beam Method (CBBM)

The critical strain energy release rates calculated by using previous methods (CCM, MBT, MCC) are depended on precise visual crack length measurements during the tests which is not easy to perform experimentally. Moreover, a fracture process zone which develops ahead of crack tip in a results of plastification or micro-cracking makes the visual crack tip localization more difficult. The compliance-based beam method (CBBM) which based on the crack equivalent concept takes into account these features and is described below. This methodology depends only on the specimen’s compliance (C = δ/P) which is expressed as follows:
C = 8 a 3 E 1 B h 3
where a specimen crack length, E1 longitudinal elastic modulus, B specimen width and h is arm’s height. Root rotation effect at crack tip is also taken into account by introduction a correction Δ calculated based on formula proposed by Hashemi et al. [32]:
Δ = h E 1 11 G 12 [ 3 2 ( Γ 1 + Γ ) 2 ] ; Γ = 1.18 E 1 E 2 G 12
Here E1, E2, G12 are material properties of carbon/epoxy laminate. Additionally, in order to account for scatter of elastic properties between different specimens a corrected flexural modulus can be calculated by using Equation (7) where C0 and a0 stands for initial compliance and crack length, respectively.
E f = 8 ( a 0 + | Δ | ) 3 B h 3 C 0
The equivalent crack length ae (Equation (8)) is a function of the current compliance C, which is calculated directly from the measured load-displacement (P-δ) curve and by using the corrected flexural modulus Ef
a e = ( C E f B h 3 8 ) 1 / 3
The values of the mode I critical strain energy release rate can be calculated from the Irwin-Kies equation:
G I = P 2 2 B d C d a
which leads to
G IC = 12 P 2 a e 2 B 2 h 3 E f

4. Effect of Elastic Couplings on Composite Laminates

Coupling behavior is dependent on the form of each of the elements in the extensional [A], the coupling [B] and the bending [D] stiffness matrices. In order to validate the influence of elastic couplings on the DCB laminates, values of the Aij, Bij and Dij components of the stiffness matrices and the non-dimensional parameters Dc and Bt [9] were calculated for chosen specimens according to the classical lamination theory (CLT) by using MatLab software environment. Obtained values were collected in Table 3 and Table 4. Experimental investigation of influence of the elastic couplings on delamination onset and propagation was presented in the subsequent section.
Composite laminates exhibit coupling between the in-plane shear strains and the extension when A16 and A26≠0, and between the bending moments and the twisting curvatures when D16 and D26≠0. For the 0°//0° laminate those conditions were equal zero, on contrary to the angle-ply laminates where the extensional stiffness matrix was fully coupled which initiated coupling between the normal stress and the shear strains (in-plane normal forces produced mid-plane shear strains). Considering the coupling stiffness matrix, the B11, B12, B22 and B66 components were also different from zero which induced shearing deformations under the bending load. Nevertheless, all terms in the coupling matrices [B] were incredibly small. On the other hand, fully coupled bending stiffness matrix [D] with non-zero D16 component responsible for twisting curvature could cause that specimen subjected to the mode I DCB test tried to twist itself around it’s central axis. In addition, for all tested specimens with different delamination interfaces, the non-dimensional Dc parameter was close to zero and the bending-twisting coupling intensity parameter Bt was equal zero. It means that influence of elastic couplings on the strain energy release rate was insignificant in this case. For the laminates with specific stacking sequences where the bending-twisting and the bending-extension couplings occurred, the values of the non-dimensional parameters were much greater and reached: Dc = 0.2384 and Bt = 0.2363 for the BT specimen as well as Dc = 0.4089 and Bt = 0 for the BE specimen. In this case the elastic coupling phenomena strongly influenced on behavior of the laminates and in consequence the SERR distribution along specimen width which was proved by an author of this paper in [33,34].
Values of the Dc and Bt parameters depended not only on the specific laminate configuration but also on the fiber θ angles which was presented in Figure 2. For fiber angles less than θ = 20° and greater than θ = 60° values of the non-dimensional parameter Dc were small. The greatest values were obtained for fiber angles range between θ = 35° and θ = 45° reaching maximum for the bending-twisting laminate equal Dc = 0.2623 for θ = 38° and for the bending-extension laminate Dc = 0.4089 for θ = 45°. For the BT specimen values of the non-dimensional parameter Bt were similar to Dc whereas for the BE laminate those values were equal zero for the entire range of fiber angles.

5. Experimental Results

5.1. Critical Strain Energy Release Rate (c-SERR)

Table 5 presents the results of calculation of the mode I critical strain energy release rate obtained for specimens with different delamination interfaces as well as for the laminates exhibiting elastic couplings. Values of the critical strain energy release rates calculated by using the CBBM method were discussed below. For the angle-ply laminates maximum value of the GIC was obtained for specimen with interface 60°//60° and was equal to 0.79 N/mm. Slightly less value of the critical strain energy release rate was obtained for the 0°//60° interface (c-SERR = 0.74 N/mm). Specimens with delamination interfaces 30°//−30° and 90°//90° exhibited values of the GIC on average level equal 0.58 N/mm. Considering 0°//45°, 0°//90° and 60°//−60° interfaces, the mode I strain energy release rates were equal to GIC = 0.39 N/mm, GIC = 0.45 N/mm and GIC = 0.45 N/mm, respectively. The lowest value of the mode I c-SERR was obtained for the unidirectional laminate and was equal to 0.15 N/mm. For the interface 45°//45° the critical strain energy release rate was greater and equal 0.18 N/mm. For specimens with delamination interface 0°//θ and θ//θ (without interface 0°//90°) with the growth of θ angle the values of mode I c-SERR also increased. Comparing the coupled composite laminates the greatest value of the mode I c-SERR was reached for specimen with the bending-extension coupling which was equal 0.65 N/mm. For specimen with the bending-twisting coupling the mode I fracture toughness was about two times smaller and equal 0.38 N/mm. Values of the critical strain energy release rate obtained by using methods recommended by the ASTM Standard were on similar level. Moreover, differences between the standardized data reduction schemes and the compliance based beam method were small. However, for all tested specimens values of the GIC obtained by using the CBBM were slightly greater. Nevertheless, those discrepancies between the standardized and non-standardized methods were not significant. In addition, transformation from the plane strain state to the plane stress state with increasing values of the a0/b ratio described in [10] had no major effect on the c-SERR values.
Typical load versus displacement curves which were obtained during experimental tests were presented in Figure 3 and Figure 4.

5.2. Analysis of Delamination Initiation and Propagation

The acoustic emission behavior during propagation of delamination process could be divided into three regions, namely: first region, where any damage mechanism was not occurred in laminate which in consequence no significant AE activity could be observed, second region where delamination initiation caused abrupt acoustic emission activity and third region with different AE activities caused by propagation of delamination. In Figure 5, Figure 6 and Figure 7 relations of the load versus the cumulative AE energy and counts for the laminates with 30°//30°, 45°//45°, 60°//60° delamination interfaces as well as for the BT and the BE specimens were presented.
Considering specific period of time, it could be observed, that for all tested specimens sudden changes in load led in a similar increase in the acoustic emission rate. This could be caused of by an increase in the rate of AE events occurred in the corresponding range of time. Moreover, the rate of changes and the trend in cumulative AE counts were similar to those of cumulative acoustic energy in every particular range of time. It proved that more AE energy emission caused more events. Although some similarities in general behavior of the counts and acoustic emission energy could be observed, however, there were some detailed differences for each laminates, which could be caused by the effect of delamination interface as well as by the specific sequence of laminate layers.
In the case of the bending-twisting coupled laminate (Figure 5a) the AE graphs experienced a number of sharp jump. Similar trend of the load and AE curve during propagation could be observed also for the laminate with 60°//60° interface (Figure 5b). This might be due to fiber bridging phenomena. The differently oriented fibers at the interface might bridge and then their abrupt debonding or break could produce sudden occurrence of large amount of AE emission and many events which resulted in violent crack propagation observed during the DCB experiment. In some cases, single bridging fibers were merged into bridged bundle which were obstacle for micro-cracks propagation. The abrupt released of that bundle resulted in acoustic energy jump caused by tremendous AE events at crack front occurred in that time. This phenomenon highlighted in Figure 5a, b. Microscopic SEM investigation of delamination surface after the DCB test revealed, that both for the BT and the 60°//60° specimen such bridged bundle occurred which was presented in Figure 8. The nature of fiber bridging depended on the fiber angle at the interface. For the 60°//60° interface bridged fibers released with greater AE energy than specimens with 30°//30° and 45°//45° interfaces.
For the BE specimen (Figure 6a) the load and the AE energy jumped very frequently, however, each jump emitted less energy relative to the bending-twisting laminate. It could be caused that releasing of acoustic emission energy and the cycles of bundle bridging occurred in a shorter period of time and consequently less events happened and less AE energy was emitted. Similar situation could be observed for the angle-ply laminates with delamination interfaces 30°//30° and 45°//45°, which was presented in Figure 6b and Figure 7a. In addition, the SEM analysis revealed that for the BE specimen a non-zero value of the B11 term caused that the stress field was strongly dominated by the mode II component. As the shear component increased the toughness also increased which was related to the formation of cusps (Figure 9). In macro-scale this occurred as an surface roughness increase. Microscopically cusps were formed by the development of angled cracks, that with increasing of the mode II component extended towards the plies and coalesced. Locally, cusps could cause a cleavage fracture which led to increase in absorption of fracture energy and in consequence increase in the fracture toughness. It was confirmed during determination of the c-SERR for the bending-extension specimen for which the non-zero B11 component of coupling stiffness matrix is responsible for the shearing effect under the bending load. Results of calculations revealed that the value of the mode I critical strain energy release rate for the BE laminate was greater than for the BT specimen for which the B11 component was close to zero. Moreover, microscopic SEM investigation of the BE specimen showed a number of cusps spread on delamination plane.
Relations between the load, cumulative count and energy of AE signal obtained for the specimens with θ//−θ and 0°//θ delamination interfaces were presented in Figure 10, Figure 11 and Figure 12. The greatest values of released AE energy during delamination process were observed for the laminates with interfaces 30°//−30°, 60°//−60° and 0°//30°. For the 45°//−45° specimen a characteristic energy jump caused by sudden propagation of delamination was observed. In the case of the 0°//45° and 0°//60° specimens, a rapid growth of the AE energy occurred. In addition, it was observed that for several tested specimens the AE energy reached established level at the end of propagation whereas the cumulative counts were continuously growing. It could be cause by a steady-state crack propagation connected with occurrence of the matrix cracking with constant rate, which is considered as a major source of acoustic energy. During propagation, such phenomena as fiber-matrix debonding or fiber breakage also generated the acoustic waves but they had less influence on total counts of the cumulative events.
Apart from investigation of delamination propagation, the AE technique could be successfully utilized as non-destructive method in detecting crack of initiation, which played extremely important role in determination of the fracture toughness. In order to obtain the critical strain energy release rate it was necessarily to identify the point on the load versus displacement plot where crack started to propagate. Acoustic emission is a natural phenomenon caused by propagation of transient elastic wave resulted from abrupt energy releasing inside the material. The point called AE was the first sharp increase in the cumulative of the acoustic energy counts. For example, for the BE specimen the first visually observed crack growth was registered when load reached 65 N. It could be caused by phenomenon of bridged bundles of fibers at the crack tip which resisted against a crack growth. However, stored strain energy in specimen (GI) achieved the critical value (GIC) when load reached 38 N. Hence, this point should be considered as delamination initiation. It was confirmed by increase in the cumulative counts of acoustic emission signal which corresponded to increase of the cumulative AE energy. Delamination onset (AE point) determined as first increase of acoustic signal was pointed out on all plots as the end of green line.
Proper recognizing of initiation point is often questionable, so the ASTM Standard recommended three different criteria: visually observed crack front (VIS), increase of compliance (5%) and deviation from linearity (NL). Nevertheless, initiation of delamination at interface with different fibers orientation is different from when the delamination initiates between unidirectional plies. Although, the acoustic emission signal for the specimens with 45°//45° and 60°//60° interfaces indicated early damage initiation a crack tip was not visible. In fact, the early micro-cracking in 45° and 60° interfaces created more extensive fracture process zone ahead of the crack tip because of small crack opening displacement. On the other hand, for the specimen with interface 30°//30° as well as for the BT and BE laminate because of greater crack opening displacement caused by existence of larger local normal strains the crack tip was visible. It was correlated with sharp jump on the load curve obtained during the experiment. In this case, evolution of damage along the fibers created a smaller FPZ and much more visible crack onset. Also, the component D16 might influence on crack tip visibility. Non-zero D16 term caused that specimen were tensioned on one side and compressed on another which resulted that crack tip was visible only on one edge of tested specimen.

6. Conclusions

An experimental study was conducted on the mode I double cantilever beam carbon/epoxy composite beams exhibited different delamination interfaces and elastic couplings. The mode I critical strain energy release rates were calculated by using the standardized data reduction schemes recommended by the ASTM D5528 Standard and the compliance-based beam method which based on the crack-equivalent concept. Results showed that the mode I fracture toughness depended on fiber angle at the interface and type of elastic coupling. Moreover, discrepancies between the GIC values obtained by using four different calculation methods were relatively small. Application of the CBBM method to specimens with different delamination interfaces and specific stacking sequences turned out to be advantageous relative to the classical ones since it did not required measurements of crack length during test and takes into account the energy dissipated at the fracture process zone.
Analytical analysis showed that the values of the non-dimensional parameter Dc and Bt depended not only on the lay-ups composition but also on fiber angles in specific stacking sequences. The greatest values of the Dc were obtained for fiber angles range between θ = 35° and θ = 45° for the bending-twisting as well as for the bending-extension specimens. For the BT laminate values of Bt parameter were similar to Dc whereas for the BE specimen, values of coupling intensity parameter were equal zero for entire angles set.
The AE technique was successfully utilized in determination of delamination initiation. The first increase of cumulative counts was considered as a delamination onset.
Acoustic emission investigations and SEM observations revealed that delamination initiation and propagation depend on type of elastic coupling. In the case of the BT and 60°//60° laminates the twisting curvatures caused that fibers merged into bundles and energy released was significantly large. On the other hand, for the 30°//30° and 45°//45° specimens released AE energy was smaller. This confirmed that bridging phenomena also depended on the interface type. Moreover, for the bending-extension specimen, the B11 term of the coupling stiffness matrix caused contribution of the mode II shearing components in the form of cusps which was proved during the microscopic observations. Presence of the shearing effects could increase total value of the strain energy release rate measured during the mode I DCB test. Moreover, elastic couplings influenced on laminates behavior and decreased visibility of crack tip at one side of specimen during experimental tests.

Author Contributions

Conceptualization, J.R.; Methodology, J.R., S.S. and M.d.M.; Experimental investigation, J.R.; Analytical and formal analysis, J.R.; Validation, J.R., S.S. and M.d.M.; Writing—original draft, J.R.; Review and editing, S.S. and M.d.M. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the Ministerial Research Project No. DEC-2016/21/B/ST8/03160, financed by the Polish National Science Centre. The third author acknowledges the “Laboratório Associado de Energia, Transportes e Aeronáutica” (LAETA) for the financial support by the project UID/EMS/50022/2013.

Conflicts of Interest

The authors declare no conflict of interest.

Nomenclature

acrack length
aeqequivalent crack length
a0initial crack length
hhalf thickness of specimen
ncorrection parameter
Aextensional stiffness matrix
Aijcomponent of the extensional stiffness matrix
A1correction parameter
Bcoupling stiffness matrix
Bspecimen width
Bijcomponent of the coupling stiffness matrix
Btnon-dimensional parameter
Ccompliance
C0initial compliance
Dbending stiffness matrix
Dijcomponent of the bending stiffness matrix
Dcnon-dimensional parameter
E1longitudinal Young’s modulus
E2transversal Young’s modulus
Efflexural modulus
G12shear modulus
GICcritical strain energy release rate
Lspecimen length
Mmoment vector
Nforce vector
Papplied load
ν12Poisson’s ratio
δdisplacement
Δcorrection parameter
Γ correction parameter
List of acronyms
AEacoustic emission
BEbending-extension coupled specimen
BTbending-twisting coupled specimen
CBBMCompliance-Based Beam Method
DCBDouble Cantilever Beam
MCCModified Compliance Calibration
CCMCompliance Calibration Method
MBTModified Beam Theory
SEMscanning electron microscope

References

  1. D5528 Standard Test. Method for Mode I Interlaminar Fracture Toughness of Unidirectional Fiber-Reinforced Polymer Matrix Composites; ASTM International: West Conshohocken, PA, USA, 2013.
  2. D7905 Standard Test. Method for Determination of the Mode II Interlaminar Fracture Toughness of Unidirectional Fiber-Reinforced Polymer Matrix Composites; ASTM International: West Conshohocken, PA, USA, 2018.
  3. De Gracia, J.; Boyano, A.; Arrese, A.; Mujika, F. Analysis of the DCB test of angle-ply laminates including residual stresses. Theor. Appl. Fract. Mech. 2018, 94, 197–204. [Google Scholar] [CrossRef]
  4. Pereira, A.B.; de Morais, A.B. Mode II interlaminar fracture of glass/epoxy multidirectional laminates. Compos. Part A 2004, 35, 265–272. [Google Scholar] [CrossRef]
  5. Alfano, G.; Crisfield, M.A. Finite element interface models for the delamination analysis of laminated composites: Mechanical and computational issues. Int. J. Numer. Meth. Eng. 2001, 50, 1701–1736. [Google Scholar] [CrossRef]
  6. de Borst, R. Numerical aspects of cohesive-zone models. Eng. Fract. Mech. 2003, 70, 1743–1757. [Google Scholar] [CrossRef]
  7. Gliszczynski, A.; Samborski, S.; Wiacek, N.; Rzeczkowski, J. Mode I Interlaminar Fracture of Glass/Epoxy Unidirectional Laminates. Part II: Numerical Analysis. Materials (Basel) 2019, 12, 1607. [Google Scholar] [CrossRef] [Green Version]
  8. Herráez, M.; González, C.; Lopes, C.S. A numerical framework to analyze fracture in composite materials: From R-curves to homogenized softening laws. Int. J. Solids Struct. 2018, 134, 216–228. [Google Scholar] [CrossRef]
  9. Orifici, A.C.; Krueger, R. Benchmark assessment of automated delamination propagation capabilities in finite element codes for static loading. Finite Elem. Anal. Des. 2012, 54, 28–36. [Google Scholar] [CrossRef] [Green Version]
  10. Rybicki, E.F.; Kanninen, M.F. A finite element calculation of stress intensity factors by a modified crack closure integral. Eng. Fract. Mech. 1977, 9, 931–938. [Google Scholar] [CrossRef]
  11. Turon, A.; Dávila, C.G.; Camanho, P.P.; Costa, J. An engineering solution for mesh size effects in the simulation of delamination using cohesive zone models. Eng. Fract. Mech. 2007, 74, 1665–1682. [Google Scholar] [CrossRef]
  12. Gong, Y.; Zhang, B.; Zhao, L.; Zhang, J.; Hu, N.; Zhang, C. R-curve behaviour of the mixed-mode I/II delamination in carbon/epoxy laminates with unidirectional and multidirectional interfaces. Compos. Struct. 2019, 223, 110949. [Google Scholar] [CrossRef]
  13. Gong, Y.; Zhang, B.; Mukhopadhyay, S.; Hallett, S.R. Experimental study on delamination migration in multidirectional laminates under mode II static and fatigue loading, with comparison to mode I. Compos. Struct. 2018, 261, 683–698. [Google Scholar] [CrossRef] [Green Version]
  14. Shokrieh, M.M.; Heidari-Rarani, M. Effect of stacking sequence on R-curve behavior of glass/epoxy DCB laminates with 0°//0° crack interface. Mater. Sci. Eng. A 2011, 529, 265–269. [Google Scholar] [CrossRef]
  15. Schön, J.; Nyman, T.; Blom, A.; Ansell, H. A numerical and experimental investigation of delamination behaviour in the DCB specimen. Compos. Sci. Technol. 2000, 60, 173–184. [Google Scholar] [CrossRef]
  16. Davidson, B.D. Three Dimensional Analysis and Resulting Design Recommendations for Unidirectional and Multidirectional End-Notched Flexure Tests. J. Compos. Mater. 1995, 29, 2018–2033. [Google Scholar] [CrossRef]
  17. Olsson, R.A. Simplified improved beam analysis of the DCB specimen. Compos. Sci. Technol. 1992, 43, 329–338. [Google Scholar] [CrossRef]
  18. Shokrieh, M.M.; Heidari-Rarani, M.; Rahimi, S. Influence of curved delamination front on toughness of multidirectional DCB specimens. Compos. Struct. 2012, 94, 1359–1365. [Google Scholar] [CrossRef]
  19. York, C.B.; de Almeida, S.F.M. On Extension-Shearing Bending-Twisting coupled laminates. Compos. Struct. 2017, 164, 10–22. [Google Scholar] [CrossRef] [Green Version]
  20. Rzeczkowski, J.; Samborski, S.; Valvo, P.S. Effect of stiffness matrices terms on delamination front shape in laminates with elastic couplings. Compos. Struct. 2020, 23, 111547. [Google Scholar] [CrossRef]
  21. Jiang, Z.; Wan, S.; Keller, T.; Fang, Z.; Vassilopoulos, A.P. Vassilopoulos. Influence of curved delamination front on R-curve of DCB specimen. Compos. Struct. 2019, 223, 111311. [Google Scholar] [CrossRef]
  22. Jiang, Z.; Wan, S.; Zhong, Z.; Li, S.; Shen, K. Effect of curved delamination front on mode-I fracture toughness of adhesively bonded joints. Eng. Fract. Mech. 2015, 138, 73–91. [Google Scholar] [CrossRef]
  23. Amrutharaj, G.S.; Lam, K.Y.; Cotterell, B. Fracture process zone concept and delamination of composite laminates. Theor. Appl. Fract. Mech. 1995, 24, 57–64. [Google Scholar] [CrossRef]
  24. de Moura, M.F.S.F.; de Morais, A.B. Equivalent crack based analyses of ENF and ELS tests. Eng. Fract. Mech. 2008, 75, 2584–2596. [Google Scholar] [CrossRef]
  25. Rzeczkowski, J. An experimental analysis of the end-notched flexure composite laminates beams with elastic couplings. Contin. Mech. 2020, 98, 210. [Google Scholar]
  26. Rzeczkowski, J.; Samborski, S.; Paśnik, J. Experimental Verification of the DCB Test Configuration Applicability to Mechanically Coupled Composite Laminates. IOP Conf. Ser. Mater. Sci. Eng. 2018, 416, 12055. [Google Scholar] [CrossRef]
  27. De Moura, M.F.S.F.; Campilho, R.D.S.G.; Gonçalves, J.P.M. Crack equivalent concept applied to the fracture characterization of bonded joints under pure mode I loading. Compos. Sci. Technol. 2008, 68, 2224–2230. [Google Scholar] [CrossRef] [Green Version]
  28. Lissek, F.; Haeger, A.; Knoblauch, V.; Hloch, S.; Pude, F.; Kaufeld, M. Acoustic emission for interlaminar toughness testing of CFRP: Evaluation of the crack growth due to burst analysis. Compos. Part B Eng. 2018, 136, 55–62. [Google Scholar] [CrossRef]
  29. Nikbakht, M.; Yousefi, J.; Hosseini-Toudeshky, H.; Minak, G. Delamination evaluation of composite laminates with different interface fiber orientations using acoustic emission features and micro visualization. Compos. Part. B Eng. 2017, 113, 185–196. [Google Scholar] [CrossRef]
  30. Saeedifar, M.; Fotouhi, M.; Najafabadi, M.A.; Toudeshky, H.H.; Minak, G. Prediction of quasi-static delamination onset and growth in laminated composites by acoustic emission. Compos. Part. B Eng. 2016, 85, 113–122. [Google Scholar] [CrossRef]
  31. Kaw, A.K. Mechanics of Composite Materials, 2nd ed.; CRC, Taylor & Francis: Boca Raton, FL, USA, 2006; ISBN 9780849313431. [Google Scholar]
  32. Hashemi, S.; Kinloch, A.J.; Williams, J.G. The analysis of interlaminar fracture in uniaxial fibre-polymer composites. Proc. R. Soc. A 1990, 427, 173–199. [Google Scholar]
  33. Samborski, S. Analysis of the end-notched flexure test configuration applicability for mechanically coupled fiber reinforced composite laminates. Compos. Struct. 2017, 163, 342–349. [Google Scholar] [CrossRef]
  34. Samborski, S. Numerical analysis of the DCB test configuration applicability to mechanically coupled Fiber Reinforced Laminated Composite beams. Compos. Struct. 2016, 152, 477–487. [Google Scholar] [CrossRef]
Figure 1. Geometry of the DCB specimen.
Figure 1. Geometry of the DCB specimen.
Materials 13 05146 g001
Figure 2. Non-dimensional parameters: (a) values of the DC parameter calculated for different fiber angles θ for the bending-twisting and the bending-extension laminates, (b) values of the DC and Bt parameters calculated for different fiber angles θ for the bending-twisting laminate.
Figure 2. Non-dimensional parameters: (a) values of the DC parameter calculated for different fiber angles θ for the bending-twisting and the bending-extension laminates, (b) values of the DC and Bt parameters calculated for different fiber angles θ for the bending-twisting laminate.
Materials 13 05146 g002
Figure 3. Load-displacement curves obtained during the DCB experiments: (a) for the UD and the 0°//θ laminates, (b) for the UD and the θ//θ laminates.
Figure 3. Load-displacement curves obtained during the DCB experiments: (a) for the UD and the 0°//θ laminates, (b) for the UD and the θ//θ laminates.
Materials 13 05146 g003
Figure 4. Load-displacement curve obtained during the DCB experiments: (a) for the UD and the θ//−θ laminates, (b) for the UD and the coupled laminates.
Figure 4. Load-displacement curve obtained during the DCB experiments: (a) for the UD and the θ//−θ laminates, (b) for the UD and the coupled laminates.
Materials 13 05146 g004
Figure 5. Relation between load, cumulative count and energy of AE signal obtained during the DCB experiment: (a) for specimen exhibiting the bending-twisting elastic coupling, (b) for specimen with 60°//60° delamination interface.
Figure 5. Relation between load, cumulative count and energy of AE signal obtained during the DCB experiment: (a) for specimen exhibiting the bending-twisting elastic coupling, (b) for specimen with 60°//60° delamination interface.
Materials 13 05146 g005
Figure 6. Relation between load, cumulative count and energy of AE signal obtained during the DCB experiment: (a) for specimen exhibiting the bending-extension elastic coupling, (b) for specimen with 30°//30° delamination interface.
Figure 6. Relation between load, cumulative count and energy of AE signal obtained during the DCB experiment: (a) for specimen exhibiting the bending-extension elastic coupling, (b) for specimen with 30°//30° delamination interface.
Materials 13 05146 g006
Figure 7. Relation between load, cumulative count and energy of AE signal obtained during the DCB experiment: (a) for specimen with 45°//45° delamination interface, (b) for UD laminate.
Figure 7. Relation between load, cumulative count and energy of AE signal obtained during the DCB experiment: (a) for specimen with 45°//45° delamination interface, (b) for UD laminate.
Materials 13 05146 g007
Figure 8. Delamination surfaces after the scanning electron microscope (SEM) investigations for the specimen with 60°//60° delamination interface and for the bending-twisting (BT) laminate.
Figure 8. Delamination surfaces after the scanning electron microscope (SEM) investigations for the specimen with 60°//60° delamination interface and for the bending-twisting (BT) laminate.
Materials 13 05146 g008
Figure 9. Delamination surfaces after the SEM investigations of the bending-extension specimen.
Figure 9. Delamination surfaces after the SEM investigations of the bending-extension specimen.
Materials 13 05146 g009
Figure 10. Relation between load, cumulative count and energy of AE signal obtained during the DCB experiment: (a) for specimen with 30°//−30° delamination interface, (b) for specimen with 45°//−45° delamination interface.
Figure 10. Relation between load, cumulative count and energy of AE signal obtained during the DCB experiment: (a) for specimen with 30°//−30° delamination interface, (b) for specimen with 45°//−45° delamination interface.
Materials 13 05146 g010
Figure 11. Relation between load, cumulative count and energy of AE signal obtained during the DCB experiment: (a) for specimen with 60°//−60° delamination interface, (b) for specimen with 0°//30° delamination interface.
Figure 11. Relation between load, cumulative count and energy of AE signal obtained during the DCB experiment: (a) for specimen with 60°//−60° delamination interface, (b) for specimen with 0°//30° delamination interface.
Materials 13 05146 g011
Figure 12. Relation between load, cumulative count and energy of AE signal obtained during the DCB experiment: (a) for specimen with 0°//45° delamination interface, (b) for specimen with 0°//60° delamination interface.
Figure 12. Relation between load, cumulative count and energy of AE signal obtained during the DCB experiment: (a) for specimen with 0°//45° delamination interface, (b) for specimen with 0°//60° delamination interface.
Materials 13 05146 g012
Table 1. Laminate configurations.
Table 1. Laminate configurations.
LaminateStacking Sequence (for One Sub-laminate)
BT(45°/0°/45°/45°/0°/−45°/0°/−45°/−45°/−45°/−45°/0°/−45°/45°/0°/0°/45°/45°)
BE(45°/−45°/0°/−45°/0°/45°/90°/45°/−45°)
Table 2. Mechanical properties of carbon/epoxy laminate.
Table 2. Mechanical properties of carbon/epoxy laminate.
E1 (GPa)E2 (GPa)ν12 (-)G12 (GPa)
112.1057.4210.2703.338
Table 3. Values of the stiffness matrices components calculated for the laminates with different delamination interfaces.
Table 3. Values of the stiffness matrices components calculated for the laminates with different delamination interfaces.
InterfaceA (MPa × mm) × 105B (MPa × mm2) × 10−10D (MPa × mm3) × 105DcBt
0°//0°5.01740.08970−0.29100.003408.29460.148300.00480
0.08970.332100.00340.004500.14830.54910
000.148700−0.0011000.2458
30°//30°4.89800.14010.0888−0.29100.003408.29390.14850.00050.00480
0.14010.35060.03050.00340.004500.14850.54920.0002
0.08880.03050.199100−0.00110.00050.00020.2461
45°//45°4.81230.15700.0689−0.29100.003408.29340.14860.00040.00480
0.15700.40260.06890.00340.004500.14860.54950.0004
0.06890.06890.2160000.00110.00040.00040.2462
60°//60°4.76020.14010.0305−0.29100.003408.29310.14850.00020.00480
0.14010.48840.08880.00340.004500.14850.55000.0005
0.03050.08880.199100−0.00110.00020.00050.2461
Table 4. Values of the stiffness matrices components calculated for laminates with the bending-twisting (BT) and the bending-extension (BE) couplings.
Table 4. Values of the stiffness matrices components calculated for laminates with the bending-twisting (BT) and the bending-extension (BE) couplings.
LaminateA (MPa × mm) × 105B (MPa × mm2) × 10−11D (MPa × mm3) × 104DcBt
BT1.42570.451200.090900.09096.60612.09061.56080.23840.2363
0.45120.598900−0.090902.09062.77511.5608
000.48240.090900.18191.56081.56082.2353
LaminateA (MPa × mm) × 104B (MPa × mm2) × 103D (MPa × mm3) × 103DcBt
BE5.75062.25600−3.6104005.18033.306000.40890
2.25604.3726003.610403.30605.16060
002.4121000003.4870
Table 5. Critical strain energy release rate (c-SERR) obtained for different DCB specimens by using different calculation methods (average values and standard deviations).
Table 5. Critical strain energy release rate (c-SERR) obtained for different DCB specimens by using different calculation methods (average values and standard deviations).
InterfaceCCMMBTMCCCBBM
GIC (N/mm)Error (%)GIC (N/mm)Error (%)GIC (N/mm)Error (%)GIC (N/mm)Error (%)
0°//0°0.140.08 0.140.170.142.080.150.26
0°//30°0.182.500.192.610.191.950.216.28
0°//45°0.393.490.384.070.384.680.392.21
0°//60°0.534.330.511.160.501.930.7410.08
0°//90°0.4010.280.4213.400.4010.550.4512.29
30°//30°0.161.680.172.250.172.110.172.37
30°//−30°0.4911.540.511.920.5111.040.5616.66
45°//45°0.161.590.161.440.161.360.180.76
45°//−45°0.227.320.228.640.238.760.248.94
60°//60°0.673.100.664.110.635.260.794.64
60°//−60°0.464.210.464.870.423.590.453.75
90°//90°0.542.200.543.600.546.830.619.33
BE0.513.560.502.320.495.380.659.06
BT0.363.060.362.590.342.110.381.47
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Rzeczkowski, J.; Samborski, S.; de Moura, M. Experimental Investigation of Delamination in Composite Continuous Fiber-Reinforced Plastic Laminates with Elastic Couplings. Materials 2020, 13, 5146. https://doi.org/10.3390/ma13225146

AMA Style

Rzeczkowski J, Samborski S, de Moura M. Experimental Investigation of Delamination in Composite Continuous Fiber-Reinforced Plastic Laminates with Elastic Couplings. Materials. 2020; 13(22):5146. https://doi.org/10.3390/ma13225146

Chicago/Turabian Style

Rzeczkowski, Jakub, Sylwester Samborski, and Marcelo de Moura. 2020. "Experimental Investigation of Delamination in Composite Continuous Fiber-Reinforced Plastic Laminates with Elastic Couplings" Materials 13, no. 22: 5146. https://doi.org/10.3390/ma13225146

APA Style

Rzeczkowski, J., Samborski, S., & de Moura, M. (2020). Experimental Investigation of Delamination in Composite Continuous Fiber-Reinforced Plastic Laminates with Elastic Couplings. Materials, 13(22), 5146. https://doi.org/10.3390/ma13225146

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