Optimising the Workability and Strength of Concrete Modiﬁed with Anacardium Occidentale Nutshell Ash

: Strength failure persists both in structural and mechanical analysis. One of its prominent characteristics is the adequate provision for parameters that minimise or maximise strength objectives while satisfying boundary conditions. The previous optimisation of concrete strength usually neglects mix design mechanisms induced by optimisation. Recent efforts to accurately optimise the concrete compressive strength have factored in some of these mechanisms. However, optimising concrete strength modiﬁed with high silica and alumina precursors, and crucial mix design factors are rare. Consequently, this paper optimised the concrete workability and strength, incorporating binding, water/binder ratio, binder/aggregate ratio, and curing mechanisms using the Box–Behnken design approach (BBDA). A waste material, anacardium occidentale (cashew) nutshell ash, was valorised and used at 5, 10, and 15 wt.% of cement. The composites were made, cured and tested at 14–90 d. The results revealed a high precision between the experimental slump and the optimisation slump at 97% R 2 . In addition, a 5% increase in compressive strength was obtained compared with the target compressive strength. Besides, the correlation between the model equation obtained from this study and predictions of previous studies via BBDA yielded a strong relationship.


Introduction
In recent times, the reality of material science as a subject area is still of interest considering the behavioural mechanisms of materials and the suitable modelling of their mathematical relationships [1]. A concrete mix is optimised to reduce the constituents' costs and expended time on conducting the concrete trial mixes without compromising the quality and performance of concrete under an applied loading condition [2]. Aggregates and cement paste are the primary constituents of concrete (composite material). Although the water/cement ratio governs the quality of required cement paste while the required proportion of cement paste to attain a targeted concrete strength and quality depends on the aggregates' characteristics such as types, voids, shapes, and surface area [2]. Therefore, as reported in previous studies [3][4][5], optimising a concrete mix can be conducted by adjusting the critical levels of mix factors such as binder/aggregate ratio, water/binder ratio, cementitious (binding) material contents, and binding contents ratio.
Box-Behnken design (BBD) is a statistical and analytical method of response surface that examines the effect of a set of quantitative experimental variables of factors on the response [6]. Response surface methodology (RSM), a design of experiment (DOE), is usually engaged to identify a set of vital factors (operating conditions), which generate the experimental design. Design factors such as binders' ratio, water/binder ratio, curing ages, and binder/aggregate ratio were selected as independent (continuous) variables. At the same time, the compressive strength was considered as a response variable. The data were analysed statistically via Minitab 17 software, and the regression model was developed for concrete strength as a function of mixture variables. The developed approach to optimising the mixtures of concrete proportions was based on the experimental works (involving the required characteristics of concrete performance) and the statistical analysis of data generated; this would reduce time, cost, and the number of trial mixes.

Materials' Preparation and Characterisations
Cashew nutshell was valorised using the pyrolysis method, and about 23 wt.% of cashew nutshells was obtained as ash. The results of chemical compositions of the ash and Portland limestone cement (PLC) used, via XRF analyser, Model PW-1800, are presented in Table 1. The results indicated that cashew nutshell ash (CNA) met a 70% minimum required specification of SiO 2 + Al 2 O 3 + Fe 2 O 3 , recommended by BS EN 8615-2 [21], hence suitable for use as SCM. A 42.5 R PLC was used as a binder. Moreover, the fineness, specific gravity, and specific surface area were obtained for both CNA and PLC as 8.10 and 3.15%, 2.85 and 3.15 g/cm 3 , and 505 and 375 m 2 /kg, respectively, following the procedure outlined in BS EN 196-3 [22]. The particle size distribution (PSD) of the binders, as shown in Figure 1, was evaluated using Laser diffraction, Model Beckman Coulter LS-100. The results signified that more volume and water would be needed, and a higher capacity against an alkaline environment would be yielded when CNA replaces PLC [23].   Aggregates, fine aggregate (≤4.5 mm size) and coarse aggregate (granite 12.5 mm size) were used. The specific gravity (SG) and water absorption (WA) of aggregates were determined by weighing 2 kg of each sample, immersing in clean water at 25 • C for 24 h and drying in an oven for 24 h at 110 • C [11]. The WA results indicated 0.70 and 0.80% and 2.60 and 2.65 g/cm 3 for both fine aggregates (FA) and coarse aggregates (CA). In addition the moisture content (MC) was examined using a clean container covered with a lid. A 2 kg of each sample was oven-dried at 105 • C for 24 h, removed and cooled for 45 min, and then weighed. The results showed 0.3% and 0.2% MC for both FA and CA, respectively. As outlined in BS EN 12620 [24], Grading was conducted on the aggregates used, and the results are shown in Figure 2.

