關 聰，陳 英，王 哲

（廣西民族大學，廣西 南寧 530000）

1 問題背景

回焊爐是集成電路板等電子產品生產的重要工具之一。在生產過程中電子制造行業(yè)中印刷電路板（PCB）回流焊工藝是SMT技術重要的一環(huán)，回流焊溫度曲線會實時記錄焊接過程溫度的變化。在現(xiàn)實的生產過程中，不少企業(yè)認為確?；亓骱笭t處于良好的制程狀態(tài)是質量控制的關鍵，因此需要對回焊爐溫度曲線進行持續(xù)的測量檢查，所以爐溫曲線的研究可以保證焊接質量，降低測試成本，確?；亓骱笭t具有優(yōu)良的制程能力。回焊爐截面示意圖如圖1所示。

圖1 回焊爐截面示意圖

2 目標問題

本文旨在對工具機理模型進行分析研究，為有效測量各溫區(qū)溫度值，并確定最優(yōu)傳送帶傳送速度，建立溫度變化規(guī)律數(shù)學模型，從而解決爐溫曲線相應問題。

問題一：在給定各溫區(qū)溫度設定值條件下，計算焊接區(qū)域中心溫度變化情況，并畫出相應爐溫曲線。制程界限如表1所示。

表1 制程界限

問題二：以峰值溫度為中心線的兩側超過217℃的爐溫曲線應盡量對稱，結合給出的傳送帶最大過爐速度，畫出最優(yōu)曲線、各溫區(qū)設定的溫度及傳送帶過爐溫度及相應指標。

3 模型的假設

模型的假設如下：①集成電路板的邊界長l=2.5 cm，且兩端邊界相等；②兩塊相等溫度的集成電路板邊界與間隙的最低溫度差值ΔT=170℃；③過爐速度勻速穩(wěn)定，不存在突變；④受傳熱影響，相同溫度的兩塊小溫區(qū)之間間隙中心溫度低于5℃；⑤假設受傳熱影響，第一個邊界略高于室溫為26℃；⑥不考慮電子元件材料等對溫度測量和焊接的影響；⑦假設更改小溫區(qū)溫度值后，加熱溫度最終趨于穩(wěn)態(tài)；⑧假設過爐過程中，電路板表面加熱溫度即為小溫區(qū)溫度；⑨假設調整小溫區(qū)非邊界溫度都能精確到設定值；⑩假設在調節(jié)溫區(qū)溫度時，原先速度始終為該狀態(tài)下最大過爐速度。

4 模型的建立與求解

4.1 模型的建立

4.1.1 穩(wěn)態(tài)熱傳導過程分析

4.1.1.1 導熱微分方程

滿足直角坐標系一維、穩(wěn)態(tài)、常物性、無內熱源、第一邊界條件下控制方程為：

將控制方程積分計算得出通解：

4.1.1.2 回焊爐內熱傳導穩(wěn)態(tài)溫度分布分析

回歸本題中，假設將小溫區(qū)看成一有限長的細狀熱源，回焊爐內有以下情況的穩(wěn)態(tài)熱傳導過程：爐前區(qū)域與第一小溫區(qū)之間、溫度相等的溫區(qū)之間、溫度由低到高的兩各溫區(qū)之間、溫度由高到低的溫區(qū)之間。

假設剛進入爐前初始邊界溫度為T00，在每一段有限長熱源的兩端靠近端點處都取一相同距離l位置的溫度，靠近左端點的溫度為Ti1，靠近右端點的位置為Ti2。x為小溫區(qū)與起始點距離，假設小溫區(qū)間隙為5 cm。

爐內加熱過程如圖2所示。

圖2 爐內加熱過程

根據(jù)圖2，小溫區(qū)有以下四個一維穩(wěn)態(tài)熱傳導過程：爐前區(qū)域起始溫度分布、溫度相等的相鄰小溫區(qū)區(qū)域的溫度分布、溫度升高的相鄰小溫區(qū)區(qū)域的溫度分布、溫度降低的相鄰小溫區(qū)區(qū)域的溫度分布。

過程一：爐前區(qū)域與小溫區(qū)區(qū)域溫度分布。

研究物體內溫度分布（溫度分布為線性分布）：

過程二：溫度相等的相鄰小溫區(qū)區(qū)域的溫度分布。

根據(jù)熱傳導過程，相同溫度小溫區(qū)與間隙處于熱平衡狀態(tài)，小溫區(qū)釋放能量，間隙吸收能量，則間隙中點必然存在一個溫度最低點，即為間隙中點，間隙長度為5 cm，間隙離中點距離為2.5 cm，設該點與溫區(qū)的溫度差ΔT=5℃。

利用導熱方程通解算出平衡狀態(tài)下溫度分布：

過程三：溫度升高的相鄰小溫區(qū)區(qū)域的溫度分布。

根據(jù)熱傳導過程，低溫小溫區(qū)與高溫小溫區(qū)區(qū)域存在溫度差，同時相鄰小溫區(qū)區(qū)域會達到熱平衡穩(wěn)態(tài)。

