Cyclic True Triaxial Tests on Aeolian Sand Considering Initial Shear Effect

: Under trafﬁc loading, the soil elements in subgrade are subjected to a complex 3D stress path. To investigate the cyclic behavior of desert subgrade under initial shear stress conditions, the number of cyclic true triaxial tests were implemented on aeolian sand from the Tengger Desert. A large range of initial shear stress levels and different cyclic stress paths (various combinations of cyclic major and intermediate principal stresses) were designed in the experiments. The results show that the initial shear stress level signiﬁcantly inﬂuences the cyclic response, and the response mode of aeolian sand under initial shear stress is incremental failure and elastic shakedown, while it is critical failure without initial shear stress. With the increase in initial shear stress, the permanent strain increases ﬁrst and then decreases, with the maximum permanent strain occurring at q 0 = 50 kPa. Moreover, by comparing the test data under different cyclic stress paths, it is found that the cyclic resistance decreases with the increase in the coefﬁcient of cyclic intermediate principal stress b cyc . Compared to the limiting pore pressure criterion, the conventional 5% axial strain failure criterion may overestimate the cyclic resistance, leading to unsafe evaluation and design. Therefore, by using the pore pressure criterion, the relationship between the limit pore pressure ratio and the initial stress condition was investigated and it was found that the limit pore pressure ratio decreased linearly when the initial shear stress increased.


