CN110688797B - 一种用于汽轮机叉型叶根叶片的有限元计算方法 - Google Patents

一种用于汽轮机叉型叶根叶片的有限元计算方法 Download PDF

Info

Publication number
CN110688797B
CN110688797B CN201910916811.9A CN201910916811A CN110688797B CN 110688797 B CN110688797 B CN 110688797B CN 201910916811 A CN201910916811 A CN 201910916811A CN 110688797 B CN110688797 B CN 110688797B
Authority
CN
China
Prior art keywords
blade
fork
blade root
root
type
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
CN201910916811.9A
Other languages
English (en)
Other versions
CN110688797A (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.)
Harbin Turbine Co Ltd
Original Assignee
Harbin Turbine 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 Harbin Turbine Co Ltd filed Critical Harbin Turbine Co Ltd
Priority to CN201910916811.9A priority Critical patent/CN110688797B/zh
Publication of CN110688797A publication Critical patent/CN110688797A/zh
Application granted granted Critical
Publication of CN110688797B publication Critical patent/CN110688797B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • YGENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y02TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02TCLIMATE CHANGE MITIGATION TECHNOLOGIES RELATED TO TRANSPORTATION
    • Y02T90/00Enabling technologies or technologies with a potential or indirect contribution to GHG emissions mitigation

Landscapes

  • Turbine Rotor Nozzle Sealing (AREA)

Abstract

一种用于汽轮机叉型叶根叶片的有限元计算方法,它涉及一种有限元计算方法,具体涉及一种用于汽轮机叉型叶根叶片的有限元计算方法。本发明为了解决。本发明的具体步骤为在UG‑NX软件中创建三维实体计算模型:采用循环对称算法时叉型叶根叶片和轮缘要包含两个节距的实体;在ANSA软件中对实体模型进行网格划分:将实体模型导入到ANSA中,实体内部接触表面的部分细化网格,并保持实体形状与网格一致;将inp格式网格文件导入有限元分析软件Abaqus中,根据需求定义三维叶片材料数据、运行工况,定义接触、边界和载荷等设置;提交有限元计算方案。本发明属于汽轮机制造领域。

Description

