Next Article in Journal
A Bridge-Shaped Vibration Energy Harvester with Resonance Frequency Tunability under DC Bias Electric Field
Next Article in Special Issue
Internet of Things-Based ECG and Vitals Healthcare Monitoring System
Previous Article in Journal
Cuffless Blood Pressure Monitoring: Academic Insights and Perspectives Analysis
Previous Article in Special Issue
Simulation and Test of a MEMS Arming Device for a Fuze
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Micro-Mechanical Investigation of Interfacial Debonding in Carbon Fiber-Reinforced Composites Using Extended Finite Element Method (XFEM) Approach

1
Department of Aeronautics & Astronautics, Institute of Space Technology, Islamabad 44000, Pakistan
2
School of Mechanical & Manufacturing Engineering (SMME), National University of Sciences and Technology (NUST), Islamabad 44000, Pakistan
3
Automotive Design and Crashworthiness Research, Islamabad 44000, Pakistan
4
Aeronautical Engineering Department, Military Technological College, Muscat P.O. Box 262, Oman
*
Author to whom correspondence should be addressed.
Micromachines 2022, 13(8), 1226; https://doi.org/10.3390/mi13081226
Submission received: 4 July 2022 / Revised: 25 July 2022 / Accepted: 26 July 2022 / Published: 30 July 2022
(This article belongs to the Special Issue Micro-Scale Energy Harvesting Devices)

Abstract

:
The interface debonding in carbon fiber-reinforced polymers is analyzed and evaluated using the extended finite element method (XFEM). In order to accurately evaluate the bonding properties between fibers and matrix, different tests were carried out, including the multiple tests for different orientations to study longitudinal, transversal, and shear properties of unidirectional carbon fiber-reinforced composites. Extensive experimentation has been performed in all the different groups and categories with different dimensions and parameters in order to ascertain the values of strength and the prediction of the damage to the structure. The experimental and numerical comparison provided significant trends and data to evaluate the mechanical properties of the interface. The values of stiffness and strength are compared and validated. Development of Representative Volume Element (RVE) for progressive damage model to these damage phenomena has already been performed as a feasibility study for the model, though it is not included in this particular paper. The results of this research for all the experimental and numerical sets can serve as reliable data in the microsimulation of devices and sensitive parameters that include carbon fiber-reinforced light metal matrix composites and makes a better investigative model that contributes to various conditions. It further offers an investigation of the microscopic deformation mechanisms in the composites.

1. Introduction

The emergence of a new class of fiber composite materials for use in aircraft has its origin in a number of technical and scientific developments, the starting point being the discovery last century of synthetic organic materials derived from vegetation and coal. The development of carbon fibers of high strength and stiffness creates special problems because the graphite crystal from which they are produced is extremely anisotropic. The crystals have to be aligned with the fiber axis to obtain high performance. Fiber composites are making an impact across the whole field of structural materials and their use has been growing at about 7% per annum for the last decade and a half, regardless of the economic downturns. It is, however, the requirements of high-performance aircraft for defense and civil usage, aerospace vehicles, and rockets that have fueled the growth in the development of fiber composites. The straightforward requirement for weight loss favors their use [1,2,3].
The fiber–matrix interface is a critically thin surface between fiber and matrix. A key role is placed by the interface in determining the mechanical properties of composite materials because the load applied externally has to be transmitted from matrix to fibers through the interface [4,5,6]. The bonding quality between the fibers and matrix is responsible for the strength of the fiber–matrix interface. Bonding at the fiber–matrix interface is particular for every system. A complex chemical bond between fiber and matrix is formed in the case of the carbon/epoxy composites. Composites having a strong interface have high strength and stiffness as compared to composites having a weak interface. Composites having a strong interface behave in a brittle manner with 2–3% elastic elongation-to-failure with no sufficient plastic deformation. Therefore, in some special cases, a relatively weak interface is required. The most commonly observed failure in composites is delamination. It is the separation of layers due to the weakness of the interface between laminates. Delamination may originate from low-velocity impact, strangeness in structural load paths that cause out-of-plane loads, or from heterogeneous and stacked nature, which create local out-of-plane load [7,8,9].
The development of Representative Volume Element RVE [10,11,12] will provide us with more strength and make our research even more interesting as our domain of research is wider now. In any research, it is very important to know what the previous theories, hypotheses, or findings of researchers and scientists are. When we can understand such previous work and research, we can then take it forward to make it even better. The main portion of usage of carbon fiber-reinforced composites is in the aerospace industry including military operations, space operations, and civil aircraft, where weight reduction is critical for higher speeds and increased payload [13,14,15]. In aerospace, the use of fiber-reinforced polymers has experienced steady growth. The structural integrity and durability of these early components have built up confidence in their performance and promoted developments of other structural aircraft components, resulting in an increasing number of composites being used in military aircraft. For example, the airframe of an AV-8B contains 25% by weight of carbon fiber-reinforced epoxy, while the F22 contains 25% by weight of carbon fiber-reinforced polymers. The stealth characteristics of these aircraft are because of the usage of carbon fiber-reinforced composites. In commercial aircraft, these composites began to be used along with a few secondary components. The Boeing 777 contains carbon fiber-reinforced epoxy 10% by weight which was introduced in 1995. On the other hand, the Boeing 787 Dreamliner, which was introduced in 2009, contains carbon fiber-reinforced polymers of almost 50% by weight [16].
Damage progression includes three main technology areas for obtaining a solution: instrumentation and data processing hardware, data integration, and predictive modeling. It is believed that damage progression is, in most cases, an emerging technology. The general modeling issues that we consider are service environment, considerations in material specification, and material microstructure considerations (including single phase materials and multiphase materials). First, we start with forms of damage that can occur in our selected type of material, which can be found in many forms within the same material. The generic working model we used here is any structural deviation from a perfect material state, such as material separations, imperfections, and displacement discontinuities. Damage occurs virtually at all length scales which results in non-linear, irreversible dissipative deformation mechanisms manifesting as history-dependent behavior. Usually, damage occurs at lower length scales, evolves and interacts with these scales, and then manifests with higher scales. As damages occur at the microscale, modeling relative volume elements to begin with will provide us with more knowledge on how that damage propagates in any material. Thereafter, by characterizing that material’s behavior it will then be useful to consider a monotonic stress–strain curve which shows material behavior from initial loading up to material failure. After generating stress–strain curves, we can predict the total life or the remaining life of the material [4].

