The Study of Seismic Response on Accelerated Contained Fluid

Ground-supported cylindrical tanks are strategically very important structures used to store a variety of liquids.This paper presents the theoretical background of fluid effect on tank when a fluid container is subjected to horizontal acceleration. Fluid excites the hydrodynamic (impulsive and convective) pressures, impulsive and convective (sloshing) actions. Seismic response of cylindrical fluid filling tanks fixed to rigid foundationswas calculated for variation of the tank slenderness parameter.The calculating procedure has been adopted in Eurocode 8.


Introduction
Liquid-containing cylindrical tanks are strategically very important structures used in water distribution systems and in industries for storing water for firefighting, petroleum, oil and liquefied natural gas, chemical, toxic, flammable, and radioactive liquids.These tanks are mainly of two types: ground-supported cylindrical tanks (concrete, prestressed concrete, and steel) and elevated tanks [1][2][3].
Seismic analysis of liquid-containing tanks requires special consideration and is different from analysis of buildings [4][5][6][7].The fluid inside of tanks exerts hydrodynamic forces on the tank walls and base.The knowledge of liquid hydrodynamic pressures developed during an earthquake is important for design of tanks, because stresses of the tank wall depend on distribution of loading along the wall height.The first solution structure of a problem was by Westergaard (1933), who determined the pressures on a rectangular dam subjected to horizontal acceleration.Hopkins and Rodriquez (1934) gave the first report on analytical and experimental observations of rigid rectangular tanks under a simulated horizontal earthquake excitation.Jacobsen analyzed, in 1949 and 1951, a rigid cylindrical liquid-containing tank and cylindrical pier surrounded by liquid, subjected to horizontal acceleration.Housner (1963) derived the expressions for distribution of hydrodynamic pressure on a rigid tank wall due to lateral base excitation and considered impulsive and convective components of hydrodynamic pressure.Veletsos (1984) obtained also the distribution of hydrodynamic pressure on rigid as well as flexible wall.
Historically, mechanical models were developed for tank with rigid walls.Housner [8] simplified the method of hydrodynamic analysis in system tank-liquid and introduced the concept of two dynamic components: dynamic impulsive and convective masses.Housner established spring-mass system, which simulates the impulsive and convective mode of vibration of tank-fluid system modified by Graham and Rodriquez (1952) and Epstein (1976).Wozniak and Michell (1978) generalized the Housner model for short and slender tanks.Velestos and Yang (1977) used a different approach to arrive at a similar type of mechanical model for circular rigid tanks.Haround and Housner (1981) and Velestos (1984) developed mechanical model for flexible tanks.Malhotra et al. [9] simplified the flexible model of Velestos (1984) and studied the dynamic behavior of unanchored and baseisolated liquid cylindrical storage tanks [10,11].
The seismic design of circular tanks is recommended by various codes of practices; their implementation strategy is rather varied leading to significantly different design forces in some cases [12].Of the best structural engineering design codes that tackle fluid tank systems are [12][13][14] the American Concrete Institute, ACI 350.3, Eurocode 8 [15], the Standards Association of New Zealand, NZS, and Indian code IS 1893-1984 [16].These codes address ground-supported European Committee for Standardization prepared code Eurocode 8 (2006).Part 4 [15] of this code is related to tanks, silos, and pipelines.Eurocode 8 uses approach of Veletsos (1984) to get hydrodynamic pressure distribution in circular tanks.Eurocode 8 mentions mechanical model of Veletsos and Yang (1977) as an acceptable procedure for rigid circular tanks.Models of Haroun and Housner (1981) and Veletsos (1984) are described along with the procedure of Malhotra et al. (2000) [9] for flexible circular tanks.Model of Housner (1963) is suggested for rigid rectangular tanks.
In regard to the combination procedure of the impulsive forces and the convective forces, Eurocode 8 suggests the absolute sum combination procedure of the impulsive forces and the convective forces, while the ACI, NZS, and IS suggest the Square Root Sum of the Squares method (SRSS) [14,16,17].
In this paper, analysis for the seismic design of fluid filling in unanchored rigid vertical circular tanks on-ground, fixed to the foundation, presented by Eurocode 8 is provided [5,15].

