于 今 陳 華 劉駿豪

1.重慶大學(xué)機(jī)械傳動(dòng)國家重點(diǎn)實(shí)驗(yàn)室，重慶，400044 2.重慶大學(xué)機(jī)械工程學(xué)院，重慶，400044

0 引言

液壓機(jī)械無級變速器(hydro-mechanical continuously variable transmission,HMCVT)是一種復(fù)合型無級變速器，由機(jī)械路與液壓路組成，其中機(jī)械路主要用于傳遞大功率和擴(kuò)大傳動(dòng)范圍，液壓路主要用于實(shí)現(xiàn)無級調(diào)速與傳遞小功率[1-3]。HMCVT綜合了液壓傳動(dòng)無級變速和機(jī)械傳動(dòng)效率高的兩方面特點(diǎn)，已經(jīng)應(yīng)用在部分特種車輛上[2,4-6]。同金屬帶式無級變速器控制類似，HMCVT通過對目標(biāo)速比的跟蹤控制，間接控制發(fā)動(dòng)機(jī)的轉(zhuǎn)速，使發(fā)動(dòng)機(jī)工作在最佳工作點(diǎn)，因此具備一個(gè)穩(wěn)定、有效的速比跟蹤控制策略顯得尤為重要[1,6-8]。

國內(nèi)外部分高校和研究單位對HMCVT的速比控制策略進(jìn)行了相關(guān)研究。韓國成均館大學(xué)進(jìn)行了發(fā)動(dòng)機(jī)-液壓機(jī)械傳動(dòng)最優(yōu)控制曲線研究[9]。文獻(xiàn)[1]提出了離散控制系統(tǒng)，文獻(xiàn)[6，8，10]的研究集中在模糊PID的控制。由于模糊控制不依賴于完整的精確數(shù)學(xué)模型，所以模糊控制在HMCVT這類復(fù)雜的系統(tǒng)中應(yīng)用較為廣泛[1,7,8,10],但是由于復(fù)合變速器的速比跟蹤控制較為復(fù)雜，常規(guī)模糊PID控制仍然存在速比跟蹤控制效果不佳、適應(yīng)性不強(qiáng)和動(dòng)態(tài)品質(zhì)不佳等問題[10-11]。為解決以上問題，本文設(shè)計(jì)了基于量化因子和決策因子的變論域模糊PID控制器并搭建了AMESim和Simulink聯(lián)合仿真平臺。

1 液壓機(jī)械無級變速器數(shù)學(xué)模型

1.1 基本傳動(dòng)結(jié)構(gòu)

本文所研究的HMCVT傳動(dòng)結(jié)構(gòu)簡圖見圖1。該系統(tǒng)由1個(gè)行星排、1個(gè)變排量液壓元件和1個(gè)定排量液壓元件構(gòu)成，當(dāng)系統(tǒng)工作時(shí)，由a1軸輸入的功率經(jīng)過行星排A分流，一路經(jīng)過齒輪傳動(dòng)到由變量泵-定量馬達(dá)構(gòu)成的容積調(diào)速液壓回路，一路經(jīng)過齒輪進(jìn)入機(jī)械傳動(dòng)部分，最后兩路功率經(jīng)過a2軸匯流并輸出到負(fù)載。圖1中箭頭的方向表明了HMCVT在不發(fā)生功率循環(huán)下的功率流向。通過調(diào)節(jié)變量泵的排量，進(jìn)而調(diào)節(jié)了液壓路的傳動(dòng)比，實(shí)現(xiàn)整個(gè)HMCVT的速比的改變。

圖1 液壓機(jī)械無級變速器基本結(jié)構(gòu)Fig.1 The basic structure of HMCVT

1.2 變速器基本模型

1.2.1軸系模型

根據(jù)動(dòng)力學(xué)的基本建模方法，可以得到各軸傳遞的扭矩和轉(zhuǎn)速。

T?1=z?1T?13T?2=z?2T?52

ω?3=z?1ω?4

(1)

ω?5=z?2ω?2

(2)

ω?1+Aω?3-1+A?ω?2=0

(3)

A?=-k?/(1+k?)

其中，ai?為第i?軸的表示符號；T?i?為第i?軸的輸入扭矩；T?ij?為i?軸齒輪到j(luò)?軸齒輪傳遞的扭矩；ω?i?為第i?軸的轉(zhuǎn)速；J?i?為i? 軸的轉(zhuǎn)動(dòng)慣量；B?i?為i? 軸的黏性阻尼；z?1為a3軸與a4軸的傳動(dòng)比；z?2為a2軸與a5軸的傳動(dòng)比；T?p為變量泵的輸入扭矩；T?m為馬達(dá)的輸出扭矩；T?為負(fù)載力矩。k?為行星排特性參數(shù)，其取值范圍為2.5～4[3]。聯(lián)立式(1)～式(3)可得變速器靜態(tài)速比i?H表達(dá)式：

