蔣麒麟 薛小強(qiáng) 李翔英
(①南京航空航天大學(xué)機(jī)電學(xué)院,江蘇 南京210016;②南京工程學(xué)院機(jī)械學(xué)院,江蘇南京211167)
鈑金工藝中最重要的3 個(gè)步驟是切割、沖裁和折彎,數(shù)控切割是鈑金件制造過(guò)程中一個(gè)極其重要的環(huán)節(jié)。在大批量制造環(huán)境下,對(duì)鈑金件快速、合理、高效地?cái)?shù)控切割是鈑金件后續(xù)加工工藝的重要前提,同時(shí)也是鈑金件產(chǎn)品制造過(guò)程中工藝設(shè)計(jì)、工時(shí)與成本估算的一個(gè)重要基礎(chǔ)。
鈑金件數(shù)控切割工藝一般由3 部分組成,分別是預(yù)熱打孔、切割和空行程,其利用數(shù)控切割機(jī),根據(jù)不同的零件、具體的加工工藝編寫數(shù)控代碼,對(duì)已在整塊鈑金件上排樣好的零件輪廓進(jìn)行自動(dòng)切割。鈑金件的數(shù)控切割主要使用4 種切割方法[1]:(1)氣體火焰切割;(2)等離子弧切割;(3)激光切割(如圖1);(4)高壓水射流切割(如圖2)。
工時(shí)估算是生產(chǎn)經(jīng)營(yíng)管理的主要基礎(chǔ)工作之一,它是掌握生產(chǎn)進(jìn)度情況、安排生產(chǎn)計(jì)劃、進(jìn)行成本核算的基礎(chǔ)。在生產(chǎn)調(diào)度中,要利用零件的工序工時(shí)信息確定供貨期和原材料需求計(jì)劃;在車間管理中,要利用零件的工時(shí)信息合理地分配工人和加工設(shè)備的任務(wù);在生產(chǎn)控制中,通過(guò)利用工時(shí)和材料定額來(lái)控制各加工設(shè)備資源和物料的消耗;工時(shí)估算也是產(chǎn)品成本計(jì)算與報(bào)價(jià)和工人工作業(yè)績(jī)考核的主要依據(jù)[2-3]。
基于實(shí)例的推理(case - based reasoning,簡(jiǎn)稱CBR)是人工智能領(lǐng)域中較新崛起的一種重要的基于知識(shí)的問(wèn)題求解和學(xué)習(xí)方法,自上世紀(jì)80 年代出現(xiàn)以來(lái),已被廣泛應(yīng)用于設(shè)計(jì)、工藝過(guò)程規(guī)劃、故障診斷、分析和決策支持等領(lǐng)域。在大批量定制領(lǐng)域,CBR 主要被應(yīng)用于產(chǎn)品配置和工時(shí)與成本的快速估算之中。通常,CBR 工作被描述為一個(gè)推理循環(huán),這個(gè)循環(huán)包括6個(gè)部分,分別是[4]:
(1)表示(render):對(duì)目標(biāo)問(wèn)題進(jìn)行特征描述。
(2)檢索(retrieve):從實(shí)例庫(kù)中檢索與目標(biāo)問(wèn)題最相似的實(shí)例。
(3)重用(reuse):將最相似實(shí)例的解決方案應(yīng)用于目標(biāo)問(wèn)題。
(4)修正(revise):為了適應(yīng)新問(wèn)題,需要時(shí)對(duì)目標(biāo)實(shí)例的實(shí)際數(shù)據(jù)進(jìn)行修改。
(5)評(píng)價(jià)(remark):對(duì)修正結(jié)果進(jìn)行一定的測(cè)試、評(píng)定。
(6)保留(retain):方案確定生效后,將新實(shí)例保存在實(shí)例庫(kù)中。
CBR 系統(tǒng)中的關(guān)鍵步驟是如何正確有效地表示實(shí)例以及實(shí)例特性的描述、實(shí)例庫(kù)的組織和檢索、實(shí)例的匹配與重用、實(shí)例的更新和維護(hù)。
在一個(gè)CBR 系統(tǒng)中,實(shí)例庫(kù)將隨著應(yīng)用的積累而不斷擴(kuò)大,因此,如何在一個(gè)大型數(shù)據(jù)庫(kù)中高效地檢索目標(biāo)實(shí)例是成功應(yīng)用CBR 的關(guān)鍵之一。實(shí)例的檢索策略主要有以下幾種:最近鄰策略、知識(shí)歸納策略及知識(shí)引導(dǎo)策略。目前應(yīng)用最廣泛的就是最近鄰策略(the nearest neighbor method,NNM),通常用公式表示為[5]:
式中:T表示目標(biāo)實(shí)例,S表示實(shí)例庫(kù)中的源實(shí)例,n表示每一實(shí)例所包含的特征個(gè)數(shù),i∈[1,n];f是目標(biāo)實(shí)例T和源實(shí)例S中特征的相似性函數(shù);ω 表示的是特征i的重要性權(quán)值,即權(quán)重系數(shù),。計(jì)算得到的相似度(similarity)越高,則表示匹配程度越好。
對(duì)于大批量生產(chǎn)的鈑金件,由于加工工時(shí)數(shù)據(jù)具有重用性,因此,可以采用相似實(shí)例零件的典型工時(shí)來(lái)推算產(chǎn)品族中新零件的工時(shí)。零件工時(shí)的重用性主要是由其相似性所決定的。從機(jī)械加工工時(shí)計(jì)算的角度出發(fā),零件的相似性主要包括以下內(nèi)容[6]:
(1)零件結(jié)構(gòu)的相似性。在一個(gè)零件族中,不同零件的加工工藝都是基本相同的。盡管其幾何尺寸、精度要求、重量以及細(xì)節(jié)結(jié)構(gòu)等不盡相同,但是其總體結(jié)構(gòu)具有很強(qiáng)的相似性,零件結(jié)構(gòu)的相似性帶來(lái)工藝的相似性。
(2)加工過(guò)程的相似性。對(duì)于不同切割加工設(shè)備、不同材料、不同厚度的鈑金件,其一般工步都是上料-裝夾-空行程-預(yù)熱-打孔-切割-空行程-卸夾-下料,這種周而復(fù)始、循環(huán)的、極其相似的加工過(guò)程是本文使用CBR 方法的重要原因。
(3)制造資源的相似性。零件的制造資源主要是指加工設(shè)備以及技術(shù)工人。對(duì)于多數(shù)制造企業(yè)而言,加工設(shè)備通常更新周期比較長(zhǎng),上崗技術(shù)工人的技術(shù)水平也比較穩(wěn)定,因此,制造資源也具有相似性。
(4)制造環(huán)境的相似性。制造環(huán)境主要是指影響加工工時(shí)的企業(yè)環(huán)境,如庫(kù)存地點(diǎn)、廠房布局、設(shè)備維修和后勤保障等。對(duì)于同一個(gè)企業(yè),這些情況往往變化不大,因此,零件制造環(huán)境同樣具有相似性。
根據(jù)鈑金件數(shù)控切割的特點(diǎn),結(jié)合各種主要切割方法,可確定影響切割工時(shí)的7 個(gè)主要因素,分別是:附加工時(shí)、切割輪廓總周長(zhǎng)、噴嘴空行程總長(zhǎng)度、勻速切割速度、空行程速度、單個(gè)孔預(yù)熱打孔時(shí)間和打孔點(diǎn)個(gè)數(shù)。對(duì)每一個(gè)評(píng)價(jià)指標(biāo),根據(jù)其具體的物理意義,對(duì)應(yīng)著一系列的評(píng)價(jià)值,如:對(duì)于“勻速切割速度”這一評(píng)價(jià)指標(biāo),不同的零件厚度,對(duì)應(yīng)著不同的評(píng)價(jià)值。由此,切割工時(shí)實(shí)例可表示為:
式中:T為切割工時(shí)實(shí)例;f為附加工時(shí)(包括準(zhǔn)備時(shí)間、上下料時(shí)間等);c為切割輪廓總周長(zhǎng);s為噴嘴空行程總長(zhǎng)度;υ 為勻速切割速度;υ'為空行程速度;p為單個(gè)孔預(yù)熱打孔時(shí)間;為n為打孔點(diǎn)個(gè)數(shù)。
在定量表達(dá)零件加工工時(shí)相似程度的過(guò)程中,本文使用模糊層次分析法建立模糊判斷矩陣。模糊層次分析法(fuzzy analytic hierarchy process,F(xiàn)-AHP 法)[7]是一種定量與定性相結(jié)合的多目標(biāo)決策分析方法。這一方法的核心是將決策者的經(jīng)驗(yàn)判斷給予量化,從而為決策者提供定量形式的決策依據(jù)。其主要思想為:對(duì)于一個(gè)復(fù)雜的多規(guī)則評(píng)價(jià)問(wèn)題,可將評(píng)價(jià)指標(biāo)和評(píng)價(jià)對(duì)象劃分成多個(gè)層次,對(duì)同一層次上的元素,通過(guò)指標(biāo)之間的兩兩比較,對(duì)系統(tǒng)中各指標(biāo)的重要程度給予評(píng)判,組成模糊判斷矩陣,并利用這種評(píng)判結(jié)果來(lái)綜合計(jì)算各指標(biāo)的權(quán)重系數(shù)。
模糊判斷矩陣是一個(gè)n×n的對(duì)稱矩陣,矩陣的行、列數(shù)等于評(píng)價(jià)指標(biāo)的個(gè)數(shù)。根據(jù)心理學(xué)家提出的“人區(qū)分信息等級(jí)的極限能力為7 ±2”的研究結(jié)論,F(xiàn) -AHP法在對(duì)指標(biāo)的相對(duì)重要程度進(jìn)行測(cè)量時(shí),引入了九分位的相對(duì)重要的比例標(biāo)度,標(biāo)度及其值如表1 所示。
表1 標(biāo)度及其值
在征求工程師與專家意見的基礎(chǔ)上,通過(guò)對(duì)上述影響切割工時(shí)的7 個(gè)主要評(píng)價(jià)指標(biāo)的分析、對(duì)比,確定模糊判斷矩陣A 如下:
設(shè)為列向量的第i個(gè)元素,并將判斷矩陣A第i行的幾何平均值作為:
則得向量
將向量作歸一化處理,即令
得權(quán)重向量
顯然,有
獲得工時(shí)估算主要評(píng)價(jià)指標(biāo)的權(quán)重系數(shù)之后,根據(jù)相關(guān)系數(shù)法,即夾角余弦算法,可得到兩零件加工工時(shí)的相似度。相關(guān)系數(shù)法公式為[8]:
式中:r為估算零件與實(shí)例零件的相似度值;m為評(píng)價(jià)指標(biāo)數(shù)量(本文中m=7);Wi為每個(gè)評(píng)價(jià)指標(biāo)的權(quán)重;Xc為實(shí)例零件評(píng)價(jià)指標(biāo)值;Xd為估算零件評(píng)價(jià)指標(biāo)值;設(shè)定閾值R∈(0,1),只有在估算零件與實(shí)例零件相似度值r>R時(shí),實(shí)例零件的結(jié)果值才被調(diào)用。
指數(shù)平滑法是一種在逐次實(shí)際觀察基礎(chǔ)上,全面地對(duì)某一預(yù)測(cè)模型的估計(jì)系統(tǒng)進(jìn)行修正的預(yù)測(cè)方法。該方法的特點(diǎn)是以預(yù)測(cè)對(duì)象的前期實(shí)際值和前期的預(yù)測(cè)值為基礎(chǔ),經(jīng)過(guò)修勻處理后得到本期的預(yù)測(cè)值,指數(shù)平滑法的通用遞推公式如下[9]:
式中:F為預(yù)測(cè)值;Y為觀察值(實(shí)際值);α 為平滑系數(shù),其取值范圍為0≤α≤1;t為時(shí)間序列。
從式(7)可以看出,它是一種以時(shí)間定權(quán)的加權(quán)平均預(yù)測(cè),每一遞推觀察值的加權(quán)系數(shù)[α,α(1 -α),α(1 -α)2,…]按指數(shù)規(guī)律遞減。工時(shí)估算過(guò)程中,以產(chǎn)品的相似度取代指數(shù)平滑法中的時(shí)間維作為加工工時(shí)相似排序的新方法,從而使指數(shù)平滑法能夠科學(xué)地應(yīng)用于相似產(chǎn)品的工時(shí)估算。平滑系數(shù)α 的確定可采用經(jīng)驗(yàn)準(zhǔn)則法、最小變異系數(shù)法、0.618 優(yōu)選法和自適應(yīng)選取法等多種方法[9]。由此,相似件的切割加工工時(shí)估算值為:
設(shè)估算零件切割工時(shí)估算評(píng)價(jià)指標(biāo)值分別為f=1,c=1450,s=350,υ =60,υ' =90,p=3,n=6,某實(shí)例零件相應(yīng)的值分別為f=0.5,c=1890,s=70,υ=80,υ' =100,p=3,n=15,根據(jù)式(5)中的結(jié)果,代入式(6),可得估算零件與某實(shí)例零件的相似度值為0.733(計(jì)算過(guò)程略。)
設(shè)定匹配閾值R為0.7,檢索實(shí)例庫(kù)中相似實(shí)例。若無(wú)滿足條件實(shí)例,有兩條途徑,一是依照式(9)計(jì)算切割工時(shí),二是降低閾值,重新檢索。對(duì)于第一條途徑,可通過(guò)ObjectARX 編程按步驟計(jì)算即得結(jié)果;對(duì)于第二條途徑,依照指數(shù)平滑法修改檢索到的符合條件的實(shí)例。利用指數(shù)平滑法修改實(shí)例時(shí),結(jié)合匹配閾值來(lái)選取實(shí)例(一般R≥0.6)。平滑系數(shù)的確定采用經(jīng)驗(yàn)準(zhǔn)則法,當(dāng)時(shí)間序列波動(dòng)不大時(shí),α 值一般在0.05~0.3 之間選取,時(shí)間序列波動(dòng)較大時(shí),α 值一般在0.3~0.7 之間選取[10]。根據(jù)工時(shí)估算的特點(diǎn),經(jīng)過(guò)試算,選取α=0.5。
若大于匹配閾值的實(shí)例個(gè)數(shù)為3,相似度分別為0.749,0.823,0.775,工時(shí)值分別為1.6 min,14.3 min,18.2 min,則計(jì)算得到的估算零件的加工工時(shí)為:
在估算新零件加工工時(shí)的過(guò)程中,若實(shí)例庫(kù)較小,或閾值較高,難以檢索到合適的匹配實(shí)例,則可直接提取切割工程圖中的幾何元素,按下式計(jì)算切割工時(shí):
式中:a為準(zhǔn)備工時(shí);t0為上下料與裝夾時(shí)間;c為切割輪廓總周長(zhǎng);s為噴嘴空行程總長(zhǎng)度;υ 為勻速切割速度;υ'為空行程速度;p為單個(gè)孔預(yù)熱打孔時(shí)間;n為打孔點(diǎn)個(gè)數(shù)。
對(duì)于固定的切割設(shè)備和一定材質(zhì)、厚度的鈑金件,其a、t0可通過(guò)定額工時(shí)很方便地獲得,υ、υ'、p用戶可根據(jù)具體情況確定,c、s、n可從工程圖中直接提取。圖3為基于實(shí)例推理的鈑金件數(shù)控切割工時(shí)估算流程。
通過(guò)分析鈑金件數(shù)控切割加工工藝過(guò)程,建立了鈑金件數(shù)控切割工時(shí)計(jì)算公式,將鈑金件數(shù)控切割工時(shí)表示為實(shí)例,運(yùn)用CBR 方法,檢索并修正相似實(shí)例,在多品種、大批量生產(chǎn)環(huán)境下,推理、估算鈑金件數(shù)控切割加工工時(shí)。
工時(shí)估算與定額是企業(yè)生產(chǎn)管理工作中獲得基礎(chǔ)數(shù)據(jù)的方法之一,是企業(yè)編制主生產(chǎn)計(jì)劃、合理組織制造活動(dòng)的前提,是改進(jìn)生產(chǎn)技術(shù)、挖掘勞動(dòng)潛力、提高勞動(dòng)效率的重要手段,同時(shí)也是企業(yè)內(nèi)部進(jìn)行成本估算、評(píng)價(jià)經(jīng)濟(jì)效益的重要依據(jù)。工時(shí)定額文件制定的質(zhì)量直接影響著生產(chǎn)效率、設(shè)備利用率、產(chǎn)品生產(chǎn)周期、合同交貨期、企業(yè)員工勞動(dòng)報(bào)酬等方面的內(nèi)容,對(duì)企業(yè)的發(fā)展至關(guān)重要。本文在完成鈑金件數(shù)控切割工時(shí)估算之后,考慮其他若干因素,可進(jìn)行工時(shí)定額文件的制定,這對(duì)于產(chǎn)品成本估算與快速報(bào)價(jià)具有積極的推動(dòng)作用。
[1]李亞江. 切割技術(shù)及應(yīng)用[M]. 北京:化學(xué)工業(yè)出版社,2004.
[2]朱歷新. 基于神經(jīng)網(wǎng)絡(luò)的工時(shí)定額技術(shù)研究[D]. 西安:西北工業(yè)大學(xué),2003.
[3]董豐收. 基于工藝標(biāo)準(zhǔn)化的工時(shí)定額制定技術(shù)研究[D]. 西安:西北工業(yè)大學(xué),2006.
[4]王凌志. CBR 理論、方法及其在沖模CAD 中應(yīng)用的研究[D]. 大連:大連理工大學(xué),2003.
[5]孔繁勝. 知識(shí)庫(kù)系統(tǒng)原理[M]. 杭州:浙江大學(xué)出版社,2001.
[6]楊青海,祁國(guó)寧,黃哲人,等. 基于案例推理和事物特性表的零件工時(shí)估算方法[J]. 機(jī)械工程學(xué)報(bào),2007(5):99 -105.
[7]Billaudel P,Devillez A,Villermain Lecolier G. Performance evaluation of fuzzy classification methods designed for real time application[J]. International Journal of Approximate Reasoning,1999(20):1-20.
[8]Rissland E L,Skala D B. Combining case-based and rule-based reasoning:A heuristic approach[C]. IJCAI -89 Proceeding of the 11th International Joint Conference on Artificial Intelligence,1959,1:524-530.
[9]吳立之. 工程投資估算的新方法——模糊指數(shù)平滑法[D]. 大連:大連理工大學(xué),2003.
[10]王云霞,高成沖. 實(shí)例推理技術(shù)在定制產(chǎn)品相似件成本估算中的應(yīng)用[J]. 機(jī)械設(shè)計(jì),2007(4):65 -67.