一种用于汽轮机叉型叶根叶片的有限元计算方法
技术领域
本发明涉及一种有限元计算方法,具体涉及一种用于汽轮机叉型叶根叶片的有限元计算方法,属于汽轮机制造领域。
背景技术
叶片长期处在恶劣的工作环境,故障时有发生,强度分析和结构振动特性分析是考量汽轮机叶片安全性能的重要手段之一。
叶片强度与振动特性分析主要应用有限元分析法。有限元分析法能够较为精确的判断叶片的危险区域和最大应力的位置,着重关注叶片结构容易发生损坏的部位。而在进行整圈叶片结构振动特性分析时,如果按照实际整圈叶片结构创建模型进行有限元模拟分析,需要较大的计算机资源和较长计算时间,而且前处理工作也十分繁杂。
发明内容
本发明为解决按照实际整圈叶片结构创建模型进行有限元模拟分析,需要较大的计算机资源和较长计算时间,而且前处理工作也十分繁杂的问题,进而提出一种用于汽轮机叉型叶根叶片的有限元计算方法。
本发明为解决上述问题采取的技术方案是:本发明的具体步骤如下:
步骤一、在UG-NX软件中创建三维实体计算模型:
采用循环对称算法时叉型叶根叶片和轮缘要包含两个节距的实体,叉型叶根叶片的结构形式主要包括叉型叶根叶片的围带、叉型叶根叶片工作部分、叉型叶根叶片中间体、叉型叶根、销钉和叉型轮缘,对叉型叶根叶片的围带和叉型轮缘进行切割处理;
步骤二、在ANSA软件中对实体模型进行网格划分:
将实体模型导入到ANSA中,实体内部接触表面的部分细化网格,并保持实体形状与网格一致;
步骤三、将inp格式网格文件导入有限元分析软件Abaqus中,根据需求定义三维叶片材料数据、运行工况,定义接触、边界和载荷等设置;
步骤四、提交有限元计算方案。
本发明的有益效果是:本发明用于进行叉型叶根叶片的强度分析和结构振动特性分析,得到较为精确的叶片危险区域和最大应力的位置及整圈叶片频率,作为汽轮机叶片安全性的重要指导。
附图说明
图1是汽轮机叉型叶根叶片的有限元计算模型结构示意图;
图2是汽轮机叉型叶根叶片结构示意图;
图3是汽轮机转子结构示意图;
图4是销钉结构示意图;
图1中1-叉型叶根叶片的围带,2-叉型叶根叶片工作部分,3-叉型叶根叶片中间体,4-叉型叶根,5-销钉,6-叉型轮缘,7-转子;
图2中8-围带接触面,9-叉型叶根径向接触面,10-1-第一叶根销钉接触面,10-2-第二叶根销钉接触面,11-第三叶根销钉接触面,12-叶根齿轮槽接触面;
图3中13-第一轮槽销钉接触面,14-第二轮槽销钉接触面,15-轮槽叶根槽接触面,16-第三轮槽销钉接触面,17-1-第一转子端面,17-2-第二转子端面,18-转子中心。
具体实施方式
具体实施方式一:结合图1至图4说明本实施方式,本实施方式所述一种用于汽轮机叉型叶根叶片的有限元计算方法是通过如下步骤实现的:
步骤一、在UG-NX软件中创建三维实体计算模型:
采用循环对称算法时叉型叶根叶片和轮缘要包含两个节距的实体,叉型叶根叶片的结构形式主要包括叉型叶根叶片的围带1、叉型叶根叶片工作部分2、叉型叶根叶片中间体3、叉型叶根4、销钉5和叉型轮缘6,对叉型叶根叶片的围带1和叉型轮缘6进行切割处理;
步骤二、在ANSA软件中对实体模型进行网格划分:
将实体模型导入到ANSA中,实体内部接触表面的部分细化网格,并保持实体形状与网格一致;
步骤三、将inp格式网格文件导入有限元结构分析软件Abaqus中,根据需求定义三维叶片材料数据、运行工况,定义接触、边界和载荷等设置;
步骤四、提交有限元计算结果。
具体实施方式二:结合图1至图4说明本实施方式,本实施方式所述一种用于汽轮机叉型叶根叶片的有限元计算方法的步骤二中划分网格的基本原则是:
步骤A、将叉型叶根叶片的围带1自身画成一体;
步骤B、将叉型叶根叶片工作部分2画成一体;
步骤C、将叉型叶根4与叉型叶根叶片中间体3画成一体;
步骤D、将叉型轮缘6画成一体;
将划分好的叉型叶根叶片的围带1、叉型叶片、叉型叶根4和叉型轮缘6按照实际装配情况在ANSA中装配好后,以inp文件格式导出。
具体实施方式三:结合图1至图4说明本实施方式,本实施方式所述一种用于汽轮机叉型叶根叶片的有限元计算方法的步骤三中根据需求定义三维叶片材料数据、运行工况,定义接触、边界和载荷等设置的具体步骤如下:
步骤三(一)、在有限元网格计算模型上定义三维叶片材料数据,包括:
三维叶片材料在不同温度下的密度;
三维叶片材料在不同温度下的线胀系数;
三维叶片材料在不同温度下的弹性模量;
三维叶片材料在不同温度下的塑性数据;
三维叶片材料在不同温度下的导热系数;
三维叶片材料在不同温度下的比热;
将定义好的材料赋给相应的装配体;
步骤三(二)、在有限元网格计算模型上定义前处理需要的三维叶片面集合、节点集合、单元集合、载荷分析场,具体过程为:
定义二十一个三维叶片面集合、六个节点集合、一个单元集合、一个载荷分析场;
所述二十一个三维叶片面集合为:
叶片围带与围带之间的四个接触面;
叶片叶根与叶根之间的两个接触面;
叶根与轮槽接触面之间一个叶根面;
叶片叶根与销钉之间的三个叶根面;
叶根与轮槽接触面之间的一个轮槽面;
叶片轮槽与销钉之间的三个轮槽面;
销钉与叶片叶根、轮槽之间的三个销钉面;
循环对称面、主表面、从表面;
两个叶片压力表面;
所述六个节点集合为:
轮槽进汽侧端面节点、轮槽出汽侧端面节点、轮槽转子中心节点、销钉进汽侧端面中部节点、销钉出汽侧端面中部节点、用于加载热载荷的所有模型节点;
所述一个单元集合为:
用于加载离心力的所有模型单元;
所述一个载荷分析场为:
第一蒸汽压力场;
步骤三(三)、定义三维叶片运行工况包括:
叉型叶根叶片预应力;
叉型叶根叶片起机时强度;
叉型叶根叶片蒸汽弯应力强度;
叉型叶根叶片离心力强度;
叉型叶根叶片蠕变松弛强度;
叉型叶根叶片在跳机时的强度;
叉型叶根叶片在甩负荷时的强度;
叉型叶根叶片整圈动频率;
步骤三(四)、定义三维叶片相关设置:
定义叉型叶根叶片接触设置:叉型叶根叶片围带与围带之间、叉型叶片叶根与叶根之间、叶根与轮槽之间、叶根与销钉之间、轮槽与销钉之间等的位置和摩擦系数;
定义叉型叶根叶片循环对称设置:叉型叶根叶片只数、叶片旋转中心;
定义叉型叶根叶片位移约束设置:叉型叶根叶片转子中心、转子端面;
定义叉型叶根叶片载荷设置:叉型叶根叶片围带与围带之间的过盈载荷、叶根与叶根之间的过盈载荷、叶片蒸汽压力载荷、叶片离心力载荷;
定义叉型叶根叶片温度设置:三维叶片在起机时的、叶片在稳态时的温度、叶片在跳机或甩负荷时的温度、叶片在装配时的温度。
以上所述,仅是本发明的较佳实施例而已,并非对本发明作任何形式上的限制,虽然本发明已以较佳实施例揭露如上,然而并非用以限定本发明,任何熟悉本专业的技术人员,在不脱离本发明技术方案范围内,当可利用上述揭示的技术内容做出些许更动或修饰为等同变化的等效实施例,但凡是未脱离本发明技术方案内容,依据本发明的技术实质,在本发明的精神和原则之内,对以上实施例所作的任何简单的修改、等同替换与改进等,均仍属于本发明技术方案的保护范围之内。

