采用ABAQUS的粘膠機織物拉伸力學(xué)性能仿真

劉倩楠， 劉新金

(生態(tài)紡織教育部重點實驗室(江南大學(xué))， 江蘇 無錫 214122)

拉伸破壞作為機織物最基本的力學(xué)性能，在紡織品的設(shè)計和生產(chǎn)中有著重要的參考價值。近年來，許多研究人員采用了有限元理論分析法[1]研究機織物拉伸力學(xué)性能。曹榮平[2]利用ABAQUS有限元軟件模擬了機織物建筑膜材料單軸向與雙軸向的拉伸，討論了材料屬性、Step值、紗線相互作用等因素對機織物建筑膜材料拉伸模擬結(jié)果的影響；程建芳等[3]利用有限元軟件ANSYS求解Kevlar 129紗線及織物帶有主應(yīng)力斷裂失效準(zhǔn)則的彈性正交各向異性材料模型，分析了織物的拉伸力學(xué)性質(zhì)。這類研究均是基于紗線的材料屬性和織物的結(jié)構(gòu)特征建立機織物三維細(xì)觀模型，利用有限元軟件創(chuàng)建拉伸環(huán)境模擬織物拉伸過程，求解數(shù)學(xué)模型。

由于機織物中經(jīng)緯紗的擠壓變形，單紗材料的非線性和織物拉伸時影響因素的多樣性，導(dǎo)致機織物拉伸模型求解存在一定的難度。隨著專業(yè)紡織建模軟件的開發(fā)和有限元軟件的應(yīng)用，突破了織物三維建模的難度和拉伸環(huán)境設(shè)置的局限性。本文以粘膠平紋織物為例，基于織物結(jié)構(gòu)參數(shù)，借助專業(yè)紡織建模軟件Texgen建立織物模型，利用有限元軟件ABAQUS模擬織物拉伸力學(xué)性能，并通過實驗驗證模擬數(shù)值結(jié)果。

1 實驗部分

1.1 粘膠單紗強力測試

從織物中取出帶有屈曲的單根紗線，參照GB/T 3916—1997 《紡織品 卷裝紗 單根紗線斷裂強力和斷裂伸長率的測定》[4]，使用YG020B型電子單紗強力機，設(shè)置有效加持距離為200 mm，拉伸速度為200 mm/min，分別對織造粘膠織物的經(jīng)緯單紗重復(fù)測試10次，結(jié)果取平均值，得到粘膠單紗拉伸強力參數(shù)，如表1所示。

表1 單紗拉伸強力參數(shù)Tab.1 Single yarn tensile strength parameters

1.2 粘膠織物拉伸性能測試

粘膠平紋織物規(guī)格如表2所示，按照GB/T 3923.1—2013 《紡織品 織物拉伸性能 第1部分：斷裂強力和斷裂伸長率的測定(條樣法)》說明，用條樣法將織物制作成規(guī)格為35 cm×5 cm的試樣，使用YG026D型多功能電子織物強力機[5]，設(shè)置隔距200 mm，拉伸速度100 mm/min，重復(fù)測試5次。

表2 織物規(guī)格參數(shù)Tab.2 Fabric specification parameters

2 有限元模擬

2.1 Texgen模型建立

Texgen織物仿真軟件可準(zhǔn)確地模擬紗線或織物的幾何結(jié)構(gòu)[6]，其良好的圖形用戶界面(GUI)使用戶可以在“Weave Wizard”中設(shè)置經(jīng)緯紗數(shù)量、紗線間隔、紗線寬度、紗線高度、織物厚度和經(jīng)緯紗交織形式等一系列指標(biāo)，完成模型的初步建立；在“Modeller”中設(shè)置紗線截面形狀、屈曲形態(tài)和紗線屬性等參數(shù)進一步細(xì)化模型。

本文將表3所示粘膠平紋織物的幾何結(jié)構(gòu)參數(shù)輸入Texgen軟件，得到如圖1、2所示的織物經(jīng)向和緯向拉伸模型。再將建立的模型以.stp格式輸入有限元軟件ABAQUS中進行拉伸力學(xué)性能分析[7]。

表3 織物幾何結(jié)構(gòu)參數(shù)Tab.3 Fabric geometry parameters mm

注：圖中數(shù)值單位為mm。

注：圖中數(shù)值單位為mm。

2.2 ABAQUS有限元模擬

