A geomechanics classification for the rating of railroad subgrade performance

The type of subgrade of a railroad foundation is vital to the overall performance of the track structure. With the train speed and tonnage increase, as well as environmental changes, the evaluation and influence of subgrade are even more paramount in the railroad track structure performance. A geomechanics classification for subgrade is proposed coupling the stiffness (resilient modulus) and permanent deformation behaviour evaluated by means of repeated triaxial loading tests. This classification covers from fine- to coarse-grained soils, grouped by UIC and ASTM. For this achievement, we first summarize the main models for estimating resilient modulus and permanent deformation, including the evaluation of their robustness and their sensitivity to mechanical and environmental parameters. This is followed by the procedure required to arrive at the geomechanical classification and rating, as well as a discussion of the influence of environmental factors. This work is the first attempt to obtain a new geomechanical classification that can be a useful tool in the evaluation and modelling of the foundation of railway structures.

Abstract The type of subgrade of a railroad foundation is vital to the overall performance of the track structure. With the train speed and tonnage increase, as well as environmental changes, the evaluation and influence of subgrade are even more paramount in the railroad track structure performance. A geomechanics classification for subgrade is proposed coupling the stiffness (resilient modulus) and permanent deformation behaviour evaluated by means of repeated triaxial loading tests. This classification covers from fine-to coarse-grained soils, grouped by UIC and ASTM. For this achievement, we first summarize the main models for estimating resilient modulus and permanent deformation, including the evaluation of their robustness and their sensitivity to mechanical and environmental parameters. This is followed by the procedure required to arrive at the geomechanical classification and rating, as well as a discussion of the influence of environmental factors. This work is the first attempt to obtain a new geomechanical classification that can be a useful tool in the evaluation and modelling of the foundation of railway structures. Diameter corresponding to n% of passed material in a grain size distribution curve CBR California bearing ratio

Introduction
The understanding and knowledge about the deformation and failure mechanisms of the geomaterials under repeated loadings are extremely important for the proper design and maintenance planning in the railway structures [1][2][3][4]. Typically, the subgrade layer of railway track is composed of geomaterials (mainly coarse-grained soils due to its better performance when compared to the fine-grained soils) and presents two types of deformation when submitted to cyclic loads: the recoverable (quasi-elastic, resilient or reversible) and permanent (or plastic) deformations, which have significant importance in the long-term performance of the substructure [5][6][7]. Regarding the resilient deformation, the concepts were introduced by AASHTO in 1986 [8]. This property is often used to characterize the materials in several railway layers [9][10][11]. About the permanent deformation, an accurate estimation of the amount of cumulative settlement will benefit railway structures to avoid a mediocre performance and consequently higher annual maintenance costs [12][13][14]. Overall, resilient and permanent deformations are both important parameters for the design of railway structures [15][16][17][18]. However, the complexity and time-consuming of laboratory tests, mainly for permanent deformation evaluation, conduct to the development of predictive models based on the available data from the existing laboratory and/or field results.
In this line, one of the main objectives of this work is to review the available resilient and permanent deformation models in the scope of the geomaterials published in the specialized bibliography considering the information about the conditions for its development (properties of the materials tested, compaction tests, degree of compaction, moisture content, etc.) and the main variables/factors that can influence the response of the material. This means that this analysis also includes external factors as the physical state of the soil, which is often difficult to control because it depends on other environmental aspects. A first attempt was already performed for permanent deformation models in the work developed by [19]. Consequently, this work intends to complement it with the resilient deformation models, including comparisons among some of the selected models that are based on different materials with different classifications (UIC and ASTM), properties, granulometry and physical states. By coupling both resilient-and permanent-based deformations, we rate the geomaterial classifications according to the predicted resilient and permanent deformation performances, which is believed to be a helpful tool in the design of the railway structures. It is also noted that this rating should be interpreted under normalized conditions because it depends on several properties and soil state conditions. Based on these results, a novel geomechanical classification for track subgrade/ formation that couples the resilient modulus and the permanent deformation is purposed using a similar approach developed by [20].

