Effect of Material Gradient on Stresses of Thick FGM Spherical Pressure Vessels with Exponentially-Varying Properties

Similar documents
Multivector Functions

On composite conformal mapping of an annulus to a plane with two holes

Technical Report: Bessel Filter Analysis

Advanced Physical Geodesy

Generalized Fibonacci-Lucas Sequence

CHAPTER 5 : SERIES. 5.2 The Sum of a Series Sum of Power of n Positive Integers Sum of Series of Partial Fraction Difference Method

Mapping Radius of Regular Function and Center of Convex Region. Duan Wenxi

( ) 1 Comparison Functions. α is strictly increasing since ( r) ( r ) α = for any positive real number c. = 0. It is said to belong to

Applications of the Dirac Sequences in Electrodynamics

Some Properties of the K-Jacobsthal Lucas Sequence

Modelling rheological cone-plate test conditions

Supplementary materials. Suzuki reaction: mechanistic multiplicity versus exclusive homogeneous or exclusive heterogeneous catalysis

Core Variation in the Entrance Region Flow of Herschel- Bulkley Fluid in an Annuli

KEY. Math 334 Midterm II Fall 2007 section 004 Instructor: Scott Glasgow

A NOTE ON DOMINATION PARAMETERS IN RANDOM GRAPHS

Latticed pentamode acoustic cloak (supplementary Info)

On Some Fractional Integral Operators Involving Generalized Gauss Hypergeometric Functions

DETERMINE VARIATION OF POISSON RATIOS AND THERMAL CREEP STRESSES AND STRAIN RATES IN AN ISOTROPIC DISC

EVALUATION OF SUMS INVOLVING GAUSSIAN q-binomial COEFFICIENTS WITH RATIONAL WEIGHT FUNCTIONS

Range Symmetric Matrices in Minkowski Space

THE ANALYSIS OF SOME MODELS FOR CLAIM PROCESSING IN INSURANCE COMPANIES

Chaoyang University of Technology -- Radial Consolidation -- APPLICATION OF TERZAGHI'S THEORY OF CONSOLIDATION TO PROBLEMS INVOLVING RADIAL FLOW

Lecture 24: Observability and Constructibility

ON CERTAIN CLASS OF ANALYTIC FUNCTIONS

On ARMA(1,q) models with bounded and periodically correlated solutions

Rotational symmetry applied to boundary element computation for nuclear fusion plasma

DANIEL YAQUBI, MADJID MIRZAVAZIRI AND YASIN SAEEDNEZHAD

Non-Linear Bending Analysis of Moderately Thick Functionally Graded Plates Using Generalized Differential Quadrature Method

Generalization of Horadam s Sequence

SHIFTED HARMONIC SUMS OF ORDER TWO

LOCUS OF THE CENTERS OF MEUSNIER SPHERES IN EUCLIDEAN 3-SPACE. 1. Introduction

Progression. CATsyllabus.com. CATsyllabus.com. Sequence & Series. Arithmetic Progression (A.P.) n th term of an A.P.

2012 GCE A Level H2 Maths Solution Paper Let x,

INVERSE CAUCHY PROBLEMS FOR NONLINEAR FRACTIONAL PARABOLIC EQUATIONS IN HILBERT SPACE

Calculation of Matrix Elements in the Foldy-Wouthuysen Representation

Seismic Analysis of an Axial Blower using ANSYS

FAR FIELD SOLUTION OF SH-WAVE BY CIRCULAR INCLUSION AND LINEAR CRACK

THE GRAVITATIONAL POTENTIAL OF A MULTIDIMENSIONAL SHELL

Strong Result for Level Crossings of Random Polynomials

Applicability of Four Parameter Viscoelastic Model for Longitudinal Wave Propagation in Non-Homogeneous Rods

Design and analytically full-wave validation of the invisibility. cloaks, concentrators, and field rotators created with a general

SOME ARITHMETIC PROPERTIES OF OVERPARTITION K -TUPLES

CO-ROTATIONAL MAXWELL FLUID ANALYSIS IN HELICAL SCREW RHEOMETER USING ADOMIAN DECOMPOSITION METHOD

