DESIGN RECOMMENDATION FOR
STORAGE TANKS AND THEIR SUPPORTS
WITH EMPHASIS ON SEISMIC DESIGN(2010 EDITION)

ARCHITECTURAL INSTITUTE OF JAPANForeword
The Architectural Institute of Japan set up a “Sub-Committee for Design of Storage Tanks” toprovide a recommendation of the seismic design for the various types of storage tanks incommon use throughout Japan. The Sub-Committee first published “Design Recommendationfor Storage Tanks and Their Supports” in 1984, and amended it in the 1990, 1996 and 2010
publications. This current revised recommendation provides bulk material pressures for silos
and further guidance on seismic design methods for storage tanks based on the horizontal loadcarrying capacity in the structure.
It is envisaged by publishing the English version of “Design Recommendation for Storage
Tanks and Their Supports” that the above unique design recommendation will be promoted tothe overseas countries who are concerned on the design of storage tanks and the activities of
the Architectural Institute of Japan will be introduced them too. It is addressed how thisrecommendation is in an advanced standard in terms of the theory of the restoring force
characteristics of the structure considering the Elephant Foot Bulge (EFB), the effect of the
uplifting tank and the plastic deformation of the bottom plate at the shell-to-bottom juncture in
the event of earthquake, the design spectrum for sloshing in tanks, the design pressure for silos,
and the design methods for the under-ground storage tanks as well. The body of therecommendation was completely translated into English but the translation of the commentary
was limited to the minimum necessary parts for understanding and utilising the theories andequations in the body. The sections 2.7 Timbers and 4.9 Wooden Storage Tanks are omitted asthey are not common in overseas countries. Listing of some Japanese references are omitted
and some new references are added in this English version.
It is worthy of note that the contents of this recommendation have been verified and theirefficacy confirmed by reference to data obtained from the tanks and silos damaged during
recent earthquakes in Japan, including the Hanshin-Awaji Great Earthquake in 1995 and the
Niigataken Chuetu-oki Earthquake in 2007.

i

Introduction
There are several types of storage tanks, e.g., above-ground, flat-bottomed, cylindrical tanks forthe storage of refrigerated liquefied gases, petroleum, etc., steel or concrete silos for the storageof coke, coal, grains, etc., steel, aluminium, concrete or FRP tanks including elevated tanks forthe storage of water, spherical tanks (pressure vessels) for the storage of high pressure liquefied
gases, and under-ground tanks for the storage of water and oil. The trend in recent years is forlarger tanks, and as such the seismic design for these larger storage tanks has become more
important in terms of safety and the environmental impact on society as a whole.
The failure mode of the storage tank subjected to a seismic force varies in each structural type,with the structural characteristic coefficient (Ds) being derived from the relationship between
the failure mode and the seismic energy transferred to, and accumulated in the structure. A
cylindrical steel tank is the most common form of storage tank and its normal failure mode is abuckling of the cylindrical shell, either in the so called Elephant Foot Bulge (EFB), or as
Diamond Pattern Buckling (DPB). The Ds value was originally calculated with reference toexperimental data obtained from cylindrical shell buckling, but was later re-assessed andmodified based on the restoring force characteristics of the structure after buckling. Those
phenomena at the Hanshin-Awaji Great Earthquake and the Niigataken Chuetu-oki Earthquake
were the live data to let us review the Ds value. For the EFB, which is the typical bucklingmode of a cylindrical shell storage tank for petroleum, liquefied hydrocarbon gases, etc., it
became possible to ascertain the buckling strength by experiments on a cylindrical shell byapplying an internal hydrodynamic pressure, an axial compressive force, and a shear force
simultaneously. Details of these experiments are given in Chapter 3.
The seismic design calculations for other types of storage tanks have been similarly reviewed
and amended to take into account data obtained from recent experience and experiments.
Design recommendation for sloshing phenomena in tanks has been added in this publication.
Design spectra for sloshing, spectra for long period range in other words, damping ratios forthe sloshing phenomena and pressures by the sloshing on the tank roof have been presented.
For above-ground vertical cylindrical storage tanks without any restraining element, such as
anchor bolts or straps, to prevent any overturning moment, only the bending resistance due tothe uplift of the rim of bottom plate exists. This recommendation shows how to evaluate theenergy absorption value given by plasticity of the uplifted bottom plate for unanchored tanks,
as well as the Ds value of an anchored cylindrical steel-wall tank.
As the number of smaller under-ground tanks used for the storage of water and fuel isincreasing in Japan, the Sub-committee has added them in the scope of the recommendationand provided a framework for the seismic design of under-ground tanks. The recommendation
has accordingly included a new response displacement method and a new earth pressure
calculation method, taking into account the design methods adopted by the civil engineering
fraternity.
For silo design, additional local pressure which depends on eccentricity of discharge outlet, andequations which give approximate stress produced by this pressure are given in this 2010 publication.
August 2011
Sub-committee for Design of Storage Tanks
The Architectural Institute of Japanii

COMMITTEE FOR STRUCTURES (2010)
Chairman
Secretaries
Members

Masayoshi NAKASHIMA
Hiroshi OHMORI
Hiroshi KURAMOTO(Omitted)

Kenji MIURA

SUB-COMMITTEE FOR DESIG OF STORAGE TAKS (2010)
Chairman
Secretaries
Members

Yukio NAITO
Nobuyuki KOBAYASHI
Hiroshi AKIYAMA
Hitoshi KUWAMURA
Hideo NISHIGUCHI
Yutaka YAMANAKA

Hitoshi HIROSE
Yasushi UEMATSU
Minoru KOYAMA
Katsu-ichiro HIJIKATA
Jun YOSHIDA

Toshio OKOSHI
Kouichi SHIBATA
Hiroaki MORIWORKIG GROUP FOR SEISMIC LOAD AD RESPOSE
ESTIMATIO OF STORAGE TAKS (2008)
Chairman
Secretary
Members

Yukio NAITO
Nobuyuki KOBAYASHITokio ISHIKAWA
Nobuo NAKAI
Katsu-ichiro HIJIKATA

Hitoshi HIROSE
Shinsaku ZAMA
Hideo NISHIGUCHI
Tetsuya MATSUIiii

Hitoshi SEYA

Table of Contents

1

General --------------------------------------------------------------------------------- 1
1.1 Scope ----------------------------------------------------------------------------------- 1
1.2 Notations ------------------------------------------------------------------------------- 1

2

Materials ------------------------------------------------------------------------------- 3
2.1 Scope ----------------------------------------------------------------------------------- 3
2.2 Steels ----------------------------------------------------------------------------------- 3
2.3 Stainless Steels ------------------------------------------------------------------------ 4
2.4 Concrete -------------------------------------------------------------------------------- 4
2.5 Aluminium Alloys -------------------------------------------------------------------- 4
2.6 FRP ------------------------------------------------------------------------------------- 5
2.7 Timber ---------------------------------------------------------------------------------- 6

3

Design Loads -------------------------------------------------------------------------- 9
3.1 General --------------------------------------------------------------------------------- 9
3.2 Dead Loads ---------------------------------------------------------------------------- 9
3.3 Live Loads ----------------------------------------------------------------------------- 9
3.4 Snow Loads---------------------------------------------------------------------------- 9
3.5 Wind Loads -------------------------------------------------------------------------- 10
3.6 Seismic Loads ----------------------------------------------------------------------- 10
3.7 Design of Cylindrical Metal Shell Wall Against Buckling ------------------- 23
3.8 Earth Pressures ---------------------------------------------------------------------- 41

4

Water Tanks ------------------------------------------------------------------------- 44
4.1 Scope --------------------------------------------------------------------------------- 44
4.2 Structural Design ------------------------------------------------------------------- 44
4.3 Reinforced Concrete Water Tanks ----------------------------------------------- 46
4.4 Pre-stressed Concrete Water Tanks ---------------------------------------------- 47
4.5 Steel Water Tanks ------------------------------------------------------------------ 47
4.6 Stainless Steel Water Tanks ------------------------------------------------------- 47
4.7 Aluminium Alloy Water Tanks --------------------------------------------------- 48
4.8 FRP Water Tanks ------------------------------------------------------------------- 48
4.9 Wooden Water Tanks -------------------------------------------------------------- 49

5

Silos----------------------------------------------------------------------------------- 51
5.1 Scope --------------------------------------------------------------------------------- 51
5.2 Structural Design ------------------------------------------------------------------- 52
5.3 Reinforced Concrete Silos --------------------------------------------------------- 60
5.4 Steel Silos ---------------------------------------------------------------------------- 60

6

Seismic Design of Supporting Structures for Spherical Tanks -------------------- 62
6.1 Scope --------------------------------------------------------------------------------- 62
6.2 Physical Properties and Allowable Stresses ------------------------------------- 63
6.3 Seismic Design ---------------------------------------------------------------------- 64

iv

7

Seismic Design of Above-ground, Vertical, Cylindrical Storage Tanks ---------- 76
7.1 Scope --------------------------------------------------------------------------------- 76
7.2 Seismic Loads ----------------------------------------------------------------------- 80
7.3 Evaluation of Seismic Capacity--------------------------------------------------- 95

8

Seismic Design of Under-ground Storage Tanks ---------------------------------- 100
8.1 Scope -------------------------------------------------------------------------------- 100
8.2 Seismic Design --------------------------------------------------------------------- 100
8.3 Evaluation of Seismic Capacity-------------------------------------------------- 104

Appendix A
A1
A2
A3
A4
Appendix B

Design and Calculation Examples ----------------------------------- 106
Steel Water Tower Tanks --------------------------------------------- 106
Steel Silo ---------------------------------------------------------------- 120
Supporting Structures for Spherical Tanks ----------------------------- 149
Above-ground, Vertical, Cylindrical Storage Tanks ------------------ 154
Assessment of Seismic Designs for Under-ground
Storage Tanks ----------------------------------------------------------- 160v

1.

General

1.1

ScopeThis Design Recommendation is applied to the structural design of water storage
tanks, silos, spherical storage tanks (pressure vessels), flat-bottomed, cylindricalabove-ground storage tanks and under-ground storage tanks.
-------------------------------------------------------------------------------------------------Commentary:
This Design Recommendation is applied to the structural design, mainly the seismic design, ofwater storage tanks, silos, spherical storage tanks (pressure vessels), flat-bottomed, cylindrical,
above-ground storage tanks and under-ground storage tanks. As common requirements chapter 2
calls for material specifications which are applicable to the above -mentioned tanks, and chapter 3
calls for loads and buckling designs. Chapter 3 especially highlights the modified seismic
coefficient method and the modal analysis as approved evaluation methods of seismic design basedon the reference design acceleration response spectrum and the structural characteristic coefficient
Ds. The seismic design method adopted in the Recommendation is the allowable stress method
modified with B, the ratio of the horizontal load-carrying capacity in the structure to the short-termallowable yield strength. The design for each tank is described in chapter 4 and onwards. Chapters
4 and 5 call for general requirements of structural design of tanks and their supporting structures
for water storage tanks and silos, respectively. Chapters 6 calls for requirements of seismic design
only for supporting structures of spherical storage tanks. Chapters 7 and 8 calls for requirementsof flat-bottomed, cylindrical above-ground storage tanks and under-ground storage tanks,
respectively.
Examples of design procedure for each type of tanks in chapters 4 through 8 are included inappendices.

1.2

otations
The following notations are applicable as common notations through the chapters inthis Recommendation, and each chapter includes some additional notations to be
specifically used in the chapter.
B
C
Cj
Ds
Dh
Dη
E
F
b fcr
c fcr
s fcr
g

ｈ
I
j
K

Chapter2

ratio of the horizontal load-carrying capacity in the structure to the short-termallowable yield strength
design yield shear force coefficientdesign yield shear force coefficient at the j-th natural mode
structural characteristic coefficient Dh × Dη
coefficient determined by the damping of the structure
coefficient determined by the ductility of the structure
Young’s modulus (N/mm2)
yield strength (N/ mm2)
allowable bending stress in the cylindrical wall (N/mm2)
allowable compressive stress in the cylindrical wall (N/mm2)
allowable shear stress in the cylindrical wall (N/mm2)
acceleration of gravity (= 9.8m/s2)
damping ratio
importance factor
order of natural frequency
lateral design seismic coefficient

1

KH
KV
KHS
KHB

lateral design seismic coefficient in the groundvertical design seismic coefficient in the ground
lateral design seismic coefficient at the ground surface
lateral design seismic coefficient at the bedrock surfaceｋthe maximum number of vibration modes which largely influence to seismic
responses
ｎtotal number of mass points
Qd
design yield shear force at the base of the structure (N)Qdi design shear force imposed just below the i-th mass in case that the structure isassumed the n-th lumped mass vibration system (N)RG
ground amplification factorrradius of cylindrical tank (mm)Sv
design velocity response spectrum (m/s)
Sa
design acceleration response spectrum (m/s2)
Sa1 acceleration response at the first natural period (m/s2)
Saj
acceleration response at the j-th natural period (m/s2)
Sv1 design velocity response at the first natural period in the sloshing mode (m/s)
Tthe first natural period of the structure (s)Tjj-th natural period of the structure (s)critical period determined by the ground classification (s)TG
UH(z) horizontal displacement at the depth z from the ground surface (m)UV(z) vertical displacement at the depth z from the ground surface (m)uijj-th natural mode at the i-th massW
design weight imposed on the base of the structure (N)weight of the i-th mass (N)Wi
Zs
seismic zoning factor
σb bending stress in the cylindrical wall (N/mm2)
σc average compressive stress in the cylindrical wall (N/mm2)
σh average hoop tensile stress in the cylindrical wall (N/mm2)
τ
shear stress in the cylindrical wall (N/mm2)

Chapter222.

Materials

2.1

Scope
This chapter calls for the recommended materials used in the tanks and their supports.
Allowable stresses of the materials are determined in accordance with therecommended method set up in each chapter for each type of storage tank.

2.2

SteelsSteel materials listed in Table 2.1 should apply for.
Table 2.1 Specifications for Steels
Specification
JIS G 3136
JIS G 3101
JIS G 3106
JIS G 3444
JIS G 3466
JIS G 5101
JIS G 3201
JIS G 3115
JIS G 3120

JIS G 3126
JIS B 1186
JIS Z 3211
JIS Z 3212
JIS Z 3351
JIS Z 3352
JIS G 3112
JIS G 3117
JIS G 3551
JIS G 3536
JIS G 3538
JIS G 3109

Chapter2

Type of steel Grades
Rolled Steels for Building Structure; SN400A/B/C, SN490B/C
Rolled Steels for General Structure; SS400, SS490, SS540
Rolled Steels for Welded Structure; SM400A/B/C, SM490A/B/C,
SM490YA/YB, SM520B/C, SM570
Carbon Steel Tubes for General Structural Purposes; STK400,
STK490Carbon Steel Square Pipes for General Structural Purposes;
STKR400, STKR490
Carbon Steel Castings
Carbon Steel Forgings for General Use
Steel Plates for Pressure Vessels for Intermediate Temperature
Service; SPV235, SPV315, SPV355, SPV450, SPV490
Manganese-Molybdenum and Manganese-Molybdenum-Nickel
Alloy Steel Plates Quenched and Tempered for Pressure Vessels;
SQV1A/1B, SQV 2A/2B, SQV 3A/3B
Carbon Steel Plates for Pressure Vessels for Low Temperature
Service; SLA235A/B, SLA325A/B, SLA365
Sets of High Strength Hexagon Bolt, Hexagon Nut and Plain
Washers for Friction Grip Joints
Covered Electrodes for Mild SteelCovered Electrodes for High Tensile Strength Steel
Submerged Arc Welding Solid Wires for Carbon Steel and Low
Alloy Steel
Submerged Arc Welding Fluxes for Carbon Steel and Low Alloy
Steel
Steel Bars for Concrete Reinforcement
Rerolled Steel Bars for Concrete Reinforcement
Welded Steel Wires and Bar Fabrics
Uncoated Stress-Relieved Steel Wires and Strands for
Prestressed Concrete
Hard Drawn Steel Wires for Prestressed Concrete
Steel Bars for Prestressed Concrete
JIS : Japanese Industrial Standards

3

2.3

Stainless Steels
Stainless steel materials listed in Table 2.2 should apply for.
Table 2.2 Specifications for Stainless Steels
Specification
JIS G 4304
JIS G 4305
JIS G 3601
JIS G 3459
JIS Z 3221
JIS Z 3321

Type of steel Grades
Hot Rolled Stainless Steel Plates, Sheets and Strips; SUS304,
SUS304L, SUS316, SUS316L, SUS444, SUS329 J 4 L
Cold Rolled Stainless Steel Plates, Sheets and Strips; SUS304,
SUS304L, SUS316, SUS316L, SUS444, SUS329 J 4 L
Stainless Clad Steels; SUS304, SUS304L, SUS316, SUS316L,
SUS444
Stainless Steel Pipes; SUS304TP, SUS304LTP, SUS316TP,
SUS316LTP, SUS444TP, SUS329 J 4 LTP
Stainless Steel Covered Electrodes
Stainless Steel, Wire Rods and Solid Wires

2.4

Concrete

2.4.1

Concrete materials should comply with JASS 5 (Japanese Architectural Standard
Specifications and Commentary for Reinforced Concrete Works (2003)) 5.4“Concrete Materials” except that par. 5N.13.3 “Pre-stressed Concrete” of JASS 5N“Reinforced Concrete Work at Nuclear Power Plants (2001)” should apply for prestressed concrete.

2.4.2

Quality of concrete materials should comply with Table 2.3 and below.
-

JASS 5.5 “Mixing Proportion”
JASS 5.6 “Production of Concrete”
JASS 5.7 “Transportation and Placement”
JASS 5.8 “Curing” etc.

JASS 5N.13.3 “Pre-stressed Concrete” should apply for pre-stressed concrete.
Table 2.3 Specifications for Concrete Materials

2.5

Specification

Minimum Value of Fc(N/mm2)

Normal Weight
Concrete

18

Aggregate
CoarseFine

Gravel or Crushed
Stone

Sand or Crushed
Sand

Aluminium Alloys
Aluminium alloys used for structural members should be selected from those

Chapter2

4designated in Table 2.5 which are in accordance with the specifications of JapaneseIndustrial Standards in Table 2.4.

Table 2.4 Aluminium and Aluminium Alloy Materials
Specification
JIS H 4000

Title
Aluminium and Aluminium Alloy Sheets and Plates, Strips andCoiled Sheets
Aluminium and Aluminium Alloy Rods, Bars and Wires
Aluminium and Aluminium Alloy Extruded Tubes and Colddrawn Tubes
Aluminium and Aluminium Alloy Welded Pipes and Tubes
Aluminium and Aluminium Alloy Extruded Shape
Aluminium and Aluminium Alloy Forgings
Aluminium Alloy Castings
Aluminium and Aluminium Alloy Welding Rods and Wires

JIS H 4040
JIS H 4080
JIS H 4090
JIS H 4100
JIS H 4140
JIS H 5202
JIS Z 3232

Table 2.5 Grades of Aluminium Alloy Materials
A 1060A 1070
A 1080A 1100
A 1200A 3003
A 3203A 5052A 5056
A 5083A 5086
A 5154A 5254A 5454A 5652A 6061-T4
A 6061-T6
A 6063-T5
A 6063-T6
A 7 N 01-T4
A 7 N 01-T6

A 7003-T5

2.6

FRP

2.6.1

Resins used for FRP (Fibre Reinforced Plastics) materials should be equal to or
better than the quality of the UP-G (Unsaturated Polyester Resin-Gneral) specified in
the Japanese Industrial Standard JIS K 6919 (Liquid unsaturated polyester resins forfibre reinforced plastics), and should have excellent weather and water resistantproperties.

2.6.2

Glass fibre used for FRP materials should be of the no alkali glass type in accordancewith the Japanese Industrial Standards indicated in Table 2.6. Any additional
processed glass fibres should also be manufactured using these fibres as raw
materials.
Table 2.6 Types of glass fibre for FRPSpecification
JIS R 3411
JIS R 3412
JIS R 3413
JIS R 3415
JIS R 3416
JIS R 3417

Chapter2

Type
Textile glass chopped strand mats
Textile glass rovings
Textile glass yarns
Textile glass tapes
Finished textile glass fabrics
Woven roving glass fabrics

52.6.3

FRP materials should be processed using the resins and glass fibres specified above,and by adding any necessary fillers, coloring agents, artificial fibres, etc., as required
to ensure that the performance requirements stipulated in Table 2.7 are satisfied.
Table 2.7 Physical Properties of FRP Materials
Property
Tensile strength
Flexural strength
Flexural modulus
Glass content
Barcol hardness
Water absorption

2.7

Performance> 59N/mm2> 78N/mm2> 5.9×103 N/mm2> 25wt%> 30< 1.0%

Testing Method
JIS K 7054
JIS K 7055
JIS K 7055
JIS K 7052
JIS K 7060
JIS K 6919

TimberThis section is omitted in this English version.

Commentary1.

Component Materials of FRP
FRP basically consists of a fibre reinforcement in a resin matrix. FRP used for water tanks, use acompound resin consisting of hardening agents, colouring agents, mineral fillers, etc., in an
unsaturated polyester resin. The fibre reinforcement is generally glass fibre of the no alkali glass typebundled from 50 ~ 2000 filaments of 8 ~ 13μm in diameter in a processed mat, cloth or roving state
depending upon the intended use. The resin, glass fibre, additional fillers, and other materials shouldbe equal to or better than the requirements specified in the relevant Japanese Industrial Standards
indicated in sections 2.6.1, 2.6.2 and 2.6.3 above. Additional consideration should be given to thehygiene requirements for materials used in drinking water tanks.

2.

Performance and Allowable Stress of FRP
The performance of FRP may be changed by varying the combinations of the resin, glass fibre, andother agents to meet the requirements indicated. For water tanks, the combinations should be suchthat the mechanical and physical performance criteria specified in 2.6.3 are satisfied.
An approximation of the relationship of tensile stress and strain of FRP is shown in Fig. 2.6.1. FRP
material with a tensile strength greater than 59A/mm2 is used provided the material characteristics
satisfy the minimum requirements noted in 2.6.3. The strain at break point is between 2% ~ 3%.
Allowable stress of FRP is determined for the static properties at normal temperature, taking into
account degradation within the durable period, material deviations, etc. Tables 2.6.1, 2.6.2, and2.6.3 show the current method used to determine the allowable stresses in the design of FRP watertanks. In this illustration, the static properties of FRP at normal temperature are determined as thestandard value, the limited value - considering the degradation within durable period (15years) - isthen calculated, and finally the allowable stress is calculated by dividing the limited value by a safety
factor (See Table 2.6.3.).

Chapter2

6Static properties at normal temperature

Ft

Limited value (15 years)

Stress

0.7 Ft0.318 Ft

Allowable stress - short term loading (0.7 Ft / S1)
Allowable stress - long term loading (0.7 Ft /1.5 S1)

0.212Ft

0
Fig. 2.6.1

Strain %Characteristics of tensile stress and strain of FRP and an example of allowablestress determination. In this case, S1=2.2 (Mat hand-lay-up material).

Table 2.6.1 Allowable stress based on fracture strength
Standard Value *1

Limited Value

Tensile strength

Ft

Flexural strength

Fracture Strength

Allowable Stress
Short Term*2

Long Term*3

0.7 Ft

0.7 Ft / S1*4

0.7 Ft / 1.5S1

Fb

0.6 Fb

0.6 Fb / S1

0.6 Fb / 1.5S1In-plane shear strengthFs

0.7 Fs

0.7 Fs / S1

0.7 Fs / 1.5S1Interlaminar shear strength

FIS

0.7 FIS0.7 FIS / S1

0.7 FIS / 1.5S1Transverse shear strength

FT

0.6 FT

0.6 FT / S1

0.6 FT / 1.5S1

Compressive strength

FB

0.7 FB

0.7 FB / S1

0.7 FB / 1.5S1

*1
*2
*3

*4

Standard value -static properties at normal temperature
Short-term - considering short-term loading
Long-term - considering long-term loading
The definition of S1 is given in Table 2.6.3.

For long-term loading, hydrostatic pressure, dead loads, etc., should be divided by 1.5 to take into
account creep phenomena.

Chapter2

7

Table 2.6.2 Allowable Properties of Modulus

Standard Value

Limited Value

Allowable
Properties

Tensile Modulus

Et

0.8 Et0.8 Et / S2*

Flexural Modulus

Eb

0.8 Eb

0.8 Eb / S2

In-Plane Shear Modulus

G

0.8 G

0.8 G / S2

Modulus* The definition of S2 is given in Table 2.6.3.
Table 2.6.3 Safety factor of FRP (S1, S2)
Strength Criteria of Material (S1)

Stiffness Criteria of Structure (S2)

1.72 LS

1.58 LS

LS is the coefficient of dispersion for material properties of FRP.
(In the case of fracture strength, tensile strength and flexural strength; in the case of stiffness,tensile modulus and flexural modulus.)
LS is determined by the following equation, based on material properties obtained from more than
10 test pieces.
LS = 1 /{ 1 − 3.09( σ / x )}where:xσ= average of material properties= standard deviation.

The allowable properties to be determined depends on the criteria of fracture strength or stiffness;i.e., if the criteria is based on strength of material, the allowable stress should be determined on thebase of strength, and if the criteria is based on stiffness of structure, the allowable properties
should be determined on the base of stiffness. The safety factor should be the value considering theeffect of circumstance temperature of water tanks used, and deviation of material properties.
The material properties of FRP should be determined by the following test methods.
JIS K 7051:
JIS K 7054:
JIS K 7055:
JIS K 7056:
JIS K 7057:
JIS K 7059:

Chapter2

General rules for testing methods of glass fibre reinforced plastics.
Testing method for tensile properties of glass fibre reinforced plastics.
Testing method for flexural properties of glass fibre reinforced plastics.
Testing method for compressive properties of glass fibre reinforced plastic.
Fibre-reinforced plastic composites - Determination of apparent interlaminarshear strength by short-beam method.
Testing method for in-plane shear properties of glass fibre reinforced plastics.

8

3

Design Loads

3.1

General3.1.1

Loads should be applied to the structural design of a tank according to its intended
use, size, structure type, materials, design lifetime, location and environment, inorder to assure life safety and to maintain its essential functions.

3.1.2

The applied loads should be as follows, and their combinations should be definedconsidering the actual probability of occurrence.

(i)
Dead loads(ii) Live loads(iii) Snow loads(iv) Wind loads(v) Seismic loads(vi) Impulse and suction due to content sloshing, and pressure due to content(vii) Thermal stresses(viii) Shock, e. g., by crane(ix) Fatigue loads(x) Soil and water pressures(xi) Others. e.g., load from mechanical device.
3.2

Dead Loads
Dead loads are the sum of the weights of the tank, its associated piping andequipment and other fixed appurtenances.

3.3

Live Loads

3.3.1

Live loads should be considered to be the contents of the tank, the temporary weightof personnel, and the weight of other temporary equipment not normally fixed to thetank. If such loads involve a impulsive force, they should be increased by a suitable
impulse factor.

3.3.2

Temporary live loads may be discounted if they have been taken into consideration
with other load combinations.

3.3.3

For seismic designs, the weight of the tank contents should be considered as being
divided into fixed weight (impulsive mass) and free weight (convective mass).

3.4

Snow Loads
Snow loads should be defined by considering the location, topography, environment,
density of the snow, snow accumulating period, and the shape and temperature of thetank, as defined in “Recommendations for Loads on Buildings” (AIJ- 2004).

Chapter393.5Wind loads
Wind loads should be defined by considering the shape of the tank, its structuralcharacteristics, the location, and environment, as defined in “Recommendations for
Loads on Buildings” (AIJ - 2004).

3.6

Seismic Loads

3.6.1

Seismic Loads for Above-ground Storage Tanks

3.6.1.1

Design for Impulsive Mass
Design seismic loads for above-ground storage tanks should be calculated by either
one of the following methods:
(i)
(ii)

Modified Seismic Coefficient Method, orModal Analysis

The Modified Seismic Coefficient Method should be used for the design seismic
loads of tank foundations.
3.6.1.2

Modified Seismic Coefficient Analysis
Design yield shear force, Qd, should be calculated using equations (3.1) and (3.2).
(3.1)

Qd = CW(3.2)

C = Z s IDs

S a1
gwhere:C ≥ 0.3ZsI
Notations:Qd
design yield shear force (N)C
design yield shear force coefficientWdesign weight imposed on the base of the structure, which isequal to the sum of dead weight of structure and weight of
impulsive mass of liquid content (N)Sa1
design acceleration response spectrum corresponding to thefirst natural period, given in par. 3.6.1.6 (m/s2)gacceleration of gravity = 9.8 (m/s2)
Zs
seismic zone factor, which is the value stipulated in the
Building Standards Law and the local governmentI
importance factor, given in par. 3.6.1.7
Dsstructural characteristic coefficient, which is the value
calculated in equation (3.3) and as further detailed in
Chapters 4, 5, 6, and 7.
(3.3)

Chapter3

Ds = DhDη

10where:Dhis the coefficient determined by the damping of the structureDη
is the coefficient determined by the ductility of the structure
When the structure is assumed to have the n mass vibration system, the design shear
force imposed just below the i-th mass, Qdi, should be calculated using equation (3.4).
n

∑W h

m

∑W h

m

m

(3.4)

Qdi = Qd mn=i

m

m =1

where:
is the vertical distribution of the design shear force imposed
Qdi
just below the i-th mass (N)Notations:ntotal number of masseshiheight at the i-th mass from the ground level (m)Wiweight of the i-th mass (N)
3.6.1.3

Modal Analysis
When the structure is considered as the vibration of an n mass-system, the design
yield shear force imposed just below the i-th mass, Qdi, should be calculated using
equation (3.5).
1

(3.5)

2
2
 k  n 
Qdi = ∑ C j ∑ Wm β j umj    j =1  m=i(3.6)

C j = Z s IDs

S ajgwhere:
Qd1 ≥ 0.3ZsIW
Notations:jorder of natural frequencykthe maximum number of vibration modes which largely
influences seismic responsesntotal number of masses
Wi
weight of the i-th mass (N)βjparticipation factor of the j-th natural modej-th natural mode at the i-th massuij
Cj
design yield shear force coefficient at the j-th natural mode
Saj
design acceleration response spectrum at the j-th natural
period, given in par. 3.6.1.6 (m/s2)

Chapter3113.6.1.4

Design Seismic Loads of Weight of Convective Mass
When the stored content is liquid, the design seismic loads of the weight of
convective mass should be obtained from the design velocity response spectrum, Sv1,
corresponding to the first natural period of the convective mass. See also chapter 7.

3.6.1.5

Design Seismic Loads of Tank Foundation
The sum of the lateral force, which is obtained from the dead weight multiplied by
the horizontal design seismic coefficient, K, as shown in (3.7), and the shear force
transmitted from the upper structure should be the design seismic load acting on thetank foundation.
K ≥ 0.15(3.7)
3.6.1.6

Design Response Spectra
Design velocity response spectrum, Svj, and design acceleration response spectrum,
Saj, should be obtained from the equation (3.8) or (3.9).
(3.8)

when Tj < TG

Svj = 1.56 Tj
Saj = 9.8

(3.9)

when Tj ≥ TG

Svj = 1.56 TG
Saj = 9.8 TG / Tj

Notations:Tjj-th natural period of the structure (s)critical period determined by the ground classification shown in Table 3.1
TG
below (s).
Table 3.1 Critical Period, TG (s)Classificationof Ground

Ground Conditions(1)

Type 1

Type 2(2)(3)

the ground before the Tertiary (hereinafter
called bedrock)diluviums
alluviums which are less than 10m inthickness to bedrock

alluviums which are less than 25 m in thickness
to bedrock and of which the soft layer is less than
5 m in thickness(1)

the ground other than the above(2)

the ground which is unknown in properties

Type 3

Chapter3

TG (s)

0.64

0.96

1.28

12

3.6.1.7

Importance Factor
An importance factor is determined corresponding to the classification of the seismic
design I, II, III, or IV given in Table 3.2 in consideration of significance ofearthquake damages and consequence to environment.
Table 3.2 Importance Factor, I
Classificationof Seismic
Design3.6.1.8

Description

I

I

Storage tanks which have small capacities and do
0.6 and greater
not contain hazardous materials

II

Storage tanks which have medium to large
capacities, do not contain hazardous materials, and
0.8 and greater
will not have any significant consequential effect inthe event of earthquake damage.

III

Storage tanks to which the significance ofearthquake damages is equivalent to common use 1.0 and greaterbuildings

IV

Storage tanks which contain hazardous materials andthe failure of which will constitute a secondary 1.2 and greater
disaster

Design Seismic Loads for Allowable Stress Design
The allowable design stress for storage tanks should be obtained from the design
story shear force, Qei, given in equation (3.10).
Q(3.10)
Qei = diB
Notations:
Qei
design story shear force for allowable stress design (N)B
ratio of horizontal load-carrying capacity in the structure toshort-term allowable yield strength, the value of whichshould be between 1.0 ~ 1.5 and should be determined foreach type of structure specified in the commentary to chapter
3 and the body of chapter 4 onwards.

3.6.1.9

Horizontal Load-carrying Capacity
Horizontal load-carrying capacity of storage tanks should be obtained from theallowable stresses given in par. 3.7.4.

Chapter313------------------------------------------------------------------------------------------------------------------Commentary:
Essentially seismic design as applied to storage tanks follows the requirements prescribed by the BuildingStandard Law where the importance of structures is not explicitly considered. However, as storage tanks are
more versatile than normal building structures in that they include a range of structures covering silocontaining forages and tanks for toxic materials, an importance factor, I, is introduced to accommodate thesedifferences. Whilst this chapter summarises the common features of seismic design for above-ground storage
tanks, specific design considerations - including the importance factor, is given in the chapters dedicated toeach principal tank design type.
1.

Seismic Input
The seismic energy input, ED , the principal cause for a structure’s elasto-plastic deformation, isapproximately expressed as shown in equation (3.6.1) as follows[3.13]:

ED =(3.6.1)MS v
2

2

Aotations:M
total mass of a structure
Svvelocity response spectrum
The acceleration response spectrum, Sa, prescribed by the Building Standard Law, is obtained fromequation (3.6.2):

S a = C0 Rt g

(3.6.2)Aotations:C0
standard design shear force coefficientRtvibration characteristic coefficientgacceleration of gravity
Sa and Sv are related to each other as shown in the following equation (3.6.3).

Sa =(3.6.3)2π
Sv
T

Aotations:Tnatural period for the 1st mode of a structure
Sa and Sv are indicated by broken lines in the curves shown in Figure 3.6.1, whilst the solid lines
represent the simplified curves expressed in equations (3.8) and (3.9).

10

9.80

5

Sv (m/s)

Sa (m/s2)

2.0
1.5
1.0
0.5

Fig.3.6.1 Design Response Spectra

Chapter3

14

2.

Seismic Design Loads
For the sake of determining seismic design loads, storage tanks can be simplified to a single mass
system, the restoring force characteristics of which are assumed to be of the elasto-plastic type withyield shear force Qy.
It is also assumed that the structure develops cumulative plastic deformation, δ 0 , in both positive
and negative directions under severe earthquake conditions as shown in Figure 3.6.2.

Fig.3.6.2 Damage and Cumulative Plastic Deformation
The relationship between the total energy input exerted by an earthquake and the damage of astructure is expressed by equation (3.6.4) as follows:(3.6.4)W p + We = E − Wh =Aotations:
E
Wp
We
Wh
h

E

(1 + 3h + 1.2 h )2total energy input
damage to a structure (cumulative plastic strain energy)
elastic vibrational energy
energy absorbed by damping
damping ratio

E - Wh can be approximately expressed by the following equation (3.6.5):(3.6.5)E − Wh =

MS v
2

2

Aotations:M
mass of a single-mass system
The spectra shown in Figure 3.6.1 are given for h = 0.05. E may then be expressed by equation(3.6.6):(3.6.6)E=

1.42 2 MS v
22

Wp and We are expressed by equations (3.6.7) and (3.6.8) respectively:(3.6.7)

Chapter3

W p = 2Q y δ 0

15

(3.6.8)We =Q yδ y2

Where the yield deformation, δy, and the yield shear force, Qy, are determined by equations (3.6.9)
and (3.6.10) respectively.
(3.6.9)(3.6.10)

δy =

Qy
k

