亚洲免费av电影一区二区三区,日韩爱爱视频,51精品视频一区二区三区,91视频爱爱,日韩欧美在线播放视频,中文字幕少妇AV,亚洲电影中文字幕,久久久久亚洲av成人网址,久久综合视频网站,国产在线不卡免费播放

        ?

        基于R語(yǔ)言的制造質(zhì)量控制與預(yù)測(cè)*

        2019-03-01 09:37:08喬仁杰宋庭新
        關(guān)鍵詞:花鍵軸哈特波動(dòng)

        喬仁杰,宋庭新

        (湖北工業(yè)大學(xué) 機(jī)械工程學(xué)院,武漢 430068)

        0 引言

        在現(xiàn)代制造企業(yè)中,產(chǎn)品制造過(guò)程中的質(zhì)量控制無(wú)疑是企業(yè)最核心的競(jìng)爭(zhēng)力之一。1924年休哈特(Hugh Hart)提出的控制圖理論雖然能夠很直觀(guān)地反映質(zhì)量數(shù)據(jù)的波動(dòng)情況,但卻忽略了歷史數(shù)據(jù)對(duì)當(dāng)前數(shù)據(jù)的累積影響,對(duì)于過(guò)程的小偏移檢出力較弱。所以,采取聯(lián)合休哈特累積和(CUSUM)控制圖,增加對(duì)觀(guān)測(cè)值與目標(biāo)值之差的累積和來(lái)描點(diǎn),利用累積數(shù)據(jù),可以提升對(duì)質(zhì)量狀態(tài)小偏移的檢出能力[1-2]。近年來(lái),許多學(xué)者對(duì)質(zhì)量控制理論進(jìn)行了大量研究,如龔立雄等利用信息化技術(shù)設(shè)計(jì)了一套完整的質(zhì)量信息系統(tǒng)來(lái)對(duì)質(zhì)量進(jìn)行實(shí)時(shí)監(jiān)控與分析[3],很好地解決了實(shí)時(shí)質(zhì)量數(shù)據(jù)監(jiān)控,但不能對(duì)未來(lái)質(zhì)量數(shù)據(jù)進(jìn)行預(yù)測(cè)與分析。江宇平等提出的一種基于賦值型誤差傳遞網(wǎng)絡(luò)的多工序加工質(zhì)量預(yù)測(cè)建模方法,能有效預(yù)測(cè)加工誤差[4],但這種預(yù)測(cè)方法主要針對(duì)加工過(guò)程中的誤差,且過(guò)程較繁瑣,難以普及。R語(yǔ)言作為一種專(zhuān)門(mén)為數(shù)理統(tǒng)計(jì)、數(shù)據(jù)分析和統(tǒng)計(jì)制圖開(kāi)發(fā)的開(kāi)源編程語(yǔ)言,可以快速生成質(zhì)量控制圖并對(duì)質(zhì)量數(shù)據(jù)進(jìn)行分析與預(yù)測(cè)。本文結(jié)合汽車(chē)零部件制造中花鍵軸的實(shí)例,利用R語(yǔ)言繪制聯(lián)合休哈特累積和控制圖并針對(duì)花鍵軸特點(diǎn)建模,選取Holt-Winters質(zhì)量預(yù)測(cè)模型來(lái)對(duì)未來(lái)質(zhì)量數(shù)據(jù)進(jìn)行預(yù)測(cè)。這種方法可以很好地簡(jiǎn)化質(zhì)量預(yù)測(cè)過(guò)程,提高質(zhì)量預(yù)測(cè)效率,為解決質(zhì)量預(yù)測(cè)難以在中小型企業(yè)中普及的問(wèn)題提供了有效的思路。

        1 質(zhì)量控制理論

        1.1 休哈特控制理論

        休哈特控制理論主要通過(guò)采集質(zhì)量數(shù)據(jù)樣本,計(jì)算標(biāo)準(zhǔn)差σ和管理上限UCL及下限LCL,當(dāng)加工工序處于穩(wěn)定狀態(tài)時(shí),隨機(jī)誤差具有一定的分布規(guī)律,總體上服從正態(tài)分布N(μ,σ)或近似正態(tài)分布。在正態(tài)分布的±3σ范圍內(nèi),樣品特征值出現(xiàn)在μ+3σ與μ-3σ上下限之間的概率為99.73%,超出該部分的概率僅為0.27%。其計(jì)算步驟如下:

        (1)

        (2)求樣品極差R或者標(biāo)準(zhǔn)差S,其中S代表標(biāo)準(zhǔn)差σ。

        (2)

        (3)確定控制界限。一般控制界限選取為3σ標(biāo)準(zhǔn)差,即:

        (3)

        按照計(jì)算的參數(shù)繪制休哈特控制圖(見(jiàn)圖1),圖中折線(xiàn)代表產(chǎn)品質(zhì)量的波動(dòng),在管制上限UCL與管制下限LCL內(nèi)的波動(dòng)屬于可控狀態(tài),當(dāng)超出則說(shuō)明失控。

        圖1 休哈特控制圖

        1.2 累積和控制理論

        累積和控制圖是基于與歷史數(shù)據(jù)比較的結(jié)果,將多次數(shù)據(jù)的波動(dòng)進(jìn)行整合來(lái)放大波動(dòng)效果,使累積和控制圖對(duì)于均值微小的偏移過(guò)程更加敏感。其一般步驟如下:

        (1)求樣本標(biāo)準(zhǔn)偏差σLT;

        (4)

        (2)求累積和圖,其中Ku、Kl表示上限值和下限值,上偏差mszu與下偏差mszl參數(shù)(一般取1σ);

        (5)

        (6)

        (7)

        繪制的累積和控制圖中有上累積線(xiàn)與下累積線(xiàn)兩條折線(xiàn),兩條不同的繪圖點(diǎn)讓累積和控制圖可以同時(shí)看到累積和上下偏移目標(biāo)的情況,刻畫(huà)了微小偏移的波動(dòng)(見(jiàn)圖2)。

        圖2 累積和控制圖

        通過(guò)累積和控制圖發(fā)現(xiàn)大約在10個(gè)數(shù)據(jù)前上累積和與下累積和基本都趨于零,說(shuō)明此時(shí)數(shù)據(jù)比較平穩(wěn),而之后上累積大幅的向上偏移并超過(guò)了UCL控制上線(xiàn),此時(shí)說(shuō)明數(shù)據(jù)整體向上偏移異常,系統(tǒng)處于失控狀態(tài)。所以,通過(guò)聯(lián)合休哈特與累積和控制圖既可以直觀(guān)地反應(yīng)數(shù)據(jù)波動(dòng)情況,又可以刻畫(huà)歷史數(shù)據(jù)的波動(dòng)影響及對(duì)微小波動(dòng)的敏感性,提高質(zhì)量控制的效率及準(zhǔn)確率。

        2 案例研究

        在重型汽車(chē)車(chē)橋中,傳動(dòng)軸聯(lián)結(jié)一般采用漸開(kāi)線(xiàn)花鍵聯(lián)結(jié),因此花鍵的加工質(zhì)量非常關(guān)鍵[5]。某汽車(chē)零部件企業(yè)生產(chǎn)的花鍵軸產(chǎn)品如圖3所示,其外花鍵測(cè)量棒間距均值要求為μ0=45.550mm,測(cè)量20組外花鍵跨棒距參數(shù)D作為統(tǒng)計(jì)樣本。

        圖3 花鍵軸與外花鍵跨棒距D

        取樣本量為20,測(cè)得質(zhì)量數(shù)據(jù)矩陣為:

        在R語(yǔ)言的集成開(kāi)發(fā)環(huán)境RStudio中進(jìn)行編程和統(tǒng)計(jì)制圖。使用R語(yǔ)言的sd()函數(shù)求出數(shù)據(jù)矩陣中的標(biāo)準(zhǔn)差σ = 0.00590517,采用3σ的管控系數(shù),則管控中心線(xiàn)CL =μ0= 45.550 mm,上、下控制限為:UCL = CL+3σ = 45.53228 mm,LCL = CL-3σ = 45.56772 mm。再通過(guò)plot()函數(shù)將收集的樣本數(shù)據(jù)進(jìn)行描點(diǎn)及繪制,并附加abline()函數(shù)繪制控制圖的管控線(xiàn)UCL和LCL,繪制的控制圖如圖4所示。

        圖4 在RStudio中繪制休哈特控制圖

        圖4中的數(shù)據(jù)均處于管控上線(xiàn)UCL與管控下線(xiàn)LCL的范圍之中,說(shuō)明此時(shí)數(shù)據(jù)沒(méi)有明顯超出控制,但從圖4中發(fā)現(xiàn)數(shù)據(jù)存在向上偏移的趨勢(shì)。因此聯(lián)合使用累積和CUSUM控制,如圖5所示。

        圖5 聯(lián)合累積和CUSUM控制圖

        從圖5中可以發(fā)現(xiàn)11號(hào)樣本之前的累積上下線(xiàn)都處于0左右,說(shuō)明數(shù)據(jù)沒(méi)有明顯的偏移,而11號(hào)之后下累積線(xiàn)一直為0,上累積線(xiàn)大幅度向上偏移,說(shuō)明此時(shí)數(shù)據(jù)整體向上偏移,從第19個(gè)數(shù)據(jù)開(kāi)始,數(shù)據(jù)整體向上偏移超過(guò)控制上線(xiàn),說(shuō)明此時(shí)質(zhì)量數(shù)據(jù)處于失控狀態(tài),需要重新檢查系統(tǒng)找出問(wèn)題原因,之后再進(jìn)行系統(tǒng)評(píng)估。

        3 質(zhì)量預(yù)測(cè)

        3.1 Holt-Winters指數(shù)平滑法

        目前,在質(zhì)量預(yù)測(cè)中最常用的時(shí)間序列分析方法是指數(shù)平滑法,該方法通過(guò)對(duì)不同時(shí)間序列數(shù)據(jù)的權(quán)重分配以及從近期到遠(yuǎn)期數(shù)據(jù)逐漸賦予收斂至零的權(quán)數(shù),來(lái)合理利用各個(gè)時(shí)間段上的數(shù)據(jù),從而對(duì)未來(lái)質(zhì)量數(shù)據(jù)進(jìn)行預(yù)測(cè)[6-8]。不同次數(shù)的指數(shù)平滑對(duì)應(yīng)不同預(yù)測(cè)模型。一次指數(shù)平滑適用于無(wú)明顯趨勢(shì)和季節(jié)性的時(shí)間序列;二次指數(shù)平滑在第一次的基礎(chǔ)上再進(jìn)行指數(shù)平滑,增加了趨勢(shì)參數(shù)但沒(méi)有季節(jié)性的序列;三次指數(shù)平滑在兩次指數(shù)平滑算法上增加季節(jié)參數(shù)P來(lái)對(duì)趨勢(shì)進(jìn)行預(yù)測(cè),更能體現(xiàn)周期性變化的規(guī)律。三次指數(shù)平滑也叫Holt-Winters指數(shù)平滑法。Holt-Winters指數(shù)平滑法分為累加和累乘兩種,累加是根據(jù)線(xiàn)性時(shí)間趨勢(shì)與加法模型序列,累乘是根據(jù)線(xiàn)性時(shí)間趨勢(shì)與乘法模型序列。由于不同季度不合格品是個(gè)數(shù)的對(duì)比,比如夏季比冬季的不合格品數(shù)多30件,所以選用累加模型。而累乘一般適用于大量數(shù)據(jù)的百分比比較[9-10]。

        基于累加的Holt-Winters指數(shù)平滑公式為:

        si=α(xi-pi-k)+(1-α)(si-1-ti-1)

        (8)

        ti=β(si-si-1)+(1-β)ti-1

        (9)

        pi=γ(xi-si)+(1-γ)pi-k(k=周期)

        (10)

        其累加Holt-Winter預(yù)測(cè)公式為:

        xi+h=si+hti+pi-k+(hmodk)

        (11)

        其中,公式(8)中的α是平滑參數(shù),t用來(lái)保留平滑趨勢(shì),p為季節(jié)參數(shù),si則是i個(gè)數(shù)據(jù)的平滑值,取值為[0,1]。α越接近0,說(shuō)明數(shù)據(jù)平滑后越接近i個(gè)數(shù)據(jù)的平均值,數(shù)據(jù)越平滑;α越接近1說(shuō)明平滑后的值越接近當(dāng)前時(shí)間的數(shù)據(jù)值,數(shù)據(jù)越不平滑。α的值通常通過(guò)具體情況進(jìn)行嘗試來(lái)確定最佳值。在式(8)~式(11)中,α,β,γ的值均應(yīng)該多次測(cè)試與實(shí)驗(yàn),使其最佳值位于[0,1]之間。s,t,p的初始值的選取對(duì)于整體算法的影響較小,通常取s0=x0,t0=x1-x1,在累加中p值默認(rèn)為0[11-12]。

        3.2 Holt-Winters預(yù)測(cè)過(guò)程

        在R語(yǔ)言中預(yù)置有Holt-Winters模型包,通過(guò)調(diào)用包內(nèi)的函數(shù)并配置相關(guān)參數(shù)就能建立Holt-Winters模型。將采集到的上述花鍵軸跨棒距數(shù)據(jù)D導(dǎo)入RStudio編譯器中,利用R語(yǔ)言的ts()函數(shù)對(duì)所取數(shù)據(jù)進(jìn)行時(shí)間序列化,取10為底的log對(duì)數(shù),命名為測(cè)量距d,增加數(shù)據(jù)波動(dòng)的敏感性。再調(diào)用plot()函數(shù)打印時(shí)間序列化后的數(shù)據(jù)圖如圖6所示。

        圖6 時(shí)間序列化圖

        從圖6可以看出,經(jīng)過(guò)時(shí)間序列化后的圖形趨勢(shì)變化更加明顯。時(shí)間序列t表示在時(shí)間序列下的新坐標(biāo),1.0表示第一個(gè)時(shí)間點(diǎn)獲取的質(zhì)量數(shù)據(jù),隨后表示每經(jīng)過(guò)0.1個(gè)時(shí)間段獲取的新的數(shù)據(jù),這樣20個(gè)樣本數(shù)據(jù)分別與時(shí)間序列1.0~2.9之間的時(shí)間點(diǎn)對(duì)應(yīng)。

        接著建立Holt-Winters模型,利用R語(yǔ)言中的Holt-Winters()函數(shù),將時(shí)間序列轉(zhuǎn)化后的數(shù)據(jù)建模,將數(shù)據(jù)過(guò)濾為Holt-Winters模型下的趨勢(shì)圖(見(jiàn)圖7),其中無(wú)圓點(diǎn)的折線(xiàn)表示原數(shù)據(jù)經(jīng)過(guò)Holt-Winters三次指數(shù)平滑計(jì)算后的觀(guān)察值,而圓點(diǎn)折線(xiàn)表示對(duì)現(xiàn)有數(shù)據(jù)進(jìn)行Holt-Winters過(guò)濾后的趨勢(shì)變化。

        圖7 Holt-Winters模型圖

        根據(jù)圖7中的Holt-Winters模型發(fā)現(xiàn)數(shù)據(jù)趨勢(shì)與原數(shù)據(jù)基本吻合,再調(diào)用R語(yǔ)言包中的forecast()函數(shù)進(jìn)行數(shù)據(jù)預(yù)測(cè),其中需要設(shè)置的參數(shù)h代表預(yù)測(cè)多少個(gè)數(shù)據(jù)(這里取h=8代表預(yù)測(cè)8個(gè)數(shù)據(jù)),預(yù)測(cè)結(jié)果見(jiàn)表1。

        表1 Holt-Winters預(yù)測(cè)結(jié)果

        表1中,時(shí)間序列3.0對(duì)應(yīng)的預(yù)測(cè)數(shù)據(jù)3.819082便是對(duì)第21個(gè)未來(lái)樣本的預(yù)測(cè)值,以此類(lèi)推。Lo 80與Hi 80表示預(yù)測(cè)的一級(jí)波動(dòng)區(qū)間即80%的置信區(qū)間為(Lo 80, Ho 80),同理Lo95與Hi95表示二級(jí)波動(dòng)區(qū)間,其范圍比一級(jí)波動(dòng)區(qū)間更大,可信度更高但數(shù)據(jù)預(yù)測(cè)更不準(zhǔn)確,波動(dòng)區(qū)間為實(shí)際數(shù)據(jù)波動(dòng)提供參考。最后根據(jù)質(zhì)量控制理論計(jì)算出來(lái)的UCL與LCL管控上下線(xiàn)取對(duì)數(shù)后的新UCL與LCL上下線(xiàn)進(jìn)行質(zhì)量管控,通過(guò)abline()函數(shù)附加在原圖上,再將預(yù)測(cè)值與原數(shù)據(jù)值覆蓋一并繪圖,如圖8所示。

        圖8 預(yù)測(cè)數(shù)據(jù)及波動(dòng)區(qū)間

        根據(jù)圖8發(fā)現(xiàn)當(dāng)前的數(shù)據(jù)(圓點(diǎn)折線(xiàn))并沒(méi)有超出UCL與LCL管控線(xiàn),但是根據(jù)預(yù)測(cè)的結(jié)果(無(wú)圓點(diǎn)折線(xiàn))說(shuō)明:預(yù)測(cè)值越來(lái)越向UCL管控上線(xiàn)偏移,最終將會(huì)超出管控上線(xiàn),未來(lái)可能出現(xiàn)質(zhì)量問(wèn)題。所以此時(shí)應(yīng)該提前預(yù)警,對(duì)可能存在的質(zhì)量問(wèn)題做出檢測(cè)與排查,保證未來(lái)質(zhì)量也能處于可控狀態(tài),避免質(zhì)量問(wèn)題的發(fā)生。

        4 結(jié)論

        本文基于R語(yǔ)言結(jié)合聯(lián)合休哈特累積和控制理論,對(duì)花鍵軸制造過(guò)程進(jìn)行質(zhì)量控制,并高效地建立Holt-Winters模型進(jìn)行質(zhì)量預(yù)測(cè)。對(duì)未來(lái)可能發(fā)生的質(zhì)量問(wèn)題進(jìn)行預(yù)警與排查,將質(zhì)量問(wèn)題解決在萌芽之中。本文質(zhì)量控制與預(yù)測(cè)方法成本低廉,便捷高效,適合向中小型制造企業(yè)進(jìn)行推廣與借鑒。下一步,還可以深入研究R語(yǔ)言與質(zhì)量管理信息系統(tǒng)(QMIS)的集成問(wèn)題,將質(zhì)量控制理論和預(yù)測(cè)方法融入到QMIS的開(kāi)發(fā)中,以提高企業(yè)質(zhì)量管理水平和質(zhì)量信息化目標(biāo),為實(shí)際生產(chǎn)提供有益的管理工具和手段。

        猜你喜歡
        花鍵軸哈特波動(dòng)
        葉片式擺動(dòng)液壓馬達(dá)花鍵軸的45鋼屈服極限分析與強(qiáng)度校核
        羊肉價(jià)回穩(wěn) 后期不會(huì)大幅波動(dòng)
        C35鋼花鍵軸開(kāi)裂原因分析
        微風(fēng)里優(yōu)美地波動(dòng)
        2019年國(guó)內(nèi)外油價(jià)或?qū)⒉▌?dòng)加劇
        花鍵軸扭轉(zhuǎn)變形及破壞分析
        干濕法SO2排放波動(dòng)對(duì)比及分析
        哈特《法律的概念》導(dǎo)讀
        哈特的語(yǔ)義學(xué)
        船王挑選接班人
        亚洲爆乳精品无码一区二区三区| 最近亚洲精品中文字幕| 国内国外日产一区二区| 亚洲一区二区三区四区精品在线| 在线涩涩免费观看国产精品| 九九精品视频在线观看| 中文字幕精品亚洲二区| 99久久精品一区二区国产| 狠狠噜狠狠狠狠丁香五月| 99久久免费精品高清特色大片 | 又爽又黄又无遮挡的激情视频| 99久久综合九九亚洲| 在线观看免费的黄片小视频| 无码av专区丝袜专区| 青青草原精品99久久精品66| 日韩精品一区二区三区四区| 亚洲性av少妇中文字幕| 99麻豆久久久国产精品免费| 少妇邻居内射在线| 国产粉嫩嫩00在线正在播放| 国产精品国产自产拍高清| 亚洲妇女自偷自偷图片| 欧美中文字幕在线| 国产精品成人黄色大片| 国产av无码专区亚洲av麻豆| 天天躁日日躁狠狠躁av| 婷婷激情六月| sm免费人成虐漫画网站| 久久97久久97精品免视看| 色爱区综合激情五月综合小说 | 一区二区三区少妇熟女高潮 | 日韩狼人精品在线观看| 森中文字幕一区二区三区免费| 西西人体444www大胆无码视频| 国产国语对白一区二区三区| 99久久婷婷国产精品综合网站| 中国老熟妇506070| 成人伊人亚洲人综合网站222| 国产亚洲精品视频在线| 搡女人真爽免费视频大全| 亚洲色偷偷综合亚洲av伊人|