A family of variational principles of minima for the plasticity, the… — やさしい解説

あなたは、複雑なシステムが時間の経過とともにどのように振る舞うかを予測しようとしていると想像してください。例えば、熱によって曲がる金属製の梁、二つの粗い表面が擦れ合う様子、あるいはガラスの中に広がる亀裂などです。通常、科学者はこれらの問題をステップ・バイ・ステップで、つまり「今現在の位置」に基づいて「次の位置」を計算しながら、山を一歩ずつ登るように解こうとします。

この論文は、より「一括（all-at-once）」的な考え方を提案しています。ステップ・バイ・ステップで登るのではなく、始まりから終わりまでの「全行程」を一つの統一された経路として捉え、あらゆる可能な経路の中から「最善の」ものを見つけ出すという方法です。

以下に、簡単な比喩を用いて、この論文のアイデアを解説します。

1. 大きなアイデア：「映画」対「スナップショット」

ほとんどの工学計算は、一連のスナップショットを撮るようなものです。1秒時点の状態を計算し、次に2秒、その次に3秒……という具合です。
著者であるG. de Saxcéは、「映画」的なアプローチを提案しています。彼は**変分原理（Variational Principle）**を提示します。これは次のようなルールです。「このシステムの歴史を描いたあらゆる可能な映画の中で、自然界は特定の『コスト』を最小化する映画だけを選択する」。

もし、この「コスト」をゼロにする経路を見つけることができれば、それはシステムの真の物理的挙動を見つけたことになります。

2. ツールキット：二つの幾何学

この「映画」のルールを構築するために、著者は二種類の異なる幾何学を組み合わせています。

可逆的な部分（シンプレクティック幾何学 / Symplectic Geometry）: これは、摩擦のない振り子の揺れのような、物理学の「完璧な」部分を扱います。それは、エネルギーが保存される摩擦のないアイススケートリンクのようなものです。
不可逆的な部分（凸解析 / Convex Analysis）: これは、摩擦や塑性変形（金属が曲がったままになる現象）、あるいは亀裂のように、エネルギーが失われる「乱れた」部分を扱います。ここでは物事が「粘りけ」を持ったり「粗く」なったりします。

この論文の主なトリックは、これら二つを組み合わせることです。著者は、システムを「可逆的なエンジン（バネのようなもの）」と「散逸的なブレーキ（摩擦のようなもの）」を持つものとして扱い、これらが全タイムラインにわたって完璧にバランスを取る数学的公式を見つけ出します。

3. 「BEN」原理：完璧な経路を見つける

この論文の核心は、Brezis-Ekeland-Nayroles (BEN) と呼ばれる有名なアイデアの拡張です。

比喩: あなたが、重い砂袋（摩擦）を引きずりながら、地点Aから地点Bまでボールを転がす最もスムーズな経路を探していると想像してください。
論文の主張: あなたが想像したあらゆる経路の「粗さ」を計算する、特定の数学的公式（「汎関数」）が存在します。
- もし、自然界が取らないであろう経路を推測した場合、その公式は正の数（ペナルティ）を返します。
- もし、自然界が実際に辿る経路を推測した場合、その公式はゼロを返します。
- したがって、この問題を解くには、この公式をゼロにする経路を見つけるだけでよいのです。

4. これは何を解決するのか？

著者は、標準的な数学がしばしば苦戦する三つの困難な領域において、この「映画」的アプローチが有効であることを示しています。

塑性（金属の曲がり）: 紙クリップを曲げると、それは元の形に戻りません。この論文は、「ゼロ・コスト」のルールを用いることで、曲がるプロセス全体をステップ・バイ・ステップではなく、一度に計算する方法を示しています。
摩擦接触（表面の擦れ）: 二つの粗い表面が接触すると、それらは複雑に固着したり滑ったりします。論文では、これらを単純な「滑らかな」形状に無理やり当てはめることなく、**「バイポテンシャル（Bipotential）」**と呼ばれるツール（二面的なマップのようなもの）を使用して、この固着・滑りの挙動を記述しています。
破壊（ガラスの亀裂）: これは最も劇的な例です。亀裂が発生するとき、それは通常、特定の方向に飛び進みます。
- 問題: 古い手法では、計算が「明示的（explicit）」で小さな誤差に敏感すぎるため、亀裂が間違った方向に進むと予測してしまうことがよくありました。
- 論文の解決策: 「インプリシット（implicit/暗黙的）」な計算（ステップ全体を一度に見る方法）を用いたこの「映画」的アプローチを用いることで、著者のモデルは亀裂の経路をより正確に予測します。これは、亀裂が特定の角度で「キンク（折れ曲がり）」を生じるという実世界の実験結果と一致します。

