Contribution to resilient and permanent deformation investigation of unbound granular materials with different geological origins from Rio Grande do Sul, Brazil

Abstract This article evaluates the resilient modulus and permanent deformation of granular materials of different lithological origins widely used as a pavements base layer in south Brazil. For this, a single particle size distribution was determined for the materials that were subjected to physical, chemical, mechanical characterizations, especially resilient modulus and permanent deformation by repeated load triaxial tests. It was noticed that the denser materials had a higher resilient modulus generated by increase in the sample's stiffness. For permanent deformation this tendency has not been maintained for all materials. Therefore, the granulation and structure of the materials can influence long-term tests. The Guimarães’ model has proven to be adequate for the sample evaluation. For the shakedown research, samples showed accommodation and creep shakedown. The samples that presented accommodation had an increase in the resilient modulus after permanent deformation, while those that presented creep increased or decreased resilient modulus according to the material origin.


Introduction
Granular materials, used as a base layer for flexible pavements, influence the structure performance as a whole.In this layer, when the material receives and supports the stresses of traffic, successively returning to its original state throughout its service life, it is said that this material has a high resilient modulus (RM), preventing the formation of fatigue cracks on the pavement surface.Concomitantly, part of the deformations caused by traffic action are accumulated, called permanent deformations (PD) which, added to the deformations of other layers, are manifested on the pavement surface as rutting (Huang, 2004;Cerni et al., 2012;Erlingsson et al., 2017).
Shakedown theory is used to characterize soils and unbound granular materials employed in pavements (Werkmeister et al., 2001;Werkmeister, 2003;Werkmeister, 2006;Wang & Yu, 2013;Gu et al., 2017;Alnedawi et al., 2019a;Nazzal et al., 2020).According to the theory, the unbound granular materials (UGM) can present three different ranges: range A or plastic shakedown, after a finite number of cicles, the acumulation of permanent deformations reaches a constant; range B or creep shakedown, the permanent strain rate decreases after a number of cycles, but the resilient strains are not constant and the stiffness decrease; range C or incremental colapse, the permanent deformation increases with the cycles.The UGMs can fail by shear or overstressing.
The main parameters that affect the elastic and plastic deformability of granular layers are the active stresses, the reorientation of the main stresses, the history of stresses, number, duration and frequency of loads, degree of compaction, moisture content, particle size distribution, maximum size of aggregates, fines content, type of aggregate and particle shape (Collins & Boulbibane, 2000;Lekarp et al., 2000;Lekarp & Isacsson, 2001;Song & Ooi, 2010;Xiao et al., 2019;Soliman & Shalaby, 2015).Among these factors, for the present study, the lithological origin of the aggregates that compound the granular material is emphasized.

Abstract
This article evaluates the resilient modulus and permanent deformation of granular materials of different lithological origins widely used as a pavements base layer in south Brazil.For this, a single particle size distribution was determined for the materials that were subjected to physical, chemical, mechanical characterizations, especially resilient modulus and permanent deformation by repeated load triaxial tests.It was noticed that the denser materials had a higher resilient modulus generated by increase in the sample's stiffness.For permanent deformation this tendency has not been maintained for all materials.Therefore, the granulation and structure of the materials can influence long-term tests.The Guimarães' model has proven to be adequate for the sample evaluation.For the shakedown research, samples showed accommodation and creep shakedown.The samples that presented accommodation had an increase in the resilient modulus after permanent deformation, while those that presented creep increased or decreased resilient modulus according to the material origin.
Contribution to resilient and permanent deformation investigation of unbound granular materials with different geological origins from Rio Grande do Sul, Brazil According to Xiao et al., (2019), Ba et al. (2015), Alnedawi et al. (2019b) and Alnedawi et al. (2021), granular materials composed of different lithological origins show different deformability behavior, even when the granulometric curve is similar among them.Lima et al. (2017) have analyzed two quarries of the same lithological origin and for the same granulometric curve obtained similar results in terms of resilience and permanent deformation.Nazzal et al. (2020) have conducted RM and PD tests of granular granite, sandstone and limestone materials.The mixtures' granulometric curves were similar, although the behavior in terms of resilience, permanent deformation and shakedown are considerably different.
The objective of this study is to evaluate the resilient and permanent deformation behavior of three unbound granular materials from different lithological origins.
The characterization of the materials that compose the UGM and the results of tests and mathematical modeling of RM and PD and shakedown research are presented.The relevance of this study is related to the extensive use of granular bases regionally, due the abundance and variability of rock masses in the territory of the state of Rio Grande do Sul.Therefore, these materials should be characterized in terms of deformability.Besides, this study contributes to the database of the Brazilian M-E Design Method (MeDiNa), which is being implemented.

