Effect of Anisotropy on the Resilient Behaviour of a Granular Material in Low Traffic Pavement

Granular materials are often used in pavement structures. The influence of anisotropy on the mechanical behaviour of granular materials is very important. The coupled effects of water content and fine content usually lead to more complex anisotropic behaviour. With a repeated load triaxial test (RLTT), it is possible to measure the anisotropic deformation behaviour of granular materials. This article initially presents an experimental study of the resilient repeated load response of a compacted clayey natural sand with three fine contents and different water contents. Based on anisotropic behaviour, the non-linear resilient model (Boyce model) is improved by the radial anisotropy coefficient γ3 instead of the axial anisotropy coefficient γ1. The results from both approaches (γ1 and γ3) are compared with the measured volumetric and deviatoric responses. These results confirm the capacity of the improved model to capture the general trend of the experiments. Finally, finite element calculations are performed with CAST3M in order to validate the improvement of the modified Boyce model (from γ1 to γ3). The modelling results indicate that the modified Boyce model with γ3 is more widely available in different water contents and different fine contents for this granular material. Besides, based on the results, the coupled effects of water content and fine content on the deflection of the structures can also be observed.


Introduction
Low traffic pavements with a thin bituminous surfacing, granular base, and sub-base layers represent approximately 60% of the road network in France.
In these pavements, granular layers play an important role in the overall performance of the structure, especially in deformation behaviours, which are characterized by a long-term mechanical process and a short-term mechanical process. The short-term behaviour concerns the instantaneous deformation of the granular layers, as excited during the repeated traffic loading. The long-term behaviour concerns the accumulation of permanent deformations, which is not studied in this work.
In pavements, the granular layers are compacted to achieve maximum dry density. During the compaction process, these layers become anisotropic due to the vertical loads applied. As a result, the granular layers become more rigid vertically than horizontally.
Several studies [9][10][11][12][13][14][15][16][17] have been carried out to study anisotropy in granular materials: Casagrande and Carillo [9] have made a distinction for granular materials between initial anisotropy and the anisotropy induced by applied stress. The latter was defined as the ratio between axial and radial deformation under isotropic loading [12]. The initial anisotropy is mainly due to the shape and disposition of the grains, and the geometrical arrangement. Karasahin and Dawson [17] stated that it was possible to measure the anisotropy of resilient deformation of unbound granular material with a repeated load triaxial test.
Oda and Sudoo [13] indicated that anisotropy induced by applied stress was mainly due to the development of plastic deformation. The plastic deformation, meantime, can be influenced by the coupled effects of water content and fine content or suction [18]. As a result, anisotropy induced in granular layers becomes more complex.
Besides, the granular materials usually exhibit a complex non-linear response under repeated loading. A reliable description of resilient behaviour using experimental methods with or without soil suction control is time-consuming, and needs sophisticated equipment as well as trained personnel [19]. Hence, in recent years, various elasticity models have been proposed to predict the resilient behaviour based on resilient modulus and Poisson's ratio [20][21][22] or volumetric and shear stress-strain relationships [23,24]. A comparison of the most widely used elasticity models in pavement design can be found in the following Ref. [25].
Among these models, the Boyce model [23] presents the advantage of being simple (three parameters), and correctly simulating the effect of the stress on the resilient behaviour (dilatancy for high values of ∆q/∆p). As mentioned above, the anisotropy of the granular materials is increasingly being recognized as a property that must be modelled to adequately describe the pavement's mechanical behaviour. Hence, Hornych et al. [24] improved the Boyce model to include the anisotropy, and the new model showed a good agreement between measured and modelling results. However, the effects of water content and fine content on the anisotropy were not considered in this improvement.
In this context, this paper initially presents the resilient behaviour of a clayey sand (Missillac sand) with different water contents and three fine contents through the RLTT results. The results allow us to improve the Boyce model by considering the effects of water content and fine content on the anisotropic behaviour of the granular materials. A finite element modelling will be additionally performed by taking into account the proposed approach to calculate the deflection of a low traffic pavement constituted of granular materials with three different fine contents at different water contents.

