This study examined the influence mechanism of temperature on the interfacial shear strength (IFSS) between carbon fiber (CF) and epoxy resin (EP) matrices under various thermal loads using experimental and numerical simulation methods. To evaluate the change in IFSS as a function of the increase in temperature, a microbond test was performed under controlled temperature environment from 23°C to 150°C. The experimental results showed that IFSS values of CF/EP reduce significantly when the temperature reaches near glass transition temperature. To interpret the effect of thermal loads on IFSS, a thermal-mechanical coupling finite element model was used to simulate the process of fiber pull-out from EP. The results revealed that temperature dependence of IFSS is linked to modulus of the matrix as well as to the coefficients of thermal expansion of the fiber and matrix.
Fundamental Research Funds for the Central Universitiesxjj2017601. Introduction
The interface between fiber and matrix plays an important role in fiber-reinforced composite materials by affecting directly the load transfer efficiency between the fiber and the matrix as well as the mechanical properties of the composites [1, 2]. The interfacial shear strength (IFSS) is a key parameter that influences the adhesion performance between the fiber and matrix [3, 4]. However, because of the heat instability of EP [5, 6] and the significant difference between the fiber and the EP matrix [7], IFSS becomes unstable as it is affected by the curing process as well as service temperature [8].
The coefficients of thermal expansion (CLTE) of EP matrices are almost 100 times higher than the longitudinal thermal expansion coefficient of CF and 10 times greater than the transverse thermal expansion coefficient of CF by absolute value. These differences cause residual thermal stress at the interface during the cooling process [9–11]. As a result, tremendous efforts have been devoted to overcoming these problems. For example, Sockalingam et al. [12] developed a finite element modeling methodology to simulate the residual thermal stress with large CLTE mismatch between the fiber and the matrix caused by the cooling process. Also, Di Landro and Pegoraro [13] investigated the influence of the processing temperature and cooling rates on the residual stress in polymer composite matrices. They proposed that the residual thermal stress is mainly caused by the mismatch of CLTE between the fiber and matrix that importantly contributes to the fiber-matrix adhesion. Choo et al. [14] studied thermal path-dependency in the stress evolution of continuous-fiber-reinforced composites by comparing cooling and reheat cycles. Their results showed that the effect of time-dependent deformation becomes more significant as the fiber content increases and cooling/heating rate reduces. However, most of the above-mentioned studies mainly focus on the influence of the curing process of IFSS; only a handful of studies were devoted to investigating cured CF/EP service under elevated temperature work conditions.
In fact, IFSS also changed with serviced temperature. For example, Thomason and Yang [15] explored the temperature dependence of IFSS in glass fiber polypropylene composites (GF-PP) and noticed that IFSS in GF-PP is inversely dependent on the testing temperature with a major increase in glass transition region of the PP matrix. In their other paper [16], they found that the magnitude of the thermal residual stress due to mismatch in the thermal expansion coefficients of fiber and matrix was insufficient to explain the magnitude of the system IFSS, and they suggest that the residual stress combined with static adhesion could be the major contributor to the apparent interfacial adhesion in glass fiber-epoxy system, but they did not consider the effects of the temperature of resin on IFSS. In order to develop improved process ability and equipment operating at high temperature, Wang et al. [17] studied the thermal stability of carbon nanotube/epoxy nanocomposite by mechanical and interfacial evaluation at room temperature and 150°C. The results show that addition of low concentration CNT could improve the thermal stability of epoxy resin in turn hence the IFSS under high temperature. Also, the similar conclusions had been reached in similar studies [18]. So, when cured CF/EP service is under changing temperature work conditions, both the mismatch of CLTE between the fiber and matrix and the dependence-temperature of resin will cause a change in IFSS, hence affecting the CFRP performance. As a result, studying the mechanism of the influence of temperature on the interfacial properties of cured composites is of great importance.
The aim of this paper is to investigate the mechanism due to the influence of temperature on the interfacial properties of cured CFRP. In order to get the temperature dependence of IFSS, a microbond test was performed under controlled temperature; thermomechanical property tests were performed to get the Tg and the CLTE of resin matrix; tensile test was performed to measure the value of modulus changed with temperature. Furthermore, a FEM was developed to simulate the process of the microbond test under thermal loads to facilitate interpretation of the experiment observations. The FE analysis was carried out following three steps: (i) application of thermal preload and large mismatch in CTE between CF and EP to analyze the residual thermal stress due to the cooling process, (ii) modeling of the thermal stress in CF/EP under thermal preload through changes from room to the working temperature, and (iii) modeling the process of single fiber pull-out from the microdroplet using the prestress of both residual thermal and thermal stresses.
2. Experimental2.1. Materials and Test Sample
The EPs named TF1408 were obtained, respectively, from Hengshen Co. Ltd. Carbon fibers T300 were purchased from Toray Industries, Inc. The single CF was bonded on a thin stalloy adhered through microdroplet uncured EPs (Figure 1). All the samples were cured at 135°C and kept for 90 minutes in an oven, with the same heating rate and cooling rate of 3°C/min.
Temperature dependence HM410-microdroplet test equipment.
2.2. Microbond Testing
The temperature dependence of CF/EP IFSS was measured by HM410-microdroplet test equipment from TOHEI SANGYO, Japan (Figure 1). The temperature was automatically controlled and monitored through a digital temperature programmer. Five tests were performed at room temperature (RT), 60, 90, and 120°C.
The required force to pull out the fiber of cured epoxy resin is measured for each specimen and the interfacial shear strength τIFSS was estimated by [19, 20](1)τIFSS=Fmaxπ·df·le,where Fmax is the maximum value of pull-out force, df is the diameter of single fiber, and le is the embedded length of microdroplet. However, since at the moment when F=Fmax the crack length is not zero and the measured Fmax value includes the contributions of interfacial adhesion and friction, the value of τIFSS calculated by (1) is the apparent IFSS.
The diameter of CF utilized in this study was 7 μm. For a big effect of specimen dimensions on the pull-out force and IFSS, the embedded length of EP microdroplet was maintained at 70–90 μm and the contact angle at 40–45°C (Figure 1). The load-displacement curve from each test was recorded to obtain the maximum force (Fmax) and a typical example is shown in Figure 2, where Fd is kink force [21]; when the pull-out force reaches this value, the fiber begins to debond from the matrix; Fmax is the peak force caused by debonding force and the frictional load in debonding regions; Fb is the tail force caused by frictional load when the embedded fiber totally debond from the matrix.
Load-displacement curve for the microdroplet test.
2.3. Thermomechanical Property Tests
The thermomechanical properties of EPs were tested by dynamic mechanical analysis (DMA) and thermal expansion analysis [22]. A dynamic mechanical analyzer (DMA Q800) was used for measurement of the glass transition temperature (Tg) as well as modulus versus temperature curve of both cured EPs. DMA tests were conducted using three-point bending modes (ASTM D5023-15) as it can provide the most analogous approach to normal mechanical testing method [23]. According to ASTM D5023-15, the dimensions utilized for DMA tests were (60 × 13 × 3 mm) and temperature scanning measurements were performed from 20°C to 150°C at a heating rate of 3°C/min and an oscillation frequency of 1 Hz. To ensure that measured temperature range includes the temperature of Tg, the real temperature sweep test was performed from 0°C to 200°C. According to ASTM E289, the coefficient of the linear thermal expansion (CLTE) of EPs was measured using a dilatometer (DIL402C) from 20°C to 150°C at a heating rate 5°C/min for CLTE test specimens of (25 × 5 × 5 mm).
2.4. Tensile Testing of Epoxy Resin
The moduli changed with different temperature were tested by this test. The temperature dependence of resin under tensile loading was determined using a universal testing machine (CMT4304; MTS, USA) according to the ASTM D638-10 “Standard Test Method for Tensile Properties of Plastics.” The dumbbell-shaped samples (type I) were stretched using an extensometer at a rate of 2 mm/min and stress-strain curves were recorded for ambient temperatures ranging from RT to 120°C (at 30°C intervals) within the climatic chamber. For each condition, 5 tests were performed.
3. Numerical Simulations Procedures3.1. Finite Element Modeling
Commercial finite element software Abaqus@ (version 6.10) was utilized to carry out FE calculations. The baseline FE model of the microbond test consists of a single CF with a diameter of 7 μm and an epoxy matrix droplet with embedded length of 80 μm made by eight-node fully integrated solid elements. The matrix was modeled as an isotropic elastic material with properties shown in Table 1. Here, Young’s modulus of resin was the slope of stress-strain curves under linear-elastic range and measured by tensile tests (Figure 3); the thermal expansion of resin under different temperatures was measured by CLTE tests (Figure 4). The fiber was modeled as an orthotropic elastic material with properties shown in Table 2 [24]. The knife blade was modeled as a rigid body and owing to its symmetry, only a quarter of the fiber and droplet were modeled (Figure 5).
Matrix properties at room temperature used for the numerical simulations.
Test temperature
RT
60°C
90°C
120°C
Young’s modulus, Em (MPa)
3960
3068
2510
1590
Poisson’s ratio
0.37
0.37
0.37
0.37
Coefficient of thermal expansion (10-6/°C)
58
60
65
75
Carbon fiber properties at room temperature used for the numerical simulations.
Elastic properties
Carbon fiber
Longitudinal stiffness, E11(MPa)
233130
Transverse stiffness, E22 (MPa)
23110
Out-of-plane stiffness, E33 (MPa)
23110
Poisson’s ratio, υ12
0.2
Poisson’s ratio, υ13
0.2
Poisson’s ratio, υ23
0.4
Shear modulus, G12 (MPa)
8960
Shear modulus, G13 (MPa)
8960
Shear modulus, G23 (MPa)
8270
Coefficient of thermal expansion α1 (10-6/°C)
−0.54
Coefficient of thermal expansion α2 (10-6/°C)
10.10
Representative stress-strain curves of the epoxy resin under different temperature conditions.
Coefficient of thermal expansion of epoxy resin versus temperature.
Quarter symmetric FEM of the microdroplet test.
3.2. Interfacial Behavior
The fiber-matrix interphase is subjected to the cohesive force and the residual compressive radial force, which results in a sliding friction force in the process of interfacial debonding. So in this study, the interfacial was model with cohesive element and coulomb friction contact behavior.
The zero thickness cohesive elements with a mixed-mode traction-separation law [25] (Figure 6) was used to model the process of interface debonding. The failure model of traction-separation used in this study is quadratic stress failure criterion [26] and a Benzeggagh-Kenane (BK) criterion [27] to evaluate, respectively, both the initial damage and crack propagation. The quadratic stress failure criterion was expressed through (2)maxσn,0Nmax2+σtTmax2+σsSmax2=1,where σn, σs, and σt are the stress components of, respectively, pure normal direction, first shear direction, and second shear direction. Nmax, Smax, and Tmax are the interfacial strengths of, respectively, pure normal direction, first shear direction, and second shear direction. The stiffness coefficient of cohesive element is estimated using the peak stress and the maximum relative experimental displacement by the results of microbond test. The BK criterion for crack propagation can be evaluated through [27] (3)GIC+GIIC-GICGII+GIIIGI+GII+GIIIη=GC,where Gi and GiC(i=I,II,III) are the corresponding strain energy release rates and fracture toughness under, respectively, pure mode I, mode II, and mode III loadings. GC=GI+GII+GIII is the fracture toughness, and η is the BK exponent. Fracture toughness GiC is estimated using the kink force Fd issued from microdroplet test at stress-free temperature through [28](4)GiC=rC33s2Fdπr2+αA-αmΔT2C33s2,where Fd is the kink force measured from microbond test, r is the radius of the fiber, and αA and αm are, respectively, the axial coefficient for thermal expansion of the fiber and the matrix. ΔT is the difference between the test and stress-free temperatures and C33s is given by (5)C33s=121EA+V1V2Em,where EA and Em are, respectively, the fiber axial modulus and matrix tensile modulus. V1 and V2 are volume fractions of the fiber and the matrix, respectively. In this study, the traction law for mode I is assumed to be the same as modes II and III.
Mixed-mode traction-separation law.
A coulomb friction is used to model the friction caused by residual compressive radial force. The coefficient of friction during progressive debonding is estimated by the ratio of tail force Fb tested by the results of microbond test (Figure 2) and the radial force caused by the residual thermal stress.
3.3. Residual Thermal Stress Analysis
The residual stresses due to thermal shrinkage during cooling down from the cure temperature to room temperature were simulated by application of thermal preload to the FE model. A thermal preload was applied using Tg as stress-free temperature and the operating temperature as room temperature (23°C). Tg of the epoxy resin was commonly used to define the stress-free temperature in calculation of the residual stress [29]. The model of the knife blade was suppressed in this step and instead an axisymmetric loading was applied to the matrix droplet and fiber. The mixed-mode damage and crack initiation loads were set on the fiber-matrix interphase.
3.4. Thermal Stress Analysis
The thermal stress analysis due to thermal preload was simulated under a prestress of residual thermal stress. A thermal preload was set on the fiber and matrix from room to working temperatures and Young’s modulus of the matrix was modified to the corresponding temperature. The model of the knife blade was suppressed in this step too and axisymmetric loading was applied to the matrix droplet and the fiber. A mixed-mode damage and crack initiation loads were set on the fiber-matrix interphase.
3.5. Microdroplet Debonding Analysis
The process of single fiber pull-out from the microdroplet was modeled under the prestress of residual thermal and thermal stresses. The knife blade was activated and fixed boundary conditions were applied to the knife blade. A contact property was set between the knife blade and droplet and constant displacement loading was applied to the fiber along the axial direction. An axisymmetric loading was applied to the matrix droplet and fiber and mixed-mode damage and crack initiation loads were set on the fiber-matrix interphase.
4. Results and Discussion
Figure 7 shows the force-displacement curve of microbond test under different temperature; the results show that the peak forces are decreased with tested temperature increasing. Correspondently, the apparent IFSS also decreased with tested temperature increasing. Figure 8 illustrates results from the thermomechanical characterization of EPs obtained through DMA and apparent IFSS as a function of temperature. It can clearly be seen that the measured apparent IFSS for CF/EP is significantly affected by the thermal load, where apparent IFSS decreases substantially near Tg. Also, the change point of apparent IFSS decrease with temperature (about 120–130°C from Figure 8) is approximately equal to Tg of EPs (about 123–131°C from Figure 8).
Force-displacement curves under different tested temperature by experiment results.
EPs matrix storage/loss modulus and apparent IFSS versus temperature.
Figure 9 depicts results from the process of fiber pull-out performed at three different temperatures, which further confirmed the effect of temperature on the apparent IFSS. It is clear that failure modes of the fiber pull-out significantly change as temperature increases. At processing temperatures below Tg, increasing the pull-out force induces a crack that propagates along the interface between the fiber and the matrix, where the droplet will be in complete debonding with very little deformation (Figure 9(a)). When the processing temperature reaches Tg, the stiffness of EP reduces as the pull-out force increases. This induces the droplet to completely debond with larger deformation effect (Figure 9(b)). At processing temperature above Tg, modulus of EP becomes particularly weak due to the relaxation transition in the polymer. As the apparent IFSS is larger than EP, failure will be induced in EP (Figure 9(c)).
The process of fiber outperformed at different temperatures: (a) temperature < Tg; (b) temperature near to Tg; and (c) temperature > Tg.
However, the above experiment results just proved the effect of temperature on the fiber pull-out peak force and apparent IFSS, the effect mechanism of the thermal stress and modulus on the apparent IFSS cannot be observed from the experiment results. So, it is necessary to analyze this question by finite element analysis. For the different fracture toughness, peak stress, stiffness coefficient of cohesive element, and the different modulus and CLET of resin matrix under 23, 60, 90, and 120°C were predicted by models with the same traction laws as shown in Figure 10. The values of apparent IFSS calculated by the peak force of the experiment and simulations are shown in Figure 11. It is clear from Figures 8 and 10 that, although there is still error, the trend of force-displacement curves of FE models is in good agreement with the experimental results. And the value of apparent IFSS by FE analysis is in the error range of experimental data.
Simulated force-displacement curves under different temperature.
Comparison of apparent IFSS between the simulated and experimental data.
The results for the same fracture toughness, peak stress, stiffness coefficient of cohesive element with different modulus of EP, and thermal stress under 60 and 90°C are gathered in Figures 12(a) and 12(b). It can be observed from Figure 12 that the thermal stress caused by the mismatch of CLTE under higher temperature can increase the peak force; on the contrary, the peak forces are decreasing with the modulus of resin matrix under higher temperature. So, the temperature dependence of CF/EP IFSS was also linked to the thermal stress caused by the mismatch of CLTE and the modulus of resin matrix. Due to the mismatch of CLTE, the IFSS is increased with temperature increasing and is decreased in Young’s modulus with temperature increasing. For the trend of IFSS is decreased with the temperature increasing in this study, the effect of the thermal stress caused by the mismatch of CLTE is lower than the modulus of matrix with temperature increasing.
Simulated force-displacement curves for the effect of CLET and EP’s modulus (a) under 60°C and (b) under 90°C.
5. Conclusions
The effects of temperature on the interfacial properties of CF/EP were investigated through a microbond test under controlled temperature environment. The experimental results revealed that apparent IFSS for CF/EP at temperatures below that of glass transition of EP does not show a significant decrease. However, this becomes significantly reduced at higher temperatures than that of glass transition. To interpret the experimental results, a series of FE calculations were carried out to simulate the process of microbond test as a function of the thermal loads. The results show that the thermal stress caused by the mismatch of CLTE under higher temperature can increase the peak force and apparent IFSS; oo the contrary, the peak forces are decreasing with the modulus of resin matrix under higher temperature. The effect of the thermal stress caused by the mismatch of CLTE is lower than the modulus of matrix with temperature increasing.
Conflicts of Interest
The authors declare that there are no conflicts of interest regarding the publication of this paper.
Acknowledgments
The authors acknowledge the financial support by “the Fundamental Research Funds for the Central Universities” (no. xjj201760).
MintyR.ThomasonJ.PetersenH.2015LiH.WangY.ZhangC.ZhangB.Effects of thermal histories on interfacial properties of carbon fiber/polyamide 6 composites: Thickness, modulus, adhesion and shear strength20168531392-s2.0-8496086851410.1016/j.compositesa.2016.03.007BeggsK. M.ServinisL.GengenbachT. R.HusonM. G.FoxB. L.HendersonL. C.A systematic study of carbon fibre surface grafting via in situ diazonium generation for improved interfacial shear strength in epoxy matrix composites201511831382-s2.0-8493961506010.1016/j.compscitech.2015.08.001LiuB.LiuZ.WangX.ZhangG.LongS.YangJ.Interfacial shear strength of carbon fiber reinforced polyphenylene sulfide measured by the microbond test20133247247302-s2.0-8487656937310.1016/j.polymertesting.2013.03.020ZabihiO.KhayyamH.FoxB. L.NaebeM.Enhanced thermal stability and lifetime of epoxy nanocomposites using covalently functionalized clay: Experimental and modelling2015393226922782-s2.0-8492436332010.1039/c4nj01768d10.1039/C4NJ01768DXuG.-R.XuM.-J.LiB.Synthesis and characterization of a novel epoxy resin based on cyclotriphosphazene and its thermal degradation and flammability performance20141092402482-s2.0-8490731816710.1016/j.polymdegradstab.2014.07.020PanY.PelegriA. A.Influence of matrix plasticity and residual thermal stress on interfacial debonding of a single fiber composite2010511291422-s2.0-7795792293110.2140/jomms.2010.5.129GaoS.-L.KimJ.-K.Cooling rate influences in carbon fibre/PEEK composites. Part II: Interlaminar fracture toughness20013267637742-s2.0-003537242810.1016/S1359-835X(00)00188-3LitherlandA. E.Le'HT.Transcrystallized interphase in thermoplastic composites. II: Influence of interfacial stress, cooling rate, fibre properties and polymer molecular weight1992WangY.HahnT. H.AFM characterization of the interfacial properties of carbon fiber reinforced polymer composites subjected to hygrothermal treatments2007671921012-s2.0-3375073204210.1016/j.compscitech.2006.03.030ZhuM.LiM.WuQ.GuY.LiY.ZhangZ.Effect of processing temperature on the micro- and macro-interfacial properties of carbon fiber/epoxy composites20142154434532-s2.0-8489638910610.1080/15685543.2014.877270SockalingamS.DeyM.GillespieJ. W.Jr.KeefeM.Finite element analysis of the microdroplet test method using cohesive zone model of the fiber/matrix interface2014562392472-s2.0-8489652505310.1016/j.compositesa.2013.10.021Di LandroL.PegoraroM.Evaluation of residual stresses and adhesion in polymer composites19962798478532-s2.0-002971749510.1016/1359-835X(96)00046-2ChooH.BourkeM. A. M.DaymondM. R.A finite-element analysis of the inelastic relaxation of thermal residual stress in continuous-fiber-reinforced composites20016112175717722-s2.0-003545059110.1016/S0266-3538(01)00075-6ThomasonJ. L.YangL.Temperature dependence of the interfacial shear strength in glass-fibre polypropylene composites20117113160016052-s2.0-8005191713110.1016/j.compscitech.2011.07.006ThomasonJ. L.YangL.Temperature dependence of the interfacial shear strength in glass-fibre epoxy composites2014967122-s2.0-8489708972210.1016/j.compscitech.2014.03.009WangZ.-J.KwonD.-J.ChoiJ.-Y.ShinP.-S.YiJ.-W.ByunJ.-H.LeeH.-I.ParkJ.-K.DevriesK. L.ParkJ.-M.Inherent and interfacial evaluations of carbon nanotubes/epoxy composites and single carbon fiber at different temperatures2016911111182-s2.0-8495817426510.1016/j.compositesb.2015.12.036GuruK.MishraS. B.ShuklaK. K.Effect of temperature and functionalization on the interfacial properties of CNT reinforced nanocomposites201534959652-s2.0-8493770078110.1016/j.apsusc.2015.04.196ZhangF.-H.WangR.-G.HeX.-D.WangC.RenL.-N.Interfacial shearing strength and reinforcing mechanisms of an epoxy composite reinforced using a carbon nanotube/carbon fiber hybrid20094413357435772-s2.0-6764998293010.1007/s10853-009-3484-xPisanovaE.ZhandarovS.MäderE.AhmadI.YoungR. J.Three techniques of interfacial bond strength estimation from direct observation of crack initiation and propagation in polymer-fibre systems2001323-44354432-s2.0-003528156110.1016/S1359-835X(00)00054-3ZhandarovS.MäderE.Characterization of fiber/matrix interface strength: Applicability of different tests, approaches and parameters20056511491602-s2.0-1364427794710.1016/j.compscitech.2004.07.003WatanabeH.YamadaN.OkajiM.Linear thermal expansion coefficient of silicon from 293 to 1000 K20042512212362-s2.0-354305134310.1023/B:IJOT.0000022336.83719.43DengS.HouM.YeL.Temperature-dependent elastic moduli of epoxies measured by DMA and their correlations to mechanical testing data20072668038132-s2.0-3454769163510.1016/j.polymertesting.2007.05.003HyerM. W.WhiteS. R.1998McGraw-HillSystemsD.2011CamanhoP. P.DávilaC. G.De MouraM. F.Numerical simulation of mixed-mode progressive delamination in composite materials20033716141514382-s2.0-004252722310.1177/0021998303034505BenzeggaghM. L.KenaneM.Measurement of mixed-mode delamination fracture toughness of unidirectional glass/epoxy composites with mixed-mode bending apparatus19965644394492-s2.0-003011180910.1016/0266-3538(96)00005-XScheerR. J.NairnJ. A.A Comparison of Several Fracture Mechanics Methods for Measuring Interfacial Toughness with Microbond Tests1995531-245682-s2.0-000250063610.1080/00218469508014371JiaY.YanW.LiuH.-Y.Carbon fibre pullout under the influence of residual thermal stresses in polymer matrix composites20126279862-s2.0-8486169275310.1016/j.commatsci.2012.05.019