Next Article in Journal
A Novel Hydroforming Process by Combining Internal and External Pressures for High-Strength Steel Wheel Rims
Previous Article in Journal
Photoelasticity for Stress Concentration Analysis in Dentistry and Medicine
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Investigation and Optimisation of the Rheological Properties of Magnesium Potassium Phosphate Cement with Response Surface Methodology

1
College of Materials Science and Engineering, Chongqing University, 174 Shazheng Street, Shapingba, Chongqing 400044, China
2
School of Architecture and Planning, Yunnan University, Kunming 650051, China
3
Department of Civil, Environmental and Geomatic Engineering, University College London, Gower Street, London WC1E 6BT, UK
4
School of Civil Engineering, University of Leeds, Woodhouse Lane, Leeds LS2 9JT, UK
*
Authors to whom correspondence should be addressed.
Materials 2022, 15(19), 6815; https://doi.org/10.3390/ma15196815
Submission received: 12 May 2022 / Revised: 12 September 2022 / Accepted: 21 September 2022 / Published: 30 September 2022
(This article belongs to the Section Construction and Building Materials)

Abstract

:
Magnesium phosphate cement (MPC) is a promising alternative cement. However, the rheological property of this new binder is still to be explored. In this study, Response Surface Methodology (RSM) was adopted with Central Composite Design (CCD) to establish mathematical models describing the rheological characteristics of MPC in terms of initial mini slump (Y1), mini-slump loss (Y2), yield stress (Y3) and plastic viscosity (Y4), as a function of three independent variables, namely, water-to-solid ratio (W/S ratio, X1), MgO to MKP ratio (M/P ratio, X2) and borax dosage (X3). The results show that the M/P ratio and borax dosage could significantly affect the yield stress and mini-slump loss of MPC, while the W/S ratio was the significant coefficient influencing plastic viscosity and initial mini slump. The numerical optimised values of X1, X2 and X3 were 0.280, 7.528 and 0.170, respectively, and an MPC paste with desirable rheological characteristics (Y1 161.858 mm, Y2 11.282, Y3 0.680 Pa, Y4 0.263 Pa·s) with the highest desirability of 0.867 can be obtained.

1. Introduction

Magnesium phosphate cement (MPC), also known as chemically bonded phosphate ceramic, is an alternative clinker-free binder which has attracted increasing attention worldwide [1,2,3,4,5]. Compared to the traditional cement [6], MPC possesses many superior properties, such as super-fast setting, high early strength, strong bonding, low drying shrinkage and better durability [1,2,3,4,5]. These properties make MPC remarkably popular in the fast repair, strengthening and rehabilitation of concrete infrastructures such as pavements, highways and airport runways [7,8]. In recent decades, MPC has also demonstrated huge potential for stabilising/solidifying heavy metals, immobilising nuclear wastes, as well as 3D-printing [9,10]. As a promising alternative cement, the hydration and hardening mechanisms of MPC have been extensively studied in the past [11,12,13,14]. In general, the hydration of MPC is considered to be a through-solution acid-base exothermic reaction between dead-burnt magnesia (MgO) and an acid phosphate source in the presence of water to form magnesium phosphate hydrates [11,12,13,14]. Although monoammonium phosphate (NH4H2PO4, MAP) has widely been employed as the source of phosphate in earlier studies in the literature, monopotassium phosphate (KH2PO4, MKP) has increasingly been employed in recent years owing to its lower solubility than ammonium phosphates, which can slow down the violent reaction and reduce the heat released. Additionally, using MKP also has the advantage of avoiding the emission of unpleasant ammonium gas [14]. The main hydration product of the MgO–KH2PO4 system is MgKPO4·6H2O (k-struvite), and its composition and structure are reported to be analogous to those of NH4MgPO4·6H2O (struvite) formed in the MgO–NH4H2PO4 system [5,14,15]. Different from the poorly ordered calcium silicate hydrate (C-S-H) gel formed in PC, the magnesium phosphate hydrates formed in MPC are mainly crystalline phases, and thus, its volume stability is believed to be better than PC-based materials [16,17].
Although the fundamentals regarding MPC—such as its hydration mechanisms, the properties of the hydrates, factors affecting the properties, as well as its potential applications—have been extensively investigated [13,14,18,19], its rheological properties are rarely reported and largely unknown. Whilst the rheological properties of conventional PC-based cementitious materials have been well established and commonly described using a Bingham model defined by yield stress (τ) and plastic viscosity (μ) [20], studying the rheological properties of MPC could be rather challenging. This could partly be due to its super-fast setting, making it extremely difficult to characterise its fresh properties, such as rheology. Nonetheless, to efficiently utilise MPC-based materials in industrial applications, their rheological behaviour needs to be understood and optimised. This will also significantly benefit the transporting, placing, and finishing processes in engineering practice. Apart from its super-fast setting, the two main constituents of MPC could also cause complex issues to the fresh properties of MPC. This is because, even though the water-to-solid (W/S) ratio is widely established and expected to affect the rheology of MPC system, the MgO to MKP (M/P) ratio can also cause uncertainties to the properties of MPC. For example, it has been shown that the compressive strength of MPC is not a simple monotonic function of M/P. Instead, at a given W/S, there is always an optimal M/P for the highest compressive strength [8], which is different from what has been widely anticipated, i.e., the higher the M/P, the higher the compressive strength. In terms of rheology, whether such kinds of complex also exist is still unknown in the literature. Furthermore, except for the MgO and phosphate salt, a ‘third’ part, termed as retarder, is usually incorporated in the manufacture of MPC to modulate the setting time. This retarding agent, normally boric acid (H3BO3) or borate (Na2B4O7·5H2O), could introduce further uncertainty and complexity to the fresh properties of MPC [12,21]. Due to the above issues, performance-based design approaches have been widely adopted by different researchers to obtain the desired strength properties. Although this trial-error based method has also been extensively used in the laboratory to formulate the MPC with targeted fresh properties, it is not only time-consuming, but also resource-intensive. Moreover, even though this performance-based approach could allow a mix with good performance to be identified, it cannot ensure a mix with the optimal performance being always achieved. Therefore, Response Surface Methodology (RSM) is attempted as a potential mathematical tool in this paper to resolve the above complexes in order to optimise the rheological properties of MPC based on systematic numerical analyses.
Response Surface Methodology (RSM) is a mathematical and statistical approach using empirical modelling to explore the relationships between independent variables (explanatory factors) and the responses of interests, as well as the interactions between the independent variables [22]. Compared to the standard experimental methods, RSM makes a research strategy systematic and efficient to correlate the interactive effects of numerous parameters simultaneously with the help of statistical procedures. The Central Composite Design (CCD) is the most used surface response statistical design method, which runs much fewer experiments than a full factorial design but provides almost the same information [23]. Hence, the RSM jointly used with the statistic design method (e.g., CCD) is recognised as a significant approach in the design and optimisation process, which has also recently been effectively employed to investigate the interactive effects of independent variables in many areas [24,25,26]. In the field of cement and concrete research, RSM demonstrated great potential for experimental optimisation, particularly in optimising the rheological parameters, hydration and mechanical properties of cementitious materials [27,28,29,30]. Recently, Hou et al. [31] presented their research on optimising the setting time and early-age compressive strength of MPC cement with RSM methodology for rapid repair applications. Using water-to-binder ratio, magnesia-to-phosphate ratio and borax content as the variables, and the setting time and early-age compressive strength as the responses in RSM, a statistical model was successfully developed by Hou et al. [31], showing a good potential of applying RSM for the design and optimisation of MPC mix. However, the rheological parameters (e.g., mini slump, mini-slump loss, yield stress and plastic viscosity) of MPC were not investigated in their study.
Building upon the aforementioned information, it is anticipated that RSM-CCD methodology could be used to explore the rheological property of this new MPC materials. Hence, in this study, the RSM methodology was adopted to investigate the rheological behaviour of MPC material and also to optimise the mix proportion in terms of water/solid ratio (W/S), MgO/MKP ratio (M/P) and borax dosage, with initial mini slump, mini-slump loss, yield stress (τ, Pa) and plastic viscosity (μ, Pa·s) considered as responses. The CCD design was, therefore, conducted to develop a three variables (factors) (n = 3) experiment matrix with 20 runs in total. In addition, the analysis of variance (ANOVA) was conducted to assess the significance and adequacy of the regression models attained. Finally, to supplement the mathematical analyses, the early-age (i.e., 1 d and 7 d) compressive strength of the MPC pastes was investigated to verify its application in real engineering practices.

