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

This paper presents a complete analytical solution for the stress field inside a homogeneous, linearly elastic sphere subjected to a concentrated normal surface load, derived via elastodynamic equations and spherical harmonic expansions, with results generalized to arbitrary loading positions through rotational transformations and superposition.

The Gist

Imagine you have a perfectly round, solid rubber ball. Now, imagine someone presses a single, sharp finger right into the top of that ball with a sudden, hard poke. What happens inside the ball? Does the squish stay right under your finger, or does it ripple through the whole thing?

This paper is like a very detailed, mathematical recipe for answering that exact question. The authors, Yosuke Mori and his team, have figured out a way to calculate exactly how the stress (the internal "squeezing" and "stretching") moves and settles inside a solid ball when it gets poked at a single point.

Here is the breakdown of their work in plain language:

1. The Problem: The "Perfect" Poke

In the real world, if you poke a ball, the force spreads out. But in physics, it's hard to describe a "perfect" poke because it's infinitely small and infinitely strong at one spot. Previous math solutions worked for infinite spaces (like a giant block of rubber that goes on forever) or flat surfaces, but they struggled with a finite ball with a curved edge.

The authors wanted to solve this specific puzzle: What is the exact stress pattern inside a solid ball when a concentrated load hits the surface?

2. The Method: Listening to the Ball's "Vibrations"

Instead of just looking at the ball sitting still, the authors started by imagining the ball as a dynamic system. They treated the poke as a sudden event that sends waves rippling through the material, like dropping a pebble in a pond.

• The Waves: When you poke the ball, two types of waves shoot out:
  • P-waves (Compression waves): Like a sound wave, these squeeze the material together and move fast.
  • S-waves (Shear waves): These wiggle the material side-to-side and move slower.
• The Math Tool: They used a fancy mathematical technique called "spherical harmonics." Think of this like breaking a complex, messy sound (the stress field) into a set of pure musical notes. By figuring out the volume and pitch of each "note," they could rebuild the entire stress picture.

3. The Result: A Complete Map

The paper provides a "closed-form" solution. In simple terms, this means they didn't just give a computer code to guess the answer; they wrote down the exact mathematical formula for every single point inside the ball.

• The Static Picture: If you wait long enough for all the waves to settle down, you get a "static" picture. The authors found that the stress is incredibly high right under the poke and spreads out in a specific, predictable pattern. Interestingly, they found that the stress doesn't just stay in a straight line; it spreads out in all directions, creating a unique 3D pattern that is different from what happens in flat, 2D materials.
• The Dynamic Picture: They also showed what happens while the waves are moving. You can actually see the P-waves racing ahead, followed by the slower S-waves, and even a special wave that skims along the surface (like a ripple on a pond).

4. Why This Matters (According to the Paper)

The authors mention that this math is crucial for 3D photoelasticity.

• The Analogy: Imagine putting the ball in a special light. When you poke it, the stress inside makes the light bend and create colorful patterns (fringes), like a rainbow inside the ball.
• The Connection: Scientists use these rainbow patterns to figure out how strong the material is. However, to read the rainbow correctly, you need a perfect theoretical map of what the stress should look like. This paper provides that map. It allows researchers to check if their experiments or computer simulations are accurate by comparing their results against this "gold standard" math.

5. The "Superposition" Trick

The paper also explains how to handle more than one poke. If you poke the ball in four different places at once, you don't need to start from scratch. Because the math is linear, you can just take the solution for one poke, rotate it to match the new location, and add them all together. It's like mixing different colors of paint; you can predict the final color by knowing exactly how each individual color behaves.

Summary

In short, this paper gives us the ultimate "instruction manual" for understanding how a solid ball reacts when poked. It moves from the chaotic moment of the impact (the waves) to the calm, settled state (the static stress), providing a precise mathematical map that helps scientists verify their experiments and understand how stress concentrates inside 3D objects.

Technical Summary

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

Problem Formulation
This study addresses the boundary-value problem of determining the complete stress and displacement fields within a homogeneous, isotropic, linearly elastic solid sphere of radius $R$ subjected to a concentrated normal load on its surface. While classical solutions exist for infinite or semi-infinite domains (e.g., Kelvin and Boussinesq), exact analytical solutions for bounded three-dimensional bodies with curved boundaries remain scarce. The specific configuration involves a concentrated compressive stress $\sigma_0$ applied at a point on the surface (initially the north pole). The authors note that such a single load technically violates global force balance, inducing rigid-body motion; however, the analysis focuses strictly on the internal stress distribution, treating the rigid-body translation as a separate $n=1$ mode that does not affect the internal elastic deformation.

