Articles

Modeling and parameters optimization of biocomposite using box-Behnken response surface methodology

M. Mittal *
Central University of Haryana, India
K. Phutela
Ch. Bansilal Govt. Polytechnic, India

Modeling and parameters optimization of biocomposite using box-Behnken response surface methodology

Journal of applied research and technology, vol. 21, no. 6, pp. 991-1018, 2023

Universidad Nacional Autónoma de México, Instituto de Ciencias Aplicadas y Tecnología

Received: 24 June 2022

Accepted: 22 September 2022

Published: 31 December 2023

Abstract: The primary objective of this work includes modeling and optimization of the mechanical properties of natural fiber biocomposites using three-factor, three-level Box-Behnken design (BBD). In this context, the effect of three independent performance parameters; pineapple leaf fiber (PALF) content, fiber length, and polyethylene-grafted-maleic anhydride (MAPE) compatibilizer load have been investigated on the mechanical properties of PALF/HDPE/MAPE biocomposite. The sequential model sum of squares, lack of fit, and normal probability plots showed a good agreement in between the experimental results and those predicted by mathematical models (95% confidence level). The optimization results obtained in Design-Expert software revealed that the most optimal value of tensile strength, tensile modulus, flexural strength, flexural modulus, and impact strength as 32.35 MPa, 1475 MPa, 49.21 MPa, 1659.04 MPa, and 58.24 J/m respectively, at fiber length of 13.67 mm, PALF content of 16.84 wt.%, and MAPE load of 2.95 wt.%. To verify the mathematical models, validation tests were also performed which showed that the response surface methodology (RSM) based BBD and ANOVA tools are adequate for analytically evaluating the performance of biocomposites.

Keywords: Box-Behnken design (BBD), response surface methodology (RSM), pineapple leaf fiber (PALF), high-density polyethylene (HDPE), mechanical properties, polyethylene-grafted-maleic anhydride (MAPE).

1. Introduction

Environmental and economic concerns are stimulating research in the development of eco-friendly composites for several engineering applications. In this regard, researchers put their attention towards the utilization of low-cost lignocellulosic fibers, such as bamboo, jute, sisal, PALF, coir, ramie, etc. as reinforcements in petroleum-based thermoplastics. Owing to their exceptional features such as lightweight, renewability, biodegradability, easy availability, high specific strength & stiffness, good acoustic, thermal, and electrical insulation make them as promising candidates compared to inorganic fibers (Ali et al., 2017; Ali et al., 2018; Basu et al., 2017; Bodirlau et al., 2013; Holt et al., 2014; Khalil et al., 2007; Latif et al., 2015; Li et al., 2016; Qaiss et al., 2014; Salasinska et al., 2016; Shi et al., 2013).

Among all the natural cellulosic fibers, the PALF has significant potential to be used as a reinforcing agent in fiber or powder form. This is because of the presence of high cellulose content (70-85%), low-density, and low microfibrillar angle. Moreover, PALF is a waste product of pineapple cultivation and therefore relatively inexpensive which can be employed for industrial purposes (Arib et al., 2004; Mishra et al., 2001; Pavithran et al., 1987). The fibers derived from pineapple leaves have excellent characteristics for use in automotive, building & construction, packaging, and general engineering applications. According to the statistical database (2017) of ‘Food and Agricultural Organization’, the total production of pineapple fruit in the entire world is 25.8 million tons (FAOSTAT, 2017).

Studies showed that the PALF could be used effectively for polypropylene (PP), polyethylene (PE), thermoplastic starch, and rubber reinforcement (Kengkhetkit & Amornsakchai, 2012). In addition, its carbon footprint is much lower than that of other cultivated natural and synthetic fibers (Kengkhetkit & Amornsakchai, 2014). Although cellulosic fibers have number of advantageous features, they exhibit some problems such as hydrophilic nature, limited thermal stability, and poor dispersion characteristics within the non-polar thermoplastic matrix (Bledzki et al., 1998; Cantero et al., 2003; Kazayawoko et al., 1999; Raj et al.,1989). These limitations can be remedied by several methods like grafting functional moieties onto the fibers, introduction of coupling agents, or pre- treatment of fibers with suitable agents. There have been numerous studies conducted for improving the interfacial adhesion between natural fibers and the matrix.

Table 1 reported that the chemical treatment has profound influence on the mechanical properties of fiber-reinforced composites (FRCs). It was observed that the use of a coupling agent containing maleic anhydride, polyethylene-grafted-maleic anhydride (MAPE) has not been addressed significantly. Therefore, in the present communication, the effect of maleic anhydride coupling agent on mechanical properties of PALF/HDPE/MAPE biocomposite was investigated. A review of the existing literature (Table 2) revealed that the fiber length and loading are the other critical parameters that can affect the thermo-mechanical properties of biocomposites. It was observed from literature review that most of the work includes one parameter effect at a time while keeping other at a fixed level. The interaction effect of various parameters has not been discussed so far. Therefore, to solve this problem, the design of experiment (DOE) method was used (Ashenai Ghasemi et al., 2016; Mhalla et al., 2017; Rostamiyan et al., 2014; Rostamiyan, Fereidoon, Mashhadzadeh et al., 2015; Rostamiyan, Fereidoon, Rezaeiashtiyani et al., 2015; Rostamiyan, Fereidoon, Nakhaei et al., 2017; Subasinghe et al., 2016; Ghasemi et al., 2016).

Table 1
Previous research work on chemical treatment of fibers.
Previous research work on chemical treatment of fibers.

Table 2
Literature studies on the effect of fiber length and content on properties of biocomposites.
Literature studies on the effect of fiber length and content on properties of biocomposites.

Several optimization techniques such as RSM, Taguchi method, full factorial, fractional factorial, ANN, fuzzy logic, and GA are available which have significant potential to optimize the performance parameters (Mohamed et al., 2015). Furthermore, these techniques have significant potential to develop models that can predict and established the relationship between different inputs and response variables. RSM is a promising analytical tool to determine the significance of interactions and square terms of parameters, 3D response surface generation, and optimize the parameters. A three-level BBD is economical and popular in industrial research for modeling and optimizing the parameters to satisfy the defined desirable response variables. Therefore, in this study, the Box-Behnken design (BBD) which is a subset of RSM was employed as a DOE method for optimizing the mechanical properties of HDPE/PALF/MAPE biocomposite. The individual and simultaneous effects of the fiber length, PALF loading, and MAPE compatibilizer content on tensile, flexural, and impact properties were investigated.

