Mathematical Problems in Engineering

Mathematical Problems in Engineering / 2019 / Article

Research Article | Open Access

Volume 2019 |Article ID 6190159 |

Sangmok Lee, Byungmin Kim, Young-Joo Lee, "Seismic Fragility Analysis of Steel Liquid Storage Tanks Using Earthquake Ground Motions Recorded in Korea", Mathematical Problems in Engineering, vol. 2019, Article ID 6190159, 15 pages, 2019.

Seismic Fragility Analysis of Steel Liquid Storage Tanks Using Earthquake Ground Motions Recorded in Korea

Academic Editor: Francisco J. Montáns
Received09 Apr 2019
Revised22 Jun 2019
Accepted27 Jun 2019
Published15 Jul 2019


Liquid-containing storage tanks are important structures in industrial complexes. Because earthquake damages to liquid storage tanks can cause structural collapse, fires, and hazardous material leaks, there have been continuous efforts to mitigate earthquake damages using seismic fragility analysis. In this regard, this study focuses on the seismic responses and fragility of liquid storage tanks. First, the characteristics of earthquake ground motions are a critical factor influencing the seismic fragility of structures; thus, this study employs real earthquake records observed in the target area, southeastern Korea, with the earthquake characteristics estimated based on the ratio of peak ground acceleration to peak ground velocity. When a liquid storage tank oscillates during an earthquake, additional forces can impact the tank wall owing to hydrodynamic pressures. Therefore, this study presents a sophisticated finite element (FE) model that reflects the hydrodynamic effect of an oscillating liquid. Another advantage of such an FE model is that detailed structural responses of the entire wall shells can be estimated; this is not possible in simplified lumped mass or surrogate models. Lastly, probabilistic seismic demand models are derived for three critical limit states: elastic buckling, elephant’s foot buckling, and steel yielding. Using the real earthquake ground motion records, constructed FE model, and limit states, a seismic fragility analysis is performed for a typical anchored steel liquid storage tank in Korea. In addition, for comparison purposes, a ring-stiffened model is investigated to derive a seismic fragility curve. The results of the seismic fragility assessment show that elastic buckling is the most vulnerable damage state. In contrast, elephant’s foot buckling and steel yielding indicate relatively severe damage levels. Furthermore, it is observed that ring stiffeners decrease the elastic buckling damage, although there is no practical effect on elephant’s foot buckling and steel yielding in all ground motion intensities.

1. Introduction

Liquid storage tanks are a major component of industrial complexes. When a storage tank containing liquid is subjected to seismic loads, it can be damaged and collapsed, potentially causing the leakage of hazardous materials and major fires. According to Persson and Lönnermark [1], the 1964 Niigata earthquake in Japan led to fires that resulted in the loss of 97 tanks containing 1.1 million barrels of crude. It was also reported that 29 tanks were damaged or leaked owing to the impact of the 2003 Hokkaido earthquake. In Turkey, three naphtha tanks and a crude oil tower were damaged owing to the 1999 İzmit earthquake.

According to statistical data obtained from the Korea Industrial Complex Corporation [2], there are 1189 industrial complexes in Korea; in particular, there are 436 industrial complexes in southeastern Korea, as of 2018. In addition, earthquakes have frequently occurred in the southeastern part of Korea in recent years, exposing these many industrial complexes to seismic hazards. Thus, the importance of structural risk assessment has become increasingly important, and there have been continuous efforts to mitigate earthquake damages using seismic fragility analyses [3, 4]. This study concentrates on the structural fragility of liquid storage tanks.

The characteristics of earthquake ground motions are also a critical factor in the seismic fragility of structures. Previous studies [36] investigated the seismic fragility of structures with various types of earthquake ground motions to determine the effects of ground motion characteristics. The analysis results of these studies showed that earthquake damages to structures varied with different ground motion properties. In addition, for an accurate seismic fragility analysis of structures, numerous studies on seismic fragility assessments using analytical models [310] have utilized real ground motion records from world-wide well-known earthquakes. However, these records may not reflect the characteristics of a specific region. In this study, real ground motion records observed in the target area (i.e., southeastern Korea) are intensively collected and used in seismic fragility estimation for a more realistic simulation.

When a steel liquid storage tank is exposed to an earthquake, additional forces can impact the tank wall owing to hydrodynamic pressures. The hydrodynamic response of the liquid in a tank can be described by two distinct modes: convective and impulsive [11]. The upper part of the liquid, having convective motion (often called the sloshing), exhibits a long period oscillation, whereas the lower part of the liquid moves together with the tank wall as though rigidly attached to the tank. Veletsos [12] investigated the seismic effects of a flexible liquid storage tank and found that hydrodynamic forces could have a greater influence on a flexible tank than on a rigid one.

Dynamic buckling mainly appears in an oscillating liquid-filled tank owing to the liquid-tank interaction, and many researchers have explored the buckling of steel tanks subjected to seismic loads with numerical models [9, 1315]. These studies have demonstrated the need for a sophisticated model of structural analysis that can account for hydrodynamic pressures when analyzing the seismic response and fragility of a liquid storage tank.

