Submit manuscript...
MOJ
eISSN: 2572-8520

Civil Engineering

Research Article Volume 3 Issue 1

Dynamic Buckling of Cylindrical Storage Tanks during Earthquake Excitations

Chikhi Abderrahmen, Djermane Mohamed

Reliability Laboratory of Materials and Structures (FIMAS), Tahri Mohamed University, Algeria

Correspondence: Chikhi Abderrahmen, Reliability Laboratory of Materials and Structures (FIMAS), Tahri Mohamed University, Algeria

Received: May 26, 2017 | Published: July 5, 2017

Citation: Abderrahmen C, Mohamed D (2017) Dynamic Buckling of Cylindrical Storage Tanks during Earthquake Excitations. MOJ Civil Eng. 2017;3(1):227-232. DOI: 10.15406/mojce.2017.03.00058

Download PDF

Abstract

The behavior of storage tanks' analysis in seismic areas is of major importance because of the strategic nature of these works. The steel cylindrical tanks are the most susceptible to damage due to dynamic buckling during earthquakes. In this study, three criteria are used to estimate the critical peak ground acceleration caused the tank instability. The liquid inside the tank was modeled using specific Ansys's finite elements and fluid-structure interaction. The calculation includes modal and time history analysis, including material and geometric non-linearity. The result values are compared with standard code previsions as well as the results of previous numerical research, and show the need to improve code provisions.

Keywords: dynamic buckling, tanks, earthquakes, finite element, fluid-structure interaction, instability criteria

Introduction

Cylindrical tank is one of the strategic structures in the daily human life. These facilities are used to store petroleum products, water, oil and chemicals, etc. According to various reports and observations on the structural behavior of the reservoirs during the recent earthquakes, steel tanks are more susceptible to damage than others. Among those reports which cited1 damage of petroleum steel tanks due to the earthquakes of 1933 Long Beach, 1952 Kern County, 1964 Alaska, 1964 Niigata, 1966 Parkfield, 1971 San Fernando, 1978 Miyagi prefecture, 1979 Imperial County, 1983 Coalinga, 1994 Northridge, and 1999 Kocaeli earthquakes,2 Shibata,3 Kono,4 Manos & Clough,5,6 Boumerdès & Sezen.7 Evidence of damage to silos Korso during Boumerdès earthquake7 not provided in seismic zone. Types of failure reported for these structures: diamond or elephant’s foot buckling, uplift of their bases, pipe damage, etc. Among these negative phenomena, dynamic buckling of tank walls remains the most common and more dangerous one.1

The Elephant Foot buckling (EFB)-which is an outward bulge located just above the tank base – results from the combined action of vertical compressive stress, exceeding the critical stress, and hoop tension close to the yield limit.8 The elephant foot buckling bulge usually extends completely around the bottom of tanks due to the reverse in the direction of the seismic excitation.8 The diamond buckling is an elastic instability phenomenon due to the presence of high axial compressive stresses.8–14 The first numerical research estimated the critical PGA that causes the dynamic buckling is due to Virella et al.15 Using the technique of the added mass modeled by ABAQUS code, they examined the PGA which induces the buckling of a set of anchored cylindrical tanks. Followed by Djermane et al.8 presented the contribution to evaluate the PGA values and compared with standard code previsions by using a numerical model. This finite element developed locally is exposed.8 More recently, Sukhvarsh & Mark14 have presented stability static and dynamic analysis of cylindrical tanks. The fluid – structure interaction (FSI) is modeled in ANSYS finite element code. Unlike the most research papers cited above this work uses three different stability criteria to estimate the critical PGA.

Numerical model

In this work, the tank used has been treated by several authors.8,15 The tall tank has a height to a radius ratio (H/R) of 2. In the three-dimensional tank, most developed models use shell element for the wall, the roof and the tank bottom, beam elements for the roof rafters, and fluid element for liquid containing cylindrical tank. As mentioned earlier, the model of the fluid and the shell structure was done using ANSYS software. In the numerical analysis, we used two shell elements for tank structure, the SHELL63 element used for modal analysis and the SHELL181 for time history analysis. This last element is suited for non-linear application. The liquids inside this structure are modeled by FLUID80 element. The roof rafters are modeled by BEAM4 element. Material properties of the studied tank are presented in Table 1, Figure 1 shows the different shell thickness of the tank model.

