Next Article in Journal
Stress–Strain–Time Description and Analysis of Frozen–Thawed Silty Clay under Low Stress Level
Previous Article in Journal
Rock Mechanical Laboratory Testing of Thebes Limestone Formation (Member I), Valley of the Kings, Luxor, Egypt
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Conversion of Triaxial Compression Strain–Time Curves from Stepwise Loading to Respective Loading

1
School of Civil Engineering, Harbin Institute of Technology, Harbin 150090, China
2
Engineering Research Institute of Appraisal and Strengthening of Shandong Jianzhu University Co., Ltd., Jinan 250014, China
*
Author to whom correspondence should be addressed.
Geotechnics 2022, 2(4), 855-870; https://doi.org/10.3390/geotechnics2040041
Submission received: 10 August 2022 / Revised: 22 September 2022 / Accepted: 27 September 2022 / Published: 30 September 2022
(This article belongs to the Topic Support Theory and Technology of Geotechnical Engineering)

Abstract

:
Numerous researchers of soil creep behavior adopt stepwise loading (SL) rather than respective loading (RL) to perform the triaxial creep tests. However, a complete continuous strain–time curve of SL needs to be converted into assumed curve clusters supposing obtained under RL before the deformation data are used to develop creep constitutive models. Classical methods realize the conversion mainly by focusing on the creep deformation parts and classifying them into linear and nonlinear compositions. Mostly, the linear parts are simply superposed while the nonlinear parts are complex to consider and so are neglected. Moreover, classical methods are not sufficiently valid to eliminate the stress history effect on the conversion. Here, a new method is proposed to achieve the conversion without neglecting the stress history effect. The method rebuilds the triaxial creep test mathematically and physically, adhering to the revising of energy. The method treats the tested deformation in its entirety, instead of distinguishing it into elastic, visco-elastic, plastic and creep (linear and nonlinear) deformation to convert respectively. The comparison among actual measured SL and RL strain–time curves and the curves converted by the new method proves the stress history effect should not be neglected. The higher the vertical load level, the larger the discrepancy between the RL and SL strain–time curve, and the disparity becomes larger with time. The new method highlights the necessity of considering the stress history effect in analysis and design for higher accuracy. The comparisons illustrate the conversion method at least produces more satisfactory results for clayey soil. Primarily examined, at the later stages of loading, the disparity in strain between the converted RL and measured RL decreases by 52.5%~53.5% compared with strain between the measured SL and measured RL.

1. Introduction

Time-dependent stress–strain behavior is a basic research theme of geotechnical engineering. To obtain realistic solutions for time-dependent engineering problems, it is essential to use a constitutive model that accounts for the time dependency of the stress–strain–strength properties of soils. The research on soil creep behavior and properties is a significant branch of soil rheology research. So far, considerable creep constitutive models have been developed from laboratory tests and then adopted in numerical simulation and practical engineering [1]. Therefore, laboratory tests potentially affect the constitutive models.

1.1. Loading Mode

In laboratory triaxial rheological tests, there are two modes to exert the vertical axial stress loads ( σ 1 σ 3 ) : respective loading (RL) and stepwise loading (SL).
SL is exerting a series of loads in turn to only one soil specimen, with each load level sustained until the specimen is deformed extremely slowly towards a relatively steady state. SL produces the stress–strain–time data from a single specimen that experienced a series of loads.
RL is exerting certain constant axial loads to soil specimens one to one. Tests in this loading mode measure and record the strain–time behavior of soil specimens accurately from an initial unpressed state to its failure or unacceptable deformation state.
Usually, RL is considered closer to the concept of research on soil mechanical behavior over time with constant stress load, while in building engineering, the fact is that the total load is not exerted on the soil at one time but gradually, which is coordinated with SL. Even if soil creep behavior is considered from the moment that a building is put into service, there still exists a stress history effect because the soil withstanding forces from the upper building has already been pressed to be dense. Furthermore, if there are not enough rheometers to carry out several tests simultaneously, RL will be time-consuming [2]. Analyzing the data obtained from RL is affected by soil specimens being discrete in physical properties as well as loading conditions (environmental vibrations, temperature, etc.) being slightly different. Therefore, SL has been chosen for numerous creep research studies of soil and rock for a long time up to now [3,4,5,6,7,8,9,10,11,12,13,14].

1.2. Strain–Time Curves and Creep Curves

Here, firstly, the ‘strain–time curve’ is distinguished from the ‘creep curve’. The former contains all the total strain–time points recorded from the moment that the stress load reaches the preset value. The ‘strain’ is composed of elastic, visco-elastic, plastic, and creep strain, while only the latter contains the creep strain–time part. Nowadays, whether the effective stress of the soil specimen is stable is used to judge whether the creep deformation has begun. At an early stage, when the effective stress cannot be measured accurately, there was no unified criterion for the judgement of elastoplastic and creep strain, and some scholars adopt one hour from when the load stress became stable as the time demarcation point for the two types of strain [15,16,17]. There also exist some research works which treat ‘strain–time curve’ as ‘creep curve’ [18], however, here, we do not follow this view.
Though the strain–time curves of SL have practical significance, as mentioned in Section 1.1, researchers prefer the strain–time curves of RL to develop a constitutive model, avoiding the soil stress history effect [19]. Figure 1a shows the single stair-shaped creep curves of RL and Figure 1b shows the multi-stair-shaped creep curves of SL. Though the curves of SL contain the stress history effect of former load levels to latter load levels, researchers began to convert the tested creep deformation curves of SL into assumed curves of RL from an early stage [2,3,20].
The problem is as follows: in the gradually uplifted multi-stair-shaped deformation curve of SL, the latter sections are actually affected by the former sections. There exist micro fissures in the specimen caused by former loads when the latter loads are to be exerted and the fissures exacerbate deformation [21]. How to convert the curves to get rid of the stress history effect remains a problem [2].

1.3. Classical Conversion Methods for Creep Curves

1.3.1. The Coordinate Translation Method (CTM)

