高面板堆石壩面板地震動力反應(yīng)特性大型 振動臺模型試驗研究

楊正權(quán) 劉小生 周國斌 劉啟旺 趙劍明 楊玉生

摘要： 混凝土面板是面板堆石壩整個防滲體系的核心結(jié)構(gòu)，其結(jié)構(gòu)安全性是整個大壩系統(tǒng)正常運(yùn)行的基本保障。強(qiáng)震作用下，面板產(chǎn)生動應(yīng)力應(yīng)變反應(yīng)與靜力作用效應(yīng)相疊加，可能會對面板的結(jié)構(gòu)安全性產(chǎn)生重大影響，故面板的地震動力反應(yīng)特性一直都是大壩抗震設(shè)計關(guān)注的重點問題。通過高面板堆石壩大型地震模擬振動臺模型試驗，研究面板在地震作用下的各向動應(yīng)變反應(yīng)過程及其極值空間分布規(guī)律，再根據(jù)面板動應(yīng)變反應(yīng)過程和材料力學(xué)理論分析面板的動應(yīng)力反應(yīng)過程及其極值空間分布特性，并利用振動臺模型試驗成果推求研究原型壩混凝土面板的動力反應(yīng)特性。研究表明：振動臺模型試驗實測面板三向應(yīng)變反應(yīng)過程及其推求的主應(yīng)變、主應(yīng)力過程符合對面板動力反應(yīng)特性的一般認(rèn)識，振動臺模型試驗可作為研究面板動力反應(yīng)特性的重要手段；沿豎直方向上，面板的大、小主應(yīng)變極值均出現(xiàn)在0.6倍壩高左右的高程上；和應(yīng)變分布規(guī)律基本對應(yīng)，面板的大、小主應(yīng)力和最大剪應(yīng)力極值同樣分布在0.6倍壩高左右的高程上；面板在地震中所受動拉應(yīng)力水平總體低于動壓應(yīng)力。

關(guān)鍵詞： 地震反應(yīng)； 高面板堆石壩； 混凝土面板； 大型振動臺模型試驗

中圖分類號： TV312； TV641.4文獻(xiàn)標(biāo)志碼： A文章編號： 1004-4523（2018）02-0209-10

DOI：10.16385/j.cnki.issn.1004-4523.2018.02.003

引 言

面板堆石壩是隨著近年來碾壓堆石壩壩體施工技術(shù)快速進(jìn)步而迅速發(fā)展起來的一種土石壩壩型。面板堆石壩由于壩體體型相對“單薄”而使得其土石方量較其他土石壩壩型少得多，防滲體材料為人工材料混凝土，不受當(dāng)?shù)胤罎B材料儲量的限制，且面板系統(tǒng)的體量小、防滲效果好，因而受到了現(xiàn)代壩工領(lǐng)域的青睞，國內(nèi)外已有大量的面板堆石壩建成，并有數(shù)量眾多的高或者超高面板堆石壩正在建設(shè)[1-2]。面板堆石壩的工程優(yōu)勢主要是在于面板防滲體系的使用，也正是為此，面板的安全性和防滲效果對整個面板堆石壩結(jié)構(gòu)體系的安全性和效益發(fā)揮起到了決定性的作用。尤其是在地震工況下，地震動作用所產(chǎn)生的面板動應(yīng)力和靜力作用產(chǎn)生的應(yīng)力相疊加，很大可能使得面板的受力狀態(tài)趨于惡化，對整個混凝土面板系統(tǒng)的安全性是一個巨大挑戰(zhàn)。靜力作用下，國內(nèi)外已有數(shù)個堆石壩混凝土面板發(fā)生破壞而影響工程安全運(yùn)行的案例，如巴西的Campos Novos面板堆石壩和Barra Grande面板堆石壩發(fā)生了面板垂直縫擠壓破壞，中國的三板溪面板堆石壩和布西面板堆石壩發(fā)生了水平向施工縫擠壓破壞等[3]。也有由于地震作用而引起面板局部破壞的工程案例，如2008年汶川地震中的紫坪鋪面板堆石壩，面板發(fā)生了局部錯臺和擠壓破壞[4]。面板在地震作用下的破壞性態(tài)受其地震動力反應(yīng)性態(tài)所控制，因而面板系統(tǒng)在地震作用下的地震動力反應(yīng)特性，就成為了工程技術(shù)人員關(guān)心的問題。

