pith. sign in

arxiv: 2601.13478 · v2 · submitted 2026-01-20 · ⚛️ physics.class-ph · cond-mat.soft· physics.app-ph

Inverse Reconstruction of Moving Contact Loads on an Elastic Half-Space Using Prescribed Surface Displacement

Satoshi Takada , Yosuke Mori , Shintaro Hokada This is my paper

Pith reviewed 2026-05-16 13:09 UTC · model grok-4.3

classification ⚛️ physics.class-ph cond-mat.softphysics.app-ph
keywords moving contact loadselastic half-spaceinverse reconstructionGreen's functionsFourier inversioncontact pressureelastodynamicsMach number
0
0 comments X

The pith

A Fourier-domain inversion reconstructs the unknown moving contact traction from any prescribed surface displacement on an elastic half-space.

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

The authors first derive closed-form Green's functions for the full displacement and stress fields produced by a single point load that travels at constant speed across the surface of a two-dimensional elastic half-space, keeping the Mach-number dependence explicit. These solutions are then superposed to pose an inverse problem: given only the surface displacement profile, recover the distributed surface traction that produced it. The inversion is performed directly in Fourier space with a regularization term that stabilizes the high-frequency components, avoiding any iterative forward solves. When the method is applied to a rigid wheel whose geometry fixes the displacement inside the contact patch, the recovered pressure is smooth and symmetric while the interior stresses take analytic form involving dilogarithms; the principal-stress difference shows fringe-like contours whose forward-backward asymmetry grows with load speed.

Core claim

The surface traction that generates a prescribed displacement under constant-velocity moving contact is recovered exactly by inverting the Green's-function relation in Fourier space; regularization suppresses noise without iteration, and the resulting subsurface stresses for a wheel contact are expressible in closed form with dilogarithms whose principal-stress patterns become increasingly asymmetric as the Mach number rises.

What carries the argument

Analytical Green's functions for a constant-velocity point load on a linear isotropic elastic half-space, inverted via regularized Fourier transform to obtain the contact pressure directly from surface displacement.

If this is right

  • Subsurface stress fields are obtained in closed analytic form once the surface traction is recovered, without separate numerical integration.
  • For wheel-ground contact the reconstructed pressure distribution remains symmetric inside the patch even though the load is moving.
  • The difference between principal stresses forms spatial patterns analogous to photoelastic fringes whose asymmetry increases with Mach number.
  • Any prescribed surface displacement profile can be mapped to its generating traction and interior stresses by the same direct procedure.

Where Pith is reading between the lines

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

  • The same Green's-function inversion could be extended to three-dimensional contacts by replacing the two-dimensional kernels with their three-dimensional counterparts.
  • Testing the reconstruction on a static load whose exact solution is known would quantify how much regularization distorts sharp pressure edges.
  • The appearance of dilogarithms in the stress expressions suggests possible links to other moving-load problems whose integrals reduce to the same special functions.
  • Measured surface displacements from real rolling contacts could be fed into the procedure to infer in-situ contact pressures without assuming a particular pressure law.

Load-bearing premise

The half-space is linear and isotropic, the load velocity is exactly constant, and the regularization parameter recovers the true traction without appreciable distortion.

What would settle it

Apply the inversion to a known static rigid-cylinder indentation whose exact pressure distribution is available from Hertz theory; any large deviation between recovered and exact pressure would falsify the method.

read the original abstract

This study investigates the elastic response of a two-dimensional semi-infinite medium subjected to a moving surface load with a prescribed displacement profile. As a fundamental step, we derive analytical Green's functions for the displacement and stress fields generated by a point load traveling at a constant velocity along the surface, explicitly incorporating elastodynamic effects through Mach number dependence. These moving-load solutions serve as building blocks for constructing more general loading scenarios via linear superposition. Based on Green's functions, an inverse problem is formulated to reconstruct the unknown surface traction responsible for a given surface displacement. The inverse analysis is performed through a Fourier-domain inversion with regularization, which enables a direct and computationally efficient determination of the contact pressure without iterative forward simulations. This framework is applied to a rigid wheel-ground contact problem, where the imposed displacement is dictated by the wheel geometry. The reconstructed surface traction exhibits a smooth, symmetric distribution within the contact region, while the resulting subsurface stress fields are obtained in closed analytical form and involve dilogarithm functions. The principal stress difference reveals characteristic spatial patterns similar to photoelastic fringes, and their asymmetry increases with the Mach number, reflecting the dynamic nature of the moving contact.

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 / 1 minor

