删除或更新信息,请邮件至freekaoyan#163.com(#换成@)

考虑流体影响的水轮机叶轮固有频率

本站小编 Free考研考试/2020-03-23

李鹤, 王天任, 韩萍, 闻邦椿
东北大学 机械工程与自动化学院, 辽宁 沈阳 110819
收稿日期: 2015-06-04
基金项目: 国家自然科学基金资助项目(51175071); 中央高校基本科研业务费专项资金资助项目(N120203001,N130803001); 国家重大科学仪器设备开发专项(2013YQ470765).
作者简介: 李 鹤 (1975-),男,河南方城人,东北大学教授,博士生导师;
闻邦椿(1930-),男,浙江温岭人,东北大学教授,博士生导师,中国科学院院士。

摘要: 将工作时水轮机叶轮简化为浸在流体中的旋转圆盘,建立了流体中圆盘振动微分方程,并考虑边界条件对其固有频率的影响.分别研究了不同圆盘转速、不同流体深度以及不同流体密度对固有频率的影响.研究结果表明,与真空中圆盘的固有频率相比,沉浸在流体中的圆盘的固有频率降低.当浸在流体中的圆盘发生旋转时,圆盘的模态由驻波变为行波,不同的旋转方向将对应两个不同的固有频率值.当流体深度减小时,固有频率值降低;流体密度越大,固有频率越小.
关键词:水轮机叶轮固有频率流体深度流体密度
Natural Frequency of Turbine Impeller Considering Fluid Effects
LI He, WAGN Tian-ren, HAN Ping, WEN Bang-chun
School of Mechanical Engineering & Automation, Northeastern University, Shenyang 110819,China.
Corresponding author: LI He, E-mail: hli@mail.neu.edu.cn
Abstract: The turbine impeller was simplified to a disk immersed in fluid, and the vibration differential equation of the disk was established with considering the effect of boundary conditions on natural frequency of the disk. The influence of different disk speed, fluid depth and fluid density on the natural frequency was studied. The results show that the natural frequency of the disk immersed in fluid decreased comparing with the natural frequency of the disk in vacuum. When the disk in fluid is rotated, the mode of the disk is changed from standing wave to traveling wave. Different rotation directions correspond to two different natural frequencies. As the fluid depth reduces, the value of natural frequency decreases. The higher the fluid density is, the lower the natural frequency is.
Key Words: turbineimpellernatural frequencyfluid depthfluid density
水轮机叶轮的振动对水轮机的性能和叶轮的寿命有着重要的影响.尤其是当激振频率与叶轮的固有频率相等时,将发生严重共振,进而造成叶轮工作时的不稳定或者疲劳失效等问题[1].为了避免出现共振现象[2],必须对叶轮的固有频率进行准确的估计.但是考虑到水轮机工作时叶轮是旋转的,并且由于叶轮和流体及静止叶片之间存在相互作用,影响叶轮的固有频率的因素变得非常复杂.近些年来,采用有限元方法来解决流固耦合问题较为流行[3],该方法能够比较有效地计算结构的固有频率,但是它需要更多的计算机资源,而且只能在设计后期,叶轮的结构设计完成后才能计算其固有频率.文献[4-8]等采用简化模型来计算结构的固有频率,该方法简单直观,而且可以在叶轮的初步设计阶段快速估计叶轮的固有频率.文献[4]研究了流体中圆盘的固有频率,但是它没有考虑边界条件的影响.
文献[5]研究了简支梁在流体中的振动问题,文献[6]研究了圆柱体在流体中的动力学特性,文献[7]研究了流体中悬臂板的振动,但是均与水轮机叶轮模型相差较大.文献[8]虽然研究了无限流体域中圆盘的振动问题,但是没有考虑流体旋转的影响等等.
本文以水轮机和叶轮结构为背景,对其固有频率进行估计.首先将工作时的叶轮简化为浸在流体中的圆盘模型;然后建立流体中圆盘的振动微分方程,最后分析了不同圆盘转速和边界条件对其固有频率的影响.
1 理论模型1.1 空气中圆盘的固有频率基于薄板理论和小变形,并且忽略转动惯量和剪切变形的影响[9],建立空气中圆盘的振动微分方程:
(1)
考虑到中心固定,周边自由度的边界条件,因此,环形盘的固有频率为
(2)
式中:ρDhD分别为圆盘的密度和厚度;rout为环形盘的外径;λnm为常数,可由文献[9]中查到;D为圆盘的弯曲刚度,,其中,Eμ分别为材料的弹性模量和泊松比.
1.2 流体中圆盘的固有频率如图 1所示,考虑圆盘在充满流体的圆柱形容器中振动,容器的壁是刚性的,流体是无黏性的,并且认为圆盘相对于流体以恒定的转速ω旋转.因此有如下关系:
(3)
其中,θφ分别为圆盘和流体相对于静止参考系的夹角.
图 1(Fig. 1)
图 1 叶轮简化模型Fig.1 Simplified impeller model

