龐 輝， 劉 凡， 王 延

(西安理工大學(xué)機(jī)械與精密儀器工程學(xué)院 西安,710048)

引 言

車輛懸架系統(tǒng)的主要功能是吸收由路面不平度引起的振動響應(yīng)，同時保證輪胎與路面的不間斷接觸以提高車輛行駛平順性和操縱穩(wěn)定性[1-2]。由于半主動懸架(semi-active suspension，簡稱SAS)兼顧被動懸架(passive suspension，簡稱PS)的穩(wěn)定性和主動懸架的自適應(yīng)能力，能夠根據(jù)行駛工況實(shí)時調(diào)整阻尼力大小，實(shí)現(xiàn)與主動懸架相近的減振效果，同時無需過多的能量消耗，因而受到了研究人員廣泛關(guān)注，磁流變減振器(magneto-rheological，簡稱MR)作為非線性元件被廣泛應(yīng)用于半主動懸架系統(tǒng)中[3]。但是因?yàn)槠浞蔷€性以及磁滯特性，該減振器力學(xué)模型難以精確描述[4-5]。Stanway等[6]建立了磁流變減振器的Bingham模型，該模型能很好地描述磁流變減振器阻尼力的時域特性及阻尼力與位移的關(guān)系特性，但無法表示在速度較小的區(qū)域內(nèi)阻尼力與速度的關(guān)系。Spencer等[7]提出了修正的Bouc-Wen模型，它能較好地描述磁流變減振器的非線性特性，但模型中有過多的參數(shù)需要優(yōu)化辨識，同時由于模型復(fù)雜，在實(shí)際應(yīng)用過程中會導(dǎo)致實(shí)時性降低，時滯增大。文獻(xiàn)[8]針對參數(shù)化建模存在的問題建立了一種基于自適應(yīng)神經(jīng)模糊的磁流變減振器非參數(shù)化模型，同時驗(yàn)證該模型可以高度逼近磁流變減振器實(shí)驗(yàn)?zāi)Ｐ汀?/p>

針對半主動懸架系統(tǒng)控制器的設(shè)計(jì)與優(yōu)化，研究人員開展了大量研究，提出諸如H∞控制[9-10]，滑模變結(jié)構(gòu)控制[11]，自適應(yīng)反推控制[12]及基于T-S模糊模型的控制[13]等。其中，T-S模糊控制由于結(jié)構(gòu)簡單且控制策略不依賴于模型，近年來受到廣泛關(guān)注。例如，文獻(xiàn)[14]針對主動懸架系統(tǒng)控制執(zhí)行器的延遲和故障不確定性等問題，在滿足安全性能約束要求并保證系統(tǒng)漸近穩(wěn)定性和魯棒性前提下，設(shè)計(jì)了一種模糊魯棒H∞控制器。文獻(xiàn)[15]針對非線性半車懸架系統(tǒng)建立了T-S模糊模型，設(shè)計(jì)了半車主動懸架系統(tǒng)的優(yōu)化模糊控制器。文獻(xiàn)[16]在考慮車身振動及控制時滯條件下，基于T-S模糊模型設(shè)計(jì)了一種穩(wěn)定的模糊H∞控制器，實(shí)現(xiàn)了懸架系統(tǒng)的漸進(jìn)穩(wěn)定以及預(yù)期的H∞性能。上述文獻(xiàn)所設(shè)計(jì)的T-S模糊控制器，對懸架性能均有一定程度的改善，但傳統(tǒng)模糊控制存在過分依賴專家經(jīng)驗(yàn)的問題，因此導(dǎo)致所設(shè)計(jì)控制器控制精度難以保證。針對模糊控制存在的不足，文獻(xiàn)[17]提出了變論域模糊控制的思想，在規(guī)則不變的前提下，論域隨著誤差變小而收縮(亦可隨著誤差增大而擴(kuò)展)。論域的這種伸縮變化通過調(diào)整伸縮因子的大小來實(shí)現(xiàn)，目前對于伸縮因子的描述沒有統(tǒng)一形式。如文獻(xiàn)[18-19]采用函數(shù)模型來選擇伸縮因子大小，文獻(xiàn)[20-21]采用模糊模型來描述伸縮因子變化，但是兩種描述都存在一定不足，即函數(shù)模型難以在控制過程中精確描述伸縮因子變化，而采用模糊模型會導(dǎo)致伸縮因子的選擇完全依賴于經(jīng)驗(yàn)公式，產(chǎn)生與模糊控制同樣的問題。文獻(xiàn)[22]針對1/4車輛半主動懸架提出一種基于模糊神經(jīng)網(wǎng)絡(luò)的變論域控制策略，并取得了較好的控制效果，但并未研究分析該控制策略對于車身俯仰及其它運(yùn)動的控制有效性。