2. Materials and Methods

2.1. Materials

Magnesium oxide (MgO) provided by Richard Baker Harrison Ltd., Manchester, UK, was a Dead Burned Magnesite (DBM) calcined at about 1750 °C with a purity of 90% and a mesh size of 200 μm. Monopotassium phosphate (KH2PO4, MKP) was a food-grade MKP with a specified purity >99%, from Prayon UK. Sodium tetraborate decahydrate (borax, Na2B4O7·10H2O), supplied by Sigma-Aldrich US, was employed as the retarder. It was ACS reagent grade with assay ≥99.5% and a pH of 9.15–9.20.

2.2. Design and Sample Preparation

The work process contains three steps: (1) selecting targeted performance parameters and key influencing factors; (2) building and fitting the statistical models for each performance parameter; (3) optimising the mix proportion to achieve the aimed performance. The overall framework for this study is illustrated in Figure 1.
A typical CCD design involves 2n factorial points, 2n axial points, and nc centre points, in which n represents the number of factors/variables, and the nc refers to the number of repetitions of the centre point. The repetition of the central point is to estimate the experimental internal error and, thus, to improve the reliability of the models. In this study, a 23 CCD design with three variables (n = 3) and two extreme levels (coded as +α and −α) was applied to develop the mathematical equation and quantify the rheological properties of MPC in terms of initial mini slump, mini-slump loss, yield stress and plastic viscosity. The three variables investigated in this study were water/solid ratio (W/S) (X1), MgO/MKP ratio (M/P) (X2) and borax dosage (X3) (by the weight of MgO). A total of 20 runs, including eight (2n = 23) factorial points, six (2n = 2 × 3) axial points and six (replicates) centre points, were carried out. The CCD design and data analysis associated with the analysis of variance (ANOVA) and RSM optimisation were conducted using Design Expert 10 (Stat-Ease, Minneapolis, MN, USA) software. Table 1 below reports the actual values of the factors, and Figure 2 visually illustrates the distribution of the coded points for the designed matrix. It should be noted that the ranges/levels of the variables were determined based on the results from our previous work [32]. In total, 20 MPC pastes were designed by the CCD method and the details of the coded value of mixing portions are given in Table 2.
All 20 MPC pastes, as shown in Table 2, were fabricated in a high-shear mixer. The borax was first mixed with water and followed by adding MKP and MgO. The total mixing time was 3 min.

2.3. Test Methods

2.3.1. Rheological Properties Test

Immediately after mixing, the pastes were transferred into the container (92 mm diameter) of a modified and calibrated rheometer, which uses a helical impeller to establish the relationship between shear stress and shear rate. The results were recorded by the software and the yield stress and plastic viscosity were fitted by the Bingham model. The mini-slump test was conducted with a conical mould with a lower inner diameter of 38.1 mm, an upper inner diameter of 19 mm and a height of 57.1 mm. The mini slump was determined at 15 min and 45 min after mixing. The mini-slump loss was then calculated as the difference between the values of mini slump at 15 min and 45 min.

2.3.2. Compressive Strength Test

The 25 × 25 × 25 mm3 PVC moulds were used to cast cubes to determine the compressive strength of the paste. Immediately after mixing, the moulds were filled with the paste and a steel scraper was used to finish the specimen surface. The moulds were covered with water-saturated hessian and stored in a constant temperature room at 20 ± 1 °C for 24 h. The specimens were then demoulded, and each cube was wrapped with a water-saturated hessian before being sealed in an air-tight plastic bag. The specimens were again stored in the same constant-temperature room until being tested.

2.4. Data Analysis

The data were processed by Design Expert software. The analysis of variance (ANOVA) with Fisher’s F-test was conducted to check the significance and adequacy of the models established. The desirability function (DF) approach was used to establish the optimum criteria based on multi-variables. The general procedure of this approach is to convert each response Y k into an individual desirability function as   d k = f Y k   0     d k 1 . If the response ( Y k ) meets the target of the optimisation goal,   d k   is then valued as 1; while the response is beyond the acceptable limit, then d k is 0.

3. Results and Discussion

3.1. Model Adequacy Analysis

The experimental data in terms of the four responses (i.e., Y1, Y2, Y3 and Y4) obtained from all 20 mixes are presented in Table 3 and analysed by Design Expert 10. The best-fitting surface response model for describing the mini slump (Y1) is generated, which is suggested to be a linear relation, whilst the models of the other three responses, i.e., mini-slump loss (Y2), yield stress (Y3), and plastic viscosity (Y4), are two-factor interactions (2FI). The regression equations in the coded value attained correlating to the responses and the three independent variables are shown in Equations (1)–(4) for initial mini slump, mini-slump loss, yield stress and plastic viscosity, respectively. The analysis of variance (ANOVA) with Fisher’s F-test was then conducted to assess the significance and adequacy of these models, which are reported and discussed in separate sections (Section 3.1.1, Section 3.1.2, Section 3.1.3, Section 3.1.4 and Section 3.1.5) below. Since the repeatability of the responses at the central points was used for estimating the error of the models, the analysis of the results of the four responses at central points was also conducted, and the results are presented in Table 4.
Y1 = 153.60 + 9.48X1 + 0.14X2 − 1.24X3
Y2 = 38.75 − 9.84X1 + 28.68X2 − 13.55X3 + 3.87X1X2 − 5.62X1X3 − 23.88X2X3
Y3 = 9.61 − 9.67X1 + 11.60X2 − 10.41X3 − 14.25X1X2 + 14.33X1X3 − 15.23X2X3
Y4 = 0.73 − 0.43X1 − 0.04X2 + 0.20X3 + 0.24X1X2 − 0.37X1X3 + 0.52X2X3

