# Application of Effective Stress Model to Analysis of Liquefaction and Seismic Performance of an Earth Dam in China

Authors: Changqing Qi; Wei Lu; Jimin Wu; Xing Liu
Journal: Mathematical Problems in Engineering
 (2015)

---

## Abstract

Earthquake-induced liquefaction is one of the major causes of catastrophic earth dam failure. In order to assess the liquefaction potential and analyze the seismic performance of an earth dam in Fujian, Southeastern China, the in situ shear wave velocity test was firstly carried out. Results indicate that the gravelly filling is a type of liquefiable soil at present seismic setting. Then the effective stress model was adopted to thoroughly simulate the response of the soil to a proposed earthquake. Numerical result generally coincides with that of the empirical judgment based on in situ test. Negative excess pore pressure developed in the upper part of the saturated gravelly filling and positive excess pore pressure developed in the lower part. The excess pore pressure ratio increases with depth until it reaches a maximum value of 0.45. The displacement of the saturated gravelly soil is relatively small and tolerable. Results show that the saturated gravelly filling cannot reach a fully liquefied state. The dam is overall stable under the proposed earthquake.


---

## Body

## 1. Introduction

Earth dams usually have better seismic performance during earthquakes. According to statistics, seldom earth dams have been totally out of service after earthquakes in the past few decades in China [1]. But when the dam contains or is situated on liquefiable materials, earthquake-induced liquefaction may cause considerable reduction in stiffness and strength of soil, resulting in dam failure [2]. A number of dam failures or damages have been reported due to seismically induced liquefaction. The most classic example is the lower San Fernando dam during the 1971 San Fernando earthquake. Liquefaction induced flow slide on the upstream side of the dam nearly caused the dam to be out of service [3]. The liquefaction slide of Baihe dam of Miyun Reservoir during the 1976 Tangshan earthquake was another representative example. The estimated 0.15 million m3 volumetric slide in the upstream part of the dam aroused great panic at that time [4]. Several dam failures in Chilean [5], Japan [6], and India [7] were also reported causing great damage.In sight of its immense economic damage and loss of life, earth dam failures due to liquefaction have drawn great concern in the past half century. Several types of approaches have been developed to study this problem. The empirical relationships based on tested indexes are commonly used methods at present [8–10]. The empirical relationship methods can give an overall and quick judgment on liquefaction potential of gravelly soil. But for thoroughly understanding the response of coarse soil to cyclic shear loads, numerical modeling technique is better [11–13]. The seismic response of soil has a direct relation to the progressive build-up of pore pressure during an earthquake. The increasing pore pressure and decreasing effective stress control the resistance of the soil to deformation [14]. Thus, assessment on progressive degradation of soil strength is important in dynamic soil liquefaction analysis.The effective stress method is a useful tool in modeling the progressive loss of soil strength caused by development of pore pressure. Liyanapathirana and Poulos [15] summarized four main categories of liquefaction models based on effective stress analysis method, which are (1) models based on plasticity theory; (2) stress path methods; (3) correlations between pore pressure response and plastic volume change tendency; and (4) use of experimentally observed undrained pore pressure response. The model of the third category, which is often referred to as the Finn model, is explicit and has a lesser number of parameters. The model can take into account stiffness and strength degradation due to pore pressure development and was used in this paper for dam liquefaction and seismic performance assessment.


## 2. Methodology

The strength and stiffness of soil are primarily governed by effective stress, and so it is desirable to evaluate seismic response of soil in terms of effective stress. For saturated granular material, adopted numerical model should reflect the variation of pore pressure and thus can predict the effective stress level.For saturated soil under undrained conditions, Martin et al. [16] suggested that the pore pressure increment is related to the change of plastic volumetric strain of the soil skeleton:(1)Δu=-1CbΔεvd,where Δu is pore pressure increment, Cb is bulk compressibility, and Δεvd is the plastic volumetric strain increment. The referred model can adequately reflect the response of pore pressure until liquefaction triggering point [17] and thus is an effective tool for assessment of soil liquefaction risk.The plastic volumetric strain increment could be obtained by various constitutive theories. In this paper, the expression presented by Byrne [18] was adopted. The plastic volumetric strain increment was expressed as follows:(2)Δεvdγ=C1exp⁡-C2εvdγ,where γ is the shear strain in the current cycle, εvd is the accumulated volumetric strain from prior cycles, and C1 and C2 are constants that depend on the relative density Dr:(3)C1=7600Dr-2.5,C2=0.4C1.In this paper, the model was incorporated into the finite difference computer programFlac3D to perform a nonlinear fully coupled dynamic analysis. Flac3D is based on a continuum finite difference discretization using the Lagrangian approach [19]. The equations of motion are utilized to obtain the velocities and displacements when dynamic load is excreted. The equation of motion can be expressed as(4)ρ∂vi∂t=∂σij∂xi+ρbi,where ρ is material density, t is time domain, xi is coordinate vector, σij is stress tensor, and bi is body force.For a fully nonlinear method, any given function can be used in dynamic analysis ofFlac3D. The general constitutive equation is expressed as(5)σij˙=Mσij,εij˙,κ,where σij˙ is stress rate tensor, εij˙ is strain rate tensor, κ is the parameter taking the loading history into account, and M is the given functional expression.The Mohr-Coulomb elastic-perfectly plastic constitutive relation is the commonly used constitutive models for soil. In order to adapt the model for dynamic analysis, several modifications have been achieved. In this paper, the modification defined by Puebla et al. [20] was used. The secant shear modulus G and bulk modulus B were considered to be stress-dependent and given as follows:(6)G=kg·Pa·σm′/Pan,B=kb·Pa·σm′/Pam,where kg and kb are shear and bulk modulus numbers, n and m are modulus exponents, σm′ is the mean effective stress, and Pa is atmospheric pressure.


## 3. Statement of Dam Conditions

Dongzhen reservoir, located about 6.0 km upstream Putian City in Fujian Province, Southeastern China, has a normal storage capacity of 435 million m3. The reservoir has a comprehensive function of flood control, irrigation, and power generation. The water retaining dam has an irregular geometry (Figure 1). The longitudinal profile of the dam has an asymmetric U-shape which is steep on the right flank and gentle on the left (Figure 2(a)). The dam is a core-wall earth dam with a maximum height of 58.6 m. The normal dammed water level is about 8.1 m to the dam crest (Figure 2(b)). The core wall made of lean clay and the filling is gravelly soil. A layer of rock blocks revetment was placed on the surface to protect the slope. The thickness of the blocks revetment on the upstream surface is about 2.0 m.Figure 1
Plane view of Dongzhen Reservoir dam.Figure 2
Longitudinal (a) and cross (b) section of Dongzhen Reservoir dam. A1–A8 denote the location for excess pore pressure ratio analysis.
(a)
(b)The dam was built in the 1950s, when there was no seismic design standard in China. The study area is located in the coast of Taiwan Channel. Caused by the movement of the Pacific Plate, a series of NNE and NW fault zones can be found in this area (Figure3). The area is currently within the region with seismic intensity of 7 based on Chinese Standard [21]. According to the seismic safety analysis completed by Fujian Provincial Institute of Geological Engineering Investigation, the equivalent earthquake magnitude which affects the dam area is Ms 6.5 and the epicenter distance is about 50 km. The peak bedrock acceleration of the maximum credible earthquake is about 158 gal. Owing to the high public interest and the large amount of potentially affected persons, questions about the seismic stability of the dam are important to study in detail.Figure 3
Geological map of the study area.


## 4. Shear Wave Velocity Test

