1. 引言
长轴液下泵(以下简称液下泵)广泛用于自来水公司、污水处理厂、电厂、钢铁厂、矿山等工矿企业,适用于输送带颗粒、高粘度、强酸、碱、盐、强氧化剂等多种腐蚀性介质。其转动轴为悬臂式,最长可达十几米,结构简单,占地面积小 [1]。
由于液下泵为细长轴,在运行过程中细长轴容易弯曲变形和扭曲变形等,经常影响用户的正常使用。因此,在设计、研发阶段准确计算泵轴的强度与刚度对提升产品的可靠性有重要的现实意义。国内有学者对此展开了一些研究,夏斌 [2] 基于转子动力学对液下硫磺泵进行研究,表明在加载了不平衡量的谐响应分析中,轴承与叶轮都在前四阶固有频率下产生共振现象。熊珍兵 [3] 采用有限元法对立式液下泵静强度和刚度进行了分析,发现设计选型上存在一点的缺陷。崔常鹏 [4] 采用模态叠加法对液下泵的固有频率和各阶振型进行了分析,得到轴系磨损严重的原因是轴和轴承之间的间隙过大导致的。本文针对某型液下泵为研究对象,采用ANSYS WORKBENCH软件对其轴系统进行强度、刚度及模态分析,给出该型泵设计过程中的合理优化思路。
2. 液下泵模型
液下泵65~310型为单级单吸悬臂式立式离心结构,符合《API610-2010第11版石油石化及天然气工业用离心泵》VS5型,VS5刚性轴设计(一般轴都很粗)不带滑动轴承支撑点,也不要冲洗,液下深度不会太长(一般不超过1.5米),双轴承支撑,两轴承之间的距离长,轴肩少;65为排出口直径,310为叶轮外径;液下深度:740~1340 mm,流量:30~160 m3/h,扬程10~70 m,转速960~2050 r/min,结构和实物如图1所示。
Figure 1. Structural and physical drawings
图1. 结构和实物图
3. 液下泵强度分析
3.1. 轴向力
泵运转过程中转子上的轴向力使得转子轴向移动,该力有盖板力(指向叶轮吸入口方向)、动反力(指向叶轮后面)、轴台等引起的轴向 [5] (方向视情况而定)和转子重量引起的轴向力等组成。为了研究方便,轴向力简化公式为:
(1)
式中:k为系数,当比转速ns = 30~100时,k = 0.6;当比转速ns = 100~200时,k = 0.7,当比转速ns = 240~280时,k = 0.8。ρ为流体密度。H为扬程(m)。Rm为叶轮密封环半径(m)。Rh为叶轮有效轮毂半径(m)。有上式计算得到叶轮轴向力大约为1335 N。
3.2. 径向力
液下泵具有螺旋形压水室,在运作中会产生作用在叶轮上的径向力,使轴产生定向的挠度。螺旋形压水室设计工况下叶轮周围压水室中的速度和压力是均匀、轴对称的,理论上无径向力。但是当泵运行在非设计工况时叶轮表面压产生径向力 [6] [7]。径向力是垂直于旋转轴平面的合力。目前尚无精确的计算公式,只能按照下面的经验公式估算:
(2)
式中:Kr为实验系数,按照Steponoff公式计算。H为泵扬程(m)。D2为叶轮外径(m)。B2为叶轮包括盖板的出口宽度(m)。按上式计算得到叶轮径向力为1379 N。
3.3. 强度分析
1) 结构几何模型
通过Pro/e软件对泵轴进行三维建模,在WORKBENCH [8] [9] 软件中进行有限元网格划分、载荷约束施加,计算轴在工作状态下的结构应力及形变量,校核轴的强度是否满足要求。轴材料:2205双相不锈钢。为了计算的准确性,将泵轴按加工图三维建模,如图2所示,旋转轴为X轴。
2) 网格划分
网格划分与计算目标的匹配程度、网格的质量好坏,决定了后期有限元计算的质量,是有限元分析前处理的重中之重 [10]。液下泵轴结构相对简单,采用几何自动生成法,根据物理模型自动离散生成有限元网格模型见图3。对泵轴局部进行了加密处理,网格最大尺度控制在10 mm。
3) 载荷
泵运转过程中为:角速度:214.57 rad/s;轴向力:1335.6 N;径向力:1379.44 N,叶轮转子质量按最大重量加载,对轴承位置固定约束。
4) 结果分析
图4和图5为静应力强度与变形量。计算结果表明,泵运行过程中的最大静应力为28.5 Mpa,远小于材料的屈服极限。最大静应力位置在靠近调心滚子轴承的轴肩处,最大变形量为0.17255 × 10−3 mm,该处变形量可忽略不计。
4. 刚度分析
4.1. 挠度分析
泵轴除了校核强度以外要进行刚度计算 [11]。刚度计算包括轴的挠度(即轴运转中挠度小于叶轮与泵体之间的半径间隙)和轴的临界转速。泵运转过程中的挠度等于转子自重引起的静挠度,加上残余不平衡质量的离心惯性力引起的动挠度。
挠度是弯曲变形时横截面形心沿与轴线垂直方向的线位移,图6为液下泵轴挠度,可以看出Y方向的最大挠度值0.26 mm,位于轴头安装叶轮的位置,一般闭式叶轮的口环间隙为0.5 mm左右,液下轴长1340 mm时叶轮有可能与泵体发生刮擦,叶轮口环磨损加剧等现象,显然泵轴刚度稍微偏小,轴应加粗或在靠近叶轮位置布置支点;从图7可以看出,在轴悬臂位置的挠度值程线性变化,最小值位置调心滚子轴承的轴肩处。
Figure 7. Deflection Variation in Y Direction
图7. Y方向的挠度变化
4.2. 临界转速
当泵轴的转速增大到某一个转速时,轴运行变得不稳定,产生强烈的反复变形和振动。当转速继续增大时该现象反复出现,此现象为共振,产生共振的转速称为轴的临界转速,有第一、第二临界转速。传统的计算方法有:叠加法、能量法、逐次逼近法、图解法等 [12]。本文采用有限元法计算临界转速,ANSYS WORKBENCH模块充分考虑了系统质量、刚度、阻尼以及陀螺效应等因素。
4.2.1. 动力学方程
根据转子动力学和有限元理论,液下泵弹性结构系统的动力学通用方程为:
(3)
式中:M是结构整体的质量矩阵;C是结构整体的阻尼矩阵;K是结构整体的刚度矩阵;分别为结构整体的节点加速度矢量、速度矢量和位移矢量;为结构整体的外载荷矢量。对于结构系统无外部载荷作用,即计算泵轴的固有频率时,一般是无阻尼的自由振动,所以,系统初始状态开始的无阻尼自由振动,运动方程为:
(4)
液下泵轴的自由振动为简谐振动,及位移为正弦函数:式中:为模态固有频率。即上式为:
(5)
求解该该方程即可得到泵轴的各阶固有频率和振型。
4.2.2. 模态与振型
液下泵轴旋转频率为214.75 Hz,叶片数为6。由Modal计算可知轴的前6阶临界转速分别为200.91 Hz、202.2 Hz、597.5 Hz、598.96 Hz、910.08 Hz、1022.3 Hz,表1中为对应的主振型。可以看出轴的第一、第二临界转速均小于泵的转速2050 r/min。由于该型泵的转速范围960 r/min~2050 r/min,当泵运行在该范围的某一个转速时将会发生共振现象。泵的转速为第二临界转速的1.06倍,第三临界转速为泵的运行转速的2.69倍。该型液下泵运行的安全转速范围为:1.3倍第一临界转速与0.7倍第二临界转速之间,这样才能有效避免共振,不影响运转稳定。从轴临界转速计算可以看出,该型液下泵轴在设计过程中存在明显的缺陷,转速范围中存在某一共振点,因此在选型过程中应避开。
Table 1. Critical speed and vibration modes
表1. 临界转速与振型
5. 结论
本文采用ANSYS WORKBENCH有限元计算了液下泵轴的力学特性,详细分析了轴强度和刚度,给产品的选型和优化提供了相应的改型措施。研究表明,改型泵轴强度满足材料屈服应力,但叶轮和泵体之间的口环间隙有可能满足不了要求,第一、第二临界转速小于泵转速2050 r/min,存在共振区。