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

    基于彈塑性區(qū)間有限元的邊坡穩(wěn)定性分析

    2016-05-17 03:33:24喻和平袁明明陳熙源張家德

    喻和平,袁明明,張 聰,陳熙源,張家德

    (長(zhǎng)沙理工大學(xué) 水利工程學(xué)院,湖南 長(zhǎng)沙 410114)

    ?

    基于彈塑性區(qū)間有限元的邊坡穩(wěn)定性分析

    喻和平,袁明明,張聰,陳熙源,張家德

    (長(zhǎng)沙理工大學(xué) 水利工程學(xué)院,湖南 長(zhǎng)沙 410114)

    摘要:區(qū)間分析法是工程中對(duì)不確定問題進(jìn)行分析的方法之一?,F(xiàn)有區(qū)間分析法大多采用區(qū)間參數(shù)攝動(dòng)法求解,其存在明顯的局限性,特別對(duì)于彈塑性問題需采用增量法計(jì)算時(shí),很難得到有效的解答。在邊坡穩(wěn)定性計(jì)算中采用區(qū)間變量來描述工程中的不確定性,將區(qū)間數(shù)學(xué)與彈塑性有限元相結(jié)合,建立了彈塑性區(qū)間有限元,推導(dǎo)了基于彈塑性區(qū)間有限元方法的邊坡穩(wěn)定性區(qū)間安全系數(shù)計(jì)算公式,這一方法能有效的避免區(qū)間參數(shù)攝動(dòng)法的局限性。通過工程實(shí)例,驗(yàn)證了該方法的可行性和有效性。在此基礎(chǔ)上探討了巖土體參數(shù)對(duì)區(qū)間安全系數(shù)的影響。結(jié)果表明,內(nèi)摩擦角對(duì)邊坡安全系數(shù)的影響明顯,重度、黏聚力影響較大,彈性模量和泊松比影響小。

    關(guān)鍵詞:區(qū)間分析;彈塑性區(qū)間有限元;邊坡穩(wěn)定;區(qū)間安全系數(shù)

    區(qū)間分析法(區(qū)間數(shù)學(xué))是20世紀(jì)60年代后期發(fā)展起來的較活躍的計(jì)算數(shù)學(xué)分支,是一種不確定性分析方法。Moore R E在1966年奠定了區(qū)間數(shù)學(xué)的理論基礎(chǔ)。我國(guó)王德人和陳塑寰[1-2]等人發(fā)展了其理論。近幾十年來,區(qū)間有限元法已在結(jié)構(gòu)分析和巖土工程領(lǐng)域得到了廣泛的研究與應(yīng)用。如:Rao S S和Mullen等[3-6]結(jié)合工程中的諸多不確定性,把區(qū)間方法引入到工程中,建立了區(qū)間有限元;邱志平等[7]針對(duì)區(qū)間數(shù)學(xué)擴(kuò)張問題,把攝動(dòng)方法引入到區(qū)間有限元結(jié)構(gòu)分析中,提出了基于攝動(dòng)法的區(qū)間有限元;邵國(guó)建等[8]基于單元的子區(qū)間攝動(dòng)法求解線彈性區(qū)間有限元,得到了結(jié)構(gòu)抗滑區(qū)間安全系數(shù);黃耀英等[9]采用區(qū)間參數(shù)攝動(dòng)法,推導(dǎo)了彈塑性區(qū)間有限元計(jì)算公式,并應(yīng)用于結(jié)構(gòu)分析中;喻和平等[10-11]總結(jié)了區(qū)間有限元控制方程的求解方法,并推導(dǎo)了基于單元的彈性區(qū)間有限元計(jì)算公式和基于應(yīng)力水平的邊坡區(qū)間安全系數(shù)。現(xiàn)階段,相關(guān)學(xué)者雖取得了一定成就,但以上區(qū)間參數(shù)攝動(dòng)法的求解存在較明顯的局限性。自由度較少時(shí),計(jì)算速度快,能得到有效的解答。而自由度較多,特別是彈塑性問題需采用增量法進(jìn)行計(jì)算時(shí),不僅計(jì)算效率和精度需進(jìn)一步研究,且很難得到有效的解答,如邊坡等問題。

    為解決以上的缺陷,在邊坡計(jì)算中采用區(qū)間變量來描述工程中的不確定性,將區(qū)間數(shù)學(xué)與彈塑性有限元相結(jié)合,建立了彈塑性區(qū)間有限元,推導(dǎo)出基于彈塑性區(qū)間有限元方法的邊坡穩(wěn)定性區(qū)間安全系數(shù)計(jì)算公式。該方法是在假定滑動(dòng)面的基礎(chǔ)上,對(duì)參數(shù)在區(qū)間范圍內(nèi)根據(jù)一定步長(zhǎng)進(jìn)行搜索賦值,然后在確定性有限元的基礎(chǔ)上形成彈塑性區(qū)間有限元求解,相比區(qū)間攝動(dòng)法、彈性區(qū)間有限元法,該方法計(jì)算的安全系數(shù)區(qū)間值更為精確。本文以某邊坡為例,驗(yàn)證其效果,以期為工程實(shí)際提供參考。

    1彈塑性問題的區(qū)間有限元

    采用增量法對(duì)邊坡進(jìn)行穩(wěn)定性分析,區(qū)間有限元控制方程如下式(1):

    KIΔUI=ΔRI

    (1)

    假設(shè)材料滿足D-P屈服準(zhǔn)則,則該屈服函數(shù)如式(2):

    (2)

    式中:α、k為材料常數(shù),與材料的內(nèi)摩擦角和黏聚力有關(guān);I1、J2為應(yīng)力第一不變量和應(yīng)力偏量的第二不變量。由于材料參數(shù)的獲取存在不確定性因素,這種不確定性通過區(qū)間數(shù)來表示。因此,各材料參數(shù)都被看作區(qū)間變量,則屈服函數(shù)也是區(qū)間函數(shù),其函數(shù)如下式(3):

    (3)

    2邊坡區(qū)間安全系數(shù)

    圖1單元體的平面應(yīng)力狀態(tài)

    (4)

    (5)

    根據(jù)邊坡安全系數(shù)定義,邊坡安全系數(shù)為潛在滑動(dòng)面上各段抗剪力之和與其剪切力之和的比值。則邊坡區(qū)間安全系數(shù)計(jì)算公見式(6)~式(8)。

    (6)

    (7)

    (8)

    3程序編制

    針對(duì)區(qū)間數(shù)學(xué)擴(kuò)張問題,本文利用確定性的有限元程序[12-14],通過參數(shù)區(qū)間組合法來確定邊坡的區(qū)間安全系數(shù)。主要步驟為:(1) 依次對(duì)巖土體內(nèi)摩擦角,黏聚力,重度按一定步長(zhǎng)進(jìn)行賦值,在各區(qū)間參數(shù)范圍內(nèi)得到盡量多的參數(shù)值,然后對(duì)賦值的內(nèi)摩擦角,黏聚力,重度進(jìn)行隨機(jī)組合形成盡量多的參數(shù)組合;(2) 對(duì)形成的參數(shù)組合進(jìn)行確定性的彈塑性有限元分析,得到安全系數(shù)的確定值;(3) 對(duì)各組計(jì)算結(jié)果進(jìn)行比較,尋找安全系數(shù)的最大最小值,作為邊坡區(qū)間安全系數(shù)的上下值。

    4實(shí)例分析

    4.1工程概況與模型建立

    某邊坡為一順層基巖滑坡,經(jīng)局部解體而成。坡區(qū)分布的第四系主要為崩、坡積層及堆積層,物質(zhì)組成復(fù)雜,主要分三層:礫質(zhì)黏土夾少量碎石;黏土夾灰?guī)r大塊石、塊石、碎石等;砂巖碎石、碎屑土層,物質(zhì)成分為灰黃、淺黃色、灰綠色、雜色黏土夾碎石碎屑[15]。邊坡計(jì)算剖面圖和有限元計(jì)算模型見圖2、圖3,從上到下依次為覆蓋層、滑帶層和基巖,各層巖土體材料計(jì)算參數(shù)均值見表1,計(jì)算時(shí)左右兩邊約束水平位移,下端約束水平和垂直位移,上端自由。

    圖2 邊坡計(jì)算剖面圖

    圖3 邊坡有限元計(jì)算模型

    4.2敏感性分析

    為研究該邊坡巖土體參數(shù)對(duì)邊坡穩(wěn)定性的影響程度,在假設(shè)其它參數(shù)為均值的基礎(chǔ)上,分別通過滑動(dòng)帶巖土體彈性模量、泊松比、重度、黏聚力及內(nèi)摩擦角的區(qū)間參數(shù)變異系數(shù)(CV=(標(biāo)準(zhǔn)差÷基準(zhǔn)值)×100%)的變化對(duì)安全系數(shù)的影響程度,來反映巖土體參數(shù)對(duì)邊坡穩(wěn)定性的影響。初始計(jì)算參數(shù)值為各巖體參數(shù)的平均值,結(jié)果見圖4~圖8。

    (υ=[0.31,0.49])

    (φ=[12°,21°])

    (E=[10,19] MPa)

    (c=[15,24] kPa)

    圖7黏聚力與安全系數(shù)的變異系數(shù)關(guān)系曲線圖

    由圖4~圖8可知,泊松比、內(nèi)摩擦角、彈性模量、黏聚力與重度的變異系數(shù)對(duì)安全系數(shù)的敏感度,依次為0.00780,0.89141,0.00455,0.05002,0.05968。表明內(nèi)摩擦角對(duì)邊坡安全系數(shù)的影響明顯,重度、黏聚力影響較大,彈性模量和泊松比影響小。因此為減小計(jì)算量,只對(duì)重度、黏聚力和內(nèi)摩擦角采用區(qū)間計(jì)算,對(duì)彈性模量和泊松比采取確定性計(jì)算。

    (γ=[16.5,25.5] kN/m3)

    圖8重度與安全系數(shù)的變異系數(shù)關(guān)系曲線圖

    4.3安全系數(shù)區(qū)間計(jì)算

    邊坡采用平面四邊形單元,通過編制的彈塑性區(qū)間有限元程序進(jìn)行計(jì)算,計(jì)算區(qū)間參數(shù)見表2,計(jì)算結(jié)果見表3、表4。由表3、表4可知,區(qū)間安全系數(shù)為[1.26,1.46],表明該邊坡是基本穩(wěn)定的。區(qū)間安全系數(shù)端點(diǎn)值分別對(duì)應(yīng)著各區(qū)間參數(shù)的端點(diǎn)值,巖土體的各參數(shù)對(duì)安全系數(shù)的影響是單調(diào)的,說明端點(diǎn)組合法在邊坡區(qū)間穩(wěn)定性計(jì)算中存在一定的合理性。區(qū)間有限元法所得的區(qū)間安全系數(shù)的下限值對(duì)應(yīng)著滑帶及以上巖土體材料內(nèi)摩擦角和黏聚力的下限值、重度的上限值,滑帶以下材料內(nèi)摩擦角和黏聚力的上限值、重度的下限值。區(qū)間安全系數(shù)的上限值對(duì)應(yīng)著滑帶及以上巖土體材料內(nèi)摩擦角和黏聚力的上限值、重度的下限值,滑帶以下材料內(nèi)摩擦角和黏聚力的下限值、重度的上限值。以滑帶為分界線,滑帶以上材料的內(nèi)摩擦角和黏聚力對(duì)安全系數(shù)的影響為單調(diào)增加,重度對(duì)安全系數(shù)的影響單調(diào)減小,滑帶以下材料的內(nèi)摩擦角和黏聚力對(duì)安全系數(shù)的影響為單調(diào)減小,重度對(duì)安全系數(shù)的影響單調(diào)增加。

    表2 彈塑性區(qū)間有限元計(jì)算參數(shù)

    表3 區(qū)間安全系數(shù)下端值

    表4 區(qū)間安全系數(shù)上端值

    5結(jié)論

    通過分析了區(qū)間參數(shù)攝動(dòng)法求解過程中存在的局限性問題,基于確定性有限元的求解方法,編寫了彈塑性區(qū)間有限元計(jì)算程序,采用有限元計(jì)算結(jié)果求解了邊坡區(qū)間安全系數(shù),并對(duì)不確定性的邊坡穩(wěn)定性問題進(jìn)行了評(píng)價(jià)。其主要結(jié)論如下:

    (1) 基于區(qū)間數(shù)學(xué)分析法,結(jié)合有限元計(jì)算,形成區(qū)間有限元法,不失為一種可行且有效的邊坡穩(wěn)定性計(jì)算方法。工程實(shí)例驗(yàn)證表明,該方法是合理、可行、有效的。

    (2) 針對(duì)區(qū)間分析法分析過程中核心—區(qū)間擴(kuò)張問題,借助于確定性的有限元求解方法求解區(qū)間有限元控制方程,使得求解結(jié)果更精確。

    (3) 參數(shù)敏感性分析表明,內(nèi)摩擦角、重度、黏聚力、彈性模量及泊松比對(duì)邊坡穩(wěn)定性均有影響。其中,內(nèi)摩擦角對(duì)邊坡安全系數(shù)的影響明顯,重度、黏聚力影響較大,彈性模量和泊松比影響小。

    參考文獻(xiàn):

    [1]王德生,張連生,鄧乃揚(yáng),非線性方程的區(qū)間算法[M].上海:上??茖W(xué)技術(shù)出版,1987:200-258.

    [2]陳塑寰.結(jié)構(gòu)動(dòng)態(tài)設(shè)計(jì)的矩陣攝動(dòng)理論[M].北京:科學(xué)出版社,1999:1-60.

    [3]Rao S S, Sawyer P J. Fuzzy finite element approach for analysis of imprecisely defined systems[J]. Aiaa Journal, 1995,33(12):2364-2370.

    [4]Rao S S, Berke L. Analysis of uncertain structural systems using interval analysis[J].Aiaa Journal,1997,35(4):727-735.

    [5]Mullen R L, Muhanna R L. Bounds of structural response for all possible loading Combinations[J]. Journal of Structural Engineering, 1999,125(1):98-106.

    [6]Muhanna R L, Mullen R L. Uncertainty mechanics problems-interval-based approach[J]. Journal of Engineering Mechanics, ASCE, 2001,127(6):557-566.

    [7]Qiu Z, Elishakoff I. Anti-optimization structures with Large Uncertain but nonrandom parameters via interval analysis[J]. Computer Methods in Applied Mechanics and Engineering, 1998,152(S3/S4):361-372.

    [8]邵國(guó)建,蘇靜波.區(qū)間有限元方法及其在抗滑穩(wěn)定性分析中的應(yīng)用[J].應(yīng)用數(shù)學(xué)和力學(xué),2007,28(4):471-478.

    [9]黃耀英,吳中如,蘇靜波.基于區(qū)間參數(shù)攝動(dòng)法的三維塑性區(qū)間有限元探討[J].計(jì)算力學(xué)學(xué)報(bào),2008,25(S1):6-10.

    [10]喻和平.區(qū)間分析理論及其在邊坡工程中的應(yīng)用[D].南京:河海大學(xué),2006:25-36.

    [11]喻和平,徐衛(wèi)亞.區(qū)間有限元在邊坡穩(wěn)定分析中的應(yīng)用[J].水利水電技術(shù),2005,36(6):33-35.

    [12]于生飛,陳征宙,張明瑞,等.基于區(qū)間不確定分析方法的邊坡穩(wěn)定性分析[J].工程地質(zhì)學(xué)報(bào),2012,20(2):228-233.

    [13]楊輝.邊坡穩(wěn)定分析的有限元極限平衡法原理及程序?qū)崿F(xiàn)[D].北京:北京交通大學(xué),2014:30-45.

    [14]陳曦,劉春杰.有限元強(qiáng)度折減法中安全系數(shù)的搜索算法[J].巖土工程學(xué)報(bào),2010,32(9):1443-1447.

    [15]周建軍,徐瑞春,柳景華.大巖淌邊坡降雨-應(yīng)力-非飽和滲流研究[J].長(zhǎng)江科學(xué)院院報(bào),2008,25(6):78-81.

    Slope Stability Analysis Based on Elastic-plastic Interval Finite Element

    YU Heping, YUAN Mingming, ZHANG Cong, CHEN Xiyuan, ZHANG Jiade

    (SchoolofHydraulicEngineering,ChangshaUniversityofScienceandTechnology,Changsha,Hunan410114,China)

    Abstract:Interval analysis is one of the methods to analyze the uncertain problem in engineering. The interval parameter perturbation method is used to solve the existing interval analysis method, which has obvious limitations. It is difficult to get effective solutions for the elastic plastic problem using incremental method. Interval variables are used to describe the uncertainty in the slope. In this research the interval mathematics and elastic-plastic finite element method is adopted together to establish the elastic-plastic interval finite element method. The calculation formula of Interval safety factor of slope stability is derived based on the elastic-plastic interval finite element method. This method can effectively avoid the limitation of interval parameter perturbation method. The feasibility and effectiveness of the proposed method were then verified by engineering examples. The influence of the parameters of the rock and soil on the safety factor was discussed. The results show that the inner friction angle has sensitive impact on the safety factor of the slope. The effect of elastic modulus and Poisson's ratio on the safety factor is smaller than that of cohesion and unit weight.

    Keywords:interval analysis; elastic-plastic interval finite element method; slope stability; interval safety factor

    中圖分類號(hào):TU352.1

    文獻(xiàn)標(biāo)識(shí)碼:A

    文章編號(hào):1672—1144(2016)02—0132—04

    作者簡(jiǎn)介:喻和平(1972—),女,湖南長(zhǎng)沙人,副教授,主要從事數(shù)值計(jì)算與分析方面的教學(xué)與科研工作。E-mail:413702004@qq.com通訊作者:袁明明(1990—),男,湖南郴州人,碩士研究生,研究方向?yàn)閿?shù)值計(jì)算與分析。E-mail:634432395@qq.com。

    基金項(xiàng)目:湖南省教育廳資助項(xiàng)目“區(qū)間分析理論及其在邊坡工程中的應(yīng)用”(09C086)

    收稿日期:2015-11-26修稿日期:2015-12-30

    DOI:10.3969/j.issn.1672-1144.2016.02.026

    久久久午夜欧美精品| 欧美成人免费av一区二区三区| 精品一区二区三区av网在线观看| 黄色一级大片看看| 国产一区二区在线观看日韩| 男女之事视频高清在线观看| 老司机福利观看| 欧美日本视频| 久久欧美精品欧美久久欧美| 嫩草影院精品99| 在线免费十八禁| 国产色爽女视频免费观看| 看免费成人av毛片| 91av网一区二区| 亚洲性夜色夜夜综合| 人妻少妇偷人精品九色| 波多野结衣巨乳人妻| 人妻夜夜爽99麻豆av| 免费av毛片视频| 久久6这里有精品| 国产69精品久久久久777片| 午夜视频国产福利| 午夜福利在线在线| 亚洲最大成人av| 乱系列少妇在线播放| 国产色婷婷99| 大型黄色视频在线免费观看| or卡值多少钱| 久久久国产成人精品二区| 可以在线观看毛片的网站| 精品一区二区三区视频在线观看免费| 亚洲av免费在线观看| 日日干狠狠操夜夜爽| 夜夜夜夜夜久久久久| 赤兔流量卡办理| 国产女主播在线喷水免费视频网站 | 亚洲av中文av极速乱 | 人妻制服诱惑在线中文字幕| 91午夜精品亚洲一区二区三区 | 国产免费男女视频| 国产一区二区在线观看日韩| a在线观看视频网站| 色尼玛亚洲综合影院| 伦精品一区二区三区| 久久精品国产亚洲av涩爱 | 久久久久久伊人网av| 欧美日韩精品成人综合77777| 午夜a级毛片| 国产精品,欧美在线| 国产一区二区三区在线臀色熟女| 欧美日韩乱码在线| 啪啪无遮挡十八禁网站| 99riav亚洲国产免费| 麻豆国产av国片精品| 日本撒尿小便嘘嘘汇集6| 色尼玛亚洲综合影院| 免费在线观看日本一区| 国产精品亚洲美女久久久| 国产三级中文精品| 超碰av人人做人人爽久久| 在线看三级毛片| 色综合婷婷激情| 欧美+亚洲+日韩+国产| 亚洲国产欧美人成| 一级av片app| 欧美激情国产日韩精品一区| 18禁黄网站禁片免费观看直播| 亚洲av一区综合| 欧美丝袜亚洲另类 | 一个人观看的视频www高清免费观看| 色av中文字幕| 国产熟女欧美一区二区| 他把我摸到了高潮在线观看| 精品人妻熟女av久视频| av专区在线播放| 老熟妇仑乱视频hdxx| 琪琪午夜伦伦电影理论片6080| 人妻少妇偷人精品九色| 日韩在线高清观看一区二区三区 | 1024手机看黄色片| 免费av观看视频| 又粗又爽又猛毛片免费看| a级毛片免费高清观看在线播放| 人妻少妇偷人精品九色| 特级一级黄色大片| 搡老岳熟女国产| 99久久无色码亚洲精品果冻| 精品久久久久久久久久免费视频| 亚洲精品亚洲一区二区| 免费看美女性在线毛片视频| 亚洲国产精品久久男人天堂| 永久网站在线| 露出奶头的视频| 韩国av在线不卡| 在线播放无遮挡| 日本与韩国留学比较| 国产伦人伦偷精品视频| 级片在线观看| 中文字幕av成人在线电影| 琪琪午夜伦伦电影理论片6080| 久久精品国产亚洲网站| 色综合色国产| 国产精品人妻久久久影院| 欧美成人一区二区免费高清观看| 真实男女啪啪啪动态图| 日韩人妻高清精品专区| 欧美日韩乱码在线| 亚洲天堂国产精品一区在线| 国产精品综合久久久久久久免费| 中文字幕熟女人妻在线| 国产精品久久电影中文字幕| 国产一级毛片七仙女欲春2| 国产国拍精品亚洲av在线观看| 好男人在线观看高清免费视频| 免费在线观看日本一区| 亚洲精品日韩av片在线观看| 一边摸一边抽搐一进一小说| 精品久久久噜噜| 成人毛片a级毛片在线播放| 国国产精品蜜臀av免费| 国产精品久久电影中文字幕| 成年女人看的毛片在线观看| or卡值多少钱| 级片在线观看| 久久欧美精品欧美久久欧美| 成人无遮挡网站| 免费观看精品视频网站| 成人三级黄色视频| 夜夜爽天天搞| av在线老鸭窝| 免费看a级黄色片| 91在线精品国自产拍蜜月| 99热这里只有精品一区| 亚洲成人免费电影在线观看| 国产成人福利小说| 如何舔出高潮| 熟女电影av网| 色哟哟·www| 桃色一区二区三区在线观看| 小说图片视频综合网站| 丰满乱子伦码专区| 免费黄网站久久成人精品| 99视频精品全部免费 在线| 亚洲国产欧洲综合997久久,| 又爽又黄无遮挡网站| 白带黄色成豆腐渣| 国产黄a三级三级三级人| 成人鲁丝片一二三区免费| 少妇高潮的动态图| 国产男人的电影天堂91| 国产v大片淫在线免费观看| 网址你懂的国产日韩在线| 亚洲熟妇中文字幕五十中出| 亚洲av不卡在线观看| 国产精品久久久久久久电影| 99九九线精品视频在线观看视频| 亚洲国产精品久久男人天堂| 国产av在哪里看| 久久久久久久精品吃奶| 日韩一区二区视频免费看| av中文乱码字幕在线| 女人十人毛片免费观看3o分钟| av在线老鸭窝| 国产午夜精品久久久久久一区二区三区 | 少妇的逼好多水| 日本撒尿小便嘘嘘汇集6| 国产白丝娇喘喷水9色精品| 日韩欧美三级三区| 如何舔出高潮| 97碰自拍视频| 亚洲国产精品成人综合色| 又爽又黄a免费视频| 一区二区三区高清视频在线| 亚洲精品影视一区二区三区av| 啦啦啦韩国在线观看视频| av天堂在线播放| 久久久久久久精品吃奶| 人人妻人人看人人澡| 午夜福利欧美成人| 日韩av在线大香蕉| 日本一二三区视频观看| 欧美黑人巨大hd| 日韩高清综合在线| x7x7x7水蜜桃| 久久久精品大字幕| 女同久久另类99精品国产91| 少妇丰满av| 两个人视频免费观看高清| 女人被狂操c到高潮| 特级一级黄色大片| 午夜福利18| or卡值多少钱| 久久精品国产亚洲网站| 在线国产一区二区在线| 亚洲专区中文字幕在线| АⅤ资源中文在线天堂| 国产人妻一区二区三区在| 久久久久久久亚洲中文字幕| 熟女电影av网| 精品乱码久久久久久99久播| 中出人妻视频一区二区| ponron亚洲| 亚洲国产精品成人综合色| 国产探花在线观看一区二区| 国产高清激情床上av| 午夜精品久久久久久毛片777| 亚洲最大成人中文| 色综合站精品国产| 日韩欧美在线二视频| 精品人妻偷拍中文字幕| 午夜老司机福利剧场| 赤兔流量卡办理| 久久精品久久久久久噜噜老黄 | av在线观看视频网站免费| 嫩草影院入口| 中文资源天堂在线| 日本精品一区二区三区蜜桃| 婷婷丁香在线五月| 此物有八面人人有两片| 久久久久九九精品影院| 国内精品一区二区在线观看| 黄色视频,在线免费观看| 国产综合懂色| 久久精品国产自在天天线| 能在线免费观看的黄片| 欧美区成人在线视频| 美女xxoo啪啪120秒动态图| 中文在线观看免费www的网站| 成人国产麻豆网| 少妇猛男粗大的猛烈进出视频 | 联通29元200g的流量卡| 99久久精品热视频| 精品日产1卡2卡| 我的老师免费观看完整版| 校园人妻丝袜中文字幕| 日日摸夜夜添夜夜添小说| 在线观看一区二区三区| 国产精品一区www在线观看 | 一级黄片播放器| 欧美最新免费一区二区三区| 两个人视频免费观看高清| 国产欧美日韩精品一区二区| 欧美黑人欧美精品刺激| 欧美国产日韩亚洲一区| 精品人妻偷拍中文字幕| 男人和女人高潮做爰伦理| 黄色配什么色好看| 亚洲性夜色夜夜综合| 直男gayav资源| 欧美日本亚洲视频在线播放| 可以在线观看的亚洲视频| 一个人看的www免费观看视频| 琪琪午夜伦伦电影理论片6080| 又爽又黄a免费视频| 亚洲真实伦在线观看| 久久九九热精品免费| 久久精品国产亚洲av涩爱 | 美女高潮喷水抽搐中文字幕| 亚洲美女搞黄在线观看 | 少妇裸体淫交视频免费看高清| 精品久久久久久久久久免费视频| 一进一出好大好爽视频| 成年女人毛片免费观看观看9| 国产精品一区二区性色av| 亚洲欧美日韩卡通动漫| 欧美绝顶高潮抽搐喷水| 成人特级av手机在线观看| 国产在视频线在精品| 村上凉子中文字幕在线| 最近在线观看免费完整版| 色综合色国产| 一区二区三区四区激情视频 | 精品久久久久久久久久免费视频| 久久久久久久久大av| 干丝袜人妻中文字幕| 国产爱豆传媒在线观看| 国产淫片久久久久久久久| 亚洲av二区三区四区| a级毛片a级免费在线| 观看免费一级毛片| 精品久久久久久久久亚洲 | 国产精品av视频在线免费观看| 最后的刺客免费高清国语| 久久国产精品人妻蜜桃| 国产私拍福利视频在线观看| 精品无人区乱码1区二区| 日日摸夜夜添夜夜添av毛片 | 欧美+日韩+精品| 国产av不卡久久| 久久久久久久久久久丰满 | x7x7x7水蜜桃| 国产三级中文精品| 狂野欧美白嫩少妇大欣赏| 亚洲中文字幕一区二区三区有码在线看| 欧美黑人巨大hd| 老司机深夜福利视频在线观看| 九九热线精品视视频播放| 九色国产91popny在线| 久久精品国产自在天天线| 99热6这里只有精品| 永久网站在线| 嫁个100分男人电影在线观看| 99久久九九国产精品国产免费| 日韩欧美免费精品| 国产伦一二天堂av在线观看| 中文字幕熟女人妻在线| netflix在线观看网站| 欧美精品国产亚洲| 精品久久久久久久久久免费视频| 久久这里只有精品中国| 国产亚洲av嫩草精品影院| 久久亚洲精品不卡| 91av网一区二区| 最近最新免费中文字幕在线| 欧美高清成人免费视频www| 日韩精品有码人妻一区| 噜噜噜噜噜久久久久久91| 久久久久性生活片| 91在线观看av| 3wmmmm亚洲av在线观看| 99久久精品国产国产毛片| 国产色爽女视频免费观看| 午夜福利欧美成人| 午夜激情福利司机影院| 亚洲精品一卡2卡三卡4卡5卡| 美女xxoo啪啪120秒动态图| 国产aⅴ精品一区二区三区波| 免费在线观看影片大全网站| av视频在线观看入口| 精品久久久久久久人妻蜜臀av| 国产伦精品一区二区三区视频9| 桃红色精品国产亚洲av| 国产精品一区二区免费欧美| 嫩草影视91久久| 日韩欧美精品免费久久| 精品国内亚洲2022精品成人| 22中文网久久字幕| 欧美成人性av电影在线观看| 老师上课跳d突然被开到最大视频| 九九在线视频观看精品| 国产免费av片在线观看野外av| 欧美xxxx黑人xx丫x性爽| 亚洲三级黄色毛片| 久久久久久久久大av| 国产亚洲欧美98| 亚洲成人久久性| 性插视频无遮挡在线免费观看| 啦啦啦韩国在线观看视频| 人妻夜夜爽99麻豆av| 最近视频中文字幕2019在线8| 亚洲国产精品久久男人天堂| 亚洲成人久久性| 女人十人毛片免费观看3o分钟| 婷婷丁香在线五月| 欧洲精品卡2卡3卡4卡5卡区| 变态另类成人亚洲欧美熟女| 亚洲无线在线观看| 一级av片app| 久久精品91蜜桃| 哪里可以看免费的av片| 精品一区二区三区视频在线| 亚洲中文日韩欧美视频| 丰满的人妻完整版| 很黄的视频免费| 婷婷亚洲欧美| 精华霜和精华液先用哪个| 啪啪无遮挡十八禁网站| 成人综合一区亚洲| 欧美bdsm另类| 亚洲人与动物交配视频| 免费无遮挡裸体视频| 美女高潮喷水抽搐中文字幕| 蜜桃亚洲精品一区二区三区| av.在线天堂| 在线观看美女被高潮喷水网站| 久久久久久久亚洲中文字幕| 欧美又色又爽又黄视频| 精品久久久久久久末码| 精品一区二区三区av网在线观看| 18禁黄网站禁片免费观看直播| 成熟少妇高潮喷水视频| 在线a可以看的网站| 精品免费久久久久久久清纯| 日本在线视频免费播放| 91久久精品电影网| or卡值多少钱| 日韩大尺度精品在线看网址| av国产免费在线观看| 国产人妻一区二区三区在| 在线a可以看的网站| 欧美一区二区亚洲| 国产精品久久久久久久久免| 男人舔奶头视频| 18禁黄网站禁片午夜丰满| 一区二区三区激情视频| 人妻少妇偷人精品九色| 免费看a级黄色片| 久久久国产成人精品二区| 美女免费视频网站| 免费黄网站久久成人精品| 午夜福利欧美成人| 一本精品99久久精品77| 日本精品一区二区三区蜜桃| 国产在线男女| 18禁在线播放成人免费| 伊人久久精品亚洲午夜| 搡老妇女老女人老熟妇| 日韩精品中文字幕看吧| 一区二区三区免费毛片| 色哟哟哟哟哟哟| 熟女人妻精品中文字幕| 色在线成人网| 一卡2卡三卡四卡精品乱码亚洲| 婷婷六月久久综合丁香| 精品人妻一区二区三区麻豆 | 日本色播在线视频| 亚洲avbb在线观看| 亚洲美女视频黄频| 国产又黄又爽又无遮挡在线| 99热只有精品国产| 欧美潮喷喷水| 日韩欧美免费精品| 亚洲精品一卡2卡三卡4卡5卡| 嫁个100分男人电影在线观看| 三级国产精品欧美在线观看| 狂野欧美白嫩少妇大欣赏| 国产三级在线视频| 人妻久久中文字幕网| eeuss影院久久| 国产av一区在线观看免费| 深夜精品福利| 在线国产一区二区在线| а√天堂www在线а√下载| 成人国产综合亚洲| 日本一二三区视频观看| 嫩草影视91久久| 又爽又黄无遮挡网站| 精品国产三级普通话版| 天堂网av新在线| 亚洲无线在线观看| 人妻制服诱惑在线中文字幕| 国模一区二区三区四区视频| 91久久精品电影网| 久久精品国产99精品国产亚洲性色| 亚洲成人中文字幕在线播放| 午夜久久久久精精品| 亚洲自偷自拍三级| 在线播放无遮挡| 欧美激情在线99| 熟女电影av网| 午夜激情福利司机影院| 男人舔女人下体高潮全视频| 亚洲在线自拍视频| 亚洲性夜色夜夜综合| 成人国产综合亚洲| 午夜精品在线福利| 嫩草影视91久久| 91精品国产九色| 99热网站在线观看| 精品国内亚洲2022精品成人| 免费观看人在逋| 欧美最新免费一区二区三区| 国产真实乱freesex| 日本与韩国留学比较| 免费黄网站久久成人精品| 亚洲av免费在线观看| 人人妻人人看人人澡| 国产高清有码在线观看视频| 国产视频一区二区在线看| 1024手机看黄色片| 精品久久久久久,| 搡女人真爽免费视频火全软件 | 国产av麻豆久久久久久久| 欧美一区二区亚洲| 在线a可以看的网站| 国产淫片久久久久久久久| 免费看a级黄色片| 日韩欧美在线二视频| 免费高清视频大片| av在线亚洲专区| 性插视频无遮挡在线免费观看| 精品一区二区三区人妻视频| 性欧美人与动物交配| 免费av毛片视频| 99久久精品一区二区三区| 一夜夜www| 亚洲成人免费电影在线观看| 久久精品国产亚洲网站| 国产免费av片在线观看野外av| 联通29元200g的流量卡| 69人妻影院| 精品不卡国产一区二区三区| 人妻少妇偷人精品九色| 3wmmmm亚洲av在线观看| 成人综合一区亚洲| 一夜夜www| 亚洲国产高清在线一区二区三| 亚洲人成伊人成综合网2020| 九九久久精品国产亚洲av麻豆| 日日夜夜操网爽| 亚洲av一区综合| 又黄又爽又免费观看的视频| 国产欧美日韩一区二区精品| 欧美日本亚洲视频在线播放| 中文在线观看免费www的网站| 综合色av麻豆| 欧美激情国产日韩精品一区| 男插女下体视频免费在线播放| 久久热精品热| 老司机福利观看| 欧美精品国产亚洲| 国产精品乱码一区二三区的特点| 亚洲精品影视一区二区三区av| 成人毛片a级毛片在线播放| 亚洲精品456在线播放app | 精品久久久久久久人妻蜜臀av| 欧美zozozo另类| 国产极品精品免费视频能看的| 日本黄色片子视频| 香蕉av资源在线| 亚洲av中文字字幕乱码综合| 久久亚洲真实| 日韩人妻高清精品专区| 日韩欧美在线二视频| 麻豆国产97在线/欧美| 欧美绝顶高潮抽搐喷水| 在现免费观看毛片| 亚洲精品成人久久久久久| 男女那种视频在线观看| 一个人免费在线观看电影| 在线观看免费视频日本深夜| 99热这里只有是精品在线观看| 久久久久久伊人网av| 亚洲成av人片在线播放无| 在线观看午夜福利视频| 人人妻人人澡欧美一区二区| 久久国产精品人妻蜜桃| 欧美日韩乱码在线| 国内精品美女久久久久久| 老司机午夜福利在线观看视频| 波多野结衣高清作品| 日韩,欧美,国产一区二区三区 | 亚洲精品日韩av片在线观看| 亚洲欧美日韩高清在线视频| av视频在线观看入口| 99久久久亚洲精品蜜臀av| 亚洲av电影不卡..在线观看| 国产一区二区激情短视频| 国产乱人伦免费视频| 午夜激情欧美在线| 日日撸夜夜添| 欧美激情国产日韩精品一区| 精品欧美国产一区二区三| 亚洲欧美日韩无卡精品| 啦啦啦观看免费观看视频高清| av视频在线观看入口| 日本免费一区二区三区高清不卡| 久久精品国产亚洲av天美| 狂野欧美激情性xxxx在线观看| 波野结衣二区三区在线| xxxwww97欧美| 国产精品人妻久久久影院| 亚洲欧美日韩东京热| 简卡轻食公司| 国产蜜桃级精品一区二区三区| 小蜜桃在线观看免费完整版高清| 日本 欧美在线| 日本五十路高清| 又黄又爽又刺激的免费视频.| 18禁在线播放成人免费| 日本-黄色视频高清免费观看| aaaaa片日本免费| 91麻豆精品激情在线观看国产| 亚洲,欧美,日韩| 久9热在线精品视频| 国内揄拍国产精品人妻在线| 国产亚洲精品综合一区在线观看| 欧美xxxx黑人xx丫x性爽| 日本与韩国留学比较| 91精品国产九色| 国产精品国产三级国产av玫瑰| 久久精品国产亚洲av涩爱 | 春色校园在线视频观看| 2021天堂中文幕一二区在线观| 亚洲在线观看片| 一区二区三区四区激情视频 | 三级国产精品欧美在线观看| 亚洲久久久久久中文字幕| 亚洲av日韩精品久久久久久密| 亚洲人成网站在线播放欧美日韩| 长腿黑丝高跟| 婷婷精品国产亚洲av| 两人在一起打扑克的视频| 国产男人的电影天堂91| 国产黄片美女视频| 成年免费大片在线观看| 美女高潮喷水抽搐中文字幕| 99久久九九国产精品国产免费| 非洲黑人性xxxx精品又粗又长| 99热只有精品国产| 婷婷亚洲欧美| 日本欧美国产在线视频| 此物有八面人人有两片| 亚洲av电影不卡..在线观看| 亚洲一级一片aⅴ在线观看| 免费观看精品视频网站| 欧美性猛交黑人性爽| 亚洲性久久影院| 欧美在线一区亚洲| 日本黄色视频三级网站网址| 国产日本99.免费观看| 99久久中文字幕三级久久日本| 亚洲成av人片在线播放无| 亚洲va在线va天堂va国产| 嫩草影视91久久|