李玲玲 張士暖 李志剛 賀鵬舉
?
基于粗糙集理論和生命初態(tài)信息的繼電器壽命預(yù)測(cè)方法
李玲玲 張士暖 李志剛 賀鵬舉
(河北工業(yè)大學(xué)電磁場(chǎng)與電器可靠性省部共建重點(diǎn)實(shí)驗(yàn)室 天津 300130)
工作壽命是表征產(chǎn)品可靠性的重要特征量。對(duì)于繼電器這類長(zhǎng)壽命產(chǎn)品,以壽命為指標(biāo)進(jìn)行產(chǎn)品篩選的方法目前還未得到深入研究。近來通過對(duì)繼電器壽命試驗(yàn)數(shù)據(jù)進(jìn)行分析發(fā)現(xiàn):繼電器個(gè)體的壽命與其生命初期的性態(tài)有關(guān),初期性態(tài)較好的個(gè)體其壽命也較長(zhǎng),反之則壽命較短。本文基于這一發(fā)現(xiàn)提出了繼電器的個(gè)體壽命預(yù)測(cè)與篩選方法,即首先通過壽命試驗(yàn)獲取一組樣本在生命初期的性能參數(shù)值與壽命值,然后利用粗糙集理論中的屬性約簡(jiǎn)提取一組與壽命相關(guān)的性能參數(shù),并建立一組壽命決策規(guī)則。此后對(duì)于任一同型產(chǎn)品,只要測(cè)出其生命初態(tài)信息即可根據(jù)這組規(guī)則預(yù)測(cè)其壽命,進(jìn)而完成以壽命為指標(biāo)的繼電器產(chǎn)品篩選。本文最后進(jìn)行了實(shí)例分析,結(jié)果表明了該方法的有效性。
繼電器 生命初期特征 壽命預(yù)測(cè) 產(chǎn)品篩選 粗糙集理論
工作壽命是表征工業(yè)產(chǎn)品可靠性的重要特征量。對(duì)于繼電器而言,不同行業(yè)、系統(tǒng)和設(shè)備對(duì)其工作壽命有著不同要求,因此很多時(shí)候需要以壽命為指標(biāo)進(jìn)行繼電器的篩選。例如在車輛系統(tǒng)中,揚(yáng)聲器繼電器要求百萬次左右的工作壽命,前燈繼電器50萬次,而其他一般的繼電器達(dá)3~10萬次即能滿足需要[1]。特別地,當(dāng)繼電器服役于某些對(duì)可靠性要求極高、一旦投入使用即不易修復(fù)甚至根本無法修復(fù)的系統(tǒng)(例如太空站、海上石油平臺(tái)和人造衛(wèi)星等)時(shí),每個(gè)繼電器個(gè)體都被要求擁有比其所服役的系統(tǒng)更長(zhǎng)的工作壽命。為此不得不進(jìn)行歷時(shí)較長(zhǎng)、條件嚴(yán)苛的篩選試驗(yàn),從而占用大量人力物力并產(chǎn)生高昂費(fèi)用,況且這種歷時(shí)較長(zhǎng)的篩選試驗(yàn)本身也會(huì)減少繼電器個(gè)體的剩余壽命,使其有效服役期縮短。隨著我國(guó)太空技術(shù)與軍工技術(shù)等的迅猛發(fā)展,當(dāng)前迫切需要一種以壽命為指標(biāo)對(duì)包括繼電器在內(nèi)的電子電器進(jìn)行篩選的方法,且該方法不會(huì)過多壓縮其有效服役期。
以壽命為指標(biāo)的產(chǎn)品篩選首先需要對(duì)產(chǎn)品壽命做出預(yù)估、預(yù)測(cè)。當(dāng)前,產(chǎn)品壽命預(yù)測(cè)的方法可分為兩類:一類是基于抽樣檢測(cè)理論,根據(jù)樣本的壽命數(shù)據(jù)推測(cè)整體的壽命分布以及與壽命有關(guān)的可靠性特征量(如平均壽命、中位壽命等)。這類方法適用于批產(chǎn)品而非個(gè)體的壽命預(yù)測(cè),對(duì)于以壽命為指標(biāo)的產(chǎn)品篩選意義不大,因?yàn)榧词故峭慌蔚漠a(chǎn)品,不同個(gè)體的壽命也是有差別的,表征整體性能的壽命指標(biāo)并不能充分表達(dá)個(gè)體的情況。另一類方法是根據(jù)產(chǎn)品個(gè)體的性能參數(shù)變化特征推測(cè)其壽命,此類方法建立在產(chǎn)品個(gè)體的性能退化機(jī)理之上,可以體現(xiàn)同一批次下不同產(chǎn)品個(gè)體的壽命差別,故其適用對(duì)象為產(chǎn)品個(gè)體。
在上述第二類方法中,當(dāng)前的研究大多集中于基于神經(jīng)網(wǎng)絡(luò)[2]、回歸分析[3,4]、小波和ARMA分析[5]等數(shù)學(xué)方法建模的預(yù)測(cè)。這類數(shù)值型預(yù)測(cè)方法所需數(shù)據(jù)量較大,計(jì)算過程復(fù)雜。事實(shí)上,在為太空站這類對(duì)可靠性要求極高的系統(tǒng)篩選基礎(chǔ)電子器件(例如繼電器)時(shí),并不需要確切掌握每個(gè)器件的工作壽命,而只要求這些器件的壽命長(zhǎng)于系統(tǒng)自身的服役期,這種情況下,用于器件篩選的壽命預(yù)測(cè)應(yīng)側(cè)重于定性而非定量。另外,較復(fù)雜的自控系統(tǒng)所含繼電器成百上千,例如某型軍用運(yùn)輸機(jī)采用了300多只航天繼電器,某大型人造衛(wèi)星采用的繼電器數(shù)量達(dá)1 500多只[6],如果每個(gè)繼電器都基于數(shù)值型預(yù)測(cè)方法進(jìn)行篩選,則計(jì)算量太過龐大。簡(jiǎn)言之,以系統(tǒng)壽命為篩選依據(jù)的基礎(chǔ)器件壽命預(yù)測(cè)不便于采用上述數(shù)值型預(yù)測(cè)方法。
通過對(duì)長(zhǎng)期所積累的不同型號(hào)繼電器的壽命試驗(yàn)數(shù)據(jù)進(jìn)行分析,一個(gè)規(guī)律近來被發(fā)現(xiàn),即:繼電器個(gè)體的壽命與其生命初期的性態(tài)有關(guān),初期性態(tài)較好的個(gè)體其壽命也較長(zhǎng),反之則壽命較短。本文基于這一發(fā)現(xiàn),運(yùn)用粗糙集理論建立了繼電器產(chǎn)品初期性能參數(shù)和壽命的決策表,然后根據(jù)所建立的決策表和繼電器的初期性能參數(shù)信息即可預(yù)測(cè)繼電器個(gè)體的壽命等級(jí)(即壽命的大致范圍),最終實(shí)現(xiàn)以壽命為指標(biāo)的繼電器產(chǎn)品篩選。
繼電器的觸點(diǎn)是最易侵蝕破壞的部位,其故障率占總故障率的80%以上[7]。繼電器通電后,其觸點(diǎn)并不是馬上就達(dá)到穩(wěn)定閉合狀態(tài),而是要經(jīng)歷彈跳和動(dòng)態(tài)接觸這兩個(gè)過程。在觸點(diǎn)的彈跳過程中,彈跳時(shí)間是刻畫這一過程的重要參數(shù);而在動(dòng)態(tài)接觸過程中,動(dòng)態(tài)波動(dòng)時(shí)間和動(dòng)態(tài)接觸電阻則是刻畫該過程的重要參數(shù);觸點(diǎn)穩(wěn)定吸合后的性能則由靜態(tài)接觸電阻刻畫。動(dòng)態(tài)接觸電阻是指觸點(diǎn)閉合或斷開時(shí),隨觸點(diǎn)壓力大小而變化的接觸電阻;靜態(tài)接觸電阻是指觸點(diǎn)穩(wěn)定閉合狀態(tài)下的觸點(diǎn)間電阻。因動(dòng)態(tài)接觸電阻、靜態(tài)接觸電阻不易直接測(cè)量,故通過測(cè)量觸點(diǎn)間的動(dòng)態(tài)接觸壓降峰值(即動(dòng)態(tài)峰值壓降)、靜態(tài)接觸壓降來分別替代。
本文選取了12只某型號(hào)電磁繼電器作為試品,分三組進(jìn)行了可靠性壽命試驗(yàn)(每只試品選擇一對(duì)觸點(diǎn)),同步監(jiān)測(cè)并記錄了這些試品隨觸點(diǎn)動(dòng)作次數(shù)變化的四個(gè)性能參數(shù)(彈跳時(shí)間、動(dòng)態(tài)波動(dòng)時(shí)間、動(dòng)態(tài)峰值壓降和靜態(tài)接觸壓降)。12只試品的壽命實(shí)測(cè)值和估算值見表1,這里注釋如下:
(1)試品壽命由其觸點(diǎn)動(dòng)作次數(shù)表示,“次”表示繼電器試品正常工作時(shí)其觸點(diǎn)通、斷各一次。
(2)三組試品分別在-20℃、20℃、55℃三種環(huán)境溫度下采用定時(shí)截尾方式各自進(jìn)行試驗(yàn),定時(shí)分別為330萬次、492萬次和140萬次。至試驗(yàn)截止時(shí)僅第5、6號(hào)試品失效,壽命分別為232萬次和242萬次。失效類型前者為靜態(tài)接觸電阻大幅增加而后者為觸點(diǎn)黏連,其余試品均未失效。
(3)試驗(yàn)中,試品的線圈加額定電壓,但觸點(diǎn)并非帶額定負(fù)載,而是帶載一個(gè)微弱電信號(hào),因此表中的“壽命”并非嚴(yán)格意義上的電壽命或機(jī)械壽命。但因觸點(diǎn)實(shí)際所帶負(fù)載遠(yuǎn)小于額定值,故表中試品的壽命值應(yīng)更接近于其機(jī)械壽命。
(4)本文方法須對(duì)試品做全壽命試驗(yàn),以獲得每個(gè)試品的壽命值。但受試驗(yàn)條件所限,目前僅對(duì)該型號(hào)繼電器的12個(gè)試品分組進(jìn)行了定時(shí)截尾試驗(yàn),未獲得全部試品工作壽命的實(shí)測(cè)值。為填補(bǔ)缺失的信息,本文采用小波包變換與RBF神經(jīng)網(wǎng)絡(luò)相結(jié)合的方法[8]對(duì)各試品的各項(xiàng)性能參數(shù)進(jìn)行了預(yù)測(cè),一旦發(fā)現(xiàn)試品的任一性能參數(shù)預(yù)測(cè)值超過限定值則視其為失效,由此獲得了各試品壽命的估算值,見表1,相關(guān)的數(shù)據(jù)處理過程從略。
表1 三種溫度下繼電器試品的壽命實(shí)測(cè)值和估算值
Tab.1 The actual and estimated life of the tested relays at the temperature of -20℃, 20℃ and 55℃
由表1數(shù)據(jù)可知:即使在同一環(huán)境溫度下,各試品的壽命也會(huì)出現(xiàn)較大差別。例如試品5、試品7雖然都工作于20℃,但前者通斷232萬次即失效而后者工作至492萬次仍完好,估算其壽命可達(dá)到950萬次。經(jīng)對(duì)試驗(yàn)數(shù)據(jù)進(jìn)行全面分析發(fā)現(xiàn),試品在壽命指標(biāo)上的優(yōu)劣于其生命初期的性能參數(shù)上即有所體現(xiàn),或者說,試品在生命初期的性態(tài)好壞預(yù)示著其工作壽命的長(zhǎng)短。圖1、圖2分別為20℃下試品5和試品7在生命初期關(guān)于動(dòng)態(tài)波動(dòng)時(shí)間、靜態(tài)接觸壓降這兩項(xiàng)性能參數(shù)的對(duì)比。
由圖1、圖2兩圖易見,壽命較短的試品5在生命初期的動(dòng)態(tài)波動(dòng)時(shí)間較大,在生命最初時(shí)刻的靜態(tài)接觸壓降也較大;而壽命值相差較大的試品5和試品7在生命最初時(shí)刻之后的一段時(shí)間內(nèi),其靜態(tài)接觸壓降卻相差不大。
總結(jié)表1、圖1和圖2,可知繼電器的個(gè)體壽命與工作溫度以及生命初期的性態(tài)有關(guān),初期性態(tài)較好的個(gè)體其壽命也較長(zhǎng),反之則壽命較短。然而,繼電器的個(gè)體壽命受多個(gè)因素影響,并非所有的影響因素都是同等重要的,甚至有些是不必要的、冗余的,這些不必要的因素不但會(huì)增加預(yù)測(cè)的運(yùn)算量,還會(huì)干擾預(yù)測(cè)結(jié)果。對(duì)此,可借助于粗糙集理論刪減冗余信息。
粗糙集理論由波蘭學(xué)者Z.Pawlak創(chuàng)建于1982年,其主要優(yōu)勢(shì)是不需任何先驗(yàn)知識(shí)與信息,被作為一種數(shù)據(jù)挖掘和知識(shí)發(fā)現(xiàn)的有力工具而成功應(yīng)用在模式識(shí)別、機(jī)器學(xué)習(xí)和過程控制等領(lǐng)域。
在粗糙集理論[9]中,知識(shí)被認(rèn)為是一種對(duì)研究論域?qū)ο缶哂蟹诸惸芰Φ男畔ⅰVR(shí)表達(dá)系統(tǒng)也稱為信息系統(tǒng),多元組表示的信息系統(tǒng)如下:=<,,,>,其中為研究對(duì)象的非空有限集,稱為論域;為屬性的非空有限集,,為條件屬性,為決策屬性;是屬性值的集合,,V為屬性的值域;是一個(gè)信息函數(shù),它為每個(gè)對(duì)象的每個(gè)屬性賦予一個(gè)信息值:,,。
設(shè)是上的一個(gè)等價(jià)關(guān)系,/表示的所有等價(jià)類組成的集合。若且,則(即中所有等價(jià)類關(guān)系的交集)也是一個(gè)等價(jià)關(guān)系,稱為上的不可區(qū)分關(guān)系,記為ind()。這樣,表示與等價(jià)關(guān)系族相關(guān)的知識(shí)。
定義1 上近似、下近似與正域
稱為集合的正域。
定義2 知識(shí)的約簡(jiǎn)與核
約簡(jiǎn)是與信息系統(tǒng)的屬性全集具有相同基本集的最小屬性子集。令信息系統(tǒng)=<,,,>,設(shè),如果是獨(dú)立的,且ind()=ind(),則是的一個(gè)約簡(jiǎn)。
屬性的核是信息系統(tǒng)的屬性集的所有必要屬性構(gòu)成的集合,記為core(),常稱為的核。記red()為的所有約簡(jiǎn)集合,則有
即的所有約簡(jiǎn)的交集構(gòu)成的核。
定義3 決策規(guī)則可信度
基于粗糙集理論的預(yù)測(cè)過程可劃分為兩大步驟,即決策規(guī)則的建立與匹配,對(duì)應(yīng)于知識(shí)工程中的知識(shí)獲取與使用。以下對(duì)這兩個(gè)步驟進(jìn)行分述。
決策規(guī)則的建立(或稱獲取)方法如下:
(1)建立繼電器壽命預(yù)測(cè)的信息系統(tǒng)。將繼電器壽命試驗(yàn)中的相關(guān)數(shù)據(jù)視為信息系統(tǒng)=<,,,>,其中,表1中的試品為研究對(duì)象;屬性集=∪,其中條件屬性由試驗(yàn)條件以及用于描述繼電器觸點(diǎn)性態(tài)的四個(gè)性能參數(shù)組成,且有={環(huán)境溫度,彈跳時(shí)間,動(dòng)態(tài)波動(dòng)時(shí)間,動(dòng)態(tài)峰值壓降,靜態(tài)接觸壓降}={1,2,3,4,5},繼電器試品壽命作為決策屬性集,有={繼電器壽命};屬性值的集合通過試驗(yàn)測(cè)試得到。
信息系統(tǒng)的數(shù)據(jù)經(jīng)常是以關(guān)系表的形式給出的。關(guān)系表的每一行對(duì)應(yīng)一個(gè)具體的研究對(duì)象,每一列代表某個(gè)屬性不同的屬性值。在粗糙集理論中,稱這個(gè)關(guān)系表為決策表。繼電器壽命預(yù)測(cè)中任意一個(gè)試品的不同屬性值構(gòu)成了決策表中的一行,全部試品的同一屬性值構(gòu)成了決策表中的一列。運(yùn)用粗糙集理論處理決策表時(shí),要求決策表中的值用離散數(shù)據(jù)表示。然而繼電器產(chǎn)品壽命預(yù)測(cè)問題中的條件屬性2~5和決策屬性的值域都是實(shí)數(shù)域上的一個(gè)連續(xù)區(qū)間,故需對(duì)這些屬性的值進(jìn)行離散處理。常用的離散化方法包括經(jīng)驗(yàn)法、等距離散化、等頻離散化、基于布爾邏輯與粗糙集相結(jié)合的離散化、基于屬性重要性的離散化方法等[10]。對(duì)屬性值離散化后形成原始決策表。
(2)繼電器壽命決策表的約簡(jiǎn)。原始決策表中含有5個(gè)條件屬性,但并非所有的條件屬性都是同等重要的,甚至有些是不必要的、冗余的,不加刪減地選用所有條件屬性不但會(huì)增加計(jì)算量,還可能引入無關(guān)的隨機(jī)信息影響預(yù)測(cè)結(jié)果,因此在導(dǎo)出決策規(guī)則之前需先進(jìn)行屬性約簡(jiǎn)。常用的約簡(jiǎn)算法有基于區(qū)分矩陣的算法[11]、基于屬性依賴度的算法[12]、基于分類質(zhì)量的算法[13]和基于信息熵的算法[14]等。不同算法的特點(diǎn)不同,約簡(jiǎn)結(jié)果也不盡相同,但約簡(jiǎn)所得的核是一致的。經(jīng)過屬性約簡(jiǎn),可得到對(duì)繼電器壽命最有表征力的最簡(jiǎn)屬性集,該集合保留了與原決策信息系統(tǒng)具有相同決策能力的條件屬性。實(shí)際預(yù)測(cè)中僅考慮約簡(jiǎn)后的條件屬性即可。
(3)繼電器壽命決策規(guī)則的建立。最簡(jiǎn)決策表只是原有知識(shí)表達(dá)系統(tǒng)的精簡(jiǎn)表達(dá),并沒有形成有效的決策規(guī)則。產(chǎn)生決策規(guī)則通常有兩種方式:一種方法是將約簡(jiǎn)后的決策表直接用產(chǎn)生式規(guī)則表示出來,并給出每條規(guī)則的可信度,這樣得到的決策規(guī)則即是標(biāo)準(zhǔn)決策規(guī)則;另一種為了提髙決策規(guī)則與新對(duì)象匹配的可能性,在生成標(biāo)準(zhǔn)決策規(guī)則后,采用近似規(guī)則生成算法增加規(guī)則數(shù)量[15]。近似決策規(guī)則是提高決策規(guī)則對(duì)新對(duì)象的匹配能力的一種方法,然而在分類系統(tǒng)的實(shí)際應(yīng)用中,仍然會(huì)經(jīng)常發(fā)生新對(duì)象找不到匹配規(guī)則的現(xiàn)象,此時(shí)需要恰當(dāng)?shù)囊?guī)則匹配算法實(shí)現(xiàn)決策規(guī)則與被預(yù)測(cè)對(duì)象的匹配。
對(duì)繼電器產(chǎn)品的壽命預(yù)測(cè),實(shí)際上是根據(jù)被預(yù)測(cè)對(duì)象的條件屬性進(jìn)行規(guī)則匹配,把獲得匹配的規(guī)則結(jié)論中關(guān)于繼電器壽命的值賦給待預(yù)測(cè)產(chǎn)品,以此作為該產(chǎn)品的壽命預(yù)測(cè)值,這里,繼電器的壽命預(yù)測(cè)值是一個(gè)區(qū)間值而非一個(gè)確定的實(shí)數(shù)。
被預(yù)測(cè)對(duì)象和決策規(guī)則的匹配分為以下幾種情形(類似問題在其他一些文獻(xiàn)[16]中也被討論過):①對(duì)象與一條規(guī)則相匹配;②對(duì)象與多條規(guī)則相匹配,且結(jié)論相同;③對(duì)象與多條規(guī)則相匹配,但結(jié)論不同;④對(duì)象與任一規(guī)則均不匹配。
情形①和②的預(yù)測(cè)結(jié)果是確定的;情形③規(guī)則之間有沖突,即條件一致,結(jié)論不同,因此預(yù)測(cè)結(jié)果有一定的可信度,但并不完全確定??捎猛镀狈?、對(duì)象跟蹤法等方法[17]解決;在情形④中,不存在與被預(yù)測(cè)對(duì)象相匹配的規(guī)則。目前解決此類問題的主要方法有部分匹配法和相近規(guī)則法。部分匹配法是將被預(yù)測(cè)對(duì)象的部分屬性與決策規(guī)則進(jìn)行部分匹配,而相近規(guī)則法則是通過某種算法從決策規(guī)則中找出與預(yù)測(cè)對(duì)象屬性最接近的決策規(guī)則。此外還有其他一些規(guī)則匹配法,例如文獻(xiàn)[18]使用了一種距離度量法來匹配規(guī)則。具體方法為:
設(shè)被預(yù)測(cè)對(duì)象為,條件屬性值為{1(),2(),…,y()},規(guī)則集中規(guī)則的條件屬性值為{1(),2(),…,y()}。統(tǒng)計(jì)規(guī)則集中所有規(guī)則的條件屬性值與的條件屬性值不匹配的數(shù)目,將該數(shù)目最小的規(guī)則取出,構(gòu)成候選規(guī)則集。若候選規(guī)則集中只有1條規(guī)則,以該規(guī)則為匹配規(guī)則,否則按式(5)度量各候選規(guī)則與的距離,即相似程度,以距離最近的規(guī)則為匹配規(guī)則。
式中,ymax和ymin分別是y最大值和最小值;為決策規(guī)則中條件屬性的個(gè)數(shù)。
確定規(guī)則后,將的決策屬性值劃歸該規(guī)則的決策類別,即由規(guī)則的結(jié)論確定的壽命預(yù)測(cè)結(jié)果。
3.3.1 繼電器壽命決策表的構(gòu)建與約簡(jiǎn)
從第1節(jié)所述12組繼電器壽命數(shù)據(jù)中,隨機(jī)選用其中的9組(用1~9表示)建立壽命決策規(guī)則,剩下的3組(用10~12表示)用來驗(yàn)證規(guī)則的有效性。提取1~9的試驗(yàn)最終壽命數(shù)據(jù)及其生命中最初10萬次內(nèi)的性能參數(shù)數(shù)據(jù),并對(duì)生命中最初10萬次內(nèi)的性能參數(shù)數(shù)據(jù)做平均值統(tǒng)計(jì),得到9個(gè)試品的性態(tài)參數(shù)與壽命信息表,見表2。
表2 繼電器性能參數(shù)與壽命信息表
Tab.2 Information table of parameter performance and life of relays
其中,決策屬性(繼電器壽命)的離散可根據(jù)產(chǎn)品的實(shí)際篩選需求進(jìn)行,本實(shí)例中將其劃分為1、2、3三個(gè)等級(jí)(也可以認(rèn)為是三個(gè)模糊語言值,例如“小、中、大”),條件屬性則根據(jù)數(shù)據(jù)特點(diǎn)采用經(jīng)驗(yàn)法進(jìn)行離散。具體離散方法如下:
對(duì)于環(huán)境溫度,“1”表示其值屬于[L, 0℃這一區(qū)間;類似地,“2、3”分別表示其值屬于[0,50℃和[50,H]℃區(qū)間。其中,L、H分別為繼電器可正常工作的環(huán)境溫度下限、上限,這兩個(gè)值一般由產(chǎn)品生產(chǎn)廠家提供,也可以通過實(shí)驗(yàn)獲得。
對(duì)于彈跳時(shí)間,“1、2、3”分別表示其值屬于[0,0.5ms、[0.5,0.8ms、[0.8,+∞ms區(qū)間。
對(duì)于動(dòng)態(tài)波動(dòng)時(shí)間,“1、2、3”分別表示其值屬于[0,0.6ms、[0.6,1.0ms、[1.0,+∞ms區(qū)間。
對(duì)于動(dòng)態(tài)峰值壓降,“1、2、3”分別表示其值屬于[0,600mV、[600,1 000mV、[1 000,+∞mV區(qū)間。
對(duì)于靜態(tài)接觸壓降,“1、2、3”分別表示其值屬于[0,400mV、[400,800mV、[800,+∞]mV區(qū)間。
對(duì)于繼電器的壽命,“1、2、3”分別表示其值屬于[0,600萬次、[600,900萬次、[900,+∞]萬次 區(qū)間。
數(shù)據(jù)處理過程如下:
(1)繼電器壽命決策表的建立。按照上述離散規(guī)則,對(duì)表2進(jìn)行離散操作,并按條件屬性和決策屬性分類,其中1,2,…,5均為條件屬性,分別表示環(huán)境溫度、彈跳時(shí)間、動(dòng)態(tài)波動(dòng)時(shí)間、動(dòng)態(tài)峰值壓降、靜態(tài)接觸壓降;為決策屬性,表示工作壽命。如此便得到繼電器的壽命決策表,見表3。
表3 繼電器壽命決策表
Tab.3 Life decision table of relays
(2)繼電器壽命決策表的約簡(jiǎn)。求取與繼電器壽命有關(guān)的最簡(jiǎn)屬性集,可基于粗糙集理論中的屬性約簡(jiǎn)來完成。以下是對(duì)表3所示的繼電器壽命決策表的屬性約簡(jiǎn)過程,即條件屬性“”對(duì)決策屬性“繼電器壽命”的約簡(jiǎn)。具體約簡(jiǎn)過程如下:
記和之間依賴度為,=說明屬性集={溫度, 彈跳時(shí)間,動(dòng)態(tài)波動(dòng)時(shí)間,動(dòng)態(tài)峰值壓降,靜態(tài)接觸壓降}對(duì)繼電器壽命的分類是充分集合,完全依賴于。
由以上計(jì)算可知,1~4分別相對(duì)于不可單獨(dú)約簡(jiǎn),5相對(duì)于可單獨(dú)約簡(jiǎn)。因此,得到的核為={環(huán)境溫度、彈跳時(shí)間、動(dòng)態(tài)波動(dòng)時(shí)間、動(dòng)態(tài)峰值壓降},這也是約簡(jiǎn)所得的最簡(jiǎn)屬性集。
3.3.2 繼電器壽命決策規(guī)則的獲取
通過前述的屬性約簡(jiǎn)及求核,可以得到簡(jiǎn)化的壽命決策表見表4。
表4 繼電器壽命決策規(guī)則表
Tab.4 Life decision rules table of relays
根據(jù)表4可得到如下一組決策規(guī)則:
Rule1:(1,1)∧(2,2)∧(3,1)∧(3,3)→(,2)
Rule2:(1,1)∧(2,1)∧(3,2)∧(3,3)→(,2)
Rule3:(1,2)∧(2,2)∧(3,2)∧(3,1)→(,1)
Rule4:(1,2)∧(2,1)∧(3,1)∧(3,1)→(,1)
Rule5:(1,2)∧(2,1)∧(3,1)∧(3,2)→(,3)
Rule6:(1,2)∧(2,2)∧(3,1)∧(3,1)→(,3)
Rule7:(1,3)∧(2,2)∧(3,1)∧(3,3)→(,1)
Rule8:(1,3)∧(2,1)∧(3,2)∧(3,3)→(,1)
Rule9:(1,1)∧(2,3)∧(3,3)∧(3,3)→(,2)
根據(jù)式(4)求得本實(shí)例中每條決策規(guī)則的可信度均為1,說明所得決策規(guī)則都是確定性規(guī)則。實(shí)際應(yīng)用中可將1作為決策規(guī)則的初始可信度,在出現(xiàn)與決策規(guī)則相沖突(條件屬性值相同,決策屬性值不同)的實(shí)例后,可根據(jù)沖突的頻率、基于概率論方法調(diào)整規(guī)則可信度(通常為小于1的正實(shí)數(shù))。當(dāng)規(guī)則可信度下降到一定閾值(比如0.5),此條規(guī)則就不適宜再直接使用,此時(shí)需要重新選擇樣本進(jìn)行規(guī)則提取,從而實(shí)現(xiàn)對(duì)規(guī)則組的修正、補(bǔ)充。
3.3.3 壽命預(yù)測(cè)及決策規(guī)則有效性分析
根據(jù)壽命決策規(guī)則對(duì)繼電器1012進(jìn)行壽命預(yù)測(cè),首先提取1012的生命初期性能參數(shù)數(shù)據(jù),并做平均值統(tǒng)計(jì),再利用前述的離散化方法進(jìn)行離散分類,然后匹配決策規(guī)則進(jìn)行壽命預(yù)測(cè)。被預(yù)測(cè)繼電器的生命初期性能參數(shù)信息見表5。
表5 待預(yù)測(cè)繼電器的工作壽命與性能參數(shù)
Tab.5 The working life and performance parameters of the predicted relays
繼電器編號(hào)數(shù)據(jù)溫度/℃彈跳時(shí)間/ms動(dòng)態(tài)波動(dòng)時(shí)間/ms動(dòng)態(tài)峰值壓降/mV工作壽命/萬次 x10實(shí)際值-200.530.481 035750 離散值12132 x11實(shí)際值550.590.491 218540 離散值32131 x12實(shí)際值550.401.001 230370 離散值31331
對(duì)照3.3.2建立的壽命決策規(guī)則,繼電器10匹配第一條規(guī)則,并且可信度為1,對(duì)應(yīng)預(yù)測(cè)壽命等級(jí)為2,與實(shí)際壽命結(jié)果一致(10實(shí)際壽命為750萬次);繼電器11匹配第七條規(guī)則,并且可信度為1,對(duì)應(yīng)預(yù)測(cè)壽命等級(jí)為1,與實(shí)際壽命結(jié)果一致(11實(shí)際壽命為540萬次);然而,決策規(guī)則中沒有與繼電器12相匹配的規(guī)則,這里參照文獻(xiàn)[18]中的距離度量法進(jìn)行規(guī)則匹配。根據(jù)上述方法,發(fā)現(xiàn)繼電器12的候選規(guī)則只有一條規(guī)則,即Rule8,故以該規(guī)則為匹配規(guī)則。Rule8對(duì)應(yīng)的壽命等級(jí)為1,所以可推測(cè)繼電器12的壽命等級(jí)為1。這一推測(cè)與繼電器實(shí)際壽命370萬次一致。
通過分析繼電器的壽命試驗(yàn)數(shù)據(jù),發(fā)現(xiàn)繼電器個(gè)體的壽命與其生命初期的性態(tài)有著某種關(guān)系。本文基于這一發(fā)現(xiàn),提出了一種基于粗糙集理論的繼電器壽命預(yù)測(cè)方法,預(yù)測(cè)的精細(xì)程度取決于壽命的離散分類等級(jí)數(shù)。離散分類等級(jí)數(shù)越高,每一等級(jí)所表征的區(qū)間越窄,預(yù)測(cè)結(jié)果也就越精細(xì)。當(dāng)僅需判斷繼電器個(gè)體是否能達(dá)到期望的壽命值時(shí),壽命的離散分類等級(jí)數(shù)為2即可,分別代表“達(dá)到”和“未達(dá)到”,這種簡(jiǎn)單易行的方法為以壽命為指標(biāo)的產(chǎn)品篩選帶來極大便利。
與當(dāng)前常用的繼電器產(chǎn)品壽命預(yù)測(cè)和產(chǎn)品篩選方法相比,本文方法的特點(diǎn)體現(xiàn)在:①簡(jiǎn)單、快捷、運(yùn)算量小、所依賴的信息量少,僅需測(cè)出繼電器個(gè)體在生命初期的性能參數(shù)數(shù)據(jù)。②壽命預(yù)測(cè)的對(duì)象為繼電器產(chǎn)品個(gè)體而非批產(chǎn)品,因此非常適用于以工作壽命為指標(biāo)的繼電器產(chǎn)品等級(jí)劃分與篩選。③由于僅需對(duì)繼電器在其生命初期的性能參數(shù)進(jìn)行測(cè)試,故在篩選過程中不會(huì)過多壓縮其有效服役期,同時(shí)也大大降低了產(chǎn)品篩選的成本。
該方法亦可推廣應(yīng)用于其他產(chǎn)品的壽命預(yù)測(cè)與產(chǎn)品篩選,只要該產(chǎn)品的壽命與其生命初期特征具有相關(guān)性。
[1] 張又衡. 車輛繼電器壽命綜合實(shí)驗(yàn)臺(tái)[D]. 長(zhǎng)春: 吉林大學(xué), 2007.
[2] 張菲菲, 李志剛. 基于BP神經(jīng)網(wǎng)絡(luò)的繼電器剩余壽命預(yù)測(cè)[J]. 低壓電器, 2012(1): 11-14.Zhang Feifei, Li Zhigang. Remaining lifetime prediction of relay based on BP neural network[J].Low Voltage Apparatus, 2012(1): 11-14.
[3] 李華, 孫東旺, 賀鵬舉, 等. 基于超程時(shí)間回歸模型的繼電器壽命預(yù)測(cè)方法[J]. 電工技術(shù)學(xué)報(bào), 2013, 28(S2): 414-417. Li Hua, Sun Dongwang, He Pengju, et al. The method of relay life prediction based on the regression model of super-path time[J]. Transactions of China Electrotech- nical Society, 2013, 28(S2): 414-417.
[4] Ye Xuerong, Yu Qiong, Zhai Guofu. Reliability assessment for electromagnetic relay based on time parameters degradation[C]//2010 11th International Conference on Electronic Packaging Technology & High Density Packaging, Xi’an, 2010: 1269-1272.
[5] Yu Qiong, Qi Ming, Wang Shujuan, et al. Research on life prediction based on wavelet transform and ARMA model for space relay[C]//4th IEEE Conference on Industrial Electronics and Applications, Xi’an, 2009: 1275-1280.
[6] 任立. 航天繼電器可靠性壽命試驗(yàn)分析系統(tǒng)的研究[D]. 哈爾濱: 哈爾濱工業(yè)大學(xué), 2009.
[7] 王淑娟, 余瓊, 翟國(guó)富. 電磁繼電器接觸失效機(jī)理判別方法[J]. 電工技術(shù)學(xué)報(bào), 2010, 25(8): 38-44.Wang Shujuan, Yu Qiong, Zhai Guofu. Discrimination method of contact failure mechanisms for electro- magnetic apparatus[J]. Transactions of China Electro- technical Society, 2010, 25(8): 38-44.
[8] 馮興杰, 潘文欣, 盧楠. 基于小波包的RBF神經(jīng)網(wǎng)絡(luò)網(wǎng)絡(luò)流量混沌預(yù)測(cè)[J]. 計(jì)算機(jī)工程與設(shè)計(jì), 2012, 33(5): 1681-1686.Feng Xingjie, Pan Wenxin, Lu nan. Research on chaotic prediction method of network traffic based on wavelet packet[J]. Computer Engineering and Design, 2012, 33(5): 1681-1686.
[9] Liu Guilong. Rough set theory based on two universal sets and its applications[J]. Knowledge-Based Systems, 2010, 23(2): 110-115.
[10] 高爽, 冬雷, 高陽, 等. 基于粗糙集理論的中長(zhǎng)期風(fēng)速預(yù)測(cè)[J].中國(guó)電機(jī)工程學(xué)報(bào), 2012, 32(1):32-37.Gao Shuang, Dong Lei, Gao Yang, et al. Mid-long term wind speed prediction based on rough set theory[J]. Proceedings of the CSEE, 2012, 32(1): 32-37.
[11] Yanyan Yang, Degang Chen, Ze Dong. Novel algo- rithms of attribute reduction with variable precision rough set model[J]. Neurocomputing, 2014, 139(2): 336-344.
[12] Kai Zheng, Jie Hu, Zhenfei Zhan. An enhancement for heuristic attribute reduction algorithm in rough set[J]. Expert Systems with Applications, 2014, 41(15): 6748-6754.
[13] 高曉紅. 基于分類質(zhì)量的決策系統(tǒng)屬性約簡(jiǎn)新算法[J]. 計(jì)算機(jī)與現(xiàn)代化, 2010(1): 19-22. Gao Xiaohong. Novel attribute reduction algorithm of decision system based on classification quality[J]. Computer and Modernization, 2010(1): 19-22.
[14] Jing Tian, Quan Wang, Yu Bing. A rough set algo- rithm for attribute reduction via mutual information and conditional entropy[C]//Fuzzy Systems and Know- ledge Discovery(FSKD), Shenyang, 2013: 567-571.
[15] Wang Yun, Wu Zhiqiang. Generating algorithm of approximate decision rules and its application[C]// Granular Computing(GrC), San Jose, CA, 2010: 830- 833.
[16] Qiong Wu, Xin Pan. A novel remote sensing classify- cation rule extraction method based on discrete rough set[C]//Fuzzy Systems and Knowledge Discovery (FSKD), Shanghai, 2011: 330-334.
[17] 楊傳振, 朱玉全, 陳耿. 一種基于粗糙集屬性約簡(jiǎn)的多分類器集成方法[J]. 計(jì)算機(jī)應(yīng)用研究, 2012, 29(5): 1648-1650. Yang Chuanzhen, Zhu Yuquan, Chen Geng. New classifier ensemble method based on rough set attributereduction[J]. Application Research of Computers, 2012, 29(5): 1648-1650.
[18] 龍丹麗, 黎靜華, 韋化. 粗糙集法解多環(huán)境因素影響的母線負(fù)荷預(yù)測(cè)問題[J]. 電網(wǎng)技術(shù), 2013, 37(5): 1335-1340. Long Danli, Li Jinghua, Wei Hua. Solution of multi environmental factor-influenced bus load forecasting by rough set method[J]. Power System Technology, 2012, 37(5): 1335-1340.
The Life Prediction Method of Relay Based on Rough Set Theory and Relay’s Initial Life Information
(Province-Ministry Joint Key Laboratory of Electromagnetic Field and Electrical Apparatus Reliability Hebei University of Technology Tianjin 300130 China)
Working life is an important feature of product reliability. However the productscreening method using life as the sole index hasn’t been thoroughly studied for long life products such as relays. After the life test data of relays was analyzed,a fact was foundrecently that there are some relations between relay’s initial working state and its working life, i.e., better performance of the initial state, longer life expectancy of the individual relay, and vice versa. Based on the fact, a method of the individual relay’s life predicting and screening was presented. Firstly, obtained the initial state performance parameters values and life values of the sample relays through life test. Secondly, extracted a set ofparameters associated with life based on the method of attribute reduction from the rough set theory, then established a set of life decision rules based on the extracted parameters. After that, for any one of the relay whose model same as the samples, its life could be predicted only according to its initial state information and the rules, then the screening using life as the sole index could be achieved. At last, the availability of the above method was proved by a practical example of relay.
Relay, initial state information, life prediction, product screening, rough set theory
TM581.3
李玲玲 女,1968年生,博士,教授,博士生導(dǎo)師,研究方向?yàn)殡娖骺煽啃?、電力系統(tǒng)及其自動(dòng)化。
E-mail: lilingling@hebut.edu.cn
李志剛 男,1958年生,博士,教授,博士生導(dǎo)師,研究方向?yàn)殡娮与娖?、電器可靠性及其檢測(cè)技術(shù)。
E-mail: zgli@hebut.edu.cn(通信作者)
2014-07-28 改稿日期 2014-09-19
國(guó)家自然科學(xué)基金(61072100、51377044),河北省自然科學(xué)基金(E2014202230)和河北省高校創(chuàng)新團(tuán)隊(duì)領(lǐng)軍人才培育計(jì)劃(LJRC003)資助項(xiàng)目。