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

    基于多頻率數(shù)據(jù)融合的Lamb波損傷定位研究

    2017-10-26 02:47:22楊曉華劉學(xué)君馬廣婷
    聲學(xué)技術(shù) 2017年2期
    關(guān)鍵詞:時(shí)頻換能器時(shí)域

    楊曉華,劉學(xué)君,馬廣婷

    ?

    基于多頻率數(shù)據(jù)融合的Lamb波損傷定位研究

    楊曉華1,劉學(xué)君1,馬廣婷2

    (1. 海軍航空工程學(xué)院青島校區(qū)航空機(jī)械系,山東青島266041;2. 91899部隊(duì),遼寧葫蘆島125001)

    確定損傷位置是結(jié)構(gòu)健康監(jiān)測(cè)中的重要環(huán)節(jié),而時(shí)間延遲是確定損傷位置的關(guān)鍵問題。針對(duì)Lamb波的頻散效應(yīng)對(duì)定位精度的影響,提出了一種基于Morlet小波時(shí)頻分析和多頻率數(shù)據(jù)融合的損傷定位方法,首先采用小波分析的手段,提取多組特定頻率下的信號(hào),實(shí)現(xiàn)多頻率成分下波達(dá)時(shí)刻的準(zhǔn)確判定,進(jìn)而采用加權(quán)平均處理的方法獲得損傷位置到兩換能器的距離之和,最后基于橢圓定位法確定了損傷的位置。試驗(yàn)結(jié)果表明,該方法能夠有效減少Lamb波頻散效應(yīng)的影響,相較時(shí)域峰值包絡(luò)的方法,可有效提高定位的精度,同時(shí)增強(qiáng)了系統(tǒng)抗干擾的能力。

    Lamb波;Morlet小波;橢圓定位法;數(shù)據(jù)融合;時(shí)間延遲

    0 引言

    近年來,基于導(dǎo)波的飛機(jī)結(jié)構(gòu)健康監(jiān)測(cè)技術(shù)引起了研究人員的廣泛關(guān)注。Lamb波是在板類結(jié)構(gòu)中傳播的導(dǎo)波,可以由壓電換能器、磁致伸縮換能器等激發(fā)和感知,具有易激發(fā)、傳播距離遠(yuǎn)及監(jiān)測(cè)范圍大等特點(diǎn),是對(duì)薄板類結(jié)構(gòu)如飛機(jī)蒙皮等進(jìn)行在線監(jiān)測(cè)的有效手段。結(jié)構(gòu)健康監(jiān)測(cè)的主要任務(wù)有:判斷是否存在損傷、定位損傷的位置、評(píng)估損傷的類型及程度、分析預(yù)測(cè)結(jié)構(gòu)壽命及評(píng)估結(jié)構(gòu)安全性等,而定位損傷的位置是結(jié)構(gòu)健康監(jiān)測(cè)中的重要環(huán)節(jié)[1-5]。

    彭鴿等[12]利用信號(hào)傳播時(shí)間(Time of fight, Tof),對(duì)健康信號(hào)與含損傷信號(hào)的差信號(hào)進(jìn)行分析,采用橢圓定位法進(jìn)行了損傷定位研究;LU Y等[13]采用時(shí)域信號(hào)包絡(luò)法,確定了時(shí)間延遲,采用橢圓定位法初步確定損傷位置,進(jìn)而采用全加的數(shù)據(jù)融合方法對(duì)多傳感路徑進(jìn)行綜合定位,并成像處理。上述方法均是基于對(duì)信號(hào)時(shí)域波形的分析,但由于Lamb波的頻散特性,即不同頻率下的Lamb波群速度不同,可能會(huì)導(dǎo)致信號(hào)時(shí)域波包的變形,造成時(shí)間延遲的計(jì)算不準(zhǔn)確,進(jìn)而導(dǎo)致定位失準(zhǔn)。為消除頻散效應(yīng)的影響,嚴(yán)宏等[14]基于小波分析及概率統(tǒng)計(jì)理論,利用健康信號(hào)與差信號(hào)的能量差作為損傷指標(biāo),提出了一種損傷存在概率成像算法?;跁r(shí)頻變換的方法對(duì)換能器接收的信號(hào)進(jìn)行處理,可有效消除Lamb波的頻散影響,進(jìn)而獲得較為準(zhǔn)確的時(shí)間延遲,最終可提高定位精度。

    本文采用Morlet小波分別對(duì)健康信號(hào)和損傷信號(hào)進(jìn)行時(shí)頻分析,通過提取激勵(lì)信號(hào)中心頻率和其附近頻率的時(shí)域信號(hào),確定多個(gè)頻率下的時(shí)間延遲,結(jié)合對(duì)應(yīng)頻率下的相速度,確定板中離散點(diǎn)到換能器對(duì)的距離之和,對(duì)多個(gè)頻率下的距離之和進(jìn)行加權(quán)處理,進(jìn)而采用橢圓定位法確定損傷位置。

    1 Lamb波傳播理論

    Lamb波是板狀結(jié)構(gòu)中產(chǎn)生的平面應(yīng)變波,且在板的上、下表面應(yīng)力為零,可劃分為對(duì)稱模態(tài)和反對(duì)稱模態(tài)。Lamb波的群速度及相速度可根據(jù)Rayleigh-Lamb波方程得到[1]:

    其中,=2為板厚。

    針對(duì)LY12CZ鋁合金材料,其彈性模量為73 MPa,泊松比為0.33,密度為2 780 kg/m3,橫波波速為3 142 m/s,縱波波速為6 237 m/s,采用二分法,通過MATLAB編程對(duì)頻散曲線進(jìn)行了數(shù)值求解,群速度和相速度頻散曲線如圖1和圖2所示。

    圖1 群速度頻散曲線

    圖2 相速度頻散曲線

    2 橢圓定位法

    式中:是損傷到兩換能器的距離之和。損傷必然發(fā)生在以兩換能器為焦點(diǎn),且到達(dá)兩焦點(diǎn)距離之和為g的橢圓上。因一對(duì)換能器只能確定一個(gè)橢圓,無法實(shí)現(xiàn)損傷定位,故采用三個(gè)以上的換能器進(jìn)行定位,多個(gè)橢圓的交點(diǎn)即為損傷位置。

    圖3 橢圓定位法示意圖

    Fig.3 The schematic of Ellipse localization method

    3 基于Morlet小波時(shí)頻分析及數(shù)據(jù)融合的橢圓定位法

    3.1 基于Morlet小波變換的時(shí)間延遲確定方法

    目前常用的時(shí)頻分析方法有短時(shí)傅里葉變換(Short-Time Fourier Transform,STFT)、Wigner-Ville分布和小波變換等。短時(shí)傅里葉變換是采用固定的窗函數(shù),窗函數(shù)一旦確定,則時(shí)頻變換的分辨率就已確定,且不能兼顧時(shí)間和頻率的分辨率,僅適合于緩變信號(hào)的分析;Wigner-Ville分布相較STFT,具有較好的分辨率,但當(dāng)信號(hào)含有多個(gè)成分時(shí),會(huì)存在交叉干擾項(xiàng),引入虛假的能量分布;小波變換則是基于STFT發(fā)展起來的,其采用了小波基函數(shù),可以通過伸縮和平移運(yùn)算對(duì)信號(hào)進(jìn)行多尺度細(xì)化,可提供一個(gè)隨頻率變化的“時(shí)間-頻率”窗口,是進(jìn)行信號(hào)時(shí)頻分析的理想手段。所以本文采用小波變換作為對(duì)換能器采集信號(hào)進(jìn)行時(shí)頻分析的工具。

    本文采用Morlet函數(shù)作為分析小波,其表達(dá)式為[15]:

    通過Morlet小波變換可得到信號(hào)幅值隨時(shí)間和頻率的變化關(guān)系。由于文中采用的頻率無法滿足因材料的非線性引起的二次諧波的條件,故可認(rèn)為Lamb波在板中傳播時(shí),僅存在激勵(lì)信號(hào)的頻率成分[16]。則考慮提取激勵(lì)信號(hào)頻率0及其附近的多組頻率所對(duì)應(yīng)的時(shí)域信號(hào),取特定頻率時(shí)域信號(hào)中第一個(gè)波峰值為信號(hào)傳播時(shí)間,通過比較健康信號(hào)及差信號(hào)的信號(hào)傳播時(shí)間,即可確定對(duì)應(yīng)頻率下的時(shí)間延遲,如圖4所示。結(jié)合前文計(jì)算得到的頻散曲線,可計(jì)算得到240 kHz頻率下的Lamb波相速度p,根據(jù)信號(hào)傳播時(shí)間及相速度可計(jì)算出損傷到兩換能器的距離之和。

    圖4 240 kHz頻率下的時(shí)域信號(hào)

    3.2 損傷位置的數(shù)據(jù)融合方法

    本文采用加權(quán)平均法作為數(shù)據(jù)融合的工具,則有[17]:

    其中,為權(quán)重值;為融合頻率的個(gè)數(shù)。

    因本文采用的激勵(lì)信號(hào)是漢寧窗調(diào)制的5個(gè)周期正弦信號(hào):

    圖5 激勵(lì)信號(hào)的時(shí)域及頻域圖

    對(duì)權(quán)值進(jìn)行歸一化處理,得到歸一化的權(quán)值:

    圖6 權(quán)重確定示意圖

    得到權(quán)值之后,即可依據(jù)公式(9),對(duì)所獲得的多頻率下的損傷到兩換能器的距離L進(jìn)行融合計(jì)算,得到可靠性較高的損傷到該組換能器的距離和,并以此法獲得其他組的定位信息,進(jìn)而采用橢圓定位法,即可確定損傷的位置。

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

    采用尺寸為700 mm×300 mm×3 mm的LY12CZ鋁合金薄板為試驗(yàn)對(duì)象。通過查材料手冊(cè)得到該平板的材料參數(shù),并通過計(jì)算得到其橫波及縱波速度,見表1。試驗(yàn)件示意圖見圖7。以板的中心位置為零點(diǎn),橫向?yàn)檩S,縱向?yàn)檩S,分別在坐標(biāo)值為(-150,80)、(150,80),(150,-80)和(-150,-80)位置處粘貼PZT壓電片,其中坐標(biāo)單位為mm,并進(jìn)行編號(hào)。在板(-50,40)位置處,設(shè)置直徑為5 mm的通孔,作為損傷狀態(tài)。分別在板健康和損傷狀態(tài)下進(jìn)行試驗(yàn)監(jiān)測(cè)。采用NI數(shù)字采集設(shè)備(Pxi-5105示波器卡、Pix-5412采集卡),HSA 4052功率放大器,試驗(yàn)設(shè)備如圖8所示。激勵(lì)信號(hào)為中心頻率250 kHz的漢寧窗調(diào)制的5周期正弦信號(hào)。信號(hào)采集頻率為1 MHz,信號(hào)采集點(diǎn)數(shù)為10 000。

    表1 材料參數(shù)及橫縱波速

    圖7 試驗(yàn)件示意圖

    圖8 試驗(yàn)設(shè)備

    為驗(yàn)證本文方法的優(yōu)越性,分別由換能器1和2激勵(lì),其他換能器接收,分別在健康和損傷狀態(tài)下獲得12組信號(hào),采用信號(hào)時(shí)域包絡(luò)法確定時(shí)間延遲。以傳播路徑1-3為例,對(duì)信號(hào)進(jìn)行濾波處理以減少噪聲影響,得到的健康信號(hào)、損傷信號(hào)和差信號(hào)及其信號(hào)包絡(luò)圖如圖9所示。

    (a) 健康信號(hào)

    (b) 損傷信號(hào)

    (c) 差信號(hào)

    下面采用Morlet小波變換和多頻率數(shù)據(jù)融合的方法對(duì)損傷進(jìn)行橢圓定位。激勵(lì)信號(hào)的中心頻率為250 kHz,-3 dB帶寬為70 kHz,取230、240、250、260和270 kHz共5組頻率作為需要提取的頻率。根據(jù)公式(12)確定了這5組頻率下?lián)p傷到換能器對(duì)的距離之和的權(quán)重,同時(shí)根據(jù)頻散曲線得到了這5組頻率下對(duì)應(yīng)的相速度,見表3。

    表2 基于信號(hào)包絡(luò)閾值法確定時(shí)間延遲

    圖10 基于時(shí)域信號(hào)包絡(luò)閾值法的損傷定位

    表3 多頻率下的相速度、權(quán)重值及時(shí)間延遲

    以換能器路徑2-3為例,首先分別對(duì)健康信號(hào)和差信號(hào)進(jìn)行Morlet小波變換,得到信號(hào)的時(shí)頻圖(見圖11),繼而分別提取上述5組頻率下的時(shí)域信號(hào)(250 kHz下的時(shí)域圖見圖12),確定了時(shí)間延遲(見表3),結(jié)合權(quán)重值,最終確定損傷到換能器2和3的距離之和。

    基于本文所提方法,對(duì)多條傳感路徑進(jìn)行分析,采用多頻率下時(shí)域信號(hào)確定時(shí)間延遲,見表4。以此并進(jìn)行加權(quán)平均處理,最終獲得多個(gè)。采用橢圓定位法進(jìn)行損傷定位,并成像,經(jīng)閾值化處理,最終確定損傷位置為(-43.8,43.2),誤差為6.98 mm,如圖13所示。

    (a) 健康信號(hào)小波時(shí)頻圖

    (b) 差信號(hào)小波時(shí)頻圖

    圖11 小波時(shí)頻圖

    Fig.11 Wavelet time-frequency distribution

    圖12 250 kHz下時(shí)域信號(hào)圖

    表4 基于小波變換和多頻率融合方法獲得的路徑數(shù)據(jù)

    圖13 基于Morlet小波變換和多頻率數(shù)據(jù)融合的損傷橢圓定位

    5 結(jié)論

    本文針對(duì)Lamb波的頻散特性對(duì)時(shí)間延遲確定的影響,提出了基于小波分析多頻率數(shù)據(jù)融合的橢圓損傷定位方法。該方法基于Morlet小波分析提取了多頻率下信號(hào)隨時(shí)間的變化信息,獲得了多頻率下的時(shí)間延遲數(shù)據(jù),根據(jù)時(shí)間延遲數(shù)據(jù),對(duì)損傷到兩換能器的距離之和進(jìn)行了數(shù)據(jù)融合,進(jìn)而基于橢圓定位法進(jìn)行了損傷定位。通過試驗(yàn)驗(yàn)證,對(duì)本文所提方法與常用的信號(hào)時(shí)域包絡(luò)法進(jìn)行了比較。結(jié)果表明,該方法可有效避免因Lamb波頻散效應(yīng)引起的時(shí)域包絡(luò)法定位不準(zhǔn)的現(xiàn)象,相較信號(hào)時(shí)域包絡(luò)法有較好的定位精度,融合了多頻率下的傳感數(shù)據(jù),提高了定位系統(tǒng)的魯棒性,可減少因系統(tǒng)或人為因素造成的誤差。

    [1] 袁慎芳. 結(jié)構(gòu)健康監(jiān)控[M]. 北京: 國防工業(yè)出版社,2007. YUAN Shenfang. Structural health monitoring[M]. Beijing: National Defence Industry Press, 2007.

    [2] 王強(qiáng), 袁慎芳, 陳小惠, 等. 主動(dòng)Lamb波合成波陣面損傷成像監(jiān)測(cè)方法[J]. 儀器儀表學(xué)報(bào), 2011, 32(11): 2468-2474. WANG Qiang, YUAN Shenfang, CHEN Xiaohui, et al. Active Lamb wave synthetic wave front damage imaging monitoring method[J]. Chinese Journal of Scientific Instrument, 2011, 32(11): 2468-2474.

    [3] 王建強(qiáng),余龍,張宇, 等. 基于最大相似性的Lamb波損傷信號(hào)分解及試驗(yàn)研究[J]. 振動(dòng)與沖擊, 2014, 33(24): 200-204. WANG Jianqiang, YU Long, ZHANG Yu, et al. Damage signal decomposition of Lamb wave and tests based on similarity[J]. Journal of Vibration and Shock, 2014, 33(24): 200-204.

    [4] Alessandro Perelli, Luca De Marchi, Alessandro Marzani, et al. Frequency warped cross-wavelet multi-resolution analysis of guided waves for impact localization[J]. Signal Processing, 2014, 96(5): 51-56.

    [5] Ben B S, Ben B A, Vikramb K A, et al. Damage identi?cation in composite materials using ultrasonic based Lamb wave method[J]. Measurement, 2013, 46(2): 904-912.

    [6] 劉增華, 徐營贊, 何存富, 等. 板狀結(jié)構(gòu)中基于Lamb波單模態(tài)的缺陷成像試驗(yàn)研究[J]. 工程力學(xué), 2014, 31(4): 232-238. LIU Zenghua, XUE Yingzan, HE Cunfu. experimental study on defect imaging based on single Lamb wave mode in plate-like structures[J]. Engineering Mechanics, 2014, 31(4): 232-238.

    [7] 徐穎娣, 袁慎芳, 彭鴿. 二維結(jié)構(gòu)損傷的主動(dòng)Lamb波定位技術(shù)研究[J]. 航空學(xué)報(bào), 2004, 25(5): 446-479. XU Yingdi, YUAN Shenfang, PENG Ge. Study on two- dimensional damage location in structure based on active Lamb wave detection technique[J]. Acta Aeronautica Et Astronautica Sinica. 2004, 25(5): 446-479.

    [8] 于勇凌, 張海燕, 馮國瑞, 等. 蘭姆波檢測(cè)信號(hào)的匹配追蹤時(shí)頻分析[J]. 聲學(xué)學(xué)報(bào), 2013, 38(5): 576-582. YU Yongling, ZHANG Haiyan, FENG Guorui, et al. Matching pursuit time-frequency analysis of Lamb wave detection signals[J]. Acta Acustica, 2013, 38(5): 576-582.

    [9] Fenza A D, Sorrentino A, Vitiello P. Application of Artificial Neural Networks and Probability Ellipse methods for damage detection using Lamb waves[J]. Composite Structures, 2015, 133: 390-403.

    [10] Pai P F, Deng H, Sundaresan M J. Time- frequency characterization of Lamb waves for material evaluation and damage inspection of plates[J]. Mechanical Systems and Signal Processing(S 0888-3270), 2015, 62-63: 183-206.

    [11] GE Luyao, WANG Xinwei, JIN Chunhua. Numerical modeling of PZT-induced Lamb wave-based crack detection in plate-like structures[J]. Wave Motion(S 0165-2125), 2014, 51(6): 867-885.

    [12] 彭鴿, 袁慎芳, 徐穎娣. 基于主動(dòng)Lamb波和小波變換的二維結(jié)構(gòu)損傷定位研究[J]. 振動(dòng)工程學(xué)報(bào), 2004, 17(4): 488-493. PENG Ge, YUAN Shenfang, XU Yingdi. Damage location of two-dimensional structure based on wavelet transform and active monitoring technology of lamb waves[J]. Journal of Vibration Engineering, 2004, 17(4): 488-493.

    [13] LU Y, YE L, WANG D, et al. Conjunctive and compromised data fusion schemes for identification of multiple notches in an aluminium plate using Lamb wave signals.[J]. IEEE Transactions on Ultrasonics Ferroelectrics & Frequency Control(S 0885-3010), 2010, 57(9): 2005-2016.

    [14] 嚴(yán)宏, 周麗. 基于損傷存在概率成像方法的復(fù)合材料結(jié)構(gòu)損傷識(shí)別[J]. 振動(dòng)與沖擊, 2012, 31(13): 76-80. YAN Hong, ZHOU Li. Damage identification of composite structures based on reconstruction algorithm for probabilistic inspection of damage[J]. Journal of Vibration and Shock, 2012, 31(13): 76-80.

    [15] Boggess Albert, Narcowich Francis J. 小波與傅里葉分析基礎(chǔ)[M]. 北京: 電子工業(yè)出版社, 2013. Boggess Albert, Narcowich Francis J. Wavelets with Fourier analysis[M]. Beijing: Publishing House of Electronics Industry, 2013.

    [16] 鄧明晰, Price D C, Scott D A. 蘭姆波非線性效應(yīng)的實(shí)驗(yàn)觀察[J]. 聲學(xué)學(xué)報(bào), 2005, 31(1): 1-7. DENG Mingxi, D.C.Price, D.A.Scott. experimental observation of nonlinear effects of Lamb waves[J]. Acta Acustica, 2005, 31(1): 1-7.

    [17] 劉同明, 夏祖勛, 解洪成. 數(shù)據(jù)融合技術(shù)及應(yīng)用[M]. 北京: 國防工業(yè)出版社, 1998: 1-12.LIU Tongming, XIA Zuxun, XIE Hongcheng. Data fusion technology and its application[M]. Beijing: National Defense Industry Press, 1998.

    Research on multi-frequency data fusion based Lamb wave damage localization

    YANG Xiao-hua1, LIU Xue-jun1, MA Guang-ting2

    (1. Department of Aviation Mechanism, Qingdao Branch, Naval Aeronautical and Astronautical University, Qingdao 266041, Shandong, China;2. Unit 91899 of PLA, Huludao 125001, Liaoning, China)

    Determination of damage location is an important part of structural health monitoring, and time delay is the key parameter to determine the damage location. A damage localization method based on Morlet wavelet transform and data fusion is put forward to reduce the influence of Lamb wave dispersion on the localization accuracy. Travel time information at several special frequencies is extracted, so the signal arrival time is determined. The distance of damage to sensors is obtained based on weighted processing. Damage localization is determined based on ellipse localization method. Simulation results show that the method can reduce the influence of dispersion effect, improve the accuracy of localization and enhance the system anti-interference ability.

    Lamb waves; Morlet wavelet; ellipse localization method; data fusion; time delay

    V214.8 TB332

    A

    1000-3630(2017)-02-0133-07

    10.16300/j.cnki.1000-3630.2017.02.007

    2016-05-20;

    2016-08-13

    總裝預(yù)研項(xiàng)目(143092015)

    楊曉華(1964-), 男, 江蘇啟東人, 博士, 教授, 博士生導(dǎo)師, 研究方向?yàn)轱w機(jī)結(jié)構(gòu)壽命評(píng)定。

    劉學(xué)君, E-mail: liuxuejun_2007@163.com

    猜你喜歡
    時(shí)頻換能器時(shí)域
    基于時(shí)域信號(hào)的三電平逆變器復(fù)合故障診斷
    一種推挽式驅(qū)動(dòng)Ⅶ型彎張換能器的探討
    IV型彎張換能器輻射阻抗的等效球模型計(jì)算
    基于極大似然準(zhǔn)則與滾動(dòng)時(shí)域估計(jì)的自適應(yīng)UKF算法
    基于時(shí)域逆濾波的寬帶脈沖聲生成技術(shù)
    超磁致伸縮復(fù)合棒換能器研究
    基于高頻發(fā)射換能器的功放設(shè)計(jì)
    基于時(shí)域波形特征的輸電線雷擊識(shí)別
    基于時(shí)頻分析的逆合成孔徑雷達(dá)成像技術(shù)
    對(duì)采樣數(shù)據(jù)序列進(jìn)行時(shí)頻分解法的改進(jìn)
    亚洲aⅴ乱码一区二区在线播放 | 欧美黑人欧美精品刺激| 精品高清国产在线一区| 免费久久久久久久精品成人欧美视频| 国产精品美女特级片免费视频播放器 | 中文字幕人妻熟女乱码| 日韩三级视频一区二区三区| 久久精品国产亚洲av香蕉五月| 黄色成人免费大全| 岛国在线观看网站| 欧美精品啪啪一区二区三区| 国产aⅴ精品一区二区三区波| 久99久视频精品免费| 久久人妻熟女aⅴ| 99国产精品免费福利视频| 国产男靠女视频免费网站| 麻豆成人av在线观看| 黄色a级毛片大全视频| 夜夜夜夜夜久久久久| 国产色视频综合| 嫩草影院精品99| 天堂动漫精品| av视频免费观看在线观看| 精品人妻1区二区| 丝袜在线中文字幕| 一本大道久久a久久精品| 9191精品国产免费久久| 精品国内亚洲2022精品成人| 久久午夜综合久久蜜桃| 亚洲成a人片在线一区二区| 人成视频在线观看免费观看| 国产精品爽爽va在线观看网站 | 亚洲一区中文字幕在线| av电影中文网址| netflix在线观看网站| 午夜久久久久精精品| 又紧又爽又黄一区二区| 久久久久精品国产欧美久久久| 最新在线观看一区二区三区| 婷婷丁香在线五月| 成人国语在线视频| 亚洲国产中文字幕在线视频| 老司机午夜福利在线观看视频| 久久久久久亚洲精品国产蜜桃av| 最新在线观看一区二区三区| 可以在线观看的亚洲视频| 嫁个100分男人电影在线观看| 满18在线观看网站| 亚洲精品久久国产高清桃花| 制服丝袜大香蕉在线| 长腿黑丝高跟| 国产午夜精品久久久久久| av视频免费观看在线观看| 日本五十路高清| 此物有八面人人有两片| 天天躁夜夜躁狠狠躁躁| 国产视频一区二区在线看| 激情在线观看视频在线高清| 久久久久国产精品人妻aⅴ院| 在线国产一区二区在线| 91大片在线观看| 国产成人一区二区三区免费视频网站| 1024视频免费在线观看| 国产精华一区二区三区| 一级作爱视频免费观看| 国产欧美日韩综合在线一区二区| a在线观看视频网站| 午夜免费激情av| 大香蕉久久成人网| www.www免费av| 一级作爱视频免费观看| 在线观看午夜福利视频| 国产精品 国内视频| 中亚洲国语对白在线视频| a级毛片在线看网站| 日日干狠狠操夜夜爽| 欧美激情高清一区二区三区| 丝袜美腿诱惑在线| 叶爱在线成人免费视频播放| 变态另类丝袜制服| 香蕉丝袜av| 99国产精品一区二区蜜桃av| 日韩欧美免费精品| 久久亚洲真实| a级毛片在线看网站| 一本久久中文字幕| 18美女黄网站色大片免费观看| 一边摸一边抽搐一进一小说| 免费看十八禁软件| 免费在线观看完整版高清| 免费少妇av软件| www.999成人在线观看| 国产精品九九99| 十八禁网站免费在线| 欧美黄色片欧美黄色片| 亚洲自拍偷在线| 悠悠久久av| 精品不卡国产一区二区三区| 每晚都被弄得嗷嗷叫到高潮| 精品电影一区二区在线| 国产午夜福利久久久久久| 看黄色毛片网站| 淫妇啪啪啪对白视频| 日韩欧美在线二视频| 久久久久久久久免费视频了| 级片在线观看| 老司机午夜十八禁免费视频| 一进一出抽搐动态| 可以在线观看的亚洲视频| 日韩大码丰满熟妇| 国产三级黄色录像| 亚洲激情在线av| 亚洲中文av在线| 国产精品日韩av在线免费观看 | 亚洲免费av在线视频| 成年人黄色毛片网站| 久久午夜亚洲精品久久| 精品无人区乱码1区二区| 老司机午夜福利在线观看视频| 制服诱惑二区| 亚洲人成电影观看| 免费女性裸体啪啪无遮挡网站| 亚洲中文av在线| 色综合站精品国产| 国产一区二区在线av高清观看| 啦啦啦 在线观看视频| 国产免费av片在线观看野外av| 宅男免费午夜| 一进一出好大好爽视频| 日本精品一区二区三区蜜桃| 丁香欧美五月| 美女免费视频网站| 人人妻人人爽人人添夜夜欢视频| 波多野结衣一区麻豆| 一本久久中文字幕| 制服诱惑二区| 久久香蕉激情| 精品国产一区二区久久| 女警被强在线播放| 久久精品国产99精品国产亚洲性色 | 又大又爽又粗| 青草久久国产| 午夜老司机福利片| 中文字幕色久视频| 九色国产91popny在线| 一级毛片高清免费大全| 91在线观看av| 亚洲精品一卡2卡三卡4卡5卡| 手机成人av网站| 搞女人的毛片| 日韩大尺度精品在线看网址 | 精品一品国产午夜福利视频| 久久精品国产清高在天天线| 国产午夜精品久久久久久| 性色av乱码一区二区三区2| 18美女黄网站色大片免费观看| 国产精华一区二区三区| 午夜老司机福利片| 性色av乱码一区二区三区2| 亚洲成av人片免费观看| 亚洲成国产人片在线观看| 欧美中文日本在线观看视频| 老司机午夜福利在线观看视频| www.精华液| 在线国产一区二区在线| 女性被躁到高潮视频| 精品福利观看| 精品国产超薄肉色丝袜足j| 在线观看免费视频日本深夜| 老熟妇仑乱视频hdxx| 欧美+亚洲+日韩+国产| 日韩成人在线观看一区二区三区| 啦啦啦免费观看视频1| 美女高潮喷水抽搐中文字幕| 免费在线观看完整版高清| 制服丝袜大香蕉在线| 视频在线观看一区二区三区| 午夜福利视频1000在线观看 | 亚洲av五月六月丁香网| av免费在线观看网站| 一进一出抽搐gif免费好疼| 亚洲国产看品久久| 国产精品,欧美在线| 黑丝袜美女国产一区| 国产伦人伦偷精品视频| 中国美女看黄片| 97超级碰碰碰精品色视频在线观看| 人妻久久中文字幕网| 久久久久久久久免费视频了| 国产成+人综合+亚洲专区| 天天躁夜夜躁狠狠躁躁| 国产成人系列免费观看| 欧美老熟妇乱子伦牲交| 国产一区二区激情短视频| 日韩高清综合在线| tocl精华| 亚洲成国产人片在线观看| 可以在线观看的亚洲视频| 丰满人妻熟妇乱又伦精品不卡| 免费一级毛片在线播放高清视频 | 国产精品二区激情视频| 中文字幕人妻丝袜一区二区| 人成视频在线观看免费观看| 国产亚洲欧美精品永久| 看免费av毛片| 成年版毛片免费区| 日本欧美视频一区| 欧美激情 高清一区二区三区| 十八禁网站免费在线| 国产成人精品无人区| 久99久视频精品免费| 人人澡人人妻人| 成人国产综合亚洲| 一区二区三区高清视频在线| 亚洲第一av免费看| 叶爱在线成人免费视频播放| 亚洲人成电影观看| 久久草成人影院| 中文字幕久久专区| 咕卡用的链子| av天堂在线播放| 亚洲片人在线观看| 成人亚洲精品一区在线观看| 日韩大码丰满熟妇| 狂野欧美激情性xxxx| 精品高清国产在线一区| 国产精品免费视频内射| 欧美av亚洲av综合av国产av| 国产精品日韩av在线免费观看 | 国产人伦9x9x在线观看| 看免费av毛片| 国产成+人综合+亚洲专区| 中文字幕精品免费在线观看视频| 精品少妇一区二区三区视频日本电影| 亚洲国产高清在线一区二区三 | 亚洲午夜理论影院| 亚洲最大成人中文| 久久久久亚洲av毛片大全| 久久精品aⅴ一区二区三区四区| 两个人免费观看高清视频| 美女高潮到喷水免费观看| 亚洲精华国产精华精| 夜夜夜夜夜久久久久| 多毛熟女@视频| 一区二区三区国产精品乱码| 黄网站色视频无遮挡免费观看| 亚洲自拍偷在线| 黄色视频不卡| 给我免费播放毛片高清在线观看| 在线观看免费视频网站a站| 男人的好看免费观看在线视频 | av电影中文网址| 欧美精品啪啪一区二区三区| 激情视频va一区二区三区| 亚洲国产精品999在线| 久久精品国产清高在天天线| 国产一区二区三区在线臀色熟女| 午夜福利在线观看吧| 亚洲精品国产一区二区精华液| 看免费av毛片| 亚洲中文字幕日韩| 欧美激情 高清一区二区三区| 搡老妇女老女人老熟妇| 热99re8久久精品国产| 欧美最黄视频在线播放免费| 国产一级毛片七仙女欲春2 | 美女免费视频网站| 长腿黑丝高跟| 日韩高清综合在线| 黄色丝袜av网址大全| 国产精品免费视频内射| 亚洲片人在线观看| 免费看美女性在线毛片视频| 国产亚洲欧美98| 一区二区三区国产精品乱码| 黄色成人免费大全| 男女做爰动态图高潮gif福利片 | 在线观看日韩欧美| 桃色一区二区三区在线观看| 看免费av毛片| 国产成人欧美| 十八禁网站免费在线| 午夜福利视频1000在线观看 | 波多野结衣高清无吗| 在线永久观看黄色视频| 亚洲七黄色美女视频| 非洲黑人性xxxx精品又粗又长| 脱女人内裤的视频| 亚洲中文字幕一区二区三区有码在线看 | 50天的宝宝边吃奶边哭怎么回事| 又紧又爽又黄一区二区| 欧美黑人精品巨大| 亚洲精品中文字幕一二三四区| 热99re8久久精品国产| netflix在线观看网站| 国产又爽黄色视频| 欧美日韩乱码在线| 99久久精品国产亚洲精品| 亚洲 欧美一区二区三区| 久久香蕉精品热| 久久精品国产亚洲av香蕉五月| 女性生殖器流出的白浆| 亚洲精华国产精华精| 色av中文字幕| 欧美乱色亚洲激情| 午夜精品国产一区二区电影| 99在线人妻在线中文字幕| 亚洲欧美一区二区三区黑人| 99久久精品国产亚洲精品| 人妻丰满熟妇av一区二区三区| 午夜两性在线视频| 国产真人三级小视频在线观看| 国产黄a三级三级三级人| 一本久久中文字幕| 欧美丝袜亚洲另类 | 欧美激情高清一区二区三区| 一级作爱视频免费观看| 久久人人爽av亚洲精品天堂| 免费无遮挡裸体视频| 天天躁夜夜躁狠狠躁躁| 女性生殖器流出的白浆| 桃红色精品国产亚洲av| 日韩一卡2卡3卡4卡2021年| www.熟女人妻精品国产| 不卡av一区二区三区| 亚洲精品一区av在线观看| 亚洲久久久国产精品| 色播亚洲综合网| 人妻久久中文字幕网| 手机成人av网站| 亚洲国产精品sss在线观看| 日韩视频一区二区在线观看| 可以在线观看的亚洲视频| 一级毛片女人18水好多| 国产欧美日韩一区二区三区在线| 在线观看免费视频网站a站| 亚洲国产欧美一区二区综合| 国产亚洲欧美精品永久| 妹子高潮喷水视频| 国产精品久久久av美女十八| 成在线人永久免费视频| 成人三级做爰电影| 不卡av一区二区三区| 男女做爰动态图高潮gif福利片 | 欧美日韩一级在线毛片| 9色porny在线观看| 久久久久久免费高清国产稀缺| 亚洲av成人av| 嫁个100分男人电影在线观看| 人人澡人人妻人| 麻豆久久精品国产亚洲av| 久久久水蜜桃国产精品网| 国产精品秋霞免费鲁丝片| 一区二区三区国产精品乱码| 波多野结衣高清无吗| 极品教师在线免费播放| 人人妻人人澡人人看| 亚洲中文日韩欧美视频| 国产男靠女视频免费网站| 国产黄a三级三级三级人| 久久久精品国产亚洲av高清涩受| 啦啦啦观看免费观看视频高清 | 久久人妻熟女aⅴ| 色综合欧美亚洲国产小说| 老熟妇乱子伦视频在线观看| 中文字幕人成人乱码亚洲影| www.www免费av| 免费观看精品视频网站| 乱人伦中国视频| 欧美在线一区亚洲| 一卡2卡三卡四卡精品乱码亚洲| 777久久人妻少妇嫩草av网站| 1024香蕉在线观看| 桃红色精品国产亚洲av| 制服人妻中文乱码| 欧美日韩黄片免| 亚洲无线在线观看| 桃色一区二区三区在线观看| 亚洲一区二区三区不卡视频| 男女下面进入的视频免费午夜 | 久久亚洲精品不卡| 电影成人av| 国产精品一区二区在线不卡| 午夜精品国产一区二区电影| 成人亚洲精品一区在线观看| 亚洲免费av在线视频| 成人国产综合亚洲| 日韩精品免费视频一区二区三区| 久久精品国产亚洲av高清一级| 国产人伦9x9x在线观看| 99国产精品一区二区三区| 国产一区在线观看成人免费| 欧美成人一区二区免费高清观看 | 亚洲 欧美一区二区三区| 婷婷精品国产亚洲av在线| 日韩av在线大香蕉| 精品一区二区三区av网在线观看| 91字幕亚洲| 亚洲精华国产精华精| 国产99久久九九免费精品| 婷婷精品国产亚洲av在线| 69av精品久久久久久| 高潮久久久久久久久久久不卡| 久久久国产欧美日韩av| 亚洲成av人片免费观看| 欧美日本视频| 少妇被粗大的猛进出69影院| 亚洲一卡2卡3卡4卡5卡精品中文| 久久精品国产综合久久久| 黄色女人牲交| 18禁黄网站禁片午夜丰满| 99国产精品一区二区三区| 婷婷六月久久综合丁香| 很黄的视频免费| 欧美日韩精品网址| 成在线人永久免费视频| 久久久久国产一级毛片高清牌| 老司机深夜福利视频在线观看| 男女床上黄色一级片免费看| 两人在一起打扑克的视频| 国产精品亚洲一级av第二区| 亚洲七黄色美女视频| 国内精品久久久久精免费| 久久狼人影院| 欧美人与性动交α欧美精品济南到| 精品国产乱子伦一区二区三区| av天堂在线播放| 在线观看午夜福利视频| 男男h啪啪无遮挡| 一级a爱片免费观看的视频| 久久久久国内视频| 国产成人免费无遮挡视频| 啦啦啦免费观看视频1| 18美女黄网站色大片免费观看| 波多野结衣巨乳人妻| 国产精品久久久人人做人人爽| bbb黄色大片| 侵犯人妻中文字幕一二三四区| 黄色视频,在线免费观看| 中文字幕最新亚洲高清| 此物有八面人人有两片| 国产成人影院久久av| 黑人操中国人逼视频| 免费av毛片视频| 国产区一区二久久| 亚洲无线在线观看| 最近最新免费中文字幕在线| 一级毛片高清免费大全| 免费av毛片视频| 97超级碰碰碰精品色视频在线观看| 欧美成人午夜精品| 黄频高清免费视频| 91精品国产国语对白视频| 男女床上黄色一级片免费看| 色尼玛亚洲综合影院| 美女午夜性视频免费| a在线观看视频网站| 热re99久久国产66热| 两个人免费观看高清视频| 在线观看免费午夜福利视频| 十八禁网站免费在线| 日本vs欧美在线观看视频| 麻豆成人av在线观看| 极品教师在线免费播放| 久久人人爽av亚洲精品天堂| 国产又色又爽无遮挡免费看| 久久九九热精品免费| 国产精品免费视频内射| 久久中文字幕人妻熟女| 看免费av毛片| 在线观看一区二区三区| 久久久精品国产亚洲av高清涩受| 亚洲 欧美 日韩 在线 免费| 69av精品久久久久久| 国产精品野战在线观看| 国产av一区二区精品久久| 99久久综合精品五月天人人| 国产高清视频在线播放一区| 18禁裸乳无遮挡免费网站照片 | 在线观看免费午夜福利视频| 国产精品免费视频内射| videosex国产| 制服人妻中文乱码| 欧美日本亚洲视频在线播放| 午夜精品久久久久久毛片777| 一级a爱视频在线免费观看| 国内毛片毛片毛片毛片毛片| 中文字幕久久专区| 日本 av在线| 97超级碰碰碰精品色视频在线观看| 国产一级毛片七仙女欲春2 | 免费高清在线观看日韩| 亚洲欧洲精品一区二区精品久久久| 亚洲少妇的诱惑av| 国产欧美日韩一区二区精品| 激情在线观看视频在线高清| 90打野战视频偷拍视频| 怎么达到女性高潮| 国产免费av片在线观看野外av| 首页视频小说图片口味搜索| ponron亚洲| 香蕉国产在线看| 在线十欧美十亚洲十日本专区| 视频区欧美日本亚洲| 九色国产91popny在线| 午夜影院日韩av| 男人操女人黄网站| 久久国产精品影院| 免费观看人在逋| 亚洲午夜理论影院| 美女免费视频网站| 国产麻豆成人av免费视频| 麻豆久久精品国产亚洲av| av网站免费在线观看视频| 天天一区二区日本电影三级 | 91老司机精品| 他把我摸到了高潮在线观看| 亚洲一区高清亚洲精品| 国产精品影院久久| 久久久久国内视频| 亚洲 欧美一区二区三区| 在线免费观看的www视频| 亚洲男人的天堂狠狠| 无遮挡黄片免费观看| 国产精品美女特级片免费视频播放器 | 午夜免费鲁丝| 精品人妻在线不人妻| 女人精品久久久久毛片| 国产欧美日韩一区二区三区在线| 国产色视频综合| 精品不卡国产一区二区三区| 国产精品久久久久久人妻精品电影| 少妇熟女aⅴ在线视频| 国产av精品麻豆| 国产亚洲精品综合一区在线观看 | 久久热在线av| 丰满的人妻完整版| 欧洲精品卡2卡3卡4卡5卡区| 在线视频色国产色| 女人精品久久久久毛片| 亚洲专区中文字幕在线| 午夜福利在线观看吧| 国产精品免费视频内射| 国产精品久久视频播放| 亚洲自偷自拍图片 自拍| 90打野战视频偷拍视频| 窝窝影院91人妻| 啦啦啦 在线观看视频| 欧美久久黑人一区二区| 久久中文看片网| 日韩欧美国产在线观看| 国产蜜桃级精品一区二区三区| 日本三级黄在线观看| 国产欧美日韩一区二区三| 法律面前人人平等表现在哪些方面| 美女午夜性视频免费| 老司机午夜十八禁免费视频| 12—13女人毛片做爰片一| 久99久视频精品免费| 波多野结衣高清无吗| av片东京热男人的天堂| 久久人人精品亚洲av| 亚洲免费av在线视频| 波多野结衣av一区二区av| 精品国产一区二区久久| 999久久久精品免费观看国产| 99riav亚洲国产免费| 不卡一级毛片| 给我免费播放毛片高清在线观看| 黄色毛片三级朝国网站| 久久精品影院6| 啦啦啦观看免费观看视频高清 | 法律面前人人平等表现在哪些方面| 中文亚洲av片在线观看爽| 电影成人av| 日韩精品中文字幕看吧| 69av精品久久久久久| 日韩精品中文字幕看吧| 最近最新免费中文字幕在线| 国产乱人伦免费视频| 国产精品自产拍在线观看55亚洲| 国产精品亚洲美女久久久| АⅤ资源中文在线天堂| 成人精品一区二区免费| 久久中文字幕人妻熟女| 免费在线观看视频国产中文字幕亚洲| 人成视频在线观看免费观看| 国产精品乱码一区二三区的特点 | 国产伦一二天堂av在线观看| 中文字幕另类日韩欧美亚洲嫩草| 天天一区二区日本电影三级 | 香蕉久久夜色| 日本a在线网址| 亚洲精品一卡2卡三卡4卡5卡| 又大又爽又粗| 每晚都被弄得嗷嗷叫到高潮| 精品无人区乱码1区二区| 亚洲欧美日韩无卡精品| 美女扒开内裤让男人捅视频| 国产蜜桃级精品一区二区三区| 中文字幕久久专区| 午夜福利免费观看在线| 国产精品亚洲美女久久久| 久久久久久久久久久久大奶| 亚洲熟女毛片儿| 久久香蕉激情| 亚洲精品一区av在线观看| 亚洲第一欧美日韩一区二区三区| 色综合亚洲欧美另类图片| 久久人妻av系列| 亚洲av日韩精品久久久久久密| 久久精品国产综合久久久| 51午夜福利影视在线观看| 亚洲专区国产一区二区| 麻豆国产av国片精品|