花軍 林曉亮 陳光偉

(東北林業(yè)大學(xué)，哈爾濱，150040)

?

磨片精磨區(qū)纖維軌跡模型的建立及磨齒作用頻數(shù)的求解1)

花軍 林曉亮 陳光偉

(東北林業(yè)大學(xué)，哈爾濱，150040)

通過(guò)對(duì)熱磨機(jī)磨片精磨區(qū)內(nèi)纖維運(yùn)動(dòng)的分析，闡述了纖維的流體特性；根據(jù)流體力學(xué)的運(yùn)動(dòng)學(xué)原理求解磨片精磨區(qū)間隙內(nèi)纖維流體的運(yùn)動(dòng)軌跡模型，并建立纖維流體在精磨區(qū)內(nèi)的軌跡方程；結(jié)合磨片的結(jié)構(gòu)參數(shù)，運(yùn)用軌跡方程推導(dǎo)磨齒作用頻數(shù)的數(shù)學(xué)表達(dá)式，并通過(guò)計(jì)算實(shí)際磨片的磨齒作用頻數(shù)和統(tǒng)計(jì)纖維的研磨質(zhì)量數(shù)據(jù)探究二者之間的影響關(guān)系。結(jié)果表明：纖維流體在精磨區(qū)內(nèi)的流動(dòng)屬于庫(kù)特流動(dòng)，并符合Naiver提出的線性滑移模型；其運(yùn)動(dòng)軌跡呈螺旋形，主要受磨片角速度、磨片間隙、磨片齒傾角的影響；磨齒作用頻數(shù)受磨片結(jié)構(gòu)參數(shù)影響，并與纖維的研磨質(zhì)量具有相關(guān)性，合理的磨齒作用頻數(shù)能夠保證較高的纖維質(zhì)量。

纖維分離；精磨；流體運(yùn)動(dòng)；軌跡模型；頻數(shù)

熱磨法是纖維板生產(chǎn)中纖維制備的重要方法，該法利用熱磨機(jī)靜、動(dòng)磨片對(duì)纖維原料進(jìn)行研磨，使纖維原料在研磨過(guò)程中發(fā)生解離。目前，國(guó)內(nèi)外對(duì)于纖維解離的研究主要是從兩個(gè)方面進(jìn)行的。一方面是分析木材原料的形態(tài)、含水率和纖維質(zhì)量、設(shè)備能耗之間的關(guān)系，這類研究主要以安德里茨等纖維分離設(shè)備制造企業(yè)為代表[1]；另一方面是建立磨齒對(duì)纖維作用的力學(xué)模型，通過(guò)模型分析和實(shí)驗(yàn)驗(yàn)證得到磨齒的齒形參數(shù)與纖維分離強(qiáng)度、能量消耗等的聯(lián)系[2-3]。纖維在磨片內(nèi)研磨時(shí)的運(yùn)動(dòng)速度和運(yùn)動(dòng)軌跡決定著纖維的研磨時(shí)間，研磨時(shí)間的長(zhǎng)短影響著纖維受磨齒作用頻數(shù)的高低，從而決定纖維質(zhì)量的優(yōu)劣，而目前對(duì)于纖維在磨片內(nèi)部運(yùn)動(dòng)狀態(tài)的研究較少。

在磨片的精磨區(qū)間隙內(nèi)，從宏觀來(lái)看，纖維的運(yùn)動(dòng)特征具有統(tǒng)計(jì)意義上的平均性和一致性，若只考慮纖維的整體運(yùn)動(dòng)性質(zhì)而不考慮纖維單體的微觀運(yùn)動(dòng)，則其總體行為模式與流體極其相似[4]。纖維的含水率較高，運(yùn)動(dòng)的纖維內(nèi)部以及纖維與磨片之間存在的較大黏性切應(yīng)力和較強(qiáng)的黏附性，這是黏性流體運(yùn)動(dòng)的特點(diǎn)和標(biāo)志。大量的纖維在精磨區(qū)內(nèi)處于強(qiáng)烈的擠壓狀態(tài)，纖維單體之間相互擠壓和粘連，使整個(gè)纖維流體由連續(xù)介質(zhì)組成，因此，纖維既具有流動(dòng)性又具有連續(xù)性。筆者將運(yùn)用流體運(yùn)動(dòng)學(xué)理論對(duì)纖維流體進(jìn)行運(yùn)動(dòng)分析，并建立對(duì)運(yùn)動(dòng)軌跡模型。利用建立的數(shù)學(xué)模型推導(dǎo)出纖維流體流經(jīng)精磨區(qū)時(shí)受磨片、磨齒作用頻數(shù)的計(jì)算公式，通過(guò)分析幾種磨片工作時(shí)磨齒作用頻數(shù)與纖維質(zhì)量的統(tǒng)計(jì)數(shù)據(jù)研究作用頻數(shù)與纖維質(zhì)量的關(guān)系。

