In this paper, transient responses of the Brezina concrete arch dam, located in the west of Algeria, are determined using the finite elements commercial package, Ansys. A 3D model has been created to study the rock–fluid–structure interaction phenomenon. The foundation rock is modeled as a mass rock using solid finite elements, and the reservoir water is modeled using fluid finite element for the 3D model; the length and width of the foundation rock and water reservoir, along the global
Keywords
- Concrete dam
- Coupling equations
- Contact elements
- Transient analysis
One of the main topics in earthquake engineering is the dynamic analysis of dam–reservoir systems. Research on the dam–foundation interaction problem is available extensively in the literature [1,2,3,4,5,6,7,8].
There are several extensive methods that have been applied in past decades. Numerical procedures that include the interaction between several domains with different properties—concrete dam, foundation rock, water, bottom sediments, and bank of the reservoir—have been developed using the finite element method, the boundary element method, and various combinations of both the methods [9,10,11,12, 21, 22].
Evaluation of the important hydrodynamic forces that develop on the upstream face of a large dam during severe transient excitations has been the subject of numerous studies, starting with Westergaard's classical work in 1933. Westergaard explained the physical behavior of dam–reservoir interaction for 2D coupled systems [13].
According to the literature, two approaches may be used to model the dam–reservoir interaction phenomenon: the coupling method and the contact elements method [1].
As part of the comprehensive study undertaken on the dynamic response of the Brezina arch dam, located in El Beyadh, Algeria, and following recent works [14, 19], in this paper, the influence of the reservoir water on the transient response of the Brezina arch dam is investigated. This article is the continuity of Berrabah
The discretization of the structural dynamic equation taking into account both the dam and foundation rock subjected to ground motion, using the finite-element approach, can be given as:
The hydrodynamic pressure equation discretization is given as:
The assembled form of both Eqs (1) and (2) can be given as:
Eq. (3) is the second-order linear differential equation having unsymmetrical matrices and it can be solved using direct integration methods.
The object of the present work is the Brezina concrete arch dam, shown in Figure 1. It would be helpful to recall that this article is a continuity of Berrabah
Figure 1
Brezina concrete arch dam [ref. Wikipedia.org]

Geometrical characteristics of Brezina dam
60 | 78.5 | At crest: 5 |
At foundation level: 36.3 |
The dam–reservoir interaction system is investigated using the finite element commercial package, Ansys, with a mapped meshing [14] 3D finite element model.
In order to determine the fineness of meshing [13], a sensitivity study was performed with respect to the modal and static responses.
The 3D model features 16,252 quadratic solid elements (SOLID186) and 3,000 fluid elements (FLUID80).
It is important to note that coupling equations are used to model the fluid – “rock foundation” and/or fluid–“concrete dam” contact; however, for the “concrete dam”–“rock foundation” interface, “perfect contact” available in the Ansys code is used to ensure continuity and linearity.
Table 2 summarizes the material properties of the Brezina concrete arch dam, rock–foundation and reservoir water.
Material properties of Brezina arch dam, foundation rock, and reservoir water
Concrete dam | 28.5e+09 | 0.2 | 2,500 | |
Foundation rock | 14.5e+09 | 0.25 | 2,100 | |
Reservoir water | Compressibility modulus (N/m2) | Viscosity | 0.49 | 1,000 |
2.068E+09 | 0.001 |
As shown in Figure 2, for the 3D model, the length and width of the foundation rock and water reservoir, along the global
Figure 2
3D FE model of “fluid–dam–foundation rock” system and applied BCs. BCs, boundary conditions

