pith. sign in

arxiv: 1907.06990 · v1 · pith:BG4BKZ6Inew · submitted 2019-07-13 · 💻 cs.CE

A stabilized total-Lagrangian SPH method for large deformation and failure in geomaterials

Md Rushdie Ibne Islam , Chong Peng This is my paper

Pith reviewed 2026-05-24 21:32 UTC · model grok-4.3

classification 💻 cs.CE
keywords smoothed particle hydrodynamicstotal Lagrangian formulationgeomaterialslarge deformationhourglass controltensile instabilityfailure modelingplastic flow
0
0 comments X

The pith

Stabilized total-Lagrangian SPH models large deformation and plastic flow in geomaterials while producing smooth accurate stresses.

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

The paper introduces a total-Lagrangian smoothed particle hydrodynamics formulation that is free of tensile instability by construction. It adds a stiffness-based hourglass control to remove rank-deficiency modes and periodically resets the reference configuration so the method can continue through post-failure flow. Conventional Eulerian SPH requires multiple extra stabilizations for the same tasks; the new combination is shown to run without them and to return smoother stress fields. A reader would care because reliable simulation of geomaterial failure under large strain is needed for geotechnical engineering and hazard assessment.

Core claim

The stabilized TLSPH formulation is inherently free of tensile instability; a stiffness-based hourglass control cures rank-deficiency; and periodic reference-configuration updates allow the method to follow large plastic deformations and post-failure flows in geomaterials, yielding accurate and smooth stress results without the auxiliary treatments required by conventional Eulerian SPH.

What carries the argument

Total-Lagrangian SPH kernel with stiffness-based hourglass control and periodic reference-configuration updates.

If this is right

  • TLSPH requires fewer auxiliary numerical treatments than CESPH for the same class of problems.
  • Stress fields remain smooth and accurate even after the onset of plastic flow and failure.
  • The method can continue stable simulation once material has undergone large post-failure motion.
  • Hourglass control and configuration updates are shown to be essential for maintaining robustness.

Where Pith is reading between the lines

These are editorial extensions of the paper, not claims the author makes directly.

  • The same stabilization pattern could be tested on other particle methods that suffer rank deficiency.
  • If the reference updates can be made fully automatic, the approach might become suitable for long-duration run-out simulations.
  • Because stresses are reported smoother, the method may reduce the need for post-processing filtering in engineering output.

Load-bearing premise

The stiffness-based hourglass algorithm removes rank-deficiency modes without creating new artifacts or lowering stress accuracy, and the periodic reference updates preserve the stability of the total-Lagrangian description once flow begins.

What would settle it

A benchmark large-strain shear or collapse test in which the stabilized TLSPH produces visibly oscillatory or non-physical stress fields at the same resolution where conventional SPH with its usual stabilizations remains smooth.

Figures

Figures reproduced from arXiv: 1907.06990 by Chong Peng, Md Rushdie Ibne Islam.

