土木建筑与环境工程  2015, Vol. 37 Issue (6): 24-31   PDF    
重组竹的单轴与纯剪应力-应变关系
盛宝璐, 周爱萍, 黄东升    
南京林业大学 土木工程学院,南京 210037
收稿日期:2015-07-30
基金项目:国家自然科学基金(51378263);林业科学技术成果国家级推广项目([2015]21号)
作者简介:盛宝璐(1990-),博士生, 主要从事现代竹木结构研究,(E-mail)baolu.sheng@yahoo.com
周爱萍(通信作者),副教授,博士,(E-mail)zaping2007@163.com
摘要:重组竹是将竹丝束平行组坯、经高压胶合而成的一种生物质复合材料,是一种极具潜势的建筑结构材料。研究重组竹的基本力学性能和应力应变关系,是建立此类材料本构关系和进行重组竹结构非线性分析的基础。将重组竹理想化为横向各向同性复合材料,通过试验,给出了重组竹各主轴方向的单轴与各主平面的纯剪力学参数,建立了各种应力状态下的应力应变关系。结果表明,重组竹力学性能优于常用的结构用木材,且变异性较小。重组竹顺纹受拉强度约是顺纹受压强度的2倍;横纹受拉强度远远低于横纹受压强度;横切面内的剪切模量及强度远远低于另外两个方向,且横纹剪切强度是顺纹剪切强度的3倍。重组竹的应力应变关系和破坏模式与纤维参与受力程度密切相关。顺纹受拉时,拉应力完全由纤维承担,破坏表现为纤维的脆性拉断,强度最高,应力应变为完全线性关系;其他应力状态下,破坏均发生在基体或纤维-基体界面,若裂纹的扩展受到纤维限制,破坏呈渐进性,强度较低,应力应变曲线由早期的线性关系转入后期的非线性关系;当裂纹的扩展未受到纤维限制,破坏强度最低,应力应变呈线性关系。
关键词重组竹    应力应变关系    破坏模式    横观各向同性复合材料    
Stress-strain relationship of parallel strand bamboo under uniaxial or pure shear load
Sheng Baolu, Zhou Aiping, Huang Dongsheng    
College of Civil Engineering, Nanjing Forestry University, Nanjing 210037, P.R. China
Received: 2015-07-30
Foundation item: National Natural Science Foundation of China (No.51378263);National Promotion Program of Forestry Science and Technology Achievements (No.[2015]21)
Author brief: Sheng Baolu (1990-), PhD candidate, main research interest:mechanics performance analysis of bamboo/wood structure, (E-mail)baolu.sheng@yahoo.com.
Zhou Aiping (corresponding author), associate professor, PhD, (E-mail)zaping2007@163.com.
Abstract: Parallel Strand Bamboo (PSB) is a biocomposite composed of long narrow parallel bamboo strands which are adhesively bonded under high pressure. It has more and more attractive structural applications in building and construction engineering. It's important to well understand the stress-strain relationship to develop the constitutive law and conduct the nonlinear analysis of PSB structures. The PSB was treated as an transversely isotropic composite in the experiment and the uniaxial parameters in each main material axis and the pure shearing parameters in each main material plane were proposed, and the corresponding stress-strain relationships of each stress state were also established. The results show that compared with common used woods in construction engineering, PSB has higher strengths with less variability. Strength of tension parallel to grain is nearly as twice as that of compression parallel to grain. In perpendicular to grain direction, the strength of tension is much lower than the that of compression. Shearing in transverse-to-grain plane presents lowest modulus and strength than those shearing in other two directions. The shearing strength in perpendicular to grain is as about 3 times as that of shearing in parallel to grain direction. The stress-strain relationships and the failure modes of PSB are significantly depended on the way in which the fiber participated. In parallel to grain direction, tensile damage almost entirely contributes to the broken of fibers, which shows the highest strength and brittle behavior among all stress states. In other cases, when the expending of failure cracks are restricted by fibers, the damage presents progressive process and higher strength, and the stress-strain relationships exhibit linearity in the earliear life while turn to nonlinearity in the later life of the specimens. When damages take place in matrix or in fiber-matrix interface without fibers involved in, the material shows lower strength, and the stress-strain curves present linear and brittle behavior.
Key Words: Parallel strand bamboo    stress-strain relationship    failure mode    transversely isotropic composite    

