李恒通，石瀟杰，嚴 勃，石 峰，郭 鑫

(1.北京礦冶科技集團有限公司，北京 100160；2.北礦機電科技有限責任公司，北京 100160)

目前在地下礦開采領域，自主駕駛地下鏟運機的研究越來越廣泛，發(fā)動機模型是后續(xù)發(fā)動機與變矩器共同工作點、檔位控制、車速控制研究的基礎，發(fā)動機模型的建立對于后續(xù)開展自主行駛方面的研究具有重要意義。

發(fā)動機模型主要有外特性模型、負載特性模型、調速特性模型、萬有特性模型[1]。本文先建立了發(fā)動機調速特性模型，之后根據調速特性模型建立了負載特性模型。

1 建立發(fā)動機模型的方法

發(fā)動機的理論特性模型一般根據發(fā)動機各個參數，結合發(fā)動機各個活塞缸的燃燒特性模型，進行模擬，進而得到其理論特性模型。

但發(fā)動機在實際工作中，由于發(fā)動機的實際噴油量、空氣溫度及負載等各個因素的影響，發(fā)動機的理論模型與實際的工作狀態(tài)相差較大，故實際當中，都是以發(fā)動機在臺架試驗上作出的試驗數據為準[2]。

發(fā)動機在臺架上的試驗數據是穩(wěn)態(tài)下的發(fā)動機工作狀態(tài)，而實際當中，發(fā)動機大多數時候工作在非穩(wěn)態(tài)狀態(tài)下，故建立發(fā)動機動態(tài)工作特性模型，才能準確的反映發(fā)動機的實際工作狀態(tài)。

在我國，試驗條件較國外差，建立發(fā)動機的特性模型一般是對穩(wěn)態(tài)的發(fā)動機數據進行一定的修正得到發(fā)動機動態(tài)工作特性。

目前建立發(fā)動機特性數學模型，主要有兩種方法：有多項式擬合法[3]和函數疊加法。

1.1 多項式擬合法

記錄臺架試驗中發(fā)動機在不同油門開度下轉矩、轉速的數據，以不同的曲線用最小二乘法去擬合每一組油門開度下的轉矩、轉速數據，得到n條k次曲線，即發(fā)動機的特性模型

每條曲線如式(1)示：

Me=αt0+αt1n+αt2n2+…+αtknk

(1)

其中：Me—發(fā)動機轉矩，N·m；n—發(fā)動機轉速，r/min；αt0,αt1,αt2,αtk—擬合系數；k—擬合次數，一般取2～5。

根據在不同油門開度a1，a2，…，ak下的m組發(fā)動機轉速轉矩的數據，建立發(fā)動機的特性模型。

方法如下：

(1)每組油門開度下的轉速轉矩數據用一條曲線去擬合，則對每一開度都有Me=f(n)。

(2)從最低轉速到最高轉速按間隔取出m組轉速值，nj=nmin+int(j(nmax-nmin)/(m-1))，j=0，通過擬合求出m組轉速轉矩數據所對應的多項式，將多項式的系數提取出來構成的矩陣如下：

其中每行的數據即是相應的多項式的系數，多項式如式(2)：

Me=αt0+αt1n+αt2n2+…+αtknk-1

(2)

逐一將各油門開度下的轉矩與轉速擬合，便可得到如下公式：

AimNmk=Me

(3)

其中：

1.2 疊加函數法

將發(fā)動機調速器的作用考慮，用函數去擬合發(fā)動機各個油門開度下的調速特性曲線[4]。

函數為兩個變量——油門開度和轉速，由三個不同類型的函數組成，M1為常函數，M2為正弦函數，M3為雙曲線函數，其在反映發(fā)動機的調速特性的同時還能反映發(fā)動機的外特性，同時能保證整個曲線段內可導[5]。即：

Me(ne,α)=M1+M2(ne)-M3(ne,α)

(4)

(5)

式中：Mem—發(fā)動機的最大轉矩，N·m；Me0—額定工作點轉矩，N·m；b—修正系數，修正調速段的斜率。

(6)

式中：ne0—發(fā)動機的額定轉速，r/min；nem—發(fā)動機的最大工作轉速，r/min。

雙曲線M3保證調速段的轉速的快速降低：

(7)

式中：b—調速段修正系數，b值越大，調速段斜率越大；ner1—發(fā)動機的最低怠速，r/min；nerz—發(fā)動機的最高轉速，r/min。

多項式擬合法適合于在臺架上測得各個油門開度下發(fā)動機的轉速轉矩數據，之后再根據多項式擬合法建立發(fā)動機模型；疊加函數法適用于已知發(fā)動機某個油門開度下的轉速轉矩數據，之后用這個油門開度下數據建立各個油門開度下的發(fā)動機模型，一般地下鏟運機出廠時會附帶發(fā)動機的外特性曲線模型即100%油門開度下的發(fā)動機輸出數據，故本次使用疊加函數法建立發(fā)動機模型。

2 發(fā)動機調速特性模型

