ISRN Mechanical Engineering

VolumeΒ 2011Β (2011), Article IDΒ 291409, 9 pages

http://dx.doi.org/10.5402/2011/291409

## Mechanical and Thermal Stresses in a FGPM Hollow Cylinder due to Radially Symmetric Loads

^{1}Islamic Azad University, South Tehran Branch, Tehran, Iran^{2}Department of Mechanical Engineering, Academy of Sciences, Amirkabir University of Technology, Tehran, Iran

Received 25 May 2011; Accepted 21 June 2011

Academic Editor: J.Β Seok

Copyright Β© 2011 M. Jabbari et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

#### Abstract

The general solution of steady-state on one-dimensional Axisymmetric mechanical and thermal stresses for a hollow thick made of cylinder Functionally Graded porous material is developed. Temperature, as functions of the radial direction with general thermal and mechanical boundary-conditions on the inside and outside surfaces. A standard method is used to solve a nonhomogenous system of partial differential Navier equations with nonconstant coefficients, using complex Fourier series, rather power functions method and solve the heat conduction. The material properties, except poisson's ratio, are assumed to depend on the variable , and they are expressed as power functions of .

#### 1. Introduction

Poroelasticity is a theory that models the interaction of deformation and fluid flow in a fluid-saturated porous medium. The deformation of the medium influences the flow of the fluid and vice versa. The theory was proposed by Biot [1, 2] as a theoretical extension of soil consolidation models developed to calculate the settlement of structures placed on fluid-saturated porous soils. The historical development of the theory is sketched by De Boer [3]. The theory has been widely applied to geotechnical problems beyond soil consolidation, most notably problems in rock mechanics. Detournay and Cheng [4] survey both these methods with special attention to rock mechanics. These include familiar analytical methods (displacement potentials, method of singularities) and computational methods (finite element and boundary element). Sandhu and Wilson [5] are acknowledged for pioneering the application of finite element techniques to poroelasticity. Detournay and Cheng [6] presented fundamentals of poroelasticity.

Abousleiman and Ekbote [7] presented the analytical solutions for inclined hollow cylinder in a transversely isotropic material subjected to thermal and stress perturbations, and they systematically evaluated the effect of the anisotropy of the poromechanical material parameters as well as thermal material properties on stress and porous pressure distributions. Chen [8] presented and analyzed the problems of linear thermo elasticity in a transversely isotropic hollow cylinder of finite length by a direct power series approximation through the application of the Lanczos-Chebyshev method. Bai [9] presented then derived an analytical method solving the responses of a saturated porous media subjected to cyclic thermal loading by the Laplace transform and the Gauss-Lengender method of Laplace transform inversion. Wang and Papamichos [10, 11] presented analytical solution for the temperature, pore pressure, and stresses around a cylindrical well bore and a spherical cavity subjected to a constant fluid flow rate by coupling the conductive heat transfer with the pore-fluid flow. Ghassemi and Tao [12] presented influence of coupled chemo-poro-thermoelastic processes on pore pressure and stress distributions around a wellbore in swelling shale. Wirth and sobey [13] presented an axisymmetric and fully 3-D poroelastic model forth evolution of hydrocephalus. Yang and Zhang [14] presented poroelastic wave equation including the Biot/squirt mechanism and the solid/fluid coupling anisotropy. Arora and Tomar [15] presented the elastic waves along a cylindrical borehole in a poroelastic medium saturated by two immiscible fluids. Hamiel et al. [16] presented the coupled evolution of damage and porosity in poroelastic media theory and applications to the deformation of porous rocks. Ghassemi [17] presented stress and pore prepressure distribution around a pressurized, cooled crack in hollow permeability rock. Youssef [18] theory of generalized porothermoelasticity was presented. Jourine et al. [19] presented modeling poroelastic hollow cylinder experiments with realistic boundary conditions.

