蔡彥 楊允出 錢江瑞

摘?要：通過分析織物內(nèi)部組織結構，將組織循環(huán)內(nèi)紗線分布按規(guī)律分6類，包括交織點紗線重疊區(qū)域、經(jīng)（緯）紗交替穿越區(qū)域、經(jīng)（緯）紗懸浮區(qū)域、孔隙區(qū)域。通過截面積相等的方式將不規(guī)則紗線等效轉化為規(guī)則的幾何體。通過分析組織循環(huán)內(nèi)各區(qū)域紗線與空氣分布情況，結合復合材料等效熱導率串并聯(lián)模型，建立各微元的等效熱導率數(shù)值模型。在已建立的模型基礎上進行算例分析，得到在維持其他紗線參數(shù)不變情況下，經(jīng)緯紗高度增加，織物整體等效熱導率增加;經(jīng)緯紗間距增加，織物整體等效熱導率降低;織物交織結構越復雜，懸浮區(qū)域數(shù)量占比越大，織物等效熱導率降低，其中懸浮區(qū)域數(shù)量為0的平紋織物等效熱導率相對最大。

關鍵詞：熱導率;機織物;數(shù)學模型;結構參數(shù)

中圖分類號：TS101.1

文獻標志碼：A

文章編號：1009-265X（2021）02-0043-07

作者簡介：蔡彥（1995-），男，浙江湖州人，碩士研究生，主要從事服裝熱舒適性方面的研究。

Abstract：The yarn distribution in weaving unit is classified into six categories according to the rule by analyzing the internal structure of fabric， including overlapping area of yarn at interlacing point， alternating crossing area of warp （weft） yarn， suspension area of warp （weft） yarn and pore area. The irregular yarns are transformed into regular geometries by means of equal cross-sectional area. Based on the analysis of yarn and air distribution in each area of the weaving unit and the series-parallel model of composite equivalent thermal conductivity， numerical models of equivalent thermal conductivity for each element were established. Analysis of examples based on the established models indicates that with other yarn parameters remaining unchanged， the overall equivalent thermal conductivity of fabrics increases as the yarn height increases; the overall equivalent thermal conductivity decreases as the yarn spacing increases; the equivalent thermal conductivity decreases as the complexity of the fabric interweaving structure increases and the proportion of suspended area increases. The equivalent thermal conductivity of plain weave fabrics with 0 suspended areas is relatively largest.

Key words：thermal conductivity; woven fabric; mathematical modeling; structural parameters

織物的傳熱性能是影響服裝熱舒適性的一個重要因素，在功能性服裝的設計及織物服用舒適性優(yōu)化中受到了廣泛關注?？椢?，無論是機織還是針織，都由包含多孔材質(zhì)紗線和孔隙的重復單元構成，其中機織物主要是由經(jīng)緯紗交織而成。為了探究機織物內(nèi)部傳熱機理，需對內(nèi)部紗線及空隙的幾何結構、各級材料物理屬性與織物的傳熱性能的關系進行深入分析、建模。許多學者基于統(tǒng)計模型、神經(jīng)網(wǎng)絡模型、有限元方法，分形模型、物理建模等開展了關于纖維集合體或織物的內(nèi)部結構對其傳熱性能影響、預測研究。吳海軍等[1]利用相關性分析法，對不同結構參數(shù)下織物的傳熱性能進行了研究。李麗等[2]基于兩相復合介質(zhì)的等效熱導率模型，運用順變平面熱源法對環(huán)氧樹脂-纖維集合體進行測試，得出其等效熱導率，并驗證了其各向異性。肖紅[3]通過構建瞬態(tài)以及穩(wěn)態(tài)條件下織物傳熱方程，并通過制備樣品測試其傳熱性能，驗證了所建立方程的適用性。也有通過統(tǒng)計以及神經(jīng)網(wǎng)絡模型來分析各變量的對傳熱性能、熱阻的影響。Bhattacharjee等[4]將各類織物參數(shù)進行組合，通過響應面統(tǒng)計分析法，探究各參數(shù)集體效應對熱阻的影響，由于大部分織物結構參數(shù)為相互耦合關聯(lián)，因此很難獲得單一參數(shù)與熱阻傳熱性能的明確關系。Bhattacharjee等[5]通過建立人工神經(jīng)網(wǎng)絡模型的方式來對qamx值和熱阻進行模擬預測，雖然得到了預期的結果，但前期需要大量的數(shù)據(jù)對模型進行訓練，才能得到更為準確的結果。也有學者通過對傳熱現(xiàn)象的數(shù)學分析，進行簡化處理總結成物理模型，采用物理模型來對傳熱原理進行表征。Ismail等[6]給出了可預測平紋織物熱導率的理論模型，且該模型也可用于表征纖維在織物、紗線、纖維狀態(tài)下的等效熱導率。Zheng等[7]通過電鏡圖像實測機織物的組織結構參數(shù)和紗線交織規(guī)律，利用Texgen建立了玻璃纖維平紋織物幾何模型，利用有限元軟件開展了織物熱傳遞數(shù)值模擬研究。Zhu等[8]通過建立孔隙結構與多孔介質(zhì)結構參數(shù)函數(shù)的方式，探究了含不飽和水分的纖維狀多孔介質(zhì)的有效導熱模型，預計此模型可有助于評價紡織品的熱舒適性。