用平均半径r0简化方程(1),并将流体对圆盘的作用考虑为压力的影响[5].因此有
(4)
式中:D*与圆盘的材料和几何形状有关,它与D具有相同的量纲;p为对圆盘的动压力;r0=rout分别为圆盘的内径和外径.
方程(4)的解为
(5)
流体的速度势应满足拉普拉斯方程[4]
(6)
考虑流体在刚性表面的速度为零,流体在圆盘表面的速度为圆盘的振动速度,因此有
(7)
(8)
利用方程(5)和边界条件(7),(8)可以得到流体的速度势为
(9)
流体对圆盘的动压力为[10]
(10)
将方程(10)带入到方程(4)中,可以得到
(11)
方程(11)为流体中旋转圆盘的固有频率的计算方法,其中n为节径,n的符号代表行波的旋转方向.“+”代表行波的旋转方向与圆盘的旋转方向相同,“-”代表两者方向相反.ρF为流体的密度,如令ρF=0,得到
(12)
方程(12)可以用来计算空气中圆盘的固有频率,因此D*可以由方程(2)进行校验.
2 固有频率分析及讨论2.1 流体圆盘系统参数圆盘材料为不锈钢,内径din=40 mm,外径dout=400 mm,ρD=7 800 kg/m3E=2×1011Pa,hD=5 mm,μ=0.3.流体为水,流体深度H=30 mm,ρF=1 000 kg/m3.
2.2 圆盘转速对固有频率的影响利用已知参数进行仿真,在改变圆盘转速的情况下,分别计算空气中圆盘的固有频率和沉浸在流体中的圆盘固有频率,计算结果如表 1表 2所示.
表 1(Table 1)
表 1 空气中圆盘的固有频率Table 1 Natural frequency of disk in air
节径n圆盘转速/Hz
036
-2171.32171.32171.32
+2171.32171.32171.32
-3385.47385.47385.47
+3385.47385.47385.47


表 1 空气中圆盘的固有频率 Table 1 Natural frequency of disk in air

表 2(Table 2)
表 2 水中圆盘的固有频率Table 2 Natural frequency of disk in water
节径n圆盘转速/Hz
0123456
-2118.31119.35120.38121.41122.43123.43124.43
+2118.31117.26116.20115.13114.06112.97111.88
-3304.06305.19306.31307.43308.53309.64310.73
+3304.06302.92301.78300.63299.47298.30297.13


表 2 水中圆盘的固有频率 Table 2 Natural frequency of disk in water

通过表 1可知,当n=±2时,圆盘的固有频率值相同,说明在空气中圆盘的固有频率为单值,圆盘的模态为驻波.当圆盘的转速改变时,圆盘的同阶固有频率值仍然相同,说明旋转不影响空气中圆盘的固有频率.当n=±3时结论相同.通过表 2可知,圆盘在流体中旋转的情况下,n=+2和n=-2出现了两个不同的固有频率值.因此,当圆盘浸到流体中时,圆盘的模态由驻波变为行波,从而对应两个不同的固有频率值.当n取负值时,圆盘的固有频率值随着圆盘的转速增加而变大;n取正值时,圆盘的固有频率值随着圆盘的转速增加而减小.对于n=±3得出的结论仍然相同.对比表 1表 2,可知当圆盘沉浸在流体中时,圆盘的固有频率值降低.综上,旋转对于沉浸在流体中圆盘的固有频率影响是非常重要的.
2.3 流体深度对固有频率的影响流体深度影响水轮机结构中尾水管的设计参数,因此本文研究了在不同流体深度时,圆盘固有频率的变化.通过确定圆盘转速在3 Hz时,改变流体深度H来研究沉浸在流体中圆盘的固有频率特性.图 2a图 2b分别显示了流体深度H在30~100 mm的二阶、三阶固有频率的变化.其中纵坐标被正则为无量纲常数,fH代表流体深度为H时圆盘的固有频率,f0.03代表流体深度为0.03 m时圆盘的固有频率.
图 2(Fig. 2)
图 2 不同深度流体中的固有频率Fig.2 Natural frequency in different depth fluids (a)—二阶; (b)—三阶.

