ESPFEM2D: A MATLAB 2D explicit smoothed particle finite element method code for geotechnical large deformation analysis

Abstract

The Smoothed Particle Finite Element Method (SPFEM) has gained popularity as an effective numerical method for modelling geotechnical problems involving large deformations. To promote the research and application of SPFEM in geotechnical engineering, we present ESPFEM2D, an open-source two-dimensional SPFEM solver developed using MATLAB. ESPFEM2D discretizes the problem domain into computable particle clouds and generates the finite element mesh using Delaunay triangulation and the \( \alpha \)-shape technique to resolve mesh distortion issues. Additionally, it incorporates a nodal integration technique based on strain smoothing, effectively eliminating defects associated with the state variable mapping after remeshing. Furthermore, the solver adopts a simple yet robust approach to prevent the rank-deficiency problem due to under-integration by using only nodes as integration points. The Drucker-Prager model is adopted to describe the soil’s constitutive behavior as a demonstration. Implemented in MATLAB, this open-source solver ensures easy accessibility and readability for researchers interested in utilizing SPFEM. ESPFEM2D can be easily extended and effectively coupled with other existing codes, enabling its application to simulate a wide range of large geomechanical deformation problems. Through rigorous validation using four numerical examples, namely the oscillation of an elastic cantilever beam, non-cohesive soil collapse, cohesive soil collapse, and slope stability analysis, the accuracy, effectiveness and stability of this open-source solver have been thoroughly confirmed.