2. Material Models

Material properties of chosen fiber-reinforced composites are given in Table 1, Table 2 and Table 3 which are standard properties. The homogenized values are also calculated for the interfacial debonding study.
Here Em represents elastic modulus and υ as poisson ratios whereas XT and Xc are tensile and compressive strengths, respectively.

3. Experimental Study

3.1. Experimental Standards and Specimen

Carbon Fiber-Reinforced Composites (CFRCs) are modeled for the conditions of the multiscale XFEM framework, followed by the development of an enrichment scheme of equations. The proposed micro-mechanical model for the three-point bending test is validated by comparison with experimental data. The study provided a significant benchmark for the prediction of the XFEM progressive damage [14,18,19] and a framework for the mechanical behavior of carbon fiber-reinforced composites. A three-point bending test is performed to determine the bending stress, flexural stress, and the flexural strain in the composite materials. The three-point bending test is carried out in the same machine described previously, the test samples are prepared according to ASTM standards. The dimension design of transverse fiber tow tensile specimens refers to the standard sample size recommended by ASTM D638, as shown in Figure 1 processed through Northwestern Polytechnical University, China. The die is cut along the symmetrical middle line of the resin casting die for laying the fiber tows. In order to facilitate demolding, the following setup is used as a die mechanism in this experiment.
The experiment adopts a WDW-100 micro-controlled electronic universal testing machine as shown in Figure 2, processed through Northwestern Polytechnical University, China. The engineering constants of epoxy and carbon fiber are already discussed. The accuracy of load reading should be higher than 1% of the measured value. The fixture system can ensure that the center line of the sample is consistent with the center axis of the testing machine. The tension tester is equipped with an extensometer, which is connected to a continuous recording device to automatically record the elongation of the specimen in the clamp. The testing process refers to ASTM 3039. processed through Northwestern Polytechnical University, China.
Understanding the complicated damage mechanism of composites is still challenging. With the widespread use of CFRP, research into the damage mechanisms of CFRP structures to CFRP development in engineering structures is of major concern. Latest improvements focus on the interlaminar interface in modeling delamination propagation. Multi-layered delamination [8,20] is usually the result of damage in composite laminates where it grows up and transforms into numerous interlaminar interfaces.

3.2. Experimental Investigation of Interfacial Bonding and Mechanical Properties

To accurately evaluate the bonding properties between fibers and matrix, the tensile tests were carried out, including the tensile tests of transverse fiber orientation and the shear tests of 45-degree fiber orientations. Figure 3 shows the samples ready for the test and Figure 4a,b represents the scheme for an interfacial crack within the domains Ω1 and Ω2 for the corresponding radii re.
The mechanical properties of unidirectional lamina and adhesive are presented in Table 4. Considering the symmetric and boundary conditions of the panel, within the model.

3.3. Testing of Transverse Bond Strength at Interface for CFRCs

In the middle of the resin stretching spline die, a slot is opened, the fiber tow is laid perpendicular to the direction of the sample stretching, and the transverse tow stretching sample prepared by co-curing after pouring resin is stretched to fracture under uniform loading of the appropriate mechanical device [21]. Tensile strength, modulus of elasticity, and strain at maximum load are calculated according to the tensile stress–strain curve. Previous studies have shown that the tensile strength of transverse tow tensile specimens can reflect the strength of the fiber/resin–interface bonding. The transversal damage has been particularly studied in one of the recent studies of the author. The tested specimens are shown in Figure 5.

3.4. Interfacial Shear Strength Test of CFRPs

The interfacial shear tow tension specimen prepared by co-curing after pouring resin is stretched to fracture under uniform loading of the appropriate mechanical device. According to the tensile stress–strain curve, the interfacial shear strength and shear elastic modulus are calculated. Dimension design of tensile specimens of oblique 450 fiber arrangements refers to the standard sample size recommended by ASTM D638 [22] processed through Northwestern Polytechnical University, China. Along the symmetrical center line of the resin casting die, the die is cut diagonally 450 to lay the fiber arrangements. The other methods of operation are the same as those of transverse fiber tow tensile specimens. The drawing sample die is shown in Figure 6. The final curing and polishing of qualified resin matrix splines are shown in Figure 7.