3.1.1. Analysis of Variance (ANOVA) of the Response Model of Initial Mini-Slump

The analysis of variance of the model of initial mini slump (Equation (1)) is shown in Table 5. The obtained F-value of 9.67 and p-value of 0.0007 (<0.05) implies that the regression model was highly significant at a 5% significance level, and there was only a 0.07% chance that an F-value this large could occur due to noise. The goodness of fit of the model was assessed by the coefficient of determination R2 and equalled to 0.6446, which indicates that about 64.46% of the variation in this system was attributed to the independent variables, and only about 35% of the variation could not be explained by the model. In addition, the adequate precision value measures the signal-to-noise ratio and a value greater than 4 is desirable. The value obtained from the current model was 10.8699, suggesting that the adequate signal was obtained, and this model can be used to navigate the design space. Moreover, the significance of each of the coefficients can be checked by the p-value, as listed in Table 5, and p < 0.05 could suggest that the term is significant. From our results, evidently, the X1 is the significant model term. This suggests that the factor W/S (X1) is the most significant term for the linear model of mini slump, which significantly affects the initial mini slump of the MPC paste.

3.1.2. Model Describing Mini-Slump Loss

Similarly, the ANOVA analysis with Fisher’s F-test was conducted to check the significance and adequacy of the model of mini-slump loss (Equation (2)), which is presented in Table 6. The F-value and p-value of the 2FI model were 5.18 and 0.0063, respectively, revealing that the model was significant. Moreover, the coefficient of determination (goodness of fit) was R2 = 0.7051, which indicates that 70.51% of the variation in this system can be attributed to the independent variables. In addition, the value of adequate precision was 7.6423, which suggests that the model can be validly used in the examined range of experiments. Moreover, by checking the p-value of all the terms, the factors of X2 (0.0011) as well as X2X3 (0.0195) were significant model terms, which suggests that M/P (X2) significantly affects the mini-slump loss of MPC (95% confidence limit), and there are significant interactions between M/P and borax dosage (X2X3). On the other hand, the factor borax dosage (X3) was still within the 90% confidence limit, which also somehow showed the influence on the mini-slump loss of MPC.

3.1.3. Model Describing Yield Stress

The statistical significance and adequacy of the response model of yield stress (Equation (3)) were checked, and the results are shown in Table 7. The response model was significant, as evidenced by the F-value and p-value which were 5.13 and 0.0066, respectively. Moreover, the coefficient of determination R2 equalled to 0.7030 (goodness of fit), indicating that 70.30% of the variation in this system can be attributed to the independent variables, and only about 30% of the variation cannot be explained by the model. The lower value of R2 could be due to the indirect response (yield stress) from the factors. The value of adequate precision obtained from the current model was 9.6694 (greater than 4), suggesting that the adequate signal was obtained and this model is valid to navigate the design space. Moreover, by checking and comparing the p-value listed in Table 7, it is evident that all the factors and interactions, excepting W/S (X1), were identified as significant model terms at a 95% confidence limit. Based on the statistical analysis, the M/P (X2) and borax dosage (X3) could significantly affect the yield stress of the MPC materials, and there are significant interactions between W/S and M/P (X1X2), W/S and borax dosage (X1X3) and M/P and borax dosage (X2X3). Furthermore, the p-value of W/S (X1) was 0.0655, which is still acceptable within the allowance of the 90% confidence limit.

3.1.4. Model Describing Plastic Viscosity

The ANOVA analysis for the model of plastic viscosity (Equation (4)) is shown in Table 8. The current model, with a p-value of 0.0979 (which is less than 0.1), could still be considered as significant at a 10% significance level. The determination coefficient R2 obtained was 0.5151, which suggests that this model can interpret 51.51% of the total variation. Moreover, the adequate signal was retrieved, as evidenced by the adequate precision value of 6.1537. More importantly, the X1 was observed to be the only significant coefficient in the current model, which implies that the W/S ratio played an important role in adjusting the viscosity of the MPC materials.

3.1.5. Summary and Further Improvement

In this section, in addition to the establishment of four prediction models with the analysis factors and their interactions, an ANOVA was carried out, and the influence of the main factors was assessed along with the interaction. For all the models, the p-values were less than 0.1, indicating that these four models are significant at a 10% significance level. Furthermore, based on the analysis of each factor and their interactions, the most significant terms for each model were identified, as discussed in Section 3.2.5. However, two improvements could be further applied: (1) the results indicated strong interactions between the factors affecting mini-slump loss and yield stress, which could be further analysed; (2) the values of R2 of the four models were below 0.85 and the lowest value was 0.5151, indicating that there is some room to improve these models.

3.2. Effect of Variables on the Response of the Model

In this section, based on the statistical analysis of the four predicted models, the effects of different variables and their interactions in each predicted model are discussed.

3.2.1. Effect on Initial Mini-Slump

The 3D response surface plots of the initial mini slump with an interactive relationship between the factors (W/S, M/P, borax dosage) are shown in Figure 3a–c. Based on Equation 1, the initial mini-slump response increased with the increase in the W/S ratio and M/P ratio, whereas it decreased with the increase in borax dosage. Moreover, the higher coefficient of the W/S, which was approximately 67.1 and 7.6 times larger than those of M/P and borax dosage, respectively, indicates that the W/S significantly affected the initial mini slump of MPC. As shown in Figure 3a, the initial mini slump increased with the increase in W/S, whereas the change of the M/P imposed little effect on it. For example, the response of the predicted mini slump was 153.60 mm when the values of W/S and M/P were set as 0.24 and 7, respectively. However, when increasing the W/S to 0.28 while maintaining M/P at 7, the predicted initial mini slump increased to 163.07 mm, while it only increased to 153.73 mm if increasing M/P to 9.97 but keeping the W/S at 0.24. On the other hand, the effects of W/S and borax dosage on the initial mini slump are presented in Figure 3b. Similarly, the initial mini slump increased as the W/S increased towards its high level. In contrast, the initial mini slump decreased along with increasing borax dosage. In addition, as illustrated in Figure 3c, the initial mini slump increased with decreasing borax dosage, while the change in M/P provided little effect on it.

3.2.2. Effect on Mini-Slump Loss

