陳 旺,肖克平,葉 振,徐 乾,胡周達(dá)
(1.粵電集團(tuán)貴州有限公司,貴陽(yáng) 550081;2.電子科技大學(xué) 自動(dòng)化工程學(xué)院,成都 611731 )
作為汽車(chē)傳動(dòng)系統(tǒng)中的關(guān)鍵部件,變速器承擔(dān)了將來(lái)自于發(fā)動(dòng)機(jī)引擎的力矩平滑、可靠地傳輸給驅(qū)動(dòng)車(chē)輪,進(jìn)而推動(dòng)汽車(chē)行駛的功能。變速器性能的優(yōu)劣直接關(guān)系到整車(chē)的動(dòng)力性、舒適性以及燃油經(jīng)濟(jì)性,是汽車(chē)系統(tǒng)的核心部件。一直以來(lái),變速器都是汽車(chē)傳動(dòng)系統(tǒng)中的研究熱點(diǎn)領(lǐng)域,得到眾多學(xué)者的密切關(guān)注。作為研究變速器性能的一種重要手段,根據(jù)動(dòng)力學(xué)模型和機(jī)械傳動(dòng)原理,對(duì)變速器的工作過(guò)程進(jìn)行合精確建模,是變速器工作機(jī)理、控制算法研究的根本條件[1-6]。
文獻(xiàn)[1]針對(duì)工程類(lèi)車(chē)輛的變速器換擋規(guī)律進(jìn)行研究,根據(jù)工程車(chē)輛的特點(diǎn),設(shè)計(jì)了最優(yōu)變速器控制規(guī)律。文獻(xiàn)[2]設(shè)計(jì)了一種模塊化變速器以及配套換擋液壓系統(tǒng),并對(duì)新模塊仿真建模研究。文獻(xiàn)[3] 圍繞雙參數(shù)換檔規(guī)律的標(biāo)定過(guò)程展開(kāi)研究,對(duì)整車(chē)傳動(dòng)系統(tǒng)進(jìn)行了動(dòng)態(tài)建模;文獻(xiàn)[4]針對(duì)換擋過(guò)程中零速區(qū)域進(jìn)行了研究和數(shù)值模擬問(wèn)題。這些仿真模型從不同的角度,針對(duì)變速器的自動(dòng)控制問(wèn)題進(jìn)行了深入的研究,為后續(xù)相關(guān)問(wèn)題的研究奠定了基礎(chǔ)。隨著新能源技術(shù)的發(fā)展,新能源汽車(chē)[7-8]、新能源發(fā)電[9-10]賦予了傳統(tǒng)發(fā)動(dòng)機(jī)新的應(yīng)用場(chǎng)景,也為自動(dòng)變速器研究成果提供了新的應(yīng)用土壤。由自動(dòng)變速器的功能可知,變速器是將發(fā)動(dòng)機(jī)引擎輸出的動(dòng)力,轉(zhuǎn)變?yōu)檫m應(yīng)外界動(dòng)力需求的裝置,其主要功能就是動(dòng)力的適配控制,因此,文獻(xiàn)[7-10]涉及到的研究問(wèn)題,與自動(dòng)變速器領(lǐng)域研究的問(wèn)題有異曲同工之處,很多研究方法和研究成果可以相互借鑒。
基于上述考慮,文章以汽車(chē)的自動(dòng)變速箱仿真系統(tǒng)設(shè)計(jì)為例,利用有限狀態(tài)工具箱Stateflow,快速搭建了汽車(chē)自動(dòng)變速箱仿真系統(tǒng),該方法摒棄了大量專(zhuān)業(yè)問(wèn)題的數(shù)學(xué)建模過(guò)程,將發(fā)動(dòng)機(jī)動(dòng)力輸出作為特定事件狀態(tài)對(duì)待,將變速器不同的狀態(tài)看作是多狀態(tài)事件轉(zhuǎn)移過(guò)程,利用Stateflow狀態(tài)機(jī)工具,實(shí)現(xiàn)了各個(gè)狀態(tài)的順利遷移和仿真模型的實(shí)現(xiàn),該方法的應(yīng)用可為相關(guān)問(wèn)題的建模仿真提供借鑒意義。
汽車(chē)自動(dòng)變速器的機(jī)械結(jié)構(gòu)比較復(fù)雜,但基本都是由液力變矩器、行星齒輪機(jī)構(gòu)以及電子控制結(jié)構(gòu)組成,其結(jié)構(gòu)如圖1所示。

