Composite Interfaces, Vol. 14, No 1, pp 49-62(2007) Alsoavailableonlinewww.brillnl/ci Influence of inhomogeneous interphase on thermal stresses in fiber-reinforced composites L.H. YOU* College of Mechanical Engineering, Chongqing University, Chongqing City 400044, China and NCCA, Bournemouth University, Talbot Campus, BH12 5BB, UK Received 14 September 2005: accepted 5 January 2006 Abstract-The effects of an inhomogeneous interphase on thermal response of fiber-reinforced composites are investigated in this paper. The work is based on a linear variation of Youngs modulus and thermal expansion coefficient of the interphase and the assumption of generalized plane strain. An accurate analytical approach is presented to determine thermal stresses in composites reinforced with isotropic fibers containing an inhomogeneous interphase. How the inhomogeneous interphase, and Youngs modulus and thermal expansion coefficient of the matrix affect thermal response of the composites is examined. It is found that the inhomogeneous interphase causes different stress distribution from that of the homogeneous interphase. Raising Youngs modulus and thermal expansion coefficient of the matrix obviously increases the maximum radial, circumferential and axial stresses in all the constituents of the composites Keywords: Thermal stresses: fiber-reinforced composites; varying interphase: Youngs modulus thermal expansion coefficient; generalized plane st 1 INTRODUCTION Due to different material properties of the constituents, different deformations and stresses will occur in these constituents of fiber-reinforced composites even ubjected to the same temperature loads. These thermal deformations and stresse may cause some undesirable outcomes such as fiber breaks, debonds and cracks at interface, and plastic deformations of interfacial layers and matrix [1-3]. Therefore thermal analysis is always an active topic in composite mechanics There have been quite a number of numerical and analytical approaches deal- ing with thermal analysis of fiber-reinforced composites with a homogeneous inter- ase E-mail: Lyou@bournemouth. ac uk
Composite Interfaces, Vol. 14, No. 1, pp. 49–62 (2007) VSP 2007. Also available online - www.brill.nl/ci Influence of inhomogeneous interphase on thermal stresses in fiber-reinforced composites L. H. YOU ∗ College of Mechanical Engineering, Chongqing University, Chongqing City, 400044, China and NCCA, Bournemouth University, Talbot Campus, BH12 5BB, UK Received 14 September 2005; accepted 5 January 2006 Abstract—The effects of an inhomogeneous interphase on thermal response of fiber-reinforced composites are investigated in this paper. The work is based on a linear variation of Young’s modulus and thermal expansion coefficient of the interphase and the assumption of generalized plane strain. An accurate analytical approach is presented to determine thermal stresses in composites reinforced with isotropic fibers containing an inhomogeneous interphase. How the inhomogeneous interphase, and Young’s modulus and thermal expansion coefficient of the matrix affect thermal response of the composites is examined. It is found that the inhomogeneous interphase causes different stress distribution from that of the homogeneous interphase. Raising Young’s modulus and thermal expansion coefficient of the matrix obviously increases the maximum radial, circumferential and axial stresses in all the constituents of the composites. Keywords: Thermal stresses; fiber-reinforced composites; varying interphase; Young’s modulus; thermal expansion coefficient; generalized plane strain. 1. INTRODUCTION Due to different material properties of the constituents, different deformations and stresses will occur in these constituents of fiber-reinforced composites even subjected to the same temperature loads. These thermal deformations and stresses may cause some undesirable outcomes such as fiber breaks, debonds and cracks at interface, and plastic deformations of interfacial layers and matrix [1–3]. Therefore, thermal analysis is always an active topic in composite mechanics. There have been quite a number of numerical and analytical approaches dealing with thermal analysis of fiber-reinforced composites with a homogeneous interphase. ∗E-mail: Lyou@bournemouth.ac.uk
LH. You Ishikawa et al. presented a collocation method to calculate thermal expansion coefficients of unidirectional composites [4]. Using a three-dimensional hexagonal array finite element model, Bowles and Tompkins determined the longitudinal and transverse thermal expansion coefficients of different graphite fiber-reinforced resin, metal, and ceramic matrix composites [5]. With plane stress, plane strain and generalized plane strain finite element models, Muller and Schmauder examined the residual stresses at the fiber-matrix interface of the composites with regular cubic and hexagonal fiber arrangements [6]. Treating the geometrical structure of fiber reinforced composites as perfectly periodic, uniformly spaced fiber arrangements in square and hexagonal cells, as well as different cells in which either 30 or 60 fibers are randomly placed in the ductile matrix, Nakamura and Suresh investigated the combined effects of thermal residual stresses and fiber spatial distribution on the deformation of a 6061 aluminum alloy by using the finite element code ABAQUS [7]. Employing a 3D prism for the hexagonal fiber pattern, Abedian and Szyszkowski discussed stress singularities on the free surface at the fiber/matrix nterface with the ANSYS finite element code [8]. Recently, Choo et al. carried out a finite element analysis of the inelastic relaxation of the time-and path-dependent thermal residual stresses in continuous fiber-reinforced composites [9]. Chen and Liu presented a boundary element method with thin-body capabilities to mod multiple cells of fiber-reinforced composites [10] Numerical methods are effective in coping with complicated situations of thermal analysis of the composites. However, they are expensive. Therefore, some analytical models were developed to simplify the analysis Treating continuous fiber-reinforced composites as a cylinder consisting of a fiber and a concentric binder shell, Hashin and Rosen predicted the elastic moduli of the composites with hexagonal and random arrays of fibers [ll]. Using the rule of mixtures to calculate the properties of the outer cylinder, Mikata and Taya gave a four concentric cylinder model which consists of innermost transversely isotropic fiber, middle isotropic coating and matrix, and outermost transversely isotropic composite, and determined the elastic stress field in coated continuous fiber-reinforced composites subjected to thermomechanical loading [12]. Using the composite cylinder assemblage model of Hashin [13] to determine the properties of the outer cylinder, Kuntz et al. [14] also investigated the effects of fiber coating on the residual stresses in the same composites as those investigated by Mikata and Taya [12]. Applying the assumption of the generalized plane strain, Avery nd Herakovich developed an analytical method to determine the deformations and stresses in the composites reinforced with circumferentially orthotropic, radially orthotropic and transversely isotropic fibers [15]. Considering some composites may contain a number of coatings, Warwick and Clyne extended the model of Mikata and Taya [12] to that consisting of an arbitrary number of component cylinders, and predicted the thermal stresses in SiC monofilament systems [16] Sutcu presented a simple recursive algorithm to determine the material properties of a uniaxially aligned composite which contains an arbitrary number of coatings
50 L. H. You Ishikawa et al. presented a collocation method to calculate thermal expansion coefficients of unidirectional composites [4]. Using a three-dimensional hexagonal array finite element model, Bowles and Tompkins determined the longitudinal and transverse thermal expansion coefficients of different graphite fiber-reinforced resin, metal, and ceramic matrix composites [5]. With plane stress, plane strain and generalized plane strain finite element models, Müller and Schmauder examined the residual stresses at the fiber–matrix interface of the composites with regular cubic and hexagonal fiber arrangements [6]. Treating the geometrical structure of fiberreinforced composites as perfectly periodic, uniformly spaced fiber arrangements in square and hexagonal cells, as well as different cells in which either 30 or 60 fibers are randomly placed in the ductile matrix, Nakamura and Suresh investigated the combined effects of thermal residual stresses and fiber spatial distribution on the deformation of a 6061 aluminum alloy by using the finite element code ABAQUS [7]. Employing a 3D prism for the hexagonal fiber pattern, Abedian and Szyszkowski discussed stress singularities on the free surface at the fiber/matrix interface with the ANSYS finite element code [8]. Recently, Choo et al. carried out a finite element analysis of the inelastic relaxation of the time- and path-dependent thermal residual stresses in continuous fiber-reinforced composites [9]. Chen and Liu presented a boundary element method with thin-body capabilities to model multiple cells of fiber-reinforced composites [10]. Numerical methods are effective in coping with complicated situations of thermal analysis of the composites. However, they are expensive. Therefore, some analytical models were developed to simplify the analysis. Treating continuous fiber-reinforced composites as a cylinder consisting of a fiber and a concentric binder shell, Hashin and Rosen predicted the elastic moduli of the composites with hexagonal and random arrays of fibers [11]. Using the rule of mixtures to calculate the properties of the outer cylinder, Mikata and Taya gave a four concentric cylinder model which consists of innermost transversely isotropic fiber, middle isotropic coating and matrix, and outermost transversely isotropic composite, and determined the elastic stress field in coated continuous fiber-reinforced composites subjected to thermomechanical loading [12]. Using the composite cylinder assemblage model of Hashin [13] to determine the properties of the outer cylinder, Kuntz et al. [14] also investigated the effects of fiber coating on the residual stresses in the same composites as those investigated by Mikata and Taya [12]. Applying the assumption of the generalized plane strain, Avery and Herakovich developed an analytical method to determine the deformations and stresses in the composites reinforced with circumferentially orthotropic, radially orthotropic and transversely isotropic fibers [15]. Considering some composites may contain a number of coatings, Warwick and Clyne extended the model of Mikata and Taya [12] to that consisting of an arbitrary number of component cylinders, and predicted the thermal stresses in SiC monofilament systems [16]. Sutcu presented a simple recursive algorithm to determine the material properties of a uniaxially aligned composite which contains an arbitrary number of coatings
Influence of inhomogeneous interphase on composites on its fibers [17]. You gave an efficient method to tackle thermal analysis of composites containing anisotropic continuous fibers and an arbitrary number of interfacial layers [18] The above analytical research studies did not consider the plasticity of composites When thermal deformations of composites are large enough, some components of composites may become plastic. In order to address this problem, Zhang et al. used a two concentric cylinder model together with the increment theory of plasticity to derive two partial differential equations governing the deformation behavior of composites and numerically solved these two partial differential equations [19] Adopting Tresca's yield criterion and linear strain-hardening material model, You and Long discussed the effects of material properties of interfacial layer on thermal response of fiber-reinforced composites [20]. This work was further extended to the composites reinforced with anisotropic fibers [21] Compared to the amount of work that has been published on the composites with a homogeneous interphase, there are far fewer publications dealing with thermal or thermomechanical analysis of the composites with an inhomogeneous interphase Taking the variation of Youngs modulus of the interphase to be constant, or sub ject to power variation, reciprocal variation and cubic variation, respectively, Ja yaraman and Reifsnider investigated thermal residual stresses in fiber-reinforced composites with an inhomogeneous interphase [22]. Using the same variations to describe the thermal expansion coefficient of the interphase, they also discussed the effect of the interphase on residual thermal stresses in the composites [23]. In- troducing adhesion parameters to describe the imperfect adhesions between main constituents and taking into account the inhomogeneous interphase, Kakavas et al. determined th e coeficient of thermal expansion of fiber-reinforced compos- The introduction of the continuously varying material properties of the inhomoge- neous interphase greatly increases the complexity to achieve an analytical solution of the thermal or thermomechanical analysis of the composites. In this paper, an effort will be made to address this issue. By employing a linear variation of mater ial properties of the inhomogeneous interphase, a second order ordinary differential equation will be derived from the assumption of the generalized plane strain and the basic equations of axisymmetric problems in elasticity. The equation is further changed into a hypergeometric equation by variable substitution whose accurate analytical solution is represented with the hypergeometric function. The bound ary conditions of the composites and the continuous conditions of deformations and stresses between the constituents are employed to determine all the unknown constants. The obtained analytical solution is used to examine the effects of the inhomogeneous interphase, material properties of the matrix on thermal stresses in the composites
Influence of inhomogeneous interphase on composites 51 on its fibers [17]. You gave an efficient method to tackle thermal analysis of composites containing anisotropic continuous fibers and an arbitrary number of interfacial layers [18]. The above analytical research studies did not consider the plasticity of composites. When thermal deformations of composites are large enough, some components of composites may become plastic. In order to address this problem, Zhang et al. used a two concentric cylinder model together with the increment theory of plasticity to derive two partial differential equations governing the deformation behavior of composites and numerically solved these two partial differential equations [19]. Adopting Tresca’s yield criterion and linear strain-hardening material model, You and Long discussed the effects of material properties of interfacial layer on thermal response of fiber-reinforced composites [20]. This work was further extended to the composites reinforced with anisotropic fibers [21]. Compared to the amount of work that has been published on the composites with a homogeneous interphase, there are far fewer publications dealing with thermal or thermomechanical analysis of the composites with an inhomogeneous interphase. Taking the variation of Young’s modulus of the interphase to be constant, or subject to power variation, reciprocal variation and cubic variation, respectively, Jayaraman and Reifsnider investigated thermal residual stresses in fiber-reinforced composites with an inhomogeneous interphase [22]. Using the same variations to describe the thermal expansion coefficient of the interphase, they also discussed the effect of the interphase on residual thermal stresses in the composites [23]. Introducing adhesion parameters to describe the imperfect adhesions between main constituents and taking into account the inhomogeneous interphase, Kakavas et al. determined the coefficients of thermal expansion of fiber-reinforced composites [24]. The introduction of the continuously varying material properties of the inhomogeneous interphase greatly increases the complexity to achieve an analytical solution of the thermal or thermomechanical analysis of the composites. In this paper, an effort will be made to address this issue. By employing a linear variation of material properties of the inhomogeneous interphase, a second order ordinary differential equation will be derived from the assumption of the generalized plane strain and the basic equations of axisymmetric problems in elasticity. The equation is further changed into a hypergeometric equation by variable substitution whose accurate analytical solution is represented with the hypergeometric function. The boundary conditions of the composites and the continuous conditions of deformations and stresses between the constituents are employed to determine all the unknown constants. The obtained analytical solution is used to examine the effects of the inhomogeneous interphase, material properties of the matrix on thermal stresses in the composites
LH. You 2. ANALYTICAL FORMULAE OF DEFORMATIONS AND STRESSES Fiber-reinforced composites can be simplified as a three-phase model consisting of the inner fiber, middle interphase and outer matrix whose outer radii are rf ri and rm, respectively. In order to obtain an accurate analytical solution, all the constituents of the composites are taken to be isotropic. The variations of Youngs modulus and thermal expansion coefficient of the interphase are taken to be linear Since the Poissons ratio of the interphase has a weak influence on thermal stresses in the composites, it is assumed to be constant, i. e. does not vary along the radial direction. Using the subscript i to indicate the interphase, the variations of Youngs modulus and thermal expansion coefficient of the interphase can be written as Ei= eo+ eir Mi=ao +ar According to the following continuous conditions of material properties at the interfaces between the fiber and interphase and between the interphase and matrix r, E= et f Ei= em, ai=a the unknown constants in equation(1) are found to be r E ri E Subjected to a uniform temperature change, the deformations of the composites can be decomposed into elastic and thermal ones, respectively. Using Hooke's law in elasticity to describe the elastic deformations, the total radial, circumferential and axial strains can be formulated below E[t,-w(+0)+a△T [o-v(o+a)+a△T a2-v1(+a)]+a1△T According to the assumption of the generalized plane strain, the axial strain in all the constituents of the composites is equal to a uniform strain Ez0. Substituting it and
52 L. H. You 2. ANALYTICAL FORMULAE OF DEFORMATIONS AND STRESSES Fiber-reinforced composites can be simplified as a three-phase model consisting of the inner fiber, middle interphase and outer matrix whose outer radii are rf , ri and rm, respectively. In order to obtain an accurate analytical solution, all the constituents of the composites are taken to be isotropic. The variations of Young’s modulus and thermal expansion coefficient of the interphase are taken to be linear. Since the Poisson’s ratio of the interphase has a weak influence on thermal stresses in the composites, it is assumed to be constant, i.e. does not vary along the radial direction. Using the subscript i to indicate the interphase, the variations of Young’s modulus and thermal expansion coefficient of the interphase can be written as Ei = E0 + E1r, αi = α0 + α1r. (1) According to the following continuous conditions of material properties at the interfaces between the fiber and interphase and between the interphase and matrix r = rf , Ei = Ef , αi = αf , r = ri, Ei = Em, αi = αm, (2) the unknown constants in equation (1) are found to be E0 = riEf − rf Em ri − rf , E1 = Em − Ef ri − rf , (3) α0 = riαf − rf αm ri − rf , α1 = αm − αf ri − rf . Subjected to a uniform temperature change, the deformations of the composites can be decomposed into elastic and thermal ones, respectively. Using Hooke’s law in elasticity to describe the elastic deformations, the total radial, circumferential and axial strains can be formulated below εr = 1 Ei σr − νi(σθ + σz) + αiT, εθ = 1 Ei σθ − νi(σr + σz) + αiT, (4) εz = 1 Ei σz − νi(σθ + σr) + αiT. According to the assumption of the generalized plane strain, the axial strain in all the constituents of the composites is equal to a uniform strain εz0. Substituting it and
Influence of inhomogeneous interphase on composites the following equilibrium equation +σr-0=0 into the third of the equations(4), the axial stress can be written as 0=2,她+E1(-a△7) The substitution of equations (5)and(6) into the first two of equation(4)leads to the total radial and circumferential strains as follows ,=(--2)0-(1+的 vE20+(1+)a;△T, (1-u-202)0,+(1-v2)r d-/-le:0+(1+)x;△T Substituting the radial and circumferential strains in the above equation into the following deformation compatibility equation and making use of equation (1), a second order ordinary differential equation describing the deformation of the composites is achieved below (Eoter)r+(3Eo+2Er) E1 (E0+E1r)2=0.(9) The general solution of the above equation consists of the solution of its homoge neous equation and the particular solution of equation(9). Dividing equation(9)by Eo and denoting no= E1/Eo, the homogeneous form of equation(9)can be written (1+m0d2+(3+2nor) noo =0 (10) If we introduce a new variable z= -nor, equation(10) can be changed into the following hypergeometric equation -0, The solution of the above hypergeometric equation can be represented with the hypergeometric function which has the form of [25]
Influence of inhomogeneous interphase on composites 53 the following equilibrium equation r dσr dr + σr − σθ = 0, (5) into the third of the equations (4), the axial stress can be written as σz = 2νiσr + νir dσr dr + Ei(εz0 − αiT ). (6) The substitution of equations (5) and (6) into the first two of equation (4) leads to the total radial and circumferential strains as follows εr = 1 Ei (1 − νi − 2ν2 i )σr − νi(1 + νi)r dσr dr − νiεz0 + (1 + νi)αiT, (7) εθ = 1 Ei (1 − νi − 2ν2 i )σr + (1 − ν2 i )r dσr dr − νiεz0 + (1 + νi)αiT. Substituting the radial and circumferential strains in the above equation into the following deformation compatibility equation r dεθ dr + εθ − εr = 0, (8) and making use of equation (1), a second order ordinary differential equation describing the deformation of the composites is achieved below (E0+E1r)r d2σr dr2 +(3E0+2E1r) dσr dr −1 − 2νi 1 − νi E1σr+ α1T 1 − νi (E0+E1r)2 = 0. (9) The general solution of the above equation consists of the solution of its homogeneous equation and the particular solution of equation (9). Dividing equation (9) by E0 and denoting n0 = E1/E0, the homogeneous form of equation (9) can be written as (1 + n0r)r d2σr dr2 + (3 + 2n0r) dσr dr − 1 − 2νi 1 − νi n0σr = 0. (10) If we introduce a new variable z = −n0r, equation (10) can be changed into the following hypergeometric equation (1 − z)z d2σr dz2 + (3 − 2z) dσr dz + 1 − 2νi 1 − νi σr = 0. (11) The solution of the above hypergeometric equation can be represented with the hypergeometric function which has the form of [25]