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

    土石混合體邊坡模型隨機構(gòu)建與穩(wěn)定性分析

    2021-09-28 09:31:24張宗堂楊期君張巨峰李文菁
    華南地質(zhì) 2021年2期
    關(guān)鍵詞:混合體土石塊石

    賀 勇, 張宗堂,楊期君,張巨峰,李文菁

    (1. 長沙市規(guī)劃勘測設(shè)計研究院,長沙 410000; 2. 湖南科技大學(xué)巖土工程穩(wěn)定控制與健康監(jiān)測湖南省重點實驗室,湖南 湘潭 411201; 3. 湖南科技大學(xué)資源環(huán)境與安全工程學(xué)院,湖南 湘潭 411201)

    土石混合體(soil rock mixture)是指由具有一定尺寸且強度較高的塊石、強度較低的土體及孔隙構(gòu)成的極端不均勻的松散巖土體,是一種在自然界廣泛存在的典型非均質(zhì)巖土材料[1-5];Medley也稱之為bimrock(block-in-matric rock),意思是“由較細(xì)結(jié)構(gòu)的粘結(jié)基質(zhì)與具有巖土工程意義的塊體組成”,其中“具有巖土工程意義的塊體”表示塊體與弱基質(zhì)之間存在足夠的力學(xué)對比,并且塊體的尺寸分布和體積比例會影響工程規(guī)模下的巖體性質(zhì)[6,7]。在巖土工程領(lǐng)域,土石混合體材料十分常見,例如:川藏公路林芝至八宿段全長427 km范圍內(nèi),地質(zhì)災(zāi)害的主要物質(zhì)是土石混合體[8];浙江省上(虞)三(門)高速公路K92-K95段的邊坡為土石混合體,滑坡治理費用高達(dá)億元[9];攀西地區(qū)發(fā)生的816個滑坡中,土石混合體邊坡就有500個,占61.3%[3,10]。

    考慮到這些土石混合體的巨大空間與力學(xué)變異性,以往學(xué)者們?yōu)榱撕喕嬎?,常常選擇使用較弱強度和變形特征好模擬的土壤來開展工作[11-12],忽略了塊體對整體土石混合體強度的貢獻(xiàn)。然而,正如許多文獻(xiàn)所報道的一樣,這種簡化可能導(dǎo)致較大的計算誤差。也有實踐表明,由于土石混合體本身的多向性及不連續(xù)性等特征,使用傳統(tǒng)的均質(zhì)巖土材料假定應(yīng)用于土石混合體巖土構(gòu)筑物的分析存在明顯不足[13]。在目前的研究中缺乏與自然界真實情況符合程度較高的土石混合體模型。

    在過去的幾十年里,人們進(jìn)行了大量的研究,以建立正確描述土石混合體特征的系統(tǒng)方法。學(xué)者們大都致力于選擇適當(dāng)?shù)膹姸群妥冃螀?shù),并進(jìn)行恰當(dāng)?shù)臄?shù)值模擬,以便在復(fù)雜地層中正確地開展工程活動。在土石混合體室內(nèi)試驗結(jié)果的基礎(chǔ)上,一些學(xué)者提出了初步的強度準(zhǔn)則:假定土石混合體是均勻各向同性的且具有等效的力學(xué)性能,則土石混合體的強度可根據(jù)其塊石含量和土體強度參數(shù)確定[14]。還有一些學(xué)者分析了在用確定性方法分析非均質(zhì)地層中邊坡穩(wěn)定性時塊石的影響[15],他們強調(diào)了這種土石混合體建模方法產(chǎn)生的特殊結(jié)果,與分析均質(zhì)土體時得到的結(jié)果有顯著不同。

    基于前人的研究成果,本文提出了一種土石混合體的隨機生成算法。首先采用數(shù)字圖像處理來獲取自然界的真實塊石輪廓,然后通過MATLAB編程來生成具有特定級配和塊石含量土石混合體邊坡模型。最后分別基于FLAC2D與OPTUM G2對其進(jìn)行穩(wěn)定性分析;為了提高結(jié)果的統(tǒng)計有效性,對于每種含石量的邊坡(含石量0%的邊坡除外)進(jìn)行了5次平行試驗,并對兩種方法的計算結(jié)果進(jìn)行了對比分析。

    1 土石混合體邊坡模型的構(gòu)建

    1.1 參數(shù)說明

    本文旨在探討含石量對土石混合體邊坡穩(wěn)定性的影響,為此,進(jìn)行了大量的二維邊坡穩(wěn)定性分析。對于簡單的邊坡模型,分別采用了有限差分法(FDM)和極限分析有限元法(FELA)對邊坡模型進(jìn)行分析。其中,所有邊坡模型的坡度為45°,邊坡尺寸如圖1所示,級配均一致(圖2),但塊體含量不同。分別計算了塊石含量為0%(均值土體)、20%、40%和60%的土石混合體邊坡模型。由于它們是二維模型,塊體的面積即被認(rèn)為等于塊石的體積。

    圖1 邊坡尺寸示意圖(單位:m)Fig. 1 Schematic diagram of slope size (unit: m)

    圖2 級配曲線Fig. 2 Grading curve

    為了考慮塊石在空間和尺寸上的變化,作者編寫了一個MATLAB程序,用于在邊坡區(qū)域中生成隨機位置和大小的塊石顆粒。對于每個給定的級配和含石量,共隨機構(gòu)建了5個不同的邊坡模型,即共進(jìn)行了5次穩(wěn)定性分析。設(shè)置5組平行試驗的目的是消除隨機誤差,以增加結(jié)果的可靠性。

    此外,還分析了含石量為0%的邊坡,以評估地質(zhì)勘探人員在僅根據(jù)土體的強度和變形特性進(jìn)行設(shè)計時可能造成的潛在誤差。表1顯示了分析中使用的塊石和土體參數(shù)[1-2,16-18],兩種材料的力學(xué)參數(shù)取值符合要求(Erock/Ematrix>2和tanφrock/tanφmatrix>2)[6,19]。假定它們遵循莫爾-庫侖破壞準(zhǔn)則,并具有理想彈塑性行為;根據(jù)Lindquist和Kalender準(zhǔn)則[20],假設(shè)塊石為均勻各向同性材料,并通過不同方法對其進(jìn)行了數(shù)值模擬??偣矊?2個土石混合體邊坡分別采用FLAC2D(有限差分法)和OPTUM G2(極限分析有限元)進(jìn)行了穩(wěn)定性分析,評估并比較了安全系數(shù)。

    表1 參數(shù)取值Table 1 Parameter values

    2 邊坡模型隨機構(gòu)建

    2.1 數(shù)字圖像處理獲取塊石顆粒

    通過對真實塊石顆粒進(jìn)行拍照并進(jìn)行數(shù)字圖像處理,來提取塊石顆粒的輪廓坐標(biāo)。采用數(shù)碼相機對顆粒進(jìn)行拍照,相機型號為:佳能600D,高級攝影系統(tǒng)(APS)框架和單鏡頭反射(SLR)數(shù)碼相機,1800萬像素APS-C CMOS傳感器(尺寸:22.3*14.9 mm)和EXPEED 2圖像處理芯片,1.6倍鏡頭轉(zhuǎn)換系數(shù);鏡頭型號為:EF-S 18-135 mm,67 mm鏡頭;得到的原始塊石照片如圖3a所示。首先,利用MATLAB中的rgb2gray函數(shù)對塊石顆粒的二維照片進(jìn)行灰度處理(圖3b);然后,利用MATLAB中的fspecial和imfilter函數(shù)對圖像進(jìn)行高斯模糊平滑處理(圖3c),降低了圖像的噪聲;使用MATLAB中的im2bw函數(shù)對照片進(jìn)行二值化處理(圖3d)后,可以看到二值化處理后顆粒區(qū)域為黑色,背景區(qū)域為白色;但是顆粒上仍然有一些部分是白色,這是因為顆粒上面的顏色分布不均勻,表面有少量白色區(qū)域造成的,不過這并不會影響輪廓的提取。最后,使用MATLAB中的bwboundaries函數(shù)搜索邊界像素來提取顆粒的輪廓(圖3e)。

    需要注意的是,輪廓提取后(圖3e)得到的顆粒坐標(biāo)是像素坐標(biāo),意味著一個顆粒的輪廓點有幾千上萬個,這樣的塊石顆粒直接用于數(shù)值計算是不現(xiàn)實的,因為輪廓點太多會造成計算十分耗時。因此,需要對塊石的輪廓點進(jìn)行簡化處理。常用的簡化方法是等角度隨機采樣(圖3e到圖3f),詳細(xì)過程如圖4所示。當(dāng)采樣角度為θ°時,粒子可簡化為N(N=360/θ)個點(在本研究中θ=10°,N=36),最終得到簡化后的顆粒輪廓如圖3f所示。等角度隨機采樣的詳細(xì)過程參見文獻(xiàn)[21]。

    圖3 數(shù)字圖像處理獲取塊石輪廓Fig. 3 Block stone outline acquired by digital image processing

    圖4 等角度采樣示意圖Fig. 4 Schematic diagram of equal angle sampling

    2.2 MATLAB編程構(gòu)建邊坡模型

    在考慮塊石存在的情況下,利用MATLAB編寫程序進(jìn)行蒙特卡羅模擬,根據(jù)給定的塊石含量和級配曲線,在邊坡區(qū)域內(nèi)生成任意直徑和位置的真實塊石顆粒。對于每個不同的含石量(僅20%、40%和60%,因為0%代表均值土體),隨機構(gòu)建了5個土石混合體邊坡模型,即為了實現(xiàn)結(jié)果的統(tǒng)計有效性,進(jìn)行了5次穩(wěn)定性分析。

    當(dāng)計算得到的含石量(計算為所有塊體的面積總和除以邊坡模型的面積)對應(yīng)于輸入中定義的含石量時,使用“while”循環(huán)停止提取,允許誤差為±2.5%。

    為了最大限度地提高代碼性能,加快計算速度,塊石的放置是根據(jù)塊石粒徑的大小從大到小進(jìn)行的。MATLAB代碼滿足的其他要求是:(1)塊石不能相互穿透,否則就沒有物理意義;(2)塊石不能與它們所在的斜坡模型的外部邊界相交。實際上,任何交叉點都會導(dǎo)致塊體體積的部分損失,從而導(dǎo)致對塊石含量的低估,這將與規(guī)定的塊石含量不符。為了達(dá)到這些目的,將兩個塊石之間以及塊石與外部斜坡邊界之間的最小距離設(shè)置為5 cm。最終構(gòu)建的土石混合體模型如圖5所示(含石量40%,級配如圖2所示)。

    圖5 含石量為40%的土石混合體邊坡模型示例Fig. 5 Example of soil rock mixture slope model with 40% stone content

    MATLAB代碼的主要輸出結(jié)果是一個文本文件,其中包含每個塊石的輪廓坐標(biāo)。為了導(dǎo)入FLAC2D和OPTUM G2,所有MATLAB輸出文件都轉(zhuǎn)換成腳本文件,以便在AutoCAD軟件中可視化,然后以DXF格式保存。

    3 土石混合體邊坡模型穩(wěn)定性分析

    分別采用有限差分法(FDM)與極限分析有限元法(FELA),對含石量為0%(均值土體)、20%、40%和60%的土石混合體進(jìn)行穩(wěn)定性分析。除均值土體外,每個含石量隨機構(gòu)建了5個邊坡進(jìn)行分析,共16個邊坡,即采用兩種方法一共進(jìn)行了32次穩(wěn)定性分析。

    3.1 有限差分法(FDM)分析

    基于軟件FLAC2D對土石混合體邊坡進(jìn)行有限差分法分析。首先將前文中得到的DXF格式的邊坡模型導(dǎo)入MIDAS GTS中進(jìn)行網(wǎng)格劃分,然后編寫一個MIDAS2FLAC2D的接口程序?qū)⒛P蛯?dǎo)入到FLAC2D中。分別對不同含石量的邊坡采用FLAC2D內(nèi)置的有限差分強度折減法進(jìn)行穩(wěn)定性分析,其計算結(jié)果如圖6所示,計算所得的安全系數(shù)如表2和圖7所示。

    由圖6可以看到,均質(zhì)土體的塑性區(qū)為一條類圓弧滑動帶,而土石混合體的塑性區(qū)則較為復(fù)雜。具體來說,滑動帶沿著塊石周圍進(jìn)行擴展,呈現(xiàn)出“繞石效應(yīng)”。這是由于塊石的強度相對土體較大,滑動帶無法貫穿塊石,因此只能繞過塊石進(jìn)行擴展。同時可以看到,隨著含石量的增大,邊坡的滑動帶變得更為破碎和曲折,這與徐文杰等[1]的報導(dǎo)一致。從安全系數(shù)計算結(jié)果上來看,雖然在相同含石量下,不同組邊坡計算所得的安全系數(shù)存在一定差異,但是總體上,邊坡的安全系數(shù)隨著含石量的增大而增大,這一點可以在安全系數(shù)平均值上(表2和圖7)體現(xiàn)出來。同時可以發(fā)現(xiàn),隨著含石量的增大,5組邊坡的安全系數(shù)的標(biāo)準(zhǔn)差也隨之增大,說明含石量越大,塊石分布的隨機性對邊坡安全系數(shù)的影響越大。

    圖6 有限差分法(其中一組平行試驗)分析結(jié)果Fig. 6 Analysis results of finite difference method (one group of parallel tests)

    圖7 邊坡的安全系數(shù)平均值Fig. 7 Average safety factor of slope

    表2 有限差分法安全系數(shù)計算結(jié)果Table 2 Calculation results of safety factor by finite difference method

    3.2 極限分析有限元(FELA)分析

    基于軟件OPTUM G2對土石混合體邊坡及逆行極限分析有限元分析(取上限解)。將前文中得到的DXF文件導(dǎo)入G2中進(jìn)行網(wǎng)格劃分,然后賦予材料參數(shù)(表1)與邊界條件(G2中的“標(biāo)準(zhǔn)邊界條件”)對不同含石量的土石混合體邊坡進(jìn)行穩(wěn)定性分析,其計算結(jié)果如圖8所示,計算所得的安全系數(shù)如表3和圖7所示。

    由圖8可以看到,F(xiàn)ELA分析得到的邊坡的破壞模式與圖6中(FDM分析)邊坡的破壞模式基本一致,即均質(zhì)土體的塑性區(qū)為一條類圓弧滑動帶,而土石混合體的塑性區(qū)則較為復(fù)雜,且隨著含石量的增大,土石混合體邊坡的滑動帶變得更加破碎。由表3和圖7可以看出,邊坡的安全系數(shù)隨含石量的增大而增大,且含石量越大,5組邊坡的安全系數(shù)的標(biāo)準(zhǔn)差也越大。

    表3 極限分析有限元法安全系數(shù)計算結(jié)果Table 3 Calculation results of safety factor of limit analysis finite element method

    圖8 極限分析有限元法(其中一組平行試驗)分析結(jié)果Fig. 8 Analysis results of limit analysis finite element method (one group of parallel tests)

    3.3 FDM與FELA的比較

    由上述兩種計算方式對比可以看到,采用不同方法進(jìn)行分析的土石混合體邊坡的破壞模式相差不大,但是安全系數(shù)存在一定的差異。具體來說:FELA方法計算得到的邊坡安全系數(shù)的平均值會略大于FDM方法計算得到的安全系數(shù)的平均值,且差值隨著含石量的增大而增大(圖7)。同時,對于同一個含石量,F(xiàn)ELA計算得到的邊坡安全系數(shù)的標(biāo)準(zhǔn)差要大于FDM;這是由于在采用FELA方法進(jìn)行分析時,我們采用的是極限分析有限元的上限解。所以,F(xiàn)ELA計算出的結(jié)果總是偏安全的,即計算出的邊坡安全系數(shù)偏大。

    4 結(jié)論

    本文通過數(shù)字圖像處理獲取了自然界真實顆粒的數(shù)字化模型,然后基于MATLAB進(jìn)行特定級配與含石量的土石混合體邊坡模型的構(gòu)建,最后分別采用FDM與FELA方法對不同含石量的土石混合體邊坡進(jìn)行了穩(wěn)定性分析。結(jié)果表明:

    (1)采用本文提出的土石混合體隨機生成方法,將塊石和土體分開考慮進(jìn)行計算,可有效評估土石混合體邊坡穩(wěn)定性。

    (2)含石量對邊坡穩(wěn)定性影響較大,具體表現(xiàn)為:隨著含石量的增大,土石混合體邊坡的滑動帶變得更為破碎和曲折;邊坡的安全系數(shù)與5組邊坡安全系數(shù)的標(biāo)準(zhǔn)差隨著含石量的增大而增大。

    (3)對于土石混合體邊坡,F(xiàn)ELA方法計算得到的邊坡安全系數(shù)的平均值大于FDM方法計算的平均值。

    猜你喜歡
    混合體土石塊石
    強風(fēng)化地層上土石圍堰的截流結(jié)構(gòu)設(shè)計及施工
    不同粒徑組合塊石群水下漂移數(shù)值模擬
    沉管隧道DCM區(qū)地基塊石拋填振密施工技術(shù)
    地基土中基床塊石沉降變形研究
    水道港口(2020年5期)2020-12-09 07:37:46
    基于蒙特卡洛隨機采樣的土石混合體數(shù)值模型構(gòu)建方法
    市政道路高填方區(qū)域土石混填技術(shù)研究
    壓實度對筑堤用土石混合料力學(xué)特性影響研究
    治淮(2018年1期)2018-02-09 09:15:01
    不同水環(huán)境下土石混填地基水平推剪試驗研究
    克制的宣泄——鋼琴家普列特涅夫:浪漫主義、現(xiàn)代審美和個人風(fēng)格的混合體
    人民音樂(2016年1期)2016-11-07 10:02:53
    獨領(lǐng)風(fēng)騷的夏日
    亚洲精品国产区一区二| 成在线人永久免费视频| 99riav亚洲国产免费| 日韩欧美国产一区二区入口| 久热爱精品视频在线9| 色老头精品视频在线观看| 黄色成人免费大全| 在线观看免费视频日本深夜| 亚洲五月色婷婷综合| 精品国产超薄肉色丝袜足j| 丝袜在线中文字幕| 欧美日本亚洲视频在线播放| 久久狼人影院| 日韩欧美一区视频在线观看| 久久国产精品男人的天堂亚洲| 久久精品国产清高在天天线| 色综合站精品国产| 一级片免费观看大全| 国产成人精品无人区| 亚洲五月天丁香| 久久人人精品亚洲av| 国产伦一二天堂av在线观看| 精品国产一区二区久久| 久久久精品欧美日韩精品| 久久精品人人爽人人爽视色| 最好的美女福利视频网| 巨乳人妻的诱惑在线观看| 国产蜜桃级精品一区二区三区| 99精品欧美一区二区三区四区| 色综合站精品国产| 男女下面插进去视频免费观看| 电影成人av| 久久亚洲真实| 欧美日韩一级在线毛片| 可以在线观看的亚洲视频| 两人在一起打扑克的视频| 欧美另类亚洲清纯唯美| 国产精品一区二区免费欧美| 国产成人啪精品午夜网站| 91精品国产国语对白视频| 无限看片的www在线观看| 久久欧美精品欧美久久欧美| 香蕉丝袜av| x7x7x7水蜜桃| 中文字幕精品免费在线观看视频| 人人妻,人人澡人人爽秒播| 欧美黑人欧美精品刺激| 啪啪无遮挡十八禁网站| 大陆偷拍与自拍| 美女国产高潮福利片在线看| 韩国av一区二区三区四区| 亚洲av成人一区二区三| 亚洲av五月六月丁香网| 激情视频va一区二区三区| 午夜影院日韩av| 亚洲片人在线观看| 久久精品aⅴ一区二区三区四区| 日韩精品免费视频一区二区三区| 免费在线观看完整版高清| 香蕉丝袜av| www.www免费av| 一进一出抽搐动态| 欧美另类亚洲清纯唯美| 啦啦啦韩国在线观看视频| 亚洲国产欧美网| 国产乱人伦免费视频| 日韩免费av在线播放| 美女免费视频网站| 日本精品一区二区三区蜜桃| 日日干狠狠操夜夜爽| 青草久久国产| 最新在线观看一区二区三区| 9色porny在线观看| 国产熟女xx| 日本五十路高清| 国产精品自产拍在线观看55亚洲| 自拍欧美九色日韩亚洲蝌蚪91| 欧美性长视频在线观看| 黄色女人牲交| 女警被强在线播放| 精品高清国产在线一区| 国产精品久久久av美女十八| 日本a在线网址| 人妻久久中文字幕网| 岛国在线观看网站| 丁香欧美五月| 免费看十八禁软件| 黑人欧美特级aaaaaa片| 亚洲成人免费电影在线观看| 亚洲久久久国产精品| 一区在线观看完整版| 日日夜夜操网爽| 午夜精品国产一区二区电影| 天堂√8在线中文| 中文字幕人妻熟女乱码| av电影中文网址| 亚洲成人精品中文字幕电影| 日本a在线网址| 日韩欧美在线二视频| 亚洲专区国产一区二区| 亚洲国产看品久久| 久久精品影院6| 久久久国产精品麻豆| 日韩欧美三级三区| 欧美激情极品国产一区二区三区| 大型av网站在线播放| 中文字幕人妻熟女乱码| √禁漫天堂资源中文www| 久久久久久久久中文| 久久国产精品影院| 侵犯人妻中文字幕一二三四区| 久久欧美精品欧美久久欧美| 亚洲久久久国产精品| 中国美女看黄片| 桃色一区二区三区在线观看| 9191精品国产免费久久| 日韩免费av在线播放| 亚洲精品粉嫩美女一区| 又大又爽又粗| 色av中文字幕| 国产亚洲精品综合一区在线观看 | 国产一区二区三区在线臀色熟女| 亚洲狠狠婷婷综合久久图片| 色播亚洲综合网| 免费在线观看亚洲国产| 黑人欧美特级aaaaaa片| 老司机在亚洲福利影院| 91成人精品电影| 亚洲精品美女久久av网站| 波多野结衣一区麻豆| 久久久久精品国产欧美久久久| 国产伦一二天堂av在线观看| 波多野结衣高清无吗| 岛国视频午夜一区免费看| 色精品久久人妻99蜜桃| 大码成人一级视频| 韩国av一区二区三区四区| 亚洲专区中文字幕在线| 午夜免费观看网址| 亚洲五月天丁香| 成人亚洲精品av一区二区| 亚洲一区中文字幕在线| 天天躁狠狠躁夜夜躁狠狠躁| 国产精品九九99| 欧美成狂野欧美在线观看| 国产单亲对白刺激| 老熟妇乱子伦视频在线观看| 宅男免费午夜| 亚洲免费av在线视频| 亚洲午夜精品一区,二区,三区| 久久婷婷人人爽人人干人人爱 | 国产一区二区三区视频了| 亚洲视频免费观看视频| 精品一区二区三区四区五区乱码| 老司机深夜福利视频在线观看| 精品国产美女av久久久久小说| 后天国语完整版免费观看| 亚洲中文日韩欧美视频| 国产欧美日韩一区二区三| 黄片大片在线免费观看| 国产伦一二天堂av在线观看| 午夜福利一区二区在线看| 99在线视频只有这里精品首页| 色播在线永久视频| 久久中文看片网| bbb黄色大片| 精品高清国产在线一区| 亚洲精品在线美女| 少妇熟女aⅴ在线视频| 精品欧美一区二区三区在线| 窝窝影院91人妻| 每晚都被弄得嗷嗷叫到高潮| 欧美黑人欧美精品刺激| 国产精华一区二区三区| 欧美另类亚洲清纯唯美| 在线十欧美十亚洲十日本专区| 国产三级在线视频| 国产精品精品国产色婷婷| 国产国语露脸激情在线看| 淫妇啪啪啪对白视频| 亚洲国产高清在线一区二区三 | 在线观看免费视频日本深夜| 少妇被粗大的猛进出69影院| 色播亚洲综合网| 亚洲第一电影网av| 一边摸一边做爽爽视频免费| 日韩精品免费视频一区二区三区| 亚洲五月天丁香| 亚洲国产精品久久男人天堂| 丁香六月欧美| 久久久久久久久久久久大奶| 国产91精品成人一区二区三区| 久久久久久久午夜电影| 日本免费a在线| 亚洲免费av在线视频| 美女国产高潮福利片在线看| 99热只有精品国产| 亚洲一码二码三码区别大吗| 久久久久国内视频| 大码成人一级视频| 一级毛片女人18水好多| 一级a爱视频在线免费观看| 美女午夜性视频免费| 狠狠狠狠99中文字幕| 日本vs欧美在线观看视频| 成人欧美大片| 男人舔女人的私密视频| 成人免费观看视频高清| 久久香蕉激情| 亚洲精品av麻豆狂野| 亚洲专区中文字幕在线| 美女高潮到喷水免费观看| 国产成人一区二区三区免费视频网站| 色综合婷婷激情| 精品人妻1区二区| 视频区欧美日本亚洲| 亚洲人成电影免费在线| 男男h啪啪无遮挡| 午夜免费成人在线视频| 国产av精品麻豆| 97超级碰碰碰精品色视频在线观看| 99久久99久久久精品蜜桃| 人人澡人人妻人| 午夜福利成人在线免费观看| av福利片在线| 99国产精品99久久久久| 久久久久亚洲av毛片大全| 国产午夜精品久久久久久| 波多野结衣一区麻豆| 真人做人爱边吃奶动态| 亚洲男人天堂网一区| 麻豆av在线久日| 久久热在线av| 亚洲电影在线观看av| 亚洲精品美女久久av网站| 在线观看免费午夜福利视频| 啦啦啦观看免费观看视频高清 | 亚洲中文av在线| 无人区码免费观看不卡| 91九色精品人成在线观看| 大型黄色视频在线免费观看| 嫩草影院精品99| 嫁个100分男人电影在线观看| 两个人看的免费小视频| avwww免费| 国产麻豆成人av免费视频| 精品少妇一区二区三区视频日本电影| 性色av乱码一区二区三区2| 我的亚洲天堂| 国产精品久久久人人做人人爽| 操美女的视频在线观看| 999久久久国产精品视频| 亚洲av第一区精品v没综合| 99精品在免费线老司机午夜| 色在线成人网| 女警被强在线播放| 每晚都被弄得嗷嗷叫到高潮| 1024视频免费在线观看| 757午夜福利合集在线观看| 变态另类丝袜制服| 国产精品一区二区在线不卡| 欧美在线黄色| 国产亚洲欧美98| 国产99白浆流出| 久久亚洲真实| 久久久久国产精品人妻aⅴ院| 久久久久久大精品| 人妻丰满熟妇av一区二区三区| 欧美在线黄色| 精品久久蜜臀av无| 亚洲精品一区av在线观看| 国产熟女午夜一区二区三区| 精品不卡国产一区二区三区| 欧美成人性av电影在线观看| 99国产精品99久久久久| 19禁男女啪啪无遮挡网站| 久久午夜综合久久蜜桃| av欧美777| 国产麻豆69| 日本黄色视频三级网站网址| 狠狠狠狠99中文字幕| 亚洲av片天天在线观看| 亚洲精品粉嫩美女一区| 精品国内亚洲2022精品成人| 久久精品国产亚洲av香蕉五月| 性欧美人与动物交配| 亚洲成人久久性| 国产精品久久视频播放| 黄色成人免费大全| 亚洲精华国产精华精| 日韩国内少妇激情av| 19禁男女啪啪无遮挡网站| 国产精品一区二区在线不卡| 男女下面进入的视频免费午夜 | 国产精品永久免费网站| 亚洲精品久久成人aⅴ小说| 首页视频小说图片口味搜索| 欧美一区二区精品小视频在线| 黄片大片在线免费观看| 黄网站色视频无遮挡免费观看| av网站免费在线观看视频| 中文字幕精品免费在线观看视频| videosex国产| 一级毛片精品| 久久精品影院6| 制服人妻中文乱码| 日韩欧美一区二区三区在线观看| 一进一出抽搐动态| 久久人人97超碰香蕉20202| a在线观看视频网站| 在线十欧美十亚洲十日本专区| 天天添夜夜摸| 午夜日韩欧美国产| 国产精品自产拍在线观看55亚洲| 中文字幕av电影在线播放| 日本精品一区二区三区蜜桃| 午夜久久久久精精品| 在线观看免费视频网站a站| 成年女人毛片免费观看观看9| 97碰自拍视频| x7x7x7水蜜桃| 国产三级在线视频| 国产精品亚洲av一区麻豆| 日韩免费av在线播放| 麻豆国产av国片精品| 一边摸一边做爽爽视频免费| 色精品久久人妻99蜜桃| 国产亚洲av高清不卡| 国产精品日韩av在线免费观看 | 1024视频免费在线观看| 99国产精品一区二区蜜桃av| 在线观看日韩欧美| 国产伦一二天堂av在线观看| 免费av毛片视频| 欧美色视频一区免费| 啦啦啦韩国在线观看视频| 国产亚洲av高清不卡| 男女床上黄色一级片免费看| 免费不卡黄色视频| 精品不卡国产一区二区三区| 久久久国产精品麻豆| svipshipincom国产片| 亚洲av成人不卡在线观看播放网| 国产成人精品无人区| 精品久久久久久久人妻蜜臀av | 美女免费视频网站| 欧美成人一区二区免费高清观看 | 久久精品人人爽人人爽视色| 精品乱码久久久久久99久播| 18禁观看日本| 欧美av亚洲av综合av国产av| 老汉色av国产亚洲站长工具| 91国产中文字幕| 婷婷丁香在线五月| 亚洲专区字幕在线| 午夜两性在线视频| 最新在线观看一区二区三区| x7x7x7水蜜桃| 最近最新中文字幕大全电影3 | 欧美日韩福利视频一区二区| 99国产精品99久久久久| 99re在线观看精品视频| 亚洲av电影在线进入| 一区二区三区国产精品乱码| 欧美绝顶高潮抽搐喷水| 巨乳人妻的诱惑在线观看| 制服丝袜大香蕉在线| 婷婷丁香在线五月| 在线观看免费视频日本深夜| 变态另类丝袜制服| 精品国产国语对白av| 天天躁狠狠躁夜夜躁狠狠躁| 午夜福利在线观看吧| 黄色a级毛片大全视频| 一区二区三区激情视频| 精品久久久久久久毛片微露脸| 亚洲人成伊人成综合网2020| 欧美乱色亚洲激情| 亚洲片人在线观看| 成年人黄色毛片网站| 日日爽夜夜爽网站| 婷婷精品国产亚洲av在线| 免费看美女性在线毛片视频| 午夜老司机福利片| 国产熟女xx| 日本欧美视频一区| 老鸭窝网址在线观看| 夜夜躁狠狠躁天天躁| 亚洲国产欧美日韩在线播放| 亚洲狠狠婷婷综合久久图片| 亚洲国产精品久久男人天堂| 老司机靠b影院| 啪啪无遮挡十八禁网站| 国产日韩一区二区三区精品不卡| 久久精品91蜜桃| 久久久水蜜桃国产精品网| 免费在线观看亚洲国产| 女性被躁到高潮视频| 他把我摸到了高潮在线观看| a级毛片在线看网站| 黑人欧美特级aaaaaa片| 黑人巨大精品欧美一区二区蜜桃| 日韩欧美免费精品| 国产精品 欧美亚洲| 亚洲精品中文字幕在线视频| 久久精品国产清高在天天线| 国产极品粉嫩免费观看在线| 亚洲七黄色美女视频| 国产成人欧美在线观看| 亚洲精品一卡2卡三卡4卡5卡| 人人妻人人爽人人添夜夜欢视频| 亚洲人成电影免费在线| 亚洲va日本ⅴa欧美va伊人久久| 日韩欧美一区视频在线观看| 成人手机av| 真人做人爱边吃奶动态| 中文字幕人成人乱码亚洲影| 亚洲七黄色美女视频| 精品久久久久久久久久免费视频| 亚洲全国av大片| 人成视频在线观看免费观看| 国产精品一区二区免费欧美| 日韩欧美国产一区二区入口| 亚洲aⅴ乱码一区二区在线播放 | 午夜福利欧美成人| 91国产中文字幕| 女人爽到高潮嗷嗷叫在线视频| 在线观看免费视频网站a站| 制服诱惑二区| 老熟妇乱子伦视频在线观看| 亚洲av成人av| 亚洲免费av在线视频| 欧美国产精品va在线观看不卡| 国产一卡二卡三卡精品| 美女 人体艺术 gogo| 岛国在线观看网站| 国产欧美日韩精品亚洲av| 国产成人欧美在线观看| 成人特级黄色片久久久久久久| 亚洲狠狠婷婷综合久久图片| 国产成人啪精品午夜网站| 老熟妇仑乱视频hdxx| 黄色视频,在线免费观看| 午夜久久久久精精品| 亚洲av成人不卡在线观看播放网| 日日干狠狠操夜夜爽| 9色porny在线观看| 国产xxxxx性猛交| 一区二区三区精品91| 亚洲一区二区三区不卡视频| 久久天堂一区二区三区四区| 可以在线观看毛片的网站| 91字幕亚洲| av天堂在线播放| 69av精品久久久久久| 国产99久久九九免费精品| 亚洲国产精品999在线| 男人的好看免费观看在线视频 | 午夜影院日韩av| 色精品久久人妻99蜜桃| 久久精品国产综合久久久| 一进一出好大好爽视频| 别揉我奶头~嗯~啊~动态视频| 岛国在线观看网站| 90打野战视频偷拍视频| 亚洲精品国产色婷婷电影| 妹子高潮喷水视频| 国产精品香港三级国产av潘金莲| 成人手机av| 激情在线观看视频在线高清| 成人免费观看视频高清| 成熟少妇高潮喷水视频| 成人av一区二区三区在线看| a在线观看视频网站| 国产免费男女视频| 国产亚洲欧美在线一区二区| 国产精品一区二区在线不卡| 国产亚洲精品第一综合不卡| 麻豆国产av国片精品| 欧美 亚洲 国产 日韩一| 日韩欧美一区二区三区在线观看| av福利片在线| 黑人巨大精品欧美一区二区mp4| 一夜夜www| 国产精品野战在线观看| 日韩欧美一区二区三区在线观看| av片东京热男人的天堂| 亚洲精品中文字幕在线视频| 99国产精品一区二区蜜桃av| 黑人操中国人逼视频| 国产精品,欧美在线| 亚洲成a人片在线一区二区| 精品电影一区二区在线| 制服诱惑二区| 欧美日韩精品网址| 免费看美女性在线毛片视频| 丁香六月欧美| 欧美黄色片欧美黄色片| 在线播放国产精品三级| 免费看a级黄色片| 热99re8久久精品国产| 久久精品亚洲熟妇少妇任你| 91成年电影在线观看| 国产精品一区二区在线不卡| 男人舔女人下体高潮全视频| 夜夜夜夜夜久久久久| 亚洲专区中文字幕在线| 丝袜在线中文字幕| 视频在线观看一区二区三区| 搞女人的毛片| 淫秽高清视频在线观看| 人妻丰满熟妇av一区二区三区| 男人舔女人下体高潮全视频| 91麻豆av在线| 婷婷精品国产亚洲av在线| 99香蕉大伊视频| 身体一侧抽搐| 久久人妻熟女aⅴ| 级片在线观看| 97人妻精品一区二区三区麻豆 | 久久久国产欧美日韩av| 久9热在线精品视频| 亚洲国产精品合色在线| 国产伦人伦偷精品视频| 亚洲色图综合在线观看| 久热这里只有精品99| 亚洲第一欧美日韩一区二区三区| 久久久国产欧美日韩av| 女人精品久久久久毛片| 国产男靠女视频免费网站| 美女 人体艺术 gogo| 大型黄色视频在线免费观看| 午夜老司机福利片| 日韩欧美国产一区二区入口| 国产麻豆成人av免费视频| 韩国精品一区二区三区| 在线观看日韩欧美| 久久婷婷成人综合色麻豆| 黄片大片在线免费观看| 九色亚洲精品在线播放| 18禁美女被吸乳视频| 国产伦一二天堂av在线观看| 久久久精品欧美日韩精品| 国产熟女xx| 午夜亚洲福利在线播放| 99国产精品一区二区三区| 亚洲七黄色美女视频| 两性午夜刺激爽爽歪歪视频在线观看 | 男女午夜视频在线观看| 国产精品秋霞免费鲁丝片| ponron亚洲| 在线观看免费午夜福利视频| 成人特级黄色片久久久久久久| 禁无遮挡网站| 少妇粗大呻吟视频| 女性被躁到高潮视频| 免费女性裸体啪啪无遮挡网站| 咕卡用的链子| 欧美日韩乱码在线| 国产成+人综合+亚洲专区| 欧美中文日本在线观看视频| 午夜福利高清视频| 免费在线观看黄色视频的| 亚洲熟女毛片儿| 这个男人来自地球电影免费观看| 欧美国产精品va在线观看不卡| 亚洲国产精品成人综合色| 激情在线观看视频在线高清| 亚洲国产精品成人综合色| 美女国产高潮福利片在线看| 国产精品av久久久久免费| 麻豆一二三区av精品| 亚洲精品国产精品久久久不卡| 欧美最黄视频在线播放免费| 国产一区二区激情短视频| 亚洲 欧美一区二区三区| 午夜福利免费观看在线| 在线国产一区二区在线| 一个人免费在线观看的高清视频| 亚洲av成人一区二区三| 19禁男女啪啪无遮挡网站| 亚洲av成人一区二区三| 叶爱在线成人免费视频播放| 女人爽到高潮嗷嗷叫在线视频| 日本免费a在线| www日本在线高清视频| 欧美成人性av电影在线观看| 成人国产综合亚洲| 亚洲黑人精品在线| 老司机在亚洲福利影院| 久久精品成人免费网站| 国产午夜福利久久久久久| 一区在线观看完整版| 亚洲中文字幕日韩| 69av精品久久久久久| 欧美性长视频在线观看| 夜夜爽天天搞| 一边摸一边抽搐一进一出视频| 亚洲av第一区精品v没综合| 国产精品国产高清国产av| 亚洲熟妇中文字幕五十中出| 午夜福利成人在线免费观看| 国产免费av片在线观看野外av| 欧美一级毛片孕妇| 最好的美女福利视频网| 免费在线观看黄色视频的| 国产男靠女视频免费网站| 国产成人欧美在线观看| 婷婷六月久久综合丁香| 黑人欧美特级aaaaaa片| 国产黄a三级三级三级人| 精品高清国产在线一区| 久久久精品欧美日韩精品| 亚洲最大成人中文| 久久亚洲真实|