麥麥提明·依比布拉

(新疆維吾爾自治區(qū)塔里木河流域喀什管理局，新疆 喀什 844700)

水工鋼閘門是水工建筑物中的重要擋水結(jié)構(gòu)之一[1]，也是整個(gè)水利工程布設(shè)中一個(gè)至關(guān)重要的組成部分[2]，閘門的安全運(yùn)行影響著水電站效益的發(fā)揮以及下游人民群眾生命財(cái)產(chǎn)的安全[3]，在水利工程中發(fā)揮著重要作用[4]。閘門在長期運(yùn)行過程中，一些構(gòu)件會產(chǎn)生變形、銹蝕、焊縫、裂紋等缺陷[5]。出現(xiàn)變形缺陷后，不僅會使閘門的整體應(yīng)力重新分布，還可能會在某些重要部位產(chǎn)生應(yīng)力集中現(xiàn)象，從而影響閘門整體的安全運(yùn)行[6]。目前，針對閘門靜動力特性仿真分析開展了大量研究[7-8]，為工程設(shè)計(jì)及安全運(yùn)行提供寶貴經(jīng)驗(yàn)，而有關(guān)存在變形缺陷閘門有限元仿真分析的研究較少，其在運(yùn)行過程中的靜力特性分布規(guī)律有待深入研究。有鑒于此，從保障國民經(jīng)濟(jì)安全與水閘工程安全運(yùn)行的視角出發(fā)，有必要開展閘門構(gòu)件出現(xiàn)變形缺陷時(shí)應(yīng)力應(yīng)變的數(shù)值及分布影響的計(jì)算分析，以期揭示存在變形缺陷對平面閘門靜應(yīng)力特性的影響規(guī)律，計(jì)算結(jié)果可為實(shí)際工程中出現(xiàn)類似狀況提供參考。

1 計(jì)算原理及方法

水工鋼閘門正常運(yùn)行時(shí)，閘門材料均處于線彈性階段，其各構(gòu)件產(chǎn)生形變均較小，因此計(jì)算時(shí)為線彈性薄板問題，采用的基本方程[9]為：

閘門薄板廣義上的應(yīng)變矩陣{χ}如下所示：

(1)

應(yīng)力的矩陣{M}表示為：

(2)

由于薄板為各方向同性質(zhì)的均勻材質(zhì)，考慮水工鋼閘門材料的實(shí)際情況，視各板結(jié)構(gòu)均為彈性變形，式(1)、式(2)中矩陣{χ}、{M}的幾何關(guān)系可表示為：

{M}=[D]{χ}

(3)

式(3)中的[D]矩陣為水工鋼閘門材料的彈性矩陣，是由各構(gòu)件材料的彈性模量E與泊松比ν共同決定，屬于材料的固有屬性，矩陣[D]的計(jì)算公式可以表示為：

(4)

2 模型構(gòu)建

水工平面鋼閘門是由面板、主橫梁、水平次梁、縱梁以及滑輪等組件構(gòu)成，面板由主橫梁、水平次梁、縱梁以及邊梁形成的梁格支承，將面板焊接在梁格上。某鋼閘門尺寸6.5 m×4.3 m(寬×高)，板梁結(jié)構(gòu)。計(jì)算工況選取上游無水，下游水位61.5 m，閘門作用水頭13 m(底檻高程48.5 m)的設(shè)計(jì)工況。

開展閘門有限元分析時(shí)，常見主要有3類建模方式：板梁組合結(jié)構(gòu)建模、部分空間薄壁結(jié)構(gòu)建模以及完整空間薄壁結(jié)構(gòu)建模。由于完整空間薄壁結(jié)構(gòu)建模方式最接近于實(shí)際水工閘門的運(yùn)行工況，其將水工閘門各構(gòu)件用板殼單元離散，閘門仿真結(jié)果精確度更高、可靠性更強(qiáng)。有鑒于此，此次仿真計(jì)算采用的為第3種完整空間薄壁結(jié)構(gòu)建模方式。