Seismic fragility analyses for liquid storage tanks have been performed in recent years. Buratti and Tavano [9] derived seismic fragility curves for the buckling of the tank wall using an added-mass approximate model for liquid modeling. Paolacci et al. [16] concentrated on the seismic fragility of a liquid storage tank in terms of the buckling of the tank wall with simplified numerical models. Phan et al. [10] used lumped mass models and considered the buckling failure of the tank wall for seismic fragility assessment. However, the aforementioned studies considered the liquid-tank interaction using simplified tank models in their seismic fragility analyses. Lee et al. [6] developed an advanced FE model to take into account hydrodynamic pressures, but the buckling failure was not considered. In addition, these studies utilized worldwide well-known earthquake records in their fragility assessment.

In this study, the seismic risk of a typical anchored cylindrical storage tank is investigated, and seismic fragility curves are derived. Real ground motion data reflecting the characteristics of the target region are collected and used. To consider hydrodynamic effect, a sophisticated FE model is constructed by referring to an example model used in a previous study [17]. Using the FE model, the structural responses are analyzed by nonlinear time-history analyses. The elastic buckling, elephant’s foot buckling, and steel yielding of the tank wall are introduced as critical limit states of the tank for the fragility analysis.

2. Earthquake Ground Motions

In recent years, earthquake occurrences have increased in the Korean Peninsula. Figure 1 presents the number of earthquake occurrences from 2008 to 2017 in Korea. As shown in the figure, the number of earthquake occurrences increased significantly in 2016 and 2017, and the earthquake events (≥ local magnitude, 3.0) were also prevalent in this period.

Table 7 summarizes 52 earthquakes records (≥3.0) that occurred in southeastern Korea from 2008 to 2017, observed in the Ulsan Seismological Observatory. These records were obtained from the National Earthquake Comprehensive Information System (NECIS) database [19], which is operated by the Korea Meteorological Administration (KMA). Among these 52 earthquakes, the strongest one was the 12 September 2016 5.8 Gyeongju earthquake and the second strongest one was the 15 November 2017 5.4 Pohang earthquake, both of which triggered numerous aftershocks.

2.1. Characteristics of Earthquake Ground Motions

Characteristics of earthquake ground motions are highly dependent on various factors, including earthquake magnitudes, fault types, soil conditions, and source-to-site distances [20]. Particularly, it is well known that Korea is situated in the seismically stable region [21] and soil and rock conditions are different from other regions [22], resulting in unique earthquake ground motion characteristics in Korea. Hence, it is necessary to ascertain the characteristics of earthquake ground motions in the area under consideration prior to the assessment of a structure’s seismic fragility. The ratio of the peak ground acceleration (PGA) to the peak ground velocity (PGV), which is commonly called the ratio, is an informative measure for understanding the characteristics of earthquake ground motions [2325]. The ratios are typically classified as low, intermediate, and high. Previous studies [36] have developed seismic fragility curves of structures with different ratios and emphasized the importance of the characteristics of earthquake ground motions on the seismic fragility. This is related to the relationship of the earthquake magnitudes with the epicentral distances, frequency properties, and durations of ground motions. For example, earthquake ground motions included in the high category mostly are associated with relatively low- to moderate-magnitude earthquakes and have high frequency components and short durations, compared to those in the low category [2325].

In this study, the 52 earthquake accelerograms (Table 7) are analyzed. The gain factor is used to convert a binary expression (i.e., count) for the voltage intensity into acceleration units [26]. The gain factor varies with the types of accelerographs and data-loggers; the values used here are 2.15 × 106 count/m/s2 (before August 2016) and 1.71 × 106 count/m/s2 (after August 2016). This study applied the time-domain recursive filter [27] to the raw acceleration time series. The accelerograms are recorded in three directions: two horizontal components, east-west (EW) and north-south (NS), and one vertical component (Z). For each direction, the PGA and PGV values are estimated. The 52 earthquake ground motions are then categorized by the ratios as low ( < 0.8 /m/s), intermediate (0.8 /m/s ≤ ≤ 1.2 /m/s), and high (1.2 /m/s < ). The three ranges of ratios correspond to the three zonal combinations (i.e., < , = , and > ) in National Building Code of Canada (NBCC) (NRCC, 1985) [28]. and denote the acceleration-related and velocity-related seismic zones, respectively. These seismic zones are based on results from statistical analyses of past earthquakes in Canada [29]. Figure 2 shows the results of the classification. In the figure, the left graph presents the results in the arithmetic scale, and a graph of the common logarithm scale (base 10) is plotted on the right to show the distribution of the results more clearly. It is found in the figure that all of the ground motions are included in the high level, indicating that the earthquake ground motions that occurred in southeastern Korea have high-frequency and short-duration motion features, as explained above. Thus, it is important to derive seismic fragility curves based on real earthquake ground records which can present such ground characteristics.

2.2. Input Ground Motions for Seismic Fragility Analysis

Practical seismic design guidelines, such as ASCE 7-05 [30], Eurocode 8-1 [31], and FEMA 368 [32], recommend a proper number of input ground motions. The provisions of seismic design codes were summarized in a review paper by Beyer and Bommer [33], suggesting the use of at least seven records for the seismic analysis and design of a structure. Thus, for a seismic fragility analysis, this study uses seven accelerogram data out of the 52 earthquake accelerograms for nonlinear time-history analyses. Using seven real earthquake ground motions not only enables reflection of the dynamic properties (i.e., high ratio) of the region, but also considers the uncertainty of input ground motions.

Among the earthquake events presented in Table 7, the seven earthquakes with the largest magnitudes are selected for the nonlinear time-history analyses: 5.8, 5.4, 5.1, 5.0, 4.5, 4.3, and 4.0. Figure 13 shows the ground acceleration time-histories adopted in the analyses. Both the horizontal components (EW and NS) of each earthquake ground motion were simultaneously assigned to the structure.

