輸棉通道位置對轉(zhuǎn)杯紡纖維運動的影響

鄧茜茜 楊瑞華

摘要： 為了研究輸棉通道位置對纖維在紡紗通道中運動的影響，文章以直徑42 mm轉(zhuǎn)杯為例，采用Ansys軟件中ICEM-CFD建立三維抽氣式轉(zhuǎn)杯紡紗通道模型，在Fluent商用軟件中計算連續(xù)相氣體流場，根據(jù)連續(xù)相計算結(jié)果擬合離散相纖維運動軌跡。保持邊界條件不變，比較輸棉通道位置不同時（Y=13 mm、Y=9.75 mm、Y=6.5 mm）轉(zhuǎn)杯紡紗通道內(nèi)氣流分布特征及纖維運動軌跡的差異。試驗結(jié)果表明：Y=13 mm（工況a）時，輸棉通道內(nèi)氣流速度與壓強變化梯度最大，有利于伸直纖維;工況b時，轉(zhuǎn)杯凝聚槽負壓最大，纖維進入轉(zhuǎn)杯能較好地凝聚匯集，有利于纖維緊密凝聚成紗;輸棉通道位置過低（工況c），轉(zhuǎn)杯內(nèi)部渦流分布不均，不利于轉(zhuǎn)杯高速平穩(wěn)轉(zhuǎn)動。

關(guān)鍵詞： 轉(zhuǎn)杯紡;Fluent;輸棉通道;氣流;數(shù)值模擬

Abstract： In order to study the influence of the fiber transfer channel position on the motion of the fiber in the spinning channel， a 42 mm-diameter rotor was used as an example to establish a three-dimensional pumping rotor spinning channel model with ICEM-CFD in Ansys software. Continuous phase gas flow field was calculated in the Fluent commercial software. The discrete phase fiber motion trajectory was fitted according to the continuous phase calculation result. The boundary conditions remained unchanged， and the airflow distribution characteristics and differences of fiber motion trajectory in the rotor spinning channel were compared when the position of the fiber transport channel was different（Y=13 mm， Y=9.75 mm， Y=6.5 mm）. The test results showed that： Y=13 mm（condition a）， the airflow velocity and pressure gradient in the fiber transport channel were the largest， which is conducive to straightening the fiber. Under working condition b， the negative pressure of the rotor was the largest， and the fiber in the rotor could be well converged， which is beneficial to the close cohesion of the fibers into the yarn. When the position of the fiber transport channel was too low（condition c）， and the vortex distribution inside the rotor was uneven， which is not conducive to the high-speed and smooth rotation of the rotor.

Key words： rotor spinning; Fluent; fiber transfer channel; airflow; numerical simulation

轉(zhuǎn)杯紡作為新型紡紗的一種，具有高速度、大卷裝、廢落棉適紡性良好等特點。轉(zhuǎn)杯紡紗時，由分梳輥對纖維須條進行開松梳理，纖維隨紡紗器內(nèi)負壓形成的氣流經(jīng)過輸棉通道進入轉(zhuǎn)杯，在凝聚槽內(nèi)匯集并合形成纖維條，經(jīng)過加捻形成紗線從引紗管引出[1]。近年來，由于計算機技術(shù)的快速發(fā)展，基于計算流體動力學(xué)的數(shù)值計算方法廣泛應(yīng)用于紡織類領(lǐng)域，國內(nèi)外學(xué)者通過數(shù)值計算對紡織加工技術(shù)進行研究，以優(yōu)化改進紡紗方法。轉(zhuǎn)杯紡在封閉空間內(nèi)的高速強負壓氣流場中完成成紗過程，研究成紗過程中的氣流場特征和纖維運動形態(tài)是分析紗線結(jié)構(gòu)性能和改進成紗過程的重要手段。Kong等[2]通過數(shù)值模擬方法探索了分梳輥位置對轉(zhuǎn)杯紡氣流場的影響，分析了輸棉通道中的纖維形態(tài);武傳宇等[3]采用Fluent軟件對轉(zhuǎn)杯內(nèi)氣流場進行分析，指出轉(zhuǎn)杯滑移角小于27°最有利于纖維凝聚成紗;肖美娜等[4]研究了轉(zhuǎn)杯紡紗通道內(nèi)氣流的運動，分析了轉(zhuǎn)杯的轉(zhuǎn)速、轉(zhuǎn)動方向等因素對紡紗通道內(nèi)氣流流動特性的影響;林惠婷等[5]分析了輸棉通道內(nèi)氣流流動特征及彎鉤纖維運動形態(tài)，指出輸棉通道內(nèi)纖維傾向于向中心處運動。國內(nèi)外學(xué)者們研究了工藝參數(shù)對紡紗通道氣流場的影響，而關(guān)于纖維在紡紗通道內(nèi)運動規(guī)律的相關(guān)文章報道較少。

