某水庫潛孔弧形閘門失穩(wěn)事故分析

張懷仁，古小七，洪 偉，盛旭軍，鄭 莉，孔垂雨

(1.水利部水工金屬結(jié)構(gòu)質(zhì)量檢驗測試中心，河南 鄭州 450044；2.水利部綜合事業(yè)局，北京 100053)

2019年9月中旬，建成于20世紀70年代并于2002年進行除險加固的陜西某水庫潛孔弧形工作閘門發(fā)生失穩(wěn)事故(以下簡稱事故工作閘門)。事故經(jīng)過為：水庫進入主汛期后遭遇連續(xù)強降雨，庫水位快速上漲并超過汛限水位，為保證水庫安全，按上級防汛指令提閘泄洪，操作人員按“啟門”按鈕，閘門開度至50cm時運行正常。排沙5min后按“閉門”按鈕關(guān)閉閘門時，當(dāng)閘門閉門開度至3cm時，開度儀數(shù)據(jù)未變化，液壓站壓力表顯示管道壓力正常，此時，一聲巨響從閘門底部傳來，閘室鋼蓋板巨烈震動，開度儀無規(guī)律跳動，按“緊急停止”按鈕無變化，閘室檢修平臺檢修口有水涌出并迅速上升至液壓缸體。此次事故雖未造成人員傷亡，但作為肩負著下游一個約30萬人口城市生活供水的唯一水源地，給當(dāng)?shù)厝嗣袢罕娚a(chǎn)生活造成重大影響。

我國是世界上水庫數(shù)量最多的國家，但受客觀因素影響，病險水庫數(shù)量大，安全隱患多是當(dāng)前我國大壩安全管理中最突出的問題，金屬結(jié)構(gòu)和機電設(shè)備不安全是我國水庫安全的主要問題之一[1-4]。水工金屬結(jié)構(gòu)是水利樞紐正常有效運行的重要組成部分，其功能主要是攔截水流、控制流量、調(diào)節(jié)水位、排放漂浮物等[5-6]?；⌒伍l門相較于水工金屬結(jié)構(gòu)中其他型式的閘門具有啟閉省力、運轉(zhuǎn)可靠、受力合理、流態(tài)良好等諸多優(yōu)點，是當(dāng)前水工金屬結(jié)構(gòu)中應(yīng)用最廣泛的門型之一。水利工程建設(shè)者們通過大量的工程設(shè)計和建造經(jīng)驗避免閘門失穩(wěn)事故的發(fā)生，但國內(nèi)外依然存在一些弧形閘門失穩(wěn)案例，如陜西某渠首沖刷閘、遼寧渾河攔河閘、美國阿肯色河航運系統(tǒng)的閘門遭到破壞等[7-9]，閘門失穩(wěn)不僅會危及下游人民群眾的生命財產(chǎn)安全，往往還會造成巨大的經(jīng)濟損失，據(jù)不完全統(tǒng)計，近年來我國已有近30起弧形閘門失穩(wěn)破壞的實例[10]。因此，對該事故工作閘門進行了現(xiàn)場檢測與復(fù)核計算，分析事故原因，為其他相關(guān)工程安全運行提供借鑒。

1 事故工作閘門檢測數(shù)據(jù)

為分析事故原因，檢測人員對該事故閘門及配套的液壓啟閉機進行了巡視檢查、外觀檢測、腐蝕檢測、材質(zhì)檢測、安裝幾何位置測量等現(xiàn)場檢測。經(jīng)現(xiàn)場檢測發(fā)現(xiàn)，該事故工作閘門左右側(cè)支臂在靠近支鉸軸約0.9m處出現(xiàn)扭曲變形。檢測人員對事故工作閘門支臂材質(zhì)(牌號為Q235B，型號為20A槽鋼)進行化學(xué)成分分析和力學(xué)性能檢驗，檢驗結(jié)論為“合格”(檢測結(jié)果見表1)，對閘門安裝空間位置(未受可見性破壞部分)進行三維坐標測量，檢測結(jié)果顯示事故工作閘門左右兩支鉸軸的同軸度為16.9mm。

表1 事故工作閘門化學(xué)成分與力學(xué)性能檢測成果表

2 事故工作閘門復(fù)核計算數(shù)據(jù)

有限單元法是通過對連續(xù)結(jié)構(gòu)進行有限數(shù)目的離散化，有效地分析復(fù)雜結(jié)構(gòu)、處理復(fù)雜問題的數(shù)值計算分析方法。其原理是把連續(xù)的整體結(jié)構(gòu)劃分為有限數(shù)量的單元體，將單元富有代表性的點(一般為端點)設(shè)置為結(jié)點，從而使得相鄰單元具有連續(xù)性，構(gòu)成一個整體來代替原有的連續(xù)型結(jié)構(gòu)[11-12]。

