# Predicting the settlement of coarse granular materials under vertical loading

Juan Carlos Quezada, Gilles Saussine, Pierre Breul, Farhang Radjai

► To cite this version:

Juan Carlos Quezada, Gilles Saussine, Pierre Breul, Farhang Radjai. Predicting the settlement of coarse granular materials under vertical loading. Scientific Reports, 2015, 4, pp.5707. 10.1038/srep05707 . hal-02093091

HAL Id: hal-02093091

https://hal.science/hal-02093091

Submitted on 8 Apr 2019

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 Attribution 4.0 International License
OPEN

SUBJECT AREAS:

CIVIL ENGINEERING

FLUIDS

APPLIED PHYSICS

Received
21 March 2014

Accepted
27 June 2014

Published
16 July 2014

Correspondence and
requests for materials
should be addressed to
J.C.Q. (juan-carlos.quezada@univ-mont2.fr)

# Predicting the settlement of coarse granular materials under vertical loading

Juan Carlos Quezada¹, Gilles Saussine², Pierre Breul³ & Farhang Radjai⁴

¹SNCF PSIG, 6 Avenue François Mitterrand, 93574 La Plaine Saint Denis Cedex, France, ²SNCF Innovation and Research, 40 Avenue des Terroirs de France, F-75611 Paris Cedex 12, France, ³Institut Pascal, CNRS-Université Clermont II, F-63171 Aubière, France, ⁴Université Montpellier 2, CNRS, LMGC, F-34095 Montpellier, France.

Granular materials are widely used in industrial processes despite their complex and poorly understood mechanical behaviour both in static and dynamic regimes. A prototypical example is the settlement and compaction of a granular bed under vibrational loading. The elementary mechanisms of this process are still unclear and there is presently no established theory or methodology to predict the settlement and its statistical variability. By means of a parametric study, carried out on a full-scale track, and a critical analysis of density relaxation laws, we introduce a novel settlement model in coarse granular materials under cyclic loading. Our extensive experimental data indicate that the settlement process is governed by three independent parameters strongly correlated with the vibration intensity and initial packing fraction. We show that the mean settlement is well predicted by the model with its parameter values extracted from experimental data.

Granular materials are both pressure-dependent and density-dependent materials and exhibit a broad range of intricate behaviours due to their discrete nature, dissipative interactions and generic structural disorder¹. The packing fraction may vary as a result of particle rearrangements induced by shearing or vibrations and it leads to dramatic changes in the structure and mechanical response of a granular material²⁻⁸. A long-time logarithmic relaxation law of the packing fraction is systematically observed in experiments⁹,¹⁰. In simple compaction models, this behaviour is attributed to the exponentially increasing time for the particles to reach a new configuration of lower packing fraction. The case of settlement under cyclic loading has, however, been much less investigated. The settlement of granular bed occurs due to both compaction and side-wise spreading. An important industrial example is the railway ballast, which undergoes gradual settlement under the static and dynamic overloads induced by train traffics¹¹⁻¹⁴. The readjustment of differential settlements requires costly operations on fast-train railways. For this reason, an improved understanding of the parameters governing the settlement process is a critical technological challenge for new developments in this field.

In this paper, we show that the total settlement $\tau_N$ under vertical cyclic loading is governed by a logarithmic relaxation law as a function of the number N of cycles:

$$ \tau_N = \tau_\infty \left( 1 - \frac{1}{1 + B \ln \left( 1 + \frac{N}{N_0} \right)} \right) \quad (1) $$

where the three fitting parameters $\tau_\infty$, B and $N_0$ can be evaluated from the loading parameters, namely the frequency $\omega$ (related to the train speed for ballast) and initial packing fraction of the material. Our experimental correlations between model and loading parameters show consistently that $\tau_\infty$ and B depend on the dimensionless loading intensity $\Gamma = (A\omega^2)/(pd^2/m + g)$, where A is the vibration amplitude, p is the confining pressure (under the sleeper for ballast), d is the average particle diameter, m is the average particle mass and g is gravitational acceleration. We also find that the parameter $N_0$ is linked to the initial packing fraction of the material. In fact, this parameter controls the initial settlement rate, and it was systematically determined by means of a light penetrometer in our experiments on ballast material.
<figure><img src="image_2.png" /><figcaption>Figure 1 | Full-scale track model used in the present work.</figcaption></figure>

