CN106777635A - 一种在风机参数及电机运转参数已知条件下的空冷岛风机桥架优化方法 - Google Patents

一种在风机参数及电机运转参数已知条件下的空冷岛风机桥架优化方法 Download PDF

Info

Publication number
CN106777635A
CN106777635A CN201611117756.XA CN201611117756A CN106777635A CN 106777635 A CN106777635 A CN 106777635A CN 201611117756 A CN201611117756 A CN 201611117756A CN 106777635 A CN106777635 A CN 106777635A
Authority
CN
China
Prior art keywords
crane span
span structure
fan
vibration
formula
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.)
Granted
Application number
CN201611117756.XA
Other languages
English (en)
Other versions
CN106777635B (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.)
Beijing University of Technology
Original Assignee
Beijing University of Technology
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 Beijing University of Technology filed Critical Beijing University of Technology
Priority to CN201611117756.XA priority Critical patent/CN106777635B/zh
Publication of CN106777635A publication Critical patent/CN106777635A/zh
Application granted granted Critical
Publication of CN106777635B publication Critical patent/CN106777635B/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
    • 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
    • G06F30/00Computer-aided design [CAD]
    • G06F30/10Geometric CAD
    • G06F30/17Mechanical parametric or variational design

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Geometry (AREA)
  • General Physics & Mathematics (AREA)
  • Evolutionary Computation (AREA)
  • General Engineering & Computer Science (AREA)
  • Computer Hardware Design (AREA)
  • Computational Mathematics (AREA)
  • Mathematical Analysis (AREA)
  • Mathematical Optimization (AREA)
  • Pure & Applied Mathematics (AREA)
  • Structures Of Non-Positive Displacement Pumps (AREA)

Abstract

一种在风机参数及电机运转参数已知条件下的空冷岛风机桥架优化方法,基于风机桥架的振动情况进行实地测试,找到风机桥架的振型及振动源;分析风机旋转对风机桥架产生的激振力,并找到对风机桥架有影响的激振力为轴向力;通过将风机桥架的一阶振振型、稳态响应与相应简支梁的一阶振型、稳态响应作对比,从而将风机桥架简化为简支梁模型;通过轴向力计算公式及稳态响应范围,建立架设费用最低的桥架设计方法。通过现场测试并应用MATLAB将作出频响曲线确定桥架振动源,通过对振源的分析找到桥架的主要激振力,然后应用ANSYS仿真分析将桥架模型简化,最后得出一种基于风机轴向力公式、桥架的振动幅值范围的桥架设计方法。

Description

