Effects of Sand/Fly Ash and the Water/Solid Ratio on the Mechanical Properties of Engineered Geopolymer Composite and Mix Design Optimization

: This paper presents the results of an experimental study that investigated the effects of two parameters: sand/fly ash (S/FA) ratio and water to geopolymer-solid (W/GS) ratio on the engineered geopolymer composite. The trial mix designs were optimized using the response surface method. These parameters influence the properties of the fresh and hardened geopolymer matrix, such as slump flow, compressive strength, flexural strength, elastic modulus, flexural toughness, ductility index and drying shrinkage. The optimizing process was conducted by developing statistical models using the response surface methodology (RSM) technique. The developed models were statistically validated and could be used to determine the desired response of engineered geopolymer composite (EGC) with a significance level of more than 95%. In this study, the optimized values of the S/FA ratio and W/GS ratio were obtained as 0.341701 and 0.225184, respectively. To validate the optimized S/FA ratio and W/GS ratio, an experimental study was performed, and a difference of less than 5% was found between predicted and experimental results.


Introduction
The growing demand for ordinary portland cement (OPC) due to the extensive use of concrete by the global construction industry has created many concerns regarding sustainability and the environment. Currently, concrete is the world's most extensive consumable material after water. In contrast, cement production contributes approximately 10% to the global CO2 emissions and consumes about 15% of the total industrial energy produced, which is causing a significant burden on environmental and sustainability measures [1]. In the last few years, researchers have been making efforts to develop sustainable and environmentally-friendly solutions for producing an extensive amount of concrete for the global construction industry. In such efforts, geopolymer binders are being investigated as an innovative as well as a green alternative to traditional concrete [2]. Geopolymer binders are formed when a precursor (a powder which contains a high content of alumina and silica) interacts with an alkaline solution in the presence of heat; then, the system undergoes a chain reaction [3]. Many types of industrial byproducts, which have faced criticism regarding their safe disposal, are being used as a precursor to geopolymer binders [4]. Among the various types of industrial byproducts, fly ash is the most commonly used precursor, which is a readily available byproduct produced in coal-fired power plants. Fly ash is usually an amorphous powder that is rich in alumina and silica content, which can meet the requirements of sustainability and a green environment [5]. Many existing research studies have reported that the fly-ash-based geopolymer concrete has shown superior mechanical properties, high durability, and refined microstructures as compared to traditional cement concrete [6].
The superior compressive strength achievements of geopolymer binders create an opportunity for their application in the construction industry [7]. However, the brittleness, low tensile strength, small ultimate strain, and susceptibility to cracking of geopolymer binders are considered significant weaknesses which hinder their application in complex construction products [8]. In an effort to overcome the reported weaknesses, short and randomly oriented fibers were mixed in the matrix as an intrinsic reinforcement [9]. The short fibers in the cementing matrix can bridge the cracks, which restricts the crack propagation by enhancing the tensile strength, ductility index, and post-crack performance [10]. Many studies focused on the development of fiber-reinforced concrete (FRC) have used a different type of microfibers, such as steel [11], glass [12], basalt [13] and polyvinyl alcohol (PVA) [14]. All such types of fibers have shown various advantages and disadvantages as an intrinsic reinforcement cement matrix. For example, steel fibers (SF) improve the tensile strength and freezethaw resistance of concrete [15]. However, concerns around durability [16], cost and effects on workability [17] have been identified as the disadvantages of steel fibers. Research findings have shown that glass fibers were not stable in the alkaline environment [18]. Basalt fibers as an intrinsic reinforcement caused a softening behavior after the first crack appeared [19]. Based on various research studies, it has been found that PVA fibers are the most favorable type of fibers for achieving a high performance of FRC because of their stability in the alkaline environment, possessing high tensile strength, and because they have a significant value of elastic modulus compared to other types of fibers [20]. In some studies, researchers observed an enhancement in the mechanical properties and post-crack performance of concrete using PVA fibers [21,22].
It has been observed that in a conventional fiber reinforced concrete (FRC), the load-carrying capacity decreases after the initiation of the first crack [23]. However, advancements made in the field of FRC have resulted in the development of an engineered cementitious composite (ECC) [23]. This is classified as a high-performance FRC which is capable of achieving multiple cracking patterns with an increase in the applied load after the first crack appears [24]. Multiple cracking patterns indicate high ductility, energy absorption capacity, crack growth resistance, and toughness of ECC under flexural and tensile stresses [25]. A range of achievements in the mechanical and post-crack properties of ECC, such as compressive strength , flexural strength (10-30 MPa), and ultimate strain (1%-8%), make it a promising material for retrofitting and structural repair jobs [26].
So far, many efforts have been dedicated to improving the mechanical properties of conventional fiber-reinforced geopolymer concrete (FRGC) exhibiting softening behavior. The authors of [27] improved the flexural strength of geopolymer by adding micro steel fibers, while the authors of [28] enhanced compressive strength using basalt fiber in geopolymer. The authors of [29] observed improvements in split tensile strength and flexural strengths by adding glass fibers in geopolymer concrete. However, there is not much research work reported in the literature regarding the development of high-performance fiber-reinforced geopolymer concrete or known as engineered geopolymer composite (EGC). The mechanical and post-crack properties of EGC are considered similar, as discussed in the literature for engineered cementitious composites (ECC) [30].
Similar to that of ECC, the performance of EGC depends on multiple factors, such as slump flow of fresh mixture, compressive and flexural strength requirements, elastic modulus, flexural toughness, ductility index and drying shrinkage [30], therefore, it is a challenging task for a concrete expert to develop an appropriate mix design for achieving an optimized performance of EGC. The optimization is referred to as a systematic way of improving the reliability of the products. The Taguchi method [31], Plackett-Burman [32] and RSM [33] are reported as the most commonly-used optimization techniques. However, RSM is considered an advanced technique for systematically designing experiments, obtaining statistically significant models, and achieving optimum conditions for several desirable responses simultaneously [34,35]. The RSM technique is a widely accepted modeling and optimizing tool in the concrete industry for achieving the highest responses by acquiring an optimized mixture composition. Ghafari et al. (2014) [36] developed RSM-based statistical models for predicting the influence of hybrid fibers on the performance of self-compacting and ultra-high performance concrete. Similarly, another research reported the optimized effects of fiber volume fraction in steel fiber-reinforced concrete for achieving maximum ductility [37]. In contrast, Al-Alaily (2016) [38] applied the RSM technique to optimize the chloride permeability and mechanical properties of metakaolin-based concrete. The available literature shows that the RSM technique is successfully applied for the optimization of OPC-based concrete mixes. However, limited literature is available on its application in the field of geopolymer technology. An experimental study discussed the optimized mixture components of a lightweight geopolymer composite by minimizing the curing temperature, binder content, and curing time for achieving maximum compressive strength [39]. The EGC properties are primarily influenced by several parameters, such as NaOH molarity, Na2SiO3/NaOH ratio, curing temperature, sand to fly ash ratio, and water to geopolymer solids ratio; therefore, it is essential to determine the optimum level of dosages of the mix ingredients for achieving the most desired performance of the matrix [30].
In a previous experimental program,   [20] applied RSM to optimize the effect of the primary polymeric variables (NaOH molarity, Na2SiO3/NaOH and curing temperature) to maximize the mechanical and post-crack performance of EGC. For a further improvement in achieving the optimum mix design of EGC, the principal aim of this study was to apply the RSM technique to optimize the effects of the sand/fly ash ratio and water/geopolymer-solid ratio on the EGC properties. The effects were investigated on the slump flow, compressive strength, flexural strength, flexural toughness, ductility index, elastic modulus, and drying shrinkage. In the final phase, an experimental program was designed for validating the optimum parameters as determined from the RSM application.

Research Significance
Conventional concrete made of Portland cement is facing much criticism due to it using a massive amount of natural resources and contributing to a substantial amount in global greenhouse gas (GHG) emissions. Therefore, for the last few years, concrete researchers have been working on alternative solutions that should be environmentally friendly and sustainable. In such efforts, geopolymer binders are becoming a favorite research topic. However, there are many concerns and issues to be resolved before they take a significant share in concrete application. This paper is significant in that it is focused on the optimization of engineered geopolymer mixture design for achieving the desired structural responses.

Material Properties
For this experimental program, fly ash (FA) was obtained from the Manjung coal-fired power station located in Perak, Malaysia [40]. The supplied fly ash was tested to verify the specified characteristics; the morphological study was performed using a field emission scanning electron microscope (FESEM). Figure 1a shows that the provided fly ash complied with the specified morphology. Table 1 shows the elemental composition and surface area of FA samples as obtained using the X-ray fluorescence (XRF) and Brunauer Emmet Taller (BET) techniques. Based on the characterization test results, the supplied fly ash is classified as high-calcium as per ASTM 618-10 (CaO > 10%; SiO2 + Al2O3 + Fe2O3 > 70%). Table 2 shows the gradation of microsilica sand that was supplied by a local supplier in the state of Perak, Malaysia. Previous research suggested that the sodium-based alkaline solution performs better in achieving compressive strength, flexural strength, and deflective capacity of EGC [41]; 99% pure sodium hydroxide (NaOH) in combination with water glass, also called sodium silicate (Na2SiO3), was used for the synthesis of EGC. Its chemical characteristics are given in Table 3. The Kurary group, Japan, provided polyvinyl alcohol (PVA) fibers; their mechanical and physical properties are given in Table 4, and the FESEM image is shown in Figure 1b.

