木材的水分吸著和干縮濕脹

摘要：較系統(tǒng)地介紹了木材的水分吸著和干縮濕脹從靜態(tài)到動態(tài)研究的發(fā)展歷程，并歸納了有關恒定溫、濕條件和濕度循環(huán)條件下及濕度循環(huán)過程中3個階段木材水分吸著和干縮濕脹的研究成果及存在的問題，同時對今后需要研究的課題進行了討論。旨在幫助木材科學工作者對該領域的研究進行較全面的了解，并希望他們能夠從中得到啟迪和借鑒，把這方面的研究工作提到一個更新、更高的層次上。

關鍵詞：木材；水分吸著；干縮濕脹；靜態(tài)條件；動態(tài)條件

中圖分類號：S781.33 文獻標識碼：A 文章編號：0439-8114（2013）21-5121-05

Studies on Moisture Sorption and Hygroexpansion of Wood from Static Condition to Dynamic Condition

MA Er-ni

（College of Materials Science and Technology， Beijing Forestry University， Beijing 100083， China）

Abstract： This paper introduced systematically the progress of research on moisture sorption and hygroexpansion of wood from static condition to dynamic condition， and summarized the previous investigations and problems in studying moisture sorption and hygroexpansion of wood under three stages of constant temperature and humidity， cyclic humidity condition and humidity circulating process， respectively. Moreover， some prospects for future study were proposed. The purpose of this paper is to provide an overall understanding about studies in this field， and further to promote the research onto an advance level.

Key words： wood； moisture sorption； hygroexpansion； static condition； dynamic condition

木材是一種天然吸濕性材料，其水分吸著、干縮濕脹等各種物理性質都會受到環(huán)境溫、濕度的影響。由于日常及季節(jié)性溫、濕度的變化，木材總是在不斷地進行吸濕和解吸的同時，尺寸也發(fā)生相應的改變，造成木材在加工和使用時常常出現(xiàn)翹曲、開裂等變形問題。例如，夏季木制門窗不易關閉、冬季木地板又容易開裂等。此外，木材作為結構件利用的許多情況下，還會發(fā)生諸如機械吸濕蠕變等特異行為[1]。因此，對木材的水分吸著和干縮濕脹進行系統(tǒng)而深入的研究尤為重要。

然而，以往對于木材水分吸著與尺寸變形的研究主要是針對恒定的溫、濕條件（靜態(tài)條件，Static condition）進行的。這種手段可以簡化實驗條件和理論分析方法，但是過于理想化。在實際中，木材的加工和使用環(huán)境的溫、濕條件難以保持恒定（動態(tài)條件，Dynamic condition），因而靜態(tài)條件下吸著與變形的研究往往對木制品在日常使用中的指導意義非常有限。為了滿足實際應用的需要，動態(tài)溫、濕條件下的吸著與變形也就成了木材物理學領域研究的必然發(fā)展趨勢[2]。本研究對有關靜態(tài)條件及動態(tài)條件下木材水分吸著和干縮濕脹的研究成果進行了比較系統(tǒng)的歸納，同時提出了一些關于今后開展這方面研究工作的構想，從而幫助木材科學工作者對該領域的研究概況進行較全面的了解，旨在把這方面的研究工作提到一個更新、更高的層次。

1 恒定溫、濕條件下木材的水分吸著和干縮濕脹

有關木材水分吸著和干縮濕脹的研究可以大致分為3個階段。早期階段的研究如上所述，主要是針對恒定的溫、濕條件進行的。在該階段中，對水分吸著的研究基本上是圍繞著木材的吸濕能力展開的，并主要從兩個方面入手：一方面是木材吸濕的大小，即木材在一定溫度、相對濕度條件下所具有的平衡含水率大小，它體現(xiàn)在對木材等溫吸附曲線（Sorption isotherm）的考察上。從20世紀20年代起，先后有學者提出了不同的水分吸著理論（Moisture sorption theory），實現(xiàn)了對木材等溫吸附曲線的定量分析，至今較為著名且常用的有Brunauer-emmett-teller（BET）水分吸著理論[3]、Hailwood-horrobin水分吸著理論[4]及Dent水分吸著理論[5]。Simpson[6]對早期的BET、Hailwood-horrobin等10種水分吸著理論進行了比較，指出這些水分吸著理論都可以表示為如式（1）所示的拋物線函數(shù)形式，其中Hailwood-horrobin理論對實驗數(shù)據(jù)的模擬最為準確。