In order to study the seismic stability of the dam, the first step is to evaluate the liquefaction potential of gravelly soils [11]. The usually available field test methods are the standard penetration test (SPT), the cone penetration test (CPT), in situ shear wave velocity measurement (Vs), and the Becker penetration test (BPT) [22]. Among these approaches, the shear wave velocity measurement is more feasible in gravelly soil field because this type of soil is difficult to penetrate [10, 22, 23]. A three-component geophone setup was placed in a 30 m deep borehole to obtain the shear wave velocity at different depths. According to Chinese Standard [24], the soils can be determined as nonliquefiable when the measured shear wave velocity is greater than calculated limit shear wave velocity. The limit shear wave velocity can be calculated from the following empirical equation:(7)Vst=291·Kh·Z·rd,where Vst is the limit shear wave velocity, Kh is the peak acceleration coefficient, the value is 0.1 for the field with seismic intensity of 7, Z refers to the depth, and rd can be calculated as follows:(8)rd=1.0-0.01ZwhenZ=0~10m,rd=1.1-0.02ZwhenZ=10~20m,rd=0.9-0.01ZwhenZ=20~30m.The measured and calculated results are listed in Table1. It can be concluded that the upstream gravelly filling has the possibility of liquefaction below the depth of 12.0 m.Table 1
Liquefaction judgment using shear wave velocity test.
Number
Depth (m)
Measured shear wave velocity (m/s)
Limit shear wave velocity (m/s)
Primary liquefaction estimation
1
2–4
315
157
Nonliquefiable
2
4–6
292
201
Nonliquefiable
3
6–8
305
235
Nonliquefiable
4
8–10
315
263
Nonliquefiable
5
10–12
296
286
Nonliquefiable
6
12–14
300
304
Liquefiable
7
14–16
285
319
Liquefiable
8
16–18
314
331
Liquefiable
9
18–20
287
340
Liquefiable
10
20–22
302
350
Liquefiable
11
22–24
303
361
Liquefiable
12
24–26
304
371
Liquefiable
13
26–28
304
380
Liquefiable
14
28–30
288
387
Liquefiable


## 5. Liquefaction Assessment and Deformation Analysis

### 5.1. Numerical Modeling and Parameters

The detailed seismic response of Dongzhen Reservoir dam was simulated using softwareFlac3D. A three-dimensional modeling of the dam was established (Figure 4). x-axis of the model was set along the river, y-axis was set perpendicular to the river center line, and the positive z-axis was set upward. A field seismic wave provided by Fujian Provincial Institute of Geological Engineering Investigation was used in this paper for dam seismic response analysis. The seismic wave was selected from historical seismic wave database considering similar field condition and potential influence of the epicenter. The seismic acceleration was recorded from the Ms 6.6 Imperial Valley Earthquake in 1979 at an epicenter distance of 57.6 km. The peak seismic acceleration is 189 gal. The horizontal acceleration time history was shown in Figure 5. The viscous absorb boundary, developed by Lysmer and Kuhlemeyer [25], was used to absorb the unbalanced energy at the boundary. In view of the large permeability difference between the clay and gravelly soil, the clay core wall was considered as impermeable layer. The phreatic water level was used and normal retaining water level of the reservoir was 8.1 m below the dam crest. The hydrostatic pressure caused by the retaining water was applied to the upper stream surface of the dam.Figure 4
Three-dimensional finite difference mesh of Dongzhen Reservoir dam.Figure 5
Input earthquake acceleration record.The analysis was conducted in two steps. First, the static analysis was carried out. The Mohr-Coulomb model with stress-dependent materials properties was used for all the parts of the dam. The materials properties of the dam are given in Table2. In the second step, the effective stress model was applied to the upstream gravelly soil, and the Mohr-Coulomb model was still used for the rest parts. The relative density (Dr) of the upstream gravelly soil used in the Byrne model was obtained by field test and the mean value is 67.5%.Table 2
Properties of the soils used in the numerical modeling.
Material
Unit weight/kN/m3
Shear modulus numbers
Shear modulus exponent
Bulk modulus numbers
Bulk modulus exponents
Cohesion/kPa
Friction angle/°
Permeability coefficient cm/s
Clay
19.5
490
0.50
1470
0.50
48.0
24.0
7.50
×
10
-
6
Gravel
20.0
1200
0.63
3600
0.63
0.0
37.0
5.22
×
10
-
2
Block
22.0
1050
0.69
3150
0.69
0.0
45.0
0.55In order to illustrate the liquefaction degree of the gravelly soil, the excess pore pressure ratioRu was defined and denoted as(9)Ru=Ueσm0′,where Ue is the excess pore pressure during the earthquake and σm0′ is the mean effective stress in the static condition. Ru=1.0 represents a fully liquefied state and Ru=0.0 represents a static condition.