2.1 計(jì)算坐標(biāo)

計(jì)算坐標(biāo)系的定義如下：X軸為順河向方向，Y軸為豎向，以鋼閘門主橫梁軸向?yàn)閆方向。

2.2 網(wǎng)格剖分

平面鋼閘門是一種復(fù)雜的空間薄壁結(jié)構(gòu)體系，其材質(zhì)均為不同厚度的鋼板，因而建模時(shí)需選用不同的單元來模擬各個(gè)部件。由于該平面閘門面板、主橫梁、小橫梁、縱梁以及邊梁鋼板的厚度均遠(yuǎn)小于閘門長、寬尺寸，為準(zhǔn)確得到水工平面鋼閘門各構(gòu)件的仿真力學(xué)特性，將其離散為板殼單元，而滑輪則用實(shí)體單元來模擬。

在劃分網(wǎng)格時(shí)應(yīng)該避免出現(xiàn)重合的單元節(jié)點(diǎn)，且在單元連接處需有共同的節(jié)點(diǎn)，以保證力矩傳遞的連續(xù)性與變形的協(xié)調(diào)性。該鋼閘門模型的節(jié)點(diǎn)總數(shù)為893 286個(gè)，單元總數(shù)為96 899個(gè)。

為便于統(tǒng)計(jì)計(jì)算結(jié)果，將平面鋼閘門主橫梁自上向下依次編號為1#-6#，縱梁(不含邊梁)自左至右依次編號為1#-5#，小橫梁(包含頂、底梁)自上向下依次編號為1#-9#，支承滑塊自上向下依次編號為1#-6#。平面閘門共有3節(jié)門葉，每節(jié)門葉共有上下2根主橫梁，通過對未變形情況閘門有限元應(yīng)力計(jì)算結(jié)果分析發(fā)現(xiàn)，每節(jié)門葉下主橫梁的折算應(yīng)力要大于上主橫梁。由于平面閘門載荷分布和結(jié)構(gòu)位移的對稱性，為計(jì)算方便選取每節(jié)門葉的下主橫梁，定義對稱部分的彎曲情況，即2#，4#，6#主橫梁對半部分的3個(gè)梁格腹板，標(biāo)記為9個(gè)區(qū)域。閘門劃分網(wǎng)格后整體的有限元模型見圖1。

圖1 閘門劃分網(wǎng)格后整體有限元模型

2.3 荷載與邊界條件

開展變形平面閘門靜力分析時(shí)，考慮結(jié)構(gòu)的自重和水荷載的作用。為使模型的幾何性質(zhì)不變，視水工平面鋼閘門模型底部面板的中間節(jié)點(diǎn)在閘門寬度方向位移為零，并在滑塊上施加X方向的水平約束。

2.4 材料參數(shù)

根據(jù)規(guī)定[10]，大型工程的閘門和啟閉機(jī)安全運(yùn)行30年以上時(shí)，其時(shí)間調(diào)整系數(shù)為0.9。參考規(guī)定，將容許應(yīng)力調(diào)整系數(shù)取為0.9。根據(jù)實(shí)際工程的設(shè)計(jì)圖紙資料，平面閘門面板、小橫梁、縱梁等構(gòu)件所用鋼材為Q345鋼，厚度h≤16 mm，其容許應(yīng)力[σ]=230 MPa，調(diào)整后為207 MPa，[τ]=135 MPa，調(diào)整后為121.5 MPa；閘門主橫梁、邊梁等構(gòu)件所用鋼材為Q345鋼，16 mm

3 計(jì)算結(jié)果分析

考慮閘門不同位置和不同程度變形情況下，閘門構(gòu)件的應(yīng)力和形變的分布會發(fā)生不同的變化，變化值方向規(guī)定為正增負(fù)減，對比分析與正常工況下變形對閘門的靜力方面的影響。

