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

    乘用車動力傳動系瞬態(tài)工況扭振分析與改進(jìn)

    2017-07-06 11:02:58楊志堅喻桂華
    關(guān)鍵詞:減振器飛輪瞬態(tài)

    袁 旺,楊志堅,喻桂華

    (華南理工大學(xué) 機(jī)械與汽車工程學(xué)院, 廣州 510641)

    ?

    乘用車動力傳動系瞬態(tài)工況扭振分析與改進(jìn)

    袁 旺,楊志堅,喻桂華

    (華南理工大學(xué) 機(jī)械與汽車工程學(xué)院, 廣州 510641)

    考慮離合器扭轉(zhuǎn)減振器非線性剛度、齒輪嚙合時變剛度和齒輪側(cè)隙,建立了乘用車傳動系3擋集中參數(shù)扭振模型。用用LMS virtual.lab和Simulink聯(lián)合仿真模擬乘用車反復(fù)點踩油門-松油門的瞬態(tài)工況,飛輪、輸入軸瞬態(tài)轉(zhuǎn)速變化的仿真與實驗曲線對比驗證了模型的有效性。仿真分析發(fā)現(xiàn):小剛度的離合器扭轉(zhuǎn)減振器可以有效改善傳動系的瞬態(tài)性能。模型中換用雙質(zhì)量飛輪,仿真結(jié)果表明:在反復(fù)點踩油門-松油門的瞬態(tài)工況下,飛輪和輸入軸轉(zhuǎn)速波動幅值分別減小了14和38 r/min,輸入軸瞬時轉(zhuǎn)速中的沖擊變化量幾乎完全消除。

    動力傳動系;扭振;瞬態(tài)工況;雙質(zhì)量飛輪

    在瞬態(tài)工況,如啟停、急加速和點踩油門等時,發(fā)動機(jī)激勵會發(fā)生突變,這種突變?nèi)菀准顒恿鲃酉狄缘碗A扭轉(zhuǎn)固有頻率振動。同時,由于傳動系各部件之間存在間隙,如齒輪側(cè)隙、花鍵間隙等,在瞬態(tài)工況下容易引起敲擊、異響等振動噪聲問題。另外,離合器分段非線性剛度、齒輪副嚙合時變剛度等非線性因素也是造成瞬態(tài)行駛工況下汽車行駛不平順、乘坐舒適性惡化的原因。因此,研究并改進(jìn)汽車瞬態(tài)行駛工況下動力傳動系扭振性能對提升整車NVH性能具有重要意義。

    相比穩(wěn)定工況,有關(guān)傳動系瞬態(tài)工況扭振性能的研究文獻(xiàn)相對較少。韓海蘭[1]利用ADAMS/Driveline建立某微型客車動力傳動系模型,并與其他子系統(tǒng)通過合理的通訊器連接,組裝成微型客車整車虛擬樣機(jī)模型,實現(xiàn)了水平路面起步、加速、減速3種非穩(wěn)定工況的仿真。王東[2]在使用AMEsim建立整車動力傳動系扭振模型的基礎(chǔ)上,對緩加速工況的扭振性能進(jìn)行仿真,其結(jié)果與實驗測試一致性良好。閆明剛等[3]以某型混合動力公交車傳動系為對象,考慮傳動系的非線性特性,應(yīng)用AMESim建立系統(tǒng)扭轉(zhuǎn)振動仿真模型,計算了在包含啟動、低速和高速運行以及停機(jī)過程的完整工作循環(huán)下系統(tǒng)的扭振響應(yīng)。Raja Mazuir、Raja Ahsan Shah等[4]建模研究了某四驅(qū)車2擋快速點踩油門、松油門(tip-in/out)時的瞬態(tài)動力學(xué)性能。Fu Sheng-ping等[5]基于集中質(zhì)量法建立了車輛動力傳動系扭轉(zhuǎn)振動模型,設(shè)計了基于遺傳算法動態(tài)優(yōu)化理論的模型優(yōu)化程序,優(yōu)化結(jié)果表明:該方法有效衰減了系統(tǒng)的扭轉(zhuǎn)振動。A Farshidianfar等[6]將傳動系等效為部分集中、部分分布的混合模型,在驗證模型有效性的基礎(chǔ)上研究了瞬態(tài)工況下經(jīng)常會出現(xiàn)的叮當(dāng)聲(clonk)和拖曳(shuffle)現(xiàn)象。M Gnanakumarr等[7]利用剛?cè)狁詈夏P脱芯苛嗽诎l(fā)動機(jī)扭矩突然增大時由于齒輪間隙存在而引起的齒輪敲擊現(xiàn)象,結(jié)果表明:傳動系扭振中對應(yīng)驅(qū)動軸結(jié)構(gòu)模態(tài)的高頻響應(yīng)是由變速箱齒隙造成的沖擊激發(fā)的。

    以上文獻(xiàn)針對瞬態(tài)工況傳動系整體性能進(jìn)行研究,但較少關(guān)注系統(tǒng)具體參數(shù)對傳動系瞬態(tài)性能的影響。本文基于集中參數(shù)扭振模型的建模理論,考慮傳動系某些部件的時變特性和非線性特性,建立了某乘用車3擋行駛時的動力傳動系的扭振模型,利用LMS virtual.lab和Simulink聯(lián)合仿真研究了反復(fù)點踩油門-松油門(tip-in/out)工況下系統(tǒng)的瞬態(tài)動力學(xué)性能,在對比實驗測試結(jié)果驗證模型有效性的基礎(chǔ)上著重討論了離合器剛度特性對傳動系瞬態(tài)性能的影響,并設(shè)法尋找改進(jìn)傳動系瞬態(tài)性能的方法。

    1 傳動系模型建立及系統(tǒng)模態(tài)分析

    某乘用車動力傳動系結(jié)構(gòu)如圖1所示,發(fā)動機(jī)為四缸四沖程,前輪驅(qū)動,采用一款5擋手動變速箱。

    圖1 某車型動力傳動系結(jié)構(gòu)示意圖

    1.1 傳動系集中參數(shù)模型的建立

    依據(jù)集中參數(shù)扭振模型建立的理論,根據(jù)動力傳動系各軸(曲軸、變速器輸入軸、輸出軸、半軸等)轉(zhuǎn)動慣量和傳遞力矩的分布情況,分別等效為由離散的轉(zhuǎn)動慣量、無質(zhì)量的扭轉(zhuǎn)彈簧和阻尼器組成的力學(xué)模型[8]。3擋扭振模型如圖2所示,模型各部分對應(yīng)的物理模型及其參數(shù)如表1所示,模型中各部分阻尼的選取主要參考文獻(xiàn)[9]。在使用LMS virtual.lab的Motion模塊建立仿真模型時,為了更準(zhǔn)確地體現(xiàn)系統(tǒng)非穩(wěn)定工況下的特性,考慮了傳動系某些部件的時變特性和非線性特性:

    1) 離合器扭轉(zhuǎn)減振器的非線性剛度

    在行駛過程中,離合器扭轉(zhuǎn)減振器的剛度在扭轉(zhuǎn)角達(dá)到一定值時會發(fā)生變化,尤其在非穩(wěn)定工況下(如起步、急加速和急減速等),由于離合器扭轉(zhuǎn)角變化范圍較大,離合器剛度會發(fā)生突變,這種突變極易引起傳動系的沖擊和異響。根據(jù)廠家提供的數(shù)據(jù),所用離合器剛度特性為

    (1)

    式中:θ為離合器扭轉(zhuǎn)減震器主從動部分相對扭轉(zhuǎn)角(°);Kθ為離合器扭轉(zhuǎn)剛度(N·m/rad)。

    2) 承載齒輪副的嚙合時變剛度以及齒輪側(cè)隙

    在圖2所示的模型中,考慮了變速箱與主減速器中的承載齒輪副,I3、I4之間為3擋齒輪副,I5、I6之間為主減速器齒輪副。嚙合時變剛度是齒輪嚙合過程中參與嚙合的齒數(shù)變化引起的,在模型中主要體現(xiàn)在齒輪重合度這一參數(shù)上,取值為2.5。在齒輪嚙合傳動時,為了在齒廓之間形成潤滑油膜,避免因輪齒摩檫發(fā)熱膨脹而卡死,齒廓之間必須留有間隙,此間隙稱為側(cè)隙。在Motion中,將齒輪作為力元素來模擬,可以比較方便地考慮齒隙。側(cè)隙引起的剛度變化[10]為

    (2)

    式中:KC是嚙合剛度;xr為嚙合齒輪節(jié)圓處周向相對位移;km為單對齒輪嚙合時線性剛度值;a為節(jié)圓處側(cè)隙。3擋齒輪副中心距為80 mm。根據(jù)文獻(xiàn)[11]提供的經(jīng)驗值,側(cè)隙取0.1 mm。主減速器齒輪副中心距為127 mm,側(cè)隙取0.17 mm。

    圖2 模型示意圖

    轉(zhuǎn)動慣量/(kg·mm2)扭轉(zhuǎn)剛度/(N·m·rad-1)曲軸飛輪總成、活塞連桿機(jī)構(gòu)、離合器主動件及離合器從動盤的大部分I1159279.4離合器扭轉(zhuǎn)減震器K1,2兩段剛度離合器從動盤中與輸入軸相連的部分及輸入軸前端部分I23337.2三擋時輸入軸扭轉(zhuǎn)剛度K2,311044變速箱輸入軸及其帶動齒輪I31321.12三擋齒輪副嚙合扭轉(zhuǎn)剛度K3,4時變剛度(齒輪重合度2.5)變速箱輸出軸與輸入軸嚙合部分I42576.63三擋時輸出軸扭轉(zhuǎn)剛度K4,5118530變速箱輸出軸與差減嚙合部分I5902.7主減齒輪副嚙合扭轉(zhuǎn)剛度K5,6時變剛度(齒輪重合度2.5)差減總成及半軸之半I649654半軸扭轉(zhuǎn)剛度K6,712662.38輪胎及半軸之半I73000000輪胎扭轉(zhuǎn)剛度K7,817437.74整車平動質(zhì)量等效轉(zhuǎn)動慣量I8120000000

    1.2 傳動系扭轉(zhuǎn)模態(tài)分析

    為了便于計算系統(tǒng)固有模態(tài),對模型進(jìn)行進(jìn)一步簡化。由于嚙合齒輪副剛度很大,將被動部分通過速比轉(zhuǎn)化到主動部分,將二者視為一體,即:

    (3)

    (4)

    其中:n1為3擋齒輪副傳動比;n2為主減速器齒輪副傳動比。以曲軸轉(zhuǎn)速為基準(zhǔn)將系統(tǒng)各部分轉(zhuǎn)動慣量及扭轉(zhuǎn)剛度按照動能和彈性變形能相等的原則進(jìn)行轉(zhuǎn)化[12],同時忽略系統(tǒng)的阻尼,得到更簡單的6自由度集中參數(shù)模型。利用Matlab編程計算各階固有頻率,具體頻率值如表2所示。

    表2 系統(tǒng)前5階固有頻率

    2 傳動系瞬態(tài)性能分析

    2.1 tip-in/out工況的仿真模擬

    tip-in/out工況是指在汽車某一擋位下,踩油門加速至一定車速,然后快速松油門,減速至一定車速后再快速點踩油門,如此反復(fù)。對于擁擠道路,由于頻繁變道需要經(jīng)常性地急加速和急減速,tip-in/out工況便是反映這種情況。

    實驗在跑道上進(jìn)行,在飛輪齒盤、變速箱輸入軸固聯(lián)齒(2擋)和輸出軸固聯(lián)齒(3擋)處動力總成外殼打孔,利用霍爾傳感器測試飛輪及輸入、輸出軸處的轉(zhuǎn)速脈沖信號。在3擋40 km/h車速附近反復(fù)點踩油門-松油門(tip-in/out),持續(xù)時間約25 s,各測點轉(zhuǎn)速變化曲線如圖3所示。飛輪轉(zhuǎn)速大致在1 780~2 060 r/min內(nèi)變化。在仿真時,使用LMS virtual.lab和Simulink進(jìn)行聯(lián)合仿真,Simulink 的控制模型如圖4所示。

    圖3 實驗轉(zhuǎn)速變化

    圖4 Simulink控制模型

    以發(fā)動機(jī)轉(zhuǎn)速,即圖2中I1的轉(zhuǎn)速n作為Simulink的輸入變量,發(fā)動機(jī)輸出扭矩Ts作為Simulink的輸出變量施加在I1上,輸出扭矩由發(fā)動機(jī)扭矩外特性曲線乘以一個隨時間變化的系數(shù)k實現(xiàn)。該車型采用一款1.3T直列四缸四沖程發(fā)動機(jī),全負(fù)荷下平均輸出轉(zhuǎn)矩Te與轉(zhuǎn)速n的關(guān)系為

    Te=-4.56e-12n4+6.172e-08n3-

    3.042e-4n2+0.638 9n-254.2

    (5)

    式中:Te單位為N·m;n單位為r/min。系數(shù)k的選取使得仿真飛輪轉(zhuǎn)速隨時間變化曲線與實驗轉(zhuǎn)速曲線一致。系數(shù)k隨時間變化的曲線如圖5所示。

    圖5 仿發(fā)動機(jī)輸出扭矩系數(shù)k隨時間變化的曲線

    與實驗相對應(yīng),該系數(shù)的變化周期為7 s,在1個周期(0~7 s)內(nèi)系數(shù)k與時間的關(guān)系為:

    (6)

    因此,仿真時Simulink的輸出扭矩為

    Ts=Te·k

    (7)

    由式(5)~(7)可以確定仿真時Simulink輸入變量n與輸出變量Ts之間的關(guān)系。利用經(jīng)驗公式(8)(9)計算得到輪胎滾動阻力矩Tf、行駛空氣阻力等效阻力矩TW,并分別施加在圖2中車輪I7、整車平動質(zhì)量I8處[13]。

    (8)

    (9)

    式中:u為車速(km/h);A為汽車行駛時迎風(fēng)面的面積;C為風(fēng)阻系數(shù),CA一般在0.6~0.9之間,本研究取0.75;r是輪胎滾動半徑,為0.317 m。

    仿真時間設(shè)置為25 s,得到圖2中飛輪I1、輸入軸I3和輸出軸I5的轉(zhuǎn)速信號。圖6、7和8分別為飛輪、輸入軸和輸出軸處轉(zhuǎn)速曲線的實驗、仿真結(jié)果對比。從圖中可以看出:實驗和仿真各測點轉(zhuǎn)速變化基本一致。由于k的值是通過仿真調(diào)試得到的,與實驗實際情況存在一定誤差,同時仿真模型各參數(shù)的值也存在一定的誤差,所以仿真轉(zhuǎn)速和實驗轉(zhuǎn)速存在一定的差別,但差距較小。仿真與實驗得到的飛輪最大轉(zhuǎn)速差為42.3 r/min,最大相對誤差為2.35%。

    圖6 實驗和仿真得到的飛輪轉(zhuǎn)速

    圖7 實驗和仿真得到的輸入軸轉(zhuǎn)速

    圖8 實驗和仿真得到的輸出軸轉(zhuǎn)速

    2.2 tip-in/out傳動系轉(zhuǎn)速特性分析

    在實現(xiàn)3擋tip-in/out工況模擬的基礎(chǔ)上,著重分析此工況下動力傳動系的轉(zhuǎn)速變化特性。圖9是松油門(tip-out)時飛輪和輸入軸轉(zhuǎn)速曲線放大圖,從圖中可以看出:在快速放松油門后的0.5 s內(nèi),二者轉(zhuǎn)速均出現(xiàn)較大的不規(guī)則波動,最大波動值達(dá)到88 r/min,相比飛輪,輸入軸轉(zhuǎn)速存在更高頻的不規(guī)則波動;0.5 s之后進(jìn)入正常減速階段,此時飛輪輸入軸轉(zhuǎn)速曲線基本重合。從圖9中看出:7.894~8.383 s的0.489 s內(nèi),轉(zhuǎn)速波動3個周期,可以計算出波動頻率約為6.135 Hz,與本文1.2節(jié)計算得到傳動系1階固有頻率5.710 Hz接近。同理,圖10是tip-in時飛輪和輸入軸轉(zhuǎn)速曲線放大圖,從圖中可以看出:快速點踩油門后,二者轉(zhuǎn)速迅速從1 784 r/min增大到1 867 r/min,又減小到1 791 r/min,且輸入軸轉(zhuǎn)速出現(xiàn)明顯的高頻波動;之后進(jìn)入正常加速階段,飛輪與輸入軸轉(zhuǎn)速基本重合,轉(zhuǎn)速波動頻率約6.135 Hz,仍然與計算得到傳動系1階固有頻率5.710 Hz接近。

    在實驗測試3擋點踩油門(tip-in)和松油門(tip-out)工況時,車內(nèi)乘員主觀上明顯感受到車身前后竄動,同時聽到變速箱有明顯的“嗒”音異響。

    圖9 實驗tip-out轉(zhuǎn)速放大圖

    圖10 實驗tip-in轉(zhuǎn)速放大圖

    圖11是仿真tip-out飛輪和輸入軸轉(zhuǎn)速曲線放大圖,從圖中可以看出:在快速放松油門后的0.5 s內(nèi),轉(zhuǎn)速出現(xiàn)較大的不規(guī)則波動,最大波動量達(dá)到72 r/min,相比飛輪,輸入軸轉(zhuǎn)速存在更高頻率的不規(guī)則波動;0.5 s之后進(jìn)入正常減速階段,此時轉(zhuǎn)速曲線波動頻率約為5.61 Hz,與本文1.2節(jié)計算得到傳動系1階固有頻率5.710 Hz接近。圖12是仿真tip-in時飛輪和輸入軸轉(zhuǎn)速曲線放大圖,從圖中可以看出:仿真輸入力矩突然增大后,飛輪及輸入軸轉(zhuǎn)速迅速從1 778 r/min增大到1 812 r/min,又減小到1 764 r/min,且輸入軸轉(zhuǎn)速出現(xiàn)明顯的高頻波動;之后進(jìn)入正常加速階段,轉(zhuǎn)速波動頻率約5.59 Hz,仍然與計算得到傳動系1階固有頻率5.710 Hz接近。

    以上分析表明,在tip-in/out工況下,仿真飛輪、輸入軸轉(zhuǎn)速變化的瞬態(tài)特性與實驗測試結(jié)果比較接近,充分證明了仿真模型和仿真方法的有效性。

    圖11 仿真tip-out轉(zhuǎn)速放大圖

    圖12 仿真tip-in轉(zhuǎn)速放大圖

    在tip-in/out工況下,發(fā)動機(jī)激勵的頻率成分主要集中在低頻,容易激起傳動系第1階固有振動(頻率在10 Hz以下),而在實際行駛過程中,由式(5)(8)和(9)可知:汽車發(fā)動機(jī)激勵、滾動阻力和行駛阻力均為車速(發(fā)動機(jī)轉(zhuǎn)速)的函數(shù),因此第1階固有振動的轉(zhuǎn)速波動必然會引起發(fā)動機(jī)激勵、輪胎滾動阻力和空氣阻力的同頻率波動。在特定的輪胎滑移率范圍內(nèi),力的波動又會反作用于轉(zhuǎn)速,加劇系統(tǒng)以該頻率的轉(zhuǎn)速波動,從原理上講是一種具有自激振動性質(zhì)的扭轉(zhuǎn)振動[14-15]。由于傳動系中大阻尼部件,如離合器扭轉(zhuǎn)減振器、輪胎等的存在,這種自激振動可以得到有效的衰減,而不至于造成系統(tǒng)振動幅值的發(fā)散。在主觀上,這種具有自激性質(zhì)的扭轉(zhuǎn)振動表現(xiàn)為車身的前后竄動。在上述tip-in/out工況下,乘員所感受到的車身前后竄動現(xiàn)象正是由此引起的。

    2.3 離合器剛度特性對瞬態(tài)性能的影響

    由于實驗測試時存在測點處齒數(shù)的限制(飛輪齒數(shù)為112、輸入軸齒數(shù)為22、輸出軸齒數(shù)為41),霍爾傳感器測得的轉(zhuǎn)速信號采樣頻率較低,圖9、圖10中的實驗轉(zhuǎn)速很難清楚地反映出點踩油門-松油門瞬間0.5 s內(nèi)轉(zhuǎn)速的變化情況。相比之下,仿真可以達(dá)到很高的采樣精度。圖11、12中轉(zhuǎn)速的采樣頻率為5 000 Hz,可以清楚地分辨出0.5 s甚至更短時間內(nèi)瞬時轉(zhuǎn)速的變化情況。因此,本文通過對仿真結(jié)果的進(jìn)一步分析找到tip-in/out時輸入軸轉(zhuǎn)速發(fā)生高頻波動和沖擊變化的原因。圖13是仿真tip-out時飛輪、輸入軸轉(zhuǎn)速在6.8~7.3 s時間段內(nèi)的曲線,圖14是仿真離合器扭轉(zhuǎn)減振器主從動部分相對扭轉(zhuǎn)角變化曲線,圖15為該乘用車離合器的相對扭角-扭矩特性曲線。從圖13中可以看出:t=6.934 s時飛輪與輸入軸轉(zhuǎn)速曲線開始出現(xiàn)明顯的分離。圖14顯示此時離合器相對扭轉(zhuǎn)角為7.49°。由圖15可知:相對扭角為7.49°時,離合器扭轉(zhuǎn)減振器的扭轉(zhuǎn)剛度發(fā)生突變,變化量的絕對值為2 043.74 N·m/rad。同理,在t分別為7.035、7.234 s時,輸入軸轉(zhuǎn)速均出現(xiàn)拐點和沖擊變化,此時離合器相對扭轉(zhuǎn)角分別為-2.975°、-2.998°,與-2.99°接近。由圖15可知,相對扭轉(zhuǎn)角為-2.99°時,離合器扭轉(zhuǎn)剛度同樣發(fā)生了突變,變化量絕對值為2 043.74 N·m/rad。由以上分析可以初步判斷離合器扭轉(zhuǎn)剛度的突然變化是造成輸入軸轉(zhuǎn)速突變的主要原因。

    圖13 在 tip-out下仿真得到的轉(zhuǎn)速放大圖

    圖14 在 tip-out工況下仿真得到的離合器相對扭轉(zhuǎn)角

    圖15 離合器扭轉(zhuǎn)角-扭矩曲線

    為了進(jìn)一步明確離合器扭轉(zhuǎn)減振器剛度特性對tip-in/out工況下傳動系性能的影響,為該乘用車匹配另一款小剛度大轉(zhuǎn)角離合器。新離合器扭轉(zhuǎn)減振器的剛度特性如下:

    (10)

    式中:θ為離合器扭轉(zhuǎn)減震器主從動部分相對扭轉(zhuǎn)角(°);Kθ為減振器扭轉(zhuǎn)剛度(N·m/rad)。相比原離合器,該離合器具有扭轉(zhuǎn)剛度小、轉(zhuǎn)角大的特點。仿真研究此時傳動系瞬態(tài)性能,仿真的其余參數(shù)不變。圖16、17為tip-out工況下瞬時大小剛度離合器仿真結(jié)果對比。圖18為小剛度離合器tip-out工況下瞬時轉(zhuǎn)速仿真曲線。圖19為仿真得到的小剛度離合器相對扭轉(zhuǎn)角。

    從圖16可以看出:對于不同剛度特性的離合器,在tip-out工況下瞬時飛輪轉(zhuǎn)速波動的幅值變化很小。觀察圖17可知:換用小剛度離合器后,輸入軸轉(zhuǎn)速波動的幅值有所減小。

    圖16 仿真得到的大小剛度離合器飛輪轉(zhuǎn)速對比

    圖17 在tip-out工況下仿真得到的大小剛度離合器輸入軸轉(zhuǎn)速對比

    由圖18、19可以看出:在t=6.942 s時,離合器扭轉(zhuǎn)角為8.499°,接近于8.5°,此時離合器扭轉(zhuǎn)減振器剛度發(fā)生突變,變化絕對值為504.2 N·m/rad,飛輪轉(zhuǎn)速與輸入軸轉(zhuǎn)速的差別開始變大。由圖19看出:在t分別為7.057、7.18 s時,離合器扭轉(zhuǎn)角分別為-2.997°、-3.002°,與-3°接近,離合器扭轉(zhuǎn)剛度同樣存在絕對值為504.2 N·m/rad的突變。而圖18表明這2個時刻輸入軸轉(zhuǎn)速曲線均出現(xiàn)拐點和沖擊變化。由圖17可以看出:與原先離合器相比,小剛度大轉(zhuǎn)角離合器剛度突變值由之前的2 043.74 N·m/rad變?yōu)?04.2 N·m/rad,輸入軸轉(zhuǎn)速曲線中的沖擊變化量也明顯減小。由此進(jìn)一步確認(rèn)離合器扭轉(zhuǎn)剛度的突然變化是造成輸入軸轉(zhuǎn)速突變的主要原因,而輸入軸轉(zhuǎn)速的突變必然會引起變速箱齒輪轉(zhuǎn)速的突變。由于齒隙的存在,齒輪轉(zhuǎn)速的突變造成了嚙合齒輪間的敲擊,進(jìn)而產(chǎn)生了上述實驗測試時聽到的變速箱“嗒”音異響。

    通過本節(jié)的分析可知:離合器扭轉(zhuǎn)減振器扭轉(zhuǎn)剛度對動力傳動系瞬態(tài)性能有較大的影響,小剛度離合器可以有效地降低在tip-in/out工況下輸入軸轉(zhuǎn)速的波動幅值;對于多段剛度特性的離合器,相鄰兩段剛度差的絕對值越小,在tip-in/out工況下瞬時輸入軸轉(zhuǎn)速的沖擊變化量也就越小,有利于改善此工況下變速箱敲擊引起的“嗒”音等問題。因此,適當(dāng)減小離合器扭轉(zhuǎn)減振器剛度,降低相鄰剛度段之間扭轉(zhuǎn)剛度的差值對于改善傳動系瞬態(tài)性能是很有效的。

    圖18 在tip-ou工況下仿真得到的小剛度離合器飛輪、輸入軸轉(zhuǎn)速

    圖19 仿真得到的小剛度離合器相對扭轉(zhuǎn)角

    3 傳動系瞬態(tài)性能改進(jìn)

    3.1 改進(jìn)方法

    由以上分析可知,傳動系瞬態(tài)性能受離合器剛度特性的影響很大,小剛度大轉(zhuǎn)角離合器可以顯著改善動力傳動系的瞬態(tài)性能。但是,離合器從動盤式扭轉(zhuǎn)減振器由于自身結(jié)構(gòu)的限制,相對轉(zhuǎn)角不可能很大,為了有效傳遞發(fā)動機(jī)力矩滿足動力性要求,扭轉(zhuǎn)減振器剛度一般較大。而雙質(zhì)量飛輪(DMF)扭轉(zhuǎn)減振器的彈性元件可以布置在較大的空間內(nèi),相對扭轉(zhuǎn)角大,扭轉(zhuǎn)剛度可以設(shè)計得很小,且便于設(shè)計更多級剛度改善瞬態(tài)減振性能[16]。因此,考慮使用雙質(zhì)量飛輪改善傳動系瞬態(tài)特性。

    在仿真模型中換用雙質(zhì)量飛輪,其初級慣量與圖2模型中轉(zhuǎn)動慣量I1一樣,次級慣量對應(yīng)圖2模型的轉(zhuǎn)動慣量I2,具體數(shù)值如表3所示。而雙質(zhì)量飛輪扭轉(zhuǎn)減振器扭轉(zhuǎn)剛度對應(yīng)圖2模型中K1,2,扭轉(zhuǎn)阻尼對應(yīng)C1,2,模型其余參數(shù)值不變。

    表3 DMF模型對應(yīng)轉(zhuǎn)動慣量值

    根據(jù)廠家提供的數(shù)據(jù),該雙質(zhì)量飛輪扭轉(zhuǎn)減振器剛度特性如下:

    (11)

    阻尼C1,2仍然使用前面的值,用相同的方法計算DMF工作在最大剛度,即1 048.513 N·m/rad時系統(tǒng)各階固有頻率,頻率值在表4中列出。可以看出,換用DMF后,第2、3、4階固有頻率顯著減小。

    表4 換用DMF系統(tǒng)前五階固有頻率

    3.2 改進(jìn)效果驗證

    換用DMF后,仿真與之前相同的工況,結(jié)果如圖20、21所示。從圖20可以看出:換用DMF后,飛輪在tip-out工況下轉(zhuǎn)速波動幅值最大值由61 r/min減小到47 r/min。從圖21可以看出:換用DMF后,輸入軸zai tip-out工況下轉(zhuǎn)速波動幅值明顯減小,最大幅值由58 r/min減小到20 r/min。由于輸入軸轉(zhuǎn)速波動的減小,主觀感受車身前后竄動的幅度也會相應(yīng)的減小,有利于提升tip-in/out工況下的乘坐舒適性。進(jìn)一步分析圖21發(fā)現(xiàn):輸入軸轉(zhuǎn)速曲線中的拐點和沖擊變化量完全被消除,極大改善了tip-in/out工況下動力傳動系瞬態(tài)性能。而在主觀評價方面,換用DMF后,在tip-in/out工況下變速箱“嗒”音異響基本消除。

    圖20 在tip-out工況下仿真得到的飛輪轉(zhuǎn)速放大圖

    圖21 在tip-out工況下仿真得到的輸入軸轉(zhuǎn)速放大圖

    4 結(jié)論

    1) 考慮離合器扭轉(zhuǎn)減振器非線性剛度、齒輪嚙合時變剛度和齒輪側(cè)隙,建立了乘用車動力傳動系3擋集中參數(shù)扭振模型。

    2) LMS virtual.lab和Simulin聯(lián)合仿真模擬乘用車反復(fù)點踩油門-松油門的瞬態(tài)工況,通過仿真與實驗得到的飛輪、輸入軸瞬態(tài)轉(zhuǎn)速變化曲線對比驗證了模型的有效性。仿真分析發(fā)現(xiàn):小剛度的離合器扭轉(zhuǎn)減振器可以有效改善傳動系的瞬態(tài)性能。

    3) 在模型中換用雙質(zhì)量飛輪。仿真結(jié)果表明:在反復(fù)點踩油門-松油門的瞬態(tài)工況下,飛輪在tip-out工況下轉(zhuǎn)速波動幅值最大值由61 r/min減小到47 r/min,輸入軸在tip-out工況下轉(zhuǎn)速波動幅值明顯減小,最大幅值由58 r/min減小到20 r/min,轉(zhuǎn)速曲線中的拐點和沖擊變化量完全被消除。

    [1] 韓海蘭.微型客車動力傳動系統(tǒng)扭振及其對整車振動影響的研究[D].長春:吉林大學(xué),2012.

    [2] 王東.動力傳動系扭振激勵的整車振動研究[D].成都:西南交通大學(xué),2015.

    [3] 閆明剛,張勇,侯之超,等.混合動力客車傳動系扭振響應(yīng)及其影響因素分析[C]//全國結(jié)構(gòu)工程學(xué)術(shù)會議.2013:223-227.

    [4] SHAH B R A,MAZUIR R.Modeling of 4WD vehicle driveability during tip-in/tip-out events[C]//Proceedings of Iavsd Symposium.2011.

    [5] FU S,LI S,LUO N,et al.Dynamic optimization of tracked vehicle power train based on torsional vibration analysis[J].Advances in Mechanical Engineering,2016,8(5):115-119.

    [6] FARSHIDIANFAR A,EBRAHIMI M,BARTLETT H.Hybrid modeling and simulation of the torsional vibration of vehicle driveline systems[J].Proceedings of the Institution of Mechanical Engineers Part D Journal of Automobile Engineering,2001,215(2):217-229.

    [7] GNANAKUMARR M,THEODOSSIADES S,RAHNEJAT H,et al.Impact-induced vibration in vehicular driveline systems:Theoretical and experimental investigations[J].Proceedings of the Institution of Mechanical Engineers-Part K,2005,219(1):1-12.

    [8] 肖啟瑞,陳黎明,黃晉.重型牽引車傳動系動力學(xué)特性[J].重慶理工大學(xué)學(xué)報(自然科學(xué)),2013,27(11):23-27.

    [9] 湯志華.汽車動力傳動系扭振分析軟件開發(fā)及仿真研究[D].成都:西南交通大學(xué),2006.

    [10]CROWTHER A R,ZHANG N.Torsional finite elements and nonlinear numerical modeling in vehicle powertrain dynamics[J].Journal of Sound & Vibration,2005,284(3/4/5):825-849.

    [11]李玉龍.外嚙合齒輪泵困油機(jī)理、模型及試驗研究[D].合肥:合肥工業(yè)大學(xué),2009.

    [12]何煦.某型轎車傳動系統(tǒng)扭轉(zhuǎn)振動的研究[D].上海:上海交通大學(xué),2008.

    [13]毛恒.雙質(zhì)量飛輪系統(tǒng)扭振試驗方法研究[D].武漢:武漢理工大學(xué),2010.

    [14]鄭聯(lián)珠,劉明樹,張友坤,等.汽車傳動系自激扭振機(jī)理研究[J].汽車工程,2001,23(6):407-410.

    [15]趙克剛,李健,姚偉浩.車輛傳動系統(tǒng)固有特性的識別[J].重慶理工大學(xué)學(xué)報(自然科學(xué)),2014,28(4):5-9.

    [16]李偉,史文庫.雙質(zhì)量飛輪(DMF)的研究綜述[J].噪聲與振動控制,2008,28(5):1-5.

    (責(zé)任編輯 劉 舸)

    Study on Torsional Vibration Characteristics of Passenger Vehicle Power Transmission System During Transient Process

    YUAN Wang,YANG Zhi-jian,YU Gui-hua

    (School of Mechanical and Automotive Engineering, South China University of Technology, Guangzhou 510641, China)

    A lumped parameter model of passenger vehicle in third gear was established in which some factors such as nonlinear stiffness of clutch, the backlash and time-varying meshing stiffness of gear pairs were considered. The model was verified by comparing rotate speed signals of experiments and united simulation by LMS virtual.lab and Matlab/Simulink under tip-in/out operating condition in third gear. Simulation research results indicated that the clutch torsional vibration damper with the lower stiffness can improve transient performance of transmission system obviously. Equipping model with dual mass flywheel and the simulation results demonstrated the tip-in/out transient rotate speed fluctuation amplitude of flywheel and input shaft decreased by 14 and 38 r/min respectively and the impact of input shaft transient speed curve were eliminated almost completely.

    power transmission system; torsional vibration; transient process; dual mass flywheel

    2017-02-16

    袁旺(1992—),男,山西汾陽人,碩士研究生,主要研究從事動力傳動系扭振方面的研究,E-mail:yw0210301@163.com;通訊作者 楊志堅(1982—),男,湖南新化人,博士,副教授,主要從事機(jī)械故障診斷與信號處理、振動與控制方面的研究,E-mail:yangzhj@scut.edu.cn。

    袁旺,楊志堅,喻桂華.s[J].重慶理工大學(xué)學(xué)報(自然科學(xué)),2017(6):31-40.

    format:YUAN Wang,YANG Zhi-jian,YU Gui-hua.Study on Torsional Vibration Characteristics of Passenger Vehicle Power Transmission System During Transient Process[J].Journal of Chongqing University of Technology(Natural Science),2017(6):31-40.

    10.3969/j.issn.1674-8425(z).2017.06.005

    TB53;U46

    A

    1674-8425(2017)06-0031-10

    猜你喜歡
    減振器飛輪瞬態(tài)
    高壓感應(yīng)電動機(jī)斷電重啟時的瞬態(tài)仿真
    飛輪座注射模設(shè)計
    模具制造(2019年7期)2019-09-25 07:30:00
    輪峰推出兩款飛輪新產(chǎn)品
    十億像素瞬態(tài)成像系統(tǒng)實時圖像拼接
    汽車減振器與磁流變材料分析與運算
    基于瞬態(tài)流場計算的滑動軸承靜平衡位置求解
    Word Formation in English for Science and Technology
    DC/DC變換器中的瞬態(tài)特性分析
    推輪子的人
    意林(2014年1期)2014-07-05 05:54:04
    饋能式磁流變減振器自供電特性研究
    汽車科技(2014年6期)2014-03-11 17:45:36
    日本撒尿小便嘘嘘汇集6| 免费av毛片视频| 日韩欧美一区二区三区在线观看| 色哟哟·www| 精品人妻1区二区| 国产乱人视频| 国产精品,欧美在线| 精品午夜福利在线看| 90打野战视频偷拍视频| 精品免费久久久久久久清纯| 日韩高清综合在线| 老司机午夜十八禁免费视频| 尤物成人国产欧美一区二区三区| 国产亚洲精品av在线| 国产乱人伦免费视频| 我的老师免费观看完整版| 久久久久久九九精品二区国产| 自拍偷自拍亚洲精品老妇| 国产欧美日韩一区二区三| 国产三级黄色录像| 欧美成人a在线观看| 18禁黄网站禁片午夜丰满| 日本在线视频免费播放| 国产精品亚洲美女久久久| 在线观看免费视频日本深夜| 神马国产精品三级电影在线观看| 美女免费视频网站| 亚洲av成人av| 99riav亚洲国产免费| 国产日本99.免费观看| 天堂影院成人在线观看| 国内精品一区二区在线观看| 亚洲av中文字字幕乱码综合| 日本黄色视频三级网站网址| 国产极品精品免费视频能看的| 亚洲精品一区av在线观看| 美女黄网站色视频| 日本a在线网址| 99热这里只有是精品在线观看 | 偷拍熟女少妇极品色| 少妇人妻一区二区三区视频| 国产老妇女一区| 亚洲人成网站高清观看| 夜夜夜夜夜久久久久| 色播亚洲综合网| 一进一出好大好爽视频| 精品人妻1区二区| 熟妇人妻久久中文字幕3abv| 亚洲av电影不卡..在线观看| а√天堂www在线а√下载| 亚洲 国产 在线| 真人做人爱边吃奶动态| 亚洲av不卡在线观看| 亚洲一区高清亚洲精品| 级片在线观看| 一个人看的www免费观看视频| 国产三级中文精品| 欧美日本亚洲视频在线播放| a级毛片a级免费在线| 久久午夜福利片| 99热这里只有是精品在线观看 | 成人三级黄色视频| 日韩精品青青久久久久久| 亚洲 国产 在线| 亚洲成人免费电影在线观看| 欧美高清成人免费视频www| 国产亚洲欧美98| 一个人免费在线观看的高清视频| 亚洲va日本ⅴa欧美va伊人久久| 亚洲人成伊人成综合网2020| 国产精品三级大全| 国产免费av片在线观看野外av| 国产午夜福利久久久久久| 亚洲成av人片在线播放无| 99久久精品一区二区三区| 国产精品98久久久久久宅男小说| av国产免费在线观看| 少妇的逼好多水| 欧美在线一区亚洲| 亚洲av电影不卡..在线观看| 久久久久国产精品人妻aⅴ院| 欧美+亚洲+日韩+国产| 宅男免费午夜| 精品久久久久久久末码| av黄色大香蕉| 看十八女毛片水多多多| 久久热精品热| 亚洲熟妇熟女久久| 久久九九热精品免费| 中出人妻视频一区二区| 91狼人影院| 精品一区二区三区人妻视频| 高清日韩中文字幕在线| 男人的好看免费观看在线视频| 久久草成人影院| 人妻久久中文字幕网| av在线老鸭窝| 三级男女做爰猛烈吃奶摸视频| 精品久久久久久,| 久久精品影院6| www.www免费av| 国产v大片淫在线免费观看| 午夜日韩欧美国产| 中文字幕av成人在线电影| 久久久久久久久大av| 色播亚洲综合网| 蜜桃亚洲精品一区二区三区| 搞女人的毛片| 亚洲精品日韩av片在线观看| 最近中文字幕高清免费大全6 | 国产91精品成人一区二区三区| 我要搜黄色片| 国产日本99.免费观看| 脱女人内裤的视频| 99久久久亚洲精品蜜臀av| 宅男免费午夜| 亚洲欧美日韩高清专用| 国产精品电影一区二区三区| 国内少妇人妻偷人精品xxx网站| 亚洲综合色惰| 色精品久久人妻99蜜桃| 日本三级黄在线观看| 色5月婷婷丁香| 91久久精品电影网| a级毛片免费高清观看在线播放| 久久性视频一级片| 可以在线观看的亚洲视频| 国产大屁股一区二区在线视频| 久久久久国产精品人妻aⅴ院| 自拍偷自拍亚洲精品老妇| 国产亚洲欧美在线一区二区| 1024手机看黄色片| 欧美最黄视频在线播放免费| 国产高潮美女av| 国产精品免费一区二区三区在线| 亚洲最大成人手机在线| 观看免费一级毛片| 国产精品久久久久久久电影| 亚洲自偷自拍三级| 十八禁人妻一区二区| 国产又黄又爽又无遮挡在线| 一个人看视频在线观看www免费| 动漫黄色视频在线观看| 中文字幕人成人乱码亚洲影| 日本 欧美在线| 99久国产av精品| 国产探花在线观看一区二区| 全区人妻精品视频| 久久人人精品亚洲av| 宅男免费午夜| 中国美女看黄片| 精品久久久久久成人av| 黄色日韩在线| 一进一出好大好爽视频| 国产精品久久视频播放| 国产精品久久久久久久久免 | 久久久久国产精品人妻aⅴ院| 人人妻人人澡欧美一区二区| 乱人视频在线观看| 欧美色欧美亚洲另类二区| av专区在线播放| 国产极品精品免费视频能看的| 黄片小视频在线播放| 一个人观看的视频www高清免费观看| 桃色一区二区三区在线观看| 国产av麻豆久久久久久久| 国产极品精品免费视频能看的| 免费av毛片视频| 日本 欧美在线| 亚洲av电影不卡..在线观看| 成人国产综合亚洲| eeuss影院久久| 极品教师在线免费播放| 18美女黄网站色大片免费观看| 欧洲精品卡2卡3卡4卡5卡区| 欧美日韩综合久久久久久 | 熟妇人妻久久中文字幕3abv| 黄色日韩在线| 国语自产精品视频在线第100页| 十八禁国产超污无遮挡网站| 好男人在线观看高清免费视频| 久久九九热精品免费| 欧美乱妇无乱码| 国产亚洲精品久久久com| 高清日韩中文字幕在线| 少妇被粗大猛烈的视频| 婷婷色综合大香蕉| 每晚都被弄得嗷嗷叫到高潮| 欧美绝顶高潮抽搐喷水| 黄色日韩在线| 男人舔奶头视频| 久久热精品热| 真人做人爱边吃奶动态| av福利片在线观看| 国产精品影院久久| 亚洲av免费在线观看| 一个人观看的视频www高清免费观看| 精品人妻熟女av久视频| 国产精品久久电影中文字幕| 国产激情偷乱视频一区二区| 真实男女啪啪啪动态图| 午夜日韩欧美国产| 一进一出抽搐gif免费好疼| 久久久成人免费电影| 男女床上黄色一级片免费看| 亚洲欧美日韩卡通动漫| 久久人妻av系列| 免费在线观看影片大全网站| 我要搜黄色片| 蜜桃亚洲精品一区二区三区| 内射极品少妇av片p| 热99在线观看视频| 成人无遮挡网站| 国产aⅴ精品一区二区三区波| 国产高清视频在线播放一区| 我的女老师完整版在线观看| 男人的好看免费观看在线视频| 老女人水多毛片| 免费无遮挡裸体视频| 亚洲真实伦在线观看| 午夜福利在线观看免费完整高清在 | 九色成人免费人妻av| 久久久久久久久大av| 美女大奶头视频| 亚洲欧美日韩无卡精品| 波野结衣二区三区在线| 日韩 亚洲 欧美在线| 十八禁人妻一区二区| 热99在线观看视频| 亚洲av美国av| 一二三四社区在线视频社区8| 国产真实乱freesex| 免费看日本二区| 欧美中文日本在线观看视频| 乱码一卡2卡4卡精品| 国产老妇女一区| 人妻夜夜爽99麻豆av| 久久久久久久久久黄片| 午夜激情欧美在线| 欧美黄色淫秽网站| 又紧又爽又黄一区二区| 变态另类丝袜制服| 亚洲av一区综合| 首页视频小说图片口味搜索| 国产视频一区二区在线看| 日本精品一区二区三区蜜桃| 日本一二三区视频观看| 亚洲三级黄色毛片| 国产野战对白在线观看| 久久午夜福利片| 免费观看精品视频网站| 美女大奶头视频| 色精品久久人妻99蜜桃| 亚洲三级黄色毛片| 日本五十路高清| 一级毛片久久久久久久久女| 日本黄色片子视频| 亚洲,欧美,日韩| eeuss影院久久| 禁无遮挡网站| 琪琪午夜伦伦电影理论片6080| 国产成人影院久久av| 中文字幕高清在线视频| 一区福利在线观看| 国产精品三级大全| 男插女下体视频免费在线播放| 国产又黄又爽又无遮挡在线| 成年女人看的毛片在线观看| 人人妻,人人澡人人爽秒播| 久久久精品欧美日韩精品| 亚洲精品粉嫩美女一区| 嫩草影院精品99| 亚洲av不卡在线观看| 黄色一级大片看看| 国产精品国产高清国产av| 成人毛片a级毛片在线播放| 亚洲av电影在线进入| 麻豆av噜噜一区二区三区| 国产免费男女视频| 亚洲av.av天堂| 国产一区二区激情短视频| 高清毛片免费观看视频网站| 熟女人妻精品中文字幕| a级毛片a级免费在线| 99久久久亚洲精品蜜臀av| 少妇被粗大猛烈的视频| 日本精品一区二区三区蜜桃| 一区二区三区免费毛片| 中出人妻视频一区二区| 欧美激情在线99| 天堂动漫精品| 亚洲欧美日韩高清专用| 亚洲成人久久爱视频| 国产在视频线在精品| 成人美女网站在线观看视频| 人人妻人人看人人澡| 精品99又大又爽又粗少妇毛片 | 美女xxoo啪啪120秒动态图 | 极品教师在线视频| 中文在线观看免费www的网站| 亚洲无线观看免费| 91午夜精品亚洲一区二区三区 | 亚洲七黄色美女视频| 男插女下体视频免费在线播放| 亚洲美女搞黄在线观看 | 色av中文字幕| 国产蜜桃级精品一区二区三区| 国产精品野战在线观看| 亚洲电影在线观看av| 国产精品久久久久久亚洲av鲁大| 黄片小视频在线播放| 在现免费观看毛片| 午夜久久久久精精品| 99久久九九国产精品国产免费| 久久精品夜夜夜夜夜久久蜜豆| 精品一区二区免费观看| 亚洲七黄色美女视频| 国产成人aa在线观看| 久久午夜亚洲精品久久| 一个人观看的视频www高清免费观看| 亚洲成av人片免费观看| 无人区码免费观看不卡| 校园春色视频在线观看| 内地一区二区视频在线| 夜夜躁狠狠躁天天躁| 男女做爰动态图高潮gif福利片| 美女高潮喷水抽搐中文字幕| 色播亚洲综合网| 成人亚洲精品av一区二区| av福利片在线观看| 黄片小视频在线播放| 国产精品美女特级片免费视频播放器| 亚洲成av人片在线播放无| 97人妻精品一区二区三区麻豆| 少妇人妻精品综合一区二区 | 女生性感内裤真人,穿戴方法视频| 夜夜躁狠狠躁天天躁| 我的女老师完整版在线观看| 黄色丝袜av网址大全| 精品国产三级普通话版| 蜜桃久久精品国产亚洲av| 亚洲国产精品999在线| 日本黄色片子视频| 美女大奶头视频| 国产成人欧美在线观看| 99国产极品粉嫩在线观看| 一级黄片播放器| 此物有八面人人有两片| 97超级碰碰碰精品色视频在线观看| ponron亚洲| 久99久视频精品免费| 91久久精品电影网| 亚洲自拍偷在线| 国产欧美日韩一区二区精品| 精品一区二区三区人妻视频| 久9热在线精品视频| 麻豆久久精品国产亚洲av| 亚洲黑人精品在线| 在线观看免费视频日本深夜| 给我免费播放毛片高清在线观看| 我的老师免费观看完整版| 真人一进一出gif抽搐免费| 成人特级黄色片久久久久久久| 天天躁日日操中文字幕| 三级男女做爰猛烈吃奶摸视频| 久久天躁狠狠躁夜夜2o2o| av天堂中文字幕网| a在线观看视频网站| 伦理电影大哥的女人| 69人妻影院| 国内精品一区二区在线观看| 中亚洲国语对白在线视频| 又黄又爽又免费观看的视频| 99久久无色码亚洲精品果冻| 变态另类成人亚洲欧美熟女| 黄色视频,在线免费观看| 亚洲av中文字字幕乱码综合| 美女xxoo啪啪120秒动态图 | 国产三级在线视频| 亚洲av熟女| 日日摸夜夜添夜夜添小说| 日韩欧美国产一区二区入口| 热99re8久久精品国产| 国产精品久久久久久久电影| 最后的刺客免费高清国语| 亚洲无线观看免费| 嫩草影院新地址| 国产伦精品一区二区三区四那| 久久人人爽人人爽人人片va | 亚洲av二区三区四区| 欧美zozozo另类| 中文字幕av在线有码专区| 乱人视频在线观看| 亚洲午夜理论影院| 最近中文字幕高清免费大全6 | 国产在视频线在精品| 国产精品一及| 欧美一级a爱片免费观看看| 亚洲18禁久久av| 国产淫片久久久久久久久 | 伦理电影大哥的女人| 亚洲性夜色夜夜综合| 欧美日韩亚洲国产一区二区在线观看| 91av网一区二区| 色综合欧美亚洲国产小说| 国产乱人视频| 国产69精品久久久久777片| av在线蜜桃| 别揉我奶头~嗯~啊~动态视频| 国产精品综合久久久久久久免费| 国产美女午夜福利| 中文字幕av成人在线电影| 好男人在线观看高清免费视频| 国产高潮美女av| 男人的好看免费观看在线视频| 国产午夜福利久久久久久| 性插视频无遮挡在线免费观看| 欧美激情久久久久久爽电影| 国产伦人伦偷精品视频| 亚洲无线观看免费| 偷拍熟女少妇极品色| 欧美绝顶高潮抽搐喷水| 亚洲狠狠婷婷综合久久图片| 91狼人影院| 九九热线精品视视频播放| 国产精品不卡视频一区二区 | a级一级毛片免费在线观看| 九九久久精品国产亚洲av麻豆| 午夜久久久久精精品| 欧美bdsm另类| 久久精品国产亚洲av天美| 搡老岳熟女国产| 蜜桃久久精品国产亚洲av| 97人妻精品一区二区三区麻豆| 亚洲美女搞黄在线观看 | 一个人看视频在线观看www免费| 精品久久国产蜜桃| 真人做人爱边吃奶动态| 免费av毛片视频| 色尼玛亚洲综合影院| 在线免费观看的www视频| 亚洲美女搞黄在线观看 | 欧美高清性xxxxhd video| 久久人人精品亚洲av| 色综合欧美亚洲国产小说| 亚洲自拍偷在线| 成人亚洲精品av一区二区| 色吧在线观看| 午夜免费男女啪啪视频观看 | 黄色女人牲交| 成人一区二区视频在线观看| 观看美女的网站| 婷婷精品国产亚洲av| 美女cb高潮喷水在线观看| 91久久精品电影网| 国产69精品久久久久777片| 热99在线观看视频| 久久伊人香网站| 啪啪无遮挡十八禁网站| 精品日产1卡2卡| 欧美日本亚洲视频在线播放| 熟女电影av网| 搡女人真爽免费视频火全软件 | 成人特级av手机在线观看| 亚洲精品456在线播放app | 99久久成人亚洲精品观看| 不卡一级毛片| 琪琪午夜伦伦电影理论片6080| 淫妇啪啪啪对白视频| 18禁黄网站禁片免费观看直播| 欧美在线一区亚洲| 免费人成在线观看视频色| 亚洲精品成人久久久久久| 桃红色精品国产亚洲av| 色5月婷婷丁香| 十八禁网站免费在线| 国产精华一区二区三区| 人人妻人人看人人澡| 高清日韩中文字幕在线| 久久久国产成人精品二区| 不卡一级毛片| 亚洲欧美清纯卡通| 午夜免费男女啪啪视频观看 | 最近最新中文字幕大全电影3| 色综合亚洲欧美另类图片| 亚洲第一电影网av| 久久欧美精品欧美久久欧美| av在线观看视频网站免费| 毛片一级片免费看久久久久 | 精品久久久久久久久久久久久| 亚洲午夜理论影院| 亚洲内射少妇av| 国产精品亚洲美女久久久| 一区二区三区免费毛片| 免费人成在线观看视频色| 深爱激情五月婷婷| 在线天堂最新版资源| 99热6这里只有精品| 搡女人真爽免费视频火全软件 | 丁香欧美五月| 男女那种视频在线观看| 男人舔女人下体高潮全视频| 精品一区二区三区av网在线观看| 我要看日韩黄色一级片| 亚洲欧美日韩东京热| 欧美成狂野欧美在线观看| 又紧又爽又黄一区二区| 久久午夜亚洲精品久久| 村上凉子中文字幕在线| 夜夜夜夜夜久久久久| 亚洲综合色惰| 少妇的逼水好多| 中文字幕熟女人妻在线| 一级黄色大片毛片| 搡老妇女老女人老熟妇| 欧美区成人在线视频| 国内精品久久久久久久电影| 精品人妻1区二区| av黄色大香蕉| 深夜精品福利| 国内精品久久久久精免费| 天美传媒精品一区二区| 老熟妇乱子伦视频在线观看| 精品一区二区三区视频在线| 18禁黄网站禁片免费观看直播| 欧美日韩亚洲国产一区二区在线观看| 亚洲,欧美精品.| 老师上课跳d突然被开到最大视频 久久午夜综合久久蜜桃 | 波多野结衣高清作品| netflix在线观看网站| 一区二区三区四区激情视频 | 97超视频在线观看视频| 欧洲精品卡2卡3卡4卡5卡区| 长腿黑丝高跟| av在线老鸭窝| 一级毛片久久久久久久久女| 亚洲人成网站在线播放欧美日韩| 亚洲欧美日韩高清专用| 成人三级黄色视频| 久久国产精品人妻蜜桃| or卡值多少钱| 又紧又爽又黄一区二区| 天天躁日日操中文字幕| 永久网站在线| 色av中文字幕| 久久国产乱子免费精品| 级片在线观看| 黄色视频,在线免费观看| 日韩欧美在线乱码| 日韩欧美精品免费久久 | 午夜亚洲福利在线播放| 午夜两性在线视频| 男人的好看免费观看在线视频| 精品久久久久久,| 十八禁国产超污无遮挡网站| 精品久久国产蜜桃| av福利片在线观看| 可以在线观看的亚洲视频| 久久久久亚洲av毛片大全| 久久久久精品国产欧美久久久| 精品一区二区免费观看| 亚洲一区二区三区不卡视频| 国产精品嫩草影院av在线观看 | 人人妻,人人澡人人爽秒播| 老师上课跳d突然被开到最大视频 久久午夜综合久久蜜桃 | 日本一本二区三区精品| 亚洲七黄色美女视频| 日本免费一区二区三区高清不卡| 直男gayav资源| or卡值多少钱| 午夜福利在线观看免费完整高清在 | 夜夜看夜夜爽夜夜摸| 蜜桃亚洲精品一区二区三区| 国产av不卡久久| 国产精品一区二区三区四区免费观看 | 少妇的逼好多水| 婷婷精品国产亚洲av| 非洲黑人性xxxx精品又粗又长| 免费电影在线观看免费观看| 一个人看视频在线观看www免费| 成人av在线播放网站| 我的老师免费观看完整版| 色综合欧美亚洲国产小说| 欧美日韩中文字幕国产精品一区二区三区| 亚洲中文日韩欧美视频| 国产精华一区二区三区| 国产不卡一卡二| 中文字幕人妻熟人妻熟丝袜美| 美女xxoo啪啪120秒动态图 | 两个人视频免费观看高清| 久久久久久九九精品二区国产| 十八禁人妻一区二区| 精品久久久久久,| 日日摸夜夜添夜夜添小说| 在线观看美女被高潮喷水网站 | 欧美精品国产亚洲| 激情在线观看视频在线高清| 亚洲人成伊人成综合网2020| 一a级毛片在线观看| 亚洲中文日韩欧美视频| 欧美乱妇无乱码| 91午夜精品亚洲一区二区三区 | 国产精品亚洲美女久久久| 日韩亚洲欧美综合| 欧美最新免费一区二区三区 | 俺也久久电影网| 丁香欧美五月| 18+在线观看网站| 国产精品自产拍在线观看55亚洲| 麻豆久久精品国产亚洲av| 在线播放国产精品三级| 欧美黄色片欧美黄色片| 亚洲av五月六月丁香网| 欧美3d第一页| 香蕉av资源在线|