# Stress Distribution of CF/EP Laminated Composites under Supercritical Conditions

Authors: Haihong Huang; Zhenwen Li; Huanbo Cheng; Yanzhen Yin
Journal: International Journal of Polymer Science
 (2015)

---

## Abstract

Enormous amounts of wastes have been produced due to extensive use of carbon fiber/epoxy resin (CF/EP) composites. The fact that the supercritical fluid can be used to recycle these composites efficiently has attracted widespread concerns. A three-dimensional model of CF/EP laminates considering the interfacial layers was established. The internal stress distribution of laminates was simulated based on a heat transfer model; and the change of shear stress with supercritical temperature and pressure was investigated. The results show that the shear stress concentration was located in the interfacial layers; the maximum shear stress can be expressed by a curve of convex parabola to the temperature; and the most serious damage occurred in interfacial layers when temperature approached the glass-transition temperature of resin.


---

## Body

## 1. Introduction

Carbon fiber/epoxy resin (CF/EP) composites are used in a wide range of applications in industries such as aerospace, automotive and renewable energy, due to their interesting combination of properties, strength, durability, high strength-to-weight ratios, and resistance. The aircraft industry was an impressive example, with the Boeing 787 “Dreamliner” and Airbus A350 having up to 50% of their weight in carbon fiber reinforced plastics (CFRP) [1]. And the structural parts commonly made of composite materials in automobiles include composite modular front end, tail doors, side doors, construction of body, chassis, interiors, and seating. The average CFRP weight content of a vehicle in Europe was around 24 kg [2]. However, their inherent heterogeneous nature of the matrix and reinforcement leads to the thorny problem of waste disposal. With the advancement of sustainable development strategy in the 21st century, many countries encourage and support clean production, thereby reducing discharge of solid wastes. The worldwide productivity of carbon fiber reached approximately 100,000 tons, while the sales amount reached 50,000 tons in 2013. The increasing use of CFRP also leads to an increasing amount of wastes. In 2010, the global waste carbon fiber products have reached 20,000 tons, while the recycle rate is merely 10%. These wastes contained carbon fibers more than 10,000 tons [3]. It is urgent to explore and industrialize the recycle and reuse of the CF/EP wastes. CF/EP composites have accounted for more than 90% of CFRP based on the investigation [4]. Therefore, the environmental harmless recycle and reuse of CF/EP composites have great environmental significance and economic value.Three technologies have been proposed to recycle CF/EP composites: mechanical recycling, thermal processing, and chemical recycling [5]. Chemical recycling has been investigated using supercritical fluids as a reactive-extraction media, since they possess an interesting combination of properties such as low viscosity, high mass transport coefficients, high diffusivity, and salvation power [6]. The University of Nottingham developed a supercritical solvent decomposition method, and the recycled carbon fiber retained its original strength (up to 98%) [7, 8]. When the potassium hydroxide (KOH) was used as catalyst, yield of resin elimination reached 95%. Okajima et al. studied the decomposition of epoxy resin in the near-critical and supercritical water, in which 380°C and 25 MPa were determined preliminarily to be the optimal processing conditions [9]. Liu et al. found that supercritical water was a realistic option for the depolymerisation of epoxy resin and its composites [10].Finite element method (FEM) has been proved to be feasible and effective to simulate the stress field distribution of composites. Lu et al. acquired thermal-stress distribution of CFRP during thermal cycling using micromechanical finite element model [11]. Gao et al. carried out numerical simulation of internal thermal stress of unidirectional CF/EP composites using two-dimensional finite element model and found that internal thermal stress of the interfacial region was obviously concentrated [12]. Li et al. simulated the stress distribution for printed circuit board (PCB) in the supercritical CO2 environment and found that the internal delamination was determined by shear stress and peel stress [13]. However, the stress distribution of CF/EP composites under supercritical conditions has not been studied in the context of recycling.This paper selects CF/EP laminated composites as the research object. An equivalent three-dimensional model was established, and the model of heat transfer was proposed under supercritical condition. Thermal-stress coupling analysis was completed based on FEM, in which the supercritical fluid temperature and pressure were set as the boundary conditions of load. The internal stress distribution of laminates was analyzed. With changes of supercritical environment, that is, temperature and pressure, the variation of maximum shear stress was studied. The advantage of this simulation is that it is possible to identify the impact of different supercritical conditions on structural failure of composites, which provides useful guidelines for designing the recycling process of CF/EP composites.


