基于相關(guān)系數(shù)法的柴油機(jī)動(dòng)態(tài)特征提取

冀樹(shù)德 張興剛 孫學(xué)東 閆瑞琦 高海濤

（中國(guó)北方發(fā)動(dòng)機(jī)研究所）

基于相關(guān)系數(shù)法的柴油機(jī)動(dòng)態(tài)特征提取

冀樹(shù)德 張興剛 孫學(xué)東 閆瑞琦 高海濤

（中國(guó)北方發(fā)動(dòng)機(jī)研究所）

以某12缸柴油機(jī)裝車后行駛的砂石路面為基礎(chǔ)，針對(duì)不同駕駛員在車輛實(shí)際行駛中進(jìn)行的7次加速過(guò)程，利用高速采集設(shè)備同步采集了柴油機(jī)輸出端的轉(zhuǎn)速和扭矩信號(hào)，并以動(dòng)作響應(yīng)快的轉(zhuǎn)速信號(hào)為主要參數(shù)設(shè)計(jì)了典型的特征波形。利用相關(guān)系數(shù)法分析了7次加速信號(hào)及不同特征波形信號(hào)間的相關(guān)性，獲得了砂石路行駛條件下柴油機(jī)動(dòng)態(tài)特征信號(hào)。所提取的特征波形信號(hào)圖形可表明，采用相關(guān)系數(shù)法分析典型動(dòng)態(tài)特征信號(hào)是可行的。

1 前言

在柴油機(jī)的設(shè)計(jì)研發(fā)過(guò)程中，其性能驗(yàn)證主要以穩(wěn)態(tài)工況為主，而柴油機(jī)裝車后的實(shí)際使用過(guò)程中60%以上是動(dòng)態(tài)運(yùn)行工況，因此，基于柴油機(jī)的實(shí)際運(yùn)行情況，提取可用于臺(tái)架研究的動(dòng)態(tài)運(yùn)轉(zhuǎn)特征，摸索柴油機(jī)的動(dòng)態(tài)性能，獲取重要的試驗(yàn)數(shù)據(jù)和資料，是柴油機(jī)產(chǎn)品設(shè)計(jì)及研發(fā)環(huán)節(jié)的重要補(bǔ)充。

對(duì)于柴油機(jī)動(dòng)態(tài)特性的研究，目前廣泛開(kāi)展的是以排放研究為主的瞬態(tài)試驗(yàn)循環(huán)[1，4～9]。對(duì)于柴油機(jī)典型動(dòng)態(tài)特征的提取，主要是試驗(yàn)為主的統(tǒng)計(jì)方法研究。相關(guān)系數(shù)法作為統(tǒng)計(jì)學(xué)的一種重要方法，廣泛應(yīng)用于實(shí)測(cè)響應(yīng)空間的負(fù)荷動(dòng)態(tài)特征分類提取、多傳感器信號(hào)動(dòng)態(tài)特征分析等方面[2，3，5]。

本文以某柴油機(jī)裝車后行駛的砂石路面為基礎(chǔ)，針對(duì)不同駕駛員在車輛實(shí)際行駛中進(jìn)行的7次同樣方式的加速過(guò)程，利用相關(guān)系數(shù)法分析、提取柴油機(jī)運(yùn)轉(zhuǎn)中的典型動(dòng)態(tài)工況，以研究柴油機(jī)動(dòng)態(tài)性能。

2 相關(guān)系數(shù)法的基本原理

相關(guān)系數(shù)反映了變量間相關(guān)關(guān)系的密切程度。相關(guān)系數(shù)按積差方法計(jì)算，以兩變量與各自平均值的離差為基礎(chǔ)，2個(gè)離差相乘來(lái)反映兩變量間的相關(guān)程度。通常相關(guān)系數(shù)法用以度量不同變量間存在相關(guān)關(guān)系的可能性。本文以相同變量——柴油機(jī)轉(zhuǎn)速為基礎(chǔ)，分析不同時(shí)間段轉(zhuǎn)速歷程的相似性，以獲得柴油機(jī)動(dòng)態(tài)運(yùn)行條件下的典型運(yùn)行工況。

相關(guān)系數(shù)的具體算法如下:

式中:x、y為變量；n為變量個(gè)數(shù)；σx和σy分別為變量x和y的方差；σxy為變量x和y的協(xié)方差；γ為相關(guān)系數(shù)。

3 柴油機(jī)參數(shù)及載荷數(shù)據(jù)獲取

