李洪丞 馬國喜 羅 蓉 羅志勇

重慶郵電大學(xué)先進制造工程學(xué)院，重慶，400065

0 引言

研究表明制造業(yè)加工過程能源效率通常小于30%[1-2]，GUTOWSKI等[3]指出生產(chǎn)線能源效率只有14.8%。目前，中國已成為全球最大的鋁型材生產(chǎn)國，2015年我國鋁加工材料達到3億噸，鋁合金擠壓材料產(chǎn)量達到1.4億噸[4-5],鋁擠壓企業(yè)超過820家，其中4.5～250 MN的擠壓機超過5 800臺，占世界總臺數(shù)的85%以上[6]。熱擠壓是鋁及鋁合金加工的核心工藝之一，具有高效、無切屑、材料利用率高等優(yōu)點，但目前能耗強度偏高。我國主要鋁加工廠擠壓材生產(chǎn)的能耗是歐洲的1.43倍，溫室效應(yīng)是歐洲的2.33倍[7]。

隨著擠壓成形技術(shù)在交通運輸、航空航天等領(lǐng)域應(yīng)用越來越廣泛，擠壓成形工藝能耗問題受到越來越多關(guān)注。INGARAO等[8-9]基于生命周期清單，對比分析了熱擠壓成形與傳統(tǒng)車削工藝加工同一種軸對稱鋁合金零件的能源效率，并分析了零件批量、尺寸變化對單產(chǎn)能耗的影響。BUIS等[10]建立了擠壓成形工藝生命周期清單模型，分析了不同工藝參量與加工條件對能耗性能的影響。JEONG等[11]針對分流擠壓的鋁棒加熱效率問題，提出鋁棒勻質(zhì)熱處理與鋁棒預(yù)熱過程融合法，并評估了該方法能耗性能情況。王安哲等[12]對銅管熱擠壓與冷擠壓工藝能耗、碳排放進行了對比評估，并分析了關(guān)鍵工藝參數(shù)對擠壓工藝能耗及碳排放的影響。石建業(yè)[13]分析了電阻加熱爐與工頻感應(yīng)加熱爐在鋁及鋁合金熱擠壓中的加熱效果和能源消耗情況，并指出工頻感應(yīng)加熱具有較好的節(jié)能效果。然而，擠壓成形過程屬于復(fù)雜的熱力耦合過程，對擠壓過程熱平衡要求嚴格，而且擠壓成形過程受擠壓速度、溫度、鋁棒長度、擠壓模具幾何尺寸、材料性能等多種參數(shù)影響，現(xiàn)有研究更強調(diào)對擠壓設(shè)備或擠壓工藝能耗的宏觀統(tǒng)計分析或局部改進技術(shù)的研究，對擠壓成形能耗優(yōu)化、控制方法研究仍較匱乏。

本文以降低擠壓成形比能耗為目標，建立熱平衡約束條件等價方程，以鋁棒起始溫度、鋁棒長度、擠壓速度、擠壓筒溫度、模角等為控制變量，以熱導(dǎo)率、材料強度系數(shù)為擾動變量，集成利用實驗設(shè)計(DOE)、方差分析等方法構(gòu)建熱平衡約束的擠壓成形比能耗魯棒優(yōu)化模型，在保證擠壓穩(wěn)定的情況下，降低擠壓成形能耗。

1 擠壓成形過程能耗與熱流平衡特征

熱擠壓是高度非線性的熱力耦合過程，能源消耗用于滿足金屬塑性變形條件，以實現(xiàn)物質(zhì)轉(zhuǎn)移，達到工件形狀、尺寸變化以及改善組織、性能的要求，金屬成形過程如圖1所示。

圖1 擠壓成形過程能耗與熱流平衡Fig.1 Energy balance and heat flow balance of heat extrusion forming

從能耗輸入視角來看，鋁合金熱擠壓前先對鋁棒、模具、擠壓筒進行加熱，電能經(jīng)電動機、液壓系統(tǒng)提供擠壓力，擠壓力會用于材料變形做功和克服摩擦做功。根據(jù)能量守恒定律，擠壓成形能耗可以表達為

(1)

式中，P為擠壓力；vin為擠壓速度；Wcf為金屬流動過程中克服與擠壓筒摩擦力所做的功;Wdf為金屬流動過程中克服模具摩擦力做的功；Wform為金屬變形做的功。

在材料塑性變形過程中，從力平衡的角度看，變形區(qū)入口球面的應(yīng)力是與包括擠壓模錐面或死區(qū)界面上的摩擦力在內(nèi)的其他應(yīng)力相平衡的，因此金屬變形功Wform包括了模錐面或死區(qū)界面上的摩擦功消耗。為了計算該部分的功，本文將根據(jù)文獻[14]提出的擠壓變形球面坐標系(圖2)，基于圓棒擠壓功計算模型，采用系數(shù)修正法獲取變形功。

