运动物体在颗粒物质中的动力学过程及最大穿透深度仿真研究.pdfVIP

物理学报 Acta Phys. Sin. Vol. 62, No. 13 (2013) 130204 运动物体在颗粒物质中的动力学过程及最大 穿透深度仿真研究* † 何克晶 张金成 周晓强 ( 华南理工大学计算机科学与工程学院, 广州 510641 ) ( 2012 年12 月24 日收到; 2013 年3 月3 日收到修改稿) 利用离散元法仿真了运动物体在颗粒物质中的三维动力学过程, 仿真采用周期边界条件, 并考虑了重力、接触 力、阻尼力、摩擦力的影响. 将仿真结果和相关的三维实验结果进行了对比, 两者符合较好. 仿真结果表明穿透深度 与运动物体的冲击速度、运动物体质量、颗粒介质床的密度均有关系. 运动物体质量越大, 速度越快, 则穿透越深, 而且穿透深度和质量呈线性关系. 仿真过程较为真实地再现了小颗粒的飞溅现象. 关键词: 颗粒物质, 动力学过程, 仿真, 离散元法 PACS: 02.70.Ns, 05.10.a, 45.70.n, 81.05.Rm DOI: 10.7498/aps.62.130204 分结构变化. 传统的实验手段和方法难以观测到 1 引言 颗粒物质内部的结构变化以及各个颗粒的动力学 过程. 而采用计算机仿真让研究者较为容易观测 颗粒物质普遍存在于自然界、工业生产和日 物质内部的结构变化并跟踪每个单独颗粒的动力 常生活中, 对自然界和生活生产起着重要的影响. 学过程. 而运动物体对颗粒物质的冲击也广泛存在, 比如 因为颗粒介质固有的离散性, 以及冲击时造成 陨石撞击, 沙尘暴, 弹丸入射等. 冲击现象的普遍 的大形变和结构破坏, 传统的连续介质力学模型和 存在及重大意义, 引起了众多研究者的广泛兴趣. 方法(如有限元方法) 已难以应用. 在离散元方法7 Ciamarra 等对运动物体在准二维颗粒床中动力学 中, 系统由许多独立的颗粒组成. 这些颗粒的运动 1 过程进行了实验研究 , 实验发现阻力正比于冲击 和位移都被单独模拟. 离散元方法中各离散元是独 初速度. Lohse 等研究了运动物体在三维颗粒床中 23 立的, 并且用软接触的方法来建模接触行为, 符合 的动力学过程 , 实验表明阻力与阻力系数及物 颗粒介质的物理本质, 非常适合用于仿真颗粒物质 体进入的深度相关. Khaldoun 等研究了在河流的流 89 的动力学过程 . 沙等湿颗粒体系中, 液体黏滞性的影响以及下沉深 4 刘凯欣等提出了三维连结型离散模型, 并模拟 度与物体的剪切应力之间的关系 . 张权义、Hou 等通过高速摄像的跟踪, 研究了重力场中球形下落 了连续介质转变为非连续介质的力学过程, 其结果 物体对松散颗粒床