楊恩惠 初曦 邱華

摘要： 為了提出一種用于優(yōu)化、評估、預(yù)測針織物導(dǎo)熱性能的方法，文章探索了緯編針織物熱量傳遞的有限元仿真模擬，尤其是夏季汗衫面料，因為其導(dǎo)熱性能與人體舒適性密切相關(guān)。首先使用織物厚度儀和超景深顯微鏡獲取織物結(jié)構(gòu)參數(shù)，通過插入三次B樣條曲線建立緯平針織物的三維幾何模型。然后將單元線圈模型導(dǎo)入有限元分析軟件STARCCM+中，通過網(wǎng)格劃分、邊界設(shè)置、迭代計算步驟進行數(shù)值計算，并和實驗值進行了比較分析，誤差在4%以內(nèi)。表明此方法具有一定的實用性，可用于預(yù)測緯編針織物的導(dǎo)熱性能。

關(guān)鍵詞： 緯編針織物;導(dǎo)熱性能;三維模型;有限元仿真;網(wǎng)格劃分

中圖分類號： TS186.2文獻標(biāo)志碼： A文章編號： 10017003（2020）01003106

引用頁碼： 011106DOI： 10.3969/j.issn.10017003.2020.01.006

Finite element simulation of thermal conductivity of the weft knitted fabric

YANG Enhui， CHU Xi， QIU Hua

Abstract： In order to propose a method for optimizing， evaluating and predicting the thermal conductivity of knitted fabrics， the finite element simulation of heat transfer of weft knitted fabrics was explored in this paper. Especially for summer Tshirt fabric， its thermal conductivity is closely related to the comfort of human body. Firstly， fabric structure parameters were obtained by means of fabric thickness meter and ultradepthoffield microscope. And threedimensional geometric model of weftflat knitted fabric was established by inserting cubic Bspline curve. Then， the element loop model was imported into the finite element analysis software STARCCM+. The numerical calculation was carried out by mesh division， boundary setting， and iterative calculation steps. Compared with the experimental data， the error was less than 4%. It is shown that this method is practical and can be used to predict the thermal conductivity of weft knitted fabrics.

Key words： weft knitted fabric; thermal conductivity; threedimensional model; finite element simulation; mesh division

服裝織物的熱舒適性一般指人體在不同環(huán)境條件下或者不同活動條件下，外界環(huán)境與內(nèi)在皮膚產(chǎn)生熱量的交換，達到平衡時使人感到舒適的特性。實驗條件下，衡量織物熱舒適性能的指標(biāo)一般采用克羅值、熱阻或傳熱系數(shù)等[1]。

在研究織物熱傳導(dǎo)的過程中，以往大多研究其影響因素，宋蘭霞等[2]測試不同經(jīng)緯密的織物來研究結(jié)構(gòu)密度、厚度、含氣量對織物導(dǎo)熱性能的影響;吳海軍等[3]選擇不同組織結(jié)構(gòu)的織物，利用線性分析織物厚度、平均浮長、覆蓋系數(shù)對織物熱傳遞性能的影響。劉讓同等[4]采用織物疊加的方式探討厚度、容重等對熱阻的影響，建立了擬合回歸曲線，其二者組合參數(shù)時相關(guān)度R2=0.9872。王婷婷等[5]采用熱線法測試幾種高性能織物在不同溫度環(huán)境下的導(dǎo)熱性，并利用最小二乘法擬合出相關(guān)性曲線。后來隨著模擬軟件的發(fā)展，吳佳佳等[6]采用有限元軟件ABAQUS模擬平紋機織物的熱傳遞過程。戈夢夢[7]利用數(shù)值模擬計算探討出平紋機織物的熱阻與織物厚度和孔隙率有關(guān)。Phillip Gibson[8]建立圓筒狀的織物模型，研究環(huán)境風(fēng)速、溫度和濕度對傳熱性質(zhì)的影響。

