pith. sign in

arxiv: 2511.19402 · v2 · pith:FKT5RKUFnew · submitted 2025-11-24 · ❄️ cond-mat.soft

Evolution of the contact between rough viscoelastic solids after decreasing loads: memory erasure and monotonic increase

Pith reviewed 2026-05-25 07:48 UTC · model grok-4.3

classification ❄️ cond-mat.soft
keywords viscoelastic contactrough surfacescontact area evolutionunloadingmemory effectsfractional viscoelasticitylogarithmic aging
0
0 comments X

The pith

Linear viscoelastic models erase memory upon unloading and cannot produce decreasing contact area after load reduction.

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

Experiments show that after decreasing the normal load on rough PMMA contacts, the true contact area decreases and retains long-term memory of the prior state. Fractional viscoelastic models, adapted to capture wide relaxation spectra and logarithmic aging under constant load, instead erase that memory immediately: the area evolves exactly as if the reduced load had always been applied. A general proof establishes that no linear viscoelastic model, fractional or otherwise, can produce a decreasing phase in contact area after unloading. The authors therefore conclude that additional local internal variables are required to reproduce both the observed memory and the post-unload reduction.

Core claim

By adapting existing contact theories and numerical methods to fractional viscoelasticity, the models reproduce logarithmic aging at constant load yet show that memory of the prior contact state is erased upon unloading, with the area behaving as if the lower load had always been present. None of the simulations exhibit a decreasing regime after load reduction; this monotonic increase holds for all linear viscoelastic models despite their ability to capture aging. Additional local internal variables are therefore required to explain both long-term contact memory and contact area reduction after unloading.

What carries the argument

Adaptation of rough contact theories and simulations to fractional viscoelastic constitutive relations, which provide a broad relaxation spectrum while remaining linear.

If this is right

  • Logarithmic aging under constant load is captured by fractional models.
  • Memory of the contact state prior to unloading is erased for all linear viscoelastic models.
  • Contact area increases monotonically after unloading in any linear viscoelastic model.
  • Both the long-term memory and the decreasing phase observed in experiments require additional local internal variables.

Where Pith is reading between the lines

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

  • Friction force predictions that rely solely on linear viscoelastic contact area will miss history-dependent changes after load reduction.
  • Similar memory-erasure behavior may appear in related problems such as viscoelastic wear or adhesion under varying loads.
  • Implementing specific nonlinear or internal-variable constitutive models offers a direct route to test whether they recover the experimental decreasing phase.

Load-bearing premise

That the real area of contact evolution for linear viscoelastic materials is fully represented by adapting existing contact theories and numerical methods to fractional viscoelasticity.

What would settle it

A linear viscoelastic rough-contact simulation or analytical solution that produces a decreasing contact-area phase after a sudden load reduction.

read the original abstract

The real area of contact governs, in part, the magnitude of the friction force, yet its time evolution in rough viscoelastic interfaces remains incompletely understood. In experiments of contact between polymethylmethacrylate blocks under decreasing normal loads, Dillavou and Rubinstein have shown that the true contact area exhibits, after unloading, a decreasing phase and long-term memory of the contact state prior to unloading. It is however unclear what modeling ingredients are necessary to reproduce these two features. Here, we investigate these effects using fractional viscoelastic rough contact models. By adapting existing contact theories and numerical simulation methods to fractional viscoelasticity, which induces a wide relaxation spectrum, we reproduce logarithmic aging under constant load, but show that memory of the contact state is erased upon unloading. Indeed, the contact area behaves as if it had always experienced the reduced load, even on short time-scales, contrasting with the response of a standard linear solid. Moreover, none of our results show a decreasing regime of the contact area after unload: we ultimately prove that this is the case for all linear viscoelastic models -- despite capturing logarithmic aging -- leading to the conclusion that additional local internal variables are required to explain both long-term contact memory and contact area reduction after unloading.

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 investigates the time evolution of real contact area between rough viscoelastic solids under decreasing normal loads. Using adaptations of existing rough-contact theories and solvers to fractional viscoelasticity (which captures logarithmic aging under constant load), the authors find that upon unloading the contact area behaves as if the interface had always experienced only the reduced load, erasing prior memory on all timescales. They report no decreasing regime in contact area after unloading in their simulations and provide a general proof that this holds for all linear viscoelastic models, concluding that additional local internal variables are required to explain both long-term memory and post-unload area reduction seen in experiments.

Significance. If the general proof is valid, the result is significant for contact mechanics and tribology: it shows that linear viscoelasticity, even when extended to broad relaxation spectra via fractional operators, is insufficient to capture key memory effects observed experimentally, thereby identifying a clear modeling gap and motivating the inclusion of local internal variables. The explicit general proof (rather than case-by-case numerics) is a methodological strength.

