A practical framework for performance-based reliability analysis of subway stations based on a faultestructure combined system

It is necessary to pay particular attention to the uncertainties that exist in an engineering problem to reduce the risk of seismic damage of infrastructures against natural hazards.Moreover,certain structural performance levels should be satisfied during strong earthquakes.However,these performance levels have been only well described for aboveground structures.This study investigates the main uncertainties involved in the performance-based seismic analysis of a multi-story subway station.More than 100 pulse-like and no pulse-like ground motions have been selected.In this regard,an effective framework is presented,based on a set of nonlinear static and dynamic analyses performed by OpenSees code.The probabilistic seismic demand models for computing the free-field shear strain of soil and racking ratio of structure are proposed.These models result in less variability compared with existing relations,and make it possible to evaluate a wider range of uncertainties through reliability analysis in Rtx software using the Monte Carlo sampling method.This work is performed for three different structural performance levels(denoted as PL1ePL3).It is demonstrated that the error terms related to the magnitude and location of earthquake excitations and also the corresponding attenuation relationships have been the most important parameters.Therefore,using a faultestructure model would be inevitable for the reliability analysis of subway stations.It is found that the higher performance level(i.e.PL3)has more sensitivity to random variables than the others.In this condition,the pulse-like ground motions have a major contribution to the vulnerability of subway stations.

Seismic response comparison and sensitivity analysis of pile foundation in liquefiable and non-liquefiable soils

Case history investigations have shown that pile foundations are more critically damaged in liquefiable soils than non-liquefiable soils.This study examines the differences in seismic response of pile foundations in liquefiable and non-liquefiable soils and their sensitivity to numerical model parameters.A two-dimensional finite element(FE)model is developed to simulate the experiment of a single pile foundation centrifuge in liquefiable soil subjected to earthquake motions and is validated against real-world test results.The differences in soil-pile seismic response of liquefiable and non-liquefiable soils are explored.Specifically,the first-order second-moment method(FOSM)is used for sensitivity analysis of the seismic response.The results show significant differences in seismic response for a soil-pile system between liquefiable and non-liquefiable soil.The seismic responses are found to be significantly larger in liquefiable soil than in non-liquefiable soil.Moreover,the pile bending moment was mainly affected by the kinematic effect in liquefiable soil,while the inertial effect was more significant in non-liquefiable soil.The controlling parameters of seismic response were PGA,soil density,and friction angle in liquefiable soil,while the pile bending moment was mainly controlled by PGA,the friction angle of soil,and shear modulus of loose sand in non-liquefiable soil.

Fractal analysis of major faults and fractal dimension of lineaments in the Indo-Gangetic Plain on a regional scale

The Indo-Gangetic Plain(IGP)is one of the most seismically vulnerable areas due to its proximity to the Himalayas.Geographic information system(GIS)-based seismic characterization of the IGP was performed based on the degree of deformation and fractal dimension.The zone between the Main Boundary Thrust(MBT)and the Main Central Thrust(MCT)in the Himalayan Mountain Range(HMR)experienced large variations in earthquake magnitude,which were identified by Number-Size(NS)fractal modeling.The central IGP zone experienced only moderate to low mainshock levels.Fractal analysis of earthquake epicenters reveals a large scattering of earthquake epicenters in the HMR and central IGP zones.Similarly,the fault fractal analysis identifies the HMR,central IGP,and south-western IGP zones as having more faults.Overall,the seismicity of the study region is strong in the central IGP,south-western IGP,and HMR zones,moderate in the western and southern IGP,and low in the northern,eastern,and south-eastern IGP zones.

Seismic fragility analysis of clay-pile-pier systems considering the optimization of ground motion intensity measures

