Dynamic behaviors of concrete gravity dam against combined blast wave and bubble pulsation of underwater explosion

Abstract

Underwater explosion process mainly includes the preceding blast wave and the succeeding bubble pulsation, while the effect of bubble pulsation on hydraulic structures is always neglected. This paper aims to numerically study the dynamic behaviors of prototype concrete gravity dam against underwater explosion, considering the combined effects of blast wave and bubble pulsation. Firstly, a finite element analysis (FEA) approach for predicting the underwater explosion loadings, as well as the damage evolution and dynamic responses of dam against underwater explosions was proposed, in which the hydrostatic pressure of reservoir and the dead load of dam were considered. Based on Cole's formula and the existing underwater explosion tests, the reliability of FEA approach, including the mesh sizes, material models and the corresponding parameters, as well as the mapping, Fluid-Structure Interaction and contact algorithms, was comprehensively validated. Secondly, to realize the simulation of the prototype gravity dam impeded by the validated small Lagrange mesh size, a mesh transition strategy was proposed based on the equivalence of dam damage modes and total eroded element mass. Finally, the dynamic behaviors of a typical 120 m-height concrete gravity dam against underwater explosion, i.e., 1 t of TNT equivalence detonated at a depth of 35 m and standoff distances of 4–16 m, were examined. The dimensionless distance λR, i.e., the ratio of the standoff distance to the maximum bubble radius, was introduced to quantitatively evaluate the underwater explosion loadings and dynamic behaviors of the dam. It indicates that, the blast wave firstly leads to the local compressive damage on the upstream surface, as well as the tensile damage on the downstream folded slope and dam heel. The succeeding free vibration and bubble pulsation will further aggravate the damage degree, resulting in the cracking of dam neck and heel. Especially for the scenarios when λR ≤ 1.22, the overall sliding of dam is highly prone to occur.