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

    基于滑模迭代學(xué)習(xí)的發(fā)動(dòng)機(jī)怠速扭振模擬技術(shù)*

    2016-08-03 00:42:38李文禮石曉輝鄒喜紅王晶晶
    振動(dòng)、測試與診斷 2016年2期
    關(guān)鍵詞:滑模控制

    李文禮, 石曉輝, 柯 堅(jiān), 鄧 斌, 施 全, 鄒喜紅, 王晶晶

    (1.重慶理工大學(xué)汽車零部件先進(jìn)制造技術(shù)教育部重點(diǎn)實(shí)驗(yàn)室 重慶,400054) (2.西南交通大學(xué)機(jī)械工程學(xué)院 成都,610031) (3.重慶市科學(xué)技術(shù)研究院重慶軌道交通及汽車(摩托車)零部件研究中心 重慶,400054)

    ?

    基于滑模迭代學(xué)習(xí)的發(fā)動(dòng)機(jī)怠速扭振模擬技術(shù)*

    李文禮1,石曉輝1,柯堅(jiān)2,鄧斌2,施全1,鄒喜紅1,王晶晶3

    (1.重慶理工大學(xué)汽車零部件先進(jìn)制造技術(shù)教育部重點(diǎn)實(shí)驗(yàn)室重慶,400054) (2.西南交通大學(xué)機(jī)械工程學(xué)院成都,610031) (3.重慶市科學(xué)技術(shù)研究院重慶軌道交通及汽車(摩托車)零部件研究中心重慶,400054)

    摘要為了解決永磁同步電機(jī)在室內(nèi)傳動(dòng)系試驗(yàn)臺(tái)架上模擬發(fā)動(dòng)機(jī)扭振高動(dòng)態(tài)響應(yīng)特性的問題,首先建立了由平均扭矩、往復(fù)慣量引起的激振和燃燒壓力引起的激振信號(hào)構(gòu)成的發(fā)動(dòng)機(jī)實(shí)時(shí)扭振模型,具有數(shù)學(xué)計(jì)算量小、動(dòng)行速度快的特點(diǎn),能滿足高動(dòng)態(tài)控制響應(yīng)的要求;然后,設(shè)計(jì)了一種滑模變結(jié)構(gòu)控制(sliding mode control,簡稱SM)和迭代學(xué)習(xí)控制(iterative learing control,簡稱ILC)相結(jié)合的控制器,利用變結(jié)構(gòu)控制對(duì)系統(tǒng)參數(shù)和外部干擾具有不變性及迭代學(xué)習(xí)不依賴于系統(tǒng)精確數(shù)學(xué)模型的優(yōu)勢,提高了控制系統(tǒng)的穩(wěn)定性和動(dòng)態(tài)響應(yīng);最后,在室內(nèi)傳動(dòng)系臺(tái)架上進(jìn)行了模擬試驗(yàn)。試驗(yàn)結(jié)果驗(yàn)證了發(fā)動(dòng)機(jī)模型和控制算法在傳動(dòng)系試驗(yàn)臺(tái)架上模擬發(fā)動(dòng)機(jī)扭振的有效性和可行性。

    關(guān)鍵詞發(fā)動(dòng)機(jī)模型; 滑??刂疲?迭代學(xué)習(xí); 扭矩脈動(dòng)模擬

    引言

    傳動(dòng)系作為汽車的重要組成部分,隨著設(shè)計(jì)品質(zhì)要求的不斷提高,特別是自動(dòng)變速器的開發(fā)以及傳動(dòng)系性能(noise,vibration,harshness,簡稱NVH)測試,要求汽車傳動(dòng)系試驗(yàn)臺(tái)能夠?yàn)楸辉嚥考峁┡c實(shí)車運(yùn)行一樣的動(dòng)態(tài)工況。在傳動(dòng)系開發(fā)的初始階段,通常情況下與其匹配的發(fā)動(dòng)機(jī)也在開發(fā)測試中,因此很難獲得發(fā)動(dòng)機(jī)作為驅(qū)動(dòng)在室內(nèi)臺(tái)架上對(duì)傳動(dòng)系進(jìn)行耐久和NVH等測試。利用高動(dòng)態(tài)測功機(jī)代替發(fā)動(dòng)機(jī)進(jìn)行室內(nèi)臺(tái)架傳動(dòng)系試驗(yàn)具有安全環(huán)保、試驗(yàn)方便等諸多優(yōu)點(diǎn)[1]。

    發(fā)動(dòng)機(jī)扭矩的周期性波動(dòng)會(huì)引起傳動(dòng)系的扭轉(zhuǎn)振動(dòng),從而產(chǎn)生噪聲,降低零部件的使用壽命。如果變速器齒輪設(shè)計(jì)不當(dāng),發(fā)動(dòng)機(jī)怠速工況下的周期扭矩脈動(dòng)也會(huì)引起變速器的敲擊聲,因此如果要在室內(nèi)臺(tái)架上復(fù)現(xiàn)因發(fā)動(dòng)機(jī)周期扭轉(zhuǎn)振動(dòng)引起的NVH問題,就需要一個(gè)高精度的扭振模型(模擬活塞連桿的慣量扭矩及各缸的燃燒扭矩)和高性能的模擬器。發(fā)動(dòng)機(jī)扭振模擬系統(tǒng)包括一個(gè)低慣量高動(dòng)態(tài)特性的電力測功機(jī)及一個(gè)高頻率響應(yīng)的扭矩模擬器。為了獲得高逼真的扭矩脈動(dòng)波形,采用高頻液壓控制器和自適應(yīng)閉環(huán)控制的方法[2],可對(duì)3~10缸發(fā)動(dòng)機(jī)進(jìn)行模擬并獲得了較好的效果,但其液壓控制系統(tǒng)維護(hù)較為困難。文獻(xiàn)[3]利用直流電機(jī)動(dòng)態(tài)跟蹤航空渦輪發(fā)動(dòng)機(jī)在扭矩載荷擾動(dòng)下的轉(zhuǎn)速軌跡,分析了模擬器的動(dòng)態(tài)特性和幅頻特性。文獻(xiàn)[4]從發(fā)動(dòng)機(jī)電子控制單元(electronic control unit,簡稱ECU)獲得平均扭矩,利用測功機(jī)對(duì)發(fā)動(dòng)機(jī)的動(dòng)態(tài)特性進(jìn)行了模擬,對(duì)比分析了測功機(jī)和被模擬發(fā)動(dòng)機(jī)在相同工況下的轉(zhuǎn)速和扭矩特性,但該文獻(xiàn)并未詳細(xì)描述發(fā)動(dòng)機(jī)扭矩脈動(dòng)的模擬技術(shù)。為了在傳動(dòng)系臺(tái)架上復(fù)現(xiàn)因發(fā)動(dòng)機(jī)在怠速工況下的周期扭矩脈動(dòng)引起的傳動(dòng)系NVH問題(如變速器的敲擊聲),筆者研制了發(fā)動(dòng)機(jī)扭振模擬系統(tǒng),通過使用高動(dòng)態(tài)電力測功機(jī)作為驅(qū)動(dòng)系統(tǒng)來模擬發(fā)動(dòng)機(jī)在不同缸數(shù)、不同曲軸和飛輪轉(zhuǎn)動(dòng)慣量下的扭振特性,利用滑模變結(jié)構(gòu)迭代學(xué)習(xí)控制方法提高了控制系統(tǒng)的穩(wěn)定性和動(dòng)態(tài)響應(yīng)。為了模擬實(shí)際發(fā)動(dòng)機(jī)運(yùn)行的周期扭矩脈動(dòng)工況,電力測功機(jī)需要滿足以下要求:

    1) 快速響應(yīng);

    2) 較小的轉(zhuǎn)動(dòng)慣量;

    3) 電力測功機(jī)功率、最大轉(zhuǎn)速和最大轉(zhuǎn)矩等能夠滿足被模擬發(fā)動(dòng)機(jī)的要求;

    4) 為了滿足試驗(yàn)臺(tái)的通用性,所建模型能夠?qū)Σ煌愋偷陌l(fā)動(dòng)機(jī)進(jìn)行模擬,如不同缸數(shù)、不同曲軸和飛輪轉(zhuǎn)動(dòng)慣量、壓縮比等;

    5) 模擬系統(tǒng)要能夠快速跟蹤發(fā)動(dòng)機(jī)模型輸出信號(hào),很多控制方法是以建立精確的控制系統(tǒng)模型為前提的,而獲得精確的臺(tái)架控制模型較為困難。

    筆者使用基于滑模變結(jié)構(gòu)的迭代學(xué)習(xí)控制方法,不需要精確的控制數(shù)學(xué)模型,通過迭代學(xué)習(xí)即可達(dá)到較高的扭振曲線跟隨精度。

    1發(fā)動(dòng)機(jī)模型分析

    發(fā)動(dòng)機(jī)扭矩周期脈動(dòng)頻率很高,如一個(gè)4缸2沖程發(fā)動(dòng)機(jī)怠速轉(zhuǎn)速為900 r/min,則其點(diǎn)火頻率為30 Hz。經(jīng)測試發(fā)現(xiàn),一個(gè)完整的周期扭振波形至少需要12個(gè)點(diǎn)才能完整地再現(xiàn)出來,所以要對(duì)轉(zhuǎn)速為900 r/min的發(fā)動(dòng)機(jī)周期扭矩振動(dòng)進(jìn)行模擬,其通信頻率或控制周期至少要達(dá)到30×12=360 Hz,如果發(fā)動(dòng)機(jī)轉(zhuǎn)速為6000 r/min,則其通信頻率或控制周期應(yīng)達(dá)到2 400 Hz,發(fā)動(dòng)機(jī)模型的計(jì)算速度應(yīng)滿足實(shí)時(shí)控制要求。精確的發(fā)動(dòng)機(jī)模型在很多文獻(xiàn)中都有論述,如黑箱模型[5]和基于物理結(jié)構(gòu)的模型[6]等,由于這些模型需要大量的發(fā)動(dòng)機(jī)開發(fā)數(shù)據(jù)或詳細(xì)的發(fā)動(dòng)機(jī)設(shè)計(jì)參數(shù),這些參數(shù)在傳動(dòng)系試驗(yàn)中不易獲得,而且傳統(tǒng)的物理模型在模擬發(fā)動(dòng)機(jī)燃燒特性時(shí)復(fù)雜的數(shù)學(xué)計(jì)算很難滿足高動(dòng)態(tài)控制響應(yīng)的要求,因此這些模型不適用于傳動(dòng)系臺(tái)架的動(dòng)態(tài)模擬。

    忽略摩擦扭矩,發(fā)動(dòng)機(jī)扭矩即曲軸的轉(zhuǎn)動(dòng)扭矩,包括曲軸、連桿、活塞等機(jī)械系統(tǒng)產(chǎn)生的慣性扭矩和化學(xué)燃燒產(chǎn)生的燃燒扭矩[7-8]。為了使臺(tái)架具有高動(dòng)態(tài)的模擬效果,首先利用一個(gè)發(fā)動(dòng)機(jī)“轉(zhuǎn)速-扭矩”發(fā)動(dòng)機(jī)外特性圖為控制器提供平均扭矩給定信號(hào),同時(shí)利用數(shù)學(xué)推導(dǎo)和經(jīng)驗(yàn)公式來獲得發(fā)動(dòng)機(jī)慣量扭矩和燃燒扭矩,即在對(duì)發(fā)動(dòng)機(jī)進(jìn)行模擬過程中,把扭矩信號(hào)分離為平均扭矩信號(hào)和扭振信號(hào),而扭振信號(hào)又分解為發(fā)動(dòng)機(jī)往復(fù)慣量引起的激振和燃燒壓力引起的激振。這里的平均扭矩是各缸曲軸、連桿、活塞等機(jī)械系統(tǒng)產(chǎn)生的慣性扭矩和化學(xué)燃燒產(chǎn)生的燃燒扭矩相互疊加平均后的扭矩,以下提到的慣性扭矩和燃燒扭矩為實(shí)際發(fā)動(dòng)機(jī)扭矩去除平均扭矩后的波形,其疊加后的值為零。這種發(fā)動(dòng)機(jī)模型與傳統(tǒng)的發(fā)動(dòng)機(jī)模型相比,數(shù)學(xué)計(jì)算量小,能夠較好滿足高動(dòng)態(tài)控制響應(yīng)的要求。

    由于發(fā)動(dòng)機(jī)平均扭矩信號(hào)可以通過傳統(tǒng)方法獲得,疊加在平均扭矩信號(hào)上的扭矩脈動(dòng)信號(hào)的平均值應(yīng)為零,如圖1所示。

    圖1 發(fā)動(dòng)機(jī)扭矩組成Fig.1 De-construction of the engine torque pulse into different components

    總的發(fā)動(dòng)機(jī)扭矩[7-8]可以用下式表示

    (1)

    其中:Ten為總的發(fā)動(dòng)機(jī)扭矩;TM為平均扭矩;TI為往復(fù)慣量引起的扭矩;TC為燃燒壓力引起的扭矩。

    1.1發(fā)動(dòng)機(jī)平均扭矩模擬

    發(fā)動(dòng)機(jī)平均扭矩可由多種方法得到,包括經(jīng)驗(yàn)公式法、曲線擬合法和查表法,它是以油門開度和發(fā)動(dòng)機(jī)轉(zhuǎn)速為輸入,以發(fā)動(dòng)機(jī)輸出轉(zhuǎn)矩為輸出的模型,圖2為發(fā)動(dòng)機(jī)外特性曲線圖。

    圖2 發(fā)動(dòng)機(jī)外特性曲線圖Fig.2 Engine external characteristic diagram

    1.2往復(fù)慣量引起的扭矩TI

    發(fā)動(dòng)機(jī)運(yùn)動(dòng)部件的離心慣性力,在一定轉(zhuǎn)速下其值大小不變,而且作用方向始終通過回轉(zhuǎn)中心,因此它不會(huì)引起扭轉(zhuǎn)振動(dòng)。往復(fù)慣性力如同氣缸內(nèi)氣體壓力一樣,通過連桿作用在曲柄銷上,亦即在曲軸上產(chǎn)生周期性變化的力矩,從而引起軸系的扭振,發(fā)動(dòng)機(jī)活塞結(jié)構(gòu)如圖3所示。

    圖3 發(fā)動(dòng)機(jī)活塞結(jié)構(gòu)Fig.3 Engine piston structure diagram

    由往復(fù)運(yùn)動(dòng)部件所產(chǎn)生的往復(fù)慣性力為

    (2)

    其中:m為活塞和連桿的總質(zhì)量;x為活塞行程。

    x可以由式(3)來表示

    (3)

    其中:L為連桿長度;r為曲柄半徑;α為曲柄轉(zhuǎn)角;β為連桿與曲柄銷中心與曲軸旋轉(zhuǎn)中心線的夾角;λ為曲柄半徑與連桿長度比。

    根據(jù)運(yùn)動(dòng)學(xué)可知,活塞加速度近似公式

    (4)

    其中:ωc為曲柄旋轉(zhuǎn)角速度。

    因此,往復(fù)慣性力為

    (5)

    由往復(fù)慣性力所產(chǎn)生的力矩為

    (6)

    1.3燃燒壓力引起的扭矩TC

    燃燒扭矩是氣缸內(nèi)氣體燃燒膨脹對(duì)活塞引起的力傳遞到曲軸和飛輪上,燃燒壓力[9]由式(7)得出

    (7)

    其中:Pcomp(αi)為進(jìn)氣支管壓力;Vd(αi)為氣缸容積;V(α)為壓縮過程中氣缸容積(瞬時(shí)容積);k為膨脹系數(shù)。

    瞬時(shí)氣缸工作容積[10]為

    (8)

    曲柄受到的單缸燃燒扭矩為活塞受到的燃燒力乘以活塞行程隨曲柄轉(zhuǎn)角的變化率[9],可以由式(9)得出

    (9)

    上述單缸狀態(tài)下發(fā)動(dòng)機(jī)的往復(fù)慣量力和燃燒壓力引起的扭矩算法模型,對(duì)于多缸發(fā)動(dòng)機(jī),求取出單缸的扭矩后,根據(jù)缸數(shù)和發(fā)火順序進(jìn)行疊加即可獲得發(fā)動(dòng)機(jī)運(yùn)行時(shí)總的信號(hào)輸出。4缸發(fā)動(dòng)機(jī)的燃燒扭矩信號(hào)疊加如圖4所示。

    圖4 發(fā)動(dòng)機(jī)燃燒扭矩Fig.4 The engine combustion torque

    2臺(tái)架系統(tǒng)數(shù)學(xué)模型

    2.1臺(tái)架動(dòng)力學(xué)模型

    汽車傳動(dòng)系試驗(yàn)臺(tái)結(jié)構(gòu)如圖5和圖6所示,被試變速器輸入端由高帶寬低慣量158 kW三相交流永磁同步電機(jī)驅(qū)動(dòng)(圖6中位置1),輸出端連接高動(dòng)態(tài)交流異步電機(jī)(圖6中位置2,3為變速器輸出軸與異步電機(jī)連接的軸和軸承座,異步電機(jī)在NVH室外),同時(shí)在被試變速器(圖6中位置4)輸入輸出端分別連接轉(zhuǎn)速扭矩傳感器用來采集測試數(shù)據(jù),電機(jī)的轉(zhuǎn)速閉環(huán)回饋信號(hào)由安裝在電機(jī)上的編碼器獲得。發(fā)動(dòng)機(jī)模型控制器與永磁交流電機(jī)的通訊速率是決定臺(tái)架模擬動(dòng)態(tài)特性的關(guān)鍵因素,為了滿足快速響應(yīng)的閉環(huán)控制要求,在進(jìn)行測試之前,上位機(jī)把直接編輯好的發(fā)動(dòng)機(jī)模型下載到永磁同步電機(jī)驅(qū)動(dòng)器中,同時(shí)發(fā)動(dòng)機(jī)模型的轉(zhuǎn)速閉環(huán)信號(hào)也由編碼器直接回饋給驅(qū)動(dòng)橋控制模塊,這樣避免了利用第三方控制器接收回饋信號(hào)再傳輸給電機(jī)驅(qū)動(dòng)器方法帶來的延時(shí),提高了控制響應(yīng)速度。

    圖5 汽車傳動(dòng)系試驗(yàn)臺(tái)Fig.5 Schematic diagram of test bench for vehicle drivetrain

    圖6 汽車傳動(dòng)系試驗(yàn)臺(tái)實(shí)物圖Fig.6 Diagram of vehicle drivetrain test bench

    根據(jù)圖5中系統(tǒng)結(jié)構(gòu),驅(qū)動(dòng)電機(jī)帶動(dòng)離合器、變速器和加載測功機(jī)運(yùn)動(dòng),忽略編碼器慣量等因素的影響,臺(tái)架系統(tǒng)由牛頓定律得

    (10)

    其中:Te為驅(qū)動(dòng)電機(jī)輸出扭矩;TL為變速器輸出端扭矩;Je為驅(qū)動(dòng)電機(jī)轉(zhuǎn)動(dòng)慣量;Jt為各擋位下變速器轉(zhuǎn)動(dòng)慣量(文中均以三擋為例);Jc為聯(lián)軸器、半軸和傳感器轉(zhuǎn)動(dòng)慣量;Jd為交流異步電機(jī)轉(zhuǎn)動(dòng)慣量;B為臺(tái)架摩擦阻尼系數(shù)。

    進(jìn)行拉氏變換得臺(tái)架系統(tǒng)傳遞函數(shù)

    (11)

    2.2模擬系統(tǒng)控制模型

    模擬系統(tǒng)采用低慣量永磁同步電機(jī)(pemanent magnet synchronous motor,簡稱PMSM)作為模擬執(zhí)行器件,基于電流解耦的矢量控制,構(gòu)建電流、速度雙內(nèi)環(huán)。為了提高系統(tǒng)內(nèi)環(huán)的抗干擾能力,采用前饋補(bǔ)償?shù)姆绞健9P者直接將模擬系統(tǒng)的轉(zhuǎn)矩傳感器信號(hào)加入到電磁轉(zhuǎn)矩內(nèi)環(huán)的給定中,所以無需設(shè)計(jì)負(fù)載扭矩傳感器。對(duì)永磁同步電機(jī)數(shù)學(xué)模型作如下處理:a.假設(shè)轉(zhuǎn)子永磁磁場在氣隙空間分布為正弦波,定子電樞繞組中的感應(yīng)電動(dòng)勢也為正弦波;b.忽略電子鐵心飽和,認(rèn)為磁路為線性,電感參數(shù)不變;c.不計(jì)鐵心渦流與磁滯損耗;d.轉(zhuǎn)子上無阻尼繞組。根據(jù)文獻(xiàn)[11],采用id=0的PMSM轉(zhuǎn)子磁場定向控制,電壓方程如下

    (12)

    PMSM轉(zhuǎn)矩方程為

    (13)

    PMSM運(yùn)動(dòng)方程為

    (14)

    其中:ud,uq,id,iq,Ld,Lq分別為定子電壓、電流和電感在d,q軸上的分量;R為定子電阻;J為轉(zhuǎn)動(dòng)慣量;Te,TL分別為電磁轉(zhuǎn)矩和負(fù)載轉(zhuǎn)矩;ψa為永磁體與定子交鏈磁鏈;ωe=npω,ωe為電機(jī)電角速度,ω為電機(jī)機(jī)械角速度,np為極對(duì)數(shù);Bm為摩擦阻尼系數(shù)。

    3滑模迭代學(xué)習(xí)控制

    迭代學(xué)習(xí)控制可保證系統(tǒng)在重復(fù)運(yùn)行方向的收斂性,而變結(jié)構(gòu)控制也能保證控制系統(tǒng)在時(shí)間軸方向的收斂性。迭代學(xué)習(xí)顯著的特點(diǎn)是控制算法非常簡單,控制精度很高,理論上可以達(dá)到任意精度跟蹤給定,但主要問題之一是魯棒性問題,而變結(jié)構(gòu)控制對(duì)于系統(tǒng)參數(shù)和外部干擾具有不變性,因而具有強(qiáng)魯棒性。如果能夠設(shè)計(jì)合理的控制律來減小系統(tǒng)抖動(dòng),把滑模變結(jié)構(gòu)算法引入到迭代學(xué)習(xí)控制中能夠提高系統(tǒng)迭代的收斂速度[12-13]。

    3.1滑模變結(jié)構(gòu)控制

    定義永磁同步加載電機(jī)電磁轉(zhuǎn)矩跟蹤誤差為

    (15)

    取加載系統(tǒng)的狀態(tài)變量為

    (16)

    分別對(duì)系統(tǒng)狀態(tài)變量求導(dǎo)得

    (17)

    (18)

    設(shè)計(jì)轉(zhuǎn)矩滑模面

    (19)

    其中:c為正增益。

    對(duì)式(19)s求導(dǎo)得

    (20)

    我國高為炳院士提出的指數(shù)趨律方法不僅能較好地減弱滑模抖振,而且控制量u的求取也比較簡單直觀,其形式如下

    (21)

    由式(21)得控制量iq的表達(dá)式為

    (22)

    建立Lyapunov函數(shù)

    (23)

    由Lyapunov穩(wěn)定性理論可知,要使滑模變結(jié)構(gòu)控制的系統(tǒng)穩(wěn)定需要滿足如下條件

    (24)

    3.2迭代學(xué)習(xí)控制

    (25)

    其中:iq(k)為第k個(gè)迭代周期的控制量;iq(k+1)為第k個(gè)迭代周期的輸出量,并且作為下一個(gè)周期被控對(duì)象的輸入存放在記憶存儲(chǔ)器中;λ為控制器增益;e(k)為系統(tǒng)誤差。

    e(k)可表達(dá)為

    (26)

    從式(13)中電機(jī)扭矩與電流的關(guān)系可得

    (27)

    為了滿足收斂性,系統(tǒng)誤差應(yīng)滿足下式

    (28)

    (29)

    4試驗(yàn)

    本研究試驗(yàn)數(shù)據(jù)為模擬發(fā)動(dòng)機(jī)怠速下的曲線跟隨特性分析,文中曲線均進(jìn)行了平滑濾波處理。圖7中Ten(SM-ILC)為滑模迭代學(xué)習(xí)控制下迭代5次后發(fā)動(dòng)機(jī)瞬態(tài)扭矩控制效果,Ten(r)為模擬目標(biāo)信號(hào)。圖10為其迭代過程中的最大誤差曲線。在初始的5次迭代中由于學(xué)習(xí)參考數(shù)據(jù)不足,系統(tǒng)的跟蹤精度不高。但是從圖8~10中可以看出,模擬系統(tǒng)在學(xué)習(xí)15次后的最大跟蹤誤差為4.9%,而在常規(guī)PID控制下的最大跟蹤誤差為7.6%,滑模迭代學(xué)習(xí)控制方法明顯優(yōu)于常規(guī)PID控制。由圖10可以看出,隨著學(xué)習(xí)次數(shù)的增加,系統(tǒng)控制精度有所提高,但這要犧牲更多的迭代學(xué)習(xí)時(shí)間??傮w上可以看出,利用電力測功機(jī)在汽車傳動(dòng)系試驗(yàn)臺(tái)上進(jìn)行發(fā)動(dòng)機(jī)模擬具有可行性,獲得系統(tǒng)高逼真模擬精度的關(guān)鍵技術(shù)在于建立準(zhǔn)確且能滿足快速響應(yīng)和實(shí)時(shí)控制的發(fā)動(dòng)機(jī)模型以及高動(dòng)態(tài)的硬件控制系統(tǒng)和軟件控制算法。

    圖7 滑模迭代學(xué)習(xí)控制(迭代5次)Fig.7 The curves of the SM-ILC control (5 iterations)

    圖8 滑模迭代學(xué)習(xí)控制與常規(guī)PID控制比較(迭代15次)Fig.8 The curves of the SM-ILC control and PID control (15 iterations)

    圖9 滑模迭代學(xué)習(xí)控制與常規(guī)PID控制比較(局部放大后)Fig.9 The curves of the SM-ILC control and PID control

    圖10 迭代滑模控制迭代15次最大誤差曲線Fig.10 The maximum error curves of the SM-ILC control (15 iterations)

    5結(jié)論

    1) 在測功機(jī)控制中把發(fā)動(dòng)機(jī)扭矩分為平均扭矩、慣性扭矩和燃燒扭矩,得到的發(fā)動(dòng)機(jī)動(dòng)態(tài)扭矩信號(hào)更適合高動(dòng)態(tài)控制響應(yīng)。

    2) 設(shè)計(jì)了滑模變結(jié)構(gòu)迭代學(xué)習(xí)控制器,在該控制算法下進(jìn)行發(fā)動(dòng)機(jī)怠速模擬試驗(yàn),滑模迭代學(xué)習(xí)控制下迭代5次和15次后對(duì)比可以看出,隨著學(xué)習(xí)次數(shù)的增加,系統(tǒng)控制精度有所提高,但這要犧牲更多的迭代學(xué)習(xí)時(shí)間。模擬系統(tǒng)在學(xué)習(xí)15次后的最大跟蹤誤差為4.9%,而在常規(guī)PID控制下的最大跟蹤誤差為7.6%,滑模迭代學(xué)習(xí)控制方法優(yōu)于常規(guī)PID控制。

    3) 從試驗(yàn)數(shù)據(jù)可以看出,筆者采用的發(fā)動(dòng)機(jī)模型和控制算法在汽車傳動(dòng)系試驗(yàn)臺(tái)架上模擬發(fā)動(dòng)機(jī)動(dòng)態(tài)特性是可行的。在系統(tǒng)延時(shí)和滯后的影響下,如何提高系統(tǒng)控制精度是模擬的關(guān)鍵問題,為了達(dá)到更逼真的模擬效果,除了使用動(dòng)態(tài)性能更高的硬件外,還需進(jìn)一步優(yōu)化和改進(jìn)控制算法。

    參考文獻(xiàn)

    [1]劉力康,聞居博,董逾.汽車傳動(dòng)系動(dòng)態(tài)性能試驗(yàn)臺(tái)驅(qū)動(dòng)模擬方案的控制與仿真[J].傳動(dòng)技術(shù),2011,25(30):26-32.

    Liu Likang, Wen Jubo, Dong Yu.The control and simulation of dynamic drive emulation on the automotive transmission dynamometer[J].Drive System Echnique,2011,25(30):26-32. (in Chinese)

    [2]Kaatz S, Abe T, Vanhaaften W, et al. The ford motor company transmission NVH test cell[C]∥Noise & Vibration Conference and Exhibition. Michigan,USA:SAE Internatinal,2003: 681-690.

    [3]Corbett M, Lamm P, McNichols J, et al. Effects of transient power extraction on an integratated hardware-in-the-loop aircraft/propulsion/power system[C]∥Power Systems Conference. Washington,USA:SAE Internatinal,2008: 292-296.

    [4]Newberger N, Nevius T A, Lasota P, et al. Virtual engine dynamometer in service life testing of transmissions:a comparison between real engine and electric dynamometers as prime movers in validation test rigs. extending dynamometer performance for virtual engine simulation[C]∥SAE 2012 World Congress and Exhibition. USA:SAE Internatinal, 2010:919-930.

    [5]Luh G C, Rizzoni G. Identification of a nonlinear MIMO IC engine model during I/M240 driving cycle for on-board diagnosis[C]∥American Control Conference. Maryland, USA:IEEE, 1994:1581-1584.

    [6]Mauer G F. Modeling and experimental validation of torsional crankshaft dynamics [C]∥ International Congress and Exposition.Michigan,USA:SAE Internatinal, 1994:49-57.

    [7]Ayana E, Plahn P, Wejrzanowski K. Active torque cancellation for transmitted vibration reduction of low cylinder count engines[C]∥Vehicle Power and Propulsion Conference. Dearborn, USA:IEEE, 2011: 2971-2977.

    [8]Member G R. Estimate of indicated torque from crankshaft speed fluctuations: a model for the dynamics of the IC engine[C]∥ IEEE Transactions on Vehicular Technology. USA:IEEE Vehicular Technology Society,1989: 168-179.

    [9]Liu Feilong, Gehan A, Amaratuga J, et al. An experimental study on engine dynamics model based in-cylinder pressure estimation[C]∥SAE 2012 World Congress & Exhibition. United States:SAE Internatinal,2012:896-900.

    [10]王浩潔,周斌,展靖華,等.基于MATLAB的示功圖模擬研究[J].內(nèi)燃機(jī),2009,6(6):17-24.

    Wang Haojie, Zhou Bin, Zhan Jinghua, et al. A numerical simulation study of indicator diagrams based on MATLAB[J].Internal Combusion Engines, 2009,6(6):17-24. (in Chinese)

    [11]張曉光,趙克,孫力,等.磁同步電機(jī)滑模變結(jié)構(gòu)調(diào)速系統(tǒng)動(dòng)態(tài)品質(zhì)控制[J].中國電機(jī)工程學(xué)報(bào),2011,31(15):47-52.

    Zhang Xiaoguang, Zhao Ke, Sun Li, et al. Sliding mode control of permanent magnet synchronous motor based on a novel exponential reaching law[J].Proceedings of the CSEE,2011,31(15):47-52. (in Chinese)

    [12]張宏偉,余發(fā)山,卜旭輝,等. 基于魯棒迭代學(xué)習(xí)的永磁直線電機(jī)控制[J].電機(jī)與控制學(xué)報(bào),2012,16(6):81-86.

    Zhang Hongwei, Yu Fashan, Bu Xuhui, et al. Robust iterative learning control for permanent magnet linear motor[J]. Eleetric Machines and Control,2012,16(6): 81-86. (in Chinese)

    [13]何熊熊,秦貞華,張端.基于邊界層的不確定機(jī)器人自適應(yīng)迭代學(xué)習(xí)控制[J].控制理論與應(yīng)用,2012,29(8):1090-1093.

    He Xiongxiong, Qin Zhehua, Zhang Duan. Adaptive iterative learning control for uncertain robot based on boundary layer[J]. Control Theory and Applications, 2012,29(8): 1090-1093. (in Chinese)

    E-mail: liwenli999@163.com

    doi:10.16450/j.cnki.issn.1004-6801.2016.02.025

    收稿日期:2015-07-20;修回日期:2015-08-23

    中圖分類號(hào)TH132.46

    第一作者簡介:李文禮,男,1983年2月生,博士生。主要研究方向?yàn)槠囋囼?yàn)檢測技術(shù)與裝備開發(fā)。曾發(fā)表《關(guān)聯(lián)用戶的汽車傳動(dòng)系載荷譜室內(nèi)臺(tái)架試驗(yàn)編制方法》(《機(jī)械工程學(xué)報(bào)》2014年第50卷第20期)等論文。

    *國家自然科學(xué)基金資助項(xiàng)目(51205432);重慶市科技攻關(guān)計(jì)劃資助項(xiàng)目(cstc2012gg-yyjsB30002)

    猜你喜歡
    滑模控制
    基于干擾觀測器的PID滑模變結(jié)構(gòu)控制
    基于多算法的ROV遠(yuǎn)程協(xié)同控制器設(shè)計(jì)與實(shí)現(xiàn)
    改進(jìn)的動(dòng)態(tài)面船舶航跡跟蹤控制研究
    中國科技博覽(2016年5期)2016-04-23 05:45:45
    油氣懸架車身高度控制研究
    三相電壓型PWM整流器新型雙閉環(huán)控制策略研究
    科技視界(2015年9期)2015-04-07 09:50:54
    航空兵器(2014年4期)2015-01-15 02:39:59
    P—模糊—PI滑??刂葡到y(tǒng)設(shè)計(jì)
    適用于制導(dǎo)控制一體化的反饋線性化滑??刂品椒?/a>
    日韩在线高清观看一区二区三区| 熟妇人妻不卡中文字幕| 成人亚洲精品一区在线观看| 国产亚洲一区二区精品| 最后的刺客免费高清国语| 日本猛色少妇xxxxx猛交久久| 99热网站在线观看| 亚洲av国产av综合av卡| 伦精品一区二区三区| 久久精品国产a三级三级三级| 亚洲不卡免费看| 又粗又硬又长又爽又黄的视频| 亚洲国产av影院在线观看| 日韩成人伦理影院| 国产欧美日韩一区二区三区在线 | 丰满饥渴人妻一区二区三| 午夜老司机福利剧场| 国精品久久久久久国模美| 久久久久人妻精品一区果冻| 观看av在线不卡| 亚洲丝袜综合中文字幕| 在线亚洲精品国产二区图片欧美 | 亚洲人成网站在线观看播放| 国产精品国产三级国产专区5o| 最近中文字幕高清免费大全6| 日韩av在线免费看完整版不卡| 国产精品麻豆人妻色哟哟久久| 午夜精品国产一区二区电影| 亚洲综合精品二区| 观看美女的网站| 国产成人91sexporn| 亚洲欧美一区二区三区黑人 | 亚洲怡红院男人天堂| 精品久久蜜臀av无| 全区人妻精品视频| 久久久精品免费免费高清| 日本欧美视频一区| 在线亚洲精品国产二区图片欧美 | 九九在线视频观看精品| 国产在线视频一区二区| 在线观看人妻少妇| 毛片一级片免费看久久久久| 日本猛色少妇xxxxx猛交久久| 成年美女黄网站色视频大全免费 | 菩萨蛮人人尽说江南好唐韦庄| 日韩熟女老妇一区二区性免费视频| 两个人免费观看高清视频| 欧美激情国产日韩精品一区| 性高湖久久久久久久久免费观看| 丝瓜视频免费看黄片| 91精品国产九色| 在线免费观看不下载黄p国产| 久久97久久精品| 亚洲内射少妇av| 国产熟女午夜一区二区三区 | 人妻人人澡人人爽人人| 啦啦啦视频在线资源免费观看| 丰满少妇做爰视频| 国产成人a∨麻豆精品| 国语对白做爰xxxⅹ性视频网站| 亚洲,一卡二卡三卡| a级毛片免费高清观看在线播放| 丰满乱子伦码专区| 高清av免费在线| 中文精品一卡2卡3卡4更新| 七月丁香在线播放| 极品人妻少妇av视频| 一级片'在线观看视频| 国产高清国产精品国产三级| 久久久久人妻精品一区果冻| 色哟哟·www| 亚洲欧洲精品一区二区精品久久久 | 久久ye,这里只有精品| 久久人人爽人人爽人人片va| 多毛熟女@视频| 18禁在线播放成人免费| av在线播放精品| 国产精品久久久久久av不卡| av电影中文网址| 久久国内精品自在自线图片| 国产 精品1| 欧美丝袜亚洲另类| 国产爽快片一区二区三区| 久久精品国产亚洲av天美| av线在线观看网站| 黑人巨大精品欧美一区二区蜜桃 | 精品少妇久久久久久888优播| 18禁在线播放成人免费| 国产视频首页在线观看| 亚洲内射少妇av| 国产精品三级大全| 精品一区在线观看国产| 综合色丁香网| 最后的刺客免费高清国语| 免费日韩欧美在线观看| 少妇 在线观看| 内地一区二区视频在线| 人妻系列 视频| a级毛色黄片| a级毛色黄片| 亚洲国产欧美日韩在线播放| 美女视频免费永久观看网站| 国产精品蜜桃在线观看| 少妇熟女欧美另类| a级毛色黄片| 狠狠精品人妻久久久久久综合| 免费黄网站久久成人精品| 国产成人精品一,二区| 成人综合一区亚洲| 赤兔流量卡办理| 免费看不卡的av| 18+在线观看网站| 欧美+日韩+精品| 免费黄网站久久成人精品| 99久久人妻综合| 在线观看免费视频网站a站| 韩国av在线不卡| 国产精品女同一区二区软件| 啦啦啦在线观看免费高清www| 免费av中文字幕在线| 国产精品欧美亚洲77777| a级毛片免费高清观看在线播放| 观看美女的网站| 黄片播放在线免费| 国产男女内射视频| 18禁在线播放成人免费| 日日爽夜夜爽网站| 婷婷色综合www| 亚洲精品成人av观看孕妇| 丰满饥渴人妻一区二区三| 国产视频内射| 精品亚洲成国产av| 日韩一区二区视频免费看| 欧美激情国产日韩精品一区| 插逼视频在线观看| 人体艺术视频欧美日本| 午夜福利网站1000一区二区三区| 街头女战士在线观看网站| 777米奇影视久久| 新久久久久国产一级毛片| 国产亚洲av片在线观看秒播厂| 精品一区二区三卡| 亚洲av在线观看美女高潮| 国产亚洲精品第一综合不卡 | 男女无遮挡免费网站观看| 尾随美女入室| 久久女婷五月综合色啪小说| 永久免费av网站大全| 超色免费av| 99热国产这里只有精品6| 国产欧美日韩综合在线一区二区| 少妇丰满av| 一本—道久久a久久精品蜜桃钙片| 久久午夜综合久久蜜桃| 久久97久久精品| 亚洲精品成人av观看孕妇| 精品亚洲成国产av| kizo精华| 亚洲成色77777| av网站免费在线观看视频| freevideosex欧美| av免费观看日本| 亚洲国产av新网站| 黑人猛操日本美女一级片| 亚洲情色 制服丝袜| 美女福利国产在线| 国产精品熟女久久久久浪| 下体分泌物呈黄色| 好男人视频免费观看在线| 日韩亚洲欧美综合| 成年人午夜在线观看视频| 天堂俺去俺来也www色官网| 亚洲美女搞黄在线观看| 免费看光身美女| 精品熟女少妇av免费看| 这个男人来自地球电影免费观看 | av播播在线观看一区| 18禁裸乳无遮挡动漫免费视频| 亚洲国产精品一区三区| 一级毛片我不卡| 韩国高清视频一区二区三区| 极品少妇高潮喷水抽搐| 亚洲,欧美,日韩| 久久人妻熟女aⅴ| 久久久午夜欧美精品| 久久精品国产自在天天线| 边亲边吃奶的免费视频| 男女无遮挡免费网站观看| 精品人妻在线不人妻| 夫妻性生交免费视频一级片| 在线观看免费日韩欧美大片 | 激情五月婷婷亚洲| 日韩伦理黄色片| 日韩人妻高清精品专区| 街头女战士在线观看网站| 国产精品免费大片| 午夜激情久久久久久久| 哪个播放器可以免费观看大片| 成人黄色视频免费在线看| 一本色道久久久久久精品综合| 亚洲国产精品一区三区| 老司机亚洲免费影院| 欧美成人午夜免费资源| 在线播放无遮挡| 99视频精品全部免费 在线| 91午夜精品亚洲一区二区三区| 日本午夜av视频| 国产男人的电影天堂91| 男女高潮啪啪啪动态图| 亚洲国产av影院在线观看| 成人影院久久| 你懂的网址亚洲精品在线观看| 18+在线观看网站| 午夜免费观看性视频| 国产高清不卡午夜福利| 亚洲九九香蕉| 国产精品一区二区在线观看99| 黄片大片在线免费观看| 一区二区三区精品91| 男人操女人黄网站| 99久久国产精品久久久| 高潮久久久久久久久久久不卡| 丰满迷人的少妇在线观看| 午夜免费鲁丝| 成人黄色视频免费在线看| 亚洲人成电影免费在线| 日韩一卡2卡3卡4卡2021年| 最黄视频免费看| 精品国产一区二区三区四区第35| 嫩草影视91久久| 欧美 日韩 精品 国产| 日韩欧美一区视频在线观看| 黄片小视频在线播放| 国产亚洲欧美在线一区二区| 久久中文字幕人妻熟女| 99精品在免费线老司机午夜| 国产成人免费无遮挡视频| 9191精品国产免费久久| 亚洲av成人一区二区三| 亚洲三区欧美一区| 久久久久久久国产电影| 亚洲少妇的诱惑av| 一本一本久久a久久精品综合妖精| 少妇裸体淫交视频免费看高清 | 精品国产乱子伦一区二区三区| 国产男女内射视频| 欧美国产精品va在线观看不卡| 亚洲va日本ⅴa欧美va伊人久久| 757午夜福利合集在线观看| 欧美性长视频在线观看| 久久久久久久大尺度免费视频| 亚洲va日本ⅴa欧美va伊人久久| 9热在线视频观看99| 久热爱精品视频在线9| 变态另类成人亚洲欧美熟女 | 香蕉久久夜色| 国产av国产精品国产| 久久99热这里只频精品6学生| 一级毛片女人18水好多| 亚洲精品成人av观看孕妇| 亚洲专区国产一区二区| 一级黄色大片毛片| 欧美成人午夜精品| 成人国产一区最新在线观看| 丁香欧美五月| 国产伦人伦偷精品视频| 19禁男女啪啪无遮挡网站| 免费高清在线观看日韩| 国产一区二区激情短视频| 久久久久精品人妻al黑| 无限看片的www在线观看| 91大片在线观看| 淫妇啪啪啪对白视频| 午夜日韩欧美国产| 成人精品一区二区免费| a级毛片黄视频| 一区二区三区乱码不卡18| 国产成人av激情在线播放| 天堂8中文在线网| 超碰97精品在线观看| 女人被躁到高潮嗷嗷叫费观| 激情视频va一区二区三区| 97在线人人人人妻| 久久中文字幕一级| 国产熟女午夜一区二区三区| 欧美日韩一级在线毛片| 少妇的丰满在线观看| 国产熟女午夜一区二区三区| 国产免费视频播放在线视频| www.熟女人妻精品国产| 久久ye,这里只有精品| 啦啦啦中文免费视频观看日本| 国产亚洲精品一区二区www | 亚洲精品在线观看二区| 五月天丁香电影| av电影中文网址| 成人特级黄色片久久久久久久 | 国产成人欧美在线观看 | 日韩欧美三级三区| 日韩欧美一区二区三区在线观看 | 成人精品一区二区免费| 欧美日韩国产mv在线观看视频| 欧美精品亚洲一区二区| 久热爱精品视频在线9| 美女视频免费永久观看网站| 午夜福利视频在线观看免费| 99精品欧美一区二区三区四区| 久久久国产成人免费| 99九九在线精品视频| 日本五十路高清| 欧美激情高清一区二区三区| 精品午夜福利视频在线观看一区 | 亚洲国产毛片av蜜桃av| 亚洲精华国产精华精| 久久热在线av| 国产高清视频在线播放一区| 中文字幕另类日韩欧美亚洲嫩草| 老司机福利观看| 一边摸一边做爽爽视频免费| 一区福利在线观看| 久久久久国内视频| 正在播放国产对白刺激| 天堂俺去俺来也www色官网| 高清视频免费观看一区二区| 一进一出抽搐动态| 99国产精品99久久久久| 巨乳人妻的诱惑在线观看| 高清黄色对白视频在线免费看| 欧美国产精品一级二级三级| 两个人看的免费小视频| 飞空精品影院首页| 久久久久网色| 乱人伦中国视频| 成在线人永久免费视频| 制服诱惑二区| 夫妻午夜视频| 女人精品久久久久毛片| 又黄又粗又硬又大视频| 天堂俺去俺来也www色官网| 亚洲少妇的诱惑av| 老汉色av国产亚洲站长工具| 久久天躁狠狠躁夜夜2o2o| 欧美久久黑人一区二区| 在线观看66精品国产| 成人影院久久| 窝窝影院91人妻| 夫妻午夜视频| 丰满人妻熟妇乱又伦精品不卡| 成年版毛片免费区| 美女高潮喷水抽搐中文字幕| www.999成人在线观看| 精品少妇黑人巨大在线播放| 狠狠婷婷综合久久久久久88av| 三上悠亚av全集在线观看| 国产又色又爽无遮挡免费看| 久久久国产精品麻豆| 国产一区二区三区综合在线观看| 午夜福利在线观看吧| 久久天躁狠狠躁夜夜2o2o| 亚洲av美国av| 最近最新免费中文字幕在线| 黄色怎么调成土黄色| 欧美 亚洲 国产 日韩一| av欧美777| 国产亚洲欧美在线一区二区| 一区二区三区精品91| 热99re8久久精品国产| 人人妻人人添人人爽欧美一区卜| 手机成人av网站| 在线观看一区二区三区激情| 一级,二级,三级黄色视频| 欧美一级毛片孕妇| 咕卡用的链子| 国产精品久久久久久精品电影小说| 悠悠久久av| 久久毛片免费看一区二区三区| 大陆偷拍与自拍| 国产欧美日韩综合在线一区二区| 色综合欧美亚洲国产小说| 动漫黄色视频在线观看| 精品国产一区二区三区四区第35| 亚洲国产精品一区二区三区在线| 久久久久久久久免费视频了| 日本五十路高清| 黄片大片在线免费观看| 久久精品国产综合久久久| 中文字幕制服av| 99re在线观看精品视频| 亚洲色图 男人天堂 中文字幕| 法律面前人人平等表现在哪些方面| 王馨瑶露胸无遮挡在线观看| 日本黄色视频三级网站网址 | 美女福利国产在线| 99九九在线精品视频| 日韩欧美一区视频在线观看| 肉色欧美久久久久久久蜜桃| 日本vs欧美在线观看视频| 黄色丝袜av网址大全| 精品久久久精品久久久| 久久久久久久国产电影| 久久精品人人爽人人爽视色| 亚洲综合色网址| 欧美人与性动交α欧美软件| 在线亚洲精品国产二区图片欧美| 国产一区二区三区综合在线观看| 久久精品国产99精品国产亚洲性色 | 亚洲一码二码三码区别大吗| 国产黄频视频在线观看| 亚洲色图 男人天堂 中文字幕| 啦啦啦视频在线资源免费观看| 亚洲午夜精品一区,二区,三区| 精品一品国产午夜福利视频| 高清黄色对白视频在线免费看| 亚洲熟女精品中文字幕| 757午夜福利合集在线观看| 99re在线观看精品视频| 男女之事视频高清在线观看| www.熟女人妻精品国产| 国产精品98久久久久久宅男小说| 欧美成人午夜精品| 曰老女人黄片| av欧美777| 天天添夜夜摸| 国产成人av激情在线播放| 欧美亚洲日本最大视频资源| 超碰97精品在线观看| 国产伦人伦偷精品视频| 亚洲国产av新网站| 国产精品一区二区在线不卡| 999精品在线视频| 亚洲欧美精品综合一区二区三区| 飞空精品影院首页| 精品人妻熟女毛片av久久网站| 欧美乱妇无乱码| 亚洲精品久久午夜乱码| 国产精品麻豆人妻色哟哟久久| 亚洲精品美女久久久久99蜜臀| 最黄视频免费看| 大片电影免费在线观看免费| 国产精品99久久99久久久不卡| 热99久久久久精品小说推荐| 黄片大片在线免费观看| 久久青草综合色| 黑丝袜美女国产一区| 日韩欧美免费精品| 国产成人欧美在线观看 | 91av网站免费观看| 国产精品成人在线| 99精品在免费线老司机午夜| 美女视频免费永久观看网站| 狠狠狠狠99中文字幕| 亚洲五月婷婷丁香| 男女高潮啪啪啪动态图| 窝窝影院91人妻| 亚洲三区欧美一区| 99在线人妻在线中文字幕 | 侵犯人妻中文字幕一二三四区| 精品亚洲成国产av| 亚洲国产精品一区二区三区在线| 麻豆国产av国片精品| 国产精品一区二区在线不卡| 免费在线观看影片大全网站| 性色av乱码一区二区三区2| 不卡一级毛片| 99国产精品免费福利视频| 久久久久久久久久久久大奶| av网站在线播放免费| 亚洲专区国产一区二区| 久久精品国产亚洲av香蕉五月 | 国产av国产精品国产| 一区二区日韩欧美中文字幕| 国产精品久久久久久人妻精品电影 | 99riav亚洲国产免费| 99国产精品一区二区三区| 麻豆乱淫一区二区| 国产一区二区三区综合在线观看| 精品熟女少妇八av免费久了| 国产精品亚洲一级av第二区| 啦啦啦 在线观看视频| 80岁老熟妇乱子伦牲交| 一区在线观看完整版| 精品少妇内射三级| 老汉色∧v一级毛片| 亚洲熟妇熟女久久| 国产精品成人在线| 脱女人内裤的视频| 亚洲精品一二三| 汤姆久久久久久久影院中文字幕| 久久久久久久久免费视频了| 国产精品98久久久久久宅男小说| 亚洲精品中文字幕一二三四区 | 91九色精品人成在线观看| 999久久久精品免费观看国产| 亚洲七黄色美女视频| 免费av中文字幕在线| 丁香欧美五月| 少妇裸体淫交视频免费看高清 | 91国产中文字幕| 黄色片一级片一级黄色片| 一区二区三区激情视频| 欧美激情 高清一区二区三区| 欧美日韩成人在线一区二区| 亚洲人成电影免费在线| 一区二区三区乱码不卡18| 男女午夜视频在线观看| 久久国产精品大桥未久av| 777米奇影视久久| 激情视频va一区二区三区| kizo精华| 亚洲熟女精品中文字幕| 国产人伦9x9x在线观看| 涩涩av久久男人的天堂| 国产成人精品无人区| 国产aⅴ精品一区二区三区波| 激情视频va一区二区三区| 亚洲精品自拍成人| 91精品国产国语对白视频| 97在线人人人人妻| 国产国语露脸激情在线看| 黑人操中国人逼视频| av视频免费观看在线观看| 人妻久久中文字幕网| 亚洲精品中文字幕在线视频| 亚洲,欧美精品.| 成年人黄色毛片网站| 中文亚洲av片在线观看爽 | 建设人人有责人人尽责人人享有的| 精品乱码久久久久久99久播| 欧美日韩视频精品一区| 黄色怎么调成土黄色| 亚洲男人天堂网一区| 国产成人精品久久二区二区91| 人人澡人人妻人| 国产日韩欧美亚洲二区| 色精品久久人妻99蜜桃| 夜夜夜夜夜久久久久| 在线观看www视频免费| 免费在线观看完整版高清| 桃花免费在线播放| 色婷婷久久久亚洲欧美| 肉色欧美久久久久久久蜜桃| 色综合欧美亚洲国产小说| 飞空精品影院首页| 制服诱惑二区| 国产精品久久久久成人av| 亚洲视频免费观看视频| 伦理电影免费视频| 亚洲精品中文字幕在线视频| 香蕉丝袜av| 一级毛片女人18水好多| 国产精品麻豆人妻色哟哟久久| 久久久久国内视频| 亚洲av美国av| 国产单亲对白刺激| 国产精品久久久av美女十八| 亚洲九九香蕉| 免费在线观看日本一区| 精品人妻熟女毛片av久久网站| 成年女人毛片免费观看观看9 | 在线观看舔阴道视频| 巨乳人妻的诱惑在线观看| 老司机午夜福利在线观看视频 | 国产精品久久久久成人av| 午夜福利视频在线观看免费| 一级,二级,三级黄色视频| 亚洲欧美精品综合一区二区三区| 日韩视频在线欧美| 成在线人永久免费视频| 久久午夜综合久久蜜桃| 亚洲欧美色中文字幕在线| 亚洲av欧美aⅴ国产| 欧美精品亚洲一区二区| 亚洲中文av在线| 亚洲成人手机| 精品一区二区三区视频在线观看免费 | 亚洲男人天堂网一区| 久久精品国产99精品国产亚洲性色 | 好男人电影高清在线观看| 日韩欧美免费精品| 窝窝影院91人妻| 考比视频在线观看| 菩萨蛮人人尽说江南好唐韦庄| 亚洲熟女毛片儿| 久久精品国产亚洲av高清一级| 两性午夜刺激爽爽歪歪视频在线观看 | 夫妻午夜视频| 成人国产一区最新在线观看| 国产精品秋霞免费鲁丝片| 在线观看免费视频日本深夜| 欧美成狂野欧美在线观看| 亚洲一区中文字幕在线| www.999成人在线观看| 日日爽夜夜爽网站| 精品一品国产午夜福利视频| 操美女的视频在线观看| 丝袜美足系列| 国产一区有黄有色的免费视频| 久久久久国内视频| 欧美亚洲日本最大视频资源| 亚洲视频免费观看视频| 两人在一起打扑克的视频| 天天躁狠狠躁夜夜躁狠狠躁| 中文字幕另类日韩欧美亚洲嫩草| e午夜精品久久久久久久| 欧美精品啪啪一区二区三区| 99久久99久久久精品蜜桃| 久久精品国产亚洲av香蕉五月 | 久久精品亚洲av国产电影网| 香蕉丝袜av| bbb黄色大片| 日日爽夜夜爽网站| 热re99久久精品国产66热6| 极品教师在线免费播放| 女人爽到高潮嗷嗷叫在线视频| 久久狼人影院|