樹盤干縮異向性引起應變的測算及分析

付宗營，趙景堯，蔡英春（東北林業大學材料科學與工程學院，哈爾濱150040）

樹盤干縮異向性引起應變的測算及分析

付宗營，趙景堯，蔡英春
（東北林業大學材料科學與工程學院，哈爾濱150040）

研究探討了干燥過程中樹盤受抑制干縮應變、自由干縮應變、彈性應變、黏彈性蠕變應變以及機械吸附蠕變的圖像解析測算法；運用該方法測算了白樺樹盤常規緩慢干燥（含水率分布均勻）過程中干縮異向性引起的弦向各應變，分析了干燥過程中不同含水率階段的應力狀態及應力與各應變的關系。結果表明：應變的圖像解析測算法可滿足精度要求；樹盤干燥至fiber saturation point(FSP)以下，弦向首先受拉伸應力作用，隨著干燥的進行，拉伸應力轉變為壓縮應力；應力方向與各應變對應關系不同，與黏彈性蠕變應變無明顯對應關系，與機械吸附蠕變基本對應。

圖像解析法；樹盤；干縮異向性；干燥應力；弦向應變

1 前言

木材干燥過程中，應力的產生和發展是干燥缺陷產生的主要原因，制約著木材的高效利用。因此，其檢測與研究一直是木材行業攻關的重點。

早在20世紀50年代McMillen[1]就提出利用分層切片法測定木材干燥應力，一直沿用至今。Skaar[2]提出聲發射間接測量法，通過測量干燥時木材輕度開裂而釋放出的彈性波大小和頻率來判斷木材干燥應力的狀況。隨著研究理論的逐步完善以及相關儀器的發展，出現了光學法、應變片法、非接觸法等諸多方法，但由于使用條件、精度等局限未能推廣。長久以來，切片法作為木材干燥過程中應力、應變檢測的常規方法，一直備受研究者們青睞。李大綱等[3]用切片法對楊木高溫干燥過程中各應變進行過測定，并對整個干燥過程中各應變變化特點進行了分析。Erickson等[4]用切片法研究了紅櫟木干燥及平衡處理前后的機械吸附蠕變。切片法的主要缺點是，采用千分尺或游標卡尺等測量工具接觸式測量切片尺寸，大大增加了人為誤差；除干縮應變外，其他各應變數量級都小于10-2，測量工具的精度不能滿足要求；切片翹曲變形后復位困難，影響其長度測量精度。近年來，數字圖像等現代手段逐漸用于木材干燥應力的測量，圖像解析法就是其中之一。

對于樹盤干燥過程中應變的研究，主要集中于受抑制干縮應變（總應變）。Kang W等[5，6]對高頻真空干燥過程中落葉松樹盤的弦向總應變沿徑向分布情況進行過研究。Li等[7]曾在上述研究基礎上，研究了高頻真空干燥過程中落葉松樹盤承受外部壓載荷時弦向應變的分布和變化。目前為止，關于樹盤干燥過程中各種不同類型應變的區分檢測則未見報道。因此本研究旨在運用圖像解析法測算白樺樹盤常規緩慢干燥過程中干縮異向性所引起的弦向各應變，對干燥過程中不同含水率(MC)階段的應力狀態進行分析，并探討應力與各應變的關系。

2 材料與方法

2.1 材料

22年生人工林速生白樺（Betula platyphylla Suk）樹盤若干，平均徑級25 mm，刨光后厚度30 mm，初含水率60%。從中選擇無明顯缺陷的樹盤10塊，其中4塊用于確定初含水率和目標含水率，剩余6塊分兩組作為應變測量試件。實驗在GDS-100恒溫恒濕干燥箱內進行，采用干球溫度恒定，相對濕度階段性遞減的緩慢干燥基準，以確保干燥過程中木材內部含水率分布均勻，從而消除由含水率分布不均所產生的差異干縮應力，只研究由木材干縮異向性所產生的內應力。

2.2 樹盤干燥過程中各種應變的圖像解析測算法

2.2.1 應變切片尺寸的圖像解析法檢測

圖像解析測算法是一種基于圖像分析，非接觸式的點距測量方法，其主要利用高分辨率的數碼設備在木材干燥的不同含水率階段將被測對象連同標尺共同采集圖像，而后將圖像導入ImageJ等軟件對其實際長度進行測量分析。

按圖1所示進行應變切片（30 mm（弦向）× 10 mm（徑向））劃線及長度測點標定，干燥至不同階段對樹盤稱重，并對劃線的表面連同標尺采集圖像；當其干燥至目標含水率時，稱重、采集圖像后沿劃線劈解為應變切片，為避免其翹曲變形，采集圖像時需要將應變切片連同其他部分復位，并用彈力繩捆綁或用其他卡具夾持。使用自行制作的數字圖像采集裝置采集圖像，即將相機固定在鐵架裝置上，鏡頭距測試平面20 cm，標尺垂直于鏡頭水平擺放在測試平面上，同時在測試平面進行劃線定位，以保證每次采集圖像試件放置于同一位置。

