• <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)
    日韩精品免费视频一区二区三区 | 日韩制服丝袜自拍偷拍| 精品少妇久久久久久888优播| av免费观看日本| 十八禁网站网址无遮挡| 中文字幕制服av| 国产1区2区3区精品| 日日撸夜夜添| 交换朋友夫妻互换小说| 久久亚洲国产成人精品v| 香蕉精品网在线| 美女福利国产在线| 国产精品蜜桃在线观看| 国产1区2区3区精品| 精品久久蜜臀av无| 中文天堂在线官网| 免费看不卡的av| 婷婷色av中文字幕| 51国产日韩欧美| 久久97久久精品| 久久久久精品久久久久真实原创| 一级毛片黄色毛片免费观看视频| 丝袜喷水一区| 97精品久久久久久久久久精品| 在线观看美女被高潮喷水网站| 日韩精品有码人妻一区| 精品国产露脸久久av麻豆| 成人午夜精彩视频在线观看| 天堂俺去俺来也www色官网| 亚洲欧美日韩卡通动漫| 久久久久视频综合| 交换朋友夫妻互换小说| 国产一区二区三区综合在线观看 | 欧美xxⅹ黑人| 秋霞伦理黄片| 久久99热6这里只有精品| 久久这里只有精品19| 色网站视频免费| 91成人精品电影| www日本在线高清视频| 97在线视频观看| 七月丁香在线播放| 寂寞人妻少妇视频99o| 国产精品国产av在线观看| 蜜桃国产av成人99| 9191精品国产免费久久| 啦啦啦在线观看免费高清www| 亚洲精品,欧美精品| 国产精品.久久久| 日韩精品免费视频一区二区三区 | 精品一区二区三卡| 久久久久久久大尺度免费视频| 国产老妇伦熟女老妇高清| 视频中文字幕在线观看| 成人国语在线视频| 国产极品天堂在线| 久久久国产精品麻豆| 国产又爽黄色视频| 欧美丝袜亚洲另类| 精品国产国语对白av| 国产成人精品婷婷| 日本与韩国留学比较| 久久久久久久久久久免费av| tube8黄色片| 久久久久国产精品人妻一区二区| 天堂8中文在线网| 亚洲av成人精品一二三区| 久久人妻熟女aⅴ| 国产xxxxx性猛交| 99热全是精品| 日本爱情动作片www.在线观看| 日韩中字成人| 熟女电影av网| 久久精品aⅴ一区二区三区四区 | 亚洲成国产人片在线观看| 日韩人妻精品一区2区三区| 在线精品无人区一区二区三| av黄色大香蕉| 日韩 亚洲 欧美在线| 看免费av毛片| 少妇的丰满在线观看| 国产精品国产三级专区第一集| 亚洲国产日韩一区二区| 乱码一卡2卡4卡精品| 在线免费观看不下载黄p国产| 卡戴珊不雅视频在线播放| 少妇的逼好多水| 69精品国产乱码久久久| 制服丝袜香蕉在线| 国产在线视频一区二区| 精品视频人人做人人爽| 国产一级毛片在线| 亚洲国产色片| 精品少妇内射三级| 搡老乐熟女国产| 高清欧美精品videossex| 视频中文字幕在线观看| 国产国拍精品亚洲av在线观看| 菩萨蛮人人尽说江南好唐韦庄| 91aial.com中文字幕在线观看| 香蕉丝袜av| 久久av网站| 亚洲精品一二三| 免费少妇av软件| 日韩制服丝袜自拍偷拍| 99热6这里只有精品| 日韩免费高清中文字幕av| 美女xxoo啪啪120秒动态图| 日日撸夜夜添| av播播在线观看一区| 免费女性裸体啪啪无遮挡网站| 深夜精品福利| 日韩电影二区| 女的被弄到高潮叫床怎么办| 在线免费观看不下载黄p国产| 9191精品国产免费久久| 国产精品一区二区在线不卡| 亚洲精品乱久久久久久| 69精品国产乱码久久久| 久久免费观看电影| 国产福利在线免费观看视频| 99精国产麻豆久久婷婷| 国产欧美亚洲国产| 午夜福利,免费看| 国产精品人妻久久久久久| 在线观看免费视频网站a站| 一级爰片在线观看| 老熟女久久久| 三上悠亚av全集在线观看| a级毛色黄片| 国产精品久久久久久精品电影小说| 91精品国产国语对白视频| 美女xxoo啪啪120秒动态图| 国产国拍精品亚洲av在线观看| 国产麻豆69| 日本av免费视频播放| 母亲3免费完整高清在线观看 | 亚洲精品中文字幕在线视频| 一级黄片播放器| 国产 精品1| 美女国产高潮福利片在线看| 久久国内精品自在自线图片| 王馨瑶露胸无遮挡在线观看| 最近中文字幕高清免费大全6| 久久影院123| 伦理电影免费视频| 亚洲欧洲日产国产| 男女国产视频网站| 一级a做视频免费观看| 国产一区亚洲一区在线观看| 菩萨蛮人人尽说江南好唐韦庄| 制服诱惑二区| 亚洲激情五月婷婷啪啪| 亚洲国产精品一区二区三区在线| 波野结衣二区三区在线| 亚洲精品美女久久久久99蜜臀 | 狠狠精品人妻久久久久久综合| 韩国av在线不卡| 欧美日韩av久久| 亚洲四区av| 女性被躁到高潮视频| 又黄又粗又硬又大视频| 丝袜美足系列| 日日爽夜夜爽网站| 久久久久人妻精品一区果冻| 色哟哟·www| 18禁裸乳无遮挡动漫免费视频| 草草在线视频免费看| av免费观看日本| 一区二区三区精品91| 国产永久视频网站| 男女免费视频国产| 你懂的网址亚洲精品在线观看| 日本-黄色视频高清免费观看| 熟女电影av网| 男女下面插进去视频免费观看 | 国产日韩一区二区三区精品不卡| 国产一区亚洲一区在线观看| 精品久久久久久电影网| 国产精品秋霞免费鲁丝片| 桃花免费在线播放| 免费在线观看黄色视频的| 中文欧美无线码| 超色免费av| 2022亚洲国产成人精品| 91午夜精品亚洲一区二区三区| 伊人亚洲综合成人网| 国产极品天堂在线| 天天操日日干夜夜撸| 在线看a的网站| 久久99一区二区三区| av网站免费在线观看视频| av网站免费在线观看视频| 亚洲av中文av极速乱| 这个男人来自地球电影免费观看 | 日本爱情动作片www.在线观看| 一二三四中文在线观看免费高清| 热re99久久国产66热| 国产精品久久久久成人av| 亚洲国产精品成人久久小说| 国产黄色视频一区二区在线观看| 久久99热6这里只有精品| 久热久热在线精品观看| 久久久久久人妻| 97在线人人人人妻| 国产高清三级在线| 巨乳人妻的诱惑在线观看| 精品国产一区二区三区四区第35| 久久国产精品男人的天堂亚洲 | 日韩制服骚丝袜av| 极品人妻少妇av视频| 各种免费的搞黄视频| 夫妻午夜视频| 亚洲情色 制服丝袜| 久久 成人 亚洲| 国产在视频线精品| 日韩 亚洲 欧美在线| 超碰97精品在线观看| 老熟女久久久| 久久精品熟女亚洲av麻豆精品| 亚洲,欧美精品.| 国产成人精品在线电影| 国产福利在线免费观看视频| 亚洲av国产av综合av卡| 五月伊人婷婷丁香| 亚洲av.av天堂| 久久久国产一区二区| 亚洲欧洲精品一区二区精品久久久 | 看十八女毛片水多多多| 亚洲成人av在线免费| 国产免费一级a男人的天堂| 免费高清在线观看视频在线观看| 日韩中字成人| 18禁裸乳无遮挡动漫免费视频| 夫妻午夜视频| 久久久久久久精品精品| 咕卡用的链子| 最近的中文字幕免费完整| 中文字幕免费在线视频6| 国产日韩欧美在线精品| 中国美白少妇内射xxxbb| 免费高清在线观看日韩| 国产精品 国内视频| 国产极品粉嫩免费观看在线| 日韩中文字幕视频在线看片| 男女高潮啪啪啪动态图| 91在线精品国自产拍蜜月| 亚洲欧美清纯卡通| 人妻系列 视频| 伊人久久国产一区二区| 秋霞伦理黄片| 成人亚洲精品一区在线观看| av网站免费在线观看视频| 99香蕉大伊视频| 欧美精品国产亚洲| 综合色丁香网| 国语对白做爰xxxⅹ性视频网站| a 毛片基地| 内地一区二区视频在线| 久热这里只有精品99| 伦理电影大哥的女人| 精品福利永久在线观看| 哪个播放器可以免费观看大片| 七月丁香在线播放| 伦理电影免费视频| 国产精品久久久久久av不卡| 美女国产高潮福利片在线看| 国产精品久久久久久精品古装| 丝袜喷水一区| 亚洲国产精品999| 99国产综合亚洲精品| 亚洲精品视频女| 亚洲伊人久久精品综合| 国产亚洲午夜精品一区二区久久| 最近中文字幕高清免费大全6| 我要看黄色一级片免费的| 在线免费观看不下载黄p国产| 亚洲精品国产av蜜桃| 国产xxxxx性猛交| 99re6热这里在线精品视频| 国产av精品麻豆| 午夜日本视频在线| 欧美人与性动交α欧美精品济南到 | 少妇的逼好多水| 91在线精品国自产拍蜜月| 女性被躁到高潮视频| 中文字幕制服av| 亚洲精品一二三| 国产激情久久老熟女| 久久 成人 亚洲| 在线免费观看不下载黄p国产| 亚洲伊人色综图| 国产毛片在线视频| 免费观看av网站的网址| 日本色播在线视频| 一边亲一边摸免费视频| 亚洲,一卡二卡三卡| 国产片特级美女逼逼视频| 黄色 视频免费看| 欧美97在线视频| av有码第一页| 国产永久视频网站| 最新中文字幕久久久久| 男女国产视频网站| 咕卡用的链子| 国产69精品久久久久777片| 久久久久久久亚洲中文字幕| 久久精品国产亚洲av涩爱| 国产福利在线免费观看视频| 国产麻豆69| 国产成人av激情在线播放| 国产 精品1| 一个人免费看片子| 黄色一级大片看看| 免费高清在线观看视频在线观看| xxx大片免费视频| 啦啦啦视频在线资源免费观看| 免费黄色在线免费观看| 在线观看人妻少妇| 国产成人精品福利久久| 免费高清在线观看视频在线观看| 久久久国产一区二区| 一本久久精品| 18禁观看日本| 欧美+日韩+精品| 2022亚洲国产成人精品| 色哟哟·www| 亚洲四区av| 十八禁高潮呻吟视频| 国产在线视频一区二区| 一级,二级,三级黄色视频| 爱豆传媒免费全集在线观看| 久久国产亚洲av麻豆专区| 麻豆乱淫一区二区| 巨乳人妻的诱惑在线观看| 亚洲中文av在线| 男女边摸边吃奶| 一本—道久久a久久精品蜜桃钙片| 日本av手机在线免费观看| 欧美精品亚洲一区二区| 乱人伦中国视频| 国产在线视频一区二区| a级毛片在线看网站| 日日撸夜夜添| 久久青草综合色| 成人二区视频| 欧美变态另类bdsm刘玥| 人人妻人人澡人人看| 男女高潮啪啪啪动态图| 成人二区视频| 国产免费一级a男人的天堂| 久久午夜福利片| 欧美老熟妇乱子伦牲交| 精品福利永久在线观看| 色网站视频免费| 一区二区av电影网| 久久国内精品自在自线图片| 久久久久久人妻| 精品国产露脸久久av麻豆| 国产精品三级大全| 一级黄片播放器| 熟女av电影| 制服诱惑二区| 最近手机中文字幕大全| av线在线观看网站| 日韩一本色道免费dvd| 免费观看无遮挡的男女| 成人综合一区亚洲| 久久久久精品久久久久真实原创| 99视频精品全部免费 在线| 亚洲婷婷狠狠爱综合网| 国产色婷婷99| 久久久久精品久久久久真实原创| 日韩 亚洲 欧美在线| 精品福利永久在线观看| 精品少妇内射三级| 另类亚洲欧美激情| 成人国产av品久久久| 久久午夜综合久久蜜桃| av女优亚洲男人天堂| 久久亚洲国产成人精品v| 免费av中文字幕在线| 丁香六月天网| 久久久久国产网址| 免费观看性生交大片5| 涩涩av久久男人的天堂| 国产黄色免费在线视频| 在线 av 中文字幕| 亚洲av日韩在线播放| av在线播放精品| 两个人免费观看高清视频| 免费人妻精品一区二区三区视频| 99热网站在线观看| 又黄又爽又刺激的免费视频.| 国产精品无大码| 热99国产精品久久久久久7| a级片在线免费高清观看视频| 日本猛色少妇xxxxx猛交久久| 国产在线一区二区三区精| 亚洲av国产av综合av卡| 咕卡用的链子| 一个人免费看片子| 亚洲精品色激情综合| 在线亚洲精品国产二区图片欧美| 国产永久视频网站| videossex国产| 久久毛片免费看一区二区三区| 精品国产一区二区久久| 性色av一级| 在线天堂最新版资源| 精品第一国产精品| 日韩欧美一区视频在线观看| 日韩伦理黄色片| 精品一区二区三区视频在线| 免费看不卡的av| 国产黄频视频在线观看| 99久久精品国产国产毛片| 久久久精品免费免费高清| 这个男人来自地球电影免费观看 | 大话2 男鬼变身卡| 久久久久久久精品精品| 久久婷婷青草| 亚洲国产av新网站| 国产日韩欧美视频二区| 中文字幕另类日韩欧美亚洲嫩草| 久久99热这里只频精品6学生| 欧美成人精品欧美一级黄| 免费av不卡在线播放| 如日韩欧美国产精品一区二区三区| 久久精品久久久久久噜噜老黄| 夫妻性生交免费视频一级片| 又大又黄又爽视频免费| av免费观看日本| 伦理电影免费视频| 美女国产高潮福利片在线看| 国产又爽黄色视频| 精品一区二区三区四区五区乱码 | 一级a做视频免费观看| 人妻系列 视频| 熟女电影av网| 一区二区三区四区激情视频| 久久久久久人妻| 午夜视频国产福利| 在线看a的网站| 精品一区二区三区四区五区乱码 | 免费高清在线观看日韩| 成人黄色视频免费在线看| 欧美性感艳星| 欧美3d第一页| av在线app专区| 久久久久久久久久人人人人人人| 99精国产麻豆久久婷婷| 26uuu在线亚洲综合色| 亚洲精品久久午夜乱码| 欧美丝袜亚洲另类| 岛国毛片在线播放| av播播在线观看一区| 久久人人爽av亚洲精品天堂| 中文字幕亚洲精品专区| 超色免费av| 2021少妇久久久久久久久久久| 黄色配什么色好看| 精品第一国产精品| 另类亚洲欧美激情| 国产一区二区三区av在线| 亚洲精品日本国产第一区| 亚洲精品久久成人aⅴ小说| 草草在线视频免费看| 性色av一级| 日韩av不卡免费在线播放| 亚洲精品国产av成人精品| 久久久久精品人妻al黑| 少妇的逼好多水| 国产日韩欧美亚洲二区| av在线播放精品| 两个人免费观看高清视频| 国产欧美另类精品又又久久亚洲欧美| 欧美 日韩 精品 国产| 22中文网久久字幕| 美女xxoo啪啪120秒动态图| 国产淫语在线视频| 菩萨蛮人人尽说江南好唐韦庄| 久久人人爽av亚洲精品天堂| 在线观看www视频免费| 满18在线观看网站| 大片免费播放器 马上看| 久久久久久久久久成人| 精品一区在线观看国产| 国产成人精品一,二区| 大香蕉久久成人网| 美女国产高潮福利片在线看| 狠狠婷婷综合久久久久久88av| 国产精品久久久久久av不卡| 亚洲,欧美精品.| 欧美精品av麻豆av| 美女中出高潮动态图| www.色视频.com| 中文字幕亚洲精品专区| 久久99精品国语久久久| 啦啦啦视频在线资源免费观看| 欧美激情 高清一区二区三区| 久久久久视频综合| 亚洲精品成人av观看孕妇| 亚洲一区二区三区欧美精品| 国产精品久久久久成人av| 成人亚洲精品一区在线观看| 观看美女的网站| 考比视频在线观看| 中文欧美无线码| 黄色视频在线播放观看不卡| 午夜福利在线观看免费完整高清在| 国产激情久久老熟女| 在线观看人妻少妇| 日本wwww免费看| 久久久久久人人人人人| 91国产中文字幕| 欧美精品亚洲一区二区| 国产白丝娇喘喷水9色精品| 精品少妇久久久久久888优播| 国产精品不卡视频一区二区| 男女高潮啪啪啪动态图| 青春草视频在线免费观看| 日本vs欧美在线观看视频| 女的被弄到高潮叫床怎么办| 亚洲精品日韩在线中文字幕| av女优亚洲男人天堂| 久久久久久久久久人人人人人人| 国产免费一区二区三区四区乱码| 久久女婷五月综合色啪小说| 高清不卡的av网站| 亚洲精品美女久久av网站| 久久综合国产亚洲精品| 国产午夜精品一二区理论片| 2018国产大陆天天弄谢| av在线老鸭窝| 午夜激情久久久久久久| 国产免费又黄又爽又色| 各种免费的搞黄视频| 精品国产国语对白av| 欧美激情国产日韩精品一区| 国产亚洲欧美精品永久| 久久国产精品男人的天堂亚洲 | 欧美 日韩 精品 国产| 91aial.com中文字幕在线观看| 日韩一区二区三区影片| 亚洲激情五月婷婷啪啪| a级毛片黄视频| 男女免费视频国产| 春色校园在线视频观看| 大陆偷拍与自拍| 国产精品三级大全| 欧美变态另类bdsm刘玥| 国产精品欧美亚洲77777| 日韩精品免费视频一区二区三区 | 七月丁香在线播放| 丝袜喷水一区| 一本色道久久久久久精品综合| 亚洲国产精品国产精品| 少妇人妻久久综合中文| 高清在线视频一区二区三区| 国产 一区精品| 婷婷成人精品国产| 婷婷色综合大香蕉| 国产精品国产av在线观看| 男女免费视频国产| 91精品伊人久久大香线蕉| 亚洲欧美成人精品一区二区| 少妇的丰满在线观看| 亚洲精品国产av蜜桃| 亚洲熟女精品中文字幕| 天堂俺去俺来也www色官网| 亚洲激情五月婷婷啪啪| www日本在线高清视频| 精品熟女少妇av免费看| 国产在视频线精品| 九九在线视频观看精品| 99九九在线精品视频| 90打野战视频偷拍视频| 国产综合精华液| 欧美+日韩+精品| 黄片播放在线免费| 欧美+日韩+精品| 国产综合精华液| 欧美激情国产日韩精品一区| 精品第一国产精品| 蜜桃在线观看..| 欧美少妇被猛烈插入视频| 日韩一区二区视频免费看| 国产精品免费大片| 亚洲精品久久成人aⅴ小说| 2018国产大陆天天弄谢| 国产日韩欧美在线精品| 国产精品人妻久久久影院| 在线观看一区二区三区激情| 各种免费的搞黄视频| 亚洲第一av免费看| 大片电影免费在线观看免费| 国产精品久久久av美女十八| 大码成人一级视频| 男女无遮挡免费网站观看| 大码成人一级视频| 精品视频人人做人人爽| 久久人人爽人人爽人人片va| 美女中出高潮动态图| 色哟哟·www| 久久久欧美国产精品| 国产日韩一区二区三区精品不卡| 亚洲,欧美精品.| 黄色怎么调成土黄色| www.av在线官网国产| 国产精品久久久久久av不卡| 夫妻性生交免费视频一级片| 美国免费a级毛片| 一级毛片黄色毛片免费观看视频| 国产亚洲精品久久久com|