Design Proportions and Concrete Strength Specifications
BS EN 206 [25] was used to design the trial mix proportions, considering the materials' physical characteristics. Portland limestone cement (PLC) was replaced with CNA at 5, 10, and 15 wt. % proportions because a 15% maximum substitution of PLC with CNA exhibited an optimum compressive strength for structural applications [26,27]. Concrete grades M 25, M 30, and M 40 targeted the respective concrete strengths 25, 30, and 40 MPa at 28 days curing were adopted due to their general uses in the building construction. Ultimately, the trial mix design results are presented in Tables 2-4 for M 25,  M 30, and M 40, respectively.

Design Proportions and Concrete Strength Specifications
BS EN 206 [25] was used to design the trial mix proportions, considering the materials' physical characteristics. Portland limestone cement (PLC) was replaced with CNA at 5, 10, and 15 wt.% proportions because a 15% maximum substitution of PLC with CNA exhibited an optimum compressive strength for structural applications [26,27]. Concrete grades M 25, M 30, and M 40 targeted the respective concrete strengths 25, 30, and 40 MPa at 28 days curing were adopted due to their general uses in the building construction. Ultimately, the trial mix design results are presented in Tables 2-4

Sample Preparations and Experimental Tests
The mixing was carried out for about 15 min at 25 ± 3 • C temperature and 65 ± 5% relative humidity. The slump and density tests were performed on the concrete specimens outlined in BS EN 12350-2 [28] and BS EN 12390-7 [29], respectively. The slump cone mould was filled with fresh concrete in three equal layers and compacted with 35 tamping rod strokes. The mould was vertically lifted upward, and then the height difference between the collapsed concrete and the slump cone was measured. For the compressive strengths, mixes were prepared using the standard cube 150 mm 3 and crushed at 14-90 d following the BS EN 12390-3 [30]'s procedure. The cubes were removed from the moulds after 24 h and immersed in a curing tank for 14-90 d at 25 ± 3 • C temperatures and 65 ± 5% relative humidity. The force was applied to the concrete cube in a constant regime under a loading rate of 0.6 MPa per second until failure occurs. The values at failure were taken, and the mean of 3 crushed samples was used for the statistical analysis.

Box-Behnken Design (BBD)
The BBD was applied to evaluate the interaction between the selected continuous variables and the response variable. Minitab statistical software was used, and the BBD was created using the four (4) different continuous variables, as shown in Table 5. Thus, the experimental run number was estimated using the illustration presented in Equation (1) [6]. However, the optimisation of workability (slump) was statistically evaluated using the selected continuous variables, as shown in Table 5, except curing age (D), because slump is determined immediately after the mixing process, requiring no curing period.
where N is the number of experimental runs, v is the total continuous variables, and c p is the total applied central point. Following Equation (1), the number of the experimental run based on four (4) continuous variables and three (3) applied central points, as shown in Table 5, was twenty-seven (27). Besides, the full quadratic term, as illustrated in Equation (2), accurately optimises the relationship between the continuous variables and the response variable, hence yielding a high precision [6,8].
where Y is the response variable ( f c ), a 0 is model coefficient constant, x i and x j are continuous variables, a i is the linear coefficient, a ii : square coefficient, and a ij is the interaction coefficient.

Optimisation of Variables
The four continuous variables were minimised while maximising the response variable (compressive strength) via the optimisation concept. Based on this concept, the response variable was converted into a composite desirability function (CD) over the range of 0-1,  [6,7]. Consequently, the closer the composite desirability to 1, the better the optimisation [6,7].
From Equation (3), CD = 1 if the response value (compressive strength) is at its goal or target. Besides, CD = 0 if the response value is outside an acceptable region. The desirability is a composite based on the expressions illustrated in Equations (4)-(6) to maximise, minimise, and to make the response to be as close as possible to the target, respectively [6]. Therefore, Equations (4)-(6) were used to validate the optimisation study results.
where CD is the composite desirability function (0-1), r is the response, L, U, and T are lower, upper, and target values/limits, respectively, and w is the weight.

Statistical Analysis and Fitting of the Strength Model
In analysing the experimental data, analysis of variance (ANOVA) was used to evaluate the four continuous variables, as presented in Table 5 and considered for the development of compressive strength and the statistical model's fitting using the full quadratic terms. Degree of freedom (DF) is one of the statistical terms in ANOVA, which defines the number of final values (n) in the statistical analysis that is freely bound to change. The relationship is presented in Equation (7).
An adjusted sum of square (Adj SS) is the adjusted squared distance, as illustrated in Equation (8) for summation of all data points (n), between each observed data (Y) and the sample mean (Z) [2]. Moreover, the adjusted mean square (Adj MS), as illustrated in Equation (9), represents Adj SS ratio to the DF.
F-value signifies a ratio of Adj MS of the continuous factor to the Adj MS of the error. In the statistical analysis, an error (residual) connotes the statistical difference between the observed variables and the adopted model's predicted values. Thus, a factor with a significant effect yields a higher F-value [2]. P-value represents a significance level compared with the α-value of 0.05. Considering the null hypothesis, if P-value ≤ α-value, there is a significant effect; if P-value > α-value, then no significant impact [2].

Slump, Density, and Compressive Strengths
The slump results for M 25-40 are presented in Figure 3. The results indicated that the slump decreased as CNA content increased. About 13-33%, 14-36%, and 8-33% decrease in slump were observed for concrete grades M 25, M 30, and M 40, respectively, as the CNA content increased from 5-15%. The reason could be associated with CNA's physical characteristics, finer particle size, and lower SG compared with PLC. Consequently, water demand increased owing to more volume as CNA replaced PLC; this resulted in more cohesive concrete, hence reducing the concrete's workability [23]. The incorporation of finer pozzolans (cashew nutshell and rice husk ashes), as reported in previous studies, reduced the concrete slump in its fresh state [26,31]. Thus, a concrete modifying with CNA could be a workable mix of 15% substitution because an optimum slump of 150 mm for normal workable concrete specified by BS EN 12350-2 [28] was satisfied.

= (9)
F-value signifies a ratio of Adj MS of the continuous factor to the Adj MS of the error. In the statistical analysis, an error (residual) connotes the statistical difference between the observed variables and the adopted model's predicted values. Thus, a factor with a significant effect yields a higher F-value [2]. P-value represents a significance level compared with the α-value of 0.05. Considering the null hypothesis, if P-value ≤ α-value, there is a significant effect; if P-˃ value α-value, then no significant impact [2].

Slump, Density, and Compressive Strengths
The slump results for M 25-40 are presented in Figure 3. The results indicated that the slump decreased as CNA content increased. About 13-33%, 14-36%, and 8-33% decrease in slump were observed for concrete grades M 25, M 30, and M 40, respectively, as the CNA content increased from 5-15%. The reason could be associated with CNA's physical characteristics, finer particle size, and lower SG compared with PLC. Consequently, water demand increased owing to more volume as CNA replaced PLC; this resulted in more cohesive concrete, hence reducing the concrete's workability [23]. The incorporation of finer pozzolans (cashew nutshell and rice husk ashes), as reported in previous studies, reduced the concrete slump in its fresh state [26,31]. Thus, a concrete modifying with CNA could be a workable mix of 15% substitution because an optimum slump of 150 mm for normal workable concrete specified by BS EN 12350-2 [28] was satisfied.   Figure 4, the density decreased with increasing CNA content. Comparing with the control sample (S0), there was about 3-10% decrease in density as wt. % of CNA increased from 5-15% for all concrete grades. The decline in density with increasing CNA content could be attributed to the SG of CNA, which was 10% lesser than that of PLC; this reduced the in-   Figure 4, the density decreased with increasing CNA content. Comparing with the control sample (S0), there was about 3-10% decrease in density as wt.% of CNA increased from 5-15% for all concrete grades. The decline in density with increasing CNA content could be attributed to the SG of CNA, which was 10% lesser than that of PLC; this reduced the interfacial particle of PLC, hence decreasing the density. These results also supported the findings from previous studies [32,33], which reported that an increase in pozzolan content (CNA in this case) in the cement-pozzolan based concrete reduces the density of the tiny air bubbles of the pozzolan in the cementitious matrix. It was also evident from Figure 4 that the density increased as the concrete grades increased from M 25-M 40 due to denser, less porous, and low water content of the concrete mix [23]. Therefore, the partial substitution of PLC by CNA reduces the dead loads, contributing to the lightweight of a building.
As shown in Figure 5a-c, the concrete's compressive strengths (fc) increased with increasing CNA content. As the percentage substitution of CNA content increased from 5-15%, there was about 3-6% increase in strength for M 25-M 40 at 28 d curing than the control sample (100% PLC). The reason for strength increment could be attributed to the substantial contents of SiO 2 and Al 2 O 3 of CNA compared with PLC, which reacted with the cement's hydrating agents and resulted in the C-S-H and C-A-H; these are concrete's strengthening gel and bond, respectively [23]. Therefore, CNA could be used at 5-15% replacement level for structural applications because all concrete grades fulfilled the target strengths specified by BS EN 1992-1-1 [34] at 28 d curing. findings from previous studies [32,33], which reported that an increase in pozzolan content (CNA in this case) in the cement-pozzolan based concrete reduces the density of the tiny air bubbles of the pozzolan in the cementitious matrix. It was also evident from Figure 4 that the density increased as the concrete grades increased from M 25-M 40 due to denser, less porous, and low water content of the concrete mix [23]. Therefore, the partial substitution of PLC by CNA reduces the dead loads, contributing to the lightweight of a building. As shown in Figure 5a-c, the concrete's compressive strengths (fc) increased with increasing CNA content. As the percentage substitution of CNA content increased from 5-15%, there was about 3-6% increase in strength for M 25-M 40 at 28 d curing than the control sample (100% PLC). The reason for strength increment could be attributed to the substantial contents of SiO2 and Al2O3 of CNA compared with PLC, which reacted with the cement's hydrating agents and resulted in the C-S-H and C-A-H; these are concrete's strengthening gel and bond, respectively [23]. Therefore, CNA could be used at 5-15% replacement level for structural applications because all concrete grades fulfilled the target strengths specified by BS EN 1992-1-1 [34] at 28 d curing.

Design of Experiment (DOE)
Following the creation of response surface design in Minitab 17 using BBD with four continuous variables, as presented in Table 5, the statistical results are presented in Table  6. The fitted model analysis results to a response surface design are presented in Table 7. The linearly squared terms yielded the best precision and significant effect compared to the full quadratic terms, hence used for the modelling analysis. However, in modelling, both A 2 and B 2 were adjudged as nonsignificant terms and removed from the model. ANOVA examined the model's accuracy on continuous variables' influence on the compressive strength (fc).

Design of Experiment (DOE)
Following the creation of response surface design in Minitab 17 using BBD with four continuous variables, as presented in Table 5, the statistical results are presented in Table 6. The fitted model analysis results to a response surface design are presented in Table 7. The linearly squared terms yielded the best precision and significant effect compared to the full quadratic terms, hence used for the modelling analysis. However, in modelling, both A 2 and B 2 were adjudged as nonsignificant terms and removed from the model. ANOVA examined the model's accuracy on continuous variables' influence on the compressive strength (fc).
As revealed in Table 7, all terms, except A and C 2 , exhibited significant concrete strength effects because P-value was less than α-value (0.05). Consequently, the regression model between the response variable (fc) and the continuous variables (A, B, C, and D) is developed, as illustrated in Equation (10)  fits yielded standard distance (S) and coefficient of determination (R 2 ) as 0.279352 and 100%, respectively; this signified that the model is 100% fit to optimise all designed variables. f c = 1.61 − 16.11A + 121.7B − 9.37 C + 0.6379D − 10.0C 2 − 0.004938D 2 (10)  In practice, a balanced or nearly balanced design with a large number of observations does not significantly influence the residuals (the difference between the observed and the fitted response variable) if it departs moderately from a straight line or normality [6]. Thus, normally distributed residuals from the analysis are required for a balanced design. As shown in Figure 6, the residuals generally followed a straight-line pattern, hence no evidence of non-normality or unknown variables in the model. Furthermore, in a designed experiment, the order of observations influences the response variable if the residuals fluctuate in a random pattern around the centre line [6]. The result of versus order for the response variable, as shown in Figure 7, exhibited a randomly scattered pattern about zero, hence indicating that no evidence exists that the error terms are correlated with one another.  In optimising the workability (slump), linear term exhibited the best and significant effects compared with linearly squared, linearly interacted, and full quadratic terms. However, during the statistical analysis, the binder/aggregate ratio (B) could not be estimated and was removed. The binder/aggregate ratio (B) could not contribute to the response because it was statistically insignificant [8]. Moreover, concrete's slump is majorly influenced by binders' physical characteristics and the water/binder ratio [23]. Thus, the model summary is illustrated in Equation (11), yielding R 2 as 97%. This signifies that the model can optimise and predict the slump (S) of fresh concrete incorporating CNA for all designed factors at 95% confidence and prediction bounds for the concrete grades,   In optimising the workability (slump), linear term exhibited the best and significant effects compared with linearly squared, linearly interacted, and full quadratic terms. However, during the statistical analysis, the binder/aggregate ratio (B) could not be estimated and was removed. The binder/aggregate ratio (B) could not contribute to the response because it was statistically insignificant [8]. Moreover, concrete's slump is majorly influenced by binders' physical characteristics and the water/binder ratio [23]. Thus, the model summary is illustrated in Equation (11), yielding R 2 as 97%. This signifies that the model can optimise and predict the slump (S) of fresh concrete incorporating CNA for all designed factors at 95% confidence and prediction bounds for the concrete grades, In optimising the workability (slump), linear term exhibited the best and significant effects compared with linearly squared, linearly interacted, and full quadratic terms. However, during the statistical analysis, the binder/aggregate ratio (B) could not be estimated and was removed. The binder/aggregate ratio (B) could not contribute to the response because it was statistically insignificant [8]. Moreover, concrete's slump is majorly influenced by binders' physical characteristics and the water/binder ratio [23]. Thus, the model summary is illustrated in Equation (11), yielding R 2 as 97%. This signifies that the model can optimise and predict the slump (S) of fresh concrete incorporating CNA for all designed factors at 95% confidence and prediction bounds for the concrete grades, water/binder ratio

Main Effects of Continuous Variables on fc
The main effects of continuous variables on the response variable (fc) are presented in Figure 9 following the fitted model generation based on the stored model. It was evident from Figure 8 that the continuous variables influenced the response variable because the plotted lines were not horizontal [6]. An increase in curing age (D), binder/aggregate ratio (B), cementitious content (C), as shown in Figure 9, improved the concrete strength (fc). However, the main effects of curing age on the response (fc) signified that the strength could be optimised at 65 d against 90 d, hence minimising the curing time by about 28%. This result supported the findings of similar study [5] where the compressive strengths of geopolymer concrete incorporating corncob ash as SCM were optimised at 74 d and 83 d curing for concrete grades M 30 and M 40, respectively. Moreover, the concrete strength increased at a low water/binder ratio (A) compared with a high water/binder ratio due to the hydration process rate. The low water/binder content accelerates the setting times, compacts the matrix mechanism, and minimises the matrix pores, resulting in higher strength performance [32]. As reported in the previous studies [32,35], it was confirmed that higher cementitious and binder/aggregate contents with low water/binder ratio result in a higher compressive strength due to mechanical interlocking and cracking-arresting capacities of the cement-aggregate bond.

Main Effects of Continuous Variables on fc
The main effects of continuous variables on the response variable (fc) are presented in Figure 9 following the fitted model generation based on the stored model. It was evident from Figure 8 that the continuous variables influenced the response variable because the plotted lines were not horizontal [6]. An increase in curing age (D), binder/aggregate ratio (B), cementitious content (C), as shown in Figure 9, improved the concrete strength (fc). However, the main effects of curing age on the response (fc) signified that the strength could be optimised at 65 d against 90 d, hence minimising the curing time by about 28%. This result supported the findings of similar study [5] where the compressive strengths of geopolymer concrete incorporating corncob ash as SCM were optimised at 74 d and 83 d curing for concrete grades M 30 and M 40, respectively. Moreover, the concrete strength increased at a low water/binder ratio (A) compared with a high water/binder ratio due to the hydration process rate. The low water/binder content accelerates the setting times, compacts the matrix mechanism, and minimises the matrix pores, resulting in higher strength performance [32]. As reported in the previous studies [32,35], it was confirmed that higher cementitious and binder/aggregate contents with low water/binder ratio result in a higher compressive strength due to mechanical interlocking and cracking-arresting capacities of the cement-aggregate bond.

