Plasticity as Directional Stationarity: Yielding, Flow, and Hardening from One Functional

技術要約：方向性定常性としての塑性

問題提起
従来の塑性理論は、通常、弾性則、降伏条件、流動則、硬化関係、および負荷・除荷条件という、個別の構成要素の集合体として整理されている。この枠組みは計算機への応用において有効であることが証明されているが、共通のエネルギー的および方向的な起源を持つ関係を分離して扱っている。本論文は、これらの構成則を単一の生成対象から導出する統一的な変分言明の必要性に対処し、関連流動（associated flow）と非関連流動（non-associated flow）、混合硬化、勾配効果、および速度依存性の扱いを統一することを目的としている。

手法
本論文は、一連の片側的な許容塑性パス上で評価される単一のスカラー汎関数 $\Pi$ の**方向性定常性（directional stationarity）**に基づく定式化を提案する。

• スカラー汎関数: 定式化の核となるのは、弾性エネルギー $\phi$、塑性状態エネルギー $\psi_p$、および活動抵抗 $W$ からなる領域 $\Omega$ 上で定義される汎関数である。
  $$\Pi[u, \varepsilon^p, \lambda] = \int_\Omega [\phi(\varepsilon - \varepsilon^p) + \psi_p(\varepsilon^p) + W(\lambda)] dx - \ell(u)$$
• 方向性定常性: 降伏則や流動則を個別に課す代わりに、本論文では、汎関数の一次変分を許容される塑性接錐（tangent cone）に制限する。許容される変分は、非負の活動増分 $\delta\lambda_a \geq 0$ と方向ベクトル $M_{Ia}$ によってパラメータ化される。
• 構成則の導出:
  • 一次変分により、応力 $\sigma$、内部変数力 $X_I$、および抵抗 $R$ が得られる。
  • この変分を許容される接錐に制限することで、各メカニズムに対する縮退方向力（reduced directional force） $F_a$ が定義される：
    $$F_a = \sum_I X_I \cdot M_{Ia} - R_a$$
  • 非負の錐に対する方向性定常性の条件は、直接的にクーン・タッカー（Kuhn–Tucker）条件を導く：
    $$F_a \leq 0, \quad \delta\lambda_a \geq 0, \quad F_a \delta\lambda_a = 0$$
  • これらの条件は、降伏境界（$F_a=0$）、弾性不等式（$F_a < 0$）、および負荷・除荷の相補条件を同時に定義する。

主な貢献

1. 統一された変分フレームワーク: 本論文は、降伏関数が「単位あたりの許容塑性活動に対する負の一次変分」と識別できることを示している。このアプローチは、一つの汎関数から応力、内部力、および抵抗の微分を生成し、許容される接錐の集合が不可逆的なメカニズムを定義する。
2. 関連流動および非関連流動:
   • 関連流動: 許容される接錐が、許容される塑性計量（支持関数の劣微分）の支持方向として選択される場合に自然に生じる。この場合、流動方向は降伏面の法線方向となる。
   • 非関連流動: 接方向を、結果として得られる降伏境界に対して法線ではないものとして割り当てることで対応する。本論文は、方向的な仕事が正である限り、正の再スケーリングを介した、規則的な非関連方向の局所的な表現を提供している。
3. 自己相似的な関連ファミリー: 自己相似性解析により、正の同次的な応力ゲージが広範な降伏面ファミリーを生成することが明らかになる。本論文は、単一のゲージ $\rho$ が面の形状と流動方向を固定し、抵抗項 $R$ がスケール（等方的な膨張）を制御し、塑性状態エネルギー $\psi_p$ が中心（運動学的な移動）を制御することを示しており、幾何学的形状、サイズ、および移動の役割を分離している。
4. 一般化された硬化および勾配: 本フレームワークは以下を自然に組み込む：
   • 混合硬化: 等方硬化、運動学的な硬化、および結合硬化は、抵抗 $W$ と塑性状態エネルギー $\psi_p$ の異なる微分から導出される。
   • 勾配塑性: 抵抗汎関数に勾配項（例：$W(\lambda, \nabla\lambda)$）を含めることで、非局所的な抵抗項と自然な微視的境界条件を一次変分から直接導出できる。
   • マルチメカニズム塑性: 独立した活動変数を用いることで、潜在的な硬化結合を伴う複数の独立した活性化メカニズム（例：結晶滑り系や多表面モデル）のモデリングが可能となる。
5. 粘塑性への拡張: 粘性ポテンシャル $V(\dot{\lambda})$ を導入することで、速度依存的な発展へと拡張される。粘性スケールがゼロに近づくにつれ、レート非依存の極限が標準的なKKTグラフへと収束する。

結果と解析例
本論文は、4つの閉じた形式の解析解を通じて、本定式化を検証している：

• 多閾値ねじり: 円筒管のねじりにおける空間的に分離した活性環（annuli）の形成を示し、複数の閾値と硬化係数がどのように異なる活性半径を作り出すかを説明する。
• 運動学的硬化を伴う不均一棒: 先細りの棒に関する明示的な解を提供し、単調負荷下での背応力場と塑性ひずみがどのように空間的に分布するかを示す。
• ヒル型アニュラス（円環）: 厚肉円環の内部圧力下における圧力感受的な限界状態を再現し、圧力依存的な強度境界（パラメータ $\alpha$ によって制御）と塑性ダイラタンシー（パラメータ $\beta$ によって制御）を明確に分離する。
• 弾塑性キャビティ膨張: 円筒状空洞の膨張を解き、強度と流動キネマティクスのデカップリングを示す。この解は、応力場は降伏面の形状（$\alpha$）に依存する一方で、変位場と塑性ダイラタンシーは流動方向（$\beta$）に依存することを示しており、非関連流動のモデリングを可能にしている。

意義と主張
本論文は、この方向性定常性定式化が、伝統的に分離されてきた関係に対して共通の変分言明を提供すると主張している。降伏関数を負の一次変分として識別することで、このフレームワークは構成則の導出を統一する。また、降伏面の幾何学的形状をそのスケールおよび移動から明示的に分離しており、アドホックな修正なしに非関連流動を扱うための厳密なメカニズムを提供している。

著者は、この定式化が混合硬化、複数の塑性活動、勾配抵抗、および速度依存性を単一の構造内でカバーしていると述べている。解析例は補完的なチェックとして機能し、既知の限界状態と、正規流動および非正規流動に関連する明確に異なる変位場を正しく再現していることを確認している。結論として、最も直接的な拡張には、特定の非関連および勾配モデルのウェル・ポーズド性（well-posedness）の確立、一貫した有限要素線形化の開発、および周期的・速度依存的な実験に対するモデルの較正が含まれる。