# Application of a Viscoelastic Model to Creep Settlement of High-Fill Embankments

Authors: Liang Jia; Guangli Huang
Journal: Advances in Civil Engineering
 (2019)

---

## Abstract

In order to predict the creep settlement of high-fill embankments, the time-dependent viscoelastic model of Poynting–Thomson (the standard linear solid) has been chosen to represent the creep behavior of soils. In the present study, the hereditary integral was applied to calculate the strain while the load increase is varied with time. Calculation expressions of the creep settlement of an embankment during and after construction were obtained under one-dimensional compression conditions. Using this approach, the three parameters of every layer can be determined and adjusted to accommodate in situ monitoring data. The calculated results agreed well with those from the field, which imply that the method proposed in this paper can give a precise prediction of creep settlement of high-fill embankments.


---

## Body

## 1. Introduction

Creep settlement is a main research issue for highway deterioration and other high-fill projects [1–5]. The analysis of embankment creep settlements is the basis for the correct understanding and evaluation of the subgrade stability. In consequence, more research should be carried out into the calculation of creep settlement of embankment. Existing methods related to creep settlement research can be divided into three groups: (1) Empirical models, for example, the Asaoka method [6, 7], the hyperbolic method [8], the parabola method [9], the exponential method, and in situ tests [10, 11], in those methods, the monitored creep settlement of the embankment after construction was adopted to predict the creep settlement in the future. (2) Constitutive models: Bjerrum [12] pointed out that the influence of creep on clays was equivalent to make the clay overconsolidated, which was obtained through a large number of consolidation tests and engineering settlement data analysis, and Yao et al. [13] established the unified hardening (UH) model (a model for over-consolidated clays), which considered the time effects and analyzed the influence of creep on over-consolidated clays. Hashiguchi [14] put forward the concept of a subloading yield surface. Asaoka et al. [15] suggested the concept of superloading and subloading yield surfaces for over-consolidated soils. And Wang et al. [16] proposed an advanced elastic-viscoplastic constitutive model to predict time-dependent behavior of the over-consolidated clay based on the subloading surface model and relative overstress relation. (3) Numerical methods: Loganathan et al. [17] put forward a calculation method named “field deformation analysis” (FDA), which separates and quantifies settlement components. Li et al. [18] used the modified secant modulus method to calculate settlements, an approach that is characterized by the fact that it is not influenced by the initial void ratio. But the empirical models pay less attention on the creep mechanism of clays, and most of the constitutive models are elastoplastic constitutive models, which cannot describe the time-dependent deformation characteristic of over-consolidated clays. Experimental results indicate that the creep deformation of high-filled embankment is affected by both time and stress history. The rheological model can reflect the time-dependent stress-strain behavior of clays, and in order to simulate the layered filling process of the high-filled embankment during construction, the stress applied on the rheological model is according to the loading curve, so the rheological model used in the creep settlement calculation of the embankment can not only reflect the time-dependent stress-strain behavior but also the stress history of the embankment. The application of some relatively simple rheological models has been studied in some cases [19–22]. Morro [23] concluded that a spring-dashpot model (i.e., the Maxwell model in parallel with a dashpot) showed a constant strain rate during the secondary stage. Justo and Durand [24] obtained settlement expressions with the standard linear solid rheological model for one-dimensional settlements of an embankment during and after construction. In this paper, the Poynting–Thomson model (a linear spring in parallel with a Maxwell model) was applied to represent the time-dependent deformation of soil, and the theory of the hereditary integral used to calculate the cumulative strain is explained in detail, which reflects the strain at any time factor T and depends upon the entire stress history [25, 26].As an illustration of the creep settlement performance of a high-filled embankment, the case of the dam-like embankment of the Lanzhou–Yongjing highway was investigated using the proposed new method.


## 2. Theory

For the model of Poynting–Thomson (Figure1),(1)σ=σM+σH,(2)ε=εM=εH,where(3)σH=EHεH=EHε.Figure 1
Viscoelastic model.The constitutive equation of the Maxwell model in Figure1 is(4)ηEMσM˙+σM=ηεM˙,where(5)σM˙=σ˙−σH˙,(6)εM˙=ε˙.Substituting equations (5) and (6) into equation (4) yields(7)ηEMσ˙−σH˙+σM=ηεM˙.The calculation expression ofσM is obtained from equation (7):(8)σM=ηεM˙−ηEMσ˙−σH˙,and substituting the expressions of σM and equation (3) into equation (1) yields(9)σ−EHε=ηε˙−ηEMσ˙−σH˙=ηε˙−ηEMσ˙−EHε˙.The constitutive equation of the Poynting–Thomson model is(10)ηEMσ˙+σ=ηEMEM+EHε˙+EHε.In a creep test, the pressure maintains a constant valueσ0 when t=0; equation (10) becomes(11)σ=ηEMEM+EHε˙+EHε.The Poynting–Thomson model (Figure1) first produces the instantaneous elastic deformation ε0=σ0/EM+EH under the constant pressure σ0. As time increases, the deformation continues to increase such that when t⟶∞, ε⟶σ0/EM, and the solution of equation (11) is(12)εT=σ0EH1−EMEM+EHe−EMEH/EM+EHηt.Which otherwise can be written as(13)εt=σ0EH1−Re−T,where(14)R=EMEM+EH.R is called the creep factor.(15)Tr=EMEHEM+EHη.T
=
t
/
T
r is called the time factor.And equation (13) can also be written as(16)εt=σ0JT,where(17)JT=1EH1−Re−T.The mechanical meaning ofJT is the strain under unit force, which can easily be measured in practical tests. JT is called the creep compliance. After applying the stress σ0 when T=0, the strain of the model is σ0JT. Then, when another stress Δσ′ is applied when T=T′, there is another strain Δσ′JT−T′, which is added to σ0JT when T=T′, which is proportional to Δσ′. This strain is also a function of T, but only adds εt when T=T′ and is still dependent upon the same creep compliance:(18)εT=σ0JT+Δσ′JT−T′.IfΔσ″ is applied when T=T″, then εT should increase Δσ″JT−T″. The effect of each Δσ on εT can be superposed in this way, which is called the Boltzmann superposition principle [27]. If σ changes when T=0, the relationship between applied stress and the time factor T will be replaced by a stepped polyline, as shown in Figure 2. When dT′⟶0, the sum becomes integral:(19)εT=∫0TJT−T′dσT′.Figure 2
Sketch of the loading curve.If the applied load varies with time, according to Figure3, because σ0=0, the strain at time factor T will be given by the hereditary integral:(20)εT=∫0T1EH1−Re−T−T′dσdT′dT′.Figure 3
Variation of loading with time.The solution of the hereditary integral is(21)εT=1EHσ−Re−TI1,where(22)I1σ,T=∫0TeT′dσdT′dT′.The strain fromT to Tc (Figure 3) can be described as(23)ΔεT,Tc=εTc−εTT≤Tc.Δ
ε
T
,
T
c is obtained from equation (21):(24)ΔεT,Tc=1EHσc−σ−Re−TcI1σ,Tc−e−TI1σ,T.And whenT≥Tc, the strain from Tc to T is (Figure 3)(25)ΔεTc,T=εTc−εT.From equation (21), it is supposed that the stress maintains a constant value when T≥Tc (Figure 3):(26)ΔεTc,T=εTc+DT−εTc=REHe−TcI1σ,Tc1−e−DT.


## 3. Application to a High-Filled Embankment

If it is assumed that the height of an embankment varies with time according to Figure4, the stress at any point depends on the height of fill above it.Figure 4
Variation of fill height with time.


### 3.1. Settlement of the Embankment during Construction

#### 3.1.1. Degree of Settlement of the Embankment at HeightzI during Construction

The degree of settlement of the embankment at heightzI during construction can be described as follows:(27)SzI=∫0zIΔεTI−T,Tc−Tdz.Δ
ε can be obtained from equation (21):(28)ΔεEHγ=ZI−Re−TcI1z,Tc−I1z,TI−e−TII1z,TI−I1z,T,where(29)I1z,TI=∫0TIeT′dzdT′dT′,(30)I1z,Tc=∫0TceT′dzdT′dT′=∫0TIeT′dzdT′dT′+∫TITceT′dzdT′dT′.The integrals listed above must be calculated by numerical methods. In order to do this, the loading curve is divided intoN intervals; from T=0 up to T=TI, the curve is divided into I intervals, and from T=TI up to T=Tc, the curve is divided into N−I intervals.Substituting equation (28) into equation (27), the settlement of a point at height zI during construction is obtained from(31)SzIEHγ=ZIzI−Ra−e−TII2z,TI,where(32)a=e−TczII1TI,Tc+I2z,TI,(33)I1TI,Tc=∫TITceTdzdTdT=∑i=INeTiΔzi,(34)zi=zI+1−zI−12,zI−1=zI,zN+1=zN,(35)I2z,TI=∫0TIzeTdzdTdT=∑i=1IzieTiΔzi,zI+1=zN.If it is assumed that the height of the embankment increases linearly with time, and in order to simplify the deduction process, we assumed that the slope of the loading curve is 1, and then equations (33) and (35) will become(36)I1TI,Tc=∫TITceTdzdTdT=∫TITceTdT=eTc−eTI,(37)I2z,TI=∫0TIzeTdzdTdT=∫0TITeTdT=TIeTI−eTI+1.Substituting equations (36) and (37) into equation (31) gives(38)SzIEHγ=ZIzI−R1+e−H−e−ZI−e−zI.


