程 立，劉耀儒，呂慶超，楊 強(qiáng)

(1.水電水利規(guī)劃設(shè)計(jì)總院，北京 100120；2.清華大學(xué)水沙科學(xué)與水利水電工程國(guó)家重點(diǎn)試驗(yàn)室，北京 100084；3.中國(guó)電建集團(tuán)西北勘測(cè)設(shè)計(jì)研究院有限公司，陜西 西安 710065)

0 引 言

特高拱壩是水電開發(fā)中重要壩型。近年來，隨著我國(guó)錦屏一級(jí)、溪洛渡等7座壩高超過200 m級(jí)的特高拱壩先后建成并蓄水運(yùn)行，為特高拱壩建設(shè)積累了豐富資料和寶貴經(jīng)驗(yàn)。而在工程實(shí)踐中，卸荷松弛等勘測(cè)設(shè)計(jì)期未充分重視且常規(guī)方法難以定量分析的變形破壞問題日益突出。自然條件下的巖體結(jié)構(gòu)經(jīng)過漫長(zhǎng)的地質(zhì)演化過程而處于天然的平衡狀態(tài)，F(xiàn)airhurst等稱之為預(yù)先存在的平衡(Pre-existing equilibrium)[1]。開挖、蓄水等劇烈的工程擾動(dòng)破壞了這種平衡，擾動(dòng)巖體出現(xiàn)與時(shí)間有關(guān)的應(yīng)力和變形重分布過程，該過程可稱為非平衡演化，直接影響特高拱壩等巖體工程結(jié)構(gòu)安全建設(shè)及長(zhǎng)期穩(wěn)定運(yùn)行[2]。

近年來，在巖體工程實(shí)踐中，非平衡演化及其對(duì)工程整體穩(wěn)定性影響研究的緊迫性日益突出。如，貫穿錦屏一級(jí)水電站“開挖-澆筑-蓄水”全周期的左岸高邊坡時(shí)效變形[3]；小灣、二灘、白鶴灘等水電站特高拱壩建基面開挖產(chǎn)生嚴(yán)重的卸荷松弛現(xiàn)象[4]；金川礦區(qū)水平巷道變形破壞[2]以及深挖隧洞洞壁大變形、巖爆等，均可認(rèn)為是工程實(shí)踐中與時(shí)間有關(guān)的變形破壞現(xiàn)象。洞室開挖的新奧法施工原理(NATM)是巖體結(jié)構(gòu)非平衡演化的很好例證，陳宗基院士指出：“新奧法原理是一般流變結(jié)構(gòu)演化規(guī)律甚至地球動(dòng)力學(xué)演化規(guī)律在隧道工程背景下的縮影”[2]。

目前，特高拱壩工程實(shí)踐中穩(wěn)定性數(shù)值分析主要基于有限元法，但其核心思想來源于極限分析，即超載法和強(qiáng)度折減法。主要的判據(jù)可分為2類：①依賴于算法的穩(wěn)定性判據(jù)，如殘余力大小[5]等；②經(jīng)驗(yàn)穩(wěn)定性判據(jù)，如塑性區(qū)貫通[6]等。但這些穩(wěn)定性判據(jù)缺乏嚴(yán)格的理論基礎(chǔ)。根本原因在于極限分析法采用時(shí)間無(wú)關(guān)的觀點(diǎn)研究結(jié)構(gòu)穩(wěn)定性，認(rèn)為結(jié)構(gòu)在給定荷載作用下只有穩(wěn)定或失穩(wěn)2種可能的最終狀態(tài)，忽略了結(jié)構(gòu)達(dá)到最終狀態(tài)的時(shí)效演化過程，而嚴(yán)格的穩(wěn)定性理論應(yīng)基于動(dòng)態(tài)演化觀點(diǎn)。此外，動(dòng)態(tài)演化觀點(diǎn)對(duì)結(jié)構(gòu)長(zhǎng)期穩(wěn)定性具有重要意義，對(duì)此極限分析方法很難適用。

