[論文レビュー] 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 Mechanics参考文献 4被引用数 3

ひとこと要約

この論文は、硬いインクルージョンを有する密に詰められた周期的複合材料の有効弾性率に対するFlaherty-Kellerの公式について、プライマル・デュアル変分原理を用いて、厳密な数学的証明を提供する。主な貢献は、インクルージョン間の狭いギャップにおける応力集中を捉える特異なひずみ関数の構築であり、これにより、間隔 $\epsilon$ に対する逆平方根依存性を示す漸近公式が得られ、一致漸近解析とエネルギー推定を用いて正確な定数が導出される。

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.

研究の動機と目的

• 密に詰められた周期的複合材料に硬いインクルージョンを有する場合の、有効弾性率に対するFlaherty-Kellerの漸近公式の数学的厳密な証明を提供すること。
• これまでヒューリスティックかつ数値的に導出されてきたFlaherty-Keller公式の、厳密な正当化の空白を埋めること。
• インクルージョン間隔 $\epsilon$ が 0 に近づく際の有効モジュラスの鋭い漸近的挙動を確立すること。
• 特異なひずみ関数が、硬いインクルージョン間の狭い領域における応力集中をどのように捉えるかを示すこと。

提案手法

• 有効エネルギーの上界と下界を導出するためのプライマル・デュアル変分原理の適用。
• 狭いギャップ近傍における漸近展開に基づく、Keller型のテスト関数をプライマル原理に構築すること。
• デュアル原理のテスト関数として、ひずみ核（Lamé系の解）から導出された特異関数の使用。
• Kelvin行列とLamé系の基本解を用いて、インクルージョン付近で正確に発散するひずみ核を生成すること。
• 境界プロファイル $f(y)$ のテイラー展開を用いた漸近解析により、$\epsilon \to 0$ の極限におけるエネルギー積分を推定し、$O(\epsilon^{-1/2})$ スケーリングを導出すること。
• 発散定理と境界積分の制御を用いたエネルギー推定により、誤差項を抑え、$O(1)$ の補正項を確立すること。

実験結果

リサーチクエスチョン

• RQ1 インクルージョン間隔 $\epsilon$ が 0 に近づく際、密に詰められた硬いインクルージョンを有する周期的複合材料における有効弾性率の正確な漸近的挙動は何か？
• RQ2 これまでヒューリスティックに導出されてきたFlaherty-Keller公式は、変分法を用いてどのように厳密に正当化できるか？
• RQ3 Lamé系の解である特異なひずみ関数は、硬いインクルージョン間の狭いギャップにおける応力集中をどのように捉えるか？
• RQ4 極端な材料対比と特異な応力挙動を示す弾性問題に、プライマル・デュアル変分原理を効果的に適用できるか？
• RQ5 有効モジュラスは、インクルージョンの幾何学的形状および最も狭い点における曲率 $\kappa_0$ に対して、どのように依存するか？

主な発見

• 有効縦方向モジュラス $\mathcal{E}_1$ は、$\epsilon \to 0$ のとき $\mathcal{E}_1 = (\lambda + 2\mu)\frac{\pi}{\sqrt{\kappa_0}} \frac{1}{\sqrt{\epsilon}} + O(1)$ を満たし、係数はLamé定数と曲率に依存する。
• 有効横方向モジュラス $\mathcal{E}_2$ は、$\mathcal{E}_2 = \mu\frac{\pi}{\sqrt{\kappa_0}} \frac{1}{\sqrt{\epsilon}} + O(1)$ を満たし、導電率の場合と同様に $\epsilon^{-1/2}$ スケーリングを示す。
• デュアル変分原理に用いられた特異関数は、ひずみ核の線形結合として構築され、Lamé系を正確に満たしており、応力集中の正確な制御が可能になる。
• プライマル原理からの上界とデュアル原理からの下界が、両方とも同じ $\epsilon^{-1/2}$ 渐近的次数を示しており、公式の鋭さが確認される。
• 誤差項は $\epsilon \to 0$ のとき一様に有界であり、$O(1)$ の補正項を有するため、漸近公式が極限においても有効であることが保証される。
• 本証明により、弾性理論におけるFlaherty-Keller公式の最初の厳密な導出が確立され、複合材料理論における長年の未解決問題が解決された。