2. Experimental

2.1. Materials

High-density polyethylene [HDPE, density=950 kg/m3, melt flow index (MFI)=12.3 g/10 min] and the coupling agent, polyethylene maleic anhydride (MAPE, Mw 26000) was obtained from M/s Solvay Chemicals Ltd. Kerala, India. The pineapple leaf fiber (PALF) with diameter ranging in between 60 to 100 µm was collected from M/s Go Green products, Chennai, India. The collected fibers were thoroughly washed and then dried in an oven at 60oC to a final moisture content of 5%. The elemental composition and SEM micrograph of the studied PALF is shown in Table 3 and Figure 1, respectively.

Table 3
Elemental composition of PALF.
Elemental composition of PALF.

Surface morphology of untreated PALF.
Figure 1
Surface morphology of untreated PALF.

2.2. Fabrication of biocomposite samples

Biocomposite samples of each formulation (Table 4) were obtained through melt-blending in a twin-screw extruder (Brabender plasticorder, Germany), which has a mixing chamber of 69 cm3 volumetric capacity. To improve the interfacial bonding strength between fiber and matrix resin, a pre-impregnation technique was utilized. The obtained pre-pregs were cut into 6 mm long pellets and further blended with pure HDPE to obtain the desired composition. The mixing was conducted at 190oC and 35 RPM for 10 min. The pre-mixed pellets were molded in TS-270 injection molding machine (Windsor, India) under the temperature setting of injection hopper and nozzle at 170oC and 180oC, respectively. Specimens were conditioned in a laboratory atmosphere of 23±5oC and 50% RH prior to testing.

Table 4
Biocomposite samples according to box-Behnken experimental design.
Biocomposite samples according to box-Behnken experimental design.

2.3. Mechanical tests

Tensile and flexural (three-point bending) tests were conducted on a calibrated Tabletop Tinius Olsen Horizon H50KS, a Universal testing machine as per ASTM D638 and D790 standard, respectively. A digital pendulum impact tester was employed to determine impact strength. Four specimens for each composite sample were evaluated and their average values were reported.

2.4. Response surface method

RSM is a mathematical and statistical technique that can assess the effect of individual parameters and the interaction of parameters on response variables. It has significant potential in modeling and analysis of problems with a small number of experimental data points (Myres et al., 2016). In RSM, a three-level Box Behnken design (BBD) was used to study the effect of linear, quadratic, cubic, and cross product models of three performance parameters and to develop an experimental design matrix. The design was built in Design-Expert software (version 6.0.8) with consideration of three critical variables i.e., fiber length (A), fiber loading (B), and MAPE content (C). The range and levels of the selected parameters are shown in Table 5. In this work, the second order polynomial equation was used for fitting the experimental data and find out the relevant model terms.

y = b o + Σ b i x i + Σ b i i x i 2 + Σ b i j x i x j (1)

Table 5
Actual and coded experimental variables in Box-Behnken design.
Actual and coded experimental variables in Box-Behnken design.

Where y represents the predicted response, bo , regression equation constant; bi , linear coefficient; bii , square coefficient of each parameter; and bij , first order interaction coefficient. To estimate the desired responses with reliable measurements, a total number of 17 experiments were conducted with 12 factorial points. Furthermore, the experimental sequence was randomized to minimize the effect of uncontrollable parameters. The regression (R2), adjusted (R2), predicted (R2), pure error sum of squares, lack of fit, response plots, and adequate precision were used in the determination of the robustness of the developed model.

3. Results and discussion

3.1. ANOVA and regression models

The experimental test results were used to generate the mathematical models for each response variable. Among different mathematical models (linear, quadratic, cubic and 2FI), the quadratic was selected according to three different tests - the sequential model sum of squares, adequacy, and lack-of-fit (Table 6). The sum of square and mean square in ANOVA results (Tables 7-11) also revealed the adequacy of quadratic model for fitting the experimental data. The second-order equations (ii-vi) express the overall predictive model in terms of variables:

T e n s i l e   s t r e n g t h = 32.19 - 0.15 * A + 2.86 * B + 0.96 * C - 2.88 * A 2 - 6.62 * B 2 - 1.63 * C 2 + 0.24 * A * B + 0.097 * A * C + 0.62 * B * C (ii)

T e n s i l e   m o d u l u s = 1475.93 - 31.27 * A + 35.60 * B + 30.64 * C - 23.18 * A 2 - 58.43 * B 2 - 37.96 * C 2 + 22.17 * A * B + 15.76 * A * C + 36.35 * B * C (iii)

F l e x u r a l   s t r e n g t h = 48.98 - 1.33 * A + 6.15 * B + 5.32 * C - 6.41 * A 2 - 7.05 * B 2 - 4.79 * C 2 + 1.24 * A * B + 0.81 * A * C + 3.62 * B * C (iv)

F l e x u r a l   m o d u l u s = 1661.91 - 16.93 * A + 33.41 * B + 17.10 * C - 35.05 * A 2 - 14.22 * B 2 - 8.46 * C 2 - 19.37 * A * B + 6.22 * A * C + 8.01 * B * C (v)

I m p a c t   s t r e n g t h = 58.70 + 3.95 * A - 5.56 * B - 4.99 * C + 1.58 * A 2 + 1.85 * B 2 - 1.04 * C 2 + 0.55 * A * B - 3.26 * A * C + 2.54 * B * C (vi)

Table 6
Summary of the models.
Summary of the models.

Table 7
Initial ANOVA results and statistical parameters for tensile strength.
Initial ANOVA results and statistical parameters for tensile strength.

Table 8
ANOVA results and statistical parameters for tensile modulus.
ANOVA results and statistical parameters for tensile modulus.

Table 9
ANOVA results and statistical parameters for flexural strength.
ANOVA results and statistical parameters for flexural strength.

Table 10
Initial ANOVA results and statistical parameters for flexural modulus.
Initial ANOVA results and statistical parameters for flexural modulus.

Table 11
Initial ANOVA results and statistical parameters for impact strength.
Initial ANOVA results and statistical parameters for impact strength.

