Computation of stresses in jammed packings modeled with Tresca friction

技術要約：Tresca摩擦を用いたジャミング・パッキングにおける応力計算

問題提起
本論文は、石造建築や高密度な粒状体のような、剛体多角形セルからなるジャミング・パッキング（jammed packings）内の応力場の計算について取り組んでいる。貫入を抑制するために弾性変形を利用することが多い従来の離散要素法（DEM）とは異なり、本研究はセルが完全に剛体であり、弾性的に変形しない系に焦点を当てている。このような「ジャミング」状態において、機械的な平衡は、有意な弾性歪みなしに、接触力と摩擦のみを通じて維持される。具体的な課題は、力の変数が平衡制約の数よりも多い静定不可能な（hyperstatic）システムにおいて、内部応力分布を決定することにある。著者らは、接線方向の滑りの閾値が固定されていることを特徴とするCoulomb摩擦とは異なる、Tresca摩擦則を用いて接触物理をモデル化している。

手法
提案される手法は、以下の3段階の計算フレームワークに従う。

1. 原始エネルギー最小化（Primal Energy Minimization）:
   著者らは、全摩擦エネルギーを最小化する制約付き最小化問題を定式化している。エネルギー汎関数は、すべての内部界面におけるTresca摩擦仕事の総和で構成され、$s_T |[u]_e \cdot t_e|$ と定義される。ここで、$s_T$ は摩擦係数、$[u]_e$ は変位ジャンプ、$t_e$ は接線ベクトルである。この最小化は以下の条件に従う：

• 非貫入制約： 変位ジャンプの法線成分は非負でなければならない（$[u]_e \cdot n_e \geq 0$）。
• 閉じ込め： エネルギーには弾性復元力が存在しないため、セルが無限遠へ滑り出すのを防ぐために、ノイマン境界条件（$g$）が適用される。
• 並進不変性： 重心の位置を固定するために、グローバルな制約が加えられる。
  セルは、並進の自由度のみを持つ剛体としてモデル化されており、回転の自由度は明示的に無視されている。この簡略化は、ジャミング状態において、回転に対する抵抗（ダイラタンシーを必要とする）は、滑りに対する抵抗よりも数桁高いという観察に基づいている。

2. 双対最大化（Dual Maximization）:
   原始問題が変位（安定したジャミング状態では実質的にゼロ）を求めるものであるため、著者らは双対最大化問題を導出する。双対変数は、内部エッジに作用する力に対応する。この双対問題を解くことで、界面における法線応力および接線応力が得られる。適切な閉じ込め項に関する仮定の下で、原始エネルギーの最小化子の存在および双対問題の解の存在が確立されている。双対定式化は、可能な応力状態のハイパスタティック・多様体から特定の解を選択する役割を果たす。

3. 応力再構成（Stress Reconstruction）:
   双対問題の出力は、エッジにおける離散的な応力値を提供する。各セル内の連続的な応力場 $\sigma$ を再構成するために、著者らは最低次（lowest-order）のRaviart–Thomas ($RT_0$) 有限要素を用いる。これらの要素は $H(\text{div})$ 適合であり、その自由度はエッジ重心における法線応力に対応する。再構成には、各セル内での局所的な混合問題を解くプロセスが含まれ、以下の条件を満たす応力テンソル $\sigma_c$ を求める：

• $\text{div}(\sigma_c) = 0$ （体積荷重がない場合の平衡）。
• 境界上の法線応力が、双対解から得られた値と一致すること。
• 応力テンソルの対称性（$\sigma = \sigma^T$）は、ラグランジュ乗数を用いて弱形式で課される。これは、基礎となる離散モデルがモーメント平衡方程式を無視しているため、再構成された応力が残留トルクを表す歪対称成分を持つ可能性があることに起因する。

主な貢献

• 剛体摩擦モデル： 本論文は、標準的なDEMのアプローチとは異なり、弾性パラメータを用いずに剛体摩擦（Tresca基準）のみに依存する、ジャミング・パッキングのためのモデルを導入している。
• ハイパスタティシティ（静定不可能性）への対処： 本研究は、ジャミング・パッキングの静的決定不全性に明示的に取り組んでいる。最適化アプローチが、摩擦エネルギー汎関数を最小化することによって、自己応力状態の空間から一意の解を選択することを実証している。
• 一貫した応力再構成： 離散的な界面力から一貫した応力場を生成するために、$RT_0$ 要素を用いた新しい再構成技術を提示し、各セル内での局的な力平衡を保証している。
• 回転の無視： 著者らは、圧縮されたパッキングでは並進滑りが支配的であることを理由に、回転の自由度を無視する簡略化されたモデルを提案しており、これが非対称な応力テンソルをもたらすことを認めている。

数値結果
本手法は、様々な2D多角形メッシュ（具体的にはボロノイ・テセレーション）を用いてテストされ、以下の知見が得られた。

• 検証： 一様な境界条件を用いたテストにおいて、再構成された応力場は期待される一様な応力状態と一致し、手法の一貫性が確認された。
• せん断テスト： せん断荷重シナリオにおいて、再構成された応力テンソルは非対称性（$\sigma_{12} \neq \sigma_{21}$）を示した。著者らはこれを、離散モデルにおいてモーメント平衡が強制されていないことに起因するとし、歪対称成分が残留トルクを説明すると述べている。
• 均質圧縮： 一様な圧縮境界条件にもかかわらず、ソルバーは非均質な応力場を生成した。著者らはこれを、システムのハイパスタティシティの結果であると説明している。すなわち、最適化が零空間（null space）から特定のエネルギー最小化状態を選択するため、自然に「力鎖（force chains）」（高応力の疎なネットワーク）や低応力またはゼロ応力の領域が形成される。
• 圧縮中の引張： 結果は、全般的な圧縮下においてもセル内に引張（正の主応力）が存在することを示した。著者らは、これを数値的なアーティファクトではなく、摩擦によるインターロッキング（噛み合わせ）やブリッジング効果から生じる物理現象であると特定している。

意義と主張
本論文は、提案された手法が、剛体摩擦の物理とエネルギー最小化を活用することで、ジャミング・パッキングにおける応力を計算するための堅牢かつ一貫した方法を提供すると主張している。著者らは、弾性ベースのDEMと比較して、本アプローチは弾性変形が無視できる領域において、ジャミング系の基礎となる物理をより良く捉えられると論じている。

結果の意義は、摩擦エネルギー（$L_1$ 型の目的関数）を最小化することが、自然に疎な力ネットワークへと導くことであり、これが粒状材料における力鎖の実験的観察と一致している点にある。著者らは、その限界についても謙虚に認めている。すなわち、回転平衡の無視が非対称な応力テンソルをもたらすこと、および、Trescale摩擦（固定閾値）の使用が、より複雑で非凸な準変分不等式（Quasi-Variational Inequalities）を必要とするCoulomb摩擦と比較して簡略化されていることである。著者らは、将来の研究において、完全な機械的一貫性を達成するために、回転の自由度やより現実的な摩擦則を含むモデルへと拡張できることを示唆している。