EvaluationofHydrodynamicPressureDistribution inReservoir of Concrete Gravity Dam under Vertical Vibration Using an Analytical Solution

Fluid-structure interaction causes a hydrodynamic force, which can be exerted to the dam and affects its response. )e effect of vertical excitation of groundmotion on dynamic behavior of concrete gravity dam is themost important because of the interaction between foundation and reservoir. So, the foundation-reservoir interaction should be taken into account in designing concrete dams. In most studies, the effects of the vertical component of vibration have been ignored. While in vertical vibration, due to the interaction of the reservoir and the foundation, a significant hydrodynamic pressure is produced in the tank, which increases the dam response. In this study, the hydrodynamic pressure wave propagation in the reservoir of a concrete gravity dam caused by interaction with the foundation under vertical vibration is investigated using an analytical method. To achieve an analytical solution, the reservoir is assumed to be rectangular, and a harmonic load is vertically applied on the system from the foundation. Considering the acoustic nature of the reservoir fluid under harmonic vibration, a new method using the separation of variables method has been used for solution of hydrodynamic wave equation.)e results show a significant effect of the vertical component of earthquake on the amount of induced pressure distributed in the reservoir, which has been omitted in most previous studies. Obtained results of the proposed model can be extended to more complicated models in terms of different loading and geometrical conditions.