where A, B, and C are fiber length, fiber loading, and MAPE content, respectively. The significance of each model term was explained in terms of probability value (p-value). To eliminate non-effective terms and make regression model best fitted, ANOVA was repeated, and the results are presented in Tables 12-14. A fitted regression model with statistical significance is presented in the following equations:

T e n s i l e   s t r e n g t h = 32.19 - 0.15 * A + 2.86 * B + 0.96 * C - 2.88 * A 2 - 6.62 * B 2 - 1.63 * C 2 (vii)

F l e x u r a l   m o d u l u s = 1658.35 - 16.93 * A + 33.41 * B + 17.10 * C - 35.49 * A 2 - 14.66 * B 2 - 19.37 * A * B (viii)

I m p a c t   s t r e n g t h = 58.27 + 3.95 * A - 5.56 * B - 4.99 * C + 1.53 * A 2 + 1.79 * B 2 - 3.26 * A * C + 2.54 * B * C (ix)

Table 12
Final ANOVA results and statistical parameters for tensile strength.
Final ANOVA results and statistical parameters for tensile strength.

Table 13
Final ANOVA results and statistical parameters for flexural modulus.
Final ANOVA results and statistical parameters for flexural modulus.

Table 14
Final ANOVA results and statistical parameters for impact strength.
Final ANOVA results and statistical parameters for impact strength.

validation of regression models at 95% confidence interval was explained in terms of p, F, and R2 values. The insignificant “Lack-of-Fit” and p-values less than 0.05 confirmed the rejection of null-hypothesis. The model F (107.22) and p (<0.0001) values for tensile strength, F (560.03) and p (<0.0001) values for tensile modulus, F (475.55) and p (<0.0001) values for flexural strength, F (36.64) and p (<0.0001) values for flexural modulus, and F (53.86) and p (<0.0001) values for impact strength showed that the quadratic model was significant. The regression model was further analyzed by evaluating R2, adjusted R2, and predicted R2 which indicates the proportion of total variation in response variable predicted by model. The higher correlation coefficients confirm the suitability and correctness of the model. The adjusted R2 can be used to prevent probability error while predicted R2 indicates how well the model predicts responses for new observations. The adjusted R2 and predicted R2 are in reasonable agreement with the values of 0.9755 and 0.9433 for tensile strength, 0.9968 and 0.9816 for tensile modulus, 0.9963 and 0.9820 for flexural strength, 0.9304 and 0.8446 for flexural modulus, and 0.9586 and 0.9256 for impact strength, respectively. The R2 value for response variables implying that a high correlation exists between observed and predicted values. In addition, the adequacy of the model was confirmed using normal probabili-ty plots of the residuals, predicted versus actual values, and residuals versus predicted values plots. Figures 2-4 revealed that the error points are normally distributed along a straight line, implying that the model is adequate, and it represents the experimental data. The correlation between predicted response values and the actual values is shown in Figures 5-7, presenting uniformly distributed data points around the mean of the response variables. The linear regression fit is obtained with an R2 values of 0.9847, 0.9986, 0.9984, 0.9565, and 0.9767 for the tensile strength, tensile modulus, flexural strength, flexural modulus, and impact strength respectively, indicating that the model is accurately describing the experimental observations. In the end, the final plots of residuals versus predicted values (Figures 8-10) concluded that the residuals scattering for all mechanical responses were not significant which confirms that the proposed model was suitable. The maximum error (Tables 15-17) lies in between predicted and measured values were 4.10%, 0.26%, 1.51%, 1.01%, and 4.02% for the tensile strength, tensile modulus, flexural strength, flexural modulus, and impact strength, respectively.

Normal probability plot of the residuals for (a) tensile strength and (b) tensile modulus.
Figure 2
Normal probability plot of the residuals for (a) tensile strength and (b) tensile modulus.

Normal probability plot of the residuals for (a) flexural strength and (b) flexural modulus.
Figure 3
Normal probability plot of the residuals for (a) flexural strength and (b) flexural modulus.

Normal probability plot of the residuals for impact strength.
Figure 4
Normal probability plot of the residuals for impact strength.

Plots of predicted versus actual values for (a) tensile strength and (b) tensile modulus.
Figure 5
Plots of predicted versus actual values for (a) tensile strength and (b) tensile modulus.

Plots of predicted versus actual values for (a) flexural strength and (b) flexural modulus.
Figure 6
Plots of predicted versus actual values for (a) flexural strength and (b) flexural modulus.

Plots of predicted versus actual values for impact strength.
Figure 7
Plots of predicted versus actual values for impact strength.

Residuals versus predicted values plot for (a) tensile strength and (b) tensile modulus.
Figure 8
Residuals versus predicted values plot for (a) tensile strength and (b) tensile modulus.

Residuals versus predicted values plot for (a) flexural strength and (b) flexural modulus.
Figure 9
Residuals versus predicted values plot for (a) flexural strength and (b) flexural modulus.

Residuals versus predicted values plot for impact strength.
Figure 10
Residuals versus predicted values plot for impact strength.

Table 15
Comparison between measured and predicted values for tensile strength and modulus.
Comparison between measured and predicted values for tensile strength and modulus.

Table 16
Comparison between measured and predicted values for flexural strength and modulus.
Comparison between measured and predicted values for flexural strength and modulus.

Table 17
Comparison between measured and predicted values for impact strength.
Comparison between measured and predicted values for impact strength.

3.2. Effect of design parameters on tensile properties of biocomposites

Figure 11 (a)-(d) illustrate the 3D surface response plots for tensile strength and modulus as a function of fiber length and content, keeping the MAPE load fixed at 1 and 5 wt.%. The tensile strength was increased to a maximum value and then decreased slowly with the increase of fiber length. However, the tensile modulus was varied inversely with fiber length, and it was due to the severe fiber breakage and fibrillation of longer PALF during composite compounding. Irrespective of fiber length, the incorporation of PALF beyond a specified limit did not produce any significant reinforcing effect. This was because of the agglomeration and poor interfacial adhesion of fibers within the polymer matrix.

3D surface plots for the effect of fiber length and fiber loading on tensile strength (a) 1 wt.% MAPE, (b) 5 wt.% MAPE, and tensile modulus (c) 1 wt.% MAPE, (d) 5 wt.% MAPE.
Figure 11
3D surface plots for the effect of fiber length and fiber loading on tensile strength (a) 1 wt.% MAPE, (b) 5 wt.% MAPE, and tensile modulus (c) 1 wt.% MAPE, (d) 5 wt.% MAPE.

