Abstract：The system equation for 2D soil consolidation is given and the corresponding program is developed with element-free Galerkin method based on the equilibrium differential equation and the continuous equation of soil，combined with the initial and boundary conditions. In addition，the essential boundary conditions are enforced with penalty method. This method forms the shape functions using the moving least-squares(MLS) polynomial interpolation grounded on a set of arbitrarily distributed points in an influencing domain. Time discretization adopts the Crank-Nicolson integral scheme. The accuracy and feasibility of the EFG method for consolidation problem are verified by comparing the calculation results of a numerical example with those from the analytical method and FEM.
