CN112947050A - 一种燃气轮机仿真动态性能控制律参数确定方法 - Google Patents

一种燃气轮机仿真动态性能控制律参数确定方法 Download PDF

Info

Publication number
CN112947050A
CN112947050A CN202110121594.1A CN202110121594A CN112947050A CN 112947050 A CN112947050 A CN 112947050A CN 202110121594 A CN202110121594 A CN 202110121594A CN 112947050 A CN112947050 A CN 112947050A
Authority
CN
China
Prior art keywords
gas turbine
dynamic performance
control law
parameters
performance control
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
Application number
CN202110121594.1A
Other languages
English (en)
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.)
AECC Shenyang Engine Research Institute
Original Assignee
AECC Shenyang Engine Research Institute
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 AECC Shenyang Engine Research Institute filed Critical AECC Shenyang Engine Research Institute
Priority to CN202110121594.1A priority Critical patent/CN112947050A/zh
Publication of CN112947050A publication Critical patent/CN112947050A/zh
Pending legal-status Critical Current

Links

Images

Classifications

    • G—PHYSICS
    • G05—CONTROLLING; REGULATING
    • G05B—CONTROL OR REGULATING SYSTEMS IN GENERAL; FUNCTIONAL ELEMENTS OF SUCH SYSTEMS; MONITORING OR TESTING ARRANGEMENTS FOR SUCH SYSTEMS OR ELEMENTS
    • G05B11/00—Automatic controllers
    • G05B11/01—Automatic controllers electric
    • G05B11/36—Automatic controllers electric with provision for obtaining particular characteristics, e.g. proportional, integral, differential
    • G05B11/42—Automatic controllers electric with provision for obtaining particular characteristics, e.g. proportional, integral, differential for obtaining a characteristic which is both proportional and time-dependent, e.g. P. I., P. I. D.

Landscapes

  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Engineering & Computer Science (AREA)
  • Automation & Control Theory (AREA)
  • Feedback Control In General (AREA)

Abstract

本申请属于燃气轮机动态性能控制律参数确定技术领域，具体涉及一种燃气轮机仿真动态性能控制律参数确定方法，包括：基于设定的燃气轮机动态性能控制律参数种群，对应计算得到燃气轮机动态性能参数；基于计算得到的燃气轮机动态性能参数，对应计算得到与实际燃气轮机动态性能参数间的适应度；依据选择概率淘汰较差适应度对应的燃气轮机动态性能控制律参数种群，对剩余的燃气轮机动态性能控制律参数种群进行交叉、变异，产生燃气轮机动态性能控制律参数子种群，以子种群更新设定的燃气轮机动态性能控制律参数种群，直至达到遗传代数；选取剩余燃气轮机动态性能控制律参数种群中的燃气轮机动态性能控制律参数。

Description