Claims (1)

1.一种用于汽轮机叉型叶根叶片的有限元计算方法,其特征在于:所述一种用于汽轮机叉型叶根叶片的有限元计算方法是通过如下步骤实现的:
步骤一、在UG-NX软件中创建三维实体计算模型:
采用循环对称算法时叉型叶根叶片和轮缘要包含两个节距的实体,叉型叶根叶片的结构形式包括叉型叶根叶片的围带(1)、叉型叶根叶片工作部分(2)、叉型叶根叶片中间体(3)、叉型叶根(4)、销钉(5)和叉型轮缘(6),对叉型叶根叶片的围带(1)和叉型轮缘(6)进行切割处理;
步骤二、在ANSA软件中对实体模型进行网格划分:
将实体模型导入到ANSA中,实体内部接触表面的部分细化网格,并保持实体形状与网格一致;
步骤三、将inp格式网格文件导入有限元分析软件Abaqus中,根据需求定义三维叶片材料数据、运行工况,定义接触、边界和载荷设置;
步骤四、提交有限元计算方案;
步骤二中划分网格的基本原则是:
步骤A、将叉型叶根叶片的围带(1)自身画成一体;
步骤B、将叉型叶根叶片工作部分(2)画成一体;
步骤C、将叉型叶根(4)与叉型叶根叶片中间体(3)画成一体;
步骤D、将叉型轮缘(6)画成一体;
将划分好的叉型叶根叶片的围带(1)、叉型叶片、叉型叶根(4)和叉型轮缘(6)按照实际装配情况在ANSA中装配好后,以inp文件格式导出;
步骤三中根据需求定义三维叶片材料数据、运行工况,定义接触、边界和载荷设置的具体步骤如下:
步骤三(一)、在有限元网格计算模型上定义三维叶片材料数据,包括:
三维叶片材料在不同温度下的密度;
三维叶片材料在不同温度下的线胀系数;
三维叶片材料在不同温度下的弹性模量;
三维叶片材料在不同温度下的塑性数据;
三维叶片材料在不同温度下的导热系数;
三维叶片材料在不同温度下的比热;
将定义好的材料赋给相应的装配体;
步骤三(二)、在有限元网格计算模型上定义前处理需要的三维叶片面集合、节点集合、单元集合、载荷分析场,具体过程为:
定义二十一个三维叶片面集合、六个节点集合、一个单元集合、一个载荷分析场;
所述二十一个三维叶片面集合为:
叶片围带与围带之间的四个接触面;
叶片叶根与叶根之间的两个接触面;
叶根与轮槽接触面之间一个叶根面;
叶片叶根与销钉之间的三个叶根面;
叶根与轮槽接触面之间的一个轮槽面;
叶片轮槽与销钉之间的三个轮槽面;
销钉与叶片叶根、轮槽之间的三个销钉面;
循环对称面、主表面、从表面;
两个叶片压力表面;
所述六个节点集合为:
轮槽进汽侧端面节点、轮槽出汽侧端面节点、轮槽转子中心节点、销钉进汽侧端面中部节点、销钉出汽侧端面中部节点、用于加载热载荷的所有模型节点;
所述一个单元集合为:
用于加载离心力的所有模型单元;
所述一个载荷分析场为:
第一蒸汽压力场;
步骤三(三)、定义三维叶片运行工况包括:
叉型叶根叶片预应力;
叉型叶根叶片起机时强度;
叉型叶根叶片蒸汽弯应力强度;
叉型叶根叶片离心力强度;
叉型叶根叶片蠕变松弛强度;
叉型叶根叶片在跳机时的强度;
叉型叶根叶片在甩负荷时的强度;
叉型叶根叶片整圈动频率;
步骤三(四)、定义三维叶片相关设置:
定义叉型叶根叶片接触设置:叉型叶根叶片围带与围带之间、叉型叶片叶根与叶根之间、叶根与轮槽之间、叶根与销钉之间、轮槽与销钉之间的位置和摩擦系数;
定义叉型叶根叶片循环对称设置:叉型叶根叶片只数、叶片旋转中心;
定义叉型叶根叶片位移约束设置:叉型叶根叶片转子中心、转子端面;
定义叉型叶根叶片载荷设置:叉型叶根叶片围带与围带之间的过盈载荷、叶根与叶根之间的过盈载荷、叶片蒸汽压力载荷、叶片离心力载荷;
定义叉型叶根叶片温度设置:三维叶片在起机时的、叶片在稳态时的温度、叶片在跳机或甩负荷时的温度、叶片在装配时的温度。
CN201910916811.9A 2019-09-26 2019-09-26 一种用于汽轮机叉型叶根叶片的有限元计算方法 Active CN110688797B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201910916811.9A CN110688797B (zh) 2019-09-26 2019-09-26 一种用于汽轮机叉型叶根叶片的有限元计算方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201910916811.9A CN110688797B (zh) 2019-09-26 2019-09-26 一种用于汽轮机叉型叶根叶片的有限元计算方法