由于土石料動力變形和強(qiáng)度特性較為復(fù)雜，堆石壩體動力分析問題一直是大壩抗震研究領(lǐng)域的難點，而恰恰面板的地震動力反應(yīng)特性很大程度上受堆石壩體的地震動力反應(yīng)特性所控制。近年來，隨著土動力學(xué)和結(jié)構(gòu)動力分析理論、方法的迅速發(fā)展，面板堆石壩（堆石壩體和面板）的地震動力反應(yīng)特性研究也取得了較大進(jìn)步，對于幫助研究人員認(rèn)識面板的地震動力反應(yīng)特性很有意義，例如面板的動應(yīng)變應(yīng)力反應(yīng)過程、極值水平及分布特性和面板抗震措施等[5-6]。然而對于這些理論和分析方法，普遍缺少原型觀測和模型試驗等的實例驗證。對于土石壩的震害調(diào)查，歷來都是大壩抗震研究人員十分重視的工作，例如壩體在實際地震中的地震加速度反應(yīng)[7]、地震殘余變形[8]和大壩局部結(jié)構(gòu)動力破壞性態(tài)[9]等。對面板的震害調(diào)查，則往往局限于震后面板狀態(tài)的觀察、分析[10-11]，對于面板在地震中的動力反應(yīng)過程少見有相關(guān)文獻(xiàn)報道。制作大比尺面板堆石壩模型，在大型振動臺上進(jìn)行地震模擬振動試驗，是研究面板地震動力反應(yīng)特性的另一重要手段。由于模型壩尺寸較大，模型壩（包括面板）的地震動力反應(yīng)測試結(jié)果可靠性高，大型振動臺模型試驗成果對認(rèn)識面板的地震動力反應(yīng)特性十分有益。此外，除利用振動臺模型試驗成果直接分析面板的地震動力反應(yīng)特性外，還可以利用試驗成果對已有面板堆石壩動力分析理論、方法進(jìn)行驗證，對于改進(jìn)這些理論、方法很有幫助。

本文進(jìn)行某實際高面板堆石壩大型振動臺模型試驗，研究大壩混凝土面板的地震動力反應(yīng)特性。研究了面板的三向動應(yīng)變反應(yīng)過程及其峰值在豎直方向上的分布特點，并根據(jù)材料力學(xué)理論推求不同測點的主應(yīng)變和主應(yīng)力反應(yīng)過程，分析其峰值在豎直方向上的分布規(guī)律，再根據(jù)面板堆石壩振動臺模型試驗相似理論，推求了原型壩面板的地震動應(yīng)力反應(yīng)特性。最后，對面板堆石壩面板地震動力反應(yīng)特性大型振動臺模型試驗成果的應(yīng)用，也進(jìn)行了一些有益的討論。

第2期楊正權(quán)，等： 高面板堆石壩面板地震動力反應(yīng)特性大型振動臺模型試驗研究振 動 工 程 學(xué) 報第31卷 1 面板堆石壩大型振動臺模型試驗

某實際水利水電工程的攔河大壩為鋼筋混凝土面板堆石壩，最大壩高156 m，壩頂軸線長635 m，壩頂寬12 m。堆石壩體上游壩坡坡比為1∶1.4；下游坡設(shè)兩級馬道，第2級馬道以上壩坡坡比為1∶1.5，以下壩坡坡比為1∶1.4。壩體內(nèi)部設(shè)墊層區(qū)、反濾區(qū)、過渡區(qū)、主堆石區(qū)和次堆石區(qū)等不同材料分區(qū)。為研究大壩壩體在地震作用下的地震動力反應(yīng)特性和深入分析混凝土面板的地震動力反應(yīng)特性，進(jìn)行了大壩的大型地震模型振動臺模型試驗，本文主要基于大壩模型試驗中面板動力反應(yīng)特性相關(guān)工作展開論述。

1.1 試驗設(shè)備

試驗在三向6自由度大型高性能地震模擬振動臺上進(jìn)行。振動臺臺面尺寸為6 m×6 m，為國內(nèi)臺面尺寸最大的振動臺。振動臺不僅尺寸大，且性能優(yōu)良，振動臺主要性能參數(shù)為：最大載重為80 t；工作頻率為0.1～50 Hz；最大加速度為±1.5g（x），±1.0g（y），±0.8g（z）；最大速度為±1.0 m/s（x），±1.2 m/s（y），±0.8 m/s（z）；最大位移為±15 cm（x），±25 cm（y），±10 cm（z）；最大傾覆力矩為180 t·m。

1.2 模型壩設(shè)計

基于文獻(xiàn)[10]提出的重力場中土石壩振動臺模型試驗相似律進(jìn)行振動臺模型試驗設(shè)計。模型壩壩高1 m，原型壩壩高156 m，所以模型試驗的幾何比尺Cl為156；模型壩和原型壩的控制填筑密度分別為1.93 g/cm3和2.16 g/cm3，即密度相似常數(shù)Cρ為1.12；模型壩填筑料和原型壩堆石料最大動剪模量系數(shù)分別為1459和3159，故動變形模量相似常數(shù)Cc為2.17。通過3個控制相似常數(shù)Cl，Cρ和Cc，以及模型試驗相似律，就可以確定模型試驗其他所有項目的相似常數(shù)，具體如表1所示。

認(rèn)為河谷基巖為剛性體，選擇原型壩13個控制斷面按幾何比尺縮尺作為模型壩的控制斷面，各斷面間用三角形拼接出基巖河谷地形，用鋼筋混凝土來模擬巖石地基。將原型堆石料按混合級配法縮尺后作為模型壩填筑料，模型料最大控制粒徑為20 mm。

模型壩面板用人工配合材料來模擬水泥的作用，用細(xì)軟鐵絲網(wǎng)來模擬鋼筋的作用。根據(jù)模型試驗相似律，計算得到面板動彈模量相似比為28.73（相似律為C1/2ρC1/2lCc）。原型壩面板C30混凝土動彈模量取為32.5 GPa，換算得到模型壩面板的動彈模量為1.13 GPa。根據(jù)模型試驗對面板動彈模量相似的要求，采用人工配合材料制作面板，配合料由水泥、砂子、珍珠巖、水和面粉摻合料等組成，根據(jù)文獻(xiàn)[10]給出的方法來確定各材料的配合比例。

1.3 模型壩制作

為了控制好模型壩的填筑密度及其均勻性，整個壩體按7層稱料，并逐層擊實填筑而成，沿順河和橫河兩個方向控制填筑表面的平整度。堆石壩體填筑完成后，對上游壩坡進(jìn)行嚴(yán)格的整平，再將細(xì)鐵絲網(wǎng)分區(qū)固定于壩坡表面后分區(qū)施工面板，面板施工完成后進(jìn)行2周的晾干養(yǎng)護(hù)，圖1為制作完成的面板堆石壩三維整體模型。表1 模型試驗相似律及相應(yīng)相似常數(shù)

Tab.1 Similarity law of the model test and the corresponding similarity parameters

符號項目相似律〖〗模型試驗相似常數(shù)L壩體幾何尺寸Cl156.0ρ堆石體/面板密度Cρ1.12C堆石體模量系數(shù)Cc = Cp/ Cm2.17G堆石體剪切模量/面板楊氏模量CG= C1/2lCc C1/2ρ28.68σ堆石體/面板應(yīng)力Cσ= C1/2lCc C1/2ρ28.68ε堆石體/面板應(yīng)變Cε= C1/2ρ C1/2l/Cc6.09u堆石體位移Cu= C1/2ρ C3/2l/Cc949.82ù堆石體速度Cù= C1/4ρ C3/4l/C1/2c30.82ü堆石體加速度Cü=Cg=11ξ堆石體阻尼比Cξ=11C′堆石體有效凝聚力CC′ = CρCl174.72φ′堆石體有效摩擦角Cφ′=11t時間Ct= C1/4ρ C3/4l/C1/2c30.82f頻率Cf= C-1/4ρC-3/4lC1/2c0.032

圖1 面板堆石壩振動試驗?zāi)Ｐ蛪?/p>

Fig.1 Model of concrete faced rock-fill dam

地震作用下，混凝土面板的動應(yīng)變反應(yīng)通過粘貼在面板表面的應(yīng)變花進(jìn)行量測，應(yīng)變花由互成45°角的三向應(yīng)變片組成。所有應(yīng)變片為同一批次生產(chǎn)，規(guī)格相同，電阻為120 Ω，絲柵5×3 mm。應(yīng)變片粘貼時，先用砂紙將混凝土表面磨平拋光，再用記號筆在預(yù)設(shè)位置劃線定位，然后用專用膠水將應(yīng)變片粘貼在預(yù)設(shè)位置。在位于河床中央的主控制斷面（最大橫斷面）上沿豎直方向布置6個應(yīng)變花，自壩頂向壩踵布置于前述各填筑層的交界面高程上。應(yīng)變花構(gòu)造及其在測試斷面上的布置情況如圖2（a）所示，三向應(yīng)變片的布置方向分別是x向（橫河水平向），y向（順壩坡方向）和z向（前兩者的等角分線方向），粘貼在面板表面應(yīng)變花的實際狀態(tài)如圖2（b）所示。圖2 應(yīng)變花構(gòu)造及其在測試斷面上的布設(shè)情況（單位： mm）

Fig.2 Strain rosette and its′ layout on control cross-section of model dam （Unit： mm）1.4 輸入地震動

在大壩振動臺模型試驗中，輸入的地震波包括壓縮場地波、規(guī)范波和天然波等，輸入方式包括順河水平單向輸入、水平豎直雙向輸入和三向輸入。本文選取更能反映壩址區(qū)場地地震動特性的場地波作為分析對象，同時為了便于對比分析，使得研究問題盡量簡單化，考慮不同方向地震動輸入對大壩壩體和面板地震動力反應(yīng)影響的差異，分析實例選取大壩結(jié)構(gòu)振動狀態(tài)相對簡單的順河水平單向輸入情況，預(yù)設(shè)輸入加速度峰值分別為0.1g，0.2g和0.3g。圖3和4分別給出了預(yù)設(shè)加速度峰值為0.1g時，臺面實測輸入地震波加速度時程（實測峰值為0.099g）和相應(yīng)的傅里葉頻譜圖。

圖3 臺面實測輸入壓縮場地波加速度時程示例（圖中縱軸加速度單位為g）

