**"Process optimization of large-size bamboo bundle laminated veneer lumber (BLVL) by Box-Behnken Design,"**

*BioRes.*13(1), 1401-1412.

#### Abstract

This work focuses on optimization of the laminated lap-joint lengthening technology that is used to produce large-size bamboo bundle laminated veneer lumber (BLVL). A three-factor Box-Behnken design was developed in which lap-joint length (x1), board density (x2), and thickness of lap veneer (x3) were the three factors. Multi-objective optimization of response surface model was used to obtain 17 optimum Pareto solutions by a genetic algorithms method. The mechanical properties of BLVL predicted using the model had a strong correlation with the experimental values (R2 = 0.925 for the elastic modulus (MOE), R2 = 0.972 for the modulus of rupture (MOR), R2 = 0.973 for the shearing strength (SS)). The interaction of the x1 and x3 factors had a significant effect on MOE. The MOR and shearing SS were significantly influenced by the interaction of x2 and x3 factors. The optimum conditions for maximizing the mechanical properties of BLVL lap-joint lengthening process were established at x1 = 16.10 mm, x2 = 1.01 g/cm3, and x3 = 7.00 mm. A large-size of BLVL with a length of 14.1 m was produced with the above conditions. Strong mechanical properties and dimensional stability were observed.

Download PDF

#### Full Article

RETRACTED: June 11, 2020. The authors retract their study due to errors in data analysis that change the conclusions of this paper.

**Process Optimization of Large-Size Bamboo Bundle Laminated Veneer Lumber (BLVL) by Box-Behnken Design**

Xiaoyi Niu, Jiuyin Pang, Hanzhong Cai, Shan Li, Lei Le, and Junhua Wu*

This work focuses on optimization of the laminated lap-joint lengthening technology that is used to produce large-size bamboo bundle laminated veneer lumber (BLVL). A three-factor Box-Behnken design was developed in which lap-joint length (*x*_{1}), board density (*x*_{2}), and thickness of lap veneer (*x*_{3}) were the three factors. Multi-objective optimization of response surface model was used to obtain 17 optimum Pareto solutions by a genetic algorithms method. The mechanical properties of BLVL predicted using the model had a strong correlation with the experimental values (R^{2 }= 0.925 for the elastic modulus (MOE), R^{2 }= 0.972 for the modulus of rupture (MOR), R^{2 }= 0.973 for the shearing strength (SS)). The interaction of the *x*_{1} and *x*_{3} factors had a significant effect on MOE. The MOR and shearing SS were significantly influenced by the interaction of *x*_{2} and *x*_{3} factors. The optimum conditions for maximizing the mechanical properties of BLVL lap-joint lengthening process were established at *x*_{1 }= 16.10 mm, *x*_{2 }= 1.01 g/cm^{3}, and *x*_{3 }= 7.00 mm. A large-size of BLVL with a length of 14.1 m was produced with the above conditions. Strong mechanical properties and dimensional stability were observed.

*Keywords: Bamboo bundles; Lap-joint; Laminated veneer lumber (LVL); Box-Behnken design; **Response surface methodology; Mechanical properties*

*Contact information: Key Laboratory of Wood Material Science and Engineering of Jilin Province, Beihua University, Jilin, Jilin Province, China,132013; *Corresponding author: 54086442@qq.com*

**INTRODUCTION**

The development of energy efficient buildings made from green building materials has the potential to initiate a paradigm shift leading to a huge demand for environmentally friendly construction materials. It follows that large-size bamboo-based engineered products will become a valid alternative to traditional construction materials due to their abundance as a raw material, good structural integrity, aesthetic pleasance, effective seismic resistance, and ease of manufacturing (Nugroho and Ando 2001; Lugt *et al*. 2005; Huang *et al*. 2015; Yu *et al*. 2015). A variety of bamboo-based products, including bamboo-strip laminates, mat and curtain ply bamboo, glue-laminated bamboo, and bamboo scrimber have been developed and successfully applied in buildings, bridges, floors, transportation, and other fields (Hu and Pizzi 2013; Sharma *et al*. 2015; Chen *et al*. 2016; Yanjun *et al*. 2016). However, the dimensions of artificial bamboo boards developed to this point have been small with sizes of 1.22 m to 2.44 m, which severely hinders further applications in larger-span structural applications. Components such as double beam bamboo bundle laminated veneer lumber have limited use because of the scarcity of continuous hot pressing machinery (Li *et al*. 2014; Chen *et al*. 2016).

