ＤＨ型透平压缩机转子动力学性能计算方法技术

本发明专利技术涉及ＤＨ型透平压缩机转子动力学性能计算方法。具体计算步骤为：（１）选择分析对象和计算类型，输入轴承－转子系统的参数，完成转子系统的离散化等前处理操作，（２）调用轴承计算程序，计算各轴承的动、静特性参数，根据选定的分析对象和计算类型，组建考虑齿轮啮合的平行轴系的转子动力学方程，（３）调用求解函数计算动力学微分方程的特征值和特征向量，经有量纲化转换为动力学特性结果数据，完成系统稳定性、临界转速和强迫振动响应及振型的计算和分析。本发明专利技术在坚实的理论研究基础上开发了能用于工程实际的多平行轴系动力学性能计算分析手段，所计算的转子系统中可包含各类型支撑，较好的实现了轴承计算与转子系统动力学分析的结合。

【技术实现步骤摘要】

本专利技术涉及一种转子动力学性能计算方法，特别是一种DH型透平压缩机转子动力学性能计算方法。
技术介绍
齿轮耦合的轴承-转子系统中，多个轴系通过齿轮耦合联系在一起，系统具备由齿轮传动引起的新特性，横向振动通过齿轮传递后将引起转子产生扭转振动，即发生弯扭耦合振动。存在齿轮啮合的多平行轴系在风机、压缩机、增速器等机器中广泛存在，齿轮耦合的轴承-转子系统的典型应用之一就是双轴型透平压缩机，亦称DH型透平压缩机，其轴承-转子系统如图1所示。 齿轮的啮合作用使原本相互独立的多个轴承-转子系统联接在一起，从而使各轴承-转子系统的动力特性相互影响，整个系统的动力特性与单个子系统的动力特性大不一样。即使理想齿轮传动，齿坯刚度、轮齿刚度及其啮合刚度的时变性将使整个系统的动力特性产生变化，且齿轮啮合时存在的齿侧间隙、各种传动误差(包括齿轮各种制造误差、安装误差、热变形、磨损等)也将影响系统的动力学特性。关于齿轮啮合的研究，已经由最初的刚性啮合模型、弹性啮合模型发展到现在的世纪齿轮啮合刚度模型，数学模型与实际系统间的误差越来越小。具有齿轮耦合的轴承-转子系统最显著的特点就是弯扭耦合振动，弯扭耦合不仅影响系统的固有频率、不平衡响应，还影响系统的临界转速和稳定性。该类系统的设计过去一直采用转子单独设计的设计方法，忽略了各个转子之间的耦合关系，这种设计方法只能保证各个转子单独使用时，能够稳定运行，并不能保证耦合后还能正常运行。
技术实现思路
本专利技术的目的在于提供一种DH型透平压缩机转子动力学性能计算方法，有计算软件，计算快速，极大的方便使用者进行压缩机、风机等流体机械及增速器、减速器等具有齿轮传动的平行轴系的转子系统动力学分析，免除设计者需借助专门知识才能操作的要求，直接通过界面操作就可以进行转子动力学分析。 为了揭示一个系统的动力学行为，就必须建立一个符合实际系统的数学模型，建立多平行轴系统模型的关键是如何考虑齿轮啮合这一因素。首先建立轴承-转子系统的线性动力学模型，然后根据齿轮参数确定实际齿轮啮合刚度、啮合阻尼，并考虑齿面摩擦力、齿侧间隙、啮合脱齿及齿背啮合等情况及斜齿圆柱齿轮传动时轴向力对系统动力特性的影响，通过齿轮间的啮合力，将多个转子-轴承系统联系起来，最终建立齿轮耦合的轴承-转子系统的非线性动力学模型。针对系统动力学模型的求解，可完成系统动力稳定性、临界转速和动力响应等动力学特性分析，为系统动力学性能改进提供数值依据。 主要研究方法 1、转子系统的离散 一个实际转子是质量连续分布的弹性体，其上安装有叶轮、齿轮、联轴器、质量盘等零部件。可用集总参数法将转子离散成n个无质量的弹性轴段和n+1无弹性的集总质量，如图2所示。 2、齿轮耦合的轴承-转子系统的数学模型 (1)系统弯曲振动方程 取离散后第j轴段分析，图3为转子在xoz平面和在yoz平面内的力和力矩分布，取第j个截面在xoz平面和yoz平面中的状态变量为 T，由力、力矩及变形平衡可得第j轴段的运动方程 该单元上，外力矩、外力在x、y方向上的表达式为 式中kxx、kxy、kyx、kyy-轴承油膜刚度 dxx、dxy、dyx、dyy-轴承油膜阻尼 Mck、Nck-齿轮啮合力的轴向分量产生的弯矩 Pcx、Pcy-外加激励力、不平衡力、齿轮啮合力、重力或其它控制力 最终可得用位移、转角及其导数(x，y， ψ，......)表示的第j个单元力及力矩的平衡方程 考虑自由端部边界条件 经无量纲化，可得到转子-轴承系统的弯曲振动方程 式中，[M]、[C]、[K]、{F}分别为系统的总质量矩阵、阻尼矩阵、刚度矩阵和总激振力向量。位移向量{X}＝(X1，Φ1，…，Xj，Φj，…，Xn，Φn)T (2)扭转运动方程 图4为转子扭转振动的力学模型，其中Jzj为第j个轴段的集总极转动惯量，kθj为第j个轴段的扭转刚度，G为材料的剪切弹性模量，dj和lj分别为第j个轴段的直径和长度，θj为扭转角位移，Mp为齿轮啮合力产生的扭矩。 第j个单元的扭转振动方程为 经无量纲化扭转振动方程化为 其中 第一个和最后一个节点的扭转振动方程为 将转子各单元扭转振动方程集总，可得到转子扭转振动方程 (3)齿轮啮合力 图5为平行转子-轴承系统中，齿轮传动受力分析图，中间大齿轮为主动轮，驱动两侧平行配置的从动齿轮。在自由状态下，齿轮中心在o1、o2、o3，啮合线与齿轮中心连线间的夹角分别为α12、α13。当传递动力或受到外部作用力时，转子产生变形，齿轮中心移动到o′1、o′2、o′3，此时，啮合线方向发生变化，啮合线与齿轮中心连线的夹角分别为α′12、α′13，并且，齿轮1、2中心、齿轮1、3中心连线与x轴方向的夹角分别为β12、β13。 当不考虑齿面间的摩擦力时，可求得齿轮间的啮合力如下 式中， 为齿轮所在节点的位移向量，上标1、2、3表示转子类别，下标g、q、p表示齿轮在该转子上所处的节点位置，r1、r2、r3为齿轮的基圆半径， 分别为齿轮1、2及齿轮1、3之间的平均啮合刚度，c12、c13为齿轮的啮合阻尼，且 式中m1，m2，m3-分别为三个齿轮的等效质量 ξ-齿轮啮合阻尼比，可取0.05。 在齿轮实际啮合过程中，啮合点沿啮合线作周期性变化，所以，啮合力轴向分量的位置也在不断变化，但变化距离相对齿轮尺寸来说很小，故可近似认为啮合力始终作用在节点处。齿轮间的啮合力对相互啮合的一对齿轮而言，实际上是一对作用力和反作用力，将该作用力分别沿坐标轴方向投影可得 式中Fxi、Fyi、Fzi、Mxi、Myi、Ti为由于齿轮相互啮合而分别作用在齿轮1、齿轮2和齿轮3上的作用力和力矩，β为齿轮螺旋角，ri为齿轮的基圆半径(i＝1，2，3)。 忽略啮合线方向变化对系统的影响，并进行无量纲化处理可得 可得齿轮节点处齿轮啮合所产生的作用力和作用力矩 而在非齿轮节点处，Pc＝0，Mc＝0，Mp＝0。将各个转子的弯曲振动方程和扭转运动方程进行组装，然后将三个转子的运动方程再组装在一起，将上式代入后，可得到系统总的运动方程式 式中M-系统总质量矩阵；C-系统总阻尼矩阵；K-系统总刚度矩阵；F-系统总激振力列矢量。 3、系统的数学分析方法 (1)固有振动分析 当系统没有承受任何外界激励力时，系统固有振动方程为 设解具有一般形式X＝X0eλT，λ＝-U+iV，代入上式后，可得(Mλ2+Cλ+K)X0＝0 式中λ是特征值，其实部反映该阶模态对应的模态阻尼的大小，虚部为该阶模态下的阻尼固有频率；X0是特征向量，反映该阶模态下的振型，轴承-转子系统固有振动分析问题就转换为式(2-18)的二次特征值和特征向量问题 (2)系统强迫振动分析 在齿轮耦合的转子-轴承系统中，由于本文档来自技高网...

