The meshless method with Galerkin weak form and radial basis point interpolation method (RPIM) shows good performance in solving partial differential equation (PDE)， but it is difficult to improve the computational efficiency and accuracy at the same time. In order to improve the efficiency of the meshless method， this paper sets up a domain based on the background grid， and uses the same node during the process of integral interpolation within the domain， which can reduce some matrix computations and computing time of the RPIM method. In terms of accuracy， this paper proposes a meshless method of hybrid stress， i.e.， using Hellinger-Reissner (H-R) variational principle to derive the solving equations and using meshless method to solve the problem. Numerical example demonstrates that the new method has the higher computational accuracy and the faster computational speed when calculating two-dimensional solid mechanics， thus it has a high practical application value.