Modeling irregular boundaries using isoparametric elements in the Material Point Method

Ezra Y. S. Tjung,Shyamini Kularathna,Krishna Kumar,Kenichi Soga

DOI: https://doi.org/10.48550/arXiv.1909.13389

2019-09-30

Abstract:The Material Point Method (MPM) is a hybrid Eulerian-Lagrangian approach capable of simulating large deformation problems of history-dependent materials. While the MPM can represent complex and evolving material domains by using Lagrangian points, boundary conditions are often applied to the Eulerian nodes of the background mesh nodes. Hence, the use of a structured mesh may become prohibitively restrictive for modeling complex boundaries such as a landslide topography. We study the suitability of unstructured background mesh with isoparametric elements to model irregular boundaries in the MPM. An inverse mapping algorithm is used to transform the material points from the global coordinates to the local natural coordinates. Dirichlet velocity and frictional boundary conditions are applied in the local coordinate system at each boundary node. This approach of modeling complex boundary conditions is validated by modeling the dynamics of a gravity-driven rigid block sliding on an inclined plane. This method is later applied to a flume test of controlled debris flow on an inclined plane conducted by the United States Geological Survey (USGS).

Geophysics

What problem does this paper attempt to address?

This paper aims to solve the problem of modeling irregular boundaries in the Material Point Method (MPM). Traditional MPM methods use rectangular structured grids to represent the background grid, which is very limited when dealing with complex terrains such as geological disasters like landslides. Because the boundary conditions of these natural terrains are usually not aligned with the global coordinate axes, directly applying traditional methods will lead to inaccurate results and even produce unrealistic flow constraints.

To overcome this challenge, the author proposes a new method, that is, using unstructured quadrilateral isoparametric elements to adapt to irregular boundaries. The keys of this method are as follows:

1. **Boundary Conditions in the Local Coordinate System**: Since the constraints on the irregular boundary surface are usually not aligned with the global coordinate axes, it is necessary to define the boundary conditions of velocity and friction in the local coordinate system of each node. This involves the application of the transformation matrix from the global coordinate system to the local coordinate system.

2. **Inverse Mapping Algorithm**: In MPM, the positions of material points are used as integration points, rather than the Gaussian integration points used in the Finite Element Method (FEM). This means that the properties of material points (such as mass, momentum, stress, etc.) need to be mapped from the natural coordinate system to the corresponding node properties. For irregular elements, this inverse mapping is not directly defined, so the author adopts the general solution proposed by Zhao et al. (1999) and solves this problem in matrix form.

Through this method, the author verifies its effectiveness in simulating the sliding of a gravity - driven rigid body along an inclined plane and the dry sand flow in a control experiment. The results show that this method can accurately capture the dynamic behaviors of rigid bodies and granular flows and is highly consistent with the experimental results. This provides a solid foundation for the wider application of irregular boundary conditions in MPM in the future, especially when simulating complex geological phenomena such as landslides.