The variation in tensile properties with fiber length and compatibilizer content is shown in Figure 12 (a)-(d). The tensile strength and stiffness were increased to maximum values and then decreased with the increase of MAPE content. It was due to the excess removal of lignin compound in highly concentrated solution. Figure 13 (c)-(d) reveals the significant effect of fiber loading and MAPE content on tensile strength and modulus. The compatibilizer effect was pronounced only under high fiber loading condition (24 wt.%) which results in an increase of tensile modulus with the increase of MAPE content. It was due to the formation of rigid networks which constraint the movement of polymeric chains. The 3.6 wt.% MAPE compatibilized composite having 17% fiber content with 13 mm fiber length exhibits maximum tensile strength (Table 18).

3D surface plots for the effect of fiber length and MAPE on tensile strength (a) 8 wt.% PALF, (b) 24 wt.% PALF, and tensile modulus (c) 8 wt.% PALF, (d) 24 wt.% PALF.
Figure 12
3D surface plots for the effect of fiber length and MAPE on tensile strength (a) 8 wt.% PALF, (b) 24 wt.% PALF, and tensile modulus (c) 8 wt.% PALF, (d) 24 wt.% PALF.

3D surface plots for the effect of fiber loading and MAPE on tensile strength (a) 6 mm fiber length, (b) 20 mm fiber length, and tensile modulus (c) 6 mm fiber length, (d) 20 mm fiber length.
Figure 13
3D surface plots for the effect of fiber loading and MAPE on tensile strength (a) 6 mm fiber length, (b) 20 mm fiber length, and tensile modulus (c) 6 mm fiber length, (d) 20 mm fiber length.

Table 18
Optimum values of design parameters for maximum mechanical properties.
Optimum values of design parameters for maximum mechanical properties.

The main effect plot of each parameter on tensile properties of biocomposites is shown in Figure 14. Among all three critical variables, the fiber loading has highest degree of influence on tensile properties. The strength and stiffness of biocomposites were increased to maximum values and then decreased with the increase of fiber content.

Main effects plot for (a) tensile strength and (b) tensile modulus.
Figure 14
Main effects plot for (a) tensile strength and (b) tensile modulus.

3.3. Effect of design parameters on flexural properties of biocomposites

The interaction effect of fiber length and content on flexural properties of biocomposites is shown in Figure 15 (a)-(d). Like tensile strength, the flexural strength as well as flexural modulus were increased to maximum values and then gradually decreased with the increase of fiber length. The flexural modulus of biocomposites was increased with the increase of fiber content and it was due to the presence of hydroxyl groups in PALF which facilitate in binding of polymeric chains and the chains entanglement.

3D surface plots for the effect of fiber length and fiber loading on flexural strength (a) 1 wt.% MAPE, (b) 5 wt.% MAPE, and flexural modulus (c) 1 wt.% MAPE, (d) 5 wt.% MAPE.
Figure 15
3D surface plots for the effect of fiber length and fiber loading on flexural strength (a) 1 wt.% MAPE, (b) 5 wt.% MAPE, and flexural modulus (c) 1 wt.% MAPE, (d) 5 wt.% MAPE.

Figure 16 (a)-(d) reveals the combined effect of fiber length and MAPE content on flexural properties of biocomposites. It was observed that for all fiber length, the flexural strength and modulus were increased with the increase of compatibilizer load, and it was due to the formation of covalent linkages in between anhydride and hydroxyl (-OH) groups of fibers. The interaction effect of fiber loading and MAPE content on flexural properties of biocomposites is shown in Figure 17 (a)-(d). It is worth to noted that the reinforcement efficiency of PALF was significantly increased after the addition of compatibilizer. The interaction between fiber loading and MAPE content significantly affects the flexural strength (p-value BC<0.05). The 4.5 wt.% MAPE compatibilized composite having 22% fiber content with 14 mm length exhibits maximum flexural strength (Table 18).

3D surface plots for the effect of fiber length and MAPE on flexural strength (a) 8 wt.% PALF, (b) 24 wt.% PALF, and flexural modulus (c) 8 wt.% PALF, (d) 24 wt.% PALF.
Figure 16
3D surface plots for the effect of fiber length and MAPE on flexural strength (a) 8 wt.% PALF, (b) 24 wt.% PALF, and flexural modulus (c) 8 wt.% PALF, (d) 24 wt.% PALF.

3D surface plots for the effect of fiber loading and MAPE on flexural strength (a) 6 mm fiber length, (b) 20 mm fiber length, and flexural modulus (c) 6 mm fiber length, (d) 20 mm fiber length.
Figure 17
3D surface plots for the effect of fiber loading and MAPE on flexural strength (a) 6 mm fiber length, (b) 20 mm fiber length, and flexural modulus (c) 6 mm fiber length, (d) 20 mm fiber length.

Figure 18 illustrates the main effect plot of each individual parameter on flexural strength and modulus. It was observed that the fiber volume content has highest degree influence on the flexural behavior of developed composites. The addition of fibers beyond a limiting volume could not increase the strength of composite. However, the flexural modulus was increased progressively with the increase of fiber volume fraction. The fiber aspect ratio (length to diameter ratio) is a crucial parameter which dictates the reinforcing efficiency of fiber. The bending strength and stiffness were increased to maximum values and then decrease with the increase of fiber length. This kind of behavior was found in kenaf fiber-filled poly (butylene succinate) biocomposites (Thirmizir et al., 2011).

Main effects plot for (a) flexural strength and (b) flexural modulus.
Figure 18
Main effects plot for (a) flexural strength and (b) flexural modulus.

3.4. Effect of design parameters on impact strength of biocomposites

Figure 19 (a)-(b) illustrates the effect of fiber length and MAPE content on impact strength of biocomposites, keeping the fiber content was fixed at 8 and 24 wt.%. It was observed that the composite of higher aspect ratio (L/D) yields better impact toughness than the composite of low L/D ratio. Moreover, the low aspect ratio results in stress concentration and poor dispersion of fibers in HDPE matrix (Gamstedt et al., 2007). It is worth noted that the impact fracture strength of composites was decreased with the increase of compatibilizer content. This was attributed to the predominant action of fiber-fracture than the fiber pull-out.