The 3D response surface plots visualizing the effect of the three factors on the mini-slump loss of MPC paste are plotted in Figure 4. Generally, the mini-slump loss response decreased with the increase in W/S and borax dosage, while it increased with the increase in M/P. Moreover, the influences in the order of decreasing magnitude are M/P, borax dosage and W/S. For example, the coefficient of M/P showed an approximately 2.11 times greater effect on the mini-slump loss than the borax dosage (+28.68 vs. −13.55 in Equation (2)). As shown in Figure 4a, while keeping the borax dosage at the middle level (0.16), it is apparent that mini-slump loss decreased with decreasing M/P, whereas it decreased with increasing W/S. Moreover, the combined effects of W/S and borax dosage are shown in Figure 4b. As expected, the mini-slump loss decreased with increasing W/S and borax dosage. Similarly, as shown in Figure 4c, the mini-slump loss decreased with the decrease in M/P and increase in borax dosage. For example, when M/P was 7 and borax dosage was 0.16, the predicted mini-slump loss was 38.75. However, when M/P decreased to 4.03 while maintaining borax dosage at 0.16, the predicted mini-slump loss decreased to 10.23; meanwhile, it also decreased to 25.20 if increasing the borax dosage to 0.17 but keeping the M/P at 7.

3.2.3. Effect on Yield Stress

The 3D response surface plots visualizing the effect of the three factors on the yield stress of MPC paste are plotted as the response surfaces in Figure 5. As shown in Equation (3), the yield stress response decreased with the increase in W/S and borax dosage, whereas it increased with the increase in M/P. Moreover, the influences in the order of decreasing magnitude are M/P ratio, borax dosage and W/S. For example, the coefficient of M/P showed an approximately 1.2 times greater effect on the yield stress than that of W/S (+11.60 vs. −9.67). As shown in Figure 5a, while keeping the borax dosage at the middle level (0.16), it is evident that the yield stress increased with increasing M/P, whereas it decreased with the increase in W/S. When W/S was 0.24 and M/P was 7, the response of the predicted yield stress value was 9.61 Pa. However, when W/S decreased to 0.20 while maintaining M/P at 7, the predicted yield stress value increased to 19.28 Pa, while it increased to 21.21 Pa if increasing M/P to 9.97 but keeping the W/S at 0.24. The combined effects of W/S and borax dosage are shown in Figure 5b. As expected, the yield stress decreased with increasing W/S and borax dosage. Similarly, as shown in Figure 5c, the yield stress increased with the increase in M/P, while it decreased with increasing borax dosage. There could be interactions between the three factors in relation to the yield stress. For example, the proportion of water could not only affect the fluidity, but also change the hydration process of the MPC [33]. Therefore, further studies still need to be carried out to clarify their roles.

3.2.4. Effect on Plastic Viscosity

The 3D response surface plots based on Equation 4 are shown in Figure 6a–c. In general, the plastic viscosity response decreased with the increase in the W/S ratio and M/P ratio, while it increased with the increase in borax dosage. Moreover, it is obvious that the coefficient of the W/S (−0.43) was approximately 9.8 and 1.2 times higher than that of M/P and borax dosage, respectively, indicating that the W/S can significantly affect the plastic viscosity of MPC. The dependence of plastic viscosity on W/S and M/P for the set of borax dosage at its central level (0.16) is shown in Figure 6a. As can be seen, the plastic viscosity decreased with the increase in W/S, whereas the change in M/P imposed little effect on plastic viscosity. For example, the response of the predicted plastic viscosity was 0.73 Pa·s when the value of W/S and M/P were set as 0.24 and 7, respectively. However, when decreasing the W/S to 0.20 while maintaining the M/P at 7, the predicted plastic viscosity increased to 1.16 Pa·s. However, it only increased to 0.77 Pa·s if decreasing the M/P to 4.03 but keeping the W/S at 0.24. This phenomenon may be attributed to the fact that the change in W/S and MgO proportion could affect the cohesion of the fluids. Additionally, the effects of W/S and borax dosage on the plastic viscosity are presented in Figure 6b. The plastic viscosity decreased as the W/S increased towards its high level, while it increased along with increasing borax dosage. As illustrated in Figure 6c, the plastic viscosity increased with increasing borax dosage, while the change in M/P contributed little to plastic viscosity.

3.2.5. Discussion

To compare the effects of each variable on the rheological behaviour of MPC, the coefficients of all regression models are summarised in Table 9. Since the p-value < 0.05 would indicate that the term is significant at 95% confidence, the significant terms for each response in Table 9 are bold.
As shown in Table 9, the W/S ratio mainly affected the initial mini slump and plastic viscosity. The increase in W/S led to the increase in mini slump, but the decrease in plastic viscosity. However, it should be noticed that, although yield stress was not the significant term at 95% confidence, it was still above the 90% confidence level, indicating that it still can be influenced by W/S. In a cementitious system, Van der Waals interactions have been recognised to dominate all colloidal interactions and, therefore, dictate the inter-particle distance [34]. The increased W/S reduces the concentration of the particles which, therefore, enlarges the inter-particle distance. Since the inter-particles force is inversely proportional to the square of inter-particle distance [35], the increased inter-particle distance due to the high W/S can reduce the inter-particle force. This, in turn, can reduce the formation of the solid network through particle–particle bonds [36]. Consequently, higher mini slump and lower yield stress can be achieved. Moreover, W/S is also identified as the most significant factor influencing the plastic viscosity of MPC. This could be attributed to the fact that the plastic viscosity can somehow reflect the flocculation status in the paste and the resistance to the flow, which in turn, depends on the solid volume fraction and the packed density [37]. The alteration of W/S changed the solid volume, consequently changing the plastic viscosity of MPC.
In terms of the M/P ratio, Table 9 demonstrates that it can significantly affect the mini-slump loss and yield stress. The higher M/P ratio shows a higher mini-slump loss and yield stress. It is well known that the mini-slump loss is closely related to the setting of MPC, whilst the yield stress is closely related to the hydration process. Although it is generally accepted that the main hydration product of MKPC is k-struvite (MgKPO4·6H2O), which is formed from the reaction between MgO and KH2PO4 [38] through solution, the different ratios between MgO and KH2PO4 can lead to a different reaction process and products [13]. Specifically, the increased portion of MgO could lead to a higher pH of MPC, which promotes the dissolution of the ions available for reaction [39], and more reaction products could form and precipitate on the particle surface [14]. Consequently, the hydration process could be accelerated and setting time shortened, resulting in a higher mini-slump loss. Furthermore, the produced negatively charged k-struvite or other intermedium products [40,41] from hydration may interact with the positively charged MgO particles [42], which could increase the yield stress by forming the electrostatic attraction force among the particles.
Similarly, borax dosage, another factor which is related to the hydration process, affected the mini-slump loss and yield stress. In contrast to the effect of M/P ratio, the increase in borax dosage reduced both mini-slump loss and yield stress. The addition of borax retarders delayed the setting of MPC [43]. As discussed previously, since the mini-slump loss is linked to the setting of MPC, the increase in the borax dosage reduced the mini-slump loss of MPC by delaying the setting of MPC. Moreover, the addition of borax retards the MPC hydration by formatting the positively charged MgB(OH)4+ complex [43]. Therefore, the electrostatic repulsion between the positively charged complex and MgO could benefit the dispersion of the particles, which can further reduce the yield stress of the MPC paste.
Furthermore, it should be noted that, based on statistical analysis the interaction between any two factors (X1X2, X1X3, X2X3) can still affect the yield stress of the MPC paste. Besides, the interaction between M/P and borax dosage (X2X3) exhibited a greater influence on mini-slump loss. Such interactions lead to a more complex situation in analysing the fresh property of MPC. Therefore, although the influence of each variable on the fresh property of MPC was discussed in detail, the interactions between these variables need to be further investigated.