針對(duì)有限元法缺乏嚴(yán)格的結(jié)構(gòu)穩(wěn)定性判據(jù)這一問題，楊強(qiáng)等[7-9]提出變形加固理論(Deformation Reinforcement Theory)。變形加固理論的發(fā)展經(jīng)歷了2個(gè)階段：早期是以彈塑性理論框架和結(jié)構(gòu)失穩(wěn)定義[7-9]為基礎(chǔ)；后期基于穩(wěn)定演化為基礎(chǔ)的時(shí)效變形加固理論[10]，最終形成基于Duvaut-Lions粘塑性模型[11]的非算法依賴的結(jié)構(gòu)穩(wěn)定性理論[12]。并且可以進(jìn)一步證明，時(shí)效變形加固理論的材料本構(gòu)關(guān)系可由Duvaut-Lions模型拓展到非線性、非平衡態(tài)熱力學(xué)框架Rice正則結(jié)構(gòu)，從而適應(yīng)裂隙巖體的多尺度特點(diǎn)[12]。

1 巖體結(jié)構(gòu)非平衡演化及控制理論

傳統(tǒng)的變形加固理論在數(shù)值上采用超出屈服面的殘余不平衡節(jié)點(diǎn)力表征結(jié)構(gòu)的非平衡狀態(tài)[13-14]，這種靜態(tài)分析忽略了巖體結(jié)構(gòu)變形的時(shí)間效應(yīng)。為克服上述缺點(diǎn)，楊強(qiáng)、劉耀儒等基于Duvaut-Lions粘塑性模型建立了考慮時(shí)效的變形加固理論[10-12]。

1.1 Duvant-Lions粘塑性模型

Duvaut-Lious模型采用率形式本構(gòu)關(guān)系，在考慮關(guān)聯(lián)流動(dòng)法則情況時(shí)，一個(gè)變形體在小變形情況時(shí)彈-粘塑性材料率形式的應(yīng)力-應(yīng)變關(guān)系如下

(1)

(2)

1.2 時(shí)效的最小塑性余能原理

考慮時(shí)效的變形加固理論中，定義E為結(jié)構(gòu)的塑性余能范數(shù)，可用下式表示

(3)

(4)

式(4)表明，在恒定外部條件下，理想粘塑性結(jié)構(gòu)始終朝向塑性余能范數(shù)減小的方向變形，直到塑性余能范數(shù)達(dá)到最小值，此后塑性余能范數(shù)保持恒定，此即為最小塑性余能原理。

對(duì)于任意時(shí)刻，塑性余能范數(shù)E越大，則表明過應(yīng)力越大，結(jié)構(gòu)離平衡態(tài)越遠(yuǎn)。換句話說，塑性余能度量了結(jié)構(gòu)當(dāng)前非平衡演化狀態(tài)到平衡態(tài)的距離(E滿足Euclid距離的所有條件)。同時(shí)，式(3)表明，E是半正定的，當(dāng)且僅當(dāng)結(jié)構(gòu)達(dá)到平衡態(tài)時(shí)(E=0)，也就是結(jié)構(gòu)最終停止變形，達(dá)到穩(wěn)定。式(4)表明，在結(jié)構(gòu)變形過程中，E隨時(shí)間是不增加的，即結(jié)構(gòu)始終不會(huì)背離平衡態(tài)演化。結(jié)構(gòu)最終發(fā)揮最大的自承載力抵抗外荷載，使其自身盡量趨于平衡態(tài)。

塑性余能范數(shù)隨時(shí)間演化見圖1。圖中，t為時(shí)間；Emin為塑性余能范數(shù)E的最小值；tls為塑性余能范數(shù)達(dá)到最小值時(shí)的時(shí)刻。從圖1可知，結(jié)構(gòu)演化分為2個(gè)階段：第一階段稱為自平衡演化階段；第二階段稱為恒定演化階段。2個(gè)演化階段的分界點(diǎn)定義為極限恒定演化狀態(tài)，其對(duì)應(yīng)時(shí)刻即為tls。

圖1 塑性余能范數(shù)隨時(shí)間演化示意

(5)

基于最小塑性余能原理，塑性余能范數(shù)E和不平衡力分布U分別實(shí)現(xiàn)了結(jié)構(gòu)非平衡演化的穩(wěn)定評(píng)價(jià)與失穩(wěn)控制。該原理直觀表述為：結(jié)構(gòu)通過自身變形達(dá)到自承載力最大化、所需加固力最小化的狀態(tài)，不平衡力即是結(jié)構(gòu)的最優(yōu)加固力。