Water

Density

1000 kg/m3

Bulkmodulus

2.0684e9 Pascal

Viscosity

1.13e-3 N.S/m2

Steel

Density

7850 kg/m3

Poisson ratio

0.3

Elasticity modulus

20.67e10 Pascal

Yield stress

2.5e8 Pascal

Tangent modulus

1.45e9 Pascal

Table 1 Material properties

Figure 1 Tank model with different shell thickness.

Fluid-structure interaction

The fluid-tank interface is modeled by coupling equation; the coincident nodes at the common areas between the fluid element and the shell elements are attached in the normal direction (Figure 2).2

Figure 2 Details of the fluid-structure model.

Boundary conditions

The tank studied in this work is assumed to have a completely fixed-base (all the nodes at the base are assumed to be fixed).

Dynamic analysis

Modal analysis

Modal analysis is used to determine vibration characteristics of this model. The important parameters in the design of a structure for dynamic loading conditions are the natural frequencies and mode shapes.2 The equation of free vibrations for undamped system is given by:

[ M ] { u ..   }   +   [ K ] { u }   =   0 MathType@MTEF@5@5@+= feaagKart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaqcfa4aamWaaO qaaKqzGeaeaaaaaaaaa8qacaWGnbaak8aacaGLBbGaayzxaaqcfa4a aiWaaOqaaKqzGeWdbiaadwhajuaGdaahaaWcbeqaaKqzGeGaaiOlai aac6caaaGaaeiiaaGcpaGaay5Eaiaaw2haaKqzGeWdbiaabccacqGH RaWkcaqGGaqcfa4damaadmaakeaajugib8qacaWGlbaak8aacaGLBb Gaayzxaaqcfa4aaiWaaOqaaKqzGeWdbiaadwhaaOWdaiaawUhacaGL 9baajugib8qacaqGGaGaeyypa0Jaaeiiaiaaicdaaaa@50EE@ (1)

where, [M] = structural mass matrix, [K]= structural stiffness matrix, {u ̈ } nodal acceleration vector, and {u}= nodal displacement vector. For a linear system, free vibration will be expressed as:

u   =   ϕ _ i   c o s ω _ i   t MathType@MTEF@5@5@+= feaagKart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaqcLbsaqaaaaa aaaaWdbiaadwhacaqGGaGaeyypa0Jaaeiiaiabew9aMjaac+facaWG PbGaaeiiaiaadogacaWGVbGaam4CamXvP5wqSX2qVrwzqf2zLnhary GqHrxyUDgaiuaacaWFwaJaeqyYdCNaai4xaiaadMgacaqGGaGaamiD aiaa=Dbmaaa@4FFA@ (2)

where, ϕ _ i MathType@MTEF@5@5@+= feaagKart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaqcLbsaqaaaaa aaaaWdbiabew9aMjaac+facaWGPbaaaa@3A3E@ eigenvector representing the mode shape of the with natural frequency, = ith natural circular frequency in radians per unit time, t = time in seconds. Substitution of Eq. (2) in Eq. (1) gives:

[ M ] + [ K ] ϕ _ i = 0 MathType@MTEF@5@5@+= feaagKart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbiqaaaTdjugibi abgMYiHlabgkHiTKqbaoaadmaakeaajugibiaad2eaaOGaay5waiaa w2faaKqzGeGaey4kaSscfa4aamWaaOqaaKqzGeGaam4saaGccaGLBb Gaayzxaaqcfa4aaaGaaOqaaaGaayPkJaqcLbsacqaHvpGzcaGGFbGa amyAaiabg2da9iaaicdaaaa@4AD0@ (3)

Nonlinear time history analysis

