Large Deformation Modeling of Trapdoor Test Using MLS-MPM

Barani, Omid R.; Maleki, Amir Mohammad

doi:10.61882/NMCE.2506.1093

Large Deformation Modeling of Trapdoor Test Using MLS-MPM

Document Type : Research

Authors

Omid R. Barani

Amir Mohammad Maleki

Department of Civil Engineering, K. N. Toosi University of Technology, Tehran, Iran.

10.61882/NMCE.2506.1093

Abstract

Moving least squares material point method (MLS-MPM) is a new method to allow for sharp separation of particles and to reduce the cell-crossing error in the material point method simulations. This paper uses MLS-MPM method to numerically simulate the trapdoor test and soil arching, a phenomenon which plays a vital role in the design of many geotechnical structures, such as pile-supported embankments, buried pipelines, tunnels, and trench excavations. Large deformation of soil is considered during the plastic phase of deformation. Non-associated Mohr-Coulomb plasticity is enhanced to capture strain hardening/softening typically seen in granular soils. The proposed method is validated against centrifuge trapdoor tests. The numerical stress fields agree with the formation of triangular shaped arch after medium relative displacements and mobilization of a prism of soil after large relative displacements of the trapdoor. It is shown that using MLS-MPM, the arching behavior of soil can be described in detail throughout the trapdoor test.

Keywords

Arching

Trapdoor test

Moving Least Squares Material Point Method

Granular media

Subjects

Geotechnical Engineering

[1] Terzaghi, K., 1936. Stress distribution in dry and saturated sand above a yielding trapdoor. Proceedings 1st Conference on Soil Mechanics and Foundations Engineering, Graduate School of Engineering. Harvard University, Cambridge, MA, 307-311.

[2] Costa, Y. D., Zornberg, J.G., Bueno, B.S., and Costa, C. L., 2009. Failure mechanisms in sand over a Deep active trapdoor. Journal Geotechnical and Geoenvironmental Engineering. 10.1061/(ASCE)GT.1943-5606.0000134, 1741-1753.

[3] Engesser, Fr., 1882. Ueber den Erdduck gegen innere Stüzwande (Tunnelwande). Deutsche Bauzeitung, 16, 91-93 (in German).

[4] Iglesia G. R., 1991. Trapdoor experiments on the centrifuge: a study of arching in geomaterials and similitude in geotechnical models. Ph.D. thesis, Massachusetts Institute of Technology, Cambridge.

[5] Dewoolkar, M. M., Santichaianant, K., and Ko, H.-Y., 2007. Centrifuge modelling of granular soil response over active circular trapdoors. Soils and Foundations 47(5),931-945.

[6] Cundall, P. A., Strack, O. D. L., 1979. A discrete numerical model for granular assemblies. Geotechnique 29, 47-65.

[7] Mora, P., Place, D., 1993. A lattice solid model for the non-linear dynamics of earthquakes. International Journal of Modern Physics. C 4, 1059-1074.

[8] Radjaï, F., Dubois, F., 2011. Discrete-element Modeling of Granular Materials. Wiley-Iste, Hoboken, New Jersey, USA.

[9] Tu, X. and Andrade, J. E., 2008. Criteria for static equilibrium in particulate mechanics computations, International Journal for Numerical Methods in Engineering 75, 1581–1606.

[10] Gray G. G., Morgan J. K., and Sanz P. F., 2014. Overview of continuum and particle dynamics methods for mechanical modeling of contractional geologic structures. Journal of Structural Geology. 59, 19–36.

[11] Lim, K.-W. Krabbenhoft, K., and Andrade, J.E., 2014. A contact dynamics approach to the Granular element method, Computer Methods in Applied Mechanics and Engineering 268, 557–573.

[12] Seyyedan, S. M., Mirghasemi, A. A., Mohammadi, S., 2021. Numerical simulation of direct shear test on granular materials composed of breakable angular particles: A DEM-XFEM approach. Powder Technology 391, 450-466.

[13] Andersen, S. and Andersen, L., 2010. Modelling of landslides with the material-point method. Computational Geosciences 14 (1), 137-147.

[14] Barani, O., and Bagherzadeh, H., 2023. Simulation of the Vaiont landslide via multi-body material point method with cohesive frictional interfaces. Scientia Iranica 30, 902-917.

[15] Zhang, W., Wu, Z., Peng, C., Li, S., Dong, Y., and Yuan, W., 2023. Modelling large-scale landslide using a GPU-accelerated 3D MPM with an efficient terrain contact algorithm. Computer and Geotechnics 158, 105411

[16] Coetzee, C., J., Vermeer, P., A., and Basson, A., H., 2005. The modelling of anchors using the material point method. International Journal for Numerical and Analytical Methods in Geomechanics 29 (9), 879-895.

[17] Fetrati, M., Galavi, V., Goodarzi, M., Mörz, T., and Kreiter, S., 2024. Simulation of cone penetrometer tests in sand using three advanced constitutive models: a comparative study. Computers and Geotechnics 176, 106683.