圖2 擠壓變形時球面坐標系及應(yīng)力應(yīng)變狀態(tài)[14]Fig.2 Spherical coordinate system and stress-strain state during extrusion deformation [14]

為實現(xiàn)塑性變形而消耗在變形區(qū)內(nèi)的功率，可由變形區(qū)入口球面上的應(yīng)力與速度分布通過積分求得：

(2)

εe=lnλ

式中，vj為擠壓變形入口處的速度,本文假設(shè)該速度等于擠壓速度vin;At為擠壓筒的橫截面積;Dt為擠壓筒內(nèi)部直徑；σk為與擠壓條件有關(guān)的金屬變形抗力;α為模角;λ為擠壓比；mz為變形區(qū)模錐面或死區(qū)界面上的摩擦因子(取值范圍為0≤mz≤1.0，無潤滑熱擠壓或產(chǎn)生死區(qū)時，mz=1.0；帶潤滑熱擠壓無死區(qū)時，mz=0.5；冷擠壓無死區(qū)時，mz=0.2～0.3[14])。

材料與擠壓筒接觸面的摩擦功Wcf可根據(jù)摩擦力做功進行計算：

隨著高校招生規(guī)模的擴大，大學(xué)畢業(yè)生就業(yè)問題備受關(guān)注，畢業(yè)生遞增，市場趨于飽和，就“用工荒，就業(yè)難”的尷尬局面：一方面部分企業(yè)找不到自身需要的復(fù)合型人才，另一方面畢業(yè)生又認為找不到自己理想的工作，很大一部分原因是學(xué)生職業(yè)素養(yǎng)的缺失，在很多方面不能滿足企業(yè)需求，基于此，提高大學(xué)生的職業(yè)素養(yǎng)成為高校育人的重要內(nèi)容。

Wcf=Fcfvj

(3)

其中摩擦力Fcf可由下式計算：

Fcf=πτDtL

(4)

式中，L為鋁錠與擠壓筒的接觸長度；τ為鋁棒與擠壓筒接觸面上的摩擦剪切應(yīng)力。

綜合式(3)與式(4)，材料與擠壓筒接觸面的摩擦功Wcf可以表示為

(5)

Wcf可以看成是接觸長度L的線性函數(shù)。在材料擠壓過程中，材料會填充到整個模具中，此時材料與模具的接觸面是固定的，材料與模具接觸面的摩擦功Wdf可以視為常數(shù)，通過下式計算：

(6)

式中,vf為產(chǎn)品出口處速度；d為產(chǎn)品直徑；Ld為定徑帶長度;md為擠壓材料與模具接觸面的摩擦因子。

ΔQ=Qcf+Qdf+Qform-Qd1-Qd2-Qd3

(7)

式中，ΔQ為凈生熱量，等溫擠壓時，理想值為0；Qcf為坯料與擠壓筒之間的摩擦生熱量；Qdf為坯料與模具之間的摩擦生熱量；Qform為金屬變形區(qū)的變形生熱量；Qd1為坯料與擠壓墊片之間的傳熱量；Qd2為坯料與擠壓筒之間的傳熱量；Qd3為坯料與模具之間的傳熱量。

式(7)中，Qd1、Qd2、Qd3均可根據(jù)材料熱傳導(dǎo)方程進行計算確定。

2 熱平衡約束等價方程

根據(jù)擠壓過程中熱流平衡機理，假設(shè)變形區(qū)內(nèi)的金屬發(fā)生塑性變形產(chǎn)生的熱全部作用在?？壮隹谔幍臄D壓產(chǎn)品上，并引起?？壮隹谔帞D壓產(chǎn)品溫度升高，即可假設(shè)?？壮隹谔幃a(chǎn)品的溫度T是坯料進入變形區(qū)時的溫度T0與變形區(qū)內(nèi)塑性變形熱導(dǎo)致的溫升ΔT之和，即

T=T0+ΔT

(8)

其中，ΔT可以看作是熱擠壓過程中塑性變形生熱導(dǎo)致的溫升，已有研究通過機理分析、經(jīng)驗等解析出了溫升預(yù)測公式[14]：

(9)

式中,Af為產(chǎn)品斷面面積，cm2；l為鋁棒長度，cm；vout為產(chǎn)品流出速度,cm/s；Q為擠壓模附近塑性變形區(qū)內(nèi)產(chǎn)生的熱量，J；K為經(jīng)由擠壓筒和擠壓模的熱流系數(shù)，J/(s·℃)；T0為初始坯料溫度加上剪切變形所產(chǎn)生的溫升，再減去擠壓模初始溫度所得的值，℃；c為金屬的質(zhì)量熱容，J/(g·℃)；ρ為金屬的密度，g/cm3。