i?H=m?/[1+A?m?-Az?1z?2]

其中，m?為變量泵與定量馬達(dá)的排量比。

1.2.2變量泵-定量馬達(dá)模型

針對本文采用的變量泵-定量馬達(dá)系統(tǒng)，可以得到定量馬達(dá)輸出轉(zhuǎn)速ω?m關(guān)于泵斜盤擺角θ?(等效為排量比m?)和輸出扭矩T?m的傳遞函數(shù):

其中，ω?n為液壓部分固有頻率；ξ?n為阻尼比;ω?p為泵的輸入轉(zhuǎn)速；D?m為液壓泵的排量；ω?m為馬達(dá)的輸出轉(zhuǎn)速；θ?為變量泵變量機(jī)構(gòu)的擺角，在本文中將擺角等效為排量比；k?p為排量梯度；β?e為油液的有效體積的彈性模量；J?t為液壓馬達(dá)及負(fù)載的轉(zhuǎn)動(dòng)慣量之和；V?0為泵和馬達(dá)的工作容積；C?t為泵和馬達(dá)總泄漏系數(shù)；B?m為黏性阻尼。

2 變論域模糊PID控制器的設(shè)計(jì)

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

變論域模糊PID控制是在常規(guī)模糊PID控制的基礎(chǔ)上，加入變論域模糊控制的思想，解決一般模糊PID在控制過程中存在的缺點(diǎn)和不足[11-12]。HMCVT變論域模糊PID控制框圖見圖2。這里采用一種基于量化因子(K?e?和K?ec?)和決策因子(K?u?)收縮的方法來對模糊論域進(jìn)行伸縮。其中PID參數(shù)模糊控制器就是根據(jù)模糊規(guī)則表來設(shè)計(jì)的，控制系統(tǒng)的關(guān)鍵是要得到優(yōu)化的K?P、K?I和K?D。首先找到一組相對較優(yōu)的初始值，量化因子模糊控制器和決策因子模糊控制器通過系統(tǒng)偏差e?和偏差變化率ec?，分別進(jìn)行各自模糊規(guī)則表的查詢；然后再對量化因子和決策因子進(jìn)行修正，進(jìn)而達(dá)到對整個(gè)模糊系統(tǒng)進(jìn)行論域伸縮的目的；最后得到了一組PID參數(shù)的修正值ΔK?P、ΔK?I和ΔK?D，其中PID參數(shù)的優(yōu)化值、初始值和修正值存在如下關(guān)系：

圖2 變論域模糊PID控制器結(jié)構(gòu)框圖Fig.2 The structure of variable universe fuzzy PID controller

(4)

2.2 常規(guī)模糊PID控制器設(shè)計(jì)

以復(fù)合變速器的實(shí)際速比與目標(biāo)速比的偏差e?和偏差變化率ec?作為輸入變量，將它們的變化范圍定義為模糊集上的論域,e?、ec?對應(yīng)的模糊子集為{NB,NM,NS,ZO,PS,PM,PB}。定義3個(gè)輸出語言變量，分別是比例修正值ΔK?P、積分修正值ΔK?I、微分修正值ΔK?D，對應(yīng)的模糊子集分別如下式：

(5)

其中NB、NM、NS、ZO、PS、PM、PB分別代表負(fù)大、負(fù)中、負(fù)小、零、正小、正中、正大。

模糊規(guī)則是根據(jù)專家經(jīng)驗(yàn)知識庫得到的，根據(jù)K?P、K?I和K?D3個(gè)參數(shù)對系統(tǒng)動(dòng)態(tài)特性的影響，可以得到比例、積分、微分常數(shù)修正值所建立的模糊規(guī)則表[11-12]。表1給出了ΔK?P的模糊修正表(ΔK?I和ΔK?D同理可得)。

利用MATLAB中的Fuzzy工具箱，選擇控制器類型為Mamdani型，定義輸入輸出變量隸屬函數(shù)的類型為三角型，在Fuzzy編輯器中輸入49條模糊規(guī)則，設(shè)定非模糊化方法為重心平均法，可得ΔK?P隨系統(tǒng)偏差e?和偏差變化率ec?的變化規(guī)律(ΔK?I和ΔK?D同理可得)。

表1 ΔKP的模糊修正表

2.3 基于量化因子和決策因子變論域的實(shí)現(xiàn)

