CN106503375B - 一种基于cn群理论确定汽轮机转子临界转速的方法及系统 - Google Patents

一种基于cn群理论确定汽轮机转子临界转速的方法及系统 Download PDF

Info

Publication number
CN106503375B
CN106503375B CN201610966781.9A CN201610966781A CN106503375B CN 106503375 B CN106503375 B CN 106503375B CN 201610966781 A CN201610966781 A CN 201610966781A CN 106503375 B CN106503375 B CN 106503375B
Authority
CN
China
Prior art keywords
turbine rotor
basic sector
rotating speed
steam turbine
dynamic
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.)
Expired - Fee Related
Application number
CN201610966781.9A
Other languages
English (en)
Other versions
CN106503375A (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.)
Shandong University
Original Assignee
Shandong 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 Shandong University filed Critical Shandong University
Priority to CN201610966781.9A priority Critical patent/CN106503375B/zh
Publication of CN106503375A publication Critical patent/CN106503375A/zh
Application granted granted Critical
Publication of CN106503375B publication Critical patent/CN106503375B/zh
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

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
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F2111/00—Details relating to CAD techniques
    • G06F2111/10—Numerical modelling

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

本发明公开了一种基于CN群理论确定汽轮机转子临界转速的方法及系统，通过建立汽轮机转子的3D几何模型，选取3D几何模型的一部分作为基本扇区；建立基本扇区的局部坐标系，对基本扇区进行离散，计算离散后的各单元的动特性参数，根据CN群理论和基本扇区的动特性参数，建立汽轮机转子的动力学模型，根据动力学模型计算得到汽轮机转子的各阶固有频率，确定临界转速。本发明操作流程清晰、计算效率高、结果准确，适用于工程设计、试验及运行中快速找准汽轮机转子临界转速。

Description