Introduction
In practical engineering, soil is usually in a complex three-dimensional static and dynamic state, so it is of great practical significance to study the mechanical behavior of soil under three-dimensional stress state.Slope, foundation pit, and subgrade in the construction stage are mainly under static loading, while the long-term use and maintenance should consider the combining action of dynamic and static loads, especially for the subgrade under long-term traffic loading.For safety design, a full understanding of the static liquefaction, cyclic liquefaction, and shakedown of sandy soils are very important.Those phenomena have been studied through geotechnical investigations, model tests, and laboratory experiments [1][2][3].Conventional triaxial tests are used extensively to measure the static and cyclic characteristics of soil because of the convenient operation [4][5][6].However, in the conventional triaxial test, the cylindrical specimen is under an axisymmetric stress state, as a static test can only control the principal stress in two directions independently, and a dynamic test can only apply one or two cyclic loads, which is quite different from the actual three-dimensional stress state [7,8].Therefore, some researchers began to carry out cyclic true triaxial tests [9,10], and the deformation characteristics of soil under two cyclic principal stresses were studied.
During the process of natural deposition, gravitational force leads to anisotropic consolidation and results in initial shear stress q 0 .The cyclic triaxial tests were conducted on anisotropic consolidated samples to measure the dynamic behavior of sand under nonzero q 0 , wherein the cyclic stress was superimposed on the initial shear stress.Tests results illustrate that there is a significant difference in the cyclic response of sand between nonzero and zero initial shear stress [11][12][13].Lee and Seed [14] implemented a large number of cyclic triaxial tests, and they indicate that the existence of q 0 may improve the liquefaction resistance of sand.However, Ma [15] shows that the effect of q 0 on liquefaction resistance depends on the level of q 0 , and the liquefaction resistance is increased by a lower q 0 , while it is decreased by a higher q 0 .Vaid and Chern [16,17] also obtain similar conclusions, and it is shown that the resistance of liquefaction might show an increase or decrease trend with the existence of q 0 , which depends on the relative density of the sand.A possible explanation for the above conflicting observations by Harder and Boulanger [18] is that the q 0 plays different roles of in loose sand and dense sand.For dense sand, the q 0 may improve its liquefaction resistance, while for loose sand, it becomes damaged when the q 0 is higher.Recently, a unified relationship between state parameter ψ [19] and liquefaction resistance of Toyoura sand was established by Yang and Sze [20,21], which considered the influence of initial density and confining pressure.However, most of the existing studies use axisymmetric triaxial tests or hollow torsional shear tests to investigate the cyclic behaviors of sand under initial shear stress conditions, where the cyclic loading is confined to one-dimensional loading and the test range is limited.
The true triaxial apparatus can load independently on three orthogonal planes of a cubic soil specimen, and the static and dynamic characteristics of soil can be measured in true three-dimensional space [22][23][24].Three orthogonal principal stress loading proportions are controlled by the coefficient of intermediate principal stress (b = (σ 1 − σ 3 )/(σ 2 − σ 3 )) [25,26].Based on a series of true triaxial tests, Li et al. [25,26] studied the mechanical behavior of saturated aeolian sand under true triaxial conditions, and found that the b had a significant influence on stress-strain, pore pressure, and critical state.Shi et al. [22] investigated the influence of b on the stress-strain relationship and strength behavior of coarse-grained soil, and a shape function g(b) on π-plane was introduced.However, the cyclic tests under true triaxial conditions are rarely studied, only Wang et al. [10] implemented the number of axial and lateral cyclic triaxial tests to study the deformation characteristics of gravel under different lateral cyclic conditions.The results suggest that the intermediate principal stress has a significant influence on the strain behavior of the specimen.
As the subgrade soils are subjected to large numbers of load (10 4 -10 5 cycles) at a stress level lower than its shear strength [27], the deformation characteristics of soil under long-term cyclic loading are particularly important.The shakedown theory was proposed by Melan [28] firstly, and it has been extensively used to describe the cyclic response of subgrade soils under traffic loading.For pavement design [29][30][31] and railway track foundation design [32][33][34], the critical dynamic stress is used as the design load to distinguish the stable state from the unstable state of the pavement or track foundation.Mamou et al. [32][33][34] observed that once the critical dynamic stress was exceeded, the rate of stiffness degradation, excess pore pressure generation, and axial strain accumulation grew significantly, leading to the eventual failure of the soil.Based on a series of cyclic triaxial tests under different stress paths, the permanent deformation model with the loading cycles can be classified into three types: plastic shakedown, plastic creep, and incremental failure, according to the value of ε p 5000 -ε p 3000 that are the permanent strain obtained at the 5000th and 3000th loading cycles during tests.Nevertheless, some researchers [35][36][37] state that this method is problematic because the classification results based on this method do not match the real test results.
Moreover, a large number of tests found that the characteristic of pore pressure in the specimens could not reach the initial mean stress, and there may be a limiting value far below the initial mean stress.Therefore, the effective stress can never reach zero [38,39].Hyodo et al. [40] reveal that the limiting value of pore water pressure decreases with the increase in initial shear stress.It is also found that the direction of initial shear stress also influences the development of pore pressure.To quantitatively describe the development of pore pressure under cyclic loading, some methods have been put forward on the basis of cyclic triaxial tests.Based on Seed's pioneering work [41], an empirical model was developed and the relationship between cycle number and pore pressure was derived.The model has been used extensively for the prediction of the pore pressure development with the number of loading cycles [42,43].For example, based on the above model, Pan et al. [44][45][46] proposed a modified pore pressure accumulation model, which can be used to describe the development of pore pressure under different initial shear stress conditions.Although there are many models that can be applied to predict the development of cyclic pore pressure, the following discussion should be considered.All the above studies are based on the cyclic triaxial test under an axisymmetric state, without considering the influence of three-dimensional cyclic stress path on pore pressure.Recent work by Zhang et al. [9] observed that the limiting values of pore pressure are highly related to the level of cyclic intermediate principal stress.
To investigate the cyclic behavior of desert subgrade under initial shear stress conditions, a large number of cyclic true triaxial tests were implemented on aeolian sand from the Tengger Desert.At the same time, the construction of the Wuma highway was used as the background, which is located at the southern edge of the Tengger Desert.The effect of the interaction between the two cyclic principal stresses on the cyclic behavior of anisotropic consolidated specimens was taken into account.The differences in strain and pore pressure development under different initial shear stress levels were emphasized, in which the factors of cyclic shear stress amplitude q ampl , and the coefficient of cyclic intermediate principal stress b cyc were considered.The applicability of strain failure criterion and limiting pore pressure criterion under three-dimensional stress conditions were compared.

Test Apparatus and Specimen Preparation
The automated true triaxial apparatus developed by GDS Corporation in England was employed in this study.It can perform many functions, including isotropic and anisotropic consolidation, as well as various loading modes (see [47] for details).The most prominent merits of this device are: (1) the cyclic major principal stress (vertical direction) and the cyclic intermediate principal stress (horizontal direction) can be applied simultaneously; (2) the rigid-flexible mixed boundary control is adopted, the major and intermediate principal stresses are applied by metal plates, and the minor principal stress σ 3 is commonly provided by confining pressure, loaded with water (flexible loading).
Aeolian sand was used in all tests, which is a uniform fine sand taken from a depth of 1 m below the southern edge of the Tengger Desert [48].The primary properties of the aeolian sand are shown in Table 1.Dry deposition method was used to prepare the cubic specimens of 75 mm × 75 mm × 150 mm, which can simulate the formation process of aeolian sand and make the specimens have an initial fabric closer to the field environment (see [47] for details).The aeolian sand specimens used in this test were all medium-dense sand with a relative density of 50%.The specimens were saturated with back pressure.All the specimens were completely saturated in this study, and the Skempton's B value was higher than 0.96.

