Next Article in Journal
Mixtures of Lubricants and Ecological Refrigerants under Starved Lubrication Conditions
Previous Article in Journal
Response of Apple Flesh to Compression under the Quasi-Static and Impact Loading Conditions
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Experimental and Numerical Research of Delamination Process in CFRP Laminates with Bending-Twisting Elastic Couplings

by
Jakub Rzeczkowski
1,* and
Sylwester Samborski
2
1
Faculty of Technology Fundamentals, Lublin University of Technology, Nadbystrzycka 38, 20-618 Lublin, Poland
2
Faculty of Mechanical Engineering, Lublin University of Technology, Nadbystrzycka 36, 20-618 Lublin, Poland
*
Author to whom correspondence should be addressed.
Materials 2022, 15(21), 7745; https://doi.org/10.3390/ma15217745
Submission received: 15 October 2022 / Revised: 29 October 2022 / Accepted: 1 November 2022 / Published: 3 November 2022

Abstract

:
This paper aims at experimental and numerical research of delamination process in carbon/epoxy composite laminates with different fiber orientation angles in stacking sequence exhibiting the bending–twisting elastic couplings. Experimental specimens were subjected to the double cantilever beam (DCB) tests according to the ASTM D5528 regulations. Values of the mode I strain energy release rates were calculated by using three different data reduction schemes: the modified beam theory, the compliance calibration method and the modified compliance calibration. Determination of delamination initiation point was conducted in twofold way: by visual observation of crack tip using high resolution camera and by utilization of the acoustic emission technique. Numerical analyss were prepared in Abaqus/CAE Software environment by using the virtual crack closure technique (VCCT). The numerical beam model consisted of SC8R continuum shell elements. Obtained outcomes revealed that extensive fiber bridging phenomenon occurring during delamination process pronouncedly affected propagation values of strain energy release rate (GIprop) and numerically obtained load–displacement curves. Nevertheless, in initial stage of delamination, results obtained by using the VCCT were in agreement with experimental data. The greatest value of the mode I fracture toughness equal 0.56 N/mm was obtained for the BT45 laminate.

1. Introduction

Delamination is a common mode of failure in composite laminates. This phenomenon can be induced easily inside composite structures or near laminate edges where three-dimensional stress state occurs. In the current literature, there is considerable amount of papers devoted to research of delamination process, both numerically and experimentally [1,2,3,4,5,6,7,8,9]. According to the fracture mechanics approach, delamination process can be analyzed by comparing fracture toughness represented by the energy release rate (ERR) at delamination interface. When unidirectional (UD) laminates are tested, a pure mode of the ERR can be determined by using procedures described by the ASTM and ISO Standards [10,11]. On the other hand, when multidirectional (MD) laminates with specific stacking sequences are involved, a mix-mode condition can be present. Hence, due to dependency of mix-modes on interface fracture toughness, decomposition of energy release rate on the mode I, II and III components is necessary [12]. Moreover, it becomes extremely important during testing laminates exhibiting elastic couplings [13,14,15] in which in-plane shear deformations and twisting curvatures cause irregularities in strain energy releases distributions along specimens width [16], as well as curved delamination front [17,18,19,20]. Elastic coupling phenomena can be described by using the classical lamination theory (CLT). In accordance with CLT, resultant moments and forces with respect to laminate’s mid-plane can be expressed by the constitutive equation given below (Equation (1)):
[ N M ] = [ A B B D ] [ ε κ ]
where N is resultant force, M is resultant moment, A is in-plane stiffness matrix, B is coupling matrix, D is flexural stiffness matrix, ε is strains and κ is curvatures. Type of elastic coupling depends on layer configuration in composite stacking sequence. When laminate is symmetric, all components in the B matrix are equal to zero, hence the effects of elastic couplings vanish. On the other hand, when the B16 and B26 elements of coupling stiffness matrix are non-zero, then couplings between in-plane forces and twisting curvatures are induced in laminate. Moreover, non-zero components D16 and D26 introduce coupling between bending moments and twisting curvatures. These dependencies were extensively investigated by Davidson et al. [21]. They prepared finite element analysis of ERR distributions in anti-symmetric and symmetric laminates subjected to the double cantilever beam and the end-notched flexure configurations. Obtained results showed that asymmetry in ERR widthwise distributions as well as differences between minimum and maximum values of the energy release rate at crack front were the function of non-dimensional Dc and Bt parameters that were described in detail in Section 2. They also suggested that laminates should be characterized by Dc value less or equal to 0.25 to ensure uniform ERR distributions and crack front shape. In the available literature, there is a limited number of papers devoted to delamination problems in composite laminates with elastic couplings. In [22], the researchers presented numerical study of the double cantilever beam laminates exhibiting different types of elastic couplings. Results of simulations revealed that the ply stacking sequence in DCB model had a significant effect on the strain energy release rate distributions along specimens’ width. Similar results were drawn by Samborski [23]. He prepared numerical finite element analysis of the mode I SERR distributions at different boundary conditions (BCs). He concluded that initial BCs had essential influence on both the strain energy release rate widthwise distributions and the crack front shape. In [24], J. De Gracia et al. conducted analytical analysis and experimental tests of anti-symmetric and symmetric laminates. Regarding the former case, they obtained higher interlaminar fracture toughness values. It was attributed to non-uniform crack tip that was induced by the bending–twisting elastic coupling. As is stated by Sun et al. [25], curved crack front can lead to errors during experimental determination of delamination resistance. Values of the critical strain energy release rates (GC) obtained by using classical beam theory could be overestimated in this case. As underlined above, due to the fact that elastic couplings are extremely complex phenomena, their influence on delamination process should be more extensively investigated. Therefore, the current paper is devoted to experimental and numerical research of the mode I double cantilever beam laminates with elastic coupling. Novelty in the current paper is evaluation of influence of various fiber orientation angles in stacking sequence exhibiting the bending–twisting elastic couplings on the initial delamination resistance, as well as on the propagation values of the mode I strain energy release rate presented in the form of R-curves. Effect of fiber bridging phenomena was also investigated. Application of the acoustic emission technique in conjunction with visual recording of propagated crack tip by using high resolution camera in real time is also a new approach. This method was not used before and allows to detect delamination initiation point with high precision.

