Next Article in Journal
Redox-Stable and Multicolor Electrochromic Polyamides with Four Triarylamine Cores in the Repeating Unit
Previous Article in Journal
Injection Pultrusion of Glass-Reinforced Epoxy: Cure Kinetics, Rheology, and Force Analysis
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Three-Dimensional Thermo–Chemo–Mechanical Coupled Curing Analysis for the Filament Wound Composite Shell

Zhejiang-Italy Joint Lab for Smart Materials and Advanced Structures, School of Mechanical Engineering and Mechanics, Ningbo University, Ningbo 315211, China
*
Authors to whom correspondence should be addressed.
Polymers 2024, 16(12), 1643; https://doi.org/10.3390/polym16121643
Submission received: 27 April 2024 / Revised: 1 June 2024 / Accepted: 6 June 2024 / Published: 10 June 2024
(This article belongs to the Section Polymer Composites and Nanocomposites)

Abstract

:
Carbon fiber resin-based composite materials are widely employed in the manufacturing of composite shells. During the curing process, the temperature gradients and cure degree gradients make it easy to generate thermal strains in both carbon fibers and resin, with the resin experiencing cure shrinkage strain due to the curing reaction, ultimately leading to residual stresses and strains. In this paper, a three-dimensional thermo–chemo–mechanical coupled curing model of the composite shell was established based on a resin test, and the changes of temperature, curing degree, residual stress, and strain during the solidification of the composite shell were investigated. First, the curing property parameters and elastic modulus of HCM-2184 resin were obtained through a curing dynamic test and a tensile test. Then, considering the heat release and shrinkage reaction of solidification, a coupled thermo–chemo–mechanical curing model was developed with the CHILE ( α ) elastic model, and the curing process of the composite shell was simulated numerically. The results show that the resin used in the test belongs to the autocatalytic reaction. For thin composite shells, the heat accumulation inside the shell during curing is not obvious. During the curing process, the curing shrinkage behavior of the resin is an important factor for the generation of residual stress and residual strain.

1. Introduction

Composite materials, which have been widely applied in recent years, typically comprise two or more materials with distinct physical or chemical properties at the macroscopic scale. These materials not only retain the advantages of their constituent components, but also exhibit superior comprehensive performance compared to single materials at the macro level. Compared to traditional metallic materials and alloys, carbon fiber resin-based composites possess advantages such as low density, high-temperature resistance, corrosion resistance, high strength, high modulus, fatigue resistance, and a low thermal expansion coefficient [1,2,3,4], making them extensively used in aerospace and other fields [5]. For instance, at the same volume, ensuring strength and safety while reducing mass can effectively decrease the load. Additionally, in high-temperature and high-humidity environments, the material properties remain substantially unchanged, aligning with the lightweight manufacturing requirements of modern solid rocket engines, hence justifying their use in engine shell fabrication. The resin, serving as the matrix of composite materials, exhibits curing behavior, a chemical change where under external heating, internal groups within the resin undergo cross-linking reactions, releasing heat and transitioning from a viscous flow state to a glassy state [6], with an increase in viscosity and volume shrinkage. During the curing process of composite material shells, the presence of thermal expansion and contraction between carbon fibers and resin, resin cure shrinkage [7,8,9], and mismatches in the thermal expansion coefficients between different parts of the shell and mold materials [10] lead to thermo–chemo–mechanical coupled residual stresses and strains during the curing and molding process. After curing, some residual stresses remain within the shell, causing deformation that deviates from the designed dimensions, affecting the molding quality of shell components and their performance indicators. For example, in the strength prediction of the shell, the residual stresses and strains formed during curing can be considered to be prestress and prestrain. Therefore, the analysis of shell curing, as a factor affecting the performance of molded shells, necessitates the study of the evolution and distribution patterns of residual stresses and strains generated during curing.
The analysis of the structural quality of composite material components during curing molding essentially involves the study of changes in internal stresses and strains within the components caused by the combined effects of carbon fiber and resin thermal expansion, as well as resin curing, under the segmented processes of heating, holding, and cooling during the curing process. Currently, numerous researchers have employed experimental and numerical simulation approaches to research the curing process of composite materials based on thermosetting resins [11,12,13,14,15]. The study of resin cure kinetics forms a crucial foundation for investigating the curing deformation of composite materials. Differential Scanning Calorimetry (DSC) is commonly used to study resin curing behavior. Lee et al. [16] obtained the kinetic parameters of 3501-6 resin by DSC. Zanjanijam [17] studied the curing kinetics of phenolic resin/coal coke composites.
The curing process significantly impacts the quality of curing molding of composite materials. On the one hand, temperature affects the curing progression of the resin, leading to phenomena of uneven curing. On the other hand, the properties of resin [18,19] also change with temperature variations, such as thermal conductivity, density, and specific heat capacity. Pedersen et al. [20] investigated the thermal diffusivity of epoxy resin during the curing process, and discovered that thermal conductivity and specific heat capacity have functional relationships with the degree of cure. Johnston [21] studied the density and various thermophysical parameters of AS4/8552 composite materials.
The simulation of composite materials’ curing has undergone a developmental process [22,23,24]. Cheung [25] integrated thermal conduction with resin cure kinetics to study the temperature field distribution of composite laminate boards during the curing process. Abdelal [26] utilized thermodynamic parameters that change with temperature to construct a thermo–chemical coupled model, calculating the temperature and cure degree distribution of composite material components. To accurately simulate and predict the evolution of stress and strain during the curing process of composite material components, it is essential to establish their constitutive model throughout the entire process, and a substantial number of scholars have researched the constitutive models of composite materials. Among them, the main categories are linear elastic and viscoelastic constitutive models. Bogetti and Gillespie [27] proposed a curing-hardened instantaneous linear elastic model, assuming that the resin behaves linearly elastically at each incremental step during the curing process, and considering that the elastic modulus has a functional relationship with the degree of cure, hence referred to as the CHILE ( α ) model. However, the above methods did not account for the viscoelasticity of composite materials and the stress relaxation it causes. In order to more accurately characterize the stress relaxation behavior of composite materials, researchers introduced viscoelasticity into the calculations for the curing of composite materials [28,29,30], with Ding et al. [31,32] establishing a thermo–viscoelastic constitutive model considering stress relaxation based on a generalized Maxwell model with n Maxwell elements. The thermo–chemo–mechanical coupled curing process of laminates was numerically simulated. Nonetheless, due to the significant difficulties in measuring the viscoelastic parameters of composite materials and the increased complexity of numerical calculations involving viscoelastic constitutive equations, the CHILE ( α ) model, by comparison, is relatively simple and can more accurately predict the stress evolution during the curing of composite materials, and is thus widely adopted by scholars.
As aforementioned, the fully coupled analysis for thermo–chemo–mechanical curing is very important. While at present, there are few numerical simulation studies on the application of specific types of resins in the curing of filament wound composite shells, this article addresses the issue of residual stresses and strains leading to structural deformation during the curing process of composite material shells. The test was conducted using HCM-2184 resin. Differential scanning calorimetry (DSC) was used to test the curing kinetics of the resin, and the curing properties of the resin were obtained. A tensile test was carried out to study the tensile properties of resin-casting parts. Based on the consideration of resin-curing heat release, curing shrinkage, and anisotropy of composite materials, a thermo–chemo–mechanical coupled curing model of the filament wound composite shell was established with the CHILE ( α ) constitutive model. The evolutionary distribution of temperature, curing degree, and residual strain during curing molding was studied, which provided a reference for the quality evaluation of curing molding.

