A stable implicit nodal integration-based particle finite element method (N-PFEM) for modelling saturated soil dynamics

In this study, we present a novel nodal integration-based particle finite element method (N-PFEM) designed for the dynamic analysis of saturated soils. Our approach incorporates the nodal integration technique into a generalised Hellinger-Reissner (HR) variational principle, creating an implicit PFEM formulation. To mitigate the volumetric locking issue in low-order elements, we employ a node-based strain smoothing technique. By discretising field variables at the centre of smoothing cells, we achieve nodal integration over cells, eliminating the need for sophisticated mapping operations after re-meshing in the PFEM. We express the discretised governing equations as a min-max optimisation problem, which is further reformulated as a standard second-order cone programming (SOCP) problem. Stresses, pore water pressure, and displacements are simultaneously determined using the advanced primal-dual interior point method. Consequently, our numerical model offers improved accuracy for stresses and pore water pressure compared to the displacement-based PFEM formulation. Numerical experiments demonstrate that the N-PFEM efficiently captures both transient and long-term hydro-mechanical behaviour of saturated soils with high accuracy, obviating the need for stabilisation or regularisation techniques commonly employed in other nodal integration-based PFEM approaches. This work holds significant implications for the development of robust and accurate numerical tools for studying saturated soil dynamics.