Test Program
Isotropic consolidation was performed on saturated specimens first, and then the initial shear stress was applied.The samples were anisotropically consolidated by applying initial shear stress while keeping σ 2 = σ 3 , considering that the depositional plane of specimens is transverse isotropy.The initial shear stress q 0 can be calculated as follows: Figure 1 shows the consolidation and cyclic stress paths of the specimen in the cyclic true triaxial test.In a cyclic test, sine wave load with constant amplitude is applied, that is, only compressional stress without extensional stress, and there is no phase difference between the cyclic major and intermediate principal stresses.Different planar cyclic stress paths are applied by controlling the proportion of cyclic major principal stress ∆σ 1 and intermediate principal stress ∆σ 2 .A parameter named the coefficient of cyclic intermediate principal stress b cyc is used to represent the interaction between ∆σ 1 and ∆σ 2 , as follows: where , and σ ampl 3 are the amplitudes of cyclic major, intermediate, and minor principal stresses, respectively.The minor principal stress σ 3 is kept constant during the cycles (i.e., σ ampl 3 = 0) because considering that σ ampl 3 is small enough, it can be neglected in the stress field induced by traffic loading [27].Therefore, the slope of the stress path in the σ 1 -σ 2 plane is exactly the value of b cyc .Under three-dimensional stress state, the cyclic stress state of soil is described by the cyclic shear stress amplitude q ampl , which is expressed as follows.
Sustainability 2022, 14, 16730 5 of 14 For comparison, the initial effective confining pressure of all specimens were 100 kPa.In cyclic shear, the number of cycles for cyclic loading was set to 10,000 (the test terminates prematurely if the specimens fail before that).The frequency was 1 Hz, which was regarded as representative for real traffic [49].Considering the high loading speed of traffic loading, all tests were carried out under undrained condition.Table 2 summarizes For comparison, the initial effective confining pressure of all specimens were 100 kPa.In cyclic shear, the number of cycles for cyclic loading was set to 10,000 (the test terminates prematurely if the specimens fail before that).The frequency was 1 Hz, which was regarded as representative for real traffic [49].Considering the high loading speed of traffic loading, all tests were carried out under undrained condition.Table 2 summarizes all of the test conditions, which include 36 cyclic true triaxial tests under various combinations of q 0 , b cyc, and q ampl .Four values of q 0 (0, 50, 100, 150 kPa), which cover both isotropically and anisotropically consolidated samples, were considered in this study.Three values of b cyc (0, 0.4, and 0.8) were employed to represent different cyclic stress paths.When b cyc = 0, it is actually the one-dimensional cyclic loading that is only loaded in the vertical direction of the specimen deposition plane.By comparing the experimental results of b cyc = 0.4 and b cyc = 0.8, the different cyclic behavior of planar and one-way cyclic loading can be measured.q ampl values ranging from 20 kPa to 40 kPa were considered for studying the cyclic behavior of aeolian sand in different cyclic response modes.

Strains Development under One-Dimensional Cyclic Loading
Figure 2 shows the typical strain development under one-dimensional cyclic loading, including q 0 = 0 kPa, lower q 0, and higher q 0 .In each figure, the time-history curves of axial strains (ε 1 ) and horizontal strains (ε 2 ) are presented.Three different cyclic response modes can be observed in the test: critical failure, incremental failure, and elastic shakedown.

