Research Article  Open Access
Seismic Fragility Analysis of Steel Liquid Storage Tanks Using Earthquake Ground Motions Recorded in Korea
Abstract
Liquidcontaining 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 ringstiffened 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 [3–6] 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 [3–10] have utilized real ground motion records from worldwide wellknown 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 liquidfilled tank owing to the liquidtank interaction, and many researchers have explored the buckling of steel tanks subjected to seismic loads with numerical models [9, 13–15]. 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 addedmass 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 liquidtank 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 wellknown 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 timehistory 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 sourcetosite 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 [23–25]. The ratios are typically classified as low, intermediate, and high. Previous studies [3–6] 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 moderatemagnitude earthquakes and have high frequency components and short durations, compared to those in the low category [23–25].
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 dataloggers; the values used here are 2.15 × 10^{6} count/m/s^{2} (before August 2016) and 1.71 × 10^{6} count/m/s^{2} (after August 2016). This study applied the timedomain recursive filter [27] to the raw acceleration time series. The accelerograms are recorded in three directions: two horizontal components, eastwest (EW) and northsouth (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 accelerationrelated and velocityrelated 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 highfrequency and shortduration 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.
(a)
(b)
2.2. Input Ground Motions for Seismic Fragility Analysis
Practical seismic design guidelines, such as ASCE 705 [30], Eurocode 81 [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 timehistory 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 timehistory analyses: 5.8, 5.4, 5.1, 5.0, 4.5, 4.3, and 4.0. Figure 13 shows the ground acceleration timehistories 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 MieGrü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.

For comparison purposes, a ringstiffened 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 ringstiffened 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 crosssections 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.
(a)
(b)
(a)
(b)
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 timehistory 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.

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 timehistory 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, 42–46].
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 timehistory 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 timehistory 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 TimeHistory Analyses
Using each of the input ground motions, nonlinear timehistory 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.

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 vonMises stress, respectively. Figure 7 shows the PSDMs of the unstiffened model using a loglog scale. Table 4 presents the regression coefficients of the PSDMs (i.e., and in (3)).

(a)
(b)
(c)
As mentioned above, dynamic buckling is a wellknown 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 liquidoscillating steel tank are generally classified as elastic and elastoplastic buckling. The former includes diamondshaped 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 timehistory analyses by employing an input ground motion of 4.0 and 0.4 . Figure 8 depicts the deformed shapes of the unstiffened and ringstiffened 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 ringstiffened model. More detailed results are quantitatively discussed with the radial displacement of the tank wall in the following section.
(a)
(b)
Figure 9 compares the vonMises 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.
(a)
(b)
Figure 10 shows the stress distributions of the unstiffened and ringstiffened tank models; the maximum vonMises 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.
(a)
(b)
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 (NGAWest2) ground motion prediction equations (GMPEs) [51] yields median PGAs + one standard deviation up to 0.7 for a 5.4 earthquake when a sourcetosite 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 BudianskyRoth 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 pseudoequilibrium 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 pseudoequilibrium paths to be 0.0152 m and 0.0104 m for unstiffened and ringstiffened tanks, respectively.
(a)
(b)
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 PGAstress plane, and then the intersection point is estimated [40]. The compressive meridional stress demand is analyzed from a nonlinear timehistory analysis. The buckling stress capacity is computed using the empirical equation [52], which is also adopted in Eurocodes 84 [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 timehistory.
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].

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.
(a)
(b)
(c)
(d)
(e)
(f)
(g)
When comparing the failure probabilities of the ringstiffened 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.

 
Selected input ground motions for nonlinear timehistory analyses. 
The table shows that the ringstiffener 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 ringstiffener can significantly reduce the seismic failure risk of liquidcontaining 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 timehistory analyses. All components of the ground motions are included in high levels representing highfrequency ground motions, short durations, and smalltomedium 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 thinwalled 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 ringstiffened tank to elastic buckling decreases up to 0.147 compared with the unstiffened one, which indicates that a ringstiffener can significantly reduce the seismic failure risk of liquidcontaining 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 ringstiffened model were larger than those of the unstiffened model for PGA levels larger than 0.34 PGA.
Appendix
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 timehistory 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.
Acknowledgments
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 (19SCIPB14695902) from Construction technology research program funded by Ministry of Land, Infrastructure and Transport of Korean government.
References
 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
 “Korea Industrial Complex Corporation,” 2018. http://www.kicox.or.kr/. View at: Google Scholar
 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
 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
 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
 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
 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
 Y. Pang, X. Wu, G. Shen, and W. Yuan, “Seismic fragility analysis of cablestayed bridges considering different sources of uncertainties,” Journal of Bridge Engineering, vol. 19, no. 4, Article ID 04013015, 2013. View at: Publisher Site  Google Scholar
 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
 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
 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
 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
 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
 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
 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,” ThinWalled Structures, vol. 112, pp. 173–183, 2017. View at: Publisher Site  Google Scholar
 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
 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
 P. F. Adams, “Highstrength steels for plastic design,” Engineering Journal, AISC, vol. 3, no. 4, p. 150, 1966. View at: Google Scholar
 “National Earthquake Comprehensive Information System,” 2017. http://necis.kma.go.kr/. View at: Google Scholar
 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 nearfault ground motions,” Soil Dynamics and Earthquake Engineering, vol. 80, pp. 102–118, 2016. View at: Publisher Site  Google Scholar
 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
 H. Jung, Y. Jang, and B. G. Jo, “Uppercrust shearwave 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
 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
 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
 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
 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
 H. Kanamori, P. Maechling, and E. Hauksson, “Continuous monitoring of groundmotion parameters,” Bulletin of the Seismological Society of America, vol. 89, no. 1, pp. 311–316, 1999. View at: Google Scholar
 National Research Council Canada (NRCC), National Building Code of Canada, Ottawa, ON, Canada, 1985.
 M. Paz, International Handbook of Earthquake Engineering: Codes, Programs, and Examples, Chapman & Hall, New York, NY, USA, 1994.
 American Society of Civil Engineers (ASCE), ASCE/SEI 705 Minimum Design Loads for Buildings and Other Structures, Reston, VA, USA, 2005.
 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.
 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.
 K. Beyer and J. J. Bommer, “Selection and scaling of real accelerograms for bidirectional 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
 European StrongMotion Database. 2017. http://isesd.hi.is/.
 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
 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
 American Petroleum Institute (API), API Standard 650  Welded Tanks for Oil Storage, Washington, DC, USA, 2013.
 T. Sun, E. Azzuni, and S. Guzey, “Stability of opentopped 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
 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
 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
 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
 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
 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
 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
 J. Jeon, R. DesRoches, T. Kim, and E. Choi, “Geometric parameters affecting seismic fragilities of curved multiframe concrete boxgirder bridges with integral abutments,” Engineering Structures, vol. 122, pp. 121–143, 2016. View at: Publisher Site  Google Scholar
 S. Mangalathu, J.S. Jeon, J. E. Padgett, and R. DesRoches, “Performancebased grouping methods of bridge classes for regional seismic risk assessment: Application of ANOVA, ANCOVA, and nonparametric approaches,” Earthquake Engineering & Structural Dynamics, vol. 46, no. 14, pp. 2587–2602, 2017. View at: Publisher Site  Google Scholar
 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
 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
 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
 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, NASATN1510, USA, 1962. View at: Google Scholar
 Y. Bozorgnia, N. A. Abrahamson, L. A. Atik et al., “NGAWest2 research project,” Earthquake Spectra, vol. 30, no. 3, pp. 973–987, 2014. View at: Publisher Site  Google Scholar
 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
 European Committee for Standardization (CEN), Eurocode 8: Design of Structures for Earthquake Resistance  Part 4: Silos, Tanks And Pipelines, Brussels, Belgium, 2006.
 Joint Committee on Structural Safety. 2017. http://www.jcss.byg.dtu.dk/.
Copyright
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.