1 纖維流體運(yùn)動(dòng)速度求解

1.1 纖維流體運(yùn)動(dòng)控制方程

研磨過(guò)程中磨片間隙設(shè)置為定值，纖維在間隙內(nèi)的體積不會(huì)改變，視為不可壓縮流體?；谇把詫?duì)于纖維流體流動(dòng)性和連續(xù)性的論述，根據(jù)黏性流體運(yùn)動(dòng)學(xué)理論，纖維流體的軌跡方程將由不可壓縮流體流動(dòng)連續(xù)性方程和黏性流體運(yùn)動(dòng)的微分方程導(dǎo)出[5]。

不可壓縮流體流動(dòng)連續(xù)性方程為

(1)

黏性流體運(yùn)動(dòng)的微分方程為

(2)

式中：X、Y、Z為外部作用于單位質(zhì)量流體的體積力沿x、y、z方向的分量；為哈密爾頓算子，k。

1.2 磨片精磨區(qū)間隙內(nèi)纖維運(yùn)動(dòng)速度分析

根據(jù)速度矢量疊加原理，將纖維進(jìn)入磨片精磨區(qū)后的運(yùn)動(dòng)速度分解為繞z軸的周向運(yùn)動(dòng)和沿磨片磨齒的擴(kuò)散運(yùn)動(dòng)，如圖1所示。設(shè)精磨區(qū)內(nèi)某纖維質(zhì)點(diǎn)M的周向速度為u，擴(kuò)散速度為w，二者的合成速度為v，M位置的極徑值為r，M與磨片中心O的連線與x軸夾角為θ，磨齒的齒傾角為φ，磨片的精磨區(qū)內(nèi)徑大小為r0。

圖1 纖維質(zhì)點(diǎn)M的速度分解

1.2.1 纖維的周向運(yùn)動(dòng)

將磨片精磨區(qū)沿著圓周方向展開成直線，使動(dòng)、靜磨片轉(zhuǎn)化成平行平板結(jié)構(gòu)，采用流體平行平板間縫隙流動(dòng)模型進(jìn)行分析[6]。此時(shí)，纖維的周向運(yùn)動(dòng)轉(zhuǎn)換成直線運(yùn)動(dòng)，如圖2所示，若設(shè)纖維質(zhì)點(diǎn)M周向運(yùn)動(dòng)的角速度為ω，則纖維質(zhì)點(diǎn)M沿各坐標(biāo)軸的速度為

(3)

根據(jù)式(1)和式(3)，有?vx/?x=0，由于精磨區(qū)內(nèi)纖維的體積力相比磨片對(duì)纖維的作用力來(lái)說(shuō)很小，忽略體積力的影響。則式(2)可簡(jiǎn)化為

(4)

纖維周向運(yùn)動(dòng)的圓周軌跡上壓強(qiáng)可認(rèn)為處處相等，即?p/?x=0，得

(5)

因此式(4)進(jìn)一步簡(jiǎn)化為關(guān)于z軸坐標(biāo)值的函數(shù)，即

(6)

設(shè)動(dòng)磨片轉(zhuǎn)動(dòng)的角速度為ω0，動(dòng)、靜磨片之間的間隙為hg，由邊界條件

z=0,vz=0;

z=hg,vz=ω0r。

