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

    沉樁擠土效應(yīng)對(duì)樁周土體徑向位移數(shù)值分析

    2018-09-20 10:51:44李家華單恒年
    水道港口 2018年4期
    關(guān)鍵詞:沉樁摩擦角徑向

    李家華,單恒年,陳 祥

    (中交第四航務(wù)工程勘察設(shè)計(jì)院有限公司,廣州 510230)

    早在20世紀(jì)70年代,國(guó)內(nèi)外學(xué)者就對(duì)沉樁擠土效應(yīng)的研究做了大量的工作。歸納起來(lái),可總結(jié)為以下幾種方法:圓孔擴(kuò)張理論、應(yīng)變路徑法、有限元分析法。目前對(duì)沉樁擠土效應(yīng)的分析和解答大部分都是基于經(jīng)典平面圓孔擴(kuò)張理論[3],在巖土工程領(lǐng)域中被廣泛地運(yùn)用于沉樁引起的應(yīng)力、孔壓和土體位移、原位測(cè)試(如旁壓試驗(yàn)和靜力觸探)、爆破形成的彈坑等工程問(wèn)題的分析中,另外常用的方法還有:應(yīng)變路徑法、試驗(yàn)方法、有限元數(shù)值模擬法等。沉樁過(guò)程樁土相互作用由于其涉及到幾何非線(xiàn)性、接觸非線(xiàn)性、材料非線(xiàn)性、大變形、本構(gòu)模型等一系列的問(wèn)題,僅僅憑借著理論公式的推導(dǎo)和試驗(yàn)研究不足以滿(mǎn)足越來(lái)越精確的工程要求,以及實(shí)際施工過(guò)程中的復(fù)雜性。有限單元法(Finite Element Method)的出現(xiàn),使有限元分析(Finite Element Analysis)成為了解決這些非線(xiàn)性復(fù)雜問(wèn)題的一個(gè)有效途徑[4]。采用ABAQUS有限元分析軟件,首先建立合理的有限元模型,并與文獻(xiàn)中的離心試驗(yàn)數(shù)據(jù)進(jìn)行比較,證明模型的正確性和合理性,通過(guò)有限元模擬,研究沉樁過(guò)程中樁周土體的徑向位移變化規(guī)律。

    1 沉樁有限元分析模型構(gòu)建

    有限元模型土體的尺寸及參數(shù)參照離心試驗(yàn)來(lái)設(shè)置。假設(shè)樁徑為0.3 m,樁長(zhǎng)為5.0 m,樁靴角度為60°,樁靴采用光滑過(guò)渡。為了減小土體邊界條件的影響,將土體的分析區(qū)域長(zhǎng)度方向設(shè)為兩倍樁長(zhǎng)即10 m,寬度方向設(shè)為一倍樁長(zhǎng)即5 m,同時(shí)為了防止土體擠壓變形穿越軸線(xiàn),要預(yù)留剛性小管間隙0.001 m。土體采用修正劍橋模型,模型近似尺寸設(shè)定為40,其他土體重要參數(shù)如表1所示。

    表1 離心模型試驗(yàn)土體參數(shù)Tab.1 Soil parameters of centrifuge test data

    根據(jù)已有條件建立模型,1/2軸對(duì)稱(chēng)有限元模型及網(wǎng)格劃分情況,模型采用Abaqus/Standard中的大變形計(jì)算,將沉樁簡(jiǎn)化為平面軸對(duì)稱(chēng)問(wèn)題,通過(guò)控制分析步對(duì)樁頂施加向下位移進(jìn)行模擬沉樁,最后從后處理中導(dǎo)出模擬結(jié)果,如圖1所示。

    圖1 樁土模型(軸對(duì)稱(chēng)1/2模型)Fig.1 Pile-soil model (axisymmetric 1/2 model)

    1.1 土體本構(gòu)模型

    沉樁會(huì)使土體產(chǎn)生彈塑性區(qū),即要考慮彈性區(qū)的本構(gòu)模型也要考慮塑性區(qū)的本構(gòu)模型。彈性部分,采用多孔介質(zhì)彈性模型,它是一種非線(xiàn)性的各向同性彈性模型,可設(shè)置滲透系數(shù),與孔壓?jiǎn)卧黄鹗褂?。塑性部分,采用臨界狀態(tài)塑性模型(修正劍橋模型)。它能與彈性部分多孔介質(zhì)彈性模型很好的聯(lián)合起來(lái)使用[5],可以定義初始應(yīng)力,求解孔隙水壓力的變化。模型采用相關(guān)聯(lián)流動(dòng)法則,塑性勢(shì)面與屈服面相同,如果初始應(yīng)力狀態(tài)點(diǎn)落在屈服面外側(cè),Abaqus會(huì)自動(dòng)調(diào)整初始屈服面的位置。

    1.2 荷載邊界條件

    在初始分析步ini(Gesostatic分析步)中設(shè)置Body Force(體力)為浮重度-9 kN/m3。對(duì)于單樁貫入的問(wèn)題,樁的貫入通過(guò)樁頂參考點(diǎn)RP施加向下的位移來(lái)模擬,樁頂位移可在初始應(yīng)力平衡后的后續(xù)分析步中施加,在初始分析步中,U1、U2及UR3均設(shè)置為0。另外,在土體中心線(xiàn)處的邊界條件設(shè)置也需要特別注意。本模型簡(jiǎn)化為平面軸對(duì)稱(chēng)問(wèn)題來(lái)求解,因此不能發(fā)生穿越中心線(xiàn)的情況發(fā)生,但是遠(yuǎn)離樁中心線(xiàn)可以發(fā)生擠土效應(yīng),簡(jiǎn)單起見(jiàn),采用一個(gè)很實(shí)用的方法來(lái)解決這個(gè)問(wèn)題:在樁的底部加上一個(gè)直徑為1 mm的剛性小管,沿著小管的法線(xiàn)方向,不能承受拉力,可模擬土體開(kāi)裂,阻止土體穿越中心線(xiàn)的位移產(chǎn)生。同時(shí),還需要定義土體底面的水平和豎向位移均為0,定義土體右側(cè)的水平位移為0,土體左側(cè)的邊界條件由接觸條件來(lái)控制,無(wú)須設(shè)置,土體的邊界條件在初始分析步ini中定義生效。

    1.3 平衡初始應(yīng)力

    定義初始應(yīng)力需要注意兩個(gè)方面:首先是平衡條件,為了得到初始位移為0的狀態(tài),應(yīng)力場(chǎng)形成的等效節(jié)點(diǎn)荷載要和外荷載相平衡;其次是屈服條件,初始應(yīng)力場(chǎng)若通過(guò)高斯點(diǎn)的應(yīng)力狀態(tài)來(lái)施加,往往會(huì)出現(xiàn)一些高斯點(diǎn)的應(yīng)力狀態(tài)落在屈服面之外,雖然在后續(xù)的計(jì)算當(dāng)中,超出屈服面的應(yīng)力可以通過(guò)應(yīng)力來(lái)調(diào)整,但畢竟是不正確的,需要消耗大量的迭代,甚至可能出現(xiàn)不收斂的情況?;谝陨蟽煞矫鏃l件,常用的平衡初始應(yīng)力方法是[6]:首先將重力荷載施加于土體,并施加相應(yīng)的實(shí)際工程荷載和邊界條件,計(jì)算在重力載荷下的應(yīng)力場(chǎng),再將此應(yīng)力場(chǎng)定義為初始應(yīng)力場(chǎng)與原重力荷載一起施加在有限元模型上,這樣就可以得到滿(mǎn)足平衡條件且不違背屈服準(zhǔn)則的初始應(yīng)力場(chǎng),保證各節(jié)點(diǎn)的初始位移接近為0。若地表是水平的,則可以極大的減少定義初始應(yīng)力的工作量。

    1.4 接觸問(wèn)題

    由于接觸面的不連續(xù)約束及力學(xué)模型本身的非線(xiàn)性,在數(shù)值模擬中常常會(huì)出現(xiàn)不收斂的情況。Abaqus采用的主-從面接觸法,一個(gè)完整的接觸模擬首先需要定義主控面(master surface)和從屬面(slave surface),然后再定義它們的接觸狀態(tài),接觸屬性等。由于本模型絕大部分涉及到接觸問(wèn)題,不允許土體穿透中心線(xiàn),所以采用點(diǎn)對(duì)面的離散法,更易收斂。

    1.5 網(wǎng)格ALE技術(shù)

    圖2 深度2.0 m處樁側(cè)土體水平位移試驗(yàn)值 與模型值比較圖Fig.2 Comparison of horizontal displacement test value and model value of pile soil at depth 2.0 m

    靜壓樁沉入的過(guò)程中,在樁尖及樁端處存在著大應(yīng)變、大變形,網(wǎng)格將會(huì)扭曲的很明顯,隨著樁的沉入,邊界條件也在不斷的改變,位移的模式實(shí)際上是跳躍的,需要重新去確定邊界問(wèn)題,比較難收斂。在有限元軟件Abaqus中,提供了一系列自適應(yīng)的技術(shù),如任意拉格朗日-歐拉法自適應(yīng)網(wǎng)格劃分技術(shù)(Arbitrary Lagrangian Eulerian adaptiove meshing)等用來(lái)處理高度扭曲的網(wǎng)格。ALE分析步可以分為兩個(gè)步驟:(1)創(chuàng)建一個(gè)新的網(wǎng)格;(2)將舊網(wǎng)格的變量信息,其中包括應(yīng)力場(chǎng)、應(yīng)變場(chǎng)與動(dòng)變量等,傳輸?shù)叫戮W(wǎng)格上。這樣做可以使網(wǎng)格與物質(zhì)相分離,即使發(fā)生了網(wǎng)格大變形,也能保證在分析計(jì)算中網(wǎng)格的高質(zhì)量。ALE技術(shù)與原網(wǎng)格有一定的相關(guān)性,所以要保證初始網(wǎng)格才能提高ALE技術(shù)的效率。由于靠近中心線(xiàn)的土體區(qū)域是高應(yīng)變區(qū),網(wǎng)格的扭曲很?chē)?yán)重,所以在分析步中設(shè)置ALE區(qū)域位于離軸線(xiàn)2 m的矩形區(qū)域,同時(shí)分別設(shè)置自適應(yīng)技術(shù)里的頻率控制與對(duì)數(shù)控制。頻率控制(Frequency)主要是對(duì)每個(gè)分析步時(shí)間步長(zhǎng)里網(wǎng)格的重劃分次數(shù)控制,對(duì)于拉格朗日問(wèn)題,默認(rèn)參數(shù)為10,對(duì)歐拉問(wèn)題,默認(rèn)值為1,本模型選用1;參數(shù)控制(Mesh Sweeps)是對(duì)整個(gè)自適應(yīng)重劃分的每個(gè)子過(guò)程的控制,需要在新網(wǎng)格的基礎(chǔ)上再進(jìn)行sweep,本模型選用1。使用ALE技術(shù)會(huì)對(duì)輸出結(jié)果帶來(lái)一定的影響,在未使用ALE網(wǎng)格劃分技術(shù)時(shí),物質(zhì)點(diǎn)與網(wǎng)格點(diǎn)在計(jì)算過(guò)程中是保持相互重合的,使用完之后,網(wǎng)格的節(jié)點(diǎn)的位移將不能表示物質(zhì)點(diǎn)的位移狀態(tài)變化,而是物質(zhì)點(diǎn)移動(dòng)與網(wǎng)格重劃分共同影響的,也就是說(shuō)節(jié)點(diǎn)處的位移將不代表物質(zhì)點(diǎn)真實(shí)的位移,所有位移的云圖變化將不再是準(zhǔn)確的,而歷史變量隨時(shí)間的變化也都是無(wú)效的。

    1.6 驗(yàn)證有限元模型

    在有限元模型中對(duì)樁頂參考點(diǎn)施加向下的位移5 m,時(shí)間步設(shè)為500 s,近似模擬1 cm/s的沉樁速率,與離心模擬試驗(yàn)沉樁速率相符,由于設(shè)備有限,離心試驗(yàn)[2]只對(duì)土體位移和孔壓進(jìn)行了研究,完成后的模型數(shù)據(jù)與離心試驗(yàn)數(shù)據(jù)比較,2.0 m深處樁側(cè)土體水平位移試驗(yàn)值與模型值分布如圖2所示。

    可以看出,模型值與試驗(yàn)值具有類(lèi)似的分布規(guī)律,且誤差在允許范圍內(nèi),證明了有限元模型的合理性,可以借助該模型進(jìn)行進(jìn)一步的研究。

    2 沉樁過(guò)程土體徑向位移變化規(guī)律有限元分析

    2.1 沉樁過(guò)程土體徑向位移變化規(guī)律

    利用有限元模型模擬沉樁過(guò)程,分析沉樁至不同深度時(shí)土體徑向位移變化規(guī)律,如圖3所示。

    3-a 沉樁過(guò)程中土體表面徑向位移變化3-b 沉樁過(guò)程徑向距離r=2 d(d為樁直徑)處深度方向徑向位移變化圖3 沉樁過(guò)程土體徑向位移模型數(shù)值變化圖Fig.3 Numerical simulation of soil radial displacement model during pile driving

    由圖3分析可知:(1)地表的土體被樁向外擠出,徑向位移Ur且隨徑向距離r的增大呈對(duì)數(shù)衰減,且靠近樁的地方,擠土效應(yīng)更加明顯,衰減越快。隨著樁體的沉入,土表面徑向位移增加不大,保持對(duì)數(shù)衰減的趨勢(shì);(2)沿深度方向,徑向位移有減小的趨勢(shì),但在樁端以上0.5~1.5 m處達(dá)到最大值。Jim-HungHwang[7]等人在1994年通過(guò)觀測(cè)發(fā)現(xiàn)這個(gè)距離為10d,原因在于本文取的是均質(zhì)土層,而實(shí)際土質(zhì)并非是均勻的,并且在沉樁的過(guò)程中,土體也在不斷的重組,位移也在發(fā)生著變化。

    2.2 土體參數(shù)對(duì)土體徑向位移影響

    表2 沉樁過(guò)程不同參數(shù)下樁周土體徑向位移值Tab.2 Radial displacement of soil around pile in different parameters of pile driving process m

    整埋數(shù)據(jù),可得到沉樁過(guò)程樁周土體的徑向位移變化規(guī)律,如圖4。

    4-a 泊松比υ對(duì)徑向位移的影響 4-b 試驗(yàn)參數(shù)λ對(duì)徑向位移的影響

    4-c 試驗(yàn)參數(shù)κ對(duì)徑向位移的影響 4-d 有效內(nèi)摩擦角φ′對(duì)徑向位移的影響圖4 沉樁過(guò)程不同土體參數(shù)下樁周土體徑向位移值變化規(guī)律Fig.4 Variation law of radial displacement of soil around the pile under different soil parameters during the process of pile driving

    分析可得:(1)在距離樁入土較近的區(qū)域,土體徑向位移受泊松比影響不大,但是隨著徑向距離的增大,泊松比開(kāi)始影響土體的徑向位移,泊松比越大,徑向位移也越大;(2)試驗(yàn)參數(shù)λ、κ是劍橋模型中各向等壓加荷曲線(xiàn)和卸荷曲線(xiàn)的斜率,當(dāng)λ增大時(shí),徑向位移以一定的幅度增大,當(dāng)κ增大時(shí),則徑向位移相對(duì)減小,遠(yuǎn)離樁向的土體徑向位移基本不受影響;(3)有效內(nèi)摩擦角φ′在修正劍橋模型中存在的形式較為復(fù)雜,所以影響規(guī)律也不是單調(diào)的增加或減少,在徑向距離較小時(shí)有效內(nèi)摩擦角的改變影響不大,在0.4~0.8 m的范圍內(nèi),有效內(nèi)摩擦角的增大,徑向位移相應(yīng)減小,而在0.8 m以外有效內(nèi)摩擦角增大,徑向位移相應(yīng)增大。

    3 結(jié)語(yǔ)

    關(guān)于樁的擠土效應(yīng)的研究一直以來(lái)都是學(xué)者們關(guān)注的焦點(diǎn),從理論推導(dǎo)到試驗(yàn)研究再到現(xiàn)在的數(shù)值模擬,仍處在不斷的發(fā)展中。關(guān)于沉樁擠土效應(yīng)機(jī)理的研究更是因其特殊性而成為樁土相互作用研究的重中之重,結(jié)合離心模擬試驗(yàn)的結(jié)果,利用Abaqus有限元分析軟件建立沉樁模型,將模型計(jì)算得到徑向位移值與試驗(yàn)檢測(cè)出來(lái)的水平位移值進(jìn)行比較,驗(yàn)證模型的合理性,同時(shí)對(duì)沉樁過(guò)程樁側(cè)土體的變形特性進(jìn)行分析,為施工設(shè)計(jì)提供借鑒,具有很大的理論意義和實(shí)際意義。本文將沉樁簡(jiǎn)化為平面軸對(duì)稱(chēng)問(wèn)題,利用Abaqus有限元軟件建立沉樁模型,通過(guò)控制分析步對(duì)樁頂施加向下位移進(jìn)行模擬,最后從后處理中導(dǎo)出模擬結(jié)果。沉樁過(guò)程中,樁周土體被擠開(kāi),徑向位移比較大的區(qū)域主要集中在樁附近,并隨著與樁中心距離增加而對(duì)數(shù)衰減,到一定范圍之后則沒(méi)有什么影響。在沿沉樁的深度方向,徑向位移并不是線(xiàn)性增大,而是在樁底以上一定距離處達(dá)到最大值,這主要是受到土體重組的影響。另外,本文也研究了土質(zhì)參數(shù)對(duì)沉樁時(shí)土體徑向位移的影響,它與有限元模型中所采用的土體模型密切相關(guān),在不同區(qū)域范圍里,具有各自變化的規(guī)律,其中,泊松比系數(shù)在徑向距離大于1 m時(shí)才開(kāi)始顯著影響徑向位移的大小;試驗(yàn)參數(shù)λ增大時(shí),土體徑向位移也隨之增大;試驗(yàn)參數(shù)κ則相反;有效內(nèi)摩擦角φ′對(duì)土體徑向位移的影響規(guī)律較為復(fù)雜,在0.4~0.8 m的范圍內(nèi),有效內(nèi)摩擦角的增大,徑向位移相應(yīng)減小,而在0.8 m以外有效內(nèi)摩擦角增大,徑向位移相應(yīng)增大。

    猜你喜歡
    沉樁摩擦角徑向
    應(yīng)用摩擦角,巧解動(dòng)力學(xué)問(wèn)題
    靜壓法沉樁對(duì)周邊環(huán)境影響及質(zhì)量控制
    房建工程混凝土預(yù)制管樁靜壓沉樁施工技術(shù)
    淺探徑向連接體的圓周運(yùn)動(dòng)
    RN上一類(lèi)Kirchhoff型方程徑向?qū)ΨQ(chēng)正解的存在性
    基于PID+前饋的3MN徑向鍛造機(jī)控制系統(tǒng)的研究
    砂土層中鋼管樁高頻振動(dòng)貫入效率模型試驗(yàn)
    一類(lèi)無(wú)窮下級(jí)整函數(shù)的Julia集的徑向分布
    借助摩擦角 快解勻速運(yùn)動(dòng)問(wèn)題
    自動(dòng)沉樁技術(shù)在光伏電站建設(shè)中的應(yīng)用研究
    日韩,欧美,国产一区二区三区| 热re99久久精品国产66热6| av电影中文网址| 久久99热这里只频精品6学生| 少妇人妻精品综合一区二区| 最近2019中文字幕mv第一页| 久久久久久免费高清国产稀缺| 久久久久久人人人人人| 男人舔女人的私密视频| 国产无遮挡羞羞视频在线观看| 亚洲情色 制服丝袜| 婷婷色麻豆天堂久久| 人人妻人人爽人人添夜夜欢视频| 夜夜骑夜夜射夜夜干| 搡老乐熟女国产| 日本av免费视频播放| 天天躁狠狠躁夜夜躁狠狠躁| 丁香六月欧美| 亚洲成av片中文字幕在线观看| 欧美成人精品欧美一级黄| 日本wwww免费看| 一级毛片黄色毛片免费观看视频| 黑人巨大精品欧美一区二区蜜桃| av片东京热男人的天堂| 亚洲色图 男人天堂 中文字幕| 国产福利在线免费观看视频| 午夜福利一区二区在线看| 老司机靠b影院| 麻豆精品久久久久久蜜桃| netflix在线观看网站| 国产亚洲欧美精品永久| 午夜日本视频在线| 电影成人av| 国产免费现黄频在线看| 国产精品二区激情视频| 18禁国产床啪视频网站| 考比视频在线观看| 看非洲黑人一级黄片| 一级片免费观看大全| 精品少妇黑人巨大在线播放| 十八禁高潮呻吟视频| 欧美精品亚洲一区二区| 亚洲精品久久午夜乱码| bbb黄色大片| 成人午夜精彩视频在线观看| 秋霞伦理黄片| 久久国产精品大桥未久av| 成人午夜精彩视频在线观看| 久久婷婷青草| 制服人妻中文乱码| 免费观看a级毛片全部| 亚洲av欧美aⅴ国产| 免费高清在线观看日韩| 一本一本久久a久久精品综合妖精| 男女边吃奶边做爰视频| 搡老乐熟女国产| 国产精品欧美亚洲77777| 日日撸夜夜添| 国产成人91sexporn| 激情五月婷婷亚洲| 王馨瑶露胸无遮挡在线观看| 看免费av毛片| 五月天丁香电影| 如何舔出高潮| 精品亚洲成国产av| 欧美激情 高清一区二区三区| 飞空精品影院首页| 丝袜脚勾引网站| 美女国产高潮福利片在线看| 在现免费观看毛片| 在线看a的网站| 欧美人与性动交α欧美软件| 香蕉国产在线看| 国产精品一二三区在线看| 少妇 在线观看| 2021少妇久久久久久久久久久| 咕卡用的链子| 国精品久久久久久国模美| 多毛熟女@视频| www.熟女人妻精品国产| 国产男人的电影天堂91| 精品第一国产精品| 两个人看的免费小视频| av又黄又爽大尺度在线免费看| 日韩大码丰满熟妇| 人妻一区二区av| 丁香六月天网| 搡老岳熟女国产| 亚洲av在线观看美女高潮| 黄色 视频免费看| 亚洲欧美一区二区三区国产| 丝袜美足系列| 91精品国产国语对白视频| 国产成人免费无遮挡视频| av卡一久久| 丝袜美足系列| 人妻人人澡人人爽人人| 校园人妻丝袜中文字幕| 啦啦啦中文免费视频观看日本| 国产日韩一区二区三区精品不卡| 日韩大片免费观看网站| 一级片'在线观看视频| 大码成人一级视频| 激情五月婷婷亚洲| 日本爱情动作片www.在线观看| 丁香六月欧美| 中文字幕精品免费在线观看视频| 夜夜骑夜夜射夜夜干| 中文字幕精品免费在线观看视频| 欧美 亚洲 国产 日韩一| 国产精品嫩草影院av在线观看| 亚洲精品av麻豆狂野| 国产成人免费观看mmmm| 国产欧美日韩综合在线一区二区| 操美女的视频在线观看| 亚洲精品国产一区二区精华液| 中文字幕最新亚洲高清| 午夜久久久在线观看| 黄色怎么调成土黄色| 王馨瑶露胸无遮挡在线观看| 一二三四在线观看免费中文在| 纯流量卡能插随身wifi吗| 在线观看人妻少妇| 国产伦理片在线播放av一区| 99re6热这里在线精品视频| 19禁男女啪啪无遮挡网站| 国产福利在线免费观看视频| 99热国产这里只有精品6| 精品国产一区二区久久| 美女扒开内裤让男人捅视频| 亚洲伊人久久精品综合| 国产 一区精品| 久久99精品国语久久久| 最近最新中文字幕大全免费视频 | 最新在线观看一区二区三区 | 99re6热这里在线精品视频| 在线观看www视频免费| 国产片内射在线| 天天躁日日躁夜夜躁夜夜| 欧美成人精品欧美一级黄| 激情五月婷婷亚洲| 欧美精品av麻豆av| 久久人人爽av亚洲精品天堂| 一本久久精品| 秋霞在线观看毛片| 亚洲成人一二三区av| 午夜影院在线不卡| 久久鲁丝午夜福利片| 国产免费福利视频在线观看| 日韩人妻精品一区2区三区| 精品一区二区免费观看| 桃花免费在线播放| 自拍欧美九色日韩亚洲蝌蚪91| 黄片播放在线免费| 卡戴珊不雅视频在线播放| 久热爱精品视频在线9| 嫩草影视91久久| 在线观看一区二区三区激情| 国产一卡二卡三卡精品 | 国产亚洲午夜精品一区二区久久| 精品午夜福利在线看| 黄色怎么调成土黄色| 丝袜美足系列| 美国免费a级毛片| 老汉色∧v一级毛片| 免费看不卡的av| 男女免费视频国产| 久久久精品国产亚洲av高清涩受| 欧美日韩亚洲国产一区二区在线观看 | 欧美日韩视频精品一区| 爱豆传媒免费全集在线观看| 久久久久久久精品精品| 亚洲精品日本国产第一区| 嫩草影视91久久| 国产成人午夜福利电影在线观看| 国产成人精品福利久久| 久久久久久久久久久久大奶| 少妇人妻 视频| av天堂久久9| a级毛片黄视频| 亚洲第一青青草原| 国产成人91sexporn| 男人爽女人下面视频在线观看| 免费看不卡的av| 中文字幕人妻熟女乱码| 国产精品熟女久久久久浪| 午夜免费观看性视频| 国产精品无大码| 王馨瑶露胸无遮挡在线观看| 国产又色又爽无遮挡免| 亚洲国产av影院在线观看| 久久久久网色| 国产精品.久久久| 国产成人免费观看mmmm| 激情视频va一区二区三区| 亚洲,欧美,日韩| 日韩精品有码人妻一区| 午夜激情av网站| 免费在线观看完整版高清| 99精品久久久久人妻精品| 午夜91福利影院| 精品国产一区二区三区四区第35| 国产熟女午夜一区二区三区| 亚洲国产最新在线播放| 午夜免费观看性视频| av福利片在线| 97精品久久久久久久久久精品| 亚洲人成网站在线观看播放| 欧美精品一区二区大全| 久久这里只有精品19| 嫩草影院入口| 久久国产精品男人的天堂亚洲| 视频在线观看一区二区三区| 9191精品国产免费久久| 99re6热这里在线精品视频| 99热网站在线观看| 亚洲欧美激情在线| 精品国产露脸久久av麻豆| 亚洲国产中文字幕在线视频| 熟妇人妻不卡中文字幕| 色精品久久人妻99蜜桃| 日日摸夜夜添夜夜爱| 午夜福利视频在线观看免费| 久久99精品国语久久久| 国产精品熟女久久久久浪| 99re6热这里在线精品视频| 一级毛片 在线播放| a 毛片基地| 久久久久精品人妻al黑| 美女扒开内裤让男人捅视频| 国产乱人偷精品视频| 男的添女的下面高潮视频| 9色porny在线观看| 日本91视频免费播放| 大香蕉久久成人网| 一本大道久久a久久精品| 18禁国产床啪视频网站| 亚洲国产欧美在线一区| 纯流量卡能插随身wifi吗| 男女下面插进去视频免费观看| 国产精品女同一区二区软件| 免费少妇av软件| 久久精品亚洲av国产电影网| 亚洲国产精品成人久久小说| 黑丝袜美女国产一区| 久久ye,这里只有精品| 亚洲av成人不卡在线观看播放网 | 人人妻人人爽人人添夜夜欢视频| 一级片'在线观看视频| 欧美 日韩 精品 国产| 男女之事视频高清在线观看 | 卡戴珊不雅视频在线播放| 成年动漫av网址| 国产欧美日韩综合在线一区二区| 一边摸一边做爽爽视频免费| 99精国产麻豆久久婷婷| 日日爽夜夜爽网站| 宅男免费午夜| 美女扒开内裤让男人捅视频| 午夜精品国产一区二区电影| 视频在线观看一区二区三区| 国产精品嫩草影院av在线观看| 国产一区二区三区综合在线观看| 国产亚洲午夜精品一区二区久久| 亚洲伊人色综图| 久久天躁狠狠躁夜夜2o2o | 久久人人97超碰香蕉20202| 制服诱惑二区| 国产av国产精品国产| 亚洲精品国产av蜜桃| 国产精品.久久久| 一级毛片 在线播放| a级毛片黄视频| 80岁老熟妇乱子伦牲交| 这个男人来自地球电影免费观看 | 免费在线观看完整版高清| 最黄视频免费看| 精品卡一卡二卡四卡免费| 97在线人人人人妻| 大话2 男鬼变身卡| 国产激情久久老熟女| 下体分泌物呈黄色| 人妻一区二区av| 午夜av观看不卡| 久久97久久精品| 秋霞在线观看毛片| 人妻一区二区av| 国产精品偷伦视频观看了| 中文字幕色久视频| 欧美 日韩 精品 国产| 国产毛片在线视频| 精品国产一区二区久久| 国产精品国产三级国产专区5o| 精品国产乱码久久久久久男人| www.精华液| 午夜精品国产一区二区电影| 亚洲美女黄色视频免费看| 高清欧美精品videossex| 高清在线视频一区二区三区| 亚洲熟女精品中文字幕| 久久久久久免费高清国产稀缺| 亚洲av日韩精品久久久久久密 | 成人漫画全彩无遮挡| 国产亚洲av片在线观看秒播厂| 亚洲美女黄色视频免费看| av女优亚洲男人天堂| 久久婷婷青草| e午夜精品久久久久久久| 性色av一级| 亚洲av日韩在线播放| 丝袜美足系列| 精品酒店卫生间| 制服人妻中文乱码| 国产精品久久久人人做人人爽| 操美女的视频在线观看| 悠悠久久av| 婷婷色麻豆天堂久久| 欧美国产精品va在线观看不卡| 免费高清在线观看日韩| 国产成人a∨麻豆精品| 蜜桃在线观看..| 色综合欧美亚洲国产小说| 熟妇人妻不卡中文字幕| 亚洲一卡2卡3卡4卡5卡精品中文| 免费看av在线观看网站| 亚洲精品日韩在线中文字幕| 男人爽女人下面视频在线观看| 精品国产乱码久久久久久男人| 日日爽夜夜爽网站| 欧美日本中文国产一区发布| 日韩成人av中文字幕在线观看| 欧美变态另类bdsm刘玥| av网站在线播放免费| 国产日韩欧美亚洲二区| 日韩制服骚丝袜av| 丝瓜视频免费看黄片| 哪个播放器可以免费观看大片| 国产精品嫩草影院av在线观看| 久久国产精品大桥未久av| 久久久久精品国产欧美久久久 | 在线观看免费日韩欧美大片| 97人妻天天添夜夜摸| 精品国产露脸久久av麻豆| 中文字幕av电影在线播放| 热99久久久久精品小说推荐| 男女高潮啪啪啪动态图| 久久久精品免费免费高清| 国产精品人妻久久久影院| 久久久久久人妻| 大香蕉久久成人网| 最近2019中文字幕mv第一页| 欧美激情高清一区二区三区 | 叶爱在线成人免费视频播放| 国产成人免费无遮挡视频| 国产一区二区三区综合在线观看| 午夜免费鲁丝| 不卡视频在线观看欧美| 久久热在线av| 成人午夜精彩视频在线观看| 国产成人啪精品午夜网站| 丝袜喷水一区| 爱豆传媒免费全集在线观看| 日韩视频在线欧美| 国产精品人妻久久久影院| 黑人巨大精品欧美一区二区蜜桃| 亚洲美女视频黄频| 999久久久国产精品视频| 亚洲精品自拍成人| 在现免费观看毛片| 国产成人精品久久二区二区91 | 国产精品.久久久| 国产99久久九九免费精品| 热99久久久久精品小说推荐| 搡老乐熟女国产| 日韩一区二区视频免费看| 国产成人免费观看mmmm| 免费黄频网站在线观看国产| 亚洲成国产人片在线观看| 伊人久久大香线蕉亚洲五| 69精品国产乱码久久久| 久久久久久免费高清国产稀缺| 亚洲人成电影观看| 婷婷成人精品国产| 国产精品久久久久久人妻精品电影 | 老熟女久久久| 免费不卡黄色视频| 久久久久久人妻| 亚洲一卡2卡3卡4卡5卡精品中文| 18禁动态无遮挡网站| 超色免费av| a 毛片基地| 天天操日日干夜夜撸| 色婷婷久久久亚洲欧美| 亚洲av欧美aⅴ国产| 欧美精品亚洲一区二区| 又大又爽又粗| 成年人午夜在线观看视频| 国产成人精品久久久久久| 女人爽到高潮嗷嗷叫在线视频| 欧美日韩亚洲综合一区二区三区_| 王馨瑶露胸无遮挡在线观看| 精品国产乱码久久久久久男人| 一级毛片我不卡| 亚洲成色77777| 午夜福利视频精品| 香蕉丝袜av| 一本大道久久a久久精品| 在线天堂最新版资源| 亚洲国产欧美网| 国产精品.久久久| 亚洲婷婷狠狠爱综合网| 色94色欧美一区二区| 丝袜美足系列| 精品国产露脸久久av麻豆| 精品久久久精品久久久| 香蕉丝袜av| 中文字幕精品免费在线观看视频| 欧美xxⅹ黑人| 成年人午夜在线观看视频| 天天操日日干夜夜撸| 午夜激情av网站| 男女免费视频国产| a 毛片基地| 久久免费观看电影| 亚洲欧美日韩另类电影网站| 永久免费av网站大全| 国产精品久久久久久精品电影小说| 天天操日日干夜夜撸| 777久久人妻少妇嫩草av网站| 纵有疾风起免费观看全集完整版| 亚洲一区中文字幕在线| 亚洲av国产av综合av卡| 婷婷色av中文字幕| 91aial.com中文字幕在线观看| 亚洲精品在线美女| 亚洲一卡2卡3卡4卡5卡精品中文| 高清在线视频一区二区三区| av网站在线播放免费| 国产午夜精品一二区理论片| 欧美少妇被猛烈插入视频| 久久天堂一区二区三区四区| 黑人欧美特级aaaaaa片| 十八禁网站网址无遮挡| 色播在线永久视频| 亚洲一码二码三码区别大吗| 久久久久精品性色| 麻豆av在线久日| 麻豆精品久久久久久蜜桃| 欧美精品一区二区免费开放| 中文字幕制服av| 欧美日韩福利视频一区二区| 男女午夜视频在线观看| 2018国产大陆天天弄谢| 亚洲国产看品久久| 日本欧美视频一区| 在线天堂最新版资源| 日本爱情动作片www.在线观看| 亚洲,欧美精品.| 久久久久久久久免费视频了| 女性生殖器流出的白浆| 99九九在线精品视频| 免费女性裸体啪啪无遮挡网站| 久久久久网色| 国产精品国产三级专区第一集| 另类精品久久| 日本午夜av视频| 亚洲av成人不卡在线观看播放网 | 日韩av不卡免费在线播放| 黄网站色视频无遮挡免费观看| 男人操女人黄网站| 久久久久久人妻| 日本色播在线视频| 自线自在国产av| 永久免费av网站大全| 亚洲欧美精品综合一区二区三区| 欧美日韩成人在线一区二区| 国产精品久久久久久久久免| 日韩伦理黄色片| 狠狠精品人妻久久久久久综合| 99久久精品国产亚洲精品| 如何舔出高潮| 一级毛片 在线播放| 精品国产一区二区三区四区第35| 男女国产视频网站| 高清视频免费观看一区二区| 我的亚洲天堂| 国产精品香港三级国产av潘金莲 | 精品一区二区三区av网在线观看 | 自拍欧美九色日韩亚洲蝌蚪91| av视频免费观看在线观看| 亚洲av电影在线观看一区二区三区| 亚洲男人天堂网一区| 午夜福利视频精品| 大话2 男鬼变身卡| 秋霞在线观看毛片| 视频区图区小说| 国产深夜福利视频在线观看| 女人精品久久久久毛片| 国产免费又黄又爽又色| 伦理电影免费视频| 看十八女毛片水多多多| 曰老女人黄片| 麻豆av在线久日| 亚洲精品在线美女| 亚洲免费av在线视频| 男人添女人高潮全过程视频| 中文字幕最新亚洲高清| 午夜福利在线免费观看网站| 欧美精品亚洲一区二区| 日韩伦理黄色片| 国产一区有黄有色的免费视频| 成人影院久久| 性高湖久久久久久久久免费观看| 国产av国产精品国产| 男女之事视频高清在线观看 | 狠狠精品人妻久久久久久综合| 菩萨蛮人人尽说江南好唐韦庄| 欧美在线黄色| 精品少妇久久久久久888优播| 丝袜美腿诱惑在线| 亚洲欧美激情在线| 操美女的视频在线观看| 国产成人午夜福利电影在线观看| 丝袜美足系列| 色精品久久人妻99蜜桃| 免费黄网站久久成人精品| 久久久国产一区二区| 考比视频在线观看| 毛片一级片免费看久久久久| 天天躁夜夜躁狠狠久久av| 亚洲国产欧美日韩在线播放| 亚洲国产欧美在线一区| 十八禁网站网址无遮挡| 免费日韩欧美在线观看| 一区在线观看完整版| 在线观看一区二区三区激情| 国产高清国产精品国产三级| 久久久久久久久久久免费av| a级毛片黄视频| 99香蕉大伊视频| 美女午夜性视频免费| 亚洲精品久久午夜乱码| 男女免费视频国产| 欧美成人精品欧美一级黄| 亚洲国产欧美在线一区| 国产在线视频一区二区| 最近中文字幕高清免费大全6| 夫妻性生交免费视频一级片| www日本在线高清视频| 一区在线观看完整版| 久久国产精品男人的天堂亚洲| 天美传媒精品一区二区| 99国产综合亚洲精品| 老司机影院成人| 69精品国产乱码久久久| 国产精品香港三级国产av潘金莲 | 国产精品蜜桃在线观看| 欧美日韩亚洲国产一区二区在线观看 | 国产乱人偷精品视频| 久久ye,这里只有精品| 色网站视频免费| 亚洲三区欧美一区| 亚洲成国产人片在线观看| 中文字幕人妻丝袜一区二区 | 最新的欧美精品一区二区| 女人高潮潮喷娇喘18禁视频| 大片电影免费在线观看免费| 九草在线视频观看| 国产老妇伦熟女老妇高清| 少妇被粗大猛烈的视频| 久久精品aⅴ一区二区三区四区| av国产精品久久久久影院| 交换朋友夫妻互换小说| 久久久久久久国产电影| 丰满乱子伦码专区| 99国产综合亚洲精品| 亚洲精品在线美女| 成人三级做爰电影| av视频免费观看在线观看| 欧美日韩国产mv在线观看视频| 欧美精品人与动牲交sv欧美| 99久久人妻综合| 在线观看国产h片| av有码第一页| 亚洲五月色婷婷综合| 一级毛片 在线播放| 欧美黑人欧美精品刺激| 国产片特级美女逼逼视频| 国产精品国产三级专区第一集| tube8黄色片| 啦啦啦中文免费视频观看日本| 成年动漫av网址| av卡一久久| 免费人妻精品一区二区三区视频| 妹子高潮喷水视频| 18禁裸乳无遮挡动漫免费视频| 不卡视频在线观看欧美| 97人妻天天添夜夜摸| 国产成人91sexporn| 男女床上黄色一级片免费看| a级毛片黄视频| 91国产中文字幕| 午夜影院在线不卡| 国产1区2区3区精品| 国产一区二区三区综合在线观看| 午夜激情av网站| 亚洲av男天堂| 深夜精品福利| 亚洲精品日本国产第一区| 亚洲一卡2卡3卡4卡5卡精品中文| 欧美变态另类bdsm刘玥| 又黄又粗又硬又大视频| 在线观看免费午夜福利视频| kizo精华| 精品国产国语对白av|