# Radially Symmetric Steady State Thermal and Mechanical Stresses of a Poro FGM Hollow Sphere

Authors: M. Jabbari; S. Karampour; M. R. Eslami
Journal: ISRN Mechanical Engineering
 (2011)

---

## Abstract

A general solution for the one-dimensional steady-state thermal and mechanical stresses in a hollow thick sphere made of porous functionally graded material  (FGPM) is presented. The temperature distribution is assumed to be a function of radius, with general thermal and mechanical boundary conditions on the inside and outside surfaces of the sphere. The material properties, except Poisson's ratio, are assumed to vary along the radiusr according to a power law function. The analytical solution of the heat conduction equation and the Navier equation lead to the temperature profile, radial displacement, radial stress, and hoop stress as a function of radial direction.


---

## Body

## 1. Introduction

Functionally graded materials (FGMs) are made of a mixture with arbitrary composition of two different materials, and the volume fraction of each material changes continuously and gradually. The FGMs concept is applicable to many industrial fields such as chemical plants, electronics, and biomaterials [1]. Thick hollow sphere analysis made of FGM under mechanical and thermal loads and in asymmetric and two-dimensional (r,θ) state was conducted investigating navier equations and using legendre polynomials [2]. The analytical solution for the stresses in spheres and cylinders made of functionally graded materials are given by Lutz and Zimmerman [3, 4]. They considered thick spheres and cylinders under radial thermal loads, where radially graded materials with linear composition of the constituent materials were considered. Obata and Noda [5] studied one-dimensional steady thermal stresses in a functionally graded circular hollow cylinder and a hollow sphere using the perturbation method. The transient thermal stresses in a plate made of FGM are presented by the same authors [6]. By introducing the theory of laminated composites, Ootao and Tanigawa derived the three-dimensional transient thermal stresses of a nonhomogeneous hollow sphere with a rotating heat source [7]. Jabbari et al. presented the analytical solution of one and two-dimensional steady state thermoelastic problems of the FGM cylinder [8]. Two-dimensional nonaxisymmetric transient mechanical and thermal stresses in a thick hollow cylinder is presented by Jabbari et al. [9]. Shao and Wang presented three-dimensional solution to obtain stress fields in an FGM cylindrical plate with a finite length which is mechanical and thermal load tests and showed the results an a graph [10]. The transient thermal stress problem in a hollow sphere with homogeneous and isotropic properties is analytically solved by Cheung et al. [11]. Their assumed thermal boundary conditions are varied along the γ direction and the problem is solved with the potential function method. Takeuti and Tanigawa have employed the potential function method to obtain the analytical solution of a homogeneous spherical vessel with rotating heat source [12]. They considered the general form of thermal boundary conditions and analyzed the problem in transient condition with the potential functions method.Porous spheres of nanometer to micrometer dimensions are being pursued with great interest because of several possible technical applications in catalysis, drug delivery systems, separation techniques, photonics, as well as piezoelectric and other dielectric devices [13, 14]. The study of the thermomechanical response of fluid saturated porous materials is important for several branches of engineering [15–19]. Some of the important cases are the disposal of high-level radioactive waste, the extraction of oil or geothermal energy, the storage of hot fluids, and road subgrade or furnace foundations, which are usually subjected to cyclic changes of temperature. Wang and Papamichos [20, 21] presented analytical solutions for the temperature, pore pressure and stresses around a cylindrical wellbore and a spherical cavity subjected to a constant temperature change and a constant fluid flow rate by coupling the conductive heat transfer with the pore-fluid flow. Using the so-called heat source function method. Kurashige [22] analyzed the heat and fluid flow problem of a spherical cavity subjected to a constant heat flux in an impermeable medium.Inspite of conducted studies on spherical and cylindrical vessel made of FGM to obtain mechanical displacements and mechanical and thermal stresses, there is not any study on composition of poro and FGM materials. This study investigates the effect of mechanical and thermal stresses on one-dimensional steady state in the poro FGM hollow sphere that fluid trapped in the poro medium which is located in undrained conditions.


## 2. Analysis