發(fā)動機在工作當中，載荷時刻在變化，發(fā)動機轉速會隨著載荷的變化而變化，但有時負荷的微小變化會引起轉速的大幅度變化；尤其是發(fā)動機在最大油量位置，阻力矩突然減小時，柴油機會因超速“飛車”；當發(fā)動機在小油量位置，當轉矩增大，轉速降低時，發(fā)動機有可能熄火。為防止發(fā)動機高速時“飛車”，低速時熄火，故在發(fā)動機上一般都安裝有調速器以提高發(fā)動機對隨機的外載荷的適應能力[6]。

本次研究根據的對象為兩立方地下鏟運機，根據行走和鏟裝的工況下的功率計算，結合井下特點，進行匹配，選擇了康明斯QSB C130型電噴水冷增壓發(fā)動機，最大功率為97 kW@2 500 r/min，最大扭矩為600 Nm@1 500 r/min，發(fā)動機允許最小怠速為600 r/min，實際發(fā)動機怠速為800 r/min。

康明斯QSB C130型發(fā)動機安裝的是全程式調速器，可以對發(fā)動機各個油門開度下的任意轉速起調節(jié)作用，使發(fā)動機在整個工作轉速區(qū)間穩(wěn)定運轉。發(fā)動機的原始數據如下表1：

表1 發(fā)動機外特性

將帶有全程式調速器的QSB4.5-C130發(fā)動機的參數：Mem=600 N·m；Meo=370 N·m；b=5.5；nem=1 500 r/min；neo=2 500 r/min；ner1=600 r/min；ner2=2 700 r/min代入，得出發(fā)動機穩(wěn)態(tài)調速數學模型為：

將α等于不同值時，即油門開度位于不同位置時，發(fā)動機的特性曲線圖如圖2所示。

圖1 疊加函數法擬合圖Fig.1 Fitting curves of superposition function method

令α=100%，得出柴油機外特性模型：

通過與實際臺架測得的數據點進行比對，在實際工作區(qū)間(一般為1 500 r/min以上)，鏟運機外特性模型非常接近實際工作點，擬合效果好。

圖2 外特性曲線圖Fig.2 Engine external characteristic curve

3 負載特性模型

發(fā)動機的外特性是發(fā)動機廠商給出的，是在發(fā)動機未帶任何負載下的輸出轉矩和轉速的數據[4,7]。

發(fā)動機在實際工作中，允許的最大轉速會比額定轉速大10%左右，同時在實際計算發(fā)動機所能提供給車輛的動力時，應該扣除發(fā)動機自身附件所消耗的發(fā)動機功率，同時應該扣除車輛上制動油泵、轉向油泵以及工作油泵在卸荷狀態(tài)下消耗的發(fā)動機的功率[8]。

Mej=Me(1-α)

(8)

Mej-Mf=Mb

(9)

ne=nb

(10)

Pej=P-Pf

(11)

(12)

式中:Me—發(fā)動機扭矩，N·m；α—發(fā)動機所帶的附件對發(fā)動機功率的消耗，一般取10%；Mej—發(fā)動機凈扭矩，N·m；Mf—發(fā)動機附件所消耗的發(fā)動機轉矩，N·m；Mb—泵輪的輸入力矩，N·m；Pe—發(fā)動機凈輸出功率，kW；P—發(fā)動機額定功率，kW；Pf—發(fā)動機附件所消耗的發(fā)動機功率，kW；ne—發(fā)動機的轉速，r/min；nb—泵輪轉速，r/min。

鏟運機一般行駛狀態(tài)下，各個油泵如工作泵、制動泵及轉向泵均是卸荷狀態(tài)壓差為ΔP=1.31 MPa。

則總排量q=q1+q2+q3=74.6 mL/r。

式中：q1—工作泵排量，32.3 mL/r；q2—轉向油排量，32.3 mL/r；q3—制動泵排量，10.0 mL/r。

其消耗的發(fā)動機扭矩之和為：

(13)

式中：ip—變矩器油泵輸出機構速比取0.951；ηp—油泵機械效率，取0.9；

雙變系統(tǒng)油泵流量為79.5 L/min，ΔP=2.2 MPa，所消耗的發(fā)動機的輸出扭矩Mf2為：

式中：neo—發(fā)動機一般工作轉速；Q—雙變系統(tǒng)油泵流量，79.5 L/min。

Mf=Mf1+Mf2=18.2+14.1=32.3 N·m

發(fā)動機實際輸出扭矩為：

Mej=Me(1-α)=0.9Me

Mb=Mej-Mf=0.9Me-32.3

故可得發(fā)動機的負載特性模型如下：

(14)

具體結果：

圖3 發(fā)動機負載特性曲線Fig.3 Engine load characteristic curves

圖4 各油門下發(fā)動機負載特性曲線圖Fig.4 Engine load characteristic curves at different throttling levels

4 結論

本文通過分析現(xiàn)有的發(fā)動機建模方法，采用了疊加函數法建立了地下鏟運機發(fā)動機調速系統(tǒng)特性模型，通過與發(fā)動機臺架試驗數據對比，誤差較小，結果可信。之后根據調速系統(tǒng)特性模型同時考慮油泵等輔件的功率消耗，建立了發(fā)動機負載特性模型。上述模型可作為研究地下鏟運機動力傳動系統(tǒng)模型的基礎，為后續(xù)地下鏟運機自主行駛方面的研究提供一定的借鑒。