h/m=A+Bh-Ch2 （1）

式中，h為相對濕度，m為含水率，A、B、C為常數(shù)，其大小因水分吸著理論而異。

木材吸濕能力的另一方面體現(xiàn)在吸濕的快慢上，它屬于水分吸著動力學（Sorption kinetics）的研究范疇。在該問題上，研究者們常以菲克第二定律（Fick’s second law）為基礎來構筑水分吸著動力方程，即將吸著過程中木材內(nèi)部的水分遷移單純地視為非穩(wěn)態(tài)（Non-steady state）條件下水分子的擴散傳遞，認為擴散控制著木材內(nèi)部水分的運動，從而水分吸著可作為擴散現(xiàn)象被準確地模擬。但是，菲克微分方程的求解過程比較復雜，并且不少學者對菲克第二定律在木材-水分系統(tǒng)的適用性提出了質疑，指出木材不是一個理想系統(tǒng)，不符合菲克定律的適用條件[7]；木材的水分吸著過程還應該包含由后期的慢速吸著階段引起的非菲克擴散行為[8]。因此，一些研究者從其他角度出發(fā)對木材的水分吸著動力學進行了討論，如Nakano[9，10]基于化學自催化反應提出的木材非穩(wěn)態(tài)吸濕方程，Zhang等[11]報道的用于預測三維試材吸濕的動力學方程，以及Ma等[12]通過聯(lián)合表面吸著理論和Bradley水分吸著理論獲得的木材水分吸著動力學方程等。

而對于木材的干縮濕脹，早期階段的研究在其影響因素（含水率、溫度等）、各向異性及其機理等方面積累了大量的重要數(shù)據(jù)，馬爾妮等[13]對這方面的研究工作進行了比較系統(tǒng)的歸納。

應該指出的是，恒定的溫、濕條件與實際情況不符，并且早期階段對木材吸著與變形的研究多以最大干縮或濕脹為基礎，由于木材加工和使用環(huán)境的溫、濕度變化通常只在小范圍內(nèi)進行，其含水率也就很少在絕干態(tài)和飽和態(tài)之間跨越。因此，靜態(tài)條件下木材水分吸著和干縮濕脹的研究往往對木材的干燥更有意義，而并不能為木制品的實際應用提供科學可靠的信息。例如，Panshin等[14]早在1980年就指出，椴木在從生材到絕干的干燥過程中表現(xiàn)出很高的干縮值，基于該值的判斷，椴木是不適宜應用的，但實際上，椴木在日常環(huán)境中的變形大小只有預測的1/3～1/2，是非常適合日常使用的。

2 濕度循環(huán)條件下木材的水分吸著和干縮濕脹

鑒于恒定溫、濕條件下吸著與變形的研究對木制品實際應用的局限性，木材水分吸著和干縮濕脹方面的研究迎來了第二階段的發(fā)展，開始考察濕度循環(huán)的條件（溫度恒定）。首先在這方面開展研究的是Harris[15]，他在26 ℃條件下，將木材先后放入相對濕度為65%和95%的環(huán)境中經(jīng)歷3個循環(huán)，每個相對濕度的放置時間為24 h，并在各個循環(huán)的開始和結束階段測定了木材的含水率和弦向尺寸變化。隨后，各國木材物理學者紛紛對濕度循環(huán)條件下木材的吸著與變形性質展開了研究。Stevens[16]提出用參量“移動（Movement）”來表述經(jīng)過干燥的木材在大氣濕度變化時所發(fā)生的尺寸變化，并指出“移動”比生材干燥時所產(chǎn)生的尺寸變化小得多，因而木材的“移動”特性才是反映木制品在日常使用中尺寸穩(wěn)定性的重要指標。