對量化因子(K?e?和K?ec?)和決策因子(K?u?)伸縮實(shí)際就是對輸入輸出論域的一種伸縮。量化因子K?e?和K?ec?增加，系統(tǒng)輸入論域收縮，輸入e?和ec?對系統(tǒng)的作用就增大；反之，輸入e?和ec?對系統(tǒng)的作用就減小。決策因子K?u?增加，系統(tǒng)輸出論域擴(kuò)張，模糊輸出變量u?對系統(tǒng)的作用增大；反之，u?對系統(tǒng)的作用減小。對e?和ec?輸入論域的伸縮只與自身變化有關(guān)，設(shè)輸入e?和ec?論域的伸縮因子分別為α?和β?，它們對應(yīng)的模糊子集為{B,M,S,ZE}，得到輸入論域伸縮因子的模糊規(guī)則如表2所示,其中B,M,S,ZE分別表示正大、正中、正小和零。

表2 伸縮因子α和β的模糊規(guī)則表

輸出論域的伸縮需要根據(jù)反映系統(tǒng)狀態(tài)的e?和ec?共同確定，當(dāng)e?和ec?都較大且符號相同時(shí)，說明系統(tǒng)實(shí)際速比與目標(biāo)速比相差較大，并且在逐漸偏離目標(biāo)速比，使得誤差進(jìn)一步擴(kuò)大，此時(shí)控制器需要增加輸出，擴(kuò)大輸出論域以便讓系統(tǒng)盡快回到目標(biāo)速比附近。當(dāng)e?和ec?都較大且符號相反時(shí)，說明系統(tǒng)實(shí)際速比與目標(biāo)速比相差較大，但是正在以較快的速度回到目標(biāo)速比附近，此時(shí)應(yīng)該保持當(dāng)前論域以防止系統(tǒng)出現(xiàn)振蕩等問題。當(dāng)e?較小、ec?較大時(shí)，說明實(shí)際速比與目標(biāo)速比相差不大，但是將要以很快的速度偏離目標(biāo)速比，此時(shí)應(yīng)該適當(dāng)擴(kuò)大輸出論域，增加輸出量，使系統(tǒng)不會(huì)偏離目標(biāo)值。設(shè)系統(tǒng)輸出論域的伸縮因子為γ?，對應(yīng)模糊子集γ?={ZE、VS、S、SB、M、B、VB}，得到輸出論域的伸縮因子γ?的伸縮模糊規(guī)則如表3所示，ZE、VS、S、SB、M、B、VB分別表示零、超小、正小、中小、正中、正大、超大。

表3 伸縮因子γ的模糊規(guī)則表

3 基于AMESim-MATLAB的聯(lián)合仿真

3.1 液壓機(jī)械無級變速器聯(lián)合仿真模型

目前HMCVT系統(tǒng)模型的建立主要有3種工具，分別是功率鍵合圖、MATLAB/Simulink和AMESim。功率鍵合圖和AMESim都是基于功率流的建模工具，但是AMESim為用戶提供了圖形化建模的方式，可以讓用戶從繁雜的數(shù)學(xué)建模中解放出來[13]。MATLAB/Simulink對HMCVT建模需要對其簡化，這樣反而降低了模型精度，但是Simulink為用戶提供了實(shí)用的控制算法工具箱[10]。綜上分析，采用基于AMESim和MATLAB聯(lián)合仿真方式。HMCVT模糊PID控制的AMESim和MATLAB聯(lián)合仿真模型如圖3所示，AMESim模型主要由發(fā)動(dòng)機(jī)模型、行星排、泵-馬達(dá)系統(tǒng)和負(fù)載模型構(gòu)成。負(fù)載模型模擬了汽車行駛的過程，仿真模型的基本參數(shù)設(shè)置如表4所示。變論域模糊PID控制器和模糊PID仿真模型如圖4、圖5所示。

3.2 液壓機(jī)械無級變速器速比跟蹤仿真

為了驗(yàn)證模型的調(diào)速特性的正確性，建立了圖6所示的HMCVT實(shí)驗(yàn)平臺。該實(shí)驗(yàn)平臺主要由動(dòng)力輸入裝置、HMCVT和加載裝置等構(gòu)成，通過調(diào)節(jié)排量比的大小對復(fù)合變速器輸出轉(zhuǎn)速進(jìn)行檢測，其中實(shí)驗(yàn)平臺的基本參數(shù)與仿真模型一致。排量比為0～1時(shí)，傳動(dòng)比變化的理論、仿真和實(shí)驗(yàn)曲線如圖7所示。從圖中可以看出，理論結(jié)果比仿真和實(shí)驗(yàn)結(jié)果偏高是因?yàn)槔碚撚?jì)算時(shí)忽略了液壓損失和機(jī)械損失，但是仿真速比變化基本與理論和實(shí)驗(yàn)速比變化吻合，由此驗(yàn)證了仿真模型調(diào)速特性的正確性。

圖3 基于AMESim液壓機(jī)械無級變速器聯(lián)合仿真模型Fig.3 The simulation model of HMCVT based on AMESim

表4 仿真模型基本參數(shù)