#### 3.1.2. Settlement of a Finite Layer (zI,zII) during Construction

The strain from time factorTII up to Tc is(39)Δε=ΔεTII−T,Tc−T.SubstitutingTI by TII and equation (39) into equation (21) gives(40)ΔεEHγ=H−zII−Re−TcI1z,Tc−I1z,T−e−TIII1z,TII−I1z,T.The compression of a finite layer is(41)SzI,zII=∫zIzIIΔεdz.Substituting equation (40) into equation (41) gives(42)SzI,zIIEHγ=H−zIIΔz−Re−TcΔzI1TII,Tc+bΔze−Tc−e−TII,where(43)Δz=zII−zI,(44)b=I2TI,TII−zII1TI,TIIΔz,(45)I2TI,TII=∫TITIIzeTdzdTdT=∑i=Ii=IIzieTiΔzi,(46)I1TI,TII=∫TITIIeTdzdTdT=∑i=Ii=IIeTiΔzi,(47)Δzi=zi+1−zi−12,zII+1=zII,zII+1=zII.


### 3.2. Settlement of the Embankment after Construction

#### 3.2.1. Settlement of an Embankment at HeightzI after Construction

The settlement of an embankment at heightzI after construction can be described as follows:(48)ΔSzI=∫0zIΔεTc,Tc+DTdz.Substituting equation (26) into equation (48) gives(49)ΔSzIEHγR=a1−e−DT,where the value of a is equal to equation (32).At the top of the embankment,(50)I1Tc,Tc=0,(51)ΔSHEHγR=e−Tc1−e−DTI2z,Tc.If the height of the embankment varies linearly with time, with a slope of 1, then equations (49) and (51) will become(52)ΔSzIEHγR=zI−e−zI+e−H1−e−DT,(53)ΔSHEHγR=H+e−H−11−e−DT.


#### 3.2.2. Settlement of a Finite LayerzI,zII after Construction

The strain of a finite layerzI,zII after construction is(54)Δε=ΔSzI,zIIΔz=ΔSzII−ΔSzIΔz,The settlement of a finite layerzI,zII after construction is(55)ΔSzI,zIIEHγR=e−Tc1−e−DTΔzI1TII,Tc+b,where(56)I1TII,Tc=∫TIITceTdzdTdT=∑i=IIi=NeTiΔzi,where the value of b is equal to equation (44).


### 3.3. Parameter Determination

The three parameters of the model are the elastic modulusEH, the creep ratio R, and the time factor Tr.If settlement is measured from different heights of the embankment, from equation (55), the time factor Tr is obtained from(57)ΔSTc,T2ΔSTc,T1=1−e−DT21−e−DT1=1−e−Δt2/Tr1−e−Δt1/Tr,where Δt is the time increment after construction.The solution of equation (57) with different time increments will permit an adjustment for Tr.Substituting equation (55) into equation (42), the values of R from different layers are obtained:(58)R=ZIIe−TcI1TII,Tc+bSzI,zII/ΔSzI,zII1−e−DT+1−be−TII.Substituting the value ofTr and R into equation (41) gives(59)EH=γSzI,zIIZII−Re−TcI1TII,Tc+be−Tc−e−TII.


## 3.1. Settlement of the Embankment during Construction

### 3.1.1. Degree of Settlement of the Embankment at HeightzI during Construction