Consider a thick hollow sphere of inside radiusa and outside radius b made of poro FGM. The sphere material is graded through the radial r-direction. Thus, the material properties are functions of r. Let u be the displacement component along the radial direction. The strain-displacement relations are
(1)εrr=u,r,εθθ=εϕϕ=ur.
The stress-strain relations one-dimensional in the poro FGM hollow sphere with fluid trapped in the poro pure are
(2)σrr=C11εrr+2C12εθθ-γpδrr-zrT(r),σθθ=σϕϕ=C12εrr+(C22+C23)εθθ-γpδθθ-zθT(r),
where
(3)zr=C11αr+2C12αθ,zθ=C12αr+(C22+C23)αθ,αr=αθ,
where σij and εij(i,j=r,θ) are the stress and strain tensors, T(r) is the temperature distribution determined from the heat conduction equation, αi is the coefficient of thermal expansion of effective stress, and Cij, γ, p, and δii(i=r,θ) are elastic constants, Biots coefficient of effective stress, poro pressure, and kronecker delta.Equations in a porous material are(4)p=M(ζ-γ∈),M=ku-kγ2,B=ku-kγku,ku=k[1+γ2kf(1-γ)(γ-ϕp)kf+ϕpkf].
For fluid in undrained condition,
(5)ζ=0,p=-γM∈=-γM(εrr+εθθ+εϕϕ),εθθ=εϕϕ,
where ζ, ∈, M, ku, k, kf, B, and ϕpare variation of fluid content, volumetric strain, Biots modulus, undrained bulk modulus, drained bulk modulus, bulk modulus of fluid, compressibility coefficient, and porosity.The equilibrium equation in the radial direction, disregarding the body force and the inertia term, is(6)σrr,r+2r(σrr-σθθ)=0.
Substituting (5) into (2) leads to
(7)σrr=c1*εrr+c2*εθθ-zrT(r),σθθ=σϕϕ=c3*εrr+c4*εθθ-zθT(r),
where
(8)c1*=C11+γ2m,c2*=2(C12+γ2m),c3*=C12+γ2m,c4*=C22+C23+2γ2m.
The sphere’s material is assumed to be described with a power law function of the radial direction as
(9)ci*=c¯i*rm,zi=z̅ir2m,
where m is the power law index of the material.Using relations (1)–(9), the Navier equation in term of the displacement is
(10)u,rr+1r((m+2)c¯1*+c¯2*-2c¯3*C¯1*)u,r+1r2((m+1)c¯2*-2c¯4*c¯1*)=z¯r(c¯1*(m+1+2L)rm-2+2c¯2*(m+1-L))r2m-1,
where
(11)L=z¯θz¯r.


## 3. Heat Conduction Problem

The heat conduction equation in the steady-state condition for the one-dimensional problem in spherical coordinates and the thermal boundary conditions for a poro FGM hollow sphere are given, respectively, as(12)1r2(r2k(r)T′(r))′=0,c11T′(a)+c12T(a)=f1,c21T′(b)+c22T(b)=f2,
where k=k(r) is the thermal conduction coefficient and cij are either thermal conduction coefficient k, or convection coefficient h, depending on the type of thermal boundary conditions. The terms f1 and f2 are known constants on the inside and outside radii. It is assumed that the thermal conduction coefficient k(r) is a power function of r as
(13)k=k(r)rm.
Using (13), the heat conduction equation becomes
(14)1r2[rm+2T′(r)]′=0.
Integrating (12) twice yields
(15)T(r)=c1r-(m+1)+c2.
Using the boundary conditions (16) and (17) to determine the constants c1 and c2 yields(16)c1=c22f1-c12f2c12((m+1)c21b-(m+2)-c22b-(m+1))-c22((m+1)c11a-(m+1)-c11a-(m+1)),(17)c2=f1((m+1)c21b-(m+2)c22b-(m+1))-f2((m+1)c11a-(m+2)-c12a-(m+1))c12((m+1)c21b-(m+2)c22b-(m+1))-c22((m+1)c11a-(m+2)-c12a-(m+1)).


## 4. Solution of the Navier Equation