一种在风机参数及电机运转参数已知条件下的空冷岛风机桥 架优化方法
技术领域
本发明涉及到一种空冷岛风机桥架结构设计的方法,属于结构性能研究领域,具体涉及通过找到空冷岛风机桥架主要振动源,并简化风机桥架结构,最后设计一种在风机参数以及电机运转参数已知的条件下,使建筑开销最小的桥架设计方法。
背景技术
火力发电厂的水冷技术需要消耗大量的水资源,导致电厂建设在选址上受到极大的限制。由于北方地区水资源匮乏,近年来,水冷技术逐渐被空冷技术所替代。空冷技术对水资源依赖性较小。该技术是采用风机强制对流的方式对管束内乏汽进行冷却。由于需要冷却的乏汽量大,所以需要的空冷风机不仅多,而且规格特别大。一般的600MW机组单台机需要56个风机单元,单个风机直径达10米左右。现有的风机安装在风机桥架上,在风机运行过程中,会引起风机桥架振动。过大的振动会导致风机、电机及减速机等的破坏。2008年怀安电厂空冷岛的风机桥架振动过大,导致风机桥架上振动经常超标,后期不得不对风机桥架进行加强,总计改造费用数百万元,且延误了电厂发电工作。
加固桥架是减小振动最有效的方法,但是如果加大桥架的支撑梁,桥架的重量必然增加,而支撑桥架的钢结构以及支撑柱必然需要加固,这就直接增大了花销。所以有必要找到一个在风机以及电机运转情况已知的条件下,使开销最小的桥架设计方法。
发明内容
本发明旨在提供一种在风机,电机运转情况均已知的条件下,使开销最小的桥架设计方法。该方法的主要特点是通过对风机桥架振动源及干扰力分析后,将桥架简化为简支梁模型。然后求出简化模型的振动稳态响应。最后根据轴向力公式推导出风机桥架的优化设计方法。
本发明是采用以下技术手段实现的:
S1、基于风机桥架的振动情况进行实地测试,找到风机桥架的振型及振动源;
S2、分析风机旋转对风机桥架产生的激振力,并找到对风机桥架有影响的激振力为轴向力;
S3、通过将风机桥架的一阶振振型、稳态响应与相应简支梁的一阶振型、稳态响应作对比,从而将风机桥架简化为简支梁模型;
S4、通过轴向力计算公式及稳态响应范围,建立架设费用最低的桥架设计方法。
本发明的特点是通过现场测试并应用MATLAB将作出频响曲线确定桥架振动源,通过对振源的分析找到桥架的主要激振力,然后应用ANSYS仿真分析将桥架模型简化,最后得出一种基于风机轴向力公式、桥架的振动幅值范围的桥架设计方法。
附图说明
图1各工况频响曲线。
图2简化模型图。
图3一阶振型。
具体实施方式
步骤(1)实地测试找到桥架振型;
为了探究桥架的振型,做两组实验分别测试水平方向与垂直方向:
第一组实验,在桥架水平方向,中间位置及两端分别放置一个速度传感器,测试风机运转100%工况下,桥架水平方向振动情况。桥架传感器布置情况及测试结果中,在水平方向上,桥架中间位置和桥架端部同步振动,振速在1mm/s附近。
第二组实验,在桥架垂直方向,桥架外侧中部及桥架中部放置两个速度传感器,测试风机运转100%工况下,桥架垂直方向振动情况。桥架传感器布置情况及测试结果中,在垂直方向上,桥架外侧中部与桥架中部同步振动,振速开始时已超10mm/s,而后衰弱。
通过这两组实验,发现桥架振动时,有破坏性的振动在垂直方向,即以垂直方向一阶振动为主。
步骤(2)通过现场测试找到振动源;
通过测试风机在不同转速下的桥架振动情况,分析频响曲线,找到振动源。图1为matlab处理的电机不同转速下桥架的频响曲线。将60%,75%,80%,85%,以及100%转速下电机的频率,风机经过减速器之后的旋转频率,以及桥架的振动峰值频率总结如表1。
表1
桥架振动峰值与风机旋转频率之间的比值均在5附近,而风机的叶片数恰好为5片。所以我们可以确定桥架的激振源为风机。激振频率为风机扇叶数与风机运行频率的乘积。步骤(3)分析引起桥架振动的激振力
风机旋转对桥架的作用力主要有以下几种:
由转子质量偏心引起的扰力;由扇叶平面外振动引起的扰力;由空气反作用力引起的扭矩;由空气反作用引起的竖直向下的轴向力。
其中转子的偏心扰力为离心力,方向始终指向风机轴径向,即水平方向;扇叶振动引起的扰力会使风机轴产生弯曲应变;空气反作用力引起的扭矩方向始终沿着风机旋转的反方向;而空气反作用力引起的轴向力始终沿着风机轴向。所以对于桥架的振动方向有影响的扰力只有由空气反作用力引起的轴向力。
步骤(4)简化桥架模型
将桥架模型简化为简支梁模型,如图2。
该模型在中点位置的挠度为
该梁的刚度为
所以该模型的固有频率
其中P为集中载荷;L为梁长度;E为弹性模量;I为截面惯性矩;M为总重量;
桥架实际固有频率由仿真分析得到为ω′,则若误差在5%以内,且振型一致,则说明该简化模型成立。
举例如下:
根据某电厂的桥架及风机参数进行验证,材料为Q235,规格HM488*300,桥架长度11.236m。桥架参数如下:
E=206GPa;m=129×11.236=1449.4kg;
Ix=7.14×10-4m4;Iy=8.12×10-5m4;ρ=7.85×103kg/m3
A=164.4×10-4m2;l=11.236m;EIx=1.47×108;EIy=1.67×107
将数据代入(式2)得:
某电厂风机单重1338kg,减速机重1080kg,电机重1098kg,风筒1900kg,风罩2860kg。所以位于桥架中间位置风机、电机、减速机等总重M=8276kg;风机桥架总重M′=1449.4×2=2898.8kg。
将数据代入式(3)得
用ANSYS模拟实际桥架与风机情况振型如图3:
得到一阶固有频率f2=4.83869Hz
实际与简化模型计算结果误差为所以将桥架简化为此简支梁模型。
步骤(5)计算模型稳态响应及轴向力,得出桥架设计方法。
假设激振力P为风机扰力的20%,如果激振力方程
P(t)=F0cosωt (4)
则运动方程为
假设xp(t)=Xcos(ωt-φ) (6)
其中X,φ分别为幅值和相位。
对式(6)求一阶,二阶导数,有
将式(6)、式(7)、式(8)代入式(5)中得
X[(K-Mω2)cos(ωt-φ)-Cωsin(ωt-φ)]=F0cosωt (9)
用三角函数和差公式
cos(ωt-φ)=cosωtcosφ+sinωtsinφ
sin(ωt-φ)=sinωtcosφ-cosωtsinφ
化简得
X[(K-Mω2)cosφ+Cωsinφ]=F0
X[(K-Mω2)sinφ-Cωcosφ]=0
解得
所以系统的稳态响应为
由式(3),将式(10)、(12)、(13),化简为
其中:X为幅值;F0为激振力;
风机旋转时进出口气流的动量改变量为产生轴向力的冲量,假设质量为dm的圆柱气团被吹出风机时速度为v,轴向力为F,即:
Fdt=vdm (17)
由于dm=ρvdA;Q=vA
所以式(17)变为:
由于气体在风机风筒中会有损耗,取矫正系数k=0.9;
最后得出
由于风机流量与风机转速成正比关系,设比例系数为α,则
其中Qm、nm为满转的流量、转速;
式(19)变为
式(21)即为轴向力公式。其中
ρ为气体密度;
k为校正系数;
Qm为风机满转流量;
nm为风机满转转速;
n为风机转速;
A为风机截面积;
激振力
由式(16),式(22)得,
令vmax≤am/s (25)
在单层钢结构中,阻尼比可取0.05
所以
将式(2)、(27)代入式(24),得
综上,即可根据此公式进行桥架的选型。

