Pith. sign in

REVIEW 4 major objections 5 minor 21 references

Elastoplasticity with softening as a state-dependent sweeping process: non-uniqueness of solutions and emergence of shear bands in lattices of springs

T0 review · 4 major / 5 minor · reviewed 2026-08-05 · deepseek-v4-flash

Pith's one-line read The full Lattice Spring Model with softening, hardening, and perfectly plastic springs is equivalent to a state-dependent sweeping process.

desk verdict A genuinely useful sweeping-process formulation for lattice elastoplasticity with softening, but the paper's main numerical claim rests on an unproven existence step that the author explicitly defers. read the letter →

arxiv 2508.16333 v2 pith:MKT4YZQ7 submitted 2025-08-22 math.OC cond-mat.softmath-phmath.MP

classification math.OCcond-mat.softmath-phmath.MP MSC 74C0574N3047J2047J26
keywords state-dependentsweepingprocesselastoplasticitywithsofteningnon-uniquenessofsolutionsshearbandstrainlocalizationlatticespringmodelimplicitcatch-upalgorithmplasticitybifurcation
verification ladder T0 review T1 audit T2 compute T3 formal

The pith

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

The reading

This paper tries to show that the quasi-static evolution of a lattice of elasto-plastic springs—including the mathematically troublesome softening case—can be written exactly as a state-dependent sweeping process, a differential inclusion with a moving convex set. If the equivalence is right, the same numerical and analytical machinery applies to hardening, perfect plasticity, and softening, with the difference between them reduced to a parameter. The paper demonstrates the payoff by analytically exhibiting coexisting solutions in a two-spring softening system and by simulating shear-band formation with strain localization in regular lattices. This matters because softening plasticity is normally ill-posed: stress-strain relations lose monotonicity, models become mesh-dependent, and solutions can branch. A single finite-dimensional formulation that still produces well-defined evolutions up to spring failure would give a rigorous route into ductile-fracture initiation.

What carries the argument

The load-bearing object is the state-dependent moving set C(y, q): for each spring, a convex polytope of admissible stresses whose upper and lower thresholds shift with the damage variable q and with the applied displacement load. The sweeping-process inclusion dot-y in -N_C(y,q)(y), together with the linear compatibility and equilibrium constraints, encodes the whole lattice. The numerical scheme is the implicit catch-up map, whose fixed points at each time step are the discrete solutions; stability of these fixed points under iteration is what the paper uses to explain bifurcation and branching in softening.

What would settle it

Take the two-spring softening system of Section 5, set both springs to softening parameters with no complete failure, and choose the initial condition on the distributed saddle solution satisfying (106)-(107). Run the implicit catch-up iterations with a small time step; if for any step before failure the projection map has no fixed point or the iteration diverges or cycles, then the claimed sweeping-process formulation does not produce the evolution asserted for softening.

Watch

Extended reading notes

Core claim

At the center is a finite-dimensional equivalence (Theorem 4.1): the full displacement-stress-damage system (LSM1)-(LSM6) for a lattice built of springs of any mixture—hardening, perfectly plastic, softening—holds if and only if the pair (y, q) of reduced stress and damage variables solves the state-dependent sweeping process (58)-(59), with the moving set C(y, q) the polytope of admissible stresses. The equivalence is two-sided, which is stronger than earlier sweeping-process formulations for perfectly plastic lattices. For the toy lattice of two springs in series, the paper solves the sweeping process explicitly and finds coexisting solutions when both springs soften: a distributed solutio

Load-bearing premise

The implied existence of a well-defined solution for softening lattices: at every time step the implicit catch-up map must have a fixed point that the iteration can actually reach, even though the standard existence condition behind the cited theorem fails exactly for softening springs.

Editorial extensions

If this is right

  • Changes in the plastic regime become parameter changes: beta crossing 1 turns a hardening spring into a softening one, so one solver and one theory cover all three cases.
  • Softening lattices admit multiple coexisting evolutions; the catch-up iteration's fixed points classify them as stable or saddle, giving a route to predict which branch a load path selects.
  • Shear bands with strain localization appear as emergent, symmetry-breaking events in symmetric lattices, and their orientation is controlled by lattice geometry and defects, so the model can reproduce ductile-fracture precursors.
  • Evolution can be computed reliably only up to the first spring's complete failure; past that point the sweeping process has no admissible continuation, matching the idea of material separation.
  • The two-spring toy model provides a minimal analytically solvable case where the nonsmooth bifurcation underlying softening can be studied rigorously.

