胡少偉

（南京水利科學(xué)研究院材料結(jié)構(gòu)研究所，210024，南京）

在大壩設(shè)計和安全評價中考慮地震作用的影響大約始于1930年前后，當(dāng)時主要按擬靜力法進行大壩的抗震設(shè)防，將地震作用以一個等價的水平作用力進行代替，大多數(shù)壩的設(shè)計地震加速度采用0.1 g，動水壓力按Westergaard的近似公式進行計算。1970年前后，大壩的抗震設(shè)計有了新的發(fā)展，這期間發(fā)生的多次地震引起了水庫震害，如我國新豐江大壩的震害、美國San Fernando地震引起的Pacoima拱壩震害，促使線彈性動力分析方法在大壩抗震設(shè)計中逐步得到應(yīng)用。2000年前后，由于結(jié)構(gòu)抗震設(shè)計技術(shù)的發(fā)展，同時經(jīng)過多年的實踐，研究者和設(shè)計者們均認(rèn)識到，彈性動力分析給出的峰值應(yīng)力不能提供足夠的信息以判斷不同地震所引起的震害的嚴(yán)重性和對大壩安全的影響，同時對大壩抗震安全起關(guān)鍵作用的混凝土容許拉應(yīng)力，各國甚至各個設(shè)計單位都沒有取得共識。本文通過對比中美兩國在混凝土壩抗震設(shè)計方面的理論，分析兩國在抗震設(shè)計方面的差異，以期為完善我國混凝土壩抗震設(shè)計理論提供有益的借鑒。

一、大壩抗震設(shè)防標(biāo)準(zhǔn)

美國大壩采用運行基本地震OBE （Operating Basis Earthquake，相應(yīng)概率水準(zhǔn)為100年基準(zhǔn)期超越概率50%）和最大設(shè)計地震MDE（Maximum Design Earthquake）或最大可信地震 MCE（Maximum Credible Earthquake）兩級設(shè)防，并提出了相應(yīng)于不同設(shè)防水準(zhǔn)時大壩的性能要求。其中，OBE是從工程運行角度提出來的，對應(yīng)于運行期工程保護期望水平的最大地震，這時主要考慮地震引發(fā)的結(jié)構(gòu)損壞、機械破損和經(jīng)濟損失，設(shè)防要求是震后結(jié)構(gòu)易修復(fù)、設(shè)備可繼續(xù)運行；MDE或MCE主要是從避免引發(fā)嚴(yán)重次生災(zāi)害的角度提出來的，設(shè)防要求為大壩不發(fā)生災(zāi)難性破壞，如不致使庫水下泄失控。當(dāng)大壩失事可能導(dǎo)致危及人身安全的嚴(yán)重后果時取MDE為MCE，否則MDE一般小于MCE。MCE為壩址區(qū)可能發(fā)生的最大地震，并假定發(fā)生在離壩址最近的斷層點上。

我國的大壩采用最大設(shè)計地震MDE一級設(shè)防。對于有利地段，建筑物一般采用基本烈度設(shè)防 （相當(dāng)于50年基準(zhǔn)期超越概率10%），對于甲類設(shè)防的大壩在基本烈度基礎(chǔ)上提高一度設(shè)防（相當(dāng)于100年基準(zhǔn)期超越概率2%）。性能要求為如有局部損壞，經(jīng)一般處理后仍可正常運行。100年基準(zhǔn)期超越概率2%的設(shè)防水準(zhǔn)接近國外一些國家提出的最大可信地震（MCE）的水平，而其性能目標(biāo)又與運行基本地震（OBE）的要求相近。

因此有必要開展我國大壩抗震設(shè)計分級設(shè)防的合理性和可行性研究，必要時提出各級設(shè)防水準(zhǔn)及相應(yīng)的性能目標(biāo)和配套的分析方法，使大壩的抗震設(shè)計更為合理，既保障了其安全性，同時又符合經(jīng)濟性原則。另外對于修建在抗震不利地段或城市上游附近的大壩以及300 m級超高壩的抗震設(shè)防標(biāo)準(zhǔn)也需開展專門研究。

