[Paper Review] A proof of the Flaherty-Keller formula on the effective property of densely packed elastic composites

Hyeonbae Kang, Sanghyeon Yu|arXiv (Cornell University)|Jul 7, 2017

Composite Material Mechanics4 references3 citations

TL;DR

This paper provides a rigorous mathematical proof of the Flaherty-Keller formula for the effective elastic moduli of densely packed periodic composites with hard inclusions, using a primal-dual variational principle. The key contribution is the construction of singular strain functions that capture stress concentration in the narrow gap between inclusions, yielding asymptotic formulas showing inverse square-root dependence on the inter-inclusion distance $\epsilon$, with precise constants derived via matched asymptotic analysis and energy estimates.

ABSTRACT

We prove in a mathematically rigorous way the asymptotic formula of Flaherty and Keller on the effective property of densely packed periodic elastic composites with hard inclusions. The proof is based on the primal-dual variational principle, where the upper bound is derived by using the Keller-type test functions and the lower bound by singular functions made of nuclei of strain. Singular functions are solutions of the Lamé system and capture precisely singular behavior of the stress in the narrow region between two adjacent hard inclusions.

Motivation & Objective

• To provide a mathematically rigorous proof of the Flaherty-Keller asymptotic formula for effective elastic moduli in densely packed periodic composites with hard inclusions.
• To address the long-standing gap in rigorous justification of the Flaherty-Keller formula, which had previously been derived heuristically and numerically.
• To establish sharp asymptotic behavior of the effective moduli as the distance $\epsilon$ between inclusions tends to zero.
• To demonstrate the role of singular strain functions in capturing stress concentration in narrow inter-inclusion regions.

Proposed method

• Application of the primal-dual variational principle to derive upper and lower bounds on the effective energy.
• Construction of Keller-type test functions for the primal principle, based on asymptotic expansions near the narrow gap.
• Use of singular functions derived from nuclei of strain—solutions to the Lamé system—as test functions for the dual principle.
• Employment of the Kelvin matrix and fundamental solutions to the Lamé system to generate strain nuclei with precise blow-up behavior near inclusions.
• Asymptotic analysis using Taylor expansion of the boundary profile $f(y)$ to estimate energy integrals in the limit $\epsilon \to 0$, yielding $O(\epsilon^{-1/2})$ scaling.
• Energy estimates via divergence theorem and boundary integral control to bound error terms and establish $O(1)$ corrections.

Experimental results

Research questions

• RQ1 What is the precise asymptotic behavior of the effective elastic moduli in a periodic composite with densely packed hard inclusions as the inter-inclusion distance $\epsilon$ tends to zero?
• RQ2 How can the Flaherty-Keller formula, previously derived heuristically, be rigorously justified using variational methods?
• RQ3 What role do singular strain functions—solutions to the Lamé system—play in capturing stress concentration in narrow gaps between hard inclusions?
• RQ4 Can the primal-dual variational principle be effectively applied to elasticity problems with extreme material contrasts and singular stress behavior?
• RQ5 What is the exact dependence of the effective moduli on the geometry of the inclusions and the curvature $\kappa_0$ at the narrowest point?

Key findings

• The effective longitudinal modulus $\mathcal{E}_1$ satisfies $\mathcal{E}_1 = (\lambda + 2\mu)\frac{\pi}{\sqrt{\kappa_0}} \frac{1}{\sqrt{\epsilon}} + O(1)$ as $\epsilon \to 0$, with the prefactor depending on Lamé constants and curvature.
• The effective transverse modulus $\mathcal{E}_2$ satisfies $\mathcal{E}_2 = \mu\frac{\pi}{\sqrt{\kappa_0}} \frac{1}{\sqrt{\epsilon}} + O(1)$, showing the same $\epsilon^{-1/2}$ scaling as in the conductivity case.
• The singular functions used in the dual variational principle are constructed as linear combinations of strain nuclei and solve the Lamé system exactly, enabling precise control of stress concentration.
• The upper bound from the primal principle and the lower bound from the dual principle both yield the same $\epsilon^{-1/2}$ asymptotic order, confirming the sharpness of the formula.
• The error terms are uniformly bounded as $\epsilon \to 0$, with $O(1)$ corrections, ensuring the asymptotic formula is valid in the limit.
• The proof establishes the first rigorous derivation of the Flaherty-Keller formula in elasticity, resolving a long-standing open problem in composite materials theory.