3. FE Model of Steel Liquid Storage Tank

This study aims to derive the seismic fragility of a steel tank containing liquids, which is an important facility of industrial complexes, considering the hydrodynamic forces caused by oscillating liquids. Steel tanks can be classified to two groups, anchored and unanchored tanks, and it has been reported that their dynamic behaviors are quite different [35, 36]. This study examines the seismic fragility of an anchored steel liquid storage tank, which was relatively less well studied in previous studies but still an important facility of industrial complexes. In addition, two tank models, with and without ring stiffeners, are investigated for comparison purposes.

An FE model is constructed using a commercial FE software package, ABAQUS. The steel tank considered here is a typical steel tank and was introduced as a numerical example in a previous study [17]. It is assumed that the steel tank is made of ASTM A36 [18], whose material nonlinearity is presented in Figure 3, and the steel tank contains liquid (gasoline). Table 1 presents the material properties of the steel and liquid. The steel and liquid are modeled using shell elements (S4) and continuum elements (C3D8), respectively. To account for hydrodynamic effects in an oscillating tank, a dynamic viscosity (0.006 Pa·s) and linear (shock velocity) − (particle velocity) Hugoniot form of the Mie-Grüneisen equation of state with a velocity of sound of 1250 m/s in liquid are applied. Grüneisen’s gamma () and the linear Hugoniot slope coefficient () are assumed to be zero.

MaterialDensity (kg/m3)Young’s modulus (GPa)Poisson’s ratio

Liquid (Gasoline)680--

For comparison purposes, a ring-stiffened FE model is also constructed. In general, ring stiffeners are employed to improve resistance to buckling, and this study aimed to quantify the increase of the tank’s buckling resistance against seismic loads. The ring-stiffened model has the same shape and material properties as the unstiffened model, and two ring stiffeners are added at the top and middle of the tank wall. Both stiffeners have the same cross-sections and dimensions. The dimensions of the ring stiffeners are properly designed by referring to the API 650 design code [37] and Sun et al. [38]. Figure 4 presents the geometric details of the liquid storage tank and ring stiffener, and Figure 5 shows the FE models in ABAQUS with and without ring stiffeners.

To validate the constructed FE model, its seismic responses in terms of peak sloshing liquid height are compared with those of a previous study [17], in which seven representative earthquake ground motions were used. Table 2 provides the details of the seven earthquakes. Nonlinear time-history analyses are performed using the seven ground motions and the constructed ABAQUS FE model. Figure 6 shows a comparison with the analysis results with those of Di Carluccio et al. [17]. The maximum estimated difference is 0.50 m (i.e., 4.19%), which is within an allowable range.

CodeEarthquake eventDateCountryStation namePGA ()

000196xaMontenegro04/15/1979YugoslaviaPetrovac-Hotel Olivia0.45
006334xaSouth Iceland (aftershock)06/21/2000IcelandSolheimar0.42
000199yaMontenegro04/15/1979YugoslaviaBar-Skupstina Opstine0.36
000535yaErzincan03/13/1992TurkeyErzincan-Meteorologij Mundurlugu0.51
006263yaSouth Iceland06/17/2000IcelandKaldarholt0.51
006328yaSouth Iceland (aftershock)06/21/2000IcelandKaldarholt0.39
006334yaSouth Iceland (aftershock)06/21/2000IcelandSolheimar0.72

4. Seismic Fragility Assessment

In a fragility curve, the failure probability () generally defines the conditional probability of the seismic demand () and exceeding capacity () given an intensity measure (), as in (1) [39]. Thus, as stated by Phan and Paolacci and Bakalis et al. [10, 40], it is important to select in seismic fragility assessment. This study aims to derive the seismic fragility curves of liquid storage tanks using real earthquake ground motions that can reflect the characteristics. In previous studies [3, 4], the characteristics of collected ground motions were analyzed based on PGA/PGV () ratios, and PGA is selected as in this study.

As the first step to calculate a failure probability, nonlinear time-history analyses are performed with varying values. Based on the results, probabilistic seismic demand models (PSDMs) are constructed by linear regression of the seismic responses and corresponding values [41]. The PSDMs have been broadly utilized to derive the analytical seismic fragility curves of various structures [7, 10, 16, 39, 4246].

For the linear regression, it is assumed that the seismic demand () is expressed by the following approximation [47]:where and are the parameters that are determined from the linear regression. By taking the natural logarithms of both sides of (2), the following linear equation is obtained:Using (3) and the results of the nonlinear time-history analyses with varying IM values, the unknown parameters (k and p) can be determined by a linear regression analysis. The logarithmic dispersion of the demand given IM, , can be calculated by (Nielson and DesRoches [42] and Mangalathu et al. [43]):where is the peak demand value obtained from the nonlinear time-history analyses and is the total number of ground motions.

Assuming lognormal distributions for both the demand and capacity, the failure probability in (1) is specifically expressed as follows:where and are the median values of the demand and capacity, respectively; is the logarithmic dispersion of the capacity given IM; and is the standard normal cumulative distribution function.

4.1. Nonlinear Time-History Analyses

