CN111814379A - 一种用于核电机组低压末级长叶片的有限元分析方法 - Google Patents
一种用于核电机组低压末级长叶片的有限元分析方法 Download PDFInfo
- Publication number
- CN111814379A CN111814379A CN202010693274.9A CN202010693274A CN111814379A CN 111814379 A CN111814379 A CN 111814379A CN 202010693274 A CN202010693274 A CN 202010693274A CN 111814379 A CN111814379 A CN 111814379A
- Authority
- CN
- China
- Prior art keywords
- low
- long blade
- stage long
- pressure
- last
- 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.)
- Withdrawn
Links
- 238000004458 analytical method Methods 0.000 title claims abstract description 47
- 238000000034 method Methods 0.000 claims abstract description 11
- 230000001105 regulatory effect Effects 0.000 claims abstract description 4
- 230000001419 dependent effect Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
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
- 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)
- Geometry (AREA)
- Theoretical Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Evolutionary Computation (AREA)
- Computer Hardware Design (AREA)
- General Engineering & Computer Science (AREA)
- Pure & Applied Mathematics (AREA)
- Mathematical Optimization (AREA)
- Mathematical Analysis (AREA)
- Computational Mathematics (AREA)
- Turbine Rotor Nozzle Sealing (AREA)
Abstract
一种用于核电机组低压末级长叶片的有限元分析方法,属于核电汽轮机组长叶片分析领域。解决了现有技术中常规汽轮机组的有限元的分析方法处理流程复杂、精度低及无法适用于核电机组的问题。本发明首先构建三维几何模型,对构建的三维几何模型划分网格,定义被划分成N个多面体网格的三维几何模型中叶片的载荷、与叶片所连接的转子的边界条件和叶片与转子间的接触面,利用定义后的信息及预知的叶片和转子的材料属性进行有限元分析,获得低压末级长叶片的强度和振动频率,再根据振动频率和调节后的施加叶片上的蒸汽压力进行动应力运算,获得低压末级长叶片的动应力。主要用于对核电机组低压末级长叶片进行强度、振动频率和动应力的分析。
Description
技术领域
本发明属于核电汽轮机组长叶片分析领域,具体涉及一种用于核电机组低压末级长叶片的有限元分析方法。
背景技术
叶片是汽轮机的心脏,是汽轮机最重要的部分,而低压末级长叶片是整个机组中最为重要的一级叶片。随着核电机组的不断增加,对于核电机组低压末级长叶片的有限元分析也需要提高准确性,一致性。
现有叶片有限元的分析方法无法应用于核电机组的低压末级长叶片,难以满足当今核电汽轮机的设计需求,且现有适用于常规汽轮机组的有限元的分析方法处理流程复杂,精度低;因此,如何提供一种分析精度高,且适用于核电机组的低压末级长叶片的有限元分析方法亟需解决。
发明内容
本发明目的是为了解决现有技术中常规汽轮机组的有限元的分析方法处理流程复杂、精度低及无法适用于核电机组的问题,提供了一种用于核电机组低压末级长叶片的有限元分析方法。该分析方法主要用于对核电机组低压末级长叶片的强度、振动频率和动应力进行分析。
一种用于核电机组低压末级长叶片的有限元分析方法,该有限元分析方法包括如下步骤:
步骤一、将低压末级长叶片及其所对应的转子部分作为整体,构建三维几何模型;
步骤二、将三维几何模型进行切割,划分成N个多面体网格;其中,N为大于20万且小于50万的整数;
步骤三、定义被划分成N个多面体网格的三维几何模型中低压末级长叶片的载荷、与低压末级长叶片所连接的转子的边界条件和低压末级长叶片与转子间的接触面;根据低压末级长叶片的载荷、与低压末级长叶片所连接的转子的边界条件、低压末级长叶片与转子间接触面、预知的低压末级长叶片的材料属性和预知的与低压末级长叶片所连接的转子的材料属性进行有限元分析,获得低压末级长叶片的强度和振动频率;
步骤四、调节施加在低压末级长叶片上的蒸汽压力,根据低压末级长叶片振动频率和调节后的施加在低压末级长叶片上的蒸汽压力进行动应力运算,获得低压末级长叶片的动应力,从而完成对低压末级长叶片的有限元分析。
优选的是,步骤二中,多面体网格为六面体和/或四面体网格。
优选的是,步骤三中,所定义的低压末级长叶片的载荷包括:施加在低压末级长叶片所对应的多面体网格上所有节点的转速、施加在低压末级长叶片上的工作温度及施加在低压末级长叶片上的蒸汽压力;
所述节点为每个多面体网格上面与面之间的交点。
优选的是,步骤三中,根据低压末级长叶片的载荷、与低压末级长叶片所连接的转子的边界条件、低压末级长叶片与转子间接触面、预知的低压末级长叶片的材料属性和预知的与低压末级长叶片所连接的转子的材料属性进行有限元分析,获得低压末级长叶片的强度和振动频率的过程采用CAE软件实现。
优选的是,步骤三中,材料属性包括不同温度下所对应的密度、弹性模量和线胀系数。
优选的是,步骤三中,与低压末级长叶片所连接的转子的边界条件为转子的径向约束、转子轴向约束和转子旋转方向的约束。
优选的是,当多面体网格为六面体网格时,低压末级长叶片上受蒸汽压力部分所对应的每个六面体网格的高度、长度和宽度的比值为2:1:1。
优选的是,步骤四中,调节后的施加在低压末级长叶片上的蒸汽压力为步骤三所定义的施加在低压末级长叶片上的蒸汽压力的0.04倍。
本发明带来的有益效果是,本发明所述的有限元分析方法分析过程简单,便于实现,主要用于对核电机组低压末级长叶片的强度、振动频率和动应力进行分析,分析的过程首先构建三维几何模型,对构建的三维几何模型划分网格,定义被划分成N个多面体网格的三维几何模型中低压末级长叶片的载荷、与低压末级长叶片所连接的转子的边界条件和低压末级长叶片与转子间的接触面,利用定义后的信息及预知的低压末级长叶片和转子的材料属性进行有限元分析,获得低压末级长叶片的强度和振动频率,再根据振动频率和调节后的施加在低压末级长叶片上的蒸汽压力进行动应力运算,获得低压末级长叶片的动应力,从而完成对低压末级长叶片的有限元分析。通过上述有限元分析过程,利用定义的参量和已知的参量可准确分析出各个工况下的叶片的强度、振动频率和动应力。
具体实施方式
下面是对本发明实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例仅仅是本发明一部分实施例,而不是全部的实施例。基于本发明中的实施例,本领域普通技术人员在没有作出创造性劳动的前提下所获得的所有其他实施例,都属于本发明保护的范围。
需要说明的是,在不冲突的情况下,本发明中的实施例及实施例中的特征可以相互组合。
本实施方式所述的一种用于核电机组低压末级长叶片的有限元分析方法,该有限元分析方法包括如下步骤:
步骤一、将低压末级长叶片及其所对应的转子部分作为整体,构建三维几何模型;
步骤二、将三维几何模型进行切割,划分成N个多面体网格;其中,N为大于20万且小于50万的整数;
步骤三、定义被划分成N个多面体网格的三维几何模型中低压末级长叶片的载荷、与低压末级长叶片所连接的转子的边界条件和低压末级长叶片与转子间的接触面;根据低压末级长叶片的载荷、与低压末级长叶片所连接的转子的边界条件、低压末级长叶片与转子间接触面、预知的低压末级长叶片的材料属性和预知的与低压末级长叶片所连接的转子的材料属性进行有限元分析,获得低压末级长叶片的强度和振动频率;
步骤四、调节施加在低压末级长叶片上的蒸汽压力,根据低压末级长叶片振动频率和调节后的施加在低压末级长叶片上的蒸汽压力进行动应力运算,获得低压末级长叶片的动应力,从而完成对低压末级长叶片的有限元分析。
本实施方式所述的本发明所述的有限元分析方法先构建三维几何模型,对构建的三维几何模型划分网格,定义被划分成N个多面体网格的三维几何模型中低压末级长叶片的载荷、与低压末级长叶片所连接的转子的边界条件和低压末级长叶片与转子间的接触面,利用定义后的信息及预知的低压末级长叶片和转子的材料属性进行有限元分析,获得低压末级长叶片的强度和振动频率,再根据振动频率和调节后的施加在低压末级长叶片上的蒸汽压力进行动应力运算,获得低压末级长叶片的动应力,从而完成对低压末级长叶片的有限元分析。该分析方法主要用于对核电机组低压末级长叶片的强度、振动频率和动应力进行分析,分析过程简单,便于实现,通过上述有限元分析过程,利用定义的参量和已知的参量可准确分析出各个工况下的叶片的强度、振动频率和动应力。
进一步的,步骤二中,多面体网格为六面体和/或四面体网格。
本优选实施方式中,多面体网格为六面体网格为最优。
更进一步的,步骤三中,所定义的低压末级长叶片的载荷包括:施加在低压末级长叶片所对应的多面体网格上所有节点的转速、施加在低压末级长叶片上的工作温度及施加在低压末级长叶片上的蒸汽压力;
所述节点为每个多面体网格上面与面之间的交点。
本优选实施方式中,通过对低压末级长叶片所对应的多面体网格上所有节点的转速及施加在低压末级长叶片上的工作温度和蒸汽压力进行相关条件的预设和定义,用以获得高精度的低压末级长叶片的强度、振动频率和动应力信息。
更进一步的,步骤三中,根据低压末级长叶片的载荷、与低压末级长叶片所连接的转子的边界条件、低压末级长叶片与转子间接触面、预知的低压末级长叶片的材料属性和预知的与低压末级长叶片所连接的转子的材料属性进行有限元分析,获得低压末级长叶片的强度和振动频率的过程采用CAE软件实现。
本优选实施方式中,CAE软件为现有的软件,操作简单便于实现。
更进一步的,步骤三中,材料属性包括不同温度下所对应的密度、弹性模量和线胀系数。
更进一步的,步骤三中,与低压末级长叶片所连接的转子的边界条件为转子的径向约束、转子轴向约束和转子旋转方向的约束。
更进一步的,当多面体网格为六面体网格时,低压末级长叶片上受蒸汽压力部分所对应的每个六面体网格的高度、长度和宽度的比值为2:1:1。
更进一步的,步骤四中,调节后的施加在低压末级长叶片上的蒸汽压力为步骤三所定义的施加在低压末级长叶片上的蒸汽压力的0.04倍。
虽然在本文中参照了特定的实施方式来描述本发明,但是应该理解的是,这些实施例仅仅是本发明的原理和应用的示例。因此应该理解的是,可以对示例性的实施例进行许多修改,并且可以设计出其他的布置,只要不偏离所附权利要求所限定的本发明的精神和范围。应该理解的是,可以通过不同于原始权利要求所描述的方式来结合不同的从属权利要求和本文中所述的特征。还可以理解的是,结合单独实施例所描述的特征可以使用在其他所述实施例中。
Claims (8)
1.一种用于核电机组低压末级长叶片的有限元分析方法,其特征在于,该有限元分析方法包括如下步骤:
步骤一、将低压末级长叶片及其所对应的转子部分作为整体,构建三维几何模型;
步骤二、将三维几何模型进行切割,划分成N个多面体网格;其中,N为大于20万且小于50万的整数;
步骤三、定义被划分成N个多面体网格的三维几何模型中低压末级长叶片的载荷、与低压末级长叶片所连接的转子的边界条件和低压末级长叶片与转子间的接触面;根据低压末级长叶片的载荷、与低压末级长叶片所连接的转子的边界条件、低压末级长叶片与转子间接触面、预知的低压末级长叶片的材料属性和预知的与低压末级长叶片所连接的转子的材料属性进行有限元分析,获得低压末级长叶片的强度和振动频率;
步骤四、调节施加在低压末级长叶片上的蒸汽压力,根据低压末级长叶片振动频率和调节后的施加在低压末级长叶片上的蒸汽压力进行动应力运算,获得低压末级长叶片的动应力,从而完成对低压末级长叶片的有限元分析。
2.根据权利要求1所述的一种用于核电机组低压末级长叶片的有限元分析方法,其特征在于,步骤二中,多面体网格为六面体和/或四面体网格。
3.根据权利要求1所述的一种用于核电机组低压末级长叶片的有限元分析方法,其特征在于,步骤三中,所定义的低压末级长叶片的载荷包括:施加在低压末级长叶片所对应的多面体网格上所有节点的转速、施加在低压末级长叶片上的工作温度及施加在低压末级长叶片上的蒸汽压力;
所述节点为每个多面体网格上面与面之间的交点。
4.根据权利要求1所述的一种用于核电机组低压末级长叶片的有限元分析方法,其特征在于,步骤三中,根据低压末级长叶片的载荷、与低压末级长叶片所连接的转子的边界条件、低压末级长叶片与转子间接触面、预知的低压末级长叶片的材料属性和预知的与低压末级长叶片所连接的转子的材料属性进行有限元分析,获得低压末级长叶片的强度和振动频率的过程采用CAE软件实现。
5.根据权利要求1所述的一种用于核电机组低压末级长叶片的有限元分析方法,其特征在于,步骤三中,材料属性包括不同温度下所对应的密度、弹性模量和线胀系数。
6.根据权利要求1所述的一种用于核电机组低压末级长叶片的有限元分析方法,其特征在于,步骤三中,与低压末级长叶片所连接的转子的边界条件为转子的径向约束、转子轴向约束和转子旋转方向的约束。
7.根据权利要求2所述的一种用于核电机组低压末级长叶片的有限元分析方法,其特征在于,当多面体网格为六面体网格时,低压末级长叶片上受蒸汽压力部分所对应的每个六面体网格的高度、长度和宽度的比值为2:1:1。
8.根据权利要求3所述的一种用于核电机组低压末级长叶片的有限元分析方法,其特征在于,步骤四中,调节后的施加在低压末级长叶片上的蒸汽压力为步骤三所定义的施加在低压末级长叶片上的蒸汽压力的0.04倍。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010693274.9A CN111814379A (zh) | 2020-07-17 | 2020-07-17 | 一种用于核电机组低压末级长叶片的有限元分析方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010693274.9A CN111814379A (zh) | 2020-07-17 | 2020-07-17 | 一种用于核电机组低压末级长叶片的有限元分析方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN111814379A true CN111814379A (zh) | 2020-10-23 |
Family
ID=72864955
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202010693274.9A Withdrawn CN111814379A (zh) | 2020-07-17 | 2020-07-17 | 一种用于核电机组低压末级长叶片的有限元分析方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN111814379A (zh) |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN112580164A (zh) * | 2020-12-14 | 2021-03-30 | 哈尔滨汽轮机厂有限责任公司 | 核电汽轮机低压长叶片设计方法 |
CN112699505A (zh) * | 2020-12-28 | 2021-04-23 | 哈尔滨汽轮机厂有限责任公司 | 一种用于核电机组低压缸长叶片的动应力有限元计算方法 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101046159A (zh) * | 2006-03-31 | 2007-10-03 | 通用电气公司 | 用于降低涡轮机叶片中的应力的方法和设备 |
CN108170943A (zh) * | 2017-12-26 | 2018-06-15 | 哈尔滨汽轮机厂有限责任公司 | 基于Python语言的汽轮机三维叶片设计中有限元前处理方法 |
CN110032814A (zh) * | 2019-04-18 | 2019-07-19 | 哈尔滨汽轮机厂有限责任公司 | 一种汽轮机t型叶根预扭叶片的有限元分析方法 |
CN110210159A (zh) * | 2019-06-11 | 2019-09-06 | 哈尔滨汽轮机厂有限责任公司 | 汽轮机末叶片锁紧装置处强度的有限元计算方法 |
-
2020
- 2020-07-17 CN CN202010693274.9A patent/CN111814379A/zh not_active Withdrawn
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101046159A (zh) * | 2006-03-31 | 2007-10-03 | 通用电气公司 | 用于降低涡轮机叶片中的应力的方法和设备 |
CN108170943A (zh) * | 2017-12-26 | 2018-06-15 | 哈尔滨汽轮机厂有限责任公司 | 基于Python语言的汽轮机三维叶片设计中有限元前处理方法 |
CN110032814A (zh) * | 2019-04-18 | 2019-07-19 | 哈尔滨汽轮机厂有限责任公司 | 一种汽轮机t型叶根预扭叶片的有限元分析方法 |
CN110210159A (zh) * | 2019-06-11 | 2019-09-06 | 哈尔滨汽轮机厂有限责任公司 | 汽轮机末叶片锁紧装置处强度的有限元计算方法 |
Non-Patent Citations (3)
Title |
---|
QIN FEI等: "Dynamic stress analysis of air-cooling steam turbine blade", 《JOURNAL OF BEIJING UNIVERSITY OF TECHNOLOGY》, pages 577 - 581 * |
袁海峰: "叶轮叶片振动模态分析与实验研究", 《中国优秀硕士学位论文全文数据库》, no. 12, pages 1 - 66 * |
马义良等: "火电机组"热电解耦"后低压末级叶片动应力分析", 《汽轮机技术》, vol. 60, no. 5, pages 343 - 345 * |
Cited By (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN112580164A (zh) * | 2020-12-14 | 2021-03-30 | 哈尔滨汽轮机厂有限责任公司 | 核电汽轮机低压长叶片设计方法 |
CN112580164B (zh) * | 2020-12-14 | 2022-07-22 | 哈尔滨汽轮机厂有限责任公司 | 核电汽轮机低压长叶片设计方法 |
CN112699505A (zh) * | 2020-12-28 | 2021-04-23 | 哈尔滨汽轮机厂有限责任公司 | 一种用于核电机组低压缸长叶片的动应力有限元计算方法 |
CN112699505B (zh) * | 2020-12-28 | 2022-11-25 | 哈尔滨汽轮机厂有限责任公司 | 一种用于核电机组低压缸长叶片的动应力有限元计算方法 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN111814379A (zh) | 一种用于核电机组低压末级长叶片的有限元分析方法 | |
Andrews | A finite element analysis of bending stresses induced in external and internal involute spur gears | |
CN106777783B (zh) | 一种航空发动机叶片裂纹预测方法 | |
Schoenenborn et al. | Contribution to free and forced vibration analysis of an intentionally mistuned blisk | |
Marchukov et al. | Improving of the working process of axial compressors of gas turbine engines by using an optimization method | |
CN110909433A (zh) | 燃气轮机压气机转子枞树型榫头-榫槽连接结构优化方法 | |
Carnegie et al. | The effects of shear deformation and rotary inertia on the lateral frequencies of cantilever beams in bending | |
CN112699505B (zh) | 一种用于核电机组低压缸长叶片的动应力有限元计算方法 | |
Delgado et al. | Reconstruction methodology of a Francis runner blade using numerical tools | |
Sladojevic et al. | Investigation of the influence of aerodynamic coupling on response levels of mistuned bladed discs with weak structural coupling | |
CN110688797B (zh) | 一种用于汽轮机叉型叶根叶片的有限元计算方法 | |
CN106777479A (zh) | 基于梁理论的涡轮叶片非线性蠕变分析方法 | |
Pinto et al. | Computational analysis of flow in turbines | |
CN115510713A (zh) | 汽轮机超低负荷叶片动应力三维计算系统及其计算方法 | |
CN106503375B (zh) | 一种基于cn群理论确定汽轮机转子临界转速的方法及系统 | |
CN114564864A (zh) | 一种风机轮毂强度分析方法、系统、设备及存储介质 | |
Boccini et al. | Towards structural topology optimization of rotor blisks | |
CN115859536B (zh) | 压气机转子叶片非同步振动锁频数值模拟方法 | |
Glanville | Investigation into core compressor tip leakage modelling techniques using a 3D viscous solver | |
Fruth et al. | Influence of the Blade Count Ratio on Aerodynamic Forcing: Part II—High Pressure Transonic Turbine | |
Yi et al. | Scheme optimization for a turbine blade under multiple working conditions based on the entropy weight vague set | |
Sladojevic et al. | Forced response variation of aerodynamically and structurally mistuned turbo-machinery rotors | |
Thapa et al. | Computational methods in research of hydraulic turbines operating in challenging environments | |
Meng et al. | Coupling simulation of the fast startup of a centrifugal pump with cavitation in a closed-loop pipeline system | |
Niu et al. | Research on Modal Analysis of Rotor Shaft System of Hydro-generator |
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 | ||
WW01 | Invention patent application withdrawn after publication |
Application publication date: 20201023 |
|
WW01 | Invention patent application withdrawn after publication |