Abstract:The inversion of a mass matrix is often involved in solving partial differential equations by finite element method. Mass lumping can be used to improve the computational efficiency by diagonalizing mass matrix with special numerical integration. In this paper, a fully discrete hybrid stress quadrilateral finite element method for the wave propagation in viscoelastic solid media is presented. To realize the diagonalization of the mass matrix, the Gauss-Lobatto numerical integration with displacement interpolation nodes as quadrature nodes is used. A numerical examples is given to verify the performance of the method.