CN115879247B - 一种基于系统辨识的轮盘关键部位应力计算方法 - Google Patents
一种基于系统辨识的轮盘关键部位应力计算方法 Download PDFInfo
- Publication number
- CN115879247B CN115879247B CN202310188132.0A CN202310188132A CN115879247B CN 115879247 B CN115879247 B CN 115879247B CN 202310188132 A CN202310188132 A CN 202310188132A CN 115879247 B CN115879247 B CN 115879247B
- Authority
- CN
- China
- Prior art keywords
- stress
- wheel disc
- system identification
- parameters
- load spectrum
- 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
Images
Classifications
-
- 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
- Testing Of Devices, Machine Parts, Or Other Structures Thereof (AREA)
Abstract
本发明提供了一种基于系统辨识的轮盘关键部位应力计算方法,包括以下步骤:获取航空发动机的载荷谱参数以及轮盘观测应力;对载荷谱参数数据和轮盘观测应力数据分组,形成训练样本和验证样本;选择初始模型,并基于训练样本对初始模型训练,构建基于系统辨识的应力计算模型;采用基于系统辨识的应力计算模型对轮盘关键部位应力计算。本发明采用构建的基于系统辨识的应力模型进行计算轮盘关键部位应力,能够快速得到轮盘关键部位的应力历程,进而开展轮盘寿命消耗分析,满足轮盘寿命监视的需求。
Description
技术领域
本发明涉及航空发动机轮盘强度设计技术领域,具体为一种基于系统辨识的轮盘关键部位应力计算方法。
背景技术
发动机的轮盘作为限寿件,其一旦发生失效可能导致危害性后果,一般需要在规定使用周期内进行检查,且必须在所有可用寿命消耗完之前进行更换,为了保证安全性,先进的航空发动机一般都配有寿命监视系统。
航空发动机健康管理系统中寿命监视系统是机载系统的一部分,要求通过跟踪时变的压力、温度和转子转速的变化累积使用寿命,并驱动定寿算法即关键部件的物理失效实时模型对限寿件的寿命进行监视。
目前,限寿件的寿命通常采用有限元方法计算其关键部位的应力值后分析获得,但是该方法的计算过程较为复杂、耗时较长,不能满足寿命监视对于快速计算的需求。
鉴于此,提出了一种基于系统辨识的轮盘关键部位应力计算方法,进而对发动机中某些限寿件的寿命监视。
发明内容
本发明的目的在于提供一种简单、高效的基于系统辨识的轮盘关键部位应力计算方法。
实现发明目的的技术方案如下:一种基于系统辨识的轮盘关键部位应力计算方法,包括以下步骤:
S1、获取航空发动机的载荷谱参数以及轮盘观测应力;
S2、对载荷谱参数数据和轮盘观测应力数据分组,形成训练样本和验证样本;
S3、选择初始模型,并基于训练样本对初始模型训练,构建基于系统辨识的应力计算模型;
S4、采用基于系统辨识的应力计算模型对轮盘关键部位应力计算。
在一个实施例中,步骤S1中,航空发动机载荷谱参数和轮盘观测应力的获取方法,包括:
S11、基于试验方法或仿真分析方法,获取航空发动机各时间历程下不同载荷谱的初始参数;
S12、对发动机轮盘的引气流路进行物理特性及相关性分析,从初始参数筛选影响发动机轮盘关键部位应力的主要参数;
S13、对主要参数归一化处理,获得各时间历程下的载荷谱参数;
S14、依据有限元仿真分析方法,对各时间历程下载荷谱参数进行温度场分析和应力分析,获取各时间历程下的轮盘观测应力。
在一个实施例中,步骤中S12中,主要参数包括发动机的高度、马赫数、转速、油门杆角度、各截面温度力、各截面压力,以及喷管喉道面积反馈值。
在一个实施例中,步骤S2中,训练样本包括起飞、中间、停车中任意一个时间历程或多个时间历程的载荷谱参数和轮盘观测应力;
所述验证样本包括由起飞、中间、停车形成的全时间历程的载荷谱参数和轮盘观测应力。
在一个实施例中,步骤S3中,选择初始模型,并基于训练样本对初始模型训练,构建基于系统辨识的应力计算模型,包括:
S31、基于轮盘观测应力-时间历程的非线性特征,选择初始模型为非线性ARX模型;
S32、用训练样本对非线性ARX模型训练,获取非线性ARX模型的模型阶次和非线性算法,构建基于系统辨识的应力计算模型。
在一个实施例中,步骤S32中,所述非线性算法为小波网络算法,所述模型阶次为4阶,即前4个应力点对当前应力有影响,每个载荷谱参数的前4个数据点对当前应力有影响,不考虑样本点对应力的延迟影响。
在一个实施例中,轮盘关键部位包括轮盘中心孔、偏心孔、倒圆、榫槽中的任意一种。
在一个实施例中,步骤S3中,还包括:
采用验证样本,对构建的基于系统辨识的应力计算模型进行验证。
与现有技术相比,本发明的有益效果是:本发明设计的基于系统辨识的轮盘关键部位应力计算方法,能够快速得到轮盘关键部位的应力历程,进而开展轮盘寿命消耗分析,满足轮盘寿命监视的需求。
附图说明
为了更清楚地说明本发明实施例技术方案,下面将对实施例描述中所需要使用的附图作简单地介绍。
图1为本发明基于系统辨识的轮盘关键部位应力计算方法的流程图;
图2为具体实施方式中构建基于系统辨识的应力计算模型的流程图。
具体实施方式
下面结合具体实施例来进一步描述本发明,本发明的优点和特点将会随着描述而更为清楚。但这些实施例仅是范例性的,并不对本发明的范围构成任何限制。本领域技术人员应该理解的是,在不偏离本发明的精神和范围下可以对本发明技术方案的细节和形式进行修改或替换,但这些修改和替换均落入本发明的保护范围内。
本具体实施方式提供了一种基于系统辨识的轮盘关键部位应力计算方法,参见图1所示,轮盘关键部位应力计算方法包括以下步骤:
S1、获取航空发动机的载荷谱参数以及轮盘观测应力。
本步骤中,航空发动机载荷谱参数和轮盘观测应力的获取方法,包括:
S11、基于试验方法或仿真分析方法,获取航空发动机各时间历程下不同载荷谱的初始参数;
其中,初始参数包括发动机的高度、马赫数、转速、油门杆角度、各截面温度力、各截面压力,以及喷管喉道面积反馈值,还包括加力接通状态、机动过载参数、机动角速度等。
S12、对发动机轮盘的引气流路进行物理特性及相关性分析,从初始参数筛选影响发动机轮盘关键部位应力的主要参数;
其中,主要参数包括发动机的高度、马赫数、转速、油门杆角度、各截面温度力、各截面压力,以及喷管喉道面积反馈值。
S13、对主要参数归一化处理,获得各时间历程下的载荷谱参数;
其中,主要参数采用现有已知的归一化处理方法进行处理。
S14、依据有限元仿真分析方法,对各时间历程下载荷谱参数进行温度场分析和应力分析,获取各时间历程下的轮盘观测应力。
S2、对载荷谱参数数据和轮盘观测应力数据分组,形成训练样本和验证样本。
本步骤中,训练样本包括起飞、中间、停车中任意一个时间历程或多个时间历程的载荷谱参数和轮盘观测应力;
所述验证样本包括由起飞、中间、停车形成的全时间历程的载荷谱参数和轮盘观测应力。
S3、选择初始模型,并基于训练样本对初始模型训练,构建基于系统辨识的应力计算模型。
本步骤中,选择初始模型,并基于训练样本对初始模型训练,构建基于系统辨识的应力计算模型,包括:
S31、基于轮盘观测应力-时间历程的非线性特征,选择初始模型为非线性ARX模型;
S32、用训练样本对非线性ARX模型训练,获取非线性ARX模型的模型阶次和非线性算法,构建基于系统辨识的应力计算模型。
在一个改进实施例中,步骤S32中,所述非线性算法为小波网络算法,所述模型阶次为4阶,即前4个应力点对当前应力有影响,每个载荷谱参数的前4个数据点对当前应力有影响,不考虑样本点对应力的延迟影响。
在另一个改进实施例中,还包括对经过步骤S31~S32构建的基于系统辨识的应力计算模型进行验证评估:
S33、采用验证样本,对构建的基于系统辨识的应力计算模型进行验证;
本步骤中,基于系统辨识的应力计算模型验证使用包括按照输入参数分阶段比对模型,下一阶段的初始值使用前一阶段的终值以保证连续性,输出应力值,将输出的应力值与观测应力比较,两者误差需在±5%以内。
S4、采用基于系统辨识的应力计算模型对轮盘关键部位应力计算。
本具体实施方式中,上述轮盘关键部位应力计算方法能够对轮盘关键部位中包括但不限于轮盘中心孔、偏心孔、倒圆、榫槽的应力进行计算。
本发明提供的基于系统辨识的轮盘关键部位应力计算方法,具有以下优点:能够快速得到轮盘关键部位的应力历程,进而开展轮盘寿命消耗分析,满足轮盘寿命监视的需求。
本具体实施方式结合附图2对上述基于系统辨识的应力计算模型的构建方法进行说明:
1.通过试验方法或仿真分析方法,获取航空发动机起飞、中间、停车三个时间历程下不同载荷谱的初始参数。
初始参数包括发动机的高度、马赫数、转速、油门杆角度、发动机各截面的温度、各截面压力、喷管喉道面积反馈值、加力接通状态、机动过载参数、机动角速度等共14个参数,初始参数参加下表1所示:
表1:初始参数
2.对上述初始参数进行物理特性及相关性分析,从初始参数筛选影响发动机轮盘关键部位应力的主要参数。
主要参数包括发动机的高度、马赫数、转速、油门杆角度、各截面温度力、各截面压力,以及喷管喉道面积反馈值等共8个参数,具体参见下表2所示:
表2:主要参数
3.对上述主要参数进行归一化处理,得到航空发动机起飞、中间、停车三个时间历程下的载荷谱参数,作为构建基于系统辨识的应力计算模型的样本自变量;
依据有限元仿真分析方法,对航空发动机起飞、中间、停车三个时间历程下的载荷谱参数进行温度场分析和应力分析,获取各时间历程下的轮盘观测应力,作为构建基于系统辨识的应力计算模型的样本因变量。
4.对航空发动机起飞、中间、停车三个时间历程下的载荷谱参数和轮盘观测应力进行分组,形成训练样本和验证样本;
其中训练样本包括起飞、中间、停车中任意一个时间历程或多个时间历程的载荷谱参数和轮盘观测应力;
所述验证样本包括由起飞、中间、停车形成的全时间历程的载荷谱参数和轮盘观测应力。
5.根据轮盘观测应力-时间历程的非线性特征,选择非线性ARX模型作为初始模型,用训练样本对其进行训练,用验证样本对训练后的模型进行验证优化,得到最终的模型阶次和非线性算法,完成基于系统辨识的应力计算模型的构建;
具体的,模型阶次为4阶,即前4个应力点及每个载荷谱参数的前4个数据点对当前应力有影响,不考虑样本点对应力的延迟影响。非线性算法选择为小波网络;
在此需要说明的是,本具体实施方式中不对模型阶次和非线性算法进行限定,例如模型阶次也可以选择使用4阶以外的其他阶次,非线性算法也可以使用其他算法。
6.发动机工作时,获取起飞时间历程或中间时间历程或停车时间历程的用于轮盘中心孔应力计算的载荷谱参数,输入模型中后输出获得轮盘中心孔应力。
以上所述仅为本发明的较佳实施例而已,并不用以限制本发明,凡在本发明的精神和原则之内,所作的任何修改、等同替换、改进等,均应包含在本发明的保护范围之内。
此外,应当理解,虽然本说明书按照实施方式加以描述,但并非每个实施方式仅包含一个独立的技术方案,说明书的这种叙述方式仅仅是为清楚起见,本领域技术人员应当将说明书作为一个整体,各实施例中的技术方案也可以经适当组合,形成本领域技术人员可以理解的其他实施方式。
Claims (7)
1.一种基于系统辨识的轮盘关键部位应力计算方法,其特征在于,包括以下步骤:
S1、获取航空发动机的载荷谱参数以及轮盘观测应力;
S2、对载荷谱参数数据和轮盘观测应力数据分组,形成训练样本和验证样本;
S31、基于轮盘观测应力-时间历程的非线性特征,选择初始模型为非线性ARX模型;
S32、用训练样本对非线性ARX模型训练,获取非线性ARX模型的模型阶次和非线性算法,构建基于系统辨识的应力计算模型;
S4、采用基于系统辨识的应力计算模型对轮盘关键部位应力计算。
2.根据权利要求1所述的轮盘关键部位应力计算方法,其特征在于,步骤S1中,航空发动机载荷谱参数和轮盘观测应力的获取方法,包括:
S11、基于试验方法或仿真分析方法,获取航空发动机各时间历程下不同载荷谱的初始参数;
S12、对发动机轮盘的引气流路进行物理特性及相关性分析,从初始参数筛选影响发动机轮盘关键部位应力的主要参数;
S13、对主要参数归一化处理,获得各时间历程下的载荷谱参数;
S14、依据有限元仿真分析方法,对各时间历程下载荷谱参数进行温度场分析和应力分析,获取各时间历程下的轮盘观测应力。
3.根据权利要求2所述的轮盘关键部位应力计算方法,其特征在于,步骤中S12中,主要参数包括发动机的高度、马赫数、转速、油门杆角度、各截面温度力、各截面压力,以及喷管喉道面积反馈值。
4.根据权利要求2所述的轮盘关键部位应力计算方法,其特征在于,步骤S2中,训练样本包括起飞、中间、停车中任意一个时间历程或多个时间历程的载荷谱参数和轮盘观测应力;
所述验证样本包括由起飞、中间、停车形成的全时间历程的载荷谱参数和轮盘观测应力。
5.根据权利要求1所述的轮盘关键部位应力计算方法,其特征在于,步骤S32中,所述非线性算法为小波网络算法,所述模型阶次为4阶。
6.根据权利要求1~5任一项所述的轮盘关键部位应力计算方法,其特征在于,轮盘关键部位包括轮盘中心孔、偏心孔、倒圆、榫槽中的任意一种。
7.根据权利要求6所述的轮盘关键部位应力计算方法,其特征在于,步骤S3中,还包括:
采用验证样本,对构建的基于系统辨识的应力计算模型进行验证。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202310188132.0A CN115879247B (zh) | 2023-03-02 | 2023-03-02 | 一种基于系统辨识的轮盘关键部位应力计算方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202310188132.0A CN115879247B (zh) | 2023-03-02 | 2023-03-02 | 一种基于系统辨识的轮盘关键部位应力计算方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN115879247A CN115879247A (zh) | 2023-03-31 |
CN115879247B true CN115879247B (zh) | 2023-06-30 |
Family
ID=85761770
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202310188132.0A Active CN115879247B (zh) | 2023-03-02 | 2023-03-02 | 一种基于系统辨识的轮盘关键部位应力计算方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN115879247B (zh) |
Families Citing this family (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN116542114B (zh) * | 2023-07-04 | 2023-09-19 | 中国航发四川燃气涡轮研究院 | 基于总体参数的压气机转子温度应力变形分析方法及装置 |
CN117725802B (zh) * | 2024-02-07 | 2024-04-16 | 中国航发四川燃气涡轮研究院 | 航空发动机主轴疲劳试验标准循环载荷谱构建方法及系统 |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102411641A (zh) * | 2010-09-25 | 2012-04-11 | 集美大学 | 一种基于数值模拟的汽车覆盖件回弹误差补偿方法 |
CN110697438A (zh) * | 2017-09-19 | 2020-01-17 | 中国计量大学 | 基于神经网络的直落失重式物料下料机控制器 |
CN111859729A (zh) * | 2020-06-04 | 2020-10-30 | 北京航空航天大学 | 考虑多弹丸随机分布的喷丸模型对轮盘寿命的计算方法 |
Family Cites Families (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20070005527A1 (en) * | 2005-06-06 | 2007-01-04 | Honeywell International, Inc. | Model reduction system and method for component lifing |
US8116990B2 (en) * | 2007-10-19 | 2012-02-14 | Ashok Koul | Method and system for real-time prognosis analysis and usage based residual life assessment of turbine engine components and display |
CN109507882B (zh) * | 2018-11-20 | 2021-05-07 | 中南大学 | 一种基于rbf-arx模型的快速鲁棒预测控制方法 |
CN115563818B (zh) * | 2022-12-05 | 2023-04-28 | 中国航发四川燃气涡轮研究院 | 一种考虑瞬态历程温度影响的轮盘疲劳寿命设计方法 |
-
2023
- 2023-03-02 CN CN202310188132.0A patent/CN115879247B/zh active Active
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102411641A (zh) * | 2010-09-25 | 2012-04-11 | 集美大学 | 一种基于数值模拟的汽车覆盖件回弹误差补偿方法 |
CN110697438A (zh) * | 2017-09-19 | 2020-01-17 | 中国计量大学 | 基于神经网络的直落失重式物料下料机控制器 |
CN111859729A (zh) * | 2020-06-04 | 2020-10-30 | 北京航空航天大学 | 考虑多弹丸随机分布的喷丸模型对轮盘寿命的计算方法 |
Also Published As
Publication number | Publication date |
---|---|
CN115879247A (zh) | 2023-03-31 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN115879247B (zh) | 一种基于系统辨识的轮盘关键部位应力计算方法 | |
Li et al. | A Wiener-process-model-based method for remaining useful life prediction considering unit-to-unit variability | |
EP2495631B1 (en) | A system for analysis of turbo machinery | |
Jaw | Recent advancements in aircraft engine health management (EHM) technologies and recommendations for the next step | |
US8116990B2 (en) | Method and system for real-time prognosis analysis and usage based residual life assessment of turbine engine components and display | |
JP6069498B2 (ja) | 機械部品の寿命消費を決定する方法及びシステム | |
US20060241910A1 (en) | Methods and systems for diagnosing machinery | |
US7451122B2 (en) | Empirical design of experiments using neural network models | |
US20080183444A1 (en) | Modeling and monitoring method and system | |
EP1768007A1 (en) | Monitoring a system having degrading components | |
CN109073508B (zh) | 测试涡轮单元期间的故障诊断 | |
CN108204341B (zh) | 风电场运行状态的识别方法和装置 | |
US20160123175A1 (en) | Hybrid model based detection of compressor stall | |
CH702625B1 (en) | System and method for monitoring a gas turbine. | |
JP2016145574A (ja) | ガスタービン燃焼プロファイル監視 | |
CN111581763B (zh) | 航空发动机气路故障诊断结果评价方法 | |
CN111551383A (zh) | 一种基于异质多传感器的机械状态监测方法及系统 | |
Hale et al. | Design of built-in tests for active fault detection and isolation of discrete faults | |
US8224473B2 (en) | Apparatus and method for monitoring an industrial process | |
WO2019135747A1 (en) | Probabilistic life evaluation algorithm for gas turbine engine components | |
JP2017082776A (ja) | ターボ機械を運転するリスクを判定するためのシステム及び方法 | |
US20220060494A1 (en) | Cyber monitor segmented processing for control systems | |
US20150276548A1 (en) | Condition monitoring and analytics for machines | |
Estébanez et al. | An integrated aerospace requirement setting and risk analysis tool for life cycle cost reduction and system design improvement | |
Mironov et al. | The Role of Advanced Technologies of Vibration Diagnostics to Provide Efficiency of Helicopter Life Cycle |
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 |