圖1 汽車(chē)變速器結(jié)構(gòu)圖
圖1中,液力變矩器整個(gè)系統(tǒng)的核心組成部分,該裝置由泵輪、渦輪以及葉輪組成,通過(guò)相互作用,實(shí)現(xiàn)機(jī)械能向液體能的轉(zhuǎn)換,達(dá)到實(shí)現(xiàn)傳遞動(dòng)力的目的。該裝置的工作原理為:泵輪由發(fā)動(dòng)機(jī)曲軸驅(qū)動(dòng)旋轉(zhuǎn),導(dǎo)致泵輪葉片間的工作油液在離心力的作用下高速?lài)娝椭翜u輪,渦輪在該液體的沖擊下開(kāi)始旋轉(zhuǎn)運(yùn)動(dòng);此時(shí),高速液流不但沿著渦輪葉片通道內(nèi)流動(dòng),還在渦輪的牽連下旋轉(zhuǎn)運(yùn)動(dòng)。高速旋轉(zhuǎn)使渦輪葉片產(chǎn)生慣性,該慣性力與沖擊力一起作用,使渦輪獲得更大的輸出力矩。由渦輪流出的液體接著進(jìn)入不能旋轉(zhuǎn)的導(dǎo)輪,在導(dǎo)輪葉片的作重新返回泵輪,完成工作液體的一次循環(huán)。由上述工作過(guò)程可知,當(dāng)泵輪轉(zhuǎn)速不變時(shí),輸出軸轉(zhuǎn)矩具有跟隨渦輪轉(zhuǎn)速變化而變化的能力,因此液力變矩器是一種能夠自動(dòng)適應(yīng)汽車(chē)行駛阻力而自動(dòng)變矩出力的機(jī)械裝置,本文用RTQ表示表示這種裝置的變扭比。
由于受到空間限制,變矩器雖能增加發(fā)動(dòng)機(jī)扭矩,但其變矩比不大,不能適應(yīng)大范圍的變速,不足以滿(mǎn)足正常車(chē)輛需求,因此需要借助其它手段擴(kuò)充其變速范圍。行星齒輪變速結(jié)構(gòu)就是一種可以大范圍變速的傳動(dòng)裝置,其傳動(dòng)比RTR可以通過(guò)不同齒輪的咬合進(jìn)行改變,非常適合于動(dòng)力換檔或自動(dòng)換檔。
電子控制機(jī)構(gòu)通過(guò)外部傳感器(主要是車(chē)速和節(jié)氣門(mén)開(kāi)度信號(hào))采集車(chē)輛的實(shí)時(shí)行駛狀態(tài),并對(duì)這些信號(hào)進(jìn)行自動(dòng)處理,根據(jù)預(yù)先設(shè)置好的換檔條件,自動(dòng)確定檔位與換擋點(diǎn),并將相應(yīng)指令輸出至電磁閥線圈,通過(guò)對(duì)電磁閥線圈通斷的控制,實(shí)現(xiàn)油路的自動(dòng)切換,繼而達(dá)到自動(dòng)換檔的目的。
根據(jù)自動(dòng)變速器的工作原理,對(duì)其變速過(guò)程建模如下:
IeiNe=Te-Ti
(1)
其中:Te為發(fā)動(dòng)機(jī)的輸出扭矩,該扭矩通過(guò)曲軸與泵輪連接,帶動(dòng)泵輪旋轉(zhuǎn),發(fā)動(dòng)機(jī)的轉(zhuǎn)速用Ne表示,由于發(fā)動(dòng)機(jī)與泵輪之間,存在曲軸等連接裝置,導(dǎo)致泵輪的轉(zhuǎn)速與發(fā)動(dòng)機(jī)轉(zhuǎn)速之間不等,其大小由IeiNe表示,其中l(wèi)ei為轉(zhuǎn)動(dòng)慣量系數(shù),Ti為泵輪的轉(zhuǎn)矩,其值可由(2)式計(jì)算得出。
(2)
式(2)中,K為發(fā)動(dòng)機(jī)的能容系數(shù),通常由發(fā)動(dòng)機(jī)葉片的角度,個(gè)數(shù)等因素決定。
發(fā)動(dòng)機(jī)將動(dòng)力傳至泵輪,泵輪再通過(guò)液體經(jīng)動(dòng)力傳導(dǎo)至渦輪,帶動(dòng)渦輪旋轉(zhuǎn),定義渦輪的轉(zhuǎn)速與發(fā)動(dòng)機(jī)輸入的轉(zhuǎn)速為變速器的變扭比,其值可由式(3)計(jì)算,Nin為渦輪的轉(zhuǎn)速。
(3)
動(dòng)力經(jīng)液力變速器傳至行星齒輪處,行星齒輪輸出為:
Tout=RTRTin
(4)
Nout=RTRNin
(5)
其中:Tin,Tout分別為齒輪機(jī)構(gòu)的輸入、輸出扭矩;Nin,Nout分別為齒輪機(jī)構(gòu)的輸入、輸出轉(zhuǎn)矩;RTR為齒輪傳動(dòng)機(jī)構(gòu)的變速比。
IvNw=Rfd(Tout-Tload)
(6)
其中:Iv為車(chē)輛的慣量,Nw為車(chē)輛的角速度,Rfd為最終轉(zhuǎn)速比;Tload是與車(chē)速、空氣阻力以及剎車(chē)狀態(tài)相關(guān)的阻力矩,其由式(7)構(gòu)成:
Tload=Rload0+Rload1V2+Tbrake
(7)
Rload0分別為地面對(duì)車(chē)輛的摩擦力,Rload1為空氣對(duì)車(chē)輛阻力,其大小與車(chē)速相關(guān)。
式(1)~(7)描述了發(fā)動(dòng)機(jī)動(dòng)力傳導(dǎo)至驅(qū)動(dòng)輪,推動(dòng)汽車(chē)行駛的基本過(guò)程。實(shí)際的工作過(guò)程與此類(lèi)似,但涉及得因素更多,為了簡(jiǎn)化建模過(guò)程,僅用式(1)~(7)代表整個(gè)過(guò)程。
除此以外,在實(shí)際的建模中,還有一些關(guān)于發(fā)動(dòng)機(jī)的數(shù)據(jù)可供參考,這些數(shù)據(jù)是發(fā)動(dòng)機(jī)在出廠時(shí),利用臺(tái)架測(cè)試得到的一些數(shù)據(jù)。雖然在實(shí)際的應(yīng)用過(guò)程中,隨著使用條件的不同,這些測(cè)試數(shù)據(jù)會(huì)有一些變動(dòng),但其變化規(guī)律大致一致,在本文中,參照某發(fā)動(dòng)機(jī)廠商提供的變速切換點(diǎn)圖(圖2)和發(fā)動(dòng)機(jī)工作特性(圖3)作為建模仿真的數(shù)據(jù)。

