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

    施工期臺(tái)風(fēng)浪作用下沉箱碰撞破壞分析

    2018-09-12 11:28:44孫熙平李越松
    關(guān)鍵詞:沉箱隔板夾角

    陳 宇,孫熙平,李越松

    (交通運(yùn)輸部天津水運(yùn)工程科學(xué)研究所水工構(gòu)造物檢測、診斷與加固技術(shù)交通行業(yè)重點(diǎn)實(shí)驗(yàn)室,天津 300456)

    我國擁有漫長的海岸線,包括東部的渤海、黃海、東海以及南部的南海。隨著近幾年國家經(jīng)濟(jì)的快速發(fā)展,沿海開發(fā)碼頭、人工島等大型工程越來越多,體量越來越大,而我國又是受臺(tái)風(fēng)襲擊較為嚴(yán)重的國家之一,尤其東海、南海沿岸受熱帶氣旋引起的臺(tái)風(fēng)災(zāi)害嚴(yán)重。如2016年10月我國南部發(fā)生的“莎莉嘉”臺(tái)風(fēng)導(dǎo)致正在施工的沉箱護(hù)岸結(jié)構(gòu)發(fā)生嚴(yán)重破壞。經(jīng)調(diào)查發(fā)現(xiàn),施工期沉箱在臺(tái)風(fēng)浪作用下發(fā)生碰撞破壞是其主要破壞形式。因此,研究臺(tái)風(fēng)浪作用下沉箱碰撞破壞模式及破壞機(jī)理對(duì)提高施工期沉箱抗浪抗撞性能具有重要指導(dǎo)意義。

    關(guān)于臺(tái)風(fēng)浪作用下沉箱的碰撞破壞分析,首先需要確定作用在沉箱結(jié)構(gòu)上的波浪力,現(xiàn)有的Morsion方程[1-2]和繞射波浪理論[3-4]只針對(duì)柱體波浪力進(jìn)行求解,無法求解直立式沉箱結(jié)構(gòu)的波浪力,而《港口與航道水文規(guī)范》[5]雖然給出了直立結(jié)構(gòu)上的波浪力分布,但只能確定沉箱上的波浪力大小,要想得到沉箱碰撞過程的運(yùn)動(dòng)趨勢,必須獲得作用在沉箱上的波浪力時(shí)程。其次,結(jié)構(gòu)間發(fā)生碰撞的時(shí)間短,應(yīng)力波會(huì)在結(jié)構(gòu)中傳播,導(dǎo)致結(jié)構(gòu)發(fā)生非線性變形[6],因此需建立較為精細(xì)的有限元模型,準(zhǔn)確模擬瞬時(shí)碰撞導(dǎo)致結(jié)構(gòu)的破壞效應(yīng)[7]。傳統(tǒng)模擬技術(shù)往往需要實(shí)現(xiàn)波浪模型和結(jié)構(gòu)模型一體化建模,并實(shí)現(xiàn)結(jié)構(gòu)間碰撞效應(yīng)的模擬,計(jì)算量大,計(jì)算成本高,或因采用大尺寸有限元網(wǎng)格而降低計(jì)算精度[8]。

    本文提出一種將臺(tái)風(fēng)浪作用下施工期沉箱碰撞反應(yīng)分析分兩步走的數(shù)值模擬方法。第一步,利用高精度流體模擬軟件Flow-3D的FAVOR網(wǎng)格技術(shù)建立流體網(wǎng)格,采用連續(xù)性方程和Navier-Stokes方程作為控制方程,通過GMRES(Generalizedminimum RESidual)隱式算法求解壓力-速度耦合方程組,獲得波浪作用于沉箱結(jié)構(gòu)的波浪力時(shí)程;第二步,建立沉箱結(jié)構(gòu)精細(xì)化模型,并將第一步獲得的波浪力時(shí)程施加于沉箱結(jié)構(gòu),利用LS-DYNA顯式求解器[9]分析波浪作用下沉箱結(jié)構(gòu)間的碰撞反應(yīng)和破壞過程。將該方法應(yīng)用于“莎莉嘉”臺(tái)風(fēng)引起施工期沉箱的碰撞破壞分析,研究了臺(tái)風(fēng)引起的波浪作用方向與沉箱結(jié)構(gòu)碰撞破壞模式之間的關(guān)系。

    1 “莎莉嘉”臺(tái)風(fēng)引起的施工期沉箱碰撞破壞

    1.1 工程概況

    如意島填海工程位于??谑忻捞m區(qū)北部外海域,利用海域面積為715.75 萬m2,填海工程護(hù)岸采用直立式沉箱結(jié)構(gòu)形式。沉箱并排布置于拋石基床上,2個(gè)沉箱平行安裝間距為20 cm,沉箱尺寸如圖1所示。設(shè)計(jì)中沉箱混凝土為C40,鋼筋為HRB400,鋼筋配筋率為0.8%,工程設(shè)計(jì)水深為7.35m。

    圖1 沉箱尺寸(單位:mm)Fig.1 Caisson sizes(unit:mm)

    1.2 事故概況

    2016年10月18日上午第21號(hào)臺(tái)風(fēng)“莎莉嘉”在海南省萬寧市和樂鎮(zhèn)沿海登陸,登陸時(shí)中心附近最大風(fēng)力14級(jí),中心最低氣壓955 hPa。由于無法獲得臺(tái)風(fēng)發(fā)生時(shí)所在水域的波浪信息,根據(jù)文獻(xiàn)[10] 的臺(tái)風(fēng)特征參數(shù)及臺(tái)風(fēng)風(fēng)速預(yù)報(bào)臺(tái)風(fēng)浪波高的回歸參數(shù),可得14級(jí)臺(tái)風(fēng)作用為主導(dǎo)因素的臺(tái)風(fēng)浪的有效波高為5.0m,平均周期為8.8 s。

    圖2 沉箱內(nèi)隔板裂縫Fig.2 Cracks of caisson interior walls

    2016年11月交通運(yùn)輸部天津水運(yùn)工程科學(xué)研究所交通水運(yùn)工程試驗(yàn)檢測中心人員對(duì)受損沉箱進(jìn)行了現(xiàn)場檢測,沉箱間碰撞引起的破壞主要表現(xiàn)為4個(gè)內(nèi)隔板Oe,Of,Og,Oh(位置編號(hào)見圖1)出現(xiàn)不同程度裂縫,而沉箱4個(gè)外壁未發(fā)現(xiàn)明顯的裂縫破壞。圖2為現(xiàn)場檢查發(fā)現(xiàn)沉箱內(nèi)隔板裂縫情況。圖2(a)為沉箱內(nèi)隔板Oh發(fā)生豎向裂縫,圖2(b)為沉箱內(nèi)隔板Og發(fā)生豎向裂縫。同時(shí)現(xiàn)場發(fā)現(xiàn)裂縫主要集中于內(nèi)隔板中部和端部節(jié)點(diǎn)處。

    2 沉箱結(jié)構(gòu)的波浪力時(shí)程

    圖3 水槽-結(jié)構(gòu)三維計(jì)算模型Fig.3 3D numericalmodel for flume-structure

    本文主要研究海南如意島工程的護(hù)岸沉箱結(jié)構(gòu)在“莎莉嘉”臺(tái)風(fēng)引起的臺(tái)風(fēng)浪作用下結(jié)構(gòu)所受波浪力作用時(shí)程。高精度流體模擬軟件Flow-3D的造波技術(shù)[11]為用戶提供了三維波浪力對(duì)沉箱結(jié)構(gòu)作用的求解。由于無法獲得事故發(fā)生時(shí)波浪作用方向與沉箱x向(沿沉箱前后壁方向)的夾角,分析時(shí)取夾角在0°~90°范圍內(nèi)每隔30°為1種工況,當(dāng)夾角為90°時(shí)波浪力作用方向與沉箱并排方向垂直,導(dǎo)致并排放置的沉箱無法發(fā)生碰撞,因此僅考慮3種工況進(jìn)行分析,分別是波浪作用方向與x向的夾角為0°,30°,60°。以波浪作用方向與x方向夾角為0°為例闡述模型的建立過程。應(yīng)用Flow-3D軟件建立三維數(shù)值水槽-結(jié)構(gòu)模型,水槽長120m,寬10m,高12m,沉箱前后壁(x向)與水槽入射方向夾角為0°,三維水槽-結(jié)構(gòu)模型如圖3所示。

    圖4 邊界條件示意Fig.4 Sketch of boundary conditions

    邊界條件及初始條件:水槽邊界條件設(shè)置如圖4所示,水槽前端設(shè)置造波(Wave)邊界,末端設(shè)置靜水壓力(Specified pressure)邊界,水槽底面設(shè)為壁面(Wall)邊界;其余各面均設(shè)為對(duì)稱(Symmetry)邊界。初始時(shí)刻整個(gè)水槽中的初始速度為0,壓強(qiáng)為沿z向的靜水壓強(qiáng),采用靜水水面初始化流場。

    圖5 造波圖Fig.5 Wave generation

    根據(jù)事故概況中推算的“莎莉嘉”臺(tái)風(fēng)浪的有效波高為5.0m,平均周期為8.8 s以及設(shè)計(jì)資料的水深為7.35m,將波高、周期和水深輸入Flow-3D造波系統(tǒng),可以得到水槽造波如圖5所示。

    圖6 波浪發(fā)生越浪現(xiàn)象Fig.6 Overtopping phenomenon

    圖7 不同波浪作用與x向夾角的波浪力時(shí)程曲線Fig.7 Wave force time-history curves with different angles between wave direction and x direction

    計(jì)算時(shí)間取5個(gè)波浪周期,計(jì)算過程中波浪作用于沉箱上發(fā)生越浪現(xiàn)象,如圖6所示。通過記錄作用于沉箱結(jié)構(gòu)面上的波壓力時(shí)程可獲得沉箱所受波浪力。同理,建立了波浪作用方向與x方向夾角分別為30°,60°的三維水槽-結(jié)構(gòu)模型,通過記錄波壓力時(shí)程曲線獲得3種工況下作用于沉箱上的波浪力,如圖7所示。從圖7可以看出,波浪作用方向與x方向夾角為0°,30°,60°的沉箱結(jié)構(gòu)上所受波浪力最大值分別為2 710,3 427和3 613 kN,最大波浪力發(fā)生在第1個(gè)波峰作用時(shí)刻,而后由于波浪發(fā)生越浪現(xiàn)象,流體對(duì)結(jié)構(gòu)波浪作用力呈現(xiàn)明顯的非規(guī)則沖擊效應(yīng);同時(shí)可以看出,隨著波浪作用方向與x向夾角的增大,波浪力逐漸增大,其原因是隨著夾角的增大,迎水面面積增大從而導(dǎo)致波浪力增大。

    3 “莎莉嘉”臺(tái)風(fēng)引起沉箱碰撞破壞分析

    3.1 沉箱精細(xì)化有限元模型

    利用LS-DYNA軟件建立2個(gè)沉箱的計(jì)算模型,沉箱底面設(shè)置剛性地基簡化模擬拋石基床基礎(chǔ),如圖8所示。采用鋼筋混凝土分離式模型來模擬鋼筋混凝土沉箱結(jié)構(gòu),沉箱精細(xì)化模型中混凝土采用常應(yīng)力積分實(shí)體單元,鋼筋采用Belytschko梁單元,鋼筋和混凝土之間采用共節(jié)點(diǎn)的方式模擬,不考慮鋼筋混凝土之間的粘結(jié)滑移效應(yīng),鋼筋混凝土沉箱有限元模型如圖9所示。為了減小計(jì)算成本,設(shè)定沉箱B為剛體,重點(diǎn)分析的沉箱A混凝土材料采用損傷本構(gòu)模型,鋼筋采用線性強(qiáng)化本構(gòu)模型,模型單元尺寸為70mm。數(shù)值分析發(fā)現(xiàn),進(jìn)一步細(xì)化網(wǎng)格只能夠有限地提高計(jì)算精度,但計(jì)算成本成倍增加,甚至?xí)霈F(xiàn)計(jì)算機(jī)的溢出。

    圖8 沉箱整體模型Fig.8 Model for caissons

    圖9 沉箱精細(xì)化有限元模型Fig.9 Refined finite elementmodel for caisson

    沉箱精細(xì)化模型中混凝土材料采用LS-DYNA中的MAT_CONCRETE_DAMAGE_REL3模型,該模型可同時(shí)考慮材料的靜水壓力作用、材料的應(yīng)變率效應(yīng)、應(yīng)力強(qiáng)化效應(yīng)和材料的殘余強(qiáng)度,曾經(jīng)被成功地用于預(yù)測混凝土在沖擊荷載作用下單向、雙向和三向拉壓狀態(tài)下的響應(yīng)[12]。同時(shí)混凝土材料在碰撞等高加載率荷載作用下,具有顯著的應(yīng)變率效應(yīng),其強(qiáng)度可顯著提高[13-14],材料的應(yīng)變率效應(yīng)一般通過動(dòng)力增強(qiáng)系數(shù)(Dynamic Increase Factor)來考慮。

    混凝土強(qiáng)度的動(dòng)力增強(qiáng)系數(shù)采用K & C模型,其抗拉強(qiáng)度的TDIF定義[15]為

    (1)

    混凝土抗拉強(qiáng)度的CDIF定義為

    (2)

    表1 材料參數(shù)Tab.1 Material parameters

    鋼筋采用LS-DYNA中的MAT_PLASTIC_KINEMATIC模型,其強(qiáng)度的FDIF定義[16]為

    (3)

    α=0.074-0.040fy/414

    (4)

    式中:fy為鋼筋屈服強(qiáng)度(MPa)?;炷梁弯摻畹牟牧蠀?shù)見表1。

    2個(gè)沉箱的碰撞過程需采用接觸算法進(jìn)行模擬,對(duì)沉箱與剛性地基間的相互作用以及沉箱間的碰撞反應(yīng)采用LS-DYNA中CONTACT_AUTOMATIC_SURFACE_TO_SURFACE,沉箱與地基之間的摩擦系數(shù)按文獻(xiàn)[17] 取0.6,沉箱間混凝土的摩擦系數(shù)按文獻(xiàn)[18] 取0.5。

    3.2 沉箱碰撞破壞分析

    通過記錄作用于沉箱的波浪力時(shí)程曲線,將所記錄的波浪力時(shí)程曲線編寫成LS-DYNA的模型輸入文件,并分別按波浪力與x向夾角為0°,30°,60°共3種工況施加于所建立的沉箱B重心上,模擬波浪入射方向與x向成不同夾角時(shí)沉箱結(jié)構(gòu)的碰撞反應(yīng)。

    采用上述數(shù)值模擬方法模擬了3種工況下沉箱的碰撞破壞模式。圖10給出了波浪作用方向與x方向夾角分別為0°,30°,60°對(duì)應(yīng)3種工況a,b,c下沉箱的碰撞破壞模式??梢姡?種不同的波浪與沉箱作用方向引起沉箱碰撞破壞模式與破壞程度不盡相同。在工況a中,沉箱B與沉箱A沿0°發(fā)生碰撞,沉箱外壁AD和BC發(fā)生了嚴(yán)重破壞,在外壁節(jié)點(diǎn)處和中部出現(xiàn)了嚴(yán)重的裂縫破損。同時(shí),內(nèi)隔板Oe和Og發(fā)生了嚴(yán)重破壞,內(nèi)隔板Of,Oh未發(fā)現(xiàn)明顯裂縫破損,內(nèi)隔板g,e,O節(jié)點(diǎn)處和板中部出現(xiàn)了嚴(yán)重的裂縫破損。在工況b中,沉箱B與沉箱A沿30°發(fā)生斜向碰撞,沉箱外壁AD和BC發(fā)生了較為嚴(yán)重的破壞,內(nèi)隔板Oe,Og發(fā)生了嚴(yán)重的裂縫破壞,內(nèi)隔板Of,Oh也發(fā)生了輕微的裂縫破損,其破壞程度較內(nèi)隔板Oe和Og輕,內(nèi)隔板g,e,O節(jié)點(diǎn)處和板中部出現(xiàn)了嚴(yán)重的裂縫破損,可以看出內(nèi)隔板中部和節(jié)點(diǎn)處的破損現(xiàn)象與圖2(b)所示破壞模式相近。在工況c中,沉箱B與沉箱A沿60°發(fā)生斜向碰撞,沉箱4個(gè)外壁的破損程度較小,4個(gè)內(nèi)隔板的破壞模式與工況b相似,但內(nèi)隔板Oe和Og破損程度較工況b輕,內(nèi)隔板Of和Oh的破損程度與工況a和b相比最嚴(yán)重。出現(xiàn)該現(xiàn)象的原因是波浪作用方向與x向的夾角就越小,沉箱B與沉箱A發(fā)生碰撞時(shí)的夾角就越小,此時(shí)內(nèi)隔板Oe和Og所受x方向慣性力越大,易引起內(nèi)隔板中部和端部發(fā)生破壞,而此時(shí)內(nèi)隔板Of和Oh所受y方向的慣性力越小,越不容易發(fā)生破損。

    綜上所述,當(dāng)波浪作用方向與x向夾角為60°時(shí),沉箱破壞模式與現(xiàn)場檢查的結(jié)果相近,即沉箱外壁未發(fā)生明顯破壞,內(nèi)隔板發(fā)生不同程度的裂縫破損。值得注意的是,現(xiàn)場無法獲得實(shí)際波浪作用與沉箱夾角數(shù)據(jù),假定波浪作用方向與x向夾角為0°,30°,60°共3種工況進(jìn)行計(jì)算,數(shù)值分析表明夾角為60°時(shí)其模擬結(jié)果與現(xiàn)場情況相近,但并不代表實(shí)際波浪作用方向與沉箱夾角為60°。

    圖10 沉箱B與A沿不同角度碰撞的破壞模式Fig.10 Pounding damagemode of caissons with different angles between caissons B and A

    4 結(jié) 語

    提出了一種基于臺(tái)風(fēng)浪作用下施工期沉箱碰撞反應(yīng)分析分兩步走的數(shù)值模擬方法,并將該方法應(yīng)用于“莎莉嘉”臺(tái)風(fēng)引起施工期沉箱碰撞破壞分析研究。通過數(shù)值分析,得出如下結(jié)論:

    (1)所提出的方法可以模擬施工期沉箱在臺(tái)風(fēng)浪作用下的碰撞反應(yīng),不同波浪作用方向引起沉箱發(fā)生不同角度的碰撞會(huì)導(dǎo)致沉箱不同的碰撞破壞模式。

    (2)所模擬的3種工況中波浪作用方向與x向夾角成60°時(shí),沉箱碰撞破壞模式與“莎莉嘉”臺(tái)風(fēng)作用期間的現(xiàn)場情況更為接近,沉箱外壁未發(fā)生明顯破壞,內(nèi)隔板出現(xiàn)不同程度的裂縫破損。

    猜你喜歡
    沉箱隔板夾角
    阿爾及利亞BETHIOUA礦業(yè)碼頭工程沉箱安裝施工技術(shù)
    基于有限元分析的隔板工藝支撐技術(shù)研究
    沉箱出運(yùn)安裝施工技術(shù)分析
    壓力容器隔板的一種設(shè)計(jì)方法
    探究鐘表上的夾角
    求解異面直線夾角問題的兩個(gè)路徑
    小型傾倒沉箱打撈扶正施工技術(shù)
    橫隔板參數(shù)對(duì)裝配式箱梁寬橋受力的影響
    任意夾角交叉封閉邊界內(nèi)平面流線計(jì)算及應(yīng)用
    雙頻帶隔板極化器
    電子測試(2017年15期)2017-12-18 07:18:51
    男女无遮挡免费网站观看| 80岁老熟妇乱子伦牲交| 精品少妇黑人巨大在线播放| 亚洲精品亚洲一区二区| 又大又黄又爽视频免费| 色视频在线一区二区三区| 色婷婷久久久亚洲欧美| 女性被躁到高潮视频| 亚洲av欧美aⅴ国产| 日本免费在线观看一区| 国产精品一区二区在线观看99| 久久女婷五月综合色啪小说| 在线观看一区二区三区激情| 午夜免费观看性视频| 国产无遮挡羞羞视频在线观看| 99久久精品国产国产毛片| 久久97久久精品| 国产免费视频播放在线视频| 国产免费视频播放在线视频| 日韩欧美精品免费久久| 高清不卡的av网站| 中文字幕av电影在线播放| 久久国产亚洲av麻豆专区| 精品少妇内射三级| 在现免费观看毛片| 少妇熟女欧美另类| 国产深夜福利视频在线观看| 亚洲欧美色中文字幕在线| 欧美变态另类bdsm刘玥| 欧美激情极品国产一区二区三区 | 乱码一卡2卡4卡精品| 中文字幕久久专区| 国产伦理片在线播放av一区| 黑丝袜美女国产一区| 日本色播在线视频| 91在线精品国自产拍蜜月| 热re99久久精品国产66热6| 久久久欧美国产精品| a级毛片在线看网站| 日本猛色少妇xxxxx猛交久久| 2021少妇久久久久久久久久久| 18禁动态无遮挡网站| 丝袜喷水一区| 国产精品不卡视频一区二区| 精品酒店卫生间| 高清黄色对白视频在线免费看| 欧美日韩成人在线一区二区| 99热网站在线观看| 大码成人一级视频| 伦理电影免费视频| 国产免费视频播放在线视频| 亚洲国产色片| 蜜桃国产av成人99| 午夜福利视频在线观看免费| a级片在线免费高清观看视频| 久久久久久久久久久免费av| 两个人的视频大全免费| 哪个播放器可以免费观看大片| 一级毛片 在线播放| av天堂久久9| 麻豆精品久久久久久蜜桃| 熟女av电影| 国产极品天堂在线| 色视频在线一区二区三区| 国产精品久久久久久精品古装| 人人妻人人爽人人添夜夜欢视频| 亚洲天堂av无毛| 啦啦啦中文免费视频观看日本| 国产精品久久久久久久久免| 欧美激情极品国产一区二区三区 | 啦啦啦中文免费视频观看日本| 香蕉精品网在线| 国产一区二区三区综合在线观看 | 国产成人免费无遮挡视频| 晚上一个人看的免费电影| 不卡视频在线观看欧美| 成人国产av品久久久| 成年人免费黄色播放视频| 日本午夜av视频| 久久97久久精品| 亚洲美女搞黄在线观看| 成年美女黄网站色视频大全免费 | 中国美白少妇内射xxxbb| 日产精品乱码卡一卡2卡三| 久久影院123| 三级国产精品片| 国产成人一区二区在线| 国产伦精品一区二区三区视频9| 丝袜美足系列| 免费黄频网站在线观看国产| 边亲边吃奶的免费视频| 观看美女的网站| 18禁动态无遮挡网站| 亚洲av二区三区四区| 日日撸夜夜添| 国产精品嫩草影院av在线观看| 国产免费又黄又爽又色| 国产极品天堂在线| 18禁在线无遮挡免费观看视频| 精品亚洲乱码少妇综合久久| 在线看a的网站| 色哟哟·www| 亚洲五月色婷婷综合| 国产精品不卡视频一区二区| 看非洲黑人一级黄片| 一区二区三区免费毛片| 一个人免费看片子| 欧美亚洲 丝袜 人妻 在线| 久久人人爽av亚洲精品天堂| 亚洲人成网站在线播| 久久久久久久久久人人人人人人| 国产精品人妻久久久影院| 国产老妇伦熟女老妇高清| 久久久精品免费免费高清| 午夜日本视频在线| 一二三四中文在线观看免费高清| 亚洲精品一区蜜桃| 大话2 男鬼变身卡| 久久青草综合色| 亚洲美女搞黄在线观看| 欧美精品一区二区大全| 免费观看无遮挡的男女| 精品一区二区三卡| 你懂的网址亚洲精品在线观看| 日韩不卡一区二区三区视频在线| 少妇丰满av| 亚洲av男天堂| 国产精品人妻久久久影院| 午夜免费鲁丝| 国产精品99久久99久久久不卡 | 国产高清三级在线| 伊人久久精品亚洲午夜| a级毛片在线看网站| 日韩中文字幕视频在线看片| av卡一久久| 国产一级毛片在线| 两个人免费观看高清视频| 国产av码专区亚洲av| av福利片在线| 亚洲国产精品一区三区| 亚洲欧美一区二区三区国产| a级片在线免费高清观看视频| 亚洲av国产av综合av卡| 国产一区亚洲一区在线观看| 丰满迷人的少妇在线观看| 亚洲国产最新在线播放| 三上悠亚av全集在线观看| 国产不卡av网站在线观看| 久久99一区二区三区| 一级毛片 在线播放| 嫩草影院入口| 国产一级毛片在线| 国产成人午夜福利电影在线观看| 高清黄色对白视频在线免费看| 亚洲综合色惰| 卡戴珊不雅视频在线播放| 99热6这里只有精品| 国产色婷婷99| 国产极品天堂在线| 毛片一级片免费看久久久久| 美女中出高潮动态图| 欧美日韩成人在线一区二区| 丰满少妇做爰视频| 久久久久久久久久成人| 国产伦精品一区二区三区视频9| 老司机影院成人| 一本一本久久a久久精品综合妖精 国产伦在线观看视频一区 | 蜜桃在线观看..| 婷婷色综合大香蕉| 亚洲精品日本国产第一区| 一本一本久久a久久精品综合妖精 国产伦在线观看视频一区 | 亚洲性久久影院| 日本wwww免费看| 色婷婷久久久亚洲欧美| 免费大片黄手机在线观看| 久久久久视频综合| 如何舔出高潮| 亚洲av日韩在线播放| 青春草亚洲视频在线观看| 99热这里只有精品一区| 视频在线观看一区二区三区| av卡一久久| 在线观看免费视频网站a站| av播播在线观看一区| 精品少妇内射三级| 两个人免费观看高清视频| 久久久久精品久久久久真实原创| 黄色视频在线播放观看不卡| 女人久久www免费人成看片| 午夜激情久久久久久久| 内地一区二区视频在线| 亚洲激情五月婷婷啪啪| 午夜福利影视在线免费观看| 久久久久久久久久成人| 国产精品免费大片| 高清在线视频一区二区三区| 免费人妻精品一区二区三区视频| 日韩欧美精品免费久久| 熟女av电影| 国产成人免费观看mmmm| 丰满饥渴人妻一区二区三| 大陆偷拍与自拍| 丰满乱子伦码专区| 午夜福利视频在线观看免费| 欧美日韩视频精品一区| 亚洲国产精品成人久久小说| 黄色一级大片看看| 久久午夜综合久久蜜桃| 欧美一级a爱片免费观看看| 男人操女人黄网站| 最近2019中文字幕mv第一页| 美女脱内裤让男人舔精品视频| 日韩强制内射视频| av黄色大香蕉| 国产一区二区在线观看日韩| 好男人视频免费观看在线| 丝瓜视频免费看黄片| 在线精品无人区一区二区三| 国产精品国产三级专区第一集| 精品一区在线观看国产| 久久精品国产鲁丝片午夜精品| 天天躁夜夜躁狠狠久久av| 91精品伊人久久大香线蕉| 一区二区日韩欧美中文字幕 | 亚洲精品乱久久久久久| 国产精品久久久久久久久免| 日韩一本色道免费dvd| 久久ye,这里只有精品| 啦啦啦啦在线视频资源| 国产精品一区二区在线观看99| 亚洲av男天堂| av电影中文网址| 久久久久久久久久久丰满| 波野结衣二区三区在线| 日韩不卡一区二区三区视频在线| 99久久综合免费| 涩涩av久久男人的天堂| 插逼视频在线观看| 色吧在线观看| 最近手机中文字幕大全| 一个人免费看片子| 三级国产精品欧美在线观看| 日本猛色少妇xxxxx猛交久久| 午夜激情福利司机影院| 亚洲无线观看免费| 制服诱惑二区| 啦啦啦视频在线资源免费观看| 久久久久久久久久久丰满| 亚洲精品久久久久久婷婷小说| 最近手机中文字幕大全| 国产成人精品在线电影| 国精品久久久久久国模美| 亚洲精品久久成人aⅴ小说 | 国产精品久久久久久久久免| 亚洲一区二区三区欧美精品| 亚洲精品日本国产第一区| videosex国产| 亚洲色图综合在线观看| 一本一本久久a久久精品综合妖精 国产伦在线观看视频一区 | 又黄又爽又刺激的免费视频.| 亚洲精品视频女| 精品国产国语对白av| 91精品三级在线观看| 一级片'在线观看视频| 亚洲图色成人| 如日韩欧美国产精品一区二区三区 | 最后的刺客免费高清国语| 日日摸夜夜添夜夜添av毛片| 高清av免费在线| 亚洲美女搞黄在线观看| 国产一区二区在线观看日韩| 欧美3d第一页| 日韩av在线免费看完整版不卡| 一区二区三区乱码不卡18| 精品人妻熟女毛片av久久网站| 精品久久久久久久久av| 欧美一级a爱片免费观看看| 草草在线视频免费看| 国产精品国产三级国产专区5o| 亚洲欧美一区二区三区黑人 | 天美传媒精品一区二区| 99久久精品一区二区三区| 亚洲伊人久久精品综合| 超碰97精品在线观看| 国产永久视频网站| 曰老女人黄片| 成人亚洲欧美一区二区av| 在线观看美女被高潮喷水网站| 日韩一区二区视频免费看| 99热国产这里只有精品6| 免费观看的影片在线观看| 蜜桃在线观看..| 国产一区二区在线观看av| 亚洲欧美中文字幕日韩二区| 国产又色又爽无遮挡免| 免费不卡的大黄色大毛片视频在线观看| 欧美日韩亚洲高清精品| 黄片无遮挡物在线观看| 99热全是精品| 欧美精品一区二区大全| av女优亚洲男人天堂| 大陆偷拍与自拍| 黄色毛片三级朝国网站| av有码第一页| 国产黄片视频在线免费观看| 亚洲精品一二三| 亚洲欧洲日产国产| 午夜影院在线不卡| 色吧在线观看| av在线老鸭窝| 九九久久精品国产亚洲av麻豆| 成人二区视频| 在线观看免费日韩欧美大片 | 国产熟女午夜一区二区三区 | 国产一区二区在线观看av| 香蕉精品网在线| 最新的欧美精品一区二区| a级毛片黄视频| 夜夜看夜夜爽夜夜摸| 精品国产国语对白av| 国产精品一区二区三区四区免费观看| 中文字幕av电影在线播放| 在线精品无人区一区二区三| 丰满乱子伦码专区| 亚洲国产欧美日韩在线播放| videosex国产| 久久久久视频综合| 国产乱来视频区| 18+在线观看网站| 欧美日韩视频精品一区| 人妻人人澡人人爽人人| 秋霞伦理黄片| 国产精品一国产av| av女优亚洲男人天堂| 欧美三级亚洲精品| 免费日韩欧美在线观看| 国产精品久久久久久久电影| 日韩免费高清中文字幕av| 免费高清在线观看日韩| 精品卡一卡二卡四卡免费| 女性生殖器流出的白浆| 亚洲婷婷狠狠爱综合网| 亚洲人成77777在线视频| 卡戴珊不雅视频在线播放| 日韩,欧美,国产一区二区三区| 在线观看国产h片| 久久这里有精品视频免费| 大码成人一级视频| 在线观看免费视频网站a站| 久久久国产一区二区| 久久久久久久亚洲中文字幕| 18禁裸乳无遮挡动漫免费视频| 91成人精品电影| 国产一区二区在线观看av| 亚洲精品乱码久久久v下载方式| 久久精品久久久久久噜噜老黄| 18禁动态无遮挡网站| 国产熟女午夜一区二区三区 | 女人精品久久久久毛片| 国产高清三级在线| 高清av免费在线| 免费人成在线观看视频色| 国产国拍精品亚洲av在线观看| 国产一区有黄有色的免费视频| 国产成人91sexporn| 伊人久久国产一区二区| 国产成人91sexporn| av黄色大香蕉| 亚洲一级一片aⅴ在线观看| 丁香六月天网| 中国三级夫妇交换| 日韩电影二区| 国产精品欧美亚洲77777| 日韩欧美一区视频在线观看| 制服人妻中文乱码| 免费看光身美女| 国产黄片视频在线免费观看| 国产白丝娇喘喷水9色精品| 在线观看免费高清a一片| 久久久国产精品麻豆| 亚洲成人手机| 欧美日韩综合久久久久久| 美女内射精品一级片tv| 全区人妻精品视频| 曰老女人黄片| 亚洲一区二区三区欧美精品| 婷婷色av中文字幕| 观看av在线不卡| 亚洲人成网站在线播| 国产精品秋霞免费鲁丝片| 自线自在国产av| 日本vs欧美在线观看视频| 午夜福利视频在线观看免费| 2018国产大陆天天弄谢| 国产精品人妻久久久久久| 欧美性感艳星| av有码第一页| 能在线免费看毛片的网站| 黑丝袜美女国产一区| 婷婷色综合www| 满18在线观看网站| 秋霞在线观看毛片| 成人黄色视频免费在线看| 三上悠亚av全集在线观看| 日本wwww免费看| 精品亚洲成国产av| 国产男女内射视频| 熟女av电影| 久久99热这里只频精品6学生| 亚洲精品美女久久av网站| 黄色毛片三级朝国网站| 你懂的网址亚洲精品在线观看| 在现免费观看毛片| 18在线观看网站| 日本黄色片子视频| 日韩不卡一区二区三区视频在线| tube8黄色片| 欧美精品国产亚洲| a级片在线免费高清观看视频| 国产成人av激情在线播放 | 久久ye,这里只有精品| 国产毛片在线视频| 中文字幕av电影在线播放| 国产男女超爽视频在线观看| 在线 av 中文字幕| 只有这里有精品99| 中文字幕亚洲精品专区| av免费在线看不卡| 久热这里只有精品99| 国产视频首页在线观看| 欧美激情极品国产一区二区三区 | av在线app专区| 一本—道久久a久久精品蜜桃钙片| 国产成人精品一,二区| 97在线视频观看| 人妻一区二区av| 一区二区三区四区激情视频| 大陆偷拍与自拍| 五月天丁香电影| 婷婷色麻豆天堂久久| 欧美丝袜亚洲另类| 极品人妻少妇av视频| 两个人免费观看高清视频| 激情五月婷婷亚洲| 免费看光身美女| 国产免费一区二区三区四区乱码| 纵有疾风起免费观看全集完整版| 性色av一级| 最黄视频免费看| 国产在线一区二区三区精| 大码成人一级视频| 国产伦理片在线播放av一区| 国产成人免费无遮挡视频| 你懂的网址亚洲精品在线观看| 国产精品一区www在线观看| 91国产中文字幕| 九九爱精品视频在线观看| 国产欧美另类精品又又久久亚洲欧美| 91久久精品电影网| 王馨瑶露胸无遮挡在线观看| 五月开心婷婷网| 久久精品久久精品一区二区三区| 丝袜喷水一区| 亚洲伊人久久精品综合| 亚洲图色成人| 日韩精品有码人妻一区| 欧美日韩亚洲高清精品| 观看av在线不卡| 啦啦啦啦在线视频资源| 美女国产视频在线观看| 久久久精品免费免费高清| 天美传媒精品一区二区| 你懂的网址亚洲精品在线观看| 国产成人aa在线观看| 91aial.com中文字幕在线观看| 九九爱精品视频在线观看| 777米奇影视久久| av专区在线播放| xxxhd国产人妻xxx| 久久人妻熟女aⅴ| 亚洲欧美一区二区三区国产| 久久精品国产自在天天线| 欧美精品高潮呻吟av久久| 狂野欧美激情性xxxx在线观看| 伊人久久精品亚洲午夜| 2021少妇久久久久久久久久久| 国产极品天堂在线| 精品亚洲乱码少妇综合久久| 欧美另类一区| 欧美日本中文国产一区发布| 大香蕉97超碰在线| 欧美精品一区二区免费开放| 老女人水多毛片| 国产亚洲最大av| 成人免费观看视频高清| a 毛片基地| 国产熟女午夜一区二区三区 | 国产精品一区二区在线不卡| .国产精品久久| 日韩 亚洲 欧美在线| av免费观看日本| 欧美精品高潮呻吟av久久| 亚洲av电影在线观看一区二区三区| 国产精品一区二区三区四区免费观看| 女性被躁到高潮视频| 国产一区二区三区av在线| 中文字幕人妻熟人妻熟丝袜美| a级毛片免费高清观看在线播放| av有码第一页| 精品人妻熟女av久视频| 少妇被粗大猛烈的视频| 日韩精品免费视频一区二区三区 | 国产精品久久久久成人av| 亚洲一区二区三区欧美精品| 久久久久久久久久久丰满| 欧美少妇被猛烈插入视频| 女人久久www免费人成看片| 美女国产高潮福利片在线看| 少妇被粗大猛烈的视频| 久久婷婷青草| 卡戴珊不雅视频在线播放| 亚洲av欧美aⅴ国产| 午夜精品国产一区二区电影| 精品久久久久久久久av| 性高湖久久久久久久久免费观看| 日韩亚洲欧美综合| 亚洲熟女精品中文字幕| 成人综合一区亚洲| 亚洲精品av麻豆狂野| 国产成人精品久久久久久| 啦啦啦啦在线视频资源| 一本色道久久久久久精品综合| 在线播放无遮挡| 性色avwww在线观看| 超色免费av| 中文字幕人妻熟人妻熟丝袜美| 中文字幕最新亚洲高清| 国产有黄有色有爽视频| 欧美丝袜亚洲另类| 国产精品人妻久久久影院| 久久精品国产亚洲av涩爱| 一本一本综合久久| 日韩中字成人| 久热这里只有精品99| 丝袜美足系列| 国产精品国产三级专区第一集| 精品人妻一区二区三区麻豆| 欧美性感艳星| 寂寞人妻少妇视频99o| 欧美 日韩 精品 国产| 免费久久久久久久精品成人欧美视频 | 亚洲三级黄色毛片| 一二三四中文在线观看免费高清| 男女啪啪激烈高潮av片| av在线app专区| 青青草视频在线视频观看| 大话2 男鬼变身卡| 狠狠精品人妻久久久久久综合| 夫妻午夜视频| 最近中文字幕2019免费版| 99热这里只有是精品在线观看| 亚洲精品乱久久久久久| 日本黄色日本黄色录像| 黄色视频在线播放观看不卡| 国产一区二区在线观看av| 日韩成人伦理影院| 亚洲四区av| 26uuu在线亚洲综合色| 欧美变态另类bdsm刘玥| 人妻人人澡人人爽人人| 男人操女人黄网站| 一区二区日韩欧美中文字幕 | 在线免费观看不下载黄p国产| 美女国产视频在线观看| 国产精品99久久久久久久久| 亚洲国产色片| 国产一区二区在线观看av| 亚洲国产精品999| 在线 av 中文字幕| 精品国产一区二区久久| 亚洲无线观看免费| a级毛片黄视频| 国产伦精品一区二区三区视频9| 欧美激情 高清一区二区三区| 日本欧美国产在线视频| 视频区图区小说| 热99久久久久精品小说推荐| 我要看黄色一级片免费的| av福利片在线| 两个人免费观看高清视频| 国语对白做爰xxxⅹ性视频网站| 国产精品一国产av| 一区二区日韩欧美中文字幕 | 国产日韩欧美视频二区| 日韩中文字幕视频在线看片| 五月开心婷婷网| 久久人人爽人人片av| 桃花免费在线播放| 久久精品国产亚洲av涩爱| 九色亚洲精品在线播放| 成人国产av品久久久| 91精品国产国语对白视频| 久久久午夜欧美精品| 日本欧美视频一区| 精品少妇内射三级| 欧美日韩在线观看h| 久久精品熟女亚洲av麻豆精品| 中文字幕久久专区| 精品国产乱码久久久久久小说| 亚洲精品日本国产第一区| 国产黄色免费在线视频| 午夜精品国产一区二区电影| 精品久久蜜臀av无| 亚洲欧美色中文字幕在线| 丰满乱子伦码专区|