基于工藝場(chǎng)景數(shù)據(jù)的熱處理工藝能耗估算方法

李瑞興，張應(yīng)中

（大連理工大學(xué)機(jī)械工程學(xué)院，遼寧大連 116023）

0 引言

熱處理加工是通過(guò)對(duì)零件進(jìn)行加熱、保溫和冷卻等來(lái)改變零件材料的內(nèi)部組織，以達(dá)到滿足零件力學(xué)性能要求的一種工藝過(guò)程，因此熱處理工藝需要消耗大量的熱能。據(jù)統(tǒng)計(jì)，我國(guó)每年熱處理消耗電能總量達(dá)200 億kW·h[1]。同時(shí)在機(jī)械產(chǎn)品制造過(guò)程中，熱處理是一道十分重要、不可或缺的工藝過(guò)程，熱處理工藝能耗成為機(jī)械產(chǎn)品制造影響環(huán)境負(fù)荷的重要因素之一。

隨著環(huán)境保護(hù)意識(shí)的增強(qiáng)，綠色制造理念和技術(shù)、產(chǎn)品全生命周期評(píng)價(jià)（Life Cycle Assessment, LCA）方法被廣泛采用，LCA方法的基礎(chǔ)工作就是全生命周期清單（Life Cycle Inventory，LCI）分析，LCI分析需要產(chǎn)品全生命周期中的各個(gè)單元過(guò)程的環(huán)境資源負(fù)荷數(shù)據(jù)[2]。熱處理工藝作為產(chǎn)品制造過(guò)程中的一個(gè)重要部分，其能耗數(shù)據(jù)是評(píng)價(jià)工藝經(jīng)濟(jì)性、環(huán)境友好性的重要清單數(shù)據(jù)；然而由于熱處理工藝的離散性、對(duì)熱處理設(shè)備和工藝環(huán)境的依賴性等特性，目前非常缺乏熱處理工藝能耗數(shù)據(jù)。國(guó)內(nèi)外目前主要開(kāi)展關(guān)于機(jī)床切削工藝和設(shè)備能耗預(yù)測(cè)的研究：何彥等[3]從數(shù)控機(jī)床能量源多、加工任務(wù)及加工參數(shù)動(dòng)態(tài)變化等特點(diǎn)出發(fā)，建立了一種數(shù)控機(jī)床多能量源的動(dòng)態(tài)能耗模型與仿真方法；Kabata等[4]以能耗方程和設(shè)備利用系數(shù)為基礎(chǔ)，建立了一套適用于生產(chǎn)或設(shè)計(jì)階段的熱處理工藝能耗計(jì)算方法；林利紅等[5]為解決理論建模復(fù)雜且不具備實(shí)時(shí)性的問(wèn)題，提出了一種基于數(shù)據(jù)驅(qū)動(dòng)的電阻爐多參數(shù)能耗預(yù)測(cè)方法，其參數(shù)主要有溫度、電流、功率等；Carlsson等[6]提出一個(gè)采用統(tǒng)計(jì)建模預(yù)測(cè)電弧爐電能消耗方法。

但是在LCI分析建模中，分析者不一定是熱處理工程師，熱處理工藝也不一定發(fā)生在分析者所在的企業(yè)，很難具備熱處理工藝相關(guān)知識(shí)和數(shù)據(jù)，一般通過(guò)查詢清單數(shù)據(jù)庫(kù)獲取工藝過(guò)程能耗數(shù)據(jù)。熱處理工藝繁多，不同的工藝類型、工件材料和熱處理設(shè)備等工藝場(chǎng)景下能耗相差很大，完全采集每一種工藝和工藝場(chǎng)景數(shù)據(jù)下的能耗數(shù)據(jù)是十分困難的，也是不現(xiàn)實(shí)的，需要一個(gè)基于工藝場(chǎng)景數(shù)據(jù)的參數(shù)化模型，對(duì)熱處理能耗進(jìn)行估算。本文面向LCI分析建模需求，提出基于工藝場(chǎng)景數(shù)據(jù)知識(shí)的熱處理工藝能耗估算方法，通過(guò)輸入熱處理基本工藝場(chǎng)景數(shù)據(jù)，系統(tǒng)在知識(shí)驅(qū)動(dòng)下自動(dòng)檢索相關(guān)數(shù)據(jù)，估算出其能耗數(shù)據(jù)，可以為熱處理工藝LCI清單分析提供數(shù)據(jù)支撐。

1 熱處理工藝場(chǎng)景數(shù)據(jù)概念模型

工藝場(chǎng)景是指工藝在其生命周期中所發(fā)生活動(dòng)的特征描述，是工藝過(guò)程輸入信息的集合。熱處理工藝場(chǎng)景信息主要包括熱處理工藝類型、工件材料、熱處理爐設(shè)備、加熱方式、工藝方案、功能單位等。其中：熱處理工藝類型和工件材料決定了熱處理工藝的溫度時(shí)間曲線，對(duì)熱處理工藝的能耗起到了決定性的作用；熱處理爐設(shè)備的加熱性能和保溫性能及工藝方案對(duì)熱處理工藝的能源消耗也有重要影響。

從上述分析可以看出，熱處理工藝場(chǎng)景信息離散且相互關(guān)聯(lián)，包含豐富的語(yǔ)義信息和熱處理工藝知識(shí)；此外，熱處理工藝場(chǎng)景動(dòng)態(tài)變化，涉及的材料、熱處理設(shè)備、工藝方案眾多，不同企業(yè)的工藝場(chǎng)景都不同，要完全收集所有場(chǎng)景信息是非常困難的，需要構(gòu)建一個(gè)開(kāi)放和可共享的熱處理工藝場(chǎng)景信息模型，能夠完整一致和用戶可自定義地表示熱處理工藝場(chǎng)景信息和知識(shí)，在該信息模型基礎(chǔ)上，用戶根據(jù)自己特定的工藝場(chǎng)景添加和修改數(shù)據(jù)和知識(shí)。

本文采用本體技術(shù)提出一個(gè)面向熱處理能耗計(jì)算的工藝場(chǎng)景信息概念模型，如圖1所示。

圖1 熱處理工藝場(chǎng)景信息概念本體模型

熱處理工藝場(chǎng)景信息主要概念類如下：1）熱處理工藝類。如前面所說(shuō)的，要獲得不同力學(xué)性能，需要采用不同的熱處理工藝，因此在熱處理工藝類下派生出多個(gè)工藝子類，能耗影響較大的熱處理工藝主要是淬火、正火、退火和回火等。退火又可以分為完全退火、不完全退火、等溫退火、球化退火和去應(yīng)力退火等，從而形成一個(gè)熱處理工藝層次結(jié)構(gòu)。2）熱處理工件類。熱處理工件結(jié)構(gòu)形狀對(duì)工件加熱和散熱有一定影響，相似結(jié)構(gòu)基本上具有相似的散熱特性。本文將工件類子類定義毛坯類、軸套類、叉架類和箱體類等，同時(shí)每個(gè)工件由金屬材料組成，與工件金屬材料類通過(guò)“工件材料”關(guān)系連接。3）熱處理工件材料類。熱處理工件材料直接決定熱處理加熱溫度、保溫時(shí)間和冷卻方式，對(duì)能耗有直接影響。對(duì)材料進(jìn)行分類是必要的，按照碳鋼、合金鋼和鑄鐵等三大類主要金屬材料建立層次化的熱處理工件材料類。每一類材料聯(lián)系金屬材料數(shù)據(jù)庫(kù)表，表中存儲(chǔ)較詳細(xì)的金屬材料信息。4）熱處理爐設(shè)備類。熱處理爐設(shè)備決定了熱處理加熱方式，對(duì)熱處理能耗有直接影響。熱處理爐設(shè)備通常分為電阻爐設(shè)備和天然氣設(shè)備。電阻爐設(shè)備又分為臺(tái)車爐、井式爐等。5）加熱工藝方案類。熱處理加熱工藝方案是指加熱工件的方案，例如滿負(fù)荷（裝滿加熱爐）、半負(fù)荷、分區(qū)加熱、工裝等。

