周文婷，司玉鵬，何宏舟，王榮杰

(1.集美大學(xué)，福建省能源清潔利用與開(kāi)發(fā)重點(diǎn)實(shí)驗(yàn)室，廈門(mén)，361021；2.集美大學(xué)，福建省清潔燃燒與能源高效利用工程技術(shù)研究中心，廈門(mén)，361021；3.集美大學(xué)，廈門(mén)，361021)

0 序言

回焊爐爐溫曲線反映了焊接過(guò)程中焊點(diǎn)的溫度變化情況[1]，影響爐溫曲線調(diào)節(jié)的因素一般有爐溫曲線峰值溫度、升溫速率、電路板上的溫度分布、預(yù)熱時(shí)間等[2-5].若回焊爐各小溫區(qū)溫度調(diào)節(jié)不良，會(huì)導(dǎo)致電路板變色變形、空洞、潤(rùn)濕不良等，影響產(chǎn)品焊接質(zhì)量和可靠性[6]，因此可以通過(guò)對(duì)爐溫曲線的優(yōu)化設(shè)計(jì)來(lái)提高焊接質(zhì)量.

針對(duì)爐溫曲線優(yōu)化設(shè)計(jì)，李楠等人[7]在RSS(Ramp-Soak-Spike)升溫-保溫-峰值曲線基礎(chǔ)上，利用牛頓冷卻定律建立傳熱模型，運(yùn)用數(shù)值計(jì)算的方法進(jìn)行求解，得到通過(guò)控制傳送帶速度和溫區(qū)溫度會(huì)提高焊接質(zhì)量的結(jié)論，獲得了模擬爐溫曲線，但該方法未考慮溫區(qū)間隙對(duì)回焊爐溫度場(chǎng)的影響，會(huì)導(dǎo)致最終爐溫曲線在相鄰溫區(qū)過(guò)渡時(shí)仿真效果較差.王明泉等人[8]基于傅里葉定律建立一維偏微分熱傳導(dǎo)模型，運(yùn)用有限差分法求解爐溫曲線優(yōu)化參數(shù)，研究如何通過(guò)控制溫區(qū)溫度提高產(chǎn)品質(zhì)量，為實(shí)際生產(chǎn)提供理論支持，但該求解方法的計(jì)算精度還有待提高.方灝航等人[9]利用Fourier 定律和牛頓冷卻定律，建立爐溫曲線差分模型，運(yùn)用模擬退火算法全局尋優(yōu)，獲得高經(jīng)濟(jì)效益基礎(chǔ)上的最大過(guò)爐速度和各溫區(qū)的最優(yōu)穩(wěn)定溫度，然而模擬退火算法可能無(wú)法得到全局最優(yōu)解，影響計(jì)算結(jié)果.Lei 等人[10]建立熱傳導(dǎo)方程模型，使用經(jīng)典粒子群算法計(jì)算熱導(dǎo)率并進(jìn)行誤差分析，利用蒙特卡羅仿真(Monte Carlo simulation)[11]得到傳送帶的最大速度，但該算法易于陷入局部最優(yōu)解，對(duì)仿真結(jié)果產(chǎn)生影響.祝勝光等人[12]提出一種小型智能數(shù)控回焊爐溫度控制器的設(shè)計(jì)和實(shí)現(xiàn)方法，采用PID 算法控制溫度，縮短調(diào)節(jié)時(shí)間，提高溫控精度.綜上可知，目前在研究爐溫曲線模型參數(shù)時(shí)大多采用數(shù)值解法和傳統(tǒng)優(yōu)化算法，計(jì)算繁瑣，計(jì)算結(jié)果的精確性和收斂性還有待提高；在爐溫曲線優(yōu)化設(shè)計(jì)時(shí)，均采用單一目標(biāo)進(jìn)行優(yōu)化，沒(méi)有綜合考慮其它因素影響，整體實(shí)用性較差.因此，設(shè)計(jì)一種有效的方法來(lái)解決多目標(biāo)爐溫曲線優(yōu)化問(wèn)題十分必要.

