[論文レビュー] Interplay of rearrangements, strain, and local structure during avalanche propagation

Ge Zhang, Sean A. Ridout+1 more|arXiv (Cornell University)|Sep 23, 2020

Sports Dynamics and Biomechanics参考文献 21被引用数 5

ひとこと要約

本論文は、無熱的で準静的せん断によるジャム化した柔らかく丸い粒子系における、局所的粒子構造、再配置、およびアバランチ中におけるひずみの関係を結ぶ構造的・弾塑的フレームワークを提案する。機械学習で得た「軟らかさ」を用いて再配置を予測し、個々の再配置からのひずみ場を分析することで、遠方ひずみが $ r^{-2} $ で減少することを示し、再配置の大きさと塑性ひずみに基づく理論的予測と一致する。これは、不規則な固体におけるアバランチ伝播の洗練されたモデルの妥当性を裏付けている。

ABSTRACT

Jammed soft disks exhibit avalanches of particle rearrangements under quasistatic shear. We follow the avalanches using steepest descent to decompose them into individual localized rearrangements. We characterize the local structural environment of each particle by a machine-learned quantity, softness, designed to be highly correlated with rearrangements, and analyze the interplay between softness, rearrangements and strain. Our findings form the foundation of an augmented elastoplastic model that includes local structure.

研究の動機と目的

• 不規則な固体におけるアバランチ事象の間で、局所的構造、粒子再配置、弾性ひずみの相互作用を理解すること。
• 局所的構造の不均一性を弾塑的モデルの塑性に組み込む予測可能なフレームワークを構築すること。
• 再配置が長距離ひずみ場に与える影響を定量化し、ひずみが局所的構造的環境をどのように歪めるかを明らかにすること。
• 数値シミュレーションを用いて、局所的塑性イベントからの遠方ひずみ減衰の理論的予測を検証すること。

提案手法

• 準静的せん断中における個々の再配置を解明するため、勾配降下法によるエネルギー最小化を用い、過減衰した緩和過程を捉える。
• 隣接粒子ベクトルの変化から計算される非アフィン変位 $ D^{2}_{\text{min}} $ を用いて、粒子再配置を定義する。
• 局所的ペア相関関数に基づく機械学習で得た軟らかさ指標を用い、観測された再配置と相関を最大化するように訓練する。
• 有限時間における弾塑的方程式の解を用いて、個々の再配置からのひずみ場を測定し、体積ひずみ成分と歪み成分を含む。
• 観測されたひずみ減衰をべき関数 $ \sim r^{-2} $ にフィットさせ、数値結果を点源モデルからの解析的予測と比較する。
• 有限持続時間の塑性イベントをモデル化するため、ガウス型時間畳み込みを適用し、数値ひずみプロファイルとの一致を向上させる。

実験結果

リサーチクエスチョン

• RQ1 機械学習で得た軟らかさで定量化される局所的構造の不均一性は、アバランチ中の粒子再配置とどの程度相関するか？
• RQ2 個々の再配置が長距離弾性ひずみ場をどの程度生成するか。また、その場が距離とともにどのように減衰するか？
• RQ3 塑性イベントの空間的・時間的特性が、系全体の観測されるひずみ場にどのように影響を与えるか？
• RQ4 点源としての塑性源からのひずみ場の理論的予測が、実際のアバランチの数値シミュレーションと定量的に一致するか？

主な発見

• 遠方の歪み成分 $ \tilde{\epsilon} $ と $ xy $-ひずみ $ \epsilon_{xy} $ は、いずれも $ r^{-2} $ で減衰し、点源弾塑的モデルからの理論的予測と一致する。
• $ \tilde{\epsilon} $ の測定された比例係数は $ \tilde{c} = 0.03 $ であり、再配置の大きさと平均的塑性ひずみに基づく理論的推定値 $ 0.031 $ に非常に近い。
• $ \epsilon_{xy} $ ひずみプロファイルは、数値シミュレーションにおいて $ r^{-2} $ よりわずかに遅く減衰しており、同時に発生する再配置や有限サイズ効果の干渉を示唆している。
• $ k $-ひずみ成分（せん断応力に関連）は、幅 $ \alpha^{-1/2} \approx 127.83 $ のガウス型時間プロファイルと畳み込むことで理論モデルとよく一致し、塑性イベントの持続時間が長いことを示している。
• 圧縮源からの体積ひずみの理論モデルは、有限時間での数値結果と一致するが、せん断源モデルでは有限時間の畳み込みが必要である。
• 再配置部位における平均的な局所フィット歪み成分は $ \tilde{\epsilon} = 3.6 \times 10^{-3} $、$ \epsilon_{xy} = 1.8 \times 10^{-4} $ であり、導出された比例係数と整合的である。