基于貝葉斯網(wǎng)絡(luò)集成的船用中高速發(fā)動機(jī)磨損故障診斷模型

王永堅， 戴樂陽， 宋佳聲

(集美大學(xué) 輪機(jī)工程學(xué)院，福建 廈門 361021)

船用中高速發(fā)動機(jī)運行轉(zhuǎn)速高、使用條件惡劣，運動副磨損故障是其最常出現(xiàn)的故障之一。磨損故障的判別和診斷主要通過監(jiān)測與發(fā)動機(jī)磨損緊密相關(guān)的各類參數(shù)和檢測曲軸箱在用潤滑油理化/鐵譜/光譜等性能參數(shù)和圖像得到的原始數(shù)據(jù)，獲取相關(guān)的故障信息，并從這些信息中提取故障征狀進(jìn)行相應(yīng)處理，最后根據(jù)故障征狀和其他診斷信息定位并隔離故障，形成故障診斷與維修決策來實現(xiàn)，但磨損故障特征存在模糊性和不完備性等問題。[1]為解決該問題，可搭建基于專家知識與貝葉斯網(wǎng)絡(luò)(Bayesian Networks)集成的船用中高速發(fā)動機(jī)磨損故障診斷模型，在該模型中輸入相關(guān)監(jiān)測參數(shù)信息，及時對運行中的發(fā)動機(jī)運動副的磨損狀況及其技術(shù)狀態(tài)進(jìn)行監(jiān)控，判斷其故障磨損部位和磨損程度，給出故障預(yù)警和指導(dǎo)意見，評判發(fā)動機(jī)的“健康”狀態(tài)，確保發(fā)動機(jī)可靠、安全運行。

1 船用中高速發(fā)動機(jī)運動副摩擦學(xué)系統(tǒng)及其與監(jiān)測技術(shù)的關(guān)系

船用中高速發(fā)動機(jī)是由許多部件和運動副組成的復(fù)雜、相互關(guān)聯(lián)的結(jié)構(gòu)系統(tǒng)，不同類型的摩擦副是發(fā)動機(jī)運動副的物理實現(xiàn)。摩擦副及其潤滑控制、狀態(tài)監(jiān)測和補(bǔ)償控制系統(tǒng)組成船用發(fā)動機(jī)摩擦學(xué)系統(tǒng)?；钊h(huán)側(cè)面與缸套內(nèi)表面、連桿大端銷與大端軸承、主軸頸與主軸承及傳動齒輪組等是發(fā)動機(jī)重要的摩擦運動副。[2]圖1為船用中高速發(fā)動機(jī)運動副摩擦學(xué)物理模型。在該模型中，不同類型的摩擦運動副的輸入和輸出都伴隨著摩擦學(xué)行為，發(fā)動機(jī)摩擦學(xué)系統(tǒng)的構(gòu)成和運行機(jī)理賦予系統(tǒng)顯著的系統(tǒng)依賴性和時間依賴性[3]：

1) 系統(tǒng)的特性是所有摩擦運動副及與其關(guān)聯(lián)的工作介質(zhì)或環(huán)境的綜合特性。

2) 系統(tǒng)的特性與發(fā)動機(jī)的運行工況及其所處的運行環(huán)境密切相關(guān)。隨著發(fā)動機(jī)運轉(zhuǎn)時間的不斷增加，系統(tǒng)中構(gòu)成元素的狀態(tài)也在不斷變化。

在運動副的磨損監(jiān)測中，磨損故障與曲軸箱在用潤滑油監(jiān)測參數(shù)、運行時熱工參數(shù)和振動等有著密切關(guān)系。通過監(jiān)測這些參數(shù)，可獲取與發(fā)動機(jī)磨損故障相關(guān)的狀態(tài)特征信息，這些信息不能通過簡單的疊加產(chǎn)生更充分的故障信息，而需對各種監(jiān)測數(shù)據(jù)進(jìn)行分類、融合，并與已有的監(jiān)測診斷知識(經(jīng)驗)相關(guān)聯(lián)，進(jìn)行推理和匹配，只有如此才能給出發(fā)動機(jī)摩擦學(xué)狀態(tài)的判別結(jié)論，即通過集成各種監(jiān)測技術(shù)和融合各種特征參數(shù)來提高發(fā)動機(jī)摩擦學(xué)故障診斷精度。[4]圖2為船用中高速發(fā)動機(jī)缸套活塞組件常見的磨損故障及其與曲軸箱在用潤滑油、熱工參數(shù)和振動等監(jiān)測技術(shù)的關(guān)系圖。