## Results

**Loading system.** We carried out several experimental tests on a full-scale track model composed of a granular bed of ballast particles (Fig. 1). The bed was subjected to a sinusoidally varying vertical overload up to $10^4$ cycles for a broad range of amplitudes and frequencies. In the same way, the initial state of the granular bed was varied and a light dynamic penetrometer was used to characterize its packing fraction before each loading test$^{15–17}$. This device provides the cone penetration strength, which is a parameter strongly dependent on the packing fraction of a granular sample. The principle of the light penetrometer is to push into the material, by manual beating with a hammer, a conic tip of cross-section $A_{tip}=2 \text{ cm}^2$. For each blow, the depth of penetration and the input energy are measured to calculate the dynamic cone strength $q_d$ with the corresponding depth using the so-called Dutch formula:

$$q_d = \frac{E_K M}{A_{tip} \times e M + P} \quad (2)$$

where M is the mass of the hammer, P is the weight of the driven parts during impact (impact head, road and tip), $E_K$ is the kinetic energy transmitted to the system during the impact and e is the penetration per blow. A detailed description of the model and testing apparatus can be found in the Methods section below.

The cumulative settlement may qualitatively be described as occurring in three stages as displayed in Fig. 2a. The first stage is the initial fast compaction of the granular bed during the first 200 cycles. Its quasi-linear evolution reflects the uneven distribution of voids in the material. The largest voids disappear as a result of local particle rearrangements, leading thus to the overall increase of packing fraction and settlement. The second stage lasts from 200 to 6000

<figure><img src="image_3.png" /><figcaption>Figure 2 | (a) Settlement of ballast material as a function of the number of cycles with its three stages; (b) Comparison between the experimental settlement data obtained from a cyclic loading test with frequency of 3.3 Hz and applied load of 194 × $10^3$N, and several ballast settlement models.</figcaption></figure>

<figure><img src="image_4.png" /></figure>
<figure><img src="image_5.png" /><figcaption>**Table 1 | Cross-correlation coefficients between input and control parameters**</figcaption></figure>

<table><thead><tr><th></th><th>Γ</th><th>q<sub>d</sub>/ρ</th><th>τ<sub>∞</sub>/H<sub>0</sub></th><th>B</th><th>N<sub>0</sub></th><th>H<sub>0</sub></th></tr></thead><tbody><tr><td>Γ</td><td>1</td><td>0.03</td><td>0.91</td><td>-0.67</td><td>0.03</td><td>0.06</td></tr><tr><td>q<sub>d</sub>/ρ</td><td></td><td>1</td><td>0.02</td><td>-0.07</td><td>0.98</td><td>0.02</td></tr><tr><td>τ<sub>∞</sub>/H<sub>0</sub></td><td></td><td></td><td>1</td><td>-0.64</td><td>0.01</td><td>0.04</td></tr><tr><td>B</td><td></td><td></td><td></td><td>1</td><td>-0.07</td><td>-0.32</td></tr><tr><td>N<sub>0</sub></td><td></td><td></td><td></td><td></td><td>1</td><td>0.03</td></tr><tr><td>H<sub>0</sub></td><td></td><td></td><td></td><td></td><td></td><td>1</td></tr></tbody></table>

cycles with a nonlinear increase of settlement. This “intermediate” stage involves collective particle rearrangements since large voids are exhausted in the first stage. Hence, large-scale particle movements and longer excursions are necessary for increasingly smaller available void space¹,⁹. This slow relaxation is quite similar to that observed in glasses, spin glasses and flux lattices¹⁸⁻²⁰. The last stage of settlement is the long-term behaviour observed beyond 6000 cycles. In this part of the curve, the settlement has a slow linear evolution with the number of cycles. No collective particle rearrangements are observed so that the slow evolution should be attributed to “rare events” such as sudden micro-sliding induced by erosion or grain crushing. Beyond this schematic description, we are interested here in a functional form capable of predicting the settlement for a high number of cycles with parameters physically related to and measurable from packing fraction and loading conditions.

