熊 威，劉 智，王 姣，唐少龍，萬(wàn)小強(qiáng)

（江西省水利科學(xué)研究院，江西省水工安全工程技術(shù)研究中心，江西 南昌 330029）

0 引 言

理論計(jì)算是水工結(jié)構(gòu)設(shè)計(jì)階段獲取結(jié)構(gòu)應(yīng)力、變形和安全度的主要途徑。物理模型試驗(yàn)技術(shù)是驗(yàn)證新材料和新設(shè)計(jì)方法可靠性的實(shí)踐基礎(chǔ)[1]。與前兩種方法相比，數(shù)值模擬技術(shù)不僅能夠仿真不同荷載作用下的應(yīng)力變形，同時(shí)還能驗(yàn)證新理論方法的可行性，在研究結(jié)構(gòu)應(yīng)力變形方面更為全面與靈活[2，3]。近幾十年來(lái)，隨著人工智能和計(jì)算機(jī)技術(shù)的快速發(fā)展，各種數(shù)值模擬方法逐漸被引入到結(jié)構(gòu)應(yīng)力變形的研究中[4]。周偉等[5]采用三維彈塑性有限元方法，對(duì)溪洛渡高拱壩的壩體壩基系統(tǒng)進(jìn)行了漸進(jìn)破壞過(guò)程計(jì)算分析，取得了預(yù)期成果。陳欣等[6]利用三維非線性有限元程序?qū)﹀\屏拱壩進(jìn)行了破壞仿真分析，并探討了結(jié)構(gòu)失穩(wěn)的判據(jù)。李火坤等[7]基于FLOW-3D軟件建立了堤防決口水流模擬的三維數(shù)值模型，得到封堵過(guò)程中決口附近水位場(chǎng)和流速場(chǎng)分布規(guī)律，取得了較好的效果。有限元仿真方法作為數(shù)值模擬方法中的一種，已被廣泛應(yīng)用于多個(gè)領(lǐng)域，其對(duì)于結(jié)構(gòu)應(yīng)力變形狀態(tài)研究也已經(jīng)相當(dāng)成熟[8]，研究水平已經(jīng)能夠一定程度地對(duì)結(jié)構(gòu)的破壞過(guò)程進(jìn)行仿真，與試驗(yàn)及實(shí)際工程現(xiàn)象進(jìn)行定量比較，可更直觀的反映結(jié)構(gòu)的破壞過(guò)程。

本文以某閘壩滑移破壞實(shí)例為研究對(duì)象，結(jié)合現(xiàn)場(chǎng)查勘、設(shè)計(jì)和施工資料建立有限元仿真模型，通過(guò)折減破壞部位接觸面上的抗剪斷參數(shù)，仿真結(jié)構(gòu)破壞過(guò)程，并探討其破壞原因。

1 工程概況

某閘壩工程是一座以改善城市水環(huán)境、增加旅游景觀為主，兼有河道整治、改善生態(tài)環(huán)境和發(fā)電等綜合效益的水利樞紐工程。工程主要由溢流壩段（含高堰泄水閘壩、低堰泄水閘壩）和廠房壩段組成，全長(zhǎng)300.58m。其中，高堰泄水閘壩長(zhǎng)226.50m，由22孔泄水閘和3個(gè)閘墩組成，單孔凈寬10.00m。閘壩為混凝土結(jié)構(gòu)，采用鋼質(zhì)翻板閘門(mén)擋水，溢流堰為寬頂堰，堰頂高程179.50m（黃海高程，下同），采用底流消能，消力池池長(zhǎng)10.00m，底板厚0.80m，在消力池末端設(shè)置消力坎，坎厚1.00m，高堰泄水閘壩結(jié)構(gòu)剖面如圖1所示。

本次滑移破壞發(fā)生在高堰泄水閘壩段（樁號(hào)0+074.08～樁號(hào)0+154.08），工程仍有低堰泄水閘壩段、發(fā)電廠房壩段以及部分高堰泄水閘壩段處于正常運(yùn)行狀態(tài)。閘壩結(jié)構(gòu)滑移破壞現(xiàn)場(chǎng)情況如照片1～2所示。

圖1 高堰泄水閘壩結(jié)構(gòu)剖面圖

照片1 閘壩滑移現(xiàn)場(chǎng)情況1

照片2 閘壩滑移現(xiàn)場(chǎng)情況2

2 仿真分析

2.1 仿真模型