2 熱處理工藝能耗估算理論模型

2.1 熱處理工藝能耗影響因素

一個(gè)工藝過(guò)程通常是一個(gè)能量輸入與輸出的平衡過(guò)程。要構(gòu)建熱處理工藝能耗估算模型首先從能量輸出構(gòu)成和特性、熱處理工藝中的能量流動(dòng)因素進(jìn)行分析。

熱處理的工作設(shè)備一般比較簡(jiǎn)單，其能量主要來(lái)自電阻爐的電能消耗或者燃?xì)鉅t的燃?xì)庀模疚闹饕芯侩娮锠t的能量消耗。一般情況下，熱處理爐會(huì)根據(jù)設(shè)定的熱處理溫度時(shí)間工藝曲線對(duì)熱處理爐溫度進(jìn)行調(diào)控，電熱元件發(fā)出熱量，借輻射與對(duì)流作用將熱量傳給加熱的工件，從而使工件加熱到規(guī)定的溫度。

熱處理工藝能量輸出主要包括加熱工件的消耗的熱量（有效熱量）、工藝過(guò)程中的各項(xiàng)損失熱量?jī)纱蟛糠郑饕蛩厝鐖D2所示。主要影響因素[7]包括加熱工件能耗Q工件、加熱輔助構(gòu)件（支撐架、爐底板、工夾具及料盤等）能耗Q輔助件、熱處理爐內(nèi)氣體能耗Q氣體、通過(guò)爐壁的散熱能耗Q散熱、砌體蓄熱能耗Q蓄熱及其它熱損失能耗Q其它。

圖2 熱處理能耗影響因素示意圖

2.2 熱處理工藝能耗估算理論模型

根據(jù)上述熱處理工藝影響因素分析和熱能量輸入和輸出平衡原則，可以推導(dǎo)出如下熱處理工藝能耗估算理論模型[7]：

1）周期作業(yè)的熱處理工藝。周期作業(yè)的熱處理工藝中，在加熱階段的能量輸出是加熱工件、輔助構(gòu)件所需能量和砌體蓄熱的能量。熱處理過(guò)程中實(shí)際的蓄熱能量與熱處理冷卻階段、裝卸階段和停爐期造成的爐體降溫程度有關(guān)。其能量輸入和輸出平衡方程為

上述熱平衡方程中各項(xiàng)耗能計(jì)算如下：

a.加熱工件能耗。加熱工件能耗是熱處理過(guò)程中工件從初始溫度達(dá)到一定溫度所吸收能量。此部分能耗為工件熱處理過(guò)程中需要達(dá)到內(nèi)外溫度一致、顯微組織轉(zhuǎn)變完全所需要消耗的能量，是熱處理過(guò)程中唯一的有效能耗。可以通過(guò)工件初始和終了溫度、工件質(zhì)量及其比熱容計(jì)算得到。可以采用積分相對(duì)精確計(jì)算，其計(jì)算公式為

式中：m為工件總質(zhì)量，kg；t0和t1分別為工件加熱的初始和終了溫度，℃；c（t）為工件比熱容和溫度的函數(shù)[8]，kJ/（kg·℃）。

b.加熱輔助構(gòu)件能耗。在退火熱處理過(guò)程中，退火爐中的料框、工夾具、支撐架、爐底板及料盤等輔助構(gòu)件和工件一樣會(huì)吸收能量、升高溫度。加熱輔助構(gòu)件能耗通過(guò)其質(zhì)量、比熱容和溫度通過(guò)式（1）計(jì)算。

c.加熱控制氣體能耗。為了使工件表面不發(fā)生氧化脫碳現(xiàn)象或?qū)ぜM(jìn)行化學(xué)熱處理，需要向爐內(nèi)通入可進(jìn)行控制成分的氣氛。控制氣體在熱處理爐中升溫會(huì)吸收能量，加熱控制氣體的耗能由輸入控制氣體的體積、控制氣體比熱容和保溫溫度決定，其計(jì)算公式為

