基于混合威布爾分布的加工中心可靠性評(píng)估
引言威布爾分布是可靠性評(píng)估常用的模型，一般用于建模具有 單一故障模式的故障數(shù)據(jù)，但加工中心是復(fù)雜的機(jī)械系統(tǒng)，其 包含著眾多的零部件，每個(gè)系統(tǒng)故障可能是在多種故障機(jī)制共 同作用下發(fā)生的，其故障數(shù)據(jù)在威布爾概率紙上表現(xiàn)為曲線， 如果用原始的威布爾分布曲線來描述，會(huì)出現(xiàn)較大的誤差[1]。 利用各種改進(jìn)的威布爾模型可以較好地解決這類問題，如混合 威布爾分布[2-3]?；旌蟇eibull分布應(yīng)用的一個(gè)難點(diǎn)是參數(shù)估計(jì)問題。常用來 進(jìn)行參數(shù)估計(jì)的方法有作圖法、回歸分析法、矩估計(jì)法、最小 二乘法、極大似然估計(jì)法[4-6] (Maximum Likelihood Estimation,MLE)等。在實(shí)際應(yīng)用中，極大似然估計(jì)是最常用的 方法，也是最重要的方法。但應(yīng)用極大似然估計(jì)等傳統(tǒng)方法來進(jìn)行混合Weibull模型參數(shù)估計(jì)，一般要求解聯(lián)立的超越方程 組，相當(dāng)復(fù)雜。另外在小子樣故障數(shù)據(jù)的情況下，極大似然估 計(jì)的結(jié)果可能會(huì)產(chǎn)生比較大的偏差[7-8]。文獻(xiàn)[9]采用混合威布爾 分布描述了 FIAT汽車零部件可靠性及其故障發(fā)生規(guī)律，并且用 極大似然方法估計(jì)其分布參數(shù)，但可以看到，其故障數(shù)據(jù)為大 樣本故障數(shù)據(jù)。文獻(xiàn)[10]采用非線性的最小二乘法求解混合威布 爾模型的各參數(shù)，其結(jié)果依賴于參數(shù)初值，對(duì)于兩個(gè)子分布數(shù) 據(jù)融合度比較大時(shí)，其結(jié)果精度不高?；谪惾~斯（Bayes)理論的可靠性評(píng)估方法綜合了驗(yàn)前信 息和樣本信息，是解決小子樣故障數(shù)據(jù)較好的一種方法[11]。 Bayes理論的基本原理是利用驗(yàn)前分布和樣本信息來計(jì)算后驗(yàn) 分布，從而估算變量的點(diǎn)估計(jì)和置信區(qū)間，并進(jìn)一步推導(dǎo)其他 相關(guān)可靠性特征量的估計(jì)值。很多文獻(xiàn)將Bayes應(yīng)用到單威布 爾分布的參數(shù)求解中，并取得了很好的效果，但由于混合威布爾分布存在兩個(gè)子分布，用同一組歷史故障數(shù)據(jù)不可能求出兩 個(gè)先驗(yàn)分布的參數(shù)值，這是將Bayes應(yīng)用到混合威布爾分布中 的第一個(gè)難點(diǎn)；混合威布爾分布無共軛先驗(yàn)分布，將Bayes理 論應(yīng)用到混合威布爾分布參數(shù)求解時(shí)會(huì)產(chǎn)生很大的計(jì)算量，這 是將Bayes應(yīng)用到混合威布爾分布中的第二個(gè)難點(diǎn)。文獻(xiàn)[12] 將Bayes方法應(yīng)用于具有共軛先驗(yàn)分布的混合指數(shù)分布的參數(shù) 求解中，使得其計(jì)算量大大降低，而在復(fù)雜機(jī)械系統(tǒng)中，混合 威布爾分布是比混合指數(shù)分布更常用的一種模型；另外，文獻(xiàn)[12] 在確定先驗(yàn)分布各參數(shù)時(shí)，通過專家經(jīng)驗(yàn)法直接確定各先 驗(yàn)分布的參數(shù)值，先驗(yàn)信息的的可靠性與否很大程度上決定了 參數(shù)后驗(yàn)估計(jì)值的精確程度[13]，因此，經(jīng)驗(yàn)法確定先驗(yàn)分布失 去了 Bayes估計(jì)充分利用、挖掘先驗(yàn)信息的優(yōu)勢(shì)，引入了人為 不確定因素。針對(duì)上述問題，本文給出了一種將Bayes理論應(yīng)用于混合 威布爾分布的方法。該方法首先通過故障特征屬性的概念，用 模糊聚類的方法將可靠性試驗(yàn)樣本故障數(shù)據(jù)分為兩個(gè)子樣本， 以作為混合威布爾分布中的兩個(gè)子故障數(shù)據(jù)，然后通過各子故 障的大量歷史故障數(shù)據(jù)，確定混合威布爾分布的兩個(gè)形狀參數(shù)， 然后將混合威布爾分布轉(zhuǎn)化為具有共軛先驗(yàn)分布的混合指數(shù)分 布。再通過各子故障數(shù)據(jù)的大量歷史故障數(shù)據(jù)求得混合指數(shù)分 布的Bayes先驗(yàn)分布，利用Bayes方法得到混合指數(shù)分布的后 驗(yàn)分布并求得其混合指數(shù)分布的參數(shù)值，最后將其還原為混合 威布爾分布的各參數(shù)值，求得加工中心的可靠性評(píng)估結(jié)果。2基于故障特征屬性的故障數(shù)據(jù)模糊聚類分析本文以13臺(tái)某型號(hào)臥式加工中心定時(shí)截尾試驗(yàn)數(shù)據(jù)進(jìn)行 評(píng)估，截尾時(shí)間^=1000辦，截尾時(shí)間內(nèi)共有9臺(tái)加工中心出 現(xiàn)故障，故障集F={刀庫亂刀，刀具外部冷卻故障，主軸換刀 故障，主軸振動(dòng)異響，Y軸換刀不到位，Y軸振動(dòng)異響，B軸 回轉(zhuǎn)精度降低，托架交換速度過快，托架交換停止}，故障時(shí) 間為 r={165,207,254,307,366,436,511，591，778}。為將故障數(shù)據(jù)劃分為兩類具有相似故障機(jī)制的子故障數(shù) 據(jù)，本文提出故障特征屬性的概念，通過故障樹分析（Failure TreeAnalysis，F(xiàn)TA)建立故障的故障特征屬性集合，并根據(jù)各 故障特征屬性與故障應(yīng)力的關(guān)系，得出各故障對(duì)故障應(yīng)力的評(píng) 價(jià)值，以此作為故障信息序列，進(jìn)行模糊聚類分析，將故障數(shù) 據(jù)進(jìn)行分類，得到具有相似故障機(jī)制的兩個(gè)子故障數(shù)據(jù)。2.1故障特征屬性和故障信息序列故障通常采用故障應(yīng)力、故障機(jī)制和故障模式來表征。故 障應(yīng)力是引起故障的物理?xiàng)l件，一般來說精工機(jī)床所有的零件 都是處在多種應(yīng)力存在的復(fù)雜的物理環(huán)境中的，因此故障的發(fā) 生也是多應(yīng)力共同作用的結(jié)果；故障機(jī)制是故障應(yīng)力發(fā)生作用， 直至引起故障的動(dòng)態(tài)或靜態(tài)過程；故障模式是作為故障機(jī)制的 結(jié)果而產(chǎn)生的故障狀態(tài)及現(xiàn)象，也是加工中心使用現(xiàn)場(chǎng)所記錄 的主要故障信息。由于機(jī)械系統(tǒng)的復(fù)雜性，這三要素間往往存 在多種組合關(guān)系，即如圖1所示。對(duì)故障機(jī)制影響******的因素 是故障應(yīng)力，相同的故障應(yīng)力作用過程相似。本文采用故障應(yīng)力的相似性表征故障機(jī)制的相似性，用各 個(gè)故障對(duì)加工中心工作時(shí)存在的所有故障應(yīng)力的模糊評(píng)判結(jié)果 作為故障信息序列，進(jìn)而通過模糊聚類的方法，實(shí)現(xiàn)故障數(shù)據(jù) 的聚類。建立故障信息序列的過程也就是分析不同的故障與所 有故障應(yīng)力之間關(guān)系的過程，為表征故障的特性，以便于深入 分析故障與故障應(yīng)力的關(guān)系，本文定義故障特征屬性的概念如 下：故障特征屬性：導(dǎo)致故障發(fā)生的隨機(jī)故障事件或最小隨機(jī) 故障事件的集合。隨機(jī)故障事件是在故障應(yīng)力作用下發(fā)生的，具有一定隨機(jī) 性的事件，其發(fā)生不依賴于其它故障事件，但會(huì)導(dǎo)致其它故障 事件的發(fā)生，因此隨機(jī)故障事件也就是故障樹分析（Failure Tree Analysis，F(xiàn)TA)中的底事件，而故障特征屬性就是故障樹的最 小割集。每個(gè)最小割集中表示一個(gè)隨機(jī)故障事件，隨機(jī)故障事 件的發(fā)生導(dǎo)致系統(tǒng)故障的發(fā)生，因此隨機(jī)故障事件是故障應(yīng)力 作用的直接對(duì)象，是故障機(jī)制的發(fā)生載體，故以隨機(jī)故障事件 表征系統(tǒng)故障能更容易得到系統(tǒng)故障與故障應(yīng)力之間的關(guān)系。4結(jié)論本文針對(duì)復(fù)雜機(jī)械系統(tǒng)具有多故障模式以及小子樣可靠性 試驗(yàn)故障數(shù)據(jù)的特點(diǎn)，用混合威布爾分布模型和Bayes理論進(jìn) 行可靠性評(píng)估。為解決混合威布爾模型無共軛先驗(yàn)分布問題，提出故障特征屬性的概念，用模糊聚類的方法實(shí)現(xiàn)故障數(shù)據(jù)的 分類，得到混合威布爾模型的形狀參數(shù)，將混合威布爾分布轉(zhuǎn) 化為具有共軛先驗(yàn)分布的混合指數(shù)分布。利用故障分類后相應(yīng) 的歷史故障數(shù)據(jù)求解混合模型的先驗(yàn)分布也是本文的創(chuàng)新點(diǎn)， 這提高了求解混合模型先驗(yàn)分布的準(zhǔn)確性。 ?每個(gè)最小割集中表示一個(gè)隨機(jī)故障事件，隨機(jī)故障事 件的發(fā)生導(dǎo)致系統(tǒng)故障的發(fā)生，因此隨機(jī)故障事件是故障應(yīng)力 作用的直接對(duì)象，是故障機(jī)制的發(fā)生載體，故以隨機(jī)故障事件 表征系統(tǒng)故障能更容易得到系統(tǒng)故障與故障應(yīng)力之間的關(guān)系。本文由海天精工整理發(fā)表文章均來自網(wǎng)絡(luò)僅供學(xué)習(xí)參考，轉(zhuǎn)載請(qǐng)注明！