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

    機械變速器齒輪敲擊現(xiàn)象建模及其影響因素

    2016-12-26 03:20:38吳虎威吳光強同濟大學汽車學院上海20804東京大學生產(chǎn)技術(shù)研究所東京538505
    同濟大學學報(自然科學版) 2016年1期
    關(guān)鍵詞:振動模型

    吳虎威,吳光強,2(.同濟大學汽車學院,上海20804;2.東京大學生產(chǎn)技術(shù)研究所,東京53 8505)

    ?

    機械變速器齒輪敲擊現(xiàn)象建模及其影響因素

    吳虎威1,吳光強1,2
    (1.同濟大學汽車學院,上海201804;2.東京大學生產(chǎn)技術(shù)研究所,東京153 8505)

    利用車輛基本參數(shù)和實車傳動系統(tǒng)各部件的等效轉(zhuǎn)動慣量、扭轉(zhuǎn)剛度和阻尼扭轉(zhuǎn)特性參數(shù),建立了包含LuGre輪胎模型在內(nèi)的傳動系統(tǒng)非線性扭轉(zhuǎn)振動和車輛縱向平移運動耦合模型.利用耦合模型再現(xiàn)了車輛處于一擋加速工況下各非承載齒輪對的敲擊情況,并分析了各非承載齒輪對敲擊強度的主要影響因素.結(jié)果表明,合理設(shè)計各擋位齒輪對齒側(cè)間隙和等效轉(zhuǎn)動慣量大小可有效抑制機械變速器齒輪敲擊現(xiàn)象.

    LuGre輪胎模型;非線性扭轉(zhuǎn)振動;耦合模型;齒側(cè)間隙;等效轉(zhuǎn)動慣量

    隨著汽車工業(yè)的蓬勃發(fā)展,用戶對于汽車乘坐的舒適性要求越來越高,汽車振動噪聲品質(zhì)越來越受到重視.變速器作為汽車動力傳動系統(tǒng)的重要組

    國內(nèi)外已有關(guān)于變速器齒輪敲擊問題的研究成果.在仿真建模方面,Wang等[2-3]較早利用基于集中質(zhì)量的系統(tǒng)動力學方法針對變速器系統(tǒng)建模.模型中包括承載齒輪對和非承載齒輪對,將單個齒輪等效為集中轉(zhuǎn)動慣量,齒輪間軸段等效為扭轉(zhuǎn)剛度和扭轉(zhuǎn)阻尼,仿真分析了變速器系統(tǒng)不同非承載齒輪對敲擊強度.對于傳動系統(tǒng)模型進行較大簡化,將半軸、輪胎和整車質(zhì)量等效為較大的單一集中轉(zhuǎn)動慣量,并且假定非承載齒輪對所受阻滯力矩為定值,這與實際情況不符.Robinette等[4-5]針對裝有某款手動變速器的傳動系統(tǒng),采用集中質(zhì)量的系統(tǒng)動力學建模方法,建立了空套齒輪的阻滯力矩力學模型,并通過與仿真模型對應(yīng)的試驗裝置驗證了仿真模型的有效性.但該模型利用時域內(nèi)階梯型遞增的轉(zhuǎn)矩輸入模擬整個傳動系統(tǒng)的動力輸入,而且模型中未考慮輪胎特性和整車平動質(zhì)量,無法真實還原實車在特定工況下的敲擊現(xiàn)象.馬小英等[6]利用AVL-Timing Drive商業(yè)軟件搭建變速器系統(tǒng)模型,在齒輪阻滯力矩中考慮攪油損失力矩和軸承摩擦力矩,仿真分析了變速器在怠速工況、五檔同步器不同安裝位置(輸入軸或者輸出軸)下各擋位齒輪對嚙合力大小,提出五檔同步器安裝在輸入軸上可減小怠速工況下敲擊噪聲,但未能對齒輪敲擊強度主要影響因素進行系統(tǒng)分析.

    本文在前人研究工作的基礎(chǔ)之上,采用基于集中質(zhì)量的系統(tǒng)動力學建模方法,利用車輛基本參數(shù)和傳動系統(tǒng)各部件的等效轉(zhuǎn)動慣量、扭轉(zhuǎn)剛度和阻尼扭轉(zhuǎn)振動特性參數(shù),建立了傳動系統(tǒng)非線性扭轉(zhuǎn)振動和車輛縱向受力耦合模型.利用該耦合模型,再現(xiàn)了變速器工作在一擋工況時各非承載齒輪對齒輪敲擊現(xiàn)象,并探討了非承載齒輪對齒側(cè)間隙大小和被動齒輪等效轉(zhuǎn)動慣量大小對敲擊強度的影響規(guī)律,為進一步研究齒輪敲擊現(xiàn)象奠定了理論基礎(chǔ).

    1 動力傳動系統(tǒng)部件建模

    1.1 發(fā)動機建模

    在動力傳動系統(tǒng)中,發(fā)動機是系統(tǒng)主要激勵源.對于四缸四沖程發(fā)動機而言,發(fā)動機不平衡往復慣性力及氣體壓力的點火階次造成發(fā)動機輸出轉(zhuǎn)矩波動.發(fā)動機輸出力矩Te為平均輸出力矩Tm和波動力矩Tp之和[7],如下所示:

    式中:Tpn為n次諧波分量的幅值,N·m;φn為n次諧波分量的相位,rad;N為建模中考慮的諧波分量的個數(shù),一般取N=5;t為時間,s;ωe為發(fā)動機曲軸轉(zhuǎn)動角頻率,rad·s-1;ωp為四缸四沖程發(fā)動機點火頻率,rad·s-1.

    1.2 離合器扭轉(zhuǎn)減振器特性建模

    離合器扭轉(zhuǎn)減振器作為動力傳動系統(tǒng)的重要部件,用來避免車輛傳動系統(tǒng)共振和緩和所受的沖擊載荷,其結(jié)構(gòu)原理如圖1所示.當離合器主、從動摩擦盤完全接合后,傳動系統(tǒng)輸入扭矩波動或者負載力矩變化時,從動盤體和從動盤轂之間將產(chǎn)生相對往復擺動,依靠兩者之間減振器彈簧的彈性和阻尼特性改變傳動系統(tǒng)扭轉(zhuǎn)振動特性,消耗扭轉(zhuǎn)振動能量.

    圖1 扭轉(zhuǎn)減振器結(jié)構(gòu)原理圖Fig.1 Schematic diagram of torsional vibration damper

    變速器齒輪敲擊屬于動力傳動系統(tǒng)扭轉(zhuǎn)振動范疇,對于實車出現(xiàn)變速器齒輪敲擊現(xiàn)象,多通過合理調(diào)節(jié)離合器扭轉(zhuǎn)減振器特性參數(shù)實現(xiàn)對變速器齒輪敲擊振動噪聲的有效控制[8],因此離合器扭轉(zhuǎn)減振特性的建模至關(guān)重要,本文選取實車采用的某款離合器扭轉(zhuǎn)減振器.試驗測得的扭轉(zhuǎn)特性數(shù)據(jù)經(jīng)過曲線擬合后得到的扭轉(zhuǎn)減振器特性曲線如圖2所示[].

    圖2 考慮遲滯特性的離合器扭轉(zhuǎn)減振器特性曲線Fig.2 Characteristics curve of clutch torsional vibration damper considering hysteresis characteristics

    1.3 考慮齒側(cè)間隙的單對齒輪對建模

    對于變速器中任意一對非承載齒輪副的力學模型可簡化為如圖3所示.小齒輪為敲擊主動齒輪,大齒輪為敲擊被動齒輪,此處將主動齒輪時域內(nèi)的運動輸入作為單對齒輪系統(tǒng)的激勵,因此圖中主動齒輪處不施加驅(qū)動力矩.單對齒輪系統(tǒng)滿足如下力學關(guān)系[10]:

    式中:X為沿基圓公切線上的相對位移,m;θp、θg分別為主動齒輪、被動齒輪扭轉(zhuǎn)角度,rad;Rp、Rg分別為主動齒輪、被動齒輪基圓半徑,m;Jg為被動齒輪等效轉(zhuǎn)動慣量,kg·m2;km為齒輪副平均嚙合剛度,N·m-1;cm為齒輪副平均嚙合阻尼,N·s·m-1;TD為敲擊被動齒輪所受阻滯力矩,N·m;Fm為齒輪對敲擊力,N.其中,齒側(cè)間隙函數(shù)f(X)滿足如下關(guān)系:

    式中:b為齒輪副基圓嚙合線上的齒側(cè)間隙,m.

    圖3 單對齒輪副敲擊簡化模型Fig.3 Simplified model for rattling gear pair

    1.4 空套齒輪拖曳力矩計算

    由第1.3節(jié)可知,合理計算敲擊被動齒輪(空套齒輪)的拖曳力矩可準確反映單對齒輪對敲擊力學模型.對于不同擋位的空套齒輪,阻滯力矩TD可能包括滾子軸承損失力矩Trb、油膜剪切損失力矩Tsh和攪油損失力矩Tch.

    對于滾子軸承損失力矩Trb,可通過下式進行估計[]:

    式中:nr為空套齒輪和軸段之間(軸承內(nèi)外圈)相對轉(zhuǎn)速,r·min-1;dm為軸承內(nèi)外圈平均直徑,m;f0為潤滑因子,與潤滑油類型和潤滑方式有關(guān);ν為潤滑油運動黏度,m2·s-1.

    油膜剪切損失力矩Tsh產(chǎn)生于空套齒輪與相鄰軸段、同步器齒環(huán)之間,假設(shè)其間油膜流動為庫埃特流動方式,則此時油膜剪切損失力矩[5]

    式中:η為潤滑油絕對黏度,N·s·m-2;Lg為齒輪長度,m;R為輪齒節(jié)圓半徑,m;Δn為空套齒輪與相鄰軸段、同步器齒環(huán)之間的相對轉(zhuǎn)速,r·min-1;j為軸承徑向間隙,m.

    對于攪油齒輪,其浸油幾何關(guān)系簡圖如圖4所示,攪油損失力矩[12]

    式中:ρ為潤滑油密度,kg·m-3;ωg為齒輪旋轉(zhuǎn)角速度,rad·s-1;Sm為齒輪浸油面積,m2;Cm為攪油損失力矩系數(shù).

    圖4 攪油齒輪浸油幾何關(guān)系簡圖Fig.4 Geometry relationship of churning oil-immersed gear

    齒輪浸油面積Sm包括齒輪外側(cè)面浸油面積Sm,l和齒輪輪齒浸油面積Sm,t之和,即

    其中,

    式中:B為齒輪輪齒齒寬,m;ht為齒輪輪齒高度,m;α為齒輪輪齒壓力角,rad;Z為齒輪輪齒齒數(shù);h為齒輪浸油深度,m.

    對于攪油損失力矩系數(shù)C,計算式為m

    式中:V0為箱內(nèi)油體積,L;Fr為弗勞德數(shù);Re為雷諾數(shù);g為重力加速度,m·s-2.

    1.5 軸段扭轉(zhuǎn)剛度計算

    對于變速器輸入、輸出軸,按照各擋位齒輪所在位置分成多段,等截面軸段扭轉(zhuǎn)剛度

    式中:E為軸段材料拉壓彈性模量,GPa;IP為軸段截面極慣性矩,m4;l為等截面軸段長度,m;μ為軸段材料泊松比.

    實際簡化時,將變速器輸入軸或者輸出軸相鄰齒輪間的軸段等效為一個扭轉(zhuǎn)剛度,此時需要等效的軸段多為多截面軸段,此時多截面等效剛度

    式中:ki為第i個軸段等效剛度,N·m·rad-1;NL為軸段個數(shù).

    1.6 基于分布參數(shù)的LuGre輪胎模型

    本文采用基于分布參數(shù)的LuGre輪胎模型,其受力及運動簡圖如圖5所示[13].

    圖5 基于分布參數(shù)的LuGre輪胎模型Fig.5 LuGre tire model based on distributed parameters

    設(shè)輪胎上刷毛的平均變形量為z,輪胎相對地面的速度為vr,則毛刷變形的微分方程如下所示[14]:

    式中:σ0為縱向等效剛度系數(shù),m-1;σ1為縱向等效阻尼系數(shù),s·m-1;σ2為摩擦面相對黏性阻尼系數(shù),s·m-1;μc為動摩擦因數(shù);μs為靜摩擦因數(shù);vs為斯特列貝格相對速度,m·s-1;λ為斯特列貝格效應(yīng)指數(shù);Fz為地面支持力,N;fn(ζ)為沿輪胎印記縱向壓力分布密度函數(shù);n為接地印記沿著ζ軸方向的積分個數(shù);s為輪胎滑移率;r為輪胎滾動半徑,m.

    2 手動變速器結(jié)構(gòu)形式及建模思路

    本文以裝備于某車輛的五擋手動變速器為例,變速器主要部件包括輸入軸、輸出軸、齒輪對和同步器總成等,采用兩軸式設(shè)計(不包含倒擋軸),包括5個前進擋和1個倒擋,其具體結(jié)構(gòu)形式如圖7中所示.當變速器掛入某擋位時,相應(yīng)的齒輪對傳遞動力,稱為承載齒輪對,其余不傳遞動力的齒輪對為非承載齒輪對.

    手動變速器齒輪敲擊現(xiàn)象是動力傳動系統(tǒng)扭轉(zhuǎn)振動綜合作用的結(jié)果,針對裝備手動變速器的前置前驅(qū)車輛進行動力傳動系統(tǒng)建模,其總體建模思路如圖6所示,車輛動力傳動系統(tǒng)非線性扭轉(zhuǎn)振動分為主線扭轉(zhuǎn)振動和敲擊扭轉(zhuǎn)振動.研究學者已經(jīng)廣泛認為空套齒輪的敲擊扭轉(zhuǎn)振動對于主線扭轉(zhuǎn)振動的影響很?。?-4],因此對于各非承載齒輪對,將主線扭轉(zhuǎn)振動上的敲擊主動齒輪運動作為敲擊扭轉(zhuǎn)振動系統(tǒng)的激勵,而不考慮被動齒輪對于主動齒輪運動的影響.

    圖6 整車動力傳動系統(tǒng)建模結(jié)構(gòu)圖Fig.6 Structure of vehicle powertrain system modelling

    3 傳動系統(tǒng)扭轉(zhuǎn)振動力學模型

    以車輛行駛在水平道路上、工作擋位為一擋工況,仿真分析手動變速器齒輪敲擊現(xiàn)象,建?;谝韵录僭O(shè):

    (1)為離合器完全接合后仿真工況,不涉及離合器主、從動盤接合的動力傳動過程.對主線扭轉(zhuǎn)振動上的變速器和主減速器內(nèi)傳遞動力齒輪對,輸入軸二擋齒輪,輸出軸三、四、五擋齒輪和同步器總成,忽略攪油損失力矩或者油膜剪切損失力矩對其扭轉(zhuǎn)振動的影響.

    (2)忽略變速器輸入軸和輸出軸軸承、主減速器軸承動態(tài)力學特性對于變速器內(nèi)部齒輪系統(tǒng)扭轉(zhuǎn)振動特性的影響,忽略差速器特性和左右半軸不同扭轉(zhuǎn)剛度對傳動系統(tǒng)扭轉(zhuǎn)振動的影響,忽略車輛前后軸左、右車輪垂直載荷的差異,忽略變速器內(nèi)部部件系統(tǒng)運動帶來的潤滑油溫升.

    3.1 主線扭轉(zhuǎn)振動建模

    對于主線扭轉(zhuǎn)振動,建立如圖7所示的傳動系統(tǒng)扭轉(zhuǎn)振動力學模型和如圖8所示的整車縱向運動和受力分析模型.

    圖7 動力傳動系統(tǒng)扭轉(zhuǎn)振動模型Fig.7 Driveline torsional vibration model

    圖8 整車縱向運動及受力分析圖Fig.8 Vehicle longitudinal motion and force analysis

    圖7中,Jf為飛輪和離合器主從動摩擦盤轉(zhuǎn)動慣量,kg·m2;JC為離合器從動盤轂轉(zhuǎn)動慣量,kg· m2;JP1~JP5分別為輸入軸一~五擋齒輪轉(zhuǎn)動慣量,kg·m2;JG1S1為輸出軸一擋齒輪和相應(yīng)同步器總成轉(zhuǎn)動慣量之和,kg·m2;JP0為輸入軸倒擋齒輪轉(zhuǎn)動慣量,kg·m2;JG2~JG5分別為輸出軸二~五擋齒輪轉(zhuǎn)動慣量,kg·m2;JS2為三、四擋同步器總成轉(zhuǎn)動慣量,kg·m2;JS3為五擋同步器總成轉(zhuǎn)動慣量,kg· m2;JFD1為主減速器主動齒輪轉(zhuǎn)動慣量,kg·m2;JFD2為主減速器被動齒輪轉(zhuǎn)動慣量,kg·m2;JT為差速器、半軸和車輪的等效轉(zhuǎn)動慣量,kg·m2;k1~k5、c1~c5分別為一~五檔齒輪副平均嚙合剛度和嚙合阻尼,單位分別為N·m-1和N·s·m-1;kf、cf分別為主減速器齒輪副平均嚙合剛度和嚙合阻尼,單位分別為N·m-1和N·s·m-1;b1~b5分別為一~五擋齒輪副齒側(cè)間隙,m;bf為主減速器齒輪副齒側(cè)間隙,m;TD2~TD5分別為二~五擋被動齒輪所受阻滯力矩,N·m;Tc為離合器扭轉(zhuǎn)減振器傳遞扭矩,N·m;Tff為輪胎所受地面作用力矩,N·m.

    式中:CD為空氣阻力系數(shù);A為車輛迎風面積,m2;ρa為空氣密度,kg·m-3.

    此時,在車輛縱向加速度影響下前、后軸的載荷分配為

    根據(jù)牛頓第二定律,建立分析變速器齒輪敲擊主線扭轉(zhuǎn)振動的數(shù)學模型,即

    式中:J為慣量矩陣;K為剛度矩陣;C為阻尼矩陣;θ為(角)位移矩陣;T為作用扭矩(力)矩陣;θi(i=f,C,…,T)為對應(yīng)轉(zhuǎn)動慣量Ji(i=f,C,…,T)的角位移,rad;m為車輛整備質(zhì)量,kg;xs為車輛平移運動位移,m;Fxf為驅(qū)動前輪所受地面縱向力,N;Fxr為從動后輪所受地面縱向力,N;Fw為車輛所受空氣阻力,N.

    圖8中,xv為相對速度,在無風時即汽車的行駛速度,m·s-1;Wg為車輛質(zhì)心;La、Lb分別為車輛質(zhì)心至前、后軸距離,m;H為車輛質(zhì)心高度,m;G為車輛整備重力,N;Tt為驅(qū)動前輪等效驅(qū)動力矩,N· m;Tjw1為驅(qū)動前輪慣性力矩,N·m;Fzf為車輛前軸垂直載荷,N;Tjw2為從動后輪慣性力矩,N·m;Tfr為從動后輪受到地面作用力矩,N·m;Fzr為車輛后軸垂直載荷,N.

    如圖8中所示,車輛所受空氣阻力

    式中:a為車輛運動加速度,m·s-2;g為重力加速度,m·s-2.

    對于驅(qū)動前輪進行受力分析,扭轉(zhuǎn)振動方程為

    式中:k、c分別為主減速器被動齒輪和半軸、車輪等效轉(zhuǎn)動慣量間的等效扭轉(zhuǎn)剛度和扭轉(zhuǎn)阻尼,單位分別為N·m·rad-1和N·m·s·rad-1;f為滾動阻力系數(shù).

    對于從動后輪進行受力分析,扭轉(zhuǎn)振動方程為

    式中:JTr為從動后輪的等效轉(zhuǎn)動慣量,kg·m2;θTr為后輪轉(zhuǎn)動角位移,rad.輪胎模型采用基于分布參數(shù)的LuGre模型,地面縱向力Fxf、Fxr利用第1.6節(jié)中式(20)~(24)得到.

    3.2 敲擊扭轉(zhuǎn)振動建模

    對于敲擊扭轉(zhuǎn)振動,各非承載齒輪對沿基圓公切線上的相對位移大小分別為

    式中:RP2、RG2分別為輸入、輸出軸二擋齒輪基圓半徑,m;RP3、RG3分別為輸入、輸出軸三擋齒輪基圓半徑,m;RP4、RG4分別為輸入、輸出軸四擋齒輪基圓半徑,m;RP5、RG5分別為輸入、輸出軸五擋齒輪基圓半徑,m.

    對于二擋齒輪副,有

    式中:Fm2為二擋齒輪對敲擊力,N;Trb2為二擋被動齒輪所受軸承損失力矩,N·m;Tsh2為二擋被動齒輪所受油膜剪切損失力矩,N·m;Tch2為二擋被動齒輪所受攪油損失力矩,N·m.

    對于三擋齒輪副,有

    式中:Fm3為三擋齒輪對敲擊力,N;Trb3為三擋被動齒輪所受軸承損失力矩,N·m.

    對于四擋齒輪副,有

    式中:Fm4為四擋齒輪對敲擊力,N;Trb4為四擋被動齒輪所受軸承損失力矩,N·m.

    對于五擋齒輪副,有

    式中:Fm5為五擋齒輪對敲擊力,N;Trb5為五擋被動齒輪所受軸承損失力矩,N·m.

    4 齒輪敲擊結(jié)果及其影響因素分析

    根據(jù)整車廠家反饋,當該款變速器工作在一擋、車輛處于離合器完全接合后起步加速的初始工況下,實車存在較嚴重的齒輪敲擊現(xiàn)象.以下仿真再現(xiàn)該初始工況下各非承載齒輪對敲擊現(xiàn)象,并分析主要影響因素對于敲擊強度的影響規(guī)律.

    4.1 實車狀態(tài)及各擋位敲擊情況再現(xiàn)

    通過仿真再現(xiàn)實車工況下的發(fā)動機飛輪轉(zhuǎn)速和車速變化,分別如圖9和10所示,此時對應(yīng)的各非承載齒輪對敲擊力大小如圖11所示.

    圖9 發(fā)動機飛輪轉(zhuǎn)速隨時間變化曲線Fig.9 Engine flywheel speed curve in time domain

    圖10 車速隨時間變化曲線Fig.10 Vehicle speed curve in time domain

    圖11 各非承載齒輪對敲擊力隨時間變化曲線Fig.11 Unloaded gear pairs rattle force in time domain

    由圖9和10可知,此時發(fā)動機轉(zhuǎn)速和車速處于加速過程,發(fā)動機飛輪轉(zhuǎn)速存在波動,此轉(zhuǎn)速波動通過離合器扭轉(zhuǎn)減振器傳至變速器輸入軸而引起變速器內(nèi)發(fā)生齒輪敲擊現(xiàn)象.同時,由圖11可知,對于二擋齒輪對,隨著發(fā)動機轉(zhuǎn)速升高,敲擊強度有減小趨勢,主要是因為二擋被動齒輪隨轉(zhuǎn)速升高,攪油阻滯力矩增加,導致敲擊強度減小,再現(xiàn)工況末端二擋敲擊力大小趨于穩(wěn)定;對于三擋、四擋和五擋,隨著發(fā)動機轉(zhuǎn)速升高,敲擊強度增大,再現(xiàn)工況末端三擋、四擋和五擋敲擊力變化不大.

    圖12 二擋齒輪對間隙影響敲擊力圖Fig.12 Rattle force of the second gear pair influenced by gear backlash

    4.2 齒輪對齒側(cè)間隙大小對于敲擊強度影響規(guī)律分析

    由第1.3節(jié)單對齒輪對敲擊模型可知,齒側(cè)間隙作為數(shù)學模型中主要參數(shù),直接影響敲擊結(jié)果.由該款手動變速器廠家反饋,各擋位齒輪對齒側(cè)間隙范圍為0.06~0.12mm.現(xiàn)分析不同齒側(cè)間隙大小對于各非承載齒輪對敲擊強度影響規(guī)律,結(jié)果如圖12~15所示.

    圖13 三擋齒輪對間隙影響敲擊力圖Fig.13 Rattle force of the third gear pair influenced by gear backlash

    當齒側(cè)間隙b取較小值時,各敲擊齒輪對敲擊頻次較高,結(jié)合式(4)~(8)可知,當齒側(cè)間隙越小,敲擊被動齒輪一旦與主動齒輪脫嚙,敲擊被動齒輪齒背就更容易與下一個主動齒輪齒面產(chǎn)生敲擊,導致敲擊的頻次增多.當齒側(cè)間隙由小到大時,敲擊力大小總體呈現(xiàn)出先增大后減小的變化趨勢,這是因為當齒側(cè)間隙增加時,敲擊被動齒輪所受到的阻滯力矩作用時間相對較長,導致敲擊前主被動齒輪相對速度增大,敲擊力反而增大,當齒側(cè)間隙繼續(xù)增加時,阻滯力矩作用時間更長,敲擊被動齒輪未能發(fā)生齒背敲擊而繼續(xù)發(fā)生齒面敲擊.三擋較四擋、五擋輪齒敲擊的頻次高,這是因為在齒側(cè)間隙相同的情況下,變速器處于一擋工況時,三擋、四擋和五擋角速度大小基本一致,齒輪對數(shù)多的三擋齒輪產(chǎn)生敲擊的可能性越大.通過以上分析可知,合理設(shè)置各擋位齒側(cè)間隙對于有效抑制齒輪敲擊現(xiàn)象具有至關(guān)重要的作用,但是一旦實車工況下出現(xiàn)敲擊現(xiàn)象,變速器內(nèi)其他部件參數(shù)已經(jīng)確定,很難調(diào)整齒側(cè)間隙以抑制敲擊現(xiàn)象,因此應(yīng)在變速器正向設(shè)計中考慮多方 面因素確定合適的各擋位齒側(cè)間隙大小.

    圖14 四擋齒輪對間隙影響敲擊力圖Fig.14 Rattle force of the fourth gear pair influenced by gear backlash

    圖15 五擋齒輪對間隙影響敲擊力圖Fig.15 Rattle force of the fifth gear pair influenced by gear backlash

    4.3 非承載齒輪轉(zhuǎn)動慣量對于敲擊強度影響規(guī)律分析

    圖16 二擋被動齒輪轉(zhuǎn)動慣量影響敲擊力圖Fig.16 Rattle force of the second gear pair influenced by driven gear torsional inertia

    由第1.3節(jié)單對齒輪對敲擊模型可知,被動齒輪轉(zhuǎn)動慣量作為數(shù)學模型中主要參數(shù),也可直接影響敲擊結(jié)果.以下分析非承載齒輪對被動齒輪轉(zhuǎn)動慣量大小對于敲擊強度的影響規(guī)律,結(jié)果如圖16~19所示.

    圖17 三擋被動齒輪轉(zhuǎn)動慣量影響敲擊力圖Fig.17 Rattle force of the third gear pair influenced by driven gear torsional inertia

    圖18 四擋被動齒輪轉(zhuǎn)動慣量影響敲擊力圖Fig.18 Rattle force of the fourth gear pair influenced by driven gear torsional inertia

    圖19 五擋被動齒輪轉(zhuǎn)動慣量影響敲擊力圖Fig.19 Rattle force of the fifth gear pair influenced by driven gear torsional inertia

    由圖16結(jié)果可知,對于二擋齒輪對,被動齒輪轉(zhuǎn)動慣量JG2由6.05×10-4kg·m2變化至7.25× 10-4kg·m2時,敲擊頻次逐漸增加,然而對于敲擊力大小影響規(guī)律不明顯.由圖17結(jié)果可知,對于三擋齒輪對,當被動齒輪轉(zhuǎn)動慣量JP3由1.35×10-4kg·m2變化至2.55×10-4kg·m2,敲擊頻次逐漸增加,對于敲擊力大小影響規(guī)律不明顯.由圖18和19結(jié)果可知,對于給定的四擋和五擋被動齒輪轉(zhuǎn)動慣量變化范圍,隨著被動齒輪轉(zhuǎn)動慣量的增加,敲擊力和敲擊頻次均增強.各擋敲擊頻次隨被動齒輪轉(zhuǎn)動慣量增加而增加,是因為當主動齒輪與空套被動齒輪發(fā)生敲擊后,當空套被動齒輪轉(zhuǎn)動慣量增加時,空套被動齒輪周向轉(zhuǎn)動慣量增加,導致敲擊被動齒輪發(fā)生敲擊彈開,然后很快再次與主動齒輪發(fā)生敲擊現(xiàn)象,從而使得敲擊頻次增高.該結(jié)果與文獻[4]利用臺架試驗得到的結(jié)論保持一致.因此,在變速器主動齒輪設(shè)計過程中,應(yīng)充分考慮各擋位的轉(zhuǎn)動慣量設(shè)置(反映在各擋位齒輪具體結(jié)構(gòu)上).

    5 結(jié)論

    (1)本文利用車輛動力傳動系統(tǒng)非線性扭轉(zhuǎn)振動和縱向平移運動耦合模型,再現(xiàn)了車輛處于一擋加速工況時機械變速器各非承載齒輪對的敲擊情況.結(jié)果表明,隨著發(fā)動機轉(zhuǎn)速升高,二擋非承載齒輪對敲擊力逐漸減小,最后趨于穩(wěn)定,然而三擋、四擋和五擋非承載齒輪對敲擊力逐漸增大,最后趨于穩(wěn)定.

    (2)分析了非承載齒輪對齒側(cè)間隙大小和被動齒輪等效轉(zhuǎn)動慣量大小對敲擊強度的影響規(guī)律.結(jié)果表明,當齒側(cè)間隙位于0.06~0.12mm范圍時,隨著齒側(cè)間隙增加,各擋位齒輪對敲擊力先增加后減小,齒輪對敲擊頻次單調(diào)減小,并且三擋較四擋、五擋輪齒敲擊的頻次高.對于給定范圍內(nèi)的各擋位被動齒輪轉(zhuǎn)動慣量,隨著齒輪等效轉(zhuǎn)動慣量增加,各擋位齒輪對敲擊頻次單調(diào)增加,然而對于敲擊力大小,二擋和三擋影響規(guī)律不明顯,四擋和五擋敲擊力單調(diào)增加.

    [1]Naunheimer H,Bertsche B,Ryborz J,et al.Automotive transmissions:fundamentals,selection,design and application [M].Berlin:Springer Science &Business Media,2010.

    [2]Wang M Y,Manoj R.汽車手動變速箱的齒輪敲擊模擬和分析[J].傳動技術(shù),2002,16(4):27.Wang M Y,Manoj R.Gear rattle modelling and analysis for automotive manual transmissions[J].Drive System Technique,2002,16(4):27.

    [3]Wang M Y,Zhao W,Manoj R.Numerical modelling and analysis of automotive transmission rattle[J].Journal of Vibration and Control,2002,8(7):921.

    [4]Robinette D,Beikmann R S,Piorkowski P,et al. Characterizing the onset of manual transmission gear rattle.Part I:experimental results[J].SAE International Journal of Passenger Cars:Mechanical Systems,2009,2(1):1352.

    [5]Robinette D,Beikmann R S,Piorkowski P,et al.Characterizing the onset of manual transmission gear rattle.Part II:analytical results[J].SAE International Journal of Passenger Cars:Mechanical Systems,2009,2(1):1365.

    [6]馬小英,彭國民,余波.汽車手動變速器怠速敲擊噪聲研究及優(yōu)化[J].汽車工程學報,2011,1(Suppl.1):44.MA Xiaoying,PENG Guomin,YU Bo.Gear rattle research and analysis for automotive manual transmissions[J].Chinese Journal of Automotive Engineering,2011,1(Suppl.1):44.

    [7]Yoon J Y,Singh R.Effect of the multi-staged clutch damper characteristics on the transmission gear rattle under two engine conditions[J].Proceedings of the Institution of Mechanical Engineers,Part D:Journal of Automobile Engineering,2013,227(9):1273.

    [8]Prasad J S,Damodar N C,Naidu T S.Clutch hysteresis maximization for elimination of gear rattle in a passenger bus [R].Maharashtra:SAE Technical Paper,2013.

    [9]欒文博,吳光強.基于六自由度傳動系模型的自動變速器起步仿真[J].農(nóng)業(yè)機械學報,2013,44(4):1.LUAN Wenbo,WU Guangqiang.Launch simulation for AMT based on 6-DOF model of vehicle powertrain[J].Transactions of the Chinese Society for Agricultural Machinery,2013,44 (4):1.

    [10]Ma X,Zhan Z,Peng G,et al.Rattle simulation analysis of a manual gearbox[C]//Proceedings of the FISITA 2012 World Automotive Congress.Berlin:Springer Berlin Heidelberg,2013:219-233.

    [11]劉澤九.滾動軸承應(yīng)用手冊[M].北京:機械工業(yè)出版社,2006.LIU Zejiu.Rolling bearing application manual[M].Beijing:China Machine Press,2006.

    [12]Changenet C,Oviedo-Marlot X,Velex P.Power loss predictions in geared transmissions using thermal networksapplications to a six-speed manual gearbox[J].Journal of Mechanical Design,2006,128(3):618.

    [13]Canudas de Wit C,Tsiotras P.Dynamic tire friction models for vehicle traction control[C]//Proceedings of 38th IEEE Conference of Decision and Control.Phoenix:IEEE,1999,4 (3):3746-3751.

    [14]鄧堃,李開軍,何樂,等.用于車輛緊急制動仿真的動態(tài)輪胎模型[J].農(nóng)業(yè)機械學報,2007,38(12):11.DENG Kun,LI Kaijun,HE Le,et al.Dynamic tire model for vehicle emergency braking simulation[J].Transactions of the Chinese Society for Agricultural Machinery,2007,38(12):11.

    Modelling of Manual Transmission Gear Rattle Phenomenon and Its Impact Factors

    WU Huwei1,WU Guangqiang1,2
    (1.School of Automotive Studies,Tongji University,Shanghai 201804,China;2.Institute of Industrial Science,the University of Tokyo,Tokyo 153-8505,Japan)

    Based on basic parameters of vehicle and rotational parameters of real driveline system components,e.g.torsional inertia,stiffness and damping,the coupling model of nonlinear torsional vibration of vehicle powertrain and longitudinal translation was established.With the coupling model,each unloaded gear pair rattling phenomenon was reproduced when the vehicle in gear one was accelerating.Then,gear rattle sensitivity was analyzed for different impact factors.It is shown that appropriate gear backlash and torsional inertia could inhibit the gear rattle phenomenon of manual transmission.

    LuGre tire model;nonlinear torsional vibration;coupling model;gear backlash;equivalent torsional inertia成部分,其噪聲包括變速器嘯叫噪聲、齒輪敲擊噪聲、換擋噪聲、軸承噪聲和Clonk噪聲,是汽車振動噪聲的主要來源[1].對于變速器齒輪系統(tǒng)來說,根據(jù)不同的振動狀態(tài),噪聲主要分為齒輪敲擊噪聲和嘯叫噪聲2類,其中齒輪敲擊現(xiàn)象是發(fā)生在變速器中常嚙合非承載齒輪副.當變速器工作在某擋位時,常嚙合非承載齒輪對中存在單個齒輪空套在變速器軸上的情況,由于輪齒間存在齒側(cè)間隙,若非承載齒輪副周向運動位移差幅值超過齒側(cè)間隙大小,則將產(chǎn)生敲擊現(xiàn)象.變速器齒輪敲擊噪聲嚴重影響乘員乘坐舒適性,直接影響車輛品牌形象,因此變速器敲擊現(xiàn)象產(chǎn)生機理研究至關(guān)重要.

    U463.21

    A

    0253-374X(2016)01-0136-10

    10.11908/j.issn.0253-374x.2016.01.020

    2015-01-17

    國家自然科學基金(51175379)

    吳虎威(1987—),男,博士生,主要研究方向為車輛動力傳動系統(tǒng)振動噪聲.E-mail:1133054@#edu.cn

    吳光強(1965—),男,教授,博士生導師,工學博士,主要研究方向為車輛現(xiàn)代化設(shè)計理論及方法以及車輛動力學及其控制.E-mail:wuguangqiang@#edu.cn

    猜你喜歡
    振動模型
    一半模型
    振動的思考
    科學大眾(2023年17期)2023-10-26 07:39:14
    噴水推進高速艇尾部振動響應(yīng)分析
    重要模型『一線三等角』
    This “Singing Highway”plays music
    重尾非線性自回歸模型自加權(quán)M-估計的漸近分布
    振動攪拌 震動創(chuàng)新
    中國公路(2017年18期)2018-01-23 03:00:38
    中立型Emden-Fowler微分方程的振動性
    3D打印中的模型分割與打包
    FLUKA幾何模型到CAD幾何模型轉(zhuǎn)換方法初步研究
    免费黄频网站在线观看国产| www.色视频.com| 精品久久久久久久久av| 曰老女人黄片| 大又大粗又爽又黄少妇毛片口| 国产精品一区二区三区四区免费观看| 亚洲欧美精品专区久久| 中国国产av一级| 曰老女人黄片| 国产精品国产三级国产专区5o| 午夜影院在线不卡| 成人毛片60女人毛片免费| 国产极品粉嫩免费观看在线 | 观看av在线不卡| 男女国产视频网站| 国产日韩欧美视频二区| 精品少妇黑人巨大在线播放| 日本午夜av视频| 青春草亚洲视频在线观看| 国产精品嫩草影院av在线观看| 精品亚洲乱码少妇综合久久| 免费观看av网站的网址| 亚洲伊人久久精品综合| 国产在视频线精品| 国产免费又黄又爽又色| 成人毛片a级毛片在线播放| 18禁裸乳无遮挡动漫免费视频| 精品卡一卡二卡四卡免费| 中文字幕制服av| 十八禁网站网址无遮挡 | 久久久久久久国产电影| 大香蕉久久网| 美女国产视频在线观看| 99久国产av精品国产电影| 国产亚洲av片在线观看秒播厂| 久久国产精品大桥未久av | 国产高清三级在线| 亚洲国产精品一区二区三区在线| 日本猛色少妇xxxxx猛交久久| 亚洲四区av| 在线观看人妻少妇| 人妻夜夜爽99麻豆av| 97超视频在线观看视频| 丝袜脚勾引网站| 欧美精品高潮呻吟av久久| 最近中文字幕2019免费版| 少妇人妻久久综合中文| 欧美 亚洲 国产 日韩一| 乱码一卡2卡4卡精品| 欧美区成人在线视频| 日韩制服骚丝袜av| 亚洲四区av| 国产日韩欧美在线精品| 久久免费观看电影| 在线观看美女被高潮喷水网站| 丰满饥渴人妻一区二区三| 欧美老熟妇乱子伦牲交| 日产精品乱码卡一卡2卡三| av一本久久久久| 国产淫片久久久久久久久| 日韩视频在线欧美| 午夜福利网站1000一区二区三区| 国产高清国产精品国产三级| 18禁动态无遮挡网站| 一个人看视频在线观看www免费| 日本与韩国留学比较| 9色porny在线观看| 久久久久久久国产电影| 永久网站在线| 伦理电影免费视频| 如何舔出高潮| 一级毛片黄色毛片免费观看视频| 97在线视频观看| 狂野欧美白嫩少妇大欣赏| 国产在线男女| 国产 一区精品| 久久影院123| 欧美 亚洲 国产 日韩一| 少妇 在线观看| 国产视频内射| 一本久久精品| 人妻人人澡人人爽人人| 婷婷色综合www| 最近的中文字幕免费完整| 亚洲美女搞黄在线观看| 赤兔流量卡办理| 欧美日韩综合久久久久久| 亚洲真实伦在线观看| 一本—道久久a久久精品蜜桃钙片| 国产综合精华液| 亚洲精品日本国产第一区| 晚上一个人看的免费电影| 男人添女人高潮全过程视频| 久久久国产精品麻豆| 亚洲av中文av极速乱| 亚洲精品国产成人久久av| 麻豆乱淫一区二区| 亚洲经典国产精华液单| 九色成人免费人妻av| 水蜜桃什么品种好| 熟女av电影| 深夜a级毛片| 国产黄片视频在线免费观看| 一个人免费看片子| 男人舔奶头视频| 一级二级三级毛片免费看| 亚洲精品456在线播放app| 久久久久国产网址| 国产男人的电影天堂91| 色视频www国产| 久久久久国产网址| 一级片'在线观看视频| 久久精品夜色国产| 三级经典国产精品| 国产一区亚洲一区在线观看| 在线观看www视频免费| 亚洲在久久综合| 日本wwww免费看| 国语对白做爰xxxⅹ性视频网站| 精品久久久噜噜| 久久久久久久久久久久大奶| 啦啦啦中文免费视频观看日本| 亚洲人成网站在线播| 最近中文字幕2019免费版| 777米奇影视久久| 亚洲人与动物交配视频| 色94色欧美一区二区| 女人精品久久久久毛片| 国产精品国产三级国产专区5o| 欧美日本中文国产一区发布| 各种免费的搞黄视频| 久久久久精品性色| 久久久久久久久久成人| av免费观看日本| 精品卡一卡二卡四卡免费| 午夜影院在线不卡| 哪个播放器可以免费观看大片| 亚洲成人手机| 欧美一级a爱片免费观看看| 亚洲精品亚洲一区二区| 在线观看免费高清a一片| 一级,二级,三级黄色视频| av天堂久久9| 晚上一个人看的免费电影| av天堂中文字幕网| 亚洲精品国产色婷婷电影| 99热6这里只有精品| 人妻制服诱惑在线中文字幕| 啦啦啦在线观看免费高清www| 免费在线观看成人毛片| 国产在线免费精品| 蜜臀久久99精品久久宅男| 国产男女内射视频| 嘟嘟电影网在线观看| 久久婷婷青草| 亚洲美女黄色视频免费看| 大片免费播放器 马上看| 日韩视频在线欧美| 国产极品天堂在线| 日韩欧美精品免费久久| 另类亚洲欧美激情| 亚洲国产精品专区欧美| 少妇人妻精品综合一区二区| 黑丝袜美女国产一区| 国产亚洲欧美精品永久| 亚洲国产成人一精品久久久| 男的添女的下面高潮视频| 欧美性感艳星| 亚洲av成人精品一区久久| 最新的欧美精品一区二区| 久久国产亚洲av麻豆专区| 欧美一级a爱片免费观看看| 国产亚洲5aaaaa淫片| 美女福利国产在线| 一级毛片aaaaaa免费看小| 欧美精品一区二区免费开放| 国产老妇伦熟女老妇高清| 日日摸夜夜添夜夜添av毛片| 色婷婷av一区二区三区视频| 色视频www国产| 欧美日韩综合久久久久久| 欧美国产精品一级二级三级 | 精品久久久久久电影网| 水蜜桃什么品种好| 国产高清不卡午夜福利| 亚洲精品国产成人久久av| 久久久久久久国产电影| 日日啪夜夜撸| 99久久精品国产国产毛片| 亚洲丝袜综合中文字幕| 亚洲av欧美aⅴ国产| 美女国产视频在线观看| 国产精品久久久久久精品古装| 欧美区成人在线视频| 国产高清国产精品国产三级| 一级爰片在线观看| 国内揄拍国产精品人妻在线| 51国产日韩欧美| 大香蕉久久网| 亚洲欧洲精品一区二区精品久久久 | 狂野欧美激情性xxxx在线观看| 在线亚洲精品国产二区图片欧美 | 一级片'在线观看视频| 成人二区视频| 高清午夜精品一区二区三区| 91在线精品国自产拍蜜月| 亚洲av成人精品一区久久| av播播在线观看一区| 午夜老司机福利剧场| 国产精品不卡视频一区二区| 又粗又硬又长又爽又黄的视频| 香蕉精品网在线| 高清黄色对白视频在线免费看 | 全区人妻精品视频| 亚洲欧美日韩卡通动漫| 亚洲不卡免费看| 国产精品久久久久久久电影| 美女视频免费永久观看网站| 欧美+日韩+精品| av.在线天堂| 2018国产大陆天天弄谢| 少妇 在线观看| 亚洲丝袜综合中文字幕| 日本-黄色视频高清免费观看| 精品久久久久久久久av| 婷婷色av中文字幕| 99re6热这里在线精品视频| 欧美日韩在线观看h| 一级二级三级毛片免费看| 一级av片app| 国产熟女欧美一区二区| 国产国拍精品亚洲av在线观看| 午夜激情福利司机影院| 十八禁网站网址无遮挡 | 国产免费视频播放在线视频| 麻豆精品久久久久久蜜桃| 80岁老熟妇乱子伦牲交| 国产熟女欧美一区二区| 九色成人免费人妻av| 国产亚洲一区二区精品| 国产伦理片在线播放av一区| av播播在线观看一区| 多毛熟女@视频| 在现免费观看毛片| 精品亚洲成a人片在线观看| 插逼视频在线观看| 免费不卡的大黄色大毛片视频在线观看| 最近手机中文字幕大全| 精品一区二区三卡| 亚洲av欧美aⅴ国产| 欧美少妇被猛烈插入视频| 麻豆精品久久久久久蜜桃| 不卡视频在线观看欧美| 视频区图区小说| 一级毛片我不卡| 成人美女网站在线观看视频| 国产永久视频网站| 视频中文字幕在线观看| 精品99又大又爽又粗少妇毛片| a级毛片在线看网站| 王馨瑶露胸无遮挡在线观看| 人妻制服诱惑在线中文字幕| 中国国产av一级| 夜夜爽夜夜爽视频| 日本vs欧美在线观看视频 | 国产成人免费观看mmmm| 自线自在国产av| 国产黄频视频在线观看| 亚洲欧美日韩另类电影网站| 久久热精品热| 下体分泌物呈黄色| 国产成人精品久久久久久| 卡戴珊不雅视频在线播放| 免费看不卡的av| 亚洲av福利一区| 久久精品国产a三级三级三级| 日本av免费视频播放| 丝袜脚勾引网站| 高清午夜精品一区二区三区| 日韩欧美精品免费久久| av天堂中文字幕网| 99久久综合免费| 一个人看视频在线观看www免费| 免费av中文字幕在线| 国产亚洲av片在线观看秒播厂| 五月伊人婷婷丁香| 国产免费视频播放在线视频| 免费看不卡的av| 啦啦啦在线观看免费高清www| 国产在线一区二区三区精| 国产一区二区三区av在线| 亚洲欧美一区二区三区国产| 国产成人精品无人区| videossex国产| 国产午夜精品一二区理论片| 国产精品免费大片| 男人添女人高潮全过程视频| 国产成人精品婷婷| 日韩视频在线欧美| 中文精品一卡2卡3卡4更新| 十八禁高潮呻吟视频 | 一级爰片在线观看| 亚洲精品乱久久久久久| 久久久国产欧美日韩av| 国产成人91sexporn| 纯流量卡能插随身wifi吗| 久久99精品国语久久久| 青春草国产在线视频| 国产永久视频网站| 国产成人aa在线观看| 在线看a的网站| 亚洲国产最新在线播放| 国产高清不卡午夜福利| av卡一久久| 亚洲欧美日韩东京热| 丰满迷人的少妇在线观看| 一级毛片aaaaaa免费看小| 亚州av有码| 亚洲精品一区蜜桃| 亚洲国产精品专区欧美| videos熟女内射| 亚洲av成人精品一二三区| 男人狂女人下面高潮的视频| 成人黄色视频免费在线看| 国产日韩欧美视频二区| 日韩av在线免费看完整版不卡| 五月天丁香电影| 免费观看av网站的网址| 国产精品熟女久久久久浪| 国产熟女欧美一区二区| 丰满少妇做爰视频| 美女脱内裤让男人舔精品视频| 男男h啪啪无遮挡| 国产探花极品一区二区| 亚洲情色 制服丝袜| 少妇的逼水好多| 欧美精品国产亚洲| 成年女人在线观看亚洲视频| 精品久久国产蜜桃| 成人午夜精彩视频在线观看| 狂野欧美激情性xxxx在线观看| 日韩欧美一区视频在线观看 | 亚洲国产精品999| 成人特级av手机在线观看| 午夜免费观看性视频| 国产亚洲欧美精品永久| 亚洲丝袜综合中文字幕| 多毛熟女@视频| 日韩精品免费视频一区二区三区 | 国产精品一区二区性色av| 免费观看的影片在线观看| 99视频精品全部免费 在线| 卡戴珊不雅视频在线播放| 岛国毛片在线播放| 嫩草影院入口| 人体艺术视频欧美日本| 一本大道久久a久久精品| 精品国产露脸久久av麻豆| 欧美日韩国产mv在线观看视频| 波野结衣二区三区在线| 一本久久精品| 大话2 男鬼变身卡| 乱系列少妇在线播放| 男的添女的下面高潮视频| 下体分泌物呈黄色| 久久久午夜欧美精品| 亚洲欧美日韩另类电影网站| 国产精品成人在线| 晚上一个人看的免费电影| 日本午夜av视频| 极品教师在线视频| 少妇人妻 视频| 亚洲美女搞黄在线观看| av免费在线看不卡| 国产一区有黄有色的免费视频| 一本—道久久a久久精品蜜桃钙片| 91成人精品电影| 18禁裸乳无遮挡动漫免费视频| av在线老鸭窝| 午夜久久久在线观看| 精品视频人人做人人爽| 免费av中文字幕在线| 另类亚洲欧美激情| 欧美日韩视频高清一区二区三区二| 国产精品99久久久久久久久| 校园人妻丝袜中文字幕| 黄色一级大片看看| av国产久精品久网站免费入址| 久久女婷五月综合色啪小说| 亚洲欧美清纯卡通| 欧美精品一区二区大全| 日韩一区二区三区影片| 免费观看a级毛片全部| 女的被弄到高潮叫床怎么办| 久久久久精品久久久久真实原创| 国产日韩欧美亚洲二区| 自拍欧美九色日韩亚洲蝌蚪91 | 嫩草影院新地址| 人妻制服诱惑在线中文字幕| 精品国产国语对白av| 国产熟女欧美一区二区| 曰老女人黄片| 久久99精品国语久久久| 亚洲国产成人一精品久久久| 久久久久久伊人网av| 国产色爽女视频免费观看| 亚洲av福利一区| 精品一品国产午夜福利视频| av一本久久久久| 人妻一区二区av| av国产精品久久久久影院| 女性生殖器流出的白浆| 国产精品久久久久成人av| 国产免费一级a男人的天堂| 简卡轻食公司| 免费人妻精品一区二区三区视频| 国产精品伦人一区二区| 亚洲av成人精品一二三区| 精品亚洲乱码少妇综合久久| 又粗又硬又长又爽又黄的视频| 日韩欧美一区视频在线观看 | av福利片在线观看| 老司机亚洲免费影院| 91久久精品电影网| 国产乱来视频区| 国产熟女午夜一区二区三区 | 色网站视频免费| 日本欧美国产在线视频| 91久久精品电影网| 天堂俺去俺来也www色官网| 亚洲精品日韩av片在线观看| 97在线人人人人妻| 亚洲内射少妇av| 亚洲精品一二三| 久久青草综合色| 亚洲av二区三区四区| 高清欧美精品videossex| 中文天堂在线官网| 免费看av在线观看网站| 精品一区二区免费观看| 精品久久久噜噜| 精品一区在线观看国产| 久久鲁丝午夜福利片| 亚洲婷婷狠狠爱综合网| 高清黄色对白视频在线免费看 | 亚洲av成人精品一二三区| 久久久国产欧美日韩av| av一本久久久久| 免费看光身美女| 精品国产国语对白av| 亚洲av日韩在线播放| 一级黄片播放器| 日本色播在线视频| 日本91视频免费播放| 一级a做视频免费观看| 国产成人一区二区在线| 国产精品久久久久久av不卡| 我的老师免费观看完整版| 99热全是精品| 少妇人妻 视频| 亚洲av成人精品一区久久| 丰满乱子伦码专区| 人妻制服诱惑在线中文字幕| 国产亚洲精品久久久com| 亚洲精品乱码久久久久久按摩| 国产精品伦人一区二区| 成人国产av品久久久| 日韩熟女老妇一区二区性免费视频| 一本久久精品| freevideosex欧美| 九九爱精品视频在线观看| 各种免费的搞黄视频| 久久毛片免费看一区二区三区| 91在线精品国自产拍蜜月| 亚洲人与动物交配视频| 国产真实伦视频高清在线观看| av线在线观看网站| 午夜久久久在线观看| 欧美激情国产日韩精品一区| 成人无遮挡网站| 国产女主播在线喷水免费视频网站| 韩国高清视频一区二区三区| 精品久久久久久电影网| 人妻 亚洲 视频| 亚洲情色 制服丝袜| 免费av不卡在线播放| 成年人午夜在线观看视频| av播播在线观看一区| 日日啪夜夜爽| 亚洲精品国产色婷婷电影| 少妇高潮的动态图| 日本欧美国产在线视频| 亚洲欧美精品专区久久| 丁香六月天网| 成人18禁高潮啪啪吃奶动态图 | 精品国产露脸久久av麻豆| 乱系列少妇在线播放| 精华霜和精华液先用哪个| av在线老鸭窝| 亚洲情色 制服丝袜| 国产男女超爽视频在线观看| 国产成人aa在线观看| 3wmmmm亚洲av在线观看| 亚洲三级黄色毛片| 最近中文字幕高清免费大全6| 国产有黄有色有爽视频| 三级国产精品片| 99久久精品热视频| 国产av一区二区精品久久| 啦啦啦视频在线资源免费观看| 男女免费视频国产| 成人二区视频| 久久久久视频综合| 伦精品一区二区三区| 人妻人人澡人人爽人人| 只有这里有精品99| 91成人精品电影| 久久影院123| 欧美日韩综合久久久久久| 97精品久久久久久久久久精品| 国产极品粉嫩免费观看在线 | 免费黄色在线免费观看| 丝瓜视频免费看黄片| 少妇人妻一区二区三区视频| 成年av动漫网址| 一级毛片 在线播放| 中文字幕久久专区| 激情五月婷婷亚洲| 国产成人一区二区在线| 一二三四中文在线观看免费高清| .国产精品久久| 一本一本综合久久| 日韩人妻高清精品专区| 热re99久久精品国产66热6| 久久久久精品久久久久真实原创| 中文在线观看免费www的网站| 国产亚洲91精品色在线| 免费看不卡的av| 久久97久久精品| 国产熟女欧美一区二区| 国内精品宾馆在线| 国产成人精品无人区| 最近2019中文字幕mv第一页| 久久久久久久久久久久大奶| 国产成人精品婷婷| 91精品一卡2卡3卡4卡| 日韩成人av中文字幕在线观看| 欧美最新免费一区二区三区| 成人无遮挡网站| 寂寞人妻少妇视频99o| 日韩亚洲欧美综合| 22中文网久久字幕| 性色av一级| 少妇丰满av| tube8黄色片| xxx大片免费视频| a级片在线免费高清观看视频| 久久99热这里只频精品6学生| 精品国产一区二区三区久久久樱花| 一本大道久久a久久精品| 精品人妻偷拍中文字幕| 水蜜桃什么品种好| 国产黄频视频在线观看| 丝袜在线中文字幕| 日本wwww免费看| 又爽又黄a免费视频| 少妇的逼好多水| 国产精品99久久久久久久久| 久久97久久精品| 日本-黄色视频高清免费观看| a 毛片基地| 欧美日韩视频精品一区| 久久精品国产亚洲av天美| 97超视频在线观看视频| 黄色欧美视频在线观看| av一本久久久久| 欧美激情国产日韩精品一区| 毛片一级片免费看久久久久| 人体艺术视频欧美日本| 男人爽女人下面视频在线观看| freevideosex欧美| 如何舔出高潮| 亚洲av.av天堂| 日韩一区二区视频免费看| 美女xxoo啪啪120秒动态图| 伦理电影大哥的女人| 美女cb高潮喷水在线观看| 久久久亚洲精品成人影院| 欧美3d第一页| 99久久精品热视频| 中文字幕久久专区| 丰满迷人的少妇在线观看| 女人久久www免费人成看片| 中文字幕精品免费在线观看视频 | 51国产日韩欧美| 国产极品粉嫩免费观看在线 | 亚洲图色成人| 又粗又硬又长又爽又黄的视频| 18+在线观看网站| 午夜免费鲁丝| 免费少妇av软件| 日韩av不卡免费在线播放| 久久99热6这里只有精品| 亚洲国产av新网站| av视频免费观看在线观看| 国产黄片视频在线免费观看| 国产伦理片在线播放av一区| 六月丁香七月| 22中文网久久字幕| 久久精品久久久久久久性| 免费观看a级毛片全部| 九九爱精品视频在线观看| 99久久中文字幕三级久久日本| 涩涩av久久男人的天堂| 亚洲av福利一区|