CN114912191A - 一种导叶可调的涡扇发动机部件级机载动态模型设计方法 - Google Patents
一种导叶可调的涡扇发动机部件级机载动态模型设计方法 Download PDFInfo
- Publication number
- CN114912191A CN114912191A CN202210424287.5A CN202210424287A CN114912191A CN 114912191 A CN114912191 A CN 114912191A CN 202210424287 A CN202210424287 A CN 202210424287A CN 114912191 A CN114912191 A CN 114912191A
- Authority
- CN
- China
- Prior art keywords
- model
- lpv
- pressure ratio
- designing
- turbofan engine
- 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
- 238000000034 method Methods 0.000 title claims abstract description 28
- 238000004364 calculation method Methods 0.000 claims abstract description 25
- 239000000446 fuel Substances 0.000 claims abstract description 12
- 238000012937 correction Methods 0.000 claims abstract description 8
- 239000011159 matrix material Substances 0.000 claims description 24
- 238000013461 design Methods 0.000 claims description 10
- 239000000295 fuel oil Substances 0.000 claims description 8
- 238000005457 optimization Methods 0.000 claims description 6
- 230000000694 effects Effects 0.000 claims description 5
- 238000004422 calculation algorithm Methods 0.000 claims description 4
- 238000012546 transfer Methods 0.000 claims description 3
- 230000003416 augmentation Effects 0.000 abstract 1
- 238000004088 simulation Methods 0.000 description 5
- 238000000638 solvent extraction Methods 0.000 description 5
- PCTMTFRHKVHKIS-BMFZQQSSSA-N (1s,3r,4e,6e,8e,10e,12e,14e,16e,18s,19r,20r,21s,25r,27r,30r,31r,33s,35r,37s,38r)-3-[(2r,3s,4s,5s,6r)-4-amino-3,5-dihydroxy-6-methyloxan-2-yl]oxy-19,25,27,30,31,33,35,37-octahydroxy-18,20,21-trimethyl-23-oxo-22,39-dioxabicyclo[33.3.1]nonatriaconta-4,6,8,10 Chemical compound C1C=C2C[C@@H](OS(O)(=O)=O)CC[C@]2(C)[C@@H]2[C@@H]1[C@@H]1CC[C@H]([C@H](C)CCCC(C)C)[C@@]1(C)CC2.O[C@H]1[C@@H](N)[C@H](O)[C@@H](C)O[C@H]1O[C@H]1/C=C/C=C/C=C/C=C/C=C/C=C/C=C/[C@H](C)[C@@H](O)[C@@H](C)[C@H](C)OC(=O)C[C@H](O)C[C@H](O)CC[C@@H](O)[C@H](O)C[C@H](O)C[C@](O)(C[C@H](O)[C@H]2C(O)=O)O[C@H]2C1 PCTMTFRHKVHKIS-BMFZQQSSSA-N 0.000 description 3
- 238000013178 mathematical model Methods 0.000 description 3
- 238000006243 chemical reaction Methods 0.000 description 2
- 238000002485 combustion reaction Methods 0.000 description 2
- 238000010586 diagram Methods 0.000 description 2
- 238000005516 engineering process Methods 0.000 description 2
- 239000003921 oil Substances 0.000 description 2
- 238000011160 research Methods 0.000 description 2
- 150000001875 compounds Chemical class 0.000 description 1
- 238000010276 construction Methods 0.000 description 1
- 238000013500 data storage Methods 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
- 238000012360 testing method Methods 0.000 description 1
- 238000012795 verification 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/15—Vehicle, aircraft or watercraft design
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
- G06F17/16—Matrix or vector computation, e.g. matrix-matrix or matrix-vector multiplication, matrix factorization
-
- 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/02—Reliability analysis or reliability optimisation; Failure analysis, e.g. worst case scenario performance, failure mode and effects analysis [FMEA]
-
- 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)
- General Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Geometry (AREA)
- Mathematical Optimization (AREA)
- General Engineering & Computer Science (AREA)
- Computational Mathematics (AREA)
- Mathematical Analysis (AREA)
- Pure & Applied Mathematics (AREA)
- Evolutionary Computation (AREA)
- Computer Hardware Design (AREA)
- Mathematical Physics (AREA)
- Data Mining & Analysis (AREA)
- Algebra (AREA)
- Computing Systems (AREA)
- Databases & Information Systems (AREA)
- Software Systems (AREA)
- Automation & Control Theory (AREA)
- Aviation & Aerospace Engineering (AREA)
- Feedback Control In General (AREA)
- Combined Controls Of Internal Combustion Engines (AREA)
Abstract
本发明公开一种导叶可调的涡扇发动机部件级机载动态模型设计方法，结合航空发动机气动热力学物理建模与LPV建模思想，建立含加力的发动机的机载动态实时模型。该方法包括：在原有部件级模型的基础上，引入压气机导叶特性修正；根据发动机开加力时的部件级模型的数据，将加力燃油量作为控制量之一，建立加力接通后关于压比状态参数的含加力宽工况范围LPV模型；设计引入并优化状态参数惯性环节计算策略，获得包线内不同飞行状态平稳切换的宽工况范围机载动态模型。本发明将导叶特性修正、建立机载加力燃烧室模型、状态参数惯性环节计算策略用于建立含加力的涡扇发动机机载动态实时模型，可获得在包线内慢车以上全状态的涡扇发动机机载实时计算模型。
Description
技术领域
本发明涉及航空发动机建模与仿真领域，尤其涉及导叶可调的含加力涡扇发动机机载动态模型设计。
背景技术
现代航空发动机是一个极其复杂的非线性气动热力学系统，在研发过程中，需要大量的人员和资金投入。航空发动机数学模型是控制系统设计研究的基础，其中包括三种航空发动机数学模型：部件级模型、状态变量模型和LPV模型。
部件级模型：精度较高且能够全面反映发动机各截面工作状态。然而，由于其计算时需要反复的流量迭代计算，因此实时性较差。在控制算法实时仿真验证时，既要求数学模型有较高的精度，又要求其有较好的实时性。
状态变量模型：利用非线性部件级模型局部线性建立的模型，目前广泛应用于发动机多变控制器设计。建立航空发动机状态变量模型的方法有偏导数法、拟合法和最小二乘法。
LPV模型：由于状态变量模型是小范围的线性模型，难以满足发动机大范围的精度要求，因此需要应用LPV建模方法。LPV模型通过可测量的外部实时参数作为调度参数对系统的动态特征进行描述，LPV控制器直接采用线性鲁棒控制理论进行设计，保证所设计的控制器具有鲁棒稳定性。LPV系统在工程应用与理论研究上均具有很高的潜在价值，因此，近年来LPV建模技术受到了控制学界的高度重视。
加力燃烧一般在发动机低压涡轮后、喷管前的加力燃烧室内完成。加力燃烧室可以在特定的情况下，大幅提升发动机的推力，从而使飞行器获得短暂优势，作用不容忽视。本发明将涡扇发动机各个部件模型以及建立的关于转速和压比的LPV模型，提出一种导叶可调的含加力的涡扇发动机机载实时模型的构建方法，在精度损失较小的情况下提高发动机模型的实时性。
发明内容
本发明所要解决的技术问题是针对背景技术的缺陷，提供一种具有高精度的导叶可调的机载动态实时模型，解决原有的非线性部件级模型实时性不足，线性化模型误差较大的问题。
本发明为解决上述技术问题采用以下技术方案：
步骤A)，在原有部件级模型的基础上，引入压气机导叶特性修正；设计LPV形式的非线性共同工作方程组的无迭代求解算法，根据发动机开加力时的部件级模型的数据，将加力燃油量作为控制量之一，建立加力接通后关于压比状态参数的含加力宽工况范围LPV模型；
步骤B)，设计引入并优化状态参数惯性环节计算策略，由此获得包线内不同飞行状态平稳切换的宽工况范围机载动态模型。
作为本发明一种导叶可调的涡扇发动机部件级机载动态模型设计方法的进一步优化方案，步骤A)的具体步骤如下：
步骤A1)，在原有部件级模型的基础上，将压气机前三级与后六级分开建模，引入压气机导叶特性修正。
步骤A2)，根据发动机开加力时的部件级模型的数据，将加力燃油量作为控制量之一，建立加力接通后关于压比状态参数的含加力状态变量模型。
步骤A3)，在不同的进口温度，分别在不同的喉道面积A8与不同的导叶角角度Afa下，将在此喉道面积和导叶角角度时的不同转速下的状态变量模型组合成LPV模型。
作为一种含加力的涡扇发动机机载动态实时模型的构建方法，提出了一种减小误差的优化方案步骤B)的具体步骤如下：
步骤B1)，建立基于等温线的飞行包线划分方案。根据飞行包线划分方案，优化选取若干个T2下的包线点建立加力接通后关于压比状态参数的含加力宽工况范围LPV模型；
步骤B2)，针对全包线内飞行点，设计引入并优化状态参数惯性环节计算策略。由LPV模型求解获得换算转速与压比，提供给气路各部件气动热力学模型完成计算，由此获得状态平稳切换的宽工况范围机载动态模型。
进一步的，所述步骤B1)中飞行包线划分方案为：在包线范围内，以进口温度T2每隔20K建立一条等温线，由此获得包线内含加力宽工况范围LPV模型。
进一步的，所述步骤B2)中具体步骤如下：
步骤B2.1)，根据涡扇发动机所处的进口总温T2，计算T2与已建立LPV模型的T2的差值，当T2与已建立LPV模型的T2相差均不超过20K时，可选择差值最小的LPV模型。选择相应等温线下多项式系数，计算系数矩阵中的元素，进而得到当前高压转速nH下的各系数矩阵A，B，C，D，根据当前A8进行插值计算当前喉道截面面积A8下的转速压比。
步骤B2.2)，在原有模型的基础上，设计引入并优化状态参数惯性环节计算策略。状态参数惯性环节计算策略如图1(b)所示。
步骤B2.3)，将求出的压比与转速再送回非线性模型中，提供给气路各部件气动热力学模型完成计算，由此获得状态平稳切换的宽工况范围机载动态模型。
进一步的，步骤B2.2)中优化策略的具体步骤如下：在切换系数矩阵过程中，判断输入参数的改变是否需要更改矩阵系数，如不需要则继续使用上一时刻的矩阵系数，如需要更换矩阵系数，则重复计算若干次以达到稳态效果；在得到转速与压比后，引入转速与压比的状态参数惯性环节，减小模型在切换矩阵过程中的输出误差，一阶惯性环节的传递函数表示如下：
其中T为时间常数，s表示复频域。
本发明采用以上技术方案与现有方案相比，具有以下技术效果：
(1)本发明所使用的部件级模型，是在原有部件级模型的基础上，将压气机前六级与后三级分开建模，引入压气机导叶特性修正；
(2)本发明所使用的机载动态实时模型，将加力燃油量作为控制量之一，建立了机载加力燃烧室模型；
(3)本发明在LPV模型求解后，设计引入并优化状态参数惯性环节计算策略，能有效减小系数矩阵更换导致的误差。
附图说明
图1中，(a)是含加力的涡扇发动机机载动态实时模型原理简图，(b)是对优化计算策略的详细说明；
图2是基于等温线的包线划分计划；
图3是H＝0km,Ma＝0,Tt2＝288.15K时归一化后的燃油流量Wf、加力燃油流量Wfa和尾喷管喉道界面面积A8变化；
图4是H＝0km,Ma＝0,Tt2＝288.15K时输出参数对比；
图5是H＝11km,Ma＝0.8,Tt2＝244.38K时归一化后的燃油流量Wf、加力燃油流量Wfa和尾喷管喉道界面面积A8变化；
图6是H＝11km,Ma＝0.8,Tt2＝244.38K时输出参数对比；
图7是H＝11km,Ma＝1.5,Tt2＝314.3K时归一化后的燃油流量Wf、加力燃油流量Wfa和尾喷管喉道界面面积A8变化；
图8是H＝11km,Ma＝1.5,Tt2＝314.3K时输出参数对比。
具体实施方式
下面结合附图对本发明的技术方案做进一步的详细说明。
本发明的思路是在原有部件级模型的基础上，引入压气机导叶特性修正；根据发动机开加力时的部件级模型的数据，将加力燃油量作为控制量之一，建立加力接通后关于压比状态参数的宽工况范围LPV模型；根据包线内飞行点确定相应LPV模型，设计引入并优化状态参数惯性环节计算策略，由此获得状态平稳切换的宽工况范围机载动态模型。
本发明的具体实施方式以某型含加力的涡扇发动机部件级机载动态实时模型设计为例，图1(a)是涡扇发动机机载动态模型设计原理图，本发明说明的一种导叶可调的涡扇发动机部件级机载动态模型设计方法，具体包括以下步骤：
步骤A)，在原有部件级模型的基础上，引入压气机导叶特性修正；设计LPV形式的非线性共同工作方程组的无迭代求解算法，根据发动机开加力时的部件级模型的数据，将加力燃油量作为控制量之一，建立加力接通后关于压比状态参数的含加力宽工况范围LPV模型；
步骤A1)，在原有部件级模型的基础上，将压气机前三级与后六级分开建模，引入压气机导叶特性修正。
当导叶角小于10°时，则可使用导叶角进行插值计算：
W＝W1-(W1-W2)·Afa/10
eff＝eff1-(eff1-eff2)·Afa/10
W、eff为此导叶角下流量与效率，W1、eff1为利用第一个特性图插值得出的流量与效率，W2、eff2为利用第二个特性图插值得出的流量与效率。
步骤A2)，根据发动机开加力时的部件级模型的数据，将加力燃油量作为控制量之一，建立加力接通后关于压比状态参数的含加力状态变量模型。
发动机在稳态点(x0,u0,y0)的状态变量模型为：
Δy＝CΔx+DΔu
式中，状态量x为高低压转子转速NH和NL，控制量u为加力燃烧室供油量Wfa和尾喷管喉道面积A8，输出量y为风扇压比πfan、低压压气机压比πlcomp、高压压气机压比πhcomp、高压涡轮压比πhturb和低压涡轮压比πlturb，Δx＝x-x0，Δy＝y-y0。
参数相似归一化如下：
式中，PNL、PNH、PWfa、PA8分别为归一化后的高压转子转速、低压转子转速、加力燃烧室供油量、尾喷管喉道面积，T2为发动机进口总温，下标ds表示发动机设计点参数，用相似归一化参数表示的发动机状态变量模型如下：
式中，
分别表示矩阵A、B、C、D，利用小扰动法获得此工作点的系数矩阵A、B、C、D的初始解。
与直接计算不含加力时的状态变量模型相比，其中一个控制量均为A8,但在计算含加力燃烧室的状态变量模型时，主燃烧室燃油量Wf此时需要达到最大值，因此主燃烧室燃油量Wf默认为最大值。在使用加力燃烧室时，不采用Wf为控制量，另一个控制量选取加力燃烧室供油量Wfa；
步骤A3)，在不同的进口温度，分别在不同的喉道面积A8与不同的导叶角角度Afa下，将大量的状态变量模型组合成LPV模型：
发动机不同的喉道面积A8、导叶角afa和不同高压转速下的状态变量模型构成转速与压比的LPV模型，再利用相似换算理论将该模型在全包线内进行拓展，将发动机的LPV模型描述为：
式中，x为状态量，y为输出量，下标cor表示相似换算到地面点参数。
为了减少数据存储的空间，对数量为k的ABCD矩阵中各元素进行3阶多项式拟合：
通过拟合法多次试取获得
从而建立宽工况范围开加力时的关于压比的LPV动态实时模型。
步骤B)，设计引入并优化状态参数惯性环节计算策略，由此获得包线内不同飞行状态平稳切换的宽工况范围机载动态模型。
步骤B1)，建立基于等温线的飞行包线划分方案，飞行包线划分方案如图(2)所示。根据飞行包线划分方案，优化选取若干个T2下的包线点建立加力接通后关于压比状态参数的含加力宽工况范围LPV模型；
步骤B1)中飞行包线划分方案为：在包线范围内，以进口温度T2每隔20K建立一条等温线，飞行包线划分方案如图2所示，由此获得包线内含加力宽工况范围LPV模型。
步骤B2)，针对全包线内飞行点，设计引入并优化状态参数惯性环节计算策略。由LPV模型求解获得换算转速与压比，提供给气路各部件气动热力学模型完成计算，由此获得状态平稳切换的宽工况范围机载动态模型。
步骤B2.1)，利用含加力的关于压比状态参数LPV动态实时模型记录对应不同工况下非线性模型迭代收敛所解得的五个压比。
根据涡扇发动机所处的进口总温T2，当T2与已建立LPV模型的T2相差均不超过20K时，可使用相邻温度差距更小的LPV模型。选择相应等温线下多项式系数，计算系数矩阵中的元素，进而得到当前高压转速nH下的各系数矩阵A，B，C，D，根据当前A8进行插值计算当前喉道截面面积A8下的转速压比。若当前喉道面积为A8，当A8i≤A8≤A8i+1时，根据A8i+1和A8i分别计算在当前Afa和Wfa下两组系数矩阵中的元素，进而得到当前高压转速nH下的各系数矩阵A，B，C，D，通过LPV模型求解转速与各部件压比的计算如下；
式中，n＝[nL,nH]T，π＝[πFan,πLcomp,πHcomp,πHTurb,πLTurb]T为五个旋转部件的压比，u＝[Wfa,A8]T，下标0为稳态点。分别根据A8i和A8i+1计算得到在此Wfa下的两组转速压比。
根据当前A8进行插值计算当前喉道截面面积A8下的转速压比。
步骤B2.2)，由于修正前部件级动态实时模型输出参数有抖动误差，因此在原有模型的基础上，设计引入并优化状态参数惯性环节计算策略。
优化策略如图1(b)所示。优化策略具体为：在切换系数矩阵过程中，判断输入参数的改变是否需要更改矩阵系数，如不需要则继续使用上一时刻的矩阵系数，如需要更换矩阵系数，则重复计算100次以达到稳态效果。在得到转速与压比后，引入转速与压比的状态参数惯性环节，减小模型在切换矩阵过程中的输出误差，一阶惯性环节的传递函数表示如下：
其中T为时间常数；
步骤B2.3)，将求出的压比与转速再送回非线性模型中，与气路各部件气动热力学模型组成部件级机载动态模型，使得含加力的非线性模型可以在不迭代的情况下在不同工作点获取发动机稳态下的关键输出参数P25,T25,Ps3,P6,T6，若误差不超过2％，则可获得状态平稳切换的宽工况范围机载动态模型。
为了验证本发明所设计的导叶可调的涡扇发动机机载动态模型设计的有效性，在MATLAB环境下进行了包线内导叶可调的机载动态模型输出的数字仿真。
首先，以进口温度Tt2分别等于[244.38K、259K、273K、288.15K、314.3K、352.9K]作为覆盖全包线内的六条等温线，包线划分方案如图2所示。
在这六条等温线下分别计算：A8＝[0.191，0.211，0.218，0.266，0.314，0.362，0.41]，afa＝[0,0.156,0.279,0.351,0.369,0.4,1.4,2.4,3.4,4.4,5.4,6.4,7.4,8,9,10]下不同高压转速下含加力时的涡扇发动机数据建立关于转速和压比的状态变量模型即系数矩阵ABCD，并对不同高压转速下系数矩阵的对应元素进行3次多项式拟合，获取开加力模式下不同A8、afa和不同高压转速下ABCD矩阵元素多项式拟合系数。
在(H＝0km，Ma＝0)、(H＝11km，Ma＝0.8)、(H＝11km，Ma＝1.5)三个工作点下进行仿真试验，燃油流量Wf、加力燃油流量Wfa和尾喷管喉道界面面积A8变化如图3、图5、图7所示，导叶角Afa均为0保持不变，归一化后的机载动态实时模型与动态模型的输出对比如图4、图6、图8所示，最大稳态误差和最大动态误差如表1、2所示，在表1和表2中，最大稳态误差均未超过1％，最大稳态误差均未超过5％，结果精度较高。耗时对比如表3所示，在表3中，机载动态实时模型与部件级模型相比耗时均减少80％以上。综合以上仿真结果，该方法实现了在精度损失较小的情况下获得较高实时性能的目的。
表1最大稳态误差(％)
表2最大动态误差(％)
表3耗时对比
需要指出的是，以上所述，仅为本发明的具体实施方式，但本发明的保护范围并不局限于此，任何熟悉本技术领域的技术人员在本发明揭露的技术范围内，可轻易想到的变化和替换，都应涵盖在本发明的保护范围内。因此，本发明的保护范围应以所述权利要求的保护范围为准。
Claims (6)
1.一种导叶可调的涡扇发动机部件级机载动态模型设计方法，其特征在于，包括以下步骤：
步骤A)，在原有部件级模型的基础上，引入压气机导叶特性修正；设计LPV形式的非线性共同工作方程组的无迭代求解算法，根据发动机开加力时的部件级模型的数据，将加力燃油量作为控制量之一，建立加力接通后关于压比状态参数的含加力宽工况范围LPV模型；
步骤B)，设计引入并优化状态参数惯性环节计算策略，由此获得包线内不同飞行状态平稳切换的宽工况范围机载动态模型。
2.如权利要求1所示的一种导叶可调的涡扇发动机部件级机载动态模型设计方法，其特征在于，所述步骤A)的具体步骤如下：
步骤A1)，在原有部件级模型的基础上，将压气机前三级与后六级分开建模，引入压气机导叶特性修正；
步骤A2)，根据发动机开加力时的部件级模型的数据，将加力燃油量作为控制量之一，建立加力接通后关于压比状态参数的含加力状态变量模型；
步骤A3)，在不同的进口温度，分别在不同的喉道面积A8与不同的导叶角角度Afa下，将在此喉道面积和导叶角角度时的不同转速下的状态变量模型组合成LPV模型。
3.如权利要求1所示的一种导叶可调的涡扇发动机部件级机载动态模型设计方法，其特征在于，所述步骤B)的具体步骤如下：
步骤B1)，建立基于等温线的飞行包线划分方案，根据飞行包线划分方案，优化选取若干个进口总温T2下的包线点建立加力接通后关于压比状态参数的含加力宽工况范围LPV模型；
步骤B2)，针对全包线内飞行点，设计引入并优化状态参数惯性环节计算策略；由LPV模型求解获得换算转速与压比，提供给气路各部件气动热力学模型完成计算，由此获得状态平稳切换的宽工况范围机载动态模型。
4.如权利要求3所示的一种导叶可调的涡扇发动机部件级机载动态模型设计方法，其特征在于，所述步骤B1)中飞行包线划分方案为：在包线范围内，以进口温度T2每隔20K建立一条等温线，依据等温线建立包线内含加力宽工况范围LPV模型。
5.如权利要求3所示的一种导叶可调的涡扇发动机部件级机载动态模型设计方法，其特征在于，所述步骤B2)的具体步骤如下：
步骤B2.1)，根据涡扇发动机所处的进口总温T2，计算T2与已建立LPV模型的T2的差值，当T2与已建立LPV模型的T2相差均不超过20K时，可选择差值最小的LPV模型；选择相应等温线下多项式系数，计算系数矩阵中的元素，进而得到当前高压转速nH下的各系数矩阵A，B，C，D，根据当前A8进行插值计算当前喉道截面面积A8下的转速压比；
步骤B2.2)，在原有模型的基础上，设计引入并优化状态参数惯性环节计算策略；
步骤B2.3)，将求出的压比与转速再送回非线性模型中，提供给气路各部件气动热力学模型完成计算，由此获得状态平稳切换的宽工况范围机载动态模型。
6.如权利要求5所示的一种导叶可调的涡扇发动机部件级机载动态模型设计方法，其特征在于，步骤B2.2)中优化策略的具体步骤如下：在切换系数矩阵过程中，判断输入参数的改变是否需要更改矩阵系数，如不需要则继续使用上一时刻的矩阵系数，如需要更换矩阵系数，则重复计算若干次以达到稳态效果；在得到转速与压比后，引入转速与压比的状态参数惯性环节，减小模型在切换矩阵过程中的输出误差，一阶惯性环节的传递函数表示如下：
其中T为时间常数，s表示复频域。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202210424287.5A CN114912191A (zh) | 2022-04-21 | 2022-04-21 | 一种导叶可调的涡扇发动机部件级机载动态模型设计方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202210424287.5A CN114912191A (zh) | 2022-04-21 | 2022-04-21 | 一种导叶可调的涡扇发动机部件级机载动态模型设计方法 |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN114912191A true CN114912191A (zh) | 2022-08-16 |
Family
ID=82764634
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202210424287.5A Pending CN114912191A (zh) | 2022-04-21 | 2022-04-21 | 一种导叶可调的涡扇发动机部件级机载动态模型设计方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN114912191A (zh) |
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN115680903A (zh) * | 2022-08-31 | 2023-02-03 | 中国航发四川燃气涡轮研究院 | 可回收旁路引气控制方法 |
Citations (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN107122537A (zh) * | 2017-04-21 | 2017-09-01 | 中国航发沈阳发动机研究所 | 一种加力燃烧室总余气系数设计方法 |
| CN110362960A (zh) * | 2019-08-09 | 2019-10-22 | 哈尔滨工业大学 | 基于多胞折合平衡流形展开模型的航空发动机系统辨识方法 |
| CN111475889A (zh) * | 2020-03-30 | 2020-07-31 | 西安交通大学 | 一种物理与数据驱动的航空发动机数字孪生建模方法 |
| CN112733266A (zh) * | 2020-12-31 | 2021-04-30 | 华南理工大学 | 一种动力总成悬置系统模糊可靠性优化方法 |
| CN114297786A (zh) * | 2021-11-08 | 2022-04-08 | 南京航空航天大学 | 一种涡扇发动机机载动态模型全包线精度扩展的方法 |
| CN115017625A (zh) * | 2022-06-28 | 2022-09-06 | 南京航空航天大学 | 基于基准特性线的变循环发动机变几何部件特性建模方法 |
-
2022
- 2022-04-21 CN CN202210424287.5A patent/CN114912191A/zh active Pending
Patent Citations (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN107122537A (zh) * | 2017-04-21 | 2017-09-01 | 中国航发沈阳发动机研究所 | 一种加力燃烧室总余气系数设计方法 |
| CN110362960A (zh) * | 2019-08-09 | 2019-10-22 | 哈尔滨工业大学 | 基于多胞折合平衡流形展开模型的航空发动机系统辨识方法 |
| CN111475889A (zh) * | 2020-03-30 | 2020-07-31 | 西安交通大学 | 一种物理与数据驱动的航空发动机数字孪生建模方法 |
| CN112733266A (zh) * | 2020-12-31 | 2021-04-30 | 华南理工大学 | 一种动力总成悬置系统模糊可靠性优化方法 |
| CN114297786A (zh) * | 2021-11-08 | 2022-04-08 | 南京航空航天大学 | 一种涡扇发动机机载动态模型全包线精度扩展的方法 |
| CN115017625A (zh) * | 2022-06-28 | 2022-09-06 | 南京航空航天大学 | 基于基准特性线的变循环发动机变几何部件特性建模方法 |
Non-Patent Citations (2)
| Title |
|---|
| CANCAN WANG: ""Research on LPV-based model of a turbofan engine"", 《2015 12TH INTERNATIONAL BHURBAN CONFERENCE ON APPLIED SCIENCES AND TECHNOLOGY (IBCAST)》, 12 March 2015 (2015-03-12), pages 141 - 145 * |
| 蒋子松: ""涡扇发动机建模及控制计划研究"", 《中国优秀硕士学位论文全文数据库 工程科技II辑》, no. 2021, 15 July 2021 (2021-07-15), pages 031 - 189 * |
Cited By (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN115680903A (zh) * | 2022-08-31 | 2023-02-03 | 中国航发四川燃气涡轮研究院 | 可回收旁路引气控制方法 |
| CN115680903B (zh) * | 2022-08-31 | 2024-05-03 | 中国航发四川燃气涡轮研究院 | 可回收旁路引气控制方法 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN108828947B (zh) | 一种航空发动机含时滞的不确定性模糊动态模型建模方法 | |
| CN106951634B (zh) | 一种航空发动机鲁棒跟踪控制器设计方法 | |
| CN111680357B (zh) | 一种变循环发动机机载实时模型的部件级无迭代构建方法 | |
| Nada | Performance characterization of different configurations of gas turbine engines | |
| CN110502840B (zh) | 航空发动机气路参数在线预测方法 | |
| CN111079290B (zh) | 一种基于改进多元宇宙算法的发动机模型修正方法 | |
| CN109460628B (zh) | 一种进气道与发动机共同工作的流量匹配评估方法 | |
| CN114017201B (zh) | 一种航空燃气涡扇发动机喷口控制方法及装置 | |
| CN112257256B (zh) | 一种基于稳态数据的发动机简化动态模型设计方法 | |
| CN111914367B (zh) | 一种航空发动机部件级模型 | |
| CN114912191A (zh) | 一种导叶可调的涡扇发动机部件级机载动态模型设计方法 | |
| CN115217635B (zh) | 一种涡扇发动机全包线自适应加速控制方法 | |
| CN110321586B (zh) | 一种航空发动机偏离设计点工作状态迭代求解的取值方法 | |
| CN114048554A (zh) | 一种航空发动机三维匹配迭代方法 | |
| CN112668162A (zh) | 一种基于惯性滑模的航空发动机建模方法 | |
| CN113656907B (zh) | 一种航空发动机三维稳态仿真匹配迭代方法 | |
| Zhang et al. | An integrated modeling approach for variable cycle engine performance analysis | |
| Hao et al. | A new design method for mode transition control law of variable cycle engine | |
| CN114297786A (zh) | 一种涡扇发动机机载动态模型全包线精度扩展的方法 | |
| CN110362960B (zh) | 基于多胞折合平衡流形展开模型的航空发动机系统辨识方法 | |
| Zheng et al. | Optimal matching control schedule research on an energy system | |
| Kim | Application of machine learning and its effectiveness in performance model adaptation for a turbofan engine | |
| CN111852663A (zh) | 变循环发动机降保守性鲁棒增益调度控制器 | |
| CN114017200B (zh) | 一种涡扇发动机进口总压重构方法及装置 | |
| Saleh et al. | Modeling and simulation of a double spool turbofan engine using simulink® |
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 |