Reading between the lines

Editorial extensions of the paper, not claims the author makes directly.

  • If the promised local existence theorem for softening is established, the discrete model could serve as a regularized testbed in which shear-band width is controlled by lattice spacing, allowing mesh-dependence questions to be studied by refinement.
  • The observed relation between fixed-point stability and energy minima suggests a conjecture the paper does not prove: each iteration of the catch-up map decreases an underlying energy toward a local minimum; a cycling example would refute it, and a proof would unify variational and sweeping-process approaches.
  • The degeneracy in the plasticity-modulus formula means many parameter pairs produce the same observable modulus, so calibration from macroscopic stress-strain curves alone cannot identify them; extra physical or atomistic information would be needed to predict damage evolution and failure location.
  • A measure-valued state-dependent sweeping process for continuous media, which the paper sketches as conceivable, would complement the existing continuum variational theory for softening and could yield two-sided error estimates and optimal-control formulations for ductile-fracture initiation.
Share X Bluesky LinkedIn Reddit HN

Signed reviews

No signed human review yet.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, and a circularity audit.

Referee Report

4 major / 5 minor

Summary. The paper develops a state-dependent sweeping process formulation for quasi-static elasto-plastic lattice spring models whose springs may harden, soften, or be perfectly plastic. After analyzing a single spring (§3), it formulates the full lattice equations (LSM1)–(LSM6) and proves in Theorem 4.1 an equivalence between these equations and the state-dependent sweeping process (58)–(59), including reconstruction formulas (60)–(65). It then introduces implicit catch-up algorithms (§4.3–4.4), analyzes a two-spring toy model (§5) where non-uniqueness is asserted in Proposition 5.1, and presents numerical simulations of rectangular and triangular lattices (§6) showing shear-band localization in softening springs. The author is explicit that for softening and perfect plasticity the standard existence theorem (Theorem 1.2) does not apply because condition (93) fails, and that a local existence theorem is future work (§7).

Significance. If the central equivalence and the existence gap were closed, the paper would provide a unified finite-dimensional framework for rate-independent plasticity with softening, connecting the mechanical problem to a state-dependent sweeping process and enabling a rigorous study of non-uniqueness and shear-band bifurcations. The two-sided equivalence proof of Theorem 4.1 is detailed and self-contained, the toy-model analysis addresses a genuine non-uniqueness phenomenon, and the numerical experiments are accompanied by available code and videos. The derivation does not appear to presuppose its conclusions; the main weakness is not circularity but a missing existence theorem in the softening regime. As it stands, the advertised claim that the sweeping process 'solves' the evolution with softening outruns the proved results.

major comments (4)
  1. [§4.2, Remark 4.5, §7] The paper's central claim is that the state-dependent sweeping process solves the evolution for arbitrary mixtures including softening, but no existence theorem is proved in this regime. Theorem 1.2 requires the Lipschitz condition (93), and Remark 4.5 explicitly notes that (93) cannot hold when β_i ≤ 1 (softening/perfectly plastic springs). Section 7 defers a local existence result to future work. Since Theorem 4.1 only establishes equivalence for a given solution of (58)–(59), and Algorithm 1/2 is defined by fixed points of (91)/(99), the Section 6 simulations are currently not known to be solutions of the time-continuous problem. Please either prove existence of a fixed point for each time step under a suitable non-degeneracy assumption, or explicitly restrict the claims to the discrete fixed-point problem and state the absence of existence as a limitation in the abstract and conclusi
  2. [§5.3, Proposition 5.1] The analytical non-uniqueness result is only sketched. The proof says 'To save space we will omit complete derivations' and reports that equations were solved manually. Formulas (108)–(109) and the case analysis are central to the claim of coexisting solutions and to the bifurcation picture in Figure 8. Please provide a complete derivation (or put it in a supplementary file) with the inequality checks (110) and all case distinctions, rather than a sketch.
  3. [§6, Algorithm 2] The numerical experiments do not report the time-step size, the tolerance used in Algorithm 2, the maximum number of iterations, or the initial-guess policy; Remark 4.6 states that iteration may diverge when there are no stable fixed points. For the shear-band runs in Figs. 10–16, please report for every run that the convergence criterion was met at each time step, and provide data underlying the iteration-count plots. This is necessary to distinguish genuine solutions from algorithm failure and to support the claim that modifying the initial guess yields a second solution (Fig. 3 vs Fig. 10d).
  4. [§4.2, Theorem 4.1 proof] In the sweep-to-LSM direction the proof says 'we do not consider the state of complete failure' and uses implication (41). The theorem statement, however, does not state this exclusion. If the equivalence is intended to cover complete-failure states, the proof is incomplete; if not, the statement and the abstract's 'arbitrary placement' should be qualified. This is not merely cosmetic, since complete failure is reached in the softening simulations (Remark 3.1).