Farmer[17]則進一步提出了采用以下步驟來評估不同樹種窯干木材的“移動”特性：在25 ℃條件下使木材試樣先后于相對濕度為90%和60%的環(huán)境中達到平衡，記錄這兩個濕度環(huán)境下木材的平衡含水率和弦、徑向的尺寸變化，并以相對濕度為90%的尺寸為基數(shù)用百分數(shù)的形式表示出來，即可根據(jù)弦、徑向“移動”之和將木材分為“較小移動（Small movement，小于3.0%）”、“中等移動（Medium movement，3.0%～4.5%）”或“較大移動（Large movement，大于4.5%）”。此外，他還介紹了許多商用鋸材弦、徑向移動的大小。Wu等[18]以25 ℃、35%～55%～75%～85%～93%～75%～35%的相對濕度為一個濕度循環(huán)，考察了循環(huán)周期為12個月和6個月條件下定向刨花板的吸著行為。研究應用Nelson吸著模型通過非線性回歸技術對實驗數(shù)據(jù)進行了分析，并從第一個吸著循環(huán)確定了該模型參數(shù)的大小。研究表明，采用首個循環(huán)數(shù)據(jù)計算得到的模型參數(shù)可適用于對各個循環(huán)過程的模擬。在此基礎上，Wu等[19]探討了該濕度循環(huán)條件下定向刨花板的厚度膨脹情況，包括膨脹滯后（Swelling hysteresis）、變形速率及殘余變形等問題。García等[20]研究了水分循環(huán)處理對木材吸濕能力及尺寸穩(wěn)定性的影響。他們使木材經(jīng)歷5個濕-干水分循環(huán)（相對濕度90%至無水，50 ℃），每個水分條件下放置3 d，并分別在循環(huán)前后的3個不同濕度條件下（相對濕度57.6%、84.2%、90.2%，25 ℃）測定了木材的含水率和弦、徑向尺寸大小。研究發(fā)現(xiàn)，水分循環(huán)處理將造成木材吸濕性及其相應尺寸變化的下降，他們將這種現(xiàn)象分別稱為吸濕老化（Hygroscopic ageing）和變形老化（Dimensional ageing）。此外，木材通常所表現(xiàn)出的含水率與尺寸變化的線性關系在經(jīng)過水分循環(huán)處理后已觀察不到。木材的變形老化在Fan等[21]關于濕度循環(huán)環(huán)境中水泥刨花板的尺寸穩(wěn)定性研究中也得到了確認，這表明水分循環(huán)處理有利于木材尺寸穩(wěn)定性的提高。

可以看出，第二階段的研究雖然涉及了循環(huán)的濕度條件，但其測量均是在循環(huán)前后的水分平衡態(tài)下進行的，缺少對整個循環(huán)過程的考察，對動態(tài)環(huán)境中木材吸著與變形的了解也就十分有限。因此，即使該階段取得了一些標志性的成果，也只能認為它是一個過渡階段，而這個過渡階段的存在可能部分是由于實驗條件的限制，因為對整個動態(tài)過程的把握需要木材含水率、尺寸的實時測量技術。此外，該階段研究中的濕度幾乎都是以矩形波的形式循環(huán)變化的，而研究表明，日常的溫、濕度是近似呈正弦波變化的[22]，這也是該階段的研究中脫離實際情況的第二個不足之處。

3 濕度循環(huán)過程中木材的水分吸著和干縮濕脹

