Probabilistic evaluation of liquefaction potential under earthquake loading

This paper presents a procedure to perform the risk analysis for ground failure by liquefaction. The first part of this study describes the differential equation of a smooth hysteretic model to characterize the behavior of the soil under random loading. The parameters of the proposed model to represent the experimental relationship are discussed. The second part of this study is to develop a method to calculate the probability that a specified volume of soil will liquefy at a given depth in the deposit. The liquefaction is defined as the result of cumulative damage caused by seismic loading. The fatigue life of soil can be determined on the basis of the N---S relationship and Miner's cumulative damage law. The rain-flow method is used to count the number of cycles of stress response of the soil deposit. Finally, the probability of liquefaction is obtained by integration over all the possible ground motion and the fragility curves of liquefaction potential. The sensitivity of the reliability against liquefaction to soil system parameters is also examined.     

Analysis of the spatial variation of seismic waves and ground movements from smart-1 array data

Specially designed arrays of strong motion seismographs located near earthquake sources are required for engineering studies of the near-source properties and the spatial variation of seismic waves. The SMART-1 array in Taiwan provides good records for this type of study. Careful study of the observed strong motion data permits the identification of wave types, directions and apparent wave velocities. In this paper, a principal direction ratio R (f,α) is defined; this indicates the principal direction of the motion (along a nearly straight line) within the range 0 < R < 1. Vertical motion of the ground is also included in this study. Orbit spectrum analysis is used to verify the identification of wave directions and wave types. The spatial variation of seismic waves along the principal direction is studied. From frequency-domain analysis, mathematical models of the spatial variation of ground displacement are developed using a wave-number spectrum and the cross-spectral density function between two spatial coordinates; these models in turn can provide two alternative models for the random vibration analysis of extensive structures subject to multiple point seismic excitation. The SMART-1 array data gathered during the January 29, 1981 earthquake also are used to demonstrate calculation of the ground strains and differential movements of the array site. From time-domain analysis, the spatial variation of seismic waves is defined for ground motion along the identified principal direction. The time variation of evolutionary spectra characterized by frequency-dependent parameters is used for this formulation. The SMART-1 array data again form the basis for discussion of the spatial variation of model parameters.     

A review of the effects of particle types on oil-suspended particulate matter aggregate formation

Oil-suspended particulate matter aggregate (OSA) can form naturally when oil and particles interact. The interaction between oil and suspended particulate matter makes oil less sticky, and facilitates its dispersion in the water column. The high oil-water surface contact enhances the biodegradation of oil and thus increases the efficiency of remediation processes. There are many factors that affect OSA formation, but, particle type is one of the most important. Because different particle types have different physical, chemical, and biological properties, their interactions with oil differ greatly. Particle properties such as interlayer spaces, hydrophobicity, surface charges, polarity, organic content, and size affect the interactions between materials and oil. These different interactions determine the type, buoyancy, size, and stability of OSA that forms, thus determining its fate in the environment. This review provides a current understanding of (1) OSA formation mechanisms, (2) sources and classes of marine materials, (3) oil-particle interactions, (4) material properties and their effects on oil interaction, and (5) future research needs.     

Seismic Hazard Assessment for the Taiwan Region on the Basis of Recent Strong-Motion Data and Prognostic Zonation of Future Earthquakes

The paper describes an integrated approach to seismic hazard assessment, which was applied for the Taiwan region. First, empirical modelsfor ground motion estimation in the region were obtained on the basisof records from recent (1993-1999) earthquakes. The databaseincludes strong-motion data collected during the recent Chi-Chiearthquake (M=7.6, 21 September 1999) and large (M=6.8)aftershocks. The ground-motion database was also used for evaluationof generalised site amplification functions for typical soil classes(B, C and D). Second, the theoretical seismic catalogue (2001–2050)for the Taiwan region had been calculated using the 4D-model(location, depth, time) for dynamic deformation of the Earth' crustand 5D-model (location, depth, time, magnitude) for seismic process.The models were developed on the basis of available geophysical andgeodynamic data that include regional seismic catalogue. Third, theregion & site & time-dependent seismic analysis, which is basedon schemes of probable earthquake zones evaluated from the theoreticalcatalogue, regional ground motion models, and local site responsecharacteristics, has been performed. The seismic hazard maps arecompiled in terms of Peak Ground Acceleration (PGA) and ResponseSpectra (RS) amplitudes. The maps show distribution of amplitudesthat will not be exceeded with certain probability in condition oftypical soil classes during all possible earthquakes that may occur inthe region during time period of 2003–2025. The approach allowsintroducing new parameter that describes dependency of seismichazard on time, so-called 'period of maximum hazard'. Theparameter shows the period, during which every considered sitewill be subjected by the maximum value of ground motioncharacteristic (PGA or RS).     

Two‐stage damage detection algorithms of structure using modal parameters identified from recursive subspace identification

Structural damage assessment under external loading, such as earthquake excitation, is an important issue in structural safety evaluation. In this regard, appropriate data analysis and feature extraction techniques are required to interpret the measured data and to identify the state of the structure and, if possible, to detect the damage. In this study, the recursive subspace identification with Bona‐fide LQ renewing algorithm (RSI‐BonaFide‐Oblique) incorporated with moving window technique is utilized to identify modal parameters such as natural frequencies, damping ratios, and mode shapes at each instant of time during the strong earthquake excitation. From which the least square stiffness method (LSSM) combined with the model updating technique, called efficient model correction method (EMCM), is used to estimate the first‐stage system stiffness matrix using the simplified model from the previously identified modal parameters (nominal model). In the second stage, 2 different damage assessment algorithms related to the nominal system stiffness matrix were derived. First, the model updating technique, called EMCM, is applied to correct the nominal model by the newly identified modal parameters during the strong motion. Second, the element damage index can be calculated using element damage index method (EDIM) to quantify the damage extent in each element. Verification of the proposed methods through the shaking table test data of 2 different types of structures and a building earthquake response data is demonstrated to specify its corresponding damage location, the time of occurrence during the excitation, and the percentage of stiffness reduction.     