This method was too simple to be first thought of to perform the conversion: each beginning moment of a stair loading is regarded as the initial moment of the creep curve under the corresponding RL level, and the subsequent time is counted from this moment, that is, the time coordinate is dragged backward. This method is based on the Boltzmann linear superposition method [22] for the linearly viscoelastic body of rock. Here, the soil is assumed to be linearly rheological [23]. The linearly rheological property means that though the constitutive relationship ( σ ε ) varies with time, at the same moment, the ( σ ε ) relationship always remains linear [24].
The CTM holds that the total creep deformation of the specimen is the linear superposition of the creep deformation increments caused by each load level, and the deformation effect of each load level on the specimen is independent. However, in most cases, the soil is not linearly rheological and the creep deformation is actually affected by the former load levels. What is more, the linear superposition requires that the time of exerting each load level must be consistent, and the creep deformation rate must be close to zero so that the next load level can be exerted. In most creep tests, the time needed for the creep deformation to be stable increases with the load level, so it is difficult to determine each loading duration in advance.

1.3.2. Chen’s Method (Tan’s Method)

Tjong-kie Tan and Wen-fa Kang [3,20] came up with the experimental skill and drawing method. They tested the torsional creep behavior of sandstone and the mode of multi-stepwise moment loading was
M = M 0 H ( t ) + M 1 H ( t θ 1 ) + M 2 H ( t θ 1 ) + M 2 H ( t θ n )
which had been proposed by Tjong-kie Tan earlier in 1964. Towards the multi-stair-shaped creep deformation curve, a geometrical method is used, shown in Figure 2, where b1~b7 are the creep deformation effects caused by the moment increment (921.25 − 612.50) KG-CM and c1~c7 are the creep deformation effects caused by the moment increment (1254.90 − 921.50) KG-CM. b1~b7 are considered as the differences between the actual measured creep deformation and the assumed developed creep deformation along a1~a7 by extending the creep curve under 612.50 KG-CM. The same applies to c1~c7 under the moment increment (1254.90 − 921.50) KG-CM. Then, the creep curve of RL under 921.25 KG-CM can be obtained by a1 + b1, a2 + b2, a3 + b3…a7 + b7 and the curve under 1254.90 KG-CM can be drawn by a1 + b1 + c1, a2 + b2 + c2, a3 + b3 + c3…a7 + b7 + c7. Chen’s method requires the loading duration of each load level to be the same. Some scholars once explained the theory of this method, also related to the Boltzmann linear superposition method [18]. It has been proved that Chen’s method is better than the CTM [25]. So far, Chen’s method has been considered reasonable for clay hollow tube torsional creep tests and solid cylindrical sandstone torsion tests [20], while valid examination of Chen’s method on triaxial creep tests itself has not been directly found by the authors.

1.3.3. Sun’s Conversion Method

Jun Sun considers soil to be non-linearly rheological [2]. Suppose the nonlinear creep response is
ε T ( t ) = σ J 1 ( t ) + σ 2 J 2 ( t , t ) + σ 3 J 3 ( t , t , t ) +
where J 1 ( t ) is linear creep compliance of soil, and J 2 ( t , t ) , J 3 ( t , t , t ) … are a series of nonlinear material functions, all of which are determined by the practical properties and states of the geotechnical materials.
Towards the stress increment Δ σ i at time t i , the creep response is
ε i ( t ) = Δ σ i J 1 ( t t i ) + ( Δ σ i ) 2 J 2 ( t t i , t t i ) + ( Δ σ i ) 3 J 3 ( t t i , t t i , t t i ) +
Different from the linearly rheological hypothesis, Equation (3) reflects the effects of higher-order stress terms. The total response ε T ( t ) at any time t covers both linear and nonlinear effects. The nonlinear effects exist in the contribution of the interactions among the higher-order increments of stress and each stress level.
That is
ε i ( t ) = Δ σ 0 J 1 ( t ) + Δ σ 1 J 1 ( t t 1 ) + Δ σ 2 J 1 ( t t 2 ) + · · · + ( Δ σ 1 ) 2 J 2 ( t , t ) + Δ σ 0 J 2 ( t , t t 1 ) + Δ σ 1 Δ σ 0 J 2 ( t t 1 , t ) + ( Δ σ 1 ) 2 J 2 ( t t 1 , t t 1 ) + Δ σ 0 Δ σ 2 J 2 ( t , t t 2 ) + · · · + ( Δ σ 0 ) 3 J 3 ( t , t , t ) + ( Δ σ 0 ) 2 Δ σ 1 J 3 ( t , t , t t 1 ) + + Δ σ 0 Δ σ 1 Δ σ 2 J 3 ( t , t t 1 , t t 2 ) +
Therefore, the total strain under multi-load levels can be expressed as:
ε T ( t ) =   Δ σ i J 1 ( t t i ) +     Δ σ i Δ σ j J 2 ( t t i , t t j ) +       Δ σ i Δ σ j Δ σ k ( t t i , t t j , t t k ) +
If stress changes continuously with time, the sum form of Equation (5) will become the integrated form. When J 2 , J 3 … is omitted, Equation (5) becomes the Boltzmann linear superposition method. Although Equation (5) is presented in an analytic form, it is difficult for practical application, because the nonlinear material functions J 2 ( t , t ) , J 3 ( t , t , t ) … must be determined by additional experiments.
Comparing the CTM, Chen’s method, and Sun’s conversion method, the former two methods cannot effectively consider the non-linearity of strain–time properties and eliminate the stress history effect, while Sun’s method requires extra tests which consume more time and effort. In this manuscript, the loading mode effect on the triaxial compression strain–time curves is discussed and attempts are made to eliminate it, so that the creep constitutive model of higher accuracy can be proposed and calibrated based on the curves.

2. Method—The Revise Energy Method (REM) for the Conversion of the Strain–Time Curves

