CN103471854B - 一种航空发动机整机振动特性分析方法 - Google Patents
一种航空发动机整机振动特性分析方法 Download PDFInfo
- Publication number
- CN103471854B CN103471854B CN201310445186.7A CN201310445186A CN103471854B CN 103471854 B CN103471854 B CN 103471854B CN 201310445186 A CN201310445186 A CN 201310445186A CN 103471854 B CN103471854 B CN 103471854B
- Authority
- CN
- China
- Prior art keywords
- rotor
- stator
- engine
- analysis
- complete machine
- 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
Landscapes
- Testing Of Devices, Machine Parts, Or Other Structures Thereof (AREA)
Abstract
一种航空发动机整机振动特性分析方法,将发动机整机***隔离为转子***、静子***,通过支点***将二者耦合起来,采用转、静子***动力学特性联合仿真;将发动机整机***隔离为转子—支承***和静子承力***两个研究对象,然后应用有限元分析平台仿真分析转子—支承***在不平衡激振力和碰摩力作用下的振动响应及动态轴承支反力,得到静子***的激励输入;最后,将得到的动态轴承支反力和碰摩力加载于静子承力***,仿真分析其振动响应特性,从而得到发动机整机振动特性。本发明的优点:减少航空发动机设计中的故障隐患,并为出现的发动机振动故障提供决策依据和技术支撑。
Description
技术领域
本发明涉及航空发动机整机振动分析领域,特别涉及了一种航空发动机整机振动特性分析方法。
背景技术
目前对于发动机整机振动特性研究主要局限于转子***独立的特性分析,得到的结果与发动机实际工况相符性差,对于转、静子***特性联合研究相对较少,对工程的指导性不强。
随着控制航空发动机振动故障的需求迫切性越来越强,对分析结果的准确性、真实性和全面性要求越来越高,迫切需要针对发动机整机***进行全三维转静子动力学特性联合仿真分析,指导发动机的设计、装配和排振,保证发动机工作的可靠性,减少振动故障的发生。
发明内容
本发明的目的是为发动机生产中振动故障的控制和外场使用可靠性提供坚实的技术保障,特提供了一种航空发动机整机振动特性分析方法。
本发明提供了一种航空发动机整机振动特性分析方法,其特征在于:将发动机整机***隔离为转子***、静子***,通过支点***将二者耦合起来,采用转、静子***动力学特性联合仿真;
首先对发动机整机***进行力学分析,转子***是振动源,静子***为响应体,转子振动能量以支反力的形式通过支点传递到静子***,静子***的响应就可反映整机振动特性。把一个复杂的几乎无法分析的庞大***分隔为可分析的子***,并通过动力学特性联合仿真将各***有机的结合在一起,真实反映整机的实际状态,实现发动机整机振动特性的分析;
发动机整机包括转子—支承***和静子承力***。转子在高速旋转时,产生激振力,转子和静子承力***通过支承轴承相互耦合,转子产生的激振力以轴承支反力的形式作用于静子承力***。此外,在某些特殊工况下,转静件会发生碰摩,产生碰摩力。
首先,将发动机整机***隔离为转子—支承***和静子承力***两个研究对象,然后应用有限元分析平台仿真分析转子—支承***在不平衡激振力和碰摩力作用下的振动响应及动态轴承支反力,得到静子***的激励输入。最后,将得到的动态轴承支反力和碰摩力加载于静子承力***,仿真分析其振动响应特性,从而得到发动机整机振动特性。
航空发动机整机振动特性分析方法,具体步骤:
1)首先应用三维CAD软件建立转子—支承***和静子承力***的三维几何模型,然后利用有限元工具对模型进行网格划分、零件连接方式处理、边界条件施加等,建立转子—支承***、静子承力***的有限元模型。
2)应用多自由度强迫振动理论,对转子—支承***有限元模型进行验证;进行静子承力***模态实验,根据模态实验结果对有限元模型进行验证和修正。
3)建立高压涡轮-静子机匣的局部碰摩模型,进行非线性接触分析,得到碰摩力随径向相对位移的变化关系(图2)。
4)应用有限元软件的转子动力学计算分析功能,仿真计算转子—支承***在不平衡量作用下的振动响应。根据发动机试车中的典型故障,调整不平衡量大小、分布以及有无碰摩的组合工况,得到不同工况下的动态轴承支反力。
5)将不同工况下的动态轴承支反力和碰摩力施加在静子承力***上,进行静子机匣的振动响应特性计算。图3显示了标准工况下某一振动测点(0~4)阶范围的振动谱图阶次,1阶对应高压转子模型的频率振动,0.61阶对应低压转子模型的频率振动。图4是仿真计算得到的对应某一振动测点中水速度频谱图,图中显示在45Hz、75Hz、109Hz、148Hz、245Hz、300Hz频率点出现了较高的响应峰值,其中148Hz、245Hz分别对应低压转子和高压转子的工作频率,而45Hz、75Hz、109Hz、300Hz是静子机匣自由振动的固有频率。测试频谱图与计算频谱图在148Hz、245Hz、300Hz频率点的响应峰值有很好的吻合性。
6)将静子机匣的振动响应仿真结果与典型发动机试车中响应的实测数据进行对比分析,以验证所建的模型和计算方法的正确性。进一步分析机匣关键测点的振动响应值对不平衡量大小、分布及碰摩的敏感程度,得到发动机的整机振动特性。
整机振动故障严重困扰发动机生产及外场使用安全,必须对发动机整机振动特性和机理进行深入研究,而传统分析方法无法真实模拟发动机实际结构和工况,分析结果的参考性差。而整机全三维仿真则使模型过于庞大,无法实现计算。因此,将发动机整机***隔离为转子***、静子***,通过支点***将二者耦合起来,采用转、静子***动力学特性联合仿真,解决该问题。
本发明的优点:
目前航空发动机整机振动特性分析基本都是基于高度简化的一维或二维转子动力学模型,随着航空发动机技术的发展,对整机振动特性分析的需求也越来越强,传统的分析方法局限性越来越明显。开展全三维转、静子动力学特性联合仿真分析进行整机振动特性研究可有效解决该问题,具有十分重大的意义,本发明的实施有助于提升航空发动机从业人员的技术素质,减少航空发动机设计中的故障隐患,并为出现的发动机振动故障提供决策依据和技术支撑。
附图说明
下面结合附图及实施方式对本发明作进一步详细的说明:
图1为技术方案流程简图;
图2为碰摩力随径向相对位移的变化关系图;
图3为振动测点的试验速度频谱图;
图4为振动测点的仿真计算速度频谱图。
具体实施方式
实施例
本实施例提供了一种航空发动机整机振动特性分析方法,其特征在于:将发动机整机***隔离为转子***、静子***,通过支点***将二者耦合起来,采用转、静子***动力学特性联合仿真;
首先对发动机整机***进行力学分析,转子***是振动源,静子***为响应体,转子振动能量以支反力的形式通过支点传递到静子***,静子***的响应就可反映整机振动特性。把一个复杂的几乎无法分析的庞大***分隔为可分析的子***,并通过动力学特性联合仿真将各***有机的结合在一起,真实反映整机的实际状态,实现发动机整机振动特性的分析;
发动机整机包括转子—支承***和静子承力***。转子在高速旋转时,产生激振力,转子和静子承力***通过支承轴承相互耦合,转子产生的激振力以轴承支反力的形式作用于静子承力***。此外,在某些特殊工况下,转静件会发生碰摩,产生碰摩力。
首先,将发动机整机***隔离为转子—支承***和静子承力***两个研究对象,然后应用有限元分析平台仿真分析转子—支承***在不平衡激振力和碰摩力作用下的振动响应及动态轴承支反力,得到静子***的激励输入。最后,将得到的动态轴承支反力和碰摩力加载于静子承力***,仿真分析其振动响应特性,从而得到发动机整机振动特性。
航空发动机整机振动特性分析方法,具体步骤:
1)首先应用三维CAD软件建立转子—支承***和静子承力***的三维几何模型,然后利用有限元工具对模型进行网格划分、零件连接方式处理、边界条件施加等,建立转子—支承***、静子承力***的有限元模型。
2)应用多自由度强迫振动理论,对转子—支承***有限元模型进行验证;进行静子承力***模态实验,根据模态实验结果对有限元模型进行验证和修正。
3)建立高压涡轮-静子机匣的局部碰摩模型,进行非线性接触分析,得到碰摩力随径向相对位移的变化关系(图2)。
4)应用有限元软件的转子动力学计算分析功能,仿真计算转子—支承***在不平衡量作用下的振动响应。根据发动机试车中的典型故障,调整不平衡量大小、分布以及有无碰摩的组合工况,得到不同工况下的动态轴承支反力。
5)将不同工况下的动态轴承支反力和碰摩力施加在静子承力***上,进行静子机匣的振动响应特性计算。图3显示了标准工况下某一振动测点(0~4)阶范围的振动谱图阶次,1阶对应高压转子模型的频率振动,0.61阶对应低压转子模型的频率振动。图4是仿真计算得到的对应某一振动测点中水速度频谱图,图中显示在45Hz、75Hz、109Hz、148Hz、245Hz、300Hz频率点出现了较高的响应峰值,其中148Hz、245Hz分别对应低压转子和高压转子的工作频率,而45Hz、75Hz、109Hz、300Hz是静子机匣自由振动的固有频率。测试频谱图与计算频谱图在148Hz、245Hz、300Hz频率点的响应峰值有很好的吻合性。
6)将静子机匣的振动响应仿真结果与典型发动机试车中响应的实测数据进行对比分析,以验证所建的模型和计算方法的正确性。进一步分析机匣关键测点的振动响应值对不平衡量大小、分布及碰摩的敏感程度,得到发动机的整机振动特性。
整机振动故障严重困扰发动机生产及外场使用安全,必须对发动机整机振动特性和机理进行深入研究,而传统分析方法无法真实模拟发动机实际结构和工况,分析结果的参考性差。而整机全三维仿真则使模型过于庞大,无法实现计算。因此,将发动机整机***隔离为转子***、静子***,通过支点***将二者耦合起来,采用转、静子***动力学特性联合仿真,解决该问题。
Claims (2)
1.一种航空发动机整机振动特性分析方法,其特征在于:将发动机整机***隔离为转子***、静子***,通过支点***将二者耦合起来,采用转、静子***动力学特性联合仿真;
首先对发动机整机***进行力学分析,转子***是振动源,静子***为响应体,转子振动能量以支反力的形式通过支点传递到静子***,静子***的响应就可反映整机振动特性;把一个复杂的几乎无法分析的庞大***分隔为可分析的子***,并通过动力学特性联合仿真将各***有机的结合在一起,真实反映整机的实际状态,实现发动机整机振动特性的分析;
发动机整机包括转子***和静子***;转子在高速旋转时,产生激振力,转子和静子***通过支承轴承相互耦合,转子产生的激振力以轴承支反力的形式作用于静子***;此外,在某些特殊工况下,转静件会发生碰摩,产生碰摩力;
首先,将发动机整机***隔离为转子***和静子***两个研究对象,然后应用有限元分析平台仿真分析转子***在不平衡激振力和碰摩力作用下的振动响应及动态轴承支反力,得到静子***的激励输入;最后,将得到的动态轴承支反力和碰摩力加载于静子***,仿真分析其振动响应特性,从而得到发动机整机振动特性。
2.按照权利要求1所述的航空发动机整机振动特性分析方法,其特征在于:航空发动机整机振动特性分析方法,具体步骤:
1)首先应用三维CAD软件建立转子***和静子***的三维几何模型,然后利用有限元工具对模型进行网格划分、零件连接方式处理、边界条件施加,建立转子***、静子***的有限元模型;
2)应用多自由度强迫振动理论,对转子***有限元模型进行验证;进行静子***模态实验,根据模态实验结果对静子***有限元模型进行验证和修正;
3)建立高压涡轮-静子机匣的局部碰摩模型,进行非线性接触分析,得到碰摩力随径向相对位移的变化关系;
4)应用有限元软件的转子动力学计算分析功能,仿真计算转子***在不平衡量作用下的振动响应;根据发动机试车中的典型故障,调整不平衡量大小、分布以及有无碰摩的组合工况,得到不同工况下的动态轴承支反力;
5)将不同工况下的动态轴承支反力和碰摩力施加在静子***上,进行静子机匣的振动响应特性计算;标准工况下某一振动测点0~4阶范围的振动谱图阶次,1阶对应高压转子模型的频率振动,0.61阶对应低压转子模型的频率振动;
6)将静子机匣的振动响应仿真结果与典型发动机试车中响应的实测数据进行对比分析,以验证所建的模型和计算方法的正确性;进一步分析机匣关键测点的振动响应值对不平衡量大小、分布及碰摩的敏感程度,得到发动机的整机振动特性。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201310445186.7A CN103471854B (zh) | 2013-09-26 | 2013-09-26 | 一种航空发动机整机振动特性分析方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201310445186.7A CN103471854B (zh) | 2013-09-26 | 2013-09-26 | 一种航空发动机整机振动特性分析方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN103471854A CN103471854A (zh) | 2013-12-25 |
CN103471854B true CN103471854B (zh) | 2016-04-06 |
Family
ID=49796796
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201310445186.7A Active CN103471854B (zh) | 2013-09-26 | 2013-09-26 | 一种航空发动机整机振动特性分析方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN103471854B (zh) |
Families Citing this family (20)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103912315B (zh) * | 2014-04-14 | 2015-06-17 | 西北工业大学 | 一种航空发动机转子结构动力学的设计方法 |
CN104634445B (zh) * | 2015-01-27 | 2017-06-23 | 西安交通大学 | 一种双转子航空发动机时变振动贡献量检测方法及*** |
CN106649902B (zh) * | 2015-10-30 | 2019-08-23 | 湖北航天技术研究院总体设计所 | 一种飞行器舱段固有振动特性快速评估及优化方法 |
CN105512409B (zh) * | 2015-12-11 | 2018-08-28 | 中国北方发动机研究所(天津) | 一种发动机零部件优化设计方法 |
CN105571867B (zh) * | 2016-01-21 | 2018-10-02 | 苏州东菱振动试验仪器有限公司 | 一种航空发动机振动控制实验平台 |
CN106295074B (zh) * | 2016-08-30 | 2019-02-22 | 湖北航天技术研究院总体设计所 | 一种运载火箭舱段振动响应特性快速分析及优化方法 |
CN106777518B (zh) * | 2016-11-24 | 2019-10-18 | 北京电子工程总体研究所 | 一种针对姿轨控发动机推力的激励源辨识方法 |
CN107844673A (zh) * | 2017-12-14 | 2018-03-27 | 中国航发沈阳发动机研究所 | 一种航空发动机整机三维气动仿真方法 |
CN109117506B (zh) * | 2018-07-12 | 2020-07-17 | 北京航空航天大学 | 一种针对通用叶盘结构的嵌入式压电分支阻尼优化设计方法 |
CN109374300B (zh) * | 2018-12-13 | 2020-09-25 | 西安理工大学 | 一种基于ump的碰摩转子-轴承***振动特性分析方法 |
CN109684711B (zh) * | 2018-12-20 | 2023-09-22 | 中国人民解放军海军航空大学 | 一种涡轴发动机气动连接多转子耦合振动分析方法 |
CN110175392B (zh) * | 2019-05-22 | 2022-12-13 | 中国民航大学 | 基于otpa法和实物建模的航空发动机振动传递路径分析方法 |
CN110955987B (zh) * | 2019-06-20 | 2023-02-21 | 中南大学 | 一种基于***动力学的航空机轮振动分析方法 |
CN110543694A (zh) * | 2019-08-13 | 2019-12-06 | 中国航发贵阳发动机设计研究所 | 一种航空发动机辅助吊挂拉杆振动计算方法 |
CN111473859B (zh) * | 2020-05-21 | 2022-07-29 | 中国航发湖南动力机械研究所 | 整机振动限制值的制定方法 |
CN112525533A (zh) * | 2020-10-30 | 2021-03-19 | 中国航发沈阳黎明航空发动机有限责任公司 | 一种航空发动机滚珠轴承接触角的在线检测方法 |
CN112417620B (zh) * | 2020-11-26 | 2024-03-01 | 西安长庆科技工程有限责任公司 | 一种橇装式柱塞泵一体化集成装置的振动仿真分析方法 |
CN115077919B (zh) * | 2022-05-28 | 2024-04-12 | 西北工业大学 | 一种适于航空发动机机载的整机振动评估方法 |
CN116306143B (zh) * | 2023-03-17 | 2023-09-22 | 中国航空发动机研究院 | 双转子***不平衡故障仿真方法及电子设备 |
CN116907863B (zh) * | 2023-09-12 | 2023-11-28 | 中国航发沈阳发动机研究所 | 一种航空发动机振动特征变化与构件损伤监视和预判方法 |
Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN201277910Y (zh) * | 2008-10-22 | 2009-07-22 | 沈阳黎明航空发动机(集团)有限责任公司 | 一种航空发动机试车参数测量*** |
CN101599104A (zh) * | 2009-07-16 | 2009-12-09 | 北京航空航天大学 | 一种航空涡轮发动机叶片颤振边界的模拟方法 |
CN102507205A (zh) * | 2011-10-27 | 2012-06-20 | 西北工业大学 | 一种检测航空发动机风扇叶片颤振故障的方法 |
CN102928232A (zh) * | 2012-11-21 | 2013-02-13 | 中国民用航空飞行学院 | 一种航空发动机整机性能衰退趋势预测方法 |
CN103091108A (zh) * | 2013-01-17 | 2013-05-08 | 商旭升 | 一种用于航空发动机的振动检测*** |
CN103149029A (zh) * | 2013-01-16 | 2013-06-12 | 南京航空航天大学 | 利用倒频谱识别航空发动机转静碰摩部位的方法 |
-
2013
- 2013-09-26 CN CN201310445186.7A patent/CN103471854B/zh active Active
Patent Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN201277910Y (zh) * | 2008-10-22 | 2009-07-22 | 沈阳黎明航空发动机(集团)有限责任公司 | 一种航空发动机试车参数测量*** |
CN101599104A (zh) * | 2009-07-16 | 2009-12-09 | 北京航空航天大学 | 一种航空涡轮发动机叶片颤振边界的模拟方法 |
CN102507205A (zh) * | 2011-10-27 | 2012-06-20 | 西北工业大学 | 一种检测航空发动机风扇叶片颤振故障的方法 |
CN102928232A (zh) * | 2012-11-21 | 2013-02-13 | 中国民用航空飞行学院 | 一种航空发动机整机性能衰退趋势预测方法 |
CN103149029A (zh) * | 2013-01-16 | 2013-06-12 | 南京航空航天大学 | 利用倒频谱识别航空发动机转静碰摩部位的方法 |
CN103091108A (zh) * | 2013-01-17 | 2013-05-08 | 商旭升 | 一种用于航空发动机的振动检测*** |
Non-Patent Citations (2)
Title |
---|
航空发动机振动监测研究;孙海东等;《机械设计与制造》;20070228(第2期);第127-129页 * |
航空发动机整机有限元模型转子动力学分析;陈萌等;《北京航空航天大学学报》;20070930;第33卷(第9期);第1013-1016页 * |
Also Published As
Publication number | Publication date |
---|---|
CN103471854A (zh) | 2013-12-25 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN103471854B (zh) | 一种航空发动机整机振动特性分析方法 | |
Zeng et al. | The generalized Hamiltonian model for the shafting transient analysis of the hydro turbine generating sets | |
CN102564750B (zh) | 一种水平轴风力机叶片模态测试结果修正方法 | |
Meng et al. | An investigation on dynamic characteristics of a gas turbine rotor using an improved transfer matrix method | |
Wang et al. | Nonlinear dynamic response and stability of a rod fastening rotor with internal damping effect | |
Li et al. | Rotating blade-casing rubbing simulation considering casing flexibility | |
CN103942091A (zh) | Matlab自定义模型和psasp联合仿真的励磁***仿真方法及*** | |
CN102915388A (zh) | 一种面向对象的转子动力学***非线性非因果建模仿真方法 | |
Ma et al. | Unbalance vibration characteristics and sensitivity analysis of the dual-rotor system in aeroengines | |
CN103207924A (zh) | 风力发电机组机械暂态与电气暂态混合仿真***及方法 | |
CN104992039A (zh) | 一种圆锥面油压过盈安装过程中压力的计算方法 | |
Yin et al. | Reproduction of five degree-of-freedom loads for wind turbine using equispaced electro-hydraulic actuators | |
CN102629283A (zh) | 一种转动部件对挠性动力学影响的仿真分析方法 | |
CN104807606A (zh) | 一种matlab-stm32混合动力试验***及其试验方法 | |
Senousy et al. | Electromagnetics, structural harmonics and acoustics coupled simulation on the stator of an electric motor | |
CN103064404B (zh) | 汽车起重机的功率匹配控制仿真测试*** | |
CN104361172A (zh) | 次同步谐振的机电联合仿真方法及*** | |
CN104570767A (zh) | 一种磁流变减振器控制算法仿真试验*** | |
CN103177168A (zh) | 一种齿轮齿面摩擦因数的反求方法 | |
Feng et al. | Research on hardware-in-the-loop simulation system of hydraulic excavators | |
Xinhua et al. | Development of a rapid design system for aerial work truck subframe with UG secondary development framework | |
CN204331459U (zh) | 磁流变减振器控制算法仿真试验*** | |
Dong et al. | A precise FE model of a spur gear set considering eccentric error for quasi-static analysis | |
Qu et al. | Effect of the AERO-engine mounting stiffness on the whole engine coupling vibration | |
Hao et al. | Rigid-flexible coupling dynamics simulation of planetary gear transmission based on MFBD |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant | ||
CP01 | Change in the name or title of a patent holder |
Address after: 110043 Dong TA street, Dadong District, Shenyang, Liaoning Province, No. 6 Patentee after: Chinese Hangfa Shenyang Liming Aero engine limited liability company Address before: 110043 Dong TA street, Dadong District, Shenyang, Liaoning Province, No. 6 Patentee before: Liming Aeroplane Engine (Group) Co., Ltd., Shenyang City |
|
CP01 | Change in the name or title of a patent holder |