pith. sign in

arxiv: 2604.18917 · v1 · submitted 2026-04-20 · 💻 cs.CE · cs.NA· math.NA· physics.comp-ph

An Implicit Compact-Kernel Material Point Method for Computational Solid Mechanics

Qirui Fu , Yupeng Jiang , Minchen Li This is my paper

Pith reviewed 2026-05-10 02:37 UTC · model grok-4.3

classification 💻 cs.CE cs.NAmath.NAphysics.comp-ph
keywords material point methodcompact kernelimplicit integrationcomputational solid mechanicscell-crossing instabilitycontact accuracylarge deformation
0
0 comments X

The pith

Implicit compact-kernel MPM combines compact support with the smoothness needed for stable large-deformation solid mechanics

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

The paper develops an implicit formulation of the compact-kernel material point method and tests it on benchmark problems including cantilever bending, Hertzian contact, narrow-clearance free fall, and colliding hyperelastic rings. It establishes that the method keeps the locality advantages of compact support while delivering the smoothness required for robust implicit large-deformation runs. A sympathetic reader would care because conventional MPM kernels force trade-offs between suppressing cell-crossing noise, controlling numerical diffusion, and achieving accurate contact without artificial gaps. The reported comparisons show reduced stress noise relative to linear MPM and better contact locality relative to quadratic B-spline MPM while holding overall accuracy comparable.

Core claim

The authors present an implicit CK-MPM and demonstrate through the benchmarks that it retains the advantages of compact support while preserving the smoothness required for robust large-deformation simulation. Compared with linear MPM it reduces cell-crossing-induced stress noise and excessive numerical dissipation; compared with quadratic B-spline MPM it improves contact locality and reduces artificial contact gaps and early-contact artifacts while maintaining comparable overall smoothness and accuracy. These results indicate that CK-MPM provides a viable implicit MPM framework for computational mechanics.

What carries the argument

The compact kernel formulated inside an implicit time-stepping scheme that governs particle-to-grid and grid-to-particle transfers

If this is right

  • Implicit CK-MPM reduces cell-crossing-induced stress noise and excessive numerical dissipation relative to linear MPM.
  • Implicit CK-MPM improves contact locality and reduces artificial contact gaps and early-contact artifacts relative to quadratic B-spline MPM.
  • Implicit CK-MPM maintains comparable overall smoothness and accuracy for large-deformation cases.
  • CK-MPM supplies a workable implicit framework for computational solid mechanics simulations.

Where Pith is reading between the lines

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

  • The formulation could be tested on frictional or multi-body contact problems to check whether the locality gains persist under more complex interaction rules.
  • Lower diffusion might permit coarser background grids in practice, producing computational savings once the kernel is fixed.
  • The same compact kernel could be inserted into other implicit integrators or coupled with different constitutive models without changing the core transfer logic.

Load-bearing premise

The compact kernel can be placed inside an implicit time-stepping scheme without introducing new instabilities or demanding extra stabilization parameters that would erase the reported gains in locality and reduced diffusion.

What would settle it

A set of the same benchmark runs in which implicit CK-MPM either develops instabilities that require new stabilization parameters or fails to show the stated reductions in cell-crossing noise and contact artifacts would disprove the central claim.

Figures

Figures reproduced from arXiv: 2604.18917 by Minchen Li, Qirui Fu, Yupeng Jiang.