3D surface plots for the effect of fiber length and MAPE on impact strength
Figure 19
3D surface plots for the effect of fiber length and MAPE on impact strength

In fiber composites having strong interfacial adhesion, fiber fracture is more common and dissipates lesser energy than the fiber pull-out (Wambua et al., 2003). The overall impact toughness of a fiber-reinforced composite (FRP) depends on the nature of the constituent elements, internal structure and geometry of the composite, fiber morphology, chemical composition, and interfacial adhesion between filler element and the matrix.

The interaction effect of fiber loading and MAPE content on impact strength of biocomposites is shown in Figure 20 (a)-(b). It was observed that the impact strength of composites was decreased with the increase of PALF content. This implies that the loading of PALF results in transition from ductile to brittle behavior. The important toughening mechanism in FRP is crack bridging by fibers associated with frictional sliding during fiber pull-out. It is worth noted that the composite shows poor mechanical strength at high fiber content (24 wt%)

3D surface plots for the effect of fiber loading and MAPE on impact strength (a) 6 mm fiber length, (b) 20 mm fiber length.
Figure 20
3D surface plots for the effect of fiber loading and MAPE on impact strength (a) 6 mm fiber length, (b) 20 mm fiber length.

The 1 wt.% MAPE compatibilized composite having 9% fiber content with 14 mm length exhibits maximum impact strength (Table 18).

The main effect plot of each individual parameter on impact strength of biocomposites is shown in Figure 21. The PALF content has highest degree of influence on impact strength. Moreover, the F value of fiber content (140.29) at 0.05 probability level confirms its highest significance on impact strength. The decrement in impact strength with increase of PALF content may be due to the agglomeration of fibers. Table 14 reported that the compatibilizer content is another crucial parameter with a significant F value (112.98). It was observed that the impact strength was decreased with the increase of MAPE content. The stronger interfacial adhesion results in the major occurrence of fiber fracture than that of fiber pull-out under loading condition.

Main effects plot for impact strength.
Figure 21
Main effects plot for impact strength.

3.5. Optimization of response variables

The numerical optimization method was employed to generate the optimal condition for each response variable. The goal was to maximize strength and modulus (Table 19); therefore, the target was set at highest value obtained from experimental results. The pre-defined goal, importance level, optimum and desirable values of each parameter are shown in Figures 22-23. To validate the Box-Behnken design model, a minimum of four samples were developed at the optimum condition of performance parameters. It was observed that the obtained experimental results (Table 20) are within a 5% difference from predicted values which proves the adequacy, reliability, and significance of current model.

Table 19
Optimization criteria and importance level
Optimization criteria and importance level

Optimum values of input and response variables.
Figure 22
Optimum values of input and response variables.

Desirability values of input and response variables.
Figure 23
Desirability values of input and response variables.

Table 20
Experimental results at optimum values of performance parameters.
Experimental results at optimum values of performance parameters.

4. Conclusions

The primary objective of this work was to optimize the mechanical properties of biocomposites using Box-Behnken RSM technique. The effects of fiber length, fiber content, and compatibilizer load on mechanical properties of PALF/HDPE biocomposites as summarized below.

Acknowledgments

The authors extend their gratitude to the Delhi Technological University, Delhi (India) for supporting this research work.

References

Ali, A., Yu, L., Liu, H., Khalid, S., Meng, L., & Chen, L. (2017). Preparation and characterization of starch‐based composite films reinforced by corn and wheat hulls.Journal of Applied Polymer Science,134(32), 45159. https://doi.org/10.1002/app.45159

Ali, A., Xie, F., Yu, L., Liu, H., Meng, L., Khalid, S., & Chen, L. (2018). Preparation and characterization of starch-based composite films reinfoced by polysaccharide-based crystals.Composites Part B: Engineering,133, 122-128. https://doi.org/10.1016/j.compositesb.2017.09.017

Alvarez, V. A., Kenny, J. M., & Vázquez, A. (2004). Creep behavior of biocomposites based on sisal fiber reinforced cellulose derivatives/starch blends.Polymer Composites,25(3), 280-288. https://doi.org/10.1002/pc.20022

Anand, P., & Anbumalar, V. (2014). Mechanical properties of cellulose-filled epoxy hybrid composites reinforced with alkali-treated hemp fiber.39(1), 46-55.

Arib, R. M. N., Sapuan, S. M., Hamdan, M. A. M. M., Paridah, M. T., & Zaman, H. M. D. K. (2004). A literature review of pineapple fibre reinforced polymer composites.Polymers and Polymer Composites,12(4), 341-348. https://doi.org/10.1177/096739110401200408

Ashenai Ghasemi, F., Daneshpayeh, S., Ghasemi, I., & Ayaz, M. (2016). An investigation on the Young’s modulus and impact strength of nanocomposites based on polypropylene/linear low-density polyethylene/titan dioxide (PP/LLDPE/TiO 2) using response surface methodology.Polymer Bulletin,73, 1741-1760. https://doi.org/10.1007/s00289-015-1574-2

Basu, A., Kundu, S., Sana, S., Halder, A., Abdullah, M. F., Datta, S., & Mukherjee, A. (2017). Edible nano-bio-composite film cargo device for food packaging applications.Food Packaging and Shelf Life,11, 98-105. https://doi.org/10.1016/j.fpsl.2017.01.011

Bledzki, A. K., Gassan, J., & Theis, S. (1998). Wood-filled thermoplastic composites.Mechanics of composite materials,34(6), 563-568. https://doi.org/10.1007/BF02254666

Bodirlau, R., Teaca, C. A., & Spiridon, I. (2013). Influence of natural fillers on the properties of starch-based biocomposite films.Composites Part B: Engineering ,44(1), 575-583. https://doi.org/10.1016/j.compositesb.2012.02.039

Cantero, G., Arbelaiz, A., Llano-Ponte, R., & Mondragon, I. (2003). Effects of fibre treatment on wettability and mechanical behaviour of flax/polypropylene composites.Composites science and technology,63(9), 1247-1254. https://doi.org/10.1016/S0266-3538(03)00094-0

Chen, H. C., Chen, T. Y., & Hsu, C. H. (2006). Effects of wood particle size and mixing ratios of HDPE on the properties of the composites. European Journal of Wood and Wood Products, 64(3), 172-177.