Functionally graded materials (FGMs) are heterogeneous materials in which the elastic and thermal properties change from one surface to the other, gradually and continuously. The material is constructed by smoothly changing materials. Since ceramic has good heat resistance and metal has high strength, ceramic-Metal FGM may work at super high-temperature or under high-temperature difference field. In effect, the governing equation for the temperature and stress distributions are coordinate dependent as the material properties are functions of position. Classical method of analysis is to combine the equilibrium equations with the stress-strain and strain-displacement relations to arrive at the governing equation in terms of the displacement components called the Navier equation. There are some analytical thermal and stress calculations for functionally graded material in the one-dimensional case for thick cylinders and spheres [20, 21]. The authors have considered the nonhomogeneous material properties as liner function of . Jabbari et al. [22] studied a general solution for mechanical and thermal stresses in a functionally graded hollow cylinder due to nonaxisymmetric steady-state load. They applied separation of variables and complex Fourier series to solve the heat conduction and Navier equation. Poultangari et al. [23] presented Functionally graded hollow spheres under non-axisymmetric thermomechanical loads. Shariyat et al. [24] presented nonlinear transient thermal stress and elastic wave propagation analyses of thick temperature-dependent FGM cylinders, using a second-order point-collocation method. LΓΌ et al. [25] presented elastic mechanical behavior of nanoscaled FGM films incorporating surface energies. Afsar and Sekine [26] presented inverse problems of material distributions for prescribed apparent fracture toughness in FGM coatings around a circular hole in infinite elastic media. Zhang and Zhou [27] presented a theoretical analysis of FGM thin plates based on physical neutral surface. Fazelzadeh and Hosseini [28] presented aerothermoelastic behavior of supersonic rotating thin-walled beams made of functionally graded materials. Ootao and Tanigawa [29] presented the transient thermoelastic problem of functionally graded thick strip due to nonuniform heat supply. They obtained the exact solution for the two-dimensional temperature change in a transient state, and thermal stress of a simple supported strip under the state of plane strain. Jabbari et al. [30] presented and studied the mechanical and thermal stresses in functionally graded hollow cylinder due to radially symmetric loads. They assumed the temperature distribution to be a function of radial direction. They applied a method to solve the heat conduction and Navier equations. Farid et al. [31] presented three-dimensional temperature dependent-free vibration analysis of functionally graded material curved panels resting on two-parameter elastic foundation using a hybrid semianalytic, differential quadrature method. Bagri and Eslami [32] presented Generalized coupled thermoelasticity of functionally graded annular disk considering the Lord-Shulman theory. Shariat and Eslami [33] presented buckling of thick functionally graded plates under mechanical and thermal loads. Jabbari et al. [34] studied an axisymmetric mechanical and thermal stresses in thick short length functionally graded material cylinder. They applied the separation of variables and complex Fourier series to solve the heat conduction and Navier equation. Thieme et al. [35] presented titanium powder sintering for preparation of a porous FGM destined as a skeletal replacement implant.

In this work, a direct method of solution of the Navier equations presented which does not have limitation of the potential function method as to handle the general type of mechanical and thermal under one-dimensional steady-state temperature distribution with general type of thermal and mechanical boundary conditions is considered. The functionally graded porous material properties of the cylinder are assumed to be expressed by power functions in . The Naviear equation terms of displacements are derived and solved analytically by the direct method, so any boundary conditions for stresses and displacements can be satisfied.

#### 2. Heat Conduction Problem

Consider a hollow circular cylinder of inner radius , outer radius made of functionally graded porous material (FGPM) respectively. Axisymmetric cylindrical coordinates are considered along the radial direction. The cylinderβs material graded through the direction, thus the material properties are porous and functions of . The first law of thermodynamics for energy equation in the steady-state condition for the FGPM on dimensional cylinder is: where is temperature distribution, is the thermal conduction coefficient, and symbol (/) denotes derivative with respect to .

The thermal boundary is assumed as where (,) denotes partial derivative, and are the constant thermal parameters related to conduction and convection coefficients. We assume that nonhomogeneous thermal conduction coefficient is power function of as , where and material parameter. Using the definition for the material properties, the temperature equation becomes Integrating (4) twice yields Using the boundary conditions (2) to determine the constants and yields where constants to are given in Appendix A.

