Local point interpolation method is a new type meshless method. Because its shape functions for displacements and pore water pressures possess delta function property

the essential boundary conditions can be imposed easily. After the basic theory is introduced and the main equations for solving Biot’s consolidation problem are developed

the consolidation and settlement of one multi-layer soft ground subjected to heaped loadings are calculated with it. Furthermore

the famous Duncan-Chang model is adopted in consideration of nonlinear properties of soil

and the uniform point distribution mode is selected. At last

comparisons with the results of one-dimensional consolidation method and FEM show that this method is more accurate and satisfying.