Dabade, B. M., Reddy, G. R., Rajesham, S., & Kiran, C. U. (2006). Effect of fiber length and fiber weight ratio on tensile properties of sun hemp and palmyra fiber reinforced polyester composites.Journal of reinforced plastics and composites,25(16), 1733-1738. https://doi.org/10.1177/07316844060684

Dash, B. N., Rana, A. K., Mishra, S. C., Mishra, H. K., Nayak, S. K., & Tripathy, S. S. (2000). Novel low-cost jute-polyester composite. II. SEM observation of the fractured surfaces.Polymer-Plastics Technology and Engineering,39(2), 333-350. https://doi.org/10.1081/PPT-100100033

Das, G., & Biswas, S. (2016). Effect of fiber parameters on physical, mechanical and water absorption behaviour of coir fiber-epoxy composites.Journal of reinforced plastics and composites ,35(8), 628-637. https://doi.org/10.1088/1757-899X/115/1/012012

Devi, L. U., Bhagawan, S. S., & Thomas, S. (1997). Mechanical properties of pineapple leaf fiber‐reinforced polyester composites.Journal of Applied Polymer Science ,64(9), 1739-1748. https://doi.org/10.1002/(SICI)10974628(19970531)64:9<1739::AID-APP10>3.0.CO;2-T

El-Shekeil, Y. A., Sapuan, S. M., Abdan, K., & Zainudin, E. S. (2012). Influence of fiber content on the mechanical and thermal properties of Kenaf fiber reinforced thermoplastic polyurethane composites.Materials & Design,40, 299-303. https://doi.org/10.1016/j.matdes.2012.04.003

El-Shekeil, Y. A., Sapuan, S. M., Jawaid, M., & Al-Shuja’a, O. M. (2014). Influence of fiber content on mechanical, morphological and thermal properties of kenaf fibers reinforced poly (vinyl chloride)/thermoplastic polyurethane poly-blend composites.Materials & Design,58, 130-135. https://doi.org/10.1016/j.matdes.2014.01.047

FAOSTAT, F. (2017). Food and agriculture organization of the United Nations (FAO).FAO Statistic Database. http://www.fao.org/statistics/databases/en/

Fung, K. L., Xing, X. S., Li, R. K. Y., Tjong, S. C., & Mai, Y. W. (2003). An investigation on the processing of sisal fibre reinforced polypropylene composites.Composites science and technology ,63(9), 1255-1258. https://doi.org/10.1016/S0266-3538(03)00095-2

Fung, K. L., Li, R. K. Y., & Tjong, S. C. (2002). Interface modification on the properties of sisal fiber‐reinforced polypropylene composites.Journal of Applied Polymer Science ,85(1), 169-176. https://doi.org/10.1002/app.10584

García‐Hernández, E., Licea‐Claveríe, A., Zizumbo, A., Alvarez‐Castillo, A., & Herrera‐Franco, P. J. (2004). Improvement of the interfacial compatibility between sugar cane bagasse fibers and polystyrene for composites.Polymer composites,25(2), 134-145. https://doi.org/10.1002/pc.20011

Gamstedt, E. K., Nygård, P., & Lindström, M. (2007). Transfer of knowledge from papermaking to manufacture of composite materials. In Proceedings of the 3rd wood fibre polymer composites international symposium, Bordeaux (p. 12).

Ghasemi, F. A., Ghasemi, I., Menbari, S., Ayaz, M., & Ashori, A. (2016). Optimization of mechanical properties of polypropylene/talc/graphene composites using response surface methodology.Polymer Testing,53, 283-292. https://doi.org/10.1016/j.polymertesting.2016.06.012

Heinicke, R. M., & Gortner, W. A. (1957). Stem bromelain-a new protease preparation from pineapple plants.Economic Botany,11(3), 225-234. https://doi.org/10.1007/BF02860437

Holt, G. A., Chow, P., Wanjura, J. D., Pelletier, M. G., & Wedegaertner, T. C. (2014). Evaluation of thermal treatments to improve physical and mechanical properties of bio-composites made from cotton byproducts and other agricultural fibers.Industrial Crops and Products,52, 627-632. https://doi.org/10.1016/j.indcrop.2013.11.003

Idicula, M., Neelakantan, N. R., Oommen, Z., Joseph, K., & Thomas, S. (2005). A study of the mechanical properties of randomly oriented short banana and sisal hybrid fiber reinforced polyester composites.Journal of Applied Polymer Science ,96(5), 1699-1709. https://doi.org/10.1002/app.21636

Ismail, H., Shuhelmy, S., & Edyham, M. R. (2002). The effects of a silane coupling agent on curing characteristics and mechanical properties of bamboo fibre filled natural rubber composites.European Polymer Journal,38(1), 39-47. https://doi.org/10.1016/S0014-3057(01)00113-6

Jandas, P. J., Mohanty, S., Nayak, S. K., & Srivastava, H. (2011). Effect of surface treatments of banana fiber on mechanical, thermal, and biodegradability properties of PLA/banana fiber biocomposites.Polymer Composites,32(11), 1689-1700. https://doi.org/10.1002/pc.21165

John, M. J., Francis, B., Varughese, K. T., & Thomas, S. (2008). Effect of chemical modification on properties of hybrid fiber biocomposites.Composites Part A: Applied Science and Manufacturing,39(2), 352-363. https://doi.org/10.1016/j.compositesa.2007.10.002

Joseph, K., Thomas, S., & Pavithran, C. (1996). Effect of chemical treatment on the tensile properties of short sisal fibre-reinforced polyethylene composites.Polymer,37(23), 5139-5149. https://doi.org/10.1016/0032-3861(96)00144-9

Kazayawoko, M., Balatinecz, J. J., & Matuana, L. M. (1999). Surface modification and adhesion mechanisms in woodfiber-polypropylene composites.Journal of materials science,34, 6189-6199. https://doi.org/10.1023/A:1004790409158

Kalapakdee, A., & Amornsakchai, T. (2014). Mechanical properties of preferentially aligned short pineapple leaf fiber reinforced thermoplastic elastomer: Effects of fiber content and matrix orientation.Polymer testing,37, 36-44. https://doi.org/10.1016/j.polymertesting.2014.04.008

Kengkhetkit, N., & Amornsakchai, T. (2012). Utilisation of pineapple leaf waste for plastic reinforcement: 1. A novel extraction method for short pineapple leaf fiber.Industrial Crops and Products,40, 55-61. https://doi.org/10.1016/j.indcrop.2012.02.037