Introduction
Fluid-structure interaction is one of the most extensive discussions concerning the behavior of fluid, and numerous research works are performed on this topic. Fluid-structure interaction occurs when there is a relative motion between the two systems. e dynamic response of systems under the effect of the fluid-structure interaction is of great importance.
During the ground motion, the interaction of the reservoir with dam and foundation can produce a hydrodynamic pressure in the reservoir. Hence, besides the inertia force due to the ground movement, the hydrodynamic force on the upstream face of the dam is one of the main factors exerting a load on the dam structure in the event of a vibration. Neglecting the effects of hydrodynamic force will result inaccurate responses, which can bring some challenges for the designed structure. In other words, taking into account the fluid-structure interaction in the dynamic analysis of concrete dams will increase the accuracy of the responses obtained from the analysis and better estimation of hydrodynamic forces and the stresses created in the dam body.
e first model concerning the investigation of the interaction effect in the seismic behavior of dams was proposed by Westergaard in 1933 [1]. In this model, the damreservoir interaction under the effect of horizontal vibrations was proposed considering the two-dimensional model. In the Westergaard model, the dam has been assumed to be rigid with a semilimited reservoir and a constant depth, and the fluid pressure distribution in the dam-reservoir interface was extracted by the analytical method. Westergaard showed that the induced force because of interaction is proportional to the acceleration caused by the horizontal excitation and can be estimated by a parabolic distribution on the dam height. His method is called the added mass approximation. An important point in Westergaard solution is that the maximum energy caused by the seismic motion is related to the frequencies that are smaller than the natural frequency of the reservoir. Several research studies show that the hydrodynamic pressure distribution caused by an earthquake on the dam body is almost parabolic.
In the Westergaard added mass model, the effect of compressibility of water has been neglected. However, the accurate fluid structure analyses have shown that the compressibility of water has a considerable effect on the hydrodynamic pressure and dynamic behavior of concrete dam. In the dam-reservoir system, the effects of water compressibility are very important when the fundamental frequency of dam is close to the natural excitation frequency of the reservoir. e reason for less attention to the compressibility effects of water is that the investigation of fluidstructure interaction is only necessary for problems related to the wave propagation, and the propagation of the hydrodynamic waves only happens in the compressible flow. If the compressibility effects are neglected, an incompressible model can be used in which the propagation of hydrodynamic waves is neglected.
Chopra investigated the effects of fluid compressibility on the fluid structure interaction and induced force [2]. His model was similar to the Westergaard model. Chopra concluded that in a case of neglecting of water compressibility, a significant error is appeared in the analysis results. In most of the studies, the incompressibility assumption has been used for water because in this case, the modeling of the dam-reservoir system will be simplified. erefore, in this condition, there is no need to solve complicated integrals.
In the following, Chopra studied the effect of dam flexibility on the interaction forces by modeling the dam as a system which composes of mass, damper, and spring [3]. He extracted the parameters of the dam model by focusing on the first special mode of a triangle-shaped dam. He showed that natural frequencies of the dam-reservoir coupled system are different with than that of the uncoupled system. Chopra found that if the ratio of fundamental frequency of dam to natural frequency of the reservoir is higher than two, the compressibility effect of the water can be neglected. In the Chopra model, dam was only under the effect of horizontal acceleration of the ground motion. In 1970, Chopra analyzed the dam-reservoir model to calculate the dam response considering the interaction effect [4]. He showed that by considering the effects of the dam-reservoir interaction, the response of dam significantly increases in comparison with a case in which the effect of interaction is neglected.
ere are few works that addressed the investigation of dynamic pressure distribution in a reservoir under the effect of the vertical component. As an example, in 1971, Nath proposed a model in which the dam was considered as a cantilever beam which has been exposed under the shear and bending effects [5]. By applying the vertical component of the earthquake, the effect of vertical vibration of ground motion on the dam-reservoir interaction was investigated as a 1D model using the finite element method. Nath showed that although the hydrodynamic pressure produced on the dam due to the vertical vibration in a specific frequency is considerably higher than the pressure caused by the horizontal movement under the same frequency, but the displacement response of the structure is higher for the horizontal movement.
Chakrabarti and Chopra approximated the response of the dam under vertical earthquake component by the deformations in the fundamental mode of vibration considering the flexibility of bed rock using the analytical model. Obtained results illustrated that the hydrodynamic forces induced by vertical vibration are affected by dam-reservoir interaction and depend on the flexibility of the dam. ey showed that the vertical component of ground motion induces significant hydrodynamic forces in the horizontal direction on a vertical upstream face. So, it is very important to consider the effect of the vertical component in seismic analysis of concrete gravity dams [6].
In 1988, Humar and Jablonski presented a model in the frequency domain to investigate the hydrodynamic pressure distribution in a reservoir under the effect of vertical component of the earthquake using the boundary element method [7]. ey show that the amount of hydrodynamic pressure produced in the reservoir under the effect of the vertical component of the earthquake is more than that of the horizontal component.
In recent years, some research studies have been performed on the dynamic behavior of the concrete gravity dam using the analytical models. Ghaemian et al. investigated the nonlinear seismic behavior of a concrete gravity dam by considering the heterogeneity and displacement of the foundation [8]. e obtained results from the model show that the heterogeneity of the foundation can affect the system response and causes the rapid occurrence of a crack in the dam body during the earthquake. Ziaolhagh et al. analyzed the dam and foundation system by using the finite element method with triangle elements and 21-point Gaussian integration of the dam system [9]. ey assumed that the dam is flexible, and the system frequencies and mode shape were extracted using the finite element method. ey compared the results of 8-node and 21-node elements and showed that the accuracy of the finite element model with 21-node element is higher. Ribeiro and Pedroso provided a review on the analysis method for concrete gravity dams under the effect of seismic load [10]. ey studied the limits of the previous methods by the finite element method and extracted the dam response using the semistatic, semidynamic and dynamic methods.
Some research has examined the effects of heat transfer on solid contact. Some recent studies have used nanofluids and nanoparticles to more accurately study the partial specifications of water and the state of fluid temperature changes which can be referred to studies by Zahir Shah et al. However, according to the conditions of the dam-reservoir, the thermal changes of the water are not important, the water density is assumed to be constant, and the foundation is impenetrable.
Studying of the previous research studies illustrate that the investigation of the hydrodynamic pressure distribution caused by the foundation-reservoir interaction due to the vertical vibration of the ground motion in the concrete gravity dams has been less evaluated. In other words, the importance of attention to vertical vibration using the numerical and analytical methods is among the topic which has received less attention in previous research.
So, in this study, the hydrodynamic pressure propagation in the reservoir of a concrete gravity dam caused by reservoir and foundation interaction under the vertical excitation is investigated using an analytical model. In the analytical method, some loading and geometrical simplifications may be performed on the model, but because of the response accuracy and the ability to perform the parametric study of variables, it can be a capable practice in studying the interaction effects and hydrodynamic behavior of reservoir. Also, it should be noted that the numerical methods require the coding and high computational efforts, and they are used when obtaining the analytical solution for the problem is impossible.

