摘要:由于机油泵内部流体流动的复杂性,进行数值计算比较困难。本文对摆线转子式机油泵流场进行仿真分析,得出机油泵流量和压力波动曲线等。并对机油泵进行了空化现象预测,取得比较理想的效果。仿真结果与理论分析比较验证仿真结果的正确性,为机油泵设计提供一个指导。
关键词:转子泵 Pumplinx 流量脉动 空化
前言
机油泵是发动机润滑系统的心脏,直接影响着发动机的整机性能,尤其是发动机的使用寿命和节能效果[1]。转子式机油泵结构紧凑,外形尺寸小,重量轻,吸油真空度较大,泵油量大,供油均匀性好,成本低,在中小型发动机上应用广泛。转子式机油泵由壳体、内转子、外转子和泵盖等组成。内转子用键或销子固定在转子轴上,由曲轴齿轮直接或间接驱动,外转子由内转子带动在壳体孔中旋转。转子泵的运动方式和内啮合齿轮副一样。内转子A和外转子B中心的偏心距为e,内转子带动外转子一起沿同一方向转动。内转子比外转子少一个齿,这样就保证了内外转子同向不同步地旋转[2]。如图1所示。
机油泵流体仿真分析比较复杂,由于机油泵腔体内的液体在高速运转时体积不断变化,所以要采用动网格划分,对计算机的性能有较高的要求,并且难以保证精度。通用CFD软件对技术人员的要求较高,缺少经验的技术人员往往不能得到理想的结果。Pumplinx软件很好的解决了这些问题。
图1
1计算模型[3]
流体流动所遵循的物理定律,是建立流体运动方程的依据。流体满足质量守恒定律、动量守恒、能量守恒等,另外需要考虑不同的流态等。
(1)
(2)
式中,表示控制体,和分别表示控制体表面以及表面的外部矢量。是液体的密度。和表示液体的流速和表面流速。P表示压力,表示体力。剪应力系数有粘度和速度决定。可以表示为
(3)
(i=1,2,3)是速度在三个坐标上的分量。表示为
(4)
Pumplinx空化模型是基于两相流动思想,采用Rayleigh Plesset方程求解气泡变化动态过程[4]。
(5)
其中是气体分布率,为湍流施密特数。和是气体生成率和凝结率。可以通过下式计算。
(6)
(7)
其中,取0.02,取0.01。
2转子式机油泵流场数值模拟
2.1转子式机油泵的流场建模
根据转子泵的参数建立流场的三维流场模型,本文机油泵考虑侧盖间隙的泄漏情况。采用UG软件进行机油泵的物理建模,并通过布尔运算得到机油泵的流场部分的三维模型。在Pumplinx中导入泵的STL文件,进行网格划分。采用自动网格生成器,专有的几何等角自适应二元树算法,既CAB算法。CAB算法在由封闭表面构成的体域生成迪卡尔六面体网格。在靠近几何边界,CAB自动调整网格来适应几何曲面和几何边界线。为了适应关键性的几何特征,CAB 通过不断的分裂网格来自动的调整网格大小,这是利用最小的网格分辨细节特征的最有效方法。另外,泵转子的网格采用自动结构化网格生成器生成。
采用同位网格下的有限体积法求解瞬态三维的粘性流体力学方程。CFD算法采用了适合低速液体的、以速度和压力为基本变量的SIMPLEC方法。采用准确直观的动网格方法来处理泵的运动,初始时对动静流道分别建摸,在计算过程中,动流道部分自动处理网格的移动和变形。在每个时间步,动静区域通过全隐式滑移界面技术MGI自动在新位置重新联结起来[5]。将摆线泵的面几何分多个面。摆线泵几何中的外齿轮的内侧面、内齿轮的外侧面、内、外齿轮的上、下底壳面需为单独的面。根据面进行网格划分,首先对流体进行网格划分,流体网格是动网格,是CFD分析的关键。出入口则采用二叉树网格法划分。图2为流场模型,图3为网格划分后的模型。
图2 实体模型
图3 网格划分
2.2添加边界条件求解
根据机油泵实际工作状况,添加边界条件,主要有机油泵入口压力,油液基本参数(油液密度、粘度、油温等)。参数值如表1、2所示。在模型中添加空化和湍流分析模型,并在出口处添加监测点。
表1油液基本参数
油液参数
数值
油密度
动力粘度
液体弹性模量
含气体密度
饱和压
表2 边界条件
边界
值
进口压力
出口压力
温度
空气含量(体积)
(入口)
3.结果分析
3.1流量和出口压力脉动
图4给出了机油泵在1000RPM、4000RPM和6000RPM的瞬时流量图,从图中可以看出,当机油泵运转稳定后,流量会稳定在一个平均值,并在一个平均值附近存在一定的波动。并且,随着转速的增加,流量波动愈大。根据流量计算公式
其中,内转子齿每转扫过的面积,n为机油泵转速。计算得到不同转速下机油泵的理论流量。与仿真结果进行对比,如表3
表3 理论流量与数值计算结果比较
数
值
转
速
理论值
(L/min)
仿真值
(L/min)
1000RPM
18.9
17.66
4000RPM
75.6
69.78
6000RPM
113.4
92.76
图4 不同转速下流量变化图
从表中看出,仿真结果比理论计算稍低,这是因为,实际工作中机油泵存在径向和轴向间隙,存在着一定的泄露。所以从流量曲线变化趋势看出,转速越高,模拟值和理论值相对误差越大,因为转速越高油液压力越大,泄漏量也就越大。
流量脉动可以有流量脉动率评价,其公式为
(8)
式中,为齿轮泵的转速,为小齿轮泵的齿数。从式中看出,流量脉动率随转速的增加、齿数的增多而增大。本文,通过仿真分析了同一机油泵,即齿数一定不同转速情况下的流量脉动情况,。从图5中可以看出,随着转速的增加,流量波动越大。与理论分析完全相符。
图5瞬时流量波动
图6 进口凸台瞬时压力波动
图7 出口凸台处瞬时压力波动
流量脉动伴随相应的压力波动,从而产生噪声振动。由进出口监测点可以得到进口和出口的压力波动曲线,从图中看出,压力稳定后,在转速较高时出口压力波动较大。在进口凸台,转速增加监测点处压力平均值下降,因为转速高进油腔真空度增加。图7出口压力波动看出,转速增加压力波动幅值增加明显。如何降低内啮合摆线转子泵工作时的脉动以及由此产生的噪声和振动,一直是设计者追求的目标。杨元模等[6]采用了在出油口添加一个稳压阀,通过一个较大弹力范围的弹簧调整,达到了稳压效果。并通过实验验证,增加稳压阀后,泵的脉动频幅不及原来的1/3,从而大为降低由此产生的噪声振动。
3.2空化现象分析