Materials
To evaluate the influence of the lithological origin on the resilient behavior, permanent deformation and shakedown of the base course layer, three different types of mineral aggregates were selected.It is noteworthy that such materials came from different geomorphological provinces, including practically all lithological formations in the state of Rio Grande do Sul, Brazil.For this purpose, São Juvenal  The petrographic thin section slides provided the microscopic description of the textural, structural and mineralogical characteristics of the rocks, in addition to the physical and mechanical characterization tests.Table 1 presents the mineralogical composition, rock description and the physical and mechanical indexes.
Based on Table 1, aggregates mechanical performance is consistent with their different types of formation, since SJ and DPA are defined as volcanic igneous rocks, with aphanitic texture, and SBS as plutonic igneous rock, with phaneritic texture.Although plutonic igneous rocks have satisfactory mechanical resistance due to the relative homogeneity of rock bodies in addition to the mineralogical composition that holds minerals of high hardness, the high granulation of their minerals promotes points of weakness in the rock, increasing the occurrence of micro fractures and consequently decreasing the material mechanical strength (Curtis Neto et al., 2018;Back et al., 2021;Adomaki et al., 2021).This fact justifies the superior mechanical performance of SJ and DPA, since they have aphanitic texture and/or fine granulation, thus presenting better distributions of mechanical efforts.
The presence of large percentages of alkali feldspar corroborates to the high abrasive loss of SBS, which contains minerals of high hardness (e.g., feldspar: 6 and quartz: 7) with the rocky matrix and show low tenacity due to the granulation.Similarly, the presence of foliation also influences the SBS mechanical performance since it adds horizontal weakness planes that tend to generate a higher percentage of lamellar aggregates in the crushing process and these particles tend to break in the compaction process (Wojahn et al., 2021).In addition to the abrasive loss, there is high loss due to impact and crushing.
Contribution to resilient and permanent deformation investigation of unbound granular materials with different geological origins from Rio Grande do Sul, Brazil

Experimental program
The granulometric curve adopted for the evaluated UGM is in accordance with a highway road standard (DNER, 1998b).The material from the aforementioned rock deposits went through separation, sieving and mixing procedures in order to obtain the specified UGM.A single granulometric curve was used for the three unbound granular materials, which was included within the limits of the most used specifications for the granular base layers of Brazil's southern region by the federal highway agencies (range C) (DNIT, 2006) and state highway agencies (maximum size 3/4") (DAER, 1991).Figure 2 shows the selected particle size distribution.
The samples were subjected to the compaction test to determine the maximum dry density (MDD) and the optimum moisture content (OMC).The compaction process follows the procedures described in the standards of resilient modulus (DNIT, 2018a) and permanent deformation (DNIT, 2018b), similar to the technical procedures adopted by American Association of State Highway and Transportation Officials (AASHTO, 2004), Australia (AGPT, 2006) and the European Union (BSI, 2004).For this, a three-part cylindrical steel mold with dimensions of 200 mm in height and 100 mm in diameter was used.The compaction energy used was equivalent to that of the modified energy.Samples were considered valid for RM and PD tests if the variation was not superior than ± 1% in relation to the MDD and OMC.
The repeated load triaxial equipment used is located in the Department of Transportation from Universidade Federal de Santa Maria (UFSM).In the RM test, firstly 1,500 cycles are applied for conditioning the sample; then the test proposes the application of a 100 load cycles for each of 18 pairs of confining stress and deviator stress (DNIT, 2018a).The test was performed in triplicate and with a loading frequency of 1 Hz. Figure 3 presents a sample of UGM compacted in the conditions mentioned above, being positioned in the equipment for a subsequent RLT test.
After the resilient modulus test, the test data were adjusted to mathematical models that were representative of the paving materials mechanical behavior.The mathematical models that can describe the behavior of granular materials regarding resilience, pointed out in the technical literature and cited in this study, are in Table 2.For multiple nonlinear analysis, we used the software Statistica v.10.0.228.2, for linear analysis we used the software Microsoft Office Excel 2013.In order to evaluate the plasticization to which the granular materials are subjected, the permanent deformation tests were carried out (DNIT, 2018b).Again, the RLT was used to obtain the experimental data.In this test, 100 conditioning cycles are applied for conditioning phase with the 30×30 kPa stresses pair.It is followed by the application of 150,000 charging cycles of each stress pair.For this study, the test protocol was performed according to Lima et al. (2019), with six samples and application of the following single stage stresses (confining versus deviator): 40×40 kPa, 40×120 kPa, 80×80 kPa, 80×240 kPa, 120×240 kPa and 120×360 kPa.
Among the models for predicting permanent deformation of granular materials, the models of Barksdale (1972), Monismith et al. (1975), Uzan (1985) and Tseng & Lytton (1989) stand out.However, for the coherent characterization of tropical materials and their particularities, the model that seems more adequate is the Guimarães' model (Guimarães et al., 2018;Lima et al., 2020), a model used in this research and addressed in Brazilian regulations (Equation 8).The parameters ( ) Where: ( )

