999精品在线视频,手机成人午夜在线视频,久久不卡国产精品无码,中日无码在线观看,成人av手机在线观看,日韩精品亚洲一区中文字幕,亚洲av无码人妻,四虎国产在线观看 ?

四列角接觸球軸承振動(dòng)特性研究

2017-05-17 02:11:08鄧四二盧羽佳華顯偉張言偉黃運(yùn)生
振動(dòng)與沖擊 2017年9期
關(guān)鍵詞:振動(dòng)

鄧四二, 盧羽佳, 華顯偉, 張言偉, 黃運(yùn)生

(1.河南科技大學(xué) 機(jī)電工程學(xué)院,河南 洛陽(yáng) 471003; 2.遼寧重大裝備制造協(xié)同創(chuàng)新中心,遼寧 大連 116024)

四列角接觸球軸承振動(dòng)特性研究

鄧四二1,2, 盧羽佳1, 華顯偉1, 張言偉1, 黃運(yùn)生1

(1.河南科技大學(xué) 機(jī)電工程學(xué)院,河南 洛陽(yáng) 471003; 2.遼寧重大裝備制造協(xié)同創(chuàng)新中心,遼寧 大連 116024)

基于滾動(dòng)軸承動(dòng)力學(xué)理論,建立了四列角接觸球軸承的動(dòng)力學(xué)分析模型,并以某型號(hào)四列角接觸球軸承為例,對(duì)不同結(jié)構(gòu)參數(shù)與工況參數(shù)下的軸承振動(dòng)特性進(jìn)行理論分析。結(jié)果表明:對(duì)于四列角接觸球軸承,根據(jù)使用工況選擇大、小球列不同的內(nèi)、外溝曲率半徑系數(shù)與初始接觸角更有利于軸承的減振降噪;對(duì)軸承施加一定的軸向預(yù)緊量可有效減小軸承振動(dòng);存在較為合理的軸向載荷、傾覆力矩及內(nèi)圈轉(zhuǎn)速范圍使四列角接觸球軸承在使用時(shí)的振動(dòng)較小。

四列角接觸球軸承;動(dòng)力學(xué);結(jié)構(gòu)參數(shù);振動(dòng)

傳統(tǒng)的輪轂軸承單元多采用雙列角接觸球軸承與雙列圓錐滾子軸承,新型的四列角接觸球軸承作為輪轂軸承單元使用具有承載能力強(qiáng)、摩擦功耗小,制造成本低等優(yōu)勢(shì),可以替代雙列圓錐滾子軸承用于中型及重型汽車中[1]。作為輪轂軸承,四列角接觸球軸承的振動(dòng)特性將直接影響到汽車的安全性與舒適性。

軸承振動(dòng)的產(chǎn)生機(jī)理非常復(fù)雜,涉及工藝誤差、潤(rùn)滑劑、結(jié)構(gòu)參數(shù)及工況參數(shù)等諸多因素[2]。解決軸承振動(dòng)的基本方向?yàn)楦纳戚S承的制造工藝,通過控制軸承在加工過程中產(chǎn)生的圓度、粗糙度及波紋度等誤差實(shí)現(xiàn)軸承減振降噪的目的[3-7]。隨著制造工藝的發(fā)展,當(dāng)軸承精度超過某一界限時(shí),軸承的振動(dòng)水平不再由其制造精度決定,即使不存在幾何誤差,軸承仍會(huì)產(chǎn)生振動(dòng)。目前有關(guān)軸承結(jié)構(gòu)參數(shù)對(duì)軸承振動(dòng)特性影響的研究主要集中于單列球軸承的鋼球個(gè)數(shù)及徑向游隙等:Aktürk等[8-9]分別建立了三自由度與二自由度的角接觸球軸承支承的轉(zhuǎn)子系統(tǒng)振動(dòng)數(shù)學(xué)模型,對(duì)鋼球數(shù)、預(yù)載荷與主軸的振動(dòng)關(guān)系進(jìn)行了研究;Bai等[10-11]分別建立了角接觸球軸承與雙列調(diào)心球軸承的動(dòng)力學(xué)分析模型,開展了軸承徑向游隙、徑向力對(duì)軸承振動(dòng)特性的影響研究;殷玉楓等[12-13]分別通過聲學(xué)理論與滾動(dòng)軸承動(dòng)力學(xué)理論研究了深溝球軸承的溝曲率半徑系數(shù)、徑向游隙等結(jié)構(gòu)參數(shù)對(duì)軸承振動(dòng)噪聲的影響。汽車輪轂軸承在工作時(shí)的受力情況比較復(fù)雜,四列角接觸球軸承作為新型的汽車輪轂軸承單元結(jié)構(gòu),目前有關(guān)其振動(dòng)特性理論的研究尚未見。鑒于此,本文基于滾動(dòng)軸承動(dòng)力學(xué)分析理論,建立了四列角接觸球軸承的動(dòng)力學(xué)分析模型,以某型號(hào)四列角接觸球軸承為例,對(duì)軸承的工況和結(jié)構(gòu)等參數(shù)與軸承振動(dòng)特性的關(guān)系進(jìn)行分析,研究結(jié)果可為四列角接觸球軸承的設(shè)計(jì)提供一定的理論依據(jù)。

