Displacement general solutions in strain gradient elasticity: review and analysis

This paper reviews and extends general displacement solutions for static isotropic strain gradient elasticity by demonstrating that all classical elasticity representations can be generalized through Helmholtz decomposition, while establishing their consistency with prior work and proving the completeness of the resulting stress functions.

The Gist

Imagine you are trying to predict how a piece of rubber, a block of metal, or even a tiny grain of sand will squish, stretch, or bend when you push on it.

For over a century, scientists have had a perfect "recipe book" for this. It's called Classical Elasticity. Think of it like a standard map. If you know the shape of the object and where you are pushing, this map tells you exactly how every point inside moves. It works great for big things like bridges, cars, and buildings.

But, in the last few decades, we've started building things that are incredibly small—like tiny sensors, nanomaterials, or microscopic cracks in a bridge. When things get this small, the standard map stops working. It's like trying to use a map of a whole country to navigate a single city street; you miss all the tiny details.

This is where Strain Gradient Elasticity (SGE) comes in. It's an upgraded, high-definition map that accounts for the "texture" of the material at a microscopic level. It realizes that when you push a tiny object, the change in how it stretches (the gradient) matters just as much as the stretch itself.

The Problem: Too Many Maps, No Compass

The problem with this new, high-definition map (SGE) is that it's complicated. Over the last 60 years, different scientists have tried to write down the "recipes" (mathematical formulas) to solve these problems. They came up with about ten different ways to write the solution, named after the people who invented them (like Mindlin, Papkovich, Neuber, Boussinesq, etc.).

It was like having ten different GPS apps for the same city. They all claimed to get you to the destination, but they spoke different languages, used different symbols, and no one knew for sure if they were all actually saying the same thing. Some were too messy to use, others were too simple to be accurate.

The Solution: The "Universal Translator"

This paper is like a team of expert cartographers who decided to clean up the mess. The authors (Solyaev, Hamouda, and Sherbakov) did three main things:

1. They built a "Universal Translator"
They discovered that all these different, complicated recipes are actually just variations of the old, simple recipes we already knew, plus a little "extra sauce."

• The Analogy: Imagine the old recipes are a basic soup (Classical Elasticity). The new SGE recipes are that same soup, but with a special "gradient spice" added to make it taste right for tiny ingredients.
• The Breakthrough: They showed that you don't need to learn ten new languages. You can take any old, trusted recipe, add this specific "gradient spice" (using a mathematical tool called a Helmholtz decomposition), and you instantly get the correct answer for the tiny world.

2. They simplified the "Mindlin" recipe
One of the most famous recipes, created by a scientist named Mindlin in 1964, was incredibly complex. It required calculating derivatives (math slopes) up to the fifth order—basically, doing the math equivalent of climbing a mountain five times just to get a cup of coffee.
The authors showed how to rewrite Mindlin's recipe so it looks almost exactly like the simple, classic recipes we've used for 100 years. They proved that the "Mindlin" way and the "Papkovich-Neuber" way are actually the same thing, just dressed in different clothes. This connects two huge families of solutions that had been living in parallel for decades without realizing they were related.

3. They proved the map is complete
In math, you have to prove that your recipe book covers every possible scenario. If you leave a hole, your map is useless. The authors proved that their new, simplified way of writing these solutions is complete. It means no matter what weird shape or force you throw at a tiny material, this method can find the answer.

Why Does This Matter?

You might ask, "Who cares about tiny math recipes?"

• For Engineers: It helps them design better micro-chips, stronger nanomaterials, and safer medical implants.
• For Scientists: It stops them from wasting time trying to invent new, complicated formulas when a simple one already exists.
• For the Future: As we move toward the "nanoworld" (building things atom by atom), having a clear, unified way to predict how materials behave is essential.

The Bottom Line

Think of this paper as the unification of a chaotic library. The authors took a shelf full of confusing, overlapping, and overly complicated books on how tiny materials move, and they wrote a new index card.

That card says: "Don't panic. You don't need to learn ten new languages. Just take the old, simple language you already know, add a tiny bit of 'gradient spice,' and you can solve any problem in the microscopic world."

They didn't just find new answers; they showed us that the answers were hiding in plain sight all along, waiting for someone to connect the dots.

Technical Summary

1. Problem Statement

The paper addresses the theoretical challenge of constructing general solutions for displacement fields in the static equilibrium equations of isotropic Strain Gradient Elasticity (SGE).

• Context: Classical elasticity theory possesses well-established general solutions (e.g., Boussinesq-Galerkin, Papkovich-Neuber, Love, Lamé) that allow for the derivation of analytical solutions for various boundary value problems.
• Gap: While SGE (specifically the Mindlin-Toupin theory) introduces higher-order equilibrium equations containing two length scale parameters ($l_1, l_2$) to capture size effects, a unified and comprehensive framework for general solutions in SGE has been fragmented. Existing SGE solutions (e.g., by Mindlin, Lurie et al., Charalambopoulos et al.) often differ in complexity, the number of potentials used, and their explicit separation of classical and gradient effects.
• Objective: To review existing SGE general solutions, derive new representations, establish rigorous mathematical relationships between them, and prove their completeness.

2. Methodology

The authors employ a rigorous mathematical approach based on operator formalism and Helmholtz decomposition:

