CN110929397B - 基于零转速特性线的燃气轮机涡轮低转速特性线外推方法 - Google Patents
基于零转速特性线的燃气轮机涡轮低转速特性线外推方法 Download PDFInfo
- Publication number
- CN110929397B CN110929397B CN201911144121.2A CN201911144121A CN110929397B CN 110929397 B CN110929397 B CN 110929397B CN 201911144121 A CN201911144121 A CN 201911144121A CN 110929397 B CN110929397 B CN 110929397B
- Authority
- CN
- China
- Prior art keywords
- function
- turbine
- rotation
- characteristic line
- zero
- 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
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F18/00—Pattern recognition
- G06F18/20—Analysing
- G06F18/21—Design or setup of recognition systems or techniques; Extraction of features in feature space; Blind source separation
- G06F18/214—Generating training patterns; Bootstrap methods, e.g. bagging or boosting
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F18/00—Pattern recognition
- G06F18/20—Analysing
- G06F18/24—Classification techniques
- G06F18/241—Classification techniques relating to the classification model, e.g. parametric or non-parametric approaches
- G06F18/2411—Classification techniques relating to the classification model, e.g. parametric or non-parametric approaches based on the proximity to a decision surface, e.g. support vector machines
-
- 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)
- Data Mining & Analysis (AREA)
- Theoretical Computer Science (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Bioinformatics & Cheminformatics (AREA)
- Bioinformatics & Computational Biology (AREA)
- Artificial Intelligence (AREA)
- Evolutionary Biology (AREA)
- Evolutionary Computation (AREA)
- Physics & Mathematics (AREA)
- General Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Life Sciences & Earth Sciences (AREA)
- Control Of Turbines (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
本发明公开了一种基于零转速特性线的燃气轮机涡轮低转速特性线外推方法。它包括如下步骤,步骤一:对涡轮特性线进行描述改进;步骤二:将涡轮零转速特性数据和已知涡轮高转速特征数据作为样本数据;步骤三:结合样本数据,利用支持向量机拟合的方法对燃气轮机涡轮低转速特性线进行外推。本发明具有提高涡轮特性线外推精度的优点。
Description
技术领域
本发明涉及燃气轮机涡轮技术领域,更具体地说它是一种基于零转速线的燃气轮机涡轮低转速特性线外推方法。
背景技术
在燃气轮机的起动过程建模研究中,一个非常复杂的问题是燃气轮机的低转速部件特性难以获取,主要原因是由于低转速的部件特性受外界随机因素影响较大,并且低转速区域部件内部流场的气流速度和气流攻角与设计值差别较大,因此在试验和数值计算中都较难获取。考虑到燃气轮机高转速下的部件特性较容易获取,因此,燃气轮机低转速部件特性的获取通常是采用依据高转速部件特性数据进行拟合外推的方法。
国内对燃气轮机部件特性的描述一般采用传统描述方式,即用压比、转速、效率、流量四个变量来表征部件特性。但由于传统描述方式得到的特性线具有非线性,因而,由传统描述方式结合高转速部件特征数据进行外推而得到的低转速部件特性并不能达到满意精度。刘建勋博士、李应红院士(2007年2月《发动机起动过程部件特性的一种改进描述》,燃气涡轮试验与研究)提出了一种新的部件特性描述方式,用换算扭矩代替效率,并用出口条件表示的换算转速代替常规的换算转速,使得改进后的部件特性线性程度更高。但是根据改进后的部件特性描述方式,根据已知的高转速涡轮部件特性数据结合支持向量机的外推方法得到的低转速部件特性数据精度不高。
因此,现亟需开发一种提高涡轮特性线外推精度的燃气轮机涡轮低转速特性线改进方法。
发明内容
本发明的目的是为了提供一种基于零转速特性线的燃气轮机涡轮低转速特性线外推方法,在采用支持向量机和涡轮特性线转换描述方法相结合的基础上,通过在支持向量机的训练集中加入零转速线的特征数据,提高涡轮特性线外推精度。
为了实现上述目的,本发明的技术方案为:基于零转速特性线的燃气轮机涡轮低转速特性线外推方法,其特征在于:包括如下步骤,
步骤一:对涡轮特性线进行描述改进;
步骤二:将涡轮零转速特性数据和已知涡轮高转速特征数据作为样本数据;
步骤三:结合样本数据,利用支持向量机拟合的方法对燃气轮机涡轮低转速特性线进行外推。
在上述技术方案中,在步骤一中,对涡轮特性线进行描述的改进,涡轮特性线用换算转速ntnp和膨胀比πt来表示换算流量Gtnp和换算扭矩Tqnp,改进后的涡轮特性线描述为:
式(1-1)中:ntnp为换算转速;
πt为膨胀比;
Gtnp为换算流量;
Tqnp为换算扭矩;
f1,f2为两个未知的函数。
在上述技术方案中,在步骤二中,由于涡轮零转速线必然过压比为1且流量为0的点,结合涡轮高转速的特性数据,将零转速特征点作为训练数据的一部分。
在上述技术方案中,在步骤三中,在进行拟合时,首先通过非线性变换x→Ψ(x)将输入空间映射到一个高维的特征空间,然后依据统计学习理论中的结构风险最小化准则,在特征空间进行最优线性逼近,即:
b是某未知量,类似线性函数的截距;
w=(w1,…,wN)是确定了这个空间中的一个超平面的权值向量;
式(1-4)中,|| ||表示向量模;
C为平衡因子;
L为惩罚函数;
J为风险函数;
xi,yi表示第i个样本数据;
w=(w1,…,wN)是确定了这个空间中的一个超平面的权值向量;
惩罚函数L通常取线性ε不敏感损失函数,定义为:
在上述技术方案中,将所述式(1-4)的最小风险函数等价为如下的二次规划问题:
式(1-6)和(1-7)中,核函数K(xi,yi)=(ψ(xi)·ψ(yi))是满足Mercer条件的任意对称函数;
xi,yi为样本数据;
J为风险函数;
i为样本点序列;
ε为容忍误差;
最后,支持向量机的输出为:
K(xi,yi)=(ψ(xi)·ψ(yi))是满足Mercer条件的任意对称函数;
b是某未知量,类似线性函数的截距。
上述的特性线描述改进、支持向量机外推均为已知技术。
上述的KKT条件为现有技术,KKT条件是解决最优化问题时用到的一种方法;最优化问题通常是指对于给定的某一函数,求其在指定作用域上的全局最小值。KKT条件和拉格朗日乘子均是求解最优化问题的方法,不同之处在于英语的情形不同。
本发明具有如优点:
本发明针对燃气轮机涡轮特性线外推精度不高的问题,在采用支持向量机和涡轮特性线转换描述方法相结合的基础上,通过在支持向量机的训练集中加入零转速线的特征数据,从而提高涡轮特性线外推精度。
附图说明
图1为两种外推方法对涡轮低转速特性外推的对比图。
图1中,仿真流量为支持向量机训练集没有增加零转速线的外推结果;基于零转速线的仿真流量为支持向量机训练集增加零转速线的外推结果。
图1中,涡轮试验流量表示试验测得的流量数据;
本发明采用支持向量机和部件特性转换描述相结合的方法,对涡轮低转速特性线进行外推,申请人发现得到的涡轮低转速部件特性线精确度不高的原因是支持向量机训练集中低转速部件特性数据不足。因此本发明通过在支持向量机的训练集中加入零转速线的特征数据,从而扩充了支持向量机训练集在涡轮低转速的特征数据,使涡轮低转速特性线外推精度更高。
具体实施方式
下面结合附图详细说明本发明的实施情况,但它们并不构成对本发明的限定,仅作举例而已。同时通过说明使本发明的优点更加清楚和容易理解。
参阅附图可知:基于零转速特性线的燃气轮机涡轮低转速特性线外推方法,包括如下步骤,
步骤一:对涡轮特性线进行描述进行改进;
步骤二:将涡轮零转速特性数据和已知涡轮高转速特征数据作为样本数据;
步骤三:方程左边的变量用y表示,方程右边函数参数组成的向量用x表示,使y与x的关系满足支持向量机拟合;
步骤四:在向量机进行拟合时,首先通过非线性变换将输入空间映射到一个高维的特征空间,然后依据统计学习理论中的结构风险最小化准则,在特征空间进行最优线性逼近;
步骤五:将最小风险函数等价为二次规划问题;
步骤六:结合样本数据,利用支持向量机拟合的方法对燃气轮机涡轮低转速特性线进行外推。
进一步地,在步骤一中,涡轮特性线一般用换算转速ntnp和膨胀比πt来表示换算流量Gtnp和效率ηt,对涡轮特性线进行描述的改进,用换算扭矩来Tqnp代替效率,改进后的涡轮特性线描述为:
式(1-1)中:ntnp为换算转速;
πt为膨胀比;
Gtnp为换算流量;
Tqnp为换算扭矩;
f1,f2为两个未知的函数。
进一步地,在步骤二中,由于压气机和涡轮零转速线必然过压比为1且流量为0的点,结合涡轮高转速的样本数据,将零转速特征点作为训练数据的一部分。
进一步地,在步骤三中,方程左边的变量用y表示,方程右边函数参数组成的向量用x表示,支持向量机拟合的思想为:针对任意一个未知非线性函数y=f(x),以及对其采样得到的样本集合寻求函数/>使得函数f和/>之间的广义距离/>最小,即满足如下(1-2)所示的方程式,
L是惩罚函数;
P为联合概率;
f表示某未知函数;
x,y均表示样本数据。
进一步地,在步骤四中,在进行拟合时,首先通过非线性变换x→Ψ(x)将输入空间映射到一个高维的特征空间,然后依据统计学习理论中的结构风险最小化准则,在特征空间进行最优线性逼近,即:
b是某未知量,类似线性函数的截距;
w=(w1,…,wN)是确定了这个空间中的一个超平面的权值向量;
式(1-4)中,|| ||表示向量模;
C为平衡因子;
L为惩罚函数;
J为风险函数;
xi,yi表示第i个样本数据;
w=(w1,…,wN)是确定了这个空间中的一个超平面的权值向量;
惩罚函数L通常取线性ε不敏感损失函数,定义为:
更进一步地,在步骤五中,所述式(1-4)的最小风险函数等价为如下的二次规划问题:
其中,式(1-6)为目标函数;式(1-7)为约束条件;
式(1-6)和(1-7)中,核函数K(xi,yi)=(ψ(xi)·ψ(yi))是满足Mercer条件的任意对称函数,如多项式核[(u·v)+1]d,径向基核exp(-||u-v||2/2σ2)等(其中,u表示xi,v表示yi),xi,yi为样本数据;
J为风险函数;
i为样本点序列;
ε为容忍误差;
最后,支持向量机的输出为:
K(xi,yi)=(ψ(xi)·ψ(yi))是满足Mercer条件的任意对称函数;
b是某未知量,类似线性函数的截距。
实施例
选取涡轮转换转速为1、0.9431、0.8431、0.7349、0.6257时的特征数据,组成训练集样本(xi,yi),其中xi为[ntnp,πt]T,yi为Gtnp。选取径向基核宽度σ为0.1,平衡因子C为100,不敏感因子为0.0001,对涡轮转换转速为0.554的折合流量进行外推。
训练集样本如下表1所示:
表1训练集样本
结果如图1仿真流量所示,在外推燃气轮机部件特性时需要考虑转速线与零转速特性线的相似性,因此需要将部件的零转速特性考虑在内。
由于压气机和涡轮零转速线必然过压比为1且流量为0的点,结合涡轮高转速的特性数据,将零转速特征点作为训练数据的一部分,重新进行支持向量机训练,训练结果如图1基于零转速线的仿真流量所示。
由试验测得的涡轮转换转速为0.554的折合流量数据如图1涡轮试验流量所示。在增加零转速特征点前,外推得到的涡轮低转速特性线与实验得到的特性线的相对误差最大为10%,增加零转速特征点后,特性线的相对误差最大为0.57%,以上分析表明,零转速特征点数据的加入明显提升了外推效果。
其它未说明的部分均属于现有技术。
Claims (3)
1.基于零转速特性线的燃气轮机涡轮低转速特性线外推方法,其特征在于:包括如下步骤,
步骤一:对涡轮特性线进行描述改进;
在步骤一中,对涡轮特性线进行描述的改进,涡轮特性线用换算转速ntnp和膨胀比πt来表示换算流量Gtnp和换算扭矩Tqnp,改进后的涡轮特性线描述为:
式(1-1)中:ntnp为换算转速;
πt为膨胀比;
Gtnp为换算流量;
Tqnp为换算扭矩;
f1,f2为两个未知的函数;
步骤二:将涡轮零转速特性数据和已知涡轮高转速特征数据作为样本数据;
步骤三:结合样本数据,利用支持向量机拟合的方法对燃气轮机涡轮低转速特性线进行外推;
在步骤三中,在进行拟合时,首先通过非线性变换x→Ψ(x)将输入空间映射到一个高维的特征空间,然后依据统计学习理论中的结构风险最小化准则,在特征空间进行最优线性逼近,即:
b是某未知量,类似线性函数的截距;
w=(w1,…,wN)是确定了这个空间中的一个超平面的权值向量;
式(1-4)中,||||表示向量模;
C为平衡因子;
L为惩罚函数;
J为风险函数;
xi,yi表示第i个样本数据;
w=(w1,…,wN)是确定了这个空间中的一个超平面的权值向量;
惩罚函数L通常取线性ε不敏感损失函数,定义为:
2.根据权利要求1所述的基于零转速特性线的燃气轮机涡轮低转速特性线外推方法,其特征在于:在步骤二中,由于涡轮零转速线必然过压比为1且流量为0的点,结合涡轮高转速的特性数据,将零转速特征点作为训练数据的一部分。
3.根据权利要求2所述的基于零转速特性线的燃气轮机涡轮低转速特性线外推方法,其特征在于:将所述式(1-4)的最小风险函数等价为如下的二次规划问题:
式(1-6)和(1-7)中,核函数K(xi,yi)=(ψ(xi)·ψ(yi))是满足Mercer条件的任意对称函数;
xi,yi为样本数据;
J为风险函数;
i为样本点序列;
ε为容忍误差;
最后,支持向量机的输出为:
K(xi,yi)=(ψ(xi)·ψ(yi))是满足Mercer条件的任意对称函数;
b是某未知量,类似线性函数的截距。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201911144121.2A CN110929397B (zh) | 2019-11-20 | 2019-11-20 | 基于零转速特性线的燃气轮机涡轮低转速特性线外推方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201911144121.2A CN110929397B (zh) | 2019-11-20 | 2019-11-20 | 基于零转速特性线的燃气轮机涡轮低转速特性线外推方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN110929397A CN110929397A (zh) | 2020-03-27 |
CN110929397B true CN110929397B (zh) | 2023-06-30 |
Family
ID=69851376
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201911144121.2A Active CN110929397B (zh) | 2019-11-20 | 2019-11-20 | 基于零转速特性线的燃气轮机涡轮低转速特性线外推方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN110929397B (zh) |
Families Citing this family (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN115688554B (zh) * | 2022-08-24 | 2024-04-30 | 南京航空航天大学 | 一种涡轴发动机旋转部件低转速特性外推及修正方法 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105930662A (zh) * | 2016-04-26 | 2016-09-07 | 中国科学院工程热物理研究所 | 一种压气机低转速特性外推方法 |
CN106569977A (zh) * | 2016-09-26 | 2017-04-19 | 南京航空航天大学 | 涡轴发动机起动过程建模方法 |
CN107882636A (zh) * | 2017-10-26 | 2018-04-06 | 中国航发湖南动力机械研究所 | 航空发动机压气机低转速部件特性扩展方法和装置 |
CN108304652A (zh) * | 2018-02-01 | 2018-07-20 | 中国人民解放军海军航空大学 | 一种基于改进qpso-elm的涡轴发动机起动过程模型辨识方法 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20180107176A1 (en) * | 2016-10-18 | 2018-04-19 | United Technologies Corporation | System and method for simulating an operation of a gas turbine engine |
-
2019
- 2019-11-20 CN CN201911144121.2A patent/CN110929397B/zh active Active
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105930662A (zh) * | 2016-04-26 | 2016-09-07 | 中国科学院工程热物理研究所 | 一种压气机低转速特性外推方法 |
CN106569977A (zh) * | 2016-09-26 | 2017-04-19 | 南京航空航天大学 | 涡轴发动机起动过程建模方法 |
CN107882636A (zh) * | 2017-10-26 | 2018-04-06 | 中国航发湖南动力机械研究所 | 航空发动机压气机低转速部件特性扩展方法和装置 |
CN108304652A (zh) * | 2018-02-01 | 2018-07-20 | 中国人民解放军海军航空大学 | 一种基于改进qpso-elm的涡轴发动机起动过程模型辨识方法 |
Non-Patent Citations (1)
Title |
---|
周超等.《两种低转速特性外推方法的对比》.《工程与试验》.2019,(第undefined期),第11-14页. * |
Also Published As
Publication number | Publication date |
---|---|
CN110929397A (zh) | 2020-03-27 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN105512429B (zh) | 一种三轴燃气轮机的总体方案计算方法 | |
CN109684597B (zh) | 一种用于压气机全工况特性拓展的方法 | |
CN110929357A (zh) | 一种高性能舰船燃机压气机气动设计方法 | |
CN109992880B (zh) | 一种燃气轮机性能适配方法、设备以及系统 | |
CN111027148B (zh) | 损失落后角模型自动标定及工业轴流压缩机性能计算方法 | |
CN112270139B (zh) | 基于母型库的燃料电池离心压缩机气动优化设计方法 | |
CN110929397B (zh) | 基于零转速特性线的燃气轮机涡轮低转速特性线外推方法 | |
CN105787151A (zh) | 一种2.5维陶瓷基复合材料结构可靠性稳健优化设计方法 | |
CN114254460A (zh) | 基于图卷积神经网络的透平机械气动鲁棒性优化方法 | |
CN109902382B (zh) | 一种大中型混流式水轮机比转速的统计计算方法 | |
CN110705079B (zh) | 一种基于模拟退火算法的离心式压缩机结构优化方法 | |
Wisler | Core compressor exit stage study. Volume 1: Blading design | |
CN116663202B (zh) | 一种多级轴流压气机性能仿真模型的校核方法 | |
Fei et al. | Application of new empirical models based on mathematical statistics in the through-flow analysis | |
CN116595874A (zh) | 叶轮机械性能预测模型参数优化方法及装置、存储介质 | |
Calvert et al. | Design, test and analysis of a high-pressure-ratio transonic fan | |
CN116702511A (zh) | 一种可调导叶落后角的计算方法 | |
KR102082909B1 (ko) | 풍력 발전기의 후류모델 보정방법 | |
Chu et al. | Improved fast model migration method for centrifugal compressor based on bayesian algorithm and Gaussian process model | |
Lee et al. | Multi-disciplinary design optimization and performance evaluation of a single stage transonic axial compressor | |
CN109948299A (zh) | 一种压气机叶片中弧线计算方法 | |
Venter et al. | Benchmarking the Performance of the Actuator-Disk Method for Low-Pressure Axial Flow Fan Simulation | |
Purushothaman et al. | Cold blade profile generation methodology for compressor rotor blades using FSI approach | |
Maksiuta et al. | Deviation Angle in a Turbine Nozzle Cascade With Convergent Meridional Shape of Flow Path | |
Liu et al. | Study on component characteristic correction algorithm of turbofan engine |
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 |