本文對轉(zhuǎn)杯紡紗通道中纖維的運動特性進行分析，應(yīng)用Fluent軟件對不同輸棉通道位置時轉(zhuǎn)杯紡紗通道內(nèi)的氣流分布進行模擬計算;采用DPM模型，將氣流視為連續(xù)相，纖維簡化為離散顆粒，探索轉(zhuǎn)杯紡紗過程中纖維分布及運動規(guī)律，為優(yōu)化轉(zhuǎn)杯紡紗器結(jié)構(gòu)提供理論依據(jù)。

1 成紗器三維幾何模型

研究思路為先建立轉(zhuǎn)杯紡成紗器的動力學(xué)模型，計算氣體連續(xù)相流場;再結(jié)合流場變量計算每顆粒子的受力情況，獲得粒子的速度并追蹤每一顆粒子的軌道，研究纖維在轉(zhuǎn)杯紡紗器中的運動軌跡。

在轉(zhuǎn)杯紡紗過程中，抽氣通道將轉(zhuǎn)杯內(nèi)空氣抽走，轉(zhuǎn)杯內(nèi)形成負壓，輸棉通道與引紗管作為補氣口，輸棉通道入口進入的空氣將纖維有序地輸送至轉(zhuǎn)杯。在此背景下，建立抽氣式轉(zhuǎn)杯紡紗機的紡紗通道模型，研究抽氣式轉(zhuǎn)杯紡紗通道的氣流分布特性。根據(jù)Fluent軟件特性，將轉(zhuǎn)杯紡紗器簡化，除去接觸不到氣流的壁面，在三維坐標系下建立用于仿真模擬的三維幾何模型，如圖1所示。其中轉(zhuǎn)杯直徑D=42 mm，滑移角ɑ=22°，坐標原點位于轉(zhuǎn)杯底部中心，在X=7 mm處建立Y-Z方向的輸棉通道豎直截面A-A，在Y=3.5 mm處建立X-Z方向的凝聚槽截面B-B。

輸棉通道中心線與轉(zhuǎn)杯壁面相交，設(shè)交點為F點，F(xiàn)點的縱坐標Y值表示F點與轉(zhuǎn)杯底部的距離，該距離隨輸棉通道出口與轉(zhuǎn)杯相對位置的變化而發(fā)生改變。

保持抽風(fēng)機功率、轉(zhuǎn)杯形狀、轉(zhuǎn)杯轉(zhuǎn)速、輸棉通道長度和傾斜角等參數(shù)不變，研究輸棉通道出口位置不同時轉(zhuǎn)杯紡紗器中氣流分布規(guī)律及纖維運動規(guī)律的變化，設(shè)計方案見表1。

2 轉(zhuǎn)杯紡紗器的數(shù)值模擬

2.1 數(shù)學(xué)模型

2.1.1 連續(xù)相計算模型

先計算連續(xù)相氣流分布規(guī)律。為簡化研究，不考慮熱交換，假定轉(zhuǎn)杯內(nèi)高速氣流為可壓縮黏性流動，氣體流場為湍流模型[5]，流動滿足質(zhì)量守恒方程和動量守恒方程[6-7]。

式中：uk代表氣流在xk方向的速度分量;ρ為氣體密度，P為氣體壓強;Re為雷諾數(shù);τij為牛頓流體黏性應(yīng)力張量。其中，μ為氣體動力黏性系數(shù);δij為Komecker delta函數(shù);Gk為由平均速度梯度引起的湍流動能K的產(chǎn)生項;Gb為浮力引起的湍流動能的產(chǎn)生項;YM為可壓縮湍流中脈動膨脹對總耗散率的影響;σk和σε分別為湍流動能k與耗散能ε對應(yīng)的普朗特數(shù);C1ε、C2ε、C3ε為常數(shù)。

2.1.2 離散相計算模型