2. Resin Test

2.1. Preparation of Resin Matrix

Material: HCM-2184 epoxy resin and curing agent, XTEND832 semi-permanent release agent. First, the HCM-2184 epoxy resin was mixed with the curing agent at the mass ratio 100:31. Then, through the high speed stirring and mixing evenly, the resin was vacuumed in the vacuum oven at room temperature for an hour. Finally, when no bubbles appeared, the resin was removed, and subsequent tests were carried out.

2.2. Non-Isothermal DSC Test

In order to study the curing reaction kinetics of HCM-2184 epoxy resin, according to the curing characteristics of HCM-2184 epoxy resin, Netzsch DSC-200F3 non-isothermal DSC tests were carried out. The 5 mg mass sample was placed in an aluminum pot, heated to 300 °C at the heating rate of 5.0, 10.0, 15.0, and 20.0 °C/min, and tested under nitrogen atmosphere.

2.3. Tensile Test

When calculating the curing process of composite materials, it is necessary to input the elastic modulus parameters of the resin. Therefore, tensile tests on the resin are required. The tensile test specimens of the resin are prepared according to the current national standard in China [33], with specific dimensions, as shown in Figure 1.

2.4. Experimental Results and Discussions

The DSC heat flow diagram of epoxy resin at different heating rates (5 °C/min–20 °C/min) is shown in Figure 2. It can be seen that there is a single heat release peak at different heating rates, and the peak heat flow and peak temperature increases with the increase in heating rate, and the heat release peak moves towards the higher temperature direction with the increase in heating rate, indicating that the resin has a heat transfer lag phenomenon.
According to the results of Figure 1, the characteristic parameters such as maximum heat release temperature (peak temperature), peak heat flow, and total heat release of curing reaction ( Δ H c u r e , calculated by the area under the peak heat release) are summarized in Table 1, and the figures in brackets are standard deviations. The results of each heating rate were averaged by repeated tests (at least three times), and the total curing heat release calculated at different heating rates was averaged to obtain the final reaction heat of the resin.
To solve the cure kinetics equation of the resin, it is necessary to determine the resin’s activation energy through mathematical fitting. Activation energy is defined as the minimum energy required for a chemical reaction to occur: the smaller the activation energy, the more readily the reaction takes place. Commonly used equations for calculating activation energy include the Kissinger equation and the Ozawa equation [34]. This article employs both equations to calculate the activation energy of the resin, and the average of the two results is taken as the final activation energy.
The Kissinger equation is expressed as follows:
ln ( τ T P 2 ) = ln ( A R E a ) E a R 1 T P ,
In the equation, R represents the ideal gas constant, generally taken as 8.314 J/(mol·K); A denotes the frequency factor, with the unit S 1 . E a is the activation energy of the cure reaction being determined, τ is the heating rate, and T p is the peak temperature obtained from tests under different heating rates.
The Ozawa equation is expressed as follows:
d ln τ d 1 T p = 1.052 E a R ,
In Equation (2), the physical significance of each parameter is the same as that in the Kissinger equation.
Firstly, the activation energy is determined based on the Kissinger equation. Using the data from Table 1, a plot is constructed with 1000 / T p as the x-axis and ln ( τ / T p 2 ) as the y-axis. After linear fitting, the result is presented in Figure 3a. Using the Ozawa equation to determine the activation energy, a plot is created with 1000 / T p as the x-axis and ln ( τ ) as the y-axis. After linear fitting, the result is shown in Figure 4.
In Figure 3a, after linear fitting of the HCM-2184 resin, the slope of the line is 5.554, with R 2 = 0.998 indicating a good fit. The reaction activation energy is calculated according to Equation (1), resulting in E a = 46.18 kJ/mol. In Figure 3b, the slope of the line is 6.472, with R 2 = 0.998. The reaction activation energy is calculated according to Equation (2), E a = 51.15 kJ/mol.
Taking the average of the results calculated from the two equations, the activation energy for HCM-2184 resin is determined to be 48.67 kJ/mol.
When calculating the degree of cure ( α ) during the curing process of thermosetting resin, it can be represented by the ratio of the heat of reaction at the current time to the total heat of reaction, as shown in Equation (3). Among them, the reaction heat at the current moment ( Δ H T ) is obtained by integrating the heat flow curve measured using DSC in Figure 1 over time, and the total heat released during the curing process is shown in Table 1:
α = Δ H T Δ H c u r e ,
The variation trend of the curing degree of HCM-2184 resin with the increase in temperature under different heating rates is obtained, as shown in Figure 4.
By further deriving the degree of cure with respect to time, the relationship between the cure reaction rate of the resin and time at different heating rates is obtained. Curves depicting the degree of cure and cure rate as functions of time under four different heating rates are plotted, as shown in Figure 5. In the figure, the times at which the cure rate reaches its peak are marked vertically, and the corresponding degree of cure intervals are marked horizontally. It can be observed that, at different heating rates, the resin reaches its maximum cure reaction rate within a degree of cure range of 0.4–0.6, indicating that the cure reaction is most intense within this interval.
The N-order reaction model [35] is applicable to the resin system without autocatalytic behavior in the curing reaction, and its expression is shown in Equation (4):
d α d t = K ( T ) ( 1 α ) n ,
where n is the order of reaction associated with the material, K is a function of temperature, and α is the degree of cure.
The curing kinetic model expression of epoxy resin system with autocatalytic characteristics is shown in Equation (5):
d α d t = K ( T ) α m ( 1 α ) n ,
The relationship between the cure reaction rate and the degree of cure for epoxy resin at different heating rates is shown in Figure 6.
As shown in Figure 6, the cure reaction rate is related to the heating rate, with the cure rate significantly accelerating as the heating rate increases. At four different heating rates, the cure reaction rate exhibits a parabolic trend of first increasing and then decreasing, rather than reaching the maximum value at the outset. The interval at which the cure reaction rate peaks is within the range of 0.4–0.6, largely unaffected by the heating rate. This indicates that the basic model of the resin’s cure reaction is independent of the heating rate, thereby suggesting that the cure kinetics model for the HCM-2184 resin system belongs to an autocatalytic reaction [35].
The cure rate curves obtained from the experiments were subjected to nonlinear fitting based on the autocatalytic reaction model. Studies indicate that the extremum of the resin’s autocatalytic reaction generally occurs at a degree of cure around 0.4. Therefore, to fit the curves more accurately, a degree of cure of 0.4 was used as a demarcation point for segmental fitting. Figure 7 presents a comparison between the experimental and fitting results of the cure rate of HCM-2184 resin at different heating rates. It is observable that, at all four heating rates, the experimental and fitted curves are in basic agreement, indicating that fitting each curve individually yields good results. Subsequently, by averaging the characteristic parameters of the cure kinetics equations fitted at different heating rates, a universal formula is obtained as follows:
d α d t = 3802 exp ( 48670 8.314 T ) α 0.539 ( 1 α ) 1.002 , α 0.4 3891 exp ( 48670 8.314 T ) α 0.404 ( 1 α ) 1.264 , α > 0.4 ,
Figure 8 presents the comparison results of the curves, revealing that, although there is a certain discrepancy between the averaged curves and the experimental curves, they still essentially coincide. At the relatively lower heating rates of 5 °C/min and 10 °C/min, the fitted curves and experimental curves almost overlap. At the higher heating rates of 15 °C/min and 20 °C/min, there is a noticeable difference in the maximum values of the cure rates of the fitted curves, and yet the general trend of the cure rate remains fundamentally the same.
For the tensile test, the force-displacement curve measured using the test is shown in Figure 9a, and the stress–strain curve is drawn with the linear section. The result after linear fitting is shown in Figure 9b. The slope—that is, the elastic modulus of the resin after curing—is 1.97 GPa, which is used as the material parameter for the subsequent curing simulation.