事故工作閘門主要由面板、主橫梁、次梁、主縱梁、橫梁后翼板、上支臂、下支臂、支臂間連接梁、裙板及支鉸等金屬結(jié)構(gòu)件組成。閘門孔口尺寸為2.0m×2.4m(寬×高)，設(shè)計水頭為41.33m，校核水頭為46.34m，面板總水壓力為2810.8kN，閘門為單節(jié)弧形門葉，雙支臂，門葉頂部設(shè)置單吊點吊耳，閘門門葉及吊耳、支臂等主要結(jié)構(gòu)材料為Q235B鋼(支臂為槽鋼)，閘門主提升啟閉機為320kN單吊點液壓式啟閉機，閘門工作方式為動水啟閉。

2.1 三維實體模型

根據(jù)弧形閘門結(jié)構(gòu)設(shè)計，在SolidWorks中建立事故工作閘門三維實體模型并導(dǎo)入ANSYS進行網(wǎng)格劃分及邊界條件設(shè)定。閘門整體采用實體單元進行網(wǎng)格劃分，單元劃分按閘門結(jié)構(gòu)布置上的特點采用自動網(wǎng)格劃分，由ANSYS網(wǎng)格處理模塊面板、主橫梁腹板翼板、主縱梁腹板翼板、吊耳板、支臂、支臂裙板等規(guī)則構(gòu)件離散為八節(jié)點六面體實體單元，部分不規(guī)則板件(如支鉸座、活動支鉸及支鉸軸等)由網(wǎng)格劃分器自動離散為Solid187四節(jié)點四面體單元。實體之間接觸的部分在前處理系統(tǒng)中通過印記面(Imprint Face)的形式，使模型在劃分網(wǎng)格的時候，對共享印記面的部分劃分為接近相同的網(wǎng)格，從而使總裝配體模型中各個子模型的相鄰節(jié)點實現(xiàn)共享或在接觸面中實現(xiàn)較好的綁定，這樣可避免在有限元計算過程中，通過接觸而相鄰的兩個模型因接觸面節(jié)點不一致而出現(xiàn)局部單元剛性過大的情況，通過劃分較為細密的網(wǎng)格來保證較高的計算精度(三維模型如圖1所示)。

圖1 事故工作閘門全實體三維模型

在ANSYS環(huán)境下，閘門模型坐標系設(shè)定如下：坐標原點設(shè)置在弧形閘門面板底部中間結(jié)合處中點部位，x軸正向為垂直水流向，z軸正向指向水流下游方向，y軸正向指向重力反方向，x軸、y軸、z軸符合笛卡爾坐標系右手螺旋法則。網(wǎng)格劃分綜合考慮所校核閘門幾何尺寸以及計算機計算資源與計算精度的平衡性，設(shè)置網(wǎng)格大小為20mm。劃分完成的網(wǎng)格質(zhì)量，大部分網(wǎng)格的長寬比(Aspect Ratio)為1～10.8，雅各比比率(Jacobian Ratio)為1～6.3，翹曲系數(shù)和平行度偏差在合理范圍內(nèi)，網(wǎng)格質(zhì)量較好。

2.2 約束邊界條件

事故工作閘門正常擋水工況下的載荷和約束施加方式如下：在支鉸座底面施加固定約束(Fixed Support)用來模擬支鉸座在邊墻上的固定；弧形門葉左右兩端側(cè)水封部位施加無摩擦支撐(Frictionless Support)用來模擬門葉側(cè)導(dǎo)輪對門葉的x方向約束作用，同時避免門葉向其他方向出現(xiàn)位移；閘門吊耳處施加y方向遠程約束(Remote Displacement)以模擬閘門在閉門時的約束作用；水壓載荷作用在閘門面板正面區(qū)域處。

弧形閘門起吊工況下的約束施加方式如下：在支鉸座底面施加固定約束用來模擬支鉸座在邊墻上的固定；弧形門葉左右兩端側(cè)水封部位施加無摩擦支撐用來模擬門葉側(cè)導(dǎo)輪對門葉的x方向約束作用，同時避免門葉向其他方向出現(xiàn)位移；水壓載荷作用在閘門面板正面區(qū)域處，吊耳孔施加320kN的軸承載荷(啟閉力)，方向y方向。

有限元計算模型各工況的約束邊界條件如圖2—3所示。

圖2 正常工作載荷約束邊界條件的設(shè)置示意圖

圖3 起吊工況載荷約束邊界條件的設(shè)置示意圖

