Stress Analysis of a Square Elastic Body Under Biaxial Loading Using Airy Stress Functions

This study presents an analytical investigation of stress distributions in square elastic bodies under uniaxial and biaxial compressive loads by deriving closed-form Airy stress function solutions that satisfy the biharmonic equation and boundary conditions, showing strong agreement with experimental photoelastic data.

The Gist

Imagine you have a perfect, square piece of rubber. Now, imagine someone pushes down hard on the top and bottom edges, or perhaps pushes on all four sides at once. What happens inside that rubber? Does the pressure spread out evenly, or does it get squished in weird spots?

This paper is like a mathematical crystal ball that lets us see exactly what's happening inside that square rubber block without actually cutting it open or using expensive computer simulations. The authors, researchers from Tokyo University of Agriculture and Technology, used a classic mathematical tool called the Airy Stress Function to solve this puzzle.

Here is the breakdown of their work in plain English:

The Problem: Squares are Tricky

Scientists have known for a long time how to calculate stress in round objects (like a coin being squeezed). It's like solving a puzzle with a circular frame; the math flows smoothly. But when the shape is a square, the math gets messy. The corners and straight edges make it very hard to find a perfect, exact formula. Usually, engineers have to rely on computer programs (like Finite Element Analysis) that give approximate answers.

This paper says: "Let's find the exact answer for a square."

The Method: The "Stress Recipe"

To solve this, the authors used a special mathematical recipe (the Airy Stress Function). Think of this recipe as a master key that automatically balances all the forces inside the material so they don't fly apart.

1. Breaking it Down: They took the complex pressure pushing on the edges and broke it down into a series of simple waves (like ripples on a pond).
2. The Infinite Sum: They wrote a formula that adds up thousands of these tiny waves to build the total stress picture.
3. The Tuning Knob: They had to adjust the "volume" of each wave (mathematical coefficients) until the pressure on the edges matched exactly what they wanted (either a hard push or a smooth squeeze).

The Results: What They Found

1. The "Easy Mode" Check:
First, they tested their math on a simple case: pushing evenly on all sides. As expected, the stress inside was perfectly uniform. This proved their "recipe" worked correctly.

2. The "Squeeze" Test (Uniaxial Loading):
Next, they simulated pushing down only on the top and bottom (like a Brazilian nut test).

• The Surprise: In a round disk, the tension (pulling apart) in the middle is perfectly straight and even. But in a square, the authors found that the stress near the top and bottom isn't flat. Because the square has corners and flat sides, the material resists the squeeze differently, creating a "dip" or a localized change in stress right where the force is applied.
• The Proof: They compared their math to real-world photos of stressed plastic (called photoelasticity) and computer simulations. Their mathematical "crystal ball" matched the real-world photos almost perfectly.

3. The "Double Squeeze" (Biaxial Loading):
Finally, they looked at what happens when you push on the top/bottom and the left/right at the same time.

• They found that the stress inside becomes a complex mix of the two pushes. Depending on where you look inside the square, the "difference" between the strongest and weakest stress changes. It's like mixing two different colors of paint; the result depends on exactly where you are in the mix.

Why This Matters (According to the Paper)

The authors aren't claiming this will cure diseases or build new bridges tomorrow. Instead, they are providing a gold-standard reference.

• The Benchmark: Just as a ruler is needed to check if a tape measure is accurate, this exact mathematical solution is needed to check if computer simulations are working correctly.
• The Insight: It reveals hidden details about how square materials behave that round-object math misses. It shows that the shape of the object (square vs. circle) actually changes how the stress flows right under your fingers.

In short, this paper gives us a precise, exact map of the invisible forces inside a square block of material, proving that even in a simple shape, the physics can be surprisingly complex and unique.

Technical Summary

Technical Summary: Stress Analysis of a Square Elastic Body Under Biaxial Loading Using Airy Stress Functions

