王維志

（湖北汽車工業(yè)學(xué)院 機(jī)械工程系，湖北 十堰442002）

在機(jī)械加工制造領(lǐng)域，人們一直用過程能力指數(shù)Cpk評(píng)價(jià)生產(chǎn)過程滿足產(chǎn)品技術(shù)規(guī)格要求的能力［1－2］。 這里提出一個(gè)問題，Cpk評(píng)價(jià)的是不是生產(chǎn)過程持續(xù)滿足產(chǎn)品技術(shù)規(guī)格要求的能力，結(jié)論顯然是否定的。只有當(dāng)生產(chǎn)過程是穩(wěn)定的，或者說是始終處于統(tǒng)計(jì)穩(wěn)態(tài)，結(jié)論才是肯定的，而Cpk本身并不能保證這一點(diǎn)。因此，需要另外一種方法，計(jì)算Cpk之前確認(rèn)生產(chǎn)過程是統(tǒng)計(jì)穩(wěn)態(tài)的，即計(jì)算所使用的數(shù)據(jù)均來自同一個(gè)正態(tài)總體。一個(gè)生產(chǎn)過程是不是會(huì)持續(xù)保持穩(wěn)定，回答也是否定的，特別對(duì)于機(jī)械加工類的生產(chǎn)過程，刀具的磨損、不均勻的熱變形等因素會(huì)很快導(dǎo)致過程的失穩(wěn) （統(tǒng)計(jì)不再穩(wěn)定）。因此，現(xiàn)代生產(chǎn)質(zhì)量管理應(yīng)當(dāng)實(shí)時(shí)監(jiān)測(cè)過程的運(yùn)行，當(dāng)過程出現(xiàn)失穩(wěn)或有失穩(wěn)趨勢(shì)發(fā)生時(shí)，及時(shí)發(fā)出報(bào)警，通過調(diào)節(jié)使過程恢復(fù)統(tǒng)計(jì)穩(wěn)態(tài)。

與逐件檢驗(yàn)產(chǎn)品質(zhì)量的控制方法不同，過程統(tǒng)計(jì)控制研究的目的在于保持過程處于統(tǒng)計(jì)穩(wěn)態(tài)和及時(shí)發(fā)現(xiàn)過程運(yùn)行的變化趨勢(shì)，從而防止不合格產(chǎn)品的出現(xiàn)，實(shí)現(xiàn)預(yù)防為主的產(chǎn)品質(zhì)量管理目標(biāo)。

1 過程統(tǒng)計(jì)控制應(yīng)用現(xiàn)狀

目前生產(chǎn)管理過程中，也有一些推斷統(tǒng)計(jì)穩(wěn)態(tài)的方法［3－5］，比較常用的有 Shewhart控制圖、累積和控制圖 （CUSUM）和指數(shù)加權(quán)移動(dòng)平均控制圖（EWMA）等。其中Shewhart控制圖應(yīng)用最為普遍，計(jì)算簡(jiǎn)單，而且有許多商品化的軟件可以應(yīng)用，缺點(diǎn)是研判準(zhǔn)則比較復(fù)雜，過程均值有小的漂移時(shí)，檢出效果比較差，對(duì)于樣本數(shù)據(jù)量比較少的研究性過程分析以及定量過程分析也都存在一定的不足。累積和控制圖（CUSUM）和指數(shù)加權(quán)移動(dòng)平均控制圖（EWMA）等統(tǒng)計(jì)分析工具，雖然它們已被證明在檢 測(cè) 小 的 漂 移 時(shí) 效 果 不 錯(cuò) 外［5－7］。 也 都 存 在 和Shewhart控制圖同樣的其它缺點(diǎn)。

