Next Article in Journal
Spinel rGO Wrapped CoV2O4 Nanocomposite as a Novel Anode Material for Sodium-Ion Batteries
Next Article in Special Issue
Influence of Nanoparticles on Thermal and Electrical Conductivity of Composites
Previous Article in Journal
Thermosensitive Hydrogel for Encapsulation and Controlled Release of Biocontrol Agents to Prevent Peanut Aflatoxin Contamination
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Experimental and Numerical Assessment of Fibre Bridging Toughening Effects on the Compressive Behaviour of Delaminated Composite Plates

1
Department of Engineering, University of Campania “L. Vanvitelli”-via Roma, 81031 Aversa, Italy
2
CIRA Italian Aerospace Research Centre, Via Maiorise s/n, 81043 Capua, Italy
3
Institute of Polymers, Composites and Biomaterials, CNR–Research National Council of Italy, Granatello 80055 Portici, Naples, Italy
*
Author to whom correspondence should be addressed.
Submission received: 14 February 2020 / Revised: 26 February 2020 / Accepted: 27 February 2020 / Published: 3 March 2020
(This article belongs to the Special Issue Damage Mechanics of Polymer Composites)

Abstract

:
Increasing the Mode I inter-laminar fracture toughness of composite laminates can contribute to slowing down delamination growth phenomena, which can be considered one of the most critical damage mechanisms in composite structures. Actually, the Mode I interlaminar fracture toughness (GIc) in fibre-reinforced composite materials has been found to considerably increase with the crack length when the fibre bridging phenomenon takes place. Hence, in this paper, the fibre bridging phenomenon has been considered as a natural toughening mechanism able to replace embedded metallic or composite reinforcements, currently used to increase tolerance to inter-laminar damage. An experimental/numerical study on the influence of delamination growth on the compressive behaviour of fibre-reinforced composites characterised by high sensitivity to the fibre bridging phenomenon has been performed. Coupons, made of material systems characterised by a variable toughness related to a high sensitivity to the fibre bridging phenomenon and containing artificial through-the-width delaminations, were subjected to a compressive mechanical test and compared to coupons made of standard material system with constant toughness. Out-of-plane displacements and strains were monitored during the compression test by means of strain gauges and digital image correlation to assess the influence of fibre bridging on delamination buckling, delamination growth and on the global buckling of the specimens, including buckling shape changes. Experimental data were combined with a numerical study, performed by means of a virtual crack closure technique based procedure, named SMart Time XB – Fibre Bridging (SMXB-FB), able to mimic the crack bridging effect on the toughness properties of the material system. The combination of numerical results and experimental data has allowed the deformations and the buckling shape changes to be correlated to the onset and evolution of damage and, hence, contributes to improving the knowledge on the interaction of the failure mechanisms in the investigated composite specimens.

1. Introduction

Fibre bridging is a toughening mechanism which is able to delay delamination growth in composite structures [1,2,3]. Indeed, delaminations are among the most dangerous failure mechanisms for composite material structures, due to their frequent unstable growth, which can lead to the sudden and premature collapse of structures [4,5,6,7,8]. Actually, fracture toughness is a non-isotropic matrix property able to slow down damage progression in the weaker (out-of-plane) direction of fibre-reinforced composite laminates. Hence, improving this material property becomes mandatory to widen the possible applications of fibre-reinforced composite laminates.
To improve the fracture toughness, toughening resins and thickness reinforcements (stitching, pinning or orthogonal weaving) are generally adopted [9,10,11,12,13]. Such reinforcements, even if able to slow down and sometimes arrest the delamination advance, can strongly reduce the in-plane mechanical properties of the fibre-reinforced composites. Indeed, in the stitching process a needle is used to pierce the laminate, in order to insert a high tensile strength yarn. Literature study [14] demonstrates that stitching can cause breakage, kink and spread of the fibres. Furthermore, the formation of resin-rich regions, porosity and resin cracks can occur. Likewise, the presence of the z-pins creates resin pockets and skew of the fibres, leading to degradation of the in-plane strength of the composite [15]. Another highly effective method to improve the fracture toughness is the addition of nanoparticles in the polymeric matrix [16,17,18,19,20,21], which improves also tensile strength, stiffness and thermal properties. On the other hand, some disadvantages are associated with nanoparticles, such as viscosity increase, dispersion problems and sedimentation [22]. Hence, attention should be paid to alternative toughening approaches based on the modification of the manufacturing process of composites that are able to tune the sensitivity to natural toughening phenomena, such as the fibre bridging.
In order to experimentally assess the effects of toughness variations on the mechanical behaviour of composite components, Non-Destructive Testing (NDT) methods can be used, which, replacing the more uncertain “destructive sample testing”, allow the evolution of inter-laminar damage under loading conditions to be evaluated. Among others, Digital Image Correlation (DIC), which is a technique based on the acquisition of surface digital images at different loading stages, allows the strain and the displacement distributions all over the specimen to be evaluated. The analysis and comparison of the acquired images can give an idea of the damage evolution in the inspected structures [23,24,25,26]. The main advantage of the DIC is the ease of use, since significantly less equipment is needed for DIC inspections compared to other NDT techniques, such as the ones based on ultrasound [27,28] and thermography [29,30].
In this paper, a numerical/experimental study, on the compressive behaviour of delaminated composite plates, is introduced. The main aim is to investigate the effects of modifications to the manufacturing process of composites for tuning the sensitivity to fibre bridging. Two material systems have been considered: the former is a material highly sensitive to fibre bridging characterised by a toughness increase with crack length; the latter is a toughened material with constant toughness. Epoxy resin/carbon fibre coupons, characterized by an artificial through-the-width delamination, have been experimentally tested under compression to assess their inter-laminar damage behaviour. The specimens have been equipped with strain gauges to investigate the delamination buckling and growth, as a result of the compressive load. Moreover, a DIC system has been employed for the monitoring of displacements and deformation fields during the test.
A numerical study has also been carried out and the results have been used, together with the experimental data, to provide an interpretation of the effects of the Mode I fracture toughness on delamination phenomena. The adopted numerical tool, validated in [31], is an innovative and robust Virtual Crack Closure Technique (VCCT)-based procedure, named SMXB-FB, able to simulate the interlaminar damage in composite material structure, taking into account the bridging of the fibres, as well as the resistance curve behaviour (variation of toughness with crack length).
In Section 2, the manufacturing of the investigated plates is analysed in detail and information is provided on the experimental procedures behind the compressive mechanical tests. In Section 3, the SMXB-FB numerical tool and the finite element model are introduced. Finally, in Section 4, the experimental data and the numerical results are presented, compared and discussed in order to improve the knowledge about the interactions between fibre bridging and buckling behaviour of the investigated plates.

