Abstract
An improved version of the solid–liquid coupled material point method (MPM) is proposed to simulate ground collapses with fluidization involving transition processes from soil structures to flowing mixture. A water-saturated soil is assumed based on porous media theory, and the physical quantities of the soil and water phases are assigned to two separate sets of particles (material points). The main contribution of this study is the introduction of the fractional step projection method for the time discretization of the momentum equation of the water phase on the assumption of incompressibility. Thanks to this, the proposed solid–liquid coupled MPM is capable of suppressing the pressure oscillations caused by the weak incompressibility of water, which is commonly assumed in the previous studies, and of representing the wide range of behavior of the soil–water mixture at relatively low computational cost. Also, B-spline basis functions are utilized for the spatial discretization to suppress the cell-crossing errors caused by particles crossing element (cell) boundaries. Several numerical tests are conducted to examine the performance of the proposed method that inherits the beneficial features of MPM and demonstrate the capability of reproducing a model experiment of wave collision to sandpile that exhibits the water flow-induced fluidization process of soil involving scouring, transportation and sedimentation.
Similar content being viewed by others
References
Abe K, Soga K, Bandara S (2014) Material point method for coupled hydromechanical problems. J Geotech Geoenviron Eng 140(3):04013033
Bandara S, Ferrari A, Laloui L (2016) Modelling landslides in unsaturated slopes subjected to rainfall infiltration using material point method. Int J Numer Anal Methods Geomech 40(9):1358–1380
Bandara S, Soga K (2015) Coupling of soil deformation and pore fluid flow using material point method. Comput Geotech 63:199–214
Bardenhagen SG, Kober EM (2004) The generalized interpolation material point method. Technical Report 6
Biot MA (1956) Theory of propagation of elastic waves in a fluid-saturated porous solid. II. Higher frequency range. J Acoust Soc Am 28(2):179–191
Bishop AW, Morgenstern N (1960) Stability coefficients for earth slopes. Géotechnique 10(4):129–153
Blanc T, Pastor M (2012) A stabilized fractional step, Runge-Kutta Taylor SPH algorithm for coupled problems in geomechanics. Comput Methods Appl Mech Eng 221–222:41–53
Bowen RM (1982) Compressible porous media models by use of the theory of mixtures. In J Eng Sci 20(6):697–735
Bui HH, Nguyen GD (2017) A coupled fluid-solid SPH approach to modelling flow through deformable porous media. Int J Solids Struct 125:244–264
Bui HH, Sako K, Fukagawa R (2007) Numerical simulation of soil–water interaction using smoothed particle hydrodynamics (SPH) method. J Terramech 44(5):339–346
Chorin AJ (1968) Numerical solution of the Navier–Stokes equations. Math Comput 22(104):745–745
Codina R (2001) Pressure stability in fractional step finite element methods for incompressible flows. J Comput Phys 170(1):112–140
Dunatunga S, Kamrin K (2015) Continuum modelling and simulation of granular flows through their many phases. J Fluid Mech 779:483–513
Ehlers W, Graf T, Ammann M (2004) Deformation and localization analysis of partially saturated soil. Comput Methods Appl Mech Eng 193(27–29):2885–2910
Gao M, Tampubolon AP, Jiang C, Sifakis E (2017) An adaptive generalized interpolation material point method for simulating elastoplastic materials. ACM Trans Graph 36(6):1–12
Griffiths DV, Lane PA (2001) Slope stability analysis by finite elements. Géotechnique 51(7):653–654
Harlow FH, Welch JE (1965) Numerical calculation of time-dependent viscous incompressible flow of fluid with free surface. Phys Fluids 8(12):2182
Higo Y, Oka F, Kimoto S, Morinaka Y, Goto Y, Chen Z (2011) A coupled MPM-FDM analysis method for multi-phase elasto-plastic soils. Soils Found 50(4):515–532
Hughes TJ, Cottrell JA, Bazilevs Y (2005) Isogeometric analysis : CAD, finite elements, NURBS, exact geometry and mesh refinement. Comput Methods Appl Mech Eng 194(39–41):4135–4195
Jassim I, Stolle D, Vermeer P (2013) Two-phase dynamic analysis by material point method. Int J Numer Anal Methods Geomech 37(15):2502–2522
Brackbill JU, Ruppel HM (1986) FLIP: a method for adaptively zoned, particle-in-cell calculations of fluid flows in two dimensions. J Comput Phys 65(2):314–343
Kularathna S, Soga K (2017) Comparison of two projection methods for modeling incompressible flows in MPM. J Hydrodyn 29(3):405–412
Kularathna S, Soga K (2017) Implicit formulation of material point method for analysis of incompressible materials. Comput Methods Appl Mech Eng 313:673–686
Laloui L, Klubertanz G, Vulliet L (2003) Solid–liquid-air coupling in multiphase porous media. Int J Numer Anal Methods Geomech 27(3):183–206
Larese A, Rossi R, Oñate E, Idelsohn SR (2012) A coupled PFEM-Eulerian approach for the solution of porous FSI problems. Comput Mech 50(6):805–819
Li X, Han X, Pastor M (2003) An iterative stabilized fractional step algorithm for finite element analysis in saturated soil dynamics. Comput Methods Appl Mech Eng 192(35–36):3845–3859
Li X, Zhang X, Han X, Sheng DC (2010) An iterative pressure-stabilized fractional step algorithm in saturated soil dynamics. Int J Numer Anal Methods Geomech 34(7):733–753
Mackenzie-Helnwein P, Arduino P, Shin W, Moore JA, Miller GR (2010) Modeling strategies for multiphase drag interactions using the material point method. Int J Numer Methods Eng 83(3):295–322
Morgenstern N (1963) Stability charts for earth slopes during rapid drawdown. Géotechnique 13(2):121–131
Nithiarasu P, Zienkiewicz OC (2006) Analysis of an explicit and matrix free fractional step method for incompressible flows. Comput Methods Appl Mech Eng 195(41–43):5537–5551
Pradhana TA, Teran J, Gast T, Klár G, Fu C, Jiang C (2017) Multi-species simulation of porous sand and water mixtures an improved return mapping for sand plasticity that prevents volume gain artifacts in the traditional Drucker–Prager model. ACM Trans Graph 36(4):105
Sadeghirad A, Brannon RM, Burghardt J (2011) A convected particle domain interpolation technique to extend applicability of the material point method for problems involving massive deformations. Int J Numer Methods Eng 86(12):1435–1456
Sadeghirad A, Brannon RM, Guilkey JE (2013) Second-order convected particle domain interpolation (CPDI2) with enrichment for weak discontinuities at material interfaces. Int J Numer Methods Eng 95(11):928–952
Schrefler BA, Scotta R (2001) A fully coupled dynamic model for two-phase fluid flow in deformable porous media. Comput Methods Appl Mech Eng 190(24–25):3223–3246
Simo JC (1992) Algorithms for static and dynamic multiplicative plasticity that preserve the classical return mapping schemes of the infinitesimal theory. Comput Methods Appl Mech Eng 99(1):61–112
Simo JC, Ortiz M (1985) A unified approach to finite deformation elastoplastic analysis based on the use of hyperelastic constitutive equations. Comput Methods Appl Mech Eng 49(2):221–245
Siriaksorn T, Chi SW, Foster C, Mahdavi A (2018) u-p semi-Lagrangian reproducing kernel formulation for landslide modeling. Int J Numer Anal Methods Geomech 42(2):231–255
Steffen M, Kirby RM, Berzins M (2008) Analysis and reduction of quadrature errors in the material point method (MPM). Int J Numer Anal Methods Geomech 76(6):922–948
Stomakhin A, Schroeder C, Jiang C, Chai L, Teran J, Selle A (2014) Augmented MPM for phase-change and varied materials. ACM Trans Graph 33(4):1–11
Sulsky D, Chen Z, Schreyer HL (1994) A particle method for history-dependent materials. Comput Methods Appl Mech Eng 118(1–2):179–196
Sulsky D, Zhou SJ, Schreyer HL (1995) Application of a particle-in-cell method to solid mechanics. Comput Phys Commun 87(1–2):236–252
Tielen R, Wobbes E, Möller M, Beuth L (2017) A high order material point method. Proced Eng 175:265–272
Wei H, Chen JS, Hillman M (2016) A stabilized nodally integrated meshfree formulation for fully coupled hydro-mechanical analysis of fluid-saturated porous media. Comput Fluids 141:105–115
Yerro A, Alonso E, Pinyol N (2015) The material point method for unsaturated soils. Géotechnique 65(3):201–217
Zhang DZ, Zou Q, VanderHeyden WB, Ma X (2008) Material point method applied to multiphase flows. J Comput Phys 227(6):3159–3173
Zhang F, Zhang X, Sze KY, Lian Y, Liu Y (2017) Incompressible material point method for free surface flow. J Comput Phys 330:92–110
Zhang HW, Wang KP, Chen Z (2009) Material point method for dynamic analysis of saturated porous media under external contact/impact of solid bodies. Comput Methods Appl Mech Eng 198(17–20):1456–1472
Acknowledgements
This work is supported by JSPS KAKENHI Grant Number JP16H02137, “Joint Usage/Research Center for Interdisciplinary Large-scale Information Infrastructures” and “High Performance Computing Infrastructure” in Japan (Project ID: jh180065-NAH). Also, this study partly utilizes computational resources under Collaborative Research Project for Enhancing Performance of Programming provided by Academic Center for Computing and Media Studies, Kyoto University.
Author information
Authors and Affiliations
Corresponding author
Additional information
Publisher's Note
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Rights and permissions
About this article
Cite this article
Yamaguchi, Y., Takase, S., Moriguchi, S. et al. Solid–liquid coupled material point method for simulation of ground collapse with fluidization. Comp. Part. Mech. 7, 209–223 (2020). https://doi.org/10.1007/s40571-019-00249-w
Received:
Revised:
Accepted:
Published:
Issue Date:
DOI: https://doi.org/10.1007/s40571-019-00249-w