1 軸承動(dòng)力學(xué)分析模型

四列角接觸球軸承受載情況見圖1,為便于對(duì)四列角接觸球軸承進(jìn)行動(dòng)力學(xué)分析,作如下假設(shè):軸承零件工作表面具有理想的幾何形狀,其形心與質(zhì)心重合;軸承外圈固定,兩內(nèi)圈通過鎖扣連結(jié)后以轉(zhuǎn)速n繞X軸旋轉(zhuǎn);保持架由鋼球引導(dǎo);軸承承受徑向載荷、軸向載荷及力矩載荷;坐標(biāo)系{O;X,Y,Z}為軸承慣性坐標(biāo)系。

(a)

(b)

1.1 軸承各元件間相互作用力數(shù)學(xué)模型

圖1中,設(shè)左側(cè)大鋼球所在的列為第1列、左側(cè)小鋼球所在的列為第2列、右側(cè)小鋼球所在的列為第3列、右側(cè)大鋼球所在的列為第4列,則第k列第j個(gè)鋼球的受力情況見圖2;位于第k列的保持架受力情況見圖3;軸承內(nèi)圈的受力情況見圖4。圖2~圖4中相關(guān)符號(hào)的含義見表1,其值求解法參見文獻(xiàn)[14]。

(a)左側(cè)鋼球x-y平面(b)左側(cè)鋼球y-z平面(c)左側(cè)鋼球x-z平面(d)右側(cè)鋼球x-y平面(e)右側(cè)鋼球y-z平面(f)右側(cè)鋼球x-z平面

圖2 鋼球受力示意圖

Fig.2Schematicdiagramofballforces

1.2 軸承動(dòng)力學(xué)微分方程

根據(jù)軸承各零件的受力情況可建立四列角接觸球軸承各零件的動(dòng)力學(xué)微分方程組。

表1 符號(hào)說明

圖3 保持架受力示意圖

圖4 內(nèi)圈受力示意圖

鋼球的動(dòng)力學(xué)微分方程組為

(1)

(2)

(3)

(4)

(5)

(6)

保持架動(dòng)力學(xué)微分方程組為

(7)

(8)

(9)

內(nèi)圈動(dòng)力學(xué)微分方程組為

(10)

(11)

(13)

(14)

2 軸承振動(dòng)特性分析

采用GSTIFF(Gear Stiff)變步長(zhǎng)積分算法對(duì)建立的四列角接觸球軸承動(dòng)力學(xué)微分方程式(1)~式(14)進(jìn)行求解,以某型號(hào)四列角接觸球軸承為例,對(duì)軸承的工況和結(jié)構(gòu)等參數(shù)與軸承振動(dòng)特性的關(guān)系進(jìn)行分析。軸承主參數(shù)如表2所示,內(nèi)、外圈及滾動(dòng)體的材料為GCr15,保持架的材料為尼龍材料PA66。

表2 四列角接觸球軸承結(jié)構(gòu)參數(shù)

本文采用軸承振動(dòng)加速度級(jí)法[15]得到的軸承運(yùn)轉(zhuǎn)分貝值來評(píng)價(jià)軸承的振動(dòng),分貝值定義為

DB=20lg(a/a0)

(15)

式中:a為某一頻帶范圍內(nèi)的軸承振動(dòng)加速度均方根值;a0為參考加速度,a0=9.81×10-3m/s2。

2.1 結(jié)構(gòu)參數(shù)對(duì)軸承振動(dòng)的影響

