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

    基于自適應(yīng)變分模態(tài)分解的擴展工況傳遞路徑分析方法

    2021-08-10 08:31:18張俊紅賈宏杰周啟迪朱小龍林杰威
    公路交通科技 2021年7期
    關(guān)鍵詞:頻響巴氏模態(tài)

    張俊紅,賈宏杰,周啟迪,朱小龍,林杰威,3

    (1. 天津大學 內(nèi)燃機燃燒學國家重點實驗室,天津 300072; 2. 天津仁愛學院,天津 301636;3. 濰柴動力股份有限公司,山東 濰坊 261016)

    0 引言

    汽車作為一個復雜的機械系統(tǒng),在運行過程由于多個激勵點經(jīng)不同的傳遞路徑合成導致汽車內(nèi)部產(chǎn)生聲振問題。在車輛研發(fā)過程中,為優(yōu)化整車NVH(Noise, Vibration, Harshness)性能,需要綜合考慮所有激勵和傳遞路徑。傳遞路經(jīng)分析[1]可有效解決上述問題。經(jīng)過20多年的發(fā)展,在TPA的基礎(chǔ)上衍生出了多種分析方法,如工況傳遞路徑分析(Operational Transfer Path Analysis,OPA)[2]和擴展工況傳遞路徑分析(Operational-X Transfer Path analysis,OPAX)等[3]。其中,OPAX[4]兼顧TPA的高精度和OPA的高效率,被廣泛應(yīng)用于工程和科研領(lǐng)域。

    OPAX進行參數(shù)化建模的過程,需要采集大量的工況數(shù)據(jù)[5-6]。因測試環(huán)境存在偶發(fā)性的沖擊和外界干擾噪聲,導致采集數(shù)據(jù)摻雜一定的噪聲信號,有效的弱特征振動信號可能會被掩蓋,使OPAX分析精度出現(xiàn)較大的偏差,無法準確獲得各條傳遞路徑的貢獻度。因此,為剔除干擾信號,保證有效信號的識別精度,需要對采集到的信號進行去噪處理。

    目前常用的去噪方法包含小波變換(WT)、經(jīng)驗?zāi)B(tài)分解(EMD)、變分模態(tài)分解(VMD)等[7-9]。WT基于積分理論,在分解信號的過程中過分依賴于小波基的選擇,不具有自適應(yīng)性[10]。EMD彌補了WT的不足,可自適應(yīng)的將復雜的信號分解為多個本征模態(tài)函數(shù)(Intrinsic Mode Function,IMF),但EMD在理論上缺乏嚴格的證明,分解信號過程中存在包絡(luò)誤差、模態(tài)混淆等問題[11]。文獻[12]在傳統(tǒng)維納濾波理論的基礎(chǔ)上提出了變分模態(tài)分解(Variational Mode Decomposition,VMD),該方法具有堅實的數(shù)學基礎(chǔ),在處理非線性非平穩(wěn)信號時具有很大的優(yōu)勢。然而,VMD方法中的模態(tài)個數(shù)K和懲罰因子α需要人為預設(shè),并且缺乏標準。

    針對VMD參數(shù)優(yōu)化問題,湯周杰等[13]通過融合魯棒獨立分量優(yōu)化VMD中的模態(tài)個數(shù)K,但沒有考慮懲罰因子α的影響。鄢小安等[14]利用瞬時能量分配-置換熵作為目標函數(shù)分別對VMD中的參數(shù)K和α進行優(yōu)化,然而忽略了K和α之間的耦合影響,可能會引起局部最優(yōu)結(jié)果。肖懷碩等[15]采用傅立葉頻譜包絡(luò)作為優(yōu)化目標函數(shù)對參數(shù)K和α同時優(yōu)化,但該目標函數(shù)對非穩(wěn)態(tài)信號適應(yīng)性較差。唐貴基等[16]以Shannon熵作為適應(yīng)度函數(shù),采用粒子群優(yōu)化算法對VMD進行優(yōu)化。然而以Shannon熵為適應(yīng)度函數(shù)的優(yōu)化結(jié)果魯棒性較差。多尺度模糊熵可以衡量時間序列在不同尺度因子下的復雜程度和維數(shù)變化時產(chǎn)生新信息的概率大小,考慮到采集信號中包含具有規(guī)律性的振動信號和隨機性很強的干擾信號,將多尺度模糊熵作為判斷信號特征的指標,用于振動信號的分解去噪。PSO具有較高的計算效率,但是全局搜索能力較差,容易陷入局部最優(yōu)解[17]。模擬退火算法[18]具有概率突跳能力,使粒子在停止進化的位置重新產(chǎn)生粒子進行迭代計算,可有效避免局部最優(yōu)。

    為了削減OPAX在信號采集過程干擾噪聲的影響,提高分析精度,本研究提出一種基于自適應(yīng)變分模態(tài)分解和巴氏距離的優(yōu)化OPAX方法。該方法考慮到多尺度模糊熵能夠較好地表征非穩(wěn)態(tài)復雜信號,將其作為適應(yīng)度函數(shù);采用模擬退火粒子群算法(Simulated Annealing Particle Swarm Optimization, SAPSO)進行信號自適應(yīng)變分模態(tài)分解,最后通過巴氏距離計算原始信號和分解信號概率密度函數(shù)的相似性進行相關(guān)模態(tài)篩選,實現(xiàn)信號去噪。

    1 AVMD-BD和OPAX方法原理

    1.1 AVMD-BD降噪算法

    1.1.1 AVMD算法

    文獻[12]已詳細介紹了VMD算法的原理,本研究不再贅述。VMD方法在處理信號時需要預先設(shè)定模態(tài)個數(shù)K和懲罰參數(shù)α,參數(shù)選擇不準確會導致信號產(chǎn)生過分解、欠分解和模態(tài)混疊的問題。針對VMD參數(shù)優(yōu)化的問題,由于PSO算法存在著全局搜索能力較差、易陷入局部最優(yōu)解、精度較差等問題,借鑒模擬退火算法中的概率突跳能力,使粒子跳出已搜索的最優(yōu)區(qū)域,在更大空間中開展搜索。

    基于此,本研究綜合考慮SAPSO算法全局搜索優(yōu)勢,以及多尺度模糊熵在提取復雜非穩(wěn)態(tài)信號的優(yōu)勢[19-20],提出一種以最小IMF多尺度模糊熵為適應(yīng)度值的SAPSO搜索算法進行VMD自適應(yīng)參數(shù)尋優(yōu)。

    粒子群算法的進化方程為:

    (1)

    (2)

    為提高粒子群算法的搜索效率,使得在搜索空間中隨機產(chǎn)生的粒子j以較大的概率處于最優(yōu)點附近,采用模擬退火算法生成。

    (3)

    式中,η為擾動幅值參數(shù);rand為隨機變量。

    AVMD算法中SAPSO的參數(shù)設(shè)置如表1所示。其中M為粒子的種群規(guī)模;N為迭代次數(shù);w為慣性權(quán)重;T為初始溫度;q為退火系數(shù);c1和c2為學習因子;V為粒子移動的速度范圍。

    表1 SAPSO參數(shù)設(shè)置Tab.1 SAPSO parameter setting

    AVMD算法自適應(yīng)獲取最優(yōu)參數(shù)流程見圖1,具體步驟如下:

    圖1 AVMD算法優(yōu)化流程Fig.1 Flowchart of optimization by AVMD algorithm

    (1)確定SAPSO中的所有初始化參數(shù)。

    (2)計算每個粒子的適應(yīng)度。

    (3)對比粒子的適應(yīng)度和最好位置pg,如較好,設(shè)置為當前最好位置。

    (5)根據(jù)(3)和(4)對其他粒子的速度和位置進行優(yōu)化。

    (6)滿足停止條件輸出最優(yōu)結(jié)果。

    1.1.2 相關(guān)模態(tài)的選擇

    概率密度函數(shù)(Probability Density Function,PDF)[21]可以反映不同信號分布之間的區(qū)別。本研究采用核密度估計的方法獲取每個分解信號和原始信號的PDF。計算信號之間的相似性作為區(qū)分相關(guān)模態(tài)的依據(jù)。巴氏距離[22]是2個統(tǒng)計樣本之間的重疊量的距離測量,用于計算2個樣本之間的相對接近程度。巴氏距離的定義如下:

    在定義域X中,概率分布Q和P的巴氏距離的定義為:

    DB(P,Q)=-ln[BC(P,Q)]。

    (4)

    對于離散概率分布:

    (5)

    對于連續(xù)概率分布:

    (6)

    可以看出,DB越小表示概率分布越接近,即模態(tài)分量與輸入信號更相關(guān)。

    輸入信號和模態(tài)分量之間相似性的定義為:

    L(i)=distance[pdf(x(t)),pdf(IMFi(t))]。

    (7)

    通過評估2個相鄰IMF之間距離的斜率來確定相關(guān)模式。當斜率顯著增加時,表明該IMF之后發(fā)生的相似性迅速降低。定義θ為2個相鄰的IMF與輸入信號之間的最大斜率,表達式為:

    θ=max|L(i+1)-L(i)|,i=1, 2,…,n-1。

    (8)

    如在IMFm和IMFm+1之間產(chǎn)生最大斜率,則通過式(9)獲得濾波后的信號y(t)。

    (9)

    1.2 OPAX法基本原理

    傳統(tǒng)的TPA通常使用逆矩陣法求解工作載荷,其理論公式為:

    (10)

    式中,yk(ω)為目標點k的總貢獻量;ω為頻率;n為傳遞路徑數(shù)量;Hik(ω)和Hik(ω)分別為不耦合的主動端與響應(yīng)端之間的結(jié)構(gòu)和聲學頻率響應(yīng)函數(shù)矩陣;Fik(ω)和Qjk(ω)分別為作用在主動端的結(jié)構(gòu)載荷和聲學載荷。

    OPAX方法是在傳統(tǒng)TPA方法基礎(chǔ)上演化而來的,與傳統(tǒng)TPA方法的主要區(qū)別在于使用了參數(shù)化模型進行載荷識別,同時借鑒了OPA方法的優(yōu)點,大量采用實測工況數(shù)據(jù)進行建模,減少了FRF的測試工作量,在保證結(jié)果精度的同時節(jié)省了建模時間。

    在整車系統(tǒng)中,動力總成作為激勵源為主動側(cè),車身結(jié)構(gòu)為被動側(cè),中間的耦合元件為懸置,起到支撐及隔離振動的作用,目標點是某點在某自由度上的振動。對于目標點的振動響應(yīng),被動側(cè)在耦合點的每個自由度到響應(yīng)點均會形成1條獨立的傳遞路徑。考慮縱向、側(cè)向和垂向3個方向自由度。任一方向的頻響函數(shù)可以表示為:

    (11)

    式中,Vi(ω)為結(jié)構(gòu)表面的振動加速度;Q(ω)為麥克風處的體積速度。

    已知某傳遞路徑上的系統(tǒng)頻響函數(shù)和工作載荷,該路徑對目標位置噪聲的貢獻yki(ω)可以表示為:

    yki(ω)=Hki(ω)Fi(ω),

    (12)

    式中,k為響應(yīng)點的表示符號;Hki(ω)為第i條傳遞路徑上的頻響函數(shù);Fi(ω)為第i條路徑上的工作載荷。

    Vi(ω)直接影響頻響函數(shù)的計算精度,進而影響傳遞路徑分析的計算結(jié)果。

    如果有n條傳遞路徑,則目標點的總響應(yīng)可認為是多條路徑分貢獻量的線性疊加,可以表示為:

    (13)

    式中yk(ω)為目標點k的響應(yīng)。

    2 仿真試驗與分析

    考慮到振動信號中包含多種成分,同時測試環(huán)境中存在偶發(fā)性的沖擊和外界干擾噪聲。本研究仿真試驗由2部分組成:(1)驗證AVMD算法的準確性。由于高斯白噪聲中包含多種不確定分量,無法驗證AVMD分解結(jié)果的準確性,因此模擬信號x(t)由包含明確特征的多分量和非高斯信號組成。(2)驗證AVMD-BD去噪效果的有效性。在信號x(t)中增加高斯白噪聲模擬外界干擾噪聲的影響。

    2.1 AVMD算法驗證

    為驗證AVMD算法的準確性,模擬信號x(t)由4類信號疊加而成。

    針對下式的仿真信號進行分解:

    x(t)=x1(t)+x2(t)+x3(t)+x4(t),

    (14)

    采用PSO-VMD和SAPSO-VMD這2種方法對上述信號進行分解,計算得到的最優(yōu)參數(shù)組合為[K=5,α=436]和[K=4,α=674],從優(yōu)化參數(shù)組合中可以看出,PSO-VMD的K=5,而原始信號中只有4組信號,表明該方法出現(xiàn)過分解現(xiàn)象。SAPSO-VMD可以準確分解4個模擬信號。上述分析結(jié)果表明,SAPSO-VMD可以有效分離混合信號,為后續(xù)通過巴氏距離挑選相關(guān)模態(tài)奠定良好的基礎(chǔ)。

    2.2 AVMD-BD算法驗證

    用于計算頻響函數(shù)的信號中若摻雜干擾噪聲會嚴重影響OPAX的計算精度。為驗證AVMD-BD算法在信號去噪效果中的有效性,本研究在式(14)信號的基礎(chǔ)上疊加高斯白噪聲來進行模擬仿真。模擬信號時域結(jié)果如圖2所示。

    圖2 模擬信號仿真Fig.2 Analog signal simulation

    用AVMD對上述模擬信號進行分解,計算得到最佳的參數(shù)優(yōu)化組合為[K=16,α=3 052]。

    為了有效地剔除干擾參量,實現(xiàn)信號的精確去噪,在采用巴氏距離進行特征模態(tài)選擇時,首先利用核密度估計的方法獲得輸入信號與每個模態(tài)分量的pdf。其次根據(jù)式(4),得到信號之間的巴氏距離,如圖3所示。可以看出,巴氏距離將分解信號分為2組,每組中的xi(t)與輸入信號之間的距離相對集中,從x3(t)到x4(t)的過程中,斜率顯著增加,表明x4(t)后發(fā)生的相似性迅速減少。另外,從x4(t)到x16(t)與原始信號之間的距離相對較遠,均可以認為是遠離原始信號的非相關(guān)模態(tài)。x1(t)到x3(t)即為相關(guān)模態(tài),通過重構(gòu)相關(guān)模態(tài),從而得到去噪后的模擬信號。

    圖3 各分解信號與模擬含噪信號之間的巴氏距離Fig.3 BD between each decomposed signal and analog noisy signal

    3 AVMD-BD-OPAX效果驗證

    3.1 OPAX模型建立及工況數(shù)據(jù)采集

    本研究選取的試驗對象為某國產(chǎn)輕型客車,裝載一款4缸柴油發(fā)動機,發(fā)動機排量為2.0 L,額定乘員為21人。動力總成懸置共有3個,傳動軸懸置有2個,排氣懸置有4個,后輪鋼板彈簧懸置左右共有6個。整車動態(tài)試驗在90 km/h勻速工況進行測試,利用數(shù)采系統(tǒng)(LMS,SCADAS 05)和LMS TEST.LAB軟件對數(shù)據(jù)進行采集與分析。試驗中在懸置主被動端分別安裝三向加速度傳感器(PCB,356A26)進行振動信號采集,在車輪、發(fā)動機進排氣端和駕駛員右耳布置傳聲器(PCB,378B02)進行噪聲信號采集,傳遞路徑分析模型如圖4所示,包括15×3+7=52條傳遞路徑。傳感器安裝位置如圖5所示。

    圖4 傳遞路徑分析模型示意圖Fig.4 Schematic diagram of transfer path analysis model

    圖5 測點布置示意圖Fig.5 Schematic diagram of measuring point arrangement

    3.2 頻響函數(shù)的測量

    采用聲學互易法進行頻響函數(shù)的測量。將體積聲源作為輸入信號放置在車內(nèi)駕駛員位置,測取整車37個測點處的信號。對去噪后的信號利用式(11)計算相應(yīng)傳遞路徑的頻響函數(shù)。

    3.3 AVMD-BD-OPAX結(jié)果

    圖6為車輛90 km/h勻速行駛工況下駕駛員右耳噪聲的平均功率譜??梢钥闯?,噪聲能量主要集中在0~300 Hz的中低頻區(qū)域。因此主要關(guān)注該區(qū)域的優(yōu)化結(jié)果。

    圖6 駕駛員右耳噪聲平均功率譜Fig.6 Average power spectrum of driver’s right ear noise

    將目標點處的實測信號作為參考值,對比驗證本研究所提方法的準確性和優(yōu)越性。圖7為原始OPAX計算值和AVMD-BD優(yōu)化后OPAX計算值??梢悦黠@看出,優(yōu)化前后所得結(jié)果的走勢與實測值基本一致,在采用AVMD-BD降噪后所得結(jié)果與實測值的吻合程度更高,在40~280 Hz頻率段內(nèi)尤為明顯。

    圖7 原OPAX方法與改進后方法的計算結(jié)果Fig.7 Calculation results of original OPAX method and improved method

    為進一步驗證AVMD-BD方法在振動信號去噪處理中的有效性,以信噪比SNR和均方根誤差RMSE來作為評判的依據(jù)。其中均方根誤差越小,信噪比越大,表征信號的去噪效果越好。

    (15)

    (16)

    式中,f′(k)為去噪信號;f(k)為原始信號;k為信號長度。

    表2列出了AVMD-BD降噪優(yōu)化后相對于實測信號的SNR和RMSE??梢悦黠@看出,在采用AVMD-BD降噪優(yōu)化后,信噪比提升了71.2%,均方根誤差減小了66.9%。說明本研究所提的降噪優(yōu)化方案可有效提高OPAX方法的計算精度。

    表2 AVMD-BD降噪后SNR和RMSE比較Tab.2 Comparison of SNR and RMSE after AVMD-BD denoising

    在40~280 Hz內(nèi)振動響應(yīng)較大,且有8個比較明顯的峰值,這也是工程研究人員最為關(guān)心的問題。

    表3中列出了在8處峰值頻率處,計算值和實測值之間的相對誤差。結(jié)果顯示經(jīng)AVMD-BD優(yōu)化后的OPAX誤差均控制在5%以內(nèi),表明優(yōu)化后OPAX計算結(jié)果的精確性得到有效提高。

    表3 2種方法的峰值相對誤差比較Tab.3 Comparison of peak relative errors between 2 methods

    4 結(jié)論

    本研究針對采用擴展傳遞路徑分析(OPAX)進行頻響函數(shù)計算時,由于存在干擾信號導致頻響函數(shù)計算準確性低進而影響OPAX計算精度的問題,提出了一種基于自適應(yīng)變分模態(tài)分解和巴氏距離的優(yōu)化OPAX方法。該方法提高了傳遞路徑分析的精確性。具體結(jié)論如下:

    (1)提出了一種基于模擬退火粒子群算法的自適應(yīng)變分模態(tài)分解方法,該方法考慮到多尺度模糊熵能夠較好地表征非穩(wěn)態(tài)復雜信號,將其作為適應(yīng)度函數(shù);利用模擬退火算法的概率突跳能力,可有效避免局部最優(yōu)。通過仿真信號分析,將其與PSO-VMD進行對比,結(jié)果表明,AVMD可以有效分離混合信號中的各個模態(tài),其分解結(jié)果優(yōu)于PSO-VMD。

    (2)采用巴氏距離可以有效識別信號分量與有效信號的相關(guān)性,實現(xiàn)高信噪比的信號重構(gòu)。重構(gòu)信號的走勢與原信號基本一致,且波形更加光滑,降噪效果明顯。

    (3)提出了一種基于AVMD-BD的OPAX優(yōu)化方法。試驗分析結(jié)果表明,優(yōu)化后的計算值更加接近實測值,優(yōu)化后OPAX計算值的信噪比提高了71.2%,均方根誤差減小了66.9%,在峰值頻率處的誤差均控制在5%以內(nèi)。

    猜你喜歡
    頻響巴氏模態(tài)
    釋放巴氏新小綏螨可滿足對蘋果全爪螨的防治需求
    基于分塊化頻響函數(shù)曲率比的砌體房屋模型損傷識別研究
    美團外賣哥
    巴氏殺菌水牛奶在不同儲存條件下微生物增長規(guī)律的研究
    頻響函數(shù)殘差法在有限元模型修正中的應(yīng)用
    頻響阻抗法診斷變壓器繞組變形
    巴氏醋桿菌核酸修復酶UvrA對大腸桿菌耐受性的影響
    國內(nèi)多模態(tài)教學研究回顧與展望
    基于HHT和Prony算法的電力系統(tǒng)低頻振蕩模態(tài)識別
    由單個模態(tài)構(gòu)造對稱簡支梁的抗彎剛度
    計算物理(2014年2期)2014-03-11 17:01:39
    91在线观看av| 色哟哟哟哟哟哟| 男人舔女人的私密视频| 国产黄a三级三级三级人| 禁无遮挡网站| 欧洲精品卡2卡3卡4卡5卡区| 国产免费男女视频| 成人三级黄色视频| 久久中文字幕一级| 午夜老司机福利片| 91在线观看av| 美女大奶头视频| 久久久国产精品麻豆| 男女视频在线观看网站免费 | 又大又爽又粗| 亚洲男人的天堂狠狠| 欧美乱色亚洲激情| 午夜福利高清视频| 波多野结衣高清作品| 久久精品91无色码中文字幕| 日本三级黄在线观看| 日本熟妇午夜| www日本黄色视频网| 欧美性猛交╳xxx乱大交人| 可以在线观看的亚洲视频| 黑人欧美特级aaaaaa片| 999精品在线视频| 欧美黄色片欧美黄色片| 久久性视频一级片| 中文字幕久久专区| 亚洲国产欧美日韩在线播放| 视频区欧美日本亚洲| 欧美人与性动交α欧美精品济南到| 久久久久精品国产欧美久久久| www日本在线高清视频| 免费看a级黄色片| 久久久国产精品麻豆| 黄色片一级片一级黄色片| 精品第一国产精品| 每晚都被弄得嗷嗷叫到高潮| 视频区欧美日本亚洲| 女警被强在线播放| 91九色精品人成在线观看| www.自偷自拍.com| 久久人妻福利社区极品人妻图片| 精品一区二区三区四区五区乱码| 91九色精品人成在线观看| 精品人妻1区二区| 90打野战视频偷拍视频| 欧美一级毛片孕妇| 两性午夜刺激爽爽歪歪视频在线观看 | 国产激情偷乱视频一区二区| 此物有八面人人有两片| 叶爱在线成人免费视频播放| 精品日产1卡2卡| 中文字幕精品免费在线观看视频| 精品久久久久久,| 午夜亚洲福利在线播放| 成人亚洲精品av一区二区| 免费在线观看视频国产中文字幕亚洲| 亚洲精品久久成人aⅴ小说| 久久久久久久精品吃奶| 日本三级黄在线观看| 亚洲专区国产一区二区| 亚洲国产欧洲综合997久久, | 午夜免费成人在线视频| 精品免费久久久久久久清纯| 熟女少妇亚洲综合色aaa.| 亚洲va日本ⅴa欧美va伊人久久| 村上凉子中文字幕在线| 99国产综合亚洲精品| 欧美色欧美亚洲另类二区| 麻豆一二三区av精品| 欧美一区二区精品小视频在线| 长腿黑丝高跟| 精品熟女少妇八av免费久了| 91九色精品人成在线观看| 村上凉子中文字幕在线| a级毛片在线看网站| 国产精品久久电影中文字幕| 成人特级黄色片久久久久久久| 欧美激情 高清一区二区三区| 99riav亚洲国产免费| 在线观看午夜福利视频| 九色国产91popny在线| 亚洲电影在线观看av| 在线免费观看的www视频| 在线观看午夜福利视频| 亚洲av中文字字幕乱码综合 | 91麻豆av在线| 女同久久另类99精品国产91| 亚洲精品粉嫩美女一区| 在线看三级毛片| 女人爽到高潮嗷嗷叫在线视频| 男女午夜视频在线观看| av视频在线观看入口| 亚洲av电影不卡..在线观看| 国产精品98久久久久久宅男小说| 日本精品一区二区三区蜜桃| 午夜老司机福利片| 手机成人av网站| 日韩欧美一区视频在线观看| 青草久久国产| 色播在线永久视频| 18禁黄网站禁片免费观看直播| 满18在线观看网站| 九色国产91popny在线| 日韩欧美 国产精品| 亚洲久久久国产精品| 可以在线观看的亚洲视频| 91在线观看av| 日韩欧美一区二区三区在线观看| 一级作爱视频免费观看| 不卡av一区二区三区| 国产成人系列免费观看| 久久久久九九精品影院| 窝窝影院91人妻| 嫩草影视91久久| 丝袜美腿诱惑在线| 国内久久婷婷六月综合欲色啪| 欧美日韩精品网址| 91麻豆精品激情在线观看国产| 在线观看免费日韩欧美大片| 色播在线永久视频| 国产蜜桃级精品一区二区三区| 男女视频在线观看网站免费 | 久久久久久亚洲精品国产蜜桃av| 久久人妻av系列| 日本成人三级电影网站| 啦啦啦免费观看视频1| 啪啪无遮挡十八禁网站| 国产精品一区二区免费欧美| 国产av在哪里看| 满18在线观看网站| 亚洲在线自拍视频| 亚洲美女黄片视频| 黄片大片在线免费观看| 国产亚洲精品av在线| 亚洲第一av免费看| 亚洲色图 男人天堂 中文字幕| 国产激情偷乱视频一区二区| 中文字幕人成人乱码亚洲影| 欧美成狂野欧美在线观看| 日日爽夜夜爽网站| 搡老岳熟女国产| a级毛片在线看网站| 一级毛片女人18水好多| 久久人人精品亚洲av| 亚洲aⅴ乱码一区二区在线播放 | 男人舔女人的私密视频| 色精品久久人妻99蜜桃| 精品高清国产在线一区| 在线观看一区二区三区| 一级作爱视频免费观看| 国产乱人伦免费视频| 岛国视频午夜一区免费看| 免费在线观看黄色视频的| 精品高清国产在线一区| 一进一出好大好爽视频| 在线看三级毛片| 精品福利观看| 久久久久久久午夜电影| 色av中文字幕| 老熟妇仑乱视频hdxx| 国产精品精品国产色婷婷| 人妻久久中文字幕网| 亚洲电影在线观看av| 午夜福利免费观看在线| 成人18禁在线播放| 999久久久国产精品视频| 欧美日韩乱码在线| 欧美一级a爱片免费观看看 | 一个人观看的视频www高清免费观看 | 叶爱在线成人免费视频播放| 国产又黄又爽又无遮挡在线| 亚洲自拍偷在线| 欧美zozozo另类| 国产在线精品亚洲第一网站| 国产又色又爽无遮挡免费看| 麻豆国产av国片精品| 亚洲av片天天在线观看| 宅男免费午夜| 非洲黑人性xxxx精品又粗又长| 日本熟妇午夜| 女性生殖器流出的白浆| 国产高清videossex| 精品一区二区三区视频在线观看免费| 欧美性猛交╳xxx乱大交人| www.熟女人妻精品国产| 亚洲精品一区av在线观看| 18美女黄网站色大片免费观看| 亚洲 欧美一区二区三区| 欧美av亚洲av综合av国产av| 桃红色精品国产亚洲av| 亚洲va日本ⅴa欧美va伊人久久| 免费在线观看完整版高清| 无人区码免费观看不卡| 国产亚洲av高清不卡| 黄色视频,在线免费观看| 成人18禁在线播放| 成人三级黄色视频| 午夜老司机福利片| 亚洲精品在线美女| 搡老妇女老女人老熟妇| 中文字幕av电影在线播放| av在线播放免费不卡| 99精品久久久久人妻精品| 超碰成人久久| 成人免费观看视频高清| 国产成人av激情在线播放| 久久人人精品亚洲av| 亚洲精品久久国产高清桃花| 国产免费男女视频| 淫秽高清视频在线观看| 最新美女视频免费是黄的| 99国产极品粉嫩在线观看| av超薄肉色丝袜交足视频| 一级毛片精品| 草草在线视频免费看| 亚洲人成电影免费在线| 一边摸一边做爽爽视频免费| 亚洲av第一区精品v没综合| 精品久久久久久久久久免费视频| 国产又色又爽无遮挡免费看| 最新在线观看一区二区三区| 一区二区三区高清视频在线| 亚洲国产欧美一区二区综合| 国产精品亚洲美女久久久| 国产精品,欧美在线| av有码第一页| 熟女少妇亚洲综合色aaa.| 91字幕亚洲| 亚洲国产精品成人综合色| 中文字幕久久专区| 18禁国产床啪视频网站| 成人国语在线视频| 日韩精品青青久久久久久| 1024视频免费在线观看| 两性午夜刺激爽爽歪歪视频在线观看 | 天天添夜夜摸| 欧美久久黑人一区二区| 国内揄拍国产精品人妻在线 | 亚洲精品在线美女| 中国美女看黄片| 成人永久免费在线观看视频| 美女高潮到喷水免费观看| 精品久久久久久久久久久久久 | 男人的好看免费观看在线视频 | 亚洲国产精品sss在线观看| 国产欧美日韩一区二区精品| 久久国产精品影院| 欧美成人性av电影在线观看| 一级毛片女人18水好多| 欧美黄色片欧美黄色片| 国产人伦9x9x在线观看| 久热这里只有精品99| 成年免费大片在线观看| 免费在线观看影片大全网站| 日韩高清综合在线| 成年版毛片免费区| 狠狠狠狠99中文字幕| 亚洲精品av麻豆狂野| 韩国精品一区二区三区| 久久香蕉国产精品| 欧美一级毛片孕妇| 村上凉子中文字幕在线| 免费在线观看亚洲国产| 国产爱豆传媒在线观看 | 人人澡人人妻人| 丁香欧美五月| 日本一区二区免费在线视频| 亚洲中文字幕一区二区三区有码在线看 | 在线视频色国产色| 99热6这里只有精品| 国产亚洲精品久久久久久毛片| 欧美日韩乱码在线| 美国免费a级毛片| 亚洲精品美女久久久久99蜜臀| 校园春色视频在线观看| 天天躁狠狠躁夜夜躁狠狠躁| av福利片在线| 一区二区三区高清视频在线| 最近在线观看免费完整版| 日韩精品中文字幕看吧| 日韩欧美免费精品| 亚洲精品av麻豆狂野| 黄片小视频在线播放| 精品日产1卡2卡| 国产男靠女视频免费网站| 久久午夜亚洲精品久久| 久久久久久免费高清国产稀缺| 欧美zozozo另类| 久久久久免费精品人妻一区二区 | 伦理电影免费视频| 两个人免费观看高清视频| 久久这里只有精品19| 在线观看www视频免费| 免费人成视频x8x8入口观看| 欧美精品啪啪一区二区三区| 日本a在线网址| 9191精品国产免费久久| 免费高清在线观看日韩| 亚洲 欧美一区二区三区| 国产一卡二卡三卡精品| 欧美大码av| 午夜福利欧美成人| 99久久综合精品五月天人人| 婷婷精品国产亚洲av在线| 精品欧美国产一区二区三| 国产精品免费视频内射| 麻豆av在线久日| 美女国产高潮福利片在线看| 51午夜福利影视在线观看| 国产1区2区3区精品| 男女床上黄色一级片免费看| 国产成人av激情在线播放| 免费av毛片视频| 成熟少妇高潮喷水视频| 男女下面进入的视频免费午夜 | 日韩欧美国产一区二区入口| 国产极品粉嫩免费观看在线| www.自偷自拍.com| 精品一区二区三区四区五区乱码| 午夜福利在线在线| xxxwww97欧美| 国产av在哪里看| 国产麻豆成人av免费视频| 亚洲国产欧美日韩在线播放| av福利片在线| 久久久久久九九精品二区国产 | 中文字幕久久专区| 99在线人妻在线中文字幕| 人人澡人人妻人| 国产精品,欧美在线| 高清在线国产一区| 午夜福利一区二区在线看| 91麻豆av在线| av在线播放免费不卡| 精品欧美一区二区三区在线| 一级a爱片免费观看的视频| 在线看三级毛片| 亚洲av成人不卡在线观看播放网| 无人区码免费观看不卡| 日本五十路高清| 18禁黄网站禁片午夜丰满| 欧美亚洲日本最大视频资源| 国产av又大| 精品一区二区三区视频在线观看免费| 大香蕉久久成人网| 在线观看免费视频日本深夜| 国产成人影院久久av| 婷婷精品国产亚洲av| 一个人观看的视频www高清免费观看 | 后天国语完整版免费观看| 亚洲精品在线美女| 亚洲一卡2卡3卡4卡5卡精品中文| 亚洲全国av大片| 欧美人与性动交α欧美精品济南到| 国产一级毛片七仙女欲春2 | 国产真人三级小视频在线观看| 免费看日本二区| 国产亚洲精品第一综合不卡| 久久久久国产一级毛片高清牌| 午夜福利18| xxxwww97欧美| 搡老岳熟女国产| 精品欧美一区二区三区在线| 欧美黑人精品巨大| www.熟女人妻精品国产| 黄片小视频在线播放| 九色国产91popny在线| 99热只有精品国产| 国产成人av激情在线播放| 国产精品二区激情视频| 免费观看精品视频网站| 成人亚洲精品一区在线观看| 91在线观看av| 亚洲七黄色美女视频| 熟妇人妻久久中文字幕3abv| 欧美成人性av电影在线观看| 免费无遮挡裸体视频| 欧美zozozo另类| 成年版毛片免费区| a在线观看视频网站| 国产亚洲av高清不卡| 午夜免费激情av| 亚洲成人免费电影在线观看| 国产熟女午夜一区二区三区| 日本免费一区二区三区高清不卡| 美女 人体艺术 gogo| 韩国精品一区二区三区| 欧美中文日本在线观看视频| 久久亚洲精品不卡| 精品久久久久久久久久久久久 | 欧美日韩亚洲国产一区二区在线观看| 少妇被粗大的猛进出69影院| 成人18禁高潮啪啪吃奶动态图| 国产私拍福利视频在线观看| 黄片播放在线免费| 欧美日韩中文字幕国产精品一区二区三区| 亚洲av成人一区二区三| 久久中文字幕一级| 手机成人av网站| 日韩大尺度精品在线看网址| 在线天堂中文资源库| 麻豆一二三区av精品| 日韩视频一区二区在线观看| 熟妇人妻久久中文字幕3abv| 欧美黑人巨大hd| 亚洲欧美精品综合一区二区三区| 午夜免费激情av| 久久婷婷成人综合色麻豆| 免费在线观看黄色视频的| 国产精品影院久久| 国产成人欧美在线观看| 久久婷婷人人爽人人干人人爱| 丰满人妻熟妇乱又伦精品不卡| 黄片小视频在线播放| 高潮久久久久久久久久久不卡| 免费高清视频大片| 黄片播放在线免费| 亚洲精品粉嫩美女一区| 亚洲一区二区三区不卡视频| 国产成人精品无人区| 亚洲专区国产一区二区| 最近在线观看免费完整版| 天天躁夜夜躁狠狠躁躁| 桃色一区二区三区在线观看| 国产一区在线观看成人免费| 国产av一区在线观看免费| 亚洲第一青青草原| 真人做人爱边吃奶动态| 久久人妻福利社区极品人妻图片| 777久久人妻少妇嫩草av网站| 久久久久九九精品影院| 99久久久亚洲精品蜜臀av| bbb黄色大片| av天堂在线播放| 99久久综合精品五月天人人| 国产私拍福利视频在线观看| 国产激情偷乱视频一区二区| 一夜夜www| 亚洲av五月六月丁香网| 午夜免费激情av| 一进一出抽搐gif免费好疼| 一进一出抽搐动态| 精华霜和精华液先用哪个| 欧洲精品卡2卡3卡4卡5卡区| 午夜激情av网站| 午夜福利在线观看吧| 国产视频一区二区在线看| 最近在线观看免费完整版| 人人妻,人人澡人人爽秒播| 欧美一级毛片孕妇| 黄片小视频在线播放| 免费一级毛片在线播放高清视频| 午夜亚洲福利在线播放| 午夜激情av网站| 香蕉国产在线看| 久久欧美精品欧美久久欧美| 熟女少妇亚洲综合色aaa.| 黄片小视频在线播放| 国产激情偷乱视频一区二区| 日日干狠狠操夜夜爽| 黄色成人免费大全| 欧美日韩一级在线毛片| 亚洲自拍偷在线| 欧美国产日韩亚洲一区| 国产午夜精品久久久久久| 欧美日本视频| 午夜福利在线观看吧| 12—13女人毛片做爰片一| 国产97色在线日韩免费| 两个人免费观看高清视频| 啦啦啦观看免费观看视频高清| 成在线人永久免费视频| 一二三四在线观看免费中文在| 久久九九热精品免费| 日韩欧美国产一区二区入口| 欧美中文综合在线视频| 国产精品一区二区三区四区久久 | 亚洲第一青青草原| 亚洲av成人一区二区三| a级毛片在线看网站| 日韩欧美 国产精品| 欧美中文日本在线观看视频| 亚洲,欧美精品.| 久久 成人 亚洲| 少妇裸体淫交视频免费看高清 | 男人舔女人的私密视频| 精品日产1卡2卡| 99精品欧美一区二区三区四区| 男男h啪啪无遮挡| 欧美zozozo另类| 国产99白浆流出| 成年免费大片在线观看| 国产一卡二卡三卡精品| 国产亚洲精品av在线| 最近最新免费中文字幕在线| 91字幕亚洲| 欧美另类亚洲清纯唯美| 一二三四社区在线视频社区8| 一本综合久久免费| 亚洲精品久久国产高清桃花| 久久久久亚洲av毛片大全| 成熟少妇高潮喷水视频| 亚洲精品美女久久av网站| 欧美 亚洲 国产 日韩一| 一区二区三区精品91| 国产男靠女视频免费网站| 久久久久精品国产欧美久久久| 欧美日韩中文字幕国产精品一区二区三区| 欧美不卡视频在线免费观看 | 热re99久久国产66热| 欧美成人性av电影在线观看| 国产成人欧美| 久久精品国产99精品国产亚洲性色| 精品国产美女av久久久久小说| 男男h啪啪无遮挡| 久久亚洲真实| 久久久久久人人人人人| 婷婷丁香在线五月| 亚洲精品美女久久av网站| 亚洲最大成人中文| 欧美乱妇无乱码| 天天躁狠狠躁夜夜躁狠狠躁| 亚洲在线自拍视频| 亚洲精品中文字幕一二三四区| 成人特级黄色片久久久久久久| 国产亚洲欧美98| 国产欧美日韩一区二区三| 国产1区2区3区精品| 男女视频在线观看网站免费 | 日韩视频一区二区在线观看| 精品少妇一区二区三区视频日本电影| 亚洲欧美激情综合另类| 国产私拍福利视频在线观看| 久久精品亚洲精品国产色婷小说| 亚洲成av人片免费观看| 久久国产乱子伦精品免费另类| 亚洲人成网站高清观看| 中文字幕另类日韩欧美亚洲嫩草| 夜夜看夜夜爽夜夜摸| av片东京热男人的天堂| 国产一卡二卡三卡精品| 成人永久免费在线观看视频| 制服诱惑二区| 国产不卡一卡二| 欧美最黄视频在线播放免费| 啪啪无遮挡十八禁网站| 日韩欧美国产一区二区入口| 亚洲欧洲精品一区二区精品久久久| 欧美性猛交黑人性爽| 亚洲五月色婷婷综合| e午夜精品久久久久久久| 久久精品国产清高在天天线| 国产亚洲精品av在线| 亚洲五月婷婷丁香| 亚洲美女黄片视频| 亚洲欧美日韩高清在线视频| 亚洲aⅴ乱码一区二区在线播放 | 黄色视频,在线免费观看| 亚洲三区欧美一区| 老司机午夜福利在线观看视频| АⅤ资源中文在线天堂| 岛国在线观看网站| 国产麻豆成人av免费视频| 色综合亚洲欧美另类图片| 女人爽到高潮嗷嗷叫在线视频| 亚洲成av人片免费观看| 男人舔女人下体高潮全视频| 看黄色毛片网站| 99久久无色码亚洲精品果冻| 91字幕亚洲| 中文在线观看免费www的网站 | 午夜福利成人在线免费观看| 亚洲av日韩精品久久久久久密| 女同久久另类99精品国产91| 国产成人欧美在线观看| 亚洲精品粉嫩美女一区| 少妇 在线观看| xxxwww97欧美| 一二三四社区在线视频社区8| 亚洲午夜精品一区,二区,三区| 中文字幕精品亚洲无线码一区 | √禁漫天堂资源中文www| 嫩草影院精品99| 欧美日韩黄片免| 国产精品久久电影中文字幕| 天天添夜夜摸| 老汉色∧v一级毛片| 国产久久久一区二区三区| 一个人免费在线观看的高清视频| 亚洲精品在线美女| 99热这里只有精品一区 | 美女免费视频网站| 成人一区二区视频在线观看| 母亲3免费完整高清在线观看| 亚洲精品一区av在线观看| 免费观看精品视频网站| 午夜福利在线观看吧| 亚洲国产毛片av蜜桃av| 啦啦啦韩国在线观看视频| 欧美成人一区二区免费高清观看 | 国产精品影院久久| 国产成人av激情在线播放| 男女视频在线观看网站免费 | 亚洲精品中文字幕一二三四区| 伊人久久大香线蕉亚洲五| 国产精品一区二区免费欧美| 亚洲国产精品成人综合色| 男女午夜视频在线观看| 成人一区二区视频在线观看| 日韩大码丰满熟妇| 97人妻精品一区二区三区麻豆 |