劉倩嵐,雷億輝,肖育江,董亞玲
(吉首大學數學與統(tǒng)計學院,湖南 吉首 416000)
在高溫環(huán)境下,作業(yè)人員需要穿著特定的熱防護服,以避免熱源對人體造成的傷害,防止中暑、燒傷、灼傷等.因此,明確熱防護服的導熱規(guī)律、清楚高溫環(huán)境中工作時間的極限,從而建立科學、準確的熱防護服傳熱模型,制作出性能良好的專業(yè)防護服,可最大限度地保護作業(yè)人員的生命安全.學者對熱防護服內部熱量傳遞機理進行了研究,并建立了傳熱模型[1-6].但是,這些模型都需要求解微分方程,而在實際應用中,受邊界條件的約束,微分方程的求解十分不易.于是,筆者擬通過蒙特卡洛算法求出溫度分布模型中的相關參數,以提高建模效率.
本研究實驗數據來源于2018年高教社杯全國大學生數模競賽A題[7].將體內溫度控制在37 ℃的假人放置在實驗室的恒定高溫環(huán)境中,研究高溫環(huán)境下配備高溫作業(yè)專用服裝的人體皮膚表層的溫度變化情況.
根據實驗數據可知,導熱過程中人體皮膚表層的溫度隨時間的變化而變化,且隨著時間的推移逐漸趨于定值,由此可認為防護服的熱傳導是非周期性且非穩(wěn)態(tài)的.基于此給出以下假設:
(H1)熱防護服的紡織材料一、二、三層及空隙層分布均勻.
(H2)環(huán)境通過熱防護服及空隙向假人導熱是非周期性且非穩(wěn)態(tài)的.
(H3)熱量傳播的過程只有傳導沒有對流.
(H4)熱傳遞僅考慮熱傳導的影響,忽略汗液、水汽的影響,即不考慮濕傳遞.除熱量傳遞外,忽略環(huán)境其他因素的影響.
熱防護服的多層紡織面料導熱,可轉化為通過熱防護服一維橫截平板內的導熱.在第3類邊界條件下,由于熱防護服的特征厚度很小,邊界表面單位面積上的換熱熱阻遠大于材料內部單位導熱面積上的導熱熱阻,因此畢渥數Bi趨近于0,即
其中:d為材料層厚度;h為表面?zhèn)鳠嵯禂担沪藶閷嵯禂?進一步可知,在整個導熱過程中橫截平板各點的溫度基本一致,即同一層的紡織面料質量與熱容量都匯總到一點上,因此環(huán)境、熱防護服、假人皮膚表層之間的導熱分析可通過集中參數法進行.環(huán)境通過3層熱防護服紡織材料,以及熱防護服與人體間的空隙層(第Ⅳ層)將熱量傳遞至假人皮膚表層(圖1),具體傳熱機制如下:(1)在初始時刻,熱防護服3層材料溫度相等,為 37 ℃,不進行熱傳導;進入高溫環(huán)境后,熱量先從環(huán)境傳遞至熱防護服的第Ⅰ層,第Ⅰ層吸收熱量而升溫,再將熱量傳遞至第Ⅱ層,以此類推,熱量逐層傳遞.(2)置于高溫環(huán)境時,每一層材料的屬性不同,其吸收和放出的熱量也不同.環(huán)境傳遞至第Ⅰ層的熱量總大于第Ⅰ層傳遞至第Ⅱ層的熱量,以此類推.由此可知,熱量傳遞量是逐層遞減的.(3)接觸高溫環(huán)境后,第Ⅰ層材料外側流入熱量,隨著材料外側溫度升高,流入熱量逐漸減少;同時,材料內側溫度梯度增加,流出熱量逐漸增加,與流入熱量差距減少;當流出熱量等于流入熱量時,傳熱進入新的穩(wěn)態(tài)階段,達到平衡.其他材料層以此類推.
圖1 傳熱示意
由傅里葉導熱定律可得材料的熱流量計算公式為
(1)
其中:ΔΦ為導熱熱流量;ΔT為材料改變的溫度;A為熱傳遞面積.由(1)式可知,導熱率增大時當層材料熱流流量增加,因此降低材料的導熱率,可有效提高防護服的隔熱性能.
利用Matlab軟件將5 400 s的實驗數據轉換為直觀圖形(圖2).從圖2可以看出,假人皮膚表層溫度變化整體呈S型,變化趨勢符合阻滯增長模型的特點.為了準確建立相應的阻滯增長模型,對溫度變化數據進行具體分析:前15秒內,假人皮膚表層溫度變化不明顯;第16秒開始,假人皮膚表層溫度迅速上升;第427秒后,假人皮膚表層溫度上升速率由快變慢;第1 646秒后,假人皮膚表層溫度趨于恒定,保持在48.08 ℃.
圖2 0~5 400 s的假人皮膚表層溫度變化
根據材料間熱傳導機制可知,每層材料當前溫度T(t)等于前一刻溫度T(t-1)與傳入熱量而升高的溫度ΔT之和.ΔT的計算公式為
其中:ΔQ為熱量傳遞量;c為材料的比熱容;ρ為材料的密度.
熱傳導率表示單位截面和單位長度的材料在單位溫差和單位時間內直接傳導的熱量.引入熱傳導率,可得各層材料的熱量傳遞公式為
(2)
其中S表示材料的熱傳遞面積.
由此,可以初步建立如下溫度分布模型:
(3)
其中:T外為環(huán)境溫度;Tj(t)為各層材料的當前溫度;ΔTj(t)為每平方米各層材料改變的溫度.
由于第Ⅳ層材料是空氣,因此第Ⅳ層的溫度函數可以近似為假人皮膚表層的溫度函數.初始條件為 37 ℃的情況下,模型(3)的擬合數據曲線與實驗數據曲線的對比結果如圖3所示.
圖3 初步擬合對比結果
從圖3可以看出,當溫度恒定時,由模型(3)得到的溫度為75 ℃,遠高于實驗數據,且擬合曲線變化速率過快.由此可知,模型(3)擬合效果較差,需要對其進行修正.
考慮到人體散熱機制、溫差等因素對熱傳導的阻滯作用,對第Ⅳ層材料的溫度函數引入阻滯因子.由于阻滯值不明確,溫度函數的擬合效果不明顯,因此在阻滯因子中進一步嵌入控制參數來調整阻滯作用.
(4)
為了提高控制參數α的計算效率,使得均方差達到最小,這里采用蒙特卡洛算法進行求解.蒙特卡洛算法也稱為計算機隨機模擬方法,是指使用隨機數(或更常見的偽隨機數)來解決計算問題的方法.設計蒙特卡洛方法的初衷是用于解決物理數值模擬問題,后來隨著計算機的快速發(fā)展,這一方法廣泛應用在函數值極小化、計算幾何、組合計數等方面.如今,蒙特卡洛算法是求解科學、工程和科學技術領域大量應用問題的常用數值方法[8].利用蒙特卡洛算法進行求解,Matlab運行結果為α≈0.014 8.
利用(2)式,可得每單位各層材料熱傳遞到下一層同樣單位的材料能夠吸收的熱量(表1).
表1各層材料吸收的熱量
將ΔQj(j=1,2,3,4),Ts,α代入模型(4),可得假人皮膚表層溫度的擬合值.模型(4)的擬合數據曲線與實驗數據曲線的對比結果如圖4所示,部分時刻擬合數據與實驗數據的對比結果見表2.
圖4 修正后的擬合對比結果
表2 部分時刻的擬合數據與實驗數據的對比結果
從圖4和表2可知,擬合數據與實驗數據基本保持一致,說明模型(4)的擬合效果較好.擬合數據與實驗數據有一定偏差的原因是,材料初始溫度選取的是37 ℃,而現實環(huán)境中一開始材料溫度并未達到37 ℃.盡管存在擬合偏差,通過模型(4)設計的熱防護服依然能滿足現實要求.這是因為,由圖4可以看出,從第299秒開始擬合數據就略高于實驗數據,表明基于模型(4)設計的熱防護服實際達到的標準更高,即當熱防護服在該模型下滿足高溫作業(yè)要求時,現實中也一定能滿足防護要求.
結合熱傳導原理和傅里葉定律分析了防護服傳熱機制及材料層溫度、熱量的變化情況,建立了初步溫度分布模型.考慮到人體散熱機制等因素對熱傳導的阻滯作用,引入含控制參數的阻滯因子來修正初步溫度分布模型,并采用蒙特卡洛算法借助Matlab軟件求解相關參數,得到修正溫度分布模型.數值仿真實驗結果表明,修正溫度分布模型的擬合數據與實驗數據基本保持一致,說明該模型的擬合效果較好.