以上關于織物傳熱性能的研究大都著重于傳熱方式以及各參數(shù)耦合條件下織物傳熱性能的變化情況。結合已有研究，擬從織物內(nèi)部紗線幾何參數(shù)與等效熱導率的影響出發(fā)對織物傳熱進行探究。以機織物為研究對象，通過對織物組織結構及紗線幾何形態(tài)的分析，對織物組織內(nèi)部的紗線排列按規(guī)律進行分類;對每一類微元進行傳熱分析，結合復合材料等效熱導率串并聯(lián)模型，建立便于計算的等效熱導率數(shù)值模型。并采用建立的數(shù)值模型，再從織物內(nèi)紗線高度、間距及交織規(guī)律三個影響條件下對熱導率變化情況進行算例分析。

1?織物等效導熱系數(shù)建模

1.1?機織物組織循環(huán)的結構特征及區(qū)域劃分

對于機織物結構形態(tài)的描述通常有以下幾類經(jīng)典幾何模型：圓（橢圓）形、跑道模型、透鏡形[9]。本文以基本組織機織物為研究對象，對組織循環(huán)單元的紗線交織結構規(guī)律進行分析和區(qū)域細分。以2上1下左斜紋為例，將該組織循環(huán)劃分為交織點紗線重疊區(qū)域，紗線交替穿越區(qū)域、紗線懸浮區(qū)域和孔隙區(qū)域四類（見圖1、圖2）。由圖1可知，該斜紋組織循環(huán)共被劃分為36個基本微元區(qū)域，其中包括9個交織點紗線重疊區(qū)域（經(jīng)組織點6個、緯組織點3個），12個紗線交替穿越區(qū)域（經(jīng)紗交替穿越區(qū)6個、緯紗交替穿越區(qū)6個），6個紗線懸浮區(qū)域（經(jīng)紗懸浮區(qū)域3個、緯紗懸浮區(qū)域3個），9個孔隙區(qū)域。為了后續(xù)構建數(shù)學模型，進一步設定該組織的相關結構參數(shù)，包括經(jīng)、緯紗的寬度a1、a2，經(jīng)、緯紗的高度b1、b2，相鄰兩經(jīng)紗的間距c1，相鄰兩緯紗的間距c2。

橢圓織物幾何模型中，經(jīng)緯紗線的截面均設為橢圓形，經(jīng)緯紗線的中心軸線在交替區(qū)域采用了正弦函數(shù)曲線。為了簡化計算，紗線中心軸線的曲線形狀轉換為折線形狀（見圖2），將經(jīng)緯紗線的橢圓截面形狀采用寬度不變、面積相等的原則轉換為矩形形狀（見圖3）。

等效轉化后的經(jīng)、緯紗線高度b1′，b2′取值分別見式（1）、式（2）。

1.2?基于并聯(lián)模型的機織物組織循環(huán)各區(qū)域熱導率計算

織物內(nèi)部熱量傳遞的方式有熱傳導、熱對流與熱輻射，大多數(shù)條件下，熱對流與熱輻射傳遞的熱量相較于熱傳導可忽略，因此本文中對于織物內(nèi)部傳熱的分析只考慮熱傳導?？椢锝M織循環(huán)內(nèi)各微元的等效熱導率采用復合材料的串并聯(lián)模型進行計算。當熱量傳導的方向與兩種材料交接面平行時，為兩相并聯(lián)模型，總熱量等于分別流過兩種材料熱量之和，此時等效熱導率計算見式（3）。當傳熱方向與兩種材料交接面垂直時，熱量依次通過兩種材料，為兩相串聯(lián)模型，此時等效熱導率計算見式（4）。

1.2.1?交織點紗線重疊區(qū)域的熱導率計算