Main Effects of Continuous Variables on fc
The main effects of continuous variables on the response variable (fc) are presented in Figure 9 following the fitted model generation based on the stored model. It was evident from Figure 8 that the continuous variables influenced the response variable because the plotted lines were not horizontal [6]. An increase in curing age (D), binder/aggregate ratio (B), cementitious content (C), as shown in Figure 9, improved the concrete strength (fc). However, the main effects of curing age on the response (fc) signified that the strength could be optimised at 65 d against 90 d, hence minimising the curing time by about 28%. This result supported the findings of similar study [5] where the compressive strengths of geopolymer concrete incorporating corncob ash as SCM were optimised at 74 d and 83 d curing for concrete grades M 30 and M 40, respectively. Moreover, the concrete strength increased at a low water/binder ratio (A) compared with a high water/binder ratio due to the hydration process rate. The low water/binder content accelerates the setting times, compacts the matrix mechanism, and minimises the matrix pores, resulting in higher strength performance [32]. As reported in the previous studies [32,35], it was confirmed that higher cementitious and binder/aggregate contents with low water/binder ratio result in a higher compressive strength due to mechanical interlocking and cracking-arresting capacities of the cement-aggregate bond.   Figure 9. The main effect of continuous variables on response value. A is water/binder ratio, B is binder/aggregate ratio, C is CNA/(CNA+PLC), and D is curing age (d).

Operating Conditions of Continuous Variables on fc
Operating conditions illustrate the interaction between the continuous variables and the response variable and visualise the effects in both 2-D contour lines and 3-D surface, holding other continuous variables constant. The contour indicates the highest impact on the operating condition with the darkest green colour. Notwithstanding, any alteration to the continuous variables would change the shape patterns of the working conditions. Figure 10a,b present the visualised effects on fc in 2-D contour lines and 3-D surface, respectively. As shown in Figure 10, it was evident that the strength increased with increasing B but at a decreasing A, holding both C and D at 10% and 52 d constant, respectively. Therefore, a response (fc) greater than 52 MPa, could be attained at 0.31 and 0.42 for B and A, respectively.

Operating Conditions of Continuous Variables on fc
Operating conditions illustrate the interaction between the continuous variables and the response variable and visualise the effects in both 2-D contour lines and 3-D surface, holding other continuous variables constant. The contour indicates the highest impact on the operating condition with the darkest green colour. Notwithstanding, any alteration to the continuous variables would change the shape patterns of the working conditions. Figure 10a,b present the visualised effects on fc in 2-D contour lines and 3-D surface, respectively. As shown in Figure 10, it was evident that the strength increased with increasing B but at a decreasing A, holding both C and D at 10% and 52 d constant, respectively. Therefore, a response (fc) greater than 52 MPa, could be attained at 0.31 and 0.42 for B and A, respectively.

Figure 10. Effects of B and A on fc (MPa) (a) 2-D contour lines (b) 3-D surface.
A is water/binder ratio, B is binder/aggregate ratio, C is CNA/(CNA + PLC), and D is curing age (d).

Influence of D and A on fc
The visualised effects of D and A on fc are presented in Figure 11a,b for 2-D contour lines and 3-D surface, respectively. It was revealed that the response (fc) could be optimised (˃45 MPa) at about 65 d of curing age (D) and a low water/binder ratio ≤ 0.42, holding both B and C at 0.25 and 10% constant, respectively.

Influence of D and A on fc
The visualised effects of D and A on fc are presented in Figure 11a,b for 2-D contour lines and 3-D surface, respectively. It was revealed that the response (fc) could be optimised (¿45 MPa) at about 65 d of curing age (D) and a low water/binder ratio ≤ 0.42, holding both B and C at 0.25 and 10% constant, respectively.

Operating Conditions of Continuous Variables on fc
Operating conditions illustrate the interaction between the continuous variables and the response variable and visualise the effects in both 2-D contour lines and 3-D surface, holding other continuous variables constant. The contour indicates the highest impact on the operating condition with the darkest green colour. Notwithstanding, any alteration to the continuous variables would change the shape patterns of the working conditions. Figure 10a,b present the visualised effects on fc in 2-D contour lines and 3-D surface, respectively. As shown in Figure 10, it was evident that the strength increased with increasing B but at a decreasing A, holding both C and D at 10% and 52 d constant, respectively. Therefore, a response (fc) greater than 52 MPa, could be attained at 0.31 and 0.42 for B and A, respectively.

Influence of D and A on fc
The visualised effects of D and A on fc are presented in Figure 11a,b for 2-D contour lines and 3-D surface, respectively. It was revealed that the response (fc) could be optimised (˃45 MPa) at about 65 d of curing age (D) and a low water/binder ratio ≤ 0.42, holding both B and C at 0.25 and 10% constant, respectively.

Influence of D and B on fc
As shown in Figure 12a,b, the response (fc) increased with increasing D and B, holding A and C at 0.52 and 10%, respectively. However, the interaction revealed that the response could be optimised (¿50 MPa) at about 65 d of curing age and 0.31 B. As shown in Figure 12a,b, the response (fc) increased with increasing D and B, holding A and C at 0.52 and 10%, respectively. However, the interaction revealed that the response could be optimised (˃50 MPa) at about 65 d of curing age and 0.31 B.

Optimisation
Based on the model stored, Figure 13 presents the optimisation results for the response (fc) following the desirability functions' combination into single composite desirability (CD). The vertical brown and horizontal blue lines, as shown in Figure 13, signified a current setting and a current response (fc), respectively. The current set showed A, B, C, and D as 0.42, 0.31, 0.15, and 70 d, respectively, yielding an optimised response (fc) of 54 MPa at 1.000 composite desirability (CD). Therefore, the strength optimisation results supported a relationship, as illustrated in Equation (10), such that CD = 1 because the current response, fc (54 MPa), is higher than the target/upper strength (51.48 MPa) obtained from the experimental works. Figure 13. Optimisation plot for the response (fc, MPa). A is water/binder ratio, B is binder/aggregate ratio, C is CNA/(CNA+PLC), and D is curing age (d).

Optimisation
Based on the model stored, Figure 13 presents the optimisation results for the response (fc) following the desirability functions' combination into single composite desirability (CD). The vertical brown and horizontal blue lines, as shown in Figure 13, signified a current setting and a current response (fc), respectively. The current set showed A, B, C, and D as 0.42, 0.31, 0.15, and 70 d, respectively, yielding an optimised response (fc) of 54 MPa at 1.000 composite desirability (CD). Therefore, the strength optimisation results supported a relationship, as illustrated in Equation (10), such that CD = 1 because the current response, fc (54 MPa), is higher than the target/upper strength (51.48 MPa) obtained from the experimental works.

Influence of D and B on fc
As shown in Figure 12a,b, the response (fc) increased with increasing D and B, holding A and C at 0.52 and 10%, respectively. However, the interaction revealed that the response could be optimised (˃50 MPa) at about 65 d of curing age and 0.31 B. ure 12. Effects of D and B on fc (MPa) (a) 2-D contour lines (b) 3-D surface. A is water/binder ratio, B is bindaggregate ratio, C is CNA/(CNA + PLC), and D is curing age (d).

Optimisation
Based on the model stored, Figure 13 presents the optimisation results for the response (fc) following the desirability functions' combination into single composite desirability (CD). The vertical brown and horizontal blue lines, as shown in Figure 13, signified a current setting and a current response (fc), respectively. The current set showed A, B, C, and D as 0.42, 0.31, 0.15, and 70 d, respectively, yielding an optimised response (fc) of 54 MPa at 1.000 composite desirability (CD). Therefore, the strength optimisation results supported a relationship, as illustrated in Equation (10), such that CD = 1 because the current response, fc (54 MPa), is higher than the target/upper strength (51.48 MPa) obtained from the experimental works. Figure 13. Optimisation plot for the response (fc, MPa). A is water/binder ratio, B is binder/aggregate ratio, C is CNA/(CNA+PLC), and D is curing age (d). Figure 13. Optimisation plot for the response (fc, MPa). A is water/binder ratio, B is binder/aggregate ratio, C is CNA/(CNA+PLC), and D is curing age (d). Figure 14 compares the developed model with the earlier studies to establish and prove the model's accuracy, precision, and actual forecasting capability toward the concrete modifying with SCMs. As indicated in Figure 14, Kumar and Prasad [36] replaced fly ash and silica fume at 6-15 wt.% of cement using a design water/binder and binder/aggregate ratios of 0.4-0.5 and 0.21-0.22, respectively. Moreover, Duana et al. [37] substituted both metakaolin and blast furnace slag at 10-20 wt.% of cement using water/binder and binder/aggregate ratios of 0.50 and 0.20, respectively. Compared the experimental compressive strengths (Efc) of these previous studies [36,37] with the developed model, as illustrated in Equation (10), the predictive results (Pfc) yielded a high precision with 62-88% R 2 , hence indicating a good relationship. Therefore, this developed model established the reliability to predict and optimise the concrete compressive strength, mainly modified and produced with SCMs.

