压缩机复杂管路压力脉动及管道振动研究.docx

NuclearPowerEngineeringAug8月文章编号：0258-0926(2008)04-0079-05压缩机复杂管路压力脉动及管道振动研究徐斌，冯全科，余小玲（西安交通大学能源动力学院，西安，710049）摘要：围绕往复式压缩机管道系统的振动及往复式压缩机的管道压力脉动问题，依据平面波动理论，采用转移矩阵和刚度矩阵计算出复杂管路的气柱固有频率和压力脉动。借助于有限元方法的离散思想，建立了往复式压缩机管道振动及应力分析的数学模型，提出了恰当的边界条件，利用基于有限元的管道分析软件CAESARⅡ对模型进行求解，获得了管道系统的振动模态结果。对比试验结果与计算结果发现，利用一维平面波动方程可以比较准确地计算出往复式压缩机管路的气柱固有频率、压力脉动。关键词：往复式压缩机；复杂管路；压力脉动；管道振动中图分类号：TH457文献标识码：A阻尼系数与管道中气体流量成正比。根据气体连续方程、运动方程、波动方程可得到关于脉动压力pt的偏微分方程：1引言往复式压缩机广泛用于石油、化工、钢铁和冶金等行业。管道的强烈振动对压缩机性能及可靠性有很大的影响，压缩机管道的振动问题也越来越引起人们的重视。生产中遇到的压缩机装置振动绝大多数是气流脉动引起的[1]，而往复式压缩机所产生的气流脉动最为剧烈，对性能、工作的影响最为明显[2]。研究气流脉动的产生机理，建立合理的气流脉动数学模型，准确地进行气流脉动的预测，掌握控制气流脉动的有效方法，意义十分重大。近几年，国内外对气流脉动的研究工作主要集中于以下几个方面：完善转移矩阵，建立和完善管系单元的数学模型[3]；运用有限元的方法对气流脉动进行分析，并引入一些如CFD等新的方法[4～10]。但是这些计算气流脉动的方法都以整个管路为整体来研究，在对复杂管路进行计算时缺乏灵活性和通用性，不易于迅速解决工业实际中不同压缩机管道系统的问题。本文基于平面波动理论，考虑声学反射、流体粘性和摩擦情况，分段计算了复杂的往复式压缩机级间管路的气柱固有频率和气流脉动，并在往复式压缩机管路系统实验台进行了实验验证。4μ?3p2?2p?p?pt?t+Rt?ct=0(1)f?t?x2?x2?t3ρ?t0式中，p为管道中压力的脉动值，Pa；μ为介质t的粘性系数，(N·s)/m2；ρ0为密度的平均值，kg/m3；R′为摩擦力的脉动值，N；c为声速，m/s。f用分离变量的方法求解式(1)，得：??x??**pt=Aexp?jk?t?a??+?????x???B*exp?jk?t+??(2-1)a????1??x???**u=?Aexp?jk?t????tρ0SEi?a??????x???Bexpjkt+*????(2-2)?a??????R?ω2?1?if??ω???2k=(3)??4μω2??c?1+i2?3ρ0c??λQ=02基本理论考虑了气体内部粘性影响，并假设管壁摩擦R(4)fdS收稿日期：2007-07-26；修回日期：2007-10-2980核动力工程Vol.29.No.4.20081λ=(5)1.14+2lgd?2E=ck(6)ωS式中，ω为脉动的圆频率，Hz；Q0为管道中气体平均流量，m3/s；d为管道直径，m；S为管道的截面积，m2；Δ为管道的粗糙度。式(2)是经典平面波动方程的一般解。此二式中的复数常数A*、B*由管道的端点条件确定。根据声学原理可知，式(2-1)的第一项代表沿x轴正向传波的行波，第二项代表沿x轴负向传波的反射波[11]，管内压力场为叠加形成的驻波。所以，激发源处的压力由压缩机活塞运动引起的压力波动和反射回来的压力波叠加而成。这里，考虑极图1往复式压缩机管路系统Fig.1PipingSystemforReciprocatingCompressors表1往复式压缩机主要技术规格TechnicalSpecificationforReciprocatingCompressorTable1值情况，两列波的相位相同。令u*=2u，代入式1t(2-1)、式(2-2)，求得A*和B*。nf=mi(8)c603复杂管系的气流脉动计算用转移矩阵法计算气柱固有频率。依据平面波动原理，1点和4点的压力脉动p*和质量流量ξ*能够用转移矩阵表示式中，n为压缩机转速，r/min；m为简谐阶次，m=1，2，3，??，n。激发频率的计算结果见表2。当激发频率等于0.8～1.2倍气柱固有频率时产生气柱共振。表2管系的激发频率Table2ExcitationFrequencyofPipingSystem?p*??p*???=[M4?3][M3?2][M2?1]??41(7)ξξ??4*????1*??式中，[M]为各种元件的转移矩阵，详见文献[1]。4管道系统的振动建模及计算为了验证本文使用的计算方法和计算分析