2. Specimen Manufacturing and Experiments Details

Epoxy resin/carbon fibre specimens with a 50-mm wide trough-the-width artificial delamination were manufactured according to the geometrical description given in Figure 1. Coupons are made of 24 plies (each ply is 0.1875 mm thick) with a stacking sequence of [0,90]6s. A rectangular Teflon film was inserted between the 4th and 5th ply, as shown in Figure 1, to create an artificial delamination. Tabs on both sides of the samples were placed to allow a distribution as uniform as possible of the load from the test machine to the sample.
Two different curing processes have been used to manufacture the epoxy resin/carbon fibre plates. As a result, two sets of coupons have been obtained with the same mechanical properties but different Mode I fracture toughness:
  • set of coupons characterised by a constant GIc of 510 J/m2 (toughened material);
  • set of coupons characterized by a variable GIc, from 243 J/m2 to 456 J/m2 depending on delamination size (material sensitive to fibre bridging).
The resistance curves for the two sets of coupons, representing the GIc as a function of the crack length a, are reported in Figure 2, together with the other mechanical properties.
Three samples were cut from each plate, resulting in three specimens made of toughened material and three specimens made of material sensitive to fibre bridging. In Figure 3, the stacking sequence and the position of delamination along the thickness, between the fourth and fifth ply, are introduced.
The tested coupons were labelled according to their fracture toughness characteristics, i.e., toughened material (TOUGH) and material sensitive to fibre bridging (MFB), as listed in Table 1, where the tests matrix is reported.
According to Table 1, the TOUGH#1, TOUGH#2, MFB#1 and MFB#2 specimens were equipped with two strain gauge rosettes (SG1 and SG2), located, as shown in Figure 4a, on the thicker sub-laminate side, to measure the deformations in the 0° and 90° fibre directions. The remaining samples, TOUGH#3 and MFB#3, were prepared for the DIC by painting in white and by randomly sprinkling with black pigments the thinner sub-laminate surfaces. These two last specimens were also equipped with one 0° fibre direction strain gauge (SG1bis), placed, as shown in Figure 4b, on their thicker side.
Images of the specimens are shown in Figure 5. Indeed, in Figure 5a the specimen prepared for strain monitoring with strain gauges is shown, while in Figure 5b the specimen prepared for DIC is presented.
Before the compression tests, the samples were checked for manufacturing defects using ultrasonic inspections, performed by means of an ultrasonic device with an array of 64 sensors. No relevant defects have been found apart from some small voids at the edges of the Teflon film, probably due to the roughness of the Teflon film self. Figure 6 shows the C-scans and the B-scans of the TOUGH#1 and MFB#1 samples, as an example of the performed non-destructive ultrasonic checks.
A hydraulic testing machine was used to perform compressive mechanical tests. The tabs of the specimens were clamped in the hydraulic grip and a controlled compressive displacement was applied with a rate of 0.5 mm/min. A picture of the specimen mounted in the testing rig is shown in Figure 7. After the installation of each specimen, the strain gauges were linked to the data acquisition systems and calibrated. Consequently, the specimens were loaded with a small force before the test initiation in order to check the strain gauge measurements and verify the correct alignment and the proper introduction of the compressive load.
According to the test matrix of Table 1, DIC inspections were performed on TOUGH#3 and MFB#3 samples. DIC is a no-contact measurement technique based on numerical processing of digital images for the analysis of displacement and deformation fields. Through two cameras, several images are acquired, as well as an initial image representing the body in the reference condition, which is typically the undeformed configuration. The subsequent data analysis operations compare the reference image with those acquired with the body in the different deformed conditions, in order to calculate the relative displacements and deformations. In more detail, the area of interest is divided by the processing software into small square regions, named a subset. In each subset, the displacement of the central point is calculated for different time instants.
Since calculations depend on the subset correspondence within the acquired images, it is clear that the choice of suitable patterns becomes of main relevance. The special pattern used for DIC analysis is called a speckle pattern, which is typically composed of dark speckles of uniform size, distributed on a white background, to maximize contrast, as schematically shown in Figure 8.
Proper correlation criteria have been defined to find the correspondence between each subset in the reference image and in the images acquired during the deformed states (cross correlation criteria and sum-squared differences-based criteria [32,33]).
In Figure 9, the DIC instrumentation used during the tests is shown. The digital images were acquired every 3.5 s for the entire duration of the tests.