## 2. Equivalent Model and Heat Transfer Process Of CF/EP Laminates

The equivalent three-dimensional model considering the interfacial layers is the key for the finite element simulation. Traditional modeling holds that the two layers of CF/EP laminates are connected by resin matrix, and the interlaminar properties depend largely on properties of matrix. However, the existence of interface between fiber and resin is neglected. Interface is a very important microcosmic structure, which not only connects fiber and matrix but also transfers loads from matrix to fiber. So it is very crucial to consider the existence of the interface. Interface is not a single geometric surface between reinforcement and matrix but the transitional region from matrix surface to reinforcement surface. According to the hybrid model, a single fiber is selected as a basic cell, called representative volume element (RVE), as shown in Figure1(a) [14]. Based on the assumption that the transitional region is equivalent to the interface layer, an equivalent three-dimensional model of CF/EP laminates was established. The simplified laminates model was symmetric, which was laid by resin layer, interfacial layer, and fiber layer in turn. Equivalent three-dimensional model of CF/EP laminates was shown in Figure 1(b).(a) Representative volume element (RVE); (b) equivalent model of laminates.
(a)
(b)Heat transfer model of a lamina was shown in Figure2(a). It was assumed that the size of lamina is infinite along x direction. According to Fourier laws of heat conduction, the heat flux of laminate Q1 can be calculated by (1) [15]. Consider
(1)Q1=-kAdtdx=-kAt2-t1l=t1-t2l/kA,
where l is the single layer thickness of CF/EP, k is the thermal conductivity, A is the heat transfer area, t1 and t2 are the surface temperature, respectively, and l·k-1·A-1 is the thermal resistance.(a) Steady-state conduction of a lamina; (b) steady-state conduction of laminates.
(a)
(b)The heat flux of laminates with different thickness and physical properties was calculated. The size and thermophysical properties of laminates were shown in Figure2(b). The heat was transferred from the left to right in the model. Based on the series connection of thermal resistance, the heat flux of laminates Q2 can be calculated by
(2)Q2=t1-tn∑i=1nli/kiA,
where n is the layer number of laminates and t1 and tn are the surface temperature for two sides of laminates.Supercritical fluid was closed to the surfaces of laminates under supercritical conditions. In this case, thermal resistance of supercritical fluid should be taken in consideration. So the heat flux of laminatesQ3 can be calculated by the following:
(3)Q3=tf1-tfn1/h1A+∑i=1nli/kiA+1/hnA,
where tf1 and tfn are the supercritical fluid temperature for the two sides of laminates and h1 and hn are the surface convective heat transfer coefficient.


## 3. Thermal-Stress Coupling Analysis

Thermal-stress coupling analysis of the CF/EP laminates under supercritical conditions was simulated by ABAQUS. Temperature was calculated by transient thermal analysis to obtain the node temperature distribution firstly. And then, taking temperature field as thermal loads, the stress was calculated by structural analysis. The simulation was executed based on the following assumptions: (1) the change of convective heat transfer coefficient along with temperature and pressure was ignored; namely, the surface convective heat transfer coefficient was constant; (2) the reverse-coupling effect of deformation process to temperature field was ignored; namely, the process was regarded as one-way sequentially coupled analysis.


### 3.1. Finite Element Model