式中：V為氣體用量；t0和t1分別為氣體入爐前溫度和工作溫度；C為氣體在t0和t1溫度范圍內(nèi)的平均比熱容。加熱爐內(nèi)氣體能耗主要在控制氛圍爐中考慮，由于此部分能耗過(guò)小，在非控制氛圍爐和真空爐中一般忽略不計(jì)。

d.通過(guò)爐壁的散熱損失能耗。根據(jù)熱力學(xué)原理散失到外界的能量通過(guò)三層傳熱構(gòu)成：爐氣對(duì)爐墻內(nèi)表面的綜合傳熱（對(duì)流傳熱和輻射傳熱）；爐墻內(nèi)部的傳導(dǎo)傳熱；爐墻外表面對(duì)空氣的綜合傳熱。爐壁散熱可以用熱流密度表示，即單位面積和單位時(shí)間通過(guò)的熱能。已知爐壁兩側(cè)溫度各為t1、t2，爐壁的厚度為s，導(dǎo)熱系數(shù)為λ，α為爐墻外表面對(duì)空氣的綜合給熱系數(shù)，則對(duì)于n層爐墻的傳熱過(guò)程，可以得到爐壁散熱能耗為

e.砌體蓄熱能耗。砌體蓄熱能耗為爐子從室溫加熱至工作溫度并且達(dá)到穩(wěn)定狀態(tài)時(shí)爐襯本身吸收的熱量。采用隔熱和保溫雙層爐襯結(jié)構(gòu)的砌體蓄熱能耗損失計(jì)算公式[7]為

f.其它能耗。此項(xiàng)耗能包括未考慮到的各種熱損失及一些不易精確計(jì)算的熱損失，如爐襯密封不嚴(yán)、爐子長(zhǎng)期使用后保溫材料隔熱性能降低和爐子密封性降低等造成的熱損失、加熱元件的熱短路、電氣控制和工件傳送能耗等。此項(xiàng)耗能具有較大的不確定性，本文取上述各項(xiàng)能耗總和的10%～20%。

3 基于工藝場(chǎng)景數(shù)據(jù)的熱處理能耗估算

3.1 熱處理工藝場(chǎng)景數(shù)據(jù)與能耗估算關(guān)聯(lián)模型

實(shí)際上熱處理工藝場(chǎng)景數(shù)據(jù)基本上決定了熱處理工藝的能耗，要實(shí)施參數(shù)化的工藝能耗估算，構(gòu)建熱處理工藝場(chǎng)景數(shù)據(jù)與能耗估算的關(guān)聯(lián)模型是必要的。基于上述分析，本文提出一個(gè)圖3所示的熱處理工藝場(chǎng)景數(shù)據(jù)與能耗估算的關(guān)聯(lián)模型。

圖3 熱處理工藝場(chǎng)景數(shù)據(jù)與能耗估算的關(guān)聯(lián)模型

從圖3可以看出，熱處理工藝類型和熱處理工件材料決定了熱處理溫度時(shí)間曲線，工件材料決定了不同溫度下比熱容，再結(jié)合工件質(zhì)量就確定出工件加熱的能耗；散熱能耗、蓄熱能耗和其它能耗主要由加熱爐設(shè)備確定；輔助能耗和氣體能耗主要由熱處理工藝類型、工藝方案和加熱爐設(shè)備確定。下面主要闡述工件加熱消耗的熱能估算。

3.2 工件加熱能耗計(jì)算

根據(jù)上述熱處理能耗估算的理論模型和關(guān)聯(lián)模型，工件加熱所需的能耗完全可以從熱處理類型、熱處理工件信息計(jì)算出來(lái)。主要計(jì)算步驟如下。

3.2.1 熱處理工藝溫度時(shí)間曲線

熱處理溫度時(shí)間曲線實(shí)際就是熱處理工藝曲線，標(biāo)明了熱處理加熱溫度、保溫時(shí)間、冷卻速度變化狀況[9]。圖4給出一個(gè)熱處理工藝的溫度時(shí)間曲線示意圖。

