A Geometric Theory of Surface Elasticity and Anelasticity

In this paper we formulate a geometric theory of elasticity and anelasticity for bodies containing material surfaces with their own elastic energies and distributed surface eigenstrains. Bulk elasticity is written in the language of Riemannian geometry, and the framework is extended to material surfaces by using the differential geometry of hypersurfaces in Riemannian manifolds. Within this setting, surface kinematics, surface strain measures, surface material metric, and the induced second fundamental form follow naturally from the embedding of the material surface in the material manifold. The classical theory of surface elasticity of Gurtin and Murdoch (1975) is revisited and reformulated in this geometric framework, and then extended to anelastic bodies with anelastic material surfaces. Constitutive equations for isotropic and anisotropic material surfaces are formulated systematically, and bulk and surface anelasticity are introduced by replacing the elastic metrics with their anelastic counterparts. The balance laws are derived variationally using the Lagrange-d’Alembert principle. These include the bulk balance of linear momentum together with the surface balance of linear momentum, whose normal component gives a generalized Laplace’s law. As an application, we obtain the complete solution for a spherical incompressible isotropic solid ball containing a cavity filled with a compressible hyperelastic fluid, where the cavity boundary is an anelastic material surface with distributed surface eigenstrains. The analytical and numerical results quantify the effects of surface and fluid eigenstrains on the pressure-stretch response and residual stress.

💡 Research Summary

The paper presents a unified geometric framework for describing elasticity and anelasticity in bodies that contain material surfaces possessing their own elastic energies and distributed eigenstrains. By formulating bulk elasticity in the language of Riemannian geometry and extending this language to material surfaces through the differential geometry of hypersurfaces, the authors obtain a natural, coordinate‑free description of surface kinematics, strain measures, material metrics, and the induced second fundamental form.

The classical Gurtin‑Murdoch surface elasticity theory (1975) is revisited and recast in this geometric setting. The bulk deformation is described by a smooth map χ from a reference configuration to the current configuration, with the deformation gradient F and the bulk material metric C = Fᵀ G F. A material surface Σ is treated as an embedded 2‑dimensional hypersurface of the 3‑dimensional material manifold. Its surface deformation gradient 𝔽, surface material metric g = 𝔽ᵀ G 𝔽, and second fundamental form b arise directly from the embedding.

Constitutive relations are built from strain‑energy densities ψ_bulk(C) for the bulk and ψ_surf(g,b) for the surface. By applying the Lagrange‑d’Alembert principle, the authors derive the bulk balance of linear momentum together with a surface balance of linear momentum. The normal component of the surface balance reduces to a generalized Laplace law that couples surface tension, curvature, and surface eigenstrains.

A major contribution is the systematic inclusion of anelastic effects. Both bulk and surface are endowed with a reference material metric (G₀, g₀) and a current metric (G, g). The difference defines anelastic strain tensors (eigenstrains) that generate residual stresses τ = ∂ψ/∂E. Thus, bulk and surface residual stresses emerge naturally from metric incompatibility, without the need for ad‑hoc incompatibility tensors. The framework accommodates isotropic and anisotropic surface constitutive laws, and it can be specialized to purely membrane‑type surfaces (no curvature dependence) as done in the present work.

The theory is illustrated with a fully nonlinear problem: an incompressible isotropic solid sphere of radius R containing a concentric spherical cavity of radius a. The cavity boundary is a material surface endowed with its own anelastic metric (surface eigenstrain ε_s) and isotropic surface elastic constants (μ_s, λ_s). The cavity is filled with a compressible hyperelastic fluid characterized by a bulk strain‑energy ψ_f(V) and a natural volume ratio η_f. By exploiting spherical symmetry, the authors reduce the governing equations to a set of ordinary differential equations for the radial stretch λ(r). Boundary conditions include an external pressure applied at r = R and continuity of traction across the cavity surface. Analytic solutions are obtained for the small‑deformation limit, while a Newton‑Raphson scheme provides numerical solutions for large deformations.

Results show that positive surface eigenstrains (surface expansion) stiffen the overall response, raising the pressure required for a given stretch, whereas negative eigenstrains (surface contraction) soften the response. The fluid’s natural volume ratio also influences the pressure–stretch curve: a fluid that prefers a larger volume than the cavity imposes additional internal pressure, leading to non‑uniform residual stresses throughout the solid. The generalized Laplace law captures the combined effect of surface tension, curvature, and surface eigenstrain on the normal stress balance.

In conclusion, the authors have built a mathematically rigorous, coordinate‑free theory that treats bulk and surface elasticity and anelasticity on equal footing. By identifying the material metric as the fundamental descriptor of both elastic and anelastic deformations, the framework reproduces classical results while extending them to situations with residual stresses and eigenstrains on surfaces. The paper opens pathways for future work on surface bending effects, visco‑anelastic surface behavior, and multi‑cavity or multi‑surface composites, where the interplay of bulk and surface metrics will be essential for accurate mechanical predictions.