Figure 1
Figure 1. Figure 1: Setup of the soil block for large deformation analysis [PITH_FULL_IMAGE:figures/full_fig_p009_1.png] view at source ↗
Figure 2
Figure 2. Figure 2: Conventional SPH and removal of tensile instability with di [PITH_FULL_IMAGE:figures/full_fig_p010_2.png] view at source ↗
Figure 3
Figure 3. Figure 3: Final deformed shape free from tensile instability with conventional SPH ( [PITH_FULL_IMAGE:figures/full_fig_p010_3.png] view at source ↗
Figure 4
Figure 4. Figure 4: Comparison of the stabilised deformed shape with conventional SPH ( [PITH_FULL_IMAGE:figures/full_fig_p011_4.png] view at source ↗
Figure 5
Figure 5. Figure 5: Comparison of the deformed shape and accumulated plastic strain at di [PITH_FULL_IMAGE:figures/full_fig_p011_5.png] view at source ↗
Figure 6
Figure 6. Figure 6: Effect of update of configurations in TLSPH To investigate the influence of the threshold k, three simulations using TLSPH are performed with k = 2, 4, and 8. The results are also given in [PITH_FULL_IMAGE:figures/full_fig_p012_6.png] view at source ↗
Figure 7
Figure 7. Figure 7: Comparison of the effect of hourglass control force at time step with TLSPH TLSPH with hourglass control TLSPH without hourglass control [PITH_FULL_IMAGE:figures/full_fig_p013_7.png] view at source ↗
Figure 8
Figure 8. Figure 8: Comparison of the stabilised deformed shape in TLSPH with and without hourglass control [PITH_FULL_IMAGE:figures/full_fig_p013_8.png] view at source ↗
Figure 9
Figure 9. Figure 9: Setup of the slope for factor of safety analysis [PITH_FULL_IMAGE:figures/full_fig_p014_9.png] view at source ↗
Figure 10
Figure 10. Figure 10: Distribution of initial stress The material constants and the numerical settings are given in [PITH_FULL_IMAGE:figures/full_fig_p014_10.png] view at source ↗
Figure 11
Figure 11. Figure 11: Displacement over time at different factor of safety [PITH_FULL_IMAGE:figures/full_fig_p015_11.png] view at source ↗
Figure 12
Figure 12. Figure 12: Displacement of the slope at different time steps with fs = 1.8 [PITH_FULL_IMAGE:figures/full_fig_p015_12.png] view at source ↗
Figure 13
Figure 13. Figure 13: Plastic strain of the slope at different time steps with fs = 1.8 15 [PITH_FULL_IMAGE:figures/full_fig_p015_13.png] view at source ↗
Figure 14
Figure 14. Figure 14: Stress (σzz) distribution of the slope at different time steps with fs = 1.8 The developments of displacement and shear band of the slope under reduction factor f = 1.8 are given in [PITH_FULL_IMAGE:figures/full_fig_p016_14.png] view at source ↗
Figure 15
Figure 15. Figure 15: Comparison of the initial configuration (black line) and the deformed slope after 6 s: (a) [PITH_FULL_IMAGE:figures/full_fig_p016_15.png] view at source ↗
read the original abstract

Conventional smoothed particle hydrodynamics based on Eulerian kernels (CESPH) is widely-used in large deformation analysis in geomaterials. Despite being popular, it suffers from tensile instability and rank-deficiency; thus, it needs several numerical treatments to be stable. In this work, we present a stabilized total-Lagrangian SPH method (TLSPH), which is inherently free of tensile instability. A stiffness-based hourglass control algorithm is employed to cure the hourglass mode caused by rank-deficiency. Periodic update of reference configuration is used in simulations to allow TLSPH to model large deformation and post-failure flow in geomaterials. Several numerical examples are presented to show the performance of the stabilized TLSPH method. The comparison between TLSPH and CESPH are discussed. The influences of hourglass control and configuration update are also discussed and shown in the numerical examples. It is found that the presented stabilized TLSPH is robust and can model large deformation and plastic flows in geomaterials. Particularly, the stabilized TLSPH delivers accurate and smooth stress results.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit. Tearing a paper down is the easy half of reading it; the pith above is the substance, this is the friction.

Referee Report

3 major / 2 minor

Summary. The manuscript proposes a stabilized total-Lagrangian SPH (TLSPH) formulation for large-deformation and failure problems in geomaterials. It asserts that the total-Lagrangian kernel choice inherently removes tensile instability, that a stiffness-based hourglass control cures rank-deficiency modes, and that periodic reference-configuration updates enable continued simulation of post-failure plastic flow while retaining the stability advantages of the TL approach. Several numerical examples are presented to demonstrate robustness, accuracy, and smoothness of the resulting stress fields relative to conventional Eulerian SPH.

Significance. If the quantitative accuracy and stability claims are substantiated, the work would supply a useful alternative to the multiple ad-hoc corrections typically required for Eulerian SPH in geomechanics. The direct change of reference configuration rather than additional stabilization layers is a constructive feature.

major comments (3)
  1. [Numerical Examples] Numerical Examples section: the central claim that the stabilized TLSPH “delivers accurate and smooth stress results” is unsupported by any L2-norm, convergence rate, or other quantitative error measure against reference solutions, especially after reference-configuration updates; visual inspection alone cannot establish the accuracy asserted in the abstract and conclusion.
  2. [Method (hourglass control)] Hourglass-control subsection: the stiffness-based algorithm is introduced without a supporting stability analysis or bound showing that the added term eliminates rank-deficiency modes without degrading stress accuracy or introducing new artifacts; this is load-bearing for the robustness claim.
  3. [Method (configuration update)] Reference-configuration-update subsection: no analysis or test is given demonstrating that periodic updates preserve the tensile-instability-free property of the total-Lagrangian formulation once deformation becomes extreme; the weakest assumption identified in the stress-test note therefore remains unaddressed.
minor comments (2)
  1. [Abstract] Abstract: the statement of findings should reference the specific figures or tables that contain the TLSPH–CESPH comparisons.
  2. [Method] Notation: the definition and units of the hourglass stiffness coefficient should be stated explicitly when first introduced.

Simulated Author's Rebuttal

3 responses · 0 unresolved

We thank the referee for the constructive and detailed comments. The suggestions highlight areas where the manuscript can be strengthened with additional quantitative support and discussion. We address each major comment below and indicate the revisions we will make.

read point-by-point responses

  1. Referee: [Numerical Examples] Numerical Examples section: the central claim that the stabilized TLSPH “delivers accurate and smooth stress results” is unsupported by any L2-norm, convergence rate, or other quantitative error measure against reference solutions, especially after reference-configuration updates; visual inspection alone cannot establish the accuracy asserted in the abstract and conclusion.

    Authors: We agree that quantitative error measures would provide stronger substantiation. In the revised manuscript we will add L2-norm errors for displacement and von Mises stress against analytical or FEM reference solutions in the elastic and small-deformation benchmarks. For the large-deformation and post-failure examples we will report quantitative comparisons of final geometry and average stress values against results obtained with established Eulerian SPH and MPM implementations. These additions will directly address accuracy after reference-configuration updates. revision: yes

  2. Referee: [Method (hourglass control)] Hourglass-control subsection: the stiffness-based algorithm is introduced without a supporting stability analysis or bound showing that the added term eliminates rank-deficiency modes without degrading stress accuracy or introducing new artifacts; this is load-bearing for the robustness claim.

    Authors: The stiffness-based hourglass term is constructed to penalize only the zero-energy modes while leaving the consistent part of the stress field unchanged; its form is taken from established finite-element hourglass control and adapted to the total-Lagrangian SPH discretization. Effectiveness is demonstrated by the numerical examples in which the control suppresses oscillations that appear when the term is omitted. A full von Neumann analysis for the nonlinear elasto-plastic case lies outside the present scope, but we will expand the subsection with a short derivation showing that the added stiffness acts only on the rank-deficient kernel and does not alter the consistent stress evaluation. revision: partial

  3. Referee: [Method (configuration update)] Reference-configuration-update subsection: no analysis or test is given demonstrating that periodic updates preserve the tensile-instability-free property of the total-Lagrangian formulation once deformation becomes extreme; the weakest assumption identified in the stress-test note therefore remains unaddressed.

    Authors: Because each reference-configuration update creates a new total-Lagrangian segment in which the kernel is again defined with respect to the current reference, the tensile-instability-free property is retained locally. We will add a dedicated numerical test that tracks particle clustering and negative-pressure indicators across multiple updates in a large simple-shear problem. We will also note in the text that a general mathematical proof for arbitrary update intervals remains an open question and is left for future work. revision: yes

Circularity Check

0 steps flagged

No circularity: method is direct formulation with external validation via examples

full rationale

The paper introduces a total-Lagrangian SPH formulation that is stated to be inherently free of tensile instability, augmented by a stiffness-based hourglass control and periodic reference updates. These are presented as explicit algorithmic choices to address known CESPH issues, followed by numerical examples for validation. No equations, predictions, or central claims reduce by construction to fitted parameters, self-definitions, or self-citation chains. The derivation chain consists of standard SPH kernel changes plus standard stabilization techniques, with performance assessed against independent benchmarks rather than internal tautologies. This is self-contained against external numerical tests.

Axiom & Free-Parameter Ledger

1 free parameters · 1 axioms · 0 invented entities

Abstract-only access prevents identification of most implementation details; the hourglass stiffness coefficient is expected to be a tunable parameter, and the claim that total-Lagrangian kernels are inherently tensile-instability-free is treated as a domain assumption.

free parameters (1)
  • hourglass stiffness coefficient
    Parameter in the stiffness-based hourglass control algorithm used to suppress rank-deficiency modes; its specific value or fitting procedure is not stated in the abstract.
axioms (1)
  • domain assumption Total-Lagrangian SPH formulation is inherently free of tensile instability
    Directly stated in the abstract as a property of the TLSPH approach.

pith-pipeline@v0.9.0 · 5712 in / 1210 out tokens · 27907 ms · 2026-05-24T21:32:26.135700+00:00 · methodology

discussion (0)

Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.

Reference graph

Works this paper leans on

33 extracted references · 33 canonical work pages

  1. [1]

    J. J. Monaghan, Smoothed particle hydrodynamics and its diverse applications, Annual Review of Fluid Mechanics 44 (2012) 323–346

    work page 2012

  2. [2]

    J. S. Chen, M. Hillman, S. W. Chi, Meshfree methods: progress made after 20 years, Journal of Engineering Mechanics 143 (4) (2017) 04017001

    work page 2017

  3. [3]

    S. G. Bardenhagen, E. M. Kober, The generalized interpolation material point method, Computer Modeling in Engineering and Sciences 5 (6) (2004) 477–496

    work page 2004

  4. [4]

    Maeda, H

    K. Maeda, H. Sakai, M. Sakai, Development of seepage failure analysis method of ground with smoothed particle hydrodynamics, JSCE Structural Engineering/Earthquake Engineering 23 (2) (2006) 307s–319s

    work page 2006

  5. [5]

    Naili, T

    M. Naili, T. Matsushima, Y . Yamada, A 2D smoothed particle hydrodynamics method for liquefaction induced lateral spreading analysis, JSCE Journal of Applied Mechanics 8 (2005) 591–599

    work page 2005

  6. [6]

    H. H. Bui, R. Fukagawa, K. Sako, S. Ohno, Lagrangian meshfree particles method (SPH) for large deformation and failure flows of geo- material using elastic–plastic soil constitutive model, International Journal for Numerical and Analytical Methods in Geomechanics 32 (12) (2008) 1537–1570

    work page 2008

  7. [7]

    H. H. Bui, R. Fukagawa, K. Sako, J. C. Wells, Slope stability analysis and discontinuous slope failure simulation by elasto-plastic smoothed particle hydrodynamics (SPH), G´eotechnique 61 (7) (2011) 565–574

    work page 2011

  8. [8]

    H. H. Bui, R. Fukagawa, An improved SPH method for saturated soils and its application to investigate the mechanisms of embankment failure: Case of hydrostatic pore-water pressure, International Journal for Numerical and Analytical Methods in Geomechanics 37 (1) (2013) 31–50

    work page 2013

  9. [9]

    H. H. Bui, J. K. Kodikara, A. Bouazza, A. Haque, P. G. Ranjith, A novel computational approach for large deformation and post-failure analyses of segmental retaining wall systems, International Journal for Numerical and Analytical Methods in Geomechanics 38 (13) (2014) 1321–1340

    work page 2014

  10. [10]

    Chambon, R

    G. Chambon, R. Bouvarel, D. Laigle, M. Naaim, Numerical simulations of granular free-surface flows using smoothed particle hydrodynam- ics, Journal of Non-Newtonian Fluid Mechanics 166 (12-13) (2011) 698–712

    work page 2011

  11. [11]

    C. Peng, X. G. Guo, W. Wu, Y . Q. Wang, Unified modelling of granular media with smoothed particle hydrodynamics, Acta Geotechnica 11 (6) (2016) 1231–1247

    work page 2016

  12. [12]

    J. Wang, D. Chan, Frictional contact algorithms in SPH for the simulation of soil–structure interaction, International Journal for Numerical and Analytical Methods in Geomechanics 38 (7) (2014) 747–770

    work page 2014

  13. [13]

    L. Zhan, C. Peng, B. Y . Zhang, W. Wu, Three-dimensional modeling of granular flow impact on rigid and deformable structures, Computers and Geotechnics 112 (2019) 257–271

    work page 2019

  14. [14]

    G. C. Ganzenm ¨uller, M. Sauer, M. May, S. Hiermaier, Hourglass control for smooth particle hydrodynamics removes tensile and rank- deficiency instabilities, The European Physical Journal Special Topics 225 (2) (2016) 385–395

    work page 2016

  15. [15]

    Swegle, D

    J. Swegle, D. Hicks, S. Attaway, Smoothed particle hydrodynamics stability analysis, Journal of Computational Physics 116 (1) (1995) 123–134

    work page 1995

  16. [16]

    C. T. Dyka, P. W. Randles, R. P. Ingel, Stress points for tension instability in SPH, International Journal for Numerical Methods in Engineering 40 (13) (1997) 2325–2341. 17

    work page 1997

  17. [17]

    J. J. Monaghan, SPH without a tensile instability, Journal of Computational Physics 159 (2) (2000) 290–311

    work page 2000

  18. [18]

    J. P. Gray, J. J. Monaghan, R. P. Swift, SPH elastic dynamics, Computer Methods in Applied Mechanics and Engineering 190 (49-50) (2001) 6641–6662

    work page 2001

  19. [19]

    C. Peng, S. Wang, W. Wu, H. S. Yu, C. Wang, J. Y . Chen, LOQUAT: an open-source GPU-accelerated SPH solver for geotechnical modeling, Acta Geotechnica (2019) 1–19

    work page 2019

  20. [20]

    Belytschko, Y

    T. Belytschko, Y . Guo, W. Kam Liu, S. Ping Xiao, A unified stability analysis of meshless particle methods, International Journal for Numerical Methods in Engineering 48 (9) (2000) 1359–1400

    work page 2000

  21. [21]

    Bonet, S

    J. Bonet, S. Kulasegaram, Alternative total Lagrangian formulations for corrected smooth particle hydrodynamics (CSPH) methods in large strain dynamic problems, Revue Europ´eenne des ´El´ements Finis 11 (7-8) (2002) 893–912

    work page 2002

  22. [22]

    Vignjevic, J

    R. Vignjevic, J. R. Reveles, J. Campbell, SPH in a total Lagrangian formalism, Computer Modeling in Engineering & Sciences 14 (3) (2006) 181

    work page 2006

  23. [23]

    G. C. Ganzenm ¨uller, An hourglass control algorithm for Lagrangian smooth particle hydrodynamics, Computer Methods in Applied Mechan- ics and Engineering 286 (2015) 87–106

    work page 2015

  24. [24]

    L. Zhan, C. Peng, B. Zhang, W. Wu, A stabilized TL–WC SPH approach with GPU acceleration for three–dimensional fluid–structure interaction, Journal of Fluids and Structures 86 (2019) 329–353

    work page 2019

  25. [25]

    Wendland, Piecewise polynomial, positive definite and compactly supported radial functions of minimal degree, Advances in computa- tional Mathematics 4 (1) (1995) 389–396

    H. Wendland, Piecewise polynomial, positive definite and compactly supported radial functions of minimal degree, Advances in computa- tional Mathematics 4 (1) (1995) 389–396

    work page 1995

  26. [26]

    Monaghan, R

    J. Monaghan, R. A. Gingold, Shock simulation by the particle method sph, Journal of Computational Physics 52 (2) (1983) 374–389

    work page 1983

  27. [27]

    C. Peng, W. Wu, H. S. Yu, C. Wang, A SPH approach for large deformation analysis with hypoplastic constitutive model, Acta Geotechnica 10 (6) (2015) 703–717

    work page 2015

  28. [28]

    Belytschko, S

    T. Belytschko, S. Xiao, Stability analysis of particle methods with corrected derivatives, Computers & Mathematics with Applications 43 (3-5) (2002) 329–350

    work page 2002

  29. [29]

    Rabczuk, T

    T. Rabczuk, T. Belytschko, S. Xiao, Stable particle methods based on Lagrangian kernels, Computer Methods in Applied Mechanics and Engineering 193 (12-14) (2004) 1035–1063

    work page 2004

  30. [30]

    De Vuyst, R

    T. De Vuyst, R. Vignjevic, Total lagrangian sph modelling of necking and fracture in electromagnetically driven rings, International Journal of Fracture 180 (1) (2013) 53–70

    work page 2013

  31. [31]

    Leroch, M

    S. Leroch, M. Varga, S. Eder, A. Vernes, M. R. Ripoll, G. Ganzenm¨uller, Smooth particle hydrodynamics simulation of damage induced by a spherical indenter scratching a viscoplastic material, International Journal of Solids and Structures 81 (2016) 188–202

    work page 2016

  32. [32]

    C. Dyka, R. Ingel, An approach for tension instability in smoothed particle hydrodynamics (sph), Computers & structures 57 (4) (1995) 573–580

    work page 1995

  33. [33]

    D. V . Griffiths, P. A. Lane, Slope stability analysis by finite elements, Geotechnique 49 (3) (1999) 387–403. 18

    work page 1999