Fig.3 Example of actual measured acceleration time-history curves by the shaking table （the unit is ″g″ in figure）

圖4 臺面實測輸入壓縮場地波傅里葉譜示例

Fig.4 Example of actual measured Fourier spectra curves by the shaking table2 面板動應(yīng)變反應(yīng)

將應(yīng)變片實測的應(yīng)變時程曲線經(jīng)過歧義點剔除和漂零處理后，可以得到各測點實際三向動應(yīng)變反應(yīng)時程。模型壩在壓縮場地波作用下，處在大壩中間高程的3號測點（0.6倍壩高）的三向應(yīng)變反應(yīng)時程和相應(yīng)傅里葉譜如圖5和6所示。圖5中縱軸應(yīng)變單位為微應(yīng)變με，1με=10-6；對應(yīng)變方向的符號規(guī)定與材料力學(xué)的一般要求一致，即拉為正壓為負(fù)。對比圖3和4可知，和輸入加速度時程相對應(yīng)，應(yīng)變反應(yīng)比較大的時段出現(xiàn)在0.6 s前后，臺面加速度時程和應(yīng)變反應(yīng)時程在波形上也有很好的對應(yīng)性。從振動頻譜圖也可以看出，面板振動的頻譜特性和地震動輸入有較好的對應(yīng)關(guān)系。測點三向應(yīng)變反應(yīng)

圖5 應(yīng)變花實測測點三向應(yīng)變反應(yīng)時程（3號測點應(yīng)變花）

Fig.5 Dynamic strain response time-history of concrete measured by 3# strain rosette圖6 應(yīng)變花實測測點三向應(yīng)變反應(yīng)傅立葉譜（3號測點應(yīng)變花）

Fig.6 Dynamic strain response Fourier spectra curves of concrete measured by 3# strain rosette

圖7 面板三向應(yīng)變極值（拉壓雙向）沿高程分布的對比

Fig.7 Comparison of the concrete dynamic strain response peak values （three direction） on vertical direction

在拉壓方向上基本呈對稱分布，這表明面板的動應(yīng)變反應(yīng)也基本是在靜力狀態(tài)的基礎(chǔ)上接近對稱的反復(fù)拉壓反應(yīng)。從時程曲線還可以看出，測點的壓應(yīng)變峰值要高于拉應(yīng)變峰值。

對以上應(yīng)變反應(yīng)時稱取雙向（拉和壓）極值，可以得到各向應(yīng)變反應(yīng)極值在豎直方向上的分布情況，如圖7所示。從量值上看，隨著輸入地震動強(qiáng)度的增加，面板的動應(yīng)變反應(yīng)幅度增加；面板三向動應(yīng)變反應(yīng)中，順河向量值最大，45°方向次之，橫河向最小，這是由于振動臺輸入順河水平單向地震波對壩體進(jìn)行激振，壩體在順河方向（和順坡方向基本一致）上的振動幅度最大，面板的動力反應(yīng)在此方向上也就最強(qiáng)烈。從分布特性上看，橫河向應(yīng)變反應(yīng)的極值點出現(xiàn)在壩體頂端，且沿壩踵向壩頂是逐漸增大的，面板順坡向和45°方向上的應(yīng)變反應(yīng)最大值均出現(xiàn)在靠近壩體中部的0.6倍壩高高程上，壩踵和壩頂處面板應(yīng)變反應(yīng)的幅度最?。回Q直方向上，三向應(yīng)變拉壓雙向分布曲線形狀基本對稱，但是壓應(yīng)變水平要明顯高于拉應(yīng)變水平。應(yīng)變反應(yīng)極值在豎直方向上的分布特性，主要是受面板在地震過程中的加速度反應(yīng)和變形過程所控制。

根據(jù)材料力學(xué)理論，線彈性混凝土材料單元的大小主應(yīng)變（ε1，ε3）同雙向正交（x，y方向）及其45°夾角方向（z方向）應(yīng)變間的關(guān)系可以用下式來表示γxy=2ε45°-εx-εy

ε1=12[（εx+εy）+（εx-εy）2+γ2xy]

ε3=12[（εx+εy）-（εx-εy）2+γ2xy]（1）式中 γxy為最大剪應(yīng)變。

對前述測點的三向應(yīng)變時程按公式（1）可以計算得到測點的大小主應(yīng)變時程。圖8所示為模型壩在壓縮場地波作用下（臺面輸入加速度峰值0.099g），3號測點處面板的大小主應(yīng)變反應(yīng)時程，其縱軸應(yīng)變單位同圖5。從測點大小主應(yīng)變時程曲線的形態(tài)可以看出，大主應(yīng)變在大部分時段內(nèi)為拉應(yīng)變，而小主應(yīng)變在大部分時段為壓應(yīng)變，大小主應(yīng)變曲線沿時間軸的分布形態(tài)基本對稱，小主應(yīng)變絕對值高于大主應(yīng)變。