從圖4可以看出，不同熱處理工藝加溫到不同的溫度。例如完全退火工藝將亞共析鋼加熱至Ac3（轉(zhuǎn)變?yōu)閵W氏體的終了溫度）以上20～30 ℃，保溫足夠時(shí)間，奧氏體化后，隨爐緩慢冷卻，從而接近平衡的組織，而不完全退火將亞共析鋼在Ac1～Ac3（Ac1為珠光體向奧氏體轉(zhuǎn)變的溫度）之間或過(guò)共析鋼在Ac1～Accm（Accm為溶入奧氏體的終了溫度）之間兩相區(qū)加熱。因此獲取材料的Ac1、Ac3、Accm等臨界溫度是很重要的。

研究表明，這些臨界溫度與材料中的碳元素、合金元素和雜質(zhì)元素等相關(guān)[9]。在處理工件鋼/鑄鐵材料中，總會(huì)含有一定數(shù)量的碳元素、合金元素和雜質(zhì)元素等，這些元素所占的比例決定了鋼/鑄鐵的力學(xué)特性。熱處理過(guò)程中，需要知道熱處理工件材料晶體變化的臨界溫度，例如Ac1線，即超過(guò)Ac1（共析）線，材料奧氏體會(huì)分解為珠光體。

通過(guò)實(shí)驗(yàn)研究得到了金屬材料的臨界溫度的經(jīng)驗(yàn)計(jì)算公式[10]，根據(jù)金屬材料中碳元素、合金元素及雜質(zhì)元素等的質(zhì)量百分比計(jì)算材料的Ac1、Ac3、Accm等臨界溫度：

通過(guò)臨界溫度計(jì)算公式可以方便地獲取材料的臨界溫度，是熱處理加熱溫度確定的基礎(chǔ)。但是通過(guò)經(jīng)驗(yàn)公式得到的臨界溫度可能存在誤差，基于實(shí)際測(cè)量的材料臨界溫度構(gòu)建材料臨界溫度數(shù)據(jù)庫(kù)，輸入材料名稱即可查詢獲取材料臨界溫度。

3.2.2 熱處理工件材料比熱容曲線

圖5 碳質(zhì)量分?jǐn)?shù)為0.23%的碳素鋼在不同溫度下的平均比熱容

3.2.3 基于三次樣條插值的數(shù)值積分計(jì)算工件加熱能耗

按照2.2節(jié)給出的加熱工件能耗計(jì)算公式，工件能耗計(jì)算是一個(gè)從室溫T0到加熱終了溫度T1之間，工件質(zhì)量乘以材料的比熱容溫度函數(shù)的積分。為了方便計(jì)算機(jī)求解，本文采用基于三次樣條插值的數(shù)值積分計(jì)算工件加熱能耗。

對(duì)于給定的列表函數(shù)c(t):{(ti，yi)}，（T0=t1＜t2＜，…，＜tn=T1）（i=1，2，…，n），采用三彎矩法構(gòu)造比熱容三次樣條插值函數(shù)。設(shè)hj=tj+1-tj，（j=1，2，…，n-1），則有如下計(jì)算式[11]：

上述公式將對(duì)比熱容積分計(jì)算離散轉(zhuǎn)化為三次樣條插值的求和運(yùn)算，方便了計(jì)算機(jī)編程的實(shí)現(xiàn)。

4 估算計(jì)算編程實(shí)現(xiàn)與采集數(shù)據(jù)驗(yàn)證

4.1 估算計(jì)算編程實(shí)現(xiàn)

本文基于上述熱處理工藝場(chǎng)景數(shù)據(jù)概念模型及熱處理工藝能耗計(jì)算的理論模型，采用Qt編程實(shí)現(xiàn)一個(gè)基于工藝場(chǎng)景參數(shù)的熱處理工藝能耗估算軟件系統(tǒng)，該系統(tǒng)主要包括如下3個(gè)模塊。

