Jtam-A4.dvi JOURNAL OF THEORETICAL AND APPLIED MECHANICS 55, 4, pp. 1155-1165, Warsaw 2017 DOI: 10.15632/jtam-pl.55.4.1155 THERMAL CREEP STRESS AND STRAIN ANALYSIS IN A NON-HOMOGENEOUS SPHERICAL SHELL Pankaj Thakur Department of Mathematics, Faculty of Science and Technology, ICFAI University Baddi, Solan, India e-mail: pankaj thakur15@yahoo.co.in; dr pankajthakur@yahoo.com Satya Bir Singh Department of Mathematics, Punjabi University Patiala, Punjab, India; e-mail: sbsingh69@yahoo.com D.S. Pathania Department of Mathematics, Guru Nanak Dev Engineering College, Ludhiana, Punjab, India e-mail: despathania@yahoo.com Gaurav Verma Research Scholar, IKG Punjab Technical University Kapurthala, Punjab, India; e-amil: gkdon85@gmail.com The purpose of this paper is to present study of thermal creep stress and strain rates in a non-homogeneous spherical shell byusingSeth’s transition theory. Seth’s transition theory is applied to the problem of creep stresses and strain rates in the non-homogeneous spherical shell under steady-state temperature. Neither the yield criterion nor the associated flow rule is assumed here. With the introduction of thermal effect, values of circumferential stress decrease at the external surface as well as internal surface of the spherical shell. It means that the temperature dependentmaterialsminimize the possibility of fracture at the internal surface of the spherical shell. The model proposed in this paper is used commonly as a design of chemical and oil plants, industrial gases and stream turbines, high speed structures involving aerodynamic heating. Keywords: stress, strain rates, thermal, spherical shell, non-homogeneous 1. Introduction Spherical shell structures have found widespread use in modern technology such as design of chemical and oil plants, accumulator shells, pressure vessel for industrial gases or media trans- portation of high-pressurized fluids and piping of nuclear containment, high speed structures involving aerodynamic heating, submerged undersea structures, earth sheltered domes, and the like. These spherical systems are effective from the perspectives of both structural and architec- tural design. Inmany of these cases, the spherical shells have to operate under severemechanical and thermal loads causing significant creep and thus reducing its service life. The collapse or da- mage is initiatedbycreep, shrinkageandthermal effects, or fromtheir interactionwith structures that both experience or do not experience environmental degradation. Consequently, demand for strengthening and upgrading the existing concrete structures, because of damage caused by long-termeffects andexcessive structural deformations, hasbeen recognized.However, before the application of costly strengthening techniques, understanding of nonlinear long-term behaviour of the existing and new spherical shells is essential, and the development of suitable and reliable theoretical approaches for their analysis and safety assessment is required. Creep effects gene- rally increase deformations of a shell structure even under room temperatures, and are usually only considered to affect behaviour at the service ability limit states. Therefore, the analysis of 1156 P. Thakur et al. long term steady state creep deformations of shells is very important in these applications (Ha- med et al., 2010; Kashkoli andNejad, 2014). Due to the occurrence of these creep deformations, non-homogeneous materials are widely used in the engineering applications. Non-homogeneous materials are a specific class of composite materials known as functionally graded materials (FGM) in which constituents are graded in one or more direction with continuous variation to achieve desired properties. The smooth grading of the constituents result in better thermal properties, higher fracture toughness, improved residual stress distribution and reduced stress intensity factors. These properties allow non-homogeneous structures towithstandhigh pressure under elevated thermal environment.Therefore, the analysis of non-homogeneity in the spherical shell through a mathematical model by taking one and all the complexities into consideration is the major concern of researchers (Kar and Panda, 2016). Some degree of non-homogeneity is present in a wide class of materials such as hot rolled metals, magnesium and aluminum alloys. Non-homogenity can also be introduced by a certain external field which is a thermal gradient material as the elastic moduli of the materials vary with temperature (Olszak, 1960). Penny (1967) obtained the effects of creep in spherical shells by analysis similar to the corresponding elastic ones described here. Miller (1995) evaluated solutions for stresses and displacements in a thick spherical shell subjected to internal and external pressure loads. You et al. (2005) pre- sented an accurate model to carry out elastic analysis of thick-walled spherical pressure vessels subjected to internal pressure. Kellogg andKing (1997) developed a finite element model of co- nvection in a spherical axisymmetric shell that we use to simulate upwelling thermal plumes in the mantle. Thakur (2011) analyzed creep transition stresses of a thick isotropic spherical shell by finitesimal deformation under steady state of temperature and internal pressure by using Seth’s transition theory. Seth’s transition theory does not acquire any assumptions like the yield condition, incompressibility condition, and thus poses and solves a more general problem from which cases pertaining to the above assumptions can be worked out. This theory utilizes the concept of a generalized strain measure and asymptotic solution at critical points or turning points of differential equations defining the deformed field and has been successfully applied to a large number of problems (Seth, 1962, 1966; Thakur, 2011, 2014; Thakur et al., 2016, 2017). Seth (1962) defined the concept of generalized strain measures as eii = A eii ∫ 0 ( 1−2 A eii ) n 2 −1 d A eii= 1 n [ 1− ( 1−2 A eii ) n 2 ] i =1,2,3 (1.1) wheren is themeasure and A eii are theAlmansi finite strain components. Forn =−2,−1,0,1,2 it gives theCauchy,GreenHencky, Swainger andAlmansimeasures respectively.Non-homogeneity in a spherical shell has been taken as the compressibility of the material as C = C0r −k (1.2) where a ¬ r ¬ b, C0 and k are real constants. 2. Governing equations We consider a spherical shell whose internal and external radii are a and b, respectively, and is subjected to uniform internal pressure pi of gradually increasingmagnitude and temperature Θ0 applied to the internal surface r = a as shown in Fig. 1. The components of displacement in spherical co-ordinates are given by (Seth, 1962, 1966) u = r(1−β) v =0 w =0 (2.1) Thermal creep stress and strain analysis in a non-homogeneous spherical shell 1157 Fig. 1. Geometry of the spherical shell where u, v, w are displacement components, β is a position function depending on r. The generalized components of strain are given by Seth (1966) err = 1 n [1− (rβ′+β)n] eθθ = 1 n (1−βn)= eϕϕ erθ = eθϕ = eϕr =0 (2.2) where n is the measure and β′ = dβ/dr. Stress-strain relation. The stress-strain relations for a thermo-elastic isotropicmaterial are given by (Parkus, 1976) Tij = λδijI1+2µeij − ξΘδijTij i,j =1,2,3 (2.3) where Tij are stress components and eij is strain component, λ and µ are Lame’s constants, I1 = ekk is the first strain invariant, δij is Kronecker’s delta, ξ = α(3λ +2µ), α being the coefficient of thermal expansion, and Θ is temperature. Further, Θ has to satisfy ∇2Θ =0 (2.4) Substituting the strain components from Eq. (2.2) in Eq. (2.3), the stresses are obtained as Trr = 2µ n [1− (rβ′+β)n]+ λ n [3− (rβ′+β)n−2βn]−ξΘ Tθθ = Tϕϕ = 2µ n (1−βn)+ λ n [3− (rβ′+β)n−2βn]−ξΘ Trθ = Tθϕ = Tϕr =0 (2.5) Equation of equilibrium. The radial equilibrium of an element of the spherical shell requires dTrr dr + 2 r (Trr−Tθθ)= 0 (2.6) where Trr and Tθθ are the radial and hoop stresses, respectively. Boundary conditions. The temperature satisfying Laplace equation (2.4) with boundary condi- tions Θ = Θ0 ∧ Trr =−pi at r = a Θ =0 ∧ Trr =0 at r = b (2.7) where Θ0 is constant, is given by (Parkus, 1976) Θ = Θ0 log(r/b) log(a/b) (2.8) 1158 P. Thakur et al. Critical points or turning points. Using Eq. (2.5) in Eq. (2.6), we get a non-linear differential equation in β as nP(P +1)n−1βn+1 dP dβ = (µ′ µ − C′ C ) {3−2C −βn[2(1−C)+(1+P)n]}−2C′r(1−βn) −nβnP[2(1−C)+(1+P)n]+2Cβn[1− (1+P)n]− nCΘ0 2µβn ( ξ +rξ′ log r b ) (2.9) where Θ0 = Θ0/ log(a/b), C =2µ/λ+2µ and rβ ′ = βP (P is a function of β and β is a function of r only). The transition or turning points of β in Eq. (2.9)) are P →−1 and P →±∞. 3. Analytical solution For finding thermal creep stresses and strain rates, the transition function is taken through the principal stress difference (see Seth, 1962, 1966; Thakur, 2011, 2014; Thakur et al., 2016, 2017) at the transition point P →−1.We define the transition function ψ as ψ = Trr−Tθθ = 2µβn n [1− (P +1)n] (3.1) where ψ is a function of r only, and ψ is the dimension. Taking logarithmic differentiation of Eq. (3.1) with respect to r and substituting the value of dP/dβ from Eq. (2.9), we get d logψ dr = np r + µ′ µ − 1 rβn[1− (1+P)n] [ r (µ′ µ − C′ C ) {3−2C −βn[2(1−C)+(1+P)n]} −2rC′(1−βn)−nβnP[2(1−C)+(1+P)n] +2Cβn[1− (1+P)n]− nCΘ 2µ ( ξ+rξ′ log r b ) ] (3.2) Taking asymptotic value of Eq. (3.2) at P →−1, we get d dr (logψ)= 3µ′ µ − 2C′ C − 3n r +X (3.3) where X = 2(n−1)C r − 2Cµ′ µ + 2C′ βn − (µ′ µ − C′ C )3−2C βn + nCΘ0 2µrβn ( ξ+rξ′ log r b ) Integrating equation (3.3), we get ψ = A µ3 C2r3n exp(h) (3.4) whereh = ∫ X dr andA is a constant of integration,which canbedetermined fromtheboundary condition. From Eq. (3.1) and Eq. (3.4), we have Trr−Tθθ = A 2rµ3 2C2r3n+1 exp(h)= ArH 2 (3.5) where H = 2µ3 C2r3n+1 exp(h) Thermal creep stress and strain analysis in a non-homogeneous spherical shell 1159 Substituting Eq. (3.5) in Eq. (2.6) and integrating, we get Trr = B −A ∫ H dr (3.6) where B is a constant of integration, which can be determined from the boundary condition and asymptotic value of β as P →−1 is D/r, with D being a constant. Using boundary condition Eq. (2.7) in Eq. (3.6), we get A =− pi b ∫ a H dr B =− pi b ∫ a H dr ∫ H dr (3.7) where pi is pressure at the inner surface of the spherical shell. Using the integration constants A and B in Eq. (3.6), we get Trr = pi b ∫ a H dr b ∫ r H dr (3.8) Substituting Eq. (3.8) into Eq. (3.5), we get Tθθ = Tϕϕ = pi b ∫ a H dr ( b ∫ r H dr+ rH 2 ) (3.9) We introduce non-homogeneity in the spherical shell due to variable compressibility as given in Eq. (1.2), then Eq. (3.8) and Eq. (3.9) become Trr = pi b ∫ a H dr b ∫ r H1 dr Tθθ = Tϕϕ = pi b ∫ a H1 dr ( b ∫ r H1 dr+ rH1 2 ) (3.10) where Trr, Tθθ are radial and circumferential stresses, and H1 = 2µ3 C2r3n+1 exp(h1)= r−(3n+k+1) 4(1−C0r −k)3 exp(h1) h1 =− 2(n−1) k C0r −k− 2kC0r n−k Dn(n−k) + kC0 Dn ∫ rn−k−1(3−2C0)r −k 1−C0(bR) −k dr+2log(1−C0r −k) + αnΘ0 Dn ∫ (1−C0r −k) ( 3+ C0r −k 1−C0r −k − kC0r −k log(r/b) 1−C0r −k ) rn−1 dr C′ =−kC0r −k−1 µ ′ µ = C′ C(1−C) Equations (3.10)1,2 give thermal creep stresses for a spherical shell made of a non-homogeneous material under steady-state temperature. We introduce the following non-dimensional compo- nents as: R = r/b, R0 = a/b, σrr = τrr/pi, σθθ = τθθ/pi and αΘ0 = Θ1. Equations (3.10)1,2 in non-dimensional form become σrr = 1 ∫ R H2 dR 1 ∫ R0 H2 dr σθθ = σϕϕ = 1 1 ∫ R0 H2 dR ( 1 ∫ R H2 dR+ RH2 2 ) (3.11) 1160 P. Thakur et al. where H2 = (bR)−(3n+k+1)C0λ 3 4(1−C0b−kR−k)3 exp(h2) h2 =− 2(n−1) k C0(bR) −k− 2kC0(bR) n−k Dn(n−k) + kC0b n−k Dn ∫ Rn−k−1(3−2C0)R −k 1−C0(bR) −k dR+2log(1−C0b −kR−k) + nΘ1 Dn lnR0 ∫ (1−C0b −kR−k) ( 3+ C0(bR) −k 1−C0b −kR−k − kC0(bR) −k logR 1−C0b −kR−k ) Rn−1 dR Particular case: In the absence of temperature gradient (i.e. Θ1 =0), Eqs. (3.11), become σrr = 1 ∫ R H∗2 dR 1 ∫ R0 H∗2 dR σθθ = σϕϕ = 1 1 ∫ R0 H∗2 dR ( 1 ∫ R H∗2 dR+ RH∗2 2 ) (3.12) where H∗2 = (bR)−(3n+k+1)C0λ 3 4(1−C0b−kR−k)3 exp(h∗2) h∗2 =− 2(n−1) k C0(bR) −k− 2kC0(bR) n−k Dn(n−k) + kC0b n−k Dn ∫ Rn−k−1(3−2C0)R −k 1−C0(bR)−k dR+2log(1−C0b −kR−k) 4. Estimation of creep parameters When the creep sets in, the strains shouldbe replaced by strain rates, then stress-strain relations are given (Sokolnikoff, 1946; Parkus, 1976) eij = 1+ν E Tij − ν E δijT +αΘ (4.1) where eij is the strain component and T = Tii is thefirst stress invariant and ν =(1−C)/(2−C) is Poisson’s ratio. Differentiating Eq. (2.2) with respect to time t, we get ėθθ =−β n−1β̇ (4.2) For Swainger measure (i.e. n =1), Eq. (4.2) become ε̇θθ = β̇ (4.3) where ε̇θθ is the Swainger strain measure. FromEq. (3.1), the transition value β is given at the transition point P →−1 by β = ( n 2µ ) 1 n (Trr−Tθθ) 1 n (4.4) Using Eqs. (4.2)-(4.4) in Eq. (4.1), we get ε̇rr = m(σrr−νσθθ+αΘ) ε̇θθ = m(σθθ−νσrr+αΘ) ε̇ϕϕ =−m[ν(σrr+σθθ)+αΘ] (4.5) where ε̇rr, ε̇θθ and ε̇ϕϕ are strain rates, σrr, σθθ are stress components and m = [n(σr − σθ)(1 + ν)] 1 n −1. These are the constitutive equations used by Odquist (1974) for finding the creep stresses, provided we put n =1/N and N is the measure. Thermal creep stress and strain analysis in a non-homogeneous spherical shell 1161 5. Numerical results and discussion For calculating creep stresses and strain rates on the basis of the above analysis, the following values have been taken, ν = 0.5 (incompressible material C0 = 0), ν = 0.42857 (compressible material C0 =0.25) and 0.333 (compressiblematerial C0 =0.50), n =1/3, 1/5, 1/7 (i.e. N =3, 5, 7), thermal expansion coefficient α = 5.0 · 10−5degF−1 for Methyl Methacrylate (Levitsky andShaffer, 1975) andΘ1 = αΘ0 =0and 0.5, D =1. In classical theory, themeasure N is equal to 1/n. The definite integrals in Eqs. (3.11) have been evaluated by using Simpson’s rule. From Figs. 2-4, curves are presented between stresses along the radii ratio R = r/b in the spherical shell made of compressible as well as incompressible materials for k =−1, 0, 1. It can be seen form Figs. 2 and 3 that the circumferential stresses are maximum at the external surface for n =1/7 and k =−1, 0 for a compressible material as compared to the incompressible material. From Fig. 4, the circumferential stress is maximum at the internal surface for non-homogeneity k =1. The non-homogeneity increases the values of circumferential stress (i.e. k =−1, 0), but reverse the result for k =1. Fig. 2. Creep stresses in a non-homogeneous (k =−1) spherical shell along the radii ratio R = r/b; (a) without temperature, (b) for temperature Θ1 =0.5 1162 P. Thakur et al. Fig. 3. Creep stresses in a non-homogeneous (k =0) spherical shell along the radii ratio R = r/b; (a) without temperature, (b) for temperature Θ1 =0.5 With the introduction of a temperature gradient, the values of circumferential stress are de- creased at the external surface aswell as internal surface of the spherical shell for different values of non-homogeneity. It means that temperature dependentmaterials minimize the possibility of fracture at the internal surface of the spherical shell. From (Fig. 5), curves are producedbetween strain rates along the radii ratio R = r/b for the spherical shell made of compressible material C =0.25, i.e. saturated clay for k =−1, 0, 1. It can be seen that the strain rates are maximum at the external surface for k = −1, 0 and reverse in the case k = 1. With the introduction of a temperature gradient, the strain rates decrease at the internal surface as well as the external surface. It means that the temperature dependentmaterials minimize the possibility of fracture at the internal surface of the spherical shell. Acknowledgment Theauthors aregrateful to the referee forhis critical commentswhich led to a significant improvement of the paper. References 1. Hamed E., Bradford M.A., Gilbert R.I., 2010, Nonlinear long-term behaviour of spherical shallow thin-walled concrete shells of revolution, International Journal of Solids and Structures, 47, 2, 204-215 Thermal creep stress and strain analysis in a non-homogeneous spherical shell 1163 Fig. 4. Creep stresses in a non-homogeneous (k =1) spherical shell along the radii ratio R = r/b; (a) without temperature, (b) for temperature Θ1 =0.5 2. Kar V.R., Panda S.K., 2016, Nonlinear thermo-mechanical deformation behaviour of P-FGM shallow spherical shell panel,Chinese Journal of Aeronautics, 29, 1, 173-183 3. KashkoliM.D.,NejadM.Z., 2014,Effect of heatfluxoncreep stresses of thickwalledcylindrical pressure vessels, Journal of Applied Research and Technology, 12, 3, 585-597 4. Kellogg L.H., King S.D., 1997, The effect of temperature dependent viscosity on the structure of new plumes in the mantle: Results of a finite element model in a spherical, axisymmetric shell, Earth and Planetary Science Letters, 148, 1/2, 13-26 5. Levitsky M., Shaffer, B.W., 1975, Residual thermal stresses in a solid sphere from a thermo- setting material, Journal of Applied Mechanics, Transactions of ASME, 42, 3, 651-655 6. Miller G.K., 1995, Stresses in a spherical pressure vessel undergoing creep and dimensional changes, International Journal of Solids and Structures, 32, 14, 2077-2093 7. Odquist F.K.G., 1974, Mathematical Theory of Creep and Creep Rupture, Clarendon Press, Oxford 8. Olszak W., 1960, Non-homogeneity in elasticity and plasticity, ZAMM, 40, 10/11, 522-523 9. Parkus H., 1976,Thermo-Elasticity, Springer-VerlagWien, NewYork, USA 1164 P. Thakur et al. Fig. 5. Strain rates in a non-homogeneous spherical shell along the radii ratio R = r/b for C =0.25; (a) without temperature Θ1 =0, (b) for temperature Θ1 =0.5 10. Penny R.K., 1967, The creep of spherical shells containing discontinuities, International Journal of Mechanical Sciences, 9, 6, 373-388 11. Seth B.R., 1962, Transition theory of elastic-plastic deformation, creep and relaxation, Nature, 195, 896-897, DOI: 10.1038/195896a0 12. SethB.R., 1966,Measure concept inmechanics, International Journal of Non-LinearMechanics, 1, 1, 35-40 13. Sokolnikoff I.S., 1946,Mathematical Theory of Elasticity, 1st edition,Mc-GrawHill BookCom- pany, Inc., NewYork, 60-76 14. ThakurP., 2011, Creep transition stresses of a thick isotropic spherical shell by finitesimal defor- mation under steady state of temperature and internal pressure,Thermal Science,15, 2, S157-S165 15. Thakur P., 2014, Steady thermal stress and strain rates in a circular cylinder with non- homogeneous compressibility subjected to thermal load,Thermal Science, 18, 1, S81-S92 16. Thakur P., Singh S.B., Kaur J., 2016, Thermal creep stresses and strain rates in a circular disc with shaft having variable density,Engineering Computation, 33, 3, 698-712 Thermal creep stress and strain analysis in a non-homogeneous spherical shell 1165 17. ThakurP., GuptaN., Singh S.B., 2017,Creep strain rates analysis in cylinder under tempera- ture gradient for different material, Engineering Computations, 34, 3, DOI: 10.1108/EC-05-2016- -0159 18. You L.H., Zhang J.J., You X.Y., 2005, Elastic analysis of internally pressurized thick-walled sphericalpressurevessels of functionallygradedmaterials, International Journal of PressureVessels and Piping, 82, 5, 347-354 Manuscript received December 2, 2016; accepted for print April 15, 2017