Collapse of a nonductile concrete frame: Shaking table tests

Post‐earthquake reconnaissance has reported the vulnerability of older reinforced concrete (RC) columns lacking details for ductile response. Research was undertaken to investigate the full‐range structural hysteretic behavior of older RC columns. A two‐dimensional specimen frame, composed of nonductile and ductile columns to allow for load redistribution, was subjected to a unidirectional base motion on a shaking table until global collapse was observed. The test demonstrates two types of column failure, including flexure‐shear and pure flexural failure. Test data are compared with various simplified assessment models commonly used by practicing engineers and researchers to identify older buildings that are at high risk of structural collapse during severe earthquake events. Comparison suggests that ASCE/SEI 41‐06 produces very conservative estimates on load–deformation relations of flexure‐shear columns, while the recently proposed ASCE/SEI 41‐06 update imposes significant modifications on the predictive curve, so that improved accuracy has been achieved. Copyright © 2008 John Wiley & Sons, Ltd.     

Application of subspace identification technique to long‐term seismic response monitoring of structures

This paper addresses the issue of structural system identification using earthquake‐induced structural response. The proposed methodology is based on the subspace identification algorithm to perform identification of structural dynamic characteristics using input–output seismic response data. Incorporated with subspace identification algorithm, a scheme to remove spurious modes is also used to identify real system poles. The efficiency of the proposed method is shown by the analysis of all measurement data from all measurement directly. The recorded seismic response data of three structures (one 7‐story RC building, one midisolation building, and one isolated bridge), under Taiwan Strong Motion Instrumentation Program, are analyzed during the past 15 years. The results present the variation of the identified fundamental modal frequencies and damping ratios from all the recorded seismic events that these three structures had encountered during their service life. Seismic assessment of the structures from the identified system dynamic characteristics during the period of their service is discussed. Copyright © 2014 John Wiley & Sons, Ltd.     

Tracking modal parameters from building seismic response data using recursive subspace identification algorithm

Tracking modal parameters and estimating the current structural state of a building from seismic response measurements, particularly during strong earthquake excitations, can provide useful information for building safety assessment and the adaptive control of a structure. Therefore, online or recursive identification techniques need to be developed and implemented for building seismic response monitoring. This paper develops and examines different methods to track modal parameters from building seismic response data. The methods include recursive data‐driven subspace identification (RSI‐DATA) using Givens rotation algorithm, and RSI‐DATA using Bona fide algorithm. The question on how well the results of RSI‐DATA reflect the real condition is investigated and verified with a bilinear SDOF simulation study. Time‐varying modal parameters of a four‐story reinforced concrete school building are identified based on a series of earthquake excitations, including several seismic events, large and small. Discussions on the different methods' ability to track the time‐varying modal parameters are presented. The variation of the identified building modal frequencies and damping ratios from a series of event‐by‐event seismic responses, particularly before and after retrofitting of the building is also discussed. Copyright © 2017 John Wiley & Sons, Ltd.     

Empirical model for estimating Fourier amplitude spectra of ground acceleration in Taiwan region

A collection of ground‐motion recordings (1070 acceleration records) of moderate (5.1⩽M_L⩽6.5) earthquakes obtained during the execution of the Taiwan Strong Motion Instrumentation Program (TSMIP) since 1991 was used to study source scaling model and attenuation relations for a wide range of earthquake magnitudes and distances and to verify the models developed recently for the Taiwan region. The results of the analysis reveal that the acceleration spectra of the most significant part of the records, starting from S‐wave arrival, can be modelled accurately using the Brune's ω‐squared source model with magnitude‐dependent stress parameter Δσ, that should be determined using the recently proposed regional relationships between magnitude (M_L) and seismic moment (M₀) and between M₀ and Δσ. The anelastic attenuation Q of spectral amplitudes with distance may be described as Q=225 ƒ^1.1 both for deep (depth more than 35 km) and shallow earthquakes. The source scaling and attenuation models allow a satisfactory prediction of the peak ground acceleration for magnitudes 5.1⩽M⩽6.5 and distances up to about 200 km in the Taiwan region, and may be useful for seismic hazard assessment. Copyright © 2000 John Wiley & Sons, Ltd.     

Seismic hazard and uncertainty analysis of the Taiwan area

The objectives of this paper are (1) to obtain estimates on the effect of uncertainties of the hazard model, and (2) to evaluate the seismic hazard in Taiwan for structural analysis and design purposes. The seismic hazard in the Taiwan area is presented in terms of an iso-acceleration map. Such a map is developed for return periods of peak ground acceleration of 225 years and 475 years. The contour map of b-values and mean occurence rates for this region is also presented. Uncertainty analyses of model parameters in hazard analysis are concentrated on the analysis of dispersion of PGA values and the probabilistic modeling of stationary and nonstationary Poisson models of occurrences. Th e overall results are considered to be conservative since for most uncertainty analyses the more conservative values are used.