木材熱物理參數(shù)理論表達式的推導及其計算值與實測值比較

林 銘 ,楊慶賢 ,饒久平 ,謝擁群 ,廖益強 ,魏 微

(福建農(nóng)林大學材料工程學院 ,福建福州350002)

木材熱物理參數(shù)理論表達式的推導及其計算值與實測值比較

林 銘 ,楊慶賢 ,饒久平 ,謝擁群 ,廖益強 ,魏 微

(福建農(nóng)林大學材料工程學院 ,福建福州350002)

利用木材熱物理參數(shù)3個理論表達式分別計算28種木材的比熱容、導熱系數(shù)和導溫系數(shù) ,并分別與同條件下的實測值相比較.結果表明:28種木材理論計算的平均誤差分別為2.6%、5.7%、6.7%；最大誤差分別為8.0%、12.8%、12.5%.木材熱物理常數(shù)3個理論表達式都不是通常依據(jù)試驗數(shù)據(jù)擬合成的經(jīng)驗方程式 ,而是依據(jù)統(tǒng)計力學和物理學原理 ,經(jīng)過數(shù)學推導得到的理論表達式 ,使之學術上具有內在的理論意義.

木材；熱物理參數(shù)；理論表達式；統(tǒng)計熱力學；類比推理法

木材是三大建筑工程材料之一.木材的比熱容、橫紋導熱系數(shù)和導溫系數(shù)是木材熱工程應用和木材改性、防腐等熱加工處理以及木材科學研究中的重要熱物理參數(shù).由于木材解剖結構的復雜性、化學組成的多樣性、物理性質的各向異性等復雜因素 ,導致從機理上研究木材熱物理參數(shù)有極大困難.因此 ,國內外學者一般都是從試驗上直接研究木材的熱物理參數(shù) ,從中獲得試驗數(shù)據(jù)后 ,或以列表形式、或通過數(shù)學擬合成各種經(jīng)驗方程式[1－13]等表示木材的熱物理參數(shù).

本文將應用統(tǒng)計熱力學和類比法[14－18],推導木材的比熱容、橫紋導熱系數(shù)和導溫系數(shù)3種熱物理參數(shù)的理論表達式.應用這3個理論表達式分別計算28種木材的比熱容、導熱系數(shù)和導溫系數(shù) ,同時分別應用Dunlap[12]和Кириллов[13]比熱經(jīng)驗方程式以及Maclean[4]和Gu-Audrey[11]橫紋導熱系數(shù)公式作對比計算 ,最后均以試驗數(shù)據(jù)進行驗證.

1 木材比熱容理論表達式的推導

木材的各種物理性質主要取決于組成木材細胞壁的物質結構[19].已知纖維素是構成木材細胞壁的主要物質.纖維素大分子的基本化學結構單元是D-吡喃葡萄糖(C6H10O5) ,以β-1 ,4—糖苷鍵聯(lián)接而成的沒有分支的鏈狀分子 ,其化學結構式如圖1所示.

圖1 木材纖維素的化學結構式Fig.1 Chemical structure of wood cellulose

應用X射線衍射 ,結果表明:一切天然纖維素具有相同的晶體結構.纖維素晶體含有結晶區(qū)和無定形區(qū)[20],結晶區(qū)的分子鏈排列最有規(guī)則 ,呈現(xiàn)出清晰的X射線衍射圖.無定形區(qū)的分子鏈排列的規(guī)則性較差 ,但也并非完全無序 ,仍有相當?shù)囊?guī)則性 ,取向大致與纖維軸平行.結晶區(qū)與無定形區(qū)之間沒有明顯界限 ,而是逐漸過渡的.

按照經(jīng)典統(tǒng)計熱力學[21],一個含有N個分子 ,每個分子由S個原子組成的分子晶體 ,可以有3S組不同的簡正振型和頻率.其中每一組振型都有較高和較低兩個頻率與之對應.頻率較高的稱為光學振型 ,頻率較低的稱為聲學振型.由N個聲學振型組成的每一組貢獻出一個德拜項 ,而由N個光學振型組成的每一組貢獻一個愛因斯坦項.于是 ,由一摩爾分子組成的晶體 ,其摩爾熱容為

對于由吡喃環(huán)構成的木材纖維素來說 ,要確定吡喃環(huán)質心的位置 ,需要3個平動自由度 ,又因吡喃環(huán)呈環(huán)狀結構 ,不能看作質點 ,還需3個轉動自由度.這6個自由度滿足聲學色散關系 ,它們各自貢獻一個德拜項.這也意味著(1)式中的a＝2.

