A blended approach for evolving phase fields using peridynamics: Cyclic loading in quasi-brittle fracture

この論文は、塑性変形と損傷に伴う弾性係数の低下を考慮した新しい二点履歴依存型位相場モデルを非局所的な構成則と組み合わせることで、準脆性材料の疲労破壊や損傷を熱力学法則に則って予測するメッシュフリー手法を提案し、コンクリートの三点曲げ試験など実験結果と定量的・定性的に一致する結果を示したものである。

要約

この論文は、「コンクリートやレンガのような、少しひび割れやすい素材（準脆性材料）が、繰り返し力を受けてどう壊れるか」を予測する新しい計算方法について書かれています。

専門用語を抜きにして、日常の言葉と面白い例え話を使って説明しましょう。

1. 従来の方法との違い：「2 つの方程式」vs「1 つの方程式」

これまでのコンピュータシミュレーションでは、材料の「変形（動き）」と「壊れ方（き裂）」を予測するために、**2 つの異なるルール（方程式）**を同時に動かしていました。

• ルール A: 物体がどう動くか（ニュートンの運動法則）。
• ルール B: ひび割れがどう広がるか（別の特別なルール）。

これらは別々に計算され、最後に合体させられるので、計算が複雑で、時に矛盾が起きることがありました。

この論文の新しい方法は？
「全部、1 つのルール（ニュートンの運動法則）だけで説明しちゃおう！」という大胆なアイデアです。
ひび割れは、特別なルールで追加するのではなく、「材料の動きそのもの」から自然に生まれてくる現象として扱います。

例え話：
従来の方法は、「車の運転（動き）」と「車の故障（壊れ方）」を、それぞれ別の運転手と整備士が別々のマニュアルで管理しているようなものです。

新しい方法は、**「運転手がハンドルを切った瞬間、自動的にタイヤがすり減り、最終的にタイヤが外れるまでを、運転手一人の動きだけで全て説明する」**ようなものです。動きさえ正しければ、壊れ方も自然に決まります。

2. 核心となるアイデア：「過去の記憶」と「2 点の絆」

この方法の最大の特徴は、**「材料の歴史（記憶）」**を重視している点です。

• 2 点の絆（ボンド）：
  材料を無数の小さな点の集まりだと考えます。隣り合う 2 つの点の間に、目に見えない「ゴム紐（ボンド）」が張られていると想像してください。
• 記憶（ヒストリー）：
  このゴム紐は、「過去にどれくらい引っ張られたか」を覚えています。
  • 一度強く引っ張られて伸びきると、元には戻りません（塑性変形）。
  • さらに引っ張り続けると、ゴム紐の強さが弱くなり（損傷）、最後はパキッと切れます（破壊）。

この「ゴム紐の強さ」は、**「フェーズフィールド（状態の指標）」**という数値で表されます。

• 1.0: 新品のゴム紐（元気）。
• 0.5: 伸びて弱くなったゴム紐（損傷中）。
• 0.0: 完全に切れたゴム紐（ひび割れ発生）。

例え話：
材料を**「巨大な蜘蛛の巣」だと想像してください。
蜘蛛の糸（ゴム紐）は、過去にどれだけ虫がぶつかったか（負荷）を記憶しています。
一度ぶつかった糸は、次にぶつかった時にすぐに切れてしまいます。
このシミュレーションは、「風（外力）が吹いた時に、どの糸が伸びて、どの糸が切れて、最終的にどんな穴（ひび割れ）が開くか」**を、糸の動きだけで計算しちゃうのです。

3. 「引っ張り」と「押し付け」の使い分け

コンクリートのような素材は、**「引っ張られると弱い」けれど、「押し付けられると強い」**という特徴があります。
このモデルは、その違いを完璧に再現しています。

• 引っ張り（引張力）： 糸が伸びて弱くなり、最終的に切れます（ひび割れ）。
• 押し付け（圧縮力）： 糸が切れても、隣り合う点が押し合うと、また「弾力」で支えようとします。つまり、ひび割れた部分でも、押し付けば再び力を伝え合えます。

例え話：
壊れた陶器の皿を想像してください。

• 引っ張ると： 割れた部分が離れて、二度と元に戻りません（ひび割れ）。
• 上から押すと： 割れた部分が互いに押し合い、一時的に支え合います。
  このシミュレーションは、**「割れた皿でも、押せば支えられる」**という現実の挙動を、計算の中で自然に再現しています。

4. なぜこれがすごいのか？（実験との一致）

