Abstarct: On the basis of plane elasticity theory (PET), the governing equations for axisymmetric thick cylindrical shells made of nonhomogeneous functionally graded piezoelectric materials (FGPMs) subjected to electrical, thermal and mechanical loading are derived. The material properties are assumed to vary nonlinearly in the radial direction. The temperature distribution is assumed to be a function of radius. The obtained result (radial stress, meridional stress, axial stress ,Von Mises stress, radial displacement and electrical potential) are presented in the form of graph. The problem was also solved, using the finite element method (FEM), the results of which were compared with those of the analytical method. It was observed that the two method’s results are in good agreement.