图 2a可知,当流体深度小于0.08 m时,对于二阶固有频率值明显的上升趋势;当流体深度大于0.08 m时,固有频率值变化缓慢,最终趋于稳定.对于图 2b,流体深度小于0.06 m时,三阶固有频率值明显变大;当流体深度大于0.06 m时,固有频率值逐渐趋于稳定.
图 2可知,随着流体深度的增加,沉浸在流体中圆盘的各阶固有频率值起初均有所增加,但是当深度达到一定值时,圆盘的各阶固有频率值逐渐趋于稳定.而且模态阶数越大时,圆盘固有频率趋于稳定的速度越快.
2.4 流体密度对固有频率的影响考虑水轮机叶轮经常工作在不同流体等复杂工作环境,因此研究不同流体的密度对叶轮固有频率的影响是非常重要的.本文研究了在圆盘转速为3 Hz时,叶轮在三种不同种类的流体中,随着深度H改变固有频率的变化.三种流体和其密度分别为:水,ρ=1 000 kg·m-3;植物油,ρ=906.4 kg·m-3;丙酮,ρ丙酮=778.0 kg·m-3.利用已知参数,编写计算机程序,仿真结果如图 3所示.其中图 3a图 3b分别显示了浸在三种不同密度流体中圆盘的二阶、三阶固有频率值随着深度H的变化关系.
图 3(Fig. 3)
图 3 不同密度流体中的固有频率Fig.3 Natural frequency in fluid with different densities (a)—二阶; (b)—三阶.

图 3a,3b可知,发现对于不同的流体密度,圆盘的各阶固有频率值相差较大.当流体深度相同时,圆盘在三种液体中的固有频率有以下关系,水中固有频率值最小,丙酮中固有频率值最大,植物油中的固有频率介于两者之间.可见,流体密度对圆盘的固有频率具有很大影响,当流体深度相同时,流体的密度越大沉浸在流体中的圆盘固有频率值越小.
3 结论1) 空气中圆盘的固有频率值不受圆盘转速的影响,此时圆盘的模态为驻波.浸在流体中的圆盘固有频率相比于空气中的固有频率降低.并且当沉浸在流体中的圆盘发生旋转时,圆盘的模态由驻波变为行波,同时对应两个不同的固有频率值.可见对于浸在流体中的圆盘,旋转将对其固有频率产生很大的影响.
2) 不同的流体深度对圆盘的固有频率也将产生一定的影响,随着流体深度的增加,沉浸在流体中圆盘的固有频率起初有所增加,但是当深度达到一定值时,圆盘的固有频率逐渐趋于稳定.
3) 流体密度也是影响流体中圆盘固有频率的一项重要因素.结论表明,在其他条件相同时,沉浸在流体中圆盘的固有频率随着流体密度的变大而变小.
参考文献
[1]Eduard E, Carme V. Failure investigation of a large pump-turbine runner[J].Engineering Failure Analysis, 2012, 2(23) : 27–34.(0)
[2]Eduard E, Carme V. Condition monitoring of pump-turbines:new challenges[J].Measurement, 2015, 67(2) : 151–163.(0)
[3]陈锋, 王春江, 周岱. 流固耦合理论与算法评述[J].空间结构, 2012, 14(8) : 55–63.
( Chen Feng, Wang Chun-jiang, Zhou Dai. Review of theory and numerical methods of fluid-structure interaction[J].Spatial Structures, 2012, 14(8) : 55–63.)(0)
[4]Kubota Y, Ohashi H. A study on the natural frequencies of hydraulic pumps[J].Computers & Fluids, 1991, 13(10) : 589–593.(0)
[5]徐博侯, 王清. 简支梁在无限大不可压缩流体和声学流体中的固有振动问题[J].振动工程学报, 1990, 3(1) : 10–17.
( Xu Bo-hou, Wang Qing. Natural vibration of simple beam in infinite incompressible and acoustic fluid[J].Journal of Vibration Engineering, 1990, 3(1) : 10–17.)(0)
[6]袁振伟, 褚福磊, 王三保. 横向振动圆柱体在流体中的动力学特性[J].工程力学, 2006, 23(8) : 49–52.
( Yuan Zhen-wei, Chu Fu-lei, Wang San-bao. Dynamic characteristics of a lateral vibrating cylinder immersed in unbounded fluids[J].Engineering Mechanics, 2006, 23(8) : 49–52.)(0)
[7]Rodriguez C G, Flores P, Pierart F G. Capability of structural-acoustical FSI numerical model to predict natural frequencies of submerged structures with nearby rigid surfaces[J].Computers & Fluids, 2012, 64(12) : 117–126.(0)
[8]Askari E, Jeong K H, Amabili M. Hydroelastic vibration of circular plates immersed in a liquid-filled container with free surface[J].Journal of Sound and Vibration, 2013, 3(32) : 3064–3085.(0)
[9]屈维德. 机械振动手册[M].北京: 机械工业出版社, 1992: 201-218.
( Qu Wei-de. Manual of mechanical vibration[M].Beijing: China Machine Press, 1992: 201-218.)(0)
[10]Harris C M. Shock and vibration handbook[M].New York: McGraw-Hill, 1995: 804.(0)

