全 文 ：第 51 卷 第 4 期
2 0 1 5 年 4 月
林 业 科 学
SCIENTIA SILVAE SINICAE
Vol. 51，No. 4
Apr．，2 0 1 5
doi:10．11707 / j．1001-7488．20150417
收稿日期: 2014 － 05 － 21; 修回日期: 2014 － 06 － 23。
基金项目: 国家自然科学基金项目(31270595)。
* 蔡英春为通讯作者。
有限差分法逆求木材导热系数*
赵景尧 付宗营 宦思琪 蔡英春
(东北林业大学材料科学与工程学院 哈尔滨 150040)
摘 要: 【目的】 提出一种基于有限差分法逆求木材导热系数的方法，以期弥补传统测量法设备复杂、价格高
昂的不足，为后续建立导热系数的回归方程提供可靠数据。【方法】 运用一维热传导控制方程描述木材升温过程
中内部温度变化。以兴安落叶松弦切锯材为对象，通过试验获得锯材沿厚度方向的温度数据 (温度检测: 使用
NEC Ｒemote Scanner Jr． DC3100 多点信号巡检仪通过埋入锯材的 T 形热电偶获得)，采用有限差分逆求法，合理化
边界条件后编程求解其在不同含水率、不同温度下的径向导热系数(其中，控制方程离散后的导热系数差分矩阵采
用追赶法求解，所有差分方程均在 Matlab2010b 软件编程并运行)，探讨并分析其随含水率、温度的变化规律。【结
果】1) 兴安落叶松锯材径向导热系数计算值沿厚度方向存在一定波动性，但其平均值0. 106 1 W·m － 1 K － 1 (标准差
0. 010 8)符合实际要求，且与理论计算值(0. 110 9，0. 125 2W·m － 1 K － 1 )较为接近; 2) 含水率、温度对导热系数影响
显著，且前者影响高于后者(含水率与温度的 F －检验分别为 126. 942 1，99. 008 3); 含水率、温度共同作用对导热
系数亦存在显著影响(交互作用的 F － 检验为 164. 297 5); 导热系数随含水率的升高而增大，随温度的升高亦增
大; 3) 木材材性与内部含水率分布对导热系数的影响较大。【结论】通过试验获得木材内部可靠的温度分布与变
化数据后，运用有限差分逆求法可快速、准确获得内部与测温点相对应位置的导热系数，尽管计算值存在一定波动
性，但其平均值与理论计算值相吻合，说明运用该方法测算木材导热系数是可行的。相比传统导热系数测量方法，
该方法最大优势在于经济且不受试样尺寸限制; 同时，可以测算试样内部任意层位置导热系数。今后为提高测算
精度，木材内部由于水分迁移产生的热量变化与材性差异应考虑使用该方法; 同时，为推广此方法，将程序可视化
亦是今后研究的方向。
关键词: 含水率; 温度; 导热系数; 有限差分逆求法
中图分类号: S781. 3 文献标识码: A 文章编号: 1001 － 7488(2015)04 － 0134 － 07
Inverse Determination of Thermal Conductivity of
Wood Using Finite Difference Method
Zhao Jingyao Fu Zongying Huan Siqi Cai Yingchun
( School of Material Science and Engineering，Northeast Forestry University Harbin 150040)
Abstract: 【Objective】To develop a precise，quick and economic method for determination the thermal conductivity of
wood，this paper proposed a new numeral method based on the inverse finite difference method，which can overcome the
drawback of traditional methods such as depending on relatively complex and expensive instrumentations． This will provide
reliable data for the development of thermal conductivity’s regression equation． 【Method】A one-dimensional governing
equation of the heat conduction was used to describe the temperature variation inside the lumber during the heating
process． The timbers obtained through tangential cut in larch (Larix gmelinii) plantation，with dimensions of 600 mm in
length，220 mm in width，and 40 mm in thickness，and moisture contents of 28． 3%，41． 2%，62． 3%，were used in the
experiments． Temperature along the thickness direction was obtained via experiments． The detection of temperature was
obtained using NEC Ｒemote Scanner Jr． DC3100 with T-type in wood． With inverse finite difference method and
rationalization of boundary condition，the determination of thermal conductivity in radial direction with different moisture
content and temperature was achieved． Difference matrix of thermal conductivity from the discrete control equations could
be solved by tridiagonal matrix algorithm (TDMA) ． All difference equation were programmed and run in Matlab 2010b．
At the end，the variation law of thermal conductivity with the changing of moisture content and temperature was discussed
第 4 期 赵景尧等: 有限差分法逆求木材导热系数
and analyzed． 【Ｒesult】The results showed that: 1) although some fluctuation along the thickness direction of thermal
conductivity in larch plantation was observed，the average value (0． 106 1 W·m － 1 K － 1，standard deviation ( SD) was
0. 010 8) agreed with theoretical calculating value (0． 110 9，0． 125 2W·m － 1K － 1 ); 2) MC and T had apparent effect on
the thermal conductivity and the former influences slightly higher than the latter ( FMC -test is 126． 942 1，F T -test is
99. 008 3) ． The interaction of MC and T also had apparent effect on the thermal conductivity (FMC × T -test is 164. 297 5) ．
The thermal conductivity increased with increasing MC and T; 3 ) Wood properties and the distribution of MC had
apparent effect on the thermal conductivity． 【Conclusion】Ｒeliable temperature data in wood was obtained via
experiments． The thermal conductivity in the position corresponding to temperature measuring points was obtained quickly
and accurately with IFDM． Although some fluctuation of thermal conductivity was observed，the average value agreed with
theoretical calculating value． It indicted that the determination of thermal conductivity with finite difference method is
feasible． Compared with traditional method，the main advantage of our method is economic and unconstrained in size of
specimen． The thermal conductivity of any layer of wood can be obtained． To improve precision，heat changes with
moisture migration and the difference of wood properties should be considered． To expand the method， program
visualization should also be researched in the future．
Key words: moisture content; temperature; thermal conductivity; inverse finite different method
导热系数是木材热学性质(热容量、热传导、热
扩散与热膨胀)中最重要的物理参数之一(高瑞堂
等，1985; Yang，2001)，也是木材热处理(如木材
干燥、人造板热压等)过程中的重要工艺参数，尤其
在木材干燥方面，直接影响干燥过程中木材升温速
率、内部温度分布以及干燥热质迁移模型的准确性
等(俞昌铭，2011; Walter et al．，2011)。因此，准确
测算木材导热系数，对了解干燥过程中木材内部温
度变化规律及科学合理地实施干燥工艺意义重大。
木材是一种天然高分子多孔材料，其内部的复
杂化学成分和结构给学者研究木材导热性质带来极
大困难。侯祝强(1992)、陈瑞英等(2005)、Hunt 等
(2006)、俞自涛等(2007)、蔡从中等(2009)、林铭等
(2013)采用不同理论与方法对木材内部结构模型
化，根据类比原理与数学分析法建立了导热系数的
理论表达式并探讨了导热系数的影响因素。实际上
导热系数受含水率、温度影响较大，仅由理论导出的
平均导热系数来研究木材升温速率与内部温度变化
必将与实际产生较大误差。因此，建立导热系数依
据含水率、温度变化的回归方程更具实际意义
( MacLean， 1941; Kollmann et al．， 1968; Siau，
1984)，而导热系数的精确测算则是建立回归方程
的重要前提。
导热系数的试验测量方法主要有稳态与瞬态法
(俞自涛，2005; Dupleix et al．，2013)，近年来，有限
差分逆求法也被广泛应用于木材导热系数测算的研
究中(Yeung et al．，1996; Chen et al．，1996; Cai et
al．，2006)。与传统商用平板导热仪器相比，有限差
分法快速、简单、经济，且不受试样尺寸限制; 同时，
可以测算试样内部任意层位置的导热系数。其主要
思想是: 在特定含水率下，测量升温过程中木材内
部不同位置温度随时间的变化，据各层温度变化运
用有限差分法逆求导热系数。但上述相关研究仅探
讨了逆求某一含水率时单一温度段的各层间导热系
数，没有涉及温度、含水率对导热系数的影响; 且边
界条件不合理，如 Cai 等(2006)运用此方法测算冰
冻材的导热系数，计算过程中边界节点方程仅 1 阶
展开，导致边界计算精度与内部不匹配，影响导热系
数的求解精度。
鉴于此，本文以兴安落叶松( Larix gmelinii)木
材为对象，采用有限差分逆求法，合理化边界条件后
编程求解其在不同含水率、不同温度段变化时的径
向导热系数，并分析导热系数随含水率、温度变化的
规律，以期为后续建立导热系数随含水率、温度变化
的二元回归方程提供基础。
1 材料与方法
试材为采自黑龙江省东方红林业局的人工林兴
安落叶松。初选 3 根不同初含水率(约 60%，40%，
30% )，无开裂、腐朽、变形等明显缺陷，规格为
4 200 mm(长) × 225 mm (宽) × 42 mm (厚 ) 的锯
材，4 面刨光并锯解成规格为 600 mm(长) × 220 mm
(宽) × 40 mm(厚)的 3 组弦切板试样(图 1a，每组 6
块，共计 18 块 )，各组试样平均含水率分别为
MC#1 = 28. 3%，MC#2 = 41. 2%，MC#3 = 62. 3%。用环
氧树脂和铝箔隔热卷材对各试样端面及侧面封闭后
再用保鲜膜全面密封，以保证试验时热量只沿厚度方
向传递并防止水分流失 (试验前后水分流失不足
531
林 业 科 学 51 卷
1. 2% )。由于温度巡检仪测点数量限制，每次仅对 2
块试样进行试验(每组在相同条件下进行 3 次)，其余
试样放入冷库中备用。3 组不同平均含水率、初始温
度为室温的试样分别在 50，70，90 ℃的环境温度(湿
度均为 100%，风速 2 m·s － 1)下试验并测算。为减小
材性对试验结果的影响，每次试验结束待试样温度降
到室温后，再置入下一水平温度环境中试验。试验
前，在试样侧面用直径为 1 mm 的钻头钻孔(具体位
置见图 1b)并称重，然后将 T 形热电偶埋入试样中，
再将试样放入温度已经稳定为设定值的 DS-408 型恒
温恒湿箱中加热升温。使用 NEC Ｒemote Scanner Jr．
DC3100 多点信号巡检仪通过埋入木材的 T 形热电偶
检测材温(每分钟记录 1 次)。自冷库中取出的试样，
先室温解冻并稳定至室温后再进行上述操作。试验
结束后，对试样按 GB /T 1933—2009《木材密度测定
方法》测量基本密度，采用绝干称重法检测含水率。
试样初始温度(室温)为 16 ～ 18 ℃。
2 模型建立
2. 1 试样加热升温过程中逆求导热系数的控制
方程
模型建立前，提出合理的假设条件是必要的。
由于试样厚度小于宽度、远小于长度，且其端口和侧
面都进行了隔热处理，因此，假定温度仅沿厚度方向
传递; 同时，假定木材结构、物理化学特性在厚度方
向上是对称的，厚度方向上的中心层面即为绝热面;
假定试样初始温度分布均匀; 假定升温阶段试样含
水率保持恒定，即试样表面只存在热交换; 假定试
样外部环境不随其热量变化而变化，即环境所有状
态参数是恒定的。
试样内部温度控制方程 (逆求导热系数方程)
如下:
图 1 锯解示意与热电偶测温部位分布
Fig． 1 Sawing diagram and thermocouple arrangement
ρd cp
T
τ
= 
x
(λ T
x
)。 (1)
式中: λ 为导热系数 (W·m － 1 k － 1 );T 为木材温度
(℃ ); τ 为时间( s); ρd 为木材基本密度(kg·m
－ 3);
cp 为木材比热( J·kg
－ 1K － 1)，cp = 2 000 + 8 ． 71MC +
4 ． 98T (Steinhagen et al．，1988); x 为沿木材厚度方
向从表面到内部任意一点的距离(mm)。
显然，在式(1)中，当基本密度 ρd 、比热 cp 、各
层温度 T /x 、不同时刻温度 T /τ 变化已知的情
况下，导热系数可解。
初始条件:
T = T(0，x)。 (2)
边界条件:
木材表面处( x = L，L 为木材厚度的一半)的边
界条件:
－ λ T
x
= h T s － T( )e 。 (3)
式中: h 为表面传热系数 (W·m － 2K － 1 )，由于表面
传热系数受风速影响较大(本研究风速恒定)，为
此，假定 3 组温度水平下表面传热系数均为 h = 9. 6
W·m － 2K － 1 ( Pordage et al．，1999); T s 为表面温度
(℃ ); T e 为环境温度(℃ );
木材厚度中心层( x = 0)的边界条件:
T
x
= 0。 (4)
2. 2 逆求导热系数方程的有限差分法
运用有限差分法对式(1)(逆求导热系数方程)
进行求解。首先，对空间区域进行离散化，将区域
［0，L］分成 J 个相同的单元，节点依次为 j = 1，…，
J，共 J 个，空间步长记作 Δx，xj = jΔx(图 2); 将时间
631
第 4 期 赵景尧等: 有限差分法逆求木材导热系数
区域分成 I 个单元，节点依次为 i = 1，…，I，共 I 个，
时间步长记作 Δτ，τ i = iΔτ。其次，写出所有方程的
有限差分格式。
图 2 试样与离散示意
Fig． 2 Diagram of specimen and discretization
a．试样 Specimen; b．沿径向空间离散
Space discretization in radial direction．
1) 内部节点差分方程(0 ＜ j ＜ J):
4ρcp
Δ( )x 2
Δτ T
i+1
j － T( )
i
j = λ
i
j －1(T
i
j －1 － T
i
j +1) +
4λ ij(T
i
j +1 － 2T
i
j + T
i
j －1) +
λ ij +1(T
i
j +1 － T
i
j －1)。 (5)
2) 为保证模型精度，表面节点采用 2 阶差分方
程( j = J):
ρcp
Δ( )x 2
Δτ
(Ti+1J － T
i
J) = λ
i
j －1(T
i
j －1 － T
i
j) +
h(T e － T
i
j)Δx。 (6)
3) 心层节点差分方程( j = 1):
ρcp
Δ( )x 2
Δτ
(Ti+1j － T
i
j) = 2λ
i
j(T
i
j +1 － T
i
j)。 (7)
令矩阵 d = (λ j) J ×1，矩阵 b = ( bj) J ×1，矩阵a =
(aj，k) J × J，以上有限差分方程即可表示为 ad = b，如
下式所示:
a1，1 0 0 … 0
a2，1 a2，2 a2，3 … 0
 a3，2  a3，4 
