• <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ì)重載鐵路路基加固效果的影響
    熟妇人妻不卡中文字幕| 最近中文字幕高清免费大全6| 女性生殖器流出的白浆| 一级爰片在线观看| 免费少妇av软件| 国产日韩欧美在线精品| 18+在线观看网站| 精品人妻熟女毛片av久久网站| 欧美激情国产日韩精品一区| 人人妻人人爽人人添夜夜欢视频| 免费av不卡在线播放| 免费高清在线观看日韩| 国产在线免费精品| 国产成人精品无人区| 精品国产一区二区久久| 搡老乐熟女国产| 久久午夜福利片| 在线 av 中文字幕| 欧美xxxx性猛交bbbb| 国产老妇伦熟女老妇高清| 久久精品国产亚洲av天美| 午夜影院在线不卡| 五月开心婷婷网| 日韩人妻精品一区2区三区| av福利片在线| 全区人妻精品视频| 久久精品熟女亚洲av麻豆精品| 国产在线视频一区二区| 免费看光身美女| 男女无遮挡免费网站观看| 深夜精品福利| 最近最新中文字幕大全免费视频 | 在线观看www视频免费| 自拍欧美九色日韩亚洲蝌蚪91| 夫妻午夜视频| 卡戴珊不雅视频在线播放| 久久精品夜色国产| 久久精品国产鲁丝片午夜精品| 宅男免费午夜| 日本黄大片高清| 亚洲色图 男人天堂 中文字幕 | 中文字幕av电影在线播放| 又黄又粗又硬又大视频| 熟女人妻精品中文字幕| 欧美 亚洲 国产 日韩一| 毛片一级片免费看久久久久| 午夜91福利影院| 我要看黄色一级片免费的| 亚洲精品国产av蜜桃| 国产亚洲av片在线观看秒播厂| 日韩在线高清观看一区二区三区| 久久久久精品人妻al黑| 国产欧美亚洲国产| 亚洲欧美成人综合另类久久久| 精品熟女少妇av免费看| 热99久久久久精品小说推荐| 日韩一区二区三区影片| 激情视频va一区二区三区| 黄色视频在线播放观看不卡| 久久久久久久大尺度免费视频| 18禁动态无遮挡网站| 日韩精品免费视频一区二区三区 | 黄色视频在线播放观看不卡| 女人被躁到高潮嗷嗷叫费观| a级片在线免费高清观看视频| 国产熟女欧美一区二区| 男人舔女人的私密视频| h视频一区二区三区| 亚洲精品美女久久av网站| 少妇高潮的动态图| 久久久久精品人妻al黑| 99视频精品全部免费 在线| 国产男女超爽视频在线观看| 满18在线观看网站| 久久精品国产自在天天线| 国产高清国产精品国产三级| 国产免费一区二区三区四区乱码| 最近的中文字幕免费完整| 国产国语露脸激情在线看| 一区在线观看完整版| 国产黄色免费在线视频| 国产精品成人在线| 日本午夜av视频| 国产国语露脸激情在线看| 日韩三级伦理在线观看| 成人国语在线视频| 黑人猛操日本美女一级片| 两个人看的免费小视频| 一二三四中文在线观看免费高清| 久久综合国产亚洲精品| 老司机亚洲免费影院| 日日摸夜夜添夜夜爱| 嫩草影院入口| 秋霞在线观看毛片| 婷婷成人精品国产| 国产精品蜜桃在线观看| 丝袜在线中文字幕| 亚洲av国产av综合av卡| 97超碰精品成人国产| 自线自在国产av| 国产在线一区二区三区精| 建设人人有责人人尽责人人享有的| 9热在线视频观看99| 亚洲精品第二区| 少妇熟女欧美另类| 欧美精品亚洲一区二区| 在线 av 中文字幕| 亚洲五月色婷婷综合| 国产无遮挡羞羞视频在线观看| 18禁观看日本| 免费黄频网站在线观看国产| 中国美白少妇内射xxxbb| 亚洲国产精品一区二区三区在线| 久久午夜综合久久蜜桃| 天美传媒精品一区二区| 亚洲国产成人一精品久久久| 永久网站在线| 亚洲av电影在线观看一区二区三区| 国产精品久久久久久精品古装| 国产一区二区三区综合在线观看 | 三级国产精品片| 十八禁高潮呻吟视频| 国产在线视频一区二区| 一区二区av电影网| 观看美女的网站| 亚洲精品456在线播放app| 黄色配什么色好看| av又黄又爽大尺度在线免费看| a 毛片基地| 国产不卡av网站在线观看| 国产精品偷伦视频观看了| av国产精品久久久久影院| 免费黄网站久久成人精品| 国产有黄有色有爽视频| 少妇的逼好多水| 一级黄片播放器| 久久精品国产亚洲av涩爱| 国产极品粉嫩免费观看在线| www日本在线高清视频| h视频一区二区三区| 久久99精品国语久久久| 最后的刺客免费高清国语| 国产淫语在线视频| 69精品国产乱码久久久| 精品一区二区三区视频在线| 亚洲在久久综合| 亚洲欧洲国产日韩| 国产 精品1| 国产免费又黄又爽又色| 91aial.com中文字幕在线观看| 黄色毛片三级朝国网站| 久久婷婷青草| 成人国产麻豆网| 国内精品宾馆在线| 久久久欧美国产精品| av播播在线观看一区| 91在线精品国自产拍蜜月| 极品少妇高潮喷水抽搐| 美女大奶头黄色视频| 亚洲情色 制服丝袜| 日本午夜av视频| 美国免费a级毛片| 欧美97在线视频| 汤姆久久久久久久影院中文字幕| 日韩欧美精品免费久久| 亚洲精华国产精华液的使用体验| 久久精品国产自在天天线| 一本大道久久a久久精品| 午夜视频国产福利| 精品久久久久久电影网| 51国产日韩欧美| 国产一区二区三区综合在线观看 | 少妇精品久久久久久久| 欧美精品人与动牲交sv欧美| 久久韩国三级中文字幕| 日韩视频在线欧美| 欧美日韩精品成人综合77777| 亚洲综合精品二区| 亚洲性久久影院| 亚洲美女视频黄频| 国产精品国产三级国产专区5o| 久久久久久久久久久免费av| 在线观看国产h片| 国产欧美日韩一区二区三区在线| 日韩人妻精品一区2区三区| 亚洲精品av麻豆狂野| 天天影视国产精品| av播播在线观看一区| 日韩欧美精品免费久久| 亚洲精品一二三| 亚洲国产av影院在线观看| 下体分泌物呈黄色| 国产一区二区激情短视频 | 高清欧美精品videossex| 蜜桃在线观看..| 又大又黄又爽视频免费| 少妇熟女欧美另类| 亚洲av日韩在线播放| 天堂8中文在线网| 亚洲国产精品国产精品| 久久99蜜桃精品久久| 亚洲国产最新在线播放| 免费人成在线观看视频色| 国产精品国产三级专区第一集| 精品少妇久久久久久888优播| 少妇熟女欧美另类| 免费在线观看完整版高清| 在线看a的网站| 国产精品人妻久久久影院| 91精品伊人久久大香线蕉| 中文字幕免费在线视频6| 少妇被粗大的猛进出69影院 | 大香蕉97超碰在线| 另类精品久久| 少妇猛男粗大的猛烈进出视频| 中文乱码字字幕精品一区二区三区| 美女视频免费永久观看网站| kizo精华| 国产免费福利视频在线观看| 国产麻豆69| 在线观看美女被高潮喷水网站| 日韩 亚洲 欧美在线| a级毛片黄视频| 精品亚洲乱码少妇综合久久| 亚洲,欧美,日韩| 91aial.com中文字幕在线观看| 亚洲高清免费不卡视频| 亚洲性久久影院| 少妇的丰满在线观看| 国产69精品久久久久777片| 国产有黄有色有爽视频| 黄色毛片三级朝国网站| h视频一区二区三区| 国产欧美日韩综合在线一区二区| 亚洲五月色婷婷综合| 精品国产露脸久久av麻豆| 搡老乐熟女国产| 新久久久久国产一级毛片| 一级,二级,三级黄色视频| 精品少妇内射三级| 久久av网站| 18禁裸乳无遮挡动漫免费视频| 超碰97精品在线观看| 制服丝袜香蕉在线| 亚洲婷婷狠狠爱综合网| 欧美日本中文国产一区发布| 久久久精品94久久精品| videos熟女内射| 黄片无遮挡物在线观看| 国产精品嫩草影院av在线观看| 日本黄色日本黄色录像| 久久久久人妻精品一区果冻| 久久精品熟女亚洲av麻豆精品| 国产国语露脸激情在线看| 少妇 在线观看| 99视频精品全部免费 在线| 欧美日韩视频高清一区二区三区二| 少妇精品久久久久久久| 精品久久久久久电影网| 中文字幕免费在线视频6| 亚洲精品自拍成人| 欧美人与性动交α欧美软件 | 夜夜骑夜夜射夜夜干| 久久av网站| 大片电影免费在线观看免费| 国产亚洲av片在线观看秒播厂| 日本色播在线视频| 性色av一级| 男女下面插进去视频免费观看 | 亚洲国产精品一区二区三区在线| 午夜老司机福利剧场| 久久久久久久久久成人| 国产片内射在线| 乱码一卡2卡4卡精品| 国产熟女午夜一区二区三区| 黄网站色视频无遮挡免费观看| 一级毛片黄色毛片免费观看视频| 日日爽夜夜爽网站| 亚洲性久久影院| 午夜久久久在线观看| 亚洲图色成人| 婷婷色综合大香蕉| 亚洲精品456在线播放app| 久久免费观看电影| 久久久久精品性色| 亚洲欧美一区二区三区黑人 | 18禁在线无遮挡免费观看视频| 国产亚洲av片在线观看秒播厂| 成人亚洲欧美一区二区av| 免费不卡的大黄色大毛片视频在线观看| 免费日韩欧美在线观看| 少妇精品久久久久久久| 免费av不卡在线播放| 亚洲综合色惰| 99久久精品国产国产毛片| 少妇精品久久久久久久| 高清黄色对白视频在线免费看| 亚洲内射少妇av| 青春草视频在线免费观看| 欧美日韩成人在线一区二区| 一二三四中文在线观看免费高清| 亚洲欧美色中文字幕在线| 边亲边吃奶的免费视频| 视频中文字幕在线观看| 极品少妇高潮喷水抽搐| 亚洲一码二码三码区别大吗| 国产一区二区在线观看av| 亚洲欧美清纯卡通| kizo精华| 亚洲av.av天堂| 亚洲欧美日韩另类电影网站| 看免费成人av毛片| 亚洲欧洲日产国产| 自线自在国产av| 一区二区三区精品91| 丝袜美足系列| 亚洲国产毛片av蜜桃av| 国产永久视频网站| 大片电影免费在线观看免费| 少妇熟女欧美另类| 日韩在线高清观看一区二区三区| 日韩av在线免费看完整版不卡| 国产亚洲最大av| 99热6这里只有精品| 人人妻人人爽人人添夜夜欢视频| 欧美日韩精品成人综合77777| 校园人妻丝袜中文字幕| 三上悠亚av全集在线观看| 国产一区二区激情短视频 | 国产高清不卡午夜福利| 亚洲成人手机| 日韩熟女老妇一区二区性免费视频| 国产成人av激情在线播放| 国产av一区二区精品久久| 啦啦啦在线观看免费高清www| 亚洲精品一区蜜桃| 日韩不卡一区二区三区视频在线| 久久人人97超碰香蕉20202| 国产熟女午夜一区二区三区| 国产精品蜜桃在线观看| 国产一区二区在线观看日韩| 欧美 日韩 精品 国产| 街头女战士在线观看网站| 久久人人爽人人片av| 久久久久久久久久久久大奶| 亚洲欧美一区二区三区国产| 老司机影院毛片| 多毛熟女@视频| 久久久久久久国产电影| 黑丝袜美女国产一区| 一本色道久久久久久精品综合| 又黄又粗又硬又大视频| 又黄又爽又刺激的免费视频.| 男女免费视频国产| 2018国产大陆天天弄谢| 啦啦啦啦在线视频资源| 亚洲久久久国产精品| 欧美日韩精品成人综合77777| 热re99久久国产66热| videos熟女内射| 亚洲精品,欧美精品| 日日爽夜夜爽网站| 高清在线视频一区二区三区| 永久免费av网站大全| 亚洲婷婷狠狠爱综合网| 精品国产一区二区三区久久久樱花| 嫩草影院入口| 美女脱内裤让男人舔精品视频| 好男人视频免费观看在线| 成人免费观看视频高清| 国产精品三级大全| 精品人妻一区二区三区麻豆| 制服诱惑二区| 国产精品一区www在线观看| 久久女婷五月综合色啪小说| 免费久久久久久久精品成人欧美视频 | 国产黄色免费在线视频| 国产1区2区3区精品| 天堂俺去俺来也www色官网| 国产日韩欧美亚洲二区| 高清黄色对白视频在线免费看| 人人澡人人妻人| 中文字幕亚洲精品专区| 99精国产麻豆久久婷婷| 91久久精品国产一区二区三区| 丝袜在线中文字幕| 亚洲中文av在线| 我的女老师完整版在线观看| 亚洲五月色婷婷综合| 国产精品一区二区在线观看99| 亚洲成国产人片在线观看| 一边摸一边做爽爽视频免费| 久久久久久久久久久久大奶| 人妻系列 视频| 国产精品国产三级专区第一集| 国产精品一二三区在线看| 蜜臀久久99精品久久宅男| 韩国av在线不卡| 在线观看免费日韩欧美大片| 欧美成人精品欧美一级黄| av福利片在线| 韩国精品一区二区三区 | 国产高清三级在线| 91国产中文字幕| 国产精品嫩草影院av在线观看| 黄色一级大片看看| 午夜免费鲁丝| 看十八女毛片水多多多| 亚洲精品国产av成人精品| 午夜福利在线观看免费完整高清在| 最近中文字幕高清免费大全6| 久久久久久久国产电影| 国产精品一二三区在线看| 大陆偷拍与自拍| 日本av手机在线免费观看| 18禁观看日本| 90打野战视频偷拍视频| 日本爱情动作片www.在线观看| 美女中出高潮动态图| 亚洲国产精品专区欧美| 三上悠亚av全集在线观看| 超色免费av| videossex国产| 亚洲国产日韩一区二区| 伦理电影免费视频| 久久国产精品男人的天堂亚洲 | 日本欧美国产在线视频| 国产白丝娇喘喷水9色精品| 色视频在线一区二区三区| 午夜免费观看性视频| 七月丁香在线播放| 又黄又粗又硬又大视频| 夫妻性生交免费视频一级片| 又粗又硬又长又爽又黄的视频| 在线观看人妻少妇| 26uuu在线亚洲综合色| 久久精品夜色国产| 永久免费av网站大全| 美国免费a级毛片| 国产成人精品福利久久| √禁漫天堂资源中文www| 亚洲第一av免费看| 免费观看a级毛片全部| 大香蕉97超碰在线| av在线app专区| 国产精品无大码| 久久99热这里只频精品6学生| 赤兔流量卡办理| 国产成人91sexporn| 香蕉精品网在线| 亚洲av日韩在线播放| 免费大片18禁| 男人舔女人的私密视频| 91精品国产国语对白视频| 亚洲精华国产精华液的使用体验| 国产国拍精品亚洲av在线观看| 伦精品一区二区三区| 免费人成在线观看视频色| 国产成人91sexporn| 免费黄网站久久成人精品| 国产黄色免费在线视频| 老司机影院成人| 岛国毛片在线播放| 成人影院久久| 久久久久久人妻| 亚洲精品,欧美精品| 亚洲欧美中文字幕日韩二区| 国产在线免费精品| av天堂久久9| 成年动漫av网址| xxx大片免费视频| 久久久a久久爽久久v久久| 国产麻豆69| 香蕉精品网在线| 在线观看免费日韩欧美大片| 曰老女人黄片| 天堂中文最新版在线下载| 校园人妻丝袜中文字幕| 赤兔流量卡办理| 一本—道久久a久久精品蜜桃钙片| 精品亚洲乱码少妇综合久久| 精品国产一区二区久久| 如何舔出高潮| 国产一区二区三区av在线| 99热6这里只有精品| 国产永久视频网站| 午夜福利影视在线免费观看| 久久久久久久久久人人人人人人| 国产成人精品婷婷| 久久婷婷青草| 自拍欧美九色日韩亚洲蝌蚪91| 久久久a久久爽久久v久久| 亚洲av男天堂| 亚洲图色成人| 精品国产一区二区久久| 三上悠亚av全集在线观看| 国产黄色免费在线视频| 国产 精品1| 交换朋友夫妻互换小说| av卡一久久| 久久久国产精品麻豆| 青青草视频在线视频观看| 欧美激情 高清一区二区三区| 久久午夜综合久久蜜桃| 亚洲欧美成人精品一区二区| 97在线视频观看| 国产av一区二区精品久久| 欧美另类一区| 一区二区日韩欧美中文字幕 | 亚洲国产精品一区三区| 欧美变态另类bdsm刘玥| 黑人猛操日本美女一级片| 欧美精品高潮呻吟av久久| 日日撸夜夜添| 国产乱人偷精品视频| 26uuu在线亚洲综合色| 亚洲精品中文字幕在线视频| av播播在线观看一区| 国产精品蜜桃在线观看| 又大又黄又爽视频免费| 久久久久久久国产电影| 亚洲国产精品专区欧美| 午夜影院在线不卡| 国产亚洲午夜精品一区二区久久| 日韩一区二区视频免费看| 一级毛片电影观看| 日韩电影二区| 国产一区二区在线观看av| 如何舔出高潮| av在线观看视频网站免费| 美女国产高潮福利片在线看| 欧美日韩综合久久久久久| 精品少妇久久久久久888优播| 两个人免费观看高清视频| 日韩中文字幕视频在线看片| 熟妇人妻不卡中文字幕| av在线观看视频网站免费| 久久久久网色| 老司机影院毛片| 久久 成人 亚洲| 黑丝袜美女国产一区| 高清欧美精品videossex| 制服丝袜香蕉在线| 欧美国产精品va在线观看不卡| 永久网站在线| a 毛片基地| 91午夜精品亚洲一区二区三区| 99热全是精品| 超色免费av| 啦啦啦啦在线视频资源| 欧美日韩视频精品一区| 欧美亚洲日本最大视频资源| 街头女战士在线观看网站| 我要看黄色一级片免费的| 成年人免费黄色播放视频| 国产精品一区二区在线不卡| 亚洲第一av免费看| 亚洲伊人久久精品综合| 国产精品久久久久久精品古装| 亚洲av日韩在线播放| av国产久精品久网站免费入址| 人妻一区二区av| 国产精品三级大全| 国产av一区二区精品久久| 国产视频首页在线观看| 午夜福利,免费看| 人人妻人人添人人爽欧美一区卜| 黄色一级大片看看| 在线观看免费高清a一片| 看免费成人av毛片| 国产精品女同一区二区软件| av视频免费观看在线观看| 欧美精品国产亚洲| 老女人水多毛片| 精品一区二区三区四区五区乱码 | 日韩伦理黄色片| 国产精品久久久久久久久免| 精品久久久久久电影网| 国产色婷婷99| 久久午夜综合久久蜜桃| av有码第一页| av不卡在线播放| 国产免费一级a男人的天堂| 黄网站色视频无遮挡免费观看| 国产欧美日韩综合在线一区二区| 在线观看人妻少妇| 免费av中文字幕在线| 久久精品熟女亚洲av麻豆精品| 国产片内射在线| 国产欧美亚洲国产| 婷婷色综合www| 视频在线观看一区二区三区| 纯流量卡能插随身wifi吗| 日韩大片免费观看网站| 国产又爽黄色视频| av一本久久久久| 在线天堂中文资源库| 丰满少妇做爰视频| 99国产综合亚洲精品| 男人爽女人下面视频在线观看| 亚洲av欧美aⅴ国产| 久久精品久久久久久久性| 人人妻人人澡人人看| 伊人久久国产一区二区| 亚洲综合色网址| 蜜桃在线观看..| 卡戴珊不雅视频在线播放| 国产成人精品久久久久久| 免费黄色在线免费观看| 26uuu在线亚洲综合色| 夜夜爽夜夜爽视频| 青青草视频在线视频观看| 99视频精品全部免费 在线| 最近的中文字幕免费完整| 又大又黄又爽视频免费|