On singular behaviour in a plane linear elastostatics problem

Basic Information

• Paper ID: 2409.07954
• Title: On singular behaviour in a plane linear elastostatics problem
• Authors: H. Gimperlein, M. Grinfeld, R. J. Knops, M. Slemrod
• Classification: math.AP (Mathematics - Partial Differential Equations)
• Publication Date: November 25, 2024 (v2)
• Paper Link: https://arxiv.org/abs/2409.07954

Abstract

This paper constructs a vector field analogous to those introduced separately by Artstein and Dafermos, derived from tangent vectors to a monotonically increasing single-parameter family of non-concentric circles tangent at the origin. These circles define and fill a lens-shaped region Ω, whose outer and inner boundaries are the largest and smallest circles, respectively. The cusp singularity at the origin generates geometric singularity, at which the vector field is indeterminate and exhibits non-unique limiting behavior. Through the semi-inverse method involving Airy stress functions, it is shown that this vector field corresponds to the displacement vector field of a linear plane compressible non-homogeneous isotropic elastostatics equilibrium problem on Ω with relatively rigid rotational boundary conditions, which may induce fracture or tearing at the origin. The paper finds a family of solutions where not only do the Lamé parameters satisfy the strong ellipticity condition, but the non-unique limiting behavior of displacements is also preserved.

Research Background and Motivation

Problem Formulation

1. Core Problem: Investigate the singular behavior of displacement fields in plane elastostatics problems with geometric singularities. Specifically, explore the indeterminacy and non-unique limiting behavior of the displacement field at the singular point when the boundary undergoes relative rigid rotation on a lens-shaped region (with a cusp at the origin).
2. Problem Significance:
   • From a mathematical perspective, this problem involves existence, uniqueness, and regularity of solutions to partial differential equations on domains with geometric singularities
   • From a physical perspective, singular behavior at singular points may correspond to important phenomena such as material fracture, tearing, or brittle failure
   • This research extends the work of Artstein and Dafermos on non-unique solutions, providing new interpretations within the framework of continuum mechanics
3. Limitations of Existing Methods:
   • Artstein considered a variant of this vector field in control systems, but no smooth stabilizing control exists
   • Dafermos used it as an example of non-unique non-smooth solutions in evolutionary ordinary differential equations
   • Prior work (reference 3) proved that this vector field constitutes Lagrangian trajectories of fluid particles in stable compressible flow
   • However, systematic investigation of this vector field's singularity within the elastostatics framework was lacking
4. Paper Motivation:
   • Employ a method different from Artstein and Dafermos to construct the vector field (based on tangents to circle families)
   • Interpret singular behavior within the linear plane elastostatics context
   • Find non-homogeneous material parameters satisfying the strong ellipticity condition
   • Deepen understanding of singular point properties through computation of total force, total torque, and strain energy

Core Contributions

1. Geometric Construction of Vector Field: Proposes a vector field construction method based on tangent lines to a single-parameter family of non-concentric circles, which is naturally generalizable and clearly exhibits the source of geometric singularity
2. Elastostatics Interpretation: Through the semi-inverse method and Airy stress functions, proves that this vector field corresponds to the displacement field of a non-homogeneous compressible isotropic linear elastic body in equilibrium under zero body forces
3. Strong Ellipticity Condition: Finds a family of solutions whose Lamé parameters satisfy the strong ellipticity condition $\mu(\lambda + 2\mu) > 0$ while preserving the non-unique limiting behavior of displacements
4. Singularity Analysis:
   • Proves the indeterminacy of the vector field at the origin and its non-unique limiting behavior (Proposition 4.1)
   • Computes total force (zero) and total torque (proportional to $R^2-1$) (Proposition 6.1)
   • Analyzes singularity of strain energy ($O(r^{-1})$ when $k>0$, non-singular when $k=0$)
5. Physical Interpretation: Interprets the relative rigid rotation of the boundary and singular behavior at the origin as possible fracture, tearing, or brittle failure phenomena

Methodology Details

Task Definition

Investigate the plane linear elastostatics boundary value problem on lens-shaped region $\Omega \subset \mathbb{R}^2$:

• Input: Geometric domain (defined by non-concentric circle family), boundary conditions (rigid rotation)
• Output: Displacement field $\mathbf{u}(x_1, x_2)$, stress field $\sigma_{\alpha\beta}$, non-homogeneous Lamé parameters $\lambda(x), \mu(x)$
• Constraints: Equilibrium equations $\sigma_{\alpha\beta,\beta} = 0$, constitutive relations, strong ellipticity condition

Geometry and Vector Field Construction

1. Circle Family Definition Consider a single-parameter family of non-concentric circles: $\psi_c(x) := c^{-2}\{(x_1-c)^2 + x_2^2\} = 1$ where $1 \leq c \leq R$, with centers at $(c,0)$, radius $c$, all circles tangent at the origin $(0,0)$.

2. Tangent Vector Field Unit tangent vector at point $x$ on the circle: $t_1(x_1,x_2) = \frac{x_2}{c}, \quad t_2(x_1,x_2) = \frac{c-x_1}{c}$

3. Displacement Vector Field Setting $\mathbf{u} = c\mathbf{t}$, obtain: $u_1(x_1,x_2) = x_2$$u_2(x_1,x_2) = \frac{x_2^2-x_1^2}{2x_1} = c-x_1$

4. Domain Definition

• Outer boundary: $\psi_R(x) = 1$ (radius $R$)
• Inner boundary: $\psi_1(x) = 1$ (radius $1$)
• Region $\Omega$ is lens-shaped with symmetric cusps at the origin

Vector Field Properties Analysis

1. Strain Components$e_{11} = 0, \quad e_{22} = \frac{x_2}{x_1}, \quad e_{12} = \frac{x_1^2-x_2^2}{4x_1^2}$

2. Polar Coordinate Representation In polar coordinates $(r,\theta)$, the circle equation is $r = 2c\cos\theta$, with radial and tangential components: $u_r(c,\phi) = 0, \quad u_\phi(c,\phi) = -c$ indicating the vector field is tangential to circles with constant magnitude $c$.

3. Rigid Motion Decomposition$\mathbf{u}(x_1,x_2) = c\mathbf{e}_2 + (x_2\mathbf{e}_1 - x_1\mathbf{e}_2)$ Expressed as: translation by $c$ in $x_2$ direction + unit rigid body rotation clockwise about the origin.

4. Singularity (Proposition 4.1)

• At the origin, $u_2$ given by $u_2 = (x_2^2-x_1^2)/(2x_1)$ is indeterminate
• When approaching the origin along circle $\psi_c=1$: $\lim_{x_1\to 0} u_2 = c$
• Limit value depends on which circle is approached, hence non-unique
• Jump across different circles: $[u_2] = c_1 - c_2$

Semi-Inverse Method and Airy Stress Function

1. Airy Stress Function The general solution to equilibrium equations $\sigma_{\alpha\beta,\beta}=0$ can be expressed using Airy function $\chi(x_1,x_2)$: $\sigma_{11} = -\chi_{,22}, \quad \sigma_{12} = \chi_{,12}, \quad \sigma_{22} = -\chi_{,11}$

2. Governing Equation Define parameter: $\Lambda(x_1,x_2) := \frac{e_{11}-e_{22}}{e_{12}} = \frac{\sigma_{11}-\sigma_{22}}{\sigma_{12}}$ obtaining a hyperbolic equation independent of Lamé parameters: $\chi_{,11} - \chi_{,22} - \Lambda(x_1,x_2)\chi_{,12} = 0$

3. Complex Variable Solution Introduce $z = x_1 + ix_2$, $\bar{z} = x_1 - ix_2$, transforming the equation to: $z^2\chi_{,zz} + \bar{z}^2\chi_{,\bar{z}\bar{z}} = 0$

Solving within the special class $z^2\chi_{,zz} = J(z,\bar{z}) = -\bar{z}^2\chi_{,\bar{z}\bar{z}}$, where: $J(z,\bar{z}) = z\bar{z} + k$ with $k$ a non-negative constant.

4. Airy Function Solution$\chi(z,\bar{z}) = (z\bar{z}-k)\log\frac{z}{\bar{z}}$ In polar coordinates: $\chi(r,\theta) = (r^2-k)2i\theta$

Stress, Strain and Material Parameters

Derived from the Airy function:

Stress Components: $\sigma_{11} = -2\left[2\theta + \frac{r^2+k}{r^2}\sin 2\theta\right]$$\sigma_{22} = 2\left[-2\theta + \frac{r^2+k}{r^2}\sin 2\theta\right]$$\sigma_{12} = 2\frac{r^2+k}{r^2}\cos 2\theta$

Lamé Parameters: $\mu = 4\frac{r^2+k}{r^2}\cos^2\theta$$\lambda = -4\theta\cot\theta - 4\frac{r^2+k}{r^2}\cos^2\theta$

Strong Ellipticity Condition: $\lambda + 2\mu = 4\cos\theta\left(\cos\theta - \frac{\theta}{\sin\theta}\right) + \frac{k}{c^2}$

When $k=0$, $\lambda+2\mu < 0$, violating strong ellipticity. When $k > 0$ and sufficiently large, one can ensure $\lambda+2\mu > 0$ and $\mu > 0$, satisfying strong ellipticity.

Technical Innovations

1. Systematic Geometric Construction: Constructing the vector field from circle family tangents provides greater geometric intuition and generalizability compared to Artstein and Dafermos's methods
2. Complex Variable Techniques: Cleverly employing complex variables to transform the hyperbolic PDE into a more tractable form
3. Parametric Control: Introducing parameter $k$ enables control of Lamé parameter strong ellipticity while maintaining displacement singularity
4. Non-homogeneous Materials: Proves that this problem can only be realized within the non-homogeneous material framework; homogeneous materials cannot satisfy the equilibrium equations

Experimental Setup

This paper is pure mathematical theoretical research without numerical experiments or experimental data, relying primarily on rigorous mathematical derivation and analysis.

Analysis Tools

1. Analytical Computation: Explicit expressions for strain, stress, and Lamé parameters
2. Limit Analysis: Investigation of vector field limiting behavior at the origin
3. Integral Computation: Calculation of total force, total torque, and strain energy
4. Asymptotic Analysis: Asymptotic behavior of physical quantities as $a\to 0$

Computational Cases

The paper provides detailed calculations of the following physical quantities:

1. Geometric Parameters (Table 1):

• Coordinates of intersection points between circle $\psi_c=1$ and ball $B(0,a)$ and corresponding trigonometric values
• Boundary parameterization for subsequent integral calculations

2. Boundary Traction: $F_1 = -4\theta\cos 2\theta$$F_2 = -4\theta\sin 2\theta + 2\left(1+\frac{k}{r^2}\right)$

Experimental Results

Main Results

1. Total Force (Proposition 6.1) Through integration over the punctured region $\Omega_a(c_2,c_1)$ and taking $a\to 0$:

• Total force in $x_1$ direction: $T_1 = 0$ (since $F_1$ is odd in $\theta$)
• Total force in $x_2$ direction: $T_2 = 0$

2. Total Torque (Proposition 6.1)$\Gamma = 4\pi(R^2-1)$ indicating that relative rotation between inner and outer boundaries produces non-zero torque, proportional to the difference of squared radii.

3. Strain Energy Over region $\Omega_a(1,R)$: $\lim_{a\to 0}|E_a(1,R) - \frac{2k}{a}(R-1)| = 0$

• When $k > 0$: strain energy $E_a = O(a^{-1})$, singular at the origin
• When $k = 0$: strain energy non-singular (but Lamé parameters violate strong ellipticity)

4. Strong Ellipticity

• $k=0$: $\mu > 0$ but $\lambda+2\mu < 0$, violating strong ellipticity
• $k > 0$ sufficiently large: $\mu > 0$ and $\lambda+2\mu > 0$, satisfying strong ellipticity

Important Findings

1. Preservation of Non-uniqueness: Even when strong ellipticity is satisfied, the non-unique limiting behavior of displacement at the origin persists, a fundamental characteristic of geometric singularity
2. Relationship Between Torque and Energy:
   • Zero total force indicates overall system equilibrium
   • Non-zero torque reflects relative boundary rotation
   • Existence/absence of strain energy singularity relates to material parameter strong ellipticity
3. Physical Interpretation of Singularity: Singular behavior at the origin may correspond to:
   • Material fracture or tearing
   • Brittle failure
   • Dislocations or other defects
