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

    多尺度特征組合優(yōu)化的航空液壓管路故障診斷

    2022-10-13 08:41:30薛政坤汪曦于曉光王寵劉思遠張景博
    機床與液壓 2022年18期
    關(guān)鍵詞:學習機管路尺度

    薛政坤,汪曦,于曉光,王寵,劉思遠,張景博

    (遼寧科技大學機械工程與自動化學院,遼寧鞍山 114051)

    0 前言

    航空液壓系統(tǒng)是航空機械設(shè)備重要的動力支持系統(tǒng),其下的液壓管路系統(tǒng)是航空液壓系統(tǒng)中不可或缺的一部分,擔負著能量傳輸?shù)娜蝿?wù)。液壓管路長期服役于復雜的環(huán)境中,受高壓環(huán)境下強烈振動的影響,導致管體自身疲勞以及過載和碰撞等問題出現(xiàn),極易造成積累性的損傷故障,影響整個液壓系統(tǒng)正常運作,甚至可能引發(fā)嚴重的安全事故。航空液壓管路故障一直是影響航空液壓系統(tǒng)可靠性的重要問題,因此開展航空液壓管路故障診斷研究,對提高航空液壓系統(tǒng)的可靠性具有重要的意義。

    航空液壓管路的振動機制具有一定復雜性,一方面泵源壓力脈動頻率與液壓流體諧振頻率接近時易引發(fā)耦合振動,另一方面泵源壓力脈動頻率與管路結(jié)構(gòu)固有頻率接近時易引發(fā)流固耦合振動,導致管路振動信號呈現(xiàn)非平穩(wěn)和非線性等特點。從復雜的振動信號中提取出故障特征,是實現(xiàn)航空液壓管路故障診斷的關(guān)鍵。李哲洙等將經(jīng)驗模態(tài)分解方法結(jié)合希伯爾特黃變換方法引入到航空液壓管路的振動信號分析,實現(xiàn)管路裂紋故障的診斷。王立文等針對航空液壓管路破裂泄漏等故障問題,利用仿真建模得到不同液壓泄漏故障程度的管路的固有頻率作為特征信息組成數(shù)據(jù)樣本,利用概率神經(jīng)網(wǎng)絡(luò)模型實現(xiàn)故障識別。

    為進一步提取信號中的特征信息,可以利用信號的信息熵作為特征指標來表示,它用于反映時域序列的隨機性和復雜程度。ROSTAGHI和AZAMI提出一種用于衡量時間序列復雜度的指標散布熵(Dispersion Entropy,DE),與樣本熵和排列熵相比具有更好的穩(wěn)定性和較強的計算速度。

    由于非線性非平穩(wěn)振動信號的復雜性,僅從單一尺度衡量信號的復雜度,可能會遺漏其他尺度上的重要信息,難以全面反映特征信息。因此,多尺度信息熵作為判斷不同尺度信號復雜程度的指標和特征參數(shù),被用來表征不同類型的振動信號在不同尺度下的復雜性。例如多尺度樣本熵、多尺度排列熵等多尺度信息熵作為表征信號信息的特征指標,并被引入故障診斷領(lǐng)域中。

    為從振動信號中提取更具狀態(tài)表征性的特征信息,本文作者采用能量比值法對變分模態(tài)分解的分解層數(shù)進行優(yōu)化選取。根據(jù)模態(tài)分量信號與原信號的相關(guān)系數(shù),選取最佳的模態(tài)分量信號進行模態(tài)信號重構(gòu)作為分析信號。計算重構(gòu)模態(tài)信號的優(yōu)化多尺度散布熵(Multiscale Dispersion Entropy,MDE)作為多尺度特征指標,將其熵值組合成具有表征性的特征向量集。利用經(jīng)麻雀搜索算法(Sparrow Search Algorithm,SSA)優(yōu)化的極限學習機(Extreme Learning Machine,ELM)網(wǎng)絡(luò),對特征向量集進行識別和故障狀態(tài)診斷,實現(xiàn)對航空液壓管路的故障識別及診斷。

    1 基本理論

    1.1 變分模態(tài)分解

    變分模態(tài)分解算法(Variational Mode Decomposition,VMD)是一種能夠自適應(yīng)復雜信號且能夠?qū)⒕哂形锢硪饬x的本征模態(tài)分量分解出來的信號處理方法。其主要原理是通過對待處理信號進行受條件約束的變分模型構(gòu)造,并采用迭代方式搜尋約束變分模型的最優(yōu)解,從而獲得每個本征模態(tài)函數(shù)分量的中心頻率和帶寬,進而實現(xiàn)對信號的最佳分解。其中變分模型構(gòu)造是對待分解信號()進行分解,獲得個模態(tài)分量信號()以及估計中心頻率,并對各個模態(tài)分量信號進行希爾伯特變換求取解析信號,與估計中心頻率混合,在滿足各分量信號之和等于原信號的約束條件下,建立變分模型如式(1)所示:

    (1)

    在模型求解過程中,通過引入拉格朗日乘法算子和懲罰因子,將模型轉(zhuǎn)換成不受條件約束的變分模型,如式(2)所示:

    ({},{},)=

    (2)

    利用交替方向乘法算子,通過迭代更新的方式得到變分模型的鞍點,即受約束變分模型最優(yōu)解。

    具體迭代更新過程如下:

    (3)

    (4)

    (5)

    (6)

    1.2 麻雀搜索算法

    麻雀搜索算法是一種模仿自然界中麻雀生活中覓食和躲避捕食者過程行為方式的群體尋找最優(yōu)解的智能優(yōu)化算法模型,在計算時不需要依賴復雜的梯度信息,僅需要對目標函數(shù)和適應(yīng)度函數(shù)進行設(shè)定。在求解最優(yōu)問題時,該算法具有很強的實用性和較好的魯棒性等特點。具體工作流程如圖1所示。

    圖1 麻雀搜索算法流程

    1.3 優(yōu)化多尺度特征提取

    (7)

    (8)

    (9)

    式中:均值和標準差取自時間序列。

    (10)

    式中:int為取整函數(shù);為類別數(shù)。

    (11)

    式中:為嵌入維數(shù);為時間延遲。

    (12)

    (13)

    通過分析航空液壓管路振動信號發(fā)現(xiàn),傳統(tǒng)的信號特征分析方法只從單一時間尺度方面考慮,所分析的信號特征信息不夠全面,未能較好地表征出液壓管路的狀態(tài)。因此,從多個時間尺度進行信號的特征分析,能夠找到更多隱藏的特征信息。

    多尺度散布熵MDE能夠描述信號序列在不同尺度因子下的熵值特征信息,MDE的函數(shù)表達式如式(14)所示:

    (14)

    優(yōu)化多尺度散布熵IMDE是基于多尺度散布熵MDE的一種優(yōu)化方法,它能夠避免傳統(tǒng)多尺度粗?;^程中造成的鄰近序列信息遺漏,同時起到抑制因時間尺度因子增大序列數(shù)減少而導致的熵值突變現(xiàn)象,在一定程度上能提升熵值計算的穩(wěn)定性。通過分別計算優(yōu)化粗粒序列在不同時間尺度因子下的散布熵值的均值,即可得到IMDE的函數(shù)表達式如(15)所示:

    (15)

    1.4 優(yōu)化極限學習機

    極限學習機ELM是一種單隱含層前饋神經(jīng)網(wǎng)絡(luò)的算法,通過設(shè)置隱含層神經(jīng)元個數(shù),同時隨機初始化輸入權(quán)值和隱含層閾值,就能夠確定網(wǎng)絡(luò)結(jié)構(gòu)。與傳統(tǒng)神經(jīng)網(wǎng)絡(luò)相比,ELM需要設(shè)置的參數(shù)少,在保證精度的前提下,ELM因具有良好的泛化性和學習的高效性被廣泛應(yīng)用于各個工程領(lǐng)域。

    由于極限學習機神經(jīng)網(wǎng)絡(luò)的輸入權(quán)值與隱含層閾值是隨機產(chǎn)生的,且在訓練過程中不發(fā)生改變,在一定程度上影響了性能,因此選擇合適的參數(shù)是提升極限學習機網(wǎng)絡(luò)性能的關(guān)鍵。

    為提高極限學習機的性能,利用麻雀搜索算法,將樣本訓練集與測試集的分類誤差率之和作為適應(yīng)度函數(shù),對極限學習機網(wǎng)絡(luò)的輸入權(quán)值和隱含層閾值參數(shù)進行優(yōu)化選取,具體流程如圖2所示。

    圖2 優(yōu)化極限學習機流程

    2 故障診斷流程

    多尺度特征組合優(yōu)化的航空液壓管路故障診斷方法流程如圖3所示,詳細步驟如下:

    圖3 航空液壓管路故障診斷流程

    (1)按照一定的采樣頻率,分別采集航空液壓管路正常狀態(tài)、裂紋故障狀態(tài)、凹坑故障狀態(tài)的振動信號;

    (2)分別對管路不同狀態(tài)下的振動信號進行VMD分解,利用能量比值法確定適應(yīng)該信號的VMD的最佳分解層數(shù);

    (3)根據(jù)最佳分解層數(shù)進行VMD分解,獲得模態(tài)分量信號,分別計算每個模態(tài)分量信號與原信號的相關(guān)系數(shù)值,并根據(jù)相關(guān)系數(shù)值選取合適的模態(tài)分量信號進行重構(gòu);

    (4)計算每種狀態(tài)下重構(gòu)信號的優(yōu)化多尺度散布熵,選取合適時間尺度因子下的散布熵作為故障特征指標組成數(shù)據(jù)集;

    (5)將故障特征數(shù)據(jù)集代入麻雀搜索算法優(yōu)化后的極限學習機網(wǎng)絡(luò)進行訓練,得到故障診斷模型,實現(xiàn)管路故障識別分類。

    3 信號分析處理

    3.1 仿真信號分析

    為驗證所提方法的有效性,給出如式(16)所示的仿真信號模型:

    (16)

    其中:=50 Hz、=150 Hz、=250 Hz、=350 Hz;為高斯白噪聲。采樣頻率為1 000 Hz;采樣點數(shù)為1 000。生成的仿真信號時頻域波形如圖4所示。

    圖4 仿真信號時頻域波形

    首先,采用較通用的懲罰因子=2 000作為因子值,對仿真信號進行變分模態(tài)分解。分別計算出仿真信號在為2~7之間時各模態(tài)分量的能量比,結(jié)果如圖5所示。可知:當由2變化到4時,能量比大幅度增大;當從4變化到5時,未發(fā)生大幅度增大且趨于穩(wěn)定,說明當=4時達到最佳分解個數(shù),如果繼續(xù)增大會導致信號過度分解,造成模式混疊現(xiàn)象。

    圖5 仿真信號各模態(tài)分量能量比值曲線

    通過中心頻率觀察法對能量比值法判別選取的正確性進行驗證,計算仿真信號在不同下各個模態(tài)分量的模態(tài)中心頻率并進行分析。分別計算為2~7時各模態(tài)分量的中心頻率,結(jié)果如圖6所示??芍簩Ρ葹?~3時的分解結(jié)果,能夠發(fā)現(xiàn)各模態(tài)之間的中心頻率差距過大,信號分解明顯不足;當=5時,IMF4和IMF5模態(tài)分量的中心頻率較為接近,由此推斷信號分解過程中,出現(xiàn)模態(tài)混疊現(xiàn)象。為了防止過度分解出現(xiàn)模態(tài)混疊現(xiàn)象,綜合考慮=4符合該信號模態(tài)的最佳分解個數(shù)。

    圖6 仿真信號各模態(tài)分量中心頻率

    3.2 管路實驗

    通過模擬航空液壓管路工作環(huán)境,對不同狀態(tài)類型的液壓管路分別進行實驗。采取人工植入的方式,分別對液壓直管和彎管的管體表面設(shè)置裂紋和凹坑兩種故障。

    圖7所示為航空液壓管路實驗平臺示意,主要由電氣控制系統(tǒng)、液壓動力系統(tǒng)、實驗臺和信號采集系統(tǒng)組成。

    圖7 航空液壓管路實驗平臺示意

    實驗中所涉及參數(shù):液壓系統(tǒng)壓力為10 MPa,電動機轉(zhuǎn)速為1 200 r/min,采樣頻率為7 680 Hz,使用加速度傳感器采集振動數(shù)據(jù),不同振動測試管路不同傳感器布局位置如圖8所示。

    圖8 液壓管路振動實驗

    3.3 實驗信號分析

    通過實驗獲得了液壓管路不同狀態(tài)類型的振動信號,文中以直管裂紋故障狀態(tài)類型的振動信號為例進行分析。利用變分模態(tài)分解算法處理振動信號,其中默認懲罰因子=2 000,利用能量比值法確定分解層數(shù)。圖9所示為不同時直管裂紋故障狀態(tài)振動信號的能量比值曲線,通過觀察可以選取直管裂紋故障狀態(tài)的振動信號分解層數(shù)=6。

    圖9 直管裂紋故障各模態(tài)分量能量比值曲線

    對直管裂紋的振動信號進行分解,得到模態(tài)分量信號如圖10所示。

    圖10 液壓直管裂紋故障狀態(tài)振動信號VMD分解結(jié)果

    為盡可能消除干擾信號,獲得信號中重要的信息成分,剔除與原信號相關(guān)性較低的分量,選取相關(guān)性較高的信號重構(gòu)。計算直管裂紋故障狀態(tài)下,分解后的各個模態(tài)分量信號IMF與原信號的相關(guān)系數(shù),結(jié)果如表1所示。

    表1 直管裂紋故障狀態(tài)模態(tài)分量相關(guān)系數(shù)

    模態(tài)分量信號的相關(guān)系數(shù)越大,說明它與原信號之間的相關(guān)性越高,分量信號中包含與原信號有關(guān)的信息越多。當相關(guān)系數(shù)小于1/時,表示相關(guān)性較差,因此選用與原信號相關(guān)性較高的模態(tài)分量IMF1和IMF2進行重構(gòu),得到的重構(gòu)模態(tài)分量作為進一步多尺度特征分析的對象。

    根據(jù)上述方法依次計算各項數(shù)據(jù)的重構(gòu)模態(tài)分量,并分別計算其優(yōu)化多尺度散布熵作為特征指標。其中,計算多尺度散布熵所涉及的參數(shù)嵌入維數(shù)=2、類別數(shù)=4、時間延遲=1、時間尺度=20。

    標準差可以較好地衡量離散程度,散布熵的標準差曲線如圖11所示??芍褐惫軣o故障狀態(tài)和直管裂紋故障狀態(tài)數(shù)據(jù)集的散布熵標準差較低,說明重構(gòu)信號隨著時間尺度因子的增大,整體熵值較為平穩(wěn);直管凹坑故障狀態(tài)數(shù)據(jù)集的散布熵標準差隨著的增大而遞增,說明信號整體較不穩(wěn)定,在大于8的較高時間尺度時出現(xiàn)交叉疏散現(xiàn)象。彎管無故障狀態(tài)和彎管裂紋故障狀態(tài)以及彎管凹坑故障狀態(tài)數(shù)據(jù)集的散布熵標準差曲線具有一定近似性,隨著時間尺度因子的改變逐漸產(chǎn)生小幅度變化。在選擇信號的特征指標時,特征數(shù)量太多會導致信息冗余,特征數(shù)量太少又不能完全反映信息。從整體角度觀察其多尺度散布熵的標準差曲線,經(jīng)綜合考慮,選取為1~8的時間尺度的散布熵值作為特征指標。

    圖11 不同狀態(tài)管路改進多尺度散布熵標準差曲線

    3.4 狀態(tài)識別分類

    運用文中所述方法構(gòu)建多尺度特征向量集,液壓管路直管裂紋故障、直管凹坑故障、彎管裂紋故障、彎管凹坑故障的標簽分別設(shè)置為a、b、c、d ,每種故障的振動信號各選取40組,共160組數(shù)據(jù),每組數(shù)據(jù)的長度為11 520個點。從每種狀態(tài)類型中隨機選取25組數(shù)據(jù)作為訓練樣本,共計100組,每種狀態(tài)類型剩余數(shù)據(jù)作為測試樣本。將訓練集樣本輸入到利用麻雀搜索算法優(yōu)化的極限學習機網(wǎng)絡(luò)模型進行訓練,進而實現(xiàn)故障識別及診斷。經(jīng)過多次嘗試,選取優(yōu)化算法的最佳參數(shù)為種群規(guī)模為50、最大迭代次數(shù)為100、警戒閾值為0.6。

    為驗證文中方法的優(yōu)越性,探討多尺度特征提取對識別效果的影響以及變分模態(tài)分解是否對管路狀態(tài)的識別結(jié)果產(chǎn)生影響。選用排列熵PE作為特征指標和采用支持向量機SVM網(wǎng)絡(luò)等方法作為識別模型,并選用相同振動信號數(shù)據(jù)進行對比分析,結(jié)果如表2所示。

    由表2可知:選用傳統(tǒng)多尺度散布熵MDE作為特征指標輸入到ELM網(wǎng)絡(luò)模型中的方法(VMD-MDE-ELM),分類準確率能達到95%,與單尺度特征的方法(VMD-DE-ELM)相比,分類準確率提高13.33%;利用傳統(tǒng)多尺度散布熵MDE作為特征指標輸入到極限學習機的方法(VMD-MDE-ELM),分類準確率低于優(yōu)化多尺度IMDE,說明優(yōu)化多尺度粗?;^程能夠提高網(wǎng)絡(luò)模型分類準確率。

    表2 不同方法識別結(jié)果

    采用直接計算原始振動信號的排列熵PE為特征值的極限學習機方法(PE-ELM)分類準確率僅達到73.33%,利用散布熵DE作為特征指標時(DE-ELM)的分類效果準確率能達到80%,而經(jīng)過變分模態(tài)分解處理的方法VMD-PE-ELM和VMD-DE-ELM的準確率分別能夠達到76.67%和81.67%,說明特征指標選取散布熵與排列熵更具有優(yōu)勢,驗證了變分模態(tài)分解處理對識別準確率具有一定影響。

    分別將優(yōu)化多尺度IMDE方法(VMD-IMDE- SVM)和傳統(tǒng)多尺度MDE方法(VMD-MDE-SVM)輸入到支持向量機網(wǎng)絡(luò)模型,其分類準確率分別為96.67%和93.33%,而VMD-IMDE-ELM方法能夠達到98.33%;與此同時,極限學習機網(wǎng)絡(luò)整體上運行所需時間較少。說明相同特征指標下極限學習機網(wǎng)絡(luò)比支持向量機網(wǎng)絡(luò)在運行速度和識別準確率方面擁有更顯著的優(yōu)勢。

    將優(yōu)化多尺度IMDE方法輸入到經(jīng)麻雀搜索算法優(yōu)化的極限學習機網(wǎng)絡(luò)模型中的方法VMD-IMDE-SSA-ELM的分類準確率能夠達到100%,說明對極限學習機網(wǎng)絡(luò)進行優(yōu)化,能提高其識別的綜合準確率。

    4 結(jié)論

    針對航空液壓管路故障特征難以提取、故障診斷困難的問題,提出了多尺度特征組合優(yōu)化的航空液壓管路故障診斷方法。經(jīng)過實驗驗證,所述方法能夠準確有效地實現(xiàn)對不同類型航空液壓管路的故障識別及診斷。主要結(jié)論如下:

    (1)利用能量比值法能夠?qū)ψ兎帜B(tài)分解的分解層數(shù)進行優(yōu)化選取;通過選取合適的分量信號進行重構(gòu)并進行多尺度特征提取,能夠提升網(wǎng)絡(luò)模型的識別準確率。

    (2)與傳統(tǒng)多尺度特征提取方法相比,優(yōu)化多尺度特征提取方法穩(wěn)定性更好;選擇散布熵作為特征指標具有更好的分類效果。

    (3)與支持向量機網(wǎng)絡(luò)相比,極限學習機網(wǎng)絡(luò)需要設(shè)置的參數(shù)更少且在運行速度和分類識別方面較優(yōu);利用麻雀搜索算法優(yōu)化后的極限學習機網(wǎng)絡(luò),識別準確率得到有效提升。

    猜你喜歡
    學習機管路尺度
    基于水質(zhì)變化的供熱采暖管路設(shè)計
    財產(chǎn)的五大尺度和五重應(yīng)對
    液壓管路系統(tǒng)隨機振動下疲勞分析
    硅鋼軋制過程中乳化液流量控制解耦研究及應(yīng)用
    山西冶金(2019年2期)2019-05-31 11:30:04
    極限學習機綜述
    基于極限學習機參數(shù)遷移的域適應(yīng)算法
    自動化學報(2018年2期)2018-04-12 05:46:21
    分層極限學習機在滾動軸承故障診斷中的應(yīng)用
    宇宙的尺度
    太空探索(2016年5期)2016-07-12 15:17:55
    9
    美航天服漏水或因管路堵塞
    太空探索(2014年4期)2014-07-19 10:08:58
    国产人伦9x9x在线观看| 亚洲人成77777在线视频| 欧美亚洲日本最大视频资源| 欧美中文综合在线视频| 精品国产乱码久久久久久男人| 日本免费一区二区三区高清不卡 | 国产成人精品无人区| 国产成人系列免费观看| 亚洲色图av天堂| 麻豆国产av国片精品| 视频区欧美日本亚洲| 久久久久久免费高清国产稀缺| a级毛片黄视频| 欧美日韩精品网址| 亚洲午夜精品一区,二区,三区| 岛国在线观看网站| 色婷婷久久久亚洲欧美| 老司机深夜福利视频在线观看| 女人爽到高潮嗷嗷叫在线视频| 午夜91福利影院| av电影中文网址| 无限看片的www在线观看| 国产免费男女视频| 国产精品国产高清国产av| 九色亚洲精品在线播放| 午夜视频精品福利| 国产野战对白在线观看| 精品久久久久久久毛片微露脸| 亚洲色图 男人天堂 中文字幕| 男女高潮啪啪啪动态图| 欧美不卡视频在线免费观看 | 操美女的视频在线观看| 国产亚洲精品一区二区www| 国产高清视频在线播放一区| 国产一区二区在线av高清观看| 久久久久九九精品影院| 国产在线精品亚洲第一网站| 亚洲熟女毛片儿| 亚洲精品久久成人aⅴ小说| 国产亚洲精品久久久久5区| 日本一区二区免费在线视频| 国产激情欧美一区二区| 成人永久免费在线观看视频| 宅男免费午夜| 亚洲国产精品999在线| 最好的美女福利视频网| 欧美日韩乱码在线| 中出人妻视频一区二区| 手机成人av网站| 亚洲精品一区av在线观看| 黄色视频不卡| 后天国语完整版免费观看| 久久影院123| 一二三四在线观看免费中文在| 国产不卡一卡二| 午夜亚洲福利在线播放| 亚洲第一欧美日韩一区二区三区| 又黄又爽又免费观看的视频| 80岁老熟妇乱子伦牲交| 亚洲五月天丁香| 老司机福利观看| 欧美乱色亚洲激情| 欧美黑人精品巨大| 色精品久久人妻99蜜桃| 两个人看的免费小视频| 国产无遮挡羞羞视频在线观看| 欧美激情 高清一区二区三区| 麻豆av在线久日| 日本三级黄在线观看| 亚洲午夜精品一区,二区,三区| 精品一区二区三卡| 亚洲国产欧美一区二区综合| 伦理电影免费视频| 黑人巨大精品欧美一区二区mp4| 男女做爰动态图高潮gif福利片 | 51午夜福利影视在线观看| 亚洲一区高清亚洲精品| 亚洲色图av天堂| 欧美黑人精品巨大| 欧美精品一区二区免费开放| 国产1区2区3区精品| 午夜视频精品福利| 亚洲狠狠婷婷综合久久图片| 国产97色在线日韩免费| 成年版毛片免费区| 岛国视频午夜一区免费看| 久久九九热精品免费| 人人澡人人妻人| 亚洲精品在线美女| 国产午夜精品久久久久久| 69精品国产乱码久久久| 国产熟女xx| 国产亚洲精品第一综合不卡| cao死你这个sao货| 久久久久久久久免费视频了| 午夜精品国产一区二区电影| 国产精品自产拍在线观看55亚洲| 欧美成人免费av一区二区三区| 老司机亚洲免费影院| 天天影视国产精品| 亚洲五月婷婷丁香| 色播在线永久视频| 一a级毛片在线观看| 人成视频在线观看免费观看| 精品一区二区三区av网在线观看| 国产欧美日韩一区二区三区在线| 99国产综合亚洲精品| 女人精品久久久久毛片| 老司机亚洲免费影院| avwww免费| 久久天躁狠狠躁夜夜2o2o| 国产人伦9x9x在线观看| 好男人电影高清在线观看| 国产精品自产拍在线观看55亚洲| 国产成人啪精品午夜网站| 欧美精品亚洲一区二区| 88av欧美| 99久久精品国产亚洲精品| 美国免费a级毛片| 天天躁狠狠躁夜夜躁狠狠躁| 亚洲一区二区三区欧美精品| 老司机福利观看| 国产成人精品久久二区二区免费| 日韩成人在线观看一区二区三区| 国产成年人精品一区二区 | 亚洲av成人不卡在线观看播放网| 高清黄色对白视频在线免费看| www.精华液| 亚洲中文字幕日韩| 日韩视频一区二区在线观看| www.熟女人妻精品国产| 丰满饥渴人妻一区二区三| 长腿黑丝高跟| 国产成人一区二区三区免费视频网站| 搡老熟女国产l中国老女人| 91成年电影在线观看| 欧美丝袜亚洲另类 | 黑人巨大精品欧美一区二区蜜桃| 成人亚洲精品一区在线观看| 欧洲精品卡2卡3卡4卡5卡区| 少妇的丰满在线观看| 精品高清国产在线一区| 另类亚洲欧美激情| 99久久国产精品久久久| 激情在线观看视频在线高清| 日本黄色日本黄色录像| 99久久综合精品五月天人人| 久久婷婷成人综合色麻豆| 在线免费观看的www视频| 国产精品久久久久久人妻精品电影| 18禁黄网站禁片午夜丰满| xxxhd国产人妻xxx| 动漫黄色视频在线观看| 另类亚洲欧美激情| 日本a在线网址| 日本一区二区免费在线视频| 成人18禁高潮啪啪吃奶动态图| 亚洲午夜理论影院| 亚洲成人免费电影在线观看| 亚洲成人精品中文字幕电影 | 欧美精品一区二区免费开放| 老司机福利观看| 9191精品国产免费久久| 夜夜看夜夜爽夜夜摸 | av超薄肉色丝袜交足视频| 热99re8久久精品国产| 亚洲在线自拍视频| 国产成人精品无人区| 国产亚洲欧美98| 国产精品一区二区免费欧美| 女人爽到高潮嗷嗷叫在线视频| 不卡一级毛片| 亚洲片人在线观看| 国产精品av久久久久免费| 欧美色视频一区免费| 成人影院久久| 国产精品国产高清国产av| 午夜福利,免费看| 久久久精品国产亚洲av高清涩受| 一二三四在线观看免费中文在| 长腿黑丝高跟| 夫妻午夜视频| ponron亚洲| 国产伦人伦偷精品视频| xxxhd国产人妻xxx| 日本 av在线| 成人黄色视频免费在线看| 午夜91福利影院| 婷婷六月久久综合丁香| 精品一区二区三区四区五区乱码| 国产亚洲精品久久久久久毛片| av有码第一页| 天堂影院成人在线观看| 在线视频色国产色| 久久久久国内视频| 欧美乱色亚洲激情| 国产三级黄色录像| 亚洲精品国产色婷婷电影| 老熟妇仑乱视频hdxx| 中文字幕人妻丝袜制服| 亚洲色图综合在线观看| 亚洲国产精品999在线| 嫁个100分男人电影在线观看| 真人做人爱边吃奶动态| 国产精品国产av在线观看| 欧美成狂野欧美在线观看| 91在线观看av| 老汉色av国产亚洲站长工具| 美女扒开内裤让男人捅视频| 神马国产精品三级电影在线观看 | 亚洲五月色婷婷综合| 色综合站精品国产| 久久热在线av| 精品人妻1区二区| 中亚洲国语对白在线视频| 国产1区2区3区精品| 精品一区二区三卡| 性少妇av在线| 国产一区二区三区视频了| 叶爱在线成人免费视频播放| 我的亚洲天堂| 亚洲国产精品sss在线观看 | 女同久久另类99精品国产91| 国产伦人伦偷精品视频| 一区二区三区国产精品乱码| 宅男免费午夜| 亚洲avbb在线观看| 欧美日本亚洲视频在线播放| 亚洲一卡2卡3卡4卡5卡精品中文| 一二三四在线观看免费中文在| videosex国产| 亚洲情色 制服丝袜| 在线视频色国产色| 国产欧美日韩一区二区精品| 一二三四在线观看免费中文在| 90打野战视频偷拍视频| 国产在线精品亚洲第一网站| 国产男靠女视频免费网站| 国产三级黄色录像| 天堂影院成人在线观看| 乱人伦中国视频| 成人免费观看视频高清| 免费久久久久久久精品成人欧美视频| 亚洲av成人不卡在线观看播放网| 国产黄色免费在线视频| 超碰成人久久| 国产黄a三级三级三级人| 精品午夜福利视频在线观看一区| 成人国产一区最新在线观看| 欧美乱妇无乱码| 欧美性长视频在线观看| 亚洲成av片中文字幕在线观看| av片东京热男人的天堂| 三上悠亚av全集在线观看| 日韩欧美在线二视频| 老熟妇乱子伦视频在线观看| 成在线人永久免费视频| av国产精品久久久久影院| av在线播放免费不卡| 波多野结衣一区麻豆| 亚洲国产精品合色在线| 久久人妻福利社区极品人妻图片| 一区二区三区国产精品乱码| 人人妻人人爽人人添夜夜欢视频| 高清av免费在线| 女人被躁到高潮嗷嗷叫费观| 中文欧美无线码| 日韩精品中文字幕看吧| 香蕉丝袜av| 精品国产乱子伦一区二区三区| 国产真人三级小视频在线观看| 视频区图区小说| 国产精品亚洲一级av第二区| 免费搜索国产男女视频| av电影中文网址| 亚洲精品美女久久久久99蜜臀| 亚洲欧美精品综合久久99| 大型黄色视频在线免费观看| 久9热在线精品视频| 国产精品免费视频内射| 人成视频在线观看免费观看| 久久精品亚洲熟妇少妇任你| 大陆偷拍与自拍| 免费观看精品视频网站| a级片在线免费高清观看视频| 女人被狂操c到高潮| 老熟妇仑乱视频hdxx| 日韩大码丰满熟妇| 少妇 在线观看| 中文字幕另类日韩欧美亚洲嫩草| 黑人操中国人逼视频| 88av欧美| 国产精品一区二区在线不卡| 黄色女人牲交| 国产精品久久电影中文字幕| 精品国产一区二区三区四区第35| 国产成人欧美在线观看| av天堂在线播放| 波多野结衣一区麻豆| 一级黄色大片毛片| 一个人观看的视频www高清免费观看 | 亚洲人成77777在线视频| 桃色一区二区三区在线观看| 国产片内射在线| 久热爱精品视频在线9| 国产欧美日韩一区二区三| 亚洲人成电影观看| 人人澡人人妻人| 午夜视频精品福利| 99国产综合亚洲精品| 一级片'在线观看视频| 亚洲精品一卡2卡三卡4卡5卡| 男女做爰动态图高潮gif福利片 | 十八禁网站免费在线| 国产麻豆69| 一边摸一边抽搐一进一出视频| 亚洲第一青青草原| 一进一出抽搐动态| 欧美人与性动交α欧美软件| 少妇被粗大的猛进出69影院| 最新在线观看一区二区三区| 亚洲av成人不卡在线观看播放网| 香蕉国产在线看| 国产精品免费视频内射| 亚洲五月天丁香| 国产免费男女视频| 免费女性裸体啪啪无遮挡网站| 欧美最黄视频在线播放免费 | 亚洲精品中文字幕一二三四区| 91老司机精品| 色在线成人网| 日本免费a在线| 高清av免费在线| 亚洲熟女毛片儿| 色综合站精品国产| 午夜久久久在线观看| 超碰97精品在线观看| 久久精品人人爽人人爽视色| 欧美一级毛片孕妇| 动漫黄色视频在线观看| 免费高清在线观看日韩| 国产激情久久老熟女| 麻豆久久精品国产亚洲av | 好看av亚洲va欧美ⅴa在| 精品一品国产午夜福利视频| 国产av又大| 免费在线观看视频国产中文字幕亚洲| 亚洲视频免费观看视频| 日韩欧美在线二视频| 精品一区二区三卡| 免费久久久久久久精品成人欧美视频| avwww免费| www国产在线视频色| 欧美+亚洲+日韩+国产| 久久香蕉激情| 国产高清视频在线播放一区| 69精品国产乱码久久久| 18禁观看日本| 精品一区二区三区四区五区乱码| 日韩精品青青久久久久久| 成人黄色视频免费在线看| 亚洲成人免费电影在线观看| 免费看十八禁软件| 又黄又爽又免费观看的视频| 夜夜爽天天搞| 欧美日韩av久久| 免费人成视频x8x8入口观看| 亚洲精品一卡2卡三卡4卡5卡| 国产精品久久久久久人妻精品电影| 国产精品乱码一区二三区的特点 | 桃色一区二区三区在线观看| 天天躁狠狠躁夜夜躁狠狠躁| 老司机福利观看| 午夜影院日韩av| 国产成人一区二区三区免费视频网站| 亚洲全国av大片| 国产精品成人在线| 日韩中文字幕欧美一区二区| 欧美日韩视频精品一区| 999久久久国产精品视频| 日韩有码中文字幕| 欧美一区二区精品小视频在线| 麻豆国产av国片精品| 亚洲 欧美 日韩 在线 免费| 久久久国产成人免费| 亚洲五月色婷婷综合| 久久久久国内视频| 国产无遮挡羞羞视频在线观看| 如日韩欧美国产精品一区二区三区| 亚洲三区欧美一区| 国产熟女午夜一区二区三区| 亚洲人成网站在线播放欧美日韩| 19禁男女啪啪无遮挡网站| 美女扒开内裤让男人捅视频| 久久人人97超碰香蕉20202| 久久99一区二区三区| 日韩大尺度精品在线看网址 | 男女午夜视频在线观看| 国产精品 国内视频| 可以在线观看毛片的网站| 欧美在线黄色| 欧美日韩福利视频一区二区| 欧美国产精品va在线观看不卡| 亚洲av成人av| 亚洲男人的天堂狠狠| 人人妻,人人澡人人爽秒播| 黄片大片在线免费观看| 亚洲自偷自拍图片 自拍| www.999成人在线观看| 精品欧美一区二区三区在线| 纯流量卡能插随身wifi吗| 视频在线观看一区二区三区| 夜夜看夜夜爽夜夜摸 | 国产片内射在线| 1024香蕉在线观看| 亚洲精品久久午夜乱码| 午夜福利免费观看在线| 欧美人与性动交α欧美软件| 亚洲国产欧美一区二区综合| 久久精品91蜜桃| 国产熟女午夜一区二区三区| 亚洲伊人色综图| 操美女的视频在线观看| 男人舔女人下体高潮全视频| 亚洲精品国产色婷婷电影| 99热只有精品国产| 可以在线观看毛片的网站| 国产成人精品在线电影| 女人爽到高潮嗷嗷叫在线视频| 国产伦人伦偷精品视频| 国产精品免费视频内射| 国产精品成人在线| 97碰自拍视频| 97超级碰碰碰精品色视频在线观看| 亚洲国产精品999在线| 亚洲精品国产一区二区精华液| 他把我摸到了高潮在线观看| 欧美精品啪啪一区二区三区| 亚洲三区欧美一区| 亚洲一码二码三码区别大吗| 亚洲五月色婷婷综合| 在线观看免费日韩欧美大片| 国产高清视频在线播放一区| 亚洲免费av在线视频| 天天影视国产精品| 国产高清国产精品国产三级| 国产一区二区三区综合在线观看| 国产成人一区二区三区免费视频网站| 涩涩av久久男人的天堂| www.精华液| 欧美黄色片欧美黄色片| 可以在线观看毛片的网站| 黑丝袜美女国产一区| 中文字幕精品免费在线观看视频| 电影成人av| 9热在线视频观看99| 国产精品av久久久久免费| 久久国产亚洲av麻豆专区| 宅男免费午夜| 高清av免费在线| 91成人精品电影| 欧美人与性动交α欧美精品济南到| av免费在线观看网站| 中文字幕另类日韩欧美亚洲嫩草| 99国产精品免费福利视频| 欧洲精品卡2卡3卡4卡5卡区| 久久久国产欧美日韩av| 国产欧美日韩一区二区精品| 女同久久另类99精品国产91| 欧美日韩乱码在线| 99精品在免费线老司机午夜| 一区福利在线观看| 亚洲第一欧美日韩一区二区三区| 色哟哟哟哟哟哟| 黄频高清免费视频| 一级毛片精品| 日本黄色日本黄色录像| 久久午夜综合久久蜜桃| 国产主播在线观看一区二区| 桃色一区二区三区在线观看| 午夜精品国产一区二区电影| 如日韩欧美国产精品一区二区三区| 一边摸一边抽搐一进一出视频| 亚洲av第一区精品v没综合| 欧美日韩国产mv在线观看视频| 久久久久久大精品| 在线观看免费高清a一片| 国产亚洲精品久久久久5区| 国产精品免费视频内射| 日韩有码中文字幕| 99re在线观看精品视频| 9色porny在线观看| 黄色视频,在线免费观看| 欧美日韩瑟瑟在线播放| 成人手机av| 亚洲少妇的诱惑av| 男女高潮啪啪啪动态图| 久久久久久久午夜电影 | 亚洲国产精品一区二区三区在线| 丰满饥渴人妻一区二区三| 免费在线观看黄色视频的| 日韩成人在线观看一区二区三区| 成人影院久久| 国产一区二区三区在线臀色熟女 | 久久精品91无色码中文字幕| 亚洲专区中文字幕在线| 欧美日韩亚洲国产一区二区在线观看| 亚洲熟妇中文字幕五十中出 | 国产伦人伦偷精品视频| 国产精品乱码一区二三区的特点 | 一进一出抽搐动态| 免费看a级黄色片| 中文字幕色久视频| 欧美av亚洲av综合av国产av| 免费不卡黄色视频| 成人国产一区最新在线观看| 在线观看日韩欧美| 亚洲成人久久性| 中文欧美无线码| 夜夜看夜夜爽夜夜摸 | 黑人欧美特级aaaaaa片| 国产av一区二区精品久久| 精品久久久久久成人av| 欧美在线黄色| 中文字幕色久视频| 国产精品久久久av美女十八| 亚洲一区高清亚洲精品| 女人高潮潮喷娇喘18禁视频| 国产国语露脸激情在线看| 999精品在线视频| 桃红色精品国产亚洲av| 日本黄色视频三级网站网址| 久久久久九九精品影院| 亚洲自偷自拍图片 自拍| 啪啪无遮挡十八禁网站| 欧美日韩亚洲高清精品| 久久久国产一区二区| 中文字幕人妻丝袜一区二区| 一本大道久久a久久精品| 久久久久九九精品影院| 国产野战对白在线观看| av天堂久久9| 极品教师在线免费播放| 夜夜看夜夜爽夜夜摸 | 亚洲精品中文字幕在线视频| 欧美中文日本在线观看视频| 亚洲成av片中文字幕在线观看| 国产免费男女视频| 后天国语完整版免费观看| 在线观看一区二区三区激情| 免费观看精品视频网站| 99国产精品一区二区蜜桃av| 欧美一区二区精品小视频在线| 久久久国产一区二区| 国产午夜精品久久久久久| 欧美不卡视频在线免费观看 | 国产精品久久久av美女十八| 人人妻人人添人人爽欧美一区卜| 狠狠狠狠99中文字幕| 麻豆久久精品国产亚洲av | 日本精品一区二区三区蜜桃| 热99国产精品久久久久久7| 制服诱惑二区| 亚洲自拍偷在线| 黄片播放在线免费| 精品免费久久久久久久清纯| 12—13女人毛片做爰片一| 成熟少妇高潮喷水视频| 女性生殖器流出的白浆| 久热爱精品视频在线9| 色婷婷av一区二区三区视频| 在线免费观看的www视频| 超碰成人久久| 在线观看日韩欧美| 午夜福利一区二区在线看| 色综合婷婷激情| 亚洲人成电影观看| 成人三级黄色视频| 一级毛片精品| 亚洲一区二区三区欧美精品| 午夜影院日韩av| 精品国产乱子伦一区二区三区| 国产成人精品无人区| 男人舔女人的私密视频| 老司机午夜福利在线观看视频| 久久精品aⅴ一区二区三区四区| 成人影院久久| 国产精品野战在线观看 | 国产黄色免费在线视频| 日本免费一区二区三区高清不卡 | 欧美中文日本在线观看视频| 国产精品国产av在线观看| 久久 成人 亚洲| 性色av乱码一区二区三区2| 国产色视频综合| 亚洲国产毛片av蜜桃av| 可以在线观看毛片的网站| 狠狠狠狠99中文字幕| 在线观看免费视频日本深夜| 亚洲色图 男人天堂 中文字幕| 久久久久亚洲av毛片大全| 久久精品国产99精品国产亚洲性色 | 国产色视频综合| 婷婷精品国产亚洲av在线| 可以免费在线观看a视频的电影网站| 亚洲国产毛片av蜜桃av| 亚洲精品国产色婷婷电影| 国产亚洲欧美在线一区二区| 99国产精品一区二区三区| 久久久久国内视频| 日本免费一区二区三区高清不卡 | 日韩一卡2卡3卡4卡2021年| 久久国产精品男人的天堂亚洲| 国产日韩一区二区三区精品不卡| av有码第一页|