一种汽轮机T型叶根预扭叶片的有限元分析方法与流程

本发明涉及计算机辅助工程(cae)技术领域，具体为一种汽轮机t型叶根预扭叶片的有限元分析方法。

背景技术：

叶片是汽轮机的心脏，是极为重要的零件，它直接起着将蒸汽或燃气的热能转变为机械能的作用，叶根是将叶片固定在叶轮或转鼓上的连接部分，它应保证在任何运行条件下的连接牢靠，同时力求制造简单、装配方便。

t型叶根具有结构简单，加工和装配方便等特点，随着汽轮机中t型叶根预扭叶片使用的不断增加，对叶片的有限元分析提出了高效性、经济性、准确统一性等多方面的要求。但是预扭叶片的有限元分析不但要分析叶片强度与振动频率，同时还要考核叶片的预扭特性。

现有技术中仅有对预扭叶片的强度分析方法，但是并没有针对预扭叶片的振动频率和预扭特性的分析方法，这难以满足现今汽轮机领域的需求。

技术实现要素：

本发明的目的是：针对现有技术中仅有对预扭叶片强度的分析方法，并没有针对预扭叶片的振动频率和预扭特性的分析方法的从而导致叶片的分析结果与实际存在较大差异的问题，提供一种针对预扭叶片强度、振动频率和预扭特性的分析方法。

本发明采用如下技术方案实现：一种汽轮机t型叶根预扭叶片的有限元分析方法，包括以下步骤：

步骤一：几何模型处理

首先选取单级叶片、单级填隙条以及单级叶片对应的转子的部件几何模型，然后使用三维软件制作一个扇区内的几何模型；

步骤二：cae模型处理

首先对几何模型进行循环对称部分切割，然后对切割后的模型在网格划分软件中画六面体网格模型，最后在cae计算软件中根据需求定义叶片、转子、填隙条的材料数据，并定义接触、边界和载荷等设置，随后提交作业并进行计算；

步骤三：cae后处理

分析叶片各个工况下静应力，得到叶片所受静应力强度，使用频率safe图通过危险区域避开率来考核叶片振动特性，分析各个接触面的接触面积与接触应力大小得到预扭结构的预扭特性。

进一步的，所述cae计算软件为abaque软件。

进一步的，所述三维软件为nx。

进一步的，所述网格划分软件为ansa或hypermesh。

进一步的，在所述六面体网格模型中，叶片为整体的六面体网格，填隙条和转子为一体的六面体网格。

进一步的，所述叶片的部分非关心区域用四面体过度。

进一步的，所述的材料数据为叶片、转子、填隙条在不同温度下的密度、线胀系数、弹性模量、导热系数、比热等材料属性。

本发明采用上述技术方案，具有如下有益效果：本发明对汽轮机t型叶根预扭叶片的计算进行了标准化和规范化，通过分析叶片各个工况下静应力，考核叶片所受静应力强度。使用频率safe图分析叶片频率来考核叶片振动特性，通过分析各个接触面的接触面积与接触应力大小来分析预扭结构的预扭特性，本发明能有效分析叶片在汽轮机不同运行工况下的静强度，振动特性以及预扭特性，并且能提高30％工作效率。

附图说明

图1为一个扇区内单级预扭叶片几何模型示意图。

图2为本发明单级叶片对应转子选取方案图。

图3为叶片cae模型设置过程中定义接触面的con_bla_3的示意图。

图4为叶片cae模型设置过程中定义接触面的con_bla_3-1的示意图。

图5为cae模型设置过程中定义节点集axis示意图。

图6为cae模型设置过程中定义节点集bc_rot_us_a示意图。

图7为cae模型设置过程中定义节点集bc_rot_ds_a示意图。

图8为cae模型设置过程中定义节点集bc_rot_ds_1示意图。

图9为cae模型设置过程中定义节点集bc_rot_ds_1示意图。

图10为叶片各个工况下静应力示意图1。

图11为叶片各个工况下静应力示意图2。

图12为叶片各个工况下静应力示意图3。

图13为叶片各个工况下静应力示意图4。

图14为各个接触面的接触面积与接触应力大小示意图1。

图15为各个接触面的接触面积与接触应力大小示意图2。

图16为con_bla_1接触面定义图。

具体实施方式

具体实施方式一：下面结合图1和图2说明本实施方式，在本实施方式中，一种汽轮机t型叶根预扭叶片的有限元分析方法，包括以下步骤：

步骤一：几何模型处理

首先选取单级叶片、单级填隙条以及单级叶片对应的转子的部件几何模型，然后使用三维软件制作一个扇区内的几何模型；

步骤二：cae模型处理