3. Construction of the Composite Curing Coupled Model

3.1. The Thermo–Chemo–Mechanical Coupled Model

The curing process of resin-based carbon fiber composite material is essentially the result of the joint action of the temperature applied externally and the internal heat source generated from the resin releasing heat. The external applied temperature belongs to the heat conduction, and the curing of the resin belongs to the chemical reaction. The heat conduction and curing reaction in composites lead to the generation of residual stress and strain, resulting in curing deformation. Therefore, the curing process is a thermo–chemo–mechanical coupled process. The three-dimensional Fourier heat conduction equation of the anisotropic composite and the internal heat source of the resin are shown as follows:
ρ c C p T t = k 11 2 T x 2 + k 22 2 T y 2 + k 33 2 T z 2 + q ˙ 0 ,
q ˙ 0 = ρ r V r H r d α d t ,
where ρ c is the density of the composite material, k i i denotes the heat conduction coefficient of the composite material, C p is the specific heat capacity, q ˙ 0 is the internal heat source generated during the curing process of the resin, ρ r is the density of the resin, V r is the volume fraction of the resin, H r is the total heat released after the resin is completely cured by the test, and d α / d t is the curing rate.
For the resin used in this paper, a semi-empirical phenomenological kinetic model based on macro-scale was adopted. According to the test, HCM-2184 was determined to be an autocatalytic resin, and the curing rate equation is shown as Equation (6).
The mechanical behavior during the curing process of thermo–chemo–mechanical coupled, that is, the constitutive relationship of the composite material, is described in Section 3.2.

3.2. Constitutive Relationships of the Composite Material