3. Numerical Procedure Description and Finite Elements Model Definition

The experimental research activity, described in the previous section, was numerically simulated with the aim to investigate the damage mechanisms onset, evolution and interaction when delamination grows under compression, depending on the different interlaminar toughness conditions. The robust numerical tool SMXB FB, validated in [33], was used to mimic the compressive behaviour of the abovementioned delaminated composite coupons taking into account the fibre bridging effects.

3.1. SMXB FB Numerical Tool Description

The delamination growth phenomenon is usually simulated in commercial Finite Elements Method (FEM) codes by connecting the propagation region nodes’ couple with appropriate contact elements and by defining a propagation criterion able to “kill” such interface elements when proper conditions on the Strain Energy Release Rate (SERR) are met. This method is labelled as Fail Release (FR) approach. The SERR can be calculated by using the Virtual Crack Closure Technique (VCCT) equations, described in detail in [34]. In particular, the VCCT equation for the SERR calculation in the case of 4-noded solid elements is reported in Equation (1).
G j = F j Δ u j 2 Δ A
where the subscript j identifies the opening (I), sliding (II) and tearing (III) fracture mode, Fj the force at the crack tip, Δuj the opening displacement and ΔA the delaminated area.
The SMXB-FB numerical tool, implemented in the framework of the ANSYS® FEM software [35], is based on the VCCT-FR approach and employs the linear power law, in Equation (2), as the delamination growth criterion.
E d = G I G I c + G I I G I I c + G I I I G I I I c 1
According to Equation (2), the SERR values (Gj with j = I, II, III) are compared with the corresponding critical values (Gjc with j = I, II, III), which are experimentally determined [36,37,38].
The standard VCCT-FR approach usually underestimates or overestimates the delaminated area, due to mesh size and time step dependence, as demonstrated in [39]. The SMXB-FB procedure, is capable of overcoming these dependencies and providing a correct evaluation of the delaminated area as a function of the applied load.
According to Equation (2), when the condition Ed > 1 is verified for a certain load step size, the commercial VCCT-based codes release the constraint between the pair of nodes and define the delaminated zone (ANUM), which depends on the element size (ΔAe). In such a way, an over/underestimation of the damaged area can be achieved. The VCCT-based SMXB-FB numerical procedure iteratively changes the load step size, in order to find the value of applied load for which the nominal released area (AES) exactly corresponds, when the growth criterion Ed = 1 is satisfied, to the numerically computed damaged area (ANUM). This principle is clarified in Equation (3), where N is the number of nodes couples characterized by Ed = 1 and
A E S = i = 1 N ( j = 1 3 F j i Δ u j i 2 G j c ) i = 1 N Δ A i D = f · Δ A e = A N U M
Additionally, the SMXB-FB procedure is capable of correctly defining the local coordinates system at the delamination front for the computation of SEER components, even for irregular delamination front shapes.
In the SMXB-FB numerical methodology, the ability to overcome the mesh and time step dependences is combined with the capability of taking into account the fibre bridging phenomenon by considering the materials’ resistance curves. Indeed, the growth criterion in Equation (2) has been rewritten as reported in Equation (4) (where a is the crack length).
E d = G I G I c ( a ) + G I I G I I c + G I I I G I I I c 1
According to Equation (3), when the equivalence is satisfied, the constrains between the nodes are released and the delamination front shape is adapted. As the crack advances, the Fibre Bridging (FB) modulus of the SMXB-FB numerical tool can assign to each node of the new delamination front, resulting from the propagation, a new value of GIc, according to the resistance curve. The FB modulus operating principles are schematically described in Figure 10.
The above-mentioned features, such as time step and mesh independence and the ability to mimic the bridging of the fibres, make the SMXB-FB numerical tool a robust and innovative procedure able to realistically simulate the interlaminar damage evolution in composite structures.

3.2. Finite Elements Model

The finite element model was created according to the tested coupons geometry, shown in Figure 1, and discretized with solid ANSYS® layered elements, as shown in Figure 11. In order to model the region of possible delamination propagation, the couples of nodes, belonging to the two considered sublaminates (see Figure 3) were connected by means of contact elements (see Figure 12). Such contact elements are released depending on whether the propagation criterion is satisfied; this feature is called “birth and death option”.
With the aim to simulate the experimental compressive tests boundary conditions, the panel was clamped on one edge, while compressive displacement was applied to the other edge. The tabs regions were neglected in the FEM model.

4. Numerical Results/Experimental Data Comparisons and Discussion

