# Data-Driven Prediction of Cement-Stabilized Soils Tensile Properties

Mario Castaneda-Lopez, Thomas Lenoir, Jean-Pierre Sanfratello, Luc Thorel

► To cite this version:

Mario Castaneda-Lopez, Thomas Lenoir, Jean-Pierre Sanfratello, Luc Thorel. Data-Driven Prediction of Cement-Stabilized Soils Tensile Properties. Infrastructures, 2023, 8 (10), pp.146. 10.3390/infrastructures8100146. [https://doi.org/10.3390/infrastructures8100146](https://doi.org/10.3390/infrastructures8100146)

HAL Id: hal-04262321

https://hal.science/hal-04262321v1

Submitted on 27 Oct 2023

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons CC BY 4.0 - Attribution - International License
<figure><img src="image_2.png" /></figure>

Article

# Data-Driven Prediction of Cement-Stabilized Soils Tensile Properties

Mario Castaneda-Lopez <fnref n="1" />,<fnref n="2" /> ID, Thomas Lenoir <fnref n="1" />,* ID, Jean-Pierre Sanfratello <fnref n="3" /> and Luc Thorel <fnref n="2" /> ID

**Abstract:** The indirect tensile strength of two geomaterials treated with variable cement contents, degrees of compaction and water contents were tested after several curing times. A statistical review through an analysis of variance allows for identifying the significant variables and generating prediction models. The distribution of associated uncertainties was measured. Based on these probabilistic results, numerical models were constructed using Latin Hypercube Sampling as the space filling technique. Predictions from the numerical sampling were in accordance with the experimental results. The numerical results suggest that the net gain in accuracy was not affected by the soil type. In addition, it increases rapidly as a function of the sampling size. The proposed approach is broad. It can help to highlight the physical mechanisms involved in behaviors of multi-component materials.

**Keywords:** cement-stabilized soils; indirect tensile strength; analysis of variance; response surface method; DOE; numerical experiments

<figure><img src="image_3.png" /></figure>

check for
updates

**Citation:** Castaneda-Lopez, M.; Lenoir, T.; Sanfratello, J.-P.; Thorel, L. Data-Driven Prediction of Cement-Stabilized Soils Tensile Properties. *Infrastructures* **2023**, *8*, 146. https://doi.org/10.3390/infrastructures8100146

Academic Editor: Francesca Dezi

Received: 22 August 2023

Revised: 21 September 2023

Accepted: 6 October 2023

Published: 12 October 2023

<figure><img src="image_4.png" /></figure>

**Copyright:** © 2023 by the authors. Licensee MDPI, Basel, Switzerland. This article is an open access article distributed under the terms and conditions of the Creative Commons Attribution (CC BY) license (https://creativecommons.org/licenses/by/4.0/).

## 1. Introduction

Adding small amounts of cement as a binder is a common technique to enhance the engineering and mechanical properties of in situ soils according to stress rates imposed by civil engineering projects [1,2]. This practice has the advantage of minimizing environmental impacts and reducing infrastructure cost [3], especially in terms of preserving quality natural resources. The use of in situ soil avoids the use of quarry materials [4].

Cement stabilization consists of incorporating a volume of cement-based additive to the soil. In the presence of water, pozzolanic reactions occur over short and long terms, and improve the stiffness and strength of the mix [5]. The current state of practice tends to limit the application of cement to coarse and low plastic soils [6,7]; however, satisfactory performances have been reported with fine-grained and even with expansive soils [8–11].

These treated materials are called Cement-Stabilized Soils (CSS). They constitute an intermediate class of geomaterials whose mechanical behavior is at the boundary between classical soil mechanics and rock mechanics [12,13]. Cement stabilization is widely used in civil engineering applications such as embankments, slope protection of dams, pavement of highways, building pads and foundation stabilization [5,14]. In the field of transportation geotechnics, CSSs are generally used in capping layers and can be used in a road base or subbase [13].

The amounts of hydraulic binder in CSSs commonly range between 5 and 10% per weight of dry material. However, cement manufacturing is one of the most significant sources of anthropogenic greenhouse gases, in particular with carbon dioxide (CO₂) emissions, and production requires the use of considerable energy and materials [15]. Total emissions from the cement industry contributes as much as 8% of global CO₂ emissions [16].

<footnote n="1">¹ Advanced Materials for Transportation Infrastructures Laboratory, Materials and Structures Department, University Gustave Eiffel, F-44344 Bouguenais, France; mario-alexander.castaneda-lopez@univ-eiffel.fr</footnote>

<footnote n="2">² Centrifugues for Geotechnics Laboratory, Geotechnical Engineering, Environment, Natural Hazards and Earth Sciences Department, University Gustave Eiffel, F-44344 Bouguenais, France; luc.thorel@univ-eiffel.fr</footnote>

<footnote n="3">³ DTRD, COLAS SA, 78771 Magny-les-Hameaux, France; jeanpierre.sanfratello@colas.com</footnote>

<footnote>* Correspondence: thomas.lenoir@univ-eiffel.fr</footnote>
Thus, in the field of road geotechnics, due to the amount of material involved, the optimization of dosage procedures is a key issue. Decreasing cement amounts could bring a real environmental benefit.

Treated layers of land transport infrastructures are subjected to tensile stress at their bottom. This critical mechanical characteristic is generally used for material evaluation and structure design [2,17,18]. For an easy-performing, time-saving and cost-effective method, this feature is mainly ascertained from the measurement of Indirect Tensile Strength (ITS) at the laboratory scale [19].

The mechanical properties of CSSs depend on variables such as cement content, cement type, moisture content, curing time, degree of compaction, temperature during curing, aggregate type, volume of voids and porosity/binder ratio [2,14,17–27].

To consider the possible performance variations due to construction procedures in worksites, a laboratory test program is established during the design process by varying the mixture parameters. In France, for instance, mixture proportioning design follows a Central Composite Design (CCD) based on variations around the standard proctor values of the reference mixture, i.e., Optimum Moisture Content (OMC) and Maximum Dry Density (MDD). Performance is evaluated for a fixed curing time (long term, e.g., 90 days) [28].

This sensitivity analysis only focuses on unfavorable conditions. That limits the scope and these aspects are critical for rationalizing structure design, particularly when binder amounts are less than 5%. Unfortunately, increasing the number of variables can lead to an infeasible experimental program [29] and reducing it is unacceptable and may bias the conclusions [26]. Currently, there is a need to develop models based on a reasonable number of tests that allow a straightforward assessment of tensile properties.

Design Of Experiment (DOE) techniques have become crucial for dosage evaluations and the mixing optimization of CSSs [30,31]. DOE techniques, for example, can be implemented in numerical works to assess the influence of sampling on the accuracy of prediction models [29]. Similarly, Response Surface Methodology (RSM) as part of the DOE process has become a useful tool to rationalize cemented amounts in hydraulically bound materials. RSM is often used with innovative complex materials when interactions between constituents are not well understood [19,23,26,32–36].

The first aim of our work was to investigate, at the laboratory scale, the mechanical performances of CSSs treated with a low amount of binder considering the variability of field conditions. To achieve this goal, we measured the ITS of two geo-materials treated with 3% of cement. To consider the field conditions, the dosage preparations (i.e., curing time, compaction degree, cement content and water content) were varied around a reference mixture, following a CCD design. Then, an Analysis of Variance (ANOVA) and multilinear regression were applied to the experimental results to assess and predict the mechanical performance. Response surfaces and contour plots were generated to highlight the links between preparation parameters. The second aim of the study was to propose a tool to facilitate and support the selection of a given sampling in variability studies. Based on the ascertained regression models and measured uncertainties, numerical experiments were simulated using DOE. For a given number of numerical experiments, varying from 10 up to 100, 1000 samplings were carried out. Synthetic regression models were obtained and the distances between the experimental and numerical models were measured. Through this approach, we expressed the mean differences, and thus the accuracy as a function of the sampling size.

## 2. Materials and Methods

### 2.1. Materials

#### 2.1.1. Soils

Two soils named A and B were studied. The results of the characterization are shown in Table 1 and the grain size curves are shown in Figure 1.
* Soil A

This material was cut from the covering soil of an aggregate quarry at Soignies Belgium (50°34′25.759″ N, 4°2′16.476″ E). The maximum diameter of the larger elements in this material was less than 1 mm. Through an 80 µm sieve, the cumulated undersize was 94%, which means that its mechanical behavior was controlled by its fine fraction. It is classified as a fine-grained soil [37]. The clay activity is low, A = PI/C₂µm = 0.29. Swelling potential was thus not significant.

* Soil B

This material was a quarry tailing from Chauvé, France (47°8′22.886″ N, 2°0′45.934″ W). The maximum diameter of the larger elements in this material was less than 50 mm. Through an 8 mm sieve, the cumulated undersize was 93%, and fine content, less than 80 µm, was 17%, which means that its mechanical behavior was controlled by its coarse fraction. It is classified as a coarse-grained soil [37]. Despite a high clay activity, A = 1.3, the low clay content and low PI imply a low swelling potential.

<table><caption>**Table 1.** Physical properties of soils.</caption><thead><tr><th>Properties</th><th>Reference</th><th>Soil A</th><th>Soil B</th></tr></thead><tbody><tr><td>Density of soil particles, ρ<sub>s</sub></td><td>[38]</td><td>2.59 g/cm<sup>3</sup></td><td>2.64 g/cm<sup>3</sup></td></tr><tr><td>Liquid limit, W<sub>L</sub></td><td></td><td>30%</td><td>38%</td></tr><tr><td>Plastic limit, W<sub>P</sub></td><td>[39]</td><td>17%</td><td>25%</td></tr><tr><td>Plasticity index, (PI = W<sub>L</sub> - W<sub>P</sub>)</td><td></td><td>13%</td><td>13%</td></tr><tr><td>Clay content (d < 2 µm)<sup>a</sup>, C<sub>2µm</sub></td><td></td><td>24%</td><td>6%</td></tr><tr><td>Silt (2 µm < d < 63 µm)<sup>a</sup></td><td></td><td>69%</td><td>11%</td></tr><tr><td>Sand (63 µm < d < 2 mm)<sup>a</sup></td><td>[37]</td><td>6%</td><td>40%</td></tr><tr><td>Gravel (2 mm < d < 200 mm)<sup>a</sup></td><td></td><td>1%</td><td>43%</td></tr><tr><td>AASHTO<sup>b</sup> classification</td><td>[40]</td><td>A-6(11)</td><td>A-2-6(0)</td></tr><tr><td>GTR<sup>c</sup> classification d</td><td>[41]</td><td>B5</td><td>A2</td></tr><tr><td>USCS<sup>d</sup> classification c</td><td>[42]</td><td>CL</td><td>SM-SC</td></tr></tbody></table>

### 2.1.2. CSS—Treatment and Sample Preparation

Both soils were treated with 3% of cement (weight of dry material) CEM II/A-LL 42.5 R [44]. Lime pre-treatment was not needed due to the suitable workability in both mixtures. The OMC and the MDD relationships for both CSSs were obtained following Proctor compaction tests [45], for standard compaction effort (600 kN m/m³):

* OMC = 15.5%, MDD = 1.79 g/cm³ for mixture A

* OMC = 10.0%, MDD = 1.98 g/cm³ for mixture B

These values were taken as a reference for the sample preparation. Based on the French guidelines for soil treatment, a variability study following a fractional CDD was proposed [28]. Samples were prepared at several Varying conditions (V) around the reference mixture (V3) with four varying factors: Curing Time (CT), Cement Content (CC), Degree of Compaction (DC) and Water content (W). Table 2 shows all the parameters assessed, their definition and the levels of variation used as the reference in this work. The varied factor for each V is highlighted in bold. In total, seven treatments were proposed and for each one, at least three specimens were tested. The complete variability study included a total of 107 treated specimens for A and 106 for mixture B. Additionally, four Non-Treated (NT) samples were tested for each material.

<footnote><sup>a</sup> Soil size range based on European standard NF EN ISO 14688-1 [43]; <sup>b</sup> American Association of State Highway and Transportation Officials; <sup>c</sup> French Road Earthworks Manual; <sup>d</sup> Unified Soils Classification System.</footnote>
<figure><img src="image_5.png" /><figcaption>**Figure 1.** Particle size distribution of soils A and B. Soil size range: F (fine), M (medium) and C (coarse) [37,43].</figcaption></figure>

<table><caption>**Table 2.** Reference values for sample preparation for the variability study.</caption><thead><tr><th rowspan="2">V</th><th colspan="4">Varying Factors</th></tr><tr><th>Curing Time, CT (Days)</th><th>Cement Content, CC (%) <sup>a</sup></th><th>Degree of Compaction, DC (%) <sup>b</sup></th><th>Water Content, W (%) <sup>c</sup></th></tr></thead><tbody><tr><td>V1</td><td>7, 28, 90, 180, 360</td><td>3.0</td><td>94.0</td><td>100</td></tr><tr><td>V2</td><td>7, 28, 90, 180, 360</td><td>2.0</td><td>96.0</td><td>100</td></tr><tr><td>V3 (ref.)</td><td>7, 28, 90, 180, 360</td><td>3.0</td><td>96.0</td><td>100</td></tr><tr><td>V4</td><td>7, 28, 90, 180, 360</td><td>4.0</td><td>96.0</td><td>100</td></tr><tr><td>V5</td><td>7, 28, 90, 180</td><td>3.0</td><td>96.0</td><td>90</td></tr><tr><td>V6</td><td>7, 28, 90, 180</td><td>3.0</td><td>96.0</td><td>110</td></tr><tr><td>V7</td><td>7, 28, 90, 180, 360</td><td>3.0</td><td>98.5</td><td>100</td></tr><tr><td>NT</td><td>-</td><td>0.0</td><td>96.0</td><td>100</td></tr></tbody></table>

Samples were compacted in a double static mode in one layer into a 100 mm diameter (D), 100 mm high (H) cylindrical mold and they were stored in a room with temperature regulated at 20 °C and protected with plastic films to keep constant water content. No saturation process was made before the experiments. Details of the sampling preparation are available online [46].

## 2.2. Indirect Tensile Strength

All tests were carried out on a Zwick/Roell Z150 press with a crosshead speed of 0.3 mm/min in absolute value. Applied load was measured with a force sensor TME F521TC (max. load = 10 kN). Specimens were placed between two diametrically opposite generators and then subjected to a compressive load until failure [47]. ITS (MPa) was then calculated using Equation (1),

$$ ITS = \frac{2F_r}{\pi HD} \qquad (1) $$

<footnote><sup>a</sup> CC as percentage of dry material weight, <sup>b</sup> DC as percentage of MDD, <sup>c</sup> W as percentage of OMC.</footnote>
where F<sub>r</sub> is the applied failure load (N); H and D are the thickness and diameter of the specimen (mm).

2.3. Statistical Approach

2.3.1. Analysis of Variance

Analysis of Variance, ANOVA, is a statistical tool used to identify and establish the effect of variations of a given factor on the output variable, by grouping and comparing the mean values according to the levels of the factor under study. This analysis was considered due to the large data set obtained from the experimental program, and the number of investigated variables. This tool was applied to the experimental results to determine which factors significantly affected the ITS of both CSSs. The results of NT were not included.

An ANOVA tests the hypothesis that all group means are equal (i.e., the null hy-
pothesis, H₀) against the alternative supposition that at least one group mean is different
from the others (H₁). Null hypotheses are evaluated by a Fisher's test (F-test) to compare
the variation within groups with a significance level α = 0.05 as the threshold. If the
calculated probability (p-value) is lower than α, then H₀ is rejected as it inherently implies a
significant effect and thus the factor is considered in the prediction model. The hypothesis
of no differences in the treatment means was performed, as described by the means model
(Equation (2)) [48]. It was analyzed using type II ANOVA, i.e., considering only main ef-
fects. This approach is justified because mixture parameters are independent variables [49]

$$
y_{ijklm} = \bar{\mu} + \alpha_i + \beta_j + \gamma_k + \theta_l + \epsilon_{ijklm} \quad (2)
\left\{
\begin{array}{l}
i = 1, 2, 3, 4, 5 \\
j = 1, 2, 3 \\
k = 1, 2, 3 \\
l = 1, 2, 3 \\
m = 1, 2, 3
\end{array}
\right.
$$

where:

* y<sub>ijklm</sub> is the experimental measurement of ITS when the factors CT, CC, DC and W are in the i-th, j-th, k-th and l-th level, respectively, for the m-th replicate;

* μ̄ is the overall mean effect, αᵢ is the effect of the i-th level of the CT factor, βⱼ is the effect of the j-th level of the CC factor, γₖ is the effect of the k-th level of the DC factor and θ₁ is the effect of the l-th level of the W factor;

* ε<sub>ijklm</sub> is the error (or residual) component. It is assumed normally distributed with a constant variance SD<sup>2</sup><sub>ε</sub> and zero mean (μ<sub>ε</sub> = 0).

2.3.2. Regression Model

Once significant variables were identified by the ANOVA, a multilinear regression
analysis was conducted to develop predictive equations for both treated soils. The least
square method was employed to estimate the regression coefficients. Models for both
mixtures follow Equation (3).

$$
y_e = a_0 + a_1 \log_{10} \text{CT} + a_2 \text{CC} + a_3 \text{DC} + a_4 W \quad (3)
$$

where y<sub>e</sub> = f<sub>e(A)</sub> or f<sub>e(B)</sub> are the predicted ITS based on the experimental data in MPa.
Predictor variables are CT (days), CC (%), DC (%) and W (%); a<sub>i</sub>, with i = 0, 1, 2, 3, 4, as the
regression coefficients. First-order models were selected. A logarithmic relationship was
proposed to represent the evolution of CT, because its effect on CSSs usually follows this
pattern [2,11].

An examination of the residuals of both models allowed the construction of Probabilistic Density Functions (PDFs). The variations of residuals were considered to follow a normal distribution. The representation of variables as a PDF is convenient in geotechnical problems [50].
The Response Surface Method, RSM, [51] and contour plots were used to visualize the relationships between factors and their response according to the regression models.

## 2.4. Numerical Approach

### 2.4.1. Experimental Design

The inference space defined in the statistical approach (i.e., 7 < CT (days) < 360, 2 < CC (%) < 4, 94 < DC (%) < 98.5 and 90 < W (%) < 110) is the general frame of the numerical works. Sampling followed Latin Hypercube Sampling (LHS) [52]. This space-filling technique is based on DOE methods. One particularity of this random sampling method is that replications are avoided [29,53]. Besides this, it is computationally cheap and can cope with many input variables [54]. The use of LHS has been successfully proven for data generation in geotechnical problems when random variables are considered [55].

To match the laboratory work, the number of observations, n, varies from n = 10 to n = 100. To be statistically relevant, sampling was generated k = 10³ times. Based on a numerical toolbox for DOE [56], the LHS matrix (LHS<sub>n,k</sub>) was constructed. At the k-th iteration, b<sub>n,k</sub> is the n-th sample (vector) whose coordinates, i.e., mixture parameters, were generated by LHS. Then, B<sub>n,k</sub> is a data subset (matrix) that contains n samples for the k-th iteration. Therefore, the LHS matrix corresponds to the merged B<sub>n,k</sub> subsets, as described by Equation (4).

$$ \text{LHS}_{n,k} = \begin{bmatrix} B_{10,1} & \cdots & B_{10,k} \\ \vdots & \ddots & \vdots \\ B_{n,1} & \cdots & B_{10^2,10^3} \end{bmatrix}, \quad \text{where } B_{n,k} = \begin{bmatrix} b_{10,k} \\ \vdots \\ b_{n,k} \\ \vdots \\ b_{10^2,k} \end{bmatrix} \quad \text{and } b_{n,k} = (\text{CT}, \text{CC}, \text{DC}, \text{W})_{n,k} \quad (4) $$

### 2.4.2. Numerical ITS Values and Prediction Models

Numerical ITS values were generated following two steps: firstly, by applying the experimental prediction models f<sub>e</sub>, (Equation (3)) to the B<sub>n,k</sub> subsets. Second, a random residual value generated from the PDF was added. From there, numerical prediction models, named f<sub>n,k</sub>, i.e., f<sub>n,k(A)</sub> or f<sub>n,k(B)</sub>, were obtained by fitting the numerical ITS values to the corresponding B<sub>n,k</sub> sampling subset obtained by LHS.

### 2.4.3. Model Similarity

A similarity assessment was used to compare experimental and numerical predictions [29]. A first indicator δ<sub>n,k</sub>, defined as the absolute difference between numerical and experimental predictions for the B<sub>n,k</sub> subset, was used to quantify their distance. To generalize results, this parameter is normalized by the value of the experimental prediction on the same sampling, as presented in Equation (5). The indicator, δ<sub>n,k</sub>, is expressed in kPa and the normalized value δ<sub>n,k</sub> is expressed in percent.

$$ \overline{\delta_{n,k}} = \frac{\delta_{n,k}}{f_e(B_{n,k})}, \quad \text{where } \delta_{n,k} = |f_{n,k}(B_{n,k}) - f_e(B_{n,k})| \quad (5) $$

For k iterations, both values δ<sub>n,k</sub> and δ<sub>n,k</sub> are averaged and represented by a mean value, i.e., δ<sub>n</sub> and δ<sub>n</sub>, respectively. Both indicators are considered normally distributed with a constant standard deviation SD<sub>δ</sub> for δ<sub>n</sub> and SD<sub>δ-norm</sub> for δ<sub>n</sub>.

All calculations were performed using Python version 4.1 [57]. The following scientific libraries were used: SciPy [58], NumPy [59], pandas [60], statsmodels [61], doePy [56], matplotlib [62] and seaborn [63].
### 3. Results and Discussion

#### 3.1. *Experimental Results*

Hereinafter, Standard Deviation (SD) is indicated in brackets. The analysis of the experiments is based, for each case (V), on both the mean value and the SD (Table 3). Individual results can be consulted online [46]. A first analysis of ITS results shows:

<table><caption>Table 3. Experimental ITS results: mean ITS values and SD, in MPa.</caption><thead><tr><th rowspan="2">Material</th><th rowspan="2">V</th><th colspan="2">NT</th><th colspan="2">7 Days</th><th colspan="2">28 Days</th><th colspan="2">90 Days</th><th colspan="2">180 Days</th><th colspan="2">360 Days</th></tr><tr><th>Mean</th><th>SD</th><th>Mean</th><th>SD</th><th>Mean</th><th>SD</th><th>Mean</th><th>SD</th><th>Mean</th><th>SD</th><th>Mean</th><th>SD</th></tr></thead><tbody><tr><td rowspan="7">A</td><td>V1</td><td></td><td></td><td>0.116</td><td>0.008</td><td>0.154</td><td>0.009</td><td>0.152</td><td>0.014</td><td>0.162</td><td>0.013</td><td>0.230</td><td>0.010</td></tr><tr><td>V2</td><td></td><td></td><td>0.092</td><td>0.006</td><td>0.093</td><td>0.005</td><td>0.101</td><td>0.008</td><td>0.144</td><td>0.007</td><td>0.187</td><td>0.012</td></tr><tr><td>V3</td><td></td><td></td><td>0.140</td><td>0.019</td><td>0.151</td><td>0.009</td><td>0.176</td><td>0.010</td><td>0.226</td><td>0.003</td><td>0.215</td><td>0.031</td></tr><tr><td>V4</td><td></td><td></td><td>0.141</td><td>0.009</td><td>0.203</td><td>0.029</td><td>0.292</td><td>0.022</td><td>0.318</td><td>0.041</td><td>0.354</td><td>0.030</td></tr><tr><td>V5</td><td></td><td></td><td>0.125</td><td>0.009</td><td>0.155</td><td>0.009</td><td>0.204</td><td>0.006</td><td>0.211</td><td>0.014</td><td>-</td><td>-</td></tr><tr><td>V6</td><td></td><td></td><td>0.117</td><td>0.010</td><td>0.138</td><td>0.014</td><td>0.194</td><td>0.020</td><td>0.215</td><td>0.011</td><td>-</td><td>-</td></tr><tr><td>V7</td><td></td><td></td><td>0.144</td><td>0.016</td><td>0.162</td><td>0.009</td><td>0.218</td><td>0.004</td><td>0.272</td><td>0.032</td><td>0.258</td><td>0.011</td></tr><tr><td colspan="2">NT</td><td>0.019</td><td>0.004</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td rowspan="7">B</td><td>V1</td><td></td><td></td><td>0.105</td><td>0.011</td><td>0.135</td><td>0.009</td><td>0.179</td><td>0.019</td><td>0.198</td><td>0.014</td><td>0.264</td><td>0.013</td></tr><tr><td>V2</td><td></td><td></td><td>0.085</td><td>0.004</td><td>0.107</td><td>0.012</td><td>0.157</td><td>0.007</td><td>0.195</td><td>0.001</td><td>0.225</td><td>0.004</td></tr><tr><td>V3</td><td></td><td></td><td>0.120</td><td>0.007</td><td>0.174</td><td>0.008</td><td>0.210</td><td>0.014</td><td>0.261</td><td>0.013</td><td>0.313</td><td>0.008</td></tr><tr><td>V4</td><td></td><td></td><td>0.149</td><td>0.001</td><td>0.199</td><td>0.010</td><td>0.269</td><td>0.007</td><td>0.317</td><td>0.018</td><td>0.319</td><td>0.006</td></tr><tr><td>V5</td><td></td><td></td><td>0.108</td><td>0.004</td><td>0.140</td><td>0.018</td><td>0.190</td><td>0.009</td><td>0.252</td><td>0.010</td><td>-</td><td>-</td></tr><tr><td>V6</td><td></td><td></td><td>0.101</td><td>0.002</td><td>0.137</td><td>0.005</td><td>0.210</td><td>0.009</td><td>0.230</td><td>0.021</td><td>-</td><td>-</td></tr><tr><td>V7</td><td></td><td></td><td>0.126</td><td>0.002</td><td>0.172</td><td>0.017</td><td>0.253</td><td>0.008</td><td>0.306</td><td>0.017</td><td>0.344</td><td>0.007</td></tr><tr><td colspan="2">NT</td><td>0.010</td><td>0.001</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td></tr></tbody></table>

Considering all variabilities, the mean ITS values for both CSSs were between 0.085 and 0.354 MPa. Uncertainties on measurements, in terms of the SD, range between 0.001 and 0.041 MPa. The Coefficient Of Variation (COV), defined as the ratio between the SD and the corresponding mean ITS value, was under 15% for all CSSs. The results from the experiments can be considered as statistically relevant.

Measured ITS values were relevant according to the literature. For instance, the values for ITS on CSSs have been reported between 0.050 and 0.300 MPa for gravel and limestone [21], 0.050 and 0.550 MPa for silt [14], 0.010 and 0.370 MPa for sand [25] when CC varied from 4 to 8% (7 < CT < 21 days), 3% to 9% (7 < CT < 28 days) and 1 to 12% (CT = 7 days), respectively. For cement-modified loess, when treated with three different hydraulic binders (CC = 5.5% and 7.0%) ITS varied between 0.330 MPa and 0.530 MPa for CT = 90 days [13]. For cemented sands [64], ITS ranged between 0.020 and 0.130 MPa (42 < CT < 90 days) when 2 < CC < 6%. To compare with similar materials, a review of recycled cement-treated mixtures [65] found that the average ITS values were 0.500 MPa, with 75% of the individual values lower than 0.650 MPa. In the case of lime-treated soils [66], with binder content varying from 3 to 5%, ITS values were reported varying between 0.310 and 0.800 MPa (CT = 7 and 28 days). For mixed treatment, reported ITS values varied between 0.200 and 0.700 MPa when soils were treated with 5 to 6% of hydraulic binder (7 < CT < 720 days) [67].

A noteworthy strength improvement with cement addition and curing time. For CT = 7 days, the ratio between the ITS for V3 and NT was 7 for mixture A and 12 for mixture B. For CT = 360 days the same ratios were, respectively, 11 and 31. A number of researchers have reported the key role of CC and CT in improving the ITS on CSSs [13,21,25]. Specifically, the enhancing of quarry waste by adding small amounts of cement (CC = 2%) has been described in the literature [68].
Lower and higher values of ITS, in bold, seem to be related to CC (V2 and V4) and DC (V1 and V7). To counterbalance the use of lower CC by increasing CD, for instance, statistical models have been proposed in cement-treated materials [22].

The effects of W (V5 and V6) are not clear. The existence of a moisture content threshold where variations can improve the performance of CSSs is suggested in the literature [2].

## 3.2. Statistical Approach

### 3.2.1. Analysis of Individual Responses

The means model from Equation (2) is represented in Figure 2. The overall mean effect ($\bar{\mu}$), represented as dashed lines, show that ITS were slightly higher for mixture B than for mixture A. Effects of CT ($\alpha_i$), CC ($\beta_j$), DC ($\gamma_k$) and W ($\theta_l$) are shown as solid lines and their individual effects are related to a standardized slope.

<figure><img src="image_6.png" /><figcaption>**Figure 2.** Main effect plot for means. ITS (top) and COV (bottom) of CSSs as function of the studied factors: CT (days), CC (% dry matter), DC (%MDD) and W (%OMC).</figcaption></figure>

Figure 2 shows that ITS follows a quite linear function of Log10(CT) for both materials, as usual for CSSs [2]. With the higher slopes, Figure 2 suggests that the most relevant factor for B is CT. Nevertheless, with two different slopes, the effects of CT are probably different in each mixture, varying from short- to long-term. This could illustrate different kinetics over time, i.e., that the process of cementation could be different in the two materials. At short-term (CT = 7 days), ITS<sub>A</sub> > ITS<sub>B</sub>, whereas for CT > 28 days, ITS<sub>B</sub> > ITS<sub>A</sub>.

For B, the effects of CC come after the CT ones. For A, this effect is unclear, even if it is comparable with the slope from CT. Nevertheless, in both cases, measured ITS values as a function of CC can be described as essentially linear. This behavior has already been reported in the literature [12,14,69]. It is well known that CC plays an important role to improve the cohesion in CSSs. It can be noted that for A, the slope is higher than for B.

Regarding DC, mixture B is slightly more affected than A. A linear relationship is identified for both CC and DC. Finally, with slopes close to zero, variations around the overall mean due to W seem not to be determining.
The dispersion of results can also be discussed graphically. The coefficient of variation, COV, is plotted in Figure 2 (bottom). Regardless of the studied material, COV seems to be a parameter-independent, varying around constant values.

### 3.2.2. ANOVA

The results of the ANOVA are summarized in Table 4 for each source in terms of Sum of Squares (SS), F value ($F_o$), $p$-value and percent contribution (defined as the normalized SS) [48]. As degrees of freedom are one for each factor, the Mean Square Treatment (MST) are equal to SS.

<table><caption>**Table 4.** ANOVA results for the studied factors.</caption><thead><tr><th>CSS</th><th>Source of Variation</th><th>Degrees of Freedom</th><th>SS (×10<sup>-4</sup>)</th><th>F<sub>o</sub></th><th>p-Value</th><th>Percent Contribution</th></tr></thead><tbody><tr><td rowspan="5">A</td><td>Log<sub>10</sub> CT</td><td>1</td><td>1396.7</td><td>218.7</td><td>4.1 × 10<sup>-27</sup></td><td>44.5</td></tr><tr><td>CC</td><td>1</td><td>1514.0</td><td>237.0</td><td>2.3 × 10<sup>-28</sup></td><td>48.2</td></tr><tr><td>DC</td><td>1</td><td>225.6</td><td>35.3</td><td>3.9 × 10<sup>-8</sup></td><td>7.2</td></tr><tr><td>W</td><td>1</td><td>4.7</td><td>0.74</td><td>0.39</td><td>0.1</td></tr><tr><td>Error</td><td>102</td><td>651.5</td><td></td><td></td><td></td></tr><tr><td rowspan="5">B</td><td>Log<sub>10</sub> CT</td><td>1</td><td>3381.3</td><td>770.2</td><td>4.6 × 10<sup>-49</sup></td><td>77.3</td></tr><tr><td>CC</td><td>1</td><td>677.6</td><td>154.3</td><td>4.6 × 10<sup>-22</sup></td><td>15.5</td></tr><tr><td>DC</td><td>1</td><td>311.7</td><td>71.0</td><td>2.6 × 10<sup>-13</sup></td><td>7.1</td></tr><tr><td>W</td><td>1</td><td>2.7</td><td>0.62</td><td>0.43</td><td>0.1</td></tr><tr><td>Error</td><td>101</td><td>443.4</td><td></td><td></td><td></td></tr></tbody></table>

With the highest contributions in both cases, CT and CC are the main factors explaining ITS variations. For CT, the percent contribution was higher for B than for A (77.3 vs. 44.5). For CC, the percent contribution was higher for A than for B (48.2 vs. 15.5). The contribution of DC was around 7% for both materials.

The influence of W, in terms of percent contribution, was negligible for both materials. In the same way, $p$-values related to W are higher than $\alpha = 0.05$. The F value ($F_o$), is lower than the F critical value, $F_{crit}$ calculated from the degree of freedom, error and $\alpha$. For mixture A, ($F_o(W) = 0.39 < F_{crit}(1,102; 0.05) = 3.9$) and B, ($F_o(W) = 0.43 < F_{crit}(1,101; 0.05) = 3.9$). This result shows that for the studied inference space, this factor does not have a significant effect on ITS values for either mixtures A or B. It was excluded from the prediction models. All these results are in accordance with observations from Figure 2.

Nevertheless, the results of the ANOVA clarify the role of each parameter variation, whose scope is limited to the studied range. In a similar investigation, it has been reported that W was the most significant variable on ITS for CSS. However, variations were conducted over a wider range (60 to 140% of OMC) [21].

### 3.2.3. Multilinear Regression

Multilinear regression models (Equations (6) and (7)) allowed the prediction of ITS with high accuracy for both mixtures. For A, $R^2$ is 0.842 with a number of observations, n, of 107. For B, $R^2$ = 0.923, with n = 106. The accuracy of the models can be confirmed using the observed F value, $F_{obs}$, which is higher than the F critical value. For A, $F_{obs} = 183.2 > F_{crit}(3, 103;0.05) = 2.7$ and for B, $F_{obs} = 406.3 > F_{crit}(3, 102;0.05) = 2.7$. In Equations (6) and (7), CT is the curing time (days), CC is the cement content (%), and DC is the degree of compaction, as a percentage of the MDD.

$$ ITS_A (\text{MPa}) = (-131.99 + 6.34 \log_{10} \text{CT} + 6.90 \text{CC} + 1.24 \text{DC}) \times 10^{-2} \quad (6) $$

$$ ITS_B (\text{MPa}) = (-160.45 + 9.99 \log_{10} \text{CT} + 4.79 \text{CC} + 1.55 \text{DC}) \times 10^{-2} \quad (7) $$
By comparing the regression coefficients of Equations (6) and (7), some previous analyses from Figure 2 can be quantitatively confirmed: (i) the effects of CT are stronger for mixture B than for mixture A, (ii) the influence of CC is higher in mixture A and (iii) with a similar order of magnitude, the effects of DC are comparable for both mixtures. Coefficient values are in the same range for both materials, which implies comparable values of ITS regardless of soil type.

The experimental data and the predicted ITS values from Equations (6) and (7) are compared in Figure 3. Residual values (ε), defined in Equation (2), are the difference between measured and predicted values. They can be represented as the horizontal distances from each point to the line of equality. For both models, residuals are accurately described by a normal PDF (R² of 0.96 and 0.99 for mixtures A and B, respectively) with mean values close to zero (με (A) = -1.74 × 10⁻¹⁵ MPa and με (B) = -6.35 × 10⁻¹⁴ MPa), and a constant variance value (SD²ε (A) = (0.0248 MPa)² and SD²ε (B) = (0.0205 MPa)²). Thus, the underlying assumption of normality of the residual distributions was verified.

<figure><img src="image_7.png" /><figcaption>**Figure 3.** Comparison between experimental and predicted values of ITS for mixtures (a) A and (b) B.</figcaption></figure>

### 3.2.4. Coded Multilinear Regression

In order to determine which factor has the strongest effect on ITS, dimensionless models with coded units were considered. By a scale change on input variables, the magnitudes of the model's coefficients can be directly compared as they measure the effect of changing each design factor over the same one-unit interval [-1, +1] [48]. The construction of the coded variables is presented in Table 5.

<table><caption>**Table 5.** Natural and coded factors. Coded values in square brackets.</caption><thead><tr><th>Source of Variation</th><th>Min. Level</th><th colspan="3">Intermediate Level(s)</th><th>Max. Level</th></tr></thead><tbody><tr><td>CT</td><td>Log<sub>10</sub>(7)<br/>[-1.00]</td><td>Log<sub>10</sub>(28)<br/>[-0.30]</td><td>Log<sub>10</sub>(90)<br/>[+0.30]</td><td>Log<sub>10</sub>(180)<br/>[+0.65]</td><td>Log<sub>10</sub>(360)<br/>[+1.00]</td></tr><tr><td>CC</td><td>2.0<br/>[-1.00]</td><td colspan="3">3.0<br/>[0.00]</td><td>4.0<br/>[+1.00]</td></tr><tr><td>DC</td><td>94.0<br/>[-1.00]</td><td colspan="3">96.0<br/>[+0.11]</td><td>98.5<br/>[+1.00]</td></tr></tbody></table>
Coded equations, presented below, followed the same form that the reference model described in Equation (3). In the case of mixture A, coefficients in Equation (8) show that CC has the strongest effect (6.90 × 10⁻²) on ITS, followed by CT (5.44 × 10⁻²) and DC (2.78 × 10⁻²). For mixture B, (Equation (9)), the strongest effect is due to CT (8.56 × 10⁻²), followed by CC (4.80 × 10⁻²) and DC (3.50 × 10⁻²). The relative effect of every factor agrees with the percent contribution shown in Figure 2 and the ANOVA (Table 4). In all cases, the effects are positive.

$$ \text{ITS}_{C-A}(-) = (18.46 + 5.44\log_{10}\text{CT} + 6.90\text{CC} + 2.78\text{DC}) \times 10^{-2} \quad (8) $$

$$ \text{ITS}_{C-B}(-) = (20.46 + 8.56\log_{10}\text{CT} + 4.80\text{CC} + 3.50\text{DC}) \times 10^{-2} \quad (9) $$

### 3.2.5. Response Surfaces

As the varying factors that play a significant role on ITS are limited to 3, it is possible to visualize the prediction models (Equations (6) and (7)) as response surfaces (Figure 4).

<figure><img src="image_8.png" /><figcaption>Figure 4. Response surfaces at reference DC values (as percentage of MDD): mixture (a) A and (b) B.</figcaption></figure>

Responses surfaces (Figure 4) show the classical ITS behavior of CT for a given CC and DC. In both cases, it follows the same pattern: from quick short-term stabilization, to slow long-term effects. As pointed by Firoozi et al. [5], hydration in cement occurs quickly, which allows for an immediate strength gain at a young age. A reduction of the ITS gain at later ages for a given CC and DC can be linked to the deceleration of hydration [70] and a slow pozzolanic reaction.

For our discussion about the influence of parameters, one parameter is fixed and contour plots in 2D are generated from the response surfaces (Figures 5 and 6).

In Figures 5a and 6a, contours of ITS (iso-lines) as a function of CT and CC are plotted. Gradients represent the effect of CC on ITS as function of time for DC = 96% MDD. Despite the logarithmic shape of the ITS curves, ITS variation is constant and depends on CC variations. For a 1% increase of CC (e.g., 2.5% to 3.5%), the gain on ITS is around 0.069 MPa for A and 0.048 MPa for B regardless of CT. This suggests that CC governs the number of bonds between aggregates. The lower the CC, the lower the number of bonds and thus the lower ITS. However, it also suggests that the chemical kinetics were not modified.
<figure><img src="image_9.png" /><figcaption>**Figure 5.** ITS contour plots for mixture A: (a) CT vs. CC (DC = 96% MDD), (b) CT vs. DC (CC = 3%), (c) DC vs. CC (CT = 90 days). Contour lines: ITS, in MPa.</figcaption></figure>

<figure><img src="image_10.png" /><figcaption>**Figure 6.** ITS contour plots for mixture B: (a) CT vs. CC (DC = 96% MDD), (b) CT vs. DC (CC = 3%), (c) DC vs. CC (CT = 90 days). Contour lines: ITS, in MPa.</figcaption></figure>

The coarse material B (Figure 6a) exhibits higher gradients than the clayey mixture A (Figure 5a) under the same preparation conditions. This may be linked to different stabilization mechanisms. Particularly, clay content and OMC can explain these results. The higher occurrence of water in A could explain the quick reactivity with cement in the short-term. Conversely, the occurrence of clay could affect the long-term quality of cement bonds between soil aggregates. It is accepted that the best cement-stabilization results are observed on silt and coarse-grained materials [9] whereas fine-grained soil may require a higher CC [27]. For example, for a fixed CC, at DC = 96% MDD, with CT varying from 7 to 360 days, the increase in ITS is 0.109 MPa for mixture A and 0.171 MPa for mixture B.

To achieve a required ITS value with lower quantities of cement, it is necessary to have a higher CT. For instance, in Figure 5a for material A, to obtain ITS = 0.200 MPa we consider both CT = 360 days and CC = 2.5%, or CT = 7 days and CC = 4%. This is an important point, especially in the short term, when ITS values are critical issues on practice (e.g., opening to traffic for road structures). Finally, this relationship can be used to support the choice of smaller amounts of binder in applications where the CT is not a critical factor. It also could be used, in the field, for estimating mechanical performance variations due to dosage fluctuations.

Figures 5b and 6b plot ITS contours as a function of CT and DC. Regardless of the materials, contour plots present the same pattern previously described for the interaction between CT and CC. For instance, when the DC ranges from 94.0 to 98.5% of MDD (at CC = 3%), the gain on ITS is 0.056 MPa for A and 0.070 for B, regardless of CT.
In the same way, a higher DC would lead to achieving a given value of ITS with a lower CT. Specifically, the effects of the DC are stronger for mixture B. In Figure 6b (CC = 3%), for B, for example, to obtain ITS = 0.250 MPa we can consider both CT = 65 days and DC = 98.5% MDD, or CT = 320 days and DC = 94% MDD. The mechanism might be the following: the higher the DC, the higher the contact points between aggregates in the matrix and thus the creation of more cemented bonds. With a higher importance of DC for the coarse material (mixture B), this suggestion is in accordance with previous observations [71]. These remarks support the importance of compaction control and their consequence even in long term.

Figures 5c and 6c plot ITS contours as a function of DC and CC. As the effects of both DC and CC are linear, the contour plots are straight lines. In mixture B at CT = 90 days (Figure 6c), the contour plots show that a reduction of 1% on CC (e.g., from 3 to 2%) can be compensated by increasing the DC on 3.1% (e.g., from 94% to 97.1% MDD). In practice, cement savings could be expected with good compaction.

A practical example to use these contour plots could be the following. For a given level of expected ITS performance (e.g., material A, ITS = 0.200 MPa), several combinations are possible. For instance, for regular field practices (CC = 3%, DC = 96%), this performance is achieved normally for CT = 100 days (Figure 5a). The curing time can be reduced by increasing DC for the same CC of 3% (CT = 40 days when DC = 98%, in Figure 5b). Increasing the DC also allows CC to be reduced. For a given curing time (e.g., CT = 90 days, Figure 5c), varying DC from 96% to 98% allows CC to be reduced from 3% to 2.7% to achieve the aimed performance.

### 3.3. *Numerical Approach*

An example of one numerical experiment for the mixture A is illustrated in Figure 7, for a set of n = 100 observations at the k = 10³ iteration ($B_{10^2,10^3}$). Numerical ITS values are generated and presented as a function of dosage variables. In Figure 7, the hue and marker size represent the numerical ITS values. The LHS sampling is illustrated with the points scattering. In spite of the addition of noise, the effects of dosage variables on ITS can be identified. When compared with iso-lines from the contour plot of A (Figure 5), the ITS increase due to CT, CC and DC is clear. The introduction of noise on the numerical data generates variations around the reference model.

<figure><img src="image_11.png" /><figcaption>**Figure 7.** Numerical ITS experiment for mixture A ($B_{10^2,10^3}$). Iso-lines for ITS, in MPa, are generated at: (a) DC = 96% MDD, (b) CC = 3% and (c) CT = 90 days. Marker size and color hue refer to numerical values of ITS.</figcaption></figure>

The parameter $\delta_{n,k}$ as defined in Equation (5) represents the distance between the noised numerical model and the experimental model. Once all the iterations are complete, the mean ($\delta_n$) and normalized ($\bar{\delta}_n$) indicators, and their respective standard deviations, can
be calculated. Figure 8 illustrates that for both materials, the COV of $\delta_n$ and $\bar{\delta}_n$ are nearly constant (in average, 37%) which implies a proportional reduction of the SD as function of n for both indicators.

<figure><img src="image_12.png" /><figcaption>**Figure 8.** Variation of (a) $\delta_n$, its COV, and (b) $\bar{\delta}_n$ with its COV as function of number of observations ($k = 10^3$).</figcaption></figure>

The higher the number of samples n, the lower $\delta_n$ as shown in Figure 8a. For the studied range ($10 < n < 100$), $\delta_n$ decreases from 12.1 (4.6) to 3.9 (1.4) kPa for A and from 10.0 (3.9) to 3.2 (1.1) kPa for mixture B. That is, there is a total net gain on accuracy of 8.2 kPa for A and 6.8 kPa for B. The lower values of $\delta_n$ for B than for A are related to the error of experimental residuals, which is lower for B than for A ($SD^2_{\epsilon}(B) < SD^2_{\epsilon}(A)$). The numerical results are in accordance with the experimental ones. This is clearer in Figure 8b with the normalized indicator $\bar{\delta}_n$ (Equation (5)). For the range of observations analyzed ($10 < n < 100$), $\bar{\delta}_n$ decreases from 6.8 (2.8)% to 2.0 (0.7)% for A and from 4.8 (2.1)% to 1.4 (0.5)% for B material. For both indicators and both materials, the rate of decline slows as n increases (Figure 8).

A curve that represents the gain in accuracy for both materials is presented in Figure 9. The net gain for a given n is described as the ratio between $\delta_{10} - \delta_n$ and $\delta_{10} - \delta_{100}$. For the studied range of n, regardless of the material, both curves exhibit the same increasing pattern as a function of n. Most of the improvement (82%) is achieved when n varies from 10 to 50 observations. That corresponds approximately to half of the experimental sampling. Moreover, 43% of the total net gain is obtained when n is varied from 10 to 20 observations. This implies that good accuracy in the numerical models can be rapidly achieved, e.g., with less than 50 samples by LHS. This result can be used to optimize the sampling of experimental procedures.
<figure><img src="image_13.png" /><figcaption>Figure 9. Net gain as a function of number of observations for studied mixtures.</figcaption></figure>

## 4. Conclusions

We have experimentally and numerically assessed the impact of dosage in two CSSs in a fixed inference space. Based on a conventional approach, an experimental variability study on ITS was planned for both materials. Curing time, cement content, degree of compaction and water content were the dosage variables. An ANOVA was used to identify the significant dosage variables and to build prediction models. A numerical approach was used to iteratively generate synthetic ITS values from statistical prediction models. An LHS space-filling technique was used. Experimental and numerical predictions were compared. Based on the experimental results, statistical analysis and numerical results, within the inference space, the following conclusions can be drawn:

* The addition of a few percentage points of cement significantly increased the overall performance of both sandy and clayey soils. Mean treated ITS values for the reference mixture varied from 7 up to 11 times the non-treated ITS for mixture A, and from 12 up to 31 times for mixture B.

* The content of water variations in the interval we studied (±10% OMC) showed no significant effects on ITS values for both mixtures. Once this variable was excluded, density was identified as the factor with the lowest effect. Conversely, CC and CT proved to be most significant variable on mixture A and B, respectively.

* The experimental ITS was described accurately by using dosage variables and curing time as predictors on multilinear regression models: R² of 0.84 and 0.92 for mixture A and B, respectively. The effects of CC were stronger on the clayey soil (A) whereas the influence of CT and DC showed to be higher on the sandy mixture (B).

* The ANOVA proved to be a straightforward method to construct experimental prediction models. These models can be used to interpolate conditions that were not explicitly evaluated in the experimental program. Preparation parameters are ranked regarding their effect on ITS values.

* The combined effects of dosage variables were assessed by means of contour plots. These tools can be used for mixture optimization according to technical specifications and operating conditions.
* The combined effects of dosage variables were observed on the response of surfaces and associated contour plots. These tools can be used to assess specifications for mixture optimization, and operating conditions from short to long-term.

* Long-term results can be achieved in several days by increasing CC. Reducing the binder consumption can be achieved by increasing CT or DC. Specifically, the effects of compaction were shown to be relevant even in the long-term in both materials.

* The residual values of experimental models are well described with normal PDFs. These statistical descriptions were used to generate numerical data.

* LHS as space-filling technique was a powerful method to generate efficiently numerical sampling data.

* As expected, the number of observations was critical in explaining the differences between experimental and numerical models. The results showed that regardless of material type, the complete sampling size (n~100) may not be necessary for achieving a similar degree of accuracy. A net accuracy gain of 43% was measured when the number of observations varied from 10 to 20 for both mixtures. This gain reached 82% in both cases when n varied from 10 to 50. This supports the optimization of the experimental work in terms of the reduction of sampling. Moreover, the proposed methodology associates a number of experimental observations with a given degree of uncertainty on ITS.

* Finally, the proposed methodology can be generalized to other experimental scientific fields.

**Author Contributions:** Conceptualization, T.L.; methodology, M.C.-L. and T.L.; validation, L.T.; formal analysis, M.C.-L.; investigation, M.C.-L. and T.L.; resources, T.L. and J.-P.S.; writing—original draft preparation, M.C.-L.; writing—review and editing, T.L. and L.T.; supervision, T.L. and L.T.; project administration, T.L. and J.-P.S.; funding acquisition, J.-P.S. All authors have read and agreed to the published version of the manuscript.

**Funding:** This research was funded by the French Ministry for the Ecological and Inclusive Transition and the COLAS company, grant number RP2-E20029.

**Data Availability Statement:** Castaneda-Lopez, M. A., Lenoir, T., Sanfratello, J. P., & Thorel, L. Tensile properties of cement stabilized soils: experimental data. Recherche Data Gouv. 2023. https://doi.org/10.57745/DFW5LS.

**Acknowledgments:** Authors thank the French Ministry for the Ecological and Inclusive Transition and the COLAS company for their funding. The manuscript benefited from comments and suggestions by four anonymous reviewers.

**Conflicts of Interest:** The authors declare no conflict of interest. The funder had no role in the design of the study; in the collection, analyses, or interpretation of data; in the writing of the manuscript; or in the decision to publish the results.

## References

1. Preteseille, M.; Lenoir, T. Mechanical Fatigue of a Stabilized/Treated Soil Tested with Uniaxial and Biaxial Flexural Tests. *Transp. Res. Procedia* **2016**, *14*, 1923–1929. [CrossRef]

2. Xuan, D.X.; Houben, L.J.M.; Molenaar, A.A.A.; Shui, Z.H. Mechanical properties of cement-treated aggregate material—A review. *Mater. Des.* **2012**, *33*, 496–502. [CrossRef]

3. Preteseille, M.; Lenoir, T. Structural test at the laboratory scale for the utilization of stabilized fine-grained soils in the subgrades of High Speed Rail infrastructures: Experimental aspects. *Int. J. Fatigue* **2016**, *82*, 505–513. [CrossRef]

4. Celauro, B.; Bevilacqua, A.; Lo Bosco, D.; Celauro, C. Design Procedures for Soil-Lime Stabilization for Road and Railway Embankments. Part 1—Review of Design Methods. *Procedia Soc. Behav. Sci.* **2012**, *53*, 754–763. [CrossRef]

5. Firoozi, A.A.; Guney Olgun, C.; Firoozi, A.A.; Baghini, M.S. Fundamentals of soil stabilization. *Int. J. Geo Eng.* **2017**, *8*, 26. [CrossRef]

6. Pedarla, A.; Chittoori, S.; Puppala, A. Influence of mineralogy and plasticity index on the stabilization effectiveness of expansive clays. *Transp. Res. Rec.* **2011**, *2212*, 91–99. [CrossRef]

7. Sirivitmaitrie, C.; Puppala, A.; Saride, S.; Hoyos, L. Combined lime-cement stabilization for longer life of low-volume roads. *Transp. Res. Rec.* **2011**, *2204*, 140–147. [CrossRef]

8. Chittoori, B.C.S.; Puppala, A.J.; Pedarla, A. Addressing Clay Mineralogy Effects on Performance of Chemically Stabilized Expansive Soils Subjected to Seasonal Wetting and Drying. *J. Geotech. Geoenviron. Eng.* **2018**, *144*, 1–12. [CrossRef]
9. Sariosseiri, F.; Muhunthan, B. Effect of cement treatment on geotechnical properties of some Washington State soils. *Eng. Geol.* **2009**, *104*, 119–125. [CrossRef]

10. Consoli, N.C.; Vaz Ferreira, P.M.; Tang, C.S.; Veloso Marques, S.F.; Festugato, L.; Corte, M.B. A unique relationship determining strength of silty/clayey soils—Portland cement mixes. *Soils Found.* **2016**, *56*, 1082–1088. [CrossRef]

11. Horpibulsuk, S.; Miura, N.; Nagaraj, T.S. Assessment of strength development in cement-admixed high water content clays with Abrams' law as a basis. *Geotechnique* **2003**, *53*, 439–444. [CrossRef]

12. Schnaid, F.; Prietto, P.D.M.; Consoli, N.C. Characterization of Cemented Sand in Triaxial Compression. *J. Geotech. Geoenviron. Eng.* **2001**, *127*, 857–868. [CrossRef]

13. Lenoir, T.; Dubreucq, T.; Lambert, T.; Killinger, D. Safety factor calculation of a road structure with cement-modified loess as subgrade. *Transp. Geotech.* **2021**, *30*, 100604. [CrossRef]

14. De Baldovino, J.J.A.; dos Izzo, R.L.S.; Pereira, M.D.; de Rocha, E.V.G.; Rose, J.L.; Bordignon, V.R. Equations Controlling Tensile and Compressive Strength Ratio of Sedimentary Soil-Cement Mixtures under Optimal Compaction Conditions. *J. Mater. Civ. Eng.* **2020**, *32*, 04019320. [CrossRef]

15. Feiz, R.; Ammenberg, J.; Baas, L.; Eklund, M.; Helgstrand, A.; Marshall, R. Improving the CO2 performance of cement, part II: Framework for assessing CO2 improvement measures in the cement industry. *J. Clean. Prod.* **2015**, *98*, 282–291. [CrossRef]

16. Andrew, R.M. Global CO2 emissions from cement production. *Earth Syst. Sci. Data* **2018**, *10*, 195–217. [CrossRef]

17. Wang, M.; Huston, M.T. Direct-tensile stress and strain of a cement-stabilized soil. *Highw. Res. Rec.* **1972**, *379*, 19–24.

18. Consoli, N.C.; Dalla Rosa, J.A.; Gauer, E.A.; dos Santos, V.R.; Moretto, R.L.; Corte, M.B. Key parameters for tensile and compressive strength of silt-lime mixtures. *Géotech. Lett.* **2012**, *2*, 81–85. [CrossRef]

19. MolaAbasi, H.; Khajeh, A.; Semsani, S.N.; Kordnaei, A. Prediction of Zeolite-Cemented Sand Tensile Strength by GMDH type Neural Network. *J. Adhes. Sci. Technol.* **2019**, *33*, 1611–1625. [CrossRef]

20. Anagnos, J.N.; Kennedy, T.W.; Hudson, W.R. *Evaluation and Prediction of Tensile Properties of Cement-Treated Materials*; Research Report Number 98-8; Center for Highway Research—University of Texas: Austin, TX, USA, 1970.

21. Moore, R.K.; Kennedy, T.W.; Hudson, W.R. Factors affecting the tensile strength of materials. *Highw. Res. Rec.* **1970**, *315*, 64–80.

22. Fedrigo, W.; Núñez, W.P.; Castañeda López, M.A.; Kleinert, T.R.; Ceratti, J.A.P. A study on the resilient modulus of cement-treated mixtures of RAP and aggregates using indirect tensile, triaxial and flexural tests. *Constr. Build. Mater.* **2018**, *171*, 161–169. [CrossRef]

23. Adresi, M.; Khishdari, A.; Ahmadi, A.; Rooholamini, H. Influence of high content of reclaimed asphalt on the mechanical properties of cement-treated base under critical environmental conditions. *Int. J. Pavement Eng.* **2019**, *20*, 1098–1105. [CrossRef]

24. Consoli, N.C.; Foppa, D.; Festugato, L.; Heineck, K.S. Key Parameters for Strength Control of Artificially Cemented Soils. *J. Geotech. Geoenviron. Eng.* **2007**, *133*, 197–205. [CrossRef]

25. Consoli, N.C.; Cruz, R.C.; Floss, M.F.; Festugato, L. Parameters Controlling Tensile and Compressive Strength of Artificially Cemented Sand. *J. Geotech. Geoenviron. Eng.* **2010**, *136*, 759–763. [CrossRef]

26. Zhang, J.; Deng, A.; Jaksa, M. Optimizing micaceous soil stabilization using response surface method. *J. Rock. Mech. Geotech. Eng.* **2021**, *13*, 212–220. [CrossRef]

27. Ordoñez Muñoz, Y.; Luis dos Santos Izzo, R.; Leindorf de Almeida, J.; Arrieta Baldovino, J.; Lundgren Rose, J. The role of rice husk ash, cement and polypropylene fibers on the mechanical behavior of a soil from Guabirotuba formation. *Transp. Geotech.* **2021**, *31*, 100673. [CrossRef]

28. LCPC; SETRA. *Trajetoyment Des Sols à la Chaux et/ou aux Liants Hydrauliques*, 2nd ed.; LCPC: Paris, France, 2000.

29. Ma, S.; Henry, M.; Opon, J. Exploring the experimental design and statistical modeling of cementitious composite systems using various sampling methods. *J. Adv. Concr. Technol.* **2021**, *19*, 501–518. [CrossRef]

30. Slebi-Acevedo, C.J.; Castro-Fresno, D.; Pascual-Muñoz, P.; Lastra-González, P. A combination of DOE—Multi-criteria decision making analysis applied to additive assessment in porous asphalt mixture. *Int. J. Pavement Eng.* **2022**, *23*, 2489–2502. [CrossRef]

31. Ikeagwuani, C.C.; Agunwamba, J.C.; Nwankwo, C.M.; Eneh, M. Additives optimization for expansive soil subgrade modification based on Taguchi grey relational analysis. *Int. J. Pavement Res. Technol.* **2021**, *14*, 138–152. [CrossRef]

32. Olgun, M. The effects and optimization of additives for expansive clays under freeze-thaw conditions. *Cold Reg. Sci. Technol.* **2013**, *93*, 36–46. [CrossRef]

33. Shahbazi, M.; Rowshanzamir, M.; Abtahi, S.M.; Hejazi, S.M. Optimization of carpet waste fibers and steel slag particles to reinforce expansive soil using response surface methodology. *Appl. Clay Sci.* **2017**, *142*, 185–192. [CrossRef]

34. Chen, K.; Wu, D.; Zhang, Z.; Pan, C.; Shen, X.; Xia, L.; Zang, J. Modeling and optimization of fly ash-slag-based geopolymer using response surface method and its application in soft soil stabilization. *Constr. Build. Mater.* **2022**, *315*, 125723. [CrossRef]

35. Xuan, D.X.; Houben, L.J.M.; Molenaar, A.A.A.; Shui, Z.H. Mixture optimization of cement treated demolition waste with recycled masonry and concrete. *Mater. Struct. Constr.* **2012**, *45*, 143–151. [CrossRef]

36. Longarini, N.; Crespi, P.; Zucca, M.; Giordano, N.; Silvestro, G. The advantages of fly ash use in concrete structures. *Inz. Miner.* **2014**, *15*, 141–145.

37. NF EN ISO 17892-4; Geotechnical Investigation and Testing—Laboratory Testing of Soil—Part 4: Determination of Particle Size Distribution. AFNOR: Saint-Denis, France, 2018.

38. NF EN ISO 17892-3; Geotechnical Investigation and Testing—Laboratory Testing of Soil—Part 3: Determination of Particle Density. AFNOR: Saint-Denis, France, 2015.
39. NF EN ISO 17892-12; Geotechnical Investigation and Testing—Laboratory Testing of Soil—Part 12: Determination of Liquid and Plastic Limits. AFNOR: Saint-Denis, France, 2018.

40. M 145-91; Classification of Soils and Soil-Aggregate Mixtures for Highway Construction Purposes. AASHTO: Washington, DC, USA, 2021.

41. NF P11-300; Classification of Materials for Use in the Construction of Embankments and Capping Layers of Road Infrastructures. AFNOR: Saint-Denis, France, 1992.

42. D2487-17; Standard Practice for Classification of Soils for Engineering Purposes (Unified Soil Classification System). ASTM: West Conshohocken, PA, USA, 2017.

43. NF EN ISO 14688-1; Geotechnical Investigation and Testing—Identification and Classification of Soil—Part 1: Identification and Description. AFNOR: Saint-Denis, France, 2018.

44. NF EN 197-1; Cement—Part 1: Composition, Specifications and Conformity Criteria for Common Cements. AFNOR: Saint-Denis, France, 2012; p. 42.

45. NF P94-093; Soils: Investigation and Testing—Determination of the Compaction Reference Values of a Soil Type—Standard Proctor Test—Modified Proctor Test. AFNOR: Saint-Denis, France, 2014.

46. Castaneda-Lopez, M.A.; Lenoir, T.; Sanfratello, J.P.; Thorel, L. Tensile properties of cement stabilized soils: Experimental data. *Rech. Data Gouv.* **2023**. [CrossRef]

47. NF EN 13286-42; Unbound and Hydraulically Bound Mixtures—Part 42: Test Method for the Determination of the Indirect Tensile Strength of Hydraulically Bound Mixtures. AFNOR: Saint-Denis, France, 2003.

48. Montgomery, D.C. *Design and Analysis of Experiments*, 9th ed.; Wiley: Hoboken, NJ, USA, 2017.

49. Abdallah, A.; Russo, G.; Cuisinier, O. Statistical and Predictive Analyses of the Strength Development of a Cement-Treated Clayey Soil. *Geotechnics* **2023**, *3*, 465–479. [CrossRef]

50. Johari, A.; Fooladi, H. Comparative study of stochastic slope stability analysis based on conditional and unconditional random field. *Comput. Geotech.* **2020**, *125*, 103707. [CrossRef]

51. Myers, R.H.; Montgomery, D.C.; Anderson-Cook, C.M. *Response Surface Methodology: Process and Product Optimization Using Designed Experiments*, 3rd ed.; Wiley: Hoboken, NJ, USA, 2009.

52. Viana, F.A.C. A Tutorial on Latin Hypercube Design of Experiments. *Qual. Reliab. Eng. Int.* **2016**, *32*, 1975–1985. [CrossRef]

53. Joseph, V.R. Space-filling designs for computer experiments: A review. *Qual. Eng.* **2016**, *28*, 28–35. [CrossRef]

54. Sacks, J.; Welch, W.J.; Mitchell, T.J.; Wynn, H.P. Design and analysis of computer experiments. *Stat. Sci.* **1989**, *4*, 409–423. [CrossRef]

55. Zhang, H.; Wu, S.; Zhang, Z.; Huang, S. Reliability analysis of rock slopes considering the uncertainty of joint spatial distributions. *Comput. Geotech.* **2023**, *161*, 105566. [CrossRef]

56. Sarkar, T.; Flanagan, D. doepy 0.0.1. 2023. Available online: https://pypi.org/project/doepy/ (accessed on 29 March 2023).

57. Van Rossum, G. *Python Tutorial CS-R9526*; Centrum voor Wiskunde en Informatica (CWI): Amsterdam, The Netherlands, 1995.

58. Virtanen, P.; Gommers, R.; Oliphant, T.E.; Haberland, M.; Reddy, T.; Cournapeau, D.; Burovski, E.; Peterson, P.; Weckesser, W.; Bright, J.; et al. SciPy 1.0: Fundamental algorithms for scientific computing in Python. *Nat. Methods* **2020**, *17*, 261–272. [CrossRef] [PubMed]

59. Harris, C.R.; Millman, K.J.; van der Walt, S.J.; Gommers, R.; Virtanen, P.; Cournapeau, D.; Wieser, E.; Taylor, J.; Berg, S.; Smith, N.J.; et al. Array programming with NumPy. *Nature* **2020**, *585*, 357–362. [CrossRef] [PubMed]

60. McKinney, W. Data Structures for Statistical Computing in Python. In *Proceedings of the 9th Python in Science Conference (SciPy 2010)*, Austin, TX, USA, 28 June–3 July 2010; pp. 56–61. [CrossRef]

61. Seabold, S.; Perktold, J. *Statsmodels: Econometric and Statistical Modeling with Python*. In *Proceedings of the 9th Python in Science Conference (SciPy 2010)*, Austin, TX, USA, 28 June–3 July 2010; pp. 92–96. [CrossRef]

62. Hunter, J.D. Matplotlib: A 2D graphics environment. *Comput. Sci. Eng.* **2007**, *9*, 90–95. [CrossRef]

63. Waskom, M. Seaborn: Statistical data visualization. *J. Open Source Softw.* **2021**, *6*, 3021. [CrossRef]

64. Mola-Abasi, H.; Khajeh, A.; Naderi Semsani, S. Effect of the Ratio between Porosity and SiO2 and Al2O3 on Tensile Strength of Zeolite-Cemented Sands. *J. Mater. Civ. Eng.* **2018**, *30*, 04018028. [CrossRef]

65. Fedrigo, W.; Núñez, W.P.; Visser, A.T. A review of full-depth reclamation of pavements with Portland cement: Brazil and abroad. *Constr. Build. Mater.* **2020**, *262*, 120540. [CrossRef]

66. Piratheepan, J.; Gnanendran, C.T.; Arulrajah, A. Determination of $c$ and $\varphi$ from IDT and Unconfined Compression Testing and Numerical Analysis. *J. Mater. Civ. Eng.* **2012**, *24*, 1153–1164. [CrossRef]

67. Preteseille, M. Comportement à la Fatigue des Sols Traités aux Liants Hydrauliques Dans les Plates-Formes des Structures Ferroviaires Pour L.G.V: Modélisations Numérique et Expérimentale de Leur Comportement. Ph.D. Thesis, Ecole Centrale de Nantes, Nantes, France, 2014.

68. Mwumvaneza, V.; Hou, W.; Ozer, H.; Tutumluer, E.; Al-Qadi, I.L.; Beshears, S. Characterization and stabilization of quarry byproducts for sustainable pavement applications. *Transp. Res. Rec.* **2015**, *2509*, 1–9. [CrossRef]

69. Consoli, N.C.; Fonseca, A.V.; Silva, S.; Cruz, R.C.; Fonini, A. Parameters controlling stiffness and strength of artificially cemented soils. *Géotechnique* **2012**, *62*, 177–183. [CrossRef]
70. Bullard, J.W.; Jennings, H.M.; Livingston, R.A.; Nonat, A.; Scherer, G.W.; Schweitzer, J.S.; Scrivener, K.L.; Thomas, J.J. Mechanisms of cement hydration. *Cem. Concr. Res.* **2011**, *41*, 1208–1223. [CrossRef]

71. Kezdi, A.; Rethati, L. *Soil Mechanics of Earthworks, Foundations and Highway Engineering*; Elsevier: Amsterdam, The Netherlands, 1988. [CrossRef]

**Disclaimer/Publisher's Note:** The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.