根據(jù)連續(xù)相氣流場模擬計算的結(jié)果擬合顆粒軌跡，顆粒運動滿足牛頓第二定律[9]，計算方程為：

式中：F-表示加速度項，F(xiàn)D（μ--μ-p）為每單位粒子的阻力，μ-表示流體相速度，μ-p表示粒子速度，μ表示流體分子黏度，ρ表示流體密度，ρp表示粒子密度，dp表示粒子直徑，Re為雷諾數(shù)。

2.2 邊界條件設(shè)置

氣流場數(shù)值模擬采用標準k-epsilon模型，近壁面應(yīng)用壁面函數(shù)法，應(yīng)用Simple算法進行氣流場的耦合求解[10-11]。其中邊界條件設(shè)置為：

1）進口邊界：引紗管出口相當于補氣口，不考慮分梳輥的影響，設(shè)輸棉通道入口與引紗管出口為壓力入口，為環(huán)境大氣壓。

2）出口邊界：根據(jù)紡紗工藝，將轉(zhuǎn)杯與假捻盤之間的空隙設(shè)為壓力出口，出口相對壓為-8 000 Pa，方向為垂直壁面方向。

3）固體邊界：將轉(zhuǎn)杯設(shè)為旋轉(zhuǎn)壁面，轉(zhuǎn)速為120 000 r/min，壁面皆設(shè)為無滑移壁面，遵循無滑移邊界條件。

為了簡化研究，忽略粒子質(zhì)量對運動的影響，設(shè)置粒子類型為Massless，將輸棉通道入口設(shè)置為噴射源，其余選項保持默認。

3 模擬結(jié)果與分析

3.1 流場分析

3.1.1 壓強分布規(guī)律

圖2表示不同工況時輸棉通道中心線上的壓強變化。圖3為不同工況時轉(zhuǎn)杯內(nèi)壓強變化情況。圖4為轉(zhuǎn)杯凝聚槽一周（0°～360°）各個角度處的靜壓分布，0°表示輸棉通道中心延長線與凝聚槽交匯處，角度沿順時針方向遞增。由圖2可知，輸棉通道中心線的壓強分布范圍及變化梯度基本一致，各工況下輸棉通道內(nèi)靜壓絕對值由入口到出口逐漸增大，在輸棉通道出口達到最大值。工況a壓強分布梯度最大。由圖3、圖4可知，工況a凝聚槽靜壓分布為-6 300～-7 990 Pa，工況b凝聚槽靜壓分布為-4 080～-8 030 Pa，工況c凝聚槽靜壓分布為-5 100～-7 980 Pa，且靜壓最大值均處于輸棉通道中心延長線與凝聚槽交匯處（0°位置）。工況b時0°位置靜壓值最大，為-4 080 Pa;工況c時0°位置的靜壓值居中，為-5 100 Pa;工況a時0°位置的靜壓值最小，為-6 300 Pa。

在抽氣風(fēng)機和轉(zhuǎn)杯高速旋轉(zhuǎn)的共同作用下，轉(zhuǎn)杯紡紗器內(nèi)部形成強負壓。負壓的大小對纖維的輸送與并合凝聚有重要影響。輸棉通道內(nèi)的負壓有利于纖維從分梳輥處剝離、加速運動至輸棉通道出口并向轉(zhuǎn)杯壁面滑移。進入轉(zhuǎn)杯的纖維在轉(zhuǎn)杯內(nèi)負壓的帶動下沿滑移面凝聚。轉(zhuǎn)杯凝聚槽的負壓越高，越有利于吸附纖維，使纖維在轉(zhuǎn)杯凝聚槽匯集。

3.1.2 氣流速度分布規(guī)律

圖5表示A-A截面氣流速度變化情況。圖6表示不同工況時輸棉通道中心線的氣流速度分布情況。圖7表示不同工況時凝聚槽截面（Y=3.5 mm）上氣流速度矢量變化。圖8表示不同工況時轉(zhuǎn)杯內(nèi)氣流速度矢量分布。由圖5可見，氣流從輸棉通道進入，在輸送通道出口處達到速度最大值并沖入轉(zhuǎn)杯，形成沖擊氣流渦，在高速旋轉(zhuǎn)轉(zhuǎn)杯帶動下隨轉(zhuǎn)杯旋轉(zhuǎn)方向轉(zhuǎn)移至凝聚槽。隨著輸棉通道橫截面積逐漸變小，氣流速度不斷增大，有利于纖維在輸送過程中的伸直。

