Determination of Permanent Deformations of Non-Cohesive Soils in Pavement Structures under Repeated Traffic Load

: One of the main types of distress in pavement structures is rutting, which may also reduce serviceability significantly. Most design methods typically attribute rutting to the asphalt layer alone, proposing that it can be managed by controlling vertical deformation or stress at the subgrade’s top. Furthermore, these methods frequently lack precise measurements for rut depth. On-site measurements show that the majority of permanent deformation occurs in the unbound layers beneath the asphalt; attention should be directed towards these layers. In recent literature, there are calculation methods that account for accumulating strains in soils. However, further investigation is needed regarding the effect of soil properties and the significance of the pavement cross-section. The literature is also somewhat contractionary regarding the origin of permanent deformations. In this research, the residual settlement of soils (base, subbase, and subgrade) under flexible pavement systems was analyzed due to the repeated traffic loads. Rut depths were calculated and analyzed using the High-Cycle Accumulation (HCA) model. The different behaviour in each course of the pavement system is discussed to reveal their contribution to total ruts. The effect of the grain size distribution of the subgrade was analyzed, and its significance on the rutting depths was demonstrated. Standardized pavement cross-sections with different asphalt thicknesses were evaluated, and the calculated settlements of the pavement originating from the ground during the design lifetime are also presented. It is shown that, with the same number of repetitions, the settlements calculated in each traffic load class are proportional to the thickness of the asphalt course. The contribution of the base, subbase, and subgrade courses to the total settlement is also presented for different subgrade types and traffic load classes.


Introduction
The goal when designing a pavement is to find the most economical combination of pavement material and thickness by also considering the local climatic and ground conditions in such a way that the pavement is capable of bearing the traffic over its design lifetime without restricting its functionality.One of the main types of distress in pavement structures is rutting.It was revealed in [1] that an approximately 28 mm deep rut will reduce the present serviceability index (PSI) of a relatively new pavement to an unacceptable classification.It is therefore crucial to look at the mechanism, possible reasons, as well as techniques to reduce the rutting in order to offer long service life and safe pavement structures [2], since besides serviceability, rutting is considered a major safety issue for vehicles [3].
According to a recent literature review [2], rutting is considered to be of asphalt concrete origin only.Also, some of the recent studies [4][5][6] focused on the characteristics and properties of the hot mixed asphalt to reduce rutting sensitivity.However, instigations at the AASHTO road test revealed that 68% of the permanent deformations occur not in the asphalt concrete but in the underlying base, subbase, and subgrade layers [7].In a more recent study [8], comparable findings emerged from on-site assessments, indicating that base and subbase deformations account for 73% to 90% of the overall rut depth.To reduce the rutting sensitivity of pavement structures, the research and calculation method should concentrate on the unbound soil layers instead of the asphalt concrete.
According to a recent review of rutting prediction models [9], there is no widely recognized and documented method for calculating permanent deformations in unbound granular layers (base, subbase courses, and subgrade) beneath asphalt or for determining rut depth.Ref. [10] addressed only the deformation of the granular base course, contradicting the findings of [11].In other instances [12], the material constants or parameter calibrations are not documented, making these models unsuitable for different boundary conditions.In [13,14], the permanent deformations of the unbound layers were determined using state-of-the-art FE analysis.The material constants were calibrated based on cyclic triaxial laboratory tests.However, both studies used a very low number of cycles (only up to 1000 cycles), which is significantly smaller than the millions of axle repetitions usual in pavement structures.Other calculation methods [15] used a state-of-the-art 3D-FEM analysis to quantify the rutting depth.However, only the AC-layer has a constitutive model to enable calculating permanent deformations.The underlying soils are modelled with the Drucker-Prager soil model, so that no cyclic-induced strain can appear in the soils.
In research focusing on rutting, it is crucial to prioritize methods that assess the permanent deformation of the structural layers due to cyclic loading and account for diverse boundary conditions.Since the literature is somewhat controversial regarding the origin of ruts and the calculation methodology, the residual settlements of the unbound layers in the flexible pavement systems induced by the repeated traffic loads were investigated in this paper.The accumulated strains due to cyclic loading were determined using the HCA (High-Cycle Accumulation) model [16].The theoretical basis of the HCA model is presented in Section 2.1.To calculate the rutting depth, a novel calculation model presented in [17] was used, which is suitable for addressing these problems.
This study aims to elucidate how ruts develop over the design lifetime as a result of axle repetitions.A key focus was to determine the ratio and contribution of each specific unbound pavement layer to the overall rut depth, particularly emphasizing the subgrade, given its natural variability in quality.The research further explores how the grain size distribution characterized by the coefficient of uniformity and mean grain size of the subgrade influences the total rutting of the pavement structure.Additionally, the study examines how selecting different standard cross-sections and pavement thicknesses can affect the serviceability of the pavement.The novelty of this research lies in applying the HCA model for pavement structures.In this way, the behaviour of the different structural layers over the design lifetime and depth under the pavement can be investigated.By analyzing the settlements across different subgrade types and traffic load classes, the study provides insights into the development of rutting and informs strategies for mitigating it through informed design choices.

