CN116384191A - 用于评估汽轮机叶片水蚀缺损安全性的有限元计算方法 - Google Patents
用于评估汽轮机叶片水蚀缺损安全性的有限元计算方法 Download PDFInfo
- Publication number
- CN116384191A CN116384191A CN202310366634.8A CN202310366634A CN116384191A CN 116384191 A CN116384191 A CN 116384191A CN 202310366634 A CN202310366634 A CN 202310366634A CN 116384191 A CN116384191 A CN 116384191A
- Authority
- CN
- China
- Prior art keywords
- blade
- defect
- water erosion
- erosion defect
- stress
- 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
Links
- 230000003628 erosive effect Effects 0.000 title claims abstract description 174
- 230000007547 defect Effects 0.000 title claims abstract description 171
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 title claims abstract description 138
- 238000004364 calculation method Methods 0.000 title claims abstract description 34
- 238000000034 method Methods 0.000 claims abstract description 18
- 230000003068 static effect Effects 0.000 claims description 44
- 239000012530 fluid Substances 0.000 claims description 29
- 239000000463 material Substances 0.000 claims description 29
- 230000004044 response Effects 0.000 claims description 11
- 230000008878 coupling Effects 0.000 claims description 6
- 238000010168 coupling process Methods 0.000 claims description 6
- 238000005859 coupling reaction Methods 0.000 claims description 6
- 239000007787 solid Substances 0.000 claims description 6
- 238000010586 diagram Methods 0.000 claims description 4
- 238000005452 bending Methods 0.000 claims description 3
- 125000004122 cyclic group Chemical group 0.000 claims description 3
- 238000013016 damping Methods 0.000 claims description 3
- 230000005284 excitation Effects 0.000 claims description 3
- 238000012795 verification Methods 0.000 claims description 3
- 238000012423 maintenance Methods 0.000 abstract description 6
- 238000013461 design Methods 0.000 description 3
- 230000002950 deficient Effects 0.000 description 2
- 230000007774 longterm Effects 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000008859 change Effects 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000011156 evaluation Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000008569 process Effects 0.000 description 1
- 230000009467 reduction Effects 0.000 description 1
- 238000012502 risk assessment Methods 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
- 230000007704 transition Effects 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
- G06F30/23—Design optimisation, verification or simulation using finite element methods [FEM] or finite difference methods [FDM]
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/10—Geometric CAD
- G06F30/17—Mechanical parametric or variational design
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
- G06F30/28—Design optimisation, verification or simulation using fluid dynamics, e.g. using Navier-Stokes equations or computational fluid dynamics [CFD]
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T17/00—Three dimensional [3D] modelling, e.g. data description of 3D objects
- G06T17/20—Finite element generation, e.g. wire-frame surface description, tesselation
- G06T17/205—Re-meshing
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2113/00—Details relating to the application field
- G06F2113/08—Fluids
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2119/00—Details relating to the type or aim of the analysis or the optimisation
- G06F2119/14—Force analysis or force optimisation, e.g. static or dynamic forces
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Geometry (AREA)
- General Physics & Mathematics (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- General Engineering & Computer Science (AREA)
- Mathematical Optimization (AREA)
- Mathematical Analysis (AREA)
- Pure & Applied Mathematics (AREA)
- Computing Systems (AREA)
- Fluid Mechanics (AREA)
- Algebra (AREA)
- Mathematical Physics (AREA)
- Computational Mathematics (AREA)
- Computer Graphics (AREA)
- Software Systems (AREA)
- Testing Of Devices, Machine Parts, Or Other Structures Thereof (AREA)
Abstract
用于评估汽轮机叶片水蚀缺损安全性的有限元计算方法,属于汽轮机叶片技术领域。本发明解决了现有技术中汽轮机叶片水蚀缺损极限无法量化的问题。技术要点:步骤1,确定叶片水蚀缺损区域;步骤2,建立叶片水蚀缺损模型;步骤3,对叶片水蚀缺损模型进行有限元分析,确定叶片水蚀缺损尺寸上限;步骤4,评估叶片水蚀缺损结构的安全性,确定叶片模型缺损深度和长度的上限尺寸。本发明采用叶片三维气动流场分析,以及叶片现场运行记录,建立叶片水蚀缺损模型,并基于三维有限元方法分析叶片的强度与振动安全性,量化给出叶片的水蚀缺损极限,判断了目标叶片的安全运行风险,对于叶片的安全运行、精密检修具有重要的指导意义。
Description
技术领域
本发明涉及一种叶片水蚀缺损安全性的计算方法,具体涉及一种用于评估汽轮机叶片水蚀缺损安全性的有限元计算方法,属于汽轮机叶片技术领域。
背景技术
汽轮机在运行时不仅受到离心力的作用,还要受到水蒸汽的冲击,在汽轮机低压通流中,从进口到出口随着水蚀汽温度和压力参数的不断降低,水蒸汽会出现由单相向两相转换的过程,这使得通流中部分水蒸汽变成了水滴,水滴夹杂在高速的水蒸汽中,不断冲击着叶片,在高速水滴的冲蚀下叶片缺损不断变大,直至叶片出现断裂失效。由于汽轮机的工质为水蒸汽,导致叶片的水蚀是不可避免的,在汽轮机机组大修时,对于叶片水蚀缺损的安全性没有一个量化的指标,在多数情况下,叶片能否继续运行完全依靠工程经验,这使得叶片在运行过程中存在安全性风险。因此,亟待开发一种用于评估汽轮机叶片水蚀缺损安全性的有限元计算方法,判断目标叶片的安全运行风险,用以解决叶片水蚀缺损极限无法量化的问题。
发明内容
本发明为了克服上述现有技术中汽轮机叶片水蚀缺损极限无法量化的问题,提供了一种用于评估汽轮机叶片水蚀缺损安全性的有限元计算方法,在下文中给出了关于本发明的简要概述,以便提供关于本发明的某些方面的基本理解。应当理解,这个概述并不是关于本发明的穷举性概述。它并不是意图确定本发明的关键或重要部分,也不是意图限定本发明的范围。
本发明的技术方案:
一种用于评估汽轮机叶片水蚀缺损安全性的有限元计算方法,包括以下步骤:
步骤1,确定叶片水蚀缺损区域
进行叶片全三维气动流场分析确定,或者根据叶片在现场的运行情况确定水蚀区域,记录机组每次大修时确定叶片的水蚀缺损区域,将缺损严重的区域确定为水蚀缺损区域;
步骤2,建立叶片水蚀缺损模型
根据现场的运行情况以及全三维气动流场分析结果建立叶片水蚀缺损模型,在叶片水蚀缺损区域取不同深度、不同长度建立组合模型,模型缺损深度和长度的上限尺寸为叶片的强度与振动的许用值;
步骤3,对叶片水蚀缺损模型进行有限元分析,确定叶片水蚀缺损尺寸上限
有限元分析包括静应力、动频率以及动应力分析;
其中:静应力基于有限元弹性分析,分析考核叶片的危险界面的平均应力和局部的峰值应力;
动频率基于运行转速下的预应力,分析考核叶片的共振转速是否避开危险运行区间;
动应力基于谐响应分析,将流体计算得到的压力作为激振力,加载至叶片结构表面进行响应分析,分析考核叶片的动应力是否满足叶片材料的耐振强度极限;
步骤4,评估叶片水蚀缺损结构的安全性,确定叶片模型缺损深度和长度的上限尺寸
根据步骤3中的有限元分析结果,得到的叶片静应力水蚀缺损极限,动频率水蚀缺损极限,以及动应力水蚀缺损极限,取三个极限中叶片水蚀缺损尺寸最小值,作为判断目标叶片是否超出水蚀缺损的尺寸极限,如果叶片水蚀缺损超出此尺寸极限,叶片将存在失效风险;如果叶片水蚀缺损小于此尺寸极限,叶片可继续运行。
进一步地:所述步骤1中,进行叶片全三维气动流场分析,得到叶片在不同工况下的流线和湿度,结合电厂实际运行后的叶片水蚀缺损情况,确定叶片的水蚀区域;具体地:
首先,进行叶片的流体计算几何模型建立,模型采用循环周期对称模型,几何参数包括静、动叶片的中径、叶高、根部扩张角、根部扩张角、叶片数量以及叶片转速;
然后,基于几何模型进行流体网格划分,网格采用结构化六面体网格,并进行网格无关性验证,取计算效率和精度最优的网格数量;
然后,对流体网格模型加载边界参数,流体状态为粘性、可压缩流动,流体满足物理守恒定律;流体采用水和水蒸气工质,采用SST湍流模型,壁面采用光滑和绝热的边界条件,静叶流域设置为静止域,动叶流域设置为旋速3000rpm,计算流域进口给定总压、总焓,出口给定静压,得到叶片在不同工况下的流线和湿度;
最后,结合电厂实际运行后的叶片水蚀缺损情况,确定叶片的水蚀区域,当叶片在41%以上负荷工况时,水蒸汽主要冲蚀叶片的进汽侧,且由于叶型顶部的气流速度比叶型根部大,选择叶片的叶型顶部进汽侧为水蚀缺损区域;当叶片在41%以下负荷工况时,由于流道中出现了涡流、回流现象,导致叶型根部出汽侧出现水蚀,选择叶片叶型根部出汽侧为水蚀缺损区域。
进一步地:所述步骤2中,在同一水蚀区域建立不同高度和深度的水蚀缺损模型,水蚀在高度方向缺损模型,每50mm建立一个模型;水蚀在深度方向缺损模型,每2mm建立一个模型。
进一步地:所述步骤3中,对叶片水蚀缺损模型进行有限元分析,有限元计算模型采用循环对称设置,循环对称体包括一个节距的叶片和转子轮槽,对于叶片的围带与围带之间、拉筋与拉筋之间、叶根与叶根之间、叶根与轮槽之间采用接触设置,对于叶片和转子设置材料的密度、弹性模量、泊松比;叶片和转子网格采用结构化的六面体网格,网格边长比控制在3:1,采用一阶减缩积分实体单元;
叶片静载荷包括离心力和蒸汽压力,离心力设置叶片和转子的旋转角速度,蒸汽压力基于流体分析结果,经流固耦合后加载至叶片表面,进行静应力分析;
叶片动频率计算在2000rpm~3300rpm转速之间,每隔100rpm进行动频率分析,并在坎贝尔图中绘制共振转速和共振频率;
基于叶片流体计算结果,经流固耦合给定叶片激振力,并根据叶片的阶数和节径数给定激振力的激振因子,结合叶片的结构类型,给定阻尼系数,采用谐响应方法进行叶片的动应力分析。
进一步地:静应力分析中,当叶片水蚀缺损模型的峰值应力超出叶片材料的许用值时,此时的叶片水蚀缺损结构即为静强度考核下的上限尺寸,对于叶片的静应力考核公式如下:
峰值应力σF=(σ2+3τ2)1/2 (1)
[σ]= Ki / Si (2)
峰值应力σF<[σ](3)
式中:σ为主应力;τ为剪应力;Ki为材料的屈服强度;Si为安全系数;σF为叶片峰值应力;[σ]为叶片的静强度许用值。
进一步地:动频率分析中,当叶片水蚀缺损模型的动频率超出允许的安全范围时(工作频率*6%~工作频率*3%),此时的叶片水蚀缺损结构即为动频率考核下的上限尺寸,对于叶片的动频率考核公式如下:
工作频率*6%<f01 or f02<工作频率*3%(6)
式中:f01为弯曲振动频率;f02为扭转振动频率;L为叶高;A为截面面积;ρ为叶片材料密度;E为材料弹性模量;G为剪切模量;T0为扭转刚度几何因子;I0为叶片截面极惯性矩。
进一步地:动应力分析中,当叶片水蚀缺损模型的动应力超出叶片材料的许用值时,此时的叶片水蚀缺损结构即为动应力考核下的上限尺寸,对于叶片的动应力考核公式如下:
[σa]=σ-1 (8)
本发明有益效果体现在:
相对于现有技术,本发明采用叶片三维气动流场分析,以及叶片现场运行记录,建立叶片水蚀缺损模型,并基于三维有限元方法分析叶片的强度与振动安全性,量化给出叶片的水蚀缺损极限,判断了目标叶片的安全运行风险,解决了叶片水蚀缺损极限无法量化给出的问题。采用本发明可制定叶片的长期维修策略,有效的指导电厂检修人员高效检修,避免了人为因素引起的误判,大大降低了因叶片水蚀缺损引起的运行安全性风险。本发明操作简单,通用性强,能够普遍应用于汽轮机叶片水蚀缺损的安全性风险评估。本发明对于叶片的安全运行、精密检修具有重要的指导意义。
附图说明
图1是叶片水蚀缺损安全性有限元分析评估流程图;
图2是汽轮机低压通流水蒸汽相变示意图;
图3是叶片在不同负荷工况下的流体变化趋势;
图4是气动流场分析结果(叶片湿度云图);
图5是叶片运行现场水蚀情况,a为叶片进汽侧水蚀,b为叶片出汽侧水蚀;
图6是叶片常见的水蚀缺损示意图。
具体实施方式
为使本发明的目的、技术方案和优点更加清楚明了,下面通过附图中示出的具体实施例来描述本发明。但是应该理解,这些描述只是示例性的,而并非要限制本发明的范围。此外,在以下说明中,省略了对公知结构和技术的描述,以避免不必要地混淆本发明的概念。
实施例1,结合图1-图6说明本实施例,本实施例提供了一种用于评估汽轮机叶片水蚀缺损安全性的有限元计算方法,包括以下步骤:步骤1,确定叶片水蚀缺损区域,进行叶片全三维气动流场分析确定,或者根据叶片在现场的运行情况确定水蚀区域,记录机组每次大修时确定叶片的水蚀缺损区域,将缺损严重的区域确定为水蚀缺损区域;步骤2,建立叶片水蚀缺损模型,根据现场的运行情况以及全三维气动流场分析结果建立叶片水蚀缺损模型,在叶片水蚀缺损区域取不同深度、不同长度建立组合模型,模型缺损深度和长度的上限尺寸为叶片的强度与振动的许用值;步骤3,对叶片水蚀缺损模型进行有限元分析,确定叶片水蚀缺损尺寸上限,有限元分析包括静应力、动频率以及动应力分析;其中:静应力基于有限元弹性分析,分析考核叶片的危险界面的平均应力和局部的峰值应力;动频率基于运行转速下的预应力,分析考核叶片的共振转速是否避开危险运行区间;动应力基于谐响应分析,将流体计算得到的压力作为激振力,加载至叶片结构表面进行响应分析,分析考核叶片的动应力是否满足叶片材料的耐振强度极限;步骤4,评估叶片水蚀缺损结构的安全性,确定叶片模型缺损深度和长度的上限尺寸,根据步骤3中的有限元分析结果,得到的叶片静应力水蚀缺损极限,动频率水蚀缺损极限,以及动应力水蚀缺损极限,取三个极限中叶片水蚀缺损尺寸最小值,作为判断目标叶片是否超出水蚀缺损的尺寸极限,如果叶片水蚀缺损超出此尺寸极限,叶片将存在失效风险;如果叶片水蚀缺损小于此尺寸极限,叶片可继续运行。
实施例2,结合图1-图6说明本实施例,本实施例提供了一种用于评估汽轮机叶片水蚀缺损安全性的有限元计算方法,具体地:
步骤1,确定叶片水蚀缺损区域。叶片的水蚀缺损区域主要通过两种途径来确定。
一是进行叶片全三维气动流场分析,得到叶片在不同工况下的流线(如图3)和湿度(如图4),结合电厂实际运行后的叶片水蚀缺损情况,确定叶片的水蚀区域;具体地:
首先,进行叶片的流体计算几何模型建立,模型采用循环周期对称模型,几何参数包括静、动叶片的中径、叶高、根部扩张角、根部扩张角、叶片数量以及叶片转速;
然后,基于几何模型进行流体网格划分,网格采用结构化六面体网格,并进行网格无关性验证,取计算效率和精度最优的网格数量;
然后,对流体网格模型加载边界参数,流体状态为粘性、可压缩流动,流体满足物理守恒定律;流体采用水和水蒸气工质,采用SST湍流模型,壁面采用光滑和绝热的边界条件,静叶流域设置为静止域,动叶流域设置为旋速3000rpm,计算流域进口给定总压、总焓,出口给定静压,得到叶片在不同工况下的流线和湿度;
最后,结合电厂实际运行后的叶片水蚀缺损情况,确定叶片的水蚀区域,根据图3中叶片的流线可以看到,当叶片在41%以上负荷工况时,水蒸汽主要冲蚀叶片的进汽侧,且由于叶型顶部的气流速度比叶型根部大,选择叶片的叶型顶部进汽侧为水蚀缺损区域;当叶片在41%以下负荷工况时,由于流道中出现了涡流、回流现象,导致叶型根部出汽侧出现水蚀,选择叶片叶型根部出汽侧为水蚀缺损区域。
二是根据叶片在现场的运行情况确定水蚀区域,如图5所示,记录机组每次大修时确定叶片的水蚀缺损区域,将缺损较为严重的区域进行重点关注。
步骤2,建立叶片水蚀缺损模型。
根据实际的运行情况以及全三维气动流场分析结果建立叶片水蚀缺损模型。如图6所示,水蚀缺损模型包括以下类型,一是叶片仅有叶型顶部进汽侧为主要的水蚀缺损区域,以原始设计模型为基准模型,在叶型顶部进汽侧区域取不同深度(以深度2mm为缺损增量)、不同长度(以长度50mm为缺损增量)建立组合模型,模型缺损深度和长度的上限尺寸为叶片的强度与振动的许用值。二是叶片仅有叶型根部出汽侧为水蚀缺损区域,以原始设计模型为基准模型,在叶型根部出汽侧取不同深度(以深度2mm为缺损增量)、不同长度(以长度50mm为缺损增量)建立组合模型,即在同一水蚀区域建立不同高度和深度的水蚀缺损模型,水蚀在高度方向缺损模型,每50mm建立一个模型;水蚀在深度方向缺损模型,每2mm建立一个模型。模型缺损深度和长度的上限尺寸为叶片的强度与振动的许用值。三是叶片叶型顶部进汽侧以及叶型根部出汽侧均有缺损时,以原始设计模型为基准模型,取不同深度(以深度2mm为缺损增量)、不同长度(以长度50mm为缺损增量)建立组合模型,模型缺损深度和长度的上限尺寸为叶片的强度与振动的许用值。
步骤3,对叶片水蚀缺损模型进行有限元分析,确定叶片水蚀缺损尺寸上限。
有限元分析包括静应力、动频率以及动应力分析;其中:静应力基于有限元弹性分析,分析考核叶片的危险界面的平均应力和局部的峰值应力;动频率基于运行转速下的预应力,分析考核叶片的共振转速是否避开危险运行区间;动应力基于谐响应分析,将流体计算得到的压力作为激振力,加载至叶片结构表面进行响应分析,分析考核叶片的动应力是否满足叶片材料的耐振强度极限;
对叶片水蚀缺损模型进行有限元分析,有限元计算模型采用循环对称设置,循环对称体包括一个节距的叶片和转子轮槽,对于叶片的围带与围带之间、拉筋与拉筋之间、叶根与叶根之间、叶根与轮槽之间采用接触设置。对于叶片和转子设置材料的密度、弹性模量、泊松比等参数。叶片和转子网格采用结构化的六面体网格,网格边长比控制在3:1,采用一阶减缩积分实体单元。叶片静载荷包括离心力和蒸汽压力,其中离心力设置叶片和转子的旋转角速度,蒸汽压力基于流体分析结果,经流固耦合后加载至叶片表面,进行静力学分析。叶片动频率计算在2000rpm~3300rpm转速之间,每隔500rpm进行动频率分析,并坎贝尔图中绘制共振转速和共振频率。基于叶片流体计算结果,经流固耦合给定叶片激振力,并根据叶片的阶数和节径数给定激振力的激振因子,结合叶片的结构类型,给定阻尼系数,采用谐响应方法进行叶片的动应力分析。
根据叶片水蚀缺损模型进行有限元静强度分析结果,对叶片的有限元峰值应力进行静强度考核,当叶片水蚀缺损模型的峰值应力超出叶片材料的许用值时,此时的叶片水蚀缺损结构即为静强度考核下的上限尺寸,对于叶片的静应力考核公式如下:
峰值应力σF=(σ2+3τ2)1/2 (1)
[σ]= Ki / Si (2)
峰值应力σF<[σ](3)
式中:σ为主应力;τ为剪应力;Ki为材料的屈服强度;Si为安全系数;σF为叶片峰值应力;[σ]为叶片的静强度许用值。
根据叶片水蚀缺损模型进行有限元动频率分析结果,基于叶片坎贝尔图进行动频率考核,当叶片水蚀缺损模型的动频率超出允许的安全范围时(工作频率*6%~工作频率*3%),此时的叶片水蚀缺损结构即为动频率考核下的上限尺寸,对于叶片的动频率考核公式如下:
工作频率*6%<f01 or f02<工作频率*3%(6)
式中:f01为弯曲振动频率;f02为扭转振动频率;L为叶高;A为截面面积;ρ为叶片材料密度;E为材料弹性模量;G为剪切模量;T0为扭转刚度几何因子;I0为叶片截面极惯性矩。
根据叶片水蚀缺损模型进行有限元动应力分析结果,基于叶片材料的耐振强度曲线进行动应力考核,当叶片水蚀缺损模型的动应力超出叶片材料的许用值时,此时的叶片水蚀缺损结构即为动应力考核下的上限尺寸,对于叶片的动应力考核公式如下:
[σa]=σ-1 (8)
步骤4,评估叶片水蚀缺损结构的安全性,确定叶片模型缺损深度和长度的上限尺寸。
根据叶片水蚀缺损模型有限元分析结果,得到的叶片在静强度考核下的水蚀缺损极限尺寸,动频率考核下的水蚀缺损极限尺寸,以及动强度考核下的水蚀缺损极限尺寸,取三种考核下的水蚀缺损极限尺寸下的最小值,即为叶片水蚀缺损极限尺寸,从而可以定两给出叶片的水蚀缺损极限,判定方法如下所示。当叶片水蚀缺损到达3个极限尺寸中最小值时,叶片存在失效风险。当叶片水蚀缺损均小于3个极限尺寸中最小值时,叶片可继续运行。
叶片水蚀缺损极限=min{极限尺寸静应力;极限尺寸动频率;极限尺寸动应力} (10)。
上述实施例是一种用于评估汽轮机叶片水蚀缺损安全性的有限元计算方法,该方法采用数值分析结合实际运行情况,可以定量给出叶片的水蚀缺损极限尺寸。叶片在运行时发生缺损失效,可能会引起整圈叶片的失效,严重的可将汽缸射穿,造成设备损坏、人员伤亡。采用上述实施例可制定叶片的长期维修策略,降低叶片运行风险,避免叶片故障引发的停机事故;操作简单,通用性强,能够普遍应用于具有水蚀缺损的汽轮机叶片。
以上所述的实施例,对本申请的目的、技术方案和有益效果进行了进一步详细说明,所应理解的是,以上所述仅为本申请的实施例而已,并不用于限定本申请的保护范围,凡在本申请的技术方案的基础之上,所做的任何修改、等同替换、改进等,均应包括在本申请的保护范围之内。
Claims (7)
1.用于评估汽轮机叶片水蚀缺损安全性的有限元计算方法,其特征在于,包括以下步骤:
步骤1,确定叶片水蚀缺损区域
进行叶片全三维气动流场分析确定,或者根据叶片在现场的运行情况确定水蚀区域,记录机组每次大修时确定叶片的水蚀缺损区域,将缺损严重的区域确定为水蚀缺损区域;
步骤2,建立叶片水蚀缺损模型
根据现场的运行情况以及全三维气动流场分析结果建立叶片水蚀缺损模型,在叶片水蚀缺损区域取不同深度、不同长度建立组合模型,模型缺损深度和长度的上限尺寸为叶片的强度与振动的许用值;
步骤3,对叶片水蚀缺损模型进行有限元分析,确定叶片水蚀缺损尺寸上限
有限元分析包括静应力、动频率以及动应力分析;
其中:静应力基于有限元弹性分析,分析考核叶片的危险界面的平均应力和局部的峰值应力;
动频率基于运行转速下的预应力,分析考核叶片的共振转速是否避开危险运行区间;
动应力基于谐响应分析,将流体计算得到的压力作为激振力,加载至叶片结构表面进行响应分析,分析考核叶片的动应力是否满足叶片材料的耐振强度极限;
步骤4,评估叶片水蚀缺损结构的安全性,确定叶片模型缺损深度和长度的上限尺寸根据步骤3中的有限元分析结果,得到的叶片静应力水蚀缺损极限,动频率水蚀缺损极限,以及动应力水蚀缺损极限,取三个极限中叶片水蚀缺损尺寸最小值,作为判断目标叶片是否超出水蚀缺损的尺寸极限,如果叶片水蚀缺损超出此尺寸极限,叶片将存在失效风险;如果叶片水蚀缺损小于此尺寸极限,叶片可继续运行。
2.根据权利要求1所述的用于评估汽轮机叶片水蚀缺损安全性的有限元计算方法,其特征在于:所述步骤1中,进行叶片全三维气动流场分析,得到叶片在不同工况下的流线和湿度,结合电厂实际运行后的叶片水蚀缺损情况,确定叶片的水蚀区域;具体地:
首先,进行叶片的流体计算几何模型建立,模型采用循环周期对称模型,几何参数包括静、动叶片的中径、叶高、根部扩张角、根部扩张角、叶片数量以及叶片转速;
然后,基于几何模型进行流体网格划分,网格采用结构化六面体网格,并进行网格无关性验证,取计算效率和精度最优的网格数量;
然后,对流体网格模型加载边界参数,流体状态为粘性、可压缩流动,流体满足物理守恒定律;流体采用水和水蒸气工质,采用SST湍流模型,壁面采用光滑和绝热的边界条件,静叶流域设置为静止域,动叶流域设置为旋速3000rpm,计算流域进口给定总压、总焓,出口给定静压,得到叶片在不同工况下的流线和湿度;
最后,结合电厂实际运行后的叶片水蚀缺损情况,确定叶片的水蚀区域,当叶片在41%以上负荷工况时,水蒸汽主要冲蚀叶片的进汽侧,且由于叶型顶部的气流速度比叶型根部大,选择叶片的叶型顶部进汽侧为水蚀缺损区域;当叶片在41%以下负荷工况时,由于流道中出现了涡流、回流现象,导致叶型根部出汽侧出现水蚀,选择叶片叶型根部出汽侧为水蚀缺损区域。
3.根据权利要求2所述的用于评估汽轮机叶片水蚀缺损安全性的有限元计算方法,其特征在于:所述步骤2中,在同一水蚀区域建立不同高度和深度的水蚀缺损模型,水蚀在高度方向缺损模型,每50mm建立一个模型;水蚀在深度方向缺损模型,每2mm建立一个模型。
4.根据权利要求3所述的用于评估汽轮机叶片水蚀缺损安全性的有限元计算方法,其特征在于:所述步骤3中,对叶片水蚀缺损模型进行有限元分析,有限元计算模型采用循环对称设置,循环对称体包括一个节距的叶片和转子轮槽,对于叶片的围带与围带之间、拉筋与拉筋之间、叶根与叶根之间、叶根与轮槽之间采用接触设置,对于叶片和转子设置材料的密度、弹性模量、泊松比;叶片和转子网格采用结构化的六面体网格,网格边长比控制在3:1,采用一阶减缩积分实体单元;
叶片静载荷包括离心力和蒸汽压力,离心力设置叶片和转子的旋转角速度,蒸汽压力基于流体分析结果,经流固耦合后加载至叶片表面,进行静应力分析;
叶片动频率计算在2000rpm~3300rpm转速之间,每隔100rpm进行动频率分析,并在坎贝尔图中绘制共振转速和共振频率;
基于叶片流体计算结果,经流固耦合给定叶片激振力,并根据叶片的阶数和节径数给定激振力的激振因子,结合叶片的结构类型,给定阻尼系数,采用谐响应方法进行叶片的动应力分析。
5.根据权利要求4所述的用于评估汽轮机叶片水蚀缺损安全性的有限元计算方法,其特征在于:静应力分析中,当叶片水蚀缺损模型的峰值应力超出叶片材料的许用值时,此时的叶片水蚀缺损结构即为静强度考核下的上限尺寸,对于叶片的静应力考核公式如下:
峰值应力σF=(σ2+3τ2)1/2(1)
[σ]=Ki/Si (2)
峰值应力σF<[σ] (3)
式中:σ为主应力;τ为剪应力;Ki为材料的屈服强度;Si为安全系数;σF为叶片峰值应力;[σ]为叶片的静强度许用值。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202310366634.8A CN116384191A (zh) | 2023-04-07 | 2023-04-07 | 用于评估汽轮机叶片水蚀缺损安全性的有限元计算方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202310366634.8A CN116384191A (zh) | 2023-04-07 | 2023-04-07 | 用于评估汽轮机叶片水蚀缺损安全性的有限元计算方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN116384191A true CN116384191A (zh) | 2023-07-04 |
Family
ID=86965216
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202310366634.8A Pending CN116384191A (zh) | 2023-04-07 | 2023-04-07 | 用于评估汽轮机叶片水蚀缺损安全性的有限元计算方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN116384191A (zh) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN116933447A (zh) * | 2023-09-17 | 2023-10-24 | 浙江大学高端装备研究院 | 一种评估带涂层汽轮机叶片阻尼结构的可靠性的方法 |
-
2023
- 2023-04-07 CN CN202310366634.8A patent/CN116384191A/zh active Pending
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN116933447A (zh) * | 2023-09-17 | 2023-10-24 | 浙江大学高端装备研究院 | 一种评估带涂层汽轮机叶片阻尼结构的可靠性的方法 |
CN116933447B (zh) * | 2023-09-17 | 2024-01-09 | 浙江大学高端装备研究院 | 一种评估带涂层汽轮机叶片阻尼结构的可靠性的方法 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN108829917B (zh) | 一种秸秆抛送叶轮疲劳寿命预测方法 | |
CA2893408C (en) | Rotor blade trailing edge | |
CN116384191A (zh) | 用于评估汽轮机叶片水蚀缺损安全性的有限元计算方法 | |
CN108108577B (zh) | 一种水泵叶片疲劳寿命预测及其与导叶最优间距确定方法 | |
CN112861288A (zh) | 一种流体机械校核方法 | |
Shibukawa et al. | A correlation between vibration stresses and flow features of steam turbine long blades in low load conditions | |
CN111444639A (zh) | 一种卧式半调节轴流泵叶轮叶片固定螺栓疲劳寿命预测方法 | |
CN1908382A (zh) | 一种汽轮机转子等效应力安全裕度系数在线监控方法 | |
Algolfat et al. | Dynamic responses analysis of a 5MW NREL wind turbine blade under flap-wise and edge-wise vibrations | |
Lübbe et al. | Design and Validation of a Large Steam Turbine End-Stage Blade to Meet Current and Future Market Demands | |
Veselý et al. | A case study of upgrading of 62.5 MW Pelton turbine | |
Müller et al. | A case study of the fluid structure interaction of a Francis turbine | |
Kaneko et al. | Development and verification of 3000RPM 48inch integral shroud blade for steam turbine | |
Kushner et al. | Critical Review Of Compressor Impeller Vibration Parameters For Failure Prevention. | |
Radgolchin et al. | Investigating the effects of shroud and blade thickness profiles on aeromechanical behavior and fatigue-life of 17-4PH impellers | |
CN112580164B (zh) | 核电汽轮机低压长叶片设计方法 | |
Huang et al. | Effects of tip clearance on aerodynamic damping in a linear turbine cascade | |
Jeanty et al. | Numerical simulation of cavitation phenomena in a centrifugal pump | |
Thakker et al. | Computed effect of guide vane shape on performance of impulse turbine for wave energy conversion | |
CN108915938B (zh) | 一种减轻抽水蓄能电站机组和厂房振动的电站运行方法 | |
Wallin et al. | High speed testing and numerical validation of an aggressive intermediate compressor duct | |
Tao et al. | Detached eddy simulations for high speed axial flow fuel pumps with swept and straight blade impellers | |
Kobayashi et al. | Numerical analysis of stress on pump blade by one-way coupled fluid-structure simulation | |
Tanuma et al. | A Study of Unsteady Aerodynamic Forces on Last Stage Blades in a Low-Pressure Model Steam Turbine at Very Low Operating Conditions | |
Hlous et al. | Development of a New Titanium 54’’LSB: Static and Dynamic Assessment |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination |