高 珊，王立海，王 洋

（森林持續(xù)經(jīng)營(yíng)與環(huán)境微生物工程黑龍江省重點(diǎn)實(shí)驗(yàn)室 （東北林業(yè)大學(xué)），哈爾濱 150040）

木材質(zhì)量無(wú)損檢測(cè)是提高木產(chǎn)品質(zhì)量，提高木 材利用率的一個(gè)重要途經(jīng)。應(yīng)用應(yīng)力波法進(jìn)行活立木的木材質(zhì)量檢測(cè)是準(zhǔn)確評(píng)估森林經(jīng)濟(jì)價(jià)值，有效利用森林資源的一項(xiàng)最為直接的手段之一。東北林區(qū)是我國(guó)最大的林業(yè)生產(chǎn)基地。在該林區(qū)大部分樹(shù)木的采伐作業(yè)在冬季仍持續(xù)進(jìn)行。而寒冷的冬季進(jìn)行采伐作業(yè)時(shí)，樹(shù)木內(nèi)水分呈現(xiàn)冰凍結(jié)晶狀態(tài)，這種水分存在狀態(tài)有利于樹(shù)木防御凍害。冰點(diǎn)以下的樹(shù)木其細(xì)胞內(nèi)脂肪和單寧物質(zhì)的含量均有所增加。此時(shí)細(xì)胞液變濃，原生質(zhì)膜形成擬脂層，降低了其透性，因此稱(chēng)此時(shí)的樹(shù)木處于凍結(jié)狀態(tài)［1］。木材是吸濕性極強(qiáng)的生物質(zhì)材料，活立木和剛采伐的新鮮原木都有大量的水分存在。低溫對(duì)活立木應(yīng)力波傳播的影響對(duì)木材質(zhì)量等級(jí)的評(píng)估尤為重要。但是目前，應(yīng)力波檢測(cè)技術(shù)的研究大都是在常見(jiàn)作業(yè)環(huán)境下，針對(duì)一些板材和木質(zhì)產(chǎn)品進(jìn)行應(yīng)力波傳播速度與木材的各種性能之間的關(guān)系研究，而對(duì)立木尤其是低溫下處于凍結(jié)狀態(tài)下的立木的研究的報(bào)道還很少。

Koichi Yamamoto等人（1998）和 William T.Simpson（2001）對(duì)板材和橡膠樹(shù)試驗(yàn)研究了應(yīng)力波傳播速度與含水率的關(guān)系，結(jié)果表明含水率與應(yīng)力波的傳播速度之間存在顯著的負(fù)相關(guān)關(guān)系［2－3］。2005年，van Dyk和Rice［4］進(jìn)行了一項(xiàng)實(shí)驗(yàn)研究來(lái)考察是否可以用超聲波速度來(lái)測(cè)定凍結(jié)材和非凍結(jié)材中的含水率。通過(guò)對(duì)云杉試材（38 mm×100 mm）在五種不同含水率下的徑向波速進(jìn)行測(cè)定后發(fā)現(xiàn)，凍結(jié)材和非凍結(jié)材中的平均徑向波速分別為2 173 m/s和2 077 m/s。凍結(jié)材中的徑向波速比非凍結(jié)材中相比高5%，具有統(tǒng)計(jì)顯著性。2003年，Wang和Ross［5］針對(duì)結(jié)構(gòu)材的應(yīng)力分等對(duì)凍結(jié)材的波速進(jìn)行了試驗(yàn)研究。試材為50根新鮮黃松（Ponderosa Pinus）板材（51 mm×102 mm×2.59 m），含水率均在纖維飽和點(diǎn)以上。應(yīng)力波測(cè)試分別在室溫（15.6℃）和冰窖中（－18.3℃）兩個(gè)溫度條件下進(jìn)行。測(cè)試結(jié)果表明，當(dāng)溫度從常溫降到凍結(jié)狀態(tài)時(shí)，板材的應(yīng)力波速度顯著增加。在室溫時(shí)板材應(yīng)力波的最低、最高傳播速度分別為1 367 m/s和2 978 m/s，而在 －18.3℃時(shí)分別為2 555 m/s和3 833 m/s。凍結(jié)時(shí)的應(yīng)力波傳播速度比常溫時(shí)高出49%。早在2007年我們的研究團(tuán)隊(duì)由王立海帶領(lǐng)對(duì)應(yīng)力波在凍結(jié)狀態(tài)白樺和落葉松活立木中的傳播速度進(jìn)行了測(cè)試，并將研究結(jié)果在2008和2009年發(fā)表。將試驗(yàn)結(jié)果與國(guó)內(nèi)外已有的常溫下應(yīng)力波在活立木中傳播速度的研究結(jié)果進(jìn)行比較時(shí)發(fā)現(xiàn)，應(yīng)力波在凍結(jié)狀態(tài)白樺活立木中的徑向傳播速度明顯高于常溫情況，且傳播速度受含水率影響較大，二者成負(fù)相關(guān)［6－7］。Chan等人（2010）研究溫度和含水率對(duì)200 mm×80 mm×55 mm放射松（Radiata pine）中波速和動(dòng)彈性模量的影響進(jìn)行研究。研究結(jié)果表明，含水率在纖維飽和點(diǎn)以下的板材中的聲速，雖與溫度之間存在著線性關(guān)系，但是纖維飽和點(diǎn)以上的板材，聲速在冰點(diǎn)發(fā)生突變［8］。

