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

    附屬構(gòu)件及橋面粗糙度對(duì)近流線型寬體箱梁氣動(dòng)靜力系數(shù)的影響

    2016-06-23 13:03:01張亮亮劉連杰劉書洋
    實(shí)驗(yàn)流體力學(xué) 2016年1期
    關(guān)鍵詞:迎角欄桿升力

    張亮亮, 吳 波,*, 楊 陽(yáng), 劉連杰, 劉書洋

    (1. 重慶大學(xué) 土木工程學(xué)院, 重慶 400045; 2. 山地城鎮(zhèn)建設(shè)安全與防災(zāi)協(xié)同創(chuàng)新中心, 重慶 400045)

    附屬構(gòu)件及橋面粗糙度對(duì)近流線型寬體箱梁氣動(dòng)靜力系數(shù)的影響

    張亮亮1,2, 吳 波1,2,*, 楊 陽(yáng)1,2, 劉連杰1, 劉書洋1

    (1. 重慶大學(xué) 土木工程學(xué)院, 重慶 400045; 2. 山地城鎮(zhèn)建設(shè)安全與防災(zāi)協(xié)同創(chuàng)新中心, 重慶 400045)

    針對(duì)高寬比為1/12的近流線型寬體箱梁進(jìn)行了一系列參數(shù)化的三分力測(cè)力試驗(yàn),研究了欄桿透風(fēng)率、檢修車軌道位置、橋面粗糙度對(duì)三分力系數(shù)的影響。運(yùn)用CFD工程軟件Fluent 15.0模擬二維風(fēng)場(chǎng),對(duì)三分力系數(shù)變化的可能機(jī)理進(jìn)行討論。結(jié)果表明,對(duì)于近流線型寬體箱梁:(1) 隨迎角由負(fù)到正加大,CD先降后升,CL與Cm總體趨于升高,但均不關(guān)于零迎角對(duì)稱或反對(duì)稱;(2) 欄桿透風(fēng)率的減小一方面加大了斷面擋風(fēng)面積,使尾流變寬,CD升高;另一方面擴(kuò)大了斷面上方負(fù)壓區(qū)范圍(正迎角工況),提高了上表面壓力系數(shù)(負(fù)迎角工況),使正、負(fù)迎角下的升力均增加,即CL絕對(duì)值升高;(3) 正迎角來流時(shí),檢修車軌道的內(nèi)移使斷面下方的負(fù)壓區(qū)縮小、壓力系數(shù)減小、尾流變窄,則CD、CL均降低;但負(fù)迎角來流時(shí),其位置的移動(dòng)未對(duì)分離區(qū)造成明顯影響,則CD與CL均未出現(xiàn)顯著波動(dòng);(4) 正迎角來流時(shí),橋面粗糙度的提高減小了上表面的分離角,使負(fù)壓區(qū)縮小、尾流變窄,則CD降低;且負(fù)壓區(qū)的縮小使得朝上的升力減小,即CL降低。負(fù)迎角來流時(shí),粗糙度的提高使上表面的阻力與壓力系數(shù)均增加,則CD升高、CL的絕對(duì)值升高;(5) 欄桿透風(fēng)率、檢修車軌道位置及橋面粗糙度對(duì)Cm的影響可以忽略不計(jì)。

    扁平箱梁;氣動(dòng)靜力系數(shù);節(jié)段模型試驗(yàn);橋面粗糙度;欄桿透風(fēng)率;檢修車軌道位置;繞流場(chǎng)特性

    0 引 言

    自然風(fēng)對(duì)橋梁的作用分為定常力與非定常力,定常力是靜風(fēng)荷載、抖振響應(yīng)、馳振穩(wěn)定性分析等抗風(fēng)研究的基礎(chǔ),可通過阻力系數(shù)、升力系數(shù)及力矩系數(shù)來表示[1-3]。目前,獲取靜力三分力系數(shù)的常用方法是節(jié)段模型測(cè)力實(shí)驗(yàn)。諸多研究表明[4-7],斷面形狀、附屬結(jié)構(gòu)物對(duì)靜力三分力的影響不可忽視。馬凱[6]指出,近流線型斷面的三分力系數(shù)與桁架橋、疊合梁橋明顯不同;高亮[7]研究了欄桿、檢修車軌道對(duì)帶中央開槽的分離式箱梁(高寬比1/9)三分力系數(shù)的影響;但關(guān)于附屬結(jié)構(gòu)物對(duì)整體式寬體扁平箱梁三分力系數(shù)的影響,目前鮮有專門報(bào)道。關(guān)于橋面粗糙度對(duì)斷面三分力系數(shù)影響的報(bào)道則更少。針對(duì)上述問題,以高寬比1/12的整體式寬體箱梁為對(duì)象,進(jìn)行了1∶60的節(jié)段模型風(fēng)洞實(shí)驗(yàn),系統(tǒng)地探討了迎角、欄桿透風(fēng)率、檢修車軌道位置以及橋面粗糙度對(duì)流線型斷面三分力系數(shù)的影響。

    計(jì)算流體動(dòng)力學(xué)(CFD)是數(shù)值模擬的理論基礎(chǔ),具有良好的可視化功能,可將氣流流經(jīng)橋梁斷面時(shí)的流場(chǎng)分布特性直觀地呈現(xiàn)出來[8-9]。一些學(xué)者[10-11]利用CFD識(shí)別橋梁斷面的三分力系數(shù),取得了較高的精度。文中對(duì)主梁斷面進(jìn)行了CFD二維求解,對(duì)比了斷面在不同工況下的壓強(qiáng)、流線分布特征,從微觀層面探討三分力系數(shù)隨迎角、欄桿透風(fēng)率、檢修車軌道位置以及橋面粗糙度變化的可能機(jī)理。

    1 工程概況

    擬建的重慶寸灘長(zhǎng)江大橋是重慶機(jī)場(chǎng)專用快速路工程南段的重要組成部分,包括跨江主橋和南北引橋,總長(zhǎng)1.6km,跨徑250m+880m+250m。主橋?yàn)殡p塔單跨懸索橋,主跨880m,矢跨比1/8.8。主梁為寬體扁平鋼箱梁,寬42m,中心處梁高3.5m,高寬比1/12,設(shè)2%雙向排水坡。圖1為主梁標(biāo)準(zhǔn)斷面圖。

    圖1 主梁標(biāo)準(zhǔn)斷面(單位:cm)

    2 節(jié)段模型風(fēng)洞實(shí)驗(yàn)

    2.1 實(shí)驗(yàn)方案

    主梁節(jié)段模型縮尺比1∶60,采用環(huán)氧樹脂板與優(yōu)質(zhì)木材制作,長(zhǎng)L=2.1m,寬B=0.70m,高D=0.058m。試驗(yàn)在西南交通大學(xué)工業(yè)風(fēng)洞(XNJD-1)第二試驗(yàn)段中進(jìn)行,模型由三分力專用側(cè)壁支撐固定于三分量應(yīng)變式天平,來流為均勻流,風(fēng)速U=10m/s。通過計(jì)算機(jī)控制模型姿態(tài)角來設(shè)置來流迎角,α=-12°~12°,Δα=1°;使用數(shù)據(jù)采集儀測(cè)量靜風(fēng)荷載下的阻力(FD)、升力(FL)及俯仰力矩(MZ),然后按式(1)~(3)計(jì)算風(fēng)軸下的三分力系數(shù);圖2為成橋狀態(tài)靜力三分力實(shí)驗(yàn)。

    圖2 靜力三分力系數(shù)節(jié)段模型測(cè)力試驗(yàn)

    試驗(yàn)中通過在模型上表面粘貼輕、薄的砂紙來模擬橋面粗糙度,所用砂紙的厚度僅為模型高度的0.01~0.02倍,確保砂紙平整貼合于模型表面,避免對(duì)主梁外形產(chǎn)生影響。

    2.2 實(shí)驗(yàn)工況

    除裸梁工況外,實(shí)驗(yàn)還設(shè)置了無欄桿(有檢修車軌道)及3種不同透風(fēng)率的欄桿(59.8%、45.8%、35.6%)、3處不同位置的檢修車軌道(分別距底板邊緣0、2、4cm)及3種橋面粗糙度(無砂紙、120目砂紙、80目砂紙),每類工況均測(cè)量了來流迎角α=-12°~12°下的三分力系數(shù)。

    2.3 靜力三分力系數(shù)

    主梁斷面在風(fēng)軸坐標(biāo)系下的靜力三分力系數(shù)計(jì)算式為[12]:

    阻力系數(shù):

    升力系數(shù):

    力矩系數(shù):

    式中,U表示來流風(fēng)速,ρ為空氣密度,取1.225kg/m3。FD、FL、MZ分別表示橫橋向阻力、豎向升力及俯仰力矩。B為參考寬度,取箱梁頂板寬度0.70m。D為參考高度,取箱梁中心高度0.058m。L為節(jié)段長(zhǎng)度,取2.1m。

    3 CFD數(shù)值模擬

    數(shù)值模擬軟件為Ansys Fluent 15.0,由于主梁外形沿主跨方向不變,處理為二維,流動(dòng)為定常不可壓,靜力按定常處理,不考慮溫度的影響。對(duì)應(yīng)的數(shù)學(xué)模型為雷諾時(shí)均N-S方程,湍流模型為標(biāo)準(zhǔn)k-ε模型,近壁面以標(biāo)準(zhǔn)壁面函數(shù)處理,壓力速度耦合選用

    SIMPLE算法,離散格式為二階迎風(fēng)插值[11,13]。

    在前處理軟件GAMBIT中形成計(jì)算網(wǎng)格,數(shù)值模型的縮尺比也為1:60,與節(jié)段模型實(shí)驗(yàn)相同。為減小邊界對(duì)斷面流場(chǎng)的影響,計(jì)算區(qū)域設(shè)置得足夠大[14-15]:上下邊界到斷面中心的距離為1.0m(>15D),入口距斷面中心4.2m(6B),出口距斷面中心8.4m(12B);入口為速度條件,出口為壓力出口,主梁及其附屬物外壁均為無滑移的壁面。

    斷面附近第1層網(wǎng)格高度0.5mm,附屬構(gòu)件附近的第1層網(wǎng)格高度0.1mm,使壁面Y+≈30[16]。在流場(chǎng)變化劇烈的地方劃分較密的網(wǎng)格,然后逐漸過渡,外部網(wǎng)格劃分較疏,以使網(wǎng)格更能適應(yīng)流場(chǎng)的變化,最大網(wǎng)格尺寸為0.035m。

    Fluent中以Roughness height表示粗糙度厚度Ks。對(duì)于砂粒平均分布的砂紙,Ks取粒徑平均高度。經(jīng)游標(biāo)卡尺測(cè)量,實(shí)驗(yàn)中所用的120目、80目砂紙對(duì)應(yīng)的Ks分別為0.97mm與1.28mm。將包括風(fēng)嘴在內(nèi)的模型上表面單獨(dú)劃分成一個(gè)壁面邊界,求解時(shí),設(shè)置相應(yīng)的Roughness height值即可模擬橋面粗糙。為更好分辨橋面的粗糙度顆粒,在斷面上表面劃分邊界層網(wǎng)格,第1層網(wǎng)格高度0.05mm。

    4 結(jié) 果

    4.1 CFD計(jì)算值與實(shí)驗(yàn)結(jié)果的對(duì)比

    表1列出了一些典型工況三分力系數(shù)的CFD計(jì)算值與實(shí)驗(yàn)結(jié)果的對(duì)比,可見:CFD可準(zhǔn)確反映CD、CL、Cm隨迎角、欄桿透風(fēng)率、檢修車軌道位置及橋面粗糙度的變化趨勢(shì)。其中,CD的模擬精度最高,CL次之,Cm的模擬精度稍差,但由于其值很小,接近于0,誤差尚屬于可接受的范圍[10],這與譚紅霞[11]利用標(biāo)準(zhǔn)k-ε模型取得的模擬精度相當(dāng)。

    對(duì)于裸梁工況,涉及到模擬橋面粗糙度的工況k~n的精度尚可,但稍差于不模擬粗糙度的工況i、j,可能與砂紙顆粒被近似處理為平均高度的粗糙度條帶有關(guān)。

    表1 代表性工況三分力系數(shù)的CFD計(jì)算值與實(shí)驗(yàn)值的對(duì)比Table 1 Comparison between CFD and experimental values under representative working cases

    4.2 結(jié)果分析

    4.2.1 迎角的影響

    圖4為本次試驗(yàn)條件下,成橋態(tài)10m/s風(fēng)速下的風(fēng)軸三分力系數(shù)隨迎角α的變化曲線,可見:(1)隨α加大,阻力系數(shù)CD先降后升,即小迎角范圍內(nèi)較小,正、負(fù)大迎角范圍內(nèi)更大,主要是由于大迎角下斷面的擋風(fēng)面積更大。但由于斷面上、下部的不對(duì)稱及風(fēng)嘴、橋面附屬物等的作用,CD-α曲線并不是關(guān)于0°迎角對(duì)稱,最小值位于-3°附近;(2)負(fù)迎角來流時(shí),分離區(qū)位于斷面以下,呈負(fù)壓,斷面以上呈正壓,上、下壓力差為正值,則升力FL朝下,即升力系數(shù)CL為負(fù)值;隨α趨于0°迎角,分離區(qū)縮小,上、下壓力差減小,F(xiàn)L減小,CL的絕對(duì)值減小,即CL-α曲線呈正斜率;正迎角來流時(shí),分離區(qū)位于斷面以上,上、下壓力差為負(fù)值,升力FL朝上,即升力系數(shù)CL為正值;隨α加大,分離區(qū)加大,CL逐漸升高;(3)隨α加大,力矩系數(shù)Cm趨于升高,但波動(dòng)范圍小于CL。

    圖4 三分力系數(shù)隨迎角的變化曲線(成橋態(tài),U=10m/s)

    4.2.2 欄桿透風(fēng)率的影響

    在研究欄桿透風(fēng)率的影響時(shí),4類工況的透風(fēng)率(ventilation rate)分別為100%(無欄桿)、59.8%、45.8%及35.6%。試驗(yàn)條件下,欄桿透風(fēng)率對(duì)三分力系數(shù)的影響見圖5:

    (1) 隨透風(fēng)率減小,CD逐漸升高,尤以小迎角范圍的趨勢(shì)更明顯;(2) 在正迎角范圍內(nèi),CL隨透風(fēng)率減小而升高;在負(fù)迎角范圍內(nèi),CL隨透風(fēng)率減小而降低(負(fù)值,絕對(duì)值增加);(3)Cm未隨透風(fēng)率的變化而發(fā)生明顯波動(dòng):最大幅度低于0.02。

    選取5°為代表來研究三分力系數(shù)隨欄桿透風(fēng)率變化的微觀機(jī)理。圖6為不同透風(fēng)率工況下(工況c~f),斷面上、下表面壓力系數(shù)隨無量綱距離x/B變化曲線。其中,壓力系數(shù)的計(jì)算式為:

    圖5 欄桿透風(fēng)率對(duì)三分力系數(shù)的影響(U=10m/s)

    圖6 欄桿透風(fēng)率對(duì)主梁表面壓力系數(shù)的影響(α=5°)

    (4)

    其中,pi為測(cè)點(diǎn)靜壓,p為參考點(diǎn)靜壓。

    圖7為各工況的流線分布,結(jié)合圖6,將三分力系數(shù)隨欄桿透風(fēng)率變化的原因歸結(jié)如下:

    (1) 與有欄桿的工況c、d、e相比,無欄桿的工況f擋風(fēng)面積更小,尾流寬度更小,CD明顯更低;由于缺少欄桿的阻礙,工況f的上表面壓力系數(shù)低于(負(fù)值,絕對(duì)值大于)、下表面壓力系數(shù)高于工況c、d、e,可獲得更大的升力(朝上,正值),即CL更大;隨透風(fēng)率逐漸減小,工況c、d、e上氣流再附點(diǎn)向下游移動(dòng),加之擋風(fēng)面積的增加使上表面的分離角逐漸加大,負(fù)壓分離區(qū)加高,主旋渦加大(圖7),使尾流變寬,因而CD逐漸升高;小迎角來流作用時(shí),氣流近乎垂直作用于欄桿,則透風(fēng)率降低所增加的擋風(fēng)效應(yīng)大于大迎角工況,因此CD-α曲線在小迎角范圍內(nèi)的升高趨勢(shì)更顯著。

    (2) 正迎角來流時(shí),隨欄桿透風(fēng)率減小,上表面分離區(qū)加高,負(fù)壓范圍更大,則回流作用更強(qiáng),風(fēng)速更大;由伯努利方程,風(fēng)速升高引起靜壓減小,即上表面壓力系數(shù)隨欄桿透風(fēng)率減小而減小(圖6),但由于其呈負(fù)值,實(shí)際上是絕對(duì)值增大,因此工況c、d、e的升力(朝上)逐漸增加,即CL-α曲線在正迎角范圍內(nèi)隨透風(fēng)率減小而升高;

    (3) 負(fù)迎角來流時(shí),分離區(qū)形成于斷面下方,欄桿透風(fēng)率的改變并不會(huì)對(duì)下表面壓力系數(shù)造成較大影響。而斷面以上處于正壓區(qū),透風(fēng)率的減小加強(qiáng)了對(duì)上表面氣流的阻擋,使風(fēng)速減??;由伯努利方程,壓力系數(shù)升高(圖6),則升力(朝下)更大,即CL降低(負(fù)值,絕對(duì)值增加)。

    圖7 欄桿透風(fēng)率對(duì)主梁上表面附近流線分布的影響(α=5°)

    Fig.7 Impact of ventilation rate on streamline distributions near the upper surface

    4.2.3 檢修車軌道位置的影響

    在研究檢修車軌道位置的影響時(shí),3類工況的檢修車軌道分別位于底板最外側(cè)(DT=0)、內(nèi)移2cm(DT=2cm)、內(nèi)移4cm(DT=4cm),欄桿透風(fēng)率均為59.8%。試驗(yàn)條件下,檢修車軌道位置對(duì)三分力系數(shù)的影響如圖8所示。

    圖8 檢修車軌道位置對(duì)三分力系數(shù)的影響(U=10m/s)

    (1) 隨軌道內(nèi)移,CD在正迎角范圍內(nèi)逐漸降低,但在負(fù)迎角范圍內(nèi)未出現(xiàn)顯著波動(dòng);(2)隨軌道內(nèi)移,CL在正迎角范圍內(nèi)逐漸降低,在負(fù)迎角范圍內(nèi)未發(fā)生明顯變化;(3)Cm未隨檢修車軌道的移動(dòng)而發(fā)生明顯波動(dòng)。

    在研究檢修車軌道位置的移動(dòng)對(duì)三分力系數(shù)造成影響的微觀機(jī)理時(shí),同樣選取α=5°為代表工況。圖9為不同檢修車軌道位置(工況c、g、h)的斷面上、下表面壓力系數(shù)曲線,圖10為其流線分布。將正迎角來流下,三分力系數(shù)隨檢修車軌道位置移動(dòng)而發(fā)生變化的原因歸結(jié)如下:

    (1) 正迎角來流作用下,隨軌道內(nèi)移,斷面下表面的分離點(diǎn)內(nèi)移(如圖9),即負(fù)壓分離區(qū)縮小,且氣流更緊貼斷面而行,使尾流變窄(如圖10),則CD逐漸降低;軌道內(nèi)移在縮小下方分離區(qū)的同時(shí),還減小了下表面壓力系數(shù),則斷面上、下壓力差減小,升力(朝上)減小,即CL降低;

    圖9 檢修車軌道位置對(duì)主梁表面壓力系數(shù)的影響(α=5°)

    Fig.9 Impact of vehicle track location on the surface pressure coefficient

    圖10 檢修車軌道位置對(duì)主梁下表面附近流線分布的影響(α=5°)

    Fig.10 Impact of overhaul-doll rail location on the streamline distributions near the lower surface

    (2) 負(fù)迎角來流時(shí),氣流在風(fēng)嘴與斜腹板附近即完成了分離,位于斷面以下的分離區(qū)將檢修車軌道全部包圍,軌道位置的移動(dòng)并不會(huì)對(duì)分離區(qū)和尾流產(chǎn)生較大的影響,因此CD與CL均未出現(xiàn)顯著波動(dòng)。

    4.2.4 橋面粗糙度的影響

    實(shí)驗(yàn)?zāi)P捅砻嬗森h(huán)氧樹脂板制作,其粗糙度可忽略不計(jì);此外,在模型上表面粘貼了120目、80目2種規(guī)格的砂紙來研究橋面粗糙度的影響。3類工況均在裸梁的基礎(chǔ)上進(jìn)行,以突出粗糙度的影響。試驗(yàn)條件下,橋面粗糙度對(duì)裸梁三分力系數(shù)的影響如圖11所示。

    圖11 橋面粗糙度對(duì)三分力系數(shù)的影響(裸梁,U=10m/s)

    (1) 隨橋面粗糙度加大,CD在正迎角范圍內(nèi)逐漸降低,在負(fù)迎角范圍內(nèi)逐漸升高;(2)隨橋面粗糙度加大,CL在正負(fù)迎角范圍內(nèi)均趨于降低;(3)橋面粗糙度的變化未對(duì)Cm造成明顯影響。

    在討論正迎角來流時(shí)橋面粗糙度對(duì)三分力系數(shù)影響的微觀機(jī)理時(shí),仍然選取α=5°為代表。圖12為不同橋面粗糙度工況(i、k、m)的上表面附近流線分布,將正迎角來流下三分力系數(shù)隨粗糙度變化的原因歸結(jié)如下。

    圖12 正迎角下橋面粗糙度對(duì)上表面流線的影響(裸梁,α=5°)

    Fig.12 Impact of deck roughness on the streamline distributions near the upper surface at a positive attack angle

    正迎角來流作用下,斷面上方處于負(fù)壓區(qū),粗糙度顆粒對(duì)貼近表面的氣流具有一定的吸附作用,隨粗糙度的提高,吸附作用加強(qiáng),使風(fēng)嘴處的分離角縮小,分離區(qū)變窄,旋渦縮小(如圖12),則尾流變窄,因此CD逐漸降低;同時(shí),分離區(qū)在豎直方向的變窄即上方負(fù)壓區(qū)的縮小,因此斷面所受升力(朝上)降低,即CL降低。

    負(fù)迎角下,橋面粗糙度對(duì)主梁壓力系數(shù)的影響見圖13。將負(fù)迎角來流下三分力系數(shù)隨粗糙度變化的原因歸結(jié)如下。

    負(fù)迎角來流作用時(shí),上表面處于正壓區(qū),粗糙度顆粒對(duì)氣流起阻擋作用。CFD的計(jì)算結(jié)果可定性反映這種阻擋作用:由頂板所產(chǎn)生的CD分量從工況j(Roughness height=0)的0.265增加至工況n(Roughness height=1.28mm)的0.366,增量0.101。而工況j、n的CD差值為0.105,即上表面阻力的增加是CD隨粗糙度提高而升高的主要原因;此外,粗糙度顆粒阻擋作用的加強(qiáng)使得上表面風(fēng)速略微減小,由伯努利方程可知,壓力系數(shù)升高(如圖13),則上、下壓力差更大,即升力(朝下)更大,CL降低(負(fù)值,絕對(duì)值升高)。

    圖13 負(fù)迎角下橋面粗糙度對(duì)主梁壓力系數(shù)的影響(裸梁,α=-5°)

    Fig.13 Impact of deck roughness on the surface pressure coefficient at a negative attack angle

    5 結(jié) 論

    橋梁斷面三分力系數(shù)隨迎角、欄桿透風(fēng)率、檢修車軌道位置及橋面粗糙度變化的微觀機(jī)理在于繞流場(chǎng)的變化。對(duì)于近流線型的扁平箱梁:

    (1) 隨迎角由負(fù)到正加大,阻力系數(shù)CD先降后升,升力系數(shù)CL與力矩系數(shù)Cm總體趨于升高,但均不關(guān)于0°迎角對(duì)稱或反對(duì)稱;

    (2) 隨欄桿透風(fēng)率逐漸減小,上表面氣流再附點(diǎn)朝下游移動(dòng),且斷面擋風(fēng)面積逐漸增加,則尾流變寬,CD升高;正迎角來流時(shí),欄桿透風(fēng)率的減小擴(kuò)大了斷面上方負(fù)壓區(qū)的范圍,則斷面所受朝上的升力更大,即CL升高;負(fù)迎角來流時(shí),欄桿透風(fēng)率的減小使得上表面風(fēng)速減小,壓力系數(shù)升高,則朝下的升力更大,即負(fù)值CL降低(絕對(duì)值升高);

    (3) 正迎角來流時(shí),隨檢修車軌道內(nèi)移,下表面分離點(diǎn)內(nèi)移,即負(fù)壓分離區(qū)縮小,且氣流更緊貼斷面而行,則尾流變窄,CD逐漸降低;分離、再附點(diǎn)位置的變化也使得斷面下方負(fù)壓區(qū)范圍縮小,壓力系數(shù)減小,則朝上的升力減小,CL降低;但由于負(fù)迎角來流時(shí),檢修車軌道位置的移動(dòng)未對(duì)分離區(qū)造成較大影響,因此CD與CL均未出現(xiàn)顯著波動(dòng);

    (4) 正迎角來流時(shí),粗糙度的提高使上表面分離角縮小,則尾流變窄,CD降低;分離角的縮小也使得斷面上方負(fù)壓區(qū)范圍縮小,則朝上的升力降低,即CL降低。而負(fù)迎角來流時(shí),粗糙度的提高使得頂板所受阻力逐漸增加,即CD逐漸升高;此外,粗糙度的提高還增加了上表面的壓力系數(shù),使朝下的升力增加,即負(fù)值CL降低,絕對(duì)值升高;

    (5) 本實(shí)驗(yàn)中各工況的Cm均很小(-0.18~0.16),欄桿透風(fēng)率、檢修車軌道位置及橋面粗糙度對(duì)其的影響均小于0.02的變化幅值,基本上可忽略不計(jì)。

    [1] 陳政清. 橋梁風(fēng)工程[M]. 北京: 人民交通出版社, 2005.

    Cheng Zhengqing. Wind engineering of bridge[M]. Beijing: China Communications Press, 2005.

    [2] 王騎, 廖海黎, 李明水, 等. 大振幅下薄翼和流線型箱梁的氣動(dòng)遲滯研究[J]. 實(shí)驗(yàn)流體力學(xué), 2013, 27(1): 32-37, 45.

    Wang Qi, Liao Haili, Li Mingshui, et al. Aerodynamic hysteresis of thin airfoil and streamline box girder under large amplitude oscillation[J]. Journal of Experiments in Fluid Mechanics, 2013, 27(1): 32-37, 45.

    [3] 喻梅, 廖海黎, 李明水, 等. 大跨度橋梁斜風(fēng)作用下抖振響應(yīng)現(xiàn)場(chǎng)實(shí)測(cè)及風(fēng)洞試驗(yàn)研究[J]. 實(shí)驗(yàn)流體力學(xué), 2013, 27(3): 51-55, 76.

    Yu Mei, Liao Haili, Li Mingshui, et al. Field measurement and wind tunnel test of buffeting response of long-span bridge under skew wind[J]. Journal of Experiments in Fluid Mechanics, 2013, 27(3): 51-55, 76.

    [4] 李永樂, 周昱, 葛世平, 等. 主梁斷面形狀對(duì)車-橋系統(tǒng)氣動(dòng)特性影響的風(fēng)洞試驗(yàn)研究[J]. 土木工程學(xué)報(bào), 2012, 45(7): 127-133.

    Li Yongle, Zhou Yu, Ge Shiping, et al. Wind tunnel test study of the effect of bridge deck cross section shape on aerodynamic characteristics of vehicle-bridge system[J]. China Civil Engineering Journal, 2012, 45(7): 127-133.

    [5] 王向陽(yáng), 宋宏敏, 熊銳. 橋梁斷面形式對(duì)三分力系數(shù)影響的數(shù)值研究[J]. 公路, 2012, (12): 39-42.

    Wang Xiangyang, Song Hongmin, Xiong Rui. Numerical study on the influence of cross section of three component coefficients of the bridge[J]. Highway, 2012, (12): 39-42.

    [6] 馬凱. 典型橋梁斷面靜力三分力系數(shù)研究[C]//第二十一屆全國(guó)橋梁學(xué)術(shù)會(huì)議論文集(下冊(cè)). 大連: 中國(guó)工程院土木水利與建筑工程學(xué)部, 2014: 914-923.

    Ma Kai. Study on the static three component coefficient of typical bridge section[C]//21th National Bridge Academic Conference(II). Dalian: Departmen of Civil and Construction Engineering of the Chinese Academy of Engineering, 2014: 914-923.

    [7] 高亮. 橋梁斷面三分力系數(shù)試驗(yàn)研究[D]. 西安: 長(zhǎng)安大學(xué), 2010.

    Gao Liang. Investigation into three-component force coefficients of bridge deek section[D]. Xi’an: Chang’an University, 2010.

    [8] 李永樂, 汪斌, 黃林, 等. 平板氣動(dòng)力的CFD模擬及參數(shù)研究[J]. 工程力學(xué), 2009, 26(03): 207-211.

    Li Yongle, Wang Bin, Huang Lin, et al. CFD simulation and parameter study on aerodynamic force of flat plate[J]. Engineering Mechanics, 2009, 26(3): 207-211.

    [9] 《中國(guó)公路學(xué)報(bào)》編輯部. 中國(guó)橋梁工程學(xué)術(shù)研究綜述·2014[J]. 中國(guó)公路學(xué)報(bào), 2014, 27(5): 1-96.

    Editorial Department of China Journal of Highway and Transport. Review on China’s bridge engineering research: 2014[J]. China Journal of Highway and Transport, 2014, 27(5): 1-96.

    [10] 劉鑰, 陳政清, 張志田. 箱梁斷面靜風(fēng)力系數(shù)的CFD 數(shù)值模擬[J]. 振動(dòng)與沖擊, 2010, 29(1): 133-137, 242.

    Liu Yue, Chen Zhengqing, Zhang Zhitian. CFD numerical simulation of aerostatic coefficients of a box girder section[J]. Journal of Vibration and Shock, 2010, 29(1): 133-137, 242.

    [11] 譚紅霞, 陳政清. CFD在橋梁斷面靜力三分力系數(shù)計(jì)算中的應(yīng)用[J]. 工程力學(xué), 2009, 26(11): 68-72.

    Tan Hongxia, Chen Zhengqin. Application of CFD in calculating static coefficients of bridge section[J]. Engineering Mechanics, 2009, 26(11): 68-72.

    [12] JTG/T D60-01-2004, 公路橋梁抗風(fēng)設(shè)計(jì)規(guī)范[S]. 北京: 人民交通出版社, 2004.

    JTG/T D60-01-2004, Wind resistant code for design of highway bridge[S]. Beijing: China Communications Press, 2004.

    [13] 劉鑰. 連續(xù)剛構(gòu)橋梁的氣動(dòng)性能數(shù)值模擬研究[D]. 長(zhǎng)沙: 湖南大學(xué), 2009.

    Liu Yue. The research on numerical simulation of ae rostatic performance of continuous rigid frame bridge[D]. Changsha: Hunan University, 2009.

    [14] 張學(xué)敏. 橋梁風(fēng)效應(yīng)數(shù)值模擬研究[D]. 長(zhǎng)沙: 湖南大學(xué), 2014.

    Zhang Xuemin. Numerical simulation of wind effects on bridges[D]. Changsha: Hunan University, 2014.

    [15] 曾鍇, 汪叢軍, 黃本才, 等. 計(jì)算風(fēng)工程中幾個(gè)關(guān)鍵影響因素的分析與建議[J]. 空氣動(dòng)力學(xué)學(xué)報(bào), 2007, 25(4): 504-508.

    Zeng Kai, Wang Congjun, Huang Bencai, et al. Analysis and suggestions of several key factors in CFD[J]. Acta Aerodynamica Sinica, 2007, 25(4): 504-508.

    [16] 王福軍. 計(jì)算流體動(dòng)力學(xué)分析[M]. 北京: 清華大學(xué)出版社, 2004.

    Wang Fujun. Analysis of computational fluid dynamics[M]. Beijing: Tsinghua University Press, 2004.

    (編輯:張巧蕓)

    Effects of subsidiary members and deck surface roughness on the aerodynamic coefficients of static forces on a flat box girder

    Zhang Liangliang1,2, Wu Bo1,2,*, Yang Yang1,2, Liu Lianjie1, Liu Shuyang1

    (1. School of Civil Engineering, Chongqing University, Chongqing 400045, China; 2. Key Laboratory of New Technology for Construction of Cities in Mountain Area, Chongqing 400045, China)

    Taking Cuntan Yangtze Bridge, which has a height-width ratio of 1/12, as an example, a series of parametric wind tunnel tests were conducted to study the effects of ventilation rate of railing, location of overhaul-dolly rails and deck surface-roughness on the bridge’s aerodynamic coefficients of static forces. The possible variation mechanism of the aerodynamic coefficients were also discussed via comparing the flow patterns and characteristics simulated by CFD software Fluent 15.0. The following results are obtained for a streamlined flat box girder: (1) As attack angle increases from negative to positive, the drag coefficient (CD) decreases firstly and then increases, while both the lift coefficient (CL) and moment coefficient (Cm) have an overall tendency of increasing. (2) The decrease of ventilation rate, on the one hand, results in an increase of the wind acting area, which leads to an expansion of the wake zone and a consequent increase ofCD;on the other hand, it extends the negative pressure zone in the case of positive attack angles, or increases the pressure coefficients on the upper surface in the case of negative attack angles, so that the magnitude ofCLincreases in both cases of positive and negative attack angles. (3) In the case of positive attack angles, when the overhaul-dolly rails are moved inward, the blunt degree of the cross section is reduced, narrowing both the negative pressure zone and the wake zone, and decreasing bothCDandCLconsequently. However, in the case of negative attack angles,the changes of separation and reattachment processes are not significant when the overhaul-dolly rails are moved, and therefore neitherCDnorCLhas a significant fluctuation. (4)In the case of positive attack angles,the increase of deck surface roughness narrows the separation zone over the upper face and the wake zone,leading to the decrease of bothCDandCL. In the case of negative attack angles,the increase of the deck surface roughness leads to an increase of both the resistance and pressure coefficients and the consequent increase of the absolute values of bothCDandCL. (5)The changes ofCmwith the ventilation rate of railing,the location of overhaul-dolly rails and the deck surface-roughness are generally negligible.

    flat box girder;aerodynamic coefficients of static force;wind tunnel test;railing ventilation rate;location of overhaul-dolly rail;deck surface roughness;flow characteristic

    1672-9897(2016)01-0074-07

    10.11729/syltlx20150059

    2015-05-01;

    2015-09-24

    國(guó)家自然科學(xué)基金(51578098);重慶市科委項(xiàng)目(20130844)

    ZhangLL,WuB,YangY.Effectsofsubsidiarymembersanddecksurfaceroughnessontheaerodynamiccoefficientsofstaticforcesonaflatboxgirder.JournalofExperimentsinFluidMechanics, 2016, 30(1): 74-80. 張亮亮, 吳波, 楊陽(yáng). 附屬構(gòu)件及橋面粗糙度對(duì)近流線型寬體箱梁氣動(dòng)靜力系數(shù)的影響. 實(shí)驗(yàn)流體力學(xué), 2016, 30(1): 74-80.

    U441+.2

    A

    張亮亮(1956-),男, 湖南益陽(yáng)人,教授,博士生導(dǎo)師。研究方向:橋梁力學(xué)性能分析、橋梁抗風(fēng)。通信地址:重慶市沙坪壩區(qū)沙坪壩北街83號(hào)重慶大學(xué)土木工程學(xué)院(400045)。E-mail: zll200510@126.com

    *通信作者 E-mail: wuswjtu@yeah.net

    猜你喜歡
    迎角欄桿升力
    Slowly and Quickly慢慢地,快快地
    跨過童年的欄桿
    高速列車車頂–升力翼組合體氣動(dòng)特性
    連續(xù)變迎角試驗(yàn)數(shù)據(jù)自適應(yīng)分段擬合濾波方法
    無人機(jī)升力測(cè)試裝置設(shè)計(jì)及誤差因素分析
    基于自適應(yīng)偽譜法的升力式飛行器火星進(jìn)入段快速軌跡優(yōu)化
    高速公路手動(dòng)欄桿智能化改進(jìn)方案設(shè)計(jì)
    升力式再入飛行器體襟翼姿態(tài)控制方法
    失速保護(hù)系統(tǒng)迎角零向跳變研究
    科技傳播(2014年4期)2014-12-02 01:59:42
    做好小事不容易
    美女视频免费永久观看网站| 一区在线观看完整版| 2021少妇久久久久久久久久久| 亚洲精品视频女| 欧美黑人欧美精品刺激| 国产伦理片在线播放av一区| 如何舔出高潮| 国产精品嫩草影院av在线观看| 9191精品国产免费久久| 91精品伊人久久大香线蕉| 亚洲精品国产色婷婷电影| 国产片内射在线| √禁漫天堂资源中文www| www.精华液| av一本久久久久| 国产男人的电影天堂91| 亚洲国产欧美在线一区| 亚洲精品国产一区二区精华液| 成年女人毛片免费观看观看9 | 麻豆乱淫一区二区| 十八禁人妻一区二区| 日韩大码丰满熟妇| 国产国语露脸激情在线看| 色综合欧美亚洲国产小说| 免费av中文字幕在线| 又粗又硬又长又爽又黄的视频| 大片电影免费在线观看免费| 亚洲精品乱久久久久久| 国产男女内射视频| 国产精品女同一区二区软件| 黑人欧美特级aaaaaa片| 成人漫画全彩无遮挡| 一级毛片 在线播放| 欧美精品一区二区大全| 纵有疾风起免费观看全集完整版| 男女边摸边吃奶| 99精品久久久久人妻精品| 久久久久精品久久久久真实原创| 人人妻人人澡人人看| 汤姆久久久久久久影院中文字幕| 亚洲一卡2卡3卡4卡5卡精品中文| 国产精品 国内视频| 高清黄色对白视频在线免费看| 大话2 男鬼变身卡| 美女脱内裤让男人舔精品视频| 国产毛片在线视频| 熟女av电影| 男女床上黄色一级片免费看| 大香蕉久久网| 免费观看人在逋| 国产在线免费精品| 黄片小视频在线播放| 国产色婷婷99| 妹子高潮喷水视频| 在线 av 中文字幕| 日日爽夜夜爽网站| 蜜桃在线观看..| 日韩 亚洲 欧美在线| 99久国产av精品国产电影| 一本大道久久a久久精品| 欧美 亚洲 国产 日韩一| 国产精品免费视频内射| 久久99热这里只频精品6学生| 免费不卡黄色视频| 中文欧美无线码| av.在线天堂| 大香蕉久久网| 久久毛片免费看一区二区三区| 久久久久久人妻| 2021少妇久久久久久久久久久| 91老司机精品| 桃花免费在线播放| 亚洲av国产av综合av卡| 国产欧美亚洲国产| 高清视频免费观看一区二区| 十八禁高潮呻吟视频| 国产一级毛片在线| 日日啪夜夜爽| 久久av网站| 深夜精品福利| 久久人人97超碰香蕉20202| 欧美日韩国产mv在线观看视频| xxxhd国产人妻xxx| 无限看片的www在线观看| 国产精品嫩草影院av在线观看| 日韩一区二区三区影片| 中文字幕亚洲精品专区| 国产 精品1| 老司机影院成人| 少妇 在线观看| 亚洲欧美激情在线| 精品国产超薄肉色丝袜足j| 精品卡一卡二卡四卡免费| 日韩伦理黄色片| 满18在线观看网站| 亚洲av欧美aⅴ国产| 亚洲欧美一区二区三区黑人| 精品国产露脸久久av麻豆| videosex国产| 亚洲人成网站在线观看播放| 一级爰片在线观看| 九色亚洲精品在线播放| 亚洲精品第二区| 免费日韩欧美在线观看| 免费日韩欧美在线观看| 飞空精品影院首页| 国产精品久久久久久精品古装| 国产成人欧美| 久久国产精品男人的天堂亚洲| 久久久精品免费免费高清| 国产探花极品一区二区| av有码第一页| 精品少妇黑人巨大在线播放| 老熟女久久久| 人人澡人人妻人| 日韩中文字幕视频在线看片| svipshipincom国产片| 爱豆传媒免费全集在线观看| 最新的欧美精品一区二区| 如日韩欧美国产精品一区二区三区| 亚洲av国产av综合av卡| 国产 一区精品| 亚洲精品自拍成人| 伦理电影大哥的女人| 日韩欧美一区视频在线观看| 欧美97在线视频| 亚洲精品美女久久av网站| 亚洲精品美女久久久久99蜜臀 | 在现免费观看毛片| netflix在线观看网站| 王馨瑶露胸无遮挡在线观看| 99久久99久久久精品蜜桃| 91精品伊人久久大香线蕉| 99久久综合免费| 久久97久久精品| 大码成人一级视频| 我的亚洲天堂| 久久精品亚洲熟妇少妇任你| 女人爽到高潮嗷嗷叫在线视频| 最近中文字幕高清免费大全6| 男女高潮啪啪啪动态图| 男女边吃奶边做爰视频| 日日撸夜夜添| 亚洲国产av影院在线观看| av免费观看日本| 丝袜人妻中文字幕| 久久免费观看电影| 久久久久久久久久久久大奶| 国产成人精品无人区| 男女之事视频高清在线观看 | 国产老妇伦熟女老妇高清| 国产一区二区 视频在线| 国产成人a∨麻豆精品| 一本久久精品| 建设人人有责人人尽责人人享有的| 亚洲av成人不卡在线观看播放网 | 美女中出高潮动态图| 美女中出高潮动态图| 黄色一级大片看看| 国产有黄有色有爽视频| 高清黄色对白视频在线免费看| 亚洲一码二码三码区别大吗| 国产精品无大码| 日韩中文字幕视频在线看片| 日日摸夜夜添夜夜爱| 一边摸一边做爽爽视频免费| 狂野欧美激情性bbbbbb| 亚洲色图综合在线观看| 亚洲,一卡二卡三卡| 午夜福利影视在线免费观看| 男人舔女人的私密视频| 久久精品aⅴ一区二区三区四区| 亚洲av日韩精品久久久久久密 | 免费看av在线观看网站| 韩国高清视频一区二区三区| 纯流量卡能插随身wifi吗| 久久久久久免费高清国产稀缺| 97精品久久久久久久久久精品| videos熟女内射| 老司机影院毛片| 国产 一区精品| 国产精品秋霞免费鲁丝片| 色视频在线一区二区三区| av国产久精品久网站免费入址| 欧美日韩av久久| 亚洲美女搞黄在线观看| 国产精品国产av在线观看| 老汉色av国产亚洲站长工具| 在线精品无人区一区二区三| 你懂的网址亚洲精品在线观看| 欧美日韩成人在线一区二区| 精品午夜福利在线看| 卡戴珊不雅视频在线播放| av不卡在线播放| 亚洲欧美精品综合一区二区三区| 精品少妇一区二区三区视频日本电影 | 又大又爽又粗| 成人国语在线视频| 亚洲成国产人片在线观看| 国产激情久久老熟女| 综合色丁香网| 欧美日韩视频高清一区二区三区二| 少妇的丰满在线观看| 2018国产大陆天天弄谢| 国产成人免费观看mmmm| 国产精品国产三级专区第一集| 欧美日韩一区二区视频在线观看视频在线| 中文字幕人妻丝袜制服| 香蕉丝袜av| 久久人妻熟女aⅴ| 欧美日韩成人在线一区二区| 久久天躁狠狠躁夜夜2o2o | 90打野战视频偷拍视频| 女性生殖器流出的白浆| 夫妻性生交免费视频一级片| 精品国产一区二区三区久久久樱花| 免费观看av网站的网址| 大话2 男鬼变身卡| 亚洲欧美一区二区三区黑人| 岛国毛片在线播放| 国产成人精品福利久久| 99国产综合亚洲精品| 午夜福利网站1000一区二区三区| 欧美精品人与动牲交sv欧美| 欧美精品人与动牲交sv欧美| 亚洲av国产av综合av卡| 午夜av观看不卡| av免费观看日本| 日韩伦理黄色片| 可以免费在线观看a视频的电影网站 | 伦理电影免费视频| 国产av国产精品国产| 亚洲av综合色区一区| 久久精品国产亚洲av涩爱| 麻豆av在线久日| 亚洲成国产人片在线观看| 国产精品欧美亚洲77777| 91精品国产国语对白视频| 一本久久精品| 自线自在国产av| 九草在线视频观看| 91国产中文字幕| 男女边吃奶边做爰视频| 亚洲四区av| 亚洲成人一二三区av| 精品亚洲乱码少妇综合久久| 亚洲第一青青草原| 国产精品久久久久久精品电影小说| 大片免费播放器 马上看| 高清欧美精品videossex| 欧美精品亚洲一区二区| 美国免费a级毛片| 建设人人有责人人尽责人人享有的| 国产精品久久久av美女十八| 一本一本久久a久久精品综合妖精| 男女边摸边吃奶| 国产精品 欧美亚洲| 久久久久精品性色| 国产深夜福利视频在线观看| 亚洲欧美成人精品一区二区| 国产在线视频一区二区| 午夜福利网站1000一区二区三区| 午夜av观看不卡| 美女扒开内裤让男人捅视频| 久久久久国产一级毛片高清牌| 少妇猛男粗大的猛烈进出视频| 亚洲av综合色区一区| 国产精品国产av在线观看| 亚洲七黄色美女视频| 国产 一区精品| 亚洲成国产人片在线观看| 中国三级夫妇交换| 天美传媒精品一区二区| 下体分泌物呈黄色| 久久狼人影院| 欧美精品一区二区大全| 大码成人一级视频| a级毛片黄视频| av在线老鸭窝| 美女高潮到喷水免费观看| 中文欧美无线码| 最新的欧美精品一区二区| 久久久国产一区二区| 91老司机精品| 超碰97精品在线观看| 久久久欧美国产精品| 欧美精品人与动牲交sv欧美| 欧美在线一区亚洲| 操出白浆在线播放| 日本一区二区免费在线视频| 成人18禁高潮啪啪吃奶动态图| 国产精品嫩草影院av在线观看| 亚洲图色成人| 国产探花极品一区二区| 最新的欧美精品一区二区| 亚洲国产av影院在线观看| av在线播放精品| 飞空精品影院首页| 精品一品国产午夜福利视频| 亚洲国产日韩一区二区| 国产精品久久久久久人妻精品电影 | 亚洲精品中文字幕在线视频| 在线观看一区二区三区激情| 在线精品无人区一区二区三| av网站免费在线观看视频| 满18在线观看网站| 亚洲人成电影观看| 日日撸夜夜添| 中文字幕av电影在线播放| 国产黄色视频一区二区在线观看| 亚洲av综合色区一区| 日韩电影二区| 久久精品人人爽人人爽视色| 99精品久久久久人妻精品| 多毛熟女@视频| 中文字幕精品免费在线观看视频| 建设人人有责人人尽责人人享有的| 欧美在线黄色| 最近中文字幕2019免费版| 69精品国产乱码久久久| 午夜激情久久久久久久| 成人影院久久| 水蜜桃什么品种好| 青春草视频在线免费观看| 国产高清不卡午夜福利| av片东京热男人的天堂| av卡一久久| 在线观看人妻少妇| av网站免费在线观看视频| 日本爱情动作片www.在线观看| 欧美精品亚洲一区二区| 亚洲av成人不卡在线观看播放网 | 欧美日韩亚洲国产一区二区在线观看 | 国产日韩欧美视频二区| 日韩视频在线欧美| 国产伦人伦偷精品视频| 国产一区亚洲一区在线观看| 97人妻天天添夜夜摸| 高清视频免费观看一区二区| 操出白浆在线播放| 91精品国产国语对白视频| 国产欧美日韩一区二区三区在线| 久久狼人影院| 叶爱在线成人免费视频播放| 亚洲美女视频黄频| 91精品国产国语对白视频| 亚洲国产毛片av蜜桃av| 日韩av在线免费看完整版不卡| av一本久久久久| 麻豆精品久久久久久蜜桃| 欧美av亚洲av综合av国产av | 国产精品三级大全| 97人妻天天添夜夜摸| h视频一区二区三区| 国产免费现黄频在线看| 操出白浆在线播放| 女人爽到高潮嗷嗷叫在线视频| 日韩精品免费视频一区二区三区| 这个男人来自地球电影免费观看 | 午夜福利网站1000一区二区三区| 久久这里只有精品19| 久久久久网色| 免费少妇av软件| 一级毛片黄色毛片免费观看视频| 成人国语在线视频| 国产片特级美女逼逼视频| 亚洲,欧美,日韩| 亚洲欧美成人综合另类久久久| 一本久久精品| 久久久欧美国产精品| 中文欧美无线码| 啦啦啦视频在线资源免费观看| 免费看av在线观看网站| 欧美人与善性xxx| 交换朋友夫妻互换小说| 国产成人a∨麻豆精品| 巨乳人妻的诱惑在线观看| 黄色视频在线播放观看不卡| 晚上一个人看的免费电影| 女人久久www免费人成看片| 国产精品秋霞免费鲁丝片| 男男h啪啪无遮挡| 老熟女久久久| 成年美女黄网站色视频大全免费| av一本久久久久| 国产精品久久久久久精品古装| 欧美激情极品国产一区二区三区| 母亲3免费完整高清在线观看| 黑人猛操日本美女一级片| 在线免费观看不下载黄p国产| 亚洲激情五月婷婷啪啪| 日韩精品免费视频一区二区三区| 久久久久久久久久久久大奶| 夫妻午夜视频| 黄色 视频免费看| 久久人妻熟女aⅴ| 久久99热这里只频精品6学生| 一个人免费看片子| 成人手机av| 天天躁日日躁夜夜躁夜夜| 欧美激情 高清一区二区三区| 亚洲伊人久久精品综合| 精品卡一卡二卡四卡免费| 日本91视频免费播放| 精品免费久久久久久久清纯 | 国产爽快片一区二区三区| av网站免费在线观看视频| 一区二区三区四区激情视频| 国产在线一区二区三区精| 精品国产乱码久久久久久男人| svipshipincom国产片| 男人舔女人的私密视频| 欧美在线一区亚洲| 黄色毛片三级朝国网站| 国产成人系列免费观看| 最新的欧美精品一区二区| 美女福利国产在线| 久久精品国产a三级三级三级| 亚洲成人av在线免费| 大码成人一级视频| 一二三四在线观看免费中文在| 亚洲成人国产一区在线观看 | 亚洲国产精品一区二区三区在线| 欧美少妇被猛烈插入视频| 欧美亚洲日本最大视频资源| 国产成人精品久久久久久| 黄色怎么调成土黄色| 精品视频人人做人人爽| 别揉我奶头~嗯~啊~动态视频 | 看免费成人av毛片| 国产成人欧美在线观看 | 男女下面插进去视频免费观看| 黄色怎么调成土黄色| 国产色婷婷99| 日本一区二区免费在线视频| 国产精品99久久99久久久不卡 | 久久热在线av| 少妇精品久久久久久久| 中文欧美无线码| 欧美精品亚洲一区二区| 欧美日韩亚洲综合一区二区三区_| 成人影院久久| 久久精品国产a三级三级三级| 伊人亚洲综合成人网| 亚洲综合色网址| 亚洲欧美成人综合另类久久久| 亚洲欧美清纯卡通| 久久久精品免费免费高清| 另类精品久久| 男女之事视频高清在线观看 | 久久精品国产亚洲av高清一级| 欧美精品av麻豆av| 久久久久久人妻| 国产成人精品久久二区二区91 | 国产精品久久久久成人av| 一区二区日韩欧美中文字幕| 免费观看a级毛片全部| 王馨瑶露胸无遮挡在线观看| 久久婷婷青草| 99精品久久久久人妻精品| 老司机深夜福利视频在线观看 | 日韩免费高清中文字幕av| 日韩电影二区| 观看av在线不卡| 国产精品久久久久成人av| 欧美老熟妇乱子伦牲交| 精品一区二区免费观看| 亚洲欧美一区二区三区国产| 在线天堂中文资源库| 777米奇影视久久| 建设人人有责人人尽责人人享有的| 午夜福利免费观看在线| 国产人伦9x9x在线观看| 伊人久久国产一区二区| 欧美精品人与动牲交sv欧美| 91精品伊人久久大香线蕉| 人妻 亚洲 视频| 日韩中文字幕视频在线看片| 18禁观看日本| 亚洲美女黄色视频免费看| 亚洲国产精品成人久久小说| 免费高清在线观看日韩| 黑人欧美特级aaaaaa片| 新久久久久国产一级毛片| 一区二区三区精品91| 国产精品久久久久成人av| 亚洲成av片中文字幕在线观看| av网站在线播放免费| 黑人巨大精品欧美一区二区蜜桃| 这个男人来自地球电影免费观看 | 成人亚洲精品一区在线观看| 久久99热这里只频精品6学生| 午夜免费观看性视频| 少妇人妻久久综合中文| 一二三四在线观看免费中文在| 亚洲成色77777| 久久久久人妻精品一区果冻| 久久精品熟女亚洲av麻豆精品| 久久久久精品久久久久真实原创| 国产日韩欧美视频二区| 视频在线观看一区二区三区| 欧美日韩国产mv在线观看视频| 免费在线观看完整版高清| 狂野欧美激情性xxxx| 国产老妇伦熟女老妇高清| 黄色一级大片看看| av一本久久久久| h视频一区二区三区| 精品一区二区三区四区五区乱码 | 2018国产大陆天天弄谢| 免费黄频网站在线观看国产| 日韩制服骚丝袜av| 亚洲国产成人一精品久久久| 久久久久久久久久久久大奶| 亚洲婷婷狠狠爱综合网| 黄片小视频在线播放| 一边亲一边摸免费视频| 九九爱精品视频在线观看| 色精品久久人妻99蜜桃| 午夜91福利影院| 国产成人精品福利久久| 亚洲欧洲日产国产| 丁香六月天网| 国产麻豆69| 国产欧美日韩一区二区三区在线| 国产精品成人在线| 99热全是精品| 叶爱在线成人免费视频播放| 极品人妻少妇av视频| 婷婷色av中文字幕| 电影成人av| 黄色视频在线播放观看不卡| 中文天堂在线官网| 夜夜骑夜夜射夜夜干| 国产成人午夜福利电影在线观看| 日韩制服骚丝袜av| 丰满饥渴人妻一区二区三| 国产欧美亚洲国产| 国产一区亚洲一区在线观看| 天堂俺去俺来也www色官网| 最近最新中文字幕大全免费视频 | 精品第一国产精品| 成人午夜精彩视频在线观看| 国产男女内射视频| 日韩精品免费视频一区二区三区| 亚洲av福利一区| 无遮挡黄片免费观看| 亚洲欧美成人综合另类久久久| 亚洲精品中文字幕在线视频| 午夜日本视频在线| 国产老妇伦熟女老妇高清| 狠狠婷婷综合久久久久久88av| 亚洲婷婷狠狠爱综合网| 久久女婷五月综合色啪小说| 国产成人免费无遮挡视频| 男女床上黄色一级片免费看| 在线观看国产h片| 久久国产精品男人的天堂亚洲| 亚洲美女视频黄频| 国产精品蜜桃在线观看| 欧美日韩视频精品一区| 韩国高清视频一区二区三区| 亚洲激情五月婷婷啪啪| 水蜜桃什么品种好| 亚洲av欧美aⅴ国产| 国产片特级美女逼逼视频| 男女之事视频高清在线观看 | 欧美人与性动交α欧美精品济南到| 欧美精品亚洲一区二区| 大话2 男鬼变身卡| 2021少妇久久久久久久久久久| av.在线天堂| 成年人午夜在线观看视频| 天天躁夜夜躁狠狠久久av| 婷婷色综合www| 国产熟女欧美一区二区| 久久韩国三级中文字幕| 亚洲在久久综合| 精品一区二区三区av网在线观看 | 丝袜脚勾引网站| 国产福利在线免费观看视频| 高清欧美精品videossex| 叶爱在线成人免费视频播放| 亚洲婷婷狠狠爱综合网| videosex国产| 91精品伊人久久大香线蕉| 19禁男女啪啪无遮挡网站| 久久久久久久精品精品| 99re6热这里在线精品视频| 久久久久久久久免费视频了| 欧美激情 高清一区二区三区| 夫妻午夜视频| 日韩精品免费视频一区二区三区| 777米奇影视久久| 国产色婷婷99| 色综合欧美亚洲国产小说| 欧美在线黄色| 大片电影免费在线观看免费| 美女中出高潮动态图| 国产成人啪精品午夜网站| avwww免费| 男女高潮啪啪啪动态图| 午夜福利乱码中文字幕| 大片电影免费在线观看免费| 另类亚洲欧美激情| 亚洲av在线观看美女高潮| 99香蕉大伊视频| av线在线观看网站| 少妇人妻 视频| 乱人伦中国视频| 国产97色在线日韩免费| 欧美av亚洲av综合av国产av | 欧美日韩亚洲高清精品| 精品一区二区免费观看| 美女午夜性视频免费|