Response Surface Method (RSM) Application for Experiment Design
The trial and error method has conventionally been practiced to find out the influence on the responses by changing the value of one variable at a time. This process is time-consuming, especially for multivariable systems. Therefore, it is needed to utilize a robust technique for solving multiobjective optimization problems. The response surface method (RSM) is a robust approach for optimizing complex problems through the development of statistical interaction between factors (independent variables) and responses (dependent variables). The process of optimization through RSM is initiated with the systematic designing of experiments and through the collections of responses by performing experimental programs. Finally, statistical models are developed for the prediction of each response, and optimization is performed for the desired targets [36]. The most practically adopted method for experimental design, statistical modeling, and optimization is the central composite design (CCD) of RSM [42]. The process of selection of data points through CCD using RSM is illustrated in Figure 2. The experiments are designed so that the space of independent variables must include corner points (+1 and −1), a center point, and extreme points (+ and -) [43]. Quadratic models that are used to estimate each response are developed using the 2nd-degree polynomial Equation (1) [36].
where the predicted response (slump flow, compressive strength, flexural strength, elastic modulus, flexural toughness, and ductility index of EGC) is represented by " ". Moreover, , , , and are the constant coefficient, linear coefficient, the quadratic coefficient, and interaction coefficient, respectively. Furthermore, the independent variables in Equation (1) are denoted by and . Table 5 shows the defined independent variables (factors) and their corresponding code and levels.
An experimental program was designed using commercially available Design-Expert ® software version 10 ( Table 6). The ranges of sand to fly ash and water to geopolymer solids were selected as 0 to 0.8 and 0.22 to 0.26, respectively, while other factors such as NaOH molarity, Na2SiO3/NaOH, and curing temperature were kept constant at 8 M and 2.5 and 60 °C, respectively [30]. Table 5. Boundaries of factors in response surface methodology (RSM).

Factor
Code Note: Water to geopolymer (W/GP) solids and sand/fly ash were selected according to the central composite design configurations for two factors using Design expert ® software. Geopolymer solids (solids in NaOH and Na2SiO3 solution plus fly ash), PVA fibers, NaOH molarity, and Na2SiO3/NaOH and curing temperature were selected from the available literature on EGC [30]. The dosage of the alkaline solution was selected from the literature [20]. Extra water was added to maintain desired W/GP solids. a Quantity in a mass ratio of fly ash. b Quantity in the volume fraction of the total volume of material.

