Abstract：Various interface elements are successfully applied in traditional finite element method(FEM) to model discontinuities，in which Goodman element is the most representative one. The paper discusses the feasibility of introducing Goodman element to model geological discontinuities in meshless method in full length，and points out the problems existing in some current studies. The Goodman element is presented in the framework of FEM；the displacement mode of Goodman element is designed to be compatible with the finite element along the common boundary between finite element and Goodman element. But，in the meshless method，the displacement mode of Goodman element generally is not compatible with the meshless displacement mode which is based on a number of discrete nodes，the number，however，cannot be known beforehand. The key to solve this problem is that the stiffness matrix of the Goodman element must be computed through numerical integration，then the computed matrix but not the analytical stiffness matrix in the traditional FEM is added to the general matrix of the system when the contribution of the Goodman element to the general stiffness matrix is considered. So，discontinuities can be represented by Goodman element in analytical model in implicit or explicit mode. In the framework of the natural element method(NEM)，the implementation scheme of the Goodman element，implicitly and explicitly respectively，is addressed in detail for illustration of the idea. The presented method is general and suitable for all existing meshless methods.