As introduced in Section 1.3, the classical methods for the conversion mainly focus on creep strain–time. This is because the researchers only needed the creep strain parts at that time, and the conversion mainly starts from the perspective of ‘strain’, not ‘deformation’. These methods are based on the concept of classifying strain as linear and nonlinear.
The REM will perform the conversion on deformation, treating deformation as an entirety, to obtain the RL strain–time curves corresponding to SL.
Take the example of exerting the first two load levels in turn (SL) and exerting equal loads at one time (RL). Figure 3a shows the process of exerting the first stepwise load F 1 and Figure 3b shows that of adding the second SL F 2 on the same specimen A. Figure 3c shows the RL of the same load F 1 on specimen B. Figure 3d shows the RL of the same load F 1 + F 2 on specimen C.
The SL and RL of F 1 will result in the same deformation. The actual deformation conversion should begin from the + F 2 loading on specimen A and F 1 + F 2 loading on specimen C.
Conservation of energy is usually used as a process calculation criterion from the beginning to the end. When the specimen is deformed to the stable state under the 2nd level SL, the external force work W 2 should be expressed as
W 2 = F 1 ( l 1 + l 2 ) + F 2 l 2   + t 0 t 1 0 l 1 ( σ 3 · A t ) d l d t + t 1 t 2 0 l 2 ( σ 3 · A t ) d l d t + t 0 t 1 0 l x 1 ( σ 3 · A s ) d l x d t + t 1 t 2 0 l x 2 ( σ 3 · A s ) d l x d t
To simplify the explanation of so many symbols in the manuscript, here, we extend the constant subscripts to variable subscripts and explain them. The symbols with superscripts represent similar meanings while under RL, if not explained specially.
W n is the total external force work of n vertical load levels under SL.
n is the level of vertical SL at present or to be calculated.
F i is the ith external vertical load. It is the single load under SL, not the accumulated loads, (N). i = 1, …, n.
l j   ( l x j ) is the vertical stable deformation increment (the mean radial swell increment) under the jth SL compared with the (j − 1)th loading, (mm). j = 1, …, n.
σ 3 is the confining pressure under SL, (MPa). ( σ 3 , under RL).
A t   ( A s ) is the top (side) area of the triaxial specimen under SL, (mm2). ( A t , under RL) ( A s , under RL).
t i 1 ,   t i are, respectively, the beginning time and the stable time of the ith SL, (min). t 0 = 0 .
A s ,   A t ,   σ 3 ( A s ,   A t ,   σ 3 ) change with time.
By setting the value adoption mode of the physical terms in the integral Equation (6), Equation (6) can be simplified as
W 2 = F 1 ( l 1 + l 2 ) + F 2 l 2 + σ 3 · A t 1 · l 1 + σ 3 · A t 2 · l 2 + σ 3 · A s 1 · l x 1 + σ 3 · A s 2 · l x 2
A t i ( A s i ) is the stable top (side) area of the triaxial specimen under the ith level SL, (mm2). ( A t i , under RL F i ) ( A s i , under RL F i ) (Note: the load at this moment is accumulated loads of all the former load levels).
(Note: σ 3 is actually changeable during the test process, but few experiments performed earlier chose to change it because a complex stress history will make the analysis more difficult.)
When the equal level RL test develops to the steady state, the external force work W 2 should be expressed as
W 2 = F 2 l 2   + t 0 t 2 0 l 2 ( σ 3 · A t ) d l d t + t 0 t 2 0 l x 2 ( σ 3 · A s ) d l x d t
where W n is the total external force work of n vertical load levels under RL.
n is the level of vertical RL at present or needs to be calculated. For RL, i = n.
F i is the ith external vertical load. It is the single load under RL, equaling the total quantity from 1st to ith load level under SL, (N). The same, i = n for RL.
l i   is the vertical stable deformation under the F i RL, (mm).
l x i is the mean radial swell increment under the F i RL, (mm).
t i is the time needed until the stable deformation under the F i RL, (mm). t 0 = 0 .
Usually, the conversion of strain–time (deformation–time) curves from SL to RL keeps the confining pressure the same: σ 3 = σ 3 .
Similarly to Equation (7), Equation (8) can be simplified as below
W 2 = F 2 l 2 + σ 3 · A t 2 · l 2 + σ 3 · A s 2 · l x 2
The conversion process of strain–time (deformation–time) curves is to convert the known practical deformation under SL to the unknown deformation under RL. The connection, that is, the consistency equation, could be keeping the total external force work of the ith load level the same between two loading modes:
W 2 = W 2
However, by comparing Equations (7) and (9), there apparently exists a disparity in external force work. With further research, the revised energy term Δ W 2 should be added to W 2 :
Δ W 2 = 2 F 2 l 1 + 2 [ ( F 1 + F 2 ) | θ | + σ 3 · A t 2 · | θ | ]
where θ is a disparity deformation to revise the stress history effect on soil.
θ = l 1 F 1 ( F 1 + F 2 ) ( l 1 + l 2 )
Then Equation (10) is changed into
W 2 + Δ W 2 = W 2
To obtain the l i , Equation (13) is extended by substituting Equations (7), (9) and (11) in, and l i should be the only unknown variable. The strain–time (deformation–time) under SL is already known. Clearly, all the other uncertain geometrical variables must be transformed into expressions of the known variables such as l or the unknown variable l i , so that Equation (13) can be used to derive the l i . The l x j and other geometrical variables may be measured directly if the rheometer is advanced enough [26], which will be considered to be known, making Equation (13) easier.
If the geometrical variables cannot be obtained directly, a hypothesis regarding the geometrical change of the triaxial specimen will push the calculation forward: the triaxial specimen stays in a cylinder shape when it is deformed. This hypothesis refers to the true shape change of specimens in the creep tests, as Figure 4 shows. If there is no friction between the filter paper and the specimen, and also no friction between the filter paper and the porous stone, the cylinder specimen will keep its cylinder shape [27,28]. Of course, the friction is impossible to avoid, but the l x j ,   l x i are both mean values of the radial swell increments, which makes the cylinder hypothesis as approximate as the practical drum-like shape specimen. For UU creep tests, the volume of the cylinder specimen can be assumed not to change, while for the CU and CD, the shrinkage volume of the specimen must be known for higher accuracy.
Following the geometrical hypothesis, rebuild the geometrical change of specimens. Take CD creep tests as an example.
For the 1st level SL, the geometrical Equations (14)–(17) are established
( 1 ε v 1 ) π 4 d 0 2 h 0 = π 4 d 1 2 h 1 = π 4 d 1 2 ( h 0 l 1 )
l x 1 = d 1 d 0
A s 1 = π d 1 h 1 = π d 1 ( h 0 l 1 )
A t 1 = π 4 d 1 2
While for the same quantity 1st level RL,
( 1 ε v 1 ) π 4 d 0 2 h 0 = π 4 d 1 2 h 1 = π 4 d 1 2 ( h 0 l 1 )
l x 1 = d 1 d 0
A s 1 = π d 1 h 1 = π d 1 ( h 0 l 1 )
A s 1 = π 4 d 1 2
where d 0 ( h 0 ) is the initial diameter (height) of the specimen, d i ( h i ) is the stable diameter (height) of the specimen under the ith SL. ε v i is the stable bulk strain of the ith SL. ε v i is the stable bulk strain of the ith RL.
Similarly, the geometrical equations of the 2nd level SL and the same quantity of the 2nd level RL can be obtained.
Obviously, l 1 = l 1 , l x 1 = l x 1 = d 0 ( ( 1 ε v 1 ) h 0 h 0 l 1 1 ) = d 0 ( ( 1 ε v 1 ) h 0 h 0 l 1 1 ) .
To obtain l 2 , a further calculation is needed.
Combine Equations (7), (9), (11)–(13), and the geometrical equations of 2nd level SL and RL. Let F 2 = ( F 1 + F 2 ) , and Equation (22) will be established
W 2 + Δ W 2 = ( F 1 + F 1 ) l 2 + σ 3 A t 2 l 2 + σ 3 A s 2 l x 2 = ( F 1 + F 2 ) l 2 + σ 3 · π 4 d 2 2 l 2 + σ 3 · π d 2 h 2 l x 2 = ( F 1 + F 2 ) l 2 + σ 3 · π 4 d 2 2 l 2 + σ 3 · π d 2 ( h 0 l 2 ) ( d 2 d 0 ) = ( F 1 + F 2 ) l 2 + σ 3 · π 4 d 0 2 ( ( 1 ε v 2 ) h 0 h 0 l 2 ) l 2 + σ 3 · π d 0 ( 1 ε v 2 ) h 0 h 0 l 2 ( h 0 l 2 ) ( d 0 ( 1 ε v 2 ) h 0 h 0 l 2 d 0 )
To simplify Equation (22), let ( 1 ε v 2 ) h 0 h 0 l 2 = ξ , then l 2 = ( ξ 2 1 + ε v 2 ) ξ 2 h 0 , and substitute them back into Equation (22), finally, Equation (23) is
1 4 π σ 3 d 0 2 h 0 ξ 4 + [ ( F 1 + F 2 ) h 0 + 3 4 π σ 3 d 0 2 h 0 3 4 π ε v 2 σ 3 d 0 2 h 0 ( W 2 + Δ W 2 ) ] ξ 2 ( σ 3 π d 0 2 h 0 ε v 2 π σ 3 d 0 2 h 0 ) ξ + ( F 1 + F 2 ) ( ε v 2 1 ) h 0 = 0
To solve l 2 , ξ should be solved first. Fortunately, Equation (23) is a quartic equation of one variable, which is solvable analytically. However, unless the RL is practically performed, ε v 2 will not be determined, so ε v 2 which can be measured in SL is used instead. It should also be noted,
( 1 ε v i ) π 4 d 0 2 h 0 = π 4 d i 2 h i = π 4 d i 2 ( h 0 l i ) d i = d 0 ( 1 ε v i ) h 0 h 0 l i ν = ( d 0 d i d 0 ) / ( l i h 0 ) = ( d 0 d 0 ( 1 ε v i ) h 0 h 0 l i d 0 ) / ( l i h 0 ) = h 0 l i ( 1 ( 1 ε v i ) h 0 h 0 l i )
Here, ν is the Poisson ratio of the triaxial cylinder specimen at the stable state. With the proper time interval selected, the real-time Poisson ratio will also be obtained in the same way.
On working out the ξ , l 2 is also obtained.
Repeat the calculation idea and steps until the nth loading, and Equation (23) will be extended as
1 4 π σ 3 d 0 2 h 0 ξ 4 + [ ( i = 1 n F i ) h 0 + 3 4 π σ 3 d 0 2 h 0 3 4 π ε v n σ 3 d 0 2 h 0 ( W n + Δ W n ) ] ξ 2 ( σ 3 π d 0 2 h 0 ε v n π σ 3 d 0 2 h 0 ) ξ + ( i = 1 n F i ) h 0 ( ε v n 1 ) = 0
ξ = ( 1 ε v i ) h 0 h 0 l i  
l i = ( ξ 2 1 + ε v i ) ξ 2 h 0
Δ W n = 2 i = 2 n ( F i l i 1 ) + 2 [ ( F 1 + i = 2 n F i ) | θ | + σ 3 · A t n · | θ | ]
θ = l 1 F 1 ( F 1 + i = 2 n F i ) ( l 1 + i = 2 n l i )
In solving Equation (23) or (25), there will be four solutions. According to Equation (27), the useful solution ξ is the positive real number which is slightly more than 1.0.
By now, all the calculation steps in the REM have only finished the conversion of stable deformation points of the deformation–time curves while other points at any time between loading moment and a stable state were not converted. To obtain the complete converted strain–time (deformation–time) curves, the method should be used for the whole process of SL on a specimen.
To reduce the considerable calculation amount, for each stair of a deformation–time curve, select a series of calculation points in a kind of time mode, for example, t k = 100 k , k = 0 , 1 , 2 , 3 … (min) as the temporary stable time points until the real end moment of a stair. Carry out the method again by replacing the integral end time point with t k .
Then there will be a series of quartic equations of ξ
1 4 π σ 3 d 0 2 h 0 ξ t k 4 + [ ( i = 1 n F i ) h 0 + 3 4 π σ 3 d 0 2 h 0 3 4 π ε v n , t k σ 3 d 0 2 h 0 ( W n , t k + Δ W n ) ] ξ t k 2 ( σ 3 π d 0 2 h 0 ε v n , t k π σ 3 d 0 2 h 0 ) ξ t k + ( i = 1 n F i ) ( ε v n , t k 1 ) h 0 = 0  
ξ t k = h 0 ( 1 ε v i , t k ) h 0 l t k   ,   l t k = ( ξ t k 2 1 + ε v i , t k ) ξ t k 2 h 0
where W n , t k   is the total external force work up to the moment of t k at nth stair of a deformation–time curve, (n = 1, 2, 3, … t k = 100 k , k = 1 , 2 , 3 ). The simplified W n , t k   is
W n , t k = i = 1 n = 1 ( F i · j = i n 1 l j ) + F n · l t k   + i = j = 1 n 1 ( σ 3 · A t i · l j ) + σ 3 · A t t k · l t k + i = j = 1 n 1 ( σ 3 · A s i · l x j ) + σ 3 · A s t k · l x t k
where l t k ( l x t k ) is up to the t k moment, the vertical deformation increment (the mean radial swell increment) under the SL of i = 1 n F i level, (mm).
A t t k ( A s t k ) is up to the t k moment, the top (side) area of the triaxial specimen under SL of i = 1 n F i level, (mm2).
l t k is the vertical deformation up to the t k moment under the F i RL converted from SL, (mm). l t k is also the variable that the REM concerns.
ε v n , t k ( ε v i , t k ) is up to the t k moment, the bulk strain of the specimen under nth (ith) SL.
Perform the conversion at t k moments, respectively. With the points converted, a complete deformation–time curve of RL can be drawn.

