CN112464529A - 柔性支撑汽轮发电机组运行状态下载荷分配实时计算方法 - Google Patents

柔性支撑汽轮发电机组运行状态下载荷分配实时计算方法 Download PDF

Info

Publication number
CN112464529A
CN112464529A CN202011276103.2A CN202011276103A CN112464529A CN 112464529 A CN112464529 A CN 112464529A CN 202011276103 A CN202011276103 A CN 202011276103A CN 112464529 A CN112464529 A CN 112464529A
Authority
CN
China
Prior art keywords
load distribution
bearing
calculating
steam turbine
flexible support
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.)
Pending
Application number
CN202011276103.2A
Other languages
English (en)
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 Thermal Power Research Institute Co Ltd
Original Assignee
Xian Thermal Power Research Institute 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 Xian Thermal Power Research Institute Co Ltd filed Critical Xian Thermal Power Research Institute Co Ltd
Priority to CN202011276103.2A priority Critical patent/CN112464529A/zh
Publication of CN112464529A publication Critical patent/CN112464529A/zh
Pending legal-status Critical Current

Links

Images

Classifications

    • 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)
  • Testing Of Devices, Machine Parts, Or Other Structures Thereof (AREA)

Abstract

柔性支撑汽轮发电机组运行状态下载荷分配实时计算方法,1)使用有限元法建立轴系及轴承模型,使用轴承原始标高计算刚性支撑静载荷分配,在此基础上考虑弹簧弹性及油膜厚度影响,迭代计算得到柔性基础汽轮发电机组运行状态下的载荷分配;2)计算得到运行状态下柔性支撑载荷分配与刚性支撑静载荷分配的影响系数关系曲线;3)按照曲线可在实际运行时实时计算得到当前转速下的轴系载荷分配。为更精确的计算运行状态下的载荷分配必须考虑到柔性基础本身的变形量以及轴承油膜的厚度,其关键技术就是柔性支撑、运行状态下载荷分配影响系数曲线的确定。

Description

