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

    動(dòng)載荷作用下中介軸承的非線性熱行為研究1)

    2021-03-24 06:11:48高朋侯磊陳予恕
    力學(xué)學(xué)報(bào) 2021年1期
    關(guān)鍵詞:雙穩(wěn)態(tài)游隙滾子

    高朋 侯磊 陳予恕

    (哈爾濱工業(yè)大學(xué)航天學(xué)院,哈爾濱 150001)

    引言

    隨著轉(zhuǎn)子系統(tǒng)向高速重載方向發(fā)展,對(duì)支承系統(tǒng)的動(dòng)力學(xué)性能和熱性能提出了更高的要求.中介軸承作為雙轉(zhuǎn)子系統(tǒng)中高、低壓轉(zhuǎn)子之間的重要支承和傳動(dòng)部件,其動(dòng)力學(xué)特性[1]和熱特性[2]相較于普通支承軸承更加復(fù)雜,甚至許多情況下是非線性的[3].然而,目前關(guān)于中介軸承非線性熱行為的研究幾乎處于空白狀態(tài).因此,對(duì)中介軸承非線性熱行為機(jī)理進(jìn)行深入地研究,闡明動(dòng)力學(xué)參數(shù)以及熱參數(shù)對(duì)其的影響是非常有必要的.

    迄今為止,已經(jīng)有許多學(xué)者研究了軸承的Hertz接觸分?jǐn)?shù)指數(shù)非線性以及徑向游隙等非線性因素對(duì)轉(zhuǎn)子系統(tǒng)動(dòng)力學(xué)特性的影響.Yamamoto[4]研究了軸承徑向游隙對(duì)軸承?轉(zhuǎn)子系統(tǒng)共振特性的影響,發(fā)現(xiàn)共振峰的頻率和振幅隨徑向游隙的增大而減小.Fukata 等[5]研究了球軸承在恒定徑向載荷作用下的徑向振動(dòng),發(fā)現(xiàn)軸承的Hertz 接觸分?jǐn)?shù)指數(shù)非線性以及徑向游隙等非線性因素會(huì)引起超諧共振和亞諧共振等非線性共振.Holmes 等[3]在試驗(yàn)臺(tái)上復(fù)現(xiàn)了某中型噴氣發(fā)動(dòng)機(jī)中雙轉(zhuǎn)子結(jié)構(gòu)的跳躍現(xiàn)象、亞諧共振和組合共振等非線性動(dòng)力學(xué)行為.Mevel 和Guyader[6]基于Fukata 的動(dòng)力學(xué)模型研究了輕載球軸承通往混沌的兩種不同路徑,即次諧波路徑和準(zhǔn)周期路徑.Tiwari 等[7-8]應(yīng)用改進(jìn)的諧波平衡法,從理論上模擬了軸承Hertz 接觸和徑向游隙對(duì)水平轉(zhuǎn)子非線性動(dòng)力學(xué)行為的影響,并進(jìn)行了實(shí)驗(yàn)驗(yàn)證.Ghafari等[9]提出了考慮滾動(dòng)體非線性剛度的集中質(zhì)量?阻尼?彈簧模型,研究了徑向游隙對(duì)軸承平衡點(diǎn)的影響.Bai 等[10]通過(guò)數(shù)值仿真和實(shí)驗(yàn)驗(yàn)證,建立了考慮軸承非線性因素六自由度動(dòng)力學(xué)模型,研究了對(duì)稱球軸承?轉(zhuǎn)子系統(tǒng)的亞諧共振.Zhang 等[11-12]基于文獻(xiàn)[7-8]中的分析方法HB-AFT,重點(diǎn)研究了軸承的Hertz 接觸分?jǐn)?shù)指數(shù)非線性以及徑向游隙等非線性因素對(duì)轉(zhuǎn)子系統(tǒng)的共振滯后特性的影響.然而以上研究工作中,都集中在軸承?轉(zhuǎn)子系統(tǒng)的非線性動(dòng)力學(xué)行為上,沒(méi)有考慮軸承的熱效應(yīng).

    摩擦?xí)璧K物體運(yùn)動(dòng)并以摩擦熱的形式造成能量損失.中介軸承運(yùn)行中產(chǎn)生的摩擦熱可使軸承溫度升高,可采用摩擦力矩來(lái)度量.摩擦力矩的大小很大程度上取決于潤(rùn)滑方式,對(duì)于采用油潤(rùn)滑的圓柱滾子軸承,潤(rùn)滑油將占據(jù)軸承內(nèi)部空間,從而阻礙了滾子的運(yùn)動(dòng)[13].摩擦力矩的大小與潤(rùn)滑劑的性能、填充量以及滾子運(yùn)轉(zhuǎn)速度都有關(guān)系,然而,這些因素通常是緊密耦合在一起的,很難區(qū)分開(kāi)來(lái)[14].1945 年,Palmgren 等[15]通過(guò)對(duì)各種類型和尺寸的滾動(dòng)軸承進(jìn)行大量試驗(yàn),得出了軸承摩擦力矩的經(jīng)驗(yàn)公式,Palmgren 經(jīng)驗(yàn)公式作為一種預(yù)測(cè)滾動(dòng)軸承摩擦力矩的精確方法已被廣泛認(rèn)可和使用.根據(jù)熱傳遞的基礎(chǔ)理論,Harris 等[13]建立了滾動(dòng)軸承穩(wěn)態(tài)熱傳遞能量守恒方程,理論預(yù)測(cè)了滾動(dòng)軸承主要零部件滾珠、內(nèi)圈和外圈的穩(wěn)態(tài)溫度.Winer 等[16]則更進(jìn)一步建立了圓錐滾子軸承和軸承座間的熱傳遞模型,同時(shí)給出了轉(zhuǎn)軸、軸承以及軸承座之間的熱阻計(jì)算公式.DeMul等[17-18]提出了一個(gè)五自由度軸承受力模型,同時(shí)采用矩陣法分析載荷與撓度之間的關(guān)系.Jorgensen 和Shin[19]提出了一個(gè)準(zhǔn)三維傳熱模型來(lái)預(yù)測(cè)考慮熱膨脹的主軸?軸承系統(tǒng)的穩(wěn)態(tài)溫度分布.基于Palmgren的經(jīng)驗(yàn)公式,Stein 和Tu[14]提出了一個(gè)狀態(tài)空間模型,用于預(yù)測(cè)由角接觸球軸承的熱膨脹而引起的預(yù)緊力,并分析了轉(zhuǎn)速和初始預(yù)緊力的影響.Sun 等[20]開(kāi)發(fā)了一種葉片損耗模擬方法,并建立了一個(gè)傳熱模型來(lái)估計(jì)葉片損耗以及軸承主要零部件的熱膨脹數(shù)據(jù).Takabi 和Khonsari[21]提出了油浴潤(rùn)滑的深溝球軸承非穩(wěn)態(tài)熱傳遞模型,從而用來(lái)研究軸承的瞬態(tài)溫度.艾思源等[22]聚焦于脂潤(rùn)滑的雙列圓錐滾子軸承的熱行為,并提出了一種準(zhǔn)靜態(tài)模型,用以計(jì)算軸承載荷分布和運(yùn)動(dòng)參數(shù).Than 和Huang[23]提出了一種將準(zhǔn)靜態(tài)模型與有限元相結(jié)合的方法,研究了高速主軸軸承在預(yù)載荷作用下的非線性熱行為.然而,上述文獻(xiàn)均未考慮系統(tǒng)動(dòng)力學(xué)特性對(duì)軸承熱行為的影響.

    本文提出一種將傳熱學(xué)與非線性動(dòng)力學(xué)相結(jié)合的理論方法,研究雙轉(zhuǎn)子系統(tǒng)的非線性動(dòng)力學(xué)特性對(duì)中介軸承熱行為的影響機(jī)制,并詳細(xì)分析動(dòng)力學(xué)參數(shù)和熱參數(shù)的影響規(guī)律.根據(jù)雙轉(zhuǎn)子系統(tǒng)的非線性動(dòng)力響應(yīng)定義出中介軸承的動(dòng)載荷,借助Palmgren 經(jīng)驗(yàn)公式將其代入軸承熱平衡方程,從而建立動(dòng)載荷作用下的中介軸承非穩(wěn)態(tài)熱傳遞模型.通過(guò)數(shù)值模擬研究動(dòng)載荷作用下的中介軸承熱行為,發(fā)現(xiàn)軸承的運(yùn)行溫度不僅與轉(zhuǎn)子系統(tǒng)的運(yùn)行轉(zhuǎn)速相關(guān),也與轉(zhuǎn)子系統(tǒng)的動(dòng)力學(xué)特性息息相關(guān),這也為軸承的設(shè)計(jì)提供了一種新的思路.

    1 中介軸承動(dòng)載荷的非線性行為

    簡(jiǎn)化雙盤(pán)四支承的雙轉(zhuǎn)子系統(tǒng)如圖1 所示,每個(gè)轉(zhuǎn)子有一個(gè)圓盤(pán)和一個(gè)轉(zhuǎn)軸組成,圓盤(pán)固結(jié)在轉(zhuǎn)軸上.中介軸承位于高、低壓轉(zhuǎn)子之間,其不同于通常用于支承的軸承,中介軸承外圈也會(huì)隨著高壓轉(zhuǎn)子旋轉(zhuǎn).其中l(wèi)i(i=1,2,3,4,5) 為軸的長(zhǎng)度,ki和ci(i=1,2,3) 表示彈性支撐的剛度和阻尼系數(shù),ω1和ω2(rad/s) 是低壓轉(zhuǎn)子和高壓轉(zhuǎn)子轉(zhuǎn)速.定義轉(zhuǎn)速比為λ=ω2/ω1,由于高壓轉(zhuǎn)子轉(zhuǎn)速大于低壓轉(zhuǎn)子,故λ>1 為雙轉(zhuǎn)子同向旋轉(zhuǎn),λ <?1 為雙轉(zhuǎn)子反向旋轉(zhuǎn).

    圖1 簡(jiǎn)化雙盤(pán)四支承的雙轉(zhuǎn)子系統(tǒng)Fig.1 Simple dual-rotor system with two disks and four supports

    簡(jiǎn)化雙盤(pán)四支承的雙轉(zhuǎn)子系統(tǒng)動(dòng)力學(xué)方程可采用Lagrange 第二方程推導(dǎo),具體推導(dǎo)過(guò)程見(jiàn)文獻(xiàn)[24-25],系統(tǒng)八自由度動(dòng)力學(xué)方程如下

    方程動(dòng)力學(xué)參數(shù)物理意義及數(shù)值見(jiàn)文獻(xiàn)[24-25].

    對(duì)于航空發(fā)動(dòng)機(jī)雙轉(zhuǎn)子系統(tǒng),中介軸承一般采用徑向圓柱滾子軸承.中介軸承具有間隙、Hertz 接觸引起的分?jǐn)?shù)指數(shù)以及變剛度[26]等非線性因素,會(huì)對(duì)中介軸承的彈性恢復(fù)力產(chǎn)生不可忽視的影響.中介軸承運(yùn)動(dòng)示意圖如圖2 所示.

    第k個(gè)滾子的瞬時(shí)角位置可以表示為

    其中保持架轉(zhuǎn)速為ωc=(ω1ri+ω2ro)/(ri+ro),式中ri,ro分別表示中介軸承內(nèi)、外圈半徑.

    根據(jù)幾何關(guān)系,第k個(gè)滾子與滾道間的相對(duì)變形為

    其中2δ0為中介軸承的徑向游隙.

    圖2 中介軸承運(yùn)動(dòng)示意圖Fig.2 Motion diagram of inter-shaft bearing

    豎直方向和水平方向上的中介軸承恢復(fù)力為

    其中H(·)為Heaviside 函數(shù),Kb和Nb表示中介軸承的Hertz 接觸剛度和滾子數(shù)目.

    中介軸承采用德國(guó)FAG?公司生產(chǎn)的型號(hào)為NU1020 圓柱滾子軸承,其參數(shù)列于表1.

    表1 NU1020 圓柱滾子軸承參數(shù)Table 1 Parameters of the radial cylindrical roller bearing NU1020

    中介軸承動(dòng)載荷定義為一個(gè)周期內(nèi)豎直方向和水平方向上的中介軸承彈性恢復(fù)力的有效值[27],即方均根(root-mean-square,RMS),定義式如下

    其中,T是周期,N表示周期內(nèi)的離散化點(diǎn)數(shù),和為豎直方向和水平方向彈性恢復(fù)力周期內(nèi)的平均值.

    2 中介軸承非穩(wěn)態(tài)熱傳遞方程

    2.1 中介軸承摩擦熱

    中介軸承一般采用圓柱滾子軸承,總摩擦力矩M分為3 個(gè)部分:載荷摩擦力矩Ml、黏度摩擦力矩Mν和滾子端面?擋邊摩擦力矩Mf.中介軸承不承受軸向載荷,因此可以忽略Mf,則總摩擦力矩M為

    摩擦力矩單位統(tǒng)一采用N·mm.

    本文引用中介軸承動(dòng)載荷Fb來(lái)計(jì)算,而非以往文獻(xiàn)采用的靜載荷[23],則載荷摩擦力矩Ml為

    式中,fl是一個(gè)與軸承類型和載荷相關(guān)的系數(shù),對(duì)于帶保持架的向心圓柱滾子軸承,fl=0.000 2~0.000 4,小值用于輕系列軸承,大值用于重系列軸承[13];Dm是中介軸承節(jié)圓直徑.

    黏度摩擦力矩Mν顯然與潤(rùn)滑劑黏度ν 有關(guān),采用運(yùn)動(dòng)黏度(cSt,即mm2/s)來(lái)表示.中介軸承內(nèi)、外圈均隨低、高壓轉(zhuǎn)子同步運(yùn)轉(zhuǎn),因此采用內(nèi)外圈轉(zhuǎn)速差

    替代原公式中的轉(zhuǎn)速,單位為r/min.則黏度摩擦力矩Mν為

    式中fν是一個(gè)與軸承類型和潤(rùn)滑方式相關(guān)的系數(shù),帶保持架的向心圓柱滾子軸承,脂潤(rùn)滑fν=0.6~1,油霧潤(rùn)滑fν=1.5~2.8,油浴潤(rùn)滑fν=2.2~4,小值用于輕系列軸承,大值用于重系列軸承[13].

    中介軸承總摩擦熱Q為

    中介軸承載荷摩擦熱Ql為

    中介軸承黏度摩擦熱Qν為

    摩擦熱單位統(tǒng)一采用W.

    2.2 潤(rùn)滑劑黏溫關(guān)系

    潤(rùn)滑劑黏度通常是溫度的函數(shù),由此可見(jiàn),溫度對(duì)于黏度影響很大.因此,在計(jì)算黏度摩擦熱時(shí),需要考慮潤(rùn)滑劑的黏溫關(guān)系曲線.根據(jù)文獻(xiàn)[23],潤(rùn)滑劑在30,40,50,60,70,80?C 下的運(yùn)動(dòng)黏度分別為15,10,7.8,5.9,5,4 mm2/s.

    Reynolds 黏溫模型[28]如下

    式中,γ 表示黏溫系數(shù),礦物油有γ=0.018~0.036?C?1;T0表示初始溫度;ν0表示初始溫度下的初始運(yùn)動(dòng)黏度.

    利用Reynolds 黏溫模型擬合的實(shí)驗(yàn)數(shù)據(jù)[23],擬合曲線如圖3 所示.擬合曲線的黏溫系數(shù)為γ=0.026?C?1,這表明Reynolds 黏溫模型適用于該潤(rùn)滑劑.根據(jù)擬合出的黏溫關(guān)系曲線,可以得到任意溫度下潤(rùn)滑劑的運(yùn)動(dòng)黏度.此外,潤(rùn)滑劑運(yùn)動(dòng)黏度隨溫度變化范圍很大,這進(jìn)一步證明了考慮潤(rùn)滑劑黏溫關(guān)系的必要性.

    2.3 中介軸承熱傳遞模型

    由于軸承鋼GCr15 的Biot 數(shù)通常比較小(Bi<0.1),可以采用集總參數(shù)法[29]建立中介軸承熱傳遞模型,忽略滾子、內(nèi)圈和外圈的內(nèi)部溫差,從而大大簡(jiǎn)化了建模難度.中介軸承的熱結(jié)點(diǎn)及熱傳遞網(wǎng)絡(luò)如圖4 所示,圖4(a) 給出了中介軸承傳熱過(guò)程中的6 個(gè)熱結(jié)點(diǎn),分別為中介軸承滾子、內(nèi)圈、外圈、低壓軸段、高壓軸段和潤(rùn)滑劑.圖4(b)給出了中介軸承熱傳遞網(wǎng)絡(luò),Rri,Rro,Ri,Ro為熱傳導(dǎo)的接觸熱阻,RLr,RLi,RLo,RLP,RHP為對(duì)流傳熱的熱阻.

    圖3 潤(rùn)滑劑黏溫關(guān)系曲線Fig.3 Lubricant viscosity-temperature curve

    圖4 中介軸承的熱結(jié)點(diǎn)及熱傳遞網(wǎng)絡(luò)Fig.4 Thermal nodes and heat transfer network of the inter-shaft bearing

    潤(rùn)滑劑溫度假設(shè)為T(mén)L=(Tr+Ti+To)/3,即中介軸承滾子、內(nèi)圈和外圈溫度的平均值[14].此外,總摩擦熱假設(shè)為平均分配[14],即分配給滾子的摩擦熱Qr=0.5Qt,分配給內(nèi)圈的摩擦熱Qi=0.25Qt,分配給外圈的摩擦熱Qo=0.25Qt.對(duì)每個(gè)熱結(jié)點(diǎn)列寫(xiě)熱平衡方程,即可得到中介軸承熱傳遞控制方程

    式中,t為時(shí)間,csteel表示鋼的比熱容,mr,mi,mo,mLP,mHP分別是滾子、內(nèi)圈、外圈、低壓軸段和高壓軸段的質(zhì)量;內(nèi)圈?低壓軸段熱阻Ri=ln(di/dL)/(2πksteelB),外圈?高壓軸段熱阻Ro=ln(dH/do)/(2πksteelB);滾子?內(nèi)圈熱阻Rri=滾子?外圈熱阻Rro=式中nb為平均受力滾子數(shù)[30],Rone=1.13/(ksteel為單個(gè)滾子與內(nèi)圈或外圈的熱阻,其中A表示滾子與內(nèi)圈或外圈接觸面面積,Pe?表示修正的Peclet 數(shù)[31].

    而對(duì)流傳熱熱阻的計(jì)算通常簡(jiǎn)化為Nusselt 數(shù)的計(jì)算,通過(guò)下式

    式中,h為表面?zhèn)鳠嵯禂?shù),k為流體熱導(dǎo)率,L為特征長(zhǎng)度,A為對(duì)流傳熱面積.

    潤(rùn)滑劑?滾子的Nusselt 數(shù)[32]準(zhǔn)則方程為

    式中,Re=VL/ν 為Reynolds 數(shù),Pr=ν/α 為Prandtl數(shù),該準(zhǔn)則方程適用于10?1<Re<105.

    潤(rùn)滑劑?內(nèi)圈或外圈的Nusselt 數(shù)[33]準(zhǔn)則方程為

    式中Ta=為T(mén)aylor 數(shù),其中δio=(do?di)/2表示中介軸承內(nèi)外圈的間距,r表示滾道半徑.

    高壓軸段或低壓軸段?環(huán)境Nusselt 數(shù)[34]準(zhǔn)則方程為

    3 結(jié)果與討論

    3.1 中介軸承動(dòng)載荷非線性行為

    由于中介軸承彈性恢復(fù)力存在不可忽視的非線性因素,導(dǎo)致雙轉(zhuǎn)子系統(tǒng)動(dòng)力學(xué)方程也是非線性的.本文采用四階Runge-Kutta 法求解系統(tǒng)動(dòng)力學(xué)方程(1)~(8),得到了雙轉(zhuǎn)子系統(tǒng)的非線性動(dòng)力學(xué)響應(yīng),根據(jù)中介軸承彈性恢復(fù)力表達(dá)式(9),中介軸承彈性恢復(fù)力同時(shí)也就得到了.圖5 表示當(dāng)轉(zhuǎn)速ω1=715 rad/s時(shí)豎直方向水平方向上的彈性恢復(fù)力的時(shí)間歷程曲線,雙轉(zhuǎn)子系統(tǒng)參數(shù):轉(zhuǎn)速比λ=1.2,低壓轉(zhuǎn)子偏心距e1=3 μm 以及高壓轉(zhuǎn)子偏心距e2=2 μm.豎直方向和水平方向的中介軸承彈性恢復(fù)力都是時(shí)變的,其中水平方向彈性恢復(fù)力Fy在0 附近周期性變化,而豎直方向彈性恢復(fù)力Fx在600 N 附近周期性變化,這是因?yàn)橄到y(tǒng)在豎直方向存在重力.中介軸承彈性恢復(fù)力在豎直方向和水平方向均是時(shí)變的,因此,當(dāng)雙轉(zhuǎn)子系統(tǒng)以特定轉(zhuǎn)速運(yùn)行時(shí),采用彈性恢復(fù)力的形式來(lái)描述中介軸承所承受的載荷顯然是比較復(fù)雜的,這樣并不利于中介軸承受力情況的描述和分析.

    圖5 中介軸承彈性恢復(fù)力Fig.5 Restoring forces of the inter-shaft bearing

    圖6 表示動(dòng)載荷以及振幅隨轉(zhuǎn)速變化曲線.升速曲線表示轉(zhuǎn)速由低到高的過(guò)程,而降速曲線表示轉(zhuǎn)速由高到低的過(guò)程.與中介軸承彈性恢復(fù)力相比,當(dāng)雙轉(zhuǎn)子系統(tǒng)以特定轉(zhuǎn)速運(yùn)行時(shí),中介軸承動(dòng)載荷變成了一個(gè)恒定的常數(shù),因此,利用動(dòng)載荷來(lái)描述中介軸承的受力狀況顯然比時(shí)變的彈性恢復(fù)力更加簡(jiǎn)單方便.

    圖6 動(dòng)載荷以及振幅隨轉(zhuǎn)速變化曲線Fig.6 The dynamic load and the vibration amplitude against the rotation speed

    由圖6(b) 中雙轉(zhuǎn)子系統(tǒng)幅頻曲線可以看出,在ω1=570~900 rad/s 轉(zhuǎn)速區(qū)間內(nèi),系統(tǒng)存在兩個(gè)共振區(qū),分別由高壓轉(zhuǎn)子高低壓轉(zhuǎn)子的不平衡激勵(lì)引起.升速過(guò)程中,振幅急劇增大,直至轉(zhuǎn)速達(dá)到和振幅驟降,即發(fā)生了跳躍現(xiàn)象;降速過(guò)程中,當(dāng)轉(zhuǎn)速達(dá)到振幅驟增,即發(fā)生了跳躍現(xiàn)象.當(dāng)轉(zhuǎn)速處于ω1∈和ω1∈升速曲線和降速曲線沒(méi)有重合在一起,即出現(xiàn)雙穩(wěn)態(tài)現(xiàn)象.從數(shù)學(xué)上,這意味著系統(tǒng)動(dòng)力方程(1)~(8)在兩個(gè)雙穩(wěn)態(tài)區(qū)間和內(nèi)存在兩個(gè)穩(wěn)定解,系統(tǒng)動(dòng)力方程收斂到哪一個(gè)穩(wěn)定解取決于系統(tǒng)運(yùn)動(dòng)的初始狀態(tài).升速過(guò)程和降速過(guò)程中,系統(tǒng)在4 個(gè)跳躍點(diǎn)Adown,Aup,Bdown和Bup的運(yùn)動(dòng)初始狀態(tài)不同,因此升速曲線和降速曲線在雙穩(wěn)態(tài)區(qū)間不重合.由圖6(a)可以看出,中介軸承動(dòng)載荷變化規(guī)律與系統(tǒng)振幅變化規(guī)律一致.動(dòng)載荷在4 個(gè)跳躍點(diǎn)Adown,Aup,Bdown和Bup處也發(fā)生跳躍現(xiàn)象,而在兩個(gè)雙穩(wěn)態(tài)區(qū)間和內(nèi)也發(fā)生雙穩(wěn)態(tài)現(xiàn)象.總而言之,雙轉(zhuǎn)子系統(tǒng)動(dòng)力學(xué)建模中,考慮了中介軸承徑向游隙、Hertz 接觸的分?jǐn)?shù)指數(shù)關(guān)系以及變剛度等非線性因素,同時(shí)根據(jù)系統(tǒng)的非線性動(dòng)力學(xué)影響定義了中介軸承的動(dòng)載荷,因此,動(dòng)載荷會(huì)發(fā)生跳躍和雙穩(wěn)態(tài)等非線性現(xiàn)象[35],這也說(shuō)明中介軸承動(dòng)載荷能夠在一定程度上反映系統(tǒng)的非線性動(dòng)力學(xué)特性.

    3.2 中介軸承非線性熱行為

    中介軸承熱傳遞控制方程(18)~(22) 本質(zhì)上是一階常微分方程組,本文采用四階Runge-Kutta 法進(jìn)行數(shù)值求解,計(jì)算精度默認(rèn)為<10?6,初始溫度等于環(huán)境溫度T∞=20?C,計(jì)算時(shí)間設(shè)置為10 000 s,計(jì)算時(shí)間應(yīng)足夠大,確保溫度收斂于穩(wěn)態(tài)溫度.參數(shù)取值如下:轉(zhuǎn)速比λ=1.2、低壓轉(zhuǎn)子偏心距e1=3 μm、高壓轉(zhuǎn)子偏心距e2=2μm.

    當(dāng)中介軸承處于穩(wěn)態(tài)熱傳遞狀態(tài)時(shí),其滾子、內(nèi)圈以及外圈溫度隨低壓轉(zhuǎn)子轉(zhuǎn)速變化曲線如圖7(a)所示,圖7(b)為中介軸承線性化時(shí)滾子、內(nèi)圈以及外圈溫度隨低壓轉(zhuǎn)子轉(zhuǎn)速變化曲線,中介軸承線性化方法及過(guò)程參見(jiàn)文獻(xiàn)[25].由圖7(a)可以看到,無(wú)論是升速曲線還是降速曲線,都有滾子溫度Tr高于內(nèi)圈溫度Ti高于外圈溫度To,即Tr>Ti>To.同時(shí)發(fā)現(xiàn),升速曲線和降速曲線的Tr,Ti和To隨轉(zhuǎn)速的變化規(guī)律基本相同,因此,后文以Tr為例對(duì)中介軸承的非線性熱行為進(jìn)行詳細(xì)分析.在共振區(qū)A和B內(nèi),升速曲線Tr急速升高,直到轉(zhuǎn)速增至出現(xiàn)跳躍現(xiàn)象,即Tr驟降;降速曲線Tr緩慢降低,直到轉(zhuǎn)速降至出現(xiàn)跳躍現(xiàn)象,即Tr驟增.當(dāng)轉(zhuǎn)速處于范圍內(nèi)時(shí),升速曲線和降速曲線的Tr不重合,即出現(xiàn)了雙穩(wěn)態(tài)現(xiàn)象.對(duì)比圖7(a)和圖7(b),可以看出,當(dāng)中介軸承線性化之后,其熱行為的跳躍和雙穩(wěn)態(tài)現(xiàn)象均會(huì)消失不見(jiàn)!

    圖7 滾子、內(nèi)圈和外圈的穩(wěn)態(tài)溫度隨轉(zhuǎn)速變化曲線(實(shí)線為升速曲線,虛線為降速曲線)Fig.7 Steady-state temperatures of rollers,inner race and outer race against rotation speed(solid lines represent run-up curves,dotted lines represent run-down curves)

    為了便于后文進(jìn)一步分析中介軸承的非線性熱行為,現(xiàn)定義相關(guān)符號(hào)和參數(shù):Adown,Aup,Bdown和Bup表示“跳躍點(diǎn)”;表示“跳躍點(diǎn)頻率”;表示“跳躍幅值”;表示“雙穩(wěn)態(tài)區(qū)間”.

    為了進(jìn)一步揭示非線性熱行為的內(nèi)在機(jī)理,非常有必要進(jìn)行摩擦熱分析,包括總摩擦熱分析、載荷摩擦熱分析和黏度摩擦熱分析.圖8 為中介軸承總摩擦熱、載荷摩擦熱和黏度摩擦熱隨轉(zhuǎn)速的變化曲線,其中實(shí)線表示升速曲線,虛線表示降速曲線.可以看到,升速過(guò)程中,總摩擦熱Q、載荷摩擦熱Ql和黏度摩擦熱Qν在跳躍點(diǎn)Aup和Bup處均出現(xiàn)跳躍現(xiàn)象;升速過(guò)程中,三者在跳躍點(diǎn)Adown和Bdown處均出現(xiàn)跳躍現(xiàn)象.此外,在Aup和Bup處,Ql向下跳躍而Qν向上跳躍,但Ql跳躍幅值顯著大于Qν跳躍幅值,因此Q與Ql跳躍方向相同;在Adown和Bdown處,Ql向上跳躍而Qν向下跳躍,但Ql跳躍幅值顯著大于Qν跳躍幅值,因此Q與Ql跳躍方向相同.同時(shí),三者在雙穩(wěn)態(tài)區(qū)間?ωA和?ωB內(nèi)均出現(xiàn)雙穩(wěn)態(tài)現(xiàn)象.

    圖8 中介軸承摩擦熱隨轉(zhuǎn)速變化曲線Fig.8 Friction heat of the inter-shaft bearing against rotation speed

    總而言之,中介軸承溫度出現(xiàn)非線性熱行為,即跳躍和雙穩(wěn)態(tài)現(xiàn)象,直接原因是總摩擦熱出現(xiàn)非線性行為,而根本原因是雙轉(zhuǎn)子系統(tǒng)非線性動(dòng)力學(xué)特性所導(dǎo)致的中介軸承動(dòng)載荷出現(xiàn)的非線性行為.之所以產(chǎn)生這一獨(dú)特的現(xiàn)象是因?yàn)楸疚囊胫薪檩S承動(dòng)載荷代入Palmgren 經(jīng)驗(yàn)公式來(lái)計(jì)算載荷摩擦熱,如果仍舊采用靜載荷則無(wú)法發(fā)現(xiàn)中介軸承的非線性熱行為.

    3.3 轉(zhuǎn)速比的影響

    圖9 不同轉(zhuǎn)速比下滾子溫度隨轉(zhuǎn)速變化曲線Fig.9 Variation for temperature of rollers with rotation speed under different rotation speed ratio

    圖9 為不同轉(zhuǎn)速比下滾子溫度隨轉(zhuǎn)速變化曲線對(duì)比,轉(zhuǎn)速比分別取λ=1.1,λ=1.15,λ=1.2和λ=1.25.可見(jiàn),轉(zhuǎn)速比λ 不僅影響滾子溫度Tr而且影響共振區(qū)A和B內(nèi)中介軸承非線性熱行為.隨著λ 的增大,Tr顯著增大;跳躍點(diǎn)頻率基本不變,而則顯著降低;跳躍幅值均明顯增大;雙穩(wěn)態(tài)區(qū)間?ωA和?ωB長(zhǎng)度也基本不變,但?ωB略寬于?ωA.

    3.4 偏心距的影響

    首先分析低壓轉(zhuǎn)子偏心距的影響,圖10 為不同低壓轉(zhuǎn)子偏心距下滾子溫度隨轉(zhuǎn)速變化曲線對(duì)比,低壓轉(zhuǎn)子偏心距分別取e1=2μm,e1=3μm,e1=4μm和e1=5 μm.可見(jiàn),滾子溫度Tr除共振區(qū)B外幾乎重合,表明低壓轉(zhuǎn)子偏心距e1只影響共振區(qū)B.隨著e1的增大,跳躍點(diǎn)頻率基本不變,而則顯著增大;跳躍幅值基本不變,而則明顯增大;雙穩(wěn)態(tài)區(qū)間?ωA長(zhǎng)度基本不變,而?ωB則變窄.

    圖10 不同低壓轉(zhuǎn)子偏心距下滾子溫度隨轉(zhuǎn)速變化曲線Fig.10 Temperature and friction heat against rotation speed under different unbalances of lower pressure rotor

    圖11 不同高壓轉(zhuǎn)子偏心距下滾子溫度隨轉(zhuǎn)速變化曲線Fig.11 Temperature and friction heat against rotation speed under different unbalances of higher pressure rotor

    最后分析高壓轉(zhuǎn)子偏心距的影響,圖11 為不同高壓轉(zhuǎn)子偏心距下滾子溫度隨轉(zhuǎn)速變化曲線對(duì)比,高壓轉(zhuǎn)子偏心距分別取e2=2μm,e2=3μm,e2=4μm和e2=5 μm.可見(jiàn),滾子溫度Tr除共振區(qū)A外幾乎重合,表明高壓轉(zhuǎn)子偏心距e2只影響共振區(qū)A.隨著e2的增大,跳躍點(diǎn)頻率基本不變,而則顯著增大;跳躍幅值基本不變,而則明顯增大;雙穩(wěn)態(tài)區(qū)間?ωB長(zhǎng)度基本不變,而?ωA則變窄.

    總而言之,低壓轉(zhuǎn)子偏心距只影響共振區(qū)B,而高壓轉(zhuǎn)子偏心距只影響共振區(qū)A,隨著對(duì)應(yīng)的偏心距的增大,對(duì)應(yīng)的跳躍點(diǎn)頻率和跳躍幅度均增大,而對(duì)應(yīng)的雙穩(wěn)態(tài)區(qū)間則變窄.

    3.5 中介軸承徑向游隙的影響

    圖12 不同徑向游隙下滾子溫度隨轉(zhuǎn)速變化曲線Fig.12 Temperature and friction heat against rotation speed under different radial clearance

    圖12 為不同徑向游隙下滾子溫度隨轉(zhuǎn)速變化曲線對(duì)比,徑向游隙分別取δ0=3 μm,δ0=5 μm,δ0=8 μm 和δ0=10 μm.可見(jiàn),徑向游隙δ 主要影響共振區(qū)A和B內(nèi)中介軸承非線性熱行為,隨著δ 的增大,滾子溫度Tr的最大值幾乎不變;跳躍點(diǎn)頻率均略微降低;跳躍幅值均略微增大;雙穩(wěn)態(tài)區(qū)間?ωA和?ωB均明顯變寬.值得注意的是,當(dāng)δ0=3μm 時(shí),中介軸承跳躍現(xiàn)象和雙穩(wěn)態(tài)現(xiàn)象等非線性熱行為“消失了”,這意味著徑向游隙的存在是引起中介軸承溫度出現(xiàn)非線性熱行為的因素之一.

    中介軸承徑向游隙延遲了滾子與內(nèi)圈和外圈的接觸時(shí)間,這相當(dāng)于以另一種方式降低了中介軸承的剛度,從而降低了跳躍點(diǎn)頻率.此外,徑向游隙本質(zhì)上是一種非線性因素,適當(dāng)減小徑向間隙將顯著抑制中介軸承非線性熱行為,即減小跳躍幅值和縮窄雙穩(wěn)態(tài)區(qū)間,當(dāng)徑向間隙減小到一定值時(shí),中介軸承非線性熱行為甚至?xí)?

    3.6 中介軸承Hertz 接觸剛度的影響

    圖13 為不同中介軸承Hertz 接觸剛度下滾子溫度隨轉(zhuǎn)速變化曲線對(duì)比,徑向游隙分別取Kb=8Kb0,Kb=10Kb0,Kb=15Kb0和Kb=20Kb0(Kb0=107N/m10/9).可見(jiàn),中介軸承Hertz 接觸剛度Kb主要影響共振區(qū)A和B內(nèi)中介軸承非線性熱行為,隨著Kb的增大,滾子溫度Tr的最大值幾乎不變;跳躍點(diǎn)頻率均顯著增大;跳躍幅值均略微減小;雙穩(wěn)態(tài)區(qū)間?ωA和?ωB也略微變窄.

    圖13 不同中介軸承Hertz 接觸剛度下滾子溫度隨轉(zhuǎn)速變化曲線Fig.13 Temperature and friction heat against rotation speed under different Hertz contact stiffness of the inter-shaft bearing

    3.7 中介軸承滾子數(shù)目的影響

    圖14 為不同滾子數(shù)目下滾子溫度隨轉(zhuǎn)速變化曲線對(duì)比,徑向游隙分別取Nb=12,Nb=16,Nb=20和Nb=24.可見(jiàn),滾子數(shù)目Nb不僅影響滾子溫度Tr而且影響共振區(qū)A和B內(nèi)中介軸承非線性熱行為.隨著Nb的增大,Tr顯著減小;跳躍點(diǎn)頻率均顯著增大;跳躍幅值均略微減小;雙穩(wěn)態(tài)區(qū)間?ωA和?ωB也略微變窄.

    圖14 不同滾子數(shù)目下滾子溫度隨轉(zhuǎn)速變化曲線Fig.14 Temperature and friction heat against rotation speed under different roller numbers of the inter-shaft bearing

    4 結(jié)論

    本文考慮中介軸承的徑向游隙、分?jǐn)?shù)指數(shù)非線性和參數(shù)激勵(lì)等多種非線性因素,基于雙轉(zhuǎn)子系統(tǒng)非線性動(dòng)力響應(yīng)定義了中介軸承動(dòng)載荷,同時(shí)考慮了潤(rùn)滑劑的黏溫關(guān)系,根據(jù)Palmgren 經(jīng)驗(yàn)公式建立了動(dòng)載荷作用下的中介軸承熱傳遞模型,研究了動(dòng)載荷非線性行為對(duì)中介軸承溫度的影響,并進(jìn)行了詳細(xì)的參數(shù)分析.主要結(jié)論如下:

    (1)中介軸承的非線性動(dòng)載荷將導(dǎo)致中介軸承溫度出現(xiàn)非線性熱行為,即表現(xiàn)出跳躍和雙穩(wěn)態(tài)現(xiàn)象.

    (2)轉(zhuǎn)速比和滾子數(shù)目對(duì)中介軸承溫度和非線性熱行為都有重要影響,降低轉(zhuǎn)速比或增加滾子數(shù)目不僅能夠有效降低溫度,而且有利于抑制非線性熱行為.

    (3)偏心距、中介軸承徑向游隙和Hertz 接觸剛度均只影響中介軸承非線性熱行為,降低轉(zhuǎn)子偏心距、適當(dāng)縮小徑向游隙或增大Hertz 接觸剛度均有利于抑制非線性熱行為.

    本文研究表明,雙轉(zhuǎn)子系統(tǒng)的非線性動(dòng)力學(xué)特性會(huì)使中介軸承熱行為表現(xiàn)出復(fù)雜非線性行為,未來(lái)的工作將集中在非線性熱行為的實(shí)驗(yàn)驗(yàn)證上.

    猜你喜歡
    雙穩(wěn)態(tài)游隙滾子
    一維有界區(qū)域上雙穩(wěn)態(tài)方程多重正解的存在性
    六層非對(duì)稱正交雙穩(wěn)態(tài)復(fù)合材料層合板的動(dòng)態(tài)跳躍研究1)
    含彈性碰撞作用的雙級(jí)雙穩(wěn)態(tài)結(jié)構(gòu)振動(dòng)能量采集研究
    圓錐滾子軸承半凸?jié)L子的優(yōu)化研究
    哈爾濱軸承(2021年1期)2021-07-21 05:43:10
    仿真模擬在多聯(lián)推力滾子軸承研發(fā)中的應(yīng)用
    哈爾濱軸承(2021年4期)2021-03-08 01:00:52
    軸承游隙對(duì)隔爆電機(jī)的影響分析
    水平軸設(shè)計(jì)中軸承游隙的計(jì)算與選擇
    滿滾子徑向自鎖軸承的優(yōu)化設(shè)計(jì)
    軸承(2018年10期)2018-07-25 07:22:22
    基于雙穩(wěn)態(tài)的振動(dòng)能量收集系統(tǒng)的設(shè)計(jì)
    雙列球面滾子軸承力學(xué)分析及滾子受載計(jì)算方法
    h日本视频在线播放| 亚洲内射少妇av| 免费大片18禁| 非洲黑人性xxxx精品又粗又长| 欧美日韩亚洲国产一区二区在线观看| 欧美激情国产日韩精品一区| 亚洲精品在线美女| 丰满的人妻完整版| 69av精品久久久久久| 无人区码免费观看不卡| 看免费av毛片| 亚洲精品一卡2卡三卡4卡5卡| 久久久久亚洲av毛片大全| 天堂动漫精品| 久久久国产成人免费| 99热这里只有精品一区| 日韩国内少妇激情av| ponron亚洲| 日韩欧美一区二区三区在线观看| 免费搜索国产男女视频| 欧美成人免费av一区二区三区| 小蜜桃在线观看免费完整版高清| 在线观看av片永久免费下载| 欧美黑人欧美精品刺激| 久久性视频一级片| 国内精品久久久久精免费| 久久精品国产亚洲av涩爱 | 97超级碰碰碰精品色视频在线观看| 精品久久久久久久久久免费视频| 国产老妇女一区| www日本黄色视频网| 国产精品自产拍在线观看55亚洲| 国产亚洲欧美98| 赤兔流量卡办理| 精品一区二区三区视频在线观看免费| 色综合亚洲欧美另类图片| 观看美女的网站| 亚洲人与动物交配视频| 国产在线男女| 中文字幕av成人在线电影| 天堂动漫精品| 一进一出抽搐动态| 色尼玛亚洲综合影院| 免费看光身美女| 亚洲精品在线观看二区| 亚洲狠狠婷婷综合久久图片| 日本撒尿小便嘘嘘汇集6| 99久久精品一区二区三区| 亚洲电影在线观看av| 草草在线视频免费看| 日日干狠狠操夜夜爽| 男插女下体视频免费在线播放| 色av中文字幕| .国产精品久久| 亚洲精品成人久久久久久| 中文字幕熟女人妻在线| 国产人妻一区二区三区在| 九色国产91popny在线| 内射极品少妇av片p| 嫩草影视91久久| 天堂√8在线中文| 国模一区二区三区四区视频| 又爽又黄a免费视频| 激情在线观看视频在线高清| 亚洲美女视频黄频| 日韩高清综合在线| 色尼玛亚洲综合影院| 成人国产综合亚洲| 日韩欧美一区二区三区在线观看| 色精品久久人妻99蜜桃| 尤物成人国产欧美一区二区三区| 成人特级黄色片久久久久久久| www.熟女人妻精品国产| 日本黄色视频三级网站网址| 国产精品影院久久| 亚洲国产高清在线一区二区三| 亚洲avbb在线观看| 精品久久久久久久久久久久久| 中文资源天堂在线| 老司机深夜福利视频在线观看| 精品不卡国产一区二区三区| 人人妻,人人澡人人爽秒播| 成人国产一区最新在线观看| www.色视频.com| 亚洲激情在线av| 在线免费观看的www视频| 亚洲自偷自拍三级| 日韩高清综合在线| 老师上课跳d突然被开到最大视频 久久午夜综合久久蜜桃 | 日韩欧美国产一区二区入口| 亚洲性夜色夜夜综合| 国产成人aa在线观看| 人妻丰满熟妇av一区二区三区| 一区二区三区激情视频| 精品乱码久久久久久99久播| 久久久国产成人免费| 中文字幕久久专区| 不卡一级毛片| 最好的美女福利视频网| 中文字幕av在线有码专区| 琪琪午夜伦伦电影理论片6080| 长腿黑丝高跟| 午夜精品久久久久久毛片777| 乱码一卡2卡4卡精品| 久久国产精品影院| 久久中文看片网| 高清毛片免费观看视频网站| 国产大屁股一区二区在线视频| 亚洲经典国产精华液单 | 日本撒尿小便嘘嘘汇集6| 国产精品人妻久久久久久| 神马国产精品三级电影在线观看| 最近中文字幕高清免费大全6 | 久久精品国产99精品国产亚洲性色| 中文资源天堂在线| 国产主播在线观看一区二区| 午夜老司机福利剧场| 中文字幕人妻熟人妻熟丝袜美| 成人国产综合亚洲| 草草在线视频免费看| 18美女黄网站色大片免费观看| 成人亚洲精品av一区二区| 757午夜福利合集在线观看| 美女免费视频网站| 色精品久久人妻99蜜桃| 成人美女网站在线观看视频| 波多野结衣高清无吗| 噜噜噜噜噜久久久久久91| 非洲黑人性xxxx精品又粗又长| 国产精品不卡视频一区二区 | 久久草成人影院| 天天一区二区日本电影三级| 婷婷色综合大香蕉| 久久午夜福利片| 国产大屁股一区二区在线视频| 精品国产亚洲在线| 丰满乱子伦码专区| 麻豆av噜噜一区二区三区| 成人永久免费在线观看视频| 男女做爰动态图高潮gif福利片| 在线观看午夜福利视频| av在线老鸭窝| 成人特级黄色片久久久久久久| 精品一区二区三区人妻视频| 噜噜噜噜噜久久久久久91| 亚洲美女视频黄频| 国产高清激情床上av| 丰满的人妻完整版| 亚洲自拍偷在线| 色吧在线观看| 欧美激情在线99| 日韩欧美免费精品| 午夜福利视频1000在线观看| 欧美3d第一页| 国产在线精品亚洲第一网站| 成人永久免费在线观看视频| 不卡一级毛片| 久久国产乱子伦精品免费另类| 久久99热这里只有精品18| 亚洲精品日韩av片在线观看| 97人妻精品一区二区三区麻豆| 国产精品一及| 国产精品免费一区二区三区在线| 变态另类丝袜制服| 国产精品久久久久久人妻精品电影| 欧美3d第一页| 午夜两性在线视频| 淫秽高清视频在线观看| 欧美另类亚洲清纯唯美| av在线蜜桃| 午夜福利免费观看在线| 亚洲不卡免费看| 亚洲人成网站在线播放欧美日韩| 日韩欧美三级三区| 一级作爱视频免费观看| 99国产精品一区二区三区| 男女下面进入的视频免费午夜| 久久久久国内视频| 看黄色毛片网站| 国产午夜精品久久久久久一区二区三区 | 国产麻豆成人av免费视频| 国产伦一二天堂av在线观看| 少妇被粗大猛烈的视频| 久久久国产成人免费| 在线播放无遮挡| 丁香欧美五月| 亚洲精品影视一区二区三区av| 久久人妻av系列| 每晚都被弄得嗷嗷叫到高潮| av视频在线观看入口| 观看免费一级毛片| 日本 欧美在线| 国产色婷婷99| 最好的美女福利视频网| 禁无遮挡网站| 搞女人的毛片| 搡女人真爽免费视频火全软件 | 国产91精品成人一区二区三区| 真人一进一出gif抽搐免费| 色综合婷婷激情| 少妇被粗大猛烈的视频| 麻豆成人午夜福利视频| 国产精品乱码一区二三区的特点| 精品一区二区免费观看| 欧美xxxx性猛交bbbb| 国产免费男女视频| 男女床上黄色一级片免费看| 真人一进一出gif抽搐免费| 日本在线视频免费播放| 免费人成在线观看视频色| 美女被艹到高潮喷水动态| 搞女人的毛片| 99久久九九国产精品国产免费| 国产亚洲精品久久久久久毛片| av天堂在线播放| 国产欧美日韩一区二区三| 69人妻影院| 欧美色欧美亚洲另类二区| 国产精品电影一区二区三区| 亚洲第一区二区三区不卡| 男人舔女人下体高潮全视频| 我要看日韩黄色一级片| 午夜激情欧美在线| 看免费av毛片| 午夜福利在线观看免费完整高清在 | 99久久九九国产精品国产免费| 亚洲国产精品999在线| 99久久精品热视频| 性色av乱码一区二区三区2| 在线免费观看不下载黄p国产 | 精品人妻偷拍中文字幕| 自拍偷自拍亚洲精品老妇| 女生性感内裤真人,穿戴方法视频| 国产欧美日韩一区二区精品| 国产又黄又爽又无遮挡在线| 欧美区成人在线视频| 草草在线视频免费看| 欧美激情在线99| 亚洲av电影不卡..在线观看| 麻豆一二三区av精品| 精品国产亚洲在线| 丰满人妻一区二区三区视频av| 国产精品久久视频播放| 亚洲av电影不卡..在线观看| 老司机深夜福利视频在线观看| 国产野战对白在线观看| 男女视频在线观看网站免费| 99久久精品一区二区三区| eeuss影院久久| 成人av在线播放网站| 老熟妇乱子伦视频在线观看| 最近视频中文字幕2019在线8| av在线老鸭窝| 少妇被粗大猛烈的视频| 色视频www国产| 亚洲国产高清在线一区二区三| 亚洲av一区综合| 男女那种视频在线观看| 一区福利在线观看| 首页视频小说图片口味搜索| 身体一侧抽搐| 十八禁人妻一区二区| 国产精品亚洲美女久久久| 美女高潮喷水抽搐中文字幕| 在线播放国产精品三级| 午夜福利18| 给我免费播放毛片高清在线观看| 亚洲色图av天堂| 90打野战视频偷拍视频| 精品人妻一区二区三区麻豆 | 国产三级黄色录像| 一区二区三区四区激情视频 | 欧美在线一区亚洲| 亚洲国产精品久久男人天堂| 免费av不卡在线播放| 少妇人妻一区二区三区视频| 久久久久久久久久黄片| 欧美乱妇无乱码| 国产在线精品亚洲第一网站| 麻豆久久精品国产亚洲av| 欧美黑人欧美精品刺激| 久久99热这里只有精品18| 性色avwww在线观看| 国产精品一区二区免费欧美| 亚洲色图av天堂| 亚洲 国产 在线| 深爱激情五月婷婷| 国产精品久久电影中文字幕| 亚洲经典国产精华液单 | 亚洲男人的天堂狠狠| 欧美激情国产日韩精品一区| 日韩成人在线观看一区二区三区| 亚洲av第一区精品v没综合| 亚洲av电影不卡..在线观看| 国产精品98久久久久久宅男小说| 欧美午夜高清在线| 青草久久国产| 国产欧美日韩精品一区二区| 国产精品综合久久久久久久免费| 琪琪午夜伦伦电影理论片6080| 欧美黄色淫秽网站| 亚洲国产精品合色在线| 一个人免费在线观看的高清视频| 亚洲三级黄色毛片| 亚洲18禁久久av| 校园春色视频在线观看| av在线观看视频网站免费| 夜夜爽天天搞| 在线十欧美十亚洲十日本专区| АⅤ资源中文在线天堂| 国内精品久久久久精免费| 精品一区二区免费观看| 精品久久久久久久久久久久久| 少妇丰满av| 麻豆av噜噜一区二区三区| 又粗又爽又猛毛片免费看| 欧美色视频一区免费| 国产视频一区二区在线看| 亚洲狠狠婷婷综合久久图片| 久久精品国产99精品国产亚洲性色| 国产综合懂色| 全区人妻精品视频| 精华霜和精华液先用哪个| 成人国产综合亚洲| 亚洲精品一区av在线观看| 美女xxoo啪啪120秒动态图 | 日韩免费av在线播放| 3wmmmm亚洲av在线观看| 亚洲欧美日韩卡通动漫| 蜜桃久久精品国产亚洲av| 久久精品国产亚洲av涩爱 | 亚洲,欧美精品.| 久9热在线精品视频| 99国产精品一区二区三区| 五月玫瑰六月丁香| bbb黄色大片| 又黄又爽又免费观看的视频| 深夜精品福利| 最新在线观看一区二区三区| av福利片在线观看| or卡值多少钱| 国产在视频线在精品| 色在线成人网| 国产成年人精品一区二区| 国模一区二区三区四区视频| 免费观看精品视频网站| 国产乱人伦免费视频| a级毛片a级免费在线| www.熟女人妻精品国产| 亚洲中文字幕日韩| 成人特级av手机在线观看| av黄色大香蕉| 国产高潮美女av| 制服丝袜大香蕉在线| 久久久久久久亚洲中文字幕 | 熟女人妻精品中文字幕| 国产精品人妻久久久久久| 久久国产乱子免费精品| 欧美性猛交╳xxx乱大交人| 91麻豆av在线| 久久久久久九九精品二区国产| 亚洲七黄色美女视频| 久久人妻av系列| 深夜a级毛片| 男人和女人高潮做爰伦理| 欧美日韩瑟瑟在线播放| 哪里可以看免费的av片| 很黄的视频免费| 婷婷精品国产亚洲av| 欧美在线黄色| 亚洲五月婷婷丁香| 亚洲最大成人手机在线| 久久精品国产亚洲av涩爱 | 五月玫瑰六月丁香| 亚洲男人的天堂狠狠| 国产精品av视频在线免费观看| 别揉我奶头 嗯啊视频| 国产精品久久久久久亚洲av鲁大| aaaaa片日本免费| 国产精品久久久久久久电影| 日韩亚洲欧美综合| 成人av一区二区三区在线看| 99久久无色码亚洲精品果冻| 国产亚洲精品久久久久久毛片| 岛国在线免费视频观看| 中文字幕精品亚洲无线码一区| 色综合婷婷激情| 日本黄色视频三级网站网址| 五月伊人婷婷丁香| 亚洲欧美日韩高清专用| 久久人妻av系列| 亚洲在线观看片| 久久久久久九九精品二区国产| 亚洲av免费在线观看| 性色avwww在线观看| 午夜久久久久精精品| 国产精品1区2区在线观看.| 综合色av麻豆| 成年人黄色毛片网站| 精品一区二区三区av网在线观看| 日本黄色视频三级网站网址| 国产精品一及| 久久99热6这里只有精品| 18+在线观看网站| 精品人妻偷拍中文字幕| 少妇的逼好多水| 最新在线观看一区二区三区| 久久精品国产亚洲av涩爱 | 亚洲人成伊人成综合网2020| www.熟女人妻精品国产| 99久久九九国产精品国产免费| 亚洲综合色惰| 天美传媒精品一区二区| 在现免费观看毛片| 国产综合懂色| 午夜福利高清视频| 亚洲自偷自拍三级| 亚洲人成电影免费在线| 国产aⅴ精品一区二区三区波| 婷婷亚洲欧美| 黄片小视频在线播放| 国产伦人伦偷精品视频| 国产人妻一区二区三区在| 日日摸夜夜添夜夜添av毛片 | 国产精品久久久久久久电影| 床上黄色一级片| 亚洲人成伊人成综合网2020| 久久国产精品人妻蜜桃| 色哟哟哟哟哟哟| 精品一区二区免费观看| 成年人黄色毛片网站| 人妻制服诱惑在线中文字幕| 一区福利在线观看| a级毛片免费高清观看在线播放| 内地一区二区视频在线| 国产精品1区2区在线观看.| 国产精品免费一区二区三区在线| 欧美3d第一页| 亚洲va日本ⅴa欧美va伊人久久| 日本黄色视频三级网站网址| 国产成人福利小说| 国产探花在线观看一区二区| 黄色视频,在线免费观看| 国产久久久一区二区三区| 国模一区二区三区四区视频| 欧美成狂野欧美在线观看| 尤物成人国产欧美一区二区三区| 国产精品影院久久| 国产私拍福利视频在线观看| 又爽又黄无遮挡网站| 亚洲一区二区三区不卡视频| 久久精品国产清高在天天线| 国产成年人精品一区二区| 啪啪无遮挡十八禁网站| av黄色大香蕉| 婷婷精品国产亚洲av在线| 久久久久久大精品| 亚洲av熟女| 久久久久精品国产欧美久久久| 亚洲在线自拍视频| 日本黄色视频三级网站网址| 俄罗斯特黄特色一大片| 亚洲av一区综合| 首页视频小说图片口味搜索| 精品福利观看| 亚洲专区中文字幕在线| 亚洲av免费在线观看| 国产成人啪精品午夜网站| 亚洲中文字幕一区二区三区有码在线看| av欧美777| 国产69精品久久久久777片| 国产精品不卡视频一区二区 | 亚洲熟妇熟女久久| 欧美性感艳星| 国产成人影院久久av| 成年人黄色毛片网站| 色尼玛亚洲综合影院| 亚洲国产高清在线一区二区三| 亚洲精品亚洲一区二区| 亚洲,欧美精品.| 久久久久久久精品吃奶| 岛国在线免费视频观看| 国产伦人伦偷精品视频| 国内精品久久久久久久电影| 观看免费一级毛片| 日韩国内少妇激情av| 日韩精品青青久久久久久| 久久人人爽人人爽人人片va | 两性午夜刺激爽爽歪歪视频在线观看| 欧美激情久久久久久爽电影| av天堂在线播放| 亚洲av日韩精品久久久久久密| 给我免费播放毛片高清在线观看| www.熟女人妻精品国产| 老熟妇乱子伦视频在线观看| 日韩欧美精品v在线| 成人特级av手机在线观看| 18禁黄网站禁片午夜丰满| 久久伊人香网站| 丁香六月欧美| 亚洲精品在线观看二区| 欧美一区二区国产精品久久精品| 欧美成人免费av一区二区三区| 国产极品精品免费视频能看的| 亚洲精品成人久久久久久| 亚洲最大成人av| 国产人妻一区二区三区在| 美女被艹到高潮喷水动态| 亚洲精品在线观看二区| 亚洲av美国av| 观看美女的网站| 亚洲成a人片在线一区二区| 亚洲在线观看片| 亚洲一区二区三区色噜噜| 欧美黑人欧美精品刺激| 久久午夜亚洲精品久久| 国产成人av教育| 99精品在免费线老司机午夜| 国产精品国产高清国产av| 成人特级黄色片久久久久久久| 国产精品av视频在线免费观看| 国产国拍精品亚洲av在线观看| 欧美一区二区亚洲| 女人被狂操c到高潮| 国产免费男女视频| 欧美bdsm另类| 免费搜索国产男女视频| 亚洲人成网站在线播| 白带黄色成豆腐渣| 草草在线视频免费看| 女生性感内裤真人,穿戴方法视频| 青草久久国产| 国产精品久久电影中文字幕| 成人高潮视频无遮挡免费网站| 日韩成人在线观看一区二区三区| 在线观看美女被高潮喷水网站 | or卡值多少钱| 99热这里只有是精品在线观看 | 国产老妇女一区| 免费av毛片视频| 丝袜美腿在线中文| 色尼玛亚洲综合影院| 99视频精品全部免费 在线| 国产精品av视频在线免费观看| 午夜福利在线观看吧| 国产又黄又爽又无遮挡在线| 国产亚洲精品久久久久久毛片| 久久久久久久久久黄片| 精品一区二区三区视频在线| 少妇的逼水好多| 亚洲国产精品sss在线观看| 亚洲精品粉嫩美女一区| 亚洲欧美日韩高清专用| 在线国产一区二区在线| 国产午夜精品久久久久久一区二区三区 | 婷婷精品国产亚洲av在线| 午夜精品一区二区三区免费看| 免费av不卡在线播放| 日韩欧美 国产精品| 国产精品一及| 18美女黄网站色大片免费观看| 怎么达到女性高潮| 成人av一区二区三区在线看| 国产一区二区在线av高清观看| 很黄的视频免费| 天天一区二区日本电影三级| 亚洲欧美日韩无卡精品| 老鸭窝网址在线观看| 99精品在免费线老司机午夜| 国产亚洲精品久久久久久毛片| 国产一区二区在线av高清观看| 热99re8久久精品国产| netflix在线观看网站| av欧美777| 青草久久国产| 人人妻人人看人人澡| 精品人妻视频免费看| 午夜老司机福利剧场| 亚洲av中文字字幕乱码综合| 国产午夜福利久久久久久| АⅤ资源中文在线天堂| 99热这里只有是精品在线观看 | 亚洲最大成人中文| 最近视频中文字幕2019在线8| 老司机福利观看| 免费av不卡在线播放| 国产av不卡久久| bbb黄色大片| 99久久九九国产精品国产免费| 国产男靠女视频免费网站| 简卡轻食公司| 一本综合久久免费| 男女做爰动态图高潮gif福利片| 俺也久久电影网| 看黄色毛片网站| 国产精品av视频在线免费观看| 国产一区二区三区在线臀色熟女| 午夜激情欧美在线| 1000部很黄的大片| 淫秽高清视频在线观看| 欧美日韩亚洲国产一区二区在线观看| 婷婷精品国产亚洲av| 亚洲熟妇熟女久久| 精品日产1卡2卡| 久久精品91蜜桃| 亚洲精品色激情综合| av黄色大香蕉| 久久欧美精品欧美久久欧美| 久久6这里有精品| 99热精品在线国产| 国产爱豆传媒在线观看| 欧美日韩福利视频一区二区| 特级一级黄色大片| 色综合亚洲欧美另类图片| 亚洲av免费高清在线观看|