圖6分別為三種工況時輸棉通道中心線上氣流速度的分布情況，其中橫坐標代表輸棉通道中心線位置，0代表輸棉通道中心線入口點，100代表輸棉通道中心線出口點;縱坐標表示輸棉通道中心線處的氣流速度值。工況a時，輸棉通道出口距離轉(zhuǎn)杯底部的距離（Y值）為13 mm，輸棉通道內(nèi)氣流速度為10.7～118 m/s。工況b時，Y值為9.75 mm，輸棉通道內(nèi)氣流速度為10.2～115 m/s。工況c時，Y值為6.5 mm，輸棉通道內(nèi)氣流速度為12.87～116.15 m/s。

模擬結(jié)果表明，當轉(zhuǎn)杯轉(zhuǎn)速和輸棉通道入口邊界條件均一致時，輸棉通道內(nèi)氣流分布特性與輸棉通道位置有關(guān)。圖6中，曲線a斜率變化最大，即與工況b和c相比較，工況a時氣流變化最大。輸棉通道內(nèi)氣流速度差異大有利于纖維平行伸直度的提高，即工況a最有利于彎鉤纖維及卷曲纖維的伸直。

在實際生產(chǎn)中可以根據(jù)原料特性，例如當原料纖維中彎鉤較多時，可以參考工況a，設(shè)定合適的機件參數(shù)。

由圖7、圖8可知，各工況下，轉(zhuǎn)杯內(nèi)渦流速度分布范圍不同。Y值越小，氣流從輸棉通道轉(zhuǎn)移至凝聚槽時的速度越大，轉(zhuǎn)杯內(nèi)兩股渦流強度差異越大。工況a時，兩股反向渦流邊界在轉(zhuǎn)杯中心處，凝聚槽處氣流速度最大值為47 m/s;工況b時右側(cè)渦流范圍減小，凝聚槽最大速度為56 m/s;工況c時渦流邊界靠近凝聚槽，氣流速度最大值為68 m/s。Y值越小，輸棉通道出口與凝聚槽的距離越小，氣流從輸棉通道出口經(jīng)過轉(zhuǎn)杯壁面，氣流轉(zhuǎn)移至凝聚槽的過程中能量損耗越少。工況a、b時，氣流從輸棉通道出口向轉(zhuǎn)杯轉(zhuǎn)移的過程中與轉(zhuǎn)杯壁面接觸，氣流運動至凝聚槽時速度減小;工況c時，氣流由輸棉通道出口直接進入轉(zhuǎn)杯，凝聚槽處氣流速度最大。

3.2 纖維運動分析

忽略纖維自身的質(zhì)量及纖維與纖維之間、纖維與轉(zhuǎn)杯壁面的相互作用力，假定纖維是從輸棉通道入口自由釋放的運動粒子，運用纖維-氣流兩相流模型模擬轉(zhuǎn)杯紡紗器中纖維的運動狀態(tài)，擬合纖維在轉(zhuǎn)杯紡紗通道的運動曲線。

圖9表示纖維進入輸棉通道入口的位置。在三種工況的輸棉通道入口平面上選取三個位置，分別為靠近輸棉通道上壁面的位置Ⅰ，輸棉通道中心處的位置Ⅱ，靠近輸棉通道下壁面的位置Ⅲ。

圖10—圖12分別表示各工況時不同位置的纖維在紡紗通道運動的軌跡。

工況a時，纖維流入轉(zhuǎn)杯時與轉(zhuǎn)杯壁面接觸的位置離轉(zhuǎn)杯出口較近，從圖10（a）可以看出，輸棉通道靠近壁面處的纖維運動的最大運動速度小于90.68 m/s，纖維離開輸棉通道后受轉(zhuǎn)杯出口負壓的影響，從轉(zhuǎn)杯出口流出形成落棉;從圖10（b）可以看出，輸棉通道中心處的纖維運動速度最大值約120 m/s，纖維進入轉(zhuǎn)杯后沿著轉(zhuǎn)杯滑移面向凝聚槽處匯集;從圖10（c）可以看出，輸棉通道下壁面處的纖維運動速度最大值約為50 m/s，纖維與轉(zhuǎn)杯壁面碰撞后隨B-B面所示氣流沿穿過轉(zhuǎn)杯底部運動至轉(zhuǎn)杯180°（輸棉通道出口處為0°），并在凝聚槽內(nèi)聚集。