某柴油機(jī)（主要技術(shù)參數(shù)見(jiàn)表1）裝車后以常用的砂石路面為基礎(chǔ)，由不同駕駛員駕駛車輛在此路面連續(xù)進(jìn)行7次相同方式的加速，同時(shí)采用高速采集設(shè)備實(shí)時(shí)采集、記錄柴油機(jī)輸出端的轉(zhuǎn)速和扭矩，構(gòu)成柴油機(jī)的動(dòng)態(tài)載荷數(shù)據(jù)，采樣頻率為500 Hz。

表1 柴油機(jī)技術(shù)參數(shù)

4 典型特征波形的設(shè)計(jì)

柴油機(jī)在砂石路面連續(xù)7次加速獲得的道路載荷數(shù)據(jù)見(jiàn)圖1。由圖1可看出，與扭矩信號(hào)相比，轉(zhuǎn)速信號(hào)對(duì)于駕駛操作變化、路面工況變化等有更快的響應(yīng)特性。因此，以轉(zhuǎn)速信號(hào)為基礎(chǔ)進(jìn)行柴油機(jī)基于砂石路面的運(yùn)行特征研究，通過(guò)比較和分析7次道路載荷數(shù)據(jù)，并結(jié)合路面和既定的駕駛員操作，設(shè)計(jì)圖2的波形為典型特征波形。

5 道路載荷特征提取

為進(jìn)一步分析7次加速信號(hào)間的相似性及共有特征存在的可能性，對(duì)7次加速信號(hào)進(jìn)行相關(guān)系數(shù)法分析，其結(jié)果見(jiàn)表2。由表2可知，第5次加速信號(hào)與其它信號(hào)間的最大相關(guān)系數(shù)絕對(duì)值僅為0.23，幾乎接近于0，由此推測(cè)該加速信號(hào)與其他信號(hào)間基本不存在相關(guān)關(guān)系。由于選取信號(hào)的變量相同，只是時(shí)間段不同，可進(jìn)一步推測(cè)該信號(hào)與其它信號(hào)間幾乎不存在相似性。在路面特征、車輛結(jié)構(gòu)和柴油機(jī)未發(fā)生變化的情況下，駕駛習(xí)慣的變化是導(dǎo)致其與其它信號(hào)相關(guān)關(guān)系較小的最可能原因。在所有7次加速中，第5次加速信號(hào)與其它信號(hào)間的相關(guān)系數(shù)都較小（表2中陰影部分）。為分析路面固有特征引起柴油機(jī)運(yùn)行工況或特性的變化情況，在此項(xiàng)研究中忽略第5次加速信號(hào)。在剩余的信號(hào)中，某一信號(hào)與其余信號(hào)之間的相關(guān)關(guān)系系數(shù)同時(shí)存在較小可忽略、較大不可忽略的情況，為保持分析信息的完整性，不再進(jìn)行信號(hào)的剔除處理。

表2 不同加速信號(hào)間的相關(guān)系數(shù)

從去除第5次外的其它加速信號(hào)中提取典型特征波形信號(hào)，如圖3所示。針對(duì)每種特征的波形信號(hào)進(jìn)行不同加速信號(hào)間的相關(guān)性分析，結(jié)果見(jiàn)表3～表6。由表3可知，第1種特征的第1次加速信號(hào)與第6次、第7次加速信號(hào)間的相關(guān)性較??；由表4可知，第2種特征的第4次加速信號(hào)與第3次、第6次加速信號(hào)間的相關(guān)系數(shù)較小，與第7次加速信號(hào)間的相關(guān)系數(shù)較小或負(fù)相關(guān)；由表5可知，第3種特征的第7次加速信號(hào)與其它5次加速信號(hào)間主要以負(fù)相關(guān)性為主，其值也較小；由表6可知，第4種特征的第6次加速信號(hào)與其它5次加速信號(hào)間相關(guān)性較小。表3～表6中，陰影部分表示加速信號(hào)間相關(guān)系數(shù)大于0.50的部分，即路面特征變化引起的2次加速行為具有一定的繼承性或相似性?；谝陨戏治黾皥D3的特征圖形可知，相關(guān)系數(shù)分析結(jié)果與實(shí)際圖形是相吻合的。

表3 第1種特征不同加速信號(hào)間的相關(guān)性

表4 第2種特征不同加速信號(hào)間的相關(guān)性

表5 砂石路第3種特征不同加速信號(hào)間的相關(guān)性

表6 第4種特征不同加速信號(hào)間的相關(guān)性

為了獲得特征信號(hào)間的相似性或共有特征，舍去信號(hào)中不相關(guān)或相關(guān)性小的信號(hào)，將特征中剩余部分信號(hào)計(jì)算平均值，同步扭矩的特征即構(gòu)成柴油機(jī)砂石路面條件下加速狀態(tài)的運(yùn)行特征信號(hào)，即動(dòng)態(tài)特征信號(hào)，見(jiàn)圖4。

6 結(jié)束語(yǔ)

采用相關(guān)系數(shù)法分析了車輛7次加速過(guò)程中柴油機(jī)的運(yùn)行特征，比較了加速信號(hào)之間相似性和差異性。通過(guò)7次加速信號(hào)的相關(guān)系法分析，最大程度地保留了加速信號(hào)的共有特征，也表征了柴油機(jī)實(shí)際應(yīng)用過(guò)程中的典型動(dòng)態(tài)特征，將其應(yīng)用于基于臺(tái)架試驗(yàn)驗(yàn)證的設(shè)計(jì)開(kāi)發(fā)，可拓展、完善產(chǎn)品設(shè)計(jì)驗(yàn)證方法，提高產(chǎn)品成品率。

1Constantine D.Rakopoulos，Evangelos Gikoumis.Diesel engine transient operation.London:Springer，2009.

2李欣然，林舜江，劉楊華，等.基于實(shí)測(cè)響應(yīng)空間的負(fù)荷特性分類原理與方法.中國(guó)電機(jī)工程學(xué)報(bào)，2006，26（8）:39～43.

3趙鳴，柳瑩.結(jié)構(gòu)動(dòng)態(tài)響應(yīng)多傳感器監(jiān)測(cè)信號(hào)數(shù)據(jù)融合分析.同濟(jì)大學(xué)學(xué)報(bào)（自然科學(xué)版），2008，36（5）:574-578.

4冀樹(shù)德，許世永，張偉，等.柴油機(jī)瞬態(tài)試驗(yàn)循環(huán)的線性回歸分析.車用發(fā)動(dòng)機(jī)，2010（5）:52～55.

5劉波瀾，黃英，張付軍，等.車用柴油機(jī)加速瞬態(tài)過(guò)程控制策略研究.車用發(fā)動(dòng)機(jī)，2008（6）:8～19.

6王忠恕，吳楠，許允，等.增壓直噴柴油機(jī)瞬態(tài)工況燃燒參數(shù)的變化規(guī)律.內(nèi)燃機(jī)學(xué)報(bào)，2007（5）:385～389.

7郭紅霞.相關(guān)系數(shù)及其應(yīng)用.武警工程學(xué)院學(xué)報(bào)，2010（5）: 18～21.

8吳建國(guó)，周鋐，陳棟華.目標(biāo)用戶道路譜與試驗(yàn)場(chǎng)道路譜的載荷當(dāng)量等效模擬研究.汽車技術(shù)，2007（7）:21～24.

9劉再生，霍福祥，王長(zhǎng)明，等.轎車懸架臺(tái)架多軸疲勞試驗(yàn)載荷開(kāi)發(fā).汽車技術(shù)，2012（06）:47～50.

（責(zé)任編輯文楫）

修改稿收到日期為年月日

Extraction of Dynamic Characteristic for Diesel Engine Based on Correlation Coefficient Method

Ji Shude，Zhang Xinggang，Sun Xuedong，Yan Ruiqi，Gao Haitao
（China North Engine Research Institute）

For the regularly running sandstone road of vehicle powered by a 12-cylinder diesel engine，and seven accelerations by different drivers at actual driving，we collect the output speed and torque signals of diesel engine with the high-speed acquisition unit，and use the speed signal which has rapid response as main parameter to design the typical characteristic waveform.The correlation between acceleration signals and different characteristic wave signals is analyzed with the correlation coefficient method，and finally the dynamic characteristic signal in sandstone road driving is acquired. The extracted characteristic wave signal graph shows that the correlation coefficient method is feasible to analyze the typical dynamic characteristic signal.

Diesel engine，Dynamic performance，Road load;Correlation coefficient method

柴油機(jī)動(dòng)態(tài)特性道路載荷相關(guān)系數(shù)法

U464.12+1

A

1000-3703（2014）12-0018-04