Abstarct: In this research, the governing equation of heterogeneous axisymmetric thick- walled piezoelectric sphere subjected to electrical, thermal and mechanical loading baxsed on the plane elasticity theory is presented. The material properties are assumed to vary nonlinearly in the radial direction. The temperature distribution is assumed to be a function of radius. Using FEM simulations, the analytical findings for FGPM spheres under the influence of internal pressure, electric gradient and electro-thermo-mechanical loads were compared to finite element results. It was observed that the two methods results are in good agreement. It is concluded that by changing the value of n_i, the stresses and displacement distribution of¬ FGM an be so modified that the lowest stress levels are reached.