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

    基于ABAQUS土坡降雨入滲應(yīng)力應(yīng)變及穩(wěn)定性研究

    2017-11-07 06:03:52
    水資源開(kāi)發(fā)與管理 2017年10期

    (新疆瑪納斯河流域管理局水利管理中心第三水利管理站,新疆 石河子 832000)

    基于ABAQUS土坡降雨入滲應(yīng)力應(yīng)變及穩(wěn)定性研究

    獨(dú)敏

    (新疆瑪納斯河流域管理局水利管理中心第三水利管理站,新疆 石河子 832000)

    降雨會(huì)影響土體的力學(xué)參數(shù),通過(guò)ABAQUS對(duì)雨水入滲的瞬態(tài)滲流對(duì)土坡應(yīng)力、應(yīng)變的影響進(jìn)行研究,以一理想邊坡作為案例進(jìn)行建模與計(jì)算,不考慮降雨積水的影響,認(rèn)為靜水壓力隨深度線性增加,建立CPE4P網(wǎng)格進(jìn)行計(jì)算。后處理中,對(duì)其滲透率、孔壓、橫縱向位移、偏應(yīng)力、塑性應(yīng)變以及應(yīng)力路徑等作了分析,得出一系列降雨入滲作用下應(yīng)力、應(yīng)變及穩(wěn)定性的相關(guān)結(jié)論。

    降雨;滲流;ABAQUS;失穩(wěn);邊坡;塑性區(qū)

    土體的抗剪強(qiáng)度等力學(xué)參數(shù)會(huì)隨著降雨的作用而降低[1],空隙水的壓力會(huì)隨著地下水位的抬高而升高[2- 3],造成穩(wěn)定性的下降[4]。除此之外,高強(qiáng)度、長(zhǎng)時(shí)間的降雨會(huì)導(dǎo)致暫態(tài)飽和區(qū)的出現(xiàn),以此,會(huì)導(dǎo)致空隙水壓力的升高?;诖?,很有必要對(duì)雨水入滲情況下瞬態(tài)滲流影響下的土坡應(yīng)力、應(yīng)變等力學(xué)特性進(jìn)行研究。

    降雨入滲作為—種典型的非飽和流固耦合現(xiàn)象,目前已有較多學(xué)者對(duì)降雨入滲進(jìn)行了研究,取得了較多成果[5- 10]。但研究領(lǐng)域多在降雨入滲的影響因素、滑坡、初始含水率、穩(wěn)定性等層面上的分析,而基于多參數(shù)的、多指標(biāo)的綜合分析研究較少。本文從數(shù)值模擬角度出發(fā),全面對(duì)滲透率、孔壓、橫縱向位移、偏應(yīng)力、塑性應(yīng)變以及應(yīng)力路徑等作了分析,以此更好地為降雨入滲對(duì)土坡的穩(wěn)定性研究做參考與支持。

    1 基于ABAQUS的降雨入滲計(jì)算力學(xué)原理及邊界條件

    1.1 計(jì)算原理

    土體材料的滲透系數(shù)與基質(zhì)吸力關(guān)系為

    kw=awkws/{[aw+bw(ua-uw)]cw}

    (1)

    式中kw—— 土體飽和時(shí)的滲透系數(shù),取5.0×10-6m/s(0.018m/h);

    uw、ua—— 土體中的水壓力、氣壓,因坡面與大氣接觸的原因,在式中取為0;

    aw、bw、cw—— 材料系數(shù),本例中分別取為1000、0.01、1.7。

    飽和度隨基質(zhì)吸力的關(guān)系式為

    Sr=Si+(Sn-Si)as/{[bw(ua-uw)]cw}

    (2)

    式中Sr—— 飽和度;

    Si—— 殘余飽和度,可取0.08;

    Sn—— 最大飽和度,取1;

    as、bs、cs—— 材料系數(shù),分別取為1、5.0×10-5和3.5。

    值得注意的是,ABAQUS中,有效應(yīng)力的定義為σ′=σ-(χuw+(1-χ)ua),χ簡(jiǎn)單的取為飽和度Sr,將ua賦值為0,采用基于有效應(yīng)力的本構(gòu)關(guān)系,該本構(gòu)重點(diǎn)體現(xiàn)了基質(zhì)吸力對(duì)模型的影響,這與反映黏聚力的強(qiáng)度本構(gòu)具有較大差異。

    1.2 邊界條件

    降雨入滲的物理過(guò)程十分復(fù)雜,Larson與Mein研究發(fā)現(xiàn),通過(guò)土壤允許的入滲容量fp,降雨強(qiáng)度q以及水力傳導(dǎo)系數(shù)kws參數(shù)的分析,可以揭示出降雨入滲行為如下:

    a.q

    b.fp>q>kws,雨水在該情況下全部入滲,隨著入滲度的不斷增加,fp不斷減小,但該強(qiáng)度下,降雨強(qiáng)度還未超過(guò)允許容量,因此,入滲率并不會(huì)出現(xiàn)降低的情況,其仍具有較高的數(shù)值。該情況時(shí),坡面為流量邊界。

    c.q>fp,此時(shí),降雨強(qiáng)度已經(jīng)超過(guò)了土壤的入滲容量,因此會(huì)出現(xiàn)部分降雨不入滲而形成地表徑流的情況。該情況中坡面上的土體為飽和狀態(tài)的土體,當(dāng)達(dá)到入滲容量后,降雨將不斷下降。

    本文考慮所述第二種情況,即fp>q>kws情況,文中的降雨邊界函數(shù)用降雨強(qiáng)度(單位流通量)來(lái)表示。同時(shí),不考慮地表積水的情況,使分析區(qū)域頂面均受到降雨作用,設(shè)置入滲強(qiáng)度和持續(xù)時(shí)間分別為20mm/h、66h。

    2 降雨入滲模型的建立與求解

    降雨入滲的模型來(lái)自于新疆瑪納斯河邊坡。新疆瑪納斯河位于新疆維吾爾自治區(qū)準(zhǔn)噶爾盆地南部,源出天山北麓北流注入瑪納斯湖,長(zhǎng)約450km,上游水急多峽谷,下游平原坦蕩,河曲發(fā)育。本文所建立的數(shù)學(xué)模型即來(lái)自于新疆瑪納斯河邊坡,為了便于研究,將其簡(jiǎn)化成理想邊坡,該邊坡坡角為40°,坡高為30m,其初始地下水位位于壩腳處。土體的彈性模量為10MPa,泊松比為0.3,黏聚力c′為15kPa,摩擦角φ′為30°,以此進(jìn)行降雨入滲模型的建立。

    基于上文降雨入滲計(jì)算力學(xué)原理及邊界條件的選擇,進(jìn)行有限元模型的建立。所建立的有限元模型如圖1所示,其element type選擇為CPE4P,滲透系數(shù)隨飽和度以及基質(zhì)吸力隨飽和度的變化可根據(jù)公式(1)和公式(2)得到。

    圖1 有限元網(wǎng)格模型

    利用distribution定義孔壓邊界條件,在左、右兩側(cè)水位以下的邊界上設(shè)置隨深度線性增加的靜水壓力孔,按照式(3)進(jìn)行設(shè)置:

    p=1010-Y

    (3)

    分析步的時(shí)間步長(zhǎng)為72h,以0.1作為初始的時(shí)間增量步,將孔壓設(shè)定為沿深度線性關(guān)系增加,降雨強(qiáng)度的時(shí)程賦值曲線Amp- 1見(jiàn)下表。

    時(shí)間和幅值的大小關(guān)系表

    對(duì)所研究的坡頂區(qū)域的降雨入滲邊界條件進(jìn)行設(shè)置,嵌入有效應(yīng)力有限元文件,即可進(jìn)行該模型的降雨入滲有限元分析。

    3 降雨入滲后的結(jié)果分析

    降雨發(fā)生后46h、72h的孔壓等值線云圖見(jiàn)圖2、圖3。由圖2、圖3可知:在降雨入滲發(fā)生后,孔壓的分布云圖與初始狀態(tài)發(fā)生顯著差異,在坡頂部范圍以下的吸力區(qū)發(fā)生范圍減小的現(xiàn)象,基質(zhì)吸力參數(shù)較初始狀態(tài)減?。浑S著降雨時(shí)間的增加,法向狀態(tài)延長(zhǎng),飽和度也隨之增大,孔隙水壓力參數(shù)也隨之增大。隨著降雨時(shí)間的增加,淺層土體的基質(zhì)吸力不斷減小甚至消失。

    而隨著降雨的減少直至停止,在時(shí)間不斷延長(zhǎng)的情況下,飽和度出現(xiàn)逐漸減小的情況,孔隙水壓力也隨之逐漸減小,土體淺層的基質(zhì)吸力出現(xiàn)逐漸增加的情況。

    圖2 降雨發(fā)生后46h的孔壓云圖

    圖3 降雨發(fā)生后72h的孔壓云圖

    圖4與圖5是降雨發(fā)生72h后的水平位移及沉降的等值線云圖。由圖4和圖5可知:發(fā)生在坡腳的最大水平位移是9.55mm;土坡的中部發(fā)生最大沉降現(xiàn)象,沉降位移是4.06mm。由于初始設(shè)置的分布條件,導(dǎo)致坡頂沒(méi)有出現(xiàn)最大沉降。

    圖4 降雨72h后的水平位移云圖

    圖5 降雨72h后的沉降云圖

    對(duì)坡頂進(jìn)行單獨(dú)分析,當(dāng)降雨入滲發(fā)生后,吸力出現(xiàn)降低情況,孔壓出現(xiàn)增加情況,有效應(yīng)力出現(xiàn)減小情況,發(fā)生了卸載回彈現(xiàn)象。除此之外,降雨入滲的不斷持續(xù),使得土體容重和含水率增加,進(jìn)而使得沉降和應(yīng)力均出現(xiàn)增加現(xiàn)象。

    位移矢量圖如圖6所示。從圖6可以看出,降雨入滲發(fā)生后,邊坡出現(xiàn)明顯的滑動(dòng)變形趨勢(shì),由此可知,邊坡穩(wěn)定性逐漸變低。

    圖6 降雨72h后的位移矢量圖

    圖7和圖8是降雨入滲23h及72h后的等效塑性應(yīng)變?cè)茍D。由圖7和圖8可知,坡腳淺層由于降雨入滲作用首先出現(xiàn)塑性區(qū),且塑性區(qū)沿坡面會(huì)不斷地向上擴(kuò)展。

    圖7 降雨入滲23h的等效塑性應(yīng)變?cè)茍D

    選取塑性區(qū)中心點(diǎn),給出等效塑性應(yīng)變時(shí)程曲線,如圖9所示。由圖9可見(jiàn),當(dāng)t=20h時(shí),塑性區(qū)出現(xiàn)在坡腳,隨后降雨入滲時(shí)間快速增加,當(dāng)時(shí)間增加到50h后,降雨情況開(kāi)始減小,此后,塑性區(qū)便不再向后擴(kuò)展,等效塑性應(yīng)變隨之也不再增加,保持穩(wěn)定不變。

    圖8 降雨入滲72h的等效塑性應(yīng)變?cè)茍D

    圖9 塑性中心單元等效塑性應(yīng)變時(shí)程曲線

    利用ABAQUS提供的combine函數(shù),繪出塑性區(qū)中心點(diǎn)和非塑性區(qū)中心點(diǎn)的有效應(yīng)力路徑,如圖10和11所示。

    圖10 塑性中心單元有效應(yīng)力路徑曲線

    圖11 非塑性區(qū)單元有效應(yīng)力路徑曲線

    由圖10和圖11可知,塑性中心單元和非塑性區(qū)中心單元的有效應(yīng)力路徑明顯出現(xiàn)了不同的力學(xué)特征。對(duì)于出現(xiàn)在坡腳的塑性中心單元(見(jiàn)圖10),由于降雨入滲,孔壓開(kāi)始增加,有效平均應(yīng)力出現(xiàn)減小的情況,而減小到一定閾值后,有效應(yīng)力路徑表現(xiàn)出屈服現(xiàn)象,其處于屈服面位置,此時(shí),應(yīng)力路徑向左下方逐漸下移,其沿著摩爾庫(kù)侖強(qiáng)度包線形成的屈服面朝左下方移動(dòng),當(dāng)降雨量減小后,移動(dòng)停止,此時(shí),吸力開(kāi)始增加,孔壓開(kāi)始減小,在有效應(yīng)力開(kāi)始增大后,其逐漸遠(yuǎn)離屈服面。

    而對(duì)于非塑性單元(見(jiàn)圖11),其處于土方內(nèi)部,上方的單元吸水后容重增加,導(dǎo)致塑性區(qū)單元的平均有效應(yīng)力和偏應(yīng)力都是增加的,直至降雨的后期,接近降雨結(jié)束時(shí)平均有效應(yīng)力和偏應(yīng)力都有所下降,即降雨入滲作用下主要是邊坡淺層可能會(huì)出現(xiàn)失穩(wěn)現(xiàn)象,該區(qū)域具有較高的失穩(wěn)傾向性。

    4 結(jié) 語(yǔ)

    基于ABAQUS軟件及其計(jì)算力學(xué)原理,以理想邊坡為研究對(duì)象,對(duì)土坡降雨入滲現(xiàn)象進(jìn)行了非飽和流固耦合分析,主要結(jié)論為:

    a. 在降雨入滲的作用下,位移矢量圖表明,邊坡有明顯的滑動(dòng)變形的趨勢(shì),邊坡穩(wěn)定性由此逐漸降低。

    b. 坡腳淺層會(huì)首先出現(xiàn)塑性區(qū),并沿著坡面不斷向上擴(kuò)展。

    c. 塑性點(diǎn)和非塑性點(diǎn)的有效應(yīng)力路徑表現(xiàn)出了不同的特征。

    [1] 李振宇. 利用地面核磁共振研究滑坡地下水特征及穩(wěn)定性[D]. 北京:中國(guó)地質(zhì)大學(xué), 2004.

    [2] 汪仕旭. 降雨入滲對(duì)土坡穩(wěn)定性的影響[D]. 廣州:廣東工業(yè)大學(xué), 2011.

    [3] 曾耀. 水庫(kù)庫(kù)岸老滑坡復(fù)活過(guò)程中地下水作用機(jī)制研究[D]. 北京:中國(guó)地質(zhì)大學(xué), 2010.

    [4] 季洪濤. 基于突變理論的土石壩壩坡安全穩(wěn)定分析[J]. 水利建設(shè)與管理, 2015(12): 24- 34.

    [5] 張家明. 植被發(fā)育斜坡非飽和帶土體大孔隙對(duì)降雨入滲影響研究[D]. 昆明:昆明理工大學(xué), 2013.

    [6] 周家文, 徐衛(wèi)亞, 鄧俊曄, 等. 降雨入滲條件下邊坡的穩(wěn)定性分析[J]. 水利學(xué)報(bào), 2008, 39(9):1066 - 1073.

    [7] 蔣中明, 熊小虎, 曾鈴. 基于FLAC~(3D)平臺(tái)的邊坡非飽和降雨入滲分析[J]. 巖土力學(xué), 2014, 35(30):855- 861.

    [8] 陳洪松, 邵明安, 王克林. 土壤初始含水率對(duì)坡面降雨入滲及土壤水分再分布的影響[J]. 農(nóng)業(yè)工程學(xué)報(bào), 2006, 25(1):44- 47.

    [9] 朱偉, 陳學(xué)東, 鐘小春. 降雨入滲規(guī)律的實(shí)測(cè)與分析[J]. 巖土力學(xué), 2006, 27(11):1873- 1879.

    [10] 王協(xié)群, 張有祥, 鄒維列, 等. 降雨入滲條件下非飽和路堤變形與邊坡的穩(wěn)定數(shù)值模擬[J]. 巖土力學(xué), 2010, 31(11):3640- 3644.

    StudyonstressstrainandstabilitybasedonABAQUSsoilsloperainfallinfiltration

    DU Min

    (No.3WaterConservancyManagementStation,XinjiangManasRiverBasinAuthorityWaterConservancyManagementCenter,Shihezi832000,China)

    Rainfall affects the mechanical parameters of soil, the influence of rainwater infiltration transient seepage flow on soil slope stress and strain is studied through ABAQUS. On the basis, an ideal slope is regarded as a case for modeling and calculation. The influence of rainfall and water accumulation is not considered. It is believed that hydrostatic pressure increases linearly with depth, and a CPE4P grid is established for calculation. In the post- processing, the permeability, pore pressure, horizontal longitudinal displacement, deviatoric stress, plastic strain and stress path, etc. are analyzed, and a series of related conclusions of stress, strain and stability under rainfall infiltration action are obtained.

    rainfall; seepage; ABAQUS; instability; slope; plastic zone

    10.16616/j.cnki.10- 1326/TV.2017.010.012

    TV211.2

    A

    2096-0131(2017)010-0039-05

    国产精品久久久久久久电影 | 国产人伦9x9x在线观看| 欧美日韩一级在线毛片| www国产在线视频色| 日本在线视频免费播放| 麻豆一二三区av精品| 无限看片的www在线观看| 九九在线视频观看精品| ponron亚洲| 丁香欧美五月| 国产精品九九99| 香蕉丝袜av| 午夜福利在线观看吧| 国产真实乱freesex| 欧美色欧美亚洲另类二区| 欧美另类亚洲清纯唯美| 久久久久久久久中文| 激情在线观看视频在线高清| www.自偷自拍.com| 色吧在线观看| 91在线精品国自产拍蜜月 | 啪啪无遮挡十八禁网站| 一级作爱视频免费观看| 色尼玛亚洲综合影院| 国产欧美日韩精品一区二区| 国产黄片美女视频| 亚洲 国产 在线| 国产不卡一卡二| 一二三四社区在线视频社区8| 综合色av麻豆| 国产欧美日韩一区二区精品| 99久久成人亚洲精品观看| 国产成人av激情在线播放| 亚洲国产日韩欧美精品在线观看 | 好男人电影高清在线观看| 成人三级黄色视频| 国产高清视频在线观看网站| 91在线精品国自产拍蜜月 | 国产亚洲精品久久久com| 国产美女午夜福利| 亚洲天堂国产精品一区在线| 亚洲中文日韩欧美视频| 日韩欧美在线乱码| 日韩中文字幕欧美一区二区| 亚洲精品一区av在线观看| 他把我摸到了高潮在线观看| a级毛片a级免费在线| 免费搜索国产男女视频| 午夜激情欧美在线| 高潮久久久久久久久久久不卡| 在线十欧美十亚洲十日本专区| 不卡一级毛片| 久久精品aⅴ一区二区三区四区| 99久久无色码亚洲精品果冻| 国产精品亚洲av一区麻豆| 九色国产91popny在线| 中文资源天堂在线| 国产一级毛片七仙女欲春2| 少妇人妻一区二区三区视频| 午夜福利免费观看在线| 好男人在线观看高清免费视频| 久久久久国产一级毛片高清牌| 夜夜躁狠狠躁天天躁| 亚洲va日本ⅴa欧美va伊人久久| 悠悠久久av| 男女下面进入的视频免费午夜| 97人妻精品一区二区三区麻豆| 久久人人精品亚洲av| 久久久成人免费电影| 国产一区二区三区在线臀色熟女| 一边摸一边抽搐一进一小说| 91麻豆av在线| 真人做人爱边吃奶动态| 亚洲专区国产一区二区| 在线观看一区二区三区| 88av欧美| 亚洲av第一区精品v没综合| 精品熟女少妇八av免费久了| 法律面前人人平等表现在哪些方面| 国产成人精品久久二区二区免费| 日韩免费av在线播放| 久久这里只有精品中国| 嫩草影院精品99| 国产亚洲av嫩草精品影院| 欧美黄色淫秽网站| 亚洲av成人一区二区三| 两个人的视频大全免费| 日韩国内少妇激情av| 精品一区二区三区视频在线观看免费| 18禁国产床啪视频网站| 午夜福利视频1000在线观看| 久久久久久久久免费视频了| 叶爱在线成人免费视频播放| 99热只有精品国产| tocl精华| 一级毛片高清免费大全| 久久国产乱子伦精品免费另类| 亚洲av第一区精品v没综合| 美女扒开内裤让男人捅视频| 一个人看的www免费观看视频| 黄色女人牲交| 亚洲欧美日韩高清专用| 欧美绝顶高潮抽搐喷水| 亚洲av美国av| 在线看三级毛片| 一本综合久久免费| 他把我摸到了高潮在线观看| 国产成人aa在线观看| 9191精品国产免费久久| 热99在线观看视频| 久9热在线精品视频| 免费一级毛片在线播放高清视频| 亚洲国产看品久久| 免费看十八禁软件| 日韩大尺度精品在线看网址| 男女午夜视频在线观看| 91在线精品国自产拍蜜月 | 国产免费男女视频| 1024香蕉在线观看| 精品一区二区三区四区五区乱码| 国产精品亚洲美女久久久| 97人妻精品一区二区三区麻豆| 久久亚洲真实| 伦理电影免费视频| 黄色成人免费大全| 巨乳人妻的诱惑在线观看| 久久久久久大精品| 精品一区二区三区视频在线 | 国产av麻豆久久久久久久| 小说图片视频综合网站| 免费av毛片视频| 亚洲欧美一区二区三区黑人| 亚洲成a人片在线一区二区| 欧美最黄视频在线播放免费| 日韩欧美在线乱码| 久久久久久人人人人人| 长腿黑丝高跟| 久久性视频一级片| 中文字幕久久专区| 人妻夜夜爽99麻豆av| 色综合亚洲欧美另类图片| 国产伦人伦偷精品视频| 女人被狂操c到高潮| 日日摸夜夜添夜夜添小说| 亚洲熟女毛片儿| 男人舔女人下体高潮全视频| 岛国在线免费视频观看| 色视频www国产| 国产激情欧美一区二区| 99久久久亚洲精品蜜臀av| 日韩欧美一区二区三区在线观看| 亚洲精华国产精华精| 国产单亲对白刺激| 国产成+人综合+亚洲专区| 国内精品美女久久久久久| 午夜福利在线在线| 亚洲精品色激情综合| 俄罗斯特黄特色一大片| 亚洲成a人片在线一区二区| 99热精品在线国产| 欧美+亚洲+日韩+国产| 变态另类丝袜制服| 操出白浆在线播放| 欧美高清成人免费视频www| 两个人的视频大全免费| 国产亚洲欧美在线一区二区| 18禁黄网站禁片午夜丰满| 黑人欧美特级aaaaaa片| 亚洲成人久久性| 久久天堂一区二区三区四区| 亚洲精品久久国产高清桃花| 国产精品永久免费网站| 国语自产精品视频在线第100页| 最好的美女福利视频网| 久久久精品欧美日韩精品| 亚洲色图 男人天堂 中文字幕| 最近最新中文字幕大全电影3| 成年版毛片免费区| 国产免费av片在线观看野外av| 国产高清有码在线观看视频| 少妇的丰满在线观看| 欧美一区二区精品小视频在线| 午夜亚洲福利在线播放| 伦理电影免费视频| 国产精品久久久久久亚洲av鲁大| 99国产精品一区二区蜜桃av| 色av中文字幕| 亚洲欧美日韩东京热| 男女午夜视频在线观看| 欧美性猛交╳xxx乱大交人| 狂野欧美白嫩少妇大欣赏| 亚洲 欧美一区二区三区| 黄色 视频免费看| 久久久久久久久久黄片| 制服丝袜大香蕉在线| 黄色片一级片一级黄色片| 国产又黄又爽又无遮挡在线| 999久久久精品免费观看国产| 岛国视频午夜一区免费看| 免费观看的影片在线观看| 非洲黑人性xxxx精品又粗又长| 国产1区2区3区精品| 日韩欧美 国产精品| 国产高清视频在线观看网站| 欧美激情在线99| 性欧美人与动物交配| 99久久99久久久精品蜜桃| 夜夜看夜夜爽夜夜摸| 亚洲成人久久性| 日韩欧美三级三区| 精品熟女少妇八av免费久了| 欧美国产日韩亚洲一区| 黄色日韩在线| 中文字幕人妻丝袜一区二区| 18禁黄网站禁片午夜丰满| 一区二区三区国产精品乱码| 一级a爱片免费观看的视频| 99国产综合亚洲精品| 99久久成人亚洲精品观看| 国产精品久久久久久亚洲av鲁大| 欧美乱妇无乱码| 国产伦在线观看视频一区| 免费电影在线观看免费观看| 99久国产av精品| 1000部很黄的大片| 18禁裸乳无遮挡免费网站照片| 亚洲中文日韩欧美视频| 亚洲18禁久久av| 久久久国产成人精品二区| 亚洲中文字幕日韩| 男女之事视频高清在线观看| 免费一级毛片在线播放高清视频| 丰满人妻一区二区三区视频av | 两人在一起打扑克的视频| 少妇的逼水好多| 大型黄色视频在线免费观看| 搡老妇女老女人老熟妇| 欧美大码av| 一进一出抽搐gif免费好疼| 麻豆国产97在线/欧美| 欧美日韩福利视频一区二区| 成人永久免费在线观看视频| 草草在线视频免费看| 麻豆av在线久日| 九色国产91popny在线| 欧美日韩瑟瑟在线播放| 老熟妇仑乱视频hdxx| 欧美高清成人免费视频www| 精品福利观看| 一区二区三区高清视频在线| 免费在线观看影片大全网站| 日本 av在线| 国产又色又爽无遮挡免费看| avwww免费| 久久中文字幕人妻熟女| 熟女人妻精品中文字幕| 国产 一区 欧美 日韩| 亚洲成av人片在线播放无| 精品一区二区三区视频在线观看免费| 久久亚洲精品不卡| 国产黄色小视频在线观看| 欧美黑人巨大hd| 成年女人毛片免费观看观看9| 99精品在免费线老司机午夜| av天堂中文字幕网| 欧美av亚洲av综合av国产av| 欧美黑人欧美精品刺激| 一夜夜www| 国产成人影院久久av| 免费在线观看成人毛片| 欧美成狂野欧美在线观看| 两人在一起打扑克的视频| 国产av一区在线观看免费| 成人欧美大片| 岛国在线观看网站| 悠悠久久av| 男人舔女人下体高潮全视频| 欧美午夜高清在线| 成在线人永久免费视频| 久久久久九九精品影院| 欧美黄色淫秽网站| 一级毛片女人18水好多| 搡老熟女国产l中国老女人| 99久久久亚洲精品蜜臀av| 国产精品一及| 91久久精品国产一区二区成人 | 亚洲欧美一区二区三区黑人| 91麻豆av在线| 日韩精品青青久久久久久| 精品乱码久久久久久99久播| 在线免费观看不下载黄p国产 | 亚洲成人中文字幕在线播放| 久久久国产成人免费| 特大巨黑吊av在线直播| 亚洲中文字幕一区二区三区有码在线看 | 19禁男女啪啪无遮挡网站| 日本撒尿小便嘘嘘汇集6| av福利片在线观看| 欧美日韩中文字幕国产精品一区二区三区| 99久久综合精品五月天人人| 精品国产亚洲在线| 在线观看美女被高潮喷水网站 | 美女扒开内裤让男人捅视频| 夜夜看夜夜爽夜夜摸| 母亲3免费完整高清在线观看| 特级一级黄色大片| 欧美性猛交╳xxx乱大交人| 日本三级黄在线观看| 国内精品美女久久久久久| av在线天堂中文字幕| 亚洲第一电影网av| 1000部很黄的大片| 欧美在线一区亚洲| 国产激情久久老熟女| 精品国产超薄肉色丝袜足j| 国产亚洲欧美在线一区二区| 成人三级做爰电影| 亚洲成人免费电影在线观看| 国产精品av久久久久免费| 久久久久精品国产欧美久久久| 一区二区三区高清视频在线| 亚洲国产欧美一区二区综合| 亚洲av中文字字幕乱码综合| 一a级毛片在线观看| 亚洲欧美日韩高清专用| 欧美不卡视频在线免费观看| 亚洲国产欧洲综合997久久,| 人妻久久中文字幕网| 亚洲电影在线观看av| 国产精品九九99| 老司机午夜福利在线观看视频| 成人特级黄色片久久久久久久| 国产一区二区激情短视频| 亚洲精品色激情综合| 国产成人欧美在线观看| 婷婷亚洲欧美| 亚洲av美国av| 一卡2卡三卡四卡精品乱码亚洲| 日本黄大片高清| 香蕉丝袜av| www.熟女人妻精品国产| 97碰自拍视频| 午夜福利在线观看吧| 日韩欧美在线乱码| 欧美黑人欧美精品刺激| 国产精品 欧美亚洲| 嫩草影视91久久| 色精品久久人妻99蜜桃| 黑人操中国人逼视频| 1024香蕉在线观看| 国产69精品久久久久777片 | 一夜夜www| xxx96com| 国产真人三级小视频在线观看| 免费观看的影片在线观看| 少妇熟女aⅴ在线视频| 国产av不卡久久| 一本一本综合久久| 日韩国内少妇激情av| 亚洲精品色激情综合| 亚洲五月天丁香| 精品日产1卡2卡| 欧美一区二区精品小视频在线| 啪啪无遮挡十八禁网站| 搡老熟女国产l中国老女人| 97超视频在线观看视频| 99久久综合精品五月天人人| 97碰自拍视频| 精品一区二区三区视频在线 | 久久久久国产一级毛片高清牌| 麻豆久久精品国产亚洲av| 亚洲av电影在线进入| 最新美女视频免费是黄的| 国产三级在线视频| 亚洲欧美日韩高清专用| 变态另类成人亚洲欧美熟女| 国产精品女同一区二区软件 | 不卡一级毛片| 一夜夜www| 一级毛片精品| 色综合亚洲欧美另类图片| 精品无人区乱码1区二区| 国产免费男女视频| 精品久久久久久,| 国产精品久久视频播放| 国产伦一二天堂av在线观看| 免费高清视频大片| 亚洲天堂国产精品一区在线| 日韩 欧美 亚洲 中文字幕| 国产麻豆成人av免费视频| av天堂中文字幕网| 性色av乱码一区二区三区2| 可以在线观看毛片的网站| 亚洲欧美日韩无卡精品| 在线观看免费午夜福利视频| 九九在线视频观看精品| 亚洲aⅴ乱码一区二区在线播放| 久久精品国产综合久久久| 精品人妻1区二区| 给我免费播放毛片高清在线观看| 两性午夜刺激爽爽歪歪视频在线观看| 亚洲欧美一区二区三区黑人| 欧美日韩精品网址| 99久久无色码亚洲精品果冻| 国产精品久久久久久久电影 | 国产久久久一区二区三区| 久久精品国产亚洲av香蕉五月| 国产人伦9x9x在线观看| 国产av不卡久久| 成人精品一区二区免费| 偷拍熟女少妇极品色| 欧美乱色亚洲激情| 日本a在线网址| 国产毛片a区久久久久| 视频区欧美日本亚洲| 在线观看舔阴道视频| 18禁美女被吸乳视频| 色老头精品视频在线观看| 亚洲国产精品成人综合色| 国产精品99久久99久久久不卡| a在线观看视频网站| 身体一侧抽搐| 亚洲欧美激情综合另类| 国产又色又爽无遮挡免费看| 国产又黄又爽又无遮挡在线| 国产三级黄色录像| 国产亚洲精品一区二区www| 两个人视频免费观看高清| 国产精品98久久久久久宅男小说| 嫁个100分男人电影在线观看| 国产极品精品免费视频能看的| 99精品久久久久人妻精品| 成人国产一区最新在线观看| 成人鲁丝片一二三区免费| 最新中文字幕久久久久 | 国产一区在线观看成人免费| 久久午夜综合久久蜜桃| 久久精品综合一区二区三区| 亚洲va日本ⅴa欧美va伊人久久| 嫩草影院入口| 欧美日韩黄片免| 精品熟女少妇八av免费久了| 日韩成人在线观看一区二区三区| www.熟女人妻精品国产| 国产午夜精品久久久久久| 亚洲国产欧美人成| 一级毛片高清免费大全| 91在线精品国自产拍蜜月 | 免费电影在线观看免费观看| 免费一级毛片在线播放高清视频| 国产伦人伦偷精品视频| 九九在线视频观看精品| 99久久国产精品久久久| 99re在线观看精品视频| 中文在线观看免费www的网站| www.999成人在线观看| 真人做人爱边吃奶动态| 少妇裸体淫交视频免费看高清| 手机成人av网站| 免费观看的影片在线观看| 好男人电影高清在线观看| 少妇的丰满在线观看| 丰满的人妻完整版| 天堂动漫精品| 国产伦一二天堂av在线观看| 狠狠狠狠99中文字幕| 一二三四在线观看免费中文在| 青草久久国产| 最好的美女福利视频网| 高清在线国产一区| 日韩三级视频一区二区三区| 国产亚洲精品综合一区在线观看| 亚洲av熟女| 动漫黄色视频在线观看| 桃色一区二区三区在线观看| 亚洲美女视频黄频| 中文字幕久久专区| 久久久精品欧美日韩精品| 国产精品99久久久久久久久| 国产成人精品无人区| 午夜福利在线观看免费完整高清在 | 一级a爱片免费观看的视频| 亚洲精品乱码久久久v下载方式 | 最近最新免费中文字幕在线| 久久久久久九九精品二区国产| x7x7x7水蜜桃| 久久精品夜夜夜夜夜久久蜜豆| 精品久久久久久成人av| 国产蜜桃级精品一区二区三区| 亚洲 欧美一区二区三区| 一个人免费在线观看电影 | bbb黄色大片| 少妇人妻一区二区三区视频| xxx96com| 久久中文字幕一级| 一级黄色大片毛片| 亚洲人与动物交配视频| 99re在线观看精品视频| 啦啦啦观看免费观看视频高清| 最近视频中文字幕2019在线8| 国产精品一区二区三区四区免费观看 | 一级作爱视频免费观看| 五月伊人婷婷丁香| 国产又黄又爽又无遮挡在线| a级毛片a级免费在线| 亚洲国产精品合色在线| 久久久久国内视频| 又粗又爽又猛毛片免费看| 欧美另类亚洲清纯唯美| 国产黄a三级三级三级人| 无限看片的www在线观看| 欧美日韩乱码在线| 成年版毛片免费区| 18禁黄网站禁片午夜丰满| 免费看a级黄色片| svipshipincom国产片| 琪琪午夜伦伦电影理论片6080| 色老头精品视频在线观看| 国产欧美日韩一区二区精品| 嫩草影视91久久| 校园春色视频在线观看| 19禁男女啪啪无遮挡网站| 好看av亚洲va欧美ⅴa在| 国产精品久久视频播放| 久久精品人妻少妇| 欧美丝袜亚洲另类 | 最新美女视频免费是黄的| 久久久久久大精品| 久久国产精品影院| 中文字幕精品亚洲无线码一区| 久久这里只有精品中国| 99久久精品一区二区三区| 国产精品久久久久久人妻精品电影| 国产精品1区2区在线观看.| 成人av在线播放网站| 男人舔女人下体高潮全视频| 国产1区2区3区精品| 久久亚洲真实| 国产精品久久久久久亚洲av鲁大| 级片在线观看| 香蕉丝袜av| 国产99白浆流出| 首页视频小说图片口味搜索| 波多野结衣高清作品| 欧美精品啪啪一区二区三区| 别揉我奶头~嗯~啊~动态视频| 久久人人精品亚洲av| 亚洲avbb在线观看| 久久精品aⅴ一区二区三区四区| 欧美成狂野欧美在线观看| 国产精品电影一区二区三区| 亚洲狠狠婷婷综合久久图片| 搞女人的毛片| 这个男人来自地球电影免费观看| 国内精品美女久久久久久| 国内精品久久久久久久电影| 精品久久久久久,| 亚洲男人的天堂狠狠| 九九久久精品国产亚洲av麻豆 | 欧美不卡视频在线免费观看| 一边摸一边抽搐一进一小说| 国产极品精品免费视频能看的| 在线观看免费午夜福利视频| 人人妻,人人澡人人爽秒播| 亚洲av电影不卡..在线观看| 国产91精品成人一区二区三区| 国产欧美日韩精品一区二区| 制服丝袜大香蕉在线| 麻豆成人av在线观看| 国产美女午夜福利| 美女免费视频网站| 国产v大片淫在线免费观看| 国产精品一及| 国产激情欧美一区二区| 国产精品女同一区二区软件 | 午夜福利免费观看在线| 日本熟妇午夜| 一边摸一边抽搐一进一小说| av在线蜜桃| 亚洲一区高清亚洲精品| 国产免费av片在线观看野外av| 欧美最黄视频在线播放免费| 免费观看人在逋| 亚洲欧美精品综合一区二区三区| 欧美乱妇无乱码| 国产野战对白在线观看| 国产av一区在线观看免费| 国产成年人精品一区二区| 日日夜夜操网爽| 99国产综合亚洲精品| 岛国在线观看网站| 男人和女人高潮做爰伦理| 国产av在哪里看| 日韩欧美在线二视频| avwww免费| 日日夜夜操网爽| 51午夜福利影视在线观看| 色老头精品视频在线观看| 日韩高清综合在线| 日韩免费av在线播放| 又爽又黄无遮挡网站| 香蕉av资源在线| 国产成人精品久久二区二区免费| 99国产精品一区二区蜜桃av| 国产一区二区三区视频了| 在线观看日韩欧美| 日韩欧美精品v在线| 99久久成人亚洲精品观看| 亚洲中文av在线| 成人高潮视频无遮挡免费网站| 日本一二三区视频观看| 女人高潮潮喷娇喘18禁视频| 国内久久婷婷六月综合欲色啪| 国产三级黄色录像| 757午夜福利合集在线观看|