工況b時，纖維從輸棉通道出口釋放，與轉(zhuǎn)杯壁面接觸后改變運動方向。從圖11（a）可以看出，纖維沿輸棉通道上壁面運動時的速度最大值小于29.15 m/s，纖維在轉(zhuǎn)杯內(nèi)復(fù)雜渦流的帶動下橫穿過轉(zhuǎn)杯底部，進入凝聚槽;從圖11（b）（c）可以看出，輸棉通道中心處、運動速度最大值約為116 m/s的纖維與輸棉通道下壁面處、速度最大值低于29.15 m/s的纖維運動軌跡基本一致，纖維離開輸棉通道后沿逆時針方向由轉(zhuǎn)杯壁面向凝聚槽匯集。工況b時落棉較小。

工況c時，輸棉通道位置低，氣流直接將纖維引入凝聚槽，纖維頭端與凝聚槽碰撞，易產(chǎn)生前彎鉤。從圖12（a）可以看出，纖維沿輸棉通道上壁面運動時速度較小，纖維在輸棉通道內(nèi)最大速度約為50 m/s，低速運動的纖維離開輸棉通道后速度增大，隨轉(zhuǎn)杯內(nèi)渦流運動至凝聚槽;從圖12（b）可以看出，輸棉通道中心位置運動速度高達118 m/s的纖維離開輸棉通道后，沿順時針方向運動至在凝聚槽;從圖12（c）可以看出，纖維沿輸棉通道下壁面運動時速度小于50 m/s，纖維離開輸棉通道向轉(zhuǎn)杯凝聚槽處匯集后隨轉(zhuǎn)杯運動。

纖維在高速氣流帶動下從輸棉通道入口進入轉(zhuǎn)杯，由于輸棉通道采用漸進縮小式結(jié)構(gòu)，輸棉通道內(nèi)的壓強絕對值與氣流速度沿縮減方向逐漸增大，因此總體上纖維呈加速度運動，在輸棉通道出口處達到最大速度。輸棉通道同一截面處的纖維運動速度不同，靠近輸棉通道壁面的纖維比中心部分的運動速度明顯較慢;越靠近輸棉通道中心線，纖維速度越大。當纖維從輸棉通道出口向轉(zhuǎn)杯轉(zhuǎn)移時，速度慢的纖維易受轉(zhuǎn)杯內(nèi)復(fù)雜渦狀氣流的干擾，在向轉(zhuǎn)杯凝聚槽運動的同時也隨渦流運動，形成彎鉤或者包纏纖維。

4 結(jié) 論

輸棉通道承擔(dān)了將纖維由分梳棍輸送到轉(zhuǎn)杯內(nèi)的任務(wù)，其在紡紗器內(nèi)的幾何位置對轉(zhuǎn)杯及輸棉通道內(nèi)氣流場特征和纖維運動軌跡有重要的影響。本文通過Fluent軟件對抽氣式轉(zhuǎn)杯紡紗器進行三維模擬，分析了三種輸棉通道位置對轉(zhuǎn)杯紡紗器內(nèi)氣流分布特點和纖維運動情況，得出如下結(jié)論：

1）三種工況下輸棉通道流場分布及速度分布規(guī)律基本一致。以工況a為例，氣流在輸棉通道加速運動，運動至出口處速度達到最大值。在輸棉通道與轉(zhuǎn)杯壁面交匯處有小面積的高壓區(qū)產(chǎn)生，有利于纖維流向轉(zhuǎn)杯。

2）工況b時，轉(zhuǎn)杯凝聚槽處壓強較大，有利于纖維在凝聚槽處匯集、凝聚成紗，紗線中纖維纏繞緊密。

3）工況a時，輸棉通道位置最高，輸棉通道內(nèi)壓強和速度變化梯度最大，有利于纖維平行伸直，進入轉(zhuǎn)杯的纖維能較好匯集凝聚;纖維進入轉(zhuǎn)杯時離負壓出口較近，易受負壓影響，落棉量較大。

4）輸棉通道位置對轉(zhuǎn)杯內(nèi)的渦狀氣流分布有一定的影響。工況c時，由于輸棉通道位置逐漸下降，轉(zhuǎn)杯內(nèi)渦流差異越大，不利于轉(zhuǎn)杯轉(zhuǎn)動的穩(wěn)定性，易導(dǎo)致轉(zhuǎn)杯震蕩、機械磨損，縮短轉(zhuǎn)杯壽命。

