• <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

    999久久久国产精品视频| 午夜91福利影院| 亚洲在久久综合| 中文字幕最新亚洲高清| 久久精品aⅴ一区二区三区四区 | 国产精品久久久久久久久免| 精品一区二区三区四区五区乱码 | 女的被弄到高潮叫床怎么办| 波野结衣二区三区在线| 啦啦啦视频在线资源免费观看| 免费大片黄手机在线观看| 黄色怎么调成土黄色| 成人亚洲精品一区在线观看| 色吧在线观看| 亚洲精品久久成人aⅴ小说| 亚洲欧洲国产日韩| 乱人伦中国视频| 成人国产麻豆网| 午夜福利乱码中文字幕| 亚洲欧美清纯卡通| 亚洲成人av在线免费| 老司机影院成人| 免费人妻精品一区二区三区视频| 亚洲国产毛片av蜜桃av| 国产成人免费无遮挡视频| 一级片免费观看大全| 成年av动漫网址| 亚洲欧洲精品一区二区精品久久久 | 日韩不卡一区二区三区视频在线| 亚洲,欧美精品.| 免费观看性生交大片5| 亚洲精品aⅴ在线观看| 丝袜喷水一区| 精品卡一卡二卡四卡免费| 国产一区二区激情短视频 | 国产精品国产av在线观看| 在线观看美女被高潮喷水网站| 亚洲中文av在线| 美女视频免费永久观看网站| 在线观看国产h片| 日韩av不卡免费在线播放| 乱人伦中国视频| 精品久久久精品久久久| 一区二区av电影网| 欧美激情极品国产一区二区三区| 亚洲国产欧美网| 日韩人妻精品一区2区三区| 久久精品国产亚洲av天美| 天堂8中文在线网| av网站在线播放免费| 久久久亚洲精品成人影院| 久久久久精品人妻al黑| 少妇精品久久久久久久| 国产精品熟女久久久久浪| av网站在线播放免费| 汤姆久久久久久久影院中文字幕| 精品人妻一区二区三区麻豆| 9热在线视频观看99| av福利片在线| 狠狠婷婷综合久久久久久88av| 少妇精品久久久久久久| 午夜激情av网站| 国产人伦9x9x在线观看 | 老司机亚洲免费影院| 精品午夜福利在线看| 在线观看免费日韩欧美大片| 精品亚洲成a人片在线观看| 亚洲婷婷狠狠爱综合网| 亚洲伊人久久精品综合| 国产探花极品一区二区| 天天躁日日躁夜夜躁夜夜| 中国三级夫妇交换| 狠狠婷婷综合久久久久久88av| 久久国产精品男人的天堂亚洲| 一本大道久久a久久精品| 日韩中文字幕视频在线看片| 伊人亚洲综合成人网| 亚洲精品久久成人aⅴ小说| av国产久精品久网站免费入址| 国产精品秋霞免费鲁丝片| 狠狠婷婷综合久久久久久88av| 欧美精品高潮呻吟av久久| 性色av一级| 青春草亚洲视频在线观看| 精品福利永久在线观看| 日韩,欧美,国产一区二区三区| 国产成人欧美| a级毛片黄视频| 久久99热这里只频精品6学生| 亚洲av.av天堂| 成年美女黄网站色视频大全免费| 男女无遮挡免费网站观看| 婷婷成人精品国产| 人人妻人人澡人人看| 免费高清在线观看日韩| 国产成人系列免费观看| 国产乱人伦免费视频| 黄色丝袜av网址大全| 18美女黄网站色大片免费观看| 神马国产精品三级电影在线观看 | 搡老熟女国产l中国老女人| 成人三级做爰电影| 国产精品 国内视频| 搡老岳熟女国产| 一本大道久久a久久精品| 成人18禁在线播放| 青草久久国产| 国产一区二区激情短视频| 精品人妻1区二区| 可以在线观看毛片的网站| 国产精品亚洲av一区麻豆| 中文字幕另类日韩欧美亚洲嫩草| 国产又爽黄色视频| 三上悠亚av全集在线观看| 高清欧美精品videossex| 亚洲精品久久午夜乱码| 午夜精品在线福利| 新久久久久国产一级毛片| 激情视频va一区二区三区| 夫妻午夜视频| 日韩大码丰满熟妇| 午夜福利在线免费观看网站| 成人三级黄色视频| 中文字幕av电影在线播放| 亚洲国产中文字幕在线视频| 真人一进一出gif抽搐免费| 久久人妻熟女aⅴ| 国产一区在线观看成人免费| 丁香六月欧美| 50天的宝宝边吃奶边哭怎么回事| 欧美日韩福利视频一区二区| 亚洲av第一区精品v没综合| avwww免费| 一区二区三区激情视频| 欧美激情久久久久久爽电影 | 国产视频一区二区在线看| av电影中文网址| 欧美中文综合在线视频| 国产一区二区在线av高清观看| 国产xxxxx性猛交| av免费在线观看网站| 欧洲精品卡2卡3卡4卡5卡区| 色婷婷av一区二区三区视频| 亚洲精品中文字幕一二三四区| 90打野战视频偷拍视频| 在线观看一区二区三区| 黑人巨大精品欧美一区二区蜜桃| 亚洲熟女毛片儿| 一级毛片高清免费大全| 多毛熟女@视频| 一个人观看的视频www高清免费观看 | 国产精品国产高清国产av| 亚洲欧美日韩另类电影网站| 亚洲五月色婷婷综合| 久久香蕉精品热| 日韩免费av在线播放| 久久人人精品亚洲av| 国产精品免费一区二区三区在线| 亚洲成人国产一区在线观看| 91精品三级在线观看| 婷婷六月久久综合丁香| 黑人巨大精品欧美一区二区蜜桃| 一区二区三区激情视频| 亚洲精品国产色婷婷电影| 91成年电影在线观看| 一进一出抽搐gif免费好疼 | 热re99久久国产66热| 久久草成人影院| 欧美乱妇无乱码| 大陆偷拍与自拍| 久久久国产欧美日韩av| 两性午夜刺激爽爽歪歪视频在线观看 | 久久人妻av系列| 69精品国产乱码久久久| 亚洲性夜色夜夜综合| 亚洲视频免费观看视频| 一进一出抽搐gif免费好疼 | 老司机深夜福利视频在线观看| 久久 成人 亚洲| 亚洲五月天丁香| 19禁男女啪啪无遮挡网站| 巨乳人妻的诱惑在线观看| 免费在线观看影片大全网站| ponron亚洲| 搡老熟女国产l中国老女人| 亚洲国产中文字幕在线视频| 亚洲精品中文字幕一二三四区| 黑人巨大精品欧美一区二区蜜桃| 午夜激情av网站| 人人妻人人添人人爽欧美一区卜| 久久久久精品国产欧美久久久| av超薄肉色丝袜交足视频| 99精品欧美一区二区三区四区| 欧美日韩精品网址| 国产99久久九九免费精品| 国产亚洲欧美98| 中文欧美无线码| 免费高清视频大片| 黄色女人牲交| 伊人久久大香线蕉亚洲五| 亚洲一区二区三区不卡视频| 国产日韩一区二区三区精品不卡| 日韩大码丰满熟妇| 久久中文字幕人妻熟女| 97碰自拍视频| 日本三级黄在线观看| 亚洲国产欧美一区二区综合| 午夜福利一区二区在线看| 久久国产精品影院| 两个人看的免费小视频| 国产成人一区二区三区免费视频网站| 一二三四社区在线视频社区8| 亚洲第一av免费看| 午夜精品国产一区二区电影| 国产亚洲精品一区二区www| 正在播放国产对白刺激| 精品一区二区三区四区五区乱码| 这个男人来自地球电影免费观看| 午夜福利免费观看在线| 97碰自拍视频| 久久国产精品影院| 99国产精品99久久久久| 99香蕉大伊视频| 亚洲av日韩精品久久久久久密| 在线观看一区二区三区激情| 午夜视频精品福利| 欧美乱色亚洲激情| 日韩免费高清中文字幕av| 成人影院久久| 18禁裸乳无遮挡免费网站照片 | 久久久久久久午夜电影 | 成人黄色视频免费在线看| 成人亚洲精品一区在线观看| 夜夜躁狠狠躁天天躁| 精品久久久久久成人av| 亚洲专区字幕在线| 日韩中文字幕欧美一区二区| 老熟妇仑乱视频hdxx| 男女午夜视频在线观看| 村上凉子中文字幕在线| 他把我摸到了高潮在线观看| 亚洲精品国产精品久久久不卡| 中文字幕人妻丝袜制服| 99国产精品一区二区蜜桃av| 中国美女看黄片| 欧美中文日本在线观看视频| 精品欧美一区二区三区在线| 啦啦啦在线免费观看视频4| 中文字幕精品免费在线观看视频| 999久久久精品免费观看国产| 国产男靠女视频免费网站| 搡老岳熟女国产| 亚洲狠狠婷婷综合久久图片| 在线看a的网站| 久久久国产成人精品二区 | 亚洲三区欧美一区| 国产99白浆流出| 午夜免费成人在线视频| 日本黄色日本黄色录像| 嫩草影院精品99| 黄片大片在线免费观看| 精品熟女少妇八av免费久了| 欧美中文日本在线观看视频| 成人18禁高潮啪啪吃奶动态图| av电影中文网址| 三上悠亚av全集在线观看| 精品一品国产午夜福利视频| 91精品三级在线观看| 一级,二级,三级黄色视频| 国产欧美日韩一区二区精品| 精品一区二区三区av网在线观看| 女人爽到高潮嗷嗷叫在线视频| 欧美日韩黄片免| 久久国产乱子伦精品免费另类| 亚洲久久久国产精品| 欧美成人性av电影在线观看| 成年人黄色毛片网站| 高清av免费在线| 亚洲av片天天在线观看| 色哟哟哟哟哟哟| 国产精品香港三级国产av潘金莲| 亚洲av日韩精品久久久久久密| 欧美在线一区亚洲| 久久久精品国产亚洲av高清涩受| 免费av中文字幕在线| 99精品在免费线老司机午夜| 又紧又爽又黄一区二区| a级片在线免费高清观看视频| 在线观看免费日韩欧美大片| aaaaa片日本免费| 国产1区2区3区精品| 成熟少妇高潮喷水视频| 在线观看舔阴道视频| 亚洲成人免费电影在线观看| 一区二区日韩欧美中文字幕| 中文字幕人妻丝袜制服| 成年人黄色毛片网站| 国产激情久久老熟女| 精品高清国产在线一区| 男人舔女人的私密视频| 18禁国产床啪视频网站| 欧美成人性av电影在线观看| 成人黄色视频免费在线看| 高清av免费在线| 在线观看舔阴道视频| 女人高潮潮喷娇喘18禁视频| 天堂√8在线中文| 亚洲国产欧美一区二区综合| 亚洲一区高清亚洲精品| 在线观看免费高清a一片| 不卡一级毛片| 男女床上黄色一级片免费看| 69av精品久久久久久| 精品久久久久久电影网| 黄片播放在线免费| 9热在线视频观看99| 村上凉子中文字幕在线| 欧美黄色片欧美黄色片| 人成视频在线观看免费观看| 成年人黄色毛片网站| 国产单亲对白刺激| 国产成+人综合+亚洲专区| 久久国产精品男人的天堂亚洲| 久久中文字幕一级| 国产区一区二久久| av国产精品久久久久影院| 久久久久亚洲av毛片大全| 亚洲中文av在线| 性欧美人与动物交配| 亚洲国产精品合色在线| 亚洲国产毛片av蜜桃av| 久久精品91无色码中文字幕| 高清av免费在线| 国产熟女xx| av欧美777| 18禁美女被吸乳视频| 自拍欧美九色日韩亚洲蝌蚪91| 一边摸一边抽搐一进一出视频| 国产在线观看jvid| 亚洲熟妇熟女久久| 午夜精品在线福利| 免费观看精品视频网站| 亚洲专区国产一区二区| 亚洲第一青青草原| 国产精品成人在线| 国产精品电影一区二区三区| 12—13女人毛片做爰片一| 一区福利在线观看| 淫秽高清视频在线观看| aaaaa片日本免费| 99国产精品一区二区三区| av在线播放免费不卡| 91字幕亚洲| 男女做爰动态图高潮gif福利片 | 中文字幕精品免费在线观看视频| 神马国产精品三级电影在线观看 | 久久久久久大精品| aaaaa片日本免费| 成年版毛片免费区| 久久天堂一区二区三区四区| 如日韩欧美国产精品一区二区三区| 免费观看精品视频网站| 亚洲全国av大片| 免费女性裸体啪啪无遮挡网站| 操出白浆在线播放| 午夜91福利影院| 亚洲全国av大片| 美女高潮到喷水免费观看| 黄色片一级片一级黄色片| 国产免费男女视频| 午夜福利一区二区在线看| 一区在线观看完整版| www日本在线高清视频| 午夜免费观看网址| 欧美精品一区二区免费开放| 韩国av一区二区三区四区| 深夜精品福利| 性少妇av在线| 国产激情欧美一区二区| 91成人精品电影| 欧美激情高清一区二区三区| 国产成人欧美在线观看| 欧美日韩精品网址| 国产成年人精品一区二区 | 久久香蕉精品热| 亚洲熟妇中文字幕五十中出 | 乱人伦中国视频| 电影成人av| 亚洲欧美日韩高清在线视频| 神马国产精品三级电影在线观看 | 操出白浆在线播放| 久久久久久久久中文| 黄色女人牲交| 久久久国产成人免费| 女人爽到高潮嗷嗷叫在线视频| 久久国产精品男人的天堂亚洲| 黄片大片在线免费观看| 91老司机精品| 亚洲成人免费av在线播放| 黄片大片在线免费观看| 搡老熟女国产l中国老女人| 午夜免费观看网址| aaaaa片日本免费| 欧美成人午夜精品| 香蕉丝袜av| 国产在线精品亚洲第一网站| 亚洲 国产 在线| 精品福利永久在线观看| 一级毛片精品| 99re在线观看精品视频| 真人一进一出gif抽搐免费| 色婷婷久久久亚洲欧美| 亚洲,欧美精品.| 亚洲免费av在线视频| 妹子高潮喷水视频| 天堂中文最新版在线下载| 夜夜躁狠狠躁天天躁| 国产亚洲精品第一综合不卡| 男人操女人黄网站| 老熟妇乱子伦视频在线观看| 久久精品国产清高在天天线| 桃色一区二区三区在线观看| 日韩中文字幕欧美一区二区| 一二三四社区在线视频社区8| 国产av又大| 亚洲欧美日韩高清在线视频| 亚洲美女黄片视频| 午夜精品在线福利| 又黄又粗又硬又大视频| 一区二区日韩欧美中文字幕| 欧美日韩视频精品一区| 欧美乱码精品一区二区三区| 免费少妇av软件| 午夜精品在线福利| 欧美日韩国产mv在线观看视频| 岛国视频午夜一区免费看| 最新在线观看一区二区三区| 亚洲午夜精品一区,二区,三区| 午夜福利免费观看在线| 9色porny在线观看| 国产男靠女视频免费网站| 欧美日韩视频精品一区| 国产激情欧美一区二区| 另类亚洲欧美激情| 精品午夜福利视频在线观看一区| 天天躁夜夜躁狠狠躁躁| 1024视频免费在线观看| 美国免费a级毛片| 亚洲成av片中文字幕在线观看| 精品一区二区三区四区五区乱码| 在线观看66精品国产| 免费观看人在逋| 国产aⅴ精品一区二区三区波| 丁香欧美五月| 香蕉国产在线看| 99国产精品一区二区三区| 久久精品人人爽人人爽视色| 91老司机精品| 在线永久观看黄色视频| 女人精品久久久久毛片| 免费在线观看影片大全网站| 久久精品91无色码中文字幕| 久久人妻av系列| 老司机深夜福利视频在线观看| 欧美黄色片欧美黄色片| 亚洲精品国产一区二区精华液| 天堂俺去俺来也www色官网| 亚洲精品中文字幕一二三四区| 久久 成人 亚洲| 久久天堂一区二区三区四区| 首页视频小说图片口味搜索| 男人舔女人的私密视频| 日韩高清综合在线| 欧美精品啪啪一区二区三区| 国产亚洲精品久久久久久毛片| 琪琪午夜伦伦电影理论片6080| 国产人伦9x9x在线观看| 欧美日本亚洲视频在线播放| 欧美乱妇无乱码| 一级a爱视频在线免费观看| 精品久久久久久久毛片微露脸| 一级毛片精品| 免费高清视频大片| 精品久久久久久成人av| 欧美日本亚洲视频在线播放| 一区在线观看完整版| 久久香蕉精品热| 久久久久久久精品吃奶| 黑人操中国人逼视频| 亚洲久久久国产精品| 人妻久久中文字幕网| 亚洲伊人色综图| 久久香蕉激情| 亚洲精品成人av观看孕妇| 亚洲aⅴ乱码一区二区在线播放 | 老司机福利观看| 亚洲av成人不卡在线观看播放网| 精品福利永久在线观看| av电影中文网址| 中文字幕人妻丝袜一区二区| www.熟女人妻精品国产| 亚洲成人精品中文字幕电影 | 久久国产精品人妻蜜桃| 色在线成人网| 日本撒尿小便嘘嘘汇集6| 亚洲精品中文字幕在线视频| 男女下面插进去视频免费观看| 在线观看免费高清a一片| 日本精品一区二区三区蜜桃| 日韩欧美在线二视频| 亚洲成人精品中文字幕电影 | 精品一区二区三区av网在线观看| 亚洲国产精品999在线| 制服诱惑二区| 亚洲欧美日韩另类电影网站| 国产成人精品无人区| 9191精品国产免费久久| 99热只有精品国产| 成年女人毛片免费观看观看9| 久久久久亚洲av毛片大全| 久9热在线精品视频| 亚洲人成网站在线播放欧美日韩| 日日干狠狠操夜夜爽| 咕卡用的链子| 大陆偷拍与自拍| 看片在线看免费视频| 琪琪午夜伦伦电影理论片6080| 日韩大尺度精品在线看网址 | 成人三级黄色视频| 一进一出抽搐动态| 成人av一区二区三区在线看| 一个人观看的视频www高清免费观看 | 性欧美人与动物交配| 国产无遮挡羞羞视频在线观看| 午夜精品在线福利| 亚洲成人免费电影在线观看| a级毛片在线看网站| 97人妻天天添夜夜摸| 亚洲av第一区精品v没综合| 日韩免费高清中文字幕av| videosex国产| 黄片播放在线免费| 老汉色av国产亚洲站长工具| 亚洲精品av麻豆狂野| 欧美激情极品国产一区二区三区| 人人澡人人妻人| 国产91精品成人一区二区三区| 在线播放国产精品三级| 午夜福利欧美成人| 高清在线国产一区| 欧美老熟妇乱子伦牲交| 国产高清videossex| 午夜a级毛片| 男女高潮啪啪啪动态图| 91精品三级在线观看| 亚洲精品av麻豆狂野| 国产精品av久久久久免费| 欧美 亚洲 国产 日韩一| 大型av网站在线播放| 在线观看免费视频日本深夜| 成人手机av| 日韩精品青青久久久久久| 亚洲欧美激情综合另类| 精品人妻1区二区| 国产精品秋霞免费鲁丝片| 韩国精品一区二区三区| 日本一区二区免费在线视频| www.自偷自拍.com| 极品人妻少妇av视频| 天天躁狠狠躁夜夜躁狠狠躁| 操美女的视频在线观看| 精品久久久久久,| 性欧美人与动物交配| 国产在线观看jvid| 午夜成年电影在线免费观看| 国产又色又爽无遮挡免费看| 亚洲精品国产一区二区精华液| 国产av在哪里看| 首页视频小说图片口味搜索| 在线视频色国产色| 很黄的视频免费| 最好的美女福利视频网| 日本wwww免费看| 99国产极品粉嫩在线观看| av片东京热男人的天堂| 国产在线观看jvid| 两人在一起打扑克的视频| 日韩欧美国产一区二区入口| 午夜老司机福利片| 国产深夜福利视频在线观看| 亚洲成a人片在线一区二区| 最好的美女福利视频网| 可以免费在线观看a视频的电影网站| 欧美日本亚洲视频在线播放| 亚洲一区高清亚洲精品| 天天躁狠狠躁夜夜躁狠狠躁| 成人亚洲精品av一区二区 | 午夜精品久久久久久毛片777| 一进一出抽搐gif免费好疼 | 中文亚洲av片在线观看爽| 在线观看免费日韩欧美大片| 热re99久久精品国产66热6| 国产欧美日韩综合在线一区二区| 黑人巨大精品欧美一区二区蜜桃| 免费不卡黄色视频| 久久久国产欧美日韩av| 老司机福利观看| 国产精品久久久av美女十八| 午夜成年电影在线免费观看| 91麻豆精品激情在线观看国产 | 欧美黑人精品巨大| 黑人欧美特级aaaaaa片| 国产精品日韩av在线免费观看 | 精品日产1卡2卡| 一进一出好大好爽视频|