This study's objective is to propose a full three-dimensional finite element modeling for examining the buckling behaviors of functionally graded (FG) porous materials with plates, spherical cap, and cylindrical shapes. An improved first-order shear deformation theory (IFSDT) is utilized in order to develop the three-dimensional solid-shell element. The governing equations are developed in a way that offers a parabolic distribution of the transverse shear strains through the FG porous shell thickness and a zero condition of the transverse shear stresses on the top and bottom surfaces. The given IFSDT solid-shell element has the ability to model any kind of shell structures and incorporates a three-dimensional material law. In this study, two types of porosity distribution functions are considered. The performance of the current IFSDT solid-shell element has been demonstrated by examining the buckling behavior of Functionally Graded Porous plates, spherical caps, and cylindrical shells. The effect of the porosity distribution, porosity coefficient, power-law index, geometrical design parameters and boundary conditions on the buckling behaviors of FG porous shells structures is evaluated.