設(shè)軸承在環(huán)境溫度下工作,分析軸承在內(nèi)圈轉(zhuǎn)速為900 r/min,徑向載荷為5 900 N,軸向載荷為2 400 N,力矩載荷為644.325 N·m,軸承軸向預(yù)緊量為10 mm時(shí),不同結(jié)構(gòu)參數(shù)對(duì)軸承振動(dòng)特性的影響。

2.1.1 內(nèi)溝曲率半徑系數(shù)對(duì)軸承振動(dòng)的影響

圖5為四列角接觸球軸承的大、小球列在不同初始接觸角α下,軸承的內(nèi)溝曲率半徑系數(shù)對(duì)其振動(dòng)的影響。由圖5可知,對(duì)于軸承的大、小球列,隨著內(nèi)溝曲率半徑系數(shù)以及初始接觸角的增加,軸承振動(dòng)均呈現(xiàn)先減小后增加的趨勢(shì),因此,對(duì)于給定工況條件,存在最優(yōu)的內(nèi)溝曲率半徑系數(shù)與初始接觸角,使四列角接觸球軸承的振動(dòng)最小。由圖5還可知,軸承振動(dòng)最小時(shí),大、小球列對(duì)應(yīng)的最優(yōu)內(nèi)溝曲率半徑系數(shù)與最優(yōu)初始接觸角均不同,因此,在設(shè)計(jì)四列角接觸球軸承時(shí),可以考慮針對(duì)軸承大、小球列選擇不同的內(nèi)溝曲率半徑系數(shù)與初始接觸角,以最大限度的降低軸承振動(dòng)。

圖5 內(nèi)溝曲率半徑系數(shù)對(duì)軸承振動(dòng)的影響

Fig.5 Effect of inner groove curvature radius coefficients on vibration value

2.1.2 外溝曲率半徑系數(shù)對(duì)軸承振動(dòng)的影響

圖6為四列角接觸球軸承的大、小球列在不同初始接觸角下,軸承的外溝曲率半徑系數(shù)對(duì)其振動(dòng)的影響。由圖6可知,軸承振動(dòng)隨外溝曲率半徑系數(shù)與初始接觸角的增加呈先減小后增加的趨勢(shì),因此,對(duì)于軸承的大、小球列,均存在一個(gè)使軸承振動(dòng)最小的外溝曲率半徑系數(shù)與初始接觸角。與內(nèi)溝曲率半徑系數(shù)類似,在給定工況下,大、小球列在軸承振動(dòng)最小時(shí)對(duì)應(yīng)的外溝曲率半徑系數(shù)與初始接觸角均不同,因此,在設(shè)計(jì)四列角接觸球軸承時(shí),可以考慮對(duì)軸承大、小球列選擇不同的外溝曲率半徑系數(shù)與初始接觸角,以最大限度的降低軸承振動(dòng)。

圖6 外溝曲率半徑系數(shù)對(duì)軸承振動(dòng)的影響

2.1.3 保持架兜孔間隙對(duì)軸承振動(dòng)的影響

圖7為四列角接觸球軸承在左、右兩側(cè)的兩個(gè)保持架分別為整體式設(shè)計(jì)與分離式設(shè)計(jì)時(shí),保持架兜孔間隙對(duì)軸承振動(dòng)的影響。由圖7可知,隨著保持架兜孔間隙的增加,軸承振動(dòng)呈先減小后增加的趨勢(shì),因此,存在最優(yōu)的保持架兜孔間隙使軸承的振動(dòng)最小。由圖7還可知,對(duì)于保持架采用整體式設(shè)計(jì)的軸承,隨著保持架兜孔間隙的變化,其振動(dòng)值均大于保持架采用分離式設(shè)計(jì)的軸承,因此,分離式設(shè)計(jì)的保持架更有利于軸承的減振降噪。

圖7 保持架兜孔間隙對(duì)軸承振動(dòng)的影響

2.2 工況參數(shù)對(duì)軸承振動(dòng)的影響

設(shè)軸承在環(huán)境溫度下工作,分析軸承結(jié)構(gòu)參數(shù)為初始設(shè)計(jì)值時(shí),不同工況參數(shù)對(duì)軸承振動(dòng)特性的影響。

2.2.1 徑向載荷對(duì)軸承振動(dòng)的影響

