In this paper we present a new application for a four variable refined plate theory to analyse the nonlinear cylindrical bending behavior of functionally graded plates subjected to thermomechanical loadings. This recent theory is based on the assumption that the transverse displacements consist of bending and shear components in which the bending components do not contribute toward shear forces and, likewise, the shear components do not contribute toward bending moments. The theory accounts for a quadratic variation of the transverse shear strains across the thickness, and satisfies the zero traction boundary conditions on the top and bottom surfaces of the plate without using shear correction factors. The material properties are assumed to vary continuously through the thickness of the plate according to a power-law distribution of the volume fraction of the constituents. The non-linear strain-displacement relations in the von Karman sense are used to study the effect of geometric non-linearity. The solutions are achieved by minimizing the total potential energy and the results are compared to the classical and the first-order theories reported in the literature. It can be concluded that the proposed theory is accurate and simple in solving the nonlinear cylindrical bending behavior of functionally graded plates.