相关话题/水轮机 叶轮

闂備浇宕甸崰鎰洪幋锔藉€块柨鏃傛櫕閻濊埖鎱ㄥ璇蹭壕闂佺娅曢幑鍥х暦濮椻偓椤㈡棃宕卞Δ鍐ㄥ箰闂傚倷绀侀崥瀣熆濮椻偓瀹曨垶顢曢姀銏℃闂侀潧顧€婵″洭鍩㈤弮鍫熺厵闁硅鍔曢惃娲煟閻旈绉洪柟顔肩秺瀹曟儼顦查柍閿嬫⒐缁绘繈鍩€椤掑嫬鐒洪柛鎰硶閻鈹戦埥鍡楃仯缁绢厼鐖煎畷鎴﹀箻鐎靛摜鎳濋梺閫炲苯澧撮柕鍡楀€块幊鐐哄Ψ閿曗偓閸斿懘姊洪崷顓炰壕濠⒀冮叄瀹曟垿骞樼€涙ê顎撶紓浣圭☉椤戝棛绱炴繝鍐х箚妞ゆ劑鍊曢。濂告煕閳轰緤鍔熸俊鍙夊姍瀹曞綊顢曢敐鍕攭婵犵數鍋涘Λ娆撳箰缁嬪簱鏋嶉柨鐔哄У閸嬶絿鎲稿鍥e亾濮樼厧骞橀柤娲憾椤㈡瑩鎮欓浣稿姎闂備礁鎲℃笟妤呭储閻撳寒鍤曢柛褎顨嗛悡銉︾箾閹寸儐娈樻い锔奸檮缁绘繈鍩€椤掑嫷鏁傞柛鈩冪懄閿涘繘姊洪崜鑼帥闁哥姵鐗滈弫顕€骞庨懞銉ヤ化闂佸綊鍋婇崜娑氱矓閻戞ü绻嗛柟缁樺笒椤忣偊鏌熸搴⌒i悗鐢靛帶閳诲酣骞嬮悪鍛簥
2婵犵數鍋為崹鍫曞箰缁嬫5娲Ω閳哄绋忛梺鍦劋椤ㄥ棝宕甸埀顒勬⒑閸涘﹤濮﹀ù婊呭仱瀹曟椽鏁撻悩宕囧幘闂佸搫顦悘婵嬪汲濞嗘垟鍋撶憴鍕8闁稿海鏁婚獮鍐╃鐎n偒妫冨┑鐐村灦绾板秹骞夐妶鍡欑闁瑰鍋為ˉ鏍偣娓氬﹦鎮兼俊鍙夊姇閳诲酣骞囬崜浣虹厬婵犵妲呴崹闈涒枍閿濆拋娼╅柕濞炬櫆閻撳繘鏌涢埄鍐╃缂佷胶澧楅妵鍕敃閵忊晛鍓遍梺鎸庣⊕閸旀瑩鐛€n亖鏋庨悘鐐村灊婢规洘绻涙潏鍓ф偧閺嬵亞鈧鍠栭悘婵嬪煘閹达箑纾兼繛鎴i哺閻g晫绱撴担椋庤窗闁稿锕ら锝夘敆閳ь剟鍩㈡惔鈾€鏋庨柟瀵稿濡叉挳姊绘担鍛婂暈闁荤喆鍎甸弫鍐Χ閸滀礁鏅犳繝鐢靛У濮樸劌鐣垫笟鈧幃妤呮濞戞瑥鏆堥梺绯曟櫅閸婃瓕鐏嬮梺缁橆殔閻楃偤宕楅鍌滅<閺夊牄鍔庣粻鎾绘煃閽樺妯€闁诡喕绮欏畷褰掝敃閵堝繒鐭�547闂傚倷绀佸﹢閬嶃€傛禒瀣;闁瑰墽绮埛鎺楁煕閺囨娅呴柣蹇d邯閺岋絽螖鐎n偄顏�4婵犵數鍋為崹鍫曞箰缁嬫5娲晲閸モ晝顦梺鐟邦嚟閸嬬喖鍩㈤弮鈧妵鍕疀閹惧顦遍梺绋款儐閹瑰洤鐣峰鈧、鏃堝幢濡ゅ啫骞愰梻鍌欑閸氬顭垮鈧畷顖烆敃閵忋垺娈鹃梺闈涱檧婵″洭鍩㈤弮鍫熺厵闁硅鍔曢惃娲煟閻旈攱璐$紒杈ㄥ笚瀵板嫮鈧綆浜炴禒鐓庘攽閻愯泛绱﹂柛妤勬珪娣囧﹦鈧稒蓱婵挳鎮峰▎娆戝埌濞存粓绠栭弻娑㈠箛閸忓摜鏁栭梺绯曟櫔缁绘繈寮诲☉婊呯杸闁挎繂鎳庨~鈺呮⒑閸涘﹦鎳曠紒杈ㄦ礋楠炲繘鎮╃拠鑼啋闂佸憡渚楅崹鎶剿囬埡鍛拺閻炴稈鈧厖澹曢梻浣告贡鏋紒銊ㄥ亹缁厽寰勯幇顓犲幘闂佸搫瀚换鎺旇姳閹稿簺浜滈柡鍥悘鑼偓娈垮櫘閸嬪﹤顕g捄琛℃瀻闁诡垎鍏俱倝姊绘担鐟邦嚋缂佸鍨胯棟妞ゆ挶鍨归悞鍨亜閹烘埊鏀婚悗姘炬嫹40缂傚倸鍊风粈渚€藝椤栨粎鐭撶€规洖娲ㄧ粻鏃堟煙閹屽殶闁崇粯姊归幈銊ヮ潨閳ь剛娑甸幖浣告瀬閻庯綆浜堕悢鍡涙偣閸ワ絺鍋撻崗鍛棜缂傚倷娴囨ご鎼佹偡閳哄懎钃熺€光偓閸曨偆顓洪梺鎸庣箓濞茬娀宕戦幘缁樻優閻熸瑥瀚崢褰掓⒑閸濆嫭宸濋柛搴㈠灩濡叉劕饪伴崼鐔哄帾闂佺硶鍓濆ú鏍姳缁寬闂傚倷绶氬ḿ褍螞閺冨倻鐭嗗ù锝堫嚉閻熸嫈鏃堝川椤旇瀵栭梻浣告啞娓氭宕㈤挊澶嗘瀺鐎光偓閸曨剛鍘遍梺鍝勫暞閹稿墽澹曢幐搴涗簻闁靛骏绱曢。鑼磼閺冨倸鏋涙い銏℃礋閹晠顢曢~顓烆棜濠电偠鎻紞鈧い顐㈩樀閹繝鍩€椤掑嫭鈷掗柛灞捐壘閳ь兛绮欓、娆愮節閸曨剦娼熼梺鍓插亝濞叉牜绮荤紒妯镐簻闁规崘娉涢宀勬煛娴e壊鍎忔い顓℃硶閹风娀鍨鹃崗鍛寜闂備線鈧偛鑻崢鎼佹煠閸愯尙鍩e┑锛勬暬閹瑩寮堕幋顓炴婵犳鍠楅敃鈺呭礂濞戞碍顫曢柨婵嗘偪瑜版帗鏅查柛鈩冪懅閻撳倸顪冮妶鍐ㄧ仾缁炬澘绉剁划鈺呮偄绾拌鲸鏅濋梺闈涚墕濞诧箓骞嗛敐澶嬧拺闁圭ǹ娴烽埥澶嬨亜閿旇寮柛鈹惧亾濡炪倖甯掗敃銉р偓姘炬嫹28缂傚倸鍊风欢锟犲磻婢舵劦鏁嬬憸鏃堢嵁閸愵喖绠婚柛鎾茬窔閳瑰繑绻濋姀锝嗙【妤楊亝鎸冲畷婵嬪川鐎涙ḿ鍘遍梺鍝勫暊閸嬫挾绱掔€n偅宕岀€殿噮鍋婇幃浠嬪垂椤愩垹骞戦梻浣告惈濞诧箑顫濋妸鈺傚仭闁宠桨璁插Σ鍫ユ煙閻愵剚缍戦柟鍏兼倐閺屽秷顧侀柛鎾卞姂楠炲繘鏁撻敓锟�1130缂傚倸鍊风粈渚€藝椤栨粎鐭撻柛鎾茬閸ㄦ繈鏌ㄩ悢鍝勑㈢紒鈧崒鐐寸厱婵炴垶锕弨缁樼箾閸繄鐒搁柡灞稿墲瀵板嫮鈧綆浜滈~鍥⒑闁偛鑻晶顕€鏌涙繝鍕殌閾荤偤鏌曢崼婵囧濠殿垰銈搁弻锝夊箣閻忔椿浜濈粋宥夊醇閺囩喓鍘搁梺鍛婃礋濞佳囨倶閿濆棛绡€闁逞屽墴閸╁嫰宕樼捄銊х厬婵犵妲呴崹闈涚暦闂堟侗鍟呮繝闈涚墢绾捐棄霉閿濆毥褰掑汲閺冨牊鐓冪憸婊堝礂濞戞碍顐芥慨妯跨堪閳ь剙鍊块幊鐐哄Ψ閿曗偓閸斿懘姊虹憴鍕姢闁哥喎娼¢幃锟犲即閻愬秵鐩幃褔宕奸鈧崣鏇犵磽娴d粙鍝虹紒璇插暙宀h法鎷犻崣鍌滃枔缁辨帒螣閻戝洨鐣辨い顓℃硶閹瑰嫰鎼归崜鎰剁磿缁辨帡鍨鹃崘宸純閻庢鍣崑鍕煝鎼淬劌绠涙い鎾跺Х閳诲矂姊绘担鐟邦嚋缂佸鍨剁缓浠嬪籍閸屾粌宕ラ柣搴㈢⊕宀e潡鍩㈤弮鍫熺厽婵﹩鍓氶~婕冩繝鐢靛仦濞兼瑩宕ョ€n喗鍤屽Δ锝呭暙缁犵喖鏌熼幆鐗堫棄闂佽¥鍊濋弻娑㈠Ψ閹存繂鏆㈡い鏂款槸閳规垿鎮╁ù瀣級闂侀潧鐗嗛幊搴ㄥ几濞嗘挻鈷戦悹鍥b偓铏亖闂佸憡鏌ㄦ鎼佸煝閹捐绠i柨鏃傜摂濮婂潡姊洪柅鐐茶嫰婢ь噣寮崼銉︾厱閻忕偛甯哄璺虹;闁瑰墽绮弲鎼佹煥閻曞倹瀚