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

    移動車輛荷載作用下中低速磁浮大跨度連續(xù)梁橋動力響應(yīng)分析

    2019-03-24 01:25:40朱志輝胡明勛馮典瑾郭向榮劉澍
    關(guān)鍵詞:積分法梁橋阻尼比

    朱志輝,胡明勛,馮典瑾,郭向榮,劉澍

    移動車輛荷載作用下中低速磁浮大跨度連續(xù)梁橋動力響應(yīng)分析

    朱志輝1, 2,胡明勛1,馮典瑾3,郭向榮1, 2,劉澍1, 2

    (1. 中南大學(xué) 土木工程學(xué)院,湖南 長沙 410075;2. 中南大學(xué) 高速鐵路建造技術(shù)國家工程實(shí)驗(yàn)室,湖南 長沙 410075;3. 湖南磁浮交通發(fā)展股份有限公司,湖南 長沙 410075)

    為研究中低速磁浮大跨度連續(xù)梁橋在移動列車荷載作用下的動力響應(yīng),采用移動荷載模擬磁懸浮車輛作用,采用振型分解精細(xì)積分法建立橋梁動力方程并求解;通過與現(xiàn)場試驗(yàn)結(jié)果對比,驗(yàn)證本文模型的準(zhǔn)確性和方法的正確性,對比分析不同行車速度、不同橋梁阻尼比等工況下的橋梁動力響應(yīng)。研究結(jié)果表明:橋梁動力系數(shù)隨列車速度增加而增大,但整體增加幅度不大;在列車以設(shè)計(jì)運(yùn)行速度100?km/h行駛時,橋梁阻尼比的變化對橋梁響應(yīng)影響較小;不同工況下橋梁動力系數(shù)均低于我國《中低速磁浮交通設(shè)計(jì)規(guī)范》中的限值。

    中低速磁??;大跨連續(xù)梁橋;動力響應(yīng);振型分解精細(xì)積分法;現(xiàn)場試驗(yàn)

    中低速磁浮交通作為一種新興的交通方式,因爬坡能力大、轉(zhuǎn)彎半徑小、噪聲低、污染小、環(huán)境兼容性好等優(yōu)點(diǎn)已越來越受到人們的關(guān)注。隨著我國長沙中低速磁浮交通運(yùn)營線的順利開通,國內(nèi)諸多城市亦開始規(guī)劃中低速磁浮交通的建設(shè),中低速磁浮交通具有良好的發(fā)展前景[1]。為了在城市軌道交通中提高其競爭力,磁浮線路大量采用較為輕盈的高架橋梁。磁浮車輛通過橋梁區(qū)段時,橋梁對車輛的懸浮穩(wěn)定性及動力響應(yīng)有顯著的影響,國內(nèi)外學(xué)者針對磁懸浮橋梁振動響應(yīng)問題進(jìn)行了大量的研究。滕延鋒等[2]將磁浮車輛簡化為均布二系懸掛模型,建立了磁浮車輛—軌道梁耦合振動分析模型,并對磁浮列車駛過三跨連續(xù)梁時軌道梁的振動響應(yīng)進(jìn)行了數(shù)值仿真分析;時瑾等[3]建立磁浮列車橋梁耦合動力學(xué)模型,提出求解機(jī)電強(qiáng)耦合微分方程的積分方法,并通過現(xiàn)場測試結(jié)果驗(yàn)證所提出的方法的準(zhǔn)確性。耿杰等[4-5]建立了中低速磁懸浮車橋耦合動力學(xué)模型,通過現(xiàn)場試驗(yàn)驗(yàn)證自編程序的正確性,并結(jié)合現(xiàn)場試驗(yàn)對中低速磁懸浮簡支梁橋的動力系數(shù)進(jìn)行了研究;楊平等[6]建立車?軌道豎向耦合振動分析模型,并利用該模型分析車速、車重及橋梁阻尼比等對20?m簡支梁動力系數(shù)的影響,同時結(jié)合現(xiàn)場試驗(yàn)測試結(jié)果與規(guī)范中動力系數(shù)公式計(jì)算結(jié)果進(jìn)行對比。鄧建良等[7]將磁懸浮列車荷載簡化為移動均布荷載,使用模態(tài)疊加法求解得到了簡支軌道梁的響應(yīng)解析解,并與響應(yīng)數(shù)值解和列車-橋梁耦合振動解做了對比,指出移動均布荷載可以反映軌道梁振動的特征和振動規(guī)律;SONG等[8]將磁懸浮列車簡化為均布移動荷載,并通過對精細(xì)有限元模型的數(shù)值計(jì)算驗(yàn)證了簡化方法的正確性和可行性。由于磁懸浮車輛相對于橋梁質(zhì)量較小,車橋耦合作用不顯著,因此可以采用移動荷載模擬磁懸浮列車對橋梁的動力作用。為提高移動荷載過橋時動力分析的精度和效率,張亞輝等[9]在荷載作用位置隨著時間發(fā)生變化的情形下應(yīng)用PIM方法,并與常規(guī)計(jì)算方法進(jìn)行比較,指出該方法計(jì)算精度及效率較通常的數(shù)值積分方法均得到了顯著提高;桂水榮等[10]根據(jù)模態(tài)綜合疊加技術(shù)的優(yōu)勢提出了基于PIM的車橋耦合振動模型新算法,并以移動常量力作用于簡支梁為例,對比新型算法與Newmark- β法等計(jì)算結(jié)果,指出基于PIM的新型算法求解車橋耦合動力問題時不受積分步長限制,具有快速收斂的優(yōu)勢。本文在以往研究的基礎(chǔ)上,采用移動荷載模擬磁浮車輛荷載作用,建立了精細(xì)三跨連續(xù)梁橋有限元模型,并基于得到的MDPIM(振型分解精細(xì)積分法)建立動力方程,這樣既克服了直接應(yīng)用精細(xì)積分法求解大型結(jié)構(gòu)矩陣尺度太大的困難,又保持了精細(xì)積分法高精度的優(yōu)點(diǎn)[11]。以長沙中低速磁浮瀏陽河大橋?yàn)槔?85 m+110 m+85 m),自編MATLAB計(jì)算程序分析了磁浮列車開行速度、橋梁阻尼比對磁浮橋梁豎向振動響應(yīng)的影響規(guī)律。

    1 車輛、橋梁動力學(xué)模型

    1.1 車輛荷載模型

    由于磁懸浮車輛相對于橋梁質(zhì)量較小,車橋耦合作用不顯著,所以磁懸浮列車采用一組移動荷載模擬。3車編組的磁浮列車每節(jié)車重量相同,將每個轉(zhuǎn)向架與橋梁系統(tǒng)之間電磁力的作用簡化為集中力,對于一節(jié)車而言,每2個集中力的間距為0,2節(jié)車的首尾集中荷載間距為1,最終得到磁浮列車的荷載簡化模型如圖1所示。

    圖1 磁浮列車荷載簡化模型

    1.2 橋梁動力模型

    根據(jù)橋梁不同構(gòu)件的特性,可以采用不同單元建立橋梁的有限元模型。建立橋梁模型后,Ansys可以導(dǎo)出結(jié)構(gòu)的質(zhì)量、阻尼和剛度矩陣。橋梁系統(tǒng)動力方程可寫為:

    根據(jù)正則坐標(biāo)變換可以得到:

    則等式(1)可以改寫為:

    在求解橋梁動力方程時,集中荷載作用到橋梁有限元模型兩相鄰節(jié)點(diǎn)和+1節(jié)點(diǎn)之間的單元上時,可以通過Hermite三次插值將集中荷載分配到單元兩端節(jié)點(diǎn)上。

    2 精細(xì)積分法求解流程

    精細(xì)時程分析法(Precise Integration Method, PIM)是利用矩陣指數(shù)函數(shù)可以在計(jì)算機(jī)字長范圍內(nèi)精確計(jì)算的特點(diǎn)得到很精密的解答[15]。與傳統(tǒng)的逐步積分法即中心差分法、Newmark-法和Wilson-法相比,精細(xì)積分法可以降低對步長的要求,在積分步長較大時得到較為精確的結(jié)果,從而能大大提高數(shù)值積分的計(jì)算效率。其中,積分步長影響計(jì)算的時間,積分間隔影響計(jì)算的精度,所以在計(jì)算時一般采用較大的積分步長和較小的積分間隔。該方法采用了2類的算法,相當(dāng)于在每一時間步長內(nèi)進(jìn)一步劃分成2個精細(xì)步長,大大提高了計(jì)算結(jié)果的精度;此外該方法降低了對于積分步長的要求,在積分步長較大時可以得到較為精確的結(jié)果,從而提高了數(shù)值積分的計(jì)算效率[15]{翟婉明, 2001 #150}。

    從文獻(xiàn)[16]中可以得到有關(guān)PIM計(jì)算方法的求解過程,該方法計(jì)算流程如圖2所示。

    圖2 精細(xì)積分法計(jì)算流程圖

    對于下式所示的結(jié)構(gòu)動力學(xué)方程:

    精細(xì)積分法較宜于處理一階方程,故采用Hamilton體系方法進(jìn)行方程變換[15],將上式化為一階方程:

    1) 齊次方程的求解

    齊次方程的通解可以形式地寫成

    令時間步長?=,則

    2) 指數(shù)矩陣的精細(xì)計(jì)算

    利用加法原理,指數(shù)矩陣的精細(xì)計(jì)算方法可表示為:

    為計(jì)算矩陣,式(11)可表示為:

    3) 非齊次項(xiàng)的處理

    其中:

    從而有:

    3 瀏陽河大橋現(xiàn)場試驗(yàn)

    3.1 橋梁概況

    長沙中低速磁懸浮工程從長沙南站到長沙黃花國際機(jī)場線路全長約18.643 km,設(shè)計(jì)行車速度100 km/h。全線初期設(shè)車站3座,高架區(qū)間總長度16.795 km,占86.3%。其中瀏陽河段大跨連續(xù)梁橋是一座三跨變截面預(yù)應(yīng)力混凝土雙線(直線)連續(xù)梁橋,主梁長度為279.7 m。該連續(xù)梁橋采用梁上承軌方案,承軌梁高為0.9 m,頂板和腹板厚度為22 cm;下層箱梁橫截面為單箱單室直腹板,箱梁頂面寬度8.1 m,底面寬度5.7 m,頂板厚度為40 cm,腹板厚度分別為40 cm和80 cm。箱梁在中支點(diǎn)截面中心線處梁高為8.00 m,跨中截面中心線處梁高為5.00 m。邊支座中心線至梁端0.70 m,邊支座橫橋向中心距4.50 m,中支座橫橋向中心距5.00 m。三跨連續(xù)梁橋示意圖及斷面圖如圖3所示,其中1-1為中支點(diǎn)斷面,2-2為中跨跨中斷面。

    3.2 試驗(yàn)方案

    3.2.1 地基微波雷達(dá)概況

    近年來地基微波雷達(dá)開始應(yīng)用于橋梁的撓度測量等領(lǐng)域之中,相對于其他傳統(tǒng)測量方式而言具有非接觸式遠(yuǎn)程測量、精度高、環(huán)境適應(yīng)性強(qiáng)以及支持多點(diǎn)同時測量等優(yōu)點(diǎn)[17]。

    單位:cm

    中南大學(xué)高速鐵路國家工程實(shí)驗(yàn)室和國防科技大學(xué)基于高分辨雷達(dá)遙感成像與雷達(dá)干涉原理聯(lián)合開發(fā)了地基微波雷達(dá),該設(shè)備采用電磁波遠(yuǎn)程對橋梁、高層建筑、鐵塔等結(jié)構(gòu)的微小變形進(jìn)行高速和實(shí)時的監(jiān)測,廣泛應(yīng)用于靜撓度、動撓度及斜拉橋索力等測量領(lǐng)域,是一種具有較高精確度的新型干涉雷達(dá)測試系統(tǒng)。不同于傳統(tǒng)的雷達(dá)干涉儀采用步進(jìn)頻率(SF-CW)波形,該地基微波雷達(dá)采用線性調(diào)頻連續(xù)波(LFMCW)[18]。LFMCW雷達(dá)具有低截獲概率(LPI)特性,測試精度受目標(biāo)物體運(yùn)動速度影響較小,波束抗干擾能力強(qiáng)等特點(diǎn);同時易于擴(kuò)頻和實(shí)現(xiàn)高分辨率,還具有無距離盲區(qū)和大時間帶寬的特性,能夠進(jìn)行多目標(biāo)精準(zhǔn)識別[19]。

    相位干涉測量技術(shù)是一項(xiàng)用途十分廣泛的測量技術(shù)。應(yīng)用該技術(shù),雷達(dá)主機(jī)可通過對相同檢測區(qū)域的重復(fù)觀測,根據(jù)回波信號相位差?反演生成目標(biāo)形變圖。設(shè)為雷達(dá)發(fā)射波波長,則目標(biāo)位移量為:

    式中:為電磁波波速;為地基微波雷達(dá)的工作頻率。在動態(tài)變形測試時,位移測試精度優(yōu)于0.1?mm,還可通過進(jìn)一步調(diào)節(jié)工作頻段,獲得更短的發(fā)射波波長,以提高測試精度。

    圖4 位移投影技術(shù)豎直投影方向

    3.2.2 試驗(yàn)測點(diǎn)布置

    瀏陽河大橋現(xiàn)場試驗(yàn)使用地基微波雷達(dá)系統(tǒng)測量列車過橋時,橋梁主要測點(diǎn)的豎向位移響應(yīng)。測點(diǎn)布置在連續(xù)梁每一跨的跨中位置,如圖5所示。為提高測點(diǎn)處雷達(dá)散射質(zhì)量進(jìn)而提高測量精度,在橋梁測點(diǎn)處布置了雷達(dá)標(biāo)靶。

    單位:m

    4 橋梁有限元模型及驗(yàn)證

    4.1 橋梁有限元模型

    圖6為瀏陽河大橋有限元模型。其中主梁、承軌梁及橋墩均采用BEAM188單元,主橋和承軌梁劃分單元長度為0.4?m,彈性模量均為5×1010N/m。承軌梁與主梁之間使用剛臂連接;支座采用COMBIN14彈簧單元模擬。該橋主梁、承軌梁及橋墩密度均為2 650?kg/m3,橋上二期恒載為55?kN/m,在進(jìn)行動力分析時把二恒等效為橋梁均布質(zhì)量施加到主橋梁單元之中。

    單位:m

    4.2 模態(tài)驗(yàn)證

    為了驗(yàn)證建立的橋梁有限元模型的正確性,以及自編MATLAB計(jì)算程序的可靠性,對比瀏陽河大橋自振特性的理論值和實(shí)測值。表1給出了有限元模型前4階振型和頻率(試驗(yàn)實(shí)測值只提供了前兩階頻率),圖7所示為有限元模型前4階振型圖。由表1可知,瀏陽河大橋有限元模型經(jīng)過模態(tài)求解得到的前2階頻率與試驗(yàn)實(shí)測結(jié)果接近。

    表1 橋梁頻率理論值和試驗(yàn)值

    4.3 靜載作用

    靜力試驗(yàn)按照2列車雙線布置的方式,每列磁懸浮列車采用3車編組,每節(jié)車滿載重量為33 t。列車荷載簡化為圖1所示的荷載形式,其中0為2.8?m,1為4.4?m,集中荷載為64.68?kN,靜載試驗(yàn)時,列車荷載分別沿各跨的跨中截面對稱分布,如圖8所示。

    (a) 1階振型;(b) 2階振型;(c) 3階振型;(d) 4階振型

    (a) 縱橋向靜力加載圖;(b) 橫橋向靜力加載圖

    工況1,工況2和工況3分別為列車荷載沿第1跨,第2跨和第3跨的跨中截面對稱分布,各跨的跨中截面位于圖5中測點(diǎn)a,b和c位置處,通過現(xiàn)場實(shí)測得到各測點(diǎn)位置的撓度值。對比3種工況下各個測點(diǎn)的撓度值,結(jié)果如表2所示。(“—”表示現(xiàn)場試驗(yàn)實(shí)測結(jié)果中未提供相關(guān)結(jié)果,正值表示橋梁位移向下,負(fù)值表示橋梁位移向上)。

    由表2可知,靜力作用加載時3種工況下測點(diǎn)的現(xiàn)場試驗(yàn)實(shí)測撓度值與使用自編程序計(jì)算得到的測點(diǎn)撓度值較為接近,其誤差均小于5%。

    表2 跨中撓度試驗(yàn)值與理論值

    圖9 試驗(yàn)實(shí)測值與理論值對比圖

    4.4 動載作用

    動力加載試驗(yàn)為一列磁浮列車以60?km/h的速度通過三跨連續(xù)梁橋,采用圖1所示的移動荷載模擬列車軸重,每節(jié)車載重為33 t。橋梁阻尼比取值為0.02,計(jì)算積分步長為0.024?s,計(jì)算時間為20?s。針對中跨跨中測點(diǎn)b,采用振型分解精細(xì)積分法,根據(jù)高速鐵路歐洲規(guī)范中橋梁加速度限值對應(yīng)的頻率[20],取前100階振型進(jìn)行求解得到理論結(jié)果,將試驗(yàn)實(shí)測結(jié)果與理論結(jié)果進(jìn)行對比,如圖9所示。在動力作用加載工況下通過現(xiàn)場試驗(yàn)實(shí)測得到的中間跨跨中測點(diǎn)的位移時程曲線與使用自編程序計(jì)算得到的位移時程曲線吻合很好。

    綜上所述,瀏陽河大橋有限元模型經(jīng)過模態(tài)驗(yàn)證、靜載試驗(yàn)、動載試驗(yàn)得到的結(jié)果均與試驗(yàn)實(shí)測值較為吻合,表明了瀏陽河大橋有限元模型的準(zhǔn)確性以及自編MATLAB計(jì)算程序的正確性。

    5 大跨度連續(xù)梁橋動力分析

    當(dāng)列車通過橋梁時,在動力荷載的作用下橋梁結(jié)構(gòu)的豎向動力響應(yīng)大于列車荷載靜止在橋上所引起的靜力響應(yīng),在假定車輛滿載條件下(每節(jié)車輛滿載重量為33 t),本節(jié)分別研究車速及橋梁阻尼比對橋梁豎向動力響應(yīng)的影響規(guī)律。

    5.1 車速的影響

    磁浮車輛速度變化范圍為20~120?km/h,變化間隔為20?km/h。動力分析時間積分步長為0.024?s,阻尼比為0.02。圖10和圖11所示分別為不同車速下橋梁中間跨的跨中測點(diǎn)的位移和加速度時程 曲線。

    圖10 不同速度下中跨跨中位移時程曲線

    從圖10~11中可以看出,不同速度工況下,三跨連續(xù)梁橋的位移時程曲線吻合較好,且跨中撓度峰值和加速度時程變化均不顯著,表明在運(yùn)營速度范圍內(nèi)不會引起橋梁響應(yīng)的劇烈變化。

    圖11 不同速度下中跨跨中加速度時程曲線

    在設(shè)定車速變化范圍內(nèi),三跨連續(xù)梁中跨跨中的動力系數(shù)隨速度變化的點(diǎn)線圖如圖12所示。中跨跨中動力系數(shù)大致隨速度的提高而增加,當(dāng)車速為120?km/h時,最大值僅為1.010 3。

    圖12 不同速度工況下中跨跨中動力系數(shù)

    5.2 阻尼比的影響

    為研究橋梁阻尼對沖擊效應(yīng)的影響,在設(shè)計(jì)時速100?km/h條件下,對比阻尼比介于0.02~0.05時,阻尼比對橋梁中跨跨中的動力系數(shù)的影響規(guī)律。阻尼比分別為0.02,0.03,0.04和0.05時,瀏陽河大橋中跨跨中的動力系數(shù)分別為1.002 5,1.004 2,1.005 2和1.005 9,如圖13所示。從圖中可以看出,動力系數(shù)變化非常小,阻尼比的變化對橋梁動力系數(shù)影響較小。

    圖13 不同阻尼比工況下中跨跨中動力系數(shù)

    6 結(jié)論

    1) 當(dāng)車輛以20~120?km/h的速度通過三跨連續(xù)梁橋時,橋梁跨中位移響應(yīng)沒有發(fā)生顯著變化,橋梁動力系數(shù)隨列車運(yùn)行速度的提高而有所增加,但增幅不大,表明磁浮列車在正常運(yùn)行速度范圍內(nèi),速度變化不會引起橋梁響應(yīng)的劇烈變化。

    2) 車輛以設(shè)計(jì)運(yùn)營速度100 km/h駛過三跨連續(xù)梁橋時,橋梁阻尼比的變化對連續(xù)梁橋動力響應(yīng)的影響較小。

    3) 當(dāng)車輛以20~120?km/h的速度過橋時,以及橋梁阻尼比在0.02~0.05之間變化時,動力系數(shù)均小于我國《中低速磁浮交通設(shè)計(jì)規(guī)范》中規(guī)定的動力系數(shù)最小限值1.15。

    [1] 鄧小星. 中低速磁浮車輛系統(tǒng)動力學(xué)性能研究[D]. 成都: 西南交通大學(xué), 2009. DENG Xiaoxing. Study on dynamic performance of medium and low speed maglev vehicle system[D]. Chengdu: Southwest jiaotong University, 2009.

    [2] 滕延鋒, 滕念管, 寇新建. 承受移動均布二系懸掛模型的磁浮三跨連續(xù)梁振動分析[J]. 振動與沖擊, 2008, 27(3): 120?123. TENG Yanfeng, TENG Nianguan, KOU Xinjian. Vibration analysis of maglev three-span continuous beam subjected to moving two-line suspension model[J]. Journal of Vibration and Shock, 2008, 27(3): 120?123.

    [3] 時瑾, 魏慶朝, 萬傳風(fēng), 等. 隨機(jī)不平順激勵下磁浮車輛軌道梁動力響應(yīng)[J]. 力學(xué)學(xué)報(bào), 2006, 38(6): 850?857. SHI Jin, WEI Qingchao, WAN Chuanfeng, et al. Study on dynamic responses of high-speed maglev vehicle/ guideway coupling system under random irregularity[J]. Theoretical and Applied Mechanics, 2006, 38(6): 850? 857.

    [4] 耿杰, 王黨雄, 李小珍, 等. 中低速磁浮列車-簡支梁系統(tǒng)耦合振動試驗(yàn)研究[J]. 鐵道學(xué)報(bào), 2018(2): 117? 124. GENG Jie, WANG Dangxiong, LI Xiaozhen, et al. Experimental study on coupled vibration of low-medium speed maglev train and simply supported girder system [J]. Journal of the China Railway Society, 2018(2): 117?124.

    [5] 李小珍, 洪沁燁, 耿杰, 等. 中低速磁浮列車-軌道梁豎向耦合模型與驗(yàn)證[J]. 鐵道工程學(xué)報(bào), 2015, 32(9): 103?108. LI Xiaozhen, HONG Qinye, GENG Jie, et al. The medium and low speed maglev train-track beam coupling vibration model and verification[J]. Journal of Railway Engineering Society, 2015, 32(9): 103?108.

    [6] 楊平, 劉德軍, 李小珍. 中低速磁浮簡支軌道梁動力系數(shù)研究[J]. 橋梁建設(shè), 2016, 46(4): 79?84. YANG Ping, LIU Dejun, LI Xiaozhen. Investigation of dynamic factors of low and medium speed maglev simply-supported guideway beam[J]. Bridge Construction, 2016, 46(4): 79?84.

    [7] 鄧建良, 吳定俊, 李奇. 移動均布荷載作用下簡支梁橋動力系數(shù)分析[J]. 工程力學(xué), 2013, 30(5): 56?62. DENG Jianliang, WU Dingjun, LI Qi. Impact factor analysis of simply-supported girder bridges traversed by uniformly-distributed moving loads[J]. Engineering Mechanics, 2013, 30(5): 56?62

    [8] SONG X D, WU D J, LI Q. Dynamic impact analysis of double-tower cable-stayed maglev bridges using a simple model[J]. Journal of Bridge Engineering, 2014, 19(1): 34?43.

    [9] 張亞輝, 張守云, 趙巖, 等. 橋梁受移動荷載動力響應(yīng)的一種精細(xì)積分法[J]. 計(jì)算力學(xué)學(xué)報(bào), 2006, 23(3): 290?294. ZHANG Yahui, ZHANG Shouyun, ZHAO Yan, et al. A precise integration method for bridges subjects to moving loads[J]. Chinese Journal of Computational Mechanics, 2006, 23(3): 290?294.

    [10] 桂水榮, 萬水, 陳水生. 移動車輛荷載過橋耦合振動精細(xì)積分算法[J]. 建筑科學(xué)與工程學(xué)報(bào), 2016, 33(2): 56?62. GUI Shuirong,WAN Shui,CHEN Shuisheng.Precise integration method for bridge subjected to moving vehicle load[J]. Journal of Architecture and Civil Engineering, 2016, 33(2): 56?62.

    [11] 儲德文, 王元豐. 結(jié)構(gòu)動力方程的振型分解精細(xì)積分法[J]. 鐵道學(xué)報(bào), 2003, 25(6): 89?92. CHU Dewen, WANG Yuanfeng. Mode decomposition precision time integration method for structural dynamics [J]. Journal of the China Railway Society, 2003, 25(6): 89?92.

    [12] ZHANG N, XIA H. Dynamic analysis of coupled vehicle–bridge system based on inter-system iteration method[M]. Pergamon Press, Inc, 2013.

    [13] R.克拉夫, J.彭津. 結(jié)構(gòu)動力學(xué)[M]. 北京: 高等教育出版社, 2006. Clough Ray W, Penzien J. Dynamics of structures[M]. Beijing: Higher Education Press, 2006.

    [14] Nguyen D V, Kim K D, Warnitchai P. Simulation procedure for vehicle–substructure dynamic interactions and wheel movements using linearized wheel–rail interfaces[J]. Finite Elements in Analysis & Design, 2009, 45(5): 341?356.

    [15] 鐘萬勰. 暫態(tài)歷程的精細(xì)計(jì)算方法[J]. 計(jì)算力學(xué)學(xué)報(bào), 1995, 12(1): 1?6. ZHONG Wanxie. Precise computation for transient analysis[J]. Chinese Journal of Computational Mechanics, 1995, 12(1): 1?6.

    [16] 鐘萬勰. 結(jié)構(gòu)動力方程的精細(xì)時程積分法[J]. 大連理工大學(xué)學(xué)報(bào), 1994(2): 131?136. ZHONG Wanxie. On precise time-integration method for structural dynamics[J]. Journal of Dalian University of Technology, 1994(2): 131?136.

    [17] 徐亞明, 王鵬, 周校, 等. 地基干涉雷達(dá)IBIS-S橋梁動態(tài)形變監(jiān)測研究[J]. 武漢大學(xué)學(xué)報(bào)(信息科學(xué)版), 2013, 38(7): 845?849. XU Yaming, WANG Peng, ZHOU Xiao, et al. Research on dynamics deformation monitoring of bridges using ground-based Interferometric radar IBIS-S[J]. Geomatics & Information Science of Wuhan University, 2013, 38(7): 845?849.

    [18] WANG Guochao, Munoz-Ferreras J M, GU Changzhan, et al. Application of linear-frequency-modulated continuous-wave (LFMCW) radars for tracking of vital signs[J]. IEEE Transactions on Microwave Theory & Techniques, 2014, 62(6): 1387?1399.

    [19] 陳偉民, 李存龍. 基于微波雷達(dá)的位移/距離測量技術(shù) [J]. 電子測量與儀器學(xué)報(bào), 2015, 29(9): 1251?1265. CHEN Weimin, LI Cunlong. Radar-based displacement/ distance measuring techniques[J]. Journal of Electronic Measurement & Instruction, 2015, 29(9): 1251?1265.

    [20] 李星新, 田啟賢, 鄭平偉, 等. 大跨度高速鐵路橋梁測試與分析中的技術(shù)問題[J]. 橋梁建設(shè), 2009(增2): 89? 92. LI Xingxin, TIAN Qixian, ZHENG Pingwei, et al. Technical issues in testing analysis of long span high-speed railway bridges[J]. Bridge Construction, 2009(Suppl 2): 89?92.

    Dynamic response analysis of low and medium speed maglev long-span continuous beam bridge traversed by moving train load

    ZHU Zhihui1, 2, HU Mingxun1, FENG Dianjin3, GUO Xiangrong1, 2, LIU Shu1, 2

    (1. School of Civil Engineering, Central South University, Changsha 410075, China; 2. National Engineering Laboratory for High Speed Railway Construction, Central South University, Changsha 410075, China; 3. Hunan Maglev Transportation Development Co., Ltd, Changsha 410075, China)

    In order to study the dynamic responses of the low and medium speed maglev long-span continuous beam bridge under the action of the moving train load, the action of maglev train was simulated by the moving load. The dynamic equation of bridge was established and was solved by MDPIM (Mode Decomposition Precision Time Integration Method), through comparison with field test, the accuracy of the finite element model and correctness of method were verified. At last, the dynamic response of bridge under different driving speeds and different damping ratios of the bridge were compared and analyzed. The results show that the dynamic coefficient increases as the train speed increase on the whole, but the increase is not modest. When the train is traveling at the design speed of 100?km/h, the change of the damping ratio of the bridge has little impact on the bridge response. And under different conditions the dynamic coefficient values are lower than the limit value in Code for Design of Low and Medium Speed Maglev Transit.

    low and medium speed maglev; long-span continuous beam bridge; dynamic response; mode decomposition precision time integration method; field test

    U237;U44

    A

    1672 ? 7029(2019)07? 1695 ? 09

    10.19713/j.cnki.43?1423/u.2019.07.013

    2018?10?24

    國家自然科學(xué)基金資助項(xiàng)目(51678576);湖南省科技重大專項(xiàng)項(xiàng)目(2015GK1001)

    朱志輝(1979?),男,河南潢川人,教授,博士,從事車?線?橋耦合動力學(xué)研究;E?mail:zzhh0703@163.com

    (編輯 涂鵬)

    猜你喜歡
    積分法梁橋阻尼比
    基于細(xì)觀結(jié)構(gòu)的原狀黃土動彈性模量和阻尼比試驗(yàn)研究
    地震研究(2021年1期)2021-04-13 01:05:24
    黏滯阻尼器在時程分析下的附加有效阻尼比研究
    振動與沖擊(2019年4期)2019-02-22 02:33:34
    波形分析法求解公路橋梁阻尼比的探討
    上海公路(2018年3期)2018-03-21 05:55:48
    巧用第一類換元法求解不定積分
    結(jié)構(gòu)構(gòu)件阻尼比對大跨度懸索橋地震響應(yīng)的影響
    隨機(jī)結(jié)構(gòu)地震激勵下的可靠度Gauss-legendre積分法
    π型鋼板應(yīng)用于板梁橋加固的非線性研究
    基于積分法的軸對稱拉深成形凸緣區(qū)應(yīng)力、應(yīng)變數(shù)值解
    探討不定積分分部積分法
    河南科技(2014年15期)2014-02-27 14:12:50
    基于攝動法的梁橋開裂后模態(tài)分析
    亚洲精品456在线播放app | 久久天躁狠狠躁夜夜2o2o| 给我免费播放毛片高清在线观看| 亚洲男人的天堂狠狠| 国产精品不卡视频一区二区 | 国产人妻一区二区三区在| 精品午夜福利视频在线观看一区| 欧美激情久久久久久爽电影| 国产午夜精品久久久久久一区二区三区 | 高清毛片免费观看视频网站| 国产在线精品亚洲第一网站| 国产精品国产高清国产av| 久久人人精品亚洲av| 综合色av麻豆| 国产精品亚洲av一区麻豆| 美女被艹到高潮喷水动态| 小说图片视频综合网站| 在线观看一区二区三区| 搡老妇女老女人老熟妇| 校园春色视频在线观看| 小蜜桃在线观看免费完整版高清| 在线播放无遮挡| 精品人妻一区二区三区麻豆 | 亚洲欧美清纯卡通| 亚洲精品影视一区二区三区av| 91久久精品国产一区二区成人| 日本一二三区视频观看| 一本久久中文字幕| 三级毛片av免费| 两人在一起打扑克的视频| 久久精品国产99精品国产亚洲性色| 内地一区二区视频在线| 国产高清三级在线| 小蜜桃在线观看免费完整版高清| 成人高潮视频无遮挡免费网站| 亚洲成人免费电影在线观看| 麻豆久久精品国产亚洲av| 无遮挡黄片免费观看| www.色视频.com| 午夜a级毛片| 中出人妻视频一区二区| 欧美黄色淫秽网站| 夜夜夜夜夜久久久久| 91午夜精品亚洲一区二区三区 | 免费在线观看亚洲国产| 99精品久久久久人妻精品| 男人狂女人下面高潮的视频| 欧美国产日韩亚洲一区| 国产精品久久电影中文字幕| 色尼玛亚洲综合影院| 亚洲,欧美,日韩| 亚洲电影在线观看av| 国产爱豆传媒在线观看| 日本免费一区二区三区高清不卡| 搡老岳熟女国产| 欧美精品国产亚洲| 免费观看人在逋| 国产淫片久久久久久久久 | 亚洲国产精品合色在线| 国产一区二区三区视频了| 亚洲成人精品中文字幕电影| 九九久久精品国产亚洲av麻豆| 深夜精品福利| 国产精品三级大全| 亚洲七黄色美女视频| 天天躁日日操中文字幕| 美女cb高潮喷水在线观看| 免费看美女性在线毛片视频| 俄罗斯特黄特色一大片| 午夜免费激情av| 国产精品伦人一区二区| 成熟少妇高潮喷水视频| 老师上课跳d突然被开到最大视频 久久午夜综合久久蜜桃 | 少妇的逼水好多| 看免费av毛片| 好看av亚洲va欧美ⅴa在| 两性午夜刺激爽爽歪歪视频在线观看| 激情在线观看视频在线高清| 搡老熟女国产l中国老女人| 日韩欧美在线乱码| 看免费av毛片| 国产在线精品亚洲第一网站| 精品一区二区三区av网在线观看| 日本在线视频免费播放| 级片在线观看| 人人妻人人看人人澡| 老熟妇乱子伦视频在线观看| 日本免费一区二区三区高清不卡| 亚洲av成人不卡在线观看播放网| 欧美三级亚洲精品| 一区二区三区激情视频| 精品人妻偷拍中文字幕| 最好的美女福利视频网| 亚洲精品一区av在线观看| 长腿黑丝高跟| 丰满乱子伦码专区| 最后的刺客免费高清国语| 日韩有码中文字幕| 乱码一卡2卡4卡精品| 国产高清视频在线播放一区| 村上凉子中文字幕在线| 高清日韩中文字幕在线| 我要看日韩黄色一级片| 久久性视频一级片| 亚洲欧美精品综合久久99| 欧美高清性xxxxhd video| 男女那种视频在线观看| 色综合婷婷激情| .国产精品久久| 欧美乱色亚洲激情| 亚洲午夜理论影院| 99热这里只有是精品50| 乱码一卡2卡4卡精品| 日日摸夜夜添夜夜添小说| 国内久久婷婷六月综合欲色啪| 国产色爽女视频免费观看| 国产精品乱码一区二三区的特点| 国产精品一区二区性色av| 亚洲综合色惰| 69av精品久久久久久| 精品欧美国产一区二区三| 精品久久久久久久久av| 亚洲熟妇中文字幕五十中出| 成人无遮挡网站| 校园春色视频在线观看| 搞女人的毛片| 久久人人精品亚洲av| 欧美精品国产亚洲| 亚洲欧美清纯卡通| 久久人妻av系列| 激情在线观看视频在线高清| 日本与韩国留学比较| 国产欧美日韩一区二区三| 国产国拍精品亚洲av在线观看| 我要看日韩黄色一级片| 好男人在线观看高清免费视频| 久久精品国产99精品国产亚洲性色| 免费观看精品视频网站| 一区二区三区高清视频在线| 国产精品不卡视频一区二区 | 精品日产1卡2卡| 别揉我奶头 嗯啊视频| 欧美高清成人免费视频www| 欧美+亚洲+日韩+国产| АⅤ资源中文在线天堂| 一个人观看的视频www高清免费观看| 国产成人aa在线观看| 国产成人欧美在线观看| 欧美高清性xxxxhd video| 色综合婷婷激情| 不卡一级毛片| 老司机午夜十八禁免费视频| 国产精品一区二区三区四区久久| 成人特级av手机在线观看| 国产精品久久久久久人妻精品电影| 嫩草影院入口| 性欧美人与动物交配| 亚洲精品乱码久久久v下载方式| 别揉我奶头 嗯啊视频| av黄色大香蕉| 人人妻人人看人人澡| 麻豆成人av在线观看| 757午夜福利合集在线观看| 精品一区二区三区人妻视频| 少妇人妻一区二区三区视频| 亚洲自偷自拍三级| 国产精品女同一区二区软件 | 亚洲熟妇熟女久久| 日韩欧美免费精品| 最后的刺客免费高清国语| 可以在线观看的亚洲视频| 特大巨黑吊av在线直播| 久久久精品大字幕| 91九色精品人成在线观看| 久久精品91蜜桃| 91字幕亚洲| 国产精品,欧美在线| 国产高清有码在线观看视频| 午夜精品久久久久久毛片777| 国产色爽女视频免费观看| 又紧又爽又黄一区二区| 丁香欧美五月| 俄罗斯特黄特色一大片| 三级毛片av免费| 国产欧美日韩精品一区二区| 日韩中文字幕欧美一区二区| 99久久精品热视频| 久久6这里有精品| 色视频www国产| 中文字幕高清在线视频| 成人毛片a级毛片在线播放| 久久人人爽人人爽人人片va | 别揉我奶头 嗯啊视频| 午夜福利高清视频| 亚洲精品久久国产高清桃花| 欧美日韩国产亚洲二区| 久久人人爽人人爽人人片va | 日韩欧美国产一区二区入口| 精品久久久久久久末码| 欧美激情国产日韩精品一区| 国产国拍精品亚洲av在线观看| 蜜桃亚洲精品一区二区三区| 欧美三级亚洲精品| 岛国在线免费视频观看| 99久久99久久久精品蜜桃| 热99在线观看视频| 久久久国产成人精品二区| 麻豆成人av在线观看| 中文字幕熟女人妻在线| 天堂影院成人在线观看| 久久精品国产自在天天线| 久9热在线精品视频| 天天一区二区日本电影三级| 成人特级av手机在线观看| 亚洲不卡免费看| 在线观看av片永久免费下载| 嫁个100分男人电影在线观看| 日韩欧美在线乱码| АⅤ资源中文在线天堂| 级片在线观看| 亚洲乱码一区二区免费版| 日本一二三区视频观看| 一个人看的www免费观看视频| 久久精品综合一区二区三区| 精品一区二区三区av网在线观看| www日本黄色视频网| 亚洲最大成人手机在线| 国语自产精品视频在线第100页| 国产亚洲欧美在线一区二区| 色尼玛亚洲综合影院| 99国产极品粉嫩在线观看| 18美女黄网站色大片免费观看| 赤兔流量卡办理| 久久99热这里只有精品18| 美女 人体艺术 gogo| 在线观看66精品国产| 欧美性猛交黑人性爽| 直男gayav资源| 久久欧美精品欧美久久欧美| 亚洲国产色片| 免费av不卡在线播放| 听说在线观看完整版免费高清| 国产又黄又爽又无遮挡在线| 国产高清激情床上av| 可以在线观看的亚洲视频| 亚洲专区国产一区二区| 深爱激情五月婷婷| 赤兔流量卡办理| 色吧在线观看| 精品久久久久久成人av| 色在线成人网| 亚洲欧美激情综合另类| 午夜福利高清视频| 欧美色视频一区免费| 黄色一级大片看看| 一个人免费在线观看电影| 夜夜躁狠狠躁天天躁| 婷婷色综合大香蕉| 精品人妻熟女av久视频| 亚洲第一区二区三区不卡| 精品免费久久久久久久清纯| 99久久99久久久精品蜜桃| 五月玫瑰六月丁香| 在线播放国产精品三级| 村上凉子中文字幕在线| 国产日本99.免费观看| 国产极品精品免费视频能看的| 一夜夜www| 窝窝影院91人妻| 国产高清激情床上av| 色视频www国产| 国产高清视频在线播放一区| 日韩亚洲欧美综合| 婷婷精品国产亚洲av| 国产精品女同一区二区软件 | 一区二区三区高清视频在线| av在线天堂中文字幕| 久久草成人影院| 看片在线看免费视频| 久久精品91蜜桃| 国产91精品成人一区二区三区| 最近视频中文字幕2019在线8| 99热这里只有是精品在线观看 | 国内揄拍国产精品人妻在线| av在线天堂中文字幕| 非洲黑人性xxxx精品又粗又长| 国产亚洲精品综合一区在线观看| 久久亚洲精品不卡| 国产精品久久久久久久久免 | 青草久久国产| 在线a可以看的网站| 一级毛片久久久久久久久女| 少妇的逼水好多| 精品午夜福利在线看| 中文字幕免费在线视频6| 亚洲乱码一区二区免费版| 国产成人aa在线观看| 午夜福利在线观看免费完整高清在 | 99久久精品国产亚洲精品| 1000部很黄的大片| 欧美色视频一区免费| 在线观看免费视频日本深夜| 成人午夜高清在线视频| 国产熟女xx| 久久精品综合一区二区三区| 看片在线看免费视频| 亚洲国产欧洲综合997久久,| 淫妇啪啪啪对白视频| www.色视频.com| 亚洲av不卡在线观看| 欧美日本亚洲视频在线播放| 午夜激情欧美在线| 日韩欧美国产一区二区入口| x7x7x7水蜜桃| 97碰自拍视频| 免费在线观看成人毛片| 欧美+亚洲+日韩+国产| 麻豆成人午夜福利视频| 一区二区三区高清视频在线| 成人性生交大片免费视频hd| 悠悠久久av| 久久久久九九精品影院| 在线天堂最新版资源| 91久久精品国产一区二区成人| 久久精品国产亚洲av涩爱 | 网址你懂的国产日韩在线| 国内毛片毛片毛片毛片毛片| 在线国产一区二区在线| 成人特级av手机在线观看| 成人一区二区视频在线观看| 亚洲精品粉嫩美女一区| 真人做人爱边吃奶动态| 成人特级av手机在线观看| 可以在线观看毛片的网站| 国产高清视频在线观看网站| 99热精品在线国产| 九色成人免费人妻av| 国内精品久久久久精免费| 亚洲av.av天堂| 直男gayav资源| 国产高清激情床上av| 一卡2卡三卡四卡精品乱码亚洲| 2021天堂中文幕一二区在线观| 亚洲成人久久爱视频| 亚洲 欧美 日韩 在线 免费| 一级毛片久久久久久久久女| 久久精品国产亚洲av天美| 三级男女做爰猛烈吃奶摸视频| 他把我摸到了高潮在线观看| 别揉我奶头~嗯~啊~动态视频| 亚洲最大成人中文| 亚洲激情在线av| 可以在线观看的亚洲视频| 色综合婷婷激情| 少妇人妻一区二区三区视频| 午夜福利在线观看吧| 美女被艹到高潮喷水动态| 久久精品国产自在天天线| 国产真实乱freesex| 亚洲,欧美,日韩| 亚洲国产欧洲综合997久久,| 国产成人欧美在线观看| 麻豆久久精品国产亚洲av| 成人午夜高清在线视频| 人人妻人人看人人澡| 免费电影在线观看免费观看| 老司机午夜十八禁免费视频| 成人特级黄色片久久久久久久| 一级毛片久久久久久久久女| 亚洲在线自拍视频| 精华霜和精华液先用哪个| 麻豆国产av国片精品| 精品福利观看| a级毛片a级免费在线| www.熟女人妻精品国产| 色在线成人网| 可以在线观看的亚洲视频| 亚洲电影在线观看av| 女人十人毛片免费观看3o分钟| 人妻制服诱惑在线中文字幕| 亚洲五月天丁香| 国产精品98久久久久久宅男小说| 18禁裸乳无遮挡免费网站照片| 变态另类成人亚洲欧美熟女| 毛片女人毛片| 国产乱人视频| 成年免费大片在线观看| 五月伊人婷婷丁香| 内地一区二区视频在线| 丰满人妻一区二区三区视频av| 国产大屁股一区二区在线视频| 在线免费观看的www视频| 淫妇啪啪啪对白视频| 久久人人爽人人爽人人片va | 欧美bdsm另类| 高清在线国产一区| 波多野结衣高清作品| 亚洲欧美日韩高清在线视频| 精品久久久久久久久av| 91九色精品人成在线观看| 午夜亚洲福利在线播放| 黄色日韩在线| 久久午夜亚洲精品久久| 男人舔奶头视频| 一级黄片播放器| 国产av在哪里看| 男人狂女人下面高潮的视频| 日本 av在线| 日韩欧美在线二视频| 香蕉av资源在线| 成人美女网站在线观看视频| 国产精品三级大全| 久久久久久久久久成人| 欧美成人a在线观看| 久久香蕉精品热| 亚洲av免费高清在线观看| 1000部很黄的大片| 97超级碰碰碰精品色视频在线观看| 亚洲欧美日韩无卡精品| 国内久久婷婷六月综合欲色啪| 黄片小视频在线播放| 国产在线精品亚洲第一网站| 久久久久久大精品| 日本熟妇午夜| 日韩欧美精品免费久久 | 国产亚洲精品久久久com| 婷婷精品国产亚洲av| 最近视频中文字幕2019在线8| 精品久久久久久久久久免费视频| 国产真实伦视频高清在线观看 | 亚洲最大成人手机在线| 亚洲最大成人中文| av在线天堂中文字幕| 日本一本二区三区精品| 热99在线观看视频| 欧美3d第一页| 亚洲 欧美 日韩 在线 免费| 男女视频在线观看网站免费| 国产精品亚洲美女久久久| a级毛片免费高清观看在线播放| 亚洲五月婷婷丁香| 午夜精品在线福利| 一个人看的www免费观看视频| 国产伦人伦偷精品视频| 亚洲av中文字字幕乱码综合| 国产91精品成人一区二区三区| 欧美另类亚洲清纯唯美| 亚洲国产精品999在线| 亚洲人成伊人成综合网2020| 国产极品精品免费视频能看的| 久久久久国内视频| 波野结衣二区三区在线| 国产精品久久久久久人妻精品电影| 午夜久久久久精精品| 国产精品亚洲一级av第二区| 床上黄色一级片| 级片在线观看| 免费av毛片视频| 午夜久久久久精精品| 中文亚洲av片在线观看爽| 中国美女看黄片| 黄色一级大片看看| 久久婷婷人人爽人人干人人爱| 午夜老司机福利剧场| 亚洲精品一卡2卡三卡4卡5卡| 日本a在线网址| 欧美日韩瑟瑟在线播放| 岛国在线免费视频观看| 韩国av一区二区三区四区| 97人妻精品一区二区三区麻豆| 极品教师在线视频| 在线观看一区二区三区| 久久99热6这里只有精品| 精品久久久久久,| 欧美日韩亚洲国产一区二区在线观看| 久久婷婷人人爽人人干人人爱| 精品久久久久久久久久久久久| 99久久精品热视频| 亚洲第一电影网av| 欧美午夜高清在线| 真实男女啪啪啪动态图| 久久精品国产清高在天天线| 国产大屁股一区二区在线视频| 久久久久久久久大av| 99热这里只有是精品50| 国产高清视频在线观看网站| 欧美最黄视频在线播放免费| 十八禁网站免费在线| 国产aⅴ精品一区二区三区波| 精品午夜福利在线看| 欧美潮喷喷水| 国产三级黄色录像| 亚洲国产日韩欧美精品在线观看| 淫秽高清视频在线观看| 韩国av一区二区三区四区| 午夜亚洲福利在线播放| 国产精品一区二区三区四区免费观看 | 欧美日韩亚洲国产一区二区在线观看| 欧美在线黄色| 亚洲男人的天堂狠狠| 国产白丝娇喘喷水9色精品| 他把我摸到了高潮在线观看| 国产久久久一区二区三区| 一区福利在线观看| 极品教师在线视频| 久久亚洲精品不卡| 级片在线观看| 亚洲av成人av| 国产精品电影一区二区三区| 90打野战视频偷拍视频| 搡老熟女国产l中国老女人| 制服丝袜大香蕉在线| 日韩精品中文字幕看吧| 亚洲精华国产精华精| 黄色日韩在线| 国产亚洲精品久久久com| 熟女电影av网| 久久久色成人| 成人毛片a级毛片在线播放| 午夜激情欧美在线| 免费看a级黄色片| 如何舔出高潮| 91麻豆精品激情在线观看国产| 看免费av毛片| av在线天堂中文字幕| 国产毛片a区久久久久| 日韩av在线大香蕉| 在线观看66精品国产| 波多野结衣高清作品| 久99久视频精品免费| 亚洲人成网站在线播放欧美日韩| 午夜福利在线观看吧| 九九在线视频观看精品| 久久久久九九精品影院| 亚洲 欧美 日韩 在线 免费| 成年人黄色毛片网站| 两个人的视频大全免费| 亚洲人成电影免费在线| 精品一区二区三区人妻视频| 欧美色视频一区免费| 日本在线视频免费播放| 日韩精品青青久久久久久| 欧美日韩福利视频一区二区| 成人毛片a级毛片在线播放| 成人av在线播放网站| 久久久久久九九精品二区国产| 简卡轻食公司| 91午夜精品亚洲一区二区三区 | 91久久精品国产一区二区成人| 亚洲欧美日韩东京热| 久久久久久久精品吃奶| 网址你懂的国产日韩在线| 欧美成人a在线观看| 国产精品1区2区在线观看.| 欧美成人a在线观看| 精品不卡国产一区二区三区| 国产成人影院久久av| 精品不卡国产一区二区三区| 亚洲最大成人手机在线| 不卡一级毛片| 亚洲人成电影免费在线| 亚洲专区中文字幕在线| 免费在线观看成人毛片| 婷婷精品国产亚洲av在线| 夜夜夜夜夜久久久久| 国产一区二区在线观看日韩| 欧洲精品卡2卡3卡4卡5卡区| 丰满的人妻完整版| 真实男女啪啪啪动态图| 麻豆国产97在线/欧美| 91狼人影院| 狠狠狠狠99中文字幕| 757午夜福利合集在线观看| 日本 欧美在线| 欧美激情久久久久久爽电影| 欧美一区二区亚洲| 欧美精品啪啪一区二区三区| 51午夜福利影视在线观看| 国内精品一区二区在线观看| 国内毛片毛片毛片毛片毛片| 一个人看视频在线观看www免费| www日本黄色视频网| 国内久久婷婷六月综合欲色啪| 亚洲av美国av| 真实男女啪啪啪动态图| 麻豆国产97在线/欧美| 国产亚洲精品综合一区在线观看| 性色av乱码一区二区三区2| 亚洲欧美日韩东京热| www.www免费av| 身体一侧抽搐| 欧美色欧美亚洲另类二区| 久久精品国产自在天天线| 久久久久久久亚洲中文字幕 | 一本精品99久久精品77| 亚洲国产欧洲综合997久久,| av专区在线播放| 床上黄色一级片| 国产成+人综合+亚洲专区| 男女之事视频高清在线观看| 国产精品乱码一区二三区的特点| 欧美成人a在线观看| 国产免费一级a男人的天堂| 非洲黑人性xxxx精品又粗又长| 日韩欧美精品v在线| www.熟女人妻精品国产| 亚洲熟妇中文字幕五十中出| 国产亚洲精品av在线| 久久6这里有精品| 国产精品1区2区在线观看.| 国产精品嫩草影院av在线观看 | 亚洲专区国产一区二区| 亚洲av免费高清在线观看| 亚洲乱码一区二区免费版| 日本三级黄在线观看|