この新しい方法は、以下の点で画期的です。

1. エネルギーの保存： 計算の中で「エネルギーが勝手に増えたり消えたりしない」ことを数学的に証明しています。熱力学の法則に厳密に従っています。
2. 実験との一致： コンクリートの「3 点曲げ試験（真ん中に力を加えて折る実験）」や、**「繰り返し力を加えた時のヒステリシス（力と変形のループ）」**という複雑な現象を、実験データとほぼ同じように予測できました。
3. サイズ効果の再現： 「大きな梁（はり）ほど、単位面積あたりの強度が下がる」という、コンクリート特有の不思議な現象（サイズ効果）も、特別な調整なしに自然に現れました。

例え話：
これまでの計算は、**「大きなビルと小さな模型を同じ強度だと仮定して計算」してしまいがちでした。
しかし、この新しい方法は、「ビルは重すぎて、小さなひび割れがすぐに全体を崩壊させるが、模型は丈夫だ」**という、現実の「サイズによる違い」を、計算の仕組みそのものから自然に導き出します。

まとめ

この論文は、**「材料の壊れ方を、特別なルールで無理やり追加するのではなく、材料の『動き』と『過去の記憶』から自然に導き出す」**という、シンプルで美しいアプローチを提案しています。

• 従来の方法： 動きと壊れ方を別々に管理する（複雑）。
• この論文の方法： 動きさえ正しければ、壊れ方は自動的に決まる（シンプルで自然）。

これにより、コンクリート構造物の耐久性予測や、地震後の損傷評価など、実社会で非常に役立つ予測が、より正確に、かつ計算コストを抑えて行えるようになる可能性があります。まるで、**「材料の心（歴史）を読み解くことで、未来の破損を予言する」**ような技術なのです。

技術概要

論文要約：準脆性材料における塑性変形と損傷を考慮したペリダイナミクスと位相場を融合したアプローチ

1. 背景と課題

準脆性材料（コンクリート、岩盤、煉瓦など）の破壊現象を記述する際、従来の連続体力学における古典的な破壊力学や、既存の位相場法（Phase Field Method）にはいくつかの課題がありました。

• 既存の位相場法の限界: 従来の動的破壊モデルでは、変位の進化（ニュートンの第二法則）と、位相場（損傷の進化）を記述する独立した時間発展方程式（通常はギンツブルグ・ランダウ方程式など）を連成させる必要があります。これにより、数値計算の複雑さが増し、熱力学的整合性を保つための追加的な制約条件が必要となることがあります。
• 塑性と損傷の統合: 準脆性材料では、弾性変形に加え、不可逆な塑性変形（ひずみ硬化や軟化）が破壊過程に大きく寄与します。特に、荷重の繰り返し（サイクリック荷重）による履歴現象（ヒステリシス）や、引張と圧縮における異なる挙動を、単一の理論枠組みで記述することは困難でした。
• メッシュ依存性: 従来の有限要素法では、破壊の局所化に伴うメッシュ依存性の問題が発生しやすく、メッシュフリー手法の必要性が指摘されていました。

本研究は、これらの課題を解決し、ニュートンの第二法則のみを用いて、変位、損傷、塑性ひずみを統一的に記述する新しい場理論（Field Theory）を提案することを目的としています。

2. 提案手法：融合アプローチ（Blended Approach）

著者らは、ペリダイナミクス（Peridynamics）の非局所性（Non-locality）と、位相場概念を融合させた新しい定式化を提案しました。

2.1 理論的枠組み

• ニュートンの第二法則のみの使用: 位相場の進化を記述する独立した微分方程式は導入されません。代わりに、位相場は「履歴依存の構成則（Constitutive Law）」の一部として、変位場によって自動的に決定されます。
• 2 点履歴依存位相場: 材料の劣化モデルを記述するために、2 点間の履歴に依存する位相場を導入しました。これにより、引張力と圧縮力による 2 つの独立した損傷状態を表現します。
  • 引張位相場 ($\gamma_+$): 引張損傷と塑性ひずみを記述。
  • 圧縮位相場 ($\gamma_-$): 圧縮損傷と塑性ひずみを記述。
