As a boundary-type meshless method

the hybrid boundary node method (HBNM) is based on a modified variational principle and moving least squares (MLS) approximation; so it has the advantages of both BEM and meshless methods. The dual reciprocity method (DRM) is combined with HBNM

in which the solution is divided into particular solution and general solution. The general solution is solved by means of HBNM; and the particular solution is approximated by using radial basis functions. Only randomly distributed nodal points on the bounding surface of the domain are required and it doesn’t need extra equations to compute internal parameters in the domain. The post-processing is very simple. The numerical examples show that this method has high precision and is more suitable for computation in geotechnical engineering.