陳鑫,陳富民

(西安交通大學(xué)機(jī)械制造系統(tǒng)工程國(guó)家重點(diǎn)實(shí)驗(yàn)室,710049,西安)

隨著經(jīng)濟(jì)的發(fā)展,用戶需求日益多樣化,多品種小批量生產(chǎn)已逐漸成為了制造業(yè)的主導(dǎo)模式[1]。然而,由于該種模式數(shù)據(jù)量少、產(chǎn)品參數(shù)多變等特點(diǎn),傳統(tǒng)的統(tǒng)計(jì)過(guò)程控制(SPC)方法的應(yīng)用收效甚微。因此,如何有效解決多品種小批量生產(chǎn)的質(zhì)量控制問(wèn)題成為當(dāng)前研究的熱點(diǎn)[2]。

近年來(lái),國(guó)內(nèi)外學(xué)者在此方面進(jìn)行了大量研究。部分學(xué)者將擴(kuò)充數(shù)據(jù)量作為研究重點(diǎn),如余忠華等提出了基于相似元的工序相似性分析方法,并應(yīng)用于SPC技術(shù)中[3];吳小勇引入了直覺(jué)模糊集理論來(lái)解決成組工序構(gòu)建中的相似性評(píng)判問(wèn)題,將相似工序分類成組進(jìn)行分析[4];WIEDERHOLD等通過(guò)識(shí)別相似生產(chǎn)過(guò)程并將它們的數(shù)據(jù)進(jìn)行聚類從而增加樣本量[5]。還有一些學(xué)者致力于質(zhì)量控制方法研究,如牛占文等針對(duì)多品種小批量生產(chǎn)模式下SPC應(yīng)用效果不佳的問(wèn)題,提出結(jié)合應(yīng)用通用控制圖、累計(jì)和控制圖、指數(shù)加權(quán)移動(dòng)平均值控制圖的解決方法[6];GU等分析了各種質(zhì)量控制方法的不足之處,提出了一種新型的t控制圖來(lái)監(jiān)測(cè)產(chǎn)品質(zhì)量[7]。

現(xiàn)有大部分研究只關(guān)注于歷史生產(chǎn)數(shù)據(jù)或當(dāng)前生產(chǎn)數(shù)據(jù)之一,產(chǎn)生了很大的信息浪費(fèi)。貝葉斯理論從歷史數(shù)據(jù)等先驗(yàn)知識(shí)中提取所需信息,結(jié)合現(xiàn)有數(shù)據(jù)就能有效解決多品種小批量生產(chǎn)模式下數(shù)據(jù)量少的問(wèn)題。目前,已有一些學(xué)者將貝葉斯理論應(yīng)用于中小批量生產(chǎn)的質(zhì)量控制中,取得了一定的成果[8-10]。然而,此類研究大多只利用單一產(chǎn)品的質(zhì)量數(shù)據(jù),不太適用于多品種小批量生產(chǎn)模式。

針對(duì)新模式下的工序質(zhì)量控制問(wèn)題,本文基于貝葉斯理論,提出了一種動(dòng)態(tài)質(zhì)量控制方法。首先根據(jù)成組技術(shù)和統(tǒng)計(jì)變換從歷史生產(chǎn)批次中擴(kuò)充數(shù)據(jù)量作為先驗(yàn)信息;然后融合指數(shù)加權(quán)思想,應(yīng)用共軛貝葉斯方法建模,計(jì)算相應(yīng)的工序能力指數(shù),隨著生產(chǎn)過(guò)程的進(jìn)行實(shí)時(shí)計(jì)算控制限,實(shí)現(xiàn)動(dòng)態(tài)質(zhì)量控制;最后通過(guò)仿真比較本文方法與現(xiàn)有常用方法的優(yōu)劣,結(jié)果表明,本文方法在生產(chǎn)的各個(gè)階段都有更好的控制效果和異常檢出力,可為多品種小批量生產(chǎn)模式提供一種有效的質(zhì)量控制方法。