**Best-fit model.** We tested different models proposed for the railway ballast settlement based on loading parameters such as vertical stress, subgrade stiffness and calibration constants²¹. None of these models was found to describe satisfactorily our settlement curves shown in Fig. 2b. Shenton expression fits well the beginning of the settlement curve but overestimates long-term settlement. The expressions of Sato and Hettler do not provide a better fit until 5000 cycles. Moreover, a fundamental short-coming of such phenomenological models is that they introduce calibration constants that have no obvious physical interpretation and cannot be estimated by means of independent experimental measurements.

An alternative approach is the density relaxation law for granular materials tested on model systems. This law, called “Chicago fit”, was obtained by experimental analysis of the evolution of packing fraction for mono-disperse spherical particles (glass beads) in a cylindrical tube under a series of external excitations consisting of vertical shakes or “taps” applied to the container²,³,⁹. This relaxation law is an inverse-logarithmic law:

$$ \rho(t) = \rho_f - \frac{\rho_f - \rho_0}{1 + B \ln(1 + t/t_0)} \quad (3) $$

where $ρ_f$ is the ultimate value of packing fraction, $ρ_0$ is the initial packing fraction, B is a fitting parameter and $t_0$ is a characteristic time. Equation (1) given in Introduction is nothing but a transcription of equation (3) in which the number of cycles N replaces time, $τ_N$ replaces the packing fraction, and the characteristic time is replaced by a characteristic number of cycles $N_0$. From 360 independent settlement data, the parameters for the two models were identified from the loading parameters (amplitude and frequency) of each test.

Table 1 shows the inter-correlations coefficients among the model and loading parameters. The parameters $τ_∞$ and B are clearly linked to the dimensionless acceleration $Γ$ whereas $N_0$ depends on the initial packing fraction of the granular material. The correlation between $τ_∞$ normalized by the initial thickness of the granular layer $H_0$ and $Γ$ has a high correlation coefficient $R^2 = 0.9$ (Fig. 3a), thus suggesting a linear relationship between the two parameters. Different values of amplitude and confining pressure lead to a similar behaviour with the predicted values of $τ_∞/H_0$. In the same way, for different values of A, $ω$ and p but similar values of $Γ$, the settlement curves show a quite similar behaviour, as seen in Fig. 4a. The parameters A, $ω$ are independent from each other, while $Γ$ increase as a function of both (Fig. 4b and 4c). $Γ$ appears thus to be the governing parameter of settlement.

Figure 3b shows the correlation between B and $Γ$. Despite the relative dispersion of data due to the inherent variability of the system, it is seen that B globally declines as the acceleration parameter increases. Regarding the $N_0$ parameter, a logarithmic relationship with the tip strength (normalized by p) is observed; see Fig. 3c. Indeed, increasing the initial packing fraction, in correlation with the penetration strength $q_d$, leads to higher values of $N_0$ and therefore lower initial slope or slower settlement in the first loading cycles.

In Fig. 5a the experimental data are compared with the above model of settlement. This model describes quite well the evolution of settlement up to a high number of cycles and for the three stages of settlement when the three model parameters are adjusted. In order to assess the quality of this model, we calculate the mean square error (MSE) obtained between all the experimental and predicted values of settlement. MSE is defined as

$$ \text{MSE}(\hat{\theta}|0) = \langle (\hat{\theta} - 0)^2 \rangle \quad (4) $$

where $\hat{\theta}$ is the vector of predicted values for each settlement test and $θ$ is the vector of the corresponding measured values. In Fig. 5b we have plotted the cumulative distribution function of MSE between 360 experimental tests: 40% of the settlement curves are predicted with an error below 5% and in 60% of tests the error is below 10%. This indicates that the model based on the Chicago fit provides a quite accurate estimation of the evolution of settlement in coarse granular materials and, what is more, the model parameters have a clear

<figure><img src="image_6.png" /><figcaption>**Figure 3 | Correlations between the model parameters and loading parameters.** (a) The relation between settlement normalized by the thickness of the granular bed $τ_∞/H_0$ and $Γ$; (b) The relation between B and $Γ$; (c) The relation between $N_0$ and $q_d/p$.</figcaption></figure>
<figure><img src="image_7.png" /><figcaption>**Figure 4** | (a) Settlement for different values of A, $\omega$ and $p$, but the same value of $\Gamma$; (b) Values of $\Gamma$ as a function of the amplitude for $\omega$=18.85 Hz and $p$=138 × 10³ Pa; (c) Values of $\Gamma$ as a function of $\omega$, for A=0.8 × 10⁻³ m and $p$=170 × 10³ Pa.</figcaption></figure>

