Finite element analysis of the seismic response of anchored and unanchored liquid storage tanks

The seismic response of liquid-filled cylindrical storage tanks has been investigated using finite element techniques implemented in the general purpose structural analysis computer code ANSYS. Both added mass concepts and displacement-based fluid finite elements were employed to allow for the effects of the liquid. Simplified response spectrum modal analyses of a tank making use of the axisymmetric harmonic displacement patterns of the principal modes of deformation were found to give accurate predictions of the tank behaviour with a rigidly anchored base. Time history analyses of three-dimensional finite element models of unanchored and flexibly anchored tanks, with gap conditions between the tank base and the supporting floor to allow lift-off of the base, indicated that stresses in the tank and resultant loads on the floor can be much greater than for a rigidly restrained tank. These results demonstrate the importance of carefully considering the restraint conditions when performing seismic design calculations on storage tanks.     

Seismic risk assessment of liquid storage tanks via a nonlinear surrogate model

A performance‐based earthquake engineering approach is developed for the seismic risk assessment of fixed‐roof atmospheric steel liquid storage tanks. The proposed method is based on a surrogate single‐mass model that consists of elastic beam‐column elements and nonlinear springs. Appropriate component and system‐level damage states are defined, following the identification of commonly observed modes of failure that may occur during an earthquake. Incremental dynamic analysis and simplified cloud are offered as potential approaches to derive the distribution of response parameters given the seismic intensity. A parametric investigation that engages the aforementioned analysis methods is conducted on 3 tanks of varying geometry, considering both anchored and unanchored support conditions. Special attention is paid to the elephant's foot buckling formation, by offering extensive information on its capacity and demand representation within the seismic risk assessment process. Seismic fragility curves are initially extracted for the component‐level damage states, to compare the effect of each analysis approach on the estimated performance. The subsequent generation of system‐level fragility curves reveals the issue of nonsequential damage states, whereby significant damage may abruptly appear without precursory lighter damage states.     

Developing efficient scalar and vector intensity measures for IDA capacity estimation by incorporating elastic spectral shape information

Scalar and vector intensity measures are developed for the efficient estimation of limit‐state capacities through incremental dynamic analysis (IDA) by exploiting the elastic spectral shape of individual records. IDA is a powerful analysis method that involves subjecting a structural model to several ground motion records, each scaled to multiple levels of intensity (measured by the intensity measure or IM), thus producing curves of structural response parameterized by the IM on top of which limit‐states can be defined and corresponding capacities can be calculated. When traditional IMs are used, such as the peak ground acceleration or the first‐mode spectral acceleration, the IM‐values of the capacities can display large record‐to‐record variability, forcing the use of many records to achieve reliable results. By using single optimal spectral values as well as vectors and scalar combinations of them on three multistorey buildings significant dispersion reductions are realized. Furthermore, IDA is extended to vector IMs, resulting in intricate fractile IDA surfaces. The results reveal the most influential spectral regions/periods for each limit‐state and building, illustrating the evolution of such periods as the seismic intensity and the structural response increase towards global collapse. The ordinates of the elastic spectrum and the spectral shape of each individual record are found to significantly influence the seismic performance and they are shown to provide promising candidates for highly efficient IMs. Copyright © 2005 John Wiley & Sons, Ltd.     

Wavelet‐based non‐stationary seismic rocking response of flexibly supported tanks

The non‐stationary rocking response of liquid storage tanks under seismic base excitations including soil interaction has been developed based on the wavelet domain random vibration theory. The ground motion has been characterized through statistical functionals of wavelet coefficients of the ground acceleration history. The tank–liquid–foundation system is modelled as a multi‐degree‐of‐freedom (MDOF) system with both lateral and rocking motions of vibration of the foundation. The impulsive and convective modes of vibration of the liquid in the tank have been considered. The wavelet domain coupled dynamic equations are formulated and then solved to get the expressions of instantaneous power spectral density function (PSDF) in terms of functionals of input wavelet coefficients. The moments of the instantaneous PSDF are used to obtain the stochastic responses of the tank in the form of coefficients of hydrodynamic pressure, base shear and overturning base moment for the largest expected peak responses. Parametric variations are carried out to study the effects of various governing parameters like height of liquid in the tank, height–radius ratio of the tank, ratio of total liquid mass to mass of foundation, and shear wave velocity in the soil medium, on the responses of the tank. Copyright © 2003 John Wiley & Sons, Ltd.     