The Navier equation for the radial displacementu was given in (10). Equation (10) is the nonhomogeneous Euler differential equation with the general and particular solutions. The general solution, ug, is obtained by assuming(18)ug(r)=Brη.
Substituting (18) into the homogeneous form of (10) yields
(19)η2+((m+1)c¯1*+c¯2*-2c¯3*c¯1*)η+((m+1)c¯2*-2c¯4*c¯1*)=0,
where (19) has two real roots η1 and η2 as(20)η1,2=-((m+1)c¯1*+c¯2*+2c¯3*)±((m+1)c¯1*+c¯2*+2c¯3*)2-4c¯1*((m+1)c¯2*-2c¯4*)2c¯1*.Thus, the general solution is
(21)ug(r)=B1rη1+B2rη2.
The particular solution up(r)is assumed to have the form
(22)up(r)=D1rm+D2r2m+1.Substituting (22) in (10) yields
(23)(m(m-1)+((m+2)c¯1*+c¯2*-2c¯3*c¯1*)m+((m+1)c¯2*+2c¯4*c¯1*))D1rm-2+((2m+1)(2m)+((m+2)c¯1*+c¯2*+2c¯3*c¯1*)×(2m+1)+((m+1)c¯2*+2c¯4*c¯1*))D2r2m-1=z¯r(c¯1*(m+1-2L)rm-2+2c¯2*(m+1-L)r2m-1).
Equating the coefficients of the identical powers yields(24)D1=z¯r(c¯1*(m+1)-2L)m(m-1)c¯1*+((m+2)c¯1*+c¯2*+2c¯3*)m+(m+1)c¯2*+2c¯4*,(25)D2=z¯r(2c¯2*(m+1-L))(2m+1)(2m)c¯1*+((m+2)c¯1*+c¯2*+2c¯3*)(2m+1)+((m+1)c¯2*+2c¯4*). The complete solution for u(r) is the sum of the general and particular solutions as
(26)u(r)=ug(r)+up(r).
Thus,
(27)u(r)=B1rη1+B2rη2+D1rm+D2r2m+1.
Substituting (27) into (1) and (2), the strains and stresses are obtained as
(28)εrr=B1η1rη1-1+B2η2rη2-1+D1mrm-1+D2(2m+1)r2m,εθθ=B1rη1-1+B2rη2-1+D1rm-1+D2r2m,σrr=c¯1*[S(c¯1*r-(m+1)+c¯2*)B1(η1+h)rη1-1+B2(η2+h)rη2-1+D1(m+h)rm-1+D2(2m+1+h)r2m-S(c¯1*r-(m+1)+c¯2*)],σθθ=c¯1*[B1(η1+f)rη1-1+B2(η2+f)rη2-1+D1(m+f)rm-1+D2(2m+1+f)r2m-g(c¯1*r-(m+1)+c¯2*)],
where
(29)h=c¯2*c¯1*,S=z¯rc¯1*,f=c¯4*c¯3*,g=z¯θc¯3*.
To determine the constants B1 and B2, the boundary conditions for stresses may be used. Consider the mechanical boundary conditions on the inside and outside radii as
(30)σrr(a)=-Pa,σrr(b)=-Pb.


## 5. Numerical Results and Discussion

The analytical solution obtained in the previous section may be checked for a one of example. Then, consider a thick hollow sphere of inner radiusa = 1 m, and the outer radius b = 1.2 m of poro FGM material with properties are given in Table 1. The boundary conditions for temperature are taken as T(a)=10°C and T(b)=0°C. The hollow sphere may be assumed to be under internal pressure of 50 MPa and zero external pressure (σrr(a) = −50 MPa andσrr(b) = 0 MPa). For different values of m, temperature profile, radial displacement, radial stresses, and hoop stresses along the radial direction are plotted in Figures 1–4.  Figure 1 shows that as the power law increases, the temperature decreases.  Figure 2 shows that for higher values of m, radial displacement decreases. Figure 3 represents the radial stress along the radial direction and decreases as the power law index increases. The circumferential stress versus the radial direction is shown in Figure 4. It is seen that for  m<1, the circumferential stress decreases along the radial direction. When m>1, the situation is reversed, and the circumferential stress increases along the radial direction. The curve associated with m=1 shows that the variation of circumferential stress along the radial direction is minor and is almost uniform across the radius. Effective porosity and compressibility on mechanical radial displacement in Figures 5 and 6, radial stress in Figures 7 and 8 and hoop stress in Figures 9 and 10 for m=1 are shown. As noted from Figures 5–10, distribution that with increasing porosity increases radial displacement, and all mechanical stresses decrease. Also, with increasing compressibility, radial displacement decreases, and all mechanical stresses increases.Table 1
Matrial constants.
Elastic constantsPoro constantsC11=1.867×1011γ=0.27C12=0.819×1011ku=4.1×1010C22=1.628×1011k=3.5×1010C23=0.270×1011kf=3.3×103αr=1.2×10-6B=0.45ϕp=0.01Figure 1
Radial temperature distribution.Figure 2
Radial displacement distribution.Figure 3
Radial stress distribution.Figure 4
Circumferential stress distribution.Figure 5
Effective porosity on mechanical radial displacement form=1.Figure 6
Effective compressibility on mechanical radial displacement form=1.Figure 7
Effective porosity on radial stress form=1.Figure 8
Effective compressibility on radial stress form=1.Figure 9
Effective porosity on circumferential stress form=1.Figure 10
Effective compressibility on circumferential stress form=1.



---

*Source: 305402-2011-07-28.xml*