Mixing and Casting of Samples
To optimize the water to geopolymer solid (W/GS) ratio and the sand to fly ash ratio (S/FA), nine EGC mixes, as shown in Table 6, were prepared. A total of 320 g of NaOH pellets were dissolved in the tap water for making one litter of 8 M NaOH solution. Since it is an exothermic reaction, the NaOH solution was placed in the laboratory until it cooled down. To achieve a homogeneous paste, mixing activity was done for a total of five minutes using the Hobart mixture. In the beginning, fly ash and sand samples were dry-mixed for thirty seconds in slow mode, then the alkaline solution was gradually added. After adding the alkaline solution, the PVA fibers were mixed in the mixture and the mixing operation was carried out for four minutes in the fast mode. After completing the mixing process, the matrix was poured into the standard molds (e.g., 50 mm cubes for compressive strength test) and left in the molds until the samples hardened. Upon hardening, the samples were demolded and wrapped in three folds of thin plastic wrap to avoid excessive loss of water during heat curing. All samples were placed in an oven at 60 °C for 24 h. After completion of the curing time, samples were placed in the lab at a constant temperature of 23 °C. Testing was carried out on the 28th day of the casting.

Experimental and Test Procedures
The workability of the fresh EGC mix was evaluated via a flow table test complying with the requirements of ASTM C1437 [44]. A truncated cone (70 × 100 × 60 mm) was placed in the middle of the flow table, subsequently filled with EGC paste and lifted to allow the paste to flow. After removing the cone, the paste was subjected to 25 jerks in 15 s using mechanical rising and dropping of the flow table. The slump flow is expressed as a percentage change in the flow diameter of the paste and the original base diameter. The same procedure was used to assess the workability of ECC mixes [45].
A total of three samples from all mix designs were tested against each property. A compressive strength test was performed applying the ASTM C109 standard test procedure. The samples were subjected to a gradually applied load; the load was increased at a rate of 3.0 kN/s until complete failure. An elastic modulus test was carried out on cylindrical samples (150 × 300 mm) as per the requirements given in ASTM C 469-14. Flexural strength, flexural toughness, and ductility index were estimated by conducting a three-point loading test on the prismatic samples; the experimental arrangement and sample details are given in Figure 3a, while Figure 3b shows the typical loaddeflection curve of a ductile material. The load was applied under deflection control mode; the deflection was increased at a rate of 0.15 mm/s using a universal testing machine (UTM) with a 100 kN capacity. A linear variable differential transformer (LVDT) was attached to the sample for measuring the midspan deflection. Flexural strength is the ability of a material to resist bending failure and is evaluated according to the procedure described in ASTM C1609 [46]. Flexural toughness is calculated by estimating the area under the load-deflection curve (Figure 3b), according to ASTM C1609, [46]. The ductility index is the measure of the ductility of a given material. It is defined as the ratio between deflection at peak load to the deflection at which the first crack appears on the surface of the sample. The ductility index was calculated according to the procedure proposed by Nematollahi et al. (2014) [38], using Equation (2): where and represent deflection at peak load and deflection at the first crack (Figure 3b). The drying shrinkage test was performed according to the procedure given in ASTM C157-17.

Results and Discussion
The experimental test results for the EGC properties (slump flow, compressive strength, elastic modulus, flexural strength, flexural toughness, ductility index, and drying shrinkage) of all mixtures are given in Table 7. The obtained results were analyzed for the effect of sand to fly ash (S/FA) and water to geopolymer (W/GP) solids on each EGC property using Design-Expert ® software by developing three-dimensional response surface diagrams, contour diagrams, and statistical models as given in subsections. Finally, model validation and optimization of S/FA and W/GS were performed.

