低速刚性体冲击密实散体介质的动力学实验和理论研究.pdfVIP

低速刚性体冲击密实散体介质的动力学 实验与理论研究 庞 勇1 刘才山1 （1 湍流与复杂系统国家重点实验室,北京大学,北京,100871 ） 摘要: 干沙颗粒所组成的散体介质是一类具有典型多尺度耦合效应的物理对象。如果将颗粒等效为理想质点，从经典力学的角度 理解，该类系统可表示为一个包含无穷多质点的复杂多体系统。该类系统在受到外部冲击作用，由于散体颗粒间之间受到典型 的Coulomb 摩擦效应，所反映出来的宏观性质在不同的物理阶段具有明显不同的宏观物理特征。散体性质会发生从固态到流动 状态，以及从流动状态到固态的转化过程。研究刚体冲击作用下的散体介质的动力学行为已成为当前多个学科所关注的重要科 学问题。 揭示这类复杂系统动力学过程的一个有效手段是基于实验现象发现和归纳相关的物理规律。我们首先建立了具有3米高程的 跌落台架，通过不同高度释放具有不同直径的刚性圆柱体，利用加速度计和高速摄像装置采集圆柱体冲击干沙界面时的动态加 速度和位移信息。通过大量的实验数据我们发现，刚性圆柱体在与干沙颗粒冲击作用的过程中具有典型的三个特征阶段：颗粒 流化阶段；流动加速阶段；流动减速至固化阶段。实验数据同时表明：颗粒流化阶段存在一个稳定的临近屈服强度，并且该临 界屈服强度满足土壤力学塑性滑移线理论模式。在流动加速阶段中，颗粒之间的Coulomb摩擦作用将呈现与速度相关的粘性力特 征，并且具有一个稳定的常数值。而在流动减速阶段，该粘性力将与刚体运动的速度以及静水压力密切相关。 在以上实验现象观察的基础上，我们提出了描述圆柱型刚体与散体介质不同冲击阶段的有效的简化模型。在颗粒流化阶 段，将冲击过程等效为固体介质中的冲击波效应，并根据塑性滑移线理论确定颗粒介质流化时的静水压力。在散体流动过程中， 将散体介质与刚性圆柱之间的相互作用等效为与压人深度相关的静水压力项和与速度相关的粘性阻力项。同时，该粘性阻力项 与时变的特征长度相关。并且发现，特征长度的变化规律受到一个无量纲惯性常数的控制。在以上模型的基础上，数值模拟结 果能够有效的和相关实验结果吻合。同时，该模型能够有效的描述于球壳刚体与散体介质的冲击行为，并与相关的实验结果相 吻合。 本文的研究成果对深刻理解受冲击作用下的散体介质性质具有积极的作用。如何建立散体介质流动过程中的粘性力特征 与Coulomb摩擦之间的关联效应将是我们今后的重要研究内容。 关键词: 刚体;散体;冲击;模型 1) **基金(和** 基金 （9876543 ）资助项目 2) E-mail ：py82@163.com liucs@ 引 言 散体的动力学行为在自然界中广泛存在，从大沙漠中的沙子流动到雪山上的雪崩，从工 业上各种颗粒产品的生产到农业中粮食的输送。散体的动力学行为与我们的生活息息相关。 各种散体根据单位体积内的密度比上单位体积内散体颗粒以最密实的方式排列的密度比，可 以分成密实散体和松散散体。同样散体根据其流动速度可以分成快速流动散体和低速流动散 体，还有不运动散体。各种情况下的散体表现出不同的性质。静止不动的散体在外力的作用 下没有发生流动时表现出来的性质和固体类似。流动以后的散体表现出来的性质又象流体。 散体的这些性质，很难找到一种具体模型来描述。但实际生活中与散体相关的事情很多，人 们迫切需要找到其物理规律。各个领域的科学家从各自的领域出发，找出了各种方法对本领 域内的散体的各种行为进行了研究,其中一个重要方面是刚体与散体相互作用时的动力学研 究。由于散体的动力学行为还没有成熟的理论，因此大家采用的研究方法主要是分子动力学 模拟和试验方法。分子动力学方法理论简单，但计算量大，其不能完全代替实验，但试验比 较难实现的时候常用分子动力学模拟。现在这方面的研究大部分采用的试验方法。刚体冲击 散体的动力学过程很早就有人研究。Allen 在1957[1]年对刚体高速冲击散体的动力学行为进 行了总结，他把刚体在冲击过程中的阻力分成速度相关项和常数项，与速度相关项又分成速 度的一次项和速度平方项。提出在高速的情况下刚体受到的阻力主要是速度的二次项起作 用，在低速下主要是速度的一次项起作用，高速和低速的分界点在声速。最近国际上的研究 主要是对刚体低速冲击散体的动力学研究。大家仍然把刚体受到的力分成速度相关项和速度 无关项，但速度无关项与位移相关。John R. de Bruyn[2]认为与速度的相关项为速度