In carbon fiber-reinforced resin composite materials, the carbon fibers acting as the reinforcement material only undergo thermo–mechanical changes during the curing process, and their constitutive behavior is linear elastic. In contrast, the resin, serving as the matrix, is a polymer material that undergoes chemical changes [36] during curing, and generally exhibits viscoelastic behavior. Scholars have proposed computational models for the evolution of modulus during the resin’s curing process, such as the CHILE model, the path-dependent model, and the viscoelastic generalized Maxwell model, among which the CHILE model is widely adopted. It assumes that the resin can be approximately considered to be a linear elastic material within a single time step, with the model linearly varying with the degree of cure, as shown in Equation (9).
E r = E r 0 ( 1 α m o d ) E r 0 + α m o d E r E r     α < α c 1 α c 1 α < α c 2 α α c 2 ,
where α m o d = ( α α g e l ) / ( 1 α g e l ) , α g e l is the gel point when the resin starts crosslinking. E r 0 and E r are the elastic moduli of the resin in the uncured state and when curing is completed, respectively. E r / E r 0 is usually set to be one thousand. In this case, E r = 1.97 GPa, E r 0 = 1.97 MPa. α c 1 , and α c 2 are the curing degrees corresponding to the starting and ending points of curing shrinkage, respectively.
The composite shell is composed of fiber winding, which is mainly divided into a barrel body, a pole hole, and a head. Fiber winding is divided into two types: spiral winding, and circumferential winding. The circular winding is only distributed in the body, and the angle is about 90°. In the composite shell model with spiral fiber winding, the winding angle from the barrel to the pole hole is calculated for the fiber winding layer as follows:
α ( r ) = sin 1 ( r 0 r ) ,
where r 0 is the radius at the polar hole of the shell, and r is the radius of the fiber winding layer at any position of the head section. When r   =   r 0 , the winding angle is 90°, that is, the fiber winding layer is tangent to it at the polar hole. The initial winding angle of the spiral winding layer in the barrel body is 19.5°. When approaching the pole hole, the winding angle increases to 90°. The circumferential winding layer is only distributed in the barrel body, and the angle is 90°, as shown in Figure 10a. Figure 10b shows the change of the winding angle of the spiral winding layer along the axial direction of the shell.
In the actual winding process, the direction of the fibers continuously changes with the winding angle. Therefore, it is unreasonable to place the direction of the fiber material in the global coordinate system, which requires two coordinate rotations.
During the first rotation, the entire shell is placed in a cylindrical coordinate system, with the 1, 2, and 3 axes representing the radial, circumferential, and axial directions, respectively. For the end-cap section shown in Figure 11, take the 2-axis as the rotation axis, rotate the 1–3 plane by an angle α , and the rotated plane is a 1′–3′ plane. α is the angle between the original radial axis 1 and the axis 1′ along the generatrix direction of the frustum of a cone.
The second rotation is a plane rotation, taking the 3′ axis as the rotation axis, rotating from the 1′–2 plane to the 1″–2″ plane, and the rotation angle θ is the fiber winding angle α ( r ) calculated by Equation (10). The rotation process is shown in Figure 12.
The winding shell is composed of composite materials, which are considered transversely isotropic materials. The material stiffness matrix is denoted as C , where each parameter represents the composite material properties calculated using micromechanical equations for composite materials, incorporating both carbon fibers and resin.
After the rotation, the constitutive equation incorporating angular variables is formed, with the stress–strain relationship as given in Equation (11). C represents the stiffness matrix of the composite material after coordinate transformation.
σ x σ y σ z τ x y τ x z τ y z = T σ 1 C T ε ε x ε y ε z γ x y γ x z γ y z = C ε x ε y ε z γ x y γ x z γ y z ,
The stress and strain transfer matrices are given as follows:
T σ = cos 2 θ sin 2 θ 2 sin θ cos θ 0 0 0 sin 2 θ cos 2 θ 0 2 sin θ cos θ 0 0 0 0 1 0 0 0 sin θ cos θ sin θ cos θ 0 cos 2 θ sin 2 θ 0 0 0 0 0 0 cos θ sin θ 0 0 0 0 sin θ cos θ ,
T ε = cos 2 θ sin 2 θ 0 sin θ cos θ 0 0 sin 2 θ cos 2 θ 0 sin θ cos θ 0 0 0 0 1 0 0 0 2 sin θ cos θ 2 sin θ cos θ 0 cos 2 θ sin 2 θ 0 0 0 0 0 0 cos θ sin θ 0 0 0 0 sin θ cos θ ,
On the other hand, the equation of the linear elastic constitutive equation is as follows:
σ i j = C i j k l ε k l ,
C = 1 ν 23 ν 32 E 2 E 3 Δ ν 12 + ν 13 ν 32 E 2 E 3 Δ ν 13 + ν 12 ν 23 E 2 E 3 Δ 0 0 0 ν 12 + ν 13 ν 32 E 2 E 3 Δ 1 ν 13 ν 31 E 1 E 3 Δ ν 23 + ν 21 ν 13 E 1 E 3 Δ 0 0 0 ν 13 + ν 12 ν 23 E 2 E 3 Δ ν 23 + ν 21 ν 13 E 1 E 3 Δ 1 ν 12 ν 21 E 1 E 2 Δ 0 0 0 0 0 0 G 12 0 0 0 0 0 0 G 13 0 0 0 0 0 0 G 23 ,
where
Δ = 1 ν 12 ν 21 ν 23 ν 32 ν 13 ν 31 + 2 ν 21 ν 32 ν 13 E 1 E 2 E 3 ,
After coordinate conversion, the constitutive equation including temperature and cure variables can be written as:
σ i j ( t ) = t C i j k l ( α , T , t τ ) t ε k l ( τ ) ε k l t c ( τ ) d τ ( i , j , k , l = 1 , 2 , 3 ) ,

3.3. Flowchart of the Curing Coupled Model

The curing analysis flow chart of thermo–chemo–mechanical coupled composite materials is shown in Figure 13, which is divided into curing dynamics module, heat conduction module and curing deformation module. Among them, the curing dynamics module and heat conduction module are coupled with each other, and the two modules are sequentially coupled with the curing deformation module. The curing calculation of thermo–chemo–mechanical coupling requires the writing of ABAQUS user subroutines. Different modules require different subroutines. The curing kinetics module uses UMAT definitions, the heat transfer module uses DISP, UEXPAN, and UMATHT definitions, and the curing deformation module uses UMAT definitions.

3.4. FE Model of the Composite Shell