二、混凝土材料的動態(tài)性能

美國墾務(wù)局的Raphael在西方5座混凝土壩中鉆孔取樣進行試驗，在0.05s的時間內(nèi)加載到極限強度（相當(dāng)于大壩5Hz的振動頻率），得出動態(tài)抗壓強度較靜強度平均提高31%，直接拉伸強度平均提高66%，劈拉強度平均提高45%，試驗結(jié)果有一定離散性。在美國等許多國家的設(shè)計標(biāo)準(zhǔn)中，一般參考Raphael的試驗結(jié)果，將地震作用下混凝土的動態(tài)強度較靜態(tài)值提高一個百分比。

我國現(xiàn)行《水工建筑物抗震設(shè)計規(guī)范》（DL 5093—2000）也作了類似規(guī)定。這些規(guī)定不考慮大壩和地震激勵的具體情況，一律將地震作用下混凝土的強度或彈性模量等動態(tài)特性參數(shù)較靜態(tài)常數(shù)提高一個固定百分比，作法顯得過于粗略。實際上在地震作用下，不同的大壩、不同的部位、不同的瞬時，其應(yīng)變速率的變化各不相同，從而混凝土的強度和剛度特性將隨著激勵地震波的特性而變化。其變化規(guī)律還和混凝土的原材料特性密切相關(guān)。另外，在混凝土壩的強度校核中，作用效應(yīng)綜合了地震作用和靜態(tài)作用的影響，僅用動態(tài)強度來評價大壩靜動疊加響應(yīng)也不盡合理，尤其是對于靜態(tài)應(yīng)力水平較高的大壩。目前，國內(nèi)開展了在不同靜態(tài)預(yù)載水平下大壩混凝土動態(tài)彎拉強度和彈性模量隨加載速率（應(yīng)變率）變化的全級配混凝土試驗，初步成果認(rèn)為，在無靜態(tài)預(yù)載時，相對于靜態(tài)強度的動態(tài)提高率小于規(guī)范基于濕篩小試件成果所采用的30%，約為17%；在初始靜載達到80%混凝土強度以內(nèi)時動態(tài)提高率隨初始靜載水平的提高而有所提高；初始靜載超過80%混凝土強度后動態(tài)提高率明顯下降；動態(tài)彈性模量較靜態(tài)值略有提高 （約為5%），但明顯小于彎拉強度的提高率。由于試驗試件數(shù)量較少，關(guān)于全級配大壩混凝土試件的動態(tài)抗力取值的一般規(guī)律還有待于進一步研究。

三、壩—地基相互作用

在壩體和地基的動力相互作用中，對地基影響的考慮，與靜力問題不同之處，在于它不僅限于計入地基的柔度，而且要考慮地基質(zhì)量的慣性和其因振動能量自遠處逸散而形成的所謂輻射阻尼?；炷翂闻c地基間的動力相互作用問題的研究已有40年歷史，國內(nèi)也進行了近30年的研究工作，在許多領(lǐng)域取得研究成果，分析這一問題最有效的方法是子結(jié)構(gòu)法，其中關(guān)鍵問題是建立基巖的阻抗矩陣。在重力壩的平面分析中，半平面地基阻抗矩陣可從波動方程得到理論解，對于空間問題只能借助于數(shù)值方法得到近似解。在一些拱壩動力程序中更多采用美國Clough教授提出的無質(zhì)量地基假定，但它有很大的近似性。建立基巖的阻抗矩陣時，一般先求出與壩基面相對應(yīng)的基巖動力柔度系數(shù)。我國這一領(lǐng)域的研究相對比較深入，清華大學(xué)提出了有限元—邊界元—無窮邊界元的計算模型，中國水利水電科學(xué)研究院則發(fā)展了廖振鵬的透射邊界計算模型，對有關(guān)相互作用對高拱壩地震變形和應(yīng)力的影響有了更深入的了解。各種計算模型都有一定特點，也有一定局限性，目前還比較缺乏實測資料的檢驗。大連理工大學(xué)發(fā)展了Wolf與Song的阻尼影響抽取法的計算模型以及比例邊界有限元的模型，并研究了地基剛度變化和地基不均勻性對拱壩和重力壩地震響應(yīng)的影響，獲得了一些新的認(rèn)識。研究結(jié)果發(fā)現(xiàn)采用無質(zhì)量地基的計算模型時，壩頂部的動力放大系數(shù)基本上不隨地基剛度的變化而變化，但考慮無限地基動力相互作用的影響以后，隨著地基與壩體彈性模量比值的不同，壩頂部的動力放大系數(shù)將隨之發(fā)生比較大的變化。地基愈柔性，特別是當(dāng)?shù)鼗膭偠?彈性模量)低于壩體的剛度時，地震響應(yīng)有較大幅度的降低。

