CN106844839A - 用于优化汽轮机叶片型线的方法 - Google Patents

用于优化汽轮机叶片型线的方法 Download PDF

Info

Publication number
CN106844839A
CN106844839A CN201611154462.4A CN201611154462A CN106844839A CN 106844839 A CN106844839 A CN 106844839A CN 201611154462 A CN201611154462 A CN 201611154462A CN 106844839 A CN106844839 A CN 106844839A
Authority
CN
China
Prior art keywords
entry
point
circular arc
arc
spl
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.)
Granted
Application number
CN201611154462.4A
Other languages
English (en)
Other versions
CN106844839B (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.)
China Chang Jiang Energy Group Corp
Original Assignee
China Chang Jiang Energy Group Corp
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 China Chang Jiang Energy Group Corp filed Critical China Chang Jiang Energy Group Corp
Priority to CN201611154462.4A priority Critical patent/CN106844839B/zh
Publication of CN106844839A publication Critical patent/CN106844839A/zh
Application granted granted Critical
Publication of CN106844839B publication Critical patent/CN106844839B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

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与进/出汽边圆弧延伸的圆交于点N；当圆心C2与进/出汽边圆弧圆心C1的距离Dc小于R2时，过C1作方向与矢量C2C1相同的射线，使射线与进/出汽边圆弧延伸的圆交于点N’；点N或点N’即为经过优化的背弧样条曲线与进/出汽边圆弧H的连接点；
(2)迭代步骤：利用上述步骤(1)获得的点N或点N’，作为背弧样条曲线的初始端点，再和相邻的两个拟合点S1、S2，再次执行步骤(1)中的优化方法，直至点N或点N’与上一次迭代计算的背弧样条曲线的初始端点的距离小于设定的目标值，从而得到优化的背弧样条曲线与进/出汽边圆弧H的连接点的叶片型线。
3.如权利要求1或2任一项所述的一种用于优化汽轮机叶片型线的方法，其特征在于，当所述的初始端点S0和相邻的两个拟合点S1、S2在同一条直线上时，为初始端点S0的x坐标和y坐标各加一随机值,二随机值相同或不同，随机值的范围为[-1/10·R1,1/10·R1]，R1为所述进/出汽边圆弧H的半径。
4.如权利要求1或2任一项所述的一种用于优化汽轮机叶片型线的方法，其特征在于，迭代次数至多为10000次。
5.如权利要求4所述的一种用于优化汽轮机叶片型线的方法，其特征在于，当迭代次数大于1000次时，为初始端点S0的x坐标和y坐标各加一随机值，二随机值相同或不同，随机值的范围为[-1/10·R1,1/10·R1]，R1为所述进/出汽边圆弧H的半径。
6.如权利要求5所述的一种用于优化汽轮机叶片型线的方法，其特征在于，所述的叶片型线是动叶片型线或静叶片型线。
7.如权利要求1-6任一项所述的用于优化汽轮机叶片型线的方法，其特征在于，所述的叶片型线应用于1MW以上的汽轮机动叶片型线或静叶片型线上。
CN201611154462.4A 2016-12-14 2016-12-14 用于优化汽轮机叶片型线的方法 Active CN106844839B (zh)

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 true CN106844839A (zh) 2017-06-13
CN106844839B 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)

Cited By (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109063352A (zh) * 2018-08-10 2018-12-21 上海汽轮机厂有限公司 提高汽轮机叶片型面光顺度的方法
CN109359354A (zh) * 2018-09-28 2019-02-19 中船黄埔文冲船舶有限公司 一种挡水部件设计方法
CN110966969A (zh) * 2019-12-20 2020-04-07 杭州汽轮动力集团有限公司 一种燃气轮机涡轮静叶片关键结构检测装置和处理方法
CN111275792A (zh) * 2020-02-26 2020-06-12 江苏大学 一种混流式水轮机轴面流线绘制的新方法
CN113414102A (zh) * 2021-05-13 2021-09-21 天津水泥工业设计研究院有限公司 一种可调节成品粒度级配的双转子机构及其应用

Citations (8)

* Cited by examiner, † Cited by third party
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 如皋透平叶片制造有限公司 一种透平叶片型线的设计方法

Patent Citations (8)

* Cited by examiner, † Cited by third party
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)