A thermo–chemo–mechanical coupled analysis is performed on the whole structure of the carbon fiber composite material shell. The shell consists of carbon fiber winding layers, metal fittings, rubber insulation layers and a sand core mold, as illustrated in Figure 14. In the figure, the blue part represents the metal fittings, the cyan part denotes the rubber insulation layer, the gray part is the sand core mold, and the yellow part is the composite material layer.
The specific size of the filament wound composite shell is shown in Figure 15. The two ends of the shell have symmetrical polar holes, and the size is the same, the thickness of the rubber insulation layer is 2 mm, the barrel body is 146.5 mm, the length is 400 mm, the distance from the head to the polar hole is 98 mm, the radius of the metal polar hole is 50 mm, and the thickness of the single layer of the fiber winding layer is 0.2 mm. The model was established in the finite element software ABAQUS 2017, and the fiber winding method was [±19.5°/90° ± 19.5°]. The coordinate conversion of material parameters and the calculation of the curing process were carried out by the user subroutine written by ABAQUS.
Apart from resin and carbon fibers, the model incorporates three other materials: metal fittings, the sand core mold, and the rubber insulation layer. The insulation layer is made of ethylene propylene diene monomer (EPDM) rubber, and the front and rear metal fittings are made of 30CrMnSiA material. The mechanical properties of each material are presented in Table 2, while the elastic mechanical parameters of T800 carbon fibers are shown in Table 3.
The curing process of the composite shell is shown in Figure 16, which is applied to the outer surface of the shell model. The contact type between the fiber winding layer and the insulation layer is the Tie constraint. Considering the symmetry properties of the shell model [37], the 1/36 model is adopted for calculation. The symmetric boundary conditions along the column coordinates are set on both sides of the cross-section to prevent rigid body displacement.

4. Results and Discussions

4.1. Temperature and Curing Degree Analysis

A, B, and C are selected as the observation points, which are distributed from the outside to the inside near the polar hole, as shown in Figure 17.
Figure 18a presents the temperature comparison, and Figure 18b compares the degree of cure and cure rate. In the first heating stage, the outer surface is heated, and temperature conducts from the outside inwards, resulting in internal temperatures being lower than external ones. The increase in temperature initiates resin curing, with different temperatures leading to a gradient distribution in the degree of cure, and external curing reactions occurring faster than internal ones. When the heating of the first stage is nearly complete and enters the subsequent heat preservation stage, the temperature is conducted to the interior. The resin begins to undergo intense curing reactions, with the internal degree of cure increasing. During this process, a significant amount of heat is released. In Figure 18a, it can be seen that points B and C have significantly higher temperatures than point A, which is similar to the results in the literature [38]. In addition, it is evident that the concentration of heat causes the fastest cure rate at point C, as observed in Figure 18b. However, due to the thin fiber winding layer of the shell, the phenomenon of heat concentration is not obvious. In the later stages of heating and holding—that is, after about 650 min—no significant accumulation of heat occurs as the resin cure reaction weakens. Hence, the temperature lags behind the curing process. During the cooling stage, internal temperatures are higher than external, conducting heat from the inside outwards, with temperatures gradually approaching the cure temperature.

4.2. Residual Stress and Residual Strain Analysis

During the curing process, points A, B, and C distributed from the outside to the inside of the cylindrical section of the shell, as shown in Figure 19, are selected to plot the trends of axial residual stress and strain, as illustrated in Figure 20a,b. It can be found that, in the first heating stage, the strain primarily increases due to thermal expansion. At the holding stage, as the resin curing progresses, the strain decreases due to the cure shrinkage effect. When the resin curing is near completion, the strain change is positively correlated with the temperature change, which is similar to the strain change in the literature [39]. When comparing the strain at internal and external points, due to the shell being thin, their strains are essentially the same.
For the stress situation, the overall trend of stress during the curing process increases. In each heating stage, stress decreases when the temperature rises. In the holding stage before the end of curing, the cure shrinkage of the resin causes a significant increase in stress, indicating that cure shrinkage is a significant factor in the occurrence of residual stress. Comparing the stress at internal and external points, since the strain is similar, and points A and C are in the helical winding layer with the same winding angle, the stress changes are essentially the same. However, point B is in the circumferential winding layer, and after coordinate transformation, the constitutive equation differs from the other two points, resulting in stress being greater at point B than at points A and C.

5. Conclusions

In this paper, the results of the curing deformation of filament wound composite shells are studied using experimental and numerical studies. A thermo–chemo–mechanical coupled model for the curing analysis of the composite shell is proposed, which can be used to simulate the curing process and evaluate the quality of the composite shell. The detailed conclusions are listed as follows:
  • The non-isothermal DSC test was carried out on HCM-2184 resin, and the curing reaction was analyzed to belong to the autocatalytic reaction model. The parameters of the curing kinetic formula were obtained by mathematical fitting. The elastic modulus of the resin after curing was obtained using a tensile test. These material parameters are used in the curing numerical simulation.
  • Combining the heat conduction, the curing behavior of resin, and the curing deformation of composite materials, a thermo–chemo–mechanical coupled numerical model was constructed on the basis of considering the coordinate transformation of the composite shell, and the curing process of the composite shell model was simulated.
  • The curing process of the composite shell is analyzed. The results show that the heat accumulation caused by curing heat release can be observed in the thicker parts of the shell, such as the head section. However, due to the thin shell, the curing heat generation is not obvious. During the curing process, the strain changes inside and outside the shell barrel body are close, and the stress is different due to the change in fiber winding angle. During the intensive curing stage of the resin, the strain decreases rapidly and the stress increases, indicating that the curing shrinkage behavior of the resin is an important factor in the generation of residual stress and residual strain.

Author Contributions