correlation with the packing fraction and loading parameters and can be obtained by independent tests.

## Discussion

In this work, we carried out a parametric study on a full-scale ballast track in order to test various models for the prediction of settlement in coarse granular materials under cyclic forcing. The settlement as a function of the number of cycles occurs in three stages with nearly the same functional dependence but varying with the initial state of the material and vibration amplitude. We showed that a model based on the Chicago density relaxation law provides an accurate estimation of settlement for a high number of cycles with parameters that reflect the loading parameters and initial mechanical state of the material and can be measured by independent tests. This model, together with its testing procedure, provides a reliable approach in applications where long-term settlement needs to be estimated in coarse granular materials.

Our findings raise also a fundamental question as to the predictability of the mechanical response of granular materials despite their natural variability and complex shapes of ballast grains. The behaviour is chaotic in the sense that small variations in the initial configuration of a granular packing are expected to be amplified with the number of loading cycles. This is not what we observe although fluctuations are observed between independent tests as shown in Fig. 3. This indicates that the fluctuating parts of consecutive

incremental settlements average out and the system tends to a well-defined ultimate mechanical state. In other words, the mechanical state of the granular material under vibrations is determined by its distance from the ultimate state. This is consistent with the form of equation (1) where the relevant settlement parameter appears to be $\tau_N/\tau_\infty$ with $\tau_\infty$ only depending on the dimensionless acceleration parameter $\Gamma$. This feature of compaction can be generic to all systems characterized by jamming transition and in which a maximally jammed state may be defined¹⁹,²⁰. Paradoxically, in such systems the long-term behaviour may be more accurately predicted than the short-term behaviour.

## Methods

The track model is composed of two monoblock sleepers lying on a ballast layer of thickness 0.35 m. The ballast layer is confined in a wood box of dimensions 4.11 m by 1.5 m. In this track model, the subgrade is represented by a wooden plate. The sleepers are fixed together by a rail section. A hydraulic cylinder applies the sinusoidal loading on the track. With the aim of providing a uniform load distribution, a beam is fixed to a hydraulic cylinder applying the overload on the rails. To study the influence of loading on settlement, we apply a sinusoidal load during 10⁴ cycles. The peak load is varied from 194 to 272 × 10³ N and the frequency from 3.3 Hz to 6 Hz. The cyclic loading is applied by a hydraulic cylinder on two sleepers placed on the granular bed. Before applying the cyclic loading, we characterized with several light penetration tests the initial mechanical state of the granular bed in our full-scale ballast track. In the course of cyclic loading, we measured with 4 non-contact optical sensors the evolution of settlement as a function of the number of cycles. The tests were carried out with 15 different combinations of the mean load and frequency and 6 independent initial conditions.

<figure><img src="image_8.png" /><figcaption>**Figure 5** | (a) Prediction of settlement by a logarithmic law in comparison to the experimental data for two tests (frequency of 4.5 Hz and applied load of 239 × 10³ N; Frequency of 3.3 Hz and applied load of 194 × 10³ N). The inset shows the data on a log-linear scale; (b) Cumulative distribution function of the mean square error (MSE) of model prediction compared to 360 experimental tests.</figcaption></figure>
<figure><img src="image_9.png" /></figure>

1. Jaeger, H. M., Nagel, S. R. & Behringer, R. P. Granular solids, liquids, and gases. *Rev. Mod. Phys.* **68**, 1259–1273 (1996).

2. Knight, J. B., Fandrich, C. G., Lau, C. N., Jaeger, H. M. & Nagel, S. R. Density relaxation in a vibrated granular material. *Phys. Rev. E* **51**, 3957–3962 (1995).

3. Nowak, E. R., Knight, J. B., Ben-Naim, E., Jaeger, H. M. & Nagel, S. R. Density fluctuations in vibrated granular materials. *Phys. Rev. E* **57**, 1971–1982 (1998).

4. Nicolas, M., Duru, P. & Pouliquen, O. Compaction of a granular material under cyclic shear. *Eur. Phys. J. E. Soft Matter* **3**, 309–314 (2000).