量子粒子群優(yōu)化算法(quantum particle swarm optimization，QPSO)是Sun 等人[13]基于量子力學(xué)提出的一種新的粒子群優(yōu)化算法.該算法因?qū)崿F(xiàn)簡(jiǎn)單方便、搜索高效準(zhǔn)確、迭代收斂快速、結(jié)果預(yù)測(cè)精準(zhǔn)、需要控制的參數(shù)少，全局搜索能力強(qiáng)，被廣泛用于處理各種工程實(shí)際優(yōu)化問(wèn)題[14-17].文中運(yùn)用QPSO算法，在滿足制程界限的同時(shí)，以加熱因子[18]最小和傳送帶速度最優(yōu)為目標(biāo)，對(duì)回焊爐爐溫曲線進(jìn)行優(yōu)化設(shè)計(jì)，為解決目標(biāo)沖突的回焊爐溫區(qū)溫度設(shè)置提供一種解決方案.

1 爐溫曲線模型建立

1.1 基本參數(shù)

電子元器件的焊接過(guò)程會(huì)經(jīng)歷爐前區(qū)域、爐內(nèi)區(qū)域和爐后區(qū)域.根據(jù)回焊爐內(nèi)部各區(qū)段溫度和功能的不同可分為預(yù)熱區(qū)段、恒溫區(qū)段、回流區(qū)段、冷卻區(qū)段，并分為11 個(gè)小溫區(qū).環(huán)境溫度為25 ℃，回焊爐各溫區(qū)長(zhǎng)度和溫度見(jiàn)表1.

表 1 回焊爐各溫區(qū)長(zhǎng)度和溫度Table 1 Length and temperature of each temperature zone in rewelding furnace

1.2 基本假設(shè)

由于受到模擬和實(shí)際研究難度的限制，將忽略回焊爐內(nèi)溫度場(chǎng)熱輻射的傳熱；假設(shè)電路板的材料均勻，具有各向同性，忽略其內(nèi)部溫差，溫度分布均勻；假設(shè)電路板進(jìn)入溫區(qū)瞬間其內(nèi)部溫度保持穩(wěn)定且不變；由于電路板厚度較小，其沿厚度方向的熱傳導(dǎo)將忽略.

1.3 熱傳導(dǎo)模型

通過(guò)已知各小溫區(qū)溫度可獲得回焊爐內(nèi)溫度隨時(shí)間的數(shù)學(xué)表達(dá)式.各小溫區(qū)溫度恒定不變，則對(duì)應(yīng)回焊爐內(nèi)溫度場(chǎng)的溫度恒定不變；相鄰溫區(qū)若不存在溫差，則溫區(qū)間隙溫度也恒定不變；反之，間隙溫度呈線性變化.

對(duì)于回焊爐空間內(nèi)的物體而言，可建立三維非穩(wěn)態(tài)導(dǎo)熱微分方程，如式(1)所示[19].

式中：τ為 時(shí)間；T為回焊爐內(nèi)的溫度場(chǎng)溫度；ρ為密度；c為比熱容；λ為導(dǎo)熱系數(shù)；Φ˙為單位時(shí)間單位體積中的內(nèi)熱源生成熱；x，y，z分別為該物體在空間x，y和z3 個(gè)方向上的坐標(biāo)值.

由于各間隙的邊界溫度穩(wěn)定不變，可將式(1)簡(jiǎn)化為一維穩(wěn)態(tài)無(wú)內(nèi)熱源常物性導(dǎo)熱微分方程，如式(2)所示.

對(duì)式(2)進(jìn)行微分求解，可得溫區(qū)間隙處的溫度分布，如式(3)所示.

式中：i表示溫區(qū)區(qū)段，i=0，1，2，3···；εi，βi為第i溫區(qū)的常數(shù)參數(shù)；xi為爐前到第i溫區(qū)與(i+1)溫區(qū)之間間隙的距離；Ti表示第i溫區(qū)處對(duì)應(yīng)的溫度場(chǎng)溫度.

回焊爐各小溫區(qū)間隙溫度的線性分布如式(4)所示.

式中：v表示電路板的過(guò)爐速度；表示為從爐前到離開(kāi)第i溫區(qū)區(qū)段的距離；表 示從爐前到進(jìn)入第i+1溫區(qū)區(qū)段開(kāi)始位置的距離.

1.4 電路板焊接中心溫度模型

電路板進(jìn)入回焊爐瞬間，其內(nèi)部溫度保持不變，隨后焊接中心溫度隨時(shí)間開(kāi)始發(fā)生變化，最終趨近于回焊爐的溫度場(chǎng)溫度，故可將其簡(jiǎn)化為非穩(wěn)態(tài)、有內(nèi)熱源的導(dǎo)熱微分方程，如式(5)所示.

式中：t為焊接中心溫度；a為熱擴(kuò)散系數(shù).

由于焊接中心處吸收的熱量與焊接表面上的對(duì)流換熱量相等，得到熱平衡公式，如式(6)所示.

