三維六角形編織結(jié)構(gòu)的計算機模擬

高彥濤，KO Frank，胡 紅，3

（1.東華大學 紡織學院，上海 201620；2.英屬哥倫比亞大學 材料系 ＆先進材料工程工藝實驗室，溫哥華 V6T1Z4；3.香港理工大學 紡織與服裝系，香港）

三維六角形編織結(jié)構(gòu)的計算機模擬

高彥濤1，KO Frank2，胡 紅1，3

（1.東華大學 紡織學院，上海 201620；2.英屬哥倫比亞大學 材料系 ＆先進材料工程工藝實驗室，溫哥華 V6T1Z4；3.香港理工大學 紡織與服裝系，香港）

基于三維六角形編織原理，采用編織紗線位置點的矩陣表示法，根據(jù)紗線的運動情況，得到每一步運動后的紗線位置坐標，然后采用B－spline曲線對每一根紗線進行擬合，用管狀實體代表紗線對整個編織體進行三維顯示，經(jīng)過優(yōu)化后得到編織體的仿真圖像.仿真結(jié)果表明，該模擬方法能夠?qū)碗s六角形編織結(jié)構(gòu)進行形象模擬，為探索編織體中內(nèi)部紗線排布提供有力工具.

六角形編織；B－spline曲線；模擬仿真；Matlab編程

新型三維六角形編織是一種旋轉(zhuǎn)編織方法，它能夠編織多種復雜形狀結(jié)構(gòu)件，以及利用幾十微米紗線進行微觀精細化編織［1－4］，其產(chǎn)品在生物醫(yī)用以及復合材料領域具有較大應用潛力.這種六角形編織與傳統(tǒng)的縱橫編織原理不同，因此，縱橫編織的模擬方法不適用于六角形編織.文獻［5－7］詳細闡述了基于傳統(tǒng)縱橫編織原理得到的編織體的模擬方法，為了實現(xiàn)三維模擬仿真，它們都聯(lián)合多種編程軟件，通過一些文件修改與傳遞進行反復模擬，這給模擬的靈活性帶來一定問題.

基于數(shù)學算法并利用Matlab編程，本文建立了一種針對六角形編織的模擬方法.一些新算法很容易加入模擬程序，易于進一步優(yōu)化.同時，該方法可以進一步擴展到其他編織方法的模擬仿真，也適用于單胞模型的建立，對編織結(jié)構(gòu)的研究有一定意義.

1 三維六角形編織的基本原理

旋轉(zhuǎn)編織是通過控制凸輪旋轉(zhuǎn)以及卷取機理來獲得編織結(jié)構(gòu)的方法.三維六角形編織是由Frank Ko教授發(fā)明的一種新型旋轉(zhuǎn)編織［2］.事實上，控制凸輪的運動是為了移動載紗器使它在軌道板上運動來使紗線纏繞.為了得到需要的編織結(jié)構(gòu)，在進行編織之前，必須確定攜紗器（編織紗）的初始位置、凸輪的運動和卷取速度.

圖1通過兩個相鄰六角形凸輪的3步連續(xù)運動描述了六角形編織基本原理.圖1中3個不同標號的圓圈分別代表進行編織的3根紗線，凸輪的旋轉(zhuǎn)模式由凸輪軸邊上的旋轉(zhuǎn)箭頭表示.其中，圖1（a）代表編織紗的初始位置，圖1（b）～1（d）分別表示凸輪步進一次后編織紗的位置.從圖1（b）可以看出，經(jīng)過左邊六角形凸輪逆時針旋轉(zhuǎn)后，紗線1和2分別按照逆時針移動了一個位置.依次類推，可得到紗線在底盤的運動規(guī)律.

圖1 三維六角形編織工藝示意圖Fig.1 Schematic of 3Dhexagonal braiding process

2 三維六角形編織命令文件與編織紗的矩陣表示

為獲得需要的編織結(jié)構(gòu)，必須根據(jù)編織結(jié)構(gòu)的要求確定電機運轉(zhuǎn).電機的運轉(zhuǎn)模式保存于一個.txt文件，即指令文件，該文件包含編織機上所有凸輪每一步運動信息.指令文件包括兩部分：一是機器卷取速度；二是電機運動模式.各個電機運動模式都是用數(shù)字表示，“1”代表順時針旋轉(zhuǎn)60°，“－1”代表逆時針旋轉(zhuǎn)60°，“0”代表靜止.一個用于編織管狀結(jié)構(gòu)的電機指令文件實例如圖2所示.

圖2中左側(cè)第1和2列代表電機的定位，例如，motor 1－2代表第1圈第2個電機（凸輪）.從第3列開始的數(shù)字分別代表對應電機的運動模式，例如，motor 1－2后的3個數(shù)字分別為－1，－1，0，這表示motor 1－2在前兩步編織時，逆時針旋轉(zhuǎn)60°，第3步靜止不動，以此類推.假如某電機第n步運動狀態(tài)為1，代表電機在第n步執(zhí)行順時針旋轉(zhuǎn)60°.利用此方法就可以表達所有編織過程.