Kengkhetkit, N., & Amornsakchai, T. (2014). A new approach to “Greening” plastic composites using pineapple leaf waste for performance and cost effectiveness.Materials & Design,55, 292-299. https://doi.org/10.1016/j.matdes.2013.10.005

Khalil, H. A., Issam, A. M., Shakri, M. A., Suriani, R., & Awang, A. Y. (2007). Conventional agro-composites from chemically modified fibres.Industrial Crops and Products,26(3), 315-323. https://doi.org/10.1016/j.indcrop.2007.03.010

Latif, S. S., Nahar, S., & Hasan, M. (2015). Fabrication and electrical characterization of bamboo fiber-reinforced polypropylene composite.Journal of reinforced plastics and composites , 34(3), 187-195.

Leblanc, J. L., Furtado, C. R., Leite, M. C., Visconte, L. L., & de Souza, A. M. (2007). Effect of the fiber content and plasticizer type on the rheological and mechanical properties of poly (vinyl chloride)/green coconut fiber composites.Journal of Applied Polymer Science ,106(6), 3653-3665. https://doi.org/10.1002/app.26567

Li, M. C., Wu, Q., Song, K., Cheng, H. N., Suzuki, S., & Lei, T. (2016). Chitin nanofibers as reinforcing and antimicrobial agents in carboxymethyl cellulose films: Influence of partial deacetylation.ACS Sustainable Chemistry & Engineering, 4(8), 4385-4395. https://doi.org/10.1021/acssuschemeng.6b00981

Liu, W., Drzal, L. T., Mohanty, A. K., & Misra, M. (2007). Influence of processing methods and fiber length on physical properties of kenaf fiber reinforced soy based biocomposites.Composites Part B: Engineering ,38(3), 352-359. https://doi.org/10.1016/j.compositesb.2006.05.003

Maheswari, C. U., Reddy, K. O., Muzenda, E., Shukla, M., & Rajulu, A. V. (2013). Mechanical properties and chemical resistance of short tamarind fiber/unsaturated polyester composites: Influence of fiber modification and fiber content.International Journal of Polymer Analysis and Characterization,18(7), 520-533. https://doi.org/10.1080/1023666X.2013.816073

Maffezzoli, A., Calo, E., Zurlo, S., Mele, G., Tarzia, A., & Stifani, C. (2004). Cardanol based matrix biocomposites reinforced with natural fibres.Composites science and technology ,64(6), 839-845. https://doi.org/10.1016/j.compscitech.2003.09.010

Mathew, L., Joseph, K. U., & Joseph, R. (2004). Isora fibres and their composites with natural rubber. Progress in Rubber Plastics and Recycling Technology, 20(4), 337-349.

Mhalla, M. M., Bahloul, A., & Bouraoui, C. (2017). Analytical models for predicting tensile strength and acoustic emission count of a glass fiber reinforced polyamide using response surface method.Journal of Alloys and Compounds,695, 2356-2364. https://doi.org/10.1016/j.jallcom.2016.11.108

Mishra, S., Naik, J. B., & Patil, Y. P. (2000). The compatibilising effect of maleic anhydride on swelling and mechanical properties of plant-fiber-reinforced novolac composites.Composites science and technology ,60(9), 1729-1735. https://doi.org/10.1016/S0266-3538(00)00056-7

Mishra, S., Misra, M., Tripathy, S. S., Nayak, S. K., & Mohanty, A. K. (2001). Potentiality of pineapple leaf fibre as reinforcement in PALF-polyester composite: Surface modification and mechanical performance.Journal of reinforced plastics and composites , 20(4), 321-334.

Mittal, M., & Chaudhary, R. (2018). Effect of fiber content on thermal behavior and viscoelastic properties of PALF/Epoxy and COIR/Epoxy composites. Materials Research Express, 5(12), 125305.

Mwaikambo, L. Y., & Ansell, M. P. (2002). Chemical modification of hemp, sisal, jute, and kapok fibers by alkalization.Journal of Applied Polymer Science ,84(12), 2222-2234. https://doi.org/10.1002/app.10460

Mohanty, S., Verma, S. K., Nayak, S. K., & Tripathy, S. S. (2004). Influence of fiber treatment on the performance of sisal-polypropylene composites.Journal of Applied Polymer Science ,94(3), 1336-1345. https://doi.org/10.1002/app.21161

Mohanty, S., Nayak, S. K., Verma, S. K., & Tripathy, S. S. (2004). Effect of MAPP as coupling agent on the performance of sisal-PP composites.Journal of reinforced plastics and composites ,23(18), 2047-2063. https://doi.org/10.1177/0731684404041711

Mohamed, O. A., Masood, S. H., & Bhowmik, J. L. (2015). Optimization of fused deposition modeling process parameters: a review of current research and future prospects.Advances in manufacturing, 3, 42-53. https://doi.org/10.1007/s40436-014-0097-7

Myers, R. H., Montgomery, D. C., & Anderson-Cook, C. M. (2016). Response surface methodology: process and product optimization using designed experiments. John Wiley & Sons.

Nakagaito, A. N., & Yano, H. (2008). The effect of fiber content on the mechanical and thermal expansion properties of biocomposites based on microfibrillated cellulose.Cellulose,15, 555-559. https://doi.org/10.1007/s10570-008-9212-x

Nakhaei, M. R., Mostafapour, A., & Naderi, G. (2017). Optimization of mechanical properties of PP/EPDM/clay nanocomposite fabricated by friction stir processing with response surface methodology and neural networks.Polymer Composites,38, E421-E432. https://doi.org/10.1002/pc.23942

Qaiss, A. E. K., Bouhfid, R., & Essabir, H. (2014). Natural fibers reinforced polymeric matrix: thermal, mechanical and interfacial properties.Biomass and bioenergy: processing and properties, 225-245. https://doi.org/10.1007/978-3-319-07641-6_14

Oushabi, A., Sair, S., Hassani, F. O., Abboud, Y., Tanane, O., & El Bouari, A. (2017). The effect of alkali treatment on mechanical, morphological and thermal properties of date palm fibers (DPFs): Study of the interface of DPF-Polyurethane composite.South African Journal of Chemical Engineering,23, 116-123. https://doi.org/10.1016/j.sajce.2017.04.005

