## Static and dynamic analyses of the heightening of concrete face gravel dam Limon (Peru)

**Authors:**Lyapichev Y.P.**Issue:**Vol 15, No 2 (2019)**Pages:**158-168**Section:**Dynamics of structures and buildings**URL:**http://journals.rudn.ru/structural-mechanics/article/view/21083**DOI:**http://dx.doi.org/10.22363/1815-5235-2019-15-2-158-168

#### Abstract

Aims of research. Expert validation of all proposed design solutions, development of necessary design solutions for the heightening of the Limon dam according to the ICOLD recommendations. Methods. The detailed static and seismic (dynamic) analyses of stressstrain state and seepage of concrete face rockfill dam Limon (Peru) were performed using the advanced software FLAC-3D (USA) and PLAXIS 2D (Holland), respectively. The elasto-plastic model with Mohr - Coulomb criterion with variable shear angles of gravel and pebble zones of dam materials and its foundation soils was used in the static and seismic (dynamic) analyses of the dam. The dynamic nonlinear analyses of stress-strain state of two variants of Limon dam with full reservoir under Maximum Credible Earthquake (MCE) action of the Mar-Chile Earthquake accelerogram. Results. On the base of these analyses the recommendations were developed for the project of the dam heightening from 43 up to 82 m before the initial filling of the reservoir. Expert validation of all proposed design solutions, necessary design solutions for the heightening of the Limon dam were developed according to the ICOLD recommendations.