(6)

結(jié)構(gòu)失穩(wěn)控制：阻止結(jié)構(gòu)失穩(wěn)的加固力為U，從塑性余能范數(shù)最小的角度，根據(jù)最優(yōu)不平衡力Umin確定的加固力具有全局最優(yōu)性，或者說具有最高加固效率。通過有限元分析獲得精確到節(jié)點(diǎn)的全局最優(yōu)加固力配置是變形加固理論的最大特點(diǎn)，也是其較快得到工程應(yīng)用的首要原因。

變形加固理論是基于有限元法的結(jié)構(gòu)穩(wěn)定性理論，但它并未對(duì)有限元法本身進(jìn)行實(shí)質(zhì)性更改。塑性余能與不平衡力均為有限元法固有中間變量，特別是不平衡力，它就是有限元法中的殘余力(residual force)，即增量節(jié)點(diǎn)平衡方程的右端項(xiàng)。從計(jì)算角度講，變形加固理論與傳統(tǒng)有限元結(jié)構(gòu)分析的唯一區(qū)別在于額外計(jì)算并記錄塑性余能與不平衡力這2項(xiàng)結(jié)構(gòu)響應(yīng)，它們與位移、應(yīng)力、塑性區(qū)等結(jié)構(gòu)響應(yīng)的地位是并列的，只是塑性余能與不平衡力反映了結(jié)構(gòu)失穩(wěn)的根源。

綜上，變形加固理論并非更改有限元法，而是深入挖掘有限元法易被忽略的內(nèi)涵。

1.3 與剛體極限平衡法的對(duì)比

對(duì)于單滑塊單自由度的抗滑穩(wěn)定問題，變形加固理論與剛體極限平衡法的分析結(jié)果是一致的。單滑塊穩(wěn)定分析見圖2。圖2中，T為滑動(dòng)力，R為抗滑力，N為法向力，Q為加固力。從圖2可知，假定抗滑條件滿足Mohr-Coulomb準(zhǔn)則，在法向力N的作用下，最大可能的抗滑力R=fN+cA。其中，f和c分別為摩擦系數(shù)和粘聚力，A為滑面面積。當(dāng)T=R，滑塊處于臨界穩(wěn)定狀態(tài)。當(dāng)T>R時(shí)，滑塊無(wú)法保持平衡，需要施加與滑動(dòng)方向相反的加固力Q來保持穩(wěn)定，Q=-(T-R)。

圖2 單滑塊穩(wěn)定分析

從變形加固理論的觀點(diǎn)分析，當(dāng)抗滑條件(R≤fN+cA)和平衡條件(T=R)無(wú)法同時(shí)滿足，滑塊失去穩(wěn)定。由于滑面真實(shí)抗滑力的最大值為R=fN+cA，因此上述2個(gè)條件中，平衡條件不能滿足，即T>R。此時(shí)作用在滑塊上的不平衡力為F，F(xiàn)=T-R=-Q。根據(jù)潘家錚“最大最小原理”，此時(shí)的滑動(dòng)面是唯一的滑面，失穩(wěn)時(shí)滑面將發(fā)揮最大的抗滑作用。根據(jù)牛頓第三定律，不平衡力量值正比于滑塊加速度。這就表明，滑塊失穩(wěn)時(shí)不平衡力趨于最小，而Q就是最優(yōu)的加固力。

1.4 不平衡力性質(zhì)的討論

不平衡力是彈塑性有限元計(jì)算中的一個(gè)中間變量，本質(zhì)是一種虛擬力，并不會(huì)在實(shí)際工程中出現(xiàn)。如果有限元計(jì)算中結(jié)構(gòu)出現(xiàn)不平衡力，則實(shí)際工程中會(huì)通過損傷、開裂等方式使不平衡力得到釋放。根本原因是由于在常用的基于位移有限元分析中，結(jié)構(gòu)連續(xù)性和變形協(xié)調(diào)條件被人為強(qiáng)制滿足，無(wú)法反映結(jié)構(gòu)局部開裂帶來的不連續(xù)性和非協(xié)調(diào)性。理想彈塑性材料是無(wú)損材料，出現(xiàn)不平衡力意味著結(jié)構(gòu)無(wú)法繼續(xù)滿足在連續(xù)、完整、無(wú)損的條件下收斂到一個(gè)穩(wěn)定解，此時(shí)的不平衡力可從加固和破壞2個(gè)角度理解：