Dynamic buckling and seismic fragility of anchored steel tanks by the added mass method

Buckling plays a fundamental role in the design of steel tanks because of the small thicknesses of the walls of this class of structures. The first part of the paper presents a review of this phenomenon for liquid‐containing circular cylindrical steel tanks that are fully anchored at the base, considering the different buckling modes and especially the secondary buckling occurring in the top part of the tank. A case study based on a cylindrical tank is then introduced in order to investigate various aspects of dynamic buckling. The finite element model of the case study tank is set‐up using the added mass method for fluid modelling. The influence of pre‐stress states caused by hydrostatic pressure and self‐weight on the natural periods of the structure is first studied and it is found that this influence is very small as far as the global behaviour of the tanks is considered, while it is important for local, shell‐type, vibration modes. In the following, the efficiency and sufficiency of different ground motion intensity measures is analysed by means of cloud analysis with a set of 40 recorded accelerograms. In particular, the peak ground displacement has been found being the most efficient and sufficient intensity measure so far as the maximum relative displacement of the tank walls is concerned. Finally, incremental nonlinear time‐history analyses are performed considering the case study structure under recorded earthquake ground motions in order to identify the critical buckling loads and to derive fragility curves for the buckling limit state. Copyright © 2013 John Wiley & Sons, Ltd.     

A vector‐valued ground motion intensity measure consisting of spectral acceleration and epsilon

The ‘strength’ of an earthquake ground motion is often quantified by an Intensity Measure (IM), such as peak ground acceleration or spectral acceleration at a given period. This IM is used to predict the response of a structure. In this paper an intensity measure consisting of two parameters, spectral acceleration and epsilon, is considered. The IM is termed a vector‐valued IM, as opposed to the single parameter, or scalar, IMs that are traditionally used. Epsilon (defined as a measure of the difference between the spectral acceleration of a record and the mean of a ground motion prediction equation at the given period) is found to have significant ability to predict structural response. It is shown that epsilon is an indicator of spectral shape, explaining why it is related to structural response. By incorporating this vector‐valued IM with a vector‐valued ground motion hazard, we can predict the mean annual frequency of exceeding a given value of maximum interstory drift ratio, or other such response measure. It is shown that neglecting the effect of epsilon when computing this drift hazard curve leads to conservative estimates of the response of the structure. These observations should perhaps affect record selection in the future. Copyright © 2005 John Wiley & Sons, Ltd.     

Seismic assessment of unanchored steel storage tanks by endurance time method

Liquid storage tanks are essential structures that are often located in residential and industrial areas; thus an assessment of their seismic performance is an important engineering issue. In this paper, the seismic response of unanchored steel liquid storage tanks is investigated using the endurance time (ET) dynamic analysis procedure and compared to responses obtained for anchored tanks under actual ground motions and intensifying ET records. In most cases, the results from ground motions are properly obtained with negligible differences using ET records. It is observed that uplifting of the tank base, which is closely related to the tank aspect ratio, has the greatest significance in the responses of the tank and can be predicted with reasonable accuracy by using currently available ET records.     

Non‐stationary seismic response of tanks with soil interaction by wavelets

A wavelet‐based random vibration theory has been developed for the non‐stationary seismic response of liquid storage tanks including soil interaction. The ground motion process has been characterized via estimates of statistical functionals of wavelet coefficients obtained from a single time history of ground accelerations. The tank–liquid–soil system has been modelled as a two‐degree‐of‐freedom (2‐DOF) system. The wavelet domain equations have been formulated and the wavelet coefficients of the required response state are obtained by solving two linear simultaneous algebraic equations. The explicit expression for the instantaneous power spectral density function (PSDF) in terms of the functionals of the input wavelet coefficients has been obtained. The moments of this PSDF are used to estimate the expected pseudo‐spectral acceleration (PSA) response of the tank. Parametric variations are carried out to study the effects of tank height, foundation natural frequency, shear wave velocity of soil and ratio of the mass of tank (including liquid) to the mass of foundation on the PSA responses of tanks. Copyright © 2001 John Wiley & Sons, Ltd.     

