[论文解读] Sound wave velocities in dry and lubricated granular packings packings: numerical simulations and experiments

Ivana Agnolin, Jean-Noël Roux+3 more|arXiv (Cornell University)|Jan 20, 2009

Granular flow and fluidized beds参考文献 15被引用 5

一句话总结

本研究结合数值模拟与实验，表明尽管润滑颗粒堆集的密度较低，但其声波传播速度显著高于干燥堆集，原因在于其配位数更高。模拟结果表明，制备方法——尤其是润滑——对微观结构的影响大于密度，而配位数是影响刚度和波传播的主导因素，挑战了传统认为密度单独决定颗粒材料弹性性质的假设。

ABSTRACT

Numerical simulations are used to investigate the origins of the different wave velocities measured in dense granular samples assembled with different methods. Glass bead packings are prepared in the lab either by pouring and vibrating the dry material in a container, or by mixing with a very small amount of a viscous lubricant. Lubricated samples, although less dense, exhibit significantly higher wave velocities for confining pressures in the 100 kPa range. Numerical predictions for elastic moduli agree much better with experimental results when the computational preparation of the samples mimics the laboratory one, albeit in a simplified manner. A plausible explanation to the laboratory observations is that the coordination number, which influences the material stiffness more than its density, is notably higher in lubricated packings.

研究动机与目标

• 理解尽管密度较低，为何润滑颗粒堆集的弹性波速度高于干燥堆集。
• 研究不同的样品制备方法（干倒法与润滑混合法）对配位数和刚度等微观结构特性的影响。
• 评估数值模拟是否能准确再现颗粒堆集中的实验波速测量结果。
• 确定在致密颗粒材料中，配位数还是密度是影响弹性波速度的主导因素。
• 评估简化数值模型在捕捉真实实验室制备过程对颗粒微观结构和力学响应影响方面的有效性。

提出的方法

• 采用离散元动力学方法进行数值模拟，对实验室制备方法进行简化建模：'振动'干燥（C型）和'润滑'（A型与B型）组装。
• 在围压作用下模拟机械平衡状态，随后通过接触网络计算刚度矩阵，以推导弹性模量。
• 采用Voigt-Reuss-Hill平均方法，从刚度矩阵计算有效体积模量（B）和剪切模量（G）。
• 利用标准弹性理论公式，从模量和质量密度推导纵向波速（VP）和横向波速（VS）：VP = √[(B + 4G/3)/ρ]，VS = √[G/ρ]。
• 将模拟得到的波速与不同围压（30–1000 kPa）下干燥（E1）和润滑（E2）玻璃微珠堆集的超声测量实验数据进行对比。
• 分析模拟样品中的配位数和接触网络结构，以关联微观结构与宏观弹性响应。

实验结果

研究问题

• RQ1 为何润滑颗粒堆集的声波速度高于干燥堆集，尽管其堆积密度较低？
• RQ2 样品制备方法——特别是润滑——在多大程度上影响致密颗粒堆集中配位数和弹性模量？
• RQ3 能否通过简化数值模拟实验室制备过程，再现颗粒材料中实验测得的波速趋势？
• RQ4 在颗粒堆集中，配位数是否比密度更具影响力，从而决定弹性波速度？
• RQ5 在不同围压下，数值模拟堆集的力学响应与实验结果相比如何？

主要发现

• 润滑堆集（E2）的纵向波速和横向波速比干燥堆集（E1）高出10–20%，尽管其固相分数较低（0.62 vs. ~0.64）。
• 数值模拟显示，C型样品（模拟振动干燥堆积）的波速与实验测得的干燥堆集（E1）结果匹配度远高于A型（理想润滑）或B型（非理想润滑）模拟结果。
• 润滑堆集中波速更高的原因在于其配位数更高，从而增强了刚度，其影响超过密度的影响。
• 模拟结果表明，在致密颗粒堆集中，配位数比堆积密度对弹性波速度的影响更为显著。
• 本研究证明，数值模型必须复现制备工艺，而不仅仅是最终密度，才能准确预测力学行为，因为配位数是独立的微观结构变量。
• 在低围压（~100 kPa）下，高、低配位数堆集之间的差异至关重要；而在高围压（>1 MPa）下，不同模拟类型之间的差异减小，与实验趋势一致。