• ひずみの分解: 2 点間のひずみは、弾性ひずみと不可逆（塑性）ひずみの和として定義されます。
• 構成則と履歴: 材料の応力 - ひずみ関係は、強度エンベロープ（Strength Envelope）と、損傷・塑性に基づく弾性係数の軟化を考慮したアンローディング則（Unloading laws）で記述されます。
  • 引張領域では、ピーク応力以降、塑性変形と弾性軟化が発生します。
  • 圧縮領域では、弾性挙動を維持しつつ、特定の条件下で損傷が進行します。
  • アンローディング比 ($\beta$): 材料軟化に伴う塑性ひずみと、軟化を無視した塑性ひずみの比率として定義され、脆性破壊から準脆性破壊までの挙動を制御するパラメータとなります。

2.2 数値的実装

• メッシュフリー法: ペリダイナミクスに基づき、要素や節点間の幾何学的接続を必要とせず、積分近似（クアドラチャ）に基づく離散化手法を採用しています。
• 初期境界値問題: 変位場、弾性応力、弾性ひずみ、および不可逆ひずみ場が一意に決定される、数学的に適切に定式化された（Well-posed）初期境界値問題として定式化されています。
• 収束性: 非局所性のスケール（$\epsilon$）がゼロに近づく極限において、このモデルはクラシックな線形弾性力学およびグリフィスの破壊エネルギー則に収束することが理論的に示されています。

3. 主要な貢献と理論的性質

1. 熱力学的整合性の保証: 運動方程式から直接、エネルギーバランス（運動エネルギー、弾性エネルギー、散逸エネルギー）を導出しました。損傷率（エネルギー散逸率）は常に非負であり、熱力学の法則を満たします。
2. 破壊エネルギーの導出: 平坦な亀裂に対する破壊エネルギーが、エネルギー解放率と亀裂面積の積としてモデルから直接導出され、非局所性の長さスケールに依存しないことが示されました。
3. 塑性ヒステリシスの記述: 引張と圧縮の両方における塑性変形と損傷の履歴を、単一の構成則内で自然に記述可能であり、サイクリック荷重下でのヒステリシスループを再現できます。
4. サイズ効果の予測: 材料定数のみから、構造物のサイズに依存する破壊強度（サイズ効果）を、明示的なサイズ依存ルールなしに予測可能であることを示しました。
5. パラメータの独立性: 非局所性の長さスケール（$\epsilon$）は材料定数ではなく、数値解の精度を制御するパラメータとして扱われ、強度や破壊エネルギーの較正には使用されません。

4. 数値実験結果

本研究では、コンクリート試料を対象とした複数のベンチマーク問題で数値シミュレーションを行い、実験データと比較しました。

• メッシュ収束性: 3 点曲げ試験において、ホライズンサイズ（$\epsilon$）とグリッド比（$m$）を変化させた収束解析を行い、提案モデルの精度と安定性を確認しました。
• モード I 破壊（単調・サイクリック荷重）:
  • 高強度コンクリートの 3 点曲げ試験において、荷重 - クラック開口変位（CMOD）曲線が実験結果と定量的に一致しました。
  • 繰り返し荷重下でのヒステリシス挙動を正確に再現し、塑性ひずみの局所化（プロセスゾーン）を可視化しました。
• 混合モード破壊:
  • 片持ち梁や L 字型コンクリートパネルにおけるせん断と引張が複合する破壊経路を予測し、実験で観測された亀裂経路と良好な一致を示しました。
• サイズ効果:
  • 異なる寸法のコンクリート梁（3 種類のサイズ）を用いた試験において、梁のサイズが大きくなるにつれて見かけの強度が低下する「サイズ効果」を、同一の材料定数と構成則で再現することに成功しました。これは、従来の塑性破壊理論や最大応力則では説明できない現象です。

5. 意義と結論

本研究で提案された「ペリダイナミクスと位相場を融合したアプローチ」は、準脆性材料の動的破壊、特に塑性変形を伴う複雑な履歴挙動を記述する強力なツールとなります。

• 理論的革新: 独立した位相場進化方程式を不要とし、ニュートンの第二法則と履歴依存の構成則のみで破壊を記述する枠組みは、計算効率と理論的整合性の面で画期的です。
• 実用性: 実験データ（荷重 - 変位履歴、破壊経路、サイズ効果）と高い一致を示しており、コンクリート構造物の損傷評価や破壊予測への応用が期待されます。
• 将来展望: この枠組みは、異方性材料や複合材料への拡張、およびテンソル値の位相場を用いたより一般的な状態ベースのペリダイナミクスへの発展の可能性を開いています。

要約すると、この論文は、熱力学的に整合性を持ち、メッシュフリーで、かつ塑性と損傷の複雑な相互作用を自然に扱う新しい破壊力学の計算手法を確立し、その有効性を数値的・実験的に実証した重要な研究です。