式中：Qa為焊接中心處吸收的熱量.

采用牛頓冷卻定律建立界面換熱模型，如式(7)[20]所示.

式中：h(T)為回焊爐溫度場(chǎng)內(nèi)空氣流體與電路板之間的對(duì)流換熱系數(shù)；A為換熱面積；V為體積.將式(7)化簡(jiǎn)，即

式中：θ=t(h，τ)-Ti為過(guò)余溫度.邊界條件為初始溫度，如式(9)所示.

式中：ti為進(jìn)入第i溫區(qū)時(shí)焊接中心的初始溫度；θi為進(jìn)入第i溫區(qū)時(shí)焊接中心的初始過(guò)余溫度.

綜上所述，可得到焊接中心處溫度變化表達(dá)式，如式(10)所示[21].

式中：T∞為回焊爐溫度場(chǎng)內(nèi)經(jīng)一段時(shí)間傳熱后恒定不變時(shí)溫度場(chǎng)的溫度；取，i=0，1，2···10，RCi的 值僅與回焊爐內(nèi)流體溫度有關(guān)；hi(T)表示第i溫區(qū)回焊爐溫度場(chǎng)內(nèi)空氣流體與電路板之間的對(duì)流換熱系數(shù).

2 多目標(biāo)爐溫曲線優(yōu)化模型

加熱因子是指使回焊爐爐溫曲線溫度從錫膏熔點(diǎn)溫度到峰值溫度所覆蓋的面積，即回焊爐爐溫曲線液相線以上的溫度對(duì)時(shí)間的積分[22-24]，文中采用無(wú)鉛常規(guī)錫膏，熔點(diǎn)為217 ℃.在建立焊接中心溫度變化模型的基礎(chǔ)上，以加熱因子最小和傳送帶速度最大為目標(biāo)，利用量子粒子群算法進(jìn)行多目標(biāo)優(yōu)化，進(jìn)而獲得最優(yōu)回焊爐爐溫曲線.

2.1 目標(biāo)函數(shù)1：加熱因子最小

為使加熱因子最小，建立目標(biāo)函數(shù)如式(11)所示.

式中：ta為無(wú)鉛常規(guī)錫膏熔點(diǎn)，ta=217 ℃；Scover表示加熱因子；τa表示焊接中心溫度第一次達(dá)到ta的時(shí)間；τmax表示焊接中心溫度達(dá)到峰值溫度的時(shí)間；τb表示焊接中心溫度第二次達(dá)到ta的時(shí)間.

2.2 目標(biāo)函數(shù)2：最優(yōu)傳送帶速度

傳送帶速度最優(yōu)，即在滿足爐溫曲線制程界限和工藝要求的前提下過(guò)爐速度最大，如式(12)所示.

式中：Vbest為電路板的最大過(guò)爐速度.

2.3 約束條件

約束條件包括傳送帶的過(guò)爐速度、爐溫曲線的峰值溫度和焊接中心溫度的變化速率、預(yù)熱時(shí)間、焊接區(qū)域溫度超過(guò)無(wú)鉛常規(guī)錫膏熔點(diǎn)所用時(shí)間.

實(shí)際生產(chǎn)時(shí)傳送帶速度需滿足速度限值，如式(13)所示.

電路板在焊接生產(chǎn)過(guò)程中，爐溫曲線需滿足的制程界限如表2 所示.

表 2 傳送帶速度需滿足的制程界限Table 2 Process limits for conveyor speed

3 量子粒子群算法

在量子粒子群優(yōu)化算法中，將粒子的搜索空間從經(jīng)典空間轉(zhuǎn)化為量子空間，并使用薛定諤方程的波函數(shù)對(duì)量子空間中的粒子狀態(tài)進(jìn)行表述.由于局部吸引子，粒子將最終會(huì)聚集到某個(gè)區(qū)域，其如式(14)[25]所示.

式中：o為粒子數(shù)；j為解空間維數(shù)；d為迭代次數(shù)；poj(d)為局部吸引子；Goj(d)為粒子的個(gè)體最佳位置；Zj為種群的全局最佳位置；φoj(d)是區(qū)間(0，1)上具有均勻概率密度函數(shù)的隨機(jī)數(shù).

在解空間中，粒子根據(jù)量子阱和波函數(shù)給出的密度函數(shù)向阱中心收斂，然后使用式(15)和式(16)更新所有粒子.