竹材是一种可再生、可降解的天然生物质复合材料,一般3~5年就可以成材,具有硬阔叶树材的诸多优良特性,其比强度和比刚度高于钢材,是一种理想的绿色高强建筑材料。竹结构建筑的环境负荷远小于钢材、混凝土等传统建筑。但原竹材壁薄中空,几何变异性大、含糖高易虫蛀、质地不均匀、耐久性较差,不能满足现代工程结构对构件的力学性能和几何构形要求。原竹剖割成约长2 m、宽15 mm、厚3 mm的竹蔑,经80 ℃恒温烘干至含水率低于11%,再将竹篾碾压成竹丝束,同方向平行组坯并浸渍酚醛树脂胶,通过高温热压胶合,制成重组竹[1-3]。由于在工业化制造的过程中,竹材经过筛选,剔除了原竹的缺陷,且含水率低,故重组竹力学性能均匀、变异性小、强度高、耐久性好,很少出现收缩、翘曲、开裂等现象。重组竹可以加工成不同截面的型材,适用于房屋的梁、柱及跨度较大的构件[4-6]。因此,重组竹可以满足现代建筑结构对材料的力学、环保与耐久性等方面的性能要求。研究表明,重组竹结构在节能环保、工业化生产、装配式施工等方面有着传统材料不可替代的优势,以重组竹作为结构材料建造多层甚至高层建筑是极具潜势的发展方向。

重组竹纤维沿材料纵向相互平行,沿横向随机均匀分布,其力学性能具有明显的定向性,不仅纵向(顺纤维方向)与横向(垂直纤维方向)的力学性能差异很大,即使同一方向的拉、压性质亦明显不同[7],故重组竹是一种纤维定向增强复合材料[8]。试验表明,重组竹的纵向受压应力应变关系呈明显的非线性特征,这使得重组竹结构基本构件在接近或达到承载力极限状态时,荷载变形关系亦呈明显的非线性[9]。在结构设计中考虑这种非线性特性,是概率极限状态设计方法的基本要求,采用中国现行《木结构设计规范》(GB 50005-2013)[10]中基于线性原理的方法计算重组竹构件的承载力或变形误差很大。

此外,结构服役期内,材料常常处于复杂应力状态。作为复杂应力状态下纤维定向复合材料的破坏准则,常包含多个材料常数。如,Hill[11]提出的正交各向异性材料破坏准则,包涵了6个弹性常数,由于Hill准则忽略了材料拉、压性质的差异,故不适用于重组竹材料强度分析。Hoffman[12]考虑了材料拉、压力学性能的差异,提出了如式(1)的破坏准则。

$ \begin{array}{l} {C_1}{({\sigma _2}-{\sigma _3})^2} + {C_2}{({\sigma _3}-{\sigma _1})^2} + {C_3}{({\sigma _1}-{\sigma _2})^2} + \\ {C_4}{\sigma _1} + {C_5}{\sigma _2} + {C_6}{\sigma _3} + {C_7}\tau _{12}^2 + {C_9}\tau _{12}^2 = 1 \end{array} $ (1)

式中:σiτij(i, j=1, 2, 3)分别是3个主材料方向主应力和3个主平面的剪应力;Ck(k=1, 2, …, 9)是材料常数,由3个主轴方向的单轴拉、压和3个主平面的纯剪强度确定。Hoffman准则可以作为重组竹材料的屈服条件。破坏准则建立了材料的应力状态与强度之间的关系,但为了分析结构从加载至破坏全过程响应,还须了解材料在各种应力状态下的应力应变关系。可见,确定重组竹的弹性常数、建立其非线性应力应变关系,是重组竹结构基本构件非线性分析的基础。

本文将重组竹理想化为横向各向同性复合材料,通过各主轴和主平面的单轴与纯剪试验,研究重组竹3个主方向的拉、压力学性能和3个主平面的剪切性能[13-15],建立相应的应力应变关系,为进一步建立重组竹的本构关系[16]奠定基础。

1 重组竹力学性能的描述