4. Numerical Model Estimation and Representation

The dimensional relations among its different components were still undetermined. We have one full cylindrical fiber at the center and four quarter cylindrical fibers at all corners for the sake of periodicity. So, we have a total of two carbon fibers other than carbon. Now we can say that the total volume occupied by the carbon fibers is equal to the volume of the unit cell multiplied by the fiber volume fraction which is 0.6. This is shown in Equation (1), where V f and V c represent the volume of carbon fibers and the cubical unit cell. Equation (2) simplifies the expression of the volume of fiber and the unit cell. Where a represents the length of the sides of the unit cell and r f represents the radius of the carbon fiber. Here, VFF is the volume filler fraction.
2 V f = V F F   V c
2 π × r f 2 × a = V F F   a 3
The relationship is presented in Equation (3), which gives us the length of the unit cell.
a = 2 π × r f 2 V F F
Using the same concept as the above equations, we derived the expressions given in Equations (4) and (5). The volume fraction for the carbon fibers was selected to be 0.5% of the total matrix volume which is less value compared to the carbon fibers and it will still give satisfactory results for stress-bearing in the normal direction to the carbon fibers.
V n f = 0.005 × V m
V s n f = π ( r n f ) 2 × a
N n = V n f / V s n f
where V n f and V s n f are the total volume of Nano-fibers and the volume of a single Nano-fiber, respectively. We can calculate the number of Nano-fibers in our unit cell by using the expression given in Equation (6). The r n f represents the radius we selected for our Nano-fibers. Using all of the above expressions, we calculated different dimensions of our representative volume element, which are given in Table 5. We calculated the number of Nano-fibers to be 14 and, as already discussed, we have two carbon fibers. We selected our fiber diameter to be 9 µm which is the most common diameter of the Carbon T-300 fiber, processed through Northwestern Polytechnical University, China being used these days.
The Extended Finite Element Method (XFEM) is based on the Finite Element Method (FEM) used to treat discontinuities. The objective of the research is to evaluate the bonding properties between fiber and matrix and to define a technique for comparison of mechanical properties i.e., strength and stiffness in interfacial debonding. This will help us improve the interfacial properties of composite laminates in order to avoid delamination. Implementation of XFEM [23,24,25] and simulations of different loading conditions for the proposed ECDM model and the modified CZM are tools for the prediction of the crack growth due to delamination. The simplification of the modeling of discontinuous phenomena is the main benefit of XFEM methods for different problems in materials science. In the traditional formulation of the FEM, a fracture is modeled by requiring the fracture to follow element edges. In contrast, the fracture geometry in the X-FEM does not need to be aligned with the element edges, which is a great flexibility.
In damage analysis, there are, in general, two different types of failure criteria for CFRCs material. Strength of Material Criteria based on stresses for the damage initiation and Fracture Mechanics Criteria based on energy for the damage growth. Hashin and Rotem (1973), Hashin (1980), Matzenmiller et al. (1995), and Camanho and Davila (2002) all contributed to the Abaqus anisotropic damage model. Fiber rupture in stress, fiber buckling and kinking in compression, matrix cracking under transverse tension and shearing, and matrix crushing under transverse compression and shearing are all taken into account. In Abaqus, the initiation criteria proposed by Hashin and Rotem (1973) and Hashin (1980), in which the failure surface is expressed in the effective stress space, determine the onset of damage. The response of the material is computed from Equation (7).
σ = C d ε
In Equation (7) ε is the strain and Cd is the elasticity matrix, which reflects any damage and has the form shown in Equation (8).
C d = 1 D 1 d f E 11 1 d f 1 d m v 21 E 11 0 1 d f 1 d m v 12 E 22 1 d m E 22 0 0 0 1 d s G 12 D
In Equation (8) D = 1 1 d f 1 d m v 12 v 21 , d f is the current state of fiber damage, d m is the current state of matrix damage, d s is the current state of shear damage, E 11 is the Young’s modulus in fiber direction (longitudinal direction), E 22 is the Young’s modulus in the direction perpendicular to the fibers (transverse direction), G12 is the shear Modulus and v 12 , and v 21 are Poisson ratios. The longitudinal and transverse direction of uni-directional lamina is shown in Figure 8 as direction 1 and 2, respectively.
The damage variables d f , d m , and d s were derived, from the damage variables d f t , d f c , d m   t , and d m c shown in Equations (9)–(13) which were referred to in the four failure modes discussed previously.
d f = { d f   t   i f   σ 11 0
d f = {   d f   c   i f   σ 11 < 0
d m = {   d m   t   i f   σ 22 0
d m = {   d m   c   i f   σ 22 < 0
d s = 1 1 d f   t 1 d f   c 1 d m   t 1 d m   c
The damage variables d f   t ,   d f   c ,   d m   t , and d m   c in Equations (9)–(13) are the internal damage variables in the fiber and matrix phases of the lamina, under tension or compression loadings. Damage initiation refers to the beginning of a degradation. For CFRCs, it is based on Hashin’s failure criteria, using four different damage mechanisms that is fiber tension, fiber compression, matrix tension, and matrix compression as shown in Equations (14)–(17).
Fiber Tension: σ 11 0
F f T = σ 11 X 1 T 2 + α τ 12 S 12 2 = 1
Fiber Compression: σ 11 < 0
F f C = σ 11 X 1 C 2 = 1
Matrix Tension: σ 22 0
F m T = σ 22 X 2 T 2 + α τ 12 S 12 2 = 1
Matrix Compression: σ 22 < 0
F m C = σ 22 2 S 13 2 + X 2 C 2 S 13 2 1 σ 22 X 2 C + τ 12 S 12 2 = 1
In Equations (14)–(17) σ 11 ,   σ 22 , and τ 12 are the longitudinal, transverse, and shear stresses in the lamina, X 1 T and X 1 C refers to the tensile and compression strength in the fiber direction (longitudinal tensile and compressive strength), and X 2 T   and X 2 C   refers to the tensile and compression strength in the transverse direction. S 12 and S 13 are the longitudinal and transverse shear strength, respectively. The coefficient α determines the contribution of shear stress to the fiber tensile damage initiation in the present work. The material was linearly elastic before damage initiation, based on the brittle behavior of the CFRC. Damage evolution took place after the damage was initiated. After one or more damage initiation conditions are met, the damage evolution description describes how the material degrades. Multiple damage evolution forms, one for each given damage initiation criterion, can act on a material at the same time. Damage variables that have values ranging from zero (undamaged state) to one (damaged state) influence the reduction of stiffness coefficients. For the Hashin’s damage evolution model, the data table contains the fields shown in Table 6 and a brief comparison in Table 7 with statistical segregation. You define damage stabilization for fiber-reinforced materials by entering viscosity coefficients for each of the potential failure modes. Each of the viscous coefficients should be small compared to the increment size. Viscous regularization is intended to improve convergence as the material fails.
Whilst the matrix crack propagated across the fiber–matrix interface, it also deflects alongside the interface. Fracture mechanics and shear strength methods are the best ways to deal with the interface debond problem. The fracture mechanics method deals the fiber–matrix interface debonding as a crack propagation problem, in which the interface debonding [8,26] happens because the strain energy release rate at interface reaches debonding toughness. Moreover, the shear strength method is based on maximum shear stress criteria. Sun and Singh analyzed matrix multi-cracking and fiber/matrix interfacial debonding. However, for an extended interface debond period, the fracture mechanics approach provided a good fit. The fracture mechanics approach is used to determine the fiber/matrix–interface debonding, and is given as: Applied the XFEM step enrichment to model narrow damage localization zones as given in Equation (18) for interface debond energy and Equation (19) for the interface debond length. Table 6 shows the numerical evaluation and parameter of selection for failure in different scenarios.
ζ d = F 4 π r f w f 0 L d 1 2 0 L d τ i v x L d d x
where ζd is interface debond energy, F is fiber load at matrix crack plane, v(x) is relative displacement between fiber and matrix, and w f (0) shows the fiber axial displacement at the matrix crack plane. The interface debond length L d is determined by equation.
L d = r f 2 V m E m σ V f E c τ i 1 ρ   r f 2 ρ 2 + r f V m E m E f τ i 2 E c ζ d
Table 8 shows the predicted values for the model and comparison with the literature and previous studies with the shear dominated failure criterion in Table 6 and Table 7, for all the variables and properties previously discussed. The result’s predicted values lie in a very good agreement for longitudinal strength prediction. It clearly represents the validation for the transversal and longitudinal principle for the unidirectional fibers. Where X T are tensile, and X C are compressive strengths for the set.
A 3D model is developed with geometry dimensions similar to experimental specimens which is shown in Table 4, Table 5 and Table 6. In comparison to natural specimen geometry, the geometry of the numerical model is made to appear as authentic as possible. A three-dimensional, deformable shell planar function is used to build the model. Figure 9 shows that three specimen geometry with different orientation. 7a shows the zero-degree fiber orientation and 14 plies, 9b shows the 90-degree fiber orientation with 14 plies, and 9c shows 5 plies with 0/90/0/90/0 plies orientation.
The research started with tensile testing to study the transversal damage for category A, B and C [27]. Specimen A and B are composite layups of 0° and 90°, respectively, whereas Specimen C is a cross-ply composite layup of 5 plies. One ply usually consists of two constituents, fiber and matrix, which can both be damaged individually. A laminate is made by piling multiple plies together in different orientations. The lay-up is the arrangement of the laminate that shows its ply composition with various fiber orientations. The detailed model and the specimen are shown in Table 8, Table 9 and Table 10 and the representation is shown in Figure 9a–c.

5. Major Critical Results and Discussions

Furthermore, it is followed by the comparison of stress component to the fiber (vertical) direction for reference models on beginning the crack propagation for a similar width and the simulations are shown below in Figure 10 for XFEM damage comparison for 3P/ 0˚ specimen which depicts a very reliable model in Figure 11 and Figure 12. The experimental method is the same as that of transverse fiber tow tension. XFEM damage specimen. The same is crystal clear from Table 11, Table 12 and Table 13 for the comparison of models for the stresses and overall strength, which is promising.
The experimental results in Figure 10 and Figure 11 display the damage and crack propagation as comparison and related load–displacement curve and contour plots for every carbon fiber composite ply as the crack and delamination. It is clearly proven that the transversal damage for the group A of 45 has crack initiation at the center of the specimen, which can be evaluated from the other cases. After the failure of the specimen, from load and displacement data points and analyzes the results.
The curves clearly illustrate the anisotropic brittle behavior of carbon-reinforced composites. When a load is applied in fiber direction (longitudinal direction) Figure 10 clearly illustrates that aligned continuous fibers have a much higher tensile strength.
It is possible to obtain the mechanical properties of composite material after conducting tensile tests and processing the results, as shown in Table 14. It is important that the material has high mechanical strength and strain values, and that the elastic modulus is in the right range, a representation for a sorted set of experiment is shown in Figure 12.
Different fracture modes, such as brittle matrix fracture and fiber splitting, were also observed in the test results as shown in Figure 13.
The initiated crack propagates in the upper surface that causes the initiation of damage to the adjoining layers and comparative studies have been considered as [28], the result lies in a reasonable agreement. The values had been calculated as a comparison of the predicted values and the results from the experimental outputs had been numerically calculated as given in Table 14.
The boundary conditions and the displacement load for 90-degree orientation have been analyzed and lies within the error range of less than 3%. It is observed that in the case of transverse tension or compression, the phenomenon of transverse bearing capacity of the matrix loses when the overall bearing capacity of the matrix is damaged.
At the end of the test, if the fiber bundles are intact and only the resin at the standard interval necks or breaks, which belongs to the invalid mode, the test results will be discarded. Table 15 shows the detailed results.

6. Conclusions

This experimental work is in agreement with numerical simulations which provides deep understanding of crack behavior in carbon fiber-reinforced composites, particularly for the interfacial micro-debonding. The extensive experimentation and numerical analysis has been carried out for the development of a flexible framework for the prediction of the damage computational model, techniques for fracture and delamination behaviour of the crack in carbon fiber-reinforced composites, particularly for the different modes of loads, and boundary conditions have been evaluated using XFEM. The experimentation and numerical response of the framework is in a good agreement, hence it conforms to the proposed technique as a simplified and efficient tool for beginners and professionals. This comparative study has also been carried out with previous numerically studied sets and the results lie in reasonable convergence. The results obtained with this technique are in agreement with experimental results and the error percentage is quite reasonable. This research has a unique methodology by considering different phenomenon at the same time for interfacial debonding strength. It is recommended to carry out further quantitative work in order to validate the numerical algorithms. The stress–strain behavior gives the reader a better perspective in order to investigate various affects on the composite. This study is concise and converged with the previously performed models and research for all sets of results which were grouped into multiple categories. This study has been carried out for both transversal damage and the phenomenon of delamination in composites with the multiple modes and categories A, B and C.

Author Contributions

The authors confirm contribution to the paper as follows: study conception and design: R.F.S., S.R.Q. and A.A.K.; data collection: M.U.S., O.M.Q. and A.M.; analysis and interpretation of results: R.F.S., S.R.Q., M.U.S. and A.M.; draft manuscript preparation: R.F.S., S.R.Q., A.A.K. and O.M.Q. Supervision and structure of writing: A.A.K. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Data Availability Statement

No repositories involved or shared.

Acknowledgments

I would like to acknowledge my research group specially AsfanYar Amjad, Zahid, Mahnoor, Saood, and Zeeshan for giving the article a significant share. I would like to acknowledge Hossein Raza Hamdani for providing us a smooth atmosphere of research and arrangements whenever asked. Would like to thank every member of Institute of Space Technology and Northwestern Polytechnical University for helping us out in experimentation.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Hulton, A.W.; Cavallaro, P.V. Comparing Computational and Experimental Failure of Composites Using XFEM. Am. Soc. Mech. Eng. 2016, 50633, V009T12A004. [Google Scholar] [CrossRef]
  2. Bai, R.; Bao, S.; Lei, Z.; Yan, C.; Han, X. Finite element inversion method for interfacial stress analysis of composite single-lap adhesively bonded joint based on full-field deformation. Int. J. Adhes. Adhes. 2018, 81, 48–55. [Google Scholar] [CrossRef]
  3. Swati, R.F.; Wen, L.H.; Elahi, H.; Khan, A.A.; Shad, S. Extended finite element method (XFEM) analysis of fiber reinforced composites for prediction of micro-crack propagation and delaminations in progressive damage: A review. Microsyst. Technol. 2018, 25, 747–763. [Google Scholar] [CrossRef]
  4. García, I.; Paggi, M.; Mantič, V. Fiber-size effects on the onset of fiber–matrix debonding under transverse tension: A comparison between cohesive zone and finite fracture mechanics models. Eng. Fract. Mech. 2014, 115, 96–110. [Google Scholar] [CrossRef]
  5. Gilabert, F.; Garoz, D.; Van Paepegem, W. Macro- and micro-modeling of crack propagation in encapsulation-based self-healing materials: Application of XFEM and cohesive surface techniques. Mater. Des. 2017, 130, 459–478. [Google Scholar] [CrossRef] [Green Version]
  6. Longbiao, L. Micromechanical Modeling for Tensile Behaviour of Carbon Fiber—Reinforced Ceramic—Matrix Composites. Appl. Compos. Mater. 2015, 29, 773–790. [Google Scholar] [CrossRef]
  7. Yazdani, S.; Rust, W.J.; Wriggers, P. Delamination onset and growth in composite shells. Comput. Struct. 2018, 195, 1–15. [Google Scholar] [CrossRef]
  8. Higuchi, R.; Okabe, T.; Nagashima, T. Numerical simulation of progressive damage and failure in composite laminates using XFEM/CZM coupled approach. Compos. Part A Appl. Sci. Manuf. 2017, 95, 197–207. [Google Scholar] [CrossRef]
  9. Zhao, L.; Zhi, J.; Zhang, J.; Liu, Z.; Hu, N. XFEM simulation of delamination in composite laminates. Compos. Part A Appl. Sci. Manuf. 2016, 80, 61–71. [Google Scholar] [CrossRef]
  10. El Moumen, A.; Kanit, T.; Imad, A. Numerical evaluation of the representative volume element for random composites. Eur. J. Mech. A/Solids 2021, 86, 104181. [Google Scholar] [CrossRef]
  11. Henrich, M.; Pütz, F.; Münstermann, S. A Novel Approach to Discrete Representative Volume Element Automation and Generation-DRAGen. Materials 2020, 13, 1887. [Google Scholar] [CrossRef] [Green Version]
  12. Bargmann, S.; Klusemann, B.; Markmann, J.; Schnabel, J.E.; Schneider, K.; Soyarslan, C.; Wilmers, J. Generation of 3D representative volume elements for heterogeneous materials: A review. Prog. Mater. Sci. 2018, 96, 322–384. [Google Scholar] [CrossRef]
  13. Sanjay, M.R.; Arpitha, G.R.; Yogesha, B. Study on Mechanical Properties of Natural-Glass Fibre Reinforced Polymer Hybrid Composites: A Review. Mater. Today Proc. 2015, 2, 2959–2967. [Google Scholar] [CrossRef]
  14. Swati, R.F.; Wen, L.H.; Elahi, H.; Khan, A.A.; Shad, S. Experimental and numerical investigation of transversal damage in carbon fiber reinforced composites using X-FEM analysis. J. Mech. Sci. Technol. 2018, 33, 205–211. [Google Scholar] [CrossRef]
  15. Tan, H.; Liu, C.; Huang, Y.; Geubelle, P. The cohesive law for the particle/matrix interfaces in high explosives. J. Mech. Phys. Solids 2005, 53, 1892–1917. [Google Scholar] [CrossRef]
  16. Wu, C.-P.; Liu, Y.-C. A review of semi-analytical numerical methods for laminated composite and multilayered functionally graded elastic/piezoelectric plates and shells. Compos. Struct. 2016, 147, 1–15. [Google Scholar] [CrossRef]
  17. Swati, R.F.; Hua, W.L.; Elahi, H.; Khan, A.A. Extended Finite Element Method Damage Analysis of Carbon Fiber Reinforced Composites Using Multiscale Method and Study of Crack Propagation in Mixed Mode using ABAQUS. In Proceedings of the 4th International Conference on Mechanics and Mechatronics Research, Xian, China, 22–24 May 2021. [Google Scholar]
  18. Chen, J. Predicting Progressive Delamination of Stiffened Fibre-Composite Panel and Repaired Sandwich Panel by Decohesion Models Application of decohesive model with mixed damage scale in fracture. J. Thermoplast. Compos. Mater. 2002, 15, 429–442. [Google Scholar] [CrossRef]
  19. Iqbal, N.; Ali, Y.; Haq, I.U.; Lee, S. Progressive interface debonding in composite electrodes of Li-ion batteries via mixed-mode cohesive zone model: Effects of binder characteristics. Compos. Struct. 2021, 259, 113173. [Google Scholar] [CrossRef]
  20. Wang, H.; Zhou, H.; Ji, H.; Zhang, X. Application of extended finite element method in damage progress simulation of fiber reinforced composites. Mater. Des. 2014, 55, 191–196. [Google Scholar] [CrossRef]
  21. Matveeva, A.Y.; Lomov, S.V.; Gorbatikh, L. Debonding at the fiber/matrix interface in carbon nanotube reinforced composites: Modelling investigation. Comput. Mater. Sci. 2019, 159, 412–419. [Google Scholar] [CrossRef]
  22. Almeida, J.H.S.; Souza, S.D.B.; Botelho, E.C.; Amico, S.C. Carbon fiber-reinforced epoxy filament-wound composite laminates exposed to hygrothermal conditioning. J. Mater. Sci. 2016, 51, 4697–4708. [Google Scholar] [CrossRef] [Green Version]
  23. Petrov, N.A.; Gorbatikh, L.; Lomov, S.V. A parametric study assessing performance of eXtended Finite Element Method in application to the cracking process in cross-ply composite laminates. Compos. Struct. 2018, 187, 489–497. [Google Scholar] [CrossRef]
  24. Swati, R.F.; Elahi, H.; Wen, L.H.; Khan, A.A.; Shad, S.; Mughal, M.R. Investigation of tensile and in-plane shear properties of carbon fiber reinforced composites with and without piezoelectric patches for micro-crack propagation using extended finite element method. Microsyst. Technol. 2018, 25, 2361–2370. [Google Scholar] [CrossRef] [Green Version]
  25. van der Meer, F.P.; Sluys, B.; Hallett, S.R.; Wisnom, M.R. Computational modeling of complex failure mechanisms in laminates. J. Compos. Mater. 2012, 46, 603–623. [Google Scholar] [CrossRef]
  26. Sun, Q.; Zhou, G.; Meng, Z.; Guo, H.; Chen, Z.; Liu, H.; Kang, H.; Keten, S.; Su, X. Failure criteria of unidirectional carbon fiber reinforced polymer composites informed by a computational micromechanics model. Compos. Sci. Technol. 2019, 172, 81–95. [Google Scholar] [CrossRef]
  27. Chen, J. An extended cohesive damage model with a length scale in fracture analysis of adhesively bonded joints. Eng. Fract. Mech. 2014, 119, 202–213. [Google Scholar] [CrossRef] [Green Version]
  28. Li, X.; Chen, J. An extended cohesive damage model for simulating multicrack propagation in fibre composites. Compos. Struct. 2016, 143, 1–8. [Google Scholar] [CrossRef] [Green Version]
Figure 1. Tensile sample die.
Figure 1. Tensile sample die.
Micromachines 13 01226 g001
Figure 2. Tensile sample die.
Figure 2. Tensile sample die.
Micromachines 13 01226 g002
Figure 3. Samples.
Figure 3. Samples.
Micromachines 13 01226 g003
Figure 4. Tension specimen for experimentation (a) domain and (b) estimation for radii analysis. A half of FE mesh was used.
Figure 4. Tension specimen for experimentation (a) domain and (b) estimation for radii analysis. A half of FE mesh was used.
Micromachines 13 01226 g004
Figure 5. Test results of the different specimen.
Figure 5. Test results of the different specimen.
Micromachines 13 01226 g005
Figure 6. Shear strength testing die for fiber arrangement.
Figure 6. Shear strength testing die for fiber arrangement.
Micromachines 13 01226 g006
Figure 7. Grinding qualified resin matrix spline.
Figure 7. Grinding qualified resin matrix spline.
Micromachines 13 01226 g007
Figure 8. Unidirectional CFRC lamina.
Figure 8. Unidirectional CFRC lamina.
Micromachines 13 01226 g008
Figure 9. Model geometry of tensile specimen: (a) category A, (b) category B, and (c) category C.
Figure 9. Model geometry of tensile specimen: (a) category A, (b) category B, and (c) category C.
Micromachines 13 01226 g009
Figure 10. XFEM damage comparison for 3P/0° specimen.
Figure 10. XFEM damage comparison for 3P/0° specimen.
Micromachines 13 01226 g010
Figure 11. Crack propagation in vertical direction.
Figure 11. Crack propagation in vertical direction.
Micromachines 13 01226 g011
Figure 12. Modulus of elasticity of categories.
Figure 12. Modulus of elasticity of categories.
Micromachines 13 01226 g012
Figure 13. Matrix and fiber breaking.
Figure 13. Matrix and fiber breaking.
Micromachines 13 01226 g013
Table 1. Material properties and homogenized values.
Table 1. Material properties and homogenized values.
LinearEm (GPa)υmυp
Elasticity3.73 ± 0.300.38 ± 0.010.30
Damageσft (MPa)σft (MPa)GIC (J/m2)
Model61.6 ± 4.6300 ± 30.6334.1 ± 73
Table 2. Material properties and homogenized values for the development of RVE.
Table 2. Material properties and homogenized values for the development of RVE.
MaterialElastic Modulus (GPa)Poisson’s RatioDensity (kg/m3)
Carbon Fiber
(T-300)
2310.21760
Matrix
(Cast Epoxy)
3.4250.321250
Table 3. Mechanical properties and parameters for the unidirectional set for the interfacial debonding model [17].
Table 3. Mechanical properties and parameters for the unidirectional set for the interfacial debonding model [17].
E11 (MPa)E22 (MPa)G12 (MPa)G23 (MPa)υ12XT (MPa)XC (MPa)S
Carbon Fiber22113.894.80.235282500-
Epoxy Resin3.53.51.2961.2960.3511224189.6
Table 4. Mechanical properties of carbon fiber in analysis.
Table 4. Mechanical properties of carbon fiber in analysis.
E11E22 = E33G12 = G13G23v12
245 (GPa)19.8 (GPa)29.191 (GPa)5.922 (GPa)0.28
Table 5. Different dimensions of RVE.
Table 5. Different dimensions of RVE.
ParametersSymbolValueUnits
Fiber Radius r f 4.5µm
RVE Length a 14.5621µm
Nano-Fiber Radius r n f 0.1µm
Table 6. Shear dominated failure assessment.
Table 6. Shear dominated failure assessment.
Shear Dominated FailureNU-DanielApplied Criteria [14]
σ 22 T r a n σ 22 0 τ 12 S L 2 + 2 α σ 22 S L = 1 τ 12 S L 2 + α σ 22 Y T = 1 α = Y T σ 22 T r a n τ 12 T r a n S L 1
Table 7. Comparison with criterion.
Table 7. Comparison with criterion.
TypeE22Relative Error %YT (MPa)Error %YC(MPa)Error %
Tsai Hill9.101.56%117.82.33%25410.92%
Tsai Wu9.101.56%80.530.07%26817.03%
Modified criteria9.101.56%87.523.99%18917.47%
Table 8. 0 deg specimen dimensions.
Table 8. 0 deg specimen dimensions.
Specimen AWidth (mm)Thickness (mm)Failure Load (kN)
112.12.9643.21
212.22.9452.73
312.32.8742.46
412.52.9547.48
Table 9. 90 deg specimen dimensions.
Table 9. 90 deg specimen dimensions.
Specimen BWidth (mm)Thickness (mm)Failure Load (kN)
124.992.760.872
225.212.640.974
325.382.590.914
425.142.530.987
Table 10. 0/90/0/90/0 specimen dimensions.
Table 10. 0/90/0/90/0 specimen dimensions.
Specimen CWidth (mm)Thickness (mm)Failure Load (kN)
124.53.065.62
224.72.9872.65
325.42.9075.17
Table 11. Comparison between predicted and theoretical values with reference to interface debonding.
Table 11. Comparison between predicted and theoretical values with reference to interface debonding.
ParameterLongitudinal Strength (MPa)Predicted Values (MPa)Error %
E₁₁150.8155.22.92%
XT248924192.81%
XC176917252.49%
Table 12. Comparison of theoretical and experimental strength.
Table 12. Comparison of theoretical and experimental strength.
TypesOrientationsTheoretical Strength (MPa)Experimental Strength (MPa)%Diff.
ATensile 0°1074893.716.78
BTensile 0°, 90°, 0°, 90°642593.67.50
CTensile 90°454534.317.62
Table 13. Strength comparison.
Table 13. Strength comparison.
GroupAverage Strength
MPa
Error %% Variation
0° 3P specimen
0° 3P test1382.45--13.79
RVE2112.913.4411.12
0° 3P simulation1829.38.7 13.2
+45°, −45° 3P specimen
+45°, −45° test178.5--12.40
ABAQUS RVE208.338.64.5
+45°, −45° 3P simulation191.611.224.9
(0˚, 90˚, 0˚, 90˚) 3P specimen
Laminate test (Experimental)647.111.334.55
Laminate simulation571.815.77%2.45
Table 14. Experimental results.
Table 14. Experimental results.
SpecimenFailure Load
(kN)
Failure Strength
(MPa)
Young’s Modulus E
(GPa)
Category A
143.211326.8108
252.731569.1
342.461258.09
447.481423.6
Category B
10.87213.706.44
20.97415.25
30.91414.73
40.98715.39
Category C
165.62890.171.01
272.65985.09
375.171019.20
Table 15. Summarized results.
Table 15. Summarized results.
Sample S/No.Width/mmThickness/mmArea/mm2Max Force (N)Tensile Strength (MPa)Modulus (GPa)
14.943.7218.38294.516.033.920
25.503.2017.60235.513.384.043
34.982.2811.35134.511.853.498
45.092.9314.91190.012.743.916
55.094.2121.43246.511.505.106
Average values 13.094.096
Standard deviation 1.7950.601
Discrete Coefficient 13.7014.66
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Swati, R.F.; Qureshi, S.R.; Sohail, M.U.; Munir, A.; Qureshi, O.M.; Khan, A.A. Micro-Mechanical Investigation of Interfacial Debonding in Carbon Fiber-Reinforced Composites Using Extended Finite Element Method (XFEM) Approach. Micromachines 2022, 13, 1226. https://doi.org/10.3390/mi13081226

AMA Style

Swati RF, Qureshi SR, Sohail MU, Munir A, Qureshi OM, Khan AA. Micro-Mechanical Investigation of Interfacial Debonding in Carbon Fiber-Reinforced Composites Using Extended Finite Element Method (XFEM) Approach. Micromachines. 2022; 13(8):1226. https://doi.org/10.3390/mi13081226

Chicago/Turabian Style

Swati, Raees Fida, Saad Riffat Qureshi, Muhammad Umer Sohail, Adnan Munir, Omer Masood Qureshi, and Abid Ali Khan. 2022. "Micro-Mechanical Investigation of Interfacial Debonding in Carbon Fiber-Reinforced Composites Using Extended Finite Element Method (XFEM) Approach" Micromachines 13, no. 8: 1226. https://doi.org/10.3390/mi13081226

APA Style

Swati, R. F., Qureshi, S. R., Sohail, M. U., Munir, A., Qureshi, O. M., & Khan, A. A. (2022). Micro-Mechanical Investigation of Interfacial Debonding in Carbon Fiber-Reinforced Composites Using Extended Finite Element Method (XFEM) Approach. Micromachines, 13(8), 1226. https://doi.org/10.3390/mi13081226

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