20個原子(S＝20) ,共有60個自由度.除了平動和轉動6個自由度外 ,還有54個自由度 ,它們均屬于吡喃環(huán)內各原子的振動自由度 ,包括C—C ,C—H ,C—O ,O—H各鍵的伸縮與撓曲運動.它們均屬于光學振型.各貢獻一個愛因斯坦項.此外 ,C—O鍵強度與環(huán)內鍵強度相近.但環(huán)上的鍵成環(huán)狀結構 ,看作一個整體 ,貢獻德拜項 ,而兩環(huán)之間的這個氧原子因其質量遠小于作為整體的吡喃環(huán)總質量 ,因而不能構成德拜項.因此 ,這3個自由度亦屬光學振型 ,各貢獻一個愛因斯坦項.由此可從(1)式得到絕干木材的比熱表達式

式中μ是葡萄糖基的摩爾質量.

為簡化(2)式的計算 ,可利用熱力學溫度T與攝氏溫度t之間的關系 ,把(2)式C0(T)在T0處展開為泰勒級數(shù) ,取到t的一次項 ,并代入已知數(shù)據(jù) ,整理后再利用絕干木材比熱容C0(T) ,推出含水率為W%的濕木材的比熱容C(t ,W)為

2 木材橫紋導熱系數(shù)理論表達式的推導

2.1 木材細胞結構模型

在分析木材生物細胞微觀形態(tài)的基礎上 ,作為一級近似可把細胞的微觀結構看作一個中空的、細長的、方匣子模型 ,如圖2所示.設細胞長度為d ,橫截面是邊長為b的正方形.中央是細胞腔 ,其橫截面是邊長為a的正方形 ,長度與細胞等長.腔內充滿導熱系數(shù)為λA的空氣.細胞腔四周是厚為(b－a)/a的細胞壁 ,它是由導熱系數(shù)為λC的木材細胞壁物質構成.

2.2 木材導熱的熱阻計算

由于導熱與導電具有許多相似的宏觀屬性 ,根據(jù)類比推理方法 ,可以利用歐姆電阻公式以及電阻的串聯(lián)、并聯(lián)公式來計算木材的熱阻.

現(xiàn)在考慮熱流沿圖2中的X方向流動的一維導熱問題.由圖2可以看出 ,木材在該方向的熱阻R可以看作是由細胞的上下兩個水平壁的熱阻R1,R2與細胞腔中的空氣的熱阻R3并聯(lián)之后再與細胞的前后兩個豎直壁的熱阻R4,R5串聯(lián)而成 ,于是有

由圖2可以看出 ,當熱流從后豎直壁流向前豎直壁(X方向)時 ,必須經(jīng)過上下兩個水平壁和中央空氣腔.由于空氣的導熱系數(shù)λA遠小于細胞壁物質的導熱系數(shù)λC,將導致熱流在前后兩個豎直壁中的流場呈不均勻分布狀態(tài) ,相當于熱流在前后壁中經(jīng)過的橫截面S不再是固定的恒量S＝bd ,而是變化的截面S＝S(x) ,于是熱阻的計算需要積分.為克服S(x)是未知函數(shù)給積分造成的困難 ,引入“有效截面參數(shù)e” ,用“平均”截面ebd替代原幾何截面bd.從而把不均勻熱阻當作均勻熱阻作近似計算.結果如(4)式所示.

圖2 木材細胞結構模型Fig.2 Cell structure model of wood

2.3 橫紋導熱系數(shù)理論表達式

根據(jù)導熱與導電的類比推理 ,可以推出木材的導熱系數(shù)是單位長度、單位橫截面積的木材熱阻的倒數(shù).據(jù)此 ,由(4)式并代入已知數(shù)據(jù)可導出木材橫紋導熱系數(shù)

式中V＝a2db2d ,它是細胞中的空腔體積a2d與細胞總體積b2d的比值 ,稱為木材孔隙率.木材孔隙率是未被木材的細胞壁物質和木材中的水分所占據(jù)的體積與木材總體積之比.據(jù)此這一定義可推出木材孔隙率[14 ,15]

式中ρ是為木材比重.

為確定(5)式中的待定參數(shù)e ,可利用某一已知橫紋導熱系數(shù)試驗值的木材來確定.這里以含水率為12%、密度為900 kg?m－3的母生(Homalium hainanense)的橫紋導熱系數(shù)試驗值等于0.1965 W?m－1?K－1來確定.把這些已知數(shù)據(jù)分別代入(6)式和(5)式 ,便可求出e＝0.45.再把該e值代入(5)式 ,經(jīng)適當整理后 ,可得木材橫紋導熱系數(shù)理論表達式

