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

    中央扣對大跨懸索橋顫振穩(wěn)定性的影響

    2021-05-17 17:21:16李凱韓艷蔡春聲樊中武林超

    李凱 韓艷 蔡春聲 樊中武 林超

    摘 ? 要:為研究中央扣對大跨度懸索橋顫振穩(wěn)定性的影響,以矮寨大橋為工程背景,基于大橋精細(xì)化空間桁架梁有限元模型,根據(jù)主梁整體剛度等效原則,采用懸臂梁位移法建立了大橋等效單主梁有限元模型;考慮了跨中短吊桿(無中央扣)、一對柔性中央扣、三對柔性中央扣和剛性中央扣4種不同結(jié)構(gòu)形式,計算分析了中央扣對大跨度懸索橋自振特性的影響;基于試驗獲得的顫振導(dǎo)數(shù),采用脈沖響應(yīng)函數(shù)結(jié)合Roger有理函數(shù)并利用非線性最小二乘擬合方法擬合其系數(shù)從而得到主梁斷面自激力的時域表達(dá)式,隨后利用ANSYS二次開發(fā),實現(xiàn)了大橋顫振穩(wěn)定性時域分析,研究了中央扣對顫振臨界風(fēng)速、顫振頻率及全橋三維顫振姿態(tài)的影響規(guī)律. 結(jié)果表明:不論柔性還是剛性中央扣都能夠顯著提高主梁縱飄振型的振動頻率,其對反對稱側(cè)彎和反對稱扭轉(zhuǎn)頻率的影響比正對稱大;剛性中央扣能夠大幅提高反對稱扭轉(zhuǎn)振型的頻率. 由于矮寨大橋是以一階正對稱豎彎、二階正對稱豎彎和一階正對稱扭轉(zhuǎn)相互耦合的振型發(fā)生彎扭耦合顫振,因此,中央扣對顫振臨界風(fēng)速的影響極小,但對顫振頻率與主梁三維顫振姿態(tài)有一定影響,并一定程度上有利于顫振穩(wěn)定性. 此外還發(fā)現(xiàn)當(dāng)結(jié)構(gòu)阻尼很低時,由于顫振頻率落于固有頻率分布十分密集的區(qū)域,主梁顫振狀態(tài)有復(fù)雜拍振現(xiàn)象(間歇性顫振現(xiàn)象)的出現(xiàn).

    關(guān)鍵詞:中央扣;顫振穩(wěn)定性;大跨度橋梁;自激力;時域分析

    中圖分類號:U447 ? ? ? ? ? ? ? ? ? ? ? ? ? ? 文獻(xiàn)標(biāo)志碼:A

    Influence of Central Buckle on Flutter Stability

    of Long-span Suspension Bridge

    LI Kai1,HAN Yan1?,CAI Chunsheng1,2,F(xiàn)AN Zhongwu3,LIN Chao3

    (1. School of Civil Engineering,Changsha University of Science & Technology,Changsha 410114,China;

    2. Department of Civil and Environmental Engineering,Louisiana State University,Baton Rouge LA70803,USA;

    3. Bridge Management Office of Aizhai Bridge,Hunan Expressway Group Co Ltd,Jishou 416000,China)

    Abstract:To investigate the effects of central buckles on the flutter stability of long-span suspension bridges,the Aizhai Bridge in China was selected as the engineering background. Based on a refined spatial-truss-girder model and according to the principle of stiffness equivalence of the main girder in all directions,an equivalent-single-girder finite element model was firstly ?established by using the displacement method of a cantilever beam. Subsequently,four different connection options between the main cable and the girder near the mid-span position,namely,a short suspender,one pair of flexible central buckles,three pairs of flexible central buckles and one pair of rigid central buckles,were considered and their effects on the dynamic characteristics of long span suspension bridges were studied. Then,based on the flutter derivatives obtained from wind tunnel tests,the time domain formulations of self-excited forces in the girder section,expressed in terms of convolution integrals of impulse response functions,were obtained using a nonlinear least square fitting method. Based on APDL (ANSYS Parametric Design Language) offered by ANSYS,a time-domain flutter analysis was realized. Finally,the influences of the central buckles on the critical flutter velocity,flutter frequency,and three-dimensional flutter states of the bridge were investigated. The results show that the central buckles can significantly increase the frequency of the longitudinal floating mode of the bridge and have greater influence on the frequencies of asymmetric lateral bending mode and asymmetric torsion mode than that of symmetric ones. The rigid central buckle can largely increase the frequency of asymmetric torsion mode. The central buckles have negligible impact on the critical flutter velocity because the flutter mode shape of the Aizhai Bridge is coupled with the symmetric vertical bending mode shape and the symmetric torsion mode shape. However,it has a certain impact on the flutter frequency and the three-dimensional flutter states of the bridge,which benefits the flutter stability. In addition,it is found that the phenomenon of complex beat vibration (called intermittent flutter phenomenon) appeared when the structural damping was very low,because the flutter frequency falls in an area where the natural frequency distribution is very dense.

    Key words:central buckle;flutter stability;long-span suspension bridge;self-excited force;time-domain analysis

    大跨懸索橋往往伴隨著柔度大、阻尼低等特點[1],很容易受到交通、風(fēng)及地震等動態(tài)環(huán)境荷載的影響,采取適當(dāng)?shù)拇胧﹣砜刂苿討B(tài)載荷作用下的結(jié)構(gòu)響應(yīng)是其最重要的問題之一. 此外,對于典型的不設(shè)置中央扣的懸索橋,在動力荷載作用下主纜與加勁梁的位移不同步甚至反向,時常導(dǎo)致懸索橋短吊桿彎曲及斷裂的現(xiàn)象發(fā)生. 因此,在懸索橋的主纜和加勁梁跨中處設(shè)置中央扣聯(lián)結(jié)開始廣泛應(yīng)用于大跨懸索橋中,它能夠在一定程度上防止跨中短吊桿的彎曲和斷裂,抑制結(jié)構(gòu)動態(tài)位移響應(yīng),提高結(jié)構(gòu)的力學(xué)性能. 另外,中央扣的設(shè)計和使用可以追溯到New Tacoma Bridge的建設(shè)[2],中央扣的使用能提高橋梁的抗風(fēng)穩(wěn)定性. 正是由于中央扣具有以上優(yōu)點且構(gòu)造簡單,使其成為大跨度懸索橋的一種發(fā)展趨勢,但中央扣對顫振性能的影響機(jī)制仍是不明確的,因此探究中央扣的設(shè)置方式對懸索橋顫振穩(wěn)定性的影響能夠為日后中央扣在大跨度懸索橋中的廣泛應(yīng)用提供參考價值,具有重要意義.

    近年來,我國已有多座懸索橋采用中央扣結(jié)構(gòu),相關(guān)的文獻(xiàn)報道多集中于中央扣對懸索橋固有模態(tài)和抗震性能的改變,而關(guān)于中央扣對懸索橋抗風(fēng)穩(wěn)定性的影響則鮮有報道. 徐勛等[3-4]等基于有限元分析發(fā)現(xiàn)中央扣能夠明顯提高懸索橋一階反對稱扭轉(zhuǎn)和縱飄振型的頻率. Qin等[5]也采用有限元法得到了與上述類似的結(jié)論. 胡騰飛等[6]通過矮寨大橋氣彈模型動力測試研究發(fā)現(xiàn)反對稱側(cè)彎頻率隨中央扣剛度的增大而增加,中央扣能夠提高懸索橋的整體剛度,其中以反對稱扭轉(zhuǎn)和縱向振動頻率提升最為顯著,進(jìn)一步驗證了上述基于有限元分析的部分結(jié)論. 此外,彭旺虎等[7]通過理論研究發(fā)現(xiàn)主纜彈性剛度與主纜重力剛度、加勁梁剛度之和的比值決定了第一階對稱扭轉(zhuǎn)和反對稱扭轉(zhuǎn)振型的形態(tài)和排列次序,中央扣剛度和跨中纜梁間距是影響反對稱扭轉(zhuǎn)基頻大小的主要因素. 在中央扣對懸索橋地震響應(yīng)的影響研究方面,Wang[8]和徐勛[4]等發(fā)現(xiàn)中央扣可以減小由地震激發(fā)的主梁端部的縱向位移但會增加塔頂?shù)目v向位移及塔底的彎矩. 另外,Wang等[9-10]基于車橋耦合振動分析發(fā)現(xiàn)中央扣對行車舒適性指數(shù)幾乎沒有影響. Wang等[11-12]通過抖振計算分析發(fā)現(xiàn)當(dāng)風(fēng)速為10 ~ 40 m/s時,中央扣能顯著減小主梁抖振豎向位移和主塔的縱向位移. 然而,采用時域分析方法研究中央扣對大跨度懸索橋顫振穩(wěn)定性的影響,國內(nèi)外鮮有文章報道.

    本文以矮寨大橋為工程背景,采用懸臂梁位移法建立了大橋等效單主梁有限元模型;考慮了跨中短吊桿(無中央扣)、一對柔性中央扣、三對柔性中央扣和剛性中央扣4種結(jié)構(gòu)形式,計算分析了中央扣對大跨度懸索橋自振特性的影響;基于試驗獲得的顫振導(dǎo)數(shù),采用脈沖響應(yīng)函數(shù)結(jié)合Roger有理函數(shù)并利用非線性最小二乘擬合方法擬合其系數(shù)從而得到主梁斷面自激力的時域表達(dá)式;隨后利用ANSYS二次開發(fā),實現(xiàn)了大橋顫振穩(wěn)定性時域分析,研究了中央扣對顫振臨界風(fēng)速、顫振頻率及主梁三維顫振姿態(tài)的影響規(guī)律. 研究結(jié)果可為日后中央扣在大跨度懸索橋中的廣泛應(yīng)用提供參考.

    1 ? 有限元建模及動力特性分析

    1.1 ? 工程背景

    矮寨大橋是一座跨越深山峽谷的單跨懸索橋,主纜矢跨比1 ∶ 9.6,兩主纜形心間距為27 m,主纜孔跨布置為242 m+1 176 m+116 m,加勁梁全長1 000.5 m,主跨跨中附近設(shè)置三對柔性中央扣,其中中央扣關(guān)于主梁跨中不對稱,如圖1(a)所示. 主梁為鋼桁加勁梁,桁寬27 m,桁高7.5 m. 鋼桁架主梁上采用鋼縱梁與混凝土橋面板相結(jié)合的橋面體系,橫斷面布置如圖1(b)所示.

    1.2 ? 有限元模型

    為了準(zhǔn)確獲得矮寨大橋的動力特性,首先采用有限元計算軟件ANSYS建立大橋空間桁架主梁有限元模型,其中鋼桁加勁梁與主塔采用beam188梁單元模擬,主纜和吊桿采用Link10桿單元模擬并設(shè)置為僅有拉伸剛度. 鋼-混組合橋面體系采用梁-殼單元混合有限元建模的方法建模,其中混凝土橋面板采用Shell63殼單元模擬,鋼縱梁采用Beam188梁單元模擬,同時對梁-殼單元采用MPC方法進(jìn)行耦合連接. 有限元模型及其笛卡爾坐標(biāo)系如圖2所示. 由其計算得到的矮寨大橋動力特性列于表1中. 為大幅提高顫振時域分析的計算效率并方便準(zhǔn)確地對有限元模型施加自激力荷載,基于主梁各方向整體剛度等效的原則,采用懸臂梁單位荷載位移法[13]建立了等效單主梁有限元模型(即魚骨梁模型). 空間桁架主梁模型與等效單主梁模型動力特性結(jié)果對比如表1所示,兩者主要相應(yīng)振型的頻率誤差基本都在1%以內(nèi),僅一階反對稱側(cè)彎振型頻率誤差約為1.7%,因此可以采用等效單主梁模型進(jìn)行顫振時域分析.

    1.3 ? 中央扣對動力特性的影響

    為研究不同結(jié)構(gòu)形式的中央扣對懸索橋固有模態(tài)及顫振穩(wěn)定性的影響,按跨中處主纜和加勁梁的不同連接形式定義了如表2所示的4種有限元計算模型. 其中FM-A的短吊桿、FM-B與FM-C的柔性中央扣都采用LINK10單元模擬,并設(shè)置為僅有拉伸剛度. FM-D的剛性中央扣參考文獻(xiàn)[3]中同為鋼桁懸索橋的四渡河大橋設(shè)計尺寸,如圖3所示. 剛性中央扣采用Beam188單元模擬. 同時為了考慮中央扣索夾和主纜構(gòu)成的局部位置能夠向剛性中央扣及主梁傳遞彎矩的效應(yīng),模型中對此處采用剛臂梁單元模擬,如圖4所示. 對比FM-B、FM-C與FM-D可知剛性扣與柔性扣對大跨懸索橋影響的差異. 模態(tài)分析得到的4種有限元模型的主要自振頻率及振型如表3所示.

    由表3中4種有限元模型自振頻率與振型特征的分析比較可以發(fā)現(xiàn):1)由于中央扣的傾斜拉桿能夠一定程度限制主纜與加勁梁之間的相對運(yùn)動,相當(dāng)于加勁梁在跨中位置被施加了縱向位移約束進(jìn)而顯著提高了主梁的縱向剛度,因此FM-A的縱飄振型伴隨一階反對稱豎彎振型出現(xiàn),頻率為0.102 99 Hz;而FM-B、FM-C、FM-D的縱飄振型都伴隨三階反對稱豎彎振型出現(xiàn),頻率分別為0.355 2 Hz、0.413 8 Hz、0.421 6 Hz,由此可見中央扣能夠大幅提高大跨度懸索橋縱飄振型的頻率并提高與其耦合的豎彎振型的階數(shù),且剛性中央扣的提升幅度最大. 2)相比于正對稱側(cè)彎振型的頻率,反對稱側(cè)彎的頻率受中央扣的影響明顯偏大. ?3)豎彎振型的頻率受中央扣的影響很小主要是由于中央扣僅提高了主梁跨中位置的局部豎向剛度. ?4)正對稱扭轉(zhuǎn)振型頻率受中央扣的影響較小,其中剛性中央扣能夠提高3%一階正對稱扭轉(zhuǎn)振型的頻率. 3種中央扣都能夠大幅度提高反對稱扭轉(zhuǎn)振型的頻率,以一階反對稱扭轉(zhuǎn)振型為例,相比FM-A,有中央扣結(jié)構(gòu)的FM-B、FM-C、FM-D頻率分別提高了5.9%、10.6%、14.8%,可見剛性中央扣提升幅度更大,因此剛性中央扣的抗風(fēng)穩(wěn)定性能更佳.

    2 ? 顫振穩(wěn)定性分析

    從上一節(jié)的分析中發(fā)現(xiàn)中央扣對反對稱振型的影響要明顯大于正對稱振型,但其他模態(tài)對顫振的貢獻(xiàn)需進(jìn)一步分析,且中央扣對于整個結(jié)構(gòu)的非線性尤其是幾何非線性[14]的影響是未可知的. 因此,中央扣對柔性主梁顫振行為的實際影響需進(jìn)一步采用顫振時域分析的方法定性定量地探討.

    2.1 ? 顫振時域分析方法

    橋梁結(jié)構(gòu)在均勻流中的運(yùn)動方程可以描述為:

    M[q] + C[q]+ Kq = Fae ? ? ? ?(1)

    式中:M、C和K分別為結(jié)構(gòu)質(zhì)量、阻尼和剛度矩陣; [q]、[q] 和q分別代表節(jié)點加速度向量、速度向量和位移向量;Fae表示橋梁受到的氣動荷載向量,當(dāng)進(jìn)行顫振穩(wěn)定性分析時,可只考慮氣動自激力. 采用脈沖響應(yīng)函數(shù)表達(dá)的單位長度斷面自激升力和自激扭轉(zhuǎn)(忽略側(cè)向位移的影響)分別為[15]:

    Lse(t) = Lseh(t) + Lseα(t) =

    ρU2(ILh(t-τ)h(τ)+ILα(t-τ)α(τ))dτ

    (2)

    Mse(t) = Mseh(t) + Mseα(t) =

    ρU2(IMh(t-τ)h(τ)+IMα(t-τ)α(τ))dτ

    (3)

    式中:Lsex(x = h、α)代表相應(yīng)位移產(chǎn)生的氣動升力;Msex(x = h、α)代表相應(yīng)位移產(chǎn)生的氣動扭矩;Ifx(f = L、M,x = h、α)為脈沖響應(yīng)函數(shù). 以豎向位移產(chǎn)生的氣動升力為例(忽略由第三項加速度項引起的氣動力),其最終純時域表達(dá)式為[16]:

    Lseh(t) = ρU2A1h(t) + A2

    [h] (t) +

    ρU2Ane[h] (t)dτ ?(4)

    式中:A1、A2、An、dn(n = 4,5,…,m)為用近似連續(xù)函數(shù)Roger有理函數(shù)連續(xù)表達(dá)離散的顫振導(dǎo)數(shù)時待擬合的系數(shù)[16],可采用最小二乘法擬合確定[17],本文不再贅述. 矮寨大橋0°風(fēng)攻角下節(jié)段模型顫振導(dǎo)數(shù)試驗值和用有理函數(shù)近似連續(xù)表達(dá)反算的顫振導(dǎo)數(shù)擬合值對比如圖5所示,可見試驗值與擬合值吻合得很好.

    基于自激力時域表達(dá)式通過ANSYS APDL語言編寫自激力數(shù)值計算程序,實現(xiàn)大橋顫振時域分析,其中式(4)中的卷積積分項在計算程序中需轉(zhuǎn)化為數(shù)值遞推關(guān)系式. 計算程序采用Newmark-β法求解式(1)運(yùn)動方程,通過時間步長無關(guān)性檢驗后計算時間步長取0.01 s. 最后,通過逐步提高風(fēng)速,觀察節(jié)點位移隨時間變化的特性,便可得到結(jié)構(gòu)顫振臨界風(fēng)速.

    2.2 ? 不考慮結(jié)構(gòu)阻尼

    鋼桁梁橋的結(jié)構(gòu)阻尼比一般取0.5%[18]. 然而,現(xiàn)存很多鋼橋的實測結(jié)構(gòu)阻尼比小于0.5%,甚至只有0.2%~0.3%左右. 因此,為充分研究中央扣對大跨懸索橋顫振穩(wěn)定性的影響同時比較結(jié)構(gòu)阻尼對顫振臨界風(fēng)速的影響,本文對不考慮結(jié)構(gòu)阻尼和考慮結(jié)構(gòu)阻尼進(jìn)行了分析研究.

    圖6~圖9為不考慮結(jié)構(gòu)阻尼0°攻角情況下模型FM-C跨中節(jié)點位移時程. 可知風(fēng)速74 m/s時豎向和扭轉(zhuǎn)位移整體上既無發(fā)散也無收斂的趨勢;而風(fēng)速78 m/s時豎向和扭轉(zhuǎn)位移整體上都有明顯的發(fā)散趨勢. 細(xì)化計算風(fēng)速范圍74~78 m/s時,發(fā)現(xiàn)74 m/s可被認(rèn)為是該條件下的顫振臨界風(fēng)速. 另由圖6~圖9還可以發(fā)現(xiàn)FM-C顫振狀態(tài)具有明顯多頻率參與的拍振特性,稱為間歇性顫振現(xiàn)象.

    圖10、11分別為風(fēng)速74 m/s時豎向和扭轉(zhuǎn)位移響應(yīng)的頻譜分析圖,可知豎向位移和扭轉(zhuǎn)位移響應(yīng)具有3個一致的卓越頻率. 其原因在于該橋在一階正對稱扭轉(zhuǎn)頻率附近,密集分布著一些側(cè)彎、豎彎、扭轉(zhuǎn)相互耦合的振型,如圖12所示. 當(dāng)主梁發(fā)生彎扭耦合顫振時,其顫振頻率為0.262 9 Hz,落于上述頻率分布密集區(qū),當(dāng)結(jié)構(gòu)阻尼為0時,極易牽連帶動頻率分布密集區(qū)的振型間歇性參與振動. 此外,由圖12可知FM-A、FM-B和FM-D的扭轉(zhuǎn)基頻附近頻率分布也很密集,因此它們的顫振狀態(tài)也有多頻率參與的拍振現(xiàn)象.

    表4為4種不同模型顫振臨界風(fēng)速與顫振頻率的對比,可知3種不同結(jié)構(gòu)形式的中央扣對顫振臨界風(fēng)速幾乎沒有影響.

    2.3 ? 考慮結(jié)構(gòu)阻尼

    本節(jié)首先探討結(jié)構(gòu)阻尼比對顫振臨界風(fēng)速以及上述間歇性顫振現(xiàn)象的影響. 圖13為顫振臨界風(fēng)速隨結(jié)構(gòu)阻尼比的變化趨勢,可以發(fā)現(xiàn)顫振臨界風(fēng)速并不隨結(jié)構(gòu)阻尼比成線性增長而是增長率不斷降低的趨勢. 此外,如圖14所示,當(dāng)結(jié)構(gòu)阻尼比較低時(ζ=0.05%、0.1%、0.15%),仍具有明顯的顫振拍現(xiàn)象,由此可見低阻尼情況下同樣存在顫振拍現(xiàn)象. 顫振拍現(xiàn)象隨結(jié)構(gòu)阻尼比的增長而不斷弱化,如圖15和16所示,當(dāng)結(jié)構(gòu)阻尼比為0.5%時顫振拍現(xiàn)象很快消失并達(dá)到單頻穩(wěn)幅的振動狀態(tài). 另外,可以預(yù)見,隨著跨度的增長,主梁的密頻和低阻尼特性將越明顯,因此越易產(chǎn)生間歇性顫振現(xiàn)象,這對主梁的顫振穩(wěn)定性是有利的. 這些在顫振頻率附近的模態(tài)間歇性參與顫振能不斷消耗斷面從氣流中吸收的能量,其作用相當(dāng)于一個間歇性調(diào)諧質(zhì)量阻尼器(TMD). 該TMD的阻尼效果隨著結(jié)構(gòu)阻尼的增長而不斷減小,因此可以預(yù)測其為顫振臨界風(fēng)速增長率隨結(jié)構(gòu)阻尼比增長而降低的一個重要原因,其具體效應(yīng)有待進(jìn)一步驗證,因為顫振臨界風(fēng)速隨結(jié)構(gòu)阻尼的這種變化趨勢還與氣動阻尼隨風(fēng)速的非線性變化有關(guān).

    接下來詳細(xì)探討結(jié)構(gòu)阻尼比為0.5%時,中央扣對顫振臨界風(fēng)速以及主梁三維顫振形態(tài)的影響. 圖15和圖16分別為風(fēng)速94.6 m/s時FM-C跨中節(jié)點豎向位移和扭轉(zhuǎn)位移時程,易知94.6 m/s為顫振臨界風(fēng)速,此值與風(fēng)洞試驗結(jié)果95.1 m/s[19]接近,說明程序計算結(jié)果是可信的.

    對圖15后500 s的位移時程做頻譜分析,結(jié)果如圖17所示.豎向和扭轉(zhuǎn)位移都只有一個卓越頻率0.257 8 Hz,原因是正結(jié)構(gòu)正阻尼下,頻率密集區(qū)的固有模態(tài)的位移響應(yīng)迅速衰減,或顫振頻率很難牽連帶動頻率密集區(qū)的固有模態(tài)間歇性參與振動. 表5為4種模型顫振臨界風(fēng)速與顫振頻率的對比,可知3種不同中央扣對顫振臨界風(fēng)速和顫振頻率幾乎無影響. 為探明其機(jī)理,基于主梁所有節(jié)點的位移時程提取出顫振臨界時的主梁振動形態(tài),如圖18所示. 可見豎向振動的位移最大值在1/4跨位置,而主梁一階對稱豎彎的最大值在跨中,二階正對稱豎彎的最大值在1/4跨,如圖19所示. 因此其顫振振型由一階正對稱豎彎、二階正對稱豎彎、一階正對稱扭轉(zhuǎn)振型相互耦合而成,但中央扣對這3種模態(tài)的頻率幾乎沒有影響,從而導(dǎo)致3種不同中央扣對顫振臨界風(fēng)速和顫振頻率幾乎沒有影響.

    如圖20和圖21所示,當(dāng)主梁發(fā)生彎扭耦合顫振時,其展向不同位置的位移響應(yīng)存在相位差且同一位置處豎向位移與扭轉(zhuǎn)位移也存在相位差. 主梁不同位置之間的振動相位差即可定量揭示主梁振動時的展向同步性,同步性越高,能量越容易累積,顫振穩(wěn)定性就越差,因此相位差的分析可以一定程度上揭示顫振穩(wěn)定性的機(jī)理,同時可以通過4個模型的對比分析探究中央扣對主梁顫振展向同步性的影響. 采用信號處理的方法可獲得兩個同頻信號之間的相位差,圖22為主梁其他位置與其跨中位置豎向位移相位差沿跨長的變化規(guī)律,可以發(fā)現(xiàn)模型FM-A的主梁其他位置與其跨中位置豎向位移相位差隨著二者距離的增大而增大,最大達(dá)到25°. 對比4個模型,可以發(fā)現(xiàn)關(guān)于跨中設(shè)置不對稱的中央扣(非對稱型的中央扣,如圖4所示)打亂了局部位置處相位差隨距離增大而增大的規(guī)律,并增大了主梁左邊與跨中處的豎向位移相位差,最大達(dá)40°,中央扣的不同形式造成的影響差別極小. 整體上非對稱型的中央扣有打亂主梁豎向位移展向同步性的趨勢,這對顫振穩(wěn)定性是有利的.

    圖23為主梁其他位置與其跨中位置扭轉(zhuǎn)位移相位差沿跨長的變化規(guī)律,可以發(fā)現(xiàn)主梁其他位置與其跨中位置扭轉(zhuǎn)位移相位差也隨二者距離的增大而增大,但最大僅3°,說明主梁顫振時扭轉(zhuǎn)位移展向同步性明顯高于豎向位移. 中央扣對主梁其他位置與其跨中位置扭轉(zhuǎn)位移相位差的影響微弱.

    圖24為主梁同一位置豎向位移與扭轉(zhuǎn)位移之間的相位差沿跨長的變化規(guī)律,可知FM-A跨中位置豎向位移與扭轉(zhuǎn)位移相位差最小,最小為14°,并向主梁展向兩端不斷增大. 對比FM-A、FM-B、FM-C和FM-D,非對稱型中央扣打亂了局部位置處豎向位移與扭轉(zhuǎn)位移相位差由跨中向兩邊增大的趨勢,增大了主梁左邊豎向位移與扭轉(zhuǎn)位移的相位差,最大達(dá)50°. 由此可見,中央扣對主梁的整個三維顫振姿態(tài)產(chǎn)生比較復(fù)雜的影響.

    3 ? 結(jié) ? 語

    本文通過模態(tài)分析和顫振時域分析研究了中央扣的不同結(jié)構(gòu)形式對大跨懸索橋自振頻率及顫振穩(wěn)定性的影響,得出以下結(jié)論:

    1)3種不同結(jié)構(gòu)形式的中央扣都能夠大幅提高主梁縱飄的頻率,其中剛性中央扣的增幅最大;對于側(cè)彎和扭轉(zhuǎn)振型,中央扣對反對稱振型的影響要遠(yuǎn)大于正對稱振型且剛性中央扣的影響最大,其中剛性中央扣能夠提高15%一階反對稱扭轉(zhuǎn)的頻率;此外,豎彎振型的頻率受中央扣的影響很小,主要由于中央扣僅提高了主梁跨中處的局部豎向剛度.

    2)在不考慮結(jié)構(gòu)阻尼情況下,3種不同結(jié)構(gòu)形式的中央扣對顫振臨界風(fēng)速的影響都極小. 此外,由于該橋的顫振頻率落在固有頻率分布密集的區(qū)域,因此,當(dāng)不考慮結(jié)構(gòu)阻尼或當(dāng)結(jié)構(gòu)阻尼比很低時顫振頻率很容易牽連帶動頻率密集區(qū)的固有模態(tài)間歇性參與振動,出現(xiàn)間歇性顫振現(xiàn)象.

    3)考慮結(jié)構(gòu)阻尼時,3種不同中央扣對顫振臨界風(fēng)速與顫振頻率的影響都極小,這是由于矮寨大橋以一階正對稱豎彎、二階正對稱豎彎和一階正對稱扭轉(zhuǎn)相互耦合的振型發(fā)生彎扭耦合顫振,而中央扣對這些振型頻率的影響甚微. 剛性中央扣能提高15%一階反對稱扭轉(zhuǎn)的頻率,若主梁發(fā)生反對稱豎彎和扭轉(zhuǎn)耦合的顫振,可以預(yù)見中央扣將能大幅提高顫振臨界風(fēng)速甚至改變顫振振型. 此外,通過主梁展向不同位置處顫振位移響應(yīng)之間相位差的研究,可知主梁不同位置之間振動的同步性隨著距離的增大而降低;該橋彎扭耦合振動時扭轉(zhuǎn)位移的展向同步性要遠(yuǎn)高于豎向位移;中央扣對主梁的整個三維顫振姿態(tài)產(chǎn)生比較復(fù)雜的影響,一定程度上有利于顫振穩(wěn)定性.

    主梁顫振時不同形式的中央扣其局部位置的應(yīng)力分布不同,當(dāng)結(jié)構(gòu)發(fā)生大幅振動時會對中央扣構(gòu)件產(chǎn)生很大的疲勞損傷,尤其是剛性中央扣,因此下一步可進(jìn)一步深入研究.

    參考文獻(xiàn)

    [1] ? ?劉志文,謝普仁,陳政清,等. 大跨度流線型箱梁懸索橋顫振穩(wěn)定性氣動優(yōu)化[J]. 湖南大學(xué)學(xué)報(自然科學(xué)版),2019,46(3):1—9.

    LIU Z W,XIE P R,CHEN Z Q,et al. Aerodynamic optimization of flutter stability for a long-span streamlined box girder suspension bridge[J]. Journal of Hunan University (Natural Sciences),2019,46(3):1—9.(In Chinese)

    [2] ? ?VIOLA J M,SYED S,CLENANCE J. The new Tacoma Narrows suspension bridge:Construction support and engineering[C]// Proceedings of the 2005 Structures Congress and the 2005 Forensic Engineering Symposium. New York:Structure Engineering Institute of the America Society of Civil Engineers,2005:1—12.

    [3] ? ?徐勛,強(qiáng)士中,賀栓海. 中央扣對大跨懸索橋動力特性和汽車車列激勵響應(yīng)的影響[J]. 中國公路學(xué)報,2008,21(6):57—63.

    XU X,QIANG S Z,HE S H. Influence of central buckle on dynamic behavior and response of long-span suspension bridge under vehicle group excitation[J]. ?China Journal of Highway and Transport,2008,21 (6):57—63. (In Chinese)

    [4] ? ?徐勛,強(qiáng)士中. 中央扣對大跨懸索橋動力特性和地震響應(yīng)的影響研究[J]. 鐵道學(xué)報,2010,32(4):84—91.

    XU X,QIANG S Z. Influence of central buckle on dynamic behavior and seismic response of long-span suspension bridge[J]. ?Journal of the China Railway Society,2010,32 (4):84—91. (In Chinese)

    [5] ? ?QIN F J,DI J,DAI J,et al. Influence of central buckle on dynamic behavior of long-span suspension bridge with deck-truss composite stiffening girder[J]. ?Advanced Materials Research,2013,838/841:1096— 1101.

    [6] ? ?胡騰飛,華旭剛,溫青,等. 中央扣對大跨懸索橋模態(tài)特性的影響[J]. 公路交通科技,2015,32(6):89—94.

    HU T F,HUA X G,WEN Q,et al. Influence of central buckles on modal characteristics of long-span suspension bridge[J]. ?Journal of Highway and Transportation Research and Development,2015,32 (6):89—94. (In Chinese)

    [7] ? ?彭旺虎,邵旭東. 設(shè)置中央扣懸索橋的扭轉(zhuǎn)自振分析[J]. 中國公路學(xué)報,2013,26(5):76—87.

    PENG W H,SHAO X D. Analysis of free torsional vibration of suspension bridges with center ties[J]. ?China Journal of Highway and Transport,2013,26 (5):76—87. (In Chinese)

    [8] ? ?WANG H,LI A Q,JIAO C K,et al. Influence of rigid central buckle on seismic response of long-span suspension bridges[J]. High Technology Letters,2011,17(2):214—219.

    [9] ? ?WANG D,CHEN C M,LIU Y. Central buckle influence research on dynamical characteristics of suspension bridge[J]. ?Applied Mechanics and Materials,2013,405/408:1489—1493.

    [10] WANG D,DENG Y,LIU Y. Influence of central buckle on suspension bridge dynamic characteristics and driving comfort[J]. ?Journal of Central South University,2015,22(8):3108—3115.

    [11] ?WANG H,LI A Q,ZHAO G W,et al. Non-linear buffeting response analysis of long-span suspension bridges with central buckle[J]. ?Earthquake Engineering and Engineering Vibration,2010,9(2):259—270.

    [12] ?WANG H,TAO T Y,ZHOU R,et al. Parameter sensitivity study on flutter stability of a long-span triple-tower suspension bridge[J]. ?Journal of Wind Engineering and Industrial Aerodynamics,2014,128:12—28.

    [13] ?李春光,陳政清,張志田. 大跨度桁式加勁梁懸索橋氣彈模型等效設(shè)計[J]. 振動與沖擊,2009,28(9):171—174.

    LI C G, CHEN Z Q,ZHANG Z T. Equivalent design of aeroelastic model for a long-span suspension bridge with truss stiffening girder [J]. ?Journal of Vibration and Shock,2009,28(9):171—174. (In Chinese)

    [14] ?吳長青,張志田,張偉峰. 考慮幾何非線性的橋梁后顫振極限環(huán)特性[J]. 湖南大學(xué)學(xué)報(自然科學(xué)版),2018,45(5):1—10.

    WU C Q,ZHANG Z T,ZHANG W F. Limit cycle oscillation of bridge post-flutter with geometric nonlinearities[J]. ?Journal of Hunan University(Natural Sciences),2018,45(5):1—10. (In Chinese)

    [15] ?LIN Y K,YANG J N. Multimode bridge response to wind excitations[J]. ?Journal of Engineering Mechanics,1983,109(2):586—603.

    [16] ?CHEN X,MATSUMOTO M,KAREEM A. Time domain flutter and buffeting response analysis of bridges [J]. ?Journal of Engineering Mechanics,2000,126(1):7—16.

    [17] ?張志田. 大跨度橋梁非線性抖振及其對抗風(fēng)穩(wěn)定性影響的研究[D]. ?上海:同濟(jì)大學(xué), 2004:69—71.

    ZHANG Z T. Study on the nonlinear buffeting of long-span bridges and its impact on wind stability[D]. Shanghai:Tongji University,2004:69—71. (In Chinese)

    [18] ?中國公路橋梁抗風(fēng)設(shè)計規(guī)范:JTG/T 3360-01—2018[S]. 北京:人民交通出版社,2018:35—36.

    Wind-resistant design specification for highway bridges:JTG/T 3360-01—2018[S]. Beijing:Peoples Communications Press,2018:35—36. (In Chinese)

    [19] ?陳政清. 矮寨大橋懸索橋風(fēng)洞試驗及風(fēng)振響應(yīng)分析研究總報告[R]. 長沙:湖南大學(xué)風(fēng)工程試驗中心,2013:53—55.

    CHEN Z Q. General report on wind tunnel test and wind vibration response analysis of the Aizhai Bridge[R]. Changsha:Wind Engineering Test Center of Hunan University,2013:53—55. (In Chinese)

    久久国产精品男人的天堂亚洲| 老司机午夜十八禁免费视频| 国产高清videossex| 美女福利国产在线| 精品少妇一区二区三区视频日本电影| 日本91视频免费播放| 亚洲精品美女久久av网站| 免费看十八禁软件| 一级片'在线观看视频| 91九色精品人成在线观看| 在线观看www视频免费| 久久久精品免费免费高清| 最黄视频免费看| 国产精品亚洲av一区麻豆| 日本色播在线视频| 亚洲美女黄色视频免费看| bbb黄色大片| 免费av中文字幕在线| 啦啦啦中文免费视频观看日本| 亚洲欧美日韩高清在线视频 | 欧美黑人欧美精品刺激| 精品久久久久久久毛片微露脸 | 亚洲第一av免费看| 精品一区二区三卡| 亚洲视频免费观看视频| 日韩电影二区| 成年人免费黄色播放视频| 嫁个100分男人电影在线观看 | 精品人妻熟女毛片av久久网站| 久久久亚洲精品成人影院| xxx大片免费视频| 亚洲成人手机| 黄色 视频免费看| 国产又色又爽无遮挡免| 男女免费视频国产| 日本vs欧美在线观看视频| av在线播放精品| 日韩中文字幕欧美一区二区 | 真人做人爱边吃奶动态| 老司机靠b影院| 深夜精品福利| 欧美日韩成人在线一区二区| 亚洲国产中文字幕在线视频| 热99国产精品久久久久久7| 99久久人妻综合| 黄色一级大片看看| 久9热在线精品视频| 黑人欧美特级aaaaaa片| 一级a爱视频在线免费观看| 国产真人三级小视频在线观看| 久久国产精品大桥未久av| 国产成人av教育| 国产99久久九九免费精品| 麻豆乱淫一区二区| 国产精品久久久久久精品电影小说| 一边摸一边做爽爽视频免费| 亚洲精品国产av蜜桃| 国产真人三级小视频在线观看| 少妇裸体淫交视频免费看高清 | 欧美精品av麻豆av| 丰满少妇做爰视频| 麻豆乱淫一区二区| 亚洲欧美一区二区三区国产| 国产精品99久久99久久久不卡| 咕卡用的链子| 国产精品免费视频内射| 悠悠久久av| 国产成人系列免费观看| 亚洲精品一二三| 国产男女超爽视频在线观看| 国产伦人伦偷精品视频| 制服诱惑二区| 成年女人毛片免费观看观看9 | 在线观看免费高清a一片| 免费一级毛片在线播放高清视频 | 999精品在线视频| 丝袜喷水一区| 亚洲国产精品999| 黑丝袜美女国产一区| 麻豆av在线久日| 一区二区av电影网| 久久久久精品国产欧美久久久 | 国产高清videossex| 色播在线永久视频| av电影中文网址| 又黄又粗又硬又大视频| 精品国产超薄肉色丝袜足j| 中文字幕亚洲精品专区| 亚洲成人免费av在线播放| 国产伦理片在线播放av一区| 美女脱内裤让男人舔精品视频| 久久久亚洲精品成人影院| 日本wwww免费看| 国产精品一区二区在线不卡| 99九九在线精品视频| 侵犯人妻中文字幕一二三四区| 亚洲伊人色综图| 蜜桃在线观看..| 男女之事视频高清在线观看 | 悠悠久久av| 成人免费观看视频高清| 久久这里只有精品19| 国产免费一区二区三区四区乱码| 97在线人人人人妻| 香蕉丝袜av| 可以免费在线观看a视频的电影网站| 欧美 亚洲 国产 日韩一| 国产日韩欧美亚洲二区| 国产一区有黄有色的免费视频| 岛国毛片在线播放| 亚洲欧美成人综合另类久久久| 黄色毛片三级朝国网站| 老汉色∧v一级毛片| 午夜免费成人在线视频| 久久国产精品男人的天堂亚洲| 成年人免费黄色播放视频| 99香蕉大伊视频| 黄片播放在线免费| 国产无遮挡羞羞视频在线观看| 亚洲精品美女久久久久99蜜臀 | 精品亚洲成国产av| 亚洲欧美精品自产自拍| 高清av免费在线| 亚洲五月色婷婷综合| 熟女av电影| 波多野结衣一区麻豆| 丁香六月欧美| 97在线人人人人妻| 国产成人精品无人区| 国产激情久久老熟女| 亚洲成人手机| 丝瓜视频免费看黄片| 国产精品二区激情视频| 又大又爽又粗| 18禁裸乳无遮挡动漫免费视频| 久久99热这里只频精品6学生| 中文字幕精品免费在线观看视频| 婷婷成人精品国产| 久久精品aⅴ一区二区三区四区| 久久中文字幕一级| 成人黄色视频免费在线看| 啦啦啦在线观看免费高清www| 真人做人爱边吃奶动态| 久久久久久免费高清国产稀缺| 一区福利在线观看| 在线看a的网站| 国产成人免费观看mmmm| 久久 成人 亚洲| 97在线人人人人妻| 亚洲av欧美aⅴ国产| 成人18禁高潮啪啪吃奶动态图| 国产欧美日韩一区二区三区在线| 成在线人永久免费视频| 色综合欧美亚洲国产小说| 午夜两性在线视频| 成人亚洲精品一区在线观看| 亚洲av成人精品一二三区| 黑人欧美特级aaaaaa片| 亚洲精品国产区一区二| 十八禁网站网址无遮挡| 十分钟在线观看高清视频www| 色婷婷久久久亚洲欧美| 国产成人系列免费观看| 纵有疾风起免费观看全集完整版| 国产成人91sexporn| 国精品久久久久久国模美| 97人妻天天添夜夜摸| 日本91视频免费播放| 国产成人精品在线电影| 免费在线观看黄色视频的| 久久久久久久久免费视频了| 侵犯人妻中文字幕一二三四区| avwww免费| 免费观看av网站的网址| 黄色片一级片一级黄色片| 两个人看的免费小视频| 777久久人妻少妇嫩草av网站| 国产野战对白在线观看| 一级毛片黄色毛片免费观看视频| 中文字幕最新亚洲高清| 国产深夜福利视频在线观看| 久久久久久久久久久久大奶| 久久99热这里只频精品6学生| 免费黄频网站在线观看国产| 日韩大码丰满熟妇| 91精品国产国语对白视频| 国产在线观看jvid| 久久九九热精品免费| 精品久久久精品久久久| 亚洲精品一二三| 别揉我奶头~嗯~啊~动态视频 | tube8黄色片| 十八禁人妻一区二区| 亚洲国产欧美网| 国产精品香港三级国产av潘金莲 | 亚洲欧美色中文字幕在线| 亚洲五月婷婷丁香| 国产在视频线精品| 91麻豆精品激情在线观看国产 | av电影中文网址| 1024香蕉在线观看| 99久久精品国产亚洲精品| 午夜激情av网站| 在线天堂中文资源库| 可以免费在线观看a视频的电影网站| 国产亚洲一区二区精品| 大陆偷拍与自拍| 久久午夜综合久久蜜桃| 韩国精品一区二区三区| 欧美精品人与动牲交sv欧美| 国产成人a∨麻豆精品| 精品一区二区三卡| 99热国产这里只有精品6| 老司机影院成人| 在线天堂中文资源库| 啦啦啦在线观看免费高清www| 丝袜在线中文字幕| 亚洲七黄色美女视频| 多毛熟女@视频| 日韩 亚洲 欧美在线| 免费看不卡的av| 国产精品 国内视频| 国产成人a∨麻豆精品| 80岁老熟妇乱子伦牲交| 亚洲欧美成人综合另类久久久| 九色亚洲精品在线播放| 99国产精品免费福利视频| 欧美黄色淫秽网站| 国产亚洲一区二区精品| 亚洲欧美一区二区三区黑人| 亚洲色图综合在线观看| 男女无遮挡免费网站观看| 国产成人影院久久av| 精品第一国产精品| 亚洲av综合色区一区| 亚洲久久久国产精品| 69精品国产乱码久久久| 精品国产国语对白av| 18禁国产床啪视频网站| 各种免费的搞黄视频| 热re99久久精品国产66热6| 精品国产乱码久久久久久男人| 久久国产精品人妻蜜桃| 一个人免费看片子| 观看av在线不卡| 国产片特级美女逼逼视频| 黄色一级大片看看| 中文字幕最新亚洲高清| 波多野结衣av一区二区av| 亚洲美女黄色视频免费看| 五月天丁香电影| 国产精品九九99| 午夜久久久在线观看| 欧美变态另类bdsm刘玥| 日韩制服骚丝袜av| 免费看av在线观看网站| 人成视频在线观看免费观看| 中文字幕高清在线视频| 久热这里只有精品99| 免费不卡黄色视频| 美女视频免费永久观看网站| 色精品久久人妻99蜜桃| 我要看黄色一级片免费的| 在线精品无人区一区二区三| 熟女少妇亚洲综合色aaa.| 91字幕亚洲| 欧美日韩成人在线一区二区| 老司机影院成人| 亚洲精品美女久久av网站| svipshipincom国产片| av天堂在线播放| 中文欧美无线码| 女性被躁到高潮视频| 国产在线视频一区二区| 欧美国产精品一级二级三级| 男女午夜视频在线观看| 别揉我奶头~嗯~啊~动态视频 | 亚洲成人国产一区在线观看 | 在线精品无人区一区二区三| √禁漫天堂资源中文www| 国产亚洲精品第一综合不卡| 99热国产这里只有精品6| 亚洲欧美激情在线| 最近最新中文字幕大全免费视频 | 国产成人精品久久二区二区91| 中文字幕高清在线视频| 91精品三级在线观看| 国产男人的电影天堂91| 极品少妇高潮喷水抽搐| 国产熟女午夜一区二区三区| 亚洲精品第二区| 在线看a的网站| 久久精品国产综合久久久| 亚洲精品美女久久av网站| 成人影院久久| 国产精品香港三级国产av潘金莲 | kizo精华| 国产激情久久老熟女| 午夜福利视频在线观看免费| 一区二区av电影网| 老司机靠b影院| 国产成人影院久久av| 美女国产高潮福利片在线看| 免费观看人在逋| 久久 成人 亚洲| 深夜精品福利| 国产精品一区二区在线不卡| 日韩欧美一区视频在线观看| 婷婷丁香在线五月| 免费黄频网站在线观看国产| 国产成人一区二区三区免费视频网站 | 99国产精品一区二区三区| 国产一区二区激情短视频 | 蜜桃国产av成人99| 18在线观看网站| 精品福利观看| 国产女主播在线喷水免费视频网站| 亚洲国产欧美网| www.av在线官网国产| 在线观看国产h片| 80岁老熟妇乱子伦牲交| 国产男女内射视频| 久久国产精品男人的天堂亚洲| 两人在一起打扑克的视频| 久久人人97超碰香蕉20202| 后天国语完整版免费观看| 侵犯人妻中文字幕一二三四区| 99国产精品一区二区蜜桃av | 久久性视频一级片| 国产亚洲av高清不卡| 国产女主播在线喷水免费视频网站| 亚洲精品国产av成人精品| 人人妻人人澡人人爽人人夜夜| 国产成人av激情在线播放| 蜜桃在线观看..| 丰满人妻熟妇乱又伦精品不卡| 国产精品香港三级国产av潘金莲 | 亚洲男人天堂网一区| svipshipincom国产片| avwww免费| 又粗又硬又长又爽又黄的视频| 少妇精品久久久久久久| 99久久精品国产亚洲精品| 日韩制服丝袜自拍偷拍| 99香蕉大伊视频| 国产精品久久久久久精品古装| 亚洲五月婷婷丁香| 亚洲中文字幕日韩| 亚洲一码二码三码区别大吗| 看免费av毛片| 搡老岳熟女国产| 十八禁人妻一区二区| 啦啦啦视频在线资源免费观看| 欧美在线黄色| 成人黄色视频免费在线看| 最近手机中文字幕大全| 久久精品国产亚洲av涩爱| 成人三级做爰电影| 纯流量卡能插随身wifi吗| 免费女性裸体啪啪无遮挡网站| 国产成人欧美在线观看 | 少妇 在线观看| 久久久久精品人妻al黑| 十分钟在线观看高清视频www| 国产黄色视频一区二区在线观看| 国产在线视频一区二区| 777米奇影视久久| 久久久久久人人人人人| 大话2 男鬼变身卡| 一级毛片女人18水好多 | xxx大片免费视频| 女人被躁到高潮嗷嗷叫费观| 国产成人精品久久二区二区91| 日日摸夜夜添夜夜爱| 国产精品香港三级国产av潘金莲 | 三上悠亚av全集在线观看| 日韩人妻精品一区2区三区| netflix在线观看网站| 免费看不卡的av| 国产片内射在线| 成年动漫av网址| 中文精品一卡2卡3卡4更新| 少妇精品久久久久久久| 亚洲国产欧美在线一区| 久久狼人影院| svipshipincom国产片| 国产午夜精品一二区理论片| 欧美另类一区| 国产午夜精品一二区理论片| 中国国产av一级| 国产91精品成人一区二区三区 | 国产亚洲午夜精品一区二区久久| 国产精品久久久人人做人人爽| 多毛熟女@视频| 国产午夜精品一二区理论片| 亚洲av电影在线观看一区二区三区| 国产又色又爽无遮挡免| 亚洲情色 制服丝袜| 欧美黄色片欧美黄色片| 久久国产精品人妻蜜桃| 日本黄色日本黄色录像| 精品高清国产在线一区| 叶爱在线成人免费视频播放| 夜夜骑夜夜射夜夜干| 国产免费福利视频在线观看| 咕卡用的链子| 捣出白浆h1v1| 男女国产视频网站| 欧美av亚洲av综合av国产av| 亚洲精品久久久久久婷婷小说| 9191精品国产免费久久| 国产av一区二区精品久久| 国产97色在线日韩免费| 久久久久久久大尺度免费视频| 一级片免费观看大全| 亚洲,一卡二卡三卡| 久久天堂一区二区三区四区| 亚洲av男天堂| 美女视频免费永久观看网站| 久久免费观看电影| 一本综合久久免费| 啦啦啦 在线观看视频| 成人18禁高潮啪啪吃奶动态图| 男女无遮挡免费网站观看| 日本午夜av视频| 天天躁日日躁夜夜躁夜夜| 亚洲人成77777在线视频| 99九九在线精品视频| 波野结衣二区三区在线| 午夜91福利影院| 国产成人精品久久久久久| 久久久久网色| 只有这里有精品99| 午夜精品国产一区二区电影| 性少妇av在线| 国产淫语在线视频| 久久热在线av| 视频区欧美日本亚洲| 午夜福利影视在线免费观看| 亚洲一码二码三码区别大吗| 欧美成人午夜精品| 别揉我奶头~嗯~啊~动态视频 | 成在线人永久免费视频| 国产片特级美女逼逼视频| 一级片免费观看大全| 韩国高清视频一区二区三区| bbb黄色大片| 丝袜喷水一区| 久久久久精品国产欧美久久久 | 欧美xxⅹ黑人| 国产精品一区二区精品视频观看| 只有这里有精品99| 亚洲 国产 在线| 国产成人欧美在线观看 | 美女福利国产在线| 日本a在线网址| 一级黄片播放器| 别揉我奶头~嗯~啊~动态视频 | 黄频高清免费视频| 人妻一区二区av| 久久久久精品国产欧美久久久 | 久久精品国产亚洲av涩爱| 美女福利国产在线| 国产福利在线免费观看视频| 在线观看www视频免费| 亚洲人成电影免费在线| 日本91视频免费播放| 夫妻午夜视频| 亚洲精品美女久久久久99蜜臀 | 日韩av在线免费看完整版不卡| 亚洲精品中文字幕在线视频| 极品人妻少妇av视频| 国产爽快片一区二区三区| 中文字幕人妻熟女乱码| 美女视频免费永久观看网站| 波多野结衣一区麻豆| 国产精品国产三级国产专区5o| 大型av网站在线播放| 欧美国产精品va在线观看不卡| 国产有黄有色有爽视频| 99久久99久久久精品蜜桃| 黄色视频不卡| 日韩免费高清中文字幕av| 色精品久久人妻99蜜桃| 国产爽快片一区二区三区| 久久久久国产一级毛片高清牌| kizo精华| 9191精品国产免费久久| 日韩一本色道免费dvd| 欧美 日韩 精品 国产| 在线av久久热| 欧美激情 高清一区二区三区| 亚洲中文日韩欧美视频| 大型av网站在线播放| 超碰97精品在线观看| 国产成人免费观看mmmm| 蜜桃在线观看..| 亚洲中文日韩欧美视频| 国产精品 国内视频| 少妇粗大呻吟视频| 色综合欧美亚洲国产小说| 亚洲,欧美,日韩| 精品国产国语对白av| 老鸭窝网址在线观看| 久久久久久久大尺度免费视频| 狂野欧美激情性xxxx| 欧美黄色片欧美黄色片| av电影中文网址| 国产成人影院久久av| 国产在线一区二区三区精| 夫妻午夜视频| 久久精品人人爽人人爽视色| 国产黄色视频一区二区在线观看| 又大又黄又爽视频免费| 久久亚洲精品不卡| 91精品国产国语对白视频| cao死你这个sao货| 在线亚洲精品国产二区图片欧美| 一级片'在线观看视频| 十八禁网站网址无遮挡| 免费日韩欧美在线观看| 久久九九热精品免费| 无遮挡黄片免费观看| 久久热在线av| 免费女性裸体啪啪无遮挡网站| 性色av一级| 首页视频小说图片口味搜索 | 亚洲国产欧美网| 男女边摸边吃奶| 国产深夜福利视频在线观看| 国产片内射在线| 免费高清在线观看日韩| 99热全是精品| 一级毛片 在线播放| 熟女av电影| 美国免费a级毛片| 亚洲五月婷婷丁香| 亚洲国产毛片av蜜桃av| 久久久久久久久免费视频了| 性色av乱码一区二区三区2| 国产精品一二三区在线看| 亚洲免费av在线视频| 免费在线观看完整版高清| 亚洲欧美清纯卡通| 中文欧美无线码| 国产精品久久久av美女十八| 免费在线观看日本一区| 人成视频在线观看免费观看| 伊人久久大香线蕉亚洲五| 一级毛片电影观看| 色94色欧美一区二区| 又紧又爽又黄一区二区| 亚洲国产av新网站| 18禁国产床啪视频网站| 婷婷色综合www| 欧美日韩视频精品一区| 亚洲第一av免费看| 国产av国产精品国产| 母亲3免费完整高清在线观看| 日韩制服丝袜自拍偷拍| av天堂久久9| 中文字幕另类日韩欧美亚洲嫩草| 天堂俺去俺来也www色官网| 爱豆传媒免费全集在线观看| 国产一区二区三区综合在线观看| 不卡av一区二区三区| 99精品久久久久人妻精品| 美女中出高潮动态图| 日日爽夜夜爽网站| 日韩av不卡免费在线播放| 精品人妻1区二区| 亚洲情色 制服丝袜| 男女免费视频国产| 成年av动漫网址| 中文字幕最新亚洲高清| 黄频高清免费视频| 中文字幕亚洲精品专区| 精品一区在线观看国产| 亚洲九九香蕉| 老司机靠b影院| √禁漫天堂资源中文www| 大陆偷拍与自拍| 亚洲av综合色区一区| 美女扒开内裤让男人捅视频| 夫妻性生交免费视频一级片| 久久久久久免费高清国产稀缺| 一级片'在线观看视频| 国产免费现黄频在线看| 日韩中文字幕欧美一区二区 | 亚洲第一青青草原| 99国产精品99久久久久| 精品国产乱码久久久久久男人| 亚洲精品第二区| 国产国语露脸激情在线看| 日本91视频免费播放| 丰满人妻熟妇乱又伦精品不卡| 国产欧美日韩精品亚洲av| 国产精品99久久99久久久不卡| e午夜精品久久久久久久| 十八禁网站网址无遮挡| 狂野欧美激情性bbbbbb| 亚洲成国产人片在线观看| 国产在线免费精品| 日韩一区二区三区影片| 中文精品一卡2卡3卡4更新| 看十八女毛片水多多多| 丰满饥渴人妻一区二区三| 久久久久国产一级毛片高清牌| 亚洲av日韩精品久久久久久密 | 国产女主播在线喷水免费视频网站| 亚洲av成人不卡在线观看播放网 | 婷婷色av中文字幕| 欧美中文综合在线视频| 你懂的网址亚洲精品在线观看| 三上悠亚av全集在线观看|