[18] Galavi, V., Martinelli, M., 2024. MPM simulation of the installation of an impact driven pile in dry sand and subsequent axial bearing capacity. Journal of Geotechnical and Geoenvironmental Engineering 150 (4), 04024019.

[19] Alonso, E. and Zabala, F., 2011. Progressive failure of Aznalc_ollar dam using the material point method. Geotechnique 61 (9), 795-808.

[20] Talbot, L., E., D, Given, J., Tjung, E., Y., S., Liang, Y., Chowdhury, K., Seed, R., and Soga, K., 2024. Modeling large-deformation features of the Lower San Fernando Dam failure with the material point method. Computers and Geotechnics 165(2), 105881.

[21] Hu, Y., Fang, Y., Ge, Z., Qu, Z., Zhu, Y., Pradhana, A., Jiang, C., 2018. A moving least squares material point method with displacement discontinuity and two-way rigid body coupling. ACM Transactions on Graphics. 37 (4), 1–14.

[22] Song, J., U., and Kim, H., G., 2021. An improved material point method using moving least square shape functions. Computational Particle Mechanics. 8(4), 751–766.

[23] Nakamura, K., Matsumura, S. and Mizutani, T., 2021 Particle-to-surface frictional contact algorithm for material point method using weighted least squares Computers and Geotechnics 134, 104069.

[24] Zhao, S., He, S., Li, X., Scaringi, G., Liu, Y., Deng, Y., 2024. Investigating the dynamic process of a rock avalanche through an MLS-MPM simulation incorporated with a nonlocal μ(I) rheology model. Landslides. 21, 1483–1499.

[25] Zhao, S., Liu, Y., Scaringi, G., Li, X., He, S., He, G., Zhu, L., 2025. Non-local μ(I) rheology improves landslide deposition modeling in MLS-MPM simulations. Engineering Geology. 348, 107963.

[26] Belytschko, T., Lu, Y., and Gu, L., 1994. Element-free Galerkin methods. International Journal for Numerical Methods in Engineering 37(2), 229–256.

[27] Liu, W., K., Li S. and Belytschko, T., 1997. Moving least square reproducing kernel methods Part I: Methodology and convergence. Computer Methods in Applied Mechanics and Engineering. 143(1-2), 113-154.

[28] Huerta, A., and Fernandez-Mendez, S., 2000. Enrichment and coupling of the finite element and meshless methods. International Journal for Numerical Methods in Engineering. 48(11), 1615-1636.

[29] Jiang, C., Schroeder, C., Selle, A., Teran, J., and Stomakhin, A., 2015. The affine particle-in-cell method, ACM Transactions on Graphics. 34(4), 1-10.

[30] Jiang, C., Schroeder, C., and Teran, J., 2017. An angular momentum conserving affine-particle-in-cell method, Journal of Computational Physics. 338, 137-164.

[31] Simo, J., C., 1988. A framework for finite strain elastoplasticity based on maximum plastic dissipation and the multiplicative decomposition: Part I. continuum formulation. Computer Methods in Applied Mechanics and Engineering, 66(2), 199–219.

[32] Mast, C., 2013. Modeling landslide-induced flow interactions with structures using the Material Point Method. PhD thesis. University of Washington, USA.

[33] Truesdell, C., A., 1991. A First Course in Rational Continuum Mechanics, Academic Press, San Diego, USA.

[34] Stomakhin, A., Howes, R., Schroeder, C., and Teran, J., 2012. Energetically consistent invertible elasticity. In Proceeding of ACM SIGGRAPH/Eurographics Symposium on Computer Animation, 25–32.

[35] Soga, K., Alonso, E., Yerro, A., Kumar, K. and Bandara, S., 2016. Trends in large deformation analysis of landslide mass movements with particular emphasis on the material point method. Géotechnique 66(3), 248–273.

[36] Payan, M., Fathipour, H., Hosseini, M., Chenari, R. J., and Shiau, J. S. 2022. Lower bound finite element limit analysis of geo-structures with non-associated flow rule. Computers and Geotechnics, 147, 104803.

[37] Hosseini, M., B. Yaghoobi, P. Z. Ranjbar, M. Payan, and R. J. Chenari. 2023.Stability analysis of geosynthetic-reinforced shallow foundations with a lower-bound FELA approach adopting the nonassociated plastic flow rule. International Journal of Geomechanics, 23(2), 04022280.

[38] Bazant, Z.P., Oh, B.-H., 1983. Crack band theory for fracture of concrete. Materials and Structures 16, 155–177.

[39] Iglesia, G., R., Einstein, H., H., and Whitman, R., V., 2011. Validation of centrifuge model scaling for soil systems via trapdoor tests. Journal of Geotechnical and Geoenvironmental Engineering. 137 (11), 1075-1089.

[40] Iglesia, G., R., Einstein, H., H., and Whitman, R., V., 2013. Investigation of soil arching with centrifuge tests. Journal of Geotechnical and Geoenvironmental engineering. 140(2), 04013005.

[41] Barandiaran Villegas L., B., 2018. Sinkhole development over karstic substratum. An MPM approach. Master's thesis, Universität Politecnica de Catalunya.