1983年，Chomcharn等[23]把木材干縮濕脹方面的研究引入了第三個階段，即對濕度循環(huán)過程（溫度恒定）的考察，是真正意義上的動態(tài)條件。研究先使美國椴木（Tilia americana）、黃樺（Betula alleghaniensis）、櫻桃木（Prunus serotina）的生材在25 ℃、77%相對濕度條件下達到氣干平衡，然后將它們放在相對濕度為77%～47%的正弦變化環(huán)境中經(jīng)歷若干循環(huán)，循環(huán)周期分別為5.33、10.67、16.00和25.33 h（1∶2∶3∶5），最后使試材在47%相對濕度條件下達到平衡狀態(tài)。在濕度循環(huán)前后的兩個水分平衡態(tài)和濕度循環(huán)過程中，研究測定了試材的含水率及弦、徑向尺寸變化，并以水分膨脹系數(shù)（Moisture expansion coefficient）X、濕度膨脹系數(shù)（Humidity expansion coefficient）Y和水分吸著系數(shù)（Moisture sorption coefficient）Z為參數(shù)，對試材在濕度循環(huán)前后的平衡態(tài)（靜態(tài)）和濕度循環(huán)過程中的非平衡態(tài)（動態(tài)）的橫向干縮濕脹性質進行了計算和比較。結果表明，動態(tài)水分膨脹系數(shù)（Xd）比相應的靜態(tài)值（Xs）高，它在循環(huán)過程中近似為常數(shù)，與循環(huán)周期的大小無關；動態(tài)濕度膨脹系數(shù)（Yd）則還不到相應的靜態(tài)值（Ys）的一半，它隨著循環(huán)周期的增長而增大；而動態(tài)水分吸著系數(shù)（Zd）則隨著循環(huán)次數(shù)的增加及循環(huán)周期的變短而減小。圖1[23]表示了循環(huán)周期為5.33 h條件下椴木的含水率與相對濕度的關系曲線（實線）——木材的動態(tài)吸著曲線，以及Spalt關于椴木的吸濕和解吸等溫線（虛線）——木材的靜態(tài)吸著曲線。圖1中橢圓長軸的斜率代表動態(tài)水分吸著系數(shù)（Zd）的大小，吸濕和解吸等溫線的斜率代表靜態(tài)水分吸著系數(shù)（Zs）的大小。由圖1可知，Zd比Zs小得多，Chomcharn等[23]指出這是由于循環(huán)過程中強烈的吸濕滯后作用所造成的。Chomcharn等[23]的研究表明，在實際使用環(huán)境中，木材的含水率、尺寸變化遠小于平衡狀態(tài)下的結果。因此，使用水分平衡態(tài)下的實驗數(shù)據(jù)來對使用環(huán)境中木制品含水率及尺寸穩(wěn)定性進行預測是不科學的。

Droin-Josserand等[24]提出了一個含有基于濃度變化的水分擴散系數(shù)的數(shù)學模型，該模型可以用于描述相對濕度以一定速度變化的吸著過程。但是，該研究只測試了一個水分循環(huán)，且對解吸過程的計算結果與試驗數(shù)據(jù)相差較遠。

Time[25，26]考察了25 ℃、相對濕度為54%～94%的條件下，以1 d和7 d兩種周期循環(huán)變化環(huán)境中挪威云杉（Picea abies）的含水率及吸濕滯后情況，并構筑了基于菲克定律和水蒸氣壓力驅動的水分傳輸數(shù)學模型。

J？觟nsson[27]測量了20 ℃、相對濕度為40%～80%循環(huán)過程中及自然氣候環(huán)境中膠合層積材橫紋方向上的水分及應變分布情況，該研究為動態(tài)條件下木材在某一方向上不同位置的含水率及尺寸變化的理論分析提供了可靠的實驗數(shù)據(jù)。