解得纖維質(zhì)點(diǎn)M在磨片間隙內(nèi)的周向速度為

(7)

由式(7)看出，纖維在磨片間隙內(nèi)作周向運(yùn)動(dòng)時(shí)，運(yùn)動(dòng)速度與動(dòng)磨片角速度、磨片間隙及位置半徑有關(guān)，因纖維自身的黏性還會(huì)沿間隙方向(z軸方向)產(chǎn)生速度梯度，即和動(dòng)磨盤磨片接觸的纖維層與磨片同步旋轉(zhuǎn)，與靜磨盤磨片接觸的纖維層靜止不動(dòng)，兩層之間的纖維運(yùn)動(dòng)速度呈線性分布，形成平行于磨片表面的剪切運(yùn)動(dòng)。該流動(dòng)形式在流體運(yùn)動(dòng)學(xué)中稱為庫(kù)特流動(dòng)[7]。其速度分布見圖3。

圖2 流體平行平板間縫隙流動(dòng)模型

圖3 纖維流體流動(dòng)的速度分布

1.2.2 纖維的擴(kuò)散運(yùn)動(dòng)

磨片精磨區(qū)間隙內(nèi)的纖維遇到磨齒的作用時(shí)，受齒傾角的影響，纖維流體會(huì)沿磨齒的延伸方向向外擴(kuò)散。由于周向速度與磨齒齒面的夾角為π/2+φ，則纖維的擴(kuò)散速度w為

(8)

纖維的擴(kuò)散速度表達(dá)式表明磨片精磨區(qū)間隙內(nèi)纖維的擴(kuò)散運(yùn)動(dòng)同樣符合庫(kù)特流動(dòng)，其大小同樣與動(dòng)磨片角速度、磨片間隙及位置半徑有關(guān)，并且受齒傾角的影響。由于擴(kuò)散速度是極徑值r的函數(shù)，而齒傾角又影響著纖維擴(kuò)散路徑的大小。因此齒傾角的設(shè)置決定著纖維在精磨區(qū)研磨時(shí)間的長(zhǎng)短，合理的齒傾角能保證纖維良好的研磨質(zhì)量。

1.2.3 纖維的運(yùn)動(dòng)速度及其修正

由速度合成定理，利用纖維的擴(kuò)散速度和周向速度表達(dá)式對(duì)合成速度的大小進(jìn)行解算。根據(jù)余弦定理得

(9)

解得纖維的合成速度為

(10)

從式(10)中可以看出，纖維在磨片精磨區(qū)間隙內(nèi)的速度受磨片間隙、磨片角速度及磨齒傾角的影響。這些參數(shù)共同影響了纖維的研磨時(shí)間，纖維制備工藝中通過(guò)調(diào)整各參數(shù)的大小，可以控制纖維分離的質(zhì)量。

對(duì)纖維周向速度、擴(kuò)散速度和合成速度表達(dá)式的導(dǎo)出都是建立在磨片表面為水力光滑壁的條件之上的；但是，由于磨片內(nèi)表面并不光滑且存在磨齒，纖維在磨片間隙內(nèi)的流動(dòng)將受到很大阻力，頂層纖維的運(yùn)動(dòng)無(wú)法與動(dòng)磨片同步運(yùn)動(dòng)，而呈現(xiàn)出速度滯后。另一方面，由于動(dòng)、靜磨片之間的間隙很小，動(dòng)磨片的旋轉(zhuǎn)對(duì)靜磨片表層纖維的拖拽作用不可忽略，即底層纖維運(yùn)動(dòng)速度不為零。

基于以上兩個(gè)方面的分析，可以得出與靜磨片接觸的纖維在解離過(guò)程中將會(huì)產(chǎn)生線性滑移。這種滑移又會(huì)減小z軸方向纖維運(yùn)動(dòng)的速度梯度，與Naiver提出的線性滑移模型相符[8-9]，滑移速度的大小與hg成反比，如圖4所示。