5. Philippe, P. & Bideau, D. Numerical model for granular compaction under vertical tapping. *Phys. Rev. E. Stat. Nonlin. Soft Matter Phys.* **63**, 051304 (2001).

6. Philippe, P. & Bideau, D. Compaction dynamics of a granular medium under vertical tapping. *Europhys. Lett.* **60**, 677 (2002).

7. Richard, P., Nicodemi, M., Delannay, R., Ribière, P. & Bideau, D. Slow relaxation and compaction of granular systems. *Nat Mater.* **4**, 121–128 (2005).

8. Vandewalle, N., Lumay, G., Gerasimov, O. & Ludewig, F. The influence of grain shape, friction and cohesion on granular compaction dynamics. *Eur. Phys. J. E. Soft Matter.* **22**, 241–248 (2007).

9. Ben-Naim, E., Knight, J. B. & Nowak, E. R. Slow relaxation in granular compaction. *J. Chem. Phys.* **100**, 6778 (1996).

10. Ribière, P., Richard, P., Philippe, P., Bideau, D. & Delannay, R. On the existence of stationary states during granular compaction. *Eur. Phys. J. E. Soft Matter.* **22**, 249–253 (2007).

11. Cerasi, P. & Mills, P. Insights in erosion instabilities in nonconsolidated porous media. *Phys. Rev. E* **58**, 6051–6060 (1998).

12. Foster, M., Fell, R. & Spannagle, M. The statistics of embankment dam failures and accidents. *Can. Geotech. J.* **37**, 1000–1024 (2000).

13. Lobo-Guerrero, S. & Vallejo, L. E. Discrete element method analysis of railtrack ballast degradation during cyclic loading. *Granul. Matter.* **8**, 195–204 (2006).

14. Tutumluer, E., Huang, H., Hashash, Y. & Ghaboussi, J. Aggregate shape effects on ballast tamping and railroad track lateral stability. Paper presented at the AREMA annual Conference. Louisville. Arema (2006).

15. Gourvès, R. & Barjot, R. The panda ultralight dynamic penetrometer. Paper presented at the Eleventh European Conference on Soil Mechanics and Foundation Engineering. Copenhagen. Dansk Geoteknisk Forening. (1995).

16. Chaigneau, L. Caractérisation des sols de surface à l'aide d'un pénétromètre. Ph.D. thesis, Université Blaise Pascal, Clermont-Ferrand (2001).

17. Breul, P., Benz, M., Gourvès, R. & Saussine, G. Penetration test modeling in a coarse granular medium. Paper presented at Powders and grains 2009:

Proceedings of the 6th international conference on micromechanics of granular media. Golden, Colorado. AIP Conference Proceedings. (2009).

18. Nicodemi, M., Coniglio, A. & Herrmann, H. J. Frustration and slow dynamics of granular packings. *Phys. Rev. E* **55**, 3962–3969 (1997).

19. Liu, A. J. & Nagel, S. R. Nonlinear dynamics: Jamming is not just cool any more. *Nature.* **396**, 21–22 (1998).

20. Liu, A. J. & Nagel, S. R. Jamming and rheology: constrained dynamics on microscopic and macroscopic scales (CRC Press, 2001).

21. Dahlberg, T. Some railroad settlement models - a critical review. *Proc. Inst. Mech. Eng. F J. Rail Rapid Transit.* **215**, 289–300 (2001).

## Acknowledgments

The authors would like to thank Rodolphe Potvin and the personal staff from the CEV test centre of the SNCF.

## Author contributions

J.C.Q. performed the experiments and prepared the figures. J.C.Q., G.S., P.B. and F.R. worked on data analysis, did the theoretical work and reviewed the manuscript.

## Additional information

Competing financial interests: The authors declare no competing financial interests.

How to cite this article: Quezada, J.C., Saussine, G., Breul, P. & Radjaï, F. Predicting the settlement of coarse granular materials under vertical loading. *Sci. Rep.* **4**, 5707; DOI:10.1038/srep05707 (2014).

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

This work is licensed under a Creative Commons Attribution 4.0 International License. The images or other third party material in this article are included in the article's Creative Commons license, unless indicated otherwise in the credit line; if the material is not included under the Creative Commons license, users will need to obtain permission from the license holder in order to reproduce the material. To view a copy of this license, visit http://creativecommons.org/licenses/by/4.0/