* Cited by examiner, † Cited by third party
Title
毛鹏枭等: "汽轮机静叶片内背弧型线优化设计与研究", 《机械设计与制造》 *

Cited By (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109063352A (zh) * 2018-08-10 2018-12-21 上海汽轮机厂有限公司 提高汽轮机叶片型面光顺度的方法
CN109063352B (zh) * 2018-08-10 2023-01-06 上海汽轮机厂有限公司 提高汽轮机叶片型面光顺度的方法
CN109359354A (zh) * 2018-09-28 2019-02-19 中船黄埔文冲船舶有限公司 一种挡水部件设计方法
CN110966969A (zh) * 2019-12-20 2020-04-07 杭州汽轮动力集团有限公司 一种燃气轮机涡轮静叶片关键结构检测装置和处理方法
CN110966969B (zh) * 2019-12-20 2021-05-14 杭州汽轮动力集团有限公司 一种燃气轮机涡轮静叶片关键结构检测装置和处理方法
CN111275792A (zh) * 2020-02-26 2020-06-12 江苏大学 一种混流式水轮机轴面流线绘制的新方法
CN111275792B (zh) * 2020-02-26 2023-12-15 江苏大学 一种混流式水轮机轴面流线绘制的方法
CN113414102A (zh) * 2021-05-13 2021-09-21 天津水泥工业设计研究院有限公司 一种可调节成品粒度级配的双转子机构及其应用
CN113414102B (zh) * 2021-05-13 2022-05-13 天津水泥工业设计研究院有限公司 一种可调节成品粒度级配的双转子机构及其应用

Also Published As

Publication number Publication date
CN106844839B (zh) 2020-01-31

Similar Documents

Publication Publication Date Title
CN106844839A (zh) 用于优化汽轮机叶片型线的方法
CN106527351B (zh) 一种过渡曲线长度动态确定的小线段实时插补方法
CN104181860B (zh) 数控机床s型加减速控制方法
CN106640210A (zh) 一种径‑轴流式透平膨胀机叶轮叶片三维型线的设计方法
CN103631198B (zh) 基于G2连续Bézier曲线的刀具轨迹压缩方法
CN108563872B (zh) 网格参数化方法及基于该网格参数化方法的轴流涡轮气动优化设计方法
Reutter et al. Endwall contouring and fillet design for reducing losses and homogenizing the outflow of a compressor cascade
CN102023616B (zh) 三角Bézier曲面数控精加工刀轨快速生成方法
CN103279069A (zh) 用于数控机床的三角函数二阶连续可导加减速算法
CN105550388A (zh) 钣金工艺的工艺模型生成方法及装置
CN106125674A (zh) 一种高精度实时轮廓误差估计方法
CN103744346A (zh) 一种电子凸轮曲线生成方法
CN112084589A (zh) 基于轴面流速控制的轴贯流式水轮机转轮叶片优化方法
CN106020122A (zh) 基于牛顿迭代的数控轨迹控制方法
CN104317251A (zh) 基于Obrechkoff算法的三次NURBS曲线实时插补方法
CN103902786B (zh) 一种涡轮气冷动叶片伸根段外形参数化设计方法
CN116430796A (zh) 一种基于几何裁剪的加工路径局部光顺方法
CN106874526A (zh) 叶轮机叶片的生产坐标的生成方法和装置
CN110245407A (zh) 一种汽轮机叶片扭曲规律设计方法
CN107017820A (zh) 一种开关磁阻电机减振降噪系统
CN105574221A (zh) 一种改进cst翼型参数化方法
CN105240480B (zh) 基于水滴状翼型函数的液力变矩器叶片厚度确定方法
CN103942397A (zh) 基于幂函数的修形齿轮数字化建模方法
CN103870626A (zh) 一种径轴流式透平膨胀机叶轮子午面型线设计及校核方法
CN108304606B (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