Explicit Calculation Method and HCA Model
Implicit methods that directly calculate soil deformations under repeated loading are unsuitable for a large number of cycles.This is because errors in calculation (integration) and material models accumulate with each load cycle, making these methods reliable only for N < 50 cycles [18].Hence, explicit approaches are preferred for high-cycle loading scenarios, which give the load repetition-accumulated strain (N − ε acc ) relationship on an empirical basis.One of the most extensively documented material models for evaluating permanent deformations under repeated loads, supported by numerous laboratory tests, is the High-Cycle Accumulation (HCA) model, introduced by [16] and validated in [19].The model's detailed presentation is available in [20]; here, only the essential calculation steps pertinent to this research are discussed.If cycles are applied with a low loading frequency, the inertia forces are negligible, leading to quasi-static cyclic loading; at higher frequencies, the loading becomes dynamic due to significant inertia forces.The accumulation model presented in this section based on strain amplitude, is applicable regardless of whether the loading is quasi-static or dynamic.In this research, traffic cyclic loading was considered quasi-static.
The basic principle of the HCA model is to combine the increment-based .σ − . ε advanced implicit material model and the cycle number-accumulated strain (N − ε acc ) explicit model.In the first step, state parameters such as the void ratio, the mean effective normal stress, the stress state and the cyclic strain amplitude (ε ampl ) are determined with a conventional advanced material model.Subsequently, the strain amplitude serves as input for calculating the rate of strain accumulation using the explicit model (Equation ( 1)).In drained conditions, if the mean effective normal stress is not changing and no static plastic strains develop, Equation (1) will directly give the plastic strain rate. .
where the following variables are used: . ε acc is the plastic strain rate per load cycle; f ampl is the function describing the influence of strain amplitude; . f N is the function describing the preloading history (memory of soil); f e is the function describing the influence of void ratio (density); f p is the function describing the influence of average mean pressure; f Y is the function describing the influence of the average stress ratio; f π is the function describing the influence of polarisation changes.
The permanent deformation resulting from the ∆N cycle package is determined by Equation (2) using the plastic strain rate for a constant strain amplitude, ε ampl .As soil structure and state parameters undergo long-term changes, the explicit calculation is periodically halted, and the implicit analysis is initiated in control cycles.This process ensures the continuation of explicit calculation for another ∆N cycle package, incorporating the "updated" actual value of ε ampl .
The HCA model provides the accumulation rate as a vector quantity.Thus, besides the scalar magnitude of the plastic strain, a flow rule is required to define the direction of accumulation (m), delineating the deviatoric ε q and volumetric ε v components of the strain.These components are instrumental in determining the vertical strain elements essential for settlement analysis. .
This study applied the Modified Cam Clay (MCC) flow rule to the subgrade, as suggested by findings from [21].Tests conducted on crushed stone base course revealed anisotropic behaviour under compaction, unlike sand [22], rendering the MCC flow rule inapplicable.Given that pavements commonly incorporate crushed stone in their base and subbase courses, an anisotropic flow rule based on [22] was adopted for this investigation.
The scalar component of the accumulation rate (Equation (1)) can be calculated as the product of five empirical functions: f ampl , .f N , f p , f Y , and f e representing the influence of strain amplitude, load history, mean normal stress, average stress ratio, and void ratio, respectively.These factors are described by seven material constants, along with the critical friction angle φ cc , necessary for defining the flow rule and f Y .Notably, φ cc may differ from the critical friction angle φ c obtained from monotonic CU triaxial tests.This study determined φ cc , based on recommendations from [21].
For determining the HCA parameters, a simplified calibration procedure outlined in [23,24] was employed.This procedure was developed for clean subangular quartz sand and gravel soils devoid of fines content, subject to certain constraints: the coefficient of uniformity must fall between 1.5 and 8.0, d 50 should range from 0.1 mm to 3.5 mm, and the maximum cycle number should not exceed N max = 2 × 10 6 .While the laboratory's utilization of 2 million cycles represents the upper limit of technical reliability, it surpasses the cycle counts in most literature.However, road pavements may experience a higher number of load cycles.Given the uncertainty regarding soil behaviour beyond 2 × 10 6 cycles, caution and critical evaluation are warranted for assumptions associated with N > 2 × 10 6 , as the validity of the C i parameters becomes uncertain.The HCA parameters for the crushed stone base course and subbase course were established following recommendations from [22].
The f N function of the HCA model, which describes the relationship between accumulated residual strain and cycle count, comprises both logarithmic and linear components (Equation ( 4)).If the logarithmic part of the equation is more pronounced depending on the grain size distribution and grain shape [23], the C N3 factor will be very small, and the equation will appear as a linear line in a semi-logarithmic coordinate system.If the linear part of the equation is more pronounced, the C N3 factor will be relatively large, and the equation will bend upwards in a semi-logarithmic coordinate system.

Calculation Method
In this research, a simplified laminar model was used on the analogy of [25][26][27] to calculate the settlement of the pavement.The calculation procedure is discussed in detail in [23]; in this section, only the most important steps are summarized.The used calculation model operates akin to conventional laminar models typically utilized for settlement computations, wherein the subsoil is segmented into multiple sub-layers.Firstly, stress distribution under the axle/wheel load is assessed, and the additional vertical load ∆σ 1,i in each sub-layer is determined.Subsequently, the cyclic strain amplitude ε i ampl (averaged from values at the top and base of the sub-layer) is computed at the centerline of each sublayer, employing a nonlinear approach considering mean normal stress, density, and small strain-dependent stiffness.Next, the accumulated residual strain ε i acc is calculated for each sub-layer using ε i ampl as an input parameter, and compression of each sub-layer is obtained by multiplying it with the thickness and associated cyclic flow rule: s 1,i acc = m 1 × ε i acc d i .The principles of the laminar model are summarized in Figure 1.Ref. [23] demonstrates that the simplified procedure provides sufficiently accurate and reliable results for both the resilient and the plastic residual strains.They showed that the difference between settlements given by the simplified and finite element procedures at the end of the design lifetime is approx.2-7% regardless of the subgrade type and the traffic load class.