major comments (2)
  1. [Abstract (final paragraph)] Abstract (final paragraph): the central claim that 'we ultimately prove that this is the case for all linear viscoelastic models' is load-bearing, yet the manuscript provides no outline of the proof steps, assumptions, or handling of the nonlinear Signorini condition that determines the pressure field and thus the area; without these details the extrapolation from the fractional numerics to arbitrary relaxation functions cannot be assessed.
  2. [Numerical methods (implied)] Numerical adaptation section (implied by abstract): the claim that fractional viscoelasticity is faithfully represented inside existing iterative contact solvers rests on unstated details of kernel discretization and multi-scale asperity coarse-graining; because area is a nonlinear functional of the evolving pressure, any inconsistency here would invalidate both the reported absence of a decreasing regime and the general proof.
minor comments (2)
  1. The abstract would be clearer if it named the specific fractional derivative order or relaxation function employed in the simulations.
  2. Ensure all citations to prior rough-contact solvers (e.g., the base theories being adapted) appear explicitly in the methods description.

Simulated Author's Rebuttal

2 responses · 0 unresolved

We thank the referee for their positive assessment of the work's significance and for the constructive comments. We address each major comment below and will revise the manuscript to improve clarity on the requested points.

read point-by-point responses
  1. Referee: [Abstract (final paragraph)] Abstract (final paragraph): the central claim that 'we ultimately prove that this is the case for all linear viscoelastic models' is load-bearing, yet the manuscript provides no outline of the proof steps, assumptions, or handling of the nonlinear Signorini condition that determines the pressure field and thus the area; without these details the extrapolation from the fractional numerics to arbitrary relaxation functions cannot be assessed.

    Authors: The complete proof appears in the main text, but we agree that the abstract and introduction lack an explicit outline of steps, assumptions, and treatment of the Signorini condition. The argument uses linearity of the viscoelastic convolution to show that post-unload contact area depends only on the instantaneous load via superposition on the pressure field; the unilateral constraint is incorporated by demonstrating that any admissible pressure history consistent with the reduced load yields the same area evolution. We will insert a concise outline of these elements into the abstract and add a short proof-summary subsection. revision: yes

  2. Referee: [Numerical methods (implied)] Numerical adaptation section (implied by abstract): the claim that fractional viscoelasticity is faithfully represented inside existing iterative contact solvers rests on unstated details of kernel discretization and multi-scale asperity coarse-graining; because area is a nonlinear functional of the evolving pressure, any inconsistency here would invalidate both the reported absence of a decreasing regime and the general proof.

    Authors: We accept that the numerical-methods description would benefit from explicit statements on kernel discretization (via Prony-series approximation of the fractional operator) and the multi-scale coarse-graining procedure. These choices preserve the underlying linear viscoelastic integral equation at every scale, which was verified against analytic limits. We will expand the relevant section with the discretization scheme, convergence checks, and coarse-graining rules to ensure full transparency. revision: yes

Circularity Check

0 steps flagged

No significant circularity; general proof stands independently of numerics

full rationale

The paper's central result is an explicit general proof that no post-unload decreasing contact-area regime exists for any linear viscoelastic model (despite logarithmic aging under constant load). This proof is presented as extending beyond the fractional case studied numerically, and the adaptation of existing contact solvers is described as a technical extension rather than a fitted or self-defining step. No equation reduces by construction to a prior fit, no uniqueness theorem is imported from the authors' own prior work, and no ansatz is smuggled via self-citation. The derivation chain therefore remains self-contained against external benchmarks and does not exhibit any of the enumerated circular patterns.

Axiom & Free-Parameter Ledger

0 free parameters · 1 axioms · 0 invented entities

The paper relies on the domain assumption of linear viscoelasticity and adaptations of prior contact theories without introducing new free parameters, axioms beyond standard math, or invented entities.

axioms (1)
  • domain assumption Linear viscoelasticity governs the material response at the asperity scale
    Invoked throughout the abstract when considering all linear viscoelastic models and their inability to produce the observed features.

pith-pipeline@v0.9.0 · 5752 in / 1249 out tokens · 37607 ms · 2026-05-25T07:48:52.332310+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