Introduction1 In July 2012 the Government of Lambayeque province (Peru) invited the author of this article as an international expert and member of ICOLD to perform the expert validation of design of the heightening of concrete face gravel dam (CFGD) Limon from 43 to 82 m. The dam is the main element of project “Proyecto Especial Olmos - Tinajones (PEOT)”. The hydraulic transfer scheme of 1© Lyapichev Yu.P., 2019 This work is licensed under a Creative Commons Attribution 4.0 International License the project includes the TransAndes water-transfer 26 km long tunnel now completed. The 82 m high Limon CFGD is located on the right bank of Huancabamba river in remote region of Andes wit very high seismicity. Maximum Credible Earthquake (MCE) with the return period T = 5000 years and Amax = 0.57g corresponds to ICOLD recommendations: Bulletins 148, 122, 154, 155, 167 [1-5] and was much more dangerous than adopted in 2009 Brazilian design: Amax = 0.39g, T ≈ 1000 years [6]. In the first PEOT project, developed by Hydroproject Institute (Moscow) in 1982, the variant of 82 m high Limon rockfill dam with clay core 158 DYNAMICS OF STRUCTURES AND BUILDINGS Ляпичев Ю.П. Строительная механика инженерных конструкций и сооружений. 2019. Т. 15. № 2. С. 158-168 was adopted for the one-stage dam construction. But later due to the political and financial problems in Peru the project implementation was delayed for 20 years and was resumed as a twostage construction by BOT scheme (build, operate and transfer to owner), proposed by Odebrecht construction company (Brazil). The company changed the Soviet design of one-stage 82 m high Limon traditional rockfill dam with clay core in favor of the two-stage CFGD (43 and 82 m high). The Soviet project of shore spillway remained unchanged. 1. Seepage analysis of Limon CFGD (H = 82 m) and its alluvial (40 m deep) foundation Seepage analysis was made using the software PLAXIS 2D PlaxFlow (The Netherlands) [7]. In figure 1 is presented the geometry of dam with zones of materials and its 40 m deep alluvial strata of foundation in channel section 10-10' and their permeability coefficients. In figure 2 is presented the finite element mesh of dam and its foundation in channel section 10. Figure 1. Zoning and permeability of soils of CFGD Limon H = 82 m and its foundation in channel section Figure 2. Finite element mesh of CFGD Limon H = 82 m and (40 m deep) foundation in channel section ДИНАМИКА КОНСТРУКЦИЙ И СООРУЖЕНИЙ 159 Lyapichev Yu.P. Structural Mechanics of Engineering Constructions and Buildings, 2019, 15(2), 158-168 Figure 3. Equipotential lines of the total seepage heads in rock foundation below the concrete diaphragm The results of the equipotential lines of total seepage heads in the rock foundation below the plastic concrete diaphragm are showed in figure 4, verifying the significant reduction of the total seepage and pressure heads in the rock foundation by effect of the plastic concrete diaphragm in foundation. Table 1 shows the unit seepage flows in the dam foundation for construction stages of H = 43 m and H = 82 m in sections 8-8' (in right abutment) and 10-10' (channel section) below the concrete diaphragm, in the central dam axis and below the dam toe. The relationship of unit seepage flows in section 8-8' and 10-10' shows that seepage flow in the foundation of the dam H = 82 m would be more than twice the seepage flow in the foundation of the dam H = 43 m. Figure 4. Equipotential lines of seepage pressure heads in rock foundation below the concrete diaphragm Unit seepage flows in the dam foundation for construction stages of H = 43 m and H = 82 m Table 1 Dam Section Unit seepage flows (m3/s/m) Below concrete diaphragm In axis of dam cross-section Below toe of downstream slope I Stage H = 43 m 8-8' 1.373×10-3 1.37×10-3 0.744×10-3 10-10' 1.495×10-3 1.38×10-3 0.652×10-3 II Stage H = 82 m 8-8' 2.865×10-3 2.759×10-3 1.904×10-5 10-10' 3.163×10-3 2.791×10-3 0.536×10-3 Relation QH82/QH43 8-8' 2.09 2.01 2.56 10-10' 2.12 2.02 0.82 160 DYNAMICS OF STRUCTURES AND BUILDINGS Ляпичев Ю.П. Строительная механика инженерных конструкций и сооружений. 2019. Т. 15. № 2. С. 158-168 2. Seismic (dynamic) analysis of 82 m high Limon CFGD under MCE action (Amax = 0.57 g) The main results of dynamic nonlinear analysis of stress-strain state of Limon CFGD (H = 82 m, adopted variant 2 with additional downstream rockfill zone) with full reservoir under Maximum Credible Earthquake (MCE) action of the Mar-Chile Earthquake accelerogram are given in figures 5-17. The previous dynamic analysis of variant 1 (without this downstream rockfill zone) are omitted. Another MCE of the Lima-Peru Earthquake accelerogram was considered also in the dynamic analysis, but its action was less dangerous than that of the Mar-Chile Earthquake accelerogram. In figure 5 the accelerogram of Mar-Chile Earthquake normalized to the maximum acceleration of Amax = 0.57g is shown. The Mar-Chile Earthquake with the return period T = 5000 years and Amax = 0.57g corresponds to the recommendations of the ICOLD Bulletins [1-5] and was much more dangerous than adopted in previous (2009) Brazilian design: Amax = 0.39 g, T ≈ 1000 years [6]. The static and dynamic analyses of stress-strain state of Limon CFGD (H = 43 and 82 m) were made by FLAC-3D software (USA) [8], which was estimated in ICOLD Congress (Canada, 2003) [9] as one of the best software for dynamic analyses of large rockfill dams including CFRDs. The finite element model of Limon CFGD (H = 43 and 82 m) with its foundation is shown in the figure 6. Figure 5. Accelerogram of Mar-Chile Figure 6. The finite element model of Limon CFGD Earthquake normalized to Amax = 0.57 g (H = 43 and 82 m) with its foundation Parameters of the elasto-plastic model with Mohr - Coulomb criterion for dam materials and foundation soils in static analyses of Limon CFGD (H = 43 and 82 m) are given in table 2. Parameters of Mohr-Coulomb model in static analyses of Limon CFGD (H = 43 and 82 m) Table 2 Numbers and names of zones of dam materials and foundation soils Material or soils Dry density and void ratio Parameters of deformation Parameters of shear strength of materials γdr, t/m3 n E (MPa) Angle of dilatancy (°) ν C (MPa) ψ(°) 1st stage dam (H = 43 m) 1, 3. Foundation Alluvium 2.15 0.2 108 0 0.30 0 42 2. Diaphragm Concrete 2.25 0 320 0 0.40 0.4 30 4. Plint slab Concrete 2.5 0 20000 0 0.17 1.0 60 5. Embankment zone Gravels and pebbles 2.2 0.15 168 0 0.30 0 46.5 6. Transition zone Gravels 2.15 0.2 150 10 0.33 0 42 7. Transition zone Sand 2.1 0.25 100 10 0.33 0 40 8. Concrete face Concrete 2.5 0 20000 0 0.17 1.0 60 2nd stage dam (H = 82 m) 9. Embankment zone Gravels 2.2 0.15 168 0 0.30 0 46.5 10. Embankment zone Gravels 2.2 0.15 168 0 0.30 0 46.5 11. Transition zone Gravels 2.15 0.2 150 10 0.33 0 42 12. Transition zone Sand 2.1 0.25 100 10 0.33 0 40 13. Concrete face Concrete 2.5 0 20000 0 0.17 1.0 60 14. Downstream zone with 2 berms Pebbles 2.1 0.25 150 0 0.30 0 46.5 ДИНАМИКА КОНСТРУКЦИЙ И СООРУЖЕНИЙ 161 Lyapichev Yu.P. Structural Mechanics of Engineering Constructions and Buildings, 2019, 15(2), 158-168 between the dam crest and upper alluvial layers of dam foundation. Figure 7. Scheme of CFGD Limon (H = 42 and 82 m) (adopted variant with d-s zone 14 with 2 berms) Figure 9. Distribution of seismic accelerations through the dam height (H = 82 m) using shear wedge method Figure 8. Scheme of CFGD Limon (H = 43 and 82 m) with variable shear angles of gravel and pebble zones 10-11, 15-17 Table 3 Values of shear angles of gravel and pebble zones 10-11, 15-17 depending on normal stresses Normal stresses, σn , MPa 0.2 0.5 0.8 1.0 ≥1.2 Shear angles y(°) of gravel and pebble zones 46.5 46.3 42.0 41.1 40.0 Figure 10. Factors of seismic (Fmin = 1,19 > Fperm = 1,06) and static stability (Fmin = 1,69 > Fperm = 1,25) of downstream slope Scheme of zoning of CFGD Limon (H = 82 m) with variable shear angles of gravel and pebble zones 10-11, 15-17 (figure 8) was used in the pseudostatic analyses of the downstream slope stability under action of the acceleration in dam foundation Ahor = 2/3 • Amax = 2/3 • 0.57 g = 0.38 g. The distribution of seismic accelerations through the dam height was received according to Russian seismic design norms for dams (SNiP 33-01-2003) using the shear wedge method (figure 9). Figure 10 shows results of static (the most dangerous circular surface 2) and seismic (the most dangerous circular surface 1) stability of downstream slope of Limon CFGD (H = 82 m) taking into account the variable shear angles of gravel and pebble zones 10-11, 15-17. This figure show that the minimum factor of the downstream slope stability under action of seismic loads is more that permissible as per design norms SNiP 33-01-2003 (Fmin = 1,22 > Fperm = 1,06) and corresponds to the deep circular sliding surface The comparison of results of the seismic stability analysis of the downstream slope of Limon CFGD (H = 82 m) with additional pebble zone 14 with two berms (figure 7) with results of the same analysis of the dam but without the additional zone show that the inclusion of this zone in the downstream slope provide a significant increase of the minimum factor of the downstream slope stability from 1.05 up to 1.22. Below in figures 11, 13-15 the main results of dynamic nonlinear analysis of stress-strain state of Limon CFGD (H = 43 and 82 m, variant with the additional downstream pebble zone) with full reservoir under action of MCE of the Mar-Chile Earthquake accelerogram are presented. The dam zones with the shear stress state of soils are painted in orange and zones with the tension stress state of soils are painted in blue (figure 11). Parameters of Mohr-Coulomb model used in the dynamic nonlinear analysis of Limon CFGD (H = 43 and 82 m) are given in table 4. 162 DYNAMICS OF STRUCTURES AND BUILDINGS Ляпичев Ю.П. Строительная механика инженерных конструкций и сооружений. 2019. Т. 15. № 2. С. 158-168 Figure 11. Zones of Limon CFGD (H = 43 and 82 m) with the shear and tension stress state of soils under action of SMC of the Mar-Chile Earthquake accelerogram with Amax = 0.57 g Parameters of Mohr-Coulomb model in dynamic analyses of Limon CFGD (H = 43 and 82 m) Table 4 Numbers and names of zones of dam materials and foundation soils Material or soils Dry density and void ratio Dynamic modulus of elasticity Edyn, MPa Shear modulus Gmax (MPa) Initial coefficient of damping ξ, % Reduction of parameters Gmax and ξ γdr, t/m3 n 1st stage dam (H = 43 m) 1, 3. Foundation Alluvium 2.15 0.2 1300 Gmax = 35(σm)0.5 5 see figure A 2. Diaphragm Concrete 2.25 0 1600 G = Edyn / [2(1 + ν)] 3 - 4. Plint slab Concrete 2.5 0 20000 G = Edyn / [2(1 + ν)] 2 - 5. Embankment zone Gravels and pebbles 2.2 0.15 2000 Gmax = 40(σm)0.5 5 see figure A 6. Transition zone Gravels 2.15 0.2 1000 Gmax = 22(σm)0.5 4 see figure A 7. Transition zone Sand 2.15 0.2 700 Gmax = 20(σm)0.5 4 see figure A 8. Concrete face Concrete 2. 5 0 20000 G = Edyn / [2(1 + ν)] 5 - 2nd stage dam (H = 82 m) 9. Embankment zone Gravels 2.2 0.15 2000 Gmax = 40(σm)0.5 5 see figure A 10. Embankment zone Gravels 2.2 0.15 2000 Gmax = 40(σm)0.5 5 see figure A 11. Transition zone Gravels 2.15 0.2 1000 Gmax = 22(σm)0.5 4 see figure A 12. Transition zone Sand 2.1 0.25 700 Gmax = 20(σm)0.5 4 see figure A 13. Concrete face Concrete 2.5 0 20000 G = Edyn / [2(1 + ν)] 5 - 14. Downstream zone with 2 berms Pebbles 2.15 0.2 1500 G = Edyn / [2(1 + ν)] 5 see figure A N o t e : (σm) - medium stress (effective) in kPa. Figure 12. Curves of reduction of the shear modulus G/Gmax and initial coefficient of damping ξ, % of soils of the dam and its foundation Under action of the Mar-Chile Earthquake the dam would suffer elasto-plastic deformations with large plastic displacements in the wide zone of the downstream slope (figure 13). The large plastic deformations modified the dynamic stress-strain state of the dam and its foundation (figures 13-14). The horizontal and vertical displacement in the dam after the Mar-Chile Earthquake in the upper part of the downstream slope are, respectively, 2.0 and 1.0 m; in the upper berm - 2.2 and 1.1 m; in the lower berm - 2.5 and 1.3 m and at the toe of the slope - 6.0 m and zero (figure 13). The intensity of shear deformations (figure 13) is concentrated in the narrow zone in lower part of dam downstream slope. ДИНАМИКА КОНСТРУКЦИЙ И СООРУЖЕНИЙ 163 Lyapichev Yu.P. Structural Mechanics of Engineering Constructions and Buildings, 2019, 15(2), 158-168 Figure 13. Horizontal (a) and vertical (b) displacements in Limon dam (82 m) after Mar-Chile Earthquake Figure 14. Intensity of the shear deformations in Limon dam (82 m) after the Mar-Chile Earthquake Figure 15. The time history of the residual horizontal (a) and vertical (b) displacements of the crest of Limon dam (82 m) during the Mar-Chile Earthquake The time history of the residual horizontal (a) and vertical (b) displacements of the dam crest during the Mar-Chile Earthquake is shown in figure 15. The maximum horizontal and vertical displacements of the dam crest during the Mar-Chile Earthquake are, respectively, 1.5 and 1.1 m. 164 DYNAMICS OF STRUCTURES AND BUILDINGS Ляпичев Ю.П. Строительная механика инженерных конструкций и сооружений. 2019. Т. 15. № 2. С. 158-168 3. Results of analysis of the concrete face These results are showed in figure16 and 17 (Limon dam, the 2nd stage, H = 82 m). The displacements, bending moments and longitudinal forces (axial) in the concrete face are presented for action of Maximum Credible Eartquake (MCE) of Mar-Chile after filling of the reservoir up to maximum elevation (Limon dam, the 2nd stage, H = 82 m). It’s shown that the greatest influences on the concrete face is axial compression, being of lesser bending moment value, therefore, the concrete face will be in compression state. With the dimensions of concrete face, adopted forces and moments it can determine the bearing capacity and reinforcement of thick 0.55 and 0.42 m concrete face based on the diagram of interaction force-bending. Also this reinforcement is also recommended to absorb stresses due to shrinkage and thermal changes in concrete face. Figure 16. Horizontal (a) & vertical (b) displacements of concrete face after filling of the reservoir up to maximum elevation (CFRD of the 2nd stage, H = 82 m): Figure 17. Longitudinal forces (a) and bending moments (b) in the concrete face after filling of the reservoir up to maximum elevation (CFRD of the 2nd stage, H = 82 m) LEGEND FOR FIGURES 16 AND 17: A - point on perimetral joint; B - intermedial point; C - point on the crest of dam of the 1st stage; A - point on perimetral joint; B - intermedial point; C - point on the crest of dam of the 2nd stage. Displacements in points of concrete face: Forces N & bending moments M in points of concrete face: A - Ux = 0,14 m; Uz = -0,16 m A - N = -80 t (compressión); M = -6 tm B - Ux = 0,16 m; Uz = -0,18 m B - N = -40 t (compressión); M = -2 tm C - Ux = 0,14 m; Uz = -0,14 m C - N = -50 t (compressión); M = 0 tm D - Ux = 0,08 m; Uz = -0,10 m D - N = -50 t (compressión); M = -2 tm E - Ux = 0,06 m; Uz = -0,04 m E - N = -50 t (compressión); M = -4 tm Deflections of concrete face: in point B (maximum) - 24 cm; in point C - 19 cm; in point D - 13 cm ДИНАМИКА КОНСТРУКЦИЙ И СООРУЖЕНИЙ 165 Lyapichev Yu.P. Structural Mechanics of Engineering Constructions and Buildings, 2019, 15(2), 158-168 4. Supposed designed behavior of the concrete face of Limon dam From the pattern of deformations the behavior of concrete face acts as a rigid wall to vibrations in the transverse direction of the valley, differs considerably from the rockfill and transition zone materials, the transverse response across the valley of rockfill dam may be limited by the relatively rigid concrete face. This may lead to tension stresses in the plane of concrete face, seismic forces transferred from rockfill to concrete face are limited by the frictional forces between the transition zone of rockfill and concrete face. Like all the water load is supported by the concrete face, these frictional forces are relatively high and therefore stresses in the plane of concrete face can be significant enough to cause local deformations or shears of concrete face slabs along its longitudinal joints. The dam deformations can cause the crack opening in concrete slabs and sliding along the crack surface and suffer oscillating movements. This behavior of Limon dam corresponds to many case studies of behavior of concrete face rockfill dams, in detail described and discussed in [10-15]. The main advantage of this dam with concrete face is its high resistance to erosion in case of water seepage through cracked concrete face. If the zone of materials under the concrete face has a proper particle sizes the permeability coefficient will be 10-3 cm/s and then zone will be stable against its erosion. This eliminates the risk of high leakages developed under cracked concrete face. Main conclusions 1. Horizontal and vertical displacements of the downstream slope after the Mar-Chile Earthquake are, respectively, 2.0-2.5 and 1-1.3 m. The maximum horizontal and vertical displacements of dam 3. During the forthcoming reservoir filling the strict monitoring of instruments embedded in dam body (piezometers, benchmarks, strain gauges, accelerometers) with expert control is necessary.