3.1 面 板

在閘門主橫梁5，10和20 mm彎曲變形情況下，不同變形位置閘門面板的折算應(yīng)力、變化值及變化幅度見圖2-圖4。

圖2 閘門面板折算應(yīng)力

圖3 閘門面板實(shí)際變化值

圖4 閘門面板折算應(yīng)力變化幅度

由計(jì)算結(jié)果及圖2-圖4可知：①各位置變形模型下，面板的折算應(yīng)力分布規(guī)律趨勢也基本相同。1#-9#主橫梁腹板5 mm變形模型的面板的折算應(yīng)力變化幅值相差較小(≤3%)，且最大折算應(yīng)力出現(xiàn)的位置基本一致，均出現(xiàn)在由2#縱梁、3#縱梁、4#縱梁和3#主橫梁圍成的梁格的區(qū)域，且均處于較為靠近跨中的區(qū)域。綜上所述，平面鋼閘門主橫梁梁格腹板發(fā)生小彎曲變形會增加閘門面板的折算應(yīng)力，但對應(yīng)力分布影響較小。②面板折算應(yīng)力隨著變形程度的增加而增加。當(dāng)變形位置出現(xiàn)在4#和5#區(qū)域的主橫梁腹板上時(shí)(即4#主橫梁跨中區(qū)域的梁格腹板)，引起面板折算應(yīng)力的變化稍大于其他位置變形產(chǎn)生的影響。腹板產(chǎn)生微小變形對面板的折算應(yīng)力與應(yīng)力分布產(chǎn)生影響較小，面板的折算應(yīng)力未超出容許應(yīng)力范圍。③閘門面板折算應(yīng)力分布較均勻，未出現(xiàn)明顯的應(yīng)力集中現(xiàn)象。表征由面板、主橫梁、小橫梁、邊梁、縱梁組成的梁格體系能將水工平面鋼閘門面板受到的水壓力以較均勻的方式傳遞，可大幅提高面板的承壓能力，其面板梁格體系布置形式較合理。

3.2 主橫梁

選取變形程度最大的20 mm變形情況，以2#主橫梁為例，繪制模型主橫梁最大折算應(yīng)力折線，見圖5。

圖5 20 mm變形閘門主橫梁折算應(yīng)力

由計(jì)算結(jié)果及圖5可知：①與未變形情況的閘門相比，2#主橫梁1#-3#位置梁格腹板存在彎曲變形主橫梁的折算應(yīng)力變化幅度較小(≤±5%)，表征2#主橫梁上梁格變形對主橫梁產(chǎn)生的折算應(yīng)力變化幅度影響較小，最小折算應(yīng)力出現(xiàn)的位置均為主橫梁跨中線上，最大折算應(yīng)力均出現(xiàn)在鋼閘門主橫梁后翼緣與腹板的連接區(qū)域。②變形工況下，主橫梁腹板變形位置處折算應(yīng)力增加?？拷鳈M梁跨中區(qū)域的變形對主橫梁應(yīng)力分布的影響較大，產(chǎn)生的折算應(yīng)力變化比也更大，建議實(shí)際工程中開展復(fù)核。變形發(fā)生在主橫梁上靠近邊梁的區(qū)域?qū)χ鳈M梁應(yīng)力分布影響較小。

進(jìn)一步研究2#主橫梁腹板變形對主橫梁徑向應(yīng)力、軸向應(yīng)力分布的影響。繪制主橫梁最大軸向應(yīng)力折線及變化值，見圖6。軸向力及變化值見圖7。

圖6 20 mm變形主橫梁徑向應(yīng)力及變化值

圖7 2#主橫梁變形的軸向應(yīng)力及變化值