Qy = MgC
Aotations:yield shear force
Qyyield deformation
δyspring constant of a singl-mass systemkCyield shear coefficient

The natural period of a single-mass system, T, is expressed by equation (3.6.11).
(3.6.11)M 
T = 2π  
 k 

1 2

C is obtained by substituting (3.6.7) and (3.6.8) into (3.6.4) and using (3.6.9) - (3.6.11) as shown inequation (3.6.12) as follows:(3.6.12)C=(3.6.13)

η=

2π S v
1 + 4η 1 + 3h + 1.2 h T g
1

1.42

δ0
δy

Aotations:η

averaged cumulative plastic deformation ratio as shown inequation (3.6.13).

Knowing that Sa = 2πSv /T , and applying Zs = I = 1.0, Ds in (3.3) is determined by comparing(3.2) and (3.6.12) using equations (3.6.14) and (3.6.15) as follows:(1)(3.6.14)Dh =(3.6.15)

Dη =

1.42
1 + 3h + 1.2 h
1
1 + 4η

Value of Dh
Large diameter storage tanks, such as oil storage tanks, are subject to an interaction between thesoil and the structure. Assuming that a structure in direct contact with the ground is a continuous
shear stratum, the one dimensional wave propagation theory yields an estimate of damping ratio,h, due to energy dissipation through the soil - structure interaction as illustrated in equation(3.6.16).
(3.6.16)h=

2 ρ 0V0π ρ 1V1

Aotations:

Chapter316ρ0 , ρ1densities of structure and the ground, respectively (kg/m3)

V0 , V1

shear wave velocities of structure and the ground,respectively (m/s)ρ 0 and V0 are expressed by equations. (3.6.17) and (3.6.18).
(3.6.17)

ρ0 =

M
HA0

(3.6.18)

V0 =

4H
T

Aotations:
M
H
A0
T

mass of structure (kg)height of structure (m)bottom area of structure (m2)fundamental period of structure (s)ρ 0V0 is approximately expressed by equation (3.6.19).
(3.6.19)

ρ 0V0 =

4M
A0T

Assuming ρ1 = 2.4 × 103 kg/m3 and V1 = 150 m/s for Type 3 Classification of Ground with somemargin, ρ1V1 is expressed by equation (3.6.20).
(3.6.20)

ρ1V1 =

360 3
⋅ 10
ah

Where ah is shown in Table 3.6.1.
Damping ratio, h, is then expressed by equation (3.6.21) by substituting equations (3.6.19) and(3.6.20) into equation (3.6.16).
(3.6.21)h=

7ah M⋅ 10− 6
A0T

If a large –scale cylindrical LPG storage tank on soft soil (Type 3) is given for an example, andassuming the diameter of the tank D=60m, liquid height H=25m, T=0.32s, A0=3,120m2 and
M=2.12 × 107kg, ρ 0V0 is estimated as ρ 0V0 = 8.5 × 10 4 kg /( s ⋅ m 2 ) . The damping ratio is thendetermined as, h=0.15, from equation (3.6.21).
These values are consistent with the damage noted after the Miyagiken-oki and other earthquakes,
and are confirmed by the analytical results of a three-dimensional thin-layer element model.
Equation (3.6.21) is applicable to flat tanks such as large-scale cylindrical oil storage tanks build
on soft ground. For non-flat tanks, the damping ratio is smaller.
The ρ1 and V1 for equation (3.6.16) can also be directly measured by velocity logging etc., ifnecessary.
Table.3.6.1 Value of α hClassification of GroundαhType 10.6Type 20.75

Type 3

1.0* The classification, Type 1~3, is defined in Table 3.1.

Chapter317(2) Value of Dη
As there are many variations in the design of storage tanks, estimating approximate values of η isaccordingly dependent on the information provided. It is therefore essential that sufficient details
are provided about the structural design to enable a reasonable estimation to be made. As atentative measure, the minimum value of η , and Dη as obtained from equation (3.6.15) should be
used (see Table 3.6.2).
Table.3.6.2 Value of Dη

Structural TypeηDηFrame Structure

1.3

0.40

1.0

0.45

RC

Wall Structure

S 0.75～0.25 0.5～0.7
RC : reinforced concrete
S : steel3.

Relationship between modified seismic coefficient method and modal analysis
A seismic load can be obtained by the modified seismic coefficient method or by modal analysis.
The advantage of the modified seismic coefficient method lies in its simplicity of application,
whereas modal analysis is a little bit more complicated, as indicated by equations (3.5) and (3.6),
and can result in a smaller seismic load than the modified seismic coefficient method.
When comparison of the design yield shear force at the base of structure is made using the modal
analysis and the modified seismic coefficient method for structures in which shear deformation isdominant, the ratio of the design yield shear force given by the modal analysis, Qd1, to the design
yield shear force given by the modified seismic coefficient method, Qd, lies between 1.0 and 0.8.
When flexural deformation is added, the ratio may be further decreased. Provided that structuralparameters can be adequately estimated, the lower seismic loads given by modal analysis may be
applied. However, in view of the uncertainty in estimating structural parameters, it is preferable
that Qd1 is limited to ensure that,(3.6.22)4

Qd1 ≥ 0.7 Qd

Design seismic load applied to liquid sloshing
The sloshing oscillation of a contained liquid is obtained from the velocity response spectrum.
Although the velocity response spectra in the range of longer periods are not yet fully clarified,equation(7.7) is considered to be a valid conservative estimate.

5.

Evaluation of seismic resistance
The seismic resistance of structures can be evaluated by the following alternative methods:(1) Horizontal load-carrying capacity design
Horizontal load-carrying capacity in each part of a structure is checked to ensure that the criteria
contained in equation (3.6.23) are satisfied.
(3.6.23)

Qyi ≥ Qdi
Aotations:
Qyi
Qdi

Chapter3

horizontal load-carrying capacity in the i-th mass
design yield shear force in the i-th mass18

(2) Allowable stress design
Similarly, the design shear force in each storey, Qei, is checked to confirm that the stresses causedby this force and the vertical loads are less than the allowable stresses.
Qei is obtained using equation (3.6.24).

(3.6.24)

Qei =

Qdi
B

Aotations:
Qei
design story shear force in the i-th mass pointB
the ratio of the horizontal load-carrying capacity to theshort-term allowable horizontal strength ( ≥ 1.0)
The first method corresponds to the ultimate strength design, whilst the second method corresponds
to the conventional allowable stress design. The second method is needed, because the first one is
not established yet for some types of storage tanks. Generally B is greater than unity, and whenapplied to ordinary building structures, the following values can be taken:
For frame structures
For wall structures

B
B

3.6.2

Seismic Load for Under-ground Storage Tanks

3.6.2.1

General==

1.2 ~1.5
1.0

The following paragraphs give details of under-ground seismic intensity and underground soil displacement as applicable to the seismic intensity and the response
displacement methods. Further details of their application can be found in Chapter 8
Seismic Design of under-ground Tanks, and Appendix B.

3.6.2.2

Under-ground seismic intensity
Horizontal seismic intensity, KH, is based on the linear relation to the depth, andshould be determined according to equation (3.11) below.
z(3.11)K H = K HS − ( K HS − K HB )H
where:
KHS
horizontal seismic intensity at the ground surface ascalculated by equation (3.12) below,
KHB
horizontal seismic intensity at the bedrock surface ascalculated by equation (3.13) below,zdepth from the ground surface (m),H
depth from the ground surface to the bedrock surface.
(3.12)(3.13)

Chapter3

K HS = Z S RG K 0
3
K HB = Z S K 0
4
where:
ZS
seismic zone factor,
RG
soil amplification factor given in Table 3.3 below,

19

K0

standard seismic intensity obtained from equation (3.14).

K 0 = 0.2
Vertical seismic intensity KV should be determined by equation (3.15).

(3.14)

KV =(3.15)

1
KH
2Table 3.3 Soil Amplification Factor RG
Ground Type

Soil Amplification Factor - RG

Type 1 Ground (diluvial deposit)

1.0Type 2 Ground (alluvium)

1.1Type 3 Ground (soft ground)

1.2

Note:

1) RG may be 0.9 for bedrock.
2) The classification, Type 1~3, is defined in Table 3.1.

3.6.2.3

Under-ground soil displacement
Horizontal and vertical displacements UH and UV of under-ground soil should bedetermined in accordance with equations (3.16) and (3.17).
(3.16)

UH =(3.17)UV =2

πSV K HB TS

2

H−z
H

1
UH
2

where:
horizontal displacement of under-ground soil at the depth z
from the ground surface (m),vertical displacement of under-ground soil at the depth z
from the ground surface (m),design velocity response spectrum for a unit seismic
intensity determined according to the surface soil natural
period TS (m/s),
natural period of the ground surface layer which isdetermined by equation (3.18), taking into account the shear
strain level of the surface soil (s)

UH
UV
SV
TS

(3.18)

TS = 1.25T0(3.19)

T0 =

n

4Hi

i =1

si

∑V

where:
Hi
Vsi

Chapter3

thickness of the stratum at the i-th layer (m),mean shear wave velocity at the i-th layer (m/s),

20

ntotal number of the layers.

However, when the bedrock surface level is close to the tank bottom level, UH may
be determined by equation (3.20)
UH =(3.20)

2

SV K HB TS cos(

πz

)

2H
π
------------------------------------------------------------------------------------------------------------------2

Commentary:1.

Under-ground seismic intensity
Under-ground seismic intensity may be calculated as a linear function of the depth using theseismic intensities at the ground and the bedrock surfaces provided the surface layers around thetank are nearly uniform. The term “bedrock” here is defined as 25 or more A-value for clay, 50 ormore A-value for sand, or 300m/sec. or more shear wave velocity .

2.

Under-ground soil displacement
Because under-ground tanks are bound by the surrounding soil, the seismic response is influenced
more by the deformation of the soil than that of their structural characteristics. Therefore,
regarding the soil deformation as forced displacement, as well as inertia force of the structure, isconsistent with actual condition.
The recommendation published in 1990 adopted a sine curve distribution for under-ground soil
displacement, where a dynamic shear deformation equation for homogeneous one-layer soil wasused with a sine curve assumption for shear wave form in the soil. Later FEM simulations
confirmed that assumed sine curve distribution gives a smaller relative displacement (deformation)
than the FEM results to tanks in deep bedrock cases, because a relatively soft soil rigidity gives
larger relative displacement in the surface region. Accordingly, it is recommended that linear
distribution of soil displacement for deep bedrock cases and sine curve distribution for shallow
bedrock cases are used in the seismic design for under-ground tanks.
Appendix B examines the seismic design of Under-ground storage tanks and gives the basis fordeep bedrock cases, whereas the basis for shallow bedrock cases is given in the following
paragraphs.
Although soil displacements during earthquakes are very complicated with the actual conditions
depending on the propagating earthquake, for design purposes it is assumed that all movements
are in shear wave. Sine wave motion and cosine ( πz / 2 H ) distribution in depth are assumed asin Figure 3.6.3.

thickness, H, of thesurface stratum

ground surface

bedrock

Figure 3.6.3 Under-ground soil displacement

Chapter3

21

If a homogeneous one-layer soil is assumed, the vibration equation for shear waves is as given inequation (3.6.25).

γt(3.6.25)

∂ 2u
∂u ∂
∂u
∂ 2u B
+
C
−
(
G
)
=
−
γ
t
∂t ∂z
∂z
∂t 2
∂t 2

where:

γtunit soil mass
horizontal displacement at depth z from the ground surface
horizontal displacement of the bedrock
damping factor
soil shear modulus

u
uB
C
G

Bedrock and the ground surface are assumed to be fixed and free respectively. Spectrum modal
analysis was performed for the first eigenvalue, assuming C = 0, then the maximum horizontal
displacement at ground surface umax is given in equation (3.6.26).

u max =

(3.6.26)

where:
umax

4

π

Sd =

2

π2

Sυ Ts

Sd

maximum horizontal displacement at the ground surface by
shear wave
the first natural period at the ground surface
displacement spectrum at the bedrockSυ

velocity spectrum at the bedrock ( = 2π S d / Ts )Ts

Equation (3.6.26) is the basis for the displacement from equation (3.20) in chapter 3.6.2.3.
Ts is given by the equation Ts=1.25T0 - see equations. (3.18) and (3.19), where a reduction in therigidity according to the strain level during earthquake is considered.

Sυ is evaluated as design velocity spectrum Sυ at bedrock surface multiplied by horizontal
seismic intensity KHB at bedrock surface. Sυ is calculated for 1g seismic intensity at bedrocksurface (see Fig. 3.6.4).

Sυ = Sυ K HB

(3.6.27)

Fig. 3.6.4 is obtained as an envelope of sixteen velocity response spectra, with soil damping ratio
0.2, of strong motion records at the second class (Type 2) grounds (alluvium) normalized to 1g
level.

1.0(m/s)

response velocity for1g seismic intensity(0.6, 1.5)

0.5
(0.1, 0.25)
0.1

natural period of ground surface layer (s)Figure 3.6.4

Chapter3

Velocity spectrum for 1g level

22

3.7

Design of Cylindrical Metal Shell Wall Against Buckling

3.7.1

Design Stress Criterion
The allowable stress design for both the elastic design and the ultimate strength
design is applied. The general design format for the design criterion is given by thefollowing two conditions, eqs. (3.21) and (3.22).
σcσ(3.21)+ b ≦1cf cr

τ

(3.22)
s

f cr

b

f cr

≦1

where:
σcaverage axial compressive stress = W/A(N/mm2)
σbmaximum axial compressive stress caused by the overturning moment=M/Z(N/mm2)
maximum shear stress caused by the lateral shear force =τ
2Q/A(N/mm2)W
axial compressive force(N)Moverturning moment(Nmm)Q
lateral shear force(N)A
sectional area(mm2)Z
elastic section modulus(mm3)allowable axial compressive stress(N/mm2)c f crallowable bending stress(N/mm2)
b f crallowable shear stress(N/mm2)s f cr
The allowable stresses defined above are composed of ones against the long-termload, ones against the short-term load exclusive of the seismic load, and ones against
the seismic load exclusively.
3.7.2

Allowable Stresses Against the Long-Term Load
The allowable stresses against the long-term axial compressive force, overturning
moment and lateral shear force are given by the following three terms (1), (2) and (3),respectively, which all depend on both σ h / F , the ratio of the average
circumferential tensile stress caused by the internal pressure σ h to the basic value fordetermining the yield stress F, and r/t, the ratio of the cylinder’s inside radius ofcurvature r to the cylinder’s wall thickness t.

(1)

Allowable axial compressive stress c f cr against the long-term vertical force(3.23)

for 0≦

σh(3.24)

for

σh

≧0.3

Chapter3F
where:

≦0.3

F
c

c

f cr = c f cr +f cr = f cr 0 ( 1 −23

0.7 f cr 0 − c f cr σ h
(
)
0.3
F

σh
F

)

σh

average circumferential tensile stress caused by the internalF
f cr

pressure (N/mm2)
basic value for determining the yield stress (N/mm2)
allowable axial compressive stress against the long-term load

f cr 0exclusive of internal pressure (N/mm2)
basic value for determining the strength of elephant-footcbulge (N/mm2)
The value of f cr 0 is given by the following equations according to the value of r/t.
rFE
(3.25)for ≦0.069
f cr 0 =
t
F
1.5
r
E
E
(3.26)
for 0.069 ≦ ≦0.807
FtF(3.27)r F 0.807 − t E 
f cr 0 = 0.267 F + 0.4 F  0.738 


r
1
E
for 0.807 ≦f cr 0 =σ cr 0 ,eFt
2.25

where:rt
EF
σ cr 0 ,ecylinder’s inside radius of curvature (mm)cylinder’s wall thickness (mm)Young’s modulus (N/mm2)
basic value for determining the yield stresses (N/mm2)
axially symmetric elastic buckling stress (N/mm2)

The value of σ cr 0, e is given by the following equation.
(3.28)

σ cr 0 ,e =where:ν

0.8
3( 1 − ν )
2

E

t
r

Poisson’s ratio

The value of c f cr is given by the following equations according to the value of r/t.
rE
≦ 0.377  tF

(3.29)

for(3.30)

Efor 0.377 
F
 

0.72

0.72

c

f cr =

F
1.5

r
E≦ ≦ 2.567 
t
F
 

0.72

0.72

rF 
 2.567 −   t E c f cr = 0.267 F + 0.4 F


2.190





Chapter324

(3.31)Efor 2.567 
F

0.72

≦

r
t

c

f cr =

1c σ cr ,e
2.25

where:
the elastic axially asymmetric buckling stress, based onc σ cr ,elower bound formula proposed by NASA[3.1].
The value of c σ cr , e is given by the following equation.
(3.32)
(2)

c

t 
r






1/ 2

1 r
 
 16  t 

σ cr ,e = 0.6 E 1 − 0.9011 − exp −



 
 
 

Allowable bending stress b f cr against the long –term overturning moment(3.33)

for 0≦

σh(3.34)

for

σh

≧0.3

≦0.3

F

f cr = b f cr +

b

F




b f cr = f cr 0 1 −0.7 f cr 0 − b f cr  σ h 
 
0.3
 F 

σh 

F 

where:
allowable bending stress against the long-time overturningb f crmoment exclusive of internal pressure (N/mm2)
f cr 0value given by eq. (3.25), (3.26) or (3.27) (N/mm2)
The value of b f cr is given by the following equations according to the value of r/t.
rE
≦ 0.274 tF

(3.35)

for(3.36)

Efor 0.274 
F
 

0.78

0.78

b

f cr =

F
1.5

r
E≦ ≦ 2.106 
t
F
 

0.78

0.78

rF 
 2.106 −   tE 

f
=
0
.
267
F
+
0
.
4
F
b cr


1.832





(3.37)

Efor 2.106 
F
 

0.78

≦

r
t

b

f cr =

1b σ cr , e
2.25

where:
the elastic bending buckling stress, based on lower boundb σ cr ,eformula proposed by NASA[3.1].
The value of b σ cr ,e is given by the following equation

Chapter3

25

t 
r


(3)





1/ 2

1 r
 
16
t



1 − 0.7311 − exp −b σ cr ,e = 0.6 E(3.38)



 
 

 

Allowable shear stress s f cr against the long-time lateral forceF− s f cr
σh
σh 
1
.
5
3
≦0.3 s f cr = s f cr +(3.39)
for 0≦
 
0.3
F
 F 
σ
(3.40)
for h ≧0.3 s f cr = F
F
1.5 3
where:
allowable shear stress against the long-time lateral loads f crexclusive of internal pressure (N/mm2)
The value of s f cr is given by the following equations according to the value of r/t.
(3.41)

E
0.204 
r
F
for ≦
0 .4
t
l
 
r

(3.42)

E
0.204 
F
for
0.4
l
 
 
r

(3.43)

E
1.446 
F
for
0 .4
l
 
r

0.81

s

f cr =

0.81

F
1.5 3
0.81

E
1.446 
r
F
≦ ≦
0 .4
t
l
 
 
r
0 .4
0.81

rl F 
1.446 −     
0.267 F 0.4 F t r E +
s f cr =


1.242
3
3




0.81

≦

r
t

s

f cr =

1s σ cr ,e
2.25

where:
ls

shear buckling wave length in the axial direction (mm)σ cr ,ethe elastic shear buckling stress, based on 95 percent
probability formula introduced by multiplying the Donnel’s
torsional buckling formula by 0.8

The value of s σ cr ,e is given by the following equation.

Chapter326

(3.44)3.7.3

3
 l  r 1 / 2  t
1 + 0.0239     s σ cr , e = 0.8
2
 r  t    l  r 1 / 2  r     r  t  

1/ 2

4.83 E

Allowable Stresses against the Short-Term Load
The allowable stresses against the short-term load exclusive of the seismic load
should be the values given by multiplying the allowable stresses against the longterm load shown in 3.7.2 by 1.5.

3.7.4

Allowable Stress against the Seismic Load
The allowable axial compressive, bending and shear stresses against the seismic
load are given by the following three terms (1), (2) and (3), respectively, which arealso depend on both σ h /F and r/t, and which should be applied to the ultimate
strength design.

(1)

Allowable axial compressive stress against the seismic load
σ(3.45)
for 0≦ h ≦0.3 c f cr = c f cr + 0.7 f crs − c f cr  σ h 
0.3
F
 F 
σ
(3.46)
for h ≧0.3 c f cr = f crs 1 − σ h F Fwhere:
allowable seismic axial compressive stress exclusive ofc f crfcrs

internal pressure (N/mm2)
basic value for determining the strength of elephant-footbulge (N/mm2)

The value of fcrs is given by the following equations according to the value of r/t.
(3.47)(3.48)(3.49)

rE≦0.069f crs = Ft
F
r
E
E
for 0.069 ≦ ≦0.807
F
t
F

for

r F

 0.807 − t E 
f crs = 0.6 F + 0.4 F  0.738 rE
for 0.807 ≦f crs = σ cr 0 , eFtwhere:

σ cr 0,e
Chapter3

value given by eq. (3.28)

27

The value of c f cr is given by the following equations according to the value of r/t.
r≦0.377  E tF

(3.50)

for(3.51)

for 0.377  E 

0.72

F

0.72

c

f cr = Fr≦ ≦2.567  E 
t
F

0.72

0.72

rF 
2
567
−
.

  
t E c f cr = 0.6 F + 0.4 F


2.190





(3.52)

for 2.567  E 

0.72

F

≦

r
t

c

f cr = c σ cr ,ewhere:
value given by eq. (3.32)
c σ cr ,e(2)Allowable bending stress against the seismic overturning momentσ(3.53)
for 0≦ h ≦0.3 b f cr = b f cr + 0.7 f crs − b f cr  σ h 
0.3
F
 F 
σ
(3.54)
for h ≧0.3 b f cr = f crs 1 − σ h F Fwhere:
allowable seismic bending stress exclusive of internalb f crf crs

pressure (N/mm2)
value given by the eq. (3.47), (3.48) or (3.49). (N/mm2)

The value of b f cr is given by the following equations according to the value of r/t.
r≦0.274  E tF

(3.55)

forfor 0.274  E 

0.78

(3.56)

for 2.106  E 

0.78

(3.57)

0.78b

f cr = F
0.78

r≦ ≦2.106  E 
t
F
F
0.78

rF 
 2.106 −   tE 
b f cr = 0.6 F + 0.4 F


1.832




F

≦

r
t

b

f cr = b σ cr ,ewhere:
value given by eq. (3.38)
b σ cr , e

Chapter328

(3)Allowable shear stress against the seismic lateral forceF− s f crσσh (3.58)
for 0≦ h ≦0.3 s f cr = s f cr + 3
 
0.3
F
 F 
σ
(3.59)
for h ≧0.3 s f cr = F
F3
where:
allowable seismic shear stress exclusive of internal pressures f cr
(N/mm2)The value of s f cr is given by the following equations according to the value of r/t.

(3.60)

E
0.204 
r
F
for ≦
0 .4
t
l
 
r

(3.61)

E
0.204 
F
For
0 .4
l
 
r

(3.62)

E
1.446 
F
for
0 .4
l
 
r

0.81

s

0.81

f cr =

F
3
0.81

E
1.446 
r
F
≦ ≦
0 .4
t
l
 
r
0 .4
0.81

rl F 
.
1
446
−

    
0. 6 F 0. 4 F t r E +
s f cr =

1.242
3
3 




0.81

≦

r
t

s

f cr = s σ cr ,ewhere:lshear buckling wave length in the axial direction (mm)valuegiven by eq. (3.44)σs cr ,e3.7.5

Remarks on the Evaluation of the Design Shear Stress(1) The design shear stress can be evaluated by the mean value within the range of theshear bucking wave length in the axial direction.
(2) The cylinder’s inside radius of curvature r and wall thickness t can be defined bytheir mean values within the range of the shear bucking wave length in the axial
direction.

Chapter329

(3) In the design of silos containing granular materials, the allowable shear stress can becalculated by setting l/r =1 and σ h =0.

Commentary:
1.

Empirical formula for buckling stresses

1.1 Case without internal pressure
The lower bound values of elastic buckling stresses which apply to cylindrical shells without
internal pressures are given by the following empirical formula.
(a)

(b)

(c)

elastic axial compressive buckling stresses (empirical formula by AASA[3.1])
1 r −
16 t    r (3.7.1)
σ
=
0
.
6
1
−
0
.
901
1
−
E
e

c cr , e /  t   
elastic flexural buckling stresses (empirical formula by AASA[3.1])
1 r −
16 t    r (3.7.2)b σ cr , e = 0.6 E 1 − 0.7311 − e /  t   
elastic shear buckling stresses [3.2] [3.3](3.7.3)

τ cr ,e =4.83E × 0.8
2

l
1 + 0.0239
r

3r   r / t   t 

l r 

r t 

Eq.(3.7.3) is given by the torsional buckling stress derived by Donnell multiplied by the reduction
factor of 0.8. These formula correspond to the 95％ reliability limit of experimental data.
Based on the elastic buckling stresses, the nonlinearity of material is introduced. In the region
where the elastic buckling stresses exceeds 60％ of the yield point stress of material ( σ y or τ y ),the buckling mode is identified to be nonlinear. In Fig.3.7.1, the relationship between the buckling
stress and the r/t ratio (the buckling curve) is shown. The buckling curve in the nonlinear range ismodified as shown in Fig.3.7.1. (r/t)I is the (r/t) ratio at which the elastic buckling stress σ cr ,ereaches 60％ of the yield stress. (r/t)II denotes the (r/t) ratio at which the buckling stress exceeds
the yield point stress. In the region of (r/t) ratios less than (r/t)II, the buckling stresses are assumed
to be equal to the yield point stress. In the region of (r/t) ratios between (r/t)II and (r/t)I, thebuckling stresses are assumed to be on a line as shown in the figure. (r/t)II was determined on thecondition that the elastic buckling stress at (r/t)II reaches 6.5 times to 7.5 times as large as the yield
point stress.

Fig.3.7.1 Buckling curve

Chapter330

(r/t)I and (r/t)II are determined according to the stress conditions as shown in Table 3.7.1.
The occurrence of buckling under combined stresses can be checked by the following criteriaσcσ+ b = 1
σ
σ
c cr
b cr

(3.7.4)
or

τ

=1

τ cr

Aotations:σcσbmean axial compressive stresscompressive stress due to bending momentτshear stressEq.(3.7.4) implies that the axial buckling and the flexural buckling are interactive, while the shear
buckling is independent of any other buckling modes. The experimental data are shown in
Fig.3.7.2 . Figure (a) shows the data obtained by Lundquist on cylindrical shells made ofduralumin all of which buckled in the elastic range[3.4]. Figure (b) indicates the test results on
steel cylindrical shells[3.3]. The ordinate indicates σ c / c σ cr + σ b / b σ cr and the abscissa indicatesτ / τ cr . The buckling criteria which correspond to Eq.(3.7.4) is shown by the solid lines. Althougha few data can not be covered by Eq.(3.7.4), Eq.(3.7.4) should be applied as a buckling criterion
under the seismic loading considering that the estimate of Ds-values are made conservatively. For
short-term loadings other than the seismic loading, the safety factor of 1.5 is introduced in theestimate of the elastic buckling stresses to meet the scatter of the experimental data.
Table 3.7.1

Limit r/t ratios

Stress Condition
Axial
Compression
Bending
Shearing

Chapter331

(a) Duralumin cylindrical shell(b) Steel cylindrical shell

Fig.3.7.2 Comparison between test results and predictions

1.2 Case with internal pressure
The effects of the internal pressure on the buckling of cylindrical shells are twofold :・The internal pressure constrains the occurrence of buckling.
・The internal pressure accelerates the yielding of shell walls and causes the elephant foot bulge.
The effect of the internal pressure is represented by the ratio of the tensile hoop stress, σ h , to theyield point stress, σ y , (σ h / σ y ) [3.5].
Referring to experimental data, it was made clear that the buckling stresses under axial
compression and bending moment can be expressed by the following empirical formula[3.6] [3.7]:
(0.7σ cr 0 − σ cr )(σ h / σ y ) 
σ cr = σ cr +0.3
for σ h / σ y > 0.3, σh σ cr = σ cr 0 1 − σy for σ h / σ y ≤ 0.3,

(3.7.5)

Aotations:

σ cr
σ cr 0buckling stress in case without internal pressurecompressive buckling stress in axi-symmetric mode

Fig.3.7.3 Influence of internal pressure on buckling stress

Chapter332

Eq.(3.7.5) is shown in Fig.(3.7.3). In the region below σ h / σ y = 0.3 , the buckling occurs ininextensional mode, whereas beyond σ h / σ y = 0.3 , the buckling occurs in elephant foot bulge mode.

σ cr 0 for cylinders which buckle elastically was found to be expressed byt 
3(1 − ν )  r 
Aotations:
Poisson's ratioν

σ cr 0, e =(3.7.6)0.8 E

2Eq.(3.7.6) is the theoretical buckling stress of axially compressed cylinders multiplied by aempirical reduction factor of 0.8. The buckling stress in the inelastic range, σ cr 0 , can be related
to σ cr 0 ,e , similarly to c σ cr ,e and b σ cr , e as follows:
for

r
E
≤ 0.069 tF
(3.7.7)

E r
for 0.807  ≤F t(3.7.8)σ cr 0 = σ yσ cr 0 = σ cr 0 ,e

E r
E
for 0.069  < < 0.807 
F t
F
σ cr 0 is linearly interpolated between the values given by Eq.(3.7.7) and Eq.(3.7.8)

Fig.3.7.4

Chapter3

Comparison between test results and predictions for axial loading

33

The applicability of Eq.(3.7.5) is demonstrated in the following:i) Axially compressed cylinders experimental data and the results of analysis are shown in
Fig.3.7.4, together with Eq.(3.7.5) [3.7]. The abscissa indicates the hoop stress ratio, σ h / σ y , andthe ordinate indicate the ratio of the buckling stress to the classical theoretical buckling stress, σ c1 .
Shell walls were perfectly clamped at both ends. Test data are indicated by the symbols of ◆ thatindicates the buckling under diamond pattern, ● that indicates the buckling under elephant footbulge mode and ▲ that is the case of simultaneous occurrence of diamond pattern and elephantfoot bulge mode.
Curves with a specific value of ∆ are theoretical values for elephant foot bulge mode. ∆ is theinitial imperfection divided by the thickness of wall. The chained line going up right-handedly is
the theoretical curve of diamond pattern of buckling derived by Almroth[3.8]. The polygonal line A
is the design curve for seismic loading given by Eq.(3.7.5). The polygonal line B is the design
curve for loading other than the seismic loading.
ii) Cylinders subjected to shear and bending[3.6]
Referring to experimental results on cylindrical shells subjected to shear and bending, it was madeclear that the shear buckling pattern becomes easily localized with a faint presence of the internalpressure.
In Fig.3.7.5, a shear buckling pattern under σ h / σ y = 0.05 is shown. Thus, it is understood thatthe internal pressure makes the substantial aspect ratio, l / r in Eq.(3.7.3), reduced drastically.
At σ h / σ y = 0.3 , the shear buckling disappears.
Considering these facts, the shear buckling stress under the internal pressure can be summarized asfollows:(3.7.9)for σ h / σ y ≤ 0.3 , τ cr = τ cr +(τ − τ )(σ / σ )
y

cr

h

y

0. 3

for σ h / σ y > 0.3 , τ cr = τ yAotations:shear yield-point stressτyFig.3.7.5

Chapter3

Buckling wave pattern of shear buckling
σ h / σ y = 0.05, r / t = 477, l / r = 1.0, σ y = 32.3kA / cm2

34

(t/cm 2 =9.8kN/cm 2 )(t/cm 2 =9.8kN/cm 2 )σm(t/cm2=9.8kN/cm2)bbb

σm(t/cm2=9.8kN/cm2)

σm(t/cm2=9.8kN/cm2)bb

σm(t/cm2=9.8kN/cm2)

σm(t/cm2=9.8kN/cm2)

Fig.3.7.6 Influence of internal pressure on buckling(in case of small internal pressure)Fig.3.7.7 Buckling strength under combined internal pressure and flexural-shear loading

Chapter335

In Fig.3.7.6, test results in the range of small values of σ h / σ y are compared with the estimate
made by Eq.(3.7.5) and Eq.(3.7.9). The ordinate indicates the bending stress at the lower edge ofcylinder. Therefore, the shear buckling stress is converted to the equivalent bending stress by thefollowing equation.
h(3.7.10)σ cr = τ crr
Aotations:hdistance between the lower edge of shell and the point wherehorizontal load is applied
Eq.(3.7.5) is shown by solid lines and Eq.(3.7.9) is shown by broken lines. As shown in the figure,cylindrical shells subjected to internal pressures are hardly influenced by shear buckling.
In Fig.3.7.7, test data under wider range of parameters are shown, predictions made by Eq.(3.7.5)
and Eq.(3.7.9) are shown by solid lines. The point shown by ○ indicates the inextentional mode ofbuckling, and the point shown by ● indicates the elephant foot bulge mode of buckling.
2.

Safety factor on buckling stresses
Safety factors for various loading conditions are set up as shown in Table 3.7.2. The level of safety
factor depends on the probability of occurrence of the specified design loading.
Table 3.7.2 Safety factor for buckling
Safety factor for buckling

Loading
conditionLong term
condition
Short term
(except
earthquake)

Linear interpolation
Linear interpolation

Earthquake
condition3.

Supplementary comment on shear buckling
When geometrical parameters and structural parameters change within a range of shear buckling
length, the mean values of them can be taken according to the usual practice of engineering.
Generally, the buckling strengths of silos containing granular contents are influenced by theinternal pressure and stiffness associated with the properties of contents, and the shear buckling
length, l / r is considered to be less than unity. Still, the exact estimate of internal pressure underseismic loading is difficult. As a tentative measure of compromise, the constraining effect ofinternal pressure is to be disregarded on the condition that l / r is fixed to be unity.

4.

Deformation characteristics after buckling
The general load-deformation curve of cylindrical shells subjected to the monotonic horizontal
force is schematically shown in Fig.3.7.8. The ordinate indicates the maximum bending stress, b σ .
The abscissa indicates the mean inclination angle, θ . In the post-buckling range, the deformation
develops with a gradual decrease of resistant force, finally approaching a horizontal asymptote.
The level of the asymptote depends on the level of the axial stress, σ 0 . When σ 0 / c σ cr becomeslarger than 0.2, the asymptote disappears, and the resistant force decreases monotonously as thedeformation increases.
When σ 0 / c σ cr remains less than 0.2, cylindrical shells exhibit a stable energy absorption capacity

Chapter336[3.3] [3.6] [3.10] [3.11].
The hysteretic behavior of cylindrical shells under repeated horizontal forces is shown inFig.(3.7.9).
The basic rule which governs the hysteretic behavior is described as follows.

・

Definitions to describe the hysteretic rule:a)
The loading path is defined by b σdθ > 0 .
b)

・
・
・
・

The unloading path is defined by b σdθ < 0 .

The b σ − θ relationship under the monotonic loading is defined to be the skeleton curve.
The unloading point is defined to be the point which rests on the skeleton curve andterminates the loading path.
The initial unloading point is defined to be the point at which the buckling starts.
The intermediate unloading point is defined to be the point which does not rest on theskeleton curve and terminates the loading path.

On the assumption that the loading which reaches the initial loading point has been already madein both positive and negative directions, the hysteretic rule is described as follows.
a)

The loading path points the previous unloading point in the same loading domain. Afterreaching the previous unloading point, the loading path traces the skeleton curve.
b) The unloading path from the unloading point points the initial unloading point in the
opposite side of loading domain. The unloading path from the intermediate unloading point
has a slope same as that of the unloading path from the previous unloading point in the sameloading domain.
These hysteretic rule has been ascertained to apply even for the case of σ h / σ y ≠ 0 .

Fig.3.7.8

Load-deformation relationship under monotonic loadingskeleton curveloading pathunloading path

parallelinitial unloading point
unloading pointintermediate unloading point

Fig.3.7.9 Hysteretic rule

Chapter337

5.

Dη -values for cylindrical structuresDη -values for structures with short natural periods which has been obtained analytically using

the load-deformation characteristics shown in Fig.3.7.9 are shown in Fig.3.7.10[3.12].
Essential points in the analysis are summarized as follows.
5.1

Dη -values in case that uplift of the bottom of the structure is not allowed by restraint of ananchoring system:
1)

2)
3)