一种燃气轮机仿真动态性能控制律参数确定方法
技术领域
本申请属于燃气轮机动态性能控制律参数确定技术领域，具体涉及一种燃气轮机仿真动态性能控制律参数确定方法。
背景技术
燃气轮机是重要的动力部件，对燃气轮机运行过程进行仿真，模拟燃气轮机实际工作过程，可方便的获取燃气轮机全流程参数，包括实际试车难以测量参数和表征燃气轮机危险性指标的参数，基于此可总结出合理的控制规律，指导燃气轮机设计与试车。
燃气轮机控制系统采用串级闭环形式，如图1所示，包括外控制回路、内控制回路，其中外控制回路为动力涡轮闭环，内控制回路为高压换算转速闭环，内控制回路、外控制回路由经典PID控制律实现，控制律参数均包括比例系数、积分系数和微分系数。
对燃气轮机运行过程进行仿真时，控制系统控制律参数准确设定是仿真结果有效的重要保证，当前，对燃气轮机运行过程进行仿真时，控制系统控制律参数多是以人工调整参数的方式进行计算得到，该种方案效率较低，且严重依靠经验，存在主观任意性，难以保证计算精度。
鉴于上述技术缺陷的存在提出本申请。
需注意的是，以上背景技术内容的公开仅用于辅助理解本发明的发明构思及技术方案，其并不必然属于本专利申请的现有技术，在没有明确的证据表明上述内容在本申请的申请日已经公开的情况下，上述背景技术不应当用于评价本申请的新颖性和创造性。
发明内容
本申请的目的是提供一种燃气轮机仿真动态性能控制律参数确定方法，以克服或减轻已知存在的至少一方面的技术缺陷。
本申请的技术方案是：
一种燃气轮机仿真动态性能控制律参数确定方法，包括：
基于设定的燃气轮机动态性能控制律参数种群，对应计算得到燃气轮机动态性能参数；
基于计算得到的燃气轮机动态性能参数，对应计算得到与实际燃气轮机动态性能参数间的适应度；
依据选择概率淘汰较差适应度对应的燃气轮机动态性能控制律参数种群，对剩余的燃气轮机动态性能控制律参数种群进行交叉、变异，产生燃气轮机动态性能控制律参数子种群，以子种群更新设定的燃气轮机动态性能控制律参数种群，直至达到遗传代数；
选取剩余燃气轮机动态性能控制律参数种群中的燃气轮机动态性能控制律参数。
根据本申请的至少一个实施例，上述的燃气轮机仿真动态性能控制律参数确定方法中，燃气轮机动态性能控制律参数种群中的燃气轮机动态性能控制律参数包括燃气轮机控制系统内控制回路的比例系数、积分系数、微分系数，以及包括燃气轮机控制系统外控制回路的比例系数、积分系数、微分系数。
根据本申请的至少一个实施例，上述的燃气轮机仿真动态性能控制律参数确定方法中，燃气轮机动态性能参数为燃气轮机的转速、总温、总压或功率。
根据本申请的至少一个实施例，上述的燃气轮机仿真动态性能控制律参数确定方法中，所述基于计算得到的燃气轮机动态性能参数，对应计算得到与实际燃气轮机动态性能参数间的适应度，具体为按照以下公式进行计算：
J＝∫t|e(t)|dt；
其中，
J为计算得到的燃气轮机动态性能参数、实际燃气轮机动态性能参
数间的适应度；
t为燃气轮机动态性能参数对应的时间；
e(t)为计算得到的燃气轮机动态性能参数、实际燃气轮机动态性能
参数在t时刻的偏差。
根据本申请的至少一个实施例，上述的燃气轮机仿真动态性能控制律参数确定方法中，所述选取剩余燃气轮机动态性能控制律参数种群中的燃气轮机动态性能控制律参数，具体为：
选取最佳适应度对应燃气轮机动态性能控制律参数种群中的燃气轮机动态性能控制律参数。
附图说明
图1是本申请实施例提供的示意图；
图2是本申请实施例提供的燃气轮机仿真动态性能控制律参数确定方法的流程图。
为了更好说明本实施例，附图某些部件会有省略、放大或缩小，并不代表实际产品的尺寸；此外，附图用于示例性说明，其中描述位置关系的用语仅限于示例性说明，不能理解为对本专利的限制。
具体实施方式
为使本申请的技术方案及其优点更加清楚，下面将结合附图对本申请的技术方案作进一步清楚、完整的详细描述，可以理解的是，此处所描述的具体实施例仅是本申请的部分实施例，其仅用于解释本申请，而非对本申请的限定。需要说明的是，为了便于描述，附图中仅示出了与本申请相关的部分，其他相关部分可参考通常设计，在不冲突的情况下，本申请中的实施例及实施例中的技术特征可以相互组合以得到新的实施例。
此外，除非另有定义，本申请描述中所使用的技术术语或者科学术语应当为本申请所属领域内一般技术人员所理解的通常含义。本申请描述中所使用的“上”、“下”、“左”、“右”、“中心”、“竖直”、“水平”、“内”、“外”等表示方位的词语仅用以表示相对的方向或者位置关系，而非暗示装置或元件必须具有特定的方位、以特定的方位构造和操作，当被描述对象的绝对位置发生改变后，其相对位置关系也可能发生相应的改变，因此不能理解为对本申请的限制。本申请描述中所使用的“第一”、“第二”、“第三”以及类似用语，仅用于描述目的，用以区分不同的组成部分，而不能够将其理解为指示或暗示相对重要性。本申请描述中所使用的“一个”、“一”或者“该”等类似词语，不应理解为对数量的绝对限制，而应理解为存在至少一个。本申请描述中所使用的“包括”或者“包含”等类似词语意指出现在该词前面的元件或者物件涵盖出现在该词后面列举的元件或者物件及其等同，而不排除其他元件或者物件。
此外，还需要说明的是，除非另有明确的规定和限定，在本申请的描述中使用的“安装”、“相连”、“连接”等类似词语应做广义理解，例如，连接可以是固定连接，也可以是可拆卸连接，或一体地连接；可以是机械连接，也可以是电连接；可以是直接相连，也可以通过中间媒介间接相连，还可以是两个元件内部的连通，领域内技术人员可根据具体情况理解其在本申请中的具体含义。
下面结合附图1至图2对本申请做进一步详细说明。
一种燃气轮机仿真动态性能控制律参数确定方法，包括：
基于设定的燃气轮机动态性能控制律参数种群，对应计算得到燃气轮机动态性能参数；
基于计算得到的燃气轮机动态性能参数，对应计算得到与实际燃气轮机动态性能参数间的适应度；
依据选择概率淘汰较差适应度对应的燃气轮机动态性能控制律参数种群，对剩余的燃气轮机动态性能控制律参数种群进行交叉、变异，产生燃气轮机动态性能控制律参数子种群，以子种群更新设定的燃气轮机动态性能控制律参数种群，直至达到遗传代数；
选取剩余燃气轮机动态性能控制律参数种群中的燃气轮机动态性能控制律参数。
对于上述实施例公开的燃气轮机仿真动态性能控制律参数确定方法，领域内技术人员可以理解的是，其基于遗传算法实施，通过多次迭代，快速搜索到与燃气轮机运行控制相符的燃气轮机动态性能控制律参数种群，选取其中的燃气轮机动态性能控制律参数对燃气轮机运行过程进行仿真，具有较高的精度。
对于上述实施例公开的燃气轮机仿真动态性能控制律参数确定方法，领域内技术人员还可以理解的是，其可基于MATLAB平台进行实施，在SIMULINK平台完成适应度函数定义，具有较高的效率。
对于上述实施例公开的燃气轮机仿真动态性能控制律参数确定方法，领域内技术人员还可以理解的是，设定的燃气轮机动态性能控制律参数种群的数量、选择概率的大小、交叉概率的大小、变异概率的大小、遗传代数的多少，可由相关技术人员在应用本申请时根据具体实际进行确定。
在一些可选的实施例中，上述的燃气轮机仿真动态性能控制律参数确定方法中，燃气轮机动态性能控制律参数种群中的燃气轮机动态性能控制律参数包括燃气轮机控制系统内控制回路的比例系数、积分系数、微分系数，以及包括燃气轮机控制系统外控制回路的比例系数、积分系数、微分系数。
在一些可选的实施例中，上述的燃气轮机仿真动态性能控制律参数确定方法中，燃气轮机动态性能参数为燃气轮机的转速、总温、总压或功率。
在一些可选的实施例中，上述的燃气轮机仿真动态性能控制律参数确定方法中，所述基于计算得到的燃气轮机动态性能参数，对应计算得到与实际燃气轮机动态性能参数间的适应度，具体为按照以下公式进行计算：
J＝∫t|e(t)|dt；
其中，
J为计算得到的燃气轮机动态性能参数、实际燃气轮机动态性能参
数间的适应度；
t为燃气轮机动态性能参数对应的时间；
e(t)为计算得到的燃气轮机动态性能参数、实际燃气轮机动态性能
参数在t时刻的偏差。
对于上述实施例公开的燃气轮机仿真动态性能控制律参数确定方法，领域内技术人员还可以理解的是，计算得到的燃气轮机动态性能参数、实际燃气轮机动态性能参数间的适应度J，表示以燃气轮机动态性能控制律参数种群中燃气轮机动态性能控制律参数对燃气轮机运行过程进行仿真，与实际燃气轮机运行动态性能参数的吻合程度，其数值越小表示吻合程度越高，以对应燃气轮机动态性能控制律参数种群中燃气轮机动态性能控制律参数对燃气轮机运行过程进行仿真，所得的结果有效度越高。
在一些可选的实施例中，上述的燃气轮机仿真动态性能控制律参数确定方法中，所述选取剩余燃气轮机动态性能控制律参数种群中的燃气轮机动态性能控制律参数，具体为：
选取最佳适应度对应燃气轮机动态性能控制律参数种群中的燃气轮机动态性能控制律参数。
说明书中各个实施例采用递进的方式描述，每个实施例重点说明的都是与其他实施例的不同之处，各个实施例之间相同相似部分互相参见即可。
至此，已经结合附图所示的优选实施方式描述了本申请的技术方案，领域内技术人员应该理解的是，本申请的保护范围显然不局限于这些具体实施方式，在不偏离本申请的原理的前提下，本领域技术人员可以对相关技术特征作出等同的更改或替换，这些更改或替换之后的技术方案都将落入本申请的保护范围之内。