1 試驗(yàn)材料及方法

1.1 試驗(yàn)材料

試驗(yàn)材料為東北地區(qū)有代表性的10個(gè)主要樹(shù)種，來(lái)自哈爾濱林業(yè)示范基地的人工林場(chǎng)（經(jīng)度126°，緯度45°）。每個(gè)樹(shù)種分別選擇胸徑接近、長(zhǎng)勢(shì)優(yōu)好、節(jié)疤較少的健康活立木30棵進(jìn)行應(yīng)力波傳播速度測(cè)試。5種針葉樹(shù)，有云杉（Picea asperata）、黑皮油松（Pinus tabulae formis）、落葉松（Larix gmelini）、樟子松（Pinus sylvestris）和冷杉（Abies fabric），5種闊葉樹(shù)，有楊樹(shù)（Populus pseudosimonii）、白樺（Betula platyphylla）、榆樹(shù)（Ulmus pumila）、蒙古櫟（Quercus mongolica）和水曲柳（Fraxinus mandshurica）。

1.2 試驗(yàn)設(shè)備

本試驗(yàn)應(yīng)用儀器主要有應(yīng)力波傳播記時(shí)器（FAKOPP Microsecond Timer，匈牙利進(jìn)口）。該儀器用來(lái)測(cè)量在立木距根部在1 m距離內(nèi)的應(yīng)力波的縱向與徑向傳播時(shí)間。儀器特點(diǎn)靈活小巧，方便攜帶，并且能量損失小，測(cè)量精度高。其規(guī)定的工作溫度范圍在（－10～50℃），但經(jīng)過(guò)預(yù)先測(cè)試證明在超過(guò)該規(guī)定范圍的5～7℃其工作的穩(wěn)定性仍可保證。含水率的測(cè)試主要采用意大利生產(chǎn)的針式木材水分計(jì)（KT－80），用來(lái)測(cè)量立木的邊材含水率。其它便攜儀器有樹(shù)木測(cè)徑儀，皮尺，鋼尺等。由于野外對(duì)立木的含水率的測(cè)量難度相對(duì)較大，因此本實(shí)驗(yàn)中通過(guò)采用幾種木材測(cè)濕儀器對(duì)立木邊材含水率進(jìn)行測(cè)量以后，對(duì)比發(fā)現(xiàn)KT－80的測(cè)量結(jié)果相對(duì)穩(wěn)定，故本實(shí)驗(yàn)采用KT－80進(jìn)行野外立木含水率的測(cè)量。

1.3 試驗(yàn)方法