事故工作閘門于2004年投入使用，運行已約16年，時間系數(shù)取0.90；閘門為弧形工作門，重要性調(diào)整系數(shù)取0.95，所以其綜合許用應(yīng)力調(diào)整系數(shù)為0.855(閘門實際許用應(yīng)力值見表2)。

表2 閘門實際容許應(yīng)力

2.3 有限元計算結(jié)果分析

(1)整體結(jié)構(gòu)分析

在約束邊界條件環(huán)境下，用ANSYS計算閘門整體結(jié)構(gòu)有限元結(jié)果為：在正常工作工況下，閘門主結(jié)構(gòu)在最大綜合von Mises應(yīng)力為282.33MPa，超過Q235B鋼局部承壓最大容許應(yīng)力為205.2MPa，超過Q235B鋼屈服應(yīng)力為235.0MPa(主結(jié)構(gòu)應(yīng)力位移數(shù)據(jù)詳見表3)，應(yīng)力最大位置在支臂部件從支鉸處上游方向第一根斜撐同上支腿交界處(圖4)，其原因是門葉在水壓荷載下發(fā)生變形時將水壓荷載傳遞至支撐臂，支腿作為承力構(gòu)件發(fā)生變形，變形過程中部分斜撐同上支腿交界處出現(xiàn)最大擠壓應(yīng)力；且上下支腿部件的應(yīng)力超過規(guī)范要求的區(qū)域面積較大，事故工作閘門整體結(jié)構(gòu)部分區(qū)域強度不足。

表3 正常工作工況下事故工作閘門主結(jié)構(gòu)應(yīng)力位移表

圖4 主結(jié)構(gòu)最大Vonmises應(yīng)力位置(單位：MPa)

(2)上下支腿構(gòu)件分析

對上支腿及下支腿進行有限元計算，上下支腿部件的最大綜合von Mises應(yīng)力為271.32MPa，超過Q235B鋼局部承壓最大容許應(yīng)力為205.2MPa，超過Q235B鋼屈服應(yīng)力為235.0MPa，上下支腿部件絕大部分區(qū)域的綜合von Mises應(yīng)力均在182.71～240MPa，很大部分區(qū)域應(yīng)力超過Q235B鋼局部承壓最大容許應(yīng)力為205.2MPa(圖5)；上下支腿部件Z向應(yīng)力基本處于壓應(yīng)力狀態(tài)，其大部分區(qū)域在-137～-205MPa(圖6)，大部分區(qū)域明顯超出了Q235B容許應(yīng)力(抗拉抗壓應(yīng)力為136.8MPa)，上下支腿部件YZ平面剪應(yīng)力在-84～-105MPa(圖7)，大部分區(qū)域剪應(yīng)力超過了Q235B抗剪容許應(yīng)力為81.2MPa，支腿的強度不足，易導(dǎo)致閘門在正常工況下失穩(wěn)。

圖5 上下支腿von Mises應(yīng)力云圖(單位：MPa)

圖6 z軸應(yīng)力云圖(單位：MPa)

圖7 yz面剪應(yīng)力云圖(單位：MPa)

3 結(jié)論與建議

(1)對事故工作閘門在正常工作工況下的復(fù)核計算顯示，閘門主結(jié)構(gòu)在最大綜合von Mises應(yīng)力為282.33MPa，超過Q235B鋼局部承壓最大容許應(yīng)力為205.2MPa，超過Q235B鋼屈服應(yīng)力為235.0MPa，應(yīng)力最大位置在支臂部件從支鉸處上游方向第一根斜撐同上支腿交界處。計算結(jié)果與現(xiàn)場檢測成果吻合，事故工作閘門存在設(shè)計不恰當(dāng)、支臂結(jié)構(gòu)設(shè)計單薄(且采用槽鋼制造)、支臂結(jié)構(gòu)強度不足等問題，在正常擋水工況下支臂失穩(wěn)發(fā)生彎扭屈，從而導(dǎo)致整扇閘門破壞。

(2)事故發(fā)生前，事故工作閘門為高水頭小開度運行，閘門開啟瞬時與水體接觸面積大，受水壓荷載影響導(dǎo)致較大的動水激勵，也是導(dǎo)致閘門失穩(wěn)因素之一。

(3)事故工作閘門左右側(cè)支鉸同軸度安裝存在較大誤差，閘門啟閉不夠平順，閘門啟閉運行時摩阻力較大，長期運行也易產(chǎn)生失穩(wěn)。

我國大部分病險水庫除險加固工作已經(jīng)完成且水庫安全狀況改善顯著，但部分水庫水工金屬結(jié)構(gòu)設(shè)備仍存有安全隱患，加強已建水庫水工金屬結(jié)構(gòu)安全評價檢測是消除隱患的重要技術(shù)手段。