Claims (2)

1.一种在风机参数及电机运转参数已知条件下的空冷岛风机桥架优化方法,其特征在于:
S1、基于风机桥架的振动情况进行实地测试,找到风机桥架的振型及振动源;
S2、分析风机旋转对风机桥架产生的激振力,并找到对风机桥架有影响的激振力为轴向力;
S3、通过将风机桥架的一阶振振型、稳态响应与相应简支梁的一阶振型、稳态响应作对比,从而将风机桥架简化为简支梁模型;
S4、通过轴向力计算公式及稳态响应范围,建立架设费用最低的桥架设计方法。
2.根据权利要求1所述的一种在风机参数及电机运转参数已知条件下的空冷岛风机桥架优化方法,其特征在于:
步骤(1)实地测试找到桥架振型;
为了探究桥架的振型,做两组实验分别测试水平方向与垂直方向:
第一组实验,在桥架水平方向,中间位置及两端分别放置一个速度传感器,测试风机运转100%工况下,桥架水平方向振动情况;桥架传感器布置情况及测试结果中,在水平方向上,桥架中间位置和桥架端部同步振动,振速在1mm/s附近;
第二组实验,在桥架垂直方向,桥架外侧中部及桥架中部放置两个速度传感器,测试风机运转100%工况下,桥架垂直方向振动情况;桥架传感器布置情况及测试结果中,在垂直方向上,桥架外侧中部与桥架中部同步振动,振速开始时已超10mm/s,而后衰弱;
通过这两组实验,发现桥架振动时,有破坏性的振动在垂直方向,即以垂直方向一阶振动为主;
步骤(2)通过现场测试找到振动源;
通过测试风机在不同转速下的桥架振动情况,分析频响曲线,找到振动源;matlab处理的电机不同转速下桥架的频响曲线;将60%,75%,80%,85%,以及100%转速下电机的频率,风机经过减速器之后的旋转频率,以及桥架的振动峰值频率总结如表1;
表1
桥架振动峰值与风机旋转频率之间的比值均在5附近,而风机的叶片数恰好为5片;所以确定桥架的激振源为风机;激振频率为风机扇叶数与风机运行频率的乘积;
步骤(3)分析引起桥架振动的激振力
风机旋转对桥架的作用力主要有以下几种:
由转子质量偏心引起的扰力;由扇叶平面外振动引起的扰力;由空气反作用力引起的扭矩;由空气反作用引起的竖直向下的轴向力;
其中转子的偏心扰力为离心力,方向始终指向风机轴径向,即水平方向;扇叶振动引起的扰力会使风机轴产生弯曲应变;空气反作用力引起的扭矩方向始终沿着风机旋转的反方向;而空气反作用力引起的轴向力始终沿着风机轴向;所以对于桥架的振动方向有影响的扰力只有由空气反作用力引起的轴向力;
步骤(4)简化桥架模型
将桥架模型简化为简支梁模型;
该模型在中点位置的挠度为
该梁的刚度为
所以该模型的固有频率
其中P为集中载荷;L为梁长度;E为弹性模量;I为截面惯性矩;M为总重量;
桥架实际固有频率由仿真分析得到为ω′,则若误差在5%以内,且振型一致,则说明该简化模型成立;
步骤(5)计算模型稳态响应及轴向力,得出桥架设计方法;
假设激振力P为风机扰力的20%,如果激振力方程:
P(t)=F0cosωt (4)
则运动方程为
假设xp(t)=X cos(ωt-φ) (6)
其中X,φ分别为幅值和相位;
对式(6)求一阶,二阶导数,有
将式(6)、式(7)、式(8)代入式(5)中得
X[(K-Mω2)cos(ωt-φ)-Cωsin(ωt-φ)]=F0cosωt (9)
用三角函数和差公式
cos(ωt-φ)=cosωtcosφ+sinωtsinφ
sin(ωt-φ)=sinωtcosφ-cosωtsinφ
化简得
X[(K-Mω2)cosφ+Cωsinφ]=F0
X[(K-Mω2)sinφ-Cωcosφ]=0
解得
所以系统的稳态响应为
由式(3),将式(10)、(12)、(13),化简为
其中:X为幅值;F0为激振力;
风机旋转时进出口气流的动量改变量为产生轴向力的冲量,假设质量为dm的圆柱气团被吹出风机时速度为v,轴向力为F,即:
Fdt=vdm (17)
由于dm=ρvdA;Q=vA
所以式(17)变为:
由于气体在风机风筒中会有损耗,取矫正系数k=0.9;
最后得出
由于风机流量与风机转速成正比关系,设比例系数为α,则
其中Qm、nm为满转的流量、转速;
式(19)变为
式(21)即为轴向力公式;ρ为气体密度;k为校正系数;Qm为风机满转流量;nm为风机满转转速;n为风机转速;A为风机截面积;
激振力
由式(16),式(22)得,
令vmax≤am/s (25)
在单层钢结构中,阻尼比可取0.05
所以
将式(2)、(27)代入式(24),得
综上,即可根据此公式进行桥架的选型。
CN201611117756.XA 2016-12-07 2016-12-07 一种在风机参数及电机运转参数已知条件下的空冷岛风机桥架优化方法 Active CN106777635B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201611117756.XA CN106777635B (zh) 2016-12-07 2016-12-07 一种在风机参数及电机运转参数已知条件下的空冷岛风机桥架优化方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201611117756.XA CN106777635B (zh) 2016-12-07 2016-12-07 一种在风机参数及电机运转参数已知条件下的空冷岛风机桥架优化方法