為了探索結(jié)構(gòu)破壞的原因，本節(jié)通過(guò)建立仿真分析模型，模擬工程實(shí)際工作性態(tài)及其滑移過(guò)程。高堰泄水閘結(jié)構(gòu)在基底及上下游方向延伸約1.5倍閘墩高度形成地基，模型共有節(jié)點(diǎn)數(shù)164 237個(gè)，單元數(shù)149 872個(gè)，模型如圖2所示。本次結(jié)構(gòu)破壞過(guò)程中鋼筋未對(duì)閘墩位移起約束作用，因此仿真模型中未設(shè)立鋼筋單元。本次仿真計(jì)算設(shè)置的主要潛在滑動(dòng)面及開(kāi)裂面包括：堰體下部C10埋石混凝土與上部C20混凝土堰面之間的結(jié)合部位、二期C30混凝土與堰面C20混凝土之間結(jié)合部位、二期C30混凝土與下部C10埋石混凝土之間的結(jié)合部位、上部1.00m厚C20混凝土堰面兩次澆筑之間的結(jié)合部位。接觸面模型如圖3所示。

圖2 高堰泄水閘壩結(jié)構(gòu)有限元仿真模型

圖3 結(jié)構(gòu)潛在滑動(dòng)面及開(kāi)裂面接觸單元模型示意圖

2.2 材料力學(xué)參數(shù)

由于缺少材料力學(xué)性能參數(shù)的相關(guān)資料，本次仿真計(jì)算材料參數(shù)按照《水工混凝土結(jié)構(gòu)設(shè)計(jì)規(guī)范》（SL191-2008）混凝土標(biāo)號(hào)建議值進(jìn)行取值，具體參數(shù)如表1所示。

本次閘壩破壞形式不同于建基面滑動(dòng)或地基深層滑動(dòng)，潛在滑動(dòng)面出現(xiàn)在混凝土分層澆筑的結(jié)合部位。由于無(wú)法對(duì)現(xiàn)場(chǎng)結(jié)構(gòu)的層間抗剪斷參數(shù)開(kāi)展原位試驗(yàn)，本次仿真結(jié)構(gòu)潛在滑動(dòng)面的層間抗剪斷參數(shù)參照《碾壓混凝土壩設(shè)計(jì)規(guī)范》（SL 314-2018）進(jìn)行估值，我國(guó)部分碾壓混凝土工程層間面抗剪斷參數(shù)如表2所示，經(jīng)綜合比較分析后，工程潛在滑動(dòng)面的層間抗剪斷參數(shù)基準(zhǔn)值采用摩擦系數(shù)和凝聚力。

2.3 邊界條件及荷載組合

本次仿真計(jì)算采用基于不連續(xù)介質(zhì)力學(xué)的混合算法，在模型水平方向施加法向約束，地基底部施加三向約束，潛在張開(kāi)、滑移塊體部位不施加約束。計(jì)算荷載包括閘門(mén)承受的靜水荷載及結(jié)構(gòu)自重荷載，破壞時(shí)壩前水位為183.23m，閘門(mén)擋水水頭3.73m，通過(guò)工程潛在滑動(dòng)面的層間抗剪斷參數(shù)變化將本次仿真分析分為3個(gè)工況，其荷載組合詳見(jiàn)表3。

表1 材料力學(xué)性能參數(shù)

表2 國(guó)內(nèi)部分工程碾壓層間面抗剪斷參數(shù)表（現(xiàn)場(chǎng)原位試驗(yàn)）

表3 仿真分析各工況荷載組合

計(jì)算分為兩個(gè)運(yùn)行塊進(jìn)行，第一塊僅地基參與計(jì)算，計(jì)算的結(jié)果位移清零，應(yīng)力不清零，作為下一塊計(jì)算時(shí)的初始狀態(tài)，即扣除地基沉降，但保持初始地應(yīng)力；第二塊地基與結(jié)構(gòu)共同參與運(yùn)算，分析結(jié)構(gòu)的工作性態(tài)。

2.4 仿真結(jié)果分析

通過(guò)改變層間抗剪斷參數(shù)，總結(jié)不同參數(shù)對(duì)工程破壞形態(tài)以及位移量值的影響，初步判斷實(shí)際破壞發(fā)生時(shí)結(jié)構(gòu)潛在滑動(dòng)面部位的力學(xué)性能指標(biāo)。在充分考慮結(jié)構(gòu)間相互制約作用的同時(shí)簡(jiǎn)化仿真計(jì)算量，提高計(jì)算效率。將溢流堰面上游部位按模型分區(qū)將節(jié)點(diǎn)打斷，閘墩之間的堰面潛在滑動(dòng)面抗剪斷參數(shù)采用一較大的固定值，不隨其它接觸面抗剪斷參數(shù)進(jìn)行折減。

（1）工況1仿真結(jié)果分析