圖2 電機指令文件Fig.2 Command file of motors

實現(xiàn)編織結(jié)構(gòu)的三維模擬需要從電機指令文件出發(fā)，依據(jù)電機指令文件獲取載紗器的初始位置，并將其存儲起來，存儲形式為一個兩維矩陣.圖3描述了編織紗線的矩陣表示形式.

圖3 編織紗的矩陣表示示意圖Fig.3 Schematic of matrix representation of braiding yarn

依據(jù)圖3中矩陣表示方法，在模擬編織過程中需要計算矩陣的縱行和橫列數(shù)目，這兩個參數(shù)在程序中分別用“nb＿row”和“nb＿column”表示，其值可由下列方程式給出：

其中：nb＿level表示編織機器中凸輪的層數(shù).

3 三維六角形編織的模擬方法

編織結(jié)構(gòu)的模擬由編織工藝決定，三維六角形編織的主要工藝參數(shù)包括凸輪運動模式、編織尺寸、面板上攜紗器布置形狀等.為了正確地模擬編織結(jié)構(gòu)首先必須得到這些工藝參數(shù)，其模擬的流程圖如圖4所示.

圖4 三維六角形編織模擬流程圖Fig.4 Flow chart of hexagonal braiding simulation

由圖4可以看出，實現(xiàn)三維六角形編織模擬的主要過程包括4部分：編織工藝參數(shù)的讀取；紗線位置的計算；創(chuàng)建B－spline曲線；模擬優(yōu)化.

3.1 初始數(shù)據(jù)的讀取

新型六角形編織是一種三維旋轉(zhuǎn)編織，具有編織復雜形狀的能力.為創(chuàng)造復雜結(jié)構(gòu)，指令文件必須分別對每個電機獨立控制.在進行編織模擬時，先讀取指令文件.基于電機指令文件格式，編寫相應程序來讀取一些基本的編織參數(shù)，程序流程如圖5所示.

3.2 計算紗線各步的位置

圖5 數(shù)據(jù)讀取流程圖Fig.5 Flow chart of data acquisition

在考慮編織運動之前，必須先定義攜紗器的初始位置.第一根紗線用數(shù)字“1”來表示，紗線總數(shù)為Nc.基于矩陣表示方法，可通過編寫程序來實現(xiàn)各個紗線位置點分布.攜紗器的初始位置確定以后，根據(jù)編織機每一步的運動情況就可以確定攜紗器在每一步的位置.載紗器位置計算流程如圖6所示.由圖6可以看出，對于編織循環(huán)的每一步而言，必須考慮所有的凸輪轉(zhuǎn)動.通過研究凸輪轉(zhuǎn)動，建立起每步運動后攜紗器在相應矩陣中的位置點.

圖6 攜紗器位置計算流程圖Fig.6 Flow chart of carrier position calculation

3.3 B－spline曲線的生成

圖6中獲得編織紗線矩陣給出了所有紗線在不同編織步的位置，基于這些數(shù)據(jù)可以畫出每一根紗線沿編織方向延伸的位置.在模擬過程中，引入?yún)?shù)步長hstep，它表示兩步編織之間的編織長度.3根紗線在空間的延伸路徑模擬如圖7所示.

圖7 編織體中紗線路徑Fig.7 Yarn paths in braid

由圖7可知，編織紗線位置點延伸是由線段組成，它們不能用來代表紗線.為了使編織紗線在編織結(jié)構(gòu)中平滑延伸，可利用B－spline曲線來擬合數(shù)據(jù)點，所得到的擬合曲線能夠合理地代表紗線.為了定義B－spline曲線，必須確定3個參數(shù)（節(jié)點ti、控制點數(shù)f、曲線次數(shù)n）.依據(jù)這些參數(shù)，B－spline曲線可做如下描述.

首先，所有確定節(jié)點按如下方式排列：

同時定義一個整數(shù)：

然后，B樣條曲線可以寫成如下形式：

其中：Pi為控制點位置向量；

通過編寫B(tài)－spline曲線模擬程序即可實現(xiàn)圖7中對應于矩陣控制點的三階樣條曲線，如圖8所示.

為了更加直觀形象地描述紗線，由一個自定義函數(shù)tubeplot來實現(xiàn)紗線的管狀實體顯示.圖8中的3根紗線利用此功能進行的管狀實體表達形式如圖9所示.

圖8 編織紗線路徑B－spline曲線擬合Fig.8 B－spline curve fitting of braiding yarn paths

圖9 紗線的管狀實體顯示Fig.9 Description of braiding yarn using solid tube

3.4 優(yōu)化過程

通過以上步驟，可以看出紗線編織后的纏繞方式，但是由于紗線之間沒有接觸，這不能代表紗線的真實位置.圖10給出了一個優(yōu)化前的編織體模擬圖實例，可以看出紗線在編織體內(nèi)沒有相互接觸.事實上，當進行編織時，由于紗線張力作用，紗線在編織體中趨向于走最短路徑從而貫穿于整個編織體長度內(nèi)，紗線會彼此接觸緊挨在一起.因此，為了更好地模擬編織體，必須實施編織結(jié)構(gòu)的模擬優(yōu)化功能.

