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

    邊坡動力破壞機(jī)理的振動臺試驗研究

    2014-03-06 05:45:12劉振平劉福海趙顯波
    關(guān)鍵詞:振動臺曲率滑動

    劉振平,劉福海,趙顯波

    (1.黑龍江工程學(xué)院 土木與建筑工程學(xué)院,哈爾濱 150050;

    2.東北電力設(shè)計院,長春130021;3.黑龍江省水利科學(xué)研究院,哈爾濱 150080)

    在土木工程和水利工程中,有眾多的自然邊坡、人工邊坡如堤壩、路基路塹等邊坡工程。中國的許多重大工程如大型水電站、南水北調(diào)西線工程、西部山區(qū)高速公路等均處于強(qiáng)震區(qū),特別是2008年汶川“5·12”地震以來,邊坡動力破壞機(jī)理引起學(xué)者們的廣泛關(guān)注,成為巖土地震工程領(lǐng)域中重要的研究課題之一。有關(guān)邊坡地震的原型觀測資料很少,所以室內(nèi)振動臺試驗就成為研究地震作用下邊坡動力問題的重要手段之一。Lin等[1]采用大型振動臺試驗研究土坡在地震作用下的響應(yīng),認(rèn)為砂土邊坡破壞面較淺,只是坡體表面的破壞。徐光興等[2-3]通過振動臺模型試驗研究了混合土邊坡的動力特性,并結(jié)合數(shù)值分析結(jié)果認(rèn)為在強(qiáng)震作用下單一均質(zhì)土層邊坡的破壞模式仍然是沿著某一弧形潛在滑動面失穩(wěn)。許強(qiáng)等[4]通過多組振動臺試驗認(rèn)為均質(zhì)土坡和塊狀巖質(zhì)斜坡的動力破壞模式都是坡頂拉裂,中下部剪切滑移破壞。陳新民等[5-6]采用振動臺試驗研究了粘土邊坡動力特性、動力反應(yīng)和宏觀變形,發(fā)現(xiàn)在動力作用下邊坡坡頂和坡腳出現(xiàn)裂縫。葉海林等[7]巖質(zhì)邊坡大型振動臺模型試驗,結(jié)果表明地震滑動面為上部拉裂縫和下部剪切滑移面形成貫通的破裂面。劉君等[8]把PIV技術(shù)應(yīng)用在土質(zhì)邊坡振動臺模型試驗中,獲取了邊坡破壞的完整過程。Wang等[9]在砂土邊坡振動臺試驗中采用PIV技術(shù)獲取了邊坡表面的位移場,探討了邊坡動力失穩(wěn)機(jī)理。Wartman等[10-11]通過多組粘土邊坡振動臺模型試驗研究了邊坡產(chǎn)生的永久位移,并與基于峰值和殘余強(qiáng)度的Newmark公式計算結(jié)果進(jìn)行了比較。Katz等[12]采用振動沙箱研究了砂土邊坡破壞的類型和頻率的關(guān)系。以上針對邊坡的動力特性、宏觀破壞現(xiàn)象及個別點的位移進(jìn)行的研究,而沒有涉及邊坡的位移場,劉君、Wang等雖用PIV技術(shù)得到了整個邊坡的位移場,但沒有對應(yīng)變場進(jìn)行研究。

    應(yīng)變是直接對含有噪聲的離散位移進(jìn)行差分計算得到的,那么即使微小的位移測量誤差也會被急劇地放大,使計算的應(yīng)變不可靠[13]。一般的數(shù)據(jù)擬合平滑方法也不能滿足應(yīng)變計算的精度要求,潘兵等[13-14]提出了局部最小二乘擬合的應(yīng)變計算方法,其原理同Savitzky-Golay平滑方法。該方法能夠在一定程度上去除原始位移場的噪聲,缺點是由于每個位移點都要一個計算窗口平滑,計算量巨大,并且對于非均勻應(yīng)變場,應(yīng)變計算窗口的大小顯著影響應(yīng)變的計算精度,目前還沒有合適的方法來選擇最佳的窗口大小。本文采用數(shù)字圖像位移和有限元數(shù)據(jù)平滑技術(shù),得到了整個邊坡模型的位移場和應(yīng)變場,同時探討了土質(zhì)邊坡的地震破壞機(jī)理。

    1 基于有限元數(shù)據(jù)平滑的應(yīng)變計算

    數(shù)字圖像位移測量技術(shù)現(xiàn)在已經(jīng)成熟,采用互相關(guān)圖像匹配技術(shù)獲取位移場后,再用有限元數(shù)據(jù)平滑方法對位移場進(jìn)行處理后,得到應(yīng)變場。

    1.1 基本原理

    基本思路是引入泛函f(F)[15-19],使該泛函取得極小值的自變函數(shù)F就是要求的光滑位移函數(shù)。

    根據(jù)試驗數(shù)據(jù)點的位置分布,把區(qū)域A劃分為若干單元,在每一個單元內(nèi)用形函數(shù)逼近光滑函數(shù)F,這個過程與有限元方法是一樣的。

    在第j個單元內(nèi),把光滑函數(shù)F表示成式中:N j為第j個單元的形函數(shù);uej為第j個單元的節(jié)點自由度。

    式(9)是有限元的標(biāo)準(zhǔn)格式,其中K=K1+λK2,λ是平滑系數(shù),U表示所有單元的節(jié)點自由度。平滑系數(shù)的確定方法在下文專門論述。

    只要確定了平滑系數(shù),就可從式(9)中解出U,再從式(2)中得到平滑位移函數(shù)u,其應(yīng)變的計算方法如下:

    有限變形條形下的格林(Green)應(yīng)變[20](以受壓為正)為

    式中:u和v分別為上述所求的x和y方向的位移函數(shù);εx和εy分別為x和y方向的正應(yīng)變;γxy為工程剪應(yīng)變;γmax為最大剪應(yīng)變。

    小變形條件下柯西(Cauchy)應(yīng)變的計算可不考慮上述應(yīng)變表達(dá)式中的二次項。

    1.2 形函數(shù)的確定

    采用三角形單元,以便適應(yīng)任意形狀的數(shù)據(jù)區(qū)域。同時為提高插值函數(shù)的光滑可導(dǎo)性,采用的是五次多項式的單元插值函數(shù)[15],其表達(dá)式為式中:X為變量項;C為系數(shù)列向量。由于沿三角形單元三條邊滿足法向?qū)?shù)連續(xù)的條件,21個系數(shù)可縮減為18個。

    1.3 平滑系數(shù)的確定

    為了能夠自動得到最佳的平滑系數(shù),Golub等[21]、Craven等[22]及 Bates等[23]提 出 了 廣 義 交 互驗證GCV方法,建立了關(guān)于λ的GCV函數(shù)X是n階滿秩陣,K2是奇異陣,其秩為r。擬總體剛度陣K1是對稱正定陣,因此式(15)成立,可以證明X是滿秩陣,因此y也能夠由式(16)式唯一確定。

    求取廣義交互驗證函數(shù)GCV的最小值,就可以得到最佳平滑系數(shù)λ,具體計算方法詳見文獻(xiàn)[24]。

    2 試驗設(shè)備及模型設(shè)計

    2.1 試驗設(shè)備

    試驗在大連理工大學(xué)工程抗震試驗室的水平與豎直雙向水下振動臺上進(jìn)行,其主要性能參數(shù)如下:數(shù)字式控制方式,臺面尺寸3 m×4 m,最大載重100 k N,工作頻率0.1~50 Hz;滿載工況下最大加速度±1.0g(水 平 向)、±0.7g(豎 向)、最 大 速 度±50 cm/s(水平向)、±35 cm/s(豎向)、最大位移±75 mm(水平向)、±50 mm(豎向)。

    2.2 試驗?zāi)P?/h3>

    在鋼制模型箱內(nèi)堆制模型進(jìn)行試驗,為便于觀察和圖像采集,一側(cè)為有機(jī)玻璃。模型壩壩高1.3 m,邊坡坡率1∶1.6,其模型如圖1所示。模型填土基本性能參數(shù)見表1。

    表1 試驗填土基本性能參數(shù)

    圖1 邊坡振動臺模型示意圖(單位:cm)

    試驗加載的是頻率10 Hz的正弦增幅波(喇叭波),如圖2所示。從0開始逐漸增加,直至邊坡完全失穩(wěn)。

    采用的圖像采集設(shè)備是CANON EOS 450D高清數(shù)碼相機(jī),最大分辨率4 272×2 848像素,配備10~22 mm廣角鏡頭,可以在距模型2 m處拍攝整個模型的圖片。存儲為JPG格式的圖片文件,最大連拍速度約3.3幀/s。振動過程中連續(xù)的采集圖像,用于數(shù)字圖像位移和應(yīng)變識別。

    圖2 輸入的加速度時程

    3 試驗結(jié)果分析

    3.1 模型的位移應(yīng)變計算流程

    開發(fā)了一套數(shù)字圖像測量分析程序,可以批處理多種格式的圖片,可以得到一系列的位移場和應(yīng)變場,還能以等值線圖、云圖和矢量圖等形式輸出各種位移場、應(yīng)變場和速度場。其計算流程和步驟如下:

    位移計算流程:圖像像素塊劃分→計算總位移→計算模型箱位移→計算相對位移(總位移減去模型箱剛體位移)→刪除計算的位移壞點→輸出位移場等后處理。

    得到位移場后,用有限元平滑位移場再微分計算應(yīng)變場。其計算流程如下:生成單元和節(jié)點信息→計算單元剛度陣和單元節(jié)點荷載列陣→集成總體剛度陣和總體節(jié)點荷載列陣→計算平滑系數(shù)→求解節(jié)點未知參數(shù)→平滑位移場微分求應(yīng)變→輸出應(yīng)變場等后處理。

    3.2 邊坡模型的破壞過程描述

    邊坡模型破壞過程基本可分3個階段,第1階段是整體變形階段(圖3(a)),當(dāng)輸入加速度較小時,整個模型以均勻沉降變形為主,同時伴有相對于模型箱的水平往復(fù)運動;第2階段是滑移變形階段(圖3(b)),當(dāng)輸入加速度大約在0.1g時,模型開始沿邊坡一側(cè)向下滑移,隨著輸入加速度的增大,滑移量也逐漸增大,在輸入加速度在0.45g時,坡頂出現(xiàn)張拉裂縫,坡體內(nèi)部也出現(xiàn)剪切裂隙,形成明顯的滑動帶,邊坡達(dá)到極限破壞狀態(tài);第3階段是破壞階段,滑動帶上的裂隙使得滑動體與邊坡相對分離,滑動體加速下滑,最后整個坡體坍塌,邊坡完全破壞。

    從圖3(c)可以看出,邊坡模型動力失穩(wěn)存在深層滑動帶,滑動面大概呈圓弧狀。

    3.3 位移分析

    圖4、圖5和圖6分別是8.42 s(對應(yīng)輸入加速度為0.45g)時即在極限破壞狀態(tài)時的水平位移、豎向位移和總位移的等值線圖,其水平位移出現(xiàn)在坡面中部,豎向位移在坡頂附近,水平位移和總位移等值線基本平行于坡面,均符合一般的邊坡破壞時的位移分布規(guī)律。

    圖3 邊坡模型位移矢量

    圖4 8.42 s時(輸入加速度0.45 g)時水平位移場(單位:mm)

    圖5 8.42 s時(輸入加速度0.45 g)時豎向位移場(單位:mm)

    圖6 8.42 s時(輸入加速度0.45 g)時總位移場(單位:mm)

    3.4 應(yīng)變分析

    從圖7和圖8最大剪應(yīng)變和水平方向應(yīng)變分布上可以看出,坡體中部到坡腳是剪切破壞,坡頂一定深度是拉剪破壞。其變形破壞模式與汶川地震觀察到的邊坡破壞現(xiàn)象相吻合[25]:土坡主要是坡頂向下一定深度內(nèi)的拉破壞,坡腳向上延伸形成剪切滑移帶,最終二者連通形成貫通的破裂面。這也與許強(qiáng)等[4]的振動臺試驗結(jié)果以及鄭穎人等[26-27]用 FLAC程序做的數(shù)值分析結(jié)果基本一致。

    圖7 8.42 s時(輸入加速度0.45 g)時水平應(yīng)變場(單位:%)

    8 8.42 s時(輸入加速度0.45 g)時最大剪應(yīng)變場(單位:%)

    3.5 位移應(yīng)變時程分析

    A點位于滑動體上,B點位于滑動體外(圖1),從圖9的位移時程可以看出,位移隨時間是逐漸增加的,沒有明顯的突變點,說明邊坡在地震作用下是漸進(jìn)式破壞的。

    圖9 位移時程曲線

    由于位移時程曲線沒有明顯的突變拐點,如何判斷模型破壞的時刻就成為一個難題。地震作用下邊坡形成滑移帶后,滑動體與坡體相對分離,由于滑移帶的阻隔,滑動體以外的壩體振動對其影響將大大減小,滑動體將比較平穩(wěn)的運動,而不會隨輸入地震波而振動。所以可以把位移時程曲線的曲率做為判斷壩坡失穩(wěn)滑移的一個物理量,這里所指的曲率是有正負(fù)之分的,正號表示是凹曲線,負(fù)號表示是凸曲線,可定義為廣義曲率。

    式中:κ為廣義曲率;y′是位移的一階導(dǎo)數(shù);y″是位移的二階導(dǎo)數(shù)。

    從A點的位移時程曲線曲率(圖10)可以看出,在8.42 s之前,曲率以零線為中心上下波動,說明在此之前沒有發(fā)生連續(xù)的滑動,而在此之后,曲率迅速趨于零,不再發(fā)生波動,該點此時已經(jīng)失穩(wěn)破壞,曲率發(fā)生突變表明滑移帶已經(jīng)貫通形成,邊坡整體失穩(wěn),達(dá)到了極限抗震狀態(tài)。位于滑動體以外的B點,曲率自始自終都處于振動狀態(tài),沒有發(fā)生失穩(wěn)破壞。由此可知,用位移時程曲線曲率做為判斷邊坡動力破壞的物理量是可行的。

    圖10 位移時程曲線的曲率

    4 結(jié) 論

    將有限元數(shù)據(jù)平滑方法引入到數(shù)字圖像位移應(yīng)變測量中是可行的,得到了邊坡振動臺模型整個試驗過程的位移場和應(yīng)變場。初步表明土質(zhì)邊坡的變形是漸進(jìn)式的,坡體中部到坡腳是剪切破壞,坡頂一定深度是拉剪破壞,破壞時有深層的圓弧狀滑動面,用位移時程的廣義曲率做為判斷邊坡動力破壞的物理量是可行的。

    [1]Lin M L,Wang K L.Seismic slope behavior in a largescale shaking table model test [J].Engineering Geology,2006,86(2/3):118-133.

    [2]徐光興,姚令侃,高召寧,等.邊坡動力特性與動力響應(yīng)的大型振動臺模型試驗研究[J].巖石力學(xué)與工程學(xué)報,2008,27(3):624-631.

    Xu G X,Yao L K,Gao Z N,et al.Large-scale shaking table model test study on dynamic characteristics and dynamic responses of slope[J].Chinese Journal of Rock Mechanics and Engineering,2008,27(3):624-631.

    [3]徐光興,姚令侃,李朝紅,等.邊坡地震動力響應(yīng)規(guī)律及地震動參數(shù)研究[J].巖土工程學(xué)報,2008,30(6):918-923.

    Xu G X,Yao L K,Li C G,et al.Dynamic response of slopes under earthquakes and influence of ground motion parameters[J].Chinese Journal of Geotechnical Engineering,2008,30(6):918-923.

    [4]許強(qiáng),陳建君,馮文凱,等.斜坡地震響應(yīng)的物理模擬試驗研究[J].四川大學(xué)學(xué)報:工程科學(xué)版,2009,41(3):266-272.

    Xu J,Chen J J,F(xiàn)eng W K,et al.Study of the seismic response of slopes by physical modeling [J].Journal of Sichuan University:Engineering Science Edition,2009,41(3):266-272.

    [5]陳新民,沈建,魏平,等.下蜀土邊坡地震穩(wěn)定性的大型振動臺試驗研究(I)-模型試驗設(shè)計[J].防災(zāi)減災(zāi)工程學(xué)報,2010,30(5):492-502.

    Chen X M,Shen J,Wei P,et al.Large-scale shaking table test of seismic stability of Xiashu loess slopes(I):design of model test[J].Journal of Disaster Prevention and Mitigation Engineering,2010,30(5):492-502.

    [6]陳新民,沈建,魏平,等.下蜀土邊坡地震穩(wěn)定性的大型振動臺試驗研究(II)-試驗結(jié)果及分析[J].防災(zāi)減災(zāi)工程學(xué)報,2010,30(6):587-594.

    Chen X M,Shen J,Wei P,et al.Large-scale shaking table test of seismic stability of Xiashu loess slopes(II):analysis of test results [J].Journal of Disaster Prevention and Mitigation Engineering,2010,30(6):587-594.

    [7]葉海林,鄭穎人,杜修力,等.邊坡動力破壞特征的振動臺模型試驗與數(shù)值分析[J].土木工程學(xué)報,2012,45(9):128-135.

    Ye H L,Zheng Y R,Du X L,et al.Shaking table model test and numerical analysis on dynamicfailure characteristics of slope[J].China Civil Engineering Journal,2012,45(9):128-135.

    [8]劉君,劉福海,孔憲京,等.PIV技術(shù)在大型振動臺模型試驗中的應(yīng)用[J].巖土工程學(xué)報,2010,32(3):368-374.

    Liu J,Liu F H,Kong X J,et al.Application of PIV in largescale shaking table model tests [J].Chinese Journal of Geotechnical Engineering,2010,32(3):368-374.

    [9]Wang K L,Lin M L.Initiation and displacement of landslide induced by earthquake-a study of shakingtable model slope test[J].Engineering Geology,2011,122:106-114.

    [10]Wartman J.Physical model studies of seismically induced deformations in slopes [D].California:University of California-Berkeley,1999.

    [11]Wartman J,Seed R B,Bray J D.Shaking table modeling of seismically induced deformations in slopes [J].Journal of Geotechnical and Geoenvironmental Engineering,2005,131(5):610-622.

    [12]Katz O,Aharonov E.Andslides in vibrating sand box:what controls types of slope failure and frequencymagnitude relations?[J].Earth and Planetary Science Letters,2006,247(3):280-294.

    [13]潘兵,謝惠民.數(shù)字圖像相關(guān)中基于位移場局部最小二乘擬合的全場應(yīng)變測量[J].光學(xué)學(xué)報,2007,27(11):1980-1986.

    Pan B,Xie H M.Full-field strain measurement based on least-square fitting of local displacement for digital image correlation method[J].Acta Optica Sinica,2007,27(11):1980-1986.

    [14]Pan B,Xie H M,Guo Z Q,et al.Full-field strain measurement using a two-dimensional Savitzky-Golay digital differentiator in digital image correlation [J].Optical Engineering,2007,46(3):1-10.

    [15]Segalman D J,Woyak D B,Rowlands R E.Smooth spline-like finite-element differentiation of full-field experimental data over arbitrary geometry [J].Experimental Mechanics,1979,19(12):429-437.

    [16]Engelstad M J,Chambless D A,Swinson W F,et al.Hybrid stress analysis of vibrating plates using holographic interferometry and finite elements [J].Experimental Mechanics,1987,27(1):23-30.

    [17]Feng Z,Rowlands R E.Continuous full-field representation and differentiation of three-dimensional experimental vector data[J].Computers and Structures,1987,26(6):979-990.

    [18]Sutton M A,Turner J L,Burck H A,et al.Full-field repersentation of discretely smapled surface deformational for displacement and strain analysis [J].Experimental Mechnaics,1991,31(2):168-177.

    [19]Freese C E,Gee L.A multilevel treatment of moiré fringe data using finite elements [J].Experimental Mechanics,1999,39(4):304-310.

    [20]陸明萬,羅學(xué)富.彈性理論基礎(chǔ)[M].北京:清華大學(xué)出版社,2001:38-43.

    [21]Golub G H,Heath M,Wahba G.Generalized cross validation as a method for choosing a good ridge parameter[J].Technometrics,1979,21(2):215-223.

    [22]Craven P,Wahba G.Smoothing noisy data with spline functions:estimating the correct degree of smoothing by the method of generalized cross-validation [J].Numerische Mathematik,1979,31:377-403.

    [23]Bates D M,Lindstrom M J,Wahba G,et al.Gcvpak-routines for generalized cross validation [J].Communications in Statistics-Simulation and Computation,1987,16(1):263-297.

    [24]Kent J T,Mohammadzadeh M.Global optimization of the generalized cross-validation criterion [J].Statistics and Computing,2000,10(3):231-236.

    [25]許強(qiáng),董秀軍.汶川地震大型滑坡成因模式[J].中國地質(zhì)大學(xué)學(xué)報:地球科學(xué),2011,36(6):1134-1142.

    Xu J,Dong X J.Genetic types of large-scale landslides induced by Wenchuan Earthquake[J].Journal of China University of Geosciences:Earth Science,2011,36(6):1134-1142.

    [26]鄭穎人,葉海林,黃潤秋.地震邊坡破壞機(jī)制及其破裂面的分析探討[J].巖石力學(xué)與工程學(xué)報,2009,28(8):1714-1723.

    Zheng Y R,Ye H L,Huang R Q.Analysis and discussion of failure mechanism and fracture surface of slope under earthquake [J].Chinese Journal of Rock Mechanics and Engineering,2009,28(8):1714-1723.

    [27]鄭穎人,葉海林,黃潤秋,等.邊坡地震穩(wěn)定性分析探討[J].地震工程與工程振動,2010,30(2):173-180.

    Zheng Y R,Ye H L,Huang R Q,et al.Study on the seismic stability analysis of a slope [J].Journal of Earthquake Engineering and Engineering Vibration,2010,30(2):173-180.

    (編輯王秀玲)

    猜你喜歡
    振動臺曲率滑動
    大曲率沉管安裝關(guān)鍵技術(shù)研究
    一類雙曲平均曲率流的對稱與整體解
    基于振動臺試驗的通信機(jī)柜地震易損性分析
    半正迷向曲率的四維Shrinking Gradient Ricci Solitons
    一種新型滑動叉拉花鍵夾具
    Big Little lies: No One Is Perfect
    大型液壓離心振動臺控制策略的仿真研究
    滑動供電系統(tǒng)在城市軌道交通中的應(yīng)用
    420 kV避雷器振動臺抗震試驗
    一種基于變換域的滑動聚束SAR調(diào)頻率估計方法
    人体艺术视频欧美日本| 亚洲成色77777| 九色亚洲精品在线播放| 七月丁香在线播放| 亚洲激情五月婷婷啪啪| 全区人妻精品视频| 国产精品人妻久久久久久| 国产深夜福利视频在线观看| 五月玫瑰六月丁香| 九九在线视频观看精品| 少妇 在线观看| √禁漫天堂资源中文www| 欧美日韩一区二区视频在线观看视频在线| 最新的欧美精品一区二区| 国产深夜福利视频在线观看| 超碰97精品在线观看| 国语对白做爰xxxⅹ性视频网站| 中文天堂在线官网| 高清视频免费观看一区二区| 日韩一区二区三区影片| 久久青草综合色| 日日摸夜夜添夜夜爱| 免费观看无遮挡的男女| 国产高清有码在线观看视频| 亚洲欧美日韩卡通动漫| 亚洲综合精品二区| 久久午夜综合久久蜜桃| 女人久久www免费人成看片| 最近的中文字幕免费完整| 国产综合精华液| 天天躁夜夜躁狠狠久久av| 亚洲美女搞黄在线观看| 精品少妇久久久久久888优播| 亚洲成人av在线免费| 日本黄色片子视频| 久久 成人 亚洲| 色哟哟·www| 亚洲伊人久久精品综合| 免费观看av网站的网址| 国产不卡av网站在线观看| 国产亚洲一区二区精品| 免费久久久久久久精品成人欧美视频 | 又大又黄又爽视频免费| 99久久精品一区二区三区| 亚洲av成人精品一二三区| 国产精品99久久久久久久久| 狂野欧美激情性bbbbbb| 寂寞人妻少妇视频99o| 精品一区二区三区视频在线| 免费播放大片免费观看视频在线观看| 男人爽女人下面视频在线观看| 97超碰精品成人国产| videossex国产| 只有这里有精品99| 精品一区二区三卡| 热re99久久国产66热| 男女免费视频国产| 日本欧美视频一区| 欧美 亚洲 国产 日韩一| 夫妻午夜视频| 亚洲熟女精品中文字幕| 日本-黄色视频高清免费观看| 丰满乱子伦码专区| 99国产综合亚洲精品| 亚洲精品乱码久久久v下载方式| 国产男人的电影天堂91| 亚洲国产日韩一区二区| 99热这里只有精品一区| 亚洲婷婷狠狠爱综合网| 一区在线观看完整版| 大香蕉久久网| 午夜激情av网站| 美女内射精品一级片tv| 一本—道久久a久久精品蜜桃钙片| 久久精品国产亚洲网站| 日韩强制内射视频| 日韩一区二区视频免费看| 熟女人妻精品中文字幕| 女人精品久久久久毛片| 国产极品天堂在线| 热99久久久久精品小说推荐| 久久久久精品久久久久真实原创| 性色avwww在线观看| 国产极品粉嫩免费观看在线 | 免费黄频网站在线观看国产| 熟女电影av网| 视频区图区小说| 国产精品国产三级专区第一集| 欧美xxxx性猛交bbbb| 婷婷成人精品国产| 精品久久久噜噜| 久久久久久久久久久免费av| 亚洲av日韩在线播放| 日本wwww免费看| 国产日韩一区二区三区精品不卡 | 亚洲不卡免费看| 十八禁网站网址无遮挡| 日韩伦理黄色片| 黄片播放在线免费| 三上悠亚av全集在线观看| 久久久精品免费免费高清| 日韩人妻高清精品专区| 99久久中文字幕三级久久日本| 最近最新中文字幕免费大全7| 中文字幕人妻丝袜制服| 久久狼人影院| 少妇熟女欧美另类| 熟妇人妻不卡中文字幕| 成人亚洲欧美一区二区av| 亚洲av男天堂| 亚洲国产最新在线播放| 免费人成在线观看视频色| 亚洲第一区二区三区不卡| 亚洲精品色激情综合| 久久国内精品自在自线图片| 少妇的逼水好多| 18在线观看网站| 亚洲情色 制服丝袜| 婷婷色av中文字幕| 老司机影院毛片| 久久精品久久久久久噜噜老黄| 久久精品久久精品一区二区三区| 国产精品蜜桃在线观看| 免费高清在线观看日韩| 国产男女内射视频| 国产无遮挡羞羞视频在线观看| 亚洲精品一二三| 在线观看免费日韩欧美大片 | 欧美激情 高清一区二区三区| 成人毛片a级毛片在线播放| 亚洲内射少妇av| 亚洲五月色婷婷综合| kizo精华| 亚洲国产av影院在线观看| 国语对白做爰xxxⅹ性视频网站| av线在线观看网站| 91久久精品电影网| 校园人妻丝袜中文字幕| 国产精品麻豆人妻色哟哟久久| 黄色一级大片看看| 97超视频在线观看视频| 午夜激情福利司机影院| 日韩熟女老妇一区二区性免费视频| 69精品国产乱码久久久| 春色校园在线视频观看| 午夜免费男女啪啪视频观看| 国产极品粉嫩免费观看在线 | 日本欧美视频一区| 国产精品偷伦视频观看了| av天堂久久9| 久久久久久久大尺度免费视频| 久久亚洲国产成人精品v| av播播在线观看一区| 中文字幕av电影在线播放| 自线自在国产av| 免费观看在线日韩| 久久影院123| 美女视频免费永久观看网站| 免费观看在线日韩| 天美传媒精品一区二区| 又黄又爽又刺激的免费视频.| 欧美日韩一区二区视频在线观看视频在线| 午夜激情久久久久久久| 亚洲av综合色区一区| 看十八女毛片水多多多| 一级片'在线观看视频| 少妇熟女欧美另类| 九九在线视频观看精品| 久久午夜综合久久蜜桃| 高清在线视频一区二区三区| 午夜日本视频在线| 王馨瑶露胸无遮挡在线观看| 国产精品人妻久久久久久| 3wmmmm亚洲av在线观看| 国产片特级美女逼逼视频| 国产精品不卡视频一区二区| 亚洲伊人久久精品综合| 亚洲欧美日韩另类电影网站| 国产高清不卡午夜福利| 国产日韩一区二区三区精品不卡 | 亚洲av中文av极速乱| 综合色丁香网| 免费人成在线观看视频色| 街头女战士在线观看网站| 久久国产亚洲av麻豆专区| 免费黄频网站在线观看国产| 国产精品偷伦视频观看了| 亚洲欧美日韩卡通动漫| 成年女人在线观看亚洲视频| 少妇被粗大猛烈的视频| 九色亚洲精品在线播放| videosex国产| 久久久精品94久久精品| 国产乱人偷精品视频| 寂寞人妻少妇视频99o| 18禁在线播放成人免费| 日本欧美视频一区| 最新中文字幕久久久久| 水蜜桃什么品种好| 免费高清在线观看日韩| 18禁在线播放成人免费| 男人添女人高潮全过程视频| 国产一区二区在线观看日韩| 黄色视频在线播放观看不卡| 亚洲人成77777在线视频| 国产男女内射视频| 国产乱人偷精品视频| 国产永久视频网站| 丰满少妇做爰视频| 亚洲,欧美,日韩| 久久久久久伊人网av| 亚洲国产精品成人久久小说| 国产欧美亚洲国产| 免费久久久久久久精品成人欧美视频 | 大话2 男鬼变身卡| 国产高清有码在线观看视频| 青春草视频在线免费观看| 久久99精品国语久久久| 女性被躁到高潮视频| 国产探花极品一区二区| 夫妻午夜视频| 97在线人人人人妻| 亚洲色图 男人天堂 中文字幕 | 99久国产av精品国产电影| 欧美激情 高清一区二区三区| 99热国产这里只有精品6| 2021少妇久久久久久久久久久| 久久毛片免费看一区二区三区| 亚洲熟女精品中文字幕| av女优亚洲男人天堂| 亚洲无线观看免费| 国产精品免费大片| 久久久久久久久久久丰满| a级毛片在线看网站| 亚洲经典国产精华液单| 777米奇影视久久| 成人二区视频| 日韩av在线免费看完整版不卡| 青春草亚洲视频在线观看| 高清在线视频一区二区三区| 亚洲天堂av无毛| 亚洲四区av| 一个人看视频在线观看www免费| 精品国产一区二区久久| 国产精品久久久久久精品古装| 中文字幕久久专区| av有码第一页| 久久久国产精品麻豆| 亚洲经典国产精华液单| 美女脱内裤让男人舔精品视频| 女性生殖器流出的白浆| 亚洲精品,欧美精品| 美女国产视频在线观看| av天堂久久9| 欧美日韩av久久| 婷婷色av中文字幕| 国产不卡av网站在线观看| 国产有黄有色有爽视频| 亚洲精品乱码久久久久久按摩| 日韩不卡一区二区三区视频在线| 超色免费av| 久久久久久伊人网av| 亚洲内射少妇av| 亚洲精品中文字幕在线视频| 男女啪啪激烈高潮av片| 国产精品熟女久久久久浪| 国产在线视频一区二区| 青春草亚洲视频在线观看| 天天影视国产精品| 国产成人免费观看mmmm| 久久久久久人妻| 一级,二级,三级黄色视频| 99久久综合免费| a级片在线免费高清观看视频| 色婷婷av一区二区三区视频| 香蕉精品网在线| 日韩人妻高清精品专区| 夜夜骑夜夜射夜夜干| 久久久久国产网址| 成人18禁高潮啪啪吃奶动态图 | 成人黄色视频免费在线看| 国产精品99久久99久久久不卡 | 最近2019中文字幕mv第一页| 永久免费av网站大全| av卡一久久| 人妻一区二区av| 大片电影免费在线观看免费| 欧美精品高潮呻吟av久久| 国产国语露脸激情在线看| 久久精品国产亚洲网站| 亚洲激情五月婷婷啪啪| 久久久久精品性色| 亚洲av日韩在线播放| 亚洲婷婷狠狠爱综合网| 亚洲中文av在线| 久久久久人妻精品一区果冻| 午夜视频国产福利| 精品久久国产蜜桃| 亚洲精品久久成人aⅴ小说 | 美女视频免费永久观看网站| 欧美日韩在线观看h| 九色亚洲精品在线播放| 国产永久视频网站| 九色亚洲精品在线播放| 天天躁夜夜躁狠狠久久av| 一区二区三区乱码不卡18| 国产一区二区三区综合在线观看 | 亚洲五月色婷婷综合| 我的老师免费观看完整版| 一区二区三区免费毛片| 女性生殖器流出的白浆| 亚洲精品久久午夜乱码| 亚洲,欧美,日韩| 亚洲精品国产色婷婷电影| 青春草视频在线免费观看| 国产成人免费观看mmmm| 国产亚洲av片在线观看秒播厂| 国产精品一国产av| 美女视频免费永久观看网站| 成人手机av| 久久久久国产精品人妻一区二区| 国产成人精品久久久久久| 精品国产国语对白av| 欧美日本中文国产一区发布| 人人妻人人添人人爽欧美一区卜| 亚洲人成网站在线观看播放| 激情五月婷婷亚洲| 午夜视频国产福利| 国产精品欧美亚洲77777| 国产一区有黄有色的免费视频| 少妇被粗大猛烈的视频| 狠狠精品人妻久久久久久综合| 婷婷色麻豆天堂久久| 曰老女人黄片| 欧美精品一区二区大全| 曰老女人黄片| 久久久亚洲精品成人影院| 国产精品熟女久久久久浪| 国产熟女欧美一区二区| 99九九在线精品视频| 九草在线视频观看| 亚洲,一卡二卡三卡| 精品人妻熟女毛片av久久网站| 男女国产视频网站| 少妇的逼水好多| 午夜久久久在线观看| 午夜影院在线不卡| 日韩三级伦理在线观看| 一级爰片在线观看| 最近最新中文字幕免费大全7| 熟妇人妻不卡中文字幕| 菩萨蛮人人尽说江南好唐韦庄| 国产无遮挡羞羞视频在线观看| 日韩人妻高清精品专区| 免费黄色在线免费观看| a级毛片在线看网站| 国产精品一区www在线观看| 午夜福利在线观看免费完整高清在| 好男人视频免费观看在线| 人人妻人人澡人人爽人人夜夜| 国产精品女同一区二区软件| 一本一本久久a久久精品综合妖精 国产伦在线观看视频一区 | 美女脱内裤让男人舔精品视频| 国产av国产精品国产| 国产一区二区三区综合在线观看 | 婷婷成人精品国产| 少妇丰满av| 欧美亚洲日本最大视频资源| 26uuu在线亚洲综合色| 国产精品久久久久久久久免| 国产一区二区在线观看日韩| 老熟女久久久| 色吧在线观看| 午夜影院在线不卡| 多毛熟女@视频| 午夜视频国产福利| 欧美性感艳星| 久久鲁丝午夜福利片| 天美传媒精品一区二区| 女人久久www免费人成看片| 国产免费视频播放在线视频| 久久久久久久久久人人人人人人| av卡一久久| 国产男女超爽视频在线观看| 亚洲欧美清纯卡通| 日本-黄色视频高清免费观看| 亚洲成人手机| 国产精品久久久久久久电影| 大香蕉97超碰在线| 欧美精品一区二区大全| 国产乱来视频区| 一区二区三区精品91| 热re99久久精品国产66热6| 国产精品偷伦视频观看了| 国产69精品久久久久777片| 在线观看免费日韩欧美大片 | 菩萨蛮人人尽说江南好唐韦庄| 国产欧美日韩一区二区三区在线 | 在线看a的网站| av女优亚洲男人天堂| 永久免费av网站大全| 国产精品99久久99久久久不卡 | 国产精品蜜桃在线观看| 日韩制服骚丝袜av| 丝袜喷水一区| 日本猛色少妇xxxxx猛交久久| 欧美97在线视频| 久久人人爽人人爽人人片va| 久热这里只有精品99| av天堂久久9| av专区在线播放| 美女国产视频在线观看| 97精品久久久久久久久久精品| 亚洲性久久影院| 亚洲精品色激情综合| 中文字幕免费在线视频6| 久久久亚洲精品成人影院| 男人爽女人下面视频在线观看| 精品国产露脸久久av麻豆| 免费观看的影片在线观看| 高清黄色对白视频在线免费看| 亚洲精品日本国产第一区| 天天躁夜夜躁狠狠久久av| 国产精品国产三级国产专区5o| 最新的欧美精品一区二区| 99热这里只有是精品在线观看| 男女啪啪激烈高潮av片| 超碰97精品在线观看| 五月伊人婷婷丁香| 国产亚洲av片在线观看秒播厂| 国产成人免费观看mmmm| 精品一品国产午夜福利视频| 欧美国产精品一级二级三级| 在线精品无人区一区二区三| 亚洲成人手机| 成年美女黄网站色视频大全免费 | 国产免费福利视频在线观看| 欧美日韩视频高清一区二区三区二| 欧美 日韩 精品 国产| 一区二区av电影网| 91成人精品电影| 国产精品人妻久久久影院| 精品视频人人做人人爽| 美女xxoo啪啪120秒动态图| 伊人久久精品亚洲午夜| 性色avwww在线观看| 午夜福利网站1000一区二区三区| 美女国产高潮福利片在线看| 这个男人来自地球电影免费观看 | 视频区图区小说| tube8黄色片| 日韩av免费高清视频| 欧美变态另类bdsm刘玥| 人人澡人人妻人| 日本黄色日本黄色录像| 亚洲高清免费不卡视频| 夜夜看夜夜爽夜夜摸| 男人爽女人下面视频在线观看| 欧美xxⅹ黑人| 韩国av在线不卡| 精品久久久久久久久亚洲| 国产日韩欧美亚洲二区| 亚洲内射少妇av| 国产毛片在线视频| 免费观看无遮挡的男女| 肉色欧美久久久久久久蜜桃| 国产成人aa在线观看| 亚洲色图 男人天堂 中文字幕 | 精品视频人人做人人爽| 丰满迷人的少妇在线观看| 国产精品.久久久| 春色校园在线视频观看| 久久久亚洲精品成人影院| 婷婷成人精品国产| 久久久久久伊人网av| 日韩一区二区视频免费看| 色视频在线一区二区三区| 人人妻人人澡人人爽人人夜夜| 人妻人人澡人人爽人人| 人成视频在线观看免费观看| 国产无遮挡羞羞视频在线观看| 精品国产乱码久久久久久小说| 久久精品国产a三级三级三级| 亚洲一级一片aⅴ在线观看| av专区在线播放| 视频中文字幕在线观看| 国产成人精品在线电影| 美女国产高潮福利片在线看| 亚洲精品日韩av片在线观看| 精品亚洲成a人片在线观看| 亚洲内射少妇av| 免费看光身美女| 亚洲欧美清纯卡通| 这个男人来自地球电影免费观看 | 日韩中文字幕视频在线看片| 亚洲国产精品999| 午夜福利视频精品| 亚洲国产精品一区二区三区在线| 亚洲精品成人av观看孕妇| 看非洲黑人一级黄片| 亚洲第一av免费看| 99九九线精品视频在线观看视频| 久久鲁丝午夜福利片| 我要看黄色一级片免费的| 亚洲精品自拍成人| 18禁动态无遮挡网站| 国产一区亚洲一区在线观看| 麻豆成人av视频| 欧美日韩国产mv在线观看视频| 三级国产精品欧美在线观看| 国产精品久久久久久精品古装| 日韩制服骚丝袜av| 80岁老熟妇乱子伦牲交| 卡戴珊不雅视频在线播放| 国产成人精品久久久久久| 欧美丝袜亚洲另类| 中国三级夫妇交换| 亚洲欧美成人精品一区二区| 99久久综合免费| 日韩av不卡免费在线播放| 99热这里只有精品一区| 少妇的逼好多水| 亚洲精品美女久久av网站| 春色校园在线视频观看| 日本与韩国留学比较| av黄色大香蕉| 久久午夜综合久久蜜桃| 国产av精品麻豆| 国产精品国产三级国产专区5o| 国产国语露脸激情在线看| 亚洲丝袜综合中文字幕| 啦啦啦视频在线资源免费观看| 全区人妻精品视频| 天堂俺去俺来也www色官网| av免费在线看不卡| 18+在线观看网站| 亚洲色图综合在线观看| 国产男女超爽视频在线观看| 久久国内精品自在自线图片| 亚洲av欧美aⅴ国产| av一本久久久久| 欧美 日韩 精品 国产| www.av在线官网国产| 一级毛片我不卡| 我的老师免费观看完整版| 九色亚洲精品在线播放| 亚洲精品久久久久久婷婷小说| 女性被躁到高潮视频| 国产黄色免费在线视频| 91精品国产九色| 国产极品粉嫩免费观看在线 | 91精品三级在线观看| 日日摸夜夜添夜夜爱| 亚洲欧美精品自产自拍| 日韩一区二区三区影片| 3wmmmm亚洲av在线观看| 国产高清三级在线| 男女高潮啪啪啪动态图| 91精品三级在线观看| 国产高清有码在线观看视频| 亚洲成人av在线免费| 国产精品久久久久成人av| 久久久亚洲精品成人影院| 成人二区视频| 日韩大片免费观看网站| 考比视频在线观看| 国产综合精华液| 免费高清在线观看日韩| 一本一本综合久久| 天天躁夜夜躁狠狠久久av| 免费看av在线观看网站| 一区二区av电影网| 中文字幕人妻熟人妻熟丝袜美| 99热国产这里只有精品6| 精品人妻一区二区三区麻豆| 黑人高潮一二区| 国产不卡av网站在线观看| 亚洲综合精品二区| 欧美人与性动交α欧美精品济南到 | 色吧在线观看| 99re6热这里在线精品视频| 亚洲国产色片| 成人国产麻豆网| 日韩在线高清观看一区二区三区| 99久久综合免费| 久久99精品国语久久久| 亚洲三级黄色毛片| 久久韩国三级中文字幕| 超色免费av| 大又大粗又爽又黄少妇毛片口| 亚洲,一卡二卡三卡| 国产永久视频网站| 国产亚洲av片在线观看秒播厂| 9色porny在线观看| 成人毛片60女人毛片免费| av线在线观看网站| 在线观看美女被高潮喷水网站| 国产有黄有色有爽视频| 女的被弄到高潮叫床怎么办| 伦理电影免费视频| 国产成人精品久久久久久| 亚洲精品乱码久久久v下载方式| 精品久久国产蜜桃| .国产精品久久| 中文乱码字字幕精品一区二区三区| 国产女主播在线喷水免费视频网站| 日本爱情动作片www.在线观看| av电影中文网址| 精品熟女少妇av免费看| 狂野欧美激情性bbbbbb| 日本与韩国留学比较| 18禁观看日本| 777米奇影视久久| 性高湖久久久久久久久免费观看| 国产片特级美女逼逼视频|