Assuming that the structure is one-mass system and equating the structural energy
absorption capacity to the energy input exerted by the earthquake, the maximum deformationof the structure is obtained.
The energy input to the structure is obtained by considering the elongation of the vibrational
period of structure.
The Dη -values is defined by(3.7.11)Dη =4)the level of earthquake which causes the buckling deformation, δ cr
the level of earthquake which causes the maximum deformation of (1 + µ )δ cr

Aotations:δµ = m −1
δ cr

nondimensionalized maximum deformation

δ crbuckling deformation

The assumed load-deformation curve in monotonic loading is one that exhibits aconspicuous degrading in strength in the post-buckling range as shown in Fig.3.7.10. q in
the figure denotes the nondimensionalized level of the stationary strength in the post buckling
range.
When we select the value of Dη to be 0.5, the pairs of values of ( q , µ ) are read from
Fig.3.7.10 as follows : (0.6, 1.0), (0.5, 2.0), and (0.4, 4.0). See Table 3.7.3.
These q − µ relationships which corresponds to Dη = 0.5 are compared with thenondimentionalized load-deformation curves obtained experimentally in Fig.3.7.11 andFig.3.7.12. The q − µ relations are written as white circles in these figures. The
experimental curves lie almost above the q − µ relations for Dη = 0.5 (see Table 3.7.3). It
implies that the actual structures are more abundant in energy absorption capacity than theanalytical model equipped with the load-deformation characteristics shown in Fig.3.7.10.
Thus, to apply the value of 0.5 for Dη -values of cylindrical structures influenced by bucklingseems reasonable. When the axial compressive stress becomes large, however, the Dη values should be kept larger than 0.5 since the strength decreases drastically in the post
buckling range.
Consequently, the value of Dη is determined as follows:

Chapter338

Fig.3.7.10 Dη -valuesfor σ 0 / c σ cr ≦0.2,(3.7.12)Dη = 0.5

for σ 0 / c σ cr >0.2,(3.7.13)Dη = 0.7

Table 3.7.3 q(=Q/Qcr) and µ (= [δ − δ cr ] / δ cr ) relationship (in case of δ 0 / δ cr ≦ 0.2 , Dη = 0.5 )
µ
q
0.6
1.0
0.5
2.0
0.4
4.0

Fig.3.7.11 Load-deformation relationshipwhich corresponds to Dη = 0.5( σ h /σ y = 0 )Chapter339

Fig.3.7.12 Load-deformation relationshipwhich corresponds to Dη = 0.5( σ h /σ y ≠ 0 )5.2

Dη -values in case that uplift of the bottom of the structure is allowed without an anchoring

system:The Dη -value defined in Eq. (3.6.15) of sec. 3.6 should be expressed as follows:(3.7.14)Dη =1

 Wp 

1 + 

 We 

The cumulative plastic strain energy Wp induced by buckling of the shell plate is as indicatedin Eq. (3.6.7). However, the elastic vibration energy We should be of the strain energy of thecylindrical shell added by the strain energy due to the uplift of the bottom of the shell. In casethat uplift of the bottom of the structure is not allowed, replacing Dη and We with Dη f and Wef,
respectively, Eq. (3.7.14) should be expressed as follows:(3.7.15)

1

Dη f =

 Wp 

1+ 
 Wef 


In case that uplift of the bottom of the structure is considered, Eq. (3.6.8) should be expressed asfollows:(3.7.16)

We =

Qyδ y
2

Whilst in case that uplift of the bottom of the structure is not considered, Wef should be expressed

as follows:(3.7.17)Wef =Q yδ y  T f 
2  Te 

2

From (3.7.14) to (3.7.17) Dη value, in case that uplift of the bottom of the structure is allowed,
should be given in Eq. (3.7.18).
(3.7.18)

Chapter3

Dη =

1
1

  1
 T f 2  2  
1 + 
−1
 Te    Dη f 2




40

Aotations:Dη f ,WefDη , We in case that uplift of the bottom of the structure is

not allowed.

Tfnatural period in case that uplift of the bottom of thestructure is not allowed.
natural period in case that uplift of the bottom of thestructure is considered.

Te
In general as Te > Tf , Dη > Dη f .
6.

B-value for cylindrical structures
In thin steel cylindrical structures, the buckling starts when the maximum stress reaches the
buckling stress, σ cr .
Therefore, the stress redistribution due to yielding is hardly considered. Thus, the value of B forsteel cylindrical shells should be set as follow:.
for r / t ≧( r / t )II ,(3.7.19)

B=1.0
Aotations:( r / t )II：refer to Table 3.7.1

3.8

Earth Pressures

3.8.1

Earth Pressure in Normal Condition

Earth pressure in normal condition should be determined according to“Recommendations for Design of Building Foundations, 2001 - Architectural
Institute of Japan ”.
------------------------------------------------------------------------------------------------------------------Commentary:
Earth pressures in the normal condition for sand and clay shall be calculated in accordance withthe “Recommendations for design of building foundations (2001)” compiled by the Architectural
Institute of Japan (AIJ).

3.8.2

Earth Pressure during Earthquakes
Earth pressures which act on the underground walls of tanks during earthquakes
may be considered to be either active or passive and should be calculated inaccordance with equations (3.63) and (3.64) for active earth pressures, and (3.65)
and (3.66) for passive earth pressures.
(a)

Active earth pressure(3.63)

PEA = ( 1 − KV )( γz +(3.64)K EA =

Chapter3

ω

)K EA

cos θ

cos 2 ( φ − θ )  sin( φ + δ ) sin( φ − θ ) 1 / 2 cos θ cos( δ + θ )1 +  cos( δ + θ )
 
 

41

2

(b) Passive earth pressure(3.65)

PEP = ( 1 − KV )( γz +(3.66)K EP =

ω
cosθ

)K EP
cos 2 ( φ − θ )  sin( φ + δ ) sin( φ − θ ) 1 / 2 cos θ cos( δ + θ )1 −  cos( δ + θ )
 
 

2

where:
PEA
KV
γ
z
ω
θ
K
(3.67)

K=
KH
KEA
φ
δ
PEP
KEP

horizontal active earth pressure during earthquakes acting
on vertical wall at depth z (N/mm2) - refer also to Figure 3.1,
vertical underground seismic intensity given in chapter
3.6.2.2,
unit weight of soil (N/mm3),
depth from the ground surface (mm),load pressure (N/mm2), - refer also to Figure 3.1,
tan-1K,
combined seismic coefficient, see equation (3.67)
KH
1 − KVhorizontal underground seismic intensity given in chapter
3.6.2.2,
coefficient of active earth pressure during earthquakes,
internal friction angle of soil,
wall friction angle,
horizontal passive earth pressure during earthquakes acting
on vertical wall at depth z (N/mm2), - refer also to figure 3.1
coefficient of passive earth pressure during earthquakes.

ω (N / mm2)

PEA , PEP (N / mm2)tank

Fig.3.1 Definition of symbols------------------------------------------------------------------------------------------------------------------Commentary:
Although the actual mechanism of soil pressure during earthquakes is complicated and unknown to someextent, Mononobe and Okabe presented their formula for earth pressure during earthquakes using
seismic intensity and Coulomb’s earth pressure theory of ultimate slip. Vertical tank wall and horizontal
ground surface are assumed in equations (3.8.1) and (3.8.2) for simplicity.

Chapter342

(a)Seismic active earth pressure resultant(3.8.1)(b)γ

ωPEA ' = (1 − KV )( +) K EA H 2
2 H cosθSeismic passive earth pressure resultant(3.8.2)γ

ωPEP ' = (1 − KV )( +) K EP H 2
2 H cosθ
where:PEA 'seismic active earth pressure resultant acting on the wall

K EA

with depth H (A/mm)
see eq. (3.64) in chap. 3.8.2

KV

PEP '

see eq. (3.15) in chap. 3.6.2.2
unit weight of soil (A/mm3)
depth from the ground surface to the tank bottom (mm)imposed loads (A/mm2)
seismic composite angle (=tan-1K)
seismic passive earth pressure resultant acting on the wall

K EP

with depth H (A/mm)
see equation (3.66) in chap. 3.8.2

γ
H
ω

θ

As shown above, Mononobe and Okabe’s formula gives the resultant acting on the tank wall for aunit width soil wedge above the slip surface. For convenience, this recommendation transforms
equations (3.8.1) and (3.8.2) to pressures per unit area of the tank wall as in equations (3.63) and(3.65). If these equations are integrated in the interval z = 0~H, equations (3.8.1) and (3.8.2) will
be obtained.
Equations (3.63) ~ (3.66) give pressures acting in the inclined direction with angle δ from thehorizontal plane. But the angle of friction between the wall and the soil may be neglected, and thepressures may be regarded as acting in the normal direction to the wall.
The equations in this chapter give ultimate earth pressure, but they have limitations in that theirbasis is static equilibrium and they depend on assumed seismic intensity. Therefore additional
methods, such as response displacement, seismic intensity, or dynamic analysis methods should
also be used in the design for seismic loads.

References:
[3.1]
[3.2]
[3.3]
[3.4]
[3.5]
[3.6]

[3.7]

[3.8]

AASA SP-8007: Buckling of Thin-walled Circular Cylinders, Aug. 1968
Donnell, L. H.: Stability of Thin-walled Tubes under Torsion, AACA, TA479, 1933
Akiyama, H., Takahashi, M., Hashimoto, S.: Buckling Tests of Steel Cylindrical Shells Subjected to
Combined Shear and Bending, Trans. of Architectural Institute of Japan, Vol. 371, 1987
Lundquist, E. E.: Strength Test of Thin-walled Duralumin Cylinders in Combined Transverse Shear andBending, AACA, TA523, 1935
Yamada, M.: Elephant's Foot Bulge of Cylindrical Steel Tank Shells JHPI, Vol. 18, Ao. 6, 1980
Akiyama, H., Takahashi, M., Aomura, S.: Buckling Tests of Steel Cylindrical Shells Subjected to
Combined Bending, Shear and Internal Pressure, Trans. of Architectural Institute of Japan, Vol. 400,
1989
Yamada, M., Tsuji, T.: Large Deflection Behavior after Elephant's Foot Bulge of Circular Cylindrical
Shells under Axial Compression and Internal Pressure
Summaries of Technical Papers of Annual Meeting, AIJ, Structure-B, pp. 1295-1296, 1988
Almroth Bo O., Brush D. O.: Postbuckling Behavior of Pressure Core-Stabilized Cylinders under Axial
Compression. AIAA J., Vol. 1, pp.2338-2441, 1963

Chapter343[3.9]

Yamada, M.: Elephant's Foot Bulge of Circular Cylindrical Shells
The Architectural Reports of the Tohoku University, Vol. 20, pp. 105-120, 1980
[3.10] Shibata, K., Kitagawa, H., Sampei, T. Misaji, K: The Vibration Characteristics and Analysis of After
Buckling in Metal Silos, Trans. of Architectural Institute of Japan, Vol. 367, 1986
[3.11] Kitagawa, H., Shibata, K.: Studies on Comparison of Earthquake - Proof of Metal and FRP Silos,
Summaries of Technical Papers of Annual Meeting, AIJ, Structure - B, 1988
[3.12] Akiyama, H.: Post Buckling Behavior of Cylindrical Structures Subjected to Earthquakes, 9th SMiRT,
Lausanne, 1987
[3.13] Kato, B., Akiyama, H.: Energy Input and Damages in Structures Subjected to Severe Earthquakes,
Trans. of Architectural Institute of Japan, Vol. 235, 1975

Chapter3444Water Tanks

4.1

Scope

4.1.1

GeneralThis chapter is applicable to the design of potable water storage tanks and theirsupports.

4.1.2

Materials
Potable water storage tanks are normally constructed using reinforced concrete, prestressed concrete, carbon steel, stainless steel, aluminium, FRP, timber, or any
combination of these materials.

4.1.3

Notations
B
hi
Qei
W

ratio of the horizontal load-carrying capacity to the short term allowablestrength.
height of the i-th mass from ground level (m).
design storey shear force for designing allowable stresses (N).
design weight imposed on the base of the structure exclusive of
convective mass (N). See chap. 4.2.4.4.

4.2

Structural Design

4.2.1

Design Principles

4.2.1.1

Potable water storage tanks should be located such that the contents are not subject
to insanitary pollution. They should also be located to enable easy access forinspection, maintenance, and repair.

4.2.1.2

The shape of a water tank should be as simple as possible, be symmetrical about its
axis, and mechanically clear of any unnecessary obstructions.

4.2.1.3

Due consideration should be given to stress concentrations at openings or
penetrations for pipework.

4.2.1.4

Suitable measures should be taken to prevent a sudden internal pressure drop causedby any uncontrolled release of water.

4.2.1.5

The foundations for water tanks should be designed such that the total forces of theoverall structure are transferred to the ground as smoothly as possible.

4.2.2

Allowable Stresses and Material Constants
Allowable stresses and material constants used in the construction of water tanks
and their support structures should be in accordance with the applicable codes andstandards, and as indicated below.

Chapter 445

(1) Long term allowable stresses of reinforcing bars in a reinforced concrete structure
should be properly defined after taking into account the allowable cracking width.
(2) Allowable stresses of stainless steel should be defined in accordance with the“Design Standard for Steel Structures – Based on Allowable Stress Concept -” - AIJ
2005. The standard F value should be taken as 0.2% strength as stipulated by JIS.
(3) The standard value F, which defines yielding stress of all welded aluminium alloy
structure, should be defined in accordance with the “Aluminum Handbook – Japan
Aluminium Association (2001)” etc.
(4) Allowable stresses and constants of FRP structures should be defined by reference tothe material characteristics obtained from tensile and flexural strength tests,(“Method of tension test for FRP” - JIS K 7054, and “Method of flexural test of
FRP” - JIS K 7055), after taking into account degradation due to aging, etc.
4.2.3

Loads

4.2.3.1

Loads for structural design should be as shown in chapter 3 in addition to theprovisions of 4.2.3.2 below and so on.

4.2.3.2

Water tanks should be considered as being full for the design calculations. However,a partial load for the tank contents may be considered when there is a dangerous
combination of loads that would increase the overall loading.

4.2.4

Seismic Design

4.2.4.1

Seismic design of water tanks should be classified below in accordance with 3.6.1.7:-

Small water tanks used at private houses
Essential water tanks for public use
OthersⅠ
Ⅲ
Ⅱ

4.2.4.2

Design seismic loads should be obtained from the simplified seismic coefficientmethod or the modal analysis which are shown in 3.6.1.2, 3.6.1.3 and 4.2.4.5 except
for 4.2.4.3.

4.2.4.3

Design yield shear forces should be determined based on the horizontal seismic
coefficient of 1.5 for the water tanks placed on the roof.

4.2.4.4

Contents in the water tank should be considered in two different masses, i.e. theimpulsive mass of the tank contents that moves in unison with the tank and the otherconvective mass of the tank contents that acts as sloshing liquid. See the commentaryof chap. 7.1 and 7.2 for more information.

4.2.4.5

Calculation of design shear forces based on Simplified Seismic Coefficient Method:

Chapter 446

(1) The design yield shear force, Qd, at the ground base of the water tank should be
calculated by the equations (3.1) by use of the simplified seismic coefficientmethod.

Where:
W = design weight imposed on the base of the structure exclusive of convective
mass (N)Ds = structural characteristic coefficient, as is given on Table 4.1
Table 4.1 Ds Value(a) Elevated Water Tanks
Structure
Metallic Structure
Reinforced Concrete Structure
Steel Reinforced Concrete Structure
Pre-stressed Concrete Structure

Type of Structure
Moment Frame
Truss/Bracing
Plate/Shell
Moment Frame

Ds
0.40
0.40
0.5~0.7
0.40

Plate

0.45

(b) Water Tanks Installed Directly on the Foundation or Ground

4.3

Structure

Ds

Other than below
Reinforced Concrete Structure
Pre-stressed Concrete Structure

0.55
0.45

Reinforced Concrete Water Tanks
The design of reinforced concrete water tanks should be as indicated in 2.2, 2.3, 2.4
and 4.2, and the “Standard for Structural Calculation of Reinforced Concrete
Structures” - AIJ 2010. However, the following provisions should also be taken into
consideration.
(1) Concrete for water tanks should be as specified in Chapter 23: Watertight Concrete JASS 5, AIJ 2009.
(2) The minimum ratio of reinforcing bars of side walls and bottom slab should be morethan 0.3% in both orthogonal directions.
(3) The minimum concrete covering of reinforcing bars should be as given in Chapter
3.11: Cover - JASS 5, AIJ 2009.
(4) Any required water proofing should be applied to the surface that is in contact withthe stored water.

Chapter 447

(5) The minimum concrete covering of reinforcing bars to the surface in contact withstored water should be 30mm, and in the case of light duty water proof finishing,e.g., mortar, or in cases where no water proof finishing is applied, the minimum
cover of reinforcing bars should be 40mm and 50mm respectively.

4.4

Pre-stressed Concrete Water Tanks
Pre-stressed concrete water tanks should be as stipulated in 2.2, 2.4, and 4.2 of thisdocument, and the “Standard for Structural Design and Construction of Pre-stressed
Concrete Structures” - AIJ 1998. Reinforced concrete components should be as
noted in 4.3. In addition, the following provisions should be taken into
consideration.
(1) For usual use conditions, cross sections of members should be designed in full prestressing condition. Stresses from water pressure should be considered as live loads.
(2) The additional loads and stresses caused by pre-stressing should also be taken into
consideration.
(3) Thermal stresses and additional stresses caused by dry shrink strain should beconsidered.

4.5

Steel Water Tanks
The design of steel water tanks should conform to the “Design Standard for Steel
Structures - Based on Allowable Stress Concept - ” 2005 edited by the Institute, andthe following items.
(1) Cylindrical shell design of the water tanks should give investigation on the bucklingof the cylindrical shell conforming to section 3.7.
(2) The water tanks and pipeline design should give full preventative measures tocorrosion.

4.6

Stainless Steel Water Tanks
The design of stainless steel water tanks should conform to the “Guideline of AntiSeismic Design and Construction of Building Equipments (2005)” edited by the
Building Center of Japan (BCJ) etc. The design of stainless steel water tanks support
section should conform to the “Design and Construction Standard on the StainlessSteel Building Structure, and the Explanation (2001)” edited by the Stainless SteelJapanese Society of Steel Construction, and the following items.
(1) Cylindrical shell design of the water tanks should give investigation on buckling of
the cylindrical shell conforming to section 3.7

Chapter 448

(2) The water tanks and pipeline design should give full preventative measures tocorrosion.
4.7

Aluminium Alloy Water Tanks
The design of aluminium alloy structure should conform to section 2.5, 4.2 and the“Commentary, Design and Calculation Examples of Technical Standard on the
Aluminium Alloy Building Structure (2003) by BCJ etc.”, and to the following
items.
(1) Cylindrical shell design of the water tanks should give investigation on buckling of
the cylindrical shell conforming to section 3.7(2) The water tanks and pipeline design should give full preventative measures tocorrosion.
(3) When the stored water is potable water, adequate lining should be applied to theinner wall surface of the tank.

4.8

FRP Water Tanks
The acceleration response of the story at which FRP water tanks is installed should
be taken into consideration on the seismic design of the tank. They are usually
installed at a higher elevation.
The pipe connections and the anchoring system should be designed with carefulness.
Both of them are known to be relatively easily suffered damages during earthquakes.
Consideration should also be given to the sloshing response of the water during
earthquake which may cause to damage the ceiling panels.
Materials used in the manufacture of FRP water tanks should meet the requirementsof paragraphs, 2.6 and 4.2, and should be designed in accordance with a recognised
code, e.g., the Japan Reinforced Plastics Society “Structural Calculation Method for
FRP Water Tanks (1996)”.
However, the following precautions should be always considered.
(1) The wall thickness of the tank should be more than 3mm except for the manhole
cover plate.
(2) The design should give consideration to the anisotropic physical properties of the
proposed FRP materials.
(3) The design should give consideration to the thermal properties of the proposedmaterials in case the tank contains the hot water.
(4) The design should give consideration to the corrosion prevention in case the tank isreinforced with metals.

Chapter 449

(5) The design should give consideration to the measure to prevent formation of algae
due to the transmission of light.
(6) Materials selected for potable water tanks should be of the ones not impair thequality of the stored water, and should be resistant to the chemicals used forcleaning the tank.
(7) The structural design should give consideration to the expected loads and thephysical characteristics of the proposed FRP materials.

4.9

Wooden Water Tanks
This section is omitted in this English version.

------------------------------------------------------------------------------------------------------------Commentary:
In addition to designing for structural integrity, the design of FRP water tanks should include
all hygiene requirements for potable water.
The following points should be given due consideration in the design of FRP water tanks:(1)

With disregard to the calculated strength requirements and the derived wall thickness
requirements, a minimum wall thickness of 3mm should be applied to all FRP watertank designs.

(2)

The orientation of glass fibre has a significant effect on the strength and stiffness of
FRP. Therefore, due consideration should be given to the anisotropy of the glass fibre
reinforcement in the design of FRP water tanks.

(3)

The effect of temperatures between - 40°C and + 40°C on the material characteristics
of FRP is minimal. As the normal temperature range of potable water is + 4°C to +30°C the effect of temperature on FRP materials is therefore insignificant. However,at temperatures above + 60°C the strength and stiffness of FRP is impaired, therefore
the design of water tanks for use at high temperatures should take into account thisloss of strength and stiffness.

(4)

In cases where a combination of FRP and other materials is necessary, e.g.,galvanised steel or stainless steel reinforcement, preventative measures should betaken to prevent corrosion of the non FRP reinforcement by contaminants or otherchemicals, e.g., chlorine, in the water .

(5)

In order to prevent the formation of algae inside water tanks, and therefore to
maintain water quality, the intensity ratio of illumination (the ratio of intensity of lightbetween the inside and outside of the water tank), should not exceed 0.1%. Careshould therefore be taken in the design and placement of manholes, air vents, roof
panels, etc., which may provide points of ingress for light.

(6)

Materials used in the construction of potable water tanks should meet the relevanthygiene standards. Particular attention should be paid to the prevention of materials
being used in the manufacturing process which may be soluble in water and which canaffect the quality of potable water. As the inside of water tanks are washed with an

Chapter 450aqueous solution of sodium hypochlorite, only materials which are known to beresistant to this cleansing agent should be selected for water tanks.
(7)

Chapter 4

FRP does not have the distinct yield point, and neither does it have the ductility
similar to metallic materials. Therefore the structural design should take into account
this deformation and strength characteristics of the FRP structure.

51

5.

Silos

5.1

Scope

5.1.1

This chapter covers the structural design of upright containers and their supports(silos) for storing granular materials.

5.1.2

Definitions
silos
containers
silo wall
hopperNote:

the generic name for upright containers used for storing granularmaterials and their supports.
tanks for storing granular materials.
the upright wall of the containers.
the inclined wall of the containers.

The term granular materials is taken to include “fine particles”.

5.1.3

The structural materials of the silos which are considered in this chapter are
reinforced concrete and steel.

5.1.4

Notations

Cd

Cf
Ci
CL
d
H
hm
Ks
lc
d:m
Pf

Chapter 5silo wallChorizontal area of the silo (m2)ratio of the horizontal load-carrying capacity of the structure to the shortterm allowable yield strength
design yield shear force coefficient at thebase of silo
overpressure factor, used to consider
increases of pressure occurring during
discharge, converting from static pressure to
design pressure
safety factor of friction force, used to longterm design friction force between silo wallassumed freeand granular materials
surface
impact factor, used to consider pressure
increase due to sudden filling, convertingfrom static pressure to design pressure
ratio of the local pressure to the design
horizontal pressure
diameter (inside) of silo (see Fig.5.1) (m)height of the silo (see Fig.5.1) (m)effective height of material above bottom ofhopper assuming top of stored material islevelled (see Fig.5.1) (m)ratio of Ph to Pv
perimeter of horizontal inside cross section ofcontainer (or hopper) (m, cm)
vertical design force per unit length of silowall horizontal section (N/mm)
friction force of unit area of silo wall surface
Fig.5.1 Silo Section
skirt

A
B

52

Ph
Pv
Pa
dPh
dPL
dPv
dPa
rw
x

α
γ
µf
φi
φr

at depth x (N/mm2)
horizontal pressure of unit area at depth x due to stored material (N/mm2)
vertical pressure of unit area at depth x due to stored material (N/mm2)
pressure normal to the surface of inclined hopper wall (N/mm2)
design horizontal pressure of unit area (N/mm2)
design local horizontal pressure of unit area (N/mm2)
design vertical pressure of unit area (N/mm2)
design normal pressure of unit area to hopper wall (N/mm2)
hydraulic radius of horizontal cross section of storage space (m)depth from surface of stored material to point in question (m)angle of hopper from horizontal (°)
weight per unit volume for stored material (N/mm3)
coefficient of friction between stored material and wall or hopper surface
internal friction angle of stored material ( ° )repose angle of stored material ( ° )5.2Structural Design

5.2.1

General5.2.1.1

Silos and their supports should be designed to contain all applicable loads taking
into account the properties of stored materials, the shape of the silos, methods ofmaterial handling, etc.

5.2.1.2

The shape of the silo should be as simple as possible, be symmetrical about its axis,
and should have structural members which are proportioned to provide adequate
strength.

5.2.1.3

The foundations for silos should be designed to support stresses from the upper
structural members of the silos and their supports.

5.2.1.4

The design should include measures to prevent dust or gas explosions andexothermic reaction of stored materials.

5.2.1.5

In the event that fumigation with insecticides is required, the silo should be air tightin accordance with the applicable regulations.

5.2.1.6

The internal surfaces of silos which are used for storing food materials or feed stock
should meet the relevant food safety and sanitary regulations.

5.2.1.7

Physical property tests using actual granular materials are expected to find weight
per unit volume γ , internal friction angle φi ; and deformation characteristics(dependency of rigidity and damping ratio on strain level)

Chapter 553

5.2.2Loads

5.2.2.1

GeneralWith the exception of loads for structural design shown in this chapter, all otherloadings should, in principle, be based on chapter 3.

(1) Silos, hopper walls, and bottom slabs should be designed to resist all applicableloads, and should include forces from stored materials, static pressures, dynamic
over-pressure and under-pressure during charge and discharge, arch, collapse of arch,
aeration and eccentric discharge.
(2) Design for group silos should take into account the varying operational scenarios.
e.g., some silos are full, others are empty, and some are in the process of being filled
or emptied.
(3) In the case that the characteristics of the stored material may be subject to change
during storage, the design should consider the most severe load conditions.
(4) The designer may use an alternative method for calculating pressures if it can beverified by experiments etc.
5.2.2.2

Pressures and friction forces of static stored materials
Static pressures exerted by stored material at rest should be calculated using thefollowing methods.

(1)

The vertical pressure, Pv, at depth x below the surface of the stored material iscalculated using equation (5.1).
(5.1)

Pv =

− µ f ⋅ K s ⋅ x / rw γ ⋅ rw ⋅ 1 − eµ f ⋅ Ks in which Ks and rw are obtained from equations (5.2) and (5.3).

(2)

(5.2)

Ks =

1 − sin φi
1 + sin φi

(5.3)

rw

=

A/lc

The horizontal pressure, Ph, of a unit area at depth x from the free surface is obtained
from equation (5.4).
(5.4)

(3)

where : K s ≧ 0.3

Ph = K s ⋅ Pv

The friction force, Pf, of a unit area of silo wall at depth x from the assumed freesurface should be computed from equation (5.5).
(5.5)

Chapter 5

Pf = µ f ⋅ Ph

54

(4)

The static unit pressure Pa normal to the surface inclined at angle α to thehorizontal at depth x below the surface of stored material is obtained form equation(5.6).
Pa = Ph ⋅ sin 2 α + Pν ⋅ cos 2 α(5.6)
5.2.2.3(1)

Design pressures and forces exerted by stored materials
Design pressures should be obtained by multiplying the static pressures by the
overpressure correction factor Cd or impact factor Ci in accordance with theappropriate equations (5.7), (5.8), or (5.9).

dPν = Ci ⋅ Pν

(5.7)

for silo wall.

dPν = Ci ⋅ Pν  for hopper, whichever is greater.
dPν = C d ⋅ Pν 
(5.8)

dPh = Cd ⋅ Ph

(5.9)

dPa = Ci ⋅ Pa 
 whichever is greater.
dPa = C d ⋅ Pa 

(2)

The impact factor, Ci, should be used to consider the material properties, thecharging method and the charging rates, and should be taken a value between the 1.0
and 2.0.

(3)

The minimum required values of overpressure factor Cd are given in Figure 5.2.

Cdfor hm/d≧4.0

hm-h0

hm

hm

hm-h0

Cd=2.0for 4.0>hm/d>1.0
Cd=4/3+hm/6d

h0

h0

Cd
1.2 × C d
1.0

for 1.0≧hm/d
Cd=1.5h0 : hopper heightdFig.5.2 Minimum Required Values of Overpressure Factor C d(4)

For small scale and 1.5≧hm/ d silos, minimum required factor may be Cd =1.0,
however Pv is taken as Pν = γ ⋅ x .

(5) The long - term local pressure shown below should be taken into account, as well as
above discussed Ci and Cd , to consider the increase and decrease of pressure byarching and collapse of the arch etc. in stored materials.

Chapter 5551) The design pressure Cd ⋅ Ph is multiplied by CL, which ranges from 0.15 to 0.4
and depends on eccentricity of the discharge outlet, and this local pressure dPL (=C L ⋅ C d ⋅ Ph ) is applied to the partial area of the silo wall of 0.1d width. A set ofpressures, compressive directions or tensile ones both, is applied to any point symmetric locations of silo wall, as shown in Figure 5.4. The eccentricity “e” isdefined in Figure 5.3.
2) CL=0.15 when e=0, and CL=0.4 when e=0.5d, the case where the outlet is just on
the wall. For other values of “e”, CL is given by CL=0.15+0.5 (e/d).
3) For cases where “e” is less than or equal to 0.1d, the local pressure for the upper
part of the silo wall may be reduced. The CL is given by CL=0.15+0.5(e/d) if theheight of the location from the bottom end of the silo walls is less than 1.5d, andtakes zero at the top. Between the 1.5d height and the top, a linear function of height
gives the CL value. See equation (5.10) and Figure 5.5.

discharge outlet

dPL
tensileecompressive
compressive

dPL
0.1d
tensilee : eccentricity

Fig.5.4 Local Pressure

Fig.5.3 Definition of Eccentricity(5.10)(6)

dPL = C L ⋅ C d ⋅ Ph C L = 0.15 + 0.5(e / d )Vertical long-term design force d:m per unit length of silo wall horizontal section, at
depth x from free surface, should be computed as shown in equation (5.11).
(5.11)

d: m = ( γ ⋅ x − Pν ) ⋅ rw ⋅ C f

where Cf is taken both 1.0 and more than 1.5, because either case may generate thelarger stress.
(7)

Vertical short-term design force d:m by the seismic force should be calculated byequation.5.12(5.12)

d:m = ( γ ⋅ x − Pν ) ⋅ rw

(8) The axial force in the direction of generator on the hopper wall horizontal section perunit length, : φ , and the tensile force in the circumferential direction, :θ , arecalculated using equations (5.13) and (5.14) respectively. See Figure 5.6.
(W + Ws ) dPv ⋅ d '(5.13):φ = h+lc ⋅ sin α
4 ⋅ sin α

Chapter 556

(5.14):θ =dPa ⋅ d '2 ⋅ sin α

Where :

Wh

weight of the contents of the hopper beneath the section (N)Wsweight of the hopper beneath the section (N)

lc
d'
α
dPv

perimeter of the hopper at the section level (mm)diameter of the hopper at the section level (mm)slope of the hopper ( o )design vertical pressure of unit area at the section level (N/mm2)dPa

design normal pressure of unit area to hopper wall at the section level(N/mm2)For design of the connection of the hopper and the silo wall, however, an additionalcase where : φ is assumed to be zero should be considered.
e/d＞0.1
e/d≦0.1

1.5d

00.15～0.4

dFig.5.5 Local Pressure Factor CL
d'

dPv

dPa

section level
generator direction

generator directionαcircumferential direction
Fig.5.6 Force Acting on the Hopper
5.2.2.4

Mixing by Stirring and Rotation
The effect of mixing materials by stirring and rotation, in which stored materials aredischarged to other silos and simultaneously charged from other silos, should beconsidered for both dynamic overpressure and impact pressure.

Chapter 557

5.2.2.5

Eccentric discharge
The design should take into account the effect of asymmetric flow, the local pressure
given by equation (5.10), caused by concentric or eccentric discharge openings. If asilo has mobile discharge openings, the maximum value of eccentricity should be
adopted.

5.2.2.6

Pressure increase due to aeration
In the case of aeration in a silo and pneumatic conveyance, the possibility ofpressure increase should be considered.

5.2.2.7

Flushing phenomenon
The possibility of flushing phenomenon should be considered where there is thepossibility that the silos could be mis-operated, or where it is likely that material
fluidity could be a contributing factor in any seismic event.

5.2.2.8

Thermal stresses
Consideration should be given to thermal stresses generated in the walls of silos inwhich hot materials are stored.

------------------------------------------------------------------------------------------------------------Commentary:
5.2.2.2

Pressures and friction forces of static stored materials
Foreign standards5.1)~5.3) may be referred to for better understanding of this chapter. Janssen’s
equation5.4) and Rankine’s earth pressure (equations (5.1) and (5.2)) are used for static pressures. To
obtain design pressures, overpressure factor Cd , impact factor Ci and local pressure dPL should beconsidered.
5.2.2.3
Design pressures and forces exerted by stored materials
1) Local pressure
Although large pressures of stored materials to silo walls had been observed in many cases, theactual effect of them on the wall stress (strain) was regarded as unclear. Recently simultaneous
measurements of wall strain and pressure were made to make clear the effect of the pressure to the strain
as well as to discriminate between bending and membrane strains 5.5), 5.6). Simulation analyses of themeasured results assuming local pressure were carried out 5.6). Dependency of the magnitude andfluctuation of the local pressure on the eccentricity of the discharge outlet was also studied 5.5). The local
pressure given by the equation (5.10) is determined with due regard to above references and other design
and accident examples.
2) Local stress coefficient
Based on parametric studies and for approximate estimations, equations below are made to
calculate bending stresses.
2-1) The maximum bending moment for a simple ring loaded with a pair of normal local force at
opposite sides is given by the equation below.
M 0 = 0.318 ⋅ 0.1d ⋅ dPL ⋅ ( d / 2)
(5.2.1)
Where M 0 is the bending moment per unit width of the ring calculated using simple ring model.
If a 3D cylindrical model with high rigidity bottom and top boundaries is considered instead of
the simple ring, the bending moment will be smaller.
2-2) For steel silos, approximate estimation of the maximum bending stress (moment) is given bymultiplying the coefficient “ k ” below into the “ M 0 ” above in some cases defined below. High rigidity
boundary cases, with stiffening beams at the boundary, roofs, hoppers e.g., are assumed.

Chapter 5

58

The local stress coefficient “ k ” is applicable when hm / d ≦ 5 , d ≦10m and t ≦ 20mm
For each plate thickness of the silo wall, the location (height) where the bending stress is maximum
should be checked.
iii) The local stress coefficient “ k ” is given by the equation below.
k = 0.1 + 0.05( hm / d − 2)
(5.2.2)2-3) The bending moment M b is given by the equation below.
M b = kM 0
(5.2.3)where M b is the bending moment per unit width of the silo wall for which the effect of rigid boundaries
is approximately considered.
2-4) For reinforced concrete silos, another local stress coefficient is defined below, similarly to steelsilos. High rigidity boundary cases are assumed.
(5.2.4)
k = 0.25 + 0.2( hm / d − 1)
In case hm / d is smaller than 1, k is fixed at 0.25. In many cases, the rigidity of the top end is
not sufficient to ensure the accuracy of the equation (5.2.4). Therefore, the bending moment for the silowall near the top end is modified. From the top surface of the stored material to 0.2 hm depth, thebending moment at the 0.2 hm depth is fixed for this region.
3) In silos with mobile discharge openings or a large number of discharge outlets, the bottom region isflat compared to single hopper cases. Therefore, large pressure to the bottom region is probable andshould be designed with due regard.
4) Large pressure and stress are expected at the connection of the silo wall and the hopper because of thehopper weight, the pressure and the frictional force on the hopper, etc. Sufficient margin is necessary in
structural design.
5) Equations (5.13) and (5.14) are derived from the membrane shell theory 5.7), and the equation (5.13)
can also be derived from equilibrium of forces. In the equation (5.14), only the normal pressure isconsidered and the effect of the frictional force is neglected.