### 5.2. Result of the Analyses

The responses of the saturated gravelly soil in terms of excess pore pressure ratio (Ru) at different depths are presented in Figure 6. The relative positions selected for illustration are shown in Figure 2. The distance of the positions is about 5.0 m and the depth of point A1 is about 2.0 m which is just below the rock block revetment layer. Figure 6 indicates that negative excess pore pressure ratio developed in the upper part of the saturated gravelly filling, which means a decrease of pore pressure and increase of effective stress. This phenomenon coincides with the dynamic behavior research of moderate dense granular material with low confining stress [26]. The increase of effective stress prevents the occurrence of potential liquefaction in this range. This conclusion also coincides with the result of in situ shear wave test. Below the depth of 12.0 m, only positive excess pore pressures were built up during the earthquake. The excess pore pressures ratios increase sharply in the time 7.5 to 10 s, corresponding to the period of strong shaking, and then level off. The permanent mean excess pore pressure ratio increases with depth and then decreases. The maximum permanent mean excess pore pressure ratio is less than 0.45, which means that upstream saturated gravelly filling cannot reach a fully liquefied state.Figure 6
Excess pore pressure ratio of upstream gravelly soil at points A1 to A8. The location of points A1 to A8 is shown in Figure2(b).
(a)
(b)
(c)
(d)
(e)
(f)
(g)
(h)Figure7 illustrates the time history of the maximum horizontal displacement in the upstream gravelly soil of the dam. We can observe that the permanent displacement is relatively small. The maximum horizontal displacement of the upstream filling occurs at middle slope with a value about 7.3 cm. The deformations are tolerable and do not have significant influence on the service function of the whole dam according to Hynes-Griffin and Franklin [27].Figure 7
Time history of maximum horizontal displacement in upstream gravelly soil.


## 5.1. Numerical Modeling and Parameters

The detailed seismic response of Dongzhen Reservoir dam was simulated using softwareFlac3D. A three-dimensional modeling of the dam was established (Figure 4). x-axis of the model was set along the river, y-axis was set perpendicular to the river center line, and the positive z-axis was set upward. A field seismic wave provided by Fujian Provincial Institute of Geological Engineering Investigation was used in this paper for dam seismic response analysis. The seismic wave was selected from historical seismic wave database considering similar field condition and potential influence of the epicenter. The seismic acceleration was recorded from the Ms 6.6 Imperial Valley Earthquake in 1979 at an epicenter distance of 57.6 km. The peak seismic acceleration is 189 gal. The horizontal acceleration time history was shown in Figure 5. The viscous absorb boundary, developed by Lysmer and Kuhlemeyer [25], was used to absorb the unbalanced energy at the boundary. In view of the large permeability difference between the clay and gravelly soil, the clay core wall was considered as impermeable layer. The phreatic water level was used and normal retaining water level of the reservoir was 8.1 m below the dam crest. The hydrostatic pressure caused by the retaining water was applied to the upper stream surface of the dam.Figure 4
Three-dimensional finite difference mesh of Dongzhen Reservoir dam.Figure 5
Input earthquake acceleration record.The analysis was conducted in two steps. First, the static analysis was carried out. The Mohr-Coulomb model with stress-dependent materials properties was used for all the parts of the dam. The materials properties of the dam are given in Table2. In the second step, the effective stress model was applied to the upstream gravelly soil, and the Mohr-Coulomb model was still used for the rest parts. The relative density (Dr) of the upstream gravelly soil used in the Byrne model was obtained by field test and the mean value is 67.5%.Table 2
Properties of the soils used in the numerical modeling.
Material
Unit weight/kN/m3
Shear modulus numbers
Shear modulus exponent
Bulk modulus numbers
Bulk modulus exponents
Cohesion/kPa
Friction angle/°
Permeability coefficient cm/s
Clay
19.5
490
0.50
1470
0.50
48.0
24.0
7.50
×
10
-
6
Gravel
20.0
1200
0.63
3600
0.63
0.0
37.0
5.22
×
10
-
2
Block
22.0
1050
0.69
3150
0.69
0.0
45.0
0.55In order to illustrate the liquefaction degree of the gravelly soil, the excess pore pressure ratioRu was defined and denoted as(9)Ru=Ueσm0′,where Ue is the excess pore pressure during the earthquake and σm0′ is the mean effective stress in the static condition. Ru=1.0 represents a fully liquefied state and Ru=0.0 represents a static condition.


