羅鳳潔 張玉林

摘? ?要：本文通過建立數(shù)學(xué)模型對(duì)三層織物材料組成的高溫作業(yè)專用服裝的設(shè)計(jì)問題進(jìn)行研究。運(yùn)用Fourier transform算法求解偏微分方程從而建立熱傳導(dǎo)模型，并利用實(shí)驗(yàn)數(shù)據(jù)擬合得到溫度的分布圖像與函數(shù)，得到最后的熱傳導(dǎo)模型公式。通過求解該模型的數(shù)值解，與實(shí)驗(yàn)數(shù)據(jù)進(jìn)行對(duì)比，進(jìn)行模型檢驗(yàn)。

關(guān)鍵詞：熱傳遞? 高溫作業(yè)專用服裝? ?Fourier transform? 數(shù)值擬合

中圖分類號(hào)：O242? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? 文獻(xiàn)標(biāo)識(shí)碼：A? ? ? ? ? ? ? ? ? ? ? ? 文章編號(hào)：1674-098X（2019）11（c）-0111-02

當(dāng)人們處于高溫環(huán)境下，人體的皮膚會(huì)受到一定程度的燒傷，因此，人們需要穿著一種專業(yè)服裝，在環(huán)境與人體之間充當(dāng)一種隔絕層進(jìn)而對(duì)人體皮膚進(jìn)行保護(hù)。人們對(duì)其經(jīng)濟(jì)實(shí)用性以及舒適性等都有了新的要求因此也進(jìn)行了相關(guān)的研究。本文運(yùn)用Fourier transform算法求解偏微分方程建立了熱傳導(dǎo)模型，并通過實(shí)驗(yàn)數(shù)據(jù)有效驗(yàn)證了模型的正確性，為今后人們研究更加高效經(jīng)濟(jì)的高溫作業(yè)熱防護(hù)服奠定一定的基礎(chǔ)。本文的研究背景以及研究數(shù)據(jù)均來源與大學(xué)生數(shù)學(xué)建模競賽2018年B題。

1? 問題提出

已知專用服裝材料的某些參數(shù)值例如密度、比熱、熱傳導(dǎo)率和厚度的值等。我們需要根據(jù)實(shí)驗(yàn)所得數(shù)據(jù)建立數(shù)學(xué)模型，并計(jì)算溫度分布。根據(jù)傳熱學(xué)與工程熱力學(xué)的相關(guān)知識(shí)可知，熱傳遞主要分為三種基本方式即熱傳導(dǎo)、熱輻射和熱對(duì)流，為了方便計(jì)算，我們暫時(shí)忽略以上因素對(duì)熱傳遞所造成的影響，主要考慮三大主因素的影響，從而根據(jù)傅里葉定律建立起熱傳導(dǎo)模型;得到溫度分布函數(shù)，進(jìn)而用理論的溫度分布函數(shù)與所得到的溫度分布函數(shù)進(jìn)行對(duì)比并檢驗(yàn)。

2? 符號(hào)定義

本文中令q為熱對(duì)流（輻射）密度;k為熱傳導(dǎo)系數(shù);σ為斯蒂芬-波爾滋蔓常數(shù);β為輻射吸收常數(shù);FL為向左輻射通量;FR為向右輻射通量;ε1為Ⅰ層左側(cè)的輻射量;h為服裝層或空隙層的厚度;ωf為纖維內(nèi)水蒸氣濃度。

3? 熱傳遞模型的建立與求解

3.1 模型建立

⑴在垂直于熱量傳導(dǎo)的方向上每個(gè)面的溫度都相同，我們稱之為等溫面;把環(huán)境到皮膚最外層表面的距離x分為n個(gè)小段，記為、、……、;每相鄰的之間的溫度差，記為（）;所以溫度梯度表示為。

⑵在考慮熱輻射的情況下，該服裝中的輻射能量即為熱能，記為在只考慮一維熱傳導(dǎo)的情況下，x點(diǎn)處向左輻射的總熱輻射量，并記熱流密度為q，由能量守恒得：。若在不考慮右熱輻射所傳遞的熱能的情況下，則在點(diǎn)x處的凈熱能量為：，其中β為輻射吸收常數(shù)即單位面積下該服裝的吸收能力。

綜合以上內(nèi)容的分析，得出問題一的模型：

（1）

其中：

3.2 模型求解

其中（1）看成兩個(gè)一階線性微分方程，它們的邊界條件為：

解得：

所以方程（1）進(jìn)一步寫為：

（2）

為解方程（3.2），我們需用到傅里葉積分變換，方程（2）變化為：

（3）

因此，方程（3）可看成一個(gè)含參量ω的一階線性非齊次常微分方程的初值問題，再利用常數(shù)變異法，再對(duì)取傅里葉逆變換，根據(jù)傅里葉變換卷積性質(zhì)得：

令：

得到：? ? ? ?（4）

其中分別表示Ⅰ層、II層、III層、IV層的溫度關(guān)于時(shí)間t、厚度x的二元函數(shù)并且。

4? 結(jié)果分析

圖1是用MATLAB軟件根據(jù)模型一的數(shù)值解模擬出來的三維圖像，反映了時(shí)間、溫度與厚度的關(guān)系。圖2反映了模型一的數(shù)值解以及附件2所給數(shù)據(jù)的對(duì)比情況，其中1線表示附件2的數(shù)據(jù)，2線表示模型一的數(shù)值解;從圖2中我們可以看出模型一的數(shù)值解與實(shí)際數(shù)據(jù)在前期接近程度較高，但隨著時(shí)間的推移，偏移量逐漸增加，因?yàn)樵谀Ｐ鸵恢袥]有對(duì)織物材料結(jié)構(gòu)參數(shù)中孔隙率、孔半徑等、熱濕傳遞以及物理參數(shù)如水汽擴(kuò)散等因素的影響，隨著時(shí)間的推移，溫度在不斷升高，而在高溫狀態(tài)下，以上因素對(duì)溫度降低的影響較大，因此，隨著溫度的升高，利用模型一算出的溫度將會(huì)比實(shí)際溫度偏低。

參考文獻(xiàn)

[1] 盧琳珍，徐定化，徐映紅.應(yīng)用三層熱防護(hù)服熱傳遞改進(jìn)模型的皮膚燒傷度預(yù)測(cè)[J].紡織學(xué)報(bào)，2018，39（1）.

[2] 上海交通大學(xué)數(shù)學(xué)系.數(shù)學(xué)物理方法[M].上海：上海交通大學(xué)出版社，2016.

[3] 潘斌.熱防護(hù)服裝熱傳遞數(shù)學(xué)建模及參數(shù)決定反問題[D].浙江理工大學(xué)，2017.

[4] 潘啟天，徐映紅，徐定華.基于熱濕耦合模型的三層 多孔織物厚度決定反問題[J].浙江理工大學(xué)學(xué)報(bào)：自然科學(xué)版，2017，37（6）：771-777.

[5] 盧琳珍.多層熱防護(hù)服裝的熱傳遞模型及參數(shù)最優(yōu)決定[D].浙江理工大學(xué)，2018.

[6] 余躍.紡織材料熱濕傳遞數(shù)學(xué)建模及其設(shè)計(jì)反問題[D].浙江理工大學(xué)，2016.

[7]蘇云，王云儀，李俊.消防服衣下空氣層熱傳遞機(jī)制研究進(jìn)展[J].紡織學(xué)報(bào)，2016，37（1）：67-172.