一种基于CN群理论确定汽轮机转子临界转速的方法及系统
技术领域
本发明涉及一种基于CN群理论确定汽轮机转子临界转速的方法及系统。
背景技术
现有的获取汽轮机转子临界转速方法有两种：计算法和实测法。实测法虽能准确测出转子的临界转速，但存在很大的局限性，如在临界转速大于工作转速的情况下，就无法实测得到转子的临界转速。计算法多采用对实际汽轮机转子的几何形状、支承刚度等过度简化的模型，导致计算值与实测值有一定的误差；而采用转子真实模型计算规模巨大，求解效率低。
关于汽轮机转子临界转速计算方法的报道，多出现于论文中，暂未检索到相关专利。如 2015年第9届国际转子动力学会议公开了《The Influence of Blade Row Dynamicson Lateral and Torsional Shaft Vibrations in Steam Turbines》，其将汽轮机转子中的叶片简化为位于叶片顶端的质量单元，相邻叶片间的连接及叶根处的连接以弹簧等效代替，该简化模型能高效地计算出叶片振动与轴振动间的相互影响规律，但无法反映出叶片质量在叶片长度及宽度方向上的分布，导致汽轮机转子临界转速的计算值与真实值间有较大误差。《动力工程学报》2014 年第6期公开了《汽轮机长叶片弯曲与轴系扭转耦合振动研究》，其基于汽轮机转子实际结构参数建立三维模型，采用全环模型进行计算，能较为准确地得到汽轮机转子的各阶临界转速，但求解规模过于巨大、效率低。
除此之外，2014年ASME涡轮技术峰会公开了《Last Stage Blade Coupled ShaftTorsional Vibration Analysis of 1000MW Steam Turbine Generator Set by aReduced 3D Finite Element Method》，其基于汽轮机转子实际结构参数建立三维模型，并采用子结构模态综合法提取重要的转子低阶振型，既保证计算的准确性又提高计算效率，但该方法对技术人员理论基础及编程能力要求较高，有一定的操作难度。
综上，现有的获取汽轮机转子临界转速方法存在实测法应用受限、计算法与真实值有一定误差或计算难度大等问题，使得计算法难以用于含多级叶片汽轮机转子。
发明内容
本发明为了解决上述问题，提出了一种基于CN群理论确定汽轮机转子临界转速的方法及系统，本方法利用汽轮机转子的循环对称性，在保证计算精度的同时降低计算规模，能较为高效准确地获得汽轮机转子的临界转速。
为了实现上述目的，本发明采用如下技术方案：
一种基于CN群理论确定汽轮机转子临界转速的方法，包括以下步骤：
(1)建立汽轮机转子的3D几何模型，选取3D几何模型的一部分作为基本扇区；
(2)建立基本扇区的局部坐标系，对基本扇区进行离散，计算离散后的各单元的动特性参数；
(3)根据CN群理论和基本扇区的动特性参数，建立汽轮机转子的动力学模型，根据动力学模型计算得到汽轮机转子的各阶固有频率，确定临界转速。
所述步骤(1)中，根据轴段、轮盘、叶片的结构参数，建立汽轮机转子的3D几何模型。
所述步骤(1)中，确定汽轮机转子中各级叶片的最大公约数N，在圆周方向取3D模型的1/N个连续的叶片形成的区域作为基本扇区。
所述步骤(2)中，采用有限单元法对基本扇区进行离散。
所述步骤(2)中，对各离散单元的动特性参数进行组装，得到基本扇区的质量矩阵、刚度矩阵和阻尼矩阵。
所述步骤(3)中，根据CN群理论和基本扇区的动特性参数，建立汽轮机转子的动力学模型，得到用基本扇区动特性参数表征的汽轮机转子的特征方程。
所述步骤(3)中，采用Newmark法求解特征方程，得到汽轮机转子的临界转速。
一种基于上述方法的系统，包括：
建模模块，被配置为根据轴段、轮盘、叶片的结构参数，建立汽轮机转子的3D几何模型，选取3D几何模型的一部分作为基本扇区；
离散模块，被配置为建立基本扇区的局部坐标系，对基本扇区进行离散，计算离散后的各单元的动特性参数；
处理模块，被配置为根据CN群理论和基本扇区的动特性参数，建立汽轮机转子的动力学模型，根据动力学模型计算得到汽轮机转子的各阶固有频率，确定临界转速。
本发明的有益效果为：
(1)本发明根据汽轮机转子的实际特征数据，建立较为精确的汽轮机转子几何模型，基于CN群的理论和基本扇区的动特性参数，得到用基本扇区动特性参数表征的较低阶的汽轮机转子特征方程，求解特征方程可得汽轮机转子的各阶固有频率，即临界转速；
(2)本发明避免了采用汽轮机转子整体计算规模巨大，采用过度简化的汽轮机转子模型计算误差大的问题，同时克服了实测法应用的局限性，具有准确度高、适用范围广的优点，可为机组的运行或高速动平衡提供可靠依据；
(3)本发明利用汽轮机转子的循环对称性，在保证计算精度的同时降低计算规模，操作流程清晰、计算效率高、结果准确，适用于工程设计、试验及运行中快速找准汽轮机转子临界转速。
附图说明
图1为本发明的基于CN群理论计算汽轮机转子临界转速的方法的流程图；
图2为本发明的一个具体实施例的汽轮机转子几何模型；
图3为本发明的一个具体实施例的汽轮机转子基本扇区划分示意图；
图4为本发明的一个具体实施例的汽轮机转子基本扇区离散模型；
其中Ⅰ-Ⅹ分别表示汽轮机转子的第Ⅰ-Ⅹ级叶片。
具体实施方式：
下面结合附图与实施例对本发明作进一步说明。
如图1所示，本发明的基于CN群理论计算汽轮机转子临界转速的方法主要通过三大模块完成：基本扇区创建模块、基本扇区动特性参数计算模块、动力学建模及求解模块。各模块之间依次进行数据处理，首先在圆周方向取汽轮机转子几何模型的1/N作基本扇区，然后用有限单元法对基本扇区离散得到其动特性参数；最后求解用基本扇区动特性参数表征的较低阶的汽轮机转子特征方程，得到所需各阶临界转速。
所述各模块进行数据处理是通过以下步骤实现的：
1)划分基本扇区
1-1)如图2所示，根据轴段、轮盘、叶片的结构参数，建立汽轮机转子的3D几何模型；
1-2)图2中的汽轮机转子在轴向关于中部平面对称，为降低求解规模，取几何模型的一半进行计算，同时各级叶片的最大公约数为20，在圆周方向取1/2模型的1/20，作为基本扇区，如第Ⅹ级包含60个叶片，则取含有3个叶片的扇区作基本扇区，如图3所示。
2)基本扇区动特性参数计算
2-1)如图4所示，建立基本扇区的局部坐标系oi-xiyizi，采用有限单元法对基本扇区进行离散；
2-2)记基本扇区两侧界面分别为l、r，两界面内部区域为g，则基本扇区的质量矩阵M、阻尼矩阵C和刚度矩阵K分别为
Figure BDA0001145836980000041
Figure BDA0001145836980000042
Figure BDA0001145836980000043
3)特征方程及固有频率求解
3-1)根据CN群理论，汽轮机转子的临界转速可通过求解用基本扇区动特性参数表征的特征方程得到，进而降低计算规模，所述特征方程为：
(λ2αp+λβp+γp)ψp＝0
其中：
Figure BDA0001145836980000044
Figure BDA0001145836980000045
Figure BDA0001145836980000046
Ψp为特征值λ对应的特征向量，p为汽轮机转子振动模态的节径数，p＝0，1，2，…，N/2，α＝2π/N；
3-2)将基本扇区的动特性参数代入所述特征方程，并采用Newmark法求解所述特征方程，可得汽轮机转子的各阶固有频率，即临界转速。
本发明基于CN群理论和有限单元法，利用汽轮机转子结构的循环对称性，用基本扇区的动特性参数来表征汽轮机转子整体的特征方程，在保证计算精度的同时能较大地提升计算效率，另外还可避免实测法在应用范围上的局限性。
上述虽然结合附图对本发明的具体实施方式进行了描述，但并非对本发明保护范围的限制，所属领域技术人员应该明白，在本发明的技术方案的基础上，本领域技术人员不需要付出创造性劳动即可做出的各种修改或变形仍在本发明的保护范围以内。