3. Experiment and Result—An Instance for the Revise Energy Method

3.1. Experiment

Here, to examine the validity of the revise energy method, Zhou’s [29] triaxial creep experimental results of the remolded saturated clayey soil are adopted. The basic physical properties are listed in Table 1.
Through triaxial CU tests, the cohesion of the remolded clayey soil is 2.37 kPa, and the internal friction angle is 16.86°. The triaxial creep test plan was designed as Table 2.
Apparently, Test 4 and 5 are in a stepwise loading path while Test 1, 2, and 3 are in their respective loading path. The specimens tested are all 61.8 mm ( d 0 ) × 125 mm ( h 0 ). The stable criteria for deformation is < 0.01   mm / d . Axial strain–time and bulk strain–time data are both extracted for conversion.

3.2. Result

The conversion axial strain–time curves of ( σ 1 σ 3 ) = 120   kPa   &   240   kPa in Test 4 are illustrated in Figure 5 and Figure 6, respectively, and they are also compared with Test 2 and Test 3. While Figure 7 compares the conversion of ( σ 1 σ 3 ) = 240   kPa in Test 5 and the respective loading result of Test 3. All the conversions are performed with bulk strain tested in stepwise loading because the conversion itself is based on not knowing the bulk strain of the respective loading.
In Figure 5, soil specimens exerted the same external stress load 120 kPa, the SL strain–time curve firstly exceeds the RL curve, while, along with time, the RL curve exceeds the SL curve. The converted strain ε 2 –time (corresponding to l 2 –time) curve exhibits a higher level than both strain–time of SL and RL. Focusing on the trend that curves develop, in general, the later stage disparity between the converted curve and RL curve becomes smaller, while the disparity between the RL and SL curves becomes larger.
In Figure 6, soil specimens exerted the same external stress load 240 kPa, the RL strain–time curve quickly exceeds the SL curve and later exceeds the converted curve. The converted strain ε 4 –time (corresponding to l 4 –time) curve exhibits a middle level between the strain–time curves of SL and RL. Focusing on the trend that curves develop, the later stage disparity among the three curves shows an inconspicuous change.
In Figure 7, under the same external stress load of 240 kPa, the RL strain–time curve quickly exceeds the SL curve and the converted curve at an early moment. The converted strain ε 2 –time (corresponding to l 2 –time) curve exhibits a middle level between the strain–time curves of SL and RL. The later stage disparity among the three curves also shows an inconspicuous change.