2 貝葉斯網(wǎng)絡(luò)

貝葉斯網(wǎng)絡(luò)又稱信念網(wǎng)絡(luò)，是一種基于網(wǎng)絡(luò)結(jié)構(gòu)的有向無環(huán)圖解概率模型(DAG模型)，是人工智能、概率理論、圖論和決策理論相結(jié)合的產(chǎn)物。[5]

由一組變量X=(X1,X2,…,Xn)構(gòu)成的貝葉斯網(wǎng)絡(luò)用二元組B(G，P)表示，其中：G為有向無環(huán)圖；P為條件概率表。圖3為一個簡單的貝葉斯網(wǎng)絡(luò)結(jié)構(gòu)圖(省略各節(jié)點的條件概率)。

貝葉斯網(wǎng)絡(luò)能將與故障診斷有關(guān)的各種信息(如故障征狀、故障模式及測試值等)用節(jié)點的形式表示，節(jié)點之間的關(guān)聯(lián)關(guān)系用網(wǎng)絡(luò)中的有向邊表示，其緊密程度用節(jié)點間的條件概率(CPT)表達(dá)。[6]將貝葉斯網(wǎng)絡(luò)及其強(qiáng)大的學(xué)習(xí)和統(tǒng)計推斷功能應(yīng)用于設(shè)備故障診斷中，不僅能有效處理故障診斷過程中的不確定問題，還能對多源異類信息進(jìn)行有效表達(dá)和融合[7]，在設(shè)備故障診斷領(lǐng)域具有較高的推廣應(yīng)用價值。

3 貝葉斯診斷模型的搭建

3.1 發(fā)動機(jī)磨損故障貝葉斯診斷網(wǎng)絡(luò)結(jié)構(gòu)的搭建

船用中高速發(fā)動機(jī)磨損故障診斷一般從發(fā)動機(jī)運動副發(fā)生磨損故障時出現(xiàn)與之關(guān)聯(lián)的典型故障征狀(如曲軸箱在用潤滑油監(jiān)測參數(shù)、運行時熱工參數(shù)、振動與噪聲及發(fā)動機(jī)性能指標(biāo)的變化)出發(fā)，挖掘運動副出現(xiàn)磨損故障時呈現(xiàn)的故障征狀與故障模式或原因(如缸套擦傷等)之間的因果關(guān)系，最終獲得故障原因和故障位置的推理。表1為船用中高速發(fā)動機(jī)磨損故障經(jīng)多名專家確認(rèn)之后的層次結(jié)構(gòu)。表1中，發(fā)動機(jī)的磨損故障診斷過程分為3個層次，即與發(fā)動機(jī)磨損故障有關(guān)的典型征狀(主要是定期采集故障監(jiān)測參數(shù))、與磨損故障征狀密切相關(guān)的故障模式和由故障模式引發(fā)的故障運動副。由于不同運動副在進(jìn)行磨損故障診斷時存在復(fù)雜性、監(jiān)測手段局限性和知識表達(dá)的不確定性等問題，均會導(dǎo)致故障征狀與故障模式之間的因果映射出現(xiàn)隨機(jī)性和不確定性；加上在磨損故障診斷中能獲得的觀測數(shù)據(jù)(特別是故障監(jiān)測樣本)非常有限，甚至是不完整、不準(zhǔn)確的，船用發(fā)動機(jī)運動副只有在出現(xiàn)嚴(yán)重磨損故障影響其正常運行時，與之關(guān)聯(lián)的典型故障征狀(如曲軸箱在用潤滑油監(jiān)測參數(shù)等)才會被監(jiān)測，致使故障樣本非常少。因此，在進(jìn)行磨損故障診斷時，通常結(jié)合專家豐富的知識(經(jīng)驗)及對發(fā)動機(jī)運動副結(jié)構(gòu)、所用材料等內(nèi)容的分析研究，通過專家會診和問詢等方式進(jìn)行故障原因分析及判別。采用提取專家知識，利用專家的因果映射知識的方式搭建發(fā)動機(jī)磨損故障貝葉斯診斷網(wǎng)絡(luò)結(jié)構(gòu)，具體步驟如下。

