Impact of soil deformation on phreatic line in earth-fill dams
Highlights
► Unsaturated seepage analysis with and without coupling deformation is conducted. ► A simplified approach for consolidation analysis of unsaturated soil is employed. ► Impact of soil deformation on phreatic line in earth-fill dam is investigated.
Introduction
When a dam is subjected to varied water levels, its seepage analysis and stress–strain analysis are usually conducted separately. Seepage analysis is firstly carried out and subsequently followed by stress–strain analysis. Based on such an approach, the impact of stress–strain properties of soil on seepage is neglected. Seepage analysis is a classic topic in soil mechanics and it can be conducted by using many numerical methods such as the routine finite difference method, the finite-volume method, the boundary-fitted coordinate transformation method, the finite element method, the numerical manifold method, the meshless method etc. (Bathe and Khoshgoftaar, 1979, Darbandi et al., 2007, Desai, 1976, Jiang et al., 2010, Jie et al., 2004, Lam and Fredlund, 1984, Li et al., 2003, Zheng et al., 2005). The soil mass below the phreatic line is under saturated conditions with positive pore water pressure, then seepage theories for saturated soil are applicable. Regarding the soil mass above the phreatic line with negative pore water pressure, seepage theories for unsaturated soil are required and the coefficient of permeability varies with negative pore water pressure (Fredlund and Rahardjo, 1993).
If the variation in water level is small, it is generally believed that the seepage is marginally affected by the stress–strain properties of soil. However, when the water level varies greatly, the impact of stress–strain properties of soil on seepage cannot be neglected and consolidation theories are required. Biot's consolidation theory has been extensively used in the analysis of saturated soil (Biot, 1941, Sandhu and Wilson, 1969). However, if negative pore water pressure exits, the consolidation theories for unsaturated soil will be more suitable.
The consolidation model coupling deformation, pore water pressure and pore air pressure was first proposed by Barden (1965). Closed formulations were derived by using continuity equations of water and gas, Darcy's law, suction state function, Bishop's effective stress equation and the relationship between porosity and effective stress. Other typical consolidation formulations were proposed by Scott (1963), Lloret and Alonso (1980), and Fredlund et al. (Fredlund and Hasan, 1979, Fredlund and Morgenstern, 1976, Fredlund and Rahardjo, 1993).
In this paper, a simplified approach for consolidation analysis of unsaturated soil suggested by Shen (2003) is used to conduct consolidation analysis of an earth-rock dam subjected to water filling. This approach is based on Bishop's effective stress (Bishop, 1959). By introducing the air drainage ratio, pore air pressure can be solved indirectly and is no longer treated as an unknown quantity in governing equations, greatly simplifying the amount of computation and the complexity of programming. It has been successfully used to analyze surface cracks on clay by Deng and Shen (Deng et al., 2003, Deng et al., 2006, Deng and Shen, 2006). Here, this approach is employed to analyze the seepage in an earth-rock dam during water filling and to study the impact of deformation of soil on the development of phreatic line.
Section snippets
Governing equations for consolidation analysis
To be consistent with the general mechanical analysis, sign convention used in elasticity mechanics is adopted in this section unless otherwise stated. Such sign convention is opposite to that generally used in soil mechanics. In order to assure the value of suction positive, suction is defined as s=uw−ua, different from the conventional definition s=ua−uw, where ua is the pore air pressure and uw is the pore water pressure.
Bishop's effective stress is adopted with its definition as follows:
Results
The simplified consolidation approach for unsaturated soils better corresponds to reality than that for saturated soils regarding numerical simulations of earth-rock dams. At the same time, computational complexity is not increased too much. So it is very feasible in carrying out numerical analysis of earth-rock dams.
Fig. 1 shows the mesh of the cross-section of a high core-wall dam. The bottom elevation of the mesh is 142 m and the top elevation is 283 m. Upstream side is on the left and the
Conclusions
In this study, a simplified approach for consolidation analysis of unsaturated soil is applied to numerical simulation of an earth-rock dam during the process of water-filling. The computational results include stress and displacement fields within the dam and the variations of pore water pressure and phreatic line. The results show that due to the coupling effect between deformation and pore water pressure, the development of pore water pressure in the core-wall of the dam is quicker than that
Acknowledgments
The supports of Natural Science Foundation of China (51039003), National Basic Research Program of China (973 Program 2010CB732103), and the State Key Laboratory of Hydroscience and Engineering (2012-KY-02) are gratefully acknowledged.
References (25)
- et al.
Modeling unconfined seepage flow using three-dimensional numerical manifold method
Journal of Hydrodynamics
(2010) - et al.
Seepage analysis based on boundary-fitted coordinate transformation method
Computers and Geotechnics
(2004) - et al.
Saturated–unsaturated transient finite element seepage model for geotechnical engineering
Advances in Water Resources
(1984) - et al.
Free surface seepage analysis based on the element-free method
Mechanical Research Communication
(2003) Consolidation of compacted and unsaturated clays
Geotechnique
(1965)- et al.
Finite element free surface seepage analysis without mesh iteration
International Journal for Numerical and Analytical Methods in Geomechnics
(1979) General theory of three-dimensional consolidation
Journal of Applied Physics
(1941)The principle of effective stress
Teknisk Ukeblad
(1959)- et al.
A moving-mesh finite-volume method to solve free-surface seepage problem in arbitrary geometries
International Journal for Numerical and Analytical Methods in Geomechnics
(2007) - et al.
Numerical simulation of crack formation process in clays during drying and wetting
Geomechanics and Geoengineering
(2006)
Numerical simulation of crack formation due to desiccation in the clay surface
Chinese Journal of Geotechnical Engineering
Cited by (9)
Treatment of transitional element with the Monte Carlo method for FEM-based seepage analysis
2013, Computers and GeotechnicsCitation Excerpt :Not surprisingly, a similar application is to conduct consolidation analysis of unsaturated soil. In fact, the numerical simulation of water-filling in an earth-rock dam reported in Ref. [27] has actually employed the Monte Carlo method to deal with the transitional element. The Monte Carlo method is used in this paper to improve the convergence of seepage analysis.
Characteristics and mechanism of inhomogeneous pore pressure of core wall in super-high rockfill dam
2022, Shuili Xuebao/Journal of Hydraulic EngineeringMultifield Coupling Numerical Simulation of the Seepage and Stability of Embankment Dams on Deep Overburden Layers
2022, Arabian Journal for Science and EngineeringAnalyses on finite earth pressure and slope safety factors
2019, Qinghua Daxue Xuebao/Journal of Tsinghua UniversityEvaluation of seismic performance of earth dams due to the level of its reservoir using finits element method
2018, International Conference on Civil, Structural and Transportation Engineering