The laminates model with size of 30 mm × 30 mm × 1.66 mm was established. As shown in Figure3, it consisted of five fiber layers, six resin layers, and ten interfacial layers. The stacking sequence of fiber layers was [0/90/0/90/0]. To simplify calculation, 1/4 of CF/EP laminated composites was analyzed due to its symmetry. Resin layers and interface layers were considered as isotropic homogeneous material, and fiber layers were considered as transversely isotropic material. Physical parameters of resin layers and fiber layers were given by reference [16], and parameters of the interfacial layers were given based on the rule of the general simple mixture [17], as shown in Table 1. The convective heat transfer coefficient of supercritical fluid was 4000 W/(m2·°C).Table 1
The property of CF/EP laminates at room temperature.
Parameter
Fiber layer
Resin layer
Interfacial layer
Thickness (mm)
0.15
0.15
0.001
Density (kg/m3)
1600
1200
1450
Thermal conduction coefficient (W/(m·°C))
k
1
=
4.6, k2=k3=0.42
0.2
2.4
Thermal expansion coefficient (/°C)
α
1
=
-
1
×
10
-
6, α2=26×10-6
60
×
10
-
6
27
×
10
-
6
Poisson ratio
ν
12
=
0.28, ν23=0.4
0.35
0.315
Elastic modulus (MPa)
E
1
=
138
×
10
3, E2=E3=11×103
4
×
10
3
7.5
×
10
3
Specific heat (J/(g·°C))
0.48
1.2
0.84Figure 3
Model of CF/EP laminates and interface stress element.


### 3.2. Analysis of Simulation Results

As shown in Figure4(a), the deformation of CF/EP laminates appeared under the thermal load, and the maximum displacement was more than two times the thickness of laminates at four corners. The thermal expansion coefficient of resin layer was far more than that of fiber layer along the length of fibers. Owing to the mismatch of thermal expansion coefficient, significantly interfacial slip deformation appeared, as shown in Figure 4(b).(a) Displacement distribution pattern of laminates at 300°C, 20 MPa; (b) displacement distribution pattern of local laminates at 300°C, 20 MPa.
(a)
(b)The shear stressτxz, τyz and peel stress σz were caused by interfacial slip and mainly contributed to the interfacial failure of laminates. When the temperature of supercritical environment was set as 300°C, the compressive stress increased significantly with the increase of pressure, but both shear stresses τxz and τyz had little change, as shown in Figure 5(a). Due to the compressive stress acting in the opposite direction of peel stress σz, the delamination of CF/EP laminates was restrained gradually. When the pressure of supercritical environment was set as 10 MPa, both shear stresses τxz and τyz increased significantly with the increase of temperature, as shown in Figure 5(b). The shear stress τyz increased from 2.9 MPa to 14.7 MPa when the temperature goes from 50°C to 325°C.The impact of (a) pressure on the maximum shear/peel stress; and (b) temperature on the maximum shear stress.
(a)
(b)The shear stress became the main resource to result in structural damage of laminates, due to the fact that the peel stress manifested as compressive stress in supercritical conditions. As shown in Figure6, the shear stress concentrated mainly on the interfacial layers, especially near the free edge, and the stress decreased rapidly from the free edge to the inside. At either side of 0° plies, shear stresses τxz were concentrated on the free edge of interfacial layer and their directions were opposite; and it was the same situation for the shear stresses τyz at either side of 90° plies.(a) Shear stressτxz of local laminates; (b) shear stress τxz of local interfacial layer; (c) shear stress τyz of local laminates; and (d) shear stress τyz of local interfacial layer.
(a)
(b)
(c)
(d)


### 3.3. Discussion

Under the supercritical environment, the mechanical properties of resin were changed by the combined effect of supercritical fluid and temperature, while the effects of pressure can be neglected since there was little change of composite properties when increasing the pressure. In addition, the Elastic modulus of resin had considerable effect on the stress. The Elastic modulus of epoxy resin plummeted about 2-3 orders of magnitude, when it transformed from glassy state to rubbery state with the increase of temperature [18]. Although the glass-transition temperature of epoxy resin was affected by the composition, curing agent, and crosslinking degree, it was about 200°C when the crosslinking degree was 70%–90% based on the published researches [19–23]. The Elastic modulus and Poisson ratio of laminates in different temperatures were shown in Table 2.Table 2
Elastic modulus and Poisson ratio of laminates at different temperature.
Temperature (°C)
50
100
150
200
250
275
300
325
Elastic modulus (MPa)
3.97
×
10
3
3.5
×
10
3
3.1
×
10
3
2.6
×
10
3
1.8
×
10
3
1.3
×
10
3
0.5
×
10
3
0.07
×
10
3
Poisson ratio
0.353
0.362
0.384
0.4
0.421
0.43
0.442
0.453The shear stress of laminates under different supercritical conditions was calculated when the pressure of supercritical environment was set as 10 MPa, as shown in Figure7. The change of Elastic modulus and Poisson ratio with the temperature was taken into consideration in the calculation. It showed that the maximum shear stress can be expressed by a curve of convex parabola to the temperature, and the peak of shear stress appeared at 250°C. With the increase of supercritical fluid temperature, the maximum shear stress increased gradually at first and then decreased rapidly when the temperature was higher than the glass-transition temperature. The interfacial shear stress concentration was weakened due to the reduction of Elastic modulus.Figure 7
The variation of maximum shear stress when considering the change of Elastic modulus and Poisson ratio with temperature.