圖8 計算求得測點大小主應(yīng)變反應(yīng)時程（3號測點應(yīng)變花）

Fig.8 Calculated principal strain response time-history of concrete （3# strain rosette）

圖9給出了大小主應(yīng)變雙向（拉壓）極值沿豎直方向上分布的對比。如圖9所示，大主應(yīng)變在受拉方向上的極值和小主應(yīng)變在受壓方向上的極值在豎直方向上分布曲線的形態(tài)基本一致，而大主應(yīng)變在受壓方向上的極值和小主應(yīng)變在受拉方向上的極值在豎直方向上分布曲線的形態(tài)也基本一致。面板所受壓應(yīng)變極值高于所受拉應(yīng)變極值，拉壓應(yīng)變的極值點均出現(xiàn)在靠近壩體中部0.6倍壩高的高程上。

圖9 面板主應(yīng)變極值（拉壓雙向）沿豎直方向分布對比

Fig.9 Comparison of the concrete principal strain response peak values on vertical direction3 面板動應(yīng)力反應(yīng)

基于材料力學(xué)理論，可根據(jù)單元體任意正交兩向應(yīng)變和最大剪應(yīng)變計算求得單元體的大小主應(yīng)力（σ1，σ3）和最大剪應(yīng)力（τmax），如下式所示σ1=12E1-μ（εx+εy）+12E1+μ（εx-εy）2+γ2xy

σ3=12E1-μ（εx+εy）-12E1+μ（εx-εy）2+γ2xy

τmax=E2（1+μ）（εx-εy）2+γ2xy（2） 如前述，模型壩面板混凝土的動彈模量E取為1.13 GPa；根據(jù)經(jīng)驗公式，確定動泊松比ν為0.16[10]，式中其他符號意義同前。

基于測點三向應(yīng)變時程，根據(jù)公式（2）可以計算得到測點的大小主應(yīng)力和最大剪應(yīng)力時程。圖10給出了模型壩在壓縮場地波作用下（臺面輸入加速度峰值0.099g），3號測點處混凝土大小主應(yīng)力和最大剪應(yīng)力反應(yīng)時程，圖中縱軸應(yīng)力單位為Pa。由圖10可知，面板所受動態(tài)大小主應(yīng)力時程曲線的相互關(guān)系和大小主應(yīng)變的相互關(guān)系一致：大主應(yīng)力大部分時段內(nèi)為拉應(yīng)力，而小主應(yīng)力大部分時段為壓應(yīng)

圖10 計算求得測點大小主應(yīng)力和最大剪應(yīng)力反應(yīng)時程（3號測點應(yīng)變花）

Fig.10 Calculated principal stress response time-history of concrete （3# strain rosette）力，大小主應(yīng)力曲線沿時間軸的分布形態(tài)基本對稱，但小主應(yīng)力的絕對值極值高于大主應(yīng)力。

對面板大小主應(yīng)力反應(yīng)時程雙向（拉和壓）取極值，最大剪應(yīng)力時程取極值，得到大小主應(yīng)力和最大剪應(yīng)力極值沿豎直方向分布的對比，如圖11所示。由圖11可知，和大小主應(yīng)變沿高程分布規(guī)律一致，大主應(yīng)力在受拉方向上的極值和小主應(yīng)力在受壓方向上的極值在豎直方向上的分布曲線形態(tài)基本一致，而大主應(yīng)力在受壓方向上的極值和小主應(yīng)力在受拉方向上的極值在豎直方向上的分布曲線形態(tài)也基本一致。面板所受動壓應(yīng)力極值高于所受動拉應(yīng)力，拉壓應(yīng)力的極值點均出現(xiàn)在靠近壩體中部0.6倍壩高左右的高程上。0.6～0.8倍壩高高程處混凝土所受動剪應(yīng)力水平最高。面板在地震過程中的動拉壓應(yīng)變反應(yīng)和動拉壓應(yīng)力反應(yīng)極值均出現(xiàn)在大壩

圖11 面板大小主應(yīng)力（拉壓雙向）和最大剪應(yīng)力極值沿高程分布對比

Fig.11 Comparison of the concrete principal stress response peak values on vertical direction0.6左右高程上的這一分布特點，和面板及壩體在地震中的動力反應(yīng)特性是相關(guān)聯(lián)的：橫斷面上，面板是一個長高比較大的梁結(jié)構(gòu)，在靜動力作用下，中間段靜動力變形（撓度）大，所引起的應(yīng)力應(yīng)變反應(yīng)自然較大；而且，面板坐落于壩坡表面，其動力反應(yīng)和壩體在地震中的動力反應(yīng)，尤其是壩坡表面的動力反應(yīng)特性也是相關(guān)聯(lián)的。

此外，試驗結(jié)束后對面板狀態(tài)進(jìn)行了仔細(xì)觀察，未發(fā)現(xiàn)有明顯的張拉裂縫和擠壓破壞，面板和周圍山體也沒有出現(xiàn)明顯的裂縫和錯位。在最大幅值為0.3g的場地地震波作用下，面板未出現(xiàn)會影響其安全性的明顯破壞。