Bamboo bundle laminated veneer lumber (BLVL) is manufactured by arranging broomed bamboo bundle fibers impregnated with phenol-formaldehyde resin (PF) in a parallel lay-up, which was then hot pressed. Bamboo bundle laminated veneer lumber as a potential construction material has a favorable uniform density, dimensional stability, and stability of mechanical performance (Chen *et al*. 2013, 2014). The effects of scarf-joint and finger-joint forms on the mechanical properties of large-scale BLVL have been compared (Yeh and Lin 2012; Deng *et al*. 2014). Zhang *et al*. (2014) evaluated the effects of the four veneer-joint forms (*i.e.*, button joint, lap joint, toe joint, and tape joint) on the mechanical properties of BLVL and found that the best veneer-joint form was the lap joint laminate. Recently, large span BLVL (length > 2.44 m) was successfully manufactured by combining the prepress densification process and the intermittent hot-press method (Sharma *et al*. 2015). According to the intermittent hot-press process, research is required on the parallel lap-joint and assembly technology of bamboo bundle veneer in order to allow for effective production of large-span BLVL. Additionally, broomed bamboo bundle fibrosis veneer is a loosely laminated reticulate sheet and is much harder than wood, which makes it inconvenient to mill. Based on the lap-joint method used in laminated veneer lumber (LVL) (Zhang *et al*. 2014), loose bamboo bundle fibers are laminated in a crisscross overlap with each other at the lap-joint. The local density and stress distribution at the lap-joint have a large variation caused by the technological parameters of the process such as the thickness of lap veneer, lap-joint length, and board density. The variance of density and stress distribution is of key importance for the mechanical behavior and dimensional stability of large-span BLVL. Therefore, the lap-joint process and lay-up method of bamboo bundle veneer would have a great influence on the mechanical properties of board, especially at the lap-joint.

Once a certain mechanical index of board needs to be satisfied for the engineering design requirements, manufacturing processes parameters should be adjusted accordingly. Typical methods are based on conducting a large amount of destructive testing, which is repetitive, labor-intensive, inefficient, blind, and a waste of resources. If the inner mathematical relationships between process parameters and the mechanical properties of a board could be developed using numerical simulation, then a predetermined regard for the lap-joint assembly pattern could effectively be obtained with the optimization process parameters. Based on this, the parameters and model are then corrected and reconfirmed by the re-experimental data, which improves the work efficiency and saves resources. Response surface methodology (RSM) is an accurate and effective tool which is a kind of statistical and mathematical method used for optimizing experimental processes (Bezerra *et al.* 2008; Jiang *et al*. 2013; Subramonian *et al*. 2015). By designing a series of experiments for adequate and reliable data for the response of variables of interest, a mathematical model of the second order response surfaces with the best fittings is developed. The best finite-parameter solutions with a maximum or minimum value of interest and the interactive effects of process parameters in 3-D contour plots could be observed (Aslan and Cebeci 2007; Lai *et al. *2013; Yuan *et al*. 2015).

There are few studies available in literature about the optimization of lap-joint assembly process for large-span bamboo bundle laminated veneer lumber (BLVL) using RSM. Therefore, the focus of this study is to optimize the laminated lap-joint lengthening technology using RSM on a three-level three-factorial Box-Behnken experimental design. The Functional relationship between the three factors (x_{1}: lap-joint length, x_{2}: board density, x_{3}: veneer thickness) on the mechanical properties of lap-joint BLVL were investigated.

**EXPERIMENTAL**

**Manufacture Process of BLVL**

Cizhu bamboo (*Neosinocalamus affinis*) bundle veneers were immersed in phenol formaldehyde (PF) resin (solid content of 12%) for 5 to 7 min and then dried to a moisture content of 10% to 12% under an ambient environment. The resin soaked veneers were layered symmetrically along their grain direction with the outer layer of bamboo bundle facing upward. The overlap joint in each bamboo bundle layer was set on a 1/6 to 1/2 position along the length of the board. The lap-joints in adjacent layers were designed with horizontal separation. A pressure, temperature, and time of 3.5 MPa, 150 °C, and 30 min, respectively, were applied to the BLVL during the hot press. The dimensions of BLVL were 300mm × 150 mm × 12.5 mm (length*width* thickness).

**Fig. 1. **Production of bamboo bundle veneer and lap-joint allocation on BLVL

**Mechanical Properties**

