SUMMARYAn element-free Galerkin (EFG) method with linear, quadratic and cubic approximations, which can exactly, in a numerical sense, pass the corresponding patch tests is proposed and is named as consistent EFG (CEFG) method. The development of this method is based on the Hu-Washizu three-field variational principle. Numerical integration schemes with corrected nodal derivatives at quadrature points are proposed according to the satisfaction of the orthogonality condition between stress and strain difference. Thus, the method is variationally consistent. The consistency of the corrected nodal derivatives and the satisfaction of patch test conditions are theoretically proved and also numerically validated. Numerical results show that the proposed CEFG method greatly improves the numerical performance of the EFG method in terms of accuracy, convergence, efficiency and stability, especially the proposed cubic CEFG method, which shows exceptional accuracy and convergence.