任何一種鋁合金材料都有一個最佳擠出溫度，實際溫度與最佳擠出溫度之間的波動量越小，擠出產(chǎn)品質(zhì)量越穩(wěn)定。假設(shè)某種產(chǎn)品最佳擠出溫度為Tbest，則產(chǎn)品擠出溫度的波動量

δT=T-Tbest

(10)

綜合以上分析，熱平衡約束條件可以定義為

δT→min(T-Tbest)

(11)

3 成形比能耗魯棒優(yōu)化模型

3.1 比能耗優(yōu)化流程

本文采用圖3所示的成形比能耗優(yōu)化模型來優(yōu)化成形能耗，該模型是建立在擠壓成形能耗平衡與熱平衡機理分析基礎(chǔ)上的，從該模型中可知，若要實現(xiàn)能效的魯棒性優(yōu)化，需進行以下幾個方面研究。

圖3 成形能耗的優(yōu)化過程Fig.3 Optimization process of forming energy

(1)成形過程能耗平衡、熱平衡機理分析。該部分主要是根據(jù)金屬擠壓變形過程，通過對擠壓力受力狀態(tài)、擠壓力影響因素等進行分析，一方面確定擠壓過程能耗平衡機制，另一方面確定擠壓過程熱生成、傳遞、轉(zhuǎn)換機理，確定熱平衡機制。

(2)熱平衡控制原則分析。根據(jù)擠壓過程中金屬的傳熱行為特征，建立擠壓過程的熱流平衡關(guān)系，明確擠壓控制變量對熱平衡及出?？桩a(chǎn)品溫度的影響規(guī)律，在此基礎(chǔ)上確定擠壓成形過程中熱平衡控制原則，即確定工藝控制變量的選擇基準。

(3)變量類型及范圍設(shè)定。根據(jù)擠壓過程特點，將關(guān)鍵工藝參數(shù)、設(shè)計參數(shù)、擾動參數(shù)定義為擠壓成形過程的控制變量集，并確定不同變量的取值范圍。本文選取模角、鋁棒長度、鋁棒起始溫度為設(shè)計變量(DV)，擠壓速度、擠壓溫度為工藝變量(PV)，而描述擠壓工藝設(shè)備、產(chǎn)品設(shè)計參數(shù)的變量定義為固定變量(FV), 對影響成形能耗但難以控制的變量集定義為擾動變量(NV)。

(4)能效優(yōu)化模型。將擠壓成形比能耗定義為優(yōu)化目標，構(gòu)建能效的優(yōu)化模型，并基于DOE構(gòu)建魯棒性優(yōu)化算法，實現(xiàn)能效的優(yōu)化，其中比能耗可采用下式計算：

(12)

F(DV,PV,NV)為單位時間內(nèi)產(chǎn)品擠出量與能耗的比值，定義為產(chǎn)品擠出速度與產(chǎn)品成形功耗的比值，單位為mm/J，該值越大表明擠壓成形能效越高。

(5)能效魯棒優(yōu)化算法的構(gòu)建。主要基于DOE構(gòu)建能效的魯棒性算法，在擠壓熱平衡約束下實現(xiàn)能效的魯棒優(yōu)化。

3.2 能效魯棒性優(yōu)化算法

根據(jù)圖1對成形能耗的分析，構(gòu)建了能效魯棒性優(yōu)化算法，如圖4所示。其中，步驟(3)為比能耗、熱約束評價指標的計算模擬?；跀D壓成形能耗分析模型以及控制變量、擾動變量的正交矩陣確定不同輸入模式下(工藝參數(shù)、設(shè)計參數(shù)不同組合方案確定了不同輸入模式)成形比能耗、溫度波動量的評價指標的均值、方差以及變異系數(shù)等。其中變異系數(shù)可采用下式計算：

η=-10lg{(1/μ2)·[1+3(σ2/μ2)]}

(13)

式中，μ為成形比能耗或溫升的均值；σ為成形比能耗或溫升的方差。

圖4 魯棒優(yōu)化算法過程Fig.4 Robust optimization algorithm

4 案例研究

4.1 實驗設(shè)計

實驗選取圖5所示的鋁合金零部件為研究對象，該零部件材料為6063鋁合金，其材料性能參數(shù)如表1所示。為了加工該鋁合金件，采用圖6所示的擠壓模具。上模尺寸為φ100 mm×50 mm，上模設(shè)有4個分流孔；下模尺寸為φ100 mm×50 mm，焊合高度為20 mm，設(shè)有一級空刀槽，空刀槽起到支撐工作帶和避免模具與型材之間的接觸減少摩擦的作用。鋁合金零部件的擠壓條件如表2所示。

