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

    多尺度能量熵與優(yōu)化極限學(xué)習(xí)機(jī)的航空液壓管路故障診斷方法

    2022-08-08 06:22:56薛政坤于曉光張小龍
    液壓與氣動(dòng) 2022年7期
    關(guān)鍵詞:特征向量管路尺度

    薛政坤,汪 曦,于曉光,王 寵,張小龍

    (遼寧科技大學(xué) 機(jī)械工程與自動(dòng)化學(xué)院,遼寧 鞍山 114051)

    引言

    航空液壓管路在工作中受振動(dòng)影響出現(xiàn)自身疲勞,以及過載和碰撞等因素易造成液壓管路的累積性損傷進(jìn)而導(dǎo)致液壓系統(tǒng)失效故障[1-3]。振動(dòng)引起的航空液壓管路故障一直是影響航空液壓系統(tǒng)可靠性的重要問題之一[4-5],因此開展航空液壓管路故障診斷研究,對(duì)提高航空液壓系統(tǒng)的可靠性具有重要意義。

    航空液壓管路在航空液壓系統(tǒng)中具有一定復(fù)雜性,導(dǎo)致管路振動(dòng)信號(hào)存在非平穩(wěn)以及非線性等特點(diǎn),對(duì)航空液壓管路的故障特征提取存在一定困難。當(dāng)前,針對(duì)處理非平穩(wěn)、非線性信號(hào)主要的分析處理方法為經(jīng)驗(yàn)?zāi)B(tài)分解(Empirical Mode Decomposition,EMD)[6],但EMD過程迭代次數(shù)較多,存在明顯的端點(diǎn)效應(yīng)。在經(jīng)驗(yàn)?zāi)B(tài)分解基礎(chǔ)上SMITH J S[7]提出了局部均值分解(Local Mean Decomposition,LMD),相比于EMD,LMD迭代次數(shù)更少,且一定程度上對(duì)端點(diǎn)效應(yīng)有一定的抑制。為了更好的提取非線性振動(dòng)信號(hào)中的特征信息,將信息熵作為特征信息的方式被引入到故障診斷領(lǐng)域中,同時(shí)得到了較廣泛的應(yīng)用[8]。師沖等[9]將EMD方法與一維密集卷積神經(jīng)網(wǎng)絡(luò)方法結(jié)合,解決了電液換向閥內(nèi)泄漏故障難以識(shí)別診斷的問題。陳敏捷等[10]運(yùn)用仿真模擬故障方式提取數(shù)據(jù),利用BP神經(jīng)網(wǎng)絡(luò)成功實(shí)現(xiàn)對(duì)液壓支架液壓系統(tǒng)的不同故障類型準(zhǔn)確診斷。

    為了從管路振動(dòng)信號(hào)中提取到更具有表征性的特征信息,更好的提高模型的故障診斷效果,本研究利用LMD對(duì)信號(hào)進(jìn)行預(yù)處理,分解出若干個(gè)PF分量信號(hào),利用相關(guān)系數(shù)-能量比準(zhǔn)則篩選出品質(zhì)最佳的PF分量信號(hào),計(jì)算最佳PF分量信號(hào)的多尺度能量熵(Multi-scale Energy Entropy,MEE),將其熵值組合成具有表征性的特征向量,并輸入到經(jīng)麻雀搜索算法(Sparrow Search Algorithm,SSA)優(yōu)化的極限學(xué)習(xí)機(jī)(Extreme Learning Machine,ELM)網(wǎng)絡(luò)學(xué)習(xí)訓(xùn)練,利用ELM對(duì)特征向量進(jìn)行分類,從而對(duì)航空液壓管路實(shí)現(xiàn)有效的故障診斷。

    1 基本理論

    1.1 LMD

    LMD[11]是一種自適應(yīng)信號(hào)處理方法,能夠?qū)⒂蓮?fù)雜的調(diào)幅調(diào)頻信號(hào)組成的信號(hào)分解為若干個(gè)PF分量信號(hào),本質(zhì)是依靠移動(dòng)平均法平滑處理待處理信號(hào)極值點(diǎn)的均值和包絡(luò)估計(jì)值,以獲得局部均值函數(shù)和包絡(luò)估計(jì)函數(shù),從而將純調(diào)頻信號(hào)和包絡(luò)子信號(hào)從待處理信號(hào)中分離出來,進(jìn)而通過迭代循環(huán)解調(diào)的方式得到純調(diào)頻信號(hào)和包絡(luò)信號(hào),其中包絡(luò)信號(hào)由迭代過程產(chǎn)生的各個(gè)包絡(luò)子信號(hào)相乘獲得,最終將純調(diào)頻信號(hào)和包絡(luò)信號(hào)相乘獲得PF分量信號(hào)。信號(hào)x(t)LMD流程圖如圖1所示[12-13],信號(hào)x(t)最終分解為N項(xiàng)PF分量和單調(diào)函數(shù)uN(t),即式(1)所示:

    (1)

    1.2 相關(guān)系數(shù)-能量比準(zhǔn)則

    信號(hào)經(jīng)LMD處理后分解成多個(gè)PF分量,由于分解后信號(hào)頻率由高頻到低頻產(chǎn)生,其幅值和頻率都會(huì)有所衰減,能量逐漸降低,且容易在低頻區(qū)產(chǎn)生虛假分量,因此提出相關(guān)系數(shù)-能量比準(zhǔn)則對(duì)各分量進(jìn)行篩選。

    能量比系數(shù)[14]表達(dá)式:

    (2)

    式中,ε——能量比系數(shù)

    EPF——各個(gè)PF分量具有的能量值

    E——總能量值

    局部均值分解形成的PF分量能量比值是不同的,單純采用能量比系數(shù)作為判定依據(jù)存在一定誤差,因此通過能量比與相關(guān)系數(shù)值相結(jié)合來綜合分析,將其定義為PF分量的Cc值,如式(3)所示:

    Cc=αC+(1-α)ε

    (3)

    式中,C——PF分量的相關(guān)系數(shù)值

    ε——PF分量的能量比系數(shù)

    α——PF分量的相關(guān)系數(shù)值對(duì)應(yīng)Cc值的權(quán)值

    1.3 MEE

    不同類型的故障產(chǎn)生的振動(dòng)信號(hào)頻率不同,不同的信號(hào)能量之間也會(huì)存在差異,因此可以通過能量熵反映信號(hào)的特征,能量熵[15]的計(jì)算公式為式(4):

    (4)

    (5)

    (6)

    式中,pi——信號(hào)的能量比值

    Ei——不同位置上具有的幅值能量

    MEE是基于能量熵的一種改進(jìn)方法,能夠描述信號(hào)的時(shí)間序列在不同尺度因子下能量化信號(hào)中包含信息的價(jià)值,具有較好的表征性。

    MEE的計(jì)算步驟如下:

    (7)

    (2) 設(shè)原始數(shù)據(jù)為長(zhǎng)度為N,計(jì)算原始數(shù)據(jù)時(shí)間序列經(jīng)粗?;玫溅觽€(gè)粗粒化序列的能量熵,構(gòu)造關(guān)于尺度因子τ的函數(shù),即可得到MEE的函數(shù)表達(dá)式:

    (8)

    1.4 SSA

    SSA是一種模仿自然界中麻雀覓食和躲避捕食者過程行為方式的群體尋找最優(yōu)解的智能優(yōu)化算法模型[17],在計(jì)算時(shí)不需要依賴復(fù)雜的梯度信息,僅需要對(duì)目標(biāo)函數(shù)和適應(yīng)度函數(shù)進(jìn)行設(shè)定,針對(duì)求解最優(yōu)問題時(shí),具有很強(qiáng)的實(shí)用性和較好的魯棒性等特點(diǎn)。

    在SSA中,麻雀?jìng)€(gè)體主要分為發(fā)現(xiàn)者、跟隨者、偵查者3種身份,且每只麻雀的身份能夠互相轉(zhuǎn)換。其中,發(fā)現(xiàn)者和跟隨者的數(shù)量維持固定比例,且具有較好的適應(yīng)度值的麻雀成為偵查者,其余的麻雀成為跟隨者。無(wú)論是發(fā)現(xiàn)者還是跟隨者,都有可能成為偵查者,但偵查者的數(shù)量是一定的。偵查者的主要工作是負(fù)責(zé)警戒偵查是否有危險(xiǎn),有危險(xiǎn)時(shí)放棄當(dāng)前位置,將食物移動(dòng)到安全位置,群體中間的麻雀則隨機(jī)移動(dòng)接近其他麻雀。

    其中,發(fā)現(xiàn)者的個(gè)體位置更新通過式(9)實(shí)現(xiàn):

    (9)

    跟隨者的個(gè)體位置更新通過式(10)實(shí)現(xiàn):

    (10)

    同時(shí)將從種群中隨機(jī)選擇一定比例個(gè)體作為偵查者,根據(jù)式(11)更新個(gè)體位置,經(jīng)過若干次迭代后完成優(yōu)化。

    (11)

    圖2 SSA執(zhí)行流程圖Fig.2 Execution flow chart of SSA

    1.5 ELM

    ELM是一種單隱含層前饋神經(jīng)網(wǎng)絡(luò)的算法[18],與傳統(tǒng)神經(jīng)網(wǎng)絡(luò)相比,需要設(shè)置的參數(shù)要求少,在保證精度的前提下,具有良好的泛化性和快速的學(xué)習(xí)速度,被廣泛應(yīng)用于諸多工程領(lǐng)域,該算法隨機(jī)初始化輸入權(quán)重和隱含層偏置,通過廣義逆矩陣?yán)碚撚?jì)算得到相應(yīng)的輸出權(quán)重。

    ELM神經(jīng)網(wǎng)絡(luò)在訓(xùn)練過程中不需要調(diào)整過多的參數(shù)值,只需確定隱含層神經(jīng)元的數(shù)量,便可以獲得唯一的最優(yōu)解。典型的ELM網(wǎng)絡(luò)結(jié)構(gòu)如圖3所示,網(wǎng)絡(luò)由輸入層和隱含層以及輸出層3部分構(gòu)成,且輸入層與隱含層以及隱含層與輸出層之間神經(jīng)元對(duì)應(yīng)逐個(gè)連接,網(wǎng)絡(luò)中有n個(gè)輸入層神經(jīng)元,l個(gè)隱含層神經(jīng)元,m個(gè)輸出層神經(jīng)元,網(wǎng)絡(luò)的輸出表達(dá)式為:

    圖3 ELM網(wǎng)絡(luò)結(jié)構(gòu)圖Fig.3 Network structure diagram of ELM

    (12)

    式中,g(x)為激活函數(shù),wi=[wi1,wi2,…,wij]表示輸入層中對(duì)應(yīng)的第i個(gè)神經(jīng)元與隱含層第j個(gè)神經(jīng)元之間的權(quán)重值;βi=[βj1,βj2,……,βjk]表示隱含層對(duì)應(yīng)的第j個(gè)神經(jīng)元與輸出層第k個(gè)神經(jīng)元之間的權(quán)重值;bl=[b1,b2,…,bl]表示第l隱含層的偏置。

    由于ELM神經(jīng)網(wǎng)絡(luò)的輸入權(quán)值與隱含層閾值通過隨機(jī)產(chǎn)生,且在訓(xùn)練過程中不發(fā)生改變,在一定程度上影響了性能[19]。為了提高極限學(xué)習(xí)機(jī)的性能,利用SSA根據(jù)樣本訓(xùn)練集與測(cè)試集的分類誤差率之和作為適應(yīng)度函數(shù),對(duì)ELM網(wǎng)絡(luò)的輸入權(quán)值和隱含層閾值參數(shù)進(jìn)行優(yōu)化選取,其流程圖如圖4所示。

    圖4 優(yōu)化ELM流程圖Fig.4 Flow chart of optimized ELM

    2 方法流程

    由于航空液壓系統(tǒng)結(jié)構(gòu)復(fù)雜,在實(shí)際工作中受到不同類型的振動(dòng)耦合,導(dǎo)致航空液壓管路的故障機(jī)理受多種因素影響,而且受外界強(qiáng)大噪聲干擾,故障特征難以提取。直接計(jì)算采集的原始振動(dòng)信號(hào)的能量熵作為故障特征不能很好的表示管路的故障特征。為了提高診斷結(jié)果,本研究提出了基于LMD與MEE的SSA與ELM相結(jié)合的航空液壓管路故障診斷方法,用于對(duì)航空液壓管路多故障狀態(tài)的分類診斷,航空液壓管路的故障診斷流程圖如圖5所示。

    圖5 航空液壓管路故障診斷流程Fig.5 Fault diagnosis process of aviation hydraulic pipeline

    基于LMD-MEE和SSA-ELM的航發(fā)液壓管路故障診斷方法如下:

    (1) 對(duì)采集來的液壓管路振動(dòng)信號(hào)進(jìn)行LMD,獲得PF分量信號(hào);

    (2) 針對(duì)PF分量信號(hào)分別計(jì)算其相關(guān)系數(shù)值和能量比值,并根據(jù)相關(guān)系數(shù)-能量比準(zhǔn)則選取最佳分量信號(hào);

    (3) 計(jì)算最佳分量信號(hào)的多尺度能量熵,選擇合適時(shí)間尺度下的能量熵值組成故障特征向量;

    (4) 將故障特征向量樣本帶入SSA優(yōu)化的ELM網(wǎng)絡(luò)訓(xùn)練,形成故障診斷模型,實(shí)現(xiàn)不同管路類型的故障識(shí)別診斷。

    3 實(shí)驗(yàn)分析

    3.1 振動(dòng)信號(hào)采集

    本實(shí)驗(yàn)使用如圖6所示的航空液壓實(shí)驗(yàn)裝置動(dòng)力系統(tǒng)進(jìn)行。實(shí)驗(yàn)裝置由電動(dòng)機(jī)、柱塞泵、冷卻器、壓力表、油箱、管路系統(tǒng)和控制系統(tǒng)等組成。

    圖6 航空液壓實(shí)驗(yàn)裝置動(dòng)力系統(tǒng)Fig.6 Power system of aviation hydraulic test device

    液壓系統(tǒng)的工作壓力為10 MPa,電動(dòng)機(jī)轉(zhuǎn)速為1200 r/min,采樣頻率為7680 Hz,使用加速度傳感器同步采集不同狀態(tài)管路的振動(dòng)數(shù)據(jù),實(shí)驗(yàn)過程中傳感器分布如圖7所示。

    圖7 液壓管路傳感器分布示意圖Fig.7 Distribution diagram of hydraulic pipeline sensor

    3.2 振動(dòng)信號(hào)分析

    對(duì)采集來的管路不同狀態(tài)下的振動(dòng)信號(hào)進(jìn)行局部均值分解,獲得若干個(gè)PF分量信號(hào)的波形圖,如圖8所示。

    圖8 液壓管路振動(dòng)信號(hào)PF分量波形圖Fig.8 PF component waveform of hydraulic pipeline vibration signal

    上述液壓管路振動(dòng)信號(hào)經(jīng)過局部均值分解,得到PF分量。運(yùn)用相關(guān)系數(shù)-能量比準(zhǔn)則,篩選其中分量得到其最佳分量。如表1所示,直管無(wú)故障的PF分量具有的最大Cc值約為0.8791,而其他分量的Cc值不大于0.5,因此選用PF2分量作為最佳分量;直管裂紋故障的PF分量具有的最大Cc值約為0.9076,而其他分量的Cc值不大于0.5,因此選用PF1分量作為最佳分量;直管凹坑故障的PF分量具有的最大Cc值約為0.8175,而其他分量的Cc值不大于0.5,因此選用PF1分量作為最佳分量;彎管無(wú)故障的PF分量具有的最大Cc值約為0.6839,而其他分量的Cc值不大于0.5,因此選用PF2分量作為最佳分量;彎管裂紋故障的PF分量具有的最大Cc值約為0.8023,而其他分量的Cc值不大于0.5,因此選用PF1分量作為最佳分量;彎管凹坑故障的PF分量具有的最大Cc值約為0.6600,而其他分量的Cc值不大于0.5,因此選用PF1分量作為最佳分量。

    表1 各PF分量相關(guān)系數(shù)值、能量比和Cc值Tab.1 Phase relation values,energy ratio and Cc values of PF components

    (續(xù)表1)

    將不同管路類型的最佳PF分量信號(hào)作為進(jìn)一步的分析信號(hào),計(jì)算不同類型液壓管路的最佳PF分量的多個(gè)時(shí)間尺度的能量熵。

    如圖9所示,為τ=1時(shí)的原始尺度能量熵S曲線,通過觀察能夠看出直管無(wú)故障與彎管無(wú)故障時(shí)能量熵值較為接近;直管裂紋故障與直管凹坑故障和彎管裂紋故障時(shí)的能量熵值差異較小,可分性較差;彎管無(wú)故障與彎管裂紋故障以及彎管凹坑故障之間能量熵值存在的差異明顯,可分性較強(qiáng)。

    圖9 原始尺度能量熵曲線Fig.9 Original scale energy entropy curve

    如圖10所示為多尺度能量熵SM曲線,隨著尺度因子τ不斷增大,液壓管路不同狀態(tài)下的能量熵之間的差值變大,但從整體角度觀察多尺度能量熵SM曲線,發(fā)現(xiàn)隨尺度因子τ增大,會(huì)出現(xiàn)一定程度的交叉干擾情況。

    圖10 多尺度能量熵曲線Fig.10 Multi-scale energy entropy curve

    由于多尺度能量熵相對(duì)于單一尺度能量熵能夠提升不同管路類型之間的可分性,因此選擇合適的時(shí)間尺度的能量熵值對(duì)特征向量的建立是非常重要的。在建立特征向量樣本數(shù)據(jù)集時(shí),特征數(shù)量太多會(huì)產(chǎn)生信息冗余,特征數(shù)量太少又不能完全反映信息。當(dāng)尺度因子τ為3~6時(shí),液壓管路不同狀態(tài)下的能量熵曲線變化較為規(guī)律且穩(wěn)定。綜合考慮,選取時(shí)間尺度3~6的能量熵值作為特征向量數(shù)據(jù)集。

    3.3 故障識(shí)別分類

    選取直管無(wú)故障、直管裂紋故障、直管凹坑故障、彎管無(wú)故障、彎管裂紋故障、彎管凹坑故障,6種不同類型的液壓管路狀態(tài)下振動(dòng)信號(hào)各40組數(shù)據(jù),共240組數(shù)據(jù)。隨機(jī)的從每種狀態(tài)中選取15組特征向量數(shù)據(jù)作為訓(xùn)練樣本,訓(xùn)練集共計(jì)90組數(shù)據(jù)。將每種狀態(tài)剩下的25組特征向量數(shù)據(jù)作為測(cè)試樣本,測(cè)試集共計(jì)150組數(shù)據(jù)。

    利用LMD對(duì)各個(gè)狀態(tài)下的振動(dòng)信號(hào)進(jìn)行分解,得到各個(gè)狀態(tài)下振動(dòng)信號(hào)對(duì)應(yīng)的若干個(gè)PF分量信號(hào),通過相關(guān)系數(shù)-能量比準(zhǔn)則篩選出各個(gè)狀態(tài)中的最佳PF分量,計(jì)算各個(gè)狀態(tài)下最佳PF分量的多尺度能量熵,從多尺度能量熵中選取同一尺度因子下的能量熵組成對(duì)應(yīng)各個(gè)狀態(tài)的特征向量數(shù)據(jù)集。

    根據(jù)上述方法,構(gòu)建特征向量樣本,輸入到經(jīng)SSA優(yōu)化的ELM網(wǎng)絡(luò)進(jìn)行訓(xùn)練,實(shí)現(xiàn)故障的識(shí)別分類。為了驗(yàn)證本研究所述方法的優(yōu)越性,選擇支持向量機(jī)網(wǎng)絡(luò)(Support Vector Machine,SVM)與之作對(duì)比,分類結(jié)果如表2所示。

    表2 故障識(shí)別分類結(jié)果Tab.2 Fault identification and classification results

    從結(jié)果來看,采用本研究所述方法LMD-MEE-SSA-ELM對(duì)狀態(tài)識(shí)別的分類準(zhǔn)確率達(dá)到100%,對(duì)航空液壓管路故障診斷識(shí)別分類具有較好的效果,與原始尺度能量熵作為特征向量樣本輸入到LMD-EE-SSA-ELM的方法相比,分類準(zhǔn)確率提高了15.33%,說明利用多尺度能量熵作特征向量相對(duì)于單一尺度原始能量熵作特征向量具有更好的識(shí)別效果。

    經(jīng)過與不同特征向量提取方法對(duì)比可以發(fā)現(xiàn),ELM網(wǎng)絡(luò)比支持向量機(jī)網(wǎng)絡(luò)在分類準(zhǔn)確度上具有更明顯的優(yōu)勢(shì),同時(shí)經(jīng)過SSA優(yōu)化后的ELM網(wǎng)絡(luò)模型的識(shí)別分類準(zhǔn)確率能夠得到較好地提升。

    4 結(jié)論

    針對(duì)航空液壓管路故障特征難以提取,故障診斷困難的問題,本研究提出了一種基于LMD與LMD和SSA優(yōu)化的ELM的航空液壓管路故障診斷方法;

    (1) 利用LMD能夠自適應(yīng)的將管路的振動(dòng)信號(hào)分解成若干個(gè)PF分量,利用相關(guān)系數(shù)-能量比準(zhǔn)則進(jìn)行篩選,有效抑制虛假分量干擾,提取含有重要信息的最佳分量;

    (2) MEE從多個(gè)時(shí)間尺度考慮了信號(hào)中能量的差異,相比能量熵作為特征向量,利用MEE作為特征向量在ELM分類中具有更高的識(shí)別準(zhǔn)確率;

    (3) 本研究所述方法能夠準(zhǔn)確有效的實(shí)現(xiàn)對(duì)不同類型航空液壓管路的故障識(shí)別診斷,同時(shí)發(fā)現(xiàn),ELM網(wǎng)絡(luò)比SVM網(wǎng)絡(luò)在分類識(shí)別準(zhǔn)確率方面擁有更好的效果,經(jīng)SSA優(yōu)化后的ELM網(wǎng)絡(luò),識(shí)別準(zhǔn)確率能夠得到有效地提升。

    猜你喜歡
    特征向量管路尺度
    二年制職教本科線性代數(shù)課程的幾何化教學(xué)設(shè)計(jì)——以特征值和特征向量為例
    基于水質(zhì)變化的供熱采暖管路設(shè)計(jì)
    克羅內(nèi)克積的特征向量
    財(cái)產(chǎn)的五大尺度和五重應(yīng)對(duì)
    液壓管路系統(tǒng)隨機(jī)振動(dòng)下疲勞分析
    硅鋼軋制過程中乳化液流量控制解耦研究及應(yīng)用
    山西冶金(2019年2期)2019-05-31 11:30:04
    一類特殊矩陣特征向量的求法
    EXCEL表格計(jì)算判斷矩陣近似特征向量在AHP法檢驗(yàn)上的應(yīng)用
    宇宙的尺度
    太空探索(2016年5期)2016-07-12 15:17:55
    9
    超碰av人人做人人爽久久| 日韩 亚洲 欧美在线| 校园春色视频在线观看| 脱女人内裤的视频| 动漫黄色视频在线观看| 他把我摸到了高潮在线观看| 欧美成人一区二区免费高清观看| 亚洲不卡免费看| 国产探花极品一区二区| 亚洲美女黄片视频| 俄罗斯特黄特色一大片| bbb黄色大片| 久久久久精品国产欧美久久久| 欧美激情国产日韩精品一区| 人人妻人人澡欧美一区二区| 麻豆国产97在线/欧美| 免费搜索国产男女视频| 亚洲第一欧美日韩一区二区三区| 国产精品亚洲美女久久久| 日韩人妻高清精品专区| 久久久久久国产a免费观看| 亚洲精品久久国产高清桃花| 久久伊人香网站| 少妇熟女aⅴ在线视频| 国产男靠女视频免费网站| 国产一区二区在线av高清观看| 十八禁人妻一区二区| 久久久国产成人精品二区| 国产精品人妻久久久久久| 一区二区三区高清视频在线| 欧美一区二区国产精品久久精品| bbb黄色大片| 美女被艹到高潮喷水动态| 国产午夜福利久久久久久| 国产免费av片在线观看野外av| 最近中文字幕高清免费大全6 | 校园春色视频在线观看| 在线观看免费视频日本深夜| 日本 av在线| 亚洲av第一区精品v没综合| 两人在一起打扑克的视频| 真人做人爱边吃奶动态| 国产av在哪里看| 午夜福利免费观看在线| 久9热在线精品视频| 国产av一区在线观看免费| 激情在线观看视频在线高清| 国产精品av视频在线免费观看| 亚洲欧美日韩卡通动漫| 国产精品伦人一区二区| 欧美绝顶高潮抽搐喷水| 亚洲国产精品sss在线观看| 露出奶头的视频| 精品一区二区三区视频在线| 免费观看精品视频网站| 亚洲av成人不卡在线观看播放网| 国产精品日韩av在线免费观看| 久99久视频精品免费| 啦啦啦韩国在线观看视频| 国产一级毛片七仙女欲春2| 亚洲,欧美精品.| 高清毛片免费观看视频网站| 免费人成在线观看视频色| 国内精品美女久久久久久| 美女高潮喷水抽搐中文字幕| 国产真实乱freesex| 看免费av毛片| 国产高清视频在线播放一区| a级毛片a级免费在线| 国产精品美女特级片免费视频播放器| 欧美潮喷喷水| 亚洲精品在线观看二区| 亚洲av一区综合| 美女高潮的动态| 国产午夜精品论理片| 亚洲精品一区av在线观看| 噜噜噜噜噜久久久久久91| 亚洲人成伊人成综合网2020| 欧美精品啪啪一区二区三区| 国产高清激情床上av| 欧美黑人巨大hd| 日韩欧美在线二视频| 国产乱人伦免费视频| 色综合亚洲欧美另类图片| 欧美xxxx黑人xx丫x性爽| 亚洲第一电影网av| 两个人的视频大全免费| 精品久久久久久久久av| 亚洲成a人片在线一区二区| 在线播放无遮挡| 国产三级在线视频| 天天躁日日操中文字幕| 九色国产91popny在线| 亚洲 欧美 日韩 在线 免费| 色综合站精品国产| 女人被狂操c到高潮| 又黄又爽又免费观看的视频| 91麻豆精品激情在线观看国产| 日韩高清综合在线| 日日夜夜操网爽| 国产精品影院久久| 亚洲国产精品成人综合色| 欧洲精品卡2卡3卡4卡5卡区| 亚洲中文字幕一区二区三区有码在线看| 哪里可以看免费的av片| 国内精品久久久久久久电影| 特级一级黄色大片| 国产一区二区在线观看日韩| 嫁个100分男人电影在线观看| 国产精品一区二区免费欧美| 亚洲成人久久性| 亚洲国产欧美人成| 免费在线观看成人毛片| 亚洲av五月六月丁香网| 欧美黑人欧美精品刺激| xxxwww97欧美| 女同久久另类99精品国产91| 国产亚洲欧美在线一区二区| 成年免费大片在线观看| 国产精品1区2区在线观看.| 我的女老师完整版在线观看| 好男人在线观看高清免费视频| 国产精品电影一区二区三区| 午夜免费成人在线视频| 亚洲精品久久国产高清桃花| 欧美高清成人免费视频www| 日本与韩国留学比较| 51午夜福利影视在线观看| 精品99又大又爽又粗少妇毛片 | 黄色女人牲交| 日韩亚洲欧美综合| 亚洲,欧美,日韩| 欧美bdsm另类| 最近在线观看免费完整版| 国产激情偷乱视频一区二区| 精品久久久久久久久亚洲 | 国内久久婷婷六月综合欲色啪| 午夜两性在线视频| 噜噜噜噜噜久久久久久91| 特大巨黑吊av在线直播| 欧美日韩黄片免| 琪琪午夜伦伦电影理论片6080| 成人鲁丝片一二三区免费| 亚洲国产精品合色在线| 亚洲综合色惰| 非洲黑人性xxxx精品又粗又长| 欧美乱色亚洲激情| 美女黄网站色视频| 国产视频内射| 国内精品美女久久久久久| 淫秽高清视频在线观看| 蜜桃久久精品国产亚洲av| 国产极品精品免费视频能看的| 色哟哟·www| 国内精品久久久久久久电影| 熟女电影av网| 床上黄色一级片| 两人在一起打扑克的视频| 国产黄片美女视频| 色综合站精品国产| 首页视频小说图片口味搜索| 亚洲av二区三区四区| 国产国拍精品亚洲av在线观看| 老司机午夜十八禁免费视频| 在线a可以看的网站| 亚洲av美国av| 麻豆国产av国片精品| 国产精品人妻久久久久久| 一a级毛片在线观看| 日本一二三区视频观看| 欧美乱色亚洲激情| 亚洲熟妇熟女久久| 丰满人妻熟妇乱又伦精品不卡| 亚洲一区二区三区色噜噜| or卡值多少钱| 黄色丝袜av网址大全| 一区二区三区高清视频在线| 最近最新免费中文字幕在线| 亚洲中文日韩欧美视频| 热99在线观看视频| 很黄的视频免费| 国产精品国产高清国产av| 国产精品一及| 亚洲欧美激情综合另类| 色5月婷婷丁香| 亚洲精品在线美女| 天美传媒精品一区二区| 亚洲不卡免费看| 国内精品一区二区在线观看| 国产精品三级大全| 丰满人妻熟妇乱又伦精品不卡| 欧美一区二区精品小视频在线| 国内毛片毛片毛片毛片毛片| 亚洲在线观看片| av中文乱码字幕在线| 久久亚洲真实| 国产激情偷乱视频一区二区| 国产精品免费一区二区三区在线| 国产成人aa在线观看| 国产精品一区二区三区四区免费观看 | 亚洲五月天丁香| 免费一级毛片在线播放高清视频| 国产免费av片在线观看野外av| 18禁黄网站禁片午夜丰满| 白带黄色成豆腐渣| 麻豆国产97在线/欧美| 成人亚洲精品av一区二区| 不卡一级毛片| 免费人成视频x8x8入口观看| 色综合婷婷激情| 婷婷六月久久综合丁香| 日本黄色片子视频| 1000部很黄的大片| 欧美日韩乱码在线| 亚洲欧美日韩高清在线视频| 国内揄拍国产精品人妻在线| 久久久国产成人免费| 黄色配什么色好看| 九色成人免费人妻av| 极品教师在线视频| 每晚都被弄得嗷嗷叫到高潮| 中亚洲国语对白在线视频| 真人做人爱边吃奶动态| 亚洲18禁久久av| 性欧美人与动物交配| 亚洲电影在线观看av| 性色av乱码一区二区三区2| 又黄又爽又刺激的免费视频.| 婷婷精品国产亚洲av在线| 欧美xxxx黑人xx丫x性爽| 国产亚洲精品av在线| 丰满人妻一区二区三区视频av| 国产日本99.免费观看| 在线观看66精品国产| 90打野战视频偷拍视频| 国产aⅴ精品一区二区三区波| 亚洲真实伦在线观看| www.www免费av| 欧美日本视频| 精品乱码久久久久久99久播| 国产v大片淫在线免费观看| ponron亚洲| 一区二区三区免费毛片| 免费看光身美女| 免费高清视频大片| 又爽又黄a免费视频| 欧美日本亚洲视频在线播放| 国产人妻一区二区三区在| 伦理电影大哥的女人| 综合色av麻豆| 啦啦啦观看免费观看视频高清| 国产真实伦视频高清在线观看 | 天堂网av新在线| av国产免费在线观看| 简卡轻食公司| 色哟哟哟哟哟哟| 久久久久久久精品吃奶| 国产精品99久久久久久久久| 亚洲最大成人av| 老女人水多毛片| 91av网一区二区| 别揉我奶头~嗯~啊~动态视频| 亚洲精品一卡2卡三卡4卡5卡| 亚洲最大成人av| 成人美女网站在线观看视频| 他把我摸到了高潮在线观看| 欧洲精品卡2卡3卡4卡5卡区| 免费看a级黄色片| 欧美乱妇无乱码| 欧美精品啪啪一区二区三区| 制服丝袜大香蕉在线| 夜夜躁狠狠躁天天躁| 热99在线观看视频| 校园春色视频在线观看| 精华霜和精华液先用哪个| 国产麻豆成人av免费视频| 国产老妇女一区| 制服丝袜大香蕉在线| 丁香六月欧美| 国产黄a三级三级三级人| 别揉我奶头 嗯啊视频| 色在线成人网| 国内精品久久久久久久电影| 成人国产一区最新在线观看| 国产欧美日韩一区二区精品| 日本黄色视频三级网站网址| 国产精品久久视频播放| h日本视频在线播放| 欧美性感艳星| 久久久久久久久中文| 欧美黄色片欧美黄色片| 久久精品国产亚洲av天美| 精品人妻视频免费看| 国产成人av教育| 亚洲电影在线观看av| 欧美黑人巨大hd| 一个人看视频在线观看www免费| 麻豆国产97在线/欧美| 国产色爽女视频免费观看| avwww免费| 又爽又黄无遮挡网站| 嫩草影视91久久| 真人一进一出gif抽搐免费| 国产伦一二天堂av在线观看| 国产一区二区三区在线臀色熟女| 天堂av国产一区二区熟女人妻| 亚洲自偷自拍三级| 在线播放国产精品三级| 久久精品国产自在天天线| 久久99热6这里只有精品| 免费在线观看影片大全网站| 给我免费播放毛片高清在线观看| 丁香六月欧美| 亚洲av成人av| 日韩欧美在线二视频| 丁香欧美五月| 国产成人欧美在线观看| 国产三级在线视频| 欧美3d第一页| 成年女人看的毛片在线观看| 一本久久中文字幕| 亚洲国产高清在线一区二区三| 天天躁日日操中文字幕| av福利片在线观看| 丝袜美腿在线中文| 国产精品国产高清国产av| 观看美女的网站| 亚洲成av人片在线播放无| 床上黄色一级片| 夜夜躁狠狠躁天天躁| 丰满人妻一区二区三区视频av| 日本在线视频免费播放| 精品国产三级普通话版| 99热这里只有精品一区| 少妇人妻一区二区三区视频| 亚洲一区高清亚洲精品| 亚洲精品在线观看二区| 亚洲成a人片在线一区二区| 97碰自拍视频| 少妇人妻精品综合一区二区 | 性插视频无遮挡在线免费观看| 天天一区二区日本电影三级| 日本a在线网址| 国产伦一二天堂av在线观看| 舔av片在线| 午夜免费男女啪啪视频观看 | 久久中文看片网| 亚洲精品乱码久久久v下载方式| 美女大奶头视频| 99久久精品国产亚洲精品| 人妻夜夜爽99麻豆av| 丰满乱子伦码专区| 国内精品久久久久精免费| 欧美黄色片欧美黄色片| 欧美成人免费av一区二区三区| 欧美日韩黄片免| 国产亚洲精品久久久com| 全区人妻精品视频| 在现免费观看毛片| 亚洲成人久久性| 国产欧美日韩一区二区精品| 精品人妻偷拍中文字幕| bbb黄色大片| 中文字幕av成人在线电影| 亚洲成人久久爱视频| 免费在线观看影片大全网站| 国产美女午夜福利| 亚洲真实伦在线观看| 99热这里只有是精品在线观看 | 欧美高清性xxxxhd video| 男人舔女人下体高潮全视频| 丰满人妻熟妇乱又伦精品不卡| 男女做爰动态图高潮gif福利片| 一卡2卡三卡四卡精品乱码亚洲| 亚洲国产精品999在线| 成年免费大片在线观看| 18禁黄网站禁片免费观看直播| 国产精品久久久久久久久免 | aaaaa片日本免费| 好男人电影高清在线观看| 亚洲成人久久爱视频| 亚洲 国产 在线| 精品无人区乱码1区二区| 大型黄色视频在线免费观看| 非洲黑人性xxxx精品又粗又长| 成人国产一区最新在线观看| 欧美丝袜亚洲另类 | 99国产综合亚洲精品| 亚洲午夜理论影院| 免费观看人在逋| 国产伦一二天堂av在线观看| 亚洲aⅴ乱码一区二区在线播放| 757午夜福利合集在线观看| 久久久久久久精品吃奶| 欧美区成人在线视频| 日本一二三区视频观看| 一区二区三区高清视频在线| 免费在线观看亚洲国产| 精华霜和精华液先用哪个| 一卡2卡三卡四卡精品乱码亚洲| 一区福利在线观看| av黄色大香蕉| 国产av不卡久久| 日日夜夜操网爽| a级毛片免费高清观看在线播放| 99久久99久久久精品蜜桃| 俄罗斯特黄特色一大片| 黄色一级大片看看| 国产成人影院久久av| 精品一区二区三区视频在线观看免费| 精品人妻熟女av久视频| 午夜福利18| 成熟少妇高潮喷水视频| 日韩欧美一区二区三区在线观看| 真实男女啪啪啪动态图| 亚洲熟妇熟女久久| 国产熟女xx| 国产高清有码在线观看视频| 黄色日韩在线| 真人一进一出gif抽搐免费| 在线观看av片永久免费下载| 99久久精品国产亚洲精品| 精品久久久久久,| 极品教师在线视频| 久久久久久久久大av| 成人美女网站在线观看视频| 综合色av麻豆| 成人高潮视频无遮挡免费网站| 女人十人毛片免费观看3o分钟| 黄色日韩在线| 欧美色视频一区免费| 欧美+亚洲+日韩+国产| 最新在线观看一区二区三区| 欧美bdsm另类| 又黄又爽又刺激的免费视频.| 国产国拍精品亚洲av在线观看| 久久久久久久久中文| 观看美女的网站| 免费无遮挡裸体视频| 嫁个100分男人电影在线观看| 免费在线观看影片大全网站| 久久99热6这里只有精品| 精品久久久久久久久亚洲 | 99视频精品全部免费 在线| 欧美色视频一区免费| 国产精品久久久久久亚洲av鲁大| 免费看日本二区| 男人舔女人下体高潮全视频| 亚洲av成人精品一区久久| 国产精品日韩av在线免费观看| 一本综合久久免费| 久久天躁狠狠躁夜夜2o2o| 一区二区三区激情视频| 一进一出抽搐动态| 国内毛片毛片毛片毛片毛片| 可以在线观看毛片的网站| 露出奶头的视频| 亚洲欧美日韩东京热| 97热精品久久久久久| 两个人的视频大全免费| 女人被狂操c到高潮| 国产伦精品一区二区三区视频9| 国产 一区 欧美 日韩| 哪里可以看免费的av片| 久久久久久久久久成人| 亚洲精品久久国产高清桃花| 午夜福利高清视频| 夜夜夜夜夜久久久久| 精品国产三级普通话版| 中文字幕久久专区| av在线天堂中文字幕| 国产白丝娇喘喷水9色精品| 中国美女看黄片| 一区二区三区高清视频在线| 国产精品伦人一区二区| 国产精品,欧美在线| 午夜视频国产福利| 日韩av在线大香蕉| 国产午夜福利久久久久久| 很黄的视频免费| 男人舔女人下体高潮全视频| 国产色婷婷99| 日本免费a在线| 国产老妇女一区| 国产成人aa在线观看| 18禁黄网站禁片免费观看直播| 中文字幕人妻熟人妻熟丝袜美| 女生性感内裤真人,穿戴方法视频| 国产精品99久久久久久久久| 国产成人影院久久av| 大型黄色视频在线免费观看| 成人亚洲精品av一区二区| 男女床上黄色一级片免费看| 亚洲成av人片在线播放无| 夜夜爽天天搞| 国产精品一区二区性色av| 日本黄色片子视频| 色吧在线观看| 欧美成人免费av一区二区三区| 色吧在线观看| 国产亚洲精品av在线| 亚洲精品影视一区二区三区av| 免费在线观看影片大全网站| 国产精品久久久久久久久免 | 内地一区二区视频在线| 有码 亚洲区| 成人高潮视频无遮挡免费网站| 给我免费播放毛片高清在线观看| 亚洲乱码一区二区免费版| 免费黄网站久久成人精品 | 两性午夜刺激爽爽歪歪视频在线观看| 高清毛片免费观看视频网站| 亚洲精品一区av在线观看| www.熟女人妻精品国产| 午夜福利免费观看在线| 免费av不卡在线播放| 国产真实乱freesex| 亚洲精品成人久久久久久| 老师上课跳d突然被开到最大视频 久久午夜综合久久蜜桃 | 欧美国产日韩亚洲一区| 有码 亚洲区| 久久久久国产精品人妻aⅴ院| 午夜精品在线福利| 国产三级中文精品| 国内久久婷婷六月综合欲色啪| 成年女人永久免费观看视频| 看十八女毛片水多多多| 成人性生交大片免费视频hd| 五月伊人婷婷丁香| 精品久久久久久久久久免费视频| 天美传媒精品一区二区| 精品一区二区三区人妻视频| 国产单亲对白刺激| 十八禁网站免费在线| 国产精品99久久久久久久久| 久久午夜亚洲精品久久| 精品久久久久久,| 99热这里只有是精品50| 深爱激情五月婷婷| 全区人妻精品视频| 久久久国产成人精品二区| 美女高潮喷水抽搐中文字幕| 两性午夜刺激爽爽歪歪视频在线观看| 真人一进一出gif抽搐免费| 亚洲成人中文字幕在线播放| 啦啦啦韩国在线观看视频| 69av精品久久久久久| 日韩免费av在线播放| 欧美极品一区二区三区四区| 亚洲国产欧美人成| av天堂中文字幕网| 久久99热6这里只有精品| 日韩欧美三级三区| 两性午夜刺激爽爽歪歪视频在线观看| 国产欧美日韩一区二区精品| 国产一区二区在线观看日韩| 天堂av国产一区二区熟女人妻| 毛片一级片免费看久久久久 | 欧美在线黄色| a级一级毛片免费在线观看| 亚洲国产精品成人综合色| 中文字幕av成人在线电影| 国产精品永久免费网站| 90打野战视频偷拍视频| 熟女电影av网| 蜜桃亚洲精品一区二区三区| www.999成人在线观看| 黄色一级大片看看| 少妇人妻精品综合一区二区 | 亚洲av不卡在线观看| 天堂√8在线中文| 免费无遮挡裸体视频| 中文字幕精品亚洲无线码一区| 国产伦精品一区二区三区视频9| 9191精品国产免费久久| 亚洲最大成人手机在线| 91av网一区二区| av专区在线播放| 国模一区二区三区四区视频| av中文乱码字幕在线| 精品久久久久久久久久免费视频| av视频在线观看入口| 欧美+亚洲+日韩+国产| 国产综合懂色| 我要搜黄色片| 国产高清视频在线播放一区| 在线播放国产精品三级| 日韩欧美免费精品| 日韩欧美 国产精品| 在线播放国产精品三级| 免费观看人在逋| 国产精品综合久久久久久久免费| 免费黄网站久久成人精品 | 变态另类成人亚洲欧美熟女| 久久久久性生活片| 97超视频在线观看视频| 免费看光身美女| 狂野欧美白嫩少妇大欣赏| 色综合亚洲欧美另类图片| 国产免费男女视频| 久久久国产成人免费| 国产免费av片在线观看野外av| 成人永久免费在线观看视频| 午夜精品在线福利| 国产高清激情床上av| 亚洲无线观看免费| 一级作爱视频免费观看| 国产伦精品一区二区三区视频9| 一本综合久久免费| 亚洲人成电影免费在线| 久久久成人免费电影| 国产伦在线观看视频一区| 成人性生交大片免费视频hd| 国产精品一及| 欧美激情久久久久久爽电影|