3.3. Desirability Functions for Numerical Optimisation

Based on the results and the analysis presented in previous sections, the models for all four responses, in terms of initial mini slump, mini-slump loss, yield stress and plastic viscosity (Equations (1)–(4)), were utilised to optimise the best fresh property of MPC paste. Factors including W/S, M/P and borax dosage are involved in the numerical optimisation with the design goal as ‘in range’. On the other hand, the goals in terms of yield stress, plastic viscosity, initial mini slump and mini-slump loss are set as ‘minimised’, ‘in range’, ‘maximised’ and ‘minimised’ to obtain the MPC system with desirable fresh properties. Considering the ease to process the fresh MPC in its future engineering application, the importance of yield stress, initial mini slump and mini-slump loss was set as 5, whereas it was set as 3 (default value) for other terms. The factors and predicted response of the target are listed in Table 10. As shown in the table, the predicted values for the four responses generated from the optimisation are 0.680 Pa (yield stress), 0.263 Pa·s (plastic viscosity), 161.858 mm (initial mini slump) and 11.282 (mini-slump loss), with the highest desirability of 0.867. Accordingly, the optimum conditions of the three variables are 0.280, 7.528 and 0.170 for the W/S, M/P and borax dosage, respectively. Thus, it is anticipated that by following the optimum recipe computed for W/S, M/P and borax dosage, a MPC mix with desirable fresh properties can be obtained. Our paper reported for the first time the optimised rheological parameters of the MPC materials by the statistic RSM-CCD methodology.

3.4. Compressive Strength

To use MPC in practice, its compressive strength should be examined. Considering that MPC develops its strength very quickly at an early age, the 1 d and 7 d compressive strength of the paste cubes were tested, and the results are reported in Table 11. The strength values were not particularly high, which is due to the high W/S ratio used in this study, as reported in Table 1. For all mixes, the compressive strength increased with age (from 1 d to 7 d), which agrees with the results reported in the literature [17]. The development of the strength, which has been widely investigated, is due to the continuous generation of the hydration product i.e., struvite. It should be noted that the strength of the Run 15 is not available, which is due to the break-down of the cubes caused by expansion. It can be noticed that the compressive strength varied with the change of the three factors, i.e., W/S Ratio (X1), M/P Ratio (X2) and borax dosage (X3), indicating the potential effects of these factors on compressive strength. The relationships between the strength and the rheologic parameters, as well as the interactions between the variables and their combined effects on the strength and rheologic properties will be investigated in future.

4. Conclusions

Based on the results obtained, the following conclusions can be drawn:
  • The RSM-CCD methodology was successfully adopted to investigate the rheological behaviour of MPC material and to optimise the mix proportion in terms of W/S, M/P and borax dosage, with initial mini slump, mini-slump loss, yield stress and plastic viscosity considered as responses.
  • The W/S ratio was identified as the significant factor (95% confidence level) affecting the plastic viscosity and the initial mini slump. The increase in W/S led to the decrease in the plastic viscosity, whereas it increased the mini slump. Moreover, the influence on the yield stress could not be ignored, since it remained at a 90% confidence level.
  • The yield stress and mini-slump loss were influenced by the M/P ratio. The increase in the M/P ratio was shown to increase the yield stress and mini-slump loss.
  • Borax dosage clearly affected the yield stress and mini-slump loss of MPC. With the increase in borax dosage, the yield stress and mini-slump loss decreased.
  • The numerical optimisation showed that the best predicted values for the four responses are 0.680 Pa (yield stress), 0.263 Pa·s (plastic viscosity), 161.858 mm (initial mini slump) and 11.282 (mini-slump loss), with the desirability of 0.867.

Author Contributions

Conceptualization, Y.Y., J.Q. and Y.B.; data curation, Y.Y. and J.R.; formal analysis, J.R., K.Y. and D.W.; funding acquisition, J.Q. and Y.B.; investigation, Y.Y.; methodology, Y.Y. and J.R.; project administration, Y.Y.; supervision, J.Q. and Y.B.; validation, Y.Y.; visualization, J.R., K.Y. and D.W.; writing—original draft, Y.Y.; writing—review and editing, J.R. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the UK Engineering and Physical Sciences Research Council (EPSRC) (grant number EP/M003159/1) and the National Natural Science Foundation of China (NSFC) (grant number 51461135003) under the EPSRC-NSFC Collaborative Research Scheme on ‘Sustainable Materials for Infrastructure’.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Data are contained within the article.

Acknowledgments