Claims (5)

1.一种基于CN群理论确定汽轮机转子临界转速的方法，其特征是：包括以下步骤：
(1)建立汽轮机转子的3D几何模型，选取3D几何模型的一部分作为基本扇区；
(2)建立基本扇区的局部坐标系，对基本扇区进行离散，计算离散后的各单元的动特性参数；
(3)根据CN群理论和基本扇区的动特性参数，建立汽轮机转子的动力学模型，根据动力学模型计算得到汽轮机转子的各阶固有频率，确定临界转速；
所述步骤(1)中，确定汽轮机转子中各级叶片的最大公约数N，在圆周方向取3D模型的1/N个连续的叶片形成的区域作为基本扇区；
所述步骤(2)中，对各离散单元的动特性参数进行组装，得到基本扇区的质量矩阵、刚度矩阵和阻尼矩阵；
所述步骤(3)中，根据CN群理论和基本扇区的动特性参数，建立汽轮机转子的动力学模型，得到用基本扇区动特性参数表征的汽轮机转子的特征方程。
2.如权利要求1所述的一种基于CN群理论确定汽轮机转子临界转速的方法，其特征是：所述步骤(1)中，根据轴段、轮盘、叶片的结构参数，建立汽轮机转子的3D几何模型。
3.如权利要求1所述的一种基于CN群理论确定汽轮机转子临界转速的方法，其特征是：所述步骤(2)中，采用有限单元法对基本扇区进行离散。
4.如权利要求1所述的一种基于CN群理论确定汽轮机转子临界转速的方法，其特征是：所述步骤(3)中，采用Newmark法求解特征方程，得到汽轮机转子的临界转速。
5.一种基于如权利要求1-4中任一项所述的方法的系统，其特征是：包括：
建模模块，被配置为根据轴段、轮盘、叶片的结构参数，建立汽轮机转子的3D几何模型，选取3D几何模型的一部分作为基本扇区；
离散模块，被配置为建立基本扇区的局部坐标系，对基本扇区进行离散，计算离散后的各单元的动特性参数；
处理模块，被配置为根据CN群理论和基本扇区的动特性参数，建立汽轮机转子的动力学模型，根据动力学模型计算得到汽轮机转子的各阶固有频率，确定临界转速。
CN201610966781.9A 2016-10-28 2016-10-28 一种基于cn群理论确定汽轮机转子临界转速的方法及系统 Expired - Fee Related CN106503375B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201610966781.9A CN106503375B (zh) 2016-10-28 2016-10-28 一种基于cn群理论确定汽轮机转子临界转速的方法及系统

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201610966781.9A CN106503375B (zh) 2016-10-28 2016-10-28 一种基于cn群理论确定汽轮机转子临界转速的方法及系统

Publications (2)

Publication Number Publication Date
CN106503375A CN106503375A (zh) 2017-03-15
CN106503375B true CN106503375B (zh) 2020-01-10