表1 船用中高速發(fā)動機(jī)磨損故障層次結(jié)構(gòu)

1)數(shù)據(jù)提?。和ㄟ^定性采訪的方式向?qū)＜姨岢鲆幌盗杏沙橄蟮骄唧w、由淺入深的與發(fā)動機(jī)磨損故障相關(guān)的探測性問題，以獲得盡可能多地涉及貝葉斯網(wǎng)絡(luò)目標(biāo)域的數(shù)據(jù)信息，并以描述性文本的形式存儲這些信息。例如：當(dāng)發(fā)動機(jī)缸套活塞組件過度磨損時，曲軸箱在用潤滑油性能指標(biāo)、發(fā)動機(jī)熱工參數(shù)和振動等會表現(xiàn)出哪些故障征狀；當(dāng)主軸瓦表面擦傷時，在用潤滑油理化/鐵譜/光譜監(jiān)測參數(shù)、發(fā)動機(jī)滑油壓力和振動與噪聲等會出現(xiàn)哪些異常變化等。類似的問題通過詢問專家，獲取與磨損故障征狀有關(guān)的描述性文本信息。

2)因果映射提取與建立：首先針對步驟1)中形成的專家知識描述性文本，識別其中的因果闡述；其次將每條因果闡述拆分為原因項、因果連接詞和結(jié)果項等3部分，形成初始的因果映射文本；最后經(jīng)必要的提煉和編譯，將專家針對發(fā)動機(jī)磨損故障診斷的因果映射轉(zhuǎn)化為編譯之后的認(rèn)知映射。

3)轉(zhuǎn)化因果映射為貝葉斯網(wǎng)絡(luò)結(jié)構(gòu)：通過分析船用中高速發(fā)動機(jī)故障征狀節(jié)點、故障模式節(jié)點和故障運動副節(jié)點間(如表1所示)的條件獨立性，區(qū)分節(jié)點間的直接因果關(guān)系和間接因果關(guān)系，挖掘節(jié)點間潛在的因果關(guān)系，并在排除循環(huán)因果關(guān)系之后轉(zhuǎn)化因果映射網(wǎng)絡(luò)為磨損故障的貝葉斯診斷網(wǎng)絡(luò)結(jié)構(gòu)?；谏鲜龇治鰴C(jī)理，結(jié)合表1中各節(jié)點的設(shè)置，搭建船用中高速發(fā)動機(jī)磨損故障貝葉斯診斷網(wǎng)絡(luò)結(jié)構(gòu)圖(見圖4)。

3.2 貝葉斯診斷網(wǎng)絡(luò)參數(shù)的學(xué)習(xí)

船用中高速發(fā)動機(jī)磨損故障貝葉斯診斷模型的搭建分為模型網(wǎng)絡(luò)結(jié)構(gòu)的構(gòu)建、網(wǎng)絡(luò)參數(shù)的學(xué)習(xí)和確定及推理算法等3部分。在搭建基于專家知識和因果映射的發(fā)動機(jī)磨損故障貝葉斯診斷網(wǎng)絡(luò)結(jié)構(gòu)之后，還需確定網(wǎng)絡(luò)結(jié)構(gòu)中各節(jié)點的先驗概率和條件概率(CPT)等數(shù)據(jù)信息。[8]由于依據(jù)專家知識(經(jīng)驗)搭建的磨損診斷網(wǎng)絡(luò)結(jié)構(gòu)所需的與發(fā)動機(jī)磨損故障有關(guān)的各種知識、經(jīng)驗及實際檢測數(shù)據(jù)、實驗室存儲的歷史數(shù)據(jù)較為豐富，故該診斷網(wǎng)絡(luò)結(jié)構(gòu)與實際情況較為接近。[9]但是，網(wǎng)絡(luò)結(jié)構(gòu)中各節(jié)點的條件概率是診斷模型的定量信息，單靠專家知識等進(jìn)行概率確定，會與實際故障情況存在一定的偏差，會對后驗概率的準(zhǔn)確性、故障模式的判別產(chǎn)生不利影響。因此，在發(fā)動機(jī)的磨損故障診斷分析中，貝葉斯診斷網(wǎng)絡(luò)參數(shù)的學(xué)習(xí)更為重要。在發(fā)動機(jī)磨損故障貝葉斯診斷網(wǎng)絡(luò)結(jié)構(gòu)中，各故障征狀節(jié)點條件概率的學(xué)習(xí)和確定采用以下2種方式：