The authors would like to acknowledge Richard Baker Harrison Ltd. for providing magnesium oxide (MgO) for this research.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Qin, J.H.; Qian, J.S.; You, C.; Fan, Y.R.; Zhen, W. Bond behavior and interfacial micro-characteristics of magnesium phosphate cement onto old concrete substrate. Constr. Build. Mater. 2018, 167, 166–176. [Google Scholar] [CrossRef]
  2. Yang, Q.; Zhu, B.; Wu, X. Characteristics and durability test of magnesium phosphate cement-based material for rapid repair of concrete. Mater. Struct. 2000, 33, 229–234. [Google Scholar] [CrossRef]
  3. Qiao, F.; Chau, C.K.; Li, Z. Property evaluation of magnesium phosphate cement mortar as patch repair material. Constr. Build. Mater. 2010, 24, 695–700. [Google Scholar] [CrossRef]
  4. Sun, W.; Zheng, Y.; Zhou, L.; Song, J.; Bai, Y. A study of the bond behavior of FRP bars in MPC seawater concrete. Adv. Struct. Eng. 2020, 24, 1110–1123. [Google Scholar] [CrossRef]
  5. Ma, H.Y.; Xu, B.; Li, Z. Magnesium potassium phosphate cement paste: Degree of reaction, porosity and pore structure. Cem. Concr. Res. 2014, 65, 96–104. [Google Scholar] [CrossRef]
  6. Mahmoud, A.A.; Elkatatny, S.; Al-Majed, A.; Ramadan, M.A. The use of graphite to improve the stability of Saudi class G oil-well cement against the carbonation process. ACS. Omega 2022, 7, 5764–5773. [Google Scholar] [CrossRef]
  7. Seehra, S.S.; Gupta, S.; Kumar, S. Rapid setting magnesium phosphate cement for quick repair of concrete pavements—Characterisation and durability aspects. Cem. Concr. Res. 1993, 23, 254–266. [Google Scholar] [CrossRef]
  8. Ma, H.Y.; Xu, B.; Liu, J.; Pei, H.; Li, Z. Effects of water content, magnesia-to-phosphate molar ratio and age on pore structure, strength and permeability of magnesium potassium phosphate cement paste. Mater. Des. 2014, 64, 497–502. [Google Scholar] [CrossRef]
  9. Lai, Z.Y.; Wang, H.T.; Hu, Y.; Yan, T.; Lu, Z.Y. Rapid Solidification of Highly Loaded High-Level Liquid Wastes with Magnesium Phosphate Cement. Ceram. Int. 2018, 45, 5050–5057. [Google Scholar] [CrossRef]
  10. Weng, Y.W.; Ruan, S.Q.; Li, M.Y.; Mo, L.W.; Unluer, C.; Tan, M.J.; Qian, S.Z. Feasibility study on sustainable magnesium potassium phosphate cement paste for 3D printing. Constr. Build. Mater. 2019, 221, 595–603. [Google Scholar] [CrossRef]
  11. Abdelrazig, B.; Sharp, J.H.; El-Jazairi, B. The Chemical Composition of Mortars Made from Magnesia-Phosphate Cement. Cem. Concr. Res. 1988, 18, 415–425. [Google Scholar] [CrossRef]
  12. You, C.; Qian, J.S.; Qin, J.H.; Wang, H.; Wang, Q.; Ye, Z. Effect of early hydration temperature on hydration product and strength development of magnesium phosphate cement (MPC). Cem. Concr. Res. 2015, 78, 179–189. [Google Scholar] [CrossRef]
  13. Xu, B.W.; Lothenbach, B.; Leemann, A.; Winnefeld, F. Reaction mechanism of magnesium potassium phosphate cement with high magnesium-to-phosphate ratio. Cem. Concr. Res. 2018, 108, 140–151. [Google Scholar] [CrossRef]
  14. Xu, B.W.; Winnefeld, F.; Kaufmann, J.; Lothenbach, B. Influence of magnesium-to-phosphate ratio and water-to-cement ratio on hydration and properties of magnesium potassium phosphate cements. Cem. Concr. Res. 2019, 123, 105781. [Google Scholar] [CrossRef]
  15. Xu, B.W.; Ma, H.; Li, Z. Influence of magnesia-to-phosphate molar ratio on microstructures, mechanical properties and thermal conductivity of magnesium potassium phosphate cement paste with large water-to-solid ratio. Cem. Concr. Res. 2015, 68, 1–9. [Google Scholar] [CrossRef]
  16. Qin, Z.H.; Ma, C.; Zheng, Z.Q.; Long, G.C.; Chen, B. Effects of metakaolin on properties and microstructure of magnesium phosphate cement. Constr. Build. Mater. 2020, 234, 117353. [Google Scholar] [CrossRef]
  17. Li, Y.; Chen, B. Factors that affect the properties of magnesium phosphate cement. Constr. Build. Mater. 2013, 47, 977–983. [Google Scholar] [CrossRef]
  18. Jia, X.W.; Li, J.M.; Wang, P.; Qian, J.S.; Tang, M.H. Preparation and mechanical properties of magnesium phosphate cement for rapid construction repair in ice and snow. Constr. Build. Mater. 2019, 229, 116927.1–116927.13. [Google Scholar] [CrossRef]
  19. Wang, D.Q.; Yue, Y.F.; Mi, T.W.; Yang, S.Y.; Bai, Y. Effect of magnesia-to-phosphate ratio on the passivation of mild steel in magnesium potassium phosphate cement. Corros. Sci. 2020, 174, 108848. [Google Scholar] [CrossRef]
  20. Banfill, P.F.G.; Teixeira, M.A.O.M.; Craik, R.J.M. Rheology and vibration of fresh concrete: Predicting the radius of action of poker vibrators from wave propagation. Cem. Concr. Res. 2011, 41, 932–941. [Google Scholar] [CrossRef]
  21. Ribeiro, D.V.; Paula, G.R.; Morelli, M.R. Effect of boric acid content on the properties of magnesium phosphate cement. Constr. Build. Mater. 2019, 214, 557–564. [Google Scholar] [CrossRef]
  22. Bezerra, M.A.; Santelli, R.E.; Oliveira, E.P.; Villar, L.S.; Escaleira, L.A. Response surface methodology (RSM) as a tool for optimization in analytical chemistry. Talanta 2008, 76, 965–977. [Google Scholar] [CrossRef] [PubMed]
  23. Ren, J.; Fang, Y.H.; Ma, Q.M.; Tan, H.B.; Wang, X.F. Effect of storage condition on basic performance of polycarboxylate superplasticiser system incorporated sodium gluconate. Constr. Build. Mater. 2019, 223, 852–862. [Google Scholar] [CrossRef]
  24. Ren, J. Superplasticiser for NaOH-Activated Slag: Competition and Instability between Superplasticiser and Alkali-Activator. Ph.D. Thesis, University College London (UCL), London, UK, 2016. [Google Scholar]
  25. Mohammed, S.J.M.; Luti, K.J.K. Response surface methodology: A review on its applications and challenges in microbial cultures. Mater. Today-Proc. 2021, 42, 2277–2284. [Google Scholar] [CrossRef]
  26. Yolmeh, M.; Jafari, S.M. Applications of Response Surface Methodology in the Food Industry Processes. Food. Bioprocess. Tech. 2017, 10, 413–433. [Google Scholar] [CrossRef]
  27. Zl, A.; Dlab, C.; Xga, B. Multi-objective optimization of gap-graded cement paste blended with supplementary cementitious materials using response surface methodology. Constr. Build. Mater. 2020, 248, 118552. [Google Scholar] [CrossRef]
  28. Ghafari, E.; Costa, H.; Júlio, E. Statistical mixture design approach for eco-efficient UHPC. Cem. Concr. Compos. 2015, 55, 17–25. [Google Scholar] [CrossRef]
  29. Sonebi, M.; Abdalqader, A.; Fayyad, T.; Perrot, A.; Bai, Y. Optimisation of rheological parameters, induced bleeding, permeability and mechanical properties of supersulfated cement grouts. Constr. Build. Mater. 2020, 262, 120078. [Google Scholar] [CrossRef]
  30. Sonebi, M.; Lachemi, M.; Hossain, K. Optimisation of rheological parameters and mechanical properties of superplasticised cement grouts containing metakaolin and viscosity modifying admixture. Constr. Build. Mater. 2013, 38, 126–138. [Google Scholar] [CrossRef]
  31. Hou, D.; Chen, D.; Wang, X.; Wu, D.; Yu, R. RSM-based modelling and optimization of magnesium phosphate cement-based rapid-repair materials. Constr. Build. Mater. 2020, 263, 120190. [Google Scholar] [CrossRef]
  32. Yue, Y.F.; Hu, W.Y.; You, C.; Qian, J.H.; McCague, C.; Jin, F.; Al-Tabbaa, A.; Mo, L.W.; Deng, M.; Bai, Y. A Statistical Investigation of the Rheological Properties of Magnesium Phosphate Cement. In Proceedings of the International Workshop on Innovation in Low-Carbon Cement & Concrete Technology, London, UK, 21–24 September 2016. [Google Scholar]
  33. Lahalle, H.; Coumes, C.C.D.; Mercier, C.; Lambertin, D.; Gauffinet, S. Influence of the w/c ratio on the hydration process of a magnesium phosphate cement and on its retardation by boric acid. Cem. Concr. Res. 2018, 109, 159–174. [Google Scholar] [CrossRef]
  34. Robert, J.; Flatt. Dispersion forces in cement suspensions. Cem. Concr. Res. 2004, 34, 399–408. [Google Scholar] [CrossRef]
  35. Perrot, A.; Lecompte, T.; Khelifi, H.; Brumaud, C.; Hot, J.; Roussel, N. Yield stress and bleeding of fresh cement pastes. Cem. Concr. Res. 2012, 42, 937–944. [Google Scholar] [CrossRef]
  36. Rubio-Hernández, F.; Adarve-Castro, A.; Velázquez-Navarro, J.; Páez-Flor, N.; Delgado-García, R. Influence of water/cement ratio, and type and concentration of chemical additives on the static and dynamic yield stresses of Portland cement paste. Constr. Build. Mater. 2020, 235, 117744. [Google Scholar] [CrossRef]
  37. Struble, L.J.; Lei, W.G. Rheological changes associated with setting of cement paste. Adv. Cem. Based. Mater. 1995, 2, 224–230. [Google Scholar] [CrossRef]
  38. Lu, Z.; Hou, D.; Ma, H.; Fan, T.; Li, Z. Effects of graphene oxide on the properties and microstructures of the magnesium potassium phosphate cement paste. Constr. Build. Mater. 2016, 119, 107–112. [Google Scholar] [CrossRef]
  39. Le Rouzic, M.; Chaussadent, T.; Platret, G.; Stefan, L. Mechanisms of k-struvite formation in magnesium phosphate cements. Cem. Concr. Res. 2017, 91, 117–122. [Google Scholar] [CrossRef]
  40. Zhang, X.; Hu, J.; Spanjers, H.; Lier, J.B.V. Struvite crystallization under a marine/brackish aquaculture condition. Bioresour. Technol. 2016, 218, 1151–1156. [Google Scholar] [CrossRef]
  41. Moseke, C.; Saratsis, V.; Gbureck, U. Injectability and mechanical properties of magnesium phosphate cements. J. Mater. Sci. Mater. Med. 2011, 22, 2591–2598. [Google Scholar] [CrossRef]
  42. Feng, P.; Zhou, S.G.; Zhao, J. Effect of cattle manure ash on workability and mechanical properties of magnesium phosphate cement. Constr. Build. Mater. 2016, 129, 79–88. [Google Scholar] [CrossRef]
  43. Lahalle, H.; Céline, C.; Mesbah, A.; Lambertin, D.; Cannes, C.; Delpech, S.; Gauffinet, S. Investigation of magnesium phosphate cement hydration in diluted suspension and its retardation by boric acid. Cem. Concr. Res. 2016, 87, 77–86. [Google Scholar] [CrossRef]