i)
ii)

5.2.3

Seismic Design

5.2.3.1

The importance factor is determined in accordance with the classification of theseismic design I, II and IV given in Table 5.1 below.
Table 5.1 Classification of Seismic Design and Essential Facility Factor
Seismic Design
Classification

Silos

Importance
Factor I

I

Small silos - the height of which is
below 8.0m

0.6 and over

II

Other
silos
not
Classification I and IV

0.8 and overIV

Silos used for the storage of hazardous
or explosive materials

includedin

1.2 and over

5.2.3.2

The design seismic loads are obtained from the modified seismic coefficient analysis
or the modal analysis which are shown in 3.6.1.2 and 3.6.1.3.

5.2.3.3

In case that the first natural period, T, is unknown, T should be taken as 0.6s for the
purposes of the silo design. In this case, the modified seismic coefficient analysis
should be used.

Chapter 559

5.2.3.4

Using the modified seismic coefficient analysis, the design yield shear force of thesilo at the ground base of the silo Qd should be taken from equations (3.1).

Notations:design yield shear force coefficient at the base of the silo
design weight imposed on the base of silo
importance factor, given in Table 5.1
a value obtained from the plastic deformation property of the structure,
given in Table 5.2
a value obtained from the radiation damping of the silo basement, given
in Table 5.3

C
W
I
Dη

Dh

Table 5.2 Dη Value
DηStructure
Reinforced concrete

0.45

Steel or aluminium alloy

0.5 ~ 0.7

Table 5.3 Dh Value
Short length of silo foundation

5.2.3.5

Area of silo foundation≦25m≦2000m2

25m ~ 50m

2000m2 ~ 4000m2≧50m≧4000m2Dh

1.0
1.0 ~ 0.9
0.9

Assessment of Seismic Capacity
The seismic capacity of the silo design should be assessed using equation.(3.10)
given in paragraph 3.6.1.8. The acceptable ratios of the horizontal load-carrying
capacity of the silo structure to the short-term allowable yield strength are:5.2.3.6

Reinforced concrete silos:

B =1.5.

Metal silos:

B =1.0.

When calculating seismic forces on silos by the modified seismic coefficient ormodal analysis, the lateral forces can be reduced due to, other than the radiation

Chapter 560

damping of silo structures, energy loss caused by the internal friction between thestored materials or between the materials and the silo walls. The weight of the stored
material can be reduced for this purpose. Similarly, when calculating seismic forces
for group silos which are not always fully loaded, a reduction of material weight canbe considered. But in neither case should the reduced weight of stored material beless than 80% of the maximum stored material weight.
5.2.3.7

Dynamic analysis considering the deformation characteristics of the stored material
In case large ground motion of the design earthquake is assumed, or for some kindsof stored material and silo configuration, effect of the deformation characteristics ofthe stored material is very large. In such cases, dynamic response analysis forseismic design has an advantage because the energy absorption by the stored
material friction and slip can be taken into consideration when the silo response andload are estimated.

5.2.3.8

Group silo
Due consideration should be given to group silo design because loads and stresses
are different from single silo cases.

5.3

Reinforced Concrete Silos
The design of reinforced concrete silos should conform to the requirements of the“Standard for Structural Calculation of Reinforced Concrete Structures (2010) Architectural Institute of Japan”, together with the provisions of the following
paragraphs.

5.3.1

Both the membrane stress and the out-of-plane stress should be given due
consideration in the design of silo walls and hoppers. Incidental stress and stresstransmission at connection of shells and plates should be taken into consideration
especially of group silos and silos with unusual configurations.

5.3.2

Design allowable crack width should be appropriately determined considering
necessary margin against deterioration of reinforcing bars and leakage.

5.3.3

The cross-sectional shape and the quality of the concrete parts and the arrangementsof steel bars should be carefully considered when using the slip form construction
method.

5.4

Steel Silos
The design of steel silos should conform to section 2.2, the “Design Standard for
Steel Structure – Based on Allowable Stress Concept - (2005) - Architectural
Institute of Japan” and the recommendations given below.

5.4.1

Design of Cylindrical Walls
The guidelines given in section 3.7.1 should be followed when examining thebuckling stress of cylindrical walls.

Chapter 5

61

The calculation of stresses on the parts of the silo during earthquakes should bebased on the requirements of section 3.7.4.
The calculation of short-term stresses on parts caused by other than earthquakes
should be based on section 3.7.3.
5.4.2

Design of Hopper Walls
The axial stress on the section of the hopper wall in the direction of generator, σφ ,and the tensile stress in the circumferential direction, σθ , are calculated using
equations (5.15) and (5.16) respectively.
(5.15)

σφ =(5.16)σθ =(Wh + Ws ) +lc ⋅ t2 ⋅ sin α

dPv ⋅ d ′≤ ft
4 ⋅ t 2 ⋅ sin αdPa ⋅ d ′≤ ft
2 ⋅ t2 ⋅ sin α

where:
Wh
Ws
t2
lc
d’
α
dPv
dPa
ft

weight of the contents of the hopper beneath the section (N)weight of the hopper beneath the section (N)wall thickness of the hopper (mm)perimeter of the hopper at the section level (mm)diameter of the hopper wall at the section level (mm)inclination of the hopper wall (degree)
vertical design pressure per unit area at the section level (N/mm2)
perpendicular design pressure per unit area on the hopper wall at thesection level (N/mm2)
allowable tensile stress (N/mm2)

5.4.3 Design of Supports
Local bending stresses should be considered in design procedure that may act on theconnections between the hopper or machines and their supports or other similar parts.
References:
[5.1]
[5.2]

[5.3]
[5.4]
[5.5]

[5.6]

International Standard : ISO 11697, Bases for Design of Structures – Loads due to Bulk Materials,
1995.6
American Concrete Institute : Standard Practice for Design and Construction of Concrete Silos andStacking Tubes for Storing Granular Materials (ACI 313-97) and Commentary—ACI 313R-97,
1997.1
Australian Standard : Loads on Bulk Solids Containers (AS 3774-1996), 1996.10
Janssen, H.A. : Versuche U&&ber Getreidedruck in Silozellen : VID Zeitschrift (Dusseldorf) V. 39, pp.
1047～1049, 1885.8:aito, Y., Miyazumi, K., Ogawa, K., :akai, :. : Overpressure Factor Based on a Measurement of a
Large - scale Coal Silo during Filling and Discharge, Proceedings of the 8th International
Conference on Bulk Materials Storage, Handling and Transportation (ICBMH’04), pp. 429 – 435,
2004.7
Ishikawa, T., Yoshida, J. : Load Evaluation Based on the Stress Measurement in a Huge Coal Silo (inJapanese), Journal of Structural and Construction Engineering (Transactions of AIJ), :o.560, pp.
221 – 228, 2002.10

Chapter 562[5.7][5.8][5.9][5.10]

S.P. Timoshenko and S.W. Kreiger : Theory of Plates and Shells, McGraw-Hill, 1959
G.E. Blight : A Comparison of Measured Pressures in Silos with Code Recommendations, Bulk
Solids Handling Vol. 8, :o. 2, pp. 145～153, 1988.4
J.E. Sadler, F. T. Johnston, M. H. Mahmoud : Designing Silo Walls for Flow Patterns, ACI
Structral Journal, March-April 1995, pp. 219～228 (Title :o. 92-S21), 1995
A.W. Jenike: Storage and Flow of Solid, Bulletin :o. 123 of the Utah Engineering Experiment
Station, University of Utah, 1964.11

Chapter 5636.

Seismic Design of Supporting Structures for Spherical Tanks

6.1

Scope

6.1.1

This chapter is applicable to the seismic design of supporting structures for sphericalsteel tanks. Other general requirements for designing structures should be inaccordance with the applicable codes and standards.

6.1.2

In general, the supporting structures should be composed of steel pipe columns withbracings made of steel pipes or solid round bars. However, other steel sections may
be used and designed in accordance with this chapter.

6.1.3

Notationsexternal diameter of the supporting pipe column (mm)weight ratio of the effective mass of the tank contents
short-term allowable lateral force of the supporting structure (kN)short-term allowable compressive strength of bracings (kN)short-term allowable tensile strength of the bracings (kN)wall thickness of the supporting pipe column (mm)usable liquid volume (m3)volume of the spherical tank (m3)design weight imposed on the base of structure (kN)dead weight of the structure (kN)weight of the stored contents (kN)δmhorizontal displacement of the column head when a supporting column
reaches the short-term allowable strength (m)δyhorizontal displacement of the column head when a bracing reaches theshort-term allowable strength (m)ξratio of the usable liquid volume to the volume of spherical tank = VI /VDηmean cumulative plastic deformation rate of the supporting structure------------------------------------------------------------------------------------------------------------d
f(ξ)
Fy
:C
:T
t1
VI
VD
W
WD
WI

Commentary:bracing
supporting pipe column

Fig.6.1.1 General View of Spherical Tank

Chapter 664Steel spherical tanks are usually supported on steel pipe columns where steel pipe bracings
reinforce the column structures as illustrated in Figure 6.1.1. The intention of this chapter isto provide seismic design recommendations for this type of structure.

6.2

Material Properties and Allowable Stresses

6.2.1

Materials for the structures should be in accordance with section 2.2 of thisrecommendation.

6.2.2

Table 6.1 gives the physical properties of steel materials.
Table 6.1 Physical Properties of Steel Materials
Materials

Young’s Modulus(N/mm2)

Shear Modulus(N/mm2)

Poisson’s Ratio

Expansion
Coefficient(1/°C)

Steel Plates, Pipeand Cast Steel

2.05 x 105

7.9 x 104

0.3

12 x 10-6

6.2.3

Allowable stresses should be determined as follows:(1) Allowable stresses for steel materials should be determined by use of the F value andin accordance with Chapter 5, Design Standard for Steel Structures – Based on
Allowable Stress Concept – (2005), Architectural Institute of Japan.
(2) F values of steel plates and pipe should be as shown in Table 6.2. F values for hardsteel wires should be taken to be the value either 70% of tensile strength or yieldstrength whichever is lesser.
Table 6.2 F-Values

Type of steel

Grade≦40mm thickRolledSteelsfor
General Structure (JISG 3101), Rolled Steels
for Welded Structure;(JIS G 3106), Rolled
Steel for Building
Structure;
(JIS
G
3136)

SS400, SM400, SN400

235

F(N/mm2)> 40mm thick > 75mm thickand ≦ 75mm and ≦ 100mm
thick
thick
215

SM490

315

295

SN490

325

295

SM490Y

343

335

325

SM520

355

335

325

SM570

399

339

STK400235

215

STK490315

295

Carbon Steel Tubes for
General
Structural
Purposes; (JIS G 3444)

Chapter 6

65

Carbon Steel Square
Pipes for General
Structural
Purposes;(JIS G 3466)

STKR400

245

215

STKR490325

295> 6mm thick and≦ > 50mm thick and > 100mm thickand ≦200mm≦100mm thick
50mm thick
thick
Steel
PlatesforPressure Vessels forIntermediate
Temperature Service;(JIS G 3115)
Mn-Mo and Mn-MoNickel Alloy Steel
Plates Quenched andTempered for Pressure
Vessels; (JIS G 3120)
Carbon Steel Plates for
Pressure Vessels for
Low
Temperature
Service; (JIS G 3126)

SPV235SPV315SPV355SPV450SPV490

235
315
355
399
427

215
295
335
399
427

SQV1A, 2A, 3A

345

SQV1B, 2B, 3B

434≦40mm thick
235

SLA235SLA325
SLA365

195
275
315
399
427

> 40mm thick
215
308
343

6.3

Seismic Design

6.3.1

Seismic Loads

6.3.1.1

Seismic loads should be determined in accordance with section 3.6 of chapter 3.

6.3.1.2

Design seismic loads should be determined by the modified seismic coefficientmethod in par. 3.6.1.2.

(1)

Design yield shear force, Qd, should be calculated by the equations (6.1) and (6.2)
below.
(6.1)

Qd = C ⋅ W(6.2)

C = Z s ⋅ I ⋅ Ds ⋅S a1
gWhere : C ≧ 0.3Z s ⋅ I
Notations:Qd
design yield shear force (kN)C
design yield shear force coefficientWdesign load imposed on the base of the structure, (kN), whichis equal to the sum of dead weight of structure and weight of
impulsive mass of liquid content = WD + f(ξ) WI (kN)WD
dead weight of the structure (kN)

Chapter 666

weight of the stored contents (kN)weight ratio of the impulsive mass of the tank contents – refer
to Fig. 6.1 below.
ratio of the usable liquid volume to the volume of sphericaltank i.e., VI /VD
usable liquid volume (m3)volume of the spherical tank (m3)design acceleration response spectrum, (m/s2), corresponding
to the first natural period, given in par. 3.6.1.6
acceleration of gravity (9.8m/s2)
seismic zone factor, which is the value as per the BuildingStandards Law, or the local government
importance factor, given in par. 3.6.1.7
structural characteristic coefficient, which is the value
specified in par. 6.3.1.2.(2)

WI
f(ξ)

ξ
VI
VD
Sa1
g
Zs
I
Ds

Fig.6.1

Weight Rtio of the Impulsive Mass of the Tank Contents(2) The structural characteristic coefficient, Ds, should be calculated using equations (6.3)
and (6.4) below.

(6.3)

Ds =

(6.4)

η=

1
1 + 4 aηδm − δ yδywhere:aη ≤ 3.0

Notations:ηmean cumulative plastic deformation ratio of the supporting
structureacorrection factor for η considered the reduction of thestrength of the bracings
a=1.0 for : T ≧ 2 : C
a=0.75 for : T＜2 : C:CChapter 6

short-term allowable compressive strength of bracings (kN)

67

:T

δm

short-term allowable tensile strength of the bracings (kN)horizontal displacement of the column head when aδysupporting column reaches the short-term allowable strength(m)horizontal displacement of the column head when a bracing
reaches the short-term allowable strength (m)6.3.2

Evaluation for Seismic Resistance of Supporting Structures
Evaluation for seismic resistance of supporting structures should be calculated byusing equation (6.5).
(6.5)

Fy ≥ Qd
where:
Fy

6.3.3

short-term allowable lateral force of the supporting structures(kN), provided that the lateral force of the compressive
bracings should be ignored where :T ≥ 2:C

Design of Supporting Structure Components
Design of supporting structure components and bases should be performed inaccordance with “Design Standard for Steel Structures – Based on Allowable Stress
Concept – (2005) ” and “Recommendations for Design of Building Foundations(2001)” both published by Architectural Institute of Japan. In addition, the following
requirements should be satisfied.
(1) The slenderness ratio of the supporting columns, should be, ≦ 60. The ratio of
diameter to thickness of the steel supporting columns is determined in accordance withequation (6.6) using F value (N/mm2).
(6.6)

d1 19,600≤
t1
F
where:
d1
t1

external diameter of the pipe supporting column (mm)thickness of the pipe supporting column (mm).

(2) The ratio of diameter to thickness of the steel pipe bracings, should be, ≤ 40.
(3) The strength of the joints in the bracings should be greater than the tensile yield forceof the bracing members.
(4) Bracings should be installed without slackness.
-------------------------------------------------------------------------------------------------------------Commentary:
As the structure discussed in this chapter is considered to be the vibration system of a mass, theseismic force is obtained by the modified seismic coefficient method given in 3.6.1.2. The design

Chapter 6

68

weight imposed on the base of the structure W is calculated as the sum of WD , and f(ξ) WI . The
weight ratio of the impulsive mass of the tank contents, f(ξ ), is determined with reference toexperimental data of the vibration of a real scale structure [6.1] [6.2].
The structural characteristic coefficient, Ds , is obtained from equation (6.3) which is based on themean cumulative plastic deformation ratio of the supporting structure, η, (see Commentary 3.6.1),
which is obtained from equation (6.4). This equation assumes that δm is the horizontal
displacement of the column head when a supporting column reaches the short-term allowablestrength, that δy is the horizontal displacement of the column head when a bracing reaches theshort-term allowable strength, and δm - δy is the mean cumulative plastic deformation. When theallowable strength of the bracing is governed by a buckling, i.e., when :T < 2:C , the energy
absorption capacity (deformation capacity) is lowered due to the lower strength of the structure
after the structure becomes plastic. Considering this lowering, equation (6.3) adopts “a” value forlowering the deformation capacity by 25%. Seismic design is evaluated by reference to equation(6.5) assuming that the short-term allowable strength Fy is considered to be the lateral retaining
strength [6.3].

1.

The calculation method of η and the design method to obtain the plastic deformation
capacity of bracings is given in the following paragraphs.

1.1 Elastic analysis of the structure (when :T < 2:C)
When the horizontal section of the supporting structure is a regular polygon with n equal sides asshown in Figure 6.3.1, the maximum stress takes place in the side which is parallel (or is theclosest to parallel) to the direction of the force. The shearing force in the side is obtained fromequation (6.3.1).

Fig.6.3.1 Modelling(6.3.1)Q=

2Fnwhere:
F
is the lateral force imposed on the supporting structure,nis the number of columns = the number of the sides.

Chapter 6

69

Although it should be sufficient that an analysis of the side which receives the in-plane shear force
is carried out, and as such can be considered to represent the total design, it should be noted thatthe adjacent bracings to the sides at point “b” have both a tensile force and a compressive force
respectively. The compressive strength of the bracing is generally smaller than the short-termallowable tensile strength, and the compressive force due to the adjacent bracing is generally
smaller than the force in the bracing being analysed as the compressive force is not in the analysed
plane. To simplify the analysis, the design force acting at point “b” is considered as being twice
the horizontal component of the tensile bracing force at point “b”
The analysis should proceed under the following conditions:i)

In Figure 6.1.1, h = ho – lr/2 as demonstrated by experimental data and where lr is obtained
from equation (6.3.2).

(6.3.2)Dd
2where:
D
diameter of spherical tankddiameter of supporting column

lr =

ii)

An elastic deformation and a local deformation of the sphere can be neglected as the rigidityof the sphere is relatively large.

iii)

To simplify the analysis, the conservative assumption is made that the sphere-to-column joint
is considered to be a fixed joint in the plane tangent to the sphere, and a pinned joint in theplane perpendicular to the tangential plane and including the column axis.

iv)

The column foot is considered to be a pinned joint.

v)

Column sections are considered uniform along their entire length.

vi)

Column deformation in the axial direction is neglected.

The column head “a”, shown in Figure 6.3.2 (a), displaces with δ, [as shown in equation (6.3.3)]
due to the lateral force Q.
(6.3.3)δ = δ1 + δ2 + δ3
where:

δ1
δ2
δ3

Chapter 6

displacement due to an elongation of the bracing
column displacement at b-point due to the horizontal component 2S ofthe bracing
relative displacement between “a”point and b-point when the joint
translations angle of the column is assumed

70

Fig.6.3.2 Stresses in Supporting Column
From the elasticity theory of bars, δ, δ1, and δ2 may be calculated in accordance with equations(6.3.4), (6.3.5), and (6.3.6).
2h 3(6.3.4)δ=(6.3.5)

δ1 =Sh1Ab E cos 2 θ ⋅ sin θ(6.3.6)

δ2 =Sh1 h2 (4h1 + 3h2 )6 EIh 33

3h h1 − h122

(δ1 + δ 2 )3Equation (6.3.7) is obtained from the above three equations:(6.3.7)

 h12 h2 3 (4h1 + 3h2 )

h1
+
S
3 
3
2
2
Ab E cos θ ⋅ sin θ 
6 EIh
(3h h1 − h1 ) 
where:Ab
cross sectional area of the bracingI
moment of inertia of the supporting column
E
Young’s modulus of the steel
2h 3

δ=

Figure 6.3.2 (b) shows a stress diagram of a supporting column with a force of 2S exerted by thebracing when the column head “a” does not displace. The reacting forces and moments at points“a”, “b”, and “c” are obtained from the following equation (6.3.8).
(6.3.8)
2

Rc =

Sh2
(3h1 + 2h2 )
h3

Ra =

Sh1
2
2
(2h1 + 6h1h2 + 3h2 )h3

M a1 =

Sh1h2(2h1 + h2 )h22

M b1 = −Sh1h2(3h1 + 2h2 )h3

M a > M b when h2 < 2h1

Chapter 6

71















Figure 6.3.2 (c) shows a stress diagram of the supporting column when the column head “a”
displaces by δ. The reacting forces and moments at points, “a”, “b”, and “c”, are obtained
from the following equation (6.3.9).

(6.3.9)

R′ =

M a2
M
R′ = a 2hh

M a2 =
M b2 =

 h12 h23 (4h1 + 3h2 )
3EI
1
6 Ihh1
+
δ=
S
2
3 
2
2
2
h
h
Ab cos θ ⋅ sin θ 
(3h h1 − h1 )  h12 h23 ( 4h1 + 3h2 )

1
(3h

2

2 
− h1 ) 

h

3

+


6 Ih1
S
2
Ab cos θ ⋅ sin θ 













The lateral force Q is calculated using equation (6.3.10).
(6.3.10)

Q = Ra + R 'The reacting force at the supporting column foot, R, is calculated using equation (6.3.11).
(6.3.11)

R = R '− Rc

The bending moments at points a, and b, Ma and Mb, are given in equations (6.3.12), and (6..3.13)
respectively.

1.2(6.3.12)Ma = Ma1 + Ma2(6.3.13)

Mb = Mb1 + Mb2

Elastic analysis of the structure (when : T ≧ 2 : C )In this case, the compressive force due to the bracing is neglected. Therefore, the horizontal force
2S acting at point “b” shown in Figure 6.3.2 is converted to S, and equations (6.3.14) ~ (6.3.17)
are given instead of equations (6.3.6) ~ (6.3.9). Equations (6.3.1) ~ (6.3.5) and equations (6.3.10)~ (6.3.13) have no modification except the definition of δ 2 , the displacement due to S in this caseinstead of 2S.
(6.3.14)(6.3.15)(6.3.16)

Chapter 6

Sh12 h23 (4h1 + 3h2 )12 EIh 3
 h12 h23 (4h1 + 3h2 )

h1
2h 3
+
δ=
S
2
3 
3
2
(3h h1 − h1 ) 
12 EIhAb E cos θ ⋅ sin θ 
Sh 2
Rc = 23 (3h1 + 2h2 )

2h

Sh
Ra = 13 (2h12 + 6h1 h2 + 3h22 ) 

2h

Sh1h2

M a1 =(2h1 + h2 )2
2hSh1h22
M b1 = −(3h1 + 2h2 ) 
3
2h


δ2 =

72

M a2 =
(6.3.17)

 h12 h23 (4h1 + 3h2 )
3EI
1
6 Ihh1
δ=
+
S
2
2
3 
2
2
h
(3h h1 − h1 ) 2h
Ab cos θ ⋅ sin θ 
 h12 h23 (4h1 + 3h2 )

1
6 Ih1
M b2 =
+
S
2
2 
3
2
(3h − h1 ) 2h
Ab cos θ ⋅ sin θ M
R' = a 2
h












1.3 Yield strength Fy and yield deformation δ y of supporting structures
When : T < 2 : C , the yield strength of supporting structures should be equal to the strength ofwhich the maximum stress element of a bracing in compression becomes yielding, i.e., S = Sy =1.5Ab fc cosθ (where fc is the allowable compressive stress) is input into equation Ra of (6.3.8), andin the equation R' of equation (6.3.9). Qy is obtained as Q = Qy=Ra+R’ in (6.3.10) and the yieldstrength (short-term allowable strength) of the supporting structure Fy is obtained from (6.3.1).
The yield deformation δ y is similarly obtained by inputting S = Sy = 1.5 Ab fc cosθ into (6.3.7).
On the other hand, when :T≧2:C, the yield strength of supporting structures should be equal to
the strength of which the maximum stress element of a bracing in tension becomes yielding, i. e.,
S=Sy=1.5 Ab ft cosθ (where ft is the allowable tensile stress) is input into equation Ra of (6.3.16),
and in the equation R’ of (6.3.17). Qy is obtained as Q=Qy=Ra+R’ in (6.3.10) and the yield strength(short- term allowable strength) of the supporting structure Fy is obtained from (6.3.1).
The yield deformation δ y is similarly obtained by inputting S=Sy=1.5 Ab ft cosθ into (6.3.15).
The spring constant k can be calculated from equation (6.3.18) as follows:(6.3.18)

1.4

k=

Fy

δy

Critical deformation of a frame δ mWhen the lateral displacement of the column head increases after the bracing reaches its elastic
limit, points “a” and “b” reach their yield point. In this situation, Ma1 and Mb1 do not increase as
the bracing is already in yielding condition, but Ma2 and Mb2 increase so that eventually the column
reaches its yield point. It is noted that point “a” at the column head reaches its yield point first as
long as h2 is relatively short as shown in equation (6.3.8).
Because the rigidity of the sphere-to-column joint, ±lr / 2 distance region from point “a”, is large,
the bending moment diagram of the column at this joint is not linear function of the height asshown in Figure 6.3.3. Figure 6.3.3 shows that the critical bending moment, eMa, should be
calculated by the Equation (6.3.19) at a distance of lr/4 below point “a”.

Chapter 673

Fig.6.3.3 Modelling of the Top Part of Supporting Column(6.3.19)elQ
Ma = Ma − r
4Values for the lateral force Q, and the bending moment eMa are obtained from equations (6.3.10),
and (6.3.12) respectively.
The yielding moment of the column My is calculated using equation (6.3.20) as follows:(6.3.20):
M y = Z σ y − A: = (WD + WI ) / n:otations:Z
section modulus of the column:
axial compressive force in the columnA
sectional area of the column
σyshort-term allowable tensile stress (yield stress)nnumber of columns
The surplus in the critical bending at the −lr / 4 distant location from point “a”, ∆M, is obtained
from equation (6.3.21) as follows:(6.3.21)∆M = My - (eMa)ywhere: (eMa)y is equal to eMa when the supporting structure reaches the yield strength Fy
The increment in the virtual bending moment at point “a” is shown in equation (6.3.22).
h  h − lr 4This equation is necessary to obtain the increment in displacement at point “a”.
(6.3.22)

Δ M a = ΔM 

The increment in displacement, ∆δ, at point “a” is given in equation (6.3.23).
(6.3.23)Δδ = δm -δy =∆M a h 2
3EI

η is obtained from (6.4) by use of (6.3.23) andδy , as follows.
(6.3.24)

Chapter 6

η = (δm -δy ) /δy74

:ote that buckling is not considered in equation (6.3.20) in the evaluation of a supporting column.
The overall length of the supporting column, h, should be considered in assessing out-of-plane
buckling as no restraint is provided at joint “b” of the bracings in terms of out-of-plane
deformation. On the other hand, in-plane buckling between points “a” and “b” need not be
considered as the length is relatively short, and if the effect of buckling is considered between
points “b” and “c”, the evaluation at point “a” by (6.3.20) governs. Out-of-plane buckling may beconsidered as Euler Buckling, and as such a suitable evaluation method is given in equation(6.3.25).
(6.3.25)

σc
fc

+ c

σb
ft

≤1

However, steel pipe columns need not checking of flexural-torsional buckling, no out-of-planesecondary moments due to deflection and axial forces. Therefore, equation (6.3.25) has little
meaning for pipe columns. As the equation (6.3.25) is based on a pipe column, it is assumed that fb= ft .
Taking the above observations into consideration, equation (6.3.20) is adopted as an acceptableevaluation method on condition that the slenderness ratio is less than 60 and on the assumptionthat the buckling length is h. The effective slenderness ratio is approximately 42 (0.7×60) for outof-plane buckling and the column will never buckle before yielding at the extremity of point “a”.
Although the above summary assumes that the bracing yields first, such design that point “a” ofthe column yields first is allowable. However, this is not a practical solution in anti-seismic design,
as the energy absorption by the plastic elongation of the bracing cannot be utilised. As the
supporting ability against lateral seismic load decreases when the bracing which is in compression
stress buckles, the yield strength of the frame will be reduced, and the restoring force
characteristics after yielding will become inferior. However, experimental data[6.4] show that insupporting structures which use pipe bracings, the lateral retaining strength is reduced byapproximately 25% of the yield strength even when a comparatively large plastic deformation takes
place in the bracing.
Typical examples of the joint design between the pipe column and bracing are shown in Figures
6.3.4 to 6.3.7 [6.5].
In Figure 6.3.4 the design force at the joint should be determined by the 130% of the yield strengthof the bracing [6.6].
In Figure 6.3.5(b) the short-term allowable force Pa1 at the bracing cross is obtained in equation(6.3.26);2(6.3.26)

Pa1 =2

29t1 F 29t1 F(:)=
sin ϕ1
sin 2θ

Where Pa1 > AρF (where Aρis the sectional area of pipe bracing), equation (6.3.27) is used.
(6.3.27)

t1 >Aρ sin ϕ1
5.5

(mm)

In case the requirements of the above equation cannot be satisfied, either the thickness of the pipe
wall should be increased at the joint, refer to Figure 6.3.6(a), or the joint should be reinforced witha rib as indicated in Figure 6.3.6(b).
In Figure 6.3.5(c) the short-term allowable force, Pa2, at the joint between the column and bracing
is obtained from equation (6.3.28) below:

Chapter 675

(6.3.28)

 d1   f (n ) 2
 ⋅ t2 F (:)
 d 2   sin ϕ 22Pa2 = 5.351 + 4.6

where:

sin ϕ 2 = cosθ
(when n ≧-0.44)

f(n)

= 1.0

f(n)

= 1.22 – 0.5 n(when n < -0.44)

n = :/AcF = ratio of the axial force in the column to the yielding force.
:otations :Fyield strength of column (refer to Table 6.2 for appropriate
value):
axial force in the column (compression is expressed with minus(-))
Ac
sectional area of the column
Pa2 should satisfy the condition given below.
(6.3.29)
Pa2 > Aρ F1
where: F1 is the yield strength of the bracings, (refer to Table 6.2 for the appropriate value)
In case the requirements of equation (6.3.29) cannot be satisfied, either the thickness of the column
wall should be increased at the joint, or the joint should be reinforced with a rib, refer to Figure
6.3.7.

gusset platerib

Fig.6.3.4 Reinforcement at the Column to Bracing Joint (Solid round bar)

Fig.6.3.5 Joint Detail of Steel Pipe Bracings

Chapter 676

Fig.6.3.6 Reinforcement at the Bracing Joint (Steel pipe)

ring stiffenerribFig.6.3.7 Reinforcement at the Column to Bracing Joint (Steel pipe)
References:
[6.1]
[6.2]
[6.3]
[6.4]
[6.5]
[6.6]

Akiyama, H. et al.: Report on Shaking Table Test of Spherical Tank, LP Gas Plant, Japan LP
Gas Plant Association, Vol. 10, (1973) (in Japanese )Akiyama, H.: Limit State Design of Spherical Tanks for LPG Storage, High Pressure Gas,The High Pressure Gas Safety Institute of Japan , Vol.13, :o.10 (1976), ( in Japanese )Architectural Institute of Japan: Lateral Retaining Strength and Deformation Characteristics
on Structures Seismic Design, Ch. 2 Steel Structures (1981), ( in Japanese )Kato, B. & Akiyama, H: Restoring Force Characteristics of Steel Frames Equipped with
Diagonal Bracings, Transactions of AIJ :o.260, AIJ (1977), ( in Japanese )Architectural Institute of Japan: Recommendation for the Design and Fabrication of Tubular
Structures in Steel, (2002), ( in Japanese )Plastic Design in Steel, A Guide and Commentary, 2nd ed. ASCE, 1971

Chapter 677

7.

Seismic Design of Above-ground, Vertical, Cylindrical Storage Tanks

7.1

Scope

7.1.1

This chapter covers the seismic design of above-ground vertical cylindrical tanks forliquid storage. General structural design are regulated by the relevant standards.

7.1.2

Notations
Ce
D
Hl
Ps
Pso
Pw
Pwo

ηs
Qy
Qdw
Qds

design shear force coefficientdiameter of the cylindrical tank (m)depth of the stored liquid (m)design dynamic pressure of the convective mass - the sloshing mass
(N/mm2)
reference dynamic pressure of the convective mass - the sloshing mass
(N/mm2)
design dynamic pressure of the impulsive mass - the effective mass of thetank contents that moves in unison with the tank shell (N/mm2)
reference dynamic pressure of the impulsive mass - the effective mass ofthe tank contents that moves in unison with the tank shell (N/mm2)
height of sloshing wave (m)design yield shear force of the storage tank (kN)design shear force of the impulsive mass vibration - the effective mass ofthe tank contents that moves in unison with the tank shell (kN)design shear force of the convective mass vibration - the effective mass ofthe first mode sloshing contents of the tank (kN)

Commentary:
There are two main types of above-ground cylindrical tanks for liquid storage; one is for oil
storage and the other is for liquefied gas. Many oil storage tanks are generally placed onearth foundations with a concrete ring at periphery. They are usually not equipped withanchoring devices. In case of oil storage tanks having floating roofs, the roof moves up anddown following the movements of the oil surface. (See Fig.7.1.1) Whilst, liquefied gas storage
tanks is usually a double container type tank which consists of an inner tank having a fixeddome roof and an outer tank also having a domed roof. They are usually placed on concreteslab foundations. The gas pressure in the inner tank, in which the liquefied gas is contained, isresisted by the inner dome roof and with anchoring devices provided at the bottom of theinner shell not to lift up the bottom plate (See Fig.7.1.2).

floating roof
wall
bottom platefoundation
Fig.7.1.1 Oil storage tank

Chapter 778

outerouterwallwall
inner
inner
wallwallbreathing tank

N2gas
foundation
Fig.7.1.2 Liquefied gas storage tank

These two different structures give different seismic response characteristics.
The rocking motion, accompanied by an uplifting of the rim of the annular plate or the bottom
plate, is induced in oil storage tanks by the overturning moment due to the horizontal inertia
force. In this case, particular attention should be paid to the design of the bottom corner(ⓐ)
of the tank as shown in Fig.7.1.3, because this corner is a structural weak point of verticalcylindrical storage tanks.

wall
a annular platebottom platecircumferential RC ringfoundation

Fig.7.1.3

Bottom corner of tankOn the other hand the rocking motion in the liquefied gas storage tank caused by theoverturning moment induces pulling forces in anchor straps or anchor bolts in place of
uplifting the annular plate. In this case, the stretch of the anchor straps or anchor bolts which
are provided at the bottom course of the tank, should be the subject of careful design.
Figure 7.1.4 shows several examples of tank damages observed after the Hanshin-Awaji
Great Earthquake in Japan in 1995 [7.1]. Figure 7.1.4(a) shows an anchor bolt that was
pulled from its concrete foundation. Figure 7.1.4(b) shows an anchor bolt whose head was
cut off. Figure 7.1.4(c) shows the bulge type buckling at the bottom course of tank wall plate,(Elephant’s Foot Bulge - E.F.B.). Figure 7.1.4(d) shows diamond pattern buckling at thesecond course of tank wall plate. Figure 7.1.4(e) shows an earthing wire that has been pulled
out when the tank bottom was lifted due to the rocking motion.
The type and extent of this damage shows the importance of seismic design for storage tanks.
Especially the rocking motion which lifts up the rim of the tank bottom plate, buckles thelower courses of tank wall, and the pulls out the anchor bolts or straps from their foundations.
The contained liquid, which exerts hydrodynamic pressure against the tank wall, is modelled
into two kinds of equivalent concentrated masses: the impulsive mass m1 that represents
impulsive dynamic pressure, and the convective mass ms that represents convective dynamic
pressure, (Housner) [7.2].

Chapter 779

(b) Cut off of anchor bolt(a) Pulled out of anchor bolt(c) Elephant foot bulge type buckling of tank wall(d) Diamond pattern buckling of tank wall(e) Pulled out of earthing wire
Fig.7.1.4 Tank damage mode at Hanshin-Awaji Earthquake

Chapter 7

80