在中國，動態(tài)條件下木材水分吸著和干縮濕脹的研究工作尚處于起步階段。Ma等[28]將杉木（Cunninghamia lanceolata）置于25 ℃，相對濕度為45%～75%線性周期變化的動態(tài)環(huán)境中，分別在1.0、2.0和3.5 h 3個周期條件下經(jīng)歷多個循環(huán)，并在該過程中測定了試材含水率和弦、徑向尺寸的變化情況，其結果與Chomcharn等[23]的研究一致。此外，Ma等[28]進一步指出，試材徑向的干縮濕脹落后于其弦向的干縮濕脹。Stevens[16]、Chomcharn等[29]也觀察到了類似的現(xiàn)象，Hart[30]認為這與水分在射線薄壁細胞中的滯留有關。吳義強等[31]考察了尾巨桉（Eucalyptus urophylla）在40 ℃、相對濕度為20%～80%循環(huán)過程中試材弦向的尺度變化情況，并探討了初含水率和循環(huán)次數(shù)對木材脹縮性的影響，研究發(fā)現(xiàn)，對于吸濕過程，循環(huán)次數(shù)的影響明顯；而對于解吸過程，初含水率的影響顯著。

近年來，Ma等[32]將西岸云杉（Picea sitchensis）置于20 ℃，相對濕度為45%～75%的正弦變化環(huán)境中，分別在1、6、24 h 3個周期條件下經(jīng)歷多個循環(huán)，并在該過程中測定了試材含水率和弦、徑向尺寸的變化情況。結果表明，試材含水率及弦、徑向尺寸也呈正弦變化趨勢，但三者都在相位上滯后于相對濕度的變化，如圖2所示。

在收集的實驗數(shù)據(jù)基礎上，研究構建了預測木材動態(tài)吸著的數(shù)學模型，該模型既考慮了木材表面的水分交換，又考慮了木材內(nèi)部的水分擴散，不但模擬效果優(yōu)異，還可用于計算試材厚度方向上水分的實時分布情況。此外，研究采用傅立葉法（Fourier analysis）對試材含水率、尺寸響應的振幅（大?。┖拖辔粶螅ㄋ俾剩┻M行了定量考察，并指出隨著循環(huán)周期的增長，含水率和弦、徑向尺寸相位滯后的趨勢有所下降，而振幅有所增加，Chomcharn等[23]認為這是由于較長的周期為木材對環(huán)境濕度的響應提供了一定的時間，使之能夠隨之跟上相對濕度的變化。在對含水率和尺寸相位滯后的比較分析中，研究發(fā)現(xiàn)木材的尺寸變化快于其含水率的變化。針對該現(xiàn)象，Ma等[33]進一步利用傅立葉法計算得到了試材厚度方向表層、中層、芯層各自含水率的相位滯后大小，并在此基礎上提出了相位滯后假說（Phase lag hypothesis），該假說指出木材厚度方向上各層在相位滯后上的差異是導致木材含水率響應慢于其尺寸響應的原因。

水分非平衡態(tài)下“木材的尺寸變化快于其含水率變化”這一有趣現(xiàn)象的發(fā)現(xiàn)及其解釋的提出十分有助于人們了解木材在日常加工和使用時的一些特異行為并豐富了干縮濕脹理論。但是，該現(xiàn)象的發(fā)現(xiàn)僅是基于一種樹種及試材厚度。眾所周知，木材中水分的擴散受到這兩個因素的影響，對于不同樹種與厚度的試材，木材的尺寸變化是否仍然快于其含水率變化，相位滯后假說是否依然成立，Ma等[32]建立的數(shù)學模型是否還能有效地模擬木材的動態(tài)吸著行為，這些都有待進一步研究證實。

4 研究工作構想

經(jīng)過研究者們20多年的不斷努力，有關動態(tài)條件下木材水分吸著和干縮濕脹的研究已取得了初步進展。但是，從以往研究中獲得的信息還遠遠不能滿足理論和實際上的需求，并且，主要存在以下幾個問題：①涉及溫度循環(huán)乃至溫、濕度共同循環(huán)作用的研究寥寥無幾；②“木材的尺寸變化快于含水率變化”的特異現(xiàn)象及相位滯后假說需要進一步進行樹種和試材厚度的驗證；③第二階段的研究提示人們，水分循環(huán)處理可以提高木材的尺寸穩(wěn)定性，這個效果可否用較短的周期實現(xiàn)，從而節(jié)約時間及能耗，便于應用；④動態(tài)條件下木材的變形的各向異性、殘余應變等對指導木制品在日常加工和使用中具有重要的現(xiàn)實意義，而這些卻尚未得到關注。另一方面，從木材吸著與變形測量中獲得的認識往往是宏觀的，為了從本質上明確上述問題，可從以下兩個層次進行探討。