圖3 預(yù)設(shè)裂紋的試件在逐步加載過程中開裂破壞

圖4 預(yù)設(shè)裂紋的試件在逐步加載過程中不平衡力矢量

(1)通過施加大小相等、方向相反的加固力來消除不平衡力，進(jìn)而結(jié)構(gòu)連續(xù)、完整、無(wú)損性可得以維持。

(2)如不施加加固力，結(jié)構(gòu)只能通過在不平衡力出現(xiàn)處產(chǎn)生損傷開裂(改變本構(gòu)關(guān)系)來消除不平衡力，以此來弱化強(qiáng)制的變形協(xié)調(diào)條件。因此，不平衡力實(shí)質(zhì)是一種材料損傷開裂程度的等效表述。

以不平衡力表征損傷類似于損傷力學(xué)中以有效應(yīng)力表征損傷，但這種方法在求解過程中只需提供常見的剛度和強(qiáng)度參數(shù)，無(wú)需額外增加損傷演化相關(guān)參數(shù)。該方法具有一階無(wú)條件收斂的特性，特別適用于大規(guī)模的復(fù)雜巖體工程計(jì)算。

圖3為單軸壓縮含預(yù)制裂紋的試件開裂破壞過程[15](單軸壓力由0逐級(jí)增加到2 MPa)。單軸壓力為1.6 MPa時(shí)有明顯的開裂，開裂首先出現(xiàn)在預(yù)制裂紋的尖端，隨著繼續(xù)加載，新增裂縫近似垂直于預(yù)制裂紋；單軸壓力為2.0 MPa時(shí)，相鄰預(yù)制裂紋之間的裂縫貫穿，試件最終破壞。圖4為對(duì)應(yīng)的不平衡力矢量圖[15]。從圖4可知，單軸壓力為1 MPa時(shí)，不平衡力首先出現(xiàn)預(yù)制裂紋的尖端；繼續(xù)增加荷載，不平衡力擴(kuò)展如圖4的粗實(shí)線，與預(yù)制裂紋的走向近似垂直。綜上，不平衡力的出現(xiàn)和擴(kuò)展可用于近似模擬新裂縫的出現(xiàn)與擴(kuò)展，不平衡力量值與裂縫開裂程度近似正相關(guān)。

2 特高拱壩建基面卸荷松弛研究

白鶴灘水電站特高拱壩壩高289 m，樞紐區(qū)地質(zhì)條件十分復(fù)雜，地應(yīng)力高，橫河向剖面見圖6。在樞紐區(qū)發(fā)育有較多的結(jié)構(gòu)面。其中，LS331、F17、LS3319、LS3318等在拱間槽出露。壩基處于高地應(yīng)力區(qū)，且側(cè)壓力系數(shù)較高。2014年12月，左岸建基面開挖至628 m高程，巡視發(fā)現(xiàn)左岸拱間槽松弛卸荷嚴(yán)重，部分?jǐn)鄬映霈F(xiàn)錯(cuò)動(dòng)現(xiàn)象，工程緊急停工，進(jìn)行卸荷松弛研究。

為精細(xì)模擬白鶴灘水電站左岸建基面開挖卸荷松弛情況，采用大范圍、小單元精細(xì)化網(wǎng)格計(jì)算。模擬了從950 m高程(天然邊坡)連續(xù)分步開挖至538 m高程共19步開挖工況，模擬了C3等主要結(jié)構(gòu)面。

圖5 白鶴灘水電站左岸壩基開挖剖面

圖6 典型剖面開挖不平衡力矢量