4. Discussions—Analysis of the Revise Energy Method

4.1. The Trend Analysis

Usually, within the comparisons of the curves, the RL curve shows a dramatic increase, this is because the RL curve begins from 0 strain at 0 min, while the SL curve, unless it is the first stair, will always begin from a non-zero value in the actual tests and so too the converted curve. In some ways, the differences in beginning points among the curves are the embodiment of the stress history effect. The slope of the RL curve at the beginning stage is greater than the RL and converted curves, corresponding to the elastic modulus variation of the soil specimens. The RL curve is obtained from a specimen in an originally non-loaded state, while the SL curve is obtained from a specimen in an already loaded state, which signifies the latter specimen possesses greater density and higher rigidity (the soil tested is strain-hardening type), and this is also the explanation of stress history effect on the material. It firmly proves again that if the SL curve is adopted for analysis or utilized without conversion, an unsatisfying discrepancy will occur.
The non-zero beginning point cannot be avoided in the conversion because the conversion should be based on already obtained results. To change this, it could be achieved by manually revising the first point of the curve to be (0,0).
At the later stage, with enough loading duration, the specimens are deformed towards a stable state, so the developing tendency of the curves is roughly similar. The tendency of the curves in Figure 5 illustrates the specimens may not have crept to become sufficiently stable.

4.2. The Disparity Analysis