1) 在缺乏足夠多的故障訓(xùn)練樣本的情況下，通過詢問專家獲取相關(guān)數(shù)據(jù)信息。為便于專家進(jìn)行概率估計，采用概率刻度法進(jìn)行賦值(見圖5)。圖5中，確定、很有可能、較為可能、50%可能、不太確定、不太可能和不可能對應(yīng)的概率依次為1.00、0.85、0.75、0.50、0.25、0.15和0。

2) 采用完整數(shù)據(jù)學(xué)習(xí)的貝葉斯方法(Bayes_update_params)進(jìn)行參數(shù)學(xué)習(xí)。貝葉斯方法參數(shù)學(xué)習(xí)由觀測前的先驗知識和觀測到的樣本數(shù)據(jù)2部分組成。[7]先驗知識包括參數(shù)先驗分布的選取和分布參數(shù)的選取規(guī)則。學(xué)習(xí)的目標(biāo)是根據(jù)概率分布的先驗分布和觀測樣本數(shù)據(jù)計算概率分布的后驗概率分布。

設(shè)概率分布向量θij=(θij,θij,…,θij)，在θij相互獨立的情況下，若無數(shù)據(jù)缺失，則后驗概率分布保持獨立，且有

(1)

式(1)中：D為學(xué)習(xí)樣本集；S為網(wǎng)絡(luò)結(jié)構(gòu)模型；θS為概率分布向量；取Dirichlet分布Dir(θij|αij1,αij2,…,αijri)作為先驗分布，結(jié)合磨損故障征狀樣本數(shù)據(jù)，可得后驗概率分布為

(2)

至此，可獲得貝葉斯網(wǎng)絡(luò)的后驗(條件)概率參數(shù)，即先驗信息和故障樣本數(shù)據(jù)信息的綜合?；谏鲜鼍W(wǎng)絡(luò)參數(shù)的學(xué)習(xí)機(jī)理，船用中高速發(fā)動機(jī)磨損故障診斷網(wǎng)絡(luò)后驗概率參數(shù)學(xué)習(xí)可由先驗概率分布和各監(jiān)測參數(shù)樣本數(shù)據(jù)來確定。

3.3 貝葉斯診斷網(wǎng)絡(luò)模型的推理算法

貝葉斯網(wǎng)絡(luò)推理算法分為精確推理算法和近視推理算法。[10]根據(jù)實際需求，采用精確推理算法中的聯(lián)合樹推理算法(engine=jtree_inf_engine(bnet))。該算法適用于網(wǎng)絡(luò)節(jié)點少、結(jié)構(gòu)簡單的網(wǎng)絡(luò)結(jié)構(gòu)，能滿足大多數(shù)貝葉斯網(wǎng)絡(luò)的推理。船用中高速發(fā)動機(jī)磨損故障貝葉斯診斷模型中各故障模式節(jié)點后驗概率的獲取及推理過程均在MATLAB軟件的仿真環(huán)境中采用FullBNT-1.07平臺推理演算。圖6為船用中高速發(fā)動機(jī)磨損故障貝葉斯診斷模型的搭建和診斷流程圖，故障診斷和推理運算過程為：

1)定義發(fā)動機(jī)磨損故障3種節(jié)點(故障征狀節(jié)點、故障模式節(jié)點和故障運動副節(jié)點)，根據(jù)專家知識構(gòu)建貝葉斯診斷網(wǎng)絡(luò)模型；在此基礎(chǔ)上，采用詢問專家與貝葉斯方法相結(jié)合的方法確認(rèn)各故障模式節(jié)點的條件概率表(CPT)。