Analyzed Pavement Structure
The aim of this study was to ascertain the accumulated residual strains in the standardized pavement cross-sections outlined in [28] resulting from the repeated passage of a standard axle.These standardized pavement cross-sections, detailed in Table 1, consist of four layers: a sand and gravelly sand subgrade course, a granular subbase course, a crushed stone base course, and an asphalt course.The subgrade comprises quartz sand and gravelly sand with subangular particles, as illustrated in Figure 2, and is assumed to possess a relative compaction of R.C. = 93%.A 200 mm thick subbase course overlays the subgrade, the material of which is extensively discussed in the thesis by [29], with its particle size distribution depicted in Figure 2 as well.It is presumed that the subbase course attains a relative compaction of R.C. = 95%.The bearing capacity (E2), modulus of the subbase material, estimated at approximately 150 MPa based on finite element calculations [30], suggests an E2,m of 60 MPa atop the subbase.Literature data on the HCA parameters of crushed stone base courses is only available for crushed stone bases that comply with Colombian Highway Specifications, which is in very close agreement with international crushed stone base course specifications, so this material was chosen in the tests of this research.It is assumed that the crushed stone has a relative compaction of R.C. = 96% and a uniform thickness of 200 mm (Table 1).Bearing capacity (E2) modulus of the crushed stone material based on finite element calculation [30] is approx.230 MPa; thus, an E2,b of 120 MPa can be expected on top of the base course, which complies with the international practice of highway construction.The base course is overlain by an asphalt course of traffic load class-dependent thickness (Table 1).

Analyzed Pavement Structure
The aim of this study was to ascertain the accumulated residual strains in the standardized pavement cross-sections outlined in [28] resulting from the repeated passage of a standard axle.These standardized pavement cross-sections, detailed in Table 1, consist of four layers: a sand and gravelly sand subgrade course, a granular subbase course, a crushed stone base course, and an asphalt course.The subgrade comprises quartz sand and gravelly sand with subangular particles, as illustrated in Figure 2, and is assumed to possess a relative compaction of R.C. = 93%.A 200 mm thick subbase course overlays the subgrade, the material of which is extensively discussed in the thesis by [29], with its particle size distribution depicted in Figure 2 as well.It is presumed that the subbase course attains a relative compaction of R.C. = 95%.The bearing capacity (E 2 ), modulus of the subbase material, estimated at approximately 150 MPa based on finite element calculations [30], suggests an E 2,m of 60 MPa atop the subbase.Literature data on the HCA parameters of crushed stone base courses is only available for crushed stone bases that comply with Colombian Highway Specifications, which is in very close agreement with international crushed stone base course specifications, so this material was chosen in the tests of this research.It is assumed that the crushed stone has a relative compaction of R.C. = 96% and a uniform thickness of 200 mm (Table 1).Bearing capacity (E 2 ) modulus of the crushed stone material based on finite element calculation [30] is approx.230 MPa; thus, an E 2,b of 120 MPa can be expected on top of the base course, which complies with the international practice of highway construction.The base course is overlain by an asphalt course of traffic load class-dependent thickness (Table 1).The first step in the traditional empirical-based pavement design approach is to determine the design traffic and select the corresponding design traffic class.The design traffic, F100, is defined as the number of passing standard axles (mass of 10 tons).Axle loads that differ from this are converted to the standard axle load using load factors that are defined depending on the vehicle class (buses, articulated buses, heavy trucks, trucks with trailers and semi-trailers).In this research, the analysis was performed using the mean value of the design traffic ranges given in Table 1.

Analyzed Soils and Their Properties
The subgrade was assumed to be composed of sand and gravelly sand soils.The reason behind this is that reliable data on HCA parameters in the literature are mostly available for these soils.Five different types of subgrades were analyzed, which cover well the spectrum of soils for which the simplified calibration procedure presented in [23,24] can be used.The effect of mean diameter, d 50 , was investigated on Subgrade 1, 2 and 3, where the value of mean diameter was changed while keeping the coefficient of uniformity constant (C U = 3).The mean diameter for the three subgrades was 0.2 mm, 2.0 mm and 0.6 mm, respectively.The effect of the coefficient of uniformity was investigated on Subgrades 3, 4 and 5, where the value of the coefficient of uniformity was varied between C U = 3.0-5.0-8.0 while keeping the mean diameter constant (d 50 = 0.60).Grain size distribution curves of the subgrades are presented in Figure 2. Because no Proctor test was available for the subgrade, ρ d max was defined in such a way that the actual void ratio should be e = 0.95e min at 100% modified Proctor compaction following the recommendation of [20].The small-strain shear stiffness, G max , and the maximum oedometer modulus, M max , were determined using the correlation in [31,32] as a function of d 50 and C u .The backbone curve was determined using the correlations of [33] as a function of C u .The material of the subbase was defined according to [29] and the parameters of the crushed stone base were taken from [22,34].The particle size distribution of the analyzed layers is shown in Figure 2, whereas their properties are summarized in Table 2.
Eng 2024, 5 1072 The first step in the traditional empirical-based pavement design approach is to determine the design traffic and select the corresponding design traffic class.The design traffic, F100, is defined as the number of passing standard axles (mass of 10 tons).Axle loads that differ from this are converted to the standard axle load using load factors that are defined depending on the vehicle class (buses, articulated buses, heavy trucks, trucks with trailers and semi-trailers).In this research, the analysis was performed using the mean value of the design traffic ranges given in Table 1.