Strains Development under One-Dimensional Cyclic Loading
Figure 2 shows the typical strain development under one-dimensional cyclic loading, including q0 = 0 kPa, lower q0, and higher q0.In each figure, the time-history curves of axial strains (ε1) and horizontal strains (ε2) are presented.Three different cyclic response modes can be observed in the test: critical failure, incremental failure, and elastic shakedown.Figure 2a shows the critical failure for the isotropically consolidated specimen.It can be seen that the strain starts to develop slowly under cyclic loading.The increase in strain is accelerated when it comes to 400 cycles, and, thereafter, tends to approach a large value, resulting in the failure of the specimen.Figure 2b illustrates the typical results of the aeolian sand with a small q0.In this figure, a typical behavior of the incremental failure can be observed, which is different from the critical failure shown in Figure 2a.The strain of the specimen develops rapidly at the beginning of loading, and the rate of strain  Figure 2a shows the critical failure for the isotropically consolidated specimen.It can be seen that the strain starts to develop slowly under cyclic loading.The increase in strain is accelerated when it comes to 400 cycles, and, thereafter, tends to approach a large value, resulting in the failure of the specimen.Figure 2b illustrates the typical results of the aeolian sand with a small q 0 .In this figure, a typical behavior of the incremental failure can be observed, which is different from the critical failure shown in Figure 2a.The strain of the specimen develops rapidly at the beginning of loading, and the rate of strain accumulation decreases gradually with the cyclic loading, as observed in Figure 2b. Figure 2c shows the elastic shakedown of the aeolian sand with a large initial shear stress.In the figure, an abrupt increase in the strain is evident in the first cycles, while, in the following cycles, the strain does not rise.The plastic deformation is no longer developed (ε p 5000 -ε p 3000 < 1%), and the whole specimen presents elastic characteristics.Although the experimental phenomena in Figure 2b, and c are similar, significantly differences still exist.For example, in the beginning cycles, the rate of strain increase in Figure 2c is faster than that in Figure 2b, whereas the rate of strain increase is almost 0 in Figure 2c at the end of loading, while in Figure 2b the strain is still increasing at a slow rate.
To make a comparison among different kinds of cyclic behavior about the cyclic resistance, Hyodo et al. [50] and Ishihara [51] proposed the failure criterion based on strain.This approach uses the 5% single-amplitude axial strain to determine the liquefaction or failure state of sand.In Figure 2, the number of cycles when the axial strain reaches 5% is also shown.It should be noted that for the elastic shakedown response mode, the deformation of the specimen tends to be stable after 1000 cycles, and the failure criterion of axial strain 5% is not reached after 10,000 cycles, as shown in Figure 2c.
In order to analyze the effect of q 0 on the deformation of specimen, Figure 3 depicts the variation of axial and horizontal strains with loading cycles under different q 0 in the b cyc = 0 test.In general, the q 0 has a significant effect on the deformation mode of the specimen, and with the increase in q 0 , the ultimate deformation mode of the specimen is elastic shakedown.
mode, the deformation of the specimen tends to be stable after 1000 cycles, and the failure criterion of axial strain 5% is not reached after 10,000 cycles, as shown in Figure 2c.
In order to analyze the effect of q0 on the deformation of specimen, Figure 3 depicts the variation of axial and horizontal strains with loading cycles under different q0 in the bcyc = 0 test.In general, the q0 has a significant effect on the deformation mode of the specimen, and with the increase in q0, the ultimate deformation mode of the specimen is elastic shakedown.Figure 3a shows a comparison between different q0 under small cyclic stress (q ampl = 20 kPa) about the cyclic response of the specimen.The strain that occurs in all tests is less than 0.1%, except for the test under q0 = 150 kPa.For the test of q ampl = 30 kPa, the deformation of the specimen is the largest in the q0 = 50 kPa test and the smallest in the q0 = 150 kPa test, as shown in Figure 3b. Figure 3c shows the cyclic response under q ampl = 40 kPa.With the increase in q0, the cyclic response of the specimen is close to critical failure, incremental failure, and elastic stability.The maximum deformation is observed in the q0 = 50 kPa test.

Strains Development under Two-Dimensional Cyclic Loading
This part mainly illustrates the test results of strains under the interaction of two cyclic principal stresses.Figure 4 shows time-history curves of axial and horizontal strain under different bcyc.Overall, with the increase in bcyc, the cyclic resistance is decreasing.Under the condition of q0 = 0 kPa, the failure mode of aeolian sand is critical failure, and under the condition of non-zero q0, the cyclic response is incremental failure and elastic shakedown.
Horizontal Strain Number Figure 3a shows a comparison between different q 0 under small cyclic stress (q ampl = 20 kPa) about the cyclic response of the specimen.The strain that occurs in all tests is less than 0.1%, except for the test under q 0 = 150 kPa.For the test of q ampl = 30 kPa, the deformation of the specimen is the largest in the q 0 = 50 kPa test and the smallest in the q 0 = 150 kPa test, as shown in Figure 3b. Figure 3c shows the cyclic response under q ampl = 40 kPa.With the increase in q 0 , the cyclic response of the specimen is close to critical failure, incremental failure, and elastic stability.The maximum deformation is observed in the q 0 = 50 kPa test.