為此，針對依據(jù)實(shí)驗(yàn)數(shù)據(jù)所建立的磁流變減振器模型及1/2車輛半主動懸架系統(tǒng)，筆者提出一種改進(jìn)變論域模糊控制器設(shè)計(jì)方法。通過建立懸架系統(tǒng)的T-S模糊控制器，設(shè)計(jì)合適的模糊神經(jīng)網(wǎng)絡(luò)結(jié)構(gòu)對T-S模糊控制器的兩個輸入變量論域進(jìn)行伸縮調(diào)整，結(jié)合模糊控制的邏輯語言和神經(jīng)網(wǎng)絡(luò)的自學(xué)習(xí)能力，實(shí)現(xiàn)對伸縮因子的精準(zhǔn)描述。最后通過仿真分析，驗(yàn)證了筆者所提出控制方法的有效性。

1 基于實(shí)驗(yàn)數(shù)據(jù)的MR減振器阻尼特性建模

文中選取路虎極光2016款SUV車所搭載的MagneRide減振器作為研究對象，其內(nèi)部簡化結(jié)構(gòu)如圖1左圖所示，圖1右圖為筆者所采用的MagneRide減振器實(shí)體。內(nèi)部結(jié)構(gòu)主要包括磁流變液、活塞、導(dǎo)電線圈、浮動塞以及彈簧。通過調(diào)整導(dǎo)電線圈電流大小改變磁流變液中磁性分子排列方式，進(jìn)而實(shí)現(xiàn)減振器阻尼力連續(xù)可調(diào)。

圖1 MR減振器實(shí)體及內(nèi)部結(jié)構(gòu)示意圖Fig.1 Photograph and inner schematic configuration of the MR shock absorber

為了精確描述MR減振器的非線性磁滯特性，學(xué)者們研究了多種力學(xué)模型。其中自適應(yīng)神經(jīng)網(wǎng)絡(luò)模糊非參數(shù)化模型在避免大量參數(shù)識別的同時，能夠較為精確描述MR減振器阻尼特性[23]。其結(jié)構(gòu)如圖2所示，其中，系統(tǒng)Ⅰ確定不同電壓下阻尼力輸出等級K；系統(tǒng)Ⅱ描述特定電壓下阻尼力與活塞位移和速度之間的非線性關(guān)系，F(xiàn)=K×f為減振器輸出阻尼力。

圖2 MR減振器自適應(yīng)神經(jīng)模糊非參數(shù)化模型Fig.2 Non-parametric model of the MR shock absorber based on adaptive neuro-fuzzy

采用如圖3所示的最大輸出力為20 kN的電液伺服疲勞試驗(yàn)機(jī)對筆者所采用的MagneRide減振器進(jìn)行實(shí)驗(yàn)測試，通過分析實(shí)驗(yàn)數(shù)據(jù)并基于MatlabAnfis工具箱對實(shí)驗(yàn)數(shù)據(jù)訓(xùn)練建立如圖3所示的MR減振器自適應(yīng)神經(jīng)模糊模型。圖4為頻率2 Hz、振幅25 mm時MR減振器輸出阻尼力與活塞位移關(guān)系曲線以及對應(yīng)模型仿真結(jié)果。

圖3 MR減振器試驗(yàn)臺架Fig.3 Testing setup of MR damper

圖4 實(shí)驗(yàn)數(shù)據(jù)與模型響應(yīng)對比Fig.4 Comparison between the identified model and the measured data

2 半主動懸架系統(tǒng)動力學(xué)模型

1/2車輛懸架模型包含了車輛動力學(xué)分析的主要特征，且結(jié)構(gòu)簡單，所以在懸架控制策略的研究中被廣泛應(yīng)用。圖5為4自由度1/2車輛懸架系統(tǒng)模型，其中：Is為轉(zhuǎn)動慣量；a，b為前后軸到質(zhì)心的距離；θ為車身俯仰角；ms為折算到半車模型上的簧載質(zhì)量；mu f，mu r為非簧載質(zhì)量；kt f，kt r為輪胎剛度系數(shù)；zs，zu f，zu r，qf，qr分別為簧載質(zhì)量、非簧載質(zhì)量和路面的絕對位移；kf，kr為懸架的剛度。模型中用線性彈簧代替彈性輪胎，忽略輪胎阻尼；減振器固有阻尼及可變阻尼力代替磁流變減振器，因此前后磁流變減振器阻尼力可表示為