Seismic Analysis of Liquid-Filled Tanks
Under earthquake loading, container walls and liquid are subjected to horizontal acceleration.The seismic load acting on wall and bottom of cylindrical tanks (Figure 1) can be divided into the following components: (i) The rigid impulsive component, caused by the inertia of the liquid, if the rigid tank moves together with the foundation (ii) The convective load component, the fluid vibration in the rigid tank (sloshing) (iii) The impulsive flexible tank shell (e.g., steel tanks) with the liquid [12].
The motion of contained fluid in unanchored rigid vertical circular tanks on-ground, fixed to the foundation tank, may by expressed as the sum of two separate contributions, called "rigid impulsive" and "convective," respectively.The "rigid impulsive" component satisfies exactly the boundary conditions at the walls and the bottom of the tank but, incorrectly, due to the presence of the waves in the dynamic response, gives zero pressure at the original position of the free surface of the fluid.The second "convective" term does not alter those boundary conditions that are already satisfied, while fulfilling the correct equilibrium condition at the free surface.Use is made of a cylindrical coordinate system, , , , with origin at the center of the tank bottom, with  vertical axis;  is the radius of fluid filling; and H is the original height of the free surface of fluid; see Figure 1.The mass density of the fluid is ,  = / is nondimensional radius, and  = / is dimensionless coordinate [15].

Rigid Impulsive Component.
The seismic loads acting on wall and bottom of cylindrical tanks (Figure 1) can be divided into the following components.
The spatial-temporal variation of the "rigid impulsive" pressure is given by the expression where ]  = (2 + 1)/2 and  = /,  1 (⋅) is the modified Bessel function of order 1, and   1 (⋅) is derivate that can be expressed in terms of modified Bessel function of orders 0 and 1   1 () =  1 ()/ =  0 () −  1 ()/.The function   gives the distribution along the height of the timedependence pressure   . is angle of circumference,  = / is slenderness parameter, and   () is the free-field ground horizontal acceleration of ground with peak value denoted by   .The result of an equivalent single-degree-of-freedom system is   with the impulsive period   .The distribution   along the height of tank wall gives the function   (see ( 2)) for  = 1 (i.e., at the wall of the tank) and cos  = 1 (i.e., in the plane of the horizontal seismic action).Figure 2 presents the vertical distribution of impulsive pressure   on the tank wall along the wall height, normalized to   , for seven values of the slenderness parameter  = /.
For a number of purposes it is useful to evaluate the horizontal resultant of the pressure at the base of the wall, the total base shear, and the moments of the pressures with respect to an axis orthogonal to the direction of the seismic motion.  is the total impulsive (bending) moment, with respect to an axis orthogonal to the direction of the seismic action motion, immediately above the tank bottom and includes only the contributions of the pressures on the walls (see (7)). *  is the total impulsive (overturning) moment immediately below the tank bottom and includes the contributions of the pressures on the walls and of those on the bottom (see (5)).
The horizontal impulsive pressure resultant, total impulsive base shear, at the base of the wall is where   is called impulsive mass and indicates the mass of the contained fluid moved together with the walls and   is given by the expression where  =  2  is total mass of the fluid contained in tank.
The quantities   are plotted in Figure 3 as function of the tank slenderness parameter  = /.
The total impulsive base moment immediately below the tank bottom is given as where and total impulsive base moment immediately above the tank bottom is where The quantities ℎ  and ℎ  * are plotted in Figure 4 as function of the tank slenderness parameter  = /.

