• <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值點偏移問題
    淺談某水電站引水隧洞的圍巖變形量測研究
    久久av网站| 老鸭窝网址在线观看| 亚洲成人免费电影在线观看 | 手机成人av网站| 色视频在线一区二区三区| 精品少妇内射三级| 国产一区二区 视频在线| 肉色欧美久久久久久久蜜桃| 亚洲图色成人| 不卡av一区二区三区| 精品少妇久久久久久888优播| 国产成人影院久久av| 久久精品国产亚洲av高清一级| 亚洲欧美精品自产自拍| 晚上一个人看的免费电影| 国产一区二区三区av在线| 国产在线观看jvid| videos熟女内射| 亚洲专区国产一区二区| 久久国产亚洲av麻豆专区| 无限看片的www在线观看| 久久久久久免费高清国产稀缺| 成人午夜精彩视频在线观看| 伊人亚洲综合成人网| 国产成人免费观看mmmm| 中文字幕精品免费在线观看视频| 在线观看人妻少妇| 免费观看人在逋| 亚洲国产av新网站| 久久久久久久久免费视频了| 丰满迷人的少妇在线观看| 婷婷色综合大香蕉| 久久人妻熟女aⅴ| 无遮挡黄片免费观看| 日韩 亚洲 欧美在线| 国精品久久久久久国模美| 国产精品人妻久久久影院| 亚洲美女黄色视频免费看| 桃花免费在线播放| 久久精品亚洲熟妇少妇任你| 极品少妇高潮喷水抽搐| 国产亚洲精品久久久久5区| 日韩制服骚丝袜av| 精品久久久精品久久久| 97在线人人人人妻| 久久人妻福利社区极品人妻图片 | 成人特级黄色片久久久久久久| 国产一区二区三区视频了| 亚洲一区中文字幕在线| 一卡2卡三卡四卡精品乱码亚洲| 日韩欧美免费精品| 国产亚洲欧美在线一区二区| 国产私拍福利视频在线观看| 久久狼人影院| 男女之事视频高清在线观看| www.精华液| 国产精品乱码一区二三区的特点| 国产高清激情床上av| 亚洲七黄色美女视频| 亚洲成国产人片在线观看| 在线国产一区二区在线| 天堂影院成人在线观看| av在线天堂中文字幕| 香蕉久久夜色| 99久久无色码亚洲精品果冻| 欧美黑人精品巨大| 正在播放国产对白刺激| 亚洲一卡2卡3卡4卡5卡精品中文| 亚洲av成人av| 亚洲精品国产区一区二| 国产成+人综合+亚洲专区| 日日爽夜夜爽网站| 99精品久久久久人妻精品| 亚洲aⅴ乱码一区二区在线播放 | 欧美精品啪啪一区二区三区| 日本一本二区三区精品| 欧美日韩一级在线毛片| 51午夜福利影视在线观看| 国产一区在线观看成人免费| 男男h啪啪无遮挡| 亚洲人成伊人成综合网2020| 丝袜美腿诱惑在线| 国产成人系列免费观看| 国产99白浆流出| 欧美日韩精品网址| 一边摸一边做爽爽视频免费| 欧美一区二区精品小视频在线| 成年版毛片免费区| 国产又爽黄色视频| 国产一卡二卡三卡精品| 美女大奶头视频| 欧美激情久久久久久爽电影| 国产精品永久免费网站| 一级作爱视频免费观看| 免费在线观看成人毛片| 国产极品粉嫩免费观看在线| 亚洲成a人片在线一区二区| 人人澡人人妻人| 女同久久另类99精品国产91| 黄色女人牲交| 欧美成狂野欧美在线观看| xxx96com| 久久欧美精品欧美久久欧美| 免费观看人在逋| 少妇的丰满在线观看| 中文资源天堂在线| 国产成人系列免费观看| av视频在线观看入口| 亚洲第一欧美日韩一区二区三区| 又紧又爽又黄一区二区| 亚洲激情在线av| 久久久久国内视频| 99国产精品99久久久久| 日韩中文字幕欧美一区二区| 97碰自拍视频| 91字幕亚洲| 日日摸夜夜添夜夜添小说| www.熟女人妻精品国产| 欧美成狂野欧美在线观看| 亚洲国产精品成人综合色| 中国美女看黄片| 高清在线国产一区| 久久久久久九九精品二区国产 | 丝袜美腿诱惑在线| 婷婷丁香在线五月| 一级毛片高清免费大全| 高潮久久久久久久久久久不卡| 国产精品亚洲av一区麻豆| 一个人观看的视频www高清免费观看 | 高清在线国产一区| 99久久99久久久精品蜜桃| 久久欧美精品欧美久久欧美| 一个人观看的视频www高清免费观看 | 巨乳人妻的诱惑在线观看| 欧美在线一区亚洲| 中文在线观看免费www的网站 | 在线永久观看黄色视频| 美女高潮到喷水免费观看| 十八禁网站免费在线| 亚洲国产精品sss在线观看| 老司机在亚洲福利影院| 一级毛片女人18水好多| 一级作爱视频免费观看| 亚洲成av人片免费观看| 欧美另类亚洲清纯唯美| 老汉色∧v一级毛片| 啦啦啦 在线观看视频| 亚洲成人久久爱视频| 午夜日韩欧美国产| 男人舔女人的私密视频| 男女做爰动态图高潮gif福利片| 久久国产乱子伦精品免费另类| 超碰成人久久| 国产高清videossex| 91大片在线观看| 成人亚洲精品av一区二区| 日本在线视频免费播放| 国产乱人伦免费视频| 18禁国产床啪视频网站| 黄色片一级片一级黄色片| 久久精品亚洲精品国产色婷小说| 他把我摸到了高潮在线观看| 美女 人体艺术 gogo| 亚洲一区二区三区不卡视频| 午夜免费观看网址| 国产欧美日韩一区二区精品| 丰满人妻熟妇乱又伦精品不卡| 国产久久久一区二区三区| 精品日产1卡2卡| 成人三级做爰电影| 国产av一区在线观看免费| 亚洲avbb在线观看| 午夜老司机福利片| 久久久国产成人免费| 国产色视频综合| 午夜精品久久久久久毛片777| 成人亚洲精品av一区二区| 女同久久另类99精品国产91| 叶爱在线成人免费视频播放| 国产av一区在线观看免费| 麻豆久久精品国产亚洲av| 91成年电影在线观看| 亚洲色图 男人天堂 中文字幕| 又大又爽又粗| 久久久国产欧美日韩av| 1024香蕉在线观看| 琪琪午夜伦伦电影理论片6080| 两性夫妻黄色片| 亚洲免费av在线视频| 日韩欧美在线二视频| 在线国产一区二区在线| 亚洲片人在线观看| 99在线视频只有这里精品首页| 欧美乱码精品一区二区三区| 男女床上黄色一级片免费看| 熟女少妇亚洲综合色aaa.| www.自偷自拍.com| 亚洲精华国产精华精| 国产精品野战在线观看| 国内揄拍国产精品人妻在线 | svipshipincom国产片| 久久香蕉精品热| 国产男靠女视频免费网站| 久久精品国产清高在天天线| 亚洲五月天丁香| 欧美一级a爱片免费观看看 | 精品免费久久久久久久清纯| 久久这里只有精品19| 精品国内亚洲2022精品成人| 男女午夜视频在线观看| 母亲3免费完整高清在线观看| 人人妻人人澡人人看| 老汉色av国产亚洲站长工具| 此物有八面人人有两片| av欧美777| 国产精品爽爽va在线观看网站 | 日本精品一区二区三区蜜桃| 久久婷婷人人爽人人干人人爱| 欧美日韩亚洲综合一区二区三区_| 欧美乱妇无乱码| 中出人妻视频一区二区| 国产激情偷乱视频一区二区| 曰老女人黄片| 亚洲国产精品999在线| 亚洲狠狠婷婷综合久久图片| 免费在线观看黄色视频的| 久久青草综合色| 亚洲欧美日韩无卡精品| 真人一进一出gif抽搐免费| 香蕉久久夜色| 在线天堂中文资源库| 久久精品国产99精品国产亚洲性色| 精品福利观看| 国内毛片毛片毛片毛片毛片| 久久青草综合色| 久久国产精品男人的天堂亚洲| 成年人黄色毛片网站| 美国免费a级毛片| 国产一卡二卡三卡精品| 特大巨黑吊av在线直播 | 精品久久久久久久末码| e午夜精品久久久久久久| 日韩一卡2卡3卡4卡2021年| 欧洲精品卡2卡3卡4卡5卡区| 丰满人妻熟妇乱又伦精品不卡| 级片在线观看| 99国产精品一区二区三区| 午夜视频精品福利| 午夜免费观看网址| 老鸭窝网址在线观看| 国产激情久久老熟女| 亚洲精品中文字幕在线视频| 国产成人精品无人区| 亚洲人成网站高清观看| 国产成人欧美| 美女午夜性视频免费| 亚洲久久久国产精品| av在线天堂中文字幕| 嫁个100分男人电影在线观看| 女警被强在线播放| 午夜免费鲁丝| 国内少妇人妻偷人精品xxx网站 | 88av欧美| 久久精品91无色码中文字幕| 窝窝影院91人妻| 日本黄色视频三级网站网址| 搡老妇女老女人老熟妇| 观看免费一级毛片| 热re99久久国产66热| 亚洲中文字幕一区二区三区有码在线看 | 黄色 视频免费看| 亚洲真实伦在线观看| 日本精品一区二区三区蜜桃| 宅男免费午夜| 国产伦人伦偷精品视频| 可以免费在线观看a视频的电影网站| 国产黄片美女视频| 欧美日韩中文字幕国产精品一区二区三区| 午夜成年电影在线免费观看| 成人三级黄色视频| 99久久精品国产亚洲精品| 欧美三级亚洲精品| 精品久久久久久久末码| 欧美激情久久久久久爽电影| 黄片小视频在线播放| 亚洲黑人精品在线| 国产亚洲欧美98| 波多野结衣巨乳人妻| 婷婷亚洲欧美| 亚洲专区中文字幕在线| 久久精品国产清高在天天线| 国产成人精品久久二区二区91| 午夜精品久久久久久毛片777| 一卡2卡三卡四卡精品乱码亚洲| 成人精品一区二区免费| 国产成人欧美| 久久精品国产清高在天天线| 青草久久国产| 日本一本二区三区精品| 午夜视频精品福利| 母亲3免费完整高清在线观看| 91在线观看av| videosex国产| 窝窝影院91人妻| 久久精品国产清高在天天线| 久久伊人香网站| 无遮挡黄片免费观看| 在线观看日韩欧美| 不卡一级毛片| 国产色视频综合| 白带黄色成豆腐渣| 国产成人精品无人区| 久久国产亚洲av麻豆专区| 亚洲自拍偷在线| 国产成人精品久久二区二区91| 久久久精品国产亚洲av高清涩受| 日本一本二区三区精品| 亚洲午夜精品一区,二区,三区| 久久欧美精品欧美久久欧美| 夜夜爽天天搞| 可以免费在线观看a视频的电影网站| 黄色成人免费大全| 久久精品国产清高在天天线| 别揉我奶头~嗯~啊~动态视频| 国产1区2区3区精品| 亚洲中文字幕一区二区三区有码在线看 | 夜夜爽天天搞| 色老头精品视频在线观看| 天天添夜夜摸| 国产精品电影一区二区三区| 国产精品二区激情视频| 丰满人妻熟妇乱又伦精品不卡| 久久国产精品男人的天堂亚洲| 18禁黄网站禁片午夜丰满| 精品国产超薄肉色丝袜足j| 久久久久久久午夜电影| 中出人妻视频一区二区| 在线观看舔阴道视频| 韩国精品一区二区三区| 1024香蕉在线观看| 中国美女看黄片| 午夜成年电影在线免费观看| 欧美精品亚洲一区二区| 免费av毛片视频| 亚洲一码二码三码区别大吗| 最新在线观看一区二区三区| av电影中文网址| 婷婷六月久久综合丁香| 女警被强在线播放| 亚洲精品中文字幕在线视频| 91国产中文字幕| 国语自产精品视频在线第100页| 国产伦人伦偷精品视频| 精品卡一卡二卡四卡免费| 两个人免费观看高清视频| 欧美色视频一区免费| 国产精品电影一区二区三区| 变态另类成人亚洲欧美熟女| 无遮挡黄片免费观看| 国产精品爽爽va在线观看网站 | 成年女人毛片免费观看观看9| www日本在线高清视频| avwww免费| 香蕉国产在线看| 曰老女人黄片| 亚洲欧美一区二区三区黑人| 欧美日本视频| 色综合站精品国产| 久久青草综合色| 亚洲精品久久成人aⅴ小说| 在线看三级毛片| 亚洲成人久久性| www.www免费av| 日本熟妇午夜| or卡值多少钱| 一级毛片高清免费大全| www.999成人在线观看| 国产久久久一区二区三区| 999精品在线视频| 亚洲精品美女久久av网站| 久久久久九九精品影院| 亚洲av电影不卡..在线观看| 99精品在免费线老司机午夜| 亚洲国产欧美网| 不卡一级毛片| www.999成人在线观看| 久9热在线精品视频| 日日爽夜夜爽网站| 国产又爽黄色视频| 97超级碰碰碰精品色视频在线观看| 国产精品免费视频内射| 夜夜爽天天搞| 老熟妇乱子伦视频在线观看| 狠狠狠狠99中文字幕| 亚洲中文日韩欧美视频| 久久国产精品男人的天堂亚洲| 一边摸一边抽搐一进一小说| 99热6这里只有精品| 欧美久久黑人一区二区| 成年版毛片免费区| 亚洲精品在线观看二区| 欧美成人性av电影在线观看| 在线av久久热| 无限看片的www在线观看| 国产不卡一卡二| 黄频高清免费视频| 午夜老司机福利片| 日日爽夜夜爽网站| 亚洲熟女毛片儿| 国产成年人精品一区二区| 久久精品91无色码中文字幕| 久久国产精品人妻蜜桃| 亚洲精品久久成人aⅴ小说| 精品熟女少妇八av免费久了| 亚洲一区二区三区不卡视频| 一级a爱片免费观看的视频| 日韩欧美免费精品| 精品欧美国产一区二区三| 黄片播放在线免费| 欧美日韩黄片免| 18禁黄网站禁片免费观看直播| 露出奶头的视频| 可以在线观看的亚洲视频| 成人一区二区视频在线观看| 亚洲三区欧美一区| 国产欧美日韩精品亚洲av| 国产视频内射| 久久精品国产综合久久久| 在线观看66精品国产| 19禁男女啪啪无遮挡网站| 两个人看的免费小视频| 一a级毛片在线观看| 国产精品美女特级片免费视频播放器 | 天堂影院成人在线观看| 一级a爱视频在线免费观看| 国产欧美日韩精品亚洲av| 91麻豆av在线| 18禁观看日本| 久久精品国产亚洲av高清一级| 亚洲精品国产一区二区精华液| 亚洲国产欧美日韩在线播放| 免费看美女性在线毛片视频| 中文字幕精品免费在线观看视频| 母亲3免费完整高清在线观看| 又紧又爽又黄一区二区| 亚洲色图 男人天堂 中文字幕| 国产精品综合久久久久久久免费| av片东京热男人的天堂| 啦啦啦免费观看视频1| 变态另类成人亚洲欧美熟女| 久久 成人 亚洲| 国产真人三级小视频在线观看| 成人手机av| 日本熟妇午夜| 久久久久久国产a免费观看| 高清在线国产一区| 热99re8久久精品国产| 老司机靠b影院| 别揉我奶头~嗯~啊~动态视频| 一边摸一边做爽爽视频免费| 日韩av在线大香蕉| 级片在线观看| 国产av在哪里看| 男人的好看免费观看在线视频 | 久久午夜亚洲精品久久| 中文字幕高清在线视频| 日韩欧美一区二区三区在线观看| 成年免费大片在线观看| 日本五十路高清| 国产成人欧美在线观看| 午夜福利一区二区在线看| 人人澡人人妻人| 亚洲精品国产精品久久久不卡| 99久久99久久久精品蜜桃| 成人亚洲精品一区在线观看| 脱女人内裤的视频| 亚洲激情在线av| 欧美一级毛片孕妇| 国产99白浆流出| 极品教师在线免费播放| 最近最新中文字幕大全电影3 | 无遮挡黄片免费观看| 久久午夜亚洲精品久久| 变态另类成人亚洲欧美熟女| a级毛片a级免费在线| 日本一本二区三区精品| 国产精品久久久久久精品电影 | 男人舔女人下体高潮全视频| 国产在线精品亚洲第一网站| 免费高清在线观看日韩| 亚洲自偷自拍图片 自拍| 色播亚洲综合网| 久久久久久免费高清国产稀缺| 99热6这里只有精品| 亚洲国产精品sss在线观看| 亚洲男人的天堂狠狠| 18禁国产床啪视频网站| 精品高清国产在线一区| 两个人视频免费观看高清| 免费无遮挡裸体视频| 免费观看人在逋| 成人一区二区视频在线观看| 婷婷六月久久综合丁香| 亚洲一区高清亚洲精品| 国产99白浆流出| 久久亚洲真实| 人人妻人人看人人澡| 日韩欧美一区二区三区在线观看| 久久久国产成人精品二区| 国产成人欧美| 亚洲五月婷婷丁香| 国产乱人伦免费视频| 最近最新中文字幕大全免费视频| 中文亚洲av片在线观看爽| 可以在线观看的亚洲视频| 啦啦啦免费观看视频1| 两个人看的免费小视频| 亚洲狠狠婷婷综合久久图片| 黄色成人免费大全| 国产国语露脸激情在线看| 18禁观看日本| 国产真实乱freesex| 999久久久精品免费观看国产| 亚洲一码二码三码区别大吗| 日日爽夜夜爽网站| 精品国产一区二区三区四区第35| aaaaa片日本免费| 日韩欧美 国产精品| 桃红色精品国产亚洲av| 精品福利观看| 国产成人精品久久二区二区91| 亚洲午夜理论影院| 亚洲精品国产精品久久久不卡| 精品国内亚洲2022精品成人| 国产精品,欧美在线| 成人欧美大片| 人成视频在线观看免费观看| 麻豆久久精品国产亚洲av| 女人被狂操c到高潮| 国产一卡二卡三卡精品| 热re99久久国产66热| 欧美黑人巨大hd| 12—13女人毛片做爰片一| 婷婷丁香在线五月| 长腿黑丝高跟| 日韩国内少妇激情av| 一本一本综合久久| 精华霜和精华液先用哪个| xxx96com| 国产伦人伦偷精品视频| 亚洲中文字幕一区二区三区有码在线看 | xxxwww97欧美| 亚洲成人久久性| 亚洲五月天丁香| 国产精品亚洲美女久久久| 国产黄a三级三级三级人| 天堂影院成人在线观看| 欧美av亚洲av综合av国产av| 一级片免费观看大全| 99国产精品99久久久久| 国产三级黄色录像| 精品久久久久久久人妻蜜臀av| 最近最新中文字幕大全免费视频| 精品国产亚洲在线| 好男人在线观看高清免费视频 | 一级作爱视频免费观看| 国产亚洲精品第一综合不卡| 一进一出抽搐动态| 一级a爱视频在线免费观看| 99riav亚洲国产免费| 日本成人三级电影网站| 丝袜人妻中文字幕| 在线十欧美十亚洲十日本专区| 成人免费观看视频高清| 高潮久久久久久久久久久不卡| 色综合亚洲欧美另类图片| 亚洲欧美精品综合久久99| 国产亚洲精品第一综合不卡| 亚洲成国产人片在线观看| 无遮挡黄片免费观看| 国产成人欧美在线观看| 亚洲成人精品中文字幕电影| 亚洲成a人片在线一区二区| 嫩草影视91久久| 亚洲成人久久爱视频| 伊人久久大香线蕉亚洲五| 午夜激情福利司机影院| 神马国产精品三级电影在线观看 | 好看av亚洲va欧美ⅴa在| 日韩大码丰满熟妇| 日韩成人在线观看一区二区三区| 国内精品久久久久精免费| 国产精品九九99| 久久精品国产综合久久久| 亚洲精品中文字幕一二三四区| 在线十欧美十亚洲十日本专区| 欧美激情高清一区二区三区| 国产亚洲精品一区二区www| 久久香蕉国产精品| 久久热在线av| 欧美黄色片欧美黄色片| 久久热在线av| 天堂√8在线中文| 欧美最黄视频在线播放免费| av有码第一页| 亚洲av中文字字幕乱码综合 | 中文字幕最新亚洲高清| 亚洲国产欧洲综合997久久, | 中亚洲国语对白在线视频| 亚洲全国av大片| 国产片内射在线| 午夜亚洲福利在线播放| 老熟妇仑乱视频hdxx| 可以免费在线观看a视频的电影网站| 少妇被粗大的猛进出69影院|