Seismic Motion Research Articles

System reliability analysis of Cable-stayed bridge using pair Copula and ICLHS with unequal weights under seismic loading

Considering the correlation between different components, quantifying the failure probability of a high-level hyper-static structure under non-stationary seismic loading is of great significance in guiding earthquake-resistant design and fortification of bridges. In this study, a high-pier and large-span railway cable-stayed bridge in the mountainous region were employed to evaluate the seismic performance. First, the basic theoretical methods of the non-stationary random seismic motion simulation, equivalent extreme value theory based on extreme value distribution (EVD), improved Latin hypercube sampling (ICLHS) with unequal weighted fractional moment estimation and Pair Copula technique are presented. Second, the finite element analysis model of the cable-stayed bridge was established based on the OpenSees platform, the seismic response samples of the bridge were calculated under three levels of seismic ground motion to determine the first passage failure probability of the main components of the bridge under the design failure modes. Finally, the Pair Copula function was introduced to estimate the seismic system reliability of the cable-stayed bridge. The optimal Copula type that can describe the correlation between the bridge components was selected through goodness-of-fit tests such as AIC, BIC, and the minimum distance method. The failure probability of the entire system was obtained by calculating the selected Copula function type under seismic loading.

Fiber-optic gyroscope for rotational seismic ground motion monitoring of the Campi Flegrei volcanic area

The real-time monitoring of densely populated areas with high seismic and volcanic risk is of crucial importance for the safety of people and infrastructures. When an earthquake occurs, the Earth surface experiences both translational and rotational motions. The latter are usually not monitored, but their measurement and characterization are essential for a full description of the ground motion. Here we present preliminary observational data of a high-sensitivity rotational sensor based on a 2-km-long fiber-optic Sagnac gyroscope, presently under construction in the middle of the Campi Flegrei Volcanic Area (Pozzuoli, Italy). We have evaluated its performance by analyzing data continuously recorded during an acquisition campaign of five months. The experimental setup was composed of a digital nine-component seismic station equipped with both a rotational sensor and conventional seismic sensors (seismometers, accelerometers, and tiltmeters). During this experiment we detected seismic noise and ground rotations wavefield induced by small to medium local earthquakes (M D <3). The prototype gyroscope shows a very promising sensitivity in the range of 5×10−7−8×10−9rad/s/Hz over the frequency bandwidth 5 mHz–50 Hz. Future upgrades and perspectives are discussed.

Historical Evolution of the Impact of Seismic Incident Angles on the Safety Assessment of Various Building Construction Typologies

In the existing building stock, typically characterised by a high degree of irregularity, the effects of earthquakes are strongly dependent on the epicentre–structure direction and the angle of incidence of the seismic motion. However, the scientific community has not yet reached a unanimous consensus on the evaluation of the effects of seismic incidence angles. Therefore, this paper conducts an extensive investigation of the international literature on current methods to consider seismic directionality, systematically reviewing more than 80 publications on this topic. Following a brief overview of the problem and an analysis of the initial developments of the multidirectionality concept of seismic input, a state-of-the-art review is presented based on the considered analysis methods, specifically response spectrum analysis, nonlinear static analysis, and nonlinear response history analysis. Moreover, the adoption of multidirectional seismic input in popular codes and standards is presented and discussed. This study provides the first comprehensive synthesis of research on the seismic incidence angles across diverse building typologies, offering crucial insights for future code revisions and highlighting significant gaps in current analytical methods and standards, thereby setting a new direction for subsequent empirical investigations. Specifically, the extensive state-of-the-art review revealed that, until now, the evaluation of the angle of incidence was primarily conducted on existing reinforced concrete buildings with a limited number of storeys, analysed with nonlinear response history analysis. This underscores the need for future research to extensively investigate the impact of the angle of incidence on other types of construction typologies.

Seismic mitigation analysis of three-dimensional base-isolated nuclear structures with soil-dependent isolation system under extreme earthquakes