Finite q-identities related to well-known theorems of Euler and Gauss. Johann Cigler

Advanced Higher Formula List

Available online at ScienceDirect. Procedia Engineering 153 (2016 ) 16 23

Lecture 2: Stress. 1. Forces Surface Forces and Body Forces

L8b - Laplacians in a circle

The Application of a Maximum Likelihood Approach to an Accelerated Life Testing with an Underlying Three- Parameter Weibull Model

PEF-5750 Estruturas Leves Ruy Marcelo de Oliveira Pauletti ARGYRIS NATURAL MEMBRANE ELEMENT THE NATURAL FORCE DENSITY METHOD

Using Difference Equations to Generalize Results for Periodic Nested Radicals

Auchmuty High School Mathematics Department Sequences & Series Notes Teacher Version

ON EUCLID S AND EULER S PROOF THAT THE NUMBER OF PRIMES IS INFINITE AND SOME APPLICATIONS

ADDITIONAL INTEGRAL TRANSFORMS

STRESS ANALYSIS OF LAMINATED COMPOSITE CYLINDERS UNDER NON-AXISYMMETRIC LOADING ABSTRACT

a) The average (mean) of the two fractions is halfway between them: b) The answer is yes. Assume without loss of generality that p < r.

Sums of Involving the Harmonic Numbers and the Binomial Coefficients

Conditional Convergence of Infinite Products

A note on random minimum length spanning trees

Some Integral Mean Estimates for Polynomials

Chapter Introduction to Finite Element Methods

ELEMENTARY AND COMPOUND EVENTS PROBABILITY

AN ANALYTICAL SOLUTION OF THE REYNOLDS EQUATION FOR THE FINITE JOURNAL BEARING AND EVALUATION OF THE LUBRICANT PRESSURE

On the Explicit Determinants and Singularities of r-circulant and Left r-circulant Matrices with Some Famous Numbers

ME 354, MECHANICS OF MATERIALS LABORATORY MECHANICAL PROPERTIES AND PERFORMANCE OF MATERIALS: TORSION TESTING*

Strong Result for Level Crossings of Random Polynomials. Dipty Rani Dhal, Dr. P. K. Mishra. Department of Mathematics, CET, BPUT, BBSR, ODISHA, INDIA

2 Governing Equations

IJRET: International Journal of Research in Engineering and Technology eissn: pissn:

physicsandmathstutor.com

Pankaj Thakur 1, Monika Sethi 1, Shivdev Shahi 2, Satya Bir Singh 2, Fadugba Sunday Emmanuel 3

VIII.3 Method of Separation of Variables. Transient Initial-Boundary Value Problems

RECIPROCAL POWER SUMS. Anthony Sofo Victoria University, Melbourne City, Australia.

Solution of a Spherically Symmetric Static Problem of General Relativity for an Elastic Solid Sphere

The Pigeonhole Principle 3.4 Binomial Coefficients

FINITE ELEMENT ANALYSIS OF A BWR FEED WATER DISTRIBUTOR UNDER EXTREME TRANSIENT PRESSURE LOAD

BINOMIAL THEOREM An expression consisting of two terms, connected by + or sign is called a

A Galerkin Finite Element Method for Two-Point Boundary Value Problems of Ordinary Differential Equations

On the Zeros of Daubechies Orthogonal and Biorthogonal Wavelets *

BINOMIAL THEOREM NCERT An expression consisting of two terms, connected by + or sign is called a

MATH /19: problems for supervision in week 08 SOLUTIONS

Inertia Effects in an Externally Pressurized Thrust Bearing Using Herschel - Bulkley Lubricants

VIII.3 Method of Separation of Variables. Transient Initial-Boundary Value Problems

An Exact Solution for Kelvin-Voigt Model Classic Coupled Thermo Viscoelasticity in Spherical Coordinates

Generalized Near Rough Probability. in Topological Spaces

FIXED POINT AND HYERS-ULAM-RASSIAS STABILITY OF A QUADRATIC FUNCTIONAL EQUATION IN BANACH SPACES