浮放储罐三维地震反应有限元分析

针对立式储罐,考虑液固耦合效应、地基与储罐结构的相互作用,采用有限元分析方法,对储罐在三维地震荷载作用下动反应进行了数值分析。分析结果表明：储罐三维地震加速度反应较一维地震加速度反应增加、提离高度明显放大、储罐轴向应力增加、基底剪力与弯矩增大。     

液压阻尼系统控制浮放立式储罐提离失稳

立式圆柱储罐在油田、石油化工等企业用途广泛，多用作存贮易燃、易爆介质，一旦遭遇震害，后果将十分严重。本文针对浮放立式储罐地震响应的提离失稳问题，采用液压阻尼系统（HDS）用以控制储罐的提离反应，建立了安装HDS储罐的力学分析模型并进行了数值计算，对安装HDS后储罐的提离反应进行了分析。仿真计算结果表明：HDS可以有效地减小提离应力，从而地震烈度可以降低1度进行储罐的提离设计。     

Modeling spatially correlated spectral accelerations at multiple periods using principal component analysis and geostatistics

Regional seismic risk assessments and quantification of portfolio losses often require simulation of spatially distributed ground motions at multiple intensity measures. For a given earthquake, distributed ground motions are characterized by spatial correlation and correlation between different intensity measures, known as cross‐correlation. This study proposes a new spatial cross‐correlation model for within‐event spectral acceleration residuals that uses a combination of principal component analysis (PCA) and geostatistics. Records from 45 earthquakes are used to investigate earthquake‐to‐earthquake trends in application of PCA to spectral acceleration residuals. Based on the findings, PCA is used to determine coefficients that linearly transform cross‐correlated residuals to independent principal components. Nested semivariogram models are then fit to empirical semivariograms to quantify the spatial correlation of principal components. The resultant PCA spatial cross‐correlation model is shown to be accurate and computationally efficient. A step‐by‐step procedure and an example are presented to illustrate the use of the predictive model for rapid simulation of spatially cross‐correlated spectral accelerations at multiple periods.     

Seismic vulnerability mitigation of liquefied gas tanks using concave sliding bearings

The aim of this paper is to evaluate the effectiveness of a concave sliding bearing system for the seismic protection of liquefied gas storage tanks through a seismic fragility analysis. An emblematic case study of elevated steel storage tanks, which collapsed during the 1999 ?zmit earthquake at Habas Pharmaceutics plant in Turkey, is studied. Firstly, a fragility analysis is conducted for the examined tank based on a lumped-mass stick model, where the nonlinear shear behaviour of support columns is taken into account by using a phenomenological model. Fragility curves in terms of an efficient intensity measure for different failure modes of structural components demonstrate the inevitable collapse of the tank mainly due to insufficient shear strength of the support columns. A seismic isolation system based on concave sliding bearings, which has been demonstrated a superior solution to seismically protect elevated tanks, is then designed and introduced into the numerical model, accounting for its non-linear behaviour. Finally, a vulnerability analysis for the isolated tank is performed, which proves a high effectiveness of the isolation system in reducing the probability of failure within an expected range of earthquake intensity levels.     

Seismic performance of wood‐frame houses in south‐western British Columbia

