Influence of Load Frequency on Cohesive Soil Respond

The mechanism of cohesive soils response to cycling loading is less investigated compared to cohesionless soils. Multiple load-unload cycles cause significant changes in the structure of cohesive soils, which result in complex behaviour under the given load. The aim of the paper was to investigate and study the influence of load frequency on cohesive soil reaction. In order to obtain results, tests were conducted using the cyclic triaxial apparatus. Three cyclic tests were carried out, each for different frequency −0.5 Hz, 1.0 Hz, 2.0 Hz and one static triaxial test. The maximal value of deviator stress qmax, used in the cyclic tests, was set to 40 kPa. Afterwards samples were unloaded to qmin = 30 kPa. Cyclic loading triaxial tests were performed in a consolidated-undrained (CU) one-way loading manner, a sinusoidal waves were used. After the cycling loading was completed, a static triaxial shear test was conducted. Changes in the cohesive soil responses depending on cycling load frequency were presented in the paper. Differences in the accumulation of plastic strains were noticed, as well as changes of degradation index values, resilient degradation index values and differences in the excess pore water pressure development.


Introduction
Analysis of the influence of cyclic load acting on a soil is a common geotechnical problem nowadays.The increasing economy creates the demand for a new transportation network.In the last decades, the repeating loading problem was limited to the issue which includes the studies on dynamic load bearing capacity of foundations, the reaction of machine foundations subjected to dynamic loads, the interaction between soil and surrounding constructions loaded by earthquake-generated seismic waves and the strength of dams and embankments during an earthquake [1][2][3].Most of the technical design standards do not include the rules concerning the behaviour of cohesive soils subjected to cyclic load [4].Multiple load-unload cycles cause significant changes in soil structure, which affects its behaviour [5,6].The most unpleasant occurrence of cyclic loading impact on soil subgrade is long-term settlement [7,8].Natural sources (water waves, earthquakes, storms, etc.) and human activity (machinery vibration, wind power plants, construction operations, traffic loads and many others) are known to be the sources of this particular type of load.Proper identification of load characteristics has a special importance when reconstructing the load conditions in the laboratory [8].
In this paper, the attention is focused on so-called 'quasi-static' loadings.This term describes the repeating loading in which the inertia forces are not considered or are negligible.The conditions of quasi-static loading must not exceed a certain value of frequency and amplitude.The quasi-static loading occurs when the acceleration of the loaded soil particles is significantly lower than the gravitational acceleration [9].
The behaviour of cohesive soils subjected to the cyclic load has not been sufficiently recognized and described [10][11][12].So far, focus of the researches starts to turn on the problem of cohesive soils under cyclic loads, when compared to non-cohesive soils, for which the literature is known to be more extensive [13][14][15][16][17][18].Due to its nature, the cohesive soil may be subjected to many load-unload cycles before the destruction of soil occurs [19,20].However, in cohesionless soils, less resistance to cyclic loads is observed, which may lead to liquefaction of the soil [21,22].
One of the interesting topics in the subject of cohesive soils cyclic loading is cyclical degradation of clays.This particular subject has been studied by Idriss et al. and others for over 40 years [19,20,[23][24][25][26].The cyclic degradation is an important problem since the results of this type of loading can be dangerous for human life.One of the examples was given by Yamanouchi and Yasuhara [27] where the case history of highway embankment was presented.The embankment was founded on soft cohesive soil, namely alluvial Ariake Clay and the subgrade exhibited large settlement after being opened do traffic.The settlement was equal to nearly 2 m and was observed on the length of 1700 m.The soft soil term in this article refers to a clayey soils with high plasticity and low stiffness which is caused usually by high water content.
The cohesive soil respond to cyclic loading depends on many factors, among which most important are confining pressure, void ratio, strain amplitude, temperature, overconsolidation ratio and anisotropic stress [28].What is more, the plasticity index and permeability of cohesive soils, which may lead undrained conditions, have a significant impact.
One of the factors that has been given less attention is frequency of loading.The triaxial tests, which were conducted on cohesive soils, led so far to a few conclusions.Generally, there exists the impact of the loading frequency on soil respond.Nevertheless, this impact tends to vanish with the number of cycles.The same effect was observed in case of deviator stress amplitude.Another interesting conclusion concerning frequency impact is when the soil is loaded at lower frequencies, the sample maintain higher cyclic strength.What is more, the tests have shown that the lower frequency is applied the larger shear strains are accumulated and the larger pore pressure is generated [25,[28][29][30][31][32].
The cyclic triaxial tests performed by Lei et al. [20] on Tianjin Binhai clay covered investigation on the frequency and cyclic stress ratio impact on the creep behaviour of soil subjected to cyclic loading.The results of the tests lead to the conclusion that cyclic loads produce greater compaction than static loads and the Tianjin Binhai clay creeps more under the repeating loading compared to static constant stress.The analysis of frequency impact has shown a presence of a specific frequency and a safe load under which, the strain of the soft clay increases very slowly.It was found that greater frequency produces more rapid soil failure.
This paper presents results of cyclic triaxial tests where impact of cyclic loading on cohesive soil response was studied.The objective of this study was to investigate the frequency impact on the strength degradation of cohesive soil samples.The primary goal was to determine the cohesive soil responses to repeating loading in terms of deformation which was described by Idriss [26] conception of degradation modulus.The secondary goal was to determine how ultimate strength after cyclic loading in different frequencies is changing.The knowledge about these effects is important for evaluation of cohesive soils strength properties for analyses of the cyclic responses of the soil deposits.

Material Characterization
The soil specimens used in this article were cohesive soils taken from the road construction project sites in Warsaw, Poland.The cohesive soil was deposited in a glacial erosion conditions.To characterize the sample, a number of test were conducted that have been described below.
In order to determine granular composition, the aerometric method was applied.Based on granular structure, determined by the aerometric method, the soil was classified as sasiCl (sandy silty clay, according to PN-EN ISO 14688:2006 [33] standard, based on international standards).The results of the gradation analysis tests are presented Figure 1.
In order to determine granular composition, the aerometric method was applied.Based on granular structure, determined by the aerometric method, the soil was classified as sasiCl (sandy silty clay, according to PN-EN ISO 14688:2006 [33] standard, based on international standards).The results of the gradation analysis tests are presented Figure 1.In the next step, the Atterberg limits were examined.Plastic limit was determined by rolling a portion of soil into a ball with 7 mm diameter.Afterwards, sample should be rolled into a thread with 3 mm diameter.When cracks were observed, sample's weight was tested and then sample was dried in 105 °C for 24 h to determine moisture content.Examined soil sample's plastic limit (PL) equals 13.71%.Liquid limit (LL) was examined using Casagrande cup and equals around 26.0% (Figure 2).The plasticity index (PI) of this soil is equal to 12.29%.Optimal moisture content and maximum dry intensity was defined in the Proctor compaction test, which has been chosen due to typical engineering application preservation.Nine samples were prepared and examined with different moisture content.Compaction in this particular test requires 0.59 J of energy to be transferred to 1.0 cm 3 of soil specimen.In order to prepare a proper sample for triaxial test, a part of soil was compacted in 5 layers, each blown 40 times with a compactor.The In the next step, the Atterberg limits were examined.Plastic limit was determined by rolling a portion of soil into a ball with 7 mm diameter.Afterwards, sample should be rolled into a thread with 3 mm diameter.When cracks were observed, sample's weight was tested and then sample was dried in 105 • C for 24 h to determine moisture content.Examined soil sample's plastic limit (PL) equals 13.71%.Liquid limit (LL) was examined using Casagrande cup and equals around 26.0% (Figure 2).The plasticity index (PI) of this soil is equal to 12.29%.
In order to determine granular composition, the aerometric method was applied.Based on granular structure, determined by the aerometric method, the soil was classified as sasiCl (sandy silty clay, according to PN-EN ISO 14688:2006 [33] standard, based on international standards).The results of the gradation analysis tests are presented Figure 1.In the next step, the Atterberg limits were examined.Plastic limit was determined by rolling a portion of soil into a ball with 7 mm diameter.Afterwards, sample should be rolled into a thread with 3 mm diameter.When cracks were observed, sample's weight was tested and then sample was dried in 105 °C for 24 h to determine moisture content.Examined soil sample's plastic limit (PL) equals 13.71%.Liquid limit (LL) was examined using Casagrande cup and equals around 26.0% (Figure 2).The plasticity index (PI) of this soil is equal to 12.29%.Optimal moisture content and maximum dry intensity was defined in the Proctor compaction test, which has been chosen due to typical engineering application preservation.Nine samples were prepared and examined with different moisture content.Compaction in this particular test requires 0.59 J of energy to be transferred to 1.0 cm 3 of soil specimen.In order to prepare a proper sample for triaxial test, a part of soil was compacted in 5 layers, each blown 40 times with a compactor.The  Optimal moisture content and maximum dry intensity was defined in the Proctor compaction test, which has been chosen due to typical engineering application preservation.Nine samples were prepared and examined with different moisture content.Compaction in this particular test requires 0.59 J of energy to be transferred to 1.0 cm 3 of soil specimen.In order to prepare a proper sample for triaxial test, a part of soil was compacted in 5 layers, each blown 40 times with a compactor.The optimum moisture content equals around 11.6%, while maximum dry density is equal to 1.97 g/cm 3 .The results of Proctor test are presented on Figure 3.
The oedometric test was also conducted in order to find the yield stress.The results presents Figure 4. Based on Casagrante's method, the yield stress or pre-consolidation stress σ' p was estimated.Three points were specified namely the minimum pre-consolidation stress σ' pMIN , the maximum pre-consolidation stress σ' pMAX and σ' p .Such procedure is important since the compacted soils have the pre-consolidation stress but the exact value is hard to find.The value of σ' pMIN is equal to 170 kPa, the value of σ' p = 270 kPa and the value of σ' p MAX is equal to 305 kPa.The conclusion can be drown, that the compacted cohesive soil in optimal moisture content with respect to Proctor's method have pre-consolidation stress and therefore, the soil tested in this article in σ' 3 = 45 kPa conditions is in pre-consolidated state.
Geosciences 2019, 9 FOR PEER REVIEW 4 of 13 optimum moisture content equals around 11.6%, while maximum dry density is equal to 1.97 g/cm 3 .The results of Proctor test are presented on Figure 3.The oedometric test was also conducted in order to find the yield stress.The results presents Figure 4. Based on Casagrante's method, the yield stress or pre-consolidation stress σ'p was estimated.Three points were specified namely the minimum pre-consolidation stress σ'pMIN, the maximum preconsolidation stress σ'pMAX and σ'p.Such procedure is important since the compacted soils have the pre-consolidation stress but the exact value is hard to find.The value of σ'pMIN is equal to 170 kPa, the value of σ'p = 270 kPa and the value of σ'p MAX is equal to 305 kPa.The conclusion can be drown, that the compacted cohesive soil in optimal moisture content with respect to Proctor's method have preconsolidation stress and therefore, the soil tested in this article in σ'3 = 45 kPa conditions is in preconsolidated state.

Triaxial Test Procedure
The tests were conducted in a cyclic triaxial apparatus for which soil specimens are required to have diameter equal 7 cm and height of 14 cm.The soil samples were preliminary compacted with respect to the Proctor method with energy of compaction equal to 0.59 J/cm 3 .Investigations were carried out with the samples fully saturated (measurement was concluded to guarantee that the degree of saturation was greater than 97% for all specimens).This type of compacted cohesive soil is optimum moisture content equals around 11.6%, while maximum dry density is equal to 1.97 g/cm 3 .The results of Proctor test are presented on Figure 3.The oedometric test was also conducted in order to find the yield stress.The results presents Figure 4. Based on Casagrante's method, the yield stress or pre-consolidation stress σ'p was estimated.Three points were specified namely the minimum pre-consolidation stress σ'pMIN, the maximum preconsolidation stress σ'pMAX and σ'p.Such procedure is important since the compacted soils have the pre-consolidation stress but the exact value is hard to find.The value of σ'pMIN is equal to 170 kPa, the value of σ'p = 270 kPa and the value of σ'p MAX is equal to 305 kPa.The conclusion can be drown, that the compacted cohesive soil in optimal moisture content with respect to Proctor's method have preconsolidation stress and therefore, the soil tested in this article in σ'3 = 45 kPa conditions is in preconsolidated state.

Triaxial Test Procedure
The tests were conducted in a cyclic triaxial apparatus for which soil specimens are required to have diameter equal 7 cm and height of 14 cm.The soil samples were preliminary compacted with respect to the Proctor method with energy of compaction equal to 0.59 J/cm 3 .Investigations were carried out with the samples fully saturated (measurement was concluded to guarantee that the degree of saturation was greater than 97% for all specimens).This type of compacted cohesive soil is

Triaxial Test Procedure
The tests were conducted in a cyclic triaxial apparatus for which soil specimens are required to have diameter equal 7 cm and height of 14 cm.The soil samples were preliminary compacted with respect to the Proctor method with energy of compaction equal to 0.59 J/cm 3 .Investigations were carried out with the samples fully saturated (measurement was concluded to guarantee that the degree of saturation was greater than 97% for all specimens).This type of compacted cohesive soil is different from the soils tested in previous works concerning soft soils.High density, low void ratio and low moisture content characterizes rather strong cohesive soils.
After the saturation stage was completed, specimens were subjected to isotropic effective confining pressure of 45 kPa, which led to specimen's consolidation.Table 1 presents the physical properties of soil the samples.The first step of the test was cyclic triaxial test.The maximal value of deviator stress q max was set to 40 kPa, afterwards samples were unloaded to q min = 30 kPa.Cyclic loading triaxial tests were performed in a consolidated-undrained (CU) one-way loading manner, a sinusoidal waves were used (average deviator stress q m was equal to around 35 kPa and deviator stress amplitude q a was equal to around 5 kPa).The load acting on a soil was applied in a vertical direction.To determine, whether the frequency of cyclic load affects the cohesive soil reaction, tests were conducted under three frequency conditions: f = 0.5 Hz, f = 1.0 Hz and f = 2.0 Hz.The number of cycles each time was N = 20,000.During the test the initial σ' 3 was equal to 45 kPa.The value of σ' 3 was set based on the test results concerning horizontal stress responses of subgrade subjected to one passage of two-axle truck.The 70 cm wheel caused in this case 690 kPa.The horizontal pressure in all three cases of the test was less than 45 kPa at the distance 300 cm from pavement edge and between 40 and 60 kPa 100 cm from the pavement edge [34].Based on this tests results we decided to set the σ' 3 value to constant 45 kPa.
The frequencies which were used for this study, are typical for common frequencies observed in human environment and in nature.For example, the typical frequencies observed in soils surrounding metro train tunnels are between 0.5 to 30 Hz.
The second step was static triaxial test which was conducted immediately after cyclic triaxial loading step.For this step additional static triaxial test was conducted in order to provide a comparison of obtained post-cyclic static test results and for calculation of cyclic stress ratio.
The strain measurements was performed with use of highly precise (1 × 10 −4 %) electro-mechanical actuator mounted to the loading beam.

The Results and Discussion
The triaxial test results analysis are divided into two subsections which covers cyclic behaviour of compacted sandy silty clay, the response to cyclic loading and the post-cyclic behaviour.This division will help to describe more precisely the impact of load frequency on cohesive soil strength and deformation characteristics.

Cyclic Triaxial Loading Test Results
Three series of cyclic triaxial tests and one static triaxial test were performed to investigate the frequency impact on cohesive soil respond.Stress-controlled loading program was applied using digital control system.During the triaxial test, the data on axial vertical strain, axial stress, pore pressure and cell pressure was recorded and stored as a raw data.On Figure 5 the results of all four tests are presented.The plot shows the stress-strain characteristics.The cyclic loading in all three cases lead to accumulation of plastic strains.Largest vertical plastic strains ε p were observed in first cycle of loading.
For frequency equal to 0.5 Hz the plastic strain was 0.0038, for 1 Hz the ε p was equal to 0.0044 and for 2.0 Hz the vertical plastic strain was equal to 0.0037.Figure 6 presents characteristics of accumulated plastic strain in following cycles.
Geosciences 2019, 9 FOR PEER REVIEW 6 of 13 to 0.0044 and for 2.0 Hz the vertical plastic strain was equal to 0.0037.Figure 6 presents characteristics of accumulated plastic strain in following cycles.The cyclic loading in further cycles (N > 2 and N < ~1000) was in logarithmic manner in all three cases.Next, the intermediate state was observed where the rate of plastic strain accumulation was steadily decreasing.The last stage (N > 5000) characterizes itself with lower rate of plastic accumulation and again with logarithmic increase of this phenomena.The presented number of cycles is comparative.The slight differences between total accumulated plastic strain is due to amount of plastic strain accumulated in the first cycle.This occurrence can be explained by difference between void ratio e which also refers to lower density in case of the sample tested in 1 Hz frequency conditions.
The respond of cohesive soil to cyclic loading can be examined in reference to the degradation index.The idea of parameters describing the soil degradation which can be defined by degradation index δ was presented by Idriss et al. [26].The δ value is defined as the ratio of dynamic secant modulus in the N-th cycle to dynamic secant modulus in the first cycle [26], shown in Equation ( 1  to 0.0044 and for 2.0 Hz the vertical plastic strain was equal to 0.0037.Figure 6 presents characteristics of accumulated plastic strain in following cycles.The cyclic loading in further cycles (N > 2 and N < ~1000) was in logarithmic manner in all three cases.Next, the intermediate state was observed where the rate of plastic strain accumulation was steadily decreasing.The last stage (N > 5000) characterizes itself with lower rate of plastic accumulation and again with logarithmic increase of this phenomena.The presented number of cycles is comparative.The slight differences between total accumulated plastic strain is due to amount of plastic strain accumulated in the first cycle.This occurrence can be explained by difference between void ratio e which also refers to lower density in case of the sample tested in 1 Hz frequency conditions.
The respond of cohesive soil to cyclic loading can be examined in reference to the degradation index.The idea of parameters describing the soil degradation which can be defined by degradation index δ was presented by Idriss et al. [26].The δ value is defined as the ratio of dynamic secant modulus in the N-th cycle to dynamic secant modulus in the first cycle [26], shown in Equation ( 1 The cyclic loading in further cycles (N > 2 and N < ~1000) was in logarithmic manner in all three cases.Next, the intermediate state was observed where the rate of plastic strain accumulation was steadily decreasing.The last stage (N > 5000) characterizes itself with lower rate of plastic accumulation and again with logarithmic increase of this phenomena.The presented number of cycles is comparative.The slight differences between total accumulated plastic strain is due to amount of plastic strain accumulated in the first cycle.This occurrence can be explained by difference between void ratio e which also refers to lower density in case of the sample tested in 1 Hz frequency conditions.
The respond of cohesive soil to cyclic loading can be examined in reference to the degradation index.The idea of parameters describing the soil degradation which can be defined by degradation index δ was presented by Idriss et al. [26].The δ value is defined as the ratio of dynamic secant modulus in the N-th cycle to dynamic secant modulus in the first cycle [26], shown in Equation (1): where, σ d is maximal axial stress (in terms of triaxial test conditions the referenced stress is maximal deviator stress q max which is equal to q max = q m + q a ) in N-th cycle and ε d is axial strain at maximal stress in N-th cycle.The degradation properties of cohesive soils under repeated loading depend on axial strain amplitude.The proposed exponential model of soil properties degradation was later examined in terms of consolidation ratio [33].
Resilient modulus M r describes the resilient respond of soil under repeating loading and can be obtained from cyclic triaxial tests.During unloading phase, the elastic strain also called resilient strain ε r can be observed which shows the exact amount of recoverable deformation in the particular cycle.The resilient strain can be calculated as a difference between total strain ε T in one cycle to permanent strain ε P in one cycle ε r = ε T − ε P .The Equation ( 2) was used to calculate M r value: where σ d,N is deviator stress in N-th cycle and ε r,N is resilient strain in N-th cycle.The idea of degradation index is expanded by authors into the resilient degradation index which refers to the first resilience characteristics and is presented by Equation (3): Results of calculations of the degradation index and resilient degradation index presents Figure 7.
where, σd is maximal axial stress (in terms of triaxial test conditions the referenced stress is maximal deviator stress qmax which is equal to qmax = qm + qa) in N-th cycle and εd is axial strain at maximal stress in N-th cycle.The degradation properties of cohesive soils under repeated loading depend on axial strain amplitude.The proposed exponential model of soil properties degradation was later examined in terms of consolidation ratio [33].
Resilient modulus Mr describes the resilient respond of soil under repeating loading and can be obtained from cyclic triaxial tests.During unloading phase, the elastic strain also called resilient strain εr can be observed which shows the exact amount of recoverable deformation in the particular cycle.The resilient strain can be calculated as a difference between total strain εT in one cycle to permanent strain εP in one cycle εr = εT − εP.The Equation ( 2) was used to calculate Mr value: where σd,N is deviator stress in N-th cycle and εr,N is resilient strain in N-th cycle.The idea of degradation index is expanded by authors into the resilient degradation index which refers to the first resilience characteristics and is presented by Equation ( 3): Results of calculations of the degradation index and resilient degradation index presents Figure 7.
(a) (b) The degradation index δ(N) can be presented as the degradation parameter t which is defined as negative quotient of common logarithm from degradation index δ(N) to common logarithm of number of cycles N, presented by Equation (4) [35]: The degradation index δ(N) can be presented as the degradation parameter t which is defined as negative quotient of common logarithm from degradation index δ(N) to common logarithm of number of cycles N, presented by Equation (4) [35]: The Young's modulus for a first cycle was equal to 9.84 MPa, 8.74 MPa and 10.11 MPa respectively for frequencies 0.5 Hz, 1 Hz and 2 Hz.The degradation index in all three cases in first 120 cycles was the same.After this, the degradation of sample loaded with frequency equal to 0.5 Hz slowed down.The degradation of samples for which f was equal to 1 Hz and 2 Hz went in the same way by 5000 cycles.After that, the degradation index of sample for which frequency was equal to 1 Hz was slightly smaller than for sample tested with f = 2 Hz.This phenomena corresponds with plastic strain accumulation presented on Figure 6.The degradation parameter t for the cycles from 1 to 120 was equal to 0.06 in cycle 120.The results of calculations show different response to cyclic loading of compacted cohesive soil.The soil sample loaded with frequency equal to 0.5 Hz accumulates fewer plastic strains during the cyclic loading.The degradation index confirms this behaviour.Samples loaded with frequency greater than 0.5 Hz accumulate more plastic strains and therefore the degradation index is smaller.
The resilient degradation index analysis shows that despite the fact the soil samples tested in frequencies equal to 1 Hz and 2 Hz accumulate more plastic strains, resilient modulus after 20,000 cycles will be greater than at the beginning of the test.This means that the stiffness of the soil samples in this frequencies will behave differently than in case of test in f equal to 0.5 Hz.The resilient modulus in first cycle was equal to 71.9 MPa, 70.7 MPa and 79.0 MPa respectively.In first few hundred cycles the resilient modulus decreased.The minimal M r values were observed in cycle around 7000 for f equal to 0.5 Hz, cycle around 550 for f equal to 1 Hz and around 250 cycle for frequency 2 Hz.Based on this test results, it can be seen that the minimal resilient modulus value occurs faster with increase in frequency.The resilient degradation and stiffness degradation will be observed in greater amount in case of low frequencies.
The degradation of resilient behaviour seems to be in opposite to plastic strain accumulation phenomena.The frequency impact on soil sample respond will be observed as less stiff and with less settlement in case of lower frequencies.
In order to understand this phenomena, the pore water pressure analysis was conducted.The Young's modulus for a first cycle was equal to 9.84 MPa, 8.74 MPa and 10.11 MPa respectively for frequencies 0.5 Hz, 1 Hz and 2 Hz.The degradation index in all three cases in first 120 cycles was the same.After this, the degradation of sample loaded with frequency equal to 0.5 Hz slowed down.The degradation of samples for which f was equal to 1 Hz and 2 Hz went in the same way by 5000 cycles.After that, the degradation index of sample for which frequency was equal to 1 Hz was slightly smaller than for sample tested with f = 2 Hz.This phenomena corresponds with plastic strain accumulation presented on Figure 6.The degradation parameter t for the cycles from 1 to 120 was equal to 0.06 in cycle 120.The results of calculations show different response to cyclic loading of compacted cohesive soil.The soil sample loaded with frequency equal to 0.5 Hz accumulates fewer plastic strains during the cyclic loading.The degradation index confirms this behaviour.Samples loaded with frequency greater than 0.5 Hz accumulate more plastic strains and therefore the degradation index is smaller.
The resilient degradation index analysis shows that despite the fact the soil samples tested in frequencies equal to 1 Hz and 2 Hz accumulate more plastic strains, resilient modulus after 20,000 cycles will be greater than at the beginning of the test.This means that the stiffness of the soil samples in this frequencies will behave differently than in case of test in f equal to 0.5 Hz.The resilient modulus in first cycle was equal to 71.9 MPa, 70.7 MPa and respectively.In first few hundred cycles the resilient modulus decreased.The minimal Mr values were observed in cycle around 7000 for f equal to 0.5 Hz, cycle around 550 for f equal to 1 Hz and around 250 cycle for frequency 2 Hz.Based on this test results, it can be seen that the minimal resilient modulus value occurs faster with increase in frequency.The resilient degradation and stiffness degradation will be observed in greater amount in case of low frequencies.
The degradation of resilient behaviour seems to be in opposite to plastic strain accumulation phenomena.The frequency impact on soil sample respond will be observed as less stiff and with less settlement in case of lower frequencies.
In order to understand this phenomena, the pore water pressure analysis was conducted.The excess pore water pressure development in all three cases had the same scenario.In the first phase, the excess pore water pressure rises to maximal value Δumax which is equal to 19.5 kPa, 27.0 kPa and 23 kPa respectively.The Δumax value was reached after around 395, 205 and 630 cycles.The excess pore water pressure development in all three cases had the same scenario.In the first phase, the excess pore water pressure rises to maximal value ∆u max which is equal to 19.5 kPa, 27.0 kPa and 23 kPa respectively.The ∆u max value was reached after around 395, 205 and 630 cycles.Further cyclic loading lead to decrease of excess pore water pressure.This behaviour lasted until cycle from 3000 to 6000.Subsequently, the excess pore water pressure reached steady characteristics and the ∆u change was limited.On the Figure 8a, the ∆u value is presented versus axial vertical strain development.Shape of this characteristics are close to the characteristic which was presented on Fig, 8b.This characteristic was presented in logarithmic scale of number of cycles.This means that, the excess pore water pressure characteristics have impact on the accumulation of the plastic strains.This can also explain why lower frequencies resulted in lower accumulation of vertical plastic strains.Lower value of excess pore water pressure in case of cyclic triaxial test in frequency 0.5 Hz results in lower accumulated plastic strain because the ∆u value was lower in comparison to the tests in higher frequencies.
Excess pore water pressure generation phenomena is also connected with degradation index characteristics.The degradation of tested soil modulus occurs due to the occurrence of the plastic strains.Therefore, the higher generation of excess pore water pressure is the grater degradation will occur.
The resilient degradation index seems to be independent from above presented mechanism.Nevertheless, the increase of the resilient modulus occurs in the zone of steady excess pore water pressure and starts just before occurrence of ∆u max value in case of the tests for f equal to 1 Hz and 2 Hz.In case of cyclic triaxial test performed in frequency equal to 0.5 Hz the increase of resilient degradation index and therefore the increase of resilient modulus starts when third stage of ∆u characteristics occurs.The conclusion can be drown that, the steady value of excess pore water pressure indicates reverse of the soil respond to cyclic loading and after the phase of softening, the hardening will occur.This also means that long term cyclic loading can improve mechanical behaviour of cohesive soils if the additional settlement is accepted.

Static Triaxial Loading Test Results
Static triaxial test was performed after the cyclic loading and one additional test was conducted for soil sample without performing the cyclic loading program.The results of the static test are presented on Figure 5.The static loading was conducted with constant velocity of displacement equal 1 mm/h.The static loading stage started just after the cyclic loading test.The static triaxial test on normally consolidated sample without pre-cyclic loading shows typical for this type of soil stress-strain characteristics.The maximal deviator stress q max during this test was hard to find due to constant increase of the stress strain characteristic.In order to find the maximal deviator stress the deviator stress-mean effective stress plot is presented on Figure 9.
Geosciences 2019, 9 FOR PEER REVIEW 9 of 13 Further cyclic loading lead to decrease of excess pore water pressure.This behaviour lasted until cycle from 3000 to 6000.Subsequently, the excess pore water pressure reached steady characteristics and the Δu change was limited.On the Figure 8a, the Δu value is presented versus axial vertical strain development.Shape of this characteristics are close to the characteristic which was presented on Fig, 8b.This characteristic was presented in logarithmic scale of number of cycles.This means that, the excess pore water pressure characteristics have impact on the accumulation of the plastic strains.This can also explain why lower frequencies resulted in lower accumulation of vertical plastic strains.Lower value of excess pore water pressure in case of cyclic triaxial test in frequency 0.5 Hz results in lower accumulated plastic strain because the Δu value was lower in comparison to the tests in higher frequencies.
Excess pore water pressure generation phenomena is also connected with degradation index characteristics.The degradation of tested soil modulus occurs due to the occurrence of the plastic strains.Therefore, the higher generation of excess pore water pressure is the grater degradation will occur.
The resilient degradation index seems to be independent from above presented mechanism.Nevertheless, the increase of the resilient modulus occurs in the zone of steady excess pore water pressure and starts just before occurrence of Δumax value in case of the tests for f equal to 1 Hz and 2 Hz.In case of cyclic triaxial test performed in frequency equal to 0.5 Hz the increase of resilient degradation index and therefore the increase of resilient modulus starts when third stage of Δu characteristics occurs.The conclusion can be drown that, the steady value of excess pore water pressure indicates reverse of the soil respond to cyclic loading and after the phase of softening, the hardening will occur.This also means that long term cyclic loading can improve mechanical behaviour of cohesive soils if the additional settlement is accepted.

Static Triaxial Loading Test Results
Static triaxial test was performed after the cyclic loading and one additional test was conducted for soil sample without performing the cyclic loading program.The results of the static test are presented on Figure 5.The static loading was conducted with constant velocity of displacement equal 1 mm/h.The static loading stage started just after the cyclic loading test.The static triaxial test on normally consolidated sample without pre-cyclic loading shows typical for this type of soil stressstrain characteristics.The maximal deviator stress qmax during this test was hard to find due to constant increase of the stress strain characteristic.In order to find the maximal deviator stress the deviator stress-mean effective stress plot is presented on Figure 9.The maximal deviator stress at which the failure occurs q max,f , based on the Figure 9 for static triaxial test, is equal to around 45 kPa.Therefore the cyclic stress ratio (CSR) which is defined as (5): where, the q max is maximal deviator stress during the cyclic loading which in this study was equal to around 40 kPa.Therefore, the CSR is equal to 0.89.The static loading of the pre-consolidated cohesive soil in terms of stress-strain characteristic has ductile characteristic.This type of behaviour is characteristic for soil which was artificially compacted.The pre-cyclic loading of the same type of soil results in change of the behaviour type.Now the brittle type of soil respond is observed.This behaviour type is common for naturally pre-consolidated soils.The behaviour change was caused by increase of the pore water pressure.The maximal deviator stress q max,f was equal to 60.4 kPa, 48.1 kPa and 47.9 kPa respectively.The q max,f value during cyclic loading increased.The reason of such behaviour is the excess pore water pressure and soil skeleton rearrangement which lead to increase of soil strength.
Besides the increase of q max,f , the change of soil mechanical parameters is not observed very well.The slope of critical state line is in all four cases more or less the same.The only exception can be seen in case of the test for 0.5 Hz but this behaviour was caused by temporary technical problems.The slope of critical state line M can be defined as ( 6): where, p' max is the mean effective stress in maximal deviator stress conditions.In case of this tests the M parameter is equal to around 1.65.The stress paths in case of the post cyclic loading show another kind of shearing mechanics than in case of standard static triaxial tests.Usually, in undrained conditions, the vertical loading is cause of pore pressure increase.The excess pore water pressure leads to decrease of effective stress which is acting on soil skeleton and leads to failure.In the case of post cyclic shear tests, the loading leads to slight decrease of the pore pressure.This process ends when rapid change of the pore pressure behaviour is observed.This phenomena is connected with soil failure and the deviator stress, as well as the effective stress, are decreasing.When stress path reaches the critical state, further stress path change of direction is observed.The soil behaviour during undrained shear process after the cyclic loading is applied shows that the cyclic loading generates the excess pore water pressure.The increase of pore pressure is the cause of preconsolidated behaviour of clay samples during shearing.

Conclusions
The cyclic loading of cohesive soil is an important topic for future constructions based on this type of subgrade.The soil under such loading conditions will respond differently when different value of frequency will be applied.Tested soil in this article was compacted in optimal moisture content.The results presented in this article show cohesive soil response to cyclic loading for three different frequencies.The test results and its analysis lead to the following conclusions: 1.
The stress-strain observed during cyclic loading triaxial test, indicates the pattern of soil behaviour.Three phases were identified which refers to different rate of plastic strain accumulation.The plastic strain accumulation characteristics was logarithmic, which means that, with increasing number of cycles the accumulated plastic strain will decrease per one cycle.This type of behaviour was reported before by Goldscheider and Gudehus [36].

2.
The degradation index and the resilient degradation index were employed in order to study the plastic strain accumulation and change in stiffness more precisely.

3.
The degradation characteristics in all cases of frequency loading were similar for first 120 cycles.The long term cyclic loading caused grater degradation in case of sample tested in f equal to 1 Hz and 2 Hz. 4.
The stiffness change analysis was conducted by application of resilient modulus M r which for first cycles of loading was between 70.7 MPa and 79.0 MPa.The resilient degradation index proposed in this article refers to quotient of resilient modulus in N-th cycle to M r in first cycle.5.
The resilient degradation analysis have shown that, despite the greater amount of accumulated plastic strains the samples tested in 1 Hz and 2 Hz will recover its stiffness in long term cyclic loading process.The recovery process was also observed in case of 0.5 Hz frequency loading sample but the process started later and the sample did not recovered its initial stiffness.6.
The excess pore water pressure analysis have shown that the degradation index is connected to the pore pressure generation.The steady value of excess pore which occurs in long term cyclic loading indicates reverse of the soil respond to cyclic loading and after the phase of softening, the hardening occurs.Pore pressure generation in long term cyclic loading of soft soil was studied by Tang et al. [37], the test results have shown, that the excess pore water pressure generation is different than for compacted cohesive soils.The cyclic loading with q a equal to 10 kPa and with σ' 3 equal to 95 kPa lead to generation of ∆u equal to 28 kPa after 2000 repetitions.The frequency in this test was equal to 0.5 Hz.The tests in the same stress conditions and with frequency equal to 2.5 Hz has led to generation of ∆u equal to 19 kPa. 7.
The characteristics of accumulated plastic strain versus number of cycles and the characteristics of excess pore water pressure development were found to be similar.The conclusion can be drown, the excess pore water pressure characteristics impact on the accumulation of the plastic strains.Such phenomena can also explain why in lower frequencies lower accumulated vertical plastic strains are observed.The results of Tang et al. [37] studies for cohesive soft soil shows different picture.The lower frequency causes higher accumulated plastic strains (0.7% and 0.96% respectively for 0.5 Hz and 2.5 Hz).This indicates that there exists the difference between compacted and soft soil response to cyclic loading of cohesive soils.8.
The soil behaviour during undrained shear process after the cyclic loading shows the increase of pore pressure is the cause of preconsolidated behaviour of clay samples during shearing.This type of behaviour differs from this observed during static triaxial test where normally consolidated soil respond was observed.

Figure 1 .
Figure 1.Soil classification (a); Grain size distribution curve; (b) Application of Feret's triangle to determine soil classification

Figure 2 .
Figure 2. Determination of liquid limit of analysed soil.

Figure 1 .
Figure 1.Soil classification (a); Grain size distribution curve; (b) Application of Feret's triangle to determine soil classification

Figure 1 .
Figure 1.Soil classification (a); Grain size distribution curve; (b) Application of Feret's triangle to determine soil classification

Figure 2 .
Figure 2. Determination of liquid limit of analysed soil.

Figure 2 .
Figure 2. Determination of liquid limit of analysed soil.

Figure 3 .
Figure 3. Soil compaction curve based on the Proctor test.

Figure 4 .
Figure 4. Results of oedometric tests on compacted in optimal moisture content clay sample.

Figure 3 .
Figure 3. Soil compaction curve based on the Proctor test.

Figure 3 .
Figure 3. Soil compaction curve based on the Proctor test.

Figure 4 .
Figure 4. Results of oedometric tests on compacted in optimal moisture content clay sample.

Figure 4 .
Figure 4. Results of oedometric tests on compacted in optimal moisture content clay sample.

Figure 5 .
Figure 5. Results of triaxial tests on cohesive soil.The stress-strain characteristic for different frequency of loading.

Figure 6 .
Figure 6.The accumulation of plastic strain characteristics in number of cycles.

Figure 5 .
Figure 5. Results of triaxial tests on cohesive soil.The stress-strain characteristic for different frequency of loading.

Figure 5 .
Figure 5. Results of triaxial tests on cohesive soil.The stress-strain characteristic for different frequency of loading.

Figure 6 .
Figure 6.The accumulation of plastic strain characteristics in number of cycles.

Figure 6 .
Figure 6.The accumulation of plastic strain characteristics in number of cycles.

Figure 7 .
Figure 7. Modulus change characteristics as the degradation index versus number of cycles: (a) degradation index; (b) resilient degradation index.

Figure 7 .
Figure 7. Modulus change characteristics as the degradation index versus number of cycles: (a) degradation index; (b) resilient degradation index.
Figure 8 presents excess pore water pressure development during the cyclic triaxial tests.

Figure 8 .
Figure 8. Excess pore water pressure during cyclic triaxial tests: (a) ∆u value versus vertical strain ε; (b) ∆u value versus number of cycles N.

Table 1 .
Physical properties of tested soil in cyclic triaxial apparatus.