为了定量的描述重组竹材料的力学性能,采用如图 1所示的三维笛卡尔坐标系,1轴代表材料纵向,即纤维方向;2轴、3轴分别为垂直于1轴的另外2个方向,即横向。因此,需9个弹性常数来描述重组竹刚度系数,即6个单轴拉、压模量和3个剪切模量,即顺纹、横纹和横切面内(本文将垂直于纤维的平面定义为横切面)的剪切模量。根据横观各向同性假定,材料在2轴与3轴方向的弹性常数相同,因此,重组竹的弹性常数缩减为7个,即4个拉、压弹性模量E11tE11cE22tE22c,和3个剪切模量G12G31G23。这里,E为弹性模量,其中,第一个下标代表所在平面的外法线方向,第二个下标表示弹性模量方向,上标t、c分别表示拉伸和压缩;G为剪切模量,其中,第一个下标表示剪切平面的外法线方向,第二个下标表示剪应力方向。同样,描述重组竹强度的参数也有7个,分别为3个主方向的拉、压强度fitfic (i=1, 2, 3, f2t=f3tf2c=f3c)和3个主平面的剪切强度S23S31S12。上述所有弹性常数和强度参数均可通过重组竹3个主方向的单轴拉、压和3个主平面纯剪试验得到。

图 1 重组竹3个材料主方向 Fig. 1 Main direction of PSB composite

2 重组竹基本力学性能试验

试验竹材取自浙江安吉5年生毛竹(phyllostachys),按照重组竹标准制作工艺制成试验材料,试件的密度和含水率分别为11.0 kN/m3和8.0%。每个试验类型的试件为30个,试件构形见表 1。单轴试验参照ASTM D143-09[17]、剪切试验参照ASTM D7078标准[18]进行。采用TDS-530多通道数据采集箱同时采集荷载、变形和应变,采集频率为1 Hz。各类型试验过程与现象详见表 2

表 1 试件的尺寸、形状和相关参数计算公式 Table 1 Sketch of experimental specimens and the calculation formulae for related parameters

表 2 试验过程与现象 Table 2 Process and phenomenon from the experiments

表中:ΔF11tF11tu分别为顺纹受拉的线性阶段荷载增量和极限荷载;Δε11t和Δε22t分别是顺纹与横纹受拉线性阶段的应变增量;b1t1是顺纹受拉试件有效部位的厚度与宽度,试验中通过对断口实际测量得到。ΔF22tF22tu分别为横纹受拉的线性阶段荷载增量和极限荷载;b2t2分别为横纹受拉试件有效部位的厚度和宽度。ΔF11cF11cp分别为顺纹受压的线性阶段荷载增量和极限荷载;Δε11c为顺纹受压线性阶段应变增量;w1d1为顺纹受压试件横截面尺寸。ΔF22cF22cu分别为横纹受压的线性阶段荷载增量和峰值荷载;Δε22c为横纹受压线性阶段应变增量;w2d2为横纹受压试件横截面的尺寸。ΔF21、ΔF23、ΔF12分别为2-1、2-3和1-2面内纯剪试验线性响应阶段的剪切荷载增量;F21uF23uF12u分别为2-1、2-3和1-2面内纯剪试验的峰值剪力;Δε45、Δε-45分别为在剪切试件的测试点与荷载方向成45°和-45°的应变增量;lδ为剪切试件V形切口处的截面宽度和厚度。

3 试验结果与分析
3.1 主要力学性能参数

统计分析表明,各类型试验的力学参数基本呈正态分布,图 2给出了横纹受压力学参数的正态检验结果。表 3为试验测得的力学参数统计分析结果。可以看出,重组竹顺纹拉、压模量差异很小,顺纹受拉强度约为顺纹受压强度的两倍;横纹受拉弹性模量约是横纹受压弹性模量的两倍,横纹受拉强度远远低于横纹受压强度;顺纹剪切弹性模量与横纹剪切弹性模量非常接近,顺纹剪切强度远远低于横纹剪切强度,而横切面内的剪切模量及强度远远低于另外2个方向。

图 2 横纹受压正态分布检验图 Fig. 2 Gaussian distribution test for compression in perpendicular-to-grain

表 3 试验力学参数 Table 3 The mean values and CVs of mechanical parameters of each type of experiments

表 4将重组竹与常用的结构用木材性能[19]进行了比较,可以看到,重组竹顺纹弹性模量高于大部分的常用的结构用木材,而强度远高于所有常用的结构用木材;重组竹材料力学参数的变异性比木材小。因此,重组竹的性质优于木材,是一种理想的建筑结构材料。

