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

    風(fēng)屏障對(duì)橋梁及列車的氣動(dòng)特性影響研究

    2022-07-29 02:35:46董國朝黃佳穎韓艷
    交通科學(xué)與工程 2022年2期
    關(guān)鍵詞:橋梁

    董國朝,黃佳穎,韓艷

    (長沙理工大學(xué) 土木工程學(xué)院,湖南 長沙 410114)

    隨著基礎(chǔ)建設(shè)的快速發(fā)展,交通運(yùn)輸體系外擴(kuò)[1],中國規(guī)劃建設(shè)的大型運(yùn)輸鐵路網(wǎng)覆蓋了各種強(qiáng)風(fēng)場地區(qū)。在高速鐵路上建設(shè)橋梁,可提高空間利用率,降低混行風(fēng)險(xiǎn),優(yōu)化建設(shè)布局,增大建筑空間,但廣泛應(yīng)用會(huì)大幅度提高列車行駛高程,降低橋上列車行駛的安全性。為保證高速列車在鐵路橋梁上的行車安全性,通??紤]設(shè)置風(fēng)屏障來降低風(fēng)速[2],但風(fēng)屏障的存在增大了橋梁的擋風(fēng)面積,對(duì)車橋系統(tǒng)的流場產(chǎn)生了干擾,影響了整體的抗風(fēng)設(shè)計(jì)[3]。因此,考慮在橋梁上加裝風(fēng)屏障時(shí),應(yīng)設(shè)置合理的風(fēng)屏障高度和透風(fēng)率,以達(dá)到最優(yōu)效果。目前,國內(nèi)外針對(duì)風(fēng)屏障高度和透風(fēng)率的設(shè)置研究主要是數(shù)值模擬和風(fēng)洞試驗(yàn)。TELENTA 等人[4]通過試驗(yàn)和數(shù)值模擬方法,研究了風(fēng)屏障擋風(fēng)角度調(diào)整與不同湍流流動(dòng)的機(jī)理,但僅考慮了風(fēng)屏障對(duì)橋梁的影響。LEE 等人[5]討論了風(fēng)屏障的設(shè)計(jì)準(zhǔn)則,在側(cè)風(fēng)作用下如何提高高速公路運(yùn)行車輛安全性的方法。施成華等人[6]通過建立三維風(fēng)-車-橋模型,進(jìn)行了動(dòng)態(tài)數(shù)值模擬。通過設(shè)置橋梁風(fēng)屏障后,對(duì)列車駛?cè)胨淼肋^程中的行車安全性和氣動(dòng)性能進(jìn)行研究。劉葉等人[7]對(duì)側(cè)風(fēng)作用下的平層公鐵橋梁-列車-風(fēng)屏障系統(tǒng)氣動(dòng)特性進(jìn)行了風(fēng)洞試驗(yàn)研究,針對(duì)風(fēng)屏障不同的透風(fēng)率和高度,對(duì)不同偏角下橋上中間列車的三分力系數(shù)進(jìn)行了測定,研究了風(fēng)屏障參數(shù)問題。周蕾等人[3]利用計(jì)算流體動(dòng)力學(xué)(computational fluid dynamics,簡稱為CFD)進(jìn)行了數(shù)值模擬和風(fēng)洞試驗(yàn),從風(fēng)屏障的結(jié)構(gòu)選型和高度等角度討論了對(duì)車橋系統(tǒng)氣動(dòng)性能的影響,但未考慮列車頭的流線型特征,未能較好地得到實(shí)際流場分布,其數(shù)值模擬的精度也不夠準(zhǔn)確。這些研究成果均表明:風(fēng)屏障提高了鐵路橋梁上列車行車安全性與舒適性。但針對(duì)風(fēng)屏障高度、透風(fēng)率等對(duì)列車不同部位影響的研究較少,從車橋整體氣動(dòng)性能優(yōu)化的角度,得到最佳風(fēng)屏障參數(shù)取值的研究更為鮮見。因此,本研究擬通過數(shù)值模擬和風(fēng)洞試驗(yàn)方法,對(duì)側(cè)風(fēng)下平層公鐵兩用橋梁的車-橋-風(fēng)屏障系統(tǒng)氣動(dòng)性能進(jìn)行研究。通過對(duì)不同的風(fēng)屏障高度和透風(fēng)率進(jìn)行數(shù)值模擬,探究風(fēng)屏障參數(shù)對(duì)實(shí)際橋梁動(dòng)力特性的影響,并與風(fēng)洞試驗(yàn)結(jié)果進(jìn)行比較,驗(yàn)證模型的準(zhǔn)確性,再通過數(shù)據(jù)包絡(luò)分析(data envelopment analysis,簡稱為DEA)計(jì)算不同風(fēng)屏障參數(shù)下對(duì)應(yīng)的車橋整體相對(duì)效率值,給出該類車橋下風(fēng)屏障參數(shù)的最佳取值,為實(shí)際工程提供借鑒。

    1 數(shù)值模擬參數(shù)設(shè)置

    1.1 模型參數(shù)

    對(duì)某平層公鐵兩用大橋和高速列車進(jìn)行了多工況的風(fēng)洞試驗(yàn),為確保準(zhǔn)確性,數(shù)值模擬中保留公路與鐵路之間設(shè)置的防護(hù)欄。數(shù)值模擬中,車橋模型與風(fēng)洞試驗(yàn)中均采用1∶43 的幾何縮尺比??紤]到風(fēng)洞試驗(yàn)布置和數(shù)值模擬中網(wǎng)格劃分的可行性,對(duì)試驗(yàn)?zāi)P瓦M(jìn)行了部分簡化,去除了列車原型的列車輪對(duì)、轉(zhuǎn)向架和受電弓等列車細(xì)部構(gòu)造,同時(shí)忽略了橋梁上檢修軌道、拉索等附屬設(shè)施。

    縮尺后的列車模型尺寸為長1 775.8 mm×寬78.6 mm×高81.4 mm,其三維模型如圖1 所示。列車的流線型車頭對(duì)車橋系統(tǒng)的流場影響不可忽略。因此,在風(fēng)洞試驗(yàn)和數(shù)值模擬中,完整保留頭車和尾車的幾何外形。列車采用頭車、中車和尾車3節(jié)列車模擬等效整個(gè)列車。外形相對(duì)比較規(guī)則的中間車輛,這種模型氣動(dòng)性較為穩(wěn)定,能較好地反映整個(gè)中間車節(jié)。

    圖1 列車簡化模型(單位:mm)Fig.1 Simplified model of train(unit:mm)

    簡化后的橋梁保留了對(duì)整體氣動(dòng)特性影響較大的橋上欄桿,設(shè)置列車底部與橋梁頂面距離為0.5 m 來等效輪對(duì)高度。車橋整體三維模型如圖2所示。

    圖2 公鐵兩用橋梁模型Fig.2 The mode of combined highway and railway bridge

    采用等效透風(fēng)率法對(duì)風(fēng)屏障進(jìn)行等效處理[4-7]。采用多塊單板組成風(fēng)屏障,如圖3 所示。其中,s為風(fēng)屏障條的間隙,h為單個(gè)風(fēng)屏障條的高度。風(fēng)屏障高度(H)定義為最頂部單板頂面距離橋面的高度,單板之間的距離(S)相等,透風(fēng)率β=5S/H×100%。

    圖3 風(fēng)屏障等效模型Fig.3 Equivalent model of wind barrier

    1.2 工況設(shè)置

    為研究不同風(fēng)屏障高度和透風(fēng)率對(duì)車橋系統(tǒng)氣動(dòng)性能的影響,本試驗(yàn)分別對(duì)風(fēng)屏障高度為0、2.0、3.0、3.5 和4.0 m,透風(fēng)率為20%、40%和60%的13 種不同工況進(jìn)行數(shù)值模擬,測定對(duì)應(yīng)的三分力系數(shù)。

    2 數(shù)值方法和網(wǎng)格劃分

    根據(jù)實(shí)際流場性質(zhì),采用三維、非定常、不可壓縮和黏性流場對(duì)風(fēng)屏障作用下的車橋系統(tǒng)進(jìn)行數(shù)值模擬。

    由于本研究的計(jì)算對(duì)象是車橋,重點(diǎn)考慮逆壓梯度分離引起的力的變化。因此,選用考慮了逆壓梯度分離的SST k-w 模型??紤]到SST k-w 湍流模型對(duì)近壁面區(qū)域網(wǎng)格細(xì)化的要求[8],一般控制Yplus值在1左右,確保達(dá)到湍流應(yīng)用條件。

    幾何模型的網(wǎng)格在ICEM CFD 中進(jìn)行,考慮計(jì)算資源和精度,本試驗(yàn)對(duì)整體結(jié)構(gòu)采用全結(jié)構(gòu)六面體進(jìn)行網(wǎng)格劃分,而不采用精度差的混合網(wǎng)格劃分方法,對(duì)列車流線型車頭進(jìn)行細(xì)化分塊,確保復(fù)雜表面網(wǎng)格的質(zhì)量。考慮到高速列車的車身表面復(fù)雜性,將整個(gè)計(jì)算域分為高速列車區(qū)域和橋梁-風(fēng)障區(qū)域,以方便對(duì)2 個(gè)區(qū)域采用不同的切塊形式,對(duì)列車車身周圍網(wǎng)格進(jìn)行了內(nèi)外O 型切分,實(shí)現(xiàn)網(wǎng)格貼體和更好的網(wǎng)格尺寸過渡,能對(duì)近壁面的第一層網(wǎng)格高度進(jìn)行控制,滿足網(wǎng)格中Yplus值要求。

    選用全結(jié)構(gòu)化網(wǎng)格的另一個(gè)優(yōu)勢是可更好地控制網(wǎng)格質(zhì)量,列車及橋梁-風(fēng)屏障部分網(wǎng)格如圖4 所示。高速列車區(qū)域的網(wǎng)格(80.96 萬)質(zhì)量最小值為0.710,橋梁風(fēng)障區(qū)域的網(wǎng)格(481.80 萬)質(zhì)量最小值為0.941,確保了數(shù)值模擬的準(zhǔn)確度。

    圖4 全結(jié)構(gòu)網(wǎng)格劃分示意Fig.4 Signal of full structure grid

    計(jì)算域的設(shè)置直接影響數(shù)值模擬結(jié)果的準(zhǔn)確性,為確保數(shù)值模擬準(zhǔn)確性與計(jì)算效率,采用的計(jì)算域如圖5所示。由于不考慮風(fēng)偏角對(duì)車橋系統(tǒng)的影響,本數(shù)值模擬采用單進(jìn)風(fēng)口,設(shè)置迎風(fēng)側(cè)流場邊界為速度入口,背風(fēng)側(cè)流場邊界為壓力出口,設(shè)置橋梁兩端為對(duì)稱邊界,其余均設(shè)置為壁面邊界。

    圖5 計(jì)算域示意(單位:m)Fig.5 Schematic diagram of computational domain(unit:m)

    3 風(fēng)洞試驗(yàn)驗(yàn)證

    為驗(yàn)證數(shù)值模擬的準(zhǔn)確性,風(fēng)洞試驗(yàn)和數(shù)值模擬選用的橋梁節(jié)段和列車模型尺寸一致,均保持與原車橋1∶43的縮尺比,在長沙理工大學(xué)的風(fēng)洞高速試驗(yàn)段進(jìn)行模型試驗(yàn)。本試驗(yàn)?zāi)M中采用的風(fēng)速為10.6 m/s 的均勻來流,湍流強(qiáng)度控制在0.5%以內(nèi)。

    模型為某平層公鐵兩用大橋,橋梁中間為雙線鐵路,線間距為4.6 m,鐵路兩邊是雙向6 車道的高速公路,公路與鐵路之間設(shè)防護(hù)欄。主梁斷面如圖6 所示,橋梁三分力通過測力天平得到數(shù)據(jù)。

    圖6 主梁斷面Fig.6 The main beam section

    列車原型為CRH2型客車,列車模型尺寸為長2 300 mm×寬78.6 mm×高81.4 mm。在橫風(fēng)條件下,考慮列車風(fēng)對(duì)橋梁和風(fēng)屏障的影響,為更好地捕捉頭車和尾車的復(fù)雜流場,保留車頭部分完整的細(xì)部特征,確保與數(shù)值模擬的一致性。列車測壓斷面和測壓孔布置如圖7所示,采用掃描閥連接進(jìn)行數(shù)據(jù)采集。

    圖7 列車測壓孔位布置Fig.7 Layout of holes for train pressure measuring

    風(fēng)洞試驗(yàn)中,采用等效透風(fēng)率法,選擇2.75 m高、40%透風(fēng)率作為試驗(yàn)的風(fēng)屏障參數(shù)取值??紤]到不同風(fēng)攻角對(duì)應(yīng)的三分力系數(shù)呈現(xiàn)不同變化趨勢,試驗(yàn)中,通過轉(zhuǎn)動(dòng)底部轉(zhuǎn)盤實(shí)現(xiàn)不同的風(fēng)攻角,從而測定-1°、0°和+1°不同風(fēng)攻角對(duì)應(yīng)的三分力系數(shù)。

    橋梁采用2臺(tái)六分量應(yīng)變天平測量氣動(dòng)力和氣動(dòng)力矩,列車采用原裝進(jìn)口美國PSI DTC Initium電子壓力掃描閥,512 個(gè)通道自由組合供選用(64通道/塊×8 塊),各壓力測量通道精度為0.06%;參考點(diǎn)的風(fēng)速由皮托管測量得到。試驗(yàn)中采樣時(shí)長為30 s,采樣頻率為330 Hz。

    數(shù)值模擬的三分力系數(shù)與風(fēng)洞試驗(yàn)結(jié)果對(duì)比如圖8所示。其中,扭矩系數(shù)隨風(fēng)攻角變化的趨勢與試驗(yàn)監(jiān)測基本一致,但由于升力矩系數(shù)在小攻角下數(shù)值較小,因此該誤差相對(duì)較大。試驗(yàn)測定的列車頭部阻力系數(shù)和升力系數(shù)與數(shù)值模擬基本吻合,最大誤差僅為8.86%(0°攻角下的列車阻力系數(shù)),試驗(yàn)結(jié)果表明:對(duì)車橋三分力的網(wǎng)格劃分和數(shù)值計(jì)算是可行的。

    圖8 三分力系數(shù)結(jié)果Fig.8 Three component coefficient results

    4 數(shù)值模擬結(jié)果分析

    4.1 風(fēng)屏障透風(fēng)率和高度對(duì)橋梁氣動(dòng)性能的影響

    不同風(fēng)屏障高度和透風(fēng)率下橋梁的三分力系數(shù)變化規(guī)律如圖9所示。

    從圖9 中可以看出,隨著風(fēng)屏障高度的增加,風(fēng)屏障增大了橋梁的擋風(fēng)面積,橋梁的阻力系數(shù)一直上升??疾旄鞣侄沃g的斜率,發(fā)現(xiàn)橋梁阻力隨風(fēng)屏障高度增大呈先逐漸快速增加后放緩,特別是在風(fēng)屏障高度為3.0~3.5 m 時(shí),阻力系數(shù)上升較為平緩,隨后在3.5~4.0 m 時(shí)上升趨勢明顯加快。透風(fēng)率對(duì)橋梁阻力系數(shù)的影響較為直觀,風(fēng)屏障高度相同的情況下,透風(fēng)率越高,阻力系數(shù)越小。

    圖9 不同風(fēng)屏障參數(shù)下的橋梁三分力系數(shù)Fig.9 Three component coefficient of bridge at different wind barrier parameters

    橋梁的升力系數(shù)在整組工況中為負(fù)值,隨著風(fēng)屏障高度的增加,先下降后上升,隨后再次下降,在風(fēng)屏障高度為3.0 m 時(shí)取極小值。由于本試驗(yàn)所有工況僅考慮0°攻角下三分力系數(shù)的變化,扭矩系數(shù)變化趨勢不大,在0值附近上下波動(dòng)。因此,扭矩系數(shù)保持在0值附近。其中,在風(fēng)屏障高度為3.5 m,透風(fēng)率為20% 時(shí),橋梁的扭矩系數(shù)與0差值最小。

    4.2 風(fēng)屏障透風(fēng)率和高度對(duì)列車氣動(dòng)性能的影響

    不同風(fēng)屏障高度和透風(fēng)率下列車的三分力系數(shù)變化規(guī)律如圖10所示。

    從圖10 中可以看出,橋面增加了風(fēng)屏障后,列車的氣動(dòng)性能得到明顯改善,阻力系數(shù)隨風(fēng)屏障高度升高而下降,3.5 m 之后已無明顯變化,這是由于此時(shí)的風(fēng)屏障高度基本與列車高度接近。在60%高透風(fēng)率下,隨著風(fēng)屏障高度的增加,列車的阻力系數(shù)在經(jīng)歷長段下降之后再次上升。這是因?yàn)楸3州^高的透風(fēng)率,增加風(fēng)屏障高度,遮蔽效果會(huì)變差,導(dǎo)致列車受到更大的橫風(fēng)干擾。升力系數(shù)與阻力系數(shù)變化趨勢類似。表明:高透風(fēng)率下風(fēng)屏障的遮蔽效果不明顯,扭矩系數(shù)隨風(fēng)屏障高度增加逐漸減??;低透風(fēng)率下風(fēng)屏障工況會(huì)有更好的防風(fēng)效果,在4.0 m時(shí)效果最好。

    圖10 不同風(fēng)屏障參數(shù)下的列車三分力系數(shù)Fig.10 Three component coefficients of bridge at different wind barrier parameters

    5 風(fēng)屏障參數(shù)優(yōu)化選型

    5.1 DEA評(píng)價(jià)模型

    DEA 算法引入一個(gè)相對(duì)效率作為參考指標(biāo)來進(jìn)行效率評(píng)價(jià),主要使用投入、產(chǎn)出指標(biāo)的決策單元(decision making unit,簡稱為DMU)對(duì)本研究中不同風(fēng)屏障參數(shù)取值進(jìn)行評(píng)價(jià)分析。通過數(shù)學(xué)規(guī)劃對(duì)各種決策單元之間的相對(duì)有效性進(jìn)行評(píng)價(jià)。由于DEA 方法對(duì)多指標(biāo)決策方案(單元)的評(píng)價(jià)具有很大優(yōu)勢,被廣泛應(yīng)用到很多領(lǐng)域[9-10]。

    假定有n個(gè)決策單元,每個(gè)決策單元有p種輸入和q種輸出形式,Xi和Yi為第i個(gè)風(fēng)屏障輸入及輸出指標(biāo)向量,決策單元的效率評(píng)價(jià)指數(shù)為:

    式中:U=(u1,u2,…,uq)T、V=(v1,v2,…,vp)T分別為輸出和輸入指標(biāo)對(duì)應(yīng)的權(quán)重系數(shù)。

    對(duì)應(yīng)的最優(yōu)化模型為:

    再經(jīng)過Charnes-Cooper變換和對(duì)偶理論的線性規(guī)劃變換,可得:

    式中:θ為第1 個(gè)風(fēng)屏障的相對(duì)效率;λi為變量;S-、S*分別為松弛變量。

    在風(fēng)洞試驗(yàn)中,得到風(fēng)屏障參數(shù)對(duì)橋梁和橋上列車的氣動(dòng)性能影響規(guī)律。考慮橋梁與列車各自三分力受風(fēng)屏障參數(shù)取值的影響不同,未得到風(fēng)屏障參數(shù)的合理取值。從橋梁設(shè)計(jì)角度出發(fā),考慮橋梁安全性能與列車安全性能,利用DEA 方法,利用數(shù)值模擬得到的數(shù)據(jù)進(jìn)行多目標(biāo)優(yōu)化評(píng)估。

    5.2 評(píng)估結(jié)果分析

    使用DEA 評(píng)價(jià)模型對(duì)不同工況下列車與橋梁的三分力系數(shù)進(jìn)行評(píng)估,結(jié)果見表1~3。

    表1 僅考慮橋梁時(shí)風(fēng)屏障參數(shù)取值的相對(duì)效率Table 1 Relative efficiency of wind barrier parameters under the consideration of bridge

    表2 僅考慮列車時(shí)風(fēng)屏障參數(shù)取值的相對(duì)效率Table 2 Relative efficiency of wind barrier parameters under the consideration of train

    表3 綜合考慮列車與橋梁時(shí)風(fēng)屏障參數(shù)取值的相對(duì)效率Table 3 Relative efficiency of wind barrier parameters under the comprehensive consideration of train and bridge

    由表1~3 可知,僅考慮橋梁安全性能時(shí),風(fēng)屏障參數(shù)的最佳取值為0.0 m,即不設(shè)置風(fēng)屏障。僅考慮列車時(shí),風(fēng)屏障參數(shù)中屏障的高度和透風(fēng)率的最佳取值分別為4.0 m 和20%。考慮車橋系統(tǒng)氣動(dòng)力系數(shù)隨風(fēng)屏障參數(shù)取值變化的影響,風(fēng)屏障高度為3.5 m、透風(fēng)率為20%時(shí),得到相對(duì)效率值最大(1.763 3)。多目標(biāo)評(píng)價(jià)的結(jié)果在單目標(biāo)評(píng)價(jià)結(jié)果(僅考慮車輛和僅考慮橋梁)之中,驗(yàn)證了多目標(biāo)評(píng)價(jià)的合理性。

    6 結(jié)論

    以某公鐵兩用橋?yàn)槔?,利用CFD 對(duì)不同風(fēng)屏障參數(shù)取值的車橋氣動(dòng)三分力進(jìn)行了數(shù)值模擬,并將所得成果與風(fēng)洞試驗(yàn)結(jié)果進(jìn)行對(duì)比,驗(yàn)證了數(shù)值模擬的準(zhǔn)確性。通過分析風(fēng)屏障高度和透風(fēng)率的變化對(duì)橋梁和橋上列車氣動(dòng)力的影響規(guī)律,得到結(jié)論為:

    1)風(fēng)屏障高度升高,透風(fēng)率減小,擋風(fēng)面積變大,橋梁的氣動(dòng)力系數(shù)逐漸增大。表明:橋梁氣動(dòng)力受風(fēng)屏障高度比透風(fēng)率的變化更為敏感,而風(fēng)屏障的存在為橋上列車提供了遮蔽效果,使橋上列車的氣動(dòng)力系數(shù)呈相反的變化趨勢。

    2)風(fēng)屏障透風(fēng)率較大時(shí),即使增加風(fēng)屏障高度,也達(dá)不到較好的防風(fēng)效果,但僅從單獨(dú)的三分力進(jìn)行評(píng)估,尚未得到合理的風(fēng)屏障參數(shù)取值。

    3)通過DEA算法,分別計(jì)算得到單獨(dú)考慮列車、橋梁和車橋整體氣動(dòng)性能下的風(fēng)屏障相對(duì)效率值。表明:風(fēng)屏障高度為3.5 m 且透風(fēng)率為20%時(shí),車橋整體的氣動(dòng)性能最優(yōu)。

    猜你喜歡
    橋梁
    一種橋梁伸縮縫防滲水裝置
    手拉手 共搭愛的橋梁
    句子也需要橋梁
    橋梁的加固研究
    江西建材(2018年4期)2018-04-10 12:37:32
    加固技術(shù)創(chuàng)新,為橋梁健康保駕護(hù)航
    中國公路(2017年11期)2017-07-31 17:56:30
    無人機(jī)在橋梁檢測中的應(yīng)用
    中國公路(2017年10期)2017-07-21 14:02:37
    高性能砼在橋梁中的應(yīng)用
    現(xiàn)代鋼橋制造對(duì)橋梁鋼的更高要求
    焊接(2016年8期)2016-02-27 13:05:15
    城鄉(xiāng)建設(shè)一體化要注重橋梁的建筑設(shè)計(jì)
    南昌54座橋梁進(jìn)行兩個(gè)月的夏季體檢
    久久久国产成人免费| 色精品久久人妻99蜜桃| 男人舔女人下体高潮全视频| 欧美绝顶高潮抽搐喷水| 国产单亲对白刺激| 成人无遮挡网站| 亚洲精品色激情综合| 两个人的视频大全免费| 国产精品一区二区免费欧美| a在线观看视频网站| 少妇的逼好多水| 2021天堂中文幕一二区在线观| 亚洲国产欧美网| 精品99又大又爽又粗少妇毛片 | 少妇丰满av| 人妻丰满熟妇av一区二区三区| 亚洲国产欧美人成| 美女大奶头视频| 动漫黄色视频在线观看| 18+在线观看网站| 成人18禁在线播放| 麻豆成人午夜福利视频| av在线天堂中文字幕| 美女大奶头视频| 国产视频一区二区在线看| 国产精品野战在线观看| 国产精品久久电影中文字幕| 18禁裸乳无遮挡免费网站照片| 亚洲无线在线观看| 男女下面进入的视频免费午夜| 久久香蕉精品热| 亚洲成av人片在线播放无| 日韩欧美国产在线观看| 成人av在线播放网站| 女生性感内裤真人,穿戴方法视频| 国内揄拍国产精品人妻在线| 老汉色∧v一级毛片| 日韩欧美在线乱码| 国产成人aa在线观看| av女优亚洲男人天堂| 日本 欧美在线| 欧美在线一区亚洲| 18禁在线播放成人免费| 中文字幕人妻丝袜一区二区| 国产蜜桃级精品一区二区三区| 搡老岳熟女国产| 3wmmmm亚洲av在线观看| 少妇人妻一区二区三区视频| 午夜福利高清视频| 变态另类成人亚洲欧美熟女| 在线播放无遮挡| 丁香六月欧美| 国产精品一区二区免费欧美| 久久久国产成人精品二区| 国产精品久久久久久人妻精品电影| 亚洲国产欧美人成| 深夜精品福利| 无遮挡黄片免费观看| 每晚都被弄得嗷嗷叫到高潮| 夜夜看夜夜爽夜夜摸| av天堂中文字幕网| 搡老熟女国产l中国老女人| 制服人妻中文乱码| 久久久久精品国产欧美久久久| 特大巨黑吊av在线直播| 99热6这里只有精品| 一本久久中文字幕| 亚洲五月天丁香| 岛国在线免费视频观看| 免费搜索国产男女视频| 精品久久久久久久毛片微露脸| 久久久久久久久大av| www日本黄色视频网| 嫩草影视91久久| 国产av在哪里看| 欧美大码av| 18禁在线播放成人免费| 琪琪午夜伦伦电影理论片6080| 叶爱在线成人免费视频播放| 九色成人免费人妻av| 亚洲精品色激情综合| 亚洲av美国av| 久久久久久久久大av| 美女黄网站色视频| 九色成人免费人妻av| 丰满的人妻完整版| 亚洲精品在线观看二区| 身体一侧抽搐| 国产亚洲精品综合一区在线观看| 国产激情偷乱视频一区二区| 成人一区二区视频在线观看| 一级毛片高清免费大全| 欧美最黄视频在线播放免费| 精品一区二区三区av网在线观看| 色尼玛亚洲综合影院| 亚洲一区高清亚洲精品| 成人av在线播放网站| 嫁个100分男人电影在线观看| 丰满乱子伦码专区| 岛国视频午夜一区免费看| av欧美777| 亚洲美女黄片视频| 午夜福利在线在线| 亚洲国产高清在线一区二区三| 亚洲成av人片在线播放无| 精品一区二区三区av网在线观看| 久久伊人香网站| 亚洲欧美精品综合久久99| 少妇的丰满在线观看| 真人一进一出gif抽搐免费| 久久久久国内视频| 国产精品女同一区二区软件 | 18禁裸乳无遮挡免费网站照片| 日本免费a在线| 乱人视频在线观看| 女生性感内裤真人,穿戴方法视频| 香蕉av资源在线| 久久亚洲精品不卡| 精品久久久久久,| 成人国产一区最新在线观看| 12—13女人毛片做爰片一| 又紧又爽又黄一区二区| 国产成+人综合+亚洲专区| 亚洲成av人片在线播放无| 一区二区三区高清视频在线| 网址你懂的国产日韩在线| 国产精品亚洲美女久久久| 在线国产一区二区在线| 国产色婷婷99| 精品久久久久久久久久久久久| 色综合欧美亚洲国产小说| 特大巨黑吊av在线直播| 国产视频内射| 女人被狂操c到高潮| 女人被狂操c到高潮| 性欧美人与动物交配| 国产成人欧美在线观看| 99久国产av精品| 成年女人毛片免费观看观看9| 欧美3d第一页| 亚洲久久久久久中文字幕| 成人18禁在线播放| 国产在视频线在精品| 婷婷丁香在线五月| 蜜桃久久精品国产亚洲av| 亚洲国产中文字幕在线视频| 蜜桃久久精品国产亚洲av| 中文字幕人妻熟人妻熟丝袜美 | 亚洲中文字幕日韩| 精品国产超薄肉色丝袜足j| 日韩中文字幕欧美一区二区| 黄色成人免费大全| 久久国产精品影院| 免费av不卡在线播放| 国产黄片美女视频| 在线天堂最新版资源| 亚洲专区中文字幕在线| 九九在线视频观看精品| 身体一侧抽搐| 成人特级黄色片久久久久久久| 不卡一级毛片| 成人国产一区最新在线观看| 国产成人系列免费观看| 少妇人妻一区二区三区视频| 国产高清视频在线观看网站| 免费高清视频大片| 变态另类丝袜制服| 狠狠狠狠99中文字幕| 俄罗斯特黄特色一大片| 中文字幕人成人乱码亚洲影| 亚洲成人免费电影在线观看| 国产视频一区二区在线看| 在线观看av片永久免费下载| 精品久久久久久,| 一a级毛片在线观看| 精品国产三级普通话版| 宅男免费午夜| 天堂√8在线中文| 亚洲无线在线观看| 欧美在线一区亚洲| 久久久久九九精品影院| 精品无人区乱码1区二区| svipshipincom国产片| 国产av一区在线观看免费| 亚洲在线自拍视频| 亚洲午夜理论影院| 国产精品一区二区三区四区免费观看 | 女人高潮潮喷娇喘18禁视频| 91av网一区二区| 日本精品一区二区三区蜜桃| 少妇的逼好多水| or卡值多少钱| 亚洲精品影视一区二区三区av| 午夜精品在线福利| 在线免费观看不下载黄p国产 | 白带黄色成豆腐渣| 国产精品98久久久久久宅男小说| 一级黄色大片毛片| 国产精品综合久久久久久久免费| 午夜福利18| 国内揄拍国产精品人妻在线| 久久精品国产亚洲av香蕉五月| 九色成人免费人妻av| 日本a在线网址| 国模一区二区三区四区视频| 欧美乱妇无乱码| 国产中年淑女户外野战色| 狂野欧美白嫩少妇大欣赏| 免费在线观看亚洲国产| 夜夜夜夜夜久久久久| 成人18禁在线播放| 白带黄色成豆腐渣| 日韩欧美国产一区二区入口| 日本与韩国留学比较| 一a级毛片在线观看| 久久精品国产亚洲av香蕉五月| 国产精品一及| 日韩欧美免费精品| 日韩欧美免费精品| 亚洲最大成人手机在线| 人妻久久中文字幕网| 麻豆一二三区av精品| 亚洲久久久久久中文字幕| 一级毛片女人18水好多| 欧美绝顶高潮抽搐喷水| 国产精品乱码一区二三区的特点| 日日摸夜夜添夜夜添小说| 国产成人啪精品午夜网站| 一卡2卡三卡四卡精品乱码亚洲| 蜜桃久久精品国产亚洲av| 啦啦啦观看免费观看视频高清| 一级毛片高清免费大全| 99久久九九国产精品国产免费| av在线蜜桃| 色综合亚洲欧美另类图片| 久久久精品欧美日韩精品| 日韩av在线大香蕉| 少妇人妻精品综合一区二区 | 亚洲av电影不卡..在线观看| 亚洲第一欧美日韩一区二区三区| 欧洲精品卡2卡3卡4卡5卡区| 国产精品电影一区二区三区| 香蕉久久夜色| 亚洲中文字幕一区二区三区有码在线看| 婷婷精品国产亚洲av在线| 日本 av在线| 少妇的逼好多水| av视频在线观看入口| 国产精品嫩草影院av在线观看 | 亚洲欧美日韩高清在线视频| 国产午夜精品论理片| 精品不卡国产一区二区三区| 日本黄大片高清| 97人妻精品一区二区三区麻豆| 中文字幕人妻丝袜一区二区| 淫秽高清视频在线观看| 日本一二三区视频观看| 99久久精品国产亚洲精品| 天堂网av新在线| 小说图片视频综合网站| 国产真实伦视频高清在线观看 | 尤物成人国产欧美一区二区三区| 欧美+日韩+精品| 亚洲中文字幕一区二区三区有码在线看| eeuss影院久久| 欧美性感艳星| 亚洲最大成人手机在线| 最近在线观看免费完整版| 欧美日韩亚洲国产一区二区在线观看| 国产午夜精品论理片| 婷婷精品国产亚洲av在线| 免费电影在线观看免费观看| 香蕉丝袜av| 一级作爱视频免费观看| av黄色大香蕉| 欧美不卡视频在线免费观看| 一级a爱片免费观看的视频| 色综合站精品国产| 亚洲成av人片免费观看| 成人亚洲精品av一区二区| 在线播放无遮挡| 国内精品美女久久久久久| 久久婷婷人人爽人人干人人爱| 久久天躁狠狠躁夜夜2o2o| 制服人妻中文乱码| 精品人妻偷拍中文字幕| 90打野战视频偷拍视频| 亚洲乱码一区二区免费版| 精品久久久久久,| 国产一区二区亚洲精品在线观看| 一区二区三区免费毛片| 欧美中文日本在线观看视频| 亚洲在线观看片| 在线a可以看的网站| 国产爱豆传媒在线观看| 搡女人真爽免费视频火全软件 | 久久久国产成人精品二区| 国内精品久久久久精免费| 精品福利观看| 亚洲美女黄片视频| 国产精品日韩av在线免费观看| 欧美黑人巨大hd| 亚洲成人久久爱视频| 中文在线观看免费www的网站| 日韩国内少妇激情av| 免费看a级黄色片| 久久久精品欧美日韩精品| 亚洲人与动物交配视频| 国产又黄又爽又无遮挡在线| 老熟妇乱子伦视频在线观看| 亚洲成人精品中文字幕电影| 听说在线观看完整版免费高清| 午夜精品久久久久久毛片777| 国产精品电影一区二区三区| 欧美日韩黄片免| 久久99热这里只有精品18| 免费看十八禁软件| 观看免费一级毛片| 又紧又爽又黄一区二区| 黄色女人牲交| 老鸭窝网址在线观看| 亚洲av不卡在线观看| 国产熟女xx| 久久久久亚洲av毛片大全| 国产精品三级大全| 一本一本综合久久| 真实男女啪啪啪动态图| 日本黄大片高清| 国产精品自产拍在线观看55亚洲| 在线看三级毛片| 在线播放无遮挡| 亚洲真实伦在线观看| 黄色视频,在线免费观看| 女人被狂操c到高潮| 熟妇人妻久久中文字幕3abv| 精品人妻1区二区| 超碰av人人做人人爽久久 | 老司机午夜十八禁免费视频| 国产午夜精品论理片| 国产精品综合久久久久久久免费| АⅤ资源中文在线天堂| 国产成人aa在线观看| 欧美高清成人免费视频www| 国产三级在线视频| 久久6这里有精品| 亚洲精品乱码久久久v下载方式 | 午夜福利免费观看在线| 好男人在线观看高清免费视频| 午夜精品一区二区三区免费看| 国产高潮美女av| 国产成人福利小说| 亚洲国产色片| 国产在视频线在精品| 亚洲专区国产一区二区| 草草在线视频免费看| 99久久无色码亚洲精品果冻| 亚洲男人的天堂狠狠| 亚洲最大成人中文| 91av网一区二区| 日日干狠狠操夜夜爽| 欧美成人免费av一区二区三区| 天天一区二区日本电影三级| 亚洲片人在线观看| 亚洲精品影视一区二区三区av| 蜜桃久久精品国产亚洲av| 国产乱人伦免费视频| 少妇丰满av| 男女视频在线观看网站免费| 精品日产1卡2卡| 特级一级黄色大片| 久久久久久国产a免费观看| 日韩亚洲欧美综合| 黑人欧美特级aaaaaa片| 欧美日韩综合久久久久久 | 亚洲无线在线观看| 免费在线观看成人毛片| 亚洲国产欧美网| 午夜免费成人在线视频| 国内毛片毛片毛片毛片毛片| 日本黄色片子视频| 一个人看视频在线观看www免费 | 别揉我奶头~嗯~啊~动态视频| 99久久精品一区二区三区| 午夜福利高清视频| 色综合站精品国产| 桃色一区二区三区在线观看| 特级一级黄色大片| 3wmmmm亚洲av在线观看| 亚洲黑人精品在线| 九色国产91popny在线| 亚洲国产高清在线一区二区三| 日本精品一区二区三区蜜桃| 欧美日韩一级在线毛片| 天堂av国产一区二区熟女人妻| 中文字幕av成人在线电影| 亚洲无线在线观看| 桃色一区二区三区在线观看| 香蕉久久夜色| 国产精品香港三级国产av潘金莲| 十八禁人妻一区二区| 国产精品电影一区二区三区| 嫩草影视91久久| 99国产精品一区二区蜜桃av| 免费在线观看日本一区| 国产真实伦视频高清在线观看 | 99久久精品一区二区三区| 999久久久精品免费观看国产| 淫秽高清视频在线观看| 日本三级黄在线观看| 国产午夜福利久久久久久| 90打野战视频偷拍视频| 久久久久精品国产欧美久久久| 一进一出抽搐gif免费好疼| 欧美性猛交╳xxx乱大交人| 校园春色视频在线观看| 在线a可以看的网站| 国产国拍精品亚洲av在线观看 | 久久久久久久亚洲中文字幕 | 色哟哟哟哟哟哟| 午夜福利免费观看在线| 人妻久久中文字幕网| 成人特级黄色片久久久久久久| 九色成人免费人妻av| 成人鲁丝片一二三区免费| 国语自产精品视频在线第100页| 国产精品,欧美在线| 亚洲人与动物交配视频| 搞女人的毛片| 99久久精品热视频| 国产麻豆成人av免费视频| 又黄又爽又免费观看的视频| 国产国拍精品亚洲av在线观看 | 国产伦精品一区二区三区四那| 精品久久久久久,| 男人的好看免费观看在线视频| 天堂av国产一区二区熟女人妻| 国产三级在线视频| 亚洲精品影视一区二区三区av| 看片在线看免费视频| 婷婷六月久久综合丁香| 一本一本综合久久| 午夜免费男女啪啪视频观看 | 亚洲精品在线观看二区| 天美传媒精品一区二区| 精品久久久久久久久久久久久| 精品一区二区三区人妻视频| 99国产综合亚洲精品| www.熟女人妻精品国产| 成人三级黄色视频| 成年人黄色毛片网站| 日韩国内少妇激情av| 日本精品一区二区三区蜜桃| 亚洲欧美精品综合久久99| 国产午夜福利久久久久久| 听说在线观看完整版免费高清| 99久久无色码亚洲精品果冻| 国产亚洲精品久久久com| 国产亚洲精品综合一区在线观看| 午夜a级毛片| 久久久国产成人免费| 成人高潮视频无遮挡免费网站| 哪里可以看免费的av片| 91麻豆av在线| 色哟哟哟哟哟哟| 日韩av在线大香蕉| 国产伦精品一区二区三区四那| 男人舔女人下体高潮全视频| 可以在线观看的亚洲视频| 亚洲av五月六月丁香网| 床上黄色一级片| 法律面前人人平等表现在哪些方面| 亚洲精品亚洲一区二区| 男女视频在线观看网站免费| 免费一级毛片在线播放高清视频| 日韩欧美一区二区三区在线观看| 熟女人妻精品中文字幕| 国产黄片美女视频| 久久久精品欧美日韩精品| 亚洲天堂国产精品一区在线| 欧美日本亚洲视频在线播放| 男人的好看免费观看在线视频| 岛国在线观看网站| 国产一区二区在线av高清观看| 亚洲avbb在线观看| 一进一出抽搐gif免费好疼| 久久久久久久精品吃奶| 国产探花在线观看一区二区| 免费av毛片视频| 亚洲欧美精品综合久久99| 免费看十八禁软件| 搡老熟女国产l中国老女人| 久久亚洲精品不卡| 精品无人区乱码1区二区| 国产午夜福利久久久久久| 最近最新免费中文字幕在线| 欧美乱码精品一区二区三区| 亚洲国产欧洲综合997久久,| 级片在线观看| 91av网一区二区| 免费在线观看亚洲国产| 亚洲一区二区三区色噜噜| 久久亚洲真实| 9191精品国产免费久久| 亚洲18禁久久av| 精品人妻一区二区三区麻豆 | 免费人成在线观看视频色| 国产精品av视频在线免费观看| 久久久久性生活片| 亚洲成人免费电影在线观看| 欧美激情在线99| 国产亚洲精品一区二区www| 波多野结衣高清作品| 国产在视频线在精品| 久久精品夜夜夜夜夜久久蜜豆| 欧美3d第一页| 国产高清视频在线播放一区| 一区二区三区免费毛片| 黄色日韩在线| 一个人免费在线观看电影| 韩国av一区二区三区四区| 亚洲无线在线观看| 午夜视频国产福利| 国产国拍精品亚洲av在线观看 | 国产亚洲精品一区二区www| 听说在线观看完整版免费高清| 国产激情偷乱视频一区二区| 国产一区二区亚洲精品在线观看| 中文字幕久久专区| www.色视频.com| 中文在线观看免费www的网站| 美女黄网站色视频| 免费人成在线观看视频色| 九九热线精品视视频播放| 在线国产一区二区在线| 亚洲国产欧美人成| 精品国内亚洲2022精品成人| 亚洲av美国av| 国产色婷婷99| 日韩欧美三级三区| 一卡2卡三卡四卡精品乱码亚洲| 最好的美女福利视频网| 亚洲国产色片| 国产午夜福利久久久久久| 母亲3免费完整高清在线观看| 国产精品亚洲美女久久久| 久久久久久久午夜电影| 久久久久久久亚洲中文字幕 | 2021天堂中文幕一二区在线观| 午夜精品久久久久久毛片777| 精品国产超薄肉色丝袜足j| 中国美女看黄片| 一级黄色大片毛片| 69av精品久久久久久| 亚洲最大成人中文| 搡老妇女老女人老熟妇| netflix在线观看网站| 国产一区二区在线av高清观看| 两个人看的免费小视频| 在线观看一区二区三区| 国产成人欧美在线观看| 人人妻人人澡欧美一区二区| 别揉我奶头~嗯~啊~动态视频| 在线观看免费视频日本深夜| 91久久精品国产一区二区成人 | 国产高清视频在线播放一区| 一区二区三区免费毛片| 国产av麻豆久久久久久久| 国产精品av视频在线免费观看| 国产美女午夜福利| 麻豆成人av在线观看| 久久精品综合一区二区三区| 人妻久久中文字幕网| 每晚都被弄得嗷嗷叫到高潮| 亚洲七黄色美女视频| 天堂网av新在线| 亚洲自拍偷在线| 高清毛片免费观看视频网站| 国产亚洲精品综合一区在线观看| 黄色女人牲交| 欧美黄色淫秽网站| 中文亚洲av片在线观看爽| 国产一区二区三区在线臀色熟女| 在线播放国产精品三级| 此物有八面人人有两片| 99riav亚洲国产免费| 欧美+日韩+精品| 欧美黄色片欧美黄色片| 少妇的丰满在线观看| 国产黄片美女视频| 国产精品三级大全| 色av中文字幕| 在线观看av片永久免费下载| 18禁美女被吸乳视频| 国产高清三级在线| 日韩欧美国产一区二区入口| av视频在线观看入口| 精品一区二区三区av网在线观看| 国产av在哪里看| 男插女下体视频免费在线播放| 亚洲国产精品合色在线| 亚洲aⅴ乱码一区二区在线播放| 日本熟妇午夜| 最近最新中文字幕大全免费视频| 无限看片的www在线观看| 18禁在线播放成人免费| 国产伦一二天堂av在线观看| 男女之事视频高清在线观看| 特级一级黄色大片| 91在线精品国自产拍蜜月 | www.999成人在线观看| 一区二区三区国产精品乱码| 一个人免费在线观看的高清视频| 国产亚洲精品av在线| 国产黄色小视频在线观看| 日本熟妇午夜|