在此對(duì)式(7)、(8)和(10)得出得速度表達(dá)式引入?yún)?shù)Bv和d進(jìn)行修正。其中：Bv為考慮磨齒阻力的修正系數(shù)；d為靜磨片表面到速度為零的虛擬點(diǎn)O′的距離，是考慮到磨片表面粗糙程度的修正系數(shù)。修正后的速度表達(dá)式為

(11)

由式(11)可知，磨片精磨區(qū)間隙內(nèi)的纖維運(yùn)動(dòng)速度符合庫(kù)特流動(dòng)模型，速度大小受磨片角速度、磨片間隙、磨齒傾角和磨片表面粗糙程度及磨齒阻力的控制，通過(guò)影響纖維的研磨時(shí)間，進(jìn)而影響纖維的分離質(zhì)量和熱磨機(jī)的能耗[10]。

圖4 纖維流體的線性滑移

2 纖維流體運(yùn)動(dòng)軌跡方程的建立

2.1 軌跡方程的建立

如圖5所示，在極坐標(biāo)下，假設(shè)纖維質(zhì)點(diǎn)M在t時(shí)刻到達(dá)M1，極徑值為r，經(jīng)過(guò)時(shí)間dt后，以速度v到達(dá)點(diǎn)M2處，極徑值為r+dr，對(duì)應(yīng)的轉(zhuǎn)角為dθ，則

(12)

圖5 纖維在磨片間隙內(nèi)的運(yùn)動(dòng)

當(dāng)dθ趨于零時(shí)，cosdθ趨于1，則式(12)簡(jiǎn)化為

(13)

解微分方程(13)得

(14)

當(dāng)t=0時(shí)，r=r0，即C0=r0代入(14)得

(15)

式(15)即為磨片精磨區(qū)間隙內(nèi)纖維運(yùn)動(dòng)軌跡的數(shù)學(xué)模型。從軌跡方程的表達(dá)式可以看出纖維流體運(yùn)動(dòng)的軌跡為螺旋線，其軌跡形狀受磨片轉(zhuǎn)速、磨片間隙和磨齒結(jié)構(gòu)及表面粗糙度影響。

2.2 修正系數(shù)的估算

式(15)中的修正系數(shù)BV和d在實(shí)際生產(chǎn)中難以測(cè)定，但可以通過(guò)生產(chǎn)設(shè)備年產(chǎn)量和熱磨機(jī)及磨片的相關(guān)參數(shù)進(jìn)行估算。

設(shè)熱磨機(jī)的年產(chǎn)量為Va，則纖維的每秒產(chǎn)出率可算，設(shè)為Vs。對(duì)于動(dòng)、靜磨片磨齒與齒槽進(jìn)行體積補(bǔ)償運(yùn)算，設(shè)磨片外徑為R0，齒高為hc，則磨片間纖維流體的出口面積為

S=2πR0(hg+hc)。

(16)

設(shè)纖維流體垂直于磨片方向的出口速度為vs，則有

(17)

此時(shí)應(yīng)取z=hg/2，則將BV和d的估算轉(zhuǎn)換為對(duì)BV(hg+d)的計(jì)算，估算式為

(18)

現(xiàn)以某國(guó)產(chǎn)磨片(內(nèi)徑1 080 mm，外徑1 372 mm)的相關(guān)參數(shù)為例，取磨片精磨區(qū)間隙內(nèi)z=hg/2的纖維流體層，利用Matlab繪制出磨片間隙內(nèi)某纖維質(zhì)點(diǎn)的運(yùn)動(dòng)軌跡如圖6所示。

圖6 纖維在磨片精磨區(qū)間隙內(nèi)的運(yùn)動(dòng)軌跡

通過(guò)模型及圖像可以看出，精磨區(qū)纖維在x-y平面上運(yùn)動(dòng)的軌跡呈螺旋形。當(dāng)z值改變，則不同纖維流體層的運(yùn)動(dòng)軌跡線發(fā)生相應(yīng)變化。通過(guò)幾何求解得到纖維的擴(kuò)散路程sw為

當(dāng)磨片規(guī)格和齒傾角一定時(shí)，sw為定值。由式(8)知，當(dāng)z值減小時(shí)，相同r值時(shí)的擴(kuò)散速度小，因此研磨時(shí)間長(zhǎng)。由于ω0為定值，則軌跡線的圈數(shù)就多。