The short-period ground motion of earthquakes induce impulsive dynamic pressure. Impulsive
dynamic pressure is composed of two kinds of dynamic pressures (Velestos)[7.3], as shown in
Fig.7.1.5. That is, the impulsive dynamic pressure in rigid tanks, pr , and the impulsive
dynamic pressure in flexible tanks, pf . The ordinates are for the height in the liquid with
z/Hl=1.0 for the surface, and D/Hl is the ratio of the tank diameter to the liquid depth. The
long-period ground motion of earthquakes induce convective dynamic pressure to tank walls
as the result of the sloshing response of the contained liquid in the tank. Moreover, thesloshing response causes large loads to both floating and fixed roofs.
The effect of long-period ground motion of earthquakes has been noticed since the 1964
:iigata Earthquake. As well as this earthquake, the 1983 :ihonkai-chubu earthquake and the2003 Tokachioki Earthquake caused large sloshing motion and damage to oil storage tanks.
Particularly, the Tokachioki Earthquake (M8.0) caused damage to oil storage tanks at
Tomakomai about 225km distant from the source region. Two tanks fired (One was a ring fire
and the other was a full surface fire.) and many floating roofs submerged. Heavy damage offloating roofs was restricted to single-deck type floating roofs. :o heavy damage was found in
the double-deck type floating roofs.

pfpr(a) Distribution in flexible tank(b) Distribution in rigid tankFig.7.1.5 Impulsive pressure distribution along tank wall

These two kinds of dynamic pressures, impulsive and convective, to the tank wall cause arocking motion which accompanies a lifting of the rim of the annular plate or bottom plate by
their overturning moments for unanchored tanks. Wozniak [7.5] presented the criteria forbuckling design of tanks introducing the increment of compressive force due to the rocking
motion with uplifting in API 650 App.E [7.6] . Kobayashi et al. [7.7] presented a mass-springmodel for calculating the seismic response analysis shown as Fig.7.1.6. The Disaster
Prevention Committee for Keihin Petrochemical Complex established, A Guide for Seismic
Design of Oil Storage Tanks which includes references to the rocking motion and bottom plate
lifting [7.8].

a(a) Uplift region(b) Rocking model
Fig.7.1.6 Rocking model of tank

Chapter 781(c) Spring mass system

This recommendation presents the seismic design criteria that takes into account both therocking motion with uplifting of bottom plate for unanchored tanks, and the elasto-plastic
stretch behaviour of anchor bolts or anchor straps for anchored tanks.

7.2

Seismic Loads

7.2.1

Seismic loads for vertical cylindrical tanks used for liquid storage should be
calculated using the procedure described in section 3.6.1. However, in case sloshing
response is taken into consideration, the velocity response spectrum shown in section
7.2.3. should be used.

7.2.2

Design dynamic pressure of the impulsive mass should be calculated using themodified seismic coefficient method in 3.6.1.2.
(1) Design dynamic pressure of the impulsive mass Pw is calculated using equations (7.1)
and (7.2). Design shear force of the impulsive mass vibration Qdw is obtained byintegrating the component in seismic force direction of Pw around the inside of tankshell.
(7.1)

Pw = Ce ·Pwo

(7.2)

Ce = Zs·I·Ds·Sa1/g
where,
Ce ≧ 0.3Zs·I
Notations:Cedesign shear force coefficientreference dynamic pressure of the impulsive mass (thePwoeffective mass of the tank contents that moves in unison withthe tank shell) (N/mm2)
seismic zone factor
Zs
I
importance factorDs
structural characteristic coefficient
Sa1
acceleration response spectrum of the first natural period
obtained from equations (3.8) or (3.9) (m/s2)gacceleration of gravity; 9.8 m/s2

(2) Structural characteristic coefficient, Ds , should be calculated using one of thefollowing two methods:a) For unanchored tanks, Ds is obtained from equations (7.3), (7.4) and (7.6),
In the case that the yield ratio of the annular plate, Yr , is smaller than 0.8, Ds is:

(7.3)

Ds =

1.42
1 + 3h + 1.2 h

⋅

1
T 
1 + 84 1 
 Te 

2

In the case that Yr is greater than 0.8, Ds is:
Chapter 782

(7.4)

Ds =1.42
1 + 3h + 1.2 h

⋅

1
T 
1 + 24 1 
 Te 

2b) For anchored tanks, Ds is obtained from equations.(7.5) and (7.6).
(7.5)

Ds =

1.42
1 + 3h + 1.2 h⋅

1
3.3 tσ y2  T1  
1+la paσ y  Te 2c) Buckling of tank wall plate should be examined with Ds in equation (7.6).
(7.6)

Ds =

1.42
1 + 3h + 1.2 h⋅

1 Tf 
1 + 3  Te 

2

Notations:hdamping ratio
T1 natural period of the storage tank when only the bottom plate
and anchor bolts (or anchor straps) deform (s)Te modified natural period of the storage tank considering
deformation of the wall plate as well as the bottom plate andanchorage (s)2σ y yield stress of the bottom plate (N/mm )2

σ y yield stress of anchor bolt or anchor strap (N/mm )la effective length of anchor bolt or anchor strap (mm)tthickness of the bottom plate (mm)pstatic pressure imposed on the bottom plate (N/mm2)
Yr the ratio of yield stress to ultimate tensile strength of annularplate
Tf fundamental natural period related to tank wall deformation
without uplifting of bottom plate.
a

7.2.3

Design seismic loads of the convective mass should be calculated using theprocedure described in section 3.6.1.
However, in case sloshing response by long-period ground motion of earthquakes is taken
into account, the velocity response spectrum shown in equation (7.7) in item (1) should be
used for periods longer than 1.28s. This spectrum depends on the damping ratio ofsloshing response shown in item (2).
(1) Velocity response spectrum for the first mode of sloshing is obtained from equation(7.7) in case the damping ratio is 0.5%.
 IS v = 2(m / s ) 1.28s ≦ T ≦11ｓ
(7.7)
 IS = 22 (m / s )11s ≦ T v T
Notations : T period (s)Chapter 783For cases of other damping ratio than 0.5%, the velocity response is varied according to
the damping ratio. However, when the sloshing response of higher mode, higher than the
1st, is calculated, 2m/s should be adopted at the minimum for the velocity response even if
less velocity is obtained in accordance with the larger damping ratio. This velocity
response includes an importance factor, I, and therefore means IS v .
The seismic zone factor, Z, should be 1.0 except there is special reason and basis. The
velocity response spectrum for cases of Z=1.0 are shown in figure 7.2.7 in commentary.
(2) The damping ratio for the 1st mode sloshing should be less than 0.1% for fixed rooftanks, 0.5% for single-deck type floating roof tanks and 1% for double-deck type floating
roof tanks. For the higher modes, the damping ratio should be less than the 1st mode
damping ratio multiplied by the frequency ratio, the ratio of the frequency of the mode to
the 1st mode frequency. The 1st mode damping ratios are shown in table 7.2.1 and thedamping ratios for higher modes are shown in figure 7.2.8 in commentary.
(3) Predicted strong motions considering precisely the source characteristics and patheffects can be used as design seismic ground motions according to the judgment of thedesigner instead of the velocity response spectrum shown in equation (7.7).
(4) Design dynamic pressure of the convective mass Ps is calculated by using equation(7.8). Design shear force of the convective mass vibration Qds for sloshing is obtained byintegrating Ps around the inside of tank shell.
(7.8)

Ps =

η s ⋅ Pso
D(5) Height of sloshing wave η s is calculated by using equation.(7.9).
(7.9)

η s = 0.802 ⋅ Z s ⋅ I ⋅ S v1

 3.682 H l 
Dtanh 
gD
Notations:

ηs
D
Pso
Hl
S v1

height of sloshing wave (m)diameter of the cylindrical tank (m)reference dynamic pressure of the convective mass (N/mm2)
depth of the stored liquid (m)design velocity response spectrum of the first natural period
in the sloshing mode as calculated by equations (7.7) (m/s)

7.2.4

Stress and Response of Floating Roofs
Concerning floating roof type tanks, overflow of contained liquid and damage offloating roofs by sloshing should be checked at design stage. When stresses in floatingroofs by sloshing response are calculated, higher modes response as well as the 1st mode
response should be taken into consideration.

7.2.5

Sloshing Pressure on Fixed Roofs
Impulsive pressure and hydrodynamic pressure by sloshing on fixed roofs should be taken
into consideration for the seismic design of fixed roof tanks.
Commentary:

7.2.2 1.

Design dynamic pressure

Chapter 7

84

The response magnification factor of tanks is approximately 2, when the equivalent damping ratioof a tank is estimated to be 10% [7.9]. The reference dynamic pressure of the impulsive mass Pwo
is calculated as follows:
The impulsive mass mf is obtained by integrating equation (7.2.1) over the whole tank wall.
(7.2.1)

Pwo =

2

pf +

(

)

1
2
p r − p f cos ϕ4Base shear SB is obtained from equations (7.2.2) and (7.2.3).
(7.2.2)m f = ∫∫ (Pwo / g ) cos ϕrdϕdz = f f ⋅ ml
where, ml is total mass of contained liquid, and ff is given in Fig.7.2.1.

(7.2.3)

S B = Ce ⋅ g ⋅ m f

Reference dynamic pressure of the convective mass Pso and natural period of sloshing response Ts
are calculated as shown in equations (7.2.4) and (7.2.5).
(7.2.4)ZH Pso = γ ⋅ g ⋅ D cosh 3.682  cosϕ cosh 3.682 l 
D
D



(7.2.5)Ts = 2πD
3.682 g tanh (3.682 H l D )where:
DHlγ2.

diameter of the cylindrical tankdepth of the stored liquiddensity of the contained liquid.

Coefficient determined by the ductility of tank Dη

2.1 Unanchored tank
As the failure modes of the unanchored tank are the yielding of the bottom corner due to uplifting,and the elephant foot bulge type buckling on the tank wall, Dη can be obtained from equation (7.6)
applying Dη f = 0.5 in equation (3.7.18) .
The following paragraphs show how to derive Dη for the uplifting tank.
The elastic strain energy We stored in the uplifting tank is as shown in equation (7.2.6),(7.2.6)We =k eδ e22where, ke is the spring constant of tank wall - including the effect of uplifting tank, and δ e is theelastic deformation limit of uplifting tank.
The relation of the uplifting resistance force, q, and the uplifting displacement, δ , per unit
circumferential length of tank wall in the response similar to that shown in Fig.7.1.6(b) is given as
the solid line in Fig.7.2.2.

Chapter 785

Fig.7.2.2 q vs δ relation

Fig.7.2.1 Effective mass ratio

The rigid plastic relation as shown by the dotted line is introduced for simplifying, the plastic
work Wp1 for one cycle of horizontal force is as shown in equation (7.2.7)(7.2.7)

W p1 = 2πrq y δ B
where:rqy

δBradius of the cylindrical tankyield force of the bottom plate
the maximum uplifting deformation,

Cumulative plastic strain energy Wp is given as twice Wp1 (Akiyama [7.10]) as shown in equation(7.2.8).
(7.2.8)W p = 2W p1 = 4πrq yδ B

When the bottom plate deforms with the uplifting resistance force, q, and liquid pressure, p, asshown in Fig.7.2.3, the equilibrium of the bottom plate, as the beam, is given in (7.2.9)(7.2.9)

EIy ′′′′ + p = 0

Deformation y is resolved in equation (7.2.10) under boundary conditions y = y ′ = 0 at x=0,y ′ = y ′′ = 0 at x=l.
(7.2.10)y=(1− px 4 24 + plx 3 9 − pl 2 x 2 12
EI)

where, EI is bending stiffness.
The uplifting resistance force q, the bending moment at the fixed end Mc and the relationshipbetween q and δ are obtained as follows:q= 2pl/3

Mc= pl2/6

δ= 9q4/128EI p3.

When Mc comes to the perfect plastic moment, Mc is equal to σ y t2/4. In this condition, the yielding
force of the bottom plate qy , the uplifting deformation for the yielding force δ y , and the radiusdirection uplifting length of bottom plate ly are obtained from equations (7.2.11).
q y = 2t 1.5 pσ y 3(7.2.11)δ y = 3tσ y2 8 Ep
l y = t 3σ y 2 p

Chapter 7

86

Mc
Mc

βφy

Fig.7.2.4 Mc vs φ relation in bottom plateFig.7.2.3 Bottom plate deformation

When the relationship between the bending moment Mc and the curvature ϕ of elasto-plastic beam
is as indicated in Fig.7.2.4, the relationship between q and δ , and Mc and δ are obtained bynumerical calculation, refer to Fig.7.2.5. In these figures, β denotes the ratio of the curvature at
the beginning of strain hardening to that at the elastic limit, and est denotes the ratio of the stiffness
at the beginning of strain hardening to that at the elastic limit.
The values β = 1.0 and est = 0.005 are applied for high tensile steel, and β = 1.0 to 11.0, and est= 0.03 are applied for mild steel. Equation.7.2.12 is a reasonable approximation of therelationship between Mc and δ , which is shown with dot-dash line in Fig.7.2.5(b).
Mδ= 1 + 32 c − 1


δy
My


Mc/My

(7.2.12)(b) Mc vs δ relation(a) q vs δ relation

Fig.7.2.5 Wall reaction and bending moment at bottom corner
As the bottom plate failure is considered to be caused when Mc reaches the fully plastic moment
My , which is obtained from the tensile strength of bottom plate material, the maximum value ofMc / My is given as(7.2.13)(Mc / My)max = σ B / σ ywhere, σ B is the tensile strength of bottom plate material and is defined as (σ B / σ y ) being equal
to 1.45 for materials of which the ratio of yield stress to ultimate tensile strength of annular plate
Yr is smaller than 80%, or 1.10 for materials of which Yr is greater than 80%.
By substituting the relationships in equation (7.2.12), the limit uplift displacement δ B is obtained
as shown in equation (7.2.14).

Chapter 787

(7.2.14)

δ B = 14σ yδ B = 4σ yfor Yr ≤ 80% for Yr ≥ 80% 

The spring constant of unit circumferential length for uplifting resistance k1 is defined as thetangential stiffness at the yield point of q- δ relation, and is shown in equation (7.2.15).
(7.2.15)

k1 =

qy

δy

=

16 pE 1.5 p
9σ yσyThe overturning moment M can be derived as the function of the rocking angle of tank θ as shownin equation (7.2.16),(7.2.16)2π∫ k rθ (1 + cosϕ ) r (1 + cosϕ ) rdϕ = K θ
K = 3∫ k r sin ϕ rdϕ = 3k π r

M=

0

1

M

2π

M

0

2

2

3

1

1

where, KM is the spring constant that represents the rocking motion of tank. The spring constant K1
as the one degree of freedom system that represents the horizontal directional motion of tank isgiven as,(7.2.17)

K1 = M (0.44 H l ) θ = 48.7r 3 k1 H l22

where 0.44Hl is the gravity centre height of impulsive mass.
Then the natural period of the storage tank when only the bottom plates deform, T1 , is derived form
equation (7.2.18).
(7.2.18)

T1 = 2π mt K 1

where, mt is the effective mass which is composed of the roof mass mr , mf and the mass of wall
plate mw.
The modified natural period of the storage tank considering deformation of both the wall plate andthe bottom plate, Te , can be written as equation (7.2.19).
(7.2.19)

Te = T f2 + T12

In this equation, Tf is the fundamental natural period related to tank wall deformation without
uplifting of bottom plate (JIS B8501[7.4] , MITI [7.9]), and is given as(7.2.20)

Tf =

2

λ

m0 π E t1 / 3

where, t1/3 is the wall thickness at Hl/3 height from the bottom, m0 is the total mass of tank which iscomposed of ml , mw and the mass of roof structure mr , and λ is given as following equation,

λ = 0.067(H l D )2 − 0.30(H l D ) + 0.46 .
As the equivalent spring constant considering deformation of both the wall plate and the bottom
plate, ke , can be derived from equation (7.2.21),(7.2.21)

Ke = K1(T1/ Te)2

the elastic deformation of the centre of gravity of the tank, δ e , is determined by using the design
yield shear force of the storage tank, Qy , which is further described in section 7.3 below, andequation (7.2.22).
(7.2.22)

Chapter 7

δ e = Q y / ke88

By substituting equations (7.2.21) and (7.2.22) into equation (7.2.6), We can be calculated, andfrom equation (7.2.14), Wp can be calculated.
Dη , in equation (3.6.15), can then be written as equation (7.2.23).
(7.2.23)1Dη =1 + W p We

Then, the coefficient determined by the ductility of the bottom plate uplifting, Dη , is obtained inaccordance with the equations (7.2.24).

1

Dη =

1 + 84(T1 Te )12

(7.2.24)Dη =1 + 24(T1 Te )2

for Yr ≤ 80%
for Yr ≥ 80%Finally, the structural characteristic coefficient Ds is determined as shown in equations (7.3) and(7.4) by multiplying Dh and Dη .
2.2

Anchored tank
The yield force of anchor straps for unit circumferential length aqy is given as equation (7.2.25)(7.2.25)

Aa ⋅a σ y − π r 2 Pia qy =

2π r

where:

σyyield stress of anchor straps

Aa
Pi

sectional area of anchor straps
internal pressure of tank.

a

From Hook’s low, the anchor strap deformation for yield point, a δ y , and the spring constant ofanchor straps for unit circumferential length, ak1 , are calculated in accordance with equations(7.2.26) and (7.2.27).

δ y = la ⋅a σ y E(7.2.26)a(7.2.27)

a k1 =

a qy
aδ y

A ⋅E
= e a
2π rlawhere:
la is the effective length of anchor bolt or anchor strap,
and e Aa = Aa (1 − π r 2 Pi / Aa ⋅a σ y ) .
When the deformation limit of the anchor is set to be δ B , the mean deformation of anchor strap δ B
is calculated as δ B = (2 + π )δ B / 2π . By referring to Fig.7.2.6, δ B is determined as one half of thebottom plate displacement at which point the second plastic hinge occurs in the bottom plate. The
white circles in Fig.7.2.5 show the second plastic hinge, and δ is approximately 7δ y . Then δ B
can be rewritten as equation (7.2.28).

Chapter 789

Fig.7.2.6 Average elongation of anchor strap(7.2.28)δB =(2 + π ) δ B = 1.08 tσ y2 = 1.08 tσ y2 ⋅a δ y
2π

la p⋅a σ yEp

By referring to equation(7.2.7), the plastic work Wp becomes
W p = 2π r ⋅a q y ⋅ δ B .
Then, the coefficient determined by the ductility of the anchor strap elongation, Dη , and the spring
constant K1 are determined in accordance with equations (7.2.29) and (7.2.30) respectively.

1

Dη =
1+
(7.2.29)2 
la p ⋅a σ y  Te 

Te = T f2 + T12
T1 = 2π(7.2.30)

3.3 tσ y2  T1 

mt
K1















K1 = 48.7 ra3 k1 H l2

Finally, the structural characteristic coefficient, Ds , is determined as equation (7.5) by multiplying
Dh and Dη .
2.3

Elephant foot type buckling of wall (E.F.B.)
When the uplifting of bottom plate is caused due to the overturning moment, the coefficient
determined by the ductility of the wall plate of the cylindrical tank, Dη is given as equation (7.2.31)
from equation (3.7.18).
(7.2.31)

Dη =

1

(

)

1 + 3 T f Te 2

Then, the structural characteristic coefficient, Ds , is determined as equation (7.6) multiplying Dhand Dη .
7.2.3
1. Damping ratio of sloshing
Concerning the damping ratios of sloshing for oil storage tanks, experimental data with actual andmodel tanks are obtained and published 7.12), 7.13). Damping ratios for the 1st mode are determined basedon mainly earthquake observations and forced vibration tests. Although there are few questions to
support the 0.1% or less damping ratio for fixed roof tanks (free surface), neither experimental nor
observational data is sufficient to support the 0.5% damping ratio for single-deck type floating roofs andthe 1.0% damping ratio for double-deck type.

Chapter 7

90

Table7.2.1 Damping Ratio for the 1st Mode Sloshing
Surface (Roof) Condition
Free Surface (Fixed Roof)Single-deck (Floating Roof)
Double- Deck (Floating Roof)

Damping Ratio
0.1%
0.5%
1.0%

Concerning the damping ratios for higher modes, some data for actual and model tanks show 5%
or more damping ratio on the one hand, but other data show that the same damping ratio as the 1st mode
gives the best simulation of the earthquake response of actual tanks.
Therefore, the damping ratios given for the higher modes have not much safety margin, and enoughmargin for the 1st mode seismic design is recommended. Figure 7.2.8 shows the upper limits of thehigher mode damping ratios.
2. Velocity response spectrum for sloshing response
Simple spectrum is regarded as preferable. Therefore, 1.0m/s uniform velocity spectrum forengineering bedrock is assumed for 5% damping response with coefficient 1.2 for importance factor,
modification coefficient 1.29 for 0.5% damping response and 1.29 for large scale wave propagation andbi-directional wave input effects. For damping ratios other them 0.5%, 1.10 /(1 + 3h + 1.2 h ) should be
multiplied as below and Figure 7.2.7.
1.10(7.2.32)
ISv = 2 ⋅[m/s]
1 + 3h + 1.2 h
Bi-directional effect is considered because tank structures are usually axisymmetric and the
maximum response in some direction usually different from given X-Y direction occurs without fail. If
the earthquake waves for X and Y directions are identical, the coefficient for this effect is 2 , but
usually it ranges from 1.1 to 1.2.
As explained above, importance factor, 1.2, is already considered. Seismic zone factor Zs is fixedat 1.0 with the exception of particular basis because there are insufficient data for long-period groundmotion.
In case velocity response spectrum for periods shorter than 1.28s is necessary, to obtain higher
mode sloshing response, the velocity response spectrum is calculated from the 9.8m/s2 uniform
acceleration response spectrum if the damping ratio is 0.5% or more. In case the damping ratio is lessthan 0.5%, the same modification as equation (7.2.32) is necessary.
damping ratio(%)

ZsISv (m/s)

2.2
2.0
1.8

free surface
single-deck
2.11
2.00
1.91

double-deck

reciprocally proportional
to period above 11s

single-deck

1.6
1.4

0

1.28

double-deck

1.0

11

0.5
0.1
0

T(s)

13

15

Fig.7.2.7

free surface
1.0

frequency ratio to
the 1st frequency

Velocity Response Spectrum for the 1st
Fig.7.2.8 Upper Limits for the Higher Mode
Mode SloshingDamping Ratios
Figure 7.2.9 shows the spectrum with typical observed and predicted long period earthquake
ground motions 7.13). Compared with the observed earthquakes and considering the bi-directional effect,
the proposed spectrum shows appropriate relations in magnitude to the observed. On the other hand,some of the predicted earthquake spectra exceed the 2.0m/s level in magnitude.

Chapter 791

3. :atural period of higher mode sloshing
The i-th natural period of free liquid surface sloshing in cylindrical tank is calculated by theequation below.
(7.2.33)

Ti =

2π D
2 gε i tanh( 2ε i

Hl)DWhere :
Ti
i-th natural period of sloshing (s)gacceleration of gravity (m/s2)
D
diameter of the cylindrical tank (m)εii-th positive root of the equation, J '1 (ε i ) = 0 .
J1 is Bessel function of the first kind. ε1 = 1.841 , ε 2 = 5.33 , ε 3 = 8.54
In case the tank has a single-deck type floating roof, the natural sloshing periods are calculated byequation (7.2.33), in other words, the effect of the roof is negligible. In double-deck type floating roof
cases, only the effect to the 1st mode is negligible.
7.2.4. Response and Stress Analysis of sloshing in a tank with a floating roof
Stress in peripheral pontoon and roof is calculated by recently proposed analysis methods andsimplified approximate calculation methods. In these analysis methods, nonlinear FEM for liquid andelastic body, velocity potential theory, calculus of variation, boundary integral, and Bessel function
expansion of the roof modes are used.

Chapter 792

5K-NET Tomakomai・NS
K-NET Tomakomai・EW
near Kashiwazaki NPP Unit-1・NS
near Kashiwazaki NPP Unit-1・EW

pSv(m/s)

4
3
2
1
0

0

5

10

T(s)

15

comparison with observed earthquakes (26 Sep. 2003 Tokachi-oki (Tomakomai)
and 16 Jul. 2007 :iigataken Chuetsu-oki (Kashiwazaki))
5

TS・TOK・NS
TS・YKL・NS
KH・SNJ・NS
KH・STY・NS
KH・YKH・NS

pSv(m/s)

4
3

2
1
0

0

5

10

T(s)

15

comparison with predicted earthquakes for Tokyo and Yokohama region
C・SAN・EW
A・NST・EW
A・SJB・EW
KK・OSA・NS
KK・WOS・EW
HS6・FKS・NS
HS6・OSK005・EW
HS6・OSKH02・NS

5

pSv(m/s)

4
3

2
1
0

0

5

10

T(s)

15

comparison with predicted earthquakes for :agoya and Osaka
Fig.7.2.9 Comparison of the Velocity Response Spectrum with Typical Observed andPredicted Long Period Earthquakes

Chapter 7937.2.5 Sloshing force acting on the fixed roof of the tank
When the sloshing is excited in the tank with the fixed roof subjected to the earthquake motion, thecontained liquid hits the fixed roof as shown in Fig.7.2.10 in the case that the static liquid surface is close to
the fixed roof. Figure 7.2.11 sketches the typical time evolution of the dynamic pressure acting on the point A ofthe fixed roof 7.14). The impulsive pressure Pi acts on the fixed roof at the moment that the sloshing liquid hits it,
and the duration of acting of the impulsive pressure is very short 7.14), 7.15). The hydrodynamic pressure Ph actson the fixed roof after impulsive pressure diminishes as the sloshed contained liquid runs along the fixed roof
7.14)
. Since the duration of acting of the hydrodynamic pressure is long, the hydrodynamic pressure can be
regarded as the quasi static pressure.

Fixed roof

ζR

A

ϕ

ζr

h

Liquid surfacerR
Sloshing surface

Fig.7.2.10

Fixed roof and sloshing configuration

Impulsive pressurehydrodynamic pressure

Pi
Ph
time
Fig.7.2.11 Impulsive and hydrodynamic pressure
The behavior of the sloshed liquid depends on the shape of the fixed roof. When the roof angle ϕ is large,
the restraint effect from the roof on sloshing is small as shown in Fig. 7.2.12. Then the sloshed liquid elevatesalong the roof with small resistance from the fixed roof, and yields the hydrodynamic pressure Ph . When theroof angle ϕ is small, the restraint effect on sloshing is large and the sloshing liquid hardly increases it’selevation at all as shown in Fig. 7.2.13. Then small hydrodynamic pressure Ph yields on the fixed roof of whichangle ϕ is small 7.14), 7.16), 7.17).
The impulsive pressure Pi acting on the point A of a cone roof tank is given by Yamamoto [7.15] as(7.2.34)

Pi =

π
2

( )

ρ cot ϕ ς&rA2[:/m2],

where ς&r is the sloshing velocity at the point A. ς r is the sloshing height at the radial distance r fromthe center of the tank 7.15). The first sloshing mode is enough to evaluate the sloshing velocity and the magnitude
Aof the impulsive pressure. According to Yamamoto [7.15], Eq. (7.2.34) is applicable to the roof angle ϕ ≥ 5 .
oChapter 794

Fig.7.2.12

Sloshing at cone roof with 25deg slope 7.17)

Fig.7.2.13

Sloshing at plate roof 7.17)

Kobayashi [7.14] presented that the impulsive pressure Pi acting on the dome and cone roof tank can be
estimated by Eq. (7.2.34), using the maximum sloshing velocity response without the fixed roof subjected to thesinusoidal input of which the period coincides with the sloshing natural period and three waves repeat.
Transforming the impulsive pressure shown in the literature [7.14] represented by the gravimetric unit to SI
unit, the following relation is obtained;(7.2.35)( )Pi = 34.97 ρ ς&rA

1 .6[:/m2].

Since the sloshing response is sinusoidal, the sloshing velocity ς&r at the point A is approximately given
as Eq. (7.2.36) using the distance h from the static liquid surface to the fixed roof;A(7.2.36)



ς&rA = ς r ω s cos sin −1h ς r 

[m/s]

where ωs is the circular natural frequency of the sloshing.
From these discussions, the impulsive pressure and the hydrodynamic pressure at a point oriented anangle θ to the exciting direction is given as follows;
1. in the case of ϕ ≥ 5o
Impulsive pressure:(7.2.37)

Chapter 7

Pi =

π
2

(

ρ cot ϕ ς&rA cosθ)2[:/m2]

95

hydrodynamic pressure:
Ph = ρg (ς r − h )cosθ(7.2.38)[:/m2]

2. in the case of ϕ < 5o
Impulsive pressure:(

Pi = 34.97 ρ ς&rA cos θ(7.2.39))

1 .6[:/m2]hydrodynamic pressure: excluded
Since the hydrodynamic pressure varies slowly according to the natural period of the sloshing, theearthquake resistance of the fixed roof can be examined by integrating the hydrodynamic pressure along thewetted surface of the roof statically.
Since the magnitude of the impulsive pressure is high and the duration of action is very short, theearthquake resistance examination of the fixed roof by the impulsive pressure should take the dynamic responseof the fixed roof into account 7.17). For example, if the fixed roof is assumed to be a single degree of freedom
system and the impulsive pressure is approximated as the half-triangle pulse shown in Fig.7.2.14, the transient
response of the system is obtained from the response to the half-triangle pulse load in the very short period ∆t.
Thus the impulse response of the fixed roof can be calculated using the load fs which is obtained by integrating
the impulsive pressure along the wetted surface of the roof.
When the equation of motion of a single degree of freedom system is given as Eq. (7.2.40);
mr &x& + cr x& + k r x = f (t ) ,(7.2.40)

where mr, cr and kr are the equivalent mass, the equivalent damping coefficient and the equivalent stiffness ofthe fixed roof, the response x to the half-triangle pulse under the initial conditions x = 0, x& = 0 is obtained asfollowing using the Duhamel integral;x (t ) =(7.2.41)
4I
 ∆t 
2 sin ω n  t −  − sin ω n (t − ∆t ) − sin ω n t  ,2 2 
k r ω n ∆t 



f

fs
1
f s ∆t = I
2

t
∆t

0

Fig.7.2.14

Pressure shape of half-triangle pulse

where ωn is the natural circular frequency of the fixed roof, I=(1/2)fs∆t is the magnitude of the impulse during
the small time interval ∆t and t>0.
In the case that ωn is sufficiently smaller than 2π /∆t, the impulse is assumed to be a delta function ofmagnitude I=(1/2)fs∆t and Eq. (7.2.41) is simplified into Eq. (7.2.42).
(7.2.42)x (t ) =f
1∆tω n ⋅ s sin ω n t .
2
kr

From Eq. (7.2.42), multiplying the magnitude of response (1 / 2)∆tω n by pressure given as Eqs. (7.2.37) and(7.2.39) gives the equivalent pressure. The resultant force yielded by the equivalent pressure is obtained from

Chapter 796integrating along the wetted surface of the fixed roof. Stresses in the fixed roof due to the resultant force canbe obtained using the resultant force as the static force.
Since the rising time of impulsive pressure is commonly shorter than 0.01s 7.17), it is recommended to set ∆t
as 0.02s on the safe side. ωn should be calculated according to the actual structural characteristics of the fixedroof. For example, if ∆t is 0.02s and ωn is 12rad/s, the magnitude of response is reduced to 12% of thatobtained on the assumption of static fs itself acting. This result is obtained, from Eq. (7.2.42), on the assumptionthat the envelope of the maximum value of the impulsive pressure obtained by Eqs. (7.2.37) or (7.2.39) acts
simultaneously over the wetted surface of the fixed roof.
The maximum impulsive pressures, actually, never act simultaneously all over the fixed roof, because thecontact point, where the impulsive pressure yields, moves as the cap of the sloshing wave moves. Though thisassumption seems too safe side to design, it is recommended to keep a margin of an unpredictable response
magnification due to such as the moving force effect.
A guide to estimate the resultant force by applying Eq. (7.2.42) is that the natural frequency of the fixedroof is less than 15Hz when ∆t is 0.02s. When the natural frequency of the fixed roof is more than 15Hz, themaximum response during 0 ≤ t ≤ ∆t can be estimated using Eq. (7.2.41) and the magnitude of response to thestatic force can be obtained. Alternatively, the magnitude of response can be estimated by interpolating linearly
between 15Hz to 45Hz of the natural frequency of the fixed roof, assuming the magnitude of response(1 / 2)∆tω n is 1.0 in the case of 15Hz and (1 / 2)∆tω n is 1.5 in the case of 45Hz.

7.3

Evaluation of Seismic Capacity

7.3.1

Checking of Seismic Capacity of Whole Tank
Design shear force of the impulsive mass vibration should satisfy equation (7.10).
(7.10)fQy ≧Qdw

Design shear force of the convective mass vibration should satisfy equation (7.11).
(7.11)sQy ≧Qds

Notations:yield shear force of the storage tank at the impulsive mass vibration (kN)fQyyield shear force of the storage tank at the convective mass vibration (kN)sQy
Qdw
design shear force of the impulsive mass vibration (kN)
Qds
design shear force of the convective mass vibration (kN)
7.3.2

Checking of Seismic Capacity for Tank
Checking of seismic capacity for tank parts should be conducted by comparing
working stresses and strains at corresponding parts in the storage tank with theallowable values. Working stresses and strains should be obtained by combining
seismic load and ordinary operation load.

-----------------------------------------------------------------------------------------------------Commentary:
1

Design shear force
The design shear force of the impulsive mass vibration Qdw is given by the equation (7.3.1)(7.3.1)

Qdw= Ds·SB=Zs·I·Ds·Sa1·mf

The design hoop stress at the bottom course of tank wall Qσhd is given by the equation (7.3.2)(7.3.2)

Chapter 7

Q σ hd= Qdw 2.5 H l t 0 + ml g πrt 097where, t0 is the wall thickness at the bottom.
The design shear force of the convective mass vibration Qds is given by the equation (7.3.3)(7.3.3)

Qds=Zs·I·Sa1·ms

The yield shear force of the storage tank at the convective mass vibration, sQy , is approximately
given as shown in equation (7.3.4)(7.3.4)2sQy =0.44 fQy .

Design yield shear force
The design yield shear force of the storage tank, Qy , is determined in accordance with equation(7.3.5).
(7.3.5)Q y = 2πr 2 q y 0.44 H lThis formula is modified as shown in equation (7.3.6), when elephant’s foot bulge type buckling isevaluated,(7.3.6)e

Q y = πr 2 ⋅ e σ cr t 0 0.44 H lwhere, e σ cr is the critical stress for the elephant foot bulge type buckling.
It can be further modified as shown in equation (7.3.7), when the yield of the anchor strap isevaluated,(7.3.7)3

a Q y = 2πr2⋅a q y 0.44 H l

Evaluation flow
Figures 7.3.1 and 7.3.2 summarise the evaluation procedure for anchored and unanchored tanks
respectively.

References:
[7.1]
[7.2]
[7.3]
[7.4]
[7.5]
[7.6]
[7.7]

[7.8]
[7.9]
[7.10][7.11]

Kobayashi, :. and Hirose, H.: Structural Damage to Storage Tanks and Their Supports, Report on
the Hanshin-Awaji Earthquake Disaster, Chp.2, Maruzen (1997) pp353-378, ( in Japanese )Housner, G.W.: Dynamic Pressures on Accelerated Containers, Bull. Seism. Soc. Amer., Vol47,(1957) pp15-35
Velestos, A. S.: Seismic Effects in Flexible Liquid Storage Tanks, Proc. 5th World Conf. on
Earthquake Engineering, Session 2B (1976)
Japanese Industrial Standards Committee: Welded Steel Tanks for Oil Storage, JIS B8501, (1985),( in Japanese )Wozniak, R. S. and Mitchel, W. W.: Basis of Seismic Design Provisions for Welded Steel Storage
Tanks, API Refinery Department 43rd Midyear Meeting, (1978)
API: Welded Steel Tanks for Oil Storage, App.E, API Standard 650, (1979)
Ishida, K and Kobayashi, :.: An Effective Method of Analyzing Rocking Motion for Unanchored
Cylindrical Tanks Including Uplift, Trans., ASME, J. Pressure Vessel Technology, Vol.110 :o.1
(1988) pp76-87
The Disaster Prevention Committee for Keihin Petrochemical Complex, A Guide for Seismic
Design of Oil Storage Tank, (1984), ( in Japanese )Anti-Earthquake Design Code for High Pressure Gas Manufacturing Facilities, MITI (1981), ( inJapanese )Akiyama, H.: :o.400, AIJ (1989), ( in Japanese )Akiyama, H. et al.: Report on Shaking Table Test of Steel Cylindrical Storage Tank, 1st to 3rdreport, J. of the High Pressure Gas Safety Institute of Japan, Vol.21, :o.7 to 9 (1984) (inJapanese )Chapter 798[7.12][7.13][7.14][7.15][7.16][7.17]