2. Tested Laminates

The DCB experimental specimens were carbon/epoxy composites with specific layup sequence consisting of 36 layers that generate the bending–twisting elastic coupling. The main dimensions of beam laminates were: total length L = 210 mm, width B = 25 mm and total thickness h = 5 mm, as presented in Figure 1. The pre-crack was defined by inserting a PTFE teflon foil (13 μm thick) at mid-plane between upper and lower sub-laminate. The length of film insert was 75 mm. To aid visual observation of crack propagation tip, both edges of laminates were coated with typewriter white fluid. In addition, a vertical solid line was marked in every 1 mm for the first 5 mm length and in every 5 mm for the remaining 30 mm.
Mechanical properties of CFRP laminates were obtained during tensile tests prepared in accordance with the ASTM D3039 [26] and the ASTM D3518 [27] Standards. Obtained results are collected in Table 1. In order to investigate the influence of fiber orientation angles on delamination process in coupled DCB laminates, three groups of specimens with stacking sequence [θ /0°/θ/θ/0°/−θ /0°/−θ /−θ/−θ/−θ/0°/−θ/θ/0°/0°/θ/θ] inducing the bending–twisting elastic couplings were chosen, for which the θ angles were {30°, 45°, 60°} (Table 2).
Effects of elastic coupling phenomena in tested laminates were assessed by calculations of non-dimensional parameters proposed by Davidson [21]. They were as follows: the bending stiffness ratio Dc and the bending–twisting coupling intensity Bt. Those parameters are related to non-uniformity in strain energy release rate SERR distributions along delamination front and can be determined as a ratio of D11, D12, D22, D16 elements in the bending stiffness matrix [D]. Equations that described non-dimensional parameters are presented below:
D c = D 12 2 D 11 D 22
B t = | D 16 | D 11

3. Experimental Procedures

Experimental tests were prepared on straight-sided DCB geometry specimens according to the ASTM D5528 regulations. Prior to experiments, piano hinges were bonded to laminates by using Loctite HY4070 epoxy glue. Subsequently, specimens were mounted in properly calibrated Shimadzu ASG-X tensile testing machine with sufficiently sensitive 1 kN load cell. All tests were conducted with constant crosshead velocity equal to 1 mm/min. The load–displacement responses were registered in real-time by plotting the P-δ charts. In order to monitor of crack tip and its extension during mode I fracture loading a Nikon D500 4K photo camera (Nikon, Tokyo, Japan) with set of objectives (Nikon Nikkor AF-S DX Micro 40 mm f/2.8G) and specialized led lamps for crack propagation monitoring was positioned on one side of the laminate. At first, crack tip was observed during pre-crack loading up to visually observed movement of cleavage from its initial position. Hence, the a0—initial delamination length value was determined. Afterward, a laminate was unloaded. During reloading procedure, front of growing crack was recorded in each 1 mm for the first 5 mm crack growth, and then, for every 5 mm up to 30 mm crack extension. In order to increase precision of defining the exact delamination initiation point, the acoustic emission (AE) technique was used [28]. The acoustic emission set Vallen AMSY-5 (Figure 2) supplemented with piezoelectric sensor mounted to tested laminate measured in real time various parameters of acoustic waves. The first growth of cumulative energy was treated as delamination onset and was defined as the AE point. Exemplary application of the acoustic emission technique combined with visual monitoring of crack to detection of delamination initiation point was depicted in Figure 3. In all cases, the AE points determined from cumulative energy plots corresponded to visually observed initial movement of propagated cleavage.