The Fukushima nuclear accident in March 2011 has raised more stringent requirements for the ability of nuclear power plants (NPPs) to withstand extreme disasters such as extreme earthquakes. To address this challenge, this study proposed a new hybrid multi-dimensional passive control system that combines three-dimensional base isolation and geotechnical seismic isolation. In this system, high-damper rubber bearings were employed to isolate horizontal seismic motions, whereas disc springs, accounting for frictional forces, were used to mitigate vertical seismic motions. Additionally, a geotechnical seismic isolation system arranged under the base mat was leveraged to further reduce the seismic response of the superstructure. Taking a large-scale reactor building in complex layered sites as an example, the effectiveness of the hybrid multi-dimensional isolation system was studied using a wave motion method. The research results indicate that when subjected to extreme earthquakes, hybrid multi-dimensional isolation exhibits superior isolation performance compared with three-dimensional base isolation, which can significantly reduce the three-dimensional dynamic response of nuclear structures without increasing the displacement of the isolation bearings. The proposed hybrid multi-dimensional isolation system offers an effective solution to the challenges faced by three-dimensional isolation systems, thereby enhancing the seismic resistance of nuclear structures against extreme earthquakes.

Time-frequency analysis of nonlinear structures under nonstationary excitation by directly solving structural dynamic equation through absolute displacement

The time-frequency evolution of seismic ground motion energy under strong earthquakes, which can be described by the evolutionary power spectral density (EPSD) matrix, will amplify the most unfavorable responses of members in nonlinear structural systems. However, the existing general finite element software fails to reduce the dimension of the EPSD matrix and then decouple time from frequency in frequency domain analysis. This failure greatly limits the simulation and calculation of multi-dimensional nonlinear structures with multiple degrees of freedom by the stochastic vibration theory. In this study, the pseudo excitation method (PEM), which can direct solve structural dynamic equations using absolute displacement, is used to derive the self-spectral and cross-spectral densities of seismic ground motions in one-dimensional and multi-dimensional cases. In MATLAB, the EPSD matrix is transformed into a four-dimensional uniform modulation excitation matrix independent of the time variable. Nonlinear iteration is conducted and the precise integration method (PIM) is employed through a transient analysis module. Furthermore, a simple one-dimensional example is presented to verify the correctness of the proposed algorithm herein. Finally, as an example, we select important parameters of the coherence model at different supports of a half-through three-span concrete-filled steel tubular (CFST) tied-arch bridge in view of real engineering conditions, and calculate the time-varying power spectrum and variance of the responses of arch ribs, arch feet and the bridge deck system. The results show that the theoretical derivation and the algorithm constructed in this paper can simulate the time-delay effect of the nonlinear structural system. This proves the feasibility of general finite element software in solving and analyzing random responses to multi-dimensional and multi-support non-stationary excitation of nonlinear structures.

Bayesian maximum entropy interpolation analysis for rapid assessment of seismic intensity using station and ground motion prediction equations

In this paper, we explored the combination of seismic station data and ground motion prediction equations (GMPE) to predict seismic intensity results by using Bayesian Maximum Entropy (BME) method. The results indicate that: 1) In earthquake analysis in Japan, soft data has predicted higher values of intensity in disaster areas. BME corrected this phenomenon, especially near the epicenter. Meanwhile, for earthquakes in the United States, BME corrected the erroneous prediction of rupture direction using soft data. 2) Compared with other spatial interpolation methods, the profile results of Japan earthquake and Turkey earthquake show that BME is more consistent with ShakeMap results than IDW and Kriging. Moreover, IDW has a low intensity anomaly zone. 3) The BME method overcomes the phenomenon that the strength evaluation results do not match the actual failure situation when the moment magnitude is small. It more accurately delineates the scope of the disaster area and enriches the post-earthquake processing of disaster area information and data. BME has a wide range of applicability, and it can still be effectively used for interpolation analysis when there is only soft data or few sites with data available.

Seismic behavior analysis of long-span cable-stayed bridge under bi-directional near-fault ground motions