Influence of subgrade in design and performance of railroad track structures
The quality and the support given by the subgrade are important aspects in the performance of the railroad track structure. The subgrade is an integrant component of the track structures (ballasted and ballastless tracks), and its properties are important in the performance of track and track quality. Indeed, this is why subgrade performance indicators are important and should be implemented into the track designs and assessments. Thus, there are numerous means of measuring and predicting subgrade performance, which include the laboratory and in situ tests. The in situ tests can be advantageous since the results are representative of the real conditions and some can be obtained quickly and more cost-effective. However, laboratory tests could be also very interesting and necessary to carry out parametric studies considering variations occurring during the service life of the structure. Indeed, recently, new methods have been developed to determine the resilient modulus, such as the performance of the cyclic lightweight deflectometer test [21]. Thus, despite the importance of the subgrade and its significant impact on the track quality and required maintenance, the evaluation of the subgrade is not straightforward since there is an extensive list of factors that can affect its short-and long-term performance: the moisture content and/ or, suction, and/ or degree of saturation, temperature (freezing and thawing), shear strength, stiffness parameters, consolidation, etc. Indeed, recent studies show the influence of soil suction on the deformation characteristics of railway formation materials [22]. Furthermore, other variables can affect the performance of the subgrade related to the train type (passenger or freight), axle load, train speed, train configuration, drainage, rail condition, tie spacing, and wheel condition [23][24][25]. Thus, these factors added to the environmental factors can affect the performance of the subgrade and also the global performance of the track structure. Indeed, the poor performance of the subgrade can lead to some issues in the railroad track structure: ballast fouling, ballast pockets, pumping of soil fines through the ballast, and slope stability failure. Thus, there are two main types of subgrade failure (mostly related to the fine subgrade soils): subgrade progressive shear failure and excessive subgrade plastic deformation which are well documented [15,17,26]. More recent works in Refs. [27] and [28] emphasize that lowand medium-plasticity soils can also be subjected to different types of failure such as fluidization and mud pumping due to excessive increase in axle loads of freight trains in recent years. Nevertheless, these phenomena and many others like liquefaction (large displacement caused by repeated loading, saturated silt, and fine sand), massive shear failure (slope stability-this failure is associated with the weight of the train, inadequate soil strain and increase of water content), consolidation settlement (the static soil stress increased due to the embankment weight and saturated fine-grained soils), frost action (the track become rough caused by periodic freezing), swelling/shrinkage, slope erosion and soil collapse caused by ground settlement are out of scope of this paper.
The design of the subgrade (design period corresponds to 50-100 years) is defined according to the type of traffic, bearing capacity of the subgrade, configuration of the track, climatic and hydro-geological conditions [15,29,30]. Therefore, the subgrade should be designed to show a good short-and long-term performance, resisting to failure and excessive deformation, respectively, induced by repeated loads. The short performance is influenced directly by the resilient modulus (that has influence on the stress levels on the subgrade) and by the strength parameters, among other factors. The stress paths and the strength parameters (such as cohesion and friction angle) have a significant effect on the development of the permanent deformation. (This topic will be analysed in more detail further in this work.) Usually, the design of the ballasted track involves the determination of the deformation and stresses at critical locations in the track (sleeper-ballast or ballast-sub-ballast contacts). Posteriorly, the magnitude of the settlements and stresses (induced by the passage of the trains) is compared with allowable values in an iterative process where the dimensions of the sleepers and/or thickness of the granular layers are adjusted. These methodologies are applied by several authors and standards: • Reference [15] developed a methodology that prevents progressive shear failure and also excessive plastic deformation where the plastic strain (e p ) and the cumulative deformation should be limited to 2% and 25 mm, respectively. • Reference [31] presented a set of references for the design and maintenance of track substructure. According to [31], there are four soil quality classes: QS0 (unsuitable soils), QS1 (poor soils), QS2 (average soils) and QS3 (good soils). This concept will be explored further in this work. From the soils' classification and thickness of the granular layers, the track is evaluated in three categories: P1, P2 and P3 (mediocre, average and good, respectively). • British Rail Method and Network Rail Code defined recommendations related to the thickness of the trackbed layers [17]. • West Japan Railways defined a limit value for bearing capacity (288 kPa). Below this value, the improvement of the foundation is required. Thus, the bearing capacity of 288 kPa equates to a compressive strength r s of approximately 112 kPa assuming a cohesion model plastic solution to a simple strip footing where r b = 2.57r s .
These methods allow preventing the development of an overstressed subgrade that can present cumulative permanent deformation (formation of ballast pocket) or progressive shear failure (subgrade squeezing) under repeated loads. Indeed, these are common failure modes of typical subgrades composed of fine-grained cohesive soils [15,26,30,32,33]. These types of problems can lead to an increase in dynamic loads and accelerate the deterioration of the track.
Regarding the ballastless track, the design theories of the ballastless track vary across the world according to the experience and background of each country [34].
Thus, two important parameters should be included in the design of railway tracks related to the short-and longterm performance: the resilient modulus (M r ) and the permanent deformation (e p ). These two parameters should be analysed together in an integrated approach since, from a simplistic point of view, the stress levels at the rail track subgrade vary with the value of its M r , which have impact on the permanent deformation. These parameters will be analysed together with the UIC classification in this work to purpose a novel geomechanical classification.