The performance of clay-pile-pier system under earthquake shaking was comprehensively examined via three-dimensional finite element analyses,in which the complex stress-strain relationships of a clay and piled pier system were depicted by a hyperbolic-hysteretic and an equivalent elastoplastic model,respectively.One hundred twenty ground motions with varying peak accelerations were considered,along with the variations in bridge superstructure mass and pile flexural rigidity.Comprehensive comparison studies suggested that peak pile-cap acceleration and peak pile-cap velocity are the optimal ground motion intensity measures for seismic responses of the pier and the pile,respectively.Furthermore,based on two optimal ground motion intensity measures and using curvature ductility to quantify different damage states,seismic fragility analyses were performed.The pier generally had no evident damage except when the bridge girder mass was equal to 960 t,which seemed to be comparatively insensitive to the varying pile flexural rigidity.In comparison,the pile was found to be more vulnerable to seismic damage and its failure probabilities tended to clearly reduce with the increment of pile flexural rigidity,while the influence of the bridge girder mass was relatively minor.

Comparison of seismic fragility analysis methods for arch dams

This study focuses on the seismic fragility analysis of arch dams.The multiple stripe analysis(MSA),cloud analysis(CLA),and incremental dynamic analysis(IDA)methods are compared.A comprehensive dam-reservoir-foundation rock system,which considers the opening of contraction joints,the nonlinearity of dam concrete and foundation rock,the radiation damping effect of semi-unbounded foundation,and the compressibility of reservoir water,is used as a numerical example.225,80,and 15 earthquake records are selected for MSA,CLA,and IDA,respectively.The results show that MSA provides satisfactory fragility analysis,while both CLA and IDA have assumptions that may lead to deviations.Therefore,MSA is the most reliable method among the three methods and is recommended for the fragility analysis of arch dams.It is also shown that the choice of demand level affects the reliability of fragility curves and the effect of the material uncertainty on the fragility of the dam is not significant.

Generalized response displacement methods for seismic analysis of underground structures with complex cross section

The response displacement method(RDM)is recommended for the seismic analysis of underground structures in the transverse direction for many codes,including bases for design of structures-seismic actions for designing geotechnical works(ISO 23469)and code for seismic design of urban rail transit structures(GB 50909-2014).However,there are some obvious limitations in the application of RDM.Springs and the shear stress of the soil could be approximately evaluated for the structures having a simple cross section,such as rectangular and circular structures.It is necessary to propose simplified seismic analysis methods for structures with complex cross sections.This paper refers to the idea of RDM and proposes three generalized response displacement methods(GRDM).In GRDM1,a part of the soil surrounding a structure is selected to generate a generalized underground structure with a rectangular cross section,and the same analysis model as RDM is applied to analyze the responses of the structure.In GRDM2,a hollow soil model without a generalized structure is used to compute the equivalent load caused by the relative displacement of the soil,and the soil-structure interaction model is applied to calculate the responses of the structure.In GRDM3,a continuous soil model is applied to compute the equivalent load caused by the relative displacement and shear stress of the soil,and the soil-structure interaction model is applied to analyze the responses of the structure,which is the same as the model used in GRDM2.The time-history analysis method(THAM)is used to evaluate the accuracy of the proposed simplified methods.Results show that the error of GRDM1 is about 20%,while the error is only 5%for GRDM2 and GRDM3.Among the three proposed methods,GRDM3 has obvious advantages regarding calculation efficiency and accuracy.Therefore,it is recommended to use GRDM3 for the seismic response analysis of underground structures that have conventional simple or complex cross sections.

Lineament analysis as a seismic precursor:The El Mayor Cucapah earthquake of April 4,2010(MW7.2),Baja California,Mexico

An earthquake called the MW7.2“El Mayor Cucapah”earthquake on April 4,2010 has been analyzed for seismic precursor.The changes in the lineament system concerning its pattern and time intervals were analyzed during the earthquake preparation period and occurrence using the automated lineament detection method.The Landsat 5 TM images were processed using LESSA and ADALGEO software obtaining similar results.The statistical analysis revealed the stress accumulation due to plate interaction during earthquake formation.The study shows that the number of extracted lineaments changes rapidly about 23 months before the earthquake,and the systems return to the initial stage after 23 months.Most lineaments coincide with the extension of the San Andreas Fault as NW direction is the dominant trend.Thus,it can be concluded that the featural changes within the Rose diagram corresponding to the different strokes direction along with oriented elongation lines as disclosed in the present study using satellite images could be identified as a mid-term and/or short-term precursors of the earthquake.However,even though the dynamism of the El Mayor Cucapah earthquake is found in the extracted lineaments,it is possible to isolate more significant earthquakes even if new ones appear near the zone.Moreover,using two algorithms for lineament detection allows for the tectonics to corroborate the obtained lineaments and dynamism.

