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

    基于小波變換的橋梁動(dòng)態(tài)稱重系統(tǒng)車軸高精度識(shí)別研究

    2017-05-09 06:14:42趙華譚承君張龍威喬?hào)|欽
    關(guān)鍵詞:小波變換

    趙華 譚承君 張龍威 喬?hào)|欽

    摘要:首先利用小波變換對(duì)一不能明顯識(shí)別車軸信息的數(shù)值仿真信號(hào)進(jìn)行處理,證明小波變換能夠高效放大車軸經(jīng)過傳感器時(shí)產(chǎn)生的不連續(xù)變化斜率,從而識(shí)別出車軸信息。然后基于實(shí)橋測(cè)試,對(duì)那些不能直接識(shí)別出車輛信息的FAD信號(hào),通過聯(lián)合控制最小Shan-non熵值和最大相關(guān)系數(shù)選取最適變換尺度和最適變換小波函數(shù)進(jìn)行小波變換。分析結(jié)果表明:對(duì)于不能直接識(shí)別出車輛信息的FAD信號(hào),小波變換也能準(zhǔn)確地識(shí)別車輛行駛速度、車軸數(shù)目以及車軸間距。小波變換可提高橋梁動(dòng)態(tài)稱重(BWIM)系統(tǒng)車軸識(shí)別的效率及精度,為將BWIM系統(tǒng)發(fā)展為超載車輛控制的有效工具提供技術(shù)支撐。

    關(guān)鍵詞:橋梁動(dòng)態(tài)稱重;車軸識(shí)別;小波變換;小波函數(shù)選?。蛔儞Q尺度

    中圖分類號(hào):U491;TN911.7 文獻(xiàn)標(biāo)識(shí)碼:A

    我國(guó)超載現(xiàn)象比比皆是,超載車輛對(duì)道路和橋梁造成的破壞不容忽視。近些年來,超載車輛導(dǎo)致橋梁垮塌的報(bào)道屢見不鮮。高效率、高精度的超載車輛監(jiān)管系統(tǒng)可控制超載車輛并監(jiān)測(cè)橋梁安全,在我國(guó)有著廣闊的應(yīng)用前景。

    動(dòng)態(tài)稱重系統(tǒng)(weigh-in-motion:WIM)近年來在國(guó)際上越來越多地應(yīng)用于對(duì)橋梁和道路交通車輛的監(jiān)測(cè)。該技術(shù)在測(cè)試車輛數(shù)據(jù)時(shí)不需要中斷交通,效率較高,是一種獲取車輛荷載數(shù)據(jù)、控制超重車輛的有效工具。傳統(tǒng)的動(dòng)態(tài)稱重系統(tǒng)(也稱為路面動(dòng)態(tài)稱重系統(tǒng):Pavement WIM)通過在路面埋置傳感器來測(cè)試通行車輛的軸重、軸距和車速等信息。該系統(tǒng)是永久性的,不能拆遷,使用壽命較短,且安裝及維護(hù)費(fèi)用較高,精度也有待提高。橋梁動(dòng)態(tài)稱重系統(tǒng)(bridge weigh-in-motion:BWIM)直接利用橋梁為載體,在線監(jiān)測(cè)橋梁在移動(dòng)車輛荷載作用下的動(dòng)態(tài)響應(yīng),并通過程序反算出車輛軸重。和傳統(tǒng)路面WlM系統(tǒng)相比,BWIM系統(tǒng)不僅可以在不阻斷交通的情況下連續(xù)不間斷地識(shí)別車輛軸重和總重,而且可以在線監(jiān)測(cè)橋梁在車輛荷載作用下的動(dòng)態(tài)響應(yīng)信號(hào),進(jìn)而獲得結(jié)構(gòu)的實(shí)際影響線以及在車輛荷載作用下的結(jié)構(gòu)荷載橫向分布參數(shù)和橋梁沖擊系數(shù)等。

    BWIM的概念最初由Moses于1979年提出。BWIM系統(tǒng)利用安裝在主梁下緣的傳感器所采集的動(dòng)態(tài)信號(hào)反算出車輛軸重。早期BWIM系統(tǒng)除了在主梁下緣安裝傳感器外,還需要在橋頭路面上埋置便攜式磁帶開關(guān)或壓電電纜來測(cè)試車輛速度、軸數(shù)、軸間距等信息,所以早期BWIM系統(tǒng)也存在耐久性問題,且安裝時(shí)需要中斷交通,更換也不易。

    近幾十年來,各國(guó)學(xué)者在原始BWIM模型基礎(chǔ)上推廣和改進(jìn),共同創(chuàng)建了新型商業(yè)BWIM系統(tǒng)。

    新型BWIM系統(tǒng)摒棄了早期BWIM系統(tǒng)所需的埋置在橋頭路面上的便攜式磁帶開關(guān)或壓電電纜,而是直接在橋面板下安裝額外的傳感器(也稱車軸探測(cè)傳感器或FAD傳感器),從而獲得車軸信息。新型BWIM系統(tǒng)的安裝、調(diào)試及數(shù)據(jù)采集全部在橋下進(jìn)行,橋面上不布置任何測(cè)試設(shè)備。相較于早期BWIM系統(tǒng),新型BWIM系統(tǒng)不僅可以在不阻斷交通的情況下連續(xù)不間斷地識(shí)別車輛軸重和總重,更具有攜帶方便、安裝和測(cè)試隱蔽、可以反復(fù)使用,且安裝、維護(hù)及使用成本低等優(yōu)點(diǎn)。

    Moses,Znidaric和Peters都曾提到車輛信息識(shí)別(軸數(shù)及軸間距)的準(zhǔn)確度是影響B(tài)WIM系統(tǒng)軸重及總重識(shí)別精度的一個(gè)主要因素,同時(shí)有效且準(zhǔn)確地識(shí)別出車軸信息也是BWIM系統(tǒng)精確識(shí)別車輛軸重的前提及基礎(chǔ)。然而試驗(yàn)證明,有時(shí)難以直接從FAD信號(hào)中準(zhǔn)確識(shí)別車輛信息,特別是在連軸(相鄰車軸間距很?。┗蛐盘?hào)動(dòng)態(tài)成分較大的情況下。Dunne和Chatterjee提出通過對(duì)FAD信號(hào)進(jìn)行小波變換可更有效地識(shí)別車輛信息,但對(duì)于小波函數(shù)選取及最適變換尺度都缺少研究。

    本文基于仿真分析及實(shí)橋測(cè)試結(jié)果(湖南省懷化市舞水五橋引橋的現(xiàn)場(chǎng)測(cè)試),在標(biāo)定車輛經(jīng)過測(cè)試橋梁時(shí),采集記錄FAD信號(hào)和橋梁橋底動(dòng)態(tài)響應(yīng)信號(hào),結(jié)合最小Shannon熵值以及最大相關(guān)性從而選取合適小波函數(shù)以及變換尺度,并利用小波技術(shù)變換FAD信號(hào)來獲得未能直接識(shí)別的車輛信息。分析結(jié)果表明:作為強(qiáng)大的信號(hào)處理工具,小波變換在BWIM系統(tǒng)中能夠有效提高對(duì)車軸信息的識(shí)別。

    1小波理論

    1982年法國(guó)工程師Jean Morlet首先提出小波變換的概念,它可以認(rèn)為是經(jīng)典傅里葉變換的延伸。兩者最大的區(qū)別在于小波變換可以同時(shí)在時(shí)域和頻域上定域,而傅里葉變換只是時(shí)域和頻域之間的轉(zhuǎn)換工具。小波變換可分為連續(xù)小波變換(Con-tinuous wavelet transform:CWT)和離散小波變換(Discrete wavelet transform:DWT)。

    1.1連續(xù)小波變換

    連續(xù)小波變換的數(shù)學(xué)定義為:

    (1)

    小波一詞意味著特定的小波函數(shù),即式(1)中的ψ(t),例如db2和rbi02。4(圖1)等。這種小波函數(shù)都是緊支的,即函數(shù)的定義域是有限的。傅里葉變換是將信號(hào)分解成一系列不同頻率的正弦波疊加,而小波變換是通過縮放和平移這些小波函數(shù)來逼近信號(hào)。顯然對(duì)于一些尖銳變化的信號(hào),用不規(guī)則的小波函數(shù)逼近要比光滑的正弦曲線好。

    經(jīng)過連續(xù)小波變換得到許多小波系數(shù),這些系數(shù)就是縮放因子(尺度)和平移(位置)的函數(shù),其中小波分解尺度與傅里葉變換中的頻率相對(duì)應(yīng),所以連續(xù)小波變換又可以定義為式(2)。它表示信號(hào)f(x)與被縮放和平移的小波函數(shù)ψ(t)之積在信號(hào)存在的整個(gè)時(shí)間段求和的結(jié)果。

    (2)

    1.2最適小波函數(shù)和變換尺度的選取

    在MATLAB小波分析工具中有13個(gè)小波函數(shù)族,包括60多個(gè)常用小波函數(shù),而對(duì)于同一個(gè)目標(biāo)信號(hào),每一個(gè)小波函數(shù)通過變換都會(huì)得到不一樣的結(jié)果。小波變換最大的挑戰(zhàn)是對(duì)于不同情況怎么選取小波函數(shù)來獲得最理想的結(jié)果。目前,對(duì)于小波函數(shù)選取的方法總體可以分為兩類:定性方法和定量方法。

    定性方法是基于小波函數(shù)的本身基本特性(正交性、對(duì)稱性、緊支性、光滑性、消失矩)或通過目標(biāo)信號(hào)與小波函數(shù)之間的形狀匹配度來選取最適小波函數(shù)。對(duì)于特定的目標(biāo)信號(hào),需要針對(duì)想要結(jié)果的特性選取合適的小波函數(shù)。例如,考慮到小波函數(shù)的緊支性、消失矩和正交性,選擇Coil4小波函數(shù)以便有效分離肌動(dòng)電流圖中的脈沖波和滋補(bǔ)成分。

    定性方法中大部分是通過肉眼來判斷,缺少客觀依據(jù),而定量方法則是通過計(jì)算出特定的數(shù)值(能量、熵值、相關(guān)系數(shù)等)以便更客觀地對(duì)比小波函數(shù)之間的差異。

    對(duì)于連續(xù)小波變換,它在每個(gè)可能的尺度上都進(jìn)行計(jì)算,得到一系列的小波變換系數(shù),其中只有部分結(jié)果是有用的,所以如何選取變換尺度至關(guān)重要。本文提出采用最小Shannon熵方法選擇最佳小波變換尺度。

    小波變換系數(shù)在尺度s下的能量為:

    (3)式中:ωt(s,i)為小波變換系數(shù);N為小波變換系數(shù)總?cè)萘俊?/p>

    小波變換系數(shù)的Shannon熵的分布如下:

    (4)

    顯然小波變換系數(shù)Shannon熵值的范圍為0≤Eentropy(s)≤log2N。僅當(dāng)所有小波變換系數(shù)都為0時(shí),Shannon熵值為0;且僅當(dāng)所有小波變換系數(shù)均相同時(shí),Shannon熵值等于log2 N。Shannon熵值的大小反映了小波變換系數(shù)概率分布的均勻性,熵值越大,概率分布越不均勻,能量就越分散;熵值越小,概率分布越均勻,能量分布就越集中。因此,對(duì)于不同尺度下一系列的小波變換系數(shù),其中某尺度對(duì)應(yīng)的小波變換系數(shù)包含少量大幅度的系數(shù)且其他的小波系數(shù)相對(duì)而言較小或者可以忽略時(shí),將產(chǎn)生最小的Shannon熵值。這一特性正適合在BWIM系統(tǒng)中利用FAD信號(hào)對(duì)車軸信息進(jìn)行識(shí)別。因?yàn)閷?duì)于FAD信號(hào),只需考慮當(dāng)車軸經(jīng)過FAD傳感器時(shí)產(chǎn)生的明顯波峰,其他時(shí)刻的振蕩均需忽略,使得車軸產(chǎn)生的峰值更加凸顯(即能量集中分布在車軸經(jīng)過傳感器時(shí)刻)。

    2數(shù)值仿真結(jié)果的小波變換

    仿真分析中模擬一輛兩軸車以20 m/s的速度經(jīng)過跨度為15 m的簡(jiǎn)支梁。車輛模型的簡(jiǎn)化車身質(zhì)量平均作用在車軸上(如圖2所示),且每個(gè)車輪自身重為7.5 kN。表1列出了計(jì)算分析中車輛和橋梁采用的詳細(xì)參數(shù)。

    圖3示出了車輛經(jīng)過橋梁時(shí)數(shù)值模擬產(chǎn)生的跨中梁底動(dòng)態(tài)應(yīng)變響應(yīng)信號(hào)(采樣頻率為6 400 Hz)。由于車橋耦合振動(dòng),很難直接從圖中獲得車軸經(jīng)過橋梁跨中時(shí)刻的準(zhǔn)確時(shí)間信息。本文擬用rbio2.4小波函數(shù)對(duì)圖3中模擬信號(hào)進(jìn)行連續(xù)小波變換,并為了選取最適合變換尺度,算出利用rbio2.4函數(shù)的連續(xù)小波變換在尺度1至64的Shannon熵值(圖4)。顯然,當(dāng)尺度為1時(shí),Shannon熵值最小,所以選取當(dāng)尺度為1時(shí)的連續(xù)小波變換系數(shù)作為變換結(jié)果(圖5)。其結(jié)果顯示的2個(gè)明顯峰值,表明所經(jīng)過車輛為一兩軸車。2個(gè)峰值所對(duì)應(yīng)的時(shí)刻分別為0.375 s和0.450 s,因此2個(gè)車軸經(jīng)過跨中點(diǎn)的時(shí)間差即為0.075 s。已知車輛行駛速度為20 m/s,即可算得峰值間距離(即軸距)為1.5 m,與車軸間距模型的設(shè)計(jì)值完全吻合。由此表明經(jīng)過小波變換后得到的波峰對(duì)應(yīng)時(shí)間點(diǎn)即為車軸經(jīng)過跨中點(diǎn)的瞬間。

    以上分析表明:對(duì)于沒有出現(xiàn)任何尖峰的跨中梁底車橋耦合振動(dòng)信號(hào),小波變換能非常有效地從這些信號(hào)中識(shí)別出各車軸經(jīng)過傳感器的準(zhǔn)確時(shí)間信息。其原因在于:這個(gè)數(shù)值模擬信號(hào)是由具有非?!凹饨恰钡撵o態(tài)信號(hào)疊加在振動(dòng)效果上所得到的。而小波變換就像是信號(hào)的放大器,能夠找到并放大這些“尖角”。

    3實(shí)橋測(cè)試

    3.1試驗(yàn)橋跨簡(jiǎn)介

    實(shí)橋測(cè)試選擇湖南省懷化市舞水五橋懷黔路側(cè)簡(jiǎn)支跨引橋作為測(cè)試跨。該橋主橋?yàn)轭A(yù)應(yīng)力混凝土連續(xù)箱梁橋,主橋跨徑布置為65 m+100 m+65 m。主橋兩側(cè)各有一跨徑為40 m引橋,為全預(yù)應(yīng)力混凝土簡(jiǎn)支T梁結(jié)構(gòu),其橫向由10片T梁組成。測(cè)試跨的傳感器布置如圖6,圖7和圖8所示。40 m引橋沿行車方向共布置5道橫隔板,如圖8所示。該橋橫向?qū)挾葹?4.0 m,橋?qū)挷贾脼椋?.0 m(人行道)+3.0 m(非機(jī)動(dòng)車道)+4×3.5 m(行車道)+3.0 m(非機(jī)動(dòng)車道)+2.0 m(人行道),具體布置如圖7所示。

    由圖6,圖7和圖8可以看到,在每個(gè)車道對(duì)應(yīng)的橋面板下都安裝了一對(duì)FAD傳感器,在每片T梁下緣都安裝了一個(gè)采集橋梁動(dòng)態(tài)應(yīng)變的傳感器(稱重傳感器)。如圖8所示,編號(hào)1-10為布置在梁底的稱重傳感器,編號(hào)11-18分別為布置在各個(gè)車道橋面板下的FAD傳感器。其中編號(hào)為11和15的傳感器對(duì)應(yīng)于行車道1,編號(hào)分別為L(zhǎng)1-FAD1和L1-FAD2;編號(hào)為12和16的傳感器對(duì)應(yīng)于行車道2,編號(hào)分別為L(zhǎng)2-FAD1和L2-FAD2;編號(hào)為13和17的傳感器對(duì)應(yīng)于行車道3,編號(hào)分別為L(zhǎng)3-FAD1和L3-FAD2;編號(hào)為14和18的傳感器對(duì)應(yīng)于行車道4,編號(hào)分別為L(zhǎng)4-FAD1和L4-FAD2。對(duì)各車道而言,F(xiàn)AD1和FAD2分別表示車輛先經(jīng)過編號(hào)為FAD1的傳感器,然后再經(jīng)過編號(hào)為FAD2的傳感器。

    在新型BWIM系統(tǒng)中,對(duì)測(cè)試的橋跨結(jié)構(gòu),要事先用標(biāo)定車對(duì)該橋跨進(jìn)行標(biāo)定試驗(yàn)。在整個(gè)標(biāo)定試驗(yàn)中采用標(biāo)定車(已知軸重和軸距)重復(fù)進(jìn)行跑車試驗(yàn),然后根據(jù)測(cè)試得到的動(dòng)態(tài)響應(yīng)來修正或計(jì)算結(jié)構(gòu)的實(shí)際影響線。在標(biāo)定后得到的影響線基礎(chǔ)上,以FAD傳感器識(shí)別的車輛信息(軸數(shù)、軸距等)為前提,根據(jù)連續(xù)采集的梁底動(dòng)態(tài)應(yīng)變計(jì)算出行駛車輛軸重及總重。

    試驗(yàn)現(xiàn)場(chǎng)傳感器的布置圖如圖9所示。試驗(yàn)時(shí)選取2輛兩軸車(分別稱為A和B)作為標(biāo)定車,同時(shí)或分別在每個(gè)車道上重復(fù)進(jìn)行跑車試驗(yàn)(如圖10所示)。試驗(yàn)中所采用的標(biāo)定車詳細(xì)信息見表2。試驗(yàn)過程中采用揚(yáng)州科動(dòng)公司生產(chǎn)的KD4001工具式應(yīng)變傳感器作為測(cè)試元件,日本TML公司生產(chǎn)的動(dòng)態(tài)應(yīng)變儀DC-204R作為數(shù)據(jù)采集儀。標(biāo)定試驗(yàn)采集頻率為200 Hz。

    3.2試驗(yàn)結(jié)果的小波變換

    3.2.1單輛車過橋時(shí)FAD信號(hào)的小波變換

    在試驗(yàn)中,當(dāng)標(biāo)定車輛行駛過橋時(shí),采集的FAD信號(hào)大部分能夠直接識(shí)別車輛信息,但有時(shí)由于車橋耦合振動(dòng)的動(dòng)態(tài)成分較大而導(dǎo)致FAD信號(hào)難以直接識(shí)別車輛信息。例如僅車輛A在車道3跑車時(shí),如圖8所示,編號(hào)為13的傳感器(L3-FAD1)和編號(hào)為17的傳感器(L3-FAD2)記錄了橋面板下緣處在車輛荷載作用下的動(dòng)態(tài)應(yīng)變信號(hào)(圖11)。從圖中可清晰地看出,標(biāo)定車經(jīng)過車道3時(shí),先經(jīng)過L3-FAD1,然后再經(jīng)過L3-FAD2,這和傳感器的布置一致。

    如圖11所示,L3-FAD1和L3-FAD2的信號(hào)都只出現(xiàn)了一個(gè)明顯的峰值。當(dāng)標(biāo)定車A的2個(gè)軸分別經(jīng)過每一個(gè)FAD傳感器時(shí),該傳感器采集的測(cè)試信號(hào)理應(yīng)會(huì)顯示出2個(gè)明顯的峰值,但是L3-FAD2的信號(hào)只顯示了一個(gè)明顯峰值信號(hào)以及在它之前的一個(gè)小波峰,而L3-FAD1的信號(hào)只呈現(xiàn)了一個(gè)明顯峰值信號(hào)及它前后許多小波峰。就L3-FAD1信號(hào)而言,程序難以直接識(shí)別出哪些波峰是由于車軸經(jīng)過該FAD傳感器時(shí)所產(chǎn)生的。

    小波變換擬用來提高FAD信號(hào)的車軸信息識(shí)別。為了更好地說明小波函數(shù)以及變換尺度的選取,以Re-verseBior(rbio)小波函數(shù)族中rbio1.3,rbio2.2,rbio2.4,rbio2.8,rbio3.3,rbio3.9,rbio4.4和rbio6.8為例。對(duì)于L3-FAD1信號(hào),由公式(3)和(4)算出上述小波函數(shù)從尺度1至64的連續(xù)小波變換的Shannon熵值分布(圖12)。由圖可知,每個(gè)小波函數(shù)所對(duì)應(yīng)小波系數(shù)的Shannon熵值曲線在尺度5至15存在最小值。以rbio2.2函數(shù)為例,當(dāng)尺度為6時(shí),Shannon熵達(dá)到最小值。圖13列出了rbio2.2在不同尺度下的連續(xù)小波變換系數(shù)??梢钥闯?,當(dāng)尺度較小時(shí),小波變換系數(shù)雜亂無序,對(duì)車軸信息識(shí)別毫無意義;當(dāng)尺度較大時(shí),小波變換系數(shù)的波峰將被稀釋(包括車軸產(chǎn)生的波峰),這也導(dǎo)致無法精確識(shí)別車軸信息。當(dāng)Shannon熵值最小時(shí)(尺度為6)呈現(xiàn)出2個(gè)明顯能識(shí)別的波峰,在此能量集中分布于2個(gè)波峰所對(duì)應(yīng)的時(shí)刻。并且當(dāng)尺度與最小Shannon熵值對(duì)應(yīng)的尺度越接近時(shí),由車軸經(jīng)過FAD傳感器產(chǎn)生的波峰越明顯。因此,也證明使Shannon熵值最小的小波變換尺度對(duì)應(yīng)的小波系數(shù)能獲得最佳車軸識(shí)別效果。

    前述已證明在最小Shannon熵時(shí),能夠得到最佳車軸信息識(shí)別效果,圖14列出上述一系列小波函數(shù)對(duì)應(yīng)最小熵值尺度下的連續(xù)小波變換系數(shù)。對(duì)于每一個(gè)小波函數(shù),盡管都處于最優(yōu)車軸識(shí)別效果狀態(tài),但不一定能夠有效地識(shí)別出車軸信息。從圖14中可以看出,rbio1.3,rbio3.3和rbio3.9的小波變換結(jié)果在車軸經(jīng)過FAD傳感器的時(shí)刻產(chǎn)生峰值,但由于波峰本身是反對(duì)稱形狀,導(dǎo)致波峰峰值減小,或者如rbio3.9結(jié)果所示,一個(gè)車軸經(jīng)過FAD時(shí)顯現(xiàn)出2個(gè)波峰,這都將導(dǎo)致車軸信息難以準(zhǔn)確識(shí)別。另外如rbio2.8和rbio6.8小波變換結(jié)果所示,在由車輛后軸產(chǎn)生的大波峰兩側(cè)會(huì)出現(xiàn)小波峰,峰值有時(shí)甚至大于由車軸產(chǎn)生的小波峰,這也將影響對(duì)車軸信息的判斷。

    事實(shí)上,只有rbio2.2,rbio2.4和rbio4.4小波變換系數(shù)顯示了2個(gè)可識(shí)別的波峰且其波峰形狀與目標(biāo)信號(hào)相似。本文提出通過計(jì)算相關(guān)系數(shù)來進(jìn)一步選取小波函數(shù)。相關(guān)系數(shù)計(jì)算公式見式(5)。

    (5)式中:x為目標(biāo)信號(hào);x為目標(biāo)信號(hào)平均值;y為所得小波變換系數(shù);y為小波變換系數(shù)的平均值。

    由公式(5)計(jì)算得到圖14中的小波變換結(jié)果與L3-FAD1信號(hào)的相關(guān)系數(shù)(圖15)。識(shí)別效果較差的rbio1.3,rbio3.3和rbio3.9小波函數(shù)所對(duì)應(yīng)的相關(guān)系數(shù)很小,不超過0.04;而rbio2.8和rbio6.8小波函數(shù)對(duì)應(yīng)的相關(guān)系數(shù)稍大,但也不超過0.13;車軸識(shí)別效果較好的rbio2.2,rbio2.4和rbio4.4小波函數(shù)所對(duì)應(yīng)的相關(guān)系數(shù)均超過0.25,且rbio2.2對(duì)應(yīng)的相關(guān)系數(shù)達(dá)到0.50。

    因此,利用小波變換來高效識(shí)別FAD車軸信息,可以通過最大相關(guān)性和最小Shannon熵值來選取最適合的小波函數(shù)及其最佳變換尺度。針對(duì)MATLAB中所有常用小波函數(shù),基于上述方法對(duì)圖11中的FAD信號(hào)分別進(jìn)行小波分析。對(duì)于L3-FAD1信號(hào),在rbio2.2且尺度為6時(shí)獲得最大相關(guān)系數(shù)以及最小熵值,而L3-FAD2信號(hào),在mexh小波函數(shù)尺度為3時(shí)獲得最大相關(guān)系數(shù)以及最小熵值。所對(duì)應(yīng)連續(xù)小波變換結(jié)果如圖16所示。這些原始信號(hào)經(jīng)過連續(xù)小波變換后,可以很明顯地識(shí)別出車軸數(shù)目,為BWIM系統(tǒng)識(shí)別車輛軸重及總重提供前提條件。

    當(dāng)車軸數(shù)目及對(duì)應(yīng)時(shí)刻被一對(duì)FAD傳感器識(shí)別后,可計(jì)算出車輛行駛速度和車軸間距,同時(shí)還可與靜態(tài)測(cè)量的軸距進(jìn)行比較。從圖16可以得到,標(biāo)定車A經(jīng)過L3-FAD1時(shí),前后2個(gè)波峰對(duì)應(yīng)的時(shí)刻分別為t1=2.115 s和t2=2.730 s,經(jīng)過L3-FAD2時(shí)2個(gè)峰值對(duì)應(yīng)的時(shí)刻分別為t'1=3.150 s和t'2=3.740 s。在已知2個(gè)FAD傳感器距離(如圖6和圖8所示,LFAD=8 m)前提下,車輛行駛速度可計(jì)算如下:

    (6)

    在求出速度V后,根據(jù)任何一個(gè)FAD信號(hào)中2個(gè)波峰之間的時(shí)間差即可計(jì)算出軸距,計(jì)算公式如下:

    A1,2=V(t'2-t'1)或者A1,2=V(t2-t1)。 (7)通過公式(6)和(7)可得到車輛的軸間距并和相應(yīng)的靜態(tài)測(cè)量結(jié)果比較,其結(jié)果均列在表3中。

    3.2.2多輛車同時(shí)過橋時(shí)FAD信號(hào)的小波變換

    為了進(jìn)一步驗(yàn)證小波變換應(yīng)用于BWIM系統(tǒng)中提高車軸識(shí)別效果的有效性,下面利用小波變換識(shí)別當(dāng)兩輛車同時(shí)行駛過橋時(shí)所產(chǎn)生的不易直接獲得的車輛信息的FAD信號(hào)。

    標(biāo)定車A和B分別同時(shí)沿車道4和車道3行駛過橋,對(duì)應(yīng)于這2個(gè)車道的FAD傳感器L3-FAD1,L3-FAD2,L4-FAD1和L4-FAD2采集的動(dòng)態(tài)應(yīng)變信號(hào)如圖13所示。由于試驗(yàn)持續(xù)時(shí)間長(zhǎng),傳感器L3-FAD2工作不穩(wěn)定,產(chǎn)生了白噪聲干擾信號(hào)(如圖17(a)所示)。

    如圖17所示,圖中L3-FAD1,L3-FAD2,L4-FAD1和L4-FAD2的信號(hào)都只顯示出一個(gè)峰值信號(hào)和一些小波峰,通過這些信號(hào)難以直接識(shí)別過橋車輛車軸信息,尤其是信號(hào)L3-FAD1和L4-FAD2,這種現(xiàn)象更明顯。運(yùn)用上述方法從常用小波函數(shù)里面找出最適小波函數(shù)分別對(duì)2個(gè)車道的FAD信號(hào)進(jìn)行連續(xù)小波變換,得到的結(jié)果如圖18所示。盡管信號(hào)L3-FAD2摻雜有較大的白噪音信號(hào),在coif1函數(shù)尺度14的小波變換下仍能高效識(shí)別出車軸信息。另外L3-FAD1,L4-FAD1和L4-FAD2信號(hào)分別在rbio2.2尺度為6,gaus2尺度為4和rbio2.2尺度為7時(shí)獲得理想結(jié)果(如圖18所示)。

    這些FAD信號(hào)在經(jīng)過連續(xù)小波變換后都表現(xiàn)出2個(gè)明顯的尖峰。對(duì)于車道3或車道4,從圖18中可明確得到車輛各車軸經(jīng)過每個(gè)傳感器的時(shí)刻。同樣,根據(jù)相對(duì)應(yīng)的一對(duì)FAD傳感器,由公式(6)和(7)可計(jì)算出車輛A和B的軸距,結(jié)果列于表4。

    4結(jié)論

    1)從數(shù)值模擬信號(hào)的小波分析結(jié)果可以看出,小波變換能夠高效率地發(fā)現(xiàn)和放大車軸經(jīng)過傳感器時(shí)產(chǎn)生的不連續(xù)變化斜率,從而準(zhǔn)確識(shí)別出過橋車輛車軸信息。

    2)實(shí)橋試驗(yàn)結(jié)果分析表明,小波變換能有效解決新型橋梁動(dòng)態(tài)稱重(BwIM)系統(tǒng)中部分FAD信號(hào)難以直接識(shí)別車軸數(shù)目、車軸距離及車輛速度的問題。在2組分析例子中,對(duì)于所采集的不理想FAD信號(hào),由小波變換所識(shí)別的車輛軸距誤差在4%以內(nèi)。

    3)在利用小波變換識(shí)別車軸信息時(shí),小波變換尺度至關(guān)重要,低尺度小波系數(shù)顯示高頻信息,而高尺度小波系數(shù)顯示低頻信息。根據(jù)最小shannon熵的方法得到,只有在使shannon熵值最小時(shí)才能獲得最利于識(shí)別車軸信息的小波變換尺度。

    4)通過對(duì)常用的幾十種小波函數(shù)對(duì)同一目標(biāo)信號(hào)進(jìn)行變換,所得的結(jié)果不一致,并且部分小波函數(shù)并不能達(dá)到預(yù)期結(jié)果,而只有選取合適的小波函數(shù)才能準(zhǔn)確地識(shí)別車軸信息。依據(jù)本文所提出的采用最大相關(guān)系數(shù)的方法,可以獲取最合適的小波函數(shù),從而高效率地識(shí)別車軸信息。

    5)在新型BWIM系統(tǒng)中,對(duì)于那些不能直接識(shí)別出車軸信息的FAD信號(hào),小波變換是一個(gè)高效而準(zhǔn)確地識(shí)別車軸信息的工具。在應(yīng)用小波變換時(shí),首先基于最小Shannon熵值計(jì)算出每個(gè)小波函數(shù)的最適變換尺度,然后求出各最適尺度下的變換系數(shù)與目標(biāo)信號(hào)的相關(guān)系數(shù),具有最大相關(guān)系數(shù)者即為最適小波函數(shù)。試驗(yàn)證明,對(duì)于難以識(shí)別車軸信息的FAD信號(hào),依據(jù)最小shannon熵值方法和最大相關(guān)系數(shù)確定的小波變換系數(shù)也能準(zhǔn)確、有效地識(shí)別車軸信息,進(jìn)而提高新型BWIM系統(tǒng)識(shí)別過橋車輛軸重及總重量的效率和精度。

    猜你喜歡
    小波變換
    最優(yōu)小波包變換在齒輪箱振動(dòng)信號(hào)去噪中的應(yīng)用研究
    負(fù)壓波法管道泄漏監(jiān)測(cè)定位系統(tǒng)實(shí)踐與探索
    卷宗(2016年10期)2017-01-21 15:48:38
    基于峰度和小波變換的超短波信號(hào)調(diào)制識(shí)別
    基于雙樹四元數(shù)小波變換的圖像降噪增強(qiáng)
    MATLAB在《數(shù)字圖像處理》課程中的輔助教學(xué)
    基于互信息和小波變換的圖像配準(zhǔn)的研究
    保持細(xì)節(jié)的Retinex紅外圖像增強(qiáng)算法
    基于小波變換的數(shù)字水印技術(shù)
    考試周刊(2016年45期)2016-06-24 13:55:17
    基于Matlab的遙感圖像IHS小波融合算法的并行化設(shè)計(jì)
    科技視界(2016年11期)2016-05-23 08:13:35
    基于小波包絡(luò)分析的滾動(dòng)軸承典型故障診斷技術(shù)研究
    商情(2016年11期)2016-04-15 07:39:54
    高清毛片免费观看视频网站| 欧美中文日本在线观看视频| 亚洲成人精品中文字幕电影| 又黄又爽又免费观看的视频| 国产精品电影一区二区三区| 国产伦一二天堂av在线观看| 精品不卡国产一区二区三区| 日本成人三级电影网站| 级片在线观看| www.999成人在线观看| 91成人精品电影| 午夜福利高清视频| 成人av一区二区三区在线看| 欧美三级亚洲精品| 欧美三级亚洲精品| 每晚都被弄得嗷嗷叫到高潮| 少妇粗大呻吟视频| 欧美+亚洲+日韩+国产| 男女午夜视频在线观看| 国产精品久久久av美女十八| 亚洲国产精品999在线| 久久亚洲精品不卡| 久久中文看片网| 美女大奶头视频| 亚洲精品中文字幕在线视频| 美女 人体艺术 gogo| 国产极品粉嫩免费观看在线| 老司机午夜十八禁免费视频| 校园春色视频在线观看| 50天的宝宝边吃奶边哭怎么回事| 在线国产一区二区在线| netflix在线观看网站| 女人爽到高潮嗷嗷叫在线视频| 亚洲在线自拍视频| 可以在线观看毛片的网站| 可以在线观看毛片的网站| 欧美一级a爱片免费观看看 | 18禁观看日本| 国产免费男女视频| 99精品在免费线老司机午夜| 精品免费久久久久久久清纯| 午夜亚洲福利在线播放| 十八禁人妻一区二区| 亚洲精品久久国产高清桃花| 国产成人影院久久av| 18禁美女被吸乳视频| 国产欧美日韩一区二区三| 欧美成人午夜精品| 人人妻人人看人人澡| √禁漫天堂资源中文www| 久久香蕉激情| 亚洲九九香蕉| 国产私拍福利视频在线观看| 韩国精品一区二区三区| 美国免费a级毛片| 黄色a级毛片大全视频| 热99re8久久精品国产| 性色av乱码一区二区三区2| 午夜免费观看网址| cao死你这个sao货| 极品教师在线免费播放| 长腿黑丝高跟| 不卡一级毛片| 亚洲一区中文字幕在线| 亚洲在线自拍视频| 特大巨黑吊av在线直播 | 国产又黄又爽又无遮挡在线| 久久久久久久久免费视频了| 老司机福利观看| xxxwww97欧美| 欧美日韩一级在线毛片| 午夜免费激情av| 亚洲欧美一区二区三区黑人| 国产伦人伦偷精品视频| 嫩草影院精品99| www日本黄色视频网| 老汉色∧v一级毛片| 中出人妻视频一区二区| 国产一区二区激情短视频| 日韩欧美一区二区三区在线观看| 精品不卡国产一区二区三区| 99久久99久久久精品蜜桃| 看黄色毛片网站| 精品久久久久久成人av| 成人手机av| 国产精品99久久99久久久不卡| 亚洲男人的天堂狠狠| 波多野结衣av一区二区av| 国内精品久久久久精免费| 国产99久久九九免费精品| 国产精品电影一区二区三区| 观看免费一级毛片| 黄色女人牲交| 国产亚洲av高清不卡| 黑丝袜美女国产一区| 草草在线视频免费看| 天天躁狠狠躁夜夜躁狠狠躁| 免费无遮挡裸体视频| 久久久久九九精品影院| 亚洲一码二码三码区别大吗| 一个人观看的视频www高清免费观看 | 欧美乱码精品一区二区三区| 一区二区三区激情视频| 欧美精品亚洲一区二区| 精品第一国产精品| 18禁国产床啪视频网站| 亚洲国产日韩欧美精品在线观看 | 亚洲,欧美精品.| 成人国产一区最新在线观看| 99国产极品粉嫩在线观看| 999久久久精品免费观看国产| 黄网站色视频无遮挡免费观看| 日本撒尿小便嘘嘘汇集6| 男人操女人黄网站| 大型av网站在线播放| 久久国产精品男人的天堂亚洲| 日韩一卡2卡3卡4卡2021年| 免费在线观看完整版高清| 国产精品久久久av美女十八| 久久久久久免费高清国产稀缺| 国产精品98久久久久久宅男小说| 国产伦在线观看视频一区| 久久精品国产亚洲av高清一级| 岛国视频午夜一区免费看| 亚洲自拍偷在线| 国产精品久久久久久亚洲av鲁大| 一区二区三区激情视频| 可以在线观看的亚洲视频| 51午夜福利影视在线观看| 国产精品久久视频播放| 很黄的视频免费| 国产激情久久老熟女| 亚洲av成人一区二区三| 99久久久亚洲精品蜜臀av| 国产成人系列免费观看| 国产精品二区激情视频| 天天躁夜夜躁狠狠躁躁| 免费看a级黄色片| 国产野战对白在线观看| 国产又色又爽无遮挡免费看| 在线观看一区二区三区| 久久久国产精品麻豆| 久久人妻福利社区极品人妻图片| 欧美在线一区亚洲| cao死你这个sao货| 在线观看免费日韩欧美大片| 给我免费播放毛片高清在线观看| 97超级碰碰碰精品色视频在线观看| 亚洲欧美激情综合另类| 成年免费大片在线观看| 男女那种视频在线观看| 欧美不卡视频在线免费观看 | 国产激情久久老熟女| 97人妻精品一区二区三区麻豆 | 在线观看免费日韩欧美大片| 国产黄色小视频在线观看| 久久香蕉激情| 国产精品久久电影中文字幕| 久久 成人 亚洲| 日日爽夜夜爽网站| 午夜福利在线观看吧| 久久香蕉激情| 精品久久久久久久久久久久久 | 2021天堂中文幕一二区在线观 | 国内精品久久久久精免费| 国产精品98久久久久久宅男小说| 国产成人系列免费观看| 90打野战视频偷拍视频| 禁无遮挡网站| www.熟女人妻精品国产| 国产伦在线观看视频一区| 一进一出抽搐gif免费好疼| 久热这里只有精品99| 日本免费a在线| 18禁黄网站禁片午夜丰满| 50天的宝宝边吃奶边哭怎么回事| 亚洲三区欧美一区| 黄片小视频在线播放| 人成视频在线观看免费观看| 精品第一国产精品| 18禁观看日本| 国产高清有码在线观看视频 | 我的亚洲天堂| 黄色 视频免费看| av在线播放免费不卡| 欧美日本亚洲视频在线播放| √禁漫天堂资源中文www| 看免费av毛片| 九色国产91popny在线| 国产熟女午夜一区二区三区| 很黄的视频免费| 麻豆久久精品国产亚洲av| 亚洲中文字幕日韩| 波多野结衣av一区二区av| 露出奶头的视频| 后天国语完整版免费观看| 午夜久久久久精精品| 亚洲色图av天堂| 欧美av亚洲av综合av国产av| 欧美日韩中文字幕国产精品一区二区三区| 亚洲五月婷婷丁香| 国产aⅴ精品一区二区三区波| 欧美一区二区精品小视频在线| 免费观看人在逋| 亚洲第一电影网av| 麻豆久久精品国产亚洲av| 久久久精品欧美日韩精品| 人妻久久中文字幕网| 国产亚洲av高清不卡| 国产精品九九99| 亚洲免费av在线视频| 欧美精品亚洲一区二区| 成人免费观看视频高清| 人妻久久中文字幕网| 亚洲国产高清在线一区二区三 | 亚洲人成电影免费在线| 免费在线观看视频国产中文字幕亚洲| 免费av毛片视频| 99精品久久久久人妻精品| 国产成人精品久久二区二区免费| 午夜福利欧美成人| 在线观看免费午夜福利视频| 亚洲国产精品合色在线| 亚洲av成人av| 老司机午夜福利在线观看视频| 国产亚洲精品久久久久久毛片| 国产三级黄色录像| 亚洲av电影不卡..在线观看| 妹子高潮喷水视频| 中文字幕另类日韩欧美亚洲嫩草| 两个人看的免费小视频| 99国产极品粉嫩在线观看| 色播亚洲综合网| 我的亚洲天堂| 欧美在线黄色| 麻豆成人午夜福利视频| 亚洲最大成人中文| 国产午夜精品久久久久久| 国产成人欧美在线观看| 欧美三级亚洲精品| 91av网站免费观看| 国产aⅴ精品一区二区三区波| 国产av在哪里看| 一级片免费观看大全| 欧美最黄视频在线播放免费| 国产精品久久久av美女十八| 男人操女人黄网站| 午夜激情av网站| 久久香蕉国产精品| 婷婷丁香在线五月| 一区二区三区激情视频| 亚洲精品中文字幕一二三四区| 午夜成年电影在线免费观看| 自线自在国产av| 免费人成视频x8x8入口观看| 免费一级毛片在线播放高清视频| 国产精品久久久人人做人人爽| 亚洲人成电影免费在线| 色老头精品视频在线观看| 久久久水蜜桃国产精品网| 欧美中文日本在线观看视频| 国产一级毛片七仙女欲春2 | 在线观看66精品国产| 久久久国产精品麻豆| 精品久久久久久成人av| 激情在线观看视频在线高清| 美女午夜性视频免费| 国内毛片毛片毛片毛片毛片| 久久久国产成人精品二区| 午夜日韩欧美国产| 亚洲 国产 在线| 亚洲精品国产精品久久久不卡| 欧美成狂野欧美在线观看| 国产aⅴ精品一区二区三区波| 手机成人av网站| АⅤ资源中文在线天堂| 国产区一区二久久| 成人特级黄色片久久久久久久| 国产精品av久久久久免费| 国内精品久久久久精免费| 美女午夜性视频免费| 又黄又粗又硬又大视频| 免费高清视频大片| 亚洲第一电影网av| 啦啦啦免费观看视频1| 久久久国产成人免费| 90打野战视频偷拍视频| 国内毛片毛片毛片毛片毛片| 国产在线精品亚洲第一网站| 国产又爽黄色视频| 非洲黑人性xxxx精品又粗又长| 亚洲成人免费电影在线观看| 欧美激情高清一区二区三区| 可以在线观看毛片的网站| 精品国产美女av久久久久小说| 日韩一卡2卡3卡4卡2021年| 黄网站色视频无遮挡免费观看| 中文字幕另类日韩欧美亚洲嫩草| 中文字幕av电影在线播放| 宅男免费午夜| 人妻久久中文字幕网| 伊人久久大香线蕉亚洲五| 18美女黄网站色大片免费观看| 一本大道久久a久久精品| 天堂√8在线中文| 国产高清videossex| 亚洲国产毛片av蜜桃av| 在线观看免费午夜福利视频| 黑人巨大精品欧美一区二区mp4| 午夜免费观看网址| 一级作爱视频免费观看| 无人区码免费观看不卡| 久久久国产精品麻豆| 午夜福利视频1000在线观看| 亚洲久久久国产精品| 一二三四社区在线视频社区8| 久久亚洲真实| 国产精品久久电影中文字幕| 亚洲午夜精品一区,二区,三区| 国产蜜桃级精品一区二区三区| 青草久久国产| 国产1区2区3区精品| www.自偷自拍.com| 久久精品国产亚洲av高清一级| 国产精品 欧美亚洲| 国产午夜福利久久久久久| 亚洲专区中文字幕在线| 亚洲五月色婷婷综合| 国产成年人精品一区二区| 女生性感内裤真人,穿戴方法视频| 性欧美人与动物交配| 黑人巨大精品欧美一区二区mp4| 黄色视频不卡| 首页视频小说图片口味搜索| 久久青草综合色| 村上凉子中文字幕在线| 制服丝袜大香蕉在线| 真人一进一出gif抽搐免费| 黄色丝袜av网址大全| 国产一级毛片七仙女欲春2 | 热99re8久久精品国产| 亚洲天堂国产精品一区在线| 老司机靠b影院| 香蕉av资源在线| 禁无遮挡网站| 亚洲国产日韩欧美精品在线观看 | 欧美日韩乱码在线| 国产片内射在线| 国产爱豆传媒在线观看 | 亚洲五月婷婷丁香| 美国免费a级毛片| 久久久久久久久久黄片| 国产麻豆成人av免费视频| 又黄又爽又免费观看的视频| 长腿黑丝高跟| 99在线人妻在线中文字幕| 啪啪无遮挡十八禁网站| 日韩大尺度精品在线看网址| 国产真人三级小视频在线观看| 巨乳人妻的诱惑在线观看| 亚洲一卡2卡3卡4卡5卡精品中文| 无限看片的www在线观看| 宅男免费午夜| e午夜精品久久久久久久| 亚洲五月色婷婷综合| 亚洲人成77777在线视频| 亚洲avbb在线观看| www.自偷自拍.com| 母亲3免费完整高清在线观看| 国产1区2区3区精品| 色在线成人网| 狠狠狠狠99中文字幕| 观看免费一级毛片| 此物有八面人人有两片| 人成视频在线观看免费观看| 精品国产乱子伦一区二区三区| 亚洲自拍偷在线| 青草久久国产| 久久精品亚洲精品国产色婷小说| 一进一出好大好爽视频| av片东京热男人的天堂| 黑人巨大精品欧美一区二区mp4| 国产一区二区三区在线臀色熟女| √禁漫天堂资源中文www| 好男人电影高清在线观看| 母亲3免费完整高清在线观看| 真人做人爱边吃奶动态| 亚洲精品久久成人aⅴ小说| 亚洲欧美激情综合另类| 日本黄色视频三级网站网址| 国产高清有码在线观看视频 | 一本一本综合久久| 9191精品国产免费久久| 亚洲国产欧洲综合997久久, | 人人妻,人人澡人人爽秒播| 午夜福利在线观看吧| 久久久久久久久中文| 男女做爰动态图高潮gif福利片| 日本一区二区免费在线视频| 色老头精品视频在线观看| 国内揄拍国产精品人妻在线 | 久久久久国产精品人妻aⅴ院| 97人妻精品一区二区三区麻豆 | 国产在线精品亚洲第一网站| 午夜福利18| 欧美成狂野欧美在线观看| 又黄又粗又硬又大视频| www.999成人在线观看| 欧美另类亚洲清纯唯美| 男人舔奶头视频| 亚洲五月婷婷丁香| 日韩国内少妇激情av| 亚洲精品粉嫩美女一区| 午夜福利成人在线免费观看| 大型黄色视频在线免费观看| 日本 欧美在线| 看黄色毛片网站| 国产高清激情床上av| 精品熟女少妇八av免费久了| 欧美日韩瑟瑟在线播放| АⅤ资源中文在线天堂| 熟女少妇亚洲综合色aaa.| 亚洲精品一卡2卡三卡4卡5卡| 久久精品亚洲精品国产色婷小说| 午夜福利一区二区在线看| 18美女黄网站色大片免费观看| 婷婷丁香在线五月| 国产精品乱码一区二三区的特点| 成人三级黄色视频| 亚洲欧美一区二区三区黑人| 久久草成人影院| 欧美国产精品va在线观看不卡| 亚洲欧美日韩无卡精品| 日日干狠狠操夜夜爽| 成熟少妇高潮喷水视频| 岛国在线观看网站| 日韩 欧美 亚洲 中文字幕| cao死你这个sao货| 欧美大码av| 99久久99久久久精品蜜桃| 国产精品一区二区三区四区久久 | 制服人妻中文乱码| 亚洲欧美日韩无卡精品| 一级片免费观看大全| 黄色丝袜av网址大全| 成人三级做爰电影| 亚洲 国产 在线| 午夜福利视频1000在线观看| 亚洲国产精品合色在线| or卡值多少钱| 日韩欧美国产在线观看| 丁香欧美五月| www国产在线视频色| 动漫黄色视频在线观看| 国产成人精品无人区| 一级毛片高清免费大全| 亚洲男人天堂网一区| 国产视频内射| 一区二区三区激情视频| 久久久久久久精品吃奶| 丁香六月欧美| 欧美黄色片欧美黄色片| 亚洲男人的天堂狠狠| 日韩欧美国产一区二区入口| 亚洲成av人片免费观看| 中文字幕精品免费在线观看视频| www.999成人在线观看| 午夜福利成人在线免费观看| 此物有八面人人有两片| 免费在线观看黄色视频的| 精品国产超薄肉色丝袜足j| 男女视频在线观看网站免费 | 啦啦啦观看免费观看视频高清| 久久精品夜夜夜夜夜久久蜜豆 | 中文资源天堂在线| 亚洲国产精品999在线| 国产麻豆成人av免费视频| 午夜激情av网站| 欧美日韩亚洲综合一区二区三区_| 久久热在线av| 中文字幕人成人乱码亚洲影| 日本撒尿小便嘘嘘汇集6| 久久婷婷成人综合色麻豆| 18美女黄网站色大片免费观看| 国产伦一二天堂av在线观看| 亚洲中文av在线| xxx96com| 午夜福利一区二区在线看| 久久精品国产99精品国产亚洲性色| 老熟妇仑乱视频hdxx| 18禁裸乳无遮挡免费网站照片 | 久久久久亚洲av毛片大全| 香蕉国产在线看| 久热爱精品视频在线9| 午夜久久久久精精品| 天天躁狠狠躁夜夜躁狠狠躁| 大香蕉久久成人网| 国产精品1区2区在线观看.| 搡老妇女老女人老熟妇| 精品久久久久久久久久免费视频| 黄色女人牲交| 免费在线观看成人毛片| 欧美成人午夜精品| 日日干狠狠操夜夜爽| 亚洲国产中文字幕在线视频| 久久婷婷人人爽人人干人人爱| 久久国产精品人妻蜜桃| 91麻豆精品激情在线观看国产| 99久久综合精品五月天人人| 亚洲第一av免费看| 中出人妻视频一区二区| 亚洲片人在线观看| 久久伊人香网站| 免费看a级黄色片| 国产成人av教育| 麻豆久久精品国产亚洲av| 欧美人与性动交α欧美精品济南到| 国产亚洲av高清不卡| cao死你这个sao货| 岛国在线观看网站| 黄色视频,在线免费观看| 91在线观看av| 两个人视频免费观看高清| 国产精品1区2区在线观看.| 在线观看午夜福利视频| 国产精品1区2区在线观看.| 免费在线观看影片大全网站| 啦啦啦 在线观看视频| 黑人欧美特级aaaaaa片| 在线永久观看黄色视频| 亚洲avbb在线观看| av欧美777| 日韩精品免费视频一区二区三区| 亚洲精品色激情综合| 婷婷丁香在线五月| 国产99白浆流出| 亚洲一区中文字幕在线| 叶爱在线成人免费视频播放| 成人精品一区二区免费| 国产av不卡久久| 最新美女视频免费是黄的| 久久精品夜夜夜夜夜久久蜜豆 | 一级毛片女人18水好多| 国产精品久久久av美女十八| 国产成人精品久久二区二区91| 99久久国产精品久久久| 999久久久精品免费观看国产| 在线观看午夜福利视频| 变态另类成人亚洲欧美熟女| 久久精品国产亚洲av高清一级| 91九色精品人成在线观看| 久久久久久久午夜电影| 麻豆一二三区av精品| 久久久国产精品麻豆| 香蕉久久夜色| 日本免费a在线| 欧美日韩瑟瑟在线播放| cao死你这个sao货| 亚洲专区国产一区二区| 最近最新免费中文字幕在线| 亚洲专区国产一区二区| 国产精品 国内视频| 日韩欧美在线二视频| 国产一级毛片七仙女欲春2 | 国产精品亚洲美女久久久| 亚洲狠狠婷婷综合久久图片| 欧美激情久久久久久爽电影| 久久天躁狠狠躁夜夜2o2o| 欧美成人午夜精品| 欧美乱妇无乱码| 国产一级毛片七仙女欲春2 | 天天躁夜夜躁狠狠躁躁| 国产日本99.免费观看| 最好的美女福利视频网| 欧美成人免费av一区二区三区| 久久婷婷人人爽人人干人人爱| 人妻丰满熟妇av一区二区三区| 国产三级黄色录像| 搞女人的毛片| 久久天躁狠狠躁夜夜2o2o| 亚洲国产高清在线一区二区三 | 啦啦啦观看免费观看视频高清| 国产精品日韩av在线免费观看| 一本综合久久免费| 欧美日韩乱码在线| 国产精品乱码一区二三区的特点| 69av精品久久久久久| 黄网站色视频无遮挡免费观看| 国产成人av教育| 波多野结衣高清无吗| 国产成人一区二区三区免费视频网站| 91字幕亚洲| e午夜精品久久久久久久| 美女大奶头视频| √禁漫天堂资源中文www| 久久天堂一区二区三区四区| 怎么达到女性高潮| 亚洲av日韩精品久久久久久密| 一本综合久久免费| 日本免费a在线| 久久国产精品影院| 嫁个100分男人电影在线观看| 亚洲精品在线美女| 午夜精品久久久久久毛片777| 婷婷精品国产亚洲av在线| 少妇的丰满在线观看| 男女下面进入的视频免费午夜 | 妹子高潮喷水视频| x7x7x7水蜜桃| 18禁黄网站禁片免费观看直播| 高清在线国产一区| 国产av不卡久久| 搡老熟女国产l中国老女人| 一级片免费观看大全| 色哟哟哟哟哟哟|