CN107895077B - 考虑盘鼓多制造因素的燃气轮机拉杆转子装配参数优化方法 - Google Patents

考虑盘鼓多制造因素的燃气轮机拉杆转子装配参数优化方法 Download PDF

Info

Publication number
CN107895077B
CN107895077B CN201711107238.4A CN201711107238A CN107895077B CN 107895077 B CN107895077 B CN 107895077B CN 201711107238 A CN201711107238 A CN 201711107238A CN 107895077 B CN107895077 B CN 107895077B
Authority
CN
China
Prior art keywords
drum
disk
rotor
gas turbine
pull rod
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
CN201711107238.4A
Other languages
English (en)
Other versions
CN107895077A (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.)
Xian Jiaotong University
Original Assignee
Xian Jiaotong University
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 Xian Jiaotong University filed Critical Xian Jiaotong University
Priority to CN201711107238.4A priority Critical patent/CN107895077B/zh
Publication of CN107895077A publication Critical patent/CN107895077A/zh
Application granted granted Critical
Publication of CN107895077B publication Critical patent/CN107895077B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00Computer-aided design [CAD]
    • G06F30/10Geometric CAD
    • G06F30/17Mechanical 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

本发明一种考虑盘鼓多制造因素的燃气轮机拉杆转子装配参数优化方法,用于对燃气轮机拉杆转子盘鼓制造误差带来的转子初始弯曲和不平衡量同时进行优化。该方法采用轴向轮廓线来近似转子外轮廓面,采用轴向轮廓线长度的离散度来表征转子的初始弯曲;同时,通过各级盘鼓的圆周跳动度和盘鼓质量得到拉杆转子部分的不平衡量分布。该方法将轴向轮廓线长度的标准差和x、y方向的不平衡量绝对值同时作为优化目标函数,以各级盘鼓的装配相位角作为优化参数,编写了非支配排序多目标优化程序进行优化,优化后轴向轮廓线长度的标准差和x、y方向的转子不平衡量绝对值下降幅度均大于99%,减小了转子装配后的初始弯曲量和初始总体不平衡量,优化效果显著。

Description