W = mgdz = mgh. We can express this potential as a function of z: V ( z) = gz. = mg k. dz dz

1. Using Einstein Summation notation, prove the identity: = A

A dual-reciprocity boundary element method for axisymmetric thermoelastodynamic deformations in functionally graded solids

p-adic Invariant Integral on Z p Associated with the Changhee s q-bernoulli Polynomials

( ) ( ) ( ) ( ) Solved Examples. JEE Main/Boards = The total number of terms in the expansion are 8.

Applied Mathematical Sciences, Vol. 2, 2008, no. 9, Parameter Estimation of Burr Type X Distribution for Grouped Data

PROGRESSION AND SERIES

Ch 3.4 Binomial Coefficients. Pascal's Identit y and Triangle. Chapter 3.2 & 3.4. South China University of Technology

ELASTIC ANALYSIS OF CIRCULAR SANDWICH PLATES WITH FGM FACE-SHEETS

LINEAR PLATE BENDING

This web appendix outlines sketch of proofs in Sections 3 5 of the paper. In this appendix we will use the following notations: c i. j=1.

SOLUTION OF THE RADIAL N-DIMENSIONAL SCHRÖDINGER EQUATION USING HOMOTOPY PERTURBATION METHOD

AIEEE 2004 (MATHEMATICS)

CRACK DETECTION IN EULER-BERNOULLI BEAMS ON ELASTIC FOUNDATION USING GENETIC ALGORITHM BASED ON DISCRETE ELEMENT TECHNIQUE

Minimal order perfect functional observers for singular linear systems

Transcription:

M. Zamai Nejad et al, Joual of Advaced Mateials ad Pocessig, Vol.2, No. 3, 204, 39-46 39 Effect of Mateial Gadiet o Stesses of Thick FGM Spheical Pessue Vessels with Expoetially-Vayig Popeties M. Zamai Nejad a, *, M. Ghaibi a a Depatmet of Mechaical Egieeig, Yasouj Uivesity, Yasouj, Ia. ARTICLE INFO Aticle histoy: Received 5 Ap. 204 Accepted 23 May 204 Available olie 3 Aug. 204 Keywods: Thick sphee Pessue vessel Expoetial Fobeius seies method (FSM) Fiite elemet method (FEM) ABSTRACT Usig the Fobeius seies method (FSM), a aalytical solutio is developed to obtai mechaical stesses of thick spheical pessue vessels made of fuctioally gaded mateials (FGMs). The cylide pessue vessel is subjected to uifom iteal pessue. The modulus of elasticity is gaded alog the adial diectio accodig to powe fuctios of the adial diectio. It is assumed that Poisso s atio is costat acoss the cylide thickess. Pimaily, displacemets ad stesses have bee obtaied as closed-fom solutio. Next, the pofiles ae plotted fo diffeet values of ihomogeeity costat alog the adial diectio. Fially, the poblem was solved by usig the fiite elemet method (FEM). The obtaied esults of fiite elemet method wee compaed with those of the aalytical method. The aalytical solutios ad the solutios caied out though the FEM showed good ageemet. The values used i this study ae abitaily chose to demostate the effect of ihomogeeity o displacemets, ad stesses distibutios.. Itoductio I mateials sciece, fuctioally gaded mateials (FGMs) ae advaced composite mateials whose mechaical popeties vay cotiuously fom oe suface to aothe at maco level. Fom the pespective of solid mechaics, FGMs ae o-homogeeous elastic mediums. Usig this popety i the FGMs, a egiee ca desig composite mateials such that ay potio of the mateials eaches the same safety level []. Tutucu ad Oztuk [2], by usig the ifiitesimal theoy of elasticity, obtaied closed-fom solutios fo stesses ad displacemets i fuctioally gaded cylidical ad spheical vessels subjected to iteal pessue. Assumig that the mateial popeties except Poisso s atio ae vaiable as powe law of adius, Nayak et al. [3] ivestigated displacemets, stais, ad stesses of the FGM spheical vessel. Usig a accuate method, elastic aalysis of iteally pessuized thick-walled spheical pessue vessels of fuctioally gaded mateials was studied by You et al. [4]. Tutucu [5] deived the solutios of FGM thick cylidical shells with expoetially-vayig popeties. Assumig that the mateial popeties ae Coespodig autho: E-mail addess: m_zamai@yu.ac.i (Mohammad Zamai Nejad).

