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

    改進(jìn)TSA降噪與平方包絡(luò)譜分析的故障特征提取

    2021-06-10 01:25:05郭遠(yuǎn)晶金曉航魏燕定楊友東
    振動(dòng)工程學(xué)報(bào) 2021年2期
    關(guān)鍵詞:相關(guān)系數(shù)故障診斷

    郭遠(yuǎn)晶 金曉航 魏燕定 楊友東

    摘要: 為實(shí)現(xiàn)旋轉(zhuǎn)機(jī)械設(shè)備故障診斷,提出基于相關(guān)性檢測(cè)的振動(dòng)信號(hào)改進(jìn)時(shí)域同步平均(Time Synchronous Averaging,TSA)降噪方法,并利用平方包絡(luò)譜提取故障特征。在采樣振動(dòng)信號(hào)中任取一段參考子信號(hào),令參考子信號(hào)從振動(dòng)信號(hào)初始時(shí)間開始,沿時(shí)間軸逐點(diǎn)滑動(dòng)至終了時(shí)間,同時(shí)計(jì)算參考子信號(hào)與其在振動(dòng)信號(hào)中遮掩子信號(hào)之間的相關(guān)系數(shù);利用步長(zhǎng)迭代算法獲取相關(guān)系數(shù)最優(yōu)閾值,選擇與參考子信號(hào)相關(guān)系數(shù)大于最優(yōu)閾值的所有遮掩子信號(hào)作為總體同步信號(hào),平均計(jì)算后獲得降噪的目標(biāo)信號(hào);對(duì)目標(biāo)信號(hào)進(jìn)行平方包絡(luò)譜分析。仿真信號(hào)、齒輪和軸承振動(dòng)信號(hào)的分析結(jié)果表明,所提出的方法能夠有效用于旋轉(zhuǎn)機(jī)械設(shè)備的振動(dòng)信號(hào)降噪與故障特征提取。

    關(guān)鍵詞: 故障診斷; 時(shí)域同步平均; 相關(guān)系數(shù); 步長(zhǎng)迭代算法; 平方包絡(luò)譜

    中圖分類號(hào): TH165+.3; TN911.7 ? ?文獻(xiàn)標(biāo)志碼: A ? ?文章編號(hào): 1004-4523(2021)02-0402-09

    DOI:10.16385/j.cnki.issn.1004-4523.2021.02.021

    引 ?言

    旋轉(zhuǎn)機(jī)械設(shè)備出現(xiàn)故障損傷時(shí),振動(dòng)信號(hào)中通常包含有規(guī)律性的故障特征,其提取對(duì)于設(shè)備故障診斷具有重要意義。但由于噪聲干擾、傳遞路徑衰減、多振源耦合等因素[1],往往需要對(duì)振動(dòng)信號(hào)進(jìn)行降噪處理,如濾波器降噪、小波閾值降噪、SVD降噪以及經(jīng)驗(yàn)?zāi)B(tài)分解降噪等,才能順利提取出故障特征。此外,在旋轉(zhuǎn)機(jī)械設(shè)備平穩(wěn)運(yùn)行工況下,適應(yīng)于振動(dòng)信號(hào)的周期性變化特征,時(shí)域同步平均(Time Synchronous Averaging,TSA)降噪方法也得到了廣泛應(yīng)用[2?7]。

    振動(dòng)信號(hào)TSA降噪的關(guān)鍵在于參與平均計(jì)算的各段振動(dòng)信號(hào)要保持同步,即無相位差。目前保持振動(dòng)信號(hào)同步的方法主要有按整數(shù)周期截取信號(hào)和時(shí)標(biāo)脈沖觸發(fā)同步采樣[8]。按整數(shù)周期截取信號(hào)的方法需要獲取振動(dòng)信號(hào)的準(zhǔn)確周期[7],但這并不是一件容易的事情,原因在于旋轉(zhuǎn)機(jī)械設(shè)備的轉(zhuǎn)速通常不能保持時(shí)刻恒定,尤其是在故障損傷狀態(tài)下,轉(zhuǎn)速的波動(dòng)使振動(dòng)信號(hào)沒有一個(gè)嚴(yán)格穩(wěn)定的周期,因此,按照一個(gè)固定周期對(duì)振動(dòng)信號(hào)進(jìn)行截?cái)啵瑫?huì)出現(xiàn)相位誤差累積效應(yīng)[9],影響降噪效果。時(shí)標(biāo)脈沖觸發(fā)同步采樣則需要在參考軸上安裝轉(zhuǎn)速計(jì)或光電編碼器,輸出時(shí)標(biāo)脈沖觸發(fā)振動(dòng)信號(hào)的同步采樣或者插值重采樣[2,5,10?11]。此類方法能夠保證采樣振動(dòng)信號(hào)的準(zhǔn)確同步,但需要額外增加硬件成本,并且在很多實(shí)際應(yīng)用場(chǎng)合中,轉(zhuǎn)速計(jì)或光電編碼器的安裝也并非易事[12],因此,此類方法的應(yīng)用也受到一定限制。

    為了獲取時(shí)域同步振動(dòng)信號(hào)參與TSA降噪,本文從旋轉(zhuǎn)機(jī)械設(shè)備上采集的振動(dòng)信號(hào)中選取一段參考子信號(hào),令參考子信號(hào)從振動(dòng)信號(hào)的初始時(shí)間開始,沿時(shí)間軸滑動(dòng)至終了時(shí)間,期間參考子信號(hào)不斷地在振動(dòng)信號(hào)中遮掩一段等長(zhǎng)度的子信號(hào),同時(shí)計(jì)算參考子信號(hào)與遮掩子信號(hào)之間的相關(guān)系數(shù);然后利用步長(zhǎng)迭代算法獲取相關(guān)系數(shù)的最優(yōu)閾值;最后選擇與參考子信號(hào)相關(guān)系數(shù)大于最優(yōu)閾值的所有遮掩子信號(hào)作為總體同步信號(hào),平均計(jì)算后獲得降噪的目標(biāo)信號(hào)。對(duì)目標(biāo)信號(hào)進(jìn)行平方包絡(luò)譜分析,提取故障特征,判別旋轉(zhuǎn)機(jī)械設(shè)備故障狀態(tài)。

    1 TSA降噪原理

    旋轉(zhuǎn)機(jī)械設(shè)備運(yùn)行時(shí)的振動(dòng)信號(hào)為y(t),采樣時(shí)間間隔為Δt,則采樣振動(dòng)信號(hào)為y(kΔt),記作y(k),k=0,1,…,N(N為振動(dòng)信號(hào)采樣數(shù)據(jù)點(diǎn)數(shù))。若y(k)由周期為NT(NT為一個(gè)時(shí)間周期內(nèi)的采樣數(shù)據(jù)點(diǎn)數(shù))的特征信號(hào)s(k)和白噪聲sn(k)組成,即

    式中 ?fr為齒輪轉(zhuǎn)頻,fm為齒輪嚙合頻率,N為沖擊特征個(gè)數(shù),An為第n個(gè)沖擊特征的最大幅值,[tn ,tn+1)表示第n個(gè)沖擊特征的持續(xù)時(shí)間范圍。

    假設(shè)一對(duì)嚙合齒輪中發(fā)生局部損傷齒輪的齒數(shù)z=20,轉(zhuǎn)頻fr=25 Hz,嚙合頻率fm=zfr=500 Hz。由于齒輪發(fā)生局部損傷的故障特征頻率與其轉(zhuǎn)頻fr一致,且考慮諧波分量sh(t)中調(diào)制信號(hào)的初相位為π/2,經(jīng)過計(jì)算可知,第1個(gè)沖擊特征出現(xiàn)的時(shí)間t1=0.03 s,持續(xù)時(shí)間(t2-t1)=1/fr=0.04 s。沖擊特征的最大幅值序列An由1.5?1.7的均勻分布隨機(jī)序列產(chǎn)生。設(shè)定信號(hào)采樣頻率fs=10 kHz,采樣時(shí)間t=12 s,疊加高斯白噪聲,使信噪比為-8 dB,由此產(chǎn)生長(zhǎng)度L=120000的仿真振動(dòng)信號(hào)y(k)(k=1,2,…,L),如圖3(a)所示。

    此處期望從y(k)中選取的參考子信號(hào)包含5個(gè)周期的沖擊特征,則參考子信號(hào)長(zhǎng)度至少為Np×(1/fr)×fs=2000。然后在y(k)中選擇一個(gè)時(shí)間點(diǎn)k=58930,對(duì)應(yīng)時(shí)間t=5.8930 s,以此開始截取一段長(zhǎng)度為2000的參考子信號(hào),如圖3(a)中間框選的信號(hào)段所示,參考子信號(hào)的時(shí)域波形如圖3(b)所示,其中的沖擊特征完全被噪聲所淹沒。參考子信號(hào)對(duì)應(yīng)的無噪波形如圖3(c)所示。

    令參考子信號(hào)從仿真振動(dòng)信號(hào)y(k)的初始時(shí)間點(diǎn)k=1,即初始時(shí)間t=0開始,如圖3(a)左側(cè)框選位置所示,沿時(shí)間軸逐點(diǎn)滑動(dòng)至終了時(shí)間點(diǎn)k=118001,即滑動(dòng)終了時(shí)間t=11.8001 s,如圖3(a)右側(cè)框選位置所示。參考子信號(hào)滑動(dòng)至每個(gè)時(shí)間點(diǎn)處遮掩一段等長(zhǎng)度的子信號(hào),計(jì)算參考子信號(hào)與各個(gè)遮掩子信號(hào)之間的相關(guān)系數(shù)。將相關(guān)系數(shù)從小到大排列,得到重排的相關(guān)系數(shù)序列ρ ?_sm (k_t)(kt=1,2,…,118001),舍棄其負(fù)值及末尾Nr=10個(gè)取值后的曲線如圖3(d)所示。

    需要說明的是,ρ_sm^*的選取帶有一定主觀性,選取的ρ_sm^*只是相對(duì)最優(yōu),并不是唯一,但只要ρ_sm^*在一定范圍內(nèi),如在本仿真實(shí)例中,只需滿足0.10≤ρ_sm^*≤0.15,就可以使改進(jìn)TSA取得良好的降噪效果。R(r)關(guān)于ρth的關(guān)系曲線則具有良好的指示意義,借助它可以方便地選取出合適的最優(yōu)閾值ρ_sm^*,用于篩選符合要求的總體同步信號(hào),實(shí)現(xiàn)改進(jìn)TSA降噪。改進(jìn)TSA降噪后目標(biāo)信號(hào)的平方包絡(luò)譜如圖3(g)所示,可見故障特征頻率25 Hz被順利提取出來。

    對(duì)比圖3(b)和3(f)可見,仿真振動(dòng)信號(hào)中的噪聲得到了較好的抑制。但對(duì)比圖3(c)和3(f)可見,降噪后的目標(biāo)信號(hào)幅值存在一定的衰減,原因?yàn)樵诳傮w同步信號(hào)中,以參考子信號(hào)為參考,雖然超前的遮掩子信號(hào)與滯后的遮掩子信號(hào)之和無相位差,但兩者在幅值上存在一定的相互抵消,因此兩者求和之后再求平均,獲得的幅值要比參考子信號(hào)小。

    為進(jìn)一步驗(yàn)證改進(jìn)TSA降噪方法的有效性,在y(k)中再任選一個(gè)時(shí)間點(diǎn)k=12689,對(duì)應(yīng)時(shí)間t=1.2689 s,以此開始截取一段長(zhǎng)度同樣為2000的參考子信號(hào),按前述相同步驟獲取降噪后的目標(biāo)信號(hào)如圖4(a)所示,其平方包絡(luò)譜如圖4(b)所示。在y(k)中改選一個(gè)時(shí)間點(diǎn)k=102371,對(duì)應(yīng)時(shí)間t=10.2371 s,按相同方法獲取的降噪后目標(biāo)信號(hào)及其平方包絡(luò)譜分別如圖5(a)和5(b)所示。通過對(duì)比可以發(fā)現(xiàn),參考子信號(hào)選取的時(shí)間點(diǎn)不同,對(duì)最終獲取目標(biāo)信號(hào)的降噪效果沒有影響,只是目標(biāo)信號(hào)的相位可能會(huì)發(fā)生改變,但這對(duì)目標(biāo)信號(hào)平方包絡(luò)譜中故障特征的提取并無影響。

    5 實(shí)際振動(dòng)信號(hào)分析

    5.1 齒輪故障振動(dòng)信號(hào)分析

    齒輪故障振動(dòng)信號(hào)取自一臺(tái)風(fēng)力發(fā)電機(jī),其齒輪箱中一個(gè)小齒輪發(fā)生輪齒斷裂損傷故障,如圖6所示。該小齒輪的齒數(shù)z=32,額定轉(zhuǎn)速n=1770 r/min,轉(zhuǎn)頻fr=29.5 Hz,故障特征頻率fg=fr=29.5 Hz,嚙合頻率fm=z·fn=960 Hz,振動(dòng)信號(hào)采樣頻率fs=97656 Hz,采樣得到的一段振動(dòng)信號(hào)如圖7(a)所示。

    由于該小齒輪的故障特征頻率為29.5 Hz,此處期望截取的參考子信號(hào)中至少包含6個(gè)故障沖擊特征,那么參考子信號(hào)的長(zhǎng)度至少為6×(1/29.5)×fs=19531。因此,在采樣的振動(dòng)信號(hào)中任選一個(gè)時(shí)間點(diǎn)k=50836,對(duì)應(yīng)時(shí)間為t=0.5206 s,以此開始截取一段長(zhǎng)度為19854的參考子信號(hào),時(shí)域波形如圖7(b)所示,其中噪聲較大,故障特征并不顯著。令參考子信號(hào)從采樣振動(dòng)信號(hào)初始時(shí)間開始,沿時(shí)間軸逐點(diǎn)滑動(dòng),計(jì)算參考子信號(hào)與各時(shí)間點(diǎn)處遮掩子信號(hào)之間的相關(guān)系數(shù)?;瑒?dòng)結(jié)束后,將相關(guān)系數(shù)從小到大排列,得到重排相關(guān)系數(shù)序列,舍棄其負(fù)值及末尾Nr=10個(gè)取值后的曲線如圖7(c)所示。

    選擇最優(yōu)閾值ρ_sm^*的搜尋區(qū)間為[0.3ρ_sm^max, ρ_sm^max],即[0.0997,0.3324],設(shè)定迭代步數(shù)M=100,迭代計(jì)算后獲得均方根誤差R(r)關(guān)于閾值變量ρth的關(guān)系曲線,如圖7(d)所示。選取R(r)在低值平穩(wěn)變化過程中最小值處的閾值ρth=0.1928作為最優(yōu)閾值ρ_sm^*,然后選擇與參考子信號(hào)相關(guān)系數(shù)ρsm>ρ_sm^*=0.1928的所有遮掩子信號(hào)作為總體同步信號(hào),平均計(jì)算后獲得降噪的目標(biāo)信號(hào),如圖7(e)所示,可見其中顯著的周期性故障沖擊特征。

    對(duì)降噪后的目標(biāo)信號(hào)進(jìn)行平方包絡(luò)譜分析,如圖7(f)所示,可見小齒輪故障特征頻率fg及其倍頻(2?13)×fg被清晰地提取出來。目標(biāo)信號(hào)的包絡(luò)譜如圖7(g)所示,故障特征頻率fg及其倍頻的譜線不如平方包絡(luò)譜清晰,因此,相比于包絡(luò)譜,平方包絡(luò)譜在振動(dòng)信號(hào)分析中具有一定的優(yōu)勢(shì)。

    作為比較,采用基于周期分段的傳統(tǒng)TSA降噪方法對(duì)采樣振動(dòng)信號(hào)進(jìn)行處理,獲得的降噪信號(hào)如圖8(a)所示,其平方包絡(luò)譜與包絡(luò)譜分別如圖8(b)和8(c)所示。從時(shí)域上講,由于齒輪故障損傷會(huì)引起轉(zhuǎn)速波動(dòng),振動(dòng)信號(hào)的周期并非嚴(yán)格恒定,按周期分段截取的各段振動(dòng)信號(hào)存在一定的相位誤差,且相位誤差具有累計(jì)效應(yīng),相比于改進(jìn)TSA降噪方法,傳統(tǒng)TSA降噪方法獲取的目標(biāo)信號(hào)降噪不夠徹底,存在較多的高頻分量,其中的周期性故障沖擊特征不太容易被識(shí)別,信號(hào)幅值也存在較大程度的衰減。從頻域上講,改進(jìn)TSA降噪方法與傳統(tǒng)TSA降噪方法對(duì)于故障特征頻率的提取效果基本相同,且平方包絡(luò)譜的故障特征頻率及其倍頻譜線要比包絡(luò)譜更加清晰可辨。

    5.2 軸承故障振動(dòng)信號(hào)分析

    軸承故障振動(dòng)信號(hào)取自Case Western Reserve University軸承數(shù)據(jù)中心。選擇試驗(yàn)臺(tái)驅(qū)動(dòng)端的NTN滾動(dòng)軸承為研究對(duì)象,其內(nèi)滾道上有一直徑為0.7112 mm(0.028"),深度為1.27 mm(0.050")的局部損傷故障。在試驗(yàn)臺(tái)運(yùn)行過程中,驅(qū)動(dòng)電機(jī)功率為2.2 kW,轉(zhuǎn)速為1730 r/min,加載電機(jī)施加3 N·m載荷,軸承內(nèi)圈的故障特征頻率fBPFI=123 Hz。振動(dòng)數(shù)據(jù)采樣頻率為12 kHz,采樣到的一段振動(dòng)信號(hào)如圖9(a)所示。

    在采樣振動(dòng)信號(hào)中任選一個(gè)時(shí)間點(diǎn)k=20185,對(duì)應(yīng)時(shí)間為t=1.6821 s,以此開始截取一段長(zhǎng)度為2000的參考子信號(hào),如圖9(b)所示,其中故障沖擊特征并不顯著。令參考子信號(hào)從采樣的振動(dòng)信號(hào)初始時(shí)間開始,沿時(shí)間軸逐點(diǎn)滑動(dòng),計(jì)算參考子信號(hào)與遮掩子信號(hào)之間的相關(guān)系數(shù)。滑動(dòng)結(jié)束后獲得參考子信號(hào)與所有遮掩子信號(hào)之間的相關(guān)系數(shù),將其從小到大排列后得到重排相關(guān)系數(shù)序列,舍棄其負(fù)值及末尾Nr=10個(gè)取值后的曲線如圖9(c)所示。

    選擇最優(yōu)閾值ρ_sm^*的搜尋區(qū)間為[0?ρ_sm^max,ρ_sm^max],即[0,0.5822],設(shè)定迭代步數(shù)M=100,迭代計(jì)算后獲得均方根誤差R(r)關(guān)于閾值變量ρth的關(guān)系曲線,如圖9(d)所示。選取R(r)從低值平穩(wěn)變化向急劇增大變化轉(zhuǎn)變的一個(gè)臨界值ρth=0.4366作為最優(yōu)閾值ρ_sm^*,然后選擇與參考子信號(hào)相關(guān)系數(shù)ρsm>ρ_sm^*=0.4366的所有遮掩子信號(hào)作為總體同步信號(hào),平均計(jì)算后獲得降噪的目標(biāo)信號(hào)如圖9(e)所示,可見噪聲得到了較好的抑制,周期性特征非常顯著。降噪后目標(biāo)信號(hào)的平方包絡(luò)譜如圖9(f)所示,明顯可見軸承內(nèi)圈的故障特征頻率fBPFI及其倍頻(2?5)×f_BPFI。降噪后目標(biāo)信號(hào)的包絡(luò)譜如圖9(g)所示,故障特征頻率fBPFI及其倍頻的譜線也基本清晰可辨,但不如平方包絡(luò)譜突出。

    作為比較,采用傳統(tǒng)TSA降噪方法對(duì)軸承的采樣振動(dòng)信號(hào)進(jìn)行處理,獲得降噪的目標(biāo)信號(hào)如圖10(a)所示。從該目標(biāo)信號(hào)中能夠大致看出信號(hào)的周期性,但難以辨別出周期性的故障沖擊特征,原因在于軸承故障損傷引起轉(zhuǎn)速波動(dòng),導(dǎo)致采樣振動(dòng)信號(hào)的周期并非嚴(yán)格恒定,按周期分段截取的各段時(shí)域同步信號(hào)在求和過程中存在較大的相位誤差累計(jì)效應(yīng)。相比于改進(jìn)TSA降噪方法,傳統(tǒng)TSA降噪方法獲取的目標(biāo)信號(hào)降噪不夠徹底,存在較多的高頻成分以及故障特征頻率的0.5倍、1.5倍、2.5倍等分?jǐn)?shù)倍頻成分,這從傳統(tǒng)TSA降噪目標(biāo)信號(hào)的平方包絡(luò)譜(如圖10(b)所示)與包絡(luò)譜(如圖10(c)所示)中也可以看出來。雖然傳統(tǒng)TSA降噪結(jié)合平方包絡(luò)譜或包絡(luò)譜也能夠提取出故障特征頻率,但由于故障特征頻率分?jǐn)?shù)倍頻成分的存在,對(duì)故障特征頻率識(shí)別有一定的干擾,影響故障判別。

    6 結(jié) ?論

    (1)提出一種基于相關(guān)性檢測(cè)的改進(jìn)TSA降噪方法,該方法在旋轉(zhuǎn)機(jī)械設(shè)備的采樣振動(dòng)信號(hào)中任選一段參考子信號(hào),令參考子信號(hào)在采樣振動(dòng)信號(hào)中逐點(diǎn)滑動(dòng),計(jì)算參考子信號(hào)與遮掩子信號(hào)之間的相關(guān)系數(shù);然后利用步長(zhǎng)迭代算法獲取相關(guān)系數(shù)的最優(yōu)閾值;最后選擇與參考子信號(hào)相關(guān)系數(shù)大于最優(yōu)閾值的所有遮掩子信號(hào)作為總體同步信號(hào),平均計(jì)算后可以獲得降噪效果良好的目標(biāo)信號(hào)。

    (2)改進(jìn)TSA降噪方法成功應(yīng)用于齒輪與軸承故障振動(dòng)信號(hào)的處理,降噪后得到的目標(biāo)信號(hào)經(jīng)過平方包絡(luò)譜分析,能夠有效提取出齒輪與軸承的故障特征頻率,從而實(shí)現(xiàn)旋轉(zhuǎn)機(jī)械設(shè)備故障診斷。

    (3)改進(jìn)TSA降噪方法對(duì)旋轉(zhuǎn)機(jī)械設(shè)備轉(zhuǎn)速的穩(wěn)定性以及振動(dòng)信號(hào)的平穩(wěn)性沒有嚴(yán)格要求,無需獲取振動(dòng)信號(hào)的準(zhǔn)確周期,避免了振動(dòng)信號(hào)按周期分段而導(dǎo)致的相位誤差累積效應(yīng);無需轉(zhuǎn)速傳感器,硬件成本低、應(yīng)用場(chǎng)合廣,具有較好的適用性。

    參考文獻(xiàn):

    [1] 熊國(guó)良,胡俊鋒,陳慧,等. 基于SK-NLM包絡(luò)的滾動(dòng)軸承故障沖擊特征增強(qiáng)[J]. 儀器儀表學(xué)報(bào),2016, 37(10):2176-2184.

    Xiong Guoliang, Hu Junfeng, Chen Hui, et al. Rolling bearing fault impact feature enhancement based on spectral kurtosis and non-local means[J]. Chinese Journal of Scientific Instrument, 2016, 37(10): 2176-2184.

    [2] Ahamed N, Pandya Y, Parey A. Spur gear tooth root crack detection using time synchronous averaging under fluctuating speed[J]. Measurement, 2014, 52: 1-11.

    [3] Ha J M, Youn B D, Oh H, et al. Autocorrelation-based time synchronous averaging for condition monitoring of planetary gearboxes in wind turbines[J]. Mechanical Systems and Signal Processing, 2015, 70-71: 161-175.

    [4] Li H, Xu F, Liu H, et al. Incipient fault information determination for rolling element bearing based on averaging reassigned wavelet scalogram[J]. Measurement, 2015, 65: 1-10.

    [5] Roy S K, Mohanty A R, Kumar C S. Fault detection in a multistage gearbox by time synchronous averaging of the instantaneous angular speed[J]. Journal of Vibration and Control, 2016, 22(2): 468-480.

    [6] Camerini V, Coppotelli G, Bendisch S, et al. Impact of pulse time uncertainty on synchronous average: Statistical analysis and relevance to rotating machinery diagnosis[J]. Mechanical Systems and Signal Processing, 2019, 129: 308-336.

    [7] 龍海軍,孫燦飛,莫固良. 直升機(jī)振動(dòng)檢測(cè)通用算法的研究與實(shí)現(xiàn)[J]. 振動(dòng)、測(cè)試與診斷,2016,36(3):524-528.

    Long Haijun, Sun Canfei, Mo Guliang. Research and implementation of general algorithm for helicopter vibration detection[J]. Journal of Vibration, Measurement & Diagnosis, 2016, 36(3):524-528.

    [8] 邵毅敏,周曉君,歐家福. 變周期信號(hào)的時(shí)域同步平均新算法[J]. 振動(dòng)工程學(xué)報(bào),2009,22(5):497-502.

    Shao Yimin, Zhou Xiaojun, Ou Jiafu. A algorithm of time domain averaging for fluctuant-period signal[J]. Journal of Vibration Engineering, 2009, 22(5):497-502.

    [9] 沈國(guó)際,陶利民,徐永成. 時(shí)域同步平均的相位誤差累積效應(yīng)研究[J]. 振動(dòng)工程學(xué)報(bào),2007,20(4):335-339.

    Shen Guoji, Tao Limin, Xu Yongcheng. Research on phase error accumulation effect of time synchronous averaging[J]. Journal of Vibration Engineering, 2007, 20(4):335-339.

    [10] 羅德?lián)P. 時(shí)域同步平均原理與應(yīng)用[J]. 振動(dòng)、測(cè)試與診斷,1999,19(3):202-207.

    Luo Deyang. Principals and applications of time domain synchronous averaging[J]. Journal of Vibration, Measurement & Diagnosis, 1999, 19(3): 202-207.

    [11] 呂琛,宋希庚,鄒積斌. 基于DSP的振動(dòng)信號(hào)階比與時(shí)域同步平均分析[J]. 振動(dòng)與沖擊,2002,21(2):53-57.

    L? Chen, SONG Xigeng, ZOU Jibin. DSP-based order domain and time domain synchronous averaging analysis of vibration signal[J]. Journal of Vibration and Shock, 2002, 21(2): 53-57.

    [12] Combet F, Gelman L. An automated methodology for performing time synchronous averaging of a gearbox signal without speed sensor[J]. Mechanical Systems and Signal Processing, 2007, 21(6): 2590-2606.

    [13] 何正嘉,劉雄,屈梁生. 信號(hào)時(shí)域平均原理和應(yīng)用[J]. 信號(hào)處理,1986,2(4):46-53.

    He Zhengjia, Liu Xiong, Qu Liangsheng. The principle and application of signal time domain average[J]. Signal Processing, 1986, 2(4): 46-53.

    [14] 郭遠(yuǎn)晶,魏燕定,周曉軍. 基于STFT時(shí)頻譜系數(shù)收縮的信號(hào)降噪方法[J]. 振動(dòng)、測(cè)試與診斷,2015,35(6):1090-1096.

    Guo Yuanjing, Wei Yanding, Zhou Xiaojun. Signal denoising method based on STFT time-frequency spectrum coefficients shrinkage[J]. Journal of Vibration, Measurement & Diagnosis, 2015, 35(6): 1090-1096.

    [15] Antoni J. Cyclic spectral analysis of rolling-element bearing signals: Facts and fictions[J]. Journal of Sound and Vibration, 2007, 304(3-5): 497-529.

    [16] Borghesani P, Pennacchi P, Chatterton S. The relationship between kurtosis- and envelope-based indexes for the diagnostic of rolling element bearings[J]. Mechanical Systems and Signal Processing, 2014, 43(1-2): 25-43.

    Fault feature extraction based on improved TSA denoising and squared envelope spectrum

    GUO Yuan-jing1, JIN Xiao-hang2, WEI Yan-ding3, YANG You-dong1

    (1. Zhijiang College, Zhejiang University of Technology, Shaoxing 312030, China;

    2. College of Mechanical Engineering, Zhejiang University of Technology, Hangzhou 310023, China;

    3. Key Laboratory of Advanced Manufacturing Technology of Zhejiang Province, Zhejiang University, Hangzhou 310027, China)

    Abstract: In order to achieve the fault diagnosis of the rotating machine equipment, a method based on improved time synchronous averaging (TSA) denoising using correlation detection and fault features extracting using squared envelope spectrum is proposed. A reference sub-signal is arbitrarily taken from the sampled vibration signal. The reference sub-signal starts from the initial time of the vibration signal and moves point by point along the time axis to the end time. Meanwhile, the correlation coefficient between the reference sub-signal and its masking sub-signal in the vibration signal is calculated. The correlation coefficient optimal threshold are obtained using the step iterative algorithm. All the masking sub-signals whose correlation coefficients are greater than the optimal threshold are selected as overall synchronous signals whose arithmetic average is calculated for the denoised target signal. The target signal analysis is achieved using squared envelope spectrum. The analysis results of the simulated signal, gear and rolling bearing fault vibration signals have shown that the proposed method is an efficient implementation for rotating machine vibration signal denoising and fault feature extraction.

    Key words: fault diagnosis; time synchronous averaging; correlation coefficient; step iterative algorithm; squared envelope spectrum

    作者簡(jiǎn)介: 郭遠(yuǎn)晶(1987-),男,博士,講師。電話:(0575)81112562;E-mail:gyjyn@126.com

    猜你喜歡
    相關(guān)系數(shù)故障診斷
    凍干機(jī)常見故障診斷與維修
    人口老齡化對(duì)我國(guó)消費(fèi)結(jié)構(gòu)影響研究
    南京市能見度變化趨勢(shì)及其影響因素
    基于量子萬有引力搜索的SVM自駕故障診斷
    因果圖定性分析法及其在故障診斷中的應(yīng)用
    秦皇島海域夜光藻種群密度與環(huán)境因子的關(guān)系
    電子鼻傳感器陣列優(yōu)化對(duì)豬肉新鮮度法的檢測(cè)
    肉類研究(2015年5期)2015-08-08 12:48:29
    最大生成樹聚類算法研究
    食用菌菌渣基礎(chǔ)特性分析
    基于LCD和排列熵的滾動(dòng)軸承故障診斷
    免费av不卡在线播放| 欧美高清成人免费视频www| 亚洲最大成人手机在线| 国产精品一及| 亚洲自拍偷在线| 欧美区成人在线视频| 91精品一卡2卡3卡4卡| 美女xxoo啪啪120秒动态图| 七月丁香在线播放| 人体艺术视频欧美日本| 直男gayav资源| 18禁在线无遮挡免费观看视频| 亚洲精品国产成人久久av| 成人毛片60女人毛片免费| 国产探花极品一区二区| 成人美女网站在线观看视频| a级毛色黄片| 欧美国产精品一级二级三级 | 亚洲成人精品中文字幕电影| 在线精品无人区一区二区三 | 国产黄a三级三级三级人| 汤姆久久久久久久影院中文字幕| 日韩成人伦理影院| 中文在线观看免费www的网站| 日韩一区二区三区影片| 亚洲va在线va天堂va国产| 天美传媒精品一区二区| 人妻夜夜爽99麻豆av| 亚洲国产成人一精品久久久| 中文欧美无线码| 国产高清三级在线| av在线播放精品| 蜜臀久久99精品久久宅男| 国产精品一二三区在线看| 久热这里只有精品99| 亚洲aⅴ乱码一区二区在线播放| 99热这里只有精品一区| 日韩制服骚丝袜av| 亚洲天堂国产精品一区在线| 色吧在线观看| tube8黄色片| 一级片'在线观看视频| 国产老妇女一区| 亚洲欧美清纯卡通| 亚洲精品乱久久久久久| av国产免费在线观看| 天天躁夜夜躁狠狠久久av| 最近中文字幕2019免费版| 欧美+日韩+精品| 久热这里只有精品99| 国产成人精品一,二区| 国产精品国产三级国产专区5o| 亚洲av一区综合| 欧美精品国产亚洲| 亚洲av成人精品一二三区| 亚洲一区二区三区欧美精品 | 精品一区二区免费观看| 亚洲欧美精品自产自拍| 一个人看视频在线观看www免费| 一级毛片aaaaaa免费看小| 免费看日本二区| 久久精品国产亚洲网站| 老司机影院成人| 综合色丁香网| 精品一区在线观看国产| 99久久精品热视频| 亚洲一区二区三区欧美精品 | 99视频精品全部免费 在线| 精品亚洲乱码少妇综合久久| 久热久热在线精品观看| 亚洲欧美一区二区三区国产| 久久久久九九精品影院| 国产男人的电影天堂91| 男的添女的下面高潮视频| av天堂中文字幕网| 十八禁网站网址无遮挡 | 亚洲综合色惰| 内射极品少妇av片p| 国产精品一区www在线观看| 婷婷色麻豆天堂久久| 日本三级黄在线观看| 91在线精品国自产拍蜜月| 日韩人妻高清精品专区| 99久久人妻综合| 在线观看美女被高潮喷水网站| 亚洲第一区二区三区不卡| 国产爽快片一区二区三区| 国精品久久久久久国模美| 男女啪啪激烈高潮av片| 精品少妇黑人巨大在线播放| 少妇的逼水好多| 不卡视频在线观看欧美| 乱码一卡2卡4卡精品| 婷婷色av中文字幕| 青春草国产在线视频| 精品一区二区三区视频在线| 免费观看的影片在线观看| 在线观看免费高清a一片| 激情 狠狠 欧美| 欧美国产精品一级二级三级 | 婷婷色综合www| 91久久精品国产一区二区成人| 听说在线观看完整版免费高清| 自拍欧美九色日韩亚洲蝌蚪91 | 欧美日韩国产mv在线观看视频 | 男女那种视频在线观看| 亚洲精品456在线播放app| 亚洲最大成人手机在线| 能在线免费看毛片的网站| 丰满少妇做爰视频| 丝袜美腿在线中文| 国产成人福利小说| 国产男人的电影天堂91| 舔av片在线| 欧美国产精品一级二级三级 | 好男人视频免费观看在线| 亚洲国产欧美在线一区| 国产精品三级大全| 国产一区亚洲一区在线观看| 日本色播在线视频| 亚洲精品亚洲一区二区| 97超碰精品成人国产| 国产成人精品一,二区| 久久精品国产鲁丝片午夜精品| 国产精品国产三级国产专区5o| 狂野欧美白嫩少妇大欣赏| 熟女人妻精品中文字幕| 午夜激情福利司机影院| 久久99热6这里只有精品| 啦啦啦中文免费视频观看日本| 久久韩国三级中文字幕| 美女内射精品一级片tv| 一级毛片黄色毛片免费观看视频| 日韩一区二区视频免费看| 国产在线男女| 一个人看的www免费观看视频| 交换朋友夫妻互换小说| 在线观看av片永久免费下载| 午夜福利在线在线| 国产免费福利视频在线观看| 色网站视频免费| 老司机影院毛片| 成人高潮视频无遮挡免费网站| 国产一区二区三区综合在线观看 | 国产精品精品国产色婷婷| 欧美成人一区二区免费高清观看| 一级黄片播放器| 日本三级黄在线观看| 久久精品国产鲁丝片午夜精品| 免费观看无遮挡的男女| 能在线免费看毛片的网站| 久久久色成人| 日本一二三区视频观看| 一二三四中文在线观看免费高清| 国产成人a∨麻豆精品| 美女内射精品一级片tv| 日本一二三区视频观看| 好男人在线观看高清免费视频| 成年人午夜在线观看视频| 麻豆成人午夜福利视频| 少妇人妻精品综合一区二区| 女人被狂操c到高潮| 中文字幕亚洲精品专区| 久久精品国产鲁丝片午夜精品| 九九爱精品视频在线观看| av黄色大香蕉| 51国产日韩欧美| 男男h啪啪无遮挡| 午夜福利视频1000在线观看| 国产精品嫩草影院av在线观看| 精品一区二区三区视频在线| 精品人妻熟女av久视频| 男插女下体视频免费在线播放| 日本一二三区视频观看| 五月天丁香电影| 一区二区三区乱码不卡18| 春色校园在线视频观看| 国产片特级美女逼逼视频| 亚洲欧美一区二区三区国产| 69人妻影院| 免费大片18禁| 精品久久久久久久末码| 亚洲综合精品二区| 亚洲精品成人av观看孕妇| 特大巨黑吊av在线直播| 乱码一卡2卡4卡精品| 日韩一区二区视频免费看| 久久鲁丝午夜福利片| 边亲边吃奶的免费视频| 99久久中文字幕三级久久日本| 欧美日韩视频高清一区二区三区二| 22中文网久久字幕| 国产探花极品一区二区| 欧美日韩综合久久久久久| 全区人妻精品视频| 22中文网久久字幕| 亚洲精品久久久久久婷婷小说| 国产黄频视频在线观看| 一区二区三区精品91| 少妇丰满av| 18禁裸乳无遮挡动漫免费视频 | 一级二级三级毛片免费看| 久久久久久久精品精品| 69av精品久久久久久| 好男人视频免费观看在线| 好男人视频免费观看在线| 能在线免费看毛片的网站| 亚洲欧美一区二区三区国产| 免费在线观看成人毛片| 亚洲人与动物交配视频| 亚洲精品456在线播放app| 欧美 日韩 精品 国产| 成人美女网站在线观看视频| av专区在线播放| 伦精品一区二区三区| 亚洲欧美日韩东京热| 国产探花极品一区二区| 精品国产一区二区三区久久久樱花 | 在线观看三级黄色| 综合色av麻豆| 国产精品久久久久久av不卡| 国产亚洲av嫩草精品影院| 一级黄片播放器| av线在线观看网站| 熟女av电影| 国产精品久久久久久久电影| 男人舔奶头视频| 男的添女的下面高潮视频| 亚洲成人av在线免费| 狂野欧美白嫩少妇大欣赏| 国产av码专区亚洲av| 夫妻午夜视频| 人人妻人人爽人人添夜夜欢视频 | 亚洲精品日韩av片在线观看| 亚洲国产欧美人成| 亚洲最大成人手机在线| 高清在线视频一区二区三区| 国产欧美亚洲国产| 最新中文字幕久久久久| 嫩草影院新地址| 国产亚洲一区二区精品| 老女人水多毛片| 精品熟女少妇av免费看| 精品午夜福利在线看| 男插女下体视频免费在线播放| 国产乱人偷精品视频| 夫妻性生交免费视频一级片| 欧美日韩精品成人综合77777| 视频中文字幕在线观看| a级毛色黄片| 免费看不卡的av| 日韩一本色道免费dvd| 欧美精品国产亚洲| 国产精品嫩草影院av在线观看| 老女人水多毛片| 欧美日韩一区二区视频在线观看视频在线 | 男的添女的下面高潮视频| 啦啦啦在线观看免费高清www| 毛片女人毛片| 精品视频人人做人人爽| 男的添女的下面高潮视频| 国产成人91sexporn| 午夜老司机福利剧场| 高清在线视频一区二区三区| h日本视频在线播放| 精品国产乱码久久久久久小说| 久久人人爽av亚洲精品天堂 | 超碰97精品在线观看| 麻豆成人av视频| 国产精品无大码| 在线观看美女被高潮喷水网站| 91精品国产九色| 看免费成人av毛片| 亚洲熟女精品中文字幕| 日日摸夜夜添夜夜爱| 日日啪夜夜爽| 免费av不卡在线播放| 18禁裸乳无遮挡动漫免费视频 | 一个人看的www免费观看视频| 男女下面进入的视频免费午夜| 亚洲国产精品成人综合色| 成年av动漫网址| 人体艺术视频欧美日本| 日韩大片免费观看网站| 日韩亚洲欧美综合| 国产色婷婷99| 欧美精品人与动牲交sv欧美| 国产亚洲午夜精品一区二区久久 | 寂寞人妻少妇视频99o| 黑人高潮一二区| 色婷婷久久久亚洲欧美| 一级毛片aaaaaa免费看小| 国产乱人视频| 91久久精品国产一区二区三区| 亚洲成人一二三区av| www.av在线官网国产| 亚洲精品乱久久久久久| 最后的刺客免费高清国语| 真实男女啪啪啪动态图| 午夜福利在线观看免费完整高清在| 各种免费的搞黄视频| 成年人午夜在线观看视频| 97超碰精品成人国产| 国产精品福利在线免费观看| 国产高清有码在线观看视频| 国产乱人偷精品视频| 日本免费在线观看一区| 又粗又硬又长又爽又黄的视频| 简卡轻食公司| 亚洲人成网站在线观看播放| 国产高潮美女av| 各种免费的搞黄视频| 少妇 在线观看| 欧美日本视频| 国产真实伦视频高清在线观看| 黄色怎么调成土黄色| 一级黄片播放器| 久久久久久久久久成人| 少妇人妻 视频| 欧美成人一区二区免费高清观看| 国产一区有黄有色的免费视频| 我的老师免费观看完整版| 美女xxoo啪啪120秒动态图| 涩涩av久久男人的天堂| 欧美精品人与动牲交sv欧美| 一区二区三区四区激情视频| 成人综合一区亚洲| 日产精品乱码卡一卡2卡三| 最近中文字幕高清免费大全6| 亚洲美女搞黄在线观看| 亚洲欧美中文字幕日韩二区| 最近最新中文字幕大全电影3| 在线观看av片永久免费下载| 青青草视频在线视频观看| 一级av片app| 亚洲经典国产精华液单| 亚洲精品色激情综合| 亚洲欧美日韩另类电影网站 | 免费黄色在线免费观看| 色5月婷婷丁香| 亚洲av成人精品一二三区| 久久国内精品自在自线图片| 午夜免费男女啪啪视频观看| 久久国产乱子免费精品| 精品久久久精品久久久| 午夜免费观看性视频| 日韩欧美精品v在线| videossex国产| 久久久a久久爽久久v久久| 欧美精品人与动牲交sv欧美| 欧美性感艳星| 欧美成人午夜免费资源| 在线天堂最新版资源| 97在线视频观看| 国产爽快片一区二区三区| 身体一侧抽搐| 成人漫画全彩无遮挡| 成年av动漫网址| 三级国产精品欧美在线观看| 偷拍熟女少妇极品色| 亚洲精品一区蜜桃| 边亲边吃奶的免费视频| 免费大片18禁| 又爽又黄a免费视频| 亚洲一区二区三区欧美精品 | 伦精品一区二区三区| 免费av毛片视频| 国产中年淑女户外野战色| 国产欧美亚洲国产| 丰满人妻一区二区三区视频av| 人妻系列 视频| 九九在线视频观看精品| 久久久精品94久久精品| 如何舔出高潮| 日本免费在线观看一区| 国国产精品蜜臀av免费| 免费看av在线观看网站| 全区人妻精品视频| 久久ye,这里只有精品| 亚洲美女搞黄在线观看| 日韩人妻高清精品专区| 久久97久久精品| 99热这里只有是精品50| 中文欧美无线码| 婷婷色综合www| 国产一区二区三区综合在线观看 | 国产在线男女| 男女啪啪激烈高潮av片| 亚洲精品久久午夜乱码| 最近最新中文字幕大全电影3| 亚洲欧美一区二区三区黑人 | 亚洲成人中文字幕在线播放| 亚洲高清免费不卡视频| 国产淫语在线视频| 国产精品偷伦视频观看了| 最后的刺客免费高清国语| 特级一级黄色大片| av又黄又爽大尺度在线免费看| 国产精品一二三区在线看| 亚洲成色77777| av一本久久久久| 久久精品国产a三级三级三级| 最近中文字幕2019免费版| 中国国产av一级| 国产一区有黄有色的免费视频| 亚洲精品视频女| 一个人观看的视频www高清免费观看| 男人添女人高潮全过程视频| 久久精品国产自在天天线| 国产中年淑女户外野战色| 欧美日韩综合久久久久久| 69人妻影院| 波多野结衣巨乳人妻| 日日撸夜夜添| 干丝袜人妻中文字幕| 欧美最新免费一区二区三区| 国国产精品蜜臀av免费| 亚洲图色成人| 国产精品国产三级专区第一集| 六月丁香七月| 免费av观看视频| 国产亚洲91精品色在线| 成人鲁丝片一二三区免费| 国产免费一级a男人的天堂| 91精品伊人久久大香线蕉| 国产又色又爽无遮挡免| 免费av毛片视频| 建设人人有责人人尽责人人享有的 | 一二三四中文在线观看免费高清| 精品一区二区三卡| 亚洲欧洲国产日韩| 国产黄色视频一区二区在线观看| 久久久久国产网址| 欧美少妇被猛烈插入视频| 久久精品国产亚洲av天美| 国产精品国产三级专区第一集| 日韩亚洲欧美综合| 波野结衣二区三区在线| 搡老乐熟女国产| 日日啪夜夜爽| 亚洲熟女精品中文字幕| 久久久久网色| 亚洲成人av在线免费| 一级二级三级毛片免费看| 亚洲三级黄色毛片| 又大又黄又爽视频免费| 亚洲精品成人av观看孕妇| 99视频精品全部免费 在线| 在线观看国产h片| 成人美女网站在线观看视频| 黄片wwwwww| 国产中年淑女户外野战色| 肉色欧美久久久久久久蜜桃 | 日韩不卡一区二区三区视频在线| 亚洲精品亚洲一区二区| 国产黄片视频在线免费观看| 可以在线观看毛片的网站| 免费av毛片视频| 91精品国产九色| 免费高清在线观看视频在线观看| 毛片女人毛片| 99久久九九国产精品国产免费| 日韩av免费高清视频| 久久久久久久久久久丰满| 神马国产精品三级电影在线观看| 欧美日韩视频高清一区二区三区二| 久久久久久久亚洲中文字幕| 久久久久精品久久久久真实原创| 国产中年淑女户外野战色| 1000部很黄的大片| 能在线免费看毛片的网站| 国产毛片a区久久久久| 国产精品女同一区二区软件| 国产大屁股一区二区在线视频| 久久精品国产自在天天线| 三级国产精品片| 91久久精品电影网| 亚洲四区av| 三级男女做爰猛烈吃奶摸视频| 久久99热这里只频精品6学生| 少妇的逼水好多| 久久人人爽人人爽人人片va| 97在线人人人人妻| 国内揄拍国产精品人妻在线| 国产中年淑女户外野战色| 交换朋友夫妻互换小说| 男插女下体视频免费在线播放| h日本视频在线播放| 久久久国产一区二区| 亚洲国产成人一精品久久久| 只有这里有精品99| 久久午夜福利片| 亚洲av中文字字幕乱码综合| 国产一区亚洲一区在线观看| 久久久久久久久久久丰满| 最新中文字幕久久久久| 18禁动态无遮挡网站| 亚洲欧美精品专区久久| 国产中年淑女户外野战色| 边亲边吃奶的免费视频| 99久久精品国产国产毛片| 国产毛片a区久久久久| 在线观看国产h片| a级一级毛片免费在线观看| 国产成人精品婷婷| 免费黄频网站在线观看国产| 亚洲伊人久久精品综合| 欧美高清性xxxxhd video| 精品人妻视频免费看| 秋霞在线观看毛片| 边亲边吃奶的免费视频| 日本欧美国产在线视频| 亚洲国产高清在线一区二区三| 国产精品久久久久久av不卡| 日韩,欧美,国产一区二区三区| 丰满乱子伦码专区| 日韩电影二区| 国产精品人妻久久久影院| 观看免费一级毛片| 亚洲av男天堂| 直男gayav资源| 国内揄拍国产精品人妻在线| 精品99又大又爽又粗少妇毛片| 国产一区二区在线观看日韩| 婷婷色av中文字幕| 少妇人妻一区二区三区视频| 91精品国产九色| 成人黄色视频免费在线看| 国产精品久久久久久久久免| 国产精品无大码| 亚洲精品456在线播放app| 亚洲成人中文字幕在线播放| 男人狂女人下面高潮的视频| 日本色播在线视频| 国产精品久久久久久久久免| 亚洲熟女精品中文字幕| 嘟嘟电影网在线观看| 欧美日韩视频精品一区| 热re99久久精品国产66热6| 国产精品成人在线| 成人鲁丝片一二三区免费| 91aial.com中文字幕在线观看| 一本一本综合久久| 麻豆久久精品国产亚洲av| 国产精品久久久久久精品电影小说 | 女的被弄到高潮叫床怎么办| 亚洲精品久久久久久婷婷小说| 免费黄色在线免费观看| 日产精品乱码卡一卡2卡三| 少妇熟女欧美另类| 亚洲一区二区三区欧美精品 | 亚洲精品久久久久久婷婷小说| 久久久久性生活片| 中文天堂在线官网| 一二三四中文在线观看免费高清| 五月开心婷婷网| 亚洲欧美日韩卡通动漫| 少妇高潮的动态图| av卡一久久| 在线天堂最新版资源| 男人和女人高潮做爰伦理| 中文天堂在线官网| 亚洲精品成人av观看孕妇| 2021天堂中文幕一二区在线观| 人妻系列 视频| 免费观看a级毛片全部| 在线亚洲精品国产二区图片欧美 | 国产亚洲91精品色在线| 天堂网av新在线| 少妇的逼好多水| 午夜免费男女啪啪视频观看| 我要看日韩黄色一级片| 亚洲性久久影院| 91在线精品国自产拍蜜月| 少妇高潮的动态图| 亚洲国产最新在线播放| 国产有黄有色有爽视频| 老司机影院成人| 久久久午夜欧美精品| 亚洲精品456在线播放app| 免费播放大片免费观看视频在线观看| 精品99又大又爽又粗少妇毛片| 99久久精品热视频| 中文在线观看免费www的网站| 亚洲精品乱久久久久久| 免费人成在线观看视频色| av一本久久久久| 内射极品少妇av片p| 好男人视频免费观看在线| 夫妻性生交免费视频一级片| 大香蕉97超碰在线| 1000部很黄的大片| 麻豆乱淫一区二区| 日韩中字成人| 高清在线视频一区二区三区| 七月丁香在线播放| 国产免费又黄又爽又色| a级毛片免费高清观看在线播放| 欧美人与善性xxx| 六月丁香七月| 国产午夜福利久久久久久| 亚洲国产色片| 在线看a的网站| 亚洲精品日本国产第一区| 91午夜精品亚洲一区二区三区| 亚洲成人一二三区av| 亚洲av成人精品一二三区| 丝袜喷水一区| 欧美少妇被猛烈插入视频| 免费看光身美女| 91精品伊人久久大香线蕉| 在线a可以看的网站| 麻豆成人午夜福利视频| 国产精品国产三级专区第一集| 国产精品久久久久久精品古装|