Creo 液体有限元分析：矩形罐体三阶液体晃动模态解析

项目过程

本项目委托方为某专用车厂的制动液储液罐配套厂家，需求是在运输颠簸工况下评估矩形罐内液体晃动响应。客户给出的原始尺寸是罐长 L=1.2 m、宽度 0.4 m、初始水深 h=0.5 m，对应液深比 h/L≈0.42。计算目标是提取前 3 阶液体晃动固有频率，并在 1 阶模态下输出自由液面波形与罐壁上的无量纲动压分布 p′/(ρgA)，以验证后续是否需要加装防晃隔板。

本项目使用 Creo Simulate 的 Fluid Modes 模块。罐体取 1/4 对称模型，沿长边方向（x 轴）施加对称边界条件，自由液面采用线性化的 Bernoulli 边界（基于小波假设），液体视为不可压缩势流，密度取 1000 kg/m³，重力加速度 9.81 m/s²。求解器为 Lanczos 法，提取前三阶非零频率，归一化幅值 A 取水深 h 以保证动压无量纲。

案例分析

计算结果与图中曲线吻合良好：当水深从 0.1 m 增加到 0.5 m 时，1 阶晃动频率 f₁ 从 0.36 Hz 单调上升至 0.75 Hz（即图中星标工作点），2 阶由 0.71 Hz 升到 1.15 Hz，3 阶由 1.10 Hz 升到 1.40 Hz，三条曲线在 h>0.5 m 后趋于水平，符合势流理论中 h/L 较大时频率饱和的特征。客户实际工作水深恰好是 0.5 m，对应 1 阶 0.75 Hz。

在自由液面波形图上，1 阶模态呈现明显的 1/2 波长驻波：罐两端 z=0.5 m 出现波峰与波谷的最大幅值（左正右负），x=0.6 m 节点处幅值为零，节点位置与无量纲坐标 x/L=0.5 完全一致。沿罐长方向的无量纲动压云图也呈现典型的反对称分布：左半罐 p′/(ρgA)≈+1.0（红），右半罐≈−1.0（蓝），节点处动压为零。动压幅值与自由表面波高同相位，最大动压出现在两端端壁中部。

总结教训

本项目最容易踩的坑是忘记使用「线性化」的自由液面边界。Creo 默认的流体模态算法基于小波假设，如果实际工况振幅超过水深 5%~10%，频率会被高估 5%~8%。本项目 h=0.5 m、激励幅值约 0.04 m（远小于水深），属于小波假设适用范围，结果可信。

另一个教训是网格密度对动压分布的影响：在 x 方向节点附近加密后，节点位置向中央偏了约 0.3%，原因是节点处动压梯度大，不足网格会”磨平”反相过渡。建议在该区域至少布置 10 个单元。

项目经验

经验一：当客户只关心 1 阶频率时，不必勉强算 3 阶；如果只是核对运输颠簸谱（典型 0.52 Hz）与 1 阶频率错开 30% 以上即可视为安全，1 阶 0.75 Hz 对应避频区间 0.530.98 Hz，远高于一般卡车悬挂激励频率（<2 Hz 的低频峰），实际不必加防晃板；如果客户工况包含船用或越野（0.2~0.4 Hz 低频大位移），则需重审。

经验二：动压云图可用于评估罐壁局部强度。本项目最大动压出现在端壁中部，对应 σ_hoop≈p′·R/t，在罐壁 t=3 mm 时约 1.6 MPa，远低于铝合金屈服，可忽略；若改为 0.5 mm 薄壁罐，则需在端壁中部加筋。

经验三：液深比 h/L 是设计防晃板高度的第一指标。本项目 h/L=0.42，加 1 道 0.15 m 防晃板可将 1 阶频率推高约 30%，对颠簸激励更友好。

有限元仿真