1 選取先驗(yàn)信息

在開(kāi)始質(zhì)量控制之前,應(yīng)該分析企業(yè)生產(chǎn)狀況,以確定需要進(jìn)行質(zhì)量控制的關(guān)鍵質(zhì)量特性等信息。針對(duì)關(guān)鍵質(zhì)量特性,根據(jù)成組技術(shù)收集過(guò)去一定時(shí)間段內(nèi)滿足要求的歷史質(zhì)量特征數(shù)據(jù),通過(guò)統(tǒng)計(jì)變換方法,使各批次歷史數(shù)據(jù)滿足“獨(dú)立同分布”的要求,然后以變換后數(shù)據(jù)作為貝葉斯建模的先驗(yàn)信息。

1.1 成組技術(shù)

按照一定的相似性準(zhǔn)則將產(chǎn)品分類成組是成組技術(shù)的核心[11]。將成組技術(shù)應(yīng)用到多品種小批量質(zhì)量控制,企業(yè)應(yīng)該側(cè)重于對(duì)產(chǎn)品質(zhì)量影響較大的關(guān)鍵工序來(lái)制定相似性分類標(biāo)準(zhǔn)。

目前,相似工序的分類方法主要有4種:目測(cè)分類、生產(chǎn)流程分析[12]、編碼分類[13]、聚類分析[14]。企業(yè)應(yīng)根據(jù)實(shí)際情況選擇適合的方法,然后從人、設(shè)備、物料、方法、環(huán)境、測(cè)量6個(gè)方面[15]制定分類準(zhǔn)則,將相似工序分類成組。對(duì)關(guān)鍵質(zhì)量特性進(jìn)行質(zhì)量控制前,收集相關(guān)關(guān)鍵工序同組內(nèi)的歷史批次質(zhì)量數(shù)據(jù),作為備選先驗(yàn)信息。

1.2 數(shù)據(jù)變換

由于各批次歷史數(shù)據(jù)來(lái)自不同品種的產(chǎn)品,不滿足“獨(dú)立同分布”要求中關(guān)于“服從同一個(gè)分布”的條件,所以需要對(duì)這些數(shù)據(jù)進(jìn)行變換,使各批次數(shù)據(jù)具有相同的分布。

進(jìn)行數(shù)據(jù)變換的準(zhǔn)則有兩點(diǎn):①消除不同過(guò)程均值所帶來(lái)的影響;②消除不同的公差帶來(lái)的影響?；谧儞Q準(zhǔn)則,可以確定數(shù)據(jù)變換方式為

(1)

2 貝葉斯建模

2.1 制作控制圖

隨著生產(chǎn)的進(jìn)行,產(chǎn)生了n個(gè)質(zhì)量特征數(shù)據(jù),經(jīng)過(guò)數(shù)據(jù)變換后為{y1,y2,…,yn}=Y,假設(shè)Y的總體分布為正態(tài)分布N(μ,σ2),則控制圖控制限為

(2)

由式(2)可知,對(duì)質(zhì)量控制的效果主要取決于對(duì)超參數(shù)μ、σ估計(jì)的準(zhǔn)確程度。利用貝葉斯方法進(jìn)行動(dòng)態(tài)建模能夠充分利用先驗(yàn)信息,隨著生產(chǎn)過(guò)程的進(jìn)行,逐漸提高對(duì)μ、σ估計(jì)的準(zhǔn)確度,從而達(dá)到更好的質(zhì)量控制效果。

(3)

式中:nj為第j批歷史數(shù)據(jù)樣本容量。

由于多品種小批量生產(chǎn)條件多變,加工時(shí)間更接近當(dāng)前生產(chǎn)的歷史批次應(yīng)當(dāng)更能反映當(dāng)前的生產(chǎn)狀況,所以根據(jù)時(shí)間序列對(duì)各批次歷史數(shù)據(jù)進(jìn)行指數(shù)加權(quán)。

