Dynamic analysis of rectangular tank using response spectra

. 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 one-sided 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 fluid-structure interaction (FSI) was the matter of interest here. Deformations and bending moments of walls were obtained.


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 1998-1 [2], in each area of Slovakia, the values of basic seismic acceleration are defined ( = 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.

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 ( ) and mode shapes are calculated.For a given direction of excitation (, , ), the participation factor of each mode shape  (1) is calculated.It expresses the contribution of each mode shape to individual global directions of excitation:  =  . ( where  -normalised th mode shape,  -mass matrix,  -vector describing the direction of excitation. To calculate the modal factor (  ) we proceed with the relation Eq. (2).Subsequently, according to the relationships given in the EN standard [2], the corresponding acceleration ( ) VIBROENGINEERING PROCEDIA.APRIL 2019, VOLUME 23 can be calculated for each mode shape and its corresponding frequency: In the listed example, 2 directions of construction excitation were used in the  and  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.

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 ( = 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 [3][4][5] and [6] are also dealing with the interaction of structure and subsoil.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  = 1.10 m/s 2 .

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.

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 results are summarized in Tables 2 and 3.For the location where the basic seismic acceleration is  = 1,10 m/s 2 (Fig. 2) and the seismic excitation is in the direction of -axis, there are deformations of the tank and bending VIBROENGINEERING PROCEDIA.APRIL 2019, VOLUME 23 moments of the tank wall shown in the Fig. 5-Fig.8. Other variables and values are listed in the Table 2.In the deformation figures, it can be seen that if the load is applied in the direction of -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 half-wave.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 -axis then maximum bending moment  is in middle of the bottom edge (near clamp) of the walls which are perpendicular to the load direction (longer walls).Maximum bending moment  is in the top corners of the tank.For the location where the basic seismic acceleration is  = 1,10 m/s 2 (Fig. 2) and the seismic excitation is in the direction of -axis, there are deformations of the tank and bending moments of the tank wall shown in the Fig. 9-Fig.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 -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 half-wave.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 -axis then maximum bending moment  is in middle of the bottom edge (near clamp) of the walls which are perpendicular to the load direction (shorter walls).Maximum bending moment  is in the middle of the top edge of shorter walls of the tank.

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.

Fig. 5 .Fig. 6 .
Fig. 5. Deformation in the direction of -axis by seismic excitation in the direction of -axis with the maximum value 4,785 mm (red surface)

Fig. 7 .Fig. 8 .
Fig. 7. Bending moments  as a result of seismic excitation in the direction of -axis with the maximum value 14,561 kNm/m (red surface) Fig. 8. Bending moments  as a result of seismic excitation in the direction of -axis with the maximum value 28,021 kNm/m (red surface)

Table 1 .
Natural frequencies and participation factors corresponding to global mode shapes Mode shape no.Frequency  Participation factor  Direction

Table 3 .
Results corresponding to excitation of structure in the direction of -axis