3 木材橫紋導溫系數(shù)理論表達式的推導

3.1 導溫系數(shù)及其定義式

對于非穩(wěn)態(tài)導熱 ,在相同的溫度梯度下 ,導熱系數(shù)λ越大的物體 ,傳遞的熱量越多 ,其內部各點溫度趨于均勻分布的時間越短；而單位體積熱容量ρc越小的物體 ,其內部各點溫度趨于均勻分布的時間越短.換句話說 ,導熱系數(shù)λ越大、單位體積熱容量ρc越小的物體 ,其內部各點溫度趨于均勻分布的能力越大.或者說 ,傳遞溫度變化的能力越大.我們把表征這種傳遞溫度變化能力的物理量稱為導溫系數(shù)a ,并把它定義為

3.2 木材橫紋導溫系數(shù)理論表達式

根據(jù)導溫系數(shù)定義式(8) ,把(3)式和(7)式代入(8)式 ,易得木材橫紋導溫系數(shù)(10－6m－2?s－1)理論表達式

4 理論計算、試驗驗證與比較

4.1 比熱容計算與驗證和比較

應用(3)式計算28種木材的比熱容 ,并與同條件下的試驗值比較(表1).結果表明:28種木材比熱容計算的平均誤差小至2.5% ,最大誤差也僅有8.0%.

應用Dunlap的經(jīng)驗方程式計算表1中28種木材的比熱容 ,其平均誤差達到9.1% ,最大誤差為15.8%.應用Кириллов經(jīng)驗方程式做同樣計算 ,其平均誤差為10.4% ,最大誤差則達到19.2%.若應用Кириллов的另一個經(jīng)驗方程式計算 ,則其平均誤差為47.8% ,最大誤差為52.9%.

由此可見:本文推導的木材比熱容理論表達式的計算誤差遠比Dunlap和Кириллов經(jīng)驗方程式計算的誤差小得多.

4.2 橫紋導熱系數(shù)計算與驗證和比較

應用(7)式聯(lián)合(6)式計算28種木材的導熱系數(shù) ,并與同條件下的試驗值比較(表1).結果表明:28種木材橫紋導熱系數(shù)理論計算平均誤差為5.7% ,最大誤差為12.8%.

同時應用Maclean經(jīng)驗方程式和Gu-Audrey公式分別計算表1中28種木材的橫紋導熱系數(shù)值 ,兩者的平均誤差分別為15.4%和14.2% ,兩者的最大誤差分別達到29.7%和28.3%.

以上比較計算結果表明:本文推導的木材橫紋導熱系數(shù)理論表達式的計算誤差遠比用Maclean經(jīng)驗方程式和Gu-Audrey公式計算的誤差小得多.

4.3 橫紋導溫系數(shù)理論計算與試驗驗證

應用(9)式計算28種木材的橫紋導溫系數(shù)并與同條件下的試驗值比較(表1).結果表明:28種木材橫紋導溫系數(shù)理論計算平均誤差為6.6% ,最大誤差為12.5%.

5 分析與結論

(1)根據(jù)統(tǒng)計熱力學原理推導的木材比熱理論表達式具有較高的理論計算精確度.經(jīng)典統(tǒng)計熱力學有史以來只應用在單質和一些簡單化合物晶體的比熱容計算.本文嘗試用在化學組成與結構均極其復雜的木材上 ,取得如此成功 ,是經(jīng)典統(tǒng)計熱力學在應用領域上的一個巨大拓展.

(2)依據(jù)導熱與導電宏觀規(guī)律的相似性 ,推出導熱體熱阻的計算公司用于計算木材導熱熱阻 ,從中導出的木材橫紋導熱系數(shù)理論表達式 ,經(jīng)試驗驗證 ,其理論計算誤差較小.

(3)根據(jù)導溫系數(shù)的定義式 ,從木材比熱容和橫紋導熱系數(shù)理論表達式推導木材橫紋導溫系數(shù)公式 ,其理論計算值與試驗值比較 ,兩者吻合程度令人滿意.

(4)本文推導木材橫紋導熱系數(shù)和導溫系數(shù)表達式的理論方法也適用于推導木材徑向和弦向導熱系數(shù)與導溫系數(shù)表達式.

(5)三個熱學參數(shù)的理論值與實驗值基本上相符 ,使本文的研究成果不僅具有實際應用價值 ,而且因他們分別是從微觀機理和宏觀原理上推導所得而具有理論價值.

[1]成俊卿.木材學[M].北京:中國林業(yè)出版社 ,1985:496 ,497.

