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

    一種將線性粘彈微分型本構(gòu)方程應(yīng)用到ABAQUS的方法

    2010-09-07 07:32:04潘曉明余俊楊釗孔娟
    關(guān)鍵詞:本構(gòu)微分管片

    潘曉明,余俊,楊釗,孔娟

    (1.同濟(jì)大學(xué)巖土及地下工程教育部重點(diǎn)實(shí)驗(yàn)室,上海200092; 2.中南大學(xué)土木建筑學(xué)院,湖南長(zhǎng)沙410075)

    一種將線性粘彈微分型本構(gòu)方程應(yīng)用到ABAQUS的方法

    潘曉明1,余俊2,楊釗1,孔娟1

    (1.同濟(jì)大學(xué)巖土及地下工程教育部重點(diǎn)實(shí)驗(yàn)室,上海200092; 2.中南大學(xué)土木建筑學(xué)院,湖南長(zhǎng)沙410075)

    以3參量Merchant模型、Poyning-Thomson模型,以及4參量Burger模型為例,推導(dǎo)3種常用微分型本構(gòu)方程的松弛剪切模量Prony級(jí)數(shù)形式,并確定其在ABAQUS中的輸入?yún)?shù),供ABAQUS調(diào)用分析.以算例分析巖土流變性態(tài),將3種模型的應(yīng)變?nèi)渥?、松弛?shù)值解與一維情形下的解析解進(jìn)行對(duì)比,結(jié)果吻合較好.對(duì)土體采用Merchant粘彈性流變模型的過(guò)江盾構(gòu)隧道進(jìn)行粘彈性數(shù)值分析,得到管片位移均趨于穩(wěn)定值,說(shuō)明所提出的方法是正確和可靠的,可以應(yīng)用于巖土工程數(shù)值計(jì)算.

    微分型本構(gòu);Lap lace變換;Prony級(jí)數(shù);蠕變;松弛;ABAQUS

    對(duì)于巖土所呈現(xiàn)的各種不同的流變特性,一般采用3種方法進(jìn)行描述.(1)通過(guò)室內(nèi)或現(xiàn)場(chǎng)試驗(yàn)和實(shí)測(cè)[1-3],建立經(jīng)驗(yàn)流變或蠕變公式;(2)建立積分型本構(gòu)方程求解流變問(wèn)題[4-6];(3)根據(jù)流變介質(zhì)模型,建立微分型本構(gòu)方程,對(duì)巖土體流變進(jìn)行分析.第3種方法直觀、形象,且便于有限元編程計(jì)算,因而得到廣泛的應(yīng)用[7-9].ABAQUS是目前國(guó)際上最為先進(jìn)的大型通用有限元計(jì)算分析軟件之一.它所提供的粘彈本構(gòu)模型是基于松弛剪切模量Prony級(jí)數(shù)形式的積分型本構(gòu)方程[10],這給傳統(tǒng)巖土體線性粘彈微分型本構(gòu)方程的應(yīng)用帶來(lái)了不便.本文研究一種將線性粘彈性微分型本構(gòu)方程轉(zhuǎn)化成松弛剪切模量的Prony級(jí)數(shù)形式的方法,以供ABAQUS調(diào)用分析.

    1 線性粘彈微分型本構(gòu)方程轉(zhuǎn)換成Prony級(jí)數(shù)形式

    巖土流變模型是以反映蠕變、應(yīng)力松弛等為前提,利用彈性元件(胡克體)、粘性元件(牛頓體),通過(guò)串并聯(lián)關(guān)系形成的各種流變模型.文中討論的3參量Merchant模型、Poyning-Thom son模型,以及4參量Burger模型,分別如圖1所示.

    圖1 流變模型Fig.1 Rheologicalmodel

    1.1 微分型本構(gòu)方程的Laplace變換[11-12]

    線性粘彈性微分型本構(gòu)關(guān)系可以表示為1.3 Merchant模型

    M erchant模型的一維本構(gòu)方程為

    通過(guò)同樣的方法,可得到Merchant模型在ABAQUS中的輸入?yún)?shù)為

    1.4 Poyn ing-Thom son模型

    Poyning-Thom son模型的一維本構(gòu)方程為

    通過(guò)同樣的方法,可得到Poyning-Thom son模型在ABAQUS中的輸入?yún)?shù)為

    圖2 三維計(jì)算模型Fig.2 Three-dimensional computationalmodel

    2 算例驗(yàn)證

    驗(yàn)證算例為一個(gè)三維C3D8單元試件,如圖2所示.件長(zhǎng)、寬、高均為1 m,試件底部固定,頂部受均布?jí)毫進(jìn)行蠕變?cè)囼?yàn)數(shù)值模擬,以及施加恒定位移進(jìn)行松弛試驗(yàn)數(shù)值模擬.材料基本參數(shù):EH為2.0 M Pa,泊松比μ為0.25,EM為2.0 M Pa,ηM為21.6 M Pa·月,EK為2.0 M Pa,ηK為21.6 M Pa·月.

    通過(guò)式(13)~(17),將微分型本構(gòu)Burger模型、M erchant模型及Poyning-Thom son模型參數(shù)轉(zhuǎn)換為可供ABAQUS調(diào)用的輸入?yún)?shù).

    (1)Burger模型的一維蠕變和一維松弛解析解為

    由此可以得到Burger模型在ABAQUS中的輸入?yún)?shù):E0=EM=2.0 M Pa,μ0=μM=0.25,g1= 0.768,g2=0.231,τ1=3.59,τ2=27.0.

    (2)Merchant模型的一維蠕變和一維松弛解析解為

    由此可以得到M erchant模型在ABAQUS中的輸入?yún)?shù):E0=EH=2.0 M Pa,μ0=μH=0.25,g1= 0.545 454 5,τ1=4.91.

    (3)Poyning-Thom son模型的一維蠕變和一維松弛解析解為

    由此可得到Poyning-Thom son模型在ABAQUS中的輸入?yún)?shù),即E0=EH+EM=4.0 M Pa,μ0= μH=0.25,g1=0.454 545 45,τ1=10.8.

    在ABAQUS中,應(yīng)力以受壓為負(fù),受拉為正.在施加荷載的中,均布荷載p為0.1 M Pa,進(jìn)行蠕變數(shù)值試驗(yàn);恒定位移u3為0.3 m,進(jìn)行松弛數(shù)值試驗(yàn).

    對(duì)于流變計(jì)算,在ABAQUS中一般設(shè)置兩個(gè)分析步.第1個(gè)分析步為Static靜力分析步,時(shí)間設(shè)置很小,如1×10-20s;第2個(gè)分析步為V isco準(zhǔn)靜態(tài)分析步,采用固定時(shí)間增量步為0.5月,計(jì)算總時(shí)間為150月.

    在均布荷載作用下,Burger模型、Merchant模型及Poyning-Thom son模型的應(yīng)變?nèi)渥償?shù)值解與一維情形下蠕變解析解,如圖3所示.從圖3可知,Burger模型能夠較好地模擬蠕變的第2階段,其軸向應(yīng)變?chǔ)烹S時(shí)間(t)的增加而增加;M erchant模型和Poyning-Thom son模型能夠較好地模擬蠕變的第1階段,軸向應(yīng)變隨時(shí)間的增加而趨于穩(wěn)定值.

    在常應(yīng)變作用下,3種模型應(yīng)力σ松弛數(shù)值解與一維松弛解析解,如圖4所示.從圖4可知,Buger模型存在應(yīng)力松弛,隨著時(shí)間的增加最終趨于零;Merchant模型和Poyning-Thom son模型具有松弛特性,隨時(shí)間的增加最終趨于穩(wěn)定值.

    從圖3,4可以看出,3種模型的蠕變及松弛的解析解與數(shù)值解吻合較好,蠕變最大相對(duì)誤差不超過(guò)1.34%,松弛最大誤差不超過(guò)1.28%.

    圖4 模型軸向應(yīng)力解析解與數(shù)值解比較 Fig.4 Comparison between analytical and numerical axial stresses

    圖3 模型軸向應(yīng)變解析解與數(shù)值解比較Fig.3 Comparison between analytical and numerical axial strains

    3 工程實(shí)例

    某過(guò)江隧道襯砌外徑外直徑為8.7 m,管片內(nèi)直徑為7.9 m,管片厚為400 mm;有限元模型尺寸長(zhǎng)60 m,寬79.38 m,拱頂埋深17.6 m,其網(wǎng)格圖如圖5所示.坐標(biāo)系采用水平向右為X方向,垂直向上為Y方向.底部邊界全約束邊界條件,左、右側(cè)邊界則采用X軸方向?qū)ΨQ約束,上部邊界采用自由邊界.

    圖5 有限元計(jì)算網(wǎng)格Fig.5 Finite elementmesh

    土體單元選用雙線性四邊形完全積分單元,單元數(shù)5 275個(gè),節(jié)點(diǎn)數(shù)5 525個(gè).管片采用實(shí)體單元,選用能考慮彎曲并可克服剪切自鎖的四邊形雙線性非協(xié)調(diào)等參單元.管片單元數(shù)168個(gè),節(jié)點(diǎn)數(shù)252個(gè).土體采用M erchant本構(gòu)模型,試驗(yàn)參數(shù)和經(jīng)轉(zhuǎn)化為ABAQUS的粘彈輸入?yún)?shù),如表1所示.

    隧洞中心處的外水壓力為0.367 4 M Pa,拱頂至拱底的側(cè)向水壓力呈梯形分布.對(duì)砂性土采用水土分算,外水壓力直接作用在管片襯砌外層周?chē)?計(jì)算時(shí)設(shè)置4個(gè)分析步:(1)設(shè)置Geostatic重力分析步.計(jì)算后將土體位移場(chǎng)清零,只保留應(yīng)力場(chǎng);(2)設(shè)置為Static靜力分析步.開(kāi)挖隧道范圍內(nèi)土體并施加管片,計(jì)算時(shí)間為1× 10-20s;(3)設(shè)置為Visco準(zhǔn)靜態(tài)分析步.施加管片重力及外水壓,計(jì)算時(shí)間為36月,時(shí)間增量為0.5月.

    在第3個(gè)分析步開(kāi)始時(shí),管片在施加重力荷載及施加外水壓作用下,瞬時(shí)產(chǎn)生的彈性變形如圖6所示.從圖6可知,由于管片接頭位置不對(duì)稱,造成水平位移不對(duì)稱.在土體采用粘彈材料,計(jì)算時(shí)間為36月時(shí),管片的位移云圖如圖7所示.從圖7可知,管片水平相對(duì)位移由24.1 mm變化為17.8 mm;管片豎向相對(duì)位移由25.4 mm變化為19.5 mm.這主要是由于土體發(fā)生蠕變作用,使得管片有向內(nèi)收斂.管片拱頂、拱底豎向位移,以及右側(cè)拱腰水平位移隨時(shí)間的變化曲線,如圖8所示.從圖8可看出,隨著時(shí)間的增加,管片在不同位置處的位移均達(dá)到穩(wěn)定的收斂值,這與一維情形下M erchant模型的蠕變特性相吻合.拱頂由瞬時(shí)彈性豎向位移-0.08 mm,經(jīng)14月穩(wěn)定至-16.2 mm;拱底豎向位移由25.3 mm變化至0.4 mm;右側(cè)拱腰水平位移由13.5 mm變化至10.7 mm.

    表1 土體計(jì)算參數(shù)表Tab.1 Computational parametersof soil

    圖6 Visco準(zhǔn)靜態(tài)分析步開(kāi)始時(shí)的管片瞬時(shí)彈性位移Fig.6 Instantaneous elastic disp lacement of segment at the beginning of step 3

    圖7 Visco準(zhǔn)靜態(tài)分析步結(jié)束時(shí)的管片粘彈位移Fig.7 Visco-elastic disp lacement of segment at the end of step 3

    4 結(jié)論

    圖8 管片位移隨時(shí)間的變化曲線Fig.8 Curves of disp lacement to time at different location of segment

    研究將線性粘彈微分型本構(gòu)方程應(yīng)用于ABAQUS分析平臺(tái)的方法.推導(dǎo)Burger模型,Merchant模型及Poyning-Thom son線性粘彈微分型本構(gòu)方程的松弛剪切模量p rony級(jí)數(shù)形式,并確定其在ABAQUS中的輸入?yún)?shù).

    通過(guò)三維蠕變和松弛的數(shù)值模擬結(jié)果與一維情形下的解析解進(jìn)行對(duì)比,表明兩者的計(jì)算結(jié)果吻合較好.最后,對(duì)土體采用Merchant粘彈性流變模型的過(guò)江盾構(gòu)隧道進(jìn)行了粘彈性數(shù)值分析,得到了管片位移均趨于穩(wěn)定值,進(jìn)一步說(shuō)明所提出的方法是正確的和可靠的,可以應(yīng)用于巖土工程數(shù)值計(jì)算.

    對(duì)于將其他線性粘彈微分型本構(gòu)方程的應(yīng)用到ABAQUS中,所提出方法具有一定的參考價(jià)值.

    [1] 趙永輝.潤(rùn)揚(yáng)長(zhǎng)江大橋北錨碇土體流變特性的試驗(yàn)研究[J].地下空間,2003,23(4):417-420.

    [2] 謝寧,孫鈞.上海地區(qū)飽和軟粘土流變特性[J].同濟(jì)大學(xué)學(xué)報(bào),1996,24(3):233-237.

    [3] 劉朝豐,鮑文魁,徐景龍.剪切模量的近似函數(shù)和精確數(shù)值積分表達(dá)式[J].上海航天,2003,20(5):26-28.

    [4] 焦春茂,趙春風(fēng),張征,等.遺傳積分型本構(gòu)方程的兩種推演方法及比較[J].地下空間與工程學(xué)報(bào),2008,3(4):42-45.

    [5] 王琛,劉浩吾.線性積分模型和微分模型的普遍流變方程[J],土木工程學(xué)報(bào),2003,36(10):99-102.

    [6] THOMAS F S.Finite-element calculation of stresses in glass parts undergoing viscous relaxation[J].J Am Ceram Soc,1999,70(2):90-95.

    [7] 孫鈞.巖土材料流變及其工程應(yīng)用[M].北京:中國(guó)建筑工業(yè)出版社,1999:1-40.

    [8] 章根德,何鮮,朱維耀.巖石介質(zhì)流變學(xué)[M].北京:科學(xué)出版社,1999:20-60.

    [9] 王智祥,劉金明,劉柏雄,等.AZ61鎂合金熱壓縮流變應(yīng)力的實(shí)驗(yàn)[J].華僑大學(xué)學(xué)報(bào):自然科學(xué)版,2009,30(3): 248-252

    [10] HKS Inc.ABAQUS user’smanual[M].Paw tucket:Hibbit,Karlsson&Sorensen Inc,2002.

    [11] 王芝銀,李云鵬.巖體流變理論及其數(shù)值模擬[M].北京:科學(xué)出版社,2008:6-30.

    [12] 陳子蔭.圍巖力學(xué)分析中的解析方法[M].北京:煤炭工業(yè)出版社,1994:303-390.

    A Method Using L inear Viscoelastic D ifferen tial Constitutive Equation in to ABAQUS

    PAN Xiao-m ing1,YU Jun2, YANG Zhao1,KONG Juan1
    (1.Key Laboratory of Geotechnical and Underground Engineering of Ministry of Education,Tongji University,Shanghai 200092,China; 2.School of Civil Engineering and A rchitecture,Central South University,Changsha 410075,China)

    Taking three-parameter Merchantmodel,Poyning-Thomsonmodel and four-parameter Burgermodel as examp le,the relaxation shear modulus w ith p rony series fo rm s is derived by three kinds of commonly used differential-type constitutive equation,their parameters is determined to input in ABAQUS.In order to analyze the rheological characteristics,the creep and relaxation numerical solutionsof threemodelsagreewell to analytical solutions.Using Merchant viscoelastic rheologicalmodel,the viscoelasto numerical calculation of cross-river shield tunnel is carried out,the stabile values of segment displacement are obtained,w hich indicate the p roposed method is co rrect and reliable,and can be app lied to geotechnical engineering numerical calculation.

    differential constitutive;Laplace transform;p rony series;creep;relax;ABAQUS

    TU 452

    A

    (責(zé)任編輯:錢(qián)筠 英文審校:方德平)

    1000-5013(2010)05-0570-06

    2009-10-11

    潘曉明(1979-),男,博士研究生,主要從事隧道及地下結(jié)構(gòu)方面的研究.E-mail:p xm155138@163.com.

    國(guó)家自然科學(xué)基金資助項(xiàng)目(40872179);中國(guó)博士后科研基金資助項(xiàng)目(20080440652)

    猜你喜歡
    本構(gòu)微分管片
    擬微分算子在Hp(ω)上的有界性
    上下解反向的脈沖微分包含解的存在性
    離心SC柱混凝土本構(gòu)模型比較研究
    管片拼裝技術(shù)研究
    盾構(gòu)管片封頂塊拼裝施工技術(shù)研究
    鋸齒形結(jié)構(gòu)面剪切流變及非線性本構(gòu)模型分析
    一種新型超固結(jié)土三維本構(gòu)模型
    借助微分探求連續(xù)函數(shù)的極值點(diǎn)
    地鐵盾構(gòu)管片受力分析及管片破損的控制措施研究
    對(duì)不定積分湊微分解法的再認(rèn)識(shí)
    国产一区亚洲一区在线观看| 少妇人妻精品综合一区二区| 性色av一级| 特大巨黑吊av在线直播| 真实男女啪啪啪动态图| 国产伦精品一区二区三区四那| 欧美激情国产日韩精品一区| 亚洲欧美精品自产自拍| av在线观看视频网站免费| 久久久成人免费电影| 国产精品国产三级国产av玫瑰| 97在线视频观看| 欧美日韩视频高清一区二区三区二| 国产精品久久久久久av不卡| 亚洲精品成人av观看孕妇| 亚洲精品国产av蜜桃| 亚洲人成网站高清观看| 特级一级黄色大片| 午夜福利网站1000一区二区三区| 国产真实伦视频高清在线观看| 亚洲伊人久久精品综合| 国国产精品蜜臀av免费| 插阴视频在线观看视频| 日韩一区二区三区影片| 99热这里只有是精品50| 国产精品一区二区三区四区免费观看| 国产乱人视频| 成人国产av品久久久| 国产日韩欧美亚洲二区| 国产伦在线观看视频一区| 老司机影院毛片| 国产永久视频网站| 在线观看一区二区三区| 欧美性猛交╳xxx乱大交人| 一区二区三区免费毛片| 婷婷色av中文字幕| 国产在线男女| 99精国产麻豆久久婷婷| 欧美亚洲 丝袜 人妻 在线| 日韩精品有码人妻一区| 久久久精品欧美日韩精品| 丝袜美腿在线中文| 免费在线观看成人毛片| 国产成人精品一,二区| 不卡视频在线观看欧美| 日本-黄色视频高清免费观看| 只有这里有精品99| 亚洲国产欧美在线一区| 六月丁香七月| 国产亚洲一区二区精品| 大又大粗又爽又黄少妇毛片口| 高清视频免费观看一区二区| 交换朋友夫妻互换小说| 18禁裸乳无遮挡动漫免费视频 | av卡一久久| 男人爽女人下面视频在线观看| 日本黄色片子视频| 可以在线观看毛片的网站| 麻豆成人av视频| 最近中文字幕2019免费版| 日日啪夜夜撸| 免费不卡的大黄色大毛片视频在线观看| 深爱激情五月婷婷| 精品人妻熟女av久视频| 国产亚洲精品久久久com| 成人鲁丝片一二三区免费| 日日啪夜夜爽| 久热这里只有精品99| 久久久久久久久久久免费av| 国产极品天堂在线| 日韩av免费高清视频| 国产av国产精品国产| av女优亚洲男人天堂| 女的被弄到高潮叫床怎么办| 欧美一级a爱片免费观看看| 亚洲国产日韩一区二区| 国产精品国产三级专区第一集| 少妇人妻久久综合中文| 超碰av人人做人人爽久久| 日本一二三区视频观看| 日本三级黄在线观看| 国产精品人妻久久久久久| 99精国产麻豆久久婷婷| 国产精品av视频在线免费观看| 久久精品综合一区二区三区| 97精品久久久久久久久久精品| 国产精品久久久久久久久免| 日韩成人伦理影院| 久久久久久久久久人人人人人人| 九九在线视频观看精品| 精品久久久久久电影网| 久久6这里有精品| 黄色视频在线播放观看不卡| 亚洲三级黄色毛片| 欧美三级亚洲精品| 亚洲欧美成人综合另类久久久| 激情五月婷婷亚洲| 麻豆成人午夜福利视频| 亚洲国产精品国产精品| 国产精品99久久久久久久久| 午夜福利在线在线| 欧美国产精品一级二级三级 | 伊人久久精品亚洲午夜| 大话2 男鬼变身卡| 亚洲综合色惰| av国产精品久久久久影院| 国产片特级美女逼逼视频| 午夜精品国产一区二区电影 | 久久久午夜欧美精品| 视频区图区小说| 亚洲内射少妇av| 亚洲高清免费不卡视频| 亚洲久久久久久中文字幕| 久久久精品免费免费高清| 在线a可以看的网站| 国产真实伦视频高清在线观看| 有码 亚洲区| 搡老乐熟女国产| 青青草视频在线视频观看| 麻豆成人午夜福利视频| 免费看av在线观看网站| 久久精品国产a三级三级三级| 白带黄色成豆腐渣| 国产真实伦视频高清在线观看| 久久久精品94久久精品| 亚洲欧美一区二区三区国产| 亚洲精品一区蜜桃| 国产精品99久久久久久久久| 久久99精品国语久久久| 亚洲国产欧美人成| 国产老妇伦熟女老妇高清| 午夜精品国产一区二区电影 | 夜夜爽夜夜爽视频| 成人亚洲精品一区在线观看 | 精品熟女少妇av免费看| 在线观看一区二区三区| 狂野欧美激情性xxxx在线观看| 国产视频首页在线观看| 精品人妻熟女av久视频| 国产欧美日韩精品一区二区| 欧美日韩一区二区视频在线观看视频在线 | 国产伦精品一区二区三区四那| 国产熟女欧美一区二区| www.色视频.com| 免费大片18禁| 国产精品国产三级国产专区5o| 老师上课跳d突然被开到最大视频| 99re6热这里在线精品视频| 亚洲精品aⅴ在线观看| 亚洲av成人精品一二三区| 欧美日韩视频高清一区二区三区二| 丝袜脚勾引网站| 直男gayav资源| 免费av观看视频| 夜夜看夜夜爽夜夜摸| 亚洲精品成人久久久久久| 精品午夜福利在线看| 黄色怎么调成土黄色| 狂野欧美白嫩少妇大欣赏| 亚洲精品一区蜜桃| 亚洲欧美日韩东京热| 汤姆久久久久久久影院中文字幕| 中文字幕av成人在线电影| 久久久久久伊人网av| 久久99热这里只频精品6学生| 国产黄色视频一区二区在线观看| 男女无遮挡免费网站观看| 在线a可以看的网站| 搡女人真爽免费视频火全软件| 99九九线精品视频在线观看视频| 超碰97精品在线观看| 日韩欧美一区视频在线观看 | 美女视频免费永久观看网站| 中国美白少妇内射xxxbb| 又黄又爽又刺激的免费视频.| 三级经典国产精品| 国产成人免费观看mmmm| 亚洲婷婷狠狠爱综合网| 新久久久久国产一级毛片| 国产精品一区二区三区四区免费观看| 免费电影在线观看免费观看| 色视频在线一区二区三区| 最近的中文字幕免费完整| 欧美97在线视频| 肉色欧美久久久久久久蜜桃 | 日本一二三区视频观看| 可以在线观看毛片的网站| 永久免费av网站大全| 18禁裸乳无遮挡动漫免费视频 | 插逼视频在线观看| 七月丁香在线播放| 国产黄片视频在线免费观看| 黄色视频在线播放观看不卡| 少妇的逼水好多| 国内少妇人妻偷人精品xxx网站| 日韩av免费高清视频| 深爱激情五月婷婷| 亚洲人成网站在线播| 国产大屁股一区二区在线视频| .国产精品久久| 国产精品秋霞免费鲁丝片| 91在线精品国自产拍蜜月| 成人鲁丝片一二三区免费| 亚洲国产精品999| 亚洲一区二区三区欧美精品 | 久久精品国产亚洲av天美| 内射极品少妇av片p| 国产精品福利在线免费观看| 国产精品国产三级国产av玫瑰| 亚洲四区av| 涩涩av久久男人的天堂| 亚洲av中文av极速乱| 在线观看国产h片| 久久久久国产精品人妻一区二区| 人人妻人人看人人澡| 久久99热这里只有精品18| 六月丁香七月| 大又大粗又爽又黄少妇毛片口| 国产午夜福利久久久久久| 18禁在线播放成人免费| 欧美 日韩 精品 国产| 国产毛片a区久久久久| 成人黄色视频免费在线看| 国产老妇伦熟女老妇高清| 国产色婷婷99| 国内少妇人妻偷人精品xxx网站| 国产成人福利小说| 熟女人妻精品中文字幕| 亚洲人成网站在线观看播放| 国产真实伦视频高清在线观看| 视频中文字幕在线观看| 狠狠精品人妻久久久久久综合| 国产欧美亚洲国产| 国产精品秋霞免费鲁丝片| 一二三四中文在线观看免费高清| 亚洲天堂国产精品一区在线| 寂寞人妻少妇视频99o| 免费不卡的大黄色大毛片视频在线观看| 一区二区三区乱码不卡18| 亚洲精品一区蜜桃| 国产女主播在线喷水免费视频网站| 大码成人一级视频| 亚洲成人久久爱视频| 秋霞伦理黄片| 免费高清在线观看视频在线观看| 久久人人爽av亚洲精品天堂 | 午夜福利视频精品| 日韩视频在线欧美| 亚洲真实伦在线观看| 大片电影免费在线观看免费| 狂野欧美激情性bbbbbb| 九草在线视频观看| 亚洲精品色激情综合| 99久久精品国产国产毛片| 国产探花在线观看一区二区| 亚洲精品aⅴ在线观看| 激情 狠狠 欧美| 麻豆精品久久久久久蜜桃| 在线播放无遮挡| 91aial.com中文字幕在线观看| 超碰av人人做人人爽久久| 中文字幕免费在线视频6| 欧美最新免费一区二区三区| 又粗又硬又长又爽又黄的视频| 99久国产av精品国产电影| 免费少妇av软件| 亚洲精品日韩在线中文字幕| 18禁在线播放成人免费| 又爽又黄无遮挡网站| av在线老鸭窝| 国内少妇人妻偷人精品xxx网站| a级毛片免费高清观看在线播放| 免费在线观看成人毛片| 夫妻性生交免费视频一级片| 亚洲av免费在线观看| 大香蕉97超碰在线| 欧美97在线视频| 日韩人妻高清精品专区| 免费观看性生交大片5| 午夜激情福利司机影院| 亚洲av在线观看美女高潮| 久久久久久久久久人人人人人人| 激情五月婷婷亚洲| av线在线观看网站| 肉色欧美久久久久久久蜜桃 | 国产欧美日韩一区二区三区在线 | 亚洲av成人精品一区久久| 亚洲av日韩在线播放| 少妇的逼好多水| 秋霞在线观看毛片| 欧美人与善性xxx| 国产成人精品婷婷| 国产大屁股一区二区在线视频| 久久热精品热| 热99国产精品久久久久久7| 午夜视频国产福利| 亚洲精品成人久久久久久| videossex国产| 亚洲人成网站高清观看| freevideosex欧美| 亚洲国产精品成人综合色| 亚洲精品一区蜜桃| 91午夜精品亚洲一区二区三区| 搡老乐熟女国产| 亚洲av不卡在线观看| .国产精品久久| 中文精品一卡2卡3卡4更新| 乱系列少妇在线播放| 国产伦在线观看视频一区| 亚洲av成人精品一二三区| 2021少妇久久久久久久久久久| 精品一区二区三卡| 亚洲真实伦在线观看| 国产免费视频播放在线视频| 精品久久久噜噜| 男女啪啪激烈高潮av片| 国产黄片视频在线免费观看| 亚洲天堂国产精品一区在线| 日韩人妻高清精品专区| 国产av国产精品国产| 国产 一区 欧美 日韩| 欧美精品一区二区大全| 黄色怎么调成土黄色| 国产乱人偷精品视频| 国产成人精品久久久久久| 国产亚洲午夜精品一区二区久久 | 男女啪啪激烈高潮av片| 国产成人a∨麻豆精品| tube8黄色片| 久久久久久九九精品二区国产| 听说在线观看完整版免费高清| 女的被弄到高潮叫床怎么办| 午夜福利在线观看免费完整高清在| 日日摸夜夜添夜夜添av毛片| 国产精品一区二区在线观看99| 亚洲国产欧美在线一区| 国产视频内射| 国产白丝娇喘喷水9色精品| xxx大片免费视频| 久久99热6这里只有精品| 免费av不卡在线播放| 黑人高潮一二区| 亚洲成人久久爱视频| 成人高潮视频无遮挡免费网站| 最近最新中文字幕免费大全7| 久久人人爽人人片av| 成人免费观看视频高清| 97超视频在线观看视频| 精品人妻一区二区三区麻豆| 午夜日本视频在线| 国产成人免费无遮挡视频| 可以在线观看毛片的网站| 日韩伦理黄色片| av在线天堂中文字幕| 亚洲精品成人av观看孕妇| 国产欧美亚洲国产| 亚洲av.av天堂| 99热这里只有是精品在线观看| 黄色配什么色好看| 一本色道久久久久久精品综合| 美女高潮的动态| 波野结衣二区三区在线| 嫩草影院精品99| 亚洲国产成人一精品久久久| 精品国产露脸久久av麻豆| 成人漫画全彩无遮挡| 在线免费观看不下载黄p国产| 亚洲四区av| 欧美精品国产亚洲| 欧美日韩视频高清一区二区三区二| 熟女人妻精品中文字幕| 亚洲久久久久久中文字幕| 亚洲人成网站在线播| a级毛色黄片| 成人一区二区视频在线观看| 美女脱内裤让男人舔精品视频| 国产欧美另类精品又又久久亚洲欧美| 最后的刺客免费高清国语| 涩涩av久久男人的天堂| 真实男女啪啪啪动态图| 欧美人与善性xxx| 久久99蜜桃精品久久| 国产综合精华液| 国产精品久久久久久av不卡| 久久久久久久大尺度免费视频| 国产探花极品一区二区| 在线观看一区二区三区| 在线 av 中文字幕| 亚洲人成网站高清观看| 欧美日韩综合久久久久久| 欧美另类一区| 亚洲成人一二三区av| 波野结衣二区三区在线| 18禁裸乳无遮挡动漫免费视频 | 欧美日韩国产mv在线观看视频 | 在线观看三级黄色| 美女cb高潮喷水在线观看| 中文字幕亚洲精品专区| 成人亚洲精品一区在线观看 | 国产亚洲av嫩草精品影院| 黄色视频在线播放观看不卡| 免费少妇av软件| 国内精品美女久久久久久| 一级a做视频免费观看| 精品一区二区三区视频在线| 色视频www国产| 波多野结衣巨乳人妻| 色综合色国产| 欧美日韩视频精品一区| 国产成人freesex在线| 亚洲精品视频女| 一级二级三级毛片免费看| 色婷婷久久久亚洲欧美| 人妻少妇偷人精品九色| 80岁老熟妇乱子伦牲交| 草草在线视频免费看| 简卡轻食公司| 亚洲成色77777| 综合色av麻豆| 国产高清三级在线| 日本-黄色视频高清免费观看| 亚洲av电影在线观看一区二区三区 | 亚洲精品国产色婷婷电影| 波多野结衣巨乳人妻| 久久人人爽av亚洲精品天堂 | 欧美日韩视频精品一区| 国产在线一区二区三区精| 免费不卡的大黄色大毛片视频在线观看| 男女那种视频在线观看| 日韩三级伦理在线观看| 婷婷色av中文字幕| 日韩人妻高清精品专区| 人人妻人人爽人人添夜夜欢视频 | 男人狂女人下面高潮的视频| 日本-黄色视频高清免费观看| 在线观看免费高清a一片| 精品久久久噜噜| 一级毛片黄色毛片免费观看视频| 毛片一级片免费看久久久久| 亚洲最大成人中文| 又粗又硬又长又爽又黄的视频| 涩涩av久久男人的天堂| 亚洲美女视频黄频| 国产片特级美女逼逼视频| 久久久午夜欧美精品| 人妻少妇偷人精品九色| 国产熟女欧美一区二区| 丝袜美腿在线中文| 午夜精品国产一区二区电影 | 国产精品久久久久久精品电影| 大片免费播放器 马上看| 91狼人影院| 少妇的逼好多水| 新久久久久国产一级毛片| 久久国产乱子免费精品| 在线免费观看不下载黄p国产| av在线蜜桃| 男男h啪啪无遮挡| 亚洲国产精品999| 国产av不卡久久| 天天一区二区日本电影三级| 色哟哟·www| av线在线观看网站| 亚洲精品日韩在线中文字幕| 国产黄色视频一区二区在线观看| 亚洲精品aⅴ在线观看| 欧美成人a在线观看| 三级男女做爰猛烈吃奶摸视频| 欧美激情在线99| 中国美白少妇内射xxxbb| 久久久精品94久久精品| 日韩一区二区三区影片| 国内少妇人妻偷人精品xxx网站| 成年人午夜在线观看视频| 国产成人精品久久久久久| 18禁裸乳无遮挡免费网站照片| av福利片在线观看| 中文精品一卡2卡3卡4更新| 精品人妻一区二区三区麻豆| 久久精品久久久久久噜噜老黄| 韩国高清视频一区二区三区| 免费高清在线观看视频在线观看| 国产成人aa在线观看| 欧美xxⅹ黑人| 男插女下体视频免费在线播放| 人妻制服诱惑在线中文字幕| 美女xxoo啪啪120秒动态图| 国产成人午夜福利电影在线观看| 一级毛片电影观看| 人妻 亚洲 视频| 少妇的逼好多水| 午夜精品一区二区三区免费看| 久久久久久久久大av| 99热这里只有精品一区| 免费观看a级毛片全部| 青春草国产在线视频| 亚洲精品aⅴ在线观看| 国产精品熟女久久久久浪| 欧美区成人在线视频| 91精品伊人久久大香线蕉| 国产日韩欧美亚洲二区| 日本wwww免费看| 中文字幕制服av| 2021天堂中文幕一二区在线观| 亚洲人与动物交配视频| 麻豆久久精品国产亚洲av| 欧美日韩国产mv在线观看视频 | 精品一区二区三区视频在线| 欧美人与善性xxx| 亚洲伊人久久精品综合| 国产乱人偷精品视频| 精品国产乱码久久久久久小说| 三级国产精品片| 午夜福利高清视频| 亚洲国产色片| 夜夜看夜夜爽夜夜摸| 欧美日韩国产mv在线观看视频 | 国产欧美日韩一区二区三区在线 | 2021天堂中文幕一二区在线观| 国产 精品1| 国产免费一级a男人的天堂| 一级毛片电影观看| 人妻 亚洲 视频| 麻豆国产97在线/欧美| 国产一级毛片在线| 边亲边吃奶的免费视频| 男人和女人高潮做爰伦理| 黄色配什么色好看| 色哟哟·www| 国产探花在线观看一区二区| 久久久久精品久久久久真实原创| 麻豆成人午夜福利视频| 永久免费av网站大全| 老司机影院成人| 一级毛片黄色毛片免费观看视频| av国产久精品久网站免费入址| 一个人看的www免费观看视频| 亚洲aⅴ乱码一区二区在线播放| 特级一级黄色大片| av在线天堂中文字幕| 免费人成在线观看视频色| 午夜激情久久久久久久| 国产色婷婷99| 777米奇影视久久| 天堂网av新在线| 18禁裸乳无遮挡动漫免费视频 | 黄色日韩在线| 在线 av 中文字幕| 日韩av不卡免费在线播放| 在线免费观看不下载黄p国产| 一个人看的www免费观看视频| 精品人妻熟女av久视频| 亚洲欧美日韩无卡精品| 天堂俺去俺来也www色官网| 草草在线视频免费看| 久久久a久久爽久久v久久| 五月玫瑰六月丁香| 久热这里只有精品99| 91在线精品国自产拍蜜月| 国产成人精品福利久久| 国产成人免费无遮挡视频| 男女啪啪激烈高潮av片| a级毛色黄片| 亚洲真实伦在线观看| videos熟女内射| 五月玫瑰六月丁香| 日韩一区二区视频免费看| 国产亚洲午夜精品一区二区久久 | av国产免费在线观看| 免费看av在线观看网站| 黄片wwwwww| 久久精品国产a三级三级三级| 久久精品久久久久久噜噜老黄| 国产一级毛片在线| 精品一区二区免费观看| 日本一本二区三区精品| 亚洲自偷自拍三级| 国产成人精品一,二区| 有码 亚洲区| 在线播放无遮挡| 我要看日韩黄色一级片| 亚洲,欧美,日韩| 日日摸夜夜添夜夜爱| 亚洲av二区三区四区| 日本一二三区视频观看| 亚洲国产高清在线一区二区三| 久久久久九九精品影院| 成年版毛片免费区| 久久97久久精品| 国产亚洲一区二区精品| 欧美xxⅹ黑人| 日本一二三区视频观看| 亚洲欧美精品专区久久| 国产免费又黄又爽又色| 日日摸夜夜添夜夜爱| 神马国产精品三级电影在线观看| videos熟女内射| 91在线精品国自产拍蜜月| 91午夜精品亚洲一区二区三区| 国产 一区 欧美 日韩| 国产av国产精品国产| 男人舔奶头视频| 国产综合懂色| 国产黄片视频在线免费观看| 91午夜精品亚洲一区二区三区| 最近最新中文字幕大全电影3| 亚洲欧美日韩卡通动漫| 日韩一本色道免费dvd| 极品少妇高潮喷水抽搐| av国产精品久久久久影院| 亚洲精品日韩在线中文字幕| 成人美女网站在线观看视频| 女人十人毛片免费观看3o分钟| 日韩av在线免费看完整版不卡| av黄色大香蕉|