Results and analysis
The results are presented in four distinct topics, in the following order: results of sample compaction; resilient modulus, permanent deformation and the relationship between void ratio and deformability.

Compaction
The results from the compaction tests in this study are shown in Figure 4. SJ basalt granular material reached higher densities, due the high specific gravity of the basalt aggregates.The DPA rhyodacite showed the highest OMC and lowest density among all samples.According to Paiva (2017), this rock presents devitrification of the rock matrix, which is replaced by clay minerals, thus increasing porosity and absorption, justifying the high optimum moisture content.
Where: RM : resilient modulus; 3 σ : confining stress; d σ : deviator stress; θ : first stress invariant (  The SBS syenogranite, on the other hand, has the lowest OMC among the mixtures due to its lower absorption, as mentioned in Table 1.

Resilient modulus
In order to evaluate the stiffness properties of the UGMs samples, after the repeated load triaxial tests the laboratory results were analyzed and submitted to mathematical modeling by means of models k-3 σ , k-d σ , k-θ, Compound, Universal, Witczak and Witczak and Uzan.
The regression parameters are presented in Table 3 with the values of the coefficient of determination (R 2 ) and the value of the linear resilient modulus, that is, the average of the RM values obtained for each stress pair after performing the mathematical modeling.For the Compound model, coefficient k 2 referring to the confining stress action has higher impact on the RM value to the detriment of coefficient k 3 related to the deviator stress.However, the fact that the coefficient k 3 has a positive value indicates that, by increasing the deviator stress, an increase in the resilient modulus occurs.The Compound model is the most discussed in this study since it is currently used in the Brazilian M-E Design Method (MeDiNa) for characterization related to the stiffness of subgrade soils and granular materials (Guimarães & Motta, 2016;Freitas et al., 2020;Lima et al., 2020;Sagrilo, 2020;Pascoal et al., 2021;Zago et al., 2021;Pascoal et al., 2023).
The analysis of the Universal model indicates that the increase of the first stress invariant collaborate to the RM of the material, by the positive value of k 2 .On the other hand, the coefficient k 3 presented a negative value due to the increase in the octahedrical shear stress, causing a lower RM value.This is similar to the Witczak's and Uzan's model, in which coefficients k 2 and k 3 are related to the same stresses as the Universal model.The model by Witczak also demonstrate that the increase of first stress invariant reflects in a higher RM, due to the positive value of k 2 for all materials; although the increase of the deviator stress does not strongly affect the RM, due to the negative and low value of k 3 .

Permanent deformation
In order to understand the behavior of the granular material subjected to damage by permanent deformation, the UGM of different lithologies were submitted to the RLT test of long duration.Subsequently, the shakedown analysis was carried out according to Werkmeister (2003), the Guimarães' model was implemented to the mixtures and loss or gain of stiffness of the samples was verified after the stress history of the PD test.
Figures 6, 7 and 8 show the graphs representing the permanent deformation tests of UGM SJ, DPA and SBS, respectively.The graphs at the left (a) show the accumulated PD versus number of cycles, so the increase rate of the curve gives evidence of stabilization or not of these deformations; the right graph line shows the behavior of the permanent deformation increase rate over the accumulation of deformations.Contribution to resilient and permanent deformation investigation of unbound granular materials with different geological origins from Rio Grande do Sul, Brazil For SJ basalt the highest permanent deformations recorded were observed in tests whose stress pairs are 120×360 kPa and 80×240 kPa, confining versus deviator, respectively, reaching values close to 2.3 mm at the end of the tests (Figure 6a).The pair 120×240 kPa presented considerably smaller deformations than the pairs 120×360 kPa and 80×240 kPa.The smallest deformations were found in the 40×40 kPa e 80×80 kPa tests.This behavior indicates that the permanent deformations magnitude is related to the ratio among 3 / d σ σ in a directly proportional relation.
Observing the increase rate graph in PD by the cumulative vertical deformation of the SJ it can be seen that only tests whose stress ratio 3 / d σ σ is equal to one entered shakedown (i.e., type A behavior) according to the parameters of Dawson & Wellner (1999) and Werkmeister (2003) with the permanent deformation increase rate in 10 -7 x10 -3 meters per load application cycle in 150000 cycles.Those tests are presented in the Figure 6b with filled markers.The other pairs presented creep shakedown, i.e., behavior B. Projecting the PD values with the Guimarães' model, the pair 40×120 kPa would reach the rate of 10 -7 x10 -3 meters per load application cycle in 155000 cycles, while the same would happen to the pair 120×240 kPa at 185000 cycles.
For the DPA rhyodacite the deformation results were similar to that observed with basalt, in relation to the pairs that deformed more or less, as well as the value of higher accumulated deformation, close to 2.5 mm (Figure 7a).The 40×40 kPa pair reached the accommodation -behavior A, as it possible to see in Figure 7b; the same happened to the pair 80×80 kPa by the projection of the Guimarães' model for this material, reaching 10 -7 x10 -3 meters per load application cycle.The other pairs present behavior B.
Syenogranite, an aggregate with a different texture from the others, presented considerably higher plastic deformability than the other analyzed materials.The highest deformation observed was approximately 5 mm, although the densities obtained in the compaction test were close to that for the rhyodacite (Figure 8a).This behavior seems to be related to this rock's inferior mechanical performance in relation to the others with respect to the laboratory tests of mechanical characterization, as previously shown in Table 1.The stress pairs that showed a tendency to shakedown with the highest number of cycles are the pairs with 3 / d σ σ is equal to one -40×40 kPa and 80×80 kPa (Figure 8b).Other pairs present type B behavior.
The materials were, by multiple nonlinear regression analysis, mathematically characterized as to their PD by the Guimarães' model.The parameters obtained for the specific permanent deformation are in Table 4. From the parameters and the coefficient of determination, it is evident that the Guimarães' model is representative of the deformability of the materials under study; also noticeable is the strong impact of the deviator stress in the PD, represented by the high values of 3 ψ , a trend already observed in Figures 6, 7  The stiffness gain of samples submitted to stress history was evaluated by the resilient modulus test after the permanent deformation test.It was possible to compare the average RM value after the test for each stress pair and the value obtained for the conventional RM test.The average RM for Compound Model after each PD test is shown in Figure 9 in a bar format, while the horizontal line represents the conventional RM test.In addition, in the center of each bar, is it shown the shakedown behavior of each sample, in this case, A or B.
It is possible to infer that there is a relationship among the lithological type, resilient behavior and shakedown.For basalt, after the PD test, all samples show a relative loss of stiffness in relation to the conventional RM, a result similar to that obtained by Guimarães (2009) for other basaltic material; it is also clear that this loss is not proportionally evidenced among all pairs.For those whose behavior towards shakedown pointed to behavior B, the loss was more significant.For rhyodacite, only the pairs that demonstrated accommodation of the particles had the stiffness increase of the samples with the history of stresses; the other pairs lost stiffness, as expressive as the higher the deformations undergone.Finally, for the syenogranite all samples gained stiffness after the long-term PD test, with the exception of the last pair.Close results were obtained by Lima et al. (2017) and Norback (2018) for UGM of similar origin.It is also noticed that the pairs that had behavior A towards the shakedown had an increase in RM which was more significant than the others.For all the materials, according to Figure 9, as the deviator stress increases for the same confining stress, the resilient modulus of the material decreases, in accordance to shakedown theory.A deeper investigation must be conducted in order to verify any relationship between this decrease and the type of aggregate.

Void ratio and deformability
Based on physical indexes and compaction parameters, it was observed a relation between the samples void ratio and its RM.This behavior seems to be dependent on the lithological origin of the UGM. Figure 10 exhibits the results of average RM for Compound model and void ratio; the filled markers represent the RM after PD and the not filled markers denote the conventional RM test.Contribution to resilient and permanent deformation investigation of unbound granular materials with different geological origins from Rio Grande do Sul, Brazil Basalt SJ presented the highest modulus, followed by SBS syenogranite and DPA rhyodacite, opposed to the void ratio achieved in the compaction.Based on this work, materials with different properties -such as abrasion, shape, texture and crushing process, even with the same particle size distribution, followed a trendline that suggest that void ratio may be a key parameter to predict RM.Currently, we are researching the relation between stress state and physical indexes and the densification, gain of stiffness and loss of void volume during the PD test.This research topic is still in progress and is not fully developed, so we limited the contribution in the present manuscript for this initial analysis.

Final considerations
This research aimed, by means of laboratory tests and mathematical modeling, to understand the deformability, especially resilient modulus and permanent deformation, of unbound granular materials.For this purpose, the granulometric curve was kept constant with the basalt, rhyodacite and syenogranite, materials representative of the Brazilian state of Rio Grande do Sul.
The stiffness of the materials, expressed by the resilient modulus, was tested and mathematically represented by models k-3 σ , k-d σ , k-θ, Compound, Universal, Witczak and Witczak and Uzan.The models that demonstrated the best fitting were those that consider the action of the confining and deviator stresses.In accordance to the literature, granular materials are mostly confining stress dependent; which was stated in this research as well.
The performance of the three UGMs regarding to RM and PD were different.For damage by permanent deformation, the trend found for the RM was not maintained.Syenogranite, for example, presented the worst behavior in view of permanent deformation, even though showed an intermediate RM.This effect may be related to the greater granulation of the minerals of the SBS rock in addition to the foliated structure, which starts to show poor mechanical behavior (see results for LA abrasion, impact value and crushing value).These mechanical characteristics result in the breaking of the particles promoted by the zones of weaknesses among the minerals, promoting a thinner granulometric curve, liable to high permanent deformations.
The Guimarães' model proved to be adequate to represent the PD of the materials under study, showing a high coefficient of determination.Some pairs seemed to show accommodation after 150000 cycles, and the model made it possible to predict the number of cycles needed.Furthermore, this study described a relationship between the conventional RM and after the PD test, lithological origin and shakedown.
Also, this study showed a relevant relation between the void ratio of the compacted UGM and its resilient modulus.The void ratio is affected by aggregate shape, texture, rock crushing processes and compaction parameters; this could be, with deeper investigation, a key parameter to predict RM.

List of symbols
quarry and Della Pasqua quarry, named respectively as SJ and DPA, are originated from the Planalto Meridional Province, the geographical coordinates (UTM; Universal Transverse Mercator coordinate system) of which are, respectively: 22J -251945.40m W 6826112.10 m S and 22J -228402.58m W 6724545.40 m S. In contrast, the third material, originated from SBS Engenharia quarry, called SBS, is located in the Escudo Sul-Rio-Grandense Province under coordinates 22J -357488.45m W 6483361.81 m S. Figure 1 shows the unbound granular material analyzed in this study.

Figure 1 .
Figure 1.General aspect of the unbound granular material from different lithologies analyzed in this work.

Figure 3 .
Figure 3. UGM sample being positioned on equipment (a, b) and viewed from the equipment set (c).
and 4 ψ were obtained by software Statistica v.10.0.228.2.The parameters of the Guimarães' model are essential to characterize and understand the deformability of soils and unbound granular materials, mainly due de advance of the Brazilian M-E Design Method (MeDiNa).
stress; a ρ : atmospheric pressure; N: number of loading cycles.

Figure 5 .
Figure 5. Examples of the statistical fitting -Compound (a) and Universal (b) Models for SJ.

Figure 9 .
Figure 9. Relationship between conventional RM and after PD -SJ, DPA and SBS.

Figure 10 .
Figure 10.Void ratio vs average RM for SJ, DPA, SBS granular materials.

Table 1 .
Mineralogical composition and rock characterization.

Table 2 .
Mathematical models of RM prediction.

Table 3 .
Parameters of mathematical modeling of RM.