4. Analytical Determination of Energy Release Rates

In order to determine the mode I critical strain energy release rates several methods of calculation were used. Namely: the modified beam theory (MBT), the compliance calibration method (CCM) and the modified compliance calibration (MCC). The classical beam theory assumed perfectly built-in DCB beam rigidly fixed at delamination front. According to this theory the mode I fracture toughness can be expressed by following equation, where P is load, δ is load–displacement point, B is specimens width and a is current crack length.
G I = 3 P δ 2 B a
In reality, the branches of double cantilever beam laminate were not perfectly clamped due to presence of initial crack region. This could cause overestimation of the GI values. To establish unwanted deviation from classical beam theory, the correction parameter Δ, obtained experimentally, should be taken into account in data reduction scheme. Hence, the MBT method could be expressed as:
G I = 3 P δ 2 B ( a + Δ )
In case of the compliance calibration method, the n correction parameter should be determined as an exponent from the slope of log (δi/Pi) versus crack length ai. Then the CCM method is described by the equation:
G I = n P δ 2 B a
The third method of calculation of the mode I interlaminar fracture toughness take into account a slope of the plot creating from function of cube root of compliance versus ratio a/h in the form of the A1 parameter. Therefore, values of the strain energy release rate could be determined from the following formula:

5. Numerical Analysis

Numerical validations of experimental tests were prepared in the Abaqus/CAE Standard commercial finite element software. All analysis were conducted on Dell T3500 (Intel Xeon W3540 8-core processor, 2.93 GHz, 12 GB RAM, RAID-5 SATA disk matrix) workstation. During numerical simulations load versus displacement plots were obtained for bending–twisting laminates with different fiber orientation angles in particular layers. Additionally, contributions of different modes at crack front were assessed by determining strain energy release rates (SERR) widthwise distributions along specimen width. The finite element model had a form of double cantilever beam that consisted of two sub-laminates bonded together. The main dimensions of numerical DCB beam were: 210 mm total length, 25 mm width and 5 mm total thickness. In addition, whole model was divided into two branches at one beam edge on the length a0 = 45 mm which served as delamination initiator. The material properties obtained during experimental tensile tests were assigned to each layer separately. Specific stacking sequences were modeled in Abaqus/CAE composite layup module and were oriented in accordance with global coordinate system. Due to simulations were prepared on laminates exhibiting the bending–twisting elastic couplings, for which shearing deformations and twisting curvatures were expected, the FE model was built in 61050 linear hexahedral SC8R continuum shell elements in total. Unlike to conventional shell, the continuum element allow finite membrane deformations and large rotations. Moreover, the SC8R include transverse shear effects, thickness change, as well as allow for a richer transverse shear stress and force prediction. In addition, the continuum shell elements are more accurate in contact modeling than conventional shell, which is considered a big advantage. The meshed bending–twisting beam model was presented in Figure 4.
In order to improve accuracy of numerical results, number of elements was increased along crack propagation path. To modeled a delamination process different methods of simulations is available such as: the cohesive zone model (CZM) [29], the continuum damage model (CDM) [30] and the virtual crack closure technique (VCCT). First method based on an exponential traction low that characterizing the behavior of the elements at the interface. The second one (CDM) utilizes stiffness degradation of adhesive elements imposed by a damage parameter. Nevertheless, for the needs of current research the origirnal virtual crack closure (VCCT) criterion, that based on linear elastic fracture mechanics, was applied. This method uses Irwin principle that assumes the change in elastic strain energy ΔU to be equal to the work Wclosure required to close the crack. The strain energy release rate is calculated by using Clapeyron’s theory. In Figure 5 presented exemplary 2D finite element representation of the VCCT method.
The mode I (GI) and the mode II (GII) strain energy release components can be determined from Equations (7) and (8), where Sjx, Sjy are forces at node j, uiy, uix, uiy*, uix* are displacements at node i and Δa is length of crack element.
G I = 1 2 Δ a S j y ( u i y u i y * )
G II = 1 2 Δ a S j x ( u i x u i x * )
For pure mode I and mode II, crack growth propagation is predicted from single mode criterion. Nevertheless, as it was presented in previous research conducted by the authors [31,32,33], coupled laminates were dominated by contribution of different modes. Therefore, pure mode criterion should be replaced by mix-modes criterion that takes into account interactions between different modes. Hence, the Reeder-law (Equation (9)) was used in numerical anaylsis.
G IC + ( G IIC G IC + G IIIC G IIC ( G III G II + G III ) ) ( G II + G III G T ) η 1
Here, the GT = GI + GII + GIII is the total energy release rate and GIC, GIIC, GIIIC are critical strain energy release rates in opening mode I, in-plane shear mode II and transverse shear mode III respectively. Detailed material and fracture parameters implemented to numerical analysis were collected in Table 3.
Initial boundary conditions for the mode I DCB numerical model were presented in Figure 6.
For this BCs, lower sub-laminate was restrained for translation in x, y and z direction (uix, uiy, uiz = 0). In order to induce crack propagation, the displacement δ in z direction was put to the upper sub-laminate (uiz = δ). Additionally, the displacement degrees of freedom in x and y direction for upper branch were also restrained.