Summary. The paper derives analytical Green's functions for the displacement and stress fields produced by a point load moving at constant velocity on a 2D elastic half-space, with explicit Mach-number dependence. These functions are superposed to solve an inverse problem that reconstructs the unknown surface traction from a prescribed surface displacement via regularized Fourier-domain division; the method is demonstrated on a rigid-wheel contact geometry, producing closed-form subsurface stresses expressed with dilogarithms and principal-stress-difference patterns that become asymmetric with increasing Mach number.

Significance. If the regularization step is placed on a firm footing and the results are validated, the work supplies a non-iterative, semi-analytical route to moving-contact traction reconstruction and closed-form subsurface fields. Such expressions are rare in elastodynamics and could serve as benchmarks or rapid design tools in tribology and dynamic contact mechanics.

major comments (2)
  1. [Inverse analysis and wheel-contact application] The Fourier inversion step (described after the Green's-function derivation) divides the displacement spectrum by the moving-load transfer function and applies an unspecified regularizer. No selection rule (discrepancy principle, L-curve, or Mach-dependent cutoff) is stated; the manuscript only notes that the parameter is chosen to produce a 'smooth, symmetric' traction. Because the forward operator is compact, an arbitrary filter can systematically alter both amplitude and support of the recovered pressure, directly affecting all subsequent dilogarithm-based stresses.
  2. [Results section] No quantitative validation is presented against the static limit (Mach = 0), against known analytical solutions for uniform moving loads, or against independent numerical simulations. Without such benchmarks or an error analysis, it is impossible to assess how much distortion the chosen regularization introduces into the reconstructed traction and the reported principal-stress patterns.
minor comments (1)
  1. [Notation and abstract] Notation for the Mach number and the regularization parameter should be introduced once and used consistently; the abstract and main text currently employ slightly different symbols for the same quantities.

Simulated Author's Rebuttal

2 responses · 0 unresolved

We thank the referee for the constructive and detailed comments. We address each major point below and will revise the manuscript to strengthen the regularization procedure and add quantitative validation.

read point-by-point responses

  1. Referee: [Inverse analysis and wheel-contact application] The Fourier inversion step (described after the Green's-function derivation) divides the displacement spectrum by the moving-load transfer function and applies an unspecified regularizer. No selection rule (discrepancy principle, L-curve, or Mach-dependent cutoff) is stated; the manuscript only notes that the parameter is chosen to produce a 'smooth, symmetric' traction. Because the forward operator is compact, an arbitrary filter can systematically alter both amplitude and support of the recovered pressure, directly affecting all subsequent dilogarithm-based stresses.

    Authors: We agree that the regularization parameter was not selected according to an objective rule in the submitted version. The original choice was made to enforce non-negativity and symmetry consistent with the static limit while damping high-frequency oscillations. In the revision we will replace this heuristic with the L-curve criterion applied in the Fourier domain, explicitly reporting the corner-point parameter for each Mach number. We will also add a brief sensitivity study showing that the reconstructed traction support and the resulting dilogarithm stresses vary by less than 5 % for parameter values within a factor of two of the L-curve optimum. revision: yes

  2. Referee: [Results section] No quantitative validation is presented against the static limit (Mach = 0), against known analytical solutions for uniform moving loads, or against independent numerical simulations. Without such benchmarks or an error analysis, it is impossible to assess how much distortion the chosen regularization introduces into the reconstructed traction and the reported principal-stress patterns.

    Authors: We acknowledge the absence of quantitative benchmarks. The revised manuscript will contain a new validation subsection that (i) recovers the classical static rigid-wheel pressure distribution and subsurface stresses as Mach → 0, (ii) compares the forward superposition of the moving-load Green’s functions against the known closed-form solution for a uniform moving pressure, and (iii) quantifies reconstruction error under additive Gaussian noise on the prescribed displacement for several Mach numbers. These comparisons will be presented both graphically and in tabular form, thereby bounding the regularization-induced distortion. revision: yes

Circularity Check

0 steps flagged

No significant circularity in the derivation chain

full rationale

The paper derives analytical Green's functions directly from the elastodynamic equations for a point load moving at constant velocity on an elastic half-space, with explicit Mach-number dependence. These functions are then used to set up the inverse problem as a direct Fourier-domain division of the prescribed displacement spectrum by the Green's function spectrum, followed by regularization to control high-wavenumber amplification. No step reduces by construction to a fitted parameter, self-referential definition, or load-bearing self-citation; the regularization choice is presented as a practical filter rather than a prediction or uniqueness theorem imported from prior author work. The overall chain remains self-contained and independent of the target reconstruction results.

Axiom & Free-Parameter Ledger

1 free parameters · 2 axioms · 0 invented entities

The framework rests on standard linear elastodynamics for a half-space plus one tunable regularization parameter whose selection rule is not specified.

free parameters (1)
  • regularization parameter
    Stabilizes the Fourier inversion; its value or selection criterion is not given in the abstract.
axioms (2)
  • domain assumption Linear isotropic elasticity governs the semi-infinite medium
    Required to derive the Green's functions from the elastodynamic equations.
  • domain assumption Load travels at constant velocity
    Enables steady-state moving-load solutions in the Mach-number formulation.

pith-pipeline@v0.9.0 · 5511 in / 1315 out tokens · 41303 ms · 2026-05-16T13:09:13.741154+00:00 · methodology

discussion (0)

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

Lean theorems connected to this paper

Citations machine-checked in the Pith Canon. Every link opens the source theorem in the public Lean library.

What do these tags mean?
matches
The paper's claim is directly supported by a theorem in the formal canon.
supports
The theorem supports part of the paper's argument, but the paper may add assumptions or extra steps.
extends
The paper goes beyond the formal theorem; the theorem is a base layer rather than the whole result.
uses
The paper appears to rely on the theorem as machinery.
contradicts
The paper's claim conflicts with a theorem or certificate in the canon.
unclear
Pith found a possible connection, but the passage is too broad, indirect, or ambiguous to say the theorem truly supports the claim.

Reference graph

Works this paper leans on

12 extracted references · 12 canonical work pages

  1. [1]

    𝛽L𝑦 𝑥2 + 𝛽L 2𝑦2 −4𝛽L𝛽T 𝛽T𝑦 𝑥2 + 𝛽T 2𝑦2 ] ,(13c) 𝒢𝑦𝑦 (2)(𝑥,𝑦)≡ − 1 𝜋𝒟 [ (1 + 𝛽T 2)2 𝛽L𝑦 𝑥2 + 𝛽L 2𝑦2 −4𝛽L𝛽T 𝛽T𝑦 𝑥2 + 𝛽T 2𝑦2 ] , (13d) 𝒢𝑥𝑦 (2)(𝑥,𝑦)≡ − 1 𝜋𝒟 2𝛽L(1 + 𝛽T 2) × [ 𝑥 𝑥2 + 𝛽L 2𝑦2 − 𝑥 𝑥2 + 𝛽T 2𝑦2]. (13e) Using these Green’s functions, the displacement and stress fields in the elastic half -space induced by an arbitrary distributed normal surface load...

    work page

  2. [2]

    [7]. In the numerical evaluations, the dilogarithm is computed using standard built-in functions in Python or other languages, and no Gibbs -type oscillations are observed in the resulting stress fields. In photoelastic experiments, interference fringes appear as contours of the principal stress difference [8, 9, 10] 𝜎1 − 𝜎2 ≡ √(𝜎𝑥𝑥 − 𝜎𝑦𝑦) 2 + 4𝜎𝑥𝑦2 . (23...

    work page

  3. [3]

    Timoshenko, S.P., and Goodier, J.N.: 'Theory of Elasticity' (McGraw Hill, 1970)

    work page 1970

  4. [4]

    , and Tong, P

    Fung, Y.C. , and Tong, P. : 'Classical and Computational Solid Mechanics' (World Scientific, 2001)

    work page 2001

  5. [5]

    Cole, J., and Huth, J.: 'Stresses Produced in a Half Plane by Moving Loads', Journal of Applied Mechanics, 1958, 25, (4), pp. 433-436

    work page 1958

  6. [6]

    Mindlin, R.D.: 'Force at a Point in the Interior of a Semi - Infinite Solid', Physics, 1936, 7, 195-202

    work page 1936

  7. [7]

    Muro, T., and O ’Brien, J.: 'Terramechanics' (CRC Press, 2021)

    work page 2021

  8. [8]

    Johnson, K.L.: 'Contact Mechanics' (Cambridge University Press, 2012)

    work page 2012

  9. [9]

    Abramowitz, M., and Stegun, I.A.: 'Handbook of Mathematical Functions: with Formulas, Graphs, and Mathematical Tables' (Dover, 1965)

    work page 1965

  10. [10]

    Coker, G., and Filon, L.N.G.: 'A Treatise on Photo - Elasticity, 2nd edn.' (Cambridge University Press, 1957)

    work page 1957

  11. [11]

    Yokoyama, Y., Mitchell, B.R., Nassiri, A., Kinsey, B.L., Korkolis, Y.P., and Tagawa, Y.: 'Integrated photoelasticity in a soft material: phase retardation, azimuthal angle, and stress - optic coefficient', Optics and Lasers in Engineering, 2023, 161, 107335

    work page 2023

  12. [12]

    : 'A novel soil stress estimation method of wheel -soil interaction using photoelasticity', Journal of Terramechanics, 2025, 120, 101076

    Nagaoka, K., and Yoshida, Y. : 'A novel soil stress estimation method of wheel -soil interaction using photoelasticity', Journal of Terramechanics, 2025, 120, 101076

    work page 2025