In this section, the experimental strain gauges reading and the DIC images and measurements are compared and integrated with the SMXB-FB numerical results, in terms of strains, out-of-plane displacements and delaminated area, as a function of the compressive load in order to better understand the influence of fracture toughness on the compressive behaviour for the analysed composite plates.
All the tested specimens showed a similar qualitative behaviour when subjected to compressive load. In order to outline the main phases during the compression tests, as an example, pictures taken during the tests are shown in Figure 13 and Figure 14, respectively, for a toughened sample and for a sample sensitive to fibre bridging.
In Figure 13a and Figure 14a the local delamination buckling onset can be seen, while in Figure 13b and Figure 14b the delamination growth onset can be noticed. Figure 13c and Figure 14c clearly show the transition to delamination buckling in the first global buckling shape, which develops in the opposite out-of-plane direction with respect to the delamination buckling. Finally, Figure 13d and Figure 14d show the transition to the second buckling shape (snap-through buckling phenomenon [40,41]) caused by a sudden delamination growth.
The fibre bridging, occurring in the coupons sensitive to fibre bridging, can be seen in Figure 14d (zoomed view).
Figure 15 introduces the experimental global compressive behaviour of the coupons in terms of load vs displacements curves. Indeed, in Figure 15 a change of slope, highlighting the global buckling phenomenon, can be observed at about 0.5 mm of applied displacement, corresponding to a compressive load of 45.2 kN in the case of toughened material and 41.2 kN in the case of material sensitive to fibre bridging. Moreover, the complete debonding state and the snap-through buckling phenomenon can be appreciated at 2.5 mm (40.4 kN) and 2 mm (39.3 kN) of applied displacement, respectively, in the case of toughened material and in the case of material sensitive to fibre bridging.
In Table 2, a summary of the loads and strains at local delamination buckling, global buckling and snap-through is provided.
From Table 2, it is clear that, as seen in the previous figure, the global buckling and the snap-through buckling are influenced by the fracture toughness characteristics of the material, while, since delamination growth has not yet started, the delamination buckling is, obviously, not influenced by the fracture toughness.
In Figure 16, the numerical strains in the loading direction as a function of the compressive load are compared to the experimental strain gauge measurements for the toughened specimens and the specimens sensitive to fibre bridging.
According to Figure 16, for both the analysed materials, the bending stiffness and the maximum achieved load are excellently predicted by the SMXB-FB routine. In Figure 16, it is also possible to notice that the snap-through phenomenon is delayed in the case of toughened material and this effect is well predicted by the numerical tool used for the simulations. Indeed, the lower toughness of material sensitive to fibre bridging at delamination growth initiation causes an acceleration of the delamination growth phenomenon with consequent earlier conditions of elastic instability leading to the snap-through event, if compared to the toughened material.
However, the ultimate strains values are slightly underestimated by the numerical code for both the material configurations. Indeed, with the displacement control method, which has been used in this work, the performed static analyses can provide incomplete or misguided information about the stability of the sample in the presence of dynamic phenomena, such as snap-through buckling. Anyway, even in the presence of dynamic snap-through buckling, the trend of numerical results is in excellent agreement with that of experimental outputs, mainly due to the almost quasi-static behaviour of the sample up to the snap-through buckling.
Such a discrepancy in the ultimate deformation values is also pointed out in Figure 17, where the strains measured by means of the DIC system at the P1 location (shown in Figure 17) are compared with the numerical strains at the same location. Indeed, at higher strain values, misalignment between numerical and experimental results can be seen because of the presence of dynamic phenomena, which were neglected in the numerical analyses.
In Figure 18, the load vs out-of-plane displacements curve, experimentally obtained by the DIC system, for the toughened material and material sensitive to fibre bridging, are compared to the numerical results. According to Figure 18, the out-of-plane displacement at the central location (Point 1) is predicted with excellent accuracy, while the experimentally determined out-of-plane displacements at Point 2 and Point 3 locations are slightly overestimated by the numerical tool. Actually, as seen from the deformed shapes, the dynamic effects due to the buckling shape change are localised at Point 2 and Point 3.
The growth of the delaminated area was analysed for the two material configurations, both experimentally and numerically during the loading process, providing interesting indications of the GIc influence on the delamination evolution. Numerical results and experimental data in terms of delamination size as a function of the compressive load are presented, compared and integrated in Figure 19, for the two analysed material configurations.
Figure 19 shows that, for both the analysed material configurations, the numerical predictions give a valuable contribution to the clarification of experimental data. Indeed, the trends of the delaminated area vs compressive load point out that, for the material sensitive to fibre bridging, the delamination growth is anticipated (32 kN) with respect to the toughened configuration (43 kN). For the material sensitive to fibre bridging configurations, after a stable delamination opening, an unstable propagation occurs, leading to the snap-through phenomenon. On the contrary, the toughened material is characterised by a highly unstable delamination growth since delamination growth initiation occurs very close to the global buckling. As already pointed out, the numerical tool adopted for the analyses is not able to account for the dynamic effects, hence it is not able to give a prediction of the delamination growth beyond the snap-through phenomenon.
As already mentioned, digital image correlation inspections were performed on two of the tested samples (TOUGH#3 and MFB#3). Figure 20 shows the out-of-plane displacement contour plot obtained by means of the DIC technique, as an example, for the case MFB#3.
As for Figure 13 and Figure 14, in Figure 20a the local delamination buckling onset can be seen, while in Figure 20b the delamination growth onset can be pointed out. Figure 20c clearly shows the transition from delamination buckling to the first global buckling shape which develops in the opposite out-of-plane direction with respect to the delamination buckling. Finally, Figure 20d shows the transition to the second buckling shape caused by a sudden and unstable delamination growth phenomenon.
Figure 21 and Figure 22 show the numerically deformed shapes with the out-of-plane displacement contour plot, evaluated at four fixed values of the applied displacement, respectively, for the sample characterized by toughened material and the sample characterized by material sensitive to fibre bridging.
Figure 21 and Figure 22 point out the delay in the delamination growth phenomenon and, hence, in the snap-through phenomenon observed for the toughened material configuration, characterised by a constantly higher GIc, with respect to the material sensitive to fibre bridging configuration. Actually, the configuration made of material sensitive to fibre bridging (MFB) is characterised by a lower starting value of the fracture toughness, if compared to the toughened material (TOUGH). This leads to an earlier, even if more stable, delamination growth in the first case (MFB) with respect to the case of toughened material (TOUGH), whose delamination growth initiates with a highly unstable behaviour at higher loads (close to the global buckling load).
Finally, in Figure 23, a comparison between experimental and deformed shapes at fixed values of the applied displacement, for a specimen characterised by sensitivity to fibre bridging, is introduced. The excellent agreement between experimental and numerical delamination shapes demonstrates the capability of the adopted numerical tool in predicting the compressive behaviour of the analysed delaminated composite plates, even in the presence of highly dynamic phenomena, such as the snap-through buckling.

5. Conclusions

An experimental/numerical study has been carried out to investigate the interlaminar damage behaviour of delaminated composite plates characterized by different fracture toughness under compressive loading conditions. Two fracture toughness configurations have been analysed: the toughened material configuration, characterised by a constantly high fracture toughness and the configuration sensitive to fibre bridging, characterised by an increasing value of fracture toughness. Actually, an experimental testing campaign was performed to analyse the compressive behaviour of the specimens and the DIC technique was used in order to measure samples displacements and deformations during the tests. Moreover, all the samples were equipped with strain gauges, in different locations, to measure deformations. To numerically simulate the mechanical behaviour of compressed specimen, a numerical tool able to consider the fibre bridging effects was adopted.
The experimental data and the numerical results have been found to be in excellent agreement in terms of deformed shapes, compressive load vs strain curves and out-of-plane displacements.
The two analysed configurations showed almost the same compressive behaviour characterised by the local delamination buckling onset, delamination growth onset, global buckling of the plate and transition from the global buckling shape to the second buckling shape (snap-through buckling phenomenon) caused by a sudden unstable delamination growth.
Differences in terms of delamination growth behaviour have been found, as expected, between the two analysed material configurations. Indeed, the coupons made of material sensitive to fibre bridging are characterised by a lower (standard) value of the fracture toughness at the beginning of the delamination growth phenomenon, which leads to an earlier, even if more stable, delamination growth, if compared to the toughened material configuration, whose delamination growth initiates with a highly unstable behaviour at higher loads (close to the global buckling load). These differences lead to a higher global buckling load for the toughened composite plates and to a delay in the snap-through phenomenon for these plates, if compared to the composite plates sensitive to fibre bridging. Actually, the snap-through phenomenon introduces dynamic effects which cannot be simulated by the tools used for the numerical prediction. Indeed, an underestimation of the experimental ultimate strains has been found; however, a tool has been found able to capture the difference between the two material configurations and the main phenomena occurring during the compression.

Author Contributions

All authors equally contributed to this work. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Conflicts of Interest

The authors declare no conflicts of interest.

References

  1. Johnson, W.S.; Mangalgiri, P.D. Investigation of fiber bridging in double cantilever beam specimens. J. Compos. Technol. Res. 1987, 9, 10–13. [Google Scholar]
  2. Nemat-Nasser, S.; Hori, M. Toughening by partial or full bridging of cracks in ceramics and fiber reinforced composites. Mech. Mater. 1987, 6, 245–269. [Google Scholar] [CrossRef]
  3. Huang, X.N.; Hull, D. Effects of fibre bridging on GIC of a unidirectional glass/epoxy composite. Compos. Sci. Technol. 1989, 35, 283–299. [Google Scholar] [CrossRef]
  4. Davim, J.P.; Reis, P. Study of delamination in drilling carbon fiber reinforced plastics (CFRP) using design experiments. Compos. Struct. 2003, 59, 481–487. [Google Scholar] [CrossRef]
  5. Schoeppner, G.A.; Abrate, S. Delamination threshold loads for low velocity impact on composite laminates. Compos. Part A Appl. Sci. Manuf. 2000, 31, 903–915. [Google Scholar] [CrossRef]
  6. Ho-Cheng, H.; Dharan, C.K.H. Delamination during drilling in composite laminates. J. Eng. Ind. 1990, 112, 236–239. [Google Scholar] [CrossRef]
  7. Mi, Y.; Crisfield, M.A.; Davies, G.A.O.; Hellweg, H.-B. Progressive delamination using interface elements. J. Compos. Mater. 1998, 32, 1246–1272. [Google Scholar] [CrossRef]
  8. Reeder, J.R.; Crews, J.H., Jr. Mixed-mode bending method for delamination testing. Aiaa J. 1990, 28, 1270–1276. [Google Scholar] [CrossRef] [Green Version]
  9. Pelstring, R.; Michael, M.; Ram, C. Stitching to improve damage tolerance of composites. Int. SAMPE Symp. Exhib. (Proc.) 1989, 34 Pt 2, 1519–1528. [Google Scholar]
  10. Dow Marvin, B.; Smith Donald, L. Damage-tolerant composite materials produced by stitching carbon fabrics. Natl. Sampe Tech. Conf. 1989, 21, 595–605. [Google Scholar]
  11. Shu, D.; Mai, Y.-W. Effect of stitching on interlaminar delamination extension in composite laminates. Compos. Sci. Technol. 1993, 49, 165–171. [Google Scholar] [CrossRef]
  12. Yasaee, M.; Lander, J.K.; Allegri, G.; Hallett, S.R. Experimental characterisation of mixed mode traction-displacement relationships for a single carbon composite Z-pin. Compos. Sci. Technol. 2014, 94, 123–131. [Google Scholar] [CrossRef] [Green Version]
  13. Dantuluri, V.; Maiti, S.; Geubelle, P.H.; Patel, R.; Kilic, H. Cohesive modeling of delamination in Z-pin reinforced composite laminates. Compos. Sci. Technol. 2007, 67, 616–631. [Google Scholar] [CrossRef]
  14. Mouritz, A.P.; Leong, K.H.; Herszberg, I. A review of the effect of stitching on the in-plane mechanical properties of fibre-reinforced polymer composites. Compos. Part A Appl. Sci. Manuf. 1997, 28, 979–991. [Google Scholar] [CrossRef]
  15. Steeves, C.A.; Fleck, N.A. In-plane properties of composite laminates with through-thickness pin reinforcement. Int. J. Solids Struct. 2006, 43, 3197–3212. [Google Scholar] [CrossRef] [Green Version]
  16. Quan, D.; Mischo, C.; Li, X.; Scarselli, G.; Ivanković, A.; Murphy, N. Improving the electrical conductivity and fracture toughness of carbon fibre/epoxy composites by interleaving MWCNT-doped thermoplastic veils. Compos. Sci. Technol. 2019, 182, 107775. [Google Scholar] [CrossRef]
  17. Rawat, P.; Singh, K.K.; Behera, R.P. Enhancement of impact properties by using multiwall carbon nanotubes as secondary reinforcement in glass/epoxy laminates. J. Test. Eval. 2020, 48, JTE20170510. [Google Scholar] [CrossRef]
  18. Lau, K.-T.; Gu, C.; Hui, D. A critical review on nanotube and nanotube/nanoclay related polymer composite materials. Compos. Part B Eng. 2006, 37, 425–436. [Google Scholar] [CrossRef]
  19. Huang, Z.-M.; Zhang, Y.-Z.; Kotaki, M.; Ramakrishna, S. A review on polymer nanofibers by electrospinning and their applications in nanocomposites. Compos. Sci. Technol. 2003, 63, 2223–2253. [Google Scholar] [CrossRef]
  20. Al-Sabagh, A.; Taha, E.; Kandil, U.; Nasr, G.-A.; Taha, M.R. Monitoring damage propagation in glass fiber composites using carbon nanofibers. Nanomaterials 2016, 6, 169. [Google Scholar] [CrossRef] [Green Version]
  21. Siddiqui, N.A.; Woo, R.S.C.; Kim, J.-K.; Leung, C.C.K.; Munir, A. Mode I interlaminar fracture behavior and mechanical properties of CFRPs with nanoclay-filled epoxy matrix. Compos. Part A Appl. Sci. Manuf. 2007, 38, 449–460. [Google Scholar] [CrossRef] [Green Version]
  22. Heidarian, M.; Shishesaz, M.R.; Kassiriha, S.M.; Nematollahi, M. Characterization of structure and corrosion resistivity of polyurethane/organoclay nanocomposite coatings prepared through an ultrasonication assisted process. Prog. Org. Coat. 2010, 68, 180–188. [Google Scholar] [CrossRef]
  23. Peters, W.H.; Ranson, W.F. Digital imaging techniques in experimental stress analysis. Opt. Eng. 1982, 21, 427–431. [Google Scholar] [CrossRef]
  24. Chu, T.C.; Ranson, W.F.; Sutton, M.A.; Peters, W.H. Applications of digital-image-correlation techniques to experimental mechanics. Exp. Mech. 1985, 25, 232–244. [Google Scholar] [CrossRef]
  25. Sirohi, R.S. Optical Methods of Measurement: Wholefield Techniques, 2nd ed.; CRC: Boca Raton, FL, USA, 2009. [Google Scholar]
  26. Bornert, M.; Brémand, F.; Doumalin, P.; Dupré, J.-C.; Fazzini, M.; Grédiac, M.; Hild, F.; Mistou, S.; Molimard, J.; Orteu, J.-J.; et al. Assessment of digital image correlation measurement errors: Methodology and results. Exp. Mech. 2009, 49, 353–370. [Google Scholar] [CrossRef] [Green Version]
  27. Tsuda, H. Ultrasound and damage detection in CFRP using fiber Bragg grating sensors. Compos. Sci. Technol. 2006, 66, 676–683. [Google Scholar] [CrossRef]
  28. Roberts, R.A. Computational prediction of micro-crack induced ultrasound attenuation in CFRP composites. J. Nondestruct. Eval. 2014, 33, 443–457. [Google Scholar] [CrossRef]
  29. Maldague, X. Nondestructive Evaluation of Materials by Infrared Thermography; Springer: New York, NY, USA, 1993. [Google Scholar]
  30. Toscano, C.; Riccio, A.; Camerlingo, F.; Meola, C. Lockin thermography to monitor propagation of delamination in CFRP composites during compression tests. In Proceedings of the 11th Quantitative Infrared Thermography, Naples, Italy, 11–14 June 2012. [Google Scholar]
  31. Riccio, A.; Russo, A.; Sellitto, A.; Raimondo, A. Development and application of a numerical procedure for the simulation of the “Fibre Bridging” phenomenon in composite structures. Compos. Struct. 2017, 168, 104–119. [Google Scholar] [CrossRef]
  32. Pan, B.; Qian, K.M.; Xie, H.M.; Asundi, A. Twodimensional digital image correlation for in-plane displacement and strain measurement: A review. Meas. Sci. Technol. 2009, 20, 062001. [Google Scholar] [CrossRef]
  33. Giachetti, A. Matching techniques to compute image motion. Image Vis. Comput. 2000, 18, 247–260. [Google Scholar] [CrossRef] [Green Version]
  34. Krueger, R. The Virtual Crack Closure Technique: History, Approach and Applications. ICASE Report No. 2002-10. Appl. Mech. Rev. 2004, 57, 109–143. [Google Scholar] [CrossRef]
  35. ANSYS®. Academic Research Mechanical, Release 18.1, Help System. In Coupled Field Analysis Guide; ANSYS, Inc.: Canonsburg, PA, USA, August 2005. [Google Scholar]
  36. ASTM D5528-13, Standard Test Method for Mode I Interlaminar Fracture Toughness of Unidirectional Fiber-Reinforced Polymer Matrix Composites; ASTM International: West Conshohocken, PA, USA, 2013.
  37. ASTM D7905/D7905M-14, 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, 2014.
  38. Li, J.; Lee, S.; Lee, E.; O’Brien, T. Evaluation of the Edge Crack Torsion (ECT) Test for Mode III Interlaminar Fracture Toughness of Laminated Composites. J. Compos. Technol. Res. 1997, 19, 174–183. [Google Scholar]
  39. Pietropaoli, E.; Riccio, A. On the robustness of finite element procedures based on Virtual Crack Closure Technique and fail release approach for delamination growth phenomena. Definition and assessment of a novel methodology. Compos. Sci. Technol. 2010, 70, 1288–1300. [Google Scholar] [CrossRef]
  40. Mehdi Farshad; Plastic Pipe Systems (Eds.) Mehdi Farshad, 4—Buckling of Plastic Pipes; Elsevier Science: Amsterdam, The Netherlands, 2006; pp. 101–133. [Google Scholar]
  41. Caulfeild, S. Composite skin snap-through buckling under transient dynamic pressure loads 2nd Joint US-Canada Conference on Composites—American Society for Composites. In Proceedings of the 26th Annual Technical Conference: Canadian Association for Composite Structures and Materials, Montreal, QC, Canada, 26–28 September 2011. [Google Scholar]
Figure 1. Geometrical description of the specimens.
Figure 1. Geometrical description of the specimens.
Polymers 12 00554 g001
Figure 2. Resistance curves and mechanical material properties.
Figure 2. Resistance curves and mechanical material properties.
Polymers 12 00554 g002
Figure 3. Specimen stacking sequence and delamination position along the thickness.
Figure 3. Specimen stacking sequence and delamination position along the thickness.
Polymers 12 00554 g003
Figure 4. Strain gauges position in (a) TOUGH#1, TOUGH#2, MFB#1 and MFB#2 specimens; (b) TOUGH#3 and MFB#3 specimens.
Figure 4. Strain gauges position in (a) TOUGH#1, TOUGH#2, MFB#1 and MFB#2 specimens; (b) TOUGH#3 and MFB#3 specimens.
Polymers 12 00554 g004
Figure 5. Samples preparation for strain monitoring (a) and DIC (b).
Figure 5. Samples preparation for strain monitoring (a) and DIC (b).
Polymers 12 00554 g005
Figure 6. Ultrasonic inspection.
Figure 6. Ultrasonic inspection.
Polymers 12 00554 g006
Figure 7. Test rig.
Figure 7. Test rig.
Polymers 12 00554 g007
Figure 8. Speckle pattern in Digital Image Correlation (DIC).
Figure 8. Speckle pattern in Digital Image Correlation (DIC).
Polymers 12 00554 g008
Figure 9. Digital Image Correlation rig.
Figure 9. Digital Image Correlation rig.
Polymers 12 00554 g009
Figure 10. Description of the Fibre Bridging (FB) modulus operation based on the R-curve.
Figure 10. Description of the Fibre Bridging (FB) modulus operation based on the R-curve.
Polymers 12 00554 g010
Figure 11. Finite elements model.
Figure 11. Finite elements model.
Polymers 12 00554 g011
Figure 12. Finite elements model: contact elements.
Figure 12. Finite elements model: contact elements.
Polymers 12 00554 g012
Figure 13. Toughened specimen deformed shape at (a) delamination onset, (b,c) intermediate delamination state and (d) complete delamination.
Figure 13. Toughened specimen deformed shape at (a) delamination onset, (b,c) intermediate delamination state and (d) complete delamination.
Polymers 12 00554 g013
Figure 14. Specimen sensitive to fibre bridging deformed shapes at (a) delamination onset, (b,c) intermediate delamination state and (d) complete delamination.
Figure 14. Specimen sensitive to fibre bridging deformed shapes at (a) delamination onset, (b,c) intermediate delamination state and (d) complete delamination.
Polymers 12 00554 g014
Figure 15. Compressive load vs applied displacements.
Figure 15. Compressive load vs applied displacements.
Polymers 12 00554 g015
Figure 16. Compressive load as a function of the measured and computed strains in the loading direction: toughened material (a) and material sensitive to fibre bridging (b).
Figure 16. Compressive load as a function of the measured and computed strains in the loading direction: toughened material (a) and material sensitive to fibre bridging (b).
Polymers 12 00554 g016
Figure 17. Compressive load as a function of the DIC measured strains: comparison with the numerical results.
Figure 17. Compressive load as a function of the DIC measured strains: comparison with the numerical results.
Polymers 12 00554 g017
Figure 18. Compressive load as a function of the DIC measured out-of-plane displacements: comparison between experimental data and numerical results.
Figure 18. Compressive load as a function of the DIC measured out-of-plane displacements: comparison between experimental data and numerical results.
Polymers 12 00554 g018
Figure 19. Delaminated area as a function of the compressive load: comparison between experimental data and numerical results.
Figure 19. Delaminated area as a function of the compressive load: comparison between experimental data and numerical results.
Polymers 12 00554 g019
Figure 20. DIC out-of-plane displacement contour plot (MFB#3 sample) at (a) delamination buckling onset, (b) delamination growth initiation, (c) global buckling and (d) snap-trough phenomenon.
Figure 20. DIC out-of-plane displacement contour plot (MFB#3 sample) at (a) delamination buckling onset, (b) delamination growth initiation, (c) global buckling and (d) snap-trough phenomenon.
Polymers 12 00554 g020
Figure 21. Numerical out-of-plane displacement contour plot for toughened material (units in mm) at (a) 0.3 mm applied displacement, (b) 0.4 mm applied displacement, (c) 0.55625 mm applied displacement, (d) 0.73125 mm applied displacement, and (e) 1.18 mm applied displacement.
Figure 21. Numerical out-of-plane displacement contour plot for toughened material (units in mm) at (a) 0.3 mm applied displacement, (b) 0.4 mm applied displacement, (c) 0.55625 mm applied displacement, (d) 0.73125 mm applied displacement, and (e) 1.18 mm applied displacement.
Polymers 12 00554 g021
Figure 22. Numerical out-of-plane displacement contour plot for material sensitive to fibre bridging (units in mm) at (a) 0.3 mm applied displacement, (b) 0.4 mm applied displacement, (c) 0.55625 mm applied displacement, and (d) 0.73125 mm applied displacement.
Figure 22. Numerical out-of-plane displacement contour plot for material sensitive to fibre bridging (units in mm) at (a) 0.3 mm applied displacement, (b) 0.4 mm applied displacement, (c) 0.55625 mm applied displacement, and (d) 0.73125 mm applied displacement.
Polymers 12 00554 g022
Figure 23. Numerical and experimental deformed shape comparison. Material sensitive to fibre bridging.
Figure 23. Numerical and experimental deformed shape comparison. Material sensitive to fibre bridging.
Polymers 12 00554 g023
Table 1. Tests matrix.
Table 1. Tests matrix.
SG1bisSG1SG2DIC
90°90°
TOUGH#1
TOUGH#2
TOUGH#3
MFB#1
MFB#2
MFB#3
Table 2. Summary of experimental results.
Table 2. Summary of experimental results.
Global Buckling Delamination Buckling Snap-Through Buckling
Load [kN]Strain [με]Load [kN]Strain [με]Load [kN]Strain [με]
TOUGH#144.37−178039.4−149040.4−9960
TOUGH#245.2−179039.1−148040.3−8700
TOUGH#345.2−178539.3−148340.3−8790
MFB#141.2−160039.5−147039.3−6830
MFB#240.9−160039.1−147539.0−6750
MFB#340.1−159039.2−147139.4−7000

Share and Cite

MDPI and ACS Style

Riccio, A.; Russo, A.; Sellitto, A.; Toscano, C.; Alfano, D.; Zarrelli, M. Experimental and Numerical Assessment of Fibre Bridging Toughening Effects on the Compressive Behaviour of Delaminated Composite Plates. Polymers 2020, 12, 554. https://0-doi-org.brum.beds.ac.uk/10.3390/polym12030554

AMA Style

Riccio A, Russo A, Sellitto A, Toscano C, Alfano D, Zarrelli M. Experimental and Numerical Assessment of Fibre Bridging Toughening Effects on the Compressive Behaviour of Delaminated Composite Plates. Polymers. 2020; 12(3):554. https://0-doi-org.brum.beds.ac.uk/10.3390/polym12030554

Chicago/Turabian Style

Riccio, Aniello, Angela Russo, Andrea Sellitto, Cinzia Toscano, Davide Alfano, and Mauro Zarrelli. 2020. "Experimental and Numerical Assessment of Fibre Bridging Toughening Effects on the Compressive Behaviour of Delaminated Composite Plates" Polymers 12, no. 3: 554. https://0-doi-org.brum.beds.ac.uk/10.3390/polym12030554

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