Family

ID=58322730

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201610966781.9A Expired - Fee Related CN106503375B (zh) 2016-10-28 2016-10-28 一种基于cn群理论确定汽轮机转子临界转速的方法及系统

Country Status (1)

Country Link
CN (1) CN106503375B (zh)

Families Citing this family (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109583063B (zh) * 2018-11-20 2023-04-18 东北大学 一种风扇转子试验模型的动力学特性相似设计方法
CN111062177B (zh) * 2018-12-29 2021-09-03 山东大学 一种基于围带阻尼的汽轮机转子系统稳定性动态优化方法
CN113761800A (zh) * 2021-09-02 2021-12-07 哈尔滨工业大学 基于临界转速对应的轴系动力学参数模型缩比设计方法

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102508978A (zh) * 2011-11-18 2012-06-20 东南大学 一种基于群论的对称杆系结构的可动性判定方法
CN103632006A (zh) * 2013-12-04 2014-03-12 中国寰球工程公司 一种多级离心泵湿转子临界转速的确定方法
CN103729547A (zh) * 2013-12-12 2014-04-16 清华大学 一种航空发动机转子盘鼓组合界面弯曲刚度计算方法

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102508978A (zh) * 2011-11-18 2012-06-20 东南大学 一种基于群论的对称杆系结构的可动性判定方法
CN103632006A (zh) * 2013-12-04 2014-03-12 中国寰球工程公司 一种多级离心泵湿转子临界转速的确定方法
CN103729547A (zh) * 2013-12-12 2014-04-16 清华大学 一种航空发动机转子盘鼓组合界面弯曲刚度计算方法

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
基于C_N群论的叶片-轴系统动态特性分析;乔鹏等;《汽轮机技术》;20170228;第13-16页 *

Also Published As

Publication number Publication date
CN106503375A (zh) 2017-03-15

Similar Documents

Publication Publication Date Title
CN106528982B (zh) 一种有拉筋和围带的干摩擦阻尼失谐叶片的振动分析方法
CN101122541B (zh) 汽轮机叶片振动试验方法及装置
CN106503375B (zh) 一种基于cn群理论确定汽轮机转子临界转速的方法及系统
CN102938003B (zh) 一种叶轮机械计入错频的气动弹性稳定性数值预测方法
CN109556765B (zh) 一种叶片非接触振动应力测量值换算方法
CN102439301A (zh) 负载测量装置及其方法与程序
Waite et al. Physical understanding and sensitivities of low pressure turbine flutter
CN103245291B (zh) 叶片类零件装配精度检测方法
CN105260563B (zh) 一种叶轮预应力模态的实体与轴对称变维度有限元分析法
CN115859536B (zh) 压气机转子叶片非同步振动锁频数值模拟方法
CN114117849A (zh) 低压涡轮叶/盘转子的叶冠阻尼减振分析方法
Oh et al. Vibration analysis of rotating cantilever beams orienting inward
CN111209639B (zh) 一种叶轮-轴承-转子系统的高效定量建模方法
Zhiquan et al. Structure dynamic analysis of a horizontal axis wind turbine system using a modal analysis method
CN117454721B (zh) 基于数字仿真实验的风电场尾流叠加效应评估方法和介质
Yangui et al. Nonlinear analysis of twisted wind turbine blade
Marten et al. Validation and comparison of a newly developed aeroelastic design code for VAWT
CN109684711A (zh) 一种涡轴发动机气动连接多转子耦合振动分析方法
CN102163263A (zh) 风机叶片振动位移及其威布尔分布拟合方法
CN105547591A (zh) 非对称转子无相位平衡法
Ibrahimbegovic et al. Long-term simulation of wind turbine structure for distributed loading describing long-term wind loads for preliminary design
CN112883555B (zh) 一种风力发电机组叶片各截面多方向载荷提取方法
CN107688705A (zh) 基于有限状态的旋翼系统诱导流场的轴向诱导速度计算方法
CN110094299A (zh) 风电机组的偏航对风自矫正方法和设备
RU2400712C1 (ru) Способ коррекции положения главной центральной оси инерции баллистического объекта

Legal Events

Date Code Title Description
C06 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
CF01 Termination of patent right due to non-payment of annual fee

Granted publication date: 20200110

CF01 Termination of patent right due to non-payment of annual fee