Using each of the input ground motions, nonlinear time-history analyses are performed with varying earthquake intensities. This parametric analysis method is known as an incremental dynamic analysis (IDA) [48]. Both horizontal components of each input ground motion are applied simultaneously. The input ground motions are scaled in such a way that the component with the higher PGA (primary) becomes the designated value, and the other component (secondary) is then also adjusted by the same scaling factor [49]. Table 3 summarizes the original and scaled PGA values for the primary and secondary components of the input ground motions.

No.MagnitudeComponentOriginal PGA ()Scaled PGA ()


As a result, the PSDMs are derived for each of the selected limit states (i.e., elastic buckling, elephant’s foot buckling, and steel yielding) using (3). Each damage state is analyzed in terms of radial displacement, compressive meridional stress, and von-Mises stress, respectively. Figure 7 shows the PSDMs of the unstiffened model using a log-log scale. Table 4 presents the regression coefficients of the PSDMs (i.e., and in (3)).

Model typeLimit state

Unstiffened tankElastic buckling0.02600.3335
Elephant’s foot buckling35.53970.1723
Steel yielding269.26610.2402
Ring-stiffened tankElastic buckling0.01550.2942
Elephant’s foot buckling38.49780.3810
Steel yielding268.21800.2307

As mentioned above, dynamic buckling is a well-known structural behavior of steel liquid storage tanks subjected to earthquake ground motions. Buratti and Tavano [9], Virella et al. [13], Djermane et al. [14], and Sobhan et al. [15] investigated the buckling stability of steel tanks with numerical models. The buckling phenomena in a liquid-oscillating steel tank are generally classified as elastic and elastoplastic buckling. The former includes diamond-shaped buckling and secondary buckling, while the latter is often called elephant’s foot buckling [9, 13]. This study deals with both elastic and elephant’s foot buckling phenomena. The former is typically found at the upper part of the steel tank, which has a larger deflection than the lower part, and the latter commonly exhibits a bulge of the near the bottom of the tank wall.

The impacts of the ring stiffeners and hydrodynamic forces owing to the oscillating liquid are investigated through nonlinear time-history analyses by employing an input ground motion of 4.0 and 0.4 . Figure 8 depicts the deformed shapes of the unstiffened and ring-stiffened models using the same deformation scale factor in ABAQUS. As shown in the figure, the deformation of the unstiffened model is observed to be much larger than that of the ring-stiffened model. More detailed results are quantitatively discussed with the radial displacement of the tank wall in the following section.

Figure 9 compares the von-Mises stress results of nonoscillating and oscillating liquid storage tanks without stiffeners. The objective of the comparative analysis is to understand the hydrodynamic effect. It is found in the figure that the oscillating liquid storage tank has a relatively higher level of stress compared to the nonoscillating tank, particularly at the lower part of the tank in which the maximum stress occurred. This phenomenon is due to the hydrodynamic forces of the impulsive component of the liquid. Hence, the steel yielding of the tank wall and the buckling phenomena are investigated in this study.

Figure 10 shows the stress distributions of the unstiffened and ring-stiffened tank models; the maximum von-Mises stress occurs at the lower part of the tank wall. It is also found that ring stiffeners have little contribution to reducing the stress, unlike deformation.

4.2. Seismic Fragility Curves

This study investigates the seismic fragility of liquid storage tanks for three limit states: elastic buckling, elephant’s foot buckling, and steel yielding. The threshold value for elastic buckling is computed based on previous literature [50]. According to the Korea Meteorological Administration archives [10], the earthquake on September 12, 2016 was the largest earthquake event and the PGA was approximately 0.43 PGA. In addition, the Next Generation Attenuation Relationships for Western US (NGA-West2) ground motion prediction equations (GMPEs) [51] yields median PGAs + one standard deviation up to 0.7 for a 5.4 earthquake when a source-to-site distance is 2 km (which is a feasible scenario for Korea). Based on these, the radial displacement is investigated up to 0.7 PGA in this numerical example.

In addition, the Budiansky-Roth criterion is widely used to obtain the dynamic buckling load for a tank [9, 13, 14]. The critical point at which the radial displacement suddenly increases as the earthquake intensity increases is regarded as the dynamic buckling load, which is identified by the pseudo-equilibrium path relating the IM (i.e., PGA in this study) and maximum radial displacement, as shown in Figure 11. The nodes on the upper part of the tank are monitored in this study, which is related to elastic buckling, and the critical points are estimated from the pseudo-equilibrium paths to be 0.0152 m and 0.0104 m for unstiffened and ring-stiffened tanks, respectively.

To determine the threshold value for elephant’s foot buckling, the buckling stress capacity () and compressive meridional stress demand () are analyzed and plotted in the PGA-stress plane, and then the intersection point is estimated [40]. The compressive meridional stress demand is analyzed from a nonlinear time-history analysis. The buckling stress capacity is computed using the empirical equation [52], which is also adopted in Eurocodes 8-4 [53]:where (=0.6/) is the ideal critical buckling stress, E the steel elastic modulus, the wall thickness, and the tank radius. is equal to ()/400 and is the steel yield strength. is the maximum interior pressure, which is the sum of the hydrostatic pressure and hydrodynamic pressure of impulsive component (). The latter pressure is estimated by where , , and are the nondimensional coordinates of radius, height, and angle, respectively. indicates the distribution of in the vertical direction, is the liquid density, and is the liquid height in the tank. is the ground acceleration time-history.