由計(jì)算結(jié)果及圖6、圖7可知：①發(fā)現(xiàn)2#主橫梁不同變形位置的模型主橫梁徑向應(yīng)力分布均呈對稱分布，主橫梁跨中區(qū)域、主橫梁腹板和后翼緣連接處區(qū)域比較大，其中最大徑向應(yīng)力均在主橫梁跨中區(qū)域，與未變形情況相比，彎曲變形對主橫梁的徑向應(yīng)力影響較小。2#主橫梁最大軸向應(yīng)力均出現(xiàn)變形梁格腹板和后翼緣連接區(qū)域，靠近跨中區(qū)域的梁格發(fā)生變形會對主橫梁上軸向應(yīng)力產(chǎn)生較大影響。由于此次仿真計(jì)算是在理想條件下開展的，考慮到實(shí)際工程中存在的焊縫未焊透、鋼板腐蝕等一系列因素，其軸向應(yīng)力可能比計(jì)算值偏大，需采取適當(dāng)措施改善應(yīng)力分布以消除安全隱患。②軸向應(yīng)力分布基本呈對稱且受變形影響較大，腹板變形位置軸向應(yīng)力增加，最大值出現(xiàn)在腹板與后翼緣連接處?？拷鳈M梁跨中區(qū)域的變形對主橫梁軸向應(yīng)力分布的影響較大，產(chǎn)生的變化比也更大，遠(yuǎn)離主橫梁跨中區(qū)域?qū)φw變形影響較小。

2#主橫梁變形情況下剪應(yīng)力具體分布為：1#-3#變形位置的主橫梁最大剪應(yīng)力與未變形情況相比變化幅度較小(≤±2%)，1#位置變形產(chǎn)生的剪應(yīng)力變化幅度較大于其他兩個(gè)位置。與未變形情況下相比，主橫梁剪應(yīng)力分布不再呈對稱分布，1#位置變形主橫梁的剪應(yīng)力分布受影響最大，在變形區(qū)域的梁格內(nèi)出現(xiàn)剪應(yīng)力急劇增大的現(xiàn)象，跨中區(qū)域的剪應(yīng)力變化值較大，最大剪應(yīng)力出現(xiàn)在變形區(qū)域的中心位置區(qū)域。

將4#和6#主橫梁20 mm變形情況下，不同變形位置主橫梁各最大應(yīng)力值繪制應(yīng)力折線，見圖8。

圖8 4#和6#主橫梁不同位置變形各應(yīng)力值

由計(jì)算結(jié)果及圖8可知，與未變形閘門相比，主橫梁腹板存在彎曲變形會引起主橫梁的折算應(yīng)力的增加，但是增加幅度均較小。4#主橫梁出現(xiàn)變形引起的變化值較大，表征閘門中下部的主橫梁對閘門的結(jié)構(gòu)穩(wěn)定性有重要影響，靠近主橫梁跨中區(qū)域的變形對應(yīng)力的影響大于靠近邊梁區(qū)域變形產(chǎn)生的影響，主橫梁應(yīng)力變化值均未超過容許應(yīng)力值。

3.3 縱 梁

