鄢 威,張 輝,張 華,江志剛,向 琴
(武漢科技大學 機械制造及自動化學院,湖北 武漢 430081)
電弧焊是現(xiàn)代工業(yè)生產(chǎn)中應用最廣泛的一種焊接方法,但其加工過程要耗費大量的資源和能源,同時產(chǎn)生大量的廢氣、煙塵、光輻射等,造成了嚴重的環(huán)境影響。據(jù)統(tǒng)計,我國現(xiàn)階段焊接用鋼量已占我國鋼材消費量的50%以上,經(jīng)焊接的鋼材總數(shù)已超3億噸,居世界第一[1]。另一方面,電弧焊機被列為國家12類高能耗產(chǎn)品之一,其單臺額定功率為10 kW~100kW,且近幾年產(chǎn)銷量均為400萬臺/年以上,居各類焊機之首,總耗電量不容小覷[2]。隨著國家節(jié)能減排戰(zhàn)略的實施及以低碳制造、綠色制造等為代表的一系列關注制造業(yè)可持續(xù)發(fā)展研究的興起,電弧焊過程的節(jié)能、節(jié)材和減排已成為目前研究的熱點問題。
隨著焊接技術的發(fā)展,如自動焊、半自動焊等技術和以激光為代表的高能束流焊接技術的發(fā)展與應用,焊接過程的能耗與環(huán)境排放已得到了較好的控制[3]。但從文獻[2]中電焊機行業(yè)的統(tǒng)計數(shù)據(jù)來看,我國手工電弧焊設備仍占很大比重,近幾年的產(chǎn)銷量均占所有焊接設備的52%左右,且仍以超過300萬臺/年的速度增長。因此,在這些設備的基礎上進行工藝優(yōu)化控制,實現(xiàn)電弧焊加工過程的電能消耗最小化和熱利用率最大化,具有重要的現(xiàn)實意義。
目前,國內外已有不少專家學者對焊接參數(shù)優(yōu)化問題展開了研究。例如:Cristiene等[4]基于逆向黃金分割技術,研究了鎢極惰性氣體保護焊的熱效優(yōu)化方法;Dey等[5]基于遺傳算法,研究了在熔滴滲透最大條件下使焊接區(qū)域最小的焊接參數(shù)優(yōu)化方法;Islam 等[6]以焊接產(chǎn)品質量為目標,利用遺傳算法和有限元分析方法建立了焊接參數(shù)優(yōu)化系統(tǒng);Feng等[7]測量了電流和電壓在一定條件下CO2氣體保護焊中的焊接速度與熱效率的變化規(guī)律,通過數(shù)據(jù)擬合方法建立了熱效最大的焊接速度優(yōu)化模型;Kumar等[8]利用基于灰度的田口方法研究了激光透射焊中以最大接頭強度和最小焊接寬度為目標的焊接參數(shù)優(yōu)化方法;韋智元等[9]利用熱彈塑性有限元法分析了焊接過程耦合殘余應力的分布規(guī)律,提出了通過調節(jié)焊接溫度消除長直焊縫殘余應力的工藝方法;王勤等[10]以最大化焊縫強度為目標,利用灰色關聯(lián)法對電弧焊工藝參數(shù)優(yōu)化問題進行了研究;羅毅等[11]以最小化焊接過程碳排放量為目標,利用廣義回歸神經(jīng)網(wǎng)絡(General Regression Neural Network,GRNN)方法對焊接工藝的參數(shù)優(yōu)化問題進行了研究;Khan等[12]通過實驗方法建立了基于能量密度的不銹鋼材料焊縫預測模型,研究了激光電弧焊工藝加工質量與能耗的多目標優(yōu)化問題;Liu等[13]通過優(yōu)化激光電弧焊的激光運行時間和運行狀態(tài),對焊接過程的成本和能耗的多目標優(yōu)化問題進行了研究。上述文獻大部分研究的都是焊接參數(shù)的單目標優(yōu)化問題,優(yōu)化目標為質量、成本、能耗和熱效率等;部分文獻雖然涉及焊接過程的多目標優(yōu)化,但是優(yōu)化目標均是在質量或成本目標的基礎上加入能耗、碳排放或熱效目標,難以同時實現(xiàn)焊接工藝參數(shù)節(jié)能和高效兩方面的優(yōu)化。
基于以上分析,本文提出一種面向節(jié)能高效的電弧焊工藝參數(shù)多目標優(yōu)化方法。在滿足電弧焊加工生產(chǎn)規(guī)范的基礎上,建立了以焊接電流、焊接速度為變量,以電能消耗和熱有效利用率為優(yōu)化目標,以焊件質量、焊接設備、焊條直徑等對工藝參數(shù)的限制為約束條件的電弧焊加工工藝參數(shù)多目標優(yōu)化數(shù)學模型,并利用自適應進化梯度小生境遺傳算法對該模型進行了求解。最后,通過某截止閥閥座密封面手工電弧焊的實例驗證了模型及方法的可行性和有效性。
電弧焊是一個復雜的物理化學過程,存在大量參數(shù)影響加工過程的能耗和熱效率,如焊接電流、電弧電壓、焊接速度、接頭形式、坡口形狀、焊絲直徑和焊機負載持續(xù)率等。在我國現(xiàn)行標準GB 10854-89和GB/T 983-1995中,焊接坡口接頭形式、形狀尺寸、焊縫形狀和焊絲直徑的選擇均有統(tǒng)一的規(guī)定,因此,這些參數(shù)可依據(jù)相關標準選擇。對焊接電流、電弧電壓和焊接速度等參數(shù),現(xiàn)行標準并沒有統(tǒng)一的要求,但在GB15579.1-2004中,提供了對焊機負載電弧電壓的計算經(jīng)驗公式,為精簡計算量,將其視為焊接電流的函數(shù)。因此,本文選擇焊接電流及焊接速度兩個參數(shù)對電弧焊加工過程的能耗和熱效進行優(yōu)化。
本文的優(yōu)化目標為電弧焊加工過程的節(jié)能與高效,分別對應降低電能消耗和提高熱有效利用率。將焊接質量、焊接設備、焊接方法等焊接參數(shù)的限定作為該優(yōu)化目標的約束條件,以避免出現(xiàn)優(yōu)化結果與焊接質量間的矛盾。
1.2.1 能耗目標函數(shù)
電弧焊加工過程的能耗主要是焊機產(chǎn)生的電能消耗,其大小與焊機的運行狀態(tài)相關。為防止加工過程中焊機過熱,GB/T 8118-2010將焊機工作時間分為若干工作周期T0,其中焊機分別處于負載加工狀態(tài)和空載散熱狀態(tài),不同加工狀態(tài)的持續(xù)時間由負載持續(xù)率決定。因此,通過計算每個工作周期中的焊機能耗,可最終建立電弧焊加工過程的能耗目標函數(shù)。
焊機在一個工作周期內的電能消耗可表示為
式中:E0為焊機在一個工作周期內的電能消耗;P1為焊機在一個工作周期內處于負載加工狀態(tài)的負載功率;P0為焊機在一個工作周期內處于空載散熱狀態(tài)的空載功率;k為焊機負載持續(xù)率;T0為工作周期的時間,國標規(guī)定為10min。
負載持續(xù)率k主要與焊機的負載電流有關:
式中:kN和IN分別為焊機的額定負載持續(xù)率和額定電流,可在焊機銘牌上查得;I2為負載電流。若計算出k≥1,則說明焊機可在此負載電流下連續(xù)工作,取k=1。
焊機負載功率
式中:η2 為焊機的功率因素;U2為電弧電壓,可記為焊接電流的函數(shù),U2=f(I2),具體計算經(jīng)驗公式可依據(jù)GB15579.1-2004查詢。
焊接時間可由焊縫長度和焊接速度計算得出,結合前文對焊機工作周期的分析,焊接時間
式中:H為焊縫長度;V為焊接速度;為焊接時間中包含的焊機工作周期數(shù)。
由式(1)~式(4)可得電弧焊加工過程能耗目標函數(shù):
1.2.2 熱效優(yōu)化目標的確定
電弧焊加工的熱效有兩種定義:①焊接熱效率ηT,指單位時間內焊件吸收的能量(有效熱功率)與電弧總功率的比值;②焊接熔化效率ηm,指單位時間內熔化焊縫金屬所需熱量(包括熔化潛熱)與電弧有效熱功率的比值[14]。由于有效熱功率不僅包含單位時間內熔化焊縫金屬所需的熱量,還包含單位時間內使焊縫金屬處于過熱狀態(tài)的熱焓和向焊縫四周傳導的熱量總和,焊接熱效率并不能反映電弧熔化金屬形成焊縫的實際熱量。有學者[15]提出熱有效利用率ηTe的概念,即通過單位時間內熔化焊縫金屬所需熱量與電弧總功率的比值計算焊接過程的熱效。熱有效利用率能更準確地反映電弧熔化焊縫金屬的實際熱量,因此本文選取熱有效利用率作為電弧焊加工過程的熱效優(yōu)化目標。
電弧焊加工過程熱有效利用率
其中焊接熱效率
式中:Q為有效熱功率,Q0為電弧功率,可根據(jù)焊接電流和焊接電壓計算,
焊接熔化效率
式中:Q1為單位時間內熔化焊縫金屬所需熱量(包括熔化潛熱);Q2為單位時間內使焊縫金屬處于過熱狀態(tài)的熱焓和向焊縫四周傳導熱量的總和;Q=Q1+Q2。
其中
式中:S為焊縫橫截面積;γ為被焊材料比重;Sm為單位重量已熔化金屬的熱焓,
式中:C為材料比熱;Tm為材料熔點;ΔSm為材料熔化的潛熱。
聯(lián)立式(6)~式(8),可得電弧焊加工過程熱效優(yōu)化目標函數(shù):
一般而言,電弧焊工藝參數(shù)的選取受到焊接設備、操作方法、焊材厚度、焊接質量等條件的限制,只能在滿足上述限制條件的范圍內進行取值。本文通過查閱相關焊接手冊、國家標準等,從焊接質量、焊接電流和焊接速度三方面建立電弧焊能耗和熱效優(yōu)化的約束。
(1)焊接質量約束 主要考慮焊接外觀質量,依據(jù)工藝手冊選取焊縫熔深、焊縫寬度和焊縫余高三個主要參數(shù)對焊接質量進行限定,即
式中:h1為焊縫熔深,b為焊縫寬度,h2為焊縫余高,h1min、h1max、bmin、bmax、h2min和h2max分別為焊縫 熔深、焊縫寬度和焊縫余高的最小值和最大值,其取值范圍依據(jù)GB150.4-2011和GB10854-89。
(2)焊接電流約束 焊接電流需要滿足兩個約束條件:①焊機本身的限制,焊機有明確的最小電流輸出Imin和最大電流輸出Imax;②焊接電流的取值必須與所選焊絲直徑相匹配,焊接手冊中規(guī)定了對應焊絲直徑的焊接電流計算系數(shù)Kmin和Kmax。設焊絲直徑為d,焊接電流約束可表示為:
(3)焊接速度約束 在焊接過程中,單位時間內輸入焊縫的熱量太大將導致焊縫處金屬過熱,太小則導致金屬不能充分熔化,均會影響焊接質量。焊接速度與焊接材料允許熱輸入qc、焊接電流和電弧電壓相關,同時也受到焊接設備的約束。焊接速度的約束可表示為:
式中:qc為焊接金屬對輸入熱量的限制,q1max為焊機最大的熱輸入。
為方便計算,將熱效轉換為其倒數(shù),即將求熱有效利用率最大值的優(yōu)化目標轉化為求其倒數(shù)最小值的優(yōu)化目標,則電弧焊節(jié)能高效多目標優(yōu)化模型為:
小生境遺傳算法是目前解決多目標優(yōu)化問題的一種常用方法,它通過設置個體相似度共享函數(shù)來調整個體適應度,達到維護種群多樣性的目的。相比于傳統(tǒng)的遺傳算法(Genetic Algorithm,GA),小生境遺傳算法具有很強的把握搜索空間的能力,更易于找出優(yōu)化問題的所有局部最優(yōu)解和全局最優(yōu)解。但由于小生境共享半徑的取值往往需要預先知道最優(yōu)解的個數(shù),人為設定若取值不當則會影響個體被選擇的概率及物種的多樣性和分布的合理性,導致計算停滯和局部最優(yōu)性能差等缺陷[15-16]。因此,本文提出一種自適應進化梯度的小生境遺傳算法(Adaptive and evolutionary Gradient Niche Genetic Algorithm,AGNGA),將進化梯度和小生境半徑均作為決策變量的一部分參與染色體的編碼,在對問題進行優(yōu)化的同時對個體的小生境半徑進行自適應調整;利用進化梯度信息調整個體向更優(yōu)解進化,并根據(jù)進化代數(shù)自適應調整實數(shù)編碼個體的交叉量和變異量,提高解的精度,加快收斂速度。AGNGA 的流程如圖1所示。
為彌補人為設定小生境半徑的不足,將進化梯度向量與小生境共享半徑均作為決策變量的一部分參與染色體的編碼,最終個體編碼串如下:
其中:xk∈[Mk,Nk]為目標函數(shù)的第k個決策變量,對應于焊接電流和焊接速度;G=[g1,g2,…,gn]為進化梯度向量,σsh為小生境共享半徑。
設目標函數(shù)為f(x),其n維梯度向量為▽f,可用目標函數(shù)的一階偏微分近似計算:
式中:ei為單位向量,δ為一足夠小的正實數(shù)。
梯度向量永遠指向函數(shù)增大的方向,沿該方向f(x)可取得極大值,即Xt+1=Xt+δ·▽f(x),沿該負方向f(x)可取得極小值,即Xt+1=Xt-δ·▽f(x)。
為防止移動量δ·▽f(x)過大造成跳過局部最優(yōu)解、得到一個較小的適應度值,可作以下調整:對種群中的每個個體進行梯度算子運算,計算種群中所有個體的共享程度,由調整后的適應度函數(shù)計算個體的適應度值,若個體適應度值小于原來的適應度值,則說明移動向量過大,需將個體沿著負梯度方向移動兩個單位,若個體的適應度值大于原來的個體適應度值,則說明該個體已經(jīng)進行了梯度的優(yōu)化計算。其具體過程如下:
根據(jù)調整后個體的適應度值大小自適應選擇,選擇算子
結合啟發(fā)式交叉算子和自適應交叉系數(shù),提出改進的自適應非均勻算術交叉算子產(chǎn)生新的個體,計算方法如下:
設變異步長為Δ,第t代的變異點xk的基因取值范圍為,則新基因值
式中:r為[0,1]范圍內的隨機數(shù),T為種群的最大進化代數(shù),b為系統(tǒng)參數(shù)。
綜上所述,小生境遺傳算法的基本步驟如下:
步驟1 算法初始化。建立初始群體,設定遺傳參數(shù)等(小生境半徑的定義域是原決策變量中最大定義域長度的一半。若表示小生境半徑的基因全部初始化為0,則隨機選取一個基因并將其值置為1)。
步驟2 完成個體適應度計算。
步驟3 執(zhí)行選擇、交叉、變異等遺傳操作。
步驟4 進行梯度運算,計算梯度向量。
步驟5 完成種群中所有個體共享度的計算。
步驟6 根據(jù)個體共享度重新計算每個個體的適應度。
步驟7 比較子代和父代個體適應度的大小。
步驟8 用適應度大的子代個體替換父代個體,形成新一代種群。若由于交叉或變異算子改變基因值導致某個個體的峰半徑基因全部為0,則隨機選取一個基因并將其值置為1。
步驟9 若滿足收斂條件,則算法終止;否則返回步驟3。
以某加工車間截止閥閥座密封面手工電弧焊過程為例,對上述模型進行驗證,閥座密封面簡圖如圖2所示。
此閥座密封面加工時選用YD-400AT3直流手工電弧焊機,其參數(shù)如表1所示。
表1 YD-400AT3電弧焊機參數(shù)
電弧焊選用D517焊條,直徑4mm,查詢工藝手冊得此焊條直徑下的焊接電流限定范圍為[120,160]A,可由表1和式(14)確定焊接電流的取值范圍;焊接采用單面坡口,查詢GB150得出熔深、焊縫寬度和余高限定范圍;閥座材質為碳素鑄鋼ZG230-450,可通過機械工程材料手冊等查詢得到計算焊接熱有效利用率時需要用到的材料物理性能參數(shù),如表2所示。
表2 焊件相關參數(shù)
將表2數(shù)據(jù)代入式(12),得到焊接速度的取值范圍。測得焊件厚度為7 mm,焊縫長度為94.2 mm,焊縫橫截面積為18mm2。將以上參數(shù)代入優(yōu)化數(shù)學模型,采用MATLAB 編寫AGNGA 程序,設群體大小為100,最大迭代次數(shù)為200,交叉概率為0.8,變異概率為0.1。焊接過程電能消耗和熱有效利用率的最終計算結果如圖3所示,對應的最優(yōu)解集如表3所示。
表3 最優(yōu)解集對應的參數(shù)
續(xù)表3
3.2.1 AGNGA 優(yōu)化結果分析
由圖3可以看出,電弧焊節(jié)能高效多目標優(yōu)化運算得到的最優(yōu)解集沿Pareto前沿面均勻分布,從而說明了本文選用的AGNGA 的可行性與有效性。同時,表3的最優(yōu)解集表明,隨著焊接電流的不斷增加,與之匹配的焊接速度并非單調遞增,而是會出現(xiàn)一定程度的波動,其最大取值發(fā)生在焊接電流為160A 時,相應的焊接能耗和熱有效利用率也不隨著焊接電流的增加而單調遞增,其中,能耗最大發(fā)生在焊接電流134.3A 時,熱有效利用率最大發(fā)生在焊接電流160A 時。該優(yōu)化解集提供了一個參數(shù)選擇空間,使得工藝人員可依據(jù)實際加工條件和需求選擇合適的電弧焊工藝參數(shù),從而減少了對工藝人員經(jīng)驗的依賴。
3.2.2 優(yōu)化結果與車間實測結果對比分析
在優(yōu)化解集中分別選取能耗最小和熱有效利用率最大時的一組最優(yōu)解進行實際驗證,并將驗證結果與當前車間參數(shù)的實際加工結果進行比較分析,對比結果如表4所示。
表4 優(yōu)化結果與實際結果對比
表4中,焊縫熔深、焊縫寬度和焊縫余高均符合GB150.4-2011和GB10854-89的要求,說明焊接質量達到了要求,表明此多目標優(yōu)化約束條件選取的合理性及解的可行性。相比于實際運行參數(shù),利用優(yōu)化解集中的焊接能耗最小時對應的參數(shù)進行加工,可減少電能消耗約23.66%,但熱有效利用率減少約5.59%;而采用熱有效利用率最大的解對應的參數(shù)加工時,熱有效利用率提高約16.55%,但焊接過程的能量消耗增加約8.45%。對比結果表明:電弧焊節(jié)能高效多目標優(yōu)化模型的最優(yōu)解集可以有效減少手工電弧焊過程的電能消耗,提高熱有效利用率,降低焊接過程的環(huán)境排放。工藝人員可依據(jù)實際需求在最優(yōu)解集中選取所需的參數(shù)進行加工,以達到節(jié)能和高效焊接的目的。
3.2.3 AGNGA 與GA 優(yōu)化結果對比分析
AGNGA 與傳統(tǒng)GA 的優(yōu)化結果對比如圖4所示。對比結果顯示,AGNGA 最優(yōu)解的個數(shù)多于傳統(tǒng)GA,且優(yōu)化結果能更好地貼合Pareto前沿面,從而說明該算法較好地保持了種群多樣性,且具有更好的全局尋優(yōu)能力和收斂速度,可找出本文提出的多目標優(yōu)化問題的局部最優(yōu)解和全局最優(yōu)解。
我國電弧焊設備量大面廣,在所有焊接設備中占很大比重,通過調節(jié)焊接工藝參數(shù)實現(xiàn)節(jié)能高效焊接,有助于減少焊接行業(yè)的總體能量消耗和環(huán)境排放,具有一定的現(xiàn)實意義?;诖?,本文針對電弧焊加工過程電能消耗和熱利用率優(yōu)化問題進行研究,建立了以焊接電流和焊接速度為自變量,以焊件質量、焊接設備、焊條直徑等工藝參數(shù)的限制為約束條件的面向節(jié)能(電能消耗最?。┡c高效(熱有效利用率最大)的電弧焊加工過程多目標優(yōu)化模型。
利用MATLAB編寫自適應進化梯度AGNGA程序,對帶約束的多目標優(yōu)化問題進行了尋優(yōu)求解,并結合某閥座密封面手工電弧焊實例對優(yōu)化解進行了驗證。結果表明,在保證焊接質量的前提下,該方法能有效減少焊接過程的電能消耗、增加熱有效利用率。同時,將本文所用算法與傳統(tǒng)GA 的計算結果進行了比較分析,結果表明,本文提出的AGNGA最優(yōu)解的分布更加合理,具有更好的全局尋優(yōu)能力和收斂速度,更易于找出優(yōu)化問題的局部和全局最優(yōu)解。
本文主要針對電弧焊工藝過程進行優(yōu)化建模,而目前以激光—電弧復合焊接為代表的高能束流焊接技術已逐步開始應用,這將是本領域發(fā)展的重要方向,如何針對新型焊接工藝建立科學合理的節(jié)能高效優(yōu)化模型,將是下一步的研究重點。
[1]LI Wushen,TANG Bogang.The problems for steel,welding features and materials development in China[C]//Proceedings of the 12th National Conference on Welding Academic.Beijing:China Machine Press,2008:1-12(in Chinese).[李午申,唐伯鋼.中國鋼材、焊接性與焊接材料發(fā)展及需要關注的問題[C]//第十二次全國焊接學術會議論文集.北京:機械工業(yè)出版社,2008:1-12.]
[2]Welding Machine Branch of China Electrical Equipment Industry Association.The 2012economic operation analysis report for electric welding machine industry[J].Electric Welding Machine,2013(8):4-11(in Chinese).[中國電器工業(yè)協(xié)會電焊機分會.2012年度電焊機行業(yè)經(jīng)濟運行分析報告[J].電焊機,2013(8):4-11.]
[3]LI Xiaoyan,WU Chuansong,LI Wushen.Study on the progress of welding science and technology in China[J].Journal of Mechanical Engineering,2012,48(6):19-31(in Chinese).[李曉延,武傳松,李午申.中國焊接制造領域學科發(fā)展研究[J].機械工程學報,2012,48(6):19-31.]
[4]GON?ALVES C V,CARVALHO S R,GUIMAR?ESG,et al.Application of optimization techniques and the enthalpy method to solve a 3D-inverse problem during a TIG welding process[J].Applied Thermal Engineering,2010,30(16):2396-2402.
[5]DEY V,PRATIHAR D K,DATTA G L,et al.Optimization of bead geometry in electron beam welding using agenetic algorithm[J].Journal of Materials Processing Technology,2009,209:1151-1157.DOI:10.1016/j.jmatproter.2008.03.019.
[6]ISLAM M,BUIJK A,RAIS-ROHANI M,et al.Simulationbased numerical optimization of arc welding process for reduced distortion in welded structures[J].Finite Elements in Analysis and Design,2014,84:54-64.
[7]FENG Jiecai,LI Liqun,CHEN Yanbin,et al.Effects of welding velocity on the impact behavior of droplets in gas metal arc welding[J].Journal of Materials Processing Technology,2012,212(11):2163-2172.
[8]KUMAR N,RUDRAPATI R,PAL P K.Multi-objective optimization in through laser transmission welding of thermoplastics using grey-based taguchi method[J].Procedia Materials Science,2014(5):2178-2187.
[9]WEI Zhiyuan,LIU Yujun,LI Rui,et al.Residual stresses evolutions of cutting and welding sequential process based on numerical simulation[J].Computer Integrated Manufacturing Systems,2013,19(9):2244-2250(in Chinese).[韋智元,劉玉君,李 瑞,等.基于數(shù)值計算的切割焊接聯(lián)系加工殘余應力演變過程[J].計算機集成制造系統(tǒng),2013,19(9):2244-2250.]
[10]WANG Qin,KUANG Lizhong,ZENG Shenbo.Grey relational analysis of the most optimal parameters of welding process based on arc signals[J].Electric Welding Machine,2010(3):75-78(in Chinese).[王 勤,匡立中,曾申波.基于電弧信號的焊接過程最優(yōu)參數(shù)的灰關聯(lián)分析[J].電焊機,2010(3):75-78.]
[11]LUO Yi,CAO Huajun,LI Hongcheng,et al.Carbon emission model and parameter optimization of CO2shielded welding based on GRNN[J].China Mechanical Engineering,2013,24(17):2398-3402(in Chinese).[羅 毅,曹華軍,李洪丞,等.基于GRNN 網(wǎng)絡的CO2氣體保護焊工藝碳排放建模與參數(shù)優(yōu)化[J].中國機械工程,2013,24(17):2398-3402.]
[12]KHAN M M A,ROMOLI L,DINI G,et al.A simplified energy-based model for laser welding of ferritic stainless steels in overlap configurations[J].CIRP Annals-Manufacturing Technology,2011,60(1):215-218.
[13]LIU Liming,HAO Xinfeng,SONG Gang.A new laser arc hybrid welding technique based on energy conservation[J].Materials Transactions,2006,47(6):1611-1614.
[14]WU Chuansong.Thermal process and molten pool shape in welding[M].Beijing:China Machine Press,2007(in Chinese).[武傳松.焊接熱過程與熔池形態(tài)[M].北京:機械工業(yè)出版社,2007.]
[15]AI Yongyi,KONG Fancong,QIN Jian,et al.The electric energy consumption study of AC and DC arc welding and CO2welding[J].Electric Welding Machine,1983(6):23-26(in Chinese).[艾雍宜,孔凡聰,秦 健,等.交直流手工電弧焊與CO2焊接能耗的研究[J].電焊機,1983(6):23-26.]
[15]WEI Wei,TAN Jianrong,F(xiàn)ENG Yixiong,et al.Multi-objective optimization method research on flexible job shop scheduling problem[J].Computer Integrated Manufacturing Systems,2009,15(8):1592-1598(in Chinese).[魏 巍,譚健榮,馮毅雄,等.柔性工作車間調度問題的多目標優(yōu)化方法研究[J].計算機集成制造系統(tǒng),2009,15(8):1592-1598.]
[16]ZHANG Yan,XIE Juanying.A new niche genetic algorithm based on improved gradient operators[J].Computer Engineering &Science,2009,31(9):71-73(in Chinese).[張 琰,謝娟英.改進梯度算子的小生境遺傳算法[J].計算機工程與科學,2009,31(9):71-73.]