Analyzed Soils and Their Properties
The subgrade was assumed to be composed of sand and gravelly sand soils.The reason behind this is that reliable data on HCA parameters in the literature are mostly available for these soils.Five different types of subgrades were analyzed, which cover well the spectrum of soils for which the simplified calibration procedure presented in [23,24] can be used.The effect of mean diameter, d50, was investigated on Subgrade 1, 2 and 3, where the value of mean diameter was changed while keeping the coefficient of uniformity constant (CU = 3).The mean diameter for the three subgrades was 0.2 mm, 2.0 mm and 0.6 mm, respectively.The effect of the coefficient of uniformity was investigated on Subgrades 3, 4 and 5, where the value of the coefficient of uniformity was varied between CU = 3.0-5.0-8.0 while keeping the mean diameter constant (d50 = 0.60).Grain size distribution curves of the subgrades are presented in Figure 2. Because no Proctor test was available for the subgrade, ρd max was defined in such a way that the actual void ratio should be e = 0.95emin at 100% modified Proctor compaction following the recommendation of [20].The small-strain shear stiffness, Gmax, and the maximum oedometer modulus, Mmax, were determined using the correlation in [31,32] as a function of d50 and Cu.The backbone curve was determined using the correlations of [33] as a function of Cu.The material of the subbase was defined according to [29] and the parameters of the crushed stone base were taken from [22,34].The particle size distribution of the analyzed layers is shown in Figure 2, whereas their properties are summarized in Table 2.The asphalt course was regarded as a uniform layer, its characteristics defined by thickness, Young's modulus, and Poisson's ratio dependent on the traffic load class.The asphalt stiffness was assumed constant (Ea = 4000 MPa, υa = 0.35), adjusted according to the equivalent temperature [35].The asphalt course was regarded as a uniform layer, its characteristics defined by thickness, Young's modulus, and Poisson's ratio dependent on the traffic load class.The asphalt stiffness was assumed constant (E a = 4000 MPa, υ a = 0.35), adjusted according to the equivalent temperature [35].
The HCA parameters for the subgrade were determined via the correlations established in [23,24] based on d 50 and C U .The input parameters for the crushed stone were determined considering the tests of [22].Because no data are available in the literature for the HCA parameters of the subbase course, the Ci parameters were partly determined by iteration considering the particle size distribution, particle shape and available laboratory test results and literature data [20,[22][23][24].The analysis utilized HCA parameters, which are outlined in Table 3 for reference.

Distribution of Strains
The accumulated residual strains caused by traffic loading along the wheel axis were computed for the base, subbase, and subgrade across five distinct subgrades and seven standard pavement cross-sections.This computation was conducted using the calculation model outlined in Section 2. The cyclic loading is constant throughout the entire loading period and is equal to the standard axle load.The numbers of cycles were determined as the average values from the design traffic ranges provided in Table 1.Due to the fact that the HCA parameters were only calibrated up to 2 million cycles, it is not known whether they are still valid above this range or not; the results in traffic load classes "E"-"K", and "R" should be treated cautiously above this threshold value.Therefore, results above 2 million cycles are shown with dashed lines in the figures.
First, let us examine the ε acc -N relationship in the most important level for "Subgrade-5" and traffic load classes "B"-"D"-"K".The most important levels were considered to be at the top of the crushed stone base (level-1), the lower part of the subbase (level-7) and the top of the subgrade (level-9).In order to present the difference between the behaviour of the subgrade and subbase, a level has also been picked from a greater depth (level-17, z = approx.2 m) where the residual accumulated strain is roughly equal to that of the subbase.The selected levels are shown in Figure 3, whereas the calculated accumulated strains are shown in Figure 4.The scaling of the axes in Figure 4 was chosen to be the same so that the results can be compared more easily; hence, the differences are more visible.
Eng 2024, 5  Figure 4 shows that the accumulated strains in the base course are very minor, they are an order of magnitude smaller than the strains in the subbase and subgrade courses.The largest strains, develop in the topmost level (level-9) of the subgrade (Figure 4c).An interesting tendency that is only observed in the subgrade can be noticed here that the rate of accumulation starts to decrease in semi-logarithmic scale after approx.650,000 cycles.The reason behind this is that due to the large residual strains (5-8%), the subgrade is compacted significantly, and the particles are rearranged In such a dense structure that further compaction is hardly achievable.This appears in the equation of the HCA model in such a way that the new and significantly smaller void ratio e resulting from the volumetric portion of the residual strain (ε v ) becomes very close to the C e parameter, which is the void ratio corresponding to the densest state.Consequently, f e becomes ~0, thus .ε acc is significantly reduced, i.e., the rate of accumulation decreases.On the other hand, the stiffness increases as the void ratio decreases, so the resilient cyclic strain amplitude ε ampl obtained by the implicit part of the calculation decreases, which also reduces the accumulation rate .ε acc through f ampl .
Comparison of level-7 and level-17 shows that strains are very similar in them until approx.10 7 cycles in traffic load classes "B"-"D"-"K", but they exhibit significantly different behaviour in the last cycle package in traffic load class "K".In the subbase, mostly the linear component of Equation ( 4) determines the strains due to the very high number of cycles, which results in a concave function in the semi-logarithmic plot.Accordingly, the effect of the N-dependent C N2 component is negligible, and thus C N3 component, which is independent of N, governs the equation.In the subgrade in level-17, however, ε ampl is already <10 −5 in the last cycle package due to the significant compaction, and as a result, further strain accumulation will not occur [20,36].Although the behaviour in levels 7 and

The Effect of Subgrade
As Figure 4 also shows, the relationship of ε acc -N is shaped by very complex processes and depends on several boundary conditions.The goal of this research is to describe the behaviour of the entire unbound part of the pavement system and thus to determine the expected residual settlements.Figure 5 shows the calculated settlements in each type of subgrade for traffic load class "C".The particle size distribution of the subgrade has a significant effect on the settlement of the whole pavement structure; the difference between the largest and smallest calculated settlements is almost five times.The reason for this is the effect of the C u and d 50 parameters describing the particle size distribution on the C Ni HCA parameters.
Let us have a look at the effect of the coefficient of uniformity.For "Subgrade-1", "-4" and "-5", the coefficient of uniformity is increased from 3.0 to 5.0 and then to 8.0, respectively, while keeping d 50 = 0.6 constant.It is clearly visible that increasing C U increases the settlement.The reason behind this is that more well-graded soils can be compacted better due to a unit of compaction energy because particles of different sizes can fit more easily into the voids between the particles.This appears in the C Ni parameters of the HCA model that the difference in, e.g., C N2 and C N3 parameters between "Subgrade-1" and "Subgrade-5" is nearly twofold and in parameter C N1 is nearly threefold.Accordingly, the rate of accumulation in more well-graded soils is high not only in the logarithmic component governing the beginning of loading but also in the linear component governing at higher cycles.The obtained result is a bit surprising at first since, in geotechnical practice, the use of soils with a low coefficient of uniformity is generally avoided because their compaction is difficult at the construction site, and their load-bearing capacity is also lower.Although this statement is true for installation, the opposite conclusion can be drawn if the subsequent compaction due to traffic is considered.
Eng 2024, 5 10 Let us have a look at the effect of the coefficient of uniformity.For "Subgrade-1", 4" and "-5", the coefficient of uniformity is increased from 3.0 to 5.0 and then to 8.0, r spectively, while keeping d50 = 0.6 constant.It is clearly visible that increasing CU increase the settlement.The reason behind this is that more well-graded soils can be compacte better due to a unit of compaction energy because particles of different sizes can fit mo easily into the voids between the particles.This appears in the CNi parameters of the HC model that the difference in, e.g., CN2 and CN3 parameters between "Subgrade-1" an "Subgrade-5" is nearly twofold and in parameter CN1 is nearly threefold.Accordingly, th rate of accumulation in more well-graded soils is high not only in the logarithmic compo nent governing the beginning of loading but also in the linear component governing higher cycles.The obtained result is a bit surprising at first since, in geotechnical practic the use of soils with a low coefficient of uniformity is generally avoided because the compaction is difficult at the construction site, and their load-bearing capacity is als lower.Although this statement is true for installation, the opposite conclusion can b drawn if the subsequent compaction due to traffic is considered.
The effect of particle diameter was analyzed in "Subgrade-2", "-3", and "-1", whe d50 changes from 2.0 mm to 0.6 mm and then to 0.2 mm, respectively, while CU = 3.0.It ca be concluded that strains decrease with increasing particle diameter.The largest strain were calculated in "Subgrade-1", which has the smallest particle size, and the smalle settlement was obtained in "Subgrade-2".There is an almost three-fold difference betwee the calculated settlements.The obtained result is in accordance with the geotechnical pra tice that soils with larger grain sizes have generally higher strength and load-bearing c pacity.This experience is described well by the CNi parameters of the HCA model.
The same tendency was observed for all other traffic load classes as for load class "C The effect of particle diameter was analyzed In "Subgrade-2", "-3", and "-1", where d 50 changes from 2.0 mm to 0.6 mm and then to 0.2 mm, respectively, while C U = 3.0.It can be concluded that strains decrease with increasing particle diameter.The largest strains were calculated in "Subgrade-1", which has the smallest particle size, and the smallest settlement was obtained in "Subgrade-2".There is an almost three-fold difference between the calculated settlements.The obtained result is in accordance with the geotechnical practice that soils with larger grain sizes have generally higher strength and load-bearing capacity.This experience is described well by the C Ni parameters of the HCA model.
The same tendency was observed for all other traffic load classes as for load class "C" that settlements increase from Subgrade 2-3-1-4 and 5.The total settlements and their ratios in the soils with the extreme values of d 50 and C U ("Subgrade-1"-"Subgrade-2" and "Subgrade-1"-"Subgrade-5", respectively) in each traffic load class are presented in Table 4.As the table shows the subgrade has a fundamental role in the total settlement.The settlement ranges from 2.1 times to 5 times the value for "Subgrade-2" with different coefficients of uniformity and from 1.7 times to 3.2 times the value for "Subgrade-2" with different mean grain sizes.The largest and smallest differences between the effect of subgrades were observed in traffic load class "D" and "R", respectively.The reason for this is that in the low-traffic load classes, large differences cannot yet emerge from the C N3 parameters due to the relatively low number of axle passes.The effect of the subgrade is the smallest in traffic load classes "K" and "R".Although it would be possible for the differences arising from the C Ni parameters to become distinctive due to the very high number of cycles in these classes, such large strains develop in the subgrade that further strain accumulation is significantly limited by the compaction (see the behaviour in Figure 4c).Accordingly, the subgrade has the greatest effect in the case of load classes "C" and "D", where the number of cycles is already high enough to trigger the differences in the soil behaviour, but the reduction in the rate of accumulation resulting from compaction only occurs in a thin layer at the top of the subgrade.Table 4.Total settlements and their ratio for subgrades with the lower and upper bound particle size distributions.

The Effect of Traffic Load Class
Figure 6 shows the residual settlements as the function of the cycle number for the different traffic load classes in the case of subgrade 3. The figure shows that the logarithmic behaviour is dominant up to approx.10 5 cycle number, which appears as a roughly linear line in the semi-logarithmic plot.With the increasing number of cycles and g A -parameter, the logarithmic behaviour becomes more and more overshadowed, and the C N2 component of the HCA equation becomes more and more insignificant.With the further increase in cycle number above 10 5 , the linear component and C N3 parameter become dominant, so the function "curves upwards", and a concave curve is obtained in a semi-logarithmic plot.
number of cycles in these classes, such large strains develop in the subgrade that further strain accumulation is significantly limited by the compaction (see the behaviour in Figure 4c).Accordingly, the subgrade has the greatest effect in the case of load classes "C" and "D", where the number of cycles is already high enough to trigger the differences in the soil behaviour, but the reduction in the rate of accumulation resulting from compaction only occurs in a thin layer at the top of the subgrade.Table 4.Total settlements and their ratio for subgrades with the lower and upper bound particle size distributions.

The Effect of Traffic Load Class
Figure 6 shows the residual settlements as the function of the cycle number for the different traffic load classes in the case of subgrade 3. The figure shows that the logarithmic behaviour is dominant up to approx.10 5 cycle number, which appears as a roughly linear line in the semi-logarithmic plot.With the increasing number of cycles and g A -parameter, the logarithmic behaviour becomes more and more overshadowed, and the CN2 component of the HCA equation becomes more and more insignificant.With the further increase in cycle number above 10 5 , the linear component and CN3 parameter become dominant, so the function "curves upwards", and a concave curve is obtained in a semi-logarithmic plot.The concept of empirical pavement design, based on traffic load classes, aims to ensure that the stresses induced by varying numbers of axle passes result in equivalent fatigue over the asphalt's entire design lifespan.However, as depicted in Figure 6, this equivalence does not hold true for settlements.Calculated settlements at the conclusion of the design lifespan escalate with higher traffic load classes.Simultaneously, settlements are observed to be smaller for the same cycle count in higher load classes.This correlation is logical since thicker asphalt courses in higher traffic load classes distribute wheel loads more efficiently, thereby reducing stresses reaching the subgrade and resulting in smaller ε ampl .Consequently, accumulated strains are lesser in better load-bearing traffic classes for equivalent cycle counts.For instance, after 65,000 axle passes, settlements at the end of the design lifespan are approximately four times greater in traffic load class "A" compared to those calculated for the same cycle count in class "R".However, by the conclusion of the design lifespan for class "R", settlements are approximately 5.4 times greater, despite enduring approximately 460 times more axle passes.
Figure 7 shows the normalized settlements according to Equation (5) in the case of "Subgrade-3".Accordingly, if the calculated settlement in the ith traffic load class is reduced by the ratio of the asphalt thickness of class "i" and the reference class "R", the normalized settlement will fall on one curve, which has a variance of s 2 = 0.83.This observation is also true for all the other subgrade types.This means that if the boundary conditions (e.g., subgrade type) and the number of cycles are the same, the settlement of the different traffic load classes will be proportional to the ratio of the thickness of the asphalt course as it appears in Equation ( 5). (5)

Ratio of Settlements
Another important aspect is to determine how each soil layer contributes to the settlement and how the settlement accumulates as a function of depth.This is illustrated in Figure 8 for traffic load class "C" and "Subgrade-1".Most of the settlement in the selected example (approx.90%) occurs in the upper approx. 1 m thick zone of the subgrade.The contribution of the base and subbase courses and the deeper zone of the subgrade is marginal to the total settlement.
Eng 2024, 5 1078 The concept of empirical pavement design, based on traffic load classes, aims to ensure that the stresses induced by varying numbers of axle passes result in equivalent fatigue over the asphalt's entire design lifespan.However, as depicted in Figure 6, this equivalence does not hold true for settlements.Calculated settlements at the conclusion of the design lifespan escalate with higher traffic load classes.Simultaneously, settlements are observed to be smaller for the same cycle count in higher load classes.This correlation is logical since thicker asphalt courses in higher traffic load classes distribute wheel loads more efficiently, thereby reducing stresses reaching the subgrade and resulting in smaller εampl.Consequently, accumulated strains are lesser in better load-bearing traffic classes for equivalent cycle counts.For instance, after 65,000 axle passes, settlements at the end of the design lifespan are approximately four times greater in traffic load class "A" compared to those calculated for the same cycle count in class "R".However, by the conclusion of the design lifespan for class "R", settlements are approximately 5.4 times greater, despite enduring approximately 460 times more axle passes.
Figure 7 shows the normalized settlements according to Equation (5) in the case of "Subgrade-3".Accordingly, if the calculated settlement in the ith traffic load class is reduced by the ratio of the asphalt thickness of class "i" and the reference class "R", the normalized settlement will fall on one curve, which has a variance of s 2 = 0.83.This observation is also true for all the other subgrade types.This means that if the boundary conditions (e.g., subgrade type) and the number of cycles are the same, the settlement of the different traffic load classes will be proportional to the ratio of the thickness of the asphalt course as it appears in Equation ( 5).

Ratio of Settlements
Another important aspect is to determine how each soil layer contributes to the settlement and how the settlement accumulates as a function of depth.This is illustrated in Figure 8 for traffic load class "C" and "Subgrade-1".Most of the settlement in the selected example (approx.90%) occurs in the upper approx. 1 m thick zone of the subgrade.The contribution of the base and subbase courses and the deeper zone of the subgrade is marginal to the total settlement.Figure 9 illustrates the percentage contribution of each soil layer to the total settlement as a function of the number of cycles for the different load classes and in the case of "Subgrade 1".It can be observed that the combined contribution of the crushed stone base course and the subbase course from load classes "A" to "E" is only approx.5% to the total settlement, and significant deformations develop mostly in the subgrade.This conclusion is true for the entire loading process; the contribution of these layers to the total settlement only changes negligibly (2-4%) with the number of cycles.Generally, in the higher load classes the percentage contribution of the crushed stone and subbase layer is higher in the initial stage of the loading.The reason for this is that the smaller stresses reaching the subgrade result in a smaller ε ampl , which yields significantly smaller deformations in the subgrade at the beginning of the loading due to the nearly quadratic relationship with Campl.For class "R", the combined contribution of the crushed stone and subbase gradually decreases from the initial approx.13% to approx.6% to 10 6 cycles, then it increases again to 15% until the end of the design lifetime.
The reason for the shift in the proportions is that the rate of accumulation in the upper coarse-grained layers is more or less constant even at very high cycle numbers, but further strain accumulation is severely limited in the upper part of the subgrade due to the high degree of compaction.Since the rate of accumulation in the subgrade decreases while it remains basically constant in the base and subbase courses, the proportions will shift towards the latter.The general shape of the proportion curves is the same in each traffic load class for the different subgrade types, but the proportions are somewhat different.The reason for this is that "Subgrade-5" is the most sensitive subgrade type, and it can be compacted the most easily, which is why the largest deformations and settlement occur in this type of soil.Since the input parameters for the crushed stone and the subbase layers were not changed, "Subgrade-5" will yield the highest contribution (Figure 10).In the case of "Subgrade-2", the effect is exactly the opposite; this soil type is the most resistant to cyclic deformations, so the settlement proportion of the crushed stone and subbase layers is higher than for other subgrade types (Figure 11).Accordingly, the extreme values occur in the case of load class "A" with "Subgrade-5" (98% of the settlement is caused by the subgrade) and load class "R" with "Subgrade-2" (75% of the settlement is caused by the subgrade).Figure 9 illustrates the percentage contribution of each soil layer to the total settlement as a function of the number of cycles for the different load classes and in the case of "Subgrade 1".It can be observed that the combined contribution of the crushed stone base course and the subbase course from load classes "A" to "E" is only approx.5% to the total settlement, and significant deformations develop mostly in the subgrade.This conclusion is true for the entire loading process; the contribution of these layers to the total settlement only changes negligibly (2-4%) with the number of cycles.Generally, in the higher load classes the percentage contribution of the crushed stone and subbase layer is higher in the initial stage of the loading.The reason for this is that the smaller stresses reaching the subgrade result In a smaller ε ampl , which yields significantly smaller deformations in the subgrade at the beginning of the loading due to the nearly quadratic relationship with C ampl .For class "R", the combined contribution of the crushed stone and subbase gradually decreases from the initial approx.13% to approx.6% to 10 6 cycles, then it increases again to 15% until the end of the design lifetime.
The reason for the shift in the proportions is that the rate of accumulation in the upper coarse-grained layers is more or less constant even at very high cycle numbers, but further strain accumulation is severely limited in the upper part of the subgrade due to the high degree of compaction.Since the rate of accumulation in the subgrade decreases while it remains basically constant in the base and subbase courses, the proportions will shift towards the latter.The general shape of the proportion curves is the same in each traffic load class for the different subgrade types, but the proportions are somewhat different.The reason for this is that "Subgrade-5" is the most sensitive subgrade type, and it can be compacted the most easily, which is why the largest deformations and settlement occur in this type of soil.Since the input parameters for the crushed stone and the subbase layers were not changed, "Subgrade-5" will yield the highest contribution (Figure 10).In the case of "Subgrade-2", the effect is exactly the opposite; this soil type is the most resistant to cyclic deformations, so the settlement proportion of the crushed stone and subbase layers is higher than for other subgrade types (Figure 11).Accordingly, the extreme values occur in the case of load class "A" with "Subgrade-5" (98% of the settlement is caused by the subgrade) and load class "R" with "Subgrade-2" (75% of the settlement is caused by the subgrade).

Conclusions
This study assessed the accumulated residual deformations of soil layers under flexible pavements using both implicit and explicit calculation methods.The research analyzed five types of subgrades with different particle size distributions and seven standardized pavement cross-sections to determine their effects on settlement.The key findings include: 1.
Influence of subgrade type: The subgrade grain size distribution characteristic significantly impacts the rutting depth, with well-graded soils showing higher accumulated strains due to better compactability.Increasing the mean particle size (d 50 ) while maintaining a constant coefficient of uniformity (C U ) positively affects settlement.

2.
Layer contribution: The subgrade accounts for the majority of the total settlement (75-98%), though its contribution decreases towards the end of the design lifetime.Most strains originate in the upper part of the subgrade, suggesting that improving this zone can most effectively reduce rutting.

3.
Traffic Load Class: By a given number of load repetitions the accumulated strains will be smaller in the cross sections with thicker asphalt layers due to a better stress distribution.The settlements, however, towards the end of design lifetime, increase with higher traffic load classes and are proportional to the thickness of the asphalt layer when boundary conditions and cycle numbers are constant.
Some of the calculated settlements, especially by higher load classes, are so significant that this probably reduces the serviceability of the pavement before the end of the design lifetime.It was observed that by limiting the vertical deviatoric stress, the settlements can be effectively reduced.This can be achieved via a better stress distribution by using thicker asphalt courses, thicker and stabilized base and subbase courses or by using rigid pavement structures, especially by higher traffic load classes.The calculation method used is well suited for parametric studies, and future research will examine additional factors affecting rutting.Future work will also determine the spatial extent of settlements, consider different national specifications, and investigate the effects of varying load cycles and seasonal temperature changes on asphalt stiffness and settlements.By understanding these factors, we can better predict and mitigate rutting in pavement structures, ultimately improving pavement performance and longevity.

Figure 1 .
Figure 1.The main calculation steps of the layered model [17].

Figure 1 .
Figure 1.The main calculation steps of the layered model [17].

Figure 3 .
Figure 3. Location of the 4 selected levels: Level-1 on top of the crushed stone base, levelmiddle of the subbase in the depth of z = 0.325 m, level-9 on top of the subgrade, level 1 subgrade in the depth of z ~ 2 m.

Figure 3 .
Figure 3. Location of the 4 selected levels: Level-1 on top of the crushed stone base, level-7 in the middle of the subbase in the depth of z = 0.325 m, level-9 on top of the subgrade, level 17 in the subgrade in the depth of z ~2 m.

Eng 2024, 5 1075
17 was very similar in the first 10 7 cycles, by the end of the design lifetime, almost half as much residual deformation developed in level-17 as in level-7.

Figure 3 .
Figure 3. Location of the 4 selected levels: Level-1 on top of the crushed stone base, level-7 in the middle of the subbase in the depth of z = 0.325 m, level-9 on top of the subgrade, level 17 in the subgrade in the depth of z ~ 2 m.

Figure 4 .
Figure 4. ε acc in the function of load cycles for four selected levels (a) Level-1 on top of the crushed stone base.(b) Level-7 (subbase).(c) Level-9 on top of the subgrade.(d) Level-17 in the subgrade.

Figure 4 .
Figure 4. ε acc in the function of load cycles for four selected levels (a) Level-1 on top of the crushed stone base.(b) Level-7 (subbase).(c) Level-9 on top of the subgrade.(d) Level-17 in the subgrade.

Figure 5 .
Figure 5.The effect of subgrade on settlements in the case of traffic load class "C".

Figure 5 .
Figure 5.The effect of subgrade on settlements in the case of traffic load class "C".

Figure 6 .
Figure 6.The relationship between settlement and cycle number in the case of "Subgrade-3".

Figure 6 .
Figure 6.The relationship between settlement and cycle number in the case of "Subgrade-3".

Figure 7 .
Figure 7. Normalized settlements with the thickness of the asphalt.

Figure 7 .
Figure 7. Normalized settlements with the thickness of the asphalt.

Figure 8 .
Figure 8. Distribution of settlement with depth in the case of traffic load class "C" and "Subgrade-1".

Figure 8 .
Figure 8. Distribution of settlement with depth in the case of traffic load class "C" and "Subgrade-1".

Figure 9 .
Figure 9.The contribution of each layer to the total settlement as a function of the number of cycles-"Subgrade-1".

Figure 9 .Figure 9 .Figure 10 .
Figure 9.The contribution of each layer to the total settlement as a function of the number of cycles-"Subgrade-1".

Figure 10 .Figure 11 .
Figure 10.The contribution of each layer to the total settlement as a function of the number of cycles-"Subgrade-5".

Figure 10 .Figure 10 .Figure 11 .
Figure 10.The contribution of each layer to the total settlement as a function of the number of cycles-"Subgrade-5".

Figure 11 .
Figure 11.The contribution of each layer to the total settlement as a function of the number of cycles-"Subgrade-2".Figure 11.The contribution of each layer to the total settlement as a function of the number of cycles-"Subgrade-2".

Table 2 .
Geotechnical parameters of the layers.

Table 3 .
HCA parameters of the layers.