首先，木材是由纖維素、半纖維素和木素3種主成分構成的天然高分子有機體。纖維素和半纖維素具有親水性，而木素具有疏水性，由于三者具有不同的化學性質和存在形式，它們對溫、濕度變化的響應必然有所不同，這造成它們在發(fā)生含水率或尺寸變化時存在一定的相位差。因此，考察木材主成分各自對溫、濕度循環(huán)作用的響應情況以及它們之間的相互作用，是了解水、熱非平衡態(tài)下木材的吸著與變形行為的關鍵。

其次，由于木材與水分的相互作用伴隨著系統(tǒng)中熱量或能量的變化，因而可以采用水分吸附熱力學的方法來研究木材-水分系統(tǒng)[34]，計算溫、濕度循環(huán)過程中木材中吸著水的熱力學參量如微分吸著熱（QL）、自由能變化（△G）和微分吸著熵（T△S）等。其中，微分吸著熱（QL）對應著水分子與木材實質分子之間的結合能；自由能變化（△G）表示水分子潤脹木材、切斷木材分子間氫鍵，從而暴露木材吸著點所需要做的功；微分吸著熵（T△S）則揭示了吸著在木材實質上的水分子排列的規(guī)則性。在此基礎上，探討動態(tài)條件下吸著水分子的吸著狀態(tài)與吸著環(huán)境，及其與木材實質分子相互作用的情況，最終在分子水平上解析非平衡態(tài)下木材干縮濕脹的機理。

參考文獻：

[1] GROSSMAN P U A. Requirements for a model that exhibits mechano-sorptive behavior[J]. Wood Sci Technol，1976，10（3）：163-168.

[2] 馬爾妮，趙廣杰.木材物理學專論[M].北京：中國林業(yè)出版社，2012.

[3] BRUNAUER S， EMMETT P H， TELLER E. Adsorption of gases in multimolecular layers [J]. J Am Chem Soc，1938，60（2）：309-319.

[4] HAILWOOD A J， HORROBIN S. Absorption of water by polymers： Analysis in terms of a simple model [J]. Trans Faraday Soc，1946，42（B）：84-102.

[5] DENT R W. A multilayer theory for gas sorption. I. Sorption of a single gas[J]. Text Res J，1977，47（2）：145-152.

[6] SIMPSON W T. Predicting equilibrium moisture content of wood by mathematical models[J]. Wood and Fiber，1973，5（1）：41-49.

[7] BABBITT J D. On the differential equations of diffusion [J]. Can J Res，1950，28A（4）： 449-474.

[8] WADS？魻 L. Unsteady-state water vapor adsorption in wood： An experimental study [J]. Wood Fiber Sci，1994，26（1）：36-50.

[9] NAKANO T. Non-steady state water adsorption of wood. Part1： A formulation for water adsorption[J]. Wood Sci Technol，1994， 28（5）：359-363.

[10] NAKANO T. Non-steady state water adsorption of wood. Part 2： Validity of the theoretical equation of water adsorption[J]. Wood Sci Technol，1994，28（6）：450-456.

[11] ZHANG M H， CAZO R， CASSENS D， et al. Water vapor adsorption in kiln-dried red oak[J]. Wood Fiber Sci，2007， 39（3）：397-403.

[12] MA E N， NAKAO T， ZHAO G J. Adsorption rate of wood during moisture sorption processes [J]. Wood Research，2009， 54（3）：13-22.

[13] 馬爾妮，趙廣杰.木材的干縮濕脹——從平衡態(tài)到非平衡態(tài)[J].北京林業(yè)大學學報，2006，28（5）：133-138.

