International Journal of Engineering and Applied Sciences (IJEAS)

ISSN: 2394-3661, Volume-3, Issue-2, February 2016

2D FEM Analysis of Slippage Phenomenon in Earth

and Rockfill Dams
Nishtha Saraswat

Abstract - Dams are gaining more attention in recent years

II. SLIPPAGE PHENOMENON
due to the rise of the environmental awareness and renewable It is observed that the peak strains which remain uniform
energy and sustainability concepts. Earth embankment dams within each constituent material change abruptly at the
are preferred over gravity dams for the ease of construction and core-shell interface. As the core consists of soft clay material
economical advantage. Despite considerable advances in the the settlements in core are greater than those in shell. In reality
field of geotechnical engineering, the occurrence of frequent
occurrence of sliding takes place between shell (rockfill) and
destructive earthquakes during the past decade and the trend
towards construction of structures of unprecedented size and of core (clay) is likely to occur. The results obtained considering
novel designs has enhanced the importance of earthquake this phenomenon are more accurate as compared to the ones
resistant structural design. Rock-fill dam is a type of earth dam which does not take into consideration the relative movement
where a compacted central clay core is supported from the rock but treats the whole dam body as a single structure. The
shells by a series of transition zones built of properly graded contact element feature of ANSYS 16.2 is used to model the
material. In recent years, rock-fill dams, especially the interface. The coefficient of friction between dry rockfill and
impervious-faced rock-fill dams (IFRD), are built all around the dry clay (D/S section interface) is 0.4 and between wet
world using asphalt or concrete as the impervious material in rockfill and wet clay (U/S section) is 0.2. The coefficient of
the upstream face of the dam. This paper analyses the slippage
friction is less in saturated conditions.
phenomenon in the core-shell interface of earth and rock-fill
dams with the aid of ANSYS16.2 package. The software was
first validated then used to examine the slippage in Tehri Dam, III. 2D FEM PLAIN STRAIN MODEL
which is located in the seismically active region of Himalayas. A As the conditions of plain strains are satisfied in 2D fem
2D fem analysis is done by modelling the dam as a linear, elastic, analysis: The cross-section of the dam was modeled using a
non-homogenous material. The slippage phenomenon is
commercially available finite element package, ANSYS 16.2
modelled using contact element feature of ANSYS. Among the
variables, the non-homogeneity of the core and shell material according to ANSYS users manual. The natural frequencies
and coefficient of friction is considered in the analysis. On and mode shapes of the 2D dam model are obtained by modal
evaluation it is found that core settlement increases where analysis. The element type used is Solid 8 node 183 plane
core-shell influence has been considered. strain solid elements which is an 8 noded structural shell,
suitable for analyzing thin to moderately thick structures. The
Index TermsCoefficient of friction, core-shell interface, element has 8 nodes with 6 degrees of freedom at each node.
Earth & Rockfill dam, settlement, slippage. The whole domain is divided into 8 x 8 meshes for all the
cases. The boundary condition is given fixed at the bottom of
I. INTRODUCTION the Dam. This condition closely resembled the field situation.
Tehri Dam is a composite earth & rock fill dam of height A total number of 1690 elements are considered. Material
260.5 m consisting of shell, riprap, filter and core constructed properties as obtained after testing are given below in Table 1.
on the Bhagirathi river in lesser Himalayas. The crest length is
570m along axis. The cross section of dam is as shown in Table1. Material Properties
Fig.1. The crest width is 20 m, flared to 25.0 m at the S.No Property Shell Core
abutments. It consists of central impervious core and a shell of
pervious nature consisting of well graded rock fill. The core 1. Shear Modulus ,G 5.883e+07 2.942e+07
of the dam is slightly inclined whereas the dam has an (N/m2)
upstream slope of 2.5:1 and downstream slope of 2:1.
2. Modulus of 1.5e+09 5.29e+08
Elasticity E (N/m2)

3. Poisons ratio 0.35 0.3

4. Moist 2.402e+04 1.961e+04
Density(N/m3)

5. Saturated 2.442e+04 2.108e+04

Density(N/m3)
Fig. 1 Cross-section of Dam
6. Dry Density(N/m3) 2.314e+04 1.814e+04

Nishtha Saraswat, Department of Earthquake Engineering, Indian

Institute of Technology, Roorkee, India, +91-70604-67059

66 www.ijeas.org
 2D FEM Analysis of Slippage Phenomenon in Earth And Rockfill Dams

VI. FREE VIBRATION ANALYSIS

The frequency analysis has been carried out in ANSYS 16.2
using the block lanczos method technique to study the
fundamental frequency. The Block Lanczos Method is a very
efficient algorithm to perform a modal analysis for large
models. It is a fast and robust algorithm and used for most
applications as the default solver. The fundamental vibration
occurs in lateral translational mode whereas the second mode
Fig.2 Discretization of 2D model is vertical translation mode. The translations are pronounced
much more near the crest of the dam. The fundamental natural
In static analysis, the submerges unit weight is taken for the periods of the model in this research have been calculated to
U/S shell elements are the water will be seeping in from the be 1.670648 s for model without considering the shell-core
end, saturated unit weight for the core and moist density to be interface. Considering a contact element between the shell
considered for the D/S shell elements. and core the fundamental natural period increased to 1.8992 s.

