李旭 史曉華 李瑞川 王建春 馬勇,2
(1.山東科技大學,青島 266590;2.山東海卓電液控制工程技術研究院,日照 276800)
主題詞:主動前輪轉(zhuǎn)向 理想變傳動比 附加轉(zhuǎn)角 滑??刂?轉(zhuǎn)向性能
主動前輪轉(zhuǎn)向系統(tǒng)是介于電動助力轉(zhuǎn)向和線控轉(zhuǎn)向之間的轉(zhuǎn)向系統(tǒng)[1],可在一定范圍內(nèi)實現(xiàn)變傳動比控制,對于提高車輛的操縱性、穩(wěn)定性和軌跡保持能力具有重要的意義[2-3]。
目前,針對主動前輪轉(zhuǎn)向的研究主要集中在可變傳動比設計和車輛穩(wěn)定性控制方面。針對可變傳動比:Tajima J首先提出車輛的理想變傳動比與穩(wěn)態(tài)橫擺角速度增益有關,為保證車輛具有良好的轉(zhuǎn)向性能,穩(wěn)態(tài)橫擺角速度增益應該是不隨車速變化的定值[4];商高高設計了基于穩(wěn)態(tài)橫擺角速度增益的理想變傳動比曲線[5];周兵通過對比分析5種擬合變傳動比曲線,證明改進型S函數(shù)與理想變傳動比曲線具有最高的接近程度[6]。以上研究僅探究了主動前輪轉(zhuǎn)向系統(tǒng)的傳動比應遵循的規(guī)律,尚未提出具體的控制策略。
針對車輛穩(wěn)定性控制,Nam K、桑楠等根據(jù)駕駛?cè)笋{駛意圖和車輛的行駛狀態(tài)分別利用H∞魯棒控制、自抗擾控制等求解了主動前輪轉(zhuǎn)向系統(tǒng)的附加轉(zhuǎn)角[7-13]。這些研究通常假定主動前輪轉(zhuǎn)向系統(tǒng)能夠理想地產(chǎn)生所需要的附加轉(zhuǎn)角,缺乏對該附加轉(zhuǎn)角應遵循的理想傳動比規(guī)律的研究[14-16]。
然而,基于理想變傳動比的附加轉(zhuǎn)角控制是決定主動前輪轉(zhuǎn)向系統(tǒng)性能優(yōu)劣的關鍵因素。王春燕提出一種理想傳動比規(guī)律下的主動前輪附加轉(zhuǎn)角閉環(huán)控制策略,但所采用的PID控制對于車輛行駛過程這一復雜時變、強非線性對象控制不佳。滑??刂祈憫杆佟Ⅳ敯粜院?、物理實現(xiàn)簡單,對非線性系統(tǒng)具有良好的控制效果,在響應速度、控制精度上均優(yōu)于傳統(tǒng)的PID控制。
基于此,本文提出一種基于理想變傳動比規(guī)律的主動前輪轉(zhuǎn)向滑??刂撇呗浴R灾鲃忧拜嗈D(zhuǎn)向系統(tǒng)為研究對象,建立七自由度車輛模型、Dugoff非線性輪胎模型;確定了固定橫擺角速度增益下的理想變傳動比規(guī)律。在此基礎上,提出主動前輪轉(zhuǎn)向附加轉(zhuǎn)角的滑??刂撇呗裕⒃陔A躍、蛇行轉(zhuǎn)向工況下對其控制效果進行了仿真分析。
考慮車輛縱向、側(cè)向、橫擺運動以及4個車輪的旋轉(zhuǎn)運動,建立非線性七自由度整車模型如圖1所示。
圖1 七自由度整車模型
結(jié)合整車受力情況,根據(jù)牛頓第二定律可得車輛非線性動力學方程為:
式中,m為整車質(zhì)量;vx、vy分別為縱向、橫向車速;r為橫擺角速度;Fxi、Fyi分別為輪胎縱向力、側(cè)向力,i=fl,fr,rl,rr分別表示左前輪、右前輪、左后輪、右后輪;δc為前輪轉(zhuǎn)角;Iz為整車繞z軸的轉(zhuǎn)動慣量;a、b分別為前、后軸到質(zhì)心的距離;tw1、tw2分別為前、后軸輪距;Itw為車輪的轉(zhuǎn)動慣量;ωi為各車輪的角速度;R為車輪半徑;Tbi為各車輪驅(qū)動力矩;Tdi為各車輪制動力矩。
與魔術公式模型相比,Dugoff非線性模型允許各輪胎有獨立的側(cè)向剛度和縱向剛度,所需輪胎模型參數(shù)較少,可以較精確地擬合輪胎的穩(wěn)態(tài)運動過程[17-19],因此本文選用Dugoff輪胎模型。
設4個輪胎具有相同的縱向剛度,同軸的左、右側(cè)輪胎具有相同的側(cè)偏剛度,縱向力、側(cè)向力計算公式為:
式中,cs為輪胎的縱向剛度;si為各輪胎的滑移率;ki為各輪胎的側(cè)偏剛度;αi為各輪胎的側(cè)偏角;Fzi為各輪胎的垂向載荷;μ為路面附著系數(shù);λi為輪胎動態(tài)參數(shù)。
各輪胎垂向載荷、滑移率、側(cè)偏角的計算公式,可由七自由度整車模型求得。
線性二自由度汽車微分方程能夠反映汽車曲線運動的最基本特征。等速行駛時,前輪角階躍輸入即可使車輛達到穩(wěn)態(tài)即等速圓周行駛。此時=0,m=0,rm為穩(wěn)態(tài)橫擺角速度。線性二自由度微分方程為:
式中,k1、k2分別為前、后輪輪胎的側(cè)偏剛度。
穩(wěn)態(tài)橫擺角速度增益為:
Tajima J提出,為保證轉(zhuǎn)向盤轉(zhuǎn)角和汽車航向角之間呈現(xiàn)出與車速無關的固定比例關系,穩(wěn)態(tài)橫擺角速度增益應為不隨車速變化的定值,將基于固定橫擺角速度增益確定的變傳動比稱為車輛的理想變傳動比[4]。
德國汽車研究所通過試驗得出[20],轎車的穩(wěn)態(tài)橫擺角速度增益范圍為0.16~0.33 s-1,本文選取Gsw=0.275 s-1。
固定穩(wěn)態(tài)橫擺角速度增益下,傳動比為:
轉(zhuǎn)向傳動比過小,轉(zhuǎn)向會過于靈敏,駕駛員微小的誤操作會導致較大的車輛響應;反之則過于遲緩,不利于車輛的換道、避障。因此,變傳動比曲線應有上、下界,下限Imin=fI(vx1),上限Imax=fI(vx2),其中vx1為下臨界車速,vx2為上臨界車速。
根據(jù)車輛仿真模型的基本參數(shù)(見表1),可確定理想變傳動比曲線如圖2所示,其規(guī)律為:
表1 主要仿真參數(shù)
圖2 理想變傳動比曲線
滑模控制策略通過將實際橫擺角速度r與其理想值rd作對比,構建滑模面迫使實際橫擺角速度跟隨理想值的變化而變化。控制器的設計基于線性二自由度模型,理想傳動比規(guī)律下,線性二自由度模型可用矩陣表示為:
式中,g為重力加速度。
以七自由度整車模型的橫擺角速度作為車輛的實際橫擺角速度,二自由度汽車模型的橫擺角速度作為理想橫擺角速度,將兩者之差作為控制誤差,并以此構建滑模面s:
式中,λ為正加權系數(shù)。
令=0,可得前輪轉(zhuǎn)角的等效輸入量δfeq:
分別為線性二自由度模型的矩陣元素。
主動前輪轉(zhuǎn)向的滑??刂坡捎蓛刹糠纸M成。一部分為等效控制律,確保系統(tǒng)能夠到達滑模面。另一部分確保系統(tǒng)在有擾動的情況下,仍能收斂于滑模面,滑模控制結(jié)構為:
式中,δf為總輸入轉(zhuǎn)角;sgn為符號函數(shù);Kr為滑模切換增益。
為了保證系統(tǒng)閉環(huán)穩(wěn)定,Kr的選擇應滿足李雅普諾夫函數(shù)必要條件[21]:
符號函數(shù)sgn在零點附近不連續(xù),導致估計值在滑模面兩側(cè)的高速切換中產(chǎn)生高頻抖振。為消除抖振,用飽和函數(shù)代替符號函數(shù):
式中,φ為邊界層厚度。
滑模面附近邊界層如圖3所示。由圖3可以看出,φ越小,系統(tǒng)的狀態(tài)誤差越小,本文取φ=0.05。
圖3 滑模面附近邊界層
滑??刂坡蔀椋?/p>
式中,δAFS為主動前輪轉(zhuǎn)向機構提供的附加轉(zhuǎn)角。
目前的主動前輪轉(zhuǎn)向滑??刂撇呗允窍葘鲃颖菼視為一固定值,代入二自由度車輛模型求出穩(wěn)態(tài)橫擺角速度,與車輛的實際橫擺角速度進行比較,求出需要提供的附加轉(zhuǎn)角,從而達到改變傳動比的目的[22-23]。本文將這種變傳動比方式稱為普通變傳動比,用于對本文所提出的基于理想變傳動比規(guī)律的主動前輪轉(zhuǎn)向滑模控制策略進行仿真對比及分析。
在MATLAB/Simulink中搭建七自由度整車、Dugoff輪胎、固定速比轉(zhuǎn)向器、滑??刂破鞯认到y(tǒng)模型進行仿真,主要仿真參數(shù)見表1。
滑??刂平Y(jié)構如圖4所示,首先設計滑??刂撇呗赃M行理想橫擺角速度跟蹤,得到附加轉(zhuǎn)角,將疊加后的輸入轉(zhuǎn)角輸入固定速比轉(zhuǎn)向器轉(zhuǎn)化為對應的車輪轉(zhuǎn)角,將其輸入到整車模型,實現(xiàn)基于理想變傳動比的主動前輪轉(zhuǎn)向滑??刂?,Simulink控制框圖如圖5所示,其中vx0為初始車速。
圖4 主動前輪轉(zhuǎn)向滑??刂平Y(jié)構
在轉(zhuǎn)向盤轉(zhuǎn)角幅值為30°,車速分別為20 km/h、80 km/h和120 km/h條件下,利用正弦輸入模擬車輛行駛過程中的換道操作,角階躍輸入模擬緊急轉(zhuǎn)向操作,進行仿真分析,仿真結(jié)果如圖6、圖7所示。
圖5 Simulink仿真控制框圖
圖6 不同車速下車輛的正弦輸入響應
由圖6、圖7可以看出:車速為20 km/h時,理想變傳動比車輛的橫擺角速度、質(zhì)心側(cè)偏角、側(cè)向加速度與普通變傳動比車輛相比,響應幅值均明顯增大,表明此時車輛轉(zhuǎn)向更為輕便;當車速增加至80 km/h和120 km/h時,與普通變傳動比車輛相比,理想變傳動比車輛的質(zhì)心側(cè)偏角、橫擺角速度、側(cè)向加速度的響應幅值均變小,表明此時車輛具有更好的轉(zhuǎn)向穩(wěn)定性能,駕駛員的路感更強。
圖7 不同車速下車輛的角階躍輸入響應
在較高車速(120 km/h)下進行轉(zhuǎn)向操作時:普通變傳動比車輛的質(zhì)心側(cè)偏角、側(cè)向加速度響應幅值出現(xiàn)較大的波動,車輛轉(zhuǎn)向穩(wěn)定性明顯變差;基于理想變傳動比控制的車輛能夠更快地達到穩(wěn)定狀態(tài),且大幅度減小響應幅值??梢娎硐胱儌鲃颖瓤刂茖τ谔岣吒咚傩旭傑囕v的轉(zhuǎn)向穩(wěn)定性具有重要意義。
本文利用七自由度車輛模型、Dugoff非線性輪胎模型、固定速比轉(zhuǎn)向器模型建立了用于主動前輪轉(zhuǎn)向分析的系統(tǒng)模型,確定了固定橫擺角速度增益下的理想變傳動比規(guī)律,并提出基于該理想變傳動比的主動前輪轉(zhuǎn)向附加轉(zhuǎn)角的滑??刂撇呗?。研究結(jié)果表明,與普通變傳動比控制相比,基于理想變傳動比規(guī)律的主動前輪轉(zhuǎn)向滑??刂撇呗钥捎行ПWC在低速時轉(zhuǎn)向車輛具有更好的轉(zhuǎn)向靈敏度,轉(zhuǎn)向更為輕便省力,高速時大幅度減小轉(zhuǎn)向車輛的響應幅值,具有更好的操縱穩(wěn)定性。基于理想變傳動比規(guī)律的車輛具有更為理想的轉(zhuǎn)向性能。
進一步研究方向:針對高速行駛的車輛在低附著路面上轉(zhuǎn)彎時極易發(fā)生側(cè)滑的情況,研究如何利用模型預測控制理論補償預測力與參考值之間的差值實現(xiàn)跟蹤控制,進一步提高車輛高速轉(zhuǎn)彎工況下的穩(wěn)定性。