Dynamic Equation of the Fluid inside the Reservoir
In the problems concerning the fluid-structure acoustic interaction, the dynamic equation of structure should be along with the Navier-Stokes and continuity equations of the fluid. To obtain the governing equation, the fluid inside the reservoir, the following assumptions are credible in the dam-reservoir system [3,6]: (i) e fluid is acoustic and can be assumed to be nonviscous with constant entropy, so that the fluid particles undergo small deformations (ii) Given the specifications of the reservoir, the water density variation is negligible (iii) e velocity values are so small that its transitional effects can be neglected (iv) e water viscosity is small, so the shear stresses can be neglected By applying the assumptions, the continuity equations can be simplified as follows: e relationship between the fluid density and the produced hydrodynamic pressure is defined by the equation of state: where k, ρ, P, V, and t denote the coefficient of fluid volume changes, fluid density, fluid pressure, fluid velocity vector, and time. From equations (1) and (2), it can be written as Neglecting the body gravity effects and based on the assumptions made, the momentum equation of the fluid is summarized as follows: e gradient of momentum equation can be written as According to equations (3) and (5), the following equation is obtained: or where C � ��� k/ρ denotes the speed of acoustic waves within the fluid. Equation (7) is the governing equation of the acoustic fluid inside the reservoir which is known as Helmholtz equation.
If the fluid is assumed to be incompressible, the Helmholtz equation becomes the Laplace equation: For solving the Helmholtz equation, the appropriate boundary conditions are required which will be addressed in the following [1][2][3][4].

Reservoir Bed Boundary.
To take into account the interaction effects, the boundary condition of the foundationreservoir interface must be extracted. If no absorption or penetration of water is occurred by the foundation, the boundary condition at the bottom of the reservoir is similar to the boundary condition at the fluid-solid interface. But, in the case of existence of water absorption at the bottom, its corresponding boundary condition can be extracted [17]. Given that the shear strength of water is negligible, in such case, the interaction in the vertical direction can only be considered where the hydrodynamic pressure waves influences the bottom of reservoir only in the vertical direction. According to the discussions, it is clear that the hydrodynamic pressure in the water P(n, t) in a one-dimensional case can be defined using the wave equation:

Mathematical Problems in Engineering
Similarly, displacement due to the interaction of reservoir and foundation, x(n, t), in the reservoir bottom can be defined as follows: where C f � ����� E f /ρ f . E f is the modulus of elasticity, and ρ f is the density of the materials of the reservoir bottom. At the bottom of the reservoir, according to the Euler equation, the vertical pressure gradient is proportional to the total acceleration of the foundation. Hence, € x(0, t) is the acceleration of reservoir bed induced by reservoir-foundation interaction. e equilibrium of particles in the reservoir bed causes the following equation: e d'Alembert solution for equation (10) is where g f is the defining function of impacted wave amplitude that moves vertically downward in the materials of the reservoir bottom. By using this relation, it can be written as and hence, where the prime sign denotes the derivative of g f respect to (n + C f t) parameter. By differentiating equation (12) and using equation (14), the following equation is obtained: or From equations (11) and (16), where q � (ρ w /ρ f C f ) is the attenuation coefficient.
If the reservoir bed is assumed horizontal, the corresponding boundary condition is defined as follows: where L rf is the boundary of the reservoir with foundation interface. For rigid reservoir bottom C f � ∞ and q � 0. Hence, the second term of equation (18) is eliminated, and the boundary condition for the rigid bottom of the reservoir is obtained:

Reservoir Free Surface Boundary.
If the effect of the gravitational waves is neglected, the boundary condition of the free surface is defined as follows: where L f is the boundary corresponding to the free surface.

Dam-Reservoir Interface Boundary.
To apply the interaction effects at the fluid-solid interface, it is necessary that the dam acceleration be related to reservoir pressure produced in the reservoir. For this purpose, by using the force balance relation in the contact surface, the boundary condition is written as follows where ρ w , L r d , a x , and (zP/zx) are the water density, damreservoir contact boundary, outward acceleration normal to the boundary, and the hydrodynamic pressure gradient along the vector normal to the interface, respectively. Given that, this study aims to analyze the reservoir under the vertical vibration; hence, it is assumed that a x � 0, and the interface boundary condition is defined as follows:

Far-End Boundary Condition.
Given the upward propagation of the hydrodynamic wave far from the reservoir, it is necessary that the pressure is completely dissipated, and there is no reflection from the far boundary. So, where L Infinite corresponds to the far end of reservoir, in which x ⟶ ∞ in analytical solution.

Analytical Solution
By using the extracted equation for the system under the fluid-structure interaction and the boundary conditions corresponding to the domains, the response of the desired system can be obtained by the analytical solution.
To achieve an analytical solution process, the reservoir system is considered as 2D with rectangular shaped. e wave equation extracted in the previous sections can be used for the reservoir. erefore, the behavior of the water inside the reservoir can be defined by the wave equation in terms of pressure variable: where C is the acoustic wave speed, S is the reservoir domain, and P is the hydrodynamic pressure created in the reservoir. Since wave equation is a second-order differential equation, the appropriate boundary conditions are required for solving.
One way to solve the wave equation in applied mathematics is to use the separation of the variables method [18].
For this purpose, it is assumed that the pressure response is defined as the product of three independent functions, one dependent on the parameter x, and the others dependent on the parameters y and t. at means, P � P x P y P t . (25) Due to the independence of x, y, and t, second derivative of P x , P y and P t can be expressed as follows: (26) e following result is obtained by placing the equations in the governing wave equation of reservoir: So, According to equation (28) and the independence of functions, to make a relationship, (P x ″ /P x ), (P y ″ /P y ), and (P t ″ /P t ) must have a constant value. erefore, because of the dynamic nature of the problem, the following conditions can be achieved: erefore, or (31) According to the mentioned equation and considering the dynamic nature of the pressure wave and due to the propagation of the wave to the far and its dissipation, the solutions to the problem can be considered as follows: By applying the far end boundary condition of the reservoir (x ⟶ ∞⇒P � 0), A 2 must be equal to zero, and the following result is obtained: (33) By using the dam-reservoir interface condition (x � 0⇒(zP/zx) � 0), it can be written as −kP y P t A 1 � 0.

(34)
For the satisfying the above condition, it is necessary that k � 0.
So, it can be concluded that in as case, a vertical vibration is applied to the reservoir bottom; as expected, the produced hydrodynamic pressure is independent of horizontal direction. By merging the constant A 1 into other coefficients, the hydrodynamic pressure response is calculated as where λ 2 � (ω 2 /C 2 ) − k 2 . e general solution of equation (37) is as follows: where C 1 and C 2 are the constants obtained by the boundary conditions. By substituting the free surface boundary condition (y � H⇒P � 0) in equation (38), Hence, C 2 is obtained as So, by applying the free surface boundary condition, the general solution is defined as Mathematical Problems in Engineering 5 or By applying the boundary condition of the foundationreservoir interface, equation (18) to equation (42), and by assuming that the harmonic vertical acceleration is a y � a y e iωt , it can be found that Considering that P t � B 1 e − iωt + B 2 e iωt , it can be written that For satisfying the above equation, B 1 constant must be equal to zero (B 1 � 0) and so In the above equation, B 1 constant has been merged with C 1 , and C 1 is obtained as follows: So, P y is obtained as e peak values of this function occurred at frequencies of (ωH/C) � (π/2), (3π/2), (5π/2), . . .. So, Obtained response is similar to result presented by Chopra (1973 and1983) [6,19]. e specific cases of the extracted function can be corresponding to the rigid foundation and incompressible fluid.

Incompressible Fluid.
If the fluid is incompressible, then ω ⟶ 0 and λ ⟶ 0. In this case, the limit of the P function is defined as follows: when λ ⟶ 0. en, Lim(Sinλ × (H − y)/λ) � (λ× (H − y)/λ) or (H − y). erefore, the compressive response is summarized as follows: which denotes the linear distribution for pressure. Obtained solution for incompressible fluid compiles with the pressure distribution equation in the vertical direction defined in hydrodynamic is ((zP/zy) � −ρa y ).
According to this relation, it can be found that in a case where the fluid is incompressible, the pressure distribution in the reservoir under the effect of the vertical acceleration is not depending on the bottom absorption and the time.