### 3.4. The Recycling of CF/EP Laminates

The simulation verifies that the CF/EP laminates can be recycled by the supercritical fluids. Under the supercritical conditions, the supercritical fluid infiltrates into laminates through initial pores, which leads to the swelling behavior of resin; and the shear stress concentration appeared because of the slip deformation. When the shear force was greater than bonding force, the bonding interface among layers would be damaged. Supercritical fluid flowed rapidly through the new channels, which not only accelerated transmission of heat but also widened the extraction channel of reaction products. Consequently, the structural and interfacial bonding damaged more seriously. However, the interfacial shear stress decreased rapidly with the decrease of Elastic modulus when the temperature was higher than glass-transition temperature, which would make the recycling of composites inefficient. Therefore, the recycling process of CF/EP laminates should be operated near the glass-transition temperature of the resin.


## 3.1. Finite Element Model

The laminates model with size of 30 mm × 30 mm × 1.66 mm was established. As shown in Figure3, it consisted of five fiber layers, six resin layers, and ten interfacial layers. The stacking sequence of fiber layers was [0/90/0/90/0]. To simplify calculation, 1/4 of CF/EP laminated composites was analyzed due to its symmetry. Resin layers and interface layers were considered as isotropic homogeneous material, and fiber layers were considered as transversely isotropic material. Physical parameters of resin layers and fiber layers were given by reference [16], and parameters of the interfacial layers were given based on the rule of the general simple mixture [17], as shown in Table 1. The convective heat transfer coefficient of supercritical fluid was 4000 W/(m2·°C).Table 1
The property of CF/EP laminates at room temperature.
Parameter
Fiber layer
Resin layer
Interfacial layer
Thickness (mm)
0.15
0.15
0.001
Density (kg/m3)
1600
1200
1450
Thermal conduction coefficient (W/(m·°C))
k
1
=
4.6, k2=k3=0.42
0.2
2.4
Thermal expansion coefficient (/°C)
α
1
=
-
1
×
10
-
6, α2=26×10-6
60
×
10
-
6
27
×
10
-
6
Poisson ratio
ν
12
=
0.28, ν23=0.4
0.35
0.315
Elastic modulus (MPa)
E
1
=
138
×
10
3, E2=E3=11×103
4
×
10
3
7.5
×
10
3
Specific heat (J/(g·°C))
0.48
1.2
0.84Figure 3
Model of CF/EP laminates and interface stress element.


## 3.2. Analysis of Simulation Results

As shown in Figure4(a), the deformation of CF/EP laminates appeared under the thermal load, and the maximum displacement was more than two times the thickness of laminates at four corners. The thermal expansion coefficient of resin layer was far more than that of fiber layer along the length of fibers. Owing to the mismatch of thermal expansion coefficient, significantly interfacial slip deformation appeared, as shown in Figure 4(b).(a) Displacement distribution pattern of laminates at 300°C, 20 MPa; (b) displacement distribution pattern of local laminates at 300°C, 20 MPa.
(a)
(b)The shear stressτxz, τyz and peel stress σz were caused by interfacial slip and mainly contributed to the interfacial failure of laminates. When the temperature of supercritical environment was set as 300°C, the compressive stress increased significantly with the increase of pressure, but both shear stresses τxz and τyz had little change, as shown in Figure 5(a). Due to the compressive stress acting in the opposite direction of peel stress σz, the delamination of CF/EP laminates was restrained gradually. When the pressure of supercritical environment was set as 10 MPa, both shear stresses τxz and τyz increased significantly with the increase of temperature, as shown in Figure 5(b). The shear stress τyz increased from 2.9 MPa to 14.7 MPa when the temperature goes from 50°C to 325°C.The impact of (a) pressure on the maximum shear/peel stress; and (b) temperature on the maximum shear stress.
(a)
(b)The shear stress became the main resource to result in structural damage of laminates, due to the fact that the peel stress manifested as compressive stress in supercritical conditions. As shown in Figure6, the shear stress concentrated mainly on the interfacial layers, especially near the free edge, and the stress decreased rapidly from the free edge to the inside. At either side of 0° plies, shear stresses τxz were concentrated on the free edge of interfacial layer and their directions were opposite; and it was the same situation for the shear stresses τyz at either side of 90° plies.(a) Shear stressτxz of local laminates; (b) shear stress τxz of local interfacial layer; (c) shear stress τyz of local laminates; and (d) shear stress τyz of local interfacial layer.
(a)
(b)
(c)
(d)