ABAQUS是一套功能相對強大的有限元分析軟件，既能分析簡單的線性靜力學(xué)問題，又能解決工程中復(fù)雜模型的高度非線性問題[8]。ABAQUS有限元分析時有ABAQUS/CAE前處理、分析計算、ABAQUS/Viewer后處理3個階段，其中分析計算階段有ABAQUS/Standard和ABAQUS/Explicit 2種分析形式，ABAQUS/Standard可分析廣泛領(lǐng)域的線性和非線性問題，但對于短暫、瞬時的動態(tài)事件的模擬，以及沖擊和高度不連續(xù)問題的求解不及ABAQUS/Explicit。粘膠平紋織物的拉伸實驗中實驗材料屬于非線性材料，且織物交織過程中經(jīng)緯紗之間存在一定的摩擦，拉伸時經(jīng)緯紗線間產(chǎn)生滑移，所以更適合使用ABAQUS/Explicit顯示分析。

在ABAQUS/CAE中導(dǎo)入Texgen創(chuàng)建的織物模型，Property模塊將紗線定義為彈塑性材料，即應(yīng)力低于屈服點時，表現(xiàn)為彈性行為，應(yīng)力大于屈服點后，發(fā)生塑性變形。參照表1數(shù)據(jù)，賦予經(jīng)緯紗這種非線性材料屬性[9]。經(jīng)緯紗非線性材料屬性基于紗線體積質(zhì)量，通過彈性模量和泊松比定義紗線彈性行為，塑性應(yīng)力和應(yīng)變定義紗線的塑性變形。

分析計算階段采用動態(tài)顯示分析，用Step模塊創(chuàng)建初始分析步后，選擇“Dynamic Explicit”為“Step1”，設(shè)置“Time Period”為1E-5，然后設(shè)置相應(yīng)的場輸出和歷程輸出，并將輸出頻率設(shè)置為200；用Interaction模塊定義織物中經(jīng)緯紗的接觸，在“Step1”下設(shè)置“General contact(Explicit)”，定義紗線與紗線間切向摩擦因數(shù)為0.17；用Load模塊設(shè)置1個“PINNED(U1=U2=U3=0)”的邊界條件，將織物一端固定，并且給另一端定義1個100 mm/min的拉伸速度；用Mesh模塊分別對經(jīng)緯紗進行網(wǎng)格劃分，設(shè)置網(wǎng)格種子尺寸為0.1，由于紗線截面近似橢圓形，采用中性軸算法用六面體對織物模型進行網(wǎng)格劃分，將其劃分為若干個單元，圖3示出粘膠平紋織物經(jīng)向拉伸試樣網(wǎng)格模型。其中圖3(b)為3(a)中黑色方格所表示的單元；在Job模塊給要分析的模型數(shù)據(jù)創(chuàng)建1個作業(yè)，進行提交計算。

圖3 織物網(wǎng)格模型Fig.3 Mesh model of fabric. (a) Mesh model of stretch in warp direction; (b) Element shown in black square of mesh model

在ABAQUS/Viewer后處理階段顯示ODB文件中的分析計算結(jié)果，包括變形前/后的模型圖、矢量/張量符號圖、材料方向圖、各種變量的分布云圖、變量的X-Y圖表、動畫等[10]。本文分別模擬了粘膠平紋織物的經(jīng)向拉伸和緯向拉伸，其應(yīng)力分布情況如圖4所示(圖中顏色越深，表示應(yīng)力越大)。

從圖4可看出，織物模型拉伸后發(fā)生變形，此時經(jīng)緯紗上均出現(xiàn)應(yīng)力分布。經(jīng)向拉伸時，經(jīng)紗上分布的應(yīng)力較緯紗明顯許多，這是由于織物經(jīng)向拉伸過程中，經(jīng)紗從屈曲狀態(tài)被拉直，然后再被拉斷，起到了主要的拉伸作用，而緯紗由于拉伸過程中紗線之間產(chǎn)生滑移、摩擦等作用(Interaction模塊賦予了紗線間0.17的摩擦因素)，因此在拉伸過程中，通過經(jīng)緯交織點的應(yīng)力傳遞，產(chǎn)生一定的應(yīng)力分布；緯向拉伸時，緯紗承受主要的拉伸作用，緯紗上分布的應(yīng)力較明顯，經(jīng)紗由于拉伸過程中紗線之間產(chǎn)生滑移、摩擦等作用，經(jīng)經(jīng)緯交織點的應(yīng)力傳遞，產(chǎn)生一定的應(yīng)力分布。圖4(b)、(d)為通過ABAQUS的“Remove Selected”選項，移除其余經(jīng)緯紗，僅剩1根經(jīng)紗和1根緯紗，從圖上可清晰地看出織物拉伸過程中經(jīng)緯紗的受力情況。經(jīng)向拉伸時，經(jīng)紗上紗線交織點處所受應(yīng)力最大，為23.78 MPa，緯紗上緯紗邊緣所受應(yīng)力最小，為6.266×10-3MPa；緯向拉伸時，緯紗上紗線交織點處所受應(yīng)力最大，為15.27 MPa，經(jīng)紗上經(jīng)紗邊緣所受應(yīng)力最小，為4.399×10-3MPa。