首先对几何模型进行循环对称部分切割，然后对切割后的模型在网格划分软件中画六面体网格模型，最后在cae计算软件中根据需求定义叶片、转子、填隙条的材料数据，并定义接触、边界和载荷等设置，随后提交作业并进行计算；

步骤三：cae后处理

分析叶片各个工况下静应力，得到叶片所受静应力强度，使用频率safe图通过危险区域避开率来考核叶片振动特性，分析各个接触面的接触面积与接触应力大小得到预扭结构的预扭特性。

实施例：

本发明可以分为如下三个部分：几何模型处理、cae模型处理和cae后处理来具体实施。

1、几何模型处理

参见图1，几何模型主要包括选择单级叶片1，单级填隙条2，以及单级叶片对应转子3的选取。本发明使用循环对称模型来计算整圈叶片，根据选择使用三维软件(nx等)制作一个扇区内的几何模型。

选择转子部分的方案如图2所示，如果是计算中间级就选取图中的1-1和2-1两个边界进行截取创建几何模型。如果是首、末级叶片，其中一条截取边界为边界1-1，另一端截取则一定选择边界3-1/4-1/5-1/6-1。

2、cae模型处理

cae模型处理主要包括网格划分，以及在软件abaqus中的根据需求定义叶片、转子、填隙条的材料数据，定义接触、边界和载荷等设置。

本发明中用的循环对称模型首先要对实体进行切割，然后对切割后模型在网格划分软件(如ansa、hypermesh等)或者其他自动划分网格平台画六面体网格。对网格要求为：叶片为整体的六面体网格，部分非关心区域可以用四面体过度，填隙条和转子部分为一体的六面体网格。

本发明中制定了接触、载荷、场、边界等条件设置的标准方法，通过abaqus的standard模型进行计算。

本发明定位为弹性计算，所以在abaqus中分别定义叶片、填隙条以及转子材料在不同温度下的密度、线胀系数、弹性模量、导热系数、比热等材料属性。

在abaqus中定义接触、载荷、场、边界等，要根据如下标准定义有限元模型中的节点集、单元集以及表面。

单元element定义：

eblade-叶片所有单元；

erotor-转子所有单元；

eall-结构所有单元。

节点node集定义：

nbladenrotornall含义与单元定义所对应，

con_shi_2-填隙条与叶根槽接触部分定义点集

cyc_sym_slave-循环对称面从面点集定义

接触面定义：

con_rot_1,con_rot_2,con_rot_3,con_rot_4,con_rot_5,con_rot_1-1,con_rot_2-1,con_rot_3-1,con_rot_4-1-转子轮槽与叶片叶根接触面。

con_bla_1,con_bla_2,con_bla_3,con_bla_4,con_bla_5,con_bla_1-1,con_bla_2-1,con_bla_3-1,con_bla_4-1-叶片叶根与轮槽接触面。

con_bla_5是叶片下表面与填隙条接触的转子部分的表面，定义分别与叶片各表面相对应，每一个接触面如con_bla_1，因为循环对称被切开,要分别定义成con_bla_1和con_bla_1-1，参考图3和图14，图16中附图标记7为con_bla_1,附图标记8为con_bla_2，附图标记9为con_bla_3，附图标记10为con_bla_4。

con_shr_1,con_shr_2-围带接触对接触面；

con_blro_1，con_blro_2-叶根接触对接触面；

con_shi_1-填隙条接触面(填隙条与叶根接触的面)；

cyc_sym_master-循环对称面主面；

bla_pres_cfd-施加蒸汽力的表面。

参考图5至图9,为了优化预扭叶片的有限元计算方法，定义约束边界需要根据标准分别定义如下节点集：

axis-转子中心处节点集；

bc_rot_us_a-转子进气侧节点集；

bc_rot_ds_1和bc_rot_ds_a-转子出气侧节点集。

根据上述定义，在汽轮机t型叶根预扭叶片的有限元计算方法中规定了模型在abaqus中的接触、载荷、场、边界约束设置。具体如下表1所示，提交作业并进行计算。

表1：有限元模型接触、载荷、场、边界设置

3、cae后处理

通过分析叶片各个工况下静应力示意图如图10至图13所示，根据具体静应力考核标准考核叶片所受静应力强度。

使用频率safe图，通过危险区域避开率来考核叶片振动特性。

通过分析各个接触面的接触面积与接触应力大小如图14和图15所示，根据预扭标准来考核预扭叶片的预扭特性。

需要注意的是，具体实施方式仅仅是对本发明技术方案的解释和说明，不能以此限定权利保护范围。凡根据本发明权利要求书和说明书所做的仅仅是局部改变的，仍应落入本发明的保护范围内。