在2010年1月份，環(huán)境溫度為－15～－10℃進(jìn)行應(yīng)力波傳播速度測(cè)試。首先用鋼卷尺分別在樹(shù)陰面距離地面50 cm和150 cm處做標(biāo)記。在這兩點(diǎn)將應(yīng)力波傳播記時(shí)器（FAKOPP Microsecond Timer）的錐形傳感器以45°角在敲入樹(shù)中（錐形傳感器應(yīng)透過(guò)樹(shù)皮）。固定傳感器后，再次用鋼卷尺精確確定兩個(gè)傳感器之間的距離，然后用錘敲擊傳感器觸發(fā)端，發(fā)送應(yīng)力波，記錄應(yīng)力波順紋理方向（縱向）的傳播時(shí)間，從而獲得該距離內(nèi)的傳播速度。一般取繞立木3個(gè)方向6～8次敲擊的平均時(shí)間作為傳播時(shí)間記錄，以測(cè)得精確的應(yīng)力波的傳播時(shí)間。再用測(cè)徑儀在徑向截面處測(cè)量樣本樹(shù)的直徑，并用Timer測(cè)量應(yīng)力波在樹(shù)木十字交叉徑向的傳播時(shí)間，獲得徑向傳播速度（如圖1與圖2所示）。再用KT－80測(cè)量含水率。由于KT－80的工作條件限制，所獲得的測(cè)量結(jié)果存在一定的不穩(wěn)定性。所以，為保證試驗(yàn)結(jié)果的準(zhǔn)確性，在KT－80的測(cè)量數(shù)據(jù)的基礎(chǔ)上，又通過(guò)采用木材試樣的烘干試驗(yàn)，并根據(jù)2種方法獲得的含水率之間的關(guān)系，對(duì)KT－80獲得的實(shí)驗(yàn)數(shù)據(jù)進(jìn)行修正。所以，本論文中所列含水率數(shù)據(jù)均為修正后的KT－80所得數(shù)值（由于篇幅限制，KT－80測(cè)得含水率的修正報(bào)告將獨(dú)立發(fā)表）。

圖1 縱向應(yīng)力波傳播速度測(cè)試示意圖Fig.1 Set-up for velocity measurements in longitudinal direction

圖2 徑向應(yīng)力波傳播速度測(cè)試示意圖Fig.2 Set-up for velocity measurements in radial direction

2 試驗(yàn)結(jié)果與分析

2.1 試驗(yàn)結(jié)果描述及分析

將試驗(yàn)測(cè)得10個(gè)樹(shù)種的應(yīng)力波縱向傳播速度、徑向傳播速度、立木含水率、胸徑等指標(biāo)數(shù)據(jù)經(jīng)過(guò)Excel進(jìn)行整理和計(jì)算，并通過(guò)SPSS統(tǒng)計(jì)軟件包完成統(tǒng)計(jì)分析。10個(gè)樹(shù)種的總體描述見(jiàn)表1和表2。凍結(jié)狀態(tài)下應(yīng)力波的縱向傳播速度總體平均為是在4 187 m·s－1，徑向傳播速度的總體平均值為2 013 m·s－1。同樣縱向應(yīng)力波的傳播速度明顯高于徑向傳播速度，幾乎是徑向度的2倍。其中應(yīng)力波在闊葉樹(shù)種中的縱向傳播速度平均為4 195 m·s－1，高于針葉樹(shù)種4 180 m·s－1;闊葉樹(shù)種中徑向傳播平均速度為2 039.56 m·s－1，高于針葉樹(shù)種1 988.09 m·s－1。因?yàn)殚熑~樹(shù)種的密度高于針葉樹(shù)種［9］，而應(yīng)力波在樹(shù)木中的傳播速度又和木材的密度成正相關(guān)［3－4］，所以應(yīng)力波在闊葉樹(shù)種中的傳播速度高于針葉樹(shù)種得結(jié)論可以解釋。在縱波試驗(yàn)中，波速又受到早材和晚材的密度、纖維的長(zhǎng)度及聲波的傳播方向的影響［10］。所以，木材內(nèi)的波速，不僅依存于密度，還依存于木材的構(gòu)造。針葉材中，波速受到纖維長(zhǎng)度的影響，在長(zhǎng)纖維材中比短纖維材中快［10］。落葉松和樟子松要比油松和云、冷杉纖維長(zhǎng)度長(zhǎng)，兩者均屬于長(zhǎng)纖維［9］，所以，應(yīng)力波在落葉松和樟子松中的傳播速度均高于油松和云、冷杉。

