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
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.
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
- 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.
Referee Report
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)
- [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.
- [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)
- The abstract would be clearer if it named the specific fractional derivative order or relaxation function employed in the simulations.
- 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
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
-
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
-
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
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
axioms (1)
- domain assumption Linear viscoelasticity governs the material response at the asperity scale
Reference graph
Works this paper leans on
-
[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)
work page 2005
-
[2]
Bachelier, Libraire, Paris (1821)
Coulomb,C.A.:ThéorieDesMachinesSimples:EnAyantÉgardAuFrottementdeLeurs Parties et À La Roideur Des Cordages. Bachelier, Libraire, Paris (1821)
-
[3]
Dieterich, J.H.: Time-dependent friction in rocks. J. Geophys. Res.77(20), 3690–3697 (1972) https://doi.org/10.1029/JB077i020p03690
-
[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]
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]
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]
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
work page 2018
-
[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
work page 2006
-
[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
work page 1960
-
[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
work page doi:10.1115/1 1966
-
[11]
Cambridge University Press, Cambridge (1985)
Johnson, K.L.: Contact Mechanics. Cambridge University Press, Cambridge (1985). https://doi.org/10.1017/CBO9781139171731 18
-
[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]
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]
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
work page 2012
-
[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]
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
work page doi:10.1111/j 1967
-
[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]
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]
World Scientific, Singapore (2022)
Mainardi, F.: Fractional Calculus and Waves in Linear Viscoelasticity: An Introduction to Mathematical Models. World Scientific, Singapore (2022)
work page 2022
-
[20]
Lakes,R.S.:ViscoelasticMaterials.CambridgeUniversityPress,Cambridge;NewYork (2009)
work page 2009
-
[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
work page 1983
-
[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]
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]
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]
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]
Vernier, P.: Homogenization of composite materials with fractional viscoelastic con- stituents. PhD thesis (2023)
work page 2023
-
[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)
work page 2002
-
[28]
Li,Z.,Brenner,R.,Frérot,L.:Supplementarycodesanddatato"Evolutionofthecontact betweenroughviscoelasticsolidsafterdecreasingloads:Memoryerasureandmonotonic increase" (2025) https://doi.org/10.5281/zenodo.17671850
-
[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]
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]
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]
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]
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]
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]
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]
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]
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]
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]
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]
Sharpe,W.N.(ed.):SpringerHandbookofExperimentalSolidMechanics.SpringerLink Bücher. Springer, Boston, MA (2008). https://doi.org/10.1007/978-0-387-30877-7
-
[41]
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]
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]
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
discussion (0)
Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.