柔性支撑汽轮发电机组运行状态下载荷分配实时计算方法
技术领域
本发明涉及旋转机械状态监测技术领域,特别涉及柔性支撑汽轮发电机组运行状态下载荷分配实时计算方法。
背景技术
目前国内大多数汽轮发电机组均使用刚性支撑基础,由于刚性支撑刚度较大,支撑本身的变化量在载荷分配计算中可忽略不计,因此其载荷分配容易计算。近年来,随着柔性基础在大容量火电、核电机组上的广泛应用,柔性基础带来了轴系载荷分配的变化,同时运行状态下载荷分配一直没有特别准确有效的测量方法。为更精确的计算运行状态下的载荷分配,辅助判断轴系运行状态,必须考虑到柔性基础本身的变形量以及轴承油膜的厚度,其关键技术就是柔性支撑、运行状态下载荷分配影响系数曲线的确定。
发明内容
为了克服以上技术问题,本发明的目的在于提供一种柔性支撑汽轮发电机组运行状态下载荷分配实时计算方法,为使用柔性基础的汽轮发电机组轴系运行状态监测提供数据支持。
为了实现上述目的,本发明采用的技术方案是:
一种柔性支撑汽轮发电机组运行状态下载荷分配实时计算方法,包括以下步骤;
1)使用有限元法建立轴系及轴承模型,使用轴承原始标高计算刚性支撑静载荷分配,在此基础上考虑弹簧弹性及油膜厚度影响,迭代计算得到柔性基础汽轮发电机组运行状态下的载荷分配;
2)计算得到运行状态下柔性支撑载荷分配与刚性支撑静载荷分配的影响系数关系曲线;
3)按照曲线可在实际运行时实时计算得到当前转速下的轴系载荷分配。
所述步骤1)具体步骤如下:
1)按照转子质量和抗弯刚度不变的等效原则,按有限元法将大型汽轮发电机组简化为集总质量模型,确定各轴承处原始标高hi (0),使用有限元法根据轴段的力平衡和变形关系得到刚性基础静止状态下的各轴承载荷分配Wi (0)(i是轴承编号);
2)按照Wi (0)的静载荷分配计算各轴承处垂直方向油膜厚度与转子转速关系曲线,确定不同转速下第i个轴承处垂直方向油膜厚度
Figure BDA0002779118010000022
(r)(i是轴承编号,r是转速,其取值从垂直方向油膜厚度与转子转速关系曲线图中读取);
3)根据Wi (0)及相应轴承处支撑刚度ki的大小确定第i个轴承处弹簧伸缩量li (1),使得:
Figure BDA0002779118010000023
4)使用计算得到的hi (1)(r)代替hi (0)后重复步骤1),得到新的载荷分配Wi (1)(r);
5)比较Wi (n)(n)与Wi (n-1)(n)之间的差别,当
Figure BDA0002779118010000021
时,使用
Figure BDA0002779118010000033
(r)和Wi (n)(r)代替
Figure BDA0002779118010000034
(r)和Wi (n-1)(r)并重复步骤1)至步骤4),当
Figure BDA0002779118010000031
时,Wi (n)(r)即柔性支撑、运行状态下的计算真实载荷。
所述步骤2)具体步骤如下:
计算运行状态下柔性支撑与刚性支撑间载荷分配的影响系数曲线:
Figure BDA0002779118010000032
所述步骤3)具体步骤如下:
针对具体机组,按照步骤1)和步骤2)计算得到其刚性支撑状态下的轴系静载荷分配,以及运行状态下柔性支撑与刚性支撑间载荷分配的影响系数曲线,即可在实际运行时实时计算得到当前转速下的轴系载荷分配。
Wi (n)(r)=Wi (0)×αi(n)。
本发明的有益效果:
1)当前载荷分配主要是依靠静态下测量计算,但运行状态与静止状态的载荷分配存在显著差别,使本发明可以计算使用柔性基础的机组运行状态下的实际载荷分配;
2)运行状态下的轴系载荷分配需要实测运行时的标高变化,测量成本巨大,且受测点安装因素影响存在误差,本发明是通过精确计算得到运行状态下的载荷分配,有效减少了测量成本;
3)本方法实现了载荷分配的实时计算,可以为机组轴系运行安全状态的检验分析提供及时的支持。
附图说明:
图1为本发明有限元法将大型汽轮发电机组简化为集总质量模型。
图2为运行状态下柔性支撑载荷分配与刚性支撑静载荷分配的影响系数关系曲线。
具体实施方式
下面结合附图对本发明作进一步详细说明。
柔性支撑汽轮发电机组运行状态下载荷分配实时计算方法分成柔性基础汽轮发电机组运行状态下的载荷分配的计算、影响系数关系曲线的确定和运行状态下轴系载荷分配的确定三个步骤,三个步骤一起共同实现准确、高效的现场动平衡。具体如下:
1、使用有限元法建立轴系及轴承模型,使用轴承原始标高计算刚性支撑静载荷分配,在此基础上考虑弹簧弹性及油膜厚度影响,迭代计算得到柔性基础汽轮发电机组运行状态下的载荷分配;
1)按照转子质量和抗弯刚度不变的等效原则,按有限元法将大型汽轮发电机组简化为集总质量模型(如附图1),确定各轴承处原始标高hi (0),使用有限元法根据轴段的力平衡和变形关系得到刚性基础静止状态下的各轴承载荷分配Wi (0)(i是轴承编号);
2)按照Wi (0)的静载荷分配计算各轴承处垂直方向油膜厚度与转子转速关系曲线,确定不同转速下第i个轴承处垂直方向油膜厚度
Figure BDA0002779118010000041
(r)(i是轴承编号,r是转速,其取值从垂直方向油膜厚度与转于转速关系曲线图中读取);
3)根据Wi (0)及相应轴承处支撑刚度ki的大小确定第i个轴承处弹簧伸缩量li (1),使得:
Figure BDA0002779118010000054
4)使用计算得到的hi (1)(r)代替hi (0)后重复步骤1),得到新的载荷分配Wi (1)(r);
5)比较Wi (n)(n)与Wi (n-1)(n)之间的差别,当
Figure BDA0002779118010000051
时,使用
Figure BDA0002779118010000055
(r)和Wi (n)(r)代替
Figure BDA0002779118010000056
(r)和Wi (n-1)(r)并重复步骤1)至步骤4),当
Figure BDA0002779118010000052
时,Wi (n)(r)即柔性支撑、运行状态下的计算真实载荷;
3、计算得到运行状态下柔性支撑载荷分配与刚性支撑静载荷分配的影响系数关系曲线(如附图2):
计算运行状态下柔性支撑与刚性支撑间载荷分配的影响系数曲线:
Figure BDA0002779118010000053
4、按照曲线可在实际运行时实时计算得到当前转速下的轴系载荷分配:
针对具体机组,按照步骤1)和步骤2)计算得到其刚性支撑状态下的轴系静载荷分配,以及运行状态下柔性支撑与刚性支撑间载荷分配的影响系数曲线,即可在实际运行时实时计算得到当前转速下的轴系载荷分配。
Wi (n)(r)=Wi (0)×αi(n)。