本試驗(yàn)中應(yīng)力波在闊葉樹(shù)種中的傳播速度水曲柳最高，其次是楊樹(shù)，然后是白樺、榆樹(shù)，蒙古櫟最小。這是因?yàn)殚熑~材波速隨纖維長(zhǎng)度的增加而減小，水曲柳隨密度增加而波速增加，其他闊葉材波速隨密度增大而減?。?］。本試驗(yàn)測(cè)得10個(gè)樹(shù)種的凍結(jié)狀態(tài)下活立木的含水率總體平均值為83.60%，在79.80%和90.70%之間波動(dòng)，其中白樺和云杉的含水率相對(duì)高些，達(dá)90%左右。且立木的陰陽(yáng)面的含水率有所不同，陽(yáng)面的含水率普遍高于陰面2～6個(gè)百分點(diǎn)，且陽(yáng)面的縱向傳播速度普遍低于陰面。

2.2 凍結(jié)狀態(tài)下應(yīng)力波傳播速度影響因素分析

應(yīng)力波在非均勻介質(zhì)和多層介質(zhì)中的傳播規(guī)律也早有研究［11－13］，但木材是各項(xiàng)異性物質(zhì)，這就為應(yīng)用應(yīng)力波法在木材中的傳播規(guī)律的探討帶來(lái)了很多困難。這里我們對(duì)凍結(jié)的活立木中應(yīng)力波的傳播速度的影響因素做數(shù)學(xué)分析。通過(guò)對(duì)凍結(jié)狀態(tài)下的縱向和徑向應(yīng)力波傳播速度進(jìn)行分布檢驗(yàn)可知，各組數(shù)據(jù)均符合正態(tài)分布，且通過(guò)K－S檢驗(yàn)。以縱向應(yīng)力波傳播速度為獨(dú)立變量（y），徑向傳播速度（VR）、含水率（MC）、胸徑（D）作為自變量進(jìn)行回歸統(tǒng)計(jì)。

（1）縱向應(yīng)力波傳播速度影響因素分析。以樟子松的數(shù)學(xué)分析結(jié)果為例，其凍態(tài)下應(yīng)力波縱向傳播速度與立木的含水率的相關(guān)性較高，二者呈現(xiàn)負(fù)相關(guān)，相關(guān)系數(shù)可達(dá)0.89，得到的回歸方程為:y=3490.37－38.64MC（這里的y為應(yīng)力波在凍結(jié)狀態(tài)下的樟子松中縱向傳播速度，MC為凍結(jié)狀態(tài)下樟子松的含水率，統(tǒng)計(jì)置信度為95%，雙尾t檢驗(yàn)顯著性概率都是均遠(yuǎn)小于0.05）。而縱向傳播速度與徑向傳播速度和胸徑的相關(guān)系數(shù)較小，低于0.4（置信度為95%，雙尾t檢驗(yàn)顯著性概率都是均遠(yuǎn)大于0.05），說(shuō)明徑向應(yīng)力波速度與胸徑對(duì)縱向應(yīng)力波傳播速度影響不顯著。冬季樟子松的應(yīng)力波縱向傳播速度與含水率的回歸散點(diǎn)如圖3所示。同時(shí)表3以白樺為例說(shuō)明應(yīng)力波縱向傳播速度、徑向速度以及邊材含水率的相關(guān)性統(tǒng)計(jì)結(jié)果。