圖1 樹盤劃線及應變切片劈解示意圖Fig.1 The cutting diagram of strain slices

2.2.2 應變切片不同狀態時的尺寸

圖2為應變切片不同狀態時的尺寸示意，主要基于Rice、Zhan等的研究[8～11]。

圖2 應變切片不同狀態時的尺寸Fig.2 The size of the strain slice at different drying states

圖2中，L0為干燥前應變切片兩測點間距離；L1為干燥至含水率為26%、18%、10%時，樹盤劈解前應變切片兩測點間距離；L2為干燥至目標含水率，沿劃線將樹盤劈解，即時的應變切片兩測點間距離；L3為將劈解后的應變切片置于恒溫恒濕箱內(介質狀態所對應的平衡含水率等于應變切片劈解時的含水率)，尺寸穩定后兩測點間距離；L4為將測完L3后的各切片用水浸泡24 h，然后汽蒸處理10 h，置于維持上述溫濕度介質條件的恒溫恒濕箱內，自由干縮至穩定狀態后應變切片兩測點間距離。

L2、L3、L4采集圖像前以及應變切片置于恒溫恒濕箱的處理過程中，應將其與其他部分復位并用彈力繩捆綁，以防止其因翹曲變形而影響測量精度。

2.2.3 干燥過程中各應變計算

受抑制干縮應變（總應變）：εR=(L0-L1)/L0；彈性應變：εE=(L1-L2)/L0；黏彈性蠕變應變：εV=(L2-L3)/ L0；機械吸附蠕變：εM=(L3-L4)/L0；自由干縮應變：εF=(L0-L4)/L0=εR+εE+εV+εM。

3 結果與分析

3.1 樹盤干燥過程中干縮異向性引起的應力狀態分析

樹盤干燥至fiber saturation point（FSP）以下，結合水的流失會使樹盤發生干縮。由于弦、徑向干縮差異(一般弦向干縮是徑向的2倍)，弦向干縮會受徑向干縮的抑制作用。該作用力由髓心向外成輻射狀分布（單位環面作用力f1）。另外，干縮異向性引起的各環層差異弦向應變（弦向自由干縮應變與實際受抑制干縮應變之差）與半徑以及弦、徑向自由干縮率的差值成正比[12]，因此，每一環層都要受到與之相鄰外層的微小壓力(單位環面壓力f2)作用。在某一環狀層上選取任意弦向微元為研究對象，受力如圖3a所示，其中F1和F2分別為作用在該微元體上f1和f2的合力，T1、T2為該微元受到的同一環層與其相鄰部分的切向拉力，該二力的合力T與F1和F2的合力F相平衡。

圖3a中，微元體在兩切向拉力下產生弦向拉伸彈性應變，隨著作用時間的增加產生黏彈性蠕變應變，當切向拉力超過比例極限后將產生機械吸附蠕變（塑性變形），超過弦向抗拉強度極限時將開裂。拉伸塑性變形的存在使得弦向收縮受到抑制，而徑向則在壓應力下收縮。隨著干燥的進行，當弦、徑向收縮暫時一致時內應力為0；此后，徑向收縮必定會超過弦向收縮，應力開始轉向，各環層由于對其內部木質徑向收縮的抑制而受到其相鄰內層的拉應力作用（單位環面作用力f1），同時，由于此前外一層比與其相鄰的內層拉伸塑性變形大，將導致前者作用于后者微小的拉伸應力（單位環面作用力為f2）。同樣，在樹盤任意環狀層上選取一弦向微元為研究對象，其受力分析如圖3b所示，F1和F2分別為作用在該微元體上f1和f2的合力；T1、T2為該微元體受到的同一環層與其相鄰部分的切向壓力，該二力的合力T與F1和F2的合力F相平衡。不同部位微元體由拉伸應力向壓縮應力轉變時的含水率狀態，取決于各環層弦向拉伸和徑向壓縮塑性變形的大小，前者與半徑成正比，后者與半徑成反比。

圖3 樹盤應力分析示意圖Fig.3 The diagram of stress analysis in wood disks

3.2 樹盤干縮異向性引起的弦向應力與各應變關系分析

含水率為26%時，樹盤各部位產生了弦向拉伸應力，在其作用下，各部位弦向產生了拉伸彈性應變、拉伸黏彈性蠕變應變，除近樹皮部位外，其他區域產生了拉伸塑性變形（機械吸附蠕變）。心材蠕變值明顯高于邊材，說明心材較邊材易產生蠕變。含水率為18%時，應力已經由拉伸狀態轉變為壓縮狀態，在壓縮應力下產生了壓縮的彈性應變，而其黏彈性蠕變依舊表現為拉伸狀態，說明應力轉向時間短，拉伸黏彈性蠕變尚未完全恢復；對于機械吸附蠕變，其整體表現為壓縮狀態，這是由于前期產生的弦向拉伸塑性變形為局部變形，在干燥至應力轉向之后，該部分變形逐漸恢復甚至產生反變形，而其他部分則產生更大弦向壓縮變形，因此整體表現為壓縮機械吸附蠕變應變。當含水率降至10%時，各部位在弦向壓縮應力作用下表現為壓縮彈性應變，但其應變值較之前略有減小。在壓縮應力作用下，黏彈性蠕變除局部位置因其材性異常而恢復、繼而成壓縮狀態外，其他部位則恢復至0值附近；而壓縮塑性變形逐漸積累，其值也相應增大。不同含水率階段各徑向部位弦向應變見圖4。

