Nonlinear Seismic Analysis of Concrete Gravity Dams Using Extended Finite Element method and Plastic Damage Model

Document Type : Original Research

Authors
M. Haghani 1
B. Navayi Neya 1
M. T. Ahmadi 2
J. Vaseghi Amiri 1
1 Babol Noshirvani University of Technology
2 Tarbiat Modares University
Abstract
The present paper investigates and compares the crack propagation in concrete gravity dams using two models of linear fracture mechanics and plasticity damage concrete. The first model is based on linear concrete behavior using the extended finite element method without considering the effect of strain softening on the crack tip while the second model is based on the nonlinear concrete behavior and the strain softening in tension with damage parameter. According to two different algorithms and based on two models, several benchmark examples are reviewed and the results compared with those reported in the literature. Then, path of the crack growth in Koyna gravity dam due to a seismic excitation of Koyna earthquake in 1962 has been performed by considering the dam-reservoir interaction.

The results show that due to low compressive stresses during analysis of concrete gravity dams, consideration of compressive nonlinear behavior has no effect on crack initiation and almost is the same for two models. However because of crack opening and closing with tapping the crack faces together in extended finite element model, the compressive stress will be more than the allowable stress of concrete. Crack initiation at downstream and upstream face occurred at angle of 90 and zero degrees respectively, which in both models, the numerical results are in agreement with the experimental model.

The crack in the extended finite element model grows faster such that the crest block of dam in this model is separated from the dam body, earlier than the concrete plastic damage model. Also the values of dam crest displacement and hydrodynamic pressure in the reservoir in extended finite element model with linear elastic fracture mechanic are more than the other model, which can be attributed to the linear and nonlinear behavior of concrete in extended finite element and concrete plastic damage model respectively. In the extended finite element model, due to using linear fracture mechanic, the maximum principal stress in the cracked elements reaches the values greater than the maximum tensile strength, but in the concrete plastic damage model as soon as the stress reaches a tension limit value, elements are damaged and the stress is reduced. In both models, the crack located at the slope change area, propagates with the downward slope from downstream dam face and connects to the crack at upstream face which is growth horizontally.

Because of laboratory sample dimension and boundary condition of dam-reservoir compared with actual manner, neither of two crack profiles covered the experimental model, accurately. But it is shown that the crack profiles in the extended finite element model are more consistent with experimental results. Finally, the results show that the crack profile are slightly different in the two models because of quasi brittle behavior of the dam concrete, which can be attributed to the small fracture process zone of the crack tip in comparison with the dimension of the concrete gravity dams such that by removing strain softness part, the error in the amount of additional computation can be neglected.

Keywords

Subjects

[1] Bhattacharjee S. S., & Leger, P. 1993: Finite element modelling of the tensile strain softening behaviour of plain concrete structures, Engineering Computations, 10(3), 205-221.
[2] Ayari ML, Saouma VE. A. 1990: fracture mechanics based seismic analysis of concrete gravity dams using discrete cracks. Eng Fract Mech, 35(1–3), 587–598.
[3] Ahmadi MT, Izadinia M, Bachmann H. 2001: A discrete crack joint model for nonlinear dynamic analysis of concrete arch dam. Comput Struct, 79(4), 403–420.
[4] Lee J. and Fenves, G.L. 1998-a: A plastic damage model for cyclic loading of concrete Structures, Journal of Engineering Mechanics ASCE, 124, 892–900.
[5] Belytschko T, Black T. 1999: Elastic crack growth in finite elements with minimal remeshing.
Int J Numer Methods Eng, 45(5), 601–620.
[6] Melenk JM, Babuška I. 1996: The partition of unity finite element method: basic theory and applications. Comput Methods Appl Mech Eng, 139(1–4), 289–314.
[7] Moes N, Dolbow J, Belytschko T. 1999: A finite element method for crack growth without remeshing. Int J Num. Meth. Eng., 46(1), 131–150.
[8] Lubarda, V.A., Kracjinvovic, D. and Mastilovic, S. 1994: Damage model for brittle elasticsolids with unequal tensile and compressive strength, Engineering Fracture Mechanics, 49, 681–697.
[9] Khoie, A. 2015: Extended Finite Element Method Theory and Applications, Wiley Series in computational mechanics.
[10] Lee, J. and Fenves, G.L. 1998-b: A plastic–damage model for earthquake analysis of dams. Earthquake Engineering and Structural Dynamics, 1(27), 937-956.
[11] Cevera, M., J. Oliver, J. and Manzoli, E. 1996: A rate-dependent isotropic damage model for the seismic analysis if coccrete dams, Earthquake engineering and structural dynamics, 25,987-1010.
[12] Calayir A. & Karaton, M. 2005: A continuum damage concrete model for earthquake analysis of concrete gravity dam–reservoir systems, Soil Dynamics and Earth. Eng., 25(11), 857-869.
[13] Zhang S., Wanga, and G., Yu X. 2013: Seismic cracking analysis of concrete gravity dams with initial cracks using the extended finite element method Engineering Structures, 56, 528–543.
[14] Omidi O. and V. Lotfi V. 2012: Plastic–Damage Analysis of Koyna Dam in Different Damping Mechanisms, Water In. 15 WCEE.
[15] Vaseghi J. 1998: Finite Element Modelling Of Cracking In Concrete Gravity Dams, Ph.D. thesis, University Of Tarbiat Modares, Iran (In Persian).
[16] Cai Q., Finite Element Modelling Of Cracking In Concrete Gravity Dams, Ph.D. thesis, University Of Pretoria, South Africa, 2007.
[17] Broek D. 1986: Elementary Engineering Fracture Mechanics, Hingham, MA, Kluwer Academic Publishers.
[18] Bazant Z. P. Pfeiffer P. A. 1987: Determination of fracture energy from size effect and brittleness number, ACI Materials Journal 84, 463–480.
[19] Bui HD. Mécanique de la rupture Fragile. Masson: New York, 1978.
[20] Kalani Sarokolayi L., Navayi Neya B., Vaseghi Amiri J. 2015: Nonlinear dynamic analysis of concrete gravity dams considering rotational component, Int J of civil Eng, 13(1), 16-29.
[21] Abaqus C.A.E User Manual. Ver. 2018.
[22] National Research Council (US). 1990: Earthquake engineering for concrete dams design, performance, and research needs Washington (DC): National Academies Press, 99–105.
[23] Hall, J.F. 1998: The dynamic and earthquake behavior of concrete dams: Review of experimental behavior and observational evidence. Soil Dynamics and Earthquake Engineering, 1(7), 57-121.
Modares Civil Engineering journal
Volume 20, Issue 3
September and October 2020
Pages 189-201