Unlike other meshless methods, the point interpolation method (PIM) adopts the polynomial approximation to construct the shape functions. The created shape functions possess the properties of the Kronecker delta function, which allows easy imposition of