4 結語

1）利用圖像解析法測算復位捆綁試件的長度可減小采用千分尺或游標卡尺測量長度過程中的人為誤差，因而計算出的樹盤弦向各應變精度更高。

2）樹盤干燥至FSP以下，首先受拉伸應力作用；隨著干燥的不斷進行，拉伸應力轉變為壓縮應力。該實驗下，應力轉向發生在含水率自26%降至18%的過程中。

3）不同部位微元體由拉伸應力向壓縮應力轉變時的含水率狀態，取決于各環層弦向拉伸和徑向壓縮塑性變形的大小，前者與半徑成正比，后者與半徑成反比。

圖4 不同含水率階段各徑向部位弦向應變Fig.4 The radial distribution of each strain at different MC

4）應力方向與各應變對應關系不同。與黏彈性蠕變應變無明顯對應關系，含水率為18%時，壓縮應力下表現為拉伸的黏彈性蠕變應變；與機械吸附蠕變基本對應，但受材性影響顯著。

[1]McMillen J M.Drying stress in red oak[J].Forest Products Journal,1955,5(2)∶71-76.

[2]Skaar C,Simpson W T,Honeycutt R M.Use of acoustic emissions to identify high levels of stress during oak lumber drying[J].Forest Products Journal.1980,30(2)∶21-22.

[3]李大綱,顧煉百.楊木高溫干燥過程中表層流變特性的研究[J].林業科學,1999,35（1）∶83-89.

[4] Robert W Erickson,Robert Seavey.Energy quantification and mechano-sorptive behavior in the kiln drying of 2.5 cm thick red oak lumber[J].Drying Technology,1992,10(5)∶1183-1206.

[5]Kang W,Lee N H,Cho J H.A radial distribution of moistures and tangential strains within a larch log cross section during radiofrequency/vacuum drying[J].Holz Roh Werkst,2004,62(1)∶59-63.

[6] Kang W,Lee N H.Relationship between radial variations in shrinkage and drying defects of tree disks[J].Journal of Wood Science,2004,50(3)∶209-216.

[7]Li C Y,Lee N H.Effect of external compressive load on tangential strain behavior in Japanese larch log cross sections during radio-frequency/vacuum drying[J].Journal of Wood Science, 2007,53(4)∶291-295.

[8] Rice R W,Youngs R L.The mechanism and development of creep during drying of red oak[J].Holz Roh Werkst,1990,48(1)∶73-79.

[9]Zhan J F,Gu J Y,Shi S Q.Rheological behavior of larch timber during conventional drying[J].Drying Technology,2009,27(10)∶1041-1050.

[10]戰劍鋒,顧繼友,蔡英春,等.木材流變學特性對板材常規干燥開裂、變形的影響[J].林業機械與木工設備,2007，35（10）∶33-36.

[11]戰劍鋒,顧繼友,艾沐野.白樺木材干燥過程橫紋流變特性的初步研究[J].林業科學,2004,40（5）∶174-179.

[12]楊亮慶.楓樺圓盤預處理及干燥特性研究[D].哈爾濱：東北林業大學，2011.

Investigation of tangential strain caused by shrinkage anisotropy using image analytical method

Fu Zongying,Zhao Jingyao，Cai Yingchun
(College of Material Science and Engineering,Northeast Forestry University,Harbin 150040,China)

This study investigated the image analytical method which was used to measure the tangential strains,including practical shrinkage strain,free shrinkage strain,elastic strain, viscoelastic creep strain,and mechano-sorptive(MS)creep,which caused by shrinkage anisotropy during slow conventional drying of wood disks.Also we analyzed the stress in different moisture content(MC)and the relations between stress and strain.The results showed that the image analytical method could meet the requirements in accurate level;below the fiber saturation point(FSP),the wood disks were initially subject to tangential tensile stress;with the decreasing MC,the tensile stress turned into compressive stress.The relations between stress and strain were complex.It had no obviously corresponding connections with the viscoelastic creep strain,and it was almost corresponding with the MS creep.

image analytical method;wood disks;shrinkage anisotropy;drying stress;tangential strain

S781.29

A

1009-1742（2014）04-0025-05

2013-12-06

國家自然科學基金面上項目（31270595）

蔡英春，1960年出生，男，內蒙古寧城縣人，研究方向為木材干燥；E-mail：ychcai@aliyun.com