交織點重疊區(qū)域是經(jīng)緯紗交織的區(qū)域，此區(qū)域中包含兩根紗線，經(jīng)緯紗見圖4。由圖4可知，紗線交織點處可視作空氣-經(jīng)紗以及空氣-緯紗串聯(lián)的集合體，因此采用等效熱導率串聯(lián)模型分別計算空氣-經(jīng)紗的等效熱導率λ1′，見式（5）;空氣-緯紗的等效熱導率λ2′，見式（6）。

圖5為紗線懸浮區(qū)域紗線分布。由圖5可知，在交織點處，經(jīng)緯紗呈串聯(lián)，結合式（4）的串聯(lián)模型，可得交織點處的等效熱導率計算方法，見式（7）。將λ1′，λ2′代入式（7）可得交織點處等效熱導率取值，如式（8）所示。

1.2.2?紗線懸浮區(qū)域熱導率計算

經(jīng)緯紗懸浮區(qū)域中只含單根的紗線，其他區(qū)域均為空氣，并且懸浮區(qū)域的兩側經(jīng)緯紗上下排列情況相同，見圖5，以經(jīng)紗為例，但又因經(jīng)緯紗幾何參數(shù)以及交織規(guī)律的不同而有所不同。由圖5可知，此區(qū)域可視為紗線與空氣的串聯(lián)，根據(jù)串聯(lián)模型可得經(jīng)紗懸浮區(qū)域等效熱導率λya的計算方法見式（9）;同理，可推得緯紗懸浮區(qū)域等效熱導率取值λxa，見式（10）。

式中：λxa、λya分別為經(jīng)紗懸浮區(qū)域、緯紗懸浮區(qū)域的等效熱導率，W/（m·K）;λ1，λ2，λa分別為經(jīng)紗、緯紗、空氣的熱導率，W/（m·K）;b1、b2為紗線的高度，mm。

1.2.3?紗線交替區(qū)域熱導率計算

紗線交替穿越區(qū)域中同懸浮區(qū)域類似，也只含一根紗線，但該紗線斜向穿越此區(qū)域，使得此區(qū)域兩側經(jīng)緯紗上下排列順序相反，見圖6，此處以經(jīng)紗為例。由圖6可知，此區(qū)域內(nèi)可視作是空氣與紗線的串聯(lián)，通過對等效轉化后紗線幾何形態(tài)的分析，運用式（4）計算，得到經(jīng)紗交替穿越區(qū)域等效熱導率λy的計算方法見式（11）;同理，推得緯紗交替穿越區(qū)域等效熱導率λx的計算方法見式（12）。

1.3?織物單元總等效熱導率計算

以上完成了對機織物組織循環(huán)內(nèi)各基本微元等效熱導率計算公式的推導。而織物組織循環(huán)是由上述6類微元進行排列組合構成的。即可將每一組織循環(huán)視為多個基本微元的并聯(lián)。因此計算織物整體等效熱導率只需分別計算得到各基本微元體積占比，以及微元數(shù)量。由圖1可知，體積占比用織物俯視圖中各基本微元的面積占比表示。以圖1中2上1下左斜紋織物為例，分別計算組織循環(huán)經(jīng)緯向的長度l1，l2以及組織循環(huán)面積s，計算方法見式（13）-式（15）。

式中：l1，l2分別為組織循環(huán)經(jīng)緯向的長度，mm;m1，m2，m3，m4分別為單根經(jīng)緯紗上各模塊的個數(shù)，s為組織循環(huán)的面積，mm2。

因此各模塊微元的體積占比分別為：

組織循環(huán)內(nèi)各微元數(shù)量與機織物的組織結構有關，結合式（4）復合材料等效熱導率的并聯(lián)模型，織物組織循環(huán)即織物的等效熱導率計算方法見式（16）。

式中：n為組織循環(huán)內(nèi)各個模塊的數(shù)量，可通過確定機織物的組織循環(huán)的交織規(guī)律來確定。以圖1中2上1下左斜紋織物為例，此時nxy取值為9、na的取值為9、nx與ny的取值為6、nxa和nya的取值為3。

2?結果與討論

由上述所建立等效熱導率數(shù)值計算模型可知，影響織物等效熱導率的主要是兩大因素：為織物內(nèi)紗線交織的方式與紗線的幾何參數(shù);織物內(nèi)部紗線以及空隙區(qū)域的空氣（或填充物）的熱導率。針對上述兩類影響條件，選取圖1中2上1下左斜紋織物作為算例進行分析，織物幾何參數(shù)見表1。從織物內(nèi)部紗線間距，紗線高度以及紗線交織規(guī)律3種影響因素條件下織物等效熱導率的變化進行分析。