3 磨片磨齒作用頻數(shù)與纖維質(zhì)量的關(guān)系

熱磨機(jī)在纖維分離過(guò)程中，磨片解離纖維的質(zhì)量和纖維流體的流動(dòng)狀態(tài)之間具有密切的關(guān)系。纖維分離是纖維原料逐次流經(jīng)磨片磨齒表面時(shí)，因磨齒形成的壓縮、剪切和摩擦等多種作用力綜合作用后，造成原料被反復(fù)壓潰、切斷和分絲的結(jié)果。纖維質(zhì)點(diǎn)M在沿著運(yùn)動(dòng)軌跡流動(dòng)的過(guò)程中，由于動(dòng)磨片的高速旋轉(zhuǎn)，M受到磨齒一定頻數(shù)的作用。若頻數(shù)適當(dāng)，則纖維分離的質(zhì)量就好，纖維流體的運(yùn)動(dòng)和所受磨齒作用的頻數(shù)與磨片齒形結(jié)構(gòu)密切相關(guān)；若磨片的齒形結(jié)構(gòu)不合理，原料在磨片內(nèi)的受力頻數(shù)或高或低，都會(huì)造成纖維分離質(zhì)量的下降，即纖維可能被過(guò)度研磨或研磨的不充分。因此，了解磨齒對(duì)纖維流體作用頻數(shù)與纖維研磨質(zhì)量的關(guān)系尤為重要。以下對(duì)磨齒作用頻數(shù)與纖維質(zhì)量關(guān)系進(jìn)行理論推導(dǎo)和實(shí)驗(yàn)研究。

由纖維質(zhì)點(diǎn)M的軌跡方程(15)，取r=R0可以計(jì)算z坐標(biāo)處纖維質(zhì)點(diǎn)M流經(jīng)精磨區(qū)到達(dá)出口所用的時(shí)間為

(19)

設(shè)纖維質(zhì)點(diǎn)M流經(jīng)精磨區(qū)所受磨齒作用的頻次為k，磨齒齒數(shù)為m，磨片轉(zhuǎn)速為n，則

(20)

將式(18)帶入式(20)得

(21)

式(21)即為磨片磨齒作用頻數(shù)k的數(shù)學(xué)表達(dá)式?？梢钥闯觯l數(shù)的大小與磨片轉(zhuǎn)速和磨片的規(guī)格、齒數(shù)、齒傾角及磨片間隙等參數(shù)有關(guān)，vs由式(17)計(jì)算。磨片尺寸和間隙的增大、磨齒齒數(shù)的增多以及齒傾角的加大都會(huì)使磨片磨齒作用頻數(shù)增加，從而加強(qiáng)對(duì)纖維的研磨作用。磨片間隙內(nèi)不同z值處的纖維研磨強(qiáng)度不同，在磨片規(guī)格、齒傾角一定時(shí)，由式(8)可知，z值越大，纖維在精磨區(qū)內(nèi)的擴(kuò)散速度越大，即被磨片研磨的時(shí)間短，經(jīng)歷的磨齒作用頻數(shù)越小，這與式(21)的表達(dá)相一致。實(shí)際生產(chǎn)中應(yīng)當(dāng)綜合考慮磨片結(jié)構(gòu)的各個(gè)參數(shù)，使纖維得到適當(dāng)?shù)淖饔妙l數(shù)，從而防止研磨過(guò)度或不充分。