表 4 重组竹与结构木材的力学性能比较 Table 4 Comparing mechanical properties of PSB with that of some softwoods in structural use

3.2 应力-应变关系

根据试验结果,重组竹各应力状态下的应力应变关系大致可以分为4类。类型1为完全线性关系,包括顺纹受拉、顺纹剪切、横切面内的剪切。类型2为2段式应力应变关系,即应力小于比例极限时,应力应变关系呈线性,超过比例极限后,应力应变曲线表现为非线性强化特征,包括横纹抗压、横纹剪切。顺纹抗压应力应变关系属于类型3,为3段式应力应变关系,即应力小于比例极限时,应力应变呈线性关系;应力介于比例极限和峰值应力之间时,应力应变关系呈非线性强化特性;超过峰值后,曲线进入非线性软化段,该类型的破坏表现出良好的延性。横纹抗拉应力应变关系属于类型4,从加载到破坏,应力应变关系都表现明显的非线性强化特征。

对于类型1的应力应变关系和其他类型应力应变关系的线性部分,可以将应力表达为关于应变的一次函数,其中,一次项的斜率为材料弹性模量,常数项为0。通过对试验曲线的非线性段特征的分析,可以采用二次函数或指数函数来描述非线性应力应变关系,则类型2与类型3的应力应变关系具体表达如下

横纹受压应力应变关系