圖10 優(yōu)化前的編織體模擬圖Fig.10 Simulation of a braid before optimization

通過改變各根紗線之間整體位置來實現(xiàn)編織結(jié)構(gòu)的模擬優(yōu)化，從第一根紗線開始，循環(huán)檢查紗線間的距離，如果不相互接觸，對B－spline曲線控制坐標乘一個系數(shù)來使紗線往編織結(jié)構(gòu)中心靠近，直到紗線之間有接觸.實現(xiàn)此功能的源代碼如下：

在紗線坐標移動過程中，紗線的形狀保持不變.優(yōu)化前后的模擬編織結(jié)構(gòu)對比如圖11所示.由圖11可以看出，經(jīng)過優(yōu)化之后的模擬圖更能反映編織體的真實結(jié)構(gòu).

圖11 優(yōu)化前后編織結(jié)構(gòu)模擬圖對比Fig.11 Comparison of braid simulation before and after optimized process

不同結(jié)構(gòu)形狀的編織體模擬圖如圖12所示.由此可以看出本文方法能夠較好地對三維六角形編織結(jié)構(gòu)體進行模擬仿真，為研究編織體的內(nèi)部結(jié)構(gòu)提供了有力工具.

圖12 不同形狀的編織體模擬圖Fig.12 Simulation of braided structures with different cross－sections

4 結(jié) 語

本文在分析三維六角形編織原理及編織過程的基礎上，提出了編織過程中編織紗線矩陣表示方法.基于編織紗線的空間點矩陣，通過B－spline曲線擬合，實現(xiàn)了對編織體中紗線運動規(guī)律的模擬，并系統(tǒng)地介紹了這種基于Matlab編程的編織結(jié)構(gòu)模擬方法.模擬得到的織物能夠真實地反映了紗線在編織結(jié)構(gòu)的運動規(guī)律，使得編織體內(nèi)部結(jié)構(gòu)可視化.

由于該模擬方法是利用基于數(shù)學方法的Matlab編程來實現(xiàn)，一些自定義新算法也可以施加在此模擬方法上，因此，該方法可以很容易擴展到其他編織方法的模擬仿真以及單胞模型的建立.

參 考 文 獻

［1］BOGDANOVICH A E，MUNGOLOV D.Recent advancements in manufacturing 3－D braided preforms and composites［C］／／Proceedings of the Acun－4International Composites Conference.Sydney，Australia，2002：61－72.

［2］KO K F，SCHREIBER F，YANG H J，et al.Recent advancements in three dimensional braiding［C］／／Proceedings of the 1st World Conference on 3D－Fabrics and Their Applications.Manchester，UK，2008：10－11.

［3］SCHREIBER F，KO K F，YANG H J，et al.Novel three dimensional braiding approach and its application［C］／／17th International Conference on Composite Materials.Edinburgh，London：IOM Communications，2009：27－31.

［4］成鈴.三維編織的計算機模擬技術［J］.天津工業(yè)大學學報，2001，21（2）：26－29.

［5］張美忠，李賀軍，李克智.三維編織復合材料矩形預制件仿真［J］.航空學報，2006，27（5）：985－988.

［6］王毅強，張立同，成來飛.三維編織體復合材料空間幾何結(jié)構(gòu)的計算機模擬［J］.航空材料學報，2008，28（2）：95－98.

［7］邵將，溫衛(wèi)東，崔海濤.三維四步法編織復合材料結(jié)構(gòu)的計算機仿真［J］.南京航空航天大學學報，2009，41（1）：36－40.

Computer Simulation of 3DHexagonal Braided Structures

GAOYan－tao1，KOFrank2，HUHong1，3
（1.College of Textiles，Donghua University，Shanghai 201620，China；2.Department of Material ＆AMPEL，University of British Columbia，Vancouver，BC V6T1Z4，Canada；3.Institute of Textiles and Clothing，The Hong Kong Polytechnic University，Hong Kong，China）

According to the principle of hexagonal braiding，yarns'positions at every step were obtained using matrix representation based on carrier movements.B－spline curve of each yarn was fitted through their spatial positions and yarns in braid were described as solid tubes.Finally，the 3D braids could be vividly simulated after optimization process.The simulation results show that the method is suitable for complicated structures，and provides a useful tool for investigating yarn's arrangements in braids.

hexagonal braiding；B－spline curve；simulation；Matlab programme

TB 332

A

1671－0444（2013）06－0785－05

2012－10－31

高彥濤（1984—），男，河南商丘人，博士，研究方向為紡織結(jié)構(gòu)復合材料.E－mail：dhytgao＠hotmail.com

胡 紅（聯(lián)系人），男，教授，E－mail：huhong－dhu＠hotmail.com