4. Isotropy Along Circles: On each circle $\psi_c=1$, the vector field is isotropic ($e_{rr}+e_{\theta\theta}=0$)

Related Work

Vector Field Origins

1. Artstein (1983): Introduced a variant of this vector field in control systems, studying cases where no smooth stabilizing control exists
2. Dafermos (2012): Applied this vector field as an example of non-unique non-smooth solutions in evolutionary ordinary differential equations, proposing "entropy rate" criteria for selecting unique solutions
3. Authors' Prior Work (2024): Proved that this vector field constitutes Lagrangian trajectories of fluid particles in stable compressible flow, with potentially non-unique trajectories

Elasticity Theory Foundations

4. Classical Textbooks:
   • Gurtin (1972): Linear elasticity theory
   • Landau & Lifshitz (1970): Elasticity theory
   • Love (1927): Monograph on mathematical elasticity
5. Strong Ellipticity: Knops & Payne (1971) on uniqueness theorems in elasticity

Singularity Theory

6. Grisvard (1985), Kondratiev (1967): Boundary value problems for elliptic equations on non-smooth domains, providing theoretical foundation for understanding geometric singularities

Unique Contributions of This Paper

• Method Innovation: Geometric construction based on circle family tangents
• Theory Deepening: Systematic investigation of singularity within elastostatics framework
• Parametric Control: Achieving balance between strong ellipticity and singularity through parameter $k$
• Complete Analysis: Comprehensive calculations of force, torque, and energy

Conclusions and Discussion

Main Conclusions

1. Existence: Constructs a class of displacement fields for non-homogeneous compressible isotropic linear elastic bodies satisfying equilibrium equations on lens-shaped regions
2. Singularity Characteristics:
   • Displacement is indeterminate at the origin with non-unique limits
   • Limit values depend on approach path (which circle)
   • This is a direct consequence of geometric singularity (cusp)
3. Strong Ellipticity: By selecting appropriate parameter $k$, one can ensure Lamé parameters satisfy the strong ellipticity condition while maintaining singular behavior
4. Physical Interpretation:
   • Relative rigid boundary rotation produces non-zero torque
   • Fracture, tearing, or brittle failure may occur at the origin
   • Strain energy singularity depends on material parameters

Limitations

1. Two-Dimensional Restriction: Current analysis limited to plane problems; three-dimensional extension requires Maxwell-Morera functions
2. Special Geometry: Lens-shaped region and circle family are special choices; while methods generalize, specific calculations become complex
3. Linear Theory: Employs linear elasticity theory; large deformation cases not considered
4. Statics: Considers only static equilibrium; dynamic behavior not addressed
5. Theoretical Analysis: Lacks numerical verification or experimental comparison

Future Directions

The paper's Section 8 proposes multiple research directions:

1. Deeper Understanding of Singularity: Analogous to Laplace equation research on sectorial domains (Grisvard, Kondratiev), recover uniqueness through appropriate growth conditions
2. Geometric Generalization:
   • Consider circle families not tangent at origin ($d_2 > 0$)
   • Investigate other curve families (ellipses, parabolas, etc.)
3. Domain Extension:
   • Extend vector field to $\Omega(0,1)$ (rigid rotation)
   • Extend to exterior region $\Omega(R,\infty)$
4. Small Hole Problem: Replace cusp with small circle $c_2=\epsilon$, study perforated domains
5. Half-Plane Problem: Let $R\to\infty$, investigate singularity distribution on half-planes
6. General Vector Fields: Apply methods to arbitrary conservative plane vector fields
7. Alternative Stress Functions: Explore other solutions to equation (5.1)
8. Helmholtz Decomposition: Express vector field as $\mathbf{u}=\nabla\Psi+\nabla\times\boldsymbol{\psi}$, connect to dislocation theory
9. Three-Dimensional Extension: Generalize to three dimensions using Maxwell-Morera functions
10. Relaxed Control: Consider $\mathbf{u}=cw(x)\mathbf{t}$ where $w(x)$ is arbitrary scalar function

In-Depth Evaluation

Strengths

1. Mathematical Rigor

• All conclusions rigorously proved
• Detailed calculations with clear logic
• Careful treatment of limiting processes (e.g., asymptotic analysis as $a\to 0$)