## 3.3. Discussion

Under the supercritical environment, the mechanical properties of resin were changed by the combined effect of supercritical fluid and temperature, while the effects of pressure can be neglected since there was little change of composite properties when increasing the pressure. In addition, the Elastic modulus of resin had considerable effect on the stress. The Elastic modulus of epoxy resin plummeted about 2-3 orders of magnitude, when it transformed from glassy state to rubbery state with the increase of temperature [18]. Although the glass-transition temperature of epoxy resin was affected by the composition, curing agent, and crosslinking degree, it was about 200°C when the crosslinking degree was 70%–90% based on the published researches [19–23]. The Elastic modulus and Poisson ratio of laminates in different temperatures were shown in Table 2.Table 2
Elastic modulus and Poisson ratio of laminates at different temperature.
Temperature (°C)
50
100
150
200
250
275
300
325
Elastic modulus (MPa)
3.97
×
10
3
3.5
×
10
3
3.1
×
10
3
2.6
×
10
3
1.8
×
10
3
1.3
×
10
3
0.5
×
10
3
0.07
×
10
3
Poisson ratio
0.353
0.362
0.384
0.4
0.421
0.43
0.442
0.453The shear stress of laminates under different supercritical conditions was calculated when the pressure of supercritical environment was set as 10 MPa, as shown in Figure7. The change of Elastic modulus and Poisson ratio with the temperature was taken into consideration in the calculation. It showed that the maximum shear stress can be expressed by a curve of convex parabola to the temperature, and the peak of shear stress appeared at 250°C. With the increase of supercritical fluid temperature, the maximum shear stress increased gradually at first and then decreased rapidly when the temperature was higher than the glass-transition temperature. The interfacial shear stress concentration was weakened due to the reduction of Elastic modulus.Figure 7
The variation of maximum shear stress when considering the change of Elastic modulus and Poisson ratio with temperature.


## 3.4. The Recycling of CF/EP Laminates

The simulation verifies that the CF/EP laminates can be recycled by the supercritical fluids. Under the supercritical conditions, the supercritical fluid infiltrates into laminates through initial pores, which leads to the swelling behavior of resin; and the shear stress concentration appeared because of the slip deformation. When the shear force was greater than bonding force, the bonding interface among layers would be damaged. Supercritical fluid flowed rapidly through the new channels, which not only accelerated transmission of heat but also widened the extraction channel of reaction products. Consequently, the structural and interfacial bonding damaged more seriously. However, the interfacial shear stress decreased rapidly with the decrease of Elastic modulus when the temperature was higher than glass-transition temperature, which would make the recycling of composites inefficient. Therefore, the recycling process of CF/EP laminates should be operated near the glass-transition temperature of the resin.


## 4. Conclusions

The internal stress distribution of laminates was simulated by FEM under supercritical conditions. The results show that the interfacial slip deformation appeared due to the mismatch of thermal expansion coefficient between resin and fiber. The shear stress concentrated on the interfacial layer near the free edge, and the maximum shear stress can be expressed by a convex parabola to the temperature. The microstructure was damaged more seriously near glass-transition temperature, and the recycling processes of CF/EP laminates under the supercritical conditions were proposed. This stress distribution was simulated only under the condition of physical environment, that is, temperature and pressure, without considering the chemical decomposition effects of supercritical fluids. The presented model and simulations were suitable for the solvents such as water, alcohols, and acetone under supercritical conditions.



---

*Source: 715753-2015-02-08.xml*