The BLVL were tested for modulus of elasticity (MOE), modulus of rupture (MOR), and shearing strength (SS) under perpendicular loading to the glue line. Seventeen coded groups were used for each test with seven replicates per group (Table 1) for a total of 357 samples prepared for testing. The testing procedures to determine the MOE, MOR, and SS were in accordance with the ASTM standards D1037 (2006), D3500 type A (2003), and D2344 (2006), respectively.

**Response Surface Methodology and Box-Behnken Experimental Design**

Response variables reflected the results of interest, independent variables were those which affected the response variables, and variables *m* and *k* were the number of response variables and independent variables, respectively. Equation 1 is an empirical expression for response variables and independent variables.

The objective was to find a functional relationship between the independent variables (*x*) and response variables (*y*), as well as what independent variable (*x*) would optimize the response variables (*y*). A second-order polynomial model was utilized according to the following equation,

where were the regression coefficients of constant, linear, quadratic, and interactions terms, respectively.

The independent variables (*x*_{1}: lap-joint length (mm), *x*_{2}: board density(g/cm^{3}), *x*_{3}: veneer thickness (mm)) were at three levels (*x*_{1}: 10 mm (-1), 20 mm (0), 30 mm (1);* x*_{2}: 0.8 g/cm^{3 }(-1), 0.95 g/cm^{3}(0), 1.1 g/cm^{3 }(1);* x*_{3}: 4.0 mm (-1), 5.5 mm (0), 7.0 mm (1)). For the three-level three-factorial Box-Behnken experimental design (BBD), a total of 17 experimental runs were made in random order (Table 1). Response variables of (*y*_{1}: MOE,* y*_{2}: MOR,* y*_{3}: SS) were expressed by means of seven measurements. The specific model for three-level three-factorial is as follows,

where represents the predicted response of interest; was the constant regression coefficients; were the linear regression coefficients; were the quadratic coefficients, and were the interaction coefficients. These coefficients were determined in the second-order model by computer simulation programming that applied the least square method. The data analyses were undertaken using the Design-Expert V.8.0.6 (Stat-Ease, Inc., Minneapolis, MN, USA).

**Table 1.** Experiment Design by BBD, Experimental and Predicted Values of Mechanical Properties of BLVL

Exp., Experimental; Pred., Predicted

**RESULTS AND DISCUSSION**

**Results of Response Surface Methodology Experiments**

The effects of three independent variables(*x*_{1}: lap-joint length, *x*_{2}: board density, and *x*_{3}: veneer thickness) with the three-coded level on the mechanical properties of BLVL (*y*_{1}: MOE,* y*_{2}: MOR, and*y*_{3}: SS) using BBD are reported in Table 1. The predicted values for the response variables matched well with the experimental data, indicating a strong correlation (Table 2). The low coefficient of variation (CV) values indicated good stability in the mechanical behavior of BLVL, which was adequate and suitable for the mathematical model as a function Eq. 3.

**Table 2.** Model Adequacy Indicators for Each Modeled Response of BLVL

After fitting the mechanical properties data to the experimental results using Eq. 3, the following quadratic response functions representing *y*_{1},* y*_{2}, and* y*_{3} were be expressed as a function of lap-joint length (*x*_{1}), board density (*x*_{2}), and veneer thickness (*x*_{3}), respectively. Adj-R^{2 }represents the adjusted value of R^{2} after removing the effect of the number of independent factors on the response variables. The coefficients in the model equations are rounded.

For the *y*_{1} MOE model equation, (4)

For *y*_{2} MOR model equation: (5)

For y3 SS model equation:

(6)

Analysis of variance for the response surface second-order model of the MOE, MOR, and SS show the values of “Prob > F” less than 0.05 indicated that the model had a strong significance (Tables 3, 4, and 5).

**Table 3. **Analysis of Variance for MOE Regression Equation

**Table 4. **Analysis of Variance for MOR Regression Equation

**Table 5. **Analysis of Variance for Shearing Strength (SS) Regression Equation

The model F-values for MOE, MOR, and SS regression equation indicated that the models were all significant. In addition, the lack of fit values of “Prob > F” were all more than 0.05, which implies the models fitted the MOE, MOR, and SS well.

**Effect of Lap-Joint Variables on the Mechanical Properties of BLVL**