Select the endpoints of all the curves to compare, the strain values are listed in Table 3.
From the comparison in Table 3, the end strain of the SL curve is generally smaller than it of the RL curve when the loading condition is the same. Comparing the different loading conditions, the larger the axial stress load is, the larger the discrepancy between the RL and SL strain–time curve is, and the disparity becomes larger with time.
As shown in Figure 8, the strain corresponding to stress ( σ 1 σ 3 ) = 240 kPa but loaded in different modes also varies. From these tests, in general, the more steps a specimen experiences, the lower the strain level is. These results can both be explained by the stress history effect and the strain hardening property of clayey soil.
As listed in Table 3, the disparity between the actual RL strain and converted strain reaches 0.38% strain. While the disparity between the actual RL and SL reaches 0.8% strain. Therefore, the revise energy method at least improves the accuracy of conversion results at the later stage of deformation with time. Compared with the original SL curves, the converted curves validly approaches the actual RL curves when the stress load becomes larger, while the remaining disparity that the REM is not able to improve may be explained by the fact that:
1
The cylinder geometrical hypothesis is not totally applicable to the test results of specimens in shape, as in Figure 4.
2
The bulk strain data used in the conversion is tested in SL, not the bulk strain tested in RL, the disparity caused by the discrepancy of bulk strain may be accumulated with the stress load increasing.
3
The initial size of the specimen may not accurately be d 0 61.8 mm × h 0 125 mm, which will not influence the RL strain–time curve decisively but will influence the conversion of REM because d0 and h0 are vital parameters for the REM.
The gravity of a triaxial specimen can be neglected in the calculation, however, the gravity of soil cannot be neglected when such a vast amount of soil is involved in the calculation, and the friction must not be neglected either. The practical effect of gravity and friction on the conversion for a large soil test remains for further study.

5. Conclusions

The comparison between actual measured SL and RL strain–time curves proves the stress history effect should not be neglected and highlights the necessity of considering the stress history effect in analysis and design for higher accuracy. For stress hardening soil, the higher the vertical load level is, the larger the discrepancy between the RL and SL strain–time curve is, and the disparity becomes larger with time. In general, to reach a certain stress level, the more steps a specimen experiences, the lower the strain level is.
Through revising the energy relations, the revise energy method rebuilds the triaxial creep test mathematically and physically, eventually achieving the conversion of strain–time curves without neglecting the stress history effect. The method treats the tested deformation as an entirety, instead of distinguishing it into elastic, visco-elastic, plastic and creep (linear and nonlinear) deformation to convert respectively.
The revise energy method does not require the SL strain–time curve to be of equal duration for each stress load level. With examination, the revise energy method is at least feasible to convert the later stage strain–time curves of triaxial creep tests of clayey soil. Through primary research, at the later stages of loading, compared with the disparity of strain between the measured SL and measured RL, the disparity of strain between the converted RL and measured RL decreases by 52.5%~53.5%.
With the revise energy method, much closer results to the practical results of RL creep tests can be obtained by converting the SL results so that the researchers can still choose SL in the soil triaxial creep test with the advantages of SL. Based on the more reliable converted strain–time curves, more accurate creep constitutive models are expected to be developed, and accordingly, engineering activities will benefit from the guidance of analysis on more practical strain–time curves and better constitutive models.

Author Contributions

Conceptualization, H.Q.; methodology, H.Q.; Software, H.Q., D.M.; validation, H.Q.; formal analysis, H.Q.; investigation, Z.N.; resources, A.T.; data curation, H.Q.; writing—original draft preparation, H.Q.; writing—review and editing, H.Q., D.M., Z.N., A.T.; visualization, H.Q.; supervision, A.T.; project administration, A.T., H.Q.; funding acquisition, A.T. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by [National Natural Science Foundation of China] grant number [41672287] and the APC was funded by [Harbin Institute of Technology]. The author Qu, H. was also funded by [China Scholarship Council] grant number [202106120229] during the preparation, submission and publication period of this research.

Data Availability Statement