Detailed seismic risk analysis of electrical substation equipment using a reliability based approach

This paper proposes a risk analysis framework for substation structures based on reliability methods.Even though several risk assessment approaches have been developed for buildings,detailed risk analysis procedures for infrastructure components have been lacking in prior studies.The proposed framework is showcased by its application to a system of interconnected structures at a power substation in Tehran.Finite element models of structures are developed and validated in accordance with previous experiments.The uncertainties in the material,mass,and geometric properties of structures are described by random variables that are input to the finite element model.An artificial ground motion model is employed to comprehensively consider uncertainty in ground motion.Monte Carlo sampling is subsequently conducted on the library of probabilistic models.The analysis resulted in the loss distribution in the life cycle of structures.Additionally,the loss associated with six earthquake scenarios having specific magnitudes and return periods is computed.The application provides insight into the most vulnerable equipment in the considered system.Furthermore,introduced risk measures can guide stakeholders to make risk-based decisions to optimize design or prioritize a retrofit of infrastructure components under conditions of uncertainty.

A Seismic Facies Analysis to Determine the Relative Age and History of the Al Idrissi Mud Volcano from Offshore Larache Located in the NW Moroccan Atlantic Margin

Formed on top of the Gulf of Cadiz, the Al Idrissi mud volcano is the shallowest and largest mud volcano in the El Arraiche mud volcano field of the northwestern Moroccan margin. The development and morphology of mud volcanoes from the El Arraiche mud volcanoes group have been studied at a large scale. However, the time interval related to their formation period still needs to be better understood. In this regard, we interpreted and analyzed the seismic facies from the 2D reflection data of the GEOMARGEN-1 campaign, which took place in 2011. The aim was to identify the seismic sequences and draw the Al Idrissi mud volcano system to determine the formation period of the Al Idriss mud volcano. And as a result, the Al Idrissi mud volcano system is made of both buried and superficial bicone and was identified along with the Upper Tortonian to Messinian-Upper Pliocene facies. As the initial mud volcano extrusive edifice, the buried bicone was formed in the Late-Messinian to Early-Pliocene period. However, the superficial bicone, as the final extrusive edifice, was included in the Late Pliocene. In this case, the timing interval between the buried and superficial bicone is equivalent to the Late-Messinian to Upper-Pliocene period. Therefore, the latter corresponds to the Al Idrissi mud volcano formation period.

Performance-based passive control analysis of adjacent structures:Optimization of Maxwell dampers

The performance-based passive control analysis of the Maxwell dampers between one 10-story and one 6-story adjacent RC frames is conducted in this work.Not only the optimal parameters but also the optimal arrangements of the Maxwell dampers are proposed based on the optimal target of making the total exceeding probability of the adjacent structures to be minimal.The applicability of the analytical expressions of the Maxwell damper damping parameters under different seismic performance targets are firstly examined and then the preferable damping parameters of the Maxwell dampers are proposed through the extensive parametric studies.Furthermore,the optimal arranging positions and optimal arranging numbers of the Maxwell dampers between the adjacent buildings are derived based on a large number of seismic fragility analyses,as well.The general arranging laws of the Maxwell dampers between the adjacent buildings are generated based on the discussion of the theoretical method through the simplified plane model.The optimal parameters and optimal arrangement of the Maxwell dampers presented make both the adjacent structures have preferable controlled effects under each seismic performance target which can satisfy the requirements of multi-performance seismic resistance of the modern seismic codes.

Early Cretaceous Thrust and Nappe Tectonics in North Qilian Shan,Northern Tibetan Plateau:Evidence from Field Mapping,Geochronology,and Deep Structural Analysis

