CN110245407A - 一种汽轮机叶片扭曲规律设计方法 - Google Patents
一种汽轮机叶片扭曲规律设计方法 Download PDFInfo
- Publication number
- CN110245407A CN110245407A CN201910486400.0A CN201910486400A CN110245407A CN 110245407 A CN110245407 A CN 110245407A CN 201910486400 A CN201910486400 A CN 201910486400A CN 110245407 A CN110245407 A CN 110245407A
- Authority
- CN
- China
- Prior art keywords
- blade profile
- coordinate
- profile
- design method
- calculate
- 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
- 238000000034 method Methods 0.000 title claims abstract description 22
- 238000004458 analytical method Methods 0.000 claims abstract description 10
- 230000003068 static effect Effects 0.000 claims abstract description 9
- 238000006243 chemical reaction Methods 0.000 claims description 4
- 230000007812 deficiency Effects 0.000 description 2
- 238000010586 diagram Methods 0.000 description 2
- 238000006073 displacement reaction Methods 0.000 description 2
- 230000000694 effects Effects 0.000 description 2
- 238000003754 machining Methods 0.000 description 2
- 230000015556 catabolic process Effects 0.000 description 1
- 238000007796 conventional method Methods 0.000 description 1
- 238000006731 degradation reaction Methods 0.000 description 1
- 238000012804 iterative process Methods 0.000 description 1
- 238000011084 recovery Methods 0.000 description 1
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 description 1
Classifications
-
- 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
-
- 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/06—Power analysis or power optimisation
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
本发明公开了一种汽轮机叶片扭曲规律设计方法,具体步骤为:1.输入最终气动叶型,并设为初始叶型;2.将叶型沿径向分成n个截面,同时将每个截面上的叶型分成前缘、尾缘、压力面、吸力面4段,并对三维叶型进行离散化处理;3.对叶型施加离心力及汽流力,进行分析,得到各节点的坐标;4.计算迭代坐标差值,计算加工叶型Ak坐标(xi,j,k,yi,j,k,zi,j,k),计算收敛,输出加工叶型Ak坐标;5.计算不收敛,在步骤4获得的加工叶型的基础上,继续对其进行大变形非线性静力分析,获得所有节点的坐标满足条件,否则继续进行迭代计算。本发明公开的一种汽轮机叶片扭曲规律设计方法,可用于所有汽轮机扭叶片的设计中,使设计更为精细化,能有效提升通流实际性能。
Description
技术领域
本发明涉及汽轮机叶片技术领域,特别是汽轮机长扭叶片。
背景技术
汽轮机是将水蒸汽的热能转化为机械能的动力设备。其用途广泛,是现代火力发电厂的主要设备,同时也可作为拖动给水泵、压缩机等工业领域的驱动设备等。叶片,是汽轮机的核心,其型线的设计水平直接影响整机效率。气动设计给出的是最终气动叶型,而结构设计图需要的是机械加工叶型。传统设计往往将最终气动设计叶型作为结构设计和机械加工的叶型。而实际上,在离心力和汽流力的作用下叶片将产生扭转恢复和伸长,它不但使气动性能恶化,而且还会影响叶片的载荷及改变叶片截面稳态应力的分布。近年来随着汽轮机叶片的排气面积不断增加,使得叶型的这种变化对气动影响不断增大。同时,随着对叶片气动性能要求不断提高,对考虑最终气动叶型和加工叶型变化的设计提出更为严格的要求,本发明正是为了解决以上问题。
发明内容
本发明方法,克服现有技术的不足,提供一种汽轮机叶片扭曲规律设计方法,能够实现汽轮机叶片从最终气动叶型向加工叶型转换,克服传统方法将气动设计叶型作为结构设计和机械加工叶型造成性能衰减的不足。
本发明方法,提供一种汽轮机叶片扭曲规律设计方法,包括以下步骤:
步骤1)输入最终气动模型,并将其设为A0;
步骤2)将最终气动叶型沿径向分成n个截面;同时将每个截面上的叶型分成前缘、尾缘、压力面、吸力面4段。对三维叶型进行离散化处理,其中前、尾缘曲率半径小的区域节点数量为a,背弧、内弧曲率半径大的区域节点数量为b。初始叶型A0离散点坐标为(xi,j,0,yi,j,0,zi,j,0);
步骤3)对叶型施加离心力及汽流力,并进行大变形非线性静力分析,得到变形后各节点的坐标(x’i,j,k,y’i,j,k,z’i,j,k)(其中i为1到n的变量,n为截面数,j为1到m的变量,m为截面节点数,k为迭代次数,初值k=1);
步骤4)计算迭代坐标差值△uk=x’i,j,k-xi,j,0,△vk=y’i,j,k-yi,j,0,△wk=z’i,j,k-zi,j,0;计算加工叶型Ak坐标(xi,j,k,yi,j,k,zi,j,k),xi,j,k=xi,j,k-1-△uk,yi,j,k=yi,j,k-1-△vk,zi,j,k=zi,j,k-1-△wk;
1<i≤n,1<j≤m,满足时,计算收敛,输出加工叶型Ak坐标(xi,j,k,yi,j,k,zi,j,k);
步骤5)计算不收敛,在步骤4获得的加工叶型Ak的基础上,继续对其进行大变形非线性静力分析,获得所有节点的坐标(x’i,j,k,y’i,j,k,z’i,j,k),其中k=k+1,满足时,计算收敛,输出加工叶型Ak坐标(xi,j,k,yi,j,k,zi,j,k),不满足 时,进行迭代计算。
作为优选:
所述的步骤2)中沿叶高叶型截面的数目n为9~11。
所述的步骤2)中叶型前、尾缘节点分布数量在10~15范围内,背弧、内弧节点分布数量在60~70范围。
所述的步骤4)中收敛条件ε的值取0.02,收敛条件ε的取值参考目前相应机床加工精度水平。
如上所述,本发明一种汽轮机叶片扭曲规律设计方法,具有以下收益效果:
本发明采用了大变形非线性分析方法,考虑了离心力、汽流力对叶型扭曲和径向伸长的影响,提出了一种汽轮机叶片扭曲规律设计方法,实现了最终气动叶型向加工叶型的转换。
附图说明
图1为本发明设计方法流程图;
图2本发明设计方法步骤2的示意图;
图3本发明设计方法步骤3的示意图。
图4为本发明叶型迭代物理过程图;
具体实施方式
下面通过实施例结合附图对本发明进一步的描述。
本实施例为某汽轮机扭叶片,叶高247.4mm,该机组运行转速为3600rpm。
参考图1,本发明介绍了汽轮机叶片从最终气动叶型向加工叶型的转换过程,一种汽轮机叶片扭曲规律设计方法,其具体步骤为:
步骤1)、输入最终气动模型,并将其设为A0
参考图2,步骤2)、将实施例叶型沿径向分成10个截面,同时将每个截面上的型线分成前缘、尾缘、压力面、吸力面4段。对三维叶型进行离散化处理,前缘节点数量为15个,尾缘节点数量为10个,背弧节点数量为70个,内弧节点数量为60个。初始叶型A0离散点坐标为(xi,j,0,yi,j,0,zi,j,0)。
参考图3,步骤3)、对实施例叶型施加转速及汽流力,进行大变形非线性静力分析,得到变形后所有节点的坐标(x’i,j,k,y’i,j,k,z’i,j,k)(其中i为1到n的变量,n为截面数,j为1到m的变量,m为截面节点数,k为迭代次数,初值k=1)。
步骤4)计算迭代坐标差值△uk=x’i,j,k-xi,j,0,△vk=y’i,j,k-yi,j,0,△wk=z’i,j,k-zi,j,0;计算加工叶型Ak坐标(xi,j,k,yi,j,k,zi,j,k),xi,j,k=xi,j,k-1-△uk,yi,j,k=yi,j,k-1-△vk,zi,j,k=zi,j,k-1-△wk;1<i≤n,1<j≤m,满足时,计算收敛,输出加工叶型Ak坐标(xi,j,k,yi,j,k,zi,j,k);更为形象的物理过程如图4所示,给出了第1次迭代过程,对A0进行大变形非线性静力分析得到A0’叶型,并获得每个节点的位移矢量将A0叶型每个节点减去位移矢量得到第一次加工叶型A1。对A1叶型进行大变形非线性静力分析,得到第一次气动叶型A1’。
步骤5)、计算不收敛,在步骤4获得的加工叶型Ak的基础上,继续对其进行大变形非线性静力分析,获得所有节点的坐标(x’i,j,k,y’i,j,k,z’i,j,k),其中k=k+1,满足时,计算收敛,输出加工叶型Ak坐标(xi,j,k,yi,j,k,zi,j,k),不满足 时,进行迭代计算。
本实施例最终迭代5次,达到收敛条件,输出加工叶型Ak坐标(xi,j,k,yi,j,k,zi,j,k)。
Claims (1)
1.一种汽轮机叶片扭曲规律设计方法,其特征在于通过特定的流程和算法,采用迭代的方法实现最终气动叶型向加工叶型的转换过程。其具体步骤为:
步骤1)、输入最终气动模型,并将其设为A0;
步骤2)、将最终气动叶型沿径向分成n个截面;同时将每个截面上的叶型分成前缘、尾缘、压力面、吸力面4段,对三维叶型进行离散化处理,其中前、尾缘曲率半径小的区域节点数量为a,背弧、内弧曲率半径大的区域节点数量为b。初始叶型A0离散点坐标为(xi,j,0,yi,j,0,zi,j,0);
步骤3)、对叶型施加离心力及汽流力,并进行大变形非线性静力分析,得到变形后各节点的坐标(x’i,j,k,y’i,j,k,z’i,j,k),所述i为1到n的变量,n为截面数,j为1到m的变量,m为截面节点数,k为迭代次数,初值k=1;
步骤4)、计算迭代坐标差值△uk=x’i,j,k-xi,j,0,△vk=y’i,j,k-yi,j,0,△wk=z’i,j,k-zi,j,0;计算加工叶型Ak坐标(xi,j,k,yi,j,k,zi,j,k),xi,j,k=xi,j,k-1-△uk,yi,j,k=yi,j,k-1-△vk,zi,j,k=zi,j,k-1-△wk;1<i≤n,1<j≤m,满足时,计算收敛,输出加工叶型Ak坐标(xi,j,k,yi,j,k,zi,j,k);
步骤5)、计算不收敛,在步骤4获得的加工叶型Ak的基础上,继续对其进行大变形非线性静力分析,获得所有节点的坐标(x’i,j,k,y’i,j,k,z’i,j,k),其中k=k+1,满足时,计算收敛,输出加工叶型Ak坐标(xi,j,k,yi,j,k,zi,j,k),不满足 时,进行迭代计算。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910486400.0A CN110245407A (zh) | 2019-06-05 | 2019-06-05 | 一种汽轮机叶片扭曲规律设计方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910486400.0A CN110245407A (zh) | 2019-06-05 | 2019-06-05 | 一种汽轮机叶片扭曲规律设计方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN110245407A true CN110245407A (zh) | 2019-09-17 |
Family
ID=67886130
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201910486400.0A Pending CN110245407A (zh) | 2019-06-05 | 2019-06-05 | 一种汽轮机叶片扭曲规律设计方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN110245407A (zh) |
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111963253A (zh) * | 2020-09-03 | 2020-11-20 | 杭州汽轮机股份有限公司 | 电站全容量给水泵汽轮机次末级动叶片 |
CN112580164A (zh) * | 2020-12-14 | 2021-03-30 | 哈尔滨汽轮机厂有限责任公司 | 核电汽轮机低压长叶片设计方法 |
CN112627901A (zh) * | 2020-12-18 | 2021-04-09 | 杭州汽轮动力集团有限公司 | 一种大负荷透平末级动叶片 |
Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109595040A (zh) * | 2018-12-21 | 2019-04-09 | 杭州汽轮动力集团有限公司 | 一种燃气轮机叶片扭曲规律设计方法 |
-
2019
- 2019-06-05 CN CN201910486400.0A patent/CN110245407A/zh active Pending
Patent Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109595040A (zh) * | 2018-12-21 | 2019-04-09 | 杭州汽轮动力集团有限公司 | 一种燃气轮机叶片扭曲规律设计方法 |
Cited By (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111963253A (zh) * | 2020-09-03 | 2020-11-20 | 杭州汽轮机股份有限公司 | 电站全容量给水泵汽轮机次末级动叶片 |
CN111963253B (zh) * | 2020-09-03 | 2022-12-16 | 杭州汽轮动力集团股份有限公司 | 电站全容量给水泵汽轮机次末级动叶片 |
CN112580164A (zh) * | 2020-12-14 | 2021-03-30 | 哈尔滨汽轮机厂有限责任公司 | 核电汽轮机低压长叶片设计方法 |
CN112580164B (zh) * | 2020-12-14 | 2022-07-22 | 哈尔滨汽轮机厂有限责任公司 | 核电汽轮机低压长叶片设计方法 |
CN112627901A (zh) * | 2020-12-18 | 2021-04-09 | 杭州汽轮动力集团有限公司 | 一种大负荷透平末级动叶片 |
CN112627901B (zh) * | 2020-12-18 | 2022-08-12 | 杭州汽轮动力集团有限公司 | 一种大负荷透平末级动叶片 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN110245407A (zh) | 一种汽轮机叶片扭曲规律设计方法 | |
CN110929357A (zh) | 一种高性能舰船燃机压气机气动设计方法 | |
CN112528575A (zh) | 一种船舶燃气轮机多级轴流压气机多维度协同设计方法 | |
CN112069630B (zh) | 压气机及叶片、叶片的二维叶型设计方法、计算机设备 | |
CN107529643B (zh) | 一种兼顾水力、强度性能的混流式转轮多工况优化方法 | |
CN108563872B (zh) | 网格参数化方法及基于该网格参数化方法的轴流涡轮气动优化设计方法 | |
CN102032213A (zh) | 一种端区叶片前缘仿生处理方法 | |
US7941300B1 (en) | Process for the design of an airfoil | |
CN114036668B (zh) | 一种基于离心叶轮中间面曲线的可变自由曲面叶片生成方法 | |
CN115186398B (zh) | 一种轴流压气机进口导叶造型关键角度参数确定方法 | |
CN109578085B (zh) | 一种通过导叶倾斜减弱涡轮动叶非定常作用力的方法 | |
CN109815624A (zh) | 一种考虑进气总压畸变影响的压气机稳定边界判断方法 | |
Yue et al. | Design and optimization of tandem arranged cascade in a transonic compressor | |
CN111444636B (zh) | 基于基元级退化特性图的轴流压气机喘振裕度计算方法 | |
CN109595040A (zh) | 一种燃气轮机叶片扭曲规律设计方法 | |
CN114186513A (zh) | 一种具有反s型前缘的轴流压气机叶片造型设计方法 | |
CN101749052B (zh) | 空冷给水泵汽轮机的末级动叶片 | |
CN103629152A (zh) | 压气机及其叶片系统的加厚及布置方法 | |
CN113565572A (zh) | 一种小展弦比的单级轴流涡轮设计方法 | |
CN106874526A (zh) | 叶轮机叶片的生产坐标的生成方法和装置 | |
CN101737091B (zh) | 空冷给水泵汽轮机末级动叶片 | |
CN101158991A (zh) | 大小叶片压气机的气动设计方法 | |
CN110287647B (zh) | 一种跨声速压气机平面叶栅激波控制的设计方法 | |
CN105240480A (zh) | 基于水滴状翼型函数的液力变矩器叶片厚度确定方法 | |
CN204851335U (zh) | 一体化叶轮 |
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 | ||
RJ01 | Rejection of invention patent application after publication |
Application publication date: 20190917 |
|
RJ01 | Rejection of invention patent application after publication |