According to various reports,1,9 the damage in the thin tank walls in seismic areas shows almost all large displacements and relatively large deformations. For these raisons, the material and geometric non-linearities are considered in this work. The elasto-plastic stress-strain curve is used for steel. The software offers several options, among which we have chosen the simplest consisting of a bilinear kinematic hardening curve, taking into account the Beauschinger effect. To obtain a satisfactory solution, the loads are applied incrementally at each time step through a series of smaller sub-step.10

The transient dynamic analysis solves the basic equation of motion:

[ M ] { u .    }   + [ C ] { u   ̇   }   +   [ K ] { u }   =   { F ( t ) }   ( 4 ) MathType@MTEF@5@5@+= feaagKart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaqcfa4aamWaaO qaaKqzGeaeaaaaaaaaa8qacaWGnbaak8aacaGLBbGaayzxaaqcfa4a aiWaaOqaaKqzGeWdbiaadwhajuaGdaahaaWcbeqaaKqzGeGaaiOlaa aacaqGGaGaaeiiaaGcpaGaay5Eaiaaw2haaKqzGeWdbiaabccacqGH RaWkjuaGpaWaamWaaOqaaKqzGeWdbiaadoeaaOWdaiaawUfacaGLDb aajuaGdaGadaGcbaqcLbsapeGaamyDaiaabccacqWIhmGjcaqGGaaa k8aacaGL7bGaayzFaaqcLbsapeGaaeiiaiabgUcaRiaabccajuaGpa WaamWaaOqaaKqzGeWdbiaadUeaaOWdaiaawUfacaGLDbaajuaGdaGa daGcbaqcLbsapeGaamyDaaGcpaGaay5Eaiaaw2haaKqzGeWdbiaabc cacqGH9aqpcaqGGaqcfa4damaacmaakeaajugib8qacaWGgbqcfa4d amaabmaakeaajugib8qacaWG0baak8aacaGLOaGaayzkaaaacaGL7b GaayzFaaqcLbsapeGaaiiOaKqba+aadaqadaGcbaqcLbsapeGaaGin aaGcpaGaayjkaiaawMcaaaaa@6A70@

where: [C] = damping matrix, { u   ̇   } MathType@MTEF@5@5@+= feaagKart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaqcfa4aaiWaaO qaaKqzGeaeaaaaaaaaa8qacaWG1bGaaeiiaiablEWaMiaabccaaOWd aiaawUhacaGL9baaaaa@3C88@ = nodal velocity vector and { F ( t ) } MathType@MTEF@5@5@+= feaagKart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaqcfa4aaiWaaO qaaKqzGeaeaaaaaaaaa8qacaWGgbqcfa4damaabmaakeaajugib8qa caWG0baak8aacaGLOaGaayzkaaaacaGL7bGaayzFaaaaaa@3E1A@ = load vector.

The earthquake records selection

The present study has used earthquake records of El-Centro and San Fernando with maximum acceleration of 0.349 g and 1.2 g respectively. A record is considered interesting in the sense of this study if the maximum amplitudes are recorded during the first seconds which could also contain the essential of the frequency content. This is the case with the El Centro 1940 earthquake. The selection of the San Fernando 1970 earthquake is due essentially to its hard PGA which caused massive damages to storage structures (Figure 3).8,17


Figure 3 Selected accelerograms: (a) San Fernando 1971, (b) El Centro 1940.

Presentation of stability criteria

Total energy-phase plane criterion (Hoff & Bruce, 1954)

The curve representing the movement is traced in a phase plane (U;U ̇) as shown in Figure 4. The stable movements are characterized by limited trajectories and do not move too much away from the solution of the static equilibrium which plays the role of a center of attraction. When the load reaches the critical value, the trajectory moves away from this pole without any oscillation around it.6

Figure 4 Phase plane diagram before and after the Pcr6.

Equation of motion criterion (Budiansky & Roth4)