圖8為四列角接觸球軸承在內(nèi)圈轉(zhuǎn)速為900 r/min,軸向載荷為0 N,力矩載荷為0 N·m,軸承軸向預(yù)緊量δ分別為5 mm,10 mm,15 mm,20 mm時(shí),軸承承受的徑向載荷對(duì)軸承振動(dòng)的影響。由圖8可知,隨著徑向載荷的增加,軸承振動(dòng)也逐漸增加;隨著軸向預(yù)緊量的增加,軸承振動(dòng)呈現(xiàn)先減小后增加的趨勢(shì),且在徑向載荷較小時(shí),軸承的最優(yōu)軸向預(yù)緊量較小,徑向載荷較大時(shí),軸承的最優(yōu)軸向預(yù)緊量也較大,因此,可以根據(jù)給定工況選擇合理的軸承軸向預(yù)緊量,以最大限度的降低軸承振動(dòng)。

圖8 徑向載荷對(duì)軸承振動(dòng)的影響

2.2.2 軸向載荷比對(duì)軸承振動(dòng)的影響

圖9為四列角接觸球軸承在內(nèi)圈轉(zhuǎn)速為900 r/min,徑向載荷為5 900 N,力矩載荷為0 N·m,軸承軸向預(yù)緊量為10 mm時(shí),軸承軸向載荷比對(duì)軸承振動(dòng)特性的影響。軸向載荷比Ca定義為軸承軸向載荷與徑向載荷之比:Ca=Fa/Fr。由圖9可知,隨著軸向載荷比的增加,軸承振動(dòng)呈先基本不變后逐漸增加的趨勢(shì)。因此,對(duì)給定結(jié)構(gòu)的軸承,存在一個(gè)合理的軸向載荷使用范圍,使軸承運(yùn)轉(zhuǎn)過程中的振動(dòng)較小。

圖9 軸向載荷比對(duì)軸承振動(dòng)的影響

2.2.3 傾覆力矩對(duì)軸承振動(dòng)的影響

圖10為四列角接觸球軸承在內(nèi)圈轉(zhuǎn)速為900 r/min,徑向載荷為5 900 N,軸向載荷為2 400 N,軸承軸向預(yù)緊量為10 mm時(shí),軸承承受的傾覆力矩對(duì)軸承振動(dòng)特性的影響。由圖10可知,隨著傾覆力矩的增加,軸承振動(dòng)呈先緩慢增加后迅速增加的趨勢(shì)。因此,對(duì)于給定結(jié)構(gòu)的軸承,存在一個(gè)較為合理的傾覆力矩使用范圍,使軸承在運(yùn)轉(zhuǎn)時(shí)的振動(dòng)較小。

圖10 傾覆力矩對(duì)軸承振動(dòng)的影響

2.2.4 內(nèi)圈轉(zhuǎn)速對(duì)軸承振動(dòng)的影響

圖11為四列角接觸球軸承在徑向載荷為5 900 N,軸向載荷為2 400 N,軸承軸向預(yù)緊量為10 mm,傾覆力矩為644.325 N·m時(shí),內(nèi)圈轉(zhuǎn)速對(duì)軸承振動(dòng)特性的影響。由圖11可知,隨著內(nèi)圈轉(zhuǎn)速的增加,軸承振動(dòng)呈先增加后減小而后繼續(xù)增加的趨勢(shì),因此,對(duì)于給定結(jié)構(gòu)的軸承,存在一個(gè)合理的轉(zhuǎn)速使用范圍。

圖11 內(nèi)圈轉(zhuǎn)速對(duì)軸承振動(dòng)的影響

3 結(jié) 論

(1) 對(duì)于四列角接觸球軸承,存在最優(yōu)的內(nèi)、外溝曲率半徑系數(shù)與初始接觸角使軸承的振動(dòng)最小;根據(jù)工況參數(shù)選擇軸承大、小球列不同的內(nèi)、外溝曲率半徑系數(shù)與初始接觸角更有利于軸承的減振降噪。

(2) 對(duì)四列角接觸球軸承施加一定的軸向預(yù)緊量可有效降低軸承的振動(dòng)值,但過大的軸向預(yù)緊量會(huì)導(dǎo)致軸承振動(dòng)迅速增大。

(3) 四列角接觸球軸承在運(yùn)轉(zhuǎn)過程中,存在較為合理的軸向載荷、傾覆力矩及內(nèi)圈轉(zhuǎn)速使用范圍使軸承的振動(dòng)較小。

致謝:感謝遼寧重大裝備制造協(xié)同創(chuàng)新中心支持。

