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

    波浪力下小尺度圓柱群樁的內(nèi)力與變形分析

    2022-06-07 13:43:20胡春林孫帥康
    公路工程 2022年2期
    關(guān)鍵詞:排樁軸力剪力

    胡春林, 孫帥康

    (武漢理工大學(xué) 土木工程與建筑學(xué)院, 湖北 武漢 430070)

    1 概述

    波浪力是海洋等深水工程下部結(jié)構(gòu)設(shè)計(jì)的控制荷載,群樁基礎(chǔ)是此類工程中應(yīng)用最為廣泛的基礎(chǔ)類型[1],波浪力作用下群樁基礎(chǔ)的受力與響應(yīng)問題一直受到廣泛的關(guān)注。目前國(guó)內(nèi)外在該領(lǐng)域的研究已有較多成果:劉逸敏[2]等分析了淹沒承臺(tái)時(shí)圓柱樁身所受波浪力的作用特性,并給出波浪力的計(jì)算方法。張力偉[3]等利用數(shù)值軟件建立三維波浪水槽,在考慮波浪與樁基耦合的基礎(chǔ)上計(jì)算了線性波浪作用下軟土地基中圓柱樁所受波浪力的數(shù)值解。鄶艷榮[4]等基于雷諾時(shí)均Navier-Stokes方程與標(biāo)準(zhǔn)k-ε湍流模型建立波浪水槽,研究了波浪入射角、KC數(shù)等參數(shù)對(duì)跨海大橋樁基所受波浪力的影響規(guī)律。徐博[5]等基于RANS和k-ω湍流方程,建立波浪和高樁承臺(tái)的耦合模型,研究了承臺(tái)所受各方向波浪力的變化規(guī)律、承臺(tái)周圍的流場(chǎng),以及凈空對(duì)波浪力的影響規(guī)律。張耀澤[6]基于線性波理論,建立三維波浪水槽和大柔度群樁基礎(chǔ)模型,在考慮雙向流固耦合的基礎(chǔ)上,分析了群樁基礎(chǔ)受到的波浪力,并將數(shù)值解與理論解進(jìn)行比對(duì)分析,研究了樁基自身柔度對(duì)其受力的影響規(guī)律。WANG[7]等通過模型試驗(yàn)分析了作用于樁身的波浪力,并結(jié)合數(shù)值模型研究了樁身彎矩的分布規(guī)律。黃雯[8]等基于CFD理論,建立三維數(shù)值波浪水槽,研究了樁距、樁位置的變動(dòng)對(duì)群樁中各基樁受力的影響規(guī)律。何海峰[9]等基于Morison方程計(jì)算樁身波浪力,研究了波浪力作用下幾種不同結(jié)構(gòu)形式群樁基礎(chǔ)的受力特性。施勇[10]等通過物理模型試驗(yàn),構(gòu)建三樁基礎(chǔ)風(fēng)電機(jī)組結(jié)構(gòu)的縮尺物理模型,系統(tǒng)研究了結(jié)構(gòu)在規(guī)則波浪、不同方向的入射波浪和隨機(jī)波浪作用下的動(dòng)力響應(yīng)。唐云[11]等通過物理模型試驗(yàn),研究了低樁承臺(tái)直立堤在不規(guī)則波浪作用下波浪力的變化規(guī)律。綜合現(xiàn)有文獻(xiàn),可以看出目前相關(guān)研究多側(cè)重于波浪力的求解與作用特性,而對(duì)波浪力作用下樁基的動(dòng)力響應(yīng)研究相對(duì)較少且主要圍繞整個(gè)上下部結(jié)構(gòu)的響應(yīng)展開,對(duì)樁基的內(nèi)力變化、分布規(guī)律與其變形特性的研究則更為少見。本文圍繞實(shí)際工況建立三維有限元數(shù)值模型,基于線性波理論計(jì)算由通航船只引發(fā)的波浪力,分析了樁基礎(chǔ)在波浪作用下產(chǎn)生的水平位移,在與監(jiān)測(cè)數(shù)據(jù)進(jìn)行對(duì)比,驗(yàn)證模型建立合理的基礎(chǔ)上,通過數(shù)值模擬方法,研究了波浪力作用下深水小尺度圓柱形群樁的內(nèi)力變化、分布規(guī)律與其變形特性,并根據(jù)研究結(jié)果提出相關(guān)建議,可為實(shí)際工程和后續(xù)研究提供參考。

    2 工況與監(jiān)測(cè)數(shù)據(jù)

    某橋梁樁基所在處水深約46 m,全橋由(5×30)m引橋和(75+4×130+75)m預(yù)應(yīng)力混凝土剛構(gòu)—連續(xù)梁組合體系主橋組成,結(jié)構(gòu)布置簡(jiǎn)圖如圖1所示?;A(chǔ)由方形承臺(tái)(10 m×10 m×4.5 m),及樁徑2.5 m的圓形鉆孔灌注樁組成,樁長(zhǎng)76 m,單個(gè)承臺(tái)下布置4根樁,縱橋向和橫橋向樁間距均為5.9 m,如圖2所示。樁身嵌入水底巖土層30m,水中自由段長(zhǎng)46 m,承臺(tái)、樁基均采用C30混凝土。水底巖土層主要包括6 m新近沉積粉質(zhì)黏土、9 m砂礫卵石、4 m強(qiáng)風(fēng)化絹云母石英片巖和15 m中風(fēng)化絹云母石英片巖。

    在樁基與承臺(tái)施工完成時(shí)開展監(jiān)測(cè)工作,波浪數(shù)據(jù),以及由其引起的樁基水平位移采用某大型通航船只經(jīng)過附近時(shí)采集的數(shù)據(jù):波浪的波高H=0.27 m,周期T=3.50 s,波長(zhǎng)L=19.11 m。在無行船經(jīng)過、無明顯引起波浪的風(fēng)時(shí),承臺(tái)頂部水平位移最大僅3 mm,基本可以忽略不計(jì)。行船通過施工現(xiàn)場(chǎng)附近約10 s,在此時(shí)間內(nèi)測(cè)得船行波主要引起樁基沿垂直于船行方向作小幅的周期性晃動(dòng),承臺(tái)頂部最大水平位移約為15 mm。

    圖1 橋梁布置簡(jiǎn)圖Figure 1 Layout sketch of bridge

    圖2 承臺(tái)-群樁布置簡(jiǎn)圖Figure 2 Layout sketch of cap-group pile

    3 數(shù)值模型的建立

    3.1 承臺(tái)-樁基-巖土層模型

    基于實(shí)際工況和監(jiān)測(cè)數(shù)據(jù),采用Midas/GTS有限元軟件建立數(shù)值模型。承臺(tái)采用實(shí)體單元,樁柱采用梁?jiǎn)卧?,兩者的彈性模量E為3×104MPa、泊松比υ為0.25、重度γ為24.5 kN/m3。

    近沉積粉質(zhì)黏土、砂礫卵石、強(qiáng)風(fēng)化絹云母石英片巖、中風(fēng)化絹云母石英片巖均采用實(shí)體單元,本構(gòu)選用摩爾-庫(kù)倫模型,4種巖土材料的彈性模量E、泊松比υ、重度γ、粘聚力c、內(nèi)摩擦角φ如表1所示。建模時(shí)對(duì)樁與土設(shè)置接觸面以模擬兩者間的相互作用,同時(shí)地基邊界尺寸的取值也要考慮樁土相互作用影響的范圍,一般根據(jù)樁徑D,水平向通常取(12~15)D,豎向取(5~10)D[12],超出此范圍即認(rèn)為樁的應(yīng)力和變形對(duì)地基土的影響已經(jīng)很小,可以忽略不計(jì)。本模型中樁徑為

    表1 巖土材料主要力學(xué)參數(shù)Table 1 Main mechanical parameters of geotechnical materials巖土層名稱E/ MPaυγ/(kN·m-3)c/kPaφ/(°)新近沉積粉質(zhì)黏土300.2716.51117砂礫卵石1500.2520.0326強(qiáng)風(fēng)化絹云母石英片巖2×1040.2821.82528中風(fēng)化絹云母石英片巖3×1040.2724.230036

    2.5 m,結(jié)合各巖土層的厚度,地基模型的尺寸確定為50 m×50 m×34 m(長(zhǎng)×寬×高),以充分考慮樁土相互作用影響的邊界范圍。

    考慮到巖土材料的破壞通常是基于壓縮、剪切變形而發(fā)生的塑性破壞,故三維模型選擇以六面體網(wǎng)格為主生成單元,對(duì)分析的穩(wěn)定性更為有利。模型中承臺(tái)、樁基和新近沉積粉質(zhì)黏土,采用混合網(wǎng)格劃分方式,砂礫卵石、強(qiáng)風(fēng)化絹云母石英片巖和中風(fēng)化絹云母石英片巖采用映射網(wǎng)格劃分方式,以最大限度地生成六面體單元為主的網(wǎng)格。建立的有限元模型和網(wǎng)格劃分如圖3所示。

    圖3 有限元模型和網(wǎng)格劃分Figure 3 Finite element model and meshing

    3.2 波浪力計(jì)算原理

    樁徑D與波長(zhǎng)L滿足D/L=2.5/19.11=0.13<0.2,可以將其定義為小尺度樁,而水深d與波長(zhǎng)L滿足d/L=46/19.11=2.41>0.2,則可選用線性波理論描述波浪運(yùn)動(dòng)[13],基于線性波理論的小尺度圓柱樁波浪力計(jì)算通常采用MORISON[14]提出的公式(以下簡(jiǎn)稱為用Morison方程),該公式將波浪力描述為由水質(zhì)點(diǎn)水平速度引起的拖曳力,以及由水質(zhì)點(diǎn)水平加速度引起的慣性力的疊加。用Morison方程求解波浪力的示意圖如圖4所示。

    圖4 波浪力求解示意圖Figure 4 Schematic diagram of wave force solution

    根據(jù)Morison方程,樁柱長(zhǎng)度增量dz上的波浪力合力為:

    (1)

    式中:dFH為dz長(zhǎng)度樁身所受波浪力的合力;ρ為水的密度;ux為水質(zhì)點(diǎn)的水平速度;t為時(shí)間;CD、CM分別為拖拽力系數(shù)和慣性力系數(shù),與樁的截面有關(guān),按照我國(guó)現(xiàn)行規(guī)范[15],分別取1.2、2.0。

    4 模型驗(yàn)證

    基于Morison方程,使用MATLAB軟件編制波浪力計(jì)算程序,將結(jié)果導(dǎo)出,并在Midas/GTS中進(jìn)行動(dòng)力加載,然后提交計(jì)算。由于通航船只經(jīng)過施工現(xiàn)場(chǎng)附近的時(shí)間約為10 s,故此處計(jì)算了4個(gè)波浪周期,即14 s內(nèi)樁基的動(dòng)力響應(yīng)。得到波浪力作用下樁基的最大位移云圖、承臺(tái)水平位移時(shí)程曲線分別如圖5、圖6所示。

    (a)順波浪傳播方向

    (b)逆波浪傳播方向

    圖6 承臺(tái)水平位移時(shí)程曲線Figure 6 Time history curve of horizontal displacement of bearing cap

    由數(shù)值模擬結(jié)果可以看出,船行波作用下承臺(tái)最大水平位移約17 mm,曲線呈現(xiàn)出三角函數(shù)式的變化,均與實(shí)際監(jiān)測(cè)結(jié)果較為接近。說明該模型基本能夠反映實(shí)際工況,數(shù)值模擬結(jié)果可靠,樁身的內(nèi)力與變形可基于該模型進(jìn)行計(jì)算分析。

    5 內(nèi)力及變形的有限元分析

    5.1 樁端軸力

    按照波浪傳播到達(dá)樁身的順序,分別將垂直于波浪傳播方向的兩排樁命名為前排樁和后排樁。根據(jù)圖6,選取承臺(tái)水平位移時(shí)程曲線出現(xiàn)波谷和波峰的8個(gè)時(shí)間點(diǎn):1.2、2.8、4.4、6.0、7.6、9.4、11.0、12.8 s,提取上述時(shí)間點(diǎn)水中自由段樁身的樁端(樁頂與樁底)軸力數(shù)據(jù),結(jié)果顯示同排樁軸力相等,各基樁的樁頂軸力與樁底軸力也相等。此處為使圖表簡(jiǎn)明清晰,同一排僅選取一根樁的數(shù)據(jù)繪出。將以上8個(gè)時(shí)間點(diǎn)2排樁的樁端軸力數(shù)據(jù)繪制成散點(diǎn)圖并用直線連接,如圖7所示。

    圖7 樁端軸力Figure 7 Axial force of pile end

    在同一時(shí)刻,前后2排樁的軸力呈現(xiàn)出大小基本相等,方向相反的規(guī)律,且變化速率保持一致。當(dāng)受到波浪力的作用時(shí),前排樁樁端首先出現(xiàn)負(fù)軸力,短時(shí)間內(nèi)處于受壓狀態(tài),后排樁樁端則出現(xiàn)正軸力,短時(shí)間內(nèi)處于受拉狀態(tài),隨后2排樁均出現(xiàn)受拉和受壓狀態(tài)的快速交替變換。盡管混凝土承臺(tái)具有很大的重量,但樁身在水平波浪力這一周期性荷載的作用下,并非一直處于受壓狀態(tài),而是表現(xiàn)為拉壓狀態(tài)的交替變換。隨著波浪力作用的持續(xù),承臺(tái)達(dá)到峰值位移時(shí),軸力也基本達(dá)到峰值。

    5.2 樁端剪力

    同樣提取上述8個(gè)時(shí)間點(diǎn)水中自由段樁身的樁端剪力數(shù)據(jù),結(jié)果顯示同排樁的剪力數(shù)據(jù)也基本保持一致,而不同排樁的樁頂剪力則存在顯著差異,對(duì)此,將這8個(gè)時(shí)間點(diǎn)2排樁的樁頂,及樁底剪力數(shù)據(jù)分別繪制成散點(diǎn)圖并用直線連接,見圖8。

    可以看出,后排樁樁頂剪力通常要比前排樁樁頂剪力大,且隨時(shí)間的變化速率更快,同一時(shí)刻前后排樁樁頂剪力相差最大可達(dá)25.4倍,總體峰值相差1.8倍,但絕大部分時(shí)刻剪力方向保持一致。其數(shù)值相差過大的主要原因在于承臺(tái)的重量,以及前后排樁波浪的相位差。波浪首先到達(dá)前排樁,使其產(chǎn)生動(dòng)力響應(yīng)的同時(shí),承臺(tái)和后排樁同樣因受到瞬時(shí)的慣性力而產(chǎn)生加速度,引起2排樁樁頂剪力的較大差異。而同一時(shí)刻2排樁的樁底剪力則基本

    (a)樁頂剪力

    (b)樁底剪力

    相等且變化速率一致。由于波浪傳播方向上不同排樁的樁頂剪力分配較為不均勻,因此在波浪荷載較大的工程環(huán)境中,如海洋工程等領(lǐng)域,需要對(duì)樁基與承臺(tái)連接部位的設(shè)計(jì)和施工予以足夠重視,如通過加密此處的鋼筋網(wǎng)絡(luò)、增加樁身嵌入承臺(tái)的深度等措施,以防止部分連接部位出現(xiàn)剪切破壞。

    5.3 樁身彎矩

    數(shù)值模擬結(jié)果顯示同排樁的彎矩也基本相同。此處,為使曲線更簡(jiǎn)潔明了,在前文所選取的8個(gè)時(shí)間點(diǎn)中選出4個(gè)時(shí)間點(diǎn),分別為承臺(tái)水平位移時(shí)程曲線出現(xiàn)首個(gè)波谷和波峰的1.2、2.8 s,以及出現(xiàn)負(fù)向和正向最大位移的7.6、9.4 s,提取上述時(shí)間點(diǎn)的前后排樁彎矩沿水中自由段樁身的分布,繪制曲線如圖9所示。

    盡管作用于不同排樁的波浪存在相位差,但由于承臺(tái)的連接與巖土層的固定作用,且深水高樁具有一定柔性,在同一時(shí)刻,前后排樁的彎矩基本相等,樁身彎矩?cái)?shù)值在水底土層部位最大,沿樁身向上逐漸減小,至樁頂位置趨于0。

    5.4 樁身變形

    由圖5、圖6可知,承臺(tái)的變形最大且主要表現(xiàn)為水平位移,同樣選取上述4個(gè)時(shí)間點(diǎn),提取對(duì)應(yīng)的數(shù)值模擬結(jié)果,結(jié)果顯示同排樁的水平位移也保持一致。將水中自由段樁身的位移數(shù)據(jù)繪制成曲線,如圖10所示。

    圖10 樁身水平位移Figure 10 Horizontal displacement of pile body

    可以看出,樁身水平位移曲線形似“倒置”的樁身彎矩曲線。同樣因?yàn)槌信_(tái)的連接、巖土層的固定,以及深水高樁的柔性,位移的同步并不局限于樁頂與承臺(tái)連接的位置,各基樁樁身相同部位的變形均基本保持一致,變形主要表現(xiàn)為側(cè)彎,并沿波浪傳播方向來回?cái)[動(dòng)。承臺(tái)的水平位移最大,樁身嵌入土層處的位移最小,其數(shù)值接近但不等于0,說明樁身在土層中也發(fā)生了微小的位移。此前大量關(guān)于波浪力作用下深水樁基的室內(nèi)模型試驗(yàn)與數(shù)值模擬研究中,僅考慮樁身在水中的自由段,而將樁基在水底處作固接處理,顯然有一定不妥之處,樁土接觸和相互作用對(duì)樁身內(nèi)力和位移的影響不能忽略且有待進(jìn)一步研究。

    6 結(jié)論

    本文主要利用有限元數(shù)值模擬,分析了水平線性波浪作用下小尺度圓柱形群樁在水中自由段樁身的樁端軸力與樁端剪力的變化規(guī)律、樁身彎矩的分布規(guī)律和樁身變形特征,主要結(jié)論如下:

    a.在波浪力的作用下,垂直于波浪傳播方向的同排樁樁端軸力、樁端剪力和樁身彎矩基本相等;在同一時(shí)刻,前后2排樁的樁端軸力大小相等方向相反,前排樁的樁頂剪力要小于后排樁,其同一時(shí)刻相差最大可達(dá)25.4倍,峰值相差1.8倍;前后排樁的樁底剪力與樁身彎矩基本保持相等。在水平波浪力的作用下,垂直于波浪傳播方向上的不同排樁樁頂剪力分配較為不均勻,在深水樁基的設(shè)計(jì)和施工中,應(yīng)對(duì)承臺(tái)和樁基的連接予以足夠重視,如加密連接處的鋼筋網(wǎng)絡(luò)、增大樁身嵌入承臺(tái)的深度等,以避免在海洋等波浪力較大的工程環(huán)境中出現(xiàn)部分連接部位的剪切破壞。

    b.群樁中各基樁樁身相同位置的變形基本保持一致,樁身變形主要表現(xiàn)為側(cè)彎,并沿波浪傳播方向作來回的周期性晃動(dòng),承臺(tái)的水平位移最大,樁身嵌入土層處的位移最小但不為0,說明樁身在巖土層中也產(chǎn)生了一定的位移,研究中不能忽略巖土層的存在和影響。

    猜你喜歡
    排樁軸力剪力
    軸力作用下軟土基坑力學(xué)場(chǎng)演化規(guī)律的原位試驗(yàn)研究*
    淺談支撐軸力伺服系統(tǒng)在深基坑中的應(yīng)用
    排樁+斜拋撐支護(hù)體系在深基坑中的應(yīng)用
    BIM技術(shù)在軟土深基坑排樁加注漿鋼管斜撐支護(hù)體系中的應(yīng)用
    基坑鋼支撐軸力監(jiān)測(cè)的優(yōu)化及實(shí)踐
    水道港口(2020年2期)2020-06-08 04:41:12
    懸臂箱形截面梁的負(fù)剪力滯效應(yīng)
    考慮截面配筋的箱梁剪力滯效應(yīng)分析
    鋼支撐軸力在安裝過程中的監(jiān)測(cè)分析
    鐵道勘察(2016年3期)2016-08-01 07:20:14
    巖石地區(qū)基坑排樁支護(hù)新型式設(shè)計(jì)
    加筋水泥土排樁加筋參數(shù)對(duì)重載鐵路路基加固效果的影響
    日韩欧美精品v在线| 午夜福利在线在线| 欧美色欧美亚洲另类二区| 亚洲av第一区精品v没综合| 精品一区二区三区人妻视频| 可以在线观看毛片的网站| 又爽又黄a免费视频| 国产在线精品亚洲第一网站| 制服丝袜大香蕉在线| 自拍偷自拍亚洲精品老妇| 两性午夜刺激爽爽歪歪视频在线观看| 18禁黄网站禁片午夜丰满| 色吧在线观看| 国产男靠女视频免费网站| 99久久九九国产精品国产免费| 国产大屁股一区二区在线视频| 一级av片app| 精品久久久久久成人av| 高清在线国产一区| 欧美一级a爱片免费观看看| a级一级毛片免费在线观看| 欧美高清成人免费视频www| 国产亚洲精品综合一区在线观看| 人人妻人人看人人澡| 亚洲人成伊人成综合网2020| 99国产综合亚洲精品| 久久久久久久久久黄片| 中文字幕人成人乱码亚洲影| 精华霜和精华液先用哪个| 久久国产精品影院| 91狼人影院| 欧美绝顶高潮抽搐喷水| 窝窝影院91人妻| 日本撒尿小便嘘嘘汇集6| 久久亚洲精品不卡| 好男人电影高清在线观看| 999久久久精品免费观看国产| 久久久久久九九精品二区国产| 久久国产乱子伦精品免费另类| 日韩中字成人| 久久久久久久久大av| 九色成人免费人妻av| а√天堂www在线а√下载| 亚洲五月天丁香| 免费人成视频x8x8入口观看| 国产一区二区亚洲精品在线观看| 国产精品精品国产色婷婷| 中国美女看黄片| 日韩 亚洲 欧美在线| 亚洲精品影视一区二区三区av| 日日干狠狠操夜夜爽| 欧美色视频一区免费| 婷婷丁香在线五月| 老师上课跳d突然被开到最大视频 久久午夜综合久久蜜桃 | 国产精品一及| 美女xxoo啪啪120秒动态图 | 成人高潮视频无遮挡免费网站| 欧美日韩瑟瑟在线播放| 亚洲精品一卡2卡三卡4卡5卡| 免费av不卡在线播放| 亚洲va日本ⅴa欧美va伊人久久| 天天一区二区日本电影三级| av在线蜜桃| 2021天堂中文幕一二区在线观| 少妇丰满av| 日本精品一区二区三区蜜桃| 成年女人永久免费观看视频| or卡值多少钱| 乱人视频在线观看| 欧美zozozo另类| 一边摸一边抽搐一进一小说| 97热精品久久久久久| 国产精品99久久久久久久久| 日本熟妇午夜| 亚洲成人久久性| 国产高清激情床上av| 男女床上黄色一级片免费看| 18禁在线播放成人免费| 欧美成人一区二区免费高清观看| 少妇人妻精品综合一区二区 | 久久久久九九精品影院| 久久久国产成人精品二区| 99精品在免费线老司机午夜| 在线十欧美十亚洲十日本专区| 国产精品一区二区三区四区免费观看 | 亚洲电影在线观看av| 成人一区二区视频在线观看| 中出人妻视频一区二区| 色综合欧美亚洲国产小说| bbb黄色大片| 久久久色成人| 丝袜美腿在线中文| 欧美激情在线99| 88av欧美| av中文乱码字幕在线| 淫妇啪啪啪对白视频| 亚州av有码| 少妇的逼好多水| 麻豆国产97在线/欧美| 69av精品久久久久久| 三级毛片av免费| 女人被狂操c到高潮| 国产精品永久免费网站| 国产av麻豆久久久久久久| 国内精品一区二区在线观看| 91久久精品电影网| 亚洲美女视频黄频| 国内毛片毛片毛片毛片毛片| 麻豆成人午夜福利视频| 国产大屁股一区二区在线视频| 久久久久精品国产欧美久久久| 听说在线观看完整版免费高清| 2021天堂中文幕一二区在线观| 欧美+亚洲+日韩+国产| 国产免费男女视频| 9191精品国产免费久久| netflix在线观看网站| 又黄又爽又刺激的免费视频.| 蜜桃久久精品国产亚洲av| 国内精品久久久久精免费| 赤兔流量卡办理| 美女被艹到高潮喷水动态| 国产乱人伦免费视频| 久久精品人妻少妇| 欧美黑人欧美精品刺激| 日本免费a在线| 国产精品美女特级片免费视频播放器| 99久久成人亚洲精品观看| 可以在线观看毛片的网站| 国产日本99.免费观看| 日本在线视频免费播放| 国产毛片a区久久久久| 在线观看舔阴道视频| 精品一区二区三区视频在线观看免费| 亚洲av.av天堂| 亚洲国产欧美人成| 18禁裸乳无遮挡免费网站照片| 校园春色视频在线观看| 九色国产91popny在线| 国产乱人伦免费视频| 国产免费一级a男人的天堂| 一a级毛片在线观看| 熟妇人妻久久中文字幕3abv| 老司机午夜十八禁免费视频| 国产淫片久久久久久久久 | 国产伦人伦偷精品视频| 亚洲七黄色美女视频| 日日干狠狠操夜夜爽| 欧美激情在线99| 一个人观看的视频www高清免费观看| 亚洲五月婷婷丁香| 国产亚洲精品久久久com| 精品无人区乱码1区二区| 午夜激情福利司机影院| 老司机午夜福利在线观看视频| 亚洲欧美日韩高清在线视频| 国产精品乱码一区二三区的特点| 成人av一区二区三区在线看| 亚洲18禁久久av| 久久久久久久精品吃奶| 亚洲欧美激情综合另类| 午夜福利成人在线免费观看| 色尼玛亚洲综合影院| a在线观看视频网站| 精品一区二区三区视频在线| 精品午夜福利视频在线观看一区| 老司机午夜福利在线观看视频| 欧美高清性xxxxhd video| 午夜久久久久精精品| 国产av麻豆久久久久久久| 麻豆成人av在线观看| 村上凉子中文字幕在线| 国产av不卡久久| 给我免费播放毛片高清在线观看| 久久午夜亚洲精品久久| 亚洲人与动物交配视频| 国产真实乱freesex| av女优亚洲男人天堂| 最新在线观看一区二区三区| 蜜桃亚洲精品一区二区三区| 亚洲人成网站在线播| 亚洲人成网站在线播| 日本一本二区三区精品| 特级一级黄色大片| 波多野结衣高清无吗| 一二三四社区在线视频社区8| 亚洲精华国产精华精| 美女高潮的动态| 亚洲不卡免费看| 啦啦啦观看免费观看视频高清| 日本三级黄在线观看| 校园春色视频在线观看| 午夜久久久久精精品| 丰满人妻熟妇乱又伦精品不卡| 欧美一级a爱片免费观看看| 真人做人爱边吃奶动态| 日本成人三级电影网站| 在线观看一区二区三区| 乱人视频在线观看| bbb黄色大片| 亚洲av成人精品一区久久| 99国产精品一区二区三区| 欧洲精品卡2卡3卡4卡5卡区| 直男gayav资源| 黄色丝袜av网址大全| 色尼玛亚洲综合影院| 十八禁国产超污无遮挡网站| 中文字幕久久专区| 国产熟女xx| 一本综合久久免费| 麻豆久久精品国产亚洲av| 美女xxoo啪啪120秒动态图 | 欧美日韩黄片免| 九九久久精品国产亚洲av麻豆| 国产av不卡久久| 偷拍熟女少妇极品色| 天天躁日日操中文字幕| 国产欧美日韩一区二区精品| 亚洲精品456在线播放app | 在线看三级毛片| 欧美另类亚洲清纯唯美| av在线老鸭窝| 国产精品一区二区三区四区免费观看 | 久久精品国产亚洲av涩爱 | 午夜免费男女啪啪视频观看 | 给我免费播放毛片高清在线观看| 老司机午夜福利在线观看视频| 99热这里只有是精品在线观看 | 免费在线观看影片大全网站| 日本与韩国留学比较| 国产av不卡久久| 国产伦人伦偷精品视频| 很黄的视频免费| 欧美潮喷喷水| 欧美日韩瑟瑟在线播放| 亚洲色图av天堂| 久久久久亚洲av毛片大全| 国产高清激情床上av| 最近视频中文字幕2019在线8| 亚洲七黄色美女视频| 午夜精品久久久久久毛片777| 亚洲精品456在线播放app | 日韩大尺度精品在线看网址| 免费无遮挡裸体视频| 人人妻人人澡欧美一区二区| 一进一出好大好爽视频| 中文字幕熟女人妻在线| 欧美激情在线99| 欧美在线一区亚洲| 国产视频内射| 免费人成在线观看视频色| 国产黄片美女视频| 亚洲,欧美精品.| 日本黄色片子视频| 久久香蕉精品热| 99在线人妻在线中文字幕| 欧美色视频一区免费| 免费黄网站久久成人精品 | 欧美绝顶高潮抽搐喷水| 国产精品1区2区在线观看.| 婷婷精品国产亚洲av在线| 色吧在线观看| a在线观看视频网站| 国产精品一区二区免费欧美| 老司机深夜福利视频在线观看| 好男人电影高清在线观看| 日韩欧美一区二区三区在线观看| 国产国拍精品亚洲av在线观看| 一区二区三区四区激情视频 | 噜噜噜噜噜久久久久久91| 午夜精品在线福利| 午夜视频国产福利| 他把我摸到了高潮在线观看| 欧美丝袜亚洲另类 | 色综合婷婷激情| 久久热精品热| 波多野结衣高清作品| 午夜精品在线福利| 欧美国产日韩亚洲一区| 啦啦啦韩国在线观看视频| 老熟妇乱子伦视频在线观看| 一区福利在线观看| 欧美一区二区亚洲| 男人和女人高潮做爰伦理| 一卡2卡三卡四卡精品乱码亚洲| 性色avwww在线观看| 国产精品美女特级片免费视频播放器| netflix在线观看网站| 中文字幕人成人乱码亚洲影| 国产精品女同一区二区软件 | 网址你懂的国产日韩在线| av天堂在线播放| 丰满人妻一区二区三区视频av| 色播亚洲综合网| 国产亚洲欧美98| 午夜日韩欧美国产| 色在线成人网| avwww免费| 国产高清有码在线观看视频| 亚洲精品456在线播放app | 亚洲欧美激情综合另类| 午夜精品在线福利| 狂野欧美白嫩少妇大欣赏| 日本撒尿小便嘘嘘汇集6| 嫩草影院新地址| 欧美色欧美亚洲另类二区| 日韩欧美在线乱码| 亚洲av成人不卡在线观看播放网| 久久久精品欧美日韩精品| 免费在线观看成人毛片| 国产高清三级在线| 国产精品一区二区三区四区免费观看 | 变态另类成人亚洲欧美熟女| 亚洲自拍偷在线| 90打野战视频偷拍视频| 在线观看午夜福利视频| 夜夜爽天天搞| 天堂网av新在线| 色哟哟哟哟哟哟| 免费高清视频大片| 热99re8久久精品国产| 在线a可以看的网站| 亚洲七黄色美女视频| 99久久99久久久精品蜜桃| 深夜精品福利| 级片在线观看| 少妇被粗大猛烈的视频| 两人在一起打扑克的视频| 国产精品一区二区性色av| 欧美成人一区二区免费高清观看| 亚洲成av人片在线播放无| 99在线人妻在线中文字幕| 又爽又黄a免费视频| 中文字幕人成人乱码亚洲影| 欧美3d第一页| 婷婷精品国产亚洲av在线| 午夜精品久久久久久毛片777| 熟女电影av网| 一区二区三区免费毛片| 激情在线观看视频在线高清| 国内毛片毛片毛片毛片毛片| 免费av观看视频| 一二三四社区在线视频社区8| 天堂影院成人在线观看| 久久精品夜夜夜夜夜久久蜜豆| 国产精品亚洲av一区麻豆| 日韩免费av在线播放| 99riav亚洲国产免费| 国产黄片美女视频| 亚洲aⅴ乱码一区二区在线播放| 一本综合久久免费| 我要看日韩黄色一级片| 在线播放无遮挡| 国产亚洲欧美在线一区二区| 人妻丰满熟妇av一区二区三区| 香蕉av资源在线| 欧美成人免费av一区二区三区| 色在线成人网| 夜夜夜夜夜久久久久| 久久热精品热| 亚洲人成电影免费在线| 日本黄大片高清| 亚洲欧美清纯卡通| 在线观看66精品国产| 久久亚洲真实| 精品福利观看| 久久6这里有精品| 一个人看视频在线观看www免费| 久久久久久久午夜电影| 国产伦精品一区二区三区四那| 久久人人爽人人爽人人片va | 亚洲经典国产精华液单 | 国产午夜福利久久久久久| 色播亚洲综合网| 成人高潮视频无遮挡免费网站| 一级毛片久久久久久久久女| 日日摸夜夜添夜夜添av毛片 | 欧美xxxx性猛交bbbb| 一区二区三区高清视频在线| 99在线视频只有这里精品首页| 在线天堂最新版资源| 一边摸一边抽搐一进一小说| 久久国产精品影院| 最后的刺客免费高清国语| 亚洲专区国产一区二区| 国模一区二区三区四区视频| 国语自产精品视频在线第100页| 日韩av在线大香蕉| 国内精品一区二区在线观看| 国产高清有码在线观看视频| 又黄又爽又免费观看的视频| 亚洲欧美日韩卡通动漫| 在线十欧美十亚洲十日本专区| 18禁在线播放成人免费| 99国产精品一区二区三区| 国产高清有码在线观看视频| 欧美激情久久久久久爽电影| 91在线精品国自产拍蜜月| 国产精品乱码一区二三区的特点| 国产精品久久久久久人妻精品电影| 黄色配什么色好看| 男女视频在线观看网站免费| 国语自产精品视频在线第100页| 午夜福利在线观看吧| 久久香蕉精品热| 午夜精品在线福利| 亚洲专区中文字幕在线| 久久久久性生活片| 美女大奶头视频| 色av中文字幕| 午夜a级毛片| 亚洲第一电影网av| 欧美黄色片欧美黄色片| 窝窝影院91人妻| 国产精品爽爽va在线观看网站| 在线播放无遮挡| 制服丝袜大香蕉在线| 一个人看视频在线观看www免费| 女人十人毛片免费观看3o分钟| 一个人看的www免费观看视频| 午夜老司机福利剧场| 精品国产亚洲在线| 可以在线观看的亚洲视频| 人人妻,人人澡人人爽秒播| 性欧美人与动物交配| 日韩中字成人| 国产在视频线在精品| 亚洲18禁久久av| 国产av麻豆久久久久久久| 嫩草影院精品99| 在线观看舔阴道视频| 亚洲七黄色美女视频| 波野结衣二区三区在线| www.www免费av| 身体一侧抽搐| 亚洲欧美激情综合另类| 少妇被粗大猛烈的视频| 日韩欧美国产在线观看| 变态另类丝袜制服| 午夜日韩欧美国产| 岛国在线免费视频观看| 麻豆成人av在线观看| 久久久色成人| 精品人妻1区二区| 99久国产av精品| 亚洲片人在线观看| 精品一区二区三区av网在线观看| 成年女人永久免费观看视频| 国产熟女xx| 色播亚洲综合网| 亚洲国产高清在线一区二区三| 国产真实乱freesex| 人人妻人人澡欧美一区二区| 亚洲久久久久久中文字幕| 91麻豆精品激情在线观看国产| 怎么达到女性高潮| 久久久精品欧美日韩精品| 成人高潮视频无遮挡免费网站| 18美女黄网站色大片免费观看| 精品久久久久久久久av| 成人精品一区二区免费| 午夜亚洲福利在线播放| 两性午夜刺激爽爽歪歪视频在线观看| 亚洲精品影视一区二区三区av| 国产一级毛片七仙女欲春2| 美女cb高潮喷水在线观看| 日本黄大片高清| 三级毛片av免费| 亚洲 欧美 日韩 在线 免费| 国产淫片久久久久久久久 | 偷拍熟女少妇极品色| 最近最新中文字幕大全电影3| 日韩中文字幕欧美一区二区| 真人一进一出gif抽搐免费| 亚洲av成人不卡在线观看播放网| 校园春色视频在线观看| 亚洲美女视频黄频| 亚洲不卡免费看| 长腿黑丝高跟| 久久久久免费精品人妻一区二区| 欧美zozozo另类| 久久精品国产99精品国产亚洲性色| 久久精品久久久久久噜噜老黄 | 嫁个100分男人电影在线观看| 搞女人的毛片| 一夜夜www| 亚洲成人免费电影在线观看| 免费在线观看成人毛片| 国产高潮美女av| 国产精品伦人一区二区| 精品一区二区免费观看| 成熟少妇高潮喷水视频| 亚洲av成人av| 91狼人影院| 一边摸一边抽搐一进一小说| 白带黄色成豆腐渣| 简卡轻食公司| 亚洲乱码一区二区免费版| 国内精品久久久久久久电影| 亚洲精品一区av在线观看| 女人被狂操c到高潮| 91麻豆精品激情在线观看国产| 亚洲成人久久爱视频| 99视频精品全部免费 在线| 国产成人av教育| 欧美bdsm另类| 99精品久久久久人妻精品| 一级a爱片免费观看的视频| 国产色爽女视频免费观看| 我要看日韩黄色一级片| 日日摸夜夜添夜夜添小说| 日日干狠狠操夜夜爽| 永久网站在线| 综合色av麻豆| 亚洲第一欧美日韩一区二区三区| 国产精品av视频在线免费观看| 欧美潮喷喷水| 精品免费久久久久久久清纯| 国产伦精品一区二区三区四那| 久久久成人免费电影| 欧美日韩亚洲国产一区二区在线观看| 欧美最新免费一区二区三区 | av在线老鸭窝| 亚洲av一区综合| 国产精品国产高清国产av| 久久久国产成人免费| 日韩欧美精品v在线| 日本黄色片子视频| 我的女老师完整版在线观看| 小蜜桃在线观看免费完整版高清| 国产高清视频在线观看网站| 欧美最黄视频在线播放免费| 99国产精品一区二区蜜桃av| 国产极品精品免费视频能看的| 久久九九热精品免费| 日本五十路高清| 老司机午夜福利在线观看视频| 久久国产精品影院| 99久久精品国产亚洲精品| 国产精品一及| 91狼人影院| 偷拍熟女少妇极品色| 成人精品一区二区免费| 久久亚洲真实| 亚洲av成人不卡在线观看播放网| 亚洲精品影视一区二区三区av| 日韩精品中文字幕看吧| av女优亚洲男人天堂| 色噜噜av男人的天堂激情| aaaaa片日本免费| 国产免费av片在线观看野外av| 午夜日韩欧美国产| 午夜福利在线观看吧| 成年女人看的毛片在线观看| 首页视频小说图片口味搜索| 丰满人妻一区二区三区视频av| 又黄又爽又刺激的免费视频.| 免费大片18禁| 日韩中文字幕欧美一区二区| 免费在线观看成人毛片| 丁香欧美五月| 美女被艹到高潮喷水动态| 欧美黑人巨大hd| 成人鲁丝片一二三区免费| 久久伊人香网站| 欧美xxxx黑人xx丫x性爽| 偷拍熟女少妇极品色| 深夜精品福利| 简卡轻食公司| 中文字幕免费在线视频6| 午夜福利高清视频| 中国美女看黄片| 简卡轻食公司| 色哟哟哟哟哟哟| 淫秽高清视频在线观看| 高潮久久久久久久久久久不卡| 麻豆成人av在线观看| 黄色视频,在线免费观看| 中亚洲国语对白在线视频| 十八禁人妻一区二区| 欧美最黄视频在线播放免费| 国产高清激情床上av| 国产野战对白在线观看| 欧美黑人欧美精品刺激| 亚洲av熟女| 久久国产精品人妻蜜桃| 亚洲精品亚洲一区二区| 精品久久久久久久久亚洲 | 国产精品亚洲美女久久久| 乱人视频在线观看| 欧美日韩亚洲国产一区二区在线观看| 亚洲精品在线美女| 久久久久久久午夜电影| 欧美高清成人免费视频www| 欧美黄色片欧美黄色片| 一卡2卡三卡四卡精品乱码亚洲| 长腿黑丝高跟| 日韩欧美国产在线观看| 国产精品日韩av在线免费观看| 亚洲av二区三区四区| 精品无人区乱码1区二区| 99国产极品粉嫩在线观看| 日本一本二区三区精品| 最近中文字幕高清免费大全6 | 少妇裸体淫交视频免费看高清| 亚洲精品一卡2卡三卡4卡5卡| 51国产日韩欧美| 香蕉av资源在线| 国产在线精品亚洲第一网站| 一夜夜www| 日韩欧美 国产精品| 男插女下体视频免费在线播放| 亚洲午夜理论影院| av在线老鸭窝| 一区福利在线观看|