## 5.2. Result of the Analyses

The responses of the saturated gravelly soil in terms of excess pore pressure ratio (Ru) at different depths are presented in Figure 6. The relative positions selected for illustration are shown in Figure 2. The distance of the positions is about 5.0 m and the depth of point A1 is about 2.0 m which is just below the rock block revetment layer. Figure 6 indicates that negative excess pore pressure ratio developed in the upper part of the saturated gravelly filling, which means a decrease of pore pressure and increase of effective stress. This phenomenon coincides with the dynamic behavior research of moderate dense granular material with low confining stress [26]. The increase of effective stress prevents the occurrence of potential liquefaction in this range. This conclusion also coincides with the result of in situ shear wave test. Below the depth of 12.0 m, only positive excess pore pressures were built up during the earthquake. The excess pore pressures ratios increase sharply in the time 7.5 to 10 s, corresponding to the period of strong shaking, and then level off. The permanent mean excess pore pressure ratio increases with depth and then decreases. The maximum permanent mean excess pore pressure ratio is less than 0.45, which means that upstream saturated gravelly filling cannot reach a fully liquefied state.Figure 6
Excess pore pressure ratio of upstream gravelly soil at points A1 to A8. The location of points A1 to A8 is shown in Figure2(b).
(a)
(b)
(c)
(d)
(e)
(f)
(g)
(h)Figure7 illustrates the time history of the maximum horizontal displacement in the upstream gravelly soil of the dam. We can observe that the permanent displacement is relatively small. The maximum horizontal displacement of the upstream filling occurs at middle slope with a value about 7.3 cm. The deformations are tolerable and do not have significant influence on the service function of the whole dam according to Hynes-Griffin and Franklin [27].Figure 7
Time history of maximum horizontal displacement in upstream gravelly soil.


## 6. Conclusion and Discussion

Soil liquefaction resulting from earthquake shaking is a major cause of damage in earth dam engineering. The field of soil liquefaction research is now only semimatured. The generally used liquefaction assessment methods include laboratory test, empirical relationships base on in situ test indexes, and numerical modeling analyses. Laboratory tests sometimes are too complicated and expensive to be used in engineering. Also, undisturbed test samples are difficult to collect and store. These two aspects restrict the wide use of laboratory tests. The empirical judgments only need several in situ test indexes. Empirical relationships were established based on numerous field cases and accordingly can give a reasonable result in liquefaction estimation. Thus, the empirical judgment is the dominated method at present. But the empirical methods can only provide overall seismic response estimation. The history of soil strength and deformation during earthquake cannot be reflected in empirical relationships. Oppositely, numerical modeling can adopt more complex geometrical, stress history and constitutive models and can easily give the results of effective stress, shear strain, and deformation at any time. So it is more feasible for detailed seismic performance analysis.Dongzhen Reservoir water retaining dam is only 6.0 km upstream the Putian City. The stability of the dam has drawn great concerns. The result of in situ shear wave velocity test demonstrates that the gravelly filling has liquefaction potential under an earthquake of intensity 7. In order to thoroughly investigate the seismic response of the dam, three-dimensional finite difference technique was adopted in this paper. The results suggest that negative excess pore pressure developed in the upper part of the saturated gravelly filling, which means that the gravelly soil is not likely to liquefy above the depth of 12.0 m. Below this depth, positive excess pore pressure developed and the maximum excess pore pressure ratio appears at the depth of about 22.0 m. The maximum excess pore pressure ratio rises to about 0.45. This means that the gravelly soil cannot reach a fully liquefied state. The displacement in upstream saturated gravelly filling is relatively small and tolerable. The deformation would not have significant impact on the overall stability of the dam. The gravelly filling still can maintain its seismic resistance.



---

*Source: 404712-2015-07-07.xml*