Conceptualization: M.Z.; methodology, M.L. and T.S.; software, L.L.; validation, M.L. and L.L.; formal analysis, L.L.; investigation, M.L. and L.L.; resources, J.D.; data curation, L.L., S.H. and Y.T.; writing—original draft preparation, L.L.; writing—review and editing, M.L. and M.Z.; visualization, L.L.; supervision, M.Z.; project administration, M.L., J.H. and M.Z.; funding acquisition, J.D. and M.Z. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the National Natural Science Foundation of China (Grant No. 11972199), the Key Program of the Natural Science Foundation of Zhejiang Province (Grant No. LZ22A020001), and the Ningbo Major Research and Development Plan Project (Grant No. 2022Z210).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

The raw/processed data required to reproduce these findings cannot be shared at this time as the data also form part of an ongoing study.

Conflicts of Interest

The authors declare no conflicts of interest.

References

  1. Dong, C. Process-induced deformation of composite T-stiffener structures. Compos. Struct. 2010, 92, 1614–1619. [Google Scholar] [CrossRef]
  2. Li, J.; Yao, X.; Liu, Y.; Cen, Z.; Kou, Z.; Dai, D. A study of the integrated composite material structures under different fabrication processing. Compos. Part A 2009, 40, 455–462. [Google Scholar] [CrossRef]
  3. Mangalgiri, P.D. Composite materials for aerospace applications. Bull. Mater. Sci. 1999, 22, 657–664. [Google Scholar] [CrossRef]
  4. Neto, E.S.B.; Chludzinski, M.; Roese, P.B.; Fonseca, J.S.O.; Amico, S.C.; Ferreira, C.A. Experimental and numerical analysis of a LLDPE/HDPE liner for a composite pressure vessel. Polym. Test. 2011, 30, 693–700. [Google Scholar] [CrossRef]
  5. Mlýnek, J.; Rahimian Koloor, S.S.; Martinec, T.; Petrů, M. Fabrication of high-quality straight-line polymer composite frame with different radius parts using fiber winding process. Polymers 2021, 13, 497. [Google Scholar] [CrossRef] [PubMed]
  6. Zhi, J.; Yang, B.; Li, Y.; Tay, T.E.; Tan, V.B.C. Multiscale thermo-mechanical analysis of cure-induced deformation in composite laminates using Direct FE2. Compos. Part A 2023, 173, 107704. [Google Scholar] [CrossRef]
  7. Khoun, L.; Centea, T.; Hubert, P. Characterization methodology of thermoset resins for the processing of composite materials—Case study: CYCOM 890RTM epoxy resin. J. Compos. Mater. 2010, 44, 1397–1415. [Google Scholar] [CrossRef]
  8. Ersoy, N.; Garstka, T.; Potter, K.; Wisnom, M.R.; Porter, D.; Stringer, G. Modelling of the spring-in phenomenon in curved parts made of a thermosetting composite. Compos. Part A 2010, 41, 410–418. [Google Scholar] [CrossRef]
  9. Wisnom, M.R.; Potter, K.D.; Ersoy, N. Shear-lag analysis of the effect of thickness on spring-in of curved composites. J. Compos. Mater. 2007, 41, 1311–1324. [Google Scholar] [CrossRef]
  10. Yuan, Z.; Yang, G.; Yang, Z.; Feng, Y.; Li, S.; Li, Y.; Song, D. Process-Induced Deformation of L-Shaped Laminates: Analysis of Tool–Part Interaction. Mech. Compos. Mater. 2021, 56, 789–804. [Google Scholar] [CrossRef]
  11. Yoo, S.H.; Han, M.G.; Hong, J.H.; Chang, S.H. Simulation of curing process of carbon/epoxy composite during autoclave degassing molding by considering phase changes of epoxy resin. Compos. Part B 2015, 77, 257–267. [Google Scholar] [CrossRef]
  12. Behzad, T.; Sain, M. Finite element modeling of polymer curing in natural fiber reinforced composites. Compos. Sci. Technol. 2007, 67, 1666–1673. [Google Scholar] [CrossRef]
  13. Brauner, C.; Soprano, P.; Herrmann, A.S.; Meiners, D. Cure-dependent thermo-chemical modelling and analysis of the manufacturing process of an aircraft composite frame. J. Compos. Mater. 2015, 49, 921–938. [Google Scholar] [CrossRef]
  14. Qiao, M.; You, B.; Xu, J.; Wang, X. Analysis of the influence of the mold on a thick thermoset composite flange during curing. Int. J. Adv. Manuf. Tech. 2015, 78, 603–612. [Google Scholar] [CrossRef]
  15. Vafayan, M.; Ghoreishy, M.H.R.; Abedini, H.; Beheshty, M.H. Development of an optimized thermal cure cycle for a complex-shape composite part using a coupled finite element/genetic algorithm technique. Iran. Polym. J. 2015, 24, 459–469. [Google Scholar] [CrossRef]
  16. Lee, W.I.; Loos, A.C.; Springer, G.S. Heat of reaction, degree of cure, and viscosity of Hercules 3501-6 resin. J. Compos. Mater. 1982, 16, 510–520. [Google Scholar] [CrossRef]
  17. Zanjanijam, A.R.; Wang, X.; Ramezani, M.; Holberg, S.; Johnson, P.A. Phenolic resin/coal char composites: Curing kinetics and thermal/mechanical performance. Polymer 2023, 281, 126103. [Google Scholar] [CrossRef]
  18. Struzziero, G.; Remy, B.; Skordos, A.A. Measurement of thermal conductivity of epoxy resins during cure. J. Appl. Polym. Sci. 2019, 136, 47015. [Google Scholar] [CrossRef]
  19. Chern, B.C.; Moon, T.J.; Howell, J.R.; Tan, W. New experimental data for enthalpy of reaction and temperature-and degree-of-cure-dependent specific heat and thermal conductivity of the Hercules 3501-6 epoxy system. J. Compos. Mater. 2002, 36, 2061–2072. [Google Scholar] [CrossRef]
  20. Friis-Pedersen, H.H.; Pedersen, J.H.; Haussler, L.; Storm, B.K. Online measurement of thermal diffusivity during cure of an epoxy composite. Polym. Test. 2006, 25, 1059–1068. [Google Scholar] [CrossRef]
  21. Johnston, A.A. An Integrated Model of the Development of Process-Induced Deformation in Autoclave Processing of Composite Structures. Ph.D. Thesis, University of British Columbia, Vancouver, BC, Canada, 1997. [Google Scholar]
  22. Pantelelis, N.; Vrouvakis, T.; Spentzas, K. Cure cycle design for composite materials using computer simulation and optimisation tools. Forsch. Ingenieurwes. 2003, 67, 254–262. [Google Scholar] [CrossRef]
  23. Loos, A.C.; Springer, G.S. Curing of epoxy matrix composites. J. Compos. Mater. 1983, 17, 135–169. [Google Scholar] [CrossRef]
  24. Yi, S.; Hilton, H.H.; Ahmad, M.F. A finite element approach for cure simulation of thermosetting matrix composites. Comput. Struct. 1997, 64, 383–388. [Google Scholar] [CrossRef]
  25. Cheung, A.; Yu, Y.; Pochiraju, K. Three-dimensional finite element simulation of curing of polymer composites. Finite Elem Anal Ded. 2004, 40, 895–912. [Google Scholar] [CrossRef]
  26. Abdelal, G.F.; Robotham, A.; Cantwell, W. Autoclave cure simulation of composite structures applying implicit and explicit FE techniques. Int. J. Mech. Mater. Des. 2013, 9, 55–63. [Google Scholar] [CrossRef]
  27. Bogetti, T.A.; Gillespie, J.W., Jr. Process-induced stress and deformation in thick-section thermoset composite laminates. J. Compos. Mater. 1992, 26, 626–660. [Google Scholar] [CrossRef]
  28. Zhang, J.T.; Zhang, M.; Li, S.X.; Pavier, M.J.; Smith, D.J. Residual stresses created during curing of a polymer matrix composite using a viscoelastic model. Compos. Sci. Technol. 2016, 130, 20–27. [Google Scholar] [CrossRef]
  29. Abouhamzeh, M.; Sinke, J.; Jansen, K.M.B.; Benedictus, R. Thermo-viscoelastic analysis of GLARE. Compos. Part B 2016, 99, 1–8. [Google Scholar] [CrossRef]
  30. Zocher, M.A.; Groves, S.E.; Allen, D.H. A three-dimensional finite element formulation for thermovis-coelastic orthotropic media. Int. J. Numer. Meth. Eng. 1997, 40, 2267–2288. [Google Scholar] [CrossRef]
  31. Ding, A.; Li, S.; Wang, J.; Zu, L. A three-dimensional thermo-viscoelastic analysis of process-induced residual stress in composite laminates. Compos. Struct. 2015, 129, 60–69. [Google Scholar] [CrossRef]
  32. Ding, A.; Li, S.; Sun, J.; Wang, J.; Zu, L. A thermo-viscoelastic model of process-induced residual stresses in composite structures with considering thermal dependence. Compos. Struct. 2016, 136, 34–43. [Google Scholar] [CrossRef]
  33. GB/T 2567-2021; Test Methods for Properties of Resin Casting Body. National Standards of People’s Republic of China: Beijing, China, 2021.
  34. Fatmi, M.; Ghebouli, B.; Ghebouli, M.A.; Chihi, T.; Hafiz, M.A. The kinetics of precipitation in Al-2.4 wt% Cu alloy by Kissinger, Ozawa, Bosswel and Matusita methods. Phys. B 2011, 406, 2277–2280. [Google Scholar] [CrossRef]
  35. Wibowo, E.S.; Park, B.D. Cure kinetics of low-molar-ratio urea-formaldehyde resins reinforced with modified nanoclay using different kinetic analysis methods. Thermochim. Acta 2020, 686, 178552. [Google Scholar] [CrossRef]
  36. Hong, L.C.; Hwang, S.J. Study of warpage due to P-V-T-C relation of EMC in IC packaging. IEEE Trans. Compon. Packag. Technol. 2004, 27, 291–295. [Google Scholar] [CrossRef]
  37. Wu, Q.; Zu, L.; Wang, P.; Liu, Z.; Zhang, Q.; Zhang, G.; Li, D. Design and fabrication of carbon-fiber-wound composite pressure vessel with HDPE liner. Int. J. Pressure Vessels Piping 2022, 200, 104851. [Google Scholar] [CrossRef]
  38. Kim, Y.K.; White, S.R. Viscoelastic analysis of processing-induced residual stresses in thick composite laminates. Mech. Compos. Mater. Struct. 1997, 4, 361–387. [Google Scholar] [CrossRef]
  39. Olivier, P.A. A note upon the development of residual curing strains in carbon/epoxy laminates. Study by thermomechanical analysis. Compos. Part A 2006, 37, 602–616. [Google Scholar] [CrossRef]
