部分装载罐车中的液体晃动是运输行业的重要安全隐患。制动过程中,液态货物向前移动,产生影响车辆稳定性、增加制动距离并可能导致侧翻事故的动态力。挡板通常安装在罐体内部以抑制这些振荡,但其效果取决于设计方案和机动操作的剧烈程度。仿真技术使工程师能够在实物样机制造之前,以虚拟方式评估挡板配置方案并量化其对晃动力的影响。shonDy 的晃动仿真能力已通过本网站专项案例研究与实验数据进行了验证。
案例描述
本案例研究对满载罐车从 80 km/h 紧急制动至完全停止的工况进行了仿真。仿真采用基于重力的方法,而非在空间中移动罐体几何模型:利用等效原理,将制动减速度表示为重力矢量随时间的变化。在随车(卡车)参考系中,液体受到等效于制动减速度的伪力,这在物理上等价于将重力矢量向前倾斜。
制动过程中的等效重力矢量为:
其中 为随时间变化的制动减速度曲线。行驶方向定义为负 x 方向 ,因此液体在制动时向 方向移动,正 x 方向位移对应向后运动。
制动曲线
制动曲线设计用于模拟满载卡车的气动制动系统,如下图所示:
在最大制动时,等效重力矢量的大小约为 11.77 m/s²(1.20 g),偏离竖直方向约 33°。

几何模型
罐体长 11.5 m,截面高度为 1.6 m。请注意,对于仿真而言,只有罐体和挡板在几何上是相关的。卡车车身仅用于指示方向(行驶方向)并提供整体设置的直观展示。
设置了两种配置进行比较:
- **无挡板:**罐体内部无障碍物,代表晃动最恶劣的工况。
- **有挡板:**内部隔板将罐体分隔成若干隔舱,减小自由晃动长度。
带挡板配置的几何模型如下图所示。

案例设置
卡车装载 10 m³ 密度为 884 kg/m³ 的液体,总流体质量为 8,840 kg。仿真采用半径为 6 mm 的粒子,共约 580 万个粒子。总仿真时长为 15 秒,涵盖完整的制动过程及随后的制动后晃动阶段。
| 参数 | 数值 |
|---|---|
| 流体密度 | 884 kg/m³ |
| 流体体积 | 10 m³ |
| 运动黏度 | 16 cSt |
| 粒子半径 | 6 mm |
| 粒子数量 | ~580 万 |
| 仿真时长 | 15 s |
| 初始速度 | 80 km/h (22.22 m/s) |
| 最大减速度 | 6.5 m/s² (0.66 g) |
结果
总览
下方视频展示了两种罐体配置的制动过程。在 1 秒的静止时间内,流体保持静止。随着制动力逐渐增大,液体涌向前壁。在无挡板的罐体中,全部流体质量以整体波的形式冲击前壁。在有挡板的配置中,每个隔舱容纳的液体体积较小,波高明显降低。
质心位移
对流体质心(CoM)在行驶方向(x)和竖直方向(y)随时间的变化进行了跟踪。位移以静止时的初始位置为参考。
车辆停止后,制动力消失,但储存在流体中的动能不会立即耗散。液体继续在罐体内来回晃动,产生制动后振荡,该振荡持续数秒,即使在车辆速度为零时也会在罐壁上产生交变力。这些振荡的幅值和衰减过程在下图中可直接观察到,表现为制动阶段结束后持续的振荡运动。
在行驶方向上,无挡板时质心最大前向位移达 2,750 mm。有挡板时降至 2,199 mm,减少 20%。有挡板配置中,制动后振荡的衰减速度也明显更快。为量化这一效果,计算了车辆停止后(t > 4.7 s)所有时间步的质心 X 方向位移标准差。标准差越大,表明振荡越强、持续时间越长。无挡板时,该值为 1,185 mm;有挡板时降至 500 mm,减少 58%。
在竖直方向上,无挡板时随着流体攀升前壁,质心最高升高 375 mm。挡板将此升高量限制在 267 mm,减少 29%。较低的竖直质心位移直接降低了制动过程中的有效重心高度,这是影响侧翻稳定性的关键因素。

质心速度
质心位置的一阶导数给出了晃动速度,表征流体质量在罐体内移动的速率。
无挡板时,行驶方向上质心速度峰值达 3.12 m/s。挡板将此峰值降至 1.44 m/s,减少 54%。有挡板配置还表现出制动后速度振荡衰减明显更快的特点,证实了内部隔板的阻尼效果。

晃动力
根据平滑后的质心加速度,以流体质量 8,840 kg 计算动态晃动力 。该力表示晃动流体施加在罐体结构上的额外动态载荷,超出固化(刚性)质量所产生的载荷。
该力的符号变化表示晃动波的反转:当流体从前壁反弹并开始向后运动时,动态晃动力改变方向。这些交变力正是使晃动具有危险性的原因。它们可以激发车辆动力学模态并降低操控稳定性。
无挡板时,行驶方向上的晃动力峰值达 76.5 kN。有挡板时,该峰值降至 33.2 kN,减少 57%。

总结
| 指标 | 无挡板 | 有挡板 | 降低幅度 |
|---|---|---|---|
| 质心最大位移 X [mm] | 2,750 | 2,199 | 20 % |
| 质心最大位移 Y [mm] | 375 | 267 | 29 % |
| 质心最大速度 X [m/s] | 3.12 | 1.44 | 54 % |
| 晃动力峰值 X [kN] | 76.5 | 33.2 | 57 % |
| 制动后振荡幅度 X [mm] | 1,185 | 500 | 58 % |
本案例研究使用 shonDy 的基于重力的方法对满载罐车的紧急制动工况进行了仿真。通过随时间修改重力矢量以表示制动减速度,仿真避免了在空间中移动整个罐体的需求,同时通过等效原理在物理上保持等价。
对两种配置进行了比较:无内部结构的罐体和有挡板的罐体。结果表明,挡板能有效降低晃动力峰值、限制质心位移,并加速制动后振荡的衰减。这些改进直接转化为更高的车辆稳定性和更低的紧急制动侧翻风险。
本文提出的仿真方法为在真实制动工况下评估罐体设计和挡板配置提供了一种实用且高效的手段,使工程师能够在无需实物样机的情况下优化货罐安全性。
渲染对比
下方视频展示了整个制动过程中两种配置的并排渲染对比。上半部分显示有挡板的罐体,下半部分显示无挡板的罐体。主要结果直接叠加显示在视频中。
案例设置
| Fluid density | 884 kg/m³ |
| Fluid volume | 10 m³ |
| Kinematic viscosity | 16 cSt |
| Particle radius | 6 mm |
| Particle count | ~5.8 million |
| Start velocity | 80 km/h |
| Max. deceleration | 6.5 m/s² (0.66 g) |