【技术保护点】
一种ＤＨ型透平压缩机转子动力学性能计算方法，其特征在于具体计算步骤如下：（１）多平行轴系统动力学分析前处理：①在压缩机、减速器、增速器等中选择计算对象，选定特征值及强迫振动响应计算或临界转速及失稳转速计算，同时可选择考虑齿轮耦合的多平行轴系统分析或其中某一根转子分析。②输入系统各转子及轴承的参数：包括轴段参数、附加质量参数、不平衡质量参数、齿轮参数及轴承类型、位置和动特性系数等。若进行某一转速下的特征值及振型分析和不平衡响应分析，输入主动轴的转速（耦合系统分析）或某一转子的转速（单根转子分析）。当进行临界转速及失稳转速计算时，输入主动轴的最低转速和最高转速（耦合系统分析）或某一转子（单根转子分析）的最低转速和最高转速，并选择转速步长。③按集总参数法，对各个转子进行离散化处理，将连续变截面的转子处理成由一系列无质量、刚度不变的轴段联结一系列质量块组成的多自由度离散系统，离散中各子段的质量按照质心不变的原则分到两端节点上，转动惯量按转动惯量不变原则简化到两边端点上，简化后轴段的等效抗弯刚度与各子段的抗弯刚度关系，按纯弯时两端截面相对转角不变求得。（２）系统动力学方程的建立：①若轴系支撑为滑动轴承，计算各轴承的载荷，调用轴承性能计算程序计算出滑动轴承的油膜刚度、阻尼系数表，每个轴承的实际工作点可根据载荷用插值方法求得，不同转速下的动特性参数包括４个刚度系数和４个阻尼系数也用插值方法求得，较好的实现了轴承分析与转子系统动力学分析的结合。②在传递矩阵法的基础上，形成转子系统的弯曲振动方程和扭转振动方程，并计入齿轮啮合力得到齿轮啮合的刚度矩阵和阻尼矩阵，最终得到系统总质量矩阵、阻尼矩阵、刚度矩阵及激振力矩阵：Ｍ＊＋Ｃ＊＋ＫＸ＝Ｆ，其中｛Ｘ｝＝（Ｘ↓［１］，Φ↓［ｌ］，…，Ｘ↓［ｊ］，Φ↓［ｊ］，…，Ｘ↓［ｎ］，Φ↓［ｎ］）↑［Ｔ］。（３）动力学方程求解及结果后处理输出：①特征值和强迫振动响应计算结果输出：可绘制给定转速下的８阶或１６阶特征值的实部（对数衰减率）、虚部（固有涡动频率）的总图和各阶振型图；可输出强迫振动响应沿轴向的分布图。当系统包括２～３根转子时，由左至右依次分段表示第１、第２、第３根转子的振型或响应，包括Ｘ向、Ｙ向弯曲模态及扭转模态。②临界转速及失稳转速计算结果输出：包括在给定转速范围内各阶特征值虚部（固有涡动频率）及实部（对数衰减率）与主动轴（耦合系统分析时）或某一轴（单根转子分析时）转速的关系曲线；以红色...

