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

    滾動(dòng)軸承故障信號(hào)在復(fù)雜路徑中的傳遞特性分析

    2021-08-21 11:12:16陳霽恒朱丹宸譚經(jīng)松
    機(jī)電工程技術(shù) 2021年12期
    關(guān)鍵詞:滾動(dòng)軸承

    陳霽恒 朱丹宸 譚經(jīng)松

    摘要:復(fù)雜傳遞路徑的存在會(huì)降低滾動(dòng)軸承故障信號(hào)采集的準(zhǔn)確性和有效性。針對(duì)此問題,建立了包含局部缺陷滾動(dòng)軸承、彈性支撐、油膜減振器和多層螺栓法蘭連接的動(dòng)力學(xué)仿真模型,研究了滾動(dòng)軸承故障信號(hào)在復(fù)雜路徑中的傳遞特點(diǎn)。采用滾動(dòng)軸承故障模擬平臺(tái)進(jìn)行滾動(dòng)軸承外圈故障實(shí)驗(yàn),對(duì)比動(dòng)力學(xué)仿真數(shù)值計(jì)算和實(shí)驗(yàn)分析的結(jié)果,驗(yàn)證了仿真結(jié)果的準(zhǔn)確性。借助典型的特征參數(shù)值,對(duì)比了不同測(cè)點(diǎn)處信號(hào)的特點(diǎn),并重點(diǎn)分析了油膜減振器的影響。研究結(jié)果表明,復(fù)雜傳遞路徑會(huì)導(dǎo)致故障信號(hào)在傳遞過程中產(chǎn)生較大的衰減,油膜減振器的存在會(huì)進(jìn)一步加劇信號(hào)的衰減,不利于滾動(dòng)軸承故障信號(hào)的有效傳遞。

    關(guān)鍵詞:復(fù)雜傳遞路徑;滾動(dòng)軸承;油膜減振器;動(dòng)力學(xué)建模

    中圖分類號(hào):TH133.3???????????? 文獻(xiàn)標(biāo)志碼:A??????? 文章編號(hào):1009-9492(2021)12-0083-07

    Research on Fault Signal of Rolling Element Bearing Influenced by Complex Transmission Path

    Chen Jiheng ,Zhu Danchen ,Tan Jingsong

    (Naval Petty Officer Academy, Bengbu, Anhui 233012, China)

    Abstract: The influence of complex transmission path will reduce the effectiveness and accuracy of the fault diagnosis of rolling element bearing. In view of this problem, the dynamic simulation model was established with the defective bearing, the oil film damper, the elastic support and the bolt flange connection structure, and the transmission characteristics of fault signals of rolling element bearing in complex transmission paths were studied. The fault simulation experiment of rolling element bearing with outer race defect was carried out by using the fault simulation platform, the accuracy of the simulation results was verified by comparing with the results of experimental analysis. With the help of typical characteristic parameters, the characteristics of signals at different measuring points were compared, and the influence of oil film damper was analyzed emphatically. The results show that the complex transmission path will lead to greater attenuation of the fault signal and the existence of oil film damper will further aggravate the attenuation of fault signal, which are not conducive to the effective transmission of the fault signal.

    Key words: complex transmission path; rolling element bearing; oil film damper; dynamic modeling

    0 引言

    滾動(dòng)軸承作為旋轉(zhuǎn)機(jī)械的核心零部件,在各行業(yè)中都得到了廣泛運(yùn)用。由于實(shí)際設(shè)備結(jié)構(gòu)較為復(fù)雜,滾動(dòng)軸承一般位于設(shè)備內(nèi)部,同時(shí)受到測(cè)試條件的制約,用于監(jiān)測(cè)滾動(dòng)軸承運(yùn)行狀態(tài)的信號(hào)大多只能從設(shè)備外側(cè)的殼體上獲取,導(dǎo)致軸承故障信號(hào)在向外傳遞的過程中需經(jīng)過復(fù)雜的傳遞路徑,對(duì)故障信號(hào)的有效傳遞產(chǎn)生較大影響,因此,研究故障沖擊信號(hào)在復(fù)雜傳遞路徑中的傳遞特點(diǎn)對(duì)提高軸承故障診斷的準(zhǔn)確性有著很強(qiáng)的現(xiàn)實(shí)意義。

    綜合國內(nèi)外對(duì)于復(fù)雜傳遞路徑動(dòng)力學(xué)特性的研究能夠看出,分析過程基本是將仿真建模分析和實(shí)驗(yàn)測(cè)試研究相結(jié)合,仿真研究主要包括數(shù)學(xué)模型仿真研究、有限元模型仿真研究。王峰等[1]建立了考慮箱體內(nèi)潤(rùn)滑油流固耦合的齒輪箱有限元模型,研究了人字齒輪傳動(dòng)系統(tǒng)的振動(dòng)傳遞特性、合理預(yù)估了箱體的結(jié)構(gòu)振動(dòng)噪聲。武星宇等和劉文璽等[2-3]利用仿真分別研究了螺旋槳激振力、對(duì)潛艇結(jié)構(gòu)振動(dòng)和聲輻射特性的影響。針對(duì)各種車輛開展的振動(dòng)傳遞特性分析能夠模擬車輛運(yùn)行過程,發(fā)現(xiàn)導(dǎo)致車輛振動(dòng)和噪聲較大的原因,為車輛的結(jié)構(gòu)優(yōu)化設(shè)計(jì)提供了技術(shù)支持[4-6]。溫華兵等[7]對(duì)加筋圓柱殼體支撐結(jié)構(gòu)傳遞特征開展了試驗(yàn)研究,結(jié)果有助于控制動(dòng)力機(jī)械設(shè)備的中低頻振動(dòng)傳遞。這些研究的目的大都是通過傳遞特性的研究,了解振動(dòng)、能量等參數(shù)在材料及設(shè)備內(nèi)部的傳遞規(guī)律,為材料性能改進(jìn)、設(shè)備結(jié)構(gòu)優(yōu)化等提供支撐。

    近年來,包含滾動(dòng)軸承的復(fù)雜支撐和傳動(dòng)結(jié)構(gòu)的振動(dòng)特性也到了廣泛研究,內(nèi)容主要包括了滾動(dòng)軸承及相關(guān)結(jié)構(gòu)的建模仿真與實(shí)驗(yàn)分析。Garad等[8]建立了一種含內(nèi)圈和外圈故障的滾動(dòng)軸承模型,考慮了故障大小、位置等的影響,對(duì)軸承的振動(dòng)特性進(jìn)行了分析。Yu 等[9]建立了一種4自由度的6201深溝球軸承數(shù)學(xué)模型,對(duì)軸承的各位置的位移、速度和加速度進(jìn)行了仿真。劉成穎等[10]建立了主軸-軸承-機(jī)殼的力學(xué)模型,將整個(gè)主軸系統(tǒng)簡(jiǎn)化為質(zhì)量-彈簧-阻尼系統(tǒng),實(shí)現(xiàn)了軸承支撐結(jié)構(gòu)的傳遞特性分析。Patel等[11]研究了深溝球軸承結(jié)構(gòu)的振動(dòng)傳遞特性,考慮了多種故障形式,對(duì)軸、基座、滾珠處的振動(dòng)信號(hào)進(jìn)行了時(shí)域和頻域分析。Xiao 等[12]構(gòu)建了一個(gè)8自由度的齒輪-軸-軸承-基座系統(tǒng),研究了振動(dòng)從齒輪到基座的傳遞特性。

    由此,本文以實(shí)驗(yàn)室滾動(dòng)軸承故障模擬平臺(tái)的軸承支撐結(jié)構(gòu)為研究對(duì)象,模型建立時(shí)考慮了滾動(dòng)軸承、彈性鼠籠支撐、油膜減振器、多層螺栓法蘭連接等結(jié)構(gòu),通過合理的模型簡(jiǎn)化,基于數(shù)值仿真計(jì)算,對(duì)滾動(dòng)軸承故障沖擊信號(hào)在復(fù)雜傳遞路徑中的傳遞情況進(jìn)行研究;重點(diǎn)研究了油膜減振器的作用;結(jié)合仿真過程和結(jié)果,進(jìn)一步通過實(shí)驗(yàn)研究軸承故障信號(hào)在不同結(jié)構(gòu)間的傳遞情況,并驗(yàn)證仿真分析的結(jié)果。

    1 帶局部故障的滾動(dòng)軸承動(dòng)力學(xué)模型

    1.1 滾動(dòng)軸承的簡(jiǎn)化動(dòng)力學(xué)模型及假設(shè)

    本文基于 Hertz 接觸理論建立了如圖1所示的滾動(dòng)軸承4自由度動(dòng)力學(xué)模型。

    該模型存在 x、y 兩個(gè)運(yùn)動(dòng)方向,假設(shè)滾動(dòng)軸承內(nèi)圈和外圈之間通過非線性接觸剛度和黏性阻尼進(jìn)行連接,忽略了滾動(dòng)體質(zhì)量、轉(zhuǎn)動(dòng)慣量以及滾動(dòng)體打滑的影響。 1.2 滾動(dòng)體與滾道間非線性接觸剛度的計(jì)算

    滾動(dòng)軸承內(nèi)圈和外圈間的總變形是滾動(dòng)體與內(nèi)、外圈之間各自變形的和,由此,滾動(dòng)體與內(nèi)外圈之間的總接觸剛度 kio 可以利用赫茲接觸理論進(jìn)行計(jì)算[12]:

    式中:kin 為滾動(dòng)體與內(nèi)圈的接觸剛度;kout 為滾動(dòng)體與外圈的接觸剛度;Eeq 為等效彈性模量,Eeq=E/(1-ν2),假設(shè)滾動(dòng)體與內(nèi)外圈材料相同,E 為彈性模量,v 為泊松比;κ為橢圓偏心參數(shù);Γ為第一類橢圓積分;Σ為第二類橢圓積分。

    則κ、Γ和Σ可以表示為[12]:

    式(4) ~(6)中,Σρin 和Σρout 分別表示滾動(dòng)體與內(nèi)、外圈接觸的曲率和,且:

    式中: d 為滾動(dòng)體的直徑;γ=dcosα/D ,D 為軸承的節(jié)徑,α為接觸角;ri 和ro分別為內(nèi)滾道和外滾道的曲率半徑。

    1.3 時(shí)變位移激勵(lì)

    本文以滾動(dòng)軸承存在外圈故障為例進(jìn)行研究,假設(shè)缺陷的長(zhǎng)度遠(yuǎn)大于缺陷的寬度,滾動(dòng)體在經(jīng)過缺陷時(shí)的下降位移小于缺陷的深度,滾動(dòng)體經(jīng)過缺陷時(shí)的示意圖如圖2所示,圖中,H 和 B 分別表示缺陷的深度和寬度,Hc 為最大的位移激勵(lì),且 Hc 可以表示為:

    由此,滾動(dòng)體經(jīng)過整個(gè)缺陷時(shí)的位移變化可以用正弦函數(shù)表示為[13]:

    式中:θ為缺陷對(duì)應(yīng)的角度;θi 為第 i 個(gè)滾動(dòng)體在任意時(shí)刻與缺陷位置的角度差;θ0為0時(shí)刻各滾動(dòng)體與缺陷位置的角度差。

    θ、θi 和θ0可以表示為[13]:

    式中:Do 和 Di 分別表示滾動(dòng)軸承外滾道和內(nèi)滾道直徑;φo 和φi 分別為外圈故障和內(nèi)圈故障情況下,缺陷位置與第一個(gè)滾動(dòng)體的角度差;ωs 為軸的轉(zhuǎn)頻;ωc 為滾動(dòng)體的公轉(zhuǎn)角速度。

    ωc 可以表示為:

    2 滾動(dòng)軸承支撐結(jié)構(gòu)

    本文以實(shí)驗(yàn)室某型試驗(yàn)臺(tái)的滾動(dòng)軸承支撐結(jié)構(gòu)為對(duì)象建立圖3所示的動(dòng)力學(xué)仿真模型。由圖可知,軸承故障信號(hào)從結(jié)構(gòu)內(nèi)部傳輸?shù)綑C(jī)匣的過程中,會(huì)受到油膜減振器、螺栓法蘭連接等結(jié)構(gòu)的影響。因此,以該結(jié)構(gòu)為研究對(duì)象進(jìn)行分析可以體現(xiàn)出傳遞路徑對(duì)于信號(hào)的干擾。

    1.油膜減振器外支撐 2.彈性鼠籠支撐 3.試驗(yàn)軸承

    4.機(jī)匣 5.外法蘭 6.底座 7.內(nèi)法蘭 8.輻板圖3 滾動(dòng)軸承支撐結(jié)構(gòu)

    2.1 擠壓油膜阻尼器模型

    本節(jié)基于流體動(dòng)壓潤(rùn)滑理論進(jìn)行油膜力的計(jì)算,并在后續(xù)支撐結(jié)構(gòu)建模時(shí)考慮油膜力的作用。在圓柱坐標(biāo)系內(nèi),依據(jù)雷諾方程,忽略流體慣性力和離心力,油膜阻尼器的壓力分布可以表示為[14]:

    式中:R 為阻尼器的半徑;h 為局部油膜的厚度;p 為壓力分布;θ為最大油膜厚度處算起的角度;z 為沿軸頸的軸線方向;μ為滑油的動(dòng)力粘度,一般來說,動(dòng)力粘度值隨溫度升高而減小; φ? 為軸頸進(jìn)動(dòng)角速度。

    根據(jù)短軸承假設(shè)和 Reynolds邊界條件,油膜力在徑向和切向可分別表示為[14]:

    式中:L 為阻尼器的軸向長(zhǎng)度;c 阻尼器的徑向間隙;ε為軸頸的偏心率;I1、I2和 I3為 Sommerfeld積分。

    I1、I2和 I3可以表示為[14~15]:

    為了方便仿真計(jì)算,需要將油膜力在直角坐標(biāo)系內(nèi)進(jìn)行表示,油膜力在x 和y 方向的分量分別用 Fx 和 Fy 進(jìn)行表示[14]:

    在直角坐標(biāo)系內(nèi),推導(dǎo)可得[14]:

    2.2 支撐結(jié)構(gòu)的動(dòng)力學(xué)建模及運(yùn)動(dòng)方程

    為了建立圖3所示滾動(dòng)軸承支撐結(jié)構(gòu)的集中質(zhì)量模

    型,首先做出一些假設(shè):(1) 將軸承內(nèi)圈與主軸考慮為一個(gè)整體,假設(shè)軸承內(nèi)圈處的集中質(zhì)量是軸承內(nèi)圈與主軸的質(zhì)量和,同時(shí),將主軸視為剛體,忽略其運(yùn)動(dòng)中的變形;(2)將軸承外圈與彈性鼠籠支撐視為一個(gè)整體計(jì)算集中質(zhì)量,將軸承外圈和油膜減振器外支撐間的連接剛度近似等于鼠籠剛度,且兩側(cè)均受到非線性油膜力的作用;(3)將結(jié)構(gòu)中的外法蘭與機(jī)匣視為一個(gè)整體;將螺栓連接處接觸面間的作用考慮為線剛度和阻尼;(4)將各部件簡(jiǎn)化為集中質(zhì)量塊,只考慮質(zhì)量而忽略剛度和阻尼;系統(tǒng)僅在x-y 坐標(biāo)系范圍內(nèi)受力。

    基于上述假設(shè),以圖3所示的軸承支撐結(jié)構(gòu)為依據(jù),本文建立了滾動(dòng)軸承支撐結(jié)構(gòu)的動(dòng)力學(xué)仿真模型,具體如圖4所示。圖中,min 、mout 分別為軸承內(nèi)、外圈處的集中質(zhì)量;m1、m2和 mh 分別為油膜減振器外支撐、內(nèi)法蘭、機(jī)匣的質(zhì)量;kio 表示軸承內(nèi)、外圈間的非線性接觸剛度,k1o 為軸承外圈與油膜減振器外支撐間的支撐剛度,k12、k2h 為油膜減振器外支撐和內(nèi)法蘭以及內(nèi)法蘭和機(jī)匣之間兩處螺栓連接結(jié)構(gòu)的接觸剛度,kh 為基座處的剛度;cio 為軸承內(nèi)、外圈間的阻尼,c12、c2h 為螺栓接觸面間的接觸阻尼,ch 為基座處的阻尼;fsx 和fsy 分別是 x 和 y 方向的油膜力,作用在軸承外圈和油膜減震器外支撐上。

    基于圖4的軸承支撐結(jié)構(gòu)集中質(zhì)量模型,考慮軸承內(nèi)、外圈之間的非線性接觸剛度和非線性油膜力的作用,系統(tǒng)的運(yùn)動(dòng)方程如式(23)所示。

    式中:g 為重力加速度,?i 為第 i 個(gè)球位置。

    ?i = +ωc t +?0

    式中:?0為第一個(gè)球與x 軸的初始夾角。

    δi 為第 i 個(gè)滾動(dòng)體的形變量,可以表示為:

    δi =(xin -xout)cos ?i +(yin -yout)sin ?i -(cr +H)式中:cr 為徑向游隙。

    由于只有當(dāng)滾動(dòng)體和滾道接觸產(chǎn)生形變,即δi>0時(shí),才存在赫茲接觸力,因此,動(dòng)力學(xué)仿真時(shí),采用變量βi 對(duì)滾動(dòng)體與滾道是否接觸進(jìn)行判斷:

    3 數(shù)值仿真分析

    本文采用 Matlab內(nèi)置的定時(shí)間步長(zhǎng)“ode45”方法對(duì)構(gòu)建的軸承支撐結(jié)構(gòu)動(dòng)力學(xué)方程進(jìn)行求解。仿真計(jì)算時(shí),由于缺陷對(duì)應(yīng)的角度較小,為了保證計(jì)算的準(zhǔn)確性,將時(shí)間步設(shè)定為Δt=1/131072 s ,求解時(shí)的絕對(duì)誤差和相對(duì)誤差分別設(shè)為1×10-7和1×10-6。

    仿真計(jì)算中,僅在 x 方向?qū)S承內(nèi)圈施加預(yù)加載荷,即 Fx=200 N ,F(xiàn)y=0 N ;假設(shè)靜止?fàn)顩r下,油膜阻尼器不存在偏心;轉(zhuǎn)速設(shè)定為1500 r/min;假設(shè)滾動(dòng)軸承存在外圈故障且故障的寬度設(shè)為0.6 mm;支撐結(jié)構(gòu)各部分在x 和y 方向的初始位移和初始速度都設(shè)為0。系統(tǒng)中的阻尼系數(shù)分別設(shè)為 cio=ch=300 N ·s/m ,c12=c2h=200 N ·s/m;剛度分別設(shè)為 k1o=2×107 N/m ,k12=k2h=kh=3×107 N/m 。仿真時(shí)所需其他參數(shù)如表1所示,且本節(jié)展示的所有仿真結(jié)果均在x 方向。

    通過計(jì)算可得,軸承內(nèi)、外圈之間的非線性接觸剛度 kio=6.32×109 N/m3/2,通過計(jì)算,6010軸承在1500 r/min 時(shí)的外圈故障特征頻率為140 Hz 。仿真計(jì)算的結(jié)果如圖5所示,本文中的包絡(luò)譜均通過對(duì)時(shí)域信號(hào)直接進(jìn)行Hilbert變換和快速傅里葉變換獲得,所得頻譜包含整個(gè)頻段,避免了帶通濾波器中心頻率和帶寬選取對(duì)結(jié)果的影響。圖 5(a)所示為軸承外圈振動(dòng)信號(hào)的時(shí)域波形,信號(hào)中存在明顯的周期性自由衰減振動(dòng),且從圖5(b)所示的包絡(luò)譜中,能夠識(shí)別出外圈故障特征頻率141 Hz 及其2~7倍頻;從圖5(c)所示的外支撐振動(dòng)信號(hào)的時(shí)域波形中,也能夠識(shí)別出周期性的沖擊成分,且圖5(d)所示的包絡(luò)譜中,軸承外圈故障特征頻率及其2~7倍頻也能得到較好識(shí)別;圖 5(e)所示為內(nèi)法蘭處振動(dòng)信號(hào)的時(shí)域波形,圖中雖然能夠識(shí)別周期性的沖擊成分,但由于外部結(jié)構(gòu)相對(duì)質(zhì)量較大,同時(shí)受到多層結(jié)構(gòu)傳遞的影響,軸承外圈故障產(chǎn)生的前一個(gè)沖擊未能得到有效衰減就出現(xiàn)了下一個(gè)沖擊,這一特點(diǎn)在圖5(g)中體現(xiàn)的更為明顯,由于時(shí)域信號(hào)中的周期性沖擊不再明顯,導(dǎo)致了圖5(h)中故障特征頻率的5~7倍頻幾乎難以識(shí)別。總的來看,從支撐結(jié)構(gòu)各部位振動(dòng)信號(hào)的包絡(luò)譜中,都能夠準(zhǔn)確提取出滾動(dòng)軸承外圈故障特征頻率及其倍頻成分,且與理論值相一致,說明了本文構(gòu)建的滾動(dòng)軸承非線性動(dòng)力學(xué)模型的準(zhǔn)確性,同時(shí)受到多層結(jié)構(gòu)傳遞的影響,軸承故障產(chǎn)生的周期性沖擊信號(hào)在靠近滾動(dòng)軸承部位要更為明顯。

    為了更直觀地展現(xiàn)軸承支撐結(jié)構(gòu)各部位振動(dòng)信號(hào)的周期性沖擊特性與能量的變化,本文選取了一些常見的參數(shù)進(jìn)行進(jìn)一步研究,包括均方根值(C1)、方根幅值(C2)、峭度值(C3)和峰值因子(C4)。信號(hào)的均方根值和方根幅值一樣都可以反映出信號(hào)能量的大小,峭度指標(biāo)可以體現(xiàn)信號(hào)中沖擊成分的大小,而峰值因子能較好的反映出軸承故障引起的瞬時(shí)沖擊振動(dòng)。對(duì)于離散信號(hào)%(n),假設(shè)其平均值和標(biāo)準(zhǔn)差分別為μ和σ,則上述4個(gè)參數(shù)可以表示為:

    圖6所示為所選4個(gè)參數(shù)隨振動(dòng)測(cè)點(diǎn)位置不同的變化曲線,由于仿真建模時(shí),從軸承外圈到機(jī)匣處的質(zhì)量不斷增大,且未考慮結(jié)構(gòu)自身振動(dòng)以及故障沖擊引起的局部共振,導(dǎo)致振動(dòng)信號(hào)從內(nèi)到外傳遞過程中,信號(hào)的均方根值和方根幅值不斷減小,信號(hào)能量呈衰減趨勢(shì)。受到多層結(jié)構(gòu)和傳遞路徑的影響,外側(cè)結(jié)構(gòu)振動(dòng)信號(hào)中的沖擊成分得不到有效衰減,軸承故障導(dǎo)致的周期性沖擊不明顯,因此,峭度值和峰度系數(shù)也呈減小趨勢(shì),這與圖5中觀察到的現(xiàn)象是一致的。

    為了重點(diǎn)研究油膜阻尼器的減振性能,本文對(duì)無阻尼器情況下系統(tǒng)各部位的響應(yīng)進(jìn)行計(jì)算,以阻尼器兩端結(jié)構(gòu)為對(duì)象,分析了油膜對(duì)其振動(dòng)的影響。圖7所示為無油膜和有油膜情況下滾動(dòng)軸承外圈和油膜減振器外支撐振動(dòng)的時(shí)域波形,為了使對(duì)比更明顯,每幅圖僅給出3個(gè)振動(dòng)衰減周期長(zhǎng)度的信號(hào)。

    根據(jù)圖7(a)所示的對(duì)比可知,相比于無油膜情況,油膜減振器使軸承外圈振動(dòng)的幅值明顯減小,且有油膜情況下,軸承故障產(chǎn)生的沖擊信號(hào)衰減更快,周期性故障沖擊更明顯;如圖7(b)所示,油膜的存在使得故障沖擊的振幅減小,沖擊信號(hào)的衰減振動(dòng)更明顯。

    對(duì)于實(shí)際設(shè)備而言,僅有溫度的改變會(huì)影響油膜的動(dòng)力粘度值(一般來說隨著溫度的升高,油膜的動(dòng)力粘度值會(huì)下降),由此,有必要研究油膜動(dòng)力粘度對(duì)油膜減振性能的影響。下面以軸承外圈處的振動(dòng)信號(hào)為研究對(duì)象,利用信號(hào)振幅、均方根值、峭度值等指標(biāo),研究動(dòng)力黏度在5×10-4~5×10-2范圍內(nèi)的油膜減振效果,結(jié)果如圖8所示。由圖可知,無油膜情況下,軸承外圈振動(dòng)信號(hào)的幅值和均方根值最大,峭度最小,此時(shí)信號(hào)能量較大且沖擊特征不明顯,隨著動(dòng)力粘度值增大,外圈振動(dòng)信號(hào)的幅值、均方根值總體呈下降趨勢(shì),峭度值呈上升趨勢(shì),說明隨著動(dòng)力粘度值的增大,油膜阻尼器的減振效果不斷增強(qiáng),振動(dòng)信號(hào)的能量不斷衰減,同時(shí),軸承故障產(chǎn)生的沖擊成分能很快衰減完畢導(dǎo)致信號(hào)中的沖擊特征明顯,峭度值較大。

    綜合來看,油膜減振器的存在改變了系統(tǒng)的阻尼特性,導(dǎo)致軸承故障信號(hào)通過油膜減振器向外傳遞時(shí)振幅產(chǎn)生明顯的的減小,且設(shè)備運(yùn)行過程中,油膜溫度的改變也會(huì)影響減振效果。

    4 實(shí)驗(yàn)研究

    由于仿真過程對(duì)整個(gè)支撐結(jié)構(gòu)包括滾動(dòng)軸承進(jìn)行了一定的簡(jiǎn)化,難以反映出結(jié)構(gòu)的全部特征。由此,為了進(jìn)一步研究復(fù)雜傳遞路徑對(duì)沖擊信號(hào)的影響,并對(duì)仿真結(jié)果進(jìn)行驗(yàn)證,本節(jié)以實(shí)驗(yàn)室的滾動(dòng)軸承故障模擬平臺(tái)為研究對(duì)象,通過實(shí)驗(yàn)對(duì)復(fù)雜傳遞路徑的特性進(jìn)行更為深入的研究。

    實(shí)驗(yàn)研究所用滾動(dòng)軸承故障模擬平臺(tái)的結(jié)構(gòu)如圖9(a)所示,主要包含軸承支撐結(jié)構(gòu)、軸向和徑向加載裝置、主軸、電機(jī)以及潤(rùn)滑系統(tǒng)等。實(shí)驗(yàn)軸承的型號(hào)選取為6010,通過激光切割的方法分別在測(cè)試軸承外圈切一個(gè)寬0.2 mm 、深0.2 mm 與軸線平行的窄縫,用以模擬軸承外圈故障,軸承的具體參數(shù)如表1所示。實(shí)驗(yàn)時(shí),軸的轉(zhuǎn)速為fr =1500 r/min ,選取圖9中的測(cè)點(diǎn)1和2進(jìn)行故障信號(hào)的測(cè)量,信號(hào)測(cè)量方向?yàn)閺较?,采樣頻率為32768 Hz。

    圖10~11所示為測(cè)點(diǎn)1和2測(cè)得的振動(dòng)信號(hào)的時(shí)域波形和包絡(luò)譜,圖 10(a)中能識(shí)別出軸承故障產(chǎn)生的周期性沖擊成分,噪聲干擾較小,且從圖10(b)的包絡(luò)譜中能明顯識(shí)別出頻率為142 Hz 的譜線(與故障特征頻率的理論值相近)及其2~6倍頻,說明該軸承存在外圈故障。測(cè)點(diǎn)2離故障位置稍遠(yuǎn),且內(nèi)、外法蘭之間通過螺栓進(jìn)行連接,加之可能存在的局部共振現(xiàn)象,導(dǎo)致測(cè)點(diǎn)2的振動(dòng)信號(hào)信噪比降低,但通過包絡(luò)分析仍能有效提取出故障特征頻率,但特征頻率對(duì)應(yīng)的譜線幅值有了明顯下降,效果略差于測(cè)點(diǎn)1,進(jìn)一步說明了復(fù)雜傳遞路徑會(huì)影響滾動(dòng)軸承故障信號(hào)的有效傳遞。通過與前一節(jié)的仿真結(jié)果進(jìn)行比較,仿真和實(shí)測(cè)信號(hào)的包絡(luò)譜中特征頻率數(shù)值相近,進(jìn)一步驗(yàn)證了滾動(dòng)軸承動(dòng)力學(xué)仿真模型的準(zhǔn)確性。

    與仿真過程類似,下面選取測(cè)點(diǎn)1和測(cè)點(diǎn)2處的信號(hào),分別計(jì)算其均方根值、方根幅值、峭度值和峰度因子,結(jié)果如表2所示。通過比較,測(cè)點(diǎn)1處信號(hào)的均方根值、方根幅值、峭度值和峰度因子均較大,說明相比于測(cè)點(diǎn)2處的信號(hào),測(cè)點(diǎn)1處信號(hào)能量值較大且沖擊特征較強(qiáng),這與仿真得到的結(jié)果一致。

    5 結(jié)束語

    (1) 建立了某軸承支撐結(jié)構(gòu)的集中質(zhì)量仿真模型,重點(diǎn)研究了當(dāng)滾動(dòng)軸承存在局部故障時(shí)離故障點(diǎn)遠(yuǎn)、近不同測(cè)點(diǎn)處的信號(hào)特點(diǎn),結(jié)果表明,受到復(fù)雜傳遞路徑的影響,離故障軸承較遠(yuǎn)測(cè)點(diǎn)處的信號(hào)能量較低且沖擊特征較弱。

    (2) 重點(diǎn)分析了油膜減振器的影響,結(jié)果表明,油膜減振器的存在會(huì)改變信號(hào)的原有特性,加速?zèng)_擊信號(hào)的衰減,且當(dāng)油膜動(dòng)力黏度值在一定范圍內(nèi)變化時(shí),較大的動(dòng)力黏度值會(huì)使得信號(hào)能量降低。

    (3) 進(jìn)行了滾動(dòng)軸承故障模擬實(shí)驗(yàn),驗(yàn)證了仿真分析結(jié)果,進(jìn)一步表明復(fù)雜傳遞路徑會(huì)加劇噪聲的干擾,減弱故障沖擊成分,致使故障特征提取的難度增大。

    參考文獻(xiàn):

    [1]王峰,方宗德, 李聲晉,等.人字齒輪系統(tǒng)振動(dòng)傳遞分析優(yōu)化與試驗(yàn)驗(yàn)證[J].機(jī)械工程學(xué)報(bào),2015,51(1):34-42.

    [2]武星宇, 紀(jì)剛,周其斗,等.螺旋槳激振力傳遞模式下的艇體振動(dòng)和聲輻射分析[J].中國艦船研究,2016, 11(6):90-96.

    [3]劉文璽,周其斗,譚路,等.尾后軸承剛度對(duì)潛艇結(jié)構(gòu)聲輻射特性的影響[J].船舶力學(xué),2018,22(2):235-247.

    [4]張琪昌,代英達(dá), 于躍斌,等.基于振動(dòng)傳遞特性分析的快捷重載貨車空載運(yùn)行時(shí)的振動(dòng)控制[J].天津大學(xué)學(xué)報(bào)(自然科學(xué)與工程技術(shù)版),2018,51(9):903-911.

    [5]褚志剛, 熊敏,楊洋,等.車內(nèi)噪聲時(shí)域傳遞路徑分析[J].振動(dòng)與沖擊,2015,34(17):161-166.

    [6] Noh H M.Contribution analysis of interior noise and floor vibra? tion in high-speed trains by operational transfer path analysis[J]. Advances in Mechanical Engineering,2017,9(8):1-14.

    [7]溫華兵,左言言,夏兆旺,等.加筋圓柱殼體支撐結(jié)構(gòu)振動(dòng)傳遞特性試驗(yàn)研究[J].船舶力學(xué),2013, 17(7):785-792.

    [8] Garad A,Shinde V J.A theoretical model of deep groove ball bear? ing for predicting the effect of localized defects on vibrations[J]. International Journal of Mechanical Engineering and Technology, 2017,8(6):760-769.

    [9] Yu G,Su M,Xia W,et al.Vibration characteristics of deep groove ball? bearing? based? on 4-DOF? mathematical? model[J].Procedia Engineering,2017(174):808-814.

    [10]劉成穎,鄭烽,王立平.高速電主軸中主軸-機(jī)殼振動(dòng)傳遞力學(xué)模型[J].清華大學(xué)學(xué)報(bào)(自然科學(xué)版),2018,58(7):671-676.

    [11] Patel V N,Tandon N,Pandey R K.A dynamic model for vibration studies of deep groove ball bearings considering single and mul? tiple? defects? in? races[J].Journal? of? Tribology,2010, 132(4):041101.

    [12] Xiao? H,Zhou X,Liu J,et? al.Vibration transmission? and? energy dissipation through the gear-shaft-bearing-housing system sub? jected? to? impulse? force? on? gear[J].Measurement,2017(102):64-79.

    [13] Liu J,Shao Y,Lim T C.Impulse vibration transmissibility charac? teristics? in? the? presence? of localized? surface? defects? in? deep groove ball bearing systems[J].Proceedings of the Institution of Mechanical Engineers,Part K:Journal of Multi-body Dynamics, 2014,228(1):62-81.

    [14]聞邦椿.高等轉(zhuǎn)子動(dòng)力學(xué)[M].北京:機(jī)械工業(yè)出版社,2000.

    [15] Booker J F.A table of the journal-bearing integral[J].Journal of Basic Engineering,1965,87(2):533-535.

    [16] Huang Y,Lin J,Liu Z,et al.A modified scale-space guiding varia? tional mode decomposition for high-speed railway bearing fault diagnosis[J].Journal of Sound and Vibration,2019(444):216-234.

    第一作者簡(jiǎn)介:陳霽恒(1974-),男,安徽安慶人,碩士,副教授,研究領(lǐng)域?yàn)榕灤瑱C(jī)械維修,已發(fā)表論文10余篇。

    (編輯:王智圣)

    猜你喜歡
    滾動(dòng)軸承
    一種三點(diǎn)接觸球軸承性能試驗(yàn)臺(tái)的研制
    面向風(fēng)力發(fā)電機(jī)組齒輪箱滾動(dòng)軸承故障診斷
    車用滾動(dòng)軸承失效原因分析及使用維護(hù)
    淺析機(jī)械維修過程中的軸承維修
    電動(dòng)機(jī)滾動(dòng)軸承裝配工藝及其方法探討
    商情(2017年1期)2017-03-22 12:06:11
    如何延長(zhǎng)拖拉機(jī)滾動(dòng)軸承的使用壽命
    一種滾動(dòng)軸承故障特征提取方法研究
    滾動(dòng)軸承裝配方法及注意事項(xiàng)
    滾動(dòng)軸承安裝使用中的幾個(gè)問題
    基于小波—EMD共振解調(diào)方法在軸承中的應(yīng)用研究
    国产91精品成人一区二区三区| 亚洲精品亚洲一区二区| 亚洲av.av天堂| 天美传媒精品一区二区| 国产成人欧美在线观看| 小蜜桃在线观看免费完整版高清| 亚洲精品在线美女| 色综合亚洲欧美另类图片| 在线观看美女被高潮喷水网站 | 日韩欧美在线乱码| 91久久精品电影网| 日韩欧美三级三区| 亚洲精华国产精华精| 国产单亲对白刺激| 国产亚洲精品综合一区在线观看| 神马国产精品三级电影在线观看| 亚洲狠狠婷婷综合久久图片| bbb黄色大片| 丰满人妻熟妇乱又伦精品不卡| 99国产精品一区二区三区| 不卡一级毛片| 天堂√8在线中文| 成人无遮挡网站| 欧美另类亚洲清纯唯美| 99久久九九国产精品国产免费| 午夜福利18| 午夜福利免费观看在线| 老司机午夜十八禁免费视频| 亚洲午夜理论影院| 午夜影院日韩av| 精品人妻1区二区| 亚州av有码| 亚洲黑人精品在线| 蜜桃久久精品国产亚洲av| 国产av一区在线观看免费| 中文字幕av在线有码专区| 99久久九九国产精品国产免费| 在线观看免费视频日本深夜| 亚洲av.av天堂| 亚洲五月天丁香| 久久草成人影院| netflix在线观看网站| 在线观看av片永久免费下载| 麻豆国产97在线/欧美| 午夜精品在线福利| 欧美一级a爱片免费观看看| 久久香蕉精品热| 久久亚洲精品不卡| 精品久久久久久久末码| 天堂影院成人在线观看| 日韩中文字幕欧美一区二区| 深夜精品福利| 日韩亚洲欧美综合| 欧美潮喷喷水| 日韩精品青青久久久久久| 国产精品三级大全| aaaaa片日本免费| 中文字幕高清在线视频| 亚洲精品亚洲一区二区| 观看美女的网站| 欧美丝袜亚洲另类 | 亚洲成人久久性| 97超视频在线观看视频| 成年人黄色毛片网站| 观看免费一级毛片| 国产极品精品免费视频能看的| 国产精品免费一区二区三区在线| 久久性视频一级片| 成人欧美大片| 精品人妻1区二区| 黄色视频,在线免费观看| 热99re8久久精品国产| 最好的美女福利视频网| 免费看美女性在线毛片视频| 99热这里只有是精品在线观看 | 日本 欧美在线| 亚洲av电影不卡..在线观看| 亚洲专区国产一区二区| 深夜精品福利| 国产中年淑女户外野战色| 国产在线男女| 色在线成人网| 久久亚洲真实| 国产精品永久免费网站| 国内精品美女久久久久久| 亚洲自偷自拍三级| 三级毛片av免费| 欧美高清成人免费视频www| 午夜福利在线在线| 国产伦人伦偷精品视频| 色综合欧美亚洲国产小说| 日本熟妇午夜| 我的女老师完整版在线观看| 欧美乱妇无乱码| 亚洲在线自拍视频| 麻豆久久精品国产亚洲av| 国产av不卡久久| 国产淫片久久久久久久久 | 久久精品国产清高在天天线| 国产人妻一区二区三区在| 国产伦精品一区二区三区四那| 在线观看av片永久免费下载| 亚洲欧美日韩卡通动漫| 国产精品一区二区免费欧美| 国产爱豆传媒在线观看| 我要看日韩黄色一级片| 亚洲av电影不卡..在线观看| 精品国产亚洲在线| 亚洲内射少妇av| 成年人黄色毛片网站| 国产精品电影一区二区三区| 精品久久久久久,| 欧美zozozo另类| 国产精品永久免费网站| 国产乱人视频| 国内毛片毛片毛片毛片毛片| 亚洲第一欧美日韩一区二区三区| 熟女电影av网| 成年女人永久免费观看视频| 免费看日本二区| 黄色一级大片看看| 波多野结衣高清作品| 亚洲国产色片| 国产午夜精品论理片| 无遮挡黄片免费观看| 91在线精品国自产拍蜜月| 欧美另类亚洲清纯唯美| 亚洲一区二区三区色噜噜| 色哟哟·www| 老司机午夜十八禁免费视频| 在线观看66精品国产| 丰满人妻一区二区三区视频av| 日韩亚洲欧美综合| 看黄色毛片网站| 亚洲成人免费电影在线观看| 午夜a级毛片| 久久精品综合一区二区三区| 给我免费播放毛片高清在线观看| 精品99又大又爽又粗少妇毛片 | 国产午夜精品论理片| 亚洲中文字幕日韩| 国产精品av视频在线免费观看| 天堂√8在线中文| 亚洲欧美日韩卡通动漫| 首页视频小说图片口味搜索| 高潮久久久久久久久久久不卡| 国内精品一区二区在线观看| 91av网一区二区| 久久人人爽人人爽人人片va | 尤物成人国产欧美一区二区三区| 黄色配什么色好看| av在线观看视频网站免费| 91麻豆精品激情在线观看国产| 欧美性猛交╳xxx乱大交人| 午夜精品在线福利| 久久久久国产精品人妻aⅴ院| 丰满的人妻完整版| 免费在线观看成人毛片| 国产白丝娇喘喷水9色精品| 欧美黄色片欧美黄色片| 国产精品久久久久久久久免 | 日本a在线网址| 级片在线观看| 免费看a级黄色片| 91久久精品国产一区二区成人| 老师上课跳d突然被开到最大视频 久久午夜综合久久蜜桃 | 老鸭窝网址在线观看| 久久久精品大字幕| 国产一区二区在线观看日韩| bbb黄色大片| 亚洲国产精品合色在线| 美女高潮喷水抽搐中文字幕| 亚洲精品456在线播放app | 美女高潮喷水抽搐中文字幕| 99热这里只有是精品在线观看 | 国产人妻一区二区三区在| 亚洲,欧美,日韩| 高清日韩中文字幕在线| 久久久国产成人精品二区| 在线a可以看的网站| 九色成人免费人妻av| 国产单亲对白刺激| 精品人妻偷拍中文字幕| 最后的刺客免费高清国语| 亚洲三级黄色毛片| 搡女人真爽免费视频火全软件 | 欧美精品国产亚洲| 亚洲一区高清亚洲精品| 免费在线观看影片大全网站| 国产黄色小视频在线观看| 国产aⅴ精品一区二区三区波| 久久热精品热| 亚洲五月天丁香| 国产激情偷乱视频一区二区| 亚洲综合色惰| 一级黄片播放器| av福利片在线观看| 一级毛片久久久久久久久女| 可以在线观看毛片的网站| 成人三级黄色视频| 麻豆国产av国片精品| 天天一区二区日本电影三级| 午夜福利免费观看在线| 国产大屁股一区二区在线视频| 久久热精品热| 男插女下体视频免费在线播放| 在线观看舔阴道视频| 亚洲激情在线av| 一级黄色大片毛片| 听说在线观看完整版免费高清| 日本在线视频免费播放| 亚洲国产精品sss在线观看| 亚洲最大成人手机在线| 99久久九九国产精品国产免费| 变态另类丝袜制服| 国产av一区在线观看免费| 国产精品久久久久久亚洲av鲁大| 97热精品久久久久久| 午夜日韩欧美国产| 欧美黑人欧美精品刺激| 国产精品久久久久久久久免 | 中国美女看黄片| 在线国产一区二区在线| 深夜精品福利| 尤物成人国产欧美一区二区三区| www日本黄色视频网| av女优亚洲男人天堂| 欧美+日韩+精品| 精华霜和精华液先用哪个| 亚洲av一区综合| 中文字幕精品亚洲无线码一区| 又爽又黄无遮挡网站| 国产亚洲精品av在线| 免费看a级黄色片| 亚洲五月婷婷丁香| 激情在线观看视频在线高清| 欧美性感艳星| 啦啦啦韩国在线观看视频| 在线观看一区二区三区| 午夜激情欧美在线| 99久久久亚洲精品蜜臀av| 亚洲aⅴ乱码一区二区在线播放| 国产精品电影一区二区三区| 亚洲av免费在线观看| 中出人妻视频一区二区| 人人妻,人人澡人人爽秒播| 国产成人影院久久av| 在线观看舔阴道视频| 久久久久久久久久黄片| 精品人妻一区二区三区麻豆 | 国产伦一二天堂av在线观看| 丰满乱子伦码专区| 脱女人内裤的视频| 亚洲,欧美精品.| 中出人妻视频一区二区| 久久人人精品亚洲av| 亚洲最大成人手机在线| 亚洲av二区三区四区| 乱人视频在线观看| 99热精品在线国产| 日本精品一区二区三区蜜桃| 国产av一区在线观看免费| 精品久久久久久,| 成人性生交大片免费视频hd| 精品不卡国产一区二区三区| 久久久久久久精品吃奶| 老司机深夜福利视频在线观看| 国产国拍精品亚洲av在线观看| 亚洲avbb在线观看| 精华霜和精华液先用哪个| 一进一出抽搐动态| 国产免费一级a男人的天堂| 亚洲五月天丁香| 国产精品99久久久久久久久| 黄色视频,在线免费观看| 人妻丰满熟妇av一区二区三区| 一区福利在线观看| 久久草成人影院| 国产三级中文精品| 亚洲人与动物交配视频| 亚洲精品在线美女| 18美女黄网站色大片免费观看| 久久久久久久久久成人| 色综合站精品国产| 国产av麻豆久久久久久久| 美女黄网站色视频| 天堂影院成人在线观看| 少妇丰满av| 夜夜夜夜夜久久久久| 亚洲va日本ⅴa欧美va伊人久久| 搡女人真爽免费视频火全软件 | 亚洲精品成人久久久久久| 午夜激情福利司机影院| 91午夜精品亚洲一区二区三区 | 99热精品在线国产| 国产精品一区二区三区四区久久| 老师上课跳d突然被开到最大视频 久久午夜综合久久蜜桃 | 午夜激情欧美在线| 91字幕亚洲| 男人的好看免费观看在线视频| 欧美不卡视频在线免费观看| 久久6这里有精品| 欧美日韩国产亚洲二区| 欧美一区二区亚洲| 91在线观看av| 亚洲狠狠婷婷综合久久图片| 中文字幕高清在线视频| 久久久久久久午夜电影| 赤兔流量卡办理| 精品午夜福利在线看| 日本在线视频免费播放| 久久99热这里只有精品18| 免费看日本二区| 国产av在哪里看| 国产人妻一区二区三区在| 亚洲国产高清在线一区二区三| 亚洲黑人精品在线| 亚洲av五月六月丁香网| 岛国在线免费视频观看| 又黄又爽又免费观看的视频| 岛国在线免费视频观看| 少妇高潮的动态图| 色播亚洲综合网| 国产成+人综合+亚洲专区| 九色成人免费人妻av| 亚洲在线观看片| 国产精品一及| 久久国产精品人妻蜜桃| 美女高潮喷水抽搐中文字幕| 午夜久久久久精精品| av欧美777| 亚洲精品456在线播放app | 日本a在线网址| 一本综合久久免费| 嫩草影院新地址| 好看av亚洲va欧美ⅴa在| 亚洲无线在线观看| 欧美高清成人免费视频www| 精品久久国产蜜桃| 中文字幕高清在线视频| 嫩草影院精品99| 丁香六月欧美| 国产精品日韩av在线免费观看| 天堂网av新在线| 亚洲av免费在线观看| 久久国产精品人妻蜜桃| www.999成人在线观看| 国产大屁股一区二区在线视频| 精品一区二区免费观看| 亚洲天堂国产精品一区在线| 亚洲av免费在线观看| 国产高清三级在线| 欧美黑人欧美精品刺激| 成人一区二区视频在线观看| 国产三级黄色录像| 国产亚洲精品av在线| 亚洲最大成人av| 少妇的逼水好多| 丰满乱子伦码专区| 91狼人影院| 国内毛片毛片毛片毛片毛片| 亚洲va日本ⅴa欧美va伊人久久| 欧美日韩乱码在线| 黄色一级大片看看| 午夜福利在线观看吧| av福利片在线观看| 青草久久国产| 成人毛片a级毛片在线播放| 久久午夜亚洲精品久久| 精品国产亚洲在线| 精品熟女少妇八av免费久了| 国产三级中文精品| 首页视频小说图片口味搜索| 3wmmmm亚洲av在线观看| 国产精品永久免费网站| 亚州av有码| 午夜精品一区二区三区免费看| 国产精品一区二区三区四区久久| 国产激情偷乱视频一区二区| 欧美性猛交╳xxx乱大交人| 亚洲无线观看免费| 一夜夜www| 日韩精品中文字幕看吧| 国产欧美日韩精品一区二区| 变态另类丝袜制服| 在线国产一区二区在线| 欧美日韩乱码在线| 国产乱人伦免费视频| 天天躁日日操中文字幕| 亚洲国产高清在线一区二区三| 免费观看的影片在线观看| 亚洲 欧美 日韩 在线 免费| 国产真实乱freesex| 国内少妇人妻偷人精品xxx网站| 熟女人妻精品中文字幕| 国产三级在线视频| 欧美另类亚洲清纯唯美| 宅男免费午夜| 啦啦啦韩国在线观看视频| 久久久久久国产a免费观看| 色综合婷婷激情| 最近视频中文字幕2019在线8| 99精品在免费线老司机午夜| 国产亚洲欧美在线一区二区| 国产午夜精品论理片| 51午夜福利影视在线观看| 99久久九九国产精品国产免费| 夜夜爽天天搞| 欧美日韩综合久久久久久 | 色哟哟·www| 午夜免费男女啪啪视频观看 | 免费av毛片视频| 校园春色视频在线观看| 亚洲精品成人久久久久久| 国产单亲对白刺激| 国产午夜精品久久久久久一区二区三区 | 国产伦在线观看视频一区| 嫩草影院精品99| 欧美日韩国产亚洲二区| 男女视频在线观看网站免费| 日韩中文字幕欧美一区二区| 熟女人妻精品中文字幕| 精品人妻视频免费看| 国产精品伦人一区二区| 他把我摸到了高潮在线观看| 天堂网av新在线| 国产精品三级大全| 亚洲,欧美,日韩| 嫩草影院新地址| 老熟妇乱子伦视频在线观看| 免费观看的影片在线观看| 成人av一区二区三区在线看| 免费在线观看成人毛片| 久久精品久久久久久噜噜老黄 | 99久久九九国产精品国产免费| 欧美xxxx黑人xx丫x性爽| 成人特级黄色片久久久久久久| 99在线视频只有这里精品首页| 男人舔奶头视频| 性欧美人与动物交配| 一a级毛片在线观看| www.熟女人妻精品国产| 久久人人爽人人爽人人片va | 长腿黑丝高跟| 国产亚洲精品综合一区在线观看| 久久6这里有精品| 九九在线视频观看精品| 简卡轻食公司| 亚洲欧美日韩无卡精品| 色尼玛亚洲综合影院| 亚洲精品色激情综合| 国产在线男女| 女人十人毛片免费观看3o分钟| 一级黄色大片毛片| 国产淫片久久久久久久久 | 99在线视频只有这里精品首页| 亚洲五月天丁香| 少妇高潮的动态图| 成年人黄色毛片网站| 97碰自拍视频| 亚洲av电影不卡..在线观看| 亚洲18禁久久av| 9191精品国产免费久久| 啪啪无遮挡十八禁网站| 国产极品精品免费视频能看的| 国产成人av教育| 18禁黄网站禁片免费观看直播| 高清在线国产一区| 女生性感内裤真人,穿戴方法视频| 国产午夜精品论理片| 老鸭窝网址在线观看| 国产一区二区三区视频了| 欧美xxxx性猛交bbbb| 婷婷色综合大香蕉| 久久久久久久亚洲中文字幕 | 国产精品乱码一区二三区的特点| 男人的好看免费观看在线视频| 亚洲 国产 在线| 91麻豆av在线| 亚洲,欧美精品.| 在线国产一区二区在线| 天天躁日日操中文字幕| 老熟妇仑乱视频hdxx| 亚洲精品456在线播放app | 中国美女看黄片| 老熟妇仑乱视频hdxx| 国产精品伦人一区二区| 高清日韩中文字幕在线| 免费av观看视频| 毛片一级片免费看久久久久 | 精品久久久久久久久亚洲 | 国产免费一级a男人的天堂| 色视频www国产| 久久婷婷人人爽人人干人人爱| 天堂网av新在线| 免费大片18禁| 全区人妻精品视频| 很黄的视频免费| 最新中文字幕久久久久| 窝窝影院91人妻| 日日干狠狠操夜夜爽| 亚洲精品在线美女| 首页视频小说图片口味搜索| 日韩欧美三级三区| 一本综合久久免费| 青草久久国产| 亚洲欧美日韩无卡精品| 精品一区二区三区av网在线观看| 制服丝袜大香蕉在线| 成人特级黄色片久久久久久久| aaaaa片日本免费| 久久久成人免费电影| 麻豆一二三区av精品| 国产精品久久久久久人妻精品电影| 国产精品一区二区性色av| 桃色一区二区三区在线观看| 免费电影在线观看免费观看| 亚洲国产日韩欧美精品在线观看| 中国美女看黄片| 国产精品久久久久久亚洲av鲁大| 国内久久婷婷六月综合欲色啪| 精品久久国产蜜桃| 国产黄片美女视频| 热99re8久久精品国产| 男女做爰动态图高潮gif福利片| 欧美最黄视频在线播放免费| 一区二区三区四区激情视频 | 热99re8久久精品国产| 日本 av在线| 欧美最黄视频在线播放免费| 女人十人毛片免费观看3o分钟| 1000部很黄的大片| 精品99又大又爽又粗少妇毛片 | 在线免费观看不下载黄p国产 | 18美女黄网站色大片免费观看| 嫩草影视91久久| 久久伊人香网站| 亚洲av日韩精品久久久久久密| 黄色一级大片看看| 国产真实乱freesex| 美女免费视频网站| 精品国内亚洲2022精品成人| 亚洲黑人精品在线| 国产免费一级a男人的天堂| 久久精品国产99精品国产亚洲性色| 深夜精品福利| 人妻制服诱惑在线中文字幕| 国产欧美日韩一区二区精品| 村上凉子中文字幕在线| 级片在线观看| 黄色视频,在线免费观看| 一级黄色大片毛片| 夜夜夜夜夜久久久久| 亚洲av免费在线观看| 国产探花极品一区二区| 国产淫片久久久久久久久 | 99热这里只有是精品在线观看 | 91九色精品人成在线观看| 嫁个100分男人电影在线观看| 精品久久久久久,| 啦啦啦观看免费观看视频高清| а√天堂www在线а√下载| 每晚都被弄得嗷嗷叫到高潮| 午夜福利成人在线免费观看| 麻豆国产97在线/欧美| 黄色日韩在线| 欧美色欧美亚洲另类二区| 亚州av有码| 亚洲国产精品合色在线| 无人区码免费观看不卡| 亚洲av.av天堂| 熟妇人妻久久中文字幕3abv| 在现免费观看毛片| 国产一区二区在线观看日韩| 少妇的逼水好多| 美女黄网站色视频| 日韩国内少妇激情av| 精品熟女少妇八av免费久了| 色综合亚洲欧美另类图片| 亚洲一区二区三区不卡视频| 麻豆国产av国片精品| 国产高清激情床上av| 91九色精品人成在线观看| 亚洲av免费高清在线观看| 午夜a级毛片| 欧美日本视频| 午夜精品在线福利| 久久精品夜夜夜夜夜久久蜜豆| 亚洲av中文字字幕乱码综合| 国产成人影院久久av| 搡老岳熟女国产| 午夜视频国产福利| 在线观看午夜福利视频| 亚洲国产精品sss在线观看| 日日摸夜夜添夜夜添小说| 色精品久久人妻99蜜桃| 村上凉子中文字幕在线| 日日夜夜操网爽| 最后的刺客免费高清国语| 亚洲精品日韩av片在线观看| 窝窝影院91人妻| 日韩高清综合在线| 五月伊人婷婷丁香| 午夜日韩欧美国产| 久久午夜亚洲精品久久| 亚洲美女视频黄频| 中文资源天堂在线| 成人性生交大片免费视频hd| 他把我摸到了高潮在线观看| 日本熟妇午夜| 一级黄片播放器| 亚洲av成人av| 黄色日韩在线| 久久草成人影院| or卡值多少钱|