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

    一種基于VMD的顫振試飛模態(tài)參數(shù)辨識(shí)方法

    2020-02-04 07:27:44王紹楠盧曉東
    航空科學(xué)技術(shù) 2020年9期
    關(guān)鍵詞:信噪比

    王紹楠 盧曉東

    摘要:顫振試飛信號(hào)具有非平穩(wěn)、模態(tài)密集、信噪比(SNR)低的特點(diǎn)。近年來(lái),新發(fā)展的變分模態(tài)分解方法在處理該類(lèi)信號(hào)時(shí)能有效避免模態(tài)混疊,且有較好抗噪能力。故將變分模態(tài)分解法和希爾伯特變換結(jié)合,提出了一種可用于顫振試飛脈沖響應(yīng)的時(shí)頻域模態(tài)參數(shù)辨識(shí)新方法,稱(chēng)為VMD-HT方法。通過(guò)仿真算例驗(yàn)證和實(shí)際顫振試飛中的應(yīng)用,表明該方法能準(zhǔn)確穩(wěn)定地辨識(shí)出脈沖信號(hào)的模態(tài)參數(shù),且對(duì)于不同部位傳感器數(shù)據(jù)的辨識(shí)結(jié)果聚集度更高,在實(shí)際型號(hào)試飛中取得了良好效果。

    關(guān)鍵詞:信噪比;顫振試飛;變分模態(tài)分解;希爾伯特變換;脈沖信號(hào)

    中圖分類(lèi)號(hào):V217文獻(xiàn)標(biāo)識(shí)碼:ADOI:10.19452/j.issn1007-5453.2020.09.008

    顫振試飛是世界公認(rèn)的一類(lèi)風(fēng)險(xiǎn)試飛科目,所有新機(jī)及結(jié)構(gòu)有重大改變的飛機(jī)都必須進(jìn)行顫振試飛[1-2]。其目的是為了獲取飛機(jī)結(jié)構(gòu)的顫振特性,包括模態(tài)參數(shù)(頻率和阻尼)、顫振類(lèi)型以及顫振邊界。而模態(tài)參數(shù)的準(zhǔn)確辨識(shí)是獲得顫振速度的關(guān)鍵。但顫振試飛數(shù)據(jù)非平穩(wěn)、模態(tài)密集、信噪比低的特點(diǎn)往往給模態(tài)參數(shù)的準(zhǔn)確辨識(shí)帶來(lái)難度[3]。目前常用的模態(tài)辨識(shí)方法有頻域法和時(shí)域法,兩者分別從頻域和時(shí)域的角度出發(fā),對(duì)系統(tǒng)進(jìn)行模態(tài)參數(shù)辨識(shí)。隨著模態(tài)辨識(shí)技術(shù)的發(fā)展,時(shí)頻聯(lián)合分析的時(shí)頻域模態(tài)辨識(shí)方法應(yīng)運(yùn)而生,其綜合了頻域和時(shí)域辨識(shí)的優(yōu)點(diǎn),逐步成為現(xiàn)代試驗(yàn)中廣泛采用的模態(tài)辨識(shí)技術(shù)[4]。

    希爾伯特-黃變換(hilbert-huang transform , HHT)是由Huang等提出的一種分析非平穩(wěn)信號(hào)的時(shí)頻聯(lián)合分析方法,由經(jīng)驗(yàn)?zāi)B(tài)分解(empirical mode decomposition, EMD)和希爾伯特變換組成[5]。而后,有不少將HHT應(yīng)用于線性系統(tǒng)模態(tài)參數(shù)辨識(shí)的研究[6]。HHT雖然是一種高效、準(zhǔn)確的辨識(shí)方法,但該方法實(shí)踐超前于理論,目前尚沒(méi)有嚴(yán)格的理論支持。其核心EMD屬于遞歸式的模態(tài)分解,采用了基于極值點(diǎn)的包絡(luò)求取方式,其包絡(luò)估計(jì)誤差經(jīng)多次遞歸分解而被放大,易出現(xiàn)模態(tài)混疊。雖然采用EEMD能對(duì)模態(tài)混疊現(xiàn)象進(jìn)行抑制,但這將大幅增加計(jì)算量。同時(shí)也存在采樣效應(yīng)和端點(diǎn)效應(yīng)的問(wèn)題。

    為了克服以上問(wèn)題,2014年,Dragomiretskiy等提出了一種完全不同于遞歸式模態(tài)分解的新方法,即變分模態(tài)分解(variational mode decomposition , VMD)[7]。其整體框架是變分問(wèn)題的求解,具有堅(jiān)實(shí)的理論基礎(chǔ)。VMD的實(shí)質(zhì)是多個(gè)自適應(yīng)維納濾波器,表現(xiàn)出更好的噪聲魯棒性。在模態(tài)分離方面,VMD可將頻率相近的兩個(gè)純諧波信號(hào)成功分離。另外,VMD的端點(diǎn)效應(yīng)也遠(yuǎn)弱于遞歸式模態(tài)分解。VMD一經(jīng)提出,即成為眾多研究者研究的熱點(diǎn)。目前,該方法已用于機(jī)械、電子、生物、能源等領(lǐng)域,尤其在機(jī)械故障診斷中應(yīng)用最廣[8-11]。眾多研究表明,VMD方法能有效避免模態(tài)混疊,有較好噪聲魯棒性,在處理信噪比低、低頻、非平穩(wěn)信號(hào)方面有很大優(yōu)勢(shì)。鑒于此,本文首次將VMD引入顫振試飛領(lǐng)域,將其和Hilbert變換結(jié)合,提出了一種可用于脈沖激勵(lì)下的顫振試飛模態(tài)參數(shù)辨識(shí)方法,稱(chēng)為VMD-HT法,并通過(guò)仿真算例和實(shí)際飛行數(shù)據(jù)驗(yàn)證了該方法的正確性、有效性和噪聲魯棒性。

    1理論介紹

    1.1 VMD理論

    綜上,VMD-HT方法進(jìn)行模態(tài)參數(shù)辨識(shí)過(guò)程可簡(jiǎn)述為:(1)對(duì)測(cè)得的脈沖響應(yīng)信號(hào)進(jìn)行頻譜分析,確定大致的模態(tài)數(shù)k;(2)通過(guò)VMD計(jì)算,將原始信號(hào)分解為k個(gè)單模態(tài)信號(hào)IMF;(3)對(duì)IMF進(jìn)行Hilbert變換,構(gòu)造解析函數(shù),引出瞬時(shí)幅值與瞬時(shí)相角;分別取對(duì)數(shù),并對(duì)數(shù)據(jù)進(jìn)行最小二乘擬合,得到相應(yīng)斜率;根據(jù)斜率計(jì)算得到模態(tài)頻率和阻尼比;(4)對(duì)每一個(gè)IMF均進(jìn)行步驟(3),即得出系統(tǒng)所有模態(tài)的頻率和阻尼比。

    2仿真算例

    為驗(yàn)證本文方法的正確性,首先構(gòu)造一個(gè)含三階模態(tài)的脈沖響應(yīng)信號(hào),為三個(gè)衰減的正弦函數(shù)的疊加,脈沖響應(yīng)形式如式(4)所述,并加入服從高斯分布的噪聲信號(hào)η,η~N(0,σ)。其中,σ為偏差量,用來(lái)控制噪聲強(qiáng)弱,這里σ= 0.2。

    圖1為該信號(hào)的時(shí)域波形。再采用VMD-HT方法辨識(shí)模態(tài)參數(shù)。頻譜分析顯示該信號(hào)有三個(gè)主要模態(tài),故模態(tài)數(shù)取3。圖2為VMD分解信號(hào)與原始信號(hào)的對(duì)比。圖3為VMD分解信號(hào)的頻譜與原始信號(hào)的對(duì)比,頻率吻合良好。由圖2可見(jiàn)也存在一定端點(diǎn)效應(yīng),為提高辨識(shí)精度,在進(jìn)行數(shù)據(jù)擬合時(shí),只選擇中部成明顯線性段的數(shù)據(jù)進(jìn)行線性擬合。圖4為第一階模態(tài)擬合曲線。

    表1中給出了模態(tài)頻率和阻尼比的辨識(shí)結(jié)果與理論值的對(duì)比,兩者吻合良好。由此,驗(yàn)證了VMD-HT方法對(duì)脈沖響應(yīng)信號(hào)模態(tài)辨識(shí)的正確性和有效性。

    此外,研究了該方法對(duì)噪聲的敏感性,分別選取σ= 0~0.7。因頻率辨識(shí)結(jié)果與理論值吻合一直很好,這里僅展示阻尼比的相對(duì)誤差(取絕對(duì)值)隨噪聲水平的變化,如圖5所示。可見(jiàn),三個(gè)模態(tài)的阻尼比相對(duì)誤差整體趨勢(shì)是隨噪聲的增加而增加。σ在0.5以?xún)?nèi)所有模態(tài)的阻尼辨識(shí)誤差不超過(guò)10%,頻率辨識(shí)誤差不超過(guò)0.2%;可見(jiàn)VMD-HT方法有較強(qiáng)的噪聲魯棒性。

    3試飛應(yīng)用

    某型飛機(jī)在顫振試飛中對(duì)飛機(jī)機(jī)翼多部位施加脈沖激勵(lì)。分別采用VMD-HT方法和顫振試飛處理脈沖數(shù)據(jù)的專(zhuān)用軟件對(duì)相同的激勵(lì)數(shù)據(jù)進(jìn)行模態(tài)參數(shù)辨識(shí)。

    圖6為機(jī)翼上某個(gè)傳感器測(cè)得的時(shí)域數(shù)據(jù)。首先,設(shè)置模態(tài)數(shù)為5,采用VMD-HT方法對(duì)該數(shù)據(jù)進(jìn)行計(jì)算。圖7為采用VMD分解出的信號(hào),可見(jiàn)符合IMF要求是u2~u5,即VMD成功分解出了其中的4階主要模態(tài)。由圖8為VMD分解出的信號(hào)頻譜與原信號(hào)頻譜的對(duì)比,可見(jiàn)有4階吻合良好。圖9展示了其中模態(tài)2(u2)的曲線擬合結(jié)果。

    對(duì)其他部位傳感器的數(shù)據(jù)進(jìn)行分析,VMD-HT方法的辨識(shí)結(jié)果與顫振專(zhuān)用軟件的辨識(shí)結(jié)果對(duì)比如圖10所示,從左到右分別為1~5階主要模態(tài)。可見(jiàn),VMD-HT方法能辨識(shí)出顫振所關(guān)心的結(jié)構(gòu)主要模態(tài),頻率與專(zhuān)用軟件結(jié)果吻合良好,阻尼比的辨識(shí)較為保守。圖11展示了這兩種方法辨識(shí)的模態(tài)2頻率和阻尼比隨傳感器位置的變化。可見(jiàn),專(zhuān)用軟件的傳統(tǒng)方法對(duì)不同傳感器數(shù)據(jù)辨識(shí)時(shí),同一模態(tài)的結(jié)果有較大差異,有時(shí)需采用取平均的辦法來(lái)確定該模態(tài)的頻率和阻尼比。而采用VMD-HT方法計(jì)算時(shí),頻率和阻尼比的辨識(shí)結(jié)果不會(huì)隨傳感器位置的變化有較大差異,結(jié)果聚集度高,更加穩(wěn)定,在實(shí)際顫振試飛數(shù)據(jù)的模態(tài)參數(shù)辨識(shí)中取得了良好效果。

    4結(jié)論

    通過(guò)研究,可以得出如下結(jié)論:

    (1)仿真算例驗(yàn)證了本文發(fā)展的VMD-HT顫振模態(tài)參數(shù)辨識(shí)方法能正確有效地辨識(shí)出脈沖響應(yīng)的頻率和阻尼,具有較好的噪聲魯棒性。

    (2)在實(shí)際顫振試飛中應(yīng)用該方法,其模態(tài)參數(shù)辨識(shí)結(jié)果不會(huì)隨傳感器位置的變化有較大差異,結(jié)果穩(wěn)定有效,可成為型號(hào)試飛中顫振試飛數(shù)據(jù)后處理分析的有力工具。

    (3)該方法的不足之處是VMD算法需要提前設(shè)置模態(tài)數(shù)。本文通過(guò)對(duì)信號(hào)先進(jìn)行頻譜分析來(lái)確定模態(tài)數(shù),有利于后續(xù)準(zhǔn)確有效地分解出IMF,是VMD-HT方法中行之有效的一步。

    參考文獻(xiàn)

    [1]張偉偉,鐘華壽,肖華,等.顫振飛行試驗(yàn)的邊界預(yù)測(cè)方法回顧與展望[J].航空學(xué)報(bào),2015,36(5): 1367-1384. Zhang Weiwei, Zhong Huashou, Xiao Hua,et al. Review and prospect of flutter boundary prediction methods for flight flutter testing[J]. Acta Aeronautica et Astronautica Sinica,2015,36(5):1367-1384.(in Chinese)

    [2]劉立坤,張春蔚,王東森.基于時(shí)序模型的顫振預(yù)測(cè)技術(shù)研究[J].航空科學(xué)技術(shù),2018,29(11):20-24. Liu Likun, Zhang Chunwei, Wang Dongsen. Flutter boundary prediction based on time series models[J]. Aeronautical Science & Technology, 2018, 29(11):20-24.(in Chinese)

    [3]秦浩,周友明.改進(jìn)頻域子空間模態(tài)分析方法及在試飛中的應(yīng)用[J].航空科學(xué)技術(shù),2017,28(5):59-62. Qin Hao, Zhou Youming. An improved modal parameter identification method with optimized frequency subspace and application in flight flutter test[J]. Aeronautical Science & Technology, 2017,28(5):59-62.(in Chinese)

    [4]李晶,曹登慶,劉紹奎,等.基于解析模態(tài)分解和希爾伯特變換的模態(tài)參數(shù)辨識(shí)新方法[J].振動(dòng)與沖擊,2016,35(1):34-39. Li Jing, Cao Dengqing, Liu Shaokui, et al. A new method for modal parameter identification based on analytical model decompositionandHilberttransformation[J]. Journalof Vibration and Shock, 2016,35(1):34-39.(in Chinese)

    [5]Huang N E,Shen Z,Long S R,et al. The empirical mode decomposition and the Hilbert spectrum for nonlinear and nonstationary time series analysis[J]. Proceedings of The Royal Society A:Mathematical,Physical and Engineering Sciences,1998,454(1971):903-995.

    [6]Yang J N,Lei Y,Pan S,et al. System identification of linear structures based on Hilbert-Huang spectral analysis. Part 1:normalmodes[J].EarthquakeEngineering&Structural Dynamics,2003,32(9):1443-1467.

    [7]Dragomiretskiy K,Zosso D. Variational mode decomposition[J]. IEEE Transactions on Signal Processing,2014,62(3):531-544.

    [8]武英杰,甄成剛,劉長(zhǎng)良.變分模態(tài)分解在風(fēng)電機(jī)組故障診斷中的應(yīng)用[J].機(jī)械傳動(dòng),2015,39(10):129-132. Wu Yingjie, Zhen Chenggang, Liu Changliang. Application of variational mode decomposition in wind power fault diagnosis[J]. Journal of Mechanical Transmission,2015, 39(10):129-132.(in Chinese)

    [9]唐貴基,王曉龍.參數(shù)優(yōu)化變分模態(tài)分解方法在滾動(dòng)軸承早期故障診斷中的應(yīng)用[J].西安交通大學(xué)學(xué)報(bào),2015,49(5):73-81. Tang Guiji, Wang Xiaolong. Parameter optimized variational mode decomposition method with application to incipient fault diagnosis of rolling bearing[J]. Journal of Xian Jiaotong University,2015,49(5):73-81.(in Chinese)

    [10]趙昕海,張術(shù)臣,李志深,等.基于VMD的故障特征信號(hào)提取方法[J].振動(dòng)測(cè)試與診斷,2018, 38(1):11-19. Zhao Xinhai, Zhang Shuchen, Li Zhishen, et al. Application of new denoising method based on VMD in fault feature extraction[J]. Journal of Vibration,Measurement&Diagnosis, 2018, 38(1):11-19.(in Chinese)

    [11]劉長(zhǎng)福,朱立達(dá),仇健,等.基于VMD和FFT的變切深側(cè)銑顫振特征提取方法[J].東北大學(xué)學(xué)報(bào)(自然科學(xué)版),2018,39(8):1153-1157. Liu Changfu, Zhu Lida, Qiu Jian, et al. Chatter feature extraction method in variable cutting depth flank milling based on VMD and FFT[J]. Journal of Northeastern University(Natural Science), 2018,39(8):1153-1157.(in Chinese)

    (責(zé)任編輯陳東曉)

    作者簡(jiǎn)介

    王紹楠(1992-)女,碩士,工程師。主要研究方向:飛行器結(jié)構(gòu)動(dòng)力學(xué)試飛。

    Tel:029-86838039E-mail:342669853@qq.com

    A Modal Parameter Identification Method for Flutter Flight Test Based on VMD

    Wang Shaonan*,Lu Xiaodong

    Chinese Flight Test Establishment,Xian 710089,China

    Abstract: The signal of flutter flight test has the characteristics of non-stationary, dense mode and low SNR. Variational mode decomposition method, which was developed in recent years, can avoid mode mixing effectively and have a good resistant ability to noise. Therefore, by combining variational mode decomposition method and Hilbert transform, a new modal parameter identification method of time-frequency domain for impulse response in flutter flight test was proposed in this paper, which is called VMD-HT method. The simulation examples and the application in the actual flutter flight test showed that this method can accurately and stably identify the modal parameters of the impulse signal, and the identification results for the sensor data in different positions have higher aggregation,which has achieved good results in the actual flutter flight test.

    Key Words: SNR; flutter flight test; variational mode decomposition; Hilbert transform; impulse signal

    猜你喜歡
    信噪比
    兩種64排GE CT冠脈成像信噪比與劑量對(duì)比分析研究
    基于經(jīng)驗(yàn)分布函數(shù)快速收斂的信噪比估計(jì)器
    一種基于2G-ALE中快速信噪比的估計(jì)算法
    無(wú)線通信中的信噪比估計(jì)算法研究
    信噪比在AR模型定階方法選擇中的研究
    自跟蹤接收機(jī)互相關(guān)法性能分析
    基于深度學(xué)習(xí)的無(wú)人機(jī)數(shù)據(jù)鏈信噪比估計(jì)算法
    低信噪比下LFMCW信號(hào)調(diào)頻參數(shù)估計(jì)
    低信噪比下基于Hough變換的前視陣列SAR稀疏三維成像
    不同信噪比下的被動(dòng)相控陣?yán)走_(dá)比幅測(cè)角方法研究
    美女大奶头视频| 一区二区三区精品91| 一进一出抽搐动态| 欧洲精品卡2卡3卡4卡5卡区| 老司机深夜福利视频在线观看| 国产亚洲精品久久久久久毛片| 人人妻,人人澡人人爽秒播| 国产蜜桃级精品一区二区三区| 露出奶头的视频| 欧美乱妇无乱码| 午夜福利成人在线免费观看| 欧美日韩瑟瑟在线播放| 午夜影院日韩av| 久久中文看片网| 90打野战视频偷拍视频| 18禁美女被吸乳视频| av视频在线观看入口| 一边摸一边抽搐一进一小说| 久久九九热精品免费| 国产色视频综合| 亚洲一卡2卡3卡4卡5卡精品中文| 免费搜索国产男女视频| 在线观看免费日韩欧美大片| 69精品国产乱码久久久| 99精品久久久久人妻精品| 欧美成狂野欧美在线观看| 久久性视频一级片| 欧美日韩福利视频一区二区| 国内精品久久久久久久电影| 久久草成人影院| 日本在线视频免费播放| 精品第一国产精品| 国产精品野战在线观看| 国产熟女xx| 手机成人av网站| 少妇粗大呻吟视频| 日本五十路高清| 老汉色av国产亚洲站长工具| 男人操女人黄网站| 黄色a级毛片大全视频| 侵犯人妻中文字幕一二三四区| 欧美+亚洲+日韩+国产| 亚洲 国产 在线| 91字幕亚洲| 国产在线观看jvid| 精品国产美女av久久久久小说| 欧美激情极品国产一区二区三区| 99在线视频只有这里精品首页| 亚洲欧洲精品一区二区精品久久久| 久久久久国产精品人妻aⅴ院| 涩涩av久久男人的天堂| 叶爱在线成人免费视频播放| 亚洲成国产人片在线观看| 亚洲av片天天在线观看| 97碰自拍视频| av在线天堂中文字幕| 非洲黑人性xxxx精品又粗又长| 免费观看精品视频网站| 12—13女人毛片做爰片一| 亚洲男人天堂网一区| 麻豆av在线久日| 亚洲免费av在线视频| 午夜久久久在线观看| 老司机午夜福利在线观看视频| 人妻久久中文字幕网| 色尼玛亚洲综合影院| 脱女人内裤的视频| 脱女人内裤的视频| 色综合婷婷激情| 9热在线视频观看99| 黄网站色视频无遮挡免费观看| 大型黄色视频在线免费观看| 美女午夜性视频免费| 在线观看免费视频网站a站| 欧美黑人欧美精品刺激| 黄色视频,在线免费观看| 亚洲av第一区精品v没综合| 一边摸一边抽搐一进一小说| 在线观看日韩欧美| 日本五十路高清| 午夜福利免费观看在线| 极品人妻少妇av视频| 真人做人爱边吃奶动态| 97碰自拍视频| 久久久久久大精品| 久久久国产成人免费| 精品日产1卡2卡| 国产成人av激情在线播放| 色综合欧美亚洲国产小说| 国产精品 欧美亚洲| 夜夜躁狠狠躁天天躁| 中文字幕久久专区| 成人永久免费在线观看视频| 亚洲午夜精品一区,二区,三区| 丝袜在线中文字幕| 久久国产乱子伦精品免费另类| 国产精品久久久人人做人人爽| 黄片小视频在线播放| 国产亚洲欧美在线一区二区| 大型av网站在线播放| 大香蕉久久成人网| aaaaa片日本免费| 成人18禁在线播放| av电影中文网址| 精品国产亚洲在线| 午夜精品国产一区二区电影| 亚洲在线自拍视频| 女性被躁到高潮视频| 精品电影一区二区在线| 十八禁网站免费在线| 免费观看人在逋| 国产色视频综合| 欧美一级毛片孕妇| 精品卡一卡二卡四卡免费| 免费在线观看影片大全网站| 88av欧美| 久久久久久久久久久久大奶| 免费在线观看亚洲国产| 黄色毛片三级朝国网站| 国产视频一区二区在线看| 日韩 欧美 亚洲 中文字幕| 日韩欧美免费精品| 国产成人影院久久av| 成人三级做爰电影| 国产伦一二天堂av在线观看| 女人被躁到高潮嗷嗷叫费观| 乱人伦中国视频| 国产一区二区三区综合在线观看| 十八禁人妻一区二区| 成人手机av| 极品人妻少妇av视频| 黄色丝袜av网址大全| 亚洲av成人一区二区三| 国产精品免费视频内射| 国产野战对白在线观看| 变态另类丝袜制服| 免费搜索国产男女视频| 国产成人一区二区三区免费视频网站| 中文字幕最新亚洲高清| 国产免费男女视频| 久久久久久人人人人人| 久久精品国产亚洲av香蕉五月| 黑丝袜美女国产一区| 如日韩欧美国产精品一区二区三区| 自拍欧美九色日韩亚洲蝌蚪91| 亚洲午夜理论影院| 久久国产精品男人的天堂亚洲| 亚洲av熟女| 在线播放国产精品三级| 老司机靠b影院| 午夜福利免费观看在线| 国产亚洲欧美精品永久| 妹子高潮喷水视频| 99国产综合亚洲精品| 人妻久久中文字幕网| 亚洲精品一卡2卡三卡4卡5卡| 精品乱码久久久久久99久播| 国产极品粉嫩免费观看在线| 女性被躁到高潮视频| 丰满的人妻完整版| 69精品国产乱码久久久| 欧美日本视频| 中文字幕人成人乱码亚洲影| 麻豆av在线久日| 高潮久久久久久久久久久不卡| 国产一区二区三区视频了| 成人免费观看视频高清| 欧美另类亚洲清纯唯美| 欧美亚洲日本最大视频资源| 黄频高清免费视频| 正在播放国产对白刺激| 国产欧美日韩一区二区三区在线| 满18在线观看网站| 日本三级黄在线观看| 天天一区二区日本电影三级 | 欧美黄色片欧美黄色片| 欧美在线一区亚洲| 99香蕉大伊视频| 午夜久久久久精精品| 99国产精品一区二区三区| 亚洲国产高清在线一区二区三 | 久久精品91蜜桃| 一区福利在线观看| 久久中文字幕人妻熟女| 制服诱惑二区| 欧美国产日韩亚洲一区| 亚洲最大成人中文| 亚洲成人精品中文字幕电影| 午夜福利免费观看在线| 法律面前人人平等表现在哪些方面| 日韩中文字幕欧美一区二区| 丝袜美腿诱惑在线| 中文字幕最新亚洲高清| 国产三级在线视频| 亚洲aⅴ乱码一区二区在线播放 | 黑人欧美特级aaaaaa片| 一级黄色大片毛片| 国产在线观看jvid| 国产欧美日韩一区二区三| 欧美日韩瑟瑟在线播放| 国产99白浆流出| 一区福利在线观看| 91成年电影在线观看| a级毛片在线看网站| 老司机午夜福利在线观看视频| 亚洲专区字幕在线| 此物有八面人人有两片| 天堂√8在线中文| 在线观看免费日韩欧美大片| 伦理电影免费视频| 久热爱精品视频在线9| 在线观看一区二区三区| 欧美激情高清一区二区三区| 无人区码免费观看不卡| 9色porny在线观看| 亚洲欧美日韩另类电影网站| 精品国产一区二区三区四区第35| 精品免费久久久久久久清纯| 免费在线观看亚洲国产| 午夜日韩欧美国产| av中文乱码字幕在线| 久久久久久国产a免费观看| 欧美人与性动交α欧美精品济南到| 美女免费视频网站| 天天躁夜夜躁狠狠躁躁| 久久精品国产亚洲av香蕉五月| av有码第一页| 在线观看日韩欧美| 亚洲色图综合在线观看| 精品国产超薄肉色丝袜足j| 在线免费观看的www视频| 在线观看午夜福利视频| 校园春色视频在线观看| 免费观看精品视频网站| 每晚都被弄得嗷嗷叫到高潮| 深夜精品福利| 国产色视频综合| 91老司机精品| 亚洲精品在线美女| 午夜福利,免费看| 大香蕉久久成人网| 欧美久久黑人一区二区| 国产一级毛片七仙女欲春2 | 高清毛片免费观看视频网站| 免费在线观看视频国产中文字幕亚洲| 国产午夜精品久久久久久| 亚洲第一欧美日韩一区二区三区| 身体一侧抽搐| 18美女黄网站色大片免费观看| 男男h啪啪无遮挡| 日韩欧美一区视频在线观看| 99在线视频只有这里精品首页| 国产精品国产高清国产av| 欧美日本中文国产一区发布| 国产精品 欧美亚洲| 波多野结衣高清无吗| 一级,二级,三级黄色视频| 国产高清videossex| 韩国av一区二区三区四区| 一a级毛片在线观看| 巨乳人妻的诱惑在线观看| 十八禁网站免费在线| 十八禁人妻一区二区| 亚洲欧美精品综合久久99| 欧美日本视频| 午夜福利18| 欧美日韩亚洲综合一区二区三区_| x7x7x7水蜜桃| 午夜日韩欧美国产| 乱人伦中国视频| 免费女性裸体啪啪无遮挡网站| 亚洲成国产人片在线观看| 亚洲自偷自拍图片 自拍| 国产在线观看jvid| а√天堂www在线а√下载| 午夜福利成人在线免费观看| 岛国在线观看网站| av视频免费观看在线观看| 桃红色精品国产亚洲av| 91老司机精品| 黑人巨大精品欧美一区二区蜜桃| 国产精品精品国产色婷婷| 国产国语露脸激情在线看| 亚洲av片天天在线观看| 男女午夜视频在线观看| 国产精品亚洲一级av第二区| 一进一出抽搐gif免费好疼| 人人妻人人爽人人添夜夜欢视频| 欧美日韩一级在线毛片| 久久人妻av系列| 精品电影一区二区在线| 日本精品一区二区三区蜜桃| 日韩有码中文字幕| 搞女人的毛片| 日韩欧美一区视频在线观看| 久久天躁狠狠躁夜夜2o2o| 亚洲av电影不卡..在线观看| 亚洲精品av麻豆狂野| 国产真人三级小视频在线观看| 欧美成狂野欧美在线观看| 中文字幕av电影在线播放| 亚洲国产毛片av蜜桃av| 日日干狠狠操夜夜爽| 成人18禁高潮啪啪吃奶动态图| 国产男靠女视频免费网站| 黄频高清免费视频| 久久久久久久久免费视频了| 无遮挡黄片免费观看| 国产精华一区二区三区| 搡老熟女国产l中国老女人| 亚洲成人国产一区在线观看| 欧美色视频一区免费| 久久精品亚洲精品国产色婷小说| 宅男免费午夜| 国产麻豆69| av免费在线观看网站| 国产亚洲欧美在线一区二区| 99久久综合精品五月天人人| 999久久久精品免费观看国产| 日韩视频一区二区在线观看| 欧美成人免费av一区二区三区| 久久精品91蜜桃| av有码第一页| 国产av一区在线观看免费| 一区二区三区高清视频在线| 国产av精品麻豆| 在线播放国产精品三级| 日本免费一区二区三区高清不卡 | 黄色视频,在线免费观看| 在线播放国产精品三级| 大型黄色视频在线免费观看| 满18在线观看网站| 我的亚洲天堂| www.自偷自拍.com| 99国产精品一区二区蜜桃av| 久久人人爽av亚洲精品天堂| 久久久久亚洲av毛片大全| 色精品久久人妻99蜜桃| 9色porny在线观看| 亚洲自拍偷在线| 欧美乱妇无乱码| 熟妇人妻久久中文字幕3abv| 国产精品免费一区二区三区在线| 在线av久久热| 妹子高潮喷水视频| 国产高清有码在线观看视频 | 成人免费观看视频高清| 国产片内射在线| 国产av在哪里看| 国产成人影院久久av| 亚洲一区中文字幕在线| av有码第一页| 国产av一区在线观看免费| 久久国产乱子伦精品免费另类| 精品欧美国产一区二区三| netflix在线观看网站| 亚洲久久久国产精品| 性欧美人与动物交配| 久久婷婷成人综合色麻豆| 最近最新免费中文字幕在线| 好看av亚洲va欧美ⅴa在| 国产成年人精品一区二区| 精品日产1卡2卡| 亚洲精品在线观看二区| 成年人黄色毛片网站| 亚洲自偷自拍图片 自拍| 国产精品爽爽va在线观看网站 | 多毛熟女@视频| 欧美+亚洲+日韩+国产| 亚洲无线在线观看| 成人18禁高潮啪啪吃奶动态图| 亚洲国产精品999在线| 99riav亚洲国产免费| 男男h啪啪无遮挡| 久久精品国产清高在天天线| 中出人妻视频一区二区| 侵犯人妻中文字幕一二三四区| 99国产精品免费福利视频| 精品国产超薄肉色丝袜足j| 精品国产美女av久久久久小说| x7x7x7水蜜桃| 亚洲一区二区三区不卡视频| 亚洲专区中文字幕在线| 两个人视频免费观看高清| 久久影院123| 欧美一区二区精品小视频在线| 国产一区二区三区在线臀色熟女| 国内精品久久久久久久电影| 国产成人精品久久二区二区免费| 国产野战对白在线观看| 亚洲精品中文字幕一二三四区| 亚洲国产中文字幕在线视频| 久久久久久人人人人人| 黑丝袜美女国产一区| 欧美色视频一区免费| 欧美最黄视频在线播放免费| 国产精品1区2区在线观看.| 亚洲精品一区av在线观看| 国产精品自产拍在线观看55亚洲| 久久精品91无色码中文字幕| 欧美国产日韩亚洲一区| 女人高潮潮喷娇喘18禁视频| 露出奶头的视频| 免费在线观看亚洲国产| 好看av亚洲va欧美ⅴa在| 精品久久久久久,| av网站免费在线观看视频| 亚洲avbb在线观看| 久热爱精品视频在线9| 日韩欧美三级三区| 69精品国产乱码久久久| 美国免费a级毛片| 波多野结衣av一区二区av| 后天国语完整版免费观看| 亚洲中文日韩欧美视频| 男女午夜视频在线观看| 亚洲精华国产精华精| 中文字幕最新亚洲高清| 国产伦一二天堂av在线观看| 久久精品国产亚洲av香蕉五月| 国产色视频综合| 免费一级毛片在线播放高清视频 | 国产一区在线观看成人免费| 亚洲第一青青草原| 99精品久久久久人妻精品| 窝窝影院91人妻| 精品一区二区三区视频在线观看免费| 日韩 欧美 亚洲 中文字幕| 免费搜索国产男女视频| 日本三级黄在线观看| 欧美日韩瑟瑟在线播放| 亚洲欧美日韩无卡精品| 免费人成视频x8x8入口观看| 欧美另类亚洲清纯唯美| 一边摸一边做爽爽视频免费| 国产xxxxx性猛交| 一区二区三区精品91| 国内精品久久久久精免费| 国产欧美日韩一区二区三区在线| 淫秽高清视频在线观看| 法律面前人人平等表现在哪些方面| 国产99白浆流出| 18禁美女被吸乳视频| 脱女人内裤的视频| 久久欧美精品欧美久久欧美| 亚洲九九香蕉| 可以在线观看的亚洲视频| 岛国在线观看网站| 女生性感内裤真人,穿戴方法视频| 女人被狂操c到高潮| 午夜福利高清视频| 丝袜人妻中文字幕| 欧美日韩瑟瑟在线播放| 日本精品一区二区三区蜜桃| 国产极品粉嫩免费观看在线| 又黄又爽又免费观看的视频| 亚洲中文日韩欧美视频| 非洲黑人性xxxx精品又粗又长| 精品一区二区三区av网在线观看| 999精品在线视频| 亚洲性夜色夜夜综合| 久久人妻av系列| 999精品在线视频| 亚洲精品国产精品久久久不卡| e午夜精品久久久久久久| 亚洲一区中文字幕在线| 99热只有精品国产| bbb黄色大片| 久久国产精品影院| 亚洲av成人av| 久久久久国产一级毛片高清牌| 一区二区三区激情视频| 亚洲av成人av| 国产欧美日韩一区二区精品| 午夜福利影视在线免费观看| 午夜免费观看网址| 日本三级黄在线观看| 麻豆久久精品国产亚洲av| 亚洲精品中文字幕一二三四区| 黑人欧美特级aaaaaa片| 97人妻精品一区二区三区麻豆 | 久久香蕉精品热| 国产欧美日韩精品亚洲av| 国产亚洲欧美98| 国产精品免费视频内射| 日韩 欧美 亚洲 中文字幕| 日韩av在线大香蕉| 中文字幕最新亚洲高清| 国产三级在线视频| 久久天躁狠狠躁夜夜2o2o| 国产aⅴ精品一区二区三区波| 日韩免费av在线播放| av视频在线观看入口| 黄色毛片三级朝国网站| 岛国视频午夜一区免费看| 电影成人av| 中文字幕高清在线视频| 少妇被粗大的猛进出69影院| 很黄的视频免费| 老熟妇乱子伦视频在线观看| 精品乱码久久久久久99久播| 亚洲一区高清亚洲精品| 一本大道久久a久久精品| 在线观看免费视频网站a站| 成年人黄色毛片网站| 一级毛片高清免费大全| 国产精品影院久久| 母亲3免费完整高清在线观看| 免费在线观看日本一区| 午夜福利18| 给我免费播放毛片高清在线观看| 久久人妻av系列| 麻豆久久精品国产亚洲av| av超薄肉色丝袜交足视频| 黑人操中国人逼视频| 宅男免费午夜| 精品日产1卡2卡| 亚洲国产精品合色在线| 精品久久久久久久毛片微露脸| 国产av又大| 97超级碰碰碰精品色视频在线观看| 欧美绝顶高潮抽搐喷水| 国产91精品成人一区二区三区| 欧美色欧美亚洲另类二区 | 精品电影一区二区在线| 国产av又大| 中文字幕精品免费在线观看视频| 亚洲激情在线av| 成年女人毛片免费观看观看9| 又黄又爽又免费观看的视频| 日本a在线网址| 黄色a级毛片大全视频| 女人高潮潮喷娇喘18禁视频| 一夜夜www| 亚洲精品国产色婷婷电影| 侵犯人妻中文字幕一二三四区| 亚洲,欧美精品.| 精品无人区乱码1区二区| 亚洲国产欧美网| 久久久久久久午夜电影| 一卡2卡三卡四卡精品乱码亚洲| 久久久久久人人人人人| a级毛片在线看网站| 人妻丰满熟妇av一区二区三区| 男女做爰动态图高潮gif福利片 | 午夜激情av网站| 久久香蕉国产精品| 伦理电影免费视频| 国产欧美日韩一区二区三| 国产区一区二久久| 午夜两性在线视频| 亚洲无线在线观看| 国产精品久久电影中文字幕| 亚洲片人在线观看| 亚洲成国产人片在线观看| 国产欧美日韩一区二区三区在线| 国产精华一区二区三区| 日韩中文字幕欧美一区二区| or卡值多少钱| 国产激情久久老熟女| 欧美+亚洲+日韩+国产| 高清黄色对白视频在线免费看| 中文字幕高清在线视频| 久久久久久大精品| 亚洲欧美日韩高清在线视频| 男人舔女人下体高潮全视频| 亚洲欧美激情综合另类| 此物有八面人人有两片| 大香蕉久久成人网| 亚洲五月色婷婷综合| 欧美在线一区亚洲| 亚洲五月色婷婷综合| 99国产精品一区二区蜜桃av| 搞女人的毛片| cao死你这个sao货| 国产亚洲精品一区二区www| 在线天堂中文资源库| 久久久久久久久中文| 黄片小视频在线播放| 十分钟在线观看高清视频www| 麻豆久久精品国产亚洲av| 欧美成人一区二区免费高清观看 | 日韩精品青青久久久久久| 国产成人一区二区三区免费视频网站| 久久精品91蜜桃| 人成视频在线观看免费观看| 搡老熟女国产l中国老女人| 在线观看免费视频网站a站| 日韩欧美国产在线观看| avwww免费| 国产精品一区二区在线不卡| 国内毛片毛片毛片毛片毛片| 不卡av一区二区三区| 国产精品美女特级片免费视频播放器 | 欧美性长视频在线观看| 国产91精品成人一区二区三区| 成人18禁在线播放| 亚洲性夜色夜夜综合| 久久久久国产精品人妻aⅴ院| 久久久精品国产亚洲av高清涩受| 天天躁狠狠躁夜夜躁狠狠躁| www国产在线视频色| 19禁男女啪啪无遮挡网站| 午夜福利影视在线免费观看| 久9热在线精品视频| 国产精品久久久久久亚洲av鲁大| 国产在线精品亚洲第一网站| 国产精品久久久人人做人人爽| 一进一出抽搐动态| 一边摸一边抽搐一进一小说| 国产成人啪精品午夜网站| 午夜免费成人在线视频| 国产成人一区二区三区免费视频网站| 亚洲自拍偷在线| 99久久国产精品久久久|