6. Results and Discussion

Experimental research of delamination process were conducted for three specimens in each type of bending–twisting laminates. Experiments focused on determination influence of fiber orientation angles in particular layers of BT specimens on initial values of mode I strain energy release rate values. Experimental tests were validated during finite element simulations prepared by using the virtual crack closure technique in Abaqus/CAE Software. Typical load–displacement plots obtained during DCB tests and FE analysis for different types of bending–twisting laminates are presented in Figure 7 and Figure 8. In the early stage, before delamination initiation moment, experimental and numerical curves were characterized by similar stiffness, especially for the BT30 and the BT45 laminates. In case of the BT60 specimens, these discrepancies were greater both for experimental and numerical P-δ curves. It could be caused by difficulties with precise determination of initial crack length obtained after the pre-cracking procedure. Although the BT60 specimens were characterized by the lowest values of non-dimensional parameters that describe the bending–twisting coupling intensity, differences in crack lengths on opposite sides of DCB beam ranged from 3 to 5 mm. It proved that initial pre-crack front was slightly non-uniform. Finally, it might lead to experimental and numerical errors. Considering definition of delamination initiation moment (determination of load and displacement point that was taken into consideration during experimental calculation of mode I fracture toughness), two different criteria were used. The first one was visual observation of cleavage movement from its initial position by using high-resolution camera. The second was measurement of different parameters of elastic waves generated by damage phenomena by using Vallen AMSY-5 acoustic emission system. Here, the first rapid growth of cumulative energy was taken into account as delamination onset. Owing to comparison of these two methods, it was found that delamination started to initiate when load–displacement curve started to deviate from linearity. This criterion is called the NL point and is clearly described in the ASTM D5528 Standard. It is also worth highlighting that because of the merging of the acoustic emission technique with monitoring of crack tip in real time, delamination onset was determined with high precision in all cases. Regarding the numerical results, it could be observed that delamination started to propagate nearly in the same point as in the experimentally load–displacement curves, in particular for bending–twisting laminates with 30° and 45° fiber orientation angles. However, for the BT60 laminate, small discrepancies between experimental and FE results were noticed. The mode I interlaminar delamination resistance values were obtained by using three different methods. Results of these calculations were collected in Table 4. The greater values of the mode I strain energy release rate determined for the BT45 laminate and were equal to approximately 0.56 N/mm. The lower values of the mode I c-SERR exhibited specimens BT30 and BT60. They were as follows: GIC = 0.51 N/mm for the BT30 and GIC = 0.41 N/mm for the BT60. For all types of laminates, discrepancies between calculated values of delamination resistance by using different methods were insignificant.
Considering delamination propagation stage, considerable differences between experimentally and numerically obtained P-δ curves can be observed. For the latter case, with growth of applied displacement load values uniformly decreased up to δ = 20 mm. For the experimentally tested DCB specimens, load trends were affected by the effect of fiber bridging. In all cases, the differences between experimental and FE plots were around 30 N. In addition, sharp load jumps were observed. This phenomenon was related to merging of singular fibers to the bundles. During debonding of such bundles, a large quantity of energy released that produced observable violent loads drop. Effect of fiber bridging was observed during experimental research for all tested laminates. One of the examples of this phenomenon is presented in Figure 9.
As it was proved in [34], extensive fiber bridging strongly affected the values of load and fracture energy during crack propagation. In this work, de Moura et al. proved that bridging phenomenon was responsible for up to 60% of energy being dissipated during delamination process. It was assessed by mechanical cutting of bridged fibers when crack propagated during opening–loading. This procedure allowed them to obtain a good agreement between the experimental and numerical load–displacement curves. With reference to current research, it could explain the differences in P-δ plots obtained for elastically coupled laminates. Propagation values of fracture toughness (R-curves) were presented in Figure 10 and Figure 11. The mode I SERR, GIprop, were obtained based on experimental data by using the modified beam theory. For initial propagation stage (0 mm < a-a0 < 10 mm), nearly for all cases, propagation values of fracture toughness were in the range between 0.4 N/mm and 0.8 N/mm. With growth of crack length, strain energy release rates significantly increased, reaching values from 0.8 N/mm up to 1.3 N/mm. In case of the BT60 laminates, for crack length between 15 mm and 25 mm, a short plateau could be observed. Here, the GIprop values were in the range between 0.9 N/mm and 1.1 N/mm. Preferable considerations about mode I propagation fracture toughness could be conducted by comparison average values of SERR obtained for BT30, BT45 and BT60 laminates. It could be noticed that for BT30 specimens, average values of GIprop started stabilizing on the level 0.85 N/mm at crack length equal to 15mm. Similar situation can be observed for BT60 laminate, but in this case, the mode I SERR values reached greater level around 1.1 N/mm. On the other hand, for BT45 specimens, propagation fracture toughness values seemed to continuously grow and any pronounced plateau was not observed. Significant differences in mode I strain energy release rate propagation values could mainly come from unstable delamination propagation caused by extensive fiber bridging.
Effect of this phenomenon was substantiated by comparing differences between total propagation crack lengths obtained experimentally as well as numerically. Note that for both cases applied displacement in z direction was 20 mm. Figure 12 presents entire crack lengths determined during finite element simulations conducted on various types of bending–twisting specimens. With regard to experimentally obtained delamination lengths, values of a were as follows: for the BT30 (a = 85 mm), for the BT45 (a = 75 mm) and for the BT60 (a = 78 mm). To sum up, in all cases, total real crack lengths were around 30 mm shorter than those obtained numerically.
The last part of numerical research covered evaluation of contributions of mix-modes in delamination process. In Table 5, the collected ratio of singular modes versus total fracture toughness is presented in percentage. It should be underlined that for all elastically coupled laminates contribution of the mode I was around 99%. The remaining 1% related to the mode II and III fracture schemes reaching the greatest values near beam model edges. Such small percentage contributions of the in-plane shear and the transverse shear modes were not expected in case of laminates with bending–twisting couplings. On the other hand, restrained displacements of upper branch in x and y directions could enforce global mode I circumstances. In reality, global mode I conditions were obtained by fixing piano hinges in tensile testing machine grips.

