Phase-field modeling of cyclic behavior in quasi-brittle materials: a micromechanics-based approach

技術要約：準脆性材料における周期挙動のフェーズフィールドモデリング：マイクロメカニクスに基づくアプローチ

問題提起
コンクリート、石積み、岩石などの準脆性材料は、繰返し荷重下で複雑な破壊メカニズムを示す。従来の疲労モデルは損傷の蓄積を扱うが、低サイクル疲労において極めて重要な、非弾性歪みの漸進的な蓄積である「ラチェット（累加歪み）」を捉えきれないことが多い。さらに、標準的なフェーズフィールドモデルは、引張と圧縮を区別するために現象論的なエネルギー分割（スペクトル分割や体積・偏差分割など）に依存することが多いが、これはマイクロクラックの開閉という根底にあるマイクロメカニズムを完全には反映していない可能性がある。課題は、マイクロクラックの摩擦滑りとマイクロクラックの進化といった微視的特徴を明示的に考慮しつつ、周期的な塑性とラチェット効果を準脆性材料のフェーズフィールドモデルに統合する、熱力学的に整合したフレームワークを開発することにある。

手法
著者らは、マイクロメカニクスに基づいたフェーズフィールド・フレームワークの拡張を提案している。本手法の核心は以下の通りである：

1. マイクロメカニクスに基づく自由エネルギー： 本モデルは、代表体積要素（RVE）内におけるマイクロクラックの「開いた状態」と「閉じた状態」を区別する。

   • 開いた状態： マイクロクラックは応力を伝達せず、材料は脆性的に振る舞い、塑性による歪みエネルギーの寄与はない。
   • 閉じた状態： マイクロクラックは摩擦滑りを起こす。自由エネルギーには、弾性歪み、塑性歪み、および新たに導入されたラチェット歪み（$\varepsilon_r$）の寄与が含まれる。ラチェット歪みは、微視的な再配列による弾性基準状態の不可逆的なシフトを表す。
   • この定式化は、明示的なエネルギー分割を用いることなく、マイクロクラックの状態に基づいて、非対称な挙動（引張 vs 圧縮）を自然に捉える。

2. 非関連塑性（Non-Associative Plasticity）とラチェット：

   • 塑性モデルには、非関連のドリュッカー・プラガー（Drucker–Prager）降伏条件を採用しており、関連塑性モデルが予測してしまう非現実的な過剰膨張を回避し、摩擦的な体積変化挙動を捉えている。
   • ラチェット歪みテンソルは内部変数として導入され、非関連のフロー則に従って進化し、塑性歪み速度によって駆動されながら、サイクルごとに漸進的に蓄積される。
   • 塑性ポテンシャルは、パラメータ $\beta_\mu$ および $\beta_K$ を通じて、偏差成分と体積成分のラチェットを独立して制御できるように定式化されている。

3. 疲労劣化：

   • 疲労は、疲労履歴変数（$F$）と劣化関数（$h(F)$）を介してモデル化される。
   • 2種類の劣化関数（漸近関数と対数関数）が実装されている。劣化関数は、閾値（$F_0$）を超えると破壊エネルギー（$G_c$）を減少させる。
   • モデルは、モードI（引張）とモードII（せん断／圧縮）の破壊エネルギーを区別し、smoothstep関数を用いてそれらの遷移を滑らかに正規化している。

4. 数値実装：

   • 偏微分方程式の系をデカップルするために、3フィールドのスタガード解法が採用されている。
   • アルゴリズムは、(1) 変位場（力学的平衡）、(2) クラックのフェーズフィールド（ヒストリー場を用いた損傷進化）、(3) 塑性およびラチェット歪み（インプリシットなリターンマッピング・アルゴリズムを使用）を反復的に解く。

主な貢献

• ラチェットの統合： 本研究の主要な貢献は、マイクロメカニクスに基づいたフェーズフィールドモデルへの、熱力学的に整合したラチェット歪みの組み込みである。これにより、非対称な繰返し荷重下での漸進的な非弾性歪みの蓄積のシミュレーションが可能となる。
• マイクロメカニクス的基礎： 本モデルは、マイクロクラックの挙動（開閉 vs 滑り）に根ざした変分構造を保持しており、純粋に現象論的なエネルギー分割に代わる、物理的根拠に基づいた手法を提供している。
  。
• ラチェットモードの独立制御： 本定式化は、体積成分と偏差成分のラチェットを独立して制御するための特定のパラメータを導入しており、複雑な繰返し応答のモデリングに対する柔軟性を提供している。
• 統一されたフレームワーク： 本モデルは、疲労損傷の進化、繰返し塑性、およびラチェットを、単一の連続体フレームワーク内で同時に捉える。

結果
モデルの検証のために、様々な条件下で数値シミュレーションが行われた：

• 単調および繰返し3点曲げ試験： モデルは、コンクリート梁の実験的な荷重－変位曲線を見事に再現した。繰返し荷重下において、クラックの発生、安定した進展、および特定のサイクル数（例では755サイクル）後の突然の破壊を捉えており、観察された疲労挙動と一致している。
• 高サイクル疲労（単一要素）： 漸近関数と対数関数の両方を用いてS-N曲線が作成された。本研究は、予測される疲労寿命が劣化関数の数学的形態に対して敏感であることを明らかにし、特に高歪みレベルにおける対数形式の予期せぬ挙動についても指摘している。
• 材料点のラチェット挙行： 応力制御および歪み制御の繰返し荷重下でのシミュレーションにより、以下の能力が示された：
  • 応力緩和： 歪み制御下において、サイクルごとに応力振幅が減少する現象。
  • ラチェット： 応力制御下において、各サイクルごとに歪みが漸進的に蓄積する現象。
  • モデルはラチェットパラメータ（$\beta_K, \beta_\mu$）に対して敏感であり、これらの値が高いほど、より大きな歪み蓄積が生じることを示した。
• 穿孔試験片（低サイクル疲労）： 中央に穴を持つ正方形の試験片に対し、力制御の引張－圧縮サイクルを適用したところ、応力集中部における局所的な損傷と非弾性歪みの蓄積をシミュレートできることが示された。結果は、1000サイクルにわたる剛性低下と顕著なラチェット歪みの蓄積を示した。

意義と主張
本論文は、提案されたフレームワークが、非弾性歪みの蓄積が支配的な破壊メカニズムである低サイクル疲労領域において、準脆性材料の繰返し挙動を調査するための強力なツールを提供すると主張している。マイクロメカニクスに根ざすことで、従来の現象論的なアプローチと比較して、破壊モード（引張 vs せん断）に対する、より物理的根拠に基づいた記述が可能となる。著者らは、本モデルが、ヒューリスティックなエネルギー分割に頼ることなく、疲労損傷の進化と繰返し塑性の相互作用、およびラチェットという重要な現象を効果的に捉えられることを強調している。本研究は、これらのモデルが実験的研究を補完し、複雑な荷重経路下での材料劣化の正確な予測を提供することで、広範な試験の必要性を減らす可能性があることを示唆している。結論として、ラチェット成分の独立した制御と、定式化の熱力学的整合性が、準脆性構造の長期的な劣化を正確にモデリングするために不可欠であるとしている。