A new perspective in linear Cauchy Elasticity: variational minimum principles for statics, dynamics, and heterogeneous materials

技術要約：線形コーシー弾性学における新たな視点：静力学、動力学、および不均質材料のための変分最小原理

問題の定義
本論文は、線形コーシー弾性の変分定式化における根本的な限界、特に不均質材料や蓄積エネルギー関数（すなわち、弾性係数 $C$ に主要対称性を欠くもの）を持たないシステムに関する課題に対処している。ハミルトンの原理のような標準的な変分原理は、弾性動力学の初期値問題（IVP）に適用する場合、重大な課題に直面する。これらの課題には、初期条件を自然に取り込むことの困難さや、変位に関する非因果的な最終時刻境界条件を要求することから生じる因果律の問題が含まれる。さらに、不定の弾性係数を持つ材料や非凸な材料の場合、原始的な系（primal system）は不良設定（ill-posed）であったり、一意の解を欠いたりすることがあり、標準的な凸変分アプローチを適用できなくなる。

手法
著者は、原始的な方程式が変分構造を持たない場合であっても、線形弾性動力学の方程式に関連付けられる一連の変分最小（凸）原理を構築する手法を開発した。その核心となる手法は以下の通りである：

1. 双対場変換（Dual Field Transformation）: 原始変数（変位 $u$、速度 $v$、対称歪 $E$、および反対称成分 $W$）を、変数変換を通じて一連の「双対」場（$\gamma, \lambda, \mu$）へと写像する。
2. 前双対汎関数（Pre-Dual Functional）の構成: 補助的なポテンシャル $H(U, \bar{U})$ を導入する。ここで $\bar{U}$ は任意の「基底状態（base state）」を表す。前双対汎関数 $\hat{S}[U, \gamma, \lambda, \mu]$ は、双対場に対してアフィンである。
3. 双対から原始への写像（DtP Map）: 原始変数 $U$ に関して前双対汎関数を最大化することにより、原始場を双対場と基底状態を用いて表す DtP 写像が導出される。
4. 凸双対汎関数: DtP 写像を前双対汎関数に代入することで、双対場のみに依存する双関数 $S[D]$ が得られる。この汎関数は、補助パラメータ $r > 0$ およびテンソル $A$ が正定値であるという条件下において、原始的な弾性係数 $C$ や密度 $\rho$ の性質に関わらず、凸（具体的には最小原理）であることが証明されている。
5. 時空領域定式化: 本手法は、ラプラス変換、周波数領域の仮定、または変数の二重化（Schwinger-Keldysh形式におけるもの）を必要とせず、実の時空領域で直接作用する。

主要な結果および理論的知見

• 最小原理の存在: 本論文は、線形弾性動力学の支配方程式（初期条件および境界条件を含む）が、凸双対汎関数 $S[D]$ のオイラー＝ラグランジュ方程式であることを確立している。
• 退化型楕円性（Degenerate Ellipticity）: 中心的な知見は、双対系のオイラー＝ラグランジュ系が、原始的な弾性動力学系が双曲型である場合でも、時空における退化型楕円系となることである。双対系の最高次微分に関連する二次形式は、正準半定値である。
• 一意性: 動力学的なケースにおいて、$C$ が正定値であり主要対称性を持つ場合、双対解の形式的な一意性が示されている。静的なケースでは、$C$ が主要対称性を欠いていても正定値であれば一意性が保持される。
• 不定の弾性係数の取り扱い: 双対系の退化型楕円性により、原始的な問題が不定の弾性係数（原始系が不良設定となる場合）を含む場合でも、最小原理の定式化が可能となる。補助ポテンシャル $H$ は選択基準として機能し、原始的な問題が非一意な解を持つ場合に、緩やかな不連続を伴う「良好な」解を回収できる可能性がある。
• 因果律: 双対場に対して最終時刻のディリクレ境界条件を課すことは、復元された原始解の因果律を損なわないと本論文は主張している。DtP 写像は双対場の時間微分を含むため、系は初期条件を満たすように最終時刻でのこれらの微分を調整することができる。
• 不均質材料: 材料の $C(x)$ や $\rho(x)$ が空間的に変化する場合でも、双対系は左辺に定数係数の最高次演算子（「比較媒体」）を持つように再定式化できる。これは、最高次演算子が解に依存する既存の定式化（分極応力法など）とは対照的である。この構造は、Moulinec-Suquet法への適応などの反復数値スキームにおいて利点を示唆している。

意義および主張
本論文は、本研究を新しい物理法則の発見ではなく、困難な偏微分方程式系を解くための数学的デバイスとして位置づけている。双対場は直接的な物理的意味を欠いていること（例：非因果的な境界条件を許容すること）が認められているが、このスキームは以下の理論的および計算的な利点を提供する：

1. 非標準的な問題への堅牢性: 本手法は、不定の弾性係数、非一意の解、あるいは「奇妙な弾性（Odd Elasticity）」を伴う弾性動力学など、原始的な系が不良設定である問題において特に価値がある。
2. 計算上の有用性: 双対汎関数の凸性は、「偽の時間（fake time）」変数を用いた勾配降下アル法による解の探索を可能にする。著者は、長時間の挙動を膨大な計算コストなしに扱うために、時空スライシングを用いた実装が可能であると述べている。
3. 均質化のポテンシャル: 不均質問題を、定数係数演算子を持つ時空境界値問題へと変換できることは、数学的な均質化スキームに対して利点をもたらすと推測される。
4. 限界: 著者は、標準的な線形弾性問題に対する実用的な効率性については控えめな姿勢をとっている。この定式化は、自由度を増加させ（原始の6から少なくとも12へ）、汎関数が厳密な強コアシビティ（coercivity）を保証するわけではないため、変分法の「大きなハンマー」（一般化ワイエルシュトラスの定理による最小化子の存在）が必ずしも適用できるとは限らない。したがって、勾配流スキームが弱解を得るための最も実用的な戦略として提示されている。

要約すると、本論文は、双対場を用いることで線形コーシー弾性動力学を凸変分問題へと変換する厳密な枠組みを提供しており、不均質材料、不定の弾性係数、および伝統的に変分的に扱うことが困難であった初期値問題を分析・解決するための新たなツールを提供している。