0  aj，j －1 aj，j a j，j +1
0 0 0 … aJ，


J J × J
λ i1
λ i2

λ iJ －1
λ

i
J J ×1
=
bi1
bi2

biJ －1
biJ － h(Te － T
i
J)Δ

x
J ×1
。 (8)
其中，矩阵 a，b 元素表示如下:
bj = ρcp
Δ( )x 2
Δτ T
i+1
j － T( )
i
j ，j = 1，J; (9)
bj = 4ρcp
Δ( )x 2
Δτ T
i+1
j － T( )
i
j ，1 ＜ j ＜ J; (10)
aj，j = 2(T
i
j +1 － T
i
j)，j = 1; (11)
aj，j －1 = T
i
j －1 － T
i
j +1，1 ＜ j ＜ J; (12)
aj，j = 4(T
i
j +1 － 2T
i
j + T
i
j －1)，1 ＜ j ＜ J; (13)
aj，j +1 = T
i
j +1 － T
i
j －1，1 ＜ j ＜ J; (14)
aj，j = T
i
j －1 － T
i
j，j = J。 (15)
显然，矩阵 a 为“三对角矩阵”，因此可通过追
赶法 ( TDMA， tridiagonal matrix algorithm ) 求 出
( Kreith et al．，2010; 王 正 林 等， 2009 )。通 过
Matlab2010b 软件对上述矩阵式编程并运行解出导
热系数(矩阵 d)。
3 结果与讨论
图 3 为表 1 所示组别 #1 (平均含水率 MC 为
28. 3%的试样分别在 50，70，90 ℃的环境温度下)升
温过程中不同深度位置(图 1)实测温度随时间的变
化曲线，其中，各时刻温度为 6 块试样同部位温度均
值。从图中可清晰看出: 升温过程中试样各层温度
随时间变化相似，即初期升温较快，随时间的延续变
化渐缓，到后期逐渐趋于一致达到目标温度。
以 T1(距表层 20 mm 的心层温度)与 T5 (距表
层 4 mm)为例，初始阶段温差逐渐增大，3 个温度水
平下均在 τ = 30 h 左右温差达到最大，分别为 8. 1，
13. 1，17. 9 ℃，随后此温差逐渐变小。
图 4 显示由矩阵关系式(8)编程计算得出的各
组别不同厚度位置( x = 4，8，12，16，20 mm)与时刻
(τ = 10，25，50，70，100 min)的导热系数(每组 6 块
试样同位置同时刻的导热系数均值)分布，显示了
导热 系 数 的 较 大 变 异 性。以 图 4 组 别 # 1 － 50
(MCave = 28. 3%，T = 50 ℃ )的导热系数分布为例，
同时刻不同位置导热系数、同位置不同时刻导热系
数都存在不同程度差异，且变化无明显规律。导热
系数这 2 项差异中的前者，主要由不同部位的密度、
构造、树脂含量等差异引起(Gu et al．，2005; 2006;
Vay et al．，2013; 成俊卿，1985); 而后者，可能是由
于试样在加热升温过程中其内部水分产生某种程度
迁移引起，而式(1)中并未考虑试样含水率分布的
变化，因而产生了误差。总之，前者反映了客观存在
即材性差异的影响，后者则显示了试验及求解误差
即模型本身的影响，也表明了后续改进研究的必要，
如考虑木材内部结构、含水率分布的影响等。
731
林 业 科 学 51 卷
图 3 试样升温过程各层温度随
时间变化曲线(组别#1)
Fig． 3 Variation of temperature as a function of time
during heating( group#1)
表 1 为试样不同平均含水率时在不同温度环境
下的平均导热系数及其标准差。从表 1 可以看出，
总平均导热系数为 0. 106 1 W·m － 1 K － 1，标准差为
0. 010 8，与按 Gu(2005)几何模型公式计算得出的
平均导热系数 0. 110 9 W·m － 1 K － 1、林铭(2013)的
理论计算导热系数 0. 125 2 W·m － 1K － 1接近，说明由
本模型计算得出的平均导热系数有效。
表 2 为含水率与温度对导热系数影响的方差分
析。从表 2 可以看出，含水率、温度对导热系数均存
在显著性影响 ( P ＜ 0. 01)，其中，含水率的影响高
(FMC = 126. 942 1 ＞ FT = 99. 008 3); 二者交互作用
亦对导热系数存在显著影响。将图 4 所显示的导热
系数按含水率和温度分组求均值(同一温度同一含
水率的所有不同时刻、不同层的导热系数均值)并
绘制图 5。
图 4 各组别下不同位置与时刻的导热系数分布
Fig． 4 Variation of thermal conductivities with the changes of depths and moments for different groups
图 5 和表 1 表明: 1) 在试样温度一定时，导热
系数随含水率的升高而增大(例如在 90 ℃环境温
度下，λ #1 － 90 = 0. 099 9 W·m
－ 1 K － 1，λ #2 － 90 = 0. 112 9
W·m － 1K － 1，λ #3 － 90 = 0. 139 8 W·m
－ 1 K － 1 )，原因为，
水分子的导热能力(λ≈0. 6 W·m － 1 K － 1 )约是木材
分子的 3 倍、空气分子的 23 倍，随着木材内水分的
增加，孔隙内空气被水分替代比例增大，因而产生上
述结果; 2) 在试样含水率恒定时，导热系数亦随温
度升高而增大(例如组别#3，即 MCave = 62. 3%状态
下，导热系数随温度升高而增大，λ #3 － 50 = 0. 113 1
W·m － 1K － 1，λ #3 － 70 = 0. 121 0 W·m
－ 1 K － 1，λ #3 － 90 =
0. 139 8 W·m － 1K － 1)，这是因为，随着温度的升高，
构成木材的分子及其内部水分子的热运动加剧，
且孔壁间的辐射能增大，导致导热系数增大; 3 )
随着含水率的增高，温度对导热系数的影响更显
著，MCave = 62. 3% 状态下，温度每升高 1 ℃，导热
系数 随 之 增 大 约 0. 067%，而 MCave = 41. 2%，
28. 3%状态下则分别增大约 0. 044%，0. 047%，可
能为温度对水分子的影响大于对空气分子的影响
所致; 4) 含水率在 28% (约为纤维饱和点)以上
831
第 4 期 赵景尧等: 有限差分法逆求木材导热系数
时，随着温度的升高，含水率对导热系数的影响程
度增大: T = 90 ℃环境下，含水率每增大 1%，导热
系数随之增大约 0. 117%，而 T = 70，50 ℃下则分
别增大约 0. 091%，0. 094% ; 低含水率范围内，温
度对含水率的影响程度则无明显作用。原因可能
为，纤维饱和点之上，温度对细胞腔内自由水分子
热运动程度的影响比对细胞壁内结合水的影响
大，致使产生本项结果。
表 1 不同温度下不同含水率时的导热系数及其标准差①
Tab． 1 Thermal conductivities at various temperatures and moisture contents
组别
Group
环境温度
T /℃
平均含水率
MCave (% )
含水率
MC(% )
基本密度
Basic density /( kg·m － 3 )
导热系数 Thermal
conductivity λ /(W·m － 1 K － 1 )
最小值 最大值 均值 标准差 均值 标准差
Min． Max． Average SD Average SD
#1 － 50 50 28. 3 25. 6 33. 0 445. 38 43. 04 0. 081 0 * 0. 010 5
#1 － 70 70 28. 3 25. 6 33. 0 445. 38 43. 04 0. 090 0 0. 008 1
#1 － 90 90 28. 3 25. 6 33. 0 445. 38 43. 04 0. 099 9 0. 011 4
#2 － 50 50 41. 2 36. 7 44. 7 380. 43 33. 45 0. 095 3 0. 011 9
#2 － 70 70 41. 2 36. 7 44. 7 380. 43 33. 45 0. 102 0 0. 010 6
#2 － 90 90 41. 2 36. 7 44. 7 380. 43 33. 45 0. 112 9 0. 010 2
#3 － 50 50 62. 3 58. 3 66. 4 421. 42 49. 33 0. 113 1 0. 013 3
#3 － 70 70 62. 3 58. 3 66. 4 421. 42 49. 33 0. 121 0 0. 010 7
#3 － 90 90 62. 3 58. 3 66. 4 421. 42 49. 33 0. 139 8＊＊ 0. 010 7
均值 Average — 43. 9 40. 2 48. 0 415. 74 41. 94 0. 106 1 0. 010 8
① λ ave = (∑
M
i = 1
∑
N
j = 1
λ i，j) /(M × N)，式中: M 为时间节点数，N 为空间节点数。M is the number of time nodes，N is the number of space nodes．
＊＊:最大平均导热系数 Maximum average thermal conductivity; * :最小平均导热系数 Minimum average thermal conductivity．
表 2 方差分析
Tab． 2 Variance analysis
误差来源
Source
自由度
Df
方差
SS
均方差
MS
F －检验
F-test
P －检验
P-test
MC 2 0. 030 72 0. 015 36 126. 942 1 0
T 2 0. 023 96 0. 011 98 99. 008 3 0
MC × T 4 0. 047 92 0. 019 88 164. 297 5 0
Error 216 0. 026 14 1. 21E － 4 — —
Total 224 0. 128 74 — — —
图 5 不同组别平均导热系数分布
Fig． 5 Average of thermal conductivities at various groups
4 结论
1) 兴安落叶松径向导热系数，尽管计算值沿厚
度方向存在一定波动性，但其平均值 0. 106 1 W·
m － 1K － 1(标准差 0. 010 8)符合实际要求，且与理论
计算值(0. 110 9，0. 125 2 W·m － 1K － 1)较为接近。
2) 含水率、温度对导热系数影响显著，前者影
响高于后者; 同时，含水率、温度共同作用对导热系
数亦存在显著影响; 导热系数随含水率的升高而增
大，随温度的升高亦增大。
3) 材性、含水率分布对导热系数的影响较大。
参 考 文 献
成俊卿 ． 1985．木材学 ． 北京: 中国林业出版社 ．
(Cheng J Q． 1985． Wood Science． Beijing． China Forestry Publishing
House． ［in Chinese］)
蔡从中，温玉锋，朱星键，等 ． 2009．木材导热系数的支持向量回归预
测 ．重庆大学学报，32(8) : 960 － 964．
(Cai C Z，Wen Y F，Zhu X J，et al． 2009． Wood thermal conductivity
prediction by using support vector regression． Journal of Chongqing
University，32(8) :960 － 964． ［in Chinese］)
陈瑞英，谢拥群，杨庆贤，等 ． 2005． 木材横纹导热系数的类比法研
究 ．林业科学，41(1) : 123 － 126．
(Chen Ｒ Y，Xie Y Q，Yang Q X，et al． 2005． Study on wood thermal
conductivity in transverse direction by analogism． Scientia Silvae
Sinicae，41(1) :123 － 126． ［in Chinese］)
高瑞堂，刘一星，李文深，等 ． 1985． 木材热学性质与温度关系的研
究 ．东北林业大学学报，13(4) : 22 － 27．
(Gao Ｒ T，Liu Y X，Li W S，et al． 1985． Study on the relationship
between the thermal properties of wood and the temperature． Journal
of Northeast Forestry University． 13(4) :22 － 27． ［in Chinese］)
侯祝强 ． 1992．木材导热系数的研究 ． 林业科学，28 (2) : 153 － 160．
931
林 业 科 学 51 卷
［in Chinese］)
(Hou Z Q． 1992． Studies on thermal conductivity of air-dried wood．
Scientia Silvae Sinicae，28(2) :153 － 160． ［in Chinese］)
林 铭，谢拥群，饶久平，等 ． 2013． 木材横纹导热系数一种新的表达
式的推导及与实测值比较 ．林业科学，49(2) : 108 － 112．
( Lin M，Xie Y Q，Ｒao J P，et al． 2013． Derivation of a new expression
for wood transverse thermal conductivity and comparison to
experimental value． Scientia Silvae Sinicae，49(2) :108 － 112． ［in
Chinese］)
王正林，龚 纯，何 倩 ． 2012．精通 MATLAB 科学计算 ． 北京: 电子
工业出版社 ．
(Wang Z L， Gong C， He Q． 2012． Master MATLAB scientific
computation． Beijing． Publishing House of Electronics Industry． ［in
Chinese］)
俞昌铭 ． 2011．多孔材料传热传质数值分析 ． 北京: 高等教育出版社
(Yu C M． 2011． Numerical analysis of heat and mass transfer for porous
materials． Beijing:Higher Education Press． ［in Chinese］)
俞自涛 ． 2005．着火前木材传热传质过程的试验和理论研究 ．杭州:浙
江大学博士学位论文 ．
(Yu Z T． 2005． Experimental and theoretical investigation on heat and
mass transfer in wood prior to fire． Hangzhou:PhD thesis of Zhejiang
University． ［in Chinese］)
俞自涛，胡亚才，田 甜，等 ． 2007． 木材横纹有效导热系数的分形模
型 ． 浙江大学学报:工学版，41(2) : 351 － 355
(Yu Z T，Hu Y C，Tian T，et al． 2007． Fractal model for predicting
effective thermal conductivity perpendicular to fibres of wood．
Journal of Zhejiang University: Engineering Science， 41 ( 2 ) :
351 － 355． ［in Chinese］)
Cai L P， Garrahan P． 2006． Inverse determination of thermal
conductivity in frozen lumber． Wood Sci Technol，40 (8 ) : 665 －
672．
Chen H T，Lin J Y，Wu C H，et al． 1996． Numerical algorithm for
estimating temperature-dependent thermal conductivity． Numerical
Heat Transfer，29(4) :509 － 522．
Dupleix A，Kusiak A，Hughes M，et al． 2013． Measuring the thermal
properties of green wood by the transient plane source ( TPS )
technique． Holzforschung，67(4) :437 － 445．
Gu H M，Zink-Sharp A． 2005． Geometric model for softwood transverse
thermal conductivity． Part I． Wood and Fiber Science，37 ( 4 ) :
699 － 711．
Gu H M，Hunt J F． 2006． Two-dimensional finite element heat transfer
model of softwood． Part II． Macro structural effects． Wood and Fiber
Science，38(4) : 599 － 608．
Hunt J F，Gu H M． 2006． Two-dimensional finite element heat transfer
model of softwood． Part I． Effective thermal conductivity． Wood and
Fiber Science，38(4) : 592 － 598．
Kreith F，Manglik Ｒ，Bohn M． 2010． Principles of heat transfer． Sixth
edition． Thomson Learning． Australia，Canada，A51 － A53．
Kollmann F F P，Cote Jr W A． 1968． Principles of Wood Science and
Technology． Vol． I． Solid Wood． Springer-Verlag，Berlin．
MacLean J D． 1941． Thermal conductivity of wood． Heating，Piping and
Air Conditioning，13:380 － 391．
Pordage L J，Langrish T A G． 1999． Simulation of the effect of air
velocity in the drying of hardwood timber． Drying Tech，17(1 /2) :
237 － 255．
Steinhagen H P，Lee H W． 1988． Enthalpy method to compute radial
heating and thawing of logs． Wood Fiber Sci，20(4) :415 － 421．
Siau J F． 1984． Transport processes in wood． Berlin，Heideberg，New
York:Spring-Verlag，245．
Sonderegger W，Hering S，Niemz P． 2011． Thermal behaviour of Norway
spruce and European beech in and between the principal anatomical
directions． Holzforschung，65(3) :369 － 375．
Vay O，Obersriebnig M，Müller U， et al． 2013． Studying thermal
conductivity of wood at cell wall level by scanning thermal
microscopy( SThM) ． Holzforschung，67(2) :155 － 159．
Yang Q X． 2001． Theoretical expression of thermal conductivity of wood．
Journal of Forestry Ｒesearch，12(1) :43 － 46．
Yeung W K， Lam T T． 1996． Second-order nite difference
approximation for inverse determination of thermal conductivity． Int J
Heat Mass Transfer，39(17) :3685 － 3693．
(责任编辑 石红青)
041