由于多品種小批量生產(chǎn)的特性,各企業(yè)所采取的檢驗(yàn)方式(全檢或抽樣檢驗(yàn))和抽樣方法不同,各批次數(shù)據(jù)量也可能不相同,故求取組間均值和組間方差時(shí)使用變樣本容量的計(jì)算公式。

上述兩方面相結(jié)合,可得歷史數(shù)據(jù)的組間均值與組間方差為

(4)

式中:k為加權(quán)系數(shù),一般取為(N+14)/14,N為平均各批次的生產(chǎn)天數(shù);m為先驗(yàn)信息的歷史數(shù)據(jù)總批次。

(2)求μ、σ2估計(jì)值。由于μ、σ2為正態(tài)分布的超參數(shù),根據(jù)共軛貝葉斯理論[16],可設(shè)σ2的先驗(yàn)分布為逆Gamma分布,記為IGa(α,β),μ的先驗(yàn)分布為正態(tài)分布,考慮到μ和σ2間有相互影響,故其共軛先驗(yàn)分布有乘積形式π(μ|σ2)π(σ2),其中

μ|σ2～N(a,b2);σ2～I(xiàn)Ga(α,β)

①估算σ2。對(duì)于σ2～I(xiàn)Ga(α,β),根據(jù)逆Gamma分布的性質(zhì),其期望和方差分別為

(5)

(6)

聯(lián)立(5)(6)兩式,可得

在得到當(dāng)前批數(shù)據(jù)Y后,根據(jù)共軛分布性質(zhì)可知,后驗(yàn)分布形式與先驗(yàn)分布相同,故可設(shè)

(7)

因此,可得σ2的貝葉斯估計(jì)為

(8)

②估算μ。對(duì)于μ|σ2～N(a,b2),由正態(tài)分布的性質(zhì),可得期望為

在得到當(dāng)前批數(shù)據(jù)Y后,根據(jù)共軛分布性質(zhì)可知

μ|σ2,Y～N(a′,b′2)

由共軛貝葉斯,可得(μ,σ2)的聯(lián)合先驗(yàn)密度函數(shù)為

π(μ,σ2)∝

因此,μ的邊緣后驗(yàn)密度

這是自由度為(2α+n)的t分布,其中

則可得μ的貝葉斯估計(jì)為

μB=E(μ|Y;μ1,μ2,…,μm)=a′=

(9)

(3)計(jì)算控制限。由于多品種小批量生產(chǎn)的數(shù)據(jù)量小,選用單值-移動(dòng)極差(X-MR)控制圖,結(jié)合實(shí)時(shí)生產(chǎn)數(shù)據(jù),使用式(8)(9),計(jì)算σ、μ的貝葉斯估計(jì)值,求單值圖的控制限

移動(dòng)極差的控制限為

(10)

(11)

2.2 工序能力分析

若當(dāng)前生產(chǎn)批次處于受控狀態(tài),應(yīng)該計(jì)算工序能力指數(shù),來(lái)判斷工序能力是否滿足要求。常用(雙邊公差)的工序能力指數(shù)計(jì)算公式如下

(12)

(13)

式中:Cp為分布中心與公差中心重合時(shí)的無(wú)偏移工序能力指數(shù);Cpk為有偏移工序能力指數(shù);μx、σx分別為原始數(shù)據(jù)的均值和標(biāo)準(zhǔn)差。根據(jù)數(shù)據(jù)變換式(1),可知

(14)

將μ=μB,σ=σB代入式(14),然后結(jié)合式(13)可得工序能力指數(shù)計(jì)算式如下

然后,根據(jù)工序能力評(píng)價(jià)標(biāo)準(zhǔn)[17]進(jìn)行評(píng)判。

3 仿 真

以某柴油機(jī)加工車間的鏜噴油器孔工序?yàn)槔?利用仿真生成數(shù)據(jù),比較本文方法與現(xiàn)有常用方法的優(yōu)劣。

表1 歷史質(zhì)量特征數(shù)據(jù)

