CN116090135B - 一种阻尼器-转子系统响应分析方法 - Google Patents
一种阻尼器-转子系统响应分析方法 Download PDFInfo
- Publication number
- CN116090135B CN116090135B CN202310218507.3A CN202310218507A CN116090135B CN 116090135 B CN116090135 B CN 116090135B CN 202310218507 A CN202310218507 A CN 202310218507A CN 116090135 B CN116090135 B CN 116090135B
- Authority
- CN
- China
- Prior art keywords
- damper
- oil film
- rotor system
- force
- rotor
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
- 230000004044 response Effects 0.000 title claims abstract description 36
- 238000004458 analytical method Methods 0.000 title claims abstract description 18
- 239000003921 oil Substances 0.000 claims description 122
- 238000000034 method Methods 0.000 claims description 32
- 238000006073 displacement reaction Methods 0.000 claims description 14
- 238000004422 calculation algorithm Methods 0.000 claims description 7
- 238000013016 damping Methods 0.000 claims description 6
- 239000010729 system oil Substances 0.000 claims description 5
- 239000011159 matrix material Substances 0.000 claims description 4
- 238000012546 transfer Methods 0.000 claims description 3
- 238000013461 design Methods 0.000 description 6
- 238000001125 extrusion Methods 0.000 description 4
- 230000008569 process Effects 0.000 description 4
- 238000012360 testing method Methods 0.000 description 4
- 238000004364 calculation method Methods 0.000 description 3
- 239000010687 lubricating oil Substances 0.000 description 3
- 230000005540 biological transmission Effects 0.000 description 2
- 230000000694 effects Effects 0.000 description 2
- 230000009467 reduction Effects 0.000 description 2
- 238000011160 research Methods 0.000 description 2
- 238000004088 simulation Methods 0.000 description 2
- 238000006467 substitution reaction Methods 0.000 description 2
- 230000009471 action Effects 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 239000012141 concentrate Substances 0.000 description 1
- 238000007796 conventional method Methods 0.000 description 1
- 238000007405 data analysis Methods 0.000 description 1
- 238000012938 design process Methods 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/10—Geometric CAD
- G06F30/17—Mechanical parametric or variational design
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2119/00—Details relating to the type or aim of the analysis or the optimisation
- G06F2119/14—Force analysis or force optimisation, e.g. static or dynamic forces
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02T—CLIMATE CHANGE MITIGATION TECHNOLOGIES RELATED TO TRANSPORTATION
- Y02T90/00—Enabling technologies or technologies with a potential or indirect contribution to GHG emissions mitigation
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Geometry (AREA)
- Theoretical Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Evolutionary Computation (AREA)
- Computer Hardware Design (AREA)
- General Engineering & Computer Science (AREA)
- Pure & Applied Mathematics (AREA)
- Mathematical Optimization (AREA)
- Mathematical Analysis (AREA)
- Computational Mathematics (AREA)
- Lubricants (AREA)
- Vibration Prevention Devices (AREA)
Abstract
本发明属于航空发动机及燃气轮机技术领域,具体涉及一种阻尼器‑转子系统响应分析方法。本发明将考虑油膜惯性力的油膜力引入阻尼器‑转子系统,建立带阻尼器的不平衡响应分析方法,分析精度更高。
Description
技术领域
本发明属于航空发动机及燃气轮机数据分析技术领域,具体涉及一种阻尼器-转子系统响应分析方法。
背景技术
支点阻尼器(见图1)是航空发动机等高性能旋转机械常用的减振结构,通过改变支点处刚度、阻尼参数,改善转子系统动力学特性,减振及降低支承外传力,从而提高轴承寿命、耐久性。
阻尼器设计应包含两个方面(见图2),阻尼器刚度、阻尼综合设计、阻尼器—转子系统动力学设计。这就要求支承阻尼器设计过程中,首先要通过雷诺方程进行油膜力仿真分析,其次建立阻尼器-转子系统动力学响应分析方法,对阻尼器—转子系统动力学特性进行分析,最后通过动力学试验验证阻尼器在转子系统的减振效果。
工程上,一般采用阻尼器短轴承半油膜经典理论对阻尼器油膜力进行分析,这样往往会对阻尼器雷诺方程进行简化,其中就包括油膜惯性。对于高频工作的挤压油膜阻尼器雷诺数较大,此时就必须考虑油膜惯性力。随着发动机研制向高压、高转速的方向发展,发现依据经典雷诺简化方程进行油膜力已无法满足设计需要,设计中发现随着转速增加,油膜惯性力对阻尼器-转子系统动力学特性影响很大。目前,相关研究集中阻尼器阻尼特性引入油膜惯性力,将考虑油膜惯性力的油膜力引入阻尼器-转子系统的研究较少,极大的影响了发动机支点阻尼器设计。
发明内容
有鉴于此,本发明提出一种阻尼器-转子系统响应分析方法,将考虑油膜惯性力的油膜力引入阻尼器-转子系统,建立带阻尼器的不平衡响应分析方法,分析精度更高。
为了实现上述技术目的,本发明所采用的具体技术方案为:
一种阻尼器-转子系统响应分析方法,包括以下步骤:
S1、基于能量近似法求解所述阻尼器-转子系统的考虑了油膜惯性力的雷诺方程,得到阻尼器-转子系统的考虑了油膜惯性力的油膜力;
S2、基于所述油膜力建立阻尼器-转子系统的动力学模型及方程;
S3、建立所述油膜力与所述阻尼器-转子系统的油膜轴颈位移关系,基于弦割法的数值迭代方法,求解阻尼器-转子系统不平衡响应。
进一步的,所述S1具体包括以下步骤:
S101、建立坐标系下的所述雷诺方程:
X方向指向所述轴颈表面的圆周方向;
Y方向指向所述轴颈表面圆周的法线方向;
Z方向为沿轴颈的轴向方向;
R为阻尼器-转子系统的挤压油膜阻尼器半径;
h为阻尼器-转子系统的油膜厚度;
p为阻尼器-转子系统的油膜压力;
S102、基于上述雷诺方程,求解出考虑油膜惯性力的油膜力:
式中:F SFD—— 油膜力;
R———挤压油膜阻尼器半径
Re———雷诺数
L———阻尼器-转子系统的阻尼器的油膜长度
C———阻尼器-转子系统的油膜阻尼器油膜间隙
进一步的,所述S2具体包括以下步骤:
S201、建立所述阻尼器-转子系统的转子-支承系统的动力学方程公式:
S202、对所述动力学方程公式求解,获得转子-支承系统的支承振动幅值的公式:
进一步的,所述S3具体包括以下步骤:
S301、通过传递矩阵法的数学分析推导出所述油膜轴颈位移关系:
其中:S SFD—— 挤压油膜阻尼器油膜轴颈位移;
F SFD—— 油膜力;
C———挤压油膜阻尼器油膜间隙
A 1———系数
A 2———系数;
定义F SFD是S SFD的函数
定义
定义迭代算法为:
进一步的,所述阻尼器-转子系统中的阻尼器为弹性环式挤压油膜阻尼器。
附图说明
为了更清楚地说明本公开实施例的技术方案,下面将对实施例中所需要使用的附图作简单地介绍,显而易见地,下面描述中的附图仅仅是本公开的一些实施例,对于本领域普通技术人员来讲,在不付出创造性劳动的前提下,还可以根据这些附图获得其它的附图。
图1为本发明背景技术中支点阻尼器的结构示意图;
图2为本发明背景技术中支点阻尼器—转子系统的结构示意图;
图3为本发明具体实施方式中一种阻尼器-转子系统响应分析方法的流程示意图;
图4为本发明具体实施方式中阻尼器坐标及有关尺寸Z向示意图;
图5为本发明具体实施方式中阻尼器坐标及有关尺寸X向示意图;
其中:1、油膜环;2、油膜轴颈。
具体实施方式
下面结合附图对本公开实施例进行详细描述。
以下通过特定的具体实例说明本公开的实施方式,本领域技术人员可由本说明书所揭露的内容轻易地了解本公开的其他优点与功效。显然,所描述的实施例仅仅是本公开一部分实施例,而不是全部的实施例。本公开还可以通过另外不同的具体实施方式加以实施或应用,本说明书中的各项细节也可以基于不同观点与应用,在没有背离本公开的精神下进行各种修饰或改变。需说明的是,在不冲突的情况下,以下实施例及实施例中的特征可以相互组合。基于本公开中的实施例,本领域普通技术人员在没有作出创造性劳动前提下所获得的所有其他实施例,都属于本公开保护的范围。
要说明的是,下文描述在所附权利要求书的范围内的实施例的各种方面。应显而易见,本文中所描述的方面可体现于广泛多种形式中,且本文中所描述的任何特定结构及/或功能仅为说明性的。基于本公开,所属领域的技术人员应了解,本文中所描述的一个方面可与任何其它方面独立地实施,且可以各种方式组合这些方面中的两者或两者以上。举例来说,可使用本文中所阐述的任何数目个方面来实施设备及/或实践方法。另外,可使用除了本文中所阐述的方面中的一或多者之外的其它结构及/或功能性实施此设备及/或实践此方法。
还需要说明的是,以下实施例中所提供的图示仅以示意方式说明本公开的基本构想,图式中仅显示与本公开中有关的组件而非按照实际实施时的组件数目、形状及尺寸绘制,其实际实施时各组件的型态、数量及比例可为一种随意的改变,且其组件布局型态也可能更为复杂。
另外,在以下描述中,提供具体细节是为了便于透彻理解实例。然而,所属领域的技术人员将理解,可在没有这些特定细节的情况下实践所述方面。
在本发明的一个实施例中,提出一种阻尼器-转子系统响应分析方法,如图3所示,包括以下步骤:
S1、基于能量近似法求解阻尼器-转子系统的考虑了油膜惯性力的雷诺方程,得到阻尼器-转子系统的考虑了油膜惯性力的油膜力;
S2、基于油膜力建立阻尼器-转子系统的动力学模型及方程;
S3、建立油膜力与阻尼器-转子系统的油膜轴颈2位移关系,基于弦割法的数值迭代方法,求解阻尼器-转子系统不平衡响应。
在本实施例中,S1具体包括以下步骤:
S101、建立坐标系下的雷诺方程:
X方向指向所述轴颈表面的圆周方向;
Y方向指向所述轴颈表面圆周的法线方向;
Z方向为沿轴颈的轴向方向;
R为阻尼器-转子系统的挤压油膜阻尼器半径;
h为阻尼器-转子系统的油膜厚度;
p为阻尼器-转子系统的油膜压力;
S102、基于上述雷诺方程,求解出考虑油膜惯性力的油膜力:
式中:F SFD—— 油膜力;
R———挤压油膜阻尼器半径
Re———雷诺数
L———阻尼器-转子系统的阻尼器的油膜长度
C———阻尼器-转子系统的油膜阻尼器油膜间隙
在本实施例中,S2具体包括以下步骤:
S201、建立阻尼器-转子系统的转子-支承系统的动力学方程公式:
S202、对动力学方程公式求解,获得转子-支承系统的支承振动幅值的公式:
在本实施例中,S3具体包括以下步骤:
S301、通过传递矩阵法的数学分析推导出油膜轴颈2位移关系:
其中:S SFD—— 挤压油膜阻尼器油膜轴颈2位移;
F SFD—— 油膜力;
C———挤压油膜阻尼器油膜间隙
A 1———系数
A 2———系数;
定义F SFD是S SFD的函数
定义
定义迭代算法为:
以下基于采用弹性环式挤压油膜阻尼器的阻尼器-转子系统进行进一步说明,阻尼器的油膜环1-油膜轴颈2及其坐标有关尺寸如图4及图5所示。
首先基于能量近似法求解考虑油膜惯性力的雷诺方程,可得到考虑油膜惯性力的油膜力。
经典雷诺方程为:
X方向指向轴颈表面的圆周方向;
Y方向指向轴颈表面圆周的法线方向;
Z方向为沿轴颈的轴向方向;
考虑油膜惯性力求解出的油膜力为:
式中:F SFD—— 油膜力;
R———挤压油膜阻尼器半径
Re———雷诺数
L———阻尼器-转子系统的阻尼器的油膜长度
C———阻尼器-转子系统的油膜阻尼器油膜间隙
其次,建立某型发动机阻尼器—转子系统动力学模型及方程。
转子-支承系统的动力学方程公式:
对上述动力学方程的第四行求解,获得支承振动幅值的公式为:
最后,通过建立油膜力与油膜轴颈2位移关系,基于弦割法的数值迭代方法,求解转子系统不平衡响应。
通过传递矩阵法的数学分析可以推导出公式(3):
其中:S SFD—— 挤压油膜阻尼器油膜轴颈2位移;
F SFD—— 油膜力;
C———挤压油膜阻尼器油膜间隙
A 1———系数
A 2———系数;
除了S SFD和F SFD,其他都为已知参数。很明显,F SFD是S SFD的函数,如下式(5)
将式(5)代入式(3)
即
S SFD是方程(7)的复数根,基于此通过弦割法求解方程(6)的复数根。
定义
定义迭代算法如下:
迭代算法是为了解方程(7)。
首先假定挤压油膜阻尼器油膜轴颈2位移的初始系数x 1和x 2,通过代入其他参数计算函数g(x)、油膜力、不平衡力和系数等。最后计算出x 3(通常是x n+1)。通过迭代算法检验结果是否收敛,如果在明确的误差范围内,这个计算过程被认为收敛,x n+1即为挤压油膜阻尼器油膜轴颈2位移。通过x n+1可求出两盘的位移。该过程中,随着转速增加,都将重复该过程。如果不满足差值,迭代作用的结果并没有超过极限误差,用n+1替代n,算法将继续直到满足误差或者说迭代作用的极限结果将超出公差值。事实上,计算结果收敛满足每一个转速,所以结果是准确的,与传统方法的结果对比如表1所示。
表 1理论、试验结果对比
由上表可知,与现有技术相比,本实施例的有益效果是:本实施例设计的一种考虑支点阻尼器油膜惯性力的转子系统响应分析方法,与发动机动力学响应进行对比,试验结果表明,考虑油膜惯性力的动力学响应相较于不考虑油膜惯性力的动力学响应更接近于试验结果,仿真精度提升30%。
以上所述,仅为本公开的具体实施方式,但本公开的保护范围并不局限于此,任何熟悉本技术领域的技术人员在本公开揭露的技术范围内,可轻易想到的变化或替换,都应涵盖在本公开的保护范围之内。因此,本公开的保护范围应以权利要求的保护范围为准。
Claims (2)
1.一种阻尼器-转子系统响应分析方法,其特征在于,包括以下步骤:
S1、基于能量近似法求解所述阻尼器-转子系统的考虑了油膜惯性力的雷诺方程,得到阻尼器-转子系统的考虑了油膜惯性力的油膜力;
S2、基于所述油膜力建立阻尼器-转子系统的动力学模型及方程;
S3、建立所述油膜力与所述阻尼器-转子系统的油膜轴颈位移关系,基于弦割法的数值迭代方法,求解阻尼器-转子系统不平衡响应;
其中:所述S1具体包括以下步骤:
S101、建立坐标系下的所述雷诺方程:
X方向指向所述轴颈表面的圆周方向;
Y方向指向所述轴颈表面圆周的法线方向;
Z方向为沿轴颈的轴向方向;
R为阻尼器-转子系统的挤压油膜阻尼器半径;
h为阻尼器-转子系统的油膜厚度;
p为阻尼器-转子系统的油膜压力;
S102、基于上述雷诺方程,求解出考虑油膜惯性力的油膜力:
式中:F SFD—— 油膜力;
Re———雷诺数
L———阻尼器-转子系统的阻尼器的油膜长度
C———阻尼器-转子系统的油膜阻尼器油膜间隙
所述S2具体包括以下步骤:
S201、建立所述阻尼器-转子系统的转子-支承系统的动力学方程公式:
S202、对所述动力学方程公式求解,获得转子-支承系统的支承振动幅值的公式:
所述S3具体包括以下步骤:
S301、通过传递矩阵法的数学分析推导出所述油膜轴颈位移关系:
其中:S SFD—— 挤压油膜阻尼器油膜轴颈位移;
F SFD—— 油膜力;
A 1———系数
A 2———系数;
定义F SFD是S SFD的函数
定义
定义迭代算法为:
2.根据权利要求1所述的阻尼器-转子系统响应分析方法,其特征在于,所述阻尼器-转子系统中的阻尼器为弹性环式挤压油膜阻尼器。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202310218507.3A CN116090135B (zh) | 2023-03-09 | 2023-03-09 | 一种阻尼器-转子系统响应分析方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202310218507.3A CN116090135B (zh) | 2023-03-09 | 2023-03-09 | 一种阻尼器-转子系统响应分析方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN116090135A CN116090135A (zh) | 2023-05-09 |
CN116090135B true CN116090135B (zh) | 2023-06-30 |
Family
ID=86212192
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202310218507.3A Active CN116090135B (zh) | 2023-03-09 | 2023-03-09 | 一种阻尼器-转子系统响应分析方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN116090135B (zh) |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
KR20140114173A (ko) * | 2013-03-18 | 2014-09-26 | 한국항공우주산업 주식회사 | 비행시험 데이터를 이용한 hqs 모델링과 튜닝을 위한 장치 및 그 제어방법 |
CN210565782U (zh) * | 2019-10-09 | 2020-05-19 | 康磊 | 一种挤压油膜阻尼轴承 |
CN114091314A (zh) * | 2022-01-24 | 2022-02-25 | 天津仁爱学院 | 一种基于磁流变阻尼器的转子系统模型的振动预测方法 |
CN115688387A (zh) * | 2022-10-13 | 2023-02-03 | 中国航发四川燃气涡轮研究院 | 一种基于轴承寿命的支点阻尼器油膜间隙选取方法 |
-
2023
- 2023-03-09 CN CN202310218507.3A patent/CN116090135B/zh active Active
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
KR20140114173A (ko) * | 2013-03-18 | 2014-09-26 | 한국항공우주산업 주식회사 | 비행시험 데이터를 이용한 hqs 모델링과 튜닝을 위한 장치 및 그 제어방법 |
CN210565782U (zh) * | 2019-10-09 | 2020-05-19 | 康磊 | 一种挤压油膜阻尼轴承 |
CN114091314A (zh) * | 2022-01-24 | 2022-02-25 | 天津仁爱学院 | 一种基于磁流变阻尼器的转子系统模型的振动预测方法 |
CN115688387A (zh) * | 2022-10-13 | 2023-02-03 | 中国航发四川燃气涡轮研究院 | 一种基于轴承寿命的支点阻尼器油膜间隙选取方法 |
Non-Patent Citations (5)
Title |
---|
Dynamical analysis of an elastic ring squeeze film damper-rotor system;Zhifei Han等;Mechanism and Machine Theory;第131卷;第406-419页 * |
基于有限元素法的挤压油膜动力特性研究;邓吟;林力;史亚杰;洪杰;;北京航空航天大学学报(第05期);第531-535页 * |
对挤压油膜阻尼器轴承和旋转机械转子―挤压油膜阻尼器轴承系统动力特性研究的回顾与展望;夏南,孟光;机械强度(第02期);第216-224,231页 * |
带挤压油膜阻尼器的转子系统动力学建模与非线性研究;刘一鸣;中国优秀硕士学位论文全文数据库(1);第C031-208页 * |
非线性挤压油膜阻尼器转子系统的非协调运动分析;陈照波,焦映厚,夏松波,黄文虎;中国机械工程(第09期);第774-776页 * |
Also Published As
Publication number | Publication date |
---|---|
CN116090135A (zh) | 2023-05-09 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Xu et al. | Vibration characteristics of bearing-rotor systems with inner ring dynamic misalignment | |
CN109829262B (zh) | 一种转子-轴承系统非线性动力学分析方法 | |
Ma et al. | Effects of eccentric phase difference between two discs on oil-film instability in a rotor–bearing system | |
Lu et al. | Nonlinear response analysis for an aero engine dual-rotor system coupled by the inter-shaft bearing | |
Li et al. | A general dynamic model coupled with EFEM and DBM of rolling bearing-rotor system | |
Rashidi et al. | Bifurcation and nonlinear dynamic analysis of a rigid rotor supported by two-lobe noncircular gas-lubricated journal bearing system | |
Pan et al. | Dynamic modeling and response analysis of rub-impact rotor system with squeeze film damper under maneuvering load | |
Wang et al. | Bifurcation and nonlinear dynamic analysis of a flexible rotor supported by relative short gas journal bearings | |
Liu et al. | Nonlinear vibration response of a complex aeroengine under the rubbing fault | |
Zhou et al. | Numerical analysis of a nonlinear double disc rotor-seal system | |
Chen et al. | Dynamic behavior of a flexible rotor system with squeeze film damper considering oil-film inertia under base motions | |
Sun et al. | Dynamic analysis of the variable stiffness support rotor system with elastic rings | |
CN105550395B (zh) | 透平机械单支撑轴系几何对中安装方法 | |
Hu et al. | An enhanced axisymmetric solid element for rotor dynamic model improvement | |
CN116090135B (zh) | 一种阻尼器-转子系统响应分析方法 | |
CN114021270A (zh) | 一种不对中球轴承支撑力计算方法 | |
Kim et al. | Mesoscale foil gas bearings for palm-sized turbomachinery: design, manufacturing, and modeling | |
Wang et al. | Effects of unbalance orientation on the dynamic characteristics of a double overhung rotor system for high-speed turbochargers | |
Al-Khazali et al. | The experimental analysis of vibration monitoring in system rotor dynamic with validate results using simulation data | |
CN115270342B (zh) | 一种剥落故障深沟球轴承多自由度动力学建模方法 | |
CN115688387A (zh) | 一种基于轴承寿命的支点阻尼器油膜间隙选取方法 | |
Ma et al. | Numerical research on rub-impact fault in a blade-rotor-casing coupling system | |
Pan et al. | Nonlinear vibration behavior of aero-engine rotor-bearing system in maneuvering flight | |
Taplak et al. | Passive balancing of a rotating mechanical system using genetic algorithm | |
Miranda et al. | Finite element method applied to the eigenvalue analysis of flexible rotors supported by journal bearings |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |