Reinforced concrete water tanks which play essential functions in the entire water supply system are typical liquid-containing structures. Cracks are prone to occur when tanks suffer from strong earthquakes and may lead to liquid leakage, making liquid-containing structures losing their strength or even completely lost. The internal liquid in the tank under the seismic actions will slosh, and the distribution of the additional hydrodynamic pressure for the tank body will be directly affected by the amplitude of the liquid sloshing. To reasonably guide the seismic design of the tank in the water supply system, the liquid sloshing of the circular tank under seismic action is the focus of this research and further obtaining the law of liquid sloshing is the purpose of this research.
There are three methods to study sloshing wave height and hydrodynamic pressure, and they are theoretical calculation, numerical simulation and experimental test [1]. The theoretical analysis of the sloshing wave height for the liquid-containing is the earliest. In 1957, professor Housner adopted the approximation method [2] to derive the formula of the maximum sloshing wave height
The formula
The calculation formula for the sloshing wave height of the oil tank under seismic action is recommended in the Chinese standard for seismic design of petrochemical steel equipment (GB/T 50761-2018). The sloshing wave height is proportional to the tank type coefficient, the seismic effect coefficient and the oil tank’s inner radius. But, the formula for estimating the water tank’s sloshing amplitude under seismic action is not given in the Chinese code of seismic design of outdoor water supply, sewerage, gas and heating engineering (GB 50032-2003). The numerical simulation method is also widely used in estimating the height of liquid sloshing waves. A numerical model based on the finite volume method is established by professor Goudarzi [6] to estimate the hydrodynamic damping caused by the vertical bounded baffle. The existence of the baffle reduces the sloshing amplitude of the liquid level through comparative analysis. Professor Wang studied the sloshing wave heights of isolated rectangular concrete tanks under unidirectional and bidirectional long-period seismic actions. [7]. The horizontal displacement and sloshing wave height increased under the bidirectional long-period seismic actions than under the unidirectional seismic actions. When other conditions remain unchanged, the larger the tank’s structural size, the larger the horizontal displacement and sloshing wave height.
The shaking table test method is usually used in conjunction with the numerical simulation method. The shaking table test’s research results are used to guide the model establishment or verify the numerical model’s validity. The time history of liquid sloshing wave height is measured by using the shaking table test of the 1:20 scale model for a large steel storage tank [1]. The maximum value of the sloshing wave height is calculated by numerical simulation is close to the value that obtained from the shaking table test. Still, the shapes of the time history curve for the two both are quite different. Although the size of the oil tank is scaled down during the shaking table test, the material properties of the steel tank and internal liquid water could not be scaled down, and the test results are still quite different from the actual situation. The liquid-containing density has a significant influence on the sloshing wave height of the liquid surface, which is analysed and confirmed by Yazdanian et al. The lower the liquid density, the larger the sloshing wave height [8]. The dynamic characteristics of the rectangular tank system are sensitive to the frequency content of ground motions considering the fluid-solid coupling and soil interaction, which has been confirmed by Safaa et al. [9]. The shaking table tests of three open rectangular tanks with a scale ratio of 1:10 are completed by Jure et al. [10]. The analysis results show that when the period of input ground motion is close to the first-order period of the liquid sloshing, the sloshing amplitude of the liquid surface is huge. The wave height and distribution characteristics of liquid sloshing are directly affected by the stiffness of the tank wall, the amount of liquid, and the types of input dynamic excitation, including period, amplitude, or duration time. The oil tank models with various liquid levels and fixed roof are not tested by Bae et al. through shaking table tests [11]. The results show that the acceleration of the tank along the height direction is increased when the amount of liquid in the tank is large. That is, the sloshing effect of the liquid is significantly different when the amount of liquid in the tank is different. The dynamic behaviour characteristics, including beam-type and oval-type vibration of a cylindrical liquid-containing tank under horizontal earthquake excitation, are investigated by shaking table tests [12]. The dynamic characteristics and sloshing wave height of the vertical storage tank are analysed by using ADINA finite element software. It is found that the radius of the storage tank and the height of the storage liquid has a more significant influence on the sloshing of the liquid surface [13]. The dynamic characteristics of liquid in sloped bottom tanks are investigated by Amiya by classifying frequency contents of six different ground motions [14].
Based on the above analyses, the ground circular tank’s dynamic characteristics in water supply systems with different capacities are analysed by using theoretical analysis and numerical simulation methods. Based on the dynamic characteristics analysis, the relationship between the amplitude of liquid sloshing and ground motion characteristics, liquid-containing sloshing frequency, tank radius, liquid-containing height, and other factors are analysed. Based on the original theoretical formula, the calculation formula of the sloshing wave height for the water tank is fitted.
The three circular tanks are named as tanks A, B, and C with 500 m3 capacity, 200 m3 capacity, and 2000 m3 capacity. The three tanks have the following structural characteristics: a bottom thickness of 0.3 m, a wall thickness of 0.25 m, a total height of 3.8 m, a maximum water storage height of 3.5 m and a reinforcement diameter of 10 mm. The inner radius of tanks A, B, and C is 6.75 m, 4.3 m and 13.5 m, respectively. For each tank, the height of the water is taken as 0%, 10%, 20%, 30%, 40%, 50%, 60%, and 70% of its maximum water storage height, respectively. The corresponding height is 0 (the state of the empty tank), 0.35 m, 0.70 m, 1.05 m, 1.40 m, 1.75 m, 2.10 m and 2.45 m. This means that there are eight different storage conditions for each tank. Since the water storage height of the tank during normal operation is about 70% of the total height, the 2.45 m height is set to the maximum storage height of the tank. In order to facilitate the analysis of the subsequent results, the ‘A-50%’ is used to indicate the condition of 500 m3 storage capacity and 1.75 m storage height, and other conditions naming are analogous.
In the process of finite element analyses for ADINA software, the tank body model is simulated as a 3D-Solid element, and the structure material is simulated by using concrete material. The reinforcement is simulated by using the truss unit. The stress-strain curves of the concrete and reinforcement are shown in Figure 1. The liquid model is simulated as a 3D-Fluid element and the element is set as the linear potential fluid element for static and modal analyses, and the subsonic potential fluid element for dynamic analysis. The water properties include a density of 1000 kg/m3, a bulk modulus of 2.3×109 Pa and a damping ratio of 0.16%.
The dynamic characteristics of liquid-solid coupling system with different water storage heights are analysed to obtain the mode period and mode shape. The partial mode periods obtained by numerical calculations are shown in Table 1. The calculation results in Table 1 show that when the tank radius is the same, the higher the water storage height, the shorter the liquid sloshing period, the longer the structural vibration period, and the more significant the liquid-solid coupling effect. When the water storage height is the same, the larger the tank radius, the longer the liquid sloshing period and the structural vibration period, and the less significant the liquid-solid coupling effect. The same order modal shape of liquid sloshing is almost the same. Most of the shapes of structural vibration modes appear in pairs due to the rotational symmetry of the circular tank structure.
Modal periods of part conditions (s).
A-0% | / | / | / | 0.026 | 0.026 | 0.023 |
A-10% | 12.195 (12.500) | 4.255 | 2.695 | 0.026 | 0.026 | 0.023 |
A-30% | 7.143 (7.299) | 2.688 | 1.873 | 0.027 | 0.026 | 0.023 |
A-50% | 5.650 (5.770) | 2.364 | 1.770 | 0.028 | 0.027 | 0.024 |
A-70% | 4.926 (5.025) | 2.268 | 1.754 | 0.030 | 0.029 | 0.027 |
B-70% | 3.333 (3.472) | 1.754 | 1.370 | 0.025 | 0.024 | 0.020 |
C-70% | 9.091 (9.615) | 3.704 | 2.632 | 0.040 | 0.039 | 0.039 |
Assume that the liquid-containing structure is a rigid body, and the liquid in the tank is non-viscous, non-rotating, incompressible and slightly vibrating. The calculation formula of the fundamental natural frequency
The theoretical formula is used to calculate the first-order period of liquid-containing sloshing under various working conditions. The values are shown in the brackets of the second column in Table 1. The numerical calculation results are similar to the theoretical calculation results. It shows that the established tank model is reliable. They are the basis for the correct analyses of the seismic time-history response of liquid sloshing for water tanks.
The dynamic characteristics of the liquid-containing structure differ from that of a general building structure. The impacts of the epicentre distance and the site type on the liquid-structure coupled system must be considered in the seismic response analysis. As such, seven natural seismic motions with amplitudes of about 100 gal during the Wenchuan earthquake’s mainshock are selected for structural analysis. A gal is a unit of acceleration (cm/s2). These seven natural seismic motions are labelled according to their collection locations and directions, as follows: HSDB-EW, JZGYF-NS, PJW-NS, BJ-EW, CC-NS, YL-NS and HX-NS. The time intervals of all seismic motions are 0.005 s, and the time lengths of all natural seismic ground motions are taken as 30 s. The predominant periods of seismic ground motions are obtained by Fourier transformation. Table 2 lists details of the seismic ground motions, and Figure 2 shows the time-history curves of the HSDB-EW and HX-NS seismic ground motions which are selected as examples.
Basic information of the selected seismic ground motions.
HSDB-EW | II class | Sichuan Province | 125.9 | −102.643 | 6.160 | 0.095 |
JZGYF-NS | II class | 263.3 | 100.224 | 18.425 | 0.168 | |
PJW-NS | II class | 81.0 | 101.149 | 10.105 | 0.297 | |
BJ-EW | III class | Shanxi Province | 513.1 | 120.292 | 15.720 | 0.611 |
CC-NS | III class | 528.6 | 107.719 | 14.980 | 1.138 | |
YL-NS | III class | 569.9 | −94.005 | 16.730 | 1.862 | |
HX-NS | III class | 599.0 | 92.090 | 13.835 | 4.096 |
Calculation formulas for the sloshing wave heights of liquid surfaces have been provided by Professor Housner [2] and they are provided in the Japanese welded steel tank standard JIS B 8501 [4] and the U.S. Atomic Energy Commission TID7024 standard [3]. Because Professor Housner’s calculation formula is practical and continues to be widely used, this formula is used in this paper to perform the theoretical calculations of the sloshing wave height. This calculation formula is shown in Eq. (3):
The result of
A-10% | 0.043 | 0.092 | 0.056 | 0.057 | 0.054 | 0.165 | 0.164 | |
0.024 | 0.052 | 0.032 | 0.032 | 0.030 | 0.093 | 0.093 | ||
0.074 | 0.112 | 0.162 | 0.154 | 0.249 | 0.345 | 0.339 | ||
A-30% | 0.055 | 0.137 | 0.143 | 0.088 | 0.268 | 0.977 | 0.394 | |
0.031 | 0.077 | 0.081 | 0.049 | 0.151 | 0.552 | 0.223 | ||
0.101 | 0.163 | 0.233 | 0.201 | 0.466 | 1.504 | 0.586 | ||
A-50% | 0.064 | 0.154 | 0.370 | 0.261 | 0.441 | 0.559 | 0.807 | |
0.036 | 0.087 | 0.209 | 0.148 | 0.249 | 0.316 | 0.456 | ||
0.109 | 0.203 | 0.487 | 0.348 | 0.589 | 0.736 | 0.946 | ||
A-70% | 0.082 | 0.148 | 0.377 | 0.344 | 0.609 | 1.020 | 0.721 | |
0.046 | 0.083 | 0.213 | 0.194 | 0.344 | 0.576 | 0.407 | ||
0.115 | 0.151 | 0.461 | 0.393 | 0.860 | 1.259 | 0.910 | ||
B-70% | 0.149 | 0.142 | 0.497 | 0.307 | 0.335 | 0.944 | 1.233 | |
0.054 | 0.051 | 0.179 | 0.111 | 0.121 | 0.340 | 0.444 | ||
0.113 | 0.114 | 0.398 | 0.220 | 0.484 | 0.943 | 1.075 | ||
C-70% | 0.062 | 0.111 | 0.074 | 0.061 | 0.104 | 0.377 | 0.261 | |
0.070 | 0.125 | 0.084 | 0.068 | 0.118 | 0.426 | 0.295 | ||
0.158 | 0.259 | 0.267 | 0.223 | 0.275 | 0.958 | 0.717 |
According to the calculated dynamic coefficient
The unidirectional horizontal seismic actions in the liquid-structure coupled tanks are analysed under X-direction input with modulated peak amplitudes of 100 gal, 200 gal and 400 gal. The distribution of the sloshing wave height is analysed by collating the calculation results. The sloshing wave height indicates the liquid surface’s sloshing degree, which affects the hydrodynamic pressure distribution directly.
The maximum and minimum values of the vertical displacement for the liquid surface are equal at each moment of the unidirectional horizontal seismic actions with different peak accelerations. When only the peak acceleration of the seismic ground motion is changed while the other conditions remain unchanged, the liquid surface’s sloshing wave height increases with increasing peak acceleration in multiples. The growth multiple is the seismic coefficient
Figure 4 shows the time histories of the maximum sloshing wave height for a liquid surface under long- and short-period seismic ground motions. For the same seismic ground motion and different water storage heights, the liquid surface’s maximum sloshing wave height has no apparent change in law. The sloshing responses at different storage heights differ due to the other seismic motion inputs. Figures 4a and 4b show that the liquid in the tank sloshes more violently under long-period seismic ground motion than under short-period seismic ground motion. The maximum sloshing wave heights of different working conditions under the seismic actions with 100 gal peak acceleration are calculated using ADINA and are listed in the
From the analysis results of the sloshing wave height, it can be observed that the plane positions of the maximum and minimum sloshing wave heights under different seismic ground motions are near the centre or the circumferential radius, and their spatial positions indicate the direction of the ground motion input. The first three-order modal shape of liquid sloshing under different conditions are similar. When the tank radius is the same, and the water storage heights are different, the sloshing shapes of the liquid surface in the maximum wave height moment are basically the same. When the water storage height is the same and the tank radii are different, the sloshing shapes of the liquid surface in the maximum wave height moment are also basically the same under long-period seismic ground motions. In addition, the position of the maximum wave height is usually at the circumferential radius, and the sloshing amplitude of the liquid surface is larger than that under short-period seismic ground motions. Figure 5 shows plots of the theoretical and ADINA calculation results for the maximum sloshing wave height of the liquid surface in each condition.
It can be seen from Figure 5 that the maximum wave height calculation results obtained using ADINA numerical simulation are larger than the theoretical calculation results regardless of the water storage height, the tank radius, or the long- and short-period seismic ground motions, but the changing trends are basically the same. Under short-period seismic motion actions, the wave heights obtained by the ADINA numerical simulation are not much different from the theoretical results. Due to the large periodic component of seismic ground motion, large-amplitude sloshings of liquid surfaces occur under long-period seismic actions. The numerical simulation results are much larger than the theoretical results, i.e. the maximum ratio between the two approaches eight. The calculation formula for the liquid-surface sloshing wave height derived by Housner assumes that the liquid-containing structure is rigid and the liquid surface exhibits a linear small-amplitude motion. However, the actual liquid sloshing will not be completely linear. For a large-amplitude motion of the liquid surface under long-period seismic actions, the nonlinear effect is enhanced. Therefore, using Eq. (3), the maximum sloshing wave of a free surface is predicted to be unsafe. The sloshing amplitude of the liquid surface under long-period seismic actions is underestimated.
The maximum sloshing wave height of the liquid surface under seismic ground motions has significant relationships with seismic coefficient, ground motion characteristic, tank radius and liquid sloshing frequency. The calculation formula of the maximum sloshing wave height for the liquid-containing structure is fitted based on the ADINA calculation results and the theoretical calculation formula. The maximum sloshing wave height data under unidirectional seismic actions with 100 gal peak acceleration are used to fit the formula with the above analyses. The
It can be seen from Figure 6 that the maximum sloshing height
Eq. (4) is the formula for calculating the maximum sloshing wave height of the liquid surface, for which the linear correlation coefficient is 0.956, and the standard deviation is 0.063.
The effect of the damping ratio on liquid sloshing is not negligible. The smaller is the damping ratio, the larger is the dynamic amplification factor of the sloshing vibration. The damping of the liquid is much smaller than that of the structure. The damping-ratio correction coefficient
Damping ratio correction coefficient
Correction coefficient |
0.44 | 0.56 | 0.78 | 1.00 | 1.18 | 1.32 | 1.53 | 1.79 |
The dynamic amplification factor corresponding to a 0.5% damping ratio is 1.79 times as shown in Table 4. The sloshing damping ratio of water measured by Tianjin University in China in a simulation test was 0.16%. The response spectra of different damping ratios recorded during several typical earthquakes are compared in the text by compile a group of earthquake engineering. The zero damping response spectrum was found to be about 1.5 to 3 times that of the damping ratio of 5%. To obtain the correction coefficient for a 0.16% damping ratio, data with damping ratios <5%, as shown in Table 4, were fitted, and the following fitting Eq. (5) obtained:
Eq. (5) calculated that the correction factor for the 0.16% damping ratio is 1.931. The zero damping ratio’s correction factor is 2.02, which is consistent with the conclusions reported in the literature.
According to the above analyses, the 5% damping ratio long-period response spectrum shown in Figure 7 can be used to design the sloshing wave height of a water storage tank under unidirectional ground motion. Simultaneously, the response spectrum value can be corrected according to the water sloshing damping ratio. That is, the damping coefficient directly increases 1.931 times based on the 5% response spectrum value. Thus, Eq. (4) becomes Eq. (6).
When using Eq. (6) to calculate the sloshing wave height of a liquid surface for a liquid-containing structure under unidirectional seismic ground motion action,
The article R7.1 in the ACI 350.3-01 code [16] gives the calculation formula of sloshing wave height dmax of circular tank subjected to seismic action, as shown in Eq. (7).
Seismic zone factor
Factor |
0.075 | 0.15 | 0.2 | 0.3 | 0.4 |
Soil profile coefficient
A | A soil profile with either: (a) a rock-like material characterised by a shear wave velocity >2500 ft/s (762 m/s), or by other suitable means of classification; or (b) medium-dense to dense or medium-stiff to stiff soil conditions where the soil depth is <200 ft (60,960 mm) | 1.0 |
B | A soil profile with predominantly medium-dense to dense or medium-stiff to stiff soil conditions, where the soil depth exceeds 200 ft (60,960 mm) | 1.2 |
C | A soil profile containing >20 ft (6096 mm) of soft to medium-stiff clay but not >40 ft (12,192 mm) of soft clay | 1.5 |
D | A soil profile containing >40 ft (12,192 mm) of soft clay characterised by a shear wave velocity <500 ft/s (152.4 m/s) | 2.0 |
Importance factor
Tanks containing hazardous materials | 1.5 |
Tanks that are intended to remain usable for emergency purposes after an earthquake or tanks that are part of lifeline systems | 1.25 |
All other tanks | 1.0 |
From Eq. (7), it can be seen that the maximum wave height of liquid sloshing in liquid-containing structures calculated by the American code takes into account the structure size of the tank, seismic map zone, site profile, structural importance and the spectral amplification effect corresponding to the first sloshing period. Eq. (7) covers various types of liquid-containing structures, and Eq. (6) established in this paper is only for reinforced concrete liquid containing tanks in water supply systems. In Eq. (7), the influence of seismic action on liquid sloshing is divided into three aspects of parameters
Calculation results of
0.034 | 0.088 | 0.134 | 0.171 | 0.211 | 0.114 | |
0.013 | 0.037 | 0.059 | 0.078 | 0.109 | 0.046 | |
2.626 | 2.360 | 2.252 | 2.191 | 1.941 | 2.476 |
It can be seen from the calculation results in Table 8 that the maximum sloshing wave height of the liquid surface obtained by both calculation methods has the same order of magnitude, and
The dynamic characteristics and seismic responses of the ground-rested circular reinforced concrete tank are researched. Based on the calculation results of the sloshing responses with different water storage heights and different tank radii, the calculation formula for the maximum sloshing wave height of a liquid surface under unidirectional horizontal seismic actions is fitted, and the obtained formula has a good estimation effect for long- and short-period seismic actions. The maximum sloshing wave height of a liquid surface is proportional to the tank radius, the ground motion amplification factor and the seismic coefficient. Combining the long-period seismic design spectrum and the damping ratio correction coefficient, the calculation method for the maximum sloshing wave height is proposed based on the seismic design response spectrum 5% damping ratio. Because of the current research comparison, the Chinese code is more conservative than the American code in calculating the sloshing wave height of liquid-containing.