基于偏微分方程的高溫作業(yè)專用服裝設計模型

高振　房越　許乾宸

摘 要：長期在高溫環(huán)境下工作會對工作人員的生理及心理狀態(tài)產生影響，引發(fā)各種意外和危機，如何設計高溫作業(yè)專業(yè)服裝，已經成為人們日益關注的焦點。文章建立衣服各層材料的一維無內熱源熱傳導偏微分方程，結合各層交界處熱傳導平衡且溫度相同的條件，建立出溫度分布模型并進行分析。

關鍵詞：熱傳導偏微分方程；差分法；服裝設計

1 問題的提出及分析

根據(jù)2018高教社杯全國大學生數(shù)學建模競賽（CUMCM）題目A題[1]，當人們進行救火、高溫搶修等一系列高溫作業(yè)時，需要穿上特制的服裝來預防灼傷。這種服裝由3層復合織物組成，與外界環(huán)境直接發(fā)生接觸的為I層，III層與人體皮膚間的空隙為IV層。現(xiàn)取一個假人，將其體溫控制在37 ℃左右，放置在高溫環(huán)境中，最后通過測量，得到所使用的假人的皮膚外側溫度。試建立數(shù)學模型分析假人的皮膚外側的溫度變化狀況，該服裝的織物材料參數(shù)值在文獻[1]附件1中已給出，當II層厚6 mm、IV層厚5 mm、周圍環(huán)境的溫度是75 ℃、工作時長90 min時，實驗測得該假人的皮膚外側溫度在附件2中給出。本文試建立模型，計算溫度分布情況。

由于專用服裝內部的溫度分布是不均勻的，故一定會發(fā)生熱量從溫度相對高的地方向溫度相對低的地方轉移的現(xiàn)象，即熱傳導。考慮到在熱傳導過程中，服裝內無熱源，故建立一維無內熱源熱傳導偏微分方程，并找出初始條件和邊界條件，建立溫度分布模型，利用差分法求解，得到溫度分布情況[2]。

2 模型的建立與求解

2.1 基于偏微分方程的溫度分布模型的建立

顯然，專用服裝內部的溫度分布是不均勻的，考慮到在熱傳導過程中，專用服裝內無熱源，故建立一維無內熱源熱傳導偏微分方程，并找出初始條件和邊界條件，利用差分法進行求解，得到溫度分布情況。

2.1.1 熱傳導偏微分方程的建立

我們只考慮垂直于服裝表面方向上的熱傳導，故建立如下一維無內熱源熱傳導偏微分方程：

2.2 基于差分法的溫度分布模型的求解

為了得到該拋物型偏微分方程的數(shù)值解，我們采用差分法進行求解。差分法的基本思想是用只含有有限個未知量的、離散的差分方程組去代替連續(xù)的偏微分方程，將差分方程組的解作為偏微分方程的近似解。原方程和邊界條件中的微分替換為差分，將偏微分方程轉化為差分方程來求解。步驟如下。

2.2.1 區(qū)域離散

對距離-時間（x-t）平面進行劃分，將平面劃分為有限個格點組成的網格，這些離散點即為網格的節(jié)點。我們將距離區(qū)間[0，L4]劃分為N1等份，步長即為，則節(jié)點坐標為xi=ih（i=1，2，…，N1）。同樣地，將總時間區(qū)間[0，T]劃分為N2等份，步長即為，則節(jié)點坐標為tk=kτ（k=1，2，…，N2）。

2.2.2 近似代替

記點（xi，tk）處的函數(shù)值為uki=（i=1，2，…，N1；k=1，2，…，N2），用差商代替微商，則一維無內熱源熱傳導偏微分方程的差分形式為：

2.2.3 逼近求解

當t=0時，計算得到u0i-1，u0i，u0i+1，將上一層求解的邊界條件作為下一層的初始條件，并通過迭代法可以求解出u1i-1，u1i，u1i+1，進而求解出uki-1，uki，uki+1（i=1，2，…，N1；k=1，2，…，N2），運用Matlab進行求解，逐層計算。求出每層內表面的溫度，令k=1，根據(jù)各層交界處的溫度關系，可求得下一層的溫度分布，最終得到一維無內熱源熱傳導方程的數(shù)值解。

3 結語

通過以上的求解，我們利用Matlab繪制出了前100 s內各層材料的內表面溫度隨時間變化的曲線，并繪制出各層溫度關于時間和距離的三維分布圖（見圖2～3）。

由上述實驗結果得出如下結論。

（1）除I層在短時間內升溫特別迅速以外，其余3層溫度都處于緩慢上升過程，且越向內部溫度變化越慢。

（2）皮膚外側的溫度在開始的一段時間保持不變，合理解釋為最薄的一層熱導率也較高，離空氣越近，在較短時間內迅速升溫并向內傳導熱量，且傳導熱量需要一定時間，故最初一段時間內IV層的溫度變化并不顯著。

（3）在時間一定時，溫度隨著層數(shù)的增加逐漸減少，在相同層數(shù)內，溫度隨著時間的增加先增長迅速后趨于穩(wěn)定，最外層最終基本接近外界環(huán)境溫度。

[參考文獻]

[1]全國大學生數(shù)學建模競賽組委會.2018高教社杯全國大學生數(shù)學建模競賽（CUMCM）題目A題[EB/OL].（2018-09-13）[2019-02-25].http：//www.mcm.edu.cn/.

[2]司守奎，孫璽菁.數(shù)學建模算法與應用[M].長沙：國防工業(yè)出版社，2011.