2)在FullBNT-1.07平臺上輸入已確認(rèn)的某次磨損故障征狀節(jié)點信息，運用平臺內(nèi)嵌的聯(lián)合樹推理引擎進(jìn)行推理運算，獲得各故障模式節(jié)點后驗(發(fā)生)概率。

3)選擇發(fā)生故障概率最大的運動副進(jìn)行維修，若故障得以排除，則診斷結(jié)束；若故障未能排除，則該運動副節(jié)點狀態(tài)為0(正常)，作為新證據(jù)融入系統(tǒng)重新進(jìn)行推理運算；以此類推，直至查找到故障問題所在，診斷結(jié)束。

4 發(fā)動機(jī)磨損故障問題的仿真和分析

基于已搭建的船用中高速發(fā)動機(jī)磨損故障貝葉斯診斷模型，以發(fā)動機(jī)具體磨損故障問題為研究對象，進(jìn)行實際案例仿真和分析，根據(jù)仿真結(jié)果驗證所搭建診斷模型的準(zhǔn)確性和分析方法的有效性。為便于進(jìn)行仿真計算，根據(jù)專家知識(經(jīng)驗)和采集的3家航運企業(yè)所屬近40艘船舶在用中高速柴油發(fā)動機(jī)(約100臺)3 a中發(fā)生磨損故障時的各種現(xiàn)場故障征狀信息，結(jié)合存儲的大量與發(fā)動機(jī)磨損故障有關(guān)的歷史數(shù)據(jù)，對仿真測試所需的各磨損故障征狀節(jié)點的先驗概率進(jìn)行統(tǒng)計，獲取相關(guān)概率信息(見表2)。

表2 發(fā)動機(jī)磨損故障征狀先驗概率

4.1 理論分析

以母型船上一臺受監(jiān)控的發(fā)電柴油原動機(jī)(濰柴6170，轉(zhuǎn)速1 000 r/min，額定功率300 kW，持續(xù)運行時間為5 882 h)為監(jiān)測對象，通過求取該發(fā)動機(jī)送檢的曲軸箱在用潤滑油理化/鐵譜/光譜檢測數(shù)據(jù)(采樣時間為5 882 h)及根據(jù)該采樣時間獲得柴油機(jī)熱工參數(shù)、性能指標(biāo)、振動等數(shù)據(jù)，并與正常使用范圍相對比，獲得磨損故障征狀信息(見表3)。從表3中可看出送檢潤滑油分析鐵譜片狀磨粒較多、直讀鐵譜Fe大小磨粒濃度較高、光譜測量B元素濃度較高、柴油機(jī)平均壓縮壓力/爆炸壓力均降低和柴油機(jī)振動較為明顯等故障征狀。針對這些故障征狀信息，結(jié)合表1、圖4和圖6對該發(fā)動機(jī)磨損故障進(jìn)行如下診斷和分析。

表3 機(jī)器狀況得分評定表

根據(jù)該發(fā)動機(jī)監(jiān)測到的與故障磨損緊密相關(guān)的典型故障征狀，結(jié)合已搭建的磨損故障貝葉斯診斷模型，獲取的故障征狀特征節(jié)點為片狀磨粒較多、B元素濃度較高、Fe大小磨粒濃度較高、缸內(nèi)壓縮/爆炸壓力降低和振動較為明顯，這些故障癥狀對應(yīng)的節(jié)點為S2、S4、S7、S9和S11。由圖4可知：此時S2、S4、S7和S9節(jié)點對應(yīng)缸套活塞組磨損過快故障模式，即節(jié)點F12發(fā)生；S4、S7和S11節(jié)點對應(yīng)缸套表面擦傷故障模式，即節(jié)點F11發(fā)生；S9和S11節(jié)點對應(yīng)大端軸瓦表面材料剝落故障模式，即節(jié)點F32發(fā)生。通過理論分析可知，出現(xiàn)F11和F12的可能性較大，且F12發(fā)生的可能性大一些。

4.2 仿真驗證