工況1時(shí)單段閘壩結(jié)構(gòu)變形示意圖見(jiàn)圖4～5，結(jié)構(gòu)各向位移量值均處于有限元范疇的微小變形階段，其中順?biāo)飨蛭灰谱畲笾禐?.55mm，垂直水流向位移最大值為8.14mm，均出現(xiàn)在閘門(mén)部位。靜水壓力作用下，閘門(mén)、閘墩及下游溢流堰面整體向下游方向移動(dòng)，其中下游溢流堰面在閘墩推力及接觸面上約束共同作用下，尾部輕微變形，此時(shí)結(jié)構(gòu)整體處于正常穩(wěn)定運(yùn)行。因此，當(dāng)層間抗剪斷參數(shù)采用基準(zhǔn)值時(shí)，潛在滑動(dòng)面及開(kāi)裂面均保持完好。

圖4 工況1時(shí)單段閘壩結(jié)構(gòu)變形示意圖

圖5 工況1時(shí)整體閘壩結(jié)構(gòu)變形示意圖

（2）工況2仿真結(jié)果分析

工況2時(shí)閘壩結(jié)構(gòu)變形示意圖見(jiàn)圖6～7，結(jié)構(gòu)各向位移量值已脫離有限元微小變形范疇，其中順?biāo)飨蛭灰谱畲笾禐?15.48mm，豎直水流位移最大值為59.03mm，均出現(xiàn)在閘門(mén)部位。靜水壓力作用下，閘門(mén)整體揚(yáng)起并向下游方向移動(dòng)，推動(dòng)閘墩與兩側(cè)溢流堰面和底部一期混凝土之間發(fā)生滑動(dòng)，下游溢流堰面在閘墩位移作用下在尾部區(qū)域開(kāi)裂，并形成與上游溢流堰面之間的貫穿裂縫，導(dǎo)致脫開(kāi)形成一塊準(zhǔn)脫離體，其破壞形態(tài)與現(xiàn)場(chǎng)實(shí)際觀測(cè)形態(tài)相似，但位移量較現(xiàn)場(chǎng)觀測(cè)數(shù)值低。因此，當(dāng)層間抗剪斷參數(shù)進(jìn)行0.5倍折減時(shí)，潛在滑動(dòng)面及開(kāi)裂面均出現(xiàn)一定程度的脫開(kāi)現(xiàn)象。

（3）工況3仿真結(jié)果分析

圖6 工況2時(shí)單段閘壩結(jié)構(gòu)變形示意圖

圖7 工況2時(shí)整體閘壩結(jié)構(gòu)變形示意圖

工況3時(shí)閘壩結(jié)構(gòu)變形示意圖見(jiàn)圖8～9，各向位移量值相比工況2進(jìn)一步增大，其中順?biāo)飨蛭灰谱畲笾禐?07.85mm，垂直水流向位移最大值為146.60mm，均出現(xiàn)在閘門(mén)部位。靜水壓力作用下閘門(mén)整體揚(yáng)起并推動(dòng)閘墩及下游溢流堰面向下游方向整體移動(dòng)，此時(shí)閘門(mén)、閘墩與下游溢流堰面形成了整體的脫離塊體。此時(shí)可初步判斷結(jié)構(gòu)已發(fā)生嚴(yán)重破壞，其破壞形態(tài)、位移量均與現(xiàn)場(chǎng)實(shí)際觀測(cè)形態(tài)相似。因此，當(dāng)層間抗剪斷參數(shù)進(jìn)行0.2倍折減時(shí)，潛在滑動(dòng)面及開(kāi)裂面均完全脫開(kāi)。

各工況下結(jié)構(gòu)變形最大位移值見(jiàn)表4。

圖8 工況3時(shí)單段閘壩結(jié)構(gòu)變形示意圖

圖9 工況3時(shí)整體閘壩結(jié)構(gòu)變形示意圖

表4 各工況下結(jié)構(gòu)變形最大位移值 mm

3 結(jié) 論

本文根據(jù)某閘壩汛期滑移破壞情況，結(jié)合現(xiàn)場(chǎng)查勘、設(shè)計(jì)和施工資料建立有限元仿真模型，將仿真分析結(jié)果與工程現(xiàn)場(chǎng)進(jìn)行對(duì)比，結(jié)果基本吻合，驗(yàn)證了有限元仿真分析方法的可靠性。

由仿真結(jié)果可知，潛在滑動(dòng)面及開(kāi)裂面的層間抗剪斷參數(shù)約在，f'=0.2～0.5，c'=0.2～0.5MPa，遠(yuǎn)低于正常值。因此，層間結(jié)合部位抗剪斷性能降低是導(dǎo)致本次閘壩出現(xiàn)滑移破壞的主要原因。