Modified Proposed Provisions For Aseismic Design of Liquid Storage Tanks: Part II - Commentary and Examples
Modified Proposed Provisions For Aseismic Design of Liquid Storage Tanks: Part II - Commentary and Examples
32-25
Jain and Medhekar3,4 had proposed provisions on aseismic design of liquid storage tanks. In Part I of this paper,
need for modification to these provisions was highlighted and a set of modified provisions is proposed. In this
part II of the paper, a detailed commentary is provided to explain the rationale of the modified provisions. Two
solved numerical examples are also included to illustrate the application of the modified provisions.
When a tank containing liquid with a free surface is subjected to horizontal earthquake ground motion, tank wall
and liquid are subjected to horizontal acceleration. The
liquid in the lower region of tank behaves like a mass
that is rigidly connected to tank wall. This mass is termed
as impulsive liquid mass, which accelerates along with
the wall and induces impulsive hydrodynamic pressure on
tank wall and similarly on base. Liquid mass in the upper
region of tank undergoes sloshing motion. This mass is
termed as convective liquid mass and it exerts convective
hydrodynamic pressure on tank wall and base. Thus, total
liquid mass gets divided into two parts, i.e., impulsive mass
and convective mass. In spring mass model of tank-liquid
system, these two liquid masses are to be suitably represented. A qualitative description of impulsive and convective hydrodynamic pressure distribution on tank wall and
base is given in Fig. 1.
Sometimes, vertical columns and shaft are present inside
the tank. These elements cause obstruction to sloshing
motion of liquid. In the presence of such obstructions,
impulsive and convective pressure distributions are likely to
change. At present, no study is available to quantify effect
of such obstructions on impulsive and convective pressures.
However, it is reasonable to expect that due to presence
of such obstructions, impulsive pressure will increase and
convective pressure will decrease.
MODIFIED PROVISIONS
Ground Supported Tank
Assistant Professor, Department of Applied Mechanics, Visvesvaraya National Institute of Technology, Nagpur 440 011, India.
Professor, Department of Civil Engineering, Indian Institute of Technology, Kanpur 208 016, India.
297
TABLE 1
EXPRESSIONS FOR PARAMETERS OF SPRING
MASS MODEL
Circular Tank
tanh 0.866 Dh
mi
=
m
0.866 Dh
hi
= 0.375 for h/D 0.75
h
0.09375
= 0.5
h/D
mi
=
m
0.866 Lh
hi
= 0.375 for h/L 0.75
h
0.09375
= 0.5
h/L
0.866 Dh
hi
=
h
2tanh 0.866 Dh
0.125 for h/D 1.33
wall. On the other hand, hi is the height at which the resultant of impulsive pressure on wall and base is located from
the bottom of tank wall. Thus, if effect of base pressure
is not considered, impulsive mass of liquid, mi will act at
a height of hi and if the effect of base pressure is considered, mi will act at hi . Heights hi and hi , are schematically
described in Figs. 1a and 1b.
Similarly, hc is the height at which resultant of convective pressure on wall is located from the bottom of tank
wall; while, hc is the height at which resultant of convective pressure on wall and base is located. Heights hc and hc
are described in Figs. 1c and 1d.
Circular and rectangular tank
The parameters of spring mass model depend on tank
geometry and were originally proposed by Housner7 . The
parameters shown in Figs. 2 and 3 (of part I of the paper6 )
are slightly different from those given by Housner7 , and
have been taken from ACI 350.38 . Expressions for these
parameters are given in Table 1.
It may be mentioned that these parameters are for tanks
with rigid walls. In the literature, spring-mass models for
tanks with flexible walls are also available (Haroun and
Housner9 and Veletsos10 ). Generally, concrete tanks are
considered as tanks with rigid wall; while steel tanks are
considered as tanks with flexible wall. Spring mass models
for tanks with flexible walls are more cumbersome to use.
Moreover, difference in the parameters (mi , mc , hi , hi , hc ,
hc and Kc ) obtained from rigid and flexible tank models
is not substantial (Jaiswal et al11 ). Hence in the present
code, parameters corresponding to tanks with rigid wall are
recommended for all types of tanks.
298
Rectangular Tank
tanh 0.866 Lh
0.866 Lh
hi
=
h
2tanh 0.866 Lh
0.125 for h/L 1.33
= 0.45 for h/L > 1.33
tanh 3.16 Lh
mc
= 0.264
h
m
L
cosh 3.16 Lh 1.0
hc
= 1
h
3.16 Lh sinh 3.16 Lh
cosh 3.16 Lh 2.01
hc
= 1
h
3.16 Lh sinh 3.16 Lh
mg
tanh2 3.16 Lh
Kc = 0.833
h
In case of tanks with variable wall thickness (particularly, steel tanks with step variation of thickness),
thickness of tank wall at one-third height from the base
should be used in the expression for impulsive mode time
period.
Expression for Ti given in this section is applicable to
only those circular tanks in which wall is rigidly attached
to base slab. In some concrete tanks, wall is not rigidly
attached to the base slab, and flexible pads are used between
the wall and the base slab (Figs. 5d to 5f, of part I paper6 ).
In such cases, flexibility of pads affects the impulsive mode
time period. Various types of flexible connections between
wall and base slab described in Fig. 5 (of part I paper6 )
are taken from ACI 350.38 , which provides more information on effect of flexible pads on impulsive mode time
period.
Ground supported rectangular tank
Eurocode 816 and Priestly et al14 also specify the same
expression for obtaining time period of rectangular tank. h
is the height of combined center of gravity of half impulsive
mass of liquid (mi /2), and mass of one wall (mw ).
For tanks without roof, deflection, d can be obtained by
assuming wall to be free at top and fixed at three edges
(Fig. 2a).
ACI 350.38 and NZS 310618 have suggested a simpler
approach for obtaining deflection, d for tanks without roof.
As per this approach, assuming that wall takes pressure
q by cantilever action, one can find the deflection, d by
considering wall strip of unit width and height, h, which is
subjected to concentrated load, P = qh (Figs. 2b and 2c).
Thus, for a tank with wall of uniform thickness, one can
obtain d as follows:
d=
Iw =
P (h)3
3EIw
(1)
1.0 t 3
12
(2)
299
Soil structure interaction has two effects: Firstly, it elongates the time period of impulsive mode and secondly
it increases the total damping of the system. Increase in
damping is mainly due to radial damping effect of soil
media. A simple but approximate approach to obtain the
time period of impulsive mode and damping of tank-soil
system is provided by Veletsos10 . This simple approach has
been used in Eurocode 816 and Priestley et al14 .
FIG. 2. DESCRIPTION OF DEFLECTION d, OF RECTANGULAR
TANK WALL
The expressions for convective mode time period of circular and rectangular tanks are taken from ACI 350.38 ,
which are based on work of Housner6 . The coefficients Cc
in the expressions for convective mode time period plotted
in Figs. 6 and 7 (of part I paper6 ) are given below:
(a) For ground supported circular tank:
2
Cc =
3.68 tanh (3.68h/D)
(b) For ground supported rectangular tank:
2
Cc =
3.16 tanh (3.16(h/L))
Convective mode time period expressions correspond to
tanks with rigid wall. It is well established that flexibility
of wall, elastic pads, and soil does not affect the convective
mode time period.
For rectangular tank, L is the inside length of tank parallel to the direction of loading, as described in Fig. 3.
300
Damping
Importance factor (I ), is meant to ensure a better seismic performance of important and critical tanks. Its value
depends on functional need, consequences of failure, and
post earthquake utility of the tank.
In this code, liquid containing tanks are put in three categories and importance factor to each category is assigned
(Table 1). Highest value of I = 1.75 is assigned to tanks
used for storing hazardous materials. Since release of these
materials can be harmful to human life, the highest value
of I is assigned to these tanks. For tanks used in water distribution systems, value of I is kept as 1.5, which is same
as value of I assigned to hospital, telephone exchange,
and fire station buildings in IS 1893 (Part 1): 200220 . Less
important tanks are assigned I = 1.0.
Response reduction factor (R), represents ratio of maximum seismic force on a structure during specified ground
motion if it were to remain elastic to the design seismic
force. Thus, actual seismic forces are reduced by a factor
R to obtain design forces. This reduction depends on overstrength, redundancy, and ductility of structure. Generally,
liquid containing tanks posses low overstrength, redundancy, and ductility as compared to buildings. In buildings,
non structural components substantially contribute to overstrength; in tanks, such non structural components are not
present. Buildings with frame type structures have high
redundancy; ground supported tanks and elevated tanks
with shaft type staging have comparatively low redundancy.
Moreover, due to presence of non structural elements like
masonry walls, energy absorbing capacity of buildings is
much higher than that of tanks. Based on these considerations, value of R for tanks needs to be lower than that for
buildings. All the international codes specify much lower
values of R for tanks than those for buildings. As an example, values of R used in IBC 200021 are shown in Table 2.
It is seen that for a building with special moment resisting
frame value of R is 8.0 whereas, for an elevated tank on
frame type staging (i.e., braced legs), value of R is 3.0. Further, it may also be noted that value of R for tanks varies
from 3.0 to 1.5.
TABLE 2
VALUES OF RESPONSE REDUCTION FACTOR
USED IN IBC 200022
Type of Structure
Building with special reinforced concrete moment resisting concrete frames
Building with intermediate reinforced concrete moment
resisting concrete frames
Building with ordinary reinforced concrete moment
resisting concrete frames
Building with special steel concentrically braced frames
Elevated tanks supported on braced/unbraced legs
Elevated tanks supported on single pedestal
Tanks supported on structural towers similar to buildings
Flat bottom ground supported anchored steel tanks
Flat bottom ground supported unanchored steel tanks
Reinforced or prestressed concrete tanks with anchored
flexible base
Reinforced or prestressed concrete tanks with reinforced
nonsliding base
Reinforced or prestressed concrete tanks with unanchored and unconstrained flexible base
R
8.0
5.0
3.0
8.0
3.0
2.0
3.0
3.0
2.5
3.0
2.0
1.5
FIG. 4. COMPARISON OF BASE SHEAR COEFFICIENT OBTAINED FROM IBC 200022 AND
(a) IS 1893 (PART 1): 2002, FOR A BUILDING WITH
SPECIAL MOMENT RESISTING FRAME. (FROM
JAISWAL et al21 )
(b) PRESENT CODE, FOR TANKS WITH HIGHEST AND
LOWEST VALUES OF R. (FROM JAISWAL et al21 )
301
Base Moment
Sa /g = 2.5
= 1.36/T
for T 0.55
= 1.67/T
for T 0.67
Damping Factor
Table 3 of IS 1893 (Part 1): 200220 gives values of multiplying factors for 0% and 2% damping, and value for 0.5%
damping is not given. One cannot linearly interpolate the
values of multiplying factors because acceleration spectrum values vary as a logarithmic function of damping28 .
In Eurocode 816 the value of multiplying factor is taken as
1.673 and as per ACI 350.38 and FEMA 36829 , this value
is 1.5.
Base Shear
Elevated Tank
This Clause gives shear at the base of staging. Base shear
at the bottom of tank wall can be obtained from Clause of
ground supported tank.
Total Base Shear
FIG. 5. VARIATION OF IMPULSIVE AND CONVECTIVE BENDING MOMENT COEFFICIENTS WITH HEIGHT
(FROM MALHOTRA31 )
Elevated tank
Structural mass ms , which includes mass of empty container and one-third mass of staging is considered to be
acting at the center of gravity of empty container. Base of
staging may be considered at the top of footing.
Total Moment
(ELx 0.3Ely )
(ELx + 0.3Ely );
(ELx 0.3Ely )
(0.3ELx + Ely );
(0.3ELx Ely )
(0.3ELx + Ely );
(0.3ELx Ely )
This clause is adapted from Priestley et al14 . Since hydrodynamic pressure varies slowly in the circumferential
303
direction, the design stresses can be obtained by considering pressure distribution to be uniform along the circumferential direction.
Linearised Pressure Distribution on Wall
(Ah )i mi
g
D/2
qc =
(Ah )c mc
g
D/2
qi
(6hi 2h)
h2
qc
bc = 2 (6hc 2h)
h
bi =
Vertical ground acceleration induces hydrodynamic pressure on wall in addition to that due to horizontal ground
acceleration. In circular tanks, this pressure is uniformly
distributed in the circumferential direction.
Hydrodynamic Pressure
Design vertical acceleration spectrum is taken as twothird of design horizontal acceleration spectrum, as per
clause 6.4.5 of IS 1893 (Part 1: 2002)20 .
To avoid complexities associated with the evaluation of
time period of vertical mode, time period of vertical mode
is assumed as 0.3 seconds for all types of tanks. However,
for ground supported circular tanks, expression for time
period of vertical mode of vibration (i.e., breathing mode)
can be obtained using expressions given in ACI 350.38 and
Eurocode 816 .
While considering the vertical acceleration, effect of
increase in weight density of tank and its content may also
be considered.
h
= Mtot gD/2
2
1
h
=
(Ah )i
D
Thus, when h/D exceeds the value indicated above, the
tank should be anchored to its foundation. The derivation
assumes that the entire liquid responds in the impulsive
mode. This approximation is reasonable for tanks with high
h/D ratios that are susceptible to overturning.
Piping
F = kh s Hs2
where kh is the dynamic coefficient of earth pressure; s is
the density of the soil; and Hs is the height of soil being
retained. This force acting at height 0.6 h above the base
should be used to increase or decrease the at-rest pressure
when wall deformations are small.
(b)
P -Delta Effect
(c)
(d)
EXAMPLES
(e)
305
TABLE 3
COMPARISON OF RESULTS OBTAINED BY MODIFIED
PROVISIONS AND IS 1893:19841
and hs = 16.3 m
Idealization of tank
Two-Mass
One-Mass
Considered
Neglected
Ratio
1. Lateral stiffness of
staging
2. Time period
Impulsive mode,
Tank empty (Ti)
Tank full (Ti )
Convective mode,
Tank full (Tc )
17,806 kN/m
55,089 kN/m
0.32
0.66 sec
0.86 sec
0.37 sec
0.57 sec
1.78
1.51
3.14 sec
0.11
0.084
0.075
0.056
1.47
1.50
0.040
212 kN
281 kN
144 kN
248 kN
1.47
1.13
5. Overturning
moment (M )
Tank empty
Tank full
4,053 kN-m
5,448 kN-m
2,762 kN-m
4,757 kN-m
1.47
1.15
Solution:
(a) Spring mass model of tank
The total liquid mass (m) is 1000 t; mass of the
tank wall (mw ) is 15.9 t; mass of roof (mt ) is 5.1 t;
h/D ratio = 0.74. From Fig. 2 (of part I paper6 )
mi /m = 0.703; mc /m = 0.309; hi / h = 0.375;
hi / h = 0.587; hc / h = 0.677; hc / h = 0.727. Thus,
mi = 703 t; mc = 309 t; hi = 3.32 m; hi = 5.19 m;
hc = 5.98 m; hc = 6.43 m.
(b) Time period
The impulsive mode period is, Ti = 0.13 sec and
convective mode period is, Tc = 3.64 sec.
(c) Design horizontal seismic coefficient
The response reduction factor, R = 2.5. Importance
factor, I = 1.5. Zone factor, Z = 0.36. Soil strata is
hard.
For impulsive mode (Ti = 0.13 sec, damping =
2%), Sa /g = 3.5
For convective mode (Tc = 3.64 sec, damping =
0.5%), Sa /g = 0.48.
Thus, (Ah )i = 0.378 and (Ah )c = 0.052
(d) Base shear
Vi = 2, 685 kN; Vc = 158 kN
V = (2685)2 + (158)2 = 2, 690 kN
The design lateral base shear acting at the tank center
of gravity is V = 2, 690 kN.
sec h 3.674
Qcb (x) = 1.125
D 3 D
D
where,
2
(Av ) =
3
Z I Sa
2R g
307
Av
(Ah )c
FIG. 8. EQUIVALENT LINEARISED CONVECTIVE PRESSURE
DISTRIBUTION
This work has been supported through a project entitled Review of Building Codes and Preparation of Commentary and Handbooks awarded to IIT Kanpur by the
Gujrat State Disaster Management Authority (GSDMA),
Gandhinagar through World Bank finances. The views and
opinions expressed therein are those of the authors and not
necessarily those of the GSDMA or the World Bank.
308
Vc
Vi
Z
ai , bi
ac , bc
d
dmax
g
h
h
hc
hi
hs
ht
hw
hc
hi
l
m
mb
mc
mi
ms
mt
mw
m
w
p
pcb
pcw
pib
piw
pv
pww
q
qi
qc
t
tb
x
REFERENCES
1.
IS 1893:1984, Indian standard criteria for earthquake resistant design of structures, Bureau of
Indian standards, New Delhi.
2.
Jain, S. K. and Sameer, S. U., A review of requirements in Indian codes for aseismic design of elevated
water tanks, Bridge and Struct. Engr., V. XXIII
No. 1, 1993, pp. 116.
3.
Jain, S. K. and Medhekar, M. S., Proposed provisions for aseismic design of liquid storage tanks:
Part I Codal provisions, J. of Struct. Engg., V. 20,
No. 3, 1993, pp. 119128.
4.
Jain, S. K. and Medhekar, M. S., Proposed provisions for aseismic design of liquid storage tanks:
Part II Commentary and examples, J. of Struct.
Engg., V. 20, No. 4, 1994, pp. 167175.
5.
6.
7.
Housner, G. W., Dynamic analysis of fluids in containers subjected to acceleration, Nuclear Reactors
and Earthquakes, Report No. TID 7024, U.S. Atomic
Energy Commission, Washington D.C. 1963a.
8.
ACI 350.3, 2001, Seismic design of liquid containing concrete structures, American Concrete Institute, Farmington Hill, MI, USA.
9.
10.
Veletsos, A. S., Seismic response and design of liquid storage tanks, Guidelines for the seismic design
309
12.
13.
14.
15.
16.
17.
18.
19.
Sameer, S. U. and Jain, S. K., Approximate methods for determination of time period of water tank
staging, The Indian Conc. J., V. 66, No. 12, 1992,
pp. 691698.
20.
21.
22.
310
IBC 2000, International Building Code International Code Council, Falls Church, Virginia, USA,
2000.
Jaiswal, O. R., Rai, D. C. and Jain, S. K., Codal
provisions on design seismic forces for liquid storage
tanks: a review, Report No. IITK-GSDMA-EQ-01V1.0, Indian Institute of Technology Kanpur, Kanpur,
2004a.
23.
24.
Rai, D. C., Retrofitting of shaft type staging for elevated tanks, earthquake spectra, ERI, V. 18 No. 4,
2002, pp. 745760.
25.
AWWA D-100, Welded steel tanks for water storage, American Water works association, Colorado,
USA, 1996.
26.
27.
28.
Newmark, N. M. and Hall, W. J., Earthquake spectra and design, Engineering monograph, Earthquake Engineering Research Institute, Berkeley,
USA, 1982.
29.
FEMA 368, 2000, NEHRP recommended provisions for seismic regulations for new buildings and
other structures, Building Seismic Safety Council,
National Institute of Building Sciences, USA.
30.
Malhotra, P. K. Wenk, T. and Wieland, M., Simple procedure for seismic analysis of liquid-storage
tanks, Structural Engineering International, 2000,
pp. 197201.
31.
32.
Munshi, J. A. and Sherman, W. C., Reinforced concrete tanks, Concrete International, Feburary, 2004,
pp. 101108.
33.
Mononobe, N. and Matsuo, H., On the determination of earth pressure during earthquakes, Proc. of
World Engg. Congress, 1929.
34.
35.
36.
37.