邊界條件為：

代入導熱方程通解算出平衡狀態(tài)下溫度分布：

過程四：溫度降低的相鄰小溫區(qū)區(qū)域的溫度分布。

同過程三，高溫小溫區(qū)與低溫小溫區(qū)區(qū)域存在溫度差，同時相鄰小溫區(qū)區(qū)域會達到熱平衡穩(wěn)態(tài)。邊界條件為：

代入導熱方程通解算出平衡狀態(tài)下溫度分布：

4.1.2 非穩(wěn)態(tài)熱傳導過程分析

已知非穩(wěn)態(tài)有內熱源的導熱微分方程：

由于焊接點的厚度為0.15 mm，可忽略不計，則電路板焊接區(qū)域溫度在統(tǒng)一瞬間各點溫度基本相等，即t僅是的一元函數(shù)，與坐標x、y、z無關，即：

從而：

界面上交換的熱量應折算成整個物體的體積熱源，即-ΦV=Ah（T-T0）。

電路板進入回焊爐內時到達某一點時初始溫度與加熱溫度的關系式為：

整理一下得出以下關系式：

已知印刷電路板（PCB）體積，受熱表面積A，物性參數(shù)ρ、c、λ，熱交換系數(shù)h均為常數(shù)，令：

一維非穩(wěn)態(tài)導熱方程變?yōu)椋?/p>

4.2 模型的求解

4.2.1 問題一的求解

首先，分析整個過程可以得出：爐內傳送過程為穩(wěn)態(tài)熱傳導過程，集成電路板的瞬時加熱過程為非穩(wěn)態(tài)的熱傳導過程。假設傳送帶過爐速度為78 cm/min，各溫區(qū)溫度的設定值分別為173℃（小溫區(qū)1～5）、198℃（小溫區(qū)6）、230℃（小溫區(qū)7）和257℃（小溫區(qū)8～9）。

其次，將求解過程分為三個步驟：①將附件中的溫度曲線和加熱集成電路板的瞬時溫度與電路板的位置關系圖做出比較關系圖；②考慮所有過程的熱量損失問題，將熱傳導的兩個過程分別進行數(shù)值計算；③最后根據(jù)不斷調試的熱量參數(shù)h，得到調試后的爐溫曲線圖。

解：根據(jù)穩(wěn)態(tài)導熱方程，已經畫出焊接區(qū)域的大概溫度曲線，如圖3所示。

圖3 焊接區(qū)域的大概溫度曲線

結合一維非穩(wěn)態(tài)導熱方程與Matlab，根據(jù)題中所給的各參數(shù)的數(shù)值：175℃（小溫區(qū)1～5）、195℃（小溫區(qū)7）、255℃（小溫區(qū)8～9）及25℃（小溫區(qū)10～11），傳送帶的過爐速度為70 cm/min，由于δ的數(shù)值很難計算，可以先假設δ=0.006，可以看出當電路板位于0～25 cm，穩(wěn)態(tài)溫度曲線與爐溫曲線的擬合程度最高。

同理，可以得出不同位置范圍內實現(xiàn)穩(wěn)態(tài)溫度曲線與爐溫曲線擬合程度最高的一維非穩(wěn)態(tài)導熱方程系數(shù)：

通過改變不同位置段的δ值，將找到擬合程度最高的穩(wěn)態(tài)溫度曲線，如圖4所示。

圖4 穩(wěn)態(tài)溫度曲線

已知電路板（PCB）在進入回焊爐之前的溫度為T0=25℃，假設爐前區(qū)域的起始溫度取T00=26℃，爐前區(qū)域的起始點為研究起始點，電路板上焊接區(qū)域中心每隔0.1 s之后，移動的位移Δx=v·Δt=78 cm/min×0.1 s=0.013 cm，對照穩(wěn)態(tài)溫度曲線找出電路板移動Δx后的位置的初始溫度T（x），則Tx-T0為始末溫度差溫度，再結合一維非穩(wěn)態(tài)導熱方程，可以算出加熱溫度ΔT，從而解得這一位置焊接區(qū)域中心的溫度，該溫度也為下一位置的焊接區(qū)域中心未加熱前的初始溫度；利用迭代方法，從而求出每隔0.1 s后對應位置上焊接區(qū)域中心的溫度。

根據(jù)問題一中所給參數(shù)的數(shù)值，變換各溫曲溫度的設定值分別為182℃（小溫區(qū)1～5）、198℃（小溫區(qū)6）、239℃（小溫區(qū)7）和257℃（小溫區(qū)8～9），畫出爐溫曲線。已知傳送帶過爐速度為78 cm/min，爐溫曲線再結合一維非穩(wěn)態(tài)導熱方程，運用迭代法，計算得出每隔0.1 s后對應位置上焊接區(qū)域中心的溫度，結合Matlab畫出焊接區(qū)域中心的溫度變化曲線。分別計算小溫區(qū)3、6、7中點及小溫區(qū)8結束處距爐前區(qū)域起始點的位移xi（i=3，6，7，8），由算出電路板到達對應位置所用的時間，再根據(jù)焊接區(qū)域中心的溫度變化曲線找出小溫區(qū)3、6、7中點及小溫區(qū)8結束處焊接區(qū)域中心的溫度，結果如表2和圖5所示。

