Thin-Walled Structures: Amir R. Shokrzadeh, Mohammad R. Sohrabi
Thin-Walled Structures: Amir R. Shokrzadeh, Mohammad R. Sohrabi
Thin-Walled Structures
journal homepage: [Link]/locate/tws
art ic l e i nf o a b s t r a c t
Article history: Steel cylindrical tanks are formed by very thin-walled shells and as typical thin-walled structures, they
Received 11 February 2016 are very sensitive to buckling under external pressures, especially when they are empty or at low liquid
Received in revised form level. Results of numerical investigations on the effects of spiral stairway on the buckling behavior of
25 May 2016
ground based steel cylindrical liquid storage tanks subjected to both wind and vacuum pressures are
Accepted 25 May 2016
Available online 1 June 2016
presented. It is concluded that by choosing appropriate circumferential position for stairway, considering
regional prevailing wind direction, designers can use the noticeable amount of added strength arising
Keywords: from stairway as a safety margin. In other words, it seems that the spiral stairway acts as an oblique
Steel cylindrical tanks stiffener on the tank wall. Contrary to the case of wind loading, the stairway has negligible effect on
Buckling
buckling resistance of tanks under vacuum pressure. However, it changes the buckling modes of tank
External pressures
significantly under vacuum loading.
Numerical investigations
Spiral stairway & 2016 Elsevier Ltd. All rights reserved.
Prevailing wind direction
[Link]
0263-8231/& 2016 Elsevier Ltd. All rights reserved.
438 A.R. Shokrzadeh, M.R. Sohrabi / Thin-Walled Structures 106 (2016) 437–447
Table 1 Table 2
Geometrical characteristics of tanks TK-a, TK-b and TK-c. Shell courses details for tank TK-2047.
Geometrical characteristics TK-a TK-b TK-c Course no.a Course thickness (mm) TK-2047 Course width (m) TK-2047
a
Courses are numbered from bottom to top (i.e. no.1 refers to bottom shell
course etc.).
2. Computational models
The tanks considered in this paper have cone roofs supported Fig. 3. Finite element model of tank (TK-a).
by a number of columns and radial beams. The tanks were as-
sumed to be fixed at the base. Since only fixed base tanks are
considered, the bottom plates are not included in the models. The
cylinder and they are placed on ring beams.
radial beams are connected at their ends directly to the tank The tank geometries chosen (TK-a, TK-b and TK-c) are similar to
those considered by Virella et al. [19]. Table 1 summarizes the
geometries of tanks TK-a, TK-b and TK-c, which are representative
of three tanks with aspect ratios of 0.40, 0.63 and 0.95,
respectively.
Fig. 2. Typical structure of rafter supported conical roof steel tanks. Fig. 4. Spiral stairway of tank TK-2047.
A.R. Shokrzadeh, M.R. Sohrabi / Thin-Walled Structures 106 (2016) 437–447 439
Fig. 5. Spiral stairway of TK-2047 details: (a) the real stairway and (b) simplified model of the stairway (beam model).
Fig. 6. Model of the tank with spiral stairway (TK-2047s) adopted for the computations.
In addition, a tank (TK-2047) in an oil plant located 20 km west above, the conical roof is supported by a number of radial rafters
of Arak (Markazi province, Iran) is considered (see Fig. 1). TK-2047 and ring beams with three groups of sections (inner rafters:
supports a conical roof with a slope equal to 1/16. As mentioned H198 99 4.5 7, outer rafters: H248 124 5 8 and ring
440 A.R. Shokrzadeh, M.R. Sohrabi / Thin-Walled Structures 106 (2016) 437–447
Table 3
Critical buckling loadsa obtained from bifurcation (eigenvalue) analysis.
1st 2nd 3rd 4th 5th 6th 7th 8th 9th 10th
TK-2047 2.357 2.358 2.359 2.360 2.366 2.366 2.371 2.372 2.386 2.387
TK-2047s 2.356 2.357 2.361 2.362 2.369 2.369 2.378 2.379 2.393 2.397
b
Relative difference (%) 0.02 0.04 0.068 0.064 0.127 0.101 0.266 0.295 0.285 0.39
TK-a 2.096 2.096 2.099 2.099 2.135 2.135 2.139 2.139 2.197 2.197
TK-as 2.097 2.097 2.112 2.115 2.14 2.142 2.173 2.178 2.215 2.229
b
Relative difference (%) 0.062 0.072 0.634 0.762 0.225 0.290 1.627 1.814 0.842 1.424
TK-b 2.999 2.999 3.010 3.010 3.021 3.022 3.045 3.046 3.100 3.100
TK-bs 3.004 3.004 3.019 3.020 3.043 3.046 3.077 3.084 3.120 3.130
b
Relative difference (%) 0.167 0.177 0.302 0.312 0.728 0.798 1.051 1.248 0.648 0.984
TK-c 2.995 2.995 3.007 3.007 3.009 3.009 3.045 3.045 3.051 3.051
TK-cs 2.999 3.0 3.013 3.013 3.033 3.034 3.062 3.064 3.098 3.103
b
Relative difference (%) 0.140 0.150 0.193 0.206 0.804 0.831 0.565 0.627 1.547 1.717
a
All buckling loads are in kPa.
λ s −λc
b
Relative difference¼ c × 100% .
λc
λcs −λc
Fig. 7. Relative difference of buckling loads ( ×100%) for lowest ten eigenvalues obtained from bifurcation analysis: (a) TK-2047 under vacuum pressure, and (b) TK-a,
λc
TK-b and TK-c under vacuum pressure.
beams: H496 199 9 14). The overall geometry is given by a element discretization of the structure. Eight-node shell elements
diameter D¼ 36.59 m, cylinder height H¼ 14.63 m, and 1.143 m with reduced integration (identified as S8R5 in ABAQUS) were
maximum elevation of the roof with respect to the cylinder. It can used for the cylinder, whereas triangular elements were employed
be seen that the aspect ratio of the tank (H/D) is quite low (0.40).
on the roof (STRI65). Rafters were modeled by beam elements. In
The tank has variable wall thicknesses and its cylindrical shell is
divided in six courses (see Fig. 2). The width and thickness of each the finite element model, columns are represented by vertical
shell course along with relevant information regarding the roof are restraints were applied at the points where they support the
summarized in Table 2. structural grid (some points on ring beams). The fixed boundary
In all the models, the tank was assumed empty, this is worst condition was applied on the bottom of the finite element model.
case for tanks subjected to wind and uniform external pressure
In all computational models, the material is assumed as elastic and
[6,18,20].
isotropic. ASTM A36 steel is assumed for all tank models, with
Fig. 8. First buckling mode of TK-as, TK-bs and TK-cs under vacuum pressure (isometric and bottom views): (a) isometric view of TK-as; (b) bottom view of TK- as;
(c) isometric view of TK-bs; (d) bottom view of TK- bs; (e) isometric view of TK-cs; (f) bottom view of TK- cs.
442 A.R. Shokrzadeh, M.R. Sohrabi / Thin-Walled Structures 106 (2016) 437–447
Fig. 10. Created tank models and their associated wind directions (the figure shows
the tanks from top view).
λ s −λc
Fig. 12. Relative difference of buckling loads ( c ×100%) : (a) TK-2047 under
λc
various wind directions, and (b) TK-a, TK-b and TK-c under various wind directions.
Table 3 and Fig. 7 show the eigenvalues and the relative dif-
ference of the first ten eigenvalues of the tanks with and without
stairway. Based on the results, generally, the critical buckling loads
obtained from eigenvalue analysis of tanks without stairway are
slightly lower than those obtained from eigenvalue analysis of
tanks with stairway, as depicted in Table 3. The relative difference
of the first ten eigenvalues of the tanks with and without stairway
varies from 0.017 to 1.8%. Moreover, the maximum difference be-
Fig. 11. Bifurcation analysis results of tanks: (a) critical loads for the lowest six tween the first eigenvalues is observed in the TK-b, which
eigenvalues of TK-2047s in different wind directions; (b) a comparative view of the corresponds to the tank with aspect ratio of 0.63.
TK-as, TK-bs and TK-cs buckling loads in different wind directions.
Fig. 8 shows the first buckling mode of tanks with stairway
under uniform pressure. Although in both cases (with, and with-
out stairway) the corresponding eigenvalues are the same (see
treads (made from galvanized grating with 32 mm thickness) is Fig. 7), the buckling modes attached to them are different. Tanks
ignored. In addition, handrail is not considered in the model without stairway show approximately the same buckling de-
(handrail stiffness is ignored). However, its weight along with formations all around the circumference of tank wall, whereas in
treads weight are taken into account by increasing density of beam tanks with stairway, the buckling waves are not formed in loca-
material. Since the stiffness of the beam stairway is smaller than tions around the stairway, as shown in Fig. 8. The deflected shape
the real one, the effect of the beam stairway on the buckling be- of tanks with spiral stairway resulting from the computations
havior (including buckling load, mode shapes etc.) can be defined (Fig. 8) is also consistent with the buckled shapes shown in Fig. 9.
as lower bound for the real stairway effect. As can be seen from Fig. 9, the buckling waves are not formed in
Furthermore, the finite element models of TK-a (H/D ¼ 0.4), TK- locations around the stairway, which confirms analytical results.
b (H/D ¼0.63) and TK-c (H/D ¼0.95) with and, without stairway Fig. 4(a) in Ref. [5] also shows a buckled tank in which the stairway
are examined (Models A, B and C of Ref. [19]) under vacuum and has locally strengthened the shell, thus causing a non-symmetric
wind pressures. The same slope angle (approximately 44°) and mode.
444 A.R. Shokrzadeh, M.R. Sohrabi / Thin-Walled Structures 106 (2016) 437–447
Fig. 13. First buckling mode of finite element models (a) TK-a; (b) TK-as-(10°); (c) TK-b; (d) TK-bs-(20°); (e) TK-c; (f) TK-cs-(35°).
λcos −λco
Fig. 14. Relative difference of buckling loads ( ×100%) for tanks with different aspect ratios.
λco
3.2. Tanks under wind pressure measured with respect to the line connecting the center of the
stairway to the tank cylinder axis perpendicularly. The wind di-
To simulate effect of wind direction on tanks with spiral stair- rection angles used for analysis are shown in Fig. 10.
way, the direction of wind incidence is considered as a variable. Clearly, the spiral stairway gives a positive reinforcement to the
For this purpose, wind directions characterized by angles ranging tank wall (see Fig. 11). The λ c=3. 05, 2. 87, 4. 00and3. 86 kPa are
between 0° and 360° are investigated. Wind direction angle (α) is the smallest eigenvalues of TK-2047, TK-a, TK-b and TK-c
A.R. Shokrzadeh, M.R. Sohrabi / Thin-Walled Structures 106 (2016) 437–447 445
Fig. 16. First buckling mode of finite element models (a) TK-ao; (b) TK-aos-(10°); (c) TK-bo; (d) TK-bos-(20°); (e) TK-co; (f) TK-cos-(35°).
perpendicularly is 15° to 25° for TK-as, 15° to 40° for TK-bs under uniform external pressure and it can cause a dramatic
and 5° to 60° for TK-cs, the buckling resistance of the conical decline of buckling deformation in locations around the
roof tanks increases 4%–13.5%, 3%–21.3% and 3%–22%, respec- stairway.
tively. These values change to 4%–12%, 12%–42% and 20%–73% for 8. Since the stiffness of the beam stairway (constructed by beam
TK-aos, TK-bos and TK-cos (open-top tanks), respectively. There- elements in numerical model) is smaller than that of real one,
fore, as the aspect ratio of the tanks increases, the stairway can the effect of the model stairway on the buckling behavior (in-
strengthen larger portion of the tank wall. In addition, by cluding buckling load, mode shapes etc.) can be defined as a
choosing appropriate circumferential position for spiral stair- lower bound for real stairway effect. Therefore, this kind of
way considering regional prevailing wind direction, designers simplification in spiral stairway modeling can be used in prac-
can use this added strength as a safety margin. tical design of steel cylindrical tanks.
5. Since the fatigue failure happens as a result of repetition of
maximum stress experienced by the specimen, by choosing
appropriate circumferential position for spiral stairway con-
References
sidering regional prevailing wind direction, number of experi-
enced maximum stress can be decreased considerably. There-
[1] P. Pavlovic, R. Folic, V. Radonjanim, M. Tatomirovic, The testing and repair of
fore, spiral stairway can improve fatigue life. steel silo, Constr. Build. Mater. 11 (5) (1997) 353–363.
6. Based on the results, the strengthening effect of the circumfer- [2] F.G. Flores, L.A. Godoy, Buckling of short tanks due to hurricanes, Eng. Struct.
ential stairway is significant and should be considered in future 20 (8) (1998) 752–760.
[3] R.C. Jaca, L.A. Godoy, Wind buckling of metal tanks during their construction,
editions of design codes.
Thin-Walled Struct. 48 (6) (2010) 453–459.
7. Contrary to the case of wind loading, the stairway has negligible [4] A. Zingoni, Liquid-containment shells of revolution: a review of recent studies
effect on buckling resistance of tanks under vacuum pressure. on strength, stability and dynamics, Thin-Walled Struct. 87 (2015) 102–114.
However, it changes the buckling modes of tanks significantly [5] L.A. Godoy, Buckling of vertical oil storage steel tanks: review of static buck-
ling studies, Thin-Walled Struct. 103 (2016) 1–21.
A.R. Shokrzadeh, M.R. Sohrabi / Thin-Walled Structures 106 (2016) 437–447 447
[6] Y. Zhao, Y. Lin, Buckling of cylindrical open-topped steel tanks under wind walled cylindrical tanks under wind load, Thin-Walled Struct. 31 (1998)
load, Thin-Walled Struct. 79 (2014) 83–94. 203–220.
[7] G. Portela, L.A. Godoy, Wind pressures and buckling of cylindrical steel tanks [16] L. Chen, J.M. Rotter, Buckling of anchored cylindrical shells of uniform thick-
with conical roof, J. Constr. Steel Res. 61 (6) (2005) 786–807. ness under wind load, Eng. Struct. 41 (2012) 199–208.
[8] G. Portela, L.A. Godoy, Wind pressures and buckling of cylindrical steel tanks [17] F. Bu, C. Qian, A rational design approach of intermediate wind girders on large
with dome roof, J. Constr. Steel Res. 61 (6) (2005) 808–824. storage tanks, Thin-Walled Struct. 92 (2015) 76–81.
[9] C. De Paor, D. Kelliher, K. Cronin, W.M.D. Wright, S.G. McSweeney, Prediction [18] C.A. Burgos, J.C. Batista-Abreu, H.D. Calabro´, R.C. Jaca, L.A. Godoy, Buckling
of vacuum-induced buckling pressures of thin-walled cylinders, Thin-Walled estimates for oil storage tanks: Effect of simplified modeling of the roof and
Struct. 55 (2012) 1–10. wind girder, Thin-Walled Struct. 91 (2015) 29–37.
[10] Y. Lin, Y. Zhao, Wind loads on fixed-roof cylindrical tanks with very low aspect [19] J.C. Virella, L.A. Godoy, L.E. Suarez, Fundamental modes of tank-liquid systems
ratio, Wind Struct. 18 (6) (2014) 651–668. under horizontal motions, Eng. Struct. 28 (2006) 1450–1461.
[11] L. Chen, J.M. Rotter, C. Doerich-Stavridis, Practical calculations for uniform [20] C. Maraveas, G.A. Balokas, K.D. Tsavdaridis, Numerical evaluation on shell
external pressure buckling in cylindrical shells with stepped walls, Thin- buckling of empty thin-walled steel tanks under wind load according to
Walled Struct. 61 (2012) 162–168. current American and European design codes, Thin-Walled Struct. 95 (2015)
[12] F.G. Flores, L.A. Godoy, Buckling of short tanks due to hurricanes, Eng. Struct. 152–160.
20 (8) (1998) 752–760. [21] ABAQUS v6.9. Simulia. Unified FEA. Dassault Systems. Warwick, Rhode Island,
[13] L.A. Godoy, Performance of storage tanks in oil facilities following hurricanes USA, 2009.
Katrina and Rita, ASCE J. Perform. Constr. Facil. 21 (6) (2007) 441–449. [22] AS/NZS. Structural design actions, Part 2 – Wind Actions, Australian/New
[14] M.A. Al-Kashif, H. Ramadan, A. Rashed, M.A. Haroun, Effect of FRP composites Zealand Standards 1170.2:2011, Sydney, 2011.
on buckling capacity of anchored steel tanks, Steel Comp. Struct. 10 (4) (2010) [23] C. De Paor, The Effect of Random Geometric Imperfections on the Buckling of
361–371. Thin Cylindrical Shells Due to External Pressure (Ph.D. thesis), National Uni-
[15] H. Schmidt, B. Binder, H. Lange, Postbuckling strength design of open thin versity of Ireland, Cork, Ireland, 2012.