Concrete Gravity Dam and Reservoir Model
e proposed analytical solution is investigated for a case study model, and the response for different parameters is examined. e geometric characteristic of the model, which is a concrete gravity dam with 100 m height of reservoir, is shown in Figure 1. A harmonic vertical acceleration of a y � a y e iωt with a y � 0.5 g has been applied to model. e water density and speed of compressive waves in the water have been considered as 1000 Kg/m 3 and 1438.66 m/s, respectively. To investigate the foundation effect, the modulus of elasticity and the density of foundation are assumed to be 3430 MPa and 2400 Kg/m 3 , respectively [20].

Results
Using the extracted analytical solution, the responses for rigid and flexible foundation were obtained and are presented in  For this purpose, the responses have been separated for the frequency less and more than the first natural frequency of reservoir (Ω � (πC/2H)), and the hydrodynamic pressure distribution in the dam-reservoir interface against the height is shown. It should be mentioned that y is the water depth of respect to the reservoir bottom in terms of meter. 6 Mathematical Problems in Engineering According to the presented curves, the trend of hydrodynamic pressure distribution under vertical vibration can be discussed. e first interesting case in Figure 2 is for the case in which ω � 0, that is related to the condition of incompressible fluid and indicates that the pressure distribution in this case is linear.
Also, by comparing Figure 2 with Figure 3 and Figure 4 with Figure 5, it can be seen that for a case in which the excitation frequency is more than the natural reservoir frequency (ω > Ω), unlike the case in which the excitation frequency is less than the natural reservoir frequency (ω < Ω), the maximum hydrodynamic pressure is not created at the bottom of the dam and occurred near the y � 0.3H.
Another interesting case is the case ω � 2.0Ω in Figures 3  and 5, where the pressure at the bottom of the reservoir is zero, which is due to the absorption of the excitation with the natural vibration of the reservoir. e analysis of pressure distribution in the figures shows that the effect of compressibility for the case in which the excitation frequencies are less than half of the natural frequency of the reservoir (ω < 0.5Ω) are not very important.
is has been mentioned in previous research for horizontal excitation.
In order to investigate the effect of excitation frequency on the hydrodynamic pressure distribution process, the hydrodynamic pressure sensitivity curve concerning the excitation frequency was extracted for both rigid and flexible foundation and isshown in Figures 6 and 7.
According to the curves, it can be found that the maximum responses of the system are corresponding to a case in which the excitation frequency is close to the natural frequencies of the reservoir (ω � (nπC/2H)).

Discussion
Also in some studies, to simplify the analysis, the foundation is assumed to be rigid; that assumption is not correct because the foundation under the dam and reservoir was flexible and its movement produces hydrodynamic pressure in the reservoir due to the interaction of the dam and the reservoir. In this study, the hydrodynamic response of a concrete gravity dam and reservoir under vertical excitation was investigated using an analytical solution. e results showed that the distribution of hydrodynamic pressure due to vertical vibration is not dependent on the longitudinal direction, and the propagation of hydrodynamic wave occurs only along the vertical direction. For an incompressible fluid, the pressure distribution along the dam is linear.
Based on the results, it can be found that for the excitation frequencies more than the natural frequency of the reservoir, unlike the excitation frequencies less than the natural frequency of the reservoir, the maximum hydrodynamic pressure does not occur in the bottom of the reservoir. e case in which ω � nΩ, the hydrodynamic pressure in the bottom of the reservoir is zero, which is due to absorption of excitation by the natural vibration of the reservoir. Also, the study of pressure distribution curves indicates that the compressibility effects for the excitation frequency are more than the case in which the frequency is equal to half the natural frequency of the reservoir. e results also show the sensitivity of the response to the loading frequency. According to the sensitivity analysis of the response to the excitation frequency, it can be seen that the maximum responses are related to the conditions where the excitation frequency is close to the natural frequencies of the reservoir.

Data Availability
e data used to support the findings of this study are included within the article.