The data presented in this study are available on request from the corresponding author. The data are not publicly available due to the original authors of the data didn’t present the data directly in number but in illustration in the cited literature.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Liingaard, M.; Augustesen, A.; Lade, P.V. Characterization of Models for Time-Dependent Behavior of Soils. Int. J. Géoméch. 2004, 4, 157–177. [Google Scholar] [CrossRef]
  2. Sun, J. Rheology of Geotechnical Materials and Its Engineering Application; China Construction Industry Press: Beijing, China, 1999. [Google Scholar]
  3. Zhai, Y.; Wang, Y.; Dong, Y. Modified Mesri Creep Modelling of Soft Clays in the Coastal Area of Tianjin (China). Teh. Vjesn.-Tech. Gaz. 2017, 24, 1113–1121. [Google Scholar] [CrossRef]
  4. Jiang, Y.Z.; Wang, R.H.; Zhu, J.B. Rheological mechanical properties of altered-rock. Mater. Res. Innov. 2015, 19, S5-349–S5-355. [Google Scholar] [CrossRef]
  5. Ou, Z.-F.; Fang, Y.-G. The Influence of Organic Matter Content on the Rheological Model Parameters of Soft Clay. Soil Mech. Found. Eng. 2017, 54, 283–288. [Google Scholar] [CrossRef]
  6. Han, B.; Fu, Q. Study on the Estimation of Rock Rheological Parameters under Multi-level Loading and Unloading Conditions. MATEC Web Conf. 2018, 213, 02003. [Google Scholar] [CrossRef]
  7. Li, K.-F.; Liu, R.; Qiu, C.-L.; Tan, R.-J. Consolidated Drained Creep Model of Soft Clay in Tianjin Coastal Areas; Springer: Singapore, 2018; pp. 157–165. [Google Scholar] [CrossRef]
  8. Wang, Z.; Wong, R.C.K. Strain-Dependent and Stress-Dependent Creep Model for a Till Subject to Triaxial Compression. Int. J. Géoméch. 2016, 16, 04015084. [Google Scholar] [CrossRef]
  9. Li, J.; Tang, Y.; Feng, W. Creep behavior of soft clay subjected to artificial freeze–thaw from multiple-scale perspectives. Acta Geotech. 2020, 15, 2849–2864. [Google Scholar] [CrossRef]
  10. Hu, K.; Shao, J.-F.; Zhu, Q.-Z.; Zhao, L.-Y.; Wang, W.; Wang, R.-B. A micro-mechanics-based elastoplastic friction-damage model for brittle rocks and its application in deformation analysis of the left bank slope of Jinping I hydropower station. Acta Geotech. 2020, 15, 3443–3460. [Google Scholar] [CrossRef]
  11. Zhang, G.; Chen, C.; Zornberg, J.G.; Morsy, A.M.; Mao, F. Interface creep behavior of grouted anchors in clayey soils: Effect of soil moisture condition. Acta Geotech. 2020, 15, 2159–2177. [Google Scholar] [CrossRef]
  12. Yao, X.; Qi, J.; Zhang, J.; Yu, F. A one-dimensional creep model for frozen soils taking temperature as an independent variable. Soils Found. 2018, 58, 627–640. [Google Scholar] [CrossRef]
  13. Bhat, D.; Kozubal, J.; Tankiewicz, M. Extended Residual-State Creep Test and Its Application for Landslide Stability Assessment. Materials 2021, 14, 1968. [Google Scholar] [CrossRef]
  14. Liu, Z.B.; Xie, S.Y.; Shao, J.F.; Conil, N. Multi-step triaxial compressive creep behaviour and induced gas permeability change of clay-rich rock. Géotechnique 2018, 68, 281–289. [Google Scholar] [CrossRef]
  15. Zhang, Y.Q. Research on Creep Model Test of Rockfill Material. Master’s Thesis, China Academy of Water Resources and Hydropower Sciences, Beijing, China, June 2013. [Google Scholar]
  16. Shen, Z.J.; Zuo, Y.M. Experimental research on rheological properties of rockfill. In Proceedings of the 6th Symposium on Soil Mechanics and Foundation Engineering of the Chinese Society of Civil Engineering, Shanghai, China, 18–22 June 1991; pp. 443–446. [Google Scholar]
  17. Cheng, Z.L.; Ding, H.S. Creep test for rockfill. Chin. J. Geotech. Eng. 2004, 26, 473–476. [Google Scholar]
  18. Zhao, L.; Gan, H.X.; Tang, W.; Tang, M.F.; Zhou, H.P. Applicability Analysis of Chen’s Method in the Research of TATB-based PBX Creep Behavior. Chin. J. Energetic Mater. 2018, 26, 608–613. [Google Scholar] [CrossRef]
  19. Zhao, X.; Wang, Q.; Gao, Y.; Yang, T.; He, L.; Lu, X. Indoor Rheological Test and Creep Model Analysis of Soft Soil in Qingyi River region in Wuhu Anhui. IOP Conf. Series: Earth Environ. Sci. 2019, 330, 032007. [Google Scholar] [CrossRef]
  20. Tan, T.K.; Kang, W.F. Locked in stresses, creep and dilatancy of rocks, and constitutive equations. Rock Mech. 1980, 13, 5–22. [Google Scholar] [CrossRef]
  21. Zhao, D.; Gao, Q.-F.; Hattab, M.; Hicher, P.-Y.; Yin, Z.-Y. Microstructural evolution of remolded clay related to creep. Transp. Geotech. 2020, 24, 100367. [Google Scholar] [CrossRef]
  22. Nie, Z. Experimental Investigation on the Long-Term Properties of Expansive Soils Subjected to Freeze-Thaw Cycling. Master’s Thesis, Harbin Institute of Technology, Harbin, China, June 2018. [Google Scholar]
  23. Tan, T.K. Basic Nature and Engineering Suggestions of Loess in Northwest China. Chin. J. Geotech. Eng. 1989, 4, 9–24. [Google Scholar]
  24. Sun, J. Rock rheological mechanics and its advance in engineering applications. Chin. J. Rock Mech. Eng. 2007, 26, 1081–1106. [Google Scholar]
  25. Zhang, X.W.; Wang, C.M.; Zhang, S.H. Comparative analysis of soft clay creep data processing method. J. Jilin Univ. (Earth Sci. Ed.) 2010, 40, 1401–1408. [Google Scholar] [CrossRef]
  26. Kishi, M.; Tani, K. Development of measuring method for axial and lateral strain distribution using CCD sensor in triaxial test. In Deformation Characteristics of Geomaterials; IS: Lyon, France, 2003; pp. 31–36. [Google Scholar] [CrossRef]
  27. Chu, J.; Lo, S.R.; Lee, I.K. Instability of Granular Soils under Strain Path Testing. J. Geotech. Eng. 1993, 119, 874–892. [Google Scholar] [CrossRef]
  28. Li, S.; Huang, X.; Wu, S.; Tian, J.; Zhang, W. Bearing capacity of stabilized soil with expansive component confined by polyvinyl chloride pipe. Constr. Build. Mater. 2018, 175, 307–320. [Google Scholar] [CrossRef]
  29. Zhou, P.J. Study on Shear Creep Characteristics and Model of Saturated Cohesive Soil. Master’s Thesis, Zhejiang University of Technology, Hangzhou, China, June 2018. [Google Scholar]