[2]高瑞堂.木材熱學性質與溫度關系的研究[J].東北林業(yè)大學學報 ,1985 ,13(4):24.

[3]鐘云志 ,鐘義山.木材導溫系數(shù)模型的研究[J].西北林學院學報 ,2003 ,18(2):84－86.

[4]MACLEAN J D.Thermal conductivity of wood[J].Heating ,Piping and Air Conditioning ,1941 ,13(6):380－391.

[5]KOOLMANN F P.Principles of Wood Science and Technology[M].Berlin:Springer-verlarg ,1968:247 ,250.

[6]KELLOGG T M ,WANGAARD F F.Variation in the cell-wall density of wood[J].Wood and Fiber Science ,1969 ,1(3):180－204.

[7]SIAU J F.Transport Processes in Wood[M].New York:Springer-verlarg ,1984:132－150.

[8]ROWLEY F B ,The heat conductivity of wood at climatic temperature differences[J].Heating ,piping ,and Air Conditioning ,1993 ,5:313－323.

[9]OZISIK N M.Basic heat transfer[M].McGraw-Hill Book Company ,New York ,1997:5－65.

[10]SULEIMAN B M ,LARFELDT J B ,LECKNER B ,et al.Thermal conductivity and diffusivity of wood[J].Wood Sci Technol ,1999 ,33(6):465－473.

[11]GU H ,ZINK-SHARP A.Geometric model for softwood transverse thermal conductivity.partⅠ[J].Wood and Fiber Science ,2005 ,37(4):699－711.

[12]DUNLAP F.The specific heat of wood[M] ,Montana:Kessinger.1912.

[13]楊慶賢.木材比熱的統(tǒng)計熱力學研究[J].化學物理學報 ,1991 ,14(6):464－468.

[14]林銘 ,謝擁群 ,楊慶賢 ,等.木材熱導率內在規(guī)律的理論研究[J].福建林學院學報 ,2004 ,24(1):25－27.

[15]林銘 ,陳瑞英 ,楊慶賢 ,等.木材弦向導熱系數(shù)的類比法研究[J].集美大學學報:自然科學版 ,2004 ,9(4):336－340.

[16]林銘 ,饒久平 ,謝擁群 ,等.木材徑向導溫系數(shù)數(shù)學模型[J].福建農(nóng)林大學學報:自然科學版 ,2012 ,41(3):320－323.

[17]林銘 ,謝擁群 ,饒久平 ,等.木材橫紋導熱系數(shù)一種新的表達式的推導及與實測值比較[J].林業(yè)科學 ,2013 ,49(2):108－111 ,49.

[18]YANG Q X.Theoretical expressions of thermal conductivity of wood[J].J Forestry Research ,2001 ,12(1):43－46.

[19]BROWNING D.L.The Chamistry of Wood[M].New York:Springer-verlarg ,1963:123－126.

[20]陳國符.植物纖維化學[M].北京:輕工業(yè)出版社 ,1984:127 ,146－148.

[21]錢學森.物理力學講義[M].北京:科學出版社 ,1962:156 ,162.

(責任編輯:吳顯達)

Derivation of a new expression for wood thermophysics parameters and comparison to experimental value

LIN Ming ,YANG Qing-xian ,RAO Jiu-ping ,XIE Yong-qun ,LIAO Yi-qiang ,WEI Wei
(College of Material Science and Engineering ,Fujian Agriculture and Forestry University ,Fuzhou ,Fujian 350002 ,China)

Theoretical equations ,derived from statistical thermodynamics and analogy ,were used to calculate specific heat ,thermal conductivity and thermal diffusivity of 28 types of wood.Then the results were compared with experimental values to verify the accu-racy of the new equations.Results showed that mean errors for above-mentioned parameters were 2.6% ,5.7%and 6.7%respective-ly ,and maximum errors were 8.0% ,12.8%and 12.5%respectively.These data were all less than results from empirical equations. Different from empirical Maclean and Gu-Audrey equations that were based on experimental results mostly ,derivation method that combines microscopic and macroscopic theories may be applied in other parameter prediction and subjects.

wood；thermal physical parameters；theoretical expression；statistical thermodynamics；analogical reasoning method

S781.37

A

1671-5470(2015)06-0646-05

10.13323/j.cnki.j.fafu(nat.sci.).2015.06.015

2015－03－11

2015－05－21

國家自然科學基金資助項目(30271057).

林銘(1963－) ,男 ,副教授 ,博士研究生.研究方向:木材科學與技術.通訊作者謝擁群(1960－) ,男 ,教授 ,博士生導師 ,研究方向:植物纖維緩沖材料的性能.Email:fjxieyq＠hotmail.com.