5. 「シンプレクティック」のひねり

著者は「シンプレクティック（Symplectic）」という高度な用語を導入しています。

簡単な説明: 物理学において「シンプレクティック」とは、位置と運動量（速度と位置）に関する情報を一体として整理する方法です。
論文の貢献: 著者はこの「シンプレクティック」な整理法を、エネルギーが失われる系（散逸系）にも適用できる形に拡張しました。通常、シンプレクティック数学は、エネルギーが保存される完璧なシステムのためのものです。著者は、この強力な数学を、摩擦や亀裂のような、乱れた現実世界のシステムに適用するための架け橋を築きました。

6. 非標準的なルールのための「バイポテンシャル」

クーロン摩擦のような物理法則は、数学の標準的な「滑らかな」ルールに従いません。それらは「非関連（non-associated）」であり、つまり、動きの方向が押されている力と完全に一致しないことを意味します。

比喩: 重い箱を押している場面を想像してください。通常、押せば押した方向に動きます。しかし、摩擦がある場合、十分に強く押すまでは箱は固着し、その後、横方向に滑り出すかもしれません。
論文のツール: 著者は**「バイポテンシャル」**を使用しています。これは、このような奇妙で非滑らかなルールを扱うことができる、特別な「翻訳機」のようなものです。これにより、物理現象が単純な直線に従わず、乱れている場合でも、「映画」の原理が機能することを可能にしています。

まとめ

この論文は新しい物理法則を発明したのではなく、既存の法則を解くための新しい方法を発明したものです。
システムの未来を一秒ごとに計算する代わりに、システムの全履歴を一度に計算する方法を提案しています。それは、正しい経路に対してはゼロになるはずの「コスト関数」を使用しています。完璧な運動の幾何学（シンプレクティック）と、乱れた損失の幾何学（凸解析）を組み合わせることで、著者は金属の曲がり、表面の擦れ、そして亀裂の成長を正確に予測する統一された枠組みを作り上げました。これは、多くの場合、従来のステップ・バイ・ステップの手法よりも優れた性能を発揮します。

技術要約：塑性、摩擦接触、および破壊力学における最小値の変分原理の族

問題提起
本論文は、動的な散逸系に対して、統一された非増分的な時空変分原理を定式化するという課題に取り組んでいる。従来の変分アプローチは、非平滑な構成則（塑性、摩擦接触、破壊など）や、散逸ポテンシャルが凸でない、あるいはフロー則が降伏面に対して法線方向ではない非関連フロー則を扱う際に、しばしば困難に直面する。また、既存の散逸系の統一フレームワーク（GENERICや速度依存のない系など）は、散逸ポテンシャルの1次同次性といった制限的な仮定に依存していることが多く、これが粘塑性や一般的な動的ケースへの適用性を制限している。著者は、これらの複雑性を一貫した幾何学的設定の中で扱うために、Brezis-Ekeland-Nayroles (BEN) 原理をシンプレクティック・フレームワークへと拡張することを目指している。

手法
手法は、凸解析とシンプレクティック幾何学の統合に基づいている。コアとなるアプローチには、以下のステップが含まれる：