Figure 1. (a): RL mode and corresponding creep strain–time curves. (b): SL mode and corresponding creep strain–time curves [2].
Figure 1. (a): RL mode and corresponding creep strain–time curves. (b): SL mode and corresponding creep strain–time curves [2].
Geotechnics 02 00041 g001
Figure 2. SL creep curve for sandstone and its conversion to RL creep curve [20]. The colors of angle of twist per cm–time lines represent the following: gray line, under 612.50 KG-CM torque; red lines, under 921.25 KG-CM torque; blue lines, under 1254.90 KG-CM torque; yellow line, under other later torque levels.
Figure 2. SL creep curve for sandstone and its conversion to RL creep curve [20]. The colors of angle of twist per cm–time lines represent the following: gray line, under 612.50 KG-CM torque; red lines, under 921.25 KG-CM torque; blue lines, under 1254.90 KG-CM torque; yellow line, under other later torque levels.
Geotechnics 02 00041 g002
Figure 3. Illustration of SL process on specimen A and RL process on specimen B and C. (a) Specimen A is deformed under F1 load. (b) Specimen A is continually deformed under another F2 load after it reached a stable deformation state under F1 load. (c) Specimen B is deformed under F1 load. (d) Specimen C is deformed under (F1 + F2) load.
Figure 3. Illustration of SL process on specimen A and RL process on specimen B and C. (a) Specimen A is deformed under F1 load. (b) Specimen A is continually deformed under another F2 load after it reached a stable deformation state under F1 load. (c) Specimen B is deformed under F1 load. (d) Specimen C is deformed under (F1 + F2) load.
Geotechnics 02 00041 g003
Figure 4. Soil specimen instances after creep tests. Specimen (a), under presupposed confining pressure σ3 = 200 kPa. Specimen (b), under presupposed confining pressure σ3 = 150 kPa. Specimen (c), under presupposed confining pressure σ3 = 100 kPa. ( σ 1 σ 3 ) exerted to the three specimens ranges from 37.98 kPa~360.23 kPa. The original sizes of all the three specimens are d0 = 61.8 mm and h0 = 125 mm.
Figure 4. Soil specimen instances after creep tests. Specimen (a), under presupposed confining pressure σ3 = 200 kPa. Specimen (b), under presupposed confining pressure σ3 = 150 kPa. Specimen (c), under presupposed confining pressure σ3 = 100 kPa. ( σ 1 σ 3 ) exerted to the three specimens ranges from 37.98 kPa~360.23 kPa. The original sizes of all the three specimens are d0 = 61.8 mm and h0 = 125 mm.
Geotechnics 02 00041 g004
Figure 5. Comparison of converted RL axial deformation and tested RL and SL axial deformation when ( σ 1 σ 3 ) = 120   kPa (Test 4, 2).
Figure 5. Comparison of converted RL axial deformation and tested RL and SL axial deformation when ( σ 1 σ 3 ) = 120   kPa (Test 4, 2).
Geotechnics 02 00041 g005
Figure 6. Comparison of converted RL axial deformation and tested RL and SL axial deformation when ( σ 1 σ 3 ) = 240   kPa (Test 4, 3).
Figure 6. Comparison of converted RL axial deformation and tested RL and SL axial deformation when ( σ 1 σ 3 ) = 240   kPa (Test 4, 3).
Geotechnics 02 00041 g006
Figure 7. Comparison of converted RL axial deformation and tested RL and SL axial deformation when ( σ 1 σ 3 ) = 240   kPa (Test 5, 3).
Figure 7. Comparison of converted RL axial deformation and tested RL and SL axial deformation when ( σ 1 σ 3 ) = 240   kPa (Test 5, 3).
Geotechnics 02 00041 g007
Figure 8. Comparison of axial strain tested in different loading modes when ( σ 1 σ 3 ) = 240   kPa (Test 3, 4, 5).
Figure 8. Comparison of axial strain tested in different loading modes when ( σ 1 σ 3 ) = 240   kPa (Test 3, 4, 5).
Geotechnics 02 00041 g008
Table 1. Basic physical properties of clayey soil.
Table 1. Basic physical properties of clayey soil.
Source
of Soil
Sampling
Depth
(m)
Specific Gravity
Gs
Moisture Content
ω (%)
Void Ratio e 0 Water Limit
w L (%)
Plastic Limit w P   (%) Density ρ (g/cm3) Particle Gradation (mm)
>0.0750.05–0.075<0.05
West Lake,
Hangzhou,
China
15–252.74280.8739.2221.8722.693.2194.1
Table 2. Triaxial creep test plan.
Table 2. Triaxial creep test plan.
TestConfining Pressure
σ 3 (kPa)
Shearing Strength
( σ 1 σ 3 ) f (kPa)
Loading Path
( σ 1 σ 3 )
140033060
2400330120
3400330240
4400330 60 120 180 240
5400330 120 240
Table 3. End strain of SL, RL, converted strain–time curve of three loading conditions.
Table 3. End strain of SL, RL, converted strain–time curve of three loading conditions.
Figure
(Loading Condition)
End Strain (%)Disparity of Strain between RL and Converted (%)Disparity of Strain between RL and SL (%)
SL CurveRL CurveConverted Curve
1 (120 kPa of Test 2 and 4)1.061.111.320.21−0.05
2 (240 kPa of Test 3 and 4)2.233.032.65−0.38−0.80
3 (240 kPa of Test 3 and 5)2.593.152.89−0.26−0.56
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Qu, H.; Mu, D.; Nie, Z.; Tang, A. Conversion of Triaxial Compression Strain–Time Curves from Stepwise Loading to Respective Loading. Geotechnics 2022, 2, 855-870. https://doi.org/10.3390/geotechnics2040041

AMA Style

Qu H, Mu D, Nie Z, Tang A. Conversion of Triaxial Compression Strain–Time Curves from Stepwise Loading to Respective Loading. Geotechnics. 2022; 2(4):855-870. https://doi.org/10.3390/geotechnics2040041

Chicago/Turabian Style

Qu, Haigang, Dianrui Mu, Zhong Nie, and Aiping Tang. 2022. "Conversion of Triaxial Compression Strain–Time Curves from Stepwise Loading to Respective Loading" Geotechnics 2, no. 4: 855-870. https://doi.org/10.3390/geotechnics2040041

APA Style

Qu, H., Mu, D., Nie, Z., & Tang, A. (2022). Conversion of Triaxial Compression Strain–Time Curves from Stepwise Loading to Respective Loading. Geotechnics, 2(4), 855-870. https://doi.org/10.3390/geotechnics2040041

Article Metrics

Back to TopTop