The degree of settlement of the embankment at heightzI during construction can be described as follows:(27)SzI=∫0zIΔεTI−T,Tc−Tdz.Δ
ε can be obtained from equation (21):(28)ΔεEHγ=ZI−Re−TcI1z,Tc−I1z,TI−e−TII1z,TI−I1z,T,where(29)I1z,TI=∫0TIeT′dzdT′dT′,(30)I1z,Tc=∫0TceT′dzdT′dT′=∫0TIeT′dzdT′dT′+∫TITceT′dzdT′dT′.The integrals listed above must be calculated by numerical methods. In order to do this, the loading curve is divided intoN intervals; from T=0 up to T=TI, the curve is divided into I intervals, and from T=TI up to T=Tc, the curve is divided into N−I intervals.Substituting equation (28) into equation (27), the settlement of a point at height zI during construction is obtained from(31)SzIEHγ=ZIzI−Ra−e−TII2z,TI,where(32)a=e−TczII1TI,Tc+I2z,TI,(33)I1TI,Tc=∫TITceTdzdTdT=∑i=INeTiΔzi,(34)zi=zI+1−zI−12,zI−1=zI,zN+1=zN,(35)I2z,TI=∫0TIzeTdzdTdT=∑i=1IzieTiΔzi,zI+1=zN.If it is assumed that the height of the embankment increases linearly with time, and in order to simplify the deduction process, we assumed that the slope of the loading curve is 1, and then equations (33) and (35) will become(36)I1TI,Tc=∫TITceTdzdTdT=∫TITceTdT=eTc−eTI,(37)I2z,TI=∫0TIzeTdzdTdT=∫0TITeTdT=TIeTI−eTI+1.Substituting equations (36) and (37) into equation (31) gives(38)SzIEHγ=ZIzI−R1+e−H−e−ZI−e−zI.


### 3.1.2. Settlement of a Finite Layer (zI,zII) during Construction

The strain from time factorTII up to Tc is(39)Δε=ΔεTII−T,Tc−T.SubstitutingTI by TII and equation (39) into equation (21) gives(40)ΔεEHγ=H−zII−Re−TcI1z,Tc−I1z,T−e−TIII1z,TII−I1z,T.The compression of a finite layer is(41)SzI,zII=∫zIzIIΔεdz.Substituting equation (40) into equation (41) gives(42)SzI,zIIEHγ=H−zIIΔz−Re−TcΔzI1TII,Tc+bΔze−Tc−e−TII,where(43)Δz=zII−zI,(44)b=I2TI,TII−zII1TI,TIIΔz,(45)I2TI,TII=∫TITIIzeTdzdTdT=∑i=Ii=IIzieTiΔzi,(46)I1TI,TII=∫TITIIeTdzdTdT=∑i=Ii=IIeTiΔzi,(47)Δzi=zi+1−zi−12,zII+1=zII,zII+1=zII.


## 3.1.1. Degree of Settlement of the Embankment at HeightzI during Construction

The degree of settlement of the embankment at heightzI during construction can be described as follows:(27)SzI=∫0zIΔεTI−T,Tc−Tdz.Δ
ε can be obtained from equation (21):(28)ΔεEHγ=ZI−Re−TcI1z,Tc−I1z,TI−e−TII1z,TI−I1z,T,where(29)I1z,TI=∫0TIeT′dzdT′dT′,(30)I1z,Tc=∫0TceT′dzdT′dT′=∫0TIeT′dzdT′dT′+∫TITceT′dzdT′dT′.The integrals listed above must be calculated by numerical methods. In order to do this, the loading curve is divided intoN intervals; from T=0 up to T=TI, the curve is divided into I intervals, and from T=TI up to T=Tc, the curve is divided into N−I intervals.Substituting equation (28) into equation (27), the settlement of a point at height zI during construction is obtained from(31)SzIEHγ=ZIzI−Ra−e−TII2z,TI,where(32)a=e−TczII1TI,Tc+I2z,TI,(33)I1TI,Tc=∫TITceTdzdTdT=∑i=INeTiΔzi,(34)zi=zI+1−zI−12,zI−1=zI,zN+1=zN,(35)I2z,TI=∫0TIzeTdzdTdT=∑i=1IzieTiΔzi,zI+1=zN.If it is assumed that the height of the embankment increases linearly with time, and in order to simplify the deduction process, we assumed that the slope of the loading curve is 1, and then equations (33) and (35) will become(36)I1TI,Tc=∫TITceTdzdTdT=∫TITceTdT=eTc−eTI,(37)I2z,TI=∫0TIzeTdzdTdT=∫0TITeTdT=TIeTI−eTI+1.Substituting equations (36) and (37) into equation (31) gives(38)SzIEHγ=ZIzI−R1+e−H−e−ZI−e−zI.