Comparison of Developed Model with the Previous Studies
prove the model's accuracy, precision, and actual forecasting capability toward the concrete modifying with SCMs. As indicated in Figure 14, Kumar and Prasad [36] replaced fly ash and silica fume at 6-15 wt. % of cement using a design water/binder and binder/aggregate ratios of 0.4-0.5 and 0.21-0.22, respectively. Moreover, Duana et al. [37] substituted both metakaolin and blast furnace slag at 10-20 wt. % of cement using water/binder and binder/aggregate ratios of 0.50 and 0.20, respectively. Compared the experimental compressive strengths (Efc) of these previous studies [36,37] with the developed model, as illustrated in Equation (10), the predictive results (Pfc) yielded a high precision with 62-88% R 2 , hence indicating a good relationship. Therefore, this developed model established the reliability to predict and optimise the concrete compressive strength, mainly modified and produced with SCMs.

Conclusions
A simplified mechanism, optimising the concrete workability and strength modified with high silica and alumina precursor, was developed based on the statistical analysis obtained through an experimentally designed work. As a result of this mechanism, the following outcomes are drawn: At 5-15 wt. % substitution of CNA, there was about 3-6% increase in concrete's compressive strengths for all concrete grades (M 25, M 30, and M 40) at 28 d curing. Moreover, the optimum values of water/binder ratio (A), binder/aggregate ratio (B), and SCM/cement ratio (C) exhibited a 5% increase in compressive strength at 28% decrease in curing age (D) compared with the target variables. This signifies a significant cost and time saving while maintaining the composite product's quality and strength performance. On the other hand, as illustrated in Equation (10), the developed model has proved to be efficient such that a good relationship with high precision was obtained when compared with the previous concrete strength incorporating SCMs. In addition the slump obtained via optimisation and that of experimental values yielded the same results, hence exhibiting a strong precision, as illustrated in Equation (11); this can be used in the prediction of normal workability of concrete modified with CNA while saving cost and time of preparing experimental works. Finally, the possibility of using less water for the optimisation of concrete strength is attainable by blending PLC with anacardium occidentale nutshell ash (cashew nutshell ash, CNA). In addition this research can be applied in the normal reinforced concrete work without vibration.

Conclusions
A simplified mechanism, optimising the concrete workability and strength modified with high silica and alumina precursor, was developed based on the statistical analysis obtained through an experimentally designed work. As a result of this mechanism, the following outcomes are drawn: At 5-15 wt.% substitution of CNA, there was about 3-6% increase in concrete's compressive strengths for all concrete grades (M 25, M 30, and M 40) at 28 d curing. Moreover, the optimum values of water/binder ratio (A), binder/aggregate ratio (B), and SCM/cement ratio (C) exhibited a 5% increase in compressive strength at 28% decrease in curing age (D) compared with the target variables. This signifies a significant cost and time saving while maintaining the composite product's quality and strength performance. On the other hand, as illustrated in Equation (10), the developed model has proved to be efficient such that a good relationship with high precision was obtained when compared with the previous concrete strength incorporating SCMs. In addition the slump obtained via optimisation and that of experimental values yielded the same results, hence exhibiting a strong precision, as illustrated in Equation (11); this can be used in the prediction of normal workability of concrete modified with CNA while saving cost and time of preparing experimental works. Finally, the possibility of using less water for the optimisation of concrete strength is attainable by blending PLC with anacardium occidentale nutshell ash (cashew nutshell ash, CNA). In addition this research can be applied in the normal reinforced concrete work without vibration.  Data Availability Statement: The data presented in this study are available on request from the corresponding author.