Figure 1. Illustration of overall experimental design for this study.
Figure 1. Illustration of overall experimental design for this study.
Materials 15 06815 g001
Figure 2. Schematic of the design points with coded values.
Figure 2. Schematic of the design points with coded values.
Materials 15 06815 g002
Figure 3. Three-dimension response surface plots of initial mini slump (Y1) in relation to: (a) W/S ratio (X1) and M/P ratio (X2); (b) W/S ratio (X1) and borax dosage (X3); (c) M/P ratio (X2) and borax dosage (X3).
Figure 3. Three-dimension response surface plots of initial mini slump (Y1) in relation to: (a) W/S ratio (X1) and M/P ratio (X2); (b) W/S ratio (X1) and borax dosage (X3); (c) M/P ratio (X2) and borax dosage (X3).
Materials 15 06815 g003
Figure 4. Three-dimension response surface plots of mini-slump loss (Y2) in relation to: (a) W/S ratio (X1) and M/P ratio (X2); (b) W/S ratio (X1) and borax dosage (X3); (c) M/P ratio (X2) and borax dosage (X3).
Figure 4. Three-dimension response surface plots of mini-slump loss (Y2) in relation to: (a) W/S ratio (X1) and M/P ratio (X2); (b) W/S ratio (X1) and borax dosage (X3); (c) M/P ratio (X2) and borax dosage (X3).
Materials 15 06815 g004
Figure 5. Three-dimension response surface plots of yield stress (Y3) in relation to: (a) W/S ratio (X1) and M/P ratio (X2); (b) W/S ratio (X1) and borax dosage (X3); (c) M/P ratio (X2) and borax dosage (X3).
Figure 5. Three-dimension response surface plots of yield stress (Y3) in relation to: (a) W/S ratio (X1) and M/P ratio (X2); (b) W/S ratio (X1) and borax dosage (X3); (c) M/P ratio (X2) and borax dosage (X3).
Materials 15 06815 g005aMaterials 15 06815 g005b
Figure 6. Three-dimension response surface plots of plastic viscosity (Y4) in relation to: (a) W/S ratio (X1) and M/P ratio (X2); (b) W/S ratio (X1) and borax dosage (X3); (c) M/P ratio (X2) and borax dosage (X3).
Figure 6. Three-dimension response surface plots of plastic viscosity (Y4) in relation to: (a) W/S ratio (X1) and M/P ratio (X2); (b) W/S ratio (X1) and borax dosage (X3); (c) M/P ratio (X2) and borax dosage (X3).
Materials 15 06815 g006aMaterials 15 06815 g006b
Table 1. Actual values for the variables used in the experimental design.
Table 1. Actual values for the variables used in the experimental design.
Independent VariablesSymbolsActual Values for the Coded Values
−α (−1.682)−10+1+α (+1.682)
W/S ratioX10.180.200.240.280.30
M/P ratioX224.0379.9712
borax dosageX30.130.140.160.170.18
Table 2. The points for 3 factors CCD design (coded value).
Table 2. The points for 3 factors CCD design (coded value).
RunsW/S Ratio (X1)M/P Ratio (X2)Borax Dosage (X3)
1000
2−11−1
3001.682
4000
5000
6−1−1−1
7−1.68200
800−1.682
91.68200
1011−1
11000
121−11
13−111
14−1−11
150−1.6820
16000
17111
18000
191−1−1
2001.6820
Table 3. The points for 3 factors CCD design (coded value) and corresponding responses.
Table 3. The points for 3 factors CCD design (coded value) and corresponding responses.
RunsVariables in Coded ValuesResponses
W/S Ratio
(X1)
M/P Ratio
(X2)
Borax Dosage
(X3)
Initial Mini Slump/mm(Y1)Mini-Slump
Loss
(Y2)
Yield Stress/Pa
(Y3)
Plastic Viscosity/Pa·s
(Y4)
1000159332.030.77
2−11−1140102122.65−1.47
3001.682156221.180.59
4000152221.760.70
5000155251.390.74
6−1−1−1142123.571.42
7−1.68200134967.383.00
800−1.6821562210.960.81
91.6820016461.170.36
1011−11751377.180.37
11000168361.460.66
121−11155100.680.33
13−111144314.202.10
14−1−11143131.840.90
150−1.6820156230.980.31
16000161311.130.69
17111158201.850.46
18000147171.400.72
191−1−116080.900.31
2001.682014710918.500.88
Table 4. Repeatability of the responses at central points.
Table 4. Repeatability of the responses at central points.
TestInitial Mini Slump
/mm
Mini-Slump LossYield Stress
/Pa
Plastic Viscosity
/Pa·s
Mean (n = 6)157.0027.331.530.71
Standard Derivation7.357.230.320.04
Standard Error3.002.950.130.02
Coefficient of Variation (%)4.6826.4520.785.44
Table 5. Analysis of variance (ANOVA) of the response model of initial mini-slump.
Table 5. Analysis of variance (ANOVA) of the response model of initial mini-slump.
SourceResponses
Initial Mini-Slump (Y1)
Sum of SquaresDFMSF-Valuep-Value (Prob > F)
Model1248.513416.179.670.0007
X11227.1011227.1028.530.0001
X20.2510.250.010.9397
X321.16121.160.490.4931
Residual688.291643.02
Lack of fit418.291138.030.700.7090
Pure Error270.00554.00
Cor Total1936.8019
R20.6446
Adeq
Precision
10.8699
Table 6. Analysis of variance (ANOVA) of the response model of mini-slump loss.
Table 6. Analysis of variance (ANOVA) of the response model of mini-slump loss.
SourceResponses
Mini-Slump Loss (Y2)
Sum of SquaresDFMSF-Valuep-Value (Prob > F)
Model19,992.1463332.025.180.0063
X11321.9011321.902.060.1753
X211,230.80111,230.8017.460.0011
X32506.0712506.073.900.0700
X1X2120.131120.130.190.6727
X1X3253.131253.130.390.5413
X2X34560.1314560.137.090.0195
Residual8361.6113643.20
Lack of fit8100.2881012.5419.370.0023
Pure Error261.33552.27
Cor Total28,353.7519
R20.7051
Adeq Precision7.6423
Table 7. Analysis of variance (ANOVA) of the response model of yield stress.
Table 7. Analysis of variance (ANOVA) of the response model of yield stress.
SourceResponses
Yield Stress (Y3)
Sum of SquaresDFMSF-Valuep-Value (Prob > F)
Model9716.0561619.345.130.0066
X11277.6611277.664.050.0655
X21836.1711836.175.820.0314
X31480.1811480.184.690.0496
X1X21624.2211624.225.140.041
X1X31642.5111642.515.20.0401
X2X31855.3211855.325.880.0307
Residual4104.8413315.76
Lack of fit4104.338513.045086.84<0.0001
Pure Error0.550.1
Cor Total13,820.8819
R20.7030
Adeq
Precision
9.6694
Table 8. Analysis of variance (ANOVA) of the response model of plastic viscosity.
Table 8. Analysis of variance (ANOVA) of the response model of plastic viscosity.
SourceResponses
Plastic Viscosity (Y4)
Sum of SquaresDFMSF-Valuep-Value (Prob > F)
Model6.8461.142.30.0979
X12.5712.575.180.0404
X20.0210.020.040.8384
X30.5710.571.150.3030
X1X20.4410.440.890.3623
X1X31.0811.082.180.1636
X2X32.1612.164.370.0569
Residual6.44130.5
Lack of fit6.4380.8533.79<0.0001
Pure Error0.0150.002
Cor Total13.2819
R20.5151
Adeq
Precision
6.1537
Table 9. Coefficient table for the variables of all responses.
Table 9. Coefficient table for the variables of all responses.
Initial Mini SlumpMini-Slump LossYield StressPlastic Viscosity
Coefficientp-ValueCoefficientp-ValueCoefficientp-ValueCoefficientp-Value
Intercept153.60 38.75 9.61 0.73
X19.48<0.0001−9.840.1753−9.670.0655−0.430.0404
X20.140.939728.680.001111.600.0314−0.040.8384
X3−1.240.4931−13.550.0700−10.410.04960.200.3030
X1X2 3.880.6727−14.250.04100.240.3623
X1X3 −5.630.541314.330.0401−0.370.1636
X2X3 −23.880.0195−15.230.03070.520.0569
Note: X1: W/S Ratio; X2: M/P Ratio; X3: borax dosage; X1X2: interaction between W/S ratio and M/P ratio; X1X3: interaction between W/S ratio and borax dosage; X2X3: interaction between M/P ratio and borax dosage.
Table 10. Characteristics of numerical optimisation.
Table 10. Characteristics of numerical optimisation.
ParametersImportanceWeightGoalPredict Value
W/S Ratio (X1)31In range0.280
M/P Ratio (X2)31In range7.528
Borax Dosage (X3)31In range0.170
Yield Stress/Pa51Minimise0.680
Plastic Viscosity/Pa·s31In range0.263
Initial Mini Slump/mm51Maximise161.858
Mini-Slump Loss51Minimise11.282
Desirability0.867
Table 11. The compressive strength of the cube specimens.
Table 11. The compressive strength of the cube specimens.
Runs1 d Compressive Strength/MPa7 d Compressive Strength/MPa
123.734.6
233.834.4
319.139.0
421.832.2
523.333.2
622.433.8
727.543.1
829.232.3
913.715.8
109.210.6
1121.030.8
1218.522.8
1321.722.9
1421.531.9
15n.a.n.a.
1621.732.8
1711.111.6
1819.927.2
1910.214.1
207.28.2
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Yue, Y.; Ren, J.; Yang, K.; Wang, D.; Qian, J.; Bai, Y. Investigation and Optimisation of the Rheological Properties of Magnesium Potassium Phosphate Cement with Response Surface Methodology. Materials 2022, 15, 6815. https://doi.org/10.3390/ma15196815

AMA Style

Yue Y, Ren J, Yang K, Wang D, Qian J, Bai Y. Investigation and Optimisation of the Rheological Properties of Magnesium Potassium Phosphate Cement with Response Surface Methodology. Materials. 2022; 15(19):6815. https://doi.org/10.3390/ma15196815

Chicago/Turabian Style

Yue, Yanfei, Jun Ren, Kai Yang, Danqian Wang, Jueshi Qian, and Yun Bai. 2022. "Investigation and Optimisation of the Rheological Properties of Magnesium Potassium Phosphate Cement with Response Surface Methodology" Materials 15, no. 19: 6815. https://doi.org/10.3390/ma15196815

APA Style

Yue, Y., Ren, J., Yang, K., Wang, D., Qian, J., & Bai, Y. (2022). Investigation and Optimisation of the Rheological Properties of Magnesium Potassium Phosphate Cement with Response Surface Methodology. Materials, 15(19), 6815. https://doi.org/10.3390/ma15196815

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop