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

    乘用車穩(wěn)定性自抗擾控制策略

    2016-12-16 11:31:37靳立強
    關(guān)鍵詞:偏角角速度質(zhì)心

    劉 剛, 靳立強, 王 熠

    (1. 吉林大學(xué) 汽車仿真與控制國家重點實驗室,吉林 長春 130025; 2. 河南工學(xué)院 自動控制系, 河南 新鄉(xiāng) 453000;3. 燕山大學(xué) 信息學(xué)院計算機教學(xué)實驗中心, 河北 秦皇島 066004)

    ?

    乘用車穩(wěn)定性自抗擾控制策略

    劉 剛1,2, 靳立強1, 王 熠3

    (1. 吉林大學(xué) 汽車仿真與控制國家重點實驗室,吉林 長春 130025; 2. 河南工學(xué)院 自動控制系, 河南 新鄉(xiāng) 453000;3. 燕山大學(xué) 信息學(xué)院計算機教學(xué)實驗中心, 河北 秦皇島 066004)

    為了提高汽車駕駛過程中的安全性能,研究傳統(tǒng)乘用車的操縱穩(wěn)定性控制策略,并采用分層控制結(jié)構(gòu)設(shè)計穩(wěn)定性控制系統(tǒng).控制策略包含參考模型、UniTire輪胎模型、參數(shù)估計模塊,控制量期望值計算模塊、自抗擾控制(ADRC)部分及直接橫擺力矩分配和滑移率控制模塊.采用動力學(xué)計算的方法估計質(zhì)心側(cè)偏角,將估算得到的質(zhì)心側(cè)偏角和傳感器采集得到的橫擺角速度作為控制變量.當(dāng)控制變量超過設(shè)定門限值時,采用ADRC計算橫擺力矩.基于UniTire輪胎模型建立主動橫擺力矩和車輪滑移率的關(guān)系確定目標(biāo)滑移率,通過變參數(shù)PID實現(xiàn)滑移率控制.應(yīng)用硬件在環(huán)仿真(HILS)平臺分別進行單移線和雙移線仿真實驗,結(jié)果表明:所提出的穩(wěn)定性控制策略能夠提高車輛的操穩(wěn)性能,使得控制變量能夠緊密跟蹤期望值.

    自抗擾控制(ADRC);車輛動力學(xué);電子穩(wěn)定控制系統(tǒng);橫擺力矩控制;硬件在環(huán)仿真(HILS)

    車輛電子穩(wěn)定系統(tǒng)(electronic stability program, ESP)是汽車主動安全技術(shù)領(lǐng)域一個重要的研究方向.ESP用于保證汽車在行駛過程中的穩(wěn)定性,即在車輛出現(xiàn)轉(zhuǎn)向不足或者轉(zhuǎn)向過多危險情況時,ESP通過發(fā)動機轉(zhuǎn)矩與制動轉(zhuǎn)矩產(chǎn)生直接橫擺力矩來糾正車輛行駛狀態(tài),從而保證車輛的穩(wěn)定狀態(tài).ESP產(chǎn)品目前只有博世、大陸等外企能夠生產(chǎn),國內(nèi)尚無相關(guān)產(chǎn)品.關(guān)于ESP系統(tǒng)的大量研究,集中在直接橫擺力矩的控制算法和橫擺力矩分配策略.目前在接橫擺力矩控制算法方面,滑??刂芠1]、模糊控制[2]、模型預(yù)測控制[3]、H∞算法[4]等已有運用.在橫擺力矩分配策略方面,Yang[5]采用二次規(guī)劃法設(shè)計橫擺力矩分配器,采用有效集法求車輪橫擺力矩分配最優(yōu)解.付皓[6]利用轉(zhuǎn)矩矢量進行力矩分配,實車試驗表明,這類方法有效提高了電動車的側(cè)向穩(wěn)定性能.Li等[7]提出了基于車載輪胎模型的橫擺力矩分配和制動壓力控制方案,并討論了輪胎模型不準(zhǔn)確和極限附著橢圓偏差2種情況下的目標(biāo)值選取方法.Chen等[8-10]將橫擺力矩分配和目標(biāo)滑移率值的確定問題轉(zhuǎn)換成了有約束的優(yōu)化問題,使用數(shù)學(xué)方法獲得橫擺力矩分配最優(yōu)解.

    本文采用分層控制策略,將車輛穩(wěn)定控制系統(tǒng)分成上、下兩層.上層控制器計算車輛穩(wěn)定所需直接橫擺力矩,考慮到車輛在復(fù)雜工況下會產(chǎn)生強非線性、模型參數(shù)擾動等控制難題,采用自抗擾控制算法來解決車輛行駛過程中的非線性和不確定性的問題.下層控制器用于橫擺力矩分配,通過研究UniTire輪胎模型建立直接橫擺力矩和輪胎目標(biāo)滑移率之間的關(guān)系,采用自適應(yīng)變參數(shù)比例積分微分(proportion integration differentiation, PID)算法實現(xiàn)滑移率控制.

    1 UniTire輪胎模型

    車輛動力學(xué)性能嚴重依賴輪胎與路面之間的作用力,因此建立高精度的輪胎模型是車輛電子穩(wěn)定研究領(lǐng)域的關(guān)鍵問題,本文采用Guo等[11]提出的UniTire輪胎模型計算車輛在復(fù)雜工況下輪胎受力.

    在UniTire模型中,定義縱向滑移率Sx和側(cè)向滑移率Sy分別為

    (1)

    (2)

    式中:Ω為車輪角速度,Re為車輪有效半徑,vx和vy分別為輪心的縱向速度和側(cè)向速度.

    在UniTire輪胎模型中,無量綱的縱向滑移率φx和無量綱側(cè)向滑移率φy分別為

    (3)

    式中:Kx為縱滑剛度,Ky為側(cè)偏剛度,Fxm和Fym分別為輪胎的縱向力極值和側(cè)向力極值,φ為無量綱總滑移率.

    由式(1)~(3)可得UniTire輪胎的側(cè)向力和縱向力分別為

    (4)

    (5)

    式中:Fz為輪胎z軸方向受力;μx為縱向滑移摩擦系數(shù);μy側(cè)向滑移摩擦系數(shù);E為曲率因子,用來調(diào)節(jié)受力曲線的彎曲度.

    考慮到車輛穩(wěn)定性控制時的加速、制動工況對車輪垂直載荷的影響,車輪載荷為

    (6)

    式中:Lf和Lr分別為車輛前后軸到車輛質(zhì)心的距離,hcg為質(zhì)心到地面的距離,L為前、后軸之間的距離,Bf為車輛寬度.

    2 ESP分層控制策略

    ESP在汽車實際運動狀態(tài)偏離駕駛員期望運動狀態(tài)時主動介入進行穩(wěn)定性控制.本文設(shè)計的電子穩(wěn)定系統(tǒng)采用分層控制策略,上層控制器負責(zé)采用自抗擾控制算法計算穩(wěn)定性控制所需的直接橫擺力矩,下層控制器負責(zé)將計算所得橫擺力矩分配到被控車輪進行滑移率控制從而實現(xiàn)車輛的穩(wěn)定性控制.上層控制器的控制變量選擇質(zhì)心側(cè)偏角和橫擺角速度,因為質(zhì)心側(cè)偏角和橫擺角速度能夠反映車輛的轉(zhuǎn)向特性,其中橫擺角速度用來判斷車輛轉(zhuǎn)向過程中是否出現(xiàn)轉(zhuǎn)向不足或轉(zhuǎn)向過多的情況,質(zhì)心側(cè)偏角用來判斷在轉(zhuǎn)向過程中是否出現(xiàn)軌跡偏離的情況,兩者共同決定了車輛的穩(wěn)定狀況.將這2個變量與其期望值分別作差,當(dāng)其中任何一個的差值大于介入門限值時,就采用自抗擾算法進行橫擺力矩計算.下層控制器首先根據(jù)車輛狀態(tài)判斷選取被控車輪,然后根據(jù)上層計算所得橫擺力矩和UniTire輪胎模型計算得出目標(biāo)滑移率,通過采用自適應(yīng)變參數(shù)PID算法實現(xiàn)滑移率控制,從而使得車輛狀態(tài)達到穩(wěn)定.如圖1所示為控制策略架構(gòu)圖.

    圖1 車輛電子穩(wěn)定性控制策略架構(gòu)圖Fig.1 Architecture of proposed vehicle electronical stability controlscheme

    3 上層控制器設(shè)計

    3.1 理想橫擺角速度和理想質(zhì)心側(cè)偏角計算

    采用車輛線性二自由度模型穩(wěn)態(tài)響應(yīng)模型計算橫擺角速度和質(zhì)心側(cè)偏角的期望值.

    通過文獻[12]所提公式計算得到橫擺角速度的期望值為

    (7)

    其中穩(wěn)定性因數(shù)為

    (8)

    式中:ωrd為汽車等速圓周運動時的橫擺角速度期望值;vx一般采用估算的方法得到,可參考文獻[6];δ為前輪側(cè)偏角,可通過方向盤轉(zhuǎn)角傳感器得到;m汽車的質(zhì)量;Cf為前輪側(cè)偏剛度,Cr為后輪的側(cè)偏剛度.

    由于式(7)未考慮路面工況的復(fù)雜,因此采用文獻[7]中利用附著系數(shù)對期望值進行限制,即橫擺角速度的期望值的取值范圍為

    (9)

    考慮到模型的瞬態(tài)響應(yīng)特性,增加一階滯后環(huán)節(jié)最終得到橫擺角速度的期望值為

    (10)

    式中:τ為目標(biāo)響應(yīng)時間.

    根據(jù)二自由度車輛模型做等速圓周運動的情況,在穩(wěn)態(tài)時,可得理想質(zhì)心側(cè)偏角為

    (11)

    式(11)是根據(jù)側(cè)向力與輪胎的側(cè)偏角成正比而得出的,并未考慮在輪胎附著接近極限值,此時側(cè)向力和輪胎側(cè)偏角已不成線性的關(guān)系.根據(jù)文獻[13-15]中提及的經(jīng)驗結(jié)論,即在路面附著系數(shù)為0.90時,質(zhì)心側(cè)偏角的期望值上線為10°;在路面附著系數(shù)為0.35的工況下,質(zhì)心側(cè)偏角期望值上限值為4°,因此修正的質(zhì)心側(cè)偏角的期望值為

    (12)

    式中:βupper為質(zhì)心側(cè)偏角的上限值:

    |βupper|=arctan (0.02μg).

    (13)

    3.2 質(zhì)心側(cè)偏角估算

    質(zhì)心側(cè)偏角在車輛穩(wěn)定性控制中多采用參數(shù)估計的方法.目前常見的質(zhì)心側(cè)偏角估算研究多采用卡爾曼濾波、非線性狀態(tài)觀測器等方法.此類方法多是以非線性模型為基礎(chǔ),算法復(fù)雜,實時性較差,在實際ESP系統(tǒng)中受限.本文利用線性模型進行質(zhì)心側(cè)偏角估算.

    考慮到在穩(wěn)定性控制中,車輛的質(zhì)心側(cè)偏角上限值僅為10°,是一個較小的角度值,將質(zhì)心側(cè)偏角的公式變化為

    (14)

    本文設(shè)計的電子穩(wěn)定系統(tǒng)需要每10 ms對質(zhì)心側(cè)偏角進行一次估算,在10 ms的時間內(nèi)車輛的縱向車速可以認為是一個常數(shù),因此式(14)的導(dǎo)數(shù)為

    (15)

    汽車理論中線性二自由度車輛模型微分方程為

    (16)

    將式(16)代入式(15)可得質(zhì)心側(cè)偏角的微分值為

    (17)

    式中:ay和ωr分別為側(cè)向加速度和橫擺角速度,兩者均可以通過汽車上的偏航率傳感器采樣得到.對式(17)求積分即可得到質(zhì)心側(cè)偏角.

    3.3 路面附著系數(shù)

    采用Fukada[16]提出的利用非線性度的附著補償系數(shù)加側(cè)向加速度的方式確定路面附著系數(shù)為

    μ=|Gy|+μe.

    (18)

    式中:μ為路面附著系數(shù);Gy=ay/g,g為重力加速度;μe為附著補償系數(shù),該系數(shù)前輪轉(zhuǎn)角的偏差值之間的關(guān)系如圖2所示.其中前輪轉(zhuǎn)角偏差為

    (19)

    式中:vc為特征車速.

    圖2 附著補償系數(shù)與前輪轉(zhuǎn)角差值的關(guān)系Fig.2 Curve of compensation factor versus front wheel angle difference

    3.4 自抗擾控制器設(shè)計

    自抗擾控制器的模型采用線性二自由度模型:

    (20)

    (21)

    車輛在行駛過程中,利用偏航率傳感器采集橫擺角速度,然后將其與橫擺度速度的期望值相減;當(dāng)差值大于系統(tǒng)設(shè)定的門限值時,基于橫擺角速度的自抗擾控制器就會介入控制.基于質(zhì)心側(cè)偏角的自抗擾控制器同理.

    3.3.1 橫擺角速度控制

    將式(21)改寫為以直接橫擺力矩ΔM為控制量的仿射模型:

    (22)

    式中:

    其中,u為輸出的直接橫擺力矩;f10為模型中的已知部分,其中的變量均可通過傳感器測量得到;f11為模型中的未知部分即車輛在運行過程中受到的各種“外擾”和“內(nèi)擾”的總和作用.

    橫擺角速度動態(tài)方程式(式(22))的擴張觀測器為

    (23)

    式中:

    其中,z11為橫擺角速度ωr的觀測器;z12為式(23)中f11的估計值;α1為濾波因子,取值范圍在(0,0.8);m1為非線性因子;k01和k02皆為輸出誤差矯正增益;fal(e,α,m)為坐標(biāo)系下在原點附近的冪次函數(shù).

    式(22)的非線性誤差反饋控制率為

    e2=ωrd-z11,

    uωr=k1fhan(e2,α2,m2).

    (24)

    在自抗擾算法中,k01、k02、k1、α2、m2、b0皆為待調(diào)節(jié)參數(shù).其中,k01、k02、k1是由系統(tǒng)的采樣步長確定的,韓京清[17]提出了與斐波那契數(shù)列相關(guān)的方法為

    (25)

    根據(jù)該方法得到觀測器系數(shù)的范圍后,利用Simulink軟件最終得到k01=800,k02=30 000,k1=18 000;參數(shù)α2為控制量增益,一般取2 000;參數(shù)m2和b0分別類似于PID控制算法中的比例增益系數(shù)和積分增益系數(shù),利用Simulink軟件優(yōu)化確定m2=600,b0=3.

    最終,基于橫擺角速度控制的直接橫擺力矩為

    ΔMωr=Iz(uωr-f10-z12).

    (26)

    3.3.2 質(zhì)心側(cè)偏角控制

    當(dāng)車輛運行時質(zhì)心側(cè)偏角與其期望值的差值達到門限值時,基于質(zhì)心側(cè)偏角的自抗擾控制器就會介入控制.針對式(18)作如下修改,將其改成以直接橫擺力矩ΔM為控制量的一階仿射模型:

    (27)

    式中:

    對基于質(zhì)心側(cè)偏角動態(tài)方程式(23)設(shè)計擴張狀態(tài)觀測器:

    (28)

    式中:z21為質(zhì)心側(cè)偏角β的觀測值,z22為f21的估計值.接著對式(26)設(shè)計非線性狀態(tài)誤差反饋率為

    (29)

    式中:k2、k03、k04、α4、m4、b1為待調(diào)試參數(shù).由式(25)得到k2=600、k03=28 000、k04=15 000、α4=1 700、m4=500、b1=2.

    基于質(zhì)心側(cè)偏角控制的直接橫擺力矩為

    ΔMβ=Iz(z22-uβ).

    (30)

    最終ESP在介入控制時上層控制器輸出的直接橫擺力矩為

    (31)

    式中:eβ-lim和erd-lim分別為質(zhì)心側(cè)偏角和橫擺角速度在ESP介入控制時的門限設(shè)定值.

    4 下層控制器設(shè)計

    4.1 直接橫擺力矩的分配策略

    車輛的轉(zhuǎn)向狀態(tài)可以通過理想橫擺角速度和實際橫擺角速度的差值進行判斷,判斷方法如圖(3)所示.圖中w為橫擺角速度偏差的門限值,當(dāng)偏差值小于門限值時,車輛正處在中性轉(zhuǎn)向狀態(tài).而橫擺角速度差值大于門限值,且當(dāng)ωr-ωrd<0時,車輛處于不足轉(zhuǎn)向;相反的情況下則處于轉(zhuǎn)向過度.

    圖3 車輛轉(zhuǎn)向狀態(tài)估算流程圖Fig.3 Flow chart of vehicle steering state estimation

    直接橫擺力矩的分配部分需要考慮2個問題,一個是被控車輪的選擇,二則是需要確定被控車輪的目標(biāo)滑移率.對于目前將直接橫擺力矩的分配轉(zhuǎn)換成有約束的優(yōu)化問題,目前通過優(yōu)化算法來求出每個車輛的最優(yōu)解[18].這種將橫擺力矩分配到4個車輪的方式可以有效地提高控制效果.但是此類方法算法復(fù)雜,實時性不太好,因此本文盡可能地將直接橫擺力矩分配到單個車輪,如果分配的縱向力大于輪胎的極限值,則可以選擇同側(cè)車輪作為輔助制動輪.當(dāng)車輛行駛過程中出現(xiàn)轉(zhuǎn)向不足時,ESP介入控制選擇內(nèi)側(cè)后輪作為主要的控制對象;當(dāng)車輛出現(xiàn)轉(zhuǎn)向過度的情況時,則選擇外側(cè)前輪作為主要制動對象,以此減小轉(zhuǎn)向過度的趨勢.被控車輪的選擇策略如表1所示.

    表1 被控車輪選擇策略

    圖4 車輛模型中的輪胎受力圖Fig.4 Force diagram of tire in vehicle motion model

    根據(jù)輪胎摩擦圓理論,當(dāng)汽車的車輪在縱向的滑移率增加Δλ,該輪胎的縱向力則會增加ΔFxfj,同時輪胎的側(cè)向力會減少ΔFyfj.如圖(4)所示為車輛模型中輪胎受力圖,當(dāng)ESP介入控制時,前輪有轉(zhuǎn)角,后輪沒有轉(zhuǎn)角,因此直接橫擺力矩分配至前輪與直接橫擺力矩分配至后輪的公式為

    (32)

    式中:δf為前輪轉(zhuǎn)角,df為車輛寬度,lf為前軸到質(zhì)心的距離.

    由式(32)可知,直接橫擺力矩的期望值由輪胎縱向力產(chǎn)生的力矩增量減去側(cè)向力的力矩增量構(gòu)成,接著我們對該式改寫構(gòu)建橫擺力矩的期望值與輪胎滑移率的關(guān)系為

    (33)

    由文獻[11]的UniTire輪胎模型得到縱向力公式和側(cè)向力(式(4)),對式(4)縱向滑移率微分后帶入得式(33),即可確定直接橫擺力矩期望值和輪胎目標(biāo)滑移率增量之間的關(guān)系為

    ΔM=g(μ,Fzi)Δλi.

    (34)

    式中:Δλi為車輪分配得到的直接橫擺力矩而引起的滑移率的變化量,μ為地面附著系數(shù),Fzi為車輪的垂直載荷.式(34)中的路面附著系數(shù)和車輪垂直載荷可以分別根據(jù)式(18)和式(6)計算得出.

    4.2 目標(biāo)滑移率控制

    通過將目標(biāo)滑移率和汽車當(dāng)前車輪的滑移率做差,即e(k)=λd(k)-λ(k),采用變參數(shù)自適應(yīng)PID計算得到車輪的目標(biāo)制動壓力變化率為

    Δp(k)=Kpe(k)+Kie(k)+Kd[e(k)-e(k-1)].

    (35)

    式中:的比例系數(shù)Kp、積分系數(shù)Ki及微分系數(shù)Kd的值均可根據(jù)滑移率差值e(k)進行選取,如表2所示,表2中E1=0.2,E2=0.5,e為滑移率差值e(k)在當(dāng)前周期的值.

    表2 自適應(yīng)比例積分微分(PID)控制參數(shù)值列表

    Tab.2 Parameters of adaptive proportion integration derivative (PID)

    ePIDe>E2300115E1

    5 硬件在環(huán)仿真

    圖5 硬件在環(huán)仿真HIC平臺及車輛電子穩(wěn)定系統(tǒng)(ESP)控制器Fig.5 Image of Hardware-in-the-loop (HIL) and electronic stability program (ESP) controller circuit

    為了驗證本文提出的控制策略,采用基于dSpace硬件在環(huán)平臺(Hardware-in-the-Loop, HIL)進行一系列的實驗.硬件在環(huán)平臺是將部分硬件系統(tǒng)與計算機仿真軟件系統(tǒng)結(jié)合起來,硬件部分能夠準(zhǔn)確反映系統(tǒng)的機械特性及液壓特性,仿真軟件能夠模擬各種路況,通過這種方式驗證算法的有效性和可靠性.如圖5所示為項目組設(shè)計的硬件在環(huán)仿真平臺和ESP控制器.

    如圖6所示為硬件在環(huán)平臺的結(jié)構(gòu)原理圖.平臺采用了dSpace公司的Simulator仿真系統(tǒng)及MicroAutobox快速原型控制器.通過計算機將Carsim整車模型下載并連接到Simulator中,再將本文基于Matlab/simulink的ESP控制策略下載至MicroAutobox.

    圖6 硬件在環(huán)平臺結(jié)構(gòu)原理圖Fig.6 Architecture of HIL stucture

    駕駛員通過方向盤、驅(qū)動踏板將方向盤轉(zhuǎn)角信號和踏板位置信號輸入Simulator中.同時Simulator根據(jù)實時工況,通過CAN總線和IO口將縱向加速度、側(cè)向加速度、橫擺角速度、輪速信號輸出至ESP控制器.MicroAutobox作為ESP控制器,接收Simulator和駕駛員輸入的車輛狀態(tài)信號和傳感器信號.MicroAutobox中的ESP控制策略根據(jù)車輛狀態(tài)及路況進行實時計算,將控制指令以PWM信號的方式發(fā)送給液壓調(diào)節(jié)單元,通過液壓調(diào)節(jié)單元中的電磁閥和液壓泵驅(qū)動制動執(zhí)行機構(gòu),實現(xiàn)車輛的主動安全控制,最后通過輪缸壓力傳感器將壓力信號送至ESP控制器以實現(xiàn)閉環(huán)控制.其中ESP控制器部分可以采用MicroAutobox,也可將Simulink的ESP控制算法生成C代碼下載至32位的單片機,將單片機作為ESP控制器.

    本文在單移線和雙移線工況上進行仿真實驗以驗證穩(wěn)定性控制策略的效果.仿真車型數(shù)據(jù)如表3所示.仿真條件設(shè)置如下:車速110 km/h,在附著系數(shù)為0.35和0.85的工況下只做轉(zhuǎn)向操作,車輛沒有加速或者制動的信號輸入,進行3種工況的仿真,分別是關(guān)閉ESP功能、基于PID算法的ESP控制器和以本文提出的算法作為控制器.

    表3 某乘用車技術(shù)參數(shù)

    如圖7所示在單移線工況下的仿真結(jié)果.圖中t表示時間,p表示輪缸壓力.如圖7(a)所示為側(cè)向加速度實驗結(jié)果,在關(guān)閉ESP功能后車輛在2.5 s時就表現(xiàn)出了失穩(wěn)狀態(tài);采用自抗擾算法控制的車輛的側(cè)向加速度的變化趨勢緊密跟隨方向盤轉(zhuǎn)角變化趨勢;PID算法控制下的車輛的側(cè)向加速度未能及時跟隨方向盤轉(zhuǎn)角變化趨勢,且在回正的過程中波動較大.如圖7(b)和圖7(c)所示分別為橫擺角速度和質(zhì)心側(cè)偏角的仿真結(jié)果.從圖中可以看出,采用本文算法控制的車輛的橫擺角速度和質(zhì)心側(cè)偏角的值緊密跟隨期望值,而采用PID算法控制的車輛對目標(biāo)值的跟隨性能較差.如圖7(d)所示為本文算法控制下的車輛4個車輪的輪缸壓力,本文提出的下層算法通過調(diào)整被控車輪的制動壓力實現(xiàn)目標(biāo)滑移率的控制.

    圖7 單移線工況下硬件在環(huán)實驗數(shù)據(jù)Fig.7 Simulation results of HIL during single lane

    如圖8所示為雙移線工況下的實驗結(jié)果.該工況下的仿真條件設(shè)置為時速110 km/h,路面附著系數(shù)為0.85的高附著路面,車輛在仿真過程中沒有加速、制動的信號輸入.如圖8(a)所示為雙移線工況下的側(cè)向加速度結(jié)果,關(guān)閉ESP功能的車輛在4 s時出現(xiàn)失穩(wěn)狀態(tài).由圖可見,本文提出的自抗擾控制的車輛的側(cè)向加速度跟隨方向盤轉(zhuǎn)角變化趨勢的性能優(yōu)于PID算法.如圖8(b)和(c)所示為橫擺角速度和質(zhì)心側(cè)偏角的實驗數(shù)據(jù),本文提出的算法緊密跟隨期望值,且將質(zhì)心側(cè)偏角的值控制在10°以內(nèi)(根據(jù)文獻[9],在附著系數(shù)為0.9的路面,質(zhì)心側(cè)偏角的上限值為10°),保證了車輛穩(wěn)定性,控制效果要優(yōu)于PID控制.

    圖8 雙移線工況下硬件在環(huán)實驗數(shù)據(jù)Fig.8 Simulation results of HIL during double lane

    6 結(jié) 語

    本文針對傳統(tǒng)乘用車提出了基于自抗擾控制的ESP控制策略.ESP控制策略采用分層控制的方法,上層進行相關(guān)參數(shù)的估計及控制量期望值的計算,當(dāng)質(zhì)心側(cè)偏角和橫擺角速度兩者中任一超過門限值后,采用自抗擾控制算法計算直接橫擺力矩的計算.下層負責(zé)直接橫擺力矩的分配,依靠直接橫擺力矩和目標(biāo)滑移率之間的數(shù)學(xué)關(guān)系進行被控車輪的選擇,再進行滑移率控制,從而提高直接橫擺力矩的控制精度.

    基于dSpace硬件在環(huán)實驗平臺的單移線和雙移線實驗表明:與PID控制算法相比,本文所提控制策略顯著增強了汽車在復(fù)雜工況下的操穩(wěn)性能;該方法能夠有效準(zhǔn)確估計質(zhì)心側(cè)偏角,且準(zhǔn)確地跟蹤控制變量的理想值.下一步研究將針對車輛狀態(tài)參數(shù)準(zhǔn)確估計,從而最終實現(xiàn)ESP產(chǎn)品化.

    [1] 林程,彭春雷,曹萬科.獨立驅(qū)動電動汽車穩(wěn)定性的滑模變結(jié)構(gòu)控制[J].汽車工程,2015, 37(2): 132-138. LIN Cheng, PENG Chun-lei, CAO Wan-ke. Sliding mode variable structure control for the stability of independent drive electric vehicle [J]. Automotive Engineering, 2015, 37(2): 132-138.

    [2] 許建,張政,李翔,等.獨立驅(qū)動電動汽車橫擺力矩的模糊控制算法[J].西安交通大學(xué)學(xué)報,2014,48(7): 83-91. XU Jian, ZHANG Zheng, LI Xiang, et al. A fuzzy control system for the direct yaw moment of 4WD electric vehicles [J]. Journal of Xi’an JiaoTong University, 2014, 48(7): 83-91.

    [3] 鄒廣才,羅禹貢,李克強.基于全輪縱向力優(yōu)化分配的4WD車輛直接橫擺力矩控制[J].農(nóng)業(yè)機械學(xué)報, 2009, 40(5): 1-9. ZHOU Guang-cai, LUO Yu-gong, LI Ke-qiang. 4WD vehicle DYC based on tire longitudinal forces optimization distribution [J].Transactions of The Chinese Society for Agricultural Machinery, 2009, 40(5): 1-9.

    [4] 玄圣夷,宋傳學(xué),靳立強,等.基于多級魯棒PID控制的汽車穩(wěn)定性控制策略[J].吉林大學(xué)學(xué)報:工學(xué)版,2010,40(1): 13-18. XUAN Sheng-yi, SONG Chuan-xue, JIN Li-qiang, et al. Multilevel robust PID control strategy for vehicle stabilitycontrol [J]. Journal of Science Jilin University: Engineering Science, 2010, 40(1): 13-18.

    [5] YANG P. Stability control strategy design and experiment of distributed electric drive vehicle [J]. Journal of Mechanical Engineering, 2013, 49(24): 128-136.

    [6] 付皓.汽車電子穩(wěn)定系系統(tǒng)質(zhì)心側(cè)偏角估計與控制策略[D].吉林:吉林大學(xué),2008: 76. FU Hao. Research on sideslip angle estimation and control algortihm for vehicle electronic stability system [D]. Jilin: Jilin University, 2008: 76.

    [7] LI L, JIA G, CHEN J, et al. A novel vehicle dynamics stability control algorithm based on the hierarchical strategy with constrain of nonlinear tyre forces [J]. Vehicle System Dynamics, 2015, 53(8): 1093-1116.

    [8] CHEN Y, HEDRICK J K, GUO K A novel direct yaw moment controller for in-wheel motor electric vehicles [J]. Vehicle System Dynamics, 2013, 51(6): 925-942.

    [9] GOODARZI A, MOHAMMADI M. Stability enhancement and fuel economy of the 4-wheel-drive hybrid electric vehicles by optimal tyre force distribution [J]. Vehicle System Dynamics, 2013, 52(4): 539-561.

    [10] SHUAI Z, ZHANG H, WANG J, et al. Lateral motion control for four-wheel-independent-drive electric vehicles using optimal torque allocation and dynamic message priority scheduling [J]. Control Engineering Practice, 2014, 24(1): 55-66.[11] GUO K, LU D, CHEN S K, et al. The UniTire model: a nonlinear and non-steady-state tyre model for vehicle dynamics simulation [J]. Vehicle System Dynamics, 2005, 43(1): 341-358.

    [12] SONG P, TOMIZUKA M, ZONG C. A novel integrated chassis controller for full drive-by-wire vehicles [J]. Vehicle System Dynamics, 2015, 53(2): 215-236.

    [13] LI L. Progress on vehicle dynamics stability control system [J]. Journal of Mechanical Engineering, 2013, 49(24): 95-110.

    [14] YAN C, WANG J. Energy-efficient control allocation with applications on planar motion control of electric ground vehicles[C] ∥ American Control Conference 2011. San Francisco: [s. n.], 2011: 2719-2724.

    [15] HE P, HORI Y, KAMACHI M,et al. Future motion control to be realized by in-wheel motored electri vehicle[C] ∥ 31st Annual Conference of IEEE. Raleigh: [s. n.], 2005: 263-283.

    [16] FUKADA Y. Slip-angle estimation for vehicle stability control [J]. Vehicle System Dynamics, 1999, 32(4/5): 375-388.

    [17] 韓京清.擴張狀態(tài)觀測器參數(shù)與斐波那契數(shù)列[J].控制工程,2008,15(2): 1-3. HAN J Q. Parameters of the extended state observer and Fibonacci sequence [J]. Control Engineering of China, 2008,15(2): 1-3.

    [18] BUCCHI F, FRENDO F. A new formulation of the understeer coefficient to relate yaw torque and vehicle handing [J]. Vehicle System Dynamics, 2016, 54(6): 1-22.

    Vehicle stability control system based on active disturbance-rejection control strategy

    LIU Gang1,2, JIN Li-qiang1, WANG Yi3

    (1.StateKeyLaboratoryofAutomotiveSimulationandControl,JilinUniversity,Changchun130025,China;2.DepartmentofAutomaticControl,HenanInstituteofTechnology,Xinxiang453000,China;3.YanshanUniversity,QinhuangIsland066004,China)

    The steering stability control strategy of traditional passager vehicle was investigated and an electronic stability control system based on hierarchical control strategy was designed to improve safety performance during vehicle driving. The control strategy insisted of reference model, UniTire model, parameter estimation model, desired parameter calculation model, active disturbance rejection control (ADRC) model, direct yaw moment distrubution and slip ratio regulation model. First of all, side slip angle was calculated using dynamics calculation method. Secondly, take the evaluated side slip angle and the yaw rate obtained by sensors as control variables. When the control variables exceeded threshold value, ADRC was adopted to calculate the yaw moment. In order to determine the target slip ratio, the relation between active yaw moment and wheel slip ratio was built up based on the UniTire model; the variable PID parameters were used to control slip ratio. Finally, hardware in loop simulation (HILS) was taken under single line and double line conditions. Results show that the proposed stability controller strategy can improve the vehicle’s controllability and stability, making the control variable track the expectation closely.

    active disturbance rejection control (ADRC); vehicle dynamics; electronic stability control system; yaw moment control; hardware in loop simulation (HILS)

    2015-12-08.

    國家自然科學(xué)基金資助項目(50907030);中國博士后科學(xué)基金資助項目(2013M54024B).

    劉剛(1981—).男,博士生.從事汽車底盤電子控制技術(shù)研究, ORCID: 0000-0001-8849-169X.E-mail: gliu14@mails.jlu.edu.cn 通信聯(lián)系人:靳立強(1976—),男,教授.ORCID: 0000-0002-0545-3872.E-mail: jinlq@jlu.edu.cn

    10.3785/j.issn.1008-973X.2016.12.007

    U 461.2

    A

    1008-973X(2016)12-2289-08

    猜你喜歡
    偏角角速度質(zhì)心
    重型半掛汽車質(zhì)量與質(zhì)心位置估計
    基于GNSS測量的天宮二號質(zhì)心確定
    翼吊長涵道發(fā)動機短艙內(nèi)偏角優(yōu)化和機理研究
    2018全國Ⅱ卷選修3-4中偏角的解法探討
    歐姆表偶然誤差分析
    圓周運動角速度測量方法賞析
    半捷聯(lián)雷達導(dǎo)引頭視線角速度提取
    基于構(gòu)架點頭角速度的軌道垂向長波不平順在線檢測
    一種海洋測高衛(wèi)星質(zhì)心在軌估計算法
    航天器工程(2014年5期)2014-03-11 16:35:53
    鐵路既有曲線整正計算方法研究
    鐵道建筑(2010年11期)2010-07-30 02:04:04
    日日爽夜夜爽网站| 长腿黑丝高跟| 啪啪无遮挡十八禁网站| 夜夜夜夜夜久久久久| 欧美久久黑人一区二区| 一级黄色大片毛片| 一a级毛片在线观看| 国产精品自产拍在线观看55亚洲| 久久久久国内视频| 精品午夜福利视频在线观看一区| 波多野结衣巨乳人妻| 身体一侧抽搐| 国产精品免费一区二区三区在线| 最近最新免费中文字幕在线| 欧洲精品卡2卡3卡4卡5卡区| 少妇 在线观看| 亚洲国产毛片av蜜桃av| 一夜夜www| 亚洲成a人片在线一区二区| 黄色毛片三级朝国网站| 亚洲av成人不卡在线观看播放网| 长腿黑丝高跟| 热99re8久久精品国产| 啦啦啦韩国在线观看视频| 黄网站色视频无遮挡免费观看| 久久国产乱子伦精品免费另类| 老汉色av国产亚洲站长工具| 99在线视频只有这里精品首页| 国产精品 欧美亚洲| 曰老女人黄片| 国产欧美日韩一区二区三| 人人妻人人澡欧美一区二区 | 嫁个100分男人电影在线观看| 久久久国产成人精品二区| 亚洲国产欧美网| 999久久久国产精品视频| 十八禁网站免费在线| 男人的好看免费观看在线视频 | 中亚洲国语对白在线视频| 午夜福利视频1000在线观看 | 91精品三级在线观看| 亚洲五月婷婷丁香| 精品国产乱子伦一区二区三区| 变态另类成人亚洲欧美熟女 | 九色国产91popny在线| 国产亚洲精品久久久久久毛片| 美女午夜性视频免费| 制服丝袜大香蕉在线| e午夜精品久久久久久久| 丝袜美足系列| 91成人精品电影| 日韩视频一区二区在线观看| 99久久综合精品五月天人人| 久久久久国产精品人妻aⅴ院| 亚洲少妇的诱惑av| 满18在线观看网站| 久久久久精品国产欧美久久久| 女性被躁到高潮视频| 精品久久久久久,| 久久中文字幕人妻熟女| 亚洲免费av在线视频| 精品久久久精品久久久| 午夜a级毛片| 男女下面插进去视频免费观看| 亚洲一区高清亚洲精品| 亚洲一区中文字幕在线| 午夜免费观看网址| a在线观看视频网站| 啪啪无遮挡十八禁网站| 午夜精品在线福利| 亚洲色图 男人天堂 中文字幕| 午夜福利视频1000在线观看 | 免费高清视频大片| 午夜精品在线福利| 在线观看舔阴道视频| 亚洲精品一区av在线观看| 国产精品av久久久久免费| 国产区一区二久久| 国产精品国产高清国产av| 男人舔女人下体高潮全视频| av网站免费在线观看视频| 亚洲avbb在线观看| 国产午夜福利久久久久久| 老汉色av国产亚洲站长工具| 欧美黄色淫秽网站| 国产精品国产高清国产av| 久久亚洲精品不卡| 波多野结衣一区麻豆| 夜夜躁狠狠躁天天躁| 亚洲国产欧美日韩在线播放| 中文字幕另类日韩欧美亚洲嫩草| 老司机福利观看| 欧美激情久久久久久爽电影 | 久久国产精品人妻蜜桃| 国产精品综合久久久久久久免费 | 日本三级黄在线观看| 在线观看66精品国产| 国产av一区二区精品久久| 国产精品久久视频播放| 大码成人一级视频| 不卡av一区二区三区| 精品人妻1区二区| 美女午夜性视频免费| 日本 欧美在线| 成人18禁在线播放| 久久精品91蜜桃| 老汉色av国产亚洲站长工具| 中文字幕色久视频| 亚洲三区欧美一区| 俄罗斯特黄特色一大片| 成人18禁在线播放| 自线自在国产av| 黄色视频,在线免费观看| 首页视频小说图片口味搜索| 少妇 在线观看| 国产精品二区激情视频| 欧美黄色淫秽网站| 侵犯人妻中文字幕一二三四区| 久久亚洲精品不卡| 18禁黄网站禁片午夜丰满| 午夜久久久久精精品| 在线十欧美十亚洲十日本专区| 精品少妇一区二区三区视频日本电影| 身体一侧抽搐| 国产不卡一卡二| 日韩三级视频一区二区三区| 高清黄色对白视频在线免费看| 久久精品91无色码中文字幕| 精品国产一区二区久久| 日本在线视频免费播放| 看黄色毛片网站| 亚洲中文字幕日韩| 天天躁狠狠躁夜夜躁狠狠躁| 99热只有精品国产| 国产一区二区三区在线臀色熟女| 国产精品久久久久久人妻精品电影| 国产成人av教育| 欧美乱色亚洲激情| 深夜精品福利| av片东京热男人的天堂| 国产精品亚洲一级av第二区| 免费高清在线观看日韩| 久久影院123| 成人av一区二区三区在线看| 欧美日韩亚洲综合一区二区三区_| 久久久久久亚洲精品国产蜜桃av| 国产精品免费视频内射| 亚洲成av片中文字幕在线观看| 国产欧美日韩一区二区精品| 亚洲国产精品合色在线| 国产一区二区激情短视频| 欧美黄色淫秽网站| 麻豆久久精品国产亚洲av| 黑丝袜美女国产一区| 亚洲第一欧美日韩一区二区三区| 日本黄色视频三级网站网址| 日本a在线网址| 黄色丝袜av网址大全| 国产99白浆流出| 成在线人永久免费视频| 欧洲精品卡2卡3卡4卡5卡区| 国产伦一二天堂av在线观看| 国产欧美日韩一区二区三| 操出白浆在线播放| 国语自产精品视频在线第100页| 久久精品国产亚洲av香蕉五月| 1024视频免费在线观看| 黑人欧美特级aaaaaa片| 亚洲第一av免费看| 亚洲一卡2卡3卡4卡5卡精品中文| 亚洲在线自拍视频| 一级a爱片免费观看的视频| 电影成人av| 免费在线观看日本一区| 丁香欧美五月| 一本久久中文字幕| 欧美 亚洲 国产 日韩一| 国产免费av片在线观看野外av| 亚洲五月婷婷丁香| 岛国在线观看网站| 久久久久久国产a免费观看| 精品午夜福利视频在线观看一区| 国产午夜福利久久久久久| 精品久久久久久久毛片微露脸| 免费不卡黄色视频| 国产精品免费视频内射| 亚洲九九香蕉| 怎么达到女性高潮| 国产高清激情床上av| 亚洲一卡2卡3卡4卡5卡精品中文| 少妇熟女aⅴ在线视频| 岛国视频午夜一区免费看| 国内精品久久久久久久电影| 成人特级黄色片久久久久久久| 国产精品一区二区三区四区久久 | 午夜福利成人在线免费观看| 老熟妇仑乱视频hdxx| 51午夜福利影视在线观看| 日韩欧美国产在线观看| 色在线成人网| 精品欧美国产一区二区三| 欧美乱色亚洲激情| 久久久久久国产a免费观看| 精品久久久久久成人av| 波多野结衣高清无吗| 一个人观看的视频www高清免费观看 | 欧美不卡视频在线免费观看 | 日韩 欧美 亚洲 中文字幕| 91国产中文字幕| 色播在线永久视频| 国产欧美日韩精品亚洲av| www国产在线视频色| 黄网站色视频无遮挡免费观看| 精品国产美女av久久久久小说| 男人操女人黄网站| 欧美激情久久久久久爽电影 | 韩国av一区二区三区四区| 黄色女人牲交| 国产成人欧美在线观看| 老熟妇仑乱视频hdxx| 久久亚洲真实| 欧美日本视频| 国产一区二区三区综合在线观看| 成年女人毛片免费观看观看9| 侵犯人妻中文字幕一二三四区| 国产一区二区在线av高清观看| 亚洲成a人片在线一区二区| 亚洲熟女毛片儿| 中文字幕另类日韩欧美亚洲嫩草| 国产在线精品亚洲第一网站| 国产一区在线观看成人免费| 国产精品秋霞免费鲁丝片| 国产乱人伦免费视频| 国产亚洲欧美在线一区二区| 亚洲情色 制服丝袜| 涩涩av久久男人的天堂| 无人区码免费观看不卡| 精品国产超薄肉色丝袜足j| 欧美日韩乱码在线| 亚洲av成人一区二区三| 变态另类丝袜制服| 久久热在线av| 色综合亚洲欧美另类图片| 两人在一起打扑克的视频| 麻豆成人av在线观看| 国内毛片毛片毛片毛片毛片| 91精品国产国语对白视频| 久久精品aⅴ一区二区三区四区| 午夜精品久久久久久毛片777| 热99re8久久精品国产| 在线观看66精品国产| 亚洲狠狠婷婷综合久久图片| 91大片在线观看| 日韩欧美国产一区二区入口| 国产熟女午夜一区二区三区| 麻豆成人av在线观看| 少妇裸体淫交视频免费看高清 | 亚洲精品粉嫩美女一区| 美女扒开内裤让男人捅视频| 亚洲国产精品sss在线观看| 在线观看免费午夜福利视频| 亚洲三区欧美一区| 免费看十八禁软件| 777久久人妻少妇嫩草av网站| 在线天堂中文资源库| 日韩欧美在线二视频| 男女午夜视频在线观看| 久久这里只有精品19| 欧美成人午夜精品| 亚洲五月天丁香| 女人爽到高潮嗷嗷叫在线视频| 日本 av在线| 三级毛片av免费| 一边摸一边抽搐一进一出视频| 日韩精品免费视频一区二区三区| 国产精品av久久久久免费| 天堂影院成人在线观看| 激情视频va一区二区三区| 成人三级黄色视频| 亚洲专区国产一区二区| 国产精品一区二区三区四区久久 | 男男h啪啪无遮挡| 久久狼人影院| 精品熟女少妇八av免费久了| 女性被躁到高潮视频| 久久中文看片网| 桃色一区二区三区在线观看| 精品久久久久久久久久免费视频| 老司机深夜福利视频在线观看| 老熟妇仑乱视频hdxx| 成在线人永久免费视频| 精品高清国产在线一区| 亚洲成av片中文字幕在线观看| 中文字幕色久视频| 最好的美女福利视频网| 午夜精品久久久久久毛片777| 欧美精品啪啪一区二区三区| 大陆偷拍与自拍| 精品一区二区三区av网在线观看| 色综合亚洲欧美另类图片| 婷婷丁香在线五月| 国产xxxxx性猛交| 老司机午夜十八禁免费视频| 国产欧美日韩一区二区三区在线| 青草久久国产| 国产精品1区2区在线观看.| 日本免费a在线| 一区二区三区国产精品乱码| videosex国产| 亚洲精品av麻豆狂野| 18禁美女被吸乳视频| 50天的宝宝边吃奶边哭怎么回事| 国产精品久久久久久亚洲av鲁大| 国产精品自产拍在线观看55亚洲| 国产在线精品亚洲第一网站| 国产91精品成人一区二区三区| 少妇的丰满在线观看| 91字幕亚洲| 香蕉国产在线看| 国产精品免费视频内射| 久久精品91无色码中文字幕| 免费无遮挡裸体视频| 免费观看人在逋| 国产精品电影一区二区三区| 丁香欧美五月| 午夜老司机福利片| 麻豆久久精品国产亚洲av| 亚洲人成电影观看| 在线观看免费视频日本深夜| 欧美成人一区二区免费高清观看 | 嫩草影院精品99| 中文字幕av电影在线播放| 亚洲久久久国产精品| 国产人伦9x9x在线观看| 久久精品国产清高在天天线| 色尼玛亚洲综合影院| 啦啦啦免费观看视频1| 成人国语在线视频| 精品久久久久久,| e午夜精品久久久久久久| 久久久久亚洲av毛片大全| 在线观看日韩欧美| 老熟妇仑乱视频hdxx| 一区在线观看完整版| 亚洲国产高清在线一区二区三 | 国产亚洲av嫩草精品影院| 亚洲欧美日韩另类电影网站| 18美女黄网站色大片免费观看| 精品卡一卡二卡四卡免费| 色综合亚洲欧美另类图片| 免费在线观看亚洲国产| www日本在线高清视频| 欧美老熟妇乱子伦牲交| 久久久精品国产亚洲av高清涩受| 午夜日韩欧美国产| 搡老岳熟女国产| 国产精品 国内视频| 亚洲情色 制服丝袜| 美女大奶头视频| 欧美成人一区二区免费高清观看 | 在线观看免费视频日本深夜| 久久久水蜜桃国产精品网| 99香蕉大伊视频| 丁香六月欧美| 欧美av亚洲av综合av国产av| 亚洲人成电影观看| 午夜福利影视在线免费观看| 亚洲片人在线观看| 制服诱惑二区| 欧美乱色亚洲激情| 宅男免费午夜| 国产欧美日韩一区二区精品| 欧美性长视频在线观看| 美女高潮喷水抽搐中文字幕| 国产精品亚洲av一区麻豆| 免费无遮挡裸体视频| 日日干狠狠操夜夜爽| 国产成人免费无遮挡视频| 男女床上黄色一级片免费看| 日日爽夜夜爽网站| 天天添夜夜摸| 精品一区二区三区av网在线观看| 97人妻精品一区二区三区麻豆 | 男人的好看免费观看在线视频 | 一区二区三区国产精品乱码| 国产av精品麻豆| 看片在线看免费视频| av有码第一页| 好看av亚洲va欧美ⅴa在| 中文字幕另类日韩欧美亚洲嫩草| 亚洲第一电影网av| 欧美黄色片欧美黄色片| 日日爽夜夜爽网站| 亚洲性夜色夜夜综合| 国产成人精品久久二区二区免费| 国产精品,欧美在线| 九色国产91popny在线| 99精品欧美一区二区三区四区| 性欧美人与动物交配| 人人妻人人澡欧美一区二区 | 国产不卡一卡二| 法律面前人人平等表现在哪些方面| 亚洲国产毛片av蜜桃av| 精品久久久久久成人av| 亚洲第一青青草原| 啦啦啦观看免费观看视频高清 | 咕卡用的链子| 狠狠狠狠99中文字幕| 久久久久国产精品人妻aⅴ院| 日日干狠狠操夜夜爽| 老司机在亚洲福利影院| 国产成人影院久久av| 在线观看免费视频日本深夜| 久久久国产成人免费| 欧美+亚洲+日韩+国产| 亚洲成a人片在线一区二区| 亚洲国产毛片av蜜桃av| 欧洲精品卡2卡3卡4卡5卡区| 99热只有精品国产| 一级毛片高清免费大全| 性色av乱码一区二区三区2| 一边摸一边抽搐一进一出视频| 国产成人精品久久二区二区91| 国产精品乱码一区二三区的特点 | 中文亚洲av片在线观看爽| tocl精华| 国产精品影院久久| 精品久久蜜臀av无| 黄色女人牲交| 看免费av毛片| av在线天堂中文字幕| 女性被躁到高潮视频| 一级,二级,三级黄色视频| 激情视频va一区二区三区| 在线观看日韩欧美| 久久国产乱子伦精品免费另类| 狂野欧美激情性xxxx| 国产主播在线观看一区二区| 淫妇啪啪啪对白视频| 两性午夜刺激爽爽歪歪视频在线观看 | 在线观看66精品国产| 欧美在线一区亚洲| 熟妇人妻久久中文字幕3abv| 亚洲av美国av| 老汉色av国产亚洲站长工具| 级片在线观看| 国产又色又爽无遮挡免费看| 99在线视频只有这里精品首页| 波多野结衣巨乳人妻| 午夜福利,免费看| av片东京热男人的天堂| 一进一出抽搐gif免费好疼| 一本久久中文字幕| 老司机在亚洲福利影院| 国产高清有码在线观看视频 | 一级黄色大片毛片| 久久中文字幕一级| 中文字幕久久专区| 19禁男女啪啪无遮挡网站| 国产成人免费无遮挡视频| 国产精品久久久久久精品电影 | 首页视频小说图片口味搜索| 国产高清视频在线播放一区| 亚洲免费av在线视频| 在线观看66精品国产| 丝袜美腿诱惑在线| 99国产精品99久久久久| 俄罗斯特黄特色一大片| 国产精品免费视频内射| 免费在线观看黄色视频的| 91成人精品电影| 桃色一区二区三区在线观看| 身体一侧抽搐| 一级片免费观看大全| 9191精品国产免费久久| 天堂动漫精品| 国产精品野战在线观看| 在线天堂中文资源库| 国产亚洲欧美精品永久| 亚洲精品在线观看二区| 欧美人与性动交α欧美精品济南到| 黄片小视频在线播放| 脱女人内裤的视频| 久久中文字幕一级| 黑人操中国人逼视频| 久久精品国产清高在天天线| 1024香蕉在线观看| 久久久精品国产亚洲av高清涩受| 欧美日本视频| 久久国产精品男人的天堂亚洲| 给我免费播放毛片高清在线观看| 久久青草综合色| 日本 av在线| 亚洲男人的天堂狠狠| 精品人妻在线不人妻| 熟女少妇亚洲综合色aaa.| 纯流量卡能插随身wifi吗| 男人舔女人的私密视频| av超薄肉色丝袜交足视频| 怎么达到女性高潮| 免费在线观看完整版高清| 国内久久婷婷六月综合欲色啪| 久久精品91蜜桃| 午夜福利在线观看吧| 久久精品亚洲熟妇少妇任你| 极品教师在线免费播放| 91九色精品人成在线观看| 久久久国产精品麻豆| 精品国内亚洲2022精品成人| 免费高清视频大片| 久久久久精品国产欧美久久久| 好男人电影高清在线观看| 人人妻人人澡人人看| 波多野结衣高清无吗| 欧美色欧美亚洲另类二区 | 黄色视频不卡| 免费在线观看日本一区| 天堂√8在线中文| 脱女人内裤的视频| 国产三级黄色录像| 色综合婷婷激情| av在线天堂中文字幕| 久久精品亚洲精品国产色婷小说| 欧美日韩中文字幕国产精品一区二区三区 | 亚洲成国产人片在线观看| 日韩精品免费视频一区二区三区| 淫妇啪啪啪对白视频| 欧美日韩乱码在线| 亚洲精品国产区一区二| 精品午夜福利视频在线观看一区| 久久天堂一区二区三区四区| 欧美乱码精品一区二区三区| 丝袜人妻中文字幕| 久久久久亚洲av毛片大全| 在线播放国产精品三级| 亚洲av电影在线进入| 亚洲专区国产一区二区| 国内精品久久久久久久电影| 亚洲激情在线av| 91九色精品人成在线观看| 亚洲人成电影观看| 精品久久久久久,| 淫妇啪啪啪对白视频| 精品免费久久久久久久清纯| 老司机靠b影院| 国产亚洲欧美精品永久| 亚洲中文av在线| 免费无遮挡裸体视频| 黄色片一级片一级黄色片| 久久婷婷人人爽人人干人人爱 | 亚洲人成77777在线视频| 一二三四在线观看免费中文在| 国产欧美日韩一区二区精品| 日韩 欧美 亚洲 中文字幕| 成人欧美大片| av超薄肉色丝袜交足视频| 国产精品久久久久久亚洲av鲁大| 国产成人啪精品午夜网站| 久久久久久国产a免费观看| 国产精品免费视频内射| av片东京热男人的天堂| 岛国视频午夜一区免费看| 亚洲va日本ⅴa欧美va伊人久久| 成人手机av| 亚洲第一av免费看| 一级毛片精品| 视频在线观看一区二区三区| 女人高潮潮喷娇喘18禁视频| 黑人欧美特级aaaaaa片| av视频在线观看入口| 日韩免费av在线播放| 亚洲国产精品久久男人天堂| 久久久精品欧美日韩精品| 亚洲国产高清在线一区二区三 | 亚洲精品国产色婷婷电影| 天天添夜夜摸| 国产一区二区三区视频了| 99久久综合精品五月天人人| 亚洲中文av在线| 亚洲欧美日韩另类电影网站| 一区在线观看完整版| 欧美黑人精品巨大| 国产熟女xx| 亚洲成av人片免费观看| 啦啦啦免费观看视频1| 麻豆一二三区av精品| 高清毛片免费观看视频网站| 97人妻精品一区二区三区麻豆 | 香蕉久久夜色| 成年人黄色毛片网站| avwww免费| 丝袜美足系列| 青草久久国产| 在线观看免费视频网站a站| 亚洲午夜理论影院| 制服诱惑二区| 男女做爰动态图高潮gif福利片 | 天堂√8在线中文| 国产精品98久久久久久宅男小说| 欧美日韩一级在线毛片| 他把我摸到了高潮在线观看| 两性午夜刺激爽爽歪歪视频在线观看 | 黄色女人牲交| 无限看片的www在线观看| 操出白浆在线播放| 欧美老熟妇乱子伦牲交| 又黄又爽又免费观看的视频| 很黄的视频免费| 美女高潮喷水抽搐中文字幕| 老熟妇乱子伦视频在线观看| 免费在线观看亚洲国产| 深夜精品福利| 国语自产精品视频在线第100页| 香蕉国产在线看| 婷婷六月久久综合丁香|