• Governing Equations: The study focuses on the equilibrium equation for isotropic SGE:
  $$\alpha(1 - l_1^2\nabla^2)\nabla\nabla\cdot \mathbf{u} - (1 - l_2^2\nabla^2)\nabla\times\nabla\times \mathbf{u} = -\frac{\mathbf{b}}{\mu}$$
  where $\mathbf{u}$ is displacement, $\mathbf{b}$ is body force, and $l_1, l_2$ are length scale parameters.
• Decomposition Strategy: The core methodology involves decomposing the displacement field into a classical part (satisfying classical equilibrium) and a gradient part (satisfying modified Helmholtz equations).
• Green's Functions: The authors utilize Green's functions for the Poisson equation ($N$) and the inhomogeneous modified Helmholtz equation ($H$) to construct particular solutions.
• Comparative Analysis: The paper systematically compares ten different solution representations, deriving transformation rules (stress function interrelations) to prove that they are mathematically equivalent.

3. Key Contributions

A. Review and Derivation of Ten General Solutions

The paper catalogs and analyzes ten distinct representations of general solutions in SGE:

1. Mindlin Solution: The original generalization of Papkovich-Neuber, noted for its complexity (requiring up to 5th-order derivatives).
2. Lurie-Belov-Volkov-Bogorodskiy (LBVB): Separates classical and gradient parts but uses redundant potentials.
3. Charalambopoulos-Gortsas-Polyzos (CGP): A simplification for the case $l_1 = l_2$.
4. Generalized Boussinesq-Galerkin (BG): Uses a single vector function satisfying a bi-Helmholtz equation; useful for proving completeness.
5. Generalized Papkovich-Neuber (PN): A major focus. The authors present a simplified form that explicitly separates classical and gradient parts using a minimal number of potentials (3 independent functions for the gradient part).
6. Generalized Ter-Mkrtychan-Naghdi-Hsu (TNH): New derivation. The authors derive two new variants of the TNH solution for SGE, utilizing integral representations with specific Green's functions.
7. Specialized Solutions: Derivations for Generalized Lamé, Love, and Boussinesq potentials are provided as specific cases of the general solutions, showing how they reduce to classical forms when $l_1, l_2 \to 0$.

B. Universal Representation (The "Universal" Decomposition)

The authors propose a universal representation (Eq. 49-51) that unifies all approaches:
$$\mathbf{u} = \mathbf{u}_c + \mathbf{u}_g - \frac{l_1^2}{\alpha\mu}\nabla\phi_b$$

• $\mathbf{u}_c$: Any classical general solution (Papkovich-Neuber, Boussinesq, etc.).
• $\mathbf{u}_g$: A gradient part defined by Helmholtz decomposition (solenoidal and irrotational parts satisfying modified Helmholtz equations).
• $\phi_b$: A correction term related to the body force potential.
  This result demonstrates that any classical elasticity solution can be directly generalized to SGE by adding a specific gradient correction.

C. Interrelations and Equivalence

The paper establishes explicit mathematical transformations between the stress functions (potentials) of different solution types:

• Mindlin $\leftrightarrow$ Generalized PN: The authors derive the exact mapping between Mindlin's complex potentials and the simplified PN potentials, resolving a long-standing lack of connection between these two major formulations.
• BG $\leftrightarrow$ TNH: A non-trivial relationship is derived between the Boussinesq-Galerkin and TNH potentials, involving arbitrary functions that satisfy fourth-order equations.
• Completeness Proof: By leveraging the known completeness of the Boussinesq-Galerkin solution and the established interrelations, the authors prove that all derived general solutions (including Mindlin's and the new TNH forms) are complete. This means they can represent any displacement field satisfying the SGE equilibrium equations.

4. Key Results

• Simplification: The Generalized Papkovich-Neuber solution is identified as the most efficient form, requiring only first-order derivatives of potentials and a minimal set of independent functions.
• Redundancy Elimination: The paper clarifies that the LBVB solution contains redundant potentials for the gradient part, which can be reduced to the form used in the Generalized PN solution.
• Size Effect Regularization: The solutions naturally incorporate size effects, regularizing singularities (e.g., in Kelvin or crack problems) that appear in classical elasticity, consistent with previous numerical and analytical findings.
• Historical Continuity: The work confirms that SGE does not require entirely new mathematical structures but rather a systematic generalization of classical elasticity operators via Helmholtz decomposition.

5. Significance

• Theoretical Unification: The paper provides the first comprehensive "dictionary" translating between the various general solutions proposed over the last 60 years of SGE research. It proves that Mindlin's, Lurie's, and Charalambopoulos's solutions are mathematically equivalent.
• Analytical Toolset: By providing the Universal Representation, the authors enable researchers to solve complex SGE boundary value problems using familiar classical solutions (like Boussinesq or Love) and simply adding the gradient correction. This significantly lowers the barrier to entry for analytical modeling in strain gradient elasticity.
• Validation and Application: These general solutions are critical for:
  • Validating numerical methods (FEM, BEM) in SGE against closed-form benchmarks.
  • Analyzing size-dependent phenomena in micromechanics, nanomaterials, and fracture mechanics.
  • Deriving asymptotic solutions for cracks and singularities where classical theories fail.
• Completeness: The rigorous proof of completeness for Mindlin's solution and the new TNH variants fills a gap in the mathematical foundation of SGE, ensuring that no physical solution is missed by these formulations.

In summary, this paper serves as a definitive reference for the mathematical structure of displacement fields in strain gradient elasticity, unifying disparate historical approaches into a coherent, complete, and computationally efficient framework.