CN107895077B - 考虑盘鼓多制造因素的燃气轮机拉杆转子装配参数优化方法 - Google Patents
考虑盘鼓多制造因素的燃气轮机拉杆转子装配参数优化方法 Download PDFInfo
- 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
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
本发明一种考虑盘鼓多制造因素的燃气轮机拉杆转子装配参数优化方法,用于对燃气轮机拉杆转子盘鼓制造误差带来的转子初始弯曲和不平衡量同时进行优化。该方法采用轴向轮廓线来近似转子外轮廓面,采用轴向轮廓线长度的离散度来表征转子的初始弯曲;同时,通过各级盘鼓的圆周跳动度和盘鼓质量得到拉杆转子部分的不平衡量分布。该方法将轴向轮廓线长度的标准差和x、y方向的不平衡量绝对值同时作为优化目标函数,以各级盘鼓的装配相位角作为优化参数,编写了非支配排序多目标优化程序进行优化,优化后轴向轮廓线长度的标准差和x、y方向的转子不平衡量绝对值下降幅度均大于99%,减小了转子装配后的初始弯曲量和初始总体不平衡量,优化效果显著。
Description
技术领域
本发明属于燃气轮机技术领域,具体涉及一种考虑盘鼓多制造因素的燃气轮机拉杆转子装配参数优化方法。
背景技术
拉杆转子(包括中心拉杆和周向均布拉杆)是重型燃气轮机中一种常见的拉杆—盘鼓组合式结构,由拉杆提供预紧力,并通过拉杆螺栓将多级盘鼓压紧而成。由于拉杆转子具有重量轻(启动快)、热膨胀性能和刚性好的优点,因此被广泛应用于燃气轮机和航空发动机。燃气轮机盘鼓制造过程中同时存在端面平行度、圆周跳动度,前者导致实际的盘鼓接触面与设计状态下的接触面不一致,存在平行度偏差,装配过程中会使拉杆转子产生初始弯曲;后者在机组运行时会产生不平衡量,在转动时产生离心激振力。由于两者同时受盘鼓的安装角度影响,在优化时需要同时考虑。通过调整盘鼓的安装相位角,同时对转子的初始弯曲和不平衡量进行优化,得到最佳的盘鼓装配角,提高燃气轮机拉杆转子的装配效率。
发明内容
本发明的目的在于提供一种考虑盘鼓多制造因素的燃气轮机拉杆转子装配参数优化方法,该方法通过盘鼓轴向轮廓线的离散度判断拉杆转子初始弯曲大小;同时由各级盘鼓圆周跳动度和质量,得到燃气轮机拉杆转子不平衡量分布。通过调整各级盘鼓安装角度,同时减小拉杆转子的总体不平衡量和初始弯曲度。
为达到上述目的,本发明采用如下技术方案来实现的:
考虑盘鼓多制造因素的燃气轮机拉杆转子装配参数优化方法,包括以下步骤:
1)以转子上的键相0°为相位起始点,通过测量燃气轮机拉杆转子各级盘鼓平行度e1,得到盘鼓宽度的最大值和最小值所在位置及其大小;
2)以转子上的键相0°为相位起始点,通过测量燃气轮机拉杆转子各级盘鼓圆周跳动度e2,得到各级盘鼓偏心距及其相位,结合盘鼓质量得到各级盘鼓不平衡量分布;
3)以各级盘鼓上最小宽度所在位置为起点,逆时针沿圆周均匀划分为n份,得到n个圆周点,结合步骤1)中所测盘鼓宽度的最大值和最小值,线性得到各圆周点处的盘鼓宽度;
4)连接各级盘鼓上对应序号的圆周点,共形成n条盘鼓轴向轮廓线,结合在步骤3)中各圆周点处盘鼓宽度,得到各条轴向轮廓线长度;
5)将n条轴向轮廓线长度的标准差σ以及各级盘鼓不平衡量p的矢量之和在x、y方向分量的绝对值作为优化目标函数,以燃气轮机各盘鼓装配相位角为优化参数,建立非支配排序多目标优化目标函数进行优化,得到各级盘鼓最佳安装角度。
本发明进一步的改进在于,步骤3)中第i个盘鼓第k个圆周点处盘鼓宽度为:
其中,[]表示取整,n为盘鼓圆周总圆周点数,l为盘鼓宽度。
本发明进一步的改进在于,步骤4)中,第k条轴向轮廓线长度L表示为:
设第i个盘鼓逆时针旋转了j个单位划分角度,则有:
将(7)代入(6)中将对应盘鼓宽度替换即可得旋转后的轴向轮廓线长度。
本发明进一步的改进在于,步骤5)中,非支配排序多目标优化目标函数为:
其中,σ表示燃气轮机拉杆转子各轴向轮廓线长度的标准差/mm,Lk表示第k条轴向轮廓线长度/mm,表示各轴向轮廓线长度算数平均值/mm,n为圆周划分圆周点数;mi、ri分别为第i级盘鼓的质量/g和偏心距/mm,px、py分别为x方向和y方向不平衡量绝对值/g·mm,为第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个圆周点的相位角为
由于盘鼓两端均为平面,可以得到盘鼓宽度最大值与最小值所在点的相位差为π(参见图6),且宽度随着角度线性变化。设有h个盘鼓,则第i个盘鼓第k圆周点处盘鼓宽度为
其中,[]表示取整。连接各个盘鼓上相同编号的点,沿圆周共形成n条轴向轮廓线,第k条轴向轮廓线长度表示为:
设第i个盘鼓逆时针旋转了j个单位划分角度,则有:
将(7)代入(6)中替换相应的盘鼓宽度即可得旋转后的轴向轮廓线长度;
3)确定燃气轮机拉杆转子盘鼓圆周跳动度及不平衡量分布
其中,px、py分别为燃气轮机拉杆转子盘鼓x方向和y方向总不平衡量的绝对值/g·mm。
4)确定优化目标函数及优化参数并进行优化
考虑到燃气轮机拉杆转子装配过程中,盘鼓相位的变化会同时影响拉杆转子的初始弯曲和不平衡量分布,采用轴向轮廓线长度的离散度来表征拉杆转子的初始弯曲,确立非支配排序多目标优化目标函数如下:
其中,σ表示燃气轮机拉杆转子各轴向轮廓线长度的标准差/mm,Lk表示第k条轴向轮廓线长度/mm,表示各轴向轮廓线长度算数平均值/mm,n为圆周划分圆周点数;mi、ri分别为第i级盘鼓的质量/g和偏心距/mm,px、py分别为x方向和y方向不平衡量绝对值/g·mm,为第i级盘鼓不平衡量的相位角,h为盘鼓数量。
表1某型燃气轮机拉杆转子优化前与优化后各级盘鼓装配相位角:
表2某型燃气轮机拉杆转子优化前与优化后目标函数值对比:
由表2中计算结果可知,优化后该燃气轮机拉杆转子轴向轮廓线标准差、x方向和y方向的不平衡量绝对值得到了明显的降低,下降幅度均大于99%,减小了转子装配后的初始弯曲量和初始总体不平衡量,优化效果显著。优化后直接得到各个盘鼓相对于起始位置的旋转角度,便于现场装配时进行调节。
5)验证本发明对实际装配过程效率的提升
在燃气轮机拉杆转子盘鼓实际装配过程中,由于盘鼓制造误差的存在,需要对盘鼓相位角进行反复的调整,且不一定能得到最满意的结果。以100组随机生成的燃气轮机拉杆转子盘鼓装配相位角为例,模拟现场装配时装配质量的不确定性,计算该100组相位角下的目标函数值并与本发明方法优化后的相位角下目标函数值进行对比,如图8所示,其中星号(*)代表随机生成的100组相位角计算的目标函数值,三角符号(△)代表优化后的相位角计算的目标函数值。可见,本发明在已知盘鼓端面平行度、圆周跳动度的情况下,可以快速对拉杆转子初始弯曲和不平衡量分布进行优化,得到各级盘鼓的最佳装配相位角,提高了装配质量和效率。
Claims (2)
1.考虑盘鼓多制造因素的燃气轮机拉杆转子装配参数优化方法,其特征在于,包括以下步骤:
1)以转子上的键相0°为相位起始点,通过测量燃气轮机拉杆转子各级盘鼓平行度e1,得到盘鼓宽度的最大值和最小值所在位置及其大小;
2)以转子上的键相0°为相位起始点,通过测量燃气轮机拉杆转子各级盘鼓圆周跳动度e2,得到各级盘鼓偏心距及其相位,结合盘鼓质量得到各级盘鼓不平衡量分布;
3)以各级盘鼓上最小宽度所在位置为起点,逆时针沿圆周均匀划分为n份,得到n个圆周点,结合步骤1)中所测盘鼓宽度的最大值和最小值,线性得到各圆周点处的盘鼓宽度;其中第i个盘鼓第k个圆周点处盘鼓宽度为:
其中,[]表示取整,n为盘鼓圆周总圆周点数,l为盘鼓宽度,li,min和li,max分别表示第i个盘鼓的宽度的最小值和最大值;
4)连接各级盘鼓上对应序号的圆周点,共形成n条盘鼓轴向轮廓线,结合在步骤3)中各圆周点处盘鼓宽度,得到各条轴向轮廓线长度;其中第k条轴向轮廓线长度L表示为:
设第i个盘鼓逆时针旋转了j个单位划分角度,则有:
将(7)代入(6)中将对应盘鼓宽度替换即可得旋转后的轴向轮廓线长度;
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)
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)
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)
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)
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 |
-
2017
- 2017-11-10 CN CN201711107238.4A patent/CN107895077B/zh active Active
Patent Citations (4)
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)
Title |
---|
有初弯的单盘不平衡转子等加速过临界时的瞬态响应分析;杨恩东 等;《有初弯的单盘不平衡转子等加速过临界时的瞬态响应分析》;20141208(第12期);30-33 * |
燃气轮机拉杆转子动力学特性研究;王少波;《中国优秀硕士学位论文全文数据库 工程科技Ⅱ辑》;20130715(第7期);C039-169 * |
重型燃气轮机转子结构及动力学特性研究综述;袁奇 等;《热力透平》;20131215;第42卷(第4期);294-301 * |
Cited By (1)
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 |