Problem Statement
This study addresses the analytical determination of stress distributions within a square-shaped, homogeneous, isotropic, and linearly elastic body subjected to concentrated and distributed compressive loads under uniaxial and biaxial conditions. While classical elasticity theory provides closed-form solutions for circular and elliptical domains (e.g., the diametrically loaded disk), analytical treatments for polygonal geometries, particularly squares, remain scarce due to the lack of radial symmetry and the complexity of applying realistic boundary conditions. Although Finite Element Methods (FEM) are routinely used for such problems, they are inherently approximate. This work seeks to fill the gap by providing exact or semi-analytical solutions for square domains, which are frequently encountered in practical engineering scenarios such as indirect tensile strength evaluations (e.g., the Brazilian test) and multi-axial material testing.

Methodology
The analysis is grounded in two-dimensional linear elasticity, utilizing the Airy stress function ($\phi$) method. This approach reduces the governing equilibrium equations to a single biharmonic equation ($\nabla^4\phi = 0$), which inherently satisfies equilibrium conditions in the absence of body forces.

1. Formulation: The problem is defined on a square domain of side length $2L$, centered at the origin. External normal stresses, $\sigma_1(y)$ and $\sigma_2(x)$, are applied to the lateral and vertical edges, respectively, with shear stresses vanishing on all boundaries (traction-free or frictionless conditions).
2. Series Expansion: To handle arbitrary stress distributions, the boundary conditions are expanded into cosine Fourier series. The Airy stress function is constructed as a separable solution involving polynomial terms and infinite series of hyperbolic and trigonometric functions ($\cosh$, $\sinh$, $\cos$).
3. Coefficient Determination: The unknown coefficients in the series expansion are determined by enforcing the boundary conditions. This process yields a coupled system of linear equations.
4. Numerical Stability: To address the divergence of coefficients for large indices ($n, m \gg 1$), the authors introduce rescaled variables and modified coefficients. This transformation ensures numerical stability, allowing the infinite system to be truncated at a finite order ($N_{max}$) and solved via matrix inversion.

Key Contributions

• Generalized Analytical Framework: The paper derives a closed-form, semi-analytical solution for stress fields in square elastic bodies under arbitrary symmetric biaxial loading, a geometry not readily addressed by classical separable coordinate systems.
• Convergence Analysis: The study demonstrates the convergence behavior of the series solution. It establishes that a truncation order of $N_{max} = 128$ provides sufficient numerical accuracy, with results for $N_{max}=128$ and $256$ overlapping with a relative error of $10^{-4}$.
• Validation: The method is validated against known solutions for uniform biaxial tension (recovering spatially uniform stress fields) and compared against experimental photoelastic fringe patterns and Finite Element Analysis (FEM) for concentrated loading cases.

Results

• Uniaxial Compression: Under concentrated uniaxial compression, the resulting stress fields exhibit strong agreement with previously observed photoelastic fringe patterns. The principal stress difference ($\Delta\sigma$) contours align with experimental interference fringes.
• Geometric Effects: A distinct finding is the behavior of the tensile stress component ($\sigma_{xx}$) along the loading line ($x=0$). Unlike circular domains where tensile stress remains constant, the square geometry induces a localized compression near the load application point due to sufficient resistance from the surrounding material. This results in a non-uniform $\sigma_{xx}$ profile near the edge, a feature not captured in circular domain solutions.
• Biaxial Compression: For biaxial loading, the solution represents a superposition of two orthogonal compression scenarios. The analysis reveals spatial variations in the principal stress difference that depend on the specific location within the domain and the ratio of the applied loads.

Significance
The authors position this work as a valuable tool for theoretical validation, benchmarking, and educational purposes. By providing exact or semi-analytical expressions, the study offers a reference standard for validating computational models (such as FEM) and interpreting experimental data in mechanical testing setups involving square or rectangular specimens. The framework is presented as a stepping stone for future extensions to rectangular, anisotropic, or layered domains, as well as dynamic and three-dimensional elasticity problems, without claiming immediate application to specific new industrial processes beyond the scope of the analyzed loading configurations.