43 extracted references · 43 canonical work pages

  1. [1]

    Lubrica- tion Engineering61, 18–22 (2005)

    Jost, H.P.: Tribology micro & macro economics: A road to economic savings. Lubrica- tion Engineering61, 18–22 (2005)

  2. [2]

    Bachelier, Libraire, Paris (1821)

    Coulomb,C.A.:ThéorieDesMachinesSimples:EnAyantÉgardAuFrottementdeLeurs Parties et À La Roideur Des Cordages. Bachelier, Libraire, Paris (1821)

  3. [3]

    Dieterich, J.H.: Time-dependent friction in rocks. J. Geophys. Res.77(20), 3690–3697 (1972) https://doi.org/10.1029/JB077i020p03690

  4. [4]

    pure and applied geophysics143(1-3), 283–302 (1994) https://doi.org/10.1007/BF00874332

    Dieterich, J.H., Kilgore, B.D.: Direct observation of frictional contacts: New insights for state-dependent properties. pure and applied geophysics143(1-3), 283–302 (1994) https://doi.org/10.1007/BF00874332

  5. [5]

    Li, Z., Pastewka, L., Szlufarska, I.: Chemical aging of large-scale randomly rough fric- tional contacts. Phys. Rev. E98(2), 023001 (2018) https://doi.org/10.1103/PhysRevE. 98.023001

  6. [6]

    PhD thesis (2020) https://doi.org/10.5075/epfl-thesis-7640

    Frérot, L.H.G.: Bridging scales in wear modeling with volume integral methods for elastic-plastic contact. PhD thesis (2020) https://doi.org/10.5075/epfl-thesis-7640

  7. [7]

    Physical Review Letters120(22), 224101 (2018) https://doi.org/10.1103/ PhysRevLett.120.224101

    Dillavou, S., Rubinstein, S.M.: Nonmonotonic Aging and Memory in a Frictional Interface. Physical Review Letters120(22), 224101 (2018) https://doi.org/10.1103/ PhysRevLett.120.224101

  8. [8]

    Advances in Physics55(3-4), 279–348 (2006) https://doi.org/10.1080/ 00018730600732186

    Baumberger, T., Caroli, C.: Solid friction from stick–slip down to pinning and aging. Advances in Physics55(3-4), 279–348 (2006) https://doi.org/10.1080/ 00018730600732186

  9. [9]

    Lee, E.H., Radok, J.R.M.: The contact problem for viscoelastic bodies. Jour- nal of Applied Mechanics27(3), 438–444 (1960) https://doi.org/10.1115/ 1.3644020 https://asmedigitalcollection.asme.org/appliedmechanics/article- pdf/27/3/438/5443475/438\_1.pdf

  10. [10]

    Journal of Applied Mechanics33(4), 845–854 (1966) https://doi.org/10.1115/1

    Ting, T.C.T.: The Contact Stresses Between a Rigid Indenter and a Viscoelastic Half- Space. Journal of Applied Mechanics33(4), 845–854 (1966) https://doi.org/10.1115/1. 3625192

  11. [11]

    Cambridge University Press, Cambridge (1985)

    Johnson, K.L.: Contact Mechanics. Cambridge University Press, Cambridge (1985). https://doi.org/10.1017/CBO9781139171731 18

  12. [12]

    The Journal of Chemical Physics115(8), 3840–3861 (2001) https://doi.org/10.1063/1.1388626

    Persson, B.N.J.: Theory of rubber friction and contact mechanics. The Journal of Chemical Physics115(8), 3840–3861 (2001) https://doi.org/10.1063/1.1388626

  13. [13]

    Papangelo, A., Ciavarella, M.: Viscoelastic normal indentation of nominally flat ran- domlyroughcontacts.InternationalJournalofMechanicalSciences211,106783(2021) https://doi.org/10.1016/j.ijmecsci.2021.106783

  14. [14]

    Proceedings oftheNationalAcademyofSciences109(6),1850–1855(2012)https://doi.org/10.1073/ pnas.1120147109

    Amir, A., Oreg, Y., Imry, Y.: On relaxations and aging of various glasses. Proceedings oftheNationalAcademyofSciences109(6),1850–1855(2012)https://doi.org/10.1073/ pnas.1120147109

  15. [15]

    The American Journal of Psychology56(2), 234 (1943) https://doi.org/10.2307/1417506 1417506

    Blair, G.W.S., Coppen, F.M.V.: The Estimation of Firmness in Soft Materials. The American Journal of Psychology56(2), 234 (1943) https://doi.org/10.2307/1417506 1417506

  16. [16]

    A., Hekker, S., Stello, D., Guti ´errez-Soto, J., Handberg, R., Huber, D., et al

    Caputo, M.: Linear Models of Dissipation whose Q is almost Frequency Independent– II. Geophysical Journal International13(5), 529–539 (1967) https://doi.org/10.1111/j. 1365-246X.1967.tb02303.x

  17. [17]

    pure and applied geophysics91(1), 134–147 (1971) https://doi.org/10.1007/BF00879562

    Caputo, M., Mainardi, F.: A new dissipation model based on memory mechanism. pure and applied geophysics91(1), 134–147 (1971) https://doi.org/10.1007/BF00879562

  18. [18]

    The European Physical Journal Special Topics193(1), 133–160 (2011) https://doi.org/10.1140/epjst/e2011-01387-1

    Mainardi, F., Spada, G.: Creep, relaxation and viscosity properties for basic fractional models in rheology. The European Physical Journal Special Topics193(1), 133–160 (2011) https://doi.org/10.1140/epjst/e2011-01387-1

  19. [19]

    World Scientific, Singapore (2022)

    Mainardi, F.: Fractional Calculus and Waves in Linear Viscoelasticity: An Introduction to Mathematical Models. World Scientific, Singapore (2022)

  20. [20]

    Lakes,R.S.:ViscoelasticMaterials.CambridgeUniversityPress,Cambridge;NewYork (2009)

  21. [21]

    Journal of Rheology27(3), 201–210 (1983) https://doi.org/10.1122/ 1.549724

    Bagley,R.L.,Torvik,P.J.:ATheoreticalBasisfortheApplicationofFractionalCalculus to Viscoelasticity. Journal of Rheology27(3), 201–210 (1983) https://doi.org/10.1122/ 1.549724

  22. [22]

    Journal of Applied Mechanics51(2), 294–298 (1984) https://doi

    Torvik,P.J.,Bagley,R.L.:OntheAppearanceoftheFractionalDerivativeintheBehav- ior of Real Materials. Journal of Applied Mechanics51(2), 294–298 (1984) https://doi. org/10.1115/1.3167615

  23. [23]

    Poly- mer Testing21(3), 325–331 (2002) https://doi.org/10.1016/S0142-9418(01)00092-7

    Hernández-Jiménez,A.,Hernández-Santiago,J.,Macias-García,A.,Sánchez-González, J.: Relaxation modulus in PMMA and PTFE fitting by fractional Maxwell model. Poly- mer Testing21(3), 325–331 (2002) https://doi.org/10.1016/S0142-9418(01)00092-7

  24. [24]

    Polymer Testing30(7), 797–799 (2011) https://doi.org/10.1016/j.polymertesting.2011

    Hu, K.-X., Zhu, K.-Q.: A note on fractional Maxwell model for PMMA and PTFE. Polymer Testing30(7), 797–799 (2011) https://doi.org/10.1016/j.polymertesting.2011. 06.008 19

  25. [25]

    International Journal of Solids and Structures67–68, 169–181 (2015) https://doi.org/10.1016/j.ijsolstr.2015.04.018

    Jalocha, D., Constantinescu, A., Neviere, R.: Revisiting the identification of general- ized Maxwell models from experimental results. International Journal of Solids and Structures67–68, 169–181 (2015) https://doi.org/10.1016/j.ijsolstr.2015.04.018

  26. [26]

    PhD thesis (2023)

    Vernier, P.: Homogenization of composite materials with fractional viscoelastic con- stituents. PhD thesis (2023)

  27. [27]

    Fractional Calculus and Applied Analysis (2002)

    Gorenflo, R., Loutchko, J., Luchko, Y.: Computation of the Mittag-Leffler function and its derivatives. Fractional Calculus and Applied Analysis (2002)

  28. [28]

    Evolutionofthecontact betweenroughviscoelasticsolidsafterdecreasingloads:Memoryerasureandmonotonic increase

    Li,Z.,Brenner,R.,Frérot,L.:Supplementarycodesanddatato"Evolutionofthecontact betweenroughviscoelasticsolidsafterdecreasingloads:Memoryerasureandmonotonic increase" (2025) https://doi.org/10.5281/zenodo.17671850

  29. [29]

    Bugnicourt,R.,Sainsot,P.,Lesaffre,N.,Lubrecht,A.A.:Transientfrictionlesscontactof aroughrigidsurfaceonaviscoelastichalf-space.TribologyInternational113,279–285 (2017) https://doi.org/10.1016/j.triboint.2017.01.032

  30. [30]

    Wear 231(2), 206–219 (1999) https://doi.org/10.1016/S0043-1648(99)00113-1

    Polonsky, I.A., Keer, L.M.: A numerical method for solving rough contact problems based on the multi-level multi-summation and conjugate gradient techniques. Wear 231(2), 206–219 (1999) https://doi.org/10.1016/S0043-1648(99)00113-1

  31. [31]

    Vollebregt,E.A.H.:Anewsolverfortheelasticnormalcontactproblemusingconjugate gradients,deflation,andanFFT-basedpreconditioner.JournalofComputationalPhysics 257, 333–351 (2014) https://doi.org/10.1016/j.jcp.2013.10.005

  32. [32]

    Journal of Tribology119(3), 481–485 (1997) https://doi.org/10.1115/1.2833523

    Stanley, H.M., Kato, T.: An FFT-Based Method for Rough Surface Contact. Journal of Tribology119(3), 481–485 (1997) https://doi.org/10.1115/1.2833523

  33. [33]

    JOSS5(51), 2121 (2020) https://doi

    Frérot, L., Anciaux, G., Rey, V., Pham-Ba, S., Molinari, J.-F.: Tamaas: A library for elastic-plastic contact of periodic rough surfaces. JOSS5(51), 2121 (2020) https://doi. org/10.21105/joss.02121

  34. [34]

    Physical Review E86(3), 035601 (2012) https://doi.org/10.1103/PhysRevE.86

    Yastrebov,V.A.,Anciaux,G.,Molinari,J.-F.:Contactbetweenrepresentativeroughsur- faces. Physical Review E86(3), 035601 (2012) https://doi.org/10.1103/PhysRevE.86. 035601

  35. [35]

    Tribology International33(1), 47–58 (2000) https://doi.org/10.1016/S0301-679X(00)00016-5

    Wu, J.-J.: Simulation of rough surfaces with fft. Tribology International33(1), 47–58 (2000) https://doi.org/10.1016/S0301-679X(00)00016-5

  36. [36]

    Tribology International 114, 161–171 (2017) https://doi.org/10.1016/j.triboint.2017.04.023

    Yastrebov, V.A., Anciaux, G., Molinari, J.-F.: On the accurate computation of the true contact-area in mechanical contact of random rough surfaces. Tribology International 114, 161–171 (2017) https://doi.org/10.1016/j.triboint.2017.04.023

  37. [37]

    International Journal of Solids and Structures52, 83–102 (2015) https://doi.org/10.1016/j.ijsolstr.2014.09.019 20

    Yastrebov,V.A.,Anciaux,G.,Molinari,J.-F.:Frominfinitesimaltofullcontactbetween rough surfaces: Evolution of the contact area. International Journal of Solids and Structures52, 83–102 (2015) https://doi.org/10.1016/j.ijsolstr.2014.09.019 20

  38. [38]

    International Journal of Mechanical Sciences52(6), 829–835 (2010) https://doi

    Greenwood, J.A.: Contact between an axisymmetric indenter and a viscoelastic half- space. International Journal of Mechanical Sciences52(6), 829–835 (2010) https://doi. org/10.1016/j.ijmecsci.2010.01.010

  39. [39]

    Physical Review Letters118(8), 085501 (2017) https://doi.org/10.1103/PhysRevLett.118.085501

    Lahini,Y.,Gottesman,O.,Amir,A.,Rubinstein,S.M.:NonmonotonicAgingandMem- ory Retention in Disordered Mechanical Systems. Physical Review Letters118(8), 085501 (2017) https://doi.org/10.1103/PhysRevLett.118.085501

  40. [40]

    Springer, Boston, MA (2008)

    Sharpe,W.N.(ed.):SpringerHandbookofExperimentalSolidMechanics.SpringerLink Bücher. Springer, Boston, MA (2008). https://doi.org/10.1007/978-0-387-30877-7

  41. [41]

    F1000Research (2021)

    Mölder, F., Jablonski, K.P., Letcher, B., Hall, M.B., Tomkins-Tinch, C.H., Sochat, V., Forster, J., Lee, S., Twardziok, S.O., Kanitz, A., Wilm, A., Holtgrewe, M., Rahmann, S., Nahnsen, S., Köster, J.: Sustainable Data Analysis with Snakemake. F1000Research (2021). https://doi.org/10.12688/f1000research.29032.2

  42. [42]

    D., 2007, @doi [Computing in Science and Engineering] 10.1109/MCSE.2007.55 , 9, 90

    Hunter, J.D.: Matplotlib: A 2D Graphics Environment. Computing in Science & Engi- neering9(3), 90–95 (2007) https://doi.org/10.1109/MCSE.2007.55

  43. [43]

    Papoulia, K.D., Panoskaltsis, V.P., Kurup, N.V., Korovajchuk, I.: Rheological represen- tationoffractionalorderviscoelasticmaterialmodels.RheolActa49(4),381–400(2010) https://doi.org/10.1007/s00397-010-0436-y 21