鄭嘯洲+孫偉

摘 要：試驗數(shù)據(jù)處理是發(fā)動機特性試驗的一個重要環(huán)節(jié)，是對發(fā)動機指標(biāo)及動力性、經(jīng)濟性能分析的基礎(chǔ)。本文利用MATLAB擬合工具箱的數(shù)據(jù)處理和繪圖功能，采用最小二乘法對某型發(fā)動機試驗數(shù)據(jù)進行擬合處理，建立其數(shù)學(xué)模型。利用所建模型能夠有效地分析發(fā)動機的各項性能，為研究發(fā)動機與變速器及整車的合理匹配提供了理論基礎(chǔ)。

關(guān)鍵詞：MATLAB；擬合工具箱；發(fā)動機；試驗數(shù)據(jù)；分析；建模

DOI：10.16640/j.cnki.37-1222/t.2017.18.001

0 引言

發(fā)動機是車輛的動力核心部件，其特性對車輛的動力性、燃油經(jīng)濟性等性能影響很大。發(fā)動機要與整車合理匹配才能發(fā)揮其自身優(yōu)勢，匹配的前提在于建立精確的發(fā)動機數(shù)學(xué)模型[1]。發(fā)動機模型要能夠全面準(zhǔn)確的反映發(fā)動機特性，并對其性能進行相應(yīng)的計算分析。

發(fā)動機主要性能指標(biāo)（如有效功率Pe、有效扭矩Ttq、燃油消耗率be等）和相關(guān)的特性參數(shù)隨運行工況變化而變化的關(guān)系，稱為發(fā)動機特性，其中與車輛密切相關(guān)的有速度特性、負荷特性、萬有特性和調(diào)速特性等[2]。本文主要根據(jù)發(fā)動機的特性試驗結(jié)果，利用MATLAB擬合工具箱進行多項式擬合，建立發(fā)動機的模型，為發(fā)動機與變速器及整車的匹配和控制提供理論和技術(shù)基礎(chǔ)。

1 發(fā)動機速度特性

柴油發(fā)動機速度特性是指在油量調(diào)節(jié)機構(gòu)（油量調(diào)節(jié)齒條、拉桿）保持在某一位置（油門開度α不變）的情況下，主要性能指標(biāo)隨發(fā)動機轉(zhuǎn)速ne的變化規(guī)律[3]。當(dāng)油量控制機構(gòu)位于標(biāo)定位置時（α=1），測得的特性為全負荷速度特性，簡稱為外特性。油量控制機構(gòu)位于低于標(biāo)定位置的其他位置時（0<α<1）的速度特性，稱為部分負荷速度特性。

2 發(fā)動機調(diào)速特性

3 發(fā)動機萬有特性

為了能在一張圖上較為全面的表示各種工況下內(nèi)燃機行的參數(shù)的變化，形成多參數(shù)的特性曲線，即為發(fā)動機的萬有特性。發(fā)動機萬有特性曲線上的等油耗曲線和等功率曲線，實際上是在不同轉(zhuǎn)速、負荷條件下燃油消耗率或功率相等點連接而成的曲線，類似于等高線，其實質(zhì)是二維曲線表示三維曲面[6]。故可用多元線性回歸的方法直接擬合得到它們的三維曲面，然后再用一系列的平行于坐標(biāo)面的平面截三維曲面，其交線就是等油耗曲線和等功率曲線，即可得到萬有特性曲線。采用多項式擬合得發(fā)動機有效功率模型和燃油消耗率模型如式（5）和式（6）所示，其三維曲面如圖3和圖4所示。

對于上面擬合得到的發(fā)動機發(fā)動機燃油消耗率和有效功率曲面，用MATLAB中的畫等高線CONTOUR語句便可得到等油耗曲線和等功率曲線。將這兩種曲線再加上外特性曲線繪制在一張圖上便可得到發(fā)動機的萬有特性模型，如圖5所示。

4 發(fā)動機轉(zhuǎn)速調(diào)節(jié)特性

對于全程調(diào)速式柴油發(fā)動機，其每一個油門開度都有一個最大有效功率點與之對應(yīng)，以發(fā)動機油門開度為橫坐標(biāo)軸，轉(zhuǎn)速為縱坐標(biāo)，依次連接這些最大有效功率點，即可得到發(fā)動機最佳動力性工作曲線。同理，在發(fā)動機萬有特性圖上將具有相同功率的最小燃油消耗率點連接起來，在發(fā)動機穩(wěn)態(tài)輸出扭矩曲面上找到與之對應(yīng)的發(fā)動機油門開度，亦能得到以油門開度為橫坐標(biāo)，輸出轉(zhuǎn)速為縱坐標(biāo)的最佳燃油經(jīng)濟性工作曲線。該型柴油發(fā)動機最佳動力性和最佳經(jīng)濟性轉(zhuǎn)速調(diào)節(jié)特性曲線如圖6所示，式（7）為兩條調(diào)速特性曲線的數(shù)學(xué)表達式。

5 小結(jié)

本文利用MATLAB擬合工具箱對發(fā)動機性能試驗數(shù)據(jù)進行擬合處理，采用最小二乘法擬合建立了描述發(fā)動機特性的數(shù)學(xué)模型，包括發(fā)動機外特性曲線、輸出轉(zhuǎn)矩模型、萬有特性曲線以及轉(zhuǎn)速調(diào)節(jié)特性曲線。該方法能夠直接利用MATLAB工具箱中現(xiàn)有的函數(shù)，快速有效地擬合發(fā)動機試驗數(shù)據(jù)，建立的發(fā)動機模型精度較高，能夠有效地描述發(fā)動機性能。同時，該方法在發(fā)動機特性計算分析、建模以及和動力傳動系統(tǒng)的合理匹配等具有一定的通用性。

參考文獻：

[1]徐立友，張明柱，周志立等.拖拉機液壓機械無級變速傳動系統(tǒng)與發(fā)動機的合理匹配[J].農(nóng)業(yè)工程學(xué)報，2006，22（09）：109-113.

[2]孫軍.車輛內(nèi)燃機原理[M].北京：高等教育出版社，2012.

[3]董敬，莊志，常思勤.汽車拖拉機發(fā)動機[M].北京：機械工業(yè)出版，2000.

[4]徐立友，周志立，張明柱.基于MATLAB的柴油機性能試驗數(shù)據(jù)的處理[J].河南科技大學(xué)學(xué)報（自然科學(xué)版），2006，27（04）：33-35，54.

[5]徐立友，周志立，張明柱等.拖拉機液壓機械無級變速傳動系統(tǒng)速比匹配策略[J].中國農(nóng)業(yè)大學(xué)學(xué)報，2006，11（04）：94-98.

[6]易雪梅，吳伶.用MATLAB語言繪制發(fā)動機萬有特性的兩種方法[J].北京汽車，2005（05）：33-35.

項目基金：軍隊科研計劃項目

作者簡介：鄭嘯洲（1993-），男，安徽黃山人，碩士研究生，研究方向：車輛工程。endprint