參考文獻：

[1]張百祥， 周慈念. 轉(zhuǎn)杯紡紗[M]. 北京： 紡織工業(yè)出版社， 1993： 61-65.

ZHANG Baixiang， ZHOU Cinian. Rotor Spinning[M]. Beijing： Textile Industry Press， 1993： 61-65.

[2]KONG L X， PLATFOOT R A. Fiber transportation in confined channel with recirculations[J]. Computers & Structures， 2000， 78（1）： 237-245.

[3]武傳宇， 楊西偉， 陳洪立， 等. 轉(zhuǎn)杯紡紗通道內(nèi)氣體三維流動的模擬[J]. 紡織學(xué)報， 2012， 33（3）： 124-128.

WU Chuanyu， YANG Xiwei， CHEN Hongli， et al. Numerical analysis of 3-D air flow in rotor spinning channel[J]. Journal of Textile Research， 2012， 33（3）： 124-128.

[4]肖美娜， 竇華書， 武傳宇， 等. 紡紗轉(zhuǎn)杯內(nèi)氣流流動特性的數(shù)值分析[J]. 紡織學(xué)報， 2014， 35（12）： 136-141.

XIAO Meina， DOU Huashu， WU Chuanyu， et al. Numerical simulations of air flow behavior in spinning rotor[J]. Journal of Textile Research， 2014， 35（12）： 136-141.

[5]林惠婷， 汪軍. 纖維在輸纖通道氣流場中運動的模擬[J]. 紡織學(xué)報， 2018， 39（2）： 55-61.

LIN Huiting， WANG Jun. Simulation on fiber motion in airflow field of transfer channel[J]. Journal of Textile Research， 2018， 39（2）： 55-61.

[6]曾泳春. 纖維在噴嘴高速氣流場中運動的研究和應(yīng)用[D]. 上海： 東華大學(xué)， 2003.

ZENG Yongchun. Study on Fiber Motion in High Speed Airflow within the Nozzle and Its Application[D]. Shanghai： Donghua University， 2003.

[7]徐曉菊. 對不可壓縮流體連續(xù)性方程的討論[J]. 科技咨詢導(dǎo)報， 2007（28）： 37-38.

XU Xiaoju. Discussion of incompressible fluid continuity equation[J]. Technology Consulting Herald， 2007（28）： 37-38.

[8]鄭永令. 流體流動狀態(tài)與伯努利方程[J]. 大學(xué)物理， 1994， 13（8）： 1-4.

ZHENG Yongling. Fluid flow state and the Bernoulli equation[J]. University Physics， 1994， 13（8）： 1-4.

[9]王福軍. 計算流體動力學(xué)分析[M]. 北京： 清華大學(xué)出版社， 2004： 7-9.

WANG Fujun. Computational Fluid Dynamics Analy-sis[M]. Beijing： Tsinghua University Press， 2004： 7-9.

[10]朱紅均， 林遠華， 謝龍漢. FLUENT流體分析及仿真實用教程[M]. 北京： 人民郵電出版社， 2010： 237-238.

ZHU Hongjun， LIN Yuanhua， XIE Longhan. FLUENT Fluid Analysis and Simulation Practical Tutorial[M]. Beijing： The Peoples Posts and Telecommunications Press， 2010： 237-238.

[11]林惠婷， 汪軍， 曾泳春. 輸棉通道幾何參數(shù)對轉(zhuǎn)杯紡氣流影響的數(shù)值研究[J]. 紡織學(xué)報， 2015， 38（2）： 98-104.

LIN Huiting， WANG Jun， ZENG Yongchun. Numerical study on effect of geometric parameters of transfer channel on airflow in rotor spinning[J]. Journal of Textile Research， 2015， 38（2）： 98-104.

[12]劉超， 楊瑞華， 王鴻博，等. 轉(zhuǎn)杯紡紗通道三維流場的數(shù)值模擬[J]. 紡織學(xué)報， 2016， 39（9）： 145-150.

LIU Chao， YANG Ruihua， WANG Hongbo， et al. Numerical simulation for 3-D flow field of rotor spinning channel[J]. Journal of Textile Research， 2016， 39（9）： 145-150.