4 面板動力反應(yīng)振動臺模型試驗成果應(yīng)用探討 根據(jù)模型試驗相似律，原型壩和模型壩間面板應(yīng)力的相似比為174.7（CρCl）。針對面板0.6倍壩高高程的3號測點，推求原型壩相應(yīng)位置處的動應(yīng)力反應(yīng)極值：在加速度峰值為0.099g，0.208g和0.286g順河水平單向場地波作用下，面板所受動壓應(yīng)力峰值分別為1.10，2.54和3.22 MPa，動拉應(yīng)力峰值分別為0.91，1.94和2.60 MPa。對比基于三維真非線性有限元分析方法的計算成果可知[10]，在量值水平相當(dāng)?shù)耐N地震波作用下，模型試驗推求原型壩面板動拉壓應(yīng)力反應(yīng)極值與實際大壩動力分析得到面板動拉壓應(yīng)力極值量值水平基本一致：主斷面0.6倍壩高處，輸入地震動峰值加速度為0.260g時，計算求得面板的動壓應(yīng)力極值和動拉應(yīng)力極值分別為2.6 MPa左右和2.5 MPa左右，模型試驗確定的面板最大拉應(yīng)力水平和計算求得的數(shù)值較接近，最大壓應(yīng)力則高于計算求得數(shù)值。同時，參考計算求得原型壩主斷面上動拉壓應(yīng)力極值沿高程分布情況可知，模型試驗和計算確定面板應(yīng)力沿高程分布一般特性也是一致的。計算確定面板拉壓應(yīng)力的極值區(qū)位于0.7倍壩高左右的高程上，但是模型試驗在此高程上沒有布置應(yīng)變片測點，前述模型試驗確定的面板拉壓應(yīng)力極值區(qū)高程略低于計算結(jié)果。

由于利用大型振動臺模型試驗進(jìn)行面板動力反應(yīng)特性研究的研究成果還較少，因此可以通過對比模型試驗成果和理論計算分析以及原型震害調(diào)查研究成果，來考察用振動臺試驗方法研究面板動力反應(yīng)特性的可行性。當(dāng)然，模型材料的材料動力特性指標(biāo)為根據(jù)經(jīng)驗公式推算而得，非試驗成果，所取動彈模量值可能和實際值有所差異，還不能從以上模型試驗和計算分析的對比關(guān)系，對模型試驗和計算研究成果的準(zhǔn)確性予以嚴(yán)格的驗證。但是，從兩者的對應(yīng)關(guān)系看，模型試驗成果和基于這些成果推求得到的原型壩相應(yīng)特性與基于一般力學(xué)理論的計算分析成果，不管是從定性規(guī)律、分布特性上看，還是從量值水平上看，兩者（模型試驗和理論分析）都有較好的對應(yīng)關(guān)系，兩者之間相互驗證，表明模型試驗成果是值得信賴的。

對于重力場中土工結(jié)構(gòu)振動臺模型試驗的相似律問題，目前仍舊處在探索階段，對其成果的適用性問題仍需進(jìn)一步的討論和驗證。由于土體動力變形特性受土體應(yīng)力狀態(tài)水平的影響，和可以近似看作是線彈性體的混凝土結(jié)構(gòu)相比，土工結(jié)構(gòu)的動力模型試驗相似律相對復(fù)雜，目前土工結(jié)構(gòu)動力模型試驗研究成果的應(yīng)用應(yīng)該還處在總體定性、局部定量的階段。但是，盡管模型應(yīng)力場和原型應(yīng)力場有差距，重力場中大型振動臺模型試驗也有其他動力模型試驗無法比擬的優(yōu)勢，由于模型相對原型結(jié)構(gòu)可以做到較大的比尺，對原型結(jié)構(gòu)和材料的模擬更加接近實際，試驗過程中的測試結(jié)果也更加可靠。而且，從目前土石壩大型振動臺模型試驗已取得的成果看，振動臺模型試驗在土石壩結(jié)構(gòu)動力特性測試，加速度反應(yīng)測試和地震動力反應(yīng)過程及破壞機(jī)理分析等方面，都和實際震害調(diào)查及理論分析成果有較好的對應(yīng)關(guān)系。這些試驗研究成果對于幫助工程技術(shù)人員認(rèn)識相關(guān)規(guī)律以及對大壩動力分析相關(guān)理論、方法的驗證和改進(jìn)，定性指導(dǎo)大壩抗震設(shè)計，還是很有益處的。從本文的研究成果也可以看到，試驗得到面板動力反應(yīng)沿大壩高程方向上的分布規(guī)律和量值水平，都是和理論分析成果可以相互驗證的，是可信的?？梢曰谡駝优_模型試驗方法，繼續(xù)深入進(jìn)行面板地震動力反應(yīng)相關(guān)影響規(guī)律的研究，并指導(dǎo)實際工程設(shè)計。

