CN106055792A - 一种超临界二氧化碳压气机与涡轮匹配动态仿真方法 - Google Patents

一种超临界二氧化碳压气机与涡轮匹配动态仿真方法 Download PDF

Info

Publication number
CN106055792A
CN106055792A CN201610378626.5A CN201610378626A CN106055792A CN 106055792 A CN106055792 A CN 106055792A CN 201610378626 A CN201610378626 A CN 201610378626A CN 106055792 A CN106055792 A CN 106055792A
Authority
CN
China
Prior art keywords
turbine
compressor
module
dynamic
model
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
CN201610378626.5A
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.)
Harbin Engineering University
Original Assignee
Harbin Engineering University
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 Harbin Engineering University filed Critical Harbin Engineering University
Priority to CN201610378626.5A priority Critical patent/CN106055792A/zh
Publication of CN106055792A publication Critical patent/CN106055792A/zh
Pending legal-status Critical Current

Links

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00Computer-aided design [CAD]
    • G06F30/30Circuit design
    • G06F30/36Circuit design at the analogue level
    • G06F30/367Design verification, e.g. using simulation, simulation program with integrated circuit emphasis [SPICE], direct methods or relaxation methods

Landscapes

  • Engineering & Computer Science (AREA)
  • Computer Hardware Design (AREA)
  • Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Microelectronics & Electronic Packaging (AREA)
  • Evolutionary Computation (AREA)
  • Geometry (AREA)
  • General Engineering & Computer Science (AREA)
  • General Physics & Mathematics (AREA)
  • Management, Administration, Business Operations System, And Electronic Commerce (AREA)

Abstract

本发明涉及一种超临界二氧化碳压气机与涡轮匹配动态仿真方法。本发明包括:步骤A):绘制压气机与涡轮共同工作网;采用加热器功率动态加减速系数法建立压气机与涡轮动态仿真模型。本发明用于解决缺乏部件特性数据、针对按照理想气体设计的传统的压气机与涡轮动态实时模型机载实时模型建模问题,仿真精度更高,稳定性更好。

Description