Effect of S/FA and W/GP Solids on the Slump-Flow (Workability) of EGC
Workability is a measure of how easy or difficult it is for a mixture to place and handle, usually expressed as the slump flow. The effect of W/GP solids and S/FA on the slump flow of EGC is illustrated with the help of a 3D response surface diagram, while the intervals of compressive strength are shown in the contour diagram as given in Figure 4. A gradual increase in composite flow was observed with the increment in W/GP solids. However, a contrary behavior was observed for S/FA which might have been because higher water content tends to reduce the internal particle friction. The spherical particle shape of the fly ash also induces the roller effect and encourages the flow of the composite. In this study, a maximum slump flow of 182% was observed when EGC synthesizes with W/GP solids and S/FA of 0.26, 0, respectively.  Furthermore, the lowest level of slump flow (85%) was observed in W/GP solids of 0.22 and S/FA of 0.8. Intervals of composite flow are given in the contour diagram, which can be used to predict the values of W/GP solids and S/FA for the desired value of composite flow. Second-order models developed by the RSM technique can be accurately utilized for the prediction and optimization of responses over first-order models [47]. In a study, Mohammed et al. (2012) [48] developed a secondorder RSM model for concrete with paper mill residuals. They concluded that the developed model could accurately predict the slump flow of the concrete for two variable systems, i.e., water to cement ratio and the residual content of the paper mill. Therefore, the second-order RSM model representing slump flow was developed to correlate between W/GP solids and S/FA (Equation (3) where and represents W/GP solids and S/FA, respectively.

Effect of W/GP Solids and S/FA on the Compressive Strength of EGC
Compressive strength is a primary mechanical property because it provides useful information for the structural design of the concrete members. From the 3D response surface diagram (Figure 5), the gradual improvement in compressive strength was observed with the reduction in W/GP solids from 0.26 to 0.244. This could have happened due to the low porosity obtained using a smaller water content in the mix [49]. In another study, Chung et al. (2017) [50] observed a similar behavior of compressive strength at a lower water to binder ratio in ECC. However, further reduction in W/GP solids (e.g., 0.22) was responsible for the sudden drop in workability, leaving macropores in EGC, and therefore, the compressive strength was abruptly decreased for the said scenario.  Additionally, an increase in S/FA caused improvements in compressive strength; a similar trend is mentioned in the available EGC literature [30]. Moreover, the optimum level of compressive strength (77.5 MPa) was observed in W/GP solids and S/FA ash of 0.22 and 0.4, respectively. Another study reported the optimized effect of fly ash and metakaolin in the concrete mixture using the RSM technique [51]. They concluded that ANOVA models developed by RSM could be accurately utilized to predict the compressive strength of concrete; therefore, in this study, the prediction model for compressive strength of EGC is developed and given in Equation (4) where and represent W/GP solids and S/FA, respectively.

Effect of W/GP Solids and Sand/fly Ash on the Modulus of Elasticity of EGC
The modulus of elasticity is an important mechanical property that describes the tendency of a material to deform elastically [52]. Figure 6 depicts the response surface and the contour for elastic modulus versus W/GP solids and S/FA. A reduction in composite elastic modulus was observed with the increment in W/GP solids; this was also reported in the literature [30]. Furthermore, the increase in S/FA tends to improve the composite elastic modulus; a similar trend with the addition of sand content was observed in conventional ECC [53]. Similarly, another research study presented an RSM prediction model for the elastic modulus of rubbercrete having crumb rubber content and water to cement ratio (W/C) as two variables of the system [54]. However, the elastic modulus of EGC can be predicted using Equation (5) where and represents W/GP solids and S/FA, respectively.

Flexural Properties
It is a well-understood fact that the flexural performance of the conventional ECC member reflects its tensile ductility; additionally, there is a good correlation between tensile strain capacity and flexural deflection capacity [55,56]. Hence, it is essential to understand the flexural characteristics of EGC for the evaluation of its performance. The deflection-hardening property is a result of the unique interaction between the fiber and EGC matrix. Furthermore, during the application of load, a crack is initiated from the defected or weak part of the beam sample. The fiber-bridging phenomenon resists the opening of this initial crack. Hence, the load is distributed to the nearest weak area of the crack, which results in another crack formation. The whole process ends with the formation of multiple microcracks at the moment zone of the beam sample, which allow it to undergo a plastic deformation with a deflection-hardening behavior [55]. As matrix-related parameters influence flexural properties, subsections focus on the effect of W/GP solids and S/FA on the flexural properties of EGC. Flexural strength versus midspan deflection curves for all mix designs is given in Figure 7. All the mixes have exhibited deflection-hardening behavior.

Effect of W/GP Solids and S/FA on Flexural Strength of EGC
Fibre-matrix interaction strongly influenced the flexural strength of fiber-reinforced composites [54]. It has been observed in previous research that increase in water content is responsible for weak fiber-matrix interaction [21]. Hence, a reduction in flexural strength was noticed with the rise in W/GP solids (Figure 8). This finding is based on the behavior observed for the water to binder ratio of ECC [57].  Furthermore, contrary to the effect of W/GP solids, improvement in flexural strength was observed with the increment in S/FA. The same behavior of sand dosage on ultimate strength properties was observed in previous studies of fiber-reinforced geopolymer concrete [30]. However, flexural strength can be predicted using the statistical model (Equation (6) where and represent W/GP solids and S/FA, respectively.

Effect of W/GP Solids and S/FA on Flexural Toughness of EGC
Flexural toughness (FT) is an indication of energy absorption capacity to be studied for the evaluation of the structural performance of EGC under hazardous dynamic loading conditions, such as earthquakes or blast. Therefore, flexural toughness provides useful information for the structural designs of such applications. As given in Figure 9, improvement in FT was observed with the initial increment in S/FA; however, a further increase in S/FA resulted in lower FT. This could have happened because of an alteration in the fiber-matrix bond at different S/FA levels. Previous studies have revealed the same findings for the sand to binder ratio of ECC [58].
Furthermore, the increase in the water to binder ratio caused a weak fiber-matrix interfacial transition zone [59] and loosened the frictional bond strength [60]. Therefore, FT was reduced with an increment in the W/GS ratio. Additionally, the optimum level of flexural toughness was observed when EGC synthesizes with W/GP solids and S/FA of 0.22 and 0.4, respectively. However, flexural toughness can be estimated for any value of W/GP solids and sand/fly ash using Equation (7):

Effect of W/GP Solids and S/FA on Ductility Index of EGC
The ductility index (DI) is a quantitative measure of the ductile nature of any material. It is a ratio of deflection at peak load ( ) and deflection at first crack load ( ) [58]. Figure 10 depicts the ductility index value with the variation of W/GP solids and S/FA. A slight improvement in DI was observed with the rise in W/GP solids up to a certain level. However, the maximum value of DI was observed at W/GP of 0.24, which is a moderate level. Furthermore, a decrease in DI was found with the increase in S/FA because the higher dosage of sand tends to increase matrix fracture toughness (Km), which causes the reduction of ductility of a material [61]. However, the DI can be predicted using Equation (8):

Effect of W/GP solids and S/FA on drying shrinkage of EGC
Reduction in length of the mortar caused by the loss of water during the drying process is known as drying shrinkage [62]. ECC does not utilize coarse aggregates; therefore, shrinkage is a significant concern [63]. The same is valid for EGC. The effect of W/GP solids and S/FA on the 28 days drying shrinkage of fiber-reinforced geopolymer is shown in Figure 11. Drying shrinkage was increased with the increment in W/GP solids. The same phenomenon for drying shrinkage was observed with the use of a high amount of water content in ECC [64].
(a) 3D response surface (b) contour diagram Furthermore, the dosage of sand reduced the value of drying shrinkage. These results followed the same pattern as that observed in ECC [65]. The ANOVA prediction model for drying shrinkage is given in Equation (9) where and represent W/GP solids and S/FA, respectively.

Validation of ANOVA Models
Analysis of variance (ANOVA) results of developed RSM models have been analyzed and statistically validated for their practical usage in the future. As given in Table 8, the relative change in adjusted R 2 and predicted R 2 of less than 0.2, a higher value of adequate precision (>4), and the low standard deviation of each response are the indicators of model acceptance [3]. Furthermore, the ANOVA result for all prediction models is shown in Table 9. Prob > F is less than 0.05, and the higher level of F-value is the indicator of the statistically significant model with a 95% confidence level [66]. Therefore, all models were statistically significant and could be used to estimate the response value for any given value of factors. To avoid repetition, plots such as predicted versus actual, normal probability, and perturbation were drawn for only one response (compressive strength). Predicted and actual values are close enough and have a good agreement with the plot, almost coinciding with a 45-degree straight line ( Figure 12). The straight-line coincident with data points of responses indicates normal distribution, as shown in Figure 13. To better understand the effect of independent variables (S/FA and W/GP solids) on compressive strength at a particular point, the perturbation plot is given in Figure 14. Steep slope and curvature indicate that the compressive strength is sensitive to both factors (W/GP solids and S/FA).

Optimization
Usually, the optimum solution of individual responses is located in a different area, and it is tough to find an optimal condition that fulfills the targeted goals of all responses concurrently. Therefore, a multiobjective optimization is utilized in this study to seek dosages of water/geopolymer solids (W/GP solids) and sand/fly ash (S/FA) satisfying the desired conditions of all responses at the same time. The RSM optimization was performed using the global desirability function, as given in Equation (10) [67].
Desirability (D) of the optimal solution is the geometric average of the desirability of all independent as well as dependent variables, where is the number of factors and responses involved in the optimization process. In this study, two factors, namely, W/GP solids and S/FA, in addition to seven responses named slump flow, compressive strength, flexural strength, elastic modulus, flexural toughness, drying shrinkage, and ductility index, were optimized simultaneously, where represents individual desirability function and ranges from 0 to 1 (0 for non-desired and 1 for desired). Furthermore, the importance of each factor or response is given as , which ranges from 1 to 5, where 1 stands for the least important, and 5 represents the most important.
A desirability value close to zero indicates a potential undesirable response and desirability close to one shows that the optimal solution is supposed to achieve the desired goal or target. A goal for responses (dependent variables) can be one of the given choices, i.e., "none", "maximize", "minimize", "target" or "in range"; however, the desired goal for factor includes "equal to" instead of the "none" choice.
A multiobjective optimization technique of response surface methodology was applied to get an optimum solution that could satisfy the desired goals of all responses concurrently. Definitions of factors and responses in multiobjective optimization problem are given in Table 10. The primary aim of this research was to maximize the responses, which include compressive strength, elastic modulus, flexural strength, flexural toughness, and ductility index, while at the same time, drying shrinkage was kept the minimum. A single solution was obtained, which satisfied the specified criteria (Table  11). Ramps of the optimized EGC mix and related optimal responses are shown in Figure 15. After selecting an optimal solution, a validation study was carried out to find the variation in results obtained by performing experiments.

Experimental Validation Study
An experimental program was designed to determine the precision of the RSM models and the consistency of the optimized condition. Mixing, casting, curing, and testing of the samples followed the same procedure as that given in Sections 2.3 and 2.4. Experimental results were found to be in close agreement with the predicted values of the responses, and the difference between them was less than 5% (Table 12). Figures 16 and 17 are added as the evidence of the ultraductile nature of optimized EGC with multiple-crack formation. By contrast, Figure 18 shows the flexural strength versus the midspan deflection curve of the optimized mix.

Conclusion
In this study, optimal dosages of water/geopolymer (W/GP) solids and sand/fly ash (S/FA) in EGC were determined by performing and analyzing the experiments through response surface methodology (RSM). Based on the findings of the investigations reported in this paper, significant conclusions are drawn as follows.
The influence of factors (W/GP solids and S/FA) on responses (slump flow, compressive strength, elastic modulus, flexural strength, flexural toughness, ductility index, and drying shrinkage) has been effectively observed with 3D response surface diagrams.
Contour diagrams can effectively be used to predict the ranges of factors for the desired EGC property (response). Moreover, ANOVA equations have been developed and statistically validated. These ANOVA equations can be used to predict desired responses for any target value of factors accurately.
The optimum dosages of S/FA and W/GP solids in EGC obtained using the multiobjective optimization technique of RSM are 0.341701 and 0.225184, respectively.
The RSM optimization technique reduces the design time and efficiently analyzes the experimental data as well as improves the performance and reliability of the existing process and product.
Author Contributions: Z.M. performed the experimental investigation, carried out data analysis, and prepared the manuscript draft. S.N. supervised the overall research program, reviewed the results, and edited the manuscript. All authors have read and agreed to the published version of the manuscript.