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

    基于精確阻尼調(diào)控的橋梁豎彎渦振Sc數(shù)影響

    2022-05-30 10:48:04華旭剛馬偉猛黃智文陳政清萬田保孫英杰劉曙光
    關(guān)鍵詞:渦振阻尼比橋梁

    華旭剛 馬偉猛 黃智文 陳政清 萬田保 孫英杰 劉曙光

    摘 要:在橋梁氣動外形和來流風(fēng)場確定的條件下,Sc數(shù)(Scruton number)是影響大跨度橋梁渦振響應(yīng)大小的關(guān)鍵因素.為了開展Sc數(shù)對大跨度橋梁豎向渦振影響的精細化研究,首 先研制了適用于橋梁節(jié)段模型風(fēng)洞試驗的永磁式板式電渦流阻尼器,可為彈性懸掛節(jié)段模型系統(tǒng)提供可連續(xù)調(diào)節(jié)的、理想線性黏滯阻尼.然后以帶風(fēng)嘴的開口斷面鋼混組合梁橋為研究對象,針對+3°、0°和-3°三個風(fēng)攻角和不同的Sc數(shù)開展了節(jié)段模型豎向渦振試驗.根據(jù)試驗結(jié)果,總結(jié)了節(jié)段模型渦振振幅的風(fēng)速變化曲線隨名義阻尼比的變化規(guī)律,分析了最大渦振振 幅隨Sc數(shù)變化規(guī)律的函數(shù)擬合,并對橋梁高階模態(tài)渦振響應(yīng)進行了預(yù)測.研究發(fā)現(xiàn)不同風(fēng)攻角下主梁豎向渦振的鎖定風(fēng)速區(qū)間都隨Sc數(shù)的增大而縮小,最大振幅都隨Sc數(shù)的增加而降 低,但變化曲線有顯著差別;通過對三組及以上不同Sc數(shù)的節(jié)段模型風(fēng)洞試驗結(jié)果進行函數(shù)擬合,可以較好地預(yù)測節(jié)段模型最大渦振振幅隨Sc數(shù)的變化規(guī)律;在模態(tài)阻尼比相同的條件下,大跨度懸索橋各階豎彎模態(tài)的渦振振幅基本相同,高階模態(tài)渦振更容易出現(xiàn)不滿足規(guī)范限值的情況.研究成果為大跨度橋梁豎彎渦振響應(yīng)的精細化預(yù)測提供了重要方法.

    關(guān)鍵詞:渦振;橋梁;節(jié)段模型;阻尼比;電渦流阻尼

    中圖分類號:U441.3文獻標志碼:A

    Effect of Sc Number on Vertical Vortex-induced Vibration of Bridges Based on Precise Damping Adjustment

    HUA Xugang1,2,MA Weimeng1,2,HUANG Zhiwen1,2?,CHEN Zhengqing1,2,WAN Tianbao3,SUN Yingjie4,LIU Shuguang4

    (1.Key Laboratory for Wind and Bridge Engineering of Hunan Province(Hunan University),Changsha 410082,China;

    2.College of Civil Engineering,Hunan University,Changsha 410082,China;

    3.China Railway Major Bridge Reconnaissance & Design Institute Co,Ltd,Wuhan 430056,China;

    4.China Railway Bridge & Tunnel Technologies Co,Ltd.,Nanjing 210061,China)

    Abstract:Under the condition that the aerodynamic shape of the bridge and the incoming wind field are deter-mined,the Sc number is one of the key factors affecting the vortex-induced response of long-span bridges.In orderto conduct a refinement study on the effect of Sc number on vertical vortex-induced vibration(VVIV)of long-span bridges,firstly a permanent magnet plate eddy current damper was developed,and it was utilized to provide linear and adjustable viscous damping for the spring-suspended sectional model system.Then wind tunnel test for a section model of a steel-concrete composite girder bridge with wind fairings was conducted under three different wind attack angles of +3°,0°and-3°,respectively,and several different Sc number.Based on the experiments,the variation of VVIV displacement with the wind velocity for different nominal damping ratio was obtained,the maximum displace-ment of the VVIV versus Sc number relationship was analyzed by curve fitting,and the higher-order modal vortex-induced vibration response of the bridge was predicted.It is found that the lock-in wind speed region as well as the maximum amplitude decrease with the increase of Sc number under all the three wind attack angles,though their varying tendencies are significantly different.The maximum amplitudes for different Sc numbers can be well pre-dicted by curve fitting of the results of wind tunnel tests for at least three different Sc numbers.On the assumption of the same modal damping ratio,the vortex-induced vibration amplitudes of all vertical bending modes of long-span suspension bridges are almost the same.As a result,the vortex-induced vibration amplitudes of higher-order modes are more difficult to meet the specification limits.The research results provide a novel method for refinement predic-tion of VVIV response of long-span bridges.

    Key words:vortex-induced vibration;bridge;sectional model;damping ratio;eddy current damping

    主梁豎向渦振是大跨度橋梁在常遇風(fēng)速下容易發(fā)生的一種風(fēng)致振動現(xiàn)象.主梁渦振響應(yīng)的鎖定風(fēng)速區(qū)間和幅值大小主要受主梁斷面氣動外形、來流 特性和結(jié)構(gòu)動力特性等因素的影響[1-3].在結(jié)構(gòu)動力特性參數(shù)中,模態(tài)阻尼比和質(zhì)量是兩個重要參數(shù),它 們對主梁渦振響應(yīng)的影響通常可歸納為無量綱參數(shù)Scruton數(shù)的影響[4].Scruton數(shù)簡稱Sc數(shù),它是結(jié)構(gòu)無量綱質(zhì)量和阻尼比的乘積,最早由 Scruton在對煙囪等圓柱體結(jié)構(gòu)渦振的研究中提出[5].對于主梁豎向渦振,文獻[6]對Sc數(shù)的原始定義進行了修正,其表達式為Sc=4πmξ/(ρBD),其中m表示主梁單位長度的質(zhì)量;ξ表示結(jié)構(gòu)模態(tài)阻尼比;ρ表示空氣密度;B和D分別表示主梁寬度和高度.本文將采用Sc數(shù)的上述定義形式.

    研究 Sc數(shù)對主梁渦振的影響是實橋渦振響應(yīng)預(yù) 測及控制的重要環(huán)節(jié).一方面,大跨度橋梁各階豎彎模態(tài)的等效質(zhì)量和阻尼比都可能存在差異,因此要 預(yù)測不同模態(tài)的渦振響應(yīng),就需要分析 Sc數(shù)的影響 規(guī)律.郭震山等[7]和王守強等[8]分別以虎門大橋和西 堠門大橋為背景對大跨度懸索橋的模態(tài)阻尼比進行了統(tǒng)計分析,結(jié)果發(fā)現(xiàn)大跨度懸索橋的豎彎模態(tài)阻尼比一般隨模態(tài)階次的增大而減小.實際上,由于大跨度橋梁阻尼比實測統(tǒng)計資料較少,目前國內(nèi)外不 同橋梁規(guī)范對阻尼比的建議值也存在較大差異.例如,對于鋼混組合纜索承重橋梁,我國抗風(fēng)設(shè)計規(guī)范建議的模態(tài)阻尼比為1%[9],而英國橋梁設(shè)計規(guī)范BD 49/01的建議值則為0.48%[10].另一方面,當采用調(diào)諧質(zhì)量減振器等機械阻尼措施進行渦振控制時,可首先確定主梁渦振響應(yīng)隨模態(tài)阻尼比或Sc數(shù)的變化規(guī)律,然后再根據(jù)所需的附加模態(tài)阻尼比或Sc數(shù)開展阻尼器參數(shù)設(shè)計.

    目前,有大量文獻報道了氣動外形和來流特性對橋梁渦振的影響規(guī)律[11-12],但有關(guān)Sc數(shù)對橋梁渦振的影響仍缺乏系統(tǒng)的研究.在試 驗研究方面,Klamo等[13]研制了電磁電渦流阻尼裝置,并依托該裝置分析了阻尼對水中圓柱體結(jié)構(gòu)橫流向渦振振幅和頻率的影響規(guī)律.Larsen[14]基于風(fēng)洞試驗研究了一個鈍體箱形斷面分別在均勻流場和紊流下的豎向渦振振幅隨Sc數(shù)的變化規(guī)律,并提出了一種渦激力的廣義范德玻振子經(jīng)驗?zāi)P?李永君等[15]以東海大橋顆珠山橋為背景研究了質(zhì)量和阻尼比對橋梁扭轉(zhuǎn) 渦振的影響,發(fā)現(xiàn)增大質(zhì)量參數(shù)或阻尼參數(shù)都可以減小扭轉(zhuǎn)渦振幅值.Marra等[16]通過節(jié)段模型風(fēng)洞試 驗研究了Sc數(shù)對寬高比為4∶1的矩形斷面豎向渦振的影響規(guī)律,并由此提出了改進的Scanlan非線性經(jīng)驗渦激力模型,但只考慮了0°風(fēng)攻角的試驗工況.在理論研究方面,為了提高橋梁渦振響應(yīng)的預(yù)測精度,近年來許多學(xué)者致力于拓展半經(jīng)驗渦激力數(shù)學(xué)模型的適用性,使之能夠在不同的Sc數(shù)下成立[17-19].在針對海洋立管等圓柱體結(jié)構(gòu)的流致渦振研究中,Van-diver[20]提出了一種新的無量綱阻尼參數(shù),并分析了它在三維渦振幅值響應(yīng)預(yù)測中的優(yōu)勢.Garcia等[21]分析了圓柱體結(jié)構(gòu)渦振中固有阻尼比對系統(tǒng)附加質(zhì)量和升力的影響.

    采用節(jié)段模型風(fēng)洞試驗和理論分析研究結(jié)構(gòu)阻尼比或Sc數(shù)對橋梁渦振響應(yīng)的影響,其關(guān)鍵在于為節(jié)段模型系統(tǒng)提供穩(wěn)定、可靠和精確可調(diào)的線性黏滯阻尼.目前橋梁節(jié)段模型風(fēng)洞試驗中主要采用油阻尼桶、纏繞鋼絲繩或電工膠帶等方式來調(diào)節(jié)節(jié)段模型懸掛系統(tǒng)的阻尼比,但上述方法都有很明顯的局限性.例如,油阻尼桶提供的剪切黏滯阻尼一般具有較好的線性特性,但其阻尼系數(shù)與剪切面的形狀和面積相關(guān),難以實現(xiàn)連續(xù)調(diào)節(jié);纏繞鋼絲繩或電工 膠帶產(chǎn)生的阻尼都在一定程度上依賴于材料之間的摩擦,因此都具有一定的非線性,也難以實現(xiàn)精確、連續(xù)調(diào)節(jié)[22].此外,若采用電磁電渦流阻尼,則可能存在電磁鐵發(fā)熱的問題.

    為了開展Sc數(shù)對大跨度橋梁豎向渦振影響的精細化研究,本文首先研制了適用于橋梁節(jié)段模型風(fēng) 洞試驗的、阻尼系數(shù)精確可調(diào)的永磁式板式電渦流阻尼器,結(jié)合電磁有限元分析和自由衰減振動試驗 驗證了阻尼器的阻尼特性.然后通過改變節(jié)段模型系統(tǒng)的附加阻尼比來改變系統(tǒng)的Sc數(shù),以帶風(fēng)嘴的開口斷面鋼混組合梁為研究對象,通過節(jié)段模型風(fēng) 洞試驗研究了實橋在+3°、0°和-3°三個風(fēng)攻角下的豎向渦振響應(yīng)隨Sc數(shù)的變化規(guī)律,并據(jù)此對橋梁高階模態(tài)渦振響應(yīng)進行了預(yù)測.

    1板式電渦流阻尼器研制及性能測試

    1.1永磁式板式電渦流阻尼器的基本構(gòu)造

    圖1給出了永磁式板式電渦流阻尼器的基本構(gòu)造.可以看到,板式電渦流阻尼器由永磁體、永磁體 背鐵、導(dǎo)體板、底座和間距調(diào)節(jié)裝置構(gòu)成.其中,導(dǎo)體 板固定在節(jié)段模型的端部,與節(jié)段模型一起沿豎向運動;底座固定在風(fēng)洞地面,其距離地面的高度可以調(diào)節(jié);永磁體分別布置在導(dǎo)體板的兩側(cè),其中沿導(dǎo)體 板運動方向的磁極交替布置,沿垂直于導(dǎo)體板運動方向的磁極相同布置.導(dǎo)體板兩側(cè)永磁體之間的空氣間隙可連續(xù)調(diào)節(jié),從而實現(xiàn)阻尼器的阻尼系數(shù)以及節(jié)段模型系統(tǒng)阻尼比的連續(xù)調(diào)節(jié).永磁體背鐵由高磁導(dǎo)率材料制作,這樣使阻尼器自身形成閉合磁路,不會對試驗儀器產(chǎn)生電磁干擾,同時還提高了電渦流阻尼器的耗能效率.永磁體采用釹鐵硼稀土永 磁體,厚度20mm;導(dǎo)體板采用鋁板,厚度6mm.

    本文只研究節(jié)段模型的豎向渦振,因此電渦流阻尼器安裝在節(jié)段模型左右兩端的中心位置,如圖2所示.當節(jié)段模型振動時,模型帶動導(dǎo)體板也做上下振動,從而使永磁體與導(dǎo)體板之間發(fā)生相對運動.于是導(dǎo)體板切割磁力線,并在其上產(chǎn)生電渦流,電渦流 產(chǎn)生的磁場與永磁體磁場相互作用,產(chǎn)生一個阻礙 導(dǎo)體板和永磁體相對運動的阻尼力,即電渦流阻尼力.從能量轉(zhuǎn)化的角度看,電渦流阻尼器把節(jié)段模型的振動能量首先轉(zhuǎn)化為導(dǎo)體板內(nèi)的電能,并最終轉(zhuǎn)化為熱能耗散掉,由此產(chǎn)生了阻尼效應(yīng).

    1.2 板式電渦流阻尼器性能

    本文以鋼混組合梁為研究對象,圖3給出了節(jié) 段模型橫斷面設(shè)計尺寸,其中模型寬度B=0.633m,高度D=0.046m,長度L=1.740m,每 延 米質(zhì)量m=11.273kg,豎向振動的固有頻率fn=2.87Hz,幾何縮尺比λL=60,試驗風(fēng)速比λU=2.1.

    在開展節(jié)段模型風(fēng)洞試驗前,通過電磁有限元分析和節(jié)段模型的自由衰減振動試驗對永磁式板式電渦流阻尼器的理想線性阻尼特性進行了驗證.在電磁有限元分析軟件 COMSOL中建立了板式電渦流阻尼器的有限元模型,分析了4種不同氣隙條件下板式電渦流阻尼器的阻尼系數(shù)隨模型振動速度的變化規(guī)律,并根據(jù)節(jié)段模型系統(tǒng)的質(zhì)量和自振頻率把 上述阻尼系數(shù)轉(zhuǎn)換為節(jié)段模型系統(tǒng)的附加結(jié)構(gòu)阻尼比.為了便于對比,不同振動速度下的附加結(jié)構(gòu)阻尼比均除以速度為0.01m/s時的附加阻尼比,以無量綱化.分析結(jié)果如圖4所示,可以看到,不管氣隙如何變化,節(jié)段模型系統(tǒng)的附加結(jié)構(gòu)阻尼比都幾乎不隨其振動速度(即導(dǎo)體鋁板的振動速度)變化,說明文中板式電渦流阻尼器的阻尼性能與模型振動狀態(tài)無關(guān).需要說明的是,根據(jù)模型固有頻率和以往試驗中鋼混組合梁斷面的渦振響應(yīng)大小,分析時保守估計了模型的最大振動速度為0.2m/s,對應(yīng)最大渦振振 幅為0.011m,約為0.24倍模型高度.

    板式電渦流阻尼器的線性阻尼特性還可以從安裝阻尼器前后節(jié)段模型系統(tǒng)的附加阻尼比隨振幅的變化上得到驗證.圖5給出了無阻尼器時由自由衰 減振動試驗得到的節(jié)段模型阻尼比和頻率隨振幅變化曲線,可以看到,節(jié)段模型系統(tǒng)的阻尼比隨節(jié)段模型振幅的增大而近似呈線性增大,但節(jié)段模型系統(tǒng)的固有頻率隨振幅變化非常小,可認為是常數(shù).值得注意的是,無阻尼器時由自由衰減振動試驗得到的阻尼包括因模型變形和構(gòu)件摩擦等因素產(chǎn)生的機械阻尼,以及無風(fēng)環(huán)境下的氣動阻尼[23].

    安裝板式電渦流阻尼器后,首先通過自由衰減 振動得到節(jié)段模型系統(tǒng)的總阻尼比隨振幅的變化曲線,再減去無阻尼器時系統(tǒng)的阻尼比(見圖5),就得到板式電渦流阻尼器提供的附加阻尼比.圖6根據(jù) 附加阻尼比隨振幅的變化結(jié)果進一步統(tǒng)計了附加阻尼比的平均值ξ0及其變動率δ隨氣隙大小的變化規(guī) 律.其中ξ0為1~10mm的振幅區(qū)間內(nèi),附加阻尼比的平均值;變動率δ=(ξmax-ξ min)/ξ0,其中ξmax和ξmin分別為1~10mm 振幅區(qū)間內(nèi),附加阻尼比的最大值和最小值.從圖6可知,不同氣隙下附加阻尼比的變動率δ都很小,大部分工況都在4%以內(nèi).因此,不同氣隙下附加阻尼比可認為是一常量,并可用平均值ξ0表示.

    需要說明的是,本文在分析阻尼或Sc數(shù)對節(jié)段模型渦振響應(yīng)的影響時,將計入阻尼器的附加阻尼、節(jié)段模型系統(tǒng)的機械阻尼和無風(fēng)環(huán)境下的氣動阻尼(即圖5中識別的系統(tǒng)阻尼),以下稱為總阻尼.實際 上,節(jié)段模型發(fā)生渦振時,應(yīng)將無風(fēng)環(huán)境下的氣動阻尼排除在影響因素之外,但要準確分離機械阻尼和無風(fēng)環(huán)境下的氣動阻尼難度較大[24-25].考慮到無風(fēng) 環(huán)境下的氣動阻尼相對附加阻尼較小,本文為了簡化分析,把無風(fēng)環(huán)境下的氣動阻尼合并到機械阻尼中,得到的渦振預(yù)測結(jié)果將偏于保守.

    為了展現(xiàn)板式電渦流阻尼器阻尼系數(shù)的連續(xù)可調(diào)功能,采用電磁有限元分析和節(jié)段模型系統(tǒng)的自由衰減振動試驗,對比分析了節(jié)段模型系統(tǒng)的附加阻尼比隨永磁體與導(dǎo)體板之間氣隙大小的變化規(guī) 律,結(jié)果如圖7所示.可以看到,隨著氣隙的減小,節(jié) 段模型系統(tǒng)的附加阻尼比逐漸增大.當氣隙大于4 cm時,由阻尼器提供的附加阻尼比幾乎為零;當氣隙減小到2 cm時,附加阻尼比可增大到2%,可完全 抑制節(jié)段模型的豎向渦振.此外可以發(fā)現(xiàn),由有限元分析預(yù)測的附加阻尼比與實測結(jié)果非常接近,證明了板式電渦流阻尼器設(shè)計方案的可行性.

    2? 阻尼及Sc數(shù)對節(jié)段模型渦振的影響

    2.1? 豎向渦振位移響應(yīng)風(fēng)速區(qū)間隨阻尼的變化

    因為機械阻尼和無風(fēng)環(huán)境下的氣動阻尼具有振 幅依賴性,所以當附加阻尼較小時,系統(tǒng)的總阻尼及對應(yīng)的Sc數(shù)也具有較明顯的振幅依賴性.橋梁渦振研究的主要目的之一是預(yù)測實橋最大渦振振幅,因 此可取各工況下節(jié)段模型最大渦振振幅對應(yīng)的總阻尼比作為該工況下的名義阻尼比,再由此計算Sc數(shù),即Sc= 4πmξ n /(ρBD),其中ξn表示節(jié)段模型名義阻尼比.

    圖8為沒有附加阻尼時,不同風(fēng)攻角下節(jié)段模型的無量綱豎向渦振振幅隨折減風(fēng)速的變化曲線,其中無量綱振幅定義為η=y/D,無量綱折減風(fēng)速定 義為Vr=U/fn B,其中U為來流風(fēng)速.可以看到,在+3?、0?和-3°來流風(fēng)攻角下節(jié)段模型系統(tǒng)的名義阻尼比均為0.1%,主梁均只有一個豎向渦振鎖定風(fēng)速區(qū)間,且不同風(fēng)攻角下的渦振鎖定風(fēng)速區(qū)間基本一致,無量綱最大渦振振幅約為0.12.

    通過板式電渦流阻尼器不斷增大節(jié)段模型系統(tǒng)的附加阻尼比,得到不同風(fēng)攻角下名義阻尼比不斷增大時,節(jié)段模型的無量綱豎向渦振振幅隨折減風(fēng)速的變化情況,如圖9所示.可以看到,不同風(fēng)攻角下,節(jié)段模型的豎向渦振最大振幅和渦振鎖定折減 風(fēng)速區(qū)間均隨著名義阻尼比的增大而逐漸減小,但變化趨勢不同.例如,對于+3? 風(fēng)攻角,即使名義阻尼比增大到1%,主梁的豎向渦振也沒有完全消失.而對于-3? 風(fēng)攻角,當系統(tǒng)的名義阻尼比略大于0.37%時,豎向渦振就基本消失.

    2.2 最大豎向渦振振幅隨Sc數(shù)的變化曲線

    根據(jù)上一節(jié)試驗結(jié)果進一步統(tǒng)計得到節(jié)段模型 豎向渦振的最大渦振振幅隨名義阻尼比及Sc數(shù)的變化曲線,分別如圖10和圖11所示.可以看到,隨著名 義阻尼比或Sc數(shù)的不斷增加,不同風(fēng)攻角下的最大渦振振幅不斷降低,但不同風(fēng)攻角下的變化規(guī)律不相同.+3? 風(fēng)攻角下的節(jié)段模型豎向渦振峰值響應(yīng)下降最慢,-3? 風(fēng)攻角下的節(jié)段模型豎向渦振峰值響應(yīng)下降最快,0? 風(fēng)攻角下的節(jié)段模型豎向渦振峰值響 應(yīng)下降速度介于前兩者之間.此外,隨著阻尼比或Sc數(shù)的增大,最大渦振振幅對阻尼比或Sc數(shù)的敏感程度逐漸降低.

    2.3最大豎向渦振振幅隨Sc數(shù)的變化曲線擬合

    利用板式電渦流阻尼器可實現(xiàn)節(jié)段模型系統(tǒng)阻尼比的連續(xù)調(diào)節(jié),得到最大渦振振幅隨阻尼比或Sc數(shù)的精確變化規(guī)律.但從實際工程的角度看,應(yīng)盡量減少振幅預(yù)測所需 Sc數(shù)的試驗工況,提高試驗效率.提高試驗效率的一個重要方法就是對最大豎向渦振振幅隨Sc數(shù)的變化曲線進行函數(shù)擬合.參考 Sc數(shù)對圓柱體渦振振幅的影響規(guī)律[26-27],采用兩種函數(shù)形 式對圖11所示試驗數(shù)據(jù)進行擬合,其中函數(shù)一為分數(shù)形式,函數(shù)二為冪指數(shù)形式,即

    式中:a、b和c表示函數(shù)一的待擬合參數(shù);a1和b1表示 函數(shù)二的待擬合參數(shù),從物理意義上看,a和a1表示當Sc數(shù)等于零時結(jié)構(gòu)可能達到的極限渦振振幅;ηmax表示主梁無量綱渦振振幅.

    不同風(fēng)攻角下函數(shù)一和函數(shù)二的參數(shù)識別結(jié)果見表1,擬合函數(shù)和試驗數(shù)據(jù)的對比如圖12所示.從圖12可看到,在試驗的Sc數(shù)范圍內(nèi),函數(shù)一和函數(shù)二均能較好地擬合主梁的最大渦振振幅隨Sc數(shù)的變化規(guī)律,但當Sc數(shù)小于試驗范圍時,函數(shù)一和函數(shù)二的擬合結(jié)果偏差逐漸增大.從表1可看到,對于函數(shù)一,3個風(fēng)攻角下的極限渦振振幅基本相同,均在0.135 左右;對于函數(shù)二,從-3°風(fēng)攻角到+3°風(fēng)攻角,對應(yīng)的極限渦振振幅逐漸減小.總體來看,函數(shù)二預(yù) 測的極限渦振振幅大于函數(shù)一.

    選擇合適的試驗方案可在保證擬合精度的前提下進一步減少工作量,下面以+3°風(fēng)攻角為例,說明 試驗方案對函數(shù)擬合精度的影響,結(jié)果如圖13(a)~(d)所示.對于函數(shù)一:方案一的擬合誤差最大,說明當Sc數(shù)小于試驗范圍時,函數(shù)一難以準確描述主梁渦振振幅隨Sc數(shù)的變化規(guī)律.相反,若用Sc數(shù)較小時的試驗數(shù)據(jù)預(yù)測高Sc數(shù)的主梁渦振響應(yīng),則相對誤差較?。ㄒ妶D13(b)).試驗方案三和試驗方案四都保留了Sc數(shù)最小和最大的試驗工況,然后再選取 Sc數(shù)介于二者之間的一個或兩個工況.從圖13(c)和圖13(d)可發(fā)現(xiàn),上述兩種方案的擬合誤差都較小,可有效預(yù)測不同 Sc數(shù)下的主梁渦振振幅.

    對比函數(shù)一和函數(shù)二的擬合結(jié)果可發(fā)現(xiàn),函數(shù)二能更精確地描述節(jié)段模型的最大渦振振幅隨Sc數(shù)的變化規(guī)律,即使要預(yù)測的Sc數(shù)小于試驗范圍,函數(shù)二的預(yù)測結(jié)果也非常準確,因此函數(shù)二將用于下文實橋渦振響應(yīng)預(yù)測研究.值得說明的是,函數(shù)二的適用性還有待針對不同的橋梁斷面開展進一步研究.

    3實橋主梁高階模態(tài)渦振響應(yīng)預(yù)測

    基于節(jié)段模型的最大渦振振幅預(yù)測公式(2),可進一步得到實橋各階豎 彎模態(tài)的渦振振幅預(yù)測公式:

    式中:ηi,max表示某個風(fēng)攻角下實橋第i階豎彎模態(tài)的最大渦振振幅;λL表示幾何縮尺比;λi,s表示第i階豎 彎模態(tài)的振型修正系數(shù);λi,c表示第i階豎彎模態(tài)的渦激力展向相關(guān)性修正系數(shù).其他參數(shù)的物理意義 如前文所示.

    需要說明的是,主梁各階豎彎模態(tài)的振型修正系數(shù)和渦激力展向相關(guān)性修正系數(shù)都與渦激力的數(shù)學(xué)模型和振型函數(shù)有關(guān),目前尚缺乏定論.本文從工 程實用的角度考慮,各階豎彎模態(tài)的振型修正系數(shù)都按 Scanlan非線性渦激力模型和簡諧波振型計算,

    即λi,s=2/3≈1.15;渦激力展向相關(guān)性修正系數(shù)按 完全相關(guān)考慮,λi,c=1.

    通過結(jié)構(gòu)有限元分析可得到大橋每延米的等效質(zhì)量,各階豎彎模態(tài)的阻尼比則參照英國抗風(fēng)規(guī)范取為0.48%,由此計算得到大橋前15階模態(tài)的Sc數(shù)如圖14所示.可以看到,大橋各階豎彎模態(tài)的Sc數(shù)基本相同,都在19左右.利用公式(3)可預(yù)測+3°風(fēng)攻角和0°風(fēng)攻角下,主梁各階豎彎模態(tài)的最大渦振振 幅,結(jié)果如圖15所示.圖中還依據(jù)我國《公路橋梁抗 風(fēng)設(shè)計規(guī)范》(JTG/T3360-01—2018)[9]給出了各階豎彎模態(tài)的渦振振幅容許值.可以看到,各階豎彎模態(tài)的最大渦振振幅基本相等,但對應(yīng)的規(guī)范容許值則隨著豎彎頻率的增加按反比例關(guān)系逐漸減小.因 此,高階模態(tài)的渦振振幅比低階模態(tài)更容易超過規(guī)范容許值,而且模態(tài)階次越高,可能超過規(guī)范容許值的程度也越大.

    綜上,對于大跨度懸索橋必須檢驗高階豎彎模態(tài)的渦振振幅是否滿足規(guī)范要求.僅對一階豎彎模態(tài)進行最大渦振振幅檢驗,一般是偏于不安全的.

    4結(jié)論

    1)永磁式板式電渦流阻尼器可為節(jié)段模型系統(tǒng) 提供連續(xù)可調(diào)的、精確的線性黏滯阻尼,為節(jié)段模型 渦振試驗的阻尼比調(diào)節(jié)提供了新方法.

    2)不同風(fēng)攻角下的主梁豎向渦振峰值都隨Sc數(shù)的增大而減小,但風(fēng)攻角不同,變化的規(guī)律也存在顯 著差異.總體來看,冪指數(shù)函數(shù)能較好地反映主梁豎向渦振峰值隨Sc數(shù)的變化規(guī)律.

    3)大跨度懸索橋各階豎彎模態(tài)的渦振振幅基本相等,高階模態(tài)的渦振響應(yīng)更容易超過規(guī)范限值,應(yīng)予以重視.

    參考文獻

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

    129-138.

    CHEN Z Q.Wind engineering of bridge[M].Beijing:China Com-munications Press,2005:129-138.(In Chinese)

    [2]EHSAN F,SCANLAN R H.Vortex-induced vibrations of flexible

    bridges[J].Journal of Engineering Mechanics,1990,116(6):1392-1411.

    [3]SIMIU E,SCANLAN R H.Wind effects on structures:fundamen-tals and applications to design [M].3nd ed.New York:Wiley,1996:217-230.

    [4]ZDRAVKOVICH M M.Scruton number:a proposal[J].Journal ofWind Engineering and Industrial Aerodynamics,1982,10(3): 263-265.

    [5]SCRUTON C.On the wind-excited oscillations of stacks,towers

    and masts [C]//Proceedings of the International Conference on the Wind Effects on Buildings and Structures.Teddington,UK,1963:798-837.

    [6]MARRA A M,MANNINI C,BARTOLI G.Van der Pol-type equa-tion for modeling vortex-induced oscillations of bridge decks[J].Journal of Wind Engineering and Industrial Aerodynamics,2011,99(6/7):776-785.

    [7]郭震山,陳艾榮,項海帆.大跨纜索承重橋梁的結(jié)構(gòu)阻尼特性[J].世界地震工程,2000,16(3):52-57.

    GUO Z S,CHEN A R,XIANG H F.Characteristics of structural damping of cable-supported bridges with long spans[J].World In-formation on Earthquake Engineering,2000,16(3):52-57.(In Chinese)

    [8]WANG S Q,ZHAO L,GE Y J.Mode-dependent structural damp-ing of suspension bridge [C]//The 2016 World Congress on Ad-vances in Civil,Environmental,and Materials Research.Jeju Is-land,Korea,2016:1-8.

    [9]公路橋梁抗風(fēng)設(shè)計規(guī)范:JTG/T3360-01—2018[S].北京:人民交通出版社,2019:55-56.

    Wind-resistant design specification for highway bridges:JTG/T3360-01—2018 [S].Beijing:?? China Communications Press,2019:55-56.(In Chinese)

    [10]Design rues for aerodynamic effects on bridges:BD 49/01[S].

    Dublin:National Roads Authority,2014.

    [11]葛耀君,趙林,許坤.大跨橋梁主梁渦激振動研究進展與思考[J].中國公路學(xué)報,2019,32(10):1-18.

    GE Y J,ZHAO L,XU K.Review and reflection on vortex-induced vibration of main girders of long-span bridges[J].China Journal of Highway and Transport,2019,32(10):1-18.(In Chinese)

    [12]許福友,丁威,姜峰,等.大跨度橋梁渦激振動研究進展與展望[J].振動與沖擊,2010,29(10):40-49.

    XU F Y,DING W,JIANG F,et al.Development and prospect of study on vortex-induced vibration of long-span bridges[J].Jour-nal of Vibration and Shock,2010,29(10):40-49.(In Chinese)

    [13]KLAMO J T,LEONARD A,ROSHKO A.The effects of damping

    on the amplitude and frequency response of a freely vibrating cylin-der in cross-flow[J].Journal of Fluids and Structures,2006,22(6/ 7):845-856.

    [14]LARSEN A.A generalized model for assessment of vortex-induced vibrations of flexible structures[J].Journal of Wind Engi-neering and Industrial Aerodynamics,1995,57(2/3):281-294.

    [15]李永君,葛耀君,杜柏松.大跨度橋梁質(zhì)量阻尼參數(shù)對渦激振動的影響[C]//第十六屆全國橋梁學(xué)術(shù)會議論文集.北京:人民交通出版社,2004:495-520.

    LI Y J,GE Y J,DU B S.Effect of mass-damping parameters on vortex-induced vibration of long-span bridges[C]//Proceedings of the16th National Conference on Bridges.Beijing:China Commu-nications Press,2004:495-520.(In Chinese)

    [16]MARRA A M,MANNINI C,BARTOLI G.Measurements and im-proved model of vortex-induced vibration for an elongated rectan-gular cylinder [J].Journal of Wind Engineering and Industrial Aerodynamics,2015,147:358-367.

    [17]ZHU L D,MENG X L,GUO Z S.Nonlinear mathematical model of

    vortex-induced vertical force on a flat closed-box bridge deck[J].Journal of Wind Engineering and Industrial Aerodynamics,2013,122:69-82.

    [18]XU K,GE Y J,ZHAO L,et al.Simulation of vortex-induced vibra-tion of long-span bridges:a nonlinear normal mode approach[J].International Journal of Structural Stability and Dynamics,2018,18(11):1850136.

    [19]ZHANG M J,WU T,XU F Y.Vortex-induced vibration of bridge

    decks:describing function-based model[J].Journal of Wind Engi-neering and Industrial Aerodynamics,2019,195:104016.

    [20]VANDIVER J K.Damping parameters for flow-induced vibration

    [J].Journal of Fluids and Structures,2012,35:105-119.

    [21]GARCIA E M H,BERNITSAS M M.Effect of damping on variable

    added mass and lift of circular cylinders in vortex-induced vibra-tions[J].Journal of Fluids and Structures,2018,80:451-472.

    [22]GAO G Z,ZHU L D.Nonlinearity of mechanical damping and stiff-ness of a spring-suspended sectional model system for wind tunnel tests[J].Journal of Sound and Vibration,2015,355:369-391.

    [23]CAO F C,GE Y J.Air-induced nonlinear damping and addedmass of vertically vibrating bridge deck section models under zero wind speed[J].Journal of Wind Engineering and Industrial Aero-dynamics,2017,169:217-231.

    [24]XU F Y,ZHANG Z B.Numerical simulation of windless-air-induced added mass and damping of vibrating bridge decks [J].Journal of Wind Engineering and Industrial Aerodynamics,2018,180:98-107.

    [25]ZHANG Z B,XU F Y.Added mass and damping effects on vibrat-ing bridge decks in still air[J].Journal of Wind Engineering and Industrial Aerodynamics,2019,191:227-238.

    [26]SARPKAYA T.A critical review of the intrinsic nature of vortex-induced vibrations[J].Journal of Fluids and Structures,2004,19(4):389-447.

    [27]SKOP R A,BALASUBRAMANIAN S.A new twist on an oldmodel for vortex-excited vibrations [J].Journal of Fluids and Structures,1997,11(4):395-412.

    猜你喜歡
    渦振阻尼比橋梁
    基于細觀結(jié)構(gòu)的原狀黃土動彈性模量和阻尼比試驗研究
    地震研究(2021年1期)2021-04-13 01:05:24
    分體式雙箱梁渦振氣動控制措施數(shù)值模擬
    板桁結(jié)合梁渦振性能及抑振措施研究
    手拉手 共搭愛的橋梁
    句子也需要橋梁
    黏滯阻尼器在時程分析下的附加有效阻尼比研究
    振動與沖擊(2019年4期)2019-02-22 02:33:34
    大跨度懸索橋渦振風(fēng)洞試驗與現(xiàn)場實測比較
    波形分析法求解公路橋梁阻尼比的探討
    上海公路(2018年3期)2018-03-21 05:55:48
    結(jié)構(gòu)構(gòu)件阻尼比對大跨度懸索橋地震響應(yīng)的影響
    高性能砼在橋梁中的應(yīng)用
    国产97色在线日韩免费| 一级作爱视频免费观看| 啦啦啦 在线观看视频| 免费看a级黄色片| 日日夜夜操网爽| 色综合欧美亚洲国产小说| 久久久久九九精品影院| 国产成人精品在线电影| 久久天堂一区二区三区四区| 巨乳人妻的诱惑在线观看| 久久久国产欧美日韩av| 老司机靠b影院| 欧美不卡视频在线免费观看 | 精品久久久久久久久久免费视频| 免费在线观看亚洲国产| 亚洲精华国产精华精| 叶爱在线成人免费视频播放| 精品国产美女av久久久久小说| av天堂在线播放| 一区二区三区激情视频| 99久久99久久久精品蜜桃| 国产aⅴ精品一区二区三区波| 香蕉丝袜av| 男人舔女人的私密视频| 校园春色视频在线观看| 久久香蕉国产精品| 九色亚洲精品在线播放| 国产人伦9x9x在线观看| 国产成人精品在线电影| 老司机靠b影院| 免费无遮挡裸体视频| 中文字幕久久专区| 亚洲第一电影网av| 女人精品久久久久毛片| 一进一出抽搐动态| 99在线视频只有这里精品首页| 18美女黄网站色大片免费观看| 一夜夜www| 日韩大尺度精品在线看网址 | 亚洲av片天天在线观看| 久久 成人 亚洲| 俄罗斯特黄特色一大片| 日韩视频一区二区在线观看| 亚洲中文字幕一区二区三区有码在线看 | 免费高清视频大片| av中文乱码字幕在线| 天天躁夜夜躁狠狠躁躁| 又黄又爽又免费观看的视频| 午夜福利影视在线免费观看| 国内精品久久久久精免费| 成人18禁在线播放| 亚洲性夜色夜夜综合| 久久国产精品男人的天堂亚洲| 国产精品免费视频内射| 久久久久久久久中文| 亚洲国产中文字幕在线视频| 大型av网站在线播放| 亚洲性夜色夜夜综合| 亚洲熟妇熟女久久| 亚洲精品久久成人aⅴ小说| 黄色 视频免费看| 欧美不卡视频在线免费观看 | 天天躁夜夜躁狠狠躁躁| 国产精品1区2区在线观看.| 可以在线观看毛片的网站| 最新美女视频免费是黄的| 真人做人爱边吃奶动态| 亚洲午夜理论影院| 中文字幕高清在线视频| 亚洲中文字幕一区二区三区有码在线看 | 国产不卡一卡二| 成人特级黄色片久久久久久久| 久久青草综合色| 在线av久久热| 午夜福利免费观看在线| 麻豆成人av在线观看| aaaaa片日本免费| 久久香蕉国产精品| 亚洲成人久久性| svipshipincom国产片| 一级作爱视频免费观看| 中文字幕人妻丝袜一区二区| 国产视频一区二区在线看| 久久午夜综合久久蜜桃| 欧美成人一区二区免费高清观看 | 国产熟女午夜一区二区三区| 老司机在亚洲福利影院| 99在线人妻在线中文字幕| 亚洲国产中文字幕在线视频| 黄频高清免费视频| 中文字幕高清在线视频| 久久这里只有精品19| 又大又爽又粗| 免费少妇av软件| 午夜激情av网站| 人妻丰满熟妇av一区二区三区| 亚洲精品国产色婷婷电影| 亚洲第一欧美日韩一区二区三区| 精品国产一区二区久久| 亚洲人成电影免费在线| 最新在线观看一区二区三区| 老司机午夜十八禁免费视频| av片东京热男人的天堂| 久久久国产成人精品二区| 女警被强在线播放| 性欧美人与动物交配| 色综合欧美亚洲国产小说| 亚洲男人的天堂狠狠| 天堂影院成人在线观看| 久久午夜综合久久蜜桃| av片东京热男人的天堂| 欧美成狂野欧美在线观看| 亚洲国产精品成人综合色| 精品国产乱码久久久久久男人| 久久人妻福利社区极品人妻图片| 久久精品影院6| 午夜影院日韩av| 婷婷精品国产亚洲av在线| 久久婷婷成人综合色麻豆| 国产乱人伦免费视频| 啦啦啦韩国在线观看视频| 色尼玛亚洲综合影院| 日本撒尿小便嘘嘘汇集6| 成年人黄色毛片网站| 男人的好看免费观看在线视频 | 在线av久久热| 久99久视频精品免费| av视频在线观看入口| 色综合站精品国产| 国产午夜精品久久久久久| 久久精品人人爽人人爽视色| 久久久久久国产a免费观看| 午夜福利免费观看在线| 久久久精品欧美日韩精品| 国产一区二区在线av高清观看| 亚洲精品粉嫩美女一区| 久久中文字幕一级| 亚洲成av片中文字幕在线观看| 国产精品影院久久| 久热爱精品视频在线9| 又黄又粗又硬又大视频| 99国产精品一区二区三区| 不卡av一区二区三区| 99久久99久久久精品蜜桃| 一区二区三区国产精品乱码| 日韩欧美在线二视频| 成人18禁在线播放| 久久人人爽av亚洲精品天堂| 免费在线观看日本一区| 日韩一卡2卡3卡4卡2021年| 亚洲人成77777在线视频| 脱女人内裤的视频| 国产99久久九九免费精品| 大香蕉久久成人网| 中文亚洲av片在线观看爽| 免费在线观看日本一区| 制服人妻中文乱码| 国产精品 国内视频| 99精品欧美一区二区三区四区| 亚洲成人精品中文字幕电影| 怎么达到女性高潮| 亚洲专区字幕在线| 曰老女人黄片| 侵犯人妻中文字幕一二三四区| 国产免费男女视频| 久久国产亚洲av麻豆专区| 757午夜福利合集在线观看| 一级黄色大片毛片| 色在线成人网| 午夜福利一区二区在线看| 久久精品人人爽人人爽视色| 91国产中文字幕| 在线av久久热| 国产激情欧美一区二区| 在线观看免费视频日本深夜| 18禁观看日本| 色在线成人网| 亚洲国产看品久久| 日本欧美视频一区| 涩涩av久久男人的天堂| 午夜视频精品福利| 12—13女人毛片做爰片一| 国产视频一区二区在线看| 国产精品98久久久久久宅男小说| 国产精品1区2区在线观看.| 一级黄色大片毛片| 国产视频一区二区在线看| 动漫黄色视频在线观看| 国产成人欧美| 悠悠久久av| 亚洲av成人不卡在线观看播放网| 日韩欧美国产一区二区入口| 欧美成人午夜精品| 久久久久久免费高清国产稀缺| 午夜福利高清视频| 亚洲全国av大片| 久久欧美精品欧美久久欧美| 亚洲精品一区av在线观看| 亚洲av成人不卡在线观看播放网| av在线播放免费不卡| 久久中文看片网| 免费高清在线观看日韩| 一个人观看的视频www高清免费观看 | 国内精品久久久久久久电影| 首页视频小说图片口味搜索| 国产成人免费无遮挡视频| 欧美国产精品va在线观看不卡| 激情在线观看视频在线高清| 老熟妇乱子伦视频在线观看| 国产乱人伦免费视频| 午夜两性在线视频| 久久人人精品亚洲av| 久久天堂一区二区三区四区| 国产精品九九99| 亚洲精品中文字幕一二三四区| a在线观看视频网站| 国产精品美女特级片免费视频播放器 | 村上凉子中文字幕在线| www.999成人在线观看| 黄频高清免费视频| 别揉我奶头~嗯~啊~动态视频| 日韩欧美在线二视频| 窝窝影院91人妻| 手机成人av网站| 自线自在国产av| 国产精品久久久久久精品电影 | 老司机在亚洲福利影院| 亚洲av成人不卡在线观看播放网| 久久香蕉精品热| 侵犯人妻中文字幕一二三四区| 国产成+人综合+亚洲专区| 国产精品野战在线观看| 国产亚洲精品一区二区www| 日本在线视频免费播放| 久久午夜综合久久蜜桃| 精品久久蜜臀av无| 日日夜夜操网爽| 在线十欧美十亚洲十日本专区| 身体一侧抽搐| 欧美老熟妇乱子伦牲交| 啪啪无遮挡十八禁网站| 首页视频小说图片口味搜索| 十分钟在线观看高清视频www| 热99re8久久精品国产| 淫妇啪啪啪对白视频| 9热在线视频观看99| 国产野战对白在线观看| 男女床上黄色一级片免费看| 男女下面插进去视频免费观看| 精品国产一区二区久久| 午夜免费鲁丝| 在线观看66精品国产| 亚洲伊人色综图| 男女下面进入的视频免费午夜 | 高清在线国产一区| 麻豆一二三区av精品| 欧美激情久久久久久爽电影 | АⅤ资源中文在线天堂| 国产一区二区三区综合在线观看| 黄色视频不卡| 91av网站免费观看| 高清毛片免费观看视频网站| 久久久国产成人免费| 最新在线观看一区二区三区| 国产精品一区二区免费欧美| 女人爽到高潮嗷嗷叫在线视频| 午夜激情av网站| 男人舔女人的私密视频| 国产xxxxx性猛交| 欧美日本中文国产一区发布| 老鸭窝网址在线观看| 老司机午夜十八禁免费视频| 九色国产91popny在线| 久久国产亚洲av麻豆专区| 欧美不卡视频在线免费观看 | 99riav亚洲国产免费| 午夜老司机福利片| 欧洲精品卡2卡3卡4卡5卡区| netflix在线观看网站| 精品久久久精品久久久| 久久婷婷成人综合色麻豆| 亚洲视频免费观看视频| 国产成人精品久久二区二区91| 亚洲成国产人片在线观看| 亚洲专区字幕在线| 美女午夜性视频免费| 男女下面插进去视频免费观看| 亚洲欧美激情综合另类| 一级毛片精品| 激情视频va一区二区三区| 亚洲视频免费观看视频| 国产亚洲av嫩草精品影院| 亚洲国产欧美网| 欧美av亚洲av综合av国产av| 国产伦人伦偷精品视频| 国产av一区在线观看免费| 久久久久久免费高清国产稀缺| 一二三四社区在线视频社区8| 国产成人精品无人区| 青草久久国产| 久久草成人影院| 久久人人爽av亚洲精品天堂| 欧美绝顶高潮抽搐喷水| ponron亚洲| 男女下面进入的视频免费午夜 | 美国免费a级毛片| 国产成人系列免费观看| 十八禁网站免费在线| av超薄肉色丝袜交足视频| 日本vs欧美在线观看视频| 最新美女视频免费是黄的| 久热爱精品视频在线9| 亚洲自拍偷在线| 免费在线观看完整版高清| 亚洲一码二码三码区别大吗| 亚洲人成电影免费在线| 成人三级做爰电影| 91麻豆av在线| 亚洲av成人一区二区三| 欧美成人午夜精品| 国产成人系列免费观看| 色综合站精品国产| 波多野结衣一区麻豆| 亚洲天堂国产精品一区在线| 岛国在线观看网站| 国产精品综合久久久久久久免费 | 国产亚洲精品久久久久5区| 99精品久久久久人妻精品| 黑人欧美特级aaaaaa片| 一二三四在线观看免费中文在| 精品乱码久久久久久99久播| 国产成人一区二区三区免费视频网站| 丝袜美足系列| 久久婷婷成人综合色麻豆| 免费在线观看日本一区| 日本欧美视频一区| 91国产中文字幕| 亚洲av片天天在线观看| 91成人精品电影| 欧美成狂野欧美在线观看| 国产激情欧美一区二区| 黄色女人牲交| 免费在线观看黄色视频的| 欧美乱色亚洲激情| 少妇被粗大的猛进出69影院| 99riav亚洲国产免费| 757午夜福利合集在线观看| 亚洲精品国产精品久久久不卡| 欧美另类亚洲清纯唯美| 午夜视频精品福利| 色播亚洲综合网| 黄色毛片三级朝国网站| 午夜免费成人在线视频| 国产精品 国内视频| 多毛熟女@视频| 大型黄色视频在线免费观看| 69av精品久久久久久| 久久青草综合色| 亚洲精品在线观看二区| 男女之事视频高清在线观看| 欧美久久黑人一区二区| 一级毛片高清免费大全| 亚洲自拍偷在线| bbb黄色大片| 亚洲一区二区三区色噜噜| 男女下面插进去视频免费观看| 韩国av一区二区三区四区| 欧美午夜高清在线| 可以免费在线观看a视频的电影网站| 免费在线观看影片大全网站| 少妇裸体淫交视频免费看高清 | 熟女少妇亚洲综合色aaa.| 超碰成人久久| 国产精品98久久久久久宅男小说| 久久久国产欧美日韩av| 亚洲欧美一区二区三区黑人| 精品一区二区三区四区五区乱码| 亚洲情色 制服丝袜| 亚洲国产精品合色在线| 18禁观看日本| 一区福利在线观看| 精品免费久久久久久久清纯| 麻豆国产av国片精品| 日韩欧美在线二视频| 午夜日韩欧美国产| 国产精品免费视频内射| 国产高清激情床上av| 久久中文字幕一级| 女人被狂操c到高潮| 精品国产乱子伦一区二区三区| 国产亚洲精品av在线| 亚洲精品中文字幕在线视频| 啦啦啦免费观看视频1| 国产精品久久久人人做人人爽| 99国产综合亚洲精品| 禁无遮挡网站| 久久这里只有精品19| 国产精品一区二区免费欧美| 成人18禁高潮啪啪吃奶动态图| 大型黄色视频在线免费观看| 免费人成视频x8x8入口观看| 亚洲 欧美一区二区三区| 一级a爱片免费观看的视频| 两个人视频免费观看高清| 91精品国产国语对白视频| 久久人人精品亚洲av| 青草久久国产| 国产成+人综合+亚洲专区| 在线免费观看的www视频| 免费在线观看黄色视频的| 国产亚洲精品久久久久久毛片| 啦啦啦 在线观看视频| 久久草成人影院| 老司机在亚洲福利影院| 国产精品九九99| 成人亚洲精品av一区二区| 在线观看免费视频日本深夜| 亚洲性夜色夜夜综合| 日韩欧美在线二视频| 人人妻,人人澡人人爽秒播| 免费看美女性在线毛片视频| 正在播放国产对白刺激| 亚洲成人精品中文字幕电影| 大型av网站在线播放| 一区二区三区国产精品乱码| 亚洲精品一区av在线观看| 欧美日韩乱码在线| 9色porny在线观看| 最好的美女福利视频网| 午夜亚洲福利在线播放| 亚洲欧美日韩高清在线视频| 亚洲欧美日韩无卡精品| 十八禁人妻一区二区| 好看av亚洲va欧美ⅴa在| 欧美日本亚洲视频在线播放| 男人的好看免费观看在线视频 | 国产一区二区三区综合在线观看| 国产欧美日韩综合在线一区二区| 欧美日韩精品网址| 亚洲精品国产精品久久久不卡| www.精华液| 757午夜福利合集在线观看| 日韩欧美一区二区三区在线观看| 美女高潮到喷水免费观看| 欧美国产日韩亚洲一区| 在线视频色国产色| 久久性视频一级片| 又黄又爽又免费观看的视频| 欧美老熟妇乱子伦牲交| 国产亚洲欧美98| 欧美亚洲日本最大视频资源| 99国产精品一区二区蜜桃av| 日日爽夜夜爽网站| 欧美老熟妇乱子伦牲交| 熟女少妇亚洲综合色aaa.| 国产又色又爽无遮挡免费看| 亚洲欧美精品综合久久99| 精品人妻1区二区| 久久九九热精品免费| 国产亚洲精品一区二区www| 男女床上黄色一级片免费看| 一卡2卡三卡四卡精品乱码亚洲| 天天躁夜夜躁狠狠躁躁| 国产精品久久视频播放| 九色亚洲精品在线播放| 国产成年人精品一区二区| 精品欧美一区二区三区在线| 黑人巨大精品欧美一区二区蜜桃| 亚洲五月婷婷丁香| 成人国产综合亚洲| 亚洲一码二码三码区别大吗| 国产99久久九九免费精品| 亚洲成人国产一区在线观看| 免费看十八禁软件| 看黄色毛片网站| 99精品欧美一区二区三区四区| 中出人妻视频一区二区| 嫩草影视91久久| 18禁美女被吸乳视频| 两性午夜刺激爽爽歪歪视频在线观看 | 这个男人来自地球电影免费观看| 我的亚洲天堂| 国产一级毛片七仙女欲春2 | 日本免费一区二区三区高清不卡 | 国产亚洲欧美精品永久| 美女国产高潮福利片在线看| 1024香蕉在线观看| 老鸭窝网址在线观看| 亚洲视频免费观看视频| 中亚洲国语对白在线视频| 淫妇啪啪啪对白视频| 一区二区三区高清视频在线| 国产av一区二区精品久久| 女警被强在线播放| 国产成人欧美| 午夜福利在线观看吧| 欧美精品啪啪一区二区三区| 91成年电影在线观看| 亚洲精品中文字幕一二三四区| 国产av一区在线观看免费| 自线自在国产av| 欧美乱妇无乱码| 一级黄色大片毛片| 国产亚洲精品久久久久5区| 熟妇人妻久久中文字幕3abv| 精品国产乱码久久久久久男人| 日日爽夜夜爽网站| 国产午夜精品久久久久久| 国产精品电影一区二区三区| 午夜a级毛片| 性欧美人与动物交配| 99国产极品粉嫩在线观看| 亚洲一卡2卡3卡4卡5卡精品中文| 国产亚洲精品久久久久久毛片| 国产精品久久视频播放| 欧美绝顶高潮抽搐喷水| 亚洲精品中文字幕一二三四区| 精品高清国产在线一区| 午夜福利免费观看在线| 色av中文字幕| 欧美av亚洲av综合av国产av| 无人区码免费观看不卡| 麻豆久久精品国产亚洲av| 国产亚洲av嫩草精品影院| 亚洲欧美激情在线| 一本久久中文字幕| 黑人巨大精品欧美一区二区mp4| 日韩精品青青久久久久久| 久久久久久国产a免费观看| 亚洲狠狠婷婷综合久久图片| 国产区一区二久久| 一级a爱片免费观看的视频| 亚洲avbb在线观看| 少妇被粗大的猛进出69影院| 久久影院123| 欧美黑人欧美精品刺激| 日日爽夜夜爽网站| 国产精品爽爽va在线观看网站 | 久久香蕉激情| 伊人久久大香线蕉亚洲五| 老司机午夜十八禁免费视频| 亚洲少妇的诱惑av| 精品一区二区三区四区五区乱码| 色在线成人网| √禁漫天堂资源中文www| 国产精品香港三级国产av潘金莲| 成人三级做爰电影| 午夜影院日韩av| 婷婷丁香在线五月| 在线观看免费视频网站a站| 法律面前人人平等表现在哪些方面| 欧美黑人精品巨大| 精品一区二区三区av网在线观看| 国产精品日韩av在线免费观看 | 香蕉丝袜av| 99国产精品一区二区三区| 久久人妻熟女aⅴ| 亚洲一区中文字幕在线| 久久人妻熟女aⅴ| 欧美 亚洲 国产 日韩一| 亚洲精品中文字幕一二三四区| 久久久久国产精品人妻aⅴ院| 男女下面插进去视频免费观看| 少妇熟女aⅴ在线视频| 欧美成人性av电影在线观看| 欧美绝顶高潮抽搐喷水| 成在线人永久免费视频| 一本综合久久免费| 丝袜美腿诱惑在线| 亚洲一区中文字幕在线| av在线天堂中文字幕| 亚洲国产精品合色在线| 大型av网站在线播放| av免费在线观看网站| 中国美女看黄片| 午夜激情av网站| 久久久久国内视频| 亚洲三区欧美一区| a在线观看视频网站| 精品一区二区三区视频在线观看免费| 老鸭窝网址在线观看| 国产精品亚洲一级av第二区| 一区二区三区精品91| 亚洲人成网站在线播放欧美日韩| 久久性视频一级片| 黑人巨大精品欧美一区二区mp4| 一本久久中文字幕| 色哟哟哟哟哟哟| 一进一出抽搐动态| 女人高潮潮喷娇喘18禁视频| 在线观看日韩欧美| 午夜免费成人在线视频| 国产精品亚洲av一区麻豆| 97碰自拍视频| 亚洲成国产人片在线观看| 成人三级黄色视频| 91大片在线观看| 国产精品一区二区三区四区久久 | 国产麻豆成人av免费视频| av天堂在线播放| 精品欧美一区二区三区在线| 在线永久观看黄色视频| 亚洲在线自拍视频| 久久香蕉国产精品| 一本综合久久免费| 人成视频在线观看免费观看| 日韩中文字幕欧美一区二区| 在线观看午夜福利视频| 身体一侧抽搐| 精品无人区乱码1区二区| 好男人电影高清在线观看| 夜夜夜夜夜久久久久| 97人妻精品一区二区三区麻豆 | 国产成人影院久久av| 一区二区三区高清视频在线|