M. Zamai Nejad et al, Joual of Advaced Mateials ad Pocessig, Vol.2, No. 3, 204, 39-46 40 vaiable i ay abitay diectio, ad Poisso s atio is costat, by usig teso aalysis, Zamai Nejad et al. [6] obtaied a complete ad cosistet 3-D set of field equatios of FGM thick shells of evolutio with abitay cuvatue ad vaiable thickess. Che ad Li [7] obtaied the appoximate solutio fo thick cylides ad spheical pessue vessels. Tutucu ad Temel [8] gaied axisymmetic displacemets ad stesses i fuctioally gaded hollow cylides, disks ad sphees subjected to uifom iteal pessue usig plae elasticity theoy ad complemetay fuctios method. Assumig that the displacemet fuctio is ukow i the goveig equatio, Che ad Li [9] obtaied displacemet ad stesses i thick-walled cylides ad sphees made of fuctioally gaded mateials. Zamai Nejad et al. [0] peseted a aalytical solutio ad a umeical solutio fo stesses ad adial displacemet of paabolic FGM solid sphees with paabolic vayig popeties subjected uifom exteal pessue. I aothe study, usig a aalytical method, Zamai Nejad et al. [] obtaied a exact solutio fo stesses ad displacemets of pessuized thick sphees made of fuctioally gaded mateial with expoetiallyvayig popeties. I this pape, by usig the Fobeius seies method, a closed-fom aalytical solutio fo displacemets ad stesses of FGM thick spheical pessue vessels with expoetial vayig popeties has bee obtaied. Fo the umeical solutio, a commecial fiite elemet pogam ANSYS 2 has bee used. 2. Expeimetal Coside a fuctioally gaded hollow sphee with ie adius a ad oute adiusb. The sphee is subjected to iteal uifom pessue P. (Fig. ) The mateial is isotopic ad fuctioally gaded while the Poisso s atio emais costat thoughout the etie sphee. Modulus of elasticity E is assumed to vay expoetial fom as follows: a β E E b a ie [] Whee E i ad β ae modulus of elasticity i ie suface ad ihomogeeity costat, espectively. I the symmetical defomatio case, the equilibium equatio fo the adial ad cicumfeetial stess compoets ( σ ad espectively) disegadig the body foces σθθ ad ietia tems takes the followig fom: dσ σ σ + 2 θθ 0 [2] d Two adial ad cicumfeetial stai compoets ( ε ad ε θθ espectively) ca be expessed as du ε [3] d u ε θθ [4] Hee, the displacemet i the -diectio is deoted by u. The stess ad stai elatios fo ohomogeeous ad isotopic spheical shell ae σ A 2B ε E [5] θθ φφ B A B σ σ + ε θθ εφφ Whee A ad B ae elated to Poisso s atio ν as ( ν) A [6] + ν 2 ν B ( )( ) ν 2 [7] ( + ν)( ν) Usig Eqs. [] to [5], the Navie equatio i tems of the adial displacemet is 2 u ( ) + 2 + h u 2( hν ) u 0 [8] Whee ν ν ν [9] β h b a Eq. [9] ca be solved by usig powe seies method with the solutio i the fom of + s u a [0] 0 Substitutig Eq. [0] ito Eq. [8] s a0 ( s )( s + 2) + ( + s )( + s + 2) a + s + h ( + s + 2ν ) a 0 [] Eq. [] has a idicial equatio ad a ecuece elatio, espectively. Sice a0 0

