[论文解读] MC-DEM: a novel simulation scheme for modeling dense granular media

Nicolas Brodu, J.A. Dijksman+1 more|arXiv (Cornell University)|Oct 23, 2014

Granular flow and fluidized bedsEngineering参考文献 16被引用 19

一句话总结

本文提出MC-DEM，一种新型的离散元方法，通过引入基于线弹性理论的形变模型，显式模拟同一颗粒上多个接触点之间的相互作用，从而提升致密颗粒堆积体中的力预测精度。该方法在力分布和标度律方面与三维实验数据高度吻合，显著优于将接触点视为独立处理的传统DEM，且无需采用有限元方法。

ABSTRACT

This article presents a new force model for performing quantitative simulations of dense granular materials. Interactions between multiple contacts (MC) on the same grain are explicitly taken into account. Our readily applicable method retains all the advantages of discrete element method (DEM) simulations and does not require the use of costly finite element methods. The new model closely reproduces our recent experimental measurements, including contact force distributions in full 3D, at all compression levels up to the experimental maximum limit of 13\%. Comparisons with traditional non-deformable spheres approach are provided, as well as with alternative models for interactions between multiple contacts. The success of our model compared to these alternatives demonstrates that interactions between multiple contacts on each grain must be included for dense granular packings.

研究动机与目标

• 解决传统DEM模拟与致密颗粒堆积体实验数据之间定量一致性差的问题。
• 建模单个颗粒上多个接触点之间的相互影响，该影响在标准DEM中被忽略。
• 开发一种计算效率高的方法，在保持DEM简单性的同时，提升致密、近阻塞体系的模拟精度。
• 通过完整的三维实验测量数据（接触力与宏观响应）验证模型。
• 证明颗粒尺度形变与跨接触力耦合对致密颗粒介质准确模拟至关重要。

提出的方法

• 通过球体上点力的解析线弹性解计算局部颗粒形变。
• 引入校正项δk→c，表示同一颗粒上其他接触点k的力在接触点c处引起的法向位移。
• 有效重叠量更新为δc + Σk δk→c，随后用于赫兹力定律F ∝ (δc + Σk δk→c)^3/2。
• 引入经验系数γ以考虑边界条件对跨接触影响的作用，提升与实验数据的拟合效果。
• 所有形变均在颗粒层面计算，无需长程相互作用或有限元方法（FEM）计算。
• 该方法可推广至其他已知形状-形变关系的颗粒体系，如乳液或泡沫。

实验结果

研究问题

• RQ1 能否通过一种考虑同一颗粒上多接触相互作用的改进DEM模型，更好地再现致密三维颗粒堆积体中的实验力分布？
• RQ2 通过弹性形变引入的跨接触力耦合，如何影响致密颗粒材料的宏观与微观响应？
• RQ3 标准DEM中独立接触点假设在多大程度上无法捕捉接近阻塞状态致密颗粒体系的真实力学行为？
• RQ4 能否通过一种简单、解析的多接触效应校正方法，替代昂贵的FEM模拟，同时保持高精度？
• RQ5 颗粒形变与力传播如何影响接触力随堆积密度和应变的标度关系？

主要发现

• MC-DEM与3DXP实验研究的三维实验数据高度一致，包括完整的三维力分布和宏观力-应变响应。
• 该模型以高保真度匹配实验标度律F ∝ f Z φ b，显著优于标准DEM，后者表现出明显偏差。
• MC-DEM模拟在卸载阶段仅呈现轻微曲率，比标准DEM更接近实验数据，尽管仍存在微小滞后现象。
• MC-DEM中最大顶部力F与实验值Rtop在小Fmin偏移范围内一致，表明具有极强的定量一致性。
• 该模型成功捕捉了由颗粒形变引发的新接触形成，这是标准DEM中缺失的关键机制。
• 经验系数γ有效补偿了边界条件效应，表明未来模型可利用颗粒层面统计量来模拟复杂约束。