四、混凝土壩抗震安全評價準(zhǔn)則

我國的水工建筑物抗震設(shè)計規(guī)范標(biāo)準(zhǔn)規(guī)定以容許應(yīng)力建立強度條件，以最大拉應(yīng)力安全系數(shù)來控制高拱壩的抗震安全性，它不容許拱壩在地震中出現(xiàn)裂縫。采用一級設(shè)防標(biāo)準(zhǔn)，對基本烈度（50年超越概率10%，重現(xiàn)期475年）為Ⅶ、Ⅷ、Ⅸ度區(qū)的場地，設(shè)計地震加速度分別取為0.1 g、0.2 g和0.4 g。對重要大壩則需將設(shè)計地震加速度的水準(zhǔn)提高到100年超越概率2%、重現(xiàn)期4 950年。

以美國為代表的一些國家采用了兩級地震設(shè)防標(biāo)準(zhǔn)。美國墾務(wù)局1970年以前大壩設(shè)計地震加速度采用0.1 g，1974年以后提出設(shè)計基準(zhǔn)地震DBE與最大可信地震MCE兩級設(shè)防的概念。美國大壩委員會1985年起草并經(jīng)國際大壩委員會1989年公布的《大壩地震系數(shù)選擇導(dǎo)則》，明確了使用安全運行地震動OBE與最大設(shè)計地震動MDE兩級設(shè)防的地震動參數(shù)選擇原則。按照這一準(zhǔn)則，在安全運行地震OBE作用時，大壩應(yīng)能保持運行功能，所受震害易于修復(fù)。故一般可進行彈性分析，并采用容許應(yīng)力準(zhǔn)則。在最大設(shè)計地震MDE作用時，要求大壩至少能保持蓄水能力。這表示可容許大壩出現(xiàn)裂縫，但不影響壩的整體穩(wěn)定，不發(fā)生潰壩。同時大壩的泄洪設(shè)備可以正常工作，震后能放空水庫。OBE一般選為100年內(nèi)超越概率50%（重現(xiàn)期145年）的地震動水平，以Housner為首的美國大壩安全委員會則建議DBE的重現(xiàn)期為200年，經(jīng)過經(jīng)濟上合理性的論證時還可適當(dāng)延長。關(guān)于MDE的概率水準(zhǔn)或重現(xiàn)期，沒有作明確確定。值得注意的是，MDE的決定一般都和大壩的失事后果相聯(lián)系，只對特別重要的壩才令MDE等于MCE。確定MCE一般有地質(zhì)構(gòu)造法和概率法兩種，國際大壩委員會的導(dǎo)則認(rèn)為，就目前的認(rèn)識水平而言，不可能明確規(guī)定必須采用哪種方法。建議同時采用兩種方法，并應(yīng)用工程經(jīng)驗進行判斷。