關(guān)于針織物的熱傳遞性能，趙超等[9]在電腦橫機上編織了幾種不同結(jié)構(gòu)參數(shù)的緯平針織物，探究彎紗深度對織物導(dǎo)熱系數(shù)的影響。孫玉釵[10]采用熱量傳遞數(shù)學(xué)模型和有限元計算來分析針織內(nèi)衣織物熱傳遞性能，將針織絮類夾層襯墊看做是均勻連續(xù)體進行研究。但實際上，針織物中紗線相互圈套，含有不規(guī)則空氣層，所以本文針對針織物結(jié)構(gòu)建立織物幾何模型，并導(dǎo)入有限元分析軟件中對熱傳遞性能進行計算，通過與實際測試結(jié)果對比分析，建立一種評估、優(yōu)化和預(yù)測針織物導(dǎo)熱性能的有效方法。

1織物模型建立

1.1織物幾何結(jié)構(gòu)

實驗樣布采用的是緯平針織物，因為緯平針組織為針織的最基本組織，且是日常汗衫的常用組織，對于研究針織物熱舒適性能有一定的代表性。樣布的材料是棉，因為棉織物輕柔、舒適，且是熱量的不良導(dǎo)體，汗衫采用棉材料為了吸濕導(dǎo)汗、柔和貼膚。所以本實驗采用三種不同紗線細(xì)度、組織密度的棉緯平針織物作為樣布，來研究針織物熱傳遞性能，織物規(guī)格如表1所示。

1.2構(gòu)建幾何模型

通過YG141D織物厚度測試儀（常州第二紡織機械有限公司）和VHX5000超景深數(shù)碼顯微鏡（基恩士公司）對織物尺寸進行測量，并且測量多次求其平均值，得到織物尺寸參數(shù)如表2所示。圖1為型值點測量位置，其中a、b、c、α、β均為在XY平面內(nèi)測量，a為點P1到P2的距離，b為P2到P3的距離，c為P4到P5的距離，α為P2P3與P3P4的夾角，β為P1P2與P2P3之間的夾角。

按照測量所得的數(shù)據(jù)可以確立P1～P4，那么P5～P8關(guān)于線段P4P5中心線對稱，所以依次確立8個型值點，然后采用非均勻有理B樣條曲線構(gòu)建紗線路徑，樣條曲線表達式為[11]：

P（u）=∑ni=oWiQiBi，k（u）∑ni=0WiBi，k（u）（1）

式中：Qi表示第i個幾何控制頂點;Wi是相對應(yīng)幾何控制頂點的權(quán)值，用來調(diào)節(jié)控制頂點的權(quán)重;k為B有理樣條曲線的冪次;Bi，k（u）為k次B樣條基函數(shù)。

本文中k取3，采用積累弦長參數(shù)化法[1213]計算樣條曲線的矢量節(jié)點，對數(shù)據(jù)點進行嚴(yán)格化規(guī)范后可得到一個參數(shù)化序列，將其分別對應(yīng)的矢量點代入到非均勻有理三次B樣條曲線函數(shù)中，并且使曲線經(jīng)過幾個特殊的插值型值點，以及滿足曲線的首尾兩端的端點與其型值點重合的條件，聯(lián)立求解即可得出曲線的控制點。然后經(jīng)過陣列重復(fù)可得到緯平針組織模型，如圖2所示。

2織物有限元模擬

2.1模型前處理

模型在導(dǎo)入到有限元分析軟件之前需要先切割成需要的結(jié)構(gòu)單元，因為既要代表性地反映整塊織物的熱傳遞性能，又要避免模型較大導(dǎo)致計算量大。所以，本實驗中將緯平針組織模型切割成單個完整線圈模型，然后將切割好的模型導(dǎo)入到有限元分析軟件STARCCM+中，運用布爾運算將固體紗線組合成一個整體，便于網(wǎng)格劃分和材質(zhì)的屬性定義，并在固體模型周圍建立空氣流體區(qū)域，然后進行網(wǎng)格劃分。網(wǎng)格是模擬計算的基礎(chǔ)，網(wǎng)格的大小和質(zhì)量影響迭代計算速度和最后收斂情況[14]。由于該模型尺寸小、表面屈曲程度較大，所以網(wǎng)格劃分時需要特別注意表面連接處的網(wǎng)格質(zhì)量。利用軟件對表面網(wǎng)格質(zhì)量進行檢查，一般情況下，檢查面穿透、面質(zhì)、接近面、自由邊、非多重邊和非多重面等幾項標(biāo)準(zhǔn)，利用自動修復(fù)和手動修復(fù)相結(jié)合，將問題網(wǎng)格一一修復(fù)。然后根據(jù)面網(wǎng)格可選擇劃分四面體網(wǎng)格、多面體網(wǎng)格或切割體網(wǎng)格，因為多面體網(wǎng)格有較多的相鄰單元，適用于表面屈曲程度大的幾何體，并且其梯度計算和流動模擬更為準(zhǔn)確，所以本實驗中采用多面體網(wǎng)格。

