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

    基于實(shí)測(cè)與規(guī)范風(fēng)譜的三塔懸索橋抖振性能對(duì)比

    2013-12-29 06:40:52陶天友李?lèi)?ài)群鄧穩(wěn)平
    關(guān)鍵詞:中塔三塔懸索橋

    王 浩 陶天友 郭 彤 李?lèi)?ài)群 鄧穩(wěn)平

    (東南大學(xué)混凝土及預(yù)應(yīng)力混凝土結(jié)構(gòu)教育部重點(diǎn)實(shí)驗(yàn)室, 南京 210096)

    隨著跨徑的增加,懸索橋結(jié)構(gòu)變得極為細(xì)柔,在風(fēng)荷載作用下的響應(yīng)逐漸占據(jù)主導(dǎo)作用[1-4].泰州長(zhǎng)江公路大橋(以下簡(jiǎn)稱泰州大橋)為主跨2×1080m的三塔雙跨鋼箱梁懸索橋,系世界首創(chuàng).與兩塔懸索橋相比,中塔的引入使得三塔懸索橋的受力特點(diǎn)與動(dòng)力特性變得更加復(fù)雜,在風(fēng)荷載激勵(lì)下的響應(yīng)應(yīng)該給予密切關(guān)注.抖振作為橋梁風(fēng)致振動(dòng)最突出的形式之一,已成為橋梁抗風(fēng)研究中的重點(diǎn)課題之一,且已開(kāi)展了大量相關(guān)研究[5-7],然而,國(guó)內(nèi)外對(duì)于三塔懸索橋抗風(fēng)性能的研究很少.鑒于大跨度三塔懸索橋新穎的結(jié)構(gòu)型式及其復(fù)雜的受力特點(diǎn),開(kāi)展其抖振研究具有重要的理論意義與工程應(yīng)用價(jià)值.

    我國(guó)地處亞洲大陸東部,地域廣大,不同的地形特點(diǎn)會(huì)引起所在地風(fēng)場(chǎng)特性的差異,因而單一的抗風(fēng)規(guī)范很難完整地詮釋全國(guó)各地的風(fēng)譜模型[8].本文以世界第一大跨度三塔懸索橋——泰州大橋?yàn)楣こ瘫尘?采用實(shí)測(cè)風(fēng)譜對(duì)該橋設(shè)計(jì)階段基于規(guī)范譜所進(jìn)行的抖振分析進(jìn)行對(duì)比驗(yàn)證.實(shí)測(cè)風(fēng)譜源自本課題組多年來(lái)在潤(rùn)揚(yáng)長(zhǎng)江大橋橋址區(qū)、蘇通長(zhǎng)江大橋橋址區(qū)的多次臺(tái)風(fēng)現(xiàn)場(chǎng)實(shí)測(cè)以及結(jié)構(gòu)健康監(jiān)測(cè)系統(tǒng)[9]數(shù)據(jù)庫(kù)所獲得的海鷗、鳳凰、卡努、麥莎、冬季強(qiáng)北風(fēng)等數(shù)據(jù)資料.根據(jù)以上實(shí)測(cè)資料,采用非線性最小二乘法對(duì)橋址區(qū)實(shí)測(cè)強(qiáng)風(fēng)紊流功率譜密度函數(shù)進(jìn)行擬合,獲得實(shí)測(cè)風(fēng)譜曲線,并分別以實(shí)測(cè)風(fēng)譜和規(guī)范譜為目標(biāo)譜模擬橋址區(qū)的三維脈動(dòng)風(fēng)場(chǎng),基于ANSYS對(duì)泰州大橋進(jìn)行非線性時(shí)域抖振響應(yīng)的對(duì)比分析,結(jié)果可為三塔懸索橋抖振分析及抗風(fēng)設(shè)計(jì)提供參考.

    1 工程背景

    泰州大橋位于江蘇省境內(nèi)長(zhǎng)江中段,橋位處江面寬約2.4km,河床呈“W”形,綜合考慮橋位處河勢(shì)水文與航道狀況以及深水岸線的充分利用,設(shè)計(jì)時(shí)主橋橋跨結(jié)構(gòu)布置為390m+1080m+1080m+390m的兩主跨三塔懸索橋(見(jiàn)圖1).加勁梁采用封閉式流線型扁平鋼箱梁,加勁梁設(shè)上斜腹板及下斜腹板構(gòu)成導(dǎo)風(fēng)嘴.鋼箱梁節(jié)段標(biāo)準(zhǔn)長(zhǎng)度為16m,中心線處梁高3.5m.主纜在設(shè)計(jì)成橋狀態(tài)矢跨比為1/9,2根主纜橫向間距為35.8m.索塔采用門(mén)式框架結(jié)構(gòu),兩邊塔為混凝土塔,中間塔為鋼塔,邊塔高178.0m,中塔高182.5m.2個(gè)邊塔在順橋向?yàn)閱沃谓Y(jié)構(gòu),而中間橋塔在順橋向則采用倒Y形結(jié)構(gòu),以增強(qiáng)縱向剛度.

    圖1 泰州大橋結(jié)構(gòu)布置圖(單位:m)

    該三塔懸索橋邊跨無(wú)懸吊加勁梁,中塔下橫梁上不設(shè)豎向支座,也不設(shè)0號(hào)索,在中塔橫梁與主梁間設(shè)橫向抗風(fēng)支座,限制主梁的橫向位移,豎向采用限位擋塊.中塔處設(shè)彈性索適當(dāng)限制縱向位移;下塔柱底部與沉井上蓋板相連;邊塔處設(shè)置豎向拉壓支座、側(cè)向抗風(fēng)支座和縱向阻尼支座.

    2 結(jié)構(gòu)有限元模型

    基于大型有限元分析軟件ANSYS平臺(tái)建立了泰州大橋的有限元計(jì)算模型,如圖2所示.在該有限元模型中,加勁梁采用魚(yú)骨梁式的計(jì)算模型,加勁梁和橋塔等構(gòu)件簡(jiǎn)化為空間梁?jiǎn)卧?加勁梁按吊桿的吊點(diǎn)進(jìn)行離散.主纜、吊桿則簡(jiǎn)化為空間桿單元,主纜按吊桿的吊點(diǎn)進(jìn)行離散,其彈性模量采用Ernst等效彈模公式進(jìn)行計(jì)算[10].主梁與吊桿之間采用剛性橫梁模擬,橫梁采用剛臂(剛度極大的無(wú)質(zhì)量梁?jiǎn)卧?模擬.分析中將主梁的質(zhì)量平均分配到各節(jié)點(diǎn)上,通過(guò)節(jié)點(diǎn)上的集中質(zhì)量慣矩來(lái)考慮主梁分布質(zhì)量慣矩的作用.根據(jù)結(jié)構(gòu)設(shè)計(jì)圖紙,耦合了主梁與邊塔在豎向、橫橋向以及繞順橋向的轉(zhuǎn)動(dòng)自由度,同時(shí)耦合了主梁與中塔在橫橋向的自由度.在邊、中塔下橫梁與加勁梁交界處設(shè)置Combin14非線性阻尼器單元,以模擬用于控制縱橋向位移的彈性拉索和縱向阻尼支座,并將同一位置設(shè)置同向阻尼器,通過(guò)改變阻尼系數(shù)來(lái)等效多個(gè)阻尼器在該方向上的輸出力.主纜和主塔底部完全固結(jié),未考慮土-樁-結(jié)構(gòu)相互作用.

    圖2 泰州大橋空間有限元計(jì)算模型

    3 泰州大橋三維脈動(dòng)風(fēng)場(chǎng)的模擬

    分別以《公路橋梁抗風(fēng)設(shè)計(jì)規(guī)范》[11]中的規(guī)范譜和實(shí)測(cè)風(fēng)譜為目標(biāo)譜,結(jié)合大跨度三塔懸索橋的結(jié)構(gòu)特點(diǎn)對(duì)基于譜特性的三維脈動(dòng)風(fēng)場(chǎng)進(jìn)行了簡(jiǎn)化模擬,并從功率譜和相關(guān)性2個(gè)方面驗(yàn)證了該模擬方法,旨在為該橋抖振響應(yīng)非線性時(shí)域分析提供準(zhǔn)確可靠的脈動(dòng)風(fēng)速時(shí)程.

    泰州大橋風(fēng)場(chǎng)可以簡(jiǎn)化成5個(gè)獨(dú)立的一維多變量隨機(jī)風(fēng)場(chǎng),如表1所示.運(yùn)用風(fēng)場(chǎng)模擬的諧波合成法,對(duì)沿跨度方向80m(鋼箱梁標(biāo)準(zhǔn)節(jié)段長(zhǎng)度的5倍)等間距分布的主梁風(fēng)場(chǎng)和沿豎向30m等間距分布的主塔風(fēng)場(chǎng)進(jìn)行了模擬.模擬點(diǎn)分布如圖3所示.

    表1 簡(jiǎn)化的5個(gè)獨(dú)立一維脈動(dòng)風(fēng)場(chǎng)

    3.1 實(shí)測(cè)風(fēng)譜模型

    2005年—2012年間,本課題組在潤(rùn)揚(yáng)大橋、蘇通大橋進(jìn)行多次臺(tái)風(fēng)實(shí)測(cè)并基于橋梁結(jié)構(gòu)健康監(jiān)測(cè)系統(tǒng)數(shù)據(jù)庫(kù)得到了海鷗、鳳凰、韋帕、卡努、麥莎、冬季強(qiáng)北風(fēng)等數(shù)據(jù)資料.基于Matlab平臺(tái),經(jīng)過(guò)數(shù)據(jù)預(yù)處理進(jìn)行了橋址區(qū)風(fēng)特性分析,運(yùn)用FFT技術(shù)得到實(shí)測(cè)風(fēng)譜的PSD(power spectrum density)函數(shù).運(yùn)用非線性最小二乘法對(duì)縱向及豎向的經(jīng)驗(yàn)譜進(jìn)行擬合以得出實(shí)測(cè)風(fēng)譜(見(jiàn)圖4).

    圖4 實(shí)測(cè)風(fēng)譜與規(guī)范風(fēng)譜對(duì)比圖

    高度Z處平均風(fēng)速U(z)的順風(fēng)向脈動(dòng)風(fēng)功率譜密度函數(shù)為

    (1)

    豎向脈動(dòng)風(fēng)功率譜密度函數(shù)為

    (2)

    式中,Su(n)為順風(fēng)向功率譜密度函數(shù);Sw(n)為豎向功率譜密度函數(shù);n為自然風(fēng)的脈動(dòng)頻率;f為莫寧坐標(biāo);u為氣流摩阻速度.

    3.2 主梁風(fēng)場(chǎng)模擬

    運(yùn)用諧波合成法[12]獲得該橋主梁28個(gè)點(diǎn)的橫橋向與豎向脈動(dòng)風(fēng)速時(shí)程,主梁橫橋向目標(biāo)譜分別采用Kaimal譜[13]和實(shí)測(cè)譜的順風(fēng)向譜,豎向采用Panofsky譜[14]和實(shí)測(cè)譜的豎向譜,脈動(dòng)風(fēng)速譜的相干函數(shù)采用Davenport相干函數(shù)[1].選用其他參數(shù)具體如下:截止頻率ωu=4π rad/s;頻率分段數(shù)N=1024;樣本時(shí)距Δt=0.25s;選用時(shí)段長(zhǎng)1800s.根據(jù)現(xiàn)行《公路橋梁抗風(fēng)設(shè)計(jì)規(guī)范》,取橋址區(qū)百年一遇設(shè)計(jì)風(fēng)速為模擬風(fēng)速,即6.6m高度處平均風(fēng)速為27.1m/s,根據(jù)風(fēng)速沿高度的指數(shù)變化規(guī)律對(duì)主梁和主塔高度處風(fēng)速進(jìn)行換算.共計(jì)生成了50個(gè)風(fēng)場(chǎng)樣本,并從功率譜和相關(guān)性2個(gè)方面對(duì)模擬結(jié)果進(jìn)行了校核.圖5給出了主梁模擬自譜與目標(biāo)譜的對(duì)比圖,圖6給出了主梁模擬風(fēng)場(chǎng)相關(guān)函數(shù)與理論值的對(duì)比圖.

    由上述功率譜和相關(guān)函數(shù)對(duì)比圖可知,模擬譜與目標(biāo)譜對(duì)應(yīng)值十分接近,模擬的相關(guān)函數(shù)與理論值也相差不大,驗(yàn)證了本文所采用的三維脈動(dòng)風(fēng)場(chǎng)模擬方法的準(zhǔn)確性與可靠性.

    3.3 主塔風(fēng)場(chǎng)模擬

    由于橋塔上各點(diǎn)的風(fēng)速不同,模擬橋塔的風(fēng)速時(shí)程時(shí),互譜密度函數(shù)需要考慮相位角.風(fēng)譜密度矩陣(復(fù)數(shù)矩陣)有可能不正定,根據(jù)對(duì)應(yīng)風(fēng)譜密度矩陣題正定性將風(fēng)頻值劃分為正定和不正定區(qū)間.在風(fēng)頻值正定區(qū)間內(nèi),可直接將風(fēng)譜密度矩陣進(jìn)行Cholesky分解;而在風(fēng)頻值不正定區(qū)間內(nèi),采用插值技術(shù)近似獲得風(fēng)譜密度矩陣分解式[15].

    主塔順風(fēng)向目標(biāo)譜分別采用Kaimal譜[13]和實(shí)測(cè)譜的順風(fēng)向譜,脈動(dòng)風(fēng)速譜采用Davenport相干函數(shù)[1].選用其他參數(shù)如下:截止頻率ωu=4π rad/s;頻率分段數(shù)N=1024;樣本時(shí)距Δt=0.25s;選用時(shí)段長(zhǎng)1800s.同樣利用上述參數(shù)生成了50個(gè)風(fēng)場(chǎng)樣本,并對(duì)模擬結(jié)果進(jìn)行了校核.圖7和圖8分別給出了主塔模擬脈動(dòng)風(fēng)速在功率譜與相關(guān)性2個(gè)方面的對(duì)比,對(duì)比分析說(shuō)明諧波合成法的正確性,可為后續(xù)大橋風(fēng)致抖振時(shí)域分析奠定基礎(chǔ).

    圖5 主梁模擬自譜與目標(biāo)譜對(duì)比

    4 泰州大橋抖振響應(yīng)分析

    4.1 基于ANSYS的三塔懸索橋抖振時(shí)域分析方法

    在上述風(fēng)場(chǎng)模擬的基礎(chǔ)上,基于ANSYS的瞬態(tài)動(dòng)力學(xué)分析功能,進(jìn)行了該橋抖振響應(yīng)的非線性時(shí)域分析,其中計(jì)入了氣動(dòng)自激力的影響,具體包括如下步驟:① 基于參數(shù)化設(shè)計(jì)語(yǔ)言APDL建立了大橋有限元計(jì)算模型;② 為了便于APDL進(jìn)行調(diào)用,將主梁斷面的靜力三分力系數(shù)、顫振導(dǎo)數(shù)等以Table方式存儲(chǔ);③ 根據(jù)模擬風(fēng)場(chǎng)中的風(fēng)速數(shù)據(jù)和顫振導(dǎo)數(shù)值確定用于模擬氣動(dòng)剛度和氣動(dòng)阻尼的Matrix27單元的參數(shù)[16],得到了便于抖振分析、考慮氣動(dòng)自激力影響的有限元模型;④ 采用風(fēng)洞實(shí)測(cè)靜力三分力系數(shù)計(jì)算靜風(fēng)荷載;⑤ 根據(jù)Devaport經(jīng)典公式計(jì)算抖振力時(shí)程;⑥ 將風(fēng)荷載施加到結(jié)構(gòu)有限元模型上,并進(jìn)行瞬態(tài)動(dòng)力學(xué)求解,分析過(guò)程中考慮幾何非線性的影響;⑦ 進(jìn)入后處理查看結(jié)構(gòu)的時(shí)程響應(yīng)分析結(jié)果,求出統(tǒng)計(jì)量(如位移響應(yīng)RMS值等).

    圖6 主梁相關(guān)函數(shù)對(duì)比(點(diǎn)1~點(diǎn)28)

    4.2 泰州大橋抖振位移響應(yīng)RMS分布

    采用上述抖振時(shí)域分析方法進(jìn)行分析,計(jì)算時(shí)長(zhǎng)取30min,時(shí)間步長(zhǎng)取0.25s,共計(jì)7200個(gè)荷載子步.根據(jù)輸入的規(guī)范譜及實(shí)測(cè)譜模擬風(fēng)場(chǎng),得到泰州大橋沿跨度方向的抖振位移響應(yīng)RMS值分布,如圖9所示.

    圖7 中塔模擬自譜與目標(biāo)譜對(duì)比

    圖8 中塔相關(guān)函數(shù)對(duì)比(點(diǎn)29~點(diǎn)35)

    圖9表明:① 泰州大橋各主跨主梁沿跨度方向側(cè)向抖振位移RMS值(LRMS)分布與兩塔懸索橋類(lèi)似,均表現(xiàn)為由跨中向兩側(cè)遞減,在邊塔或中塔處為0,此現(xiàn)象歸結(jié)于主梁在中塔與邊塔處橫向設(shè)抗風(fēng)支座,限制主梁側(cè)向位移.② 主梁沿跨度方向豎向與扭轉(zhuǎn)抖振位移RMS值(VRMS和TRMS)的分布表現(xiàn)出泰州大橋獨(dú)有特點(diǎn),主梁扭轉(zhuǎn)抖振位移RMS值分布表現(xiàn)為由跨中向兩側(cè)遞減,在邊塔處為0,中塔處大于0;主梁豎向抖振位移RMS值分布呈現(xiàn)多極值,其中邊塔處RMS值為0,RMS值最大出現(xiàn)在跨中附近,極小值出現(xiàn)在距中塔約200m處,中塔處為另一極值點(diǎn).主梁于邊塔處豎向設(shè)拉壓支座,于中塔下橫梁處設(shè)豎向限位擋塊且未設(shè)0號(hào)索,這是形成三塔懸索橋豎向與扭轉(zhuǎn)抖振位移RMS值變異性的主要原因.③ 基于實(shí)測(cè)譜的抖振位移RMS值均小于規(guī)范譜,表明基于規(guī)范譜的抖振計(jì)算結(jié)果偏于保守,采用實(shí)測(cè)譜使得泰州大橋抖振計(jì)算更加精細(xì)化.

    圖9 主梁沿跨度方向抖振位移響應(yīng)RMS值分布圖

    4.3 泰州大橋主梁抖振位移響應(yīng)功率譜對(duì)比分析

    為了分析大橋各階振型對(duì)抖振響應(yīng)的貢獻(xiàn),進(jìn)行了抖振位移響應(yīng)的功率譜分析.由于主梁所受約束由跨中向邊塔逐漸減小,跨中的風(fēng)致抖振響應(yīng)往往最大,因此選擇跨中位置為代表進(jìn)行泰州大橋抖振位移響應(yīng)功率譜分析.基于FFT技術(shù)將跨中抖振響應(yīng)時(shí)程由時(shí)域轉(zhuǎn)換到頻域,進(jìn)行不同風(fēng)譜下主梁跨中抖振位移響應(yīng)功率譜比較(見(jiàn)圖10).為了統(tǒng)一單位,圖中扭轉(zhuǎn)位移已乘以主梁寬度的1/2.

    不同風(fēng)譜下主梁跨中位移功率譜密度表明:① 采用實(shí)測(cè)譜與規(guī)范譜的主梁跨中抖振位移功率譜密度曲線的整體趨勢(shì)及峰值出現(xiàn)位置基本一致,其中低頻區(qū)表現(xiàn)不一,高頻區(qū)基本相同.② 泰州大橋豎向抖振位移響應(yīng)功率譜密度曲線的第一個(gè)峰值對(duì)應(yīng)結(jié)構(gòu)自振特性的第二階振型,為主梁第一階反對(duì)稱豎彎振型;側(cè)向位移響應(yīng)功率譜密度曲線的第一個(gè)峰值對(duì)應(yīng)結(jié)構(gòu)自振特性的第一階振型,為主梁第一階反對(duì)稱側(cè)彎振型;扭轉(zhuǎn)位移響應(yīng)功率譜密度曲線的第一個(gè)峰值對(duì)應(yīng)結(jié)構(gòu)自振特性的第十三階振型,為主梁第一階反對(duì)稱扭轉(zhuǎn)振型.可見(jiàn),各方面的第一階振型對(duì)泰州大橋主梁抖振響應(yīng)貢獻(xiàn)最大.

    5 結(jié)論

    1) 采用實(shí)測(cè)譜與規(guī)范譜的主梁跨中抖振位移功率譜密度曲線的整體趨勢(shì)及峰值出現(xiàn)位置基本一致,其中低頻區(qū)表現(xiàn)不一,高頻區(qū)基本相同.

    2) 泰州大橋主梁的豎向抖振響應(yīng)主要由結(jié)構(gòu)第一階反對(duì)稱豎彎振型控制,側(cè)向抖振響應(yīng)主要由結(jié)構(gòu)第一階反對(duì)稱側(cè)彎振型控制,扭轉(zhuǎn)抖振響應(yīng)主要由結(jié)構(gòu)第一階反對(duì)稱扭轉(zhuǎn)振型控制.各方向的第一階振型對(duì)泰州大橋主梁抖振響應(yīng)貢獻(xiàn)最大.

    3) 泰州大橋各主跨主梁沿跨度方向側(cè)向抖振位移RMS值分布與兩塔懸索橋類(lèi)似,均表現(xiàn)為由跨中向兩側(cè)遞減,在邊塔或中塔處為0,這主要是由于泰州大橋主梁在邊塔和中塔處設(shè)橫向抗風(fēng)支座,限制主梁側(cè)向位移.

    4) 在靠近中塔附近,主梁沿跨度方向豎向與扭轉(zhuǎn)抖振位移RMS值分布表現(xiàn)出泰州大橋的獨(dú)有特點(diǎn):主梁扭轉(zhuǎn)抖振位移RMS值在中塔處大于0,主梁豎向抖振位移RMS值分布呈現(xiàn)多極值.這主要是由于主梁在邊塔位置設(shè)豎向拉壓支座,在中塔下橫梁處設(shè)豎向限位擋塊且未設(shè)0號(hào)索.

    5) 就本次研究而言,采用規(guī)范譜的抖振位移響應(yīng)RMS值與采用實(shí)測(cè)譜相比偏于保守,因此從數(shù)值計(jì)算的結(jié)果來(lái)看:采用基于實(shí)測(cè)的實(shí)測(cè)譜更有利于進(jìn)行大跨度橋梁抖振精細(xì)化研究,而規(guī)范譜的計(jì)算偏于安全,可用于工程設(shè)計(jì).

    必須指出,本文的實(shí)測(cè)風(fēng)譜模型僅由蘇通和潤(rùn)揚(yáng)2座大橋處的實(shí)測(cè)風(fēng)速資料得到,且由于兩橋結(jié)構(gòu)健康監(jiān)測(cè)系統(tǒng)運(yùn)營(yíng)時(shí)間尚短,使得實(shí)測(cè)風(fēng)速資料有限,這也是文章的不足之處.隨著未來(lái)南京長(zhǎng)江四橋、泰州大橋等大跨度橋梁的相繼建成,以及監(jiān)測(cè)數(shù)據(jù)隨時(shí)間的不斷積累,本文所提出的實(shí)測(cè)風(fēng)譜模型也將進(jìn)一步完善,這也是作者未來(lái)的研究方向.

    )

    [1] Simiu E, Scanlan R H.Windeffectsonstructures[M]. New York: John Wiley & Sons, INC, 1996.

    [2] Scanlan R H. The action of flexible bridge under wind. Part 2: buffeting theory[J].JournalofSoundandVibration, 1978,160(2): 201-211.

    [3] 項(xiàng)海帆, 葛耀君, 朱樂(lè)東, 等. 現(xiàn)代橋梁抗風(fēng)理論與實(shí)踐[M]. 北京: 人民交通出版社, 2005.

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

    [5] 李永樂(lè), 廖海黎, 強(qiáng)士中. 三塔斜拉橋抖振的耦合行為研究[J]. 工程力學(xué), 2004, 21(4): 184-188.

    Li Yongle, Liao Haili, Qiang Shizhong. Coupled behavior in buffeting response of three-pylon cable-stayed bridge[J].EngineeringMechanics, 2004,21(4): 184-188.(in Chinese)

    [6] 閔志華, 孫利民, 淡丹輝. 臺(tái)風(fēng)下斜拉橋風(fēng)致振動(dòng)和動(dòng)力特性分析[J]. 同濟(jì)大學(xué)學(xué)報(bào):自然科學(xué)版, 2009, 37(9): 1139-1145.

    Min Zhihua, Sun Limin, Dan Danhui. Analysis of wind-induced response and dynamic properties of cable-stayed bridge under typhoon[J].JournalofTongjiUniversity:NaturalScience, 2009,37(9): 1139-1145.(in Chinese)

    [7] 郭薇薇, 夏禾, 徐幼麟. 風(fēng)荷載作用下大跨度懸索橋的動(dòng)力響應(yīng)及列車(chē)運(yùn)行安全分析[J]. 工程力學(xué), 2006, 23(2): 103-110.

    Guo Weiwei, Xia He, Xu Youlin. Dynamic response of long span suspension bridge and running safety of train under wind action[J].EngineeringMechanics, 2006,23(2): 103-110.(in Chinese)

    [8] 王浩, 李?lèi)?ài)群, 焦長(zhǎng)科, 等. 基于規(guī)范及實(shí)測(cè)風(fēng)譜的蘇通大橋抖振響應(yīng)對(duì)比研究[J]. 土木工程學(xué)報(bào), 2011, 44(10): 91-97.

    Wang Hao, Li Aiqun, Jiao Changke, et al. Comparable study of the buffeting response of Sutong Bridge based on specification and measured wind power spectrum[J].ChinaCivilEngineeringJournal, 2011,44(10): 91-97.(in Chinese)

    [9] 李惠, 周文松, 歐進(jìn)萍,等. 大型橋梁結(jié)構(gòu)智能健康監(jiān)測(cè)系統(tǒng)集成技術(shù)研究[J]. 土木工程學(xué)報(bào), 2006, 39(2): 46-52.

    Li Hui, Zhou Wensong, Ou Jinping, et al. A study on system integration technique of intelligent monitoring systems for soundness of long-span bridges[J].ChinaCivilEngineeringJournal, 2006,39(2): 46-52.(in Chinese)

    [10] Ernst H J. Der e-modul von seilenunterber?cksichtigung des durchhanges[J].DerBauingenieur, 1965,40(2): 52-55. (in German)

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

    [12] Shinozuka M, Jan C M. Digital simulation of random processes and its application[J].JournalofSoundandVibration, 1972,25(1): 111-128.

    [13] Kaimal J C. Spectral characteristics of surface layer turbulence[J].JournaloftheRoyalMetaorologicalSociety, 1972,98(1): 563-589.

    [14] Panofsky H A, McCormick R A. The spectrum of vertical velocity near the surface[J].JournaloftheRoyalMetaorologicalSociety, 1960,86(370): 546-564.

    [15] 賀藝華, 蔣友寶, 張建仁. 考慮相位角的脈動(dòng)風(fēng)場(chǎng)模擬[J]. 工程力學(xué), 2008, 25(12): 25-29.

    He Yihua, Jiang Youbao, Zhang Jianren. Simulation of turbulent wind field with phase angle[J].EngineeringMechanics, 2008,25(12): 25-29.(in Chinese)

    [16] 華旭剛, 陳政清, 祝志文. 在ANSYS中實(shí)現(xiàn)顫振時(shí)程分析的方法[J]. 中國(guó)公路學(xué)報(bào), 2002, 15(4): 32-34.

    Hua Xugang, Chen Zhengqing, Zhu Zhiwen. Approach of time-history analysis of flutter in ANSYS[J].ChinaJournalofHighwayandTransport, 2002,15(4): 32-34.(in Chinese)

    猜你喜歡
    中塔三塔懸索橋
    一類(lèi)模擬懸索橋的von Kármán方程的解的存在性
    懸索橋主纜除濕控制系統(tǒng)
    中外漢語(yǔ)教師教學(xué)情況的比較與思考
    不同參數(shù)下三塔懸索橋動(dòng)力特性敏感性分析
    卡拉蘇口岸:中塔貿(mào)易風(fēng)生水起
    站在大跨徑懸索橋施工技術(shù)創(chuàng)新的最前沿
    中塔企業(yè)成功簽約
    車(chē)流荷載作用下三塔懸索橋主梁動(dòng)力響應(yīng)分析
    矮寨特大懸索橋鋼桁架安裝
    三線合一、三塔懸索橋兩設(shè)計(jì)方案靜風(fēng)穩(wěn)定性對(duì)比研究
    国产探花极品一区二区| 不卡视频在线观看欧美| 国产国语露脸激情在线看| 丝袜喷水一区| 亚洲欧美激情在线| 婷婷色综合www| 中文字幕人妻熟女乱码| 肉色欧美久久久久久久蜜桃| 午夜福利一区二区在线看| 99久久人妻综合| 久久人人爽人人片av| 亚洲国产欧美网| 久久午夜综合久久蜜桃| 亚洲av成人精品一二三区| 色精品久久人妻99蜜桃| 国产av精品麻豆| 久久精品国产亚洲av涩爱| 精品人妻在线不人妻| 日日啪夜夜爽| 另类精品久久| 老司机在亚洲福利影院| 色播在线永久视频| 欧美少妇被猛烈插入视频| 日日撸夜夜添| 免费观看a级毛片全部| 国产女主播在线喷水免费视频网站| 亚洲成色77777| 高清视频免费观看一区二区| 成人午夜精彩视频在线观看| 99re6热这里在线精品视频| 亚洲精品国产区一区二| 777久久人妻少妇嫩草av网站| 午夜91福利影院| 青春草亚洲视频在线观看| 亚洲av福利一区| 极品人妻少妇av视频| 看非洲黑人一级黄片| 免费人妻精品一区二区三区视频| 亚洲一码二码三码区别大吗| 精品国产露脸久久av麻豆| 女性生殖器流出的白浆| 一本一本久久a久久精品综合妖精| 欧美在线黄色| 丝袜喷水一区| 男女午夜视频在线观看| 黑人巨大精品欧美一区二区蜜桃| 天堂俺去俺来也www色官网| 丝袜脚勾引网站| 2021少妇久久久久久久久久久| 色播在线永久视频| 精品少妇内射三级| 91精品伊人久久大香线蕉| 校园人妻丝袜中文字幕| 中国国产av一级| 国产成人啪精品午夜网站| 国产精品亚洲av一区麻豆 | 久久 成人 亚洲| 9色porny在线观看| av天堂久久9| 国产精品 欧美亚洲| 777米奇影视久久| 亚洲一级一片aⅴ在线观看| 精品久久久久久电影网| 国产精品秋霞免费鲁丝片| 欧美日韩一级在线毛片| 两性夫妻黄色片| 免费高清在线观看视频在线观看| 伊人亚洲综合成人网| 欧美精品一区二区大全| 交换朋友夫妻互换小说| 欧美在线一区亚洲| 精品福利永久在线观看| 国产成人精品久久二区二区91 | 亚洲欧美成人综合另类久久久| 欧美日韩视频精品一区| 免费高清在线观看视频在线观看| 视频在线观看一区二区三区| 久久久亚洲精品成人影院| av网站免费在线观看视频| 大话2 男鬼变身卡| 亚洲一码二码三码区别大吗| 最近中文字幕2019免费版| videosex国产| 老汉色av国产亚洲站长工具| av视频免费观看在线观看| 亚洲成人一二三区av| 午夜激情av网站| 欧美国产精品va在线观看不卡| 女性被躁到高潮视频| 女性被躁到高潮视频| 男人爽女人下面视频在线观看| 日本欧美国产在线视频| 亚洲国产精品成人久久小说| 九色亚洲精品在线播放| 久久久久精品性色| 亚洲欧美精品综合一区二区三区| 国产精品一区二区精品视频观看| 欧美日韩精品网址| 97人妻天天添夜夜摸| 最近2019中文字幕mv第一页| 久久人人爽人人片av| 亚洲三区欧美一区| 蜜桃国产av成人99| 久久久久久人人人人人| 黑人巨大精品欧美一区二区蜜桃| 91成人精品电影| 搡老岳熟女国产| 国产日韩欧美亚洲二区| 国产亚洲精品第一综合不卡| 日本91视频免费播放| 天天影视国产精品| 国产av一区二区精品久久| 亚洲国产精品999| 99re6热这里在线精品视频| 国产免费现黄频在线看| 欧美最新免费一区二区三区| 天天添夜夜摸| 男女之事视频高清在线观看 | 久久久久久人妻| 精品亚洲成国产av| 美女中出高潮动态图| 一区二区日韩欧美中文字幕| 另类亚洲欧美激情| 久久精品国产亚洲av涩爱| 一本大道久久a久久精品| 欧美人与性动交α欧美软件| 老司机影院成人| 曰老女人黄片| 国产亚洲av高清不卡| 国产精品一二三区在线看| 在线亚洲精品国产二区图片欧美| 男女无遮挡免费网站观看| 大陆偷拍与自拍| 99香蕉大伊视频| 中国国产av一级| 久久久久网色| 乱人伦中国视频| 国产精品二区激情视频| 成人毛片60女人毛片免费| 亚洲美女视频黄频| 日韩成人av中文字幕在线观看| 狂野欧美激情性bbbbbb| 国产xxxxx性猛交| 波野结衣二区三区在线| 欧美日韩视频精品一区| 亚洲欧洲精品一区二区精品久久久 | 纵有疾风起免费观看全集完整版| 99热网站在线观看| 亚洲色图综合在线观看| 最近的中文字幕免费完整| 亚洲成人免费av在线播放| 岛国毛片在线播放| 91国产中文字幕| 日本猛色少妇xxxxx猛交久久| 卡戴珊不雅视频在线播放| 天堂俺去俺来也www色官网| 日韩欧美一区视频在线观看| 国产av码专区亚洲av| 最黄视频免费看| 黑人欧美特级aaaaaa片| 777久久人妻少妇嫩草av网站| 亚洲成人免费av在线播放| 日韩熟女老妇一区二区性免费视频| 日本av免费视频播放| 男女高潮啪啪啪动态图| 欧美国产精品va在线观看不卡| 日日啪夜夜爽| 亚洲中文av在线| 亚洲天堂av无毛| 久久影院123| 老司机在亚洲福利影院| 国产免费又黄又爽又色| 色网站视频免费| 国产精品秋霞免费鲁丝片| 久久 成人 亚洲| 飞空精品影院首页| 亚洲国产精品999| 日韩大码丰满熟妇| 国产亚洲午夜精品一区二区久久| 亚洲四区av| 国产精品麻豆人妻色哟哟久久| 成人亚洲精品一区在线观看| 成人黄色视频免费在线看| 成人18禁高潮啪啪吃奶动态图| 九色亚洲精品在线播放| 国产精品香港三级国产av潘金莲 | 国产成人精品无人区| 精品午夜福利在线看| 天天躁日日躁夜夜躁夜夜| 亚洲美女搞黄在线观看| 亚洲精品美女久久av网站| 黄片无遮挡物在线观看| 国产日韩一区二区三区精品不卡| 国产精品免费视频内射| avwww免费| 51午夜福利影视在线观看| 亚洲免费av在线视频| 男男h啪啪无遮挡| 999久久久国产精品视频| 国产欧美亚洲国产| 别揉我奶头~嗯~啊~动态视频 | 欧美日韩亚洲综合一区二区三区_| 香蕉国产在线看| 十分钟在线观看高清视频www| 久久久久久久精品精品| 久久精品久久久久久久性| 久久亚洲国产成人精品v| a级毛片在线看网站| 国产精品免费大片| 一边摸一边抽搐一进一出视频| 国产精品久久久av美女十八| 狠狠婷婷综合久久久久久88av| 久久久久视频综合| 亚洲国产最新在线播放| 韩国精品一区二区三区| 人人妻人人澡人人爽人人夜夜| 美女主播在线视频| 高清不卡的av网站| 午夜精品国产一区二区电影| 宅男免费午夜| 亚洲av中文av极速乱| 丰满乱子伦码专区| 久久久久国产精品人妻一区二区| 热99久久久久精品小说推荐| 婷婷色麻豆天堂久久| 国产精品一区二区在线不卡| 日韩 亚洲 欧美在线| 亚洲成人国产一区在线观看 | 十八禁人妻一区二区| 国产精品一区二区在线不卡| 亚洲熟女毛片儿| 国产毛片在线视频| 久久精品亚洲av国产电影网| 久久综合国产亚洲精品| 韩国av在线不卡| 亚洲第一青青草原| 大香蕉久久成人网| 免费观看性生交大片5| 在线观看www视频免费| 精品人妻熟女毛片av久久网站| av免费观看日本| 久久97久久精品| 狠狠精品人妻久久久久久综合| 女性生殖器流出的白浆| 在线免费观看不下载黄p国产| 精品国产一区二区三区四区第35| 中文天堂在线官网| 欧美精品av麻豆av| 建设人人有责人人尽责人人享有的| 尾随美女入室| 成人黄色视频免费在线看| 国产av国产精品国产| 久久久久久人人人人人| 久久久久久久久久久免费av| 99久久精品国产亚洲精品| 色播在线永久视频| videos熟女内射| 中文字幕人妻丝袜制服| 2018国产大陆天天弄谢| 九色亚洲精品在线播放| 久久av网站| 欧美人与性动交α欧美精品济南到| 色播在线永久视频| 国产成人系列免费观看| 肉色欧美久久久久久久蜜桃| 国产激情久久老熟女| 在线免费观看不下载黄p国产| 欧美黑人欧美精品刺激| 久久精品久久久久久噜噜老黄| 亚洲国产欧美网| 国产成人系列免费观看| 亚洲av欧美aⅴ国产| 成人影院久久| 欧美 日韩 精品 国产| 色网站视频免费| 亚洲精品国产一区二区精华液| 最近2019中文字幕mv第一页| 午夜福利,免费看| 观看av在线不卡| 国产成人免费无遮挡视频| 亚洲国产精品国产精品| 国产av码专区亚洲av| 两性夫妻黄色片| 久久韩国三级中文字幕| 亚洲成人免费av在线播放| 欧美精品一区二区大全| 青春草亚洲视频在线观看| 狠狠精品人妻久久久久久综合| 精品国产乱码久久久久久男人| 午夜久久久在线观看| 韩国高清视频一区二区三区| 亚洲国产精品999| 日本一区二区免费在线视频| 国产精品久久久久成人av| 两个人免费观看高清视频| 国产一区二区三区av在线| 精品卡一卡二卡四卡免费| 亚洲国产看品久久| 婷婷色麻豆天堂久久| 国产免费视频播放在线视频| 少妇精品久久久久久久| 久久免费观看电影| 午夜免费鲁丝| 伊人久久国产一区二区| 国产女主播在线喷水免费视频网站| 两个人看的免费小视频| 丝瓜视频免费看黄片| 免费女性裸体啪啪无遮挡网站| 激情五月婷婷亚洲| 在线亚洲精品国产二区图片欧美| 国产精品av久久久久免费| 亚洲欧洲日产国产| 国产成人欧美| 天堂中文最新版在线下载| av一本久久久久| 亚洲国产成人一精品久久久| 国产免费福利视频在线观看| 国产精品一区二区在线不卡| www.自偷自拍.com| 精品免费久久久久久久清纯 | 亚洲第一区二区三区不卡| 国产亚洲最大av| 亚洲国产成人一精品久久久| 男女之事视频高清在线观看 | 国产极品粉嫩免费观看在线| 人妻人人澡人人爽人人| 亚洲人成网站在线观看播放| 婷婷色综合大香蕉| 欧美 亚洲 国产 日韩一| 18禁国产床啪视频网站| 成人国产麻豆网| 下体分泌物呈黄色| 啦啦啦 在线观看视频| 各种免费的搞黄视频| 亚洲三区欧美一区| 欧美日韩综合久久久久久| www.熟女人妻精品国产| 婷婷成人精品国产| 黄频高清免费视频| www.熟女人妻精品国产| www.自偷自拍.com| 久久久亚洲精品成人影院| 国产免费视频播放在线视频| 美女午夜性视频免费| 99香蕉大伊视频| 亚洲欧美中文字幕日韩二区| 99国产精品免费福利视频| 精品午夜福利在线看| 国产一级毛片在线| 国产福利在线免费观看视频| 黄色怎么调成土黄色| 国产又爽黄色视频| 少妇被粗大猛烈的视频| 国产老妇伦熟女老妇高清| 十八禁网站网址无遮挡| 免费久久久久久久精品成人欧美视频| 女人被躁到高潮嗷嗷叫费观| 永久免费av网站大全| 男男h啪啪无遮挡| 欧美变态另类bdsm刘玥| 精品国产一区二区久久| 精品少妇一区二区三区视频日本电影 | 亚洲久久久国产精品| 99热网站在线观看| 亚洲三区欧美一区| 久久人妻熟女aⅴ| 天美传媒精品一区二区| 精品国产一区二区三区四区第35| 免费观看a级毛片全部| 老司机影院成人| 九草在线视频观看| 亚洲欧洲日产国产| 久久精品国产a三级三级三级| 天天躁夜夜躁狠狠久久av| 热re99久久国产66热| 国产精品国产av在线观看| 少妇人妻久久综合中文| 建设人人有责人人尽责人人享有的| 欧美激情高清一区二区三区 | 欧美激情 高清一区二区三区| e午夜精品久久久久久久| 国产精品国产三级国产专区5o| 国产av码专区亚洲av| 黄色怎么调成土黄色| 久久精品aⅴ一区二区三区四区| bbb黄色大片| 国产精品免费视频内射| 在线天堂最新版资源| 精品一区二区三区av网在线观看 | 亚洲精品中文字幕在线视频| 美女脱内裤让男人舔精品视频| 天堂8中文在线网| 欧美av亚洲av综合av国产av | 国产精品蜜桃在线观看| 天天躁夜夜躁狠狠久久av| 精品卡一卡二卡四卡免费| 宅男免费午夜| 亚洲国产欧美在线一区| 伊人久久大香线蕉亚洲五| 国产精品国产三级专区第一集| 亚洲欧美清纯卡通| 国产精品二区激情视频| 欧美日韩亚洲国产一区二区在线观看 | 亚洲av日韩精品久久久久久密 | 哪个播放器可以免费观看大片| 国产精品99久久99久久久不卡 | av卡一久久| 黄片小视频在线播放| 国产一区亚洲一区在线观看| 999久久久国产精品视频| 又大又爽又粗| 欧美国产精品va在线观看不卡| 国产精品一区二区在线观看99| 亚洲国产精品一区三区| 巨乳人妻的诱惑在线观看| 日韩不卡一区二区三区视频在线| 亚洲熟女毛片儿| 午夜福利网站1000一区二区三区| 午夜福利视频精品| 人妻一区二区av| 2018国产大陆天天弄谢| 男女午夜视频在线观看| 中文字幕人妻熟女乱码| 亚洲精品第二区| 午夜福利网站1000一区二区三区| netflix在线观看网站| 欧美黄色片欧美黄色片| 熟女少妇亚洲综合色aaa.| 新久久久久国产一级毛片| 哪个播放器可以免费观看大片| 国产黄色视频一区二区在线观看| 99精国产麻豆久久婷婷| 国产av国产精品国产| 欧美激情 高清一区二区三区| 一本色道久久久久久精品综合| 美女视频免费永久观看网站| 丝袜美腿诱惑在线| 综合色丁香网| 日韩一区二区视频免费看| 亚洲精品国产区一区二| 日本91视频免费播放| 看非洲黑人一级黄片| 国产日韩欧美视频二区| 在线看a的网站| 午夜福利视频精品| 捣出白浆h1v1| 久久午夜综合久久蜜桃| 久久婷婷青草| 亚洲精品乱久久久久久| 99久久人妻综合| 99香蕉大伊视频| 免费看av在线观看网站| 午夜日本视频在线| 欧美 亚洲 国产 日韩一| 亚洲av电影在线进入| 久久精品久久精品一区二区三区| 黄片小视频在线播放| av一本久久久久| 国产淫语在线视频| 亚洲国产欧美在线一区| 精品亚洲乱码少妇综合久久| 观看av在线不卡| 热re99久久精品国产66热6| 日韩一本色道免费dvd| 久久精品亚洲熟妇少妇任你| 亚洲成人av在线免费| 卡戴珊不雅视频在线播放| 不卡视频在线观看欧美| 久久久久人妻精品一区果冻| 午夜福利一区二区在线看| 亚洲国产中文字幕在线视频| 亚洲少妇的诱惑av| 精品久久久久久电影网| 午夜精品国产一区二区电影| 男人操女人黄网站| 热re99久久精品国产66热6| 免费高清在线观看视频在线观看| 国产精品99久久99久久久不卡 | 无限看片的www在线观看| 久久女婷五月综合色啪小说| 久久99热这里只频精品6学生| 亚洲视频免费观看视频| 国产精品欧美亚洲77777| 久久婷婷青草| 一级片'在线观看视频| 在线观看一区二区三区激情| 国产精品亚洲av一区麻豆 | 国产成人91sexporn| 亚洲av成人不卡在线观看播放网 | 男人添女人高潮全过程视频| 亚洲精品乱久久久久久| videos熟女内射| 99久久精品国产亚洲精品| 亚洲av成人不卡在线观看播放网 | 久久精品熟女亚洲av麻豆精品| 777米奇影视久久| 巨乳人妻的诱惑在线观看| xxx大片免费视频| 一级爰片在线观看| 18禁裸乳无遮挡动漫免费视频| 乱人伦中国视频| 亚洲国产av影院在线观看| 亚洲伊人色综图| 精品一区二区免费观看| 亚洲精品一区蜜桃| 性少妇av在线| 成人国产麻豆网| 悠悠久久av| av免费观看日本| 中文字幕制服av| 视频区图区小说| 日韩制服丝袜自拍偷拍| 操出白浆在线播放| 亚洲美女搞黄在线观看| 人人妻人人澡人人看| 亚洲美女搞黄在线观看| bbb黄色大片| 韩国精品一区二区三区| 国产精品秋霞免费鲁丝片| 乱人伦中国视频| 亚洲av男天堂| 一个人免费看片子| 成人国语在线视频| 午夜影院在线不卡| 国产男人的电影天堂91| 欧美97在线视频| 中文精品一卡2卡3卡4更新| 久久精品熟女亚洲av麻豆精品| 欧美日韩亚洲国产一区二区在线观看 | 欧美日韩综合久久久久久| 一级毛片我不卡| 在线天堂最新版资源| 女性被躁到高潮视频| 纵有疾风起免费观看全集完整版| a 毛片基地| 国产精品国产三级专区第一集| 国产亚洲一区二区精品| 香蕉国产在线看| 纯流量卡能插随身wifi吗| 国产麻豆69| 午夜福利视频在线观看免费| 各种免费的搞黄视频| 亚洲国产看品久久| 黑人猛操日本美女一级片| 亚洲,欧美精品.| 久久97久久精品| 亚洲成色77777| 亚洲综合精品二区| 肉色欧美久久久久久久蜜桃| 欧美国产精品va在线观看不卡| 男的添女的下面高潮视频| 亚洲精品第二区| 国产一区二区三区av在线| 黄色 视频免费看| 欧美日韩福利视频一区二区| 亚洲伊人色综图| 日韩中文字幕视频在线看片| 日韩一区二区视频免费看| 中国国产av一级| 久久久久精品人妻al黑| 国产1区2区3区精品| 成人亚洲欧美一区二区av| 超碰成人久久| 亚洲国产欧美一区二区综合| 国产精品一二三区在线看| av福利片在线| 男女午夜视频在线观看| 成人三级做爰电影| 国产av一区二区精品久久| 国产精品女同一区二区软件| www.av在线官网国产| 天天躁日日躁夜夜躁夜夜| 国产成人精品无人区| 亚洲精品成人av观看孕妇| 国产精品无大码| 亚洲av成人不卡在线观看播放网 | 蜜桃在线观看..| 色视频在线一区二区三区| 午夜福利免费观看在线| 中文字幕精品免费在线观看视频| 亚洲av日韩在线播放| 下体分泌物呈黄色| 久久免费观看电影| 国产精品国产三级国产专区5o| 欧美日韩视频精品一区| 九九爱精品视频在线观看| 亚洲av男天堂| 久久精品aⅴ一区二区三区四区| 亚洲情色 制服丝袜| 亚洲中文av在线| 国产精品久久久久成人av| 亚洲av欧美aⅴ国产| 激情视频va一区二区三区| 又大又黄又爽视频免费| 亚洲精品一区蜜桃| 午夜激情久久久久久久| 国产伦理片在线播放av一区| 国产免费视频播放在线视频| 男女之事视频高清在线观看 | 天美传媒精品一区二区| 夫妻性生交免费视频一级片| 不卡av一区二区三区| 婷婷色麻豆天堂久久| 中文字幕制服av| 久久久久视频综合| kizo精华| 国产免费视频播放在线视频| 国产成人精品久久久久久| 国产精品免费视频内射| 哪个播放器可以免费观看大片| 91成人精品电影| 男女边吃奶边做爰视频| 亚洲欧美成人精品一区二区| 91精品国产国语对白视频| av.在线天堂|