• <tr id="yyy80"></tr>
  • <sup id="yyy80"></sup>
  • <tfoot id="yyy80"><noscript id="yyy80"></noscript></tfoot>
  • 99热精品在线国产_美女午夜性视频免费_国产精品国产高清国产av_av欧美777_自拍偷自拍亚洲精品老妇_亚洲熟女精品中文字幕_www日本黄色视频网_国产精品野战在线观看 ?

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

    2018-06-14 14:54楊正權(quán)劉小生周國斌劉啟旺趙劍明楊玉生
    振動工程學(xué)報 2018年2期
    關(guān)鍵詞:模型試驗測點極值

    楊正權(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蛪?/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

    猜你喜歡
    模型試驗測點極值
    通過函數(shù)構(gòu)造解決極值點偏移問題
    例談解答極值點偏移問題的方法
    極值點偏移問題的解法
    輸油泵的性能測試探索
    鄭州龍湖調(diào)蓄工程1#出口控制閘閘門及啟閉機(jī)布置研究
    基于監(jiān)測的空間網(wǎng)格結(jié)構(gòu)應(yīng)力應(yīng)變分析
    低路堤在車輛荷載作用下響應(yīng)的模型試驗
    阿爾塔什水利樞紐水墊塘消能方式選擇
    也談?wù)剺O值點偏移問題
    淺談某水電站引水隧洞的圍巖變形量測研究
    69精品国产乱码久久久| 99国产极品粉嫩在线观看| 欧美 亚洲 国产 日韩一| 免费在线观看日本一区| 色视频在线一区二区三区| 国精品久久久久久国模美| 50天的宝宝边吃奶边哭怎么回事| 色视频在线一区二区三区| 亚洲 欧美一区二区三区| 亚洲 欧美一区二区三区| 在线观看舔阴道视频| 日日爽夜夜爽网站| 午夜日韩欧美国产| netflix在线观看网站| 久久久精品区二区三区| 欧美变态另类bdsm刘玥| 亚洲国产精品一区二区三区在线| 欧美成人免费av一区二区三区 | 日本av手机在线免费观看| 精品福利观看| 一个人免费看片子| 少妇裸体淫交视频免费看高清 | 亚洲第一青青草原| 老鸭窝网址在线观看| 最黄视频免费看| 久久影院123| 国产精品一区二区精品视频观看| 一区二区日韩欧美中文字幕| 久久亚洲真实| 在线av久久热| 国产精品久久久久久精品电影小说| 国产男靠女视频免费网站| 在线天堂中文资源库| 18在线观看网站| 欧美日韩亚洲国产一区二区在线观看 | 亚洲黑人精品在线| 久久久久视频综合| 亚洲午夜精品一区,二区,三区| 97在线人人人人妻| 中文字幕人妻丝袜制服| 桃花免费在线播放| 又黄又粗又硬又大视频| 精品人妻1区二区| 每晚都被弄得嗷嗷叫到高潮| 变态另类成人亚洲欧美熟女 | 自线自在国产av| 亚洲成av片中文字幕在线观看| 亚洲一区二区三区欧美精品| 亚洲avbb在线观看| 久久久久网色| 国产av一区二区精品久久| 男女免费视频国产| 多毛熟女@视频| 五月天丁香电影| 在线观看人妻少妇| 中文字幕人妻丝袜制服| 国产男女内射视频| 国产精品熟女久久久久浪| videos熟女内射| 午夜福利在线免费观看网站| 99久久精品国产亚洲精品| 久久 成人 亚洲| 蜜桃国产av成人99| 亚洲精品美女久久av网站| 国产av又大| 欧美久久黑人一区二区| 欧美 日韩 精品 国产| 深夜精品福利| 国产精品美女特级片免费视频播放器 | 国产男靠女视频免费网站| 性高湖久久久久久久久免费观看| 99久久精品国产亚洲精品| netflix在线观看网站| 丝瓜视频免费看黄片| 美女高潮到喷水免费观看| 亚洲天堂av无毛| 人人澡人人妻人| 黄网站色视频无遮挡免费观看| 80岁老熟妇乱子伦牲交| 国产色视频综合| 50天的宝宝边吃奶边哭怎么回事| 亚洲av成人不卡在线观看播放网| 中文字幕人妻丝袜制服| 免费av中文字幕在线| 国产精品98久久久久久宅男小说| 深夜精品福利| 久久久精品国产亚洲av高清涩受| 大片免费播放器 马上看| 成人国产一区最新在线观看| 欧美激情高清一区二区三区| 国产精品久久电影中文字幕 | 91精品国产国语对白视频| 色尼玛亚洲综合影院| 人人妻人人澡人人看| 亚洲天堂av无毛| 女警被强在线播放| 亚洲国产中文字幕在线视频| 人人妻人人添人人爽欧美一区卜| 亚洲国产欧美网| 岛国在线观看网站| 国产一区二区三区视频了| 自线自在国产av| 久久精品aⅴ一区二区三区四区| 亚洲精品美女久久av网站| 国产精品二区激情视频| 国产精品二区激情视频| 国产免费av片在线观看野外av| 日韩精品免费视频一区二区三区| 另类亚洲欧美激情| 香蕉国产在线看| 高潮久久久久久久久久久不卡| 不卡一级毛片| 亚洲性夜色夜夜综合| 女人精品久久久久毛片| 亚洲欧美日韩另类电影网站| 国产欧美亚洲国产| 色婷婷av一区二区三区视频| 国产精品98久久久久久宅男小说| 日韩三级视频一区二区三区| 最新在线观看一区二区三区| 变态另类成人亚洲欧美熟女 | 亚洲精品久久午夜乱码| 色在线成人网| 狠狠狠狠99中文字幕| 亚洲avbb在线观看| 免费在线观看影片大全网站| 国产免费现黄频在线看| 亚洲国产欧美在线一区| 另类精品久久| 国产精品九九99| 久久久精品国产亚洲av高清涩受| 老熟妇乱子伦视频在线观看| 深夜精品福利| 热99re8久久精品国产| 午夜91福利影院| 国产野战对白在线观看| 交换朋友夫妻互换小说| 成人18禁高潮啪啪吃奶动态图| 老汉色∧v一级毛片| 女人久久www免费人成看片| 国产成人一区二区三区免费视频网站| 国产精品国产高清国产av | 久久天躁狠狠躁夜夜2o2o| 亚洲精品国产色婷婷电影| 超碰97精品在线观看| 国产午夜精品久久久久久| 99热网站在线观看| 女同久久另类99精品国产91| 国产三级黄色录像| 午夜激情久久久久久久| 亚洲精品成人av观看孕妇| 人人妻人人爽人人添夜夜欢视频| 国产97色在线日韩免费| 精品少妇一区二区三区视频日本电影| 国产亚洲欧美在线一区二区| 欧美日韩亚洲国产一区二区在线观看 | 久久国产亚洲av麻豆专区| 新久久久久国产一级毛片| 在线看a的网站| 操出白浆在线播放| 免费观看a级毛片全部| cao死你这个sao货| 超色免费av| www.熟女人妻精品国产| 欧美日韩亚洲高清精品| 精品福利永久在线观看| 黄色a级毛片大全视频| 999久久久国产精品视频| 99国产综合亚洲精品| 91麻豆av在线| 国产高清视频在线播放一区| 满18在线观看网站| 国产精品电影一区二区三区 | 亚洲av电影在线进入| 成年女人毛片免费观看观看9 | av视频免费观看在线观看| 国产欧美日韩一区二区三区在线| 亚洲 欧美一区二区三区| 精品卡一卡二卡四卡免费| 午夜日韩欧美国产| 日本黄色视频三级网站网址 | 亚洲专区国产一区二区| 少妇精品久久久久久久| 亚洲久久久国产精品| 久久精品成人免费网站| 亚洲第一欧美日韩一区二区三区 | 亚洲精品在线观看二区| 捣出白浆h1v1| 亚洲av国产av综合av卡| 午夜久久久在线观看| 老司机深夜福利视频在线观看| 久久精品国产a三级三级三级| 丝袜美腿诱惑在线| 夜夜夜夜夜久久久久| 怎么达到女性高潮| 国产成人系列免费观看| 最黄视频免费看| 亚洲va日本ⅴa欧美va伊人久久| 亚洲,欧美精品.| 视频区图区小说| 欧美激情久久久久久爽电影 | 色播在线永久视频| 国产一卡二卡三卡精品| 久久久久国产一级毛片高清牌| 麻豆乱淫一区二区| 国产亚洲精品久久久久5区| av天堂久久9| 黄色 视频免费看| 男女无遮挡免费网站观看| 日本av手机在线免费观看| 一区二区三区乱码不卡18| 女人被躁到高潮嗷嗷叫费观| 一区福利在线观看| 国产精品自产拍在线观看55亚洲 | 精品福利永久在线观看| 国产成人啪精品午夜网站| 色尼玛亚洲综合影院| 岛国在线观看网站| 免费看a级黄色片| 极品教师在线免费播放| 久久这里只有精品19| 精品少妇内射三级| 嫩草影视91久久| 国产男女超爽视频在线观看| 一区在线观看完整版| 9色porny在线观看| 亚洲精品一卡2卡三卡4卡5卡| 亚洲美女黄片视频| 国产免费现黄频在线看| 亚洲人成电影观看| 操美女的视频在线观看| 亚洲欧美一区二区三区久久| 大型黄色视频在线免费观看| 国产欧美亚洲国产| 久久久国产成人免费| 人人澡人人妻人| 不卡一级毛片| 十八禁网站网址无遮挡| 国产精品偷伦视频观看了| 久久人妻av系列| 建设人人有责人人尽责人人享有的| 日韩欧美一区二区三区在线观看 | 国产欧美日韩精品亚洲av| 久久性视频一级片| 淫妇啪啪啪对白视频| 亚洲中文字幕日韩| 日本精品一区二区三区蜜桃| 国产欧美日韩一区二区精品| 九色亚洲精品在线播放| 90打野战视频偷拍视频| 怎么达到女性高潮| 精品视频人人做人人爽| e午夜精品久久久久久久| 欧美成人午夜精品| av网站免费在线观看视频| 久久 成人 亚洲| av视频免费观看在线观看| aaaaa片日本免费| 国产亚洲av高清不卡| 久久人妻福利社区极品人妻图片| 午夜福利视频精品| 一本一本久久a久久精品综合妖精| av又黄又爽大尺度在线免费看| 亚洲国产欧美日韩在线播放| 99riav亚洲国产免费| 我要看黄色一级片免费的| 91精品三级在线观看| 人妻久久中文字幕网| 纯流量卡能插随身wifi吗| av超薄肉色丝袜交足视频| 国产成人精品久久二区二区免费| 黄片小视频在线播放| 777米奇影视久久| 18禁观看日本| 亚洲人成电影观看| 九色亚洲精品在线播放| 国产成人精品无人区| 99国产精品99久久久久| 久久精品熟女亚洲av麻豆精品| 男女边摸边吃奶| 熟女少妇亚洲综合色aaa.| 欧美人与性动交α欧美软件| av天堂在线播放| 90打野战视频偷拍视频| 精品一区二区三区四区五区乱码| 欧美乱码精品一区二区三区| 十八禁高潮呻吟视频| 少妇裸体淫交视频免费看高清 | 亚洲人成伊人成综合网2020| av国产精品久久久久影院| 好男人电影高清在线观看| 久久人妻福利社区极品人妻图片| 久久久国产一区二区| 一二三四在线观看免费中文在| 欧美黄色片欧美黄色片| 国产成人啪精品午夜网站| 涩涩av久久男人的天堂| 国产成人一区二区三区免费视频网站| 每晚都被弄得嗷嗷叫到高潮| 汤姆久久久久久久影院中文字幕| 亚洲专区国产一区二区| 久久精品国产a三级三级三级| 日韩欧美免费精品| 91麻豆av在线| 亚洲专区国产一区二区| 亚洲男人天堂网一区| 午夜精品久久久久久毛片777| 亚洲中文日韩欧美视频| 免费看a级黄色片| 国产av精品麻豆| 一级片'在线观看视频| 一二三四在线观看免费中文在| 亚洲男人天堂网一区| 在线观看免费午夜福利视频| 欧美精品啪啪一区二区三区| 国产av精品麻豆| 窝窝影院91人妻| 亚洲,欧美精品.| √禁漫天堂资源中文www| 91成人精品电影| 在线观看免费午夜福利视频| 久久久久久人人人人人| 麻豆国产av国片精品| 一边摸一边做爽爽视频免费| 亚洲精华国产精华精| 欧美国产精品va在线观看不卡| 90打野战视频偷拍视频| 色综合婷婷激情| 在线天堂中文资源库| 国产在线精品亚洲第一网站| 宅男免费午夜| 久久人妻av系列| av视频免费观看在线观看| 成人三级做爰电影| 国产精品香港三级国产av潘金莲| 亚洲精品国产一区二区精华液| 欧美av亚洲av综合av国产av| 国产97色在线日韩免费| 国产男靠女视频免费网站| 国产极品粉嫩免费观看在线| 国产亚洲av高清不卡| 国产av精品麻豆| 国产男女内射视频| 嫁个100分男人电影在线观看| 大陆偷拍与自拍| 夜夜夜夜夜久久久久| 一本久久精品| 搡老岳熟女国产| 一级片免费观看大全| 久久人人爽av亚洲精品天堂| 亚洲国产欧美日韩在线播放| 久久久水蜜桃国产精品网| 国产av又大| 中文字幕精品免费在线观看视频| 他把我摸到了高潮在线观看 | 国产精品久久久久成人av| 久久久欧美国产精品| 久久午夜亚洲精品久久| 亚洲av成人一区二区三| av视频免费观看在线观看| 午夜老司机福利片| 中文字幕精品免费在线观看视频| 精品国产一区二区三区久久久樱花| 免费一级毛片在线播放高清视频 | netflix在线观看网站| 19禁男女啪啪无遮挡网站| 久久青草综合色| 深夜精品福利| 国产av国产精品国产| 亚洲熟妇熟女久久| 精品国内亚洲2022精品成人 | 国产亚洲精品久久久久5区| 午夜福利影视在线免费观看| 亚洲av第一区精品v没综合| 日韩免费高清中文字幕av| 99国产综合亚洲精品| 欧美乱码精品一区二区三区| 亚洲性夜色夜夜综合| 亚洲欧美一区二区三区久久| 性色av乱码一区二区三区2| 中文字幕精品免费在线观看视频| 国产福利在线免费观看视频| 黄色视频在线播放观看不卡| 一级,二级,三级黄色视频| 一二三四在线观看免费中文在| 男女边摸边吃奶| 亚洲精品成人av观看孕妇| 久久久久国内视频| 不卡一级毛片| 视频区图区小说| 老熟女久久久| 少妇的丰满在线观看| 搡老岳熟女国产| 大片电影免费在线观看免费| 国产高清激情床上av| 久久久精品94久久精品| 免费少妇av软件| 日韩大码丰满熟妇| 黄色成人免费大全| 91国产中文字幕| 亚洲精品中文字幕在线视频| 国产国语露脸激情在线看| 每晚都被弄得嗷嗷叫到高潮| 国产麻豆69| 国产日韩欧美亚洲二区| 久久久久国产一级毛片高清牌| 巨乳人妻的诱惑在线观看| 2018国产大陆天天弄谢| avwww免费| 在线播放国产精品三级| 午夜福利,免费看| 黄色a级毛片大全视频| 亚洲精品美女久久av网站| 精品国内亚洲2022精品成人 | 国产成人精品无人区| 精品福利观看| 91国产中文字幕| 12—13女人毛片做爰片一| 人人妻人人澡人人看| 十八禁高潮呻吟视频| 可以免费在线观看a视频的电影网站| 在线亚洲精品国产二区图片欧美| 黄片大片在线免费观看| 热re99久久精品国产66热6| 波多野结衣一区麻豆| 在线观看www视频免费| 亚洲色图 男人天堂 中文字幕| 国产无遮挡羞羞视频在线观看| 男女之事视频高清在线观看| 精品久久久精品久久久| 亚洲精品中文字幕一二三四区 | 国产av精品麻豆| 少妇猛男粗大的猛烈进出视频| 黑丝袜美女国产一区| 两性午夜刺激爽爽歪歪视频在线观看 | 免费在线观看影片大全网站| 精品一区二区三区四区五区乱码| 亚洲第一av免费看| 成人影院久久| 岛国毛片在线播放| 欧美乱码精品一区二区三区| 交换朋友夫妻互换小说| 日本黄色视频三级网站网址 | 啦啦啦在线免费观看视频4| 一个人免费在线观看的高清视频| 亚洲自偷自拍图片 自拍| 精品福利永久在线观看| 日韩欧美国产一区二区入口| 国产一区二区三区视频了| 久久久国产精品麻豆| 最近最新中文字幕大全电影3 | tocl精华| 97人妻天天添夜夜摸| 黄频高清免费视频| 色播在线永久视频| 中文字幕人妻丝袜一区二区| 欧美中文综合在线视频| 波多野结衣av一区二区av| 在线观看免费午夜福利视频| 久久久精品区二区三区| 亚洲黑人精品在线| 久久精品亚洲熟妇少妇任你| 成人免费观看视频高清| 两个人看的免费小视频| 日韩精品免费视频一区二区三区| 在线av久久热| 97在线人人人人妻| 精品久久蜜臀av无| 1024香蕉在线观看| 精品少妇内射三级| 交换朋友夫妻互换小说| 国产无遮挡羞羞视频在线观看| 三级毛片av免费| 丁香欧美五月| 最近最新中文字幕大全电影3 | 欧美激情高清一区二区三区| 午夜激情av网站| 欧美日韩视频精品一区| 精品亚洲乱码少妇综合久久| 国产精品久久电影中文字幕 | 老司机影院毛片| 国产精品国产高清国产av | 久久热在线av| 老司机影院毛片| 曰老女人黄片| 欧美人与性动交α欧美精品济南到| 精品人妻1区二区| 欧美中文综合在线视频| 捣出白浆h1v1| 欧美 亚洲 国产 日韩一| 久久久久国内视频| 自线自在国产av| 亚洲综合色网址| 午夜福利,免费看| 亚洲av片天天在线观看| 大型av网站在线播放| 丰满饥渴人妻一区二区三| 一区二区av电影网| av在线播放免费不卡| 丰满迷人的少妇在线观看| 国产免费现黄频在线看| 精品福利观看| 精品国产乱码久久久久久小说| 国产精品麻豆人妻色哟哟久久| 色尼玛亚洲综合影院| 欧美激情极品国产一区二区三区| 不卡av一区二区三区| 免费人妻精品一区二区三区视频| 十八禁高潮呻吟视频| 多毛熟女@视频| 亚洲精品粉嫩美女一区| 亚洲va日本ⅴa欧美va伊人久久| 嫁个100分男人电影在线观看| 精品福利观看| 亚洲熟女精品中文字幕| 亚洲欧美色中文字幕在线| 久久人妻av系列| 成在线人永久免费视频| 精品少妇一区二区三区视频日本电影| 热99re8久久精品国产| 久久精品熟女亚洲av麻豆精品| 色视频在线一区二区三区| 丝袜美腿诱惑在线| 建设人人有责人人尽责人人享有的| 国产精品秋霞免费鲁丝片| 国产又色又爽无遮挡免费看| 亚洲 欧美一区二区三区| 亚洲欧美一区二区三区黑人| 少妇 在线观看| 女警被强在线播放| 国产高清激情床上av| 免费一级毛片在线播放高清视频 | 日韩中文字幕欧美一区二区| 菩萨蛮人人尽说江南好唐韦庄| 两个人免费观看高清视频| 久久 成人 亚洲| 久久久久久久精品吃奶| 中文字幕制服av| 高清av免费在线| 国产欧美日韩精品亚洲av| 久久精品亚洲熟妇少妇任你| 国产欧美日韩一区二区精品| 99re在线观看精品视频| 国精品久久久久久国模美| cao死你这个sao货| 狠狠精品人妻久久久久久综合| 欧美另类亚洲清纯唯美| 五月开心婷婷网| 亚洲午夜精品一区,二区,三区| 美女主播在线视频| h视频一区二区三区| 久久 成人 亚洲| 亚洲国产欧美一区二区综合| 黑人巨大精品欧美一区二区蜜桃| av一本久久久久| 日韩视频一区二区在线观看| 国产精品亚洲一级av第二区| 性少妇av在线| kizo精华| 91九色精品人成在线观看| 桃红色精品国产亚洲av| 亚洲成人免费av在线播放| 久久久久久久久免费视频了| 亚洲精品中文字幕一二三四区 | 国产单亲对白刺激| 精品人妻熟女毛片av久久网站| 老司机深夜福利视频在线观看| 亚洲色图 男人天堂 中文字幕| 亚洲第一av免费看| 国产av国产精品国产| 母亲3免费完整高清在线观看| 丝瓜视频免费看黄片| 亚洲欧洲精品一区二区精品久久久| 黄色视频,在线免费观看| 建设人人有责人人尽责人人享有的| 成人国语在线视频| 精品人妻熟女毛片av久久网站| a级毛片黄视频| 国产成人av教育| 亚洲视频免费观看视频| 夜夜骑夜夜射夜夜干| 亚洲天堂av无毛| 国产主播在线观看一区二区| 9191精品国产免费久久| 国产色视频综合| 精品一区二区三卡| 一区二区三区国产精品乱码| 国产黄色免费在线视频| 丝袜喷水一区| 亚洲色图综合在线观看| 久久天躁狠狠躁夜夜2o2o| 男女午夜视频在线观看| 怎么达到女性高潮| 电影成人av| 男女午夜视频在线观看| 精品少妇久久久久久888优播| 欧美精品人与动牲交sv欧美| 十八禁高潮呻吟视频| 50天的宝宝边吃奶边哭怎么回事| 交换朋友夫妻互换小说| 99国产极品粉嫩在线观看| 淫妇啪啪啪对白视频| 亚洲av欧美aⅴ国产| 丝瓜视频免费看黄片| 9色porny在线观看| 69av精品久久久久久 | 法律面前人人平等表现在哪些方面| 免费高清在线观看日韩| 不卡一级毛片| √禁漫天堂资源中文www| 亚洲九九香蕉| 久久久久久久大尺度免费视频| 久久久久精品国产欧美久久久| 肉色欧美久久久久久久蜜桃| 老熟妇仑乱视频hdxx|