arxiv: 2605.03256 · v1 · submitted 2026-05-05 · ⚛️ physics.class-ph · cond-mat.soft

Recognition: unknown

Revisiting the Stress Field Inside an Elastic Sphere Subjected to a Concentrated Load

Yosuke Mori , Kiwamu Yoshii , Satoshi Takada

Authors on Pith no claims yet

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

classification ⚛️ physics.class-ph cond-mat.soft

keywords elastic sphereconcentrated loadstress fieldanalytical solutionspherical harmonicselastodynamicsstatic limitrotational transformation

0 comments

The pith

A complete analytical solution with closed-form expressions for every stress component is derived for an elastic sphere under a concentrated surface load.

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

The paper establishes a full analytical solution for the stress field throughout the interior of a homogeneous linearly elastic sphere when a concentrated normal load acts at one point on the surface. The derivation begins with the three-dimensional elastodynamic equations, expresses the displacement using scalar and vector potentials, expands in spherical harmonics, and fixes the coefficients from the surface traction conditions. The static solution is recovered exactly as the long-time limit of the dynamic problem, yielding explicit formulas for all stress-tensor components that can be evaluated directly at any interior point. These formulas also extend to loads at arbitrary surface locations through rotational transformations and to multiple loads by superposition. A reader would care because the result supplies exact principal stresses and stress differences without discretization, replacing approximate numerical methods for stress analysis inside spherical bodies.

Core claim

We present a complete analytical solution for the stress field inside a homogeneous, linearly elastic solid sphere subjected to a concentrated normal load applied on its surface. Starting from the three-dimensional linearized elastodynamic equations, the displacement and stress fields are derived using scalar and vector potential representations combined with spherical harmonic expansions. All expansion coefficients are determined explicitly by enforcing the traction boundary conditions. The static elastic solution is obtained rigorously as the long-time limit of the dynamical formulation. Closed-form expressions for all components of the stress tensor are provided, enabling direct evaluati

What carries the argument

Scalar and vector potential representations of the displacement field combined with spherical harmonic expansions, with all coefficients fixed by surface traction boundary conditions, followed by extraction of the static fields as the long-time limit of the dynamic solution.

Load-bearing premise

The sphere is homogeneous and linearly elastic, and the static solution is recovered exactly as the long-time limit of the elastodynamic formulation.

What would settle it

Numerical evaluation of the closed-form stress expressions at an interior point, followed by verification that they satisfy the elastostatic equilibrium equations and reproduce the applied concentrated traction on the surface.

read the original abstract

We present a complete analytical solution for the stress field inside a homogeneous, inside a homogeneous, linearly elastic solid sphere subjected to a concentrated normal load applied on its surface. Starting from the three-dimensional linearized elastodynamic equations, the displacement and stress fields are derived using scalar and vector potential representations combined with spherical harmonic expansions. All expansion coefficients are determined explicitly by enforcing the traction boundary conditions. The static elastic solution is obtained rigorously as the long-time limit of the dynamical formulation. Closed-form expressions for all components of the stress tensor are provided, enabling direct evaluation of the principal stresses and their differences throughout the interior of the sphere. The analytical solution is further generalized to arbitrary loading positions by means of rotational transformations, allowing systematic treatment of multiple concentrated loads through superposition.

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.

Desk Editor's Note private letter to a colleague

This paper gives explicit closed-form stresses inside a point-loaded elastic sphere by solving the dynamic problem and taking the long-time limit, with a rotational generalization for arbitrary positions.

read the letter

The main takeaway is that the authors derive a full set of closed-form expressions for every stress component inside a homogeneous elastic sphere under a concentrated normal surface load. They begin with the three-dimensional elastodynamic equations, introduce scalar and vector potentials, expand in spherical harmonics, fix all coefficients from the traction boundary conditions, and extract the static solution as the long-time limit. They then rotate the result to handle loads at any surface point and note that superposition covers multiple loads. This produces usable analytical formulas for principal stresses and their differences throughout the interior. The explicit coefficients and the rotational step are the concrete advances over earlier static treatments, which often stopped at less complete forms or required numerical evaluation. The approach follows standard potential methods and delivers something practical for contact or soft-matter calculations where interior fields matter. The soft spot is the long-time limit step. In an undamped elastic sphere, waves reflect repeatedly from the boundary, so displacements and stresses generally keep oscillating and do not converge pointwise to the static equilibrium. The abstract calls the limit rigorous, but without damping, time-averaging, or an explicit convergence argument the static expressions could be formal rather than strictly obtained from the dynamic equations. If the manuscript supplies a proper justification or shows that the oscillating parts cancel in the stress formulas, that concern disappears; otherwise it is the part that needs checking. The rest of the derivation looks standard and free of circular fitting or invented terms. This is a paper for researchers who need analytical interior stresses in spherical geometries rather than a broad audience. A reader already working on similar boundary-value problems in elasticity will find the explicit forms worth having, even if they re-derive the limit themselves. It is grounded enough to go to referees, who can verify the coefficient algebra and the handling of the limit. I would send it out for review rather than desk-reject.

Referee Report

1 major / 1 minor

