An implicit 3D nodal integration based PFEM (N-PFEM) of natural temporal stability for dynamic analysis of granular flow and landslide problems

Abstract

The particle finite element method (PFEM) is a robust approach for modelling large deformation problems with free surface evolution. The classical PFEM, however, requires variable mapping from old to new quadrature points when adopting history-dependent material models in granular flow and landslide problems. Although the nodal integration technique circumvents this issue, it makes the PFEM temporal instable in dynamic analysis when using a displacement-based formulation. In this study, we developed a new version of a three-dimensional (3D) Nodal integration based PFEM (N-PFEM) using a mixed variational principle with the final problem resolved in mathematical programming. The proposed N-PFEM not only inherits the benefit from the nodal integration scheme that no variable mapping is required for handling history-dependent models but also is naturally temporal stable requiring no ad-hoc stabilization technique. We simulated a series of benchmark problems to demonstrate its nature of temporal stability as well as other admirable features such as the volumetric-locking free property and capability for tackling extreme configuration changes. Additionally, its application to a 3D landslide with a sensitive clay layer is shown to highlight its robustness.