隨著產(chǎn)品特點(diǎn)的不斷變化，例如質(zhì)量要求不斷提高，產(chǎn)品的批次增加，單批次的數(shù)量減少等原因，對(duì)生產(chǎn)過程統(tǒng)計(jì)特性分析也提出新的要求。例如在樣本數(shù)據(jù)量少的情況下，對(duì)較小的漂移要有比較高的檢出精度，統(tǒng)計(jì)穩(wěn)態(tài)的的研判更加的簡(jiǎn)單以及有量化的指標(biāo)等。

方差分析（Analysis of Variance，簡(jiǎn)稱ANOVA）是從觀測(cè)變量的方差入手，研究諸多控制變量中哪些變量是對(duì)觀測(cè)變量有顯著影響的變量。理論上講，控制變量可以是主動(dòng)施加的，也可以是潛在存在的。只要對(duì)觀測(cè)變量產(chǎn)生了顯著影響，就可以通過對(duì)觀測(cè)變量的方差分析找到它們。基于這樣一個(gè)原理，本文運(yùn)用ANOVA對(duì)生產(chǎn)過程進(jìn)行統(tǒng)計(jì)特性分析，建立模型，實(shí)現(xiàn)對(duì)生產(chǎn)過程的統(tǒng)計(jì)穩(wěn)態(tài)研判和連續(xù)統(tǒng)計(jì)控制。

2 磨削過程的統(tǒng)計(jì)特性分析

機(jī)械加工過程的輸出具有波動(dòng)性，導(dǎo)致波動(dòng)的原因有隨機(jī)的因素，也有非隨機(jī)的因素。完全的隨機(jī)因素引起的波動(dòng)具有正態(tài)性，分布形態(tài)為對(duì)稱結(jié)構(gòu)，范圍約為±3σ。加入非隨機(jī)因素的影響往往導(dǎo)致分布形態(tài)結(jié)構(gòu)發(fā)生改變，在實(shí)際生產(chǎn)過程中的后果可能是過程輸出超出產(chǎn)品技術(shù)規(guī)格的控制界，形成不合格品。理想的加工過程應(yīng)該力求避免非隨機(jī)影響因素的存在，事實(shí)上，完全排除也不現(xiàn)實(shí)，非隨機(jī)影響因素不顯著的情況下，可以忽略其影響，把過程的輸出看作是正態(tài)分布的，或者稱過程處于統(tǒng)計(jì)穩(wěn)態(tài)。當(dāng)非隨機(jī)因素影響顯著，過程輸出的分布發(fā)生漂移，就說生產(chǎn)過程處于統(tǒng)計(jì)不穩(wěn)態(tài)，一般來說，統(tǒng)計(jì)不穩(wěn)態(tài)的生產(chǎn)過程，其輸出超出產(chǎn)品技術(shù)規(guī)格的概率會(huì)增加，不穩(wěn)態(tài)的程度越大，超出技術(shù)規(guī)格的概率也就越大，出現(xiàn)廢品的概率也就越大。

對(duì)于大部分機(jī)械加工過程來說，在進(jìn)行正確調(diào)試和正確設(shè)置工藝參數(shù)的條件下都可以達(dá)到統(tǒng)計(jì)穩(wěn)態(tài)這樣一個(gè)要求［8－12］，不過是需要加以驗(yàn)證的。過程運(yùn)行中，統(tǒng)計(jì)穩(wěn)態(tài)只是相對(duì)的，由于不均勻熱變形、刀具磨損以及由此引起的振動(dòng)等原因會(huì)導(dǎo)致過程變得不再穩(wěn)態(tài)。

統(tǒng)計(jì)特性分析可以識(shí)別生產(chǎn)過程中的非隨機(jī)影響因素，并推斷其影響程度，一旦確認(rèn)影響比較顯著時(shí)，給出報(bào)警并可以采取措施加以消除（消除方法不在本文討論的范圍）。統(tǒng)計(jì)特性分析的優(yōu)勢(shì)還在于可以通過局部推斷總體，即利用有限的樣本推斷過程總體的變化趨勢(shì)。另外，統(tǒng)計(jì)特性分析也可以檢驗(yàn)生產(chǎn)過程是否正確設(shè)置了工藝和工藝結(jié)構(gòu)參數(shù)以及機(jī)床設(shè)備是否進(jìn)行了正確的調(diào)試。