The North Qilian Shan fold and thrust belt,located at the northern Tibetan Plateau and southern margin of the Hexi Corridor,is a key tectonic unit to decode the formation and expansion of the plateau.Previous studies emphasize the Cenozoic deformation due to the far-field response to the Indo-Asian collision,but the Mesozoic deformations are poorly constrained in this area.We conducted detailed field mapping,structural analysis,geochronology,and structural interpretation of deep seismic reflectional profiling and magnetotelluric(MT)sounding,to address the superposed results of the Mesozoic and Cenozoic deformation.The results recognized the North Qilian thrust and nappe system(NQTS),the root and the frontal belt are the North Qilian thrust(NQT),and the Yumu Shan klippe(YK),respectively.The middle belt is located between the NQT and the YK.Monzonitic granite zircon U-Pb dating from the middle belt yields an age of ca.415 Ma,which is similar to south NQT.The thrusting displacement is estimated at ca.48 km by structural interpretation of deep profiles.The timing is constrained in the early stage of the Early Cretaceous by the formation of simultaneous growth strata.We suggest that the NQTS has resulted from the far-field effect of the Lhasa-Qiangtang collision,and the Yumu Shan is uplifted by the superposed Cenozoic deformation.

Performance-based seismic design of nonstructural building components:The next frontier of earthquake engineering

With the development and implementation of performance-based earthquake engineering,harmonization of performance levels between structural and nonstructural components becomes vital. Even if the structural components of a building achieve a continuous or immediate occupancy performance level after a seismic event,failure of architectural,mechanical or electrical components can lower the performance level of the entire building system. This reduction in performance caused by the vulnerability of nonstructural components has been observed during recent earthquakes worldwide. Moreover,nonstructural damage has limited the functionality of critical facilities,such as hospitals,following major seismic events. The investment in nonstructural components and building contents is far greater than that of structural components and framing. Therefore,it is not surprising that in many past earthquakes,losses from damage to nonstructural components have exceeded losses from structural damage. Furthermore,the failure of nonstructural components can become a safety hazard or can hamper the safe movement of occupants evacuating buildings,or of rescue workers entering buildings. In comparison to structural components and systems,there is relatively limited information on the seismic design of nonstructural components. Basic research work in this area has been sparse,and the available codes and guidelines are usually,for the most part,based on past experiences,engineering judgment and intuition,rather than on objective experimental and analytical results. Often,design engineers are forced to start almost from square one after each earthquake event: to observe what went wrong and to try to prevent repetitions. This is a consequence of the empirical nature of current seismic regulations and guidelines for nonstructural components. This review paper summarizes current knowledge on the seismic design and analysis of nonstructural building components,identifying major knowledge gaps that will need to be filled by future research. Furthermore,considering recent trends in earthquake engineering,the paper explores how performance-based seismic design might be conceived for nonstructural components,drawing on recent developments made in the field of seismic design and hinting at the specific considerations required for nonstructural components.

Seismic isolation analysis of FPS bearings in spatial lattice shell structures

A theoretical model of a friction pendulum system （FPS） is introduced to examine its application for the seismic isolation of spatial lattice shell structures. An equation of motion of the lattice shell with FPS bearings is developed. Then, seismic isolation studies are performed for both double-layer and single-layer lattice shell structures under different seismic input and design parameters of the FPS. The influence of frictional coefficients and radius of the FPS on seismic performance are discussed. Based on the study, some suggestions for seismic isolation design of lattice shells with FPS bearings are given and conclusions are made which could be helpful in the application of FPS.

Performance-based methodology for assessing seismic vulnerability and capacity of buildings