H. :ishi, M. Yamada, S. Zama, K. Hatayama, K. Sekine : Experimental Study on the Sloshing
Behaviour of the Floating Roof using a Real Tank, JHPI Vol. 46 :o.1, 2008
AIJ (Architectural Institute of Japan) : Structural Response and Performance for Long Period
Seismic Ground Motions, 2007.12, (in Japanese):obuyuki Kobayashi: Impulsive Pressure Acting on The Tank Roofs Caused by Sloshing Liquid,
Proceedings of 7th World Conference on Earthquake Engineering, Vol.5, 1980, pp. 315-322
Yoshiyuki Yamamoto: in Japanese, Journal of High Pressure Institute of Japan, Vol.3, :o.1, 1965,
pp. 370-376
Chikahiro Minowa: Sloshing Impact of a Rectangular Water Tank : Water Tank Damage Causedby the Kobe Earthquake, Transaction of JSME, Ser. C, Vol.63, :o. 612, 1997, pp. 2643-2649, inJapanese
Yoshio Ochi and :obuyuki Kobayashi: Sloshing Experiment of a Cylindrical Storage Tank,
Ishikawajima-Harima Engineering Review, Vol.17, :o. 6, 1977, pp. 607-615

Chapter 799Design Spec.

1

Checking for impulsive mass
Checking of buckling capacity

Calculate K1 by Eq. (7.2.17)
Calculate e σ cr and Q σ hd by Eqs. (3.54) and (7.3.2)

Calculate Te by Eq. (7.2.19)

Calculate e Qy by Eq. (7.3.6)

Calculate Ds of buckling by Eq. (7.6)
Calculate Qdw by Eqs. (7.3.1)

1

e Q y > Qdw

No

Yes
Checking of bottom plate capacity

Calculate Ds of bottom plate by Eq. (7.3) or (7.4)

Calculate qy by Eq. (7.2.11)

Calculate Qdw by Eq. (7.3.1)

1

Calculate Qy by Eq. (7.3.5)

No

Qy > Qdw

Yes
Checking for convective mass

Calculate Qds by Eq. (7.3.3)

1

Calculate s Qy by Eq. (7.3.4)

No

s Q y > Qds

Yes
FinishFig.7.3.1 Evaluation flow for unanchored tanks

Chapter 7

100

Pressure to fixed roofby Eq.7.2.37 ~ 39Design Spec.

1

Checking for impulsive mass
Checking of buckling capacity
Calculate e σ cr and Q σ hd by Eqs. (3.54) and (7.3.2)

Calculate Ds of buckling by Eq. (7.6)
Calculate Qdw by Eqs. (7.3.1)

1

Calculate e Qy by Eq. (7.3.6)

e Q y > Qdw

No

Yes
Checking of anchor strap capacity

Calculate K1 by Eq. (7.2.30)
Calculate Te by Eq. (7.2.29)

Calculate a q y by Eq. (7.2.25)

Calculate Ds of anchor strap by Eq. (7.5)

Calculate a Q y by Eq. (7.3.7)

Calculate Qdw by Eq. (7.3.1)

1

No

a Q y > Qdw

Yes
Checking for convective mass

Calculate Qds by Eq. (7.3.3)

1

Calculate s Qy by Eq. (7.3.4)

No

s Q y > Qds

Yes
FinishFig.7.3.2 Evaluation flow for anchored tanks

Chapter 7

101

Pressure to fixed roofby Eq.7.2.37 ~ 39

8Seismic Design of Under-ground Storage Tanks

8.1

Scope
This chapter covers the seismic design of water tanks, cylindrical tanks, etc., which
are either completely buried, or have a considerable part of the outside surface
below ground level. For the purposes of this recommendation, the general structuraldesign is assumed to be regulated by the relevant codes and standards.

------------------------------------------------------------------------------------------------------------------------------Commentary:
1.

GeneralThe seismic advantage of under-ground tanks is that they are not subjected to the same
vibration amplification as above-ground structures. However, as the surrounding soil has asignificant influence on the seismic behaviour, this, together with the inertia force of thestored material, should be taken into account when designing for seismic loads.

2.

Definition of an under-ground tank
For the purposes of this recommendation, an under-ground storage tank may be defined as astorage tank which is for the most part covered by earth or similar materials and whose
seismic response depends on the surrounding soil. Figure 8.1.1 illustrate the general
concept of under-ground tanks.

Figure 8.1.1 Conceptual examples of under-ground tanks

Although under-ground tanks may be constructed from steel, concrete, wood, FRP, or othercombinations of materials, the outer part of them which is in contact with soil is generally of
concrete, e.g., under-ground steel tanks usually have outer under-ground containers ofconcrete. In view of this generally accepted method of construction, only concrete and prestressed concrete structures are considered in this chapter.

8.2

Seismic Design

8.2.1

Under-ground tanks should be constructed in satisfactory and stable ground. Softground, or inclined ground which is subject to sliding and slope failure should be
avoided, as should ground with abrupt changes in contours or materials, or atconnecting sections with different structures where significant deformation andstress can be expected. However, when construction on unstable ground isunavoidable, it is preferable that suitable soil improvement measures, piled
foundations, or other methods are adopted for higher seismic safety.

Chapter 8102

8.2.2The seismic intensity method, response displacement method or a dynamic analysis
may be used in the seismic design of under-ground tanks. A dynamic analysis ispreferable when the tank has an interface with other structures or is constructed insoil with abrupt changes of soil properties.
The design should fully consider the seismic loads resulting from the deformation of
the surrounding soil, from the inertia force of tanks and contained materials, andfrom the seismic earth pressure in the ultimate state for out-of plane stress of underground walls.

8.2.3

Soil properties for the seismic design model should include the effect of soil strainlevel during earthquakes.

8.2.4

Buoyancy should be taken into account if there is the possibility of surrounding soilliquefaction.

------------------------------------------------------------------------------------------------------------------------------Commentary:
1.

Seismic performance of under-ground tanks depend upon seismic stability and the responseof surrounding soil. Examples from recent earthquakes show that the following conditions
contributed significantly to the resulting damage of under-ground tanks:- under-ground structures built in soft ground,- under-ground structures built in ground with abrupt changes in profile or materials,- connecting section with different type of structures or incidental equipment.
Sinking, lifting, inclination and translation of under-ground structures may take place during
earthquakes in soft ground because of significant soil deformation or liquefaction. It istherefore recommended that under-ground tanks should be constructed in stable ground.
Soft ground, or inclined ground which is subject to sliding and slope failure should be
avoided, as should ground with abrupt changes in contours or materials, or at connecting
sections with different structures where significant deformation and stress can be expected.
However, when construction on unstable ground is unavoidable, it is preferable that acomprehensive geological survey is performed, and soil improvement, piled foundations, orother measures are adopted for better seismic safety, (refer to the “Recommendation for thedesign of building foundations” - Architectural Institute of Japan).

2.

Under-ground seismic intensity and soil displacement for the response displacement method
may be calculated using the guidelines and formulae contained in Chapter 3.6.2. Similarly,
seismic earth pressures may be calculated using the guidance and formulae given in Chapter
3.8.2.
As under-ground tanks are bounded by soil, their seismic behaviour will depend mainly on
the displacement and deformation of the surrounding soil.
Whilst the dynamic
characteristics of the structures themselves have less influence on their seismic response thanabove-ground tanks, the properties of the stored materials and their inertia force need to betaken into account.
The response displacement method, seismic intensity method and dynamic analysis are therecommended evaluation methods and as such are explained in the following paragraphs. In
these methods, both the surrounding soil and the structure are analysed and theirdisplacement, deformation, and inertia forces are taken into account simultaneously.

Chapter 8103

It must be noted that utility tunnels, under-ground pipes and linear sealed tunnels etc., orhorizontally long structures are excluded in this chapter. However, such kind of structures,
when associated with under-ground storage tanks, should be examined for deformation inhorizontal section, and reference should be made to the relevant standards andrecommendations.
Figure 8.2.1 below is a flow chart for the three methods,- the response displacement method uses soil springs for the ground model where
ground displacement, circumferential shear force, and inertia force are taken into
consideration,- the seismic intensity method, which normally uses a two dimensional FEM model forthe ground and gives the static under-ground seismic intensity,- and the dynamic analysis which also uses an FEM model for the ground and gives the
dynamic load (earthquake motion).
Response Displacement Method

Seismic Intensity Method

Modeling of Tank and Contents

Modeling of Tank and Contents

Modeling of Tank and Contents

Modeling of Soil (Ground Spring)

Modeling of Soil (FEM and so on)

Modeling of Soil (FEM and so on)

Ground Intensity

Input Seismic Wave

GroundDeformation

Surround
Shear
Force

Inertia
Force

Calculation of Seismic Stress

Dynamic Analysis

Calculation of Seismic StressCalculation of Seismic StressFig. 8.2.1 Method flow chart for seismic analysis methods

The following paragraphs give a conceptual outline of the three evaluation methods.
a) Response displacement
Figure 8.2.2 shows a model for applying the response displacement method to an underground tank where the predominant consideration is the shear deformation in the crosssection of the structure due to shear wave seismic motion of the ground.
The first step involves the calculation of the shear deformation for the free field using
equation (3.16) or the1-D wave propagation theory. The second step involves calculating therelative displacement of the free field to the tank bottom level soil and through to the top of
the tank level soil.
The relative displacement of the free field is applied to the soil springs around the tankstructure as shown in Fig. 8.2.2.
Circumferential shear forces are applied to the outside under-ground walls, and inertia
forces for under-ground seismic intensity, (refer to Chapter 3.6.2.2), are applied to the tankstructure and stored material.

Chapter 8104Refer also to Chapters 4, 5, and 7 for details of the inertia force of stored materials.
The response displacement method has a clear logical basis. The soil springs should
preferably be calculated using FEM model of the ground, but may also be calculated using
simpler methods. Appendix B shows dependence of the stress on soil springs, soil
deformation, circumferential shear force and inertia force.
circumferential shear force

shear
deformationof free fieldinertia force

shear
deformationof free field

circumferential shear force

Fig. 8.2.2 Structural model and forces for the response displacement methodb) Dynamic analysis
Dynamic analysis gives the most reliable result because both the under-ground tank and thesurrounding soil are analysed using an FEM model and the earthquake motion is given as adynamic load. However, it is important that the modeling region, boundary conditions, andthe input earthquake motions should carefully be chosen to ensure that the most reliable
results are obtained. (Refer to Appendix B.)c) Seismic intensity method
The seismic intensity method is similar to the dynamic analysis in that the under-ground tankand the surrounding soil are analysed using an FEM model, but the seismic intensity is given
as static load. Because the method involves a static analysis, the modeling region and theboundary conditions may be different from the dynamic analysis. (Refer to Appendix B.)
Evidence taken from seismic damage shows that the most critical parts for under-groundstructures are where there are abrupt changes in the soil conditions, structural properties, orsections connecting different types of structures. In such cases it is preferable that a detailed
analysis, such as a dynamic analysis, is performed to ensure seismic safety.
3.

In general, the seismic response of the ground is non-linear. Therefore, the soil strain levelshould be taken into consideration together with the appropriate equivalent soil stiffness
when the soil spring constants are evaluated. The equivalent soil stiffness may be calculatedusing the1-D wave propagation theory with the equivalent linear method etc. (Refer toAppendix B). Some references give a simple estimation of the modification factor for groundswith a shear wave velocity of less than 300m/sec. But as an excessive reduction of soil
stiffness reduces the safety margin, a careful estimation of the stiffness reduction is essential.

4.

Because the unit weight of the surrounding soil is greater than that of under-ground tanks orother under-ground structures, they may be subject to lifting when liquefaction of thesurrounding soil occurs. This was evident during the :iigata earthquake in 1964, and the:ihon-Kai-Chubu earthquake in 1983, when underpasses, under-ground oil tanks, watertanks and manholes to under-ground piping were damaged when they were pushed towards

Chapter 8105the surface due to the liquefaction of the surrounding soil. Therefore, due consideration
should be given to this phenomenon when there is a possibility of liquefaction. It must be
noted however, that stored materials should not be taken into consideration when evaluatingthese conditions. Figure 8.2.3 below shows the forces acting on under-ground tanks duringliquefaction.

ground surface

Figure 8.2.3 Forces acting on under-ground tanks during liquefaction.

In the event that liquefaction is possible, QS and QB are neglected. The forces shown are forunit thickness of the under-ground tank.

The safety factor FS for lifting due to liquefaction is calculated using equation (8.2.1) below.
(8.2.1)FS =WS + WB + QS + QB
US +UD

where:

8.3

WS

weight of the overburden (water included) (:/m)

WB
QS

weight of the under-ground tank (stored material and levellingconcrete included) (:/m)
shear strength of the overburden (:/m)

QBfrictional strength of the side walls and soil(:/m)

USbuoyancy on tank bottom caused by hydrostatic pressure (:/m)

UDbuoyancy caused by excess pore water pressure (:/m)

Checking of Seismic Capacity
Stresses in the structure are compared with the design limits. However, sufficientductility is necessary in the structure because the surrounding soil acts to deform all
the members.

Chapter 8106

------------------------------------------------------------------------------------------------------------------------------Commentary:
1.

Working stresses in under-ground tanks are calculated assuming the most severe combinationof seismic loads, the dead load, stored material properties, static earth pressure, hydrostaticpressure, etc., these are then compared with allowable stresses. Stresses that fall within theelastic range of the materials are recommended. The surrounding soil acts to deform thestructural members and therefore no reduction of seismic loads due to any plastic
deformation of the structure can be expected. It is therefore recommended that seismic
designs should have enough ductility to withstand all the stresses indicated above, and shouldtake into account the lessons learned from the seismic damage caused to under-groundstructures with insufficient ductility.

Chapter 8

107

APPEDIX ADESIG AD CALCULATIO EXAMPLES

APPEDIX A1STEEL WATER TOWER TAKS

A1.1

Description of water tower tankA1.1.1

Outline
The structural outline and design conditions are shown in Table A1.1.
Table A1.1 Design Conditions

Description

Outline of structure

Outline of vessel

Design Conditions

Usage
Structural type
Shape
HeightVessel diameter
Diameter of tower
H.W.L
L.W.L
Water capacity
Total volume

Standards and codes

Allowable stress
A1.2

Calculation of Tube Tower

A1.2.1

Assumed section of tube tower

Tower for water supply
Steel plate tower
Vessel – Cylindrical shell
Tower – Cylindrical column
GL + 35.35m
6.8m
1.6~3.4m
GL + 34.35m
GL + 32.0m
80m3
93.4m3a. Building Standard Law
b. Design Recommendation for Storage
Tanks and Their Supports (AIJ)c. Design Standard for Steel Structures(AIJ)
d. Standard for Structural
Calculation of RC Structures (AIJ)e. Japanese Industrial Standards (JIS)Steel – JIS G 3136 SN400B
Concrete – Fc = 21N/mm2
Reinforced bars – SD295
Refer to b. and c.

The natural period of the tower is calculated with the thickness without reducing the
1 mm of corrosion allowance. The outline of the water tower tanks is shown in
Figure A1.1.

Appendix A1

108

Fig. A1.1

A1.2.2Outline of Water Tower Tank

Assumed loads(1) Vertical Load
The total vertical load upon its base plate:

mass

Water

915.32k N (93,400 kg)

Vessel107.80 kN (11,000 kg)

Pipe (1.47kN/m)

47.04 kN (4,800 kg)

Deck (26.46kN for each)

105.84 kN (10,800 kg)

Tube weight107.80 kN (11,000 kg)
1,283.80 kN (131,000 kg)

Appendix A1

109

The sectional dimensions and modeling of the tower are shown in Table A1.2 andFigrure A1.2, respectively.
Table A1.2

Fig.A1.2
A1.2.3

Sections of Tower

Modeling of Tower

Stress Calculation(1) Modified Seismic Coefficient Method
The flexibility matrix and the mass matrix are obtained from the following equation:
Flexibility matrix,[F ] = 

0.0343 0.1228
0.1228 1.0903 

Mass matrix,0m[M ] =  1

0
m
2


Appendix A1

110∴T1 = 1.572 s
T2 = 0.106 s
B = 1.0, Z s = 1.0, I = 1.0, Ds = 0.5, Tg = 0.96,
T1 > Tg
From (3.2)∴ C = Z s ⋅ I ⋅ Ds ⋅Tg
T1

= 0.305 > 0.3Z s I

W = W1 + W 2 = 683.52 kN
Q d = C ⋅ W = 208.47 kN
Qei = Qdi / B

From (3.4)
n∑W ⋅ hmm∑W ⋅ hmf di = Qd mn=i

m

m =1

34,673 kg
56,121 kg

13,626 kg

Fig.A1.3
Appendix A1

Lumped Mass Model111

Table.A1.3

Qei Obtained by Calculations

47(2) Modal Analysis
The flexibility matrix and the mass matrix are obtained from the model shown in Fig.
A1.3(b) as follows:
0.0343 0.1228 0.1265
[F ] = 0.1228 1.0903 1.1415 
0.1265 1.1415 7.2404

δ’33=1/K2 is added to δ 33 from stiffness coefficient, where K2 is the spring
constant for the convective mass.
m1[M ] =  0
 0

0
m2
0

0
0 
m3 

Qd 1 =157.09 kN < 0.3ZsIW =205.02 kN
Stresses in modal analysis are obtained by multiplying withQei =Q di 0.3Z s IW⋅
B
Qd1

Table.A1.4 Natural Periods, u, β and C
thAppendix A1

112

0.3Z s IW
as follows:
Qd 1

Table.A1.5 Qei Obtained by calculation
mass

Table A1.6 Calculation of Qei and Mei

Table.A1.6 Calculation of Qei and Mei

modified seismic coefficient method
modal analysis method

modified lateral force
mode method
story shear

Overturning

bending moment

shear force

Fig.A1.4

Comparision of Stresses obtained using the Modified

Seismic Coefficient Method and the Modal Analysis
A1.3

Water Pressure Imposed on the Tank
Impulsive mass (water),

Q = (Qe 2 − Qe1 )Appendix A1

Wd= 69.39 kN = ∫∫ P ⋅ Rdϕ ⋅ dy
56.21

113

Convective mass (water),
Q' = 89.77 kN = ∫∫ P Rdϕ ⋅ dy
Height of wave due to sloshing; from (7.9),

ηs =0.802 ⋅ Z s ⋅ I ⋅ Sν 1
B

Sν 1 = 2.11 ×

1= 1.758 m/s
1 .2

Fig.A1.5

A1.4

D 3.682h tanh
 = 1.104 mg
 D 

Water Pressure due to Sloshing

Wind Loads(1) Wind loads
Wind loads are obtained from chapter 6 “Wind Loads” of Recommendations for
Loads on Buildings (2004), AIJ.

a)

Design wind speed
Basic wind speed

U 0 = 36 (m/s)

Height Coefficient
E H = Er ⋅ E g = 1.205 × 1.0 = 1.205
From exposure II,

Z b = 5 m, Z G = 350 m, α = 0.15

Appendix A1

114

 35.35 
Er = 1.7 ×  350 0.15

= 1.205

Return period conversion coefficientk rw = 0.63(λu − 1) I nr − 2.9λu + 3.9 = 0.952
Return periodr = 50 years, λu = U 500 / U 0 = 1.111

Design wind speed at standard heightU H = U 0 ⋅ k D ⋅ E H ⋅ k rw = 36 × 1.00 × 1.205 × 0.952 = 41.29 (m/s)b)

Design velocity pressure

qH =

1
1
× ρ × U H2 = × 1.22 × 41.29 2 = 1,040 (N/m2)
2
2

Air Density ρ = 1.22 (kg/m3)
c)

Pressure coefficient
From H / D = 35.35 / 4.6 = 7.68
Z ≦ 5m,H
C D = 1.2 × k1 ⋅ k 2 ⋅ k Z = 1.2 × 0.6 ×  
D

0.14

Z 
× 0.75 ×  b 
H 5 .0 = 1.2 × 0.6 × 7.7 0.14 × 0.75
 35.35 

2α

2× 0.15= 0.400

28.28m > Z > 5m,

 Z 
C D = 1.2 × k1 ⋅ k 2 ⋅ k Z = 1.2 × 0.6 × 7.7 0.14 × 0.75 × 
 35.35 

0.3= 0.247 × Z 0.3

Z≧28.28m,
C D = 1.2 × k1 ⋅ k 2 ⋅ k Z = 1.2 × 0.6 × 7.7 0.14 × 0.75 × 0.80.3 = 0.672
d)

Gust effect factor
exponent for mode shape β

2.500

total building mass

104420kg

M

generalized building mass

MD

λ

0.633

mode correction factor φD

0.174

turbulence scale at reference heightLH

turbulence intensity at reference height IrHAppendix A1

84644kg

115

108.6m0.158

topography factor for the standard deviationof fluctuating wind speed EI

1.000

topography factor for turbulence intensity EgI

1.000

I H = IrH ⋅ E gI

0.158

gust effect factor
natural frequency for the first mode inalong-wind direction

fD

0.713Hz

damping ratio for the first mode ine)along-wind direction ζ D

0.010

Cg

0.457

k

0.070

C’g

0.127

R

0.386

F

0.075

SD

0.404

FD

0.033

RD

2.593

vD

0.606

gD

3.604

gust effect factor GD

2.044

Wind load

WD = q H ⋅ C D ⋅ GD ⋅ A
Z≦5m

WD = 1,040 × 0.400 × 2.044 × A = 850 × A
28.28ｍ＞Z＞5m
WD = 1,040 × 0.247 × Z 0.3 × 2.044 × A = 0.525 × Z 0.3 × A
Z≧28.28m

WD = 1,040 × 0.672 × 2.044 × A = 1,429 × A(2) Wind load for vortex oscillationa)

Calculation model
In this example, the calculation model is assumed to have the shape shown in
Fig. A1.6, and the water content is assumed to be empty.

Appendix A1116Natural period
From

T1,

[F ] = 

0.0343 0.1111

 0.1111 0.7952

T1 = 0.726 s[M ] =  1
m

0

0 
m'2 

16,156 kg

13,626 kg

Fig.A1.6

Calculation Model

Table.A1.7 Sectional Forces by Wind Load in the Along Direction
qH ⋅ CD ⋅ GD

1.429

33.04

1.429

35.87

72.79

1.429

39.21

137.35

1.428

44.35

172.64

1.371

52.25

332.30

1.307

59.78

520.40

1.235

66.89

735.61

1.151

74.74

976.41

1.022

84.86

1312.74

0.850

95.57

1694.610.850

101.782124.68b) Calculation of wind load for vortex oscillation
The possibility of vortex oscillation is checked according to Chapter 6 “WindLoad” of the “Recommendations for Loads on Buildings (2004)”, AIJ as follows.

Appendix A1

117

Wr = 0.8 ⋅ ρU r2Cr

Z
A
H

in which,

U r = 5 ⋅ f L ⋅ Dm = 5 × 1.38 × 2.2 = 15.18 (m/s) : resonance wind speed
where,

f L = 1.38 : natural frequency of water tower tank (Hz)

ρ = 1.22 (kg/m3) : air density
H = 35.3 (m) : tower height
Dm = 2.20 (m) : diameter at the level 2/3 H
Cr = 4.03 : pressure coefficient in resonanceU r ⋅ Dm = 15.18 × 2.2 = 33.4 > 6.0 
0.57= 4.03 then, Cr =
ρ s ζ L = 142 × 0.02 = 20.1 > 5.0ζL

ζ L = 0.02 : damping ratio of water tower tankρs =M
37,600== 142 (kg/m3)
H ⋅ Dm ⋅ DB 35.3 × 2.2 × 3.4

M = 37,600 (kg)Then,
Wr = 0.8 ρU r2CrZ
Z
A = 0.8 × 1.22 × 15.182 × 4.03 ×
× A = 25.7 × Z × A (N)H
35.3

The cross section of the tower should be designed with the composite force ofobtained from the wind load itself and the wind load for vortex oscillation

Appendix A1

118

Table.A1.8 Forces Caused by Wind Load Vortex Oscillation
0.8 ρU r2

Section forces in windalong direction

Composite forces

0

33.04

39.41

43.20

35.87

72.79

41.66

84.64

81.32

39.21

137.3545.42

159.62101.9544.35

172.6450.81

200.50191.2352.25

332.3059.48

383.40293.5859.78

520.4067.58

597.50407.0566.89

735.6175.07

840.72529.7474.74

976.4183.17

1110.86693.9084.86

1312.7493.32

1484.85868.5995.57

1694.61

103.731904.25

1050.12101.782124.68

109.502370.02

Cr
H

19.61

0.0257

19.61
2.203

0.0257

1.57

21.18

0.0257

1.74

22.92

1.80
0.90
1.88

0.025724.80
3.60

0.0257

3.63

28.43
3.60

0.0257

3.09

31.52
3.60

2.56

0.025734.08
3.60

0.0257

2.40

36.48
4.50

2.34

0.025738.82
4.50

0.0257

1.52

40.34
4.50

0.04

0.025740.38

Table.A1.9 Forces List
Seismic(modified seismic
coefficient)

Axial
force
Full

Seismic(modal analysis)Wind inalong dir.

Empty

35.87
39.21
44.35
52.25
59.78
66.89
74.74
84.86
95.57
101.78

Appendix A1Wind for vortex
oscillation

119

72.79
137.35172.64332.30520.40735.61
976.411312.74
1694.612124.68

21.18
43.20
22.92
81.32
24.80 101.95
28.43 191.23
31.52 293.58
34.08 407.05
36.48 529.74
38.82 693.90
40.34 868.59
40.38 1050.12

Table.A1.10 Section Properties and Allowable Stresses
Long term

Table.A1.11

Design of Section

Full water tank (modified seismic coefficient method)

Full water tank (modal analysis)

Appendix A1

120

Seismic

Wind with full water tank (in along direction)

72.79
137.35172.64332.30520.40

35.87
39.21
44.35
52.25
59.78

0.203
0.869
1.092
1.538
2.045

735.61
976.411312.74
1694.612124.68

66.89
74.74
84.86
95.57
101.78

2.891
3.338
2.709
2.150
1.823

0.010
0.040
0.050
0.067
0.088
0.124
0.142
0.121
0.100
0.089

0.102
0.163
0.173
0.153
0.161
0.199
0.207
0.182
0.153
0.137

0.119
0.196
0.222
0.191
0.185
0.207
0.200
0.190
0.168
0.147

0.025
0.024
0.027
0.019
0.017
0.019
0.018
0.021
0.022
0.024

Wind with empty water tank (for vortex osc.)

84.64
159.62200.50383.40597.50

41.66
45.42
50.81
59.48
67.58

0.237
1.010
1.269
1.774
2.348

0.012
0.046
0.058
0.078
0.101

0.032
0.090
0.108
0.135
0.168

0.139
0.227
0.254
0.217
0.209

840.72 75.07
1110.86 83.17
1484.85 93.32
1904.25 103.73
2370.02 109.50

3.304
3.798
3.064
2.416
2.034

0.142
0.161
0.137
0.113
0.099

0.233
0.261
0.227
0.189
0.170

0.232
0.223
0.209
0.182
0.158

OK
OK

Appendix A1

121

0.029
0.028
0.031
0.021
0.019
0.022
0.020
0.023
0.024
0.026

0.043
0.042
0.047
0.032
0.029
0.032
0.030
0.034
0.036
0.039

APPEDIX A 2
DESIG EXAMPLESTEEL SILO

The equations referred in this chapter correspond to those given in chapters 3 and 5 of the main
body of this Recommendation.
A2.1

Specifications
Table A2.1 Specifications

Application

Silo for storing grain

Type

Welded steel silo

Configuration

A single free-standing, cylindrical shape

Height34,500 mm

Diameter

8,600 mmAngle of hopper :α45°

Thicknesses of wall, skirt and hopper

t1 = 6～15.2 mm, t2 = 6～8mm

Stored product

Grain

Storage weight

1176 × 104N

Density

For long-term design

7.35 × 10-6N/mm3

For seismic design

80% of the above

Angle of repose:

30°

Internal friction angle:

25°

Friction coefficient of tank wall:

0.3

Impact pressure coefficient:

1.5

Dynamic pressure

Wall

1.92

coefficient:

Hopper2.30(top), 1.0(bottom)

Safety factor of friction force:
Design codes

1.5- Design Recommendation for Storage Tanks andTheir Supports, AIJ- Design Standard for Steel Structures- Japanese Industrial Standards

Material

SS400(JIS G 3101)

Allowable stress

F=235N/mm2

Appendix A2

122Appendix A2
1,900

15.2

4,300

6,200

12.0

4,300
11.0

6,000

123
H=34,500

hm=30,411

26,111

10.0

6,000
8.0

6,000
6.0

6,000

2,189

8,600

Fig. A2.1 Specifications of Calculated Silo (Unit:mm)

A2.2

Stress Calculation(1) Pressures and Forces Caused by the Stored Materiala) Internal pressure in silo and axial force on tank wall:
Vertical pressure: PυPυ =From (5.1),− µ f・K ・
γ ⋅ rw s x / rw ⋅ 1 − eµ f ⋅ Ks where:
bulk density of the stored material, γ ,=internal friction angle of the stored material, ϕ i ,=radius of water pressure, rw , [as determined from (5.3)] =From (5.3),

7.35 × 10−6 (N/mm3)
25 (°)
2,150 (mm)rw = A/lc = d/4 = 8,600/4 = 2,150mm

depth from the free surface of the stored material, x , (mm)friction coefficient of the tank wall, μf ,= 0.3
Ks =From (5.2),

1 − sin ϕ i 1 − sin 25°== 0.401 + sin ϕ i 1 + sin 25°(∴

Pυ = 0.13169 1 − e −0.000558 x∴
where:dPυ =Ci =dPυ = Ci ⋅ Pυ = 1.5 ⋅ Pυ)

design vertical pressure per unit area
impact pressure factor when filling the silo with grain = 1.5
Cd = 1.92

xx=0=3,811hi=26,111

1.92

⑤
④

d=8,600

h0=4,300

hm=30,411=9,811hm / d =1.92=15,811

1.92=21,8112.30=26,111=27,611ⓐ①=29,111ⓑ

α = 45

o

1.00

③
②

30,4118,600= 3.54↓

Cd = 1.92
C 'd = 1.2Cd = 2.30 (hopper top)
C 'd = 1.0 (hopper bottom)

d2

Fig.A2.2 Minimum Value of Cd and C'dHorizontal pressure: Ph (N/mm2)
From (5.4),

Appendix A2

Ph = K s ⋅ Pυ = 0.40 ⋅ Pυ

124

∴

dPh = Cd ⋅ Ph = 1.92 ⋅ Ph
C 'd = 1.2 × Cd = 2.30 (top)for hopper

C 'd = 1.00 (bottom)
where:
dPh
design horizontal pressure per unit area on the tank wall (N/mm2)
Cd , C 'd dynamic pressure coefficient when discharging grain from the silo 1.00 ~2.3 (defined by Fig. A2.2)
Friction force per unit area on the tank wall: Pf (N/mm2)
From (5.5),

Pf = µ f ⋅ Ph = 0.3 ⋅ Phb) Normal pressure on hopper wall: Pa (N/mm2)
From (5.6),

Pa = Ph ⋅ sin 2 α + Pυ ⋅ cos 2 α

where:
αincline of the hopper wall to the horizontal plane 45 (°)∴
From (5.9),

Pa = 0.7 Pυ

dPa = C'd ⋅Pa = 1.61Pυ ~ 0.7Pυ (QC'd = 2.30 ~ 1.00)
dPa = C i ⋅ Pa = 1 .05 ⋅ Pυ (Q C i = 1 .5 )where:
dPa

design vertical pressure per unit area on the hopper wall (N/mm2)

Table A2.2(a) Internal Pressure in Silo and Axial Force on Tank Wall
VerticalPressure(N/mm2)

Horizontal Pressure (N/mm2)

Axial Force on Tank
Wall (N/mm)

x
(mm)

Pυ

dPυ

Ph

Cd

dPh

Long
Term
dm

DuringEarthquakes
dm

⑤

3,811

0.025

0.038

0.010

1.92

0.020

9.1

6.1

④

9,811

0.055

0.083

0.022

1.92

0.043

54.1

36.1

③

15,811

0.077

0.115

0.031

1.92

0.060

127.2

84.8

②

21,811

0.092

0.138

0.037

1.92

0.072

220.2

146.8

①

26,111

0.100

0.150

0.041

1.92

0.078

295.8

197.2

Calculation
Point

Dept
h

Appendix A2

125Calculation
Point

(b)

Internal Pressure and Axial Force in Hopper

Depth

Pressure

Pressure Factor

Design Pressure

x
(mm)

Ph
Pa
Pυ
2
2
(N/mm2) (N/mm ) (N/mm )

Ci
(－)

C’d
(－)

dPa
(N/mm)

dPυ
(N/mm)

top

26,111

0.100

0.041

0.071

1.50

2.30

0.163

0.150

middle27,544

0.102

0.041

0.072

1.50

1.87

0.135

0.153

bottom

28,977

0.105

0.042

0.074

1.50

1.44

0.111

0.158

c) Design generator force acting on cross section of tank wall caused by wall friction
force: dm (N/mm)
Long-term design generator force per unit length of the circumference of the tankwall, dm (N/mm).
Safety factor of the friction force between the tank wall and the granular material(long term), Cf , = 1.5 or 1.0. Then, from equation (5.11),()

d m = (γ ⋅ χ − Pυ ) ⋅ rw ⋅ C f = 7.35 × 10 −6 χ − Pυ ⋅ 2150 ⋅ (1.5)
Design generator force per unit length of the circumference of the tank wall during
earthquakes, d m ' (N/mm)()

d m ' = (γ ⋅ χ − Pυ ) ⋅ rw = 7.35 × 10 −6 χ − Pυ ⋅ 2150

Appendix A2

126

(2)

Calculation of Design Axial Force and Weighta)Design axial force = (empty-weight) + (axial force of stored material caused bywall friction force)

Table A2.3 Design Axial Force
Calculation
Point

Axial Force of Stored Material Caused
by Wall Friction Force: dm (N/mm)

Empty-weight

Cumulated
Weight

⑥

12.7

12.7

⑤

3.2

15.9

9.1

④

4.2

20.1

③

4.9

②

Empty-weight (N/mm)

Long Term

DuringEarthquakes--

Design Axial Force dm(N/mm)Long Term

DuringEarthquakes

12.7

12.7

6.1

25.0

22.0

54.1

36.1

74.2

56.2

25.0

127.2

84.8

152.2

109.8

4.5

29.5

220.2

146.8

249.7

176.3

①

9.2

38.7

295.8

197.2

334.5

235.9

⓪

1.5

40.2

475.5

388.4

-

-

b) Calculation of Weight

Table A2.4 Basic Design Weight and Weight for Calculating Seismic Force
Basic Design WeightW (kN)

Weight for Calculating
Seismic Force
Wt (kN)

SectionSectionCalculation
Point

Weight of Stored Material(kN)

Empty-weight
(kN)

⑥

344.0

346.2

277.0

690.2

690.2

621.0

621.0

⑤

87.4

2,561.8

2,049.4

2,649.2

3,339.4

2,136.8

2,757.8④

112.4

2,561.8

2,049.4

2,674.2

6,013.6

2,161.8

4,919.6③

131.1

2,561.8

2,049.4

2,692.9

8,706.5

2,180.5

7,100.1②

122.4

2,198.9

1,759.1

2,321.3

11,027.8

1,881.5

8,981.6①

290.9

1,529.9

1,223.9

1,820.8

12,848.6

1,514.8

10,496.4

γ= 0.735 x 10-2

γ = 0.735 x
10-2 x 0.8(See Fig. A2.5)

Appendix A2

127

Cumulation

Cumulation