圖2 變速器切換點(diǎn)

圖3 發(fā)動(dòng)機(jī)輸出扭矩MAP圖
作為整個(gè)系統(tǒng)的動(dòng)力輸出源,對(duì)發(fā)動(dòng)機(jī)精確建模的意義不言而喻,通常,發(fā)動(dòng)機(jī)特性曲線是一個(gè)非常重要的性能指標(biāo),它表征發(fā)動(dòng)機(jī)運(yùn)行狀態(tài)及運(yùn)行性能,由圖3可知,發(fā)動(dòng)機(jī)的輸出轉(zhuǎn)矩是發(fā)動(dòng)機(jī)轉(zhuǎn)速、油門(mén)開(kāi)度的函數(shù),在仿真中,可以根據(jù)上述數(shù)值進(jìn)行曲線擬合,得到一函數(shù)表達(dá)式,然后對(duì)該表達(dá)式進(jìn)行建模仿真,但這種方法需要消耗較多的計(jì)算資源,仿真速度較慢,且存在一定誤差。考慮到計(jì)算機(jī)仿真時(shí),均是按照一定的變量步長(zhǎng)進(jìn)行仿真,因此,在本文的仿真中,將圖3中的變量變化采取離散化處理,基于離散化的自變量,算出對(duì)應(yīng)的函數(shù)值,并將這些數(shù)據(jù)存放在一特定的表格中,在仿真時(shí),調(diào)用Simulink中提供的Look-up模塊,實(shí)現(xiàn)自變量不同的變化時(shí),自動(dòng)查找對(duì)應(yīng)的函數(shù)值,在文章的仿真中,模塊2D,1D均是采用這種方式構(gòu)建,通過(guò)這種處理方式,加快了模型的仿真進(jìn)度,提高了仿真的準(zhǔn)確性。圖4是由圖3以及式(1)實(shí)現(xiàn)的仿真模塊。