Claims (4)

1.一种柔性支撑汽轮发电机组运行状态下载荷分配实时计算方法,其特征在于,包括以下步骤;
1)使用有限元法建立轴系及轴承模型,使用轴承原始标高计算刚性支撑静载荷分配,在此基础上考虑弹簧弹性及油膜厚度影响,迭代计算得到柔性基础汽轮发电机组运行状态下的载荷分配;
2)计算得到运行状态下柔性支撑载荷分配与刚性支撑静载荷分配的影响系数关系曲线;
3)按照曲线可在实际运行时实时计算得到当前转速下的轴系载荷分配。
2.根据权利要求1所述的一种柔性支撑汽轮发电机组运行状态下载荷分配实时计算方法,其特征在于,所述步骤1)具体步骤如下:
1)按照转子质量和抗弯刚度不变的等效原则,按有限元法将大型汽轮发电机组简化为集总质量模型,确定各轴承处原始标高hi (0),使用有限元法根据轴段的力平衡和变形关系得到刚性基础静止状态下的各轴承载荷分配Wi (0)(i是轴承编号);
2)按照Wi (0)的静载荷分配计算各轴承处垂直方向油膜厚度与转子转速关系曲线,确定不同转速下第i个轴承处垂直方向油膜厚度
Figure FDA0002779115000000011
(i是轴承编号,r是转速,其取值从垂直方向油膜厚度与转子转速关系曲线图中读取);
3)根据Wi (0)及相应轴承处支撑刚度ki的大小确定第i个轴承处弹簧伸缩量li (1),使得:
Figure FDA0002779115000000021
4)使用计算得到的hi (1)(r)代替hi (0)后重复步骤1),得到新的载荷分配Wi (1)(r);
5)比较Wi (n)(n)与Wi (n-1)(n)之间的差别,当
Figure FDA0002779115000000022
时,使用
Figure FDA0002779115000000023
和Wi (n)(r)代替
Figure FDA0002779115000000024
和Wi (n-1)(r)并重复步骤1)至步骤4),当
Figure FDA0002779115000000025
时,Wi (n)(r)即柔性支撑、运行状态下的计算真实载荷。
3.根据权利要求1所述的一种柔性支撑汽轮发电机组运行状态下载荷分配实时计算方法,其特征在于,所述步骤2)具体步骤如下:
计算运行状态下柔性支撑与刚性支撑间载荷分配的影响系数曲线:
Figure FDA0002779115000000026
4.根据权利要求1所述的一种柔性支撑汽轮发电机组运行状态下载荷分配实时计算方法,其特征在于,所述步骤3)具体步骤如下:
针对具体机组,按照步骤1)和步骤2)计算得到其刚性支撑状态下的轴系静载荷分配,以及运行状态下柔性支撑与刚性支撑间载荷分配的影响系数曲线,即可在实际运行时实时计算得到当前转速下的轴系载荷分配;
Wi (n)(r)=Wi (0)×αi(n)。
CN202011276103.2A 2020-11-16 2020-11-16 柔性支撑汽轮发电机组运行状态下载荷分配实时计算方法 Pending CN112464529A (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202011276103.2A CN112464529A (zh) 2020-11-16 2020-11-16 柔性支撑汽轮发电机组运行状态下载荷分配实时计算方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202011276103.2A CN112464529A (zh) 2020-11-16 2020-11-16 柔性支撑汽轮发电机组运行状态下载荷分配实时计算方法

Publications (1)

Publication Number Publication Date
CN112464529A true CN112464529A (zh) 2021-03-09

Family

ID=74837048

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202011276103.2A Pending CN112464529A (zh) 2020-11-16 2020-11-16 柔性支撑汽轮发电机组运行状态下载荷分配实时计算方法

Country Status (1)

Country Link
CN (1) CN112464529A (zh)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN113847981A (zh) * 2021-09-16 2021-12-28 国家电网有限公司 一种基于机械特性的水电机组保护性振动监测方法

Citations (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103306302A (zh) * 2012-03-07 2013-09-18 江苏核电有限公司 一种计入弹簧基础的轴系载荷预调整方法

Patent Citations (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103306302A (zh) * 2012-03-07 2013-09-18 江苏核电有限公司 一种计入弹簧基础的轴系载荷预调整方法

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
何国安 等: "计入弹簧基础的1060MW 核电机组轴承载荷灵敏度分析", 《机械工程师》, 31 October 2010 (2010-10-31) *
刘荣强 等: "轴承标高对多跨轴系振动及稳定性的影响", 《哈尔滨工业大学学报》, 28 February 1995 (1995-02-28) *

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN113847981A (zh) * 2021-09-16 2021-12-28 国家电网有限公司 一种基于机械特性的水电机组保护性振动监测方法

Similar Documents

Publication Publication Date Title
CN109340062B (zh) 一种低风速风电机组数字双胞胎式的疲劳损伤预测方法
CN101576425B (zh) 一种汽轮机叶轮静平衡的测量方法与装置
Gong et al. Simulation investigation of wind turbine imbalance faults
CN109635438B (zh) 一种水轮机顶盖振动数值的确定方法
CN112685939B (zh) 一种基于实测的海上风力机基础疲劳损伤分析方法
CN102032854B (zh) 风电主轴承的测量装置
CN112464529A (zh) 柔性支撑汽轮发电机组运行状态下载荷分配实时计算方法
Nässelqvist et al. A methodology for protective vibration monitoring of hydropower units based on the mechanical properties
CN109684711A (zh) 一种涡轴发动机气动连接多转子耦合振动分析方法
CN114676577A (zh) 基于多物理场数字孪生模型的空气静压主轴模拟方法及系统
Chang et al. Flow field computation of the NREL S809 airfoil using various turbulence models
Cervantes et al. Porjus U9A full-scale hydropower research facility
CN116191478B (zh) 一种多个异步电机等效惯量评估和频率响应建模方法
CN215767662U (zh) 一种径向空气轴承测试装置
CN112557037B (zh) 一种风电主轴轴承模型加速寿命试验方法
Li et al. Axial Location Analysis of 1000 MW Unit Low-Pressure Rotor Rubbing Fault Based on ANSYS
CN113153453B (zh) 汽轮机末级叶片容积流量估计方法、颤振预警方法及系统和装置
CN107869959A (zh) 一种核电汽轮机低压转子膨胀量激光测量方法
Wang et al. A new fault detection strategy for wind turbine rotor imbalance based on multi-condition vibration signal analysis
Simms et al. Characterization and comparison of baseline aerodynamic performance of optimally twisted versus non-twisted HAWT blades
Dai et al. Field load testing of wind turbines based on the relational model of strain vs load
Akins et al. Pressure distributions on an operating vertical-axis wind-turbine blade element
Synnegård Dynamic modelling and analysis of vertical machines
Reinhardt et al. Laser vibrometry measurements of rotating blade vibrations
Stirling Internal load modelling of tapered-roller main bearings in wind turbines

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