實驗中所建立的三個幾何模型，在STARCCM+有限元軟件中，經(jīng)過表面修復(fù)和自動體網(wǎng)格劃分，網(wǎng)格尺寸為0.01mm，最終的網(wǎng)格數(shù)量和質(zhì)量情況分別是：32s棉織物的固體區(qū)域147899cells，流體區(qū)域286070cells;網(wǎng)格具有拓?fù)溆行裕嬗行詾?00，體積變化率在0.1～1.00為99.015%。26s棉織物的固體區(qū)域175229cells，流體區(qū)域312479cells;網(wǎng)格具有拓?fù)溆行?，面有效性?00，體積變化率在0.1～100為99.039%。21s棉織物的固體區(qū)域241105cells，流體區(qū)域464675cells;網(wǎng)格具有拓?fù)溆行?，面有效性?.00，體積變化率在0.1～1.00為99.156%。圖3（a）為劃分網(wǎng)格后的織物紗線模型，圖3（b）為劃分網(wǎng)格后透明化的整個模型。

2.2邊界條件設(shè)置

將紗線區(qū)域定義為固體域，物理模型采用三維恒密固體域，并設(shè)置分離固體焓變等屬性條件，固體材料屬性依據(jù)棉纖維材料進行界定，其中密度850.0kg/m3，導(dǎo)熱系數(shù)0.072W/（m·K），比熱值1275J/（kg·K）。將外圍空間定義為流體域，物理模型采用三維層流氣體域，并設(shè)置恒密度、分離流、分離流體焓變等屬性條件?？諝饬黧w的屬性可設(shè)置為默認(rèn)數(shù)值，密度1.184kg/m3，導(dǎo)熱系數(shù)0026W/（m·K），比熱值1003.62J/（kg·K）。邊界條件的設(shè)置如圖4所示，將流體域上表面top和下表面bottom的邊界設(shè)置為壓力出口pressure outlet，而四周的面設(shè)置為壁面，圖4模型邊界

Fig.4Model boundary不進行質(zhì)量和熱量交換。其中上表面的初始溫度設(shè)置為20℃，下表面的初始溫度設(shè)置為35℃，這是因為在實驗條件下，下表面所接觸的熱板為35℃，氣候室的溫度為20℃。

因為模型上下表面存在溫度差，所以流體中的分子運動和固體中的聲子運動可實現(xiàn)能量和質(zhì)量的傳遞[15]，傳遞過程中遵循控制方程如下[16]：

質(zhì)量守恒方程：

ρt+xi（ρui）=0（2）

動量守恒方程：

tρui+xjρuiuj=-pxi+τijxj+ρgi+Fi（3）

能量守恒方程：

tρh+xiρuih=xik+kiTxi+Sh（4）

式中：ρ為流體密度，t為時間，p為靜壓力，h為熵，k為分子傳導(dǎo)率，ki為由湍流傳遞而引起的傳導(dǎo)率，ui為流體速度沿i方向的分量，uj為流體速度沿j方向的分量，xi為i方向的坐標(biāo)，xj為j方向的坐標(biāo)，τij為應(yīng)力矢量，gi為i方向的重力分量，T為溫度，Sh為定義的體積源。

2.3后處理及分析