#### 3. Stress Analysis

Let displacement components in the radial direction. Then strain-displacement relations are and stress-strain relations of a functionally graded porous cylinder for nonaxisymmetric condition are where , (), , , , , , and are stress tensors, strain tensors, Biotβs modulus, Biotβs coefficient of effective stress, thermal expansion coefficient, lameβs coefficient, and the pore pressure, respectively, related to the Biotβs modulus, volumetric strain and the variation of fluid content.

We assume that pore-cylinder if undrained condition then () as: where: Thus, The equilibrium equation in the radial direction, disregarding body force and the inertia terms, is To obtain the equilibrium equations in terms of the displacement components for the FGPM cylinder, the functional relationship and pore of the material properties must be known. Because the cylinder material is assumed to be graded along the -direction, the modulus of elasticity and coefficient of thermal expansion are material constant

assumed to be described with the power laws as where the coefficients are described as and is the inner radius.

Using the relations (7) to (14), the Navier equations in terms of the displacement components are The Navier equation (15) is nonhomogeneous system of partial differential equations with non-constant coefficients. We assume that .

#### 4. Solution of the Navier Equation

Equation (15) is the Euler differential equation with general and particular solutions.

The general solution is assumed to have the form Substituting (16) into (15) yields Equation (17) has two roots to . Thus, the general solutions are Thus, the general solution is The particular solutions are assumed as Substituting (20) into (18) yields The complete details for solution of (21) is presented in Appendix B.

The complete solutions for is sum of the general and particular solutions and are Thus Substituting (23) into (1) and (2), the strains and stresses are obtained as

To determine the constants and , consider the boundary conditions for stresses given by

#### 5. Numerical Results and Discussion

Consider a thick hollow cylinder of inner radius = 1β(m) and outer radius = 1.2β(cm), shown properties are given in Table 1. For simplicity of analysis, we consider that the power law of material properties is the same as . To examine the proposed solution method, two example problems are considered. The example problem may have some physical interpretation.

As the example, consider a thick hollow cylinder where the inside boundary is traction free with given temperature distribution of Table 2. The outside boundary is assumed to be radially fixed with zero temperature. Therefore, the assumed boundary conditions yield of Table 2.

Figure 1 shows the variations of the temperature along the radial direction for different values of the power law index. The figure shows that as the power law index increases, the temperature decreased.

Figure 2 shows the plot of the radial displacement along the radius. The magnitude of the radial displacement is decreased as the power index is increased.

The radial and circumferential stresses are plotted along the radial direction and shown in Figures 3 and 4, and the magnitude of the radial stress is increased as is increased. The hoop stress along the radius decreases for (similar to thick cylinders made of isotropic materials), due to the acting internal pressure and zero external pressure. For , the hoop stress increases as the radius increases, since the modulus of elasticity is an increasing function of the radius. Physically, this means that the outer layers of the cylinder are biased to maintain the stress due to their higher stiffness. There is a limiting value for , where the hoop stress remains almost a constant along the radius. For low values of the ratio (Figures 7 and 8). Figures 5 and 6 show the radial and hoop thermal stresses in the cross-section of the cylinder, respectively, where the pore compressibility coefficient () is changed, the other parameters are fixed. Figures 5 and 6 show these stresses based on the pore volume fraction; () is pore volume per total volume.

Figure 9 shows the radial displacements in the cross-section of the cylinder based on the pore compressibility coefficient () changing. Figure 10 also shows these displacements based on the pore volume fraction () changing.

#### 6. Conclusions

In the present work, an attempt has been made to study the problem of general solution for the thermal and mechanical stresses in a thick FGPM hollow cylinder due to the one-dimensional axisymmetric steady-state loads. The method of solution is based on the direct method and uses power series, rather than the potential function method. The advantage of this method is its mathematical power to handle both simple and complicated mathematical function for the thermal and mechanical stresses boundary conditions. The potential function method is capable of handling complicated mathematical functions as boundary condition. The proposed method does not have the mathematical limitations to handle the general types of boundary conditions which are usually countered in the potential function method.

