专利名称：一种汽轮机高中压转子等效应力在线计算及控制方法
技术领域：
本发明涉及一种汽轮机高中压转子等效应力在线计算及控制方法，应用于汽轮机高中压转子等效应力在线控制和汽轮机寿命在线管理，属于汽轮机技术领域。
背景技术：
在汽轮机启动、停机和负荷变动的过程中，汽轮机进汽参数随时间发生变化，在汽轮机高中压转子中产生比较大的等效应力，若等效应力计算不准确，或等效应力控制不适当，过大的等效应力产生比较大的低周疲劳损伤，缩短了汽轮机高中压转子的使用寿命。现有的汽轮机高中压转子等效应力计算，有离线计算和在线计算两种方法。汽轮机高中压转子等效应力离线计算采用有限元数值计算，其优点是有限元数值计算结果准确，其缺点是计算模型复杂，等效应力计算速度慢，无法应用于汽轮机高中压转子等效应力的在线监测。汽轮机高中压转子等效应力在线计算采用等厚壁圆筒模型的解析法或差分法，其优点是计算模型简单，等效应力计算速度快，可以应用于汽轮机高中压转子等效应力在线监测，其缺点是等效应力计算结果的准确性相对比较低。
发明内容
本发明的目的是提供一种汽轮机高中压转子等效应力在线计算方法及等效应力在线控制的方法，能够在3秒钟内完成等效应力计算，计算速度快，等效应力计算相对误差小，实现汽轮机高中压转子等效应力准确性比较高的在线计算和在线控制。
为实现以上目的，本发明的技术方案是提供一种汽轮机高中压转子等效应力在线计算方法及等效应力控制方法，其方法为第一流程确定汽轮机高中压转子等效应力修正系数第一步采用简化模型计算名义等效应力σne采用现有技术的差分法离线计算汽轮机高中压转子的温度场，把有中心孔的汽轮机高中压转子简化为等厚壁圆筒模型，把无中心孔转子简化为等直径圆柱模型，采用差分法计算简化模型的温度场；采用现有技术的解析公式计算转子外表面的名义离心应力、名义蒸汽压应力和名义热应力，汽轮机高中压转子外表面名义等效应力σne的计算公式为σne=(σθt+σθp+σθth)2+(σzth)2-(σθt+σθp+σθth)(σzth)]]>式中σθt—转子外表面名义切向离心应力σθp—转子外表面名义蒸汽压应力σθth—转子外表面名义切向热应力σzth—转子外表面名义轴向热应力；第二步采用有限元模型离线计算等效应力σeqa采用现有技术，对于同型号汽轮机同一根高中压转子，按照实际结构尺寸建立两维或三维有限元模型，给定与步骤一相同的汽轮机的蒸汽温度、压力、转速和功率等的变化条件，采用数值计算方法，得出高中压转子外表面相同部位的等效应力σeqa；其特征在于，第三步确定瞬态等效应力修正系数ki对于汽轮机某一种启动过程、停机过程或负荷变动过程的某一瞬态时刻ti，计算等效应力修正系数kiki＝σeqa/σne第四步确定等效应力修正系数k的计算公式汽轮机高中压转子等效应力的修正系数k与汽轮机主蒸汽压力P0、汽轮机主蒸汽温度t0、汽轮机工作转速n、汽轮机高中压转子体积平均温度θm和汽轮机非稳态过程(启动过程、停机过程或负荷变动过程)起始时刻高中压内缸金属监测温度θ0等5个物理量有关。按照θ0≤100℃和θ0＞100℃，采用数理统计和多元回归技术，得出以下高中压转子等效应力的修正系数k的计算公式当θ0≤100℃时k＝A1+B1PD+C1t0+D1n+E1θm+F1θ0
θ0＞100℃时k＝A2+B2PD+C2t0+D2n+E2θm+F2θ0以上两个公式中A1、B1、C1、D1、E1、F1和A2、B2、C2、D2、E2、F2的取值示于表1[表1]
第二流程汽轮机高中压转子等效应力在线计算及控制第五步采用简化模型在线实时计算名义等效应力采用现有技术的等厚壁圆筒或等直径圆柱的简化模型，在线实时计算汽轮机高中压转子外表面名义等效应力σne；第六步在线实时计算等效应力修正系数k采用第四步得出的等效应力修正系数K的计算公式，采用C语言编写计算软件，运行在汽轮机的参数监视与控制的计算机上，输入在线监测的汽轮机主蒸汽温度t0、汽轮机工作转速n、汽轮机高中压转子体积平均温度θm和汽轮机非稳态过程即启动过程、停机过程或负荷变动过程，起始时刻高中压内缸金属监测温度θ0，实时计算得出汽轮机高中压转子等效应力修正系数k；
第七步在线实时计算等效应力修正值σeq使用计算机软件在线计算汽轮机高中压转子外表面等效应力的在线修正值σeqσeq＝k×σne第八步确定转子材料控制极限应力σu汽轮机高中压转子材料控制极限应力σu取工作温度下材料的屈服强度σ0.2和转子工作温度下转子材料工作1000小时的断裂强度σfr的最小值σu＝min{σ0.2；σfr}第九步确定相对应力RS汽轮机高中压转子相对应力RS为转子外表面等效应力在线修正值σeq与转子材料控制极限应力σu之比RS＝σeq/σu第十步等效应力的在线控制在线计算相对应力RS并与设定值比较，汽轮机高中压转子等效应力的控制分为以下三种情况若RS≤80％，汽轮机的主蒸汽温度变化率和负荷变化率按《汽轮机运行规程》的规定数值操作；若80％＜RS＜100％，减小汽轮机的主蒸汽温度变化率和负荷变化率，以减小汽轮机高中压转子的等效应力；若100％≤RS＜105％，控制汽轮机的主蒸汽温度变化率和负荷变动率为0，以进一步降低汽轮机高中压转子的等效应力。
第十一步等效应力的在线报警和打闸停机在线计算汽轮机高中压转子的相对应力Rs，若105％≤Rs＜125％，发出报警，30分钟后跳闸停机；若Rs≥125％，发出警报，1分钟后跳闸停机。
本发明给出一种汽轮机高中压转子等效应力在线计算方法及等效应力控制方法，可实现汽轮机高中压转子等效应力的在线计算及控制。如果汽轮机高中压转子的相对应力超过设定值，通过在线控制汽轮机主蒸汽的温度变化率、汽轮机的负荷变化率来减小汽轮机高中压转子的等效应力，以减小汽轮机高中压转子的低周疲劳寿命损伤。如果汽轮机高中压转子的相对应力达到报警值时，发出报警，并在设定的时间内停机，以减小汽轮机高中压转子的低周疲劳寿命损伤，达到控制汽轮机高中压转子等效应力并延长高中压转子使用寿命的技术效果。
本发明的优点是能够在3秒钟内完成等效应力计算，计算速度快，等效应力计算相对误差小，实现汽轮机高中压转子等效应力准确性比较高的在线计算和在线控制。
图1为本发明所采用方法的流程图；图2为本发明所采用的计算机软件框图；图3冷态启动过程的高压转子调节级叶轮后圆角部位等效应力的变化规律；图4温态启动过程的高压转子调节级叶轮后圆角部位等效应力的变化规律；图5热态启动过程的高压转子调节级叶轮后圆角部位等效应力的变化规律；图6极热态启动过程的高压转子调节级叶轮后圆角部位等效应力变化规律；图7冷态启动过程的高压转子调节级叶轮后圆角部位相对应力的变化规律。
具体实施例方式
如图1所示，本发明提供的方法由两个流程组成，本发明技术方案的第一步至第四步构成第一个流程，确定汽轮机高中压转子等效应力的修正系数k；本发明技术方案的第五步至第十一步构成第二流程，在线计算和控制汽轮机高中压转子的等效应力。以下结合附图和实施例对本发明作进一步说明。
实施例对于某型号进汽温度为566℃的超临界600MW汽轮机高压转子采用图1所示的本发明提供方法的流程图，采用C语言编写计算软件，运行在汽轮机的参数监视与控制的计算机上，如图2所示本发明方法的计算软件框图，计算得出的该汽轮机高压转子等效应力的修正值如图3至图6所示。
第一步和第五步采用现有技术的简化模型在线计算得出该汽轮机高压转子调节级叶轮外表面圆角部位在冷态启动过程、温态启动过程、热态启动过程、极热态启动过程的名义等效应力σne在图3至图6中用曲线1表示。
第二步采用现有技术的有限元模型，离线计算得出该汽轮机高压转子调节级叶轮外表面圆角部位在冷态启动过程、温态启动过程、热态启动过程、极热态启动过程的等效应力σeqa在图3至图6中用曲线2表示。
第三步在图3至图6中，同一横坐标对应曲线2的等效应力值σeqa和曲线1的名义等效应力值σne相比，得出一系列等效应力修正系数ki。
第四步按照非稳态过程(启动过程、停机过程或负荷变动过程)起始时刻高中压内缸金属监测温度θ0≤100℃和θ0＞100℃，采用数理统计和多元回归技术，得出由汽轮机在线监测的汽轮机主蒸汽压力P0、汽轮机主蒸汽温度t0、汽轮机工作转速n、汽轮机高中压转子体积平均温度θm和汽轮机非稳态过程(启动过程、停机过程或负荷变动过程)起始时刻高中压内缸金属监测温度θ0等5个物理量计算汽轮机高压转子调节级叶轮后圆角部位等效应力修正系数k的计算公式为当汽轮机冲转时刻高压内缸金属温度当θ0≤100℃时k＝1.204062-8.190503×10-3P0+4.834444×10-4t0-1.204386×10-4n+3.920135×10-5θm+3.329205×10-4θ0当汽轮机冲转时刻高压内缸金属温度当θ0＞100℃时k＝1.542736-1.049599×10-2P0-1.156659×10-3t0-3.326461×10-5n-5.762553×10-4θm+1.309901×10-3θ0第六步和第七步使用以上两个计算公式输入在线监测的汽轮机主蒸汽温度t0、汽轮机工作转速n、汽轮机高中压转子体积平均温度θm和汽轮机非稳态过程(启动过程、停机过程或负荷变动过程)起始时刻高中压内缸金属监测温度θ0，实时计算得出汽轮机高中压转子等效应力修正系数k。按照σeq＝k×σne，在线实时计算得出的该汽轮机高压转子调节级叶轮后外表面圆角部位的等效应力的修正值σeq在图3至图6中用曲线3表示。
第八步该高中压转子采用CrMoV材料，其屈服极限σ0.2与工作温度t的关系为σ0.2＝544.963-0.295t(MPa)，其断裂强度σfr＝10c，c＝[23×(273+t)-3.4498×104]/(-6.5356×103)。
第九步在线计算的该汽轮机高压转子调节级叶轮后外表面圆角部位在冷态启动过程中相对应力RS的变化规律如图7所示。
第十步和第十一步在该汽轮机高压转子相对应力RS变化规律的图7中，A区RS≤80％，按现有运行规程规定操作；B区80％＜RS＜100％，汽轮机运行人员减小汽轮机主蒸汽温度变化率和负荷变化率；C区100％≤RS＜105％，汽轮机运行人员把主蒸汽温度变化率和负荷变动率取为0；D区105％≤RS＜125％，发出报警，30分钟后跳闸停机；E区RS＞125％，发出报警，1分钟后跳闸停机。
从图3知，同有限元离线计算结果相比，采用简化模型在线计算，冷态启动过程等效应力相对误差的范围为-25.14％至10.85％；采用本发明提供方法在线计算，冷态启动过程等效应力的相对误差的范围为-9.84％至10.63％。从图4知，同有限元离线计算结果相比，采用简化模型在线计算，温态启动过程等效应力相对误差的范围为-33.38％至7.77％；采用本发明提供方法在线计算，温态启动过程等效应力的相对误差的范围为-3.25％至6.27％。从图5知，同有限元离线计算结果相比，采用简化模型在线计算，热态启动过程等效应力相对误差的范围为-34.42％至-17.48％；采用本发明提供方法在线计算，热态启动过程等效应力的相对误差的范围为-6.21％至5.04％。从图6知，同有限元离线计算结果相比，采用简化模型在线计算，极热态启动过程等效应力相对误差的范围为-33.86％至-26.20％；采用本发明提供方法在线计算，极热态启动过程最大等效应力的相对误差的范围为-5.23％至4.62％。
使用现有技术的简化模型在线计算汽轮机高压转子的等效应力，同有限元离线计算结果相比，汽轮机启动过程中高压转子外表面等效应力的相对误差的范围为-35％～11％。使用本发明提供的汽轮机高压转子等效应力的在线计算方法，同有限元离线计算结果相比，汽轮机启动过程中高压转子外表面等效应力计算相对误差的范围为-10％至11％，本发明提供的汽轮机高压转子等效应力在线计算的相对误差的绝对值约为简化模型的1/3，达到了提高汽轮机高中压转子等效应力在线计算准确性的技术效果。使用本发明提供的在线计算准确性比较高的等效应力控制汽轮机的启动过程、停机过程和负荷变动过程的汽轮机的主蒸汽温度变化率和汽轮机的负荷变化率，可以比较准确地监视和控制运行中汽轮机高中压转子的等效应力；可以实现减小汽轮机高中压转子的低周疲劳寿命损伤，达到延长汽轮机高中压转子使用寿命的技术效果，为汽轮机高中压转子寿命在线管理提供了技术依据。
权利要求
1.一种汽轮机高中压转子等效应力在线计算方法及等效应力控制方法，其方法为第一流程确定汽轮机高中压转子等效应力修正系数第一步采用简化模型计算名义等效应力σne采用现有技术的差分法离线计算汽轮机高中压转子的温度场，把有中心孔的汽轮机高中压转子简化为等厚壁圆筒模型，把无中心孔转子简化为等直径圆柱模型，采用差分法计算简化模型的温度场；采用现有技术的解析公式计算转子外表面的名义离心应力、名义蒸汽压应力和名义热应力，汽轮机高中压转子外表面名义等效应力σne的计算公式为σne=(σθt+σθp+σθth)2+(σzth)2-(σθt+σθp+σθth)(σzth)]]>式中σθt-转子外表面名义切向离心应力σθp-转子外表面名义蒸汽压应力σθth-转子外表面名义切向热应力σzth-转子外表面名义轴向热应力；第二步采用有限元模型离线计算等效应力σeqa采用现有技术，对于同型号汽轮机同一根高中压转子，按照实际结构尺寸建立两维或三维有限元模型，给定与步骤一相同的汽轮机的蒸汽温度、压力、转速和功率等的变化条件，采用数值计算方法，得出高中压转子外表面相同部位的等效应力σeqa；其特征在于，第三步确定瞬态等效应力修正系数ki对于汽轮机某一种启动过程、停机过程或负荷变动过程的某一瞬态时刻ti，计算等效应力修正系数kiki＝σeqa/σne第四步确定等效应力修正系数k的计算公式汽轮机高中压转子等效应力的修正系数k与汽轮机主蒸汽压力P0、汽轮机主蒸汽温度t0、汽轮机工作转速n、汽轮机高中压转子体积平均温度θm和汽轮机非稳态过程(启动过程、停机过程或负荷变动过程)起始时刻高中压内缸金属监测温度θ0等5个物理量有关。按照θ0≤100℃和θ0＞100℃，采用数理统计和多元回归技术，得出以下高中压转子等效应力的修正系数k的计算公式当θ0≤100℃时k＝A1+B1PD+C1t0+D1n+E1θm+F1θ0θ0＞100℃时k＝A2+B2PD+C2t0+D2n+E2θm+F2θ0以上两个公式中A1、B1、C1、D1、E1、F1和A2、B2、C2、D2、E2、F2的取值示于表1[表1]
第二流程汽轮机高中压转子等效应力在线计算及控制第五步采用简化模型在线实时计算名义等效应力采用现有技术的等厚壁圆筒或等直径圆柱的简化模型，在线实时计算汽轮机高中压转子外表面名义等效应力σne；第六步在线实时计算等效应力修正系数k采用第四步得出的等效应力修正系数K的计算公式，采用C语言编写计算软件，运行在汽轮机的参数监视与控制的计算机上，输入在线监测的汽轮机主蒸汽温度t0、汽轮机工作转速n、汽轮机高中压转子体积平均温度θm和汽轮机非稳态过程即启动过程、停机过程或负荷变动过程，起始时刻高中压内缸金属监测温度θ0，实时计算得出汽轮机高中压转子等效应力修正系数k；第七步在线实时计算等效应力修正值σeq使用计算机软件在线计算汽轮机高中压转子外表面等效应力的在线修正值σeqσeq＝k×σne第八步确定转子材料控制极限应力σu汽轮机高中压转子材料控制极限应力σu取工作温度下材料的屈服强度σ0.2和转子工作温度下转子材料工作1000小时的断裂强度σfr的最小值σu＝min{σ0.2；σfr}第九步确定相对应力RS汽轮机高中压转子相对应力RS为转子外表面等效应力在线修正值σeq与转子材料控制极限应力σu之比RS＝σeq/σu第十步等效应力的在线控制在线计算相对应力RS并与设定值比较，汽轮机高中压转子等效应力的控制分为以下三种情况若RS≤80％，汽轮机的主蒸汽温度变化率和负荷变化率按《汽轮机运行规程》的规定数值操作；若80％＜RS＜100％，减小汽轮机的主蒸汽温度变化率和负荷变化率，以减小汽轮机高中压转子的等效应力；若100％≤RS＜105％，控制汽轮机的主蒸汽温度变化率和负荷变动率为0，以进一步降低汽轮机高中压转子的等效应力。第十一步等效应力的在线报警和打闸停机在线计算汽轮机高中压转子的相对应力RS，若105％≤RS＜125％，发出报警，30分钟后跳闸停机；若RS≥125％，发出警报，1分钟后跳闸停机。
全文摘要
本发明给出一种汽轮机高中压转子等效应力在线计算方法及等效应力控制方法，可实现汽轮机高中压转子等效应力的在线计算及控制。如果汽轮机高中压转子的相对应力超过设定值，通过在线控制汽轮机主蒸汽的温度变化率、汽轮机的负荷变化率来减小汽轮机高中压转子的等效应力，以减小汽轮机高中压转子的低周疲劳寿命损伤。如果汽轮机高中压转子的相对应力达到报警值时，发出报警，并在设定的时间内停机，以减小汽轮机高中压转子的低周疲劳寿命损伤，达到控制汽轮机高中压转子等效应力并延长高中压转子使用寿命的技术效果。
文档编号F01D21/00GK1908381SQ20061003024
公开日2007年2月7日 申请日期2006年8月21日 优先权日2006年8月21日
发明者史进渊, 杨宇, 邓志成, 孙庆 申请人:上海发电设备成套设计研究院