The seismic performance of conventional wood‐frame structures in south‐western British Columbia is analytically investigated through incremental dynamic analysis by utilizing available UBC‐SAWS models, which were calibrated based on experimental test results. To define an adequate target response spectrum that is consistent with information from national seismic hazard maps, record selection/scaling based on the conditional mean spectrum (CMS) is implemented. Furthermore, to reflect complex seismic hazard contributions from different earthquake sources (i.e. crustal events, interface events, and inslab events), we construct CMS for three earthquake types, and use them to select and scale an adequate set of ground motion records for the seismic performance evaluation. We focus on the impacts of adopting different record selection criteria and of using different shear‐wall types (Houses 1–4; in terms of seismic resistance, House 1>House 2>House 3>House 4) on the nonlinear structural response. The results indicate that the record selection procedures have significant influence on the probabilistic relationship between spectral acceleration at the fundamental vibration period and maximum inter‐story drift ratio, highlighting the importance of taking into account response spectral shapes in selecting and scaling ground motion records. Subjected to ground motions corresponding to the return period of 2500 years, House 1 is expected to experience very limited extent of damage; Houses 2 and 3 may be disturbed by minor damage; whereas House 4 may suffer from major damage occasionally. Finally, we develop statistical models of the maximum inter‐story drift ratio conditioned on a seismic intensity level for wood‐frame houses, which is useful for seismic vulnerability assessment. Copyright © 2010 John Wiley & Sons, Ltd.     

Response of tanks to vertical seismic excitations

A method for analyzing the earthquake response of elastic, cylindrical liquid storage tanks under vertical excitations is presented. The method is based on superposition of the free axisymmetrical vibrational modes obtained numerically by the finite element method. The validity of these modes has been checked analytically and the formulation of the load vector has been confirmed by a static analysis. Two forms of ground excitations have been used: step functions and recorded seismic components. The radial and axial displacements are computed and the corresponding stresses are presented. Both fixed and partly fixed tanks are considered to evaluate the effect of base fixation on tank behaviour. Finally, tank response under the simultaneous action of both vertical and lateral excitations is calculated to evaluate the relative importance of the vertical component of ground acceleration on the overall seismic behaviour of liquid storage tanks.     

Intensity measures for probabilistic assessment of non‐structural components acceleration demand

A fundamental issue in the framework of seismic probabilistic risk analysis is the choice of ground motion intensity measures (IMs). Based on the floor response spectrum method, the present contribution focuses on the ability of IMs to predict non‐structural components (NSCs) horizontal acceleration demand. A large panel of IMs is examined and a new IM, namely equipment relative average spectral acceleration (E‐ASA_R), is proposed for the purpose of NSCs acceleration demand prediction. The IMs efficiency and sufficiency comparisons are based on (i) the use of a large dataset of recorded earthquake ground motions; (ii) numerical analyses performed on three‐dimensional numerical models, representing actual structural wall and frame buildings; and (iii) systematic statistical analysis of the results. From the comparative study, the herein introduced E‐ASA_R shows high efficiency with respect to the estimation of maximum floor response spectra ordinates. Such efficiency is particularly remarkable in the case of structural wall buildings. Besides, the sufficiency and the simple formulation allowing the use of existing ground motion prediction models make the E‐ASA_R a promising IMs for seismic probabilistic risk assessment. Copyright © 2015 John Wiley & Sons, Ltd.     

Sliding fragility of block‐type non‐structural components. Part 2: Restrained components

This paper focuses on seismic vulnerability assessment of restrained block‐type non‐structural components under sliding response on the basis of seismic inputs specified by current seismic codes. The general representation of restrained equipment considered in this study consists of a rigid block restrained by four post‐tensioned, symmetrically arranged cables. Two sliding‐related failure modes are considered: restraint breakage and excessive absolute acceleration. Fragility analysis is proposed as an appropriate tool to evaluate these failure modes. Sample fragility curves developed through Monte‐Carlo simulations show that the restraint breakage limit state is sensitive to the parameters of the equation of motion. For instance, fragility estimates obtained without taking into account vertical base accelerations can be significantly unconservative for relatively large values of the coefficient of friction. In contrast, the excessive absolute acceleration limit state exhibits little sensitivity to the parameters of the equation of motion. Peak absolute acceleration response is almost always equal to or greater than the horizontal peak base acceleration. Representative results suggest that reasonable response estimates for blocks located at stories other than the ground in multistorey buildings can in general be obtained by simply scaling the ground acceleration to the peak acceleration at the corresponding storey. Copyright © 2002 John Wiley & Sons, Ltd.     

Constant‐ductility response spectra for sequential earthquake and tsunami loading