對(duì)其余樹(shù)種的統(tǒng)計(jì)分析和回歸顯著性檢驗(yàn)可知，應(yīng)力波在凍結(jié)狀態(tài)下的縱向傳播速度與含水率的相關(guān)性較高，而和徑向傳播速度、胸徑的相關(guān)性較低。由表4得到10個(gè)樹(shù)種的縱向應(yīng)力波傳播速度與含水率的有效回歸方程及其相關(guān)系數(shù)統(tǒng)計(jì)可知:凍結(jié)狀態(tài)下，應(yīng)力波在活立木中的縱向傳播速度主要受含水率的影響，且影響顯著。總的趨勢(shì)是，隨著立木的含水率的降低，應(yīng)力波的縱向傳播速度明顯增高，針葉樹(shù)種中含水率與應(yīng)力波的相關(guān)系數(shù)均在0.90左右，普遍高于闊葉樹(shù)種，一般隨含水率升高1%降低50 m/s左右。圖4、圖5與圖6分別顯示了冬季水曲柳、白樺和落葉松的縱向應(yīng)力波傳播速度與含水率關(guān)系的回歸效果圖。這樣的結(jié)果可能與樹(shù)木內(nèi)部水份的存在狀態(tài)有關(guān)。在完全凍結(jié)的狀態(tài)下，木材內(nèi)的自由水固態(tài)形式存在，此時(shí)樹(shù)木的微觀成分發(fā)生變化，如細(xì)胞內(nèi)脂肪和單寧物質(zhì)都有所增加，這有對(duì)應(yīng)力波的縱向傳播可能產(chǎn)生相應(yīng)的影響，其具體影響方式還有待進(jìn)一步研究。利用不同樹(shù)種在得到凍結(jié)狀態(tài)下的含水率與應(yīng)力波傳播速度之間的關(guān)系，可為解釋凍結(jié)木材中的應(yīng)力波傳播速度高于常溫這一發(fā)現(xiàn)作為基礎(chǔ)。

表1 應(yīng)力波傳播速度在凍結(jié)狀態(tài)下統(tǒng)計(jì)信息Tab.1 Statistical information of SWPV in frozen state

表2 立木在凍結(jié)狀態(tài)下含水率和胸徑的統(tǒng)計(jì)信息Tab.2 Statistical information of MC in frozen state and DBH

表3 10個(gè)樹(shù)種在凍結(jié)狀態(tài)下縱向應(yīng)力波傳播速度和含水率有效回歸方程及其相關(guān)系數(shù)統(tǒng)計(jì)Tab.3 Regression models and correlation coefficients between longitudinal SWPV and MC of 10 species under frozen state

表4 白樺活立木在凍結(jié)狀態(tài)下縱、徑向傳播速度和含水率相關(guān)矩陣Tab.4 The correlation matrix among longitudinal SWPV，radial SWPV and MC of standing trees under frozen state

（2）徑向應(yīng)力波傳播速度影響因素分析。經(jīng)過(guò)統(tǒng)計(jì)分析可知在10個(gè)樹(shù)種的凍結(jié)活立木中，含水率與應(yīng)力波的徑向傳播速度的相關(guān)系數(shù)均在0.40以下（置信度為95%，雙尾t檢驗(yàn)顯著性概率都是均遠(yuǎn)大于0.05，說(shuō)明無(wú)明顯相關(guān)）。由此可得出結(jié)論:在凍結(jié)狀態(tài)下，應(yīng)力波的徑向傳播速度受立木的含水率的影響不明顯。而在取胸徑變化范圍較大的闊葉樹(shù)種－楊樹(shù)和針葉樹(shù)種落葉松與冷杉的試驗(yàn)中，得到徑向應(yīng)力波的傳播速度與胸徑得相關(guān)系數(shù)徘徊在0.60左右。其中楊樹(shù)的為0.69，針葉樹(shù)種冷杉和落葉松分別為0.67和0.63。此結(jié)果說(shuō)明，應(yīng)力波的徑向傳播速度與木材的胸徑可能存在一定聯(lián)系。

圖3 冬季應(yīng)力波在樟子松中縱向傳播速度與含水率關(guān)系的回歸散點(diǎn)圖Fig.3 Regression scatter diagram of relationship between longitudinal SWPV and MC of Pinus sylvestris in winter

圖4 冬季應(yīng)力波在水曲柳中縱向傳播速度與含水率關(guān)系的回歸散點(diǎn)圖Fig.4 Regression scatter diagram of relationship between longitudinal SWPV and MC of Fraxinus mandshurica in winter