考虑盘鼓多制造因素的燃气轮机拉杆转子装配参数优化方法
技术领域
本发明属于燃气轮机技术领域,具体涉及一种考虑盘鼓多制造因素的燃气轮机拉杆转子装配参数优化方法。
背景技术
拉杆转子(包括中心拉杆和周向均布拉杆)是重型燃气轮机中一种常见的拉杆—盘鼓组合式结构,由拉杆提供预紧力,并通过拉杆螺栓将多级盘鼓压紧而成。由于拉杆转子具有重量轻(启动快)、热膨胀性能和刚性好的优点,因此被广泛应用于燃气轮机和航空发动机。燃气轮机盘鼓制造过程中同时存在端面平行度、圆周跳动度,前者导致实际的盘鼓接触面与设计状态下的接触面不一致,存在平行度偏差,装配过程中会使拉杆转子产生初始弯曲;后者在机组运行时会产生不平衡量,在转动时产生离心激振力。由于两者同时受盘鼓的安装角度影响,在优化时需要同时考虑。通过调整盘鼓的安装相位角,同时对转子的初始弯曲和不平衡量进行优化,得到最佳的盘鼓装配角,提高燃气轮机拉杆转子的装配效率。
发明内容
本发明的目的在于提供一种考虑盘鼓多制造因素的燃气轮机拉杆转子装配参数优化方法,该方法通过盘鼓轴向轮廓线的离散度判断拉杆转子初始弯曲大小;同时由各级盘鼓圆周跳动度和质量,得到燃气轮机拉杆转子不平衡量分布。通过调整各级盘鼓安装角度,同时减小拉杆转子的总体不平衡量和初始弯曲度。
为达到上述目的,本发明采用如下技术方案来实现的:
考虑盘鼓多制造因素的燃气轮机拉杆转子装配参数优化方法,包括以下步骤:
1)以转子上的键相0°为相位起始点,通过测量燃气轮机拉杆转子各级盘鼓平行度e1,得到盘鼓宽度的最大值和最小值所在位置及其大小;
2)以转子上的键相0°为相位起始点,通过测量燃气轮机拉杆转子各级盘鼓圆周跳动度e2,得到各级盘鼓偏心距及其相位,结合盘鼓质量得到各级盘鼓不平衡量分布;
3)以各级盘鼓上最小宽度所在位置为起点,逆时针沿圆周均匀划分为n份,得到n个圆周点,结合步骤1)中所测盘鼓宽度的最大值和最小值,线性得到各圆周点处的盘鼓宽度;
4)连接各级盘鼓上对应序号的圆周点,共形成n条盘鼓轴向轮廓线,结合在步骤3)中各圆周点处盘鼓宽度,得到各条轴向轮廓线长度;
5)将n条轴向轮廓线长度的标准差σ以及各级盘鼓不平衡量p的矢量之和在x、y方向分量的绝对值作为优化目标函数,以燃气轮机各盘鼓装配相位角
Figure BDA0001464663220000022
为优化参数,建立非支配排序多目标优化目标函数进行优化,得到各级盘鼓最佳安装角度。
本发明进一步的改进在于,步骤3)中第i个盘鼓第k个圆周点处盘鼓宽度为:
Figure BDA0001464663220000021
其中,[]表示取整,n为盘鼓圆周总圆周点数,l为盘鼓宽度。
本发明进一步的改进在于,步骤4)中,第k条轴向轮廓线长度L表示为:
Figure BDA0001464663220000031
设第i个盘鼓逆时针旋转了j个单位划分角度,则有:
Figure BDA0001464663220000032
将(7)代入(6)中将对应盘鼓宽度替换即可得旋转后的轴向轮廓线长度。
本发明进一步的改进在于,步骤5)中,非支配排序多目标优化目标函数为:
Figure BDA0001464663220000033
其中,σ表示燃气轮机拉杆转子各轴向轮廓线长度的标准差/mm,Lk表示第k条轴向轮廓线长度/mm,
Figure BDA0001464663220000034
表示各轴向轮廓线长度算数平均值/mm,n为圆周划分圆周点数;mi、ri分别为第i级盘鼓的质量/g和偏心距/mm,px、py分别为x方向和y方向不平衡量绝对值/g·mm,
Figure BDA0001464663220000035
为第i级盘鼓不平衡量的相位角,h为盘鼓数量。
与现有技术相比,本发明的有益效果在于:
本发明提供的考虑盘鼓多制造因素的燃气轮机拉杆转子装配参数优化方法,将盘鼓制造带来的端面平行度、圆周跳动度问题同时考虑,采用轴向轮廓线的离散度大小来表征转子的初始弯曲程度,通过优化得到各级盘鼓的相位角。通过该优化方法,轴向轮廓线长度的标准差和x方向、y方向的不平衡量绝对值下降幅度均大于99%,减小了转子装配后的初始弯曲量和初始总体不平衡量,优化效果显著;以100组随机生成的盘鼓装配相位角为例,模拟装配质量的不确定性,对比该100组相位角下的目标函数值与优化后的相位角下目标函数值,结果表明在已知各盘鼓端面平行度、圆周跳动度的情况时,本发明可以快速地同时对拉杆转子初始弯曲和不平衡量分布进行优化,可以有效提高装配质量和效率,减小燃气轮机拉杆转子的总体不平衡量和初始弯曲度。
附图说明
图1是某型燃气轮机周向均布拉杆转子结构示意图,图中给出了拉杆转子盘鼓和周向均布拉杆的结构及其装配关系。
图2至图4为几种典型的燃气轮机拉杆转子初始弯曲形式。
图5是某型燃气轮机拉杆转子盘鼓部分的端面平行度测点位置示意图,图中数字代表测点编号,共14个测点。
图6是燃气轮机拉杆转子盘鼓端面平行度测量、盘鼓坐标系的建立及盘鼓周向圆周点划分示意图,其中1,2,…,n-1,n表示沿圆周均匀分割为n等份。
图7是某型燃气轮机拉杆转子盘鼓部分的圆周跳动度测点位置示意图,图中数字代表测点编号,共14个测点。
图8是随机生成的100组燃气轮机拉杆转子盘鼓装配相位角与采用本发明优化后的相位角计算的目标函数值对比图。
具体实施方式
下面结合附图对本发明做出进一步的详细描述。
图1给出了某型燃气轮机拉杆转子盘鼓和拉杆的结构及其装配关系,其中盘鼓部分通过周向均布拉杆提供预紧力。
图2至图4给出了盘鼓端面平行度影响下,三种典型的燃气轮机拉杆转子盘鼓初始弯曲形式。由图2至图4可知,盘鼓装配相位角决定了燃气轮机拉杆转子初始弯曲的程度。
结合图1至图8,本发明考虑盘鼓多制造因素的燃气轮机拉杆转子装配参数优化方法,包括以下步骤:
1)确定燃气轮机拉杆转子盘鼓端面平行度及坐标系的建立
图5给出了某型燃气轮机拉杆转子端面平行度测点位置。以盘鼓一侧接触面为基准平面,通过对每一级盘鼓进行旋转,得到盘鼓宽度最大、最小值所在位置及其差值e1,测量盘鼓最窄处宽度lmin,则盘鼓最宽处宽度lmax=lmin+e1;以盘鼓中心指向键相0°所在位置为x轴,沿逆时针方向旋转π/2作为y轴建立直角坐标系。图6右边部分是某型燃气轮机拉杆转子平行度测量方法及坐标系建立示意图。
2)盘鼓沿周向的划分及确定盘鼓宽度沿周向的分布规律
如图6左边部分所示,将盘鼓沿周向平均分为n份,则第k个圆周点相对于第1个圆周点的相位角为
Figure BDA0001464663220000051
由于盘鼓两端均为平面,可以得到盘鼓宽度最大值与最小值所在点的相位差为π(参见图6),且宽度随着角度线性变化。设有h个盘鼓,则第i个盘鼓第k圆周点处盘鼓宽度为
Figure BDA0001464663220000052
其中,[]表示取整。连接各个盘鼓上相同编号的点,沿圆周共形成n条轴向轮廓线,第k条轴向轮廓线长度表示为:
Figure BDA0001464663220000061
设第i个盘鼓逆时针旋转了j个单位划分角度,则有:
Figure BDA0001464663220000062
将(7)代入(6)中替换相应的盘鼓宽度即可得旋转后的轴向轮廓线长度;
3)确定燃气轮机拉杆转子盘鼓圆周跳动度及不平衡量分布
图7给出了某型燃气轮机拉杆转子盘鼓圆周跳动度测点位置,其测量方法与端面平行度类似。通过测量可以得到第i级盘鼓圆周跳动值ei,2,则偏心距
Figure BDA0001464663220000063
设其相位角为
Figure BDA0001464663220000064
质量为mi,则有:
Figure BDA0001464663220000065
其中,px、py分别为燃气轮机拉杆转子盘鼓x方向和y方向总不平衡量的绝对值/g·mm。
4)确定优化目标函数及优化参数并进行优化
考虑到燃气轮机拉杆转子装配过程中,盘鼓相位的变化会同时影响拉杆转子的初始弯曲和不平衡量分布,采用轴向轮廓线长度的离散度来表征拉杆转子的初始弯曲,确立非支配排序多目标优化目标函数如下:
Figure BDA0001464663220000071
其中,σ表示燃气轮机拉杆转子各轴向轮廓线长度的标准差/mm,Lk表示第k条轴向轮廓线长度/mm,
Figure BDA0001464663220000072
表示各轴向轮廓线长度算数平均值/mm,n为圆周划分圆周点数;mi、ri分别为第i级盘鼓的质量/g和偏心距/mm,px、py分别为x方向和y方向不平衡量绝对值/g·mm,
Figure BDA0001464663220000073
为第i级盘鼓不平衡量的相位角,h为盘鼓数量。
将各级盘鼓装配相位角
Figure BDA0001464663220000074
作为优化参数,建立非支配排序多目标优化目标函数进行优化,得到各级盘鼓最佳安装角度。
表1某型燃气轮机拉杆转子优化前与优化后各级盘鼓装配相位角:
Figure BDA0001464663220000075
表2某型燃气轮机拉杆转子优化前与优化后目标函数值对比:
Figure BDA0001464663220000076
Figure BDA0001464663220000081
由表2中计算结果可知,优化后该燃气轮机拉杆转子轴向轮廓线标准差、x方向和y方向的不平衡量绝对值得到了明显的降低,下降幅度均大于99%,减小了转子装配后的初始弯曲量和初始总体不平衡量,优化效果显著。优化后直接得到各个盘鼓相对于起始位置的旋转角度,便于现场装配时进行调节。
5)验证本发明对实际装配过程效率的提升
在燃气轮机拉杆转子盘鼓实际装配过程中,由于盘鼓制造误差的存在,需要对盘鼓相位角进行反复的调整,且不一定能得到最满意的结果。以100组随机生成的燃气轮机拉杆转子盘鼓装配相位角为例,模拟现场装配时装配质量的不确定性,计算该100组相位角下的目标函数值并与本发明方法优化后的相位角下目标函数值进行对比,如图8所示,其中星号(*)代表随机生成的100组相位角计算的目标函数值,三角符号(△)代表优化后的相位角计算的目标函数值。可见,本发明在已知盘鼓端面平行度、圆周跳动度的情况下,可以快速对拉杆转子初始弯曲和不平衡量分布进行优化,得到各级盘鼓的最佳装配相位角,提高了装配质量和效率。