The first and most used criterion of stability is due to Budiansky & Ruth.4 It was formulated as an engineering application of the Liapunov stability criteria. In this criterion, the time displacement curve is plotted for several values of the PGA. The PGA value corresponding to a curve which gives a ‘‘jump’’ relatively to its neighboring curves indicates the dynamic buckling critical value(PGA) cr Figure 5A.

Figure 5 Critical load: a- criteria of Budiansky & Roth; b- criteria of Ari Gur & Simonetta.

Criteria of Ari Gur & Simonetta3

This curve has been proposed originally by Budiansky & Ruth4 and generalized later by Ari-Gur & Simonetta.3 In the case of seismic excitation, this strong increase is less pronounced, and is rather replaced by a change of the slope of the pseudo-dynamic curve relating the PGA intensity and the maximum radial displacement at a fixed monitored node (several nodes are considered especially in the potential buckling zones) as shown in Figure 5B. This criterion offers an estimation which must be confirmed with the above criteria.6 The most commonly used stability criteria such as that of Budiansky-Ruth, the phase plane and energy witch ‘‘watch’’ the dynamic bifurcation are very time consuming CPU.8

Dynamic buckling of tanks: standard provisions

Based on simplified analogical models and on experimental research, different simplified design standards for anchored and unanchored tanks have been developed.8 The normal design procedure in North America for water structures is based on the specifications provided by the American Water Work Association (AWWA-D100)1 and for other fluids on the American Petroleum Industry (API) Standard 650.16 More recently, the European code EC 8-Part 49 has provided enhanced analysis methods for anchored and unanchored tanks design. Two different strategies have been used in the code practice: The simple stress limitation method adopted in many standards especially in American standards (AWWA100,1 API65016), the semi-empirical method of Rotter18 adopted by the New Zealand seismic Tank design recommendation, the European Standards ENV1998-4,9 etc.

Modal analysis results

To ensure the validity of the used model, modal analysis for the tank with FSI system was performed. The results were compared with previous research and standard code provisions. The results are summarized in the Table 3.The excellent agreements indicate the validity of the used model.

mode

Frequency

Period

Cumulative Mass Fraction

convective

1st

2

0.168578

5.932

0.26671

2nd

6

0.264308

3.7835

0.273556

3th

13

0.302594

3.3048

0.274772

4th

22

0.321403

3.1114

0.275144

Impulsive

1st

281

3.16079

0.31638

0.969402

2nd

340

6.95603

0.14376

0.970209

3rd

341

6.97019

0.14347

0.999052

4th

444

45.0829

0.0222

0.999978

Table 2 The significant convective and impulsive modes.

Ansys

Virella et al. [15]

Eurocode 8

API 650

Ratio

n Mode

Type

Timp (sec)

0.316

0.3

0,299

/

0.944

281

column

Tcon (sec)

5. 932

/

5.778

5.696

0.96

2

/

Table 3 Modal Period Comparison.9,15-16

Participation factor

Selection of significant modes is based on the criterion of effective mass participation; Figure 6 represents the actual participation factor and the number of mode. We can see that natural impulsive frequencies are greater than natural convective frequencies. Table 3 compares the obtained results versus those given by the EC8, API 650 standard and the attached masses model. This table indicates that estimated values for structural periods of vibration are reasonable. The obtained fundamental impulsive mode is a column mode type as shown in Figure 7B.

Figure 6 Modal participating mass ratio.

Figure 7 Fundamental modes' shapes of a tank model: a-the first convective mode. b- the first impulsive mode.

Dynamic buckling results

Elcentro earthquake

Figure 8 shows the pseudo equilibrium path for El Centro excitation. The discontinuity on the curve indicates that the PGA critical value (PGAcr) occurs at 0.384g. Figure 9 shows several history curves corresponding to different excitation levels. The reversal nature of earthquake excitation caused a hard difficulty to determine (PGA)cr with Budiansky-Ruth criterion type, this is clearly showed in Figure 9.8 At level 0.384, a disproportionate increase in displacements is distinguished.

Figure 8 El Centro Pseudo-dynamic path for the tank Model.


Figure 9 Time history curves before and after PGAcr for the tank Model under El Centro earthquake.