表2 焊接區(qū)域中心的溫度

圖5 焊接區(qū)域中心的溫度變化圖

4.2.2 問題二的求解

探究焊接過程，在滿足制程界限的情況下，要使以峰值溫度為中心線的兩側超過217℃的爐溫曲線應盡量對稱，需要算出合適的傳送帶過爐速度v和各溫區(qū)設定的溫度T有關，因此需要討論以下兩種情況。假設各溫區(qū)溫度的設定值分別為182℃（小溫區(qū)1～5）、203℃（小溫區(qū)6）、237℃（小溫區(qū)7）、254℃（小溫區(qū)8～9），傳送帶最大過爐速度v=79 cm/min。

第一種：探究同一次加熱情況下不同溫度對爐溫溫度的影響。

為了探究溫度的影響，采用控制變量法控制各溫區(qū)的溫度不變，恒定設置為182℃（小溫區(qū)1～5）、203℃（小溫區(qū)6）、237℃（小溫區(qū)7）、254℃（小溫區(qū)8～9），已知傳送帶的過爐速度調節(jié)范圍為65～100 cm/min，假設速度步長為0.5 cm/min，調節(jié)溫度大小。此時根據(jù)問題一求解方法結合Matlab，可以得到此時的爐溫曲線，從而得到爐溫峰值的溫度Tmax以及焊接區(qū)域中心到達峰值時刻tmax；再利用爐溫曲線上217℃對應的點得出升溫到217℃的時間ti1以及降溫到217℃的時間ti2，則算出升溫到217℃的點到爐溫峰值經歷的時間tmax-ti1以及爐溫峰值降溫到217℃的點經歷的時間ti2-tmax，算出爐溫峰值到左右217℃溫度的點之間的時間偏差Δti1=（tmax-ti1）-（ti2-tmax），則位時間偏差越接近于0，圖像越趨于對稱。其中i=1，2，3，…，70。

傳送帶的過爐速度與時間偏差關系如圖6所示。

圖6 傳送帶的過爐速度與時間偏差關系

可以直觀看出過爐速度越大，偏差值越小，所以要取符合條件的最大速度，即最大傳送帶過爐速度v=79 cm/min。

第二種：探究同一傳送帶速度下改變小溫區(qū)7、8的溫度對爐值溫度的影響。

在已知的允許的最大傳送帶過爐速度下，改變溫區(qū)7、8的溫度T7i、T8i，其他溫區(qū)的溫度為了方便計算，設置成問題二中的所給溫度。根據(jù)問題一求解方法結合Matlab，可以得到此時的爐溫曲線，類比第一種情況的方法，算出時間偏差Δt2i。

以溫區(qū)7為x軸，溫區(qū)8為y軸，時間偏差Δx為z軸，畫出三者變量的散點關系圖，如圖7所示。

圖7 散點關系圖

由圖7可以直觀地看出當溫區(qū)8溫度最高即T8=257℃，溫區(qū)7溫度最低即T7=225℃的情況下時間偏差最小，該爐溫曲線更符合題中所給要求的最優(yōu)情況。該題的指標為在不同傳送帶速度和各溫區(qū)的不同溫度下升溫到217℃的點到爐溫峰值經歷的時間tmax-ti1以及爐溫峰值降溫到217℃的點經歷的時間ti2-tmax的時間偏差，時間偏差越小，則越趨于對稱。

根據(jù)上述分析，將各溫區(qū)溫度的設定值分別為182℃（小溫區(qū)1～5）、203℃（小溫區(qū)6）、225℃（小溫區(qū)7）、257℃（小溫區(qū)8～9）。符合題意要求下允許的最大傳送帶過爐速度為v=79 cm/min時，區(qū)域最接近對稱，根據(jù)參數(shù)畫出如下爐溫曲線，如圖8所示。

圖8 焊接區(qū)域最接近對稱的爐溫曲線圖

此時曲線峰值溫度為245.9℃。

5 模型的應用和推廣

模型一中運用包含穩(wěn)態(tài)熱傳導過程與非穩(wěn)態(tài)熱傳導過程。穩(wěn)態(tài)熱傳導模型運用十分廣泛。本文穩(wěn)態(tài)熱傳導過程中主要運用一維熱傳導方程。該模型在生活中運用有密封材料的放熱、高溫作業(yè)衣物隔熱、一維穩(wěn)態(tài)傳熱過程類電流求解方法的討論、多層組件常壓隔熱性能實驗研究、純金屬型連鑄模型試驗研究，等等。

本題的整體研究方向為回焊爐爐溫曲線研究，旨在提高焊接工藝技術，延長機器使用，保障企業(yè)產品質量，對同物理原理的機器研究具有同樣適用性。