圖4 動(dòng)力輸出仿真單元
針對(duì)式(2)~(6)所表達(dá)的物理量的連接關(guān)系,建立如圖5所示的仿真模型:

圖5 變矩器仿真單元
自動(dòng)變速器在實(shí)際工作中,根據(jù)圖2所示的曲線圖,實(shí)時(shí)調(diào)整相應(yīng)電磁閥接通或斷開(kāi)離合器、制動(dòng)器等、使離合器和某檔齒輪嚙合改變動(dòng)力傳遞中的齒輪副,從而改變檔位。在本文中,選取四擋,二參數(shù)的自動(dòng)變速器作為研究控制對(duì)象,將每個(gè)擋位視為為一個(gè)狀態(tài),狀態(tài)之間的遷移根據(jù)圖2所示的條件進(jìn)行改變,文章借助Matlab中的有限狀態(tài)工具箱Stateflow工具箱實(shí)現(xiàn)上述狀態(tài)的遷移。Stateflow可以創(chuàng)建的復(fù)雜控制邏輯,并與Simulink模型有效地結(jié)合。其具體實(shí)現(xiàn)如圖6所示。

圖6 檔位變換仿真單元
作為動(dòng)力作用的最終目標(biāo),車(chē)輛行駛是檢驗(yàn)前述控制單元是否正常工作的最終體現(xiàn),該模塊的仿真模型建立,充分考慮了路面的阻力以及空氣對(duì)車(chē)輛的阻力因素,通過(guò)最終轉(zhuǎn)速的變化,檢驗(yàn)各仿真模塊是否達(dá)到了設(shè)計(jì)要求,仿真模型如圖7所示。

圖7 車(chē)輛行駛仿真單元
文章按照?qǐng)D2、圖3所示的而汽車(chē)發(fā)動(dòng)機(jī)參數(shù)進(jìn)行仿真,仿真過(guò)程中參照文獻(xiàn)[11]中的參數(shù)設(shè)置,補(bǔ)充了車(chē)的質(zhì)量,迎風(fēng)面積等必要參數(shù),具體設(shè)置如下:汽車(chē)質(zhì)量取1 260 kg,迎風(fēng)面積取2.96 m2,車(chē)輪轉(zhuǎn)動(dòng)慣量取62.57 kg·m2,飛輪轉(zhuǎn)動(dòng)慣量取1.87 kg·m2,空氣阻力系數(shù)0.25,傳動(dòng)效率取0.85,按照上述參數(shù)對(duì)仿真模型進(jìn)行初始化后,仿真正常行駛、加速超車(chē)以及剎車(chē)3種工況狀態(tài)下車(chē)輛控制的仿真效果。