In this paper, the seismic behaviors of a long-span cable-stayed bridge subjected to bi-directional near-fault ground motions are analyzed. The cable-stayed bridge with a main span of 432 m is selected as the research object and appropriately modeled using the OpenSees software. Ten sets of historical near-fault and far-field ground motions associated with the same earthquake event and site class are utilized as input excitations. Nonlinear time-history analyses are conducted to explore three key aspects: the influence of ground motion directionality, the effectiveness of the percentage rule and simplified Square-Root-of-Sum-of-Squares (SRSS) rule, and the main influential factors. The findings reveal that near-fault ground motions induce larger seismic responses compared to far-field ground motions. Moreover, the effects of ground motion directionality are more pronounced in near-fault scenarios than in far-field events. The percentage rule and simplified SRSS rule demonstrate minor significance in predicting the critical seismic responses. The bi-directional peak ground velocity (PGV) of seismic motions in the principal directions of the cable-stayed bridge may emerge as the primary factor influencing the seismic responses, making it a suitable intensity measure for bi-directional seismic fragility assessment.

Simulating non-stationary and non-Gaussian cross-correlated fields using multivariate Karhunen–Loève expansion and L-moments-based Hermite polynomial model

In practical engineering, cross-correlated random fields are often used to model structural materials or random loads containing multiple correlations. Effective and accurate simulation of these cross-correlation fields is an important prerequisite for subsequent reliability analysis and uncertainty quantification of complex systems. Therefore, this paper proposes a new simulation method for non-stationary non-Gaussian cross-correlated random fields to satisfy realistic engineering requirements. In this new method, L-moments-based Hermite polynomials model (LHPM) is extended to non-stationary non-Gaussian cross-correlated fields. Then, a transformation model from non-stationary non-Gaussian correlation matrix function (CMF) to the underlying Gaussian CMF is explicitly given. Multivariate K-L expansion is further used to approximate the underlying Gaussian cross-correlated fields, where the eigenvalue and eigenfunctions are estimated via the Nyström method with uniform weights. The approximation permits the application of few random variables to characterize the entire Gaussian cross-correlated fields. Ultimately, the simulated underlying Gaussian cross-correlated fields are mapped back to the target non-stationary non-Gaussian cross-correlated fields based on LHPM. Three typical examples, including exponential kernel, Wiener process kernel and spatially varying non-Gaussian and nonstationary seismic ground motions, are used to validate the effectiveness of the proposed method. The source code is readily available at: https://github.com/zhaozhao23/Simulation-of-non-stationary-and-non-Gaussian-cross-correlated-fields.

Seismic resistance analysis of a plywood house structure based on ABAQUS fracture mechanics

Plywood is a more affordable material for the residential construction industry due to its lower cost and easy availability. However, earthquakes are the biggest barrier for wood structures to stand for a long time, as the materials can easily deform or collapse when subjected to seismic forces. This paper aims to explore the performance of plywood structures during an active earthquake by conducting a series of tests and simulations. In this study, plywood property was explored in the four-point bending experiment to determine its stress-strain relationship, and the results were used to identify the materials young modulus, a key parameter for designing more resilient structures. In addition, a finite element model of a plywood house structure was developed in the ABAQUS application, and the seismic movements with five different velocity values were applied to determine the stress distribution and deflections of the structure. The results provide insights into the dynamic behavior of plywood structures under different seismic movements.

Video camera and seismic monitoring of water bulge explosion at Strokkur Geyser, Iceland

Eruptions of volcanoes and geysers share many fundamental similarities: for example, they are manifestations of Earth’s geothermal energy, involving the pressure-driven expulsion of fluids from the Earth’s interior. However, while volcanoes can produce spectacular lava bubbles that burst, water bubbles are rarely observed on the surface of geysers. It is still unclear why some of these low-viscosity geyser systems produce none, while others produce them regularly. There is no quantification of the size, speed, and height of these bubbles at geysers, which is the gap we fill here. Strokkur creates a water bulge in its surface pool (bulge stage). When the bulge bursts, water is ejected into the air (jet stage). The steam then continues to rise buoyantly and drift away (drift stage). Here we study the evolution of the three stages using records from video camera campaigns and a local seismic network. We find that larger bulges are associated with larger ascent velocities and cause larger jet heights. As energy is channeled into a high jet, small seismic ground motions are recorded. The bulge formation itself is barely visible seismically. Our work suggests that the 0.74±0.27 s-long bulge stage can be used as a first-order proxy for predicting eruption height. This study might also be relevant for understanding fluid dynamics in volcanic systems.