【技术特征摘要】
1.一种DH型透平压缩机转子动力学性能计算方法，其特征在于具体计算步骤如下(1)多平行轴系统动力学分析前处理①在压缩机、减速器、增速器等中选择计算对象，选定特征值及强迫振动响应计算或临界转速及失稳转速计算，同时可选择考虑齿轮耦合的多平行轴系统分析或其中某一根转子分析。②输入系统各转子及轴承的参数包括轴段参数、附加质量参数、不平衡质量参数、齿轮参数及轴承类型、位置和动特性系数等。若进行某一转速下的特征值及振型分析和不平衡响应分析，输入主动轴的转速(耦合系统分析)或某一转子的转速(单根转子分析)。当进行临界转速及失稳转速计算时，输入主动轴的最低转速和最高转速(耦合系统分析)或某一转子(单根转子分析)的最低转速和最高转速，并选择转速步长。③按集总参数法，对各个转子进行离散化处理，将连续变截面的转子处理成由一系列无质量、刚度不变的轴段联结一系列质量块组成的多自由度离散系统，离散中各子段的质量按照质心不变的原则分到两端节点上，转动惯量按转动惯量不变原则简化到两边端点上，简化后轴段的等效抗弯刚度与各子段的抗弯刚度关系，按纯弯时两端截面相对转角不变求得。(2)系统动力学方程的建立①若轴系支撑为滑动轴承，计算各轴承的载荷，调用轴承性能计算程序计算出滑动轴承的油膜刚度、阻尼系数表，每个轴承的实际工作点可根据载荷用插值方法求得，不同转速下的动特性参数包括4个刚度系数和4个阻尼系数也用插值方法求得，较好的实现了轴承分析与转子系统动力学分析的结合。②在传递矩阵法的基础上，形成转子系统的弯曲振动方程和扭转振动方程，并计入齿轮啮合力得到齿轮啮合的刚度矩阵和阻尼矩阵，最终得到系统总质量矩阵、阻尼矩阵、刚度矩阵及激振力矩阵其中{X}＝(X1，Φl，…，Xj，Φj，...

【专利技术属性】
技术研发人员：朱爱斌，丘大谋，张莉娟，张锁怀，
申请(专利权)人：朱爱斌，
类型：发明
国别省市：87[中国|西安]