（热能工程专业论文）600mw汽轮机启动优化的研究.pdf

华北电力大学硕士学位论文 摘要 本文以哈汽6 0 0 m w 汽轮机的启动过程为研究对象，分析与探讨了启动过程中 蒸汽温升率的计算方法，并在此基础上研究了蒸汽初温与转子金属温度的匹配问 题，使得汽轮机启动过程优化。同时对启动过程中的换热系数进行了计算与比较。 在启动过程中须要对应力进行监测，因此本文对汽轮机转子应力的在线计算模型进 行了研究，针对目前汽轮机转子温度场的在线求解问题，给出了克兰克尼科尔森差 分计算模型，并对6 0 0 m w 汽轮机冷态启动过程进行了仿真计算，并验证该计算模 型的合理性。该方法具有数据处理实时、稳定性好等特点。 关键词：启动，寿命分配，换热系数，克兰科一尼科尔森差分 a b s t r a c t s t a r tp r o c e s so f6 0 0 m wt u r b i n em a d eb yh a r b i nt u r b i n ec o m p a n yi so p t i m i z e di n t h i sp a p e r t h es p e e do fs t e a mt e m p e r a t u r er i s i n gi ns t a r tp r o c e s si sa n a l y z e da n d d i s c u s s e d a n da c c o r d i n gt h i s ，t h em a t c h i n go fs t e a mt e m p e r a t u r ea n dm e t a lt e m p e r a t u r e i sc a l c u l a t e di no r d e rt oo p t i m i z et h es t a r tp r o c e s s c o n v e c t i o nh o tc o e f f i c i e n to ft h es t a r t p r o c e s si sa l s oc a l c u l a t e dw i t hf o r m u l a so fd i f f e r e n tc o u n t r i e sa n dc o m p a n i e sa n dt h e y a r ec o m p a r e d t h es t r e s so ft h er o t o rs h o u l db em o n i t o r e da n dc r a n k n i c o l s o n f i n i t e d i f f e r e n c em e t h o di su s e dt oc a l c u l a t et h et e m p e r a t u r ef i e l di nt h er o t o ra r t i f i c i a l l y d u r i n gt h es u p e r c r i t i c a l6 0 0 m ws t e a mt u r b i n ec o l ds t a r t u p t h r o u g ht h ee m u l a t i o na n d t h ec o m p a r i s o nw i t ht h er e s u l to ff i n i t ee l e m e n tm e t h o d ，t h ev a l i d i t yo ft h em o d e li s v e r i f i e d t h i sm e t h o dh a st h ea d v a n t a g eo fd e a l i n gw i t hd a t af a s ta n da c c u r a t e l y , a n dt h e s t a b i l i t yo ft h ec a l c u l a t i o ni sp e r f o r m e dw e l l c a iy u ( t h e r m a lp o w e re n g i n e e r i n g ) d i r e c t e db ya s s o c i a t ep r o f z h up i n g k e yw o r d ：s t a r t u p ，l i f ed i s t r i b u t i o n ，h e a tt r a n s f e rc o e f f i c i e n t ，c r a n k n i c o l s o nf i n i t e - d i f f e r e n c e 华北电力大学硕士学位论文 摘要 本文以哈汽6 0 0 m w 汽轮机的启动过程为研究对象，分析与探讨了启动过程中 蒸汽温升率的计算方法，并在此基础上研究了蒸汽初温与转子金属温度的匹配问 题，使得汽轮机启动过程优化。同时对启动过程中的换热系数进行了计算与比较。 在启动过程中须要对应力进行监测，因此本文对汽轮机转子应力的在线计算模型进 行了研究，针对目前汽轮机转子温度场的在线求解问题，给出了克兰克尼科尔森差 分计算模型，并对6 0 0 m w 汽轮机冷态启动过程进行了仿真计算，并验证该计算模 型的合理性。该方法具有数据处理实时、稳定性好等特点。 关键词：启动，寿命分配，换热系数，克兰科一尼科尔森差分 a b s t r a c t s t a r tp r o c e s so f6 0 0 m wt u r b i n em a d eb yh a r b i nt u r b i n ec o m p a n yi so p t i m i z e di n t h i sp a p e r t h es p e e do fs t e a mt e m p e r a t u r er i s i n gi ns t a r tp r o c e s si sa n a l y z e da n d d i s c u s s e d a n da c c o r d i n gt h i s ，t h em a t c h i n go fs t e a mt e m p e r a t u r ea n dm e t a lt e m p e r a t u r e i sc a l c u l a t e di no r d e rt oo p t i m i z et h es t a r tp r o c e s s c o n v e c t i o nh o tc o e f f i c i e n to ft h es t a r t p r o c e s si sa l s oc a l c u l a t e dw i t hf o r m u l a so fd i f f e r e n tc o u n t r i e sa n dc o m p a n i e sa n dt h e y a r ec o m p a r e d t h es t r e s so ft h er o t o rs h o u l db em o n i t o r e da n dc r a n k n i c o l s o n f i n i t e d i f f e r e n c em e t h o di su s e dt oc a l c u l a t et h et e m p e r a t u r ef i e l di nt h er o t o ra r t i f i c i a l l y d u r i n gt h es u p e r c r i t i c a l6 0 0 m ws t e a mt u r b i n ec o l ds t a r t u p t h r o u g ht h ee m u l a t i o na n d t h ec o m p a r i s o nw i t ht h er e s u l to ff i n i t ee l e m e n tm e t h o d ，t h ev a l i d i t yo ft h em o d e li s v e r i f i e d t h i sm e t h o dh a st h ea d v a n t a g eo fd e a l i n gw i t hd a t af a s ta n da c c u r a t e l y , a n dt h e s t a b i l i t yo ft h ec a l c u l a t i o ni sp e r f o r m e dw e l l c a iy u ( t h e r m a lp o w e re n g i n e e r i n g ) d i r e c t e db ya s s o c i a t ep r o f z h up i n g k e yw o r d ：s t a r t u p ，l i f ed i s t r i b u t i o n ，h e a tt r a n s f e rc o e f f i c i e n t ，c r a n k n i c o l s o nf i n i t e - d i f f e r e n c e 声明尸明 本人郑重声明：此处所提交的硕士学位论文铝口咖n 簏彩锄嬲华北电 力大学攻读硕士学位期间，在导师指导下进行的研究工作和取得的研究成果。据本人所 知，除了文中特别加以标注和致谢之处外，论文中不包含其他人已经发表或撰写过的研 究成果，也不包含为获得华北电力大学或其他教育机构的学位或证书而使用过的材料。 与我一同工作的同志对本研究所做的任何贡献均已在论文中作了明确的说明并表示了 谢意。 - t ：卜一 学位论文作者签名：亟兰日期：蜊矿 关于学位论文使用授权的说明 本人完全了解华北电力大学有关保留、使用学位论文的规定，即：学校有权保管、 并向有关部门送交学位论文的原件与复印件；学校可以采用影印、缩印或其它复制手 段复制并保存学位论文；学校可允许学位论文被查阅或借阅；学校可以学术交流为 目的，复制赠送和交换学位论文；同意学校可以用不同方式在不同媒体上发表、传播学 位论文的全部或部分内容。 ( 涉密的学位论文在解密后遵守此规定) 作者签名：惠l e l 期：q2 三。堡 导师签名： 竺军 日+ 期：咯：拿：2 9 华北电力大学硕士学位论文 1 1 课题背景及意义 第一章引言帚一早亏l 函 近年来，我国的电力工业发展十分迅速，供电能力大幅度提高，电网容量不断 增大，用电结构也相应变化，电力供求之间矛盾也日益突出，电网峰谷差也日益加 剧，迫使大型火电机组频繁的参与调峰运行。而调峰过程中，机组的工作条件恶化， 机组的寿命损耗和安全性成为影响调峰运行能力的重要因素。近年来针对大中型火 电机组参与调峰运行的可行性，各种不同调峰运行方式的经济性和安全性，都进行 广泛的实验和研究，取得了一定的成果，但由于实际机组运行工况的复杂性，目前 仍有许多问题还需进一步深入研究。对机组过渡工况下的状态进行研究，提高机组 调峰运行的经济性和适应性，是当前需要解决的主要问题n 】。 现在国产大型机组，多数是以带基本负荷设计的，主辅机均难以适应大幅度调 峰运行的要求，限制了调峰运行中负荷变化的幅度和速率。机组在调峰运行的启动、 停机和变负荷过程中，各处蒸汽参数不断变化，其转子和汽缸的金属温度和应力随 之变化。对于汽缸这个厚壁部件，由于机组高压缸的设计普遍采用了双层结构，而 且汽缸壁的金属厚度较转子薄，蒸汽对汽缸内壁的换热系数也远比转子小，因而启 动时的径向温差及热应力都远比转子小，且转子长期在高温区工作，受力情况很复 杂，除热应力外，还承受着各种机械应力，因而监视转子应力情况更具必要性。参 与调峰运行的机组，在工况变化的过程中，其工作状态不断发生变化，使蒸汽与金 属之间产生剧烈变化的换热，造成部件受热不均匀，形成不均匀温度场，使汽轮机 的汽缸和转子内产生很大的热应力。这种频繁启停或大幅度负荷变动的非稳定工 况，将导致金属材料的低周疲劳损伤，缩短机组的使用寿命。汽轮机转子是工作条 件最艰苦、受力情况最复杂的汽轮机部件，其寿命基本代表了整台汽轮机组的寿命。 为了提高汽轮机组的安全性和经济性，大容量机组的转子热应力和寿命损耗的研究 已成为人们关注的焦点。只有准确了解机组在不同运行工况下的寿命，制定合理的 运行模式，才能确保火电机组的安全经济运行。 国外电网调峰机组的容量较大，到了7 0 年代后期，一些发达国家已经投运了大 批承担中间负荷的从2 0 0m w 到6 0 0 m w 的单元机组。随着世界各国电网容量的不 断增大和负荷的峰谷差值增加，调峰机组的容量将随之变大，使目前承担基本负荷 的大型高效机组也会参与调峰运行。因此越来越多的国家开始重视调峰汽轮机的设 计和调峰方式的研究，以平衡电网中日益增长的负荷峰谷差乜，。 我国5 0 、6 0 年代投运的机组，运行已超过1 0 万小时，属于老化设备，如将其 华北电力大学硕士学位论文 全部淘汰，这在我国电量急剧增加的今天几乎不可能。即使是在发电量已相当大的 发达国家，也都在努力探讨如何延长发电厂的寿命和汽轮机零部件的寿命。因此对 于这些老化设备进行剩余寿命的预测，制定更新老化部件的计划，实施适当的维修 管理，对改善这些火电机组的老化状况，延长运行寿命也是非常必要的乜1 。 1 2 国内外研究发展状况 1 2 1 国外研究现状 由于目前转子的温度和应力尚不能直接进行测量，只能通过间接方法，建立相 应的数学模型，测量相关参数，求出转子金属温度和应力的变化及寿命损耗。现在 转子应力的数学模型大多数是采用一维温度场理论解的简化式，其计算精度较低， 只能反映应力的变化趋势，而不能得到应力的精确值。若在此基础上计算转子在启 停和变负荷过程中的寿命损耗，将会产生较大的误差口1 。 国外机组寿命管理的应用在日本、美国和欧洲较为普遍。美国自6 0 年代g o l l i n 电站汽轮机失事以后，一些大的公司和研究机构如g e 、w e s t i n g h o u s e 、e p r i 等对转子的安全性更为重视；进行了深入的研究。他们将有限元等先进数值方法由 于汽轮机转子的分析计算，对转子材料的低周疲劳、高温蠕变、低温脆性和裂纹扩 展规律等诸多方面的问题进行了大量的研究，并在汽轮发电机组上安装了应力及寿 命损耗指示器以指导机组运行。 日本在汽轮机寿命管理方面也做了很多工作，除了预测可能出现裂纹的寿命外， 还对转子剩余寿命做出计算。日本的k a g a w au n i v e r s i t y 的e b a r a 等对汽轮机动叶片 采用的1 2 c r 钢和t i 6 a i 4 v 合金的疲劳特性进行研究，f u j i y a m a ，k a z u n a r i ；t a k a k i ， k e i s u k e ；n a k a t a n i ，y u j i r o 等根据统计损伤和随机损伤仿真研究，对汽轮机设备进 行寿命评估，采用先进技术设计汽轮机流通部分，以提高机组的性能和设备的可靠 性。另外，日本在无损探伤的研究方面处于世界先进水平，日立、三菱重工、东芝、 富士机电等著名大公司相继提出脆化一腐蚀法、硬度法、金属组织法、电极化法等 无损探伤方法作为改进转子寿命评估的手段豫3 1 4 1 。 德国的w i c h t m a n n ，a n d r e a s 研究了高温对汽轮机部件的蠕变损伤；z a v i s k a ， r e i c h e l 研究了汽轮机冷态启动过程中的转子温度变化，在此基础上建立了冷态启动 仿真模型；s c h e e f e r ，m ；k n o d l e r ；s c a r l i n ，b 等对电厂抗高温、高压材料进行了探 讨，一方面是发展新的材料，一方面是在已有的材料表面喷涂抗氧化性能强的图层； 以及关于机组安全经济运行方面也进行了大量的研究乜儿3 儿“1 。 1 2 2 国内研究现状 目前有关机组调峰运行过程中的热应力变化和寿命管理方面还有若干问题没 华北电力大学硕士学位论文 有彻底解决。例如在进行机组非稳态温度状态和热应力计算中蒸汽参数和换热系数 的确定，寿命预测中我国转子用钢高温疲劳曲线的确定，都有待进一步的研究和完 善。汽轮机在高温、高压和高转速的条件下工作，实际运行中参数的变动、负荷的 波动与设计工况差别很大，若用理论值和设计值来分析汽轮机的热应力和寿命损 耗，很难真实的反映机组的实际状况。只有用实测参数来进行分析计算，才能保证 其结果的真实可信。但计算中许多所需要的参数，实地的测取有一定的困难，必须 根据运行的实际情况来进行合理的处理晡q 1 。 我国从8 0 年代初开始进行转子寿命损耗预测和寿命分配研究。多年来，我国 有关研究机构、高等院校以及制造部门、电厂针对机组调峰的需要，以国产机组为 对象，研究了汽轮机主要零部件在非稳态下的温度及热应力分布、变化规律、金属 材料的疲劳特性以及部件的寿命损耗。对国产大容量机组参与调峰运行的可行性、 调峰运行的安全性和经济性、调峰幅度进行了深入的探讨，对低负荷、少汽无功和 两班制等不同的调峰方式在经济性和安全性方面进行了理论分析和实验研究，很多 单位都相继开展了汽轮机转子应力监测和寿命损耗计算的研究工作。 清华大学的黄仙、倪维斗等提出汽轮机转子热应力的复频法建模，根据转子热 应力的形成机理，通过建立“蒸汽温度变化率一转子体全温差 与“转子体全温差 一转子体内外表面热应力 的复频域数学模型，然后由控制理论方便的得到“蒸汽 温度变化率一转子体内外表面热应力的数学模型强1 。倪维斗、黄以腾等在原有汽 轮机转子热应力解析计算式的基础上，给出了递推求解热应力值的算法。利用递推 式，不仅加快了计算速度，还使得在监控过程中用来保存历史信息所需的存储空间 维持不变阳1 。支小牧、寇可新建立了有中心孔转子热应力解析递推计算法。该方法 计入了换热系数的变化，对汽温非线性变化采用分段线性叠加的分析方法，这个对 于边界介质温度变化和换热情况复杂的有中心孔的受热部件的应力计算和研究具 有普遍意义，并应用雨流法实时监测和管理转子的寿命损耗h 1 。 上海交大的周宇阳、陈汉平等利用热应力计算过程中边界条件、物性参数和初 殆条件等模型输入参数的灰色特性，提出以区间代数为工具的灰色数学与有限单元 模型相结合所形成的灰色有限元数学模型，可以充分表达模型输入参数的灰色特性 汹3 。韩刚、陈汉平针对目前国内外汽轮机组在线寿命管理系统的监测参数的选取进 行了详细的分析研究，提出了多种温度监测来提高机组寿命管理可靠性的监测模式 【8 】 东南大学的江宁、曹祖庆等提出采用仿真模型中的权函数来计算转子的热应 力，即可保证计算精度接近二维算法，又可提高计算速度。但该方法不适用于变物 性系统，而且只能根据影响最大的主要参数进行计算，不能考虑多参数的影响，限 制了该方法的应用范围。 华中科技大学的王坤、黄树红等对1 2 5 m w 汽轮机转子的启停调峰实验进行了 3 华北电力大学硕士学位论文 研究，计算了启停调峰过程中的应力谱和低周疲劳寿命损耗，探讨了如何对转子寿 命损耗进行优化控制的问题砼1 。刘华堂、李树人等对国产2 0 0 m w 汽轮机参加调峰 运行的寿命进行了分析计算。 浙江大学的陈坚红、阮伟等采用大型结构有限元分析软件，对冷态启动温升率 对转子寿命的影响进行了研究。郭乙木、鲁祖统在工业汽轮机转子材料低周疲劳实 验及对实际工况下转子瞬态温度场、应力场分析的基础上，对转子的有效寿命进行 了预测旧。1 。 华北电力大学的张光、张宝衡等提出一种计算大型汽轮机转子热应力的二维离 散模型，这种二维离散模型也可以进一步被用来计算转子轴向热膨胀和监测机组胀 差。黄仙、杨昆等对汽轮机热应力自适应模型进行了研究。崔健、安利强等对转子 热应力计算中的有限元方法的精度问题进行了探讨。苑素玲、葛永庆等应用 m a t l a b 软件中s i m u l i n k 软件包对汽轮机转子热应力进行了实时仿真，得到实 时监控下汽轮机转子热应力对温升率的仿真曲线。安利强、王璋奇提出了有限元热 应力分析的前后处理可视化的实现方法。刘静静、杨昆通过直接测出的高中压缸内 缸壁的温度来模拟高中压转子表面温度，作为第一类边界条件，导出了汽轮机转子 温度和热应力的计算公式。安江英、杨昆、张保衡对汽轮机转子长期在高温高压状 态下材质软化现象进行了研究，并建立了硬度与蠕变和应力松弛性能之间的解析关 系式n 2 1 引。 西安交大对汽轮机转子抗疲劳的模糊可靠性理论进行了分析，主要对汽轮机叶 片的疲劳寿命预测方法进行了研究。国家电力公司热工研究院也对汽轮机在变负荷 运行的动叶片和疲劳寿命特性进行了相关的研究工作。2 0 0 2 年云南电力实验研究院 与上海成套设备研究所合作研究了“汽轮机高压转子在蠕变和疲劳作用下寿命的设 计和评定体系一，为汽轮机高压转子的寿命设计、寿命评定和寿命延长奠定了坚实 的基础3 肼1 。 1 3 本文工作简述 l 以哈汽6 0 0 m w 汽轮机的启动过程为研究对象，分析与探讨了启动过程中蒸汽温 升率的计算方法，并在此基础上研究了蒸汽初温与转子金属温度的匹配问题，使得 汽轮机启动过程优化。 2 对哈汽6 0 0 m w 汽轮机转子冷态启动过程中的换热系数进行了研究。由于不同国 家和公司的换热系数计算公式不同，本文对常用的美国西屋、前苏联和阿尔斯通公 司的换热系数计算公式进行了计算和比较，并综合不同的计算结果，采用最小二乘 法对数据进行处理，得到不同部位换热系数的计算公式。 3 建立了汽轮机转子温度场在线计算模型。本文针对目前汽轮机转子温度场的在线 华北电力大学硕十学位论文 求解问题，给出了克兰克尼科尔森差分计算模型，并对6 0 0 m w 汽轮机冷态启动过 程进行了仿真计算，并验证该模型计算准确度的可靠性。 华北电力大学硕士学位论文 第二章转子寿命损耗的研究 汽轮机运行过程中，转子承受交变应力：启动过程加热转子表面承受压应力， 停机过程为拉应力。经过一定周次的交变应力循环，金属表面将出现疲劳裂纹并逐 渐扩展以致断裂。其特点是交变周期长，频率低，疲劳裂纹萌生的循环周次少，称 为低周疲劳，不仅发生在机组的启动和停机过程，在机组大负荷变化时也会发生。 另外，由于转子长期工作在高温环境下，转子也会产生高温蠕变。因此，转子通常 处在疲劳和蠕变交互耦合作用之下。 2 1 转子裂纹形成机理 金属弹塑性理论表明零部件热应力与内部温度梯度成正比，交变的温度场引起 交变的应力场，循环周期取决于机组启停或负荷变化过程时间，相对于振动等高周 波机械应力，成为低周应力。应力或应变反复作用使得材料性能发生变化，以致出 现裂纹。疲劳裂纹经历一个形成和发展过程，并以寿命损耗衡量。致裂寿命即裂纹 形成寿命，是指零件从初次投运到出现第一条宏观裂纹所经历的应力或应变循环 数。出现裂纹过程中，循环应变不同，会产生循环硬化、循环软化、循环稳定以及 硬化软化混合型n 们。 循环硬化：指在恒定的应变循环下，。循环所需的应力不断增加的现象。材料硬 度增加，屈服点和抗拉强度也提高。循环硬化是材料晶体结构中原有的和新生的晶 粒位错交互作用的结果，在有规则的原子晶格中，这种破坏是应力集中所致，应力 集中起着增加局部应力并超过剪切极限的作用。 循环软化：是指在恒定的应变循环下，循环实验所需的应力不断降低的现象。 材料的硬度、屈服点和抗拉强度也都降低了。循环软化源于微裂纹的成核和萌生。 应变能的循环输入改变了晶格中热量和内表面的能量，这样位错就趋向于集中到已 经包含有一个或多个位错的平面内，不断的输入循环，在这平面上就会聚集更多的 位错，促使在表面产生微观裂纹。 循环稳定：本阶段与微裂纹的缓慢扩展相对应，持续时间长，是评价材料抗疲 劳性能的主要阶段，约占总疲劳寿命的7 0 。 循环硬化和软化兼有的混合型：金属材料被硬化或软化过程在初始阶段发展较 快，并迅速达到饱和，一般在总疲劳寿命的2 0 一2 5 时达到稳定状态。经过循环 稳定阶段后若仍不断的循环或继续输入应变能，将使微观裂纹邻近的晶体变形和歪 曲，最终引起裂纹增长，合并成一个张开面，出现宏观裂纹，如果没有足够的面积 去支撑则将断裂。 华北电力大学硕士学位论文 汽轮机转子工作在高温高压的蒸汽环境，同时自身高速旋转，工作环境恶劣， 受力一般为复合应力。所承受的力和力矩有：叶片、围带的离心力引起的切向和径 向应力，叶片前后的压差引起的轴向推力，机组启停和负荷变化时蒸汽温度梯度引 起的热应力。从转子的受力状态看，转子的结构应力和稳定运行时产生的热应力都 很小，在裂纹形成中扮演较小的作用。汽轮机转子发生裂纹通常有两种情况：一是 机组在启停过程中交变热应力导致的低周疲劳裂纹，多发生在转子表面的应力集中 部位，如调节级前的叶轮根部，高中压转子的弹性槽等部位，热应力的大小取决于 温度的变化量和温升率，温升过快或者运行工况不稳定均能引起热应力的剧烈变 化。温度下降造成表面拉应力对寿命影响较大，是裂纹形成的主要原因。另一种是 转子制造加工时的固有缺陷，如原有微型裂纹、非金属夹杂物( 如硫化氢、碳氢化 合物) 等，通常存在于中心孔表面或靠近中心孔的转子金属内部，为裂纹形成的疲 劳源，在转子承受交变应力时会较快发展成为疲劳裂纹。 转子疲劳断裂包括裂纹的起始、扩展和断裂三个阶段。起始阶段包括微裂纹的 成核和扩展至可检验的工程裂纹两个阶段。从微观上看，金属由晶粒组成，由于晶 粒的大小及取向不同，或存在微观缺陷，如金属夹渣、刀痕或介质腐蚀等，每颗晶 粒强度不同，在交变热应力下( 高温下伴随着蠕变) 应力集中处的材料经受硬化和 软化过程，应力、应变响应变得十分复杂。裂纹的产生一般在微观级由局部塑性区 开始，形成微裂纹而后扩展成可检验的宏观微小裂纹。 低周疲劳的应力、应交响应过程复杂，转子承受的低周疲劳损伤是在应力集中 处循环塑性应变所造成的破坏。其特点是交变应力高，加载频率低，有效循环周次 1 0 4 1 0 5 ；材料经硬化、软化过程，每次循环均产生一定的塑性应变。应变损伤随 循环周次的增加按m i n e r 疲劳损伤线性法则累计： 办。革瓷1 ， 式中：一转子材料在循环应力幅仉或应变幅磊作用下的实际循环数； m 一在循环应力幅珥或应变幅岛作用下的材料失效循环数； 妒，一疲劳寿命损耗的份额。 2 2 低周疲劳寿命损耗的计算 2 2 1 m a n s o n c o f f i n 公式 转子钢的低周疲劳曲线一般用试棒在恒温及恒应变的条件下，通过对称拉压实 验得出的，未考虑载荷不对称性和因长期运行引起的蠕变和松弛。总应变幅与失效 循环数的关系，按照m a n s o n c o f f i n 公式，方程分解为弹性、塑性量部分，其函数 7 华北电力大学硕士学位论文 关系式为啪1 ： 堕2 一生2 + 堡2 ；生e ( 2 广+ g ；( 2 ) c 、7 j 、 7 式中： 一总应变范围； ，一弹性应变范围。 。一塑性应变范围； 仃：一疲劳强度系数； ：一疲劳塑性系数； 6 一疲劳强度指数； c 一疲劳塑性指数； 一裂纹周次。 对国产转子钢3 0 c r l m o l v ，不同温度下相关的参数如表2 - 1 所示。 表2 - 1 低周疲劳特性计算参数 ( 2 - 2 ) 温度o r l f | e b ， ， f c 室温 0 0 0 5 4 1 3 2 3一o 0 8 0 3 5 0 61 8 2 2 9 3 6加8 2 4 8 0 6 6 2 8 0o 0 0 2 9 9 9 0 3加0 2 1 8 2 l i 50 1 3 3 1 1 0旬5 7 3 2 9 8 4 0 0o 0 0 7 9 9 9 8 9_ 0 0 1 3 2 3 1 1 3o 1 1 0 6 0 7 8旬5 1 7 4 2 1 1 4 8 00 0 0 3 9 4 3 8 9旬0 6 1 0 5 80 9 4 8 7 6 8- 0 7 6 8 6 0 2 7 5 1 0o 0 【1 4 6 5 5 9 2司0 8 5 5 9 7 2 o 5 6 3 4 7 5_ 0 7 1 6 4 1 2 4 5 3 80 0 0 3 7 4 7 6 8 0 0 6 9 7 3 7 3o 8 1 7 5 7 6 旬7 5 5 3 1 0 4 5 6 50 0 0 3 3 0 5 0 8旬0 5 6 8 6 7 50 3 2 2 5 1 3 8- 0 6 5 5 5 3 9 2 2 2 全应变的计算 全应变的计算式为乜引： 式中： 童- 考+ 一篆+ 尝声2222 e 、放7 ( 2 - 3 ) 华北电力大学硕士学位论文 k 一循环强度系数； 万一循环应变强化指数。 不同温度下的参数如表( 2 - 2 ) 所示： 表2 - 2 应变计算参数 温度kn 室温 7 4 9 3 9 9 5o 0 5 0 4 9 3 1 7 2 8 07 3 9 2 9 0 60 0 5 4 6 9 2 0 4 4 0 08 1 8 8 5 4 6o ；6 7 0 6 0 4 4 8 06 0 3 8 1 3 50 0 6 1 7 5 8 0 3 5 1 0 5 9 2 7 2 1 9o 0 6 4 7 5 8 7 5 3 8 4 9 4 6 8 7o 0 4 5 7 8 2 8 3 5 6 54 8 6 8 0 跎0 0 4 9 4 7 5 6 2 2 30 0 - i ：t ：i n 频率修正式 c o f f i n 在考虑应力保持时间和蠕变疲劳交互作用的影响后提出频率修正式。引 入交变频率v ，则材料的塑性范围可表示为口制： 、 g 。= c n 一( h ) 声( 2 - 4 ) 式中： 芦，c 一材料常数； 七一频率作用常数； ，一循环频率。 式( 2 - 3 ) 中如果令k = 0 ，则g ，- c 畴) p 为纯蠕变，而当k = 1 时频率对寿命 没有影响，为常温疲劳。 弹性应变范围可用塑性应变范围函数的形式给出： 他一4 ( 每厂矿( 2 - 5 ) 式中 4 ，刀一材料应变硬化特征常数； 七一频率修正指数。 将式( 2 - 3 ) 和式( 2 - 4 ) 代入式( 2 - 2 ) 可得： a e , 一a n 4 1 厶+ 删d ( 2 6 ) 华北电力大学硕士学位论文 不同国家的实验结果如表2 3 所示： 表2 - 3 频率修正法实验数据 实验温 ababcd 度( ) 西屋公司 o 0 0 9 4 0 8 8 5 一o 0 9 2 o 0 3 3m 7 5 9o 0 3 45 3 8 ( 美) 路易斯研 究中心0 0 0 9 72 8 0加0 9 50 0 8 00 8 3 1o 1 0 25 3 8 ( 美) o 5 21 1 8m 0 8 8o o 1 0 4_ o 7 3o 吣85 三菱公司 ( 日) o 0 0 6 01 3 60 0 7 7o 0 0 7 7坩7 7o 0 7 75 6 0 由于上表中的数据均为试件断裂时得到的结果，式( 2 - 5 ) 中n 为致断循环周 次，因此得到的寿命为转子的全寿命，如果要求出转子的频率寿命，当初始裂纹深 度定义为0 1 3 0 3 8 m m 时，至裂循环周次为： m = n 一4 嘶( 2 7 ) 式中： m 一致裂循环周次( 疲劳寿命) ； 一致裂循环周次( 全寿命) 。 2 3 蠕变寿命损耗的计算 在机组运行过程中，转子除了承受交变应力产生低周疲劳外，在高温条件下， 钢材即使处于应力比其屈服极限低的情况下，随承载时间的增加而产生缓慢增加的 塑性变形以致断裂，称为塑性损伤。在高温条件下，即使载荷不大，钢材的寿命也 有一定的限度。因此在估算转子寿命损耗时，应考虑转子在带负荷稳定运行时高温 蠕变损伤所占百分比。 断裂力学和热弹塑性力学理论提供了估算转子裂纹扩展速度的基础。裂纹扩展 速度的决定因素是材料强度、交变应力的水平以及裂纹或缺陷的形状。为确保火电 设备的可靠性，必须对转子的老化程度进行确切的评价，但实际中很难从转子上切 割实验样品用来分析，因此转子的蠕变疲劳损伤难以得到。目前广泛采用非破坏性 华北电力大学硕士学位论文 方法来评价转子的蠕变疲劳损伤。 通过在实验室实验测得钢材取样在一定温度和一定应力下的断裂时间，可得到 转子蠕变损耗曲线。实际上。材料蠕变的破坏过程是非线性的，但有线性叠加求得 的蠕变损伤在数值上基本和非线性损伤相近。 考虑到汽轮机部件工作温度及承受的应力是变化的，根据r o b b i n s o n 损伤积累 法则，转子的蠕变损伤指数为口引： 唬；争 。 ( 2 8 ) 式中： z 一相应与一定应力、温度、材料发生蠕变断裂时间； t 一材料实际承受应力、温度的时间。 式( 2 - 5 ) 中z 由l a r s o n m i n e r 参数p 计算： p = o + 2 7 3 ) ( 2 0 + l o g 互) 式中： f 一运行温度，单位； z 一蠕变断裂时间，单位小时； 对于转子钢3 0 c r l m o v ，对应于一定的蠕变断裂应力仃， 与维氏硬度存在以下关系： p 一么( 盯) 月y + b ( 盯) ( 2 - 9 ) l a r s o n m i n e r 参数p ( 2 - 1 0 ) 彳( 歹) 一7 1 4 2 - 9 1 4 3 奉1 0 9 1 0 ( 仃) + 3 1 1 8 事1 9 9 1 0 ( 仃) 宰l o g l o ( 口) ( 2 1 1 ) 曰p ) 一- 9 1 5 3 9 + 3 8 2 7 5 0 宰1 0 曙l o ( o ) - 1 2 4 3 8 8 事l o g l o p ) 宰l 0 9 1 0 p ) ( 2 1 2 ) 由式( 2 6 ) 和( 2 - 7 ) 可以看出，工作应力一定时，随着转子硬度的降低，参 数p 值降低，相应的z = 减少，也就是说断裂时间缩短，蠕变损耗增加。所以说硬度 降低导致蠕变损耗增加，断裂时间比用新转子( 初始硬度) 预计的断裂时间缩短了， 也就是说单位时间的蠕变损耗是随着硬度降低而逐渐增加的变量。 华北电力大学硕士学位论文 第三章汽轮机启动的研究 3 16 0 0 m w 汽轮机转子寿命分配 汽轮机在启停和变负荷运行时，应控制蒸汽温度和负荷的变化率，控制汽轮机 部件的热应力，使机组的寿命损耗不超过其预分配值，在机组规定的使用年限内， 实现最佳的安全经济运行，使机组发挥最佳的经济效益，实现机组运行寿命的科学 管理。 目前通常认为汽轮机的服役年限为3 0 年，在这3 0 年的时间内，如何合理分配 汽轮机的寿命，以取得最大的经济效益是汽轮机寿命分配的出发点。 汽轮机寿命分配，取决于机组在电网中的运行方式。对于带基本负荷的机组， 汽轮机寿命的损耗主要为高温蠕变和正常检修启停所需低周疲劳对汽轮机寿命的 损耗。若年平均运行小时以7 0 0 0 h 计算，3 0 年内预计蠕变寿命损耗以及考虑不定因 素( 如负荷、蒸汽参数的波动、事故带厂用电运行等) 的损耗约占寿命的2 5 ，剩 余7 5 的寿命可分配给汽轮机的启停使用。带基本负荷的机组，终身启停次数较少。 因此，每次启停的寿命的损耗可以分配的较大，可选用较高的温升率启动和较快的 冷却法停机，以缩短机组启停过程的时间，增加机组运行时间，多发电。 对于调峰机组，除检修、维护需正常启停以外还应根据电网要求，安排一定次 数的热状态启动和一定范围的负荷变化。在分配寿命时，既要考虑汽轮机寿命的合 理损耗，又要考虑到电网的调峰需要。 我国进口大容量汽轮机技术谈判指南中规定，汽轮机的使用寿命一般为3 0 年， 在这3 0 年中，机组冷态启动1 0 0 次，温态启动7 0 0 次，热态启动3 0 0 0 次，极热态 启动1 5 0 次以及大于1 0 额定负荷突变1 2 0 0 次，其总寿命消耗不大于7 5 。 以下是不同公司6 0 0 m w 汽轮机的寿命分配表呻1 ： 表3 - 1a b b 6 0 0 m w 寿命分配 每年寿命损耗3 0 年内寿命损 运行方式寿命损耗率( ) 3 0 年使用次数 ( )耗( ) 冷态启动0 0 0 5 1 0 00 0 1 6 70 5 温态启动o 0 0 31 2 0 00 1 23 6 热态启动 0 0 0 1 34 5 0 00 1 9 55 8 5 极热态启动 o 0 0 15 0 00 0 1 6 70 5 负荷突变 0 0 0 0 2 54 5 0 00 0 3 7 51 1 2 5 合计 0 3 8 5 91 1 5 7 5 华北电力大学硕士学位论文 表3 - 2 ( 3 a6 0 0 m w 寿命分配 每年寿命损耗3 0 年内寿命损 运行方式寿命损耗率( )3 0 年使用次数 ( ) 耗( ) 冷态启动 0 0 0 11 2 0o 0 0 40 1 2 温态启动 0 0 0 21 1 0 0 0 0 7 3 3 2 2 热态启动 0 0 1 34 5 0 01 9 45 8 5 极热态启动 o 0 0 15 0 00 0 1 6 70 5 负荷突变 o 0 0 14 0 0 00 1 3 3 34 0 合计 6 5 3 2 表3 3 三菱6 0 0 m w 寿命分配 每年寿命损耗3 0 年内寿命损 运行方式 寿命损耗率( )3 0 年使用次数 ( )耗( ) 冷态启动 0 0 11 0 0o 0 31 0 温态启动o 0 1 1 2 0 0o 41 2 0 热态启动 0 0 0 9 14 5 0 01 3 6 54 0 9 5 极热态启动0 0 3 1 5 0 o 1 55 5 2 5 1 0 0 7 5 o 0 0 13 0 0 00 13 o 3 m i n 2 5 1 0 0 5 0 o 0 0 13 0 0 00 13 o 3 r a i n 8 0 1 0 0 5 0 0 0 0 13 0 0 0o 13 o 3 r a i n ，。 5 0 1 0 0 - 5 0 o 0 51 0 0o 1 6 75 0 5 m i n 3 0 - 5 0 2 0 0 0 0 13 0 o 13 o 3 r a i n 合计7 5 4 5 本文计算对象为哈汽6 0 0 m w 的汽轮机，是与三菱公司联合设计生产的，所以 在研究该汽轮机的启动时，采用三菱公司的寿命分配方式。 3 2 启动应力的计算 计算转子不稳定温度场时，可以认为转子是均匀的、各向同性且无内热源的物 体，一般来说，求解转子温度场可采用轴对称的有限元法，但这需要较大的计算机 内存、较长的时间和较烦琐的数据准备工作。根据国外相关研究机构的研究结果表 明，沿转子半径温差的计算偏差最大不超过8 。在现场实现汽轮机自动启停并进 行汽轮机的寿命管理时，要求热应力和寿命损耗计算准确快速可靠。因此可以认为， 温度场问题的一维解法适用于寿命管理，其计算误差是可以允许的口们。 1 3 华北电力大学硕士学位论文 启动时的应力主要是由温度场的不均匀造成的热应力和由叶片、叶轮以及转子 自重产生的离心力所引起的切向应力合成的。热应力分为两部分，一是热冲击造成 的热应力，二是蒸汽温升造成的热应力，总的热应力为二者的叠加。 3 2 1 热冲击应力的计算 设冲转时转子温度为气，蒸汽温度为f ，换热系数为口，将热冲击视为阶跃过程， 则仅由热冲击产生的热应力和通过如下方法求解，不考虑冲转后的汽温变化： o t ，o z t 1otv- 、 万剐( i o r + 7 8 r )a z 、 2 ， 7 f l f - 0 = t o 要k 。一导( f t s ) 石h i i j f 一 可得单独由于( f j - t o ) 的热冲击导致转子温度的变化为： ，。，。+。，一，。)：。12上；薹i聪， 据上式求得由于热( 冷) 冲击在转子外表面产生的热应力为： ( 3 - 1 ) ( 3 - 2 ) 喇。警凄南一墚涨， 净3 ) 上式可表示为q = a a t e 所的形式。4 由材料的尺寸及物性参数决定，r 为蒸汽 与金属表面的温差，口为应力的衰减系数，7 为时间。由此可见，影响冲击应力的 因素为蒸汽与金属的温差以及材料的温度。 本文对6 0 0 m w 汽轮机启动时的热冲击应力进行了计算，转子不同的初温时彳和 曰的值如表( 3 4 ) 所示： 华北电力大学硕士学位论文 表3 - 4 不同转子初温彳、曰值 转子初温( )4 口 1 0 0旬5 6 7 81 - 7 8 9 3 2 0 0- 0 5 8 8 61 6 3 5 4 3 0 0- 0 5 9 8 61 3 9 7 9 4 0 0 - 0 6 1 4 61 2 8 6 7 假定冲转温差a t = 5 0 一定，分别给定转子初温1 0 0 、2 0 0 、3 0 0 、4 0 0 进行计算，可以看到冲转时转子温度对产生冲击应力的影响( 图3 - 1 ) 。 o - 5 乍0 罡l s 葶 餐* 静 羹籀 徽 3 5 矗鑫51 霉1 $2 矗 2 _ 嚣 瓣衡讯， 图3 - 1 冲击温差一定时，转子表面冲击应力与转子初温的关系 由上图可以看出，对于相同的冲击温差( 5 0 ) ，最大的应力差值为5 m p a 左右， 材料初始温度对热冲击应力的影响不是很大，所以冲击应力主要取决于温差。 3 2 2 蒸汽温升造成的热应力 将转子视为无限长圆柱体，忽略轴向和切向的温度分布不均匀对其的影响，则 截面内温度和热应力的分布是沿轴向对称的，并随半径的改变而变化，这样的问题 就被简化为一维轴对称瞬态热传导问题，并与初始汽温一致时，若汽温随时间呈线 性变化，则解析法求解方程( 3 1 ) ，可得转子任意半径处温度随时间的变化式： f = 岛何一誓( 1 一再r 2 + p 2