CN106844839B - 用于优化汽轮机叶片型线的方法 - Google Patents
用于优化汽轮机叶片型线的方法 Download PDFInfo
- Publication number
- CN106844839B CN106844839B CN201611154462.4A CN201611154462A CN106844839B CN 106844839 B CN106844839 B CN 106844839B CN 201611154462 A CN201611154462 A CN 201611154462A CN 106844839 B CN106844839 B CN 106844839B
- Authority
- CN
- China
- Prior art keywords
- point
- arc
- steam inlet
- spline curve
- outlet edge
- 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
Links
Images
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
Landscapes
- Physics & Mathematics (AREA)
- Geometry (AREA)
- Engineering & Computer Science (AREA)
- Theoretical Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Pure & Applied Mathematics (AREA)
- Mathematical Optimization (AREA)
- Mathematical Analysis (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- General Engineering & Computer Science (AREA)
- Computational Mathematics (AREA)
- Turbine Rotor Nozzle Sealing (AREA)
Abstract
本发明涉及一种用于优化汽轮机叶片型线的方法,其实是一种优化内弧/背弧样条曲线与进/出汽边圆弧H的连接点的方法,本发明通过迭代算法用于优化内弧/背弧样条曲线与进/出汽边圆弧H的连接点,将优化的叶片型线在汽轮机叶片中应用。本发明得到的优化方法能使得叶型设计中内弧/背弧样条曲线与进/出汽边圆弧H的连接点得到优化,相对于现有技术中未得到优化的叶型,能得到更为光顺平滑的叶型曲线,从而提高汽轮机叶片部件的流体动力学性能,且具有很高的计算效率,因而极富产业价值,兼具计算机辅助设计的优点。
Description
技术领域
本发明涉及一种用于优化汽轮机叶片型线的方法,适用于汽轮机的动静叶片。
背景技术
汽轮机叶片截面积的轮廓线叫叶片型线,它通常由非均匀有理B样条曲线构成的内弧、背弧以及两端的进汽边、出汽边圆弧组成。汽轮机中的蒸汽在各级叶片通道中完成从热能转化成机械能的整个过程,因此叶片型线对汽轮机通流部分的效率影响极大。
在汽轮机叶片部件的设计中,从流体动力学性能的要求看,其型线的光顺性非常重要。叶片型线多以内、背弧样条曲线拟合点和进、出汽边圆弧圆心坐标及半径的形式存储在型线图纸文件中。因内/背弧样条曲线与进/出汽边圆弧的连接点处于两条曲线的端点位置,起着连接和过渡两段曲线的作用。由于型线设计手段和数据存储精度的限制,可能引起叶片型线在内/背弧样条曲线与进/出汽边圆弧的连接点处不连续或者不相切,从而导致叶片型线的不连续和不平滑。
发明内容
本申请的目的是为了克服上述现有技术存在的问题,提供一种用于优化汽轮机叶片型线的方法,优化后的内/背弧样条曲线与进/出汽边圆弧的连接点更平滑,从而提高汽轮机叶片部件的流体动力学性能,减少能量损失,进而提高汽轮机机组效率。
本发明是一种计算机辅助设计叶型的方法,并利用了计算机的计算能力进行迭代计算。迭代算法是用计算机解决问题的一种基本方法。它利用计算机运算速度快、适合做重复性操作的特点,让计算机对设定的指令或一定步骤进行有限次的重复执行,在每次执行这组指令或这些步骤时,都从变量的原值推出它的一个新值,从而得到接近最优解的近似解。
汽轮机的叶型具体如图1所示,其进/出汽边曲线均由圆弧构成,内弧、背弧均为非均匀有理B样条曲线,进汽边与出汽边曲线分别与内弧和背弧曲线光滑过渡。
具体地,一种用于优化汽轮机叶片型线的方法,其特征在于,该方法为一种优化内弧样条曲线与进/出汽边圆弧h的连接点的方法,该方法包括:
(1)获得优化的内弧样条曲线与进/出汽边圆弧H的连接点:取内弧样条曲线的一个初始端点S0和相邻的两个拟合点S1、S2,计算同时经过此三点的圆,其圆心为C2,半径为R2;当圆心C2与进/出汽边的圆弧圆心C1的距离Dc大于或等于R2时,作线段连接C1与C2与进/出汽边圆弧延伸的圆交于点N;当圆心C2与进/出汽边圆弧圆心C1的距离Dc小于R2时,过C1作方向与矢量C2C1相同的射线,使射线与进/出汽边圆弧延伸的圆交于点N’;点N或点N’即为经过优化的内弧样条曲线与进/出汽边圆弧H的连接点;
(2)迭代步骤:利用上述步骤(1)获得的点N或点N’,作为内弧样条曲线的初始端点,再和相邻的两个拟合点S1、S2,再次执行步骤(1)中的优化方法,直至点N或点N’与上一次迭代计算的内弧样条曲线的初始端点的距离小于设定的目标值,从而得到优化的内弧样条曲线与进/出汽边圆弧H的连接点的叶片型线。
具体地,一种用于优化汽轮机叶片型线的方法,其特征在于,该方法进一步用于优化背弧样条曲线与进/出汽边圆弧h的连接点:
(1)获得优化的背弧样条曲线与进/出汽边圆弧H的连接点:取背弧样条曲线的一个初始端点S0和相邻的两个拟合点S1、S2,计算同时经过此三点的圆,其圆心为C2,半径为R2;当圆心C2与进/出汽边的圆弧圆心C1的距离Dc大于或等于R2时,作线段连接C1与C2与进/出汽边圆弧延伸的圆交于点N;当圆心C2与进/出汽边圆弧圆心C1的距离Dc小于R2时,过C1作方向与矢量C2C1相同的射线,使射线与进/出汽边圆弧延伸的圆交于点N’;点N或点N’即为经过优化的背弧样条曲线与进/出汽边圆弧H的连接点;
(2)迭代步骤:利用上述步骤(1)获得的点N或点N’,作为背弧样条曲线的初始端点,再和相邻的两个拟合点S1、S2,再次执行步骤(1)中的优化方法,直至点N或点N’与上一次迭代计算的背弧样条曲线的初始端点的距离小于设定的目标值,从而得到优化的背弧样条曲线与进/出汽边圆弧H的连接点的叶片型线。
当所述的初始端点S0和相邻的两个拟合点S1、S2在同一条直线上时,为初始端点S0的x坐标和y坐标各加一随机值,二随机值相同或不同,随机值的范围为[-1/10·R1,1/10·R1],R1为所述进/出汽边圆弧H的半径。
迭代次数至多为10000次。
当迭代次数大于1000次时,为初始端点S0的x坐标和y坐标各加一随机值,二随机值相同或不同,随机值的范围为[-1/10·R1,1/10·R1],R1为所述进/出汽边圆弧H的半径。
所述的叶片型线是动叶片型线或静叶片型线。
所述的叶片型线应用于1MW以上的汽轮机动叶片型线或静叶片型线上。
上述方法的计算速度一般在毫秒级别或更快,在极高效率的情况下瞬间得到接近最优解的近似解。
本发明所述的用于优化汽轮机叶片型线的方法,在汽轮机叶片设计中,应用于1MW以上的汽轮机动叶片型线或静叶片型线上,因而其通用性很强,适用于各种叶型设计,具有很大的产业推广应用价值。
在此基础上,本发明得到的优化方法能使得叶型设计中内弧/背弧样条曲线与进/出汽边圆弧H的连接点得到优化,相对于现有技术中未得到优化的叶型,能得到相对而言更为光顺平滑的叶型曲线,从而提高汽轮机叶片部件的流体动力学性能,且具有很高的计算效率,因而极富产业价值,兼具计算机辅助设计的优点。
附图说明
图1为本发明的叶片型线示意图。
图2a为本发明当Dc大于或等于R2时获得优化的内弧/背弧样条曲线与进/出汽边圆弧H的连接点的方法的示意图。
图2b为本发明当Dc小于R2时获得优化的内弧/背弧样条曲线与进/出汽边圆弧H的连接点的方法的示意图。
具体实施方式
下面结合附图对本发明进行详细的说明:
图1为本发明的叶片型线示意图,如图1所示,其进/出汽边均由圆弧构成,内弧、背弧均为非均匀有理B样条曲线,进汽边与出汽边曲线分别与内弧和背弧曲线光滑过渡。
按照本发明,提供了一种获得优化的内弧样条曲线与进/出汽边圆弧H的连接点的方法。该方法具体包括以下步骤:取内弧样条曲线的一个初始端点S0和相邻的两个拟合点S1、S2,计算同时经过此三点的圆,其圆心为C2,半径为R2;当圆心C2与进/出汽边的圆弧圆心C1的距离Dc大于或等于R2时,如图2a所示,作线段连接C1与C2与进/出汽边圆弧延伸的圆交于点N;当圆心C2与进/出汽边圆弧圆心C1的距离Dc小于R2时,如图2b所示,过C1作方向与矢量C2C1相同的射线,使射线与进/出汽边圆弧延伸的圆交于点N’。点N或点N’即为经过优化的内弧样条曲线与进/出汽边圆弧H的连接点。
利用上述优化内弧样条曲线与进/出汽边圆弧H的连接点的方法获得的点N或点N’,作为内弧样条曲线的初始端点,再次执行上述获得优化的内弧样条曲线与进/出汽边圆弧H的连接点的方法,直至点N或点N’与上一次迭代计算的内弧样条曲线的初始端点的距离小于目标值,得到内弧样条曲线的经过优化的端点,从而得到优化的内弧样条曲线与进/出汽边圆弧H的连接点的叶片型线。
同理按上述一种获得优化的内弧样条曲线与进/出汽边圆弧H的连接点的方法,来获得优化的背弧样条曲线与进/出汽边圆弧H的连接点及背弧样条曲线与进/出汽边圆弧H的连接点的叶片型线。
在另一种情况下,当所述的初始端点S0和相邻的两个拟合点S1、S2在同一条直线上时,为初始端点S0的X坐标和Y坐标各加一随机值,二随机值相同或不同,随机值的范围为[-1/10·R1,1/10·R1],R1为所述进/出汽边圆弧H的半径;虑及这一情况,是为了避免无法计算导致无法生成圆弧,在两坐标各加入随机值之后,可保证计算顺利进行。
其中,迭代次数至多为10000次;当迭代次数大于1000次时,为初始端点S0的X坐标和Y坐标各加一随机值,二随机值相同或不同,随机值的范围为[-1/10·R1,1/10·R1],R1为所述进/出汽边圆弧H的半径。
所述的叶片型线为动叶片型线或静叶片型线。
所述的叶片型线应用于1MW以上的汽轮机动叶片型线或静叶片型线上。
上述方法的计算速度一般在毫秒级别或更快,在极高效率的情况下瞬间得到接近最优解的近似解。
本发明所述的用于优化汽轮机叶片型线的方法,在汽轮机叶片设计中,应用于1MW以上的汽轮机动叶片型线或静叶片型线上,因而其通用性很强,适用于各种叶型设计,具有很大的产业推广应用价值。
Claims (7)
1.一种用于优化汽轮机叶片型线的方法,其特征在于,该方法为一种优化内弧样条曲线与进/出汽边圆弧H的连接点的方法,该方法包括:
(1)获得优化的内弧样条曲线与进/出汽边圆弧H的连接点:取内弧样条曲线的一个初始端点S0和相邻的两个拟合点S1、S2,计算同时经过此三点的圆,其圆心为C2,半径为R2;当圆心C2与进/出汽边的圆弧圆心C1的距离Dc大于或等于R2时,作线段连接C1与C2与进/出汽边圆弧延伸的圆交于点N;当圆心C2与进/出汽边圆弧圆心C1的距离Dc小于R2时,过C1作方向与矢量C2C1相同的射线,使射线与进/出汽边圆弧延伸的圆交于点N’;点N或点N’即为经过优化的内弧样条曲线与进/出汽边圆弧H的连接点;
(2)迭代步骤:利用上述步骤(1)获得的点N或点N’,作为内弧样条曲线的初始端点,再和相邻的两个拟合点S1、S2,再次执行步骤(1)中的优化方法,直至点N或点N’与上一次迭代计算的内弧样条曲线的初始端点的距离小于设定的目标值,从而得到优化的内弧样条曲线与进/出汽边圆弧H的连接点的叶片型线。
2.如权利要求1所述的一种用于优化汽轮机叶片型线的方法,其特征在于,该方法进一步用于优化背弧样条曲线与进/出汽边圆弧H的连接点:
(1)获得优化的背弧样条曲线与进/出汽边圆弧H的连接点:取背弧样条曲线的一个初始端点S0’和相邻的两个拟合点S1’、S2’,计算同时经过此三点的圆,其圆心为C2’,半径为R2’;当圆心C2’与进/出汽边的圆弧圆心C1的距离Dc’大于或等于R2’时,作线段连接C1与C2’与进/出汽边圆弧延伸的圆交于点N1;当圆心C2’与进/出汽边圆弧圆心C1的距离Dc’小于R2’时,过C1作方向与矢量C2’C1相同的射线,使射线与进/出汽边圆弧延伸的圆交于点N1’;点N1或点N1’即为经过优化的背弧样条曲线与进/出汽边圆弧H的连接点;
(2)迭代步骤:利用上述步骤(1)获得的点N1或点N1’,作为背弧样条曲线的初始端点,再和相邻的两个拟合点S1’、S2’,再次执行步骤(1)中的优化方法,直至点N1或点N1’与上一次迭代计算的背弧样条曲线的初始端点的距离小于设定的目标值,从而得到优化的背弧样条曲线与进/出汽边圆弧H的连接点的叶片型线。
3.如权利要求1或2任一项所述的一种用于优化汽轮机叶片型线的方法,其特征在于,当所述的初始端点S0/ S0’和相邻的两个拟合点S1/ S1’、S2/ S2’在同一条直线上时,为初始端点S0/ S0’的x坐标和y坐标各加一随机值,二随机值相同或不同,随机值的范围为[-1/10·R1, 1/10·R1],R1为所述进/出汽边圆弧H的半径。
4.如权利要求1或2任一项所述的一种用于优化汽轮机叶片型线的方法,其特征在于,迭代次数至多为10000次。
5.如权利要求4所述的一种用于优化汽轮机叶片型线的方法,其特征在于,当迭代次数大于1000次时,为初始端点S0/ S0’的x坐标和y坐标各加一随机值,二随机值相同或不同,随机值的范围为[-1/10·R1, 1/10·R1],R1为所述进/出汽边圆弧H的半径。
6.如权利要求5所述的一种用于优化汽轮机叶片型线的方法,其特征在于,所述的叶片型线是动叶片型线或静叶片型线。
7.如权利要求1或2任一项所述的用于优化汽轮机叶片型线的方法,其特征在于,所述的叶片型线应用于1MW以上的汽轮机动叶片型线或静叶片型线上。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201611154462.4A CN106844839B (zh) | 2016-12-14 | 2016-12-14 | 用于优化汽轮机叶片型线的方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201611154462.4A CN106844839B (zh) | 2016-12-14 | 2016-12-14 | 用于优化汽轮机叶片型线的方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN106844839A CN106844839A (zh) | 2017-06-13 |
CN106844839B true CN106844839B (zh) | 2020-01-31 |
Family
ID=59140543
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201611154462.4A Active CN106844839B (zh) | 2016-12-14 | 2016-12-14 | 用于优化汽轮机叶片型线的方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN106844839B (zh) |
Families Citing this family (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109063352B (zh) * | 2018-08-10 | 2023-01-06 | 上海汽轮机厂有限公司 | 提高汽轮机叶片型面光顺度的方法 |
CN109359354B (zh) * | 2018-09-28 | 2022-12-30 | 中船黄埔文冲船舶有限公司 | 一种挡水部件设计方法 |
CN110966969B (zh) * | 2019-12-20 | 2021-05-14 | 杭州汽轮动力集团有限公司 | 一种燃气轮机涡轮静叶片关键结构检测装置和处理方法 |
CN111275792B (zh) * | 2020-02-26 | 2023-12-15 | 江苏大学 | 一种混流式水轮机轴面流线绘制的方法 |
CN113414102B (zh) * | 2021-05-13 | 2022-05-13 | 天津水泥工业设计研究院有限公司 | 一种可调节成品粒度级配的双转子机构及其应用 |
Citations (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JPH11343998A (ja) * | 1998-06-02 | 1999-12-14 | Hitachi Ltd | 軸流圧縮機 |
CA2569566A1 (fr) * | 2005-11-28 | 2007-05-28 | Snecma | Circuit de refroidissement central pour aube mobile de turbomachine |
EP1828597A1 (en) * | 2004-11-18 | 2007-09-05 | Wind Power Limited | Vertical axis turbine apparatus |
CN101349168A (zh) * | 2008-09-13 | 2009-01-21 | 东方电气集团东方汽轮机有限公司 | 空冷汽轮机末级动叶片 |
CN101725373A (zh) * | 2009-11-20 | 2010-06-09 | 东方电气集团东方汽轮机有限公司 | 空冷汽轮机末级动叶片 |
CN101813004A (zh) * | 2010-04-15 | 2010-08-25 | 上海电气电站设备有限公司 | 叶轮机械叶片中弧线计算方法 |
CN102108881A (zh) * | 2011-03-22 | 2011-06-29 | 东方电气集团东方汽轮机有限公司 | 半转速核电汽轮机用的次末级叶片 |
CN104376144A (zh) * | 2014-08-06 | 2015-02-25 | 如皋透平叶片制造有限公司 | 一种透平叶片型线的设计方法 |
-
2016
- 2016-12-14 CN CN201611154462.4A patent/CN106844839B/zh active Active
Patent Citations (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JPH11343998A (ja) * | 1998-06-02 | 1999-12-14 | Hitachi Ltd | 軸流圧縮機 |
EP1828597A1 (en) * | 2004-11-18 | 2007-09-05 | Wind Power Limited | Vertical axis turbine apparatus |
CA2569566A1 (fr) * | 2005-11-28 | 2007-05-28 | Snecma | Circuit de refroidissement central pour aube mobile de turbomachine |
CN101349168A (zh) * | 2008-09-13 | 2009-01-21 | 东方电气集团东方汽轮机有限公司 | 空冷汽轮机末级动叶片 |
CN101725373A (zh) * | 2009-11-20 | 2010-06-09 | 东方电气集团东方汽轮机有限公司 | 空冷汽轮机末级动叶片 |
CN101813004A (zh) * | 2010-04-15 | 2010-08-25 | 上海电气电站设备有限公司 | 叶轮机械叶片中弧线计算方法 |
CN102108881A (zh) * | 2011-03-22 | 2011-06-29 | 东方电气集团东方汽轮机有限公司 | 半转速核电汽轮机用的次末级叶片 |
CN104376144A (zh) * | 2014-08-06 | 2015-02-25 | 如皋透平叶片制造有限公司 | 一种透平叶片型线的设计方法 |
Non-Patent Citations (1)
Title |
---|
汽轮机静叶片内背弧型线优化设计与研究;毛鹏枭等;《机械设计与制造》;20131031(第10期);第54-56页 * |
Also Published As
Publication number | Publication date |
---|---|
CN106844839A (zh) | 2017-06-13 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN106844839B (zh) | 用于优化汽轮机叶片型线的方法 | |
CN110727995B (zh) | 叶片形状的构造方法、叶片以及计算机设备 | |
US20160342730A1 (en) | Method for designing centrifugal pump and mixed flow pump having specific speed of 150-1200 | |
CN110929357A (zh) | 一种高性能舰船燃机压气机气动设计方法 | |
JP2009144716A (ja) | ターボ機械用の多段タービンを設計する方法 | |
CN108563872B (zh) | 网格参数化方法及基于该网格参数化方法的轴流涡轮气动优化设计方法 | |
EP3613995B1 (en) | Method for designing blade of axial-flow fluid machine and blade | |
CN106774151A (zh) | 一种可避免过切的刀具半径补偿方法 | |
CN110083968B (zh) | 基于修正气封泄露量影响数值模型的压气机特性预测方法 | |
Reutter et al. | Endwall contouring and fillet design for reducing losses and homogenizing the outflow of a compressor cascade | |
Komarov et al. | Application of optimisation techniques for new high-turning axial compressor profile topology design | |
CN115994423A (zh) | 基于Bezier曲线的跨音叶型设计方法及相关装置 | |
CN106202711B (zh) | 液力变矩器流道模型的参数化方法 | |
Van den Braembussche et al. | Multidisciplinary multipoint optimization of a transonic turbocharger compressor | |
CN103902786A (zh) | 一种涡轮气冷动叶片伸根段外形参数化设计方法 | |
CN106874526A (zh) | 叶轮机叶片的生产坐标的生成方法和装置 | |
Arnone et al. | Parametric optimization of a high-lift turbine vane | |
Safari et al. | A novel combination of adaptive tools for turbomachinery airfoil shape optimization using a real-coded genetic algorithm | |
Ji et al. | A Flow Information-Based Prediction Model Applied to the Nonaxisymmetric Hub Optimization of a Centrifugal Impeller | |
CN101813004B (zh) | 叶轮机械叶片中弧线计算方法 | |
CN105787217B (zh) | 一种飞机用波纹翼型的优化设计方法 | |
Li et al. | The optimization of a centrifugal impeller based on a new multi-objective evolutionary strategy | |
Komarov et al. | OPTIMIZATION APPROAC H AND SOME RESULTS FOR 2D COMPRESSOR AIRFOIL | |
Straccia et al. | Analysis of Cavity Leakage Effects on Coupled Non-Axisymmetric Endwall-Airfoil Optimization in a Low-Speed Compressor Tandem Stator | |
CN116127628A (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 | ||
GR01 | Patent grant | ||
GR01 | Patent grant |