2.1?紗線高度對等效熱導率的影響

首先改變經(jīng)紗高度，使經(jīng)紗在原始參數(shù)基礎上以0.02 mm為一個梯度逐步增加，采用已建立的數(shù)值模型計算組織循環(huán)內(nèi)各區(qū)域等效熱導率以及整體等效熱導率，結果見圖7。

由圖7可知當經(jīng)紗高度增加時，只包含單根經(jīng)紗的懸浮區(qū)域熱導率會隨紗線高度的增加而逐步增大，經(jīng)紗交替穿越區(qū)域熱導率也同樣隨之增加，且增加幅度更為明顯，緯紗交替穿越區(qū)域與緯紗懸浮區(qū)域熱導率呈現(xiàn)下降趨勢，交織點等效熱導率未發(fā)生明顯變化，織物整體等效熱導率呈上升趨勢。這是因為紗線高度增加使得織物內(nèi)交替穿越區(qū)域與懸浮區(qū)域內(nèi)紗線所占體積占比增加，從而使得等效熱導率增加，反映到織物整體即為織物整體等效熱導率增加。

同時增加經(jīng)緯紗高度，使經(jīng)緯紗高度以同樣的梯度增加，代入已建立的數(shù)值模型計算，組織循環(huán)內(nèi)各區(qū)域模塊的等效熱導率以及整體等效熱導率計算結果見圖8。

由圖8可知，當經(jīng)緯紗高度同時增加時，交織點紗線重疊區(qū)域、經(jīng)緯紗懸浮區(qū)域等效熱導率幾乎維持恒定，交替穿越區(qū)域等效熱導率則呈現(xiàn)增長趨勢，織物整體等效熱導率呈現(xiàn)增長趨勢。這是因為經(jīng)緯紗懸浮區(qū)域中，經(jīng)緯紗高度同時增加時，其間紗線與空氣的體積占比變化幅度很小，幾乎保持不變，而經(jīng)緯紗交替穿越區(qū)域中紗線體積占比會因為織物內(nèi)部紗線高度增加而增加，導致等效熱導率計算結果增加。紗線高度的變化不會改變紗線內(nèi)各區(qū)域體積占比，因此可得紗線整體等效熱導率呈現(xiàn)增大的趨勢。

2.2?紗線間距對等效熱導率的影響

改變經(jīng)向紗線間距，使得紗線間距以0.04為一梯度逐級增加，采用已建立的數(shù)值模型計算組織循環(huán)內(nèi)各區(qū)域模塊的等效熱導率以及整體等效熱導率，計算結果見圖9。

由圖9可知，當紗線經(jīng)向間距增加時交織點區(qū)域、經(jīng)緯紗懸浮區(qū)域等效熱導率均趨于平穩(wěn)無較大變化，經(jīng)紗交替穿越區(qū)域等效熱導率呈下降的趨勢，織物整體等效熱導率也呈現(xiàn)下降趨勢。這是因為紗線間距的增加，使得紗線分布更為疏松，交替穿越區(qū)域內(nèi)紗線體積占比降低，經(jīng)紗交替穿越區(qū)域等效熱導率降低;且在組織循環(huán)內(nèi)純空氣區(qū)域、經(jīng)緯紗交替穿越區(qū)域、懸浮區(qū)域的體積占比增加，交織點區(qū)域的體積占比降低，因此織物整體熱導率呈現(xiàn)下降趨勢。

同時改變經(jīng)緯紗線間距，使兩根紗線間距以同樣的梯度增加，代入已建立的數(shù)值模型計算，組織循環(huán)內(nèi)各區(qū)域模塊的等效熱導率以及整體等效熱導率計算結果見圖10。

由圖10可知，當經(jīng)緯紗間距同時增加時，交織點處以及經(jīng)緯紗懸浮區(qū)域的等效熱導率幾乎維持恒定，交替穿越區(qū)域等效熱導率則呈現(xiàn)降低趨勢。由于經(jīng)緯紗間距增加，交替穿越區(qū)域內(nèi)紗線體積占比降低，交替穿越區(qū)域等效熱導率降低。并且經(jīng)緯紗間距的增加會使得純空氣區(qū)域、經(jīng)緯紗交替穿越區(qū)域、懸浮區(qū)域的體積占比增加。上述兩因素共同導致了織物整體等效熱導率降低。

2.3?織物交織規(guī)律對等效熱導率的影響