This paper presents a performance-based methodology for the assessment of seismic vulnerability and capacity of buildings. The vulnerability assessment methodology is based on the HAZUS methodology and the improved capacity- demand-diagram method. The spectral displacement (Sd) of performance points on a capacity curve is used to estimate the damage level of a building. The relationship between Sd and peak ground acceleration (PGA) is established, and then a new vulnerability function is expressed in terms of PGA. Furthermore, the expected value of the seismic capacity index (SCev) is provided to estimate the seismic capacity of buildings based on the probability distribution of damage levels and the corresponding seismic capacity index. The results indicate that the proposed vulnerability methodology is able to assess seismic damage of a large number of building stock directly and quickly following an earthquake. The SCev provides an effective index to measure the seismic capacity of buildings and illustrate the relationship between the seismic capacity of buildings and seismic action. The estimated result is compared with damage surveys of the cities of Dujiangyan and Jiangyou in the M8.0 Wenchuan earthquake, revealing that the methodology is acceptable for seismic risk assessment and decision making. The primary reasons for discrepancies between the estimated results and the damage surveys are discussed.

Frequency dependent polarization analysis of ambient seismic noise recorded at a broadband seismometer in the central United States

We present a new approach to polarization analysis of seismic noise recorded by three-component seismometers. It is based on statistical analysis of frequency-dependent particle motion properties determined from a large number of time windows via eigenanalysis of the 3-by-3, Hermitian, spectral covariance matrix. We applied the algorithm to continuous data recorded in 2009 by the seismic station SLM, located in central North America. A rich variety of noise sources was observed. At low frequencies （〈0.05 Hz） we observed a tilt-related signal that showed some elliptical motion in the horizontal plane. In the microseism band of 0.05-0.25 Hz, we observed Rayleigh energy arriving from the northeast, but with three distinct peaks instead of the classic single and double frequency peaks. At intermediate frequencies of 0.5-2.0 Hz, the noise was dominated by non-fundamental-mode Rayleigh energy, most likely P and Lg waves. At the highest frequencies （〉3 Hz）, Rayleigh-type energy was again dominant in the form of Rg waves created by nearby cultural activities. Analysis of the time dependence of noise power shows that a frequency range of at least 0.02-1.0 Hz （much larger than the microseism band） is sensitive to annual, meteorologically induced sources of noise.

Nonlinear earthquake response analysis and energy calculation for seismic slit shear wall structures

Based on the concept of structural passive control,a new type of slit shear wall,with improved seismic performance when compared to an ordinary solid shear wall,was proposed by the authors in 1996.The idea has been verified by a series of pseudo-static and dynamic tests.In this paper a macro numerical model is developed for the wall element and the energy dissipation device.Then,nonlinear time history analysis is carried out for a 10-story slit shear wall model tested on a shaking table.Furthermore,the seismic input energy and the individual energy dissipated by the components are calculated by a method based on Newmark-β assumptions for this shear wall model,and the advantages of this shear wall are further demonstrated by the calculation results from the viewpoint of energy.Finally,according to the seismic damage criterion on the basis of plastic accumulative energy and maximum response,the optimal analysis is carried out to select design parameters for the energy dissipation device.

Seismic response of continuous span bridges through fiber-based finite element analysis

It is widely recognized that nonlinear time-history analysis constitutes the most accurate way to simulate the response of structures subjected to strong levels of seismic excitation. This analytical method is based on sound underlying principles and has the capability to reproduce the intrinsic inelastic dynamic behavior of structures. Nonetheless, comparisons with experimental results from large-scale testing of structures are still needed, in order to ensure adequate levels of confidence in this numerical methodology. The fiber modelling approach employed in the current endeavor inherently accounts for geometric nonlinearities and material inelasticity, without a need for calibration of plastic hinges mechanisms, typical in concentrated plasticity models. The resulting combination of analysis accuracy and modelling simplicity, allows thus to overcome the perhaps not fully justifiable sense of complexity associated to nonlinear dynamic analysis. The fiber-based modelling approach is employed in the framework of a finite element program downloaded from the Intemet for seismic response analysis of framed structures. The reliability and accuracy of the program are demonstrated by numerically reproducing pseudo-dynamic tests on a four span continuous deck concrete bridge. Modelling assumptions are discussed, together with their implications on numerical results of the nonlinear time-history analyses, which were found to be in good agreement with experimental results.

Vector-intensity measure based seismic vulnerability analysis of bridge structures