A serials 3-D response surface contour plots were drawn to clarify the effect of the variables on the mechanical properties of BLVL (Figs. 2 through 4). The interactive effect of independent variables can be observed directly from the 3-D response surface contour plots. The contour plots’ oval shape (at the bottom of 3-D response surface plots) indicated the interactive effect between the response variables was significant. The contour plot representing the interaction between lap-joint length (*x*_{1}) and board density (*x*_{2}) was flat which affected the MOE of the BLVL (Fig. 2). Similarly, the elliptical shape of the contour plots in Figs. 3 and 4 indicated the interactive effect of board density (*x*_{2}) and veneer thickness (*x*_{3}) on MOR, while the board density (*x*_{2}) and veneer thickness (*x*_{3}) had a larger effect on the SS. The round contour plots in Figs. 2(a) and 4(b) and the line contour plot of Figure 3(b) showed there was a minimum interactive effect between the lap-joint length (*x*_{1}) and board density (*x*_{2}) on the MOE, lap-joint length (*x*_{1}) and veneer thickness (*x*_{3}) on the MOR, lap-joint length (*x*_{1}) and veneer thickness (*x*_{3}) on the SS, respectively.

From the variance of the MOE, MOR, and SS that the results were verified to their strong Prob > F values (Tables 3 to 5). The Prob>F values for MOE model *x*_{1}*x*_{3,} MOR model *x*_{2}*x*_{ 3,} and SS model *x*_{2}*x*_{3}were all lower than 0.05 which indicates that the results were significant at the 95% confidence level. The Prob>F values of MOE model *x*_{ 1}*x*_{ 2,}, MOR model *x*_{ 1}*x*_{ 3}, and SS model *x*_{1}*x*_{ 3}were confirmed to have no significant effect on the BLVL. It was observed from the contour plots that the BLVL had stronger mechanical properties when the board was manufactured with a middle level of lap-joint length, board density, and upper level of the veneer thickness. An appropriate lap-joint length and board density was beneficial for strong interfacial bonding between the layers of the BLVL. Based on the manufacturing conditions, bamboo bundle veneer thickness had a positive contribution to the interactive influence on mechanical properties. The main effect of these positive contributions was namely the larger the veneer thickness, the better MOR and SS obtained.

**(a) **

**(b)**

**Fig. 2.** Response surfaces plots showing the effect of lap-joint length (*x*_{1}) and (a) board density (*x*_{2}) or (b) veneer thickness (*x*_{3}) on MOE

**Fig. 3.** Response surfaces plots showing the effect of veneer thickness (*x*_{3}) and (a) board density (*x*_{2}) or (b) lap-joint length (*x*_{1}) on MOR

**Fig. 4.** Response surfaces plots showing the effect of veneer thickness (*x*_{3}) and (a) board density (*x*_{2}) or (b) lap-joint length (*x*_{1}) on SS

**Optimization of Lap-Joint Process and Its Application in Large-Span BLVL**

Multi-objective optimization was based on the response surface model by using a genetic algorithm (GA) due to its common used to generate high-quality solutions to optimization. In a GA, a population of candidate solutions of lap-joint parameters to an optimization problem is evolved toward better solutions.

**Table 6. **The 17 Pareto Solutions for Multi-Objective Optimization

*indicates selected design

Optimal parameters for obtaining the maximum object value set were determined to be control population of 30, aberration ratio of 0.1, number of iterations of 80, and crossover probability of 0.2. The 180 *E*_{b }level superior products of bamboo scrimber standard required the lower limit values of BLVL to be 18 GPa for MOE, 160 MPa for MOR, and 12 MPa for SS. The upper limit values for the 180E_{b} were 26 GPa for MOE, 180 MPa for MOR, and 19 MPa for SS. The optimized global Pareto solutions for group 17 were obtained using iterative computation with the initial values of *x*_{1},* x*_{2}, and*x*_{3} equal to 17 mm, 1.06 g/cm^{3}, and 6.0 mm (Table 6). The initial parameters of the mechanical properties of the manufactured BLVL were a MOE of 16.4 GPa, MOR of 175.6 MPa, and SS of 14.7 MPa.

The 17 groups of Pareto solutions were higher than the special requirements of the 180 *E*_{b }level superior products of bamboo scrimber standard (Table 6). In addition, all sets of 17 Pareto solutions had increased at different degrees when compared with the experimental results of BLVL processed with initial parameters. The first group selected by the computer system had the best solutions with a lap-joint length of 16.1 mm, board density of 1.01 g/cm^{3}, and veneer thickness of 7 mm, which resulted in strong response values. The desirability factor for group 1 was 0.852, which indicates that the experimental design possessed a strong feasibility (Fig. 8). Mechanical tests and aging treatment experimental results showed that the long-size BLVL possessed good strength, stiffness, and durability (Chen *et al*. 2016). The BLVL with a length of 14.1 m was developed based on the laminated lap-joint lengthening technology (Fig. 5).