一种超临界二氧化碳压气机与涡轮匹配动态仿真方法
技术领域
本发明涉及一种超临界二氧化碳压气机与涡轮匹配动态仿真方法。
背景技术
在基于超临界二氧化碳布雷顿循环的发电系统中,压气机与涡轮是最整个装置的核心部件,它们都有自己的工作特性,当把它们在组装成一个整体时,可能会因为互相制约从而偏离了自己的合理的运行区域,继而发生故障,甚至不能正常工作,所以它们匹配的好坏程度直接影响整个装置的效率。
国内外对于超临界二氧化碳作为介质的压气机与涡轮动态仿真模型由于缺乏部件特性数据,所以仿真模型都是根据理想气体设计的压气机和涡轮而建立的,而二氧化碳在超临界状态下不能看成是理想气体,而是实际气体,若按理想气体设计,在仿真实验中会出现运行不稳定,出现喘振等现象,不能满足机载需求。因此,需要研究根据实际气体而建立的压气机和涡轮动态实时模型的方法。
发明内容
本发明的目的在于提供一种用于解决缺乏部件特性数据、针对按照理想气体设计的传统的压气机与涡轮动态实时模型机载实时模型建模问题的超临界二氧化碳压气机与涡轮匹配动态仿真方法。
本发明的目的是这样实现的:
包括如下步骤:
步骤A):绘制压气机与涡轮共同工作网;
步骤Al):根据压气机和涡轮的特性曲线,计算压气机和涡轮的共同工作点;
步骤A2):主要对压气机和涡轮部件特性进行数学建模,通过计算压气机与涡轮的共同工作点,得到装置的变工况特性网;
步骤B):采用加热器功率动态加减速系数法建立压气机与涡轮动态仿真模型;
步骤B1),采用模块化的非线性建模方法来建立超临界CO2装置各部件的数学模型,模型中主要包含的模块有压气机与涡轮模块,加热器模块与换热器模块,容积模块,转子与负载模块,以及机械损失模块等。因为此次仿真重点在于观察压气机与涡轮的匹配特性;
步骤B2),利用之前搭建好的数学模型,在SIMULINK软件中分别建立各相应模块的仿真模型。主要采用变步长连续算法进行控制方程组的求解。仿真模型包括压气机仿真模型、涡轮仿真模型、加热器仿真模型、转子模块与机械损失模块仿真模;
步骤B3),采用控制加热器功率的方法,主要通过控制加热器功率,从而使涡轮入口温度改变,进而调节流量、膨胀比、输出功。转速由功率变化率决定,作为反馈信号,在用于电机发电时一般为定值。模拟超临界CO2装置沿等温比线变化时,各个参量的变化情况,从而得出部件时间的动态匹配性能。
本发明的有益效果在于:
本发明用于解决缺乏部件特性数据、针对按照理想气体设计的传统的压气机与涡轮动态实时模型机载实时模型建模问题,仿真精度更高,稳定性更好。
附图说明
图1为S-CO2装置单轴结构示意图;
图2为沿着不同温比运行时压气机与涡轮共同工作网;
图3为变工况计算流程图;
图4为整体装置仿真模型;
图5为压气机仿真模型;
图6涡轮模块仿真模型;
图7为加热器仿真模型;
图8为转子模块仿真模型;
图9为机械损失模块仿真模型。
具体实施方式
下面结合附图对本发明做进一步描述。
该方法和系统满足机载实时性要求,能为仿真实验运行稳定做保证。
本发明为解决上述技术问题采用以下技术方案
一种超临界二氧化碳压气机与涡轮匹配动态仿真模型建模方法及系统的制作方法,包括以下步骤:
步骤A),绘制压气机与涡轮共同工作网;
步骤B),采用加热器功率动态加减速系数法建立压气机与涡轮动态仿真模型。
作为本发明超临界二氧化碳压气机与涡轮匹配动态仿真模型建模方法进一步的优化方案,步骤A)的具体步骤如下:
步骤Al),根据压气机和涡轮的特性曲线,计算压气机和涡轮的共同工作点;
步骤A2),主要对压气机和涡轮部件特性进行数学建模,通过计算压气机与涡轮的共同工作点,得到装置的变工况特性网;
作为本发明一种超临界二氧化碳压气机与涡轮匹配动态仿真模型建模方法进一步的优化方案,步骤A2)主要对压气机和涡轮部件特性进行数学建模,通过计算压气机与涡轮的共同工作点,得到装置的变工况特性网的详细步骤如下:惯性微分方程与装置平衡关系方程,是本装置动态与稳态计算的数学模型。如计算平衡运行工况时(包括设计工况与各种部分负荷工况),将惯性微分方程方程转换为代数方程即可。而稳态特性的求解归结为非线性代数方程组的求解,在求解时把辅助方程方程带入,结合给定的参变量就可求出方程组的唯一解,即表征一个特定的运行工况。如另选一组参变量又可求出另外一组解,表征另一个特定工况。而参变量的个数取决于可调参数的多少,在变工况时,控制不同的输入参数便可得到对应工况下的涡轮与压气机的匹配工作点。
作为本发明一种超临界二氧化碳压气机与涡轮匹配动态仿真模型建模方法进一步的优化方案,步骤B)的具体步骤如下:
步骤B1),采用模块化的非线性建模方法来建立超临界CO2装置各部件的数学模型,模型中主要包含的模块有压气机与涡轮模块,加热器模块与换热器模块,容积模块,转子与负载模块,以及机械损失模块等。因为此次仿真重点在于观察压气机与涡轮的匹配特性。
步骤B2),利用之前搭建好的数学模型,在SIMULINK软件中分别建立各相应模块的仿真模型。主要采用变步长连续算法进行控制方程组的求解。仿真模型包括压气机仿真模型、涡轮仿真模型、加热器仿真模型、转子模块与机械损失模块仿真模型
步骤B3),采用控制加热器功率的方法,主要通过控制加热器功率,从而使涡轮入口温度改变,进而调节流量、膨胀比、输出功。转速由功率变化率决定,作为反馈信号,在用于电机发电时一般为定值。模拟超临界CO2装置沿等温比线变化时,各个参量的变化情况,从而得出部件时间的动态匹配性能。
作为本发明一种超临界二氧化碳压气机与涡轮匹配动态仿真模型建模方法进一步的优化方案,所述的超临界CO2装置整体结构选择单轴结构进行设计。
作为本发明一种超临界二氧化碳压气机与涡轮匹配动态仿真模型建模方法进一步的优化方案,步骤A1)中根据压气机和涡轮的特性曲线,计算压气机和涡轮的共同工作点时,超临界CO2装置进行动态调节时,采用节流阀调节压比的分配和流量的控制;超临界CO2装置也可以看做是一个惯性系统;计算不同温比下各个转速下的匹配点,绘制于压气机特性图中,便可得到整个装置的变工况共同工作网。
作为本发明一种超临界二氧化碳压气机与涡轮匹配动态仿真模型建模方法进一步的优化方案,步骤B1)中压气机和涡轮的工作特性可以用压比π、折合转速和折合流量以及效率η四个参数的关系来表示
图1是超临界二氧化碳装置单轴结构图,该动态仿真模型及系统的建立方法包括以下步骤:
步骤A),绘制压气机与涡轮共同工作网(图2所示);
步骤B),采用加热器功率动态加减速系数法建立压气机与涡轮动态仿真模型。
步骤Al),根据压气机和涡轮的特性曲线,计算压气机和涡轮的共同工作点;将计算所得的数据整理成以流量、转速为自变量,膨胀比、效率等为参变量的形式,并绘制于S-CO2装置涡轮特性图上。
步骤A2),主要对压气机和涡轮部件特性进行数学建模,通过计算压气机与涡轮的共同工作点,得到装置的变工况特性网;惯性微分方程与装置平衡关系方程,是本装置动态与稳态计算的数学模型。如计算平衡运行工况时(包括设计工况与各种部分负荷工况),将惯性微分方程方程转换为代数方程即可。而稳态特性的求解归结为非线性代数方程组的求解,在求解时把辅助方程方程带入,结合给定的参变量就可求出方程组的唯一解,即表征一个特定的运行工况。如另选一组参变量又可求出另外一组解,表征另一个特定工况。而参变量的个数取决于可调参数的多少,在变工况时,控制不同的输入参数便可得到对应工况下的涡轮与压气机的匹配工作点。在计算时,为了方便程序的自动化运行,可以先给定T3 *,然后假设点A进行计算,当结果大于设定误差时,返回用等转速线另外一点A1进行迭代(如图3所示),如此就可在程序中实现等温比线的绘制。
作为本发明一种超临界二氧化碳压气机与涡轮匹配动态仿真模型建模方法进一步的优化方案,步骤B)的具体步骤如下:
步骤B1),采用模块化的非线性建模方法来建立超临界CO2装置各部件的数学模型,模型中主要包含的模块有压气机与涡轮模块,加热器模块与换热器模块,容积模块,转子与负载模块,以及机械损失模块等。因为此次仿真重点在于观察压气机与涡轮的匹配特性。压气机的工作特性可以用压比π、折合转速和折合流量以及效率η四个参数的关系来表示,涡轮模块数学模型可以仿照压气机的方式建模,用膨胀比和折合转速为自变量,折合流量和效率为输出量建立相应的数学模型;加热器模块数学模型采用电加热,根据能量和质量守恒关系建立模型;负载模块数学模型采用恒速负载。
步骤B2),利用之前搭建好的数学模型,在SIMULINK软件中分别建立各相应模块的仿真模型(图4)。主要采用变步长连续算法,具体算法为龙格库塔法进行控制方程组的求解。仿真模型包括压气机仿真模型(图5)、涡轮仿真模型(图6)、加热器仿真模型(图7)、转子模块(图8)与机械损失模块仿真模型(图9)。
步骤B3),采用控制加热器功率的方法,主要通过控制加热器功率,从而使涡轮入口温度改变,进而调节流量、膨胀比、输出功。转速由功率变化率决定,作为反馈信号,在用于电机发电时一般为定值。模拟超临界CO2装置沿等温比线变化时,各个参量的变化情况,从而得出部件时间的动态匹配性能。