以20 mm變形情況為例，選取對主橫梁影響較大的跨中區(qū)域(以1#，4#，7#位置為例)縱梁繪制折算應(yīng)力及變化值，見圖9。

圖9 折算應(yīng)力及變化值

由計(jì)算結(jié)果及圖9可知，所有變形鋼閘門模型縱梁的最大折算應(yīng)力均出現(xiàn)在3#縱梁腹板與小橫梁連接處。其原因可能為此處水壓力較大以及與小橫梁相交接受到來自小橫梁的擠壓應(yīng)力，因此該處的折算應(yīng)力較大，最靠近閘門下部的折算應(yīng)力值較大。發(fā)現(xiàn)變形工況下應(yīng)力分布情況與未變形情況變化較小。表征主橫梁腹板發(fā)生微小變形對縱梁整體折算應(yīng)力有一定影響，對應(yīng)力分布影響較小。

為更直觀地觀察變形后縱梁的應(yīng)力分布趨勢，現(xiàn)將縱梁的徑向應(yīng)力和剪應(yīng)力進(jìn)行比較，繪制應(yīng)力變化折線及變化值，見圖10。

圖10 徑向應(yīng)力和剪應(yīng)力變化值對比圖

由計(jì)算結(jié)果及10可知：①未發(fā)生變形時(shí)閘門縱梁徑向力與未發(fā)生變形模型計(jì)算結(jié)果相比，1#-9#的20 mm變形模型的縱梁徑向應(yīng)力不同位置產(chǎn)生變化的幅值差值很小(≤15%)，且最大徑向應(yīng)力出現(xiàn)的位置變化不大，均出現(xiàn)在靠近跨中的縱梁與小橫梁交接的區(qū)域。此外，下縱梁的徑向應(yīng)力分布規(guī)律趨勢也基本相同。②主橫梁腹板發(fā)生變形會對縱梁的剪應(yīng)力變化產(chǎn)生一定影響，且1#主橫梁變形情況下縱梁產(chǎn)生的剪應(yīng)力變化稍大于3#，5#主橫梁變形產(chǎn)生的影響。③發(fā)現(xiàn)徑向應(yīng)力和剪應(yīng)力分布較均勻且分布規(guī)律基本相同，無明顯的應(yīng)力集中現(xiàn)象。表征主橫梁發(fā)生微小變形時(shí)，對縱梁的應(yīng)力分布及應(yīng)力值影響均較小。

綜上所述，大型平面鋼閘門主橫梁梁格腹板發(fā)生微小彎曲變形會增加縱梁徑向算應(yīng)力，但對應(yīng)力分布影響較小。主橫梁發(fā)生小形變后對縱梁的折算應(yīng)力值影響不大，對折算應(yīng)力、徑向應(yīng)力和剪應(yīng)力的分布規(guī)律影響不大，但對縱梁徑向應(yīng)力和剪應(yīng)力值影響較大。

4 結(jié) 論

1) 與未發(fā)生變形模型計(jì)算結(jié)果相比，主橫梁變形后的閘門面板的折算應(yīng)力變化差值很小，面板的折算應(yīng)力隨著變形程度的增加而增加。并且最大折算應(yīng)力出現(xiàn)的位置基本一樣，均出現(xiàn)在由2#縱梁、3#縱梁、4#縱梁和3#主橫梁圍成的較為靠近跨中的梁格的區(qū)域。此外，各位置變形模型下面板的折算應(yīng)力分布規(guī)律趨勢也基本相同。當(dāng)變形位置出現(xiàn)在4#，5#區(qū)域的主橫梁腹板上時(shí)(即3#主橫梁跨中區(qū)域的梁格腹板)，引起面板折算應(yīng)力的變化稍大于其他位置變形產(chǎn)生的影響。

2) 與未變形閘門相比，主橫梁腹板存在小彎曲變形后主橫梁的折算應(yīng)力變化幅值較小(≤±5%)，其中3#主橫梁受影響最大。變形對主橫梁上的應(yīng)力分布影響較大，具體表現(xiàn)為：腹板開始變形位置折算應(yīng)力增加，且靠近主橫梁跨中區(qū)域的變形對應(yīng)力分布的影響較大，產(chǎn)生的折算應(yīng)力變化比例也更大。發(fā)現(xiàn)主橫梁變形后主橫梁徑向應(yīng)力分布規(guī)律基本一致，均呈對稱分布，主橫梁跨中區(qū)域、主橫梁腹板和后翼緣連接處區(qū)域比較大。變形工況下，變形對主橫梁徑向應(yīng)力的變化幅度影響較大，靠近中橫梁中部的變形產(chǎn)生的影響遠(yuǎn)大于未變形情況。

3) 主橫梁存在小彎曲形變對縱梁的折算應(yīng)力值和應(yīng)力分布影響較小，但對徑向應(yīng)力值影響較大；對剪應(yīng)力分布規(guī)律影響較小，但對剪應(yīng)力值影響較大。