Experimental Study on Erosion Process of Expressway Embankment Subjected to Tsunami After Earthquake

On March 11, 2011, the Great East Japan Earthquake triggered tsunamis that reached extensive areas along Japan’s Pacific coast. There have been instances where embankments built on plains for expressways mitigated the impact of tsunami damage. In the vicinity of the Sendai-tobu highway, the presence of an embankment approximately 10 m high altered the course of the advancing tsunami, thereby preventing flooding. Establishing a multiplied defense system using road embankments necessitates understanding the deformation and collapse mechanisms of road embankments impacted by tsunamis following seismic motion. In this study, overtopping experiments were conducted by first applying seismic motion to model embankments, followed by introducing the first wave of breaking bores, and then simulating prolonged overtopping by the tsunami. The experimental findings indicated that within the embankments impacted by the tsunami, there was an immediate increase in what is presumed to be pore air pressure following the arrival of the breaking bores, followed by a rise in pore water pressure during subsequent overtopping. Moreover, embankments subjected to seismic motion exhibited accelerated erosion following the overtopping. These results imply that when embankments settle due to an earthquake, leading to relatively higher anticipated inundation depths and the potential for overtopping, it is crucial to implement measures to prevent the settlement of the crest for embankments expected to serve as part of a multiplied defense system.

Improved pseudostatic analysis of pile response to seismic ground movement: Multiple characteristic ground displacement profiles

Piles may be subjected to significant kinematic loading from ground movement under earthquakes. Pseudostatic analysis, which is based on the Winker beam model, traditionally adopts a single ground displacement profile, known as the characteristic ground displacement profile, for assessments of pile response under this loading condition. This study conducts dynamic analyses with different soil, pile-head-fixity, and seismic conditions and generates displacement, moment, and shear force response envelopes to evaluate the appropriateness of the characteristic displacement profiles of four common methods. These methods give inconsistent predictions because the maximum responses at different depths in the pile response envelopes occur at different points in time. To address this issue, we propose using multiple characteristic ground displacement profiles, which are determined by the maximum ground relative deformation or rotation. These multiple profiles, when used in pseudostatic analyses, give consistent predictions with dynamic pile response envelopes in most scenarios, thus overcoming the limitations of conventional single-profile methods.

Influence of frequency characteristics of seismic ground motion on the fragility of an RC frame structure

Ground-motion frequency characteristics have a significant impact on structural damage. A quantitative description of the impact of ground-motion spectral features on structures is an urgent topic that requires solutions. In this study, the ratio of peak displacement to peak acceleration (R d2a ) is used as a representative index of the ground-motion spectral characteristics. The seismic records collected from the Kik-net strong-motion array in Japan are grouped according to their R d2a values. By adopting the incremental dynamic analysis method, the seismic fragility curves of a six-story RC frame structure loaded by inputting seismic records with different R d2a values are calculated using the OpenSees platform. The influence of R d2a values on the structural fragility is analyzed. Ground-motion records of liquefied sites are collected, and the peak ratios of seismic records at the liquefied sites are analyzed to reveal the differences in the fragility of the six-story RC structures at liquefied and non-liquefied sites. The results of this study are expected to serve as a reference for the evaluation of structural fragility and to understand the effect of ground-motion frequency characteristics on the fragility of RC frame structures.

Seismic response of tripod suction bucket foundation for offshore wind turbine considering multidirectional earthquake loads

Tripod suction buckets offer several advantages as offshore wind turbine (OWT) foundations, including cost-effective installation, high stability, and resistance to overturning. In seismically active regions, OWTs experience horizontal and vertical directions seismic excitations. This study employed an advanced three-dimensional nonlinear finite element analysis to investigate the influence of vertical seismic motions on the dynamic behavior of tripod bucket foundations installed in nonhomogeneous clay deposits. By adopting a simple kinematic hardening model to capture the foundation-soil nonlinear dynamic interaction, parametric analyses were conducted to study the effect of vertical seismic motions, considering several earthquake types and undrained shear strengths of the clay deposit. The results show that the influence of the vertical seismic motion on the dynamic response and deformation of the nacelle depends on the dominant frequency of the vertical seismic motion and the overturning resistance capacity of the OWT in different horizontal directions. The deformation mechanisms of an OWT vary with the strengths of the clay deposits. Based on these analyses, this study offers insights into the influence of vertical seismic motions on the dynamic response of tripod suction bucket foundations.