3.1 分布參數(shù)估計(jì)誤差分析

應(yīng)用休哈特控制圖進(jìn)行質(zhì)量控制的效果主要取決于對(duì)分布參數(shù)估計(jì)的準(zhǔn)確程度,故通過(guò)仿真比較本文方法與兩種現(xiàn)有常用方法對(duì)參數(shù)μ、σ的估計(jì)情況,說(shuō)明本文方法的可行性。

針對(duì)上述柴油機(jī)加工車間的質(zhì)量控制問(wèn)題,兩種現(xiàn)有常用方法的處理過(guò)程如下。

方法1:通過(guò)成組技術(shù)、數(shù)據(jù)變換和假設(shè)檢驗(yàn)方法,利用多品種的歷史批次擴(kuò)充數(shù)據(jù)量,每批歷史數(shù)據(jù)取5個(gè),計(jì)算各批次歷史數(shù)據(jù)均值及標(biāo)準(zhǔn)差的平均值作為參數(shù)μ、σ的估計(jì)值,然后計(jì)算控制限,使用此固定控制限對(duì)工序質(zhì)量進(jìn)行控制。

方法2:利用當(dāng)前生產(chǎn)數(shù)據(jù)的均值與平均移動(dòng)極差,計(jì)算μ、σ的估計(jì)值,隨著生產(chǎn)的進(jìn)行實(shí)時(shí)計(jì)算控制限。

具體仿真過(guò)程如下:

(2)根據(jù)表1,從母體取9批數(shù)據(jù)作為歷史數(shù)據(jù);

(3)根據(jù)歷史數(shù)據(jù),按照方法1計(jì)算均值和標(biāo)準(zhǔn)差的估計(jì)值μm1、σm1,求估計(jì)誤差的絕對(duì)值(如均值估計(jì)誤差eμ,m1=|μm1-0|);

(4)根據(jù)當(dāng)前批數(shù)據(jù),以方法2計(jì)算均值和標(biāo)準(zhǔn)差估計(jì)值μm2、σm2,然后求估計(jì)誤差的絕對(duì)值;

(5)根據(jù)歷史數(shù)據(jù)和當(dāng)前批數(shù)據(jù),以本文方法計(jì)算μB、σB,求估計(jì)誤差的絕對(duì)值;

(6)從步驟(1)至(5)記為完成一次仿真,一共進(jìn)行1 000次仿真,求當(dāng)前批樣本容量為n時(shí)的平均參數(shù)估計(jì)誤差;

(7)令當(dāng)前批樣本容量n=n+1,返回步驟(1)。

隨著當(dāng)前批樣本容量的增加,3種方法對(duì)參數(shù)μ、σ的估計(jì)值平均誤差的變化見(jiàn)圖1。從圖1可以看出:方法1在生產(chǎn)初期誤差較小,在當(dāng)前批樣本容量較大時(shí)適用性較差;方法2只適用于當(dāng)前批樣本容量較大的情況;本文方法在整個(gè)生產(chǎn)階段都有更好表現(xiàn),并且隨著當(dāng)前批樣本容量的增加,估計(jì)誤差逐漸減小,較兩種現(xiàn)有常用方法有明顯優(yōu)勢(shì)。

(a)均值誤差

(b)標(biāo)準(zhǔn)差誤差圖1 分布參數(shù)估計(jì)平均誤差

由于方法1的準(zhǔn)確性依賴于歷史數(shù)據(jù)的可靠程度,當(dāng)歷史批次數(shù)據(jù)量減少或波動(dòng)加劇時(shí),方法1的控制效果顯著下降。由于多品種小批量生產(chǎn)的復(fù)雜特性可以預(yù)期,在實(shí)際生產(chǎn)中,本文方法的優(yōu)越性會(huì)更加明顯。

3.2 平均運(yùn)行鏈長(zhǎng)