圖5 冬季應(yīng)力波在白樺中縱向傳播速度與含水率關(guān)系的回歸散點(diǎn)圖Fig.5 Regression scatter diagram of relationship between longitudinal SWPV and MC of Betula platyphylla in winter

圖6 冬季應(yīng)力波在落葉松中縱向傳播速度與含水率關(guān)系的回歸散點(diǎn)圖Fig.6 Regression scatter diagram of relationship between longitudinal SWPV and MC of Larix gmelini in winter

3 結(jié)論

（1）凍結(jié)狀態(tài)下應(yīng)力波在10個(gè)活立木樹(shù)種中的縱向傳播速度總體平均值為4 187 m/s，徑向傳播速度總體平均值為2 013 m/s。

（2）凍結(jié)狀態(tài)下的應(yīng)力波縱向傳播速與含水率之間存在顯著的負(fù)相關(guān)關(guān)系，相關(guān)系數(shù)多數(shù)在0.8以上。

（3）徑向傳播速度受含水率的影響不是十分顯著，受胸徑的影響有待進(jìn)一步研究。

】

［1］陶大立，靳月華.樹(shù)木越冬傷害［M］.北京:科學(xué)出版社，2005:1－13.

［2］Koichi Y，Othman S，Rokiah H.Nondestructive detection of heart rot of acacia mangium trees in Malaysia［J］.Forest products Journal，1998，48（3）:83-86.

［3］Simpson W T，Wang X.Relationship between longitudinal stress wave transit time and moisture content of lumber during kiln-drying［J］.Forest products journal，2001，51（10）:51-54.

［4］van Dyk H，Rice R W.Ultrasonic wave velocity as a moisture indicator in frozen and unfrozen lumber［J］.Forest Products Journal，2005，55（6）:68-72.

［5］Wang X，Ross R J.Effect of freezing temperature on stress wave speed of green ponderosa pine boards［R］.Technical Report NRRI/TR-2003/10.Prepared for Carter Harvey Holt Ltd.Natural Resources Research Institute，University of Minnesota Duluth，Duluth，MN 2003:4p.

［6］王立海，高 珊，王 洋，等.應(yīng)力波在凍結(jié)狀態(tài)白樺活立木中傳播速度的試驗(yàn)［J］.東北林業(yè)大學(xué)學(xué)報(bào)，2008，36（11）:36 －38.

［7］王立海，王 洋，高 珊，等.凍結(jié)狀態(tài)下應(yīng)力波在長(zhǎng)白落葉松立木中傳播速度的研究［J］.北京林業(yè)大學(xué)學(xué)報(bào)，2009，31（3）:96－99.

［8］Chan J M，Walker J C，Raymond C A.Effects of moisture content and temperature on acoustic velocity and dynamic MOE of radiata pine sapwood boards［J］.Wood Science and Technology，2011，45（4）:609-626.

［9］渡邊治人，張勤麗譯.木材應(yīng)用基礎(chǔ)［M］.上海:上?？茖W(xué)技術(shù)出版社，1985:265－266.

［10］史伯章，尹思慈，阮錫根.聲速與順紋抗壓強(qiáng)度、含水率等物理量的關(guān)系［J］.南京林業(yè)大學(xué)學(xué)報(bào)（自然科學(xué)版），1983（3）:6－12.

［11］葛冬云，陸明萬(wàn).波在各向異性介質(zhì)中傳播規(guī)律的無(wú)網(wǎng)格法數(shù)值模擬［J］.工程力學(xué)，2004，21（5）:122 －125.

［12］劉 妍，張厚江.木質(zhì)材料力學(xué)性能無(wú)損檢測(cè)方法的研究現(xiàn)狀與趨勢(shì)［J］.森林工程，2010，26（4）:46 －49.

［13］董永香，夏昌敬.應(yīng)力波在多層介質(zhì)中傳播特性數(shù)值分析［J］.彈道學(xué)報(bào)，2004，16（3）:28 －33.