中美兩國大壩抗震設(shè)防處理方法和地震動參數(shù)的取值不同，原因在于各國國情不同，材料強度的控制標(biāo)準(zhǔn)不同，施工質(zhì)量的可靠程度不同等等。我國不容許壩在地震作用下出現(xiàn)裂縫的評價準(zhǔn)則顯得要求過高，對于復(fù)雜的高拱壩來講常常難以做到。一方面，出現(xiàn)一些裂縫可能并不影響壩的安全運行；另一方面,高拱壩是超靜定結(jié)構(gòu)，局部拉應(yīng)力過大引起的開裂是拱壩放松約束來調(diào)整內(nèi)力以適應(yīng)變化的正常反應(yīng)，在破壞之前必定存在一個不斷開裂、不斷靜定化的過程，在這個過程中拱壩的潛力未盡。國外在保持水庫蓄水能力的情況下一般容許發(fā)生一定程度震害的大壩抗震設(shè)防指導(dǎo)思想值得我們深入研究；我國規(guī)范對重要的大壩、高壩必須進行專門研究的做法值得堅持和推廣；對于重要大壩采取兩級或多級地震設(shè)防標(biāo)準(zhǔn)的研究可使大壩的抗震設(shè)計更為合理,符合安全性和經(jīng)濟性兼顧的原則。

五、結(jié) 語

大壩的抗震性評價將重點探索簡單而實用的動力響應(yīng)分析方法來取代擬靜力法，重點研究壩—地基體系相互作用、地震動輸入機制、加強混凝土材料強度隨應(yīng)變速率變化的規(guī)律分析并建立合理的計算模型等研究，進一步深入對大壩抗震評價準(zhǔn)則確定原則的研究、引入美國兩級或多級抗震設(shè)防的理念和性能設(shè)計的思想，提高我國大壩抗震設(shè)計水平。

[1]胡少偉.大壩結(jié)構(gòu)抗震安全評價與減震技術(shù)亟待加強 [J].中國水利，2008（11）.

[2]林皋.混凝土大壩抗震安全評價的發(fā)展趨向 [J].防災(zāi)減災(zāi)工程學(xué)報，2006，26（1）.

[3]陳厚群.高拱壩抗震設(shè)計研究進展[J].中國水利，2000（9）.

[4]United States Committee on Large Dams （USCOLD）.Updated Guidelines for Selecting Seismic Parameters for Dam Progects[M].Denver，United States，1999.

[5]林皋，陳健云.混凝土大壩的抗震安全評價[J].水利學(xué)報，2001（2）.

[6]Kolsky，H.，Sress waves in solids[M].New York:Dover Publications，1963.

[7]金峰，張楚漢，王光綸.結(jié)構(gòu)地基相互作用的有限元—邊界元—無限邊界元耦合模型 [J].清華大學(xué)學(xué)報，1993，33（2）.

[8]杜修力，陳厚群，侯順載.拱壩系統(tǒng)三維非線性地震波動分析[J].地震工程與工程振動，1996，16（3）.

[9]李建波，陳健云，林皋.無限地基-結(jié)構(gòu)動力相互作用分析的分區(qū)遞歸時域算法研究[J].工程力學(xué)，2005，22（5）.

[10]胡少偉，苗同臣.結(jié)構(gòu)振動理論及其應(yīng)用[M].北京：中國建筑工業(yè)出版社，005.

[11]周建平，陳觀福，黨林才.我國高壩抗震安全評價的現(xiàn)狀與挑戰(zhàn) [J].水利學(xué)報，2007（增刊）.

[12]胡少偉，陸俊.混凝土材料非線性對拱壩抗震安全評定的影響研究[J].水力發(fā)電學(xué)報，2009（5）.

[13]Hilbert， LBL， Jr.YW， Cook NGW.A new discontinuousfinite elementmethod for interaction of manydeformablebodiesin geomechanices[A].In:Pro 8th int cof comp meth，Adv.Geomech，1994，931～936.

[14]Hu Shaowei.Structural Integrity Assessment and Large Deformation Analysis Theory[M].Yellow River Conservancy Press.2006.