3 磨削過程統(tǒng)計(jì)控制建模

磨削過程連續(xù)統(tǒng)計(jì)控制一般在完成初始統(tǒng)計(jì)分析之后進(jìn)行，初始化統(tǒng)計(jì)分析確認(rèn)磨削過程正確設(shè)置了工藝和工藝結(jié)構(gòu)參數(shù)，確認(rèn)機(jī)床設(shè)備進(jìn)行了正確的調(diào)試。另外通過初始化統(tǒng)計(jì)分析以后，獲得了過程的統(tǒng)計(jì)特征參數(shù)。

從一個(gè)待觀察的磨削過程中抽取樣本，全部的觀測(cè)值記作：

式中：k為樣本數(shù)；n為樣本內(nèi)獨(dú)立產(chǎn)品觀測(cè)值的個(gè)數(shù)。假設(shè)數(shù)據(jù)來自于統(tǒng)計(jì)穩(wěn)態(tài)過程，xij服從正態(tài)或近似正態(tài)分布 N（μ，σ2）。

記各觀測(cè)值xij與總平均數(shù)的離均差平方和為ST，記作：

變換式（1）為式（2）的形式：

展開式（2）并整理得

其中：

它們分別擁有的自由度為

根據(jù) χ2分布定理［13］，在樣本數(shù)據(jù)服從正態(tài)分布的條件下，和都是隨機(jī)變量，且分別服從自由度為k1和k2的χ2分布，同樣地它們的比值也是隨機(jī)變量，且服從自由度為k1、k2的F分布。式（8）是基于現(xiàn)場(chǎng)采樣數(shù)據(jù)構(gòu)建的初始化統(tǒng)計(jì)分析模型，其值大小表示不同樣本之間的差異顯著性，是一個(gè)統(tǒng)計(jì)穩(wěn)態(tài)的量化指標(biāo)，但是僅僅依據(jù)其大小還不能推斷其差異顯著程度（過程的統(tǒng)計(jì)穩(wěn)定性），還需要有一個(gè)比對(duì)的參照系—F分布來加以擬合［13］。

在概率論和統(tǒng)計(jì)學(xué)里，F(xiàn)分布是一種連續(xù)概率分布，被廣泛應(yīng)用于似然比率檢驗(yàn)，F(xiàn)分布密度曲線是隨自由度k1、k2的變化而變化的一簇偏態(tài)曲線，F(xiàn)α稱作右尾臨界值，對(duì)應(yīng)的 Fα至+∞的概率值稱為右尾臨界概率，記作α%，α稱作顯著性水平，做假設(shè)檢驗(yàn)時(shí)，一般先設(shè)定α值，在給定自由度（k1，k2）的條件下，通過數(shù)理統(tǒng)計(jì)工具資料查出右尾臨界值Fα。如果被觀察的過程輸出分布來自正態(tài)過程，式（8）計(jì)算出來的F值落在0至Fα區(qū)間的概率為（1－α）%，作假設(shè)檢驗(yàn)一般α取值比較小，F(xiàn)值落在Fα至+∞區(qū)間的概率比較小，稱作小概率事件，一旦發(fā)生，就推斷被觀測(cè)的數(shù)據(jù)不是來自正態(tài)過程，樣本之間有顯著的差異性，過程統(tǒng)計(jì)不穩(wěn)定，誤判的風(fēng)險(xiǎn)不大于α%。