#### Appendices

#### A. Compressibility Coefficients and Pore Volume Fraction

: compressibility coefficient, sometimes called the Skempton pore pressure coefficient.

: pore volume fraction is pore per unite total volume. β and are bulk modulus of the fluid phase and bulk modulus of the poroelastic medium under the drained condition, respectively.

#### B. Constants Material

where constants to are given

#### Nomenclature

: | Inner radius |

: | Thermal constant |

: | Outer radius |

: | Thermal constant |

: | Constant temperature parameters |

: | Mechanical and thermal constants |

: | Mechanical and thermal constants |

: | Constant mechanical parameters |

: | Inner and outer temperature boundary conditions |

: | Inner and outer mechanical boundary conditions |

: | Thermal conduction coefficient |

: | Material parameter |

: | Yongβs modulus |

: | Material constant |

: | Material parameter |

: | Cylinder coordinate |

: | Cylinder temperature |

: | Coefficient of sine Fourier series |

: | Displacement components |

: | Thermal expansion coefficient |

: | Material constant |

: | Lame coefficient |

: | Poissonβs ratio |

: | Undrained Poissonβs ratio |

: | The pore pressure |

: | Biotβs modulus |

: | Biotβs coefficient of effective stress |

: | Delta carancker |

: | The variation of fluid content (undrained ) |

: | Strain tensor |

: | Volumetric strain () |

: | Stress tensor |

: | Compressibility coefficient |

: | Pore volume fraction is pore per unite total volume. |

#### References