Methodology
The analysis is formulated within the framework of three-dimensional linearized elastodynamics. The authors employ the following methodological steps:

1. Governing Equations: The Navier–Cauchy equation is used, nondimensionalized using characteristic length ($R$) and time ($R/v_L$, where $v_L$ is the longitudinal wave speed).
2. Potential Representation: The displacement field is decomposed using scalar ($\phi$) and vector ($\chi$) potentials via Helmholtz decomposition. This reduces the vector elastodynamic equation to two scalar wave equations corresponding to longitudinal (P-wave) and transverse (S-wave) modes.
3. Spectral Expansion: The solution is expanded using spherical harmonics (Legendre polynomials $P_n(\cos \theta)$) and modified spherical Bessel functions. The coefficients of these expansions are determined explicitly by enforcing traction boundary conditions at the surface ($r=R$).
4. Laplace Transform: To solve the time-dependent problem, the Laplace transform is applied to the potentials. The resulting modified Helmholtz equations are solved in the Laplace domain, and the coefficients are derived in closed form.
5. Static Limit: The static elastic solution is rigorously obtained as the long-time limit ($t \to \infty$) of the dynamical formulation using the final value theorem of the Laplace transform.
6. Generalization: The solution for a load at the pole is generalized to arbitrary loading positions $(R, \theta_0, \phi_0)$ using rotational transformations. This allows for the systematic treatment of multiple concentrated loads via superposition.
7. Dynamic Evaluation: The inverse Laplace transform is evaluated using contour deformation and the residue theorem. This identifies contributions from poles on the imaginary axis, corresponding to the static solution, P-waves, and S-waves.

Key Contributions

• Closed-Form Analytical Solution: The paper provides explicit, closed-form expressions for all components of the stress tensor ($\sigma_{\rho\rho}, \sigma_{\theta\theta}, \sigma_{\phi\phi}, \sigma_{\rho\theta}$) and displacement fields. Unlike previous series representations that often remain in implicit or highly coupled forms, all expansion coefficients are determined explicitly.
• Unified Static-Dynamic Framework: The work unifies the static and dynamic problems, deriving the static solution as a rigorous limit of the dynamical formulation rather than as a separate static derivation.
• Explicit Mode Analysis: The analysis explicitly identifies the $n=1$ spherical harmonic mode as representing rigid-body translation, which can be excluded when calculating internal stresses. This clarifies a point often treated implicitly in other analytical treatments.
• Principal Stress Difference: Closed-form expressions enable the direct calculation of principal stresses and their differences throughout the sphere's interior, a quantity critical for three-dimensional photoelasticity.

Results

• Static Stress Distribution: The principal stress difference is found to be highest near the loading point, exhibiting a divergent behavior proportional to $(1-\rho)^{-2}$ as the surface is approached. Along the loading axis, the radial stress is compressive, while the tangential stresses are tensile, reflecting lateral expansion.
• Dynamic Wave Propagation: The dynamic solution captures the propagation of P-waves and S-waves with distinct velocities. The results show the emergence of:
  • Primary P-waves and S-waves propagating concentrically from the source.
  • Rayleigh waves localized near the surface.
  • Reflected waves (PP-waves) and mode-converted waves (PS-waves) generated by boundary interactions.
  • Complex envelope structures formed by the superposition of waves emitted at different times, consistent with two-dimensional analogs.
• Multiple Loads: The superposition principle is demonstrated for configurations with multiple loads (e.g., four loads balancing each other), showing regions of stress amplification and cancellation.
• Validation: The analytical results are compared with Finite Element Method (FEM) simulations. The comparison shows good overall agreement, with minor discrepancies near the surface attributed to mesh resolution limitations in representing concentrated loads in FEM.

Significance
The paper claims that its primary significance lies in providing a rigorous theoretical reference for interpreting experimental observations in three-dimensional photoelasticity. By offering explicit closed-form expressions for the full stress tensor and principal stress differences, the solution serves as a benchmark for validating numerical methods (such as FEM) and experimental reconstruction techniques. The unified framework for both static and dynamic problems, combined with the ability to handle arbitrary loading positions through rotation, offers a systematic tool for analyzing stress concentrations in bounded spherical elastic bodies. The authors position this work as a foundational step for understanding stress singularities and wave propagation in finite elastic domains.