Claims (2)

1.考虑盘鼓多制造因素的燃气轮机拉杆转子装配参数优化方法,其特征在于,包括以下步骤:
1)以转子上的键相0°为相位起始点,通过测量燃气轮机拉杆转子各级盘鼓平行度e1,得到盘鼓宽度的最大值和最小值所在位置及其大小;
2)以转子上的键相0°为相位起始点,通过测量燃气轮机拉杆转子各级盘鼓圆周跳动度e2,得到各级盘鼓偏心距及其相位,结合盘鼓质量得到各级盘鼓不平衡量分布;
3)以各级盘鼓上最小宽度所在位置为起点,逆时针沿圆周均匀划分为n份,得到n个圆周点,结合步骤1)中所测盘鼓宽度的最大值和最小值,线性得到各圆周点处的盘鼓宽度;其中第i个盘鼓第k个圆周点处盘鼓宽度为:
Figure FDA0002273897070000011
其中,[]表示取整,n为盘鼓圆周总圆周点数,l为盘鼓宽度,li,min和li,max分别表示第i个盘鼓的宽度的最小值和最大值;
4)连接各级盘鼓上对应序号的圆周点,共形成n条盘鼓轴向轮廓线,结合在步骤3)中各圆周点处盘鼓宽度,得到各条轴向轮廓线长度;其中第k条轴向轮廓线长度L表示为:
Figure FDA0002273897070000012
设第i个盘鼓逆时针旋转了j个单位划分角度,则有:
Figure FDA0002273897070000013
将(7)代入(6)中将对应盘鼓宽度替换即可得旋转后的轴向轮廓线长度;
5)将n条轴向轮廓线长度的标准差σ以及各级盘鼓不平衡量p的矢量之和在x、y方向分量的绝对值作为优化目标函数,以燃气轮机各盘鼓装配相位角
Figure FDA0002273897070000021
为优化参数,建立非支配排序多目标优化目标函数进行优化,得到各级盘鼓最佳安装角度。
2.根据权利要求1所述的考虑盘鼓多制造因素的燃气轮机拉杆转子装配参数优化方法,其特征在于,步骤5)中,非支配排序多目标优化目标函数为:
Figure FDA0002273897070000022
其中,σ表示燃气轮机拉杆转子各轴向轮廓线长度的标准差/mm,Lk表示第k条轴向轮廓线长度/mm,
Figure FDA0002273897070000023
表示各轴向轮廓线长度算数平均值/mm,n为圆周划分圆周点数;mi、ri分别为第i级盘鼓的质量/g和偏心距/mm,px、py分别为x方向和y方向不平衡量绝对值/g·mm,
Figure FDA0002273897070000024
为第i级盘鼓不平衡量的相位角,h为盘鼓数量。
CN201711107238.4A 2017-11-10 2017-11-10 考虑盘鼓多制造因素的燃气轮机拉杆转子装配参数优化方法 Active CN107895077B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201711107238.4A CN107895077B (zh) 2017-11-10 2017-11-10 考虑盘鼓多制造因素的燃气轮机拉杆转子装配参数优化方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201711107238.4A CN107895077B (zh) 2017-11-10 2017-11-10 考虑盘鼓多制造因素的燃气轮机拉杆转子装配参数优化方法

