On representation of macroscopic crack in periodic fine-scale discrete mechanical models

技術要約：周期的な微細スケール離散機械モデルにおけるマクロな亀裂の表現について

問題提起
コンクリートのような不均質な軟化材料のマルチスケールモデリングにおいて、微細スケールの代表体積要素（RVE）に適用される境界条件（BCs）は、歪みの局所化の発生とパターンに大きな影響を与える。標準的な周期境界条件（PBCs）は、弾性または硬化材料には適しているが、軟化材料においては、局所化バンドの傾斜がモデル幾何学（例：正方形または立方体のドメイン）の周期方向と一致しない場合、人工的な延性挙動と過剰なエネルギー散逸を引き起こす。この不一致により、周期性を満たすために複数の非物理的な亀裂パターンが形成され、マクロな応答が歪められる。球形や円形モデルは方位依存性のない幾何学形状を提供するが、標準的なPBCsでは局所化バンドがその境界を横切ることを依然として妨げる。連続体モデルに対しては新しい境界条件が提案されているものの、それらが離散機械モデルに適用可能であるかは未検証である。

手法
本研究では、コンクリートのメソスケール離散粒子モデル（具体的には格子離散粒子モデル、LDPMフレームワーク）に対し、3つの新しい境界条件アプローチの適用可能性を評価する。モデルは、粒子間の接触において構成則が定義された、相互に連結された剛体粒子で構成されている。調査は、様々な荷重角における一軸引張を受ける2次元の正方形および円形モデル、ならびに選択された3次元立方体モデルを対象としている。

本研究では、以下の境界条件戦略を比較検討する：

1. 標準的な周期境界条件 (Standard Periodic BCs): 境界ノードがマクロな歪みに基づいて互いに相対的に移動するように制約される、ベースラインとなる手法。
2. 最小運動学的境界条件 (Minimal Kinematic BCs): 平均的な微細スケール歪みをゼロに強制するが、剛体運動は許容する制約。
3. 浸透経路整合境界条件 (Percolation-path-aligned BCs / Aligned BCs): 周期性を、局所化バンドの進化に合わせて回転させる。これらは、亀裂に垂直な境界セクションには強い周期性を課し、法線が一致しないセクションには弱い（積分的な）周期性を課す。
4. テセレーション境界条件 (Tessellation BCs): 周期性をシフト（オフセット $\Delta x$）させることで、隣接するドメイン内で亀裂のパターンを繰り返すことなく、局所化バンドが継続できるようにし、変形したモデルを用いて空間をテセレーション（敷き詰め）する。
5. 埋め込み亀裂境界条件 (Embedded-crack BCs / 円形モデル用): 特定の「開口ノードペア」に未知の変位ジャンプを導入することで、円形ドメインを分離させるように改良されたPBCs。

モデルは、「平滑な」表面幾何学（補助的な境界ノードを伴う）と、「粗い」表面幾何学（テセレーションがバルクと統計的に同一であるもの）の両方を利用している。シミュレーションはOpen Academic Solver (OAS) を用いて実行され、2次元正方形モデルでは20回の実現、円形および3次元立方体モデルでは10回の実現が行われた。

主な結果

• 標準的な周期境界条件: 方位依存的な軟化を生じさせることが確認された。2Dにおける $0, \pi/4, \pi/2$（または3Dにおける特定の直交方向）以外の荷重角では、人工的に長い複数の亀裂が形成されるため、エネルギー散逸が大幅に増加する（基準値の最大2.5倍）。
• 最小運動学的境界条件: 軟化モデルには不適切であることが判明した。少数のノードの脱離がエネルギー的に有利になる境界部で、偽の歪み局所化が発生し、その結果、非現実的に低いエネルギー散逸とピーク応力を招く。
• 浸透経路整合境界条件: 重大な欠陥を示した。強く制約された境界セクションと弱く制約された境界セクションの区別が、不均一な歪みを生じさせた。これにより、特に荷重角が $\pi/8$ から $3\pi/8$ の間で、複数の局所化バンド（二次亀裂）の形成や、弱く制約されたセクションでの偽の局所化が頻発した。結果はモデル内部の不均質性に強く依存していた。
• テセレーション境界条件: すべての荷重方向において、単一の明確な局所化バンドを一貫して生成した。バンドの長さは、モデルの幾何学形状と荷重角のみによって決定された。その結果、観察された散逸エネルギーの荷重方向への依存性は、境界条件のアーティファクトではなく、幾何学的要因（亀裂の長さ）に起因することが示された。2Dでは、エネルギー散逸は理論的な幾何学的因子 ($1/h$) に密接に一致した。3Dでは、離散化の影響により相関はやや弱まったが、妥当なものであった。
• 埋め込み亀裂境界条件 (円形モデル): 一貫性のない結果をもたらした。特定の「開口ノードペア」で亀裂が発生した場合には完全な局所化バンドを許容したが、その結果は内部の不均質性に左右された。亀裂が他の場所で発生した場合、モデルは標準的なPBCs下と同様の挙動を示し、複数の亀裂形成と過剰なエネルギー散逸を招いた。

意義および結論
本論文は、離散モデルにおいてテセレーション境界条件が標準的な手法よりも最も顕著な改善をもたらすと結論付けている。これらは、エネルギー散逸が境界条件の人工的な制約ではなく、局所化バンドの幾何学形状によって支配される、一貫したマクロ応答を提供する。これにより、計算された亀裂長に基づく簡便な後処理補正が可能となる。

対照的に、浸透経路整合境界条件は、複数の亀裂や偽の局所化を誘発する傾向があるため、離散モデルにおいては信頼できないとされる。埋め込み亀裂境界条件（円形モデル用）は、亀裂の発生位置に依存するため、結果が確率論的であり、ランダムな内部構造に依存するという限界がある。

著者らは、これらの結論は、均質な材料パラメータを持つ離散モデルにおける単調一軸引張に特有のものであると述べている。これらの手法の、より複雑な荷重シナリオや不均質な材料分布を持つ連続体モデルへの適用可能性は、依然として未解決の課題である。さらに、これらの新しい境界条件の実装には、局所化の開始時における亀裂法線を決定するための堅牢な手法が必要であり、それは反復解法において必ずしも容易ではない。