實(shí)際的磨削過程僅作初始統(tǒng)計(jì)分析是不夠的，磨削過程進(jìn)行中，砂輪的磨損，熱變形、振動(dòng)等影響因素隨時(shí)會(huì)使磨削過程變得失穩(wěn)或有失穩(wěn)趨勢(shì)發(fā)生，持續(xù)的統(tǒng)計(jì)特性分析是必須的。過程投入正式運(yùn)行以后，需要建立新的統(tǒng)計(jì)模型并按適當(dāng)?shù)拈g隔采集樣本，推斷磨削過程是否持續(xù)保持在統(tǒng)計(jì)穩(wěn)態(tài)。

LSD法稱作最小差數(shù)檢驗(yàn)法，是在F檢驗(yàn)的基礎(chǔ)上實(shí)現(xiàn)連續(xù)檢驗(yàn)的一種方法，它是將t檢驗(yàn)中由所求得的t之絕對(duì)值

統(tǒng)計(jì)穩(wěn)態(tài)條件下的最小顯著差數(shù)LSDα則是：

初始統(tǒng)計(jì)分析為評(píng)價(jià)后續(xù)過程是否統(tǒng)計(jì)穩(wěn)態(tài)提供了基準(zhǔn)參考依據(jù)。前承后接，有量化的統(tǒng)一標(biāo)準(zhǔn)，有較高的可靠性。LSD法又不同于每次利用2組數(shù)據(jù)進(jìn)行多個(gè)平均數(shù)兩兩比較的t檢驗(yàn)法。它解決了t檢驗(yàn)法檢驗(yàn)過程煩瑣，無統(tǒng)一的試驗(yàn)誤差且估計(jì)誤差的精確性和檢驗(yàn)的靈敏性低的問題。

4 實(shí)例驗(yàn)證

以某型號(hào)柴油發(fā)動(dòng)機(jī)零件氣門挺桿的磨削加工為例，設(shè)備為無心磨床；工藝參數(shù)略；生產(chǎn)類型為大批量生產(chǎn)；工件表面粗糙度規(guī)格要求為Ra不大于0.25μm。統(tǒng)計(jì)穩(wěn)態(tài)驗(yàn)證同時(shí)選擇工件尺寸和表面粗糙度值作為特征觀測(cè)值進(jìn)行研究，限于篇幅，本文僅列出表面粗糙度部分的原始數(shù)據(jù)和結(jié)論。需要說明的是，實(shí)際的生產(chǎn)過程也不是每一個(gè)工件尺寸都要進(jìn)行統(tǒng)計(jì)過程分析，關(guān)聯(lián)密切的尺寸只分析其中部分即可。產(chǎn)品投入正式生產(chǎn)之前，進(jìn)行生產(chǎn)過程的統(tǒng)計(jì)穩(wěn)態(tài)分析，正式生產(chǎn)后進(jìn)行連續(xù)統(tǒng)計(jì)穩(wěn)態(tài)控制。原始樣本數(shù)據(jù)完全取自企業(yè)的生產(chǎn)現(xiàn)場(chǎng)，見表1。取α分別為0.01和0.05，初始化統(tǒng)計(jì)分析結(jié)論見表2。為計(jì)算方便起見，所有的誤差數(shù)據(jù)放大一個(gè)相同的倍數(shù)。

按照一定的時(shí)間間隔連續(xù)的從后續(xù)磨削過程抽取樣本，采樣數(shù)據(jù)見表3，并按照LSD方法推斷過程的持續(xù)穩(wěn)定性。

表1 原始采樣數(shù)據(jù)

表2 初始化分析數(shù)據(jù)

表3 連續(xù)的樣本數(shù)據(jù)

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

通過上述理論推導(dǎo)和實(shí)例驗(yàn)證分析，可以歸納得出：應(yīng)用F檢驗(yàn)進(jìn)行磨削過程統(tǒng)計(jì)穩(wěn)態(tài)的假設(shè)檢驗(yàn)，事實(shí)上是一種分布假設(shè)檢驗(yàn)，以檢驗(yàn)均值的變化為主，相較于傳統(tǒng)數(shù)理統(tǒng)計(jì)理論中使用 χ2檢驗(yàn)作為分布檢驗(yàn)，更符合磨削加工的實(shí)際情景。

