A Resistivity Plate Loading Device for Assessing the Factors Affecting the Stiffness of a Cement-Stabilized Subgrade

The extent of mixing in the stabilization process and the control of the cement content (C) and water content (w) in the mixture are key to the outcome of the engineering performance of a cement-stabilized subgrade. Intelligent Compaction (IC) quality control has improved quality control and management practices during construction. Intelligent Compaction Measurement Values (ICMVs) selected to evaluate the stiffness properties of cement-stabilized soils do not directly relate to the stiffness properties of the cement-stabilized subgrade and do not consider w and C. Additional tests need to be conducted for calibration of ICMVs. In this study, our solution is the development of a resistivity plate loading test. The resistivity plate loading test features the flexibility in determining the soil stiffness, w, C, and other important factors, such as the time of test effect (hydration) (T) and dry density (ρd). To verify the accuracy of the testing method, laboratory experimental studies were conducted on cemented soils considering ρd, w, C, and T at different factor levels. Multiple response studies based on grey rational analysis (GRA) were conducted. Analysis of the input factors was performed, and their effects on the measured responses were quantified. According to the study, the ρ measured by the device was a powerful indicator of stiffness, ρd, w, C, and T, which showed that the device can be useful equipment for quality control and an advancement in the in situ testing technologies and test equipment. A statistical regression model based on the linear and linear plus interaction terms among the factors is proposed to predict the average responses.


Introduction
Soil stabilization with cement is a well-established practice for constructing rail/roadway embankments and pavement layers. Cement stabilization allows the improvement of both standard soils and substandard in situ soils to levels consistent with the construction requirements. The stiffness and performance of stabilized subgrades are related to how close the stabilized subgrade is to achieving optimal compaction properties. As a result, standard test methods, fundamental analysis, design procedures, and quality evaluation are available to obtain acceptable results [1][2][3][4][5]. However, significant challenges still exist for the accurate quality assessment of compacted soil properties under in situ conditions. Current in situ quality evaluation testing techniques include stiffness test methods [6][7][8][9][10][11][12][13]. Stiffness spot test methods are less frequently employed, as the number of test points is limited, certain test methods are time-consuming, and they cause interference with the construction operations. Intelligent compaction (IC) quality control that integrates a vibratory roller with an accelerometer-based measuring system, a high-precision global positioning system (GPS), and an on-board data acquisition system to monitor the compaction process is highly recommended [8,10]. IC quality control overcomes the limitations of traditional stiffness spot test methods. Studies suggest that stiffness spot test methods Figure 1a shows a section of a typical resistivity plate loading device and Figure 1b shows a picture of the device, which comprises mild steel, a resistivity metre, and titanium electrodes (diameter of 4.0 mm and height of 7.0 mm) properly insulated with acetal plastic. The Wenner four-probe method of ρ measurement was employed. This technique is commonly applied in geology and soil science ρ measurements. Aside from the Wenner configuration method, several other electrode configuration methods are available in the literature and are rich in theory. Refs. [36,37] discuss the merits and demerits of the various electrode configuration methods. The suitability of a particular array depends on the signalto-noise ratio, depth of investigation, geometry of the electrode, etc. The Wenner four-probe electrode arrangement was chosen, and instrumental sensitivity was not as important as in the other array geometries. The electrode spacing determines the depth of electrical current penetration and volume measurement. In a homogeneous soil, the soil volume measured is approximately πa 3 , where (a) is the electrode spacing. In an experiment, a resistivity meter is connected to the loading plate, i.e., two electrodes (A and B) for injecting current into the stabilized subgrade and two voltage electrodes (M and N) for measuring the potential difference. Equation (1) is used to convert the current and voltage to an apparent resistivity.
where K is the geometrical factor dependent on the electrode configuration; for the Wenner configuration, K = 2πa, π is 3.14, a is the electrode space, V is the voltage, and I is the current.
four-probe electrode arrangement was chosen, and instrumental sensitivity was not as important as in the other array geometries. The electrode spacing determines the depth of electrical current penetration and volume measurement. In a homogeneous soil, the soil volume measured is approximately πa 3 , where (a) is the electrode spacing. In an experiment, a resistivity meter is connected to the loading plate, i.e., two electrodes (A and B) for injecting current into the stabilized subgrade and two voltage electrodes (M and N) for measuring the potential difference. Equation (1) is used to convert the current and voltage to an apparent resistivity.
V K I    (1) where K is the geometrical factor dependent on the electrode configuration; for the Wenner configuration, K = 2πa, π is 3.14, a is the electrode space, V is the voltage, and I is the current. The ρPLT method can be used to assess a whole range of soils, including subgrade, subbase, base, and treated soils used for road and railway construction. A calibration procedure may be required to avoid errors in ρ measurements when ρPLT is conducted on soils in rocks [38]. The electrodes of the resistivity plate loading device may also suffer contact resistance problems when installed on frost, very dry soils, and soils in rock [38,39].
In an experimental setup, ρ measurement errors have two origins: (1) poor electrode contact and (2) instrumental noise. The electrode depth (7.00 mm) and diameter (4.00 mm) induce minor faults in the former. In the latter, measurement time causes minimal or no variation in the ρ measurement. The properties of the electrode geometry and configuration are similar to other laboratory measuring instruments available in the literature [40].
Acetal plastic has excellent mechanical and low moisture absorption properties. Additionally, the coefficient of friction of acetal is within the range of metals that are commonly employed for static loading plates. As a result, the acetal plastic may have minimal or no effect on the stiffness measured with this device.

Experimental Program
The experimental program involves (1) testing the basic properties of the soil, (2) determining the most significant factors that affect the compaction quality of the cementstabilized subgrade and their factor levels, and (3) selecting an orthogonal array and running the experiments based on the orthogonal array. The experiments conducted included the unconfined compression strength (UCS) test and resistivity subgrade reaction modulus test (ρK30). The ρPLT method can be used to assess a whole range of soils, including subgrade, subbase, base, and treated soils used for road and railway construction. A calibration procedure may be required to avoid errors in ρ measurements when ρPLT is conducted on soils in rocks [38]. The electrodes of the resistivity plate loading device may also suffer contact resistance problems when installed on frost, very dry soils, and soils in rock [38,39].
In an experimental setup, ρ measurement errors have two origins: (1) poor electrode contact and (2) instrumental noise. The electrode depth (7.00 mm) and diameter (4.00 mm) induce minor faults in the former. In the latter, measurement time causes minimal or no variation in the ρ measurement. The properties of the electrode geometry and configuration are similar to other laboratory measuring instruments available in the literature [40].
Acetal plastic has excellent mechanical and low moisture absorption properties. Additionally, the coefficient of friction of acetal is within the range of metals that are commonly employed for static loading plates. As a result, the acetal plastic may have minimal or no effect on the stiffness measured with this device.

Experimental Program
The experimental program involves (1) testing the basic properties of the soil, (2) determining the most significant factors that affect the compaction quality of the cementstabilized subgrade and their factor levels, and (3) selecting an orthogonal array and running the experiments based on the orthogonal array. The experiments conducted included the unconfined compression strength (UCS) test and resistivity subgrade reaction modulus test (ρK 30 ).

Soil Properties
The soil used for this study was obtained from Hebei, China. The soil is classified as A-2-4 according to the American Association of State Highway and Transportation Official (ASSHTO) and SP-SM according to the Unified Soil Classification System (USCS). Figure 2 shows the particle size distribution curve. The liquid limit is 29.52% and the plastic limit is 19.19%. The uniformity coefficient (Cu) is 48.46l, the coefficient of gradation (Cc) is 0.57, and the group index (GI) is 0.
The soil used for this study was obtained from Hebei, China. The soil is classified as A-2-4 according to the American Association of State Highway and Transportation Official (ASSHTO) and SP-SM according to the Unified Soil Classification System (USCS). Figure 2 shows the particle size distribution curve. The liquid limit is 29.52% and the plastic limit is 19.19%. The uniformity coefficient (Cu) is 48.46l, the coefficient of gradation (Cc) is 0.57, and the group index (GI) is 0. ASTM D558-96 procedures were followed to determine compaction characteristics. The test results were applied to determine the maximum dry density (MDD) and optimum moisture content (OMC). The test was conducted by adding ordinary Portland cement of grade 32.5 to the soil in proportions of 5%, 8%, and 12% by weight of the total soil dry mass. The mixture was thoroughly mixed until uniformity was achieved. Figure 3 shows the compaction characteristics. Twelve percent C was considered to give the best results. Ref. [37] recommended adding 5-12% C to soils with similar properties.  ASTM D558-96 procedures were followed to determine compaction characteristics. The test results were applied to determine the maximum dry density (MDD) and optimum moisture content (OMC). The test was conducted by adding ordinary Portland cement of grade 32.5 to the soil in proportions of 5%, 8%, and 12% by weight of the total soil dry mass. The mixture was thoroughly mixed until uniformity was achieved. Figure 3 shows the compaction characteristics. Twelve percent C was considered to give the best results. Ref. [37] recommended adding 5-12% C to soils with similar properties.
A-2-4 according to the American Association of State Highway and Transportation Official (ASSHTO) and SP-SM according to the Unified Soil Classification System (USCS). Figure 2 shows the particle size distribution curve. The liquid limit is 29.52% and the plastic limit is 19.19%. The uniformity coefficient (Cu) is 48.46l, the coefficient of gradation (Cc) is 0.57, and the group index (GI) is 0. ASTM D558-96 procedures were followed to determine compaction characteristics. The test results were applied to determine the maximum dry density (MDD) and optimum moisture content (OMC). The test was conducted by adding ordinary Portland cement of grade 32.5 to the soil in proportions of 5%, 8%, and 12% by weight of the total soil dry mass. The mixture was thoroughly mixed until uniformity was achieved. Figure 3 shows the compaction characteristics. Twelve percent C was considered to give the best results. Ref. [37] recommended adding 5-12% C to soils with similar properties.

Development of Taguchi Orthogonal Array
The ρ d , C, w, T, and other factors, such as the particle size and shape of the soil, type of cement, porosity, temperature, etc., affect the in situ strength and compaction quality control of the cement-stabilized subgrade [17,18,41]. The Taguchi method may be employed to assess the essential factors for compaction quality control and quality assurance. The Taguchi method works on the principles of an orthogonal array and provides fewer variances in results, with an optimal set of control parameters. The principle allows the characterization of complex behavior, maximizing the test coverage while minimizing the test cases to be considered and the ease of data analysis [42,43]. Table 1 presents the factor levels, and Table 2 presents details of the Taguchi orthogonal array design used in this study. The control factors were carefully selected based on previous research [1,44,45]. ρ d , w, C, and T were selected as control input parameters, and their corresponding levels were determined as shown Table 1.  The Taguchi method works best in processes and procedures where one quality characteristic is to be evaluated. However, many practical applications are a compendium of quality characteristics. Multicriteria decision-making methods are often utilized in such complex scenarios. Recommended multicriteria decision-making methods include grey relational analysis (GRA), techniques for order preferences by similarity to ideal solution (TOPSIS), the analytic hierarchy process (AHP), etc. More than one quality characteristic was assessed, turning the multiple responses into a single response problem. Multicriteria decision-making methods have been suggested to be useful, easy, and efficient when integrated with the Taguchi method [46], with the most common methodologies being GRA and TOPSIS with several modifications. The GRA method was employed in this study to combine all the multiple response values into a single response value; thus, an interrelationship between the responses based on the grey relational coefficient was obtained. The GRA method has the following advantages: the results depend on the original values of the measured responses; the calculations are simple and suitable for multiple complicated relationships between responses and fairness comparison; and dimensional attributes are ensured by normalization. The steps to employ the GRA technique are as follows:

ID Taguchi Orthogonal Array Test Results
Step 1. Define a set A of n alternatives concerning W evaluation criteria as follows . . , n and j = 1, 2, . . . , m) where b ij represents the response value of the ith alternative on the jth criterion. Step 2. If the quality characteristic expectancy maximizes the response, use the larger, the better Equation (3). Use the smaller, the better Equation (4) when the quality characteristic minimizes the response. Use Equation (5) for nominal quality characteristics.
where m is the number response; ∆ min is minimum value of ∆; ∆ max is maximum value of ∆; λ is the distinguishing coefficient which is defined in the range 0 ≤ λ ≤ 1, (in this study 0.5 was used); GC ij is the gray relational coefficient for the ith experiment and jth response; T is the target value. In this study, average response was used.

Soil Preparation
Large soil lumps were broken with a wooden mallet, oven-dried at a temperature of 105 • C, and allowed to cool. All calculations were based on the oven-dried soil mass. The required amount of soil and w were calculated, thoroughly mixed in a bowl, transferred into plastic bags, and then tightly sealed. The samples were then kept for a day for moisture equilibration. Afterwards, the required amount of C was calculated, added, and mechanically mixed until the clumps formed were less than 5% of the total mix.
Samples prepared for the model test were cast into the modified soil consolidation soil chamber unit and subjected to an impact compaction test. The weight of the material in the container was calculated. The mixture was compacted in three equal layers. In each layer, a temperature sensor that was properly insulated was installed to monitor the temperature evolution during the curing process. The temperature evolution has an important consequence on the strength development and ρ measurement [47,48]. The specimens were cured in an open dry place under sealed conditions. All samples were cast into a thin, insulated steel cylinder, and their upper surfaces were tightly sealed using a plastic cover.
Two (2) samples were prepared for the UCS test, adopting the sample plan in Table 2. The ASTM D 1633 (2014) standard requires samples with a diameter, height, and height-todiameter ratio (h/d) of 101.60 mm, 116.40 mm, and 1.15, respectively, or samples with a diameter, height, and h/d ratio of 71.10, 142.20 mm, and 2.00, respectively. In this study, samples were prepared at a height of 80.00 mm and diameter of 39.10 mm with an h/d ratio of 2.05. The height and weight of the prepared samples were within an accuracy of ±1%. The prepared samples were stored in a zip-lock polythene bag under damp conditions. The mixing, compaction process, and each test was completed in less than two (2) hours. The UCS was tested in a triaxial loading frame with a computerized data-acquisition system at a loading rate of 1.00 mm/min (ASTM D1633). Figure 4c shows the details of the test setup, which comprises a loading system (air cylinder, air compressor, air pressure regulator, and air pressure gauge), deformation measurement system (linear variable displacement transducers (LVDTs)), computerized data-acquisition system, ρPLT device, temperature sensors (embedded within the compacted soil; bottom, middle, and upper-lower with rubber on it), and a soil-containing unit (diameter of 25.20 cm; height of 27.00 cm).

Test Setup
samples were prepared at a height of 80.00 mm and diameter of 39.10 mm with an h/d ratio of 2.05. The height and weight of the prepared samples were within an accuracy of ±1%. The prepared samples were stored in a zip-lock polythene bag under damp conditions. The mixing, compaction process, and each test was completed in less than two (2) hours. The UCS was tested in a triaxial loading frame with a computerized data-acquisition system at a loading rate of 1.00 mm/min (ASTM D1633). Figure 4c shows the details of the test setup, which comprises a loading system (air cylinder, air compressor, air pressure regulator, and air pressure gauge), deformation measurement system (linear variable displacement transducers (LVDTs)), computerized data-acquisition system, ρPLT device, temperature sensors (embedded within the compacted soil; bottom, middle, and upper-lower with rubber on it), and a soil-containing unit (diameter of 25.20 cm; height of 27.00 cm). The temperature sensors utilized a Pro Led DC 12 V digital thermometer and sensor probes with a temperature range of −50 to 110 °C.

Test Setup
The soil-containing unit was constructed of a steel frame with a height of 30.0 cm and a diameter of 25.2 cm. A plastic bag was inserted into the container, and the soil was packed into and compacted using the impact method. The plastic bag prevented the walls from conducting electrical current during testing. The walls of the soil container were also such that they exerted minimal or no effects on the stresses and strains during testing.
The soil-containing unit edges exert a border effect on the infinite value of the resistivity of air. These border effects are significant challenges for laboratory ρ measurements. Solutions have been developed in [49]. The ρ measured in our study is affected by this border effect. However, we do not intend to correct this measurement. The testing conditions are similar to those likely to be experienced for in situ conditions in a typical railway The temperature sensors utilized a Pro Led DC 12 V digital thermometer and sensor probes with a temperature range of −50 to 110 • C.
The soil-containing unit was constructed of a steel frame with a height of 30.0 cm and a diameter of 25.2 cm. A plastic bag was inserted into the container, and the soil was packed into and compacted using the impact method. The plastic bag prevented the walls from conducting electrical current during testing. The walls of the soil container were also such that they exerted minimal or no effects on the stresses and strains during testing.
The soil-containing unit edges exert a border effect on the infinite value of the resistivity of air. These border effects are significant challenges for laboratory ρ measurements. Solutions have been developed in [49]. The ρ measured in our study is affected by this border effect. However, we do not intend to correct this measurement. The testing conditions are similar to those likely to be experienced for in situ conditions in a typical railway or a road embankment where the geometry may be a 2-or 3-dimensional structure surrounded by air. Accounting for these possible conditions will be complicated.

Temperature
The temperature sensors were monitored every fifteen minutes for the first two hours. Afterwards, the temperature sensors were read every two hours during the daytime until the curing periods ended.

Resistivity Subgrade Reaction Modulus Test
ρPLT was selected to assess the properties of the model samples. The subgrade reaction modulus (K 30 ) test was conducted according to the procedures of TB10621-2014. The required load was applied through the piston to the load plate from an air cylinder supported by a reaction frame. The air cylinder was regulated through the air regulator. When the required air was supplied to the air cylinder, the valve that opened to the air cylinder was closed until the next load application was opened. This step was repeated throughout the test. LVDTs with an accuracy of 0.01% of the full range (100 mm) were mounted on three (3) sides of the load plate and configured to a computer that automatically displayed the deformation readings. The average deformation reading was employed for all calculations. For accuracy, all the devices were calibrated after three (3) sample runs during testing. The equivalent K 30 was calculated with Equation (8).
where K 30 is the equivalent subgrade reaction modulus, s D = 1.25 mm multiplied by (D/30 cm); D is the diameter of the resistivity loading plate, and σ is the stress at s D . ρ was measured with an ETCR3000B digital grounding resistance soil resistivity tester. ρ was recorded five minutes after every load application, which reduces the variability of ρ.

Unconfined Compression Strength Test
After curing the samples for 28 days, the samples were immersed in water for no less than four (4) h. Afterwards, the soil samples were removed from the water, dampened, and tested. During testing, loads were applied until the load at failure was recorded. The UCS strength was calculated using Equation (9). The ASTMD 1633 recommends multiplying the samples prepared with method B by 1.10 to obtain the calculated UCS strength. Based on laboratory studies, [50] suggested multiplying samples with an h/d ratio of 2.0 by a factor of 0.86. This step was purposely applied to convert the strength of an h/d ratio of 2.00 to that of the h/d ratio of 1.15, which is commonly utilized in routine soil-cement testing. In this study, a factor of 1.10 was applied.
where q u is the UCS strength; P is the force applied; ε is axial strain; A 0 is the cross-sectional area of the specimen; and F the strength correction factor.

Compaction
The compaction characteristics of soils are evaluated regarding OMC and MDD. Figure 3 shows the compaction characteristics. The sample with 0% cement content exhibited the highest MDD compared with the other samples. The relatively low MDD of the samples with cement additives, particularly 12%, can be attributed to (1) the cement additive, which causes aggregation of the soil particles to occupy large spaces that change the effective particle size of the soil; (2) the absorption of the compaction energy by the hydration products of the sample (the OMC increases with the addition of C due to the increase in w required for cement hydration); and (3) the loss of moisture due to evaporation from heat generated during hydration [13,51]. Figure 5 shows the temperature-time plot history for the tested samples. The sample with the highest temperature exhibited a peak temperature of 13 • C (1) and a low temperature of 1.0 • C (3). An exothermic chemical reaction between cement and water occurs when they are blended. Temperature, relative humidity, type of cement, and C can affect the reaction. The relatively low temperature observed was caused by environmental factors [52]. High-temperature conditions enhance cement hydration and pozzolanic reactions in cement-stabilized soils, thus improving the strength properties. However, very high curing temperatures (20-60 • C) have caused cement hydration products to be arranged randomly, thus producing large pores in soils (crossover effect) [48]. Elsewhere, lower temperature conditions have also been reported to contribute to poor strength development [53]. The different temperature gradients affect the strength development and ρ measurement [48]. Figure 5 shows the temperature-time plot history for the tested samples. The sample with the highest temperature exhibited a peak temperature of 13 °C (1) and a low temperature of 1.0 °C (3). An exothermic chemical reaction between cement and water occurs when they are blended. Temperature, relative humidity, type of cement, and C can affect the reaction. The relatively low temperature observed was caused by environmental factors [52]. High-temperature conditions enhance cement hydration and pozzolanic reactions in cement-stabilized soils, thus improving the strength properties. However, very high curing temperatures (20-60 °C) have caused cement hydration products to be arranged randomly, thus producing large pores in soils (crossover effect) [48]. Elsewhere, lower temperature conditions have also been reported to contribute to poor strength development [53]. The different temperature gradients affect the strength development and ρ measurement [48].

Resistivity Subgrade Reaction Modulus
The study clearly showed the dependence of stiffness and ρ on ρd, w, C, and T. All ρ readings were temperature corrected using Equation (10). Figure 6a presents a typical stress-and temperature-corrected ρ, and the change in ρ with successive load application was a kind of multistage process during testing. We envisaged that the contacts between the electrodes and the soil would improve with successive loads acting on the loading plate. As forces were applied to the loading plate, comparable ρ changes were observed. The changes in ρ were analogous to the changes in the soil properties, interaction of the w, cement hydration products, and influence of the applied stress beneath the loading plate. For a given T, the higher the cement content is, the greater the number of hydration compounds, thus the higher the observed ρ [32]. Additionally, the lower w is, the higher the observed ρ. The selected ρ was calculated by establishing a linear equation between the measured ρ and the deformation. ρ is determined at the equivalent deformation used to calculate the subgrade reaction modulus. The test results are summarized in Table 2.