(1)

圖5 1/2車輛懸架系統(tǒng)模型Fig.5 The half vehicle dynamic model

如圖5所示的1/2車輛振動模型中，要求車輛相對于縱垂面完全對稱，而且左右車輪下的路面不平度完全一樣，即認(rèn)為車輛是在縱垂面上的振動，對其受力分析如下。

對于前后軸上方的垂直位移有如下關(guān)系

(2)

則懸架系統(tǒng)動力學(xué)方程為

(3)

3 變論域模糊控制系統(tǒng)設(shè)計(jì)

3.1 變論域模糊控制基本思想

變論域模糊控制的思想可以描述為：假定Xi=[-Ei,Ei](I=1,2,…,n)代表輸入變量xi(i=1,2,…,n)的模糊論域，Y=[-U,U]代表輸出變量y的模糊論域。同時{Ail}分別表示模糊語言變量Xi，Y的模糊集合[24]。具體描述如下。

1) 在規(guī)則不變的前提下，論域隨初始論域[-E,E]通過伸縮因子α(x)變換為[-a(x)E,a(x)E]，如圖6所示，a(x)為誤差變量x的連續(xù)函數(shù)。

2) 變論域的思想意味著變量Xi,Y的模糊論域隨著輸入輸出變量xi，y而變化，具體表達(dá)形式如式(4)，(5)所示。

其中：變量αi(xi)(i=1，2，…，n)，β(y)分別表示Xi，Y的伸縮因子。

另外，Xi，Y的初始論域定義為原始論域，且采用7個語言模糊子集{NB(負(fù)大)，NM(負(fù)中)，NS(負(fù)小)，ZE(零)，PS(正小)，PM(正中)，PB(正大)}。輸入輸出模糊論域伸縮變化如圖6所示。

圖6 變論域模糊控制基本思想Fig.6 The variation of contracting-expanding universe

文獻(xiàn)[25]采用模糊推理來描述伸縮因子的變化，其仿真結(jié)果表明，基于模糊推理的變論域模糊控制(variable-universe fuzzy control based on fuzzy reasoning，簡稱VUFC-FR)算法能夠有效地提高收斂速度和控制精度，但是這種控制算法伸縮因子的模糊論域劃分不夠精確，從而導(dǎo)致控制結(jié)果存在誤差。針對這一不足，筆者采用模糊神經(jīng)網(wǎng)絡(luò)結(jié)構(gòu)來控制伸縮因子的變化，結(jié)合模糊推理的語言邏輯性和神經(jīng)網(wǎng)絡(luò)的自學(xué)習(xí)調(diào)整優(yōu)勢，通過系統(tǒng)誤差分析以及神經(jīng)網(wǎng)絡(luò)的自學(xué)習(xí)調(diào)整修正伸縮因子的模糊論域劃分?；谀：窠?jīng)網(wǎng)絡(luò)的變論域模糊控制(variable-universe fuzzy control based on fuzzy neural network，簡稱VUFC-FNN)框圖如圖7所示，r和y分別代表理想的輸入和控制信號，I表示控制電流，F(xiàn)MR為減振器輸出阻尼力，e為系統(tǒng)控制誤差，并表示為e=r-y。

圖7 基于模神經(jīng)網(wǎng)絡(luò)的變論域模糊控制框圖Fig.7 The control block of VUFC based on FNN

3.2 基于T-S模糊模型的變論域控制器設(shè)計(jì)

T-S型模糊控制器可由如下方程組成：

Ri:如果xi為Ai1，x2為Ai2，…，xj為Aij，…，xn為Ain,那么

gi(X)=ai0+ai1x1+ai2x2+…+aijxj+…+ainxn

(6)

其中：Ri(i=1,2,…,R)為第i條規(guī)則；X=[x1,x2,…,xj,…,xn]T為模糊控制器的輸入矢量；xj(j=1, 2,…,R)為第j個輸入變量；Aij為模糊集合；gi(X)為第i條規(guī)則的輸出；f(X)為模糊控制器的輸出，采用中心平均解模糊方法；μi(X)為第i條規(guī)則的定義為乘積形式的滿足程度；Aij(xj)為xj對Aij的滿足程度；hij(xj)為定義在輸入論域上的隸屬度函數(shù)[26]。