Figure 1. Tensile sample size (unit: mm).
Figure 1. Tensile sample size (unit: mm).
Polymers 16 01643 g001
Figure 2. Results of non-isothermal DSC testing of HCM-2184 resin.
Figure 2. Results of non-isothermal DSC testing of HCM-2184 resin.
Polymers 16 01643 g002
Figure 3. Fitting of activation energy: (a) the Kissinger equation; (b) the Ozawa equation.
Figure 3. Fitting of activation energy: (a) the Kissinger equation; (b) the Ozawa equation.
Polymers 16 01643 g003
Figure 4. Changes in curing degree of the HCM-2184 resin at different heating rates.
Figure 4. Changes in curing degree of the HCM-2184 resin at different heating rates.
Polymers 16 01643 g004
Figure 5. Comparison of curing degree and curing rate at different heating rates: (a) 5 °C/min; (b) 10 °C/min; (c) 15 °C/min; (d) 20 °C/min.
Figure 5. Comparison of curing degree and curing rate at different heating rates: (a) 5 °C/min; (b) 10 °C/min; (c) 15 °C/min; (d) 20 °C/min.
Polymers 16 01643 g005
Figure 6. Relationship between the curing reaction rate and curing degree of the HCM-2184 resin at different heating rates.
Figure 6. Relationship between the curing reaction rate and curing degree of the HCM-2184 resin at different heating rates.
Polymers 16 01643 g006
Figure 7. Curve fitting results of the curing rate with different heating rate: (a) 5 °C/min; (b) 10 °C/min; (c) 15 °C/min; (d) 20 °C/min.
Figure 7. Curve fitting results of the curing rate with different heating rate: (a) 5 °C/min; (b) 10 °C/min; (c) 15 °C/min; (d) 20 °C/min.
Polymers 16 01643 g007
Figure 8. Average result of the curing rate fitting curve.
Figure 8. Average result of the curing rate fitting curve.
Polymers 16 01643 g008
Figure 9. (a) Force-displacement curve; (b) stress–strain curve and fitting results.
Figure 9. (a) Force-displacement curve; (b) stress–strain curve and fitting results.
Polymers 16 01643 g009
Figure 10. (a) Winding angle distribution of the spiral layer; (b) the winding angle of the spiral layer changes along the shell axial direction.
Figure 10. (a) Winding angle distribution of the spiral layer; (b) the winding angle of the spiral layer changes along the shell axial direction.
Polymers 16 01643 g010
Figure 11. Schematic of the first rotation.
Figure 11. Schematic of the first rotation.
Polymers 16 01643 g011
Figure 12. Schematic of the second rotation.
Figure 12. Schematic of the second rotation.
Polymers 16 01643 g012
Figure 13. Thermo–chemo–mechanical coupled curing flow chart.
Figure 13. Thermo–chemo–mechanical coupled curing flow chart.
Polymers 16 01643 g013
Figure 14. The model diagram of the filament wound composite shell.
Figure 14. The model diagram of the filament wound composite shell.
Polymers 16 01643 g014
Figure 15. Specific dimensions of composite gas cylinders.
Figure 15. Specific dimensions of composite gas cylinders.
Polymers 16 01643 g015
Figure 16. Curing process curve.
Figure 16. Curing process curve.
Polymers 16 01643 g016
Figure 17. The position of three points A, B, and C of the shell head segment.
Figure 17. The position of three points A, B, and C of the shell head segment.
Polymers 16 01643 g017
Figure 18. (a) Curve of temperature change; (b) curve of curing degree and curing rate.
Figure 18. (a) Curve of temperature change; (b) curve of curing degree and curing rate.
Polymers 16 01643 g018
Figure 19. The position of three points A, B, and C of the shell cylinder body.
Figure 19. The position of three points A, B, and C of the shell cylinder body.
Polymers 16 01643 g019
Figure 20. (a) Axial stress changes during curing. (b) Axial strain changes during curing.
Figure 20. (a) Axial stress changes during curing. (b) Axial strain changes during curing.
Polymers 16 01643 g020
Table 1. Characteristic parameters of resin non-isothermal DSC test.
Table 1. Characteristic parameters of resin non-isothermal DSC test.
Heating Rate
(°C/min)
Peak Temperature
(°C)
Maximum Heat
(W/g)
Total Heat Release
(J/g)
5164.85 (0.212)2.891 (0.214)1962.50 (60.104)
10184.70 (0.566)5.604 (0.173)1965.00 (74.953)
15199.40 (0.849)8.191 (0.714)1938.00 (82.024)
20204.75 (0.424)9.263 (0.651)1763.00 (67.882)
Table 2. Mechanical property parameters of shell material.
Table 2. Mechanical property parameters of shell material.
Parameter30CrMnSiA SteelCore MoldEPDM RubberUnit
E 2 × 1051 × 1067.8MPa
ν 0.330.10.47-
k 29.30.160.26w·m−1·K−1
α 1.172 × 10−51 × 10−71 × 10−5K−1
C P 52010002200J·kg−1·K−1
Table 3. Mechanics property of T800 carbon fiber.
Table 3. Mechanics property of T800 carbon fiber.
ParameterUnitValue
E 1 MPa294,000
E 2 , E 3 MPa14,000
G 12 , G 13 MPa15,000
G 23 MPa5500
ν 12 = ν 13 -0.23
ν 23 -0.23
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Lu, L.; Huan, S.; Lu, M.; Shen, T.; Tian, Y.; Hu, J.; Du, J.; Zhang, M. Three-Dimensional Thermo–Chemo–Mechanical Coupled Curing Analysis for the Filament Wound Composite Shell. Polymers 2024, 16, 1643. https://doi.org/10.3390/polym16121643

AMA Style

Lu L, Huan S, Lu M, Shen T, Tian Y, Hu J, Du J, Zhang M. Three-Dimensional Thermo–Chemo–Mechanical Coupled Curing Analysis for the Filament Wound Composite Shell. Polymers. 2024; 16(12):1643. https://doi.org/10.3390/polym16121643

Chicago/Turabian Style

Lu, Linjiao, Shengsheng Huan, Mengkai Lu, Tao Shen, Yanhui Tian, Jianying Hu, Jianke Du, and Minghua Zhang. 2024. "Three-Dimensional Thermo–Chemo–Mechanical Coupled Curing Analysis for the Filament Wound Composite Shell" Polymers 16, no. 12: 1643. https://doi.org/10.3390/polym16121643

APA Style

Lu, L., Huan, S., Lu, M., Shen, T., Tian, Y., Hu, J., Du, J., & Zhang, M. (2024). Three-Dimensional Thermo–Chemo–Mechanical Coupled Curing Analysis for the Filament Wound Composite Shell. Polymers, 16(12), 1643. https://doi.org/10.3390/polym16121643

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