式中：Xoj(d+1)為粒子o在第(d+1)次迭代中的位置；u為收縮擴(kuò)張因子，以控制量子粒子群的收斂速度，介于[0，1]之間的服從均勻分布的隨機(jī)數(shù)；Pbesdj(d)為粒子群的平均最優(yōu)位置；r為在區(qū)間(0，1)上分布的隨機(jī)數(shù)；M為種群規(guī)模.

為了避免人為因素對(duì)最優(yōu)結(jié)果選擇的影響，利用模糊隸屬函數(shù)[26]從帕累托(pareto)最優(yōu)解集中選擇合適的加熱方案.模糊隸屬函數(shù)計(jì)算的綜合滿意度值越大，說(shuō)明求解的結(jié)果越好.模糊隸屬函數(shù)如式(17)所示.

式中：φk為對(duì)所有目標(biāo)函數(shù)的總體滿意程度，φk值越高說(shuō)明解越好；表示第l個(gè)目標(biāo)的帕累托最優(yōu)解xm的滿意度；分別為第k個(gè)目標(biāo)函數(shù)的最大值和最小值；Nf為帕累托解的個(gè)數(shù)；k為目標(biāo)數(shù)；Nobj為目標(biāo)函數(shù)的數(shù)量[27].

基于量子多目標(biāo)粒子群算法的爐溫曲線優(yōu)化設(shè)計(jì)的算法流程如下.

(1)設(shè)定外部檔案最大儲(chǔ)存量Lmax、種群規(guī)模M等，根據(jù)制程界限確定可行解，在可行解中隨機(jī)初始化粒子群的位置.

(2)根據(jù)公式(11)和公式(12)計(jì)算適應(yīng)度值，將非劣解存入外部檔案中，確定粒子個(gè)體最優(yōu)位置及全局最優(yōu)位置.

(3)根據(jù)公式(15)和公式(16)計(jì)算粒子群的平均最優(yōu)位置，并更新粒子狀態(tài).

(4)判斷更新后的粒子是否滿足制程界限，若滿足，則根據(jù)式(11)和式(12)計(jì)算該粒子的適應(yīng)度值；若不滿足，則舍棄該粒子.

(5)利用pareto 原則更新非劣解集，并將其存入外部檔案.

(6)利用擁擠度排序法完成外部檔案更新.

(7)判斷是否滿足終止條件，若不滿足則返回step3重復(fù)計(jì)算直到達(dá)到最大迭代次數(shù)或滿足終止條件.

(8)利用模糊規(guī)則從pareto 解集中選擇最終方案.

4 模型求解和結(jié)果分析

根據(jù)1.3 節(jié)建立的熱傳導(dǎo)模型通過(guò)仿真模擬獲得回焊爐溫度場(chǎng)的溫度分布曲線如圖1 所示.以回焊爐基本參數(shù)為基礎(chǔ)，運(yùn)用量子粒子群算法和最小二乘法擬合求解焊接中心溫度變化模型中的待定參數(shù)獲得模擬爐溫曲線，并與測(cè)試數(shù)據(jù)進(jìn)行對(duì)比，如圖2 所示.電路板焊接中心溫度變化分布仿真結(jié)果如圖3所示.從圖2～圖3 可以很直觀的看出，模擬爐溫曲線與實(shí)測(cè)數(shù)據(jù)曲線的走勢(shì)大致相同，且收斂較好，驗(yàn)證了所建模型的合理性.

圖 1 回焊爐溫度場(chǎng)溫度分布曲線Fig. 1 Temperature distribution curve of reweld furnace temperature field

圖 2 模擬爐溫曲線與實(shí)測(cè)溫度曲線對(duì)比Fig. 2 Comparison between simulated furnace temperature curve and measured temperature curve

圖 3 電路板焊接中心溫度分布圖Fig. 3 Temperature distribution diagram of circuit board welding center

在上述建立的電路板焊接中心溫度模型的基礎(chǔ)上，為了驗(yàn)證文中提出的多目標(biāo)爐溫曲線優(yōu)化模型的有效性，另外設(shè)置3 組比較情況如下.

(1)以無(wú)鉛常規(guī)錫膏的加熱因子最小為優(yōu)化目標(biāo)，進(jìn)行爐溫曲線優(yōu)化設(shè)計(jì).

(2)以最優(yōu)過(guò)爐速度為優(yōu)化目標(biāo)，進(jìn)行爐溫曲線優(yōu)化設(shè)計(jì).