1）人機(jī)交互界面模塊。通過(guò)人機(jī)界面將工藝場(chǎng)景數(shù)據(jù)概念本體以列表框形式引導(dǎo)用戶選擇和輸入工藝場(chǎng)景數(shù)據(jù)、輸出估算計(jì)算的結(jié)果和編輯工藝場(chǎng)景數(shù)據(jù)概念本體。圖6中給出一個(gè)熱處理工藝場(chǎng)景數(shù)據(jù)選擇界面示意截圖。

圖6 工藝場(chǎng)景數(shù)據(jù)選擇界面示意圖

2）熱處理工藝場(chǎng)景數(shù)據(jù)概念本體庫(kù)和數(shù)據(jù)庫(kù)模塊。熱處理工藝場(chǎng)景數(shù)據(jù)概念本體庫(kù)以XML格式存儲(chǔ)的數(shù)據(jù)概念模型，用戶可以根據(jù)自己的需要進(jìn)行編輯修改。熱處理工藝場(chǎng)景數(shù)據(jù)庫(kù)采用微軟公司的Access 2010構(gòu)建，建立了常用熱處理金屬材料比熱容數(shù)據(jù)庫(kù)、熱處理爐設(shè)備參數(shù)數(shù)據(jù)庫(kù)等。

3）熱處理工藝能耗計(jì)算模塊。根據(jù)選擇和輸入的工藝場(chǎng)景數(shù)據(jù)，完成熱處理能耗的計(jì)算。

4.2 采集數(shù)據(jù)驗(yàn)證

為了驗(yàn)證上述參數(shù)化熱處理工藝能耗估算的有效性，本文分別在3個(gè)不同企業(yè)使用便攜式鉗式功率計(jì)采集了6組電爐加熱熱處理工藝的能耗數(shù)據(jù)，如表1所示。

表1 6組熱處理工藝場(chǎng)景下采集的能耗數(shù)據(jù)

將表2的工藝場(chǎng)景數(shù)據(jù)分別輸入到本文開(kāi)發(fā)的熱處理工藝能耗估算軟件系統(tǒng)中，得到如表2所示的能耗估算結(jié)果。

表2 6組熱處理工藝場(chǎng)景下估算的能耗數(shù)據(jù) kW·h

從表3可以看出，估算的熱處理工藝總能耗與實(shí)際采集得到的能耗誤差為3.3%～13.4%，在可以接受的誤差范圍內(nèi)，驗(yàn)證了本文方法的可行性。同時(shí)，從能耗分布來(lái)看，在熱處理工藝過(guò)程中能耗主要用于工件加熱，加熱爐能耗、輔助構(gòu)建能耗及其它不確定的能耗也是重要的影響因素。

5 結(jié)論

熱處理工藝消耗大量的能源資源，是產(chǎn)品LCI清單分析的重要一環(huán)，需要一個(gè)在工藝場(chǎng)景參數(shù)驅(qū)動(dòng)下的熱處理工藝的能耗數(shù)據(jù)。本文通過(guò)分析熱處理工藝場(chǎng)景數(shù)據(jù)信息，基于能量平衡方程，構(gòu)建熱處理工藝能耗估算理論模型；構(gòu)建熱處理工藝場(chǎng)景數(shù)據(jù)概念模型和能耗估算的關(guān)聯(lián)模型；開(kāi)發(fā)了能耗估算軟件，基于熱處理工藝場(chǎng)景數(shù)據(jù)進(jìn)行熱處理能耗估算；從不同企業(yè)采集了6組熱處理工藝場(chǎng)景數(shù)據(jù)和能耗數(shù)據(jù)，與本文提出的估算方法計(jì)算結(jié)果相比，誤差在合理范圍內(nèi)，表明該方法可以較為方便和準(zhǔn)確地預(yù)測(cè)熱處理工藝能耗。

由于熱處理工藝和設(shè)備的復(fù)雜性，不確定性因素較多，下一步將采集更多的工藝場(chǎng)景數(shù)據(jù)，從其不確定的能耗數(shù)據(jù)分析中找出規(guī)律，從而更準(zhǔn)確地估算熱處理工藝能耗。