Claims (5)

1.一种燃气轮机仿真动态性能控制律参数确定方法，其特征在于，包括：
基于设定的燃气轮机动态性能控制律参数种群，对应计算得到燃气轮机动态性能参数；
基于计算得到的燃气轮机动态性能参数，对应计算得到与实际燃气轮机动态性能参数间的适应度；
依据选择概率淘汰较差适应度对应的燃气轮机动态性能控制律参数种群，对剩余的燃气轮机动态性能控制律参数种群进行交叉、变异，产生燃气轮机动态性能控制律参数子种群，以子种群更新设定的燃气轮机动态性能控制律参数种群，直至达到遗传代数；
选取剩余燃气轮机动态性能控制律参数种群中的燃气轮机动态性能控制律参数。
2.根据权利要求1所述的燃气轮机仿真动态性能控制律参数确定方法，其特征在于，
燃气轮机动态性能控制律参数种群中的燃气轮机动态性能控制律参数包括燃气轮机控制系统内控制回路的比例系数、积分系数、微分系数，以及包括燃气轮机控制系统外控制回路的比例系数、积分系数、微分系数。
3.根据权利要求1所述的燃气轮机仿真动态性能控制律参数确定方法，其特征在于，
燃气轮机动态性能参数为燃气轮机的转速、总温、总压或功率。
4.根据权利要求1所述的燃气轮机仿真动态性能控制律参数确定方法，其特征在于，
所述基于计算得到的燃气轮机动态性能参数，对应计算得到与实际燃气轮机动态性能参数间的适应度，具体为按照以下公式进行计算：
J＝∫t|e(t)|dt；
其中，
J为计算得到的燃气轮机动态性能参数、实际燃气轮机动态性能参数间的适应度；
t为燃气轮机动态性能参数对应的时间；
e(t)为计算得到的燃气轮机动态性能参数、实际燃气轮机动态性能参数在t时刻的偏差。
5.根据权利要求1所述的燃气轮机仿真动态性能控制律参数确定方法，其特征在于，
所述选取剩余燃气轮机动态性能控制律参数种群中的燃气轮机动态性能控制律参数，具体为：
选取最佳适应度对应燃气轮机动态性能控制律参数种群中的燃气轮机动态性能控制律参数。
CN202110121594.1A 2021-01-28 2021-01-28 一种燃气轮机仿真动态性能控制律参数确定方法 Pending CN112947050A (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202110121594.1A CN112947050A (zh) 2021-01-28 2021-01-28 一种燃气轮机仿真动态性能控制律参数确定方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202110121594.1A CN112947050A (zh) 2021-01-28 2021-01-28 一种燃气轮机仿真动态性能控制律参数确定方法

Publications (1)

Publication Number Publication Date
CN112947050A true CN112947050A (zh) 2021-06-11

Family

ID=76239044

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202110121594.1A Pending CN112947050A (zh) 2021-01-28 2021-01-28 一种燃气轮机仿真动态性能控制律参数确定方法

Country Status (1)

Country Link
CN (1) CN112947050A (zh)

Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2002323904A (ja) * 2001-02-01 2002-11-08 Honda Motor Co Ltd プラントの制御装置及び内燃機関の空燃比制御装置
US20150192912A1 (en) * 2014-01-07 2015-07-09 General Electric Company Method and system for combustion mode transfer for a gas turbine engine
CN108573300A (zh) * 2018-04-18 2018-09-25 辽东学院 粒子群算法反馈自适应率参数优化方法及全阶状态观测器
CN109802609A (zh) * 2019-01-01 2019-05-24 武汉船用电力推进装置研究所（中国船舶重工集团公司第七一二研究所） 一种异步电机无速度传感器调速系统pi参数整定方法
CN109918693A (zh) * 2018-11-09 2019-06-21 西南交通大学 一种异步联网后水轮机调速器参数的分析和整定方法
CN111443599A (zh) * 2020-04-15 2020-07-24 上海电力大学 优化pid参数的核电站汽轮机转速控制方法
CN111488690A (zh) * 2020-04-16 2020-08-04 云南电网有限责任公司电力科学研究院 一种冲击式水轮发电机调速器空载调节参数优化方法

Patent Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2002323904A (ja) * 2001-02-01 2002-11-08 Honda Motor Co Ltd プラントの制御装置及び内燃機関の空燃比制御装置
US20150192912A1 (en) * 2014-01-07 2015-07-09 General Electric Company Method and system for combustion mode transfer for a gas turbine engine
CN108573300A (zh) * 2018-04-18 2018-09-25 辽东学院 粒子群算法反馈自适应率参数优化方法及全阶状态观测器
CN109918693A (zh) * 2018-11-09 2019-06-21 西南交通大学 一种异步联网后水轮机调速器参数的分析和整定方法
CN109802609A (zh) * 2019-01-01 2019-05-24 武汉船用电力推进装置研究所（中国船舶重工集团公司第七一二研究所） 一种异步电机无速度传感器调速系统pi参数整定方法
CN111443599A (zh) * 2020-04-15 2020-07-24 上海电力大学 优化pid参数的核电站汽轮机转速控制方法
CN111488690A (zh) * 2020-04-16 2020-08-04 云南电网有限责任公司电力科学研究院 一种冲击式水轮发电机调速器空载调节参数优化方法

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
李洁: "基于交叉熵算法的PID 控制器设计", 控制与决策, vol. 26, no. 5, pages 794 - 796 *
王顺 等: "基于差分进化算法的燃气轮机转速模糊PID 复合控制", 热能动力工程, vol. 32, no. 10, pages 40 - 44 *
王顺;智同生;施建强;: "自适应差分进化算法在燃气轮机转速控制中的应用", 热能动力工程, vol. 33, no. 10, pages 28 - 33 *

Similar Documents

Publication Publication Date Title
CN103970965B (zh) 燃气涡轮发动机加速寿命试验试车方法
CN105930662A (zh) 一种压气机低转速特性外推方法
CN106870955A (zh) 服务于供水管网节点需水量反演的管网监测点优化布置方法
CN103984813A (zh) 一种离心压缩机裂纹叶轮结构的振动建模与分析方法
US20130282307A1 (en) Turbo-machinery stage families tuning/calibration system and method
CN107273565A (zh) 一种基于les/lbm的离心泵内流动诱导噪声数值计算方法
CN111027148B (zh) 损失落后角模型自动标定及工业轴流压缩机性能计算方法
CN108363902A (zh) 一种致病遗传变异的精确预测方法
CN104134013A (zh) 一种风力机叶片模态分析方法
CN106295015B (zh) 一种渐开线直齿圆柱齿轮副的齿廓修形方法及与其配套的专用参数化cad系统
JP2017215832A (ja) 回転機械の運転状態を診断する診断装置及び診断方法
CN103995938A (zh) 一种压气机中介机匣端壁型线造型方法
CN105508147B (zh) 风电叶片单点疲劳加载试验弯矩匹配方法
CN112966399B (zh) 一种基于机器学习的脉管制冷机工况预测方法及系统
CN112947050A (zh) 一种燃气轮机仿真动态性能控制律参数确定方法
Schmitz et al. Novel performance prediction of a transonic 4.5 stage compressor
CN112257204B (zh) 多级压气机s2流面参数的计算方法
CN105956318A (zh) 基于改进分裂 h-k 聚类方法的风电场机群划分方法
CN109063313B (zh) 基于机器学习的列车牵引能耗计算方法
CN112949166A (zh) 一种燃气轮机稳态模型系数确定方法
CN115758603A (zh) 航空发动机涡轮叶片多工况蠕变/应力断裂寿命预测方法
CN114841091A (zh) 一种基于代理模型技术的离心通风机的气动优化设计方法
CN110955999B (zh) 轮胎三维温度场仿真方法
CN115329415A (zh) 一种大跨越输电塔的风振系数确定方法及系统
CN102521484B (zh) 基于自适应参数域遗传算法的机油泵性能曲线的估计方法

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