CN102880746B - 汽轮机转子体积平均温度计算的惯性环节系数的确定方法 - Google Patents

汽轮机转子体积平均温度计算的惯性环节系数的确定方法 Download PDF

Info

Publication number
CN102880746B
CN102880746B CN201210323783.8A CN201210323783A CN102880746B CN 102880746 B CN102880746 B CN 102880746B CN 201210323783 A CN201210323783 A CN 201210323783A CN 102880746 B CN102880746 B CN 102880746B
Authority
CN
China
Prior art keywords
temperature
inertial element
rotor
coefficient
inertial
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
Application number
CN201210323783.8A
Other languages
English (en)
Other versions
CN102880746A (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 Power Equipment Research Institute Co Ltd
Original Assignee
Shanghai Power Equipment Research Institute 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 Power Equipment Research Institute Co Ltd filed Critical Shanghai Power Equipment Research Institute Co Ltd
Priority to CN201210323783.8A priority Critical patent/CN102880746B/zh
Publication of CN102880746A publication Critical patent/CN102880746A/zh
Application granted granted Critical
Publication of CN102880746B publication Critical patent/CN102880746B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Landscapes

  • Control Of Turbines (AREA)

Abstract

本发明提供了一种汽轮机转子体积平均温度计算的惯性环节系数的确定方法，其特征在于：该方法由以下6个步骤组成：用汽轮机内缸内壁温度作为转子外表面温度；利用有限元计算获得转子温度场分布；采用3个惯性环节；确定目标函数；确定采用的惯性环节的时间常数和系数；验证惯性环节的时间常数和系数。本发明提供的方法克服了现有技术的不足，利用离线有限元计算的体积平均温度确定惯性环节的个数、时间常数和系数，提高了转子体积平均温度和转子表面的温差及热应力计算的准确性，为汽轮机转子的热应力优化控制和汽轮机变负荷运行的安全性提供了技术保障。

Description

汽轮机转子体积平均温度计算的惯性环节系数的确定方法
技术领域
本发明涉及一种汽轮机转子体积平均温度计算的惯性环节系数的确定方法，属于汽轮机技术领域。
背景技术
汽轮机在启动、停机、负荷变动时，转子的表面的蒸汽温度变化剧烈，而转子内部的温度变化相对较慢，转子的外表面和转子体积平均温度之间的温差较大。转子的外表面和转子体积平均温度之间的温差越大，外表面的热应力就越大，转子的寿命损耗也越大。通过监控转子的外表面和转子体积平均温度之间的温差，可以监控外表面的热应力。由于转子的旋转，无法安装测温元件进行直接测量温度。需要采用计算方法获得转子温度分布。根据非稳态传热原理，转子的热传递可采用经典控制理论中的惯性环节来模拟。工程上，仅采用一个惯性环节无法完成模拟，采用3个以上并列的惯性环节才能达到要求。目前，多个惯性环节的时间常数和权重系数的确定，还没有合适的方法。
发明内容
本发明要解决的技术问题是提供一种用于汽轮机转子体积平均温度计算的多个惯性环节的时间常数和权重系数的确定方法。
为了解决上述技术问题，本发明的技术方案是提供一种汽轮机转子体积平均温度计算的惯性环节系数的确定方法，其特征在于：该方法由以下6个步骤组成：
步骤1：用汽轮机内缸内壁温度作为转子外表面温度；
根据有限元计算，在冷态启动、温态启动、热态启动、极热态启动、停机、负荷变动条件下，汽轮机内缸内壁温度与转子外表面温度的差别不超过5℃，在工程上认为，用汽轮机内缸内壁温度作为转子外表面温度参与后续计算不会造成明显的误差；
步骤2：利用有限元软件，在冷态、热态、极热态启动条件，计算转子的瞬态的温度场分布，获得转子体积平均温度随时间变化的数据；
步骤3：采用3个惯性环节；
步骤4：确定目标函数objfun
要获得满足条件的惯性环节的常数和系数，需要先确定目标函数，目标函数就是为完成目标而编写的一个要求函数值最小化的函数；通过不断调整常数和系数，搜索到目标函数的最小值，完成目标，有限元方法计算的体积平均温度和惯性环节法计算的体积平均温度的差值的平方和就是一个可用的目标函数objfun；在目标函数中，根据惯性环节的迭代计算公式
y n = y n - 1 + τ a T ( x n - y n - 1 )
编写惯性环节输入输出关系；即可根据前一次的输出值yn-1和本次的输入值xn，近似计算出本次的输出值yn；其中，τa为两次计算时间间隔，T为惯性环节的时间常数，与转子材料、体积平均温度处的轴半径及转子平均温度本身有关；
T = 3600 5.783 × R 2 A 3 T m 2 + A 2 T m + A 1
其中，R为计算体积平均温度处的轴半径，A1、A2、A3为计算惯性时间的系数，与热传导率有关；Tm是迭代计算前一次获得的体积平均温度。
步骤5：确定采用的惯性环节的时间常数和系数x
在MATLAB软件中，利用基于模式的直接搜索工具箱中的搜索函数来完成惯性环节的常数和系数的确定；在MATLAB软件的命令行输入下面的语句，完成有约束问题模式搜索；
[x fval]＝patternsearch(objfun，x0，[]，[]，[]，[]，lb，ub，options)
其中，objfun为目标函数；x0为参与模式搜索的初始行向量；lb和ub分别表示变量的下边界和上边界；options是利用psoptimset函数创建的参数结构。fval为目标函数的最终值；x为获得目标函数的最终值对应的行向量；将惯性环节的时间常数、系数的初始值依次放置到x0行向量中，搜索完成，从x行向量中取出对应的惯性环节的时间常数和系数的优化值；
步骤6：验证惯性环节的时间常数和系数x
以上步骤2～步骤5是针对冷态、热态、极热态启动条件进行的，因此需要对温态启动、停机、负荷变动条件下的情况进行相应的验证；利用步骤5获得的惯性环节的时间常数和系数，采用惯性环节法，计算温态启动、停机、负荷变动条件下的体积平均温度，并与有限元计算结果进行比较；如果温度差值较大，可调整惯性环节的个数，重新从步骤3开始；如果温度差值均满足要求，汽轮机转子体积平均温度计算的惯性环节的个数、时间常数和系数完成确定。
本发明提供的一种汽轮机转子体积平均温度计算的惯性环节系数的确定方法用汽轮机转子有限元温度场计算的体积平均温度作为标准值，来确定惯性环节计算法计算汽轮机转子体积平均温度的并联的惯性环节的个数，以及每个惯性环节的时间常数和权重系数。通过通用的MATLAB软件自带的基于模式的直接搜索工具箱中的搜索函数来完成优化搜索，不需要编制复杂的优化算法，就可通过离线有限元计算的体积平均温度来确定惯性环节的个数、时间常数和系数，且获得的优化结果可通过对其它工作条件(启、停和负荷变动)下的适用性进行验证。
本发明提供的方法克服了现有技术的不足，利用离线有限元计算的体积平均温度确定惯性环节的个数、时间常数和系数，提高了转子体积平均温度和转子表面的温差及热应力计算的准确性，为汽轮机转子的热应力优化控制和汽轮机变负荷运行的安全性提供了技术保障。
附图说明
图1为本发明提供的一种汽轮机转子体积平均温度计算的惯性环节系数的确定方法的流程图；
图2为本实施例中启停工况转子外表面与内缸内壁温度变化曲线图；
图3为本实施例中转子体积平均温度惯性环节计算结构图；
具体实施方式
为使本发明更明显易懂，兹以一优选实施例，并配合附图作详细说明如下。
图1为本发明提供的一种汽轮机转子体积平均温度计算的惯性环节系数的确定方法的流程图，所述的一种汽轮机转子体积平均温度计算的惯性环节系数的确定方法通过有限元计算，获得各种启停和负荷变动条件下汽轮机转子体积平均温度的变化曲线。把汽轮机冷态、热态、极热态启动条件下的转子外表面温度和转子体积平均温度的变化曲线的数据放到优化目标函数中。优化目标函数中根据转子外部面温度变化曲线，采用并联的多个惯性环节，计算出转子体积平均温度。进而计算出对应时刻有限元计算的转子体积平均温度与惯性环节计算的转子体积平均温度的差值的平方和。基于模式的直接搜索算法搜索出优化目标函数的最小值及其相应的惯性环节的常数和系数。通过对其它启、停和负荷变动条件的验证计算，决定是否满足要求。
以某1100MW汽轮机为例，主蒸汽压力30MPa，主蒸汽温度610℃，再热蒸汽温度610℃。汽轮机转子体积平均温度计算的惯性环节系数的确定方法由以下6个步骤组成：
步骤1：汽轮机内缸内壁温度作为转子外表面温度
采用有限元方法，分别计算冷态、温态、热态、极热态、停机过程的高压内缸内表面和高压转子外表面的温度值，绘制温度变化曲线如图2。高压内缸内表面与高压转子外表面的温差不超过5℃，在后续计算可以用高压内缸内表面温度代替高压转子外表面的温度；
步骤2：利用有限元计算获得转子温度场分布
利用有限元软件分别计算冷态、热态、极热态启动过程，汽轮机转子的体积平均温度的计算结果如表1所示，表1中的数据将用于确定惯性环节法的常数和系数。
表1冷态、热态、极热态启动过程的转子体积平均温度
步骤3：确定采用的惯性环节的个数
1个惯性环节由于只有1个待定参数，灵活性不够，多个惯性环节并联使用能更好模拟转子的传热过程，采用图2所示的惯性环节法，采用3个惯性环节。其中K1为第一惯性环节的权重系数、K2为第二惯性环节的权重系数、K3为第三惯性环节的权重系数(3个并联惯性环节的权重系数之和等于1)，E2为第二惯性环节时间常数的比例系数、E3为第三惯性环节时间常数的比例系数，PT1为惯性环节。转子外表面温度Ta和时间常数T为输入量，转子体积平均温度Tm为输出量。
步骤4：确定目标函数objfun
要获得惯性环节的常数和系数，需要先确定目标函数，目标函数就是为完成目标而编写的一个要求函数值最小化的目标函数；通过不断调整常数和系数，搜索目标函数的最小值，完成目标，有限元方法计算的体积平均温度和惯性环节法计算的体积平均温度的差值的平方和就是一个可用的目标函数objfun；
需要通过搜索确定的变量有7个，分别是第一惯性环节时间常数二次项系数A3、第一惯性环节时间常数一次项系数A2、第一惯性环节时间常数的常数项A1、第二惯性环节时间常数的比例系数E2、第三惯性环节时间常数的比例系数E3、第二惯性环节的权重系数K2、第三惯性环节的权重系数K3。以这7个变量作为目标函数的输入参数；在目标函数中，根据惯性环节的迭代计算公式
y n = y n - 1 + τ a T ( x n - y n - 1 )
编写惯性环节输入输出关系；即可根据前一次的输出值yn-1和本次的输入值xn，近似计算出本次的输出值yn；其中，τa为两次计算时间间隔，T为惯性环节的时间常数，与转子材料、体积平均温度处的轴半径及转子平均温度本身有关，
T = 3600 5.783 × R 2 A 3 T m 2 + A 2 T m + A 1
R为计算体积平均温度处的轴半径；Tm是迭代计算前一次获得的体积平均温度。
利用惯性环节顺序计算冷态、热态、极热态启动过程的体积平均温度和对应有限元计算的体积平均温度的差的平方和，以累积值作为目标函数的函数输出值。
步骤5：确定采用的惯性环节的时间常数和系数x
在MATLAB软件中，利用基于模式的直接搜索工具箱中的搜索函数来完成惯性环节的常数和系数的确定；在MATLAB软件的命令行输入下面的语句，完成有约束问题模式搜索；
 options＝psoptimset(′CompletePoll′，′on′，′Cache′，′on′，′MaxIter′，4000，′PlotFcns′，psplotbestf，′Display′，′iter′)
[x fval]＝patternsearch(objfun，[50000，-20，-0.5，50.0，50.0，0.3，0.3]，[]，[]，[]，[]，[20000.0 -50.0 -1.0 1.0 1.0 0.0 0.0]，[100000 1.0 0.0 100.0 100.0 1.0 1.0]，options)
其中，options是利用psoptimset函数创建的参数结构，objfun为目标函数；
初始目标函数值为1486920；
1000次迭代计算后目标函数值为427.6765；
2000次迭代计算后目标函数值为85.0265；
3000次迭代计算后目标函数值为45.2730；
4000次迭代计算后目标函数值为37.7240；
获得的常数和系数分别是：
[26256.65625，-11.15625，0，84.7421875，7.25，0.11005859375，0.1740234375]
步骤6：验证惯性环节的时间常数和系数x
以上步骤2～步骤5是针对冷态、热态、极热态启动条件进行的，因此需要对温态启动、停机、负荷变动条件下的适用性进行相应的验证；利用步骤5获得的惯性环节的时间常数和系数，采用惯性环节法，计算得到温态启动的汽轮机转子的体积平均温度如表2所示，并与有限元计算的转子体积平均温度数据进行比较。两者的温差不超过1.0℃，表明步骤5获得的常数和系数满足精度要求。
表2惯性环节法和有限元法计算的转子温态启动体积平均温度
  时间(s)   温度(℃)(有限元法)   温度(℃)(惯性环节法)   温差(℃)
  0   436.20   436.78   0.58
  99   432.79   432.87   0.08
  198   431.30   431.22   0.08
  301   430.42   430.44   0.02
  400   429.71   429.82   0.11
  499   429.14   429.28   0.14
  701   428.30   428.43   0.13
  899   427.77   427.86   0.09
  1101   427.47   427.52   0.05
  1501   427.40   427.37   0.03
  2000   428.01   427.94   0.07
  3002   430.93   430.83   0.10
  4000   435.41   435.27   0.14
  5002   441.32   441.12   0.20
  6000   448.37   448.12   0.15
  8000   464.93   464.63   0.30
  10000   483.80   483.53   0.27
  12000   503.61   503.46   0.15
  14000   522.35   522.40   0.05
  16000   539.34   539.90   0.66
  16922   544.56   544.56   0.00
综上所述，利用本发明提供的一种汽轮机转子体积平均温度计算的惯性环节系数的确定方法，可在电站汽轮机的控制系统中，采用惯性环节法实时计算转子的体积平均温度，进而可以提高获得转子体积平均温度和转子表面的温差和热应力的准确性，为汽轮机转子的热应力优化控制和汽轮机变负荷运行的安全性提供了技术保障。

Claims (1)

1.一种汽轮机转子体积平均温度计算的惯性环节系数的确定方法，其特征在于：该方法由以下6个步骤组成： 
步骤1：用汽轮机内缸内壁温度作为转子外表面温度； 
根据有限元计算，在冷态启动、温态启动、热态启动、极热态启动、停机、负荷变动条件下，汽轮机内缸内壁温度与转子外表面温度的差别不超过5℃，在工程上认为，用汽轮机内缸内壁温度作为转子外表面温度参与后续计算不会造成明显的误差； 
步骤2：利用有限元软件，在冷态、热态、极热态启动条件，计算转子的瞬态的温度场分布，获得转子体积平均温度随时间变化的数据； 
步骤3：采用3个惯性环节； 
步骤4：确定目标函数objfun 
要获得满足条件的惯性环节的常数和系数，需要先确定目标函数，目标函数就是为完成目标而编写的一个要求函数值最小化的函数；通过不断调整常数和系数，搜索到目标函数的最小值，完成目标，有限元方法计算的体积平均温度和惯性环节法计算的体积平均温度的差值的平方和就是一个可用的目标函数objfun；在目标函数中，根据惯性环节的迭代计算公式 
编写惯性环节输入输出关系；即可根据前一次的输出值yn-1和本次的输入值xn，近似计算出本次的输出值yn；其中，τa为两次计算时间间隔，T为惯性环节的时间常数，与转子材料、体积平均温度处的轴半径及转子平均温度本身有关； 
其中，R为计算体积平均温度处的轴半径，A1、A2、A3为计算惯性时间的系数，与热传导率有关；Tm是迭代计算前一次获得的体积平均温度； 
步骤5：确定采用的惯性环节的时间常数和系数x 
在MATLAB软件中，利用基于模式的直接搜索工具箱中的搜索函数来完成惯性环节的常数和系数的确定；在MATLAB软件的命令行输入下面的语句，完 成有约束问题模式搜索； 
[x fval]＝patternsearch(objfun，x0，[]，[]，[]，[]，lb，ub，options) 
其中，objfun为目标函数；x0为参与模式搜索的初始行向量；lb和ub分别表示变量的下边界和上边界；options是利用psoptimset函数创建的参数结构；fval为目标函数的最终值；x为获得目标函数的最终值对应的行向量；将惯性环节的时间常数、系数的初始值依次放置到x0行向量中，搜索完成，从x行向量中取出对应的惯性环节的时间常数、系数的优化值； 
步骤6：验证惯性环节的时间常数和系数x 
以上步骤2～步骤5是针对冷态、热态、极热态启动条件进行的，因此需要对温态启动、停机、负荷变动条件下的情况进行相应的验证；利用步骤5获得的惯性环节的时间常数和系数，采用惯性环节法，计算温态启动、停机、负荷变动条件下的体积平均温度，并与有限元计算结果进行比较；如果温度差值均满足要求，汽轮机转子体积平均温度计算的惯性环节的个数、时间常数和系数完成确定。 
CN201210323783.8A 2012-09-04 2012-09-04 汽轮机转子体积平均温度计算的惯性环节系数的确定方法 Active CN102880746B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201210323783.8A CN102880746B (zh) 2012-09-04 2012-09-04 汽轮机转子体积平均温度计算的惯性环节系数的确定方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201210323783.8A CN102880746B (zh) 2012-09-04 2012-09-04 汽轮机转子体积平均温度计算的惯性环节系数的确定方法

Publications (2)

Publication Number Publication Date
CN102880746A CN102880746A (zh) 2013-01-16
CN102880746B true CN102880746B (zh) 2015-01-28

Family

ID=47482069

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201210323783.8A Active CN102880746B (zh) 2012-09-04 2012-09-04 汽轮机转子体积平均温度计算的惯性环节系数的确定方法

Country Status (1)

Country Link
CN (1) CN102880746B (zh)

Families Citing this family (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110087515B (zh) * 2016-12-15 2022-06-14 美德乐控股公司 用于对婴儿食品进行调温的设备
CN111428414B (zh) * 2020-04-01 2022-05-17 浙江大学 用于汽轮机服役性能监控的温度场虚拟现实实时重建方法
CN114396317B (zh) * 2021-12-01 2022-12-16 上海发电设备成套设计研究院有限责任公司 核电汽轮机多目标多维度在线联合监控方法及系统

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
SU759734A1 (ru) * 1978-05-31 1980-08-30 Teplotekh Nii Устройство для контроля прогрева ротора турбины1
JP2005227858A (ja) * 2004-02-10 2005-08-25 Ricoh Co Ltd 駆動機構設計支援装置と駆動機構設計支援方法及び記憶媒体
CN101509656A (zh) * 2008-12-17 2009-08-19 中国电力科学研究院 一种超临界直流炉综合型协调控制方法
CN101769811A (zh) * 2009-01-05 2010-07-07 华北电力科学研究院有限责任公司 一种汽轮机主汽压力测量方法和测量装置

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
SU759734A1 (ru) * 1978-05-31 1980-08-30 Teplotekh Nii Устройство для контроля прогрева ротора турбины1
JP2005227858A (ja) * 2004-02-10 2005-08-25 Ricoh Co Ltd 駆動機構設計支援装置と駆動機構設計支援方法及び記憶媒体
CN101509656A (zh) * 2008-12-17 2009-08-19 中国电力科学研究院 一种超临界直流炉综合型协调控制方法
CN101769811A (zh) * 2009-01-05 2010-07-07 华北电力科学研究院有限责任公司 一种汽轮机主汽压力测量方法和测量装置

Non-Patent Citations (5)

* Cited by examiner, † Cited by third party
Title
史进渊,杨宇,邓志成等.《大功率电站汽轮机寿命预测与可靠性设计》.《大功率电站汽轮机寿命预测与可靠性设计》.2011,全文. *
李维特,黄保海,毕仲波.《热应力理论分析及应用》.《热应力理论分析及应用》.2004,全文. *
模式搜索法在最优化问题中的应用;吴兴远;《软件导刊》;20090831;第8卷(第8期);全文 *
祝建飞,姚峻,吴建平.1 000 MW超超临界汽轮机热应力监测及自动控制.《中国电力》.2009,第42卷(第6期),全文. *
黄仙,卞双,张保衡等.一种面向在线监控的汽轮机转子热应力新模型.《现代电力》.1997,第14卷(第3期),全文. *

Also Published As

Publication number Publication date
CN102880746A (zh) 2013-01-16

Similar Documents

Publication Publication Date Title
Tian et al. Electric vehicle air conditioning system performance prediction based on artificial neural network
Badami et al. Dynamic parsimonious model and experimental validation of a gas microturbine at part-load conditions
Serrano et al. A study on the internal convection in small turbochargers. Proposal of heat transfer convective coefficients
Duan et al. Nonlinear modeling of regenerative cycle micro gas turbine
He et al. Performance prediction of an air-cooled steam condenser using UDF method
Wang et al. Hydraulic resistance identification and optimal pressure control of district heating network
Jiang et al. Data reconciliation for steam turbine on-line performance monitoring
Marchionni et al. An appraisal of proportional integral control strategies for small scale waste heat to power conversion units based on Organic Rankine Cycles
Kim et al. Program development and simulation of dynamic operation of micro gas turbines
Mehrpanahi et al. Dynamic modeling of an industrial gas turbine in loading and unloading conditions using a gray box method
CN102880746B (zh) 汽轮机转子体积平均温度计算的惯性环节系数的确定方法
JP2014181570A (ja) 蒸気タービン制御装置
Ji et al. The optimization of the start-up scheduling for a 320 MW steam turbine
Wang et al. High-power mode control for triaxial gas turbines with variable power turbine guide vanes
Verda et al. Thermoeconomic approach for the analysis of control system of energy plants
RU2727839C2 (ru) Способ и система управления машиной
RU2621422C2 (ru) Система и способ тестирования показателя работы паровой турбины
He et al. Influence from the blade installation angle of the windward axial fans on the performance of an air-cooled power plant
Zhang et al. Gain scheduling control of waste heat energy conversion systems based on an LPV (linear parameter varying) model
Zhang et al. Modeling and output feedback control of automotive air conditioning system
CN105787195B (zh) 给水加热系统中外置式蒸汽冷却器进汽流量的计算方法
Dettori et al. Nonlinear Model Predictive Control strategy for steam turbine rotor stress
CN101769811A (zh) 一种汽轮机主汽压力测量方法和测量装置
Ramos et al. Gas engine heat pump system: Experimental facility and thermal evaluation for 5 different units
Hercus et al. Robust shape optimization of uncertain dense gas flows through a plane turbine cascade

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
C14 Grant of patent or utility model
GR01 Patent grant