It is important to note that both rock foundation and reservoir water extensions are chosen after doing sensitivity analysis of the dam behavior via these two parameters. For this reason, an extension of both rock foundation and water reservoir was taken equal first as 50 m and a transient analysis was performed and results were extracted. Then these extensions were increased until the convergence of the extracted results could be obtained (the above extension values are judged as sufficient). At these extension limits, fixed supports are applied. So, at these limits, no radiation damping is present.
It is also worth mentioning that wave propagation is not taken into account, but this can be used as a perspective.
The “reservoir water”–“concrete dam” and “reservoir water”–“rock foundation” coupling is obtained as follows: at an interface, the displacement of the nodes belonging to the “reservoir water” and that of the nodes belonging to the “concrete dam” or “foundation rock” are coupled in the normal direction of the interface. The other two displacement components of these nodes, i.e., the displacement in the tangential and vertical directions, are left free (not coupled). In the studies by Fok and Chopra [10], Tan and Chopra [19], and Sani and Lotfi [20], which employ the Eulerian approach, the coupling is made between the displacements and pressures resulting in the interaction matrix. In this study the Lagrangian approach is employed; the reference configuration of the continuum, solid, or fluid, is the undeformed state; and the behavior of both the fluid and structure is expressed in terms of the displacement [19].
Two kinds of static loadings are considered in the present article: the self-weight of the dam (or the dam–rock system) and the hydrostatic pressure in the case of the full “reservoir water” system.
To perform transient analyses for the “reservoir water”–“concrete dam”–“foundation rock” model, three synthetic earthquake records (accelerograms) have been generated [18]. Figure 3 depicts the accelerogram in the
Figure 3
Acceleration time history along

As mentioned in Berrabah
For the Newmark formulation, both the final velocity basic integration equations and the final displacement basic integration equations are expressed as follows:
A time step that is too large will introduce error that affects the response of the higher modes (and hence the overall response). A time step that is too small will waste computer resources. To calculate an optimum time step, the following guidelines extracted from the ANSYS User's Manual [14] are used:
“Resolve the response frequency. The time step should be small enough to resolve the motion (response) of the structure. Since the dynamic response of a structure can be thought of as a combination of modes, the time step should be able to resolve the highest mode that contributes to the response. For the Newmark time integration scheme, it has been found that using approximately twenty points per cycle of the highest frequency of interest results in a reasonably accurate solution”. That is, if “Resolve the applied load-versus-time curve(s). The time step should be small enough to ‘follow’ the loading function. The response tends to lag the applied loads, especially for stepped loads. Stepped loads require a small ITS at the time of the step change so that the step change can be closely followed” and given as:
Resolve the wave propagation. If you are interested in wave propagation effects, the time step should be small enough to capture the wave as it travels through the elements”. This guideline is provided as a complement; wave propagation effects are not of interest in the present work. In agreement with the first two guidelines, in the present analyses, the ITS is chosen as shown in Table 3:
ITS for the two dam models
3D full dam with mass foundation | 7.566 | 0.0066 | 0.0007342 |
ITS, integration time step
Table 3 summarizes the fundamental frequencies for the studied dam case (after doing modal analyses) and the ITS to resolve both the response frequency and the applied load-versus-time curve(s).
To respect the above ITS values and to ensure a reasonable running time duration the ITS chosen are as follows:
Initial time step = 0.001 s Minimum time step = 0.0001 s Maximum time step = 0.01 s
By activating nonlinear effects represented by the option of “Large Displacement Transient” and automatic time stepping, Ansys adapts the time step (ITS) to take into account the natural frequencies and mode-shapes of the structure by optimizing the time step. The choice of the minimum and maximum time step allows to cover all the frequencies of interest of the fluid–dam–foundation rock, which vary from 0.1 Hz to 50 Hz [17].
Recall that a damping ratio of 5% is taken (using the Rayleigh method).
Transient analyses are performed for two 3D systems as well – the “concrete dam”–“foundation rock” system (“concrete dam”–“foundation rock” interaction case), and the “concrete dam”–“reservoir water”–“foundation rock” system (“concrete dam”–“reservoir water”–“foundation rock” interaction case) – to study the effects of reservoir water on the dynamic behavior of the system. Two paths are chosen to represent the results for the upstream dam face: the cantilever path and the crest path (Figure 4).
Figure 4
Cantilever and crest path for the upstream dam face

Figure 5 plots the maximum displacement in the
Figure 5
Maximum displacement in

Figure 6 depicts the maximum displacement in the
Figure 6
Maximum displacement in the

Figure 7 plots the maximum displacement in the
Figure 7
Maximum displacement in the

Figure 8 shows the effect of reservoir water on the maximum stress in the
Figure 8
Maximum stress in the