7. Conclusions

Mode I experimental and numerical research were conducted on carbon/epoxy composite laminates exhibiting the bending–twisting elastic couplings. Three groups of specimens with different fiber orientation angles were distinguished, namely BT30, BT45 and BT60. Experimental tests were prepared according to recommendations described in the ASTM D5528 Standard. Combination of visual recording of crack tip with acoustic emission technique ensured high precision of detecting the initiation delamination point. Values of critical strain energy release rate were determined by using three different methods of calculations: the direct beam theory, the compliance calibration method and the modified compliance calibration. Numerical simulations were prepared in Abaqus/CAE environment by using the virtual crack closure technique. The DCB model was built from SC8R continuum shell elements. Obtained outcomes revealed that:
  • The greatest value of the GIC equal to 0.56 N/mm was obtained for BT45 laminate. On the other hand, the lowest mode I c-SERR determined for BT60 laminate was equal to 0.41 N/mm. In case of the BT30 specimens, the mode I fracture toughness reached 0.51 N/mm. Moreover, the differences between delamination resistance values determined by using three different methods were slight.
  • For all tested laminates, an extensive fiber bridging was observed during delamination process. It influenced the unstable crack propagation and the discrepancies in propagation strain energy release rate values.
  • Numerically obtained load versus displacement plots were in agreement with experimental data only in initial stage of delamination. For propagation stage, the fiber bridging phenomenon influenced on differences between numerically and experimentally obtained P-δ curves. Nevertheless, application of the VCCT technique to simulate delamination in elastically coupled laminates was found out to be a sufficient tool, in particular for modeling the initial stage of delamination.
  • Prevailing contribution of the mode I fracture scheme in delamination process was affected by the constraints imposed by the piano hinges.