IV. MODELING OF SLIPPAGE VII. SUMMARY

Surface-surface contact has been modeled to simulate the Initial static analysis and modal analysis is being done and
core-shell slippage phenomenon. These contact elements use results are confirmed with the previous established results
a "target surface" and a "contact surface" to form a contact done on seismic analysis of Tehri dam. Further a comparative
pair. study between the dam modelled with considering the
The target surface is modeled with either TARGE169 (for slippage between core-shell interface is done and the results
2-D) (shell) of horizontal and vertical displacement, natural frequency are
The contact surface is modeled with elements, analyzed with respect to the dam in which core and shell
CONTA172.(core) elements are considered to be glued.
Contact surface has been created on both upstream and Parameter Contact Elements
downstream of core-shell interface. The various types of (Maximum values)
Used Not used
contact options available are: since the elements are not tied
X-Displacement (m) 2.456 2.388
to each other and free to slide over each other standard type of
contact is selected. Y-Displacement (cm) 6.044 5.852
Vertical - Tensile 0.77 0.65
Stress (MPa) -Compressive 8.05 8.10
Horizontal -Tensile 1.27 1.38
Stress (MPa)
-Compressive 3.57 3.51
Table 2: Results of Tehri dam

Fig. 3. Contact Elements between core-shell interfaces

V. INITIAL STATIC ANALYSIS

The first step in the analysis is to establish the in situ stress
conditions that exist before the earthquake occurs. The two
most important soil properties required for the initial static Figure 4: Vertical displacement profile along vertical
analysis are the:
Mass Density
Poissons Ratio of the materials used
Initial static analysis of the rock-fill dam is carried out due to
the self-weight load at the end of the construction. It can be
observed that there is a trend of increasing displacement from
base towards the crest until some depth below the crest it is
observed where maximum settlement has occurred and it
again decreased for the portion lying above it. On considering
the slippage phenomenon profile pattern of similar nature is
observed with higher magnitude.
Figure 5: Vertical displacement profile along horizontal

67 www.ijeas.org
 International Journal of Engineering and Applied Sciences (IJEAS)
ISSN: 2394-3661, Volume-3, Issue-2, February 2016
VIII. CONCLUSIONS
1. It is seen that near the crest of the dam, the maximum
horizontal and vertical displacement values are
more, than in the rest of the dam. This is in
agreement with the common knowledge that the top
one-third height of the dam vibrates more than other
part, and almost all displacement occurs in the top
part. The Horizontal and Vertical displacements are
within permissible limits as per I.S.:8826-1978 (1%
of the embankment height).
2. Under the assumed coefficient friction based on
various literature surveys, it is found that the core
settlement is considerably affected by the influence
of core-shell interface with the increase in excess
settlement of about 0.020-0.040 m and also
participates in the increase of natural time period.
3. The increase in the crest displacement is of the order
of 0.070 m in horizontal direction.
4. Vertical stresses have increased from 0.65 MPa to
0.77 MPa., increase in horizontal stresses is of the
order of 0.1MPa.

REFERENCES
[1] ANSYS Mechanical APDL Structural Analysis Guide
[2] Basudhar, P.K. Rao, N.S.V.K., Bhookya, M. and Dey, A. (2010), 2D
Fem Analysis of Earth and Rockfill Dams under Seismic Condition,
Fifth International Conference on Recent Advances in Geotechnical
Earthquake Engineering and Soil Dynamics and Symposium in Honor
of Professor I.M, Idriss, San Diego, California.
[3] Paul, D.K.(2000). Seismic Safety Analysis of High Rock-Fill Dam
Subjected to Severe Earthquake Motion, 1164, 12 WCEE, Auckland,
New Zealand
[4] Shakya Yuman & Suwal Rajan(2014),Sesimic Deformation Analysis
of Rockfill DamJournal of the Institute of Engineering Vol10. N0.1,
pp.199-204.

Nishtha Saraswat is currently a Masters

student in the Department of Earthquake Engineering of I.I.T Roorkee
(India) with bachelors in Civil Engineering from Harcourt Butler
Technological Institute , Kanpur( India).The Areas of interest are
Geotechnical Engineering, Structural Designing ,Soil Dynamics and slope
stability analysis. The dissertation topic of final year research is Performance
Based Design of Earth and Embankment Dams in Seismic Zone. The
research work done in bachelors was on the topic Planning and Design of '
Grand Colosseum' cricket stadium in Lucknow involving use of STAAD
PRO and AUTOCAD. Manually checking of designs in accordance with
Indian Design Codes. The internships undertaken by her were on the topics
of Analysis and Design of Large Span Steel Roof Truss as a part of summer
training at CSIR-CBRI, Roorkee. Further she also has practical experience of
onsite work while working on the project Constructional Aspects of
Modern High Rise Green Building (GAIL Jubilee Tower). She also did a
summer intern in L & T MHI Boilers Private Limited which involved
Familiarization to Structural Design of Super Critical Boiler Structure.
She has an excellent academic record and has always been top performers in
her batch during undergrad and post grad years.

68 www.ijeas.org