This paper presents a method for seismic vulnerability analysis of bridge structures based on vector-valued intensity measure （viM）, which predicts the limit-state capacities efficiently with multi-intensity measures of seismic event. Accounting for the uncertainties of the bridge model, ten single-bent overpass bridge structures are taken as samples statistically using Latin hypercube sampling approach. 200 earthquake records are chosen randomly for the uncertainties of ground motions according to the site condition of the bridges. The uncertainties of structural capacity and seismic demand are evaluated with the ratios of demand to capacity in different damage state. By comparing the relative importance of different intensity measures, Sa（T1） and Sa（T2） are chosen as viM. Then, the vector-valued fragility functions of different bridge components are developed. Finally, the system-level vulnerability of the bridge based on viM is studied with Duunett- Sobel class correlation matrix which can consider the correlation effects of different bridge components. The study indicates that an increment IMs from a scalar IM to viM results in a significant reduction in the dispersion of fragility functions and in the uncertainties in evaluating earthquake risk. The feasibility and validity of the proposed vulnerability analysis method is validated and the bridge is more vulnerable than any components.

Rockfall seismic features analysis based on in situ tests:frequency, amplitude, and duration

In the seismic event classification,determining the seismic features of rockfall is significantly important for the automatic classification of seismic events because of the huge amount of raw data recorded by seismic stations in continuous monitoring. At the same time, the rockfall seismic features are still not completely understood.This study concentrates on the rockfall frequency content, amplitude(ground velocity), seismic waveform and duration analysis, of an artificial rockfall test at Torgiovannetto(a former quarry in Central Italy). A total of 90 blocks were released in the test, and their seismic signals and moving trajectories were recorded by four tri-axial seismic stations and four cameras, respectively. In the analysis processing,all the artificial rockfall signal traces were cut separately and the seismic features were extracted individually and automatically. In this study, the relationships between a) frequency content and impacted materials, b) frequency content and the distance between block releasing position and seismic station(source-receiver distance) were discussed. As a result, we found that the frequency content of rockfall focuses on 10-60 Hz and 80-90 Hz within a source-receiver distance of 200 m, and it is well correlated with impacted material and source-receiver distance. To evaluate the difference between earthquake and rockfall, 23 clear earthquake signals recorded in a seven month-long continuous seismic monitoring, carried out with the four seismic stations, were picked out, according to the Italian national earthquakes database(INGV). On these traces we performed the same analysis as in the artificial rockfall traces, and two parameters were defined to separate rockfall events from earthquake noise. The first one, the amplitude ratio, is related to the amplitude variation of rockfall between two stations and is greater than that of earthquakes, because of the higher attenuation occurring for rockfall events, which consists in high frequencies whereas for earthquakes it consists in low frequencies. The other parameter, the shape of waveform of signal trace, showed a significant difference between rockfall and earthquake and that could be a complementary feature to discriminate between both. This analysis of artificial rockfall is a first step helpful to understand the seismic characteristics of rockfall, and useful for rockfall seismic events classification in seismic monitoring of slope.

Slope stability analysis under seismic load by vector sum analysis method

The vibration characteristics and dynamic responses of rock and soil under seismic load can be estimated with dynamic finite element method （DFEM）. Combining with the DFEM, the vector sum analysis method （VSAM） is employed in seismic stability analysis of a slope in this paper. Different from other conventional methods, the VSAM is proposed based on the vector characteristic of force and current stress state of the slope. The dynamic stress state of the slope at any moment under seismic load can he obtained by the DFEM, thus the factor of safety of the slope at any moment during earthquake can be easily obtained with the VSAM in consideration of the DFEM. Then, the global stability of the slope can be estimated on the basis of time-history curve of factor of safety and reliability theory. The VSAM is applied to a homogeneous slope under seismic load. The factor of safety of the slope is 1.30 under gravity only and the dynamic factor of safety under seismic load is 1.21. The calculating results show that the dynamic characteristics and stability state of the slope with input ground motion can be actually analyzed. It is believed that the VSAM is a feasible and practical approach to estimate the dynamic stability of slopes under seismic load.