Failure probabilities are calculated by the standard normal cumulative distribution function using (5). Table 5 presents the statistical parameters of the capacity and demand: , , and . The median value of is estimated by the PSDM. The value of for steel yielding is obtained by referring to the probabilistic model code from the Joint Committee on Structural Safety [54].

Model typeLimit state

Unstiffened tankElastic buckling0.0151 (m)0.13680.2011
Elephant’s foot buckling36.5227 (MPa)0.09700.3672
Steel yielding247.6041 (MPa)0.06990.1128
Ring-stiffened tankElastic buckling0.0103 (m)0.11020.1550
Elephant’s foot buckling31.3111 (MPa)0.21920.4450
Steel yielding247.6041 (MPa)0.06990.1157

From the PSDMs and statistical parameters, seismic fragility curves are derived with the given limit states (Figure 12). As shown in the fragility curves, the steel tank is more vulnerable to elastic buckling than elephant’s foot buckling and steel yielding; the maximum differences in the failure probabilities of elastic buckling and steel yielding and of elastic buckling and steel yielding are 0.501 at 0.58 PGA and 0.687 at 0.38 PGA for the unstiffened model, respectively. This means that a large deflection of the upper part of the tank wall, due to the elastic buckling instability, is likely to occur prior to reaching the elephant’s foot buckling capacity and yield strength of the steel.

When comparing the failure probabilities of the ring-stiffened and unstiffened models at the same earthquake intensities in terms of elastic buckling, although the ring stiffener has little influence on the seismic fragility to steel yielding, it can reduce the failure probability owing to elastic buckling. To investigate this further, the failure probabilities between 0.1 and 0.4 , which is the range corresponding to past major earthquakes in southeast Korea, are provided in Table 6.

PGA ()

Unstiffened tank0.1780.2510.3230.3910.4540.5120.5630.610
Ring-stiffened tank0.0790.1290.1860.2460.3070.3660.4230.476

PGA ()

Unstiffened tank0.6510.6880.7210.7500.7750.7980.8180.836
Ring-stiffened tank0.5260.5710.6120.6500.6840.7150.7430.767

Record No.Date & TimeEpicenter
Latitude (°)Longitude (°)

13.512/25/2017 16:1936.11129.36
23.611/20/2017 06:0536.14129.36
33.511/19/2017 23:4536.12129.36
44.311/15/2017 16:4936.12129.36
53.511/15/2017 15:0936.09129.34
63.611/15/2017 14:3236.10129.36
75.411/15/2017 14:2936.12129.36
83.104/15/2017 11:3136.11129.36
93.303/31/2017 13:4635.78129.20
103.301/06/2017 05:3135.75129.17
113.312/14/2016 17:2035.76129.17
123.312/12/2016 17:5335.76129.18
133.310/10/2016 22:5935.75129.18
143.010/02/2016 20:5335.75129.19
153.109/28/2016 16:3435.76129.18
163.509/21/2016 11:5335.75129.18
174.509/19/2016 20:3335.74129.18
183.009/14/2016 00:4835.75129.18
193.009/13/2016 14:3135.76129.18
203.209/13/2016 08:2435.76129.17
213.109/13/2016 00:3735.78129.21
223.109/12/2016 23:5235.75129.18
233.009/12/2016 20:3835.78129.20
243.609/12/2016 20:3435.78129.19
255.809/12/2016 20:3235.76129.19
263.109/12/2016 20:1035.77129.18
273.109/12/2016 19:4835.77129.19
285.109/12/2016 19:4435.77129.19
295.007/05/2016 20:3335.51129.99
303.006/03/2016 04:5336.39127.92
313.001/06/2016 20:3936.01128.07
323.211/24/2015 11:2735.67129.82
333.809/25/2014 02:2635.07129.94
343.509/23/2014 15:2735.80129.41
353.507/03/2014 21:5735.66129.76
363.003/28/2014 10:4036.47127.93
373.610/11/2013 16:0636.46129.61
383.008/13/2013 00:3636.91129.44
393.108/12/2013 04:3335.66129.75
403.502/05/2013 21:2535.77127.97
413.102/04/2013 05:3936.55128.90
423.011/22/2012 03:0435.22127.97
433.405/30/2012 02:4836.57129.57
443.202/24/2012 09:0535.20129.93
453.011/02/2011 07:2035.54130.65
463.205/29/2011 10:2235.58128.75
473.203/28/2011 13:5035.97129.95
483.302/16/2010 18:5335.63129.95
493.006/28/2009 22:0735.78128.29
504.005/02/2009 07:5836.56128.71
513.001/13/2009 19:0736.47129.64
523.512/19/2008 17:5336.49129.68

Selected input ground motions for nonlinear time-history analyses.

The table shows that the ring-stiffener reduces the failure probability of the steel tank to elastic buckling by 0.069–0.147. In particular, the maximum difference is 0.147 at 0.18 PGA. From these results, it is confirmed that a ring-stiffener can significantly reduce the seismic failure risk of liquid-containing steel tanks. In contrast, in the case of elephant’s foot buckling, the ring stiffeners have no practical effect, especially for PGA levels larger than 0.34 PGA.

5. Conclusions

This study analytically derives seismic fragility curves of anchored liquid storage tanks. Because the characteristics of earthquake ground motions are important for influencing the seismic vulnerability of structures, real earthquake ground motions recorded in Ulsan Seismological Observatory are used in nonlinear time-history analyses. All components of the ground motions are included in high levels representing high-frequency ground motions, short durations, and small-to-medium magnitude earthquakes.

