CN114201879B - 一种燃气轮机的压气机涡轮转子的角向固定对接组合方法 - Google Patents

一种燃气轮机的压气机涡轮转子的角向固定对接组合方法 Download PDF

Info

Publication number
CN114201879B
CN114201879B CN202111519726.2A CN202111519726A CN114201879B CN 114201879 B CN114201879 B CN 114201879B CN 202111519726 A CN202111519726 A CN 202111519726A CN 114201879 B CN114201879 B CN 114201879B
Authority
CN
China
Prior art keywords
compressor
rotor
turbine rotor
node
imbalance
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
CN202111519726.2A
Other languages
English (en)
Other versions
CN114201879A (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 Hangfa Gas Turbine Co ltd
Original Assignee
China Hangfa Gas Turbine Co ltd
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 Hangfa Gas Turbine Co ltd filed Critical China Hangfa Gas Turbine Co ltd
Priority to CN202111519726.2A priority Critical patent/CN114201879B/zh
Publication of CN114201879A publication Critical patent/CN114201879A/zh
Application granted granted Critical
Publication of CN114201879B publication Critical patent/CN114201879B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00Computer-aided design [CAD]
    • G06F30/20Design optimisation, verification or simulation
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00Computer-aided design [CAD]
    • G06F30/20Design optimisation, verification or simulation
    • G06F30/23Design optimisation, verification or simulation using finite element methods [FEM] or finite difference methods [FDM]
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F2119/00Details relating to the type or aim of the analysis or the optimisation
    • G06F2119/14Force analysis or force optimisation, e.g. static or dynamic forces

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Computer Hardware Design (AREA)
  • Evolutionary Computation (AREA)
  • Geometry (AREA)
  • General Engineering & Computer Science (AREA)
  • General Physics & Mathematics (AREA)
  • Turbine Rotor Nozzle Sealing (AREA)
  • Testing Of Balance (AREA)
  • Supercharger (AREA)

Abstract

一种燃气轮机的压气机涡轮转子的角向固定对接组合方法,包括以下步骤:步骤1,压气机转子动平衡的残余不平衡量记作m1,相位记作α1,涡轮转子动平衡的残余不平衡量记作m2,相位记作α2;步骤2,建立压气机涡轮转子的转子动力学有限元模型;步骤3,利用转子动力学方程计算该状态下各个节点的最大不平衡响应R,得到全部状态下的最大不平衡响应,分别为(R1,θ1)、(R2,θ2)……(Rm,θm);选出(R1,θ1)、(R2,θ2)……(Rm,θm)中最大不平衡响应R的最小值Rmin,此时Rmin对应的安装角度差记为θmin,相应的对接状态(Rmin,θmin)作为最终装配方案,可以最大程度降低燃机振动。本发明方法能够有效降低燃机振动,同时不改变当前零件状态,便于施行。

Description

一种燃气轮机的压气机涡轮转子的角向固定对接组合方法
技术领域
本发明属于燃气轮机技术领域,具体涉及一种燃气轮机的压气机涡轮转子的角向固定对接组合方法。
背景技术
组成燃机的转动部件必须进行动平衡,通过调整配重使残余不平衡量满足设计要求,有利于降低振动。比如一般情况下,压气机转子和涡轮转子分别进行动平衡,对接后再进行组合动平衡。有些燃机的压气机转子和涡轮转子对接后结构较长,受设备工作范围的长度限制,不能实现组合动平衡。目前的做法将压气机转子与涡轮转子在角向相位随机状态下对接,不进行组合动平衡而直接装机,不利于降低振动,导致燃机不同台份之间的振动差异较大,甚至某些台份的振动超标。当振动超标时,必须将燃机分解到部件状态,随机改变压气机转子与涡轮转子的角向相位后复装。这种压气机涡轮转子的角向随机对接组合方法在燃机的装配过程中未实现当前零件状态下燃机的振动最小,在燃机振动的排故过程中存在不确定性;为此我们提出一种燃气轮机的压气机涡轮转子的角向固定对接组合方法。
发明内容
本发明的目的在于提供一种燃气轮机的压气机涡轮转子的角向固定对接组合方法,在不改变零部件状态的情况下使燃机振动最小,确定压气机转子和涡轮转子的角向相位。
为了实现上述目的,本发明采用如下技术方案:
一种燃气轮机的压气机涡轮转子的角向固定对接组合方法,包括以下步骤:
步骤1,压气机转子在动平衡机上完成单独动平衡的残余不平衡量记作m1,相位记作α1,涡轮转子在动平衡机上完成单独动平衡的残余不平衡量记作m2,相位记作α2;其中相位表示的是残余不平衡量相对确定好的参考标记的角度;
步骤2,利用Matlab建立压气机涡轮转子的转子动力学有限元模型,其中转子采用梁单元建模,Node 1-Node 21代表转子的21个节点,Node 1- Node 11组成压气机转子,Node 11- Node 21组成涡轮转子;支承采用轴承单元建模,单元1、3、5分别对应低压涡轮转子的三个支承;
步骤3,将涡轮转子残余不平衡量与压气机转子残余不平衡量相位之间的任意差值用θ表示,压气机涡轮转子系统的各节点的不平衡响应用q表示,压气机涡轮转子系统各个节点的最大不平衡响应用R表示;
将压气机转子与涡轮转子角向对接,此时以压气机转子残余不平衡量的相位α1为基准,将此状态下的涡轮转子残余不平衡量与压气机转子残余不平衡量相位之间的差值记作θ1,将此状态下的压气机涡轮转子系统的各节点的不平衡响应记作q1,利用转子动力学方程计算该状态下的各个节点的最大不平衡响应R1
压气机涡轮转子系统的转子动力学方程为:
(1)
其中M 、C 、K 、G 分别为转子系统的质量、阻尼、刚度以及陀螺力矩矩阵;Ω为给定的转速;q1为压气机涡轮转子系统的各节点的不平衡响应,描述每个节点的响应,每个节点有4个自由度,包含2个平动、2个转动自由度,整个转子系统具有4n个自由度,n为压气机涡轮转子模型的节点数;Q为作用在转子系统上的残余不平衡向量:
(2)
其中Qi用来描述压气机转子位于第i个节点的残余不平衡, Qj用来描述涡轮转子位于第j个节点的残余不平衡,此时涡轮转子残余不平衡的相位相对压气机转子残余不平衡量相位之间的差值为θ1,有:
(3)
(4)
则压气机涡轮转子系统的各节点的不平衡响应q1的计算公式为:
(5)
不考虑压气机涡轮转子系统的各节点的不平衡响应q中的转动响应,则各个节点的最大不平衡响应R1
R1=max(q1)(6)
将这组状态记作(R1,θ1);
步骤4,以压气机转子为基准,涡轮转子按照逆时针方向相对压气机转子转动1个花键套齿,此状态下的涡轮转子残余不平衡量与压气机转子残余不平衡量相位之间的差值记作θ2,将此状态下的压气机涡轮转子系统的各节点的不平衡响应记作q2,重复步骤3计算该状态下各个节点的最大不平衡响应R2,将这组状态记作(R2,θ2);
步骤5,重复步骤4,得到全部状态下的最大不平衡响应,分别为(R1,θ1)、(R2,θ2)……(Rm,θm);
步骤6,选出(R1,θ1)、(R2,θ2)……(Rm,θm)中最大不平衡响应R的最小值Rmin :
Rmin =min(R1,R2,...,Rm)(7)
此时Rmin对应的安装角度差记为θmin,相应的对接状态(Rmin,θmin)作为最终装配方案。
本发明的技术效果为:
与当前的角向随机对接组合方法相比,压气机涡轮转子采用角向固定对接组合方法,基于已有的压气机转子和涡轮转子的单独动平衡结果计算最优角向,能够有效降低燃机振动,同时不改变当前零件状态,便于施行。
附图说明
图1 本发明压气机涡轮转子系统结构示意图;
图2本发明压气机涡轮转子角向固定对接组合方法流程图;
图3为压气机涡轮转子的转子系统动力学有限元模型;
图4为涡轮转子残余不平衡的相位相对压气机转子残余不平衡相位的差值θ的示意图;
1-前支点,2-压气机转子,3-中支点,4-涡轮转子,5-后支点,6-花键套齿。
具体实施方式
下面结合附图和实施例对本发明作进一步的详细说明。
一种燃气轮机的压气机涡轮转子的角向固定对接组合方法,包括以下步骤:
步骤1,如图1和图2所示,压气机转子2在动平衡机上完成单独动平衡的残余不平衡量记作m1,m1=30g.cm,相位记作α1,α1=0;涡轮转子在动平衡机上完成单独动平衡的残余不平衡量记作m2,m2=30g.cm,相位记作α2,α2=0;其中相位表示的是残余不平衡量相对确定好的参考标记的角度;
步骤2,压气机涡轮转子系统的前端设置前支点1,中部设置中支点3,后端设置后支点5,并且两者通过花键套齿6连接,花键套齿的齿数目为16,利用Matlab建立压气机涡轮转子的转子动力学有限元模型,如图3所示,其中转子采用梁单元建模,Node 1-Node 21代表转子的21个节点,Node 1- Node 11组成压气机转子2,Node 11- Node 21组成涡轮转子4;支承采用轴承单元建模,单元1、3、5分别对应低压涡轮转子的三个支承;
步骤3,如图4所示,将涡轮转子残余不平衡量与压气机转子残余不平衡量相位之间的任意差值用θ表示,压气机涡轮转子系统的各节点的不平衡响应用q表示,压气机涡轮转子系统各个节点的最大不平衡响应用R表示;
将压气机转子2与涡轮转子4角向对接,此时以压气机转子残余不平衡量的相位α1为基准,将此状态下的涡轮转子残余不平衡量与压气机转子残余不平衡量相位之间的差值记作θ1,θ1=0°,将此状态下的压气机涡轮转子系统的各节点的不平衡响应记作q1,利用转子动力学方程计算该状态下的各个节点的最大不平衡响应R1
压气机涡轮转子系统的转子动力学方程为:
(1)
其中M为转子系统的质量,C为转子系统的阻尼,K为转子系统的刚度,G为转子系统的陀螺力矩矩阵,为84×84矩阵;Ω为给定的转速,Ω=7500r/min; q1为压气机涡轮转子系统21节点的不平衡响应,描述每个节点的响应,每个节点有4个自由度,包含2个平动、2个转动自由度,整个转子系统具有84个自由度;Q为作用在转子系统上的残余不平衡向量:
(2)
其中Qi用来描述压气机转子位于第i个节点的残余不平衡,i=2,Qj用来描述涡轮转子位于第j个节点的残余不平衡,j=20,此时涡轮转子4残余不平衡的相位相对压气机转子2残余不平衡量相位之间的差值为θ1,θ1=0°有:
(3)
(4)
e为自然对数的底数,i为虚数单位,则压气机涡轮转子系统的各节点的不平衡响应q1的计算公式为:
(5)
压气机涡轮转子系统的各节点的不平衡响应q1计算结果见表1。
表1 涡轮转子4残余不平衡的相位相对压气机转子2残余不平衡量相位之间的差值θ1 = 0º时,压气机涡轮转子系统的各节点的不平衡响应q1的计算结果
不考虑压气机涡轮转子系统的各节点的不平衡响应q中的转动响应,则各个节点的最大不平衡响应R1
R1=max(q1)(6)
通过表1确定1号节点即转子系统前支点处,振动响应最大,为3.34mm/s,将这组状态记作(R1,θ1),其中R1=3.34、θ1=0º;
步骤4,以压气机转子2为基准,涡轮转子4按照逆时针方向相对压气机转子转动1个花键套齿,此状态下的涡轮转子4残余不平衡量与压气机转子残余不平衡量相位之间的差值记作θ2,θ2=22.5°,将此状态下的压气机涡轮转子系统的各节点的不平衡响应记作q2,将步骤3中的θ1换为θ2,q1换为q2,重复步骤3计算该状态下压气机涡轮转子系统的各节点的不平衡响应q2,计算结果见表2。
表2 涡轮转子4残余不平衡的相位相对压气机转子2残余不平衡量相位之间的差值θ2=22.5°时,压气机涡轮转子系统的各节点的不平衡响应q2的计算结果
通过表2确定1号节点即转子前支点处,振动响应最大,为3.32mm/s。
将这组状态记作(R2,θ2),其中R2=3.32、θ2=22.5º;
步骤5,重复步骤4,得到全部16中状态下的最大不平衡响应,分别为(R1,θ1)、(R2,θ2)……(R16,θ16),见表3;
表3 不同涡轮转子残余不平衡量与压气机转子残余不平衡量相位之间的差值θ下的最大不平衡响应
步骤6,选出(R1,θ1)、(R2,θ2)……(R16,θ16)中最大不平衡响应R的最小值Rmin :
Rmin =min(R1,R2,...,Rm)(7)
选出第9中安装状态为最小值Rmin,即θ=180°,Rmin=2.71mm/s。
确定最终装配方案为(2.71,180 º),即涡轮转子和压气机转子的残余不平衡量相位差为180 º的状态下,燃机振动最小,为2.71mm/s。

Claims (1)

1.一种燃气轮机的压气机涡轮转子的角向固定对接组合方法,其特征在于,包括以下步骤:
步骤1,压气机转子在动平衡机上完成单独动平衡的残余不平衡量记作m1,相位记作α1,涡轮转子在动平衡机上完成单独动平衡的残余不平衡量记作m2,相位记作α2;其中相位表示的是残余不平衡量相对确定好的参考标记的角度;
步骤2,建立压气机涡轮转子的转子动力学有限元模型;
步骤3,将涡轮转子残余不平衡量与压气机转子残余不平衡量相位之间的任意差值用θ表示,压气机涡轮转子系统的各节点的不平衡响应用q表示,压气机涡轮转子系统各个节点的最大不平衡响应用R表示;
将压气机转子与涡轮转子角向对接,此时以压气机转子残余不平衡量的相位α1为基准,将此状态下的涡轮转子残余不平衡量与压气机转子残余不平衡量相位之间的差值记作θ1,将此状态下的压气机涡轮转子系统的各节点的不平衡响应记作q1,利用转子动力学方程计算该状态下的各个节点的最大不平衡响应R1
压气机涡轮转子系统的转子动力学方程为:
(1)
其中M 、C 、K 、G 分别为转子系统的质量、阻尼、刚度以及陀螺力矩矩阵;Ω为给定的转速;q1为压气机涡轮转子系统的各节点的不平衡响应,描述每个节点的响应,每个节点有4个自由度,包含2个平动、2个转动自由度,整个转子系统具有4n个自由度,n为压气机涡轮转子模型的节点数;Q为作用在转子系统上的残余不平衡向量:
(2)
其中Qi用来描述压气机转子位于第i个节点的残余不平衡,Qj用来描述涡轮转子位于第j个节点的残余不平衡,此时涡轮转子残余不平衡的相位相对压气机转子残余不平衡量相位之间的差值为θ1,有:
(3)
(4)
则压气机涡轮转子系统的各节点的不平衡响应q1的计算公式为:
(5)
不考虑压气机涡轮转子系统的各节点的不平衡响应q中的转动响应,则各个节点的最大不平衡响应R1
R1=max(q1)(6)
将这组状态记作(R1,θ1);
步骤4,以压气机转子为基准,涡轮转子按照逆时针方向相对压气机转子转动1个花键套齿,此状态下的涡轮转子残余不平衡量与压气机转子残余不平衡量相位之间的差值记作θ2,将此状态下的压气机涡轮转子系统的各节点的不平衡响应记作q2,重复步骤3计算该状态下各个节点的最大不平衡响应R2,将这组状态记作(R2,θ2);
步骤5,重复步骤4,得到全部状态下的最大不平衡响应,分别为(R1,θ1)、(R2,θ2)……(Rm,θm);
步骤6,选出(R1,θ1)、(R2,θ2)……(Rm,θm)中最大不平衡响应R的最小值Rmin :
Rmin =min(R1,R2,...,Rm)(7)
此时Rmin对应的安装角度差记为θmin,相应的对接状态(Rmin,θmin)作为最终装配方案。
CN202111519726.2A 2021-12-14 2021-12-14 一种燃气轮机的压气机涡轮转子的角向固定对接组合方法 Active CN114201879B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202111519726.2A CN114201879B (zh) 2021-12-14 2021-12-14 一种燃气轮机的压气机涡轮转子的角向固定对接组合方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202111519726.2A CN114201879B (zh) 2021-12-14 2021-12-14 一种燃气轮机的压气机涡轮转子的角向固定对接组合方法

Publications (2)

Publication Number Publication Date
CN114201879A CN114201879A (zh) 2022-03-18
CN114201879B true CN114201879B (zh) 2024-04-09

Family

ID=80653078

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202111519726.2A Active CN114201879B (zh) 2021-12-14 2021-12-14 一种燃气轮机的压气机涡轮转子的角向固定对接组合方法

Country Status (1)

Country Link
CN (1) CN114201879B (zh)

Families Citing this family (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN115541116B (zh) * 2022-11-24 2023-10-10 中国航发沈阳发动机研究所 一种减小涡轮转子旋转惯性激励的控制方法

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
KR100743609B1 (ko) * 2006-06-20 2007-07-27 한국전력공사 회전기기 블레이드의 배열위치 산출방법
CN105550395A (zh) * 2015-12-04 2016-05-04 湖南科技大学 透平机械单支撑轴系几何对中安装方法
CN106224016A (zh) * 2016-08-31 2016-12-14 中国南方航空工业(集团)有限公司 涡轮转子、发动机及提高发动机转子对中可靠性的方法
CN109690024A (zh) * 2016-09-13 2019-04-26 西门子股份公司 用于平衡燃气轮机的压气机的转子的技术
CN113092123A (zh) * 2021-03-25 2021-07-09 中国航发南方工业有限公司 燃气发生器转子的装配方法

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
KR100743609B1 (ko) * 2006-06-20 2007-07-27 한국전력공사 회전기기 블레이드의 배열위치 산출방법
CN105550395A (zh) * 2015-12-04 2016-05-04 湖南科技大学 透平机械单支撑轴系几何对中安装方法
CN106224016A (zh) * 2016-08-31 2016-12-14 中国南方航空工业(集团)有限公司 涡轮转子、发动机及提高发动机转子对中可靠性的方法
CN109690024A (zh) * 2016-09-13 2019-04-26 西门子股份公司 用于平衡燃气轮机的压气机的转子的技术
CN113092123A (zh) * 2021-03-25 2021-07-09 中国航发南方工业有限公司 燃气发生器转子的装配方法

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
多转子轴系优化耦联及其不平衡响应有限元分析;高金吉, 缪红燕, 徐鸿, 王维民;振动与冲击;20050630(第02期);全文 *
面向核心机单元体独立交付的总体结构设计;潘辉;马健;林磊;宋德福;朱成;王丹丹;;燃气涡轮试验与研究;20200415(第02期);全文 *

Also Published As

Publication number Publication date
CN114201879A (zh) 2022-03-18

Similar Documents

Publication Publication Date Title
CN114201879B (zh) 一种燃气轮机的压气机涡轮转子的角向固定对接组合方法
Liu et al. Field dynamic balancing for rigid rotor-AMB system in a magnetically suspended flywheel
CN107341288B (zh) 一种通过调整轴承标高控制联合循环机组振动的优化方法
JP2010506113A (ja) 目標角度制御式クラッチ連結方法
CN110457740B (zh) 一种机械结构在基础激励下的响应特性分析方法
CN109632188B (zh) 虑及多转速状态的柔性转子动平衡双重优化方法及系统
CN106121907A (zh) 叶片配平衡方法、叶片及风力发电机
CN111339610B (zh) 一种叶轮机械转子叶片装配优化排序方法
CN111209639A (zh) 一种叶轮-轴承-转子系统的高效定量建模方法
CN113341714B (zh) 一种磁悬浮轴承转子控制系统同频干扰的抵消方法
CN112556931B (zh) 基于粒子群算法的高速轴承转子系统模态动平衡方法
CN1067764C (zh) 旋转机械动平衡的动刚度系数平衡方法及平衡装置
CN110750932A (zh) 叶盘-机匣系统的碰摩动力学特性的数字仿真模拟方法
Genta et al. DYNROT: a finite element code for rotordynamic analysis based on complex co‐ordinates
Luo et al. Stability of periodic motion on the rotor-bearing system with coupling faults of crack and rub-impact
Wang et al. The research on temperature effect for dynamic performance of high-speed permanent-magnet synchronous machines rotors
CN109424436B (zh) 用于操作一组发电机的方法和系统
Cao Transient analysis of flexible rotors with nonlinear bearings, dampers and external forces
CN113326647A (zh) 一种水轮机轴系转子动力学模态计算方法
Kamesh Oil-whirl instability in an automotive turbocharger
Lyu et al. Research on Rotor Dynamics of a 60000rpm 200KW High Speed Motor
CN117650583B (zh) 一种水电站一管多机并网多目标协调优化控制方法及系统
CN106777635B (zh) 一种在风机参数及电机运转参数已知条件下的空冷岛风机桥架优化方法
CN116090135B (zh) 一种阻尼器-转子系统响应分析方法
Wang et al. Improve the stability of rotor subjected to fluid leakage by optimum diameters design

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