It is difficult to settle essential boundary conditions in meshless method because the meshless approximation functions are not built by interpolation schemes.