將各區(qū)域模型的屬性和條件設(shè)置完成后，創(chuàng)建衍生截面用于模擬計算后的分析，并建立所需要的數(shù)據(jù)報告，以便于計算織物系統(tǒng)的熱阻、克羅值等相關(guān)數(shù)據(jù)。先進行初始化運行，然后運行迭代計算，圖5為步長為0.01s，迭代計算5000步后的溫度云圖。從監(jiān)控的殘差曲線圖上可看出，此時質(zhì)量和能量的傳遞已經(jīng)趨于穩(wěn)定，所以截取此時的溫度云圖用于可視化分析。圖5（a）（b）為織物上表面和下表面的云圖，可以看到溫度在紗線表面的分布情況，離接觸面越近的部位溫度越接近外表面邊界的溫度。圖5（c）為平行于模型XY平面，并且在Z方向上模型系統(tǒng)的1/2位置處的截面云圖，可以看到“月牙”形的溫度差異，所以在紗線交織部位，沿紗線傳遞的熱量比空氣傳遞熱量快。圖5（d）為平行于模型YZ平面，并且在X方向上模型系統(tǒng)的1/4位置處的截面云圖，可以看到沿紗線輪廓有“波浪”形的溫度變化，所以在沿Z方向傳遞的過程中，熱量在固體紗線與空氣流體交界面?zhèn)鬟f時出現(xiàn)滯后現(xiàn)象，這是因為二者的物理屬性不同。

3實驗與模擬結(jié)果對比

3.1熱舒適性實驗

測試織物導(dǎo)熱性能，采用國家標(biāo)準(zhǔn)GB/T11048—2018《紡織品生理舒適性穩(wěn)態(tài)條件下熱阻和濕阻的測定（蒸發(fā)熱板法）》，圖5溫度云圖

Fig.5Temperature cloud使用YG606G熱阻濕阻測試儀（寧波紡織機械有限公司）進行測試，測試儀器實物如圖6所示，原理如圖7所示。按照標(biāo)準(zhǔn)設(shè)置相關(guān)參數(shù)，先測試該儀器的空板熱阻Rct0，再將60cm×60cm的織物樣布放入測試板上進行測試，利用下式自動計算出織物的熱阻值Rct。

0=（Tm-Ta）·AΔH-Hc（5）

Rct=（Tm-Ta）·AH-ΔHc-Rct0（6）

式中：Rct0為儀器空板熱阻，m2·K/W;Rct為織物樣品的熱阻，m2·K/W;Tm為實驗板的溫度，℃;Ta為氣候室的溫度，℃;A為實驗板的面積，m2;H為提供給加熱面板的測試功率，W;ΔHc為熱阻測定中的加熱功率的修正量。

3.2織物系統(tǒng)熱阻計算

棉纖維紗線的導(dǎo)熱系數(shù)是已知的，但是棉紗線以不同的形態(tài)結(jié)構(gòu)織成織物后，織物系統(tǒng)的整體熱阻不是空氣和紗線熱阻的簡單疊加，并且空氣和紗線的形態(tài)均不規(guī)則。所以在計算織物系統(tǒng)熱阻時，采用的是傅里葉定律求得熱流密度，再類比電阻電流的關(guān)系推導(dǎo)出熱阻值[17]。

q=-λTx=φA0（7）

式中：q為熱流密度，W/m2;λ是導(dǎo)熱系數(shù)，W/m·℃;Tx是沿?zé)崃總鬟f方向的溫度梯度;φ為轉(zhuǎn)移總熱量，W;A0為單元模型垂直于Z方向的截面面積，m2。

織物系統(tǒng)的熱阻由式（7）推導(dǎo)而出：

R=Tb-Ttq（8）

式中：Tb為織物系統(tǒng)下表面平均溫度，℃;Tt為織物系統(tǒng)上表面平均溫度，℃。

3.3實驗值與模擬值對比