minor comments (5)
  1. [§1, §4.1.2, §6] Several typos and language issues: 'Given finctions' should be 'Given functions'; 'not not alter' in §4; 'modyfing' and 'incuded' in §6. A careful proofreading pass is needed.
  2. [§4.2, equations (67)–(74)] The notation for the moving set is very dense. A short table collecting the symbols in (55), (60), (70), (72), (73)–(74) would improve readability.
  3. [§5, Figures 8–9] The captions do not define the coordinate axes or the meaning of the wedges and red/black points. Please add explicit axis labels and a legend.
  4. [§6] The numerical setup is incomplete: no time-step size, tolerance, or iteration limits are given. These should be stated for reproducibility, especially because Algorithm 2's convergence is not guaranteed in the softening case.
  5. [§4.1] The paper relies heavily on [Gud+23] for definitions and algebraic properties. Please state explicitly in Section 4 which results are taken from that reference and which are new, so that the equivalence theorem can be checked without reading the full prior paper.

Circularity Check

0 steps flagged · score 0.0 of 10

No circular reduction: the sweeping-process formulation is a proven equivalence, and the shear-band/non-uniqueness outputs are not fitted inputs; the acknowledged missing existence theorem for softening is a completeness gap, not circularity.

full rationale

Theorem 4.1 gives a two-sided derivation: equations (LSM1)-(LSM6) are transformed into (58)-(59) by direct substitution of Hooke's law, compatibility, equilibrium and the normal-cone form of the flow rule, and the reverse direction reconstructs (LSM4)-(LSM5) and (LSM3) from the sweeping-process inclusion. Neither direction postulates the target behavior; the equivalence is algebraic and is proved in the text, with the lattice-construction material from [Gud+23] used as algorithmic/geometric background rather than as the source of the softening results. The plastic parameters (beta_i, gamma_i, etc.) are assigned before the simulations, and the non-symmetric shear bands (Fig. 10), the bifurcation of fixed points (Fig. 8, Observation 5.1), and the non-uniqueness in Proposition 5.1 are outputs of the model, not parameters fitted to those outputs. The self-citations [Gud+23, GM21] provide the lattice-rigidity and pseudoinverse identities and an external uniqueness theorem for perfect plasticity; they are not invoked to forbid alternatives or to define the softening sweeping process, so they do not constitute circularity. The paper explicitly flags that the standard existence condition (93) fails when beta_i <= 1 (Remark 4.5), that Algorithm 1 can diverge when no stable fixed point exists (Remark 4.6), and that a local existence theorem for softening is future work (Section 7); Proposition 5.1 is also only sketched. These are genuine correctness/completeness limitations of the softening claim, but they are not self-referential reductions of the derivation to its own inputs. The absence of circularity is therefore consistent with a nonzero-risk assessment on other grounds.

Assumptions & free parameters 2 free parameters · 5 assumptions · 0 invented entities

No new particles, forces, or dimensions are introduced; the 'damage variable' z is a standard internal variable in plasticity. The load-bearing postulates are the linear polytope constitutive ansatz, the small-strain rigidity assumptions, and, for softening, the unproven existence of sweeping process solutions.