建基面開挖至630 m高程，不平衡力較大的結(jié)構(gòu)面為L(zhǎng)S331(558 990 kN)、LS3319(234 680 kN)，各斷層與錯(cuò)動(dòng)帶的不平衡力和屈服區(qū)主要集中在結(jié)構(gòu)面交匯處及鄰近邊坡或開挖面出露區(qū)域。斷層F17與錯(cuò)動(dòng)帶C3-1在交匯區(qū)相互作用，引起拱間槽下游側(cè)720 m高程附近臨空面部分巖體屈服，這與現(xiàn)場(chǎng)觀測(cè)中拱間槽下游側(cè)坡沿F17產(chǎn)生裂縫相符。LS3319在下游側(cè)545～570 m高程之間有比較大的塑性區(qū)，不平衡力也集中于此處。結(jié)構(gòu)面卸荷松弛的現(xiàn)場(chǎng)監(jiān)測(cè)結(jié)果與計(jì)算結(jié)果比較吻合。變形加固理論預(yù)測(cè)成果表明，從630 m高程繼續(xù)向下開挖時(shí)應(yīng)重點(diǎn)關(guān)注LS331、LS3319、F17的松弛破壞。

圖6為從590 m高程開挖至538 m高程過程中，建基面某典型剖面的不平衡力矢量圖。計(jì)算結(jié)果也表明，從590 m高程開挖至570 m高程，LS331的不平衡力由825 690 kN激增到1 179 710 kN，而這也是陡坎成形和LS331開始出露臨空的過程。開挖到560 m高程以后，LS331的總不平衡力開始降低。因此，建議將左岸建基面570～590 m高程陡坎坡度放緩，并加強(qiáng)對(duì)LS331出露段附近勘測(cè)與監(jiān)測(cè)，該建議被建設(shè)方采納。

3 特高拱壩整體穩(wěn)定性研究

變形加固理論中采用超載法評(píng)價(jià)拱壩穩(wěn)定性，對(duì)每1個(gè)超載倍數(shù)K通過非線性迭代計(jì)算出相應(yīng)的塑性余能范數(shù)E。以K—E曲線評(píng)價(jià)拱壩穩(wěn)定性，可以近似說明巖體結(jié)構(gòu)工程在局部失穩(wěn)后的整體穩(wěn)定性情況。引入工程類比法，將某拱壩計(jì)算的K—E曲線與之前類似拱壩的K—E曲線進(jìn)行對(duì)比，可以說明拱壩整體穩(wěn)定性的情況。圖7為國(guó)內(nèi)主要高拱壩基礎(chǔ)和壩體超載倍數(shù)與塑性余能范數(shù)關(guān)系[16]。

圖7 主要高拱壩超載倍數(shù)與塑性余能范數(shù)關(guān)系

基于地質(zhì)力學(xué)模型試驗(yàn)和非線性數(shù)值模擬，可以采用3K安全系數(shù)對(duì)高拱壩進(jìn)行穩(wěn)定性評(píng)價(jià)[16]，這和基于變形加固理論的穩(wěn)定分析是對(duì)應(yīng)的：

(1)K1為起裂安全系數(shù)。變形加固理論中指出，如果外荷載大于結(jié)構(gòu)極限抗力，超出屈服面的應(yīng)力無(wú)法完全轉(zhuǎn)移，即不平衡力不能迭代到0，在不平衡力的作用下，巖體結(jié)構(gòu)進(jìn)行非平衡演化，出現(xiàn)局部起裂。巖體結(jié)構(gòu)首先出現(xiàn)不平衡力的位置就是首先起裂破壞的位置；不平衡力越大，表明此處的開裂越嚴(yán)重。

(2)K2為非線性變形安全系數(shù)。當(dāng)超載倍數(shù)超過一定值時(shí)，屈服區(qū)體積增加失去線性，認(rèn)為此時(shí)結(jié)構(gòu)整體進(jìn)入非線性狀態(tài)。在高拱壩整體進(jìn)入非線性階段之前，巖體結(jié)構(gòu)的塑性屈服區(qū)體積相比于整個(gè)結(jié)構(gòu)的體積很小，可近似認(rèn)為屈服區(qū)的體積與超載倍數(shù)為線性關(guān)系。所以，可通過超載倍數(shù)與屈服體積曲線失去線性關(guān)系作為結(jié)構(gòu)整體進(jìn)入非線性階段的判據(jù)。

(3)K3為極限承載力的安全系數(shù)。到達(dá)極限承載力后，拱壩全面進(jìn)入失穩(wěn)，不能繼續(xù)超載，結(jié)構(gòu)整體破壞；在非線性有限元計(jì)算中，表現(xiàn)為高拱壩的屈服區(qū)體積迅速增加，壩踵等薄弱區(qū)域發(fā)生較大的塑性流動(dòng)。實(shí)踐表明，地質(zhì)力學(xué)模型試驗(yàn)中不能繼續(xù)加載和拱壩塑性屈服區(qū)快速增加的2個(gè)極限承載力判據(jù)是一致的。

