罐车制动工况液体晃动流固耦合分析.docxVIP

罐车制动工况液体晃动的流固耦合分析 - 汽车 罐车制动工况液体晃动的流固耦合分析 张韬 柯龙燕 白光建 中航（上海）汽车技术有限公司 上海 200092 摘要：采用 VOF 模型对罐车制动过程中内部液体的自由晃动过程进行 了模拟，计算得出罐体各个部分的应力分布。 计算结果表明， 制动时前封头和前 四个防波板受力逐渐上升， 在达到峰值后迅速下降； 后封头和第五个防渡板受力 则呈逐渐减小的趋势。 随着与后封头距离的增大， 防波板承受最大应力也相应增 大，单个防波板承受最大应力位置处于底部孔边缘处。 关键词 ：罐车 两相流 液体晃动 流固耦合 中图分类号：U469.6+1.02 文献标识码：A  文章编号： 1004-0226(2015)08-0088-04 前言 罐车在行驶过程中由于加速、减速或转弯等运动状态的改变，会产生惯 性力，进而引起液体的晃动， 出现罐车罐体与介质的流固耦台效应。 在这种流固 耦台系统的作用下， 一方面可能会导致罐车的动力不稳定， 另一方面罐体中介质 的晃动会对罐体产生显著的压力， 严重时会造成罐体结构破坏。 罐体内一般都布 置有防波板或阻浪板， 将罐体内的液体分成多个小舱室， 在流固耦合作用下， 整 个罐体内液体的冲击力将由各防波板分担，因此需对防波板的受力过程进行分 析，以便通过在罐体内部设置合理的防晃装置降低罐车防波板损坏的可能性。 本文使用流固耦合分析软件对罐车在紧急制动工况下罐内的油料晃动 过程进行了分析计算， 研究了罐体各部分受力随时间的变化趋势以及防波板的应 力分布情况。 计算模型和计算方法 2.1 控制方程组 本文对于流体部分采用 VOF 多相流模型计算罐内气液两相非定常湍流 流动，控制方程组包括连续性方程和动量守恒方程。 2.2 罐体模型和网格划分 以某型号油罐车为研究对象， 其罐体模型如图 1 所示。罐体由前后封头、 个防波板和罐体组成。其中防波板将罐体分割为六个区域，区域之间通过两个开孔相连。罐体总长约 10.2 m ，宽约 2 3m ，高约 1.6 m 。 计算时，坐标系固定于罐体上， 原点位于罐体前封头端面底部中心位置， 砖由以罐体前进方向为正方向且平行于地面， y 轴垂直于地面，向上为正方向， 轴平行于地面，向右为正方向。模型网格划分采用混合网格结构，使用移动坐 标平面将罐体切分为 6 个小罐体和 5 个包含防波板的部分，小罐体可采用六面 体扫掠网格，其余部分采用四面体网格， 网格过渡部分保证平滑且节点一一对应， 如图 2 所示。这样可以大大减少网格总数，提高计算效率。 2.3 边界条件 罐体及防波板表面设置为无滑移壁面边界条件。多相流采用 VOF 模型， 可以精确模拟气液分界面的变化情况以模拟液体晃动过程。 初始气液交界面平行 于 xz 平面，液体和气体相对罐体速度为 0 ，气体压强为 1.013xlO 5 Pa 。湍流模型采用标准 k-a 模型和标准壁面函数模型。 罐车以 80 km/h 的初始速度正向行驶，在 O.l s 时施加 0.5 g 的反向加 速度模拟制动工况。罐车内介质为油，密度 900 kg/m3 ，粘性 0.048 kgm ．s， 充装率 80% 。实体材料均为 Q235 钢，弹性模量 2.06xlOl Pa ，泊松比为 0.250 制动速度曲线如图 3 所示。 2.4 求解器参数设置 速度和压强的耦合处理采用 PISO 方法，压力修正方程离散格式采用 Body Force Weighted 格式，对流项离散格式采用二阶迎风格式。非定常计算 中时间步长为 0.002 s 。 制动工况计算结果及分析 图 4 、5 为前后封头和各防波板压力随时间的变化曲线。从各曲线图可以看出，罐车制动后在 0.52 s 时液体各部分压力达到最大值，而后开始逐渐减小。 在制动之前，罐体各部分由于液体重力的存在而受到的压力基本相同。 在制动初期，前封头和前 4 个防波板所受压力呈现震荡式上升。 这是由于初期液 体晃动趋势较小， 在经过四次震荡后逐渐向峰值接近。 后封头和第五个防波板受 力曲线相似， 均呈现逐渐减小的趋势， 同时越靠近后封头的防波板受力越小， 这 是由于液体正向运动的惯性导致罐体前后部的受力趋势不同。 罐体受力为峰值时，液体晃动最为剧烈，随后受力减小。这是由于液体 受到罐体一侧壁面的阻挡开始做折返运动， 涌向另一侧。在整个晃动过程中前封 头受到的峰值压力最大。 图 6 、7 为初始时刻、 0.52 s 即压力达到峰值时罐体内液相分布图。图 8-11 为上述两个时刻罐体和防波板压力分布图。压力分布情况与前述分析结果 一致。罐内液体受晃动影响集中在每个腔内的前部， 因此每个腔内前部压力大于 后部压力。 0.52 s 时液体处在向前运动的动量最大处，