2. Method Innovation

• Novel and intuitive geometric construction method
• Skillful application of complex variable techniques
• Effective combination of semi-inverse method with Airy stress functions

3. Theoretical Depth

• Reveals profound connections between geometric singularity and solution non-uniqueness
• Achieves balance between strong ellipticity and singularity through parameter $k$
• Complete mechanical analysis (force, torque, energy)

4. Physical Insight

• Connects mathematical singularity to physical phenomena (fracture, tearing)
• Rigid motion decomposition provides clear physical picture
• Proves isotropy properties

5. Generalizability

• Methods applicable to other curve families
• Airy function solutions generalizable (Appendix B)
• Provides framework for related problem research

Weaknesses

1. Lack of Numerical Verification

• No numerical simulations supporting theoretical results
• No stress and strain distribution plots
• Missing quantitative analysis of parameter $k$ effects

2. Absence of Physical Experiment Comparison

• No experimental data validating fracture/tearing interpretation
• Missing comparison with actual material behavior

3. Geometric Specificity

• Lens-shaped region and circle family are special choices
• Treatment of more general geometries insufficiently discussed

4. Limited Application Scenarios

• Linear theory restricts practical applications
• Two-dimensional problems have relatively limited practical significance

5. Uniqueness Issues

• While non-uniqueness is identified, no selection criteria proposed
• Application of Dafermos's "entropy rate" criterion to elasticity problems unexplored

6. Presentation Aspects

• Complex notation in places (multiple coordinate systems)
• Lacks illustrative diagrams (only two simple figures)
• Physical interpretations could be more thorough

Impact Potential

1. Theoretical Contribution

• Provides new perspective for non-homogeneous elastic singularity research
• Enriches understanding of non-unique solutions under strong ellipticity
• Connects geometric singularity, solution properties, and physical phenomena

2. Methodological Value

• Geometric construction method can inspire related problem research
• Demonstrates application of complex variables and semi-inverse methods
• Provides technical reference for handling geometric singularities

3. Potential Applications

• Fracture mechanics: stress-strain analysis at singular points
• Materials science: non-homogeneous material design
• Dislocation theory: connections to Volterra dislocations (mentioned in Section 8)

4. Limitations

• Limited direct practical application in near term
• Requires further numerical and experimental support
• Three-dimensional extension presents challenges

Applicable Scenarios

1. Theoretical Research

• PDE singularity theory
• Non-uniqueness problems in nonlinear analysis
• Effects of geometric singularity on solutions

2. Materials Science

• Stress analysis in non-homogeneous materials
• Composite material interface problems
• Functionally graded material design

3. Fracture Mechanics

• Crack tip stress singularity
• Tearing and brittle failure mechanisms
• Multi-crack interaction

4. Mathematical Physics

• Mathematical foundations of continuum mechanics
• Non-unique solutions of conservation laws
• Entropy conditions and selection criteria

5. Educational Use

• Semi-inverse method examples in elasticity
• Complex variable applications to PDEs
• Concrete cases of geometric singularity

References

The paper cites nine main references:

1. Artstein (1983): Relaxed control stabilization in control systems
2. Dafermos (2012): Maximum dissipation in evolution equations
3. Gimperlein et al. (2024): Non-uniqueness in plane fluid flow (authors' prior work)
4. Grisvard (1985): Elliptic equations on non-smooth domains
5. Gurtin (1972): Linear elasticity theory (Encyclopedia of Physics)
6. Knops & Payne (1971): Uniqueness theorems in elasticity
7. Kondratiev (1967): Boundary value problems on conical or singular domains
8. Landau & Lifshitz (1970): Elasticity theory
9. Love (1927): Monograph on mathematical elasticity

Overall Assessment

This is a high-quality mathematical physics paper with significant contributions in both theoretical depth and methodological innovation. Through rigorous mathematical analysis, the authors reveal profound connections between geometric singularity, solution non-uniqueness, and physical phenomena. While there is room for improvement in numerical verification and practical applications, as theoretical research, this work provides important mathematical tools and physical insights for related fields. Particularly noteworthy is that the paper not only solves the specific problem but also proposes systematic research methods and multiple valuable future research directions, demonstrating good potential for academic impact.