程志新,張港
(華北水利水電大學(xué) 機(jī)械學(xué)院,河南 鄭州 450045)
弧形閘門因無門槽、啟門力小、操作方便等優(yōu)點(diǎn)而大量用于各種水道工作[1]。在水利工程運(yùn)行過程中,弧形閘門受到水流的脈動(dòng)作用會(huì)產(chǎn)生不同程度的振動(dòng),這種閘門與水流相互耦合產(chǎn)生的閘門振動(dòng)稱為流激振動(dòng),結(jié)構(gòu)會(huì)發(fā)生變形、磨損以及腐蝕的現(xiàn)象[2],這些都是引起弧形閘門失事破壞的源頭[3]。為避免其遭受破壞,可以從兩個(gè)方面入手:一是改善閘門進(jìn)出水流流態(tài),盡量避免異常水流荷載的產(chǎn)生;二是優(yōu)化其結(jié)構(gòu),改變自身固有頻率來避免水激勵(lì)力的高能區(qū)作用。但是由于改變水流流態(tài)一定會(huì)引起其它水工結(jié)構(gòu)的改變,實(shí)施起來較為困難。因此,首要研究的是弧形閘門結(jié)構(gòu)和動(dòng)力特性,以便有效減少弧形閘門的失事破壞現(xiàn)象的發(fā)生。
由于弧形閘門是復(fù)雜的空間板梁結(jié)構(gòu),也很難測(cè)定作用其上的水流脈動(dòng)壓力,所以在對(duì)弧形閘門進(jìn)行動(dòng)力分析時(shí),大都是先計(jì)算弧形閘門的自振頻率,在與作用水流的脈動(dòng)頻率相比較,來使其自振頻率遠(yuǎn)離高能脈動(dòng)頻率區(qū),確保閘門工作環(huán)境的安全[4]。
許多實(shí)驗(yàn)資料表明現(xiàn)在已經(jīng)能準(zhǔn)確得到水流脈動(dòng)頻率,水流脈動(dòng)主頻率主要在1至20Hz之間浮動(dòng),其中有48.3%在1至10Hz之間變化,超過20Hz的很少。閘門結(jié)構(gòu)在水中的振動(dòng)屬于流體彈與性理論研究的內(nèi)容,閘門的振動(dòng)是兩者協(xié)同作用的結(jié)果,所以要用彈性結(jié)構(gòu)與流體的耦聯(lián)振動(dòng)理論來分析其振動(dòng)特性。但是閘門振動(dòng)本身就是一個(gè)很復(fù)雜的問題,所以目前普遍采用理論分析的方法,計(jì)算閘門在空氣中的自振頻率,再通過水力學(xué)試驗(yàn)測(cè)定振動(dòng)擾力,即水流總脈動(dòng)壓力的幅值和頻率,以檢驗(yàn)閘門是否發(fā)生共振[5]。
由于作用于門體的水動(dòng)力載荷難以采用函數(shù)式表示,載荷作用機(jī)理尚不完全清晰,并且弧形閘門承受的動(dòng)水壓力是一個(gè)連續(xù)波動(dòng)的過程,時(shí)均壓力上下波動(dòng)的水流作用力就是脈動(dòng)壓力。因此,在對(duì)弧形閘門自振特性進(jìn)行研究時(shí),采用有限元模態(tài)計(jì)算與原型觀測(cè)試驗(yàn)相結(jié)合的方法,即簡化邊界條件,設(shè)定理論分析的工況參數(shù),在有限元模態(tài)計(jì)算的基礎(chǔ)上,通過原型觀測(cè)試驗(yàn),對(duì)理論計(jì)算分析和設(shè)計(jì)計(jì)算的成果進(jìn)行驗(yàn)證。
本文以陜西漢江蜀河水電站泄洪閘弧形工作閘門為研究對(duì)象,孔口尺寸13.0×24.3m(寬×高),弧形閘門半徑32m,底檻高程193.5m,支鉸高程217.6m,設(shè)計(jì)水頭23.8m。整體采用的是三支臂結(jié)構(gòu)、球鉸支承,使用懸掛式液壓啟閉機(jī)啟閉,啟閉機(jī)容量為2×4000kN。本文主要對(duì)弧形閘門在動(dòng)水試驗(yàn)水位216.15m下的自振頻率和模態(tài)進(jìn)行了計(jì)算分析。
2.2.1 有限元模型
采用AUTODESK公司商業(yè)軟件Inventor 2011建立弧形閘門有限元模型,對(duì)閘門的部分結(jié)構(gòu)進(jìn)行了一定簡化,例如省略了止水部件、側(cè)輪裝置、所有螺栓螺母、欄桿等,建立的模型如圖 1所示?;⌒伍l門有限元計(jì)算選取四面體單元(Solid187 element)和六面體單元(Solid186 element)混合而成的組合有限元模型,單元的劃分基本上按弧形閘門結(jié)構(gòu)布置上的特點(diǎn)采用自然離散的方式,將面板、橫梁腹板、橫梁翼板、縱梁腹板、縱梁翼板、吊耳、支臂腹板、支臂連接桿等構(gòu)件劃分為六面體單元,支鉸、軸承、鉸軸劃分為四面體單元。
圖1 弧形閘門有限元模型
2.2.2 邊界條件
(1)水壓力載荷
水壓力直接作用于面板外壁上,水體密度取 1000kg/m3[5],面板分布水壓力根據(jù)水頭按下式計(jì)算:
靜力擋水工況閘門水頭為22.65m,荷載動(dòng)力系數(shù)取1.0。水壓力寬度范圍為面板寬度 12.9m,高度范圍為面板底緣到正常高水位 22.65m,水壓力從頂部的 0逐漸增加到底部的221970Pa。
(2)閘門自重
閘門自重方向向下,由三維建模軟件自動(dòng)計(jì)算,閘門門葉重 165724kg,支臂重 198623kg,支鉸(鉸座+鉸軸)重53000kg,軸承重 2194kg,(門葉+支臂+鉸鏈+鉸軸)為419541kg。
(3)止水摩擦力
考慮側(cè)止水橡皮的摩阻力,側(cè)橡膠止水型號(hào)為 L60(SF6674),預(yù)壓縮 4mm,橡皮頭部寬 60mm,與不銹鋼摩擦系數(shù)為0.5。
通過有限元計(jì)算得到弧形閘門在動(dòng)水試驗(yàn)水位下的模態(tài)振型和參數(shù),這里取前4階。其前4階自振特性分析結(jié)果如表1所示,振型圖取前3階圖3所示。
表1 閘門結(jié)構(gòu)自振頻率 單位:Hz
模態(tài)分析結(jié)果表明,閘門全閉擋水工況的基頻在 4.4Hz左右。由圖2可以看出,弧形閘門的第一階振型以門葉扭曲為主;第二階振型以支臂部分的彎曲為主,彎曲方向向內(nèi),呈靠攏趨勢(shì);第三階振型以支臂和門葉部分同時(shí)彎曲為主,彎曲方向?yàn)閆向。
圖2 弧形閘門振型圖
通過改進(jìn)隨機(jī)子空間法得到弧形閘門的自振頻率,經(jīng)計(jì)算整理,結(jié)果如表2所示。
表2 弧形閘門自振頻率
由表2可以看出,主梁的自振頻率要大于支臂,其中支臂的自振頻率從一階到四階增加較為穩(wěn)定,而主梁的自振頻率在第三階和第四階之間有較大增加。
對(duì)弧形閘門自振頻率的有限元模態(tài)計(jì)算結(jié)果和模態(tài)測(cè)試結(jié)果進(jìn)行整理,如表3所示。
表3 弧形閘門自振頻率模態(tài)計(jì)算結(jié)果與測(cè)試結(jié)果對(duì)比表
由表3可以看出,主梁測(cè)試所得到的自振頻率值最大,閘門經(jīng)計(jì)算得到的自振頻率值最小?;⌒伍l門前三階,也就是在低階階段的自振頻率計(jì)算值與測(cè)試值契合度較高。但是弧形閘門存在制造和裝配誤差,運(yùn)行條件無法預(yù)判等問題,實(shí)際運(yùn)行工況下的邊界條件遠(yuǎn)遠(yuǎn)要比理論計(jì)算繁瑣,使得四階(含四階)以上自振頻率測(cè)試值高于計(jì)算值,但都有一致增大趨勢(shì)[6]。
本文對(duì)弧形閘門自振特性進(jìn)行了研究,提出了一種將有限元模態(tài)計(jì)算與原型觀測(cè)試驗(yàn)相結(jié)合的方法,并以蜀河水電站弧形閘門為例展開算例分析,先通過有限元模態(tài)計(jì)算得到弧形閘門的自振頻率,然后通過原型觀測(cè)試驗(yàn)中的模態(tài)測(cè)試進(jìn)行驗(yàn)證與分析。本文所提出的方法能極大滿足水工閘門工程分析的要求,計(jì)算結(jié)果為弧形閘門振動(dòng)故障診斷及預(yù)報(bào)及其結(jié)構(gòu)動(dòng)力特性的優(yōu)化設(shè)計(jì)提供了依據(jù),可供同類型閘門振動(dòng)穩(wěn)定性研究借鑒。