傳統(tǒng)的LSD法主要用于多個(gè)處理的多重比較，本文改進(jìn)后用于連續(xù)的過程統(tǒng)計(jì)控制，實(shí)現(xiàn)過程均值變化的連續(xù)監(jiān)控。

相對(duì)于傳統(tǒng)的過程統(tǒng)計(jì)控制（SPC）方法，基于ANOVA法模型實(shí)現(xiàn)了統(tǒng)計(jì)穩(wěn)態(tài)的定量分析，簡(jiǎn)化了過程失穩(wěn)的判斷條件，使判斷過程更加的客觀、簡(jiǎn)單和有效，有利于在生產(chǎn)現(xiàn)場(chǎng)推廣應(yīng)用。

［1］ 宋祥彥.解讀過程能力指數(shù) 4個(gè)基本特性［J］.標(biāo)準(zhǔn)科學(xué)，2010（1）:33－39.

［2］ 韓振魯，李長(zhǎng)河.磨削加工建模與仿真的研究進(jìn)展［J］.磨具磨料通訊，2011（1）：7－19.

［3］ 于濤，王高山.統(tǒng)計(jì)過程控制在先進(jìn)制造環(huán)境中的內(nèi)涵和作用［J］.工業(yè)工程與管理，2005（3）：51－54.

［4］ 陳健.統(tǒng)計(jì)過程控制（SPC）在產(chǎn)品質(zhì)量控制中的應(yīng)用研究［D］.南京：東南大學(xué)，2006.

［5］ 田學(xué)民，曹玉蘋.統(tǒng)計(jì)過程控制的研究現(xiàn)狀及展望［J］.石油大學(xué)學(xué)報(bào)：自然科學(xué)版，2008，32（4）：175－180.

［6］ 濮曉龍.關(guān)于累積和（CUSUM）檢驗(yàn)的改進(jìn)［J］.應(yīng)用數(shù)學(xué)學(xué)報(bào)，2003，26（2）：225－241.

［7］ 董其義.導(dǎo)彈測(cè)試系統(tǒng)統(tǒng)計(jì)過程控制方法及其實(shí)現(xiàn)［J］.海軍航空工程學(xué)院學(xué)報(bào)，2009，24（6）：666－667.

［8］ 錢夕元，荊建芬，侯旭珊.統(tǒng)計(jì)過程控制(SPC)及其應(yīng)用研究［J］.計(jì)算機(jī)工程，2004，30（19）：144－154.

［9］ 張磊，徐曉輝.砂輪動(dòng)態(tài)磨粒的概率統(tǒng)計(jì)數(shù)學(xué)模型［J］.應(yīng)用科技，2010，37（10）：17－20.

［10］ 曾偉，無心磨加工機(jī)理及工藝參數(shù)優(yōu)化研究［D］.武漢：華中科技大學(xué)，2006.

［11］ 張公緒.新編質(zhì)量管理學(xué)［M］.1版.北京：高等教育出版社，1998.

［12］ 龐麗君，尚曉峰.金屬切削原理［M］.1版.北京：國(guó)防工業(yè)出版社，2009.

［13］ 沈恒范.概率論與數(shù)理統(tǒng)計(jì)教程［M］.4版.北京：高等教育出版社,2003.

［14］修世超，李長(zhǎng)河，蔡光起.磨削加工表面粗糙度理論模型修正方法［J］.東北大學(xué)學(xué)報(bào)：自然科學(xué)版，2005，26（8）：770－773.

［15］ 王維志.無心磨削過程統(tǒng)計(jì)特性研究［J］.湖北汽車工業(yè)學(xué)院學(xué)報(bào)，2011，25（2）：35－38.