Simplified approaches for examining structural system response under sequential earthquake and tsunami loading are helpful for understanding response trends. To aid understanding, nonlinear (constant‐ductility) response spectra are developed for elastoplastic single degree of freedom systems subjected to seismic loads followed by hydrodynamic tsunami loads. The forcing function is composed of long‐duration earthquake motion concatenated with a range of tsunami hydrodynamic forces that are proportional to the pseudo‐spectral acceleration produced by the earthquake motion. The constant‐ductility spectra are thus constructed for scenarios where the loading imposed by one hazard is not dominant over the other. The spectra and basic intensity measures indicate that the amplification of response for sequential earthquake and tsunami loading over the earthquake only case is most significant for systems with long natural periods and high‐ductility capacity under seismic loading. Copyright © 2017 John Wiley & Sons, Ltd.     

Evaluation of the earthquake performance of anchored wine tanks during the San Juan,Argentina, 1977 earthquakexs

This paper details the damage sustained by a number of wine tanks during the San Juan, Argentina, 1977 earthquake. The available ground motion information is used to compare the actual tank performance with that predicted. Two of the tanks examined were located quite near the two accelerographs that recorded the ground motion in the city of San Juan for the main shock. An empirical attenuation formula which describes the intensity of the shaking was derived for this earthquake and used in a quantitative evaluation of the performance of a number of other tanks. By correlating and comparing the observed performance of the anchors and what was predicted we see the predictions are upheld by the actual damage. When the anchoring system failed during the earthquake the tanks were rendered effectively unanchored. Consequently, the state of stress that results from their base uplift was then examined in conjunction with the tank wall stability approach proposed by the author, and seen to give results consistent with the observed tank wall performance.     

The Umbria–Marche case: some suggestions for the Italian seismic zonation

Probabilistic seismic hazard maps, in terms of spectral acceleration and uniform hazard response spectra at given sites, considering local soil conditions, represent a much more complete estimate of the seismic hazard than the traditional maps in terms of peak ground acceleration or macroseismic intensity. This is particularly true when the requests of urban planners and engineers have to be met. The present analysis shows how some hazard parameters, such as the effective peak acceleration and the spectrum intensity, can well synthesise the overall information available from traditional probabilistic studies, but also suggests that soil condition is a first-order ingredient for effective seismic hazard mapping at national level. Three Italian towns, damaged by the 1997 Umbria–Marche earthquake sequence, are considered as example to demonstrate that: (1) soil condition dependent uniform hazard spectra well approximate actual spectra recorded during some events of the seismic sequence; (2) for these localities, the design spectrum of the present Italian seismic code does not seem adequate.

These considerations have induced the Italian scientific community to propose an updating of the national seismic zonation on the basis of several hazard parameters, that are described in this paper.     

Ductility‐dependent intensity measure that accounts for ground‐motion spectral shape and duration

The calculated nonlinear structural responses of a building can vary greatly, even if recorded ground motions are scaled to the same spectral acceleration at a building's fundamental period. To reduce the variation in structural response at a particular ground‐motion intensity, this paper proposes an intensity measure (IM_comb) that accounts for the combined effects of spectral acceleration, ground‐motion duration, and response spectrum shape. The intensity measure includes a new measure of spectral shape that integrates the spectrum over a period range that depends on the structure's ductility. The new IM is efficient, sufficient, scalable, transparent, and versatile. These features make it suitable for evaluating the intensities of measured and simulated ground motions. The efficiency and sufficiency of the new IM is demonstrated for the following: (i) elastic‐perfectly plastic single‐degree‐of‐freedom (SDOF) oscillators with a variety of ductility demands and periods; (ii) ductile and brittle deteriorating SDOF systems with a variety of periods; and (iii) collapse analysis for 30 previously designed frames. The efficiency is attributable to the inclusion of duration and to the ductility dependence of the spectral shape measure. For each of these systems, the transparency of the intensity measure made it possible to identify the sensitivity of structural response to the various characteristics of the ground motion. Spectral shape affected all structures, but in particular, ductile structures. Duration only affected structures with cyclic deterioration. Copyright © 2015 John Wiley & Sons, Ltd.