This study aims to investigate nonlinear dynamic behavior of sandwich beams constituted by functionally graded material (FGM) faces and functionally graded porous (FGP) core under the action of multiple moving loads. The energy equations based on a quasi-3D beam theory considering several functions of shear deformation were established for producing the nonlinear equations of motion used for describing dynamic responses of beams with various boundary conditions. By using the hybrid processes of the Newton-Raphson iteration procedure, the time-integration approach of Newmark-β, and the Gram-Schmidt-Ritz method, the new results of nonlinear analysis were explored and proposed. Several prime parameters such as the porous coefficient, sandwich thickness ratio, power law index and others that significantly affect the dynamic response of the beams were taken into investigations. It can be concluded that beams under one moving load have larger dynamic deflections compared to that under many moving loads. This is because the distance between the moving loads can produce a bending moment against the beam deformation.