[14] PANSHIN A J， DE ZEEUW C. Textbook of Wood Technology[M]. 4th Edition. New York：McGraw-Hill，1980.

[15] HARRIS J M. The dimensional stability， shrinkage intersection point and related properties of New Zealand timbers[R]. New Zealand， Wellington： Forest Research Institute，1961.

[16] STEVENS W C. The transverse shrinkage of wood [J]. Forest Prod J，1963，13（9）：386-389.

[17] FARMER R H. Handbook of Hardwoods [M]. 2nd Edition. London： Her Majesty Stationary Office，1972.

[18] WU Q L， REN Y K. Characterization of sorption behavior of oriented strandboard under long-term cyclic humidity exposure condition[J]. Wood Fiber Sci，2000，32（4）：404-418.

[19] WU Q L， LEE J N. Thickness swelling of oriented strandboard under long-term cyclic humidity exposure condition [J]. Wood Fiber Sci，2002，34（1）：125-139.

[20] GARC？魱A E L， GRIL J， DE PALACIOS P， et al. Reduction of wood hygroscopicity and associated dimensional response by repeated humidity cycles [J]. Ann For Sci，2005， 62（3）：275-284.

[21] FAN M Z， DINWOODIE J M， BONFIELD P W， et al. Dimensional instability of cement bonded particleboard. Part 2： Behaviour and its prediction under cyclic changes in RH [J]. Wood Sci Technol，2004，38（1）：53-68.

[22] SCHNIEWIND A P. Creep-rupture life of Douglas-fir under cyclic environmental conditions [J]. Wood Sci Technol，1967， 1（4）：278-288.

[23] CHOMCHARN A， SKAAR C. Dynamic sorption and hygroexpansion of wood wafers exposed to sinusoidally varying humidity[J]. Wood Sci Technol，1983，17（4）：259-277.

[24] DROIN-JOSSERAND A， TAVERDET J L， VERGNAUD J M. Modelling the absorption and desorption of moisture by wood in an atmosphere of constant and programmed relative humidity[J]. Wood Sci Technol，1988，22（4）：299-310.

[25] TIME B. Studies on hygroscopic moisture transport in Norway spruce（Picea abies）. Part 1： Sorption measurements of spruce exposed to cyclic step changes in relative humidity[J]. European Journal of Wood and Wood Products，2002，60（4）：271-276.

[26] TIME B. Studies on hygroscopic moisture transport in Norway spruce （Picea abies）. Part 2： Modelling of transient moisture transport and hysteresis in wood [J]. European Journal of Wood and Wood Products，2002，60（6）：405-410.

[27] J？魻NSSON J. Internal stresses in the cross-grain direction in glulam induced by climate variations[J]. Holzforschung，2004， 58（2）：154-159.

[28] MA E N， ZHAO G J， CAO J Z. Hygroexpansion of wood during moisture adsorption and desorption processes [J]. Forestry Studies in China，2005，7（2）：43-46.

[29] CHOMCHARN A， SKAAR C. Moisture and transverse dimensional changes during air-drying of small green hardwood wafers[J]. Wood Sci Technol，1983，17（3）：227-240.

[30] HART C A. Relative humidity， EMC， and collapse shrinkage in wood[J]. For Prod J，1984，34（11/12）：45-54.

[31] 吳義強，彭萬喜.循環(huán)解吸-吸濕處理對人工林桉木脹縮性影響的研究[J].功能材料，2007，38（A09）：3495-3497.

[32] MA E N， NAKAO T， ZHAO G J， et al. Dynamic sorption and hygroexpansion of wood subjected to cyclic relative humidity changes[J].Wood Fiber Sci，2010，42（2）：229-236.

[33] MA E N， NAKAO T， ZHAO G J. Responses of vertical sections of wood samples to cyclical relative humidity changes [J]. Wood Fiber Sci，2010，42（4）：550-552.

[34] SKAAR C. Wood-water Relations[M]. Berlin： Springer-Verlag，1988.