M. Zamai Nejad et al, Joual of Advaced Mateials ad Pocessig, Vol.2, No. 3, 204, 39-46 4 Fig.. Coss sectio of the FGM thick-walled spheical pessue vessel s + ( s )( s + 2) 0 [2] s2 2 h ( + s + 2ν ) a a [3] ( + s )( + s + 2) The idicial equatio has oots that diffe by a itege, thus oly oe of the solutios is i the fom of Eq. [0]. Expasio of the above ecuece elatio, the coefficiet a i tems of a 0 ad Gamma fuctios ae obtaied. 2 ( ) Γ ( + s + 2ν ) Γ ( s + 3) a0h a s( s + )( s + 2) Γ ( s + 2ν ) Γ ( + s) Γ ( + s + 3) [4] Fo the fist oot ( s +, a0 ), the fist solutio of Eq. [8] obtaied i the fom of ( ) + u + a s Whee 6( ) Γ ( + + 2ν ) a ( s ) h Γ ( 2 ν + ) Γ ( + ) Γ ( + 4) [5] [6] The secod solutio fo s 2 2 will be of the fom 2 2 l ( ) s u Qu + + C ( s2 ) [7] Whee 3 h ν ( )( ) ( ) ( ) ν 2ν Q lim s s s s 2 2 an s 3 d C ( s2 ) ( s s2 ) a ( s) ds Assumig that 2 ( ) Γ ( + s + 2ν ) Γ ( s + 3) ( + ) Γ ( + ν ) Γ ( + ) Γ ( + + ) L h s s s 2 s s 3 hece, [8] [9] L s2 2 ad ( ) Γ ( + 2ν ( ) ( ) Γ( ν ) h 2! 3! 2 2 d l L s s.5 2 ( ) ( ) ( ) Γ '( ) Γ '( + 2ν Γ ( ) Γ ( + ν ) ( ) ( ) ( ) ( ) 2 ds + + 2 2 Γ ' 2ν 2 Γ ' 2 Γ ' + Γ 2ν 2 Γ 2 Γ + It is obvious that Γ '( z) Ψ ( z) Γ z ( ) [20] [2] [22] Whee Ψ(z) is called Psi (Digamma) fuctio. With substitutig Eq. [22] ito Eq. [2] d l ( L ) s s2.5 2 ( ) ( 2 2 ) ds + Ψ + Ψ + ν Ψ 2ν 2 Ψ 2 Ψ + ( ) ( ) ( ) [23] C s will be as Hece, the fial fom of ( 2 ) follows ( ) Γ ( + 2ν h C ( s2 ) 2! ( 3 )! Γ( 2ν.5 + 2Ψ ( ) + Ψ ( + 2ν ( 2 2) ( 2) ( ) Ψ ν Ψ Ψ + [24] Thus, the complete solutio fo u is expessed as + u cu + c2u2 c + a ( s ) + c2 Q + a s l + ( ) ( ) + 2 2 + C ( s2 ) Substitutig Eq. [25] ito Eq. [5] AEc 2 Fa ( s ) σ + ν + ν + AEc2 Q ( + 2ν ) l ( ) + 2 3 [25] 3 + Q Fl ( ) + a ( s ) + GC ( s2 ) [26]

M. Zamai Nejad et al, Joual of Advaced Mateials ad Pocessig, Vol.2, No. 3, 204, 39-46 42 Whee F + + 2ν [27] G 2 + 2ν The costat coefficiets c ad c 2 ae foud usig the bouday coditios σ ( a) P ad σ ( b) 0 PD4 c AEi DD 4 D2D3 PD3 c2 AEi DD4 D2D3 [ ] [ ] Whee D + 2ν + Fa ( s ) a ν D 2 Q ( + 2ν ) l ( a) + 2 a 3 + + + [28] 3 Q Fl ( a) a ( s ) a GC ( s2 ) a 3 + ν + ( ) D 2 Fa s b ν D Q ( + 2ν ) l ( b) + 2 b 4 3 + Q Fl b + a s b + GC s b 3 ( ) ( ) ( 2 ) [29] We have the displacemet ad stesses equatios i the followig fom: P k u ( kq l ( ) ) [ ] AE 2 i D kd2 + 2 + kq l ( ) a ( s ) k C ( s2 ) [30] a β b a Pe σ + 2ν k Q + 2ν l ] ( ) ] ( ) ( ) ( ) ( ) D kd2 ν + 2 F kq Fl 3 + 3 2 + a s k GC s [3] a β b a Pe σ φφ + 2ν k Q + 2ν l D kd 2 ν +ν + 3 + ( ) ( ) ] F kq F l ( ) 3 ] a ( s ) k G C ( s2 ) +ν Whee D3 k D 4 F ν ( + 2) + G v ( ) + [32] [33] 3. Results ad Discussio Coside a thick spheical pessue vessel ude the iteal pessue of 40 MPa. The pessue vessel has the ie ad oute adii 30 cm ad 60 cm, espectively. I this study, it is supposed that the modulus of elasticity at the iteal adius has the value of 200 GPa ad the Poisso s atio, ν, has a costat value of 0.3. Fo a compaative study o umeical aalysis of this poblem, a geometic specime is modeled usig commecial fiite elemet code, ANSYS 2. Due to the geometical symmety i the sphee, oly a quate of the specime geomety i the fiite elemet model was cosideed. I ode to epeset the o-homogeeous specime, a 8-ode axisymmetic quadilateal elemet was used. The vaiatio i mateial popeties was implemeted by 20 layes, with each laye havig a costat value of mateial popeties, fo modelig of the FGM spheical pessue vessel. Fig. 2 shows the distibutio of elastic modulus i the adial diectio. The elastic modulus iceases as β iceases. Fig. 3 shows the distibutio of tesile adial displacemet alog the omalized adial diectio. It is obvious i this cuve that the adial displacemet deceases as β iceases at the same positio. I Fig. 4 distibutio of compessive adial stess alog omalized adial diectio is show. It is peceived that the adial stess

M. Zamai Nejad et al, Joual of Advaced Mateials ad Pocessig, Vol.2, No. 3, 204, 39-46 43 Fig. 2. Radial distibutio of modulus of elasticity Fig. 3. Radial distibutio of adial displacemet Fig. 4. Radial distibutio of adial stess iceases fo highe values of β. I Fig. 5 the tesile cicumfeetial stess alog omalized adial diectio fo diffeet values of β is plotted. Hee, it should be oted that i the same situatio, appoximately, fo a b a <, the value of the ( ) 0.5 cicumfeetial stess deceases as β iceases, wheeas fo ( a b a) > 0.5 this situatio was evesed. Futhe, alog the adial diectio, appoximately, fo β < the cicumfeetial stess deceases, while almost fo β >, the cicumfeetial stess iceases. It should also be oted that fo all the cosideed coditios, fo β the cicumfeetial stess emais almost uifom alog the adius of the sphee. The issue ca be a valuable facto fo cotollig the stess. Fo the pupose of studyig the stess distibutio alog the spheical pessue vessel adius, i Fig. 6 the vo Mises equivalet stess Fig. 5. Radial distibutio of cicumfeetial stess ( σ eq σ σφ ) is plotted i the adial diectio. The vo Mises equivalet stess deceases as the adius iceases fo all β values. I Figs. 7 to 0, the adial displacemets, adial, cicumfeetial ad vo Mises equivalet stesses values obtaied fom ANSYS commecial fiite elemets aalysis pogam ae epeseted. Fo veificatio, the Fobeius seies method esults ae compaed with those of the umeical solutio of Che ad Li method [9] (Figs ad 2 ( a 30 cm, b 60 cm, ν 0.3 )). 4. Coclusios The Fobeius seies method is a poweful techique fo obtaiig solutios of cetai diffeetial equatios which occu i applicatios. Based o basic equatios of elasticity ad usig FSM, closed-fom solutios have bee deived fo stesses ad the

M. Zamai Nejad et al, Joual of Advaced Mateials ad Pocessig, Vol.2, No. 3, 204, 39-46 44 Fig. 6. Radial distibutio of vo Mises equivalet stess Fig. 7. Radial displacemet obtaied fom ANSYS code i a FGM spheical pessue vessel ( β ) Fig. 8. Radial stess obtaied fom ANSYS code i a FGM spheical pessue vessel ( β ) Fig. 9. Cicumfeetial stess obtaied fom ANSYS code i a FGM spheical pessue vessel ( β ) Fig. 0. Vo Mises equivalet stess obtaied fom ANSYS code i a FGM spheical pessue vessel ( β ) displacemets of thick spheical pessue vessels made of fuctioally gaded mateials with expoetial vayig popeties. Followig Fig.. No dimesioal adial stess i spheical pessue vessel with a ie pessue P this, pofiles ae plotted fo diffeet values of ihomogeeity costat fo the adial displacemet, adial stess, ad cicumfeetial

M. Zamai Nejad et al, Joual of Advaced Mateials ad Pocessig, Vol.2, No. 3, 204, 39-46 45 Fig. 2. No dimesioal tagetial stess i spheical pessue vessel with a ie pessue P stess, as a fuctio of adial diectio. I this study, a umeical solutio usig a commecial fiite elemets code, ANSYS 2, is also peseted. Good ageemet was foud betwee the aalytical solutios ad the solutios caied out though fiite elemet code. The peseted esults show that the ihomogeeity costat has a sigificat ifluece o mechaical behavios of the thick spheical pessue vessels made of fuctioally gaded mateial with expoetial vayig popeties. Refeeces. M. Ghaad, M. Z. Nejad, Elastic solutio of pessuized clamped-clamped thick cylidical shells made of fuctioally gaded mateials, J. Theo. Appl. Mech., Vol. 5, 203, pp. 067-079. 2. N. Tutucu, M. Oztuk, Exact solutio fo stesses i fuctioally gaded pessue vessels, Compos. Pat B-Eg., Vol. 32, 200, pp. 683-686. 3. P. Nayak, S. C. Modal, A. Nadi, Stess, stai ad displacemet of a fuctioally gaded thick spheical vessel, It. J. Eg. Sci., Vol. 3, 200, pp. 2659-267. 4. L. H. You, J. J. Zhag, J. J. You, Elastic aalysis of iteally pessuized thick-walled spheical pessue vessels of fuctioally gaded mateials, It. J. Pes. Ves. Pip., Vol. 82, 2005, pp. 347-354. 5. N. Tutucu, Stesses i thick-walled FGM cylides with expoetially vayig popeties, Eg. Stuct., Vol. 29, 2007, pp. 2032-2035. 6. M. Z. Nejad, G. H. Rahimi, M. Ghaad, Set of field equatios fo thick shell of evolutio made of fuctioally gaded mateials i cuviliea coodiate system, Mechaika, Vol. 77, 2009, pp. 8-26. 7. Y. Z. Che, X. Y. Li, Elastic aalysis fo thick cylides ad spheical pessue vessels made of fuctioally gaded mateials, Comp. Mate. Sci., Vol. 44, 2008, pp. 58-587. 8. N. Tutucu, B. Temel, A ovel appoach to stess aalysis of pessuized FGM cylides, disks ad sphees, Compos. Stuct., Vol. 9, 2009, pp. 385-390. 9. Y. Z. Che, X. Y. Li, A alteative umeical solutio of thick-walled cylides ad sphees made of fuctioally gaded mateials, Comp. Mate. Sci., Vol. 48, 200, pp. 640-647. 0. M. Z. Nejad, M. Abedi, M. H. Lotfia, M. Ghaad, Exact ad umeical solutios fo stesses i pessuized FGM solid sphee with paabolic vayig popeties, Am. J. Sci. Res., Vol. 32, 20, pp. 82-89.. M. Z. Nejad, M. Abedi, M. H. Lotfia, M. Ghaad, A exact solutio fo stesses ad displacemets of pessuized FGM thick-walled spheical shells with expoetial-vayig popeties, J. Mech. Sci. Techol., Vol. 26, 202, pp. 408-4087.

M. Zamai Nejad et al, Joual of Advaced Mateials ad Pocessig, Vol.2, No. 3, 204, 39-46 46