[1] 謝夫勒兩合公司. 角接觸球軸承結(jié)構(gòu)的車輪軸承單元[P]. 德國(guó):80020034.1, 2006-12-20.

[2] 趙聯(lián)春. 球軸承的振動(dòng)研究[D]. 杭州:浙江大學(xué),2003.

[3] AKTüRK N. The effect of waviness on vibrations associated with ball bearings[J]. ASME Journal of Tribology, 1999, 121(4): 667-677.

[4] HARSHA S P, SANDEEP K, PRAKASH R. Nonlinear dynamic behaviors of rolling element bearings due to surface waviness[J]. Journal of Sound and Vibration, 2004, 272(3): 557-580.

[5] JANG G, JEONG S W. Vibration analysis of a rotating system due to the effect of ball bearing waviness[J]. Journal of Sound and Vibration, 2004, 269(3/4/5): 709-726.

[6] HARSHA S P. The effect of ball size variation on nonlinear vibration associated with ball bearings[J].Journal of Multi-body Dynamics, 2004, 218(4):191-210.

[7] 顧曉輝,楊紹普,劉永強(qiáng),等. 表面波紋度對(duì)滾動(dòng)軸承—轉(zhuǎn)子系統(tǒng)非線性振動(dòng)的影響[J]. 振動(dòng)與沖擊,2014,33(8):109-114.

GU Xiaohui, YANG Shaopu, LIU Yongqiang, et al. Effect of surface waveness on nonlinear virbtarion of a rotor with ball bearings[J]. Journal of Vibration and Shock, 2014, 33(8): 109-114.

[8] AKTüRK N, UNEEB M, GOHAR R. The effects of number of balls and preload on vibrations associated with ball bearings[J]. Journal of Tribology, 1997, 119(4): 747-753.

[9] HARSHA S P, SANDEEP K, PRAKASH R. Effects of preload and number of balls on nonlinear dynamic behavior of ball bearing system[J]. International Journal of Nonlinear Sciences and Numerical Simulation, 2003, 4(3): 265-278.

[10] BAI Changqing, XU Qingyu. Dynamic model of ball bearings with internal clearance and waviness[J]. Journal of Sound and Vibration, 2006, 294(1):23-48.

[11] ZHUO Yaobin, ZHOU Xiaojun, YANG Chenlong. Dynamic analysis of double-row self-aligning ball bearings due to applied loads, internal clearance, surface waviness and number of balls[J]. Journal of Sound and Vibration, 2014, 333(23):6170-6189.

[12] 殷玉楓,張建水. 結(jié)構(gòu)參數(shù)對(duì)軸承振動(dòng)噪聲的影響[J]. 噪聲與振動(dòng)控制,2014,34(1):76-81.

YIN Yufeng, ZHANG Jianshui.Effect of structural parameters on vibration noise of rolling bearings[J]. Noise and Vibration Control, 2014, 34(1):76-81.

[13] 鄧四二,孫朝陽(yáng),顧金芳,等. 低噪音深溝球軸承振動(dòng)特性研究[J]. 振動(dòng)與沖擊,2015,34(10):12-19.

DENG Sier, SUN Chaoyang, GU Jinfang, et al.Study on the vibration characteristics of low-noise deep groove ball beaings[J].Journal of Vibration and Shock, 2015, 34(10): 12-19.

[14] 鄧四二,郝建軍,滕弘飛,等. 角接觸球軸承保持架動(dòng)力學(xué)分析[J].軸承,2007(10):1-5.

DENG Sier, HAO Jianjun, TENG Hongfei, et al. Dynamic analysis on cage of angular contact ball bearings[J]. Bearing, 2007(10): 1-5.

[15] 滾動(dòng)軸承振動(dòng)(加速度)測(cè)量方法:JB/T 5314—1991[S]. 北京:機(jī)械科學(xué)研究院,1991.

Vibration characteristics of four-column angular contact ball bearings

DENG Sier1,2, LU Yujia1, HUA Xianwei1, ZHANG Yanwei1, HUANG Yunsheng1

(1. School of Machatronics Engineering, Henan University of Science and Technology, Luoyang 471003, China;2. Collaborative Innovation Center of Major Machine Manufacturing of Liaoning, Dalian 116024, China)