(3)Evaluation of Wall Stress in Circumferential Directiona) Fumigation Pressure with Long Term Design Pressureσθ =(dPh + Pi ) ⋅ d≦ ft
2 ⋅ tiwhere :σ θ : silo wall tensile stress in circumferential Direction (N/mm2)
dPh : horizontal pressure on tank wall (N/mm2)
Pi : pressure other than dPh (e. g., fumigation pressure) =500(mmAq)=0.0049(N/mm2)

ti : thickness of tank wall (mm)f t : allowable tensile stress = 157 (N/mm2)
Table A2.5 Evaluation of Wall Stress in Circumferential Direction(fumigation pressure considered)Calculation
Point
6
5
4
3
2
1
b)

ti
(mm)
6.0
8.0
10.0
11.0
12.0
15.2

d
(mm)
8,600
8,600
8,600
8,600
8,600
8,600

dPh
(N/mm2)
0.020
0.043
0.060
0.072
0.078
0.000

Pi
(N/mm2)0.00490.00490.00490.00490.00490.0000σθ2(N/mm )
17.8
25.7
27.9
30.1
29.7
0.0

ft
(N/mm2)
157
157
157
157
157
157

Local Pressure with dPh
0.1d
P

dPL
P

P

Fig. A2.3

Appendix A2

Local Pressure

128

σθ
ft
0.11
0.16
0.18
0.19
0.19
0.00

Test
Result
≦1.0
≦1.0
≦1.0
≦1.0
≦1.0
≦1.0

σ t + σ b ≦ ft
dPh ⋅ d2 ⋅ tiM
σb = bZ
M 0 = 0.318 ⋅ P ⋅ (d / 2)σt =Mb = k ⋅ M0

(N ⋅ mm/mm)(N ⋅ mm/mm)

k = 0.1 + 0.05 ⋅ (hm / d − 2) , however, if hm / d < 2 , k = 0.1 (−)=0.177
P = 0.1 ⋅ d ⋅ dPL ( N/mm)
dPL = C L ⋅ Cd ⋅ Ph

( N/mm 2 )C L = 0.15 + 0.5 ⋅ (e / d ) (−)

See Fig. A2.4

dd⑥

C L ＝ 0.04⑤

C L ＝ 0.11④

C L ＝ 0.15

hm

hm
1.5d

③

C L ＝ 0.15②

C L ＝ 0.15

Fig. A2.4

C L Value

where :

σ t : tensile stress (N/mm2)
σ b : bending stress by the local pressure (N/mm2)
M b : design bending moment by the local pressure ( N ⋅ mm/mm)
M 0 : maximum bending moment for ring model ( N ⋅ mm/mm)
k : local stress coefficient (－)

Appendix A2

129

hm : effective height of material = 30,411 (mm)P : local load (N/mm)
dPL : design local horizontal pressure per unit area (N/mm2)
dPL = C L ⋅ Cd ⋅ Ph (N/mm2)
C L : ratio of the local pressure to the design horizontal pressure (－)
e : eccentricity = 0 (mm)Table A2.6 Evaluation of Wall Stress in Circumferential Direction(local pressure considered)Calculation
Point
6
5
4
3
2
1

ti
(mm)
6.0
8.0
10.0
11.0
12.0
15.2

σt

dPh
(N/mm2)
0.020
0.043
0.060
0.072
0.078
0.0

d
(mm)
8,600
8,600
8,600
8,600
8,600
8,600

2

(N/mm )
14.3
23.1
25.8
28.1
28.0
0.0

ft
(N/mm2)
157
157
157
157
157
157

Table A2.7-1 Evaluation of Wall Stress in Circumferential Direction(local pressure considered)Calculation
Point

ti
(mm)

d
(mm)

dPh
(N/mm2)

6
5
4
3
2
1

6.0
8.0
10.0
11.0
12.0
15.2

8,600
8,600
8,600
8,600
8,600
8,600

0.020
0.043
0.060
0.072
0.078
0.0

Appendix A2

130

CL

dPL
(N/mm2)P(kN/m)

0.04
0.11
0.15
0.15
0.15

0.00080.00470.00900.01080.0117

0.69
4.04
7.74
9.29
10.06

Table A2.7-2 Evaluation of Wall Stress in Circumferential Direction(local pressure considered)

6
5
4
3
2
1
(4)

σb

Z

Mb
(kNm/m)
0.17
0.98
1.87
2.25
2.43
0.0

Calculation
Point

3

(mm /m)6,000.0
10,666.7
16,666.7
20,166.7
24,000.0
38,506.7

σt +σbσt2(N/mm )
28.3
91.9
112.2
111.6
101.3
0.0

2

(N/mm )
14.3
23.1
25.8
28.1
28.0
0.0

ft
0.27
0.73
0.88
0.89
0.83
0.00

Test
Result
≦1.0
≦1.0
≦1.0
≦1.0
≦1.0
≦1.0

Calculation of Stress during Earthquakesa) Modified Sesmic Intensity Method･Design shear force at the silo base: Qd (N)
Qd=C･Wt
C = Z s ⋅ I ⋅ Ds ⋅C ≧ 0.3･Zs･I

where
Symbols

S a1
gC: design yield shear force coefficient at the silo baseWt:

design weight imposed on the silo base = 10,496.4 × 103 (N)Zs:

seismic zone factor = 1.0 (Kanto area)I:

importance factor = 0.8

Ds

:

Dη

:

= DηDh = (0.5～0.7)･1.0 = 0.5～0.7
a constant determined by the plastic deformation of the silo = 0.5～0.7(steel structure)Dh:a constant determined by the damping capacity of the silo = 1.0 (area ofsilo basement S ＜ 2,000 m2)
acceleration response spectrum of the first natural period = 9.80 (m/s2)

Sa1 :g = 9.80 (m/ s2)
C=1.0･0.8･(0.5～0.7)･9.80/9.80=(0.4～0.56)>0.3･Zs･I=0.3･1.0･0.8=0.24
Qd = (0.4～0.56 ) ⋅ 10496.4 × 103 = (4198.56～5877.98) × 103･ Design yield shear force at each part of the structure:n∑W ⋅ hmm

Qdi = Qd ⋅ mn=i

∑W ⋅ hmm

m=1

Appendix A2131

Qdi (N)(N)

Symbols

Qdi : design horizontal yield shear force at the i-th mass when the structure isassumed to be a vibration system with n masses (N)Wi : weight of the i-th mass (N)hi :height of the i-th mass from the ground (m)n :total number of masses

Fig. A2.5･ Design shear force:Qei (N)

Qei=Qdi / B
SymbolB: ratio of the horizontal load-carrying capacity of the structure to theallowable strength during earthquakes

1.0

Table A2.8 was the calculated results of Qei.

Calculation
Point

Table A2.8 Design Shear Force/Bending Moment during Earthquakes Calculated by
Modified Seismic Intensity Method
HeightSection Height of Weight of
Massfrom theLength the i-th the i-th W = W
Point
∑ i
GroundPoint
Point
Number
(kN)
h
hs
hi
Wi
i
(m)
(m)
(m)
(kN)

⑥

34.50

6

⑤

28.50

5

④

22.50

③

n

Wi × hi∑ Wi × hi(kN-m)(kN-m)∑W ⋅ h
m

m

m =i
n

∑W ⋅ h
m

m

Qd
(kN)

Qdi
(kN)

Qei
(kN)

Pi
(kN)

Me
× 10 5 ( N ⋅ m )

m =1

( －)

34.50

621.0

621.0 21,424.5 21,424.5

0.11

461.8

6.00

28.50

2,136.8 2,757.8 60,898.8 82,323.3

0.42

1,763.4 461.8

4

6.00

22.50

2,161.8 4,919.6 48,640.5 130,963.8

0.67

2,813.1 1,763.4 1,301.6 133.5

16.50

3

6.00

16.50

2,180.5 7,100.1 35,978.3 166,942.1

0.85

3,568.8 2,813.1 1,049.7 302.3

②

10.50

2

6.00

10.50

1,881.5 8,981.6 19,755.8 186,697.9

0.95

3,988.7 3,568.8 755.7

516.4

①

6.20

1

4.30

6.20

1,514.8 10,496.4 9,391.8 196,089.7

1.00

4,198.6 4,198.6 3,988.7 419.9

687.9

⓪

0.00

4,198.6 209.9

948.3

*

6.20

461.8

27.7

The axial compressive stress ratio for equations (3.7.12) and (3.7.13) is σc /cfcr ＜ 0.2, which gives 0.5 to Ds.
3

Therefore, Qd = 4,198.6 x 10 (N) is used.

Appendix A2

132

･ Bending Moment:

Mei (N･m)

When we use the model shown in Fig. A 2.5 the bending moment can be calculated as follows:
The seismic force Pi acting on each part of the structure (each mass point) is calculated fromthe following equation using the design seismic-story shear force Qei described above:
Pi = Qei − Q e,i −1 (i = 1～n )The moment of each story(calculation point) can be calculated as follows:n

M e,i −1 = ∑ (Pm ⋅ H m ) (i = 1～n )m =iwhere the subscripts m of P, H and i of Qei, Mei mean the mass and the calculation point,respectively. Hm shows the height from the calculation point to each mass.
b) Modal Analysis･ Design yield shear force at each part of the structure:
Qdi =

n


∑  C j ⋅ ∑ Wm ⋅ β j ⋅ umj 
j =1 
m =i
k

Qdi (N)2

Cj=Zs･I･Ds･Saj /g

where
Symbols

Qdi ≧ 0.3･Zs･I･W

βj :participation factor of the j-th natural vibrationn∑ Wi ⋅ uij

β j = in=1

∑ Wi ⋅ uij 2
i =1

Symbolsuij: the j-th eigenfunction of the i-th massk: maximum number of vibration mode having a significant effect on
the seismic response = 3

j: order of natural frequency

Cj: design yield shear force coefficient for the j-th natural mode

Saj

: design acceleration response spectrum for the j-th natural period
defined in 3.6.1.6. (m/s2)

Zs :seismic zone factor = 1.0 (Kanto area)

I :essential facility factor = 0.8

Ds = Dη･Dh

SymbolsDη: a constant determined by the ductility of the structure = 0.5～0.7(steel structure)

Appendix A2

133

Dh

: a constant determined by the damping of the structure = 1.0 (areaof silo foundation S ＜ 2,000 m2)∴

Ds=(0.5～0･7)･1.0=0.5～0.7
Tj(max)=0.559<TG=1.28s

Symbols

Tj: the j-th natural period (s)TG: critical period determined by the ground class (s)Saj = 9.80 m/s2
C j = Z s ⋅ I ⋅ Ds ⋅･ Design shear force:

S ajg= 1.0 ⋅ 0.8 ⋅ (0.5～0.7 ) ⋅9.80= (0.4～0.56 ) > 0.3 ⋅ 1.0 ⋅ 0.8 = 0.24
9.80

Qei (N)

Qei=Qdi /B･ Bending Moment: Mei (N･m)
Pi=Qei－ Qe,i-1 (i =1～n)n

M e,i −1 = ∑ (Pm ⋅ H m ) (i = 1～n )m =iTable A2.9 Natural Period, Mode, and Participation Factor
j

1

2

3

i
6
5
4
3
2
1
6
5
4
3
2
1
6
5
4
3
2
1

Wi
(kN)
621.0
2,136.8
2,161.8
2,180.5
1,881.5
1,514.8621.0
2,136.8
2,161.8
2,180.5
1,881.5
1,514.8621.0
2,136.8
2,161.8
2,180.5
1,881.5
1,514.8

Cj

Tj

βj

0.40

0.511

1.606

0.40

0.125

0.894

0.40

0.064

0.464

umj

βj･ umj

1.00000.77760.54730.33350.15640.0651-1.0000-0.30070.36940.6711
0.55880.31711.0000-0.1279-0.5866
0.06700.63490.50391.60601.24880.87900.53560.25120.1046-0.8940-0.26880.33020.60000.4996
0.28350.4640-0.0593-0.27220.03110.29460.2338

Cj･Wi･βj･umj

Qij(Qij)2(kN)
398.93
1,067.37
760.09467.15189.0563.38-222.07-229.75
285.53532.32376.00171.78115.26-50.68-235.3827.13
221.72141.66(kN)
398.93
1,466.30
2,226.39
2,693.54
2,882.59
2,945.97-222.07-451.82-166.29357.03733.03904.81115.2664.58-170.80-143.67
78.05
219.71(109N2)159.1
2,150.0
4,956.8
7,255.2
8,309.3
8,678.7
49.3
204.1
27.7
127.5
537.3
818.7
13.3
4.2
29.2
20.6
6.1
48.3

* The axial compressive stress ratio for equations (3.7.12) and (3.7.13) is σc /c fcr＜ 0.2, which gives Ds= 0.5 and Cj = 0.4.

Appendix A2

134

Table A2.10 Design Shear Force/Bending Moment during Earthquakes Calculated byModal Analysis
hi (m)i

Qdi (kN)

⑥

34.5

6

470.9

⑤

28.5

5

④

22.5

③

Mei

Qei (kN)

Pi (kN)hs (m)

1,535.7

470.9

470.9

6.00

28.3

4

2,239.1

1,535.7

1,064.8

6.00

120.4

16.5

3

2,720.9

2,239.1

703.4

6.00

254.7

②

10.5

2

2,975.3

2,720.9

481.8

6.00

418.0

①

6.2

1

3,089.6

2,975.3

254.4

4.30

545.9

⓪

0.00

0

3,089.6

114.3

6.20

737.5

(5)

B

1.0

× 10 5 ( N ⋅ m)

Calculation of Stress under Wind Pressure

The construction site is assumed to be Tsurumi Port in the City of Yokohama. Considering the
ground surface conditions of the surrounding area ( inland port where mid-rise buildings lie
scattered ), the roughness of the ground surface is estimated to be class II.

For load calculation, we calculate the horizontal load for the structural frame. The design return
period (r) is assumed to be 50 years.

i) Calculation of Design Wind Velocity･ Basic wind velocity

UH

U0 = 38 m/s (Tsurumi Port, Yokohama City)･ Vertical distribution coefficient of wind velocity

E = Er･Eg

Zb = 5 m, ZG= 350 m, α= 0.15 (because the roughness of the ground surface is classII)
Er = 1.7(Z/ZG)α = 1.7 x (34.5/350)0.15=1.20
Eg = 1.0 (topography factor)∴

E = 1.20 x 1.0 = 1.20
･ Return period conversion factor of the wind velocity krw
k rw = 0.63(λU − 1) ln r − 2.9λU + 3.9
U
λU = 500
U0

Appendix A2

135

U500 : 500 - year - recurrence 10 - minute mean wind speed at 10m above ground over a flatand open terrain = 44 (m/s)
k rw = 0.93∴

Therefore, design wind velocity UH and design velocity pressure qH are
UH = U 0 K D E H k rw
= 42.4 (m/s)
qH = 0.5 ρU H2= 1,096.6 (N/m2)

Wind force coefficient for circular-planned structures
C D = 1.2k1k 2 k Z
k1 : factor for aspect ratioH= 4.01D∴
k1 = 0.73

CDk 2 : factor for surface roughness = 0.75
k Z : factor for vertical profilek Z = ( Z b / H ) 2α
= (Z / H )= 0.8

Z ≦ Zb2α2α

Z b＜Z ≦ 0.8 H

0 .8 H ≦ Z
Table A2.11 Calculation of Wind Force Coefficient

Calculation
Point

Height from the
Ground (m)

6
5
4
3
2
1

34.5
28.5
22.5
16.5
10.5
6.2

k1

(―)

0.73

Gust Effect Factor GD
GD = 1 + g D

C 'g
Cg

1 + φ 2 D RD

g D = 2 1n (600ν D ) + 1.2

Cg =Appendix A2

1
1+ = 0.457
3 + 3α 6

136

k2
(―)

kZ
(―)

0.75

1.00
0.94
0.88
0.80
0.70
0.60

Wind Force
Coefficient
CD
(―)
0.66
0.62
0.58
0.53
0.46
0.39

0.49 − 0.14α

C 'g = 2 I H

 0.63( BH / LH )0.56 
1 +
( H / B)
0.07 H / B ≧1
k=
 0.15 H / B＜1

φD ＝1.1－0.1β
β：exponent for mode = 1∴

φD ＝1.0RD
1 + RDν D = fD
RD =

πFD
4ζ D

FD =

I H2 FS D (0.57 − 0.35α + 2 R 0.053 + 0.042α
C '2g
1

R=

1 + 20F=

fDB
UH
4 f D LH
UH
5/6

2

 f D LH  
1 + 71 
 U H  


SD =

0 .9
0.5

2

 f D H   
fDB   1 + 3
1 + 6 
U H 

 U H   

φ D : mode correction factor (－)
f D : natural frequency for the first mode in along-wind direction = 1.1 (Hz)ζ D : damping ratio for the first mode in along wind direction = 0.01 (－)
H : reference height =34.5 (m)B : width of building =8.60 (m)U H : design wind speed =38.0 (m/s)
I H : turbulence intensity at reference height (－)

LH : turbulence scale at reference height (－)α : exponent of power law in wind speed profile = 0.15 (－)

Appendix A2

137

0.9

∴ SD =

0.52

1.1 × 8.6 1.1 × 34.5   1+6

 38   1 + 3 × 38   
H ∴ LH = 100 ×   = 115.0(−)
 30 
4 × 1.1 × 115.0
38.0
∴F =
= 0.051
5/ 62
1.1 × 115.0  
1 + 71 
 38.0  

1
∴R =
= 0.167
1.1 × 8.6
1 + 20
38.0

FD =

(

= 0.195

I H2 FS D 0.57 − 0.35α + 2 R 0.053 + 0.042α
C '2g

)

I H = I rz E g1
I rz : turbulence intensity at height H on the flat terrain categories

E gI : topography factor  H −α −0.05
0.1
  Z G 
I rz 
−α − 0.05
  Zb 0.
1 Z 
  G

Zb < H ≦ ZG
H ≦ Zb∴ I rz = 0.159
E
E gI = I
Eg

E I : topography factor for the standard deviation of fluctuating wind speed = 1

E g : topography factor for mean wind speed = 1∴ E gI = 1∴ I H = 0.159

C 'g = 2 I H

0.49 − 0.14α= 0.0756
 0.63( BH / LH )0.56 
1 +
( H / B)k



(

)

∴ FD =I H2 FS D 0.57 − 0.35α + 2 R 0.053 + 0.042α= 0.026
C '2g∴ RD =πFD= 2.082
4ζ D

Appendix A2138∴ν D = f DRD= 0.903
1 + RDg D = 2 ln(600ν D ) + 1.2 = 3.71
Therefore, gust effect factor G D is defined as follows.
C'g
GD = 1 + g D
1 + φ 2 D RD = 2.08
Cg
Along-wind Load WD
Table A2.12 Calculation of Along-wind Load WD

Calculation
Point

Heightfrom theGroundSurface
(m)

6
5
4
3
2
1

34.5
28.5
22.5
16.5
10.5
6.2

DesignVelocity
PressureqH(N/m2)

Wind
Force
CoefficientCD(－ )1,096.6

0.66
0.62
0.58
0.53
0.46
0.39

Gust
EffectFactorGD(－ )Projected
Area
A
(m2)

Along-Wind
Load
WD
(kN)

Wind
Shear-Force
Q
(kN)

1.52

51.67
51.70
51.72
51.73
37.08
53.51

56.84
53.43
50.00
45.70
28.43
34.78

56.84
110.27160.27205.97234.40269.18

Table A2.13 Calculation of Overturning Moment (Mw)i
Calculation
Point
⑥
⑤
④
③
②
①
⓪

Appendix A2

Height from the
Ground Surface
hi
(m)
34.5
28.5
22.5
16.5
10.5
6.2
0.0

Section Lengthhr(m)Projected Area
Aw
(m2)

Along-Wind
Load
WD
(kN)

Moment(Mw)i(kN･m)

6.0
6.0
6.0
6.0
4.3
6.2

51.67
51.70
51.72
51.73
37.08
53.51

56.84
53.43
50.00
45.70
28.43
34.78

170.50671.901,483.50
2,582.20
3,529.00
5,090.10

139

A2.3

Evaluation

A2.3.1

Evaluation of cylindrical buckling is made by the following equations:From(3.21)

σc
f cr

c

and

τ

(3.22)
s

f cr

+

σb

≤1

f cr

b

≤1

Notations :

σc
σb
τ
W
M
Q
A
Z
c fcr
b fcr
s fcr
where:
r
t
v
E
F
r/t

average compressive stress =W / A (N/mm2)
compressive bending stress =M / Z (N/ mm2)
shear stress=
2 Q / A (N/ mm2)
compressive force (N)bending moment (N･mm)
shear force (N)cross section (mm2)section modulus (mm3)allowable compressive stress (N/mm2)
allowable bending stress (N/mm2)
allowable shear stress (N/mm2)

=
=
=
=
=
=

d/2 = 4,300 (mm)6.0 ~ 15.2 (mm)0.3
2,058 x 102 (N/mm2)2.35 x 102 (N/mm2)
282.9 ~ 716.7

σh =Ph ⋅ dP= 4,300 ⋅ h
2⋅tt

SS 400(N/mm2)For the calculation of the shear stress the following values are used:
l/r = 1.0
σh = 0A2.3.1.1 The allowable long-term stress is given by the following equations:(1)-

Appendix A2Allowable compressive stress,

σh(N/mm2)≤ 0.3 → 0 ≤ σ h ≤ 70.5

In the case of0≤

From (3.23),

c f cr = c f cr +

F

c fcr

0.7 ⋅ f cr 0 − c f cr  σ h 
・

0.3
 F 

140

-

In the case of
From (3.24),

σhF> 0.3 → σ h > 70.5
c f cr = f cr 0 ⋅ 1 −

σh 

F 



where,
average tensile hoop stress caused by the internal pressure, σh , (N/mm2)
standard compressive buckling stress, fcro , (N/mm2)
c fcr

-

-

is equal toc

f cr (N/mm2), when no internal pressure existsrr
E≤ 0.069 ⋅   → ≤ 60.43
t
F
t
 

In the case ofF= 1.567 × 10 2
1.5

From (3.25),

f cr 0 =In the case of

r
E r
E
0.069 ⋅   ≤ ≤ 0.807 ⋅   → 60.43 ≤ ≤ 706.7
F
t
F
t
 
 

From (3.26), r   F  0.807 −  t  ⋅  E      = 1.655 × 10 2 − 0.145 r 
f cr 0 = 0.267 ⋅ F + 0.4 ⋅ F ⋅ 
 

0.738
t




-

-

-

r
E r
In the case of 0.807 ⋅   ≤ → 706.7 ≤t
F t

0.8 ⋅ E

From (3.27),

f cr 0 =In the case ofr
E≤ 0.377 ⋅  tFFrom (3.29),

c f cr =In the case of

From (3.30),(

2.25 3 ⋅ 1 − ν

2

0.72

→

)

r
4 r
  = 4.429 × 10  
t
t

r
≤ 49.53tF= 1.567 × 10 2
1.5

E
0.377 ⋅  
F

0.72

≤

r
E
≤ 2.567 ⋅  tF→ 49.53 ≤0.72r F 
 2.567 −   ⋅   t E c f cr = 0.267 ⋅ F + 0.4 ⋅ F ⋅ 
2.190







= 1.729 × 10 2 − 0.326 (r t )Appendix A2

0.72

141

r
≤ 337.3t-In the case ofE
2.567 ⋅  
F

0.72

≤

r
r
→ 337.3 ≤ttFrom (3.31),cf cr =

 1
0.6 ⋅ E r   r ⋅ 1 − 0.901 ⋅ 1 − exp − ⋅ 16 t   t 
2.25 









 1
r 
= 5.488 × 10 4 1 − 0.901 ⋅ 1 − exp − ⋅ 16 t 

r
 
t

Allowable bending stress, b fcr , (N/mm2)

(2)
-

-

-

≤ 0.3 → 0 ≤ σ h ≤ 0.705 × 10 2

0≤

From (3.33),

b f cr = b f cr +In the case of
From (3.34),-

σh

In the case of

σh
F

F

0.7 ⋅ f cr 0 − b f cr  σ h 
⋅

0 .3
 F 

> 0.3 → 0.705 × 10 2 < σ h

b f cr = f cr 0 ⋅ 1 −

σh 



is equal toF 

f cr when no internal pressure exists (N/mm2)

b

f cr

b

f cr is given by the following equations:b

In the case ofr
E≤ 0.274 ⋅  tFFrom (3.35),

b f cr =E
In the case of 0.274 ⋅  F

From (3.36),

0.78

≤

0.78

→

r
≤ 54.05tF= 1.567 × 10 2
1.5
r
E
≤ 2.106 ⋅  
t
F

0.78

r
≤ 415.4t

0.78
 r  F   2.106 −     t  E  b f cr = 0.267 ⋅ F + 0.4 ⋅ F ⋅ 
1
.832







r
= 1.708 × 10 2 − 0.260 t

Appendix A2→ 54.05 ≤

142

E
e) In the case of 2.106 ⋅  F

From (3.37),

b

0.78

r
r
→ 415.4 ≤tt≤

f cr =

 1
0.6 ⋅ E r 
⋅ 1 − 0.731 ⋅ 1 − exp − ⋅ 16 t 2.25 









 1
r 
= 5.488 × 10 4 ⋅ 1 − 0.731 ⋅ 1 − exp − ⋅ 16 t 

r
 
t

r
 
t

Allowable shear stress, s fcr , (N/mm2)F− s f crσ 
1 .5 3f=f+⋅ h 
From (3.39),
s cr
s cr
0.3
 F 

(3)

= s f cr (Q σ h = 0 )s--f cr is given by the following equations:

In the case ofE
0.204 ⋅  
r
F
≤
0 .4
t
l
 
r

From (3.41),s f cr =In the case ofE
0.204 ⋅  
F
0.4
l
 
r

F
1.5 ⋅ 3
0.81

0.81

→

r
≤ 49.31t= 0.905 ⋅ 10 2E
1.446 ⋅  
r
F
≤ ≤
0 .4
t
l
 
r

0.81

→ 49.31 ≤From (3.42),
0 .4
0.81

r l  F 
1.446 −   ⋅   ⋅   
0.267 ⋅ F 0.4 ⋅ F 
t  r E +⋅
s f cr =
1.242
3
3 






r
= 0.994 × 10 2 − 0.182  tAppendix A2

143r≤ 349.54t-E
1.446 ⋅  
F
0 .4
l
 
r

In the case of

From (3.43),s

f cr =

0.81

≤

r
r
→ 349.54 ≤tt

4.83 ⋅ E ⋅ 0.8
lr 
2.25 ⋅  ⋅
r
t 


2

l
r 
1 + 0.0239 ⋅  ⋅
r
t 


⋅

3

r 2
= 3.534 × 10 5 ⋅ 1 + 0.0239 ⋅  
t

r
 
t

3

r
 
t

2

where: l is the section length of the buckling (mm).

A2.3.1.2

The values of allowable short-term stresses in normal times are assumed to be
1.5 times those of allowable long-term stresses.

A2.3.1.3

Allowable stress during earthquakes is given by the following equations:(1)

Allowable compressive stress, c fcr, (N/mm2)

0≤

σh
F

≤ 0.3 → 0 ≤ σ h < 0.705 × 10 2

From (3.45),

σh
F

c

f cr = c f cr +0.7 ⋅ f crs − c f cr  σ h 


0.3
 F 

> 0.3 → 0.705 × 10 2 < σ hFrom (3.46),

Symbols:cc f crf crs･ In the case of σ 
f cr = f crs ⋅ 1 − h 
F 


:
:

f cr when no internal pressure exists (N/mm2)
basic compressive buckling stress under earthquakes(N/mm2)crr
E≤ 0.069 ⋅   → ≤ 60.43tt
FFrom (3.47), f crs = F = 2.35 × 10 2
r
E r
E
･ In the case of 0.069 ⋅   ≤ ≤ 0.807 ⋅   → 60.43 ≤ ≤ 706.7
t
F t
F
From (3.48),

Appendix A2144 r   F  0.807 −  t  ⋅  E      = 2.438 × 10 2 − 0.145  r 
f crs = 0.6 ⋅ F + 0.4 ⋅ F ⋅ 
 

0.738
t




r
E r
･ In the case of 0.807 ⋅   ≤ → 706.7 ≤
t
F t

c

0 .8 ⋅ E

f crs =

From (3.49),(

3 ⋅ 1 −ν

2

r
r
  = 9.965 /  
t
t

)

f cr is given by the following equations:･ In the case ofr
E≤ 0.377 ⋅  
t
F

0.72

→

r
≤ 49.53tc f cr = F = 2.35 × 10

From (3.50),

E
･ In the case of 0377 ⋅  
F

0.72

≤

2

r
E
≤ 2.567 ⋅  
t
F

0.72

→ 49.53 ≤r≤ 337.2t

0.72
 2.567 −  r  ⋅  F  t E c f cr = 0.6 ⋅ F + 0.4 ⋅ F ⋅ 

2.190







From (3.51),

r
= 2.511 × 10 2 − 0.327  tE･ In the case of 2.567 ⋅  
F

0.72

≤

r
r
→ 337.2 ≤
t
t



 1
r   r 



0
6
1
0
901
1
f
=
.
⋅
E
⋅
−
.
⋅
−
exp
−
⋅

c cr
 16 t   t 



 


 1
r   r 

= 1.235 × 10 5 ⋅ 1 − 0.901 ⋅ 1 − exp − ⋅  
 16 t   t 
Allowable bending stress, b f cr , (N/mm2)

From (3.52),

(2)

0≤

σh
F

≤ 0.3 → 0 ≤ σ h ≤ 0.705 × 10 2

From (3.53),

Appendix A2bf cr = b f cr +0.7 ⋅ f crs − b f cr  σ h  F 
0.3



145

･ In the case ofσhF> 0.3 → 0.705 × 10 2 < σ hFrom (3.54),b

Symbols:b･ In the case of

From (3.55),

 σ 
f cr = f crs ⋅ 1 − h F f cr :

(3)

→

0.78

≤

r
≤ 54.05t2r
E≤ 2.106 ⋅  
t
F

0.78

→ 54.05 ≤r≤ 415.4t

0.78
 r  F   2.106 −     t  E  
f
=
0
.
6
⋅
F
+
0
.
4
⋅
F
⋅


b cr
1
.832






r
 
= 2.49 × 10 2 − 0.260  
t

0.78

≤

r
r
→ 415.4 ≤
t
t



 1
r   r 



f
=
0
.
6
⋅
E
⋅
1
−
0
.
731
⋅
1
−
exp
−
⋅

b cr
 16 t   t 



 


 1
r 
= 1.235 × 10 5 ⋅ 1 − 0.731 ⋅ 1 − exp − ⋅ 16 t 

Allowable shear stress, s f cr , (N/mm2)F− s f crσ 
3From (3.58), s f cr = s f cr +⋅  h  = s f cr
0.3
 F 

･ In the case ofAppendix A2

0.78

b f cr = F = 2.35 × 10

E･ In the case of 2.160 ⋅  F

From (3.57),

f cr when no internal pressure exists (N/mm2)r
E
≤ 0.274 ⋅  
t
F

E
･ In the case of 0.274 ⋅  F

From (3.56),b

E
0.204 ⋅  
r
F
≤
0 .4
t
l
 
r

146

0.81

→

r
≤ 49.31
t

r
 
t

(Qσ h = 0)

From (3.60),s f cr =F
3

E
0.204 ⋅  F･ In the case of
0.4
l
 
r

From (3.61),= 1.357 × 10 2
0.81

E
1.446 ⋅  
r
F
≤ ≤
0 .4
t
l
 
r

0.81

r
≤ 349.54t→ 49.31 ≤

0 .4
0.81


1.446 −  r  ⋅  F  ⋅  l  
0.6 ⋅ F 0.4 ⋅ F t E
r 
+
⋅

s f cr =
1
.
242
3
3 






r
= 1.446 × 10 2 − 0.181  t･ In the case of

From (3.62),

E
1.446 ⋅  
F
0 .4
l
 
r

s

f cr =

0.81

≤

r
r
→ 349.54 ≤tt

0.8 ⋅ 4.83 ⋅ El ⋅ r
r
t 


2

l
r 
⋅ 1 + 0.0239 ⋅  ⋅
r
t 


r
= 7.952 × 105 1 + 0.0239 ⋅  tAppendix A2

147

3/ 2

3

r
 
t

r
 
t

2

A2.3.1.4(1)Evaluation Results

For Long-term Loading
Table A2.14
Working Stress

Test Results

Calculation
Point

Allowable Stressr/t
σc/c fcr +σb/ b fcr
σh
σc
σb
τ
c fcr
b fcr
s fcr
(mm) (－) (N/mm2) (N/mm2) (N/mm2) (N/mm2) (N/mm2) (N/mm2) (N/mm2)≦1⑤

6.0 716.7

7.3

22.9

32.4

14.8

4.16

0.0

0.0

④

8.0 537.5

12.1

36.8

47.8

21.2

9.27

0.0

0.0

③

10.0 430.0

13.4

49.5

62.2

28.0

15.20

0.0

0.0

②

11.0 390.9

14.6

56.0

70.6

31.5

22.67

0.0

0.0

①

12.0 358.3

14.6

62.0

78.0

35.2

27.83

0.0

0.0

⓪

15.2 282.9

0.0

80.6

97.3

48.3

31.23

0.0

0.0

t

(2)

0.18+0.0=0.18<1
0.25+0.0=0.25<1
0.31+0.0=0.31<1
0.40+0.0=0.40<1
0.45+0.0=0.45<1
0.39+0.0=0.39
<1

τ/
s fcr
≦1
0<1
0<1
0<1
0<1
0<1
0<1

For Seismic Loading (Modified Seismic Coefficient Method)
Table A2.15
Working Stress

Test Results

Calculation
Point

Allowable Stress(mm) (―) (N/mm ) (N/mm ) (N/mm ) (N/mm ) (N/mm ) (N/mm ) (N/mm )⑤

6.0 716.7

7.3

51.5

72.9

33.2

3.66

7.94

5.7

④

8.0 537.5

12.1

79.0

103.8

47.6

7.02

28.70

16.3

③ 10.0 430.0

13.4

104.6

133.3

63.0

10.97

51.98

20.8

② 11.0 390.9

14.6

117.6

144.1

71.0

16.01

80.71

24.0

① 12.0 358.3

14.6

130.4

151.6

79.2

19.63

98.55

24.6

⓪ 15.2 282.9

0.0

158.9

175.7

93.6

25.51

107.21

20.4

t

r/t

σh

c fcr
2

Appendix A2

2

b fcr

2

s fcr

2

σc

σb
2

148

σc/c fcr +σb/ b fcr

τ
2

2

≦1
0.07+0.11=
0.18<10.09+0.28=
0.37<1
0.10+0.39=
0.49<1
0.14+0.56=
0.70<1
0.15+0.65=
0.80<1
0.16+0.61=
0.77<1τ/s fcr≦1
0.17<10.34<10.33<10.34<10.31<10.22<1(3)For Seismic Loading (Modal Analysis)Table A2.16
Working Stress

Test Results

Calculation
Point

Allowable Stress

(mm)

⑤

6.0 716.7

7.3

51.5

72.9

33.2

3.66

8.11

5.8

④

8.0 537.5

12.1

79.0

103.8

47.6

7.02

25.88

14.2

③

10.0 430.0

13.4

104.6

133.3

63.0

10.97

43.80

16.6

②

11.0 390.9

14.6

117.6

144.1

71.0

16.01

65.33

18.3

①

12.0 358.3

14.6

130.4

151.6

79.2

19.63

78.21

18.3

⓪

15.2 282.9

0.0

158.9

175.7

93.6

25.51

83.38

15.0

(4)

t

r/t
(―)

σh

c fcr
2

2

b fcr

2

s fcr

2

σc

σb
2

τ/
s fcr
≦1

σc/c fcr +σb/ b fcr

τ2

2(N/mm ) (N/mm ) (N/mm ) (N/mm ) (N/mm ) (N/mm ) (N/mm )≦10.07+0.11=
0.18<10.09+0.25=
0.34<1
0.10+0.33=
0.43<1
0.14+0.45=
0.59<1
0.15+0.52=
0.67<1
0.16+0.47=
0.63<10.17<10.30<10.26<10.26<10.23<10.16<1For Short-term Loading (under Wind Pressure)
Table A2.17
Working Stress

Test Results

Calculation
Point

Allowable Stress

(mm)

⑤

6.0 716.7

7.3

34.4

48.6

22.2

