CN111611743A - 轴流式压气机特性线自适应方法 - Google Patents

轴流式压气机特性线自适应方法 Download PDF

Info

Publication number
CN111611743A
CN111611743A CN202010411884.5A CN202010411884A CN111611743A CN 111611743 A CN111611743 A CN 111611743A CN 202010411884 A CN202010411884 A CN 202010411884A CN 111611743 A CN111611743 A CN 111611743A
Authority
CN
China
Prior art keywords
compressor
axial
stage
characteristic
general
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.)
Granted
Application number
CN202010411884.5A
Other languages
English (en)
Other versions
CN111611743B (zh
Inventor
徐超
吴建平
陈健
钱江
郑隆云
应雨龙
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Shanghai Shangdian Electric Power Engineering Co ltd
Original Assignee
Shanghai Shangdian Electric Power Engineering Co ltd
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Shanghai Shangdian Electric Power Engineering Co ltd filed Critical Shanghai Shangdian Electric Power Engineering Co ltd
Priority to CN202010411884.5A priority Critical patent/CN111611743B/zh
Publication of CN111611743A publication Critical patent/CN111611743A/zh
Application granted granted Critical
Publication of CN111611743B publication Critical patent/CN111611743B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

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/25—Design optimisation, verification or simulation using particle-based methods
    • F—MECHANICAL ENGINEERING; LIGHTING; HEATING; WEAPONS; BLASTING
    • F04—POSITIVE - DISPLACEMENT MACHINES FOR LIQUIDS; PUMPS FOR LIQUIDS OR ELASTIC FLUIDS
    • F04D—NON-POSITIVE-DISPLACEMENT PUMPS
    • F04D19/00—Axial-flow pumps
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F2111/00—Details relating to CAD techniques
    • G06F2111/06—Multi-objective optimisation, e.g. Pareto optimisation using simulated annealing [SA], ant colony algorithms or genetic algorithms [GA]
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F2111/00—Details relating to CAD techniques
    • G06F2111/10—Numerical modelling
    • 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
    • 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

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • General Engineering & Computer Science (AREA)
  • Theoretical Computer Science (AREA)
  • Computer Hardware Design (AREA)
  • Evolutionary Computation (AREA)
  • Geometry (AREA)
  • General Physics & Mathematics (AREA)
  • Mechanical Engineering (AREA)
  • Structures Of Non-Positive Displacement Pumps (AREA)

Abstract

本发明涉及一种轴流式压气机特性线自适应方法，根据轴流式压气机各个级的通用级特性曲线计算公式，通过逐级叠加计算法建立多级轴流式压气机的热力模型；采集实际对象燃气轮机的气路测量参数，以热力模型计算得到的气路参数与实际对象燃气轮机的采集得到气路测量数据之间的均方根误差为目标函数，通过粒子群优化算法迭代优化热力模型中各级通用级特性曲线的形状系数向量，生成与实际对象燃气轮机的压气机特性相匹配的轴流式压气机通用特性线。从而消除同类型不同燃机之间由于制造、安装偏差而引入的不确定度；消除不同干扰及未知初始条件而引入的不确定度；部件特性线生成时逐级叠加计算假设条件所导致的误差,进而来提高压气机部件特性线的准确性。

Description

轴流式压气机特性线自适应方法
技术领域
本发明涉及一种能源与动力工程技术，特别涉及一种轴流式压气机特性线自适应方法。
背景技术
在当前的燃气轮机热力建模技术条件下，热力模型的准确性主要依赖于其部件(压气机和透平)的特性线和工质热物性计算程序的准确程度，尤其是压气机特性线的精度。这些部件特性线实际上需由发动机试车台在不同操作条件下严格的试验获得，或者通过计算流体力学(CFD)数值模拟获取。由于试车台试验费时且昂贵，发动机制造商不可能会获取每一台燃气轮机的部件特性线。因此，制造商通常只会提供给用户一套同一型号燃气轮机的部件特性线。然而对于同一型号的燃气轮机而言，由于制造和组装偏差的原因，其部件特性也会发生一定差别。此外，由于维护、改造、大修等原因，其部件特性也会发生较大的改变。因此，使用同一型号燃气轮机的同一套部件特性线来热力建模时，通常会产生一定程度的计算误差。对于用户，有时因制造商各种原因，甚至无法获得相关型号燃气轮机的部件特性线，只能通过已有的其他类型燃气轮机的部件特性线进行比例缩放后来使用，致使热力计算误差有时会难以接受。此时，逐级叠加计算法成为了一种用于生成压气机特性线的可靠而有效的手段，其计算过程主要基于一维平均半径处的连续性流动方程和一组通用的级特性曲线。对于用户来说，实际机组压气机各叶片级的几何参数通常也是未知的，因此，逐级叠加计算法通常只能采用一套从大量已有叶片级试验数据中拟合得到的通用的级特性曲线(如图1所示压气机通用级特性曲线分别阐述了压气机级相对压力系数ψ*(ψ*＝ψ/ψ0)和相对流量系数φ*(φ*＝φ/φ0＝(Ca/u)/(Ca/u)0)以及相对等熵效率η*(η*＝η/η0)和相对流量系数φ*的关系)来表征实际压气机各级的级特性,其中ψ为压力系数；φ为流量系数；Ca为气流轴向速度；u为叶片平均半径出的圆周速度；下标0为设计工况点。由于上述压气机特性线生成过程的简化条件，通过逐级叠加计算法建立的燃气轮机热力模型一定程度上存在不可避免的误差，导致部件特性线精度不够。
发明内容
本发明是针对现有技术生成压气机特性线精度不够的问题，提出了一种轴流式压气机特性线自适应方法，以实际燃气轮机的气路测量数据与热力模型计算的气路参数之间的均方根误差为目标函数，通过粒子群优化算法迭代寻优计算得到压气机各个级的级特性曲线的一组形状系数向量，生成与实际燃气轮机的压气机特性相匹配的轴流式压气机通用特性线，从而消除部件特性线生成时逐级叠加计算假设条件所导致的误差。
本发明的技术方案为：一种轴流式压气机特性线自适应方法，首先在[-0.5,1]的数值范围内随机产生一组各级通用级特性曲线的形状系数向量[SF1,SF2,...,SFm]，m表示压气机的级数；再根据轴流式压气机各个级的通用级特性曲线计算公式，通过逐级叠加计算法建立多级轴流式压气机的热力模型；采集实际对象燃气轮机的气路测量参数，以热力模型计算得到的气路参数与实际对象燃气轮机的采集得到气路测量数据之间的均方根误差为目标函数，通过粒子群优化算法迭代优化热力模型中各级通用级特性曲线的形状系数向量[SF1,SF2,...,SFm]，最终通过最优的形状系数向量生成与实际对象燃气轮机的压气机特性相匹配的轴流式压气机通用特性线。
所述各级通用级特性曲线指各叶片级的通用级特性曲线，包括亚音速级、跨音速级和超音速级。
本发明的有益效果在于：本发明轴流式压气机特性线自适应方法，消除部件特性线生成时逐级叠加计算假设条件所导致的误差，来提高压气机部件特性线的准确性，从而提高所建燃气轮机热力模型的准确性。
附图说明
图1为压气机通用的级特性曲线图；
图2为本发明轴流式压气机特性线自适应方法示意图；
图3为本发明自适应后生成的压气机流量特性通用特性线图；
图4为本发明自适应后生成的压气机效率特性通用特性线图。
具体实施方式
逐级叠加计算法(Stage-stacking)是一种用于生成压气机特性线的可靠而有效的手段，其计算过程主要基于一维平均半径处的连续性流动方程和一组通用级特性曲线。对用户来说，实际机组压气机各叶片级的几何参数通常也是未知的，因此，逐级叠加计算法通常只能采用一套从大量已有叶片级试验数据中拟合得到的通用级特性曲线来表征实际压气机各级的级特性，而忽略了各个不同级类型叶片级(如亚音速级、跨音速级和超音速级)本身具有的不同级特性。本发明以实际燃气轮机的气路测量数据与热力模型计算的气路参数之间的均方根误差为目标函数，通过粒子群优化算法迭代寻优计算得到压气机各个级的级特性曲线的形状系数向量SF＝[SF1,SF2,...,SFi,...,SFm](m表示压气机的级数)，生成与实际燃气轮机的压气机特性相匹配的轴流式压气机通用特性线，从而消除部件特性线生成时逐级叠加计算假设条件所导致的误差，来提高压气机部件特性线的准确性，从而提高所建燃气轮机热力模型的准确性。
具体的技术方案如下：
如图2所示，首先在[-0.5,1]的数值范围内随机产生一组各级通用级特性曲线的形状系数向量[SF1,SF2,...,SFm]，再根据轴流式压气机各个级的通用级特性曲线计算公式(1)与(2)通过逐级叠加计算法建立多级轴流式压气机的热力模型。采集实际对象燃气轮机的气路测量参数，以热力模型计算得到的气路参数与实际对象燃气轮机的采集得到气路测量数据之间的均方根误差为目标函数，通过粒子群优化算法迭代优化热力模型中各级通用级特性曲线的形状系数向量[SF1,SF2,...,SFm]，最终通过最优的形状系数向量生成与实际对象燃气轮机的压气机特性相匹配的轴流式压气机通用特性线(流量特性线与效率特性线)，可以消除3方面的不确定度：①消除同类型不同燃机之间由于制造、安装偏差而引入的不确定度；②消除不同干扰及未知初始条件而引入的不确定度；③部件特性线生成时逐级叠加计算假设条件所导致的误差。从而来提高压气机部件特性线的准确性，从而提高所建燃气轮机热力模型的准确性。
轴流式压气机各个级的通用级特性曲线可用如式(1)和(2)表示：
Figure BDA0002493562180000041
Figure BDA0002493562180000042
式中
Figure BDA0002493562180000043
为压气机i级相对压力系数；
Figure BDA0002493562180000044
为压气机i级相对流量系数。由级的压力系数的定义
Figure BDA0002493562180000045
和级的流量系数的定义φ＝Ca/u，通过逐级叠加计算法可得某一级的相对压力系数为
Figure BDA0002493562180000046
和该级的相对流量系数
Figure BDA0002493562180000047
这里
Figure BDA0002493562180000048
表示i级的入口比焓，
Figure BDA0002493562180000049
表示i级的出口理想比焓，Rg表示空气的气体常数，
Figure BDA00024935621800000410
表示i级的入口总温，
Figure BDA00024935621800000411
表示i级的入口总压，Gi,in表示i级的入口空气质量流量；下角标0表示设计条件下。
SFi为i级的级特性曲线的形状系数，通过调整SFi的设定值可以表征各个不同级类型(如亚音速级、跨音速级和超音速级)，这里作为待优化的变量。
其中，式(1)的具体展开式为：
Figure BDA00024935621800000412
式中，
Figure BDA00024935621800000413
为压气机第i级相对压力系数的最大值；
Figure BDA00024935621800000414
为取压气机第i级相对压力系数的最大值时对应的相对流量系数。
式(2)的具体展开式为：
当
Figure BDA0002493562180000051
时，
Figure BDA0002493562180000052
当
Figure BDA0002493562180000053
时，
Figure BDA0002493562180000054
式中，
Figure BDA0002493562180000055
为
Figure BDA0002493562180000056
取最小值时对应的相对等熵效率；
Figure BDA0002493562180000057
为
Figure BDA0002493562180000058
取最大值时对应的相对等熵效率。这里
Figure BDA0002493562180000059
取为0.04，对应的
Figure BDA00024935621800000510
取为0.20；
Figure BDA00024935621800000511
取为1.46，对应的
Figure BDA00024935621800000512
取为0.92。
基于上述方程(1)至(5)建立压气机热力模型，其中压气机热力模型的整机通用特性线可用如式(6)和(7)表示：
Gcor,rel＝f1(ncor,rel,πC,rel) (6)
ηC,rel＝f2(ncor,rel,πC,rel) (7)
式中：
Figure BDA00024935621800000513
为压气机的相对折合转速；
Figure BDA00024935621800000514
为压气机的相对折合流量；
Figure BDA00024935621800000515
为压气机的相对压比；ηC,rel＝ηC/ηC0为压气机的相对等熵效率；下角标0表示设计工况；n表示转速；
Figure BDA00024935621800000516
为压气机入口总温；Rg为空气的气体常数；G为空气的质量流量；
Figure BDA00024935621800000517
为压气机入口总压；πC为压气机压比。
再以实际燃气轮机的气路测量数据(包括对象压气机的进出口总温、总压、转速以及大气温度、压力及相对湿度)与所建压气机热力模型计算的相应气路参数之间的均方根误差为目标函数，通过粒子群优化算法迭代寻优计算得到压气机各个级的级特性曲线的形状系数向量[SF1,SF2,...,SFi,...,SFm]，生成与实际燃气轮机的压气机特性相匹配的轴流式压气机通用特性线，从而消除部件特性线生成时逐级叠加计算假设条件所导致的误差，来提高压气机部件特性线的准确性，从而提高所建燃气轮机热力模型的准确性。
其中，粒子群优化算法的初始迭代的形状系数的取值范围为[-0.5,1]，m为轴流式压气机的总级数，最终通过迭代寻优，以某型燃气轮机为例，可以生成如图3和图4所示的自适应后的压气机通用特性曲线簇。

Claims (2)

1.一种轴流式压气机特性线自适应方法，其特征在于，首先在[-0.5,1]的数值范围内随机产生一组各级通用级特性曲线的形状系数向量[SF1,SF2,...,SFm]，m表示压气机的级数；再根据轴流式压气机各个级的通用级特性曲线计算公式，通过逐级叠加计算法建立多级轴流式压气机的热力模型；采集实际对象燃气轮机的气路测量参数，以热力模型计算得到的气路参数与实际对象燃气轮机的采集得到气路测量数据之间的均方根误差为目标函数，通过粒子群优化算法迭代优化热力模型中各级通用级特性曲线的形状系数向量[SF1,SF2,...,SFm]，最终通过最优的形状系数向量生成与实际对象燃气轮机的压气机特性相匹配的轴流式压气机通用特性线。
2.根据权利要求1所述轴流式压气机特性线自适应方法，其特征在于，所述各级通用级特性曲线指各叶片级的通用级特性曲线，包括亚音速级、跨音速级和超音速级。
CN202010411884.5A 2020-05-15 2020-05-15 轴流式压气机特性线自适应方法 Active CN111611743B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202010411884.5A CN111611743B (zh) 2020-05-15 2020-05-15 轴流式压气机特性线自适应方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202010411884.5A CN111611743B (zh) 2020-05-15 2020-05-15 轴流式压气机特性线自适应方法

Publications (2)

Publication Number Publication Date
CN111611743A true CN111611743A (zh) 2020-09-01
CN111611743B CN111611743B (zh) 2023-03-28

Family

ID=72200348

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202010411884.5A Active CN111611743B (zh) 2020-05-15 2020-05-15 轴流式压气机特性线自适应方法

Country Status (1)

Country Link
CN (1) CN111611743B (zh)

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN114444331A (zh) * 2022-04-07 2022-05-06 中国航发四川燃气涡轮研究院 一种多级轴流压气机的级特性匹配方法
CN118377921A (zh) * 2024-06-21 2024-07-23 浙江大学 一种燃气轮机部件全工况特性曲线的自适应生成方法

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
DE102009033754A1 (de) * 2009-07-17 2011-01-20 Rolls-Royce Deutschland Ltd & Co Kg Axialverdichter mit einem Strömungsimpulserzeuger
CN103294892A (zh) * 2013-04-01 2013-09-11 深圳市广前电力有限公司 基于燃气轮机温控作用下的压气机在线特性分析方法
CN106126902A (zh) * 2016-06-21 2016-11-16 上海电机学院 基于粒子群优化算法辨识的燃气轮机部件特性线修正方法

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
DE102009033754A1 (de) * 2009-07-17 2011-01-20 Rolls-Royce Deutschland Ltd & Co Kg Axialverdichter mit einem Strömungsimpulserzeuger
CN103294892A (zh) * 2013-04-01 2013-09-11 深圳市广前电力有限公司 基于燃气轮机温控作用下的压气机在线特性分析方法
CN106126902A (zh) * 2016-06-21 2016-11-16 上海电机学院 基于粒子群优化算法辨识的燃气轮机部件特性线修正方法

Non-Patent Citations (4)

* Cited by examiner, † Cited by third party
Title
GUOQIANG ZHANG 等: ""Thermodynamic analysis of combined cycle under design/off-design conditions for its efficient design and operation"", 《ENERGY CONVERSION AND MANAGEMENT》 *
应雨龙 等: ""基于热力模型的燃气轮机气路故障预测诊断研究综述"", 《中国电机工程学报》 *
潘鹏飞 等: ""基于遗传算法的发动机部件级模型设计点匹配"", 《2013 年首届中国航空科学技术大会论文集》 *
董桢: ""基于PSO算法的涡轴发动机部件特性修正方法研究"", 《中国优秀硕士学位论文全文数据库 (工程科技Ⅱ辑)》 *

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN114444331A (zh) * 2022-04-07 2022-05-06 中国航发四川燃气涡轮研究院 一种多级轴流压气机的级特性匹配方法
CN118377921A (zh) * 2024-06-21 2024-07-23 浙江大学 一种燃气轮机部件全工况特性曲线的自适应生成方法

Also Published As

Publication number Publication date
CN111611743B (zh) 2023-03-28

Similar Documents

Publication Publication Date Title
CN111914362B (zh) 一种研发阶段涡扇发动机模型自适应方法
CN110717219B (zh) 航空发动机整机状态下压气机进口流量的获取方法及装置
CN108829928B (zh) 一种涡轴发动机自适应部件级仿真模型构建方法
CN109684597B (zh) 一种用于压气机全工况特性拓展的方法
CN111611743B (zh) 轴流式压气机特性线自适应方法
CN107203364B (zh) 一种用于压气机全工况特性的预测和辨识方法
CN106126902A (zh) 基于粒子群优化算法辨识的燃气轮机部件特性线修正方法
CN114577484A (zh) 一种核心机试验性能修正方法
CN110321586B (zh) 一种航空发动机偏离设计点工作状态迭代求解的取值方法
CN115903484A (zh) 基于合作博弈的航空发动机多变量鲁棒控制器优化方法
CN117725700A (zh) 基于数字孪生技术的分轴燃气轮机管理系统、方法及设备
CN114154234A (zh) 一种航空发动机建模方法、系统、存储介质
CN116663202B (zh) 一种多级轴流压气机性能仿真模型的校核方法
Li et al. Improved multiple point non-linear genetic algorithm based performance adaptation using least square method
CN114088409B (zh) 一种燃气-蒸汽联合循环单轴机组部分负荷性能试验方法
CN114818205A (zh) 一种航空发动机全生命周期叶尖间隙在线感知方法
CN114112414B (zh) 一种燃气-蒸汽联合循环分轴机组部分负荷性能试验方法
Serrano et al. Methodology to evaluate turbocharger turbine performance at high blade to jet speed ratio under quasi adiabatic conditions
CN114491417A (zh) 一种基于cdfs模态变化性能的一维输入修正方法
Zhang et al. Operational data-based adaptive improvement method of gas turbine component characteristics for performance simulation
Bruni et al. Manufacturing and Build Variations Modelling for Multi-Stage Axial Compressors: CFD Aerodynamic Predictions
Trawiński Development of flow and efficiency characteristics of an axial compressor with an analytical method including cooling air extraction and variable inlet guide vane angle
CN115688554B (zh) 一种涡轴发动机旋转部件低转速特性外推及修正方法
CN112834226B (zh) 一种航空发动机气路在线监测传感器优化配置方法
Petrovic et al. New Method for Cycle Performance Prediction Based on Detailed Compressor and Gas Turbine Flow Calculations

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