Resilient modulus
The resilient modulus is used to characterize the recoverable, reversible or quasi-elastic behaviour. In the traditional theory of elasticity, the elastic properties (which are material parameters) are defined, mostly, by the Poisson's ratio (m) and Young's modulus (E). In the analysis of the response of the material when submitted to cyclic loads, the modulus of the elasticity is replaced by the resilient modulus to consider the nonlinearity (i.e. dependence on stress level) on the performance of the material, as well as the inelastic properties of the materials, i.e. the loading and unloading of the stress-strain curve are not totally overlapped due to the dissipation of the energy, as depicted in Fig. 1. Indeed, in the case of repeated load tests with constant confining pressure, the resilient modulus and the Poisson's ratio are defined, respectively, by the following expressions: where r d is the deviatoric stress, e 1;r is the recoverable axial deformation and e 3;r is the recoverable horizontal deformation.
The test method used to determine the experimental value of the resilient modulus is described in [8] and is based on the application of deviatoric stress under constant confining stress. The response of the materials varies according to its own nature. The cohesive materials are more sensitive to the deviator stress than the confining pressure since the experimental results show that there is a decrease of the resilient modulus with the increase in the deviatoric stress. On the other hand, the granular materials show a different tendency since the resilient modulus increases with confining stress. The more recent models are dependent on important factors regarding environmental aspects such as moisture content, suction and degree of saturation [9,[35][36][37][38][39][40]. Additionally, recent work by [41] also proposes a methodology to obtain the damping from these type of tests using the standard AASHTOT-307-93 [42].
The computational modelling of resilient behaviour is a very difficult task, due to the necessity of decoding complex behaviours into simple mathematical expressions and procedures for routine analysis.

Mechanistic-empirical modelling approach
The empirical approaches can be used in the design procedures of railway structures. This means that a relationship is established between the design inputs (materials, loads, environment, geometry) and, for example, rail track failures through experience, experimentation or a combination of both. The more complex empirical approaches are based on empirical equations derived from experimentations. The relationship between rail track failure and the physical phenomena is described by empirical equations.
The mechanistic-empirical modelling approaches used to determine the resilient modulus are divided according to the type of material tested (granular or cohesive). Furthermore, the models more complex depending on the physical state of the material are also presented. Moreover, in appendix, the state of the art of resilient models developed based on laboratory tests (mostly cyclic triaxial tests) considering all types of materials (such as clays, silts, sands and gravels) is summarized.

Cohesive materials
In the cohesive materials, the fines content and its mineralogy are important factors in the analysis, which means that the resilient modulus should be dependent on the soil physical state. Indeed, in cohesive materials, certain characteristics may have more importance in the value of the resilient modulus, besides the stress state, such as the moisture content, suction and saturation degree. In the work developed by [43], the authors state that the new generation models not only provide a better fit than the older models, but they also provide a reasonable fit to the data that can capture the effects of stress state, soil type, soil structure and the soil physical state quite effectively. Furthermore, according to [44], there are prediction models that combine the effect of the stress state and matric suction on the resilient modulus. These models not only reflect the relationship between the resilient modulus of subgrade soils and the stresses but also the effects of seasonal variation of moisture content on the resilient modulus. The authors also conclude that the degrees of stress and moisture content have a significant impact on the resilient modulus of compacted cohesive soils. Moreover, the results show that the resilient modulus of the tested cohesive soils increases with the increase in effective confining pressure, matric suction, and degree of compaction and decreases with the increase in deviator stress and moisture content. Therefore, lower moisture content concomitantly with a higher degree of compaction is beneficial to the stiffness of the subgrade, mainly in those not sensitive to water content changes.
Considering the extensive bibliography about this topic, several models were developed to describe the resilient modulus and are summarized in Table 1.
Furthermore, it is also important to take into account the type of approach in terms of total and/or effective stresses. In the recent work developed by [52], the developed model can describe the phenomena of modulus development during the cyclic undrained condition and takes into account the actual values of effective stress (p'), excess pore water pressure, the loading characteristic and the position of the effective stress path to the failure line.

Granular materials
The modelling of the resilient behaviour in the case of granular materials can be performed through two different approaches: simulation of the resilient modulus considering a constant value of the Poisson's ratio and by separation of the deformation response of the material into the volumetric and shear components, which is more complex [53].
The models expressed by the resilient modulus and constant Poisson's ratio are very easy to understand, very simple and easy to implement numerically. In fact, initially, in terms of variables, the models were only dependent on the confining stress and material parameters and/or the sum of the principal stresses [54,55]. Due to the necessity to include the influence of the deviator stress, [56] modified the well-known k-h model, introducing the influence of the mean and deviatoric stress (p and q, respectively): where h is the sum of principal stresses, and k 1 , k 2 and k 3 are parameters dependent on the state of the soil and its characteristics. Posteriorly, some models were improved considering, for example, the porosity of the material and variable confining stress [57,58].
However, some of the previous models show drawbacks concerning the assumption of a constant value of the Poisson's ratio: despite the good results regarding the axial deformations, it was difficult to simulate correctly the volumetric and radial deformations. Indeed, some authors also suggested that the use of these models could lead to Poisson's ratios superiors to 0.5. However, some of these models are still used and were also improved because of their simplicity [59].
Nevertheless, in order to solve Poisson's ratio problem, other models were developed, expressing the resilient behaviour through the shear modulus (G) and bulk modulus (K). This formulation (in terms of G and K) is more comprehensive than the models only characterized by the resilient modulus and Poisson's ratio since it is dependent on more complex laboratory tests and with a more generic interpretation than a cyclic axial loading.
Indeed, Ref. [48] cited by [60] identifies three important conditions in the application and formulation of these models: • In each increment of calculus, a linear elastic behaviour is adopted. • The shear and volumetric components of stress and strain are analysed independently. • Better adjust to the experimental results, mostly when the solicitation presents a 3D character.
Reference [61] developed a well-known model dependent on bulk modulus and shear modulus: where e v and e q are respectively the volumetric and shear resilient strains: K 1 , b and G 1 are parameters of the model, with b = (1n)K 1 /6G 1 .
In this model, there is a clear division between the volumetric and shear deformation and the values of K and G are stress-dependent (p and q, respectively). Posteriorly, this model was updated by other authors [53,[62][63][64] in Table 1 Resilient moduli for cohesive materials (adapted from [45])

Models
Expressions Parameters Authors K 1 , K 2 , K 3 and K 4 Thompson and Robnett [46] Power M r ¼ kr n d k and n are parameters dependent on the type of soil and physical state; the parameter n usually has a negative value Moossazadeh and Witczak [47] Power (with the influence of the confining stress for overconsolidated saturated soils) A geomechanics classification for the rating of railroad subgrade performance 327 order to consider the anisotropy of the material, which is an important characteristic regarding the deposition and compaction conditions. Indeed, [63] and [53] purpose and apply an orthotropic version of Boyce's model introducing an anisotropic coefficient (c).

Models dependent on the physical state-suction and degree of saturation
Recently, some studies were developed in order to include the physical state into the resilient models [9,35,36,65].
In this section, a brief reference is made regarding the models dependent on the suction and degree of saturation. These models were developed based on experimental results and include suction as one of the main parameters. In the case of non-saturated soils, three important variables influence recoverable/reversible behaviour [66]: • net confining stress: r 3 -u a ; • deviatoric stress: where r 3 is the confining stress, r 1 is the axial stress, u a is the air pressure and u w is the water pressure.
Based on Uzan's model, Ref. [67] developed a model that includes the suction in the determination of the resilient modulus: where p a is atmospheric pressure; s oct is octahedral shear stress; h net is net bulk stress, h net = h-3u a , and h net is bulk stress; Dw m is the variation of the suction matrix regarding the initial suction (Du w-sat = 0); and Du w-sat represents the increase of the pore pressure in saturated conditions (w m = 0). Despite the simplicity of the model formulation (which is based on the work developed by [56]), the characterization of the suction increases the complexity and application of these types of models since this is a parameter that is difficult to characterize.
Recent studies show the influence of the water content (with special attention on the role of the suction in the interpretation of the results) and consider the possibility to apply an effective stress approach to characterize the changes in normalized resilient modulus as a function of a single parameter that includes the total stresses and suction [37]. Moreover, some studies are focused on the importance of the degree of saturation (and the controlling in the soil compaction) and its relationship with the soil structure design. In fact, the degree of saturation and also the CBR (California bearing ratio) are parameters that are easy to obtain and easy to use when compared to the suction. The authors conclude that CBR (soaked and unsoaked) and also the elastic shear modulus (representative of the stressstrain behaviour at small strains of the subgrade), unconfined compression strength, cyclic undrained shear strength, among other parameters, of unsaturated soil are controlled by q d and S r at the end of compaction [39]. More specific studies about compaction conditions and percentage of fines can be found in [68] and [69].

Main differences of the resilient modulus formulation-analysis and comparison
As previously described, the resilient behaviour is dependent on several factors such as the stress level, fines contents, particle shape, grading, density, maximum grain size, aggregate type, moisture content, stress history and number of load cycles. From all of them, two important factors stand out: the applied stress level and moisture content. The stress (in terms of confining stress, sum of principal stresses and deviator stress) is the most important factor in the analysis of the resilient behaviour since its impact on the response of the material is very significant. Indeed, most of the geomaterials (especially the granular materials) show a stress dependency under repeated loads. To replicate this nonlinearity and time dependency, the resilient response of the materials is not solved by the traditional elasticity theory. However, as mentioned previously, there are other formulations to characterize the response of the material such as the replacement of the resilient modulus and Poisson's ratio by the shear and bulk modulus. This formulation is simple and based on curve-fitting. Indeed, this approach deals better with the nonlinear response of the geomaterials (mostly in the case of granular materials) from a theoretical point of view, and it is more realistic regarding the physical meaning in a 3D stress regime. However, as mentioned previously, these models have a more complex formulation and the parameterization is also difficult to obtain through the available test data [70].
Another relevant aspect is related to the fact that some of the resilient models only fit the laboratory data used for their own development, which means that the mathematical formulation can only be used in a particular situation.
The models that include the influence of the moisture content and/or degree of saturation describe better the material conditions. In the case of the granular materials, the granulometry and the fines content (in a compacted well-graded material) have a significant influence on the response of the material since the water cannot drain freely since it is kept in the pores of the materials [71]. Indeed, when the saturation is close to the maximum, the resilient behaviour is affected [72]. Nevertheless, these models are more complex since it is difficult to characterize these parameters related to the suction and degree of saturation.
More recently, some resilient models have been developed using machine learning and statistical techniques and optimization tools as described in the work developed by [73] and [74].
Despite the several advantages of the empirical models, the main limitation is related to the confidence in the extrapolation of the analysis beyond the conditions in which they were defined. The tables presented in appendix (Tables 17-21) contain information regarding the type of soils (UIC and ASTM classification), source (authors), the mathematical models and respective variables and empirical constants. This information is quite significant in the modelling of the subgrade of the railway tracks since supports information about the materials used to define the model, as well as its physical conditions such as moisture content and dry density.

Permanent deformation
The permanent deformation is induced by repeated traffic loading and may occur in railway structures. The development of permanent deformation is the result of the accumulation of deformation throughout millions of cycles, which is a complex process and depends on several factors such as the number of load cycles (N), stress levels, loading history, the effect of the principal stress rotation, moisture content, density, aggregate type and particle size distribution. This development of permanent deformation during N loading cycles can stabilize or, in the worst situation, may lead to the ultimate collapse of the structure (excessive rutting). Regarding the moisture content, despite the extensive research about the influence of this parameter on the reversible deformation, a recent study was developed to study the effects of the compaction moisture content on the plastic behaviour of fine soils [75].
To predict permanent deformation, several approaches were proposed: numerical simulations using elastoplastic models, shakedown theory or mechanistic-empirical permanent deformation models based on laboratory tests such as cyclic triaxial tests or the hollow cylinder tests [19]. Throughout the analyses presented by these authors, more emphasis is given to the mechanistic-empirical permanent deformation models from the laboratory testing until the modelling. In fact, the work done in this paper for the resilient modulus is a mirror of the one for permanent deformation.
5 Geomechanics classification and the rating

Resilient modulus: parametric study
Considering the extended information presented in Sect. 3 and also the tables presented in appendix, a first attempt was done to compare the resilient modulus of different materials under a certain stress level. This analysis is important to comprehend how the models differ from each other in terms of formulation and the importance of certain variables.
The models were calibrated considering the same reference stress path for all materials during the calibration process. The stress path described by [76] was used. This stress path is characterized by a cyclic deviator stress of 24 kPa and a constant confining stress of 60 kPa. During the cyclic tests, the stress ratio (r d /r c ) was kept constant at 0.4, which is a representative ratio in the subgrade of a rail track full-scale model test. It is noted that other confining stresses from 60 to 210 kPa were also tested.
In this analysis, four models were applied in order to compare the resilient modulus for the coarse-grained soils, as depicted in Table 2. In order to facilitate the interpretation of Table 2, C and m are material constants and r 0 v is the effective vertical stress (more details in Tables 19 and  20)).
The first part of this work is based on the experimental results presented in [80]. The authors studied six granular materials representative of the base and sub-base materials used on flexible pavements. [80] performed rapid shear tests and repeated load tests in order to determine the shear strength parameters (cohesion and friction angle), resilient modulus, rutting potential and moisture susceptibility. The properties of the materials are presented in Table 3.
The constants of each model were determined by regression analyses. The values of these constants reflect the material properties as well as the loading cycle applied in each test that may vary from author to author. For example, in the work developed by [80], the specimens were conditioned for 1000 load repetitions while in the work developed by [9], the authors follow the procedures described in AASHTO T307-99. However, for practical purposes, the conditioning adopted for resilient modulus determination following the standard procedures (AASHTO T307-99, [40]) intends to assure that the number of cycles during resilient modulus determination will not affect the results. The results of the four models and the two materials are presented in Table 4. It is important to refer that for material 2, due to the lack of data, it was not possible to present the regression results for the model developed by [79] and MEPDG model. Analysing Table 4, in general, the models present similar r-squared values. Therefore, considering the obtained constants and the stress path selected, the models were compared in terms of the resilient modulus, as depicted in Table 5. It can be observed that the models present similar results. Material 2 presents an inferior resilient modulus when compared to the well-graded gravel, as expected. Regarding the material classified as QS3-GW, it is evident the close values between the Uzan's model, k-h and MEPDG model. After the analysis of the coarse-grained soils, the finegrained soils were evaluated. In the case of the fine-grained soils, four models (Table 6) and two materials (Table 7) were compared. The materials were tested through the triaxial test and both materials were compacted at optimum water content and 2% above the optimum. The measured resilient modulus is based on the work developed by [9].
The target moisture content and densities were selected based on AASHTO T99 test results and field-measured values [80] Well-graded gravel QS3-GW Considering the results obtained by [9] regarding the measured resilient modulus, the constants of the parameters were determined through regression analyses (Table 8). However, it is important to refer that the MEPDG model is defined in terms of total stresses, which means that the regression analysis is performed for each material on both moisture contents (opt and opt ? 2%). The remaining models are dependent on effective stresses, which means that the regression analysis is performed considering the whole set of moisture content (opt and opt ? 2%). This type of information allows analysing the sensibility of the resilient modulus regarding the moisture content. The models developed by [9,35] and [81] are dependent on the w m and v m parameters, which means that the results are presented considering M opt r and M optþ2% r , although the regression is performed with the whole set (opt and opt ? 2%).
The model developed by [81] is simpler than the remaining models, and the influence of the moisture content is indirectly included in the suction matrix in the initial mean normal stress (p 0 0 ). From multiple linear regression analyses, the empirical constants were determined, and the values are presented in Table 8.
From the constant parameters and the selected stress path, the resilient moduli of the materials were found, as depicted in Table 9.
Considering the obtained results for the fine-and coarsegrained materials, the MEPDG model [59] will be adopted in the following analysis, despite its complexity when compared to the remaining models. The main reason is related to its universal character since the model can be  Table 6 Models selected for the fine-grained soils Refs. Models The triaxial tests following the AASHTO T307-99 procedure were conducted on both materials   applied in the prediction of the resilient modulus of all types of geomaterials. This analysis shows a possible way to rank materials according to the resilient modulus, as depicted in Fig. 2. In the case of the material classified as QS2/QS3-SW/SM-ML, instead of the MEPDG model, the Uzan's model was used due to the lack of information. (It was not possible to perform a regression analysis.) However, as the previous analysis shows, the MEPDG and Uzan's models present similar results. All the materials present in Fig. 2 are at optimum moisture content, except the material classified as QS2/QS3 SW/SM-ML. In this case, the moisture content is close to the optimum conditions.

Permanent deformation: parametric study
The comparison regarding the permanent deformation models was already presented in the work developed by [19]. As in the previous exercise, in this parametric study, the authors attempted to compare different permanent deformation models available for different types of soils considering the empirical permanent deformation model developed by [76]. Thus, the results showed that two soils can be integrated into the same classification (UIC or ASTM), but the laboratory conditions may differ greatly and therefore lead to different results in terms of permanent deformation. The calibration was performed to find the best-fit for the experimental data through the parameters e p;o , a and B (which correspond to the model constants associated with the material properties). The models were calibrated considering the same stress paths (described by [76]). As in the resilient modulus analysis, the materials selected are representative of different types of materials.
In order to understand better the ranking, the materials selected for the preliminary analysis of the permanent deformation are depicted in Table 10. This table also includes information about the ASTM and UIC classifications. To complete the information about the materials, the physical properties, state conditions and strength properties of the materials are depicted in Tables 11, 12 and 13, which would be extremely important to understand the novel geomechanical classification purposed.
In the work developed by [19], the ranking process is described. The final ranking is depicted in Fig. 3.

Geomechanical classification
Considering the results presented in the previous sections, two models were selected to characterize the resilient modulus and permanent deformation of the geomaterials. Indeed, this is an attempt to relate the resilient modulus with the permanent deformation under a certain stress level. As in the previous analysis, the selected stress path is characterized by a cyclic deviator stress equal to 24 kPa and a confining constant stress equal to 60 kPa. The permanent deformation was determined considering a number of load cycles equal to 30,000. The selected materials are   Table 14. Some of the materials are the same used in Sects. 5.1 and 5.2, and others were added to increase the robustness of the analysis, the ranking and the impact of the final result of the geomechanical classification. These ''auxiliary'' data only have information about the resilient modulus or permanent deformation (they were used to define the ''limits'' of the classification). The work developed by [84] was also included but due to the significant information about both types of deformation (resilient and permanent) under different water content conditions. It is important to refer that, for these particular materials, are presented, whenever it is possible, the values of both deformations: resilient and permanent deformation.
This table also includes information about compaction conditions. The resilient modulus and permanent deformation of the selected materials (Table 14) are presented in Table 15. The materials with both values of M r and e p allow validating the mechanical classification.
From the results presented in Table 15, it is possible to define guiding/limit values associated with the UIC classification: QS1, QS2 and QS3 [31]. This work only presents a preliminary geomechanical classification. The lack of information about the models' parameters found in the bibliography (namely the models that include the suction and other complex model coefficients) and the necessity and difficulty to performed back analysis to find important  ** s and m are the parameters used to define the failure envelope: q = s ? mp parameters (as the cohesion and friction angle used to define the failure criterion) reduced the number of models and materials that could be used in this geomechanical classification. Besides, some of the models found in the bibliography tested materials classified as QS0, which implies its elimination in the geomechanical classification, since these materials require mechanical improvement. The geomechanical classification is an attempt to define a novel helpful guide to be used as a support for the modelling and design of the substructure. Indeed, this novel geomechanical classification tries to rank soils in a very similar way developed by [20] that was inspired by [40]. This work allows understating which are the acceptable values of the permanent deformation and resilient modulus according to the type of material and its classification. In this particular case, [31] classification was used as depicted in Fig. 4.
Analysing Fig. 4, a significant range of values was obtained for each classification. This was expected due to the mineralogical nature and physical properties of each material, which can have an influence on several parameters used in the determination of the resilient modulus and permanent deformation. Indeed, in the case of the materials classified as QS1 (fine soils), the type of the material (CL, CH, ML, or MH) associated with its plasticity's characteristics and consistency index (as well as the percentage of fines) can have a significant impact in the permanent deformation and resilient behaviour, as depicted in Fig. 5. For the materials classified as QS2 and QS3, the type of materials (sand or gravel) as well as its granulometry (wellor poor-graded) and the percentage of fines (in the wellgraded materials) can also influence the response of the material in terms of resilient and permanent deformations.
With this work, it was possible to define subsets associated with the properties of the soils, despite the need for more information. For example, as Fig. 5 shows, a silty soil (classified as QS1-ML) presents a superior resilient modulus when compared to clayed soils (QS1-CL). Regarding the granular materials, the subsets can be defined by the granulometry of the materials (well-or poor-graded). However, it is important to refer that the percentage of fines can also influence the obtained results. Indeed, despite the importance of the granulometry and plasticity properties, the fines content and the moisture content have also significant importance, mostly in the fine soils. As mentioned previously, the work developed by [84] shows the influence of the moisture content in the soils. Thus, [84] selected and studied the behaviour of the clay considering three compaction moisture content and dry unit weight conditions representing dry of optimum moisture content level, optimum moisture content level (corresponds to 95% of maximum dry unit weight), and wet of optimum moisture content level. From the results obtained by [84], an estimation was performed (in terms of resilient modulus and permanent deformation), and the results are presented in Figs. 6 and 7. These results are based on the properties of the clay described in Table 16. Analysing Fig. 7, the moisture content can influence both deformations, which means that can affect severely the track performance. In a more detailed analysis, it is possible to identify a reduction of around 200% in the resilient modulus from dry to wet conditions. This difference is significant, which means that the moisture content is an important parameter since influences the performance of the subgrade and corroborates the findings of [20,86] in the case of non-standard unbound granular materials for pavements. A geomechanics classification for the rating of railroad subgrade performance 335 The triaxial tests following the AASHTO T307-99 procedure were conducted on both materials

Conclusions
The capability to determine the resilient modulus and also the permanent deformation of geomaterials is imperative in the modelling, design and evaluation of the performance of railway structures. This work attempts to identify and summarize the main models to estimate the resilient modulus and also permanent deformation. At the beginning of the analysis, the most simplistic models are presented. In the case of the resilient modulus, the models should include the influence of the deviatoric and confining stresses and the constant parameters should reflect the influence of other factors as the moisture content and physical state of the material. The mechanistic-empirical models summarized in Annex were divided according to the type of material (clay, silt, sand and gravel) to better understand the model's formulation and the conditions of the material when tested in terms of its granulometry (C u and C c ) and plasticity properties (W L , W P , and IP). Tables also include some observations regarding the laboratory tests and classification of the material (UIC and ASTM). From the available mechanistic-empirical approach, it was important to evaluate the robustness of the models and their sensitivity. Thus, some models and materials were selected to perform a parametric study in terms of resilient modulus. A similar analysis was already performed for the permanent deformation in the work developed by [19]. The selection depended on the data available for each material and the variability of the geomaterials (type of soils) in terms of granulometry, percentage of fines, moisture content and plasticity properties. The analysis shows better results in the case of the coarse-grained soils. (The R-  Fig. 6 Variation of the moisture content with the permanent deformation (clay) Fig. 7 Variation of the moisture content in the clay with the resilient modulus and permanent deformation squared value is higher and the obtained values for resilient modulus are similar.) This fact can be explained by the complexity of the nature of fine-grained soils, namely the variation of its properties in terms of plasticity and fines content. At the end of this analysis, one model was selected to rating the geomechanical performance. The MEPDG model shows better results for both materials. Indeed, this choice is also justified by its universal character, which means that can be applied to all types of materials. As expected, the comparison of the results of the resilient modulus shows that the values increase with the UIC classification: higher resilient modulus is related to wellgraded materials (QS2 and QS3).
The geomechanical classification is based on the parametric studies of the resilient modulus (applying the MEPDG model) and the permanent deformation (applying Chen's model). Moreover, other materials were added to this analysis because of the available information about the moisture content. The classification presents a wide range of results in terms of permanent deformation and resilient modulus. This range is a reflection of the variation of the properties of the materials with the same classification (QS1, QS2 and QS3) that influence both types of deformations. Indeed, this work is the first attempt to obtain a novel geomechanical classification that can be a helpful tool in the modelling of the foundation of the railway structures since gather information about the UIC classification (which is a reference in the railway works) and permanent and resilient deformations. Nevertheless, several limitations are recognized in this proposal, which were discussed throughout this paper. These limitations should be overcome in future work updating this proposal with a wide range of geomaterials properly characterized (physical and mechanical characteristics) and tested under different state conditions (deviation to the optimum compaction conditions), which should include variations of water content /suction/degree of saturation and density.

Appendix
An extensive list of resilient models based on laboratory tests is presented in this appendix. Although these tables were constructed uniformly (especially regarding the nomenclature of the variables and the parameters and characteristics of the materials), the authors maintained the original symbols described in the original paper. For example, the symbol E r is used to describe the resilient modulus, which is usually represented by the symbol M r .
Some models were developed using certain SI units, and where possible, these units are described in the column ''variables and empirical constants''.
Regarding the resilient modulus tables, [87] defined w c as optimum moisture content (which is equal to W opt ) and w cr as the relation between the actual moisture content and the optimum moisture content. In this analysis, the fines are measured through the passing sieve #200 or through the percentage of fines. Furthermore, the maximum dry unit weight is also defined as maximum dry density. The constants such as A, B, C or k i are defined as fitting, regression, material and model parameters according to the original formulation defined by the author. The symbols p a and p 0 are defined as atmospheric stress/pressure and reference stress and it is different from p u which means unit pressure.
Regarding the stress variables, the symbol h is defined as bulk stress or the sum of principal stresses, and the deviator stress can be represented by the symbols q, q r and r d . The symbol q is different from q m which represents the mean value of the deviator stress (Tables 17-21).         [67] model.
The last term related to the suction parameter is normalized with the net bulk stress instead of the atmospheric pressure.
The model does not consider the effect of the pore-water pressure under saturated conditions   When the soils get saturated and matric suction is equal to zero (high moisture content) the resilient modulus is close to zero. Low moisture content values lead to high (unrealistic) resilient modulus.
The suction was determined through the paper filter method    [70].
The model is based on tests carried out in several types of coarse-grained granular materials typically used in unbound layers of road and railroad pavements.
The density was incorporated in the model by the porosity parameter. Materials classified as natural gravels.
The study also includes the crushed gravel and crushed rock is a parameter that depends on the mineralogy of the particles, shape of the grains and particle size curve (in the case of granular materials. In the case of the cohesive materials; it also depends on the percentage of fines and Atterberg limits  The effect of stress on the resilient modulus is only considered by the bulk stress (sum of the principal stresses). It is a simplified version of the universal model. This model is based on laboratory triaxial testing (constant confining pressure) [70].
The model fits well in the case of tests with constant confining pressure and shows poor results in the case of variable confining pressure [53] [92] is the mean normal stress; q is the deviator stress; e v is the volumetric strain and e s is the shear strain; K is the bulk modulus; G is the shear modulus This model uses shear-volumetric approach. This is a nonlinear elastic model that takes into account the effect of the stress path and it is expressed in terms of the bulk and shear modulus, that are stress-dependent. The model is elastic and follows the application of the theorem of reciprocity. The shear and volumetric components of stress and strain are treated separately.
The granular materials show an inelastic response and this elastic model can lead to inaccurate predictions.
The model is based on variable confining pressure. Sweere [94] used this model but keep the relationship between shear and volumetric strains with stresses independent from each other Table 21 continued Ref.
The power-law model was selected for the representation of the relationship between resilient modulus and deviator stress Table 21 continued Ref.

Equation model
Variables and empirical constants  h is the bulk stress; s oct is the octahedral shear; p a is atmospheric pressure; k 1 , k 2 , k 3 , a, b and k s are fitting parameters; b = ln(-b/a) The first term of the expression represents the environment related to the contribution of changes in moisture content