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].
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
and
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
represents the length of the sides of the unit cell and
represents the radius of the carbon fiber. Here, VFF is the volume filler fraction.
The relationship is presented in Equation (3), which gives us the length of the unit cell.
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.
where
and
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
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).
In Equation (7)
ε is the strain and
Cd is the elasticity matrix, which reflects any damage and has the form shown in Equation (8).
In Equation (8)
,
is the current state of fiber damage,
is the current state of matrix damage,
is the current state of shear damage,
is the Young’s modulus in fiber direction (longitudinal direction),
is the Young’s modulus in the direction perpendicular to the fibers (transverse direction),
G12 is the shear Modulus and
, and
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
,
, and
were derived, from the damage variables
,
,
, and
shown in Equations (9)–(13) which were referred to in the four failure modes discussed previously.
The damage variables ,,, and 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 Compression:
Matrix Compression:
In Equations (14)–(17)
,
, and
are the longitudinal, transverse, and shear stresses in the lamina,
and
refers to the tensile and compression strength in the fiber direction (longitudinal tensile and compressive strength), and
and
refers to the tensile and compression strength in the transverse direction.
and
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.
where
ζd is interface debond energy,
F is fiber load at matrix crack plane,
v(
x) is relative displacement between fiber and matrix, and
(0) shows the fiber axial displacement at the matrix crack plane. The interface debond length
is determined by equation.
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
are tensile, and
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.