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

    基于薄片有限元-無限元耦合模型的地鐵列車振動環(huán)境影響分析

    2017-08-30 12:22:29馬龍祥劉維寧蔣雅君晏啟祥
    振動與沖擊 2017年15期
    關(guān)鍵詞:有限元振動模型

    馬龍祥, 劉維寧, 蔣雅君, 晏啟祥

    (1.西南交通大學(xué) 土木工程學(xué)院,成都 610031;2. 西南交通大學(xué) 交通隧道工程教育部重點實驗室,成都 610031;3.北京交通大學(xué) 土木建筑工程學(xué)院,北京 100044)

    基于薄片有限元-無限元耦合模型的地鐵列車振動環(huán)境影響分析

    馬龍祥1,2, 劉維寧3, 蔣雅君1,2, 晏啟祥1,2

    (1.西南交通大學(xué) 土木工程學(xué)院,成都 610031;2. 西南交通大學(xué) 交通隧道工程教育部重點實驗室,成都 610031;3.北京交通大學(xué) 土木建筑工程學(xué)院,北京 100044)

    鑒于地鐵運營引發(fā)隧道及地層沿線路縱向的振動水平及特性不甚明確,應(yīng)用薄片有限元-無限元耦合模型對地鐵列車運營引發(fā)隧道基底、隧道壁及地表在垂直于線路的水平方向、鉛垂方向及線路縱向三個方向上的振動響應(yīng)進行高精度分析及比較。分析結(jié)果表明:①地鐵列車運營引發(fā)橫向、垂向及縱向三個方向6 Hz以下的低頻振動在地表均衰減得極其緩慢,且地表三個方向的振動速度、加速度響應(yīng)具有相似的頻率成分;②列車運營引起隧道壁的縱向振動響應(yīng)在1~100 Hz頻段內(nèi)較小,在該頻段內(nèi)的大部分頻率點處,其響應(yīng)甚至小于地表的縱向振動響應(yīng);③在隧道基底及隧道壁,由列車運營引發(fā)的縱向振動響應(yīng)在1~100 Hz頻段以內(nèi)顯著小于由其引發(fā)的垂向振動響應(yīng),但在地表,由列車運營引發(fā)的縱向振動響應(yīng)具有同橫向、垂向響應(yīng)相當(dāng)?shù)牧恐怠?/p>

    地鐵振動;環(huán)境響應(yīng);縱向振動響應(yīng);數(shù)值模擬;薄片有限元-無限元耦合模型

    近年來,地下鐵道作為緩解交通壓力的有效措施,在我國各大城市得到了迅猛發(fā)展。與此同時,由其運營對周圍環(huán)境產(chǎn)生的振動影響問題也引起了人們的廣泛關(guān)注[1-2]。一些學(xué)者也對該問題從分析預(yù)測方法[3-6]、現(xiàn)場振動測試[7-8]、減隔振措施性能[9]等方面或角度進行了研究。從相關(guān)研究中不難發(fā)現(xiàn),目前的絕大多數(shù)研究均只針對或關(guān)注鉛垂向及垂直于地鐵線路的水平向振動。然而,由于運行列車產(chǎn)生的振動激擾沿線路縱向分布并不均勻(由地鐵列車行進效應(yīng)、軌道不平順及離散軌道支撐等因素造成),其在引發(fā)隧道及地層產(chǎn)生上述兩個方向上振動的同時,也會引發(fā)相應(yīng)的沿線路縱向的振動。但目前針對地鐵運營引發(fā)隧道及地層縱向振動或同時關(guān)注線路三個方向上振動的研究還極為匱乏,學(xué)術(shù)及工程界對是否應(yīng)重視沿線路縱向的環(huán)境振動這一問題也尚未形成統(tǒng)一認(rèn)識。王文斌[10]通過對隧道內(nèi)自動落錘激勵引發(fā)的地表振動響應(yīng)進行研究,發(fā)現(xiàn)隧道內(nèi)的垂向激勵也會引發(fā)地表產(chǎn)生較大的縱向振動響應(yīng)。Sheng等[11]建立理論模型對垂向移動荷載作用于軌道進而引起周圍環(huán)境的振動響應(yīng)進行分析,得到了軌道上垂向移動荷載引發(fā)的振動響應(yīng)在線路三個方向上的分量處于同一數(shù)量級這一結(jié)論。王田友[12]通過現(xiàn)場實測及分析,得到了上海地鐵1號線運營引起線路近距離范圍內(nèi)(50 m范圍以內(nèi))地表縱向振動的量值和頻率成分。但由于并未涵蓋隧道及距離線路更遠地表點的縱向振動響應(yīng),他的工作也并不能全面反映地鐵縱向的環(huán)境振動水平及特性。綜上可見,地鐵運營引發(fā)隧道及地層沿線路縱向的振動水平及特性,特別是縱向振動相較于垂向及垂直于線路的水平向振動存有多大的量值及頻譜特性差異,在目前還不甚明確,尚需進行進一步研究。

    筆者于文獻[13]提出了地鐵列車運營引起環(huán)境振動預(yù)測分析的薄片有限元-無限元耦合模型,其在振動模擬中真實地考慮了列車的行進效應(yīng)、軌道不平順及離散軌道支撐等諸多因素,且妥善處理了動力計算的邊界問題,可對地鐵列車運營引起各方向上的環(huán)境振動進行有效分析。鑒于此,本文應(yīng)用薄片有限元-無限元耦合模型,對地鐵列車運營引發(fā)隧道基底、隧道壁及地表在鉛垂方向、垂直于線路的水平方向及線路縱向上的振動速度及加速度響應(yīng)進行計算及比較,全面分析地鐵列車引發(fā)環(huán)境振動在三個方向上的量值及衰減規(guī)律,并重點對沿線路縱向的環(huán)境振動水平及特性予以明確。本文的研究成果對是否應(yīng)重視地鐵運營引發(fā)線路縱向環(huán)境振動這一實際問題具有一定的指導(dǎo)意義。

    1 薄片有限元-無限元耦合模型

    本節(jié)對薄片有限元-無限元耦合模型進行簡要介紹,其詳細建立過程可參見文獻[13]。

    薄片有限元-無限元耦合模型是一個分析軌道交通振動環(huán)境影響問題的半解析半數(shù)值模型,其在頻域內(nèi)進行運算分析并首先給出頻域內(nèi)的振動響應(yīng),進而再由逆傅里葉變換求解時域內(nèi)的振動響應(yīng)。該模型考慮勻速運行列車產(chǎn)生的垂向振動激勵,使用基于無限-周期結(jié)構(gòu)理論的車軌動力耦合模型[14-15]對振源車軌系統(tǒng)的垂向振動進行分析,并以此計算出頻域內(nèi)列車經(jīng)由軌道傳遞給隧道基底的系列垂向振動激勵力(見圖1),進而再以這些振動激勵力為輸入對隧道-地層系統(tǒng)進行動力加載。在視隧道-地層結(jié)構(gòu)為沿線路縱向無限延伸的周期性結(jié)構(gòu)的基礎(chǔ)上,該模型通過利用隧道承受的系列振動激勵力在頻域內(nèi)所具有的特定周期性及疊加原理,可將線路縱向無限個振動激勵力作用引發(fā)的真實三維隧道-地層空間振動響應(yīng)問題的求解,歸結(jié)于研究線路縱向一個具有單元網(wǎng)格劃分且僅受一個振動激勵力作用的隧道-地層薄片范圍,見圖2。在橫斷面上,該薄片范圍由圍繞隧道結(jié)構(gòu)的有限單元和遠離振源近場、可提供良好邊界條件的無限單元所組成;在線路縱向,該薄片由于縱向長度很小(為一個扣件間距),只需用一個單元進行描述。即薄片有限元-無限元耦合模型只需對一個薄片范圍的隧道-地層系統(tǒng)進行計算分析,就能得到實際縱向(圖2中z方向)無限長的隧道-地層空間在系列振動激勵力作用下的振動響應(yīng)。顯然,由于大大減少了建模的自由度,該模型應(yīng)用于地鐵列車振動環(huán)境影響問題的分析,具有極高的計算效率及極強的計算能力。更為重要的是,該模型真實地考慮了列車行進效應(yīng)、軌道不平順及軌道離散支撐等因素的影響(通過輸入的振動激勵力加以考慮),規(guī)避了線路縱向的計算邊界(計算結(jié)果實際為線路縱向無限長的結(jié)果),且使用無限單元在線路橫斷面上提供了良好的動力計算邊界,因此,由其分析得到的振動響應(yīng),特別是縱向振動響應(yīng),將具有較大的參考價值。

    圖1 列車運營產(chǎn)生的振動激勵力

    圖2 列車振動環(huán)境影響分析的薄片有限元-無限元耦合模型

    Fig.2 Sliced finite element-infinite element coupling model for the analysis of train-induced environmental vibration

    2 依托工程及分析模型的建立

    以北京地鐵1號線東單站至建國門站區(qū)間隧道段為地鐵列車振動環(huán)境影響分析的依托工程。

    東單站至建國門站區(qū)間是馬蹄形隧道,其拱頂埋深為11.5 m,凈空寬4.9 m,高5.05 m,襯砌型式為復(fù)合式襯砌,初襯厚0.25 m,二襯厚0.3 m。初襯混凝土的彈性模量為28 500 MPa,泊松比0.2,密度2 450 kg/m3,材料阻尼比為0.02。二襯混凝土的彈性模量為30 000 MPa,泊松比0.2,密度2 500 kg/m3,材料阻尼比為0.02。根據(jù)地質(zhì)鉆孔資料,該區(qū)間地層可近似劃分為三層:第一層為雜填土,平均厚度為2.5 m;第二層土主要為粉質(zhì)黏土,平均厚度為18 m;第三層土為卵石和圓礫,其厚度在計算中可取無窮大。各層土的物理性質(zhì)見表1。地層及隧道橫斷面如圖3(a)所示。該地鐵區(qū)間軌道為普通DTVI2扣件軌道,運營列車為6節(jié)編組地鐵B型車,設(shè)計運營速度為60 km/h。根據(jù)上述工程條件及對稱性,取一半結(jié)構(gòu)建立薄片有限元-無限元耦合模型。模型縱向(線路延伸方向)長度按扣件間距取為0.6 m,橫斷面上有限元區(qū)域大小取深×寬為100 m×130 m,其中,隧道近場(隧道及其周邊土層)的有限元網(wǎng)格如圖3(b)所示。此外,建模時使計算區(qū)域內(nèi)有限單元尺寸不大于1 m。

    表1 東單站至建國門站區(qū)間地層參數(shù)表

    (a)橫斷面(b)隧道近場有限元網(wǎng)格

    圖3 區(qū)間橫斷面及隧道近場有限元網(wǎng)格

    Fig.3 Cross-sectional drawing of the subway section and finite element grid of the near field region of tunnel

    采用美國軌道五級譜來描述該區(qū)間的軌道不平順。應(yīng)用基于無限-周期結(jié)構(gòu)理論的車軌動力耦合模型計算隧道-地層結(jié)構(gòu)所受的頻域振動激勵力,并將其施加到建立的薄片有限元-無限元耦合模型上,以進行環(huán)境振動的分析。具體輸入到薄片有限元-無限元耦合模型中的環(huán)境振動激勵力的頻譜見圖4(該激勵力由不同頻率輪軌力誘發(fā)的多成分組成,而模型是按照輪軌力頻率依次輸入相應(yīng)成分并做相應(yīng)計算,以實現(xiàn)其的輸入,具體步驟參見文獻[13])。從圖4中可以看出環(huán)境振動激勵力的特性,即激勵力在10 Hz以下的低頻段具有較大量值,且在63 Hz附近存有峰值。此外,本文在環(huán)境振動分析時主要關(guān)注隧道結(jié)構(gòu)及地表的振動響應(yīng),具體取隧道基底、隧道壁及地表共5個振動觀測點,全面分析這5個點由地鐵列車通過所引發(fā)的鉛垂方向(后文簡寫為“垂向”)、垂直于線路的水平方向(后文簡寫為“橫向”)及線路縱向(后文簡寫為“縱向”)三個方向的速度及加速度響應(yīng)。所選取的5個振動觀測點分別為:距隧道中線1 m的隧道基底點、距基底頂面1.5 m的隧道壁點及距隧道中線40 m、80 m和100 m的3處地表點,如圖5所示。

    圖4 環(huán)境振動激勵力頻譜

    3 分析模型的準(zhǔn)確性驗證

    圖6給出了由上述模型計算得到的地鐵列車運行引起距隧道中線40 m處地表點的垂向加速度響應(yīng)模擬值與該處垂向加速度響應(yīng)實測值[6]的比較。由于相應(yīng)的現(xiàn)場實測工作只針對地表垂向振動開展了測試,因而此處僅給出垂向振動響應(yīng)模擬與實測值的對比。

    圖5 環(huán)境振動觀測點

    (a) 頻譜

    (b) 時程

    (c) 三分之一倍頻程譜

    從圖6中可以看到,模擬值與實測值無論是在頻譜、時程上,還是在1/3倍頻程譜上,均吻合得較好(模擬值和實測值在8~20 Hz頻段范圍內(nèi)之所以存有相對較大的差別,是因為測試地鐵振動時,鄰近道路上有少量的車輛通過,導(dǎo)致實測值相較于模擬值在該頻段范圍內(nèi)還多含有道路交通的貢獻),這充分驗證了上述建立的薄片有限元-無限元耦合模型在分析時的計算精度。

    4 地鐵列車振動環(huán)境影響分析

    圖7給出了地鐵列車運行引起基底、隧道壁及地表振動觀測點在橫向、垂向及縱向三個方向上的速度及加速度響應(yīng)三分之一倍頻程譜。

    對圖7中各觀測點三分之一倍頻程譜隨頻率的變化趨勢進行觀察,可以發(fā)現(xiàn):列車運營引發(fā)隧道及地表的速度及加速度響應(yīng)在63 Hz附近均會出現(xiàn)明顯峰值,這與車軌系統(tǒng)激發(fā)的(輸入薄片有限元-無限元耦合模型的)環(huán)境振動激勵力在該頻率處存在的峰值相對應(yīng)(見圖4);近激勵源的隧道壁、隧道基底的振動響應(yīng)三分之一倍頻程譜隨頻率的變化趨勢與遠激勵源的地表點有略微差異:對于隧道基底及隧道壁的三分之一倍頻程譜響應(yīng),隨著頻率的增大,在1~3 Hz頻段,響應(yīng)量值將有所降低或大致維持在相同水平,在3~63 Hz頻段,響應(yīng)量值在整體上將有所增大,而在63~100 Hz頻段,響應(yīng)量值在整體上又將出現(xiàn)下降;對于地表的三分之一倍頻程譜響應(yīng),隨著頻率的增大,在1~63 Hz頻段,響應(yīng)量值在整體上將有所增大,而同樣在63~100 Hz頻段,響應(yīng)量值在整體上又會有所降低。即,近激勵源的隧道基底、隧道壁觀測點與遠激勵源的地表點的速度、加速度振動響應(yīng)三分之一倍頻程譜在3~100 Hz頻段隨頻率變化的整體趨勢相近,但在1~3 Hz頻段隨頻率變化的趨勢有所不同。造成這種現(xiàn)象的原因為:近激勵源觀測點的振動響應(yīng)由準(zhǔn)靜態(tài)輪軌力與動態(tài)輪軌力共同決定,而遠激勵源觀測點的振動響應(yīng)受準(zhǔn)靜態(tài)輪軌力的影響很小,僅由動態(tài)輪軌力決定。加之,準(zhǔn)靜態(tài)輪軌力引起的振動響應(yīng)只集中在低頻段(此頻段與列車運行速度有關(guān),此處為3 Hz以下頻段),兩種效應(yīng)相疊加,就造成隧道基底、隧道壁在3 Hz以下頻段的振動響應(yīng)由準(zhǔn)靜態(tài)輪軌力與動態(tài)輪軌力共同決定,而它們(隧道基底、隧道壁)在3 Hz以上頻段及地表所有頻段的振動響應(yīng)都由動態(tài)輪軌力所決定,從而致使遠、近激勵源觀測點的響應(yīng)在三分之一倍頻程譜趨勢上出現(xiàn)上述差異。

    進一步觀察、比較圖7中各觀測點振動響應(yīng)的量值,并分析三個方向上環(huán)境振動的衰減規(guī)律,可以得到如下結(jié)論:

    (1) 對于隧道基底及隧道壁觀測點,橫向速度級(加速度級)在3 Hz以下低頻段可達78.6(36.4)dB, 在63 Hz時可達96.3(88.2)dB,垂向速度級(加速度級)在3 Hz以下低頻段可達86.9(49.4)dB, 在63 Hz時可達103(94.9)dB,縱向速度級(加速度級)在3 Hz以下低頻段可達75.2(39.6)dB, 在63 Hz時可達87.7(79.7)dB。

    (2) 對于地表觀測點,從距隧道中線40 m處到距隧道中線100 m處,橫向速度級(加速度級)最大量值分別為84.7(76.8)dB、75.6(67.6)dB、70.6(62.5)dB;垂向速度級(加速度級)最大量值分別為87.1(78.8)dB、74.5(64.8)dB、71.7(61.7)dB;縱向速度級(加速度級)最大量值分別為82.5(74.4)dB、74.6(65.4)dB、69.9(59.9)dB。

    (3) 對于橫向振動響應(yīng),由于隧道側(cè)壁臨空,其響應(yīng)在1~100 Hz以內(nèi)的大部分頻段最大,且地表在3~20 Hz頻段會出現(xiàn)橫向響應(yīng)大于基底的情況;對于垂向振動響應(yīng),在1~100 Hz以內(nèi)的大部分頻段,從基底到隧道壁,再到地表距隧道中線40 m、80 m及100 m處拾振點,振動響應(yīng)隨著拾振點距激勵源距離的增加,在整體上呈減小趨勢;對于縱向振動響應(yīng),隧道壁在1~100 Hz頻段以內(nèi)都將具有較小的量值,在該頻段內(nèi)的大部分頻率點處,其量值甚至小于地表的縱向響應(yīng)量值。而這種現(xiàn)象可能是由于隧道結(jié)構(gòu)的縱向變形剛度過大造成的。

    (4) 從隧道壁到地表距隧道中線40 m處,橫向振動響應(yīng)在3~6.5 Hz頻段幾乎沒有衰減,但在其他頻段會有較顯著的衰減;垂向振動響應(yīng)在5~8 Hz頻段幾乎沒有衰減,但在其他頻段會有較顯著的衰減;而縱向振動響應(yīng)在1.3~80 Hz頻段會有一定的放大。

    (5) 從地表距隧道中線40 m處到100 m處,橫向、垂向及縱向三個方向在6 Hz以上頻段的振動響應(yīng)均會有較為顯著的衰減,但三個方向在6 Hz以下的低頻段振動響應(yīng)卻均衰減得極其緩慢。

    (a) 速度響應(yīng)

    (b) 加速度響應(yīng)

    圖8給出了由地鐵列車運行引起隧道-地層系統(tǒng)不同觀測點橫向、垂向及縱向振動響應(yīng)三分之一倍頻程譜的比較。圖9給出了地表觀測點由地鐵列車運行引起的橫向、垂向及縱向振動響應(yīng)時程的比較。圖8及圖9中僅以加速度響應(yīng)為例進行比較,速度響應(yīng)的結(jié)果與此類似。

    從圖8及圖9中可以得到如下結(jié)論:

    (1) 無論是近激勵源的隧道基底、隧道壁拾振點,還是遠激勵源的地表拾振點,垂向、橫向及縱向三個方向上振動響應(yīng)的三分之一倍頻程譜隨頻率的變化趨勢在整體上是相似的;

    (2) 在隧道基底,縱向及橫向振動在1~100 Hz頻段以內(nèi)水平相當(dāng),均遠小于垂向振動水平;

    (3) 在隧道壁,橫向及垂向振動在1~100 Hz頻段以內(nèi)水平相當(dāng),均顯著大于縱向振動水平;

    (4) 對于地表垂向、橫向及縱向三個方向的振動,在5~30 Hz頻段,垂向最大,橫向次之,縱向最小,然而在包含振動最為顯著頻段的其他頻段,三個方向的振動水平相當(dāng),總體造成它們在時域量值上也僅具有較小的差別;且距線路距離越遠,地表三個方向振動水平的差別越小。

    圖10給出了典型的由地鐵列車運行引起的地表三個方向的振動速度及加速度響應(yīng)頻譜。圖中僅以距隧道中線80 m處地表觀測點的響應(yīng)為例,其它觀測點的響應(yīng)與此類似。

    (a) 隧道基底

    (b) 隧道壁

    (c) 地表距隧道中線40 m處

    (d) 地表距隧道中線80 m處

    (e) 地表距隧道中線100 m處

    (a) 地表距隧道中線40 m處

    (b) 地表距隧道中線80 m處

    (c) 地表距隧道中線100 m處

    (a) 速度頻譜

    (b) 加速度頻譜

    從圖10中可以看到,地鐵普通DTVI2扣件軌道段由列車運行引起地表的速度、加速度響應(yīng)在橫向、垂向及縱向三個方向上具有相似的頻率成分,其中速度響應(yīng)主要分布在80 Hz以下頻段,而加速度響應(yīng)主要分布在30~90 Hz頻段。

    5 結(jié) 論

    本文應(yīng)用薄片有限元-無限元耦合模型對地鐵列車運行引發(fā)隧道-地層系統(tǒng)垂向、橫向及縱向三個方向的速度及加速度響應(yīng)進行了高精度的細致分析及比較,得到了如下結(jié)論:

    (1) 地鐵列車運營引發(fā)橫向、垂向及縱向三個方向6 Hz以下的低頻振動在地表均衰減得極其緩慢;

    (2) 由列車運營引起隧道壁的縱向振動響應(yīng)在1~100 Hz頻段內(nèi)較小,在該頻段內(nèi)的大部分頻率點處,其響應(yīng)甚至小于地表的縱向振動響應(yīng);

    (3) 在隧道基底及隧道壁,由列車運營引發(fā)的縱向振動響應(yīng)在1~100 Hz頻段以內(nèi)顯著小于由其引發(fā)的垂向振動響應(yīng),但在地表,由列車運營引發(fā)的縱向振動響應(yīng)具有同橫向、垂向響應(yīng)相當(dāng)?shù)牧恐担?/p>

    (4) 地鐵列車運營引發(fā)地表的速度、加速度響應(yīng)在橫向、垂向及縱向三個方向上具有相似的頻率成分。在地鐵普通DTVI2扣件軌道段,由列車運行引起地表的速度響應(yīng)主要分布在80 Hz以下頻段,而引起的加速度響應(yīng)主要分布在30~90 Hz頻段。

    由上可見,在進行地鐵振動環(huán)境影響及控制領(lǐng)域內(nèi)的研究或設(shè)計時,地表縱向振動同橫、垂向振動一樣,也理應(yīng)受到相關(guān)人員的重視。

    [1] 劉維寧,夏禾,郭文軍. 地鐵列車振動的環(huán)境響應(yīng)[J]. 巖石力學(xué)與工程學(xué)報,1996,15(增刊1): 586-593.

    LIU Weining, XIA He, GUO Wenjun. Study of vibration effects of underground trains on surrounding environments[J]. Chinese Journal of Rock Mechanics and Engineering, 1996, 15(Sup1): 586-593.

    [2] 馬蒙,劉維寧,丁德云,等. 地鐵列車振動對精密儀器影響的預(yù)測研究[J].振動與沖擊,2011,30(3):185-190.

    MA Meng, LIU Weining, DING Deyun, et al. Prediction of influence of metro trains induced vibrations on sensitive instruments[J]. Journal of Vibration and Shock, 2011,30(3):185-190.

    [3] 劉維寧, 馬蒙, 王文斌. 地鐵列車振動環(huán)境響應(yīng)預(yù)測方法[J]. 中國鐵道科學(xué), 2013, 34(4): 110-117.

    LIU Weining, MA Meng, WANG Wenbin. Prediction method for subway train-induced environmental vibration responses[J]. Journal of China Railway Science, 2013, 34(4): 110-117.

    [4] DEGRANDE G, CLOUTEAU D, OTHMAN R, et al. A numerical model for ground-borne vibrations from underground railway traffic based on a periodic finite element-boundary element formulation[J]. Journal of Sound and Vibration, 2006, 293(3/4/5): 645-666.

    [5] 吳宗臻,劉維寧,馬龍祥,等. 地鐵浮置式軌道引起地表振動響應(yīng)解析預(yù)測模型研究[J].振動與沖擊,2014,33(17):132-137.

    WU Zongzhen, LIU Weining, MA Longxiang, et al. Analytical prediction model of ground vibration response induced by metro floating-type track[J]. Journal of Vibration and Shock, 2014,33(17):132-137.

    [6] 劉衛(wèi)豐.地鐵列車運行引起的隧道及自由場動力響應(yīng)數(shù)值預(yù)測模型研究[D]. 北京: 北京交通大學(xué), 2009.

    [7] DEGRANDE G, SCHEVENELS M, CHATTERJEE P, et al. Vibrations due to a test train at variable speeds in a deep bored tunnel embedded in London clay[J]. Journal of Sound and Vibration, 2006, 293(3/4/5): 626-644.

    [8] 栗潤德,張鴻儒,劉維寧. 地鐵引起的地面振動及其對精密儀器的影響[J]. 巖石力學(xué)與工程學(xué)報,2008, 27(1):206-214.

    LI Runde,ZHANG Hongru,LIU Weining. Metro-induced ground vibrations and their impacts on precision instrument[J]. Chinese Journal of Rock Mechanics and Engineering, 2008, 27(1):206-214.

    [9] 孫曉靜. 地鐵列車振動對環(huán)境影響的預(yù)測研究及減振措施分析[D]. 北京: 北京交通大學(xué), 2008.

    [10] 王文斌. 基于脈沖實驗的地鐵環(huán)境振動響應(yīng)傳遞函數(shù)預(yù)測方法研究[D]. 北京: 北京交通大學(xué), 2011.

    [11] SHENG X, JONES C J C, PETYT M. Ground vibration generated by a load moving along a railway track[J]. Journal of Sound and Vibration, 1999, 228(1):129-156.

    [12] 王田友. 地鐵運行所致環(huán)境振動與建筑物隔振方法研究[D]. 上海: 同濟大學(xué), 2008.

    [13] 馬龍祥. 基于無限-周期結(jié)構(gòu)理論的車軌耦合及隧道-地層振動響應(yīng)分析模型研究[D].北京: 北京交通大學(xué), 2014.

    [14] 馬龍祥,劉維寧,劉衛(wèi)豐. 移動荷載作用下周期支撐軌道結(jié)構(gòu)振動研究[J].中國鐵道科學(xué),2013, 34(1): 1-7.

    MA Longxiang, LIU Weining, LIU Weifeng. Study on vibration of periodic supported track structure under moving loads[J]. China Railway Science, 2013, 34(1): 1-7.

    [15] 馬龍祥,劉維寧,吳宗臻. 軌道結(jié)構(gòu)上輪對相互影響系數(shù)的解析求法[J]. 中南大學(xué)學(xué)報(自然科學(xué)版), 2014, 45(5): 1635-1641.

    MA Longxiang, LIU Weining, WU Zongzhen. Analytical method for solving wheelsets’ interaction coefficient on track structure[J]. Journal of Central South University (Science and Technology), 2014, 45(5): 1635-1641.

    Metro train-induced vibration influences on surrounding environments based on sliced finite element-infinite element coupled model

    MA Longxiang1,2, LIU Weining3, JIANG Yajun1,2, YAN Qixiang1,2

    (1. School of Civil Engineering, Southwest Jiaotong University, Chengdu 610031,China; 2. MOE Key Lab of Transportation Tunnel Engineering, Southwest Jiaotong University, Chengdu 610031,China; 3. School of Civil and Architectural Engineering, Beijing Jiaotong University, Beijing 100044,China)

    Because level and characteristic of metro train-induced tunnel and ground vibrations in the longitudinal direction along a subway line are not very clear, comprehensive high accuracy analyses and comparisons for metro train-induced vibration responses of tunnel base, tunnel wall and ground surface in all three directions were performed with the sliced finite element-infinite element coupled model. The results showed that metro train-induced lower frequency vibrations below 6 Hz on ground surface in the transverse direction, the vertical one and the longitudinal one are attenuated very slowly, the frequency components of ground surface vibration velocity and acceleration responses in three directions are similar; train-induced longitudinal vibration responses of the tunnel wall are smaller in the frequency range of 1-100 Hz, at most frequencies in this range, the responses of the tunnel wall are even smaller than those of the ground surface in the longitudinal direction; on the tunnel base and the tunnel wall, the train-induced vibration responses in the longitudinal direction in the frequency range of 1-100 Hz are significantly smaller than those in the vertical direction, but on the ground surface, the train-induced vibration responses in the longitudinal direction have comparative values to those in the transverse direction and vertical one.

    metro train-induced vibration; surrounding environment response; longitudinal vibration response; numerical simulation; sliced finite element-infinite element coupled model

    國家自然科學(xué)基金(51608456;51278425);中央高?;究蒲袠I(yè)務(wù)費專項資金(2682017CX003)

    2016-01-08 修改稿收到日期:2016-06-24

    馬龍祥 男,博士,講師,1988年10月生

    TU43; X827

    A

    10.13465/j.cnki.jvs.2017.15.017

    猜你喜歡
    有限元振動模型
    一半模型
    振動的思考
    重要模型『一線三等角』
    振動與頻率
    重尾非線性自回歸模型自加權(quán)M-估計的漸近分布
    中立型Emden-Fowler微分方程的振動性
    3D打印中的模型分割與打包
    磨削淬硬殘余應(yīng)力的有限元分析
    UF6振動激發(fā)態(tài)分子的振動-振動馳豫
    計算物理(2014年2期)2014-03-11 17:01:44
    基于SolidWorks的吸嘴支撐臂有限元分析
    插逼视频在线观看| 69av精品久久久久久| 春色校园在线视频观看| 中文在线观看免费www的网站| 国产毛片a区久久久久| 日韩一本色道免费dvd| av免费在线看不卡| 欧美xxxx黑人xx丫x性爽| 女的被弄到高潮叫床怎么办| 国产精品亚洲美女久久久| 亚洲激情五月婷婷啪啪| 变态另类丝袜制服| 好男人在线观看高清免费视频| 级片在线观看| 高清午夜精品一区二区三区 | 最新中文字幕久久久久| 一级毛片电影观看 | 国产一区亚洲一区在线观看| 俺也久久电影网| 别揉我奶头 嗯啊视频| 搞女人的毛片| 国产 一区精品| or卡值多少钱| 在线免费观看不下载黄p国产| 国产色婷婷99| 一区二区三区高清视频在线| 久久久久久久午夜电影| 美女 人体艺术 gogo| 联通29元200g的流量卡| 国产免费一级a男人的天堂| 国产精品,欧美在线| 久久午夜福利片| 成年av动漫网址| 黄色配什么色好看| 精品一区二区免费观看| 亚洲精品影视一区二区三区av| 亚洲av不卡在线观看| 久久久久免费精品人妻一区二区| 日本色播在线视频| av在线播放精品| 久久草成人影院| 伦精品一区二区三区| 人人妻人人澡欧美一区二区| 一区二区三区四区激情视频 | 亚洲最大成人中文| 久久这里只有精品中国| 国内精品宾馆在线| 免费一级毛片在线播放高清视频| 91久久精品电影网| 中文字幕熟女人妻在线| 国产视频内射| 日日摸夜夜添夜夜添av毛片| 18禁黄网站禁片免费观看直播| 亚洲av一区综合| 欧美3d第一页| 夜夜看夜夜爽夜夜摸| 婷婷色综合大香蕉| 日本爱情动作片www.在线观看 | 国产私拍福利视频在线观看| 亚洲国产色片| 成人精品一区二区免费| 欧美日韩在线观看h| 亚洲av一区综合| 蜜桃亚洲精品一区二区三区| 三级经典国产精品| 国产私拍福利视频在线观看| 久99久视频精品免费| 麻豆一二三区av精品| 天天躁夜夜躁狠狠久久av| 欧美一区二区国产精品久久精品| 久久精品国产99精品国产亚洲性色| 国产成人一区二区在线| 精品人妻视频免费看| 美女xxoo啪啪120秒动态图| 麻豆精品久久久久久蜜桃| 久久亚洲国产成人精品v| 神马国产精品三级电影在线观看| 成人漫画全彩无遮挡| 日本黄色片子视频| 黄色日韩在线| 天堂√8在线中文| 亚洲欧美中文字幕日韩二区| 91狼人影院| 深夜a级毛片| 一夜夜www| 精品午夜福利视频在线观看一区| 日本成人三级电影网站| 国产黄a三级三级三级人| 能在线免费观看的黄片| 一个人看视频在线观看www免费| 亚洲精品日韩av片在线观看| 性色avwww在线观看| 丝袜喷水一区| 国内精品美女久久久久久| АⅤ资源中文在线天堂| 啦啦啦啦在线视频资源| 免费av不卡在线播放| 亚洲欧美成人综合另类久久久 | 天美传媒精品一区二区| 色尼玛亚洲综合影院| 久久精品91蜜桃| 久久精品影院6| 日韩 亚洲 欧美在线| 免费av观看视频| 蜜臀久久99精品久久宅男| 午夜亚洲福利在线播放| 国内久久婷婷六月综合欲色啪| 国产亚洲精品久久久久久毛片| 在线免费观看的www视频| 亚洲av免费高清在线观看| 美女被艹到高潮喷水动态| 久久欧美精品欧美久久欧美| 国产黄色小视频在线观看| 校园春色视频在线观看| 国产成人精品久久久久久| 色哟哟·www| 在线看三级毛片| 一个人看视频在线观看www免费| 在线观看一区二区三区| 最后的刺客免费高清国语| 亚洲av成人精品一区久久| 在线观看av片永久免费下载| 狠狠狠狠99中文字幕| 亚洲图色成人| 久久99一区二区三区| 一级毛片电影观看| av线在线观看网站| 亚洲av国产av综合av卡| 六月丁香七月| 国产av一区二区精品久久| 国产精品偷伦视频观看了| 亚洲欧美日韩卡通动漫| 国产精品麻豆人妻色哟哟久久| 美女大奶头黄色视频| 亚洲国产精品一区三区| 高清在线视频一区二区三区| 欧美精品人与动牲交sv欧美| 日日爽夜夜爽网站| 成年女人在线观看亚洲视频| 婷婷色综合大香蕉| 男人和女人高潮做爰伦理| tube8黄色片| 国产成人freesex在线| 亚洲一区二区三区欧美精品| 久久精品久久久久久噜噜老黄| 91久久精品国产一区二区成人| 成人二区视频| 亚洲在久久综合| 国产真实伦视频高清在线观看| 校园人妻丝袜中文字幕| 另类精品久久| 亚洲精华国产精华液的使用体验| 亚洲欧美一区二区三区黑人 | 下体分泌物呈黄色| 日韩电影二区| 国产免费一级a男人的天堂| 一级片'在线观看视频| 99久国产av精品国产电影| 亚洲精品一区蜜桃| 国产成人午夜福利电影在线观看| 国产成人freesex在线| 天天操日日干夜夜撸| 成人国产av品久久久| 亚洲精品第二区| 成人美女网站在线观看视频| 黑丝袜美女国产一区| 国产一区二区在线观看日韩| 嫩草影院新地址| 热re99久久精品国产66热6| 久久久精品免费免费高清| 亚洲电影在线观看av| 搡女人真爽免费视频火全软件| 啦啦啦中文免费视频观看日本| 久久精品久久久久久久性| 欧美高清成人免费视频www| 3wmmmm亚洲av在线观看| 久久久久国产精品人妻一区二区| 黄色配什么色好看| 男女国产视频网站| 老女人水多毛片| 亚洲精品自拍成人| 亚洲精品国产av成人精品| 国产午夜精品一二区理论片| 看非洲黑人一级黄片| 亚洲综合精品二区| 亚洲精品国产色婷婷电影| 久久精品久久久久久久性| 美女脱内裤让男人舔精品视频| 免费av中文字幕在线| 中文字幕制服av| 国产免费一级a男人的天堂| 校园人妻丝袜中文字幕| 国产淫语在线视频| 久久精品国产自在天天线| 亚洲精品日本国产第一区| 一个人看视频在线观看www免费| 一本—道久久a久久精品蜜桃钙片| 少妇熟女欧美另类| 夫妻性生交免费视频一级片| 青青草视频在线视频观看| 在线观看人妻少妇| 国产欧美日韩综合在线一区二区 | 看免费成人av毛片| 成人特级av手机在线观看| 高清毛片免费看| 22中文网久久字幕| 一本久久精品| 国产成人freesex在线| 亚洲国产欧美在线一区| 内地一区二区视频在线| 自拍偷自拍亚洲精品老妇| 能在线免费看毛片的网站| 国产精品不卡视频一区二区| 97超视频在线观看视频| 国产无遮挡羞羞视频在线观看| 中文字幕人妻丝袜制服| 黄色毛片三级朝国网站 | 久久97久久精品| 成人综合一区亚洲| av福利片在线| 在线观看www视频免费| 在线天堂最新版资源| 9色porny在线观看| 少妇精品久久久久久久| 99久久精品热视频| 国产精品国产三级国产专区5o| 欧美日韩综合久久久久久| av一本久久久久| 日日啪夜夜爽| 黄片无遮挡物在线观看| 精品一区二区三区视频在线| 国产熟女欧美一区二区| 黄色毛片三级朝国网站 | 黄色毛片三级朝国网站 | 久久韩国三级中文字幕| 人妻一区二区av| 日韩中文字幕视频在线看片| 亚洲精品一区蜜桃| 国产一区二区在线观看日韩| 啦啦啦啦在线视频资源| 人妻夜夜爽99麻豆av| 日日摸夜夜添夜夜爱| 日韩欧美一区视频在线观看 | 欧美变态另类bdsm刘玥| 国产成人精品婷婷| 高清毛片免费看| 自拍偷自拍亚洲精品老妇| 国产精品熟女久久久久浪| 成年人免费黄色播放视频 | 国产成人精品一,二区| 国产 一区精品| 国产精品久久久久久精品电影小说| 看十八女毛片水多多多| av国产精品久久久久影院| 两个人免费观看高清视频 | 久久婷婷青草| 亚洲国产精品专区欧美| 国产亚洲一区二区精品| 一级,二级,三级黄色视频| 精品国产国语对白av| 人人妻人人爽人人添夜夜欢视频 | 久久精品国产a三级三级三级| 久久99精品国语久久久| 狂野欧美白嫩少妇大欣赏| 精品国产一区二区三区久久久樱花| 综合色丁香网| 国产伦理片在线播放av一区| 高清av免费在线| 老熟女久久久| 国内揄拍国产精品人妻在线| 91精品国产九色| 日韩av在线免费看完整版不卡| 最后的刺客免费高清国语| 精品人妻熟女毛片av久久网站| 欧美国产精品一级二级三级 | 国产日韩一区二区三区精品不卡 | 亚洲久久久国产精品| 亚洲内射少妇av| 国产日韩欧美亚洲二区| 一区二区av电影网| 久久久久久人妻| 欧美精品一区二区免费开放| 人妻系列 视频| 纵有疾风起免费观看全集完整版| 在现免费观看毛片| 夜夜看夜夜爽夜夜摸| 久久韩国三级中文字幕| 精品酒店卫生间| 久久综合国产亚洲精品| 少妇的逼水好多| 国产男人的电影天堂91| 日本欧美视频一区| 日韩大片免费观看网站| 国产亚洲av片在线观看秒播厂| 婷婷色综合大香蕉| 99热这里只有是精品在线观看| 边亲边吃奶的免费视频| 曰老女人黄片| 成人毛片60女人毛片免费| 日韩欧美精品免费久久| 色网站视频免费| 日日摸夜夜添夜夜爱| 人人妻人人爽人人添夜夜欢视频 | 欧美一级a爱片免费观看看| 草草在线视频免费看| 国产精品蜜桃在线观看| 日韩中字成人| 欧美日韩亚洲高清精品| 最黄视频免费看| 天堂俺去俺来也www色官网| 纵有疾风起免费观看全集完整版| 热99国产精品久久久久久7| 99精国产麻豆久久婷婷| 嫩草影院入口| 国产精品一区二区在线观看99| 久久久久久久精品精品| 国产一区二区三区综合在线观看 | 九色成人免费人妻av| videos熟女内射| 成人国产麻豆网| 国产精品伦人一区二区| 亚洲欧美精品专区久久| 免费看av在线观看网站| 欧美日韩在线观看h| 人妻系列 视频| 男女边摸边吃奶| 一本一本综合久久| 国产免费福利视频在线观看| 在线 av 中文字幕| 日韩欧美 国产精品| 色婷婷久久久亚洲欧美| 永久免费av网站大全| 日本欧美国产在线视频| 纵有疾风起免费观看全集完整版| 男人和女人高潮做爰伦理| 中文字幕人妻丝袜制服| 一本一本综合久久| 亚洲第一区二区三区不卡| 一本久久精品| 寂寞人妻少妇视频99o| 亚洲va在线va天堂va国产| 国产男女超爽视频在线观看| 91久久精品国产一区二区成人| 亚洲欧美日韩东京热| 天美传媒精品一区二区| 欧美激情极品国产一区二区三区 | 日韩强制内射视频| 91精品一卡2卡3卡4卡| 男女国产视频网站| 久久av网站| 搡老乐熟女国产| 国产又色又爽无遮挡免| 亚洲熟女精品中文字幕| 丝袜在线中文字幕| 久久精品久久久久久噜噜老黄| 晚上一个人看的免费电影| 日本色播在线视频| 蜜桃久久精品国产亚洲av| 婷婷色综合大香蕉| 一级av片app| 欧美日韩视频精品一区| 国产91av在线免费观看| 亚洲伊人久久精品综合| 少妇精品久久久久久久| 久久久国产一区二区| 中文字幕久久专区| 精品一区在线观看国产| 尾随美女入室| av在线播放精品| 国产成人精品福利久久| 六月丁香七月| 欧美三级亚洲精品| 国产精品久久久久久精品电影小说| 亚洲精品aⅴ在线观看| 麻豆成人午夜福利视频| 91午夜精品亚洲一区二区三区| 狂野欧美白嫩少妇大欣赏| 日韩av在线免费看完整版不卡| 少妇人妻精品综合一区二区| 国产黄频视频在线观看| 香蕉精品网在线| 国产熟女午夜一区二区三区 | 久久精品国产亚洲av涩爱| 日本爱情动作片www.在线观看| 亚洲国产日韩一区二区| 久久韩国三级中文字幕| 你懂的网址亚洲精品在线观看| 夜夜骑夜夜射夜夜干| 日本免费在线观看一区| 成年人午夜在线观看视频| 亚洲av欧美aⅴ国产| 国产精品一区www在线观看| 91午夜精品亚洲一区二区三区| 秋霞在线观看毛片| 国产熟女午夜一区二区三区 | 波野结衣二区三区在线| 国产精品久久久久成人av| 在线看a的网站| 下体分泌物呈黄色| 只有这里有精品99| 久久久久久久大尺度免费视频| 国产一区二区在线观看日韩| 国产日韩欧美亚洲二区| 一边亲一边摸免费视频| 免费看不卡的av| 欧美性感艳星| 国产精品.久久久| 精品国产国语对白av| 全区人妻精品视频| a 毛片基地| 久久久午夜欧美精品| 热re99久久精品国产66热6| 精华霜和精华液先用哪个| 亚洲美女搞黄在线观看| 亚洲国产欧美在线一区| 蜜桃在线观看..| 日韩成人av中文字幕在线观看| av有码第一页| 午夜福利影视在线免费观看| 中文天堂在线官网| 韩国av在线不卡| 国产 一区精品| 九九爱精品视频在线观看| 女人精品久久久久毛片| 日韩电影二区| 9色porny在线观看| 亚洲av二区三区四区| 一本久久精品| 一级,二级,三级黄色视频| 天堂俺去俺来也www色官网| 3wmmmm亚洲av在线观看| 免费黄频网站在线观看国产| 欧美另类一区| 亚洲国产日韩一区二区| 午夜福利影视在线免费观看| 日韩人妻高清精品专区| 美女视频免费永久观看网站| 久久久久久伊人网av| 成人综合一区亚洲| 99九九线精品视频在线观看视频| 精品酒店卫生间| 国产在视频线精品| 伊人亚洲综合成人网| 青春草国产在线视频| 欧美bdsm另类| 国产日韩欧美在线精品| 国产日韩一区二区三区精品不卡 | 我的老师免费观看完整版| 国产中年淑女户外野战色| 欧美精品一区二区免费开放| 久久精品国产鲁丝片午夜精品| 精品午夜福利在线看| 赤兔流量卡办理| 久久99热这里只频精品6学生| av天堂中文字幕网| 水蜜桃什么品种好| 精品久久久久久久久亚洲| 十八禁高潮呻吟视频 | 久久精品国产亚洲av涩爱| 亚洲高清免费不卡视频| 91久久精品电影网| 久热久热在线精品观看| 色吧在线观看| 18+在线观看网站| 亚洲欧美日韩卡通动漫| 永久网站在线| 性色av一级| 搡女人真爽免费视频火全软件| 国产成人精品一,二区| av天堂久久9| 伦理电影大哥的女人| 我要看黄色一级片免费的| 欧美最新免费一区二区三区| 毛片一级片免费看久久久久| 交换朋友夫妻互换小说| 亚洲av综合色区一区| 亚洲伊人久久精品综合| 午夜视频国产福利| 免费观看无遮挡的男女| 日本黄大片高清| 天堂俺去俺来也www色官网| 最新中文字幕久久久久| 亚洲三级黄色毛片| 嫩草影院新地址| 国产av码专区亚洲av| 少妇 在线观看| 一区在线观看完整版| 国产91av在线免费观看| 2018国产大陆天天弄谢| 亚洲无线观看免费| 嫩草影院入口| 久久免费观看电影| 午夜精品国产一区二区电影| 久久久国产一区二区| 国产精品99久久久久久久久| 亚洲va在线va天堂va国产| 亚洲一区二区三区欧美精品| 亚洲国产成人一精品久久久| 亚洲国产精品专区欧美| 另类亚洲欧美激情| 久久韩国三级中文字幕| 伊人亚洲综合成人网| 欧美日韩亚洲高清精品| 午夜激情久久久久久久| 一个人看视频在线观看www免费| 亚洲综合精品二区| 免费人妻精品一区二区三区视频| 久久精品国产鲁丝片午夜精品| 男女边吃奶边做爰视频| 国产精品.久久久| 国产精品久久久久久精品电影小说| 国产 一区精品| 久久久久久久亚洲中文字幕| 免费大片18禁| 亚洲第一区二区三区不卡| 国产精品一二三区在线看| 亚洲综合色惰| a级片在线免费高清观看视频| 大片免费播放器 马上看| 精品人妻偷拍中文字幕| 午夜精品国产一区二区电影| 精品久久久久久久久av| 女性生殖器流出的白浆| 狂野欧美激情性bbbbbb| 搡女人真爽免费视频火全软件| 久久这里有精品视频免费| 一二三四中文在线观看免费高清| 男的添女的下面高潮视频| 国产欧美日韩一区二区三区在线 | 水蜜桃什么品种好| 一区二区av电影网| 久久久欧美国产精品| av国产久精品久网站免费入址| 日韩av不卡免费在线播放| 成年av动漫网址| 国国产精品蜜臀av免费| 久久久久久久久久久免费av| 亚洲欧美一区二区三区国产| 日韩制服骚丝袜av| 熟女av电影| 插阴视频在线观看视频| 亚洲欧美精品专区久久| 免费少妇av软件| 啦啦啦中文免费视频观看日本| 国产成人精品无人区| 久久久久久久国产电影| 欧美另类一区| 大陆偷拍与自拍| 少妇熟女欧美另类| 日韩大片免费观看网站| 久久精品国产亚洲av涩爱| 在线 av 中文字幕| 国产极品粉嫩免费观看在线 | 99re6热这里在线精品视频| 亚洲欧美清纯卡通| 2022亚洲国产成人精品| 一个人免费看片子| 国产黄片视频在线免费观看| 亚洲av欧美aⅴ国产| 九九爱精品视频在线观看| 伦理电影免费视频| 大陆偷拍与自拍| 国产69精品久久久久777片| 亚州av有码| 欧美日本中文国产一区发布| 国产精品久久久久久精品电影小说| 欧美日韩在线观看h| 黄色欧美视频在线观看| 卡戴珊不雅视频在线播放| 天天躁夜夜躁狠狠久久av| 伊人久久国产一区二区| 久久综合国产亚洲精品| 成年美女黄网站色视频大全免费 | 国产成人午夜福利电影在线观看| 国语对白做爰xxxⅹ性视频网站| 亚洲美女黄色视频免费看| 国产在线一区二区三区精| 久久影院123| 大话2 男鬼变身卡| 蜜桃久久精品国产亚洲av| 国产精品一区二区在线不卡| 中文字幕人妻丝袜制服| 亚洲国产成人一精品久久久| 欧美另类一区| 国产av码专区亚洲av| 国产69精品久久久久777片| 久久人人爽人人爽人人片va| 少妇高潮的动态图| 秋霞在线观看毛片| 爱豆传媒免费全集在线观看| 国产成人freesex在线| 最近2019中文字幕mv第一页| 日日撸夜夜添| 欧美日本中文国产一区发布| 最近的中文字幕免费完整| 色吧在线观看| 成年人免费黄色播放视频 | 99久久精品一区二区三区| 9色porny在线观看| 成人亚洲欧美一区二区av| 少妇人妻 视频| 久久人人爽人人片av| av在线观看视频网站免费| 成年av动漫网址| 午夜久久久在线观看| 五月开心婷婷网| 久久久久网色| a级毛片在线看网站| 精品99又大又爽又粗少妇毛片| 啦啦啦在线观看免费高清www| 亚洲国产精品999| 婷婷色综合www| 国产黄色免费在线视频| 亚洲精品国产av蜜桃| 精品午夜福利在线看|