此外，受模型制作和試驗測試自身情況的限制，本文只在主斷面上沿壩坡方向布置了應(yīng)變花測點，測試得到了沿壩高方向面板的應(yīng)力應(yīng)變反應(yīng)分布特性。對于沿壩軸線方向，以及面板和河谷基巖連接的周邊區(qū)域的面板動應(yīng)力應(yīng)變反應(yīng)，沒有進(jìn)行相關(guān)測試，在今后的研究中有待補(bǔ)充完善。而且，受篇幅所限，本文只以單向輸入場地波為例，對振動臺模型試驗測試面板動力反應(yīng)特性的分析過程、基本研究成果和試驗成果的可靠性及適用性等進(jìn)行了分析，關(guān)于地震動類型、輸入方式和水庫蓄水狀態(tài)等相關(guān)因素的影響分析，以及面板在強(qiáng)震作用下的破壞模式和破壞規(guī)律等的研究成果，將另文予以介紹。

5 結(jié) 論

通過高面板堆石壩大型地震模擬振動臺模型試驗，研究了面板在地震作用下的各向動應(yīng)變反應(yīng)過程和空間分布規(guī)律，根據(jù)試驗測試得到面板三向動應(yīng)變反應(yīng)過程和材料力學(xué)理論，推求了面板的大小主應(yīng)變反應(yīng)過程及其極值分布規(guī)律，以及面板的動應(yīng)力反應(yīng)過程及其極值分布規(guī)律，并嘗試根據(jù)模型試驗相似律，利用振動臺模型試驗成果推求原型壩混凝土面板的動力反應(yīng)特性。主要研究結(jié)論為：

（1）振動臺模型試驗實測面板三向應(yīng)變反應(yīng)過程曲線及其推求主應(yīng)變和主應(yīng)力過程曲線符合對面板動力反應(yīng)過程的一般認(rèn)識，振動臺模型試驗可作為研究面板地震動力反應(yīng)特性的重要手段；

（2）在順河水平單向場地地震波作用下，沿豎直方向上面板的大小主應(yīng)變極值均出現(xiàn)在靠近壩體中部0.6倍壩高左右的高程上，不同高程處面板在地震過程中動拉應(yīng)變總體小于動壓應(yīng)變；

（3）和應(yīng)變分布規(guī)律類似，面板大小主應(yīng)力和最大剪應(yīng)力極值同樣分布在靠近壩體中部0.6倍壩高左右的高程上，面板所受動拉應(yīng)力低于動壓應(yīng)力；

（4）基于面板堆石壩振動臺模型試驗相似理論，根據(jù)模型試驗成果，推求得到原型大壩面板應(yīng)力應(yīng)變反應(yīng)的量值水平和分布規(guī)律，和基于一般力學(xué)理論的數(shù)值分析成果基本一致，兩者可以在一定程度上相互驗證。

參考文獻(xiàn)：

[1] 楊正權(quán)，劉小生，汪小剛，等. 高土石壩地震動力反應(yīng)特性大型振動臺模型試驗研究[J]. 水利學(xué)報， 2014，45（11）：1361—1372.

Yang Zhengquan， Liu Xiaosheng， Wang Xiaogang， et al. Study on dynamic response characteristics of high rock-fill dam in earthquake by large-scale shaking table model tests[J]. Journal of Hydraulic Engineering， 2014，45（11）：1361—1372.

[2] 楊正權(quán)，劉小生，劉啟旺，等. 高面板堆石壩地震反應(yīng)和破壞振動臺模型試驗研究[J].振動工程學(xué)報， 2015，28（6）：937—945.

Yang Zhengquan， Liu Xiaosheng， Liu Qiwang， et al. Study on shaking table model tests for dynamic response and failure of high concrete face rock-fill dam[J]. Journal of Vibration Engineering， 2015，28（6）：937—945

[3] 鄧 剛，汪小剛，溫彥鋒，等. 混凝土面板壩面板變形模式與水平向擠壓破損研究[J]. 水利學(xué)報， 2015，46（4）：396—404.

Deng Gang， Wang Xiaogang， Wen Yanfeng， et al. Study on conceptualization method of deformation pattern and horizontal breakage of face slab of concrete faced rock-fill dam[J]. Journal of Hydraulic Engineering， 2015，46（4）：396—404

[4] 楊正權(quán)，趙劍明，劉小生，等. 紫坪鋪大壩下游壩坡震后抗震加固措施大型振動臺模型試驗研究[J]. 巖土工程學(xué)報， 2015， 37（11）： 2058—2066.

Yang Zhengquan， Zhao Jianming， Liu Xiaosheng， et al. Study on large-scale shaking model test for dam slope anti-seismic and reinforcing measures of Zipingpu earth-rock fill dam after “5.12” Wenchuan Earthquake [J]. Chinese Journal of Geotechnical Engineering， 2015， 37（11）： 2058—2066.