Figure 1
Figure 1. Figure 1: Comparison of the three kernel functions used in this study. The compact kernel combines the principal advantages of the linear and quadratic kernels, possessing a compact support width of 1 and being differentiable everywhere. In MPM, the shape function needs to satisfy two properties at every material point to ensure 1st-order accurate interpolation: X i ωip = 1, (21) X i xiωip = xp. (22) 6 [PITH_FULL_I… view at source ↗
Figure 2
Figure 2. Figure 2: Particle–grid connectivity for the three kernel functions. In two dimensions, each particle interacts with 4 grid nodes for the linear kernel and 9 grid nodes for the quadratic kernel. For the compact kernel, each particle interacts with 4 nodes per grid, and two staggered grids are required to satisfy Equation 24. Based on this, CK-MPM advance the system by transferring particle information to the two gri… view at source ↗
Figure 3
Figure 3. Figure 3: Each pair of two nodes associated with the same particle corresponds to a non-zero block in the Hessian matrix. In 2D, with quadratic kernel, one node may share particles with 24 neighboring nodes at most. With compact kernel, one node only share particles with 8 neighboring nodes in the same grid and 16 nodes in the other grid. In 3D, these numbers become 124, 26, and 64, respectively, indicating a sparse… view at source ↗
Figure 4
Figure 4. Figure 4: 2D cantilever beam configuration. The beam is clamped at the left end and undergoes bending under its self￾weight. Here, L and h denote the beam length and thickness, respectively, while W and H represent the horizontal and vertical displacements of the free end. The computational setup, as shown in [PITH_FULL_IMAGE:figures/full_fig_p012_4.png] view at source ↗
Figure 5
Figure 5. Figure 5: Comparison of the bending ratio H/W versus γ obtained with MPM and CK-MPM. The results are nearly indistinguishable, showing that CK-MPM maintains numerical smoothness and predictive accuracy comparable to those of quadratic B-spline MPM. For both methods, the numerical results lie within the analytical asymptotic bounds. 4.2. Hertz Contact Between Cylinder and Rigid Plane We next consider a two-dimensiona… view at source ↗
Figure 6
Figure 6. Figure 6: The setup of Hertz contact experiment. Conceptually, it is a cylinder with infinite length put on a rigid plane. We only conduct experiment near the contact area to reduce computational costs. analytical solution. The analytical distribution of the contact pressure is given as p = −σ max yy r 1 − ( x a ) 2, for − a ≤ x ≤ a, (49) where a = 2r F R(1 − ν 2) πE , σmax yy = 2F πa . (50) The contact-pressure dis… view at source ↗
Figure 7
Figure 7. Figure 7: Pressure distribution along the contact interface predicted by CK-MPM for different grid spacings ∆x, compared with the analytical solution. As ∆x decreases, the numerical results progressively converge toward the analytical pressure distribution. solution. This observation further supports the interpretation that the improved performance of CK-MPM originates from its enhanced transfer locality, showing th… view at source ↗
Figure 8
Figure 8. Figure 8: Comparison of pressure distributions between CK-MPM and MPM. As shown in the figure, for each ∆x, the distribution obtained with CK-MPM is closer to the theoretical solution. This improvement is attributed to its smaller support radius. support is desirable because it helps preserve transfer locality and reduce artificial smoothing. Accordingly, CK￾MPM is expected to exhibit improved performance in suppres… view at source ↗
Figure 9
Figure 9. Figure 9: Relative RMSE as a function of ∆x for CK-MPM and quadratic MPM. Both methods exhibit a similar convergence trend with refinement, while CK-MPM consistently yields lower errors over the range of ∆x considered. MPM x=0.1 mm MPM x=0.05 mm MPM x=0.025 mm MPM x=0.0125 mm CK-MPM x=0.1 mm CK-MPM x=0.05 mm CK-MPM x=0.025 mm CK-MPM x=0.0125 mm 220.0 188.6 157.1 125.7 94.3 62.9 31.4 0.0 yy (KPa) [PITH_FULL_IMAGE:fi… view at source ↗
Figure 10
Figure 10. Figure 10: Global stress distributions σyy obtained with MPM and CK-MPM at different grid spacings ∆x. The two methods produce comparable global stress fields across all resolutions, indicating that the proposed method preserves the overall stress distribution while providing improved accuracy in the contact region. the benchmark is shown in [PITH_FULL_IMAGE:figures/full_fig_p017_10.png] view at source ↗
Figure 11
Figure 11. Figure 11: Configuration of the numerical diffusion experiment. (a) The values on either side of the interface are different. After transferring iterations a numerical diffusion band with width δ will appear. In this experiment we use a = 1 m. (b) The red solid line denotes the interface: particles above it are assigned φ = 1, while those below are assigned φ = −1. For visualization purposes, the grid spacing in thi… view at source ↗
Figure 12
Figure 12. Figure 12: Comparison of numerical diffusion after 50 transfer iterations for different kernels. It is obvious that the width of transition band in CK-MPM is smaller than that in other two kernels. 0.0 0.2 0.4 0.6 0.8 1.0 y (m) −1.25 −1.00 −0.75 −0.50 −0.25 0.00 0.25 0.50 0.75 1.00 1.25 φ (a) 0 20 40 60 80 100 Transferring Iteration 0.0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1.0 δ (m) Linear MPM Quadratic MPM CK-MPM (b… view at source ↗
Figure 13
Figure 13. Figure 13: (a) One-dimensional slice of particle values along x = 0.5 after 50 transfer iterations. Though compact kernel has a slightly larger supporting radius with linear kernel, it can reduce diffusion significantly compared with linear kernel. It benefits from more concentration around center area. (b) Numerical diffusion growth in our pure diffusion benchmark. CK-MPM exhibits the slowest growth rate, followed … view at source ↗
Figure 14
Figure 14. Figure 14: Schematic of the sphere-fall test through a hollow cylinder. The cylinder inner radius is set to be slightly larger than the sphere radius; therefore, the sphere should theoretically pass through the cylinder without contact-induced obstruction. In the numerical tests, CK-MPM reproduces this behavior accurately, whereas quadratic MPM fails to do so [PITH_FULL_IMAGE:figures/full_fig_p020_14.png] view at source ↗
Figure 15
Figure 15. Figure 15: Time history of the sphere center height. Owing to its larger kernel support, quadratic MPM does not correctly reproduce the passage of the sphere through the hollow cylinder. By contrast, the CK-MPM results closely follow the free￾fall trajectory, indicating that the sphere motion remains essentially unaffected by spurious contact forces throughout the process. 5.3. Impact of Two Neo-Hookean Rings We fur… view at source ↗
Figure 16
Figure 16. Figure 16: Comparison of the von Mises stress distributions. The sphere computed with MPM exhibits substantially higher stress levels than that computed with CK-MPM, reflecting the stronger spurious contact forces generated between the sphere and the hollow cylinder in MPM. is included as a baseline for transfer-induced energy dissipation, while the quadratic B-spline formulation serves as a representative smooth wi… view at source ↗
Figure 17
Figure 17. Figure 17: Initial configuration of the collision problem between two Neo-Hookean rings moving toward each other. The initial gap is set to 2.02 m. The stress fields of the two colliding rings at t = 0.9 s, 1.8 s, 2.7 s, and 3.6 s are presented in [PITH_FULL_IMAGE:figures/full_fig_p022_17.png] view at source ↗
Figure 18
Figure 18. Figure 18: Von Mises stress distributions at selected times. Compared with linear MPM, CK-MPM results in smooth stress distributions without noise. At the same time, it reduces the early-contact artifacts observed in quadratic MPM, thereby more accurately resolving the collision process. We further examine the energy conservation characteristics of the three methods. In the absence of artificial numerical dissipatio… view at source ↗
Figure 19
Figure 19. Figure 19: Evolution of the potential, kinetic, and total energies. The energy response of CK-MPM closely follows that of quadratic MPM, indicating that replacing the linear kernel with the compact kernel substantially reduces the excessive energy dissipation observed in linear MPM. 6. Conclusion In this work, we developed an implicit formulation of CK-MPM and examined its performance in several benchmark problems r… view at source ↗
read the original abstract

The numerical performance of the material point method (MPM) is strongly governed by the particle-grid kernel, which controls the trade-off among smoothness, locality, numerical diffusion, contact accuracy, and computational cost. Although wide-support smooth kernels can effectively suppress cell-crossing instability, they often introduce increased numerical diffusion, artificial contact gaps, and higher transfer cost. In contrast, the suitability of compact-kernel designs for implicit computational solid mechanics remains unclear. In this work, we develop an implicit formulation of the Compact-Kernel Material Point Method (CK-MPM) and assess its performance through benchmark problems in linear and nonlinear solid mechanics, including cantilever bending, Hertzian contact, narrow-clearance free fall, and colliding hyperelastic rings. The results show that implicit CK-MPM retains the advantages of compact support while preserving the smoothness required for robust large-deformation simulation. Compared with linear MPM, it reduces cell-crossing-induced stress noise and excessive numerical dissipation; compared with quadratic B-spline MPM, it improves contact locality and reduces artificial contact gaps and early-contact artifacts while maintaining comparable overall smoothness and accuracy. These results indicate that CK-MPM provides a viable implicit MPM framework for computational mechanics.

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

2 major / 2 minor

Summary. The manuscript develops an implicit time-stepping formulation of the Compact-Kernel Material Point Method (CK-MPM) for computational solid mechanics. It claims that this approach preserves the locality and reduced numerical diffusion of compact-support kernels while retaining the smoothness needed for stable large-deformation simulations. Through benchmarks including cantilever bending, Hertzian contact, narrow-clearance free fall, and colliding hyperelastic rings, the work reports reduced cell-crossing stress noise relative to linear MPM and improved contact locality with fewer artificial gaps relative to quadratic B-spline MPM, without introducing new stabilization parameters.

Significance. If the central claims hold, the implicit CK-MPM offers a practical advance by resolving the smoothness-locality trade-off in MPM kernels for implicit solid-mechanics simulations. The benchmarks provide concrete evidence that compact kernels can be embedded in Newton-based time integrators while maintaining robustness, which could reduce computational cost and improve accuracy in contact-dominated or large-strain problems. The parameter-free nature of the kernel design is a notable strength.

major comments (2)
  1. [Results section] Results section (Hertzian contact and colliding rings benchmarks): the reported reductions in artificial contact gaps and early-contact artifacts are described qualitatively. Quantitative metrics such as contact-force error norms, measured gap distances relative to analytical solutions, or L2 displacement errors with grid refinement would be required to substantiate the superiority over quadratic B-spline MPM and to confirm that the implicit residual evaluation preserves locality without introducing new artifacts.
  2. [Section 3] Section 3 (implicit formulation): the construction of the tangent Jacobian and residual for the compact kernel is central to the claim that no additional stabilization is needed. A brief analysis of Newton iteration counts, matrix conditioning, or solver convergence behavior across the deformation regimes would directly address whether the limited kernel support alters stability or requires implicit-specific adjustments that could offset the reported gains.
minor comments (2)
  1. [Figures] Figure captions and legends should explicitly label all compared methods (linear MPM, quadratic B-spline, CK-MPM) and indicate whether results are at the same time step or deformation level.
  2. [Introduction] The abstract and introduction would benefit from a short statement of the specific compact kernel function (e.g., its support radius and polynomial degree) to allow immediate comparison with existing literature.

Simulated Author's Rebuttal

2 responses · 0 unresolved

We thank the referee for the constructive review and the recommendation of minor revision. We address each major comment below and outline the revisions we will make to strengthen the quantitative support for our claims.

read point-by-point responses

  1. Referee: [Results section] Results section (Hertzian contact and colliding rings benchmarks): the reported reductions in artificial contact gaps and early-contact artifacts are described qualitatively. Quantitative metrics such as contact-force error norms, measured gap distances relative to analytical solutions, or L2 displacement errors with grid refinement would be required to substantiate the superiority over quadratic B-spline MPM and to confirm that the implicit residual evaluation preserves locality without introducing new artifacts.

    Authors: We agree that the current presentation of contact improvements is primarily qualitative. In the revised manuscript we will add quantitative metrics to the Hertzian contact and colliding rings benchmarks, including L2 displacement error norms under successive grid refinement, measured gap distances relative to the analytical Hertz solution, and contact-force error norms. These additions will enable direct, quantitative comparison with quadratic B-spline MPM and will confirm that the implicit residual evaluation preserves locality without new artifacts. revision: yes

  2. Referee: [Section 3] Section 3 (implicit formulation): the construction of the tangent Jacobian and residual for the compact kernel is central to the claim that no additional stabilization is needed. A brief analysis of Newton iteration counts, matrix conditioning, or solver convergence behavior across the deformation regimes would directly address whether the limited kernel support alters stability or requires implicit-specific adjustments that could offset the reported gains.

    Authors: The tangent Jacobian and residual are assembled directly from the compact kernel without auxiliary stabilization, following the standard Newton-Raphson procedure described in Section 3. To address the referee’s request, we will insert a concise convergence study (Newton iteration counts and solver residual histories) for the cantilever, Hertzian, and ring-collision benchmarks. This will show that the limited support does not degrade conditioning or require implicit-specific adjustments beyond those already present in the formulation. revision: yes

Circularity Check

0 steps flagged

No significant circularity detected; derivation is self-contained

full rationale

The paper introduces an implicit formulation of CK-MPM and supports its claims through direct numerical benchmarks (cantilever bending, Hertzian contact, narrow-clearance free fall, colliding rings) that compare against linear MPM and quadratic B-spline MPM. No load-bearing step reduces by construction to a fitted parameter, self-citation chain, or renamed input; the central performance assertions rest on independent simulation outcomes rather than tautological definitions or prior author results invoked as uniqueness theorems. The method development and validation chain remains externally falsifiable via the reported benchmark suite.

Axiom & Free-Parameter Ledger

0 free parameters · 0 axioms · 0 invented entities

Based on the abstract alone, the method builds on standard MPM particle-to-grid transfer and implicit time integration; no new free parameters, axioms, or invented entities are explicitly introduced or quantified in the provided text.

pith-pipeline@v0.9.0 · 5518 in / 1207 out tokens · 29274 ms · 2026-05-10T02:37:46.245296+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

58 extracted references · 58 canonical work pages

  1. [1]

    ACM Transactions on Graphics (TOG) , volume=

    Ck-mpm: A compact-kernel material point method , author=. ACM Transactions on Graphics (TOG) , volume=. 2025 , publisher=

    work page 2025

  2. [2]

    ACM Transactions on Graphics (TOG) , volume=

    Physical validation of simulators in Computer Graphics: A new framework dedicated to slender elastic structures and frictional contact , author=. ACM Transactions on Graphics (TOG) , volume=. 2021 , publisher=

    work page 2021

  3. [3]

    ACM SIGGRAPH 2023 Conference Proceedings , pages=

    Multi-layer thick shells , author=. ACM SIGGRAPH 2023 Conference Proceedings , pages=

    work page 2023

  4. [4]

    , author=

    Energetically Consistent Invertible Elasticity. , author=. Symposium on Computer Animation , volume=

    work page

  5. [5]

    The heavy elastica , author=

    L. The heavy elastica , author=. The London, Edinburgh, and Dublin Philosophical Magazine and Journal of Science , volume=. 1934 , publisher=

    work page 1934

  6. [6]

    The Quarterly Journal of Mechanics and Applied Mathematics , volume=

    Bending of a beam or wide strip , author=. The Quarterly Journal of Mechanics and Applied Mathematics , volume=. 1992 , publisher=

    work page 1992

  7. [7]

    Advances in applied mechanics , volume=

    Material point method after 25 years: Theory, implementation, and applications , author=. Advances in applied mechanics , volume=. 2020 , publisher=

    work page 2020

  8. [8]

    Computer Modeling in Engineering & Sciences , volume=

    The generalized interpolation material point method , author=. Computer Modeling in Engineering & Sciences , volume=. 2004 , publisher=

    work page 2004

  9. [9]

    Computer methods in applied mechanics and engineering , volume=

    The material-point method for granular materials , author=. Computer methods in applied mechanics and engineering , volume=. 2000 , publisher=

    work page 2000

  10. [10]

    Computer methods in applied mechanics and engineering , volume=

    The material point method in large strain engineering problems , author=. Computer methods in applied mechanics and engineering , volume=. 2004 , publisher=

    work page 2004

  11. [11]

    Nature materials , volume=

    Intrusion rheology in grains and other flowable materials , author=. Nature materials , volume=. 2016 , publisher=

    work page 2016

  12. [12]

    Journal of Fluid Mechanics , volume=

    A general fluid--sediment mixture model and constitutive theory validated in many flow regimes , author=. Journal of Fluid Mechanics , volume=. 2019 , publisher=

    work page 2019

  13. [13]

    Acm siggraph 2016 courses , pages=

    The material point method for simulating continuum materials , author=. Acm siggraph 2016 courses , pages=

    work page 2016

  14. [14]

    International journal of impact engineering , volume=

    Comparison study of MPM and SPH in modeling hypervelocity impact problems , author=. International journal of impact engineering , volume=. 2009 , publisher=

    work page 2009

  15. [15]

    ACM transactions on graphics , year=

    Incremental potential contact: intersection-and inversion-free, large-deformation dynamics , author=. ACM transactions on graphics , year=

    work page

  16. [16]

    Codimensional incremental potential contact

    Codimensional incremental potential contact , author=. arXiv preprint arXiv:2012.04457 , year=

    work page arXiv 2012

  17. [17]

    ACM Transactions on Graphics (TOG) , volume=

    Anisompm: Animating anisotropic damage mechanics , author=. ACM Transactions on Graphics (TOG) , volume=. 2020 , publisher=

    work page 2020

  18. [18]

    ACM Transactions on Graphics (TOG) , volume=

    CD-MPM: continuum damage material point methods for dynamic fracture animation , author=. ACM Transactions on Graphics (TOG) , volume=. 2019 , publisher=

    work page 2019

  19. [19]

    International Journal for Numerical and Analytical Methods in Geomechanics , volume=

    Evaluation of material point method for use in geotechnics , author=. International Journal for Numerical and Analytical Methods in Geomechanics , volume=. 2015 , publisher=

    work page 2015

  20. [20]

    International Journal for Numerical Methods in Engineering , volume=

    A hybrid material-point spheropolygon-element method for solid and granular material interaction , author=. International Journal for Numerical Methods in Engineering , volume=. 2020 , publisher=

    work page 2020

  21. [21]

    Acta Geotechnica , volume=

    Hybrid continuum--discrete simulation of granular impact dynamics , author=. Acta Geotechnica , volume=. 2022 , publisher=

    work page 2022

  22. [22]

    Computers and Geotechnics , volume=

    Coupled material point and level set methods for simulating soils interacting with rigid objects with complex geometry , author=. Computers and Geotechnics , volume=. 2023 , publisher=

    work page 2023

  23. [23]

    Computer Methods in Applied Mechanics and Engineering , volume=

    A barrier method for frictional contact on embedded interfaces , author=. Computer Methods in Applied Mechanics and Engineering , volume=. 2022 , publisher=

    work page 2022

  24. [24]

    Computers and Geotechnics , volume=

    Impact and erosion dynamics of inclusion-enriched dry granular flows on rigid barriers with basal clearance: Numerical insights from hybrid MP-DEM simulations , author=. Computers and Geotechnics , volume=. 2024 , publisher=

    work page 2024

  25. [25]

    Journal of Geophysical Research: Earth Surface , volume=

    Erosion and transport of dry soil bed by collisional granular flow: Insights from a combined experimental--numerical investigation , author=. Journal of Geophysical Research: Earth Surface , volume=. 2023 , publisher=

    work page 2023

  26. [26]

    Computers and Geotechnics , volume=

    Revealing the role of forests in the mobility of geophysical flows , author=. Computers and Geotechnics , volume=. 2023 , publisher=

    work page 2023

  27. [27]

    International Journal for Numerical Methods in Engineering , volume=

    An explicit phase field material point method for modeling dynamic fracture problems , author=. International Journal for Numerical Methods in Engineering , volume=. 2023 , publisher=

    work page 2023

  28. [28]

    Journal of Computational Physics , volume=

    An angular momentum conserving affine-particle-in-cell method , author=. Journal of Computational Physics , volume=. 2017 , publisher=

    work page 2017

  29. [29]

    International Journal for numerical methods in Engineering , volume=

    A convected particle domain interpolation technique to extend applicability of the material point method for problems involving massive deformations , author=. International Journal for numerical methods in Engineering , volume=. 2011 , publisher=

    work page 2011

  30. [30]

    International Journal for numerical methods in Engineering , volume=

    Second-order convected particle domain interpolation (CPDI2) with enrichment for weak discontinuities at material interfaces , author=. International Journal for numerical methods in Engineering , volume=. 2013 , publisher=

    work page 2013

  31. [31]

    ACM Transactions on Graphics (TOG) , volume=

    A moving least squares material point method with displacement discontinuity and two-way rigid body coupling , author=. ACM Transactions on Graphics (TOG) , volume=. 2018 , publisher=

    work page 2018

  32. [32]

    Computer Methods in Applied Mechanics and Engineering , volume=

    DEM-enriched contact approach for material point method , author=. Computer Methods in Applied Mechanics and Engineering , volume=. 2023 , publisher=

    work page 2023

  33. [33]

    International Journal for numerical methods in engineering , volume=

    Enhancement of the material point method using B-spline basis functions , author=. International Journal for numerical methods in engineering , volume=. 2018 , publisher=

    work page 2018

  34. [34]

    International Journal for Numerical Methods in Engineering , volume=

    Extended B-spline-based implicit material point method , author=. International Journal for Numerical Methods in Engineering , volume=. 2021 , publisher=

    work page 2021

  35. [35]

    Computational Mechanics , volume=

    A contact method for B-spline material point method with application in impact and penetration problems , author=. Computational Mechanics , volume=. 2024 , publisher=

    work page 2024

  36. [36]

    International Journal for Numerical Methods in Engineering , volume=

    Validating High-Performance Multi-GPU MPM for Debris-Fluid-Structure Interaction , author=. International Journal for Numerical Methods in Engineering , volume=. 2025 , publisher=

    work page 2025

  37. [37]

    Computational Particle Mechanics , volume=

    Nitsche-based material point method for large deformation frictional contact problems , author=. Computational Particle Mechanics , volume=. 2025 , publisher=

    work page 2025

  38. [38]

    Computer methods in applied mechanics and engineering , volume=

    ILS-MPM: An implicit level-set-based material point method for frictional particulate contact mechanics of deformable particles , author=. Computer methods in applied mechanics and engineering , volume=. 2020 , publisher=

    work page 2020

  39. [39]

    Computer Methods in Applied Mechanics and Engineering , volume=

    A mortar segment-to-segment frictional contact approach in material point method , author=. Computer Methods in Applied Mechanics and Engineering , volume=. 2024 , publisher=

    work page 2024

  40. [40]

    Computers & Structures , volume=

    B-spline based boundary conditions in the material point method , author=. Computers & Structures , volume=. 2019 , publisher=

    work page 2019

  41. [41]

    International Journal for Numerical Methods in Engineering , volume=

    An Extended B-Spline-Based Material Point Method for Contact Problems , author=. International Journal for Numerical Methods in Engineering , volume=. 2025 , publisher=

    work page 2025

  42. [42]

    Computer Methods in Applied Mechanics and Engineering , volume=

    DP-MPM: Domain partitioning material point method for evolving multi-body thermal--mechanical contacts during dynamic fracture and fragmentation , author=. Computer Methods in Applied Mechanics and Engineering , volume=. 2021 , publisher=

    work page 2021

  43. [43]

    Proceedings of the Institution of Mechanical Engineers , volume=

    One hundred years of Hertz contact , author=. Proceedings of the Institution of Mechanical Engineers , volume=. 1982 , publisher=

    work page 1982

  44. [44]

    International Journal for Numerical Methods in Engineering , volume=

    Field-gradient partitioning for fracture and frictional contact in the material point method , author=. International Journal for Numerical Methods in Engineering , volume=. 2017 , publisher=

    work page 2017

  45. [45]

    Computer physics communications , volume=

    Application of a particle-in-cell method to solid mechanics , author=. Computer physics communications , volume=. 1995 , publisher=

    work page 1995

  46. [46]

    International journal for numerical methods in engineering , volume=

    Analysis and reduction of quadrature errors in the material point method (MPM) , author=. International journal for numerical methods in engineering , volume=. 2008 , publisher=

    work page 2008

  47. [47]

    IEEE transactions on visualization and computer graphics , volume=

    Optimization integrator for large time steps , author=. IEEE transactions on visualization and computer graphics , volume=. 2015 , publisher=

    work page 2015

  48. [48]

    ACM Transactions on Graphics (TOG) , volume=

    Hierarchical optimization time integration for cfl-rate mpm stepping , author=. ACM Transactions on Graphics (TOG) , volume=. 2020 , publisher=

    work page 2020

  49. [49]

    International Journal for numerical methods in engineering , volume=

    Variational integrators and the Newmark algorithm for conservative and dissipative mechanical systems , author=. International Journal for numerical methods in engineering , volume=. 2000 , publisher=

    work page 2000

  50. [50]

    2006 , publisher=

    Numerical optimization , author=. 2006 , publisher=

    work page 2006

  51. [51]

    ACM Transactions on Graphics (TOG) , volume=

    VR-Doh: Hands-on 3D Modeling in Virtual Reality , author=. ACM Transactions on Graphics (TOG) , volume=. 2025 , publisher=

    work page 2025

  52. [52]

    ACM SIGGRAPH 2024 Conference Papers , pages=

    A dynamic duo of finite elements and material points , author=. ACM SIGGRAPH 2024 Conference Papers , pages=

    work page 2024

  53. [53]

    Computational Mechanics , volume=

    Unstructured moving least squares material point methods: a stable kernel approach with continuous gradient reconstruction on general unstructured tessellations , author=. Computational Mechanics , volume=. 2025 , publisher=

    work page 2025

  54. [54]

    International Journal for Numerical and Analytical Methods in Geomechanics , volume=

    Mapped material point method for large deformation problems with sharp gradients and its application to soil-structure interactions , author=. International Journal for Numerical and Analytical Methods in Geomechanics , volume=. 2024 , publisher=

    work page 2024

  55. [55]

    International Journal for Numerical Methods in Engineering , volume=

    Lagrangian--Eulerian multidensity topology optimization with the material point method , author=. International Journal for Numerical Methods in Engineering , volume=. 2021 , publisher=

    work page 2021

  56. [56]

    SIGGRAPH Asia 2023 Conference Papers , pages=

    Neural stress fields for reduced-order elastoplasticity and fracture , author=. SIGGRAPH Asia 2023 Conference Papers , pages=

    work page 2023

  57. [57]

    ACM Transactions on Graphics (TOG) , volume=

    A massively parallel and scalable multi-GPU material point method , author=. ACM Transactions on Graphics (TOG) , volume=. 2020 , publisher=

    work page 2020

  58. [58]

    ACM Transactions on Graphics (TOG) , volume=

    Energetically consistent inelasticity for optimization time integration , author=. ACM Transactions on Graphics (TOG) , volume=. 2022 , publisher=

    work page 2022