Hydrodynamic forces owing to an oscillating liquid in a tank are also a significant cause of structural damage to liquid storage tanks. Thus, this study develops two advanced FE models with and without ring stiffeners to consider hydrodynamic effect and analyze yielding stress acted on the tank walls. Furthermore, since the buckling phenomenon mainly occurs in thin-walled structures such as steel liquid storage tanks, the elastic buckling and elephant’s foot buckling phenomena are also investigated.

Seismic fragility curves with different limit states (i.e., elastic buckling, elephant’s foot buckling, and steel yielding) were successfully derived based on PSDMs. By using fragility curves, it is observed that the steel liquid storage tank is more vulnerable to elastic buckling than elephant’s foot buckling and steel yielding with maximum differences of failure probabilities (0.501 at 0.58 PGA and 0.687 at 0.38 PGA) for the unstiffened model. Furthermore, although there was little effect of the ring stiffener on reducing the steel yielding fragility, it is observed that the failure probabilities of the ring-stiffened tank to elastic buckling decreases up to 0.147 compared with the unstiffened one, which indicates that a ring-stiffener can significantly reduce the seismic failure risk of liquid-containing steel tanks. In contrast to the case of elastic buckling, the ring stiffeners had no practical effect on elephant’s foot buckling fragility since the failure probabilities of the ring-stiffened model were larger than those of the unstiffened model for PGA levels larger than 0.34 PGA.


Table 7 summarizes the occurrence of earthquake events in the southeastern Korean Peninsula. These records correspond to earthquake events of ≥ 3.0 in the last 10 years. There were 52 earthquake events observed at the Ulsan Seismological Observatory. The records were collected from the NECIS database. Figure 13 shows the selected input ground motions used in the time-history analyses.

Data Availability

The data used to support the findings of this study are available from the corresponding author upon request.

Conflicts of Interest

The authors declare that there are no conflicts of interest regarding the publication of this paper.