4 壩肩高邊坡的穩(wěn)定分析

大崗山水電站壩址區(qū)巖脈、擠壓破碎帶、斷層和6組節(jié)理裂隙發(fā)育。在右岸壩肩區(qū)邊坡開挖過程中，先后于2009年5月和8月份發(fā)生2次開裂，經(jīng)過緊急支護(hù)處理后，邊坡變形得到控制。然而2009年9月1日，在進(jìn)行拱肩槽邊坡開挖時(shí)，伴隨著一聲巨響，在右岸1 135 ～1 165 m高程出現(xiàn)裂縫，邊坡已開始出現(xiàn)持續(xù)的宏觀變形，導(dǎo)致工程停工近1年。由于大崗山水電站右岸邊坡斷層、巖脈密布，且初始地應(yīng)力較高，開挖所引起的邊坡變形破壞問題十分突出。

采用三維非線性有限元對(duì)開挖前后大崗山水電站右岸邊坡的整體穩(wěn)定性進(jìn)行分析，以塑性余能范數(shù)來判定整體穩(wěn)定性，以不平衡力來確定失穩(wěn)部位及所需加固力。有限元計(jì)算網(wǎng)格及模擬的主要結(jié)構(gòu)面見圖8。天然邊坡與開挖邊坡在降強(qiáng)計(jì)算中，折減系數(shù)—塑性余能范數(shù)曲線見圖9。開挖前后主要結(jié)構(gòu)面的不平衡力見表1。

圖8 計(jì)算模型及主要結(jié)構(gòu)面示意

圖9 塑性余能隨折減系數(shù)變化

表1 開挖前后邊坡主要斷層及巖脈不平衡力對(duì)比10 kN

由圖9和表1可知，開挖前后邊坡的整體安全度均在1.5左右，發(fā)生整體失穩(wěn)的可能性不大。對(duì)比開挖前后主要結(jié)構(gòu)面的不平衡力，在整個(gè)邊坡開挖及降強(qiáng)過程中，緩傾角斷層f231不平衡力均有所增加，特別是拱肩槽邊坡開挖時(shí)，需要在施工時(shí)重點(diǎn)關(guān)注。而β209不只在開挖過程中產(chǎn)生10 000 kN左右不平衡力，且在降強(qiáng)過程中不平衡力增長(zhǎng)迅速，故可以初步判斷，開挖過程中的可能破壞模式主要為β209的開裂引起的沿著緩傾角斷層f231的局部錯(cuò)動(dòng)。

5 結(jié) 語(yǔ)

處于臨界平衡態(tài)的巖體結(jié)構(gòu)在受到開挖卸荷、水荷載變化等工程擾動(dòng)后，平衡態(tài)被打破，結(jié)構(gòu)進(jìn)入非平衡演化過程。在向非平衡演化過程中，會(huì)伴隨損傷和大變形，直接影響到巖體結(jié)構(gòu)工程的長(zhǎng)期安全。時(shí)效變形加固理論基于非平衡演化規(guī)律的研究，建立了一個(gè)能用于結(jié)構(gòu)穩(wěn)定性評(píng)價(jià)的理論框架。時(shí)效塑性余能范數(shù)是評(píng)價(jià)結(jié)構(gòu)整體偏離平衡態(tài)的程度，可用于結(jié)構(gòu)穩(wěn)定性演化的評(píng)價(jià)；時(shí)效不平衡力是偏離程度在節(jié)點(diǎn)上的體現(xiàn)，分布和大小表明了結(jié)構(gòu)破壞的程度，可以用于指導(dǎo)巖體結(jié)構(gòu)工程的加固設(shè)計(jì)。該理論目前在高拱壩、高邊坡等工程的穩(wěn)定性評(píng)價(jià)和加固設(shè)計(jì)中得到了廣泛的應(yīng)用。在特高拱壩建基面開挖卸荷松弛、整體穩(wěn)定性、長(zhǎng)期運(yùn)行安全性等水工巖石力學(xué)的關(guān)鍵問題上，均得到了較好的驗(yàn)證。