Assessment of Site Effects and Numerical Modeling of Seismic Ground Motion to Support Seismic Microzonation of Dushanbe City, Tajikistan

In the territory of Dushanbe city, the capital of Tajikistan, detailed geological and geophysical data were collected during geophysical surveys in 2019–2020. The data comprise 5 microtremor array measurements, 9 seismic refraction tomography profiles, seismological data from 5 temporary seismic stations for standard spectral ratio calculations, 60 borehole datasets, and 175 ambient noise measurements. The complete dataset for Dushanbe was used to build a consistent 3D geologic model of the city with a size of 12 × 12 km2. The results of the seismological and geophysical surveys were compared and calibrated with borehole data to define the boundaries of each layer in the study area. The Leapfrog Works software was utilized to create a 3D geomodel. From the 3D geomodel, we extracted six 12 km long 2D geological cross-sections. These 2D geological cross-sections were used for 2D dynamic numerical modeling with the Universal Distinct Element Code software to calculate the local seismic response. Finally, the dynamic numerical modeling results were compared with the amplification functions obtained from the seismological and ambient noise data analysis. The 2D dynamic numerical modeling results allowed a better assessment of the site effects in the study area to support seismic microzonation and the determination of local peak ground acceleration changes in combination with regional seismic hazard maps. In addition, our results confirm the strong seismic amplification effects noted in some previous studies, which are attributed to the influence of local topographic and subsurface characteristics on seismic ground motions.

Pseudo-static Winkler springs for longitudinal underground structures subjected to shear waves

ABSTRACT Shear waves in underground longitudinal structures impose deformations that need to be accounted for in structural design. Simplified approaches for such structures typically consider the Winkler foundation model, assuming Euler-Bernoulli beam theory, which neglects shearing-induced distortions, especially significant when shear waves are a dominant component of the seismic motion. In order to overcome these limitations, Timoshenko beam models have been proposed in the literature. These approaches however depend on an appropriate determination of the ground springs. Existing analytical formulations often assume plane-strain conditions, inadequate for representing low frequencies and thus are not directly applicable for pseudo-static interaction analyses. The present paper develops analytical solutions to determine transverse and rotation Winkler springs for structures subjected shear waves. The proposed springs overcome the drawbacks of plane-strain models and can be construed as a generalisation of them. The springs are obtained as a function of the seismic wavelength and the ground-structure stiffness contrast. Results obtained are validated against solutions from the literature and numerical results from a full 3D finite-element model. A non-dimensional parametric study is also presented, that allow an expedited evaluation of ground springs for practical applications.

Developing design response spectra for Benghazi city including soil magnification effects

Earthquakes in some countries worldwide cause loss of lives and properties. In Libya, design of structures to resist earthquake forces was based on a paper work published by Mallick in 1976. In 1977, and after small changes in the earthquake zoning map of Libya, the ministry of housing, at that time, adopted it as a draft code of practice for the design of structures to resist earthquake forces. With Benghazi city undergoing significant rehabilitation and development programs, including major national projects, there is a pressing need to estimate updated probabilistic seismic hazard maps and design response spectra specific to the city. This study presented updated probabilistic seismic hazard ground motion for Benghazi city, considering different return periods and accounting for peak ground acceleration (PGA) values with various soil conditions. Proposed design response spectra for Benghazi City unveils substantial PGA amplifications in soft soil areas.

Numerical modelling of impact seismic sources using the stress glut theory