T-S模糊控制器在設(shè)計(jì)過程中有兩個關(guān)鍵性問題：a.確定控制規(guī)則數(shù)目；b.確定規(guī)則后件的系數(shù)。確定規(guī)則數(shù)目的過程，即對輸入變量劃分的過程，它需要產(chǎn)生合理的子空間。確定控制規(guī)則后件系數(shù)也就是確定區(qū)域作用量的過程。解決這兩個關(guān)鍵問題步驟如下。

1) 最優(yōu)參數(shù)。行駛平順性和操縱穩(wěn)定性是汽車動力學(xué)研究的兩個重要領(lǐng)域，懸架系統(tǒng)在協(xié)調(diào)這兩大使用性能方面起著非常重要的作用。在分析懸架諸多影響因子中，選取3個基本參數(shù)對懸架性能作評價(jià)，分別為車身加速度、懸架動撓度和輪胎動載荷。為了同時滿足車輛行駛平順性和操縱穩(wěn)定性的要求，建立如式(10)所示的目標(biāo)函數(shù)J

(10)

3.3 模糊神經(jīng)網(wǎng)絡(luò)伸縮因子設(shè)計(jì)

圖8 模糊神經(jīng)網(wǎng)絡(luò)(FNN)的結(jié)構(gòu)Fig.8 The structure of FNN

網(wǎng)格的Ⅰ，Ⅱ?qū)訉?yīng)于模糊規(guī)則“if”前提，采用高斯型函數(shù)作為隸屬函數(shù)。

網(wǎng)格的Ⅴ對應(yīng)解模糊層，采用面積重心法實(shí)現(xiàn)解模糊化。該層的輸出節(jié)點(diǎn)f(5)為調(diào)節(jié)系統(tǒng)模糊論域的伸縮因子α1，α2，即

(16)