This work was supported by the 2019 Research Fund (1.190047.01) of UNIST (Ulsan National Institute of Science and Technology). This research was also supported by a grant (19SCIP-B146959-02) from Construction technology research program funded by Ministry of Land, Infrastructure and Transport of Korean government.


  1. H. Persson and A. Lönnermark, “Tank fires,” Tech. Rep. SP Report 2004:14, SP Swedish National Testing and Research Institute, Boras, Sweden, 2004. View at: Google Scholar
  2. “Korea Industrial Complex Corporation,” 2018. View at: Google Scholar
  3. O. Kwon and A. Elnashai, “The effect of material and ground motion uncertainty on the seismic vulnerability curves of RC structure,” Engineering Structures, vol. 28, no. 2, pp. 289–303, 2006. View at: Publisher Site | Google Scholar
  4. D. Moon, Y. Lee, and S. Lee, “Fragility analysis of space reinforced concrete frame structures with structural irregularity in plan,” Journal of Structural Engineering, vol. 144, no. 8, Article ID 04018096, 2018. View at: Publisher Site | Google Scholar
  5. Y. Lee and D. Moon, “A new methodology of the development of seismic fragility curves,” Smart Structures and Systems, vol. 14, no. 5, pp. 847–867, 2014. View at: Publisher Site | Google Scholar
  6. S. Lee, Y.-J. Lee, and H.-Y. Tak, “Fragility assessment of liquid storage tanks under seismic loads by integration of two software packages,” in Proceedings of the 12th International Conference on Structural Safety and Reliability (ICOSSAR2017), TU Wien, Vienna, Austria, August 2014. View at: Google Scholar
  7. X. Chu, J. M. Ricles, and S. N. Pakzad, “Seismic fragility analysis of the smithsonian institute museum support center,” Earthquake Spectra, vol. 33, no. 1, pp. 85–108, 2017. View at: Publisher Site | Google Scholar
  8. Y. Pang, X. Wu, G. Shen, and W. Yuan, “Seismic fragility analysis of cable-stayed bridges considering different sources of uncertainties,” Journal of Bridge Engineering, vol. 19, no. 4, Article ID 04013015, 2013. View at: Publisher Site | Google Scholar
  9. N. Buratti and M. Tavano, “Dynamic buckling and seismic fragility of anchored steel tanks by the added mass method,” Earthquake Engineering & Structural Dynamics, vol. 43, no. 1, pp. 1–21, 2014. View at: Publisher Site | Google Scholar
  10. H. N. Phan, F. Paolacci, and S. Alessandri, “Fragility analysis methods for steel storage tanks in seismic prone areas,” in Proceedings of the ASME 2016 Pressure Vessels and Piping Conference, American Society of Mechanical Engineers, British Columbia, Canada, July 2016. View at: Publisher Site | Google Scholar
  11. G. W. Housner, “The dynamic behavior of water tanks,” Bulletin of Seismological Society of America, vol. 53, no. 2, pp. 381–387, 1963. View at: Google Scholar
  12. A. S. Veletsos, “Seismic effects in flexible liquid storage tanks,” in Proceedings of the 5th world conference on earthquake engineering, pp. 630–639, Roma, Italy, June 1974. View at: Google Scholar
  13. J. Virella, L. Godoy, and L. Suárez, “Dynamic buckling of anchored steel tanks subjected to horizontal earthquake excitation,” Journal of Constructional Steel Research, vol. 62, no. 6, pp. 521–531, 2006. View at: Publisher Site | Google Scholar
  14. M. Djermane, D. Zaoui, B. Labbaci, and F. Hammadi, “Dynamic buckling of steel tanks under seismic excitation: numerical evaluation of code provisions,” Engineering Structures, vol. 70, pp. 181–196, 2014. View at: Publisher Site | Google Scholar
  15. M. Sobhan, F. Rofooei, and N. K. Attari, “Buckling behavior of the anchored steel tanks under horizontal and vertical ground motions using static pushover and incremental dynamic analyses,” Thin-Walled Structures, vol. 112, pp. 173–183, 2017. View at: Publisher Site | Google Scholar
  16. F. Paolacci, H. N. Phan, D. Corritore, S. Alessandri, O. S. Bursi, and M. S. Reza, “Seismic fragility analysis of steel storage tanks,” in Proceedings of the 5th ECCOMAS Thematic Conference on Computational Methods in Structural Dynamics and Earthquake Engineering, COMPDYN 2015, pp. 2054–2065, Crete, Greece, May 2015. View at: Google Scholar
  17. A. Di Carluccio, G. Fabbrocino, and G. Manfredi, “FEM seismic analysis of steel tanks for oil storage in industrial facilities,” in Proceedings of the 14th World Conference on Earthquake Engineering, Beijing, China, October 2008. View at: Google Scholar
  18. P. F. Adams, “High-strength steels for plastic design,” Engineering Journal, AISC, vol. 3, no. 4, p. 150, 1966. View at: Google Scholar
  19. “National Earthquake Comprehensive Information System,” 2017. View at: Google Scholar
  20. T. G. Cork, J. H. Kim, G. P. Mavroeidis, J. K. Kim, B. Halldorsson, and A. S. Papageorgiou, “Effects of tectonic regime and soil conditions on the pulse period of near-fault ground motions,” Soil Dynamics and Earthquake Engineering, vol. 80, pp. 102–118, 2016. View at: Publisher Site | Google Scholar
  21. C. S. Kim, “Overview of the techonic environment in Korea with reference to HLW disposal,” Tech. Rep. KR9700364, Korea Atomic Energy Research Institute, Daejeon, South Korea, 1997. View at: Google Scholar
  22. H. Jung, Y. Jang, and B. G. Jo, “Upper-crust shear-wave velocity of South Korea constrained by explosion and earthquake data,” Bulletin of the Seismological Society of America, vol. 101, no. 6, pp. 2819–2832, 2011. View at: Publisher Site | Google Scholar
  23. T. J. Zhu, A. C. Heidebrecht, and W. K. Tso, “Effect of peak ground acceleration to velocity ratio on ductility demand of inelastic systems,” Earthquake Engineering & Structural Dynamics, vol. 16, no. 1, pp. 63–79, 1988. View at: Publisher Site | Google Scholar
  24. W. Tso, T. Zhu, and A. Heidebrecht, “Engineering implication of ground motion A/V ratio,” Soil Dynamics and Earthquake Engineering, vol. 11, no. 3, pp. 133–144, 1992. View at: Publisher Site | Google Scholar
  25. T. Sawada, K. Hirao, H. Yamamoto, and O. Tsujihara, “Relation between maximum amplitude ratio (a/v, ad/v2) and spectral parameters of earthquake ground motion,” in Proceedings of the Tenth World Conference on Earthquake Engineering, pp. 617–622, Madrid, Spain, July 1992. View at: Google Scholar
  26. D. H. Sheen, I. S. Lim, J. H. Park, and H. C. Chi, “Earthquake magnitude determination using P phase for earthquake early warning,” Journal of the Geological Society of Korea, vol. 48, no. 1, pp. 101–111, 2012 (Korean). View at: Google Scholar
  27. H. Kanamori, P. Maechling, and E. Hauksson, “Continuous monitoring of ground-motion parameters,” Bulletin of the Seismological Society of America, vol. 89, no. 1, pp. 311–316, 1999. View at: Google Scholar
  28. National Research Council Canada (NRCC), National Building Code of Canada, Ottawa, ON, Canada, 1985.
  29. M. Paz, International Handbook of Earthquake Engineering: Codes, Programs, and Examples, Chapman & Hall, New York, NY, USA, 1994.
  30. American Society of Civil Engineers (ASCE), ASCE/SEI 7-05 Minimum Design Loads for Buildings and Other Structures, Reston, VA, USA, 2005.
  31. European Committee for Standardization (CEN), Eurocode 8: Design of Structures for Earthquake Resistance - Part 1: General Rules, Seismic Actions and Rules for Buildings, Brussels, Belgium, 2004.
  32. Building Seismic Safety Council (BSSC), NEHRP Recommended Provisions for Seismic Regulations for New Buildings and Other Structures, Part 1: Provisions (FEMA 368), Washington, DC, USA, 2000.
  33. K. Beyer and J. J. Bommer, “Selection and scaling of real accelerograms for bi-directional loading: a review of current practice and code provisions,” Journal of Earthquake Engineering, vol. 11, no. 1, pp. 13–45, 2007. View at: Publisher Site | Google Scholar
  34. European Strong-Motion Database. 2017.
  35. H. N. Phan, F. Paolacci, and S. Alessandri, “Enhanced seismic fragility analysis of unanchored steel storage tanks accounting for uncertain modeling parameters,” Journal of Pressure Vessel and Technology, vol. 141, no. 1, Article ID 010903, 2019. View at: Google Scholar
  36. H. N. Phan, F. Paolacci, and P. Mongabure, “Nonlinear finite element analysis of unanchored steel liquid storage tanks subjected to seismic loadings,” in Proceedings of the ASME PVP (American Society of Mechanical Engineers, Pressure Vessels and Piping Division) 2017 Conference, vol. 8, Hawaii, HI, USA, 2017. View at: Google Scholar
  37. American Petroleum Institute (API), API Standard 650 - Welded Tanks for Oil Storage, Washington, DC, USA, 2013.
  38. T. Sun, E. Azzuni, and S. Guzey, “Stability of open-topped storage tanks with top stiffener and one intermediate stiffener subject to wind loading,” Journal of Pressure Vessel Technology, vol. 140, no. 1, Article ID 011204, 2018. View at: Google Scholar
  39. J. E. Padgett and R. DesRoches, “Methodology for the development of analytical fragility curves for retrofitted bridges,” Earthquake Engineering & Structural Dynamics, vol. 37, no. 8, pp. 1157–1174, 2008. View at: Publisher Site | Google Scholar
  40. K. Bakalis, D. Vamvatsikos, and M. Fragiadakis, “Seismic risk assessment of liquid storage tanks via a nonlinear surrogate model,” Earthquake Engineering & Structural Dynamics, vol. 46, no. 15, pp. 2851–2868, 2017. View at: Publisher Site | Google Scholar
  41. K. R. Mackie and B. Stojadinović, “Comparison of incremental dynamic, cloud, and stripe methods for computing probabilistic seismic demand models,” in Proceedings of the Structures Congress 2005: Metropolis and Beyond, pp. 1–11, New York, NY, USA, April 2005. View at: Publisher Site | Google Scholar
  42. B. G. Nielson and R. DesRoches, “Analytical seismic fragility curves for typical bridges in the central and southeastern United States,” Earthquake Spectra, vol. 23, no. 3, pp. 615–633, 2007. View at: Publisher Site | Google Scholar
  43. S. Mangalathu, E. Choi, H. C. Park, and J.-S. Jeon, “Probabilistic seismic vulnerability assessment of tall horizontally curved concrete bridges in California,” Journal of Performance of Constructed Facilities, vol. 32, no. 6, Article ID 04018080, 2018. View at: Google Scholar
  44. J.-S. Jeon, J.-H. Park, and R. Desroches, “Seismic fragility of lightly reinforced concrete frames with masonry infills,” Earthquake Engineering & Structural Dynamics, vol. 44, no. 11, pp. 1783–1803, 2015. View at: Publisher Site | Google Scholar
  45. J. Jeon, R. DesRoches, T. Kim, and E. Choi, “Geometric parameters affecting seismic fragilities of curved multi-frame concrete box-girder bridges with integral abutments,” Engineering Structures, vol. 122, pp. 121–143, 2016. View at: Publisher Site | Google Scholar
  46. S. Mangalathu, J.-S. Jeon, J. E. Padgett, and R. DesRoches, “Performance-based grouping methods of bridge classes for regional seismic risk assessment: Application of ANOVA, ANCOVA, and non-parametric approaches,” Earthquake Engineering & Structural Dynamics, vol. 46, no. 14, pp. 2587–2602, 2017. View at: Publisher Site | Google Scholar
  47. C. A. Cornell, F. Jalayer, R. O. Hamburger, and D. A. Foutch, “Probabilistic basis for 2000 SAC federal emergency management agency steel moment frame guidelines,” Journal of Structural Engineering, vol. 128, no. 4, pp. 526–533, 2002. View at: Publisher Site | Google Scholar
  48. D. Vamvatsikos and C. Allin Cornell, “Incremental dynamic analysis,” Earthquake Engineering & Structural Dynamics, vol. 31, no. 3, pp. 491–514, 2002. View at: Publisher Site | Google Scholar
  49. R. P. Dhakal, J. B. Mander, and N. Mashiko, “Identification of critical ground motions for seismic performance assessment of structures,” Earthquake Engineering & Structural Dynamics, vol. 35, no. 8, pp. 989–1008, 2006. View at: Publisher Site | Google Scholar
  50. B. Budiansky and R. S. Roth, “Axisymmetric dynamic buckling of clamped shallow spherical shells,” in Collected Papers on Instability of Shell Structures, pp. 597–606, NASA-TN-1510, USA, 1962. View at: Google Scholar
  51. Y. Bozorgnia, N. A. Abrahamson, L. A. Atik et al., “NGA-West2 research project,” Earthquake Spectra, vol. 30, no. 3, pp. 973–987, 2014. View at: Publisher Site | Google Scholar
  52. J. M. Rotter, “Local collapse of axially compressed pressurized thin steel cylinders,” Journal of Structural Engineering, vol. 116, no. 7, pp. 1955–1970, 1990. View at: Publisher Site | Google Scholar
  53. European Committee for Standardization (CEN), Eurocode 8: Design of Structures for Earthquake Resistance - Part 4: Silos, Tanks And Pipelines, Brussels, Belgium, 2006.
  54. Joint Committee on Structural Safety. 2017.

Copyright © 2019 Sangmok Lee et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

More related articles

 PDF Download Citation Citation
 Download other formatsMore
 Order printed copiesOrder

Related articles

Article of the Year Award: Outstanding research contributions of 2020, as selected by our Chief Editors. Read the winning articles.