Publications (2)

Publication Number Publication Date
CN106777635A true CN106777635A (zh) 2017-05-31
CN106777635B CN106777635B (zh) 2020-04-03

Family

ID=58881283

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201611117756.XA Active CN106777635B (zh) 2016-12-07 2016-12-07 一种在风机参数及电机运转参数已知条件下的空冷岛风机桥架优化方法

Country Status (1)

Country Link
CN (1) CN106777635B (zh)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110688718A (zh) * 2019-09-27 2020-01-14 北京工业大学 一种在梁式桥架参数及振动幅值条件下的空冷风机梁式桥架设计的方法

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101655097A (zh) * 2009-09-08 2010-02-24 北京龙源冷却技术有限公司 火电站空冷系统的风机桥架优化设计方法
US8677951B2 (en) * 2011-10-07 2014-03-25 Mahle Koenig Kommanditgesellschaft Gmbh & Co. Kg Opposed piston engine
CN104636543A (zh) * 2014-12-31 2015-05-20 哈尔滨工业大学 一种基于有限差分法的重型龙门铣床横梁重力变形预测方法

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101655097A (zh) * 2009-09-08 2010-02-24 北京龙源冷却技术有限公司 火电站空冷系统的风机桥架优化设计方法
US8677951B2 (en) * 2011-10-07 2014-03-25 Mahle Koenig Kommanditgesellschaft Gmbh & Co. Kg Opposed piston engine
CN104636543A (zh) * 2014-12-31 2015-05-20 哈尔滨工业大学 一种基于有限差分法的重型龙门铣床横梁重力变形预测方法

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110688718A (zh) * 2019-09-27 2020-01-14 北京工业大学 一种在梁式桥架参数及振动幅值条件下的空冷风机梁式桥架设计的方法
CN110688718B (zh) * 2019-09-27 2024-02-02 北京工业大学 一种在梁式桥架参数及振动幅值条件下的空冷风机梁式桥架设计的方法