圖8 正常行駛時(shí)曲線
圖8~10分別為正常行駛、加速超車(chē)以及剎車(chē)時(shí)3種工況狀態(tài)下車(chē)輛行駛時(shí)的參數(shù)曲線圖,圖中分別繪制了油門(mén)開(kāi)度、引擎轉(zhuǎn)速、輸出轉(zhuǎn)速以及汽車(chē)速度之間的關(guān)系。圖8為正常行駛時(shí),隨著油門(mén)開(kāi)度的增加,發(fā)動(dòng)機(jī)轉(zhuǎn)速也逐漸增加,當(dāng)t=2 s的時(shí)間時(shí),車(chē)速開(kāi)始上升7 km/h時(shí),發(fā)生1~2的向上換擋,換擋后轉(zhuǎn)速先是發(fā)生下降,然后逐漸增加,增加至2 000轉(zhuǎn)時(shí),發(fā)生2~3的向上換擋,此時(shí)車(chē)速達(dá)到18 km/h,從仿真曲線中顯示的換擋時(shí)的車(chē)速可知,換擋時(shí)的車(chē)速位于圖2給出的換擋車(chē)速范圍內(nèi),從而表明了仿真模型能夠真實(shí)反映實(shí)際車(chē)輛的行駛過(guò)程。

圖9 超車(chē)時(shí)行駛曲線
圖9為超速時(shí),各參數(shù)的變化過(guò)程。在沒(méi)有剎車(chē)(剎車(chē)扭矩為零)的情況下,車(chē)速為零,發(fā)動(dòng)機(jī)轉(zhuǎn)速保持在1 000轉(zhuǎn)左右,當(dāng)司機(jī)在t=0時(shí)踩60%的油門(mén)時(shí),引擎立即以?xún)杀抖嗟乃俣茸龀龇磻?yīng),這時(shí)整個(gè)變矩器的速比很低,輸出轉(zhuǎn)矩比很大,車(chē)輛加速快,發(fā)動(dòng)機(jī)和車(chē)輛都獲得速度,直到大約t=2 s,這時(shí)發(fā)生1~2向上升檔,升檔后,發(fā)動(dòng)機(jī)的速度典型地突然下降,然后恢復(fù)它的加速度逐漸上升,發(fā)生2~3和3~4的上升分別發(fā)生在4秒和8秒左右,曲線對(duì)應(yīng)的換擋速度與圖2所示的換擋范圍基本一致。

圖10 剎車(chē)時(shí)行駛曲線
圖10為剎車(chē)狀態(tài)下,汽車(chē)各參數(shù)的變化曲線。從t=0 s時(shí),汽車(chē)以油門(mén)開(kāi)度80%開(kāi)始行使,分別在車(chē)速為37 km/h,60 km/h時(shí),發(fā)生2~3,3~4檔的升檔操作,在t=5 s時(shí)對(duì)車(chē)輛施加一剎車(chē)信號(hào),車(chē)輛在剎車(chē)信號(hào)的作用下,車(chē)速由于慣性的原因,開(kāi)始下降,但發(fā)動(dòng)機(jī)在制動(dòng)力的作用下,立即從高轉(zhuǎn)速降至怠速轉(zhuǎn)速,實(shí)現(xiàn)了動(dòng)力的快速切斷。
文章首先分析了汽車(chē)傳動(dòng)裝置的工作原理,針對(duì)傳動(dòng)過(guò)程中主要物理量的變化搭建了簡(jiǎn)單的數(shù)學(xué)模型,接著基于Matlab/Simulink和Matlab/Stateflow工具箱搭建了整車(chē)動(dòng)力性能仿真模型和換檔過(guò)程的仿真模型,分別仿真模擬了正常行駛、超速以及剎車(chē)3種狀態(tài)下,模型的輸出參數(shù)特性,驗(yàn)證了模型的有效性。該方法主要通過(guò)將汽車(chē)行駛階段的不同狀態(tài),建模為多狀態(tài)事件的遷移,借助發(fā)動(dòng)機(jī)廠家提供的變速器換擋點(diǎn)圖以及發(fā)動(dòng)機(jī)工作特性圖,構(gòu)建供仿真過(guò)程調(diào)用的工作表,從而避免了繁瑣的數(shù)學(xué)建模和公式推導(dǎo),簡(jiǎn)化了模型復(fù)雜度,該方法的應(yīng)用可為相關(guān)復(fù)雜工程問(wèn)題的建模仿真提供借鑒意義。