Convective Component.
The spatial-temporal variation of the "convective" (sloshing) pressure component is given by the expression where  so and damping ratio appropriate for the sloshing of the fluid.
Only the first oscillating, or sloshing, mode and frequency of the oscillating liquid ( = 1) need to be considered in expression for design purposes in which the oscillation yields periods of the order for the usual values of  are few seconds.
Figure 5 shows the functions of sloshing pressure distributions along the height of tank for the tank slenderness parameters  = /: 0.5, 1.0, and 3.0, the first three modes.
The horizontal resultants of the convective pressure are the total base shear and the total convective moments.The total convective base shear is given as where is the nth modal convective mass and ℎ  is its height; see Figure 10.Figures 6 and 7 show the quantities   as function of the tank slenderness parameter  = /: (i) for i = 1, 2, 3, and 4 in Figure 6 (ii) only for i = 2, 3, and 4 in Figure 7.
The total convective (overturning) moment immediately below the bottom plate of the tank is where The total convective (bending) moment in tank wall immediately above the bottom plate is where The quantities ℎ  and ℎ  * as function of the tank slenderness parameter  = / are seen in Figures 8 and  9, the quantities ℎ  in Figure 8 and the quantities ℎ  * in Figure 9.
The convective component of the response may be obtained from that of oscillators having masses   , attached to the rigid tank through springs having stiffness   =  2    .The tank is subjected to the ground acceleration time-history   () and the masses responds with accelerations   ().ℎ *  or ℎ  is the level where the oscillator needs to be applied in order to give the correct value of  *  or   , respectively.

Mechanical Model.
The dynamic analysis of a liquidfilled tank may be carried out using the concept of generalized single-degree-of-freedom (SDOF) systems representing the impulsive and convective modes of vibration of the tankliquid system as shown in Figure 10.For practical applications, only the first convective mode of vibration needs to be considered in the analysis of mechanical model.

Numerical Example and Results
In this study fluid filling of ground-supported cylindrical rigid tank is analyzed.The characteristics of fluid filling are  = 2 m and R, inner radiuses of cylindrical tanks, where R depended on tank slenderness parameter  = /.For the tank slenderness parameter  = 0.3 is the inner radius of tank R 6.667 m, (i) for  = 0.5 R is 4 m, (ii) for  = 0.7  = 2.857 m, (iii) for  = 1  = 2 m, (iv) for  = 2  = 1 m, (v) for  = 3, R is given 0.667 m.
The material characteristic of fluid filling (H 2 O) is density   = 1,000 kg/m 3 .As the excitation input we consider horizontal earthquake load given by the accelerogram of the earthquake in Loma Prieta, California (18.10.1989), Figure 11.
The elastic response spectra for the Loma Prieta accelerogram in California were used for calculation of seismic response in the given examples, cylindrical fluid filling.The elastic response spectra drawn for damping 0.5%, 2%, 4%, 5%, 7% can be seen in Figure 12.The impulsive spectral accelerations are obtained from a 5% damped elastic response spectrum and the convective spectral accelerations are obtained from a 0.5% damped elastic response spectrum.
The peak values of accelerated fluid pressure onto wall are located at the bottom of the container wall and presented in Figure 13 as functions of the tank slenderness parameter .The peak values of pressures of accelerated fluid at the bottom of the tanks wall are bigger in short and large-scale tanks compared with pressures in narrow and high (slender) tanks.
The fluid pressure is zero pressure at the original position of the free surface of the fluid in the static situation.The peak values of accelerated fluid pressure onto wall at the original position of the free surface are shown in Figure 14  Dependence of ℎ  , ℎ  , and ℎ  * on the tank slenderness parameter is seen in Figure 17; they correlate well with ℎ  , ℎ  , and ℎ  * in Figures 4, 5, and 7. Figure 18 documents comparison of the total base shear V in [kN] of the fluid contained in tank as function of the tank slenderness parameter . Figure 19 shows comparison of component ratio of impulsive and convective component of base shear V on tank as functions of the tank slenderness Water surface level parameter.If the tank slenderness is smaller, then the component of convective base shear is bigger and the component of impulsive base shear is smaller.Comparison of component ratio of impulsive and convective fluid moments in tank as functions of the tank slenderness parameter  can be seen in Figures 21 and 22. Figure 21 shows comparison of total the bending moments M and Figure 22

Conclusions
The seismic analysis of liquid storage containers is highly complex problem.Therefore, the knowledge of earthquake fluid effect on solid domain of tanks and appropriate modeling of the liquid storage tanks is essential for dynamic analysis and seismic response evaluation of earthquake resistance structure/facility, tanks.The dynamic behavior for liquidcontaining structures has to use different approach as for solid structure.The weight of storage liquid varies in time because of variable liquid storage level; accordingly there exists sloshing of fluid which influences the response of the tank.
The thoroughgoing dynamic analysis of contained liquid is necessary.The knowledge of forces acting onto containers, pressures in solid of tanks, surface sloshing process, and maximal height of liquids wave during an earthquake plays essential role in reliable and durable design of earthquake resistance containers.The analysis of a coupled multiphysics system is frequently required nowadays, to understand the behavior of the system.The motion of contained fluid in tanks may by expressed as the sum of two separate contributions, impulsive and convective actions.The mechanical model represents the dynamic analysis of fluid-filled tank, using the concept of generalized SDOF systems, representing the impulsive and convective (sloshing) modes of vibration.This paper summarizes the results of a comprehensive analytical investigation concerning the seismic analysis of groundsupported cylindrical tank-fluid filling effects on the tanks, subject to a horizontal component of earthquake ground motion.The tanks were excited by acceleration of Loma Prieta earthquake in California.Basic responses of the interest were the peak values of fluid pressures, the total base shear, the total bending, and the overturning moments of fluid as functions of the tank slenderness parameter and their ratios of impulsive and convective effects.In this study fluid filling of ground-supported cylindrical rigid tank is analyzed.(iii) the convective components of hydrodynamic pressures receive the maximum values on the tank wall at the original position of the fluid free surface and they are bigger in short and large-scale tanks than in narrow high (slender) tanks; (iv) the total mass of fluid, the total base shear, and the total moments are bigger in narrow slender tanks; these results are due to bigger volume of contained liquid; the convective fluid base shear is bigger in short and large-scale tanks; (viii) the bigger values of fluid impulsive bending and overturning moments are given in narrow and slender tanks; (ix) the convective components of moments (bending and overturning) give the bigger values in short and largescale tanks.
The knowledge of fluid effect on solid domain of tanks and their interdependence is, due to earthquake, significant for safe and economic design of earthquake-resistant structure/facility, containers.

Figure 2 :Figure 3 :
Figure 2: Variation of vertical distribution function of the impulsive pressure along the wall height for seven values of the tank slenderness parameter  = /, normalized to   .

Figure 4 :
Figure 4: Ratios ℎ  /H and ℎ  * /H as functions of the parameter tank slenderness.

Figure 5 :
Figure 5: Variation of convective pressures along the height for values of the tank slenderness parameter  = /: 0.5, 1.0, and 3.0.
as function of the tank slenderness parameter.Figures 15 and 16 show (i) comparison of the total mass m in [kg] of the fluid contained in tank as function of the tank slenderness parameter , (ii) comparison of component ratio of impulsive and convective fluid masses contained in tank as functions of the tank slenderness parameter .If the tank slenderness parameter  is smaller, then the component of convective fluid mass is bigger and the component of impulsive fluid mass is smaller.

Figure 13 :Figure 14 :
Figure 13: The peak value of fluid pressure in [kPa] on the bottom of tank wall as functions of the tank slenderness .

Figure 20
Figure 20 shows comparison of the total bending moments M and overturning moments  * in [kNm] of the fluid contained in tank as functions of the tank slenderness parameter .Comparison of component ratio of impulsive and convective fluid moments in tank as functions of the tank slenderness parameter  can be seen in Figures21 and 22. Figure21shows comparison of total the bending moments M and Figure22describes comparison of total overturning moments  * .If the tank slenderness is smaller, then the components of convective fluid moments are bigger and the components of impulsive fluid moments are smaller.

Figure 15 :Figure 16 :
Figure 20 shows comparison of the total bending moments M and overturning moments  * in [kNm] of the fluid contained in tank as functions of the tank slenderness parameter .Comparison of component ratio of impulsive and convective fluid moments in tank as functions of the tank slenderness parameter  can be seen in Figures21 and 22. Figure21shows comparison of total the bending moments M and Figure22describes comparison of total overturning moments  * .If the tank slenderness is smaller, then the components of convective fluid moments are bigger and the components of impulsive fluid moments are smaller.

Figure 19 :Figure 20 :
Figure 19: Comparison of the base shear components (impulsive and convective) as functions of the tank slenderness parameter .

Figure 21 :Figure 22 :
Figure 21: Comparison of bending moment components (impulsive and convective) as functions of the tank slenderness parameter .