在實(shí)現(xiàn)伸縮因子調(diào)整的模糊神經(jīng)網(wǎng)絡(luò)FNN結(jié)構(gòu)中，輸入分量的模糊劃分個數(shù)mi是預(yù)先確定的，則需要學(xué)習(xí)的參數(shù)包括最后一層的連接權(quán)ωij(該層對應(yīng)于伸縮因子模糊控制器的輸出論域的中心值，i=1,2；j=1,2,…,m)和Ⅱ?qū)与`屬度函數(shù)的中心值cij和寬度bij(i=1,2；j=1,2,…,mi)[27]。

上述伸縮因子模糊神經(jīng)網(wǎng)絡(luò)是一種多層前饋網(wǎng)絡(luò)，因此可以采用誤差反轉(zhuǎn)的方法設(shè)計(jì)調(diào)整參數(shù)的學(xué)習(xí)算法，從而使系統(tǒng)的實(shí)際輸出逼近期望輸出。設(shè)yr和y分別表示期望輸出和實(shí)際輸出，定義誤差能量函數(shù)E=1/2(yr-y)2，則模糊神經(jīng)網(wǎng)絡(luò)的調(diào)整參數(shù)的自學(xué)習(xí)算法如式(17)所示

(17)

其中：η>0為模糊神經(jīng)網(wǎng)絡(luò)學(xué)習(xí)速率。

為了驗(yàn)證筆者提出的基于模糊神經(jīng)網(wǎng)絡(luò)的變論域模糊控制算法的有效性和穩(wěn)定性，使用Matlab/Simulink工具箱建立該控制算法的仿真模型，在隨機(jī)不平路面激勵下，選擇車身加速度、俯仰角加速度、懸架動行程及輪胎動載荷作為輸出測量變量。

4 仿真分析4.1 隨機(jī)不平路面模型

文中采用白噪聲法建立路面模型為

(18)

其中：f0為路面輸入的下截止頻率；n0為參考空間的頻率；Gq(n0)為路面不平度系數(shù)；u為車速；ω(t)為白噪聲。

Matlab/Simulink仿真所得隨機(jī)不平路面[28]如圖9所示。

圖9 B級路面激勵信號Fig.9 The B-class road excitation used for simulation

4.2 仿真驗(yàn)證

文中仿真采用的車輛磁流變半主動懸架參數(shù)見表1。

表1 1/2車輛磁流變半主動懸架系統(tǒng)模型參數(shù)

在控制仿真中，采用隨機(jī)路面激勵作為輸入信號，基于Matlab/Simulink工具箱的改進(jìn)變論域模糊控制器設(shè)計(jì)如圖10所示。

圖10 變論域模糊控制系統(tǒng)模型Fig.10 Simulation block diagrams of the enhanced VUFC

圖中：ke,kc為T-S模糊控制器比例因子；α1，α2為輸入變量伸縮因子。ANFIS采用模糊神經(jīng)網(wǎng)絡(luò)控制車身垂向加速度、車身相對速度的論域變化范圍。T-S模糊控制器輸出為控制電流I，同時I作為MR減振器模型輸入，得到控制阻尼力并作用于半主動懸架系統(tǒng)。

仿真條件為B級隨機(jī)路面，車速為20 m/s。

從圖11、圖12可以看出，相比于被動和VUFU-FR控制懸架，VUFC-FNN控制的MR半主動懸架車身加速度和俯仰角加速度優(yōu)化效果明顯； 同時在頻域響應(yīng)，車身加速度和俯仰角加速度功率譜密度也得到了較好改善。

圖11 車身加速度響應(yīng)Fig.11 Response to vehicle body acceleration

從圖13可以看出，VUFC-FNN控制算法可有效減小前后輪懸架動行程，防止懸架動行程超出限定范圍而導(dǎo)致懸架擊穿。從圖14可以看出，VUFU-FNN控制MR半主動懸架可有效改善路面作用于輪胎的沖擊載荷，提高輪胎壽命同時保證車輛操作穩(wěn)定性。圖15為前后輪MR減振器所產(chǎn)生的阻尼力。

圖12 車身俯仰角加速度響應(yīng)Fig.12 Responses to vehicle body pitch angular acceleration

圖13 懸架動行程響應(yīng)Fig.13 Responses to the suspension dynamic deflection

圖14 輪胎動載荷響應(yīng)Fig.14 Response to the tire dynamic loads

圖15 磁流變減振器輸出阻尼力Fig.15 Response to the MR damper force

為進(jìn)一步說明控制算法的有效性，表2列出了被動懸架、VUFC-FR及VUFC-FNN控制懸架在隨機(jī)不平路面擾動下車身加速度、俯仰角加速度、懸架動行程、輪胎動載荷的均方根值以及優(yōu)化的百分比。

表2 懸架性能指標(biāo)均方根值比較

Tab.2 Comparison of RMS values of suspension performance

性能PSVUFC-FRVUFC-FNN¨zs/(m·s-2)0.356 10.214 6(39.74%)0.121 7(62.82%)¨θ/(m2·s-3)0.071 20.060 5(15.02%)0.057 3(19.52%)Δyf /m0.006 70.004 3(35.82%)0.003 1(53.73%)Δyr /m0.006 20.003 8(35.82%)0.003 4(41.79%)Ftf /N241.780 8159.259 7(34.13%)87.404 0(62.85%)Ftr/N203.769 8142.435 6(30.10%)119.951 4(41.13%)

從表2可以看出，VUFC-FNN控制MR半主動懸架在隨機(jī)不平路面擾動下，各項(xiàng)性能指標(biāo)均有一定程度的提升，而且優(yōu)化效果明顯優(yōu)于VUFC-FR控制懸架，說明筆者所提出的VUFC-FNN控制算法能夠有效提高車輛行駛平順性和操作穩(wěn)定性。

5 結(jié)束語

通過建立1/2車輛MR半主動懸架系統(tǒng)動力學(xué)模型，考慮隨機(jī)路況條件的影響，研究了基于模糊神經(jīng)網(wǎng)絡(luò)的變論域模糊控制方法，并與被動懸架和基于模糊推理的變論域模糊控制進(jìn)行比較。研究結(jié)果表明，相比于被動懸架，筆者提出的VUFC-FNN控制MR半主動懸架系統(tǒng)各項(xiàng)性能指標(biāo)均有一定提升。車身加速度和俯仰角加速度分別優(yōu)化62.82%和19.52%；前后輪懸架動行程和輪胎動載荷分別改善53.73%，41.79%，62.85%和41.13%?；谀：窠?jīng)網(wǎng)絡(luò)的變論域模糊控制能適應(yīng)路面工況的變化，可有效改善車輛行駛的操作穩(wěn)定性和乘坐舒適性，且控制效果優(yōu)于基于模糊推理的變論域模糊控制。