Similarly, the effect of water presence on the maximum stress in the
Figure 9
Maximum stress in the

Figure 10 plots the effect of water presence on the maximum stress in the
Figure 10
Maximum stress in the

Figure 11 shows the effect of water presence on maximum von Mises stress.
Figure 11
Maximum von Mises stress along the concrete dam cantilever path for the two studied cases

Figure 12 plots the maximum displacement in the
Figure 12
Maximum displacement in the

Similarly, Figure 13 depicts the maximum displacement in the
Figure 13
Maximum displacement in the

Figure 14 plots the maximum displacement in the
Figure 14
Maximum displacement in the

Figure 15 shows the effect of reservoir water presence on maximum stress in the
Figure 15
Maximum stress in the

Similarly, the effect of water presence on the maximum stress in the
Figure 16
The maximum stress in the

Figure 17 plots the effect of reservoir water on the maximum stress in the
Figure 17
The maximum stress in the

Figure 18 shows the effect of reservoir water presence on the maximum von Mises stress along the crest path. As in Figure 17, the influence is small except at the end part of the path, where an increase due to reservoir water can be seen. The highest value, at 1 MPa, is very small.
Figure 18
Maximum von Mises stress along the dam upstream face path for the two studied cases

This study presents the transient responses of the “Brezina concrete dam”–“reservoir water”–“mass rock–foundation” system using 3D Ansys finite element model. The dam is subjected to an acceleration seismic record applied in the
Along the dam cantilever part:
Reservoir water increases the longitudinal displacement, the effects are negligible on the other displacement component(s), and the displacement magnitude is very small. With or without reservoir water, the longitudinal and vertical stresses are in compression. For the 3D model, reservoir water diminishes (the influence is small) von Mises stress.
Along the dam crest part:
Reservoir water increases the longitudinal displacement, while it drastically diminishes the vertical displacement. The influence on the displacement along the crest direction is small. The displacement magnitude is also very small. With or without reservoir water, the stress components are mostly in compression. Reservoir water is seen mostly to be diminishing the stresses, although, on the whole, its influence on von Mises stress is small. Again, the stress magnitude is very low.
Globally speaking, the role of reservoir water is mixed. In some cases, its presence is like a damper, which decreases the range of displacement or stresses, whereas in other cases, it increases or has a negligible effect on these quantities. Finally, it is worth emphasizing that the stresses determined are very low; thus, if the dam were to be subjected to a combination of three accelerations, it is expected that the resulting stress magnitude obtained would remain at a low or very low level.
Figure 1
![Brezina concrete arch dam [ref. Wikipedia.org]](https://sciendo-parsed-data-feed.s3.eu-central-1.amazonaws.com/6113c8f8927d846dfda25fe1/j_mme-2021-0015_fig_001.jpg?X-Amz-Algorithm=AWS4-HMAC-SHA256&X-Amz-Date=20220524T065258Z&X-Amz-SignedHeaders=host&X-Amz-Expires=18000&X-Amz-Credential=AKIA6AP2G7AKDOZOEZ7H%2F20220524%2Feu-central-1%2Fs3%2Faws4_request&X-Amz-Signature=f049580a25bcaac54099538236a005d5ceed9a5e1c8dabb477e8d8e73ff8b596)
Figure 2

Figure 3

Figure 4

Figure 5

Figure 6

Figure 7

Figure 8

Figure 9

Figure 10

Figure 11

Figure 12

Figure 13

Figure 14

Figure 15

Figure 16

Figure 17

Figure 18

Material properties of Brezina arch dam, foundation rock, and reservoir water
Concrete dam | 28.5e+09 | 0.2 | 2,500 | |
Foundation rock | 14.5e+09 | 0.25 | 2,100 | |
Reservoir water | Compressibility modulus (N/m2) | Viscosity | 0.49 | 1,000 |
2.068E+09 | 0.001 |
Geometrical characteristics of Brezina dam
60 | 78.5 | At crest: 5 |
At foundation level: 36.3 |
ITS for the two dam models
3D full dam with mass foundation | 7.566 | 0.0066 | 0.0007342 |