$ \sigma _{22}^{\rm{c}}(\varepsilon ) = \left\{ \begin{array}{l} E_{22}^{\rm{c}}(\varepsilon )\;\;\;\;\varepsilon \le \sigma _{22}^{{\rm{ce}}}\\ k{\gamma ^n}\;\;\;\;\;\sigma _{22}^{\rm{c}} \le \varepsilon \le \sigma _{22}^{{\rm{cu}}} \end{array} \right. $ (2)

横纹剪切应力-应变关系

$ {\tau _{23}}(\gamma ) = \left\{ \begin{array}{l} {G_{23}}\gamma \;\;\;\;\;\;\gamma \le \gamma _{23}^{{\rm{ce}}}\\ \mu {\gamma ^m}\;\;\;\;\;\gamma _{23}^{{\rm{ce}}} \le \gamma \le \gamma _{23}^{{\rm{cu}}} \end{array} \right. $ (3)

式中,材料常数可以通过$n=\frac{{{\rm{In}}(f_{22}^{{\rm{ce}}}/f_{22}^{{\rm{cu}}})}}{{{\rm{In}}(\varepsilon _{22}^{{\rm{ce}}}/\varepsilon _{22}^{{\rm{cu}}})}}$k=f22ce/(ε22ce)n$m=\frac{{{\rm{In}}(\tau _{23}^{{\rm{ce}}}/\tau _{23}^{{\rm{cu}}})}}{{{\rm{In}}(\gamma _{23}^{{\rm{ce}}}/\gamma _{23}^{{\rm{cu}}})}}$μ=τ22ce/γ22ce)m确定。这里,ε22cef22ceε22cuf22cu分别是横纹受压比例极限压应变、压应力和极限压应变、压应力;γ23ceτ23ceγ23cuτ23cu分别是横纹剪切比例极限剪应力、剪应变和极限剪应力、剪应变。

顺纹受压应力应变关系

$ \sigma _{11}^{\rm{c}}(\varepsilon ) = \left\{ \begin{array}{l} \;\;\;\;\;E_{11}^{\rm{c}}\varepsilon \;\;\;\;\;\;\;\;\;\;\;\;\;\;{\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} \varepsilon \le \varepsilon _{11}^{{\rm{ce}}}\\ {\lambda _1}{\varepsilon ^2} + {\lambda _2}\varepsilon + {\lambda _3}\;\;\;\;\varepsilon _{11}^{{\rm{ce}}} \le \varepsilon \le \varepsilon _{11}^{{\rm{cm}}}\\ {\beta _1}{\varepsilon ^2} + {\beta _2}\varepsilon + {\beta _3}\;\;\;\varepsilon _{11}^{{\rm{cm}}} \le \varepsilon \le \varepsilon _{11}^{{\rm{cu}}} \end{array} \right. $ (4)

式中:系数λiβi(i=1, 2, 3), 可以由下式确定

$ \begin{array}{l} {\lambda _1} = \frac{{f_{11}^{{\rm{ce}}}-f_{11}^{{\rm{cm}}}}}{{{{(\varepsilon _{11}^{{\rm{cm}}}-\varepsilon _{11}^{{\rm{ce}}})}^2}}}, {\lambda _2} = \frac{{2\varepsilon _{11}^{{\rm{cm}}}(f_{11}^{{\rm{cm}}}-f_{11}^{{\rm{ce}}})}}{{{{(\varepsilon _{11}^{{\rm{cm}}} - \varepsilon _{11}^{{\rm{ce}}})}^2}}}, {\lambda _3} = \\ \;\;\;\;\;\;\;\frac{{{{(\varepsilon _{11}^{{\rm{ce}}})}^2}f_{11}^{{\rm{cm}}} - 2\varepsilon _{11}^{{\rm{ce}}}\varepsilon _{11}^{{\rm{cm}}}f_{11}^{{\rm{cm}}} + {{(\varepsilon _{11}^{{\rm{cm}}})}^2}f_{11}^{{\rm{ce}}}}}{{{{(\varepsilon _{11}^{{\rm{cm}}} - \varepsilon _{11}^{{\rm{ce}}})}^2}}}\\ {\beta _1} = \frac{{f_{11}^{{\rm{cm}}} - f_{11}^{{\rm{cu}}}}}{{{{(\varepsilon _{11}^{{\rm{cu}}} - \varepsilon _{11}^{{\rm{cm}}})}^2}}}, {\beta _2} = \frac{{2\varepsilon _{11}^{{\rm{cu}}}(f_{11}^{{\rm{cu}}} - f_{11}^{{\rm{cm}}})}}{{{{(\varepsilon _{11}^{{\rm{cm}}} - \varepsilon _{11}^{{\rm{ce}}})}^2}}}, {\beta _3} = \\ \;\;\;\;\;\;\;\frac{{{{(\varepsilon _{11}^{{\rm{cm}}})}^2}f_{11}^{{\rm{cu}}} - 2\varepsilon _{11}^{{\rm{cm}}}\varepsilon _{11}^{{\rm{cu}}}f_{11}^{{\rm{cu}}} + {{(\varepsilon _{11}^{{\rm{cu}}})}^2}f_{11}^{{\rm{cm}}}}}{{{{(\varepsilon _{11}^{{\rm{cu}}} - \varepsilon _{11}^{{\rm{cm}}})}^2}}} \end{array} $

式中:ε11cef11ceε11cmf11cmε11cuf11cu分别是顺纹受压比例极限、屈服极限、极限压应变和压应力。

横纹受拉应力-应变关系

$ \sigma _{22}^{\rm{t}}(\varepsilon ) = \frac{{f_{22}^{{\rm{tu}}}-E_{22}^{\rm{t}}\varepsilon _{22}^{{\rm{tu}}}}}{{{{(\varepsilon _{22}^{{\rm{tu}}})}^2}}}{\varepsilon ^2} + E_{22}^{\rm{t}}\varepsilon $ (5)

式中:ε22tuf22tu是横纹受拉极限拉应变和拉应力。

图 3表示上述解析公式得到的曲线与试验曲线的比较。可以看出,两者吻合良好。

图 3 应力应变曲线曲线的比较 Fig. 3 Compare the calculating stress-strain curves to those obtained by experiments

4 结论

重组竹可以理想化为横向各向同性正交异性复合材料,重组竹的应力应变关系与破坏形态与纤维参与受力的程度有关。顺纹受拉作用主要由纤维承担,破坏表现为纤维的突然断裂,故表现了较高的强度与脆性,应力应变关系呈完全线性特征;顺纹、横纹抗压破坏时,裂纹出现在基体纤维界面上,但裂缝的扩展又被纤维限制,纤维在一定程度上参与受力,破坏过程呈渐进性特征,故它们的强度低于顺纹抗拉强度,应力应变关系在比例极限后表现出非线性特征;顺纹剪切、横纹剪切和横切面内的剪切破坏属于基体界面破坏,纤维几乎不参与受力,表现了较低的强度与脆性。本文采用指数函数或二次函数来描述重组竹的非线性应力应变关系,给出的解析表达式与试验结果吻合良好。

重组竹顺纹拉、压弹性模量相差无几,顺纹抗拉强度大约是横纹抗压强度的2倍,横纹抗拉强度低于横纹抗压强度;顺纹、横纹剪切弹性模量相近,大约是顺纹抗拉弹性模量的1/8;顺纹抗剪强度大大低于横纹抗剪强度,横切面内的抗剪弹性模量与强度都低于另外两个方向的剪切强度与剪切模量。

参考文献
[1] 张宏健, 李君, 叶喜. 竹条重组枋生产工艺的研究开发[J]. 建筑人造板, 1998(3) : 24–26.
Zhang H J, Li J, Ye X. Research on the fabricated technique of reformed bamboo[J]. Building Artificial Boards, 1998(3) : 24–26. (in Chinese)
[2] 李琴, 汪奎宏. 重组竹生产工艺的初步研究[J]. 人造板通讯, 2001(7) : 6–9.
Li Q, Wang K H. Research on the manufacture process of reformed bamboo[J]. China Wood-based Panels, 2001(7) : 6–9. (in Chinese)
[3] 叶良明, 姜志宏, 叶建华. 重组竹板材的研究[J]. 浙江林学院学报, 1991, 8(2) : 133–140.
Ye L M, Jiang Z H, Ye J H. Study on the laminated bamboo panel[J]. Journal of Zhejiang A & F University, 1991, 8(2) : 133–140. (in Chinese)
[4] Shigeyasu A D, Sun U T. Fracture properties of bamboo[J]. Composites: Part B, 2001, 32 : 451–459. DOI:10.1016/S1359-8368(01)00022-1
[5] Zhou A P, Huang D S, Li H T, et al. Hybrid approach to determine the mechanical properties of fibers and matrixes of bamboo[J]. Construction and building materials, 2012, 35 : 191–196. DOI:10.1016/j.conbuildmat.2012.03.011
[6] Zheng M H. Micromechanical prediction of ultimate strength of transversely isotropic fibrous composites[J]. Solids and Structures, 2011, 38 : 4147–4172.
[7] Amada S, Lchikawa Y, Munekata T, et al. Fiber texture and mechanical graded structure of bamboo[J]. Composites Part B, 1997, 28(1) : 13–20.
[8] Abdul H P S, Bhat I U H, Jawaid M, et al. Bamboo fiber reinforced biocomposites: A review[J]. Materials and Design, 2012, 42 : 353–368. DOI:10.1016/j.matdes.2012.06.015
[9] Yokozeki T, Ogasawara T, Ishikawa T. Effects of fiber nonlinear properties on the compressive strength prediction on unidirectional carbon-fiber composites[J]. Composite Science and Technology, 2005, 65 : 2140–2147. DOI:10.1016/j.compscitech.2005.05.005
[10] Canadian Wood Council. Wood design manual[S]. Ottawa, Ontario, Canada, 2010.
[11] Hill R. The mathematical theory of plasticity[M]. Lodon: Oxford University Press, 1950 .
[12] Hoffman O. The brittle strength of orthotropic materials[J]. J. Compos. Mater, 1967, 2 : 200–206.
[13] Naik N K, Kumar R S. Compressive strength of unidirectional composites: evaluation and comparison of prediction modes[J]. Composite Structures, 1999, 46 : 299–308. DOI:10.1016/S0263-8223(99)00098-7
[14] Shioya M, Takaku A. Estimation of fiber and interfacial shear strength by using a single-fiber composite[J]. Composites Science and Technology, 1995 : 33–39.
[15] Gundel D B, Waner F E. Experimental and theoretical assessment of the longitudinal tensile strength of unidirectional Sci-Fiber / Titanium-matrix[J]. Composite Science and Technology, 1997, 57 : 471–481. DOI:10.1016/S0266-3538(96)00163-7
[16] Matzenmiller A, Lubliner J, taylor R L. A constitutive model for anisotropic damage in fiber-composites[J]. Mechanics of Materials, 1995, 20 : 125–152. DOI:10.1016/0167-6636(94)00053-0
[17] American Society For Test Materials (ASTM). Standard test method for small clear specimen of timber, D143-09[S].
[18] American Society for Test Materials (ASTM). Standard test method for shear properties of composite materials by V-notched rail shear method. D7078/D7078M-12[S].
[19] United States Department of Agriculture, Forest Products Laboratory. Wood Handbook:Wood as an Engineering Material[M].Centennial Edition, 2010.