Resistivity Subgrade Reaction Modulus
The study clearly showed the dependence of stiffness and ρ on ρ d , w, C, and T. All ρ readings were temperature corrected using Equation (10). Figure 6a presents a typical stress-and temperature-corrected ρ, and the change in ρ with successive load application was a kind of multistage process during testing. We envisaged that the contacts between the electrodes and the soil would improve with successive loads acting on the loading plate. As forces were applied to the loading plate, comparable ρ changes were observed. The changes in ρ were analogous to the changes in the soil properties, interaction of the w, cement hydration products, and influence of the applied stress beneath the loading plate. For a given T, the higher the cement content is, the greater the number of hydration compounds, thus the higher the observed ρ [32]. Additionally, the lower w is, the higher the observed ρ. The selected ρ was calculated by establishing a linear equation between the measured ρ and the deformation. ρ is determined at the equivalent deformation used to calculate the subgrade reaction modulus. The test results are summarized in Table 2.
where α ≈ 0.025 ( • C −1 ), ρ w (t) is the temperature-compensated electrical resistance, t 0 is a fixed reference temperature, and ρ(t 0 ) is the ρ at a temperature t 0 . Figure 6a,b show the stress and resistivity curves and stress-strain curves for the tested samples. The samples exhibited different forms of stress and strain characteristics. Generally, strain hardening was observed for samples compacted at a higher w and low C. A similar observation was observed for certain samples with a low w. A linear elastic behavior was also observed for certain samples compacted at a higher ρ d .
where 0.025 is the temperature-compensated electrical resistance, 0 t is a fixed reference temperature, and   0 t  is the ρ at a temperature 0 t .
(a) (b) Figure 6. Resistivity plate loading test results (a) typical of stress-and temperature-corrected apparent resistivity curves from the concurrent resistivity subgrade reaction modulus test and (b) typical of the stress strain curves. Figure 6a,b show the stress and resistivity curves and stress-strain curves for the tested samples. The samples exhibited different forms of stress and strain characteristics. Generally, strain hardening was observed for samples compacted at a higher w and low C. A similar observation was observed for certain samples with a low w. A linear elastic behavior was also observed for certain samples compacted at a higher ρd.

Unconfined Compression Strength Test
The UCS is an index for quantifying the effectiveness of an additive [2,54] on the soil strength. The sample must not collapse, or no significant loss in strength should be observed during water immersion. The test results for UCS are presented in Table 2. The results range from 0.85 to 8.18. The highest results were associated with samples compacted at the highest ρd. The w and C content variation had a different effect on the UCS. UCS values in the range of 0.2-0.4 MPa are usually permitted for subgrade applications. However, considering the different compositions of the samples, they may have negative consequences on the performance of the compacted subgrade.

Range Analysis
Range analysis was performed in the Minitab statistical package to assess the influence of the various factors on the response variables and their ranks, as shown in Table 3. The rank based on the Delta statistics compared the relative magnitudes of their effects. The delta statistic is the highest average minus the lowest average for each factor. The ranks and average level responses are used to determine which element provides the best

Unconfined Compression Strength Test
The UCS is an index for quantifying the effectiveness of an additive [2,54] on the soil strength. The sample must not collapse, or no significant loss in strength should be observed during water immersion. The test results for UCS are presented in Table 2

Range Analysis
Range analysis was performed in the Minitab statistical package to assess the influence of the various factors on the response variables and their ranks, as shown in Table 3. The rank based on the Delta statistics compared the relative magnitudes of their effects. The delta statistic is the highest average minus the lowest average for each factor. The ranks and average level responses are used to determine which element provides the best results. All the means of the mean plots were determined using the nominal characteristics. Furthermore, all ANOVAs were performed at a confidence level of 95% and 5% significance level.

