[论文解读] Contact and Quasi-Static Impact of a Dissipationless Mechanical Model

Hiroto Kuninaka, Hisao Hayakawa|arXiv (Cornell University)|Apr 29, 2005

Adhesion, Friction, and Surface Interactions参考文献 1被引用 3

一句话总结

本文提出了一种无耗散的二维弹簧-质量模型，引入缺陷粒子和自由边界条件，以在平衡状态下再现赫兹接触力学，并在低速准静态撞击中再现相关行为。通过消除显式的能量耗散，该模型成功再现了赫兹力-形变关系以及准静态理论对恢复系数与撞击速度关系的预测，表明不可逆性可从可逆的微观动力学中产生，而无需外部阻尼。

ABSTRACT

Collisions and contacts of elastic materials are numerically and theoretically investigated. Using a two-dimensional spring-mass model with defect particles under the free boundary condition, we reproduce the Hertzian contact theory at equilibrium and the quasi-static theory for low speed impacts.

研究动机与目标

• 从纯粹的机械性、可逆的微观模型直接推导宏观弹性定律——赫兹接触与准静态撞击定律。
• 解决先前无耗散模型在接触力作用下无法达到平衡的问题。
• 在不引入显式阻尼的情况下，使数值模拟与理论预测在低速撞击下保持一致。
• 探究有限自由度系统中不可逆性的起源是否可由可逆动力学产生。

提出的方法

• 采用包含1099个质点的二维质量-弹簧晶格模型（圆盘）和1269个质点的墙模型，使用非线性弹簧，其刚度 $ k_a $ 和 $ k_b \propto k_a / R^2 $。
• 通过短程指数势 $ F \propto \exp(-a l_s) $ 建模接触力，其中 $ l_s $ 为到最近墙弹簧的距离，且 $ a = 500/R $。
• 采用四阶辛积分器进行系统演化，时间步长 $ dt = 10^{-3} R/c $，以保持能量守恒与时间反演对称性。
• 引入缺陷粒子以破坏完美周期性，抑制能量的周期性交换，从而实现向平衡态的弛豫。
• 模型采用自由边界条件，避免固定约束引发的人工振荡。
• 通过各向同性压缩与剪切的微扰方法计算弹性模量（杨氏模量 $ E $，泊松比 $ \nu $），得到 $ E = 0.773k_a $，$ \nu = 0.336 $。

实验结果

研究问题

• RQ1 纯粹机械性、无耗散的模型能否在平衡状态下再现弹性体之间的赫兹接触定律？
• RQ2 准静态理论是否能准确预测可逆系统中低速撞击的恢复系数？
• RQ3 为何先前的无耗散模型在接触力作用下无法达到平衡？如何修正此问题？
• RQ4 宏观接触与撞击现象中的不可逆性是否可从可逆的微观动力学中产生？
• RQ5 缺陷粒子与自由边界条件在实现此类模型中弛豫与平衡的过程中起到何种作用？

主要发现

• 该模型在平衡状态下成功再现赫兹接触关系 $ \delta \propto \frac{P}{\pi E^*} \left[ \ln\left( \frac{4\pi E^* R}{P} \right) - 1 - \nu \right] $，证实了无非人工耗散。
• 对于低速撞击（$ v/c \leq 0.007 $），模拟得到的恢复系数 $ e $ 与准静态理论预测高度一致，验证了该理论在该参数区间的适用性。
• 在较高速度下（$ v/c > 0.007 $），模拟与准静态理论之间的偏差源于内部振动模态的激发，导致准静态假设失效。
• 从数据中提取的特征时间 $ \tau_0 = 0.011 R/c $ 与准静态框架中的理论估计一致。
• 计算得到的弹性模量为 $ E = 0.773k_a $ 与 $ \nu = 0.336 $，与晶格的预期力学行为一致。
• 初步结果表明，重力会影响恢复系数，提示其在准静态撞击动力学中具有非平凡作用。