在MATLAB環(huán)境中，利用FullBNT-1.07平臺進(jìn)行受控柴油發(fā)動機(jī)在5 882 h時間段內(nèi)磨損故障的故障模式發(fā)生(后驗)概率推理運算，運算過程中需獲取圖1中各故障模式(F11～F41)的條件概率，采用貝葉斯參數(shù)學(xué)習(xí)與詢問專家相結(jié)合的方式獲取。表4為當(dāng)故障征狀節(jié)點S6(Pb元素濃度高)和S8(Cu元素濃度高)狀態(tài)不同(0故障征狀不發(fā)生；1故障發(fā)生)時，F(xiàn)22主軸瓦表面材料剝落的條件概率。同理，可計算出其余的故障模式節(jié)點的條件概率，限于篇幅，不再一一列出。

表4 F22的條件概率

針對受控柴油發(fā)動機(jī)已確認(rèn)的磨損征狀信息，基于已搭建的貝葉斯診斷模型、各征狀節(jié)點的先驗概率信息(如表2所示)和各故障模式的條件概率信息(如表4所示)，對受控柴油發(fā)動機(jī)磨損故障進(jìn)行仿真驗證。在FullBNT-1.07推理運算平臺上進(jìn)行仿真測試，具體實現(xiàn)過程及部分源代碼為：

N=19, dag = zeros(N,N); //定義貝葉斯DAG結(jié)構(gòu)圖

S1=1;…;S12=12;F11=13; …;F41=19;//19個網(wǎng)絡(luò)節(jié)點

dag(S1,F11)=1; //節(jié)點之間的連接關(guān)系

……

dag(S4, [F21,F22])=1;

……

discrete_nodes = 1:N; //離散節(jié)點

node_sizes = 2*ones(1,N); //節(jié)點狀態(tài)數(shù)

//手工構(gòu)造條件概率CPT表

bnet.CPD{S1}=tabular_CPD(bnet,S1, [0.03 0.97]);

bnet.CPD{S2}=tabular_CPD(bnet,S2, [0.07 0.93]);

……

bnet.CPD{F22}=tabular_CPD(bnet,F22, [0.75 0.25 0.45 0.01 0.25 0.75 0.55 0.99]);

……

bnet=mk_bnet(dag, node_sizes,'names',{'S1','S2',...'F11',…'F41'},'discrete',discrete_nodes);

figure; draw_graph(dag);

engine = jtree_inf_engine(bnet); //選擇jtree_inf_engine推理引擎

evidence = cell(1,N);

evidence{S2} =1; //輸入證據(jù)

evidence{S4} =1

……

[engine, loglike]= enter_evidence(engine, evidence);

Marg1 = marginal_nodes(engine,[F11F21F21F22F31F32F41]);

Marg1.F; //計算故障模式節(jié)點后驗概率。

仿真過程包括沒有故障征狀證據(jù)節(jié)點信息輸入和出現(xiàn)故障征狀節(jié)點信息輸入2個。當(dāng)S2、S4、S7、S9和S11節(jié)點的取值狀態(tài)為1(出現(xiàn)故障征狀)時，將這些征狀信息(evidence{S2,…,S11} =1)輸入貝葉斯診斷網(wǎng)絡(luò)模型中進(jìn)行證據(jù)推理傳播，可獲得故障證據(jù)傳播之后各故障模式節(jié)點的發(fā)生(后驗)概率，根據(jù)故障概率的高低確定磨損部位和部件，決定維修的順序。

4.3 仿真結(jié)果分析

圖7為沒有故障征狀節(jié)點信息(evidence{Si}=0故障征兆未出現(xiàn))輸入時發(fā)動機(jī)各磨損故障模式的直方圖，故障模式節(jié)點F11、F12、F21、F22、F31、F32和F41的發(fā)生概率分別為0.242 6、0.283 4、0.225 6、0.201 3、0.231 8、0.244 5和0.189 2。當(dāng)發(fā)動機(jī)磨損故障征狀信息輸入(即S2、S4、S7、S9和S11節(jié)點的取值狀態(tài))為1時，經(jīng)仿真運算，獲得各故障模式節(jié)點的概率值見表5，對應(yīng)的故障模式的直方圖見圖8。