圖5 零件尺寸Fig.5 Parameters of aluminum profile

表1 6063鋁合金零件材料性能參數(shù)Tab.1 Material parameters of 6063 aluminum alloy

圖6 擠壓模具結(jié)構(gòu)Fig.6 Parameters of extrusion dies

4.2 實驗結(jié)果分析

針對本文所給出的實驗對象的加工要求可確定優(yōu)化模型中的固定變量、控制變量、擾動變量類型及范圍分別如表3～表5所示。

表2 鋁合金零件擠壓條件Tab.2 Aluminum alloy parts extrusion condition

表3 固定變量Tab.3 Fixed variables

表4 控制變量類型及范圍Tab.4 Type and range of control variables

表5 擾動變量Tab.5 Noise variables

擾動變量包括金屬-擠壓筒摩擦因數(shù)、金屬熱導(dǎo)率、金屬極限強度系數(shù)、金屬-模具的傳熱系數(shù)等，該類變量都非常難控制，如金屬-擠壓筒摩擦因數(shù)會隨潤滑條件等隨機波動，金屬熱導(dǎo)率與金屬極限強度系數(shù)會隨每批材料質(zhì)量的差異波動等。本實驗將各個變量的經(jīng)驗值作為平均值，并假設(shè)所有變量都滿足正態(tài)分布，通過正交矩陣模擬的方法確定其方差、標準差等，然后基于方差、標準差等進一步確定擾動變量的取值范圍，結(jié)果如表6、表7所示。

表6 擾動變量的值、差值及標準差Tab.6 Values, difference, and standard deviations for NV

表7 正交試驗中擾動變量的水平值Tab.7 Level of noise variables in orthogonal test

根據(jù)表3～表7給出的各變量的取值范圍，利用Minitab17軟件可確定控制變量與擾動變量的正交矩陣，計算結(jié)果如表8、表9所示。根據(jù)上述數(shù)據(jù)，采用MATLAB平臺對魯棒優(yōu)化過程進行模擬計算，可以確定表10所示的成形比能耗的平均值、方差以及變異數(shù)。

表8 控制變量的L16(45)正交矩陣Tab.8 L16 (45) orthogonal array for control variables

表9 擾動變量的L16(43)正交矩陣Tab.9 L16(43) orthogonal array for noise variables

根據(jù)表10可以確定最佳比能耗方案，當(dāng)T0=420 ℃,l=560 mm,vin=330 mm/s,Tc=375 ℃,α=60°時，最佳比成形能耗為95.88 mm/J，溫度波動量為1.78 ℃。如圖7所示，本文利用效應(yīng)曲線圖進一步揭示了控制變量對比能耗的影響。當(dāng)T0增大時，比能耗存在一定的波動性，且存在一個最佳值，鋁棒溫度過低會導(dǎo)致金屬擠壓變形應(yīng)力增大，比能耗降低；為保證擠壓產(chǎn)品出口

表10 成形比能耗的均值、方差及變異系數(shù)Tab.10 Mean, variance, and variation coefficient for specific energy

圖7 控制變量的效應(yīng)曲線圖分析Fig.7 Effect curve of control variables

處的溫度，鋁棒溫度過高時，應(yīng)減小擠壓速度vin；此時比能耗也減少；l的增大會導(dǎo)致比能耗降低，主要是因為鋁棒長度增大，鋁棒與擠壓筒之間的摩擦阻力增大，能耗增大；Tc的增大會導(dǎo)致擠壓比能耗降低，Tc增大意味著加熱能耗的增加，但對比能耗的影響并不是很大；隨著模角α的增大，擠壓比能耗升高，主要原因是模角增大，金屬與模具的摩擦力減小，能耗需求降低，同時也增大了產(chǎn)品出口處的速度。

5 結(jié)論

本文揭示了擠壓成形過程熱平衡與能耗平衡機制，建立了熱平衡的等價約束方程。以優(yōu)化擠壓成形比能耗為目標，提出了熱平衡約束的擠壓成形比能耗魯棒優(yōu)化模型，并基于實驗設(shè)計、方差分析等方法提出比能耗魯棒算法。以某6063鋁合金零件的熱擠壓為例，確定當(dāng)初始溫度為420 ℃、鋁棒長度為560 mm、擠壓速度為330 mm/s、擠壓筒溫度為375 ℃、模角為60°時，最佳比成形能耗為95.88 mm/J，溫度波動量為1.78 ℃。而且根據(jù)各控制變量的效應(yīng)曲線，鋁棒初始溫度與擠壓速度對擠壓比能耗的影響最大。但由于生產(chǎn)實際中影響因素更復(fù)雜，需要進一步結(jié)合所提方法與工程試驗，確定最佳參量的可信度范圍，以便更好支撐生產(chǎn)實際。