Based on the dynamic theory of rolling bearings, the dynamic model of four-column angular contact ball bearings was established. Taking a certain type of four-column angular contact ball bearing as an example, the vibration characteristics of the bearing were investigated under different structural parameters and working conditions. The results showed that designing different inner and outer raceway groove curvature radius coefficients and different initial contact angles of different columns, respectively leads to lower vibration based on a specific working condition; a certain axial preload can reduce bearing vibration effectively; there are a reasonable axial load, an overturning moment and a rotating speed range to make the bearing vibration smaller.

four-column angular contact ball bearing; dynamics; structural parameters; vibration

河南省科技創(chuàng)新杰出人才(144200510020);國(guó)家自然科學(xué)基金項(xiàng)目(U1404514);廣東省省級(jí)科技計(jì)劃項(xiàng)目(2013B011301020)

2015-12-01 修改稿收到日期:2016-03-18

鄧四二 男,博士,教授,博士生導(dǎo)師,1963年生

盧羽佳 女,碩士生,1991年生

TH133.33

A

10.13465/j.cnki.jvs.2017.09.001

猜你喜歡
振動(dòng)
振動(dòng)的思考
某調(diào)相機(jī)振動(dòng)異常診斷分析與處理
振動(dòng)與頻率
This “Singing Highway”plays music
具非線性中立項(xiàng)的廣義Emden-Fowler微分方程的振動(dòng)性
中立型Emden-Fowler微分方程的振動(dòng)性
基于ANSYS的高速艇艉軸架軸系振動(dòng)響應(yīng)分析
船海工程(2015年4期)2016-01-05 15:53:26
主回路泵致聲振動(dòng)分析
UF6振動(dòng)激發(fā)態(tài)分子的振動(dòng)-振動(dòng)馳豫
帶有強(qiáng)迫項(xiàng)的高階差分方程解的振動(dòng)性
主站蜘蛛池模板: 久久九九热视频| 蜜芽一区二区国产精品| 亚洲国产在一区二区三区| 精品久久综合1区2区3区激情| 91在线精品麻豆欧美在线| 日本欧美成人免费| 国产成人91精品免费网址在线| 91亚瑟视频| 国产清纯在线一区二区WWW| 国产日韩欧美视频| 国产欧美日韩精品综合在线| 麻豆国产在线观看一区二区| 无遮挡一级毛片呦女视频| 成人精品亚洲| 色综合天天操| 日韩欧美在线观看| 国产迷奸在线看| 亚洲资源在线视频| 99草精品视频| 久久久久人妻一区精品色奶水| 91久久国产综合精品女同我| 亚洲va精品中文字幕| 亚洲av无码专区久久蜜芽| 国产乱子伦视频在线播放| 亚洲精品无码不卡在线播放| 国产成人无码综合亚洲日韩不卡| 91美女视频在线| 久久这里只有精品国产99| 国产精品成人一区二区不卡| 亚洲妓女综合网995久久| 亚洲第一极品精品无码| 老色鬼久久亚洲AV综合| 农村乱人伦一区二区| 中文字幕在线观| 青青草一区| 久久香蕉国产线| 色男人的天堂久久综合| 欧美一区福利| 国产成人精品午夜视频'| 一级毛片在线免费视频| 欧美日韩理论| 五月激情综合网| 国产成在线观看免费视频| 久久夜色精品| 无码网站免费观看| 精品伊人久久久大香线蕉欧美| 日本一区高清| 亚洲第一页在线观看| 国产成人无码Av在线播放无广告| 国产精品无码一区二区桃花视频| 69av在线| 久久精品日日躁夜夜躁欧美| 91在线高清视频| 中文字幕亚洲综久久2021| 久草网视频在线| 亚洲成a人片77777在线播放| 国产亚洲高清在线精品99| 国产AV无码专区亚洲A∨毛片| 又污又黄又无遮挡网站| 日本高清在线看免费观看| 国产不卡一级毛片视频| 女人毛片a级大学毛片免费| 国产主播喷水| 99视频在线免费| 亚洲三级成人| 国产男女XX00免费观看| 国产精品lululu在线观看| 成人在线观看一区| 久久精品无码国产一区二区三区| 伊人久久婷婷| 亚洲男人天堂网址| 国产国语一级毛片在线视频| 色综合中文综合网| 国内精自线i品一区202| 91偷拍一区| 亚洲第一黄色网| 日韩在线播放中文字幕| 女人18毛片久久| 日韩国产另类| 精品视频第一页| 中文字幕欧美日韩高清| 亚洲 欧美 偷自乱 图片|