3 結(jié)果與分析

注：圖中1、2、3、4、5表示節(jié)點。

圖5 織物拉伸力學(xué)性能對比Fig.5 Comparison of tensile properties of fabrics. (a) Fabric tension-elongation curve of stretch in warp direction; (b) Finite element simulation of tensile stress-strain curve of stretch in warp direction; (c) Fabric tension-elongation curve of stretch in weft direction; (d) Finite element simulation of tensile stress-strain curve of stretch in weft direction

采用YG026D型多功能電子織物強力機對粘膠平紋織物進行5次測試，分別得到織物經(jīng)向拉伸和緯向拉伸的平均拉伸斷裂強力、拉伸長度、伸長率等參數(shù)如表4所示。將其轉(zhuǎn)化為拉力-伸長曲線，結(jié)果如圖5所示。

織物在被拉伸的過程中屈曲的紗線首先被拉直，然后進入彈性變形階段，在這個階段紗線所受拉伸強力呈線性增長，當(dāng)達到最大斷裂強力時，斷裂強度最低的紗線開始斷裂，接著大部分紗線斷裂，部分紗線滑移致使織物被拉斷。

實際中紗線由纖維組成，織物拉伸時組成紗線的纖維間會產(chǎn)生滑移、剪切、摩擦等一系列的作用，織物的斷裂實際上是織物中紗線的斷裂，而紗線的斷裂則是紗線中纖維的滑移和斷裂。圖5(b)、(d)所示為有限元模擬織物拉伸的應(yīng)力-應(yīng)變曲線圖(所選節(jié)點分別為織物變形后經(jīng)緯紗交織點處經(jīng)紗和緯紗上節(jié)點，見圖4(a)、(c))，有限元模擬時，將紗線看做一個整體，賦予其非線性的材料屬性，使其產(chǎn)生理論化的變形，但是忽視了紗線中纖維斷裂的隨機性，達到最大應(yīng)力時開始產(chǎn)生脆性斷裂。

表4 織物拉伸強力參數(shù)Tab.4 Fabric tensile strength parameters

對比圖5(a)與圖5(b)～(d)，發(fā)現(xiàn)2類曲線的上升趨勢大致相同，并且均上升至某最大值；根據(jù)織物拉伸強力測試得到粘膠平紋織物經(jīng)向拉伸強力317.68 N、緯向拉伸強力191.87 N，經(jīng)公式σ=F/A(式中：σ為拉伸應(yīng)力，MPa；F為拉伸強力，N；A為拉伸切線方向橫截面積，mm2)轉(zhuǎn)換得到經(jīng)向拉伸應(yīng)力為24.82 MPa，緯向拉伸應(yīng)力為14.99 MPa，利用相對誤差計算公式可得到模擬與實驗的相對誤差分別為4.19%、1.87%，經(jīng)均值計算，平均誤差為3.03%，說明有限元軟件ABAQUS可用來模擬粘膠機織物拉伸力學(xué)性能，且模擬數(shù)值具有一定可靠性。

4 結(jié)束語

利用計算機軟件對織物拉伸性能進行模擬，可以使設(shè)計人員預(yù)測織物性能時，減少試紡、試織、后處理等工序，達到節(jié)約資源，降低成本，縮短周期的目的。本文借助專業(yè)紡織建模軟件Texgen建立織物幾何模型，利用有限元軟件ABAQUS模擬粘膠織物拉伸力學(xué)性能，通過所得應(yīng)力分布云圖，發(fā)現(xiàn)了織物拉伸時應(yīng)力分布規(guī)律：織物經(jīng)向拉伸時，主要是經(jīng)紗承受拉伸作用，經(jīng)緯紗交織點處所受應(yīng)力最大，緯紗受紗線間摩擦因素影響，緯紗邊緣所受應(yīng)力最小；織物緯向拉伸時，主要是緯紗承受拉伸作用，經(jīng)緯紗交織點處所受應(yīng)力最大，經(jīng)紗受紗線間摩擦因素影響，經(jīng)紗邊緣所受應(yīng)力最小。最后，通過實驗驗證了計算機模擬數(shù)值結(jié)果的可靠性。由于本文僅以粘膠織物為例，對于使用該方法預(yù)測織物性能是否具有普適性，還有待進一步論證。