現(xiàn)以我國(guó)某幾家纖維板生產(chǎn)企業(yè)的幾種磨片為例，對(duì)其研磨出的纖維采用OLYMPUSBX51型生物顯微鏡進(jìn)行測(cè)量和統(tǒng)計(jì)，利用統(tǒng)計(jì)結(jié)果和相應(yīng)的磨片參數(shù)探究磨片磨齒作用頻數(shù)與纖維質(zhì)量的關(guān)系。所用磨片配備的熱磨機(jī)電機(jī)轉(zhuǎn)速均為1 500r/min，外徑1 372mm和1 060mm的磨片的分別以年產(chǎn)量15萬(wàn)m3和10萬(wàn)m3對(duì)vs的估算[11]。頻數(shù)的z軸坐標(biāo)值統(tǒng)一取為平均間隙高度。根據(jù)《木材工業(yè)手冊(cè)》，以早材落葉松的纖維長(zhǎng)度、長(zhǎng)徑比與大青楊、山楊和小葉楊的平均纖維長(zhǎng)度、平均長(zhǎng)徑比分別作為合格纖維樣本尺寸形態(tài)的上下限參考值，即合格纖維長(zhǎng)度為1.27～3.51mm，長(zhǎng)徑比為47.0～60.5[12]。對(duì)每種纖維的長(zhǎng)度合格率和長(zhǎng)徑比合格率進(jìn)行統(tǒng)計(jì)計(jì)算，磨片的結(jié)構(gòu)參數(shù)及纖維質(zhì)量統(tǒng)計(jì)的百分比見表1。

表1 磨片結(jié)構(gòu)參數(shù)及纖維測(cè)量統(tǒng)計(jì)結(jié)果

以頻數(shù)為橫坐標(biāo)，纖維長(zhǎng)度的合格率和長(zhǎng)徑比合格率作為縱坐標(biāo)繪制曲線圖(見圖7)?？梢钥闯?，磨片磨齒作用頻數(shù)對(duì)于纖維長(zhǎng)度合格率的影響呈現(xiàn)由低到高再到低的變化趨勢(shì)，而對(duì)于長(zhǎng)徑比合格率的影響基本與纖維長(zhǎng)度影響趨勢(shì)一致。磨片磨齒作用頻數(shù)較低使纖維研磨得不夠充分，頻數(shù)較高又會(huì)研磨過(guò)度，這都會(huì)使纖維的合格率降低。因此，選擇合理的磨片結(jié)構(gòu)參數(shù)能有效提高纖維的研磨質(zhì)量。

圖7 纖維長(zhǎng)度合格率和長(zhǎng)徑比合格率與磨片磨齒作用頻數(shù)的關(guān)系曲線

4 結(jié)論

磨片精磨區(qū)間隙內(nèi)纖維流體的周向運(yùn)動(dòng)、擴(kuò)散運(yùn)動(dòng)以及合成速度均為庫(kù)特流動(dòng)，速度大小沿著z軸線性分布形成分層剪切，速度大小受磨片角速度、間隙大小及齒傾角的影響；因磨片間隙尺寸較小，又受齒形及磨片表面粗糙程度的影響，纖維流體的運(yùn)動(dòng)符合Naiver提出的線性滑移模型。

纖維運(yùn)動(dòng)的軌跡呈螺旋形，磨片角速度、間隙及磨齒齒傾角影響磨片間隙內(nèi)運(yùn)動(dòng)軌跡的形狀，不同纖維流體層內(nèi)的軌跡線形也不同。z軸坐標(biāo)值較小的纖維流體層內(nèi)，相同極徑值的纖維擴(kuò)散速度小，分離的時(shí)間長(zhǎng)，螺旋軌跡的圈數(shù)也就多。

磨片磨齒作用頻率與磨片的規(guī)格、間隙、轉(zhuǎn)速和磨齒的齒數(shù)、齒傾角等參數(shù)有關(guān)，具有大尺寸、高轉(zhuǎn)速、多齒數(shù)和較大齒傾角的磨片對(duì)纖維的作用頻數(shù)較高。z軸坐標(biāo)值相同的纖維在較大間隙內(nèi)由于擴(kuò)散速度較慢，所受磨齒作用時(shí)間長(zhǎng)，作用頻數(shù)高。纖維分離質(zhì)量和受磨片磨齒的作用頻數(shù)具有相關(guān)性，作用頻數(shù)較小，纖維研磨不充分，作用頻數(shù)較大，纖維研磨過(guò)度，都會(huì)降低纖維的合格率。