Summary. The paper claims to present a complete analytical solution for the stress field inside a homogeneous, linearly elastic solid sphere subjected to a concentrated normal load on its surface. It derives the displacement and stress fields from the three-dimensional linearized elastodynamic equations using scalar and vector potential representations combined with spherical harmonic expansions, determines all expansion coefficients explicitly from the traction boundary conditions, obtains the static solution rigorously as the long-time limit of the dynamical formulation, provides closed-form expressions for all stress tensor components, and generalizes the result to arbitrary load positions via rotational transformations and superposition.

Significance. If the long-time limit is rigorously established, the closed-form stress expressions would provide a valuable analytical benchmark for numerical elastostatic and elastodynamic codes in spherical geometries, as well as direct means to evaluate principal stresses and stress differences without discretization. The dynamical-to-static approach and the rotational generalization for multiple loads are practical strengths that could facilitate applications in contact mechanics and solid mechanics.

major comments (1)

[Abstract and long-time limit derivation] Abstract and the section deriving the static solution from the dynamical formulation: The claim that the static elastic solution is obtained 'rigorously' as the long-time limit of the undamped elastodynamic equations is load-bearing for the central result. In a bounded spherical domain without viscous damping, solutions to the elastic wave equation generally consist of persistent oscillations from repeated boundary reflections, so the pointwise limit lim_{t→∞} of the stress field does not exist in general. The manuscript must specify the mathematical device used to justify the limit (e.g., explicit damping term, time-averaging procedure, or appeal to a convergence theorem) and confirm that the resulting expressions satisfy div σ = 0 together with the given traction boundary conditions. Without this justification the static closed-form expressions cannot be regarded as proven limits of a

minor comments (1)

[Abstract] Abstract contains a duplicated phrase ('inside a homogeneous, inside a homogeneous') that should be corrected for clarity.

Simulated Author's Rebuttal

1 responses · 0 unresolved

We thank the referee for the careful reading and for recognizing the potential utility of the closed-form stress expressions as benchmarks. We address the concern about the long-time limit below.

read point-by-point responses

Referee: [Abstract and long-time limit derivation] Abstract and the section deriving the static solution from the dynamical formulation: The claim that the static elastic solution is obtained 'rigorously' as the long-time limit of the undamped elastodynamic equations is load-bearing for the central result. In a bounded spherical domain without viscous damping, solutions to the elastic wave equation generally consist of persistent oscillations from repeated boundary reflections, so the pointwise limit lim_{t→∞} of the stress field does not exist in general. The manuscript must specify the mathematical device used to justify the limit (e.g., explicit damping term, time-averaging procedure, or appeal to a convergence theorem) and confirm that the resulting expressions satisfy div σ = 0 together with the given traction boundary conditions. Without this justification the static closed-form expressions

Authors: We agree that the pointwise limit does not exist for the undamped problem in a bounded domain. The manuscript obtains the static solution by first introducing a small viscous damping term in the elastodynamic equations, solving the resulting damped initial-boundary-value problem via the same spherical-harmonic expansion, taking the long-time limit for fixed damping, and finally letting the damping coefficient tend to zero. This procedure is implicit in the coefficient determination but was not stated explicitly. We will add a dedicated paragraph in the derivation section that (i) writes the damped equations, (ii) shows that the time-dependent coefficients decay exponentially, (iii) extracts the static coefficients, and (iv) verifies by direct substitution that the resulting stress field satisfies both div σ = 0 in the interior and the prescribed traction boundary conditions on the sphere surface. revision: yes

Circularity Check

0 steps flagged

No circularity: derivation from standard elastodynamics with explicit coefficient determination from BCs

full rationale

The paper starts from the three-dimensional linearized elastodynamic equations, represents fields via scalar/vector potentials and spherical-harmonic expansions, and determines all expansion coefficients explicitly by enforcing the traction boundary conditions on the sphere surface. The static solution is obtained by taking the long-time limit of this dynamical formulation. No step reduces a claimed prediction or result to a fitted parameter, self-definition, or load-bearing self-citation by construction; the central expressions are derived directly from the governing PDEs and BCs without renaming known results or smuggling ansatzes via prior self-work. This is the most common honest non-finding for papers that solve a boundary-value problem from first principles.

Axiom & Free-Parameter Ledger

0 free parameters · 2 axioms · 0 invented entities

The central claim rests on the standard assumptions of linear elasticity and the validity of the Helmholtz decomposition into scalar and vector potentials for the elastodynamic equations, plus the mathematical legitimacy of taking the long-time limit to recover the static solution.

axioms (2)

domain assumption The solid is homogeneous and linearly elastic
Explicitly stated as the starting point for the three-dimensional linearized elastodynamic equations.
domain assumption The static solution is the long-time limit of the dynamical formulation
Used to obtain the static elastic solution from the time-dependent problem.

pith-pipeline@v0.9.0 · 5427 in / 1300 out tokens · 42379 ms · 2026-05-09T16:47:30.740293+00:00 · methodology

Review history (2 revisions) →

discussion (0)

Reference graph

Works this paper leans on

24 extracted references

[1]

