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

    基于改進(jìn)EMD的多繩摩擦提升機(jī)載荷信息特征提取

    2014-06-07 05:55:04石瑞敏楊兆建
    煤炭學(xué)報(bào) 2014年4期
    關(guān)鍵詞:端點(diǎn)提升機(jī)分量

    石瑞敏,楊兆建

    (太原理工大學(xué)機(jī)械工程學(xué)院,山西太原 030024)

    基于改進(jìn)EMD的多繩摩擦提升機(jī)載荷信息特征提取

    石瑞敏,楊兆建

    (太原理工大學(xué)機(jī)械工程學(xué)院,山西太原 030024)

    為解決多繩摩擦提升機(jī)有效提升載荷的間接識(shí)別問題,提出了一種基于振動(dòng)信號(hào)分析處理的特征提取方法。該方法以主軸裝置軸承處水平振動(dòng)信號(hào)為對(duì)象,利用改進(jìn)的經(jīng)驗(yàn)?zāi)J椒纸?EMD)方法將振動(dòng)信號(hào)分解為若干有效固有模態(tài)函數(shù)(IMFs),改善了經(jīng)典EMD方法模態(tài)混疊與端點(diǎn)效應(yīng)現(xiàn)象,選取各階IMF的能量、方差貢獻(xiàn)率與能量矩作為特征值,探討了特征值與提升機(jī)有效載荷之間的內(nèi)在聯(lián)系。結(jié)果表明,該方法中分析信號(hào)的獲取不改變提升機(jī)主軸系統(tǒng)結(jié)構(gòu),不影響提升機(jī)正常運(yùn)行,易于實(shí)現(xiàn)長(zhǎng)期在線監(jiān)測(cè),便于大量基礎(chǔ)數(shù)據(jù)的采集,為提升載荷定量識(shí)別積累了樣本數(shù)據(jù)。所選取的3個(gè)特征值從大小、權(quán)重及分布的角度較好地反映了提升載荷特征信息。

    改進(jìn)EMD;提升載荷;特征提取;振動(dòng)分析

    礦井提升設(shè)備作為煤礦生產(chǎn)中聯(lián)系井上井下的紐帶,是關(guān)系生產(chǎn)與人員安全的關(guān)鍵設(shè)備之一,一旦出現(xiàn)故障將造成巨大的經(jīng)濟(jì)損失甚至人員的傷亡。在諸類安全事故中,由于過裝或未卸凈引起的過卷、滑動(dòng)甚至墜罐,以及由于罐道劣變等引起的卡繩、松繩甚至斷繩等安全事故時(shí)有發(fā)生,而這些故障都可以直接或間接的由提升載荷反映出來。目前對(duì)提升載荷的動(dòng)態(tài)檢測(cè)方法主要有壓輪法與直接測(cè)量法,前者通過將壓輪推向工作過程中的鋼絲繩,測(cè)得提升鋼絲繩對(duì)壓輪作用的壓力從而計(jì)算得到提升載荷,在調(diào)整鋼絲繩張力差方面得到了較好的應(yīng)用效果[1]。后者將傳感器安裝于提升鋼絲繩與提升容器間[2]或安裝于提升容器器件間[3]直接測(cè)得提升載荷。兩類方法均需要對(duì)提升系統(tǒng)結(jié)構(gòu)進(jìn)行改造,對(duì)傳感器與監(jiān)測(cè)系統(tǒng)安全可靠性要求較高。

    振動(dòng)信號(hào)是旋轉(zhuǎn)機(jī)械故障特征信息的載體[4],筆者在對(duì)提升機(jī)主軸兩側(cè)軸承進(jìn)行監(jiān)測(cè)的過程中發(fā)現(xiàn),提升載荷作用于鋼絲繩,在提升過程中引發(fā)的振動(dòng)將傳導(dǎo)至滾筒乃至主軸,使得軸承處振動(dòng)信號(hào)包含著提升載荷的重要信息。與提升設(shè)備其他系統(tǒng)相比,提升載荷的特征信號(hào)比較微弱,信噪比低,具有很強(qiáng)的非線性、非平穩(wěn)及噪聲特性,傳統(tǒng)的信號(hào)處理方法在處理這類信號(hào)方面存在缺陷。如短時(shí)傅里葉變換對(duì)不同時(shí)段變化了的信號(hào)只能加相同的窗,不適應(yīng)信號(hào)頻率不斷變化的不同要求;小波分析克服了傅里葉變化窗口不可調(diào)的缺點(diǎn),提高了分辨率,但對(duì)小波窗內(nèi)的信號(hào)平穩(wěn)性提出要求,因?yàn)闆]有擺脫傅里葉變換的局限,并且存在信號(hào)能量泄露和非自適應(yīng)性的缺點(diǎn)。經(jīng)驗(yàn)?zāi)J椒纸?empirical mode decomposition, EMD)是1998年由黃鍔等提出的可對(duì)非線性非平穩(wěn)信號(hào)進(jìn)行有效分析的信號(hào)處理方法[5]。由于其良好的局部特性與自適應(yīng)性,EMD自提出以來被應(yīng)用于地震信號(hào)[6-7]、語音識(shí)別[8]、故障診斷[9-11]等各個(gè)領(lǐng)域,并取得良好的效果。但是EMD方法仍舊存在一些不足,如模態(tài)混疊、端點(diǎn)效應(yīng)、包絡(luò)線選擇及停止準(zhǔn)則等[12-14]。本文旨在對(duì)EMD方法存在的問題進(jìn)行改進(jìn),并在此基礎(chǔ)上提出一種多繩摩擦提升機(jī)載荷信息特征提取方法的新思路,為進(jìn)一步實(shí)現(xiàn)載荷定量識(shí)別提供理論依據(jù)。

    1 改進(jìn)的EMD方法

    1.1 EMD理論及存在問題

    EMD方法是一種適用于非線性、非平穩(wěn)性信號(hào)的時(shí)頻分析方法。此方法根據(jù)被分析信號(hào)本身特點(diǎn)構(gòu)造基函數(shù),自適應(yīng)的選擇頻帶,從局部時(shí)間尺度入手,確定信號(hào)在不同頻段的分辨率,將復(fù)雜信號(hào)分解為具有不同尺度特征的固有模態(tài)函數(shù)(intrinsic mode functions,IMF),因此,該方法具有良好的局部適應(yīng)性和時(shí)頻分辨率,較好地克服了小波變換存在的問題。EMD方法實(shí)質(zhì)上可看作截止頻率與帶寬隨信號(hào)變化的帶通濾波器,從數(shù)據(jù)的特征時(shí)間尺度入手對(duì)信號(hào)進(jìn)行分解,原始信號(hào)x(t)經(jīng)EMD分解后,得到有限數(shù)目IMF和一個(gè)剩余分量,即

    其中,n為分解后IMF數(shù)量;ci(t)為分解得到的n組IMF分量,各分量均滿足平穩(wěn)性條件;r(t)為分解后的剩余分量。EMD方法因其在處理非線性非平穩(wěn)性信號(hào)時(shí)表現(xiàn)出的良好特性,在實(shí)踐中有大量應(yīng)用,但是理論本身仍有以下需要完善的不足之處。

    (1)在EMD分解過程中存在模態(tài)混疊現(xiàn)象,即同一頻率信號(hào)被分解在不同的IMF中,或一個(gè)IMF中包含多個(gè)不同頻率的信號(hào)。當(dāng)原始信號(hào)中存在間隙性成分或脈沖成分時(shí),模態(tài)混疊現(xiàn)象更為明顯。

    (2)EMD分解過程中采用極值點(diǎn)構(gòu)造包絡(luò)曲線,而被分析信號(hào)端點(diǎn)處往往非極值點(diǎn),這就導(dǎo)致每個(gè)分解過程中都會(huì)在端點(diǎn)處產(chǎn)生擬合誤差,隨著分解進(jìn)行,誤差影響會(huì)逐步擴(kuò)大。尤其是分解后期的低頻IMF分量,端點(diǎn)效應(yīng)現(xiàn)象更為明顯,誤差的累積和傳播會(huì)導(dǎo)致分解結(jié)果的失真,產(chǎn)生虛假的IMF分量。

    1.2 改進(jìn)的EMD

    針對(duì)以上提出的模態(tài)混疊、端點(diǎn)效應(yīng)及虛假IMF分量的缺陷,提出采用集合經(jīng)驗(yàn)?zāi)J椒纸?ensemble empirical mode decomposition,EEMD)方法改善模態(tài)混疊現(xiàn)象,利用鏡像邊界延拓與灰色預(yù)測(cè)模型延拓方法相結(jié)合來抑制端點(diǎn)效應(yīng),最后采用相關(guān)系數(shù)法剔除假IMF分量。

    (1)EEMD算法是Zhaohua Wu和Huang提出的一種噪聲輔助的數(shù)據(jù)分析方法[15],可以在不影響EMD自適應(yīng)性的前提下減少模態(tài)混疊對(duì)EMD分解結(jié)果的影響。EEMD分解過程為:將高斯白噪聲與原始信號(hào)疊加后進(jìn)行多次EMD分解,由于高斯白噪聲具有頻率均勻分布的統(tǒng)計(jì)特性,使得不同尺度的信號(hào)成分自動(dòng)分布到合適的參考尺度上,從而降低各IMF分量的模態(tài)混疊程度;并且由于零均值噪聲的特性,經(jīng)過足夠多次平均運(yùn)算處理后所加入的高斯白噪聲將相互抵消,原始信號(hào)得到還原。

    在應(yīng)用EEMD方法時(shí),有2個(gè)重要參數(shù)需要考慮:首先,對(duì)于所加入的高斯白噪聲的幅值選取十分重要,如果幅值過小,可能不足以引起信號(hào)局部極值點(diǎn)的變化,使得白噪聲的加入對(duì)EMD無作用;如果幅值過大,則會(huì)產(chǎn)生過多虛假IMF分量,使分解誤差增大,甚至淹沒原始信號(hào)特征使分解失真。其次,隨著分解次數(shù)的增加所加入白噪聲的幅值會(huì)隨之降低,理論上,分解次數(shù)越大則分解誤差越小,但是過大的分解次數(shù)代價(jià)是計(jì)算效率的降低,故應(yīng)選取適當(dāng)值。Huang[15]建議所加入白噪聲的幅值可以由原始信號(hào)的標(biāo)準(zhǔn)差乘以一個(gè)系數(shù)(一般取0.2)來確定,信號(hào)中高頻成分多時(shí)系數(shù)適當(dāng)減小,反之則適當(dāng)增大。分解次數(shù)可通過設(shè)置分解誤差來確定。

    (2)抑制端點(diǎn)效應(yīng)多采用在原始信號(hào)兩端通過一定方法增加信號(hào)長(zhǎng)度,經(jīng)HHT變換后再從結(jié)果中除去延長(zhǎng)部分方法。目前所提出波形延拓、神經(jīng)網(wǎng)絡(luò)延拓、支持向量回歸預(yù)測(cè)及邊界極值延拓等方法,在一定程度上消除了端點(diǎn)效應(yīng)的影響,但仍存在局限性。本文提出利用灰色預(yù)測(cè)模型與鏡像延拓相結(jié)合的方法來抑制端點(diǎn)效應(yīng)。根據(jù)被分析信號(hào)的特性,考慮到其發(fā)展過程與各影響因素是一種灰色不確定性關(guān)系,選取灰色動(dòng)態(tài)預(yù)測(cè)模型GM(1,1)對(duì)信號(hào)進(jìn)行延拓,此方法具有要求樣本數(shù)據(jù)少、原理簡(jiǎn)單、計(jì)算量小、可檢驗(yàn)等優(yōu)點(diǎn)。具體延拓過程如下:

    ①以信號(hào)左端延拓為例,取靠近左端點(diǎn)處的5個(gè)極大(小)值組成原始序列,建立原始信號(hào)x(t)的灰色微分方程:

    其中,a為發(fā)展灰數(shù);b為內(nèi)生控制灰數(shù)。二者由已知數(shù)據(jù)及最小二乘法確定,在此不贅述。求出a,b后可得預(yù)測(cè)模型:

    ②先用原始序列建立的GM(1,1)模型預(yù)測(cè)得到一個(gè)值,然后把這個(gè)值加入到原始序列中,同時(shí)剔除一個(gè)最遠(yuǎn)端數(shù)據(jù),以新的原始序列重復(fù)步驟①,再預(yù)測(cè)下一個(gè)數(shù)值,循環(huán)進(jìn)行直至達(dá)到要求的預(yù)測(cè)精度。為滿足信號(hào)延拓要求,需要對(duì)數(shù)據(jù)極大值和極小值分別預(yù)測(cè)兩個(gè)值。

    ③對(duì)信號(hào)右端進(jìn)行相同處理,最后將延拓得到的端點(diǎn)值與原信號(hào)連接起來形成新的信號(hào)序列。在此基礎(chǔ)上應(yīng)用鏡像延拓法,以預(yù)測(cè)得到的各端點(diǎn)作為鏡像點(diǎn)進(jìn)行延拓,從而得到較理想的延拓?cái)?shù)據(jù)。

    (3)EMD分解后得到的各IMF分量中,真實(shí)IMF分量與原始信號(hào)在理論上滿足正交性,具有較大的相關(guān)系數(shù),而虛假IMF分量與原始信號(hào)相關(guān)系數(shù)較小。因此,可通過計(jì)算IMF分量與原始信號(hào)的相關(guān)系數(shù)來識(shí)別并剔除虛假IMF分量。

    2 改進(jìn)EMD方法的信號(hào)消噪分頻處理

    運(yùn)用改進(jìn)的HHT方法對(duì)原始信號(hào)進(jìn)行處理流程如圖1所示,其步驟包括:

    圖1 改進(jìn)的EMD處理流程Fig.1 Flowchart of improved EMD

    (1)將原始信號(hào)x(t)加入一個(gè)正態(tài)分布的白噪聲序列nj(t),使得x~(t)=x(t)+nj(t)。

    (2)將得到的信號(hào)x~(t)進(jìn)行EMD分解,過程中在求取信號(hào)極值點(diǎn)后,采用鏡像與灰色預(yù)測(cè)模型相結(jié)合的方式對(duì)信號(hào)進(jìn)行延拓,最終得到各IMF分量ci(t),i=1,2,3,…,n。

    (3)重復(fù)步驟(1)2N次,每次加入新的隨機(jī)正態(tài)分布的白噪聲序列nj(t)。

    (4)將N次分解得到的固有模態(tài)分量IMF做集成平均處理,當(dāng)N足夠大時(shí),可使添加白噪聲的影響趨近于0,進(jìn)而得到IMF的最終結(jié)果。

    (5)計(jì)算各IMF分量與原始信號(hào)的相關(guān)系數(shù),以此為依據(jù)篩選有效IMF集。

    (6)將有效IMF集重構(gòu)信號(hào),進(jìn)行Hilbert變換。

    為驗(yàn)證此方法的有效性,將3個(gè)頻率分別為3, 11,16 Hz,幅值分別為1,0.4,0.9的正弦信號(hào)疊加,如下式:x(t)=sin(2π×3t)+0.4sin(2π×11t)+ 0.9sin(2π×16t),采樣頻率為500 Hz,采樣長(zhǎng)度1 000點(diǎn),時(shí)域信號(hào)如圖2所示。

    圖2 仿真信號(hào)時(shí)域Fig.2 Time-domain chart of simulation signal

    對(duì)數(shù)據(jù)分別采用原始與改進(jìn)EMD兩種方法進(jìn)行分析。圖3為原始EMD分解結(jié)果與信號(hào)的Hilbert時(shí)頻譜圖,從圖3(a),(b)可看出,IMF1由兩種頻率的正弦信號(hào)疊加,出現(xiàn)了明顯的模態(tài)混疊現(xiàn)象,反映到圖3(c)上,11與16 Hz譜線對(duì)應(yīng)位置出現(xiàn)較大變形,未正確區(qū)分。另外從Hilbert時(shí)頻譜中可看出,信號(hào)兩端出現(xiàn)了明顯的頻率發(fā)散,并向信號(hào)內(nèi)部傳播,靠近邊界處尤為嚴(yán)重,只有信號(hào)中間部分可以反映信號(hào)的特征。采用本文改進(jìn)EMD方法對(duì)信號(hào)重新分解結(jié)果如圖4所示,從圖4(a),(b)EMD分解結(jié)果及頻譜可以看出,相近頻率11與16 Hz成功分解到兩個(gè)IMF分量中,模態(tài)混疊現(xiàn)象得以解決,從圖4(c)亦可看出,頻率發(fā)散現(xiàn)象趨勢(shì)變緩,端點(diǎn)效應(yīng)得到抑制。

    圖3 原始EMD分解結(jié)果Fig.3 Decompose results of original EMD

    圖4 改進(jìn)EMD分解結(jié)果Fig.4 Decompose results of improved EMD

    3 基于IMF能量信息的提升載荷信息特征提取

    反映提升機(jī)載荷信息的特征參數(shù)有電流、振動(dòng)或直接測(cè)量鋼絲繩張力等。當(dāng)提升機(jī)上提和下放載荷時(shí),提升容器的振動(dòng)會(huì)通過鋼絲繩作用于提升機(jī)從而反應(yīng)在包括軸承與減速器的主軸系統(tǒng)上,由于載荷變化會(huì)引起軸承水平方向振動(dòng)信號(hào)的波動(dòng),故本文拾取主軸裝置兩側(cè)軸承水平方向振動(dòng)信號(hào),采用改進(jìn)的EMD方法對(duì)振動(dòng)信號(hào)中包含的載荷特征信息進(jìn)行提取。經(jīng)改進(jìn)EMD處理后得到的有效IMF集中,各階IMF分量所含的頻率成分不同,階數(shù)越低所含信號(hào)的高頻成分則越多。各階IMF分量的能量以及在整個(gè)頻帶上的能量分布能夠刻畫出信號(hào)的整體特征,因此,本文提出以IMF分量的能量變化作為參考依據(jù)提取載荷信息特征。

    圖5為某礦副立井型號(hào)為JKM 2.8×6(I)A的多繩摩擦提升機(jī)主軸系統(tǒng),傳感器布置在軸承座頂端,采集每個(gè)提升過程中水平方向(以圖中X方向?yàn)檎?振動(dòng)信號(hào),采樣頻率500 Hz,截取勻速提升時(shí)6 s數(shù)據(jù)作為待分析信號(hào)。副立井主要負(fù)責(zé)運(yùn)送人員、設(shè)備以及提升矸石等,故載荷情況較為復(fù)雜,每次提升載荷情況差別較大,本文以1號(hào)罐提升工況選取5組數(shù)據(jù)作為示例,載荷數(shù)據(jù)見表1。

    圖5 多繩摩擦提升機(jī)主軸裝置傳感器布置Fig.5 Sensors placement on spindle device of multi-rope hoist

    表1 示例載荷數(shù)據(jù)Table 1 Load data of sampleskN

    首先采用改進(jìn)的EMD方法對(duì)被分析信號(hào)進(jìn)行消噪分頻處理,得到5階有效IMF集,如圖6所示。

    圖6 實(shí)測(cè)信號(hào)的改進(jìn)EMD分解Fig.6 Improved EMD of actual signal

    針對(duì)被測(cè)振動(dòng)信號(hào)的非平穩(wěn)性特點(diǎn),當(dāng)載荷發(fā)生變化時(shí),其振動(dòng)信號(hào)在各個(gè)頻段的能量必將發(fā)生變化,而隨著工況不同,經(jīng)改進(jìn)的EMD分解后的各有效IMF分量包含了不同頻段的信息,為了更好地提取隱藏在信號(hào)中的載荷信息,選取IMF分量能量值,方差貢獻(xiàn)率及能量矩從能量大小、權(quán)重和分布的角度對(duì)信號(hào)本質(zhì)特征進(jìn)行分析。其定義分別為

    圖7為各階IMF分量特征值,各特征值隨載荷的變化。從圖7(a)可以看到,隨著提升載荷總重的增加,第1階IMF能量值呈上升趨勢(shì),但當(dāng)提升載荷總重相近時(shí),載荷總重大但載荷差小的工況第1階IMF能量值反而小,需要從第2階及以后的IMF能量上反映真實(shí)載荷情況;圖7(b)為各階IMF分量的方差貢獻(xiàn)率,體現(xiàn)了各階能量在總能量中的權(quán)重,可以看到能量主要集中于前兩階IMF分量中;還選取了能量矩作為表征載荷信息的特征值,如圖7(c)所示,能量矩指標(biāo)考慮到了IMF能量隨時(shí)間t的分布情況,可以體現(xiàn)能量的分布特征。從上述分析可以得到,相對(duì)于不能完整快捷反映載荷信息的單一的能量指標(biāo),將可以體現(xiàn)能量的大小、權(quán)重及分布的指標(biāo)結(jié)合起來,可以更好地揭示IMF集能量特征以及與載荷的對(duì)應(yīng)關(guān)系,有利于后續(xù)的載荷識(shí)別與故障診斷。

    圖7 各階IMF分量特征值Fig.7 Characteristic parameters of IMFs

    4 結(jié) 論

    (1)多繩摩擦提升機(jī)主軸裝置軸承處水平振動(dòng)信號(hào)中蘊(yùn)含了提升載荷的一些特征信息。

    (2)提出了一種改進(jìn)的經(jīng)驗(yàn)?zāi)B(tài)分解方法,仿真結(jié)果表明此方法改善了原始EMD方法存在的模態(tài)混疊、端點(diǎn)效應(yīng)以及虛假IMF分量問題,快速準(zhǔn)確的得到了有效IMF集,可用于其他類型信號(hào)的消噪分頻處理。

    (3)利用固有模態(tài)函數(shù)的能量值、方差貢獻(xiàn)率與能量矩從能量大小、權(quán)重及分布3方面相結(jié)合作為提升載荷的主要特征,可以更好地反映提升載荷與振動(dòng)信號(hào)的對(duì)應(yīng)關(guān)系。

    (4)特征提取是提升載荷識(shí)別和診斷中的重要一環(huán),本文提出的方法可以在不改變提升主軸系統(tǒng)結(jié)構(gòu),不影響提升機(jī)正常運(yùn)行的前提下有效提取多繩摩擦提升機(jī)提升載荷特征,可作為進(jìn)一步建立提升載荷識(shí)別模型的依據(jù),為提升載荷定量識(shí)別與監(jiān)測(cè)提供了新思路。

    [1] 楊兆建,賈長(zhǎng)喜,王勤賢,等.多繩提升機(jī)鋼絲繩張力動(dòng)態(tài)測(cè)試調(diào)整研究[J].煤炭學(xué)報(bào),1995,20(5):519-524.

    Yang Zhaojian,Jia Changxi,Wang Qinxian,et al.Dynamic test and adjustment of wire-rope tension of multiple rope friction winder[J].Journal of China Coal Society,1995,20(5):519-524.

    [2] 騰孝來,肖興明,胡 明,等.多繩摩擦提升鋼絲繩張力動(dòng)態(tài)無線監(jiān)測(cè)[J].煤礦機(jī)械,2008,29(5):207-208.

    Teng Xiaolai,Xiao Xingming,Hu Ming,et al.Dynamic wireless tension monitoring of wire-rope of multiple friction winder[J].Coal Mine Machinery,2008,29(5):207-208.

    [3] 邵海燕.多繩提升機(jī)載荷動(dòng)態(tài)監(jiān)測(cè)與提升安全研究[D].青島:山東科技大學(xué),2003.

    Shao Haiyan.Study of the real-time monitoring and security of multiwire hoisting load[D].Qingdao:Shandong University of Science and Technology,2003.

    [4] 鞠萍華,秦樹人,秦 毅,等.多分辨EMD方法與頻域平均在齒輪早期故障診斷中的研究[J].振動(dòng)與沖擊,2009,28(5):97-101.

    Ju Pinghua,Qin Shuren,Qin Yi,et al.Research on earlier fault diagnosis of gear by method of multi-resolution empirical mode decomposition and frequency domain averaging[J].Journal of Vibration and Shock,2009,28(5):97-101.

    [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 of London, 1998,454(A):903-995.

    [6] Huang N E,Chern C C,Huang K,et al.A new spectral representation of earthquake data:Hilbert spectral analysis of Station TCU129, Chi-Chi,Taiwan,21 September 1999[J].Bulletin of the Seismological Society of America,2001,91(5):1310-1338.

    [7] 石春香,羅奇峰.時(shí)程信號(hào)的Hilbert-Huang變換與小波分析[J].地震學(xué)報(bào),2003,25(4):398-405.

    Shi Chunxiang,Luo Qifeng.Hilbert-Huang transform and wavelet analysis of time history signal[J].Acta Seismologica Sinica,2003, 25(4):398-405.

    [8] Aicha Bouzid,Noureddine Ellouze.Empirical mode decomposition of voiced speech signal[A].First International Symposium on Control, Communications and Signal Processing[C].2004:603-606.

    [9] 程軍圣,于德介,楊 宇,等.EMD方法在轉(zhuǎn)子局部碰摩故障診斷中的應(yīng)用[J].振動(dòng)、測(cè)試與診斷,2006,26(1):24-27.

    Cheng Junsheng,Yu Dejie,Yang Yu,et al.Application of EMD to local rub-impact fault diagnosis in rotor systems[J].Journal of Vibration Measurement&Diagnosis,2006,26(1):24-27.

    [10] Bassiunya A M,Lib Xiaoli.Flute breakage detection during end milling using Hilbert-Huang transform and smoothed nonlinear energy operator[J].International Journal of Machine Tools and Manufac-ture,2007,47(6):1011-1020.

    [11] 時(shí)培明,丁雪娟,李 庚,等.一種EMD改進(jìn)方法及其在旋轉(zhuǎn)機(jī)械故障診斷中的應(yīng)用[J].振動(dòng)與沖擊,2013,32(4):185-190.

    Shi Peiming,Ding Xuejuan,Li Geng,et al.An improved method of EMD and its application in rotating machinery fault diagnosis[J].Journal of Vibration and Shock,2013,32(4):185-190.

    [12] 李方熙,陳桂明,劉希亮,等.希爾伯特-黃變換端點(diǎn)效應(yīng)的自適應(yīng)端點(diǎn)相位正弦延拓方法[J].上海交通大學(xué)學(xué)報(bào),2013, 47(4):594-601.

    Li Fangxi,Chen Guiming,Liu Xiliang,et al.Processing method for Hilbert-Huang transform end effects self-adaptive endpoint-phase sinusoidal extension[J].Journal of Shanghai Jiaotong University, 2013,47(4):594-601.

    [13] 鐘佑明,金 濤,秦樹人.希爾伯特一黃變換中的一種新包絡(luò)線算法[J].數(shù)據(jù)采集與處理,2005,20(1):13-18.

    Zhong Youming,Jin Tao,Qin Shuren.New envelope algorithm for Hilbert-Huang transform[J].Journal of Data Acquisition&Processing,2005,20(1):13-18.

    [14] 王增才,王樹梁,任鍇勝,等.基于EEMD的提升機(jī)天輪軸承故障診斷方法[J].煤炭學(xué)報(bào),2012,37(4):659-694.

    Wang Zengcai,Wang Shuliang,Ren Kaisheng,et al.Research on the method of hoist head sheave bearing fault diagnosis based on EEMD[J].Journal of China Coal Society,2012,37(4):659-694.

    [15] Wu Zhaohua,Norden Eh Huang.Ensemble empirical mode decomposition:A noise assisted data analysis method[J].Advances in Adaptive Data Analysis,2008,1(1):1-41.

    Feature extraction for hoisting load of multiple rope friction hoist based on improved EMD

    SHI Rui-min,YANG Zhao-jian

    (School of Mechanical Engineering,Taiyuan University of Technology,Taiyuan 030024,China)

    In order to indirectly identity the payload of multi-rope friction hoist,a feature extraction method was proposed based on analysis of horizontal vibration signal which was collected from bearing caps of spindle device.Firstly, improved empirical mode decomposition(EMD)that can restrain the mode mixing phenomenon and reduce the boundary effect was used to decompose the signal into several intrinsic mode functions(IMFs),then the energy,contribution rate of variance and energy moment of IMFs were calculated as characteristic parameters.The intrinsic relationship between the characteristic parameters and the payload of hoist were discussed.The method of signal acquisition didn’t change the structure of spindle system and had no effect on normal operation of hoist.A long-term online monitoring was apt to accomplish,the basic data was easy to obtain,furthermore,a large amounts of sample datum were accumulated for quantitative identification of actual load.It is shown that the characteristic parameters proposed could sensitively reflect the load characters in terms of magnitude,weight and distribution.

    improved EMD;hoisting load;feature extraction;vibration analysis

    TD53

    A

    0253-9993(2014)04-0782-07

    石瑞敏,楊兆建.基于改進(jìn)EMD的多繩摩擦提升機(jī)載荷信息特征提取[J].煤炭學(xué)報(bào),2014,39(4):782-788.

    10.13225/j.cnki.jccs.2013.1502

    Shi Ruimin,Yang Zhaojian.Feature extraction for hoisting load of multiple rope friction hoist based on improved EMD[J].Journal of China Coal Society,2014,39(4):782-788.doi:10.13225/j.cnki.jccs.2013.1502

    2013-10-18 責(zé)任編輯:許書閣

    山西省科技重大專項(xiàng)資助項(xiàng)目(20111101040)

    石瑞敏(1983—),女,山西文水人,博士研究生。Tel:0351-6010414,E-mail:srm0018@link.tyut.edu.cn。通訊作者:楊兆建(1955—),男,河北高陽人,教授,博士生導(dǎo)師。Tel:0351-6010414,E-mail:yangzhaojian@tyut.edu.cn

    猜你喜歡
    端點(diǎn)提升機(jī)分量
    干熄焦提升機(jī)控制系統(tǒng)的改造
    山東冶金(2022年4期)2022-09-14 09:00:00
    非特征端點(diǎn)條件下PM函數(shù)的迭代根
    PLC技術(shù)在煤礦提升機(jī)控制系統(tǒng)的應(yīng)用
    帽子的分量
    一物千斤
    智族GQ(2019年9期)2019-10-28 08:16:21
    不等式求解過程中端點(diǎn)的確定
    論《哈姆雷特》中良心的分量
    分量
    參數(shù)型Marcinkiewicz積分算子及其交換子的加權(quán)端點(diǎn)估計(jì)
    基丁能雖匹配延拓法LMD端點(diǎn)效應(yīng)處理
    av在线天堂中文字幕| av福利片在线观看| 精品日产1卡2卡| 日本黄色片子视频| 一级av片app| 99久国产av精品| 中文字幕精品亚洲无线码一区| 亚洲欧美清纯卡通| 嫩草影院入口| 亚洲黑人精品在线| 精品久久久久久成人av| 国产黄片美女视频| 美女免费视频网站| 校园人妻丝袜中文字幕| 精品日产1卡2卡| 国产精品99久久久久久久久| 亚洲精品成人久久久久久| 听说在线观看完整版免费高清| 成人av在线播放网站| 精品人妻一区二区三区麻豆 | 啦啦啦韩国在线观看视频| or卡值多少钱| 日本黄色视频三级网站网址| 欧美xxxx黑人xx丫x性爽| 久久久久久久午夜电影| 亚洲av.av天堂| 亚洲成人免费电影在线观看| 精品久久久久久久久亚洲 | 在线观看舔阴道视频| 日本免费一区二区三区高清不卡| 国产av不卡久久| 午夜精品在线福利| av专区在线播放| 午夜福利在线观看免费完整高清在 | 中文字幕人妻熟人妻熟丝袜美| 成人国产综合亚洲| 色综合亚洲欧美另类图片| 久久久久久久精品吃奶| 直男gayav资源| 白带黄色成豆腐渣| 国产免费av片在线观看野外av| 99国产精品一区二区蜜桃av| 亚洲国产精品久久男人天堂| 国语自产精品视频在线第100页| 亚洲精品成人久久久久久| 男女之事视频高清在线观看| 一卡2卡三卡四卡精品乱码亚洲| 内射极品少妇av片p| 免费看a级黄色片| or卡值多少钱| 午夜福利成人在线免费观看| 一级黄色大片毛片| 在线观看一区二区三区| 日韩av在线大香蕉| 日本三级黄在线观看| 婷婷亚洲欧美| 我的女老师完整版在线观看| 99在线视频只有这里精品首页| 亚洲在线自拍视频| 亚洲 国产 在线| 成年女人永久免费观看视频| 三级毛片av免费| 天堂av国产一区二区熟女人妻| 夜夜看夜夜爽夜夜摸| 亚洲人成网站在线播放欧美日韩| 日韩欧美国产一区二区入口| 国产黄色小视频在线观看| 国产69精品久久久久777片| 亚洲精华国产精华液的使用体验 | 日本-黄色视频高清免费观看| 91麻豆精品激情在线观看国产| a在线观看视频网站| 中文字幕人妻熟人妻熟丝袜美| 中文亚洲av片在线观看爽| 色吧在线观看| 国产精品美女特级片免费视频播放器| 精品国内亚洲2022精品成人| 最好的美女福利视频网| 色播亚洲综合网| 蜜桃久久精品国产亚洲av| 伦理电影大哥的女人| 久久精品夜夜夜夜夜久久蜜豆| 99久久无色码亚洲精品果冻| 成人综合一区亚洲| av在线观看视频网站免费| 我要看日韩黄色一级片| 内地一区二区视频在线| 黄色日韩在线| 免费av毛片视频| 亚洲va在线va天堂va国产| 男女那种视频在线观看| 久久欧美精品欧美久久欧美| 久久6这里有精品| 成人亚洲精品av一区二区| 99久久九九国产精品国产免费| 最新中文字幕久久久久| 日本 欧美在线| 天堂影院成人在线观看| 国内少妇人妻偷人精品xxx网站| 成年女人永久免费观看视频| 午夜福利成人在线免费观看| 欧美最新免费一区二区三区| 欧美日韩综合久久久久久 | 国产乱人伦免费视频| 两性午夜刺激爽爽歪歪视频在线观看| 亚洲成av人片在线播放无| 日韩大尺度精品在线看网址| 蜜桃亚洲精品一区二区三区| 成人无遮挡网站| 哪里可以看免费的av片| 国产精品av视频在线免费观看| 国产一区二区三区av在线 | 日韩av在线大香蕉| 中文字幕熟女人妻在线| 日本色播在线视频| 国产高清激情床上av| 欧美激情在线99| 麻豆国产av国片精品| 日本 欧美在线| 国产精品自产拍在线观看55亚洲| 亚洲精品国产成人久久av| 狂野欧美白嫩少妇大欣赏| 亚洲精华国产精华液的使用体验 | 国产极品精品免费视频能看的| 波多野结衣高清作品| 欧美xxxx性猛交bbbb| 中国美白少妇内射xxxbb| 国产一区二区三区av在线 | 精品午夜福利在线看| 天堂√8在线中文| 熟妇人妻久久中文字幕3abv| 欧美区成人在线视频| 又粗又爽又猛毛片免费看| 亚洲精品色激情综合| 91在线精品国自产拍蜜月| 精品无人区乱码1区二区| 舔av片在线| 午夜老司机福利剧场| 又爽又黄a免费视频| 婷婷精品国产亚洲av| avwww免费| 一个人观看的视频www高清免费观看| 精品午夜福利视频在线观看一区| 国产淫片久久久久久久久| 国产黄色小视频在线观看| 极品教师在线免费播放| 色哟哟·www| 国产人妻一区二区三区在| 国产又黄又爽又无遮挡在线| 久久精品国产亚洲av涩爱 | 日韩av在线大香蕉| 一级av片app| 老师上课跳d突然被开到最大视频| 两性午夜刺激爽爽歪歪视频在线观看| 成人av在线播放网站| 很黄的视频免费| 久久久久久久久久久丰满 | 亚洲人成网站在线播放欧美日韩| 国产精品无大码| 欧美精品国产亚洲| 少妇裸体淫交视频免费看高清| 亚洲国产色片| 国产高清三级在线| 国产精品1区2区在线观看.| 国产av在哪里看| 亚洲一级一片aⅴ在线观看| 国产一区二区激情短视频| 99热这里只有是精品50| 亚洲国产欧美人成| 一卡2卡三卡四卡精品乱码亚洲| av.在线天堂| 国产蜜桃级精品一区二区三区| 国内精品宾馆在线| 一个人看视频在线观看www免费| 日日撸夜夜添| 成人鲁丝片一二三区免费| .国产精品久久| 久久精品影院6| 男人狂女人下面高潮的视频| 又紧又爽又黄一区二区| 亚洲av一区综合| 国产aⅴ精品一区二区三区波| 成人国产一区最新在线观看| 亚洲五月天丁香| 日韩大尺度精品在线看网址| 国产精品av视频在线免费观看| 国产真实乱freesex| 内射极品少妇av片p| 国产精品国产三级国产av玫瑰| 国产麻豆成人av免费视频| 色综合色国产| 亚洲av成人av| 免费黄网站久久成人精品| 日本黄大片高清| 美女xxoo啪啪120秒动态图| 亚洲人成网站在线播| 国产精品乱码一区二三区的特点| 精品国内亚洲2022精品成人| 久久亚洲真实| 亚洲成a人片在线一区二区| 午夜日韩欧美国产| 永久网站在线| 日韩欧美三级三区| 在线观看午夜福利视频| 国产激情偷乱视频一区二区| 丰满人妻一区二区三区视频av| 国产亚洲欧美98| 国产精品久久久久久久电影| 两个人视频免费观看高清| 国产在线精品亚洲第一网站| 在线观看一区二区三区| 啦啦啦啦在线视频资源| 国产女主播在线喷水免费视频网站 | 久久这里只有精品中国| 日本撒尿小便嘘嘘汇集6| 国产视频一区二区在线看| 欧美性猛交╳xxx乱大交人| 久久精品国产鲁丝片午夜精品 | 女的被弄到高潮叫床怎么办 | 性欧美人与动物交配| 熟女电影av网| 日本a在线网址| 精品久久久噜噜| 3wmmmm亚洲av在线观看| 桃色一区二区三区在线观看| 两个人的视频大全免费| 在线观看午夜福利视频| 欧美不卡视频在线免费观看| 国产精品自产拍在线观看55亚洲| 观看免费一级毛片| 男人舔女人下体高潮全视频| 久久久国产成人精品二区| 动漫黄色视频在线观看| 欧美日韩乱码在线| 亚洲精华国产精华液的使用体验 | 一区福利在线观看| 99热这里只有是精品在线观看| 精品不卡国产一区二区三区| 欧美成人性av电影在线观看| 免费av观看视频| 国产精品一及| 亚洲人成网站在线播| 午夜福利高清视频| 久久国内精品自在自线图片| 欧美日韩综合久久久久久 | 国产激情偷乱视频一区二区| 五月玫瑰六月丁香| 99久久成人亚洲精品观看| 人人妻,人人澡人人爽秒播| www.www免费av| 国产精品亚洲美女久久久| 亚洲自偷自拍三级| 国产又黄又爽又无遮挡在线| 黄色女人牲交| 99热这里只有精品一区| 日本色播在线视频| 欧美丝袜亚洲另类 | 国产日本99.免费观看| 日韩中文字幕欧美一区二区| 欧美成人一区二区免费高清观看| 成人高潮视频无遮挡免费网站| 亚洲内射少妇av| 亚洲天堂国产精品一区在线| 男女做爰动态图高潮gif福利片| 亚洲人成伊人成综合网2020| 亚洲中文字幕一区二区三区有码在线看| 亚洲七黄色美女视频| 中文字幕高清在线视频| 国产一区二区激情短视频| 午夜久久久久精精品| 精品人妻视频免费看| 国产精品一区www在线观看 | 啦啦啦啦在线视频资源| 精品人妻1区二区| 久久国内精品自在自线图片| 老熟妇乱子伦视频在线观看| 校园人妻丝袜中文字幕| 欧美人与善性xxx| 亚洲内射少妇av| 一本精品99久久精品77| avwww免费| 成人鲁丝片一二三区免费| 22中文网久久字幕| 午夜爱爱视频在线播放| 三级国产精品欧美在线观看| 给我免费播放毛片高清在线观看| 在线a可以看的网站| 天美传媒精品一区二区| 国产成人一区二区在线| 久久久精品欧美日韩精品| 啦啦啦观看免费观看视频高清| 黄色女人牲交| 亚洲第一电影网av| 精品国产三级普通话版| 蜜桃久久精品国产亚洲av| 国产真实伦视频高清在线观看 | 男人舔女人下体高潮全视频| 老司机福利观看| 自拍偷自拍亚洲精品老妇| 国产美女午夜福利| 久久热精品热| 简卡轻食公司| 日本免费一区二区三区高清不卡| 乱系列少妇在线播放| 全区人妻精品视频| 最后的刺客免费高清国语| 日韩,欧美,国产一区二区三区 | 欧美日韩综合久久久久久 | 国产伦在线观看视频一区| 久久精品综合一区二区三区| 麻豆一二三区av精品| 长腿黑丝高跟| 黄色欧美视频在线观看| 天堂√8在线中文| 国产在线精品亚洲第一网站| 88av欧美| 亚洲三级黄色毛片| 亚洲av成人av| 国产精品电影一区二区三区| 女生性感内裤真人,穿戴方法视频| 精品久久久久久久久久免费视频| 91麻豆精品激情在线观看国产| 欧洲精品卡2卡3卡4卡5卡区| 精品久久久久久久久亚洲 | 丰满人妻一区二区三区视频av| 亚洲美女搞黄在线观看 | 国产精品久久久久久av不卡| 黄色日韩在线| 日韩欧美精品免费久久| 国产精品久久久久久久电影| 亚洲四区av| 精品不卡国产一区二区三区| 男人的好看免费观看在线视频| h日本视频在线播放| 日本爱情动作片www.在线观看 | www日本黄色视频网| 亚洲一区高清亚洲精品| 美女大奶头视频| 亚洲欧美激情综合另类| 村上凉子中文字幕在线| 中出人妻视频一区二区| 美女免费视频网站| 亚洲avbb在线观看| 热99在线观看视频| 亚洲图色成人| 国产蜜桃级精品一区二区三区| 精品99又大又爽又粗少妇毛片 | 久久久久久久亚洲中文字幕| 88av欧美| 成人欧美大片| 99在线人妻在线中文字幕| 欧美一级a爱片免费观看看| 88av欧美| 国产精品电影一区二区三区| 男女下面进入的视频免费午夜| а√天堂www在线а√下载| 久久久国产成人精品二区| 韩国av一区二区三区四区| 精品国产三级普通话版| 成人性生交大片免费视频hd| 亚洲最大成人中文| 国产精品久久久久久精品电影| av女优亚洲男人天堂| 搞女人的毛片| 国产成年人精品一区二区| 成年免费大片在线观看| 亚洲中文字幕一区二区三区有码在线看| 一级毛片久久久久久久久女| 91在线观看av| 欧美一区二区亚洲| 老司机深夜福利视频在线观看| 午夜亚洲福利在线播放| 伦理电影大哥的女人| 一个人观看的视频www高清免费观看| 少妇高潮的动态图| 国产精品久久久久久av不卡| 熟女电影av网| 麻豆成人午夜福利视频| 男人舔奶头视频| 久久精品国产亚洲网站| 日本色播在线视频| 亚洲精品乱码久久久v下载方式| 一a级毛片在线观看| 蜜桃久久精品国产亚洲av| 久久婷婷人人爽人人干人人爱| 高清在线国产一区| 精品乱码久久久久久99久播| 日日夜夜操网爽| 国产黄a三级三级三级人| 99精品在免费线老司机午夜| 亚洲美女视频黄频| 国产视频内射| 日韩一区二区视频免费看| 欧美激情久久久久久爽电影| 欧美性猛交╳xxx乱大交人| 久久久久久久久中文| 一进一出好大好爽视频| 桃色一区二区三区在线观看| 人妻夜夜爽99麻豆av| 69人妻影院| 午夜日韩欧美国产| 哪里可以看免费的av片| 欧美在线一区亚洲| 不卡视频在线观看欧美| 九色成人免费人妻av| 国内久久婷婷六月综合欲色啪| 国产av在哪里看| 日本 欧美在线| av专区在线播放| 亚洲自偷自拍三级| 欧美成人免费av一区二区三区| 大型黄色视频在线免费观看| 中国美女看黄片| 伦精品一区二区三区| 亚洲熟妇熟女久久| 久久久久久久久久久丰满 | 精品久久国产蜜桃| 精品久久久久久久久久久久久| 国产 一区精品| 一级黄片播放器| 三级国产精品欧美在线观看| 少妇猛男粗大的猛烈进出视频 | 99热网站在线观看| 亚洲欧美激情综合另类| 不卡一级毛片| 久久天躁狠狠躁夜夜2o2o| 伊人久久精品亚洲午夜| 亚洲最大成人av| 麻豆一二三区av精品| 国产精品一及| 国产精品永久免费网站| 色在线成人网| 人妻少妇偷人精品九色| 日韩,欧美,国产一区二区三区 | 国产精品久久视频播放| 国产精品亚洲美女久久久| 国产亚洲91精品色在线| 又爽又黄a免费视频| 欧美bdsm另类| 国产精品一区www在线观看 | 91在线精品国自产拍蜜月| 国产伦精品一区二区三区四那| 国产一区二区三区av在线 | 国产精品,欧美在线| 久久精品国产鲁丝片午夜精品 | 国产高清三级在线| 91麻豆精品激情在线观看国产| 日日撸夜夜添| 中文字幕人妻熟人妻熟丝袜美| 久久久久久久久大av| 成人综合一区亚洲| 最近视频中文字幕2019在线8| 特级一级黄色大片| 男女视频在线观看网站免费| 亚洲国产高清在线一区二区三| 久99久视频精品免费| 村上凉子中文字幕在线| 99热这里只有精品一区| 国产精品人妻久久久影院| 啦啦啦观看免费观看视频高清| 噜噜噜噜噜久久久久久91| 国产精品乱码一区二三区的特点| 日韩高清综合在线| 12—13女人毛片做爰片一| 99热6这里只有精品| 国产 一区精品| av黄色大香蕉| 精品久久久久久久久av| netflix在线观看网站| 最近中文字幕高清免费大全6 | av视频在线观看入口| 一级黄色大片毛片| 国产精品美女特级片免费视频播放器| 少妇人妻一区二区三区视频| 久久久久久久久久久丰满 | 一本精品99久久精品77| 琪琪午夜伦伦电影理论片6080| 欧美日韩国产亚洲二区| а√天堂www在线а√下载| 九九在线视频观看精品| 非洲黑人性xxxx精品又粗又长| 日韩中字成人| 欧美一级a爱片免费观看看| 最近视频中文字幕2019在线8| 欧美不卡视频在线免费观看| 国产欧美日韩一区二区精品| av在线观看视频网站免费| 日韩欧美在线二视频| 久久人人爽人人爽人人片va| 亚洲无线观看免费| 欧美成人a在线观看| 不卡一级毛片| 久久久国产成人精品二区| 一进一出好大好爽视频| 人人妻人人看人人澡| 欧美成人性av电影在线观看| 国产一区二区三区视频了| 国产亚洲欧美98| 色综合亚洲欧美另类图片| 天堂av国产一区二区熟女人妻| 国产高清有码在线观看视频| 不卡一级毛片| 久久草成人影院| 在线免费观看的www视频| 日韩在线高清观看一区二区三区 | 亚洲男人的天堂狠狠| 欧美激情国产日韩精品一区| 色精品久久人妻99蜜桃| 亚洲av五月六月丁香网| 69人妻影院| 日韩国内少妇激情av| 观看美女的网站| 欧美潮喷喷水| 国产成年人精品一区二区| 99热这里只有是精品50| 国产亚洲91精品色在线| 特级一级黄色大片| 午夜精品在线福利| 日韩欧美国产一区二区入口| 久久久久久久久久久丰满 | 女人被狂操c到高潮| 97超视频在线观看视频| 熟妇人妻久久中文字幕3abv| 日本一二三区视频观看| 精品人妻视频免费看| 精品久久久久久久久久免费视频| 最近视频中文字幕2019在线8| 在线a可以看的网站| 亚洲人成伊人成综合网2020| 久久久久久久久大av| 亚洲午夜理论影院| 久久久久九九精品影院| 亚洲成人精品中文字幕电影| 国内久久婷婷六月综合欲色啪| 亚洲无线在线观看| 亚洲经典国产精华液单| 一级黄片播放器| 久久久久久伊人网av| 看黄色毛片网站| 中亚洲国语对白在线视频| 国产乱人视频| 日本三级黄在线观看| 国产一区二区亚洲精品在线观看| 国产精品1区2区在线观看.| 午夜激情福利司机影院| 高清日韩中文字幕在线| 亚洲,欧美,日韩| 久久九九热精品免费| 一级a爱片免费观看的视频| 亚洲av成人精品一区久久| 一级毛片久久久久久久久女| 精品一区二区三区人妻视频| 色av中文字幕| 亚洲无线在线观看| 在线天堂最新版资源| 久久99热6这里只有精品| 88av欧美| 欧美激情久久久久久爽电影| 18禁黄网站禁片午夜丰满| 黄片wwwwww| 麻豆一二三区av精品| 97热精品久久久久久| av天堂在线播放| 中文字幕人妻熟人妻熟丝袜美| 欧美性感艳星| 麻豆久久精品国产亚洲av| 91麻豆av在线| 国产成人福利小说| 在线观看66精品国产| 亚洲欧美激情综合另类| 日本精品一区二区三区蜜桃| 欧美激情久久久久久爽电影| 三级毛片av免费| 国产精品,欧美在线| 天天一区二区日本电影三级| 亚洲成人久久爱视频| 久久久久久久久大av| 国产免费男女视频| 能在线免费观看的黄片| 久久久久久久久大av| 内射极品少妇av片p| 国产精品一区www在线观看 | 少妇的逼水好多| 嫁个100分男人电影在线观看| 村上凉子中文字幕在线| 啦啦啦韩国在线观看视频| 色精品久久人妻99蜜桃| 精品久久久久久,| 国产色婷婷99| 中亚洲国语对白在线视频| 日本与韩国留学比较| 亚洲在线自拍视频| 久久亚洲真实| 午夜影院日韩av| 欧美潮喷喷水| 18禁裸乳无遮挡免费网站照片| 村上凉子中文字幕在线| 欧美潮喷喷水| 高清毛片免费观看视频网站| 欧美日韩精品成人综合77777| 国产激情偷乱视频一区二区| www.www免费av| www日本黄色视频网| 日韩欧美在线乱码| 午夜精品在线福利| 久久精品国产亚洲av香蕉五月| 国产伦在线观看视频一区| 国产精品一区www在线观看 | 国产亚洲精品av在线| 麻豆国产av国片精品| www.色视频.com| 韩国av一区二区三区四区| 在线免费十八禁| 国产精品,欧美在线| 欧美不卡视频在线免费观看|