CN107515989A - 一种计算航空发动机叶盘高温蠕变可靠性的双重响应面法 - Google Patents
一种计算航空发动机叶盘高温蠕变可靠性的双重响应面法 Download PDFInfo
- Publication number
- CN107515989A CN107515989A CN201710767673.3A CN201710767673A CN107515989A CN 107515989 A CN107515989 A CN 107515989A CN 201710767673 A CN201710767673 A CN 201710767673A CN 107515989 A CN107515989 A CN 107515989A
- Authority
- CN
- China
- Prior art keywords
- leaf dish
- double
- creep
- analysis
- reliability
- 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.)
- Pending
Links
Classifications
-
- 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
- G06F30/23—Design optimisation, verification or simulation using finite element methods [FEM] or finite difference methods [FDM]
-
- 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
- G06F2119/00—Details relating to the type or aim of the analysis or the optimisation
- G06F2119/08—Thermal analysis or thermal optimisation
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Geometry (AREA)
- General Physics & Mathematics (AREA)
- Evolutionary Computation (AREA)
- General Engineering & Computer Science (AREA)
- Computer Hardware Design (AREA)
- Computational Mathematics (AREA)
- Mathematical Analysis (AREA)
- Mathematical Optimization (AREA)
- Pure & Applied Mathematics (AREA)
- Testing Of Devices, Machine Parts, Or Other Structures Thereof (AREA)
Abstract
本发明一种计算航空发动机叶盘高温蠕变可靠性的双重响应面法。近年来可靠性理论与方法广泛应于工程结构的可靠性分析中,在叶盘高温蠕变方面的研究甚少,且求解过程随机性强,而传统的响应面法没有考虑多种失效模式的相关性且计算速度与精度不能满足需求。本发明的组成包括:确定叶盘高温蠕变分析的特性参数,建立叶盘结构有限元模型;确定叶盘高温蠕变可靠性分析的基本参数作为随机输入变量;在高温蠕变影响下对叶盘结构进行热‑结构耦合分析;采用拉丁超立方抽样技术抽取一定数量样本点构建叶盘高温蠕变可靠性分析的双重响应面法数学模型;采用双重响应面法对叶盘进行高温蠕变可靠性的分析;对双重响应面法进行有效性验证。
Description
技术领域
本发明涉及一种计算航空发动机叶盘高温蠕变可靠性的双重响应面法,属于可靠性理论技术领域。
背景技术
航空发动机是结构极其复杂的旋转机械,其工作环境十分恶劣,承受着高温、高压和高转速的工作负荷。航空发动机在工作过程中一旦发生故障,将导致整个飞行系统无法正常运行,造成严重后果。在航空发动机中涡轮叶盘转换并传递能量,对航空发动机的整体性能起着决定性作用。为了满足航空发动机高性能要求,并确保稳健地运行,因此,研究涡轮叶盘结构的可靠性对航空发动机具有重要的意义。
近年来可靠性理论与方法广泛应于工程结构的可靠性分析中,在叶盘高温蠕变方面的研究甚少,且求解过程随机性强,而传统的响应面法没有考虑多种失效模式的相关性且计算速度与精度不能满足需求。
发明内容
本发明的目的是:在计算叶盘高温蠕变可靠性时,考虑热-结构耦合场之间的相互影响以及材料属性和载荷的非线性,针对传统方法的不足,提出一种计算航空发动机叶盘高温蠕变可靠性的方法,保证精度提高计算速度与效率。
本发明提供了一种计算航空发动机叶盘高温蠕变可靠性的双重响应面法,其具体过程如下:
a.确定叶盘高温蠕变分析的特性参数,建立叶盘结构有限元模型;
b.确定叶盘高温蠕变可靠性分析的基本参数作为随机输入变量;
c.在高温蠕变影响下对叶盘结构进行热-结构耦合分析;
d.采用拉丁超立方抽样技术抽取一定数量样本点构建叶盘高温蠕变可靠性分析的双重响应面法数学模型;
e.采用双重响应面法对叶盘进行高温蠕变可靠性的分析;
f.对双重响应面法进行有效性验证。
所述的一种计算航空发动机叶盘高温蠕变可靠性的双重响应面法,所述的步骤a中,根据蠕变实验数据确定蠕变特性参数,在有限元软件中对叶盘结构进行四面体网格划分。
所述的一种计算航空发动机叶盘高温蠕变可靠性的双重响应面法,所述的步骤b中,根据发动机叶盘实际工况,确定影响叶盘高温蠕变可靠性的外部载荷因素和本身结构性能参数。
所述的一种计算航空发动机叶盘高温蠕变可靠性的双重响应面法,所述的步骤c中,热分析是通过热传导将高温燃气的热量传递叶盘结构上,利用能量守恒方程对其进行热计算;结构分析是通过插值方法将热分析结果信息传递到叶盘结构,采用有限元法进行运算,而分析结果则是运用有限元的基本方程计算获得。
所述的一种计算航空发动机叶盘高温蠕变可靠性的双重响应面法,所述的步骤d中,考虑材料属性的非线性和载荷的动态性,在叶盘总体应力和应变出现位置利用拉丁超立方抽样技术对随机输入变量进行小批量的抽样,根据抽样点计算出相应的输出响应,得到拟合方程所需要的数据点,基于最小二乘法通过MATLAB矩阵实验室工具箱对数据进行回归分析建立双重响应面数学模型。
所述的一种计算航空发动机叶盘高温蠕变可靠性的双重响应面法,所述的步骤e中,采用双重响应面法对叶盘进行高温蠕变可靠性分析,即采用MATLAB编写可靠性计算的程序并对双重响应面法数学模型进行高温蠕变可靠性的计算。
所述的一种计算航空发动机叶盘高温蠕变可靠性的双重响应面法,所述的步骤f中,对双重响应面法进行有效性性验证,即在相同的计算条件下,以公认的可靠性计算方法蒙特卡罗法为基准,和响应面法相对比。
本发明与现有技术相比具有以下有益效果:
1. 本发明考虑了高温蠕变对应力、应变的影响以及输入变量的随机性、多种失效模式共存状况下的相关性。本发明能在满足精度条件下,显著的提高计算速度与效率,减少计算成本。
2.本发明能有效的解决背景技术中所述的问题,该方法使用MATLAB中已有的分类程序通过数据接口文件实现对有限元分析软件ANSYS的直接调用,结合MATLAB便捷的编程功能和ANSYS分析功能,方便工程师使用。
附图说明
图1为叶盘高温蠕变双重响应面法的可靠性分析流程图。
具体实施方式
实施例1
一种计算航空发动机叶盘高温蠕变可靠性的双重响应面法,包括以下步骤:
a.确定叶盘高温蠕变分析的特性参数,建立叶盘结构有限元模型;
b.确定叶盘高温蠕变可靠性分析的基本参数最为随机输入变量;
c.在高温蠕变影响下对叶盘结构进行热-结构耦合分析;
d.采用拉丁超立方抽样技术抽取一定数量样本点构建叶盘高温蠕变可靠性分析的双重响应面法数学模型;
e.采用双重响应面法对叶盘进行高温蠕变可靠性的分析;
f.对双重响应面法进行有效性验证。
实施例2
根据实施例1所述的一种计算航空发动机叶盘高温蠕变可靠性的双重响应面法,所述的步骤a中,根据蠕变实验数据确定蠕变特性参数,在有限元软件中对叶盘结构进行四面体网格划分。
实施例3
根据实施例1所述的一种计算航空发动机叶盘高温蠕变可靠性的双重响应面法,所述的步骤b中,根据发动机叶盘实际工况,确定影响叶盘高温蠕变可靠性的外部载荷因素和本身结构性能参数。
实施例4
根据实施例1所述的一种计算航空发动机叶盘高温蠕变可靠性的双重响应面法,所述的步骤c中,热分析是通过热传导将高温燃气的热量传递叶盘结构上,利用能量守恒方程对其进行热计算;结构分析是通过插值方法将热分析结果信息传递到叶盘结构,采用有限元法进行运算,而分析结果则是运用有限元的基本方程计算获得。
实施例5
根据实施例1所述的一种计算航空发动机叶盘高温蠕变可靠性的双重响应面法,所述的步骤d中,考虑材料属性的非线性和载荷的动态性,在叶盘总体应力和应变出现位置利用拉丁超立方抽样技术对随机输入变量进行小批量的抽样,根据抽样点计算出相应的输出响应,得到拟合方程所需要的数据点,基于最小二乘法通过MATLAB矩阵实验室工具箱对数据进行回归分析建立双重响应面数学模型。
实施例6
根据实施例1所述的一种计算航空发动机叶盘高温蠕变可靠性的双重响应面法,所述的步骤e中,采用双重响应面法对叶盘进行高温蠕变可靠性分析,即采用MATLAB编写可靠性计算的程序并对双重响应面法数学模型进行高温蠕变可靠性的计算。
实施例7
根据实施例1所述的一种计算航空发动机叶盘高温蠕变可靠性的双重响应面法,所述的步骤f中,对双重响应面法进行有效性性验证,即在相同的计算条件下,以公认的可靠性计算方法蒙特卡罗法为基准,和响应面法相对比。
Claims (7)
1.一种计算航空发动机叶盘高温蠕变可靠性的双重响应面法,其特征在于,包括以下步骤:
确定叶盘高温蠕变分析的特性参数,建立叶盘结构有限元模型;
选取叶盘高温蠕变可靠性分析的基本参数作为随机输入变量;
在高温蠕变影响下对叶盘结构进行热-结构耦合分析;
采用拉丁超立方抽样技术抽取一定数量样本点构建叶盘高温蠕变可靠性分析的双重响应面法数学模型;
采用双重响应面法对叶盘进行高温蠕变可靠性的分析;
对双重响应面法进行有效性验证。
2.根据权利要求1所述的一种计算航空发动机叶盘高温蠕变可靠性的双重响应面法,所述的步骤a中,根据蠕变实验数据确定蠕变特性参数,在有限元软件中对叶盘结构进行四面体网格划分。
3.根据权利要求1所述的一种计算航空发动机叶盘高温蠕变可靠性的双重响应面法,所述的步骤b中,根据发动机叶盘实际工况,确定影响叶盘高温蠕变可靠性的外部载荷因素和本身结构性能参数。
4.根据权利要求1所述的一种计算航空发动机叶盘高温蠕变可靠性的双重响应面法,所述的步骤c中,热分析是通过热传导将高温燃气的热量传递叶盘结构上,利用能量守恒方程对其进行热计算;结构分析是通过插值方法将热分析结果信息传递到叶盘结构,采用有限元法进行运算,而分析结果则是运用有限元的基本方程计算获得。
5.根据权利要求1所述的一种计算航空发动机叶盘高温蠕变可靠性的双重响应面法,所述的步骤d中,考虑材料属性的非线性和载荷的动态性,在叶盘总体应力和应变出现位置利用拉丁超立方抽样技术对随机输入变量进行小批量的抽样,根据抽样点计算出相应的输出响应,得到拟合方程所需要的数据点,基于最小二乘法通过MATLAB矩阵实验室工具箱对数据进行回归分析建立双重响应面数学模型。
6.根据权利要求1所述的一种计算航空发动机叶盘高温蠕变可靠性的双重响应面法,所述的步骤e中,采用双重响应面法对叶盘进行高温蠕变可靠性分析,即采用MATLAB编写可靠性计算的程序并对双重响应面法数学模型进行高温蠕变可靠性的计算。
7.根据权利要求1所述的一种计算航空发动机叶盘高温蠕变可靠性的双重响应面法,所述的步骤f中,对双重响应面法进行有效性性验证,即在相同的计算条件下,以公认的可靠性计算方法蒙特卡罗法为基准,和响应面法相对比。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710767673.3A CN107515989A (zh) | 2017-08-31 | 2017-08-31 | 一种计算航空发动机叶盘高温蠕变可靠性的双重响应面法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710767673.3A CN107515989A (zh) | 2017-08-31 | 2017-08-31 | 一种计算航空发动机叶盘高温蠕变可靠性的双重响应面法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN107515989A true CN107515989A (zh) | 2017-12-26 |
Family
ID=60724861
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201710767673.3A Pending CN107515989A (zh) | 2017-08-31 | 2017-08-31 | 一种计算航空发动机叶盘高温蠕变可靠性的双重响应面法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN107515989A (zh) |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101104898A (zh) * | 2007-06-19 | 2008-01-16 | 中国科学院金属研究所 | 一种高热强性、高热稳定性的高温钛合金 |
CN105404756A (zh) * | 2015-12-17 | 2016-03-16 | 北京航空航天大学 | 一种涡轮叶盘结构寿命可靠性设计方法 |
CN106485032A (zh) * | 2016-11-08 | 2017-03-08 | 哈尔滨理工大学 | 一种计算叶盘振动可靠性的双重极值响应面法 |
-
2017
- 2017-08-31 CN CN201710767673.3A patent/CN107515989A/zh active Pending
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101104898A (zh) * | 2007-06-19 | 2008-01-16 | 中国科学院金属研究所 | 一种高热强性、高热稳定性的高温钛合金 |
CN105404756A (zh) * | 2015-12-17 | 2016-03-16 | 北京航空航天大学 | 一种涡轮叶盘结构寿命可靠性设计方法 |
CN106485032A (zh) * | 2016-11-08 | 2017-03-08 | 哈尔滨理工大学 | 一种计算叶盘振动可靠性的双重极值响应面法 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Bai et al. | Distributed collaborative response surface method for mechanical dynamic assembly reliability design | |
Fei et al. | Novel method and model for dynamic reliability optimal design of turbine blade deformation | |
CN102855347A (zh) | 一种排气歧管热机械强度的有限元分析方法 | |
Orozco et al. | A new methodology of thermodynamic diagnosis, using the thermoeconomic method together with an artificial neural network (ANN): A case study of an externally fired gas turbine (EFGT) | |
CN106485032A (zh) | 一种计算叶盘振动可靠性的双重极值响应面法 | |
JP2013510365A (ja) | 発電プラントの熱効率を最大化するシステム及び方法 | |
CN107729658A (zh) | 一种可靠性优化设计的模糊智能分解协调极值响应面法 | |
CN108038320A (zh) | 一种计算失谐叶盘振动可靠性的多重响应面法 | |
Zhang et al. | Multiobject reliability analysis of turbine blisk with multidiscipline under multiphysical field interaction | |
RU2621422C2 (ru) | Система и способ тестирования показателя работы паровой турбины | |
Brinkert et al. | Understanding the twin scroll turbine: flow similarity | |
CN107515989A (zh) | 一种计算航空发动机叶盘高温蠕变可靠性的双重响应面法 | |
Montero Carrero et al. | Experimental characterisation of a humidified T100 micro gas turbine | |
Budnik et al. | Exergetic cost of steam power plant operation | |
Eroglu et al. | Thermal analysis of heavy duty engine exhaust manifold using cfd | |
Bai et al. | Dynamic probabilistic analysis of stress and deformation for bladed disk assemblies of aeroengine | |
CN109388884A (zh) | 一种计算耦合叶盘疲劳寿命的广义回归极值响应面法 | |
CN108549348B (zh) | 锅炉物理场运行信息获取方法、装置和计算机设备 | |
Zhai et al. | Reliability and sensitivity analyses of HPT blade-tip radial running clearance using multiply response surface model | |
CN111523274B (zh) | 汽轮机高温承压部件数字孪生体构建方法及监控系统 | |
Ram et al. | Gas turbine power plant performance evaluation under key failures | |
Gorla et al. | Probabilistic structural and thermal analysis of a gasketed flange | |
Li et al. | Complex physical-model based dynamic system safety analysis of Aviation Piston Engine considering hybrid uncertainty of fault | |
Wang et al. | Effect of the secondary microcontact state induced by the all-factor thermal effect on bidirectional thermal contact performance | |
Ramerth et al. | A Probabilistic Secondary Flow System Design Process for Gas Turbine Engines |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
WD01 | Invention patent application deemed withdrawn after publication |
Application publication date: 20171226 |
|
WD01 | Invention patent application deemed withdrawn after publication |