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

    交通荷載對(duì)雙井子堡古城墻振動(dòng)影響的現(xiàn)場(chǎng)實(shí)測(cè)與數(shù)值分析

    2023-06-25 13:40:24喬雄劉文高駱維斌倪偉淋楊鑫黃錦聰劉錦龍
    振動(dòng)工程學(xué)報(bào) 2023年3期
    關(guān)鍵詞:數(shù)值模擬

    喬雄 劉文高 駱維斌 倪偉淋 楊鑫 黃錦聰 劉錦龍

    摘要: 由于古建筑結(jié)構(gòu)材料特殊、歷史悠久且文化內(nèi)涵豐富,其在交通振動(dòng)下的動(dòng)力響應(yīng)規(guī)律研究成為古建筑安全性研究的重要內(nèi)容。針對(duì)古城墻在道路車輛運(yùn)行影響下的振動(dòng)高程效應(yīng),以及結(jié)構(gòu)的安全性評(píng)估問題,基于多種工況開展振動(dòng)特性的現(xiàn)場(chǎng)測(cè)試和振動(dòng)高程響應(yīng)機(jī)制的理論分析,結(jié)合數(shù)值模擬方法,并根據(jù)相應(yīng)的安全判據(jù)和允許標(biāo)準(zhǔn)進(jìn)行了研究。結(jié)果表明:道路交通荷載影響下古城墻的振動(dòng)主頻率范圍為 5~25 Hz,城墻頂部振速大于底部振速,頂部的放大系數(shù)范圍為 1.08~2.21,振動(dòng)放大效應(yīng)主要受車輛載荷和行駛速度的影響;隨著模擬車輛載荷以及行駛速度的增加,城墻質(zhì)點(diǎn)峰值速度放大系數(shù) fPPV增大,振動(dòng)放大效應(yīng)明顯,fPPV處在 1.02~1.77 之間。通過對(duì) fPPV增量 ACI(Amplification Coefficient Increment)進(jìn)行分析,表明車輛載荷對(duì)振動(dòng)高程效應(yīng)的影響程度大于車輛速度;對(duì)模擬結(jié)果以及現(xiàn)場(chǎng)實(shí)測(cè)數(shù)據(jù)進(jìn)行擬合,發(fā)現(xiàn)在交通振動(dòng)激勵(lì)下城墻底部和頂部的振速存在一定的線性關(guān)系;根據(jù)現(xiàn)場(chǎng)實(shí)測(cè)數(shù)據(jù)與數(shù)值模擬結(jié)果,對(duì)古土城墻的安全性進(jìn)行了評(píng)估,提出了基于通行車輛交通管理的古土城墻防振、減振措施。

    關(guān)鍵詞: 古土城墻;交通振動(dòng);現(xiàn)場(chǎng)測(cè)試;數(shù)值模擬;高程放大效應(yīng)

    中圖分類號(hào): TU361 文獻(xiàn)標(biāo)志碼: A 文章編號(hào): 1004-4523(2023)03-0776-13

    DOI:10.16385/j.cnki.issn.1004-4523.2023.03.020

    引 言

    隨著中國(guó)交通系統(tǒng)全方位、多層次、立體化的發(fā)展,公路交通的建設(shè)使得交通網(wǎng)絡(luò)的密集度大幅度增加,但隨之而來的交通振動(dòng)問題也開始引起人們的關(guān)注。統(tǒng)計(jì)數(shù)據(jù)表明,除了工廠作業(yè)振動(dòng)和建筑工程施工振動(dòng)外,人們反應(yīng)最為強(qiáng)烈的振源即為道路交通振動(dòng)[1]。國(guó)內(nèi)外由于交通振動(dòng)對(duì)鄰近古建筑產(chǎn)生破壞的問題日益突出[2‐5],由于交通振動(dòng)傳播形式較為復(fù)雜,古建筑年代久遠(yuǎn),且結(jié)構(gòu)和材料特殊,使得如何在保障交通運(yùn)行的前提下,保證古建筑的結(jié)構(gòu)安全成為亟待解決的問題。Chikaaki[6]通過對(duì)交通振動(dòng)作用下高架鐵路橋梁的動(dòng)力響應(yīng)進(jìn)行分析,為交通振動(dòng)的研究發(fā)展奠定了基礎(chǔ)。國(guó)內(nèi)關(guān)于古建筑在交通振動(dòng)激勵(lì)下的振動(dòng)響應(yīng)研究始于 20世紀(jì) 80 年代末,對(duì)焦枝鐵路建設(shè)中的爆破振動(dòng),以及列車運(yùn)營(yíng)時(shí)產(chǎn)生的振動(dòng)動(dòng)力響應(yīng)等問題進(jìn)行了研究和分析[7]。為了能夠探究車致振動(dòng)對(duì)鄰近建筑結(jié)構(gòu)的影響,Lombaert 等[8]通過對(duì)現(xiàn)場(chǎng)實(shí)測(cè)數(shù)據(jù)與數(shù)值模擬結(jié)果進(jìn)行對(duì)比,研究了車致振動(dòng)在自由場(chǎng)中的影響因素,并對(duì)各種影響因素進(jìn)行了分析論證。Fang 等[9]通過現(xiàn)場(chǎng)實(shí)測(cè)的方法對(duì)道路車輛引起的振動(dòng)進(jìn)行了研究,探討了時(shí)間域?qū)囍抡駝?dòng)的影響規(guī)律。謝偉平等[10]通過對(duì)客運(yùn)站不同車輛、不同行駛速度以及不同行駛方式下樓板的振動(dòng)響應(yīng)進(jìn)行實(shí)測(cè),分析了以上各種因素對(duì)振動(dòng)的動(dòng)力響應(yīng)以及頻譜特性的影響,并對(duì)車輛的行駛速度閾值給出了相應(yīng)的建議。Hinzen[11]通過對(duì)地鐵振動(dòng)作用下科隆大教堂的動(dòng)力響應(yīng)進(jìn)行現(xiàn)場(chǎng)實(shí)測(cè)分析,對(duì)科隆大教堂在地鐵運(yùn)營(yíng)時(shí)、微幅振動(dòng)作用下的結(jié)構(gòu)安全進(jìn)行了評(píng)估。Vala?ková 等[12]通過波動(dòng)理論分析與數(shù)值模擬方法,探究了交通振動(dòng)影響下古建筑結(jié)構(gòu)的動(dòng)力響應(yīng),并對(duì)振動(dòng)環(huán)境下的結(jié)構(gòu)振動(dòng)進(jìn)行了預(yù)測(cè)和評(píng)估。Ma 等[13]通過現(xiàn)場(chǎng)實(shí)測(cè)與數(shù)值模擬的方法,對(duì)地鐵交通環(huán)境下鐘樓的動(dòng)力響應(yīng)進(jìn)行研究,提出了降低地表質(zhì)點(diǎn)振動(dòng)速度峰值的方法。

    當(dāng)結(jié)構(gòu)質(zhì)點(diǎn)自振頻率處于外部激勵(lì)荷載的振動(dòng)頻率范圍內(nèi)時(shí),就會(huì)出現(xiàn)上部結(jié)構(gòu)的振速大于下部振速的現(xiàn)象,振動(dòng)速度幅值在一定范圍內(nèi)會(huì)隨著高度的增加而表現(xiàn)出一定的放大效應(yīng)[14‐17]。李政等[18]通過建立軌道交通的整體模型,結(jié)合砌體結(jié)構(gòu)模態(tài)分析,發(fā)現(xiàn)列車荷載誘發(fā)了砌體結(jié)構(gòu)自振頻率的增大,從而產(chǎn)生了結(jié)構(gòu)振動(dòng)強(qiáng)度的放大效應(yīng)。鄭鑫等[19]采用傳遞矩陣法對(duì)交通振動(dòng)引起的結(jié)構(gòu)豎向位移衰減關(guān)系進(jìn)行了理論分析,得出在速度界面的波阻抗越大和覆蓋層阻尼比越小時(shí),結(jié)構(gòu)局部放大越明顯。李克飛等[20]基于現(xiàn)場(chǎng)實(shí)測(cè)數(shù)據(jù)的分析,發(fā)現(xiàn)在地鐵交通環(huán)境下古城墻動(dòng)力響應(yīng)在水平和豎直方向呈規(guī)律性變化,并以水平方向?yàn)橹?,且振?dòng)在一定范圍內(nèi)隨著高度的增加呈現(xiàn)出放大效應(yīng)。葛家琪等[21]通過現(xiàn)場(chǎng)實(shí)測(cè)數(shù)據(jù)與數(shù)值模擬相結(jié)合的方法,評(píng)價(jià)了虎門炮臺(tái)舊遺址在交通振動(dòng)環(huán)境下的安全性,并提出采用減振溝的方法來削弱交通振動(dòng)帶來的振動(dòng)效應(yīng)。樊興偉[22]對(duì)開封古城墻在交通振動(dòng)作用下的動(dòng)力響應(yīng)進(jìn)行實(shí)測(cè),發(fā)現(xiàn)在水平方向上各拾振點(diǎn)的振動(dòng)速度以及位移幅值均隨著高度的增加而增大。

    目前,基于交通振動(dòng)研究技術(shù)的發(fā)展,雖然對(duì)于交通振動(dòng)環(huán)境下的建筑結(jié)構(gòu),特別是古建筑結(jié)構(gòu)動(dòng)力響應(yīng)的研究較多,但是相關(guān)的研究仍存在一定的缺陷。例如,在研究古建筑結(jié)構(gòu)動(dòng)力響應(yīng)問題時(shí),只考慮了在既有車輛載荷和行駛速度影響下建筑結(jié)構(gòu)的安全性,只研究了結(jié)構(gòu)基礎(chǔ)部分和頂部的振動(dòng)速度變化規(guī)律,并未對(duì)結(jié)構(gòu)自下而上質(zhì)點(diǎn)的振動(dòng)高程效應(yīng)變化規(guī)律進(jìn)行分析。由于振動(dòng)傳播規(guī)律的復(fù)雜性以及結(jié)構(gòu)材料的多樣性,不同振源類型以及不同結(jié)構(gòu)屬性對(duì)于振動(dòng)的高程效應(yīng)會(huì)有很大的影響。

    通過對(duì)公路鄰近兩側(cè)的古城墻遺址在交通荷載作用下振動(dòng)特性的現(xiàn)場(chǎng)實(shí)測(cè),得到了不同車次運(yùn)行下城墻的振動(dòng)響應(yīng)情況,并對(duì)古城墻基礎(chǔ)質(zhì)點(diǎn)以及城墻頂部的振動(dòng)速率進(jìn)行測(cè)試與統(tǒng)計(jì)分析,以此探討鄰近道路車輛運(yùn)行對(duì)古城墻的振動(dòng)速度高程放大效應(yīng)的影響。利用數(shù)值模擬方法來探究不同載荷、不同速度工況下的城墻自下而上質(zhì)點(diǎn)振動(dòng)速度放大系數(shù)的變化規(guī)律,以此分析車輛載荷和車輛行駛速度對(duì)城墻質(zhì)點(diǎn)振動(dòng)速度的影響規(guī)律,為在公路交通荷載作用下,鄰近公路的古建筑振動(dòng)安全評(píng)估與防護(hù)提供理論依據(jù)和技術(shù)支撐。

    1 城墻概況與試驗(yàn)方案

    1. 1 城墻概況

    雙井子堡始建于明代,城堡有東、西兩門,坐北向南,目測(cè)城堡呈矩形,作為歷史上嘉峪關(guān)的軍事防御前哨,具有十分重要的戰(zhàn)略地位。城墻屬于市級(jí)文物保護(hù)單位,通過對(duì)古土城墻現(xiàn)狀進(jìn)行觀測(cè),確定北墻長(zhǎng) 240.0 m、西墻長(zhǎng) 123.1 m、東墻長(zhǎng) 97.5 m、南墻 長(zhǎng) 94.7 m,城墻寬度為 0.6~1.6 m,城墻高度為2.7~3.7 m。雙井子堡鄰近于 G312 國(guó)道,另有 G30(連霍高速)高速將雙井子堡分為了南北兩塊區(qū)域,北側(cè)距離老的 G312(上海‐霍爾果斯)國(guó)道較近,現(xiàn)場(chǎng)概況如圖 1 所示。

    1. 2 測(cè)點(diǎn)布置

    選擇車輛運(yùn)行振動(dòng)速度觀測(cè)系統(tǒng)時(shí),應(yīng)根據(jù)現(xiàn)場(chǎng)實(shí)際情況預(yù)估被測(cè)信號(hào)的幅值范圍和頻率分布范圍,選擇的觀測(cè)系統(tǒng)幅值范圍上限應(yīng)高于被測(cè)信號(hào)幅值上限的 20%,頻響范圍應(yīng)包含被測(cè)信號(hào)的頻率分布范圍,依據(jù)這個(gè)原則選擇的觀測(cè)系統(tǒng)能避免出現(xiàn)削波、平臺(tái)等情況。為此,選擇 TC‐4850 型測(cè)振儀、TCS‐B3 型三軸向振動(dòng)速度傳感器和計(jì)算機(jī)組成的觀測(cè)系統(tǒng)(如圖 2 所示)進(jìn)行振動(dòng)信號(hào)的采集。

    根據(jù)現(xiàn)場(chǎng)情況,共布設(shè) 3 條測(cè)線。首先在雙井子堡北區(qū)靠西側(cè)鄰近 G312 國(guó)道的城墻位置,選擇垂直于道路方向的一段城墻進(jìn)行振動(dòng)特性測(cè)試,在城墻頂部和底部分別布設(shè)測(cè)線一和測(cè)線二,每條測(cè)線上布設(shè) 3 個(gè)測(cè)點(diǎn),距離振源的水平距離分別為 2,6 和14 m。其次在雙井子堡南區(qū)靠近西側(cè)鄰近 G30 高速方向布設(shè)測(cè)線三,測(cè)線上布設(shè) 3 個(gè)測(cè)點(diǎn),其中一個(gè)測(cè)點(diǎn)位于城墻頂部,距離振源的水平距離分別為3 m,3 m 和 14 m,現(xiàn)場(chǎng)各測(cè)線布設(shè)情況如圖 3所示。

    2 試驗(yàn)數(shù)據(jù)分析與結(jié)構(gòu)安全評(píng)估

    為了避免由于現(xiàn)場(chǎng)測(cè)試環(huán)境的不確定性因素帶來的干擾,在進(jìn)行現(xiàn)場(chǎng)測(cè)試數(shù)據(jù)處理分析時(shí)進(jìn)行一定的濾波措施;為了避免引發(fā)二次災(zāi)害,在不影響測(cè)試結(jié)果準(zhǔn)確性的前提下,對(duì)重車的運(yùn)行路線及運(yùn)行時(shí)間進(jìn)行了適當(dāng)限制。通過測(cè)試結(jié)果分別分析在不同車次運(yùn)行的情況下,古城墻基礎(chǔ)質(zhì)點(diǎn)和城墻頂部振動(dòng)動(dòng)力響應(yīng)特性以及振速的傳播規(guī)律。

    根據(jù)測(cè)試儀器的特性(采樣時(shí)間為 0~20 s),截取典型測(cè)點(diǎn)振動(dòng)速度峰值附近 14 s 速度時(shí)程段,為了能夠在空間狀態(tài)下更直觀地觀測(cè)古城墻的振動(dòng)響應(yīng)情況,通過比較不同車次運(yùn)行下振動(dòng)速度峰值變化規(guī)律來探討古城墻的振動(dòng)動(dòng)力響應(yīng)情況以及對(duì)其結(jié)構(gòu)安全性進(jìn)行評(píng)估。

    2. 1 實(shí)測(cè)數(shù)據(jù)分析

    本文在現(xiàn)場(chǎng)測(cè)試時(shí)選取了具有代表性的車輛進(jìn)行測(cè)試,作為鄰近道路上行駛車輛的參數(shù)依據(jù)。依據(jù)現(xiàn)場(chǎng)測(cè)線布置情況,分別對(duì)每條測(cè)線進(jìn)行不同的行車工況振動(dòng)測(cè)試,各測(cè)試工況如表 1 所示(各車重量分別為吊車:57 t;卡車:25 t;貨車:8 t)。

    現(xiàn)場(chǎng)測(cè)試共得到 169 條振動(dòng)數(shù)據(jù),考慮到路面還有其他車輛振動(dòng)的影響,為了保證數(shù)據(jù)分析的便利性以及結(jié)論的準(zhǔn)確性,每個(gè)工況選取 4 個(gè)典型的振動(dòng)數(shù)據(jù),并對(duì)符合現(xiàn)場(chǎng)振動(dòng)情況的樣本數(shù)據(jù)進(jìn)行分析處理,城墻各測(cè)線、各測(cè)點(diǎn)以及各工況振動(dòng)速度峰值和振動(dòng)速度時(shí)程曲線如圖 4~6 所示(由于篇幅限制,只給出典型速度時(shí)程曲線,即各測(cè)線、各試驗(yàn)方案振動(dòng)速度最大值)。

    在 1 號(hào)測(cè)點(diǎn)到 2 號(hào)測(cè)點(diǎn)之間振動(dòng)幅值隨水平距離衰減很明顯,衰減梯度較大;在 2 號(hào)測(cè)點(diǎn)到 3 號(hào)測(cè)點(diǎn)之間衰減不明顯,變化相對(duì)較小,說明由于雙井子堡古城墻較厚,振動(dòng)速度衰減較快,結(jié)構(gòu)振動(dòng)響應(yīng)較明顯。測(cè)線一、二所獲得的振動(dòng)原始數(shù)據(jù)為雙井子堡北區(qū)靠西側(cè)鄰近 G312 國(guó)道的城墻位置振動(dòng)速度時(shí)程,通過對(duì)不同測(cè)線所選取的 4 組振動(dòng)速度時(shí)程曲線進(jìn)行峰值的提取,通過圖 4~6 可以發(fā)現(xiàn)在相同振源情況下,古城墻底部和城墻頂部振動(dòng)速度峰值的變化情況。結(jié)合實(shí)測(cè)數(shù)據(jù)分析表明,古城墻頂部的實(shí)測(cè)振動(dòng)速度峰值比城墻底部大,表明古城墻墻體對(duì)于振動(dòng)有一定的放大效應(yīng),振動(dòng)速度在城墻頂部和底部衰減規(guī)律基本一致。

    測(cè)線三所獲得的振動(dòng)數(shù)據(jù)為雙井子堡南區(qū)靠近西側(cè)垂直于 G30 高速的城墻位置振動(dòng)速度時(shí)程,通過對(duì)測(cè)線三實(shí)測(cè)數(shù)據(jù)進(jìn)行分析,進(jìn)一步證明了在相同振源情況下古城墻頂部的實(shí)測(cè)振動(dòng)峰值比城墻底部大。結(jié)合測(cè)線一和測(cè)線二對(duì)比可以發(fā)現(xiàn),車輛重量和行駛速度越大,所引起的質(zhì)點(diǎn)振動(dòng)速度越大,結(jié)構(gòu)振動(dòng)響應(yīng)越明顯,且當(dāng)連續(xù)過車時(shí),振動(dòng)幅值相對(duì)增加較為明顯。通過對(duì)比相同車輛不同時(shí)速以及相同時(shí)速不同車輛測(cè)試方案發(fā)現(xiàn),相比于車輛行駛速度對(duì)古城墻質(zhì)點(diǎn)振動(dòng)響應(yīng)的影響,車輛載荷對(duì)質(zhì)點(diǎn)振動(dòng)響應(yīng)的影響更為明顯。

    為了探究雙井子堡古城墻的振動(dòng)放大效應(yīng)規(guī)律,定義振動(dòng)速度峰值 PPV(Peak Particle Velocity)的放大系數(shù)為 fPPV(Factor of Peak Particle Velocity),即為城墻研究質(zhì)點(diǎn)的振動(dòng)速度峰值 PPV 與城墻底部振動(dòng)速度峰值 PPV 之比,則城墻頂部各測(cè)點(diǎn) fPPV如表2 所示 。 城墻頂部的放大比例均處在1.08~2.21 倍之間(城墻頂部比底部的振動(dòng)速率大 0.01~0.022 cm/s),說明在一定范圍內(nèi)隨著古城墻高度的增加,城墻結(jié)構(gòu)振動(dòng)響應(yīng)更加顯著。

    2. 2 結(jié)構(gòu)安全性評(píng)估

    不同于地震和爆破強(qiáng)烈而又短暫的振動(dòng),交通振動(dòng)屬于微幅振動(dòng),會(huì)持續(xù)對(duì)建筑結(jié)構(gòu)產(chǎn)生影響,鄰近道路的建筑物在長(zhǎng)期交通振動(dòng)的影響下,易產(chǎn)生墻體開裂、裂縫發(fā)展及不均勻沉降等現(xiàn)象。尤其對(duì)于鄰近道路的古建筑結(jié)構(gòu)而言,由于結(jié)構(gòu)材料特殊、歷史悠久且文化內(nèi)涵豐富,其在交通振動(dòng)影響下的安全性備受關(guān)注。

    國(guó)內(nèi)外相關(guān)標(biāo)準(zhǔn)及文獻(xiàn)[23‐26]都是將建筑結(jié)構(gòu)質(zhì)點(diǎn)振動(dòng)速度峰值作為疲勞損傷的控制標(biāo)準(zhǔn),若古建筑質(zhì)點(diǎn)振動(dòng)速度的峰值不超過相關(guān)規(guī)定控制振速的限值,則此類振動(dòng)對(duì)古建筑結(jié)構(gòu)不會(huì)產(chǎn)生直接性的破壞。

    雙井子堡為古土城墻結(jié)構(gòu),是嘉峪關(guān)市的保護(hù)文物。參考《古建筑防工業(yè)振動(dòng)技術(shù)規(guī)范》[24]中對(duì)古木、磚和石結(jié)構(gòu)振動(dòng)速度控制限值以及上述國(guó)內(nèi)外古建筑控制標(biāo)準(zhǔn)和文獻(xiàn),確定雙井子堡的容許振動(dòng)速度為 0.03 cm/s。

    通過分析各測(cè)線不同測(cè)試工況的現(xiàn)場(chǎng)實(shí)測(cè)數(shù)據(jù),結(jié)果表明由雙井子堡鄰近道路上行駛的大型車輛所產(chǎn)生的振動(dòng)速度已超過容許振動(dòng)速度峰值0.03 cm/s,對(duì)雙井子堡古城墻的安全性具有一定的影響,因此為了保證雙井子堡的完整以及城墻的穩(wěn)定,應(yīng)當(dāng)對(duì)行駛的車輛載荷以及車輛速度進(jìn)行控制。

    3 交通振動(dòng)高程響應(yīng)機(jī)制

    為了能夠更加直接地分析古城墻的振動(dòng)放大效應(yīng),將所研究的城墻均勻地分成 n 等份,將城墻主體視為質(zhì)點(diǎn) m1,所研究城墻放大效應(yīng)的某質(zhì)點(diǎn)為 mn,則可得到城墻在外荷載作用下的計(jì)算力學(xué)模型,如圖 7 所示。

    式(5)中,當(dāng)| K - θ 2 M |≠ 0 時(shí),可求得城墻的振幅 Y,將其代入式(4)即得到簡(jiǎn)諧荷載作用下城墻質(zhì)點(diǎn)任意時(shí)刻的位移。很明顯,當(dāng)簡(jiǎn)諧荷載的主頻率 θ 與城墻結(jié)構(gòu)質(zhì)點(diǎn)自振頻率 ωn 相同時(shí),就會(huì)出現(xiàn)振動(dòng)放大效應(yīng)。

    在實(shí)際計(jì)算過程中,采用圖 7 所示的力學(xué)模型計(jì)算結(jié)構(gòu)的動(dòng)力響應(yīng),理論研究模型尺寸較規(guī)整,研究質(zhì)點(diǎn)的自振頻率與主體的自振頻率基本趨于一致。雙井子堡古城墻在自然環(huán)境下出現(xiàn)了不同程度的損壞,城墻頂部的寬度明顯小于中間位置和底部位置,質(zhì)量及剛度相對(duì)較小,因此其自振頻率會(huì)相對(duì)較大。現(xiàn)場(chǎng)實(shí)測(cè)振動(dòng)主頻率主要處在 5~25 Hz 之間(如圖 8 所示,由于篇幅限制,僅給出了典型時(shí)刻幅頻曲線圖),當(dāng)所研究城墻結(jié)構(gòu)質(zhì)點(diǎn)自振頻率處于交通振動(dòng)主頻率范圍內(nèi)時(shí),就會(huì)出現(xiàn)上部結(jié)構(gòu)的振速大于下部振速的現(xiàn)象,城墻出現(xiàn)在振動(dòng)荷載激勵(lì)作用下的高程放大效應(yīng)。

    4 城墻動(dòng)力響應(yīng)的數(shù)值分析

    考慮到古建筑結(jié)構(gòu)的重要性,為了避免試驗(yàn)對(duì)城墻的進(jìn)一步破壞,并未進(jìn)行不同軸重車輛影響下城墻振動(dòng)響應(yīng)的多點(diǎn)現(xiàn)場(chǎng)試驗(yàn)。但通過現(xiàn)場(chǎng)實(shí)測(cè)數(shù)據(jù)分析發(fā)現(xiàn),城墻在交通振動(dòng)影響下的動(dòng)力響應(yīng)受車輛載荷和行駛速度的影響較為明顯。為分析不同車輛載荷和行駛速度對(duì)城墻振動(dòng)響應(yīng)的影響情況,依據(jù)《道路車輛外輪廓尺寸、軸荷及質(zhì)量限值》[30]中的規(guī)定,車輛接地載荷在 600~980 kPa 之間。但在實(shí)際行車過程中,重載車大部分有超載情況,為了能夠使得研究更符合實(shí)際,故將接地壓強(qiáng)取值調(diào)整為700~1500 kPa[31]。假設(shè)接觸壓力等于胎壓,選取胎壓分別為 800,900,100 和 1100 kPa 進(jìn)行模擬,確定在保證古城墻處于安全振動(dòng)控制標(biāo)準(zhǔn)范圍之內(nèi)的車輛行駛速度,探討古城墻在交通環(huán)境下的振動(dòng)放大效應(yīng)。

    4. 1 模型與參數(shù)

    車輪與路面的接觸實(shí)際上是較為復(fù)雜的相互作用,黃仰賢[32]認(rèn)為輪胎與地面的接觸面積由一個(gè)矩形和兩個(gè)半圓組成,為了計(jì)算方便,可將接觸面積等效成矩形,有限元模型中車輪荷載等效接觸面積為0.2 m×0.16 m,雙輪中心間距取 0.3 m。為了模擬車輛在行駛過程中對(duì)鄰近古建筑的振動(dòng)影響,采用ABAQUS USER SUBROUTINE 編寫 Fortran 程序DLOAD 進(jìn)行移動(dòng)荷載的輸入,模型采用 C3D8R 等參單元,模型的邊界條件為地基及路面結(jié)構(gòu)約束 x或 y 方向,基礎(chǔ)土體底部全部約束。三維數(shù)值模型如圖 9 所示。

    道路材料屬性為瀝青混合料路面,其中,瀝青混合料面層厚為 0.12 m,基層厚為 0.3 m,基礎(chǔ)土體厚為 4 m,研究城墻厚為 0.6 m,長(zhǎng)和高均為 3 m。通過查閱相關(guān)地基土體和瀝青路面參數(shù)的文獻(xiàn)[33‐35],并對(duì)現(xiàn)場(chǎng)土城墻取樣進(jìn)行室內(nèi)試驗(yàn),確定各材料的物理參數(shù)如表 3 所示。

    4. 2 振動(dòng)動(dòng)力響應(yīng)規(guī)律分析

    研究?jī)?nèi)容主要為行駛車輛的載荷以及速度對(duì)城墻振動(dòng)速度動(dòng)力響應(yīng)的影響,城墻高度和長(zhǎng)度均為3 m,從水平方向和豎直方向分別每隔 0.6 m 取一個(gè)研究質(zhì)點(diǎn)(豎直方向取 A1~A6,水平方向取 B1~B6,距離均為 0,0.6,1.2,1.8,2.4,3.0 m)分析交通振動(dòng)速度的變化規(guī)律。各研究質(zhì)點(diǎn)選取位置如圖10所示。

    4. 2. 1 交通振動(dòng)高程效應(yīng)的規(guī)律分析

    在車輛荷載作用下,不同時(shí)刻雙井子堡古城墻的振速云圖如圖 11 所示。

    隨著車輛荷載不斷向城墻駛近,城墻的振動(dòng)速度從底部向頂部不斷傳播 ,t=0.02 s 時(shí)(見圖 11(a)),振動(dòng)荷載開始作用到城墻上,此時(shí)振動(dòng)波速由城墻底部向頂部傳播,由于振動(dòng)速度較小,并未出現(xiàn)振動(dòng)速度的放大效應(yīng);t=0.4325 s 時(shí)(見圖 11(b)),車輛繼續(xù)向城墻附近行駛,此時(shí)可以看出城墻開始出現(xiàn)振動(dòng)放大效應(yīng),但振動(dòng)速度在豎直方向上呈現(xiàn)先減小后增大的趨勢(shì);t=0.5849 s 時(shí)(見圖 11(c)),古城墻在交通載荷的作用下出現(xiàn)完全的振動(dòng)速度放大效應(yīng),由城墻底部向頂部的振動(dòng)速度逐漸增大;t=0.826 s 時(shí)(見圖 11(d)),車輛行駛到離城墻最近位置,此時(shí)振動(dòng)速度達(dá)到峰值。

    不同高程 Z(從城墻底部到頂部)的振動(dòng)速度峰值放大系數(shù) fPPV隨著車輛行駛速度的變化結(jié)果如圖12 所示。

    通過對(duì)圖 12 所示的 4 組不同車輛載荷和車速情況下,不同高程上的 fPPV結(jié)果分析發(fā)現(xiàn),當(dāng)車輛荷載P=800 kPa 時(shí)(見圖 12(a)),從城墻底部到城墻頂部方向,fPPV先減小后增大,在 Z=0~1.8 m 位置處,fPPV 自 下 而 上 逐 漸 減 小 ,在 Z=1.8~3.0 m 位 置 處 ,fPPV 自下而上開始呈現(xiàn)出增大的趨勢(shì)(如圖 11 振動(dòng)速度云圖所示,自下而上出現(xiàn)先減小后增大的動(dòng)力響應(yīng)狀態(tài)),且隨著車輛行駛速度從 V=20 km/h 到V=50 km/h 不斷增大,fPPV逐漸增大,但整體放大系數(shù)均小于 1,說明城墻頂部振動(dòng)速度并沒有超過城墻底部;當(dāng)車輛載荷為 P=900 kPa 時(shí)(見圖 12(b)),fPPV隨高度的增加逐漸增大,當(dāng)速度為 V=40 km/h和 V=50 km/h 時(shí),分別在 Z=2.4 m 和 Z=1.8 m 位置處 fPPV 大于 1,說明隨著車輛行駛速度的增加,振動(dòng)放大效應(yīng)逐漸變得明顯 ;當(dāng)車輛載荷為 P=1000 kPa 時(shí)(見圖 12(c)),fPPV基本大于 1,且 fPPV曲線上升趨勢(shì)逐漸趨于穩(wěn)定,說明隨著車輛荷載和速度的不斷增加,放大系數(shù)逐漸增大;當(dāng)車輛載荷 P=1100 kPa 時(shí)(見圖 12(d)),fPPV均大于 1,說明在移動(dòng)載荷速度為 V=20 km/h 到 V=50 km/h 的整個(gè)歷程中,古城墻都表現(xiàn)出速度放大效應(yīng)。隨著模擬車輛載荷以及行駛速度的增加,城墻質(zhì)點(diǎn)振動(dòng)速度放大系數(shù) fPPV 隨之增大,振動(dòng)放大效應(yīng)明顯。由圖 12(a)~(d)可知,fPPV 處在 1.02~1.77 之間,與現(xiàn)場(chǎng)實(shí)測(cè)數(shù)值基本一致。

    通過對(duì)比圖 12 中車輛載荷從 800 kPa 到 1100kPa 以及車速?gòu)?20 km/h 到 50 km/h 時(shí)的 fPPV的變化規(guī)律可以發(fā)現(xiàn),在車輛荷載不變的情況下,隨著車輛行駛速度的不斷增加,城墻質(zhì)點(diǎn) fPPV逐漸增大,但當(dāng)車輛載荷和行駛速度較小時(shí),質(zhì)點(diǎn)振動(dòng)速度會(huì)出現(xiàn)先減小后增大的現(xiàn)象,說明在交通振動(dòng)環(huán)境下城墻振動(dòng)放大效應(yīng)并不是一定存在,而是當(dāng)振動(dòng)強(qiáng)度達(dá)到一定值后才發(fā)生;在車輛速度不變的情況下,隨著車輛載荷的不斷增加,城墻質(zhì)點(diǎn) fPPV逐漸增大,且車輛載荷越大,fPPV增量也隨之越大。

    定義 ACI 為古土城墻頂部 3.0 m 位置與城墻0.6 m 位置的振動(dòng)速度放大系數(shù) fPPV之差。不同車速和不同車載工況下 ACI變化結(jié)果如圖 13 所示。

    通過對(duì)圖 13 所示的不同車輛載荷和不同車速情況下 ACI 的結(jié)果分析發(fā)現(xiàn),當(dāng)車輛載荷固定的情況下(見圖 13(a)),隨著車速的增加,各載荷下的ACI 表現(xiàn)出減小的趨勢(shì),只有當(dāng) P=1100 kPa 的情況下,ACI 出現(xiàn)持續(xù)增長(zhǎng)的趨勢(shì),說明隨著車輛載荷的增加,城墻振動(dòng)放大效應(yīng)變得更加明顯;當(dāng)車速固定的情況下(見圖 13(b)),隨著車輛載荷的不斷增加,各車速下的 ACI 表現(xiàn)出持續(xù)增長(zhǎng)的趨勢(shì),僅當(dāng)V=20 km/h 的情況下,ACI 出現(xiàn)了減小的情況,再次說明隨著車輛載荷的增加,城墻振動(dòng)放大效應(yīng)逐漸明顯。

    通過上述對(duì) ACI 的分析可以確定,對(duì)振動(dòng)高程效應(yīng)的顯著性影響程度為車輛載荷大于車輛速度。

    4. 2. 2 水平方向振動(dòng)衰減規(guī)律分析

    沿城墻水平方向,從鄰近道路城墻底部到城墻最遠(yuǎn)端位置,根據(jù)數(shù)值模型所選取的研究質(zhì)點(diǎn),振動(dòng)速度峰值隨著車輛行駛速度增加的變化結(jié)果如圖14 所示。

    通過對(duì)模擬數(shù)據(jù)進(jìn)行分析可知,各工況下的振動(dòng)速度峰值衰減規(guī)律大致相同,都隨著振源距離的增加,振動(dòng)速度峰值逐漸減小。當(dāng)車輛載荷為 1000kPa 和 1100 kPa 時(shí),距離振源最近位置質(zhì)點(diǎn)的振動(dòng)速度峰值均超過安全控制標(biāo)準(zhǔn) 0.03 cm/s。當(dāng)車輛載荷低于 900 kPa 時(shí),為了保證古土城墻結(jié)構(gòu)的安全性,車輛行駛速度應(yīng)控制在 40 km/h 以下;當(dāng)車輛載荷超過 900 kPa 時(shí) ,應(yīng)對(duì)車輛進(jìn)行提前改道或分流。

    5 城墻底部與頂部振動(dòng)速度峰值的關(guān)系

    在古建筑受交通振動(dòng)影響的實(shí)際監(jiān)測(cè)過程中,由于古建筑結(jié)構(gòu)的特殊性和受保護(hù)性,不便于對(duì)其頂部以及高程方向上質(zhì)點(diǎn)的振動(dòng)信號(hào)進(jìn)行采集,故在實(shí)際現(xiàn)場(chǎng)監(jiān)測(cè)中無法實(shí)現(xiàn)結(jié)構(gòu)本身振動(dòng)傳播規(guī)律的研究分析。根據(jù)現(xiàn)場(chǎng)監(jiān)測(cè)信號(hào)以及模擬結(jié)果分析,可以用地表質(zhì)點(diǎn)速度峰值作為表征值來反映城墻高程方向上的質(zhì)點(diǎn)振速。對(duì)不同振源類型情況下城墻頂部與城墻底部質(zhì)點(diǎn)振動(dòng)速度峰值的關(guān)系進(jìn)行研究分析,通過確定兩者之間的振速關(guān)系,以此利用城墻底部質(zhì)點(diǎn)的振動(dòng)速度來表現(xiàn)出城墻頂部的振動(dòng)速度。以模擬結(jié)果為研究?jī)?nèi)容,分別分析在車輛載荷為 800,900,1000 和 1100 kPa 情況下各研究質(zhì)點(diǎn)的振速關(guān)系。城墻底部與頂部的模擬振速關(guān)系如圖15 所示。

    圖 15 中 Vt為城墻頂部質(zhì)點(diǎn)振動(dòng)速度峰值;Vb為城墻底部質(zhì)點(diǎn)振動(dòng)速度峰值。

    通過圖 15 可以看出,擬合曲線相關(guān)性系數(shù) R2分別為 0.89,0.93,0.93 和 0.93。結(jié)合圖 12(a)發(fā)現(xiàn),當(dāng)載荷為 800 kPa 時(shí),城墻并沒有出現(xiàn)明顯的振動(dòng)放大現(xiàn)象,因此可以推斷出,隨著車輛載荷的增大,城墻在高程方向上的振動(dòng)放大效應(yīng)逐漸明顯,且在該放大條件下,城墻底部和頂部質(zhì)點(diǎn)振動(dòng)速度峰值的相關(guān)性較為顯著,擬合曲線相關(guān)性系數(shù)較大。

    同時(shí),為了進(jìn)一步驗(yàn)證城墻底部與頂部在鄰近道路交通振動(dòng)激勵(lì)下的振速線性關(guān)系,結(jié)合現(xiàn)場(chǎng)測(cè)試數(shù)據(jù),對(duì)現(xiàn)場(chǎng)采集的振動(dòng)信號(hào)進(jìn)行線性擬合。根據(jù)現(xiàn)場(chǎng)測(cè)試 12 種工況,選取測(cè)線一的 1‐1#,1‐2#,1‐3#(城墻頂部)和測(cè)線二的 2‐1#,2‐2#,2‐3#(城墻底部)在 4 種工況下各個(gè)測(cè)點(diǎn)的振動(dòng)信號(hào)以及測(cè)線三 3‐1#和 3‐2#在 4 種工況下城墻底部和頂部的振動(dòng)信號(hào)進(jìn)行數(shù)據(jù)分析,統(tǒng)計(jì)分析結(jié)果如圖 16 所示。

    通過圖 16 可以看出,現(xiàn)場(chǎng)實(shí)測(cè)城墻底部和頂部振動(dòng)信號(hào)之間同樣存在一定的線性關(guān)系,且線性相關(guān)系數(shù) R2=0.92,現(xiàn)場(chǎng)實(shí)測(cè)數(shù)據(jù)與數(shù)值模擬結(jié)果線性擬合相關(guān)系數(shù)基本一致,進(jìn)一步說明了模擬結(jié)果的準(zhǔn)確性以及可靠性;同時(shí)也確定了在交通振動(dòng)激勵(lì)下城墻底部和頂部的質(zhì)點(diǎn)振動(dòng)速度確實(shí)存在一定的線性比例關(guān)系。

    6 結(jié) 論

    (1)通過對(duì)現(xiàn)場(chǎng)實(shí)測(cè)數(shù)據(jù)進(jìn)行分析發(fā)現(xiàn),載重汽車運(yùn)行過程中,城墻的頂部振速大于底部振速,頂部放大系數(shù)的范圍為 1.08~2.21;振動(dòng)頻率的總體分布范圍為 0~40 Hz,集中在 5~25 Hz,振動(dòng)主要是以低頻振動(dòng)為主,高頻振動(dòng)衰減較為明顯。

    (2)當(dāng)結(jié)構(gòu)質(zhì)點(diǎn)自振頻率處于交通振動(dòng)主頻帶范圍內(nèi)時(shí),會(huì)出現(xiàn)振動(dòng)放大效應(yīng)。隨著模擬車輛載荷以及行駛速度的增加,城墻質(zhì)點(diǎn)振動(dòng)速度放大系數(shù) fPPV隨之增大,振動(dòng)放大效應(yīng)明顯,fPPV處在 1.02~1.77 之間,與現(xiàn)場(chǎng)實(shí)測(cè)數(shù)值基本一致。通過對(duì)不同載荷和行駛速度下的 ACI 進(jìn)行分析,發(fā)現(xiàn)對(duì)振動(dòng)高程效應(yīng)的顯著性影響程度為車輛載荷大于車輛速度。

    (3)通過數(shù)值模擬結(jié)果的分析,對(duì)不同車輛加載形式下基礎(chǔ)質(zhì)點(diǎn)與城墻頂部速度峰值進(jìn)行線性擬合,發(fā)現(xiàn)城墻底部和頂部的振動(dòng)速度擬合曲線相關(guān)系數(shù) R2的范圍為 0.89~0.93,而現(xiàn)場(chǎng)實(shí)測(cè)振動(dòng)信號(hào)擬合相關(guān)性系數(shù)為 0.92,表明鄰近道路的古土城墻在車致振動(dòng)作用下,其底部和頂部的振速存在一定的線性關(guān)系,實(shí)測(cè)與數(shù)值模擬結(jié)果可以相互印證。

    (4)實(shí)測(cè)與模擬結(jié)果均表明在車致振動(dòng)影響下,古土城墻的振動(dòng)速度超過了安全控制標(biāo)準(zhǔn)。為了確保古土城墻的安全性,依據(jù)振動(dòng)響應(yīng)規(guī)律,提出了相應(yīng)的減振、防振建議:對(duì)超過 900 kPa 的車輛進(jìn)行提前改道或分流;車輛荷載低于 900 kPa 的車輛應(yīng)控制其行駛速度在 40 km/h 以下。同時(shí),當(dāng)結(jié)構(gòu)產(chǎn)生振動(dòng)放大效應(yīng)時(shí),應(yīng)當(dāng)以結(jié)構(gòu)最高處的振動(dòng)速度作為允許振速的控制標(biāo)準(zhǔn)。

    參考文獻(xiàn):

    [1] 孫曉靜 . 地鐵列車振動(dòng)對(duì)環(huán)境影響的預(yù)測(cè)研究及減振措施分析[D]. 北京:北京交通大學(xué),2008.SUN Xiaojing. Prediction of environment vibrations in‐duced by metro trains and mitigation measures analysis[D]. Beijing: Beijing Jiaotong University,2008.

    [2] Bata Miloslav. Effects on buildings of vibrations causedby traffic[J]. Building Science,1971,6(4):221-246.

    [3] LIU X R, KOU M M, LU Y M, et al. An experimen‐tal investigation on the shear mechanism of fatigue dam‐age in rock joints under pre-peak cyclic loading condition[J]. International Journal of Fatigue, 2018, 106:175-184.

    [4] 朱利明,王成龍,藍(lán)天,等 . 地鐵運(yùn)行引起的南京鼓樓振動(dòng)測(cè)試與分析[J]. 建筑結(jié)構(gòu)學(xué)報(bào),2018,39(S1):291-296.

    ZHU Liming, WANG Chenglong, LAN Tian, et al.Vibration test and analysis of Nanjing Drum Towercaused by metro operation[J]. Journal of Building Struc‐tures,2018,39(S1):291-296.

    [5] 馬蒙,劉維寧 . 我國(guó)文物建筑受列車微振動(dòng)影響研究現(xiàn)況及關(guān)鍵問題分析[J]. 噪聲與振動(dòng)控制,2019,39(4):1-6.

    MA Meng, LIU Weining. Overview and key problemanalysis of the vibration influences on historic buildingsinduced by moving trains in China[J]. Noise and Vibra‐tion Control,2019,39(4):1-6.

    [6] Chikaaki Ueda. On distribution of traffic vibration in thebuilding under overhead railway bridges[J]. Transac‐tions of the Architectural Institute of Japan Summariesof Technical Papers,1965,40:375.

    [7] 王京禮,蘇廣豐 . 焦枝鐵路對(duì)洛陽(yáng)龍門石窟的環(huán)境影響[J]. 鐵道標(biāo)準(zhǔn)設(shè)計(jì)通訊,1991(12):35-39.

    WANG Jingli, SU Guangfeng. The environmental im‐pact of Jiaozhi Railway on Longmen Grottoes in Luoy‐ang[J]. Railway Standard Design Newsletter,1991(12):35-39.

    [8] Lombaert G, Degrande G. The experimental validationof a numerical model for the prediction of the vibrationsin the free field produced by road traffic[J]. Journal ofSound and Vibration,2003,262(2):309-331.

    [9] Fang Xinfeng, Li Dong. Field test and analysis for envi‐ronmental vibration induced by road traffic load[J]. Ap‐plied Mechanics and Materials,2013,482:183-187.

    [10] 謝偉平,張輝,何衛(wèi) . 多層立體汽車客運(yùn)站車致振動(dòng)實(shí)測(cè)研究[J]. 振動(dòng)與沖擊,2020,39(13):238-245.

    XIE Weiping, ZHANG Hui, HE Wei. Field measure‐ment for vehicle induced vibration of a multi-story 3-Dbus station[J]. Journal of Vibration and Shock,2020,39(13):238-245.

    [11] Hinzen K G. Subway-induced vibrations in cologne ca‐thedral[J]. Seismological Research Letters,2014,85(3):631-638.

    [12] Vala?ková Veronika, Papán Daniel, Fri?ová Jana. The‐oretical analysis of traffic seizmicity effect on importanthistorical building in modra town[J]. Procedia Earthand Planetary Science,2015,15(1):193-198.

    [13] Ma M, Liu W N, Qian C Y, et al. Study of the train-in‐duced vibration impact on a historic Bell Tower abovetwo spatially overlapping metro lines[J]. Soil Dynamicsand Earthquake Engineering,2016,81(3):58-74.

    [14] 孫鵬昌,盧文波,雷振,等 . 單薄山體巖質(zhì)高邊坡爆破振動(dòng)響應(yīng)分析及安全控制[J]. 巖土工程學(xué)報(bào),2021,43(5):877-885.

    SUN Pengchang, LU Wenbo, LEI Zhen, et al. Blast‐ing vibration response and control of high rock slopes ofthin mountain[J]. Chinese Journal of Geotechnical Engi‐neering,2021,43(5):877-885.

    [15] 吳多華,劉亞群,李海波,等 . 地震荷載作用下順層巖體邊坡動(dòng)力放大效應(yīng)和破壞機(jī)制的振動(dòng)臺(tái)試驗(yàn)研究[J]. 巖 石 力 學(xué) 與 工 程 學(xué) 報(bào) ,2020,39(10):1945-1956.

    WU Duohua, LIU Yaqun, LI Haibo, et al. Shaking ta‐ble tests on dynamic amplification and failure mecha‐nism of layered rock slopes under seismic actions[J].Chinese Journal of Rock Mechanics and Engineering,2020,39(10):1945-1956.

    [16] LI H B, LIU Y Q, LIU L B, et al. Numerical evalua‐tion of topographic effects on seismic response of singlefaced rock slopes[J]. Bulletin of Engineering Geologyand the Environment,2019,78(3):1873-1891.

    [17] YANG G X, QI S W, WU F Q, et al. Seismic amplifi‐cation of the anti-dip rock slope and deformation charac‐teristics: a large-scale shaking table test[J]. Soil Dy‐namics and Earthquake Engineering, 2018, 115:907-916.

    [18] 李政,金浩,鄭軍 . 軌道交通引起沿線砌體結(jié)構(gòu)振動(dòng)局 部 放 大 分 析[J]. 振 動(dòng) 與 沖 擊 ,2021,40(15):263-270.

    LI Zheng, JIN Hao, ZHENG Jun. Analysis of local am‐plification of masonry structure vibration caused by railtransit[J]. Journal of Vibration and Shock,2021,40(15):263-270.

    [19] 鄭鑫,陶夏新,王福彤,等 . 軌道交通地面振動(dòng)衰減關(guān)系 中 局 部 放 大 現(xiàn) 象 形 成 機(jī) 理 研 究[J]. 振 動(dòng) 與 沖 擊 ,2014,33(3):35-40.

    ZHENG Xin, TAO Xiaxin, WANG Futong, et al.Mechanism of local amplification in attenuation ofground vibration induced by rail traffic[J]. Journal of Vi‐bration and Shock,2014,33(3):35-40.

    [20] 李克飛,劉維寧,劉衛(wèi)豐,等 . 交通振動(dòng)對(duì)鄰近古建筑的動(dòng)力影響測(cè)試分析[J]. 北京交通大學(xué)學(xué)報(bào),2011,35(1):79-83.

    LI Kefei, LIU Weining, LIU Weifeng, et al. Tests andanalysis of traffic-induced vibration effects on surround‐ing historic buildings[J]. Journal of Beijing JiaotongUniversity,2011,35(1):79-83.

    [21] 葛家琪,馬伯濤,龐玉濤,等 . 汽車所致虎門遺址振動(dòng)全過程性能化研究[J]. 振動(dòng)與沖擊,2014,33(14):95-100.

    GE Jiaqi, MA Botao, PANG Yutao, et al. Whole pro‐cess performance-based analysis on vibration of Humensite induced by automobile running[J]. Journal of Vibra‐tion and Shock,2014,33(14):95-100.

    [22] 樊興偉 . 道路交通振動(dòng)對(duì)開封古城墻的影響研究[D].鄭州:河南大學(xué),2020.

    FAN Xingwei. Study on the influence of road traffic vi‐bration on Kaifeng ancient city wall[D]. Zhengzhou:Henan University,2020.

    [23] 中國(guó)機(jī)械工業(yè)聯(lián)合會(huì) . 建筑工程容許振動(dòng)標(biāo)準(zhǔn): GB50868—2013[S]. 北京:中國(guó)計(jì)劃出版社,2013.

    China Machinery Industry Federation. Standard for al‐lowable vibration of building engineering: GB 50868—2013[S]. Beijing: China Planning Press,2013.

    [24] 中華人民共和國(guó)住房和城鄉(xiāng)建設(shè)部 . 古建筑防工業(yè)振動(dòng)技術(shù)規(guī)范: GB/T 50452—2008[S]. 北京:中國(guó)建筑工業(yè)出版社,2008.

    Ministry of Housing and Urban-Rural Development ofthe Peoples Republic of China. Technical specifica‐tions for protection of historic buildings against manmade vibration: GB/T 50452—2008[S]. Beijing: Chi‐na Architecture & Building Press,2008.

    [25] Technical Committee Schwingungsfragen in Bauwesen.Vibrations in civil engineering: DIN 4150—1983[S].Berlin: Beuth Verlag GmbH,1983.

    [26] Hinzen K G. Subway-induced vibrations in Cologne Ca‐thedral[J]. Seismological Research Letters,2014,85(3):631-638.

    [27] 包世華 . 結(jié)構(gòu)動(dòng)力學(xué)[M]. 武漢:武漢理工大學(xué)出版社,2005:49-70.BAO Shihua. Structural Dynamics[M]. Wuhan: Wu‐han University of Technology Press,2005:49-70.

    [28] 王雷,宋維源,李吉 . 框架結(jié)構(gòu)鞭梢效應(yīng)的分析與討論[J]. 黑龍江科技學(xué)院學(xué)報(bào),2004,14(3):177-179.

    WANG Lei, SONG Weiyuan, LI Ji. Analysis and dis‐cussion of the whiplash effect of frames[J]. Journal ofHeilongjiang Institute of Science and Technology,2004,14(3):177-179.

    [29] 陳明,盧文波,李鵬,等 . 巖質(zhì)邊坡爆破振動(dòng)速度的高程放大效應(yīng)研究[J]. 巖石力學(xué)與工程學(xué)報(bào),2011,30(11):2189-2195.

    CHEN Ming, LU Wenbo, LI Peng, et al. Elevationamplification effect of blasting vibration velocity in rockslope[J]. Chinese Journal of Rock Mechanics and Engi‐neering,2011,30(11):2189-2195.

    [30] 中華人民共和國(guó)國(guó)家質(zhì)量監(jiān)督檢驗(yàn)檢疫總局 . 道路車輛外廓尺寸、軸荷及質(zhì)量限值: GB 1589—2004[S]. 北京:中國(guó)標(biāo)準(zhǔn)出版社,2004.

    General Administration of Quality Supervision, Inspec‐tion and Quarantine of the Peoples Republic of China.Limits of dimensions, axle load and masses of road vehi‐cles: GB 1589—2004[S]. Beijing: China StandardsPress,2004.

    [31] 董旭鵬 . 基于材料非線性的多軸重載下瀝青路面結(jié)構(gòu)響應(yīng)分析[D]. 南京:東南大學(xué),2015.

    DONG Xupeng. Analysis of asphalt pavement responesunder multiole wheel load based on material nonlinearity[D]. Nanjing: Southeast University,2015.

    [32] 黃仰賢 . 路面分析與設(shè)計(jì)[M]. 余定選,齊誠(chéng),譯 . 北京:人民交通出版社,1998:19-21.

    HUANG Yangxian. Pavement Analysis and Design[M]. Translated by Yu Dingxuan, Qi Cheng. Beijing:China Communications Press,1998:19-21.

    [33] 張杰 . 公路瀝青路面結(jié)構(gòu)動(dòng)力響應(yīng)及車輛舒適性分析[D]. 鄭州:鄭州大學(xué),2019.

    ZHANG Jie. Dynamic response of highway asphaltpavement structure and vehicle ride comfort analysis[D]. Zhengzhou: Zhengzhou University,2019.

    [34] 孔德政 . 河西地區(qū)長(zhǎng)城土遺址抗震穩(wěn)定性研究[D]. 蘭州:蘭州交通大學(xué),2017.

    KONG Dezheng. The study of seismic stability of therammed Great Wall in Hexi area[D]. Lanzhou: Lan‐zhou Jiaotong University,2017.

    [35] 張恒 . 甘肅地區(qū)夯土墻抗震性能試驗(yàn)研究[D]. 蘭州:蘭州交通大學(xué),2021.

    ZHANG Heng. Experimental study on seismic perfor‐mance of rammed earth walls in Gansu area[D]. Lan‐zhou: Lanzhou Jiaotong University,2021.

    猜你喜歡
    數(shù)值模擬
    基于AMI的雙色注射成型模擬分析
    錐齒輪精密冷擺輾成形在“材料成型數(shù)值模擬”課程教學(xué)中的應(yīng)用
    基于氣象信息及風(fēng)場(chǎng)信息的風(fēng)機(jī)輪轂處風(fēng)速預(yù)測(cè)
    鉆孔灌注樁樁底沉渣對(duì)樁體承載特性影響的模擬分析
    西南地區(qū)氣象資料測(cè)試、預(yù)處理和加工研究報(bào)告
    科技資訊(2016年18期)2016-11-15 08:01:18
    張家灣煤礦巷道無支護(hù)條件下位移的數(shù)值模擬
    科技視界(2016年18期)2016-11-03 23:14:27
    張家灣煤礦開切眼錨桿支護(hù)參數(shù)確定的數(shù)值模擬
    科技視界(2016年18期)2016-11-03 22:57:21
    跨音速飛行中機(jī)翼水汽凝結(jié)的數(shù)值模擬研究
    科技視界(2016年18期)2016-11-03 20:38:17
    姚橋煤礦采空區(qū)CO2防滅火的數(shù)值模擬分析
    雙螺桿膨脹機(jī)的流場(chǎng)數(shù)值模擬研究
    科技視界(2016年22期)2016-10-18 14:53:19
    欧美国产精品一级二级三级| 美女视频免费永久观看网站| 免费黄频网站在线观看国产| 女人久久www免费人成看片| 超色免费av| 国产成人欧美在线观看 | 国产男靠女视频免费网站| 国产精品免费一区二区三区在线 | 淫妇啪啪啪对白视频| 午夜福利影视在线免费观看| 亚洲伊人色综图| 日韩三级视频一区二区三区| 大香蕉久久成人网| 免费久久久久久久精品成人欧美视频| 欧美日韩精品网址| 久热这里只有精品99| 中文字幕另类日韩欧美亚洲嫩草| 成人国语在线视频| 欧美 亚洲 国产 日韩一| av线在线观看网站| 一个人免费在线观看的高清视频| 亚洲人成伊人成综合网2020| 在线观看舔阴道视频| ponron亚洲| 十八禁网站免费在线| 99riav亚洲国产免费| 欧洲精品卡2卡3卡4卡5卡区| av免费在线观看网站| 国产不卡一卡二| 亚洲九九香蕉| 国产97色在线日韩免费| 国产在视频线精品| 国产精品久久久av美女十八| 日本a在线网址| e午夜精品久久久久久久| 黄色毛片三级朝国网站| 亚洲国产欧美日韩在线播放| 亚洲精品国产精品久久久不卡| 老鸭窝网址在线观看| 亚洲av电影在线进入| 啦啦啦视频在线资源免费观看| 亚洲在线自拍视频| 丁香六月欧美| 在线国产一区二区在线| 黑人操中国人逼视频| 岛国在线观看网站| 国产视频一区二区在线看| 好男人电影高清在线观看| 在线视频色国产色| 国产免费av片在线观看野外av| 99re6热这里在线精品视频| 久久人妻福利社区极品人妻图片| 精品福利永久在线观看| a在线观看视频网站| 国产极品粉嫩免费观看在线| 欧美日韩亚洲国产一区二区在线观看 | 欧美日韩视频精品一区| 欧美大码av| 亚洲五月色婷婷综合| 精品视频人人做人人爽| 亚洲一区二区三区欧美精品| 多毛熟女@视频| 激情在线观看视频在线高清 | 久久热在线av| 悠悠久久av| 深夜精品福利| 99热国产这里只有精品6| 国精品久久久久久国模美| 丝袜在线中文字幕| 露出奶头的视频| 热99re8久久精品国产| 午夜日韩欧美国产| 午夜影院日韩av| 久久中文字幕人妻熟女| 亚洲精品国产一区二区精华液| 欧美激情久久久久久爽电影 | 亚洲全国av大片| 天堂√8在线中文| 香蕉国产在线看| 无人区码免费观看不卡| 王馨瑶露胸无遮挡在线观看| 老司机在亚洲福利影院| 午夜福利欧美成人| 欧美亚洲 丝袜 人妻 在线| 最近最新中文字幕大全电影3 | 黄色视频,在线免费观看| 99热国产这里只有精品6| 精品高清国产在线一区| 国产成人系列免费观看| 一边摸一边抽搐一进一小说 | 国产午夜精品久久久久久| 69精品国产乱码久久久| 欧美日韩国产mv在线观看视频| 18禁观看日本| 少妇裸体淫交视频免费看高清 | 日韩欧美在线二视频 | 黑丝袜美女国产一区| 亚洲性夜色夜夜综合| 久热爱精品视频在线9| 男人舔女人的私密视频| 99re在线观看精品视频| 黑丝袜美女国产一区| 成人黄色视频免费在线看| 免费久久久久久久精品成人欧美视频| 亚洲午夜理论影院| 不卡一级毛片| 久久中文看片网| 一级片免费观看大全| 操出白浆在线播放| 色婷婷av一区二区三区视频| 欧美日韩精品网址| 色综合欧美亚洲国产小说| 亚洲精品av麻豆狂野| www.自偷自拍.com| 极品教师在线免费播放| √禁漫天堂资源中文www| 亚洲中文字幕日韩| 久久国产精品大桥未久av| 中文字幕色久视频| 亚洲国产精品合色在线| 国产蜜桃级精品一区二区三区 | 手机成人av网站| 国产蜜桃级精品一区二区三区 | 亚洲精品av麻豆狂野| 国产亚洲欧美在线一区二区| 亚洲男人天堂网一区| 人妻 亚洲 视频| 亚洲,欧美精品.| 日本精品一区二区三区蜜桃| a级毛片在线看网站| 久久久国产成人精品二区 | 欧美乱妇无乱码| 亚洲免费av在线视频| 亚洲欧美日韩高清在线视频| 不卡av一区二区三区| 国产亚洲精品久久久久久毛片 | 欧美乱妇无乱码| 狂野欧美激情性xxxx| 欧美另类亚洲清纯唯美| 国产男女内射视频| 色综合婷婷激情| 亚洲熟女毛片儿| 精品视频人人做人人爽| 99精国产麻豆久久婷婷| 中文字幕人妻熟女乱码| 国产精品一区二区在线不卡| 午夜激情av网站| 好男人电影高清在线观看| 久久精品国产清高在天天线| 高清黄色对白视频在线免费看| aaaaa片日本免费| 免费观看精品视频网站| 久久久精品区二区三区| 国产淫语在线视频| 51午夜福利影视在线观看| 免费观看a级毛片全部| 国产97色在线日韩免费| 亚洲成a人片在线一区二区| 99久久精品国产亚洲精品| 高清在线国产一区| 一本大道久久a久久精品| 黄网站色视频无遮挡免费观看| 成年动漫av网址| 91成年电影在线观看| 国产精品久久视频播放| 亚洲第一欧美日韩一区二区三区| 免费一级毛片在线播放高清视频 | 国产精品秋霞免费鲁丝片| 人人妻人人添人人爽欧美一区卜| 亚洲欧美日韩另类电影网站| 人妻 亚洲 视频| 亚洲一卡2卡3卡4卡5卡精品中文| 精品福利观看| 亚洲精品粉嫩美女一区| 欧美中文综合在线视频| 悠悠久久av| 国产精品av久久久久免费| 欧美日韩亚洲国产一区二区在线观看 | 日韩 欧美 亚洲 中文字幕| 久久性视频一级片| 中文字幕精品免费在线观看视频| 国产精品国产av在线观看| 欧美色视频一区免费| 一本大道久久a久久精品| 91麻豆精品激情在线观看国产 | 啦啦啦在线免费观看视频4| 国产男靠女视频免费网站| 一级,二级,三级黄色视频| 国产99久久九九免费精品| 村上凉子中文字幕在线| 久久精品亚洲av国产电影网| 国产精品国产高清国产av | 欧美日韩国产mv在线观看视频| 日韩欧美一区视频在线观看| 免费在线观看亚洲国产| 亚洲欧美激情在线| 国产有黄有色有爽视频| 亚洲精品久久午夜乱码| 国产精品久久久av美女十八| 视频区图区小说| 欧美国产精品va在线观看不卡| 欧美黄色片欧美黄色片| 自拍欧美九色日韩亚洲蝌蚪91| 一区二区三区激情视频| 日韩欧美免费精品| 一级毛片精品| 他把我摸到了高潮在线观看| 狂野欧美激情性xxxx| 超碰97精品在线观看| 国产一区二区三区在线臀色熟女 | 欧美日韩精品网址| 国产亚洲欧美在线一区二区| 国产aⅴ精品一区二区三区波| 丝瓜视频免费看黄片| 亚洲黑人精品在线| 91在线观看av| 日本黄色日本黄色录像| 在线播放国产精品三级| 亚洲国产精品sss在线观看 | 欧美黄色淫秽网站| 女人被躁到高潮嗷嗷叫费观| 日韩人妻精品一区2区三区| 少妇猛男粗大的猛烈进出视频| 精品一品国产午夜福利视频| 欧美丝袜亚洲另类 | 日日爽夜夜爽网站| 美女午夜性视频免费| 亚洲一卡2卡3卡4卡5卡精品中文| 国产视频一区二区在线看| 亚洲九九香蕉| 99国产精品99久久久久| 欧美黑人欧美精品刺激| 欧美一级毛片孕妇| 女人被躁到高潮嗷嗷叫费观| 亚洲黑人精品在线| 大陆偷拍与自拍| 自拍欧美九色日韩亚洲蝌蚪91| 飞空精品影院首页| 在线观看免费日韩欧美大片| 手机成人av网站| 国产成人精品在线电影| 亚洲 欧美一区二区三区| 国产淫语在线视频| 久久精品国产99精品国产亚洲性色 | 久久久久久久久久久久大奶| 美国免费a级毛片| 精品乱码久久久久久99久播| 美女扒开内裤让男人捅视频| 无遮挡黄片免费观看| 后天国语完整版免费观看| 欧美日韩瑟瑟在线播放| 国产亚洲欧美在线一区二区| 人人澡人人妻人| 亚洲国产精品sss在线观看 | 欧美丝袜亚洲另类 | 咕卡用的链子| 中文字幕人妻丝袜一区二区| 亚洲精品国产一区二区精华液| 日韩制服丝袜自拍偷拍| 在线观看午夜福利视频| 又黄又粗又硬又大视频| 亚洲精品粉嫩美女一区| 人人妻人人添人人爽欧美一区卜| 少妇被粗大的猛进出69影院| 精品欧美一区二区三区在线| 村上凉子中文字幕在线| 精品久久久久久电影网| 超碰97精品在线观看| 免费人成视频x8x8入口观看| 不卡av一区二区三区| 国产三级黄色录像| 精品国内亚洲2022精品成人 | 精品一区二区三区四区五区乱码| 老熟妇仑乱视频hdxx| 满18在线观看网站| 美国免费a级毛片| 99在线人妻在线中文字幕 | 日本黄色日本黄色录像| 国产在线精品亚洲第一网站| 80岁老熟妇乱子伦牲交| 国产激情久久老熟女| av天堂在线播放| 国产精品成人在线| 久久国产精品影院| 亚洲黑人精品在线| 国产精品国产高清国产av | 国产精品一区二区精品视频观看| 免费女性裸体啪啪无遮挡网站| 久久国产亚洲av麻豆专区| 久久中文看片网| 免费久久久久久久精品成人欧美视频| 国产xxxxx性猛交| 99久久国产精品久久久| 成人av一区二区三区在线看| 国产色视频综合| 国产精品1区2区在线观看. | 亚洲午夜理论影院| 捣出白浆h1v1| 国产aⅴ精品一区二区三区波| 午夜福利欧美成人| 亚洲综合色网址| 乱人伦中国视频| 亚洲一卡2卡3卡4卡5卡精品中文| 婷婷成人精品国产| 午夜精品在线福利| 精品国产超薄肉色丝袜足j| 国产精品 国内视频| 国产免费男女视频| 久久中文字幕一级| 亚洲成人手机| 成人免费观看视频高清| 国精品久久久久久国模美| 免费在线观看完整版高清| 大片电影免费在线观看免费| 国产一区在线观看成人免费| 中文字幕色久视频| 一级a爱片免费观看的视频| 很黄的视频免费| 美女视频免费永久观看网站| 欧美色视频一区免费| 国产亚洲精品久久久久久毛片 | 每晚都被弄得嗷嗷叫到高潮| 国产欧美日韩一区二区精品| 悠悠久久av| 中文欧美无线码| 国产精品一区二区精品视频观看| 亚洲国产欧美一区二区综合| 老熟妇乱子伦视频在线观看| 亚洲av欧美aⅴ国产| 不卡av一区二区三区| 俄罗斯特黄特色一大片| 亚洲一区二区三区不卡视频| 精品国产乱码久久久久久男人| 在线国产一区二区在线| 国产男女超爽视频在线观看| 久久精品国产综合久久久| 黄色视频不卡| 日日摸夜夜添夜夜添小说| 午夜福利,免费看| 成人影院久久| 久久天堂一区二区三区四区| 国产aⅴ精品一区二区三区波| 国产免费av片在线观看野外av| 校园春色视频在线观看| 亚洲国产看品久久| 欧洲精品卡2卡3卡4卡5卡区| 国产免费av片在线观看野外av| 美女扒开内裤让男人捅视频| 人妻 亚洲 视频| 国产精品免费大片| 久热这里只有精品99| 校园春色视频在线观看| 国产高清激情床上av| av网站在线播放免费| 十八禁高潮呻吟视频| 很黄的视频免费| 国产深夜福利视频在线观看| 国产麻豆69| www.自偷自拍.com| 免费不卡黄色视频| 午夜福利视频在线观看免费| 中文字幕色久视频| 久久狼人影院| 国产精品免费大片| 欧美午夜高清在线| xxx96com| 欧美日本中文国产一区发布| 亚洲色图综合在线观看| 国产免费现黄频在线看| 亚洲国产欧美网| 女人久久www免费人成看片| 中文字幕av电影在线播放| www.999成人在线观看| av天堂久久9| 在线永久观看黄色视频| 亚洲欧洲精品一区二区精品久久久| 一区二区日韩欧美中文字幕| av一本久久久久| 午夜福利一区二区在线看| 午夜成年电影在线免费观看| 国产精品国产av在线观看| 久久久久精品国产欧美久久久| 久久久精品免费免费高清| 丰满人妻熟妇乱又伦精品不卡| 国产免费av片在线观看野外av| 高清在线国产一区| 深夜精品福利| 国精品久久久久久国模美| 欧美激情极品国产一区二区三区| 亚洲av片天天在线观看| 狠狠婷婷综合久久久久久88av| 搡老熟女国产l中国老女人| 免费高清在线观看日韩| 欧美 亚洲 国产 日韩一| 男女免费视频国产| 国产精品99久久99久久久不卡| 精品福利永久在线观看| 麻豆乱淫一区二区| 国产熟女午夜一区二区三区| 中文字幕制服av| 三级毛片av免费| 亚洲av成人不卡在线观看播放网| 国产区一区二久久| 免费少妇av软件| 精品午夜福利视频在线观看一区| 久久精品国产亚洲av香蕉五月 | 久久久国产精品麻豆| 亚洲av成人一区二区三| 叶爱在线成人免费视频播放| av不卡在线播放| 狠狠狠狠99中文字幕| 老熟妇乱子伦视频在线观看| 丝袜人妻中文字幕| 国产无遮挡羞羞视频在线观看| 两性夫妻黄色片| 一进一出抽搐gif免费好疼 | 搡老熟女国产l中国老女人| 少妇裸体淫交视频免费看高清 | 97人妻天天添夜夜摸| 亚洲精品自拍成人| 国产99白浆流出| 最新的欧美精品一区二区| 中文亚洲av片在线观看爽 | 国产野战对白在线观看| 精品午夜福利视频在线观看一区| 他把我摸到了高潮在线观看| 一级毛片女人18水好多| 久久国产精品大桥未久av| 一本大道久久a久久精品| 男女下面插进去视频免费观看| 一二三四在线观看免费中文在| 免费看十八禁软件| 亚洲国产欧美日韩在线播放| 日本撒尿小便嘘嘘汇集6| 天天影视国产精品| 精品欧美一区二区三区在线| 他把我摸到了高潮在线观看| 国产99久久九九免费精品| 90打野战视频偷拍视频| 欧美日韩瑟瑟在线播放| 桃红色精品国产亚洲av| 一级毛片女人18水好多| 女人久久www免费人成看片| 91九色精品人成在线观看| 久久人妻熟女aⅴ| 美女国产高潮福利片在线看| 91国产中文字幕| 12—13女人毛片做爰片一| 精品一区二区三区四区五区乱码| 天天影视国产精品| 久久影院123| 亚洲熟妇中文字幕五十中出 | 亚洲专区中文字幕在线| 久久精品国产99精品国产亚洲性色 | svipshipincom国产片| 99热只有精品国产| 99国产精品99久久久久| 最近最新中文字幕大全免费视频| 变态另类成人亚洲欧美熟女 | 满18在线观看网站| 成人特级黄色片久久久久久久| 后天国语完整版免费观看| av线在线观看网站| 欧美另类亚洲清纯唯美| 亚洲成av片中文字幕在线观看| 亚洲成人手机| 交换朋友夫妻互换小说| 亚洲精品中文字幕在线视频| 亚洲成a人片在线一区二区| 精品视频人人做人人爽| 亚洲在线自拍视频| 极品教师在线免费播放| 国产精品1区2区在线观看. | 黑人巨大精品欧美一区二区mp4| 麻豆国产av国片精品| 黑人巨大精品欧美一区二区蜜桃| 久久久国产一区二区| 亚洲成人国产一区在线观看| 亚洲中文日韩欧美视频| 亚洲精品中文字幕一二三四区| 亚洲欧美精品综合一区二区三区| 后天国语完整版免费观看| 日本五十路高清| 熟女少妇亚洲综合色aaa.| 亚洲视频免费观看视频| 成人精品一区二区免费| 国产精品久久久久久人妻精品电影| 国产视频一区二区在线看| 日本一区二区免费在线视频| 久久狼人影院| 脱女人内裤的视频| 51午夜福利影视在线观看| 麻豆乱淫一区二区| 很黄的视频免费| 亚洲精品久久成人aⅴ小说| 怎么达到女性高潮| 狠狠婷婷综合久久久久久88av| 黑人巨大精品欧美一区二区mp4| 精品熟女少妇八av免费久了| 国产xxxxx性猛交| 久久人人爽av亚洲精品天堂| 黄色视频,在线免费观看| 国产精品一区二区精品视频观看| 操出白浆在线播放| 国产精品香港三级国产av潘金莲| 在线播放国产精品三级| 老汉色av国产亚洲站长工具| 9191精品国产免费久久| 亚洲欧美一区二区三区黑人| 亚洲欧美精品综合一区二区三区| 后天国语完整版免费观看| 日韩制服丝袜自拍偷拍| 熟女少妇亚洲综合色aaa.| 1024香蕉在线观看| 啦啦啦在线免费观看视频4| 99精国产麻豆久久婷婷| 下体分泌物呈黄色| 久久九九热精品免费| 一进一出抽搐gif免费好疼 | 国产精品欧美亚洲77777| 亚洲精品一二三| 精品少妇久久久久久888优播| 女人久久www免费人成看片| 无限看片的www在线观看| 亚洲精华国产精华精| 亚洲精品国产区一区二| 日本撒尿小便嘘嘘汇集6| 一a级毛片在线观看| 国产aⅴ精品一区二区三区波| 亚洲精品国产精品久久久不卡| 国产亚洲精品一区二区www | 亚洲精品成人av观看孕妇| 国产极品粉嫩免费观看在线| 欧美精品高潮呻吟av久久| 国产精品永久免费网站| 啦啦啦 在线观看视频| 午夜福利影视在线免费观看| 久久天躁狠狠躁夜夜2o2o| 亚洲精品在线美女| 午夜老司机福利片| 99re在线观看精品视频| 怎么达到女性高潮| 日本a在线网址| 亚洲三区欧美一区| 国产成人一区二区三区免费视频网站| 99re在线观看精品视频| 久久人人爽av亚洲精品天堂| 国产男女超爽视频在线观看| 水蜜桃什么品种好| 午夜精品久久久久久毛片777| 丁香欧美五月| 人人澡人人妻人| 亚洲性夜色夜夜综合| 水蜜桃什么品种好| 99精品久久久久人妻精品| 日韩成人在线观看一区二区三区| 黄色丝袜av网址大全| 久久久久精品人妻al黑| 欧美精品啪啪一区二区三区| 最新的欧美精品一区二区| 嫩草影视91久久| 亚洲综合色网址| 99精品欧美一区二区三区四区| 丝袜美足系列| 一边摸一边抽搐一进一小说 | 99久久精品国产亚洲精品| 日韩欧美一区视频在线观看| 一级毛片精品| 亚洲av熟女| 999精品在线视频| 高清av免费在线| 久久精品人人爽人人爽视色| 国产视频一区二区在线看| 欧美黄色淫秽网站| 日本黄色日本黄色录像| 97人妻天天添夜夜摸| 极品人妻少妇av视频| 亚洲,欧美精品.| 国产精品久久久久成人av| 亚洲一卡2卡3卡4卡5卡精品中文| 亚洲精品久久成人aⅴ小说| 午夜日韩欧美国产| 亚洲五月婷婷丁香| 国产一区有黄有色的免费视频| 99精品在免费线老司机午夜| 亚洲五月天丁香| 涩涩av久久男人的天堂| 国产在线精品亚洲第一网站| 高清黄色对白视频在线免费看| 亚洲全国av大片| 中文字幕最新亚洲高清| 天堂中文最新版在线下载| ponron亚洲| 国产高清videossex| 精品国产一区二区三区久久久樱花| 动漫黄色视频在线观看| 成熟少妇高潮喷水视频| 人人妻人人澡人人爽人人夜夜| 麻豆乱淫一区二区| 777米奇影视久久| 国产精品九九99| 一二三四在线观看免费中文在| 国产97色在线日韩免费| 亚洲熟妇熟女久久| 777久久人妻少妇嫩草av网站| 久久这里只有精品19| 一本一本久久a久久精品综合妖精| 又黄又爽又免费观看的视频| 国产麻豆69| 日本精品一区二区三区蜜桃| 欧美日韩精品网址| 国产97色在线日韩免费| 啦啦啦 在线观看视频| 在线观看免费视频日本深夜| 精品欧美一区二区三区在线|