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

    基于小波能量譜的盾構(gòu)隧道損傷間接識(shí)別

    2022-12-15 01:14:40周劍鋒
    振動(dòng)與沖擊 2022年23期
    關(guān)鍵詞:小波粗糙度盾構(gòu)

    羅 輝, 周劍鋒,2

    (1. 華中科技大學(xué) 土木與水利工程學(xué)院, 武漢 430074; 2. 中冶南方工程技術(shù)有限公司, 武漢 430223)

    我國(guó)城市軌道交通建設(shè)發(fā)展迅速,以盾構(gòu)隧道形式為主的地鐵正成為城市軌道交通運(yùn)輸?shù)闹袌?jiān)力量,地鐵的安全運(yùn)營(yíng)對(duì)城市的運(yùn)轉(zhuǎn)至關(guān)重要。因此國(guó)內(nèi)城市地鐵相繼建立了盾構(gòu)隧道健康監(jiān)測(cè)系統(tǒng)并發(fā)展了相應(yīng)的盾構(gòu)隧道損傷診斷方法,但現(xiàn)有的盾構(gòu)隧道健康監(jiān)測(cè)面臨著傳感器需求量大、監(jiān)測(cè)點(diǎn)布設(shè)復(fù)雜、經(jīng)濟(jì)成本高、全覆蓋監(jiān)測(cè)困難等問題亟需解決。

    Yang等[1]在2004年提出了橋梁間接測(cè)量法的概念,不同于直接健康監(jiān)測(cè)方法,間接測(cè)量法通過車輛動(dòng)力響應(yīng)間接識(shí)別橋梁的模態(tài)和損傷信息,具有移動(dòng)性、經(jīng)濟(jì)性、便捷性等優(yōu)點(diǎn)。國(guó)內(nèi)外學(xué)者開展了基于間接測(cè)量法的損傷識(shí)別研究,Nguyen等[2-4]通過對(duì)車橋相互作用有限元模擬得到的車輛動(dòng)力響應(yīng)進(jìn)行小波分析,通過小波分析對(duì)突變信號(hào)的敏感性識(shí)別橋梁結(jié)構(gòu)的損傷。Mei等[5]通過智能手機(jī)采集車輛通過橋梁的動(dòng)力響應(yīng),提取振動(dòng)數(shù)據(jù)的MFCCs系數(shù)作為特征,利用Kullback-Leibler散度比較特征的分布,對(duì)同一橋梁損傷分布的差異進(jìn)行量化來識(shí)別損傷。王樹棟等[6]提出了由過橋汽車的動(dòng)力響應(yīng)識(shí)別橋梁損傷的靈敏度分析方法,采用最小二乘法和正則化方法求解靈敏度方程進(jìn)行參數(shù)更新和模型修正識(shí)別了橋梁的損傷。

    近年來小波分析方法在土木工程損傷識(shí)別領(lǐng)域有著廣泛的應(yīng)用。Zhu等[7]對(duì)移動(dòng)荷載作用下含裂縫梁結(jié)構(gòu)的位移響應(yīng)進(jìn)行連續(xù)小波變換,利用小波變換系數(shù)識(shí)別了裂縫位置。丁幼亮等[8]提出了利用結(jié)構(gòu)動(dòng)力響應(yīng)的小波包能量譜進(jìn)行結(jié)構(gòu)損傷預(yù)警,具有較好的損傷敏感性和噪聲魯棒性。項(xiàng)貽強(qiáng)等[9]對(duì)測(cè)點(diǎn)損傷前后的動(dòng)態(tài)信號(hào)進(jìn)行離散小波變換,計(jì)算各測(cè)點(diǎn)的小波總能量相對(duì)變化作為損傷識(shí)別指標(biāo),通過一個(gè)兩端固支梁數(shù)值算例驗(yàn)證了損傷指標(biāo)的有效性。戰(zhàn)家旺等[10]結(jié)合離散小波變換和信息熵,提出了一種基于小波能量熵的橋梁損傷預(yù)警方法,通過鋼桁梁橋數(shù)值算例驗(yàn)證了所提損傷預(yù)警方法的有效性。這些研究表明,在信號(hào)小波變換的基礎(chǔ)之上,通過小波能量構(gòu)造損傷指標(biāo)可以提高損傷識(shí)別的準(zhǔn)確率和靈敏性。

    基于以上分析,本文提出將橋梁間接測(cè)量法的思想應(yīng)用到盾構(gòu)隧道結(jié)構(gòu)健康監(jiān)測(cè)中,以移動(dòng)車輛和盾構(gòu)隧道耦合振動(dòng)作用下車輛加速度響應(yīng)信號(hào)作為損傷信息數(shù)據(jù),采用連續(xù)小波變換的信號(hào)處理方法,基于小波能量譜構(gòu)建盾構(gòu)隧道損傷間接識(shí)別指標(biāo)。建立了移動(dòng)車輛和盾構(gòu)隧道耦合振動(dòng)有限元模型,通過數(shù)值算例驗(yàn)證了本文所提盾構(gòu)隧道損傷間接識(shí)別方法的有效性,并分析了車速、測(cè)量噪聲、粗糙度等因素的影響。在實(shí)驗(yàn)室搭建了車輛和盾構(gòu)隧道耦合振動(dòng)模型,進(jìn)行了損傷識(shí)別試驗(yàn)以驗(yàn)證本文方法的有效性。

    1 小波能量譜損傷指標(biāo)

    1.1 時(shí)間-小波能量譜

    (1)

    此時(shí),稱ψ(t)為一個(gè)基本小波或母小波。將母函數(shù)ψ(t)經(jīng)伸縮和平移后得

    (2)

    稱其為一個(gè)小波序列,其中a為尺度因子,b為平移因子。

    對(duì)于任意函數(shù)x(t)∈L2(R)的連續(xù)小波變換為

    (3)

    式中,ψ*(t)是ψ(t)的共軛。

    因?yàn)閤(t)的小波變換能量守恒,所以

    (4)

    式(4)可以寫為

    (5)

    (6)

    式中,E(b)定義為時(shí)間-小波能量譜。

    1.2 損傷識(shí)別指標(biāo)的定義與損傷識(shí)別流程

    盾構(gòu)隧道由于服役環(huán)境復(fù)雜、列車運(yùn)行密度大、環(huán)境荷載影響等多因素長(zhǎng)期作用,隧道結(jié)構(gòu)難免會(huì)出現(xiàn)一系列的病害,主要包括不均勻沉降、橫向收斂變形、滲漏水、管片錯(cuò)臺(tái)、接頭損壞、襯砌裂縫、螺栓銹蝕、螺栓松動(dòng)、混凝土劣化、混凝土脫落等。

    其中襯砌裂縫、螺栓銹蝕和松動(dòng)、混凝土劣化和脫落等病害會(huì)降低盾構(gòu)襯砌管片的局部剛度和變形能力[11]。盾構(gòu)襯砌管片產(chǎn)生裂縫部位的剛度會(huì)下降;混凝土碳化會(huì)造成襯砌管片結(jié)構(gòu)表面發(fā)生損壞,內(nèi)部鋼筋銹蝕加快從而產(chǎn)生銹蝕裂縫,導(dǎo)致結(jié)構(gòu)的剛度降低;混凝土脫落剝離使得襯砌管片局部厚度變薄,襯砌管片剛度降低;螺栓松動(dòng)銹蝕會(huì)造成盾構(gòu)隧道管片接頭部位的剛度下降。盾構(gòu)隧道襯砌結(jié)構(gòu)由于上述病害的發(fā)展造成某一截面剛度下降形成損傷,當(dāng)車輛經(jīng)過盾構(gòu)隧道損傷位置時(shí),車輛和盾構(gòu)隧道之間的接觸力發(fā)生變化,車輛動(dòng)力響應(yīng)信號(hào)就會(huì)發(fā)生突變。小波變換是檢測(cè)信號(hào)奇異性的有效手段,利用小波變換可以檢測(cè)出信號(hào)發(fā)生奇異性的位置,從而判斷結(jié)構(gòu)發(fā)生損傷的位置。

    當(dāng)信號(hào)經(jīng)過小波變換后,小波系數(shù)會(huì)在信號(hào)突變位置的多個(gè)尺度范圍內(nèi)出現(xiàn)模極大值點(diǎn),可以通過某一尺度的小波系數(shù)圖來判斷模極大值的位置。為了充分利用多個(gè)尺度范圍內(nèi)的損傷信息,提高損傷識(shí)別的準(zhǔn)確率和靈敏性,可通過式(6)計(jì)算多個(gè)尺度范圍內(nèi)的小波能量,定義損傷識(shí)別指標(biāo)如下:

    (7)

    式中:DIm表示第m個(gè)盾構(gòu)隧道單元的損傷識(shí)別指標(biāo);N為車輛經(jīng)過第m個(gè)盾構(gòu)隧道單元的采樣點(diǎn)數(shù)。

    車輛加速度響應(yīng)信號(hào)進(jìn)行連續(xù)小波變換得到小波系數(shù),計(jì)算各個(gè)時(shí)間步一定尺度范圍內(nèi)的小波能量E(n),車輛經(jīng)過第m個(gè)盾構(gòu)隧道單元時(shí)間段內(nèi)的小波能量之和即為該單元的損傷識(shí)別指標(biāo)DIm,損傷單元處的損傷指標(biāo)會(huì)遠(yuǎn)大于未損傷單元處,可由此來識(shí)別和定位損傷,具體的損傷識(shí)別流程如圖1所示。

    圖1 盾構(gòu)隧道損傷識(shí)別流程

    2 數(shù)值模型

    以車輛加速度響應(yīng)信號(hào)作為損傷信息數(shù)據(jù),基于小波能量譜構(gòu)建損傷識(shí)別指標(biāo),盾構(gòu)隧道損傷間接識(shí)別示意圖如圖2所示。

    在數(shù)值模型中,車輛簡(jiǎn)化為單自由度三參數(shù)模型,盾構(gòu)隧道模擬為連續(xù)支承于黏彈性地基上的簡(jiǎn)支梁,簡(jiǎn)支梁為考慮Rayleigh阻尼的Euler-Bernoulli梁。通過在梁和地基之間施加彈簧來模擬圍巖和盾構(gòu)隧道結(jié)構(gòu)之間的相互作用[12],移動(dòng)車輛和盾構(gòu)隧道耦合振動(dòng)的理論模型如圖3所示。

    圖3 移動(dòng)車輛和盾構(gòu)隧道耦合振動(dòng)模型

    車輛系統(tǒng)的振動(dòng)方程為

    (8)

    盾構(gòu)隧道系統(tǒng)的振動(dòng)方程為

    (9)

    盾構(gòu)隧道單元的質(zhì)量矩陣為

    (10)

    式中,ρ為盾構(gòu)隧道單位長(zhǎng)度質(zhì)量。

    盾構(gòu)隧道單元的剛度矩陣由盾構(gòu)隧道襯砌彎曲勢(shì)能引起的剛度和圍巖連續(xù)彈性支承引起的剛度組成

    (11)

    (12)

    (13)

    式(12)和(13)中,EI為盾構(gòu)隧道的抗彎剛度,kr為圍巖的剛度系數(shù)。

    盾構(gòu)隧道單元的阻尼矩陣由盾構(gòu)隧道襯砌內(nèi)摩擦引起的阻尼和圍巖支承引起的阻尼組成

    (14)

    (15)

    式中,cr為圍巖的阻尼系數(shù)。

    車輛系統(tǒng)的荷載Fv和盾構(gòu)隧道系統(tǒng)的荷載Fb計(jì)算如下

    (16)

    (17)

    在MATLAB軟件中建立移動(dòng)車輛和盾構(gòu)隧道耦合振動(dòng)的有限元模型,采用基于Newmark格式的交叉迭代算法[13]進(jìn)行求解得到車輛和盾構(gòu)隧道的動(dòng)力響應(yīng)。

    3 算例分析

    車輛的質(zhì)量為mv=3 000 kg,剛度為kv=300 kN/m,阻尼為cv=1 000 N·s/m。

    盾構(gòu)隧道的參數(shù)選取自文獻(xiàn)[14],盾構(gòu)隧道襯砌管片環(huán)由C50鋼筋混凝土制成,管片環(huán)外直徑6 000 mm,內(nèi)直徑5 400 mm,管片厚度300 mm,幅寬1 500 mm。取100環(huán)管片長(zhǎng)的區(qū)間盾構(gòu)隧道為研究對(duì)象,盾構(gòu)隧道的長(zhǎng)度為L(zhǎng)=150 m,盾構(gòu)隧道襯砌結(jié)構(gòu)抗彎剛度為EI=7.55×105MPa·m4,盾構(gòu)隧道襯砌單位長(zhǎng)度質(zhì)量為ρ=13 430 kg/m,盾構(gòu)隧道襯砌阻尼采用Rayleigh阻尼,阻尼比設(shè)為0.01,盾構(gòu)隧道圍巖單位長(zhǎng)度的剛度系數(shù)為kr=1.1×108N/m2,盾構(gòu)隧道圍巖單位長(zhǎng)度的阻尼系數(shù)為cr=1.0×104N·s/m2。

    在移動(dòng)車輛和盾構(gòu)隧道耦合振動(dòng)有限元模型中,盾構(gòu)隧道按盾構(gòu)襯砌管片寬度1.5 m劃分一個(gè)單元,共劃分為100個(gè)單元,依次編號(hào)為1~100,車輛以5 m/s的速度勻速通過盾構(gòu)隧道,采樣頻率為1 000 Hz。

    對(duì)于結(jié)構(gòu)的損傷,一般認(rèn)為它只改變單元的剛度,而不改變單元的質(zhì)量,所以可以通過假定某一單元或者某些單元不同程度的剛度折減來模擬結(jié)構(gòu)的損傷。設(shè)第i個(gè)單元發(fā)生損傷前的抗彎剛度為EIi0,剛度折減系數(shù)為α,則損傷狀態(tài)下的單元?jiǎng)偠菶Ii可以表示為

    EIi=EIi0(1-α)

    (18)

    結(jié)構(gòu)損傷識(shí)別的目標(biāo)是對(duì)損傷進(jìn)行定位和定量分析。設(shè)定如表1所示的盾構(gòu)隧道損傷工況,考慮不同位置、不同程度的盾構(gòu)隧道單損傷工況以及多損傷工況,分析采用小波能量譜的方法,通過車輛加速度響應(yīng)識(shí)別盾構(gòu)隧道損傷的效果。

    表1 損傷工況

    對(duì)工況1下的車輛加速度響應(yīng)做連續(xù)小波變換,對(duì)比了幾種不同的小波函數(shù)進(jìn)行信號(hào)處理,選擇最適合本文研究的mexh小波函數(shù),尺度為1~600,可以得到如圖4所示。

    圖4 工況1車輛加速度響應(yīng)連續(xù)小波變換

    從圖4可以看出,工況1下,車輛加速度響應(yīng)經(jīng)過連續(xù)小波變換后,在9 s附近即車輛經(jīng)過損傷單元附近時(shí),出現(xiàn)了一條明顯的亮條紋,說明了車輛加速度響應(yīng)信號(hào)發(fā)生了突變,意味著損傷的存在。由于小波分析邊界效應(yīng)的影響,因此在端部位置會(huì)出現(xiàn)明顯的亮條紋。

    當(dāng)盾構(gòu)隧道結(jié)構(gòu)出現(xiàn)損傷時(shí),含有損傷信號(hào)的車輛加速度響應(yīng)經(jīng)過連續(xù)小波變換后會(huì)在多個(gè)尺度上產(chǎn)生損傷信息,因此可以通過計(jì)算某一尺度范圍內(nèi)的小波能量,充分利用多個(gè)尺度范圍內(nèi)的損傷信息,提高損傷識(shí)別的效果。選擇尺度350~600,計(jì)算各個(gè)時(shí)間步350~600尺度內(nèi)的小波能量,計(jì)算各個(gè)盾構(gòu)隧道單元的損傷指標(biāo),可以得到損傷工況1的小波能量分布圖如圖5所示。

    圖5 工況1小波能量分布圖

    從圖5可以看出,由于小波分析邊界效應(yīng)的影響,在端部會(huì)出現(xiàn)較高的能量聚集,可以截取中間的20~80號(hào)單元區(qū)間內(nèi)的數(shù)據(jù)來分析。工況1和工況2的小波能量分布圖如圖6和圖7所示,在第30號(hào)單元和第50號(hào)單元處存在小波能量峰值,說明第30號(hào)單元和第50號(hào)單元發(fā)生了損傷,與工況1、2設(shè)定的損傷單元相吻合。工況1、2的損傷識(shí)別結(jié)果表明,通過車輛加速度響應(yīng)計(jì)算的盾構(gòu)隧道損傷指標(biāo),對(duì)于盾構(gòu)隧道不同位置處發(fā)生的單損傷均可以準(zhǔn)確識(shí)別。

    圖6 工況1小波能量分布圖

    圖7 工況2小波能量分布圖

    工況2~4是不同損傷程度的單損傷工況對(duì)比,計(jì)算盾構(gòu)隧道損傷識(shí)別指標(biāo)得到的小波能量分布圖如圖8所示,在不同程度的盾構(gòu)隧道損傷工況下,均可以準(zhǔn)確識(shí)別出第50號(hào)單元發(fā)生了損傷。隨著單元損傷程度的增大,小波能量值也不斷增大,說明小波能量值的大小可以反映盾構(gòu)隧道結(jié)構(gòu)的損傷程度。

    圖8 工況2-4小波能量分布圖

    工況5是損傷程度相同的多損傷工況,工況6是不同損傷程度的多損傷工況,計(jì)算盾構(gòu)隧道損傷識(shí)別指標(biāo)得到的小波能量分布圖如圖9和圖10所示。從圖9和圖10可以看出,在盾構(gòu)隧道發(fā)生多處損傷時(shí),通過小波能量分布圖均可以準(zhǔn)確識(shí)別出損傷發(fā)生的位置。工況5結(jié)果表明,多個(gè)單元發(fā)生相同程度的損傷時(shí),其小波能量的峰值并不相同,存在較小的差異;工況6結(jié)果表明,在盾構(gòu)隧道發(fā)生多處不同程度的損傷時(shí),根據(jù)小波能量的峰值可以判斷損傷程度的相對(duì)大小,小波能量峰值越大,發(fā)生的損傷程度越大。

    圖9 工況5小波能量分布圖

    圖10 工況6小波能量分布圖

    工況7和工況8是為了分析盾構(gòu)隧道發(fā)生損傷的單元相鄰近時(shí)的損傷識(shí)別效果,計(jì)算盾構(gòu)隧道損傷識(shí)別指標(biāo)得到的小波能量分布圖如圖11和圖12所示。從圖11可以看出,當(dāng)鄰近的單元發(fā)生相同程度的損傷時(shí),小波能量分布圖上有三個(gè)明顯的波峰,對(duì)應(yīng)著損傷發(fā)生的位置;從圖12可以看出,當(dāng)鄰近的單元發(fā)生不同程度的損傷時(shí),從小波能量分布圖上也可以看到三個(gè)比較尖銳的波峰對(duì)應(yīng)著損傷發(fā)生的位置。從圖11和圖12也可以看出,在損傷單元附近的小波能量也會(huì)出現(xiàn)一些小的波峰,這是采用的mexh小波函數(shù)的性質(zhì)引起的。發(fā)生損傷的單元對(duì)應(yīng)的小波能量的波峰是比較尖銳的,而損傷單元附近的小波能量波峰是比較平緩的,可以通過這些差異來排除損傷單元附近小波能量波峰的影響,減少損傷識(shí)別的誤判。

    圖11 工況7小波能量分布圖

    圖12 工況8小波能量分布圖

    4 參數(shù)分析

    4.1 車輛速度

    選擇表1中設(shè)置的盾構(gòu)隧道多損傷工況6,車輛的速度分別設(shè)置為2 m/s、5 m/s、10 m/s、15 m/s,車輛勻速通過盾構(gòu)隧道,采樣頻率均設(shè)置為1 000 Hz,其余參數(shù)不變。不同車速下計(jì)算盾構(gòu)隧道損傷識(shí)別指標(biāo)得到的小波能量分布圖,如圖13所示。

    從圖13(a)和圖13(b)可以看出,當(dāng)車速較小時(shí),通過小波能量分布圖可以準(zhǔn)確地識(shí)別出第35、50、70號(hào)單元發(fā)生了損傷,特別是當(dāng)車速只有2 m/s時(shí),只在損傷單元處出現(xiàn)了明顯的小波能量波峰,而其余未損傷單元處的小波能量值趨于零。當(dāng)車速增大到10 m/s時(shí),從圖13(c)可以看出,第35、50、70號(hào)單元附近存在著小波能量波峰,但相比于其他單元處的能量波峰已經(jīng)不是非常明顯,很難準(zhǔn)確識(shí)別出發(fā)生損傷的位置。而當(dāng)車速增大到15 m/s時(shí),從圖13(d)中已經(jīng)無法識(shí)別出發(fā)生損傷的位置。因此,在采用小波能量譜方法間接識(shí)別盾構(gòu)隧道的損傷時(shí),采用較低的車速可以獲得更好的損傷識(shí)別效果。

    4.2 測(cè)量噪聲

    為了研究測(cè)量噪聲對(duì)損傷識(shí)別結(jié)果的影響,在有限元模型計(jì)算得到的車輛加速度響應(yīng)中加入噪聲模擬現(xiàn)場(chǎng)測(cè)試中被噪聲污染的試驗(yàn)數(shù)據(jù)。

    am=ac+epN0σ(ac)

    (19)

    式中:ac為計(jì)算車輛加速度響應(yīng);am為加入噪聲的車輛加速度響應(yīng);ep為噪聲水平0≤ep≤1;N0為與計(jì)算車輛加速度響應(yīng)ac等長(zhǎng)度的標(biāo)準(zhǔn)正態(tài)分布隨機(jī)向量;σ(ac)為計(jì)算車輛加速度響應(yīng)ac的標(biāo)準(zhǔn)差。

    同樣選擇表1中的盾構(gòu)隧道多損傷工況6,車輛的速度為5 m/s,其余參數(shù)設(shè)置不變。在有限元模型計(jì)算得到的車輛加速度響應(yīng)中分別添加0%、1%、2%和5%水平的噪聲。不同噪聲水平下計(jì)算盾構(gòu)隧道損傷識(shí)別指標(biāo)得到的小波能量分布圖如圖14所示。

    從圖14可以看出,在車輛加速度響應(yīng)受到不同程度測(cè)量噪聲污染下,計(jì)算盾構(gòu)隧道損傷識(shí)別指標(biāo)得到的小波能量分布圖差異較小,從小波能量分布圖中均可以準(zhǔn)確識(shí)別出盾構(gòu)隧道發(fā)生損傷的位置,說明測(cè)量噪聲對(duì)損傷識(shí)別結(jié)果造成的影響較小。

    圖14 不同噪聲工況下小波能量分布圖

    4.3 軌道粗糙度

    粗糙度是間接測(cè)量法應(yīng)用的重要影響因素,根據(jù)功率譜密度函數(shù)以不同空間頻率的余弦函數(shù)疊加來模擬軌道粗糙度[15],即:

    (20)

    式中:ns,i為第i個(gè)空間頻率;di、θi為第i個(gè)余弦函數(shù)的振幅與隨機(jī)相位角。空間頻率ns選用1~100 r/m。

    隨機(jī)生成一組A級(jí)粗糙度如圖15所示。在移動(dòng)車輛和盾構(gòu)隧道耦合振動(dòng)有限元模型中添加軌道粗糙度,有無粗糙度影響下的車輛加速度響應(yīng)對(duì)比如圖16所示。從圖16可以看出,有無粗糙度下車輛加速度響應(yīng)幅值相差特別大,在考慮粗糙度影響后,車輛加速度響應(yīng)中主要是由軌道粗糙度引起的響應(yīng),而由車輛和盾構(gòu)隧道耦合振動(dòng)引起的響應(yīng)非常微小,這樣就無法從單次測(cè)量的車輛加速度響應(yīng)中基于小波能量譜構(gòu)建的損傷指標(biāo)識(shí)別盾構(gòu)隧道的損傷。

    圖15 A級(jí)粗糙度

    圖16 有無粗糙度下車輛加速度響應(yīng)對(duì)比

    為了降低粗糙度的影響,首先是要提高車輛和盾構(gòu)隧道耦合振動(dòng)引起的響應(yīng)在車輛加速度響應(yīng)中的占比。間接測(cè)量法的一個(gè)優(yōu)點(diǎn)就是檢測(cè)車輛既可以作為激勵(lì)器又可以作為接收器,但一般檢測(cè)車輛作為單自由度車輛不宜質(zhì)量過大。因此參考文獻(xiàn)[16],可以采用一輛重車牽引的方式,來對(duì)盾構(gòu)隧道進(jìn)行激勵(lì),從而加大檢測(cè)車輛和盾構(gòu)隧道耦合振動(dòng)引起的響應(yīng),重車考慮40 t的重量,輪距為5 m,檢測(cè)車輛和重車的距離為2 m,具體如圖17所示。

    圖17 重車牽引示意圖

    雖然重車激勵(lì)可以增大車輛與盾構(gòu)隧道耦合振動(dòng)引起的響應(yīng)在車輛加速度響應(yīng)中的占比,但是粗糙度引起的響應(yīng)仍然占比較大,無法通過單次測(cè)量得到的車輛加速度響應(yīng)判別損傷。由前面的分析可知,當(dāng)盾構(gòu)隧道結(jié)構(gòu)狀態(tài)發(fā)生變化時(shí),車輛加速度響應(yīng)進(jìn)行連續(xù)小波變換,計(jì)算盾構(gòu)隧道損傷識(shí)別指標(biāo)得到的小波能量分布也會(huì)發(fā)生變化。因此可以考慮在重車激勵(lì)下,利用盾構(gòu)隧道發(fā)生損傷前后計(jì)算得到的小波能量差來構(gòu)造損傷指標(biāo)??紤]粗糙度影響的盾構(gòu)隧道損傷指標(biāo)構(gòu)造流程如圖18所示。

    圖18 損傷指標(biāo)構(gòu)造流程圖

    選擇表1中盾構(gòu)隧道多損傷工況6中設(shè)置的第35、50、70號(hào)單元分別發(fā)生了5%、10%、8%的損傷,車輛的速度設(shè)置為5 m/s,其余參數(shù)設(shè)置不變。在移動(dòng)車輛和盾構(gòu)隧道耦合振動(dòng)有限元模型中添加軌道粗糙度,在車輛加速度響應(yīng)中添加0%、1%、2%和5%不同水平的噪聲。考慮粗糙度影響下,基于小波能量差的盾構(gòu)隧道損傷識(shí)別結(jié)果如圖19所示。

    從圖19(a)中可以看出,在無噪聲水平下,損傷單元位置處的小波能量差存在著明顯的波峰,其余未損傷單元處的小波能量差非常小;從圖19(b)中可以看出,在1%噪聲水平下,在損傷單元位置處的小波能量差存在著較為明顯的波峰,其余未損傷單元處的小波能量差也出現(xiàn)了一些小的波峰;從圖19(c)中可以看出,在2%噪聲水平下,在第50號(hào)和70號(hào)單元位置處的小波能量差存在著較為明顯的波峰,而第35號(hào)單元附近沒有明顯的波峰;從圖19(d)中可以看出,在5%噪聲水平下,存在著較多的波峰,已經(jīng)無法判斷出損傷單元的位置。

    通過圖19中考慮粗糙度影響和不同噪聲水平下盾構(gòu)隧道損傷識(shí)別結(jié)果的對(duì)比可知,無噪聲水平下可以準(zhǔn)確識(shí)別出損傷單元的位置,驗(yàn)證了損傷識(shí)別指標(biāo)的有效性。在低噪聲水平下,損傷單元處的小波能量差會(huì)大于未損傷單元處的小波能量差,可以通過小波能量差波峰的相對(duì)大小來判斷損傷單元的位置。隨著噪聲水平的增加,損傷程度較小處的小波能量差波峰會(huì)被淹沒,在高噪聲水平下,出現(xiàn)了眾多的小波能量差波峰,無法識(shí)別出損傷單元。

    (a) 0%噪聲

    從上述分析可知,在考慮粗糙度影響時(shí),噪聲對(duì)本方法的影響較大,這是因?yàn)樵肼暿翘砑釉谲囕v加速度響應(yīng)上的,而盾構(gòu)隧道的損傷信息是存在于車隧耦合振動(dòng)引起的響應(yīng)之中,車隧耦合振動(dòng)引起的響應(yīng)在車輛加速度響應(yīng)中的幅值占比較小,因此噪聲對(duì)車隧耦合振動(dòng)引起的響應(yīng)附加的影響就會(huì)更大,盾構(gòu)隧道的損傷信息很容易被噪聲淹沒,就無法識(shí)別出損傷。因此考慮粗糙度影響時(shí),通過損傷前后小波能量差的方法識(shí)別盾構(gòu)隧道損傷時(shí),應(yīng)保持較低的噪聲水平以獲得較好的損傷識(shí)別效果。

    5 試驗(yàn)驗(yàn)證

    本節(jié)通過實(shí)驗(yàn)室試驗(yàn)驗(yàn)證所提盾構(gòu)隧道損傷間接識(shí)別方法的有效性。使用鋼板、彈簧、車輪等制作了如圖20所示的試驗(yàn)車輛,試驗(yàn)車輛分為牽引平衡車和測(cè)量車,在牽引平衡車上放置一個(gè)500 g和兩個(gè)200 g的砝碼以加大車輛的激勵(lì)。盾構(gòu)隧道襯砌采用PVC管模擬,PVC管長(zhǎng)度4 000 mm,外徑160 mm,管壁厚3 mm,管道上隨機(jī)涂抹石灰以模擬軌道表面粗糙度,在管道上粘貼玻璃纖維扁條以防止車輛在行進(jìn)過程中滑落。盾構(gòu)隧道圍巖環(huán)境采用沙箱模擬,將管道埋入沙子中一半。電機(jī)通過細(xì)繩牽引車輛從管道一端勻速行進(jìn)到另一端,使用東華DHDAS動(dòng)態(tài)信號(hào)測(cè)試儀采集測(cè)量車豎向振動(dòng)加速度響應(yīng)信號(hào),試驗(yàn)測(cè)試現(xiàn)場(chǎng)如圖21所示。

    圖20 試驗(yàn)車輛

    圖21 試驗(yàn)現(xiàn)場(chǎng)

    本文提出的盾構(gòu)隧道損傷間接識(shí)別方法適用于盾構(gòu)隧道結(jié)構(gòu)在襯砌裂縫、混凝土劣化等病害發(fā)展下襯砌管片剛度下降的損傷。在實(shí)驗(yàn)室試驗(yàn)中多采用試驗(yàn)構(gòu)件開缺口或切割的方式使得截面剛度下降以模擬結(jié)構(gòu)損傷[17-19],因此本文試驗(yàn)將管道平均劃分為20個(gè)單元,依次編號(hào)1~20,在第8號(hào)單元上切割5 cm×5 cm的方形缺口使得截面剛度損失5%以模擬盾構(gòu)隧道襯砌損傷[20],如圖22所示。試驗(yàn)在夜間進(jìn)行以保證較低的測(cè)量噪聲,分別采集車輛在管道無損工況和有損工況下通過管道的豎向振動(dòng)加速度響應(yīng)信號(hào),采樣頻率為1 000 Hz。

    圖22 損傷設(shè)置形式

    選擇mexh小波函數(shù),尺度為1~600,對(duì)無損和有損工況下采集到的車輛加速度響應(yīng)信號(hào)做連續(xù)小波變換,分別計(jì)算無損和有損工況下尺度350~600范圍內(nèi)的小波能量譜,計(jì)算小波能量差。選擇中間4~18號(hào)單元的數(shù)據(jù),可得到損傷識(shí)別結(jié)果如圖23所示。

    從圖23中可以看出,第8號(hào)單元處的小波能量差遠(yuǎn)大于其他單元處的小波能量差,因此可以判斷第8號(hào)單元處發(fā)生了損傷,這與預(yù)設(shè)的損傷位置相吻合。由于試驗(yàn)條件限制,本文試驗(yàn)?zāi)P臀磳⒐艿廊柯裰糜谏匙又?,試?yàn)車輛未置于管道內(nèi)部行駛,但該試驗(yàn)?zāi)P涂梢苑从忱碚撃P椭熊囕v和盾構(gòu)隧道豎向的耦合振動(dòng),試驗(yàn)采集得到的是車輛豎向振動(dòng)加速度響應(yīng)信號(hào),所以通過該試驗(yàn)?zāi)P涂梢猿醪津?yàn)證本文所提盾構(gòu)隧道損傷間接識(shí)別方法的可行性。

    圖23 損傷識(shí)別結(jié)果

    6 結(jié) 論

    本文結(jié)合間接測(cè)量法和小波分析方法,提出了一種基于小波能量譜的盾構(gòu)隧道損傷間接識(shí)別方法,主要結(jié)論如下:

    (1) 不考慮粗糙度影響時(shí),對(duì)于盾構(gòu)隧道單損傷和多損傷工況,本文方法可以有效識(shí)別損傷發(fā)生的位置,損傷指標(biāo)可以反映損傷程度的相對(duì)大小。采用較低的車速可以獲得更好的損傷識(shí)別效果,具有較好的噪聲魯棒性;

    (2) 考慮粗糙度影響時(shí),提出了在重車激勵(lì)下構(gòu)建損傷前后小波能量差的損傷識(shí)別指標(biāo),數(shù)值算例結(jié)果表明,在低噪聲條件下,可以準(zhǔn)確識(shí)別盾構(gòu)隧道損傷發(fā)生的位置,并且損傷指標(biāo)可以反映損傷程度的相對(duì)大?。?/p>

    (3) 在實(shí)驗(yàn)室搭建了車輛和盾構(gòu)隧道耦合振動(dòng)模型,通過損傷識(shí)別試驗(yàn)驗(yàn)證了本文所提出的基于小波能量譜的盾構(gòu)隧道損傷間接識(shí)別方法的有效性。

    猜你喜歡
    小波粗糙度盾構(gòu)
    構(gòu)造Daubechies小波的一些注記
    基于無人機(jī)影像的巖體結(jié)構(gòu)面粗糙度獲取
    甘肅科技(2020年20期)2020-04-13 00:30:18
    基于MATLAB的小波降噪研究
    電子制作(2019年13期)2020-01-14 03:15:32
    冷沖模磨削表面粗糙度的加工試驗(yàn)與應(yīng)用
    模具制造(2019年4期)2019-06-24 03:36:48
    基于改進(jìn)的G-SVS LMS 與冗余提升小波的滾動(dòng)軸承故障診斷
    基于BP神經(jīng)網(wǎng)絡(luò)的面齒輪齒面粗糙度研究
    鋼材銹蝕率與表面三維粗糙度參數(shù)的關(guān)系
    盾構(gòu)近距離下穿房屋接收技術(shù)
    復(fù)合盾構(gòu)在縱向錨桿區(qū)的掘進(jìn)分析及實(shí)踐
    基于FPGA小波變換核的設(shè)計(jì)
    日韩欧美国产在线观看| 色av中文字幕| 亚洲国产精品成人综合色| 久久久国产精品麻豆| 自线自在国产av| 亚洲精品在线观看二区| 老熟妇仑乱视频hdxx| 国产真人三级小视频在线观看| 最好的美女福利视频网| www.999成人在线观看| 性色av乱码一区二区三区2| 高清毛片免费观看视频网站| x7x7x7水蜜桃| 亚洲男人的天堂狠狠| 久久国产精品人妻蜜桃| 久久国产精品影院| 亚洲国产高清在线一区二区三 | 黄色女人牲交| 99久久无色码亚洲精品果冻| 老司机福利观看| 国产1区2区3区精品| 一级黄色大片毛片| 日韩 欧美 亚洲 中文字幕| 国产精品免费一区二区三区在线| 亚洲一码二码三码区别大吗| 99久久国产精品久久久| 精品一区二区三区av网在线观看| 黑人巨大精品欧美一区二区mp4| 无遮挡黄片免费观看| 国产一区二区在线av高清观看| 日日干狠狠操夜夜爽| 日本一区二区免费在线视频| 91老司机精品| 色老头精品视频在线观看| 亚洲全国av大片| 欧美午夜高清在线| 久久精品成人免费网站| 国产精品永久免费网站| 欧美激情久久久久久爽电影| 啦啦啦免费观看视频1| 亚洲熟妇熟女久久| 日韩大码丰满熟妇| 国产片内射在线| 中文字幕高清在线视频| 男人舔女人的私密视频| 天堂动漫精品| 国产av又大| 欧美色欧美亚洲另类二区| xxx96com| 亚洲精品国产区一区二| 男人舔奶头视频| 99国产精品99久久久久| 天天躁狠狠躁夜夜躁狠狠躁| 久久久国产精品麻豆| 日韩欧美在线二视频| 欧美激情高清一区二区三区| 啦啦啦免费观看视频1| 色播在线永久视频| 巨乳人妻的诱惑在线观看| 热re99久久国产66热| 很黄的视频免费| 亚洲av日韩精品久久久久久密| 成熟少妇高潮喷水视频| 精品国产一区二区三区四区第35| 无限看片的www在线观看| 国产亚洲精品久久久久5区| 色老头精品视频在线观看| 女警被强在线播放| 男人舔女人下体高潮全视频| 在线播放国产精品三级| 亚洲av成人av| 女生性感内裤真人,穿戴方法视频| av在线天堂中文字幕| 久9热在线精品视频| 中文字幕最新亚洲高清| 一区二区三区激情视频| 日韩成人在线观看一区二区三区| 国产成人精品久久二区二区免费| 国产精品国产高清国产av| 亚洲av电影在线进入| 国产熟女xx| 亚洲欧美日韩高清在线视频| 桃色一区二区三区在线观看| 老鸭窝网址在线观看| 国产精品精品国产色婷婷| 一二三四社区在线视频社区8| 国产男靠女视频免费网站| 少妇 在线观看| 曰老女人黄片| 亚洲成人久久性| 88av欧美| 午夜老司机福利片| 天天一区二区日本电影三级| 1024香蕉在线观看| 在线观看免费午夜福利视频| 一本综合久久免费| 精品国产美女av久久久久小说| 久久天躁狠狠躁夜夜2o2o| √禁漫天堂资源中文www| 国产精品99久久99久久久不卡| 男女下面进入的视频免费午夜 | 午夜福利成人在线免费观看| 久久香蕉精品热| 黄片小视频在线播放| 视频在线观看一区二区三区| 脱女人内裤的视频| 99国产极品粉嫩在线观看| 嫩草影院精品99| av片东京热男人的天堂| www日本在线高清视频| 久久人人精品亚洲av| 91麻豆精品激情在线观看国产| 国产一区二区三区在线臀色熟女| 天堂动漫精品| 很黄的视频免费| 熟妇人妻久久中文字幕3abv| 99国产极品粉嫩在线观看| 在线免费观看的www视频| 男人的好看免费观看在线视频 | 精品日产1卡2卡| 亚洲va日本ⅴa欧美va伊人久久| 久久精品91蜜桃| 日韩欧美三级三区| 又大又爽又粗| 国产激情偷乱视频一区二区| 国产成人欧美| 亚洲一区二区三区色噜噜| 国产精品久久久久久精品电影 | 精品一区二区三区av网在线观看| 国产亚洲精品久久久久久毛片| 欧美黄色片欧美黄色片| 日韩大码丰满熟妇| 免费在线观看成人毛片| 午夜激情av网站| 亚洲国产精品999在线| 美女高潮到喷水免费观看| 俺也久久电影网| 色av中文字幕| 婷婷亚洲欧美| 男人舔女人下体高潮全视频| 日本免费a在线| 成熟少妇高潮喷水视频| 两性午夜刺激爽爽歪歪视频在线观看 | bbb黄色大片| 国产人伦9x9x在线观看| 亚洲三区欧美一区| 久久天堂一区二区三区四区| 午夜免费激情av| 欧美日韩亚洲国产一区二区在线观看| 中文字幕久久专区| 首页视频小说图片口味搜索| 无限看片的www在线观看| 久久久久久久久免费视频了| 欧美日韩一级在线毛片| 国产一区二区三区在线臀色熟女| 黄片小视频在线播放| 97人妻精品一区二区三区麻豆 | 男女那种视频在线观看| 欧美人与性动交α欧美精品济南到| 人人妻人人澡人人看| 日韩高清综合在线| 成人亚洲精品一区在线观看| 久久中文看片网| cao死你这个sao货| 日韩欧美在线二视频| 欧美三级亚洲精品| 亚洲精品色激情综合| av视频在线观看入口| 亚洲欧洲精品一区二区精品久久久| 美女国产高潮福利片在线看| 国产精品久久久久久精品电影 | 国产1区2区3区精品| 妹子高潮喷水视频| 女性生殖器流出的白浆| 很黄的视频免费| 亚洲精品av麻豆狂野| ponron亚洲| 亚洲精品粉嫩美女一区| 欧美成人午夜精品| 两个人免费观看高清视频| 18禁美女被吸乳视频| 非洲黑人性xxxx精品又粗又长| 免费在线观看日本一区| 免费在线观看影片大全网站| 99在线人妻在线中文字幕| 国产97色在线日韩免费| 色播在线永久视频| 国产欧美日韩一区二区三| 母亲3免费完整高清在线观看| 国产麻豆成人av免费视频| 天天一区二区日本电影三级| 国产av一区二区精品久久| 日本a在线网址| 黄色丝袜av网址大全| 又黄又爽又免费观看的视频| 国产又爽黄色视频| 人妻丰满熟妇av一区二区三区| 久久青草综合色| 老汉色∧v一级毛片| www.999成人在线观看| 久久草成人影院| 中文字幕精品免费在线观看视频| 91字幕亚洲| 人成视频在线观看免费观看| www.熟女人妻精品国产| 午夜激情av网站| 国产高清激情床上av| 免费搜索国产男女视频| 国内久久婷婷六月综合欲色啪| 变态另类成人亚洲欧美熟女| 亚洲国产欧洲综合997久久, | 亚洲中文字幕一区二区三区有码在线看 | 国产精品 欧美亚洲| 国产乱人伦免费视频| 国产主播在线观看一区二区| 可以在线观看的亚洲视频| 亚洲中文日韩欧美视频| 波多野结衣高清无吗| 亚洲国产精品成人综合色| √禁漫天堂资源中文www| 伊人久久大香线蕉亚洲五| 国产一区二区三区在线臀色熟女| 欧洲精品卡2卡3卡4卡5卡区| 久久久国产欧美日韩av| 亚洲欧美一区二区三区黑人| 亚洲av成人一区二区三| 无遮挡黄片免费观看| 一夜夜www| 成人三级做爰电影| 成人永久免费在线观看视频| 国产精品一区二区三区四区久久 | 很黄的视频免费| 麻豆一二三区av精品| 正在播放国产对白刺激| 国产成人系列免费观看| 亚洲真实伦在线观看| 欧美在线黄色| videosex国产| 视频在线观看一区二区三区| 一本一本综合久久| 97超级碰碰碰精品色视频在线观看| 国产成人啪精品午夜网站| 久久久久免费精品人妻一区二区 | 亚洲男人的天堂狠狠| 精品久久久久久久久久免费视频| 精品国内亚洲2022精品成人| 国产亚洲精品一区二区www| 黄网站色视频无遮挡免费观看| 大型黄色视频在线免费观看| 香蕉av资源在线| 两性夫妻黄色片| 日日夜夜操网爽| ponron亚洲| 色综合站精品国产| www.精华液| 亚洲熟女毛片儿| 国产av不卡久久| 久久久精品国产亚洲av高清涩受| 中文资源天堂在线| 国产真人三级小视频在线观看| 国产黄a三级三级三级人| 午夜精品在线福利| 日本一本二区三区精品| 亚洲电影在线观看av| 中文字幕人妻丝袜一区二区| 亚洲狠狠婷婷综合久久图片| 国产成人系列免费观看| 狂野欧美激情性xxxx| 黑人欧美特级aaaaaa片| 欧美大码av| 最近最新中文字幕大全电影3 | 国产精品一区二区免费欧美| 一级毛片高清免费大全| 国产伦在线观看视频一区| 亚洲精品在线美女| 91字幕亚洲| 亚洲成人精品中文字幕电影| 欧美激情久久久久久爽电影| 亚洲欧洲精品一区二区精品久久久| 亚洲国产中文字幕在线视频| 久久青草综合色| 亚洲成国产人片在线观看| 精品第一国产精品| 免费在线观看视频国产中文字幕亚洲| 免费高清视频大片| 两人在一起打扑克的视频| 久久精品91蜜桃| 一进一出好大好爽视频| 色尼玛亚洲综合影院| 亚洲国产欧美日韩在线播放| 香蕉国产在线看| 我的亚洲天堂| 国内揄拍国产精品人妻在线 | 一夜夜www| 国产久久久一区二区三区| 十八禁网站免费在线| 久久精品成人免费网站| 亚洲av第一区精品v没综合| 男女视频在线观看网站免费 | 91麻豆精品激情在线观看国产| www日本在线高清视频| 午夜视频精品福利| 日本免费一区二区三区高清不卡| 精品国产一区二区三区四区第35| 亚洲av第一区精品v没综合| 人人妻人人澡人人看| 伦理电影免费视频| 婷婷丁香在线五月| 一本一本综合久久| 老司机午夜福利在线观看视频| 亚洲avbb在线观看| 黄色丝袜av网址大全| 啦啦啦韩国在线观看视频| 天天一区二区日本电影三级| 午夜福利高清视频| 18禁黄网站禁片免费观看直播| 后天国语完整版免费观看| 久久午夜综合久久蜜桃| 成人一区二区视频在线观看| 中文字幕久久专区| 黄色成人免费大全| 一本一本综合久久| 精品无人区乱码1区二区| 最新美女视频免费是黄的| www日本在线高清视频| 亚洲国产精品sss在线观看| 国产主播在线观看一区二区| 亚洲男人天堂网一区| 成人三级黄色视频| 男女视频在线观看网站免费 | 草草在线视频免费看| 黄片大片在线免费观看| 老司机福利观看| 一级毛片高清免费大全| 侵犯人妻中文字幕一二三四区| 免费看a级黄色片| 久久久久久国产a免费观看| 国产黄a三级三级三级人| 欧美成人一区二区免费高清观看 | 日韩成人在线观看一区二区三区| 村上凉子中文字幕在线| 制服丝袜大香蕉在线| 两性夫妻黄色片| 99国产极品粉嫩在线观看| 日韩中文字幕欧美一区二区| 一级黄色大片毛片| 日本成人三级电影网站| 亚洲人成网站在线播放欧美日韩| 一二三四在线观看免费中文在| 手机成人av网站| 精品福利观看| 美国免费a级毛片| 亚洲狠狠婷婷综合久久图片| 麻豆国产av国片精品| 制服诱惑二区| 精品熟女少妇八av免费久了| 国产午夜精品久久久久久| 日韩精品免费视频一区二区三区| www日本黄色视频网| 在线十欧美十亚洲十日本专区| 久久精品国产99精品国产亚洲性色| 一级作爱视频免费观看| 人成视频在线观看免费观看| 99在线人妻在线中文字幕| 午夜影院日韩av| 亚洲成人久久性| 亚洲精品一区av在线观看| 久久99热这里只有精品18| 欧美黑人欧美精品刺激| 久久久水蜜桃国产精品网| 99精品欧美一区二区三区四区| 又黄又爽又免费观看的视频| 国产一级毛片七仙女欲春2 | 成人永久免费在线观看视频| 亚洲国产欧美网| 中文字幕久久专区| 色婷婷久久久亚洲欧美| 两性午夜刺激爽爽歪歪视频在线观看 | 99热这里只有精品一区 | 亚洲专区中文字幕在线| 悠悠久久av| 啦啦啦免费观看视频1| 国产91精品成人一区二区三区| 日韩中文字幕欧美一区二区| 婷婷精品国产亚洲av在线| 国产亚洲av高清不卡| 最好的美女福利视频网| 日韩精品免费视频一区二区三区| 十八禁网站免费在线| 中文字幕另类日韩欧美亚洲嫩草| 香蕉丝袜av| av福利片在线| 国产精品永久免费网站| 高潮久久久久久久久久久不卡| 色综合婷婷激情| 美女国产高潮福利片在线看| 久久香蕉精品热| 麻豆成人av在线观看| 少妇熟女aⅴ在线视频| 一区二区三区激情视频| 自线自在国产av| 国产精品综合久久久久久久免费| 欧美日韩黄片免| 日韩欧美国产一区二区入口| 亚洲自拍偷在线| 草草在线视频免费看| 中文字幕人妻丝袜一区二区| 日本免费a在线| 国产精品 国内视频| 午夜老司机福利片| xxxwww97欧美| 热re99久久国产66热| 国产精品98久久久久久宅男小说| 在线国产一区二区在线| 国产亚洲av高清不卡| av天堂在线播放| 国产精品久久久久久亚洲av鲁大| 欧美av亚洲av综合av国产av| 日韩欧美三级三区| 亚洲狠狠婷婷综合久久图片| 中亚洲国语对白在线视频| 一区二区三区精品91| 麻豆av在线久日| 成年人黄色毛片网站| 久久久国产成人精品二区| 黄片播放在线免费| 99精品在免费线老司机午夜| 好看av亚洲va欧美ⅴa在| 中文字幕久久专区| 国产乱人伦免费视频| 日韩成人在线观看一区二区三区| 一二三四社区在线视频社区8| 色综合站精品国产| 久久久久国内视频| 非洲黑人性xxxx精品又粗又长| 欧美国产日韩亚洲一区| 久久国产亚洲av麻豆专区| 欧美精品啪啪一区二区三区| 精品人妻1区二区| 88av欧美| 性色av乱码一区二区三区2| av福利片在线| 国产成人av激情在线播放| 又黄又爽又免费观看的视频| 天堂√8在线中文| 99re在线观看精品视频| 啦啦啦 在线观看视频| 色精品久久人妻99蜜桃| 久久婷婷成人综合色麻豆| 亚洲欧洲精品一区二区精品久久久| 欧美日本视频| 欧美+亚洲+日韩+国产| 91麻豆精品激情在线观看国产| 一区福利在线观看| 午夜福利免费观看在线| 一级作爱视频免费观看| 日本一本二区三区精品| 熟女电影av网| av在线播放免费不卡| 国产高清有码在线观看视频 | 国产高清有码在线观看视频 | 满18在线观看网站| 99riav亚洲国产免费| 欧美乱码精品一区二区三区| 此物有八面人人有两片| 日韩视频一区二区在线观看| 久久伊人香网站| 国产一级毛片七仙女欲春2 | 真人一进一出gif抽搐免费| 99精品在免费线老司机午夜| 亚洲精品中文字幕在线视频| 美女免费视频网站| 俄罗斯特黄特色一大片| 国产麻豆成人av免费视频| xxxwww97欧美| 午夜成年电影在线免费观看| 久久伊人香网站| 精品久久久久久久人妻蜜臀av| 黄网站色视频无遮挡免费观看| 精品久久久久久,| 中文字幕人妻丝袜一区二区| 真人一进一出gif抽搐免费| 日本一本二区三区精品| 色哟哟哟哟哟哟| 宅男免费午夜| 人人妻人人澡人人看| 亚洲av中文字字幕乱码综合 | 巨乳人妻的诱惑在线观看| 女生性感内裤真人,穿戴方法视频| 亚洲成av人片免费观看| 日本一本二区三区精品| 在线观看www视频免费| 在线免费观看的www视频| 免费人成视频x8x8入口观看| 久久久久国产一级毛片高清牌| 亚洲中文日韩欧美视频| 在线观看66精品国产| 99re在线观看精品视频| 免费看十八禁软件| 观看免费一级毛片| 国产高清视频在线播放一区| videosex国产| 国产成人av教育| 一本综合久久免费| 免费看a级黄色片| 少妇裸体淫交视频免费看高清 | 久久久久九九精品影院| 欧美国产日韩亚洲一区| 亚洲av成人一区二区三| 12—13女人毛片做爰片一| 在线视频色国产色| 国产精品日韩av在线免费观看| 亚洲一卡2卡3卡4卡5卡精品中文| 久久久久久久久免费视频了| 欧美日韩黄片免| 免费高清视频大片| 国产熟女午夜一区二区三区| 正在播放国产对白刺激| 老司机深夜福利视频在线观看| 日韩精品免费视频一区二区三区| 亚洲av熟女| 亚洲免费av在线视频| 国产精品野战在线观看| 精品人妻1区二区| 精品国产美女av久久久久小说| 久久精品成人免费网站| 日韩欧美国产在线观看| 亚洲精品在线美女| 国产欧美日韩精品亚洲av| 久久午夜亚洲精品久久| 我的亚洲天堂| 12—13女人毛片做爰片一| a级毛片a级免费在线| 国内精品久久久久久久电影| av电影中文网址| 久久99热这里只有精品18| 免费看十八禁软件| 久久伊人香网站| 亚洲精品中文字幕在线视频| 国产一区二区三区在线臀色熟女| 国产单亲对白刺激| av免费在线观看网站| 免费在线观看完整版高清| 精品免费久久久久久久清纯| 日本三级黄在线观看| 成人国产综合亚洲| 色综合亚洲欧美另类图片| 香蕉国产在线看| 亚洲欧美日韩高清在线视频| 好男人电影高清在线观看| 中文字幕另类日韩欧美亚洲嫩草| 99久久综合精品五月天人人| 欧美激情久久久久久爽电影| 美女免费视频网站| 在线观看免费日韩欧美大片| 国产精品一区二区精品视频观看| 精品国产乱码久久久久久男人| 波多野结衣巨乳人妻| av有码第一页| 久久 成人 亚洲| 亚洲av美国av| 最新美女视频免费是黄的| 久久精品国产亚洲av香蕉五月| 少妇的丰满在线观看| 男女午夜视频在线观看| 很黄的视频免费| 久久欧美精品欧美久久欧美| 麻豆国产av国片精品| 亚洲国产日韩欧美精品在线观看 | 在线观看免费视频日本深夜| 国产一卡二卡三卡精品| 99在线视频只有这里精品首页| 色综合婷婷激情| 99re在线观看精品视频| 俺也久久电影网| 国产久久久一区二区三区| 麻豆成人av在线观看| 国产av不卡久久| 制服诱惑二区| 色综合婷婷激情| av欧美777| 最近在线观看免费完整版| 精品高清国产在线一区| 村上凉子中文字幕在线| 久久国产精品男人的天堂亚洲| tocl精华| 在线十欧美十亚洲十日本专区| 亚洲久久久国产精品| 欧美一级毛片孕妇| 色播在线永久视频| av在线播放免费不卡| 国产精品美女特级片免费视频播放器 | 夜夜夜夜夜久久久久| 国产伦一二天堂av在线观看| 精品国产一区二区三区四区第35| 成人免费观看视频高清| 国产精品久久视频播放| 我的亚洲天堂| 日本黄色视频三级网站网址| 88av欧美| 给我免费播放毛片高清在线观看| 搡老熟女国产l中国老女人| 久久狼人影院| 99在线视频只有这里精品首页| 正在播放国产对白刺激| 成人av一区二区三区在线看| 嫩草影院精品99| 欧洲精品卡2卡3卡4卡5卡区| 一区二区三区高清视频在线| 亚洲精品粉嫩美女一区| 精品久久蜜臀av无| 亚洲人成网站在线播放欧美日韩| 成在线人永久免费视频| 高潮久久久久久久久久久不卡| av欧美777|