This increase does not correspond to a monotonic jump for the above-mentioned reasons. The phase planes' criterion illustrated in Figure 10 shows more easily in this case the instability near the (PGA) cr. The difficulty for using this criterion is, in some cases, the same as that reported for Budiansky and Ruth one, but the use of the two criteria simultaneously can be more illustrative. The nature of the obtained dynamic buckling response can be evaluated by studying the deformation and the iso-stresses around the excitation critical level see Figure 11.


Figure 10 Phase plane before and after PGAcr for the tank model under El Centro earthquake.

Figure 11 San Fernando Pseudo-dynamic path for the tank Model.

San Fernando earthquake

Figure 11 gives the PGAcr for San Fernando excitation. Using an estimation given by the pseudo dynamic path, the Budiansky–Ruth and phase plane criteria are then used to confirm the obtained value Figures 12 & 13. The dynamic buckling obtained is of the elastic type, the folds are concentrated on the upper part of the tank see Figure 14.


Figure 12 Time History curves before and after PGAcr for the tank Model A under San fernando earthquake.


Figure 13 Phase plane before and after PGAcr for the tank Model under San Fernando earthquake.

Figure 14 Von miss stress for the studied tank at time= 4.54 seconds under El Centro earthquake at PGA = 0.384g under PGAcr.

This type of buckling has been obtained by several authors Natsiavas & Babcock13 and also by Virella et al.15 Natsiavas and Babcock13 studied and characterized this type of buckling. They showed, in particular, that the dynamic pressure in the fluid exerts a negative pressure on the zone of the wall located near the fluid free surface, where the hydrostatic pressure and the wall thickness are of weak values. The resulting negative pressure (hydrostatic pressure- impulsive pressure) can cause sufficiently large compressive stresses to generate localized elastic buckling (depression buckling).

Table 4 shows the critical values under El Centro and San Fernando earthquakes at PGA of 0.384g, and 0.416g respectively. Through the obtained results see Tables 4 for the dynamic buckling of the tank under two seismic excitations, one can conclude the following points:

  1. In the El Centro earthquake, the critical values (the radial displacement at the tank wall, stresses and hydrodynamic pressure) occurred at time of 4.54 seconds, and the maximum water displacement value was late at time of 8.64 seconds.
  2. Figure 15 & 16 shows the maximum water vertical Displacement under El Centro earthquake.
  3. In the San Fernando earthquake, these critical values occurred late and are synchronized with maximum water displacement at time of 8.23 seconds. In this earthquake excitation, the peak value of water sloshing is greater than the value obtained for the El Centro excitation. This effect cannot be seen in the case of using the technique of the added masses.
  4. The maximum vertical surface displacement can contribute in the collapse of the tank roof. Table 5 compares the maximum free surface displacement obtained for spectrum analysis of the used model versus those given by the EC8 and API 650. This table indicates that estimated values from spectrum analysis are reasonable. In this analysis, we have used the EC8 spectrum curve.
  5. Table 6 compares the (PGA)cr values obtained by the numerical model versus those estimated by EC8 and API provisions. For the two earthquakes excitations, the results of the numerical model and the two codes are quite similar.

El Centro Earthquake at PGA=0.384g

San fernando Earthquake at PGA=0.416g

Displacement of tank wall (cm)

(MPa)

(MPa)

Hydrodynamic pressure (MPa)

free surface displacement (cm)

Displacement of tank wall (cm)

(MPa)

(MPa)

Hydrodynamic pressure (MPa)

2.483

139

84.15

23.14

38.14

5.043

91.2

81.5

28.7

Table 4 The critical values of the radial displacement, stresses, hydrodynamic pressure and maximum water displacement under El Centro and San fernando earthquakes

Figure 15 Maximum water vertical. Displacement at time= 8.64 seconds under El Centro earthquake at PGA = 0.384g.

Figure 16 Dynamic buckling obtained in the upper. zone of the tank wall under San fernando earthquake at PGA = 0.466g after PGAcr.

Spectrum Analysis of Used Model