free parameters (2)
  • spring constitutive parameters (E, a, b, yield stress) in numerical examples = e.g., E=2, a=0.4, b=-0.3, sigma_y0=0.4 for softening (113)
    Chosen by hand in Section 6 to place each spring in the softening, hardening or perfectly plastic regime. They are not fitted to data, but they select the regime that generates the observed shear bands.
  • lattice discretization and loading steps = e.g., 12x18 nodes, time step not specified
    The numerical demonstrations use specific lattice sizes and displacement load steps; no convergence study with respect to these is reported, leaving sensitivity of the shear band patterns unquantified.
assumptions (5)
  • standard math Convex analysis: normal cone calculus, projections, Moore-Penrose pseudoinverse properties (Lemma 2.1, Prop 4.1)
    Used throughout the derivation of the sweeping process, e.g., equations (14)-(17), (52)-(54).
  • domain assumption Constitutive ansatz: material state is (stress, internal variable) with the specific linear state-dependent polytope (22), (37)
    The paper restricts to linear isotropic hardening/softening with the particular set (22); the central claim holds only within this constitutive class.
  • domain assumption Quasi-static equilibrium and small-strain linearized geometric compatibility (LSM4)-(LSM6)
    The model is rate-independent and uses fixed reference geometry (43)-(44); large deformations are excluded.
  • ad hoc to paper Assumptions 1-4: kinematic determinacy (45)-(47), nonzero self-stress space (53), and regularity of loads (Assumption 3)
    These guarantee the algebraic reductions in Theorem 4.1; they are stated as assumptions rather than derived.
  • ad hoc to paper Existence of solutions to the state-dependent sweeping process in the softening case
    Theorem 1.2's condition (93) fails for softening (Remark 4.5). The paper assumes the implicit catch-up map has fixed points for the numerical simulations; a local existence theorem is deferred to future work (Section 7).

how reviews work

0 comments
Cite this review

Pith. "Pith review of Elastoplasticity with softening as a state-dependent sweeping process: non-uniqueness of solutions and emergence of shear bands in lattices of springs." pith.science (2026). https://pith.science/paper/MKT4YZQ7

