M039 · 离散力学

颗粒冲击链与缓冲设计

给首颗粒固定初始动能，跟踪仅受压 Hertz 接触在均匀、渐变、预压和高刚度链中的传播。沿轨道放置力探针，测量阈值到达与墙面反射脉冲，并在总质量和最大压缩预算下比较 25 个质量渐变方案。

已采集归一时间 $0$

总机械能 $0$

动能 $0$

键势能 $0$

有向能量误差 $0$

链动量 $0$

独立积分净墙面冲量 $0$

动量减初值与墙面冲量 $0$

步长减半位移差 · 最初十个时间单位 $0$

实际归一积分步长 $0$

最大积分能量误差 · 完整记录 $0$

记录与模型范围 $0$

采样墙面超额力峰值 · 完整记录 $0$

首次穿越墙面峰值的十分之一 $0$

双探针阈值传播速度 $0$

有向右墙超额冲量 $0$

链归一总质量 $0$

最大采样接触压缩 $0$

已采集质量渐变候选数 $0$

有限设计搜索状态 $0$

最佳可行末端质量比 $0$

释放 / 采集 / 核验

相同冲击能量下，质量渐变能降低墙面峰值吗？

颗粒数

N

24.00\,

渐变末端与首端质量比

m_N/m_1

1.00\,

归一 Hertz 系数

\kappa_H

1.00\,

初始接触压缩

\delta_0

0.00\,

初始动能冲击能量

E_{in}

0.02\,

归一观测时长

\tau_{end}

120.00\,

第一探针沿链位置

s_p/L

33.00\,

设计平均质量上限

\bar m_{max}

1.80\,

设计接触压缩上限

\delta_{max}

0.35\,

最大归一步长

h_{max}

0.02\,

检查已采集记录

t/t_{end}

0.00\,\%

全部量归一；可操作轨道、位移图和键盘滑块。记录在模型边界或步数预算处停止，不裁剪运动。

这是归一一维质点颗粒模型，含弹性 Hertz 接触、固定墙面，不含摩擦、塑性损失或旋转。首颗粒获得初始动能，没有额外驱动撞锤。预压参数是几何压缩而非力。接触力仅压缩或为零，不添加拉伸接触。能量使用预压平衡参考势能，固定墙不做功。探针力是超出基线的正压缩力；右墙冲量积分有向超额力，预压接触卸载时可为负。到达采用首次穿越各探针完整观测峰值十分之一的时刻。峰值与压缩限制使用记录检查点。25 候选搜索保持冲击能量、刚度、预压和时窗相同，排除不完整记录与尚未到达墙面的脉冲。较低有限时窗峰值可能来自延迟或重新分配，不代表吸收，也不认证后续安全或全局最优。接触压缩超过参考间距一半时停止模型。

先抓重点

核心问题: 较低墙面峰值来自延迟、反射，还是能量吸收？
一句话直觉: 保持输入能量相同，追踪有向墙面冲量与链中机械能，同时满足质量和压缩约束。

核心数学模型

仅受压 Hertz 接触

\delta_i=[\delta_0-r_i]_+,\qquad F_i=\kappa_H\delta_i^{3/2}

张开的接触物理接触力为零。

预压平衡参考势能

V_i=\tfrac25\kappa_H(\delta_i^{5/2}-\delta_0^{5/2})+F_0r_i,\qquad F_0=\kappa_H\delta_0^{3/2}

线性参考项核算平衡预压；物理接触仍仅受压。

固定输入能量

E_{in}=\tfrac12m_1v_1^2,\qquad m_i=m_1R^{(i-1)/(N-1)}

每个质量渐变候选都获得相同首颗粒动能。

动量与外部墙面

P-P_0=\int(F_L-F_R)\,dt

墙面冲量独立于动量读数积分。

有限约束设计目标

\min_R\max_{t\in\mathcal T}[F_R(t)-F_0]_+,\qquad \sum_i m_i\le N\bar m_{max},\quad\max_{i,t\in\mathcal T}\delta_i\le\delta_{max}

记录时刻集合与 25 个渐变比值定义有限搜索；要求墙面脉冲已到达。

常见难点

没有拉伸接触

容易误解接触弹簧会把邻近颗粒拉回。

正确理解Hertz 物理压缩仅在零接触处截断；允许接触张开。

降低峰值不等于耗散

容易误解较低墙面力峰值证明吸收冲击能量。

正确理解此模型为弹性；峰值移动或分裂时，能量仍在运动与接触中。

没有到达不能靠延迟取胜

容易误解短记录中的零墙面峰值就是最佳缓冲。

正确理解搜索排除未到达墙面的记录，只报告给定时窗内可行性。

动手实验

01
跟踪压缩脉冲
采集均匀链并检查力探针。
应该观察到： 脉冲先到达探针，再在固定墙面反射。
02
比较质量渐变
在相同冲击能量下选择重末端链。
应该观察到： 传播时间和墙面波形改变，不引入人为能量损失。
03
检查预压
比较预压接触，并直接改变探针位置。
应该观察到： 力读数为超出预压的部分；卸载可产生有向负超额冲量。
04
在双预算下设计
执行候选搜索，再收紧质量或压缩上限。
应该观察到： 可行候选同时满足两约束，并且具有测得的墙面到达。