Özturk, S. (2010). Effect of fiber loading on the mechanical properties of kenaf and fiberfrax fiber reinforced phenol-formaldehyde composites. Journal of Composite Materials. 44, 2265-2288. https://doi.org/10.1177/0021998310364265

Park, B. D., & Balatinecz, J. J. (1998). Short term flexural creep behavior of wood‐fiber/polypropylene composites.Polymer composites,19(4), 377-382. https://doi.org/10.1002/pc.10111

Pavithran, C., Mukherjee, P. S., Brahmakumar, M., & Damodaran, A. D. (1987). Impact properties of natural fibre composites.Journal of Materials Science Letters, 6, 882-884. https://doi.org/10.1007/BF01729857

Raj, R. G., Kokta, B. V., Dembele, F., & Sanschagrain, B. (1989). Compounding of cellulose fibers with polypropylene: Effect of fiber treatment on dispersion in the polymer matirx.Journal of Applied Polymer Science ,38(11), 1987-1996. https://doi.org/10.1002/app.1989.070381103

Rostamiyan, Y., Fereidoon, A., Mashhadzadeh, A. H., Ashtiyani, M. R., & Salmankhani, A. (2015). Using response surface methodology for modeling and optimizing tensile and impact strength properties of fiber orientated quaternary hybrid nano composite.Composites Part B: Engineering ,69, 304-316. https://doi.org/10.1016/j.compositesb.2014.09.031

Rostamiyan, Y., Fereidoon, A., Rezaeiashtiyani, M., Mashhadzadeh, A. H., & Salmankhani, A. (2015). Experimental and optimizing flexural strength of epoxy-based nanocomposite: effect of using nano silica and nano clay by using response surface design methodology.Materials & Design,69, 96-104. https://doi.org/10.1016/j.matdes.2014.11.062

Rostamiyan, Y., Fereidoon, A., Ghalebahman, A. G., Mashhadzadeh, A. H., & Salmankhani, A. (2015). Experimental study and optimization of damping properties of epoxy-based nanocomposite: Effect of using nanosilica and high-impact polystyrene by mixture design approach.Materials & Design(1980-2015),65, 1236-1244. https://doi.org/10.1016/j.matdes.2014.10.022

Rostamiyan, Y., Mashhadzadeh, A. H., & SalmanKhani, A. (2014). Optimization of mechanical properties of epoxy-based hybrid nanocomposite: Effect of using nano silica and high-impact polystyrene by mixture design approach.Materials & Design(1980-2015),56, 1068-1077. https://doi.org/10.1016/j.matdes.2013.11.060

Salasinska, K., Polka, M., Gloc, M., & Ryszkowska, J. (2016). Natural fiber composites: the effect of the kind and content of filler on the dimensional and fire stability of polyolefin-based composites.Polimery,61(4), 255-265 https://doi.org/10.14314/polimery.2016.255

Shi, A. M., Wang, L. J., Li, D., & Adhikari, B. (2013). Characterization of starch films containing starch nanoparticles: Part 1: Physical and mechanical properties.Carbohydrate Polymers,96(2), 593-601. https://doi.org/10.1016/j.carbpol.2012.12.042

Siddique, F. R., Khan, A. U., & Sheikh, R. A. (1984). Wood Text Abs. 1, 4196.

Subasinghe, A. D. L., Das, R., & Bhattacharyya, D. (2016). Parametric analysis of flammability performance of polypropylene/kenaf composites.Journal of materials science,51, 2101-2111. https://doi.org/10.1007/s10853-015-9520-0

Thirmizir, M. A., Ishak, Z. M., Taib, R. M., Rahim, S., & Jani, S. M. (2011). Kenaf‐bast‐fiber‐filled biodegradable poly (butylene succinate) composites: Effects of fiber loading, fiber length, and maleated poly (butylene succinate) on the flexural and impact properties.Journal of Applied Polymer Science ,122(5), 3055-3063. https://doi.org/10.1002/app.34046

Tripathy, S. S., Di Landro, L., Fontanelli, D., Marchetti, A., & Levita, G. (2000). Mechanical properties of jute fibers and interface strength with an epoxy resin.Journal of Applied Polymer Science ,75(13), 1585-1596. https://doi.org/10.1002/(SICI)1097-4628(20000328)75:13<1585::AID-APP4>3.0.CO;2-Q

Uhlig, H. (Ed.). (1998). Industrial enzymes and their applications. John Wiley & Sons.

Venkateshappa, S. C., Jayadevappa, S. Y., & Puttiah, P. K. W. (2012). Mechanical behavior of areca fiber reinforced epoxy composites.Advances in Polymer Technology,31(4), 319-330. https://doi.org/10.1002/adv.20255

Venkateshwaran, N., ElayaPerumal, A., Alavudeen, A., & Thiruchitrambalam, M. (2011). Mechanical and water absorption behaviour of banana/sisal reinforced hybrid composites.Materials & Design ,32(7), 4017-4021. https://doi.org/10.1016/j.matdes.2011.03.002

Wambua, P., Ivens, J., & Verpoest, I. (2003). Natural fibres: can they replace glass in fibre reinforced plastics?.Composites science and technology ,63(9), 1259-1264. https://doi.org/10.1016/S0266-3538(03)00096-4

Wong, S., Shanks, R., & Hodzic, A. (2004). Interfacial improvements in poly (3-hydroxybutyrate)-flax fibre composites with hydrogen bonding additives.Composites science and technology ,64(9), 1321-1330. https://doi.org/10.1016/j.compscitech.2003.10.012

Xie, X. L., Fung, K. L., Li, R. K. Y., Tjong, S. C., & Mai, Y. W. (2002). Structural and mechanical behavior of polypropylene/maleated styrene‐(ethylene‐co‐butylene)‐styrene/sisal fiber composites prepared by injection molding.Journal of Polymer Science Part B: Polymer Physics,40(12), 1214-1222. https://doi.org/10.1002/polb.10175

Notes

Peer Review under the responsibility of Universidad Nacional Autónoma de México.
Funding

The authors received no specific funding for this work.

Author notes

*Corresponding author. E-mail address: mohit.30mittal@gmail.com (M. Mittal).

Conflict of interest declaration

Conflict of interest

The authors have no conflict of interest to declare.

HTML generated from XML JATS by