• <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)
    一级毛片久久久久久久久女| 午夜久久久久精精品| 国产又色又爽无遮挡免| 男女啪啪激烈高潮av片| 高清av免费在线| 好男人视频免费观看在线| 啦啦啦啦在线视频资源| 亚洲av中文字字幕乱码综合| 国产免费福利视频在线观看| 久久久久久久久大av| 国产av不卡久久| 精品午夜福利在线看| 色噜噜av男人的天堂激情| 波多野结衣高清无吗| 亚洲电影在线观看av| 大香蕉久久网| 国产精品无大码| 综合色av麻豆| 99在线视频只有这里精品首页| 乱码一卡2卡4卡精品| 亚洲经典国产精华液单| 国产精品一及| 99久国产av精品国产电影| av黄色大香蕉| 一级黄色大片毛片| 国产精品,欧美在线| 91在线精品国自产拍蜜月| 青春草国产在线视频| 国产精品野战在线观看| 男人舔奶头视频| 99在线人妻在线中文字幕| 男插女下体视频免费在线播放| 国产久久久一区二区三区| 日韩,欧美,国产一区二区三区 | 热99re8久久精品国产| av在线老鸭窝| 午夜福利网站1000一区二区三区| 成年女人永久免费观看视频| 99热6这里只有精品| 99热这里只有是精品50| 中文天堂在线官网| 狂野欧美激情性xxxx在线观看| 久99久视频精品免费| 男女那种视频在线观看| 国产精品综合久久久久久久免费| 国产高清国产精品国产三级 | 一边摸一边抽搐一进一小说| 国产成人福利小说| 国产成人精品婷婷| 草草在线视频免费看| 欧美高清性xxxxhd video| 久久6这里有精品| 校园人妻丝袜中文字幕| 一个人观看的视频www高清免费观看| 乱系列少妇在线播放| 国产久久久一区二区三区| 一级毛片aaaaaa免费看小| 看免费成人av毛片| 免费观看在线日韩| 黄色配什么色好看| 91狼人影院| 久久亚洲国产成人精品v| 岛国毛片在线播放| 国产女主播在线喷水免费视频网站 | 校园人妻丝袜中文字幕| 极品教师在线视频| 欧美日韩综合久久久久久| 久久这里只有精品中国| 亚洲最大成人手机在线| 中文字幕av在线有码专区| 日韩一本色道免费dvd| 精品久久国产蜜桃| 韩国高清视频一区二区三区| 中国美白少妇内射xxxbb| 免费看a级黄色片| 国产精品精品国产色婷婷| 国产精品国产三级国产专区5o | 欧美成人一区二区免费高清观看| 尤物成人国产欧美一区二区三区| 国产精品一区二区性色av| 国产爱豆传媒在线观看| 久久精品影院6| 国产精品久久久久久精品电影小说 | 日韩精品有码人妻一区| 日本一二三区视频观看| 国产在线男女| 黄片无遮挡物在线观看| 国产成人a∨麻豆精品| 级片在线观看| 99视频精品全部免费 在线| 99久久中文字幕三级久久日本| 插阴视频在线观看视频| 少妇的逼水好多| 神马国产精品三级电影在线观看| 日本av手机在线免费观看| 久久草成人影院| 九草在线视频观看| 久久久久久久久大av| 日本一二三区视频观看| 精品久久久久久久人妻蜜臀av| 欧美又色又爽又黄视频| 一夜夜www| 亚洲av不卡在线观看| 欧美激情久久久久久爽电影| 国产精品av视频在线免费观看| 欧美成人一区二区免费高清观看| av在线老鸭窝| 日韩精品青青久久久久久| 亚洲av福利一区| 黄片wwwwww| 真实男女啪啪啪动态图| 99久国产av精品| 直男gayav资源| 免费看日本二区| 精品久久久噜噜| 国产高清三级在线| 又黄又爽又刺激的免费视频.| 日本黄色片子视频| av黄色大香蕉| 熟女电影av网| 国产精品爽爽va在线观看网站| 久久久久性生活片| 亚洲av二区三区四区| 高清日韩中文字幕在线| 最近2019中文字幕mv第一页| 天天一区二区日本电影三级| 日日撸夜夜添| 少妇人妻精品综合一区二区| 伦理电影大哥的女人| www.色视频.com| 久久国内精品自在自线图片| 免费av观看视频| 国产一区二区三区av在线| 男女视频在线观看网站免费| 国产日韩欧美在线精品| 男人舔奶头视频| av在线播放精品| 精品人妻一区二区三区麻豆| 国产免费福利视频在线观看| 国产黄片视频在线免费观看| 尤物成人国产欧美一区二区三区| 亚洲精品久久久久久婷婷小说 | 91精品国产九色| 国产高清不卡午夜福利| 嫩草影院精品99| 99久国产av精品| 国产av码专区亚洲av| 久久久久久久久大av| 国产真实乱freesex| 国产伦在线观看视频一区| 少妇丰满av| 亚洲国产最新在线播放| 日日摸夜夜添夜夜添av毛片| 久久久久久伊人网av| 久久久色成人| 国内精品一区二区在线观看| 国产av不卡久久| 国产69精品久久久久777片| 国产精品国产高清国产av| 久久久久久久久大av| 能在线免费观看的黄片| 伦理电影大哥的女人| 国产精品久久久久久精品电影| 毛片一级片免费看久久久久| 男人的好看免费观看在线视频| 晚上一个人看的免费电影| 女人被狂操c到高潮| 久久人妻av系列| 国产成人一区二区在线| 黄色配什么色好看| 久久精品人妻少妇| 国产亚洲午夜精品一区二区久久 | 在现免费观看毛片| 国产成人一区二区在线| 在线观看66精品国产| av女优亚洲男人天堂| av女优亚洲男人天堂| 国产视频内射| 又爽又黄无遮挡网站| 嘟嘟电影网在线观看| 午夜精品一区二区三区免费看| 小说图片视频综合网站| 97人妻精品一区二区三区麻豆| 国产亚洲最大av| 久久久国产成人免费| 免费观看精品视频网站| 神马国产精品三级电影在线观看| 最近手机中文字幕大全| 欧美人与善性xxx| 七月丁香在线播放| 国产精品蜜桃在线观看| 午夜福利高清视频| 午夜免费激情av| 日韩欧美在线乱码| 一本—道久久a久久精品蜜桃钙片 精品乱码久久久久久99久播 | 日日摸夜夜添夜夜爱| 青青草视频在线视频观看| 桃色一区二区三区在线观看| 亚洲国产日韩欧美精品在线观看| 亚洲精品久久久久久婷婷小说 | 18禁动态无遮挡网站| 成人午夜精彩视频在线观看| 久久精品91蜜桃| 亚洲自拍偷在线| 美女国产视频在线观看| 一个人观看的视频www高清免费观看| 有码 亚洲区| 亚洲四区av| 老女人水多毛片| 夜夜看夜夜爽夜夜摸| 美女黄网站色视频| 日本一二三区视频观看| 国产成人a区在线观看| 成人亚洲欧美一区二区av| 亚洲综合色惰| 男人的好看免费观看在线视频| 亚洲在线观看片| 永久免费av网站大全| 精品久久国产蜜桃| 亚洲人与动物交配视频| 亚洲三级黄色毛片| 久久久亚洲精品成人影院| 久久久久久国产a免费观看| 一本一本综合久久| 国产又黄又爽又无遮挡在线| 国产伦在线观看视频一区| 中文天堂在线官网| 精品一区二区免费观看| 偷拍熟女少妇极品色| 超碰97精品在线观看| 一级毛片aaaaaa免费看小| 国产成人精品久久久久久| 高清av免费在线| 99在线视频只有这里精品首页| 欧美高清性xxxxhd video| 日韩在线高清观看一区二区三区| 国产一级毛片七仙女欲春2| 日韩大片免费观看网站 | 国产v大片淫在线免费观看| 日韩一区二区视频免费看| 精品熟女少妇av免费看| 亚洲熟妇中文字幕五十中出| 亚洲欧美日韩高清专用| 国产精品久久电影中文字幕| 亚洲人成网站高清观看| 国产精品一区二区三区四区免费观看| 亚洲av日韩在线播放| 男女边吃奶边做爰视频| 建设人人有责人人尽责人人享有的 | 欧美日韩在线观看h| 国产一区二区亚洲精品在线观看| 国产片特级美女逼逼视频| 成人一区二区视频在线观看| 亚洲电影在线观看av| АⅤ资源中文在线天堂| 国产黄色视频一区二区在线观看 | 有码 亚洲区| 成年免费大片在线观看| www.av在线官网国产| 97在线视频观看| 老师上课跳d突然被开到最大视频| 国产亚洲5aaaaa淫片| 男人的好看免费观看在线视频| 日韩,欧美,国产一区二区三区 | 黄色配什么色好看| 日本免费a在线| 日韩一区二区视频免费看| 99热网站在线观看| 中文字幕免费在线视频6| 久久久亚洲精品成人影院| 午夜免费激情av| 岛国毛片在线播放| 日韩成人av中文字幕在线观看| 国产欧美日韩精品一区二区| 亚州av有码| 99久国产av精品| 波多野结衣高清无吗| 99久国产av精品国产电影| 看十八女毛片水多多多| 最近中文字幕高清免费大全6| 可以在线观看毛片的网站| 久久99热这里只有精品18| 建设人人有责人人尽责人人享有的 | 一区二区三区乱码不卡18| 成人鲁丝片一二三区免费| 可以在线观看毛片的网站| 日本免费在线观看一区| 亚洲熟妇中文字幕五十中出| 午夜激情欧美在线| 国产午夜精品论理片| 成人国产麻豆网| 亚洲精品影视一区二区三区av| 91精品一卡2卡3卡4卡| 国产国拍精品亚洲av在线观看| 国产在线一区二区三区精 | 亚洲国产精品成人久久小说| 免费观看在线日韩| 卡戴珊不雅视频在线播放| 国产精品99久久久久久久久| 超碰av人人做人人爽久久| 久久精品熟女亚洲av麻豆精品 | 噜噜噜噜噜久久久久久91| 国产成人精品一,二区| 热99在线观看视频| 亚洲国产成人一精品久久久| 国产精品无大码| 熟女人妻精品中文字幕| 永久免费av网站大全| 最近手机中文字幕大全| 亚洲中文字幕日韩| 春色校园在线视频观看| h日本视频在线播放| 少妇裸体淫交视频免费看高清| 国产欧美日韩精品一区二区| 亚洲欧洲日产国产| 国产伦精品一区二区三区四那| 色吧在线观看| 天堂影院成人在线观看| www.色视频.com| 成年av动漫网址| av免费观看日本| 日本黄色片子视频| 蜜臀久久99精品久久宅男| 1024手机看黄色片| 国产免费福利视频在线观看| 日韩亚洲欧美综合| 永久网站在线| 2021少妇久久久久久久久久久| av天堂中文字幕网| 18禁在线无遮挡免费观看视频| 亚洲最大成人中文| 综合色丁香网| 亚洲欧美日韩东京热| 青青草视频在线视频观看| 91精品一卡2卡3卡4卡| 欧美一区二区国产精品久久精品| 国产午夜精品一二区理论片| 激情 狠狠 欧美| 成人一区二区视频在线观看| 中文精品一卡2卡3卡4更新| 亚洲天堂国产精品一区在线| 亚洲av.av天堂| 高清毛片免费看| www日本黄色视频网| 国产亚洲一区二区精品| 少妇人妻一区二区三区视频| 三级国产精品片| 国产乱人视频| 久久久久性生活片| 高清日韩中文字幕在线| 久久精品久久久久久久性| 国产亚洲5aaaaa淫片| 精品无人区乱码1区二区| 日韩人妻高清精品专区| 国产精品一区www在线观看| 亚洲成色77777| 久久精品综合一区二区三区| 99热这里只有是精品50| 99在线人妻在线中文字幕| 国产v大片淫在线免费观看| 国产成人精品久久久久久| 免费观看的影片在线观看| 99久久九九国产精品国产免费| 观看美女的网站| 亚洲无线观看免费| 日日干狠狠操夜夜爽| 特级一级黄色大片| 久久精品久久精品一区二区三区| 不卡视频在线观看欧美| 一个人观看的视频www高清免费观看| 国产精品国产高清国产av| 一二三四中文在线观看免费高清| 看十八女毛片水多多多| 能在线免费看毛片的网站| 插逼视频在线观看| 1024手机看黄色片| 日韩一本色道免费dvd| 99久国产av精品国产电影| 亚洲欧美精品自产自拍| 1024手机看黄色片| 亚洲av中文av极速乱| 热99re8久久精品国产| 婷婷色麻豆天堂久久 | 久久亚洲精品不卡| 老师上课跳d突然被开到最大视频| 七月丁香在线播放| 免费观看精品视频网站| 淫秽高清视频在线观看| 精品少妇黑人巨大在线播放 | 免费观看在线日韩| 国产精品一区www在线观看| 久久久亚洲精品成人影院| 精品无人区乱码1区二区| 久久99精品国语久久久| 欧美成人免费av一区二区三区| 免费观看a级毛片全部| 老司机福利观看| 久久亚洲精品不卡| 国产一区有黄有色的免费视频 | av在线天堂中文字幕| 麻豆成人av视频| 最近中文字幕高清免费大全6| 国产毛片a区久久久久| 日日摸夜夜添夜夜爱| 黄色日韩在线| 男人和女人高潮做爰伦理| 国产av码专区亚洲av| 亚洲最大成人手机在线| 日韩av不卡免费在线播放| 久久精品国产亚洲网站| 黄片wwwwww| av线在线观看网站| 久久久久久久久久黄片| 秋霞伦理黄片| 少妇裸体淫交视频免费看高清| 久久久久国产网址| 美女内射精品一级片tv| 亚洲av日韩在线播放| 日韩 亚洲 欧美在线| 亚洲精品日韩在线中文字幕| 国产v大片淫在线免费观看| 国内揄拍国产精品人妻在线| 久久久久久久亚洲中文字幕| 精品久久久久久久久av| 高清视频免费观看一区二区 | 亚洲国产成人一精品久久久| 七月丁香在线播放| 久久这里有精品视频免费| 中国国产av一级| 特级一级黄色大片| 欧美成人一区二区免费高清观看| 午夜福利网站1000一区二区三区| 三级国产精品片| 亚洲精品成人久久久久久| 精品久久久久久成人av| 国产欧美另类精品又又久久亚洲欧美| 赤兔流量卡办理| 国产精品精品国产色婷婷| 亚洲欧美日韩东京热| 久久午夜福利片| 日本免费在线观看一区| 久久精品夜色国产| 日韩中字成人| 中文在线观看免费www的网站| 亚洲自偷自拍三级| 亚洲人成网站在线观看播放| 在线观看美女被高潮喷水网站| 欧美不卡视频在线免费观看| 99久国产av精品国产电影| 少妇的逼水好多| 欧美性猛交╳xxx乱大交人| 看免费成人av毛片| 插阴视频在线观看视频| ponron亚洲| 欧美丝袜亚洲另类| 亚洲欧美精品专区久久| 欧美极品一区二区三区四区| 免费电影在线观看免费观看| 国产探花在线观看一区二区| 国产视频首页在线观看| 黄色配什么色好看| 日本一二三区视频观看| 亚洲精品乱码久久久久久按摩| 麻豆av噜噜一区二区三区| 最新中文字幕久久久久| 国产成人a区在线观看| 国产亚洲最大av| 成人高潮视频无遮挡免费网站| 18禁在线无遮挡免费观看视频| 身体一侧抽搐| 99热全是精品| 精品久久久久久久久亚洲| 国产大屁股一区二区在线视频| 久久人人爽人人爽人人片va| 国产又色又爽无遮挡免| 91久久精品国产一区二区成人| 九九在线视频观看精品| 狂野欧美白嫩少妇大欣赏| 深夜a级毛片| av黄色大香蕉| 色综合站精品国产| 国产毛片a区久久久久| 亚洲四区av| 亚洲真实伦在线观看| 亚洲国产精品久久男人天堂| 亚洲内射少妇av| 精品少妇黑人巨大在线播放 | 日韩av在线免费看完整版不卡| 亚洲一区高清亚洲精品| 久久婷婷人人爽人人干人人爱| 啦啦啦韩国在线观看视频| 99九九线精品视频在线观看视频| .国产精品久久| 精品欧美国产一区二区三| 欧美日本亚洲视频在线播放| 99久久成人亚洲精品观看| 国产亚洲午夜精品一区二区久久 | 亚洲成av人片在线播放无| 日本av手机在线免费观看| 国产av码专区亚洲av| 一级黄片播放器| 少妇丰满av| 欧美97在线视频| 国产精品精品国产色婷婷| 波野结衣二区三区在线| 简卡轻食公司| 看非洲黑人一级黄片| 国产色婷婷99| 免费观看性生交大片5| 日本-黄色视频高清免费观看| 建设人人有责人人尽责人人享有的 | 国产成人a∨麻豆精品| 永久免费av网站大全| 国产精品久久久久久久久免| 免费电影在线观看免费观看| 日韩成人伦理影院| 国产精品久久视频播放| 午夜福利高清视频| 亚洲人成网站在线播| 九色成人免费人妻av| 欧美+日韩+精品| 能在线免费看毛片的网站| 神马国产精品三级电影在线观看| 亚洲精品一区蜜桃| 国国产精品蜜臀av免费| 亚洲成色77777| 欧美xxxx性猛交bbbb| 少妇丰满av| 99久久精品一区二区三区| 两个人视频免费观看高清| 1000部很黄的大片| 国产v大片淫在线免费观看| 中文精品一卡2卡3卡4更新| 不卡视频在线观看欧美| 精品欧美国产一区二区三| av免费观看日本| 欧美97在线视频| av在线播放精品| 国产女主播在线喷水免费视频网站 | 中文字幕人妻熟人妻熟丝袜美| or卡值多少钱| 全区人妻精品视频| 国产精品伦人一区二区| 国产av不卡久久| 伦精品一区二区三区| 最近手机中文字幕大全| 深夜a级毛片| 日韩视频在线欧美| 欧美zozozo另类| 亚洲精品aⅴ在线观看| 女人被狂操c到高潮| 亚洲18禁久久av| 久久这里有精品视频免费| 99九九线精品视频在线观看视频| 日韩国内少妇激情av| 精品久久久久久久末码| 99热这里只有是精品在线观看| 成人二区视频| 国产精品一及| 小蜜桃在线观看免费完整版高清| 嫩草影院新地址| 日本免费在线观看一区| 亚洲欧美一区二区三区国产| 亚洲精品成人久久久久久| 99热这里只有精品一区| 久久久久久伊人网av| 天堂网av新在线| 色综合站精品国产| 搞女人的毛片| 国产成人a∨麻豆精品| 天堂中文最新版在线下载 | 三级经典国产精品| 国内精品宾馆在线| 亚洲av中文av极速乱| 蜜桃亚洲精品一区二区三区| 丰满人妻一区二区三区视频av| 中文字幕av在线有码专区| av在线播放精品| 美女脱内裤让男人舔精品视频| 简卡轻食公司| 成人三级黄色视频| 午夜久久久久精精品| 成人特级av手机在线观看| 伦理电影大哥的女人| 好男人在线观看高清免费视频| 日韩精品青青久久久久久| 内地一区二区视频在线| 国产一区亚洲一区在线观看| 久久久久免费精品人妻一区二区| 少妇猛男粗大的猛烈进出视频 | 亚洲国产精品国产精品| 亚洲av成人精品一二三区| 在线播放国产精品三级| 国产国拍精品亚洲av在线观看| 91精品一卡2卡3卡4卡| 美女内射精品一级片tv| 成人毛片60女人毛片免费| 久久人人爽人人片av| 成人毛片a级毛片在线播放| 亚洲人与动物交配视频| 国产精品无大码| 久久精品国产鲁丝片午夜精品| 国产一区亚洲一区在线观看| 久久久久久久久中文| av视频在线观看入口| 久久精品91蜜桃| 免费在线观看成人毛片| eeuss影院久久| 男女那种视频在线观看| or卡值多少钱| 国产精品.久久久| 日韩亚洲欧美综合| 人人妻人人澡欧美一区二区| 一卡2卡三卡四卡精品乱码亚洲| 国产精品一及|