Eurocode 8

API 650/AWWA

27.52

27.43

26.86

Table 5 The maximum free surface displacement comparison

Earthquake

Numerical Model

By Codes

AWWA

Difference

EC8

Difference

El Centro

0.384g

0.399g

3.76%

0.387g

0.77%

San Fernando

0.400g

0.443g

9.70%

0.430g

6.97%

Table 6 Comparison of (PGA) cr values for the Tall tank

Conclusion

In this work, our investigation was focused on the search for the value of critical PGA that can lead to instability of the tank. We used a finite-element model considering the fluid-structure interaction, to have reliable results. We have used three different instability criteria. The results obtained by modal analysis confirm the validity of our model. The results of the instability study were compared to the values given by the international codes. The maximum free surface displacements obtained for spectrum analysis of the used model versus those given by standards codes are quite similar. The results of the instability study were compared to the values given by the international codes. It is clear from this comparison that the code results are not conservative.

Acknowledgements

None.

Conflict of interest

The author declares no conflict of interest.

References

  1. American Water Work Association. Welded steel tanks for water storage. Colorado, USA; 1997.
  2. The ANSYS Structural Software System. ANSYS INC; 2009. 12 p.
  3. Ari–Gur J, Simonetta SR. Dynamic pulse buckling of rectangular composite plates. Composites Part B: Engineering. 1997;28:301–308.
  4. Budiansky B, Roth S. Axisymmetric dynamic buckling of clamped shallow spherical shells. NASA collected papers on stability of shells structures. 1962;(TN–1510):597–606.
  5. Djermane M, Chelghoum A. Nonlinear dynamic and dynamic buckling of thin shells using a drilling finite element. International Journal of Advanced Manufacturing Technology. 2007;(1–2):175–184.
  6. Djermane M. Dynamic Buckling of Thin shells in Seismic Zones. Algeria: Doctoral thesis; 2008.
  7. Djermane M. Dynamic buckling of thin shells: evaluation of criteria. International Review of Mechanical Engineering. 2011;3(5):449–453.
  8. Djermane M. Dynamic buckling of steel tanks under seismic excitation: Numerical evaluation of code provisions. Engineering Structures. 2014;70:181–196.
  9. European Committee for Standardization. Eurocode 8: Design provisions of earthquake resistance of structures. Part 4: Silos, tanks and pipelines. Brussels; 2006.
  10. He Liu, Schubert DH. Effects of nonlinear geometric and material properties on the seismic response of fluid/tank systems. USA: University of Alaska Anchorage; 1963.
  11. Housner GW. The dynamic behavior of water tanks. Bulletin of the Seismological Society of America. 1963;2(53):381–387.
  12. Malhotra PK, Wenk T, Wieland M. Simple Procedure for Seismic Analysis of Liquid–Storage Tanks. Structural Engineering International. 3rd 2000. p. 197–201.
  13. Natsiavas S, Babcock C. Buckling at the top of a fluid–filled tank during base excitation. ASME Journal of Pressure Vessel Technology. 1987;4(109):374–380.
  14. Sukhvarsh J, Mark L. Stability analysis of cylindrical tanks under static and earthquake loading. Journal of Civil Engineering and Architecture. 2015;9:72–79.
  15. Virella JC, Godoy LA, Suárez LE. Dynamic buckling of anchored steel tanks subjected to horizontal earthquake excitation. Journal of Constructional Steel Research. 2006;(62):521–531.
  16. API–650. Welded Steel Tanks for Oil Storage. USA: American Petroleum Institute; 2005.
  17. Yaser Z. Stability of cylindrical oil storage tanks during an earthquake. International Journal of Engineering & Technology IJET–IJENS. 2013;3(13):78–83.
  18. Rotter J. Shell structures: the new european standard and current research needs. Thin–Walled Structures. 1998;(31):3–23.
Creative Commons Attribution License

©2017 Abderrahmen, et al. This is an open access article distributed under the terms of the, which permits unrestricted use, distribution, and build upon your work non-commercially.