根據(jù)式（7）和式（8）可計算出單元織物系統(tǒng)的熱阻，然后織物的克羅值約為熱阻值的6.45倍，傳熱系數(shù)為熱阻值的倒數(shù)，由此可以計算得到克羅值與熱傳系數(shù)。有限元模擬計算的結(jié)果與實際測量結(jié)果的對比如圖8所示。從圖8可看出，三種不同規(guī)格的棉針織物導(dǎo)熱性能不同，紗線較粗的織物熱阻較大，但是三種均是常用汗衫面料，總體熱阻較小，方便夏季人體熱量散失、保持涼爽。通過模擬計算所得熱阻、克羅值、熱傳系數(shù)與實際測量值的誤差在4%以內(nèi)，其中21s棉紗針織物的模擬結(jié)果誤差僅在1%左右。模擬計算出的熱阻值略小于實際測量值，分析認(rèn)為是因為在進行建模時，沒有將紗線表面的毛羽考慮在內(nèi)，而在實際測量中，表面毛羽對熱阻的測量有一定的影響。但其誤差在可接受范圍內(nèi)，表明該方法具有一定的實用性。

4結(jié)論

為了建立織物幾何模型，采用超景深數(shù)碼顯微鏡和織物厚度儀測量織物結(jié)構(gòu)尺寸，利用參數(shù)數(shù)據(jù)，插入非均勻有理B樣條曲線來構(gòu)建紗線彎曲路徑，通過掃描、陣列及切割得到單元線圈模型。將線圈模型導(dǎo)入到有限元分析軟件中分析熱傳遞性能，通過溫度云圖可清晰看到熱量在織物表面的分布情況，以及分析內(nèi)部截面的溫度分布，熱量沿固體紗線傳遞比空氣要快，但在二者交界面上出現(xiàn)滯后現(xiàn)象。最終三種不同規(guī)格的緯平針織物模擬計算結(jié)果，其熱阻值、克羅值和熱傳系數(shù)與實際測量值的誤差在4%以內(nèi)，說明通過建立針織物幾何模型對織物熱傳遞進行有限元仿真可以反映、預(yù)測織物的導(dǎo)熱性能，為優(yōu)化設(shè)計織物提供理論依據(jù)，為進一步研究針織物其他仿真，如導(dǎo)濕、電磁輻射等打下基礎(chǔ)。

參考文獻：

[1]楊勇， 范君. 淺談服裝的熱濕舒適性[J]. 黑龍江紡織， 2010（3）： 1416.

YANG Yong， FAN Jun. Discussion on the heatwet comfortability of clothing[J]. Heilongjiang Textile， 2010（3）： 1416.

[2]宋蘭霞， 陳昀. 織物結(jié)構(gòu)密度對織物導(dǎo)熱性能影響的研究[J]. 中國纖檢， 2005（3）： 2225.

SONG Lanxia， CHEN Yun. Study on the influence of fabric density on thermal conductivity of fabrics [J]. China Fiber Inspection， 2005（3）： 2225.

[3]吳海軍， 錢坤， 曹海建. 織物結(jié)構(gòu)參數(shù)對熱傳遞性能的影響[J]. 紡織學(xué)報， 2007（2）： 2123.

WU Haijun， QIAN Kun， CAO Haijian. Effect of fabric structure parameters on its heat transfer property[J]. Journal of Textile Research， 2007（2）： 2123.

[4]劉讓同， 劉淑萍， 李亮， 等. 織物結(jié)構(gòu)對機織物熱傳導(dǎo)的影響[J]. 上海紡織科技， 2017， 45（6）： 47.

LIU Rangtong， LIU Shuping， LI Liang， et al. Effect of fabric structure on heat conduction of woven fabrics[J]. Shanghai Textile Science & Technology， 2017， 45（6）： 47.

[5]王婷婷， 劉曉霞， 莊明宇. 不同溫度環(huán)境下織物導(dǎo)熱機理的探討[J]. 上海紡織科技， 2015， 43（11）： 2226.

WANG Tingting， LIU Xiaoxia， ZHUANG Mingyu. Discussion on the thermal transfer mechanism of fabrics at different temperatures [J]. Shanghai Textile Science & Technology， 2015， 43（11）： 2226.

[6]吳佳佳， 唐虹. 應(yīng)用ABAQUS的織物熱傳遞有限元分析[J]. 紡織學(xué)報， 2016， 37（9）： 3741.