分析表5和圖8可得，此時母型船受監(jiān)控柴油發(fā)動機(jī)磨損故障模式發(fā)生概率和維修檢查順序為：F12>F11>F32>F31>F41>F21>F22，F(xiàn)12(缸套活塞組磨損過快)節(jié)點和F11(缸套表面擦傷)節(jié)點發(fā)生磨損故障的概率大，且F12節(jié)點發(fā)生的概率最大，仿真結(jié)果與理論分析結(jié)果一致，驗證了所搭建的磨損故障貝葉斯診斷網(wǎng)絡(luò)模型的準(zhǔn)確性。仿真運算結(jié)果顯示，由于S9和S11故障征狀信息的出現(xiàn)，與其相關(guān)的故障模式節(jié)點F32的發(fā)生概率也比較大，此外分別與S7和S11故障征狀關(guān)聯(lián)的故障模式節(jié)點F41和F31發(fā)生故障的概率也偏大，且發(fā)生故障的概率非常接近，與理論分析結(jié)果一致，說明了所采用的磨損故障診斷分析方法的有效性。

表5 故障征狀節(jié)點輸入后各故障模式的發(fā)生概率

5 結(jié)束語

為克服船用中高速發(fā)動機(jī)運動副發(fā)生磨損故障時出現(xiàn)的故障特征不完備和模糊等問題，通過定期采集與發(fā)動機(jī)磨損故障緊密相關(guān)的曲軸箱在用潤滑油監(jiān)測參數(shù)、運行參數(shù)和熱工參數(shù)等信息獲取原始數(shù)據(jù)，結(jié)合專家知識(經(jīng)驗)搭建磨損故障貝葉斯診斷網(wǎng)絡(luò)系統(tǒng)結(jié)構(gòu)，并通過詢問專家和運用貝葉斯方法2種方法進(jìn)行網(wǎng)絡(luò)參數(shù)學(xué)習(xí)，完成基于專家知識和貝葉斯網(wǎng)絡(luò)集成的船用中高速發(fā)動機(jī)磨損故障診斷模型的搭建，通過實際故障案例的仿真和分析，驗證模型的準(zhǔn)確性和分析方法的有效性，取得令人滿意的結(jié)果，可為快速、準(zhǔn)確發(fā)現(xiàn)發(fā)動機(jī)磨損故障問題提供強(qiáng)有力的技術(shù)支持，具有較高的實際應(yīng)用價值。

[1] 宮義山, 高媛媛. 基于故障樹和貝葉斯網(wǎng)絡(luò)的故障診斷模型[J]. 沈陽工業(yè)大學(xué)學(xué)報, 2009,31(4):454-457.

[2] 滿一新. 船機(jī)維修技術(shù)[M]. 大連: 大連海事大學(xué)出版社, 1999：86-96.

[3] 楊其明, 嚴(yán)新平,賀石中. 油液監(jiān)測分析現(xiàn)場實用技術(shù)[M].北京: 機(jī)械工業(yè)出版社, 2006：126-130.

[4] 姜萬錄, 劉思遠(yuǎn). 多特征信息融合的貝葉斯網(wǎng)絡(luò)故障診斷方法研究[J]. 中國機(jī)械工程, 2010, 21(8)：940-945.

[5] 劉建偉, 黎海恩, 羅雄麟. 概率圖模型表示理論[J]. 計算機(jī)科學(xué), 2014,41(9):1-17.

[6] KOLLER D, FRIEDMAN N. Probabilistic Graphical Models: Principles and Techniques[M]. MIT Press, 2009: 1-153.

[7] 沈懷榮. 信息融合故障診斷技術(shù)[M]. 北京: 科學(xué)出版社, 2013：264-267.

[8] 劉希亮，陳桂明,李方溪,等. 基于多傳感器信息融合的貝葉斯網(wǎng)絡(luò)故障診斷方法研究及應(yīng)用[J]. 機(jī)械科學(xué)與技術(shù), 2013, 32(1): 91-95.

[9] 林洋.基于貝葉斯網(wǎng)絡(luò)船用柴油機(jī)故障診斷系統(tǒng)的研究[D]. 大連: 大連海事大學(xué), 2016.

[10] 李敏通. 柴油機(jī)振動信號特征提取與故障診斷方法研究[D]. 西安: 西北農(nóng)林科技大學(xué), 2012.