To conclude, elastic couplings strongly influenced the mode I delamination process. Therefore, it is important to recognize the effect of these phenomena on delamination under the in-plane shear and especially under the transverse shear fracture modes. It seems to be still unexplored area of mechanics of composite materials and will be an object of future research.

Author Contributions

Conceptualization, J.R. and S.S.; Formal analysis, J.R.; Investigation, J.R.; Methodology, J.R. and S.S.; Supervision, S.S.; Validation, S.S and J.R.; Visualization, J.R.; Writing—original draft, J.R.; Writing—review & editing, S.S. and J.R. All authors have read and agreed to the published version of the manuscript.

Funding

The research was financed in the framework of the project Lublin University of Technology-Regional Excellence Initiative, funded by the Polish Ministry of Science and Higher Education (contract no. 030/RID/2018/19).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Cañas, J.; Justo, J.; París, F. Evaluation of the interlaminar fracture toughness on composite materials using DCB test on symmetric and unsymmetric configurations. Compos. Struct. 2022, 297, 115944. [Google Scholar] [CrossRef]
  2. Gong, Y.; Chen, X.; Li, W.; Zhao, L.; Tao, J.; Zhang, J.; Hu, N. Delamination in carbon fiber epoxy DCB laminates with different stacking sequences: R-curve behavior and bridging traction-separation relation. Compos. Struct. 2021, 262, 113605. [Google Scholar] [CrossRef]
  3. Gong, Y.; Chen, X.; Tao, J.; Zhao, L.; Zhang, J.; Hu, N. A simple procedure for determining the mode I bridging stress of composite DCB laminates without measuring the crack opening displacement. Compos. Struct. 2020, 243, 112147. [Google Scholar] [CrossRef]
  4. Heidari-Rarani, M.; Sayedain, M. Finite element modeling strategies for 2D and 3D delamination propagation in composite DCB specimens using VCCT, CZM and XFEM approaches. Theor. Appl. Fract. Mech. 2019, 103, 102246. [Google Scholar] [CrossRef]
  5. Jia, R.; Zhao, L.; Curti, R.; Gong, X. Determination of pure mode-I fracture toughness of multidirectional composite DCB specimens. Eng. Fract. Mech. 2021, 252, 107776. [Google Scholar] [CrossRef]
  6. Rekbi, F.M.L.; Hecini, M.; Khechai, A. Experimental and numerical analysis of mode-I interlaminar fracture of composite pipes. J. Braz. Soc. Mech. Sci. Eng. 2018, 40, 185. [Google Scholar] [CrossRef]
  7. Gong, Y.; Tian, D.; Cao, T.; Zhao, L.; Zhang, J.; Hu, N.; Zhang, C. An R-curve effect-included delamination growth criterion for mixed-mode I/II delamination predictions of composite laminates. Compos. Struct. 2022, 295, 115846. [Google Scholar] [CrossRef]
  8. Hu, P.; Tao, R.; Li, X.; Lubineau, G. Decomposing the coupling damage in mode I multidirectional delamination. Compos. Sci. Technol. 2022, 229, 109684. [Google Scholar] [CrossRef]
  9. Low, K.O.; Johar, M.; Sung, A.N.; Mohd Nasir, M.N.; Rahimian Koloor, S.S.; Petrů, M.; AhmadIsrar, H.; JyeWong, K. Displacement rate effects on mixed-mode I/II delamination of laminated carbon/epoxy composites. Polym. Test. 2022, 108, 107512. [Google Scholar] [CrossRef]
  10. D30 Committee. Test Method for Mode I Interlaminar Fracture Toughness of Unidirectional Fiber-Reinforced Polymer Matrix Composites; ASTM International: West Conshohocken, PA, USA, 2013. [Google Scholar]
  11. Fibre-Reinforced Plastic Composites. Determination of Mode I Interlaminar Fracture Toughness, Gic, for Unidirectionally Reinforced Materials; BSI British Standards: London, UK, 2016. [Google Scholar]
  12. Valvo, P.S. On the calculation of energy release rate and mode mixity in delaminated laminated beams. Eng. Fract. Mech. 2016, 165, 114–139. [Google Scholar] [CrossRef]
  13. 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]
  14. Fukunaga, H.; Sekine, H. A Laminate Design for Elastic Properties of Symmetric Laminates with Extension-Shear or Bending-Twisting Coupling. J. Compos. Mater. 1994, 28, 708–731. [Google Scholar] [CrossRef]
  15. Baker, D.J. Response of Damaged and Undamaged Tailored Extension-Shear-Coupled Composite Panels. J. Aircr. 2006, 43, 517–527. [Google Scholar] [CrossRef]
  16. 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]
  17. 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]
  18. Rzeczkowski, J.; Samborski, S.; Valvo, P.S. Effect of stiffness matrices terms on delamination front shape in laminates with elastic couplings. Compos. Struct. 2020, 233, 111547. [Google Scholar] [CrossRef]
  19. Samborski, S. Prediction of delamination front’s advancement direction in the CFRP laminates with mechanical couplings subjected to different fracture toughness tests. Compos. Struct. 2018, 202, 643–650. [Google Scholar] [CrossRef]
  20. 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]
  21. Davidson, B.D.; Krüger, R.; König, M. Effect of stacking sequence on energy release rate distributions in multidirectional dcb and enf specimens. Eng. Fract. Mech. 1996, 55, 557–569. [Google Scholar] [CrossRef]
  22. Alhussein, M.; Ameer, Z.; Hamzah, A. Simulation of the effect laminate sequence on delamination mode-I with elastic couplings. Diagnostyka 2022, 23, 1–12. [Google Scholar] [CrossRef]
  23. 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]
  24. 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]
  25. Sun, C.T.; Zheng, S. Delamination characteristics of double-cantilever beam and end-notched flexure composite specimens. Compos. Sci. Technol. 1996, 56, 451–459. [Google Scholar] [CrossRef]
  26. D30 Committee. Test Method for Tensile Properties of Polymer Matrix Composite Materials; ASTM International: West Conshohocken, PA, USA, 2014. [Google Scholar]
  27. D30 Committee. Test Method for In-Plane Shear Response of Polymer Matrix Composite Materials by Tensile Test of a 45 Laminate; ASTM International: West Conshohocken, PA, USA, 2018. [Google Scholar]
  28. Paul, S.C.; Pirskawetz, S.; van Zijl, G.P.A.G.; Schmidt, W. Acoustic emission for characterising the crack propagation in strain-hardening cement-based composites (SHCC). Cem. Concr. Res. 2015, 69, 19–24. [Google Scholar] [CrossRef]
  29. Huo, X.; Luo, Q.; Li, Q.; Zheng, G.; Sun, G. On characterization of cohesive zone model (CZM) based upon digital image correlation (DIC) method. Int. J. Mech. Sci. 2022, 215, 106921. [Google Scholar] [CrossRef]
  30. Din, I.U.; Naresh, K.; Umer, R.; Khan, K.A.; Drzal, L.T.; Haq, M.; Cantwell, W.J. Processing and out-of-plane properties of composites with embedded graphene paper for EMI shielding applications. Compos. Part A Appl. Sci. Manuf. 2020, 134, 105901. [Google Scholar] [CrossRef]
  31. Rzeczkowski, J. An experimental analysis of the end-notched flexure composite laminates beams with elastic couplings. Contin. Mech. Thermodyn. 2021, 33, 2331–2343. [Google Scholar] [CrossRef]
  32. Rzeczkowski, J.; Paśnik, J.; Samborski, S. Mode III numerical analysis of composite laminates with elastic couplings in split cantilever beam configuration. Compos. Struct. 2021, 265, 113751. [Google Scholar] [CrossRef]
  33. 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. [Google Scholar] [CrossRef]
  34. De Moura MF, S.F.; Campilho, R.D.; Amaro, A.M.; Reis PN, B. Interlaminar and intralaminar fracture characterization of composites under mode I loading. Compos. Struct. 2010, 92, 144–149. [Google Scholar] [CrossRef]
