一种平行轴齿式压缩机临界负荷计算方法
【专利摘要】一种平行轴齿式压缩机临界负荷计算方法属于旋转机械转子动力学领域。具体步骤为：(1)获取计算平行轴齿式压缩机临界负荷需要的参数；(2)根据获得的参数计算不同功率下各个齿轮输出轴的动特性参数；(3)将动特性参数输入动力学方程，使用传递矩阵法计算不同功率下各个齿轮输出轴的1、2、3阶临界转速；(4)使用多项式拟合不同功率下各个齿轮输出轴的1、2、3阶临界转速，得到各阶临界转速的拟合函数；(5)根据各个齿轮输出轴的工作转速与临界转速拟合函数，判定各个齿轮输出轴是否在设计上存在临界转速并计算具体临界转速值。本发明用于平行轴齿式压缩机设计时各个齿轮输出轴是否存在临界负荷的验证，指导设计人员优化压缩机设计。
【专利说明】一种平行轴齿式压缩机临界负荷计算方法
【技术领域】
[0001] 本发明属于旋转机械转子动力学【技术领域】，涉及一种平行轴齿式压缩机临界负荷 计算方法。
【背景技术】
[0002] 平行轴齿式压缩机由于其运行效率高、结构紧凑等优点，在石油石化等流程工业 中压缩空气和二氧化碳方面得到广泛的应用。在现场实际运行过程中发现，平行轴齿式压 缩机在额定转速运行的状态下由于调整气动负荷产生振动过大的问题，这种问题引起的故 障就是临界负荷故障。我国设备故障诊断专家高金吉院士所著的机械故障诊治与自愈化一 书中针对临界负荷故障进行探讨与研究，并成功的解决了一台平行轴齿式压缩机的临界负 荷故障。
[0003] 临界负荷故障是平行轴齿式压缩机因运行条件改变而引起的一种振动过大的故 障，由于压缩机入口导叶变化而导致气动负荷改变，从而影响了压缩气体做功转子支撑轴 承的刚度变化，而轴承刚度改变会影响转子的各阶临界转速转速改变，当某一阶临界转速 在刚度变化的过程中与额定转速重合时，就会引发共振，导致机组振动异常，一旦振动过大 就会引发跳车，导致生产被迫停止。
[0004] 在平行轴齿式压缩机实际运行过程中发现，并不是所有输出轴都存在临界负荷故 障现象，通过动力学分析发现由于某些输出轴正常工作转速与各阶临界转速偏差较大，平 行轴齿式压缩机在负荷变化范围内各阶临界转速的变化并不会导致其与工作转速重合，因 此，这种平行轴齿式压缩机的输出轴并不存在发生临界负荷的可能。本发明提出的一种平 行轴齿式压缩机临界负荷计算方法，根据平行轴齿式压缩机的基本结构尺寸、个叶轮质心 位置、质量、直径转动惯量和极转动惯量，各个齿轮轴的正常运行工作转速，各齿轮的安装 角度、齿轮的安装位置、节圆直径、螺旋角及压力角以及额定运转功率计算输出轴是否存在 临界负荷，指导平行轴齿式压缩机设计人员在设计时，校验各个输出轴临界负荷问题，从而 在设计时规避各个输出轴发生临界负荷故障的可能。
【发明内容】
[0005] 本发明的目的是在于根据平行轴齿式压缩机的基本结构尺寸、个叶轮质心位置、 质量、直径转动惯量和极转动惯量，各个齿轮轴的正常运行工作转速，各齿轮的安装角度、 齿轮的安装位置、节圆直径、螺旋角及压力角以及额定运转功率等信息，对平行轴齿式压缩 机各个输出轴临界负荷进行计算，根据计算结果校验输出轴是否在运行过程中发生临界负 荷故障，从而在设计时规避临界负荷发生的可能，为平行轴齿式压缩机设计增加了设计依 据，保证了该类压缩机安全稳定的运行。
[0006] -种平行轴齿式压缩机临界负荷计算方法，其特征在于，包括：
[0007] 1)输入计算平行轴齿式压缩机临界负荷所需要的信息，包括：平行轴齿式压缩机 各个转轴的基本结构尺寸，各个叶轮质心位置、质量、直径转动惯量和极转动惯量，各个齿 轮轴的正常运行工作转速，各齿轮的安装角度、齿轮的安装位置信息、节圆直径、螺旋角及 压力角以及额定运转功率。
[0008] 2)计算平行轴齿式压缩机各个齿轮输出轴轴承在正常运行工作转速条件下不同 功率的动特性参数，包括刚度系数和阻尼系数。
[0009] 3)计算平行轴齿式压缩机各个齿轮输出轴在正常运行工作转速条件下不同功率 的1、2和3阶临界转速。
[0010] 4)使用二次多项式拟合算法拟合平行轴齿式压缩机各个输出轴正常运行工作转 速下不同功率的1、2和3阶临界转速，得到各个输出轴工作转速下不同功率1阶临界转速 拟合函数(P)、2阶临界转速拟合函数fe2 (p)、3阶临界转速拟合函数fe3 (p)，p的取值范 围为：0%彡P彡100%。
[0011] 5)计算平行轴齿式压缩机各个输出轴的1阶临界转速拟合函数fa (p)的变化 区间[fcl (0)，fcl (100) ]、2阶临界转速拟合函数fC2 (P)的变化区间[fC2 (0)，fC2 (100)]、 3阶临界转速拟合函数fC3(P)的变化区间[L(Ohfra(K)O)];若该输出轴的工作转速 Reve[fcl (0)，fcl (100)]或Reve[fC2 (0)，fC2 (100)]或Reve[fC3 (0)，fC3 (100)]，则该输 出轴存在临界负荷，临界负荷根据等式fra(P) =Rev计算得到，其中，i=I或2或3 ;若输 出轴的工作转速RevwLZn(O),./:.,(100)^^^4./:.2(0),./:.2(100)^^14./^(0^ 则该输出轴不存在临界负荷；
[0012] 计算各个齿轮轴轴承的动力特性参数，根据平行轴齿式压缩机各个输出轴额定运 转功率已，计算该功率下输出轴齿轮啮合力Fp并与输出轴的重力G进行矢量加和，得到当 前功率下各个输出轴齿轮哨合处的总受力=Fp+G;根据平行轴齿式压缩机各个输 出轴齿轮啮合处的总受力，调用轴承动特性计算程序，计算得到正常运行工作转速额定运 转功率^下的刚度系数与阻尼系数。根据上述计算方法计算正常运行工作转速下不同功 率的刚度系数与阻尼系数，功率取值范围Pe[0%，100% ]。
[0013] 计算平行轴齿式压缩机各个齿轮输出轴的临界转速，根据获得的输出轴的基本结 构尺寸建立转子模型以及各个叶轮质心位置、质量、直径转动惯量和极转动惯量，对各个转 子模型进行离散化处理，将连续的转子构造成一系列质量块组成的多自由度离散转子系 统，离散的各个子段的质量按照质量中心不变的原则分到两端节点上，转动惯量按照转动 惯量不变原则简化到两边端点上，简化后轴段的等效抗弯刚度与各个子段的抗弯刚度关 系，按照纯弯时两端截面相对转角不变求得。
[0014] 根据离散化的转子模型构建平行轴齿式压缩机各个输出轴转子系统的动力学方 程MV+Ci+心=F,式中，M为质量矩阵、C为阻尼矩阵、K为刚度矩阵，X为系统的广义坐 标矢量、F为作用在系统上的广义外力。将在正常运行工作转速条件下计算的不同功率的刚 度系数和阻尼系数带入各个输出轴转子系统的动力学方程的刚度矩阵和阻尼矩阵后，使用 传递矩阵法计算平行轴齿式压缩机各个输出轴在正常运行工作转速条件下不同功率的1、2 和3阶临界转速。
【专利附图】
【附图说明】
[0015]图1 一种平行轴齿式压缩机临界负荷计算方法流程图；
[0016] 图2平行轴齿式压缩机结构图；
[0017] 图3平行轴齿式压缩机齿轮输出轴1结构简图；
[0018] 图4平行轴齿式压缩机齿轮输出轴2结构简图；
[0019] 图5平行轴齿式压缩机齿轮输出轴3结构简图；
[0020] 图6齿轮输出轴1负荷变化与临界转速关系图；
[0021] 图7齿轮输出轴2负荷变化与临界转速关系图；
[0022] 图8齿轮输出轴3负荷变化与临界转速关系图；
【具体实施方式】
[0023] 下面结合附图和实例对本发明进行具体说明。
[0024] -种平行轴齿式压缩机临界负荷计算方法流程图如图1所示，主要包括：
[0025] (1)获取需要计算临界负荷的平行轴齿式压缩机的参数，包括：平行轴齿式压缩 机各个转轴的基本结构尺寸，各个叶轮质心位置、质量、直径转动惯量和极转动惯量，各个 齿轮轴的正常运行工作转速，各齿轮的安装角度、齿轮的安装位置信息、节圆直径、螺旋角 及压力角以及额定运转功率；
[0026] (2)根据齿轮轴齿轮安装位置、节圆直径、螺旋角及压力角以及额定运转功率，计 算各个齿轮输出轴齿轮啮合的总受力，根据输出轴的受力计算各个齿轮输出轴轴承的动特 性参数，包括刚度系数、阻尼系数；
[0027] (3)根据获得的输出轴的基本结构尺寸建立转子模型以及各个叶轮质心位置、质 量、直径转动惯量和极转动惯量，对各个转子模型进行离散化处理，建立各个齿轮输出轴的 动力学方程，使用传递矩阵法计算齿轮输出轴的1、2和3阶临界转速。
[0028] (4)计算齿轮输出轴运行功率在不同功率百分比p下的1、2、3阶临界转速，使 用二次多项式拟合算法拟合各个齿轮输出轴在不同负荷百分比下的1、2、3阶临界转速，得 到1阶临界转速的拟合函数fia(P)，2阶临界转速拟合函数fie2 (p)，3阶临界转速拟合函数 化3(口），1代表不同齿轮输出轴，1 = 1，2,3...，11;
[0029] (5)根据计算得到的临界转速拟合函数以及各个齿轮输出轴的工作转速，判别各 个齿轮输出轴是否存在临界负荷，如果存在则计算具体齿轮输出轴的临界负荷；
[0030] 本发明以某压缩机设计厂家设计的平行轴齿式压缩机为对象对该平行轴齿式压 缩机各个输出轴的临界负荷进行校验。该平行轴齿式压缩机共有5个转子，其中3个为输出 转子，如图2所示。根据平行轴齿式压缩机各个转轴的基本结构尺寸，各个叶轮质心位置、 质量、直径转动惯量和极转动惯量等信息分别建立齿轮输出轴的模型，图3为齿轮输出轴1 的结构简图、图4为齿轮输出轴2的结构简图、图5为齿轮输出轴3的结构简图。
[0031]各个齿轮的螺旋角：12度；压力角：20度。
[0032]
【权利要求】
1. 一种平行轴齿式压缩机临界负荷计算方法，其特征在于，包括： 1) 输入计算平行轴齿式压缩机临界负荷所需要的信息，包括：平行轴齿式压缩机各个 转轴的基本结构尺寸，各个叶轮质心位置、质量、直径转动惯量和极转动惯量，各个齿轮轴 的正常运行工作转速，各齿轮的安装角度、齿轮的安装位置信息、节圆直径、螺旋角及压力 角以及额定运转功率； 2) 计算平行轴齿式压缩机各个齿轮输出轴轴承在正常运行工作转速条件下不同功率 的动力特性参数，包括刚度系数和阻尼系数； 3) 计算平行轴齿式压缩机各个齿轮输出轴在正常运行工作转速条件下不同功率的1、 2和3阶临界转速； 4) 使用二次多项式拟合算法拟合平行轴齿式压缩机各个输出轴正常运行工作转速下 不同功率的1、2和3阶临界转速，得到各个输出轴工作转速下不同功率1阶临界转速拟合 函数(P)、2阶临界转速拟合函数fe2(p)、3阶临界转速拟合函数(p)，p的取值范围为： 0%^p^ 100% ； 5) 计算平行轴齿式压缩机各个输出轴的1阶临界转速拟合函数fa(p)的变化区 间[fcl (〇)，fcl (1〇〇) ]、2阶临界转速拟合函数fC2 (p)的变化区间[fC2 (0)，fC2 (100)]、 3阶临界转速拟合函数fra(p)的变化区间[&3(〇)，43(1〇〇)];若该输出轴的工作转速 RevG[fcl (〇)，fcl (1〇〇)]或RevG[fC2 (〇)，fC2 (1〇〇)]或RevG[fC3 (〇)，fC3 (1〇〇)]，则该输出 轴存在临界负荷，临界负荷根据等式fa(P) =Rev计算得到，其中，i= 1或2或3 ;若输出 轴的工作转速
则该输出轴不存在临界负荷。
2. 根据权利要求1所述的一种平行轴齿式压缩机临界负荷计算方法，其特征在于：所 述的计算各个齿轮轴轴承的动力特性参数，根据平行轴齿式压缩机各个输出轴额定运转功 率己，计算该功率下输出轴齿轮啮合力并与输出轴的重力石进行矢量加和，得到当前 功率下各个输出轴齿轮哨合处的总受力+G;根据平行轴齿式压缩机各个输出 轴齿轮啮合处的总受力，调用轴承动特性计算程序，计算得到正常运行工作转速额定运转 功率^下的刚度系数与阻尼系数；根据上述计算方法计算正常运行工作转速下不同功率的 刚度系数与阻尼系数，功率取值范围Pe[〇%，1〇〇% ]。
3. 根据权利要求1所述的一种平行轴齿式压缩机临界负荷计算方法，其特征在于：计 算平行轴齿式压缩机各个齿轮输出轴的临界转速，根据获得的输出轴的基本结构尺寸建立 转子模型以及各个叶轮质心位置、质量、直径转动惯量和极转动惯量，对各个转子模型进行 离散化处理，将连续的转子构造成一系列质量块组成的多自由度离散转子系统，离散的各 个子段的质量按照质量中心不变的原则分到两端节点上，转动惯量按照转动惯量不变原则 简化到两边端点上，简化后轴段的等效抗弯刚度与各个子段的抗弯刚度关系，按照纯弯时 两端截面相对转角不变求得； 根据离散化的转子模型构建平行轴齿式压缩机各个输出轴转子系统的动力学方程MV+Ci+iC^F，式中，M为质量矩阵、C为阻尼矩阵、K为刚度矩阵，x为系统的广义坐标 矢量、F为作用在系统上的广义外力；将在正常运行工作转速条件下计算的不同功率的刚 度系数和阻尼系数带入代入各个输出轴转子系统的动力学方程的刚度矩阵和阻尼矩阵后， 使用传递矩阵法计算平行轴齿式压缩机各个输出轴在正常运行工作转速条件下不同功率 的1、2和3阶临界转速。
【文档编号】F04C28/28GK104405642SQ201410536423
【公开日】2015年3月11日 申请日期:2014年10月12日 优先权日:2014年10月12日
【发明者】江志农, 张明, 冯坤, 马波, 高金吉 申请人:北京化工大学