### Yury P Lyapichev

JSC “Institute Hydroproject”, International Commission on Large Dams (ICOLD)
**Author for correspondence.**

Email: lyapichev@mail.ru

2 Volokolamskoe Shosse, Moscow, 125993, Russian Federation, 61 Avenue Kleber, Paris, 75016, France

Doctor of Technical Sciences, Professor, Expert of JSC “Hydroproject Institute” for foreign projects, Member of International Commission on Large Dams (ICOLD).

- ICOLD Bulletin 148. (2010). Selecting seismic parameters for large dams. Guidelines.
- ICOLD Bulletin 122. (2001). Computational procedures for dam engineering.
- ICOLD Bulletin 154. (2011). Dam safety management. Operational phase of the dam life cycle.
- ICOLD Bulletin 155. (2013). Guidelines for use of numerical models in dam engineering.
- ICOLD Bulletin 167. (2015). Regulation of dam safety. Overview of current practice.
- Informe de gestión Proyecto Especial Olmos - Tinajones (PEOT). (Junio 2009).
- PLAXIS 2D PlaxFlow. The Netherlands. https:// www.plaxis.com/product/plaxis-2d/
- FLAC-3D software. USA. http://www.itascacg.com/
- Proceedings of the ICOLD Congress. (2003). Canada.
- SPANCOLD. (1996). Proceedings of the 4th Benchmark Workshop on Numerical Analysis of Dams. Madrid, Spain.
- USCOLD. (1999). Proceedings of the 5th Benchmark Workshop on Numerical Analysis of Dams. Denver, Colorado, USA.
- ATCOLD and VERBUND. (2001). Proceedings of the 6th Benchmark Workshop on Numerical Analysis of Dams. Salzburg, Austria.
- ROCOLD. (2003). Proceedings of the 7th Benchmark Workshop on Numerical Analysis of Dams. Bucharest, Romania.
- Wuhan University. (2005). Proceedings of the 8th Benchmark Workshop on Numerical Analysis of Dams. Wuhan, China.
- Proceedings of 9th Benchmark Workshop on Numerical Analysis of Dams / VNIIG (Russian Research Institute of Hydraulic Engineering named after B.E. Vedeneev). Russia: Saint Petersburg, 2007
- COYNE ET BELLIER. (2009). Proceedings of the 10th Benchmark Workshop on Numerical Analysis of Dams. Gennevilliers, France.
- Lyapichev Yu.P., Landau Yu.A. (November - December 2011). Modern structural and technological solutions for new large dams. Hydro Review Worldwide. PennWell Corp., USA.
- Lyapichev Yu.P. (2013). Design, construction and behavior monitoring of modern high dams. Palmarium Academic, Germany, 369. (In Russ.)