由此可見，研究纖維流體在磨片間隙內(nèi)的運(yùn)動(dòng)軌跡及磨齒作用頻數(shù)與纖維研磨質(zhì)量的影響關(guān)系對(duì)獲得較高合格率的纖維制品具有重要意義。實(shí)際生產(chǎn)中選擇合理的磨片結(jié)構(gòu)參數(shù)能夠保證較高的纖維研磨質(zhì)量。

[1] 沈毅，方普新.熱磨機(jī)20年發(fā)展的回顧[J].中國(guó)人造板，2011(3)：130-133.

[2] 董繼先，職艷芳.Z形磨齒新型高濃盤磨機(jī)磨片[J].中國(guó)造紙，2008，27(11)：49-50.

[3]BERGJE，GULLIKSSONME,GRADINPA.Ontheenergyconsumptionforcrackdevelopmentinfiberwallindiscrefining-amicromechanicalapproach[J].Holzforschung，2009，63(2):204-210.

[4]SEDIGHIGILANIMS,NAVIP.Experimentalobservationsandmicromechanicalmodelingofsuccessivedamagingphenomenoninwoodcells’tensilebehavior[J].WoodScienceTechnology，2007，41 (1):69-85.

[5] 林建忠，阮曉東，陳邦國(guó)，等.流體力學(xué)[M].北京：清華大學(xué)出版社，2005：2-314.

[6]SENGERJJ,OUELLERD.Factorsaffectingtheshearforcesinhighconsistencyrefining[J].JournalofPulpandPaperScience，2002，28(11):364-369.

[7]JAMALMS.Fluidflowhandbook[M].London:McGraw-HillEducation，2002:4-45.

[8] 吳承偉，馬國(guó)軍，周平.流體流動(dòng)的邊界滑移問(wèn)題研究進(jìn)展[J].力學(xué)進(jìn)展，2008，36 (3) ：265-266.

[9] 吳承偉，馬國(guó)軍.關(guān)于流體流動(dòng)的邊界滑移[J].中國(guó)科學(xué)G輯:物理學(xué)、力學(xué)、天文學(xué)，2004，34(6)：681-690.

[10] 牛治軍,王寶金,俞敏.熱磨機(jī)的研究現(xiàn)狀與發(fā)展趨勢(shì)[J].林業(yè)機(jī)械與木工裝備，2012(12)：8-12.

[11] 農(nóng)林水產(chǎn)局林業(yè)實(shí)驗(yàn)場(chǎng)編.木材工業(yè)手冊(cè)[M].北京：中國(guó)林業(yè)出版，1991：194-195.

Establishment of Tajectory Model for Fiber in Fine Grinding Zone of Disc and Solution of the Frequency of Gear Grinding Effect//

Hua Jun, Lin Xiaoliang, Chen Guangwei(Northeast Forestry University, Harbin 150040, P. R. China)//

Journal of Northeast Forestry University，2017,45(4)：94-98.

Movement of fiber in fine grinding zone of refiner disc was analyzed and fluid properties of fiber were elaborated. Trajectory model of fiber fluid was solved by kinematics of fluid mechanics through which equation of trajectory in fine grinding zone was established. The formula of the frequency of gear grinding effect in combination with parameter of disc structure was solved through the equation of trajectory. Through calculating the frequency of gear grinding effect of real disc and compiling the statistics of the fiber quality the relationship was analyzed between them. The flow of fiber fluid is Coutte flow and accords with Linear slip model proposed by Naiver. The trajectory is helical and mainly influenced by disc speed, clearance and tooth Inclination. The frequency of gear grinding effect is influenced by parameter of disc structure, and it has a relationship with the properties of fiber. Therefore, reasonable frequency can ensure high quality of fiber.

Separating fiber; Fine grinding; Fluid movement; Trajectory model; Frequency

花軍，男，1959年5月生，東北林業(yè)大學(xué)機(jī)電工程學(xué)院，教授。E-mail：huajun81@163.com。

2016年10月7日。

TS642

1)教育部高等學(xué)校博士學(xué)科點(diǎn)專項(xiàng)科研基金(博導(dǎo)類)課題(20130062110005)。

責(zé)任編輯：戴芳天。