Since the functionally graded materials (FGMs) are used extensively as thermal barriers in many of applications. Therefore, the current article focuses on studying and presenting dynamic responses of multilayer functionally graded (FG) deep beams placed in a thermal environment that is not addressed elsewhere. The material properties of each layer are proposed to be temperature-dependent and vary continuously through the height direction based on the Power-Law function. The deep layered beam is exposed to harmonic sinusoidal load and temperature rising. In the modelling of the multilayered FG deep beam, the two-dimensional (2D) plane stress continuum model is used. Equations of motion of deep composite beam with the associated boundary conditions are presented. In the frame of finite element method (FEM), the 2D twelve-node plane element is exploited to discretize the space domain through the length-thickness plane of the beam. In the solution of the dynamic problem, Newmark average acceleration method is used to solve the time domain incrementally. The developed procedure is verified and compared, and an excellent agreement is observed. In numerical examples, effects of graduation parameter, geometrical dimension and stacking sequence of layers on the time response of deep multilayer FG beams are investigated with temperature effects.