## 3.1.2. Settlement of a Finite Layer (zI,zII) during Construction

The strain from time factorTII up to Tc is(39)Δε=ΔεTII−T,Tc−T.SubstitutingTI by TII and equation (39) into equation (21) gives(40)ΔεEHγ=H−zII−Re−TcI1z,Tc−I1z,T−e−TIII1z,TII−I1z,T.The compression of a finite layer is(41)SzI,zII=∫zIzIIΔεdz.Substituting equation (40) into equation (41) gives(42)SzI,zIIEHγ=H−zIIΔz−Re−TcΔzI1TII,Tc+bΔze−Tc−e−TII,where(43)Δz=zII−zI,(44)b=I2TI,TII−zII1TI,TIIΔz,(45)I2TI,TII=∫TITIIzeTdzdTdT=∑i=Ii=IIzieTiΔzi,(46)I1TI,TII=∫TITIIeTdzdTdT=∑i=Ii=IIeTiΔzi,(47)Δzi=zi+1−zi−12,zII+1=zII,zII+1=zII.


## 3.2. Settlement of the Embankment after Construction

### 3.2.1. Settlement of an Embankment at HeightzI after Construction

The settlement of an embankment at heightzI after construction can be described as follows:(48)ΔSzI=∫0zIΔεTc,Tc+DTdz.Substituting equation (26) into equation (48) gives(49)ΔSzIEHγR=a1−e−DT,where the value of a is equal to equation (32).At the top of the embankment,(50)I1Tc,Tc=0,(51)ΔSHEHγR=e−Tc1−e−DTI2z,Tc.If the height of the embankment varies linearly with time, with a slope of 1, then equations (49) and (51) will become(52)ΔSzIEHγR=zI−e−zI+e−H1−e−DT,(53)ΔSHEHγR=H+e−H−11−e−DT.


### 3.2.2. Settlement of a Finite LayerzI,zII after Construction

The strain of a finite layerzI,zII after construction is(54)Δε=ΔSzI,zIIΔz=ΔSzII−ΔSzIΔz,The settlement of a finite layerzI,zII after construction is(55)ΔSzI,zIIEHγR=e−Tc1−e−DTΔzI1TII,Tc+b,where(56)I1TII,Tc=∫TIITceTdzdTdT=∑i=IIi=NeTiΔzi,where the value of b is equal to equation (44).


## 3.2.1. Settlement of an Embankment at HeightzI after Construction

The settlement of an embankment at heightzI after construction can be described as follows:(48)ΔSzI=∫0zIΔεTc,Tc+DTdz.Substituting equation (26) into equation (48) gives(49)ΔSzIEHγR=a1−e−DT,where the value of a is equal to equation (32).At the top of the embankment,(50)I1Tc,Tc=0,(51)ΔSHEHγR=e−Tc1−e−DTI2z,Tc.If the height of the embankment varies linearly with time, with a slope of 1, then equations (49) and (51) will become(52)ΔSzIEHγR=zI−e−zI+e−H1−e−DT,(53)ΔSHEHγR=H+e−H−11−e−DT.


## 3.2.2. Settlement of a Finite LayerzI,zII after Construction

The strain of a finite layerzI,zII after construction is(54)Δε=ΔSzI,zIIΔz=ΔSzII−ΔSzIΔz,The settlement of a finite layerzI,zII after construction is(55)ΔSzI,zIIEHγR=e−Tc1−e−DTΔzI1TII,Tc+b,where(56)I1TII,Tc=∫TIITceTdzdTdT=∑i=IIi=NeTiΔzi,where the value of b is equal to equation (44).


## 3.3. Parameter Determination

The three parameters of the model are the elastic modulusEH, the creep ratio R, and the time factor Tr.If settlement is measured from different heights of the embankment, from equation (55), the time factor Tr is obtained from(57)ΔSTc,T2ΔSTc,T1=1−e−DT21−e−DT1=1−e−Δt2/Tr1−e−Δt1/Tr,where Δt is the time increment after construction.The solution of equation (57) with different time increments will permit an adjustment for Tr.Substituting equation (55) into equation (42), the values of R from different layers are obtained:(58)R=ZIIe−TcI1TII,Tc+bSzI,zII/ΔSzI,zII1−e−DT+1−be−TII.Substituting the value ofTr and R into equation (41) gives(59)EH=γSzI,zIIZII−Re−TcI1TII,Tc+be−Tc−e−TII.