Strains Development under Two-Dimensional Cyclic Loading
This part mainly illustrates the test results of strains under the interaction of two cyclic principal stresses.Figure 4 shows time-history curves of axial and horizontal strain under different b cyc .Overall, with the increase in b cyc , the cyclic resistance is decreasing.Under the condition of q 0 = 0 kPa, the failure mode of aeolian sand is critical failure, and under the condition of non-zero q 0 , the cyclic response is incremental failure and elastic shakedown.
The interaction of these two cyclic principal stresses has a significant effect on the development of strains, and the effect is related to both q ampl and q 0 .For the test of q 0 = 0 kPa, as b cyc increases, the number of cycles in which it reaches critical failure increases, which means that its cyclic resistance decreases, as shown in Figure 4a.For the test of q 0 = 50 kPa, the cyclic response for all tests is incremental failure, that is, the strain increases rapidly during the early stages of cycling, which is shown in Figure 4b.When b cyc = 0.8, because of the influence of the cyclic intermediate principal stress, a large axial strain is produced in the beginning cycles, resulting in a significant degradation of stiffness and the specimens reach failure in 400 cycles.For the test of q 0 50 kPa, the cyclic response of all specimens is elastic shakedown, and the values of strain components are obviously decreased, as shown in Figure 4c,d.In the case of b cyc = 0.8, at N = 5000, the axial strain when q 0 = 100 kPa is approximately 1.5 times that of q 0 = 150 kPa (7.1% vs. 4.6%).Furthermore, the existence of q 0 would change the direction of horizontal strain.For example, in the test of b cyc = 0.4 and 0.8 in Figure 4a, the horizontal strain is positive, which is the opposite of the characteristics reflected in the other figure.
To show the effects of q 0 on permanent deformation of specimen more clearly, the ε 1 p at N = 100 and 1000 under different q 0 are replotted in Figure 5.It is clear that the ε 1 p is the largest in the q 0 = 50 kPa tests, which corresponds to the b cyc = 0 tests.It can be concluded that the lower q 0 could increase the deformation of the specimen while the higher q 0 could reduce the deformation.In addition, Figure 5 shows that ε 1 p tends to increase with increasing q ampl , whereas q ampl has little influence on the trend of ε 1 p -q 0 curves.This is inconsistent with the conclusions obtained in previous studies [15][16][17].This indicates that saturated sand has different properties under true triaxial conditions.The interaction of these two cyclic principal stresses has a significant effect on the development of strains, and the effect is related to both q ampl and q0.For the test of q0 = 0 kPa, as bcyc increases, the number of cycles in which it reaches critical failure increases, which means that its cyclic resistance decreases, as shown in Figure 4a.For the test of q0 = 50 kPa, the cyclic response for all tests is incremental failure, that is, the strain increases rapidly during the early stages of cycling, which is shown in Figure 4b.When bcyc = 0.8, because of the influence of the cyclic intermediate principal stress, a large axial strain is produced in the beginning cycles, resulting in a significant degradation of stiffness and the specimens reach failure in 400 cycles.For the test of q0 ˃ 50 kPa, the cyclic response of all specimens is elastic shakedown, and the values of strain components are obviously decreased, as shown in Figure 4c,d.In the case of bcyc = 0.8, at N = 5000, the axial strain when q0 = 100 kPa is approximately 1.5 times that of q0 = 150 kPa (7.1% vs. 4.6%).Furthermore, the existence of q0 would change the direction of horizontal strain.For example, in the test of bcyc = 0.4 and 0.8 in Figure 4a, the horizontal strain is positive, which is the opposite of the characteristics reflected in the other figure.
To show the effects of q0 on permanent deformation of specimen more clearly, the ε1 p at N = 100 and 1000 under different q0 are replotted in Figure 5.It is clear that the ε1 p is the largest in the q0 = 50 kPa tests, which corresponds to the bcyc = 0 tests.It can be concluded that the lower q0 could increase the deformation of the specimen while the higher q0 could reduce the deformation.In addition, Figure 5 shows that ε1 p tends to increase with increasing q ampl , whereas q ampl has little influence on the trend of ε1 p -q0 curves.This is inconsistent with the conclusions obtained in previous studies [15][16][17].This indicates that saturated sand has different properties under true triaxial conditions.
(a) (b) Figure 6 compares ε1 p -ε2 p curves under different level of q0.The slope of the ε1 p -ε2 p curves indicates the direction of its permanent strain path.The figure shows that the ε1 p -ε2 p curves have a similar trend under different q0, with the increase in q0, the length of ε1 p -ε2 p curves decreases and the direction of ε1 p -ε2 p curves rotates counterclockwise.For the test of bcyc = 0, all of the ε1 p -ε2 p curves are below horizon, which means that ε2 p is always q ampl =20kPa q ampl =30kPa q ampl =40kPa ε p 1 (%) at N=1000 q 0 (kPa) q ampl =20kPa q ampl =30kPa q ampl =40kPa ε p 1 (%) at N=100 q 0 (kPa)  curves have a similar trend under different q 0 , with the increase in q 0 , the length of ε 1 p -ε 2 p curves decreases and the direction of ε 1 p -ε 2 p curves rotates counterclockwise.For the test of b cyc = 0, all of the ε 1 p -ε 2 p curves are below horizon, which means that ε 2 p is always negative.For the test of b cyc = 0.4 and 0.8, the ε 1 p -ε 2 p curves are above the horizon line only at q 0 = 0 kPa.Moreover, the slope of the ε 1 p -ε 2 p curves increases with the increase in b cyc , implying that the increase in b cyc leads to the increase in ε 2 p .
( Figure 6 compares ε1 p -ε2 p curves under different level of q0.The slope of the ε1 p -ε2 curves indicates the direction of its permanent strain path.The figure shows that the ε1 p ε2 p curves have a similar trend under different q0, with the increase in q0, the length of ε1 p ε2 p curves decreases and the direction of ε1 p -ε2 p curves rotates counterclockwise.For th test of bcyc = 0, all of the ε1 p -ε2 p curves are below horizon, which means that ε2 p is alway negative.For the test of bcyc = 0.4 and 0.8, the ε1 p -ε2 p curves are above the horizon line only at q0 = 0 kPa.Moreover, the slope of the ε1 p -ε2 p curves increases with the increase in bcy implying that the increase in bcyc leads to the increase in ε2 p .The above results show that the cyclic response of saturated aeolian sand depend on q ampl , bcyc, and q0. Figure 7 summarizes the cyclic response of the tests under differen stress conditions.In this figure, in addition to the three response modes for large deformations, there is a stable mode with minimal permanent strain (<0.1%).It can be seen from the figure that the cyclic response of all specimens is elastic shakedown type at q0 b cyc =0.4,q ampl =40kPa q 0 =0kPa q 0 =50kPa q 0 =100kPa q 0 =150kPa Permanent horizontal strain The above results show that the cyclic response of saturated aeolian sand depends on q ampl , b cyc, and q 0 .Figure 7 summarizes the cyclic response of the tests under different stress conditions.In this figure, in addition to the three response modes for larger deformations, there is a stable mode with minimal permanent strain (<0.1%).It can be seen from the figure that the cyclic response of all specimens is elastic shakedown type at q 0 = 150 kPa.Under the condition of q 0 = 50 kPa, the cyclic response of specimens is mainly incremental failure.In particular, all the critical failures occur at q 0 = 0 kPa.Sustainability 2022, 14, 16730 10 of 14 150 kPa.Under the condition of q0 = 50 kPa, the cyclic response of specimens is mainly incremental failure.In particular, all the critical failures occur at q0 = 0 kPa.

The Development of Pore Pressure
Under the dynamic loading, the excess pore pressure of saturated soils can be classified into two parts: transient and residual pore pressure.The transient pore pressure is the real-time pore pressure during loading and unloading, which has a relatively weaker effect on the variation of effective stress in soil.However, the residual pore

The Development of Pore Pressure
Under the dynamic loading, the excess pore pressure of saturated soils can be classified into two parts: transient and residual pore pressure.The transient pore pressure is the real-time pore pressure during loading and unloading, which has a relatively weaker effect on the variation of effective stress in soil.However, the residual pore pressure directly affects the effective stress in soil, which is the accumulated pore pressure at the end of each cycle.
Typical developments of excess pore pressure ∆u in undrained tests are presented in Figure 8. Figure 8a shows the development of ∆u without initial shear stress, the cyclic response is critical failure, and the strain development is shown in Figure 2a.It can be seen that ∆u rapidly accumulates during the beginning cycles, and ∆u reaches the limiting value of 86 kPa at 513 cycles.This is the exact opposite of the characteristic of strain development shown in Figure 2a. Figure 8b shows the development of ∆u in the existence of initial shear stress.As can be seen that ∆u rises vertically at the beginning of loading and reaches 90% of the limit pore pressure in 20 cycles.Notably, ∆u reaches the limiting value of 63 kPa at 41 cycles, whereas the strain begins to stabilize after 1000 cycles (Figure 2b).

The Development of Pore Pressure
Under the dynamic loading, the excess pore pressure of saturated soils can be classified into two parts: transient and residual pore pressure.The transient pore pressure is the real-time pore pressure during loading and unloading, which has a relatively weaker effect on the variation of effective stress in soil.However, the residual pore pressure directly affects the effective stress in soil, which is the accumulated pore pressure at the end of each cycle.
Typical developments of excess pore pressure Δu in undrained tests are presented in Figure 8. Figure 8a shows the development of Δu without initial shear stress, the cyclic response is critical failure, and the strain development is shown in Figure 2a.It can be seen that Δu rapidly accumulates during the beginning cycles, and Δu reaches the limiting value of 86 kPa at 513 cycles.This is the exact opposite of the characteristic of strain development shown in Figure 2a. Figure 8b shows the development of Δu in the existence of initial shear stress.As can be seen that Δu rises vertically at the beginning of loading and reaches 90% of the limit pore pressure in 20 cycles.Notably, Δu reaches the limiting value of 63 kPa at 41 cycles, whereas the strain begins to stabilize after 1000 cycles (Figure 2b). Figure 9 shows the relationship between residual pore pressure permanent deformation under different stress conditions, where residual pore pressure ratio r u is defined as the ratio of residual pore pressure to initial effective confining pressure.Figure 9a presents the results without initial shear stress with b cyc = 0.4.As the number of cycles increases, r u and ε 1 p increase gradually.Furthermore, r u reaches the limiting value of r u lim = 0.89 at N = 2766, with a permanent strain of approximately 2.67%.With the loading process, r u remains constant, while ε 1 p continues to increase until reaching 5% at 3099 cycles, that is, the strain-based failure criterion.Figure 9b shows the results with initial shear stress (q 0 = 50 kPa), r u and ε 1 p increase sharply in the incipient cycles, and the rate of increase in r u is significantly faster than that in Figure 9a.The limiting value of r u lim = 0.75 is approached at N = 222 with ε 1 p = 4.72%.In subsequent cycles, r u remains constant, while ε 1 p continues to increase moderately until the failure criterion is reached at 376 cycles.
process, r u remains constant, while ε1 p continues to increase until reaching 5% at 3099 cycles, that is, the strain-based failure criterion.Figure 9b shows the results with initial shear stress (q0 = 50 kPa), r u and ε1 p increase sharply in the incipient cycles, and the rate of increase in r u is significantly faster than that in Figure 9a.The limiting value of r u lim = 0.75 is approached at N = 222 with ε1 p = 4.72%.In subsequent cycles, r u remains constant, while ε1 p continues to increase moderately until the failure criterion is reached at 376 cycles.Obviously, the constant accumulation of permanent strain is associated with the development of residual pore pressure, and the limit value of residual pore pressure ratio depends on the initial shear stress.Figure 10 shows the relation between the limiting value of r u lim and q0 at q ampl = 40 kPa.It can be seen that r u lim decreases linearly with the increase in q0 and reaches its maximum at q0 = 0 kPa, indicating that the residual pore water pressure of saturated aeolian sand under isotropic consolidation is closest to the initial effective confining pressure, as critical failure dominates.In addition, it can be seen that r u lim increases with the increase in bcyc.Obviously, the constant accumulation of permanent strain is associated with the development of residual pore pressure, and the limit value of residual pore pressure ratio depends on the initial shear stress.Figure 10 shows the relation between the limiting value of r u lim and q 0 at q ampl = 40 kPa.It can be seen that r u lim decreases linearly with the increase in q 0 and reaches its maximum at q 0 = 0 kPa, indicating that the residual pore water pressure of saturated aeolian sand under isotropic consolidation is closest to the initial effective confining pressure, as critical failure dominates.In addition, it can be seen that r u lim increases with the increase in b cyc .
increases, r u and ε1 p increase gradually.Furthermore, r u reaches the limiting value of r u lim = 0.89 at N = 2766, with a permanent strain of approximately 2.67%.With the loading process, r u remains constant, while ε1 p continues to increase until reaching 5% at 3099 cycles, that is, the strain-based failure criterion.Figure 9b shows the results with initial shear stress (q0 = 50 kPa), r u and ε1 p increase sharply in the incipient cycles, and the rate of increase in r u is significantly faster than that in Figure 9a.The limiting value of r u lim = 0.75 is approached at N = 222 with ε1 p = 4.72%.In subsequent cycles, r u remains constant, while ε1 p continues to increase moderately until the failure criterion is reached at 376 cycles.
(a) (b) Obviously, the constant accumulation of permanent strain is associated with the development of residual pore pressure, and the limit value of residual pore pressure ratio depends on the initial shear stress.Figure 10 shows the relation between the limiting value of r u lim and q0 at q ampl = 40 kPa.It can be seen that r u lim decreases linearly with the increase in q0 and reaches its maximum at q0 = 0 kPa, indicating that the residual pore water pressure of saturated aeolian sand under isotropic consolidation is closest to the initial effective confining pressure, as critical failure dominates.In addition, it can be seen that r u lim increases with the increase in bcyc.Figure 11 make a comparison between the required cycles for reaching limiting pore pressure (N lim ) and that for permanent axial strain reaching 5% (N f ).When the specimens show critical failure behavior, the test data all fall on the diagonal line of the image, which indicates that N lim (pore pressure criterion) and N f (strain criterion) are close.However, in the response to incremental failure, the N lim is generally smaller than N f , which is mainly because the axial permanent deformation continues to increase after the pore pressure reaches its limit.The scenario indicates that the pore pressure criterion can be used to estimate the failure of saturated aeolian sand in a safe way.
pressure (Nlim) and that for permanent axial strain reaching 5% (Nf).When the specimens show critical failure behavior, the test data all fall on the diagonal line of the image, which indicates that Nlim (pore pressure criterion) and Nf (strain criterion) are close.However, in the response to incremental failure, the Nlim is generally smaller than Nf, which is mainly because the axial permanent deformation continues to increase after the pore pressure reaches its limit.The scenario indicates that the pore pressure criterion can be used to estimate the failure of saturated aeolian sand in a safe way.

Conclusions
This paper presents a series of cyclic true triaxial tests at various levels of initial shear stress.The combined effect of the two-dimensional cyclic stress and initial shear stress on the cyclic behavior of saturated aeolian sand is investigated.The characteristics of strain development and pore pressure are analyzed, and the experimental results of different initial stress states are compared.The main conclusions can be given as follows: The cyclic response of saturated aeolian sand in a three-dimensional stress state can be identified in three distinct types: critical failure, incremental failure, and elastic shakedown.The critical failure dominates for specimens under q0 = 0 kPa conditions, while the incremental failure and elastic shakedown are observed with the existence of a non-zero q0.
With the increase in q0, the permanent strain increases first for lower q0 and then decreases, and the maximum permanent strain occurs at q0 = 50 kPa.In addition, the increase in bcyc and q ampl lead to the growth of permanent strain.
The r u lim decreases linearly with the increase in q0, and increases moderately with the increase in bcyc.Under the critical failure response, the value of r u lim is the highest.
The pore pressure criterion is conservative in estimating the failure of saturated aeolian sand compared to the strain failure criterion (5%).The results from these two failure criteria are well-matched when q0 = 0 kPa, but they deviate when initial shear stress is in existence.

Figure 1 .
Figure 1.Stress path in cyclic true triaxial test.

Figure 1 .
Figure 1.Stress path in cyclic true triaxial test.

Figure 7 .
Figure 7. Summary of cyclic response modes under different stress conditions.

Figure 7 .
Figure 7. Summary of cyclic response modes under different stress conditions.

Figure 7 .
Figure 7. Summary of cyclic response modes under different stress conditions.

Figure 10 .
Figure 10.Relationship between limiting pore pressure ratio and initial shear stress.

Figure 11 .
Figure 11.The number of cycles to reach pore pressure (Nlim) and axial permanent strain of 5% (Nf).

Table 1 .
Primary properties of aeolian sand.

Table 2 .
Scheme of cyclic tests.