SUMMARY Meteorite impacts have proved to be a significant source of seismic signal on the Moon, and have now been recorded on Mars by InSight seismometers. Understanding how impacts produce seismic signal is key to the interpretation of this unique data, and to improve their identification in continuous seismic records. Here, we use the seismic Representation Theorem, and particularly the stress glut theory, to model the seismic motion resulting from impact cratering. The source is described by equivalent forces, some resulting from the impactor momentum transfer, and others from the stress glut, which represents the mechanical effect of plasticity and non linear processes in the source region. We condense these equivalent forces into a point-source with a time-varying single force and nine-component moment tensor. This analytical representation bridges the gap between the complex dynamics of crater formation, and the linear point-source representation classically used in seismology. Using the multiphysics modelling software HOSS, we develop a method to compute the stress glut of an impact, and the associated point-source from hypervelocity impact simulations. For a vertical and an oblique impact at 1000 m s−1, we show that the moment tensor presents a significant deviatoric component. Hence, the source is not an ideal isotropic explosion contrary to previous assumptions, and draws closer to a double couple for the oblique impact. The contribution of the point force to the seismic signal appears negligible. We verify this model by comparing two signals: (1) HOSS is coupled to SPECFEM3D to propagate the near-source signal elastically to remote seismic stations; (2) the point-source model derived from the stress-glut theory is used to generate displacements at the same distance. The comparison shows that the point-source model is accurately simulating the low-frequency impact seismic waveform, and its seismic moment is in trend with Lunar and Martian impact data. High-frequencies discrepancies exist, which are partly related to finite-source effects, but might be further explained by the difference in mathematical framework between classical seismology and HOSS’ numerical modelling.

Seismic resilience: Innovations in structural engineering for earthquake-prone areas

Abstract Objective The contemporary structural engineering notion of "seismic resilience" is to yield a public to its pre-earthquake state in precise time. The goal of our research is the OMRF (Ordinary Moment Resisting Frame), which is mid-rise building that had exposed to several earthquakes. The research examined the constructions mechanical act and seismic confrontation. Methods For hazard evaluations, the building's proneness and functionality were measured using the Seismic Resilience Index (SRI) and delicateness tops. The course had five critical phases: Selecting the goal buildings, picking and ascending a assembly of repetitive seismic ground motion (SGM) archives, emerging brittleness outsides, manufacture incremental dynamic analysis (IDA), and scheming the functionality curve using the seismic resilience index (RI) are the other steps. Findings: It was evident from the hazard evaluation, which included IDA and flimsiness surface examination, that the nominated assemblies required the structural integrity needed to endure a 15-second repeated earthquake. Novelty As the probable seismic ground acceleration raised, it was also probable to figure the variation in functionality, SRI, resilience, structural losses, and the amount of time desirable to mend from numerous presentation stages. These outcomes highpoint the worth of cutting-edge structural engineering methods for educating buildings' seismic resilience in earthquake-prone districts. More resilient configurations that can better endure and recover from seismic shocks can be attained by addressing strategy errors and improving structural presentation.

Seismic risk model for regional buildings that considers the influence of temperature and intensity measures

Probability-based seismic risk and random ground motion theories are used to develop high-precision seismic risk functions and models of regional building portfolios. Seismic intensity measures and structural demand parameters are critical in developing seismic fragility models for building clusters. Earthquake damage data indicate that ambient temperature may influence structural fragility, but the effect of temperature on seismic risk is commonly ignored. A method for quantifying seismic intensity is proposed in this study and optimized using 450,000 acceleration records detected by 15 stations during the Luding earthquake in Sichuan, China, on September 5, 2022. Considering the influence of different temperature fields, the proposed seismic intensity quantification method was used to evaluate structural damage during two destructive earthquakes in Xinjiang, China (the Aktao earthquake on September 2, 2003 and the Zhaosu earthquake on December 1, 2003). Three structural empirical seismic risk matrices and cloud maps were developed, and a vulnerability membership curve was generated using fuzzy decision-making, a membership algorithm, and nonlinear regression. A structural seismic risk index function was proposed, and vulnerability cloud diagrams and curves based on actual structural seismic vulnerability datasets were generated. Together, these analyses demonstrated the impact of temperature and intensity measures on the seismic risk to regional buildings. A model was established to evaluate the seismic risk to typical structures and was validated using historical earthquake damage data from China. This method improves the accuracy of earthquake risk and vulnerability assessments for typical building clusters given the influence of temperature. The developed seismic risk prediction model enables the evaluation of vulnerability and resilience of typical structures under low and high temperatures.