シンプレクティック分解： 状態ベクトル $z$ （自由度と運動量で構成される）の時間変化率を、可逆部分 $\dot{z}_R$ （ハミルトニアン $H$ のシンプレクティック勾配）と不可逆部分 $\dot{z}_I$ （散逸部分）に加法的に分解する。
シンプレクティック凸解析： 著者は、散逸ポテンシャルのシンプレクティック劣微分 $\partial_\omega \phi$ およびシンプレクティック・フェンシェル極 $\phi^*_\omega$ を導入する。これにより、標準的なフェンシェル不等式をシンプレクティックな文脈へと一般化し、散逸ポテンシャルが必ずしも1次同次ではない系のモデリングを可能にする。
SBEN原理： シンプレクティック・ブレシス・エケランド・ネイロールス（SBEN）原理を構築する。これは、散逸系の自然な進化経路が、許容されるすべての経路の中で特定の汎関数 $\Pi(z)$ を最小化するというものである。この汎関数は、散逸ポテンシャル、そのシンプレクティック極、およびシンプレクティック形式を組み合わせたものであり、真の物理的経路において最小値がゼロとなる。
非関連則への拡張： 非関連な構成則（クーロン摩擦など）を扱うために、著者はバイポテンシャル（双凸であり、特定の極値条件を満たす関数 $b(x,y)$ ）の概念を導入する。これはさらに、動的な非関連系を収容するためにシンプレクティック・バイポテンシャル ( $\hat{b}$ ) へと拡張される。
特定の力学への適用： フレームワークは以下に適用される：
- 標準的および粘塑性（Hillの不等式の回復）。
- シニョルニ（Signorini）の一面接触（線形相補問題への連結）。
- 有限歪塑性（加法的および乗法的分解の両方を使用）。
- 流体力学（オイラー記述におけるナビエ・ストークス流およびビンガム流）。
- 摩擦を伴う亀裂進展（亀裂の進展を、亀裂表面上のフローとしてモデリング）。

主要な貢献

統一されたシンプレクティック・フレームワーク： 本論文は、散逸ポテンシャルの制限的な1次同次性を必要とせずに、可逆的（ハミルトニアン）ダイナミクスと不可逆的（散逸的）ダイナミクスの扱いを統一する、一般的なSBEN原理を提案している。
シンプレクティック・バイポテンシャル： 非関連な構成則（これまで変分的に扱うことが困難であったクラスの問題）へ変分原理を拡張するための、新しい数学的ツールであるシンプレクティック・バイポテンシャルを定義した。
非増分定式化： これらの原理は「非増分的」である。つまり、ステップごとに解くのではなく、全時間区間 $[0, T]$ を同時に考慮し、進化経路に対するグローバルな視点を提供する。
破壊力学への適用： 本論文は、脆性破壊への適用として、亀裂の進展をフロー場を用いてモデリングし、結果として得られる変分問題を解くための暗黙的スキームを利用している。

結果

塑性と接触： SBEN原理は、塑性およびシニョルニ接触に関する古典的な構成則を、極値条件として正常に回復する。接触については、この定式化は線形相補問題（LCP）へと帰着する。
流体力学： この原理はナビエ・ストークス流およびビンガム流にも適用可能であることが示されており、無粘性限界においてベルヌーイの原理を回復する。
亀裂進展の検証： 混合モード荷重下での亀裂キンキング（kinked cracks）への適用において、本論文は、SBENによる予測（暗黙的スキームと法線則を用いたもの）を、実験データ（PMMA試料）およびRichardによる経験則と比較している。
- 明示的オイラー・スキームは、実験と比較して「悲惨な（disastrous）」予測を与えることが示されている。
- 法則（normality law）と組み合わせた暗黙的スキームは、キンク角および応力拡大係数に関して、実験データと「非常に良好な一致」を示す。
- 導出されたモードIIとモードIの破壊靭性比 ( $K_{IIc}/K_{Ic} = \sqrt{3}/2 \approx 0.86$ ) は、従来の理論的提案よりも実験範囲（0.88–0.95）によく一致している。

意義および主張
本論文は、SBEN原理が、滑らかな力学から非平滑な力学、固体から流体、そしてラグランジュ記述からオイラー記述まで、幅広いアプリケーションに適した堅牢な変分アプローチを提供すると主張している。

著者は、このアプローチが単に弱形式を導出するための手法ではなく、新しい構成則をモデリングするための「アイデアの源泉」であることを強調している。その意義は、凸解析とシンプレクティック幾何学のツールを用いて、非関連フロー則や動的な散逸系を含む単一かつ一貫した幾何学的フレームワーク内でこれらを扱える点にある。著者は謙虚に、亀裂キンキングの特定のケースにおいては法則（normality law）と局所対称性の原理が良好な結果をもたらすが、一般的なルールとしては法則が推奨され、破壊力学の問題においては明示的スキームよりも暗黙的スキームが好ましいと述べている。

本研究は、完全に新しい実験データの提示というよりも、先行研究（Brezis, Ekeland, Nayroles、および著者自身の過去の出版物を含む）の統合として提示されているが、既存の実験的ベンチマークに対してその理論的枠組みを検証している。

A family of variational principles of minima for the plasticity, the friction contact and the fracture mechanics