Subgrade Reaction Modulus
In Table 3, ρ d , T, w, and C are ranked in order of importance. The analysis of variance for the individual factor contributions and their F values of 1.74, 1.16, 0.67, and 0.62 are presented in the order of the ranks. As shown in Figure 7, a combination of 23 ρ d , 5% w, 18% C, and T of forty-eight (48) h is considered to give the best result. According to the literature, the strength of the cement-stabilized subgrade is related to the solid phase, w, C, and T [19,55]. An increase in ρ d yielded an increase in K 30 , as shown in Figure 7. An increased ρ d caused interlocking between the soil particles and cement particles, thus increasing the ability of the soil to resist the mechanical forces that act on it [56,57]. Additionally, an increase in w weakens the cohesion strength between the particles of the soil and the hydration products [18].   Table 3 shows the ranks of the input factors. The analysis of variance on the individual factor contributions and their F values according to the ranks (C, T, w, and ρd) are 3.19, 0.94, 0.46, and 0.38, respectively. C exhibits a significant effect on the ρ response. As shown in Figure 8, the line slope shows that the effect and the degree of the slope are comparative to the magnitude of the inputs. ρ increases with an increase in ρd and T and decreases with an increase in w. Several others have reported the same findings [32]. The C produces an amount of hydration products for a given T, resulting in a denser structure. The higher T  Table 3 shows the ranks of the input factors. The analysis of variance on the individual factor contributions and their F values according to the ranks (C, T, w, and ρ d ) are 3.19, 0.94, 0.46, and 0.38, respectively. C exhibits a significant effect on the ρ response. As shown in Figure 8, the line slope shows that the effect and the degree of the slope are comparative to the magnitude of the inputs. ρ increases with an increase in ρ d and T and decreases with an increase in w. Several others have reported the same findings [32]. The C produces an amount of hydration products for a given T, resulting in a denser structure. The higher T is, the higher the hydration products produced. Thus, the free water space and porosity decrease and tortuosity increases, increasing ρ. Additionally, an increase in ρ d causes interlocking and parking of the soil particles, which renders the electrical flow path tortuous, increasing ρ.

Unconfined Compressive Strength
Since all the samples were assessed at 28 days, the effect of curing was excluded. Subsequently, ρd, w, and C were evaluated to determine their impact on the UCS strength development. Table 3 shows the ranks of the input factors. The F values of the input factors according to the ranks are 5.93, 0.54, and 0.23. Figure 9 presents the main effect plots. The best strength was achieved at 1.80 ρd, 18% w, and 12% C. Several researchers using the UCS test have investigated the strength of cement-stabilized soil [44,[58][59][60]. Almost all researchers reported an increase in strength with an increase in ρd, and our results confirm this finding. The UCS strength increases with an increase in ρd and w up to 12% C, as shown in Figure 9. At 18% C, the UCS decreases, which is contrary to the results presented in [44,53,[58][59][60][61][62]. This observation can be explained by the observation that at higher C, the cement-treated soil requires enough w to undergo complete hydration [61,63]. Furthermore, at higher C, microcracks develop, easily propagating and joining under applied loads and thus reducing the UCS. This finding means that there is an optimum amount of w and C for which the number of bonds and their configuration offer the best strength development. Any deviation from these values will negatively impact the strength development.

Unconfined Compressive Strength
Since all the samples were assessed at 28 days, the effect of curing was excluded. Subsequently, ρ d , w, and C were evaluated to determine their impact on the UCS strength development. Table 3 shows the ranks of the input factors. The F values of the input factors according to the ranks are 5.93, 0.54, and 0.23. Figure 9 presents the main effect plots. The best strength was achieved at 1.80 ρ d , 18% w, and 12% C. Several researchers using the UCS test have investigated the strength of cement-stabilized soil [44,[58][59][60]. Almost all researchers reported an increase in strength with an increase in ρ d , and our results confirm this finding. The UCS strength increases with an increase in ρ d and w up to 12% C, as shown in Figure 9. At 18% C, the UCS decreases, which is contrary to the results presented in [44,53,[58][59][60][61][62]. This observation can be explained by the observation that at higher C, the cement-treated soil requires enough w to undergo complete hydration [61,63]. Furthermore, at higher C, microcracks develop, easily propagating and joining under applied loads and thus reducing the UCS. This finding means that there is an optimum amount of w and C for which the number of bonds and their configuration offer the best strength development. Any deviation from these values will negatively impact the strength development.  Figure 9. Effects of the control factors on USC.

Grey Rational Analysis
Once calculations for GC and G were finished, their values were entered in Table  ANOVA was performed to assess the effects of the experimental input parameters on t GRA grade. The results of the ANOVA are presented in Table 5. ρd, T, w, and C influenc the GRA grade values with contributions of 79.09%, 10.30%, 3.94%, and 6.67%, respe tively. A GRA near 1 in Table 4 can be considered a criterion that gives the best optim conditions for K30, ρ, and UCS.

Regression Analysis
In this study, regression analyses were employed for modelling and predicting t response variables. Two different models were initially proposed: linear interactions an linear plus interactions. The best fitting models were selected. The highest performin predictive equations obtained from the analysis are given below. The equation reveals t

Grey Rational Analysis
Once calculations for GC and G were finished, their values were entered in Table 4. ANOVA was performed to assess the effects of the experimental input parameters on the GRA grade. The results of the ANOVA are presented in Table 5. ρ d , T, w, and C influenced the GRA grade values with contributions of 79.09%, 10.30%, 3.94%, and 6.67%, respectively. A GRA near 1 in Table 4 can be considered a criterion that gives the best optimal conditions for K 30 , ρ, and UCS.

Regression Analysis
In this study, regression analyses were employed for modelling and predicting the response variables. Two different models were initially proposed: linear interactions and linear plus interactions. The best fitting models were selected. The highest performing predictive equations obtained from the analysis are given below. The equation reveals the nature of the relationships of the input parameters and their effect on the responses. Directly proportional responses are acknowledged using a + sign, while the -sign represents the opposite. Equations (11)- (14) were identified as the best-performing equations for K 30 , ρ, UCS, and GRA, respectively. The equations predicted the average K 30 , ρ, UCS, and GRA with adjusted coefficients of correlation (R 2 ) of 72.18%, 78.09%, 71.06%, and 55.63%, respectively. The individual unique predictor contributions for the equations are shown in Table 6.

Implications for Practice
Several studies have suggested that K 30 can reflect the strength, stiffness, and performance of compacted subgrade in China. However, we suggest that K 30 is the main parameter that is commonly employed but is not the only parameter that can assure com-pacted subgrade stiffness, performance, and deformation characteristics. Compacted soil typically remains in the elastic stage, away from the failure stage, as evidenced by the plate loading tests (refer to Figure 6b). As a result, using this indicator alone to capture the strength and stiffness characteristics is deemed ineffective. Additionally, the K 30 test results are affected by factors such as w, ρ d , T, and C (refer to Figure 7) for the cement-stabilized subgrade. The variability of w, ρ d , and C is reported to cause subgrade problems, such as swelling, subgrade cavity, and differential settlement.
There is no good indicator to reflect the changes in the physical and stiffness properties of the subgrade, particularly the K 30 , w, T, C, and ρ d of the cement-stabilized subgrade for construction. However, practitioners are aware of how w, T, C, and ρ d affect the test outcome of K 30 and the performance of the compacted subgrade. Nonetheless, a qualitative method to instantaneously and concurrently assess the w, T, C, and ρ d effects with K 30 has always been associated with certain problems. Our test method solves this problem and serves as an advancement in compaction measurement systems and in situ testing technologies.
The proposed regression Equations (11)-(13) are presented to predict the average responses of K 30 , ρ, and UCS, which showed high performance with R 2 values of 72.18%, 78.09%, and 71.06%, respectively. Use of the proposed equations is encouraged during the construction of cement-stabilized subgrade for similar soils utilized in this study with similar w and C requirements. Equation (14) was introduced as a unique equation that can give a combined unique identity to K 30 , UCS, and ρ, considering the input factors employed in the study with an R 2 of 55.63%. The test results presented in this study conclude that the proposed remediation was significant. Although our test method may be slow, its benefits outweigh those of commonly known methods in the literature.

Conclusions
A resistivity plate load device was developed to evaluate the properties of the cementstabilized subgrade. Experimental studies were conducted on compacted cement soils according to the Taguchi orthogonal L9 array. Multivariate stepwise regression models were developed to predict K 30 , ρ, and UCS considering ρ d , T, w, and C for compaction quality evaluation purposes. Our conclusions are as follows.
(1) The resistivity plate load device can be efficiently used to assess the w, C, ρ d , and T effects on K 30 for the cement-stabilized subgrade for construction. Therefore, the device is recommended for field use. This recommendation must be evaluated and validated in the field with a vast variability in compacted soils with variable w, C, and ρ d properties that may confirm the implications of our laboratory observation. (2) K 30 is significantly affected by C, w, ρ d , and T. It is important within subgrade construction control protocols to confirm that compacted cement-stabilized subgrade C, ρ d , T, and w are considered for a well-informed decision on K 30. (3) Different values of C, ρ d , and w have different effects on the USC, leading to variable USC strength values and confirming the need to control C, ρ d , and w during construction for enhanced performance. (4) The proposed regression models showed better performance, and therefore, are recommended to predict the average K 30 , ρ, UCS, and GRA grades during construction of the cement-stabilized subgrade for similar soils.
The statistical analysis presented in this study produced excellent results. The analysis should be treated as a simplified approach and a general guide that can be adopted for field compaction quality control for the cement-stabilized subgrade. Future studies include a field application of the test method and calibration of ICMVs for compaction management of the cement-stabilized subgrade. Notwithstanding, this study was a significant step towards providing a reliable reference for engineering practice in compaction quality control and a considerable step towards advancement in compaction measurement systems and in situ testing technologies.