[论文解读] Sub-surface granular dynamics in the context of oblique, low-velocity impacts into angular granular media

Peter Miklavcic, Hesam Askari+3 more|arXiv (Cornell University)|Jan 7, 2022

Planetary Science and Exploration参考文献 64被引用 2

一句话总结

本研究采用二维离散有限元模拟，对倾斜、低速撞击角状颗粒介质的过程进行建模，再现了弹跳、滚动和完全停止等实验行为。研究发现，提高撞击速度会增强横向能量耗散，而非导致更深的次表层形变，提出了‘表层区’（skin zone）这一永久应变区域的概念，位于表面之下。

ABSTRACT

Oblique, low-velocity impacts onto extraterrestrial terrain are an inevitable occurrence during space exploration. We conduct two-dimensional discrete simulations to model such impacts into a bed of triangular grains. Finite element method provides the basis for simulation, enabling the angular grain geometry. Our findings re-create the three classes of impact behavior previously noted from experiments: full-stop, rollout, and ricochet \citep*{Wright2020}. An application of Set Voronoi tessellation assesses packing fraction at a high resolution, revealing how grains shift relative to each other during an impact event. Calculation of Von Mises strain distributions then reveal how grains shift relative to the overall system, leading to the notion of the 'skin zone'. Intuition would suggest that the region of perturbed grains would grow deeper with higher velocity impacts, results instead show that increasing velocity may actually evoke a change in the grains' dissipative response that boosts lateral perturbation. Finally, we consider as a whole how sub-surface response could link with impactor dynamics to deepen our understanding of oblique, low-velocity impact events and help to improve mission outcomes.

研究动机与目标

• 模拟倾斜、低速撞击非球形、角状颗粒介质时的次表层颗粒动力学行为。
• 通过离散模拟再现实验观测到的撞击行为——完全停止、滚动和弹跳。
• 利用高分辨率Voronoi剖分量化局部颗粒堆积密度的变化及其在撞击过程中的演化。
• 识别并测量‘表层区’——即超过准静态屈服阈值的永久应变区域——的深度。
• 研究撞击速度和角度如何影响颗粒介质的能量耗散机制及次表层扰动。

提出的方法

• 采用离散有限元法（DFEM）对三角形颗粒进行建模，实现对角状、非球形颗粒几何形状的精确表示。
• 有限元公式支持不规则形状颗粒之间的接触检测与响应，实现真实的机械相互作用。
• 应用固定Voronoi剖分以高分辨率计算局部堆积密度，追踪撞击过程中颗粒的重排。
• 计算Von Mises应变分布，以识别永久形变区域并定义表面之下的‘表层区’。
• 在不同撞击速度（2–10 m/s）和角度（30°, 60°）下进行模拟，模拟在类地重力条件下，与实验数据进行对比。
• 未来可扩展为可视化节点应力场，以映射力网络在颗粒介质中的传播。

实验结果

研究问题

• RQ1 在倾斜、低速撞击角状颗粒介质时，次表层颗粒响应如何随撞击速度增加而变化？
• RQ2 ‘表层区’——即超过准静态屈服阈值的永久应变区域——的空间分布与深度如何？
• RQ3 局部堆积密度的变化与撞击行为类型（弹跳、滚动、完全停止）之间存在何种相关性？
• RQ4 提高撞击速度是否导致颗粒介质中更深或更横向的能量耗散？
• RQ5 撞击角度与速度在多大程度上调节次表层扰动的范围与形态？

主要发现

• 模拟成功再现了在不同撞击角度和速度下实验观测到的三种撞击行为类型——弹跳、滚动和完全停止。
• 在表面之下形成一个不对称的扰动区域，靠近表面的基底局部孔隙率较高，撞击点前方则发生压实，其范围约等于撞击坑半径。
• ‘表层区’——即超过准静态屈服阈值的永久应变区域——被量化并发现其深度随撞击速度增加而增大。
• 与直觉相反，更高速度的撞击并未使扰动区域更深，反而促进了横向能量耗散机制。
• 在更陡的角度（如60°）下，瞬时撞击坑深度（TCD）和平均表层深度（MSD）随速度的增加比在较平缓角度下更为显著。
• 节点应力场可视化揭示了力网络在颗粒介质中的传播，提示未来可进一步分析月壤中应力传输机制。