Publications (2)

Publication Number Publication Date
CN107895077A CN107895077A (zh) 2018-04-10
CN107895077B true CN107895077B (zh) 2020-03-17

Family

ID=61804973

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201711107238.4A Active CN107895077B (zh) 2017-11-10 2017-11-10 考虑盘鼓多制造因素的燃气轮机拉杆转子装配参数优化方法

Country Status (1)

Country Link
CN (1) CN107895077B (zh)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
EP3940198A1 (en) * 2020-07-15 2022-01-19 Raytheon Technologies Corporation Multi-ring spacer for gas turbine engine rotor stack assembly

Families Citing this family (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110530306B (zh) * 2019-08-27 2020-10-20 大连理工大学 一种基于实测跳动数据的典型回转体零件表征方法
CN111982402B (zh) * 2020-04-09 2021-11-05 西北工业大学 考虑初始不平衡量的航空发动机转子零件选配优化方法
CN115826407B (zh) * 2022-11-29 2024-04-09 中国航发沈阳发动机研究所 一种减小鼓筒轴旋转惯性激励的控制方法

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103970944A (zh) * 2014-04-29 2014-08-06 西安交通大学 轮盘间平面接触的燃气轮机拉杆转子预紧力设计校核方法
CN105736059A (zh) * 2016-02-01 2016-07-06 西安交通大学 带端面齿的燃气轮机拉杆转子高速动平衡能力优化设计方法
CN106383930A (zh) * 2016-08-31 2017-02-08 中国石油大学(华东) 一种尾轴承‑转子系统多重流固耦合计算方法
CN106839968A (zh) * 2016-04-15 2017-06-13 南京航空航天大学 转子空间弯曲轴线测试系统及其测试方法

Family Cites Families (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US9638037B2 (en) * 2014-11-03 2017-05-02 Pratt & Whitney Canada Corp. Method of balancing a gas turbine engine rotor

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103970944A (zh) * 2014-04-29 2014-08-06 西安交通大学 轮盘间平面接触的燃气轮机拉杆转子预紧力设计校核方法
CN105736059A (zh) * 2016-02-01 2016-07-06 西安交通大学 带端面齿的燃气轮机拉杆转子高速动平衡能力优化设计方法
CN106839968A (zh) * 2016-04-15 2017-06-13 南京航空航天大学 转子空间弯曲轴线测试系统及其测试方法
CN106383930A (zh) * 2016-08-31 2017-02-08 中国石油大学(华东) 一种尾轴承‑转子系统多重流固耦合计算方法

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
有初弯的单盘不平衡转子等加速过临界时的瞬态响应分析;杨恩东 等;《有初弯的单盘不平衡转子等加速过临界时的瞬态响应分析》;20141208(第12期);30-33 *
燃气轮机拉杆转子动力学特性研究;王少波;《中国优秀硕士学位论文全文数据库 工程科技Ⅱ辑》;20130715(第7期);C039-169 *
重型燃气轮机转子结构及动力学特性研究综述;袁奇 等;《热力透平》;20131215;第42卷(第4期);294-301 *

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
EP3940198A1 (en) * 2020-07-15 2022-01-19 Raytheon Technologies Corporation Multi-ring spacer for gas turbine engine rotor stack assembly

Also Published As

Publication number Publication date
CN107895077A (zh) 2018-04-10

Similar Documents

Publication Publication Date Title
CN107895077B (zh) 考虑盘鼓多制造因素的燃气轮机拉杆转子装配参数优化方法
US8083476B2 (en) Method of designing a multistage turbine for a turbomachine
CN111523267B (zh) 一种基于参数化有限元模型的风机主轴结构优化方法
CA2905140C (en) Method of balancing a gas turbine engine rotor
CN109598081B (zh) 基于数据降维及多二维流面的径流式透平气动优化方法
JP2008157230A (ja) 回転部材をバランスさせるためのシステム
CN106446324B (zh) 大型工业汽轮机末级扭叶片设计方法
CN110259637B (zh) 风力发电机组的叶片气动不平衡矫正方法、装置及设备
CN111460677B (zh) 一种基于几何代数理论建立转子堆叠精度预测模型的方法
CN108799199B (zh) 高负荷宽裕度弯度可调串列叶栅的设计方法
CN111104713A (zh) 一种叶-盘结构耦合振动特性分析方法
US11698320B2 (en) Rotor balancing method and apparatus
US6152698A (en) Kit of articles and method for assembling articles along a holder distance
CN116663202B (zh) 一种多级轴流压气机性能仿真模型的校核方法
Salnikov et al. A centrifugal compressor impeller: a multidisciplinary optimization to improve its mass, strength, and gas-dynamic characteristics
CN111022135B (zh) 控制发动机动力涡轮叶尖的蓖齿间隙的装配方法
RU2449180C1 (ru) Способ балансировки ротора
CN116413020A (zh) 一种可调节转子叶片静力矩测量装置及使用方法
US9482237B1 (en) Method of designing a multi-stage turbomachine compressor
CA3097789A1 (en) Rotor balancing method and apparatus
CN115828446A (zh) 一种提高转子堆叠优化质量和计算速度的方法
CN115235404A (zh) 一种重型燃气轮机转静子间隙修正方法
CN115112081A (zh) 航空发动机转子装配相位优化方法及选配方法
CN114186449B (zh) 基于梁弯曲理论的叶盘扇区间刚度评估方法
CN113361029B (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