[发明专利]一种基于最小二乘法的密封动力特性系数获取方法

说明书

本发明提供一种基于最小二乘法的密封动力特性系数获取方法，包括如下步骤：设定转子的力和位移为线性关系，忽略涡动高阶导数的影响，利用转子密封动力特性系数的刚度和阻尼矩阵得到力的位移模型；将转子密封中心O作为以x为横坐标，y为纵坐标的坐标系的原点来建立坐标系；设定转子在t时刻的偏心位置上以椭圆轨迹涡动，将偏心位置作为涡动中心，得到椭圆涡动轨迹的参数方程；根据参数方程、位移模型和最小二乘法得到x方向和y方向的残差值；根据两个方向的残差值来确定总残差平方和S的函数表达式；利用S的函数表达式对八个密封动力特性系数求偏导数，得到对八个密封动力特性系数求偏导数的方程组，并对其进行求解得到八个密封动力特性系数。

技术领域

本发明涉及一种透平机械密封动力特性系数数值获取方法，具体涉及一种基于最小二乘法的密封动力特性系数获取方法。

背景技术

由于旋转透平机械逐步向超大型、高压比、高转速的发向发展，转子的工作条件越来越恶劣，由密封引起的气流激振日益突出，对旋转透平机械的稳定运行及其产生的经济效益影响较大。研究表明，密封内流体流动是气流激振力的主要来源之一，人们通常采用八个刚度与阻尼系数来表示密封动力特性。

目前，对密封动力特性系数的识别主要分为实验法、理论法和计算流体力学法(CFD)三种方法。与计算流体力学方法相比，采用实验法在做实验过程中一方面往往会受到实验条件的限制和影响以至于难以模拟复杂工况，现如今旋转机械正往高参数条件下发展，无疑在做实验过程中很大程度上增加了难度，不利于实验的完成，且实验在流场细节等方面的测试较为困难。目前理论法中使用最为广泛的是Bulk Flow模型方法，该方法假设条件太多，引起假设与实际密封腔内气流的流动规律差别较大，从而导致求解精度较低。然而，计算流体力学法就可以很好的解决以上两种方法中所提出的问题。

随着计算机硬件的高速发展及其计算能力的不断提高，计算流体力学方法正在受到越来越多的应用。CFD方法对任意的密封几何和操作工况具有普适性，并且对密封流体流动特征提供了可视化的分析技术。目前运用数值模拟对密封动力特性识别的方法分为稳态法和瞬态法两类。稳态法中，假设转子做理想状态下的同轴旋转运动，在转子中心建立旋转坐标系，可将非稳态的运动转化为稳态。但稳态法在实际应用过程中局限性较多，如：(1)转子实际涡动轨迹较为复杂，如存在偏心、椭圆涡动等；(2)密封动力特性系数与转速、涡动频率、偏心率等都存在依赖关系；(3)对于像蜂窝式密封、孔型密封、袋式密封等较为复杂的密封壁面的模拟难以实现。瞬态法中可直接设置转子涡动运动轨迹来模拟实际流场，通过在不同方向施加激振计算密封气流力及密封动力特性系数，避免了运用稳态法时遇到的问题。但是一些学者为了简化计算采用了不合理的假设，如：(1)忽略涡动轨迹对动力特性系数的影响；(2)直接假设涡动轨迹为圆形等。以上两个假设导致从源头上不能同时识别密封的八个动力特性系数且识别精度较低，甚至导致识别结果错误。最小二乘法是一种被广泛应用的数学优化技术，并广泛运用于曲线拟合、参数求解与优化设计等领域，它通过最小化残差的平方和寻找数据的最佳函数匹配。利用最小二乘法可以简便的求得未知的数据，并使得这些求得的数据与实际数据之间误差的平方和为最小，求得拟合函数中的未知系数。

因此，需要设计一种能够解决上述问题的方法。

发明内容

本发明是为了解决上述问题而进行的，目的在于提供一种基于最小二乘法的密封动力特性系数获取方法。

本发明提供了一种基于最小二乘法的密封动力特性系数获取方法，用于对转子的密封动力特性系数进行获取，具有这样的特征，包括如下步骤：

步骤1，设定转子的力和位移为线性关系，并忽略涡动高阶导数带来的影响，利用转子的密封动力特性系数的刚度以及阻尼矩阵表示转子的力的位移模型：

步骤2，将转子的密封中心O作为以x为横坐标，y为纵坐标的坐标系的原点来建立(x，y)坐标系；