- M. A. Biot, βLe problème de la consolidation des matières argileuses sous une charge,β
*Annales de la Societe Scientifique de Bruxelles*, vol. B55, pp. 110β113, 1935. View at Google Scholar - M. A. Biot, βGeneral theory of three-dimensional consolidation,β
*Journal of Applied Physics*, vol. 12, no. 2, pp. 155β164, 1941. View at Publisher Β· View at Google Scholar Β· View at Scopus - R. De Boer, βHighlights in the historical development of the porous media theory: toward a consistent macroscopic theory,β
*Applied Mechanics Reviews*, vol. 49, no. 4, pp. 201β261, 1996. View at Google Scholar Β· View at Scopus - E. Detournay and H.-D. A. Cheng, βFundamentals of poroelasticity,β in
*Comprehensive Rock Engineering: Principles, Practice & Projects*, J. A. Hudson, Ed., pp. 113β171, Pergamon, Oxford, UK, 1993. View at Google Scholar - R. S. Sandhu and E. L. Wilson, βFinite element analysis of seepage in elastic media,β
*Journal of the Engineering Mechanics Division of the American Society of Civil Engineering*, vol. 95, no. 3, pp. 641β652, 1969. View at Google Scholar - E. Detournay and ,A. H.-D. Cheng, βFundamentals of poroelasticity,β in
*Comprehensive Rock Engineering: Principles, Practice and Projects, Vol. II, Analysis and Design Method*, C. Fairhurst, Ed., chapter 5, pp. 113β171, Pergamon, Oxford, UK, 1993. View at Google Scholar - Y. Abousleiman and S. Ekbote, βSolutions for the inclined borehole in a porothermoelastic transversely isotropic medium,β
*Journal of Applied Mechanics*, vol. 72, no. 1, pp. 102β114, 2005. View at Publisher Β· View at Google Scholar Β· View at Scopus - P. Y. P. Cheg, βAxismmetric thermal stresses in an anisotropic finite hollow cylinder,β
*Journal of Thermal Stresses*, vol. 6, no. 2–4, pp. 197β205, 1983. View at Google Scholar - B. Bai, βFluctuation responses of saturated porous media subjected to cyclic thermal loading,β
*Computers and Geotechnics*, vol. 33, no. 8, pp. 396β403, 2006. View at Publisher Β· View at Google Scholar Β· View at Scopus - Y. Wang and E. Papamichos, βAn analytical solution for conductive heat flow and the thermally induced fluid flow around a wellbore in a poroelastic medium,β
*Water Resource Reserch*, vol. 36, no. 5, pp. 3375β3384, 1994. View at Google Scholar - Y. Wang and E. Papamichos, βThermal effects on fluid flow and hydraulic fracturing from wellbores and cavities in low-permeability formations,β
*International Journal for Numerical and Analytical Methods in Geomechanics*, vol. 23, no. 15, pp. 1819β1834, 1999. View at Publisher Β· View at Google Scholar Β· View at Scopus - A. Ghassemi and Q. Tao, βInfluence of coupled chemo-poro-thermoelastic processes on pore pressure and stress distributions around a wellbore in swelling shale,β
*Journal of Petroleum Science and Engineering*, vol. 67, no. 1-2, pp. 57β64, 2009. View at Publisher Β· View at Google Scholar Β· View at Scopus - B. Wirth and I. Sobey, βAn axisymmetric and fully
_{3}D poroelastic model for the evolution of hydrocephalus,β*Mathematical Medicine and Biology*, vol. 23, no. 4, pp. 363β388, 2006. View at Publisher Β· View at Google Scholar Β· View at PubMed Β· View at Scopus - D. Yang and Z. Zhang, βPoroelastic wave equation including the Biot/squirt mechanism and the solid/fluid coupling anisotropy,β
*Wave Motion*, vol. 35, no. 3, pp. 223β245, 2002. View at Publisher Β· View at Google Scholar Β· View at Scopus - A. Arora and S. K. Tomar, βElastic waves along a cylindrical borehole in a poroelastic medium saturated by two immiscible fluids,β
*Journal of Earth System Science*, vol. 116, no. 3, pp. 225β234, 2007. View at Publisher Β· View at Google Scholar Β· View at Scopus - Y. Hamiel, V. Lyakhovsky, and A. Agnon, βCoupled evolution of damage and porosity in poroelastic media: theory and applications to deformation of porous rocks,β
*Geophysical Journal International*, vol. 156, no. 3, pp. 701β713, 2004. View at Google Scholar Β· View at Scopus - A. Ghassemi, βStress and pore prepressure disterbution around apressurized ,cooled crack in holw permeability rock,β in
*Proceedings of the 32nd Workshop on Geothermal Reservoir Engineering Stanford University*, SGP-TR-183, Stanford University, Stanford, Calif, USA, January 2007. - H. M. Youssef, βTheory of generalized porothermoelasticity,β
*International Journal of Rock Mechanics & Mining Sciences*, vol. 44, no. 2, pp. 222β227, 2007. View at Publisher Β· View at Google Scholar Β· View at Scopus - S. Jourine, P. P. Valkoo, and A. K. Kronenberg , βModelling poroelastic hollow cylinderexperiments with realistic boundary conditions,β
*International Journal for Numerical and Analytical Methods in Gwomechanics*, vol. 28, no. 12, pp. 1189β1205, 2004. View at Publisher Β· View at Google Scholar - M. P. Lutz and R. W. Zimmerman, βThermal stresses and effective thermal expansion coefficient of a functionally gradient sphere,β
*Journal of Thermal Stresses*, vol. 19, no. 1, pp. 39β54, 1996. View at Google Scholar Β· View at Scopus - R. W. Zimmerman and M. P. Lutz, βThermal stresses and thermal expansion in a uniformly heated functionally graded cylinder,β
*Journal of Thermal Stresses*, vol. 22, no. 2, pp. 177β188, 1999. View at Google Scholar Β· View at Scopus - M. Jabbari, S. Sohrabpour, and M. R. Eslami, βGeneral solution for mechanical and thermal stresses in a functionally graded hollow cylinder due to nonaxisymmetric steady-state loads,β
*Journal of Applied Mechanics*, vol. 79, pp. 493β497, 2002. View at Google Scholar - R. Poultangari, M. Jabbari, and M. R. Eslami, βFunctionally graded hollow spheres under non-axisymmetric thermo-mechanical loads,β
*International Journal of Pressure Vessels and Piping*, vol. 85, no. 5, pp. 295β305, 2008. View at Publisher Β· View at Google Scholar Β· View at Scopus - M. Shariyat, S. M. H. Lavasani, and M. Khaghani, βNonlinear transient thermal stress and elastic wave propagation analyses of thick temperature-dependent FGM cylinders, using a second-order point-collocation method,β
*Applied Mathematical Modelling*, vol. 34, no. 4, pp. 898β918, 2009. View at Publisher Β· View at Google Scholar Β· View at Scopus - C. F. Lü, W. Q. Chen, and C. W. Lim, βElastic mechanical behavior of nano-scaled FGM films incorporating surface energies,β
*Composites Science and Technology*, vol. 69, no. 7-8, pp. 1124β1130, 2009. View at Publisher Β· View at Google Scholar Β· View at Scopus - A. M. Afsar and H. Sekine, βInverse problems of material distributions for prescribed apparent fracture toughness in FGM coatings around a circular hole in infinite elastic media,β
*Composites Science and Technology*, vol. 62, no. 7-8, pp. 1063β1077, 2002. View at Publisher Β· View at Google Scholar Β· View at Scopus - D. G. Zhang and Y. H. Zhou, βA theoretical analysis of FGM thin plates based on physical neutral surface,β
*Computational Materials Science*, vol. 44, no. 2, pp. 716β720, 2008. View at Publisher Β· View at Google Scholar Β· View at Scopus - S. A. Fazelzadeh and M. Hosseini, βAerothermoelastic behavior of supersonic rotating thin-walled beams made of functionally graded materials,β
*Journal of Fluids and Structures*, vol. 23, no. 8, pp. 1251β1264, 2007. View at Publisher Β· View at Google Scholar Β· View at Scopus - Y. Ootao and Y. Tanigawa, βTransient thermoelastic problem of functionally graded thick strip due to non uniform heat supply,β
*Composite Structures*, vol. 63, no. 2, pp. 139β146, 2004. View at Google Scholar - M. Jabbari, S. Sohrabpour, and M. R. Eslami, βMechanical and thermal stresses in a functionally graded hollow cylinder due to radially symmetric loads,β
*International Journal of Pressure Vessels and Piping*, vol. 79, no. 7, pp. 493β497, 2002. View at Publisher Β· View at Google Scholar Β· View at Scopus - M. Farid, P. Zahedinejad, and P. Malekzadeh, βThree-dimensional temperature dependent free vibration analysis of functionally graded material curved panels resting on two-parameter elastic foundation using a hybrid semi-analytic, differential quadrature method,β
*Materials and Design*, vol. 31, no. 1, pp. 2β13, 2010. View at Publisher Β· View at Google Scholar Β· View at Scopus - A. Bagri and M. R. Eslami, βGeneralized coupled thermoelasticity of functionally graded annular disk considering the Lord-Shulman theory,β
*Composite Structures*, vol. 83, no. 2, pp. 168β179, 2008. View at Publisher Β· View at Google Scholar Β· View at Scopus - B. A. Shariat and M. R. Eslami, βBuckling of thick functionally graded plates under mechanical and thermal loads,β
*Composite Structures*, vol. 78, no. 3, pp. 433β439, 2007. View at Publisher Β· View at Google Scholar Β· View at Scopus - M. Jabbari, A. Bahtui, and M. R. Eslami, βAxisymmetric mechanical and thermal stresses in thick short length FGM cylinder,β
*International Journal of Pressure Vessels and Piping*, vol. 2009, no. 5, pp. 296β306, 86. View at Google Scholar - M. Thieme, K.-P. Wieters, F. Bergner et al., βTitanium powder sintering for preparation of a porous FGM destined as a skeletal replacement implant,β
*Materials Science Forum*, vol. 308–311, pp. 374β380, 1999. View at Google Scholar