CN109145505B - 一种水轮发电机组运行摆度位移峰最大峰值准确计算方法 - Google Patents
一种水轮发电机组运行摆度位移峰最大峰值准确计算方法 Download PDFInfo
- Publication number
- CN109145505B CN109145505B CN201811095036.7A CN201811095036A CN109145505B CN 109145505 B CN109145505 B CN 109145505B CN 201811095036 A CN201811095036 A CN 201811095036A CN 109145505 B CN109145505 B CN 109145505B
- Authority
- CN
- China
- Prior art keywords
- center
- swing
- axis
- rotation
- maximum
- 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/20—Design optimisation, verification or simulation
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02E—REDUCTION OF GREENHOUSE GAS [GHG] EMISSIONS, RELATED TO ENERGY GENERATION, TRANSMISSION OR DISTRIBUTION
- Y02E10/00—Energy generation through renewable energy sources
- Y02E10/20—Hydro energy
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)
- Measurement Of Mechanical Vibrations Or Ultrasonic Waves (AREA)
Abstract
本发明属于振动摆度技术领域,尤其涉及一种水轮发电机组运行摆度位移峰最大峰值准确计算方法。该方法包括如下步骤:S1、在大轴上安装多对涡流传感器;S2、获取旋转中心围绕固定参考中心相对转动角速度夹角(θ1)和大轴几何中心围绕旋转中心相对转动角速度夹角(θ2),所述固定参考中心为机组中心;S3、根据多个涡流传感器所测量的偏移量,获得大轴摆度值在X轴、Y轴的分量;S4、定义摆度矢量:S5、获得摆度位移峰峰值最大值和相位。本发明在合理假设的前提下,推导出大轴旋转中心动态变化下,大轴轴心轨迹曲线与运行摆度之间关系,给出了一种准确计算水轮机运行摆度位移峰峰值的方法。
Description
技术领域
本发明属于振动摆度技术领域,尤其涉及一种水轮发电机组运行摆度位移峰最大峰值准确计算方法。
背景技术
工程实践中,水轮发电机组大轴的稳定运行是水电机组安全稳定运行的基本要求。《GB 11348.1-1999(ISO 7919-1)旋转机械转轴径向振动的测量和评定第一部分总则》B3.2.1到B3.2.3小节给出了三种方法来计算运行摆度位移峰峰值最大值。
其中,方法1只能作为近视值使用,当同频振动占优势时,一般将过高估计S(P-P)max值,最大误差约为40%。对于圆形轨迹误差最大,当轨迹变平时误差逐渐减小,当退化成直线轨迹时误差为零。
方法2也只能作为近似值使用,当同频振动占优势时,一般将低估S(P-P)max值,最大误差约为30%。对于扁平轨迹误差最大,当轨迹向圆形变化时误差逐渐减小,当轨迹是圆形时误差为零。
方法3通过两个正交测量值求得Smax,当这两个测量值都是单一频率正弦形式时,方法3将是正确的。在其他场合下,方法3将过高估计S(P-P)max,误差和谐波振动分量的性质有关。
另外,标准假设水轮发电机组大轴的旋转中心是固定的,标准中GB 11348.1-1999(ISO 7919-1)中给出的旋转中心是时间积分平均。然而严格意义上,大轴的旋转中心是旋转的,大轴旋转中心会围绕固定中心(如固定参考中心,可选为机组中心或者水轮发电机组导瓦中心)旋转。
发明内容
(一)要解决的技术问题
针对现有存在的技术问题,本发明在合理假设的前提下,推导出大轴旋转中心动态变化下,大轴轴心轨迹曲线与运行摆度之间关系,通过定义摆度矢量,给出了一种准确计算水轮机运行摆度位移峰峰值的方法,为摆度信号分析提供了理论基础,对完善水电机组振动摆度在线监测系统,提高水电机组状态监测及故障诊断水平具有重工程应用价值。
(二)技术方案
为了达到上述目的,本发明采用的主要技术方案包括:
一种水轮发电机组运行摆度位移峰最大峰值准确计算方法,在假设同时满足如下条件的情况下:
条件1、水轮发电机组的大轴截面为绝对的圆形,且截面质量分布均匀,截面的质心与几何中心重合;
条件2、稳态运行下水轮发电机组的轴心轨迹为任意曲线,在不同转动周期下重合;
条件3、涡流传感器旋转不同角度,但是分布在同一个圆周上;
所述计算方法包括如下步骤:
S1、在所述大轴上安装多对涡流传感器,每一对所述涡流传感器中的一个电涡流传感器测量大轴沿其自身径向方向的偏移量,另一个测量大轴沿垂直自身径向方向的偏移量;
S2、获取旋转中心围绕固定参考中心相对转动角速度夹角(θ1)和大轴几何中心围绕旋转中心相对转动角速度夹角(θ2);
所述固定参考中心为机组中心;
S3、根据多个涡流传感器所测量的偏移量,依据下式,获得大轴摆度值在X轴、Y轴的分量:
式中:ex表示+X方向摆度值;ey表示+Y方向摆度值;e1表示固定参考中心到旋转中心的距离;e2表示旋转中心到大轴几何中心的距离;
S4、定义摆度矢量:
S5、获得摆度位移峰峰值最大值和相位:
式中:emaxx对应着最大摆度矢量下X轴摆度分量,emaxy对应着最大摆度矢量下Y轴摆度分量,∠Emax表示最大摆度矢量的角度。
进一步地,所述涡流传感器通过安装支架和底座支架安装在固定位置,且所述安装支架和底座支架的刚度相同。
进一步地,所述涡流传感器垂直于涡流传感器和大轴接触的接触面。
(三)有益效果
本发明的有益效果是:
1、本发明提供的水轮发电机组运行摆度位移峰最大峰值准确计算方法,通过理论分析提出了《GB 11348.1-1999(ISO 7919-1)旋转机械转轴径向振动的测量和评定第一部分总则》中给出的最大位移峰-峰值S(p-p)max,为工程实践中最大摆度测量提供了理论依据及技术指导。
2、本发明提供的水轮发电机组运行摆度位移峰最大峰值准确计算方法,在合理假设的前提下,推导出旋转中心动态变化下,大轴轴心轨迹曲线与运行摆度之间关系,通过定义摆度矢量,给出了一种准确计算水轮机运行摆度位移峰峰值的方法,给大轴运行摆度的趋势分析与评估提供了一种新方法,对掌握大轴与轴瓦之间油膜厚度提供了新途径。
附图说明
图1为本发明实施例中考虑大轴公转和自转下大轴摆度示意图;
图2为本发明实施例中自转几何中心圆形运动—公转中心椭圆运动轴心轨迹图;
图3为本发明实施例中自转几何中心圆形运动—公转中心椭圆运动+X+Y摆度波形图;
图4为本发明实施例中自转几何中心圆形运动—公转中心外八字形运动摆度矢量幅值波形。
具体实施方式
为了更好的解释本发明,以便于理解,下面结合附图,通过具体实施方式,对本发明作详细描述。
在进行推导之前,首先假设本实施例同时满足如下条件:
条件1、水轮发电机组的大轴截面为绝对的圆形,且截面质量分布均匀,截面的质心与几何中心重合。
工程实际中测量过程中要求摆度传感器安装位置应靠近精加工轴位置,同时大轴截面圆的圆度应满足标准要求,因此,可以在理论推导中假设大轴外形。
随着加工精度的提高,大轴的质量一般分布较为均匀。本发明暂且忽略质量分布对摆度计算的影响。
条件2、稳态运行下水轮发电机组的轴心轨迹为任意曲线,在不同转动周期下重合。
工程中实际的涡流传感器安装位置受到传感器支架及底座刚度的影响,随着机组固定部分一起振动,机组的轴心轨迹也存在着周期不重合的现象。机组大轴旋转中受到径向支撑瓦的影响,大轴的旋转中心存在着偏移、转动等现象,更进一步加剧轴心轨迹的重复性差。出于简化需求,本发明忽略传感器支架及底座刚度不均匀、旋转中心转动偏移等因素的影响,假设任意曲线的轴心轨迹在不同周期下均能重合。
条件3、涡流传感器旋转不同角度,但是分布在同一个圆周上。
实际中由于安装影响,传感器旋转过一定角度后重新安装,可能距离旋转中心的距离发生变化。实际上可以通过控制传感器安装工艺尽可能降低不同旋转角度下传感器安装偏差。故本发明忽略该因素的影响。
在如上假设条件1-3同时满足的情况下,本实施例的推导过程如下:
如图1所示,如果忽略旋转中心的运动,则该时刻下大轴在X方向的摆度是直线AB,在Y方向大轴的摆度是直线LN;如果考虑旋转中心的运动,则该时刻下大轴在X方向的摆度是直线C′B′,在Y方向大轴的摆度是直线M′N′。因此,在考虑旋转中心的运动的情况下,求出大轴的摆度即是求出直线C′B′和直线M′N′。
考虑到旋转中心O2到固定参考中心(可选为机组中心)O1的距离基本是道(10μm)别的,大轴几何中心O3到旋转中心O2也是道(10μm)级别的,故此时大轴在X轴和Y轴的摆度可以简化为:
其中,大轴几何中心O3围绕大轴几何中心O3的运动为自转/第1级相对转动,大轴几何中心O3围绕旋转中心O2的运动为公转/第2级相对转动,旋转中心O2围绕固定参考中心(可选为机组中心)O1的旋转为第3级相对转动,第3级相对转动运动形式未知。值得注意的是式(1)中的偏心距e1、e2其幅值均是时变的。
由式(1)可以看出,大轴截面处的摆度取决于第3级相对转动时旋转中心O2的轨迹和公转时大轴几何中心O3的轨迹,两者轴心轨迹综合作用下形成实际轴心轨迹,如图2所示。
现定义摆度矢量:
虽然式(2)与标准《GB 11348.1-1999(ISO 7919-1)旋转机械转轴径向振动的测量和评定第一部分总则》中给出的最大位移峰-峰值S(p-p)max相似,但本发明中式(2)的内涵已经发生了变化。GB 11348.1-1999(ISO 7919-1)中只是给出了使用其他测量量近似得到最大位移峰-峰值S(p-p)max,没有给出准确的计算方法。
结合式(1)、(2)推出式(3):
Δθ=θ1-θ2 (3)
最大摆度矢量:
式中,emaxx对应着最大摆度矢量下X轴摆度分量,emaxy对应着最大摆度矢量下Y轴摆度分量,∠Emax表示最大摆度矢量的角度。
因此,根据式(4)可以获知最大摆度幅值和相位。
当大轴自转表现为圆形时自转产生的摆度可以表示为式(5):
当大轴公转表现为外八字形时自转产生的摆度可以表示为式(6):
e1=(a1+b1cos2θ1) (6)
其中,设ɑ1=50μm,b1=200μm,e2=50μm,θ1=0~2π,θ2=0~2π,分别设置ɑ1为0、pi/6、pi/4,则对应的X轴、Y轴向摆度测量波形如图3,摆度矢量幅值波形如图4。
优选地,涡流传感器通过安装支架和底座支架安装在固定位置,且安装支架和底座支架的刚度相同。
忽略涡流传感器的安装支架刚度及底座支架的刚度的差异。工程实践中摆度的测量与涡流传感器的安装密切相关,实践中涡流传感器往往安装在专门定制的角钢支架上,而角钢支架通过特定胶水粘贴在油盆盖板或底座上。实际中由于粘贴牢固程度不同、油盆盖板存在分瓣等因素,底座的刚度及振动幅度也不同。为了简化问题,本发明暂时忽略该因素的影响。
优选地,涡流传感器垂直于涡流传感器和大轴接触的接触面。
忽略大轴静止和运动中涡流传感器探头表面与大轴表面的倾斜,假设涡流传感器探头表面与大轴接触面完全垂直。实际安装的传感器依赖于丰富的经验,可以保证涡流传感器探头表面与大轴截面垂直,本发明暂时不考虑该因素的影响。
综上,本实施例提供的水轮发电机组运行摆度位移峰最大峰值准确计算方法,通过理论分析提出了《GB 11348.1-1999(ISO 7919-1)旋转机械转轴径向振动的测量和评定第一部分总则》中给出的最大位移峰-峰值S(p-p)max,在合理假设的前提下,推导出旋转中心动态变化下,大轴轴心轨迹曲线与运行摆度之间关系,通过定义摆度矢量,给出了一种准确计算水轮机运行摆度位移峰峰值的方法,给大轴运行摆度的趋势分析与评估提供了一种新方法,对掌握大轴与轴瓦之间油膜厚度提供了新途径,为工程实践中最大摆度测量提供了理论依据及技术指导。。
以上结合具体实施例描述了本发明的技术原理,这些描述只是为了解释本发明的原理,不能以任何方式解释为对本发明保护范围的限制。基于此处解释,本领域的技术人员不需要付出创造性的劳动即可联想到本发明的其它具体实施方式,这些方式都将落入本发明的保护范围之内。
Claims (3)
1.一种水轮发电机组运行摆度位移峰最大峰值准确计算方法,其特征在于,
在假设同时满足如下条件的情况下:
条件1、水轮发电机组的大轴截面为绝对的圆形,且截面质量分布均匀,截面的质心与几何中心重合;
条件2、稳态运行下水轮发电机组的轴心轨迹为任意曲线,在不同转动周期下重合;
条件3、涡流传感器旋转不同角度,但是分布在同一个圆周上;
所述计算方法包括如下步骤:
S1、在所述大轴上安装多对涡流传感器,每一对所述涡流传感器中的一个电涡流传感器测量大轴沿其自身径向方向的偏移量,另一个测量大轴沿垂直自身径向方向的偏移量;
S2、获取旋转中心围绕固定参考中心相对转动角速度夹角(θ1)和大轴几何中心围绕旋转中心相对转动角速度夹角(θ2);
所述固定参考中心为机组中心;
S3、根据多个涡流传感器所测量的偏移量,依据下式,获得大轴摆度值在X轴、Y轴的分量:
式中:ex表示大轴在+X方向摆度值;ey表示大轴在+Y方向摆度值;e1表示固定参考中心到旋转中心的距离;e2表示旋转中心到大轴几何中心的距离;
S4、定义摆度矢量:
S5、获得摆度位移峰峰值最大值和相位:
式中:emaxx对应着最大摆度矢量下X轴摆度分量,emaxy对应着最大摆度矢量下Y轴摆度分量,∠Emax表示最大摆度矢量的角度。
2.根据权利要求1所述的水轮发电机组运行摆度位移峰最大峰值准确计算方法,其特征在于,所述涡流传感器通过安装支架和底座支架安装在固定位置且所述安装支架和底座支架的刚度相同。
3.根据权利要求1所述的水轮发电机组运行摆度位移峰最大峰值准确计算方法,其特征在于,所述涡流传感器垂直于涡流传感器和大轴接触的接触面。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201811095036.7A CN109145505B (zh) | 2018-09-19 | 2018-09-19 | 一种水轮发电机组运行摆度位移峰最大峰值准确计算方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201811095036.7A CN109145505B (zh) | 2018-09-19 | 2018-09-19 | 一种水轮发电机组运行摆度位移峰最大峰值准确计算方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN109145505A CN109145505A (zh) | 2019-01-04 |
CN109145505B true CN109145505B (zh) | 2022-12-02 |
Family
ID=64815243
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201811095036.7A Active CN109145505B (zh) | 2018-09-19 | 2018-09-19 | 一种水轮发电机组运行摆度位移峰最大峰值准确计算方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN109145505B (zh) |
Families Citing this family (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111157107B (zh) * | 2020-02-05 | 2021-11-30 | 山东联合电力产业发展有限公司 | 一种高稳定型发电机摆度测量装置 |
Family Cites Families (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP6371668B2 (ja) * | 2014-10-21 | 2018-08-08 | Nok株式会社 | 回転アンバランス測定装置 |
US10124230B2 (en) * | 2016-07-19 | 2018-11-13 | Blast Motion Inc. | Swing analysis method using a sweet spot trajectory |
CN106525337B (zh) * | 2016-11-30 | 2019-03-29 | 北京中元瑞讯科技有限公司 | 基于在线数据的水轮发电机组磁拉力不平衡故障分析方法 |
CN107894278A (zh) * | 2017-09-26 | 2018-04-10 | 昆明理工大学 | 一种立式水泵机组轴系振动参数的试验计算方法 |
-
2018
- 2018-09-19 CN CN201811095036.7A patent/CN109145505B/zh active Active
Also Published As
Publication number | Publication date |
---|---|
CN109145505A (zh) | 2019-01-04 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN105021349B (zh) | 一种获取转子的不平衡量的方法 | |
ATE499596T1 (de) | Verfahren und vorrichtung zum dynamischen messen der unwucht eines rotors | |
US8225657B2 (en) | Method and apparatus for quantitatively detecting unbalanced state and method for detecting clamping state of a workpiece | |
CN109145505B (zh) | 一种水轮发电机组运行摆度位移峰最大峰值准确计算方法 | |
CN111475903A (zh) | 基于多偏置误差同步补偿的大型高速回转装备多级零部件动力学特性逐级测调和分配方法 | |
CN110006590B (zh) | 获取转子的不平衡量和平衡机的不平衡量的方法 | |
CN108020409A (zh) | 一种主轴回转误差的四点动态测量与分离方法 | |
CN107607013A (zh) | 一种水轮发电机组旋转中心垂直度测量与调整的方法 | |
CN103240592A (zh) | 水轮发电机组盘车机组旋转中心精确确定的方法 | |
CN109406053B (zh) | 无去重转子动平衡方法 | |
CZ303892B6 (cs) | Zpusob urcení aktuální excentricity rotujícího rotoru a diagnostika excentricity rotujícího rotoru | |
CN110243541A (zh) | 风机在线动平衡计算方法 | |
RU2426976C2 (ru) | Способ и устройство для автоматической балансировки ротора | |
CN103934659B (zh) | 高精度多级碟片堆积转子的跳动控制与优化安装方法 | |
Alauze et al. | Active balancing of turbomachinery: application to large shaft lines | |
CN104361792B (zh) | 一种转动惯量测量装置 | |
CN105478245A (zh) | 基于主轴振动检测的双自由度精密离心机副轴动不平衡量辨识方法 | |
CN108981894A (zh) | 一种同时考虑水轮机大轴公转与自转下运行摆度计算方法 | |
CN109145512B (zh) | 水电机组稳态任意轴心轨迹涡流传感器安放角度分析方法 | |
JP4140380B2 (ja) | 動不釣合算定法及び動釣合試験装置 | |
CN109460572B (zh) | 稳态工况任意轴心轨迹下水电机组运行摆度计算方法 | |
CN109847952B (zh) | 一种基于驱动电流的双轴精密离心机回转台动平衡方法 | |
RU2432557C2 (ru) | Стенд комплексного определения массово-инерционных характеристик осесимметричных роторов | |
CN210142455U (zh) | 动平衡模拟实验装置 | |
CN109540388B (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 |