In this research, a thick hollow cylindrical shell made of bidirectional functionally graded open cell porous materials under internal thermal shock according to the classical theory of linear thermo-elasticity is examined for the first time. The cylinder is made of a porous cellular material and its porosity varies along both radial and axial directions continuously. The governing motion equations are obtained by using 2D-axisymmetric theory of linear thermo-elasticity rather than shell theories. This theory represents thickness stretching and gives more precise results. Graded finite element method is employed to model the problem. Applying this method rather than conventional FEM leads to more accurate results, especially for dynamic analyses. The cubic higher order element is used for dividing the solution domain. To obtain transient temperatures, Crank–Nicolson algorithm is used and then the Newmark procedure is used to derive time responses of displacement and stress components. The time history of displacements and stress components for different radial and axial power law exponents, porosity coefficient, boundary conditions, length-to-thickness ratio and two different porosity patterns are investigated in detail. The obtained results show that thermal-induced vibration is generally caused by hoop stress, and frequency and amplitude of vibrations and velocity of stress waves are considerably influenced by the porosity distribution, porosity coefficient and power law exponents in both directions.