This paper reformulates a time-discontinuous finite element method (TD-FEM) based on a new class of shape functions, called complex Fourier hereafter, for solving two-dimensional elastodynamic problems. These shape functions, which are derived from their corresponding radial basis functions, have some advantages such as the satisfaction of exponential and trigonometric function fields in complex space as well as the polynomial ones simultaneously, that make them a better choice than classic Lagrange shape functions, which only can satisfy polynomial function field. To investigate the validity and accuracy of the proposed method, three numerical examples are provided and the results obtained from the present method (complex Fourier-based TD-FEM) and the classic Lagrange-based TD-FEM are compared with the exact analytical solutions. According to them, using complex Fourier functions in TD-FEM leads to more accurate and stable solutions rather than those obtained from the classic TD-FEM.