**Fig. 5. **Manufacturing process of large-span BLVL from the length of 1 m to 14.1 m

**CONCLUSIONS**

- The 3-D response surface contour plots were developed in order to better understand the interactive effect of the independent variables on the mechanical properties of BLVL. A significant interactive effect between lap-joint length (
*x*_{1}) and veneer thickness (*x*_{3}) on elastic modulus (MOE) was observed. Similarly, modulus of rupture (MOR) and shearing strength (SS) were significantly affected by the board density (*x*_{2}) and veneer thickness (*x*_{3}). - Seventeen optimum Pareto solutions for determining the maximum mechanical properties of BLVL were obtained using a genetic algorithms method. This method determined that the best process combination was a lap-joint length (
*x*) of 16.10 mm, board density (_{1}*x*) of 1.01 g/cm_{2}^{3}, and veneer thickness (*x*) of 7.00 mm._{3} - The efficiency of Box-Behnken design and response surface methodology for modeling was able to optimize the fabrication of engineered bamboo-based materials, in an economical way that maximized the process parameters in the least amount of experiments as possible.

**ACKNOWLEDGEMENTS**

The authors are grateful for the financial support of the Key Technologies R&D Programme in agricultural field for the scientific and technological development of Jilin Province Funds for the Research on the super low quantity VOC aqueous isocynaote emission of furniture timber(20160204028NY). The constructive comments from the anonymous reviewers were also greatly appreciated.

**REFERENCES CITED**

ASTM D1037-12 (2006). “Standard test method for evaluating properties of wood-base fiber and particle panel materials,” ASTM International, West Conshohocken, PA.

ASTM D3500-14 (2003). “Standard test method for structural panels in tension,” ASTM International, West Conshohocken, PA.

ASTM D2344-16 (2006). “Standard test method for short-beam strength of polymer matrix composite materials and their laminates,” ASTM International, West Conshohocken, PA.

Aslan, N., and Cebeci, Y. (2007). “Application of Box-Behnken design and response surface methodology for modeling of some Turkish coals,” *Fuel* 86(1), 90-97.DOI:10.1016/j.fuel.2006.06.010

Bezerra, M. A., Santelli, R. E., Oliveira, E. P., Villar, L. S., and Escaleira, L. A. (2008). “Response surface methodology (RSM) as a tool for optimization in analytical chemistry,” *Talanta* 76(5), 965-977.DOI:10. 1016/j.talanta.2008.05.019

Chen, F., Jiang, Z., Deng, J., Wang, G., Zhang, D., Zhao, Q., Cai, L., Shi, S.Q. (2013)“Evaluation of the uniformity of density and mechanical properties of bamboo-bundle laminated veneer lumber (BLVL),” *BioResources* 9(1), 554-565. DOI:10.15376/biores.9.1.554-565

Chen, F. M., Deng, J. C., Cheng, H., Li, H. D., Jiang, Z. H., Wang, G., Zhao, Q. C., and Shi, S. Q. (2014). “Impact properties of bamboo bundle laminated veneer lumber by preprocessing densification technology,” *Journal of Wood Science* 60(6), 421-427. DOI: 10.1007/s10086-014-1424-0

Chen, F., Jiang, Z., Wang, G., Li, H. D., Smith, L. M., and Shi, S. Q. (2016). “The bending properties of bamboo bundle laminated veneer lumber (BLVL) double beams,” *Construction and Building Materials* 119, 145-151. DOI: 10.1016/ j.conbuildmat.2016.03.114

Chen, F. M., Deng, J. C., Li, X. J., Wang, G., Smith, L. M., and Shi, S. Q. (2016). “Effect of laminated structure design on the mechanical properties of bamboo-wood hybrid laminated veneer lumber,” *European Journal of Wood and Wood Products*8 3, 95-101. DOI: 10.1007/s00107-016-1080-8.

Deng, J., Li, H., Zhang, D., Cheng, F., Wang, G., and Cheng, H. (2014). “The effect of joint form and parameter values on mechanical properties of bamboo-bundle laminated veneer lumber (BLVL),” *BioResources* 9(4), 6765-6777. DOI: 10.15376/ biores.9.4.6765-6777