Publications (2)

Publication Number Publication Date
CN110688797A CN110688797A (zh) 2020-01-14
CN110688797B true CN110688797B (zh) 2023-07-28

Family

ID=69110402

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201910916811.9A Active CN110688797B (zh) 2019-09-26 2019-09-26 一种用于汽轮机叉型叶根叶片的有限元计算方法

Country Status (1)

Country Link
CN (1) CN110688797B (zh)

Families Citing this family (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN111814370A (zh) * 2020-06-24 2020-10-23 哈尔滨汽轮机厂有限责任公司 一种调节级叶片的有限元计算方法

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN105930620A (zh) * 2016-05-25 2016-09-07 济南大学 一种叶轮结构的热耦合振动特性有限元分析方法
CN107330204A (zh) * 2017-07-05 2017-11-07 沈阳鼓风机集团股份有限公司 闭式叶轮焊接工艺的有限元模拟分析方法及装置
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 哈尔滨汽轮机厂有限责任公司 汽轮机末叶片锁紧装置处强度的有限元计算方法

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN105930620A (zh) * 2016-05-25 2016-09-07 济南大学 一种叶轮结构的热耦合振动特性有限元分析方法
CN107330204A (zh) * 2017-07-05 2017-11-07 沈阳鼓风机集团股份有限公司 闭式叶轮焊接工艺的有限元模拟分析方法及装置
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 (6)

* Cited by examiner, † Cited by third party
Title
Experimental and numerical study of hollow core slabs strengthened with mounted steel bars and prestressed steel wire ropes;Xinghua Li等;《Construction and Building Materials》;20180825;第456-469页 *
基于ANSYS的汽轮机叶片参数化有限元分析;马征;《中国优秀硕士学位论文全文数据库工程科技Ⅱ辑》;20130715(第07期);C039-88 *
基于ANSYS的风力机基础结构非线性接触分析;方琳等;《黄河水利职业技术学院学报》;20180131;第30卷(第1期);第34-37页 *
某燃机动叶及压气机组件应力计算与振动分析;张强升;《中国优秀硕士学位论文全文数据库工程科技Ⅱ辑》;20130115(第01期);C039-94 *
核电汽轮机1710 mm末级长叶片整圈振动特性分析;周英;《热力透平》;20140930;第43卷(第3期);第199-202页 *
高阻尼末级长叶片调频方法研究;马义良等;《汽轮机技术》;20130630;第55卷(第3期);第192-194页 *

Also Published As

Publication number Publication date
CN110688797A (zh) 2020-01-14

Similar Documents

Publication Publication Date Title
CN110929357A (zh) 一种高性能舰船燃机压气机气动设计方法
CN108170943B (zh) 基于Python语言的汽轮机三维叶片设计中有限元前处理方法
CN112287580B (zh) 一种基于全三维数值模拟的轴流压气机喘振边界计算方法
CN110909433B (zh) 燃气轮机压气机转子枞树型榫头-榫槽连接结构优化方法
CN110688797B (zh) 一种用于汽轮机叉型叶根叶片的有限元计算方法
CN114117849A (zh) 低压涡轮叶/盘转子的叶冠阻尼减振分析方法
CN109578085B (zh) 一种通过导叶倾斜减弱涡轮动叶非定常作用力的方法
CN115392094A (zh) 一种基于热力耦合的涡轮盘结构优化方法
CN112182802B (zh) 计入随机失谐的叶轮机械气动弹性优化设计实现方法
CN111814379A (zh) 一种用于核电机组低压末级长叶片的有限元分析方法
CN111832126A (zh) 低压缸末级叶片静应力分析方法
CN115270567B (zh) 一种带凸肩的风扇转子叶片振动特性分析方法
CN112699505B (zh) 一种用于核电机组低压缸长叶片的动应力有限元计算方法
CN106777479A (zh) 基于梁理论的涡轮叶片非线性蠕变分析方法
CN115510713A (zh) 汽轮机超低负荷叶片动应力三维计算系统及其计算方法
CN111814370A (zh) 一种调节级叶片的有限元计算方法
CN112177677B (zh) 域扩展的带内环空腔的涡轮盘结构及其设计方法
Sladojevic et al. Investigation of the influence of aerodynamic coupling on response levels of mistuned bladed discs with weak structural coupling
CN106503375A (zh) 一种基于cn群理论确定汽轮机转子临界转速的方法及系统
CN112580164A (zh) 核电汽轮机低压长叶片设计方法
Hah et al. Unsteady aerodynamic flow phenomena in a transonic compressor stage
Zhou et al. Transonic compressor IGV/rotor interaction analysis including fluid structure interaction
Narkhede et al. Design and analysis of turbine blade of engine-a review
Li et al. Multiobjective optimization approach to turbomachinery blades design
Vahdati et al. An overview of time-domain computational methods for aeroelastic instabilities of multi-stage compressors

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
CB03 Change of inventor or designer information

Inventor after: Ma Yiliang

Inventor after: Gao Tieyin

Inventor after: Li Xinghua

Inventor after: Zhang Di

Inventor after: An Yang

Inventor after: Liang Tianfu

Inventor after: Li Yang

Inventor after: Wu Dujiang

Inventor after: Zhao Hongyu

Inventor after: Ma Xiaole

Inventor before: Zhang Di

Inventor before: Gao Tieyin

Inventor before: Li Xinghua

Inventor before: Ma Yiliang

Inventor before: An Yang

Inventor before: Liang Tianfu

Inventor before: Li Yang

Inventor before: Wu Dujiang

Inventor before: Zhao Hongyu

Inventor before: Ma Xiaole

CB03 Change of inventor or designer information
GR01 Patent grant
GR01 Patent grant