WU Jiajia， TANG Hong. ABAQUS based finite element analysis of heat transfer through woven fabrics [J] Journal of Textile Research， 2016， 37（9）： 3741.

[7]弋夢夢. 機織物熱傳導(dǎo)性能及其數(shù)值模擬[D]. 鄭州： 中原工學(xué)院， 2018.

YI Mengmeng. Thermal Conductivity of Woven Fabric and Its Numerical Simulation [D]. Zhengzhou： Zhongyuan Institute of Technology， 2018.

[8]GIBSON P. Modeling heat and mass transfer from fabriccovered cylinders [J]. Journal of Engineered Fabrics & Fibers， 2009， 4（1）： 18.

[9]趙超， 楊彩云. 彎紗深度與針織物結(jié)構(gòu)參數(shù)及導(dǎo)熱系數(shù)間的關(guān)系[J]. 國際紡織導(dǎo)報， 2015， 43（11）： 2224.

ZHAO Chao， YANG Caiyun. Relationship among the sinking depth and structure parameters and thermal conductive performance of knitted fabric[J]. Melliand China， 2015， 43（11）： 2224.

[10]孫玉釵. 針織物熱舒適性能研究與針織保暖產(chǎn)品設(shè)計[D]. 上海： 東華大學(xué)， 2005.

SUN Yuchai. Study on Thermal Comfort Properties of Knitted Fabrics and Design of Heat Preservation Products for Knitting [D]. Shanghai： Donghua University， 2005.

[11]RAJAB K， PIEGL L A， SMARODZINAVA V. CAD model repair using knowledgeguided NURBS [J]. Engineering with Computers， 2013， 29（4）： 477486.

[12]施法中. 計算機輔助幾何設(shè)計與非均勻有理 B 樣條[M]. 北京： 高等教育出版社， 2001： 47.

SHI Fazhong. ComputerAided Geometry Design and NonUniform Rational BSpline [M]. Beijing： Higher Education Press， 2001： 47.

[13]汝欣， 彭來湖， 呂明來， 等. 緯編針織物幾何建模及其算法[J]. 紡織學(xué)報， 2018， 39（9）： 4449.

RU Xin， PENG Laihu， L Minglai， et al. Modelling and algorithm of weft knitted fabric[J]. Journal of Textile Research， 2018， 39（9）： 4449.

[14]代文杰， 邱華， 楊恩惠， 等. 基于ANSYS CFX的織物透氣性數(shù)值計算[J]. 絲綢， 2018， 55（9）： 5156.

DAI Wenjie， QIU Hua， YANG Enhui， et al. Numerical calculation of fabric air permeability based on ANSYS CFX[J]. Journal of Silk， 2018， 55（9）： 5156.

[15]胡海霞， 孟家光. 織物熱傳導(dǎo)機理的分析和探討[J]. 北京紡織， 2005（6）： 2831.

HU Haixia， MENG Jiaguang. The analysis and discussion of heat transmission mechanism of the fabrics[J]. Beijing Textile Journal， 2005（6）： 2831.

[16]李明， 劉楠. STARCCM+11.0與流場計算[M]. 北京： 機械工業(yè)出版社， 2017： 78.

LI Ming， LIU Nan. STARCCM+ 11.0 and Flow Field Calculation [M]. Beijing： Mechanical Industry Press， 2017： 78.

[17]李瑛慧， 謝春萍， 劉新金， 等. 仿真絲織物與真絲織物的熱傳遞有限元仿真[J]. 絲綢， 2017， 54（12）： 711.

LI Yinghui， XIE Chunping， LIU Xinjin， et al. Finite element simulation of heat transfer on silklike fabrics and real silk fabrics [J].Journal of Silk， 2017， 54（12）： 711.

收稿日期： 20190403; 修回日期： 20191129

基金項目： 江蘇省產(chǎn)學(xué)研項目（BY201602214）

作者簡介： 楊恩惠（1994），女，碩士研究生，研究方向為針織物熱濕舒適性的仿真研究。通信作者：邱華，教授，qiuhua@jiangnan.edu.cn。