[論文レビュー] Deformation Measures for Granular Materials

Matthew R. Kuhn|arXiv (Cornell University)|Dec 3, 2024

Geotechnical and Geomechanical Engineering被引用数 25

ひとこと要約

論文は、2D粒状アセンブリの変形をミクロ機械的、void-cell ベースの表現として開発し、Bagi and Satake の枠組みを拡張し、2D DEM シミュレーションによって高度に非均一な変形とマイクロ帯域を観察できることを示す。

ABSTRACT

The paper presents a micromechanical representation of deformation in 2D granular materials. The representation is a generalization of K. Bagi's work and is based upon the void-cell approach of M. Satake. The general representation applies to a material region partitioned into polygonal subregions. This representation possesses a certain consistency that allows for a unique assignment of the contribution that each contact displacement makes to the average deformation of an assembly. The paper addresses construction of the particle graph and appropriate data structures for use with the Discrete Element Method. The approach is applied in a numerical simulation of a two-dimensional assembly of disks. The author presents results of the distributions of deformation and particle-group rotation, with a resolution of about a single particle diameter. Deformation was very nonuniform, even at low strains. Micro-bands, thin linear zones of intense rotation, were also observed.

研究の動機と目的

• 粒状材料の変形に対するミクロ機械的な見方を動機づける。
• 既存の空隙空間分割法をポリゴナルな部分領域へ一般化する。
• アセンブリの平均変形への各接触変位の寄与を割り当てる一貫した方法を確立する。
• Discrete Element Method (DEM) 実装に適した粒子グラフとデータ構造を開発する。
• 2Dディスク・アセンブリの数値シミュレーションで手法を実証し、変形と回転の分布を分析する。

提案手法

• 2D材料領域をポリゴナルなサブ領域（void cells）への分割として表現する。
• 各サブ領域内の頂点と辺の運動を用いて平均速度勾配を表現する。
• 辺相対速度を局所変形に関連づける一貫した Q^m 行列を定義し、再帰的構築を用いる（Eq. 11）。
• Bagi の三角形の場合を m辺ポリゴンへ一般化して局所変形を計算する（Eq. 7）。
• 平面的グラフ構築を説明する（粒子中心に頂点、接触に辺、空隙を面として）と DEM 統合のためのデータ構造（SLL、DLL、DCEL）。
• void-cell の面を構築し、トポロジの一貫性を維持するアルゴリズムを概説する（Eq. 6–11）。

Figure 1: 2D topological representations

実験結果

リサーチクエスチョン

• RQ1 多角形分割された粒状領域の変形を、多角形の頂点と辺の運動とどのように関連づけることができるか。
• RQ2 粒状アセンブリの平均変形への接触変位の寄与を一貫した、かつ一意に割り当てることができるか。
• RQ3 DEMシミュレーションで粒子グラフ/void-cellグラフを構築・更新するのに適したデータ構造とアルゴリズムは何か。
• RQ4 2D数値シミュレーションは粒状アセンブリにおける変形と回転の空間分布について何を明らかにするか。
• RQ5 一般的な荷重下で非均一場やマイクロバンドのような微視的変形特性が現れるか。

主な発見

• 2D粒状アセンブリの変形は、低ひずみでも高度に非均一である。
• シミュレートされたアセンブリには、回転が集中する薄い線状のマイクロバンドが観察される。
• 一貫した多角形ベースの変形測度は Q^m 行列を介して得られ、接触変位の局所変形への一意的寄与を可能にする。
• 本手法は、単一粒径分解能での変形と回転分布の粒子スケールの詳細な視点を提供する。
• DEM対応のデータ構造を備えた平面グラフフレームワークを提案し、接触と空隙のトポロジー追跡を効率化する。

Figure 2: Triangular face: edge, velocity, and relative velocity vectors