CN104008221A - 一种汽轮机模型的建立方法 - Google Patents
一种汽轮机模型的建立方法 Download PDFInfo
- Publication number
- CN104008221A CN104008221A CN201310390705.4A CN201310390705A CN104008221A CN 104008221 A CN104008221 A CN 104008221A CN 201310390705 A CN201310390705 A CN 201310390705A CN 104008221 A CN104008221 A CN 104008221A
- Authority
- CN
- China
- Prior art keywords
- model
- represent
- steam turbine
- regenerator
- servomotor
- 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
Links
Landscapes
- Control Of Turbines (AREA)
Abstract
一种汽轮机模型的建立方法,它涉及一种模型的建立方法,具体涉及一种汽轮机模型的建立方法。本发明为了解决解决现有汽轮机模型建立方法计算复杂,且建立的模型仿真程度不高,容易给后续生产带来损失和误差的问题。本发明的具体步骤为:建立油动机模型;建立蒸汽容积效应模型;建立转子模型;建立回热系统模型;建立再热系统模型;建立无穷大电网模型;结合油动机模型、蒸汽容积效应模型、转子模型、回热系统模型、再热系统模型、无穷大电网模型建立汽轮机模型。本发明用于汽轮机设计制造领域。
Description
技术领域
本发明涉及一种模型的建立方法,具体涉及一种汽轮机模型的建立方法。
背景技术
汽轮机是发电厂的重要设备,汽轮机在设计制造前都需要对其建立仿真模型,以此来确定最终的制造方案和定型生产规划。现有汽轮机模型建立方法计算复杂,且建立的模型仿真程度不高,容易给后续生产带来损失和误差。
发明内容
本发明为解决现有汽轮机模型建立方法计算复杂,且建立的模型仿真程度不高,容易给后续生产带来损失和误差的问题,进而提出一种汽轮机模型的建立方法。
本发明为解决上述问题采取的技术方案是:本发明的具体步骤如下:
步骤一、建立油动机模型:利用Matlab建模软件建立油动机仿真模型,其中油动机的惯性时间常数为0.084s时,对应的阀门关闭时间为0.10s,油动机的惯性时间常数为0.125s时,对应的阀门关闭时间为0.15s,油动机的惯性时间常数为0.167s时,对应的阀门关闭时间为0.20s,表征摩擦损失的自平衡系数Beita为0.03;
步骤二、建立蒸汽容积效应模型:根据公式建立蒸汽容积效应仿真模型,其中T0表示时间容积常数,l表示进气管道长度,A表示进气管道横截面积,v0表示气流在管道中的流速,ρ0表示额定工况下的气体密度;
步骤三、建立转子模型:根据公式χP(s)-χNL(s)=(β+Ta·s)χn(s)建立转子仿真模型,其中ω表示转子的角速度,ω0表示转子在额定工况下的角速度,NT0表示额定透平功率,NL表示决定于外部因素的量,NT表示透平功率表示转子的飞升时间常数,β表示自平衡系数;
步骤四、建立回热系统模型:
步骤四(一)、建立回热器的动态数学模型;
步骤四(二)、建立计及回热器热效应的汽轮机动态模型;
步骤四(三)、根据步骤四(一)和步骤四(二)建立回热系统模型;
步骤五、建立再热系统模型:
步骤六、建立无穷大电网模型;
步骤七、结合油动机模型、蒸汽容积效应模型、转子模型、回热系统模型、再热系统模型、无穷大电网模型建立汽轮机模型。
本发明的有益效果是:本发明计算简单,可迅速建立所需的汽轮机仿真模型,利用本发明所述方法建立的汽轮机仿真模型,精确度高,不会给后续设计和生产造成损失和误差。
附图说明
图1是本发明的流程框图。
具体实施方式
具体实施方式一:结合图1说明本实施方式,本实施方式所述一种汽轮机模型的建立方法的具体步骤如下:
步骤一、建立油动机模型:利用Matlab建模软件建立油动机仿真模型,其中油动机的惯性时间常数为0.084s时,对应的阀门关闭时间为0.10s,油动机的惯性时间常数为0.125s时,对应的阀门关闭时间为0.15s,油动机的惯性时间常数为0.167s时,对应的阀门关闭时间为0.20s,表征摩擦损失的自平衡系数Beita为0.03;
步骤二、建立蒸汽容积效应模型:根据公式建立蒸汽容积效应仿真模型,其中T0表示时间容积常数,l表示进气管道长度,A表示进气管道横截面积,v0表示气流在管道中的流速,ρ0表示额定工况下的气体密度;
步骤三、建立转子模型:根据公式χP(s)-χNL(s)=(β+Ta·s)χn(s)建立转子仿真模型,其中ω表示转子的角速度,ω0表示转子在额定工况下的角速度,NT0表示额定透平功率,NL表示决定于外部因素的量,NT表示透平功率Ta表示转子的飞升时间常数,β表示自平衡系数;
步骤四、建立回热系统模型:
步骤四(一)、建立回热器的动态数学模型;
步骤四(二)、建立计及回热器热效应的汽轮机动态模型;
步骤四(三)、根据步骤四(一)和步骤四(二)建立回热系统模型;
步骤五、建立再热系统模型:
步骤六、建立无穷大电网模型;
步骤七、结合油动机模型、蒸汽容积效应模型、转子模型、回热系统模型、再热系统模型、无穷大电网模型建立汽轮机模型。
具体实施方式二:结合图1说明本实施方式,本实施方式所述一种汽轮机模型的建立方法的步骤二中时间容积常数T0=0.1~0.3s。其它组成及连接关系与具体实施方式一相同。
具体实施方式三:结合图1说明本实施方式,本实施方式所述一种汽轮机模型的建立方法的步骤三中β的取值范围为0.03~0.05。其它组成及连接关系与具体实施方式一相同。
具体实施方式四:结合图1说明本实施方式,本实施方式所述一种汽轮机模型的建立方法的步骤四(二)中建立计及回热器热效应的汽轮机动态模型的具体步骤如下:
步骤A、确定静态参数K1、K2、K3:
公式①②③中αri表示第i级回热器抽汽份额系数,Dms0表示额定工况下上级组流量, 表示额定工况下下级组流量,Pri表示第i级回热器压力,Pbi表示第i级回热器抽汽点压力,βi表示回热器逆止阀门前后的压力损失系数,Pri0表示额定工况下第i级回热器压力,Pbi0表示额定工况下第i级回热器抽汽点压力,表示Pri点标幺量,表示Pbi点标幺量;
步骤B、确定回热器动态特性传递函数Wi(s):
公式④和⑤中V表示含抽汽管在内的加热器容积,Tr表示回热器蓄热时间常数对于混合式加热器Tr=30s,对于表面式加热器Tr=20s,ρ0表示额定工况下气体密度,Dr0表示额定流量,Cps表示汽侧比热系数,Cpw表示给水侧比热系数,Ms表示汽侧质量,Dw表示给水侧流量,n表示多变指数;
步骤C、将K1,K2,K3带入信号流,写成传递函数的形式并化简可得给水回热加热系统的数学模型。
其它组成及连接关系与具体实施方式一相同。
具体实施方式五:结合图1说明本实施方式,本实施方式所述一种汽轮机模型的建立方法的步骤五中再热容积时间常数为8~10秒。其它组成及连接关系与具体实施方式一相同。
Claims (5)
1.一种汽轮机模型的建立方法,其特征在于:所述一种汽轮机模型的建立方法的具体步骤如下:
步骤一、建立油动机模型:利用Matlab建模软件建立油动机仿真模型,其中油动机的惯性时间常数为0.084s时,对应的阀门关闭时间为0.10s,油动机的惯性时间常数为0.125s时,对应的阀门关闭时间为0.15s,油动机的惯性时间常数为0.167s时,对应的阀门关闭时间为0.20s,表征摩擦损失的自平衡系数Beita为0.03;
步骤二、建立蒸汽容积效应模型:根据公式建立蒸汽容积效应仿真模型,其中T0表示时间容积常数,l表示进气管道长度,A表示进气管道横截面积,v0表示气流在管道中的流速,ρ0表示额定工况下的气体密度;
步骤三、建立转子模型:根据公式χP(s)-χNL(s)=(β+Ta·s)χn(s)建立转子仿真模型,其中ω表示转子的角速度,ω0表示转子在额定工况下的角速度,NT0表示额定透平功率,NL表示决定于外部因素的量,NT表示透平功率Ta表示转子的飞升时间常数,β表示自平衡系数;
步骤四、建立回热系统模型:
步骤四(一)、建立回热器的动态数学模型;
步骤四(二)、建立计及回热器热效应的汽轮机动态模型;
步骤四(三)、根据步骤四(一)和步骤四(二)建立回热系统模型;
步骤五、建立再热系统模型:
步骤六、建立无穷大电网模型;
步骤七、结合油动机模型、蒸汽容积效应模型、转子模型、回热系统模型、再热系统模型、无穷大电网模型建立汽轮机模型。
2.根据权利要求1所述一种汽轮机模型的建立方法,其特征在于:步骤二中时间容积常数T0=0.1~0.3s。
3.根据权利要求1所述一种汽轮机模型的建立方法,其特征在于:步骤三中β的取值范围为0.03~0.05。
4.根据权利要求1所述一种汽轮机模型的建立方法,其特征在于:步骤四(二)中建立计及回热器热效应的汽轮机动态模型的具体步骤如下:
步骤A、确定静态参数K1、K2、K3:
①,
②,
③,
公式①②③中αri表示第i级回热器抽汽份额系数,Dms0表示额定工况下上级组流量,D′ms0表示额定工况下下级组流量,Pri表示第i级回热器压力,Pbi表示第i级回热器抽汽点压力,βi表示回热器逆止阀门前后的压力损失系数,Pri0表示额定工况下第i级回热器压力,Pbi0表示额定工况下第i级回热器抽汽点压力,表示Pri点标幺量,表示Pbi点标幺量;
步骤B、确定回热器动态特性传递函数Wi(s):
④,
⑤,
公式④和⑤中V表示含抽汽管在内的加热器容积,Tr表示回热器蓄热时间常数对于混合式加热器Tr=30s,对于表面式加热器Tr=20s,ρ0表示额定工况下气体密度,Dr0表示额定流量,Cps表示汽侧比热系数,Cpw表示给水侧比热系数,Ms表示汽侧质量,Dw表示给水侧流量,n表示多变指数;
步骤C、将K1,K2,K3带入信号流,写成传递函数的形式并化简可得给水回热加热 系统的数学模型。
5.根据权利要求1所述一种汽轮机模型的建立方法,其特征在于:步骤五中再热容积时间常数为8~10秒。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201310390705.4A CN104008221B (zh) | 2013-09-02 | 2013-09-02 | 一种汽轮机模型的建立方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201310390705.4A CN104008221B (zh) | 2013-09-02 | 2013-09-02 | 一种汽轮机模型的建立方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN104008221A true CN104008221A (zh) | 2014-08-27 |
CN104008221B CN104008221B (zh) | 2018-07-31 |
Family
ID=51368878
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201310390705.4A Expired - Fee Related CN104008221B (zh) | 2013-09-02 | 2013-09-02 | 一种汽轮机模型的建立方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN104008221B (zh) |
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104376157A (zh) * | 2014-11-03 | 2015-02-25 | 哈尔滨汽轮机厂有限责任公司 | 大型汽轮发电机组轴系动静态特性分析计算系统 |
CN105260548A (zh) * | 2015-10-21 | 2016-01-20 | 国网湖南节能服务有限公司 | 一种基于机组实际特性的汽轮机模型建模方法 |
CN116950962A (zh) * | 2023-09-21 | 2023-10-27 | 福建福清核电有限公司 | 油动机带载运行调节的方法及其装置、电子设备 |
Family Cites Families (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102073272B (zh) * | 2011-02-24 | 2013-05-15 | 北京恩耐特分布能源技术有限公司 | 一种分布式能源整体规划优化的系统及其优化方法 |
CN102663522B (zh) * | 2012-04-24 | 2015-01-07 | 湖南省电力公司调度通信局 | 电网在线风险评估方法 |
-
2013
- 2013-09-02 CN CN201310390705.4A patent/CN104008221B/zh not_active Expired - Fee Related
Cited By (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104376157A (zh) * | 2014-11-03 | 2015-02-25 | 哈尔滨汽轮机厂有限责任公司 | 大型汽轮发电机组轴系动静态特性分析计算系统 |
CN104376157B (zh) * | 2014-11-03 | 2017-09-19 | 哈尔滨汽轮机厂有限责任公司 | 大型汽轮发电机组轴系动静态特性分析计算系统 |
CN105260548A (zh) * | 2015-10-21 | 2016-01-20 | 国网湖南节能服务有限公司 | 一种基于机组实际特性的汽轮机模型建模方法 |
CN105260548B (zh) * | 2015-10-21 | 2018-08-17 | 国网湖南节能服务有限公司 | 一种基于机组实际特性的汽轮机模型建模方法 |
CN116950962A (zh) * | 2023-09-21 | 2023-10-27 | 福建福清核电有限公司 | 油动机带载运行调节的方法及其装置、电子设备 |
CN116950962B (zh) * | 2023-09-21 | 2024-01-12 | 福建福清核电有限公司 | 油动机带载运行调节的方法及其装置、电子设备 |
Also Published As
Publication number | Publication date |
---|---|
CN104008221B (zh) | 2018-07-31 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN104734147B (zh) | 一种综合能源系统概率能量流分析方法 | |
CN105956350B (zh) | 一种抽水蓄能机组过水系统的建模方法 | |
Al-Dabbas | A performance analysis of solar chimney thermal power systems | |
CN103336911B (zh) | 适用于多机电力系统调频分析的数学模型机理建模方法 | |
CN104009247A (zh) | 一种固体氧化物燃料电池电堆温度分布估计方法 | |
Borzuei et al. | On the performance enhancement of the three-blade Savonius wind turbine implementing opening valve | |
CN104008221A (zh) | 一种汽轮机模型的建立方法 | |
CN104775856A (zh) | 一种联合循环发电机组汽轮机快速启动控制方法和装置 | |
CN104732111B (zh) | 一种地源热泵的高效的逐时数值模拟方法 | |
Ji et al. | Heat transfer performance analysis of front-end capillary heat exchanger of a subway source heat pump system | |
CN105302948B (zh) | 单相流体网络建模方法 | |
CN109002741A (zh) | 一种压水堆核电机组一、二回路系统传递功率模拟方法及系统 | |
Wang et al. | One-way fluid-structure interaction simulation of an offshore wind turbine | |
CN105552934A (zh) | 时滞分布依赖的网络化电力系统负荷频率控制方法 | |
CN104504224A (zh) | 一种水轮机组暂态模型的仿真方法 | |
CN104950261A (zh) | 电池的硬件在环仿真测试方法及系统 | |
CN104809512B (zh) | 一种电力系统不同紧急控制措施性能指标的快速评估方法 | |
Islam et al. | Optimum blade design of pinwheel type horizontal axis wind turbine for low wind speed areas | |
CN101546179B (zh) | 一种火电机组过热器的非线性仿真装置 | |
CN104731761B (zh) | 天然气管网仿真方法和装置 | |
Bayareh | Numerical simulation of a solar chimney power plant in the southern region of Iran | |
Chaudhry et al. | Computational analysis to factor wind into the design of an architectural environment | |
Chen et al. | Optimization of two-stage closure law of turbine wicket gates and its application | |
Liu et al. | Numerical simulation and analysis of the vertical and double pipe soil-air heat exchanger | |
CN105093933A (zh) | 一种确定lpv变增益控制器的方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20180731 Termination date: 20190902 |
|
CF01 | Termination of patent right due to non-payment of annual fee |