Material Studied
The studied material is Missillac fine sand. It is an alluvium sand that comes from the quarry of Missillac in France. This soil is used as subgrade soil in low traffic pavements for full-scale pavement tests at IFSTTAR (Institut Français des Sciences et Technologies des Transports, de l'Aménagement et des Réseaux) in Nantes, France.
In this study, the natural samples have a maximum particle size of approximately 10 mm and contain three different fine contents (4.0%, 7.5%, and 15.3%, passing through the sieve 80 µm), which are named M4.0, M7.5, and M15.3. Figure 1 shows the particle size distribution for the three materials. Table 1 presents all of the characteristic parameters of these curves and the methylene blue values. Soils can also be classified based on these parameters [26,27], as reported in Table 1.  The standard Proctor compaction curves ( Figure 2) show the optimum water content as between 9.1-9.5% and a maximum dry density of approximately 2.04 Mg/m 3 for the three studied materials [28].

Repeated Load Triaxial Tests (RLTT)
Repeated load triaxial tests are performed on all three materials (M4.0, M7.5, and M15.3) with the water contents ranging between 7% and 11%, as shown in Table 2. The samples are compacted to the same dry density of 2 ± 0.06 Mg/m 3 with two different methods (vibrocompression method and vibrating hammer method), as shown in Table 2. The size of the sample is large enough to avoid the particle size effect: it has a diameter of 150 mm/160 mm, and a height of 285 mm/320 mm (Table 2).
In this study, the RLTTs are conducted with a variable confining pressure (VCP), where both the axial load and the confining pressure are cycled in phase [29]. As a result, the sample under the cyclic stress (confining pressure σ3 and axial stress σ1) can be used to simulate the pavement solicitation and measure the vertical strain (ε1) and the radial strain (ε3). To determine the resilient behaviour of the materials, the sample is first subjected to a conditioning phase that consists of applying 10 4 loading cycles to stabilize the permanent deformation: The standard Proctor compaction curves ( Figure 2) show the optimum water content as between 9.1-9.5% and a maximum dry density of approximately 2.04 Mg/m 3 for the three studied materials [28]. The standard Proctor compaction curves ( Figure 2) show the optimum water content as between 9.1-9.5% and a maximum dry density of approximately 2.04 Mg/m 3 for the three studied materials [28].

Repeated Load Triaxial Tests (RLTT)
Repeated load triaxial tests are performed on all three materials (M4.0, M7.5, and M15.3) with the water contents ranging between 7% and 11%, as shown in Table 2. The samples are compacted to the same dry density of 2 ± 0.06 Mg/m 3 with two different methods (vibrocompression method and vibrating hammer method), as shown in Table 2. The size of the sample is large enough to avoid the particle size effect: it has a diameter of 150 mm/160 mm, and a height of 285 mm/320 mm (Table 2).
In this study, the RLTTs are conducted with a variable confining pressure (VCP), where both the axial load and the confining pressure are cycled in phase [29]. As a result, the sample under the cyclic stress (confining pressure σ3 and axial stress σ1) can be used to simulate the pavement solicitation and measure the vertical strain (ε1) and the radial strain (ε3). To determine the resilient behaviour of the materials, the sample is first subjected to a conditioning phase that consists of applying 10 4 loading cycles to stabilize the permanent deformation:

Repeated Load Triaxial Tests (RLTT)
Repeated load triaxial tests are performed on all three materials (M4.0, M7.5, and M15.3) with the water contents ranging between 7% and 11%, as shown in Table 2. The samples are compacted to the same dry density of 2 ± 0.06 Mg/m 3 with two different methods (vibrocompression method and vibrating hammer method), as shown in Table 2. The size of the sample is large enough to avoid the particle size effect: it has a diameter of 150 mm/160 mm, and a height of 285 mm/320 mm ( Table 2).
In this study, the RLTTs are conducted with a variable confining pressure (VCP), where both the axial load and the confining pressure are cycled in phase [29]. As a result, the sample under the cyclic stress (confining pressure σ 3 and axial stress σ 1 ) can be used to simulate the pavement solicitation and measure the vertical strain (ε 1 ) and the radial strain (ε 3 ). To determine the resilient behaviour of the materials, the sample is first subjected to a conditioning phase that consists of applying 10 4 loading cycles to stabilize the permanent deformation:
p is the mean normal stress and q is the deviatoric stress. For a triaxial test, the p and q are defined as: where σ 1 and σ 3 are the major and minor principal stresses, respectively. In this context, the permanent axial deformation achieves the equilibrium state after 10 4 loading cycles. As a result, the study of resilient behaviour can be continued: after the conditioning phase, the samples (M4.0, M7.5, and M15.3) are subjected to a series of loads with different stress paths of ∆q/∆p = 0; 0.5; 1; 2; 3 [29,30], as shown in Table 3. Each stress path contains 100 loading and unloading cycles. The last cycle in each load sequence is used to determine the resilient behaviour. The soil resilient behaviour is generally presented by the volumetric strain (ε v ) and the shear or deviatoric strain (ε q ) varying with the mean normal stress (p). These deformations can be related to the axial strain (ε 1 ) and the radial strain (ε 3 ):     Based on the test results, it can be stated that the resilient behaviour of the soil is obviously nonlinear, and depends on the mean stress p and the stress path ∆q/∆p. The coupled effects of water content and fine content on resilient behaviour are evident: when the water content is higher, the εv r and the εq r are then higher in each stress path for each material. However, the effect of fine content on resilient deformation is related to the soil saturation state. Besides, large increments of the values of εv r and εq r and open loops can be observed for M15.3 samples at water contents of 11.1% in the stress paths of ∆q/∆p = 0; 0.5; and 1. (The same phenomenon is obtained in M15.3 samples at water contents of 10.2%, 11.0%, and 11.3%). Based on the test results, it can be stated that the resilient behaviour of the soil is obviously non-linear, and depends on the mean stress p and the stress path ∆q/∆p. The coupled effects of water content and fine content on resilient behaviour are evident: when the water content is higher, the ε v r and the ε q r are then higher in each stress path for each material. However, the effect of fine content on resilient deformation is related to the soil saturation state. Besides, large increments of the values of ε v r and ε q r and open loops can be observed for M15.3 samples at water contents of 11.1% in the stress paths of ∆q/∆p = 0; 0.5; and 1. (The same phenomenon is obtained in M15.3 samples at water contents of 10.2%, 11.0%, and 11.3%). Figure 5 presents the final deformations ε 1 and ε 3 of M15.3 samples with different water contents after an isotropic consolidation phase (∆q/∆p = 0, p = 10 kPa) before conditioning. It can be stated that there is an inherent anisotropy in these samples caused by compaction. The difference between ε 1 and ε 3 obviously decreases as water content for the samples increases, with a high degree of saturation (w >10.5%). The results of M4.0 and M7.5 are lost in the tests.  Figure 5 presents the final deformations ε1 and ε3 of M15.3 samples with different water contents after an isotropic consolidation phase (∆q/∆p = 0, p = 10 kPa) before conditioning. It can be stated that there is an inherent anisotropy in these samples caused by compaction. The difference between ε1 and ε3 obviously decreases as water content for the samples increases, with a high degree of saturation (w >10.5%). The results of M4.0 and M7.5 are lost in the tests.  Comparing the ratio value of ε1 to ε3 in Figure 5, and the ratio value of ε1,max r to ε3,max r in Figure  6b, it can be stated that the anisotropic behaviour of the granular material is very complicated and is strongly influenced by the stress history and the stress level applied. In the meantime, the comparison of the ratio value of ε1,max r to ε3,max r within Figure 6a,b shows that the anisotropic behaviour of M15.3 increases faster than M4.0 with an increase of water content from 7% to 11%.

Anisotropy
Based on these results, it can be stated that the obvious inherent and developing anisotropic behaviour during loading phases shouldn't be ignored.   Figure 5 presents the final deformations ε1 and ε3 of M15.3 samples with different water contents after an isotropic consolidation phase (∆q/∆p = 0, p = 10 kPa) before conditioning. It can be stated that there is an inherent anisotropy in these samples caused by compaction. The difference between ε1 and ε3 obviously decreases as water content for the samples increases, with a high degree of saturation (w >10.5%). The results of M4.0 and M7.5 are lost in the tests.  Figure 6 shows the maximum resilient deformations ε1,max r and ε3,max r for Missillac sand M4.0 and M15.3 in the isotropic loading stress path (∆q/∆p = 0, ∆p = 80 kPa) for the resilient phase. The anisotropy induced can be also observed obviously: the difference between ε1,max r and ε3,max r increases with water content for both materials, especially for M15.3. Comparing the ratio value of ε1 to ε3 in Figure 5, and the ratio value of ε1,max r to ε3,max r in Figure  6b, it can be stated that the anisotropic behaviour of the granular material is very complicated and is strongly influenced by the stress history and the stress level applied. In the meantime, the comparison of the ratio value of ε1,max r to ε3,max r within Figure 6a,b shows that the anisotropic behaviour of M15.3 increases faster than M4.0 with an increase of water content from 7% to 11%.
Based on these results, it can be stated that the obvious inherent and developing anisotropic behaviour during loading phases shouldn't be ignored. Comparing the ratio value of ε 1 to ε 3 in Figure 5, and the ratio value of ε 1,max r to ε 3,max r in Figure 6b, it can be stated that the anisotropic behaviour of the granular material is very complicated and is strongly influenced by the stress history and the stress level applied. In the meantime, the comparison of the ratio value of ε 1,max r to ε 3,max r within Figure 6a,b shows that the anisotropic behaviour of M15.3 increases faster than M4.0 with an increase of water content from 7% to 11%. Based on these results, it can be stated that the obvious inherent and developing anisotropic behaviour during loading phases shouldn't be ignored.

Modelling
In this section, the modified Boyce model based on anisotropy coefficients γ 1 and γ 3 will be compared in order to describe the resilient behaviour of Missillac sand.

Modified Boyce Model
Boyce [23] proposed an isotropic non-linear model for predicting the resilient volumetric deformation and resilient deviatoric deformation of granular materials, which was expressed as: where β = (1 − n)·K a /(6G a ); K a and G a are the reference bulk and shear moduli, respectively; and n is the nonlinear coefficient. Hornych et al. [24] introduced the anisotropic response of granular materials into the Boyce model through multiplying the axial stress by an anisotropy coefficient γ 1 . Thus, the new mean normal stress p* and the new deviator stress q* could be re-expressed as follows: where γ 1 is the axial anisotropy coefficient. Then, the potential function can be modified as: The new equations of resilient volumetric and deviatoric deformation could be proposed by taking the derivative of the potential function as: ε r q = 2 3 · p * n p a n−1 In fact, the anisotropy coefficient γ 1 is used to decrease the axial stress σ 1 to fit with the radial stiffness in the elastic range. However, as presented in Figures 5 and 6, the anisotropy behaviour changes sharply with water content for samples with a high fine content where the radial deformation ε 3 is dominant, and the two-time relation for ε 3 in Equation (3), should be also taken into account. As a result, we think that γ 1 could be no longer applicable for these results, and the radial anisotropy coefficient γ 3 is recommended with a broader scope. In other words, with the large difference between ε 1 and ε 3 , it is supposed that there is no longer such a relationship where γ 3 = 1/γ 1 , and the model for γ 3 is more general.
Thus, p* and q* can be redefined as follows: Similar to the above mentioned, the new resilient volumetric deformation ε v r and new resilient deviatoric deformation ε q r with anisotropy coefficient γ 3 can be also proposed by taking the derivative of the potential function as: ε r q = 2 3 · p * n p a n−1

Modelling Based on Two Approaches (γ 1 and γ 3 )
In this section, the modified Boyce model based on anisotropy coefficients γ 1 and γ 3 will be applied to describe the resilient behaviour of Missillac sand to demonstrate the improvement with γ 3 .

•
Parameter optimization for stress path of ∆q/ ∆p = 0 The resilient deformations of M15.3 samples at high water contents in the stress path of ∆q/∆p = 0 are obviously larger than other slopes (0.5, 1, 2 and 3), as shown in Figures 3f and 4f. To obtain the good correlation coefficients for these samples, a different anisotropy coefficient γ 3 * is used for the stress path of ∆q/∆p = 0 with the approach γ 3 for all three materials. In the meantime, to keep things consistent, γ 1 * is also required for the stress path of ∆q/ ∆p = 0 with the approach γ 1 .
The model parameters, as well as the correlation coefficients, are presented in Table 4 for M4.0 samples, Table 5 for M7.5 samples, and Table 6 for M15.3 samples with different water contents.  Based on the results, it can be stated that almost the same correlation coefficient values (close or higher than 0.8) can be found for M4.0 samples with both methods (Table 4). The better correlation coefficient values (higher than 0.9) are observed and almost the same for M7.5 samples with both methods (Table 5). However, a significant improvement of correlation coefficients for M15.3 samples at the high water contents (10.2%, 11.0%, 11.1%, and 11.3%) can be observed from γ 1 (γ 1 *) to γ 3 (γ 3 *) ( Table 6).
Besides, it can be noticed that the γ 1 and γ 1 * values, as well as the γ 3 and γ 3 * values, are almost constant or scatter slightly for M4.0 and M7.5 with all of the water contents (from 7.0% to 11.0%), as shown in Tables 4 and 5. For the M15.3 samples, the γ 1 and γ 1 * values, decrease and the γ 3 and γ 3 * values obviously increase as the water content increases, as shown in Table 6. Figures 7 and 8 show the example of modelling the ε v r and ε q r values for Missillac sand M15.3 at water content of w = 11.1% based on γ 1 (γ 1 *) and γ 3 (γ 3 *), respectively. As described by the correlation coefficients above, an obvious improvement of the estimated results (both ε v r and ε q r ) can be obtained.
Consequently, it can be stated that the radial anisotropy coefficient γ 3 , compared with the axial anisotropy coefficient γ 1 , can improve the prediction of the modified Boyce model for the granular materials with a high water content and a high fine content.
Consequently, it can be stated that the radial anisotropy coefficient γ3, compared with the axial anisotropy coefficient γ1, can improve the prediction of the modified Boyce model for the granular materials with a high water content and a high fine content.

Relationship between γ1 and γ3 (γ1* and γ3*)
The values of γ3 and 1/γ1 and the values of γ3* and 1/γ1* are presented in Figure 9 for M4.0, M7.5, and M15.3 samples. Generally, the results show that the relations of γ3 = 1/ γ1 and γ3* = 1/γ1* can be obtained approximately in all of the samples of M4.0 and M7.5. For M15.3 samples, the relation of γ3 = 1/ γ1 fits well, except for samples at water contents of 11.1% and 11.3%. The relation of γ3* = 1/γ1* fits only with samples in two low water contents of 8.1% and 8.6%. Consequently, it can be stated that the radial anisotropy coefficient γ3, compared with the axial anisotropy coefficient γ1, can improve the prediction of the modified Boyce model for the granular materials with a high water content and a high fine content.

Relationship between γ 1 and γ 3 (γ 1 * and γ 3 *)
The values of γ 3 and 1/γ 1 and the values of γ 3 * and 1/γ 1 * are presented in Figure 9 for M4.0, M7.5, and M15.3 samples. Generally, the results show that the relations of γ 3 = 1/ γ 1 and γ 3 * = 1/γ 1 * can be obtained approximately in all of the samples of M4.0 and M7.5. For M15.3 samples, the relation of γ 3 = 1/ γ 1 fits well, except for samples at water contents of 11.1% and 11.3%. The relation of γ 3 * = 1/γ 1 * fits only with samples in two low water contents of 8.1% and 8.6%.  The results confirm that the reciprocal relationship between the axial anisotropy coefficient and the radial anisotropy coefficient can change with the ratio of ε1 r to ε3 r after a threshold. The increase of this ratio or anisotropic behaviour is caused by the coupled effects of water content and fine content in this study. In fact, there are some other influences, such as the shape and disposition of the grains or the loading history, that can induce the anisotropic behaviour of granular materials. As a result, for more general behaviours for granular materials, it is necessary to take into account the variation of anisotropic behaviour as well.

Numerical Modelling of a Flexible Pavement Structure
The objective of this section is to compare the effects of two anisotropy coefficients, γ1 and γ3, on the deflection of a low traffic pavement structure. Finite elements modelling is performed with the finite elements code CAST3M.
The pavement structure [33] consists of a 66mm of asphalt surface layer placed on 500 mm of unbound granular base over 2.22 m of clayey sand (Missillac) subgrade (Figure 10a). Due to the symmetry, the 3D calculation is carried out on a quarter of the structure (Figure 10b). The pavement is discretized into 2200 cubical elements with 20 nodes. The applied load is a 65 kN dual wheel load, corresponding to the standard axle load used in France for pavement design. The loading area is rectangular, with a length of 0.30 m and a width of 0.18 m. The gravity and lateral stresses are first applied to the pavement structure to establish the initial in situ stress states.
The following materials are used for the calculation: • The asphalt concrete is used in the surface of pavement. The material behavior is assumed to be linear elastic. The Young elastic modulus (E) is equal to 6110 MPa, and its Poisson ratio is taken as equal to 0.35.

•
For the subgrade layer, the Missillac sands (M4.0, M7.5, and M15.3) are used with different water contents (from 7% to 11%). The elastic behaviour is described by the modified Boyce model with anisotropy coefficients γ1 and γ3, respectively. The model parameters in Tables 4-6 are used for The results confirm that the reciprocal relationship between the axial anisotropy coefficient and the radial anisotropy coefficient can change with the ratio of ε 1 r to ε 3 r after a threshold. The increase of this ratio or anisotropic behaviour is caused by the coupled effects of water content and fine content in this study. In fact, there are some other influences, such as the shape and disposition of the grains or the loading history, that can induce the anisotropic behaviour of granular materials. As a result, for more general behaviours for granular materials, it is necessary to take into account the variation of anisotropic behaviour as well.

Numerical Modelling of a Flexible Pavement Structure
The objective of this section is to compare the effects of two anisotropy coefficients, γ 1 and γ 3 , on the deflection of a low traffic pavement structure. Finite elements modelling is performed with the finite elements code CAST3M.
The pavement structure [33] consists of a 66mm of asphalt surface layer placed on 500 mm of unbound granular base over 2.22 m of clayey sand (Missillac) subgrade (Figure 10a). Due to the symmetry, the 3D calculation is carried out on a quarter of the structure (Figure 10b). The pavement is discretized into 2200 cubical elements with 20 nodes. The applied load is a 65 kN dual wheel load, corresponding to the standard axle load used in France for pavement design. The loading area is rectangular, with a length of 0.30 m and a width of 0.18 m. The gravity and lateral stresses are first applied to the pavement structure to establish the initial in situ stress states.
The following materials are used for the calculation: • The asphalt concrete is used in the surface of pavement. The material behavior is assumed to be linear elastic. The Young elastic modulus (E) is equal to 6110 MPa, and its Poisson ratio is taken as equal to 0.35.

•
For the subgrade layer, the Missillac sands (M4.0, M7.5, and M15.3) are used with different water contents (from 7% to 11%). The elastic behaviour is described by the modified Boyce model with anisotropy coefficients γ 1 and γ 3 , respectively. The model parameters in Tables 4-6 Figure 11 shows the deflections of the structures calculated by two different methods (γ1 and γ3) with three Missillac sands (M4.0, M7.5, and M15.3) at different water contents. It can be observed that the deflection values are identical with the two methods (γ1 and γ3) for all three materials at different water contents, except for the structures with M15.3 at water contents (10.2% and 11.1%, where γ1 ≠ 1/γ3). The difference of deflection between the two methods at a water content of 10.2% can be explained by the poor calculation with γ1 (Correl = 0.631). For the water content of 11.1%, the convergence cannot even reach with γ1. The simulation results show that the new proposed approach with anisotropy coefficient γ3 for the granular material has improved the modelling, especially for the case with large resilient deformation (at high water contents and high fine contents). Besides, the coupled effect of water content and fine content on the deflection of the structures can be observed (Figure 11b): the deflections of the pavement will be obviously higher when the fine content increases from 4% to 7.5%. For the pavement with a high fine content of 15.3%, the matric suction is significantly higher in the unsaturated part because of more fine particles, which leads to a lower resilient deformation. As for the saturated part, the suction is very small, and more fine  . The difference of deflection between the two methods at a water content of 10.2% can be explained by the poor calculation with γ 1 (Correl = 0.631). For the water content of 11.1%, the convergence cannot even reach with γ 1 . The simulation results show that the new proposed approach with anisotropy coefficient γ 3 for the granular material has improved the modelling, especially for the case with large resilient deformation (at high water contents and high fine contents).  Figure 11 shows the deflections of the structures calculated by two different methods (γ1 and γ3) with three Missillac sands (M4.0, M7.5, and M15.3) at different water contents. It can be observed that the deflection values are identical with the two methods (γ1 and γ3) for all three materials at different water contents, except for the structures with M15.3 at water contents (10.2% and 11.1%, where γ1 ≠ 1/γ3). The difference of deflection between the two methods at a water content of 10.2% can be explained by the poor calculation with γ1 (Correl = 0.631). For the water content of 11.1%, the convergence cannot even reach with γ1. The simulation results show that the new proposed approach with anisotropy coefficient γ3 for the granular material has improved the modelling, especially for the case with large resilient deformation (at high water contents and high fine contents). Besides, the coupled effect of water content and fine content on the deflection of the structures can be observed (Figure 11b): the deflections of the pavement will be obviously higher when the fine content increases from 4% to 7.5%. For the pavement with a high fine content of 15.3%, the matric suction is significantly higher in the unsaturated part because of more fine particles, which leads to a lower resilient deformation. As for the saturated part, the suction is very small, and more fine Besides, the coupled effect of water content and fine content on the deflection of the structures can be observed (Figure 11b): the deflections of the pavement will be obviously higher when the fine content increases from 4% to 7.5%. For the pavement with a high fine content of 15.3%, the matric suction is significantly higher in the unsaturated part because of more fine particles, which leads to a lower resilient deformation. As for the saturated part, the suction is very small, and more fine particles lead to decreasing the resistance of the resilient deformation because of the combination with water. Figure 12 also shows the deflections of the structures with three fine contents at the water content of 11% using the new proposed approach (γ 3 parameter). It can be seen that the evolution of deflection varies along the lateral position, and the effect of fine content is important.

Conclusions
This article initially presents the resilient behaviour of a natural compacted sand (Missillac sand), which is studied with RLTTs at different water contents and five stress paths varying from zero to three. The inherent anisotropy caused by compaction, and the anisotropy induced by loading stress, can be observed. The anisotropy is strongly influenced by the coupled effects of water content and fine content, which has a large increase as the water content increases (close to saturated state) for a high fine content sample.
Next, the radial anisotropy coefficient is suggested to replace the axial anisotropy coefficient in the non-linear model proposed by Boyce (1980), in order to determine the resilient behaviour of Missillac sand at three different fine contents and different water contents from 7% to 11%. The modelling of the results indicates that the radial anisotropy coefficient γ3, compared with the axial anisotropy coefficient γ3, show more general results, and thus obviously improves the accuracy of the modelling of this granular material with a high fine content close to saturated state. It also confirms the capacity of the improved model to capture the general trend of the experiments.
Finally, the finite elements modelling of a flexible pavement is performed with CAST3M to validate the improvement of the modified Boyce model (from γ1 to γ3). The modelling results present that the modified Boyce model with γ3 is more widely available in different water contents and different fine contents for this granular material. Besides, based on the results, the coupled effects of water content and fine content on the deflection of the structures can be observed apparently.
Consequently, the influence of the anisotropy on the resilient behaviour of granular materials has been seldom taken into account in modelling work. Furthermore, there is a strong link between the anisotropy behaviour and the coupled effects of water content and fine content. Understanding the complex mechanism of the anisotropic resilient deformation is helpful to propose more general resilient models or improve the existing models. This suggests that the number of tests used to predict the resilient behaviour for different moisture conditions or different fine contents can be reduced.
For the future work, this finding will be verified with different granular materials and even higher fine contents.
Acknowledgments: This work is supported by the China Scholarship Council and French national project Durée

Conclusions
This article initially presents the resilient behaviour of a natural compacted sand (Missillac sand), which is studied with RLTTs at different water contents and five stress paths varying from zero to three. The inherent anisotropy caused by compaction, and the anisotropy induced by loading stress, can be observed. The anisotropy is strongly influenced by the coupled effects of water content and fine content, which has a large increase as the water content increases (close to saturated state) for a high fine content sample.
Next, the radial anisotropy coefficient is suggested to replace the axial anisotropy coefficient in the non-linear model proposed by Boyce (1980), in order to determine the resilient behaviour of Missillac sand at three different fine contents and different water contents from 7% to 11%. The modelling of the results indicates that the radial anisotropy coefficient γ 3 , compared with the axial anisotropy coefficient γ 3 , show more general results, and thus obviously improves the accuracy of the modelling of this granular material with a high fine content close to saturated state. It also confirms the capacity of the improved model to capture the general trend of the experiments.
Finally, the finite elements modelling of a flexible pavement is performed with CAST3M to validate the improvement of the modified Boyce model (from γ 1 to γ 3 ). The modelling results present that the modified Boyce model with γ 3 is more widely available in different water contents and different fine contents for this granular material. Besides, based on the results, the coupled effects of water content and fine content on the deflection of the structures can be observed apparently.
Consequently, the influence of the anisotropy on the resilient behaviour of granular materials has been seldom taken into account in modelling work. Furthermore, there is a strong link between the anisotropy behaviour and the coupled effects of water content and fine content. Understanding the complex mechanism of the anisotropic resilient deformation is helpful to propose more general resilient models or improve the existing models. This suggests that the number of tests used to predict the resilient behaviour for different moisture conditions or different fine contents can be reduced.
For the future work, this finding will be verified with different granular materials and even higher fine contents.