Figure 1. Main dimensions of experimentally tested laminates.
Figure 1. Main dimensions of experimentally tested laminates.
Materials 15 07745 g001
Figure 2. Acoustic emission test set up.
Figure 2. Acoustic emission test set up.
Materials 15 07745 g002
Figure 3. Application of acoustic emission technique combined with visual monitoring of crack tip to detection of delamination initiation point during one of the DCB test conducted on BT45 laminate.
Figure 3. Application of acoustic emission technique combined with visual monitoring of crack tip to detection of delamination initiation point during one of the DCB test conducted on BT45 laminate.
Materials 15 07745 g003
Figure 4. Numerical model of the DCB beam.
Figure 4. Numerical model of the DCB beam.
Materials 15 07745 g004
Figure 5. Two-dimensional representation of the VCCT technique.
Figure 5. Two-dimensional representation of the VCCT technique.
Materials 15 07745 g005
Figure 6. Initial boundary conditions.
Figure 6. Initial boundary conditions.
Materials 15 07745 g006
Figure 7. Experimental and numerical load versus displacement plots obtained for: (a) the BT30 and (b) the BT45 bending–twisting laminates.
Figure 7. Experimental and numerical load versus displacement plots obtained for: (a) the BT30 and (b) the BT45 bending–twisting laminates.
Materials 15 07745 g007
Figure 8. Experimental and numerical load versus displacement plots obtained for the BT60 bending–twisting laminates.
Figure 8. Experimental and numerical load versus displacement plots obtained for the BT60 bending–twisting laminates.
Materials 15 07745 g008
Figure 9. Fiber bridging phenomenon observed during DCB experiments conducted on the BT30 and the BT45 laminates.
Figure 9. Fiber bridging phenomenon observed during DCB experiments conducted on the BT30 and the BT45 laminates.
Materials 15 07745 g009
Figure 10. R-curves obtained for (a) the BT30 and (b) the BT45 laminates.
Figure 10. R-curves obtained for (a) the BT30 and (b) the BT45 laminates.
Materials 15 07745 g010
Figure 11. R-curves obtained for the BT60 laminates.
Figure 11. R-curves obtained for the BT60 laminates.
Materials 15 07745 g011
Figure 12. Total crack lengths obtained during the FE simulations conducted on different types of BT laminates.
Figure 12. Total crack lengths obtained during the FE simulations conducted on different types of BT laminates.
Materials 15 07745 g012
Table 1. Material properties of carbon/epoxy composite laminates.
Table 1. Material properties of carbon/epoxy composite laminates.
E1 [GPa]S.D. [%]E2 [GPa]S.D. [%]ν12S.D. [%]G [GPa]S.D. [%]
112.1053.567.4212.610.278.283.3382.18
Table 2. Experimental elastically coupled laminates with different fiber orientation angles.
Table 2. Experimental elastically coupled laminates with different fiber orientation angles.
SpecimensStacking Sequence (for One Branch)Bt [-]Dc [-]
BT3030/0°/30/30/0°/−30/0°/−30/−30/−30/−30/0°/−30/30/0°/0°/30/300.22380.2208
BT4545/0°/45/45/0°/−45/0°/−45/−45/−45/−45/0°/−45/45/0°/0°/45/450.23630.2384
BT6060/0°/60/60/0°/−60/0°/−60/−60/−60/−60/0°/−60/60/0°/0°/60/600.13410.0990
Table 3. Fracture and material parameters implemented during FE simulations.
Table 3. Fracture and material parameters implemented during FE simulations.
E11
[MPa]
E22 = E33
[MPa]
G12 = G13
[MPa]
G23
[MPa]
ν12 = ν13
[-]
ν23
[-]
GIC
[N/mm]
GIIC
[N/mm]
GIIIC
[N/mm]
η
[-]
BT30112,1057421333827690.270.340.510.70.71.685
BT450.56
BT600.41
Table 4. Values of critical strain energy release rates determined for bending–twisting coupled laminates by using different data reduction schemes.
Table 4. Values of critical strain energy release rates determined for bending–twisting coupled laminates by using different data reduction schemes.
LaminateBT30BT45BT60
GIC [N/mm]GIC [N/mm]GIC [N/mm]
CCMMBTMCCCCMMBTMCCCCMMBTMCC
10.480.480.470.560.550.550.390.390.39
20.520.520.540.620.620.620.450.450.45
30.530.520.520.530.520.510.390.390.39
average0.510.510.510.570.560.560.410.410.41
Table 5. Mix-modes contributions.
Table 5. Mix-modes contributions.
BT30BT45BT60
GI/GT99 %99%99%
GII/GT0.6%0.8%0.85%
GIII/GT0.4%0.2%0.15%
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. Experimental and Numerical Research of Delamination Process in CFRP Laminates with Bending-Twisting Elastic Couplings. Materials 2022, 15, 7745. https://doi.org/10.3390/ma15217745

AMA Style

Rzeczkowski J, Samborski S. Experimental and Numerical Research of Delamination Process in CFRP Laminates with Bending-Twisting Elastic Couplings. Materials. 2022; 15(21):7745. https://doi.org/10.3390/ma15217745

Chicago/Turabian Style

Rzeczkowski, Jakub, and Sylwester Samborski. 2022. "Experimental and Numerical Research of Delamination Process in CFRP Laminates with Bending-Twisting Elastic Couplings" Materials 15, no. 21: 7745. https://doi.org/10.3390/ma15217745

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