由圖1可知，織物交織規(guī)律的變化會影響組織循環(huán)的大小以及組織循環(huán)內(nèi)各微元的數(shù)量，各微元數(shù)量的變化會影響織物內(nèi)紗線與空氣體積占比。在保持紗線參數(shù)不變的情況下，以平紋織物、2上1下左斜紋織物、2上2下方平織物為例，通過對組織結構的分析，分別計算其等效熱導率，結果見表2。

由表2可知，織物類型變化時，組織循環(huán)內(nèi)各區(qū)域微元總數(shù)量會變多，交織點和空氣區(qū)域的數(shù)量占比均保持恒定，為總數(shù)的1/4。交替穿越區(qū)域與懸浮區(qū)域會隨著組織結構的變化而發(fā)生變化，如平紋織物的經(jīng)緯向交替穿越區(qū)域數(shù)量均為4，占微元總數(shù)的1/4，且不存在懸浮區(qū)域;2上1下左斜紋織物經(jīng)緯向交替穿越區(qū)域的數(shù)量均為6，占微元總數(shù)的1/6，經(jīng)緯向懸浮區(qū)域數(shù)量為3，占微元總數(shù)的1/12;2上2下方平織物經(jīng)緯向交替穿越區(qū)域的數(shù)量均為8，占微元總數(shù)的1/8，經(jīng)緯向懸浮區(qū)域數(shù)量為8，占微元總數(shù)的1/8。

織物內(nèi)部組織類型的變化主要表現(xiàn)為交替穿越區(qū)域與懸浮區(qū)域數(shù)量的變化，由表2可知當懸浮區(qū)域數(shù)量占比增大時紗線的等效熱導率會降低，平紋織物不存在懸浮區(qū)域，因此此時紗線交織最為緊密，此時織物等效熱導率為三者中最大;2上2下方平織物懸浮區(qū)域數(shù)量占比最大為1/8，此時織物等效熱導率為三者中最小。懸浮區(qū)域等效熱導率為小于交替穿越區(qū)域的熱導率，因此懸浮區(qū)域數(shù)量占比的增大會導致織物等效熱導率降低。

2.4?織物熱導率實驗測試

基于文中對于熱導率的探究，對用于算例分析的面料進行熱導率測試，并與計算結果進行對比，驗證模型可靠性。使用平板式保溫儀測試織物的傳熱系數(shù)，通過計算得到織物熱導率。在織物上無開邊、破洞處裁取5個30 cm×30 cm的試樣，分別進行測試，并取平均值，得到平均傳熱系數(shù)為55.43 W/（m2·K），結合織物厚度0.64 mm可計算得到織物熱導率為0.035 475 W/（m·K）。采用文中建立的數(shù)值模型計算得到的等效熱導率為0.031 442 W/（m·K），兩者間誤差值大約為11.37%，在可接受范圍內(nèi)，其中誤差來源主要是在對于紗線進行等效轉化的過程中，紗線體積發(fā)生了變化，不能做到完全等效，但熱導率計算值與實驗值在量級以及數(shù)值上均較為接近，因此所建立的等效熱導率模型可用于對織物等效熱導率進行模擬分析。

3?結?論

通過對機織物組織循環(huán)進行分析，根據(jù)紗線排列方式將組織循環(huán)內(nèi)各區(qū)域分成6類微元，分別為交織點紗線重疊區(qū)域、經(jīng)（緯）紗交替穿越區(qū)域、經(jīng)（緯）紗懸浮區(qū)域、孔隙區(qū)域。通過截面積相等的方法，將不規(guī)則的紗線轉化為規(guī)則的矩形、折線形進行計算。以復合材料等效熱導率模型為基礎，對各微元進行分析，并建立等效熱導率計算模型;在此基礎上將各微元的等效熱導率模型進行整合得到織物等效熱導率模型。由模型可知，在只考慮織物內(nèi)部紗線層面的影響因素時，影響織物等效熱導率的主要因素為織物內(nèi)部紗線與空氣體積占比與紗線以及空氣（填充物）的熱導率。其中紗線的幾何參數(shù)以及交織規(guī)律的變化會影響織物組織循環(huán)內(nèi)各類微元數(shù)量以及體積占比，從而影響紗線與空氣的體積占比。保持織物內(nèi)部其余變量均保持不變，只考慮單一變量情況時，經(jīng)緯紗高度增大，會使得織物等效熱導率增加;經(jīng)緯紗間距增大，會使得織物等效熱導率降低;織物組織循環(huán)中懸浮區(qū)域數(shù)量占比增加（及交織規(guī)律變化）會使得織物等效熱導率降低，當懸浮區(qū)域數(shù)量為0是此時織物為平紋組織，等效熱導率為同等紗線參數(shù)條件下最大。