4.16

0.489

0.7

④

8.0 537.5

12.1

55.2

71.7

31.8

9.27

1.445

1.0

③

10.0 430.0

13.4

74.3

93.3

42.0

15.20

2.551

1.2

②

11.0 390.9

14.6

84.0

105.9

47.3

22.67

4.036

1.4

①

12.0 358.3

14.6

93.0

117.0

52.8

27.83

5.056

1.4

⓪

15.2 282.9

0.0

120.9

146.0

72.5

31.23

5.755

1.3

t

r/t
(―)

Appendix A2

σh

c fcr
2

2

b fcr

2

s fcr

2

σc

σb
2

σc/c fcr +σb/ b fcr

τ2

2(N/mm ) (N/mm ) (N/mm ) (N/mm ) (N/mm ) (N/mm ) (N/mm )149≦10.12+0.01=
0.13<1
0.17+0.02=
0.19<1
0.20+0.03=
0.23<10.27+0.04=
0.31<10.30+0.04=
0.34<1
0.26+0.04=
0.30<1τ/s fcr≦1
0.03<10.03<10.03<10.03<10.03<10.02<1

A2.3.2 Evaluation of Stresses in Hopper Wall
The axial stress σΦ in the longitudinal direction and the circumferential tensilestressσθ in the hopper wall are evaluated by the following equations:σφ =From (5.15),

dPv ⋅ d'≤ ft
4 ⋅ t2 ⋅ sinαlc ⋅ t2 ⋅ sinαdPa ⋅ d '≤ ftσθ =2 ⋅ t 2 ⋅ sin α

From (5.16)
Symbols:(Wh + Ws ) +

Wh
Ws
t2
lc
d’
dPν

α
dPa

:
:
:
:
:
:
:
:

weight of stored material in the hopper (N)weight of the hopper beneath the section (N)thickness of the hopper wall (mm)perimeter of the hopper wall (mm)diameter of the hopper wall (mm)design vertical pressure (N/mm2)
inclination of the hopper wall (°)
perpendicular design pressure of the hopper wall (N/mm2)

Table A2.18 Stresses in Hopper Wall2

(N/mm )

(mm)

σθ
σφ
ft
2
2
(―) (N/mm ) (N/mm ) (N/mm ) (N/mm2)

8.0 8,600 612.00 50.60 27.02 0.150

5.657

2.30

0.163

157

61.3

123.9

6.0 5,733 181.50 16.90 18.01 0.153

4.243

1.87

0.135

157

54.3

91.2

6.0 2,867 22.60

4.243

1.44

0.111

157

27.4

37.5

t2

d'

(mm) (mm)Appendix A2

Wh
(kN)

Ws
(kN)

4.20

lc
(m)

9.01

dPν

t 2 sin α

0.158

150

Cd’

dPa

2

APPEDIX A3DESIG EXAMPLESUPPORTIG STRUCTURES FOR SPHERICAL TAKS

inside diameter
manholetop landing

staircasebracing

12～supportingcolumnsstaircase

Fig.A3.1

General View of Spherical Tank

Specifications:-Classification of ground: Type 3-

Importance factor, I = 1.5 (Classification of seismic design IV)-

Structural type:pipe columns and pipe bracings-

Stored contents:liquefied butane (specific gravity 0.587)-

Sphere:D = 1 790 cm, t = 12.0~15.0 mm, material SPV 490 (F = 42.7 kN/ cm2)-

Columns:dc =φ609.6 , tc = 14, n = 12, material SPV 490-

Bracings:

db =φ267.4, tb = 7, L = 1 012, material STK 400 (F = 23.5 kN/ cm2)-ho

= 1 136.4 cm-

h1= 900 cm-l

= 463 cm

Appendix A3

151

-

θ = 1.096 rad

The dead weight of the structure WD and weight of the stored contents WI are calculated asfollows:
WD = 1 980 kN
WI = 15 550 kN
As f(ξ) is obtained as 0.7 from Figure 6.1, the design weight imposed on the base of structure
W is calculated at;
W = WD + f(ξ) WI = 12 865 kN
Those basic values to assess the retaining strength and elastic modulus of the structure are
obtained as follows:
lr =

Dd c
2

h = h0 –= 233.5 cm

lr= 1019.7cm
2h2 = h – h1 = 119.7 cm
Axial force in the column:=

WD + WIn= 1 461 kN

Sectional properties of the column:
Ac = 262 cm2
Zc = 3 810 cm3
Ic = 1.160 × 105 cm4
ic = 21.1 cm

λc = hic = 48.3 < 60

Sectional properties of the bracing:
Ac = 57.3 cm2
ib = 9.21 cm
Short-term allowable bending moment of the column:
1461 
My = FZc 1 −   = 42.7 × 3 810 × 1 −
 = 1 414 kN･m
 FA 
 42.7 × 262 c 

Short-term allowable shear force of the column:

Appendix A3152

Qca =

F
3

×Ac= 3 230 kN
2

Short-term allowable compressive force of the bracing:

λb = (L ib ) ×

1= 54.92

1.5fc = 19.7 kN/cm2, ∴ c = 1.5fcAb = 1 129 kN
Sy = c cosθ= 516 kN
Stresses in the column in Figure 6.3.2 of the Commentary are calculated as follows:
From Eq. (6.3.8),
Ma1 = 198.9S
Mb1 = – 35.7S
Ra = (M a1 − M b1 ) / h2 = 1.960 S

From Eq. (6.3.9),
Ma2 = 31.4S
Mb2 = 27.7S

R′ = M a 2 / h2 = 0.031S
Q = Ra + R′ = 1.99S
From Eq. (6.3.12) and (6.3.13),
Ma =230.3S
Mb = －8.0S
From Eq. (6.3.7),

δ = 0.00458S
The yield shear force Qy in the plain where the maximum shear stress takes place is given
as follows:
Qy =1.99Sy = 1,027 kN
From Eq. (6.3.1), the yield shear force Fy of the supporting structure is given as follows:

Appendix A3

153

Fy =

nQ y
2

= 6 161kN

The yield deformation δ y is given as follows:δ y = 0.00458S = 2.363 cm
The spring constant k is obtained from Eq. (6.3.18).
k=

Fy

δy

= 2 607kN/cm

The natural period: T = 2π

W
= 0.446 s
gk

The yield shear force coefficient: α =FyW= 0.479

The bending moment eMa in Figure 6.3.3 of the Commentary is calculated as per Eq.
(6.3.19):
eMa = Ma −lr Q
= 114.1S4∴ ( e M a )y = 114.1S y = 589 kN ⋅ m
From Eq. (6.3.21),∆ M = My – (eMa)y = 1414 – 589 = 825 kN ⋅ m
From Eq. (6.3.22),∆ Ma = ∆ M ×h= 875 kN ⋅ mh − lr 4

From Eq. (6.3.23), (6.3.24),

η=∆M a h 2/ δ y = 5.397
3EI∴ aη = 0.75 × 5.397 = 4.048 > 3.0 → 3.0
From Eq. (6.3),

Appendix A3

154

Ds =

1
1 + 4 aη= 0.278

From Eq. (6.2),
C = Z s ⋅ I ⋅ Ds ⋅S a1
= 1.0 × 1.5 × 0.278 × 1.0 = 0.417
g

Evaluation of anti-seismic characteristics is made in accordance with Eq. (6.5).
Fy = 6,161 kN
Qd = CW = 0.417×12,865 = 5,365 kN
Fy ≥ Qd

Appendix A3

155

APPEDIX A4DESIG EXAMPLABOVE-GROUD, VERTICAL, CYLIDRICAL STORAGE TAKS
A seismic capacity evaluation is conducted for a small tank which suffered elephant foot bulge
type damage during the Hanshin-Awaji Great Earthquake in Japan. The evaluation is
performed according to the flow scheme shown in Fig. 7.3.1.
A4.1

Tank size
Diameter of the cylindrical tank
Weight of the stored contents
Depth of the stored liquid
Static pressure imposed on the bottom plate
Slope angle of the conical roof
Wall thickness of the cylindrical tank- bottom course- one-third height above the bottom- thickness of the bottom annular plate
AnchorA4.2

t0 = 8 mm
t1/3 = 6 mmt
= 6 mmunanchoredDesign specification:
Yield point of bottom plate and wall plate
Young’s modulus of bottom plate and wall plate:
Poisson’s ratio
Ratio of yield stress to ultimate tensile strength
Seismic zone factor
Importance factor
Damping ratio

Critical period- (determined by the ground classification)
A4.3

D = 13.54 m
Wl = 1.905 x 104kN
Hl = 13.5 m
p = 0.0132 kN/cm2
15 o (estimated)

σ y = 23.5 kN/cm2
E = 2.06 x 104 kN/cm2
v = 0.3
Yr below 0.8 (estimated)
Zs = 1.0
I= 1.2
h = 0.1(for impulsive mass)= 0.001(for convective mass)
TG = 0.96 sec.

Evaluation for the mass of impulsive pressure
Natural period Tf is obtained using Eq. (7.2.20):
Tf

=
=
=

2

λ

Wl
πgEt1 / 3

2
1.905 × 10 40.23 π × 980 × 2.06 × 10 4 × 0.6
0.196 sec

Weight percent of the effective mass of the tank contents is given by Fig.7.2.1.

Appendix A4156

ff (D/Hl) = ff (1.00) = 0.77
1) Evaluation for buckling
Substituting values into Eq. (7.6):
Ds

=

0.71
1 + 3h + 1.2 h

1

⋅

2

=

 Tf 
0.5 1 + 3  Te 
0.71
1⋅
2
1 + 3 × 0.1 + 1.2 0.1
 0.196 
0.5 1 + 3
 0.397 
0.423 × 0.76 / 0.5

=

0.642

=

Design shear force of the impulsive mass vibration, Qdw, is given by Eq. (7.3.1) as:
Qdw

=
=
=

Zs . I . Ds . Sa1 . ff . Wl / g
1.0×1.2×0.642×980×0.77×1.905×104 /980
1.130×104 kN.

Mean hoop tensile stress in the cylindrical wall, Q σ hd , is calculated by Eq. (7.3.2)
as:

Q σ hd

=

Qdw
W
+ l
2.5H l t0 πrt 0=1.13 × 10 41.905 × 10 4+
2.5 × 1350 × 0.8 π × 677 × 0.8
4.185 + 11.20=

15.38 kN/cm 2=

Allowable bending stress in the cylindrical wall, b fcr, and basic axial compressive
buckling stress fcrs are calculated by Eqs. (3.54) and (3.28) as:
b fcr

f crs

=

fcrs (1－σhd / Qy)=

11.79(1－15.38/23.5)=

4.072 kN/cm2

=

0.8E

=
=

Appendix A4r
3 1 − ν 2   t0 
0.8 × 2.06 × 10 4()

3(1 − 0.09 ) × (677 / 0.8)
11.79 kN/cm 2

157

Yield shear force of the storage tank for buckling as shown in Eq. (7.3.6) is:
e Qy

=
=
=

π ⋅ r 2 ⋅ b f cr ⋅ t 0
0.44 H lπ (677 )2 × 4.072 × 0.8
0.44 × 1350
7.896 × 10 3 kN

It can be concluded that the buckling strength is insufficient as eQy < Qdw.
2)

Evaluation for bottom plate
The yielding force of the bottom plate qy, the uplifting deformation for the yielding
force δy, the spring constant of unit circumferential length for uplifting resistance k1,
the spring constant K1 as the one degree of freedom system that represents thehorizontal directional motion of tank and the natural period of the storage tankwhen only the bottom plate deforms, T1, are calculated from Eqs. (7.2.11), (7.2.15),(7.2.17) and (7.2.18), respectively, as:

qy

=

k1

K1

T1

=

=

( 3tσy2) / (8Ep)=(3×0.6×23.52) / (8×2.06×104×0.0132)=

0.457 cm=

qy /σy=0.273 / 0.457=

0.597 kN / cm2=

48.7r3 k1 / Hl2=

48.7×(677)3×0.597 / 13502=

4950 kN / cm

=

2π

=
=

Appendix A43
2 × 0.6 × 1.5 × 0.0132 × 23.5
3
0.273 kN/cm=

δy2t 1.5 Pσ y

f f WlgK 1

0.77 × 1.905 × 10 4980 × 4950
0.345 sec
2π

158

The modified natural period of the storage tank considering deformation of both thewall plate and the bottom plate, Te, the coefficient determined by the ductility of thebottom plate uplifting, Dη, and the structural characteristic coefficients, Ds arecalculated from Eqs. (7.2.19), (7.2.24) and (7.3), respectively, as
Te

Dη

=

T f2 + T12=

0.196 2 + 0.345 2=

0.397 sec

=

1

=
=
Ds

=
=
=

1 + 84(T1 / Te )2

1
1 + 84(0.345 / 0.397 )

2

0.125
1.42
1 + 3h + 1.2 h
0.846 × 0.125
0.105

⋅ Dη

Design shear force of the impulsive mass vibration, Qdw, and the design yield shear
force of the storage tank, Qy, are obtained from Eqs. (7.3.1) and (7.3.5) as:
Qdw

Qy

=

Z s ⋅ I ⋅ Ds ⋅ S a1 ⋅ f f ⋅ Wl / g=

1.0×1.2×0.105×1.0×0.77×1.905×104=

1848 kN=(2π r2qy) / (0.44Hl)=(2π×(677)2×0.273) / (0.44×1350)=

1324 kN.

Therefore it can be concluded that the uplifting strength is insufficient because eQy< Qdw.
A4.4

Evaluation for the mass of convective pressure
As ff is obtained from Fig.7.2.1, fs, the natural period of sloshing response Ts, andacceleration response spectrum of the first natural period Sa1 are calculated by
Eqs.(7.2.5) and (7.2.32) as

fs

Appendix A4=

1 － ff=

1 － 0.77

159

=
Ts

IS v

=

2π D
3.682 g tanh (3.682 H l / D )=2π 13543.682 × 980 tanh (3.682 × 1350 / 1354 )=

3.85 sec=

1.10
1 + 3h + 1.2 h
2.11m/s
2π
IS v ×Ts=∴ IS a1

0.23

=
=

2⋅

3.45m/s 2

Design shear force of the convective mass vibration Qds is obtained from Eq. (7.3.3)
as:
Qds

=

Z s ⋅ I ⋅ S a1 ⋅ f s ⋅ Wl / g=

1.0×345×0.23×1.905×104 / 980=

1542 kN.

The design yield shear force for buckling, sQy is obtained by Eq. (7.3.4) assQy=

0.44 eQy=

0.44×7.896×103=

3474 kN

and the design yield shear force for uplifting, sQy is obtained from Eq. (7.3.4) as:sQy=

0.44Qy=0.44×1324=

583 kN.

From these results, it can be concluded that buckling strength is sufficient as sQy >Qds, and uplifting strength is insufficient because sQy < Qds.
In conclusion, this tank has insufficient capacity for impulsive mass vibration
concerning both buckling and uplifting strength, and for convective mass vibration
concerning uplifting strength. This therefore, explains the actual seismic damage,
namely buckling type damage, to the subject tank.
A4.5

Sloshing force acting on the fixed roof

Appendix A4160

When only the 1st mode is considered, sloshing height ζ R is equal to η s .
Therefore,

ζR=

0.802 ⋅ Z s I ⋅ S v1

D
 3.682 H i tanh 

g
D 


=

0.802 ⋅ 1.0 ⋅ 2.11 ⋅13.54
 3.682 ⋅ 13.5 tanh 

9.8
 13.54 =
1.99m
Since the roof angle ϕ = 0.262 rad., the impulsive pressure Pi is, from equation(7.2.37),ππ
 15π 2
ρ cot ϕ (ζ&rA cos θ ) 2 = ⋅ 103 ⋅ cot Pi= ⋅ (1.99 ⋅ 1.63)
2
2
 180 =6.17 × 10 4 N/m 2
The hydrodynamic pressure Ph at the connection of the roof and wall is, fromequation (7.2.38),
Ph=ρ g (ζ r − hr ) cos θ = 1000 ⋅ 9.8 ⋅ (1.99 − 0)=Appendix A4

1.95 × 10 4 N/m 2

161

APPEDIX B
ASSESSMET OF SEISMIC DESIGS
FORUDER-GROUD STORAGE TAKS
B.1

Introduction
This chapter addresses that the seismic design for under-ground tanks has been revisedin the response displacement method. Some sample calculations are provided toevaluate the response displacement method comparing with other analytical methods.

B.2

Analytical Methods
The recommended analytical methods are; the response intensity method and dynamic
analysis using an FEM model other than the response displacement method. The mainchange in the response displacement method from the 1990 version [B.1] is to consider
an inertia force and surrounding shear forces in the new version. The five analyticalmethods shown in Figure B.1 are assessed in the following sections whereas Figure B.2
compares the new preferred analytical methods with the one in the previous (1990)
version of the Design Recommendation.
Dynamic Analysis
Reference [B.2]

SHAKE
Precise Response Displacement Method
Response Displacement Method

Recommendation
Simplified Response Displacement Method

SHAKE
Response Intensity Method

Precise Response Intensity Method

Recommendation
Simplified Response Intensity Method

Figure B.1
T he Previous Version

Recommending Method in the New Version

Response Displacement Me th od

Response Displacement Method

Resp onse In tensity M etho d

Dy na m ic Ana lys is

Modeling of Structure

Modeling of Ta nk and Co nte nts

M odeling of Tank and Con ten ts

M odelin g o f Tan k and Co ntents

Modeling of Soil (Ground Spring)M o d e l i n g o f S o i l ( FE M et c . )M o d e li n g o f So i l ( F E M et c . )G r o u n d I n t e n s i ty

I np u t Seis mic Wa ve

Calculation of Seismic Str ess

Calculation of Seismic Stress

Se ismic Soil Pre ssur e
Calcu lation o f
Gro und Spr in g

Gro undDe for ma tio nSurr oun ding
Shear
Force

In er tia
For ceCalculation ofGrou nd Defo rma tion

Calculation of Seismic StressCalculation of Seismic Stress

Figure B.2

Appendix B

162

Par.3.6.2.3 generally recommends the use of a reverse triangle distribution for
calculating ground displacement by the simplified response displacement method.
However, it recommends that a cosine-curve distribution is used where the bedrock islocated at the base of the under-ground tank. This appendix reports the assessment resultof the ground displacement with an assumption of a cosine-curve distribution used in the1990 version of the Design Recommendation and compares to the one with anassumption of a reverse triangle distribution (see Figure B.2).
B.3

Assessment of Structure with GroundFigure B.3 illustrates the model of an under-ground water tank with a flat slab structure,
which is to be analysed. Table B.1 shows the physical properties of the ground which,for the purposes of this example, is assumed as sandy ground.

”
Column

4,000

”

”

”

”

”

”

”

4,000

5,100

6,000
Underground
Water Level

4,000

Outer Wall

500

500

500

28,000

Water Level

”

”

”

400

”

2,000

4,000

GL

4,000

4,000

4,000

36,000Figure B.3

Table B.1
Soil

Reclaimed

Thickness

N

Vs

ρ

G0

(m)

Value

(m/s)

(kN/m3)(kN/㎡)

2.0

5

100

15.68

15,680ν

E0P(kN/㎡)(N/c ㎡)

0.33

42,140

0.98
3.92

Clay
Grain

3.0

10

170

16.66

49,0000.33

127,400

Sand

3.0

15

200

16.66

67,6200.33

176,400

6.86

Sand ～

12.0

30

250

16.66

107,800

0.49

323,400

11.76

20.0

50

300

18.62

166,600

0.49

499,800

20.58

30.0

>50

400

18.62

303,800

0.49

901,600

35.28

>50

500

18.62

470,400

0.49

－

－

GravelSand
Bed Rock

Grand Water Level ： GL－8m： S-Wave Velocity = 80.6 N0.331
Vsρ： Unit Volume WeightG0： Shear Modulus of Elastic Groundν： Poisson’s Ratio
E0： Young's Modulus of Elastic Ground： Restraint Pressure ＝ (1+2K0)/3×σv
P： Static Soil Pressure Coefficient = 0.5
K0σv： Imposed Load

Appendix B

163

B.4

Outline of Analytical Method
Figure B.4 shows a model of dynamic analysis and its assessment flow. The structurewith the surrounding ground is modeled with a two dimensional FEM model on theassumption that an energy transmission boundary is provided at the side of the model,
and there is a viscous boundary at the bottom of the model. The frequency response
analysis is adopted by use of the EL CENTRO seismic wave as an input, and theconvergence stiffness is adopted as a stiffness of the ground as the result of SHAKE[B.2].
Figure B.5 illustrates the structure’s model so that the outer wall in the model plane ismodeled as plane elements, and the outer wall perpendicular to the plane, columns,
bottom slab and top slab as beam elements.
GroundStructure
Assessment Flow
Input Wave

2E wave at the bedrock.

Energy Transmission Boundary

70ｍ

It defines El Centro Earthquake asItadjustsforthemaximum

acceleration to become 25cm/s at thesurface.
Stiffness of Ground
It sets based on the result of SHAKE.
It uses the convergence stiffness by
the result of SHAKE for the stiffnessof the ground.

36ｍ

18ｍ
Viscose Boundary

Figure B.4
Outer Wall (perpendicular to the plane)：Beam Element
Top and Bottom Slabs：Beam Elements
Columns：Beam Elements

Outer Wall (in the model plane)：Plane Element

Figure B.5

Appendix B

164

Figure B.6 shows a design flow of the response intensity method and Figure B.7 shows
its analytical model. The structure, together with the surrounding ground, is modeledwith a two dimensional FEM model on the assumption that a ground seismic coefficientis provided in the ground, and there is a static inertia force in the structure. An
assessment of the ground stiffness and the seismic ground coefficient, (the inertia force),
is made by means of both the simplified calculation based on the 1990 version of the
Design Recommendation and the precise calculation based on the result of SHAKE.
MethodAnalysis Flow
Simplified Calculation

Ground Stiffness

Modeling of Structure

Modeling of GroundCalculating Inertia Force
Horizontal Coefficient(Seismic Ground Coefficient)

Response Intensity Method

The initial stiffness is calculated considering thenon-linear characteristics of soil.
・to calculate the ground deformation inaccordance with the previous recommendation.
・to calculate the ground stiffeness considering thedegradation of rigidity due to the shear strain.

Seismic Ground Coefficient

Calculating Stiffnessof GroundPrecise Calculation

It is provided by the result of SHAKE.
・ The ground stiffeness is considered as theconvergence stiffness based on the result of
one-dimensional equivalent linear analysis.

It is provided by the result of SHAKE.
The Design Recommendation calls for;Khs = ZS・RG・K0

・It is calculated from the maximum shear force
of one-dimensional equivalent linear analysis.

Kh =2(|τimax|－|τi-1max|)
hi -1・ρi-1 + hi・ρi

Figure B.6
Boundary Condition

Structure(Horizontal direction is free.)Ground70ｍ

Inertia Force

Seismic Ground Coefficient

90ｍ

Boundary Condition (Fix)

36ｍ

90ｍ

Figure B.7
Figure B.8 shows the design flow of the response displacement method, and Figure B.9
shows its analytical model. The structure, together with the surrounding ground, ismodeled with a two dimensional FEM model on the assumption that the ground springs,(the axial and shear springs), are provided to the structure and the under-grounddeformation, ground seismic coefficient (inertia force) and the boundary shear forces areprovided statically. An assessment of the springs, seismic coefficient and shear forces ismade by means of both the simplified calculation based on the Specifications for
Highway Bridges [B.3], the 1990 version of the Design Recommendation, and theprecise calculation based on the result of SHAKE (Refer to Figure B.10).

Appendix B

165

Analytical Method

Analysis Flow

Precise Calculation

Simplified Calculation
Specifications for Highway Bridges
Ground Stiffness

Modeling of Structure

Calculating Inertia Force
Horizontal Coefficient(Seismic Ground Coefficient)Calculating Surrounding
Shear Force
Maximum Seismic
Ground Shear Force・ The ground stiffeness is considered as theconvergence stiffness based on the result of
one-dimensional equivalent linear analysis.

It is provided by the result of SHAKE.

The Design Recommendation calls for;・ Maximum ground deformation ofUH(Z)=

2

・Sv・Khb・Ts・cos（

π2

πZ
）
2H

one-dimensional equivalent linear analysis.

It is provided by the result of SHAKE.
The Design Recommendation calls for;・It is calculated from the maximum shear force
of one-dimensional equivalent linear analysis.

Khs = ZS・RG・K0 ,etc.

Kh =Surround Shear Force

Response Displacement Method

Ground Deformation

Calculating Stiffnessof Ground

Seismic Ground Coefficient

Modeling of Ground

It is provided by the result of the Static FEM
analysis.

Kv=12.8Kvo・Bv-3/4
Kh=12.8Kho・Bh-3/4
Kvo=Kho=1/30・α・Eo
Eo=28N
α=2.02(|τimax |－|τi-1max|)
hi-1・ρi-1 + hi・ρi

It is provided by the result of SHAKE.

τ= G × γ
・ It calculates the shear strain γ from theground deformation in accordance with・ Maximum ground Shear Force of
one-dimensional equivalent linear analysis.

the Design Recommendation.
・ It calculates the ground stiffness G
consideringthedegradation

of

the

stiffeness due to the shear strain.

Figure B.8
Structure
Shear Spring

6m

Axis Spring

Surrounding
Shear Force

Inertia ForceShear Spring
Axis Spring

Under-ground
Deformation

36m

Figure B.9
P

Deformation δ
Value of Ground Spring =δ

70ｍ

Unit of Force P

Boundary Condition (Fix)
90ｍ

36ｍ

Figure B.10

Appendix B

166

90ｍ

B.5

One-dimensional Seismic Motion Analysis (SHAKE)
The outline of one dimensional seismic motion analysis is shown in Figure B.11. The
stiffness and damping of the ground have strain dependent and restraint pressure
dependent characteristics [B.4]. The non-linear characteristic shown in Figure B.12 isassumed as a sandy ground.
The EL CENTRO seismic wave is adopted as the design seismic wave as being the 2E
wave at -70 m level from G.L. (Ground Level), and the maximum seismic velocity of
approximately 25 cm/s at the ground surface.
Figure B.13 shows the result of SHAKE. The maximum seismic acceleration at the
ground surface is approximately 220cm/s2, the maximum relative displacement at the
ground surface is approximately 2.6 cm to -70 m level from G.L., the maximum groundshear strain near the ground surface is approximately 4 - 6×10-4 and the ratio of theconvergence stiffness to the initial stiffness (the degradation percentage of the groundstiffness) is 0.4 - 0.8.

GL 0m

E’

F’

Non-linear Characteristic of Soil・Considering Strain Dependence of
Ground stiffness and Ground Damping
・One-dimensional Equivalent Linear Analysis

One Dimensional
Ground Model(Multiple
Reflection Theory)

Input Wave
・Design Seismic Wave (EL CENTRO)
as being 2E Wave of －70 m Level from GL・Maximum Seismic Velocity
at Ground Surface about 25 cm/s.

E

F

GL –70m

Figure B.11

Appendix B

167

1.0
Decline Ratio of
Ground Stiffness G／G0

P : Restraint Pressure(N／c ㎡)
0.8

P=49.0P=29.4P=9.8

0.6

P=19.6

P=4.9
0.4
P=0.98
0.2
0.0
10－6

10－5

10－4

10－3

10－2

10－1

Shear Strain
0.5
P : Restraint Pressure(N／c ㎡)
Damping Ratio h0.4
P=0.98
0.3
P=9.8

P=4.9

0.2
P=19.6
P=29.40.1

P=49.0
0.0 －
10 6

10－5

10－4

10－3

10－2

10－1

Shear Strain

Figure B.12
Maximum Response
Velocity
0
－20

(m)

－20

Depth

－40

－60

0

－40

－60
100

200

300

0

10

20

30

Acceleration (cm／sec2)

Velocity (cm／sec)

Maximum Relative
Displacement

Maximum
Shear StressMaximum
Shear Strain

Decline Ratio
0

－20

－20

－20

－20

(m)
Depth

－40

－60

－40

1

2

3

Displacement (cm)

0

50

100

Shear Stress (KN／cm2)

－60

0.0

168

1.0

1.5

Shear Strain (×10－3)

Figure B.13

Appendix B

－40

－40

－60

－60

0

Depth (m)

0

(m)

0

Depth (m)

0

Depth

Depth

(m)

Maximum Response
Acceleration
0

0.0

0.5

1.0

Decline Ratio

B.6

Results of Analyses

1) Ground Stiffness and Seismic Ground Coefficient
Figure B.14 shows a comparison of the ground stiffness and seismic ground coefficient
between the simplified calculation and the precise calculation. The seismic coefficient of
the precise calculation was calculated by the maximum ground shear stress of SHAKE,
and the ground stiffness was the convergence stiffness of SHAKE.
The seismic ground coefficient of the simplified calculation is based on the current
version of the Design Recommendation whilst the ground stiffness was calculated by theground deformation based on the 1990 version considering the strain dependent
characteristic of the ground stiffness (the effective strain is 65 % of the maximum
strain).
：Precise Calculation：Simplifed Calculation
Ground Stiffness E
（

N

ｍｍ2 ）

Seismic Ground Coefficient K
Reclaimed Clay

Grain Sand

2
0.05

4
0.10

6

8

10 (×105)

0.15

0.20

0.25

2.0ｍ
6.0ｍ

Structure

(8m)

Sand
12.0ｍ
Seismic Ground
Coefficient

Gravel(20m)
Seismic Ground
Coefficient
Sand

20.0ｍ

Seismic Ground Stiffness

Seismic Ground
Stiffness(40m)(14.2)(20.5)

Sand

30.0ｍ

(70m)
Bed Rock

Figure B.14

Appendix B

169

2) Ground Displacement, Surrounding Shear Force, Ground Spring.
Figure B.15 shows a comparison of the ground displacement, the surrounding shear
force and the ground springs (the axis springs) between the simplified calculation andthe precise calculation.
As the ground displacement obtained from the simplified calculation based on the 1990
version is given in the cosine-curve distribution, the magnitude of the grounddisplacement is relatively big near the bedrock, but smaller near the structure at the
ground surface layer, therefore the ground displacement obtained from the response
displacement method and the surrounding shear force calculated by this grounddisplacement are small.
The ground springs calculated by the N values obtained from the simplified calculation
following the Specifications for Highway Bridges are bigger than the ones calculated by
the convergence stiffness of SHAKE obtained from the static FEM analysis.

GroundDeformation (cm)GL－2.0ｍ
▽

0.0

0.1

0.2

0.3

SimplePrecise
Structure

GL－8.0ｍ▽

GL－2.0ｍ▽

Surrounding Shear
Force (KN)
0500

GL－2.0ｍ▽1000 1500

SimpleGround Spring(KN／mm)0500

1000 1500

Precise
SimplePrecise
Structure

Structure

GL－8.0ｍ▽

GL－8.0ｍ▽

Figure B.15

Appendix B

170

3) Comparison of the Stresses Shared in the Supporting Structure
Figure B.16 shows a comparison of the stresses shared in the frame components (i.e. theouter wall perpendicular to the model plane, columns, bottom slab and top slab) and thewall component (i.e. the outer wall in the model plane) in the supporting structure for thevarious analyses. The figure shown at the bottom of each rectangular bar indicates theratio for each analysis compared to 1.0 for the dynamic analysis.
It is noted that the stress shared in the frame is extremely small in that all analysis cases
show values of only a few percent. The total stresses obtained from the response
intensity method and the response displacement method (except for the simplified
calculation of the latter), indicate a relatively larger value which is 1.1 - 1.4 times thetotal stress of the dynamic analysis.
The value in (F.E.M Dynamic
Analysis Model

Load Shear Force of Lowest Layer

25

Wall

(kN×103)

Precise Calculation
Response Intensity

Load Force of Frame
Ground Deformation
Seismic Ground Coefficient
Surround Shear Force
Inertia Force of Structure

Response Deformation

Simple Calculation
Response Intensity

Response Deformation(2.4％)（1.9％）(2.4％)
20) shows the ratio of the load stress of the frame to the wall.

9％（1.7％）

67％
15

10

68％

64％
(－0.9％)
15％

100％

5

33％

32％

27％

1.13

1.18

0
1.00

1.43

68％

17％
0.51

Figure B.16

4) Examination about the Simplified Calculation in the Response Displacement Method.
As the total stress obtained from the simplified calculation of the response displacement
method is smaller compared to the other analysis result, an examination is carried out
with three cases in Table B.2. CASE 1 is of an examination to see the effect of theground spring, and CASE 2 and CASE 3 are examinations to see the effect of the grounddisplacement.
Figure B.17 shows the ground displacements and the surrounding shear forces of thesimplified response displacement method, the precise response displacement method andthe three cases. The analysis results of the sheared stresses at the lowest layer of thestructure are shown in Figure B.18. Although the total stress for CASE 1 indicatesapproximately 1.4 times of the simplified calculation due to the different ground spring,
it is still smaller than that of the dynamic analysis. The total stresses for CASE 2 andAppendix B

171

CASE 3 indicate 1.0 - 1.4 times of the dynamic analysis due to the increased
surrounding shear force.
Figure B.19 shows the analysis result in the case that a fall of the ground stiffness is notconsidered to calculate the surrounding shear force. All cases indicate the bigger total
stress due to the increase of the surround shear force. It is noted that the total stresses for
CASE 2 and CASE 3 indicate considerably larger when compared with the dynamic
analysis.
From the above, and with the consideration that the selected structure is a proper
example for the examination, the ground displacement of CASE2 (see Figure B.18) is
recommended as for the simplified calculation of the response displacement method.

Table B.2
CASE－1

It uses the ground spring calculated by the precise analysis.
is an opposite triangle
The ground deformation
distribution setting the same deformation at the groundsurface.
It changes the position of the bedrock into GL-40m fromGL-70m.

CASE－2

CASE－3Surrounding Shear Force (KN)

Ground Deformation to GL-8.0m (cm)▽GL－2.0ｍ 0.0

0.1

0.2

0.3

0.4

▽GL－2.0ｍ

0.5

0

250

500

1000 1250

750

CASE-2

Simple、CASE-1
CASE-3

Simple､CASE-1
Structure

CASE-3

Precise

Structure
PreciseCASE-2▽GL－8.0ｍ▽GL－8.0ｍ

Figure B.17

Appendix B

172

The value in (F.E.M Dynamic
Analysis Model(kN×103)

Simple Calculation
Response Deformation
Method) shows the ratio of the load stress of the frame to the wall.
Simple Calculation in Response Deformation Method

CASE-1

CASE-2

CASE-3（2.0％）

Load Force of Frame
Wall

Load Shear Force of Lowest Layer

Ground Deformation

25

Surround Shear Force

26％

Inertia Force of Structure（－0.1％）（1.7％）

20

34％
15（1.4％）
6％

10

100％

（－0.9％)
15％

68％

34％

58％
68％

5

60％

17％
0.51

0
1.00

6％
1.43

0.70

8％
1.07

Figure B.18
The value in (F.E.M Dynamic
Analysis Model

Response Deformation
Method) shows the ratio of the load stress of the frame to the wall.

Simple Calculation in Response Deformation Method

Simple Calculation

CASE-1

CASE-2

CASE-356180（1.0％）（kN×103）
Load Force of Frame

13％

Ground Deformation
Wall

Load Shear Force of Lowest Layer30（－0.1％）

Surround Shear Force
Inertia Force of Structure

24％

25

20

（1.7％）（1.2％）
15（－0.7％）
11％

10

5％

84％
70％

42％

100％
76％

5

53％
13％

0
1.00

0.66

0.80

3％

6％

2.88

1.52

Figure B.19
References[B.1] The Architectural Institute of Japan: Design Recommendation for Storage Tanks andTheir Supports, 1990
[B.2] Per B.Schnabel, John Lysmer and H.Bolton Seed : SHAKE-A Computer Program for
Earthquake Response Analysis of Horizontally Layered Sites. EERC 72-12(University ofCalifornia Berkeley), December.1972[B.3] Japan Road Associations：Specifications for Highway Bridges Part Ⅴ Seismic Design.
February.1990 and march. 2002
[B.4] Experimental Study on Dynamic Deformation Characteristics of Soils (Ⅱ) Toshio
Iwasaki, Fumio Tatsuoka, Yoshikazu Takagi. Public Works Research Institute Report
No.153
March.1980

Appendix B

173