Claims (1)

1.一种超临界二氧化碳压气机与涡轮匹配动态仿真方法,其特征在于,包括如下步骤:
步骤A):绘制压气机与涡轮共同工作网;
步骤Al):根据压气机和涡轮的特性曲线,计算压气机和涡轮的共同工作点;
步骤A2):主要对压气机和涡轮部件特性进行数学建模,通过计算压气机与涡轮的共同工作点,得到装置的变工况特性网;
步骤B):采用加热器功率动态加减速系数法建立压气机与涡轮动态仿真模型;
步骤B1),采用模块化的非线性建模方法来建立超临界CO2装置各部件的数学模型,模型中主要包含的模块有压气机与涡轮模块,加热器模块与换热器模块,容积模块,转子与负载模块,以及机械损失模块等因为此次仿真重点在于观察压气机与涡轮的匹配特性;
步骤B2),利用之前搭建好的数学模型,在SIMULINK软件中分别建立各相应模块的仿真模型主要采用变步长连续算法进行控制方程组的求解仿真模型包括压气机仿真模型、涡轮仿真模型、加热器仿真模型、转子模块与机械损失模块仿真模;
步骤B3),采用控制加热器功率的方法,主要通过控制加热器功率,从而使涡轮入口温度改变,进而调节流量、膨胀比、输出功转速由功率变化率决定,作为反馈信号,在用于电机发电时一般为定值模拟超临界CO2装置沿等温比线变化时,各个参量的变化情况,从而得出部件时间的动态匹配性能。
CN201610378626.5A 2016-05-31 2016-05-31 一种超临界二氧化碳压气机与涡轮匹配动态仿真方法 Pending CN106055792A (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201610378626.5A CN106055792A (zh) 2016-05-31 2016-05-31 一种超临界二氧化碳压气机与涡轮匹配动态仿真方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201610378626.5A CN106055792A (zh) 2016-05-31 2016-05-31 一种超临界二氧化碳压气机与涡轮匹配动态仿真方法

Publications (1)

Publication Number Publication Date
CN106055792A true CN106055792A (zh) 2016-10-26

Family

ID=57171972

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201610378626.5A Pending CN106055792A (zh) 2016-05-31 2016-05-31 一种超临界二氧化碳压气机与涡轮匹配动态仿真方法

Country Status (1)

Country Link
CN (1) CN106055792A (zh)

Cited By (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN111832152A (zh) * 2020-05-27 2020-10-27 东南大学 一种深冷高压储氢容器泄放过程仿真建模方法
CN111859563A (zh) * 2020-07-10 2020-10-30 西安交通大学 一种超临界二氧化碳透平试验的相似模化方法
CN113868982A (zh) * 2021-10-21 2021-12-31 山东大学 超临界二氧化碳径流式涡轮机械的数值仿真方法及系统
CN116502568A (zh) * 2023-06-28 2023-07-28 中国人民解放军国防科技大学 压气机内流特性自动化模拟的方法、装置、设备和介质

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101893515A (zh) * 2010-07-21 2010-11-24 奇瑞汽车股份有限公司 一种对发动机与涡轮增压器进行匹配的方法及装置
CN101949324A (zh) * 2010-08-17 2011-01-19 清华大学 一种内燃机涡轮增压通流匹配方法
CN102589894A (zh) * 2012-03-01 2012-07-18 南京航空航天大学 微型压气机/涡轮联合试验台及试验方法
CN105404750A (zh) * 2015-12-07 2016-03-16 南京航空航天大学 一种基于改进Broyden算法求解方程组的涡轴发动机自适应模型建立方法

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101893515A (zh) * 2010-07-21 2010-11-24 奇瑞汽车股份有限公司 一种对发动机与涡轮增压器进行匹配的方法及装置
CN101949324A (zh) * 2010-08-17 2011-01-19 清华大学 一种内燃机涡轮增压通流匹配方法
CN102589894A (zh) * 2012-03-01 2012-07-18 南京航空航天大学 微型压气机/涡轮联合试验台及试验方法
CN105404750A (zh) * 2015-12-07 2016-03-16 南京航空航天大学 一种基于改进Broyden算法求解方程组的涡轴发动机自适应模型建立方法

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
LIAO JIXIANG等: "Thermodynamic analysis of low-temperature power generation transcritical Rankine cycle with R41 and CO2", 《PROCEEDINGS OF ASME TURBO EXPO 2014: TURBINE TECHNICAL CONFERENCE AND EXPOSITION》 *
刘洋等: "固体燃料ATR涡轮/压气机匹配方法研究", 《推进技术》 *
梁墩煌等: "核反应堆系统中以超临界二氧化碳为工质的热力循环过程的建模与分析", 《厦门大学学报(自然科学版)》 *

Cited By (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN111832152A (zh) * 2020-05-27 2020-10-27 东南大学 一种深冷高压储氢容器泄放过程仿真建模方法
CN111859563A (zh) * 2020-07-10 2020-10-30 西安交通大学 一种超临界二氧化碳透平试验的相似模化方法
CN111859563B (zh) * 2020-07-10 2023-04-28 西安交通大学 一种超临界二氧化碳透平试验的相似模化方法
CN113868982A (zh) * 2021-10-21 2021-12-31 山东大学 超临界二氧化碳径流式涡轮机械的数值仿真方法及系统
CN116502568A (zh) * 2023-06-28 2023-07-28 中国人民解放军国防科技大学 压气机内流特性自动化模拟的方法、装置、设备和介质
CN116502568B (zh) * 2023-06-28 2023-09-05 中国人民解放军国防科技大学 压气机内流特性自动化模拟的方法、装置、设备和介质

Similar Documents

Publication Publication Date Title
CN106874569B (zh) 一种变几何分轴式燃气轮机实时仿真建模方法
CN101446807B (zh) 一种电力系统仿真中火电厂调速系统模型的实现方法
CN104734147B (zh) 一种综合能源系统概率能量流分析方法
CN103838920B (zh) 一种涡轴发动机动态实时模型建模方法及系统
CN106055792A (zh) 一种超临界二氧化碳压气机与涡轮匹配动态仿真方法
CN110829487B (zh) 一种电力系统的频率动态预测方法
CN107257132B (zh) 一种考虑风电弃用的电-气互联系统综合负荷削减模型构建方法
Calero et al. Compressed air energy storage system modeling for power system studies
CN101609312B (zh) 采用汽耗率修正滑压运行曲线的发电机组优化控制方法
CN107482649B (zh) 一种基于分频控制的两域互联系统负荷频率控制方法
CN105162164B (zh) 一种建立含风电并网系统的低阶动态频率响应模型的方法
CN110750882A (zh) 一种考虑频率约束的风电占比极限值解析计算方法
CN108223139B (zh) 一种分轴式燃气轮机动力涡轮前放气调节规律优化方法
CN105186511B (zh) 电池储能系统参与电网二次调频控制方法
CN111400937B (zh) 一种综合能源系统潮流计算方法
CN108875269A (zh) 一种考虑电力系统多平衡机与天然气系统慢动态特性的电-气互联系统多时段能流计算方法
CN105720592A (zh) 一种聚合空调负荷参与系统辅助二次调频的方法
CN107301268A (zh) 一种船舶燃气轮机可调静叶压气机偏转角度优化方法
CN107346889B (zh) 考虑一二次调频及最小频率偏差的负荷削减优化模型构建方法
Yee et al. Validated models for gas turbines based on thermodynamic relationships
CN110689199B (zh) 一种基于混合潮流模型的分布式协同调控方法
CN110350554A (zh) 基于混联结构的风储系统辅助电网一次调频控制方法
CN109858129A (zh) 一种关于三联供系统的燃气轮机动态仿真方法
CN112838580B (zh) 一种改进异质温控负荷聚合模型及多目标协调控制方法
CN114079293B (zh) 一种计及温控负荷与微型燃气轮机的电-气互联系统新能源功率波动平滑方法

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
RJ01 Rejection of invention patent application after publication

Application publication date: 20161026

RJ01 Rejection of invention patent application after publication