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

    橋面水馬布置對(duì)大跨度懸索橋渦激振動(dòng)的影響

    2024-01-01 00:00:00周椿檳周光偉陳昌萍
    關(guān)鍵詞:橋梁

    摘要: 為研究橋面水馬布置對(duì)橋梁渦激振動(dòng)的影響,以某大跨度懸索橋?yàn)楣こ瘫尘?,基于?jì)算流體力學(xué)(CFD)方法,研究不同風(fēng)速、橋面水馬布置方案下的橋梁斷面受力情況及其周圍流場(chǎng)的變化規(guī)律。結(jié)果表明:在橋面上布置水馬能夠改變橋梁斷面周圍流場(chǎng)的壓力和速度分布,在橋梁上部和背風(fēng)側(cè)產(chǎn)生漩渦,可觀察到明顯的漩渦脫落現(xiàn)象,與無(wú)布置水馬的橋梁相比,其發(fā)生渦激振動(dòng)的可能性大大提升;水馬與橋面邊緣的距離越大,漩渦脫落現(xiàn)象越不明顯,橋梁發(fā)生渦激振動(dòng)的可能性越低。

    關(guān)鍵詞:橋梁; 水馬布置; 漩渦脫落; 渦激振動(dòng); 數(shù)值模擬

    中圖分類號(hào): U 446.1文獻(xiàn)標(biāo)志碼: A"" 文章編號(hào): 1000-5013(2024)06-0712-09

    Influence of Bridge Deck Barriers Arrangement on "Vortex-Induced Vibration of Long-Span Suspension Bridge

    ZHOU Chunbin1, ZHOU Guangwei2, CHEN Changping1,2,3

    (1. School of Architecture and Civil Engineering, Xiamen University, Xiamen 361000, China;

    2. Fujian Provincial Key Laboratory of Wind Disaster and Wind Engineering,Xiamen University of Technology, Xiamen 361024, China;

    3. Xiamen Ocean Vocational College, Xiamen 361100, China)

    Abstract: In order to study the influence of bridge deck barriers arrangement on bridge vortex-induced vibration, taking a long-span suspension bridge as the engineering background, based on the computational fluid dynamics (CFD) method, the stress situation of bridge section and the variation of surrounding flow field are studied under different wind speeds and bridge deck barriers arrangement. The results show that the arrangement of the barriers on the bridge deck can change the pressure and speed distribution of surrounding flow field of bridge section, and the vortex is generated in the upper part and the leeward side of the bridge, obvious vortex shedding phenomenon can be observed, and compared with the bridge without the barriers, the possibility of vortex-induced vibration increases greatly. The greater the distance between the barriers and the edge of the bridge deck, the less obvious the vortex shedding phenomenon, and the lower the possibility of vortex-induced vibration of the bridge.

    Keywords:bridge; barriers arrangement; vortex shedding; vortex-induced vibration; numerical simulation

    水馬作為橋梁臨時(shí)維護(hù)的圍擋措施,廣泛應(yīng)用于橋梁的檢測(cè)、維修等市政工程中。因市政工程的特點(diǎn),水馬往往沿橋梁跨徑采取緊密排列的形式進(jìn)行布置,但這樣的水馬布置改變了橋梁的氣動(dòng)外形,增加了橋梁產(chǎn)生渦激振動(dòng)的可能性。2020年5月5日,廣東省的虎門(mén)大橋進(jìn)行市政施工時(shí),因水馬的布置改變了橋梁的氣動(dòng)外形,進(jìn)而產(chǎn)生了較大振幅的渦激振動(dòng)。渦激振動(dòng)是一種在低風(fēng)速條件下,柔性橋梁因經(jīng)過(guò)其結(jié)構(gòu)表面的流體漩渦發(fā)生規(guī)律的周期性脫落而導(dǎo)致橋梁發(fā)生限幅振動(dòng)的一種風(fēng)致振動(dòng)現(xiàn)象[1]。雖然渦激振動(dòng)不會(huì)在短時(shí)間內(nèi)對(duì)橋梁結(jié)構(gòu)造成毀滅性的破壞,但幅度較大的渦激振動(dòng)會(huì)嚴(yán)重影響橋梁上行車、行人的舒適性,加速受力構(gòu)件的疲勞破壞,導(dǎo)致橋梁出現(xiàn)結(jié)構(gòu)損傷甚至斷裂[2]。橋梁渦激振動(dòng)的產(chǎn)生不僅與橋梁本身的質(zhì)量、剛度、阻尼有關(guān)[3-4],也與橋梁截面的氣動(dòng)外形密切相關(guān)[5]。

    葛耀君等[6]總結(jié)大跨度橋梁主梁渦激振動(dòng)的研究進(jìn)展,描述渦激振動(dòng)控制措施的發(fā)展趨勢(shì),列舉主流的橋梁氣動(dòng)控制措施,對(duì)未來(lái)的研究方向做出了一定判斷。華旭剛等[7]討論大跨度懸索橋多個(gè)豎向振動(dòng)模態(tài)分布的特性,發(fā)現(xiàn)僅以一階模態(tài)進(jìn)行渦振性能評(píng)估不夠全面。黃俐等[8]針對(duì)大寬高比橋梁斷面的渦振鎖定區(qū)間進(jìn)行數(shù)值模擬研究,總結(jié)典型卡門(mén)渦街和“魚(yú)尾擺動(dòng)”兩種尾流形態(tài),發(fā)現(xiàn)“魚(yú)尾擺動(dòng)”是介于卡門(mén)渦街現(xiàn)象和平穩(wěn)流場(chǎng)之間的流場(chǎng)形態(tài)。王維民等[9]對(duì)橋梁豎彎渦激振動(dòng)情況下的行車安全舒適性進(jìn)行研究。上述研究主要針對(duì)渦激振動(dòng)特性和采用永久性氣動(dòng)外形控制措施的影響,而針對(duì)橋面水馬這樣的臨時(shí)施工措施對(duì)橋梁氣動(dòng)外形的改變而導(dǎo)致的渦激振動(dòng)的研究較為缺乏。

    鑒于此,本文采用計(jì)算流體力學(xué)(CFD)方法[10],研究低風(fēng)速情況下,不同的水馬布置情況對(duì)橋梁斷面周圍流場(chǎng)的影響,探究橋梁渦激振動(dòng)的發(fā)生機(jī)制,提出橋梁施工期間橋面水馬的合理布置方法。

    1 工程概況

    1.1 工程背景

    某橋全長(zhǎng)1 108 m(230 m+648 m+230 m),為三跨懸索橋結(jié)構(gòu)體系。全橋由鋼箱梁、索塔及吊桿組成,屬于大跨度柔性懸索橋。

    橋面設(shè)置2%的雙向坡,橋面兩側(cè)為雙向6車道,與虎門(mén)大橋主梁斷面氣動(dòng)外形相近,主梁標(biāo)準(zhǔn)斷面的寬為36.6 m,高度為2.67 m。橋面上設(shè)有路緣石,路緣石寬度為240 mm,高度為140 mm,橋梁兩側(cè)及中央?yún)^(qū)域設(shè)有欄桿。主梁標(biāo)準(zhǔn)截面的布置,如圖1所示。

    1.2 橋梁動(dòng)力特性分析

    橋梁的渦激振動(dòng)與全橋的動(dòng)力特性密切相關(guān)[11]。采用ANSYSAPDL軟件建立橋梁有限元模型[12],主梁、索塔采用BEAM4單元進(jìn)行模擬。

    有限元模型中,鋼箱梁彈性模量為0.21 TPa,泊松比為0.3,索塔彈性模量為34.5 GPa,承臺(tái)彈性模量為30.0 GPa,橋梁主纜彈性模量為0.198 TPa。

    全橋有限元模型,如圖2所示。圖2中:X,Y,Z為空間坐標(biāo)系。

    經(jīng)計(jì)算,可得該橋梁低階豎彎及扭轉(zhuǎn)模態(tài)的動(dòng)力特性,如表1所示。表1中:f為頻率。

    橋梁模態(tài)振型圖,如圖3所示。

    2 數(shù)值模擬方法

    2.1 流體計(jì)算關(guān)系

    靜三分力系數(shù)是用于描述成橋主梁斷面風(fēng)荷載特性的無(wú)量綱參數(shù),主要包括阻力系數(shù)(CD)、升力系數(shù)(CL)、扭矩系數(shù)(CM),分別與因前后表面壓強(qiáng)差而產(chǎn)生的阻力(FD)、上下表面壓強(qiáng)差產(chǎn)生的升力(FL)、升阻力合力點(diǎn)與扭心不一致產(chǎn)生的扭矩(MT)有對(duì)應(yīng)關(guān)系[13]。即

    CD=2FDρU2D,(1)

    CL=2FLρU2B,(2)

    CM=2MTρU2B2。(3)

    式(1)~(3)中:ρ為流體密度,即空氣密度1.29 kg·m-3;B為橋梁斷面的投影寬度;D為橋梁斷面的投影高度;U為來(lái)風(fēng)速。

    2.2 計(jì)算模型及工況設(shè)置

    采用Model Designer軟件建立懸索橋主梁斷面簡(jiǎn)化縮尺模型[14],模型分為鋼箱梁截面、欄桿、橋面水馬3個(gè)部分。整體流場(chǎng)計(jì)算域尺寸(寬×長(zhǎng))設(shè)置為18.3 m×54.9 m,為使橋梁周圍流場(chǎng)的模擬結(jié)果更為準(zhǔn)確,加密網(wǎng)格流場(chǎng)區(qū)域尺寸(寬×長(zhǎng))設(shè)置為4 m×8 m,水馬布置共計(jì)4排,由橋梁迎風(fēng)側(cè)至背風(fēng)側(cè)分別編號(hào)為水馬1~4,水馬斷面寬度為0.5 m,高度為1.0 m,氣動(dòng)外形簡(jiǎn)化為矩形。各方案主梁斷面模型及計(jì)算域,如圖4所示。

    橋面兩側(cè)布置水馬的方案中,每個(gè)水馬與欄桿的距離為0.5 m;橋面中央布置水馬的方案中,水馬1,4分別放置于橋面1/4和3/4處,水馬2,3與欄桿的距離均為0.5 m。

    2.3 計(jì)算方法

    計(jì)算模型建立后,采用Fluent Meshing軟件進(jìn)行網(wǎng)格劃分,采用面網(wǎng)格剖分方法,網(wǎng)格形式為非結(jié)構(gòu)化三角形網(wǎng)格,外部流場(chǎng)單元尺寸設(shè)置為2 000 mm,加密網(wǎng)格區(qū)域流場(chǎng)單元尺寸設(shè)置為500 mm,為滿足湍流模型對(duì)壁面Y+值的要求,對(duì)欄桿等細(xì)小構(gòu)件部分及主梁做邊界層網(wǎng)格膨脹處理[15]。網(wǎng)格總節(jié)點(diǎn)數(shù)為72 719,總單元數(shù)為98 574。主梁斷面周圍流場(chǎng)網(wǎng)格劃分,如圖5所示。

    為更好地得到主梁斷面周圍流場(chǎng)的參數(shù),湍流模型選用SST k-ω,流場(chǎng)入口速度根據(jù)工況設(shè)置選取,湍流強(qiáng)度為0.5%。流場(chǎng)出口定義為壓力出口,上下邊界定義為對(duì)稱邊界。設(shè)置計(jì)算步數(shù)為10 000步,計(jì)算步長(zhǎng)為0.01 s(總計(jì)算區(qū)間時(shí)長(zhǎng)為100 s)。采用SIMPLE方法處理壓力與速度的耦合關(guān)系。在計(jì)算過(guò)程中,選取主梁斷面的阻力、升力、力矩、壓強(qiáng)和速度作為計(jì)算步長(zhǎng)監(jiān)控的參數(shù),輸出結(jié)果并進(jìn)行分析。

    渦激振動(dòng)帶有自激性質(zhì),但振動(dòng)的結(jié)構(gòu)反過(guò)來(lái)會(huì)對(duì)漩渦脫落形成某種反饋?zhàn)饔?,使渦振振幅受到限制。取流場(chǎng)發(fā)展和變化趨勢(shì)較為穩(wěn)定的時(shí)間作為研究區(qū)間,即來(lái)流后的50~100 s。

    渦激振動(dòng)往往發(fā)生在低風(fēng)速情況下,因漩渦脫落的頻率接近橋梁低階振型頻率而發(fā)生較大振幅的振動(dòng)。為探究漩渦脫落頻率與風(fēng)速之間的關(guān)系,找到橋梁斷面有可能發(fā)生渦激振動(dòng)的風(fēng)速區(qū)間,設(shè)置風(fēng)速區(qū)間為1~10 m·s-1,每次計(jì)算間隔為1 m·s-1。

    通過(guò)統(tǒng)計(jì)橋梁斷面的升力卓越頻率,可得漩渦脫落頻率-風(fēng)速關(guān)系圖,如圖6所示。

    圖6中:fs為漩渦脫落頻率;v為風(fēng)速。

    由圖6可知:發(fā)生渦激振動(dòng)可能性較大的風(fēng)速區(qū)間為4~7 m·s-1,故選取該風(fēng)速參數(shù)作為速度入口條件。

    為研究水馬布置對(duì)主梁斷面周圍流場(chǎng)的影響,探究橋梁發(fā)生渦激振動(dòng)的機(jī)制,共設(shè)置6種工況。工況設(shè)置,如表2所示。

    3 實(shí)驗(yàn)結(jié)果與分析

    3.1 靜三分力系數(shù)

    為研究不同風(fēng)速情況下,不同水馬布置方案對(duì)主梁斷面靜三分力系數(shù)的影響,采取CFD方法,選用0°風(fēng)攻角對(duì)主梁斷面周圍流場(chǎng)進(jìn)行數(shù)值模擬分析,可得各工況下的阻力、升力和扭矩時(shí)程圖。設(shè)置較長(zhǎng)的時(shí)間區(qū)間,使流場(chǎng)趨于穩(wěn)定后,得到橋梁斷面的阻力、升力和扭矩,分別根據(jù)式(1)~(3)計(jì)算對(duì)應(yīng)的阻力系數(shù)、升力系數(shù)和扭矩系數(shù)。靜三分力系數(shù),如表3所示。

    由表3可知:當(dāng)橋梁斷面的氣動(dòng)外形發(fā)生改變,如在橋面的兩側(cè)或中央布置水馬時(shí),在風(fēng)速相同的情況下,靜三分力系數(shù)發(fā)生較大的改變;阻力系數(shù)從大到小排序?yàn)闃蛎鎯蓚?cè)布置水馬、橋面中央布置水馬、無(wú)水馬布置,且阻力方向均為正向;升力系數(shù)絕對(duì)值從大到小排序?yàn)闊o(wú)水馬布置、橋面兩側(cè)布置水馬、橋面中央布置水馬,且升力方向均為負(fù)向;扭矩系數(shù)從大到小排序?yàn)闃蛎嬷醒氩贾盟R、無(wú)水馬布置、橋面兩側(cè)水馬布置。因此,可推斷橋梁斷面的靜三分力系數(shù)與水馬的布置有較強(qiáng)的關(guān)聯(lián)性。

    3.2 不同工況下流場(chǎng)壓力及速度分析

    漩渦脫落導(dǎo)致的橋梁渦激振動(dòng)是一個(gè)連續(xù)的過(guò)程,因此,在數(shù)值模擬中采用瞬時(shí)計(jì)算迭代的方法計(jì)算每一步長(zhǎng)的流場(chǎng)情況,設(shè)置合理的時(shí)間步數(shù),最終疊加每一瞬時(shí)的流場(chǎng)情況,可得漩渦脫落持續(xù)時(shí)間內(nèi)流場(chǎng)的持續(xù)變化過(guò)程。為了使計(jì)算結(jié)果更加精確,將時(shí)間步長(zhǎng)設(shè)置為0.01 s。

    不同工況下流場(chǎng)模擬結(jié)果圖,如圖7所示。

    由圖7可知以下6個(gè)結(jié)論。

    1) 當(dāng)風(fēng)荷載以一定速度流經(jīng)主梁斷面時(shí),會(huì)在主梁鋼箱梁迎風(fēng)側(cè)形成以風(fēng)嘴為核心,逐漸向外擴(kuò)散的正高壓區(qū),也會(huì)在鋼箱梁截面底部的兩個(gè)端點(diǎn)處形成范圍較小的圓形負(fù)壓區(qū)。

    2) 流場(chǎng)速度沿鋼箱梁邊界梯度分布,越靠近鋼箱梁表面速度越低,并在背風(fēng)側(cè)形成相同速度的漩渦,風(fēng)速越大,主梁截面背風(fēng)側(cè)的渦量越大(圖7(b),(h))。

    3) 工況2,5(橋面兩側(cè)布置水馬)發(fā)生明顯的漩渦脫落現(xiàn)象,主梁背風(fēng)側(cè)可見(jiàn)明顯的漩渦交替脫落,鋼箱梁上端也因水馬的阻隔作用出現(xiàn)較為明顯的正負(fù)壓交替分布;工況3,6(橋面中央布置水馬)也可于主梁背風(fēng)側(cè)觀察到漩渦脫落現(xiàn)象,但相較于橋面兩側(cè)布置水馬的工況并不明顯。

    4) 氣流經(jīng)過(guò)主梁截面風(fēng)嘴時(shí)幾乎沿著風(fēng)嘴壁面運(yùn)動(dòng),氣流與壁面并未發(fā)生大角度碰撞,并在風(fēng)嘴的分流下形成上、下兩股氣流,下端的氣流因壁面邊界層的粘性,始終沿著主梁截面下部行進(jìn)至背風(fēng)側(cè);上端氣流受到了欄桿的阻擋,但由于欄桿鏤空率較大,氣流從欄桿截面間隔中穿過(guò),沿著主梁上端截面流動(dòng),與下端氣流匯聚于背風(fēng)側(cè),形成上、下兩個(gè)漩渦,二者共同作用下產(chǎn)生了一個(gè)較大的負(fù)壓區(qū)(圖7(a),(b),(g),(h))。

    5) 在橋面兩側(cè)布置水馬,氣流經(jīng)過(guò)迎風(fēng)側(cè)風(fēng)嘴被分為上、下兩股氣流后,下端氣流仍沿主梁斷面下部行進(jìn)至背風(fēng)側(cè),上端的氣流受到迎風(fēng)側(cè)水馬1的阻隔,使上部氣流沿著鋼箱梁上壁面先分離,再附著,進(jìn)而在1號(hào)水馬后側(cè)形成較大尺寸的漩渦(圖7(c),(d),(i),(j))。上述現(xiàn)象同樣發(fā)生在氣流遭遇水馬2~4的阻擋時(shí),并分別形成了尺寸大小不一的漩渦。在上、下端氣流在背風(fēng)側(cè)尾流區(qū)匯合后,因分離后的高速氣流與低速漩渦之間的氣壓差交替變化,在主梁斷面背風(fēng)側(cè)產(chǎn)生了較為明顯的漩渦脫落現(xiàn)象。從相同的時(shí)間步長(zhǎng)中脫落后的漩渦之間的距離大致相等可判斷漩渦脫落有一個(gè)較為固定的卓越頻率。

    6) 在橋面中央布置水馬,下端氣流沿著主梁斷面下部行進(jìn)至背風(fēng)側(cè)。上端氣流首先經(jīng)過(guò)欄桿斷面分流,受到水馬1阻隔后,氣流分離附著,在水馬1前后各形成一個(gè)漩渦,且上端氣流在經(jīng)過(guò)水馬2~4時(shí)也發(fā)生相同現(xiàn)象,前、后各產(chǎn)生1個(gè)大小相近的漩渦,并沿著鋼箱梁上部截面連接形成扁平狹長(zhǎng)的負(fù)壓區(qū)。上端氣流和下端氣流各自產(chǎn)生一個(gè)大小不等的漩渦,并在分離氣流的作用下發(fā)生較不顯著的漩渦脫落現(xiàn)象(圖7(e),(f),(k),(l))。

    橋面兩側(cè)布置水馬的情況能夠發(fā)生較為顯著的漩渦脫落現(xiàn)象,但工況2中,流場(chǎng)升力、扭矩卓越頻率分別為0.110,0.139 Hz,與主梁一階對(duì)稱豎彎、主梁一階正對(duì)稱扭轉(zhuǎn)模態(tài)頻率相差較大,不會(huì)導(dǎo)致橋梁發(fā)生渦激振動(dòng)。因此,取工況5為詳細(xì)研究對(duì)象,若工況5的漩渦脫落頻率與橋梁主要振型頻率相近,則會(huì)發(fā)生共振現(xiàn)象,橋梁的渦激振動(dòng)振幅也將大大增加,嚴(yán)重影響行車、行人的舒適性,甚至發(fā)生構(gòu)件的疲勞破壞。

    3.3 橋面水馬布置對(duì)橋梁斷面渦激振動(dòng)的影響

    為了進(jìn)一步研究該橋梁斷面在布置了橋面水馬后是否會(huì)產(chǎn)生渦激振動(dòng),通過(guò)對(duì)比漩渦脫落頻率與橋梁模態(tài)頻率的方法來(lái)判斷是否有共振現(xiàn)象的產(chǎn)生。

    為研究橋面水馬布置對(duì)于橋梁斷面渦激振動(dòng)的影響,選取具有代表性的工況5,6作為研究對(duì)象。漩渦脫落頻率往往與升力、扭矩的變化頻率相等,故選取升力、扭矩時(shí)程圖,通過(guò)快速傅里葉變換(FFT)得到頻譜圖與漩渦脫落的卓越頻率,進(jìn)而判斷橋梁斷面是否會(huì)發(fā)生渦激振動(dòng)。

    為了探究不同的水馬布置對(duì)升力、扭矩的影響,分別選取工況5,6的升力、扭矩時(shí)程圖為研究對(duì)象。

    當(dāng)氣流并未完全流經(jīng)橋梁斷面時(shí),橋梁斷面的升力、扭矩時(shí)程圖并不規(guī)則且偏差較大。模擬結(jié)果顯示,氣流經(jīng)入口進(jìn)入流場(chǎng)后50 s,橋梁升力、扭矩時(shí)程曲線的變化基本趨于穩(wěn)定。因此,選取來(lái)流后50~100 s作為時(shí)程圖橫坐標(biāo)。

    工況5,6的升力、扭矩時(shí)程圖及頻譜圖,如圖8,9所示。

    對(duì)比圖8,9可知:在橋梁斷面上的不同位置布置水馬,在流場(chǎng)的作用下,在主梁斷面上產(chǎn)生的升力和扭矩的大小、頻率、特征均有明顯不同。

    由圖8可知:在風(fēng)速為7 m·s-1,橋面兩側(cè)布置水馬的情況下,主梁斷面的升力、扭矩時(shí)程圖呈明顯的周期性變化,推測(cè)是漩渦脫落的頻率固定導(dǎo)致的;主梁斷面升力在-99.37~-96.94 N波動(dòng),平均值為-98.01 N,方向豎直向下;升力的卓越頻率為0.180 Hz,根據(jù)橋梁的動(dòng)力特性分析,該卓越頻率與主梁二階反對(duì)稱豎彎振型頻率0.175 76 Hz接近,工況5在該模態(tài)下有可能產(chǎn)生橋梁的豎向渦激振動(dòng);工況5下主梁斷面扭矩在102.21~169.51 N·m周期性波動(dòng),平均值為135.45 N·m,方向?yàn)轫槙r(shí)針;扭矩的卓越頻率為0.231 Hz,與橋梁的一階正對(duì)稱扭轉(zhuǎn)模態(tài)振型頻率0.226 09 Hz較接近,工況5在該模態(tài)下有可能產(chǎn)生橋梁的扭轉(zhuǎn)渦激振動(dòng)。

    由圖9可知:在風(fēng)速為7 m·s-1,橋面中央布置水馬的情況下,升力、扭矩時(shí)程圖不再呈現(xiàn)類似工況5中明顯的固定頻率周期性變化,而是在發(fā)展一定時(shí)間后趨于平穩(wěn),但在早期仍有一定的波動(dòng);主梁斷面升力在-92.71~-72.18 N不規(guī)律波動(dòng),平均值為-78.54 N,方向豎直向下;升力的卓越頻率為0.06 Hz,該卓越頻率與橋梁各模態(tài)振型相差較大,推測(cè)工況6產(chǎn)生橋梁豎向渦激振動(dòng)的可能性較??;主梁斷面扭矩在264.32~410.28 N·m不規(guī)律波動(dòng),平均值為307.62 N·m,方向?yàn)轫槙r(shí)針;扭矩存在多個(gè)卓越頻率,分別為0.02,0.12 Hz,扭矩時(shí)程曲線的卓越頻率與橋梁扭轉(zhuǎn)模態(tài)的頻率差距較大,推測(cè)產(chǎn)生橋梁扭轉(zhuǎn)渦激振動(dòng)的可能性較小。

    由此可知:橋面中央布置水馬相較橋面兩側(cè)布置水馬的升力、扭矩卓越頻率與橋梁模態(tài)頻率相差較大,在該布置方案下橋梁產(chǎn)生渦激振動(dòng)的可能性較低。

    4 結(jié)論

    1) 截面尺寸較小的水馬依然能夠顯著改變主梁斷面的氣動(dòng)外形,導(dǎo)致橋梁在一定的風(fēng)速區(qū)間內(nèi)發(fā)生渦激振動(dòng)。在橋面兩側(cè)布置水馬的情況下,流體繞過(guò)水馬這一“鈍體”,在橋梁上部及背風(fēng)側(cè)出現(xiàn)明顯的漩渦脫落現(xiàn)象。

    2) 橋面兩側(cè)布置水馬的工況下, 7 m·s-1的風(fēng)速處于橋梁發(fā)生渦激振動(dòng)的風(fēng)速區(qū)間內(nèi),推測(cè)橋梁在該工況下有較大可能發(fā)生大幅渦激振動(dòng)。

    3) 相較于橋面兩側(cè)布置水馬,橋面中央布置水馬的升力和扭矩卓越頻率與橋梁一階對(duì)稱豎彎和一階正對(duì)稱扭轉(zhuǎn)模態(tài)頻率相差較大,且該工況下周圍流場(chǎng)隨著時(shí)間跨度的增大逐漸穩(wěn)定,推測(cè)發(fā)生橋梁渦激振動(dòng)的可能性較小。

    4) 在橋梁兩側(cè)布置水馬會(huì)提高橋梁在特定風(fēng)速下發(fā)生渦激振動(dòng)的可能性。在橋梁進(jìn)行市政施工時(shí),建議避免采用水馬沿橋梁兩側(cè)擺放的布置方式。流體經(jīng)過(guò)鏤空率較大的欄桿時(shí),并未在尾部產(chǎn)生漩渦脫落現(xiàn)象,建議使用開(kāi)孔率較大的水馬替代無(wú)孔水馬,以降低橋梁渦激振動(dòng)的風(fēng)險(xiǎn)。

    參考文獻(xiàn):

    [1] FAZL E L,SCANLAN R H.Vortex-induced vibrations of flexible bridges[J].Journal of Engineering Mechanics,1990,116(6):1392-1411.DOI:10.1061/(ASCE)0733-9399(1990)116:6(1392).

    [2] 張志田,陳政清.橋梁節(jié)段與實(shí)橋渦激共振幅值的換算關(guān)系[J].土木工程學(xué)報(bào),2011,44(7):77-82.DOI:10.15951/j.tmgcxb.2011.07.009.

    [3] 祝志文,石亞光,顏爽.帶防撞欄桿扁平箱梁高階模態(tài)渦激振動(dòng)的CFD研究[J].振動(dòng)與沖擊,2021,40(10):228-234,277.DOI:10.13465/j.cnki.jvs.2021.10.029.

    [4] 陳平,陳強(qiáng)草,華旭剛,等.結(jié)構(gòu)-尾流振子耦合模型參數(shù)識(shí)別及橋梁渦激振動(dòng)預(yù)測(cè)[J].世界橋梁,2021,49(2):84-89.DOI:10.3969/j.issn.1671-7767.2021.02.014.

    [5] 樊書(shū)文,李永樂(lè).風(fēng)嘴形式對(duì)橋梁斷面三分力系數(shù)影響的數(shù)值模擬[J].四川建筑,2018,38(3):163-165.

    [6] 葛耀君,趙林,許坤.大跨橋梁主梁渦激振動(dòng)研究進(jìn)展與思考[J].中國(guó)公路學(xué)報(bào),2019,32(10):1-18.DOI:10.19721/j.cnki.1001-7372.2019.10.001.

    [7] 華旭剛,陳智文,陳政清.大跨度懸索橋的多階模態(tài)豎向渦振與控制[J].中國(guó)公路學(xué)報(bào),2019,32(10):115-124.DOI:10.19721/j.cnki.1001-7372.2019.10.011.

    [8] 黃俐,周帥,梁鵬.橋梁斷面兩個(gè)渦振鎖定區(qū)間的數(shù)值模擬研究[J].振動(dòng)與沖擊,2016,35(11):47-53,73.DOI:10.13465/j.cnki.jvs.2016.11.008.

    [9] 王維民,陳楚龍,劉葉,等.豎彎渦振作用下橋上列車行車安全舒適性研究[J].橋梁建設(shè),2023,53(增刊1):38-45.DOI:10.20051/j.issn.1003-4722.2023.S1.006.

    [10] SARWAR M W,ISHIHARA T.Numerical study on suppression of vortex-induced vibrations of box girder bridge section by aerodynamic countermeasures[J].Journal of Wind Engineering and Industrial Aerodynamics,2010,98(12):701-711.DOI:10.1007/s12205-018-0548-y.

    [11] 鮮榮,廖海黎,李明水.大跨度橋梁主梁沿跨向渦激振動(dòng)響應(yīng)計(jì)算[J].西南交通大學(xué)學(xué)報(bào),2008,43(6):740-746.DOI:10.3969/j.issn.0258-2724.2008.06.011.

    [12] 陳精一,蔡國(guó)忠.電腦輔助工程分析ANSYS使用指南[M].北京:中國(guó)鐵道出版社,2001.

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

    [14] NAKAMURA Y,OHYA Y,TSURUTA H.Experiments on vortex shedding from flat plates with square leading and trailing edges[J].Journal of Fluid Mechanics,1991,222:437-447.DOI:10.1017/S0022112091001167.

    [15] WISSINK J G.DNS of separating low Reynolds number flow in a turbine cascade with incoming wakes[J].International Journal of Heat and Fluid Flow,2003,24(4):626-635.DOI:10.1016/S0142-727X(03)00056-0.

    (責(zé)任編輯: 錢(qián)筠 ""英文審校: 方德平)

    通信作者: 陳昌萍(1971-),男,教授,博士,博士生導(dǎo)師,主要從事橋梁風(fēng)工程的研究。E-mail:cpchen@hnu.edu.cn。

    基金項(xiàng)目: 國(guó)家自然科學(xué)基金資助項(xiàng)目(52178510); 福建省科技計(jì)劃項(xiàng)目(引導(dǎo)性項(xiàng)目)(2021Y0042)https://hdxb.hqu.edu.cn/

    猜你喜歡
    橋梁
    一種橋梁伸縮縫防滲水裝置
    手拉手 共搭愛(ài)的橋梁
    句子也需要橋梁
    橋梁的加固研究
    江西建材(2018年4期)2018-04-10 12:37:32
    加固技術(shù)創(chuàng)新,為橋梁健康保駕護(hù)航
    無(wú)人機(jī)在橋梁檢測(cè)中的應(yīng)用
    高性能砼在橋梁中的應(yīng)用
    現(xiàn)代鋼橋制造對(duì)橋梁鋼的更高要求
    焊接(2016年8期)2016-02-27 13:05:15
    城鄉(xiāng)建設(shè)一體化要注重橋梁的建筑設(shè)計(jì)
    南昌54座橋梁進(jìn)行兩個(gè)月的夏季體檢
    国产亚洲av嫩草精品影院| 女性被躁到高潮视频| 麻豆成人av在线观看| 国产亚洲精品第一综合不卡| 国产精品日韩av在线免费观看| 视频在线观看一区二区三区| 精品日产1卡2卡| 国产熟女xx| 欧美亚洲日本最大视频资源| 黄片小视频在线播放| av超薄肉色丝袜交足视频| 国产成人欧美在线观看| 老司机深夜福利视频在线观看| 韩国精品一区二区三区| 成人亚洲精品一区在线观看| 国产成人av教育| 怎么达到女性高潮| 欧美一区二区精品小视频在线| 久久精品国产亚洲av香蕉五月| 亚洲七黄色美女视频| 十分钟在线观看高清视频www| 久久 成人 亚洲| 麻豆久久精品国产亚洲av| 国产高清激情床上av| 国产成人精品无人区| 国产单亲对白刺激| 中文字幕精品免费在线观看视频| 男女午夜视频在线观看| 黄频高清免费视频| 久久久久国产一级毛片高清牌| 久久国产亚洲av麻豆专区| 超碰成人久久| 精品国产超薄肉色丝袜足j| 亚洲欧美日韩无卡精品| 欧美午夜高清在线| 亚洲人成77777在线视频| 精品久久久久久久末码| 99在线视频只有这里精品首页| 桃色一区二区三区在线观看| 人妻久久中文字幕网| 最新美女视频免费是黄的| 国产激情偷乱视频一区二区| 亚洲专区字幕在线| av免费在线观看网站| 亚洲精品在线美女| 亚洲精品在线美女| 午夜福利在线在线| 久久久久亚洲av毛片大全| 一a级毛片在线观看| 国产精品免费一区二区三区在线| 美女高潮到喷水免费观看| 狠狠狠狠99中文字幕| 给我免费播放毛片高清在线观看| www国产在线视频色| 淫妇啪啪啪对白视频| 久热爱精品视频在线9| 日韩欧美一区视频在线观看| 日韩高清综合在线| 一级毛片女人18水好多| 好男人电影高清在线观看| 亚洲中文av在线| 久99久视频精品免费| 美女扒开内裤让男人捅视频| 超碰成人久久| 在线永久观看黄色视频| 成年人黄色毛片网站| 国产成人av激情在线播放| 妹子高潮喷水视频| 亚洲av日韩精品久久久久久密| 国产国语露脸激情在线看| 99国产精品99久久久久| 亚洲精品在线观看二区| 日日干狠狠操夜夜爽| 91麻豆av在线| 日本成人三级电影网站| АⅤ资源中文在线天堂| 99久久精品国产亚洲精品| 国产亚洲精品av在线| 精品人妻1区二区| 国产精品自产拍在线观看55亚洲| 久久久久国产一级毛片高清牌| 少妇熟女aⅴ在线视频| 美女免费视频网站| 免费看十八禁软件| 人妻久久中文字幕网| 欧美日韩精品网址| 男女床上黄色一级片免费看| 精品熟女少妇八av免费久了| 欧美乱色亚洲激情| 午夜a级毛片| 日本 欧美在线| 欧美国产日韩亚洲一区| 亚洲精品中文字幕在线视频| 成在线人永久免费视频| 国产成人啪精品午夜网站| 变态另类成人亚洲欧美熟女| 亚洲精品美女久久av网站| 欧美日韩一级在线毛片| 国产精品亚洲美女久久久| 久久久久亚洲av毛片大全| 欧美性猛交╳xxx乱大交人| 久久久久久亚洲精品国产蜜桃av| 色婷婷久久久亚洲欧美| 亚洲成人国产一区在线观看| 午夜老司机福利片| 欧美激情久久久久久爽电影| 亚洲成av片中文字幕在线观看| 一级作爱视频免费观看| av电影中文网址| 国产精品久久久av美女十八| 亚洲精品久久国产高清桃花| 日日爽夜夜爽网站| 色在线成人网| 一边摸一边抽搐一进一小说| 可以在线观看毛片的网站| 成人国产一区最新在线观看| 国产成人精品久久二区二区91| 国产精品久久电影中文字幕| 不卡一级毛片| 免费看美女性在线毛片视频| 高清在线国产一区| 亚洲全国av大片| or卡值多少钱| 国产久久久一区二区三区| 黄网站色视频无遮挡免费观看| 亚洲成av片中文字幕在线观看| 人人妻人人澡人人看| 国产伦在线观看视频一区| 一级a爱片免费观看的视频| a级毛片在线看网站| 久久国产精品人妻蜜桃| 国产av在哪里看| 亚洲欧美一区二区三区黑人| 精品久久久久久久久久免费视频| 99国产极品粉嫩在线观看| 亚洲av第一区精品v没综合| 日本三级黄在线观看| 亚洲熟女毛片儿| 免费高清视频大片| www国产在线视频色| 亚洲专区国产一区二区| 91成年电影在线观看| 国产v大片淫在线免费观看| 国产av在哪里看| 日韩欧美一区二区三区在线观看| 在线观看免费午夜福利视频| 欧美激情久久久久久爽电影| 欧美黑人欧美精品刺激| 女人被狂操c到高潮| 精品福利观看| 琪琪午夜伦伦电影理论片6080| 女人被狂操c到高潮| 在线观看66精品国产| 欧美黄色片欧美黄色片| 亚洲一区高清亚洲精品| 精品福利观看| 精品久久久久久久毛片微露脸| 精品午夜福利视频在线观看一区| 一区二区三区高清视频在线| 亚洲一码二码三码区别大吗| 久久青草综合色| 村上凉子中文字幕在线| 国产精华一区二区三区| 免费在线观看视频国产中文字幕亚洲| 美女高潮喷水抽搐中文字幕| 欧美黄色淫秽网站| 国产一区二区三区视频了| 亚洲性夜色夜夜综合| svipshipincom国产片| 免费在线观看影片大全网站| 大型av网站在线播放| 亚洲在线自拍视频| 搡老岳熟女国产| 在线观看免费日韩欧美大片| 免费在线观看视频国产中文字幕亚洲| 国产欧美日韩精品亚洲av| 性欧美人与动物交配| 在线观看66精品国产| 国产一区二区在线av高清观看| 精品无人区乱码1区二区| 美女国产高潮福利片在线看| 国产免费男女视频| 少妇粗大呻吟视频| 亚洲欧美一区二区三区黑人| 免费一级毛片在线播放高清视频| 久久精品夜夜夜夜夜久久蜜豆 | 琪琪午夜伦伦电影理论片6080| 免费看日本二区| 中文亚洲av片在线观看爽| 成人免费观看视频高清| 日韩欧美国产在线观看| 国产成人系列免费观看| 后天国语完整版免费观看| 国产亚洲精品av在线| 少妇被粗大的猛进出69影院| 两人在一起打扑克的视频| 黑人巨大精品欧美一区二区mp4| 精品日产1卡2卡| 亚洲精品中文字幕在线视频| 欧美zozozo另类| 久久青草综合色| 99re在线观看精品视频| 1024手机看黄色片| 又黄又粗又硬又大视频| 在线免费观看的www视频| aaaaa片日本免费| 黄色成人免费大全| a级毛片a级免费在线| www.熟女人妻精品国产| 天天躁狠狠躁夜夜躁狠狠躁| 亚洲一卡2卡3卡4卡5卡精品中文| 亚洲国产精品合色在线| 亚洲精品粉嫩美女一区| 最近最新中文字幕大全电影3 | 亚洲成av人片免费观看| 亚洲色图av天堂| 亚洲中文av在线| 亚洲第一电影网av| 一级a爱片免费观看的视频| tocl精华| 亚洲国产精品合色在线| 欧美性猛交╳xxx乱大交人| 日韩欧美国产在线观看| 嫩草影视91久久| 欧美激情高清一区二区三区| 中文字幕久久专区| 亚洲激情在线av| 亚洲免费av在线视频| 国产视频一区二区在线看| 亚洲性夜色夜夜综合| 日本一本二区三区精品| av欧美777| 日本一区二区免费在线视频| 999久久久精品免费观看国产| 日本免费a在线| 色综合欧美亚洲国产小说| 国产精品免费视频内射| 19禁男女啪啪无遮挡网站| av福利片在线| 久久精品国产综合久久久| 中文字幕人成人乱码亚洲影| 午夜精品在线福利| 白带黄色成豆腐渣| 黄频高清免费视频| 黄色毛片三级朝国网站| 日本精品一区二区三区蜜桃| 校园春色视频在线观看| 最好的美女福利视频网| 18禁黄网站禁片午夜丰满| 一级黄色大片毛片| а√天堂www在线а√下载| 欧美最黄视频在线播放免费| 国产精品香港三级国产av潘金莲| 桃色一区二区三区在线观看| 色哟哟哟哟哟哟| 亚洲aⅴ乱码一区二区在线播放 | 亚洲国产毛片av蜜桃av| 美女 人体艺术 gogo| 黄网站色视频无遮挡免费观看| 国产伦人伦偷精品视频| 国产亚洲精品综合一区在线观看 | 欧美乱妇无乱码| 日韩精品中文字幕看吧| 国产伦人伦偷精品视频| 亚洲一区二区三区不卡视频| 精品久久久久久久毛片微露脸| 日本一本二区三区精品| 国产一卡二卡三卡精品| 欧美国产日韩亚洲一区| 日韩精品免费视频一区二区三区| 久久草成人影院| 日韩免费av在线播放| 久久精品91无色码中文字幕| 激情在线观看视频在线高清| 免费女性裸体啪啪无遮挡网站| 免费观看精品视频网站| 亚洲男人天堂网一区| 精品久久久久久成人av| 99re在线观看精品视频| 老司机午夜福利在线观看视频| 好男人电影高清在线观看| 禁无遮挡网站| 日本黄色视频三级网站网址| www.www免费av| www日本在线高清视频| 琪琪午夜伦伦电影理论片6080| 久久热在线av| 国产一卡二卡三卡精品| 高清在线国产一区| 日韩欧美一区视频在线观看| 一卡2卡三卡四卡精品乱码亚洲| 亚洲av成人不卡在线观看播放网| 一级作爱视频免费观看| 成人特级黄色片久久久久久久| 大型av网站在线播放| 好看av亚洲va欧美ⅴa在| 19禁男女啪啪无遮挡网站| 最新美女视频免费是黄的| 国产精品香港三级国产av潘金莲| 成熟少妇高潮喷水视频| 搡老岳熟女国产| 国产日本99.免费观看| 久久久久久人人人人人| 国产欧美日韩一区二区精品| 999精品在线视频| 1024视频免费在线观看| 又黄又爽又免费观看的视频| 狠狠狠狠99中文字幕| 91九色精品人成在线观看| 啦啦啦 在线观看视频| 两个人看的免费小视频| 午夜免费观看网址| 天堂影院成人在线观看| 久久亚洲精品不卡| 国产精品 欧美亚洲| 美女午夜性视频免费| 欧美黄色淫秽网站| 天天躁狠狠躁夜夜躁狠狠躁| www国产在线视频色| 国产av一区二区精品久久| 日韩有码中文字幕| 亚洲全国av大片| 国产欧美日韩一区二区精品| 悠悠久久av| 高清毛片免费观看视频网站| 成人特级黄色片久久久久久久| 久久天躁狠狠躁夜夜2o2o| 99国产精品一区二区三区| 免费看美女性在线毛片视频| 一本大道久久a久久精品| or卡值多少钱| 好看av亚洲va欧美ⅴa在| 久久久久久久午夜电影| 中文字幕高清在线视频| 黄色视频不卡| 可以在线观看毛片的网站| 十八禁网站免费在线| 精品第一国产精品| 国产成年人精品一区二区| 亚洲天堂国产精品一区在线| 久久中文看片网| 天堂影院成人在线观看| 午夜免费观看网址| 国产精品一区二区免费欧美| 免费在线观看完整版高清| 成人18禁在线播放| 91麻豆精品激情在线观看国产| √禁漫天堂资源中文www| 国产不卡一卡二| 精品午夜福利视频在线观看一区| 久久狼人影院| www.熟女人妻精品国产| 亚洲成人国产一区在线观看| 精品一区二区三区av网在线观看| 亚洲精品美女久久av网站| 精品国产乱子伦一区二区三区| 宅男免费午夜| 国产在线精品亚洲第一网站| 国产av不卡久久| 国内精品久久久久精免费| 黄片小视频在线播放| 制服人妻中文乱码| 天堂动漫精品| 日日爽夜夜爽网站| 免费高清在线观看日韩| 一本大道久久a久久精品| 日韩欧美免费精品| 日韩精品中文字幕看吧| a级毛片在线看网站| 婷婷精品国产亚洲av| 琪琪午夜伦伦电影理论片6080| 精品国内亚洲2022精品成人| 亚洲av日韩精品久久久久久密| 看片在线看免费视频| 琪琪午夜伦伦电影理论片6080| 日日摸夜夜添夜夜添小说| 国产三级黄色录像| 亚洲av片天天在线观看| 久久中文字幕人妻熟女| 天天一区二区日本电影三级| 亚洲黑人精品在线| 久热爱精品视频在线9| 久久久久国内视频| 妹子高潮喷水视频| 18禁黄网站禁片午夜丰满| 97碰自拍视频| 亚洲成国产人片在线观看| 国产精品久久久av美女十八| 欧美中文综合在线视频| 18禁国产床啪视频网站| 免费在线观看日本一区| 亚洲av电影在线进入| 久久久久久久久中文| 淫秽高清视频在线观看| 成人永久免费在线观看视频| 两个人视频免费观看高清| 手机成人av网站| 99久久无色码亚洲精品果冻| 亚洲熟女毛片儿| tocl精华| 欧美又色又爽又黄视频| www国产在线视频色| 中文字幕人成人乱码亚洲影| 国产高清videossex| 免费电影在线观看免费观看| 欧美最黄视频在线播放免费| 成在线人永久免费视频| 女性被躁到高潮视频| 亚洲国产欧美一区二区综合| 亚洲三区欧美一区| 久9热在线精品视频| 亚洲专区中文字幕在线| 夜夜夜夜夜久久久久| 97人妻精品一区二区三区麻豆 | 午夜免费成人在线视频| 免费看日本二区| 日本 av在线| 久久久久精品国产欧美久久久| 亚洲激情在线av| 我的亚洲天堂| 国产精品爽爽va在线观看网站 | 国产精品九九99| 色综合婷婷激情| 免费女性裸体啪啪无遮挡网站| 成人国产综合亚洲| 禁无遮挡网站| 日韩欧美国产一区二区入口| 特大巨黑吊av在线直播 | 精品国内亚洲2022精品成人| 久久久精品欧美日韩精品| 亚洲自偷自拍图片 自拍| cao死你这个sao货| 欧美三级亚洲精品| 99热6这里只有精品| 亚洲午夜理论影院| 亚洲国产日韩欧美精品在线观看 | 国产99久久九九免费精品| 久久久久久久午夜电影| 国产精品久久久av美女十八| 成人国语在线视频| 黄色女人牲交| 久久青草综合色| 听说在线观看完整版免费高清| 中文字幕精品免费在线观看视频| 1024视频免费在线观看| 国产精品综合久久久久久久免费| 久久久久久亚洲精品国产蜜桃av| 久热这里只有精品99| 亚洲在线自拍视频| 欧美性猛交╳xxx乱大交人| 色老头精品视频在线观看| 18禁裸乳无遮挡免费网站照片 | a在线观看视频网站| 欧美日韩黄片免| x7x7x7水蜜桃| 日韩精品免费视频一区二区三区| 黄色a级毛片大全视频| 亚洲无线在线观看| 女性被躁到高潮视频| 俺也久久电影网| 亚洲专区中文字幕在线| 人妻久久中文字幕网| 日韩欧美三级三区| 91在线观看av| 国产精品久久视频播放| 亚洲专区中文字幕在线| 午夜福利成人在线免费观看| 亚洲av美国av| 好看av亚洲va欧美ⅴa在| av超薄肉色丝袜交足视频| 麻豆久久精品国产亚洲av| 97超级碰碰碰精品色视频在线观看| 一a级毛片在线观看| 男女那种视频在线观看| 露出奶头的视频| 欧美激情久久久久久爽电影| 一本一本综合久久| 久久国产乱子伦精品免费另类| 国产久久久一区二区三区| 国产av不卡久久| 无限看片的www在线观看| 中亚洲国语对白在线视频| 国产高清视频在线播放一区| 国产v大片淫在线免费观看| 黄片小视频在线播放| 国产成人精品久久二区二区免费| 夜夜爽天天搞| 久久久久久人人人人人| 精品日产1卡2卡| 日本精品一区二区三区蜜桃| 国产伦在线观看视频一区| 一边摸一边抽搐一进一小说| 亚洲av成人不卡在线观看播放网| 亚洲国产精品久久男人天堂| 哪里可以看免费的av片| 熟妇人妻久久中文字幕3abv| 在线观看免费日韩欧美大片| 欧美日本亚洲视频在线播放| 国产aⅴ精品一区二区三区波| 久久久久久大精品| 成年版毛片免费区| 国产精品美女特级片免费视频播放器 | 午夜激情福利司机影院| 久久狼人影院| 一级毛片精品| 99久久无色码亚洲精品果冻| 999久久久国产精品视频| 成熟少妇高潮喷水视频| 亚洲精品美女久久av网站| 国产免费男女视频| 一本一本综合久久| 18禁黄网站禁片免费观看直播| xxx96com| 国产伦人伦偷精品视频| 国产99白浆流出| 国产成人欧美在线观看| 亚洲成av片中文字幕在线观看| 久久精品亚洲精品国产色婷小说| 午夜免费激情av| 丰满人妻熟妇乱又伦精品不卡| 久久久水蜜桃国产精品网| 欧美绝顶高潮抽搐喷水| 亚洲色图 男人天堂 中文字幕| 日韩成人在线观看一区二区三区| 观看免费一级毛片| 中文字幕另类日韩欧美亚洲嫩草| 欧美三级亚洲精品| 久久草成人影院| 国产精品一区二区三区四区久久 | 淫秽高清视频在线观看| 一级a爱视频在线免费观看| 无人区码免费观看不卡| 国产av一区二区精品久久| 热re99久久国产66热| 性色av乱码一区二区三区2| 午夜两性在线视频| 性欧美人与动物交配| 岛国视频午夜一区免费看| 熟妇人妻久久中文字幕3abv| 欧美成狂野欧美在线观看| 欧美av亚洲av综合av国产av| 国产在线精品亚洲第一网站| 51午夜福利影视在线观看| 黄色 视频免费看| 欧美成人免费av一区二区三区| 天堂√8在线中文| 欧美激情久久久久久爽电影| 女警被强在线播放| 婷婷六月久久综合丁香| 亚洲中文日韩欧美视频| 久久精品91无色码中文字幕| 国产男靠女视频免费网站| 亚洲天堂国产精品一区在线| 中文字幕另类日韩欧美亚洲嫩草| 亚洲一区二区三区色噜噜| 精品日产1卡2卡| 国产一区二区激情短视频| 日韩精品中文字幕看吧| 人人妻人人澡人人看| 色综合亚洲欧美另类图片| 91九色精品人成在线观看| 午夜福利视频1000在线观看| 国产精品久久久av美女十八| 窝窝影院91人妻| 亚洲国产欧美一区二区综合| av中文乱码字幕在线| 成人18禁高潮啪啪吃奶动态图| 欧美日韩乱码在线| 日韩欧美 国产精品| 亚洲精品国产一区二区精华液| 欧美性猛交黑人性爽| 熟女电影av网| 首页视频小说图片口味搜索| 亚洲精品美女久久av网站| 白带黄色成豆腐渣| 在线看三级毛片| 男女床上黄色一级片免费看| 亚洲一区中文字幕在线| 非洲黑人性xxxx精品又粗又长| 国产精品,欧美在线| 国产成+人综合+亚洲专区| 欧美人与性动交α欧美精品济南到| 亚洲精品久久成人aⅴ小说| 91在线观看av| tocl精华| 两个人免费观看高清视频| 亚洲午夜理论影院| tocl精华| 两个人免费观看高清视频| 香蕉久久夜色| 中文字幕人妻熟女乱码| 一进一出好大好爽视频| 国产精品亚洲美女久久久| 亚洲一卡2卡3卡4卡5卡精品中文| 午夜福利一区二区在线看| 国产伦在线观看视频一区| 别揉我奶头~嗯~啊~动态视频| 亚洲片人在线观看| 搡老熟女国产l中国老女人| 最近在线观看免费完整版| 日韩欧美一区二区三区在线观看| 无限看片的www在线观看| 亚洲人成电影免费在线| 他把我摸到了高潮在线观看| 日韩国内少妇激情av| 99精品久久久久人妻精品| 国产精品久久久人人做人人爽| 亚洲自偷自拍图片 自拍| 国产麻豆成人av免费视频| 久久久久久久久久黄片| 久久久久国产精品人妻aⅴ院| 免费在线观看成人毛片| 十八禁网站免费在线| 亚洲专区国产一区二区| 可以免费在线观看a视频的电影网站|