Hu, J. B., and Pizzi, A. (2013). “Wood-bamboo-wood laminated composite lumber jointed by linear vibration-friction welding,” *European Journal of Wood and Wood Products* 71(5), 683-686. DOI: 10.1007/s00107-013-0714-3

Huang, D., Bian, Y., Zhou, A., and Sheng, B. (2015). “Experimental study on stress-strain relationships and failure mechanisms of parallel strand bamboo made from phyllostachys,” *Construction and Building Materials* 77, 130-138. DOI: 1 0.1016/ j.conbuildmat.2014.12.012

Jiang, S., Cai, W., and Xu, B. (2013). “Food quality improvement of soy milk made from short-time germinated soybeans,” *Foods* 2(2), 198-212.DOI: 10.3390/foods2020198.

Lai, J., Xin, C., Zhao, Y., Feng, B., He, C., Dong, Y., Fang, Y., and Wei, S. (2013). “Optimization of ultrasonic assisted extraction of antioxidants from black soybean (*Glycine max* var.) sprouts using response surface methodology,” *Molecules* 18(1), 1101-1110. DOI: 10.3390/molecules18011101

Li, H., Chen, F., Cheng, H., Deng, J., Wang, G., and Sun, F. (2014). “Large-span bamboo fiber-based composites, Part I: A prediction model based on the Lucas-Washburn equation describing the resin content of bamboo fiber impregnated with different PVAC/PF concentrations,” *BioResources *9(4), 6408-6419. DOI: 10.15376/ biores.9.4.6408-6419

Lugt, P. V. D., Dobbelsteen, A. V. D., and Janssen, J. J. A. (2005). “An environmental, economic and practical assessment of bamboo as a building material for supporting structures,” *Construction and Building Materials* 20(9), 648-656. DOI: 10.1016/j.conbuildmat2005.02.023

Nugroho, N. and Ando, N. (2001). “Development of structural composite products made from bamboo II: Fundamental properties of laminated bamboo lumber,” *Journal of Wood Science *47(3), 237-242. DOI: 10.1007/BF01171228

Sharma, B., Gatόo, A., Bock, M., and Ramage, M. (2015). “Engineered bamboo for structural applications,” *Construction and Building Materials* 81, 66-73. DOI: 10.1016/j.conbuildmat.2015.01.077

Subramonian, W., Wu, T. Y., and Chai, S. P. (2015). “An application of response surface methodology for optimizing coagulation process of raw industrial effluent using *Cassia obtusifolia* seed gum together with alum,” *Industrial Crops and Products* 70, 107-115.DOI:10.1016/j.indcrop.2015.02.026

Yanjun, L., Bin, X., Zhang, Q., and Shen Xue, J. (2016). “Present situation and the countermeasure analysis of bamboo timber processing industry in China,” *Journal of Forestry Engineering *1(1), 2-7. DOI: 10.13360/j.issn.2096

Yeh, M. C. and Lin, Y. L. (2012). “Finger joint performance of structural laminated bamboo member,” *Journal of Wood Science* 58(2), 120-127. DOI: 10.1007/s10086-011-1233-7

Yu, Y., Zhu, R., Wu, B., and Yu, W. J. (2015). “Fabrication, material properties, and application of bamboo scrimber,” *Wood Science and Technology* 49(1), 83-98. DOI: 10.1007/s00226-014-0683-7

Yuan, J., Huang, J., Wu, G., Tong, J., Xie, G., Duan, J., and Qin, M. (2015). “Multiple responses optimization of ultrasonic-assisted extraction by response surface methodology (RSM) for rapid analysis of bioactive compounds in the flower head of *Chrysanthemum morifmolium* Ramat,” *Industrial Crops and Products* 74, 192-199. DOI: 10.1016/j.indcrop.2015.04.057

Zhang, D., Wang, G. and Ren, W. (2014). “Effect of different veneer-joint forms and allocations on mechanical properties of bamboo-bundle laminated veneer lumber” *BioResources* 9(2), 2689-2695. DOI: 10.15376/biores.9.2.2689-2695

Article submitted: October 11, 2017; Peer review completed: December 29, 2017; Revised version received: December 30, 2017; Accepted: January 3, 2017; Published: January 8, 2018.

DOI: 10.15376/biores.13.1.1401-1412