Also Published As

Publication number Publication date
CN106777635B (zh) 2020-04-03

Similar Documents

Publication Publication Date Title
Fang et al. Exploiting the advantages of the centrifugal softening effect in rotational impact energy harvesting
Pavlenko et al. Application of artificial neural network for identification of bearing stiffness characteristics in rotor dynamics analysis
Zakaria et al. Experimental analysis of energy harvesting from self-induced flutter of a composite beam
CN105203132B (zh) 一种谐振式振动陀螺的输出频率检测方法
CN104657541B (zh) 一种风电发电机弹性支撑的选型方法
Tang et al. Analysis of frequency and mode shape of rotating-flexible disk-drum coupled structure with non-continuous connections
CN208294713U (zh) 一种杆塔宽频谐振取能装置
CN104458128A (zh) 一种基于动力学特性的涡轮增压器转子不平衡量控制方法
CN108897973A (zh) 一种弹簧-变截面盘-叶片系统的动力学建模方法
Ma et al. Vibration response analysis caused by rubbing between rotating blade and casing
Meng et al. Analytical and numerical study on centrifugal stiffening effect for large rotating wind turbine blade based on NREL 5 MW and WindPACT 1.5 MW models
CN106777635A (zh) 一种在风机参数及电机运转参数已知条件下的空冷岛风机桥架优化方法
Du et al. Synchronization of two co-rotating rotors coupled with a tensile-spring in a non-resonant system
Muminovic et al. Numerical and analytical analysis of elastic rotor natural frequency
CN103177171A (zh) 基于解析法的船舶推进轴系纵向强迫振动计算方法
CN104978494B (zh) 一种旋转非线性压电俘能结构中确定磁铁间距的方法
Li et al. An in-plane omnidirectional flutter piezoelectric wind energy harvester
CN110309615B (zh) 一种旋转叶片固有频率的预测方法
Jagannath Evaluation of critical speed of generator rotor with external load
Tanasheva et al. Analysis of aerodynamic characteristics of two parallel rotating cylinders
Castellani et al. Dynamic experimental and numerical analysis of loads for a horizontal axis micro wind turbine
CN104199998A (zh) 一种风力发电机组桨叶动态特性分析方法
CN110390139A (zh) 一种永磁调速器铜转子临界转速计算方法
Hou et al. An experimental study on a wind-induced vibration power generator for orientation-adaptive energy harvesting
Jauhari et al. Evaluation of dynamic behavior a machine tool spindle system through modal and unbalance response analysis

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