McGraw-Hill, New York (1970)

Timoshenko, S.P., Goodier, J.N.: Theory of Elasticity, 3 rd edn. McGraw-Hill, New York (1970)

1970
[2]

II Li near Theory, Academic Press, Cambridge (1975)

Eringen, A.C., S ¸uhubi, E.S.: Elastodynamics, Vol. II Li near Theory, Academic Press, Cambridge (1975)

1975
[3]

University Science Books, Sausalito (2002)

Aki, K., Richards, P.G.: Quantitative Seismology, 2nd Ed . University Science Books, Sausalito (2002)

2002
[4]

Fung, Y.C., Tong, P.: Classical and Computational Solid M echanics, vol. 1. W orld Sci- entiﬁc, Singapore (2001)

2001
[5]

Sternberg, E., Rosenthal, F.: The Elastic Sphere Under Co ncentrated Loads. J. Appl. Mech. 19(4), 413 (1952)

1952
[6]

Oka, Y.: Determination of the tensile stren gth of rock by a compression test of an irregular test piece

Hiramatsu, Y. Oka, Y.: Determination of the tensile stren gth of rock by a compression test of an irregular test piece. Int. J. Rock Mech. Min. Sci. G eomech. Abst. 3(2), 89 (1966)

1966
[7]

Jingu, T., Nezu, K.: Transient stress in an elastic sphere under diametrical concentrated impact loads. Bull. JSME 28(245), 2553 (1985)

1985
[8]

Sato, Y., Takada, S.: Revisiting stress propagation in a t hree-dimensional elastic sphere under diametric loading. Bull. JSME 90(933), p. 23-00262 (2024)

2024
[9]

Powder Technol

Sch¨ onert, K.: Breakage of spheres and circular discs. Powder Technol. 143–144 , 2 (2004)

2004
[10]

Ma, C.C., Hung, K.M.: Exact full-ﬁeld analysis of strain and displacement for circular disks subjected to partially distributed compressions. In t. J. Mech. Sci. 50(2), 275 (2008)

2008
[11]

Ramesh, K., Shins, K.: Stress ﬁeld equations for a disk su bjected to self-equilibrated arbitrary loads: revisited. Granul. Matter 24(2), 49 (2022)

2022
[12]

Shins, K., Ramesh, K.: Closed-form displacement ﬁeld eq uations for a disc subjected to self-equilibrated arbitrary loads. Granul. Matter 25(2), 38 (2023)

2023
[13]

W u, S.Z., Chau, K.T.: Dynamic response of an elastic sphe re under diametral impacts. Mech. Mater. 38(11), 1039 (2006)

2006
[14]

Guerrero, I., Turteltaub, M.J.: The elastic sphere unde r arbitrary concentrated surface loads. J. Elast. 2, 21 (1972)

1972
[15]

Frocht, M.M.: Photoelasticity, vol. 2. Wiley, New York ( 1941)

1941
[16]

Cambridge Uni- versity Press, Cambridge (1957)

Coker, E.G., Filon, L.N.G.: A Treatise on Photo-Elastic ity, 2nd edn. Cambridge Uni- versity Press, Cambridge (1957)

1957
[17]

Yokoyama, Y., Mitchell, B.R., Nassiri, A., Kinsey, B. L. , Korkolis, Y. P., and Tagawa, Y.: Integrated photoelasticity in a soft material: phase re tardation, zimuthal angle, and stress-optic coeﬃcient: Opt. Lasers Eng. 161, 197335 (2023)

2023
[18]

EPJ W eb Conf

Yu, P., Taghizadeh, K., Feisel, D.D., Ganguly, S., Schr¨ oter, M., Sperl, M.: Simulating the 3D photoelasticity forward problem in order to generate training images for deep learning. EPJ W eb Conf. 340, 10019 (2025)

2025
[19]

Dover Publications, New York (1965)

Abramowitz, M., Stegun, I.A.: Graphs, and Mathematical Tables. Dover Publications, New York (1965)

1965
[20]

Cambridge Univers ity Press, Cambridge (2007)

Press, W.H., Teukolsky, S.A., Vetterling, W.T., Flanne ry, B.P.: Numerical Recipes 3rd Edition: The Art of Scientiﬁc Computing. Cambridge Univers ity Press, Cambridge (2007)

2007
[21]

Jingu, T., Hisada, K., Nakahara, I., Machida, S.: Transi ent stress in a circular disk under diametrical impact loads. Bull. JSME 28(235), 13 (1985)

1985
[22]

Sato, Y., Ishikawa, H., Takada, S.; Revisiting Stress Pr opagation in a Two-Dimensional Elastic Circular Disk Under Diametric Loading. J. Elast. 156, 193 (2024)

2024
[23]

Strength Mater

Okamura, K., Sato, Y., Takada, S.: Displacement and Stre ss Analysis of an Elastic Hollow Disk: Comparison with Strength of Materials’ Predic tion. Strength Mater. 57, 102 (2025)

2025
[24]

Hughes, T. J. R.: The Finite Element Method: Linear Stati c and Dynamic Finite Ele- ment Analysis. Dover Publications (2012)

2012