在驗(yàn)證仿真模型正確性后，借助圖3的仿真模型首先找到一組控制效果較好的PID參數(shù)(ΔK?P=8，ΔK?I=40，ΔK?D=0.02)，為了便于比較本文設(shè)計(jì)的變論域模糊PID控制器與常規(guī)模糊PID控制器，使以上兩個(gè)控制器的基本參數(shù)完全相同，觀察兩者在不同目標(biāo)信號下的控制效果。

圖4 變論域模糊PID仿真模型Fig.4 The simulation model of the variableuniverse fuzzy PID

圖5 模糊PID仿真模型Fig.5 The simulation models of the fuzzy PID

圖6 液壓機(jī)械無級變速器實(shí)驗(yàn)平臺Fig.6 The experiment platform of HMCVT

圖7 理論、仿真和實(shí)驗(yàn)速比變化曲線Fig.7 Speed ratio curves of theory, simulation, experiment

3.2.1階躍信號跟蹤

以0.75的速比作為階躍輸入信號，負(fù)載不變。兩種算法速比響應(yīng)控制效果如圖8a所示，速比跟蹤誤差如圖8b所示，最大超調(diào)量和調(diào)整時(shí)間對比如表5所示。由以上數(shù)據(jù)分析可得，在階躍信號跟蹤中，變論域模糊PID控制較模糊PID控制超調(diào)量更小，調(diào)整時(shí)間更短，速比跟蹤誤差更小。

(a)速比跟蹤比較

(b)速比跟蹤誤差比較圖8 階躍信號跟蹤效果Fig.8 The follow-up control effect of step single input

控制算法最大超調(diào)量(%)調(diào)整時(shí)間(s)模糊PID4.41.3變論域模糊PID1.60.7

3.2.2斜坡信號跟蹤

斜坡信號(速比變化范圍為0～1.1)跟蹤控制效果如圖9所示。從圖中看出，兩種控制算法在斜坡上升跟蹤階段都存在跟蹤誤差，而且隨著跟蹤速比的增大而增大，在穩(wěn)態(tài)階段不存在穩(wěn)態(tài)誤差，在斜坡下降跟蹤階段存在跟蹤誤差，跟蹤速比越小跟蹤誤差就越小。模糊PID控制的跟蹤誤差在(-0.06，0.075)范圍內(nèi)，而變論域模糊PID控制的跟蹤誤差在(-0.03,0.03)范圍內(nèi)，速比跟蹤誤差約降低了50%。

(a)速比跟蹤比較

(b)速比跟蹤誤差比較圖9 斜坡信號跟蹤效果Fig.9 The follow-up control effect of ramp signal input

3.2.3正弦信號跟蹤

設(shè)定頻率為0.25 Hz, 速比變化范圍為(0.1,0.9)的正弦信號輸入，跟蹤效果如圖10所示。兩種控制算法跟蹤速比都存在一定滯后和誤差，跟蹤誤差也呈正弦趨勢變化。從圖10b可以看出，變論域模糊PID控制的跟蹤誤差小于模糊PID的跟蹤誤差，速比跟蹤誤差降低了約50%。

(a)速比跟蹤比較

(b)速比跟蹤誤差比較圖10 正弦信號跟蹤效果Fig.10 The follow-up control effect of sine signal input

綜上分析，變論域模糊PID控制器比常規(guī)模糊PID器的速比跟蹤效果更好，速比跟蹤誤差更小，具有更優(yōu)的動(dòng)態(tài)響應(yīng)特性，表現(xiàn)出更好的適應(yīng)性和動(dòng)態(tài)品質(zhì)。

4 結(jié)論

(1)本文設(shè)計(jì)了一種輸出耦合式液壓機(jī)械無級變速器，分析了液壓機(jī)械無級變速器的調(diào)速特性和系統(tǒng)數(shù)學(xué)模型，設(shè)計(jì)了一種基于量化因子和決策因子伸縮變化的變論域模糊PID控制器，分析了量化因子和決策因子在模糊PID控制中的作用，并制定了量化因子和決策因子伸縮變化的模糊規(guī)則表。

(2)搭建了基于AMESim和MATLAB的液壓機(jī)械無級變速器聯(lián)合仿真平臺，通過實(shí)驗(yàn)平臺驗(yàn)證了模型的正確性。在聯(lián)合仿真平臺上進(jìn)行了階躍信號、斜坡信號和正弦信號的速比跟蹤測試，仿真結(jié)果表明：所設(shè)計(jì)的變域模糊PID控制器比常規(guī)模糊PID控制器速比跟蹤控制效果更好，速比跟蹤誤差更小，表現(xiàn)出更優(yōu)越的動(dòng)態(tài)跟蹤品質(zhì)。研究結(jié)果為液壓機(jī)械無級變速器控制策略提供了理論參考。