M039 · 离散力学

颗粒冲击链与缓冲设计

冲击能量相同，让珠子质量逐个渐变，能降低墙受到的峰值力吗？ 在观察的这段时间内能。珠子质量逐个渐变，传到墙上的脉冲被拉宽、推迟，峰值就低了。但弹性链不吸收能量，只是把它重新分配，后面的反射照样可能把能量送到墙上。设计时还得守住总质量和压缩量的上限。

给第一颗颗粒一份固定的动能，看脉冲怎样穿过均匀链、渐变链、预压链和高刚度链：三维画面里，颗粒沿链方向运动，Hertz 接触只推不拉。在轨道上放一个力探针，测出脉冲越过阈值的到达时刻和墙面反射的脉冲；再在总质量和最大压缩量的上限内，比较 25 个质量渐变方案。

左侧真实接触力 $0$

右侧真实接触力 $0$

检查颗粒加速度 $0$

检查右侧接触 $0$

已采集归一时间 $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 压缩量不等于标记表面之间的距离，模型也不画颗粒变形后的轮廓。

先抓重点

核心问题: 墙面受力峰值变低，是因为脉冲来晚了、被反射了，还是能量被吸收了？
一句话直觉: 让每个方案的输入能量相同，跟踪带符号的墙面冲量和链里剩下的机械能，同时让总质量和压缩量都不超限。

核心数学模型

仅受压 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
在两个上限内设计
执行候选搜索，再收紧质量或压缩上限。
你会看到： 合格的方案要同时满足两个上限，而且确实测到脉冲到达墙面。