Dynamic analysis of rectangular tank using response spectra
Lenka Uhlirova^{1} , Norbert Jendželovský^{2}
^{1, 2}Faculty of Civil Engineering STU in Bratislava, Bratislava, Slovakia
^{1}Corresponding author
Vibroengineering PROCEDIA, Vol. 23, 2019, p. 99104.
https://doi.org/10.21595/vp.2019.20657
Received 13 March 2019; accepted 20 March 2019; published 25 April 2019
37th International Conference on Vibroengineering in Bratislava, Slovakia, April 2526th, 2019
. Reinforced concrete tanks are currently used to store various substances and products. Most often they serve as drinking water reservoirs, to sewage treatment or to store various technical fluids. Although rectangular tanks are less used, they have advantages, such as less sensitivity to onesided loading and better utilization of space when using a group of tanks. Rectangular tank which is the matter of interest in this paper is for the purposes of analysis considered as fulfilled with water. Dynamic analysis of the tank considering fluidstructure interaction (FSI) was the matter of interest here. Deformations and bending moments of walls were obtained.
Keywords: rectangular reinforced concrete tank, mode shape, response spectra.
1. Introduction
To solve tasks related with continuum mechanics, three numerical methods are currently used (finite element method, finite band method and boundary element method). The finite element method (FEM) is the most widely used method thanks to the development in the field of computer programming. The FEM is based on variational principles. Software ANSYS Academic was used to solve tasks included in this paper.
In this paper rectangular concrete tank was modelled included water filling. Interaction between the fluid and tank was also considered in the calculation. The subsoil category C was chosen to include subsoil effects in FSI analysis [1, 2]. According to STN EN 19981 [2], in each area of Slovakia, the values of basic seismic acceleration are defined (${a}_{gR}=$ 0,40 m/s^{2}, 0,63 m/s^{2}, 0,86 m/s^{2}, 1,10 m/s^{2}).
A response spectrum method was used to obtain results from the dynamic loading of the tank. From the obtained results, the most interesting variables are deformation of the tank and bending moments of the tank walls.
2. Response spectra method
The response spectra method is described in detail in [1]. This method is very often used for the tasks considering seismic load.
An integral part of response spectra method is the modal analysis where natural frequencies (${f}_{i}$) and mode shapes are calculated. For a given direction of excitation ($x$, $y$, $z$), the participation factor of each mode shape ${\gamma}_{i}$ (1) is calculated. It expresses the contribution of each mode shape to individual global directions of excitation:
where ${\mathbf{\Phi}}_{i}$ – normalised $i$th mode shape, $\mathbf{M}$ – mass matrix, $\mathbf{d}$ – vector describing the direction of excitation.
To calculate the modal factor (${A}_{i}$) we proceed with the relation Eq. (2). Subsequently, according to the relationships given in the EN standard [2], the corresponding acceleration (${S}_{ai}$) can be calculated for each mode shape and its corresponding frequency:
In the listed example, 2 directions of construction excitation were used in the $x$ and $y$ directions. The results of the modal analysis are listed in the Chapter 1.3. and results of spectral analysis are listed in the Chapter 1.4.
2.1. Model of the tank
The reservoir itself was modelled as rectangular, with the ground plan dimensions of 11 m×7,5 m and height of 4,5 m. Wall thickness is 200 mm and thickness of the bottom slab is 500 mm. The water column reaches a height of 4 m. The whole tank is on the ground surface (overground tank) and it is modelled as reinforced concrete structure with properties of the C25/30 concrete ($E=$ 31 GPa). Foundation slab is considered as bonded to the subsoil in the level of the terrain.
To take into account the interaction of the structure and the liquid filling, there were inserted the intermediate elements (CONTAC52 Fig. 1(c)) between the tank walls (SHELL181 Fig. 1(a)) and fluid elements (FLUID80 – Fig. 1(b)). Mentioned intermediate element (CONTAC52) takes into account the FSI effects. References [35] and [6] are also dealing with the interaction of structure and subsoil.
Fig. 1. a) Rectangular 3D shell element SHELL181 [7], b) fluid 3D element FLUID80 [7], c) contact element CONTAC52 [7]
a)
b)
c)
One of the response spectra used in the calculation is shown in Fig. 2. This is the response spectrum corresponding to the subsoil category C, area with seismic acceleration ${a}_{gR}=$ 1.10 m/s^{2}.
3. Results of the modal analysis
By the modal analysis, mode shapes, natural frequencies and participation factors were obtained.
It was necessary to separate mode shapes, where only the water excites from the global shapes, where the tank structure excites together with filling. In the Fig. 3. and Fig. 4. first two mode shapes corresponding to global excitation of the structure + filling are shown. In Table 1, natural frequency and participatory factor are listed for 6 global mode shapes.
Fig. 2. Response spectrum – type of spectrum 1
Fig. 3. Mode shape No. 1 ($f=$ 5,920 Hz, $\gamma =$ 4,42)
Fig. 4. Mode shape No. 2 ($f=$ 6,216 Hz, $\gamma =$ 10,68)
Table 1. Natural frequencies and participation factors corresponding to global mode shapes
Mode shape no.  Frequency $f$  Participation factor $\gamma $  Direction 
1  5,920 Hz  4,42  $Z$ 
2  6,216 Hz  10,68  $X$ 
3  8,699 Hz  8,51  $Y$ 
4  9,039 Hz  5,77  $Z$ 
5  12,133 Hz  0,00  – 
6  12,475 Hz  4,60  $Y$ 
4. Results of spectral analysis
After the modal analysis, spectral analysis follows. The object of spectral analysis is to obtain effect of seismic load on analysed structure. Response spectrum was chosen as a type 1 (Fig. 2), which parameters such as distance of seismic epicentre from the structure are listed in [2].
The calculation was performed for all four seismic areas in Slovakia. In each area, the values of basic seismic acceleration are different (${a}_{gR}=$ 0,40 m/s^{2}, 0,63 m/s^{2}, 0,86 m/s^{2}, 1,10 m/s^{2}).
The results are summarized in Tables 2 and 3.
Table 2. Results corresponding to excitation of structure in the direction of $x$axis
${a}_{gr}$ [m/s^{2}]  ${u}_{x}$ [mm]  ${u}_{y}$ [mm]  $u$ [mm]  ${m}_{x}$ [kNm/m]  ${m}_{y}$ [kNm/m] 
0,40  1,741  0,229  2,833  5,297  10,192 
0,63  2,741  0,360  4,461  8,340  16,049 
0,86  3,740  0,492  6,088  11,382  21,902 
1,10  4,785  0,629  7,789  14,561  28,021 
For the location where the basic seismic acceleration is ${a}_{gR}=$ 1,10 m/s^{2} (Fig. 2) and the seismic excitation is in the direction of $x$axis, there are deformations of the tank and bending moments of the tank wall shown in the Fig. 5Fig. 8. Other variables and values are listed in the Table 2.
Fig. 5. Deformation in the direction of $x$axis by seismic excitation in the direction of $x$axis with the maximum value 4,785 mm (red surface)
Fig. 6. Deformation in the direction of $y$axis by seismic excitation in the direction of $x$axis with the maximum value 0,629 mm (yellow surface)
In the deformation figures, it can be seen that if the load is applied in the direction of $x$axis, a deformation in the form of simple wave arises on the walls in the direction of this action (shorter walls). On the other hand, deformation of the perpendicular walls with respect to the load direction (longer walls) takes the form of the halfwave.
Fig. 7. Bending moments ${m}_{x}$ as a result of seismic excitation in the direction of $x$axis with the maximum value 14,561 kNm/m (red surface)
Fig. 8. Bending moments ${m}_{y}$ as a result of seismic excitation in the direction of $x$axis with the maximum value 28,021 kNm/m (red surface)
Based on the figures of bending moments of tank walls it is noteworthy to say that if the load is applied in the direction of $x$axis then maximum bending moment ${m}_{y}$ is in middle of the bottom edge (near clamp) of the walls which are perpendicular to the load direction (longer walls). Maximum bending moment ${m}_{x}$ is in the top corners of the tank.
Table 3. Results corresponding to excitation of structure in the direction of $y$axis
${a}_{gr}$ [m/s^{2}]  ${u}_{x}$ [mm]  ${u}_{y}$ [mm]  $u$ [mm]  ${m}_{x}$ [kNm/m]  ${m}_{y}$ [kNm/m] 
0,40  0,275  0,780  1,296  3,241  5,976 
0,63  0,434  1,229  2,041  5,104  9,411 
0,86  0,592  1,677  2,785  6,965  12,842 
1,10  0,757  2,145  3,563  8,911  16,430 
For the location where the basic seismic acceleration is ${a}_{gR}=$ 1,10 m/s^{2} (Fig. 2) and the seismic excitation is in the direction of $y$axis, there are deformations of the tank and bending moments of the tank wall shown in the Fig. 9Fig. 12. Other variables and values are listed in the Table 3.
In the deformation figures, it can be seen that if the load is applied in the direction of $y$axis, a deformation in the form of simple wave arises on the walls in the direction of this action (longer walls). On the other hand, deformation of the perpendicular walls with respect to the load direction (shorter walls) takes the form of the halfwave.
Fig. 9. Deformation in the direction of $x$axis by seismic excitation in the direction of $y$axis with the maximum value 0,757 mm (red surface)
Fig. 10. Deformation in the direction of $y$axis by seismic excitation in the direction of $y$axis with the maximum value 2,145 mm (red surface)
Fig. 11. Bending moments ${m}_{x}$ as a result of seismic excitation in the direction of $y$axis with the maximum value 8,911 kNm/m (red surface)
Fig. 12. Bending moments ${m}_{y}$ as a result of seismic excitation in the direction of $y$axis with the maximum value 16,430 kNm/m (red surface)
Based on the figures of bending moments of tank walls it is noteworthy to say that if the load is applied in the direction of $y$axis then maximum bending moment ${m}_{y}$ is in middle of the bottom edge (near clamp) of the walls which are perpendicular to the load direction (shorter walls). Maximum bending moment ${m}_{x}$ is in the middle of the top edge of shorter walls of the tank.
5. Conclusions
In situations where a liquid filling is a part of the structure, it should be taken into account in the analyses not only as a static load. It is appropriate to use 3D modelling with consideration of interaction between the structure and the liquid.
In dynamic analysis, modal analysis where mode shapes and natural frequencies are calculated in prior to spectral analysis. Subsequently, the location and subsoil parameters of the structure must be known. The higher the subsoil category, or the seismic acceleration value, the greater the effects of seismic load on the structure itself.
Acknowledgements
This paper was supported by Grant Agency VEGA, project No. 1/0412/18.
References
 Sokol M., Tvrdá K. Dynamics of Structures. Slovak University of Technology, Bratislava, 2011, p. 212. [Search CrossRef]
 STN EN 19981 Eurocode 8: Design of Structures for Earthquake Resistance. Part 1: General Rules, Seismic Actions and Rules for Buildings (STN 73 0036). SÚTN, Bratislava, 2005, p. 229 + NA/Z2. [Search CrossRef]
 Kotrasová K., Kormaníková E. A case study on seismic behavior of rectangular tanks considering fluidstructure interaction. International Journal of Mechanics, Vol. 10, 2016, p. 242252. [Search CrossRef]
 Kotrasová K., Kormaníková E. Effect of fluid in the moving container. Journal of Numerical analysis, Industrial and Applied Mathematics, Vol. 1, 2006, p. 19. [Search CrossRef]
 Kotrasová K., Harabinová S., Panulinová E., Kormaníková E. Seismic analysis of cylindrical liquid storage tanks considering of fluidstructuresoil interaction. Advances and Trends in Engineering Sciences and Technologies Proceedings of the International Conference on Engineering Sciences and Technologies, 2016, p. 8792. [Search CrossRef]
 Kala J., Salajka V., Hradil P. Response of water tower on wind induced vibration considering interaction of fluid and structure. 2nd International Conference on Engineering and Technology Innovation, Taiwan, 2012. [Search CrossRef]
 ANSYS Help, Theory Reference for the Mechanical APDL and Mechanical Applications. ANSYS manual. [Search CrossRef]
Cited By
INTERNATIONAL CONFERENCE OF NUMERICAL ANALYSIS AND APPLIED MATHEMATICS ICNAAM 2019
Kamila Kotrasová, Slávka Harabinová, Eva Kormaníková, Eva Panulinová

2020

MATEC Web of Conferences
Kamila Kotrasova, Eva Kormanikova, Iveta Hegedusova, J. Melcer, K. Kotrasová, I. Major, A.V. Benin

2020