## 4. Practical Verification

For the Lanzhou–Yongjing highway, a section from 24 + 160 km to 25 + 838 km was completed with a high-fill subgrade. The settlement meters were set up from 24 + 400 km to 24 + 405 km, which was semifilling and semiexcavating subgrades (Figure5(a)). On the right side of the subgrade, a steep slope was there. However, flood drainage was on the left side of the subgrade. The materials excavated from the slope were used to fill the subgrade. The maximum fill height of this section was 18.3 m. From 0 to 15 m height, the subgrade was filled with loess. From 15–16.1 m height, it was filled with a sandy pebble. From 16.1–18.3 m height, it was filled with gravel. The physical and mechanical properties of the subsoils are shown in Table 1. Single-point settlement meters were adopted to monitor the long-term settlement of the subgrade [29]. It can be seen that for both the road shoulder and driveway, 6 layers of settlement meters were embedded. The layout of settlement meters is shown in Figure 5, which were 3 m, 6 m, 9 m, 12 m, 15 m, and 18.3 m. And the embedment depth is reflected in Figure 5(b). Almost two-year settlement data are obtained from the embankment after construction. And the field-measurement settlement data during construction is obtained from [28], which monitored the settlement of each layer of the filling during the construction, and the relationship between filling height and settlement of the filling during construction is represented in Figure 6. And from the monitored data during construction in Figure 6, firstly, the value of Tr is obtained from equation (57), secondly, substituting the value of Tr and all the monitored parameters into equation (58) to calculate the value of R, thirdly, the values of Tr and R are substituted into equation (59) to get the value of EH, and finally, the creep settlement of every layer of the embankment is calculated from equations (42) and (55), respectively, which are presented in Figures 7 and 8. From the field-measured data, the Tr, R, and EH parameters for every layer were established as indicated in Table 2.Figure 5
Single-point settlement meters. (a) Layout profile. (b) Planar graph (reproduced from Jia et al. [28] (under the Creative Commons Attribution License/public domain)).
(a)
(b)Table 1
Physical and mechanical parameters of subsoil (reproduced from Jia et al. [28] (under the Creative Commons Attribution License/public domain)).
Soil type
Thickness (m)
Moisture content (%)
Gravity (kN·m3)
Void ratio
Compression coefficient (MPa−1)
Compression modulus (MPa)
Bearing capacity eigenvalues (kPa)
Loess
10.4
8.5
14.1
0.96
0.21
6.2
150
Loess
4.8
16.5
14.9
0.89
0.23
7.6
180
Pebble
3.3
21
24.8
1.15
—
—
400
Mudstone
—
—
25.3
—
—
—
600Figure 6
Relationship between filling height and settlement of the filling during construction.Figure 7
Comparison between observed and calculated settlements during construction.Figure 8
Comparison between observed and calculated settlements after construction.Table 2
Creep parameters of every layer.
Fill stratum
T
r
E
H
R
3 m∼6 m
172.8312
0.0373
0.0072
6 m∼9 m
172.8342
0.0378
0.0103
9 m∼12 m
172.8387
0.0370
0.0174
12 m∼15 m
172.8377
0.0368
0.0385Settlement values during construction from field measurement and settlement values calculated from equation (31) are presented in Figure 6, which is basically consistent with the experimental finding [30].The measured and calculated settlements after construction obtained with equation (49) are presented in Figure 7.The agreement between the measured and calculated settlements is excellent, although the measured settlement values are more irregular because this area is in a soil region that is seasonally frozen, and during the winter, soil freezing leads to the ground heave [31, 32].


## 5. Conclusions

To describe the time-dependent behavior of clay inside a high-filled embankment, the viscoelastic model of Poynting–Thomson was applied in this study in order to calculate the creep settlement of high-filled embankments. And the stress history is reflected by the hereditary integral. One-dimensional compression was assumed in this investigation. The settlement values at a point and of a finite layer during construction and after construction are obtained. The three parameters of this model,Tr, R, and EH, can be adjusted for the in situ monitoring data from every layer. By comparison of calculated results obtained by the proposed method with the field monitoring data, it was concluded that the proposed method can give a precise prediction of the creep settlement of high-filled embankments. In the future work, the other viscoelastic models which can better reflect the time-dependent creep behavior of soils will apply to amplify understanding of settlement behavior in more complex soil conditions.



---

*Source: 4627174-2019-07-28.xml*