@misc{pith2026250816333,
  author       = {Pith},
  title        = {Pith review of: Elastoplasticity with softening as a state-dependent sweeping process: non-uniqueness of solutions and emergence of shear bands in lattices of springs},
  year         = {2026},
  howpublished = {\url{https://pith.science/paper/MKT4YZQ7}},
  note         = {Machine review of arXiv:2508.16333}
}
read the original abstract

Plasticity with softening and fracture mechanics lead to ill-posed mathematical problems due to the loss of monotonicity. Multiple co-existing solutions are possible when softening elements are coupled together, and solutions cannot be continued beyond the point of complete degradation of the set of admissible stresses. We present a state-dependent sweeping process which solves the evolution of elasto-plastic Lattice Spring Models with arbitrary placement of softening, hardening and perfectly plastic springs. Using numerical simulations of regular grid lattices with softening we demonstrate the emergence of non-symmetric shear bands with strain localization. At the same time, in toy examples it is easy to analytically derive multiple co-existing solutions. These solutions correspond to fixed points in the implicit catch-up algorithm and we observe a discontinuous bifurcation with the exchange of stability of those fixed points.

Discussion (0). Continue with ORCID to comment.

Reference graph

Works this paper leans on

21 extracted references · 16 canonical work pages

  1. [14]

    Elsevier, 2008, pp

    Dislocations in Solids. Elsevier, 2008, pp. 1–42. ���: https://doi.org/10.1016/S1572-4859(07)00001-0. 61

  2. [15]

    Comparison of various models for strain-softening

    [PBT88] G. Pijaudier-Cabot, Z.P. Bažant, and M. Tabbara. “Comparison of various models for strain-softening”. In: 5.2 (Feb. 1988), pp. 141–150.���: 10.1108/eb023732. [Rd89] J.G. Rots and R. de Borst. “Analysis of concrete fracture in “direct” tension”. In:���� �� ������ �������25.12 (1989), pp. 1381–1394.���: https://doi.org/10.1016/0020- 7683(89)90107-8....

  3. [27]

    Chapter 81 - Atomistic Simulations of Dislocations in FCC Metallic Nanocrystalline Materials

    Advances in Applied Mechanics. Elsevier, 1989, pp. 83–151.���: https://doi.org/10.1016/S0065-2156(08)70195-9. [VD08] H. Van Swygenhoven and P. M. Derlet. “Chapter 81 - Atomistic Simulations of Dislocations in FCC Metallic Nanocrystalline Materials”. In:� ������� �� �� �� �� �������. Ed. by J. P. Hirth. Vol

  4. [1970]

    Structural softening, mesh dependence, and regular- isation in non-associated plastic flow

    [SB19] S.A. Sabet and R. de Borst. “Structural softening, mesh dependence, and regular- isation in non-associated plastic flow”. In:���� �� ������ ����� ������� �������� 43.13 (2019), pp. 2170–2183.���: https://doi.org/10.1002/nag.2973. [Sch06] H. C. Schnabel. “Zur Wohlgestelltheit des Gurson-Modells”. PhD thesis. Technische Universität München, 2006, p

  5. [1976]

    Some mathematical problems in perfect plasticity

    [ER04] F. Ebobisse and B. D. Reddy. “Some mathematical problems in perfect plasticity”. In:������� ������� ����� ����� ������193.48-51 (2004), pp. 5071–5094.���: 10.1016/j.cma.2004.07.002. [ET99] I. Ekeland and R. Temam.������ �������� ��� ����������� ��������. Classics in Applied Mathematics. Society for Industrial and Applied Mathematics,

  6. [1981]

    An analysis of ductile rupture in notched bars

    [NT84] A. Needleman and V. Tvergaard. “An analysis of ductile rupture in notched bars”. In: �� ����� ����� ������32.6 (1984), pp. 461–490.���: https://doi.org/10.1016/0022- 5096(84)90031-0. [OKZ13] J. Outrata, M. Kočvara, and J. Zowe.��������� �������� �� ������������ ����� ���� ���� ����������� ������������ ������� ������������ ��� ��������� �������. Non...

  7. [1984]

    Variational formulation of softening phe- nomena in fracture mechanics. The one-dimensional case

    [BDG99] A. Braides, G. Dal Maso, and A. Garroni. “Variational formulation of softening phe- nomena in fracture mechanics. The one-dimensional case”. In:����� ������� ����� �����146.1 (1999), pp. 23–58.���: 10.1007/s002050050135. [BKS04] M. Brokate, P. Krejčí, and H. Schnabel. “On uniqueness in evolution quasivariational inequalities”. In:�� ������� �����1...

  8. [1985]

    Functions of bounded deformation

    [TS81] R. Temam and G. Strang. “Functions of bounded deformation”. In:����� �������� ����� �����75.1 (1980/81), pp. 7–21.���: 10.1007/BF00284617. [Tve89] V. Tvergaard. “Material Failure by Void Growth to Coalescence”. In: ed. by J W. Hutchinson and T. Y. Wu. Vol

Show all 21 references
  1. [1991]

    A generalized 2D non-local lattice spring model for fracture simulation

    [Che+14] H. Chen, E. Lin, Y. Jiao, and Y. Liu. “A generalized 2D non-local lattice spring model for fracture simulation”. In:������� �����54.6 (Dec. 2014), pp. 1541–1558. ���: 10.1007/s00466-014-1075-4. [CIY09] C. Castaing, A. G. Ibrahim, and M. Yarou. “Some Contributions to N...

  2. [1996]

    Elastic–ideally plastic beams and Prandtl–Ishlinskii hys- teresis operators

    [KS07] P. Krejčí and J. Sprekels. “Elastic–ideally plastic beams and Prandtl–Ishlinskii hys- teresis operators”. In:������������ ������� �� ��� ������� ��������30.18 (2007), pp. 2371–2393.���: https://doi.org/10.1002/mma.892. [Lem17] A. Lemaître. “Inherent stress correlations ...

  3. [1998]

    Mathematical and Computational Methods in Plasticity

    [SMT80] G. Strang, H. Matthies, and R. Temam. “Mathematical and Computational Methods in Plasticity”. In:����������� ������� �� ��� ��������� �� ������. Ed. by S. Nemat- Nasser. Pergamon, 1980, pp. 20–28.���: https://doi.org/10.1016/B978-0-08- 024728-1.50008-2. [Str24] S. Stro...

  4. [1999]

    Deflation for semismooth equations

    [FCS20] P.E. Farrell, M. Croci, and T.M. Surowiec. “Deflation for semismooth equations”. In:������ ������� ������35.6 (2020), pp. 1248–1271.���: 10.1080/10556788.2019. 1613655. [FGM16] G.A. Francfort, A. Giacomini, and J.-J. Marigo. “The elasto-plastic exquisite corpse: a Suqu...

  5. [2002]

    Strain localization and ductile fracture in advanced high-strength steel sheets

    58 [Gru+17] G. Gruben, D. Morin, M. Langseth, and O.S. Hopperstad. “Strain localization and ductile fracture in advanced high-strength steel sheets”. In:�������� ������� �� ��������� � ��������61 (2017), pp. 315–329.���: https://doi.org/10.1016/j. euromechsol.2016.09.014. [Gud...

  6. [2008]

    Energy profile and bifurcation analysis in softening plastic- ity

    [CB04] G. Chen and G. Baker. “Energy profile and bifurcation analysis in softening plastic- ity”.In:���� ������� ����7.6(Dec.2004),pp.515–523.���:10.1260/1369433042863224. [CH91] J. Clark and D.A. Holton.� ����� ���� �� ����� ������. World Scientific,

  7. [2009]

    Snap-through and snap-back response in concrete structures and the dangers of under-integration

    [Cri86] M.A. Crisfield. “Snap-through and snap-back response in concrete structures and the dangers of under-integration”. In:���� �� ������ ������� ����22.3 (1986), pp. 751– 767.���: https://doi.org/10.1002/nme.1620220314. [Dal+08a] G. Dal Maso, A. DeSimone, M. G. Mora, and M...

  8. [2010]

    Mechanics of Distributed Cracking

    [Baž86] Z.P. Bažant. “Mechanics of Distributed Cracking”. In:����� ����� ����39.5 (May 1986), pp. 675–705.���: 10.1115/1.3143724. [BC11] H.H. Bauschke and P.L. Combettes.������ �������� ��� �������� �������� ���� ��� �� ������� ������. CMS Books in Mathematics. Springer New York,

  9. [2012]

    Mathematical analysis of strain localization

    [Jir07] M. Jirásek. “Mathematical analysis of strain localization”. In:���� ���� ����� ���� 11.7-8 (2007), pp. 977–991.���: 10.1080/17747120.2007.9692973. [JR03] M. Jirásek and S. Rolshoven. “Regularized formulations of strain-softening plas- ticity”. In:�������� ������������ ...

  10. [2014]

    Stabilization of the response of cyclically loaded lattice spring models with plasticity

    [GM21] I. Gudoshnikov and O. Makarenkov. “Stabilization of the response of cyclically loaded lattice spring models with plasticity”. In:����� ������� ������ ����� ���� 27 (2021), Paper No. S8, 43.���: 10.1051/cocv/2020043. [GPS02] H. Goldstein, C. Poole, and J. Safko.���������...

  11. [2015]

    Damage-augmentednonlocallatticeparticlemethodforfracture simulation of solids

    [ML22] CMengandYLiu.“Damage-augmentednonlocallatticeparticlemethodforfracture simulation of solids”. In:���� �� ������ �������243 (2022), p. 111561.���: https: //doi.org/10.1016/j.ijsolstr.2022.111561. [MMD03] Lambrecht M., C. Miehe, and J. Dettmar. “Energy relaxation of non-c...

  12. [2018]

    The non-uniqueness of the atomistic stress tensor and its relationship to the generalized Beltrami representation

    [Ant+25] H. Antil, R. Arndt, B.S. Mordukhovich, D. Nguyen, and C.N. Rautenberg.������� ������� �� � ����������������� �������� �������. 2025.���: 10.3934/mcrf.2025010. [AT16] N. C. Admal and E. B. Tadmor. “The non-uniqueness of the atomistic stress tensor and its relationship ...

  13. [2024]

    Sur un nouveau cadre fonctionnel pour les équations de la plasticité

    [Sup] SupplementalMaterial.,seeathttps://github.com/Ivan-Gudoshnikov/Elastoplasticity- softening-public for videos and computer programs of the simulations. [Suq78] P.-M. Suquet. “Sur un nouveau cadre fonctionnel pour les équations de la plasticité”. In:�� �� ����� ���� ����� ...

Pith tools

Reviewed August 5, 2026 · model on record in the stance chip above.