平均運(yùn)行鏈長(zhǎng)(ARL)是評(píng)價(jià)控制圖性能的重要指標(biāo)[18]。生產(chǎn)過(guò)程處于穩(wěn)態(tài)時(shí)的平均運(yùn)行鏈長(zhǎng)記為L(zhǎng)ARL,0,生產(chǎn)過(guò)程發(fā)生異常波動(dòng)時(shí)的平均運(yùn)行鏈長(zhǎng)記為L(zhǎng)ARL,1。LARL,0越大表示控制圖誤報(bào)率越低,LARL,1越小表示控制圖的漏報(bào)率越低,即異常檢出力越強(qiáng)。

據(jù)研究表明[19],當(dāng)累積和(CUSUM)控制圖的參數(shù)k=0.5、h=4.77時(shí),其受控平均運(yùn)行鏈長(zhǎng)LARL,0與常用的休哈特控制圖相當(dāng)。因此,在該參數(shù)情況下,比較累積和控制圖方法與其他休哈特控制圖相關(guān)方法的LARL,1,以比較各方法的異常檢出力。

方法3(累積和控制圖):通過(guò)成組技術(shù)、數(shù)據(jù)變換和假設(shè)檢驗(yàn)方法,利用多品種的歷史批次擴(kuò)充數(shù)據(jù)量,每批歷史數(shù)據(jù)取5個(gè),計(jì)算各批次歷史數(shù)據(jù)標(biāo)準(zhǔn)差的平均值作為分布參數(shù)σ的估計(jì)值σC,作累積和控制圖,其中參數(shù)K=0.5σC,H=4.77σC。

利用仿真方法估計(jì)從不同位置開(kāi)始發(fā)生波動(dòng)時(shí)本文方法及其他3種現(xiàn)有常用方法的LARL,1。仿真過(guò)程如下:

(3)分別應(yīng)用3種常用方法和本文方法進(jìn)行控制,計(jì)算從開(kāi)始發(fā)生波動(dòng)到控制圖報(bào)警之間的運(yùn)行鏈長(zhǎng);

(4)從步驟(1)至(3)記為完成一次仿真,一共進(jìn)行1 000次仿真;

(5)計(jì)算各方法的LARL,1;

(6)令p=p+1,返回步驟(1)。

表2列出了隨著發(fā)生波動(dòng)位置的變化,各種方法LARL,1的變化。

從表2可以看出,當(dāng)受控生產(chǎn)數(shù)據(jù)較少時(shí),方法1的LARL,1與本文方法接近,均小于方法2與方法3。隨著受控生產(chǎn)數(shù)據(jù)的增加,本文方法的優(yōu)勢(shì)漸漸呈現(xiàn),LARL,1明顯小于其他3種方法。因此,當(dāng)標(biāo)準(zhǔn)差發(fā)生波動(dòng)時(shí),本文方法有更強(qiáng)的異常檢出力。

表2 幾種方法平均運(yùn)行鏈長(zhǎng)LARL,1的比較

4 結(jié) 論

結(jié)合成組技術(shù)、指數(shù)加權(quán)及貝葉斯理論的優(yōu)點(diǎn),提出了一種動(dòng)態(tài)質(zhì)量控制方法。首先使用成組技術(shù)和數(shù)據(jù)變換方法擴(kuò)充數(shù)據(jù)量作為先驗(yàn)信息;然后結(jié)合指數(shù)加權(quán)思想及貝葉斯理論建模,并計(jì)算工序能力指數(shù);最后通過(guò)仿真驗(yàn)證了該方法的有效性。與現(xiàn)有常用方法相比,本文方法在適用范圍、控制效果、異常檢出力等方面均有提高,一定程度上解決了多品種小批量生產(chǎn)模式下的質(zhì)量控制問(wèn)題。在后期工作中,可以進(jìn)一步研究抽樣方法和分組準(zhǔn)則等因素對(duì)質(zhì)量控制的具體影響,進(jìn)一步提高多品種小批量生產(chǎn)的質(zhì)量控制效果。