李 碩,馬 順,李艷山,王 宇,陳曉龍
(中車青島四方機車車輛股份有限公司,山東 青島266580)
高速列車的懸架系統(tǒng)對列車的平穩(wěn)運行發(fā)揮著非常重要的作用[1-3]。傳統(tǒng)的列車懸架系統(tǒng)由一系懸架組成,不存在控制系統(tǒng),平穩(wěn)性相對較差。高速列車在一系懸架的基礎(chǔ)上增加二系空氣懸掛,使得列車舒適性大大增加,空氣懸架的應(yīng)用對于提高高速列車的平順性和舒適性具有巨大的作用。
目前對于高速列車懸架系統(tǒng)的研究,主要集中在懸架系統(tǒng)控制算法以及懸架結(jié)構(gòu)和參數(shù)優(yōu)化方面。如鄧力等[4]針對傳統(tǒng)半主動控制的缺陷,提出了開閉環(huán)優(yōu)化控制算法,大大提高了列車橫向運動的穩(wěn)定性。曹青松等[5]考慮了列車在運行過程中的干擾,利用標(biāo)準(zhǔn)正交基函數(shù)進(jìn)行預(yù)測控制。閆紅衛(wèi)等[6]比較了連續(xù)型和開關(guān)型天棚阻尼控制的優(yōu)缺點,提出了適用于不同場合的控制方式。于日偉等[7]建立的耦合系統(tǒng)包括座椅、車體和轉(zhuǎn)向架,并提出懸掛參數(shù)聯(lián)合優(yōu)化方法。崔慶霞等[8]和陳志輝等[9]通過準(zhǔn)零剛度隔振方法,在二系懸掛結(jié)構(gòu)中加入一對橫向線性彈簧,使車體的垂向振動明顯減小。
以上學(xué)者在研究過程中,忽略了懸架系統(tǒng)非線性對于控制效果的影響。由于二系懸架為空氣懸架,存在著許多非線性因素,這些非線性因素在一定的載荷激勵和頻域范圍內(nèi)表現(xiàn)十分突出,對控制效果具有非常大的影響,且線性控制策略并不適合于非線性系統(tǒng)[10]。因此本文著重研究二系空氣懸架的非線性特性,并通過仿真分析了不同速度以及干擾下的車體的振動情況,為后續(xù)高速列車的平順性控制研究奠定了基礎(chǔ)。
二系懸掛的核心是空氣彈簧,而空氣彈簧具有典型的非線性特性[11],采用一般的線性控制方式會產(chǎn)生比較大的控制偏差。因此本文通過研究空氣彈簧的非線性進(jìn)而對整個懸架系統(tǒng)進(jìn)行非線性控制。本文選用前進(jìn)牌1884N 型橡膠空氣彈簧,并在文獻(xiàn)[12]的研究基礎(chǔ)上,利用其實驗得出的數(shù)據(jù),進(jìn)行空氣懸架的控制研究。本文選取高速列車2自由度非線性空氣懸架作為研究對象,并假設(shè)車輪在軌道運行過程中不脫離軌道面,其模型簡圖見圖1。

圖1 懸架系統(tǒng)模型
根據(jù)牛頓第二定律,其運動學(xué)方程可描述為

其中:

式中:ms為車體質(zhì)量;mu為轉(zhuǎn)向架質(zhì)量;Fk為二系空氣懸架剛度,k1,k2,k3為空氣彈簧的剛度系數(shù);Fc為二系空氣懸架阻尼;bc1為空氣彈簧的阻尼系數(shù);u為主動懸架的控制力;Fs為系統(tǒng)受到不確定的外部擾動,因為受到的外部擾動不可能無限大,因此該外部擾動是有界的;Ft為一系懸架剛度系數(shù);Fb為一系懸架阻尼系數(shù);zs為車體的位移;zu為轉(zhuǎn)向架的位移;z0為軌道激勵輸入。
本文所采用的路面輸入為德國軌道高低不平順功率譜密度[7],其表達(dá)式如下所示:

式中:z0為高速軌道路面不平順輸入,v為列車行駛速度,Ωc、Ωr為截斷空間頻率,Av為軌道粗糙度系數(shù)。

式中:參數(shù)θ有界,即θmin≤θ≤θmax;Fs有界,即Fs≤D;且:

針對上述狀態(tài)空間方程設(shè)計控制器,取跟蹤誤差為

式中:x1r為參考軌跡信號,當(dāng)跟蹤誤差盡可能的小,會達(dá)到控制的目標(biāo)。對上式進(jìn)行求導(dǎo)可得:

選擇x2作為誤差動態(tài)的虛擬控制輸入,其理想函數(shù)為α,則兩者之間的誤差為

則公式(5)可以重新寫為

因此,該控制器的目標(biāo)是設(shè)計一個虛擬控制率α,使得跟蹤軌跡誤差e1盡可能的小,選擇第一個Lyapunov函數(shù)為

對上式進(jìn)行求導(dǎo)可得:

選擇如下的虛擬控制率:

式中:t1為一個正常數(shù)。則式(9)可以改寫為

從上式可以看出,當(dāng)e2→0 時,那么V˙1= -t1e21≤0,就可以保證e1是漸進(jìn)趨近于零的。
對e2進(jìn)行求導(dǎo)可得:


為了消除系統(tǒng)外部擾動對系統(tǒng)的影響,采用滑模控制的等效控制和切換控制的思想,選擇的控制律為

式中:t2為正實數(shù),那么有:

取系統(tǒng)的自適應(yīng)率為

則式(15)可以改寫為

由系統(tǒng)Lyapunov定律可知,隨著t→0,V2→0,則e1與e2也將趨近于0,則系統(tǒng)可以達(dá)到預(yù)定的目標(biāo)。
為了削弱系統(tǒng)抖振帶來的危害,采用飽和函數(shù)代替符號函數(shù)的做法[13],令控制律為

為了找到系統(tǒng)的零動態(tài)響應(yīng),令e1=0=0。因此可以得到:

將上式(19)代入到中,可以得到:


從上式可以看出,此2 階系統(tǒng)的所有特征根都具有負(fù)實部,具有Hurwitz 判據(jù)特征,因此該系統(tǒng)的零動態(tài)是漸近穩(wěn)定的。
由于外部擾動信號有范圍,因此本文假設(shè)為

空氣彈簧參數(shù)選取參照前文所述文獻(xiàn)[12],選取參數(shù)為0.8 MPa 下的空氣彈簧試驗擬合系數(shù),其他系統(tǒng)參數(shù)選擇參照文獻(xiàn)[14],具體數(shù)值如表1所示。軌道激勵參數(shù)選取如表2所示,低干擾譜采用的速度為300 km/h,高干擾譜采用200 km/h。控制器參數(shù)選取如表3所示,且控制器設(shè)計過程中的參考軌跡信號設(shè)置為0,即車體的參考信號為不振動,利用MATLAB/Simulink對系統(tǒng)進(jìn)行控制仿真。

表1 懸架系統(tǒng)參數(shù)選擇

表2 軌道激勵參數(shù)選擇

表3 控制器參數(shù)選擇
低軌道路面激勵模型如圖2所示,仿真時間為10秒。
圖3所示為車體的垂直振動位移仿真結(jié)果,圖4所示為車體的垂直振動加速度仿真結(jié)果,圖5所示為車體質(zhì)量的變化曲線,圖6所示為車體的垂直振動位移功率譜密度曲線,圖7所示為車體的垂直振動加速度功率譜密度曲線仿真結(jié)果。
懸架系統(tǒng)的作用在于降低車體上乘客的振動,提高列車在運行過程中的平順性和舒適性。通過圖3以及圖4可以看出,在低干擾300 km/h的路面激勵下,綜合考慮了車體質(zhì)量的變化以及外部擾動的影響,所設(shè)計的控制器能夠非常好地降低車體的垂直振動和位移,同時圖6和圖7的頻域仿真結(jié)果驗證了時域的仿真結(jié)果,說明所設(shè)計的控制器的有效性,突出了控制器在降低車輛位移和振動的巨大作用。

圖2 路面激勵輸入

圖3 車體垂直位移響應(yīng)

圖4 車體垂直振動加速度響應(yīng)

圖5 車體質(zhì)量變化

圖6 車體位移PSD

圖7 車體垂直振動加速度PSD
低軌道路面激勵模型如圖8所示。仿真時間為10 秒。圖9所示為車體的垂直振動位移仿真結(jié)果,圖10所示為車體的垂直振動加速度仿真結(jié)果,圖11所示為車體質(zhì)量的變化曲線,圖12所示為車體的垂直振動位移功率譜密度曲線,圖13所示為車體的垂直振動加速度功率譜密度曲線仿真曲線。
通過分析圖9至圖13可以發(fā)現(xiàn),控制器不僅在低干擾300 km/h 的情況下具有非常好的控制效果,在高干擾200 km/h的路面激勵下依然能夠保持非常好的控制效果。該控制器使得懸架系統(tǒng)發(fā)揮了更加明顯的作用,不僅能夠使得系統(tǒng)抗外部擾動的性能更強,還可以減小車體在運行過程中的位移以及振動。由此驗證了控制器的有效性。
本文在非線性空氣懸架的試驗基礎(chǔ)上,利用所得到的試驗數(shù)據(jù)建立2 自由度高速列車模型,進(jìn)而針對該模型進(jìn)行控制研究。在控制器設(shè)計時,既考慮到車體由于乘客數(shù)量的變化而出現(xiàn)的范圍波動,即參數(shù)的不確定性的情況,又兼顧了高速列車在運行過程中可能出現(xiàn)的擾動,因此所設(shè)計的自適應(yīng)反步控制器不僅能夠解決系統(tǒng)非線性帶來的影響,而且能夠解決系統(tǒng)參數(shù)的不確定性影響,還可以依靠切換控制的思想來解決系統(tǒng)外部擾動帶來的影響。仿真結(jié)果顯示控制器可以有效地提升懸架性能,證明了控制器的穩(wěn)定性,從而滿足了懸架系統(tǒng)的控制性能要求,使其能應(yīng)對不確定參數(shù)和擾動對系統(tǒng)的影響,達(dá)到控制懸架系統(tǒng)增強高速列車的平順性和乘坐舒適性的目的。

圖8 路面激勵輸入

圖9 車體垂直位移響應(yīng)

圖10 車體垂直振動加速度響應(yīng)

圖11 車體質(zhì)量變化

圖12 車體位移PSD

圖13 車體垂直振動加速度PSD