(3)同時(shí)以無(wú)鉛常規(guī)錫膏的加熱因子最小和最優(yōu)過(guò)爐速度為優(yōu)化目標(biāo)，進(jìn)行爐溫曲線優(yōu)化設(shè)計(jì).

3 種不同優(yōu)化設(shè)計(jì)情況對(duì)比結(jié)果如表3 所示.由表3 可以看出，文中提出的焊接中心溫度變化模型可以同時(shí)考慮最優(yōu)過(guò)爐速度和最小加熱因子.倘若僅考慮最優(yōu)過(guò)爐速度，就會(huì)導(dǎo)致加熱因子增大，焊接區(qū)域中心溫度超過(guò)無(wú)鉛常規(guī)錫膏熔點(diǎn)的時(shí)間過(guò)長(zhǎng)，降低焊接質(zhì)量；若僅考慮最小加熱因子，則會(huì)降低電路板的過(guò)爐速度，影響焊接效率.故僅考慮單一目標(biāo)優(yōu)化的代價(jià)就是影響其它優(yōu)化目標(biāo)，進(jìn)而影響最終優(yōu)化結(jié)果.因此，設(shè)計(jì)一種可以同時(shí)兼顧最優(yōu)速度和最小加熱因子兩個(gè)目標(biāo)的焊接中心溫度變化模型，獲得的最優(yōu)速度為95.55 cm/min，最小加熱因子為1 753.04.

表 3 3 種不同優(yōu)化情況結(jié)果Table 3 Results of three different optimization conditions

以傳送帶過(guò)爐速度最優(yōu)和加熱因子最小為目標(biāo)，采用量子多目標(biāo)粒子群算法進(jìn)行優(yōu)化求解，獲得各溫區(qū)的最優(yōu)爐溫曲線如圖4 所示；圖4 中A，B，C，D，E 點(diǎn)坐標(biāo)如表4 所示；各小溫區(qū)最優(yōu)溫度值如表5 所示.

圖 4 優(yōu)化后的爐溫曲線Fig. 4 Optimized furnace temperature curve

表 4 優(yōu)化后爐溫曲線內(nèi)特殊點(diǎn)坐標(biāo)值Table 4 Coordinate values of special points in the furnace temperature curve after optimization

表 5 優(yōu)化后獲得的最優(yōu)溫區(qū)溫度值Table 5 Optimal temperature value of temperature zone obtained after optimization

相同焊接質(zhì)量條件下，傳送帶速度越快越好，優(yōu)化后的速度為95.55 cm/min，接近速度界限上限.在焊接過(guò)程中，為保證焊接質(zhì)量，焊接區(qū)域中心溫度超過(guò)217 ℃的時(shí)間不宜過(guò)長(zhǎng)，峰值溫度也不宜過(guò)高，故爐溫曲線的溫度峰值越低，越接近制程界限，覆蓋面積越小，加熱因子越小，爐溫曲線越理想，優(yōu)化后的峰值溫度為240.01 ℃，接近溫度峰值界限下限.

由圖4 和表5 可知，焊接中心升溫過(guò)程中位于[150 ℃，190 ℃]溫度區(qū)間的時(shí)間為60.5 s，焊接中心溫度位于[217 ℃，tmax]溫度區(qū)間的時(shí)間為59.5 s，焊接中心溫度變化速率如式(18)所示.

5 結(jié)論

(1)以回焊爐爐溫曲線為研究對(duì)象，在滿足制程界限和實(shí)際工藝要求的前提下，以傳送帶速度最優(yōu)和以無(wú)鉛常規(guī)錫膏的加熱因子最小為目標(biāo)，運(yùn)用量子多目標(biāo)粒子群算法對(duì)爐溫曲線進(jìn)行優(yōu)化設(shè)計(jì)，獲得了最大過(guò)爐速度和各溫區(qū)的最優(yōu)溫度.

(2)運(yùn)用量子粒子群算法和最小二乘法擬合求得焊接中心溫度變化模型中的待定參數(shù)，獲得的模擬爐溫曲線與實(shí)測(cè)數(shù)據(jù)曲線走勢(shì)大致相同且收斂較好.

(3)同時(shí)兼顧最優(yōu)速度和最小加熱因子兩個(gè)目標(biāo)優(yōu)化后的電路板過(guò)爐速度為95.55 cm/min，接近速度界限上限，加熱因子最小為1 753.04；優(yōu)化后的爐溫曲線整體呈現(xiàn)先上升后下降的趨勢(shì)，存在一個(gè)峰值溫度，該峰值溫度為240.01 ℃，接近溫度峰值界限下限.