[5] 孔憲京，鄒德高，鄧學(xué)晶，等. 高土石壩綜合抗震措施及其效果的驗算[J]. 水利學(xué)報， 2006，37（12）： 1489—1495.

Kong Xianjing， Zou Degao， Deng Xuejing， et al. Comprehensive earthquake resistant measure of high earth-rock fill dams and effectiveness verification [J]. Journal of Hydraulic Engineering， 2006， 37（12）： 1489—1495.

[6] 孔憲京，張 宇，鄒德高. 超高面板堆石壩面板地震應(yīng)力改善措施研究[J]. 水利學(xué)報， 2014，45（4）：419—426.

Kong Xianjing， Zhang Yu， Zou Degao. Study on aseismic measures of high concrete-face rock-fill dam[J]. Journal of Hydraulic Engineering， 2014，45（4）：419—426.

[7] 劉振平，遲世春. 鯉魚潭大壩在集集地震中的變形分析[J]. 巖土力學(xué)， 2014，35（1）：248—254.

Liu Zhenping， Chi Shichun. Analysis of deformation of Liyutan dam after Chi-chi earthquake [J]. Rock and Soil Mechanics， 2014， 35（1）： 248—254.

[8] 陳生水，霍家平，章為民. “5.12”汶川地震對紫坪鋪混凝土面板壩的影響及原因分析[J]. 巖土工程學(xué)報， 2008，30（6）：795—801.

Chen Sheng-shui， Huo Jia-ping， Zhang Wei-min. Analysis of effects of “5.12” Wenchuan earthquake on Zipingpu concrete face rock-fill dam[J]. Chinese Journal of Geotechnical Engineering， 2008，30（6）：795—801.

[9] 趙劍明，劉小生，溫彥鋒，等. 紫坪鋪大壩汶川地震震害分析及高土石壩抗震減災(zāi)研究設(shè)想[J]. 水力發(fā)電， 2009， 35（5）： 11—14.

Zhao Jianming， Liu Xiaosheng， Wen Yanfeng， et al. Analysis of earthquake damage of the Zipingpu dam in Wenchuan earthquake and the study proposal on the anti-earthquake and disaster reduction of high earth-rock dam [J]. Water Power， 2009， 35（5）： 11—14.

[10]劉小生，王鐘寧，汪小剛，等. 面板壩大型振動臺模型試驗與動力分析[M].北京：中國水利水電出版社，2005.

Liu Xiaosheng， Wang Zhongning， Wang Xiaogang， et al. Large-scale Shaking Model Tests and Dynamic Analysis of Concrete Face Rock-fill Dam [M]. Beijing： China Water Resource and Hydro-powder Publishing House， 2005.

[11]Guan Zhicheng. Investigation of the 5.12 Wenchuan earthquake damages to the Zipingpu water control project and an assessment of it′s safety state[J]. Science in China Series E： Technological Sciences， 2009，52（4）：820—834.

Study on large-scale shaking table model tests for face slab′s dynamic response

characteristics of high concrete faced rock-fill dam

YANG Zheng-quan1，2， LIU Xiao-sheng1，2， ZHOU Guo-bin3， LIU Qi-wang1，2，

ZHAO Jian-ming1，2， YANG Yu-sheng1，2

（1.State Key Laboratory of Simulation and Regulation of Water Cycle in River Basin， Beijing 100038， China；

2.China Institute of Water Resources and Hydropower Research， Beijing 100048， China；

3.Guangdong Hydropower Planning & Design Institute， Guangzhou 510635， China）

Abstract： The concrete face slab is the key structure in the anti-seepage system of the concrete faced rock-fill dam， and its safety is the basic guarantee for normal running of the dam. Structural safety of the face slab may be affected greatly by the superposition of dynamic response and static response in the strong earthquake， and the dynamic response characteristics of the face slab has always been the focus in the areas of dam anti-seismic research. In this paper， large-scale shaking table model tests of the high concrete faced rock-fill dam are carried out to study the dynamic response characteristics of the face slab. The dynamic strain response process of the dam face slab and the spatial distribution law of peak values of strains are analyzed firstly. Then， the dynamic stress response process of the dam face slab and the spatial distribution law of peak values of stresses are studied based on the dynamic strain response process and the material mechanics theory.Finally， the dynamic response characteristics of the prototype dam's concrete face slab are investigated based on the results of shaking model tests and the similarity law. The results show that the measured dynamic strain response process of the dam face slab in the shaking model tests and the corresponding calculated dynamic principal stress & strain response process are coincident with the general recognition to the dynamic response characteristics of the dam face slab， and the dam shaking model test can be used as an important method to research dynamic response characteristics of the dam face slab. The peak values of the principal strain appear on the altitude of about 0.6 dam height on dam vertical direction， and the peak values of the principal stress also appear on the same height correspondingly. The dynamic tensile stress undertaken by the concrete face slab in earthquake is less than the dynamic compressive stress generally.

Key words： earthquake respanse； high concrete faced rock-fill dam； concrete face slab； large-scale shaking model test