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

    帶上水線拉索繞流場(chǎng)的大渦模擬研究

    2016-09-20 07:28:18杜曉慶李俊軍
    關(guān)鍵詞:水線升力風(fēng)壓

    杜曉慶, 李俊軍, 顧 明, 何 平

    (1.上海大學(xué) 土木工程系,上海 200072;2.同濟(jì)大學(xué) 土木工程防災(zāi)國(guó)家重點(diǎn)實(shí)驗(yàn)室,上海 200092)

    ?

    帶上水線拉索繞流場(chǎng)的大渦模擬研究

    杜曉慶1, 李俊軍1, 顧明2, 何平1

    (1.上海大學(xué) 土木工程系,上海 200072;2.同濟(jì)大學(xué) 土木工程防災(zāi)國(guó)家重點(diǎn)實(shí)驗(yàn)室,上海 200092)

    針對(duì)帶上水線的拉索模型,采用大渦模擬方法研究了水線和拉索的氣動(dòng)性能與繞流場(chǎng)流態(tài)之間的內(nèi)在關(guān)系.結(jié)果表明,帶水線拉索的繞流場(chǎng)呈現(xiàn)三種不同的流態(tài)結(jié)構(gòu).當(dāng)上水線位于0°~50°時(shí),水線對(duì)拉索繞流場(chǎng)的影響不大;當(dāng)上水線位于50°~70°時(shí),拉索上表面剪切層在上水線上分離后會(huì)再附著在拉索表面形成分離泡,并使拉索上表面受到很大的局部負(fù)壓作用;當(dāng)水線處在70°~90°時(shí),凸起的上水線會(huì)導(dǎo)致上側(cè)剪切層分離后遠(yuǎn)離拉索上表面.分離泡在拉索上表面周期性地出現(xiàn)和消失是斜拉索風(fēng)雨激振的主要原因.

    拉索; 風(fēng)雨激振; 大渦模擬; 上水線; 分離泡

    自1989年日本的Hikami和Shiraishi[1]首先報(bào)道了斜拉橋拉索的風(fēng)雨激振現(xiàn)象以來(lái),世界各地的許多斜拉橋拉索在風(fēng)雨共同作用下發(fā)生過(guò)這種大幅振動(dòng)[2-4].斜拉索的大幅振動(dòng)會(huì)嚴(yán)重危害拉索和斜拉橋的安全性.

    各國(guó)研究人員通過(guò)現(xiàn)場(chǎng)實(shí)測(cè)、風(fēng)洞試驗(yàn)和理論分析對(duì)拉索風(fēng)雨激振現(xiàn)象進(jìn)行較為廣泛的研究,已獲得了很多研究成果.但由于拉索風(fēng)雨激振涉及空氣、液體和固體等三種介質(zhì)的耦合作用,并且受拉索表面條件、拉索傾角、風(fēng)速、風(fēng)向和雨量等多種因素的影響,因而對(duì)風(fēng)雨激振的發(fā)生機(jī)理尚有不同觀點(diǎn).大部分研究者[1,5-8]認(rèn)為,上水線在拉索上表面的出現(xiàn)及其沿拉索表面的振蕩運(yùn)動(dòng)是拉索發(fā)生這種氣動(dòng)失穩(wěn)現(xiàn)象的主要原因,可采用基于準(zhǔn)定常假定的理論模型來(lái)解釋風(fēng)雨激振現(xiàn)象.但也有學(xué)者[2]認(rèn)為,軸向流是斜拉索發(fā)生氣動(dòng)失穩(wěn)的原因,認(rèn)為準(zhǔn)定常假定不適用風(fēng)雨激振現(xiàn)象.作者對(duì)偏/斜光圓柱和帶水線偏/斜圓柱進(jìn)行的研究表明:軸向流會(huì)對(duì)拉索的氣動(dòng)性能產(chǎn)生一定的影響,但一旦上水線在拉索上表面形成之后,上水線對(duì)拉索的氣動(dòng)性能將起決定性的作用[9].

    目前對(duì)上水線的受力情況和形成過(guò)程研究較少,上水線對(duì)拉索繞流場(chǎng)流態(tài)的作用機(jī)理尚不清楚.Du和Gu等[9]采用風(fēng)洞試驗(yàn),研究了上水線處在不同位置時(shí)上水線的表面風(fēng)壓和氣動(dòng)力作用規(guī)律.Lemaitre[10]基于滑移理論模型,研究了在風(fēng)荷載作用下拉索表面水膜形態(tài)的變化.許林汕等[11]在Lemaitre模型的基礎(chǔ)上,研究了考慮拉索運(yùn)動(dòng)影響的水膜形態(tài)變化和運(yùn)動(dòng)規(guī)律.Bi等[12]采用氣液兩相流理論和流體體積法研究了水線的形成過(guò)程.杜曉慶等[13]采用大渦模擬(LES)方法研究了作用在上水線上的風(fēng)致摩擦力特性,并分析了水線的表面風(fēng)壓和摩擦力在上水線形成過(guò)程中的作用.

    為了進(jìn)一步澄清水線和拉索氣動(dòng)性能與繞流場(chǎng)流態(tài)之間的內(nèi)在關(guān)系,本文采用大渦模擬方法,以帶上水線拉索模型的氣動(dòng)性能和繞流場(chǎng)特性為研究對(duì)象,在亞臨界區(qū)(雷諾數(shù)為1.4×105),研究了帶水線拉索的氣動(dòng)力、表面風(fēng)壓分布和繞流場(chǎng)特性隨水線位置的演變過(guò)程,研究了流態(tài)結(jié)構(gòu)與水線和拉索氣動(dòng)力之間的內(nèi)在聯(lián)系,從帶水線拉索繞流場(chǎng)的流態(tài)變化角度對(duì)拉索風(fēng)雨激振的發(fā)生機(jī)理進(jìn)行了探討.

    1 數(shù)值方法和計(jì)算模型

    1.1控制方程和亞格子模型

    大渦模擬(LES)方法中,大尺度渦通過(guò)濾波后的Navier-Stokes方程直接求解,小尺度的渦則采用亞格子尺度模型(SGS)模擬.與雷諾平均法(RANS)相比,大渦模擬方法可更好地模擬流場(chǎng)中的湍流漩渦,能捕捉到豐富的流場(chǎng)脈動(dòng)信息.因此,為了準(zhǔn)確地模擬帶水線拉索周圍的流場(chǎng)特性,本文采用大渦模擬方法進(jìn)行研究.

    經(jīng)過(guò)濾波函數(shù)的濾波,可得到大尺度渦的不可壓縮Navier-Stokes方程

    (1)

    (2)

    亞格子應(yīng)力τij采用Smagorinsky-Lilly的亞格子尺度模型,具有以下形式:

    (3)

    (4)

    式中νt是亞格子尺度的湍動(dòng)黏度,且

    (5)

    (6)

    (7)

    式中:Δi代表沿坐標(biāo)軸i方向的網(wǎng)格尺寸;CS是Smagorinsky常數(shù),本文取0.1.

    1.2計(jì)算模型

    考慮到上水線的出現(xiàn)是拉索發(fā)生風(fēng)雨激振的關(guān)鍵因素,本文的計(jì)算模型只考慮上水線的作用,假定水線的形狀和尺寸不隨水線的位置發(fā)生改變.現(xiàn)場(chǎng)實(shí)測(cè)發(fā)現(xiàn)斜拉橋拉索發(fā)生大幅振動(dòng)時(shí)雷諾數(shù)為0.5×105~2×105,因而本文的計(jì)算雷諾數(shù)取1.4×105(雷諾數(shù)根據(jù)拉索直徑D和來(lái)流風(fēng)速U0計(jì)算得到).另外,在本文的計(jì)算模型中,來(lái)流垂直作用于帶水線拉索模型,即不考慮風(fēng)向角的影響.

    圖1為帶上水線拉索的計(jì)算模型示意圖.圖中,θu為上水線的位置角,θcyl為拉索表面各點(diǎn)的位置角,CD為拉索的平均阻力系數(shù),CL為拉索的平均升力系數(shù).上水線模型的外形和尺寸與文獻(xiàn)[9]的風(fēng)洞試驗(yàn)?zāi)P拖嘟?,水線表面為圓弧曲面,圓弧的直徑d=0.13D,水線寬度為0.1D,高度為0.024D.

    圖1 帶上水線拉索的計(jì)算模型

    為了掌握上水線位置對(duì)拉索和水線氣動(dòng)性能的影響,根據(jù)上水線角度位置的不同,共計(jì)算了18個(gè)工況.其中:上水線位置角θu在0°~50°和80°~110° 之間每隔10°一個(gè)工況,50°~70°之間每隔5°一個(gè)工況,56°~64°之間每隔2°一個(gè)工況.

    計(jì)算模型如圖2所示,計(jì)算域的長(zhǎng)、寬、高分別為28D,16D和D.其中,拉索中心距入口邊界為8D,距離出口邊界為20D,距上下側(cè)面邊界為8D,拉索展向長(zhǎng)度為D.采用速度入口邊界條件,出口邊界采用自由流出邊界,計(jì)算域上下側(cè)面為對(duì)稱性邊界,拉索展向兩端采用周期性邊界條件,水線和拉索表面采用無(wú)滑移壁面邊界條件.

    本文的網(wǎng)格方案和計(jì)算參數(shù)基于文獻(xiàn)[16]的研究成果,近壁面最小網(wǎng)格為4×10-4D,量綱一時(shí)間步Δt*為 0.001(Δt*=ΔtU0/D,其中Δt為實(shí)際計(jì)算時(shí)間步),近壁面網(wǎng)格增長(zhǎng)率為1.05.另外,為了驗(yàn)證圓柱展向長(zhǎng)度對(duì)計(jì)算結(jié)果的影響,對(duì)三種不同的圓柱展向長(zhǎng)度進(jìn)行了比較(D,2D和4D).圖3為這三種不同展向長(zhǎng)度的圓柱表面平均風(fēng)壓系數(shù)Cp分布圖,圖中還給出了Cantwell和Coles[17]通過(guò)風(fēng)洞試驗(yàn)得到的平均風(fēng)壓系數(shù).從圖3可見(jiàn),三種計(jì)算值與風(fēng)洞試驗(yàn)值之間的差異主要發(fā)生在θcyl=240°~360°,即圓柱的背風(fēng)側(cè).并且隨著圓柱展向長(zhǎng)度的減小,計(jì)算值與風(fēng)洞試驗(yàn)值之間的偏差有增大的趨勢(shì).這主要是因?yàn)閳A柱展向長(zhǎng)度的減小會(huì)增強(qiáng)繞流場(chǎng)沿圓柱展向的相關(guān)性,不能很好地反映三維圓柱繞流的三維流場(chǎng)特性,從而會(huì)造成圓柱背風(fēng)面的負(fù)壓偏大.考慮到大渦模擬的計(jì)算量大,而本文的計(jì)算工況又多,為了降低計(jì)算量,本文的帶水線圓柱模型的展向長(zhǎng)度取D,與文獻(xiàn)[18]的計(jì)算模型相似.

    a 計(jì)算域

    b 近壁面網(wǎng)格

    圖3 光圓柱表面風(fēng)壓系數(shù)比較

    1.3氣動(dòng)參數(shù)定義

    拉索和水線表面的風(fēng)壓系數(shù)定義為

    (8)

    式中:p-po為當(dāng)?shù)仫L(fēng)壓和遠(yuǎn)前方上游壓力之差.

    對(duì)拉索表面的風(fēng)壓沿拉索周向進(jìn)行積分,可得到作用在拉索上的氣動(dòng)阻力和氣動(dòng)升力,對(duì)氣動(dòng)力按照式(9)和式(10)量綱一化后可得到拉索的平均阻力系數(shù)CD和平均升力系數(shù)CL(其方向定義見(jiàn)圖1)

    (9)

    (10)

    式中:FD和FL分別為作用在單位長(zhǎng)拉索上的氣動(dòng)平均阻力和氣動(dòng)平均升力.

    2 計(jì)算結(jié)果及分析

    2.1拉索的氣動(dòng)力系數(shù)

    圖4為作用在帶上水線拉索上的平均阻力系數(shù)CL和平均升力系數(shù)CD隨上水線位置角θu變化的曲線.圖中還列出了文獻(xiàn)[7,9]中的風(fēng)洞試驗(yàn)結(jié)果.

    從圖4可見(jiàn),本文采用大渦模擬計(jì)算得到的氣動(dòng)力系數(shù)曲線的變化趨勢(shì)與文獻(xiàn)[7,9]中的風(fēng)洞試驗(yàn)結(jié)果相似.隨著水線位置角的逐步增大,拉索的平均阻力系數(shù)和平均升力系數(shù)均經(jīng)歷了起伏和突變過(guò)程.與平均阻力系數(shù)相比,平均升力系數(shù)的變化則更為劇烈,經(jīng)歷了平穩(wěn)段、上升段和下降段三個(gè)階段.當(dāng)水線位于θu=0°~50°時(shí),拉索的平均升力系數(shù)接近于0,與不帶水線的光滑拉索的平均升力系數(shù)相同;當(dāng)水線位于θu=50°~70°時(shí),平均升力系數(shù)則經(jīng)歷了先迅速上升,然后突然下降的過(guò)程,最大平均升力系數(shù)達(dá)到0.55左右;而當(dāng)水線位于θu=70°~90°時(shí),拉索則受到方向相反的非零平均升力的作用,最小升力系數(shù)為-0.15.通過(guò)風(fēng)洞試驗(yàn)得到的文獻(xiàn)[7,9]的平均升力系數(shù)同樣經(jīng)歷了類似的變化過(guò)程,只是發(fā)生突變的水線位置范圍以及氣動(dòng)力系數(shù)的數(shù)值與本文大渦模擬計(jì)算結(jié)果有所差異.

    需要指出的是,與風(fēng)洞試驗(yàn)值相比,本文采用大渦模擬得到的拉索平均升力系數(shù)與試驗(yàn)值吻合較好,但平均阻力系數(shù)的計(jì)算值較試驗(yàn)值偏大.這可能與計(jì)算模型與試驗(yàn)?zāi)P退捎玫耐牧鞫?、表面粗糙度、雷諾數(shù)等的差異有關(guān).Zdravkovich[14]指出圓柱的氣動(dòng)性能對(duì)來(lái)流湍流度、雷諾數(shù)和表面粗糙度等條件非常敏感.即便文獻(xiàn)[7,9]均為風(fēng)洞試驗(yàn)的結(jié)果,兩者的氣動(dòng)力系數(shù)也有相當(dāng)大的差異.另一方面,如1.2節(jié)所述,本文計(jì)算模型采用較小的展向長(zhǎng)度會(huì)造成圓柱背風(fēng)面負(fù)壓的增大,從而會(huì)導(dǎo)致計(jì)算得到的平均阻力系數(shù)偏大.

    a 平均阻力系數(shù)

    b 平均升力系數(shù)

    2.2拉索的風(fēng)壓系數(shù)

    由上文可知,上水線處在不同位置時(shí),拉索的平均升力系數(shù)會(huì)經(jīng)歷三個(gè)不同階段,圖5~7給出了上水線位于θu=40°,58°和80°時(shí)帶水線拉索表面的平均風(fēng)壓系數(shù)和脈動(dòng)風(fēng)壓系數(shù)分布圖.這三個(gè)典型的水線位置分別處在升力系數(shù)的三個(gè)不同階段中.

    圖5 θu=40°時(shí)拉索表面風(fēng)壓系數(shù)

    Fig.5Wind pressures of cable with rivulet atθu=40°

    圖6 θu=58°時(shí)拉索表面風(fēng)壓系數(shù)

    Fig.6Wind pressures of cable with rivulet atθu=58°

    圖7 θu=80°時(shí)拉索表面風(fēng)壓系數(shù)

    Fig.7Wind pressures of cable with rivulet atθu=80°

    上水線位于θu=40°時(shí)(見(jiàn)圖5),拉索上側(cè)表面(θcyl=0°~180°)的風(fēng)壓系數(shù)與下側(cè)表面(θcyl=180°~360°)相比,主要的差別發(fā)生在上水線位置附近(θu=40°附近),其他位置的平均風(fēng)壓系數(shù)在拉索的上下表面基本對(duì)稱.脈動(dòng)風(fēng)壓系數(shù)在θcyl=80°和280°附近有兩個(gè)峰值,這兩個(gè)峰值處在拉索上下表面的分離點(diǎn)附近,與表面不帶上水線的光滑拉索的情況相近[15].這說(shuō)明水線處在這個(gè)位置時(shí),上水線不影響拉索上側(cè)剪切層的分離點(diǎn)位置,因而對(duì)拉索氣動(dòng)性能的影響不大.

    上水線位于θu=58°時(shí)(見(jiàn)圖6),拉索上下表面的平均風(fēng)壓系數(shù)和脈動(dòng)風(fēng)壓系數(shù)的形態(tài)均差別很大.對(duì)于平均風(fēng)壓系數(shù),在上水線附近拉索受到很大負(fù)壓的作用,拉索上表面負(fù)壓最大處的風(fēng)壓系數(shù)達(dá)到-3.5,遠(yuǎn)大于拉索下側(cè)受到的負(fù)壓;在拉索的背風(fēng)面,拉索上側(cè)背風(fēng)面(即θcyl=90°~180°)受到的負(fù)壓也要遠(yuǎn)大于拉索下側(cè)背風(fēng)面(即θcyl=180°~270°);這種拉索上下表面平均風(fēng)壓系數(shù)的不對(duì)稱性導(dǎo)致了帶上水線拉索會(huì)受到很大的平均升力作用.對(duì)于脈動(dòng)風(fēng)壓系數(shù),拉索上側(cè)表面在θcyl=58°和135°附近分別有兩個(gè)峰值,前一個(gè)峰值對(duì)應(yīng)于上水線位置,也是拉索上側(cè)剪切層的分離點(diǎn);而后一個(gè)峰值則對(duì)應(yīng)于上側(cè)剪切層在拉索表面再附后再次分離的分離點(diǎn),詳細(xì)的流場(chǎng)分析見(jiàn)2.3節(jié);拉索下表面的脈動(dòng)風(fēng)壓系數(shù)則在θcyl=280°附近有一個(gè)峰值,對(duì)應(yīng)于拉索下側(cè)的分離點(diǎn).

    上水線位于θu=80°時(shí)(見(jiàn)圖7),對(duì)于平均風(fēng)壓系數(shù),拉索上下表面風(fēng)壓分布的差別主要體現(xiàn)在水線位置附近,拉索表面風(fēng)壓在其他位置上基本是上下表面對(duì)稱.對(duì)于脈動(dòng)風(fēng)壓系數(shù),拉索上表面的峰值處在水線位置,即拉索上側(cè)的剪切層在上水線上發(fā)生分離;需要說(shuō)明的是,在上水線處在θu=70°~90°時(shí),上表面脈動(dòng)風(fēng)壓系數(shù)的峰值均在水線位置,這說(shuō)明凸起的上水線會(huì)控制拉索上側(cè)剪切層的分離點(diǎn),而拉索下表面的峰值均在θcyl=280°附近.

    2.3作用機(jī)理分析

    為了進(jìn)一步分析帶水線拉索氣動(dòng)性能處在三個(gè)不同階段的流場(chǎng)作用機(jī)理,圖8~10給出了水線位置角θu為40°,58°和80°時(shí)繞流場(chǎng)的瞬態(tài)流線圖、渦量圖、風(fēng)壓場(chǎng)以及拉索表面的瞬態(tài)風(fēng)壓系數(shù)圖.圖中,T1和T2分別是升力系數(shù)時(shí)程曲線上升力系數(shù)達(dá)到最小值和最大值的時(shí)刻,即對(duì)應(yīng)于升力系數(shù)時(shí)程曲線的波谷和波峰時(shí)刻.

    當(dāng)水線位于θu=40°時(shí),對(duì)于T1時(shí)刻(見(jiàn)圖8a),從渦量圖上看,拉索下側(cè)的剪切流在分離后、近尾流區(qū)形成了強(qiáng)度很大的漩渦,從而導(dǎo)致拉索下側(cè)繞流場(chǎng)處在負(fù)壓區(qū)內(nèi),拉索下側(cè)表面(θcyl=180°~360°)則承受很大的負(fù)壓,最大負(fù)壓系數(shù)達(dá)到-3.0;而在拉索的上側(cè),由于剪切流剛從拉索表面分離,還未能充分發(fā)展形成強(qiáng)度很大的漩渦,因而拉索上側(cè)表面(θcyl=0°~180°)的最大負(fù)壓系數(shù)較小,僅為-1.0左右.對(duì)于T2時(shí)刻(見(jiàn)圖8b),則拉索上下側(cè)的流場(chǎng)與T1時(shí)刻相反,拉索下側(cè)表面的最大負(fù)壓系數(shù)約為-1.0,而拉索上側(cè)表面最大負(fù)壓系數(shù)為-2.8左右.從T1和T2時(shí)刻拉索繞流場(chǎng)的瞬態(tài)流場(chǎng)信息和拉索表面的風(fēng)壓系數(shù)可見(jiàn):在升力系數(shù)時(shí)程曲線的波峰和波谷時(shí)刻,帶水線周圍流場(chǎng)和拉索上下側(cè)表面的風(fēng)壓系數(shù)基本處在一種反對(duì)稱的狀態(tài),即上水線的存在并未對(duì)拉索的繞流場(chǎng)和氣動(dòng)性能造成很大影響,因而得到的帶水線拉索的平均升力系數(shù)約為0.

    a T1時(shí)刻

    b T2時(shí)刻

    Fig.8Instantaneous streamline, vorticity, pressure field around and pressure distribution of cable with rivulet atθu=40°

    a T1時(shí)刻

    b T2時(shí)刻

    Fig.9Instantaneous streamline, vorticity, pressure field around and pressure distribution of cable with rivulet atθu=58°

    當(dāng)水線位于θu=58°時(shí),對(duì)應(yīng)于T1時(shí)刻(見(jiàn)圖9a),從渦量圖上看,拉索下側(cè)的剪切流在分離后、近尾流區(qū)形成漩渦,但漩渦中心離圓柱表面的距離比θu=40°(圖8a)時(shí)遠(yuǎn),因而拉索下側(cè)表面所受到的負(fù)壓數(shù)值也較θu=40°時(shí)??;而對(duì)于拉索的上側(cè)表面,從瞬態(tài)流線圖和渦量圖可以發(fā)現(xiàn):由于受到上水線的干擾,上側(cè)剪切層從拉索上表面的水線上分離后會(huì)再附到拉索表面,即形成 “分離泡”現(xiàn)象,使分離泡附近的拉索上表面承受很強(qiáng)的負(fù)壓作用,最大的負(fù)壓系數(shù)接近于-4.0.對(duì)于T2時(shí)刻(見(jiàn)圖9b),拉索上表面同樣發(fā)生了剪切流的分離再附現(xiàn)象,拉索上側(cè)的繞流場(chǎng)呈現(xiàn)大面積的負(fù)壓區(qū);受到分離泡和尾流漩渦的共同影響,拉索上表面也出現(xiàn)數(shù)值很大的負(fù)壓作用,最大負(fù)壓系數(shù)已超過(guò)了-4.0.從上述分析可見(jiàn),與上水線處在40°時(shí)相比,處在58°的上水線對(duì)拉索的繞流場(chǎng)和氣動(dòng)性能造成了很大影響,并且使得在一個(gè)渦脫周期內(nèi),拉索的上表面均受到數(shù)值很大的負(fù)壓作用,從而導(dǎo)致帶上水線拉索會(huì)承受較大的平均升力的作用,并且平均升力的方向指向拉索上側(cè).

    當(dāng)水線位于θu=80°時(shí),由圖10的流線圖和渦量圖可見(jiàn),T1和T2時(shí)刻拉索上表面均未發(fā)生分離剪切層的再附現(xiàn)象;T1和T2時(shí)刻拉索周圍的風(fēng)壓場(chǎng)和拉索的表面風(fēng)壓系數(shù)分布圖總體呈現(xiàn)反對(duì)稱的形態(tài),但T1時(shí)刻拉索下表面的負(fù)壓系數(shù)比T2時(shí)刻拉索上表面的負(fù)壓系數(shù)更大,從而導(dǎo)致帶水線拉索受到平均升力系數(shù)的作用,但平均升力系數(shù)的方向指向拉索的下側(cè),其方向與θu=58°時(shí)的升力方向相反.從瞬態(tài)渦量圖上看,這一反方向的平均升力的出現(xiàn)很可能是因?yàn)樯纤€對(duì)拉索上側(cè)剪切層的影響造成的.水線在80°時(shí),凸起的上水線會(huì)導(dǎo)致拉索上側(cè)的剪切層在分離后遠(yuǎn)離拉索表面,因而剪切層對(duì)拉索上表面的影響相應(yīng)減弱,拉索上側(cè)表面受到的負(fù)壓相對(duì)較小,導(dǎo)致拉索受到非零平均升力作用.

    a T1時(shí)刻

    b T2時(shí)刻

    Fig.10Instantaneous streamline, vorticity, pressure field around and pressure distribution of cable with rivulet atθu=80°

    從上文分析可知,當(dāng)水線處在θu=50°~70°時(shí),帶水線拉索的繞流場(chǎng)發(fā)生了劇烈的變化,上水線對(duì)拉索的氣動(dòng)性能也會(huì)造成嚴(yán)重的影響.對(duì)于實(shí)際的斜拉橋拉索,如果在風(fēng)雨共同作用下拉索上表面出現(xiàn)水線后,在漩渦脫落的渦激力和拉索振動(dòng)的牽連慣性力作用下,上水線會(huì)在拉索上表面作振蕩運(yùn)動(dòng),如果水線振蕩時(shí)其位置進(jìn)出在θu=50°~70°這一區(qū)域,則在拉索的上表面會(huì)周期性地出現(xiàn)分離泡現(xiàn)象,拉索上表面分離泡影響區(qū)周期性地受到很大的負(fù)壓作用,從而引發(fā)拉索發(fā)生振動(dòng).拉索的振動(dòng)會(huì)進(jìn)一步增大上水線在拉索表面的振蕩幅度,增強(qiáng)水線振蕩的規(guī)律性,進(jìn)而引發(fā)大振幅的風(fēng)雨激振的發(fā)生.

    3 結(jié)論

    采用大渦模擬方法,以帶上水線拉索模型的氣動(dòng)性能和繞流場(chǎng)特性為研究對(duì)象,在亞臨界區(qū)(雷諾數(shù)為1.4×105),研究了帶水線拉索模型的表面風(fēng)壓系數(shù)、氣動(dòng)力系數(shù)和周圍繞流場(chǎng)特性,研究了流場(chǎng)流態(tài)與帶水線拉索氣動(dòng)性能之間的內(nèi)在聯(lián)系,從帶水線拉索繞流場(chǎng)的流態(tài)變化角度對(duì)拉索風(fēng)雨激振的發(fā)生機(jī)理進(jìn)行了探討.主要結(jié)論如下:

    (1) 上水線處在θu=0°~50°時(shí),上水線不影響拉索上側(cè)剪切層的分離點(diǎn)位置,拉索的氣動(dòng)性能和繞流場(chǎng)特性基本不受上水線的影響.

    (2) 上水線處在θu=50°~70°時(shí),受到上水線的干擾,拉索上側(cè)剪切層從上水線上分離后會(huì)再附到拉索表面,形成分離泡的現(xiàn)象.在分離泡影響范圍內(nèi)的拉索表面受到很大的負(fù)壓作用,并導(dǎo)致拉索受到顯著的非零平均升力系數(shù)的作用,平均升力的方向指向拉索上側(cè).

    (3) 上水線處在θu=70°~90°時(shí),拉索上側(cè)剪切層在上水線上發(fā)生分離,上水線的位置決定拉索上側(cè)剪切層的分離點(diǎn)位置.凸起的上水線會(huì)導(dǎo)致拉索上側(cè)的剪切層在分離后遠(yuǎn)離拉索表面,因而分離剪切層對(duì)拉索氣動(dòng)性能的影響也較弱,拉索上側(cè)表面受到的負(fù)壓相對(duì)較小,帶水線拉索受到較小的非零平均升力的作用,但平均升力的方向指向拉索下側(cè).

    (4) 在漩渦脫落的渦激力和拉索振動(dòng)的牽連慣性力作用下,上水線會(huì)在拉索上表面作振蕩運(yùn)動(dòng),如果水線振蕩時(shí)其位置進(jìn)出在θu=50°~70°這一區(qū)域,則在拉索的上表面會(huì)周期性地出現(xiàn)分離泡現(xiàn)象,拉索上表面分離泡影響區(qū)周期性地受到很大的負(fù)壓作用,從而引發(fā)拉索發(fā)生振動(dòng).拉索的振動(dòng)會(huì)進(jìn)一步增大上水線在拉索表面的振蕩幅度,增強(qiáng)水線振蕩的規(guī)律性,進(jìn)而引發(fā)大振幅的風(fēng)雨激振的發(fā)生.

    [1]Hikami Y, Shiraishi N. Rain-wind induced vibrations of cables in cable stayed bridges [J]. Journal of Wind Engineering and Industrial Aerodynamics, 1988, 29(1): 409.

    [2]Matsumoto M, Shiraishi N. Rain-wind induced vibrations of cables of cable-stayed bridge [J]. Journal of Wind Engineering and Industrial Aerodynamics, 1992, 43(1): 2011.

    [3]顧明,劉慈軍,羅國(guó)強(qiáng),等. 斜拉橋拉索的風(fēng)(雨)激振及控制 [J]. 上海力學(xué),1998, 19(4): 283.

    GU Ming, LIU Cijun, LUO Guoqiang,etal. Rain-wind induced vibration of cables on cable-stayed bridges and its control [J]. Shanghai Journal of Mechanics, 1998, 19(4): 283.

    [4]Ni Y Q, Wang X Y, Chen Z Q. Field observations of rain-wind-induced cable vibration in cable-stayed Dongting Lake Bridge [J]. Journal of Wind Engineering and Industrial Aerodynamics, 2007, 95(5): 303.

    [5]Yamaguchi H. Analytical study on growth mechanism of rain vibration of cable [J]. Journal of Wind Engineering and Industrial Aerodynamics, 1990, 33(1): 73.

    [6]Gu M, Du X Q. Experimental investigation of rain-wind-induced vibration of cables in cable-stayed bridges and its mitigation [J]. Journal of Wind Engineering and Industrial Aerodynamics, 2005, 93(1): 79.

    [7]Matsumoto M, Yagi T, Saka S,etal. Steady wind force coefficients of inclined stay cables with water rivulet and their application to aerodynamics [J]. Wind and Structures, 2005, 8(2): 107.

    [8]李永樂(lè),徐幼麟,沈其民,等. 斜拉橋拉索風(fēng)-雨致振動(dòng)(I): 機(jī)理分析 [J]. 西南交通大學(xué)學(xué)報(bào),2011,46(4): 529.

    LI Yongle, XU Youlin, SHEN Qimin,etal. Rain-wind-induced vibration of cables in cable-stayed bridge (I): Mechanism analysis [J]. Journal of Southwest Jiaotong University, 2011, 46(4):529.

    [9]DU Xiaoqing, GU Ming, CHEN Suren. Aerodynamic characteristics of an inclined and yawed circular cylinder with artificial rivulet [J]. Journal of Fluids and Structures, 2013, 43: 64.

    [10]Lemaitre C. Thin water film around a cable subject to wind [J]. Journal of Wind Engineering and Industrial Aerodynamics, 2007, 95(9): 1259.

    [11]許林汕,趙林,葛耀君. 發(fā)生風(fēng)雨激振拉索表面水膜形態(tài)數(shù)值分析和實(shí)測(cè)初探 [J]. 空氣動(dòng)力學(xué)學(xué)報(bào),2011, 29(1):91.

    XU Linshan, ZHAO Lin, GE Yaojun. Numerical analysis and real-time measurement of water-film on rain-wind induced vibration cable [J]. Acta Aerodynamic Sinica, 2011, 29(1):91.

    [12]Bi J H, Lu P, Wang J,etal. Numerical simulation and analysis of the effects of water-film morphological changes on the aerodynamic lift of stay cables [J]. Journal of Fluids and Structures, 2014, 48: 376.

    [13]杜曉慶, 何平, 顧明. 斜拉索上水線氣動(dòng)性能的大渦模擬研究[J]. 土木工程學(xué)報(bào), 2014, 47(10): 12.

    DU Xiaoqing, HE Ping, GU Ming. Investigation on aerodynamic characteristics of upper rivulet on stay cable using large eddy simulation [J]. China Civil Engineering Journal, 2014, 47(10): 12.

    [14]Zdravkovich M M. Flow around circular cylinders [M]. Oxford : Oxford Science Publications, 1997.

    [15]Nishimura H, Taniike Y. Aerodynamic characteristics of fluctuating forces on a circular cylinder [J]. Journal of Wind Engineering and Industrial Aerodynamics, 2001, 89(7): 713.

    [16]Breuer M. A challenging test case for large eddy simulation: High Reynolds number circular cylinder flow [J]. International Journal of Heat and Fluid Flow, 2000, 21(5): 648.

    [17]Cantwell B, Coles D. An exprimental study of entrainment and transport in the turbulent near wake of a circular cylinder [J]. Journal of Fluids and Mechanics, 1983, 136: 321.

    [18]Behara S, Mittal S. Transition of the boundary layer on a circular cylinder in the presence of a trip [J]. Journal of Fluids and Structures, 2011, 27: 702.

    Large Eddy Simulation of Flow Around Stay Cable with Upper Rivulet

    DU Xiaoqing1, LI Junjun1, GU Ming2, HE Ping1

    (1. Department of Civil Engineering, Shanghai University, Shanghai 200072, China; 2. State Key Laboratory for Disaster Reduction in Civil Engineering, Tongji University, Shanghai 200092, China)

    To understand the relationship between the flow field characteristics and aerodynamics of the cable with the upper rivulet, the flow around and aerodynamic characteristics of a cable model with an upper rivulet were investigated using large eddy simulation (LES) in the subcritical regime. The results show that there exist three distinct flow regimes with the rivulet located at different positions. In the range ofθu=0°~50°, the rivulet has little influence on the cable. In the range ofθu=50°~70°, with the interferences of the upper rivulet, the upper-side shear layer of the cable which is detached from the rivulet can reattach on the cable surface and form a single separation bubble at the upper side of the cable, which induces an asymmetrical pressure distribution on the cable surface. In the range ofθu=70°~90°, the upper-side shear layer detaches at the rivulet, which causes a lower negative pressure on the upper-side surface compared to those of the lower-side of the cable. It is concluded that the appearance and disappearance of the single separation bubble are the underlying reasons for the occurrence of the rain-wind-induced vibration (RWIV).

    stay cable; rain-wind-induced vibration; large eddy simulation; upper rivulet; separation bubble

    2015-11-19

    上海市自然科學(xué)基金(14ZR1416000);上海市教委科研創(chuàng)新項(xiàng)目(14YZ004);國(guó)家自然科學(xué)基金(51578330)

    杜曉慶(1973—),男,副教授,工學(xué)博士,主要研究方向?yàn)榻Y(jié)構(gòu)風(fēng)工程. E-mail: dxq@shu.edu.cn

    U448.27

    A

    猜你喜歡
    水線升力風(fēng)壓
    高速列車車頂–升力翼組合體氣動(dòng)特性
    無(wú)人機(jī)升力測(cè)試裝置設(shè)計(jì)及誤差因素分析
    基于自適應(yīng)偽譜法的升力式飛行器火星進(jìn)入段快速軌跡優(yōu)化
    低風(fēng)壓架空導(dǎo)線的風(fēng)洞試驗(yàn)
    電線電纜(2017年5期)2017-10-18 00:52:03
    水線光纜
    水線光纜
    規(guī)?;怆u養(yǎng)殖場(chǎng)水線管理
    低風(fēng)壓導(dǎo)線的發(fā)展和測(cè)試
    電線電纜(2016年5期)2016-02-27 09:02:06
    升力式再入飛行器體襟翼姿態(tài)控制方法
    高層建筑風(fēng)荷載有哪些計(jì)算要求
    重慶建筑(2014年10期)2014-03-27 02:20:07
    亚洲精品自拍成人| 久久久久久久大尺度免费视频| 中文字幕免费在线视频6| 可以在线观看毛片的网站| 日韩,欧美,国产一区二区三区| 日韩人妻高清精品专区| 亚洲自偷自拍三级| 午夜精品在线福利| 欧美成人午夜免费资源| 内射极品少妇av片p| 国产av码专区亚洲av| 亚洲自拍偷在线| 亚洲国产欧美在线一区| 免费观看无遮挡的男女| 日本黄色片子视频| 成人综合一区亚洲| 亚洲色图av天堂| 日日啪夜夜撸| 国产精品爽爽va在线观看网站| 男人和女人高潮做爰伦理| 亚洲国产精品国产精品| 日韩在线高清观看一区二区三区| 久久久亚洲精品成人影院| 91av网一区二区| 日韩av在线免费看完整版不卡| 国产一区二区三区av在线| 亚洲真实伦在线观看| 日韩一本色道免费dvd| 国产女主播在线喷水免费视频网站 | 我的女老师完整版在线观看| 国产 亚洲一区二区三区 | 亚洲精品乱码久久久久久按摩| 身体一侧抽搐| 国产午夜精品一二区理论片| 在线观看人妻少妇| 欧美一级a爱片免费观看看| 国产亚洲av片在线观看秒播厂 | 日日干狠狠操夜夜爽| 男女国产视频网站| 免费少妇av软件| 91精品国产九色| 久久久久九九精品影院| .国产精品久久| 国产探花极品一区二区| 国产探花极品一区二区| 国产白丝娇喘喷水9色精品| 高清午夜精品一区二区三区| 成人综合一区亚洲| 国产v大片淫在线免费观看| 五月天丁香电影| 毛片女人毛片| 一级二级三级毛片免费看| 亚洲精品乱码久久久久久按摩| 国产永久视频网站| 又大又黄又爽视频免费| 免费看a级黄色片| freevideosex欧美| 成人毛片60女人毛片免费| 国产乱人视频| 欧美日韩一区二区视频在线观看视频在线 | 尤物成人国产欧美一区二区三区| 成年av动漫网址| 男人狂女人下面高潮的视频| 国产老妇伦熟女老妇高清| 久久6这里有精品| 国国产精品蜜臀av免费| 99久久精品热视频| 精品国产三级普通话版| 久久久久国产网址| 2021少妇久久久久久久久久久| 亚洲av电影在线观看一区二区三区 | 亚洲性久久影院| 免费看不卡的av| 人妻制服诱惑在线中文字幕| 国产成人精品福利久久| 欧美高清成人免费视频www| 少妇的逼水好多| 高清视频免费观看一区二区 | 国产午夜精品一二区理论片| 免费人成在线观看视频色| 男人狂女人下面高潮的视频| 亚洲国产成人一精品久久久| 日本三级黄在线观看| 联通29元200g的流量卡| 天堂av国产一区二区熟女人妻| 亚洲三级黄色毛片| 免费少妇av软件| 成人亚洲精品一区在线观看 | 91狼人影院| 精品国产一区二区三区久久久樱花 | 国产精品三级大全| 欧美xxxx黑人xx丫x性爽| 韩国高清视频一区二区三区| a级一级毛片免费在线观看| 国产爱豆传媒在线观看| 国产一区亚洲一区在线观看| 高清av免费在线| 2021少妇久久久久久久久久久| 中文字幕免费在线视频6| 大又大粗又爽又黄少妇毛片口| 国产免费福利视频在线观看| 久久久久久久久久黄片| 一级片'在线观看视频| 国产老妇伦熟女老妇高清| 久久精品国产亚洲av天美| 亚洲人成网站在线播| 亚洲四区av| 国产女主播在线喷水免费视频网站 | 成人无遮挡网站| 午夜福利在线在线| 夜夜看夜夜爽夜夜摸| 我的老师免费观看完整版| 亚洲精品,欧美精品| 成人毛片a级毛片在线播放| 亚洲国产精品sss在线观看| 日本午夜av视频| 日韩伦理黄色片| 我要看日韩黄色一级片| or卡值多少钱| 97超视频在线观看视频| 国产 一区精品| 成人毛片a级毛片在线播放| 国产男女超爽视频在线观看| 伊人久久精品亚洲午夜| 日韩视频在线欧美| 欧美变态另类bdsm刘玥| 欧美变态另类bdsm刘玥| 亚洲精华国产精华液的使用体验| 国产成人一区二区在线| 午夜久久久久精精品| 免费看日本二区| 欧美xxxx黑人xx丫x性爽| 青青草视频在线视频观看| 久久久久久久久久成人| 欧美xxxx性猛交bbbb| 国产成人a∨麻豆精品| 深爱激情五月婷婷| 国产精品女同一区二区软件| 一级毛片电影观看| 国产精品日韩av在线免费观看| 色哟哟·www| 久久人人爽人人爽人人片va| 国产精品人妻久久久影院| 有码 亚洲区| 久久精品久久精品一区二区三区| 啦啦啦中文免费视频观看日本| 精品久久久久久久久av| 床上黄色一级片| 乱人视频在线观看| 久久精品久久精品一区二区三区| 能在线免费看毛片的网站| 简卡轻食公司| 精品久久久久久久久久久久久| 欧美变态另类bdsm刘玥| 日韩,欧美,国产一区二区三区| 国产极品天堂在线| 搡老乐熟女国产| 精品人妻熟女av久视频| 少妇的逼水好多| 亚洲熟女精品中文字幕| 一个人免费在线观看电影| 日本av手机在线免费观看| 一区二区三区乱码不卡18| 国产日韩欧美在线精品| 永久免费av网站大全| 91久久精品电影网| 欧美精品国产亚洲| 性插视频无遮挡在线免费观看| 美女高潮的动态| 黄色配什么色好看| 久久99热这里只有精品18| 一级二级三级毛片免费看| 久久人人爽人人爽人人片va| 好男人在线观看高清免费视频| 天堂av国产一区二区熟女人妻| 久久99蜜桃精品久久| 国产亚洲最大av| 在现免费观看毛片| 熟妇人妻久久中文字幕3abv| 免费播放大片免费观看视频在线观看| 亚洲熟妇中文字幕五十中出| 在线 av 中文字幕| 国产欧美日韩精品一区二区| 免费观看a级毛片全部| 我要看日韩黄色一级片| 中文字幕免费在线视频6| 人妻制服诱惑在线中文字幕| 国产三级在线视频| 午夜福利在线观看吧| 99久久人妻综合| 80岁老熟妇乱子伦牲交| 少妇丰满av| 嫩草影院新地址| 青春草国产在线视频| 黄色配什么色好看| 久久99热这里只有精品18| 一个人观看的视频www高清免费观看| 午夜福利成人在线免费观看| 色网站视频免费| 亚洲,欧美,日韩| 日本一本二区三区精品| 在线a可以看的网站| 精品久久国产蜜桃| 亚洲一区高清亚洲精品| 大香蕉久久网| 国产成人91sexporn| 青青草视频在线视频观看| 成人av在线播放网站| 国产成人freesex在线| 久久99蜜桃精品久久| 欧美性感艳星| 内地一区二区视频在线| 大话2 男鬼变身卡| videossex国产| 国产亚洲av片在线观看秒播厂 | 免费看不卡的av| 久久国产乱子免费精品| 丰满少妇做爰视频| 日韩在线高清观看一区二区三区| 亚洲av国产av综合av卡| 亚洲精品久久午夜乱码| 国产日韩欧美在线精品| 免费高清在线观看视频在线观看| 青青草视频在线视频观看| 国产一区有黄有色的免费视频 | 成人漫画全彩无遮挡| 韩国av在线不卡| 久久这里有精品视频免费| 国产白丝娇喘喷水9色精品| 国产精品一及| 久久久久精品性色| 久99久视频精品免费| 久久久久国产网址| 亚洲精品色激情综合| 一本久久精品| 久久久久久久久久成人| 日日撸夜夜添| 最近的中文字幕免费完整| 国产成人精品福利久久| 99久久精品热视频| 成人亚洲精品av一区二区| 精品国产露脸久久av麻豆 | 大片免费播放器 马上看| 又爽又黄无遮挡网站| 国产黄片美女视频| 国产成年人精品一区二区| 麻豆成人午夜福利视频| 在线a可以看的网站| 亚州av有码| 亚洲av电影不卡..在线观看| 永久网站在线| av在线蜜桃| 日本与韩国留学比较| 最近最新中文字幕大全电影3| 国产成人精品一,二区| 欧美+日韩+精品| 三级男女做爰猛烈吃奶摸视频| 国产 一区精品| 国产精品久久久久久精品电影| 色哟哟·www| 亚洲精品国产av成人精品| 天堂俺去俺来也www色官网 | 亚洲国产最新在线播放| 午夜福利视频1000在线观看| 极品少妇高潮喷水抽搐| 国产白丝娇喘喷水9色精品| 国产真实伦视频高清在线观看| 成人午夜高清在线视频| 国产精品一二三区在线看| 国产欧美日韩精品一区二区| 国产av在哪里看| 一级片'在线观看视频| 日日干狠狠操夜夜爽| 日本黄色片子视频| 国产激情偷乱视频一区二区| 亚洲欧美日韩东京热| 91午夜精品亚洲一区二区三区| www.色视频.com| 久久6这里有精品| 天堂俺去俺来也www色官网 | 免费看av在线观看网站| 秋霞伦理黄片| 精华霜和精华液先用哪个| 成人午夜高清在线视频| 国产精品三级大全| 国产v大片淫在线免费观看| 亚洲国产精品国产精品| av国产免费在线观看| 一二三四中文在线观看免费高清| 亚洲精品一区蜜桃| av在线观看视频网站免费| 一夜夜www| 狠狠精品人妻久久久久久综合| 91午夜精品亚洲一区二区三区| 69人妻影院| 伊人久久国产一区二区| 国产有黄有色有爽视频| 淫秽高清视频在线观看| 中文字幕制服av| 日产精品乱码卡一卡2卡三| 卡戴珊不雅视频在线播放| 久久精品国产鲁丝片午夜精品| 亚洲成人精品中文字幕电影| 欧美xxxx黑人xx丫x性爽| 尤物成人国产欧美一区二区三区| 亚洲国产av新网站| 亚洲人与动物交配视频| 亚洲精华国产精华液的使用体验| 亚洲内射少妇av| 免费看日本二区| 国产一区二区三区综合在线观看 | 人妻夜夜爽99麻豆av| 色哟哟·www| 国产成人精品一,二区| 国产精品av视频在线免费观看| 国产成人aa在线观看| 日本-黄色视频高清免费观看| 一级a做视频免费观看| 蜜桃亚洲精品一区二区三区| 蜜臀久久99精品久久宅男| 欧美zozozo另类| 国产视频内射| 男女啪啪激烈高潮av片| 国产精品久久久久久av不卡| 免费少妇av软件| 亚洲乱码一区二区免费版| 嫩草影院入口| 成年女人在线观看亚洲视频 | 精品国产一区二区三区久久久樱花 | 精品亚洲乱码少妇综合久久| 大话2 男鬼变身卡| 亚洲国产最新在线播放| 国产白丝娇喘喷水9色精品| 国产亚洲午夜精品一区二区久久 | 天天躁日日操中文字幕| 久久精品熟女亚洲av麻豆精品 | 高清日韩中文字幕在线| 久久韩国三级中文字幕| 国产午夜精品一二区理论片| 高清日韩中文字幕在线| 色尼玛亚洲综合影院| av天堂中文字幕网| 一级a做视频免费观看| 成人综合一区亚洲| 五月伊人婷婷丁香| 少妇人妻精品综合一区二区| 在线a可以看的网站| 国产男人的电影天堂91| 欧美日韩亚洲高清精品| 欧美精品一区二区大全| 久久久久网色| av天堂中文字幕网| 搡老乐熟女国产| 国产又色又爽无遮挡免| 亚洲综合色惰| 精品亚洲乱码少妇综合久久| 亚洲自拍偷在线| 欧美3d第一页| 麻豆久久精品国产亚洲av| 亚洲精品久久午夜乱码| 联通29元200g的流量卡| 少妇被粗大猛烈的视频| 干丝袜人妻中文字幕| 亚洲第一区二区三区不卡| 一级二级三级毛片免费看| 中国美白少妇内射xxxbb| 国产成人91sexporn| 亚洲美女搞黄在线观看| 亚洲精品,欧美精品| av一本久久久久| 免费黄频网站在线观看国产| 只有这里有精品99| 美女高潮的动态| 大话2 男鬼变身卡| 午夜爱爱视频在线播放| 美女脱内裤让男人舔精品视频| 能在线免费看毛片的网站| 欧美高清成人免费视频www| 国产高清国产精品国产三级 | 免费高清在线观看视频在线观看| 国产精品久久视频播放| 国产有黄有色有爽视频| 国产黄色小视频在线观看| 91精品一卡2卡3卡4卡| 菩萨蛮人人尽说江南好唐韦庄| 三级国产精品片| 亚洲欧美精品专区久久| 国产片特级美女逼逼视频| 免费黄频网站在线观看国产| 亚州av有码| 国产精品伦人一区二区| 熟女电影av网| 国内精品宾馆在线| 中国美白少妇内射xxxbb| 国产欧美另类精品又又久久亚洲欧美| 国产精品嫩草影院av在线观看| 十八禁国产超污无遮挡网站| 精品久久久久久久久亚洲| 国产精品美女特级片免费视频播放器| 亚洲人成网站在线观看播放| 欧美性猛交╳xxx乱大交人| 最近最新中文字幕大全电影3| av免费在线看不卡| 亚洲av电影不卡..在线观看| 自拍偷自拍亚洲精品老妇| 免费看av在线观看网站| 国产永久视频网站| 久久人人爽人人片av| 熟女电影av网| 日韩成人av中文字幕在线观看| 精品少妇黑人巨大在线播放| 久久精品综合一区二区三区| 国产精品麻豆人妻色哟哟久久 | 色吧在线观看| 视频中文字幕在线观看| 国产亚洲最大av| 人妻系列 视频| av专区在线播放| 搞女人的毛片| 一级黄片播放器| 久久久久久久国产电影| 精品少妇黑人巨大在线播放| 最新中文字幕久久久久| av在线天堂中文字幕| 成人鲁丝片一二三区免费| 一级毛片aaaaaa免费看小| 大又大粗又爽又黄少妇毛片口| 国产午夜精品久久久久久一区二区三区| 国产精品爽爽va在线观看网站| 秋霞在线观看毛片| 丝袜喷水一区| 亚洲精品亚洲一区二区| 男女边吃奶边做爰视频| 久久精品夜色国产| 中文字幕av在线有码专区| 内射极品少妇av片p| 日韩人妻高清精品专区| 国产成人91sexporn| 一区二区三区四区激情视频| 成人欧美大片| 亚洲成人中文字幕在线播放| videossex国产| 熟妇人妻久久中文字幕3abv| 国产伦精品一区二区三区视频9| 97超视频在线观看视频| 久久精品夜色国产| 日本黄大片高清| 青春草亚洲视频在线观看| 久热久热在线精品观看| 亚洲性久久影院| 亚洲伊人久久精品综合| 高清欧美精品videossex| 99久国产av精品| 全区人妻精品视频| 欧美xxxx黑人xx丫x性爽| 日日摸夜夜添夜夜添av毛片| 国产成人a∨麻豆精品| 中文资源天堂在线| 国产免费视频播放在线视频 | 午夜福利在线观看免费完整高清在| 亚洲aⅴ乱码一区二区在线播放| 欧美高清成人免费视频www| 一边亲一边摸免费视频| 亚洲,欧美,日韩| 欧美性猛交╳xxx乱大交人| 色综合色国产| 国产精品久久久久久精品电影| 青春草视频在线免费观看| 亚洲av.av天堂| 在线观看人妻少妇| 狂野欧美激情性xxxx在线观看| 亚洲在线观看片| 我要看日韩黄色一级片| 联通29元200g的流量卡| 天美传媒精品一区二区| 国产综合懂色| 狂野欧美激情性xxxx在线观看| 日韩一区二区视频免费看| 观看美女的网站| 久久久久久久久中文| 男人舔女人下体高潮全视频| 亚洲av不卡在线观看| 国产黄色免费在线视频| 高清毛片免费看| 看非洲黑人一级黄片| 乱码一卡2卡4卡精品| 日韩不卡一区二区三区视频在线| 日韩视频在线欧美| 亚洲av成人av| 欧美+日韩+精品| av在线天堂中文字幕| 白带黄色成豆腐渣| av国产久精品久网站免费入址| 成人欧美大片| 成人午夜高清在线视频| 十八禁网站网址无遮挡 | 中文乱码字字幕精品一区二区三区 | 免费无遮挡裸体视频| 国产精品国产三级国产av玫瑰| 看黄色毛片网站| 久久久久久久午夜电影| 国产一区二区三区综合在线观看 | 亚洲国产精品sss在线观看| 午夜精品国产一区二区电影 | 亚洲色图av天堂| 乱人视频在线观看| 国产v大片淫在线免费观看| 国产午夜精品久久久久久一区二区三区| 精品亚洲乱码少妇综合久久| 国产亚洲午夜精品一区二区久久 | 狠狠精品人妻久久久久久综合| 岛国毛片在线播放| 亚洲国产精品成人综合色| 国产免费一级a男人的天堂| 欧美成人a在线观看| 国内少妇人妻偷人精品xxx网站| 大香蕉97超碰在线| 我的老师免费观看完整版| 一本久久精品| 日产精品乱码卡一卡2卡三| 欧美丝袜亚洲另类| 九九在线视频观看精品| 91在线精品国自产拍蜜月| 国内精品美女久久久久久| 午夜福利在线观看吧| 亚洲激情五月婷婷啪啪| 国产精品熟女久久久久浪| 深夜a级毛片| 女的被弄到高潮叫床怎么办| 久久久国产一区二区| 老女人水多毛片| 偷拍熟女少妇极品色| 直男gayav资源| 亚洲精品国产成人久久av| 国产精品人妻久久久久久| 午夜福利视频精品| 一个人看的www免费观看视频| 国产成人aa在线观看| 午夜日本视频在线| 别揉我奶头 嗯啊视频| 99久久中文字幕三级久久日本| 人人妻人人看人人澡| 国产高清国产精品国产三级 | 一二三四中文在线观看免费高清| 可以在线观看毛片的网站| 色播亚洲综合网| av天堂中文字幕网| 三级国产精品片| 亚洲激情五月婷婷啪啪| 日日啪夜夜爽| 七月丁香在线播放| av线在线观看网站| 国产一区亚洲一区在线观看| 秋霞伦理黄片| 色综合亚洲欧美另类图片| 亚洲av在线观看美女高潮| 三级经典国产精品| 一夜夜www| 国产高潮美女av| 国产精品熟女久久久久浪| 欧美日韩精品成人综合77777| 亚洲成人一二三区av| 亚洲av免费在线观看| 欧美成人午夜免费资源| 国产成人福利小说| 成人高潮视频无遮挡免费网站| 欧美xxxx黑人xx丫x性爽| 自拍偷自拍亚洲精品老妇| 欧美3d第一页| 2018国产大陆天天弄谢| 久久精品熟女亚洲av麻豆精品 | 欧美xxxx黑人xx丫x性爽| 插逼视频在线观看| 亚洲精品456在线播放app| 又大又黄又爽视频免费| 欧美高清成人免费视频www| 日本免费a在线| 大话2 男鬼变身卡| 97精品久久久久久久久久精品| 日韩精品有码人妻一区| 国产亚洲av片在线观看秒播厂 | 国产av不卡久久| 亚洲av不卡在线观看| 1000部很黄的大片| 免费观看无遮挡的男女| 在线a可以看的网站| 干丝袜人妻中文字幕| 97人妻精品一区二区三区麻豆| 国产欧美日韩精品一区二区| 午夜激情久久久久久久| 欧美激情国产日韩精品一区| 日本猛色少妇xxxxx猛交久久| 赤兔流量卡办理| 国产亚洲91精品色在线| 亚洲国产成人一精品久久久| 日韩制服骚丝袜av| 男人狂女人下面高潮的视频| 一区二区三区高清视频在线| 亚洲激情五月婷婷啪啪| 欧美日韩精品成人综合77777| 国产精品一区www在线观看| 日本黄大片高清| 大片免费播放器 马上看| 国产黄色小视频在线观看| 91精品一卡2卡3卡4卡| 久久久久久国产a免费观看| 午夜精品一区二区三区免费看| 五月伊人婷婷丁香| 国产在视频线在精品| 黄片wwwwww| 建设人人有责人人尽责人人享有的 | 搞女人的毛片| xxx大片免费视频| 伊人久久国产一区二区| 又大又黄又爽视频免费|