Propagation of seismic waves through liquefied soils

To predict the earthquake response of saturated porous media it is essential to correctly simulate the generation, redistribution, and dissipation of excess pore water pressure during and after earthquake shaking. To this end, a reliable numerical tool requires a dynamic, fully coupled formulation for solid–fluid interaction and a versatile constitutive model. Presented in this paper is a 3D finite element framework that has been developed and utilized for this purpose. The framework employs fully coupled dynamic field equations with a u–p–U formulation for simulation of pore fluid and solid skeleton interaction and a SANISAND constitutive model for response of solid skeleton. After a detailed verification and validation of the formulation and implementation of the developed numerical tool, it is employed in the seismic response of saturated porous media. The study includes examination of the mechanism of propagation of the earthquake-induced shear waves and liquefaction phenomenon in uniform and layered profiles of saturated sand deposits.     

Predicting seismic permanent displacement of soil walls under surcharge based on limit analysis approach

Seismic permanent displacement of the soil walls plays an important role in design of these structures. Due to the increase in growth of urban areas and the limitations in use of flat grounds, many structures are built near slopes and retaining walls. During earthquakes, these structures can apply an additional surcharge on the wall. The intensity and location of the surcharge is of considerable importance on the seismic displacements of the soil wall. In this study, by using the limit analysis and upper bound theorem, seismic permanent displacement of the soil wall under surcharge has been analyzed. Thus, a formulation is presented for calculating the yield acceleration and seismic displacement for different surcharge conditions. The effect of seismic acceleration, surcharge intensity, its location and soil properties is investigated. A parameter called the “displacement coefficient” is proposed, and is a potential modification for Newmark’s sliding-block method.     

Nonlinear modeling of cyclic response of RC beam–column joints reinforced by plain bars

Experience of previous earthquakes shows that a considerable portion of buildings reinforced with plain bars sustain relatively large damages especially at the beam–column joints where the damages are mostly caused by either diagonal shear cracks or intersectional cracks caused by bar slippage. While previous works mainly focus on shear failure mode, in this study, the emphasis is placed on slip based cracks as the dominant failure mode. A systematic procedure is introduced to predict the dominant failure mode at the joint which is based on the dimensional properties, reinforcement details, and axial and shear load at the joint. In addition, a relatively simple and efficient nonlinear model is proposed to simulate pre- and post-elastic behavior of the joints which fail under bar slippage mode. In this model, beam and column components are represented by linear elastic elements, dimensions of the joint panel are defined by rigid elements, and effect of slip is taken into account by a nonlinear rotational spring at the end of the beam. The proposed method is validated by experimental results for both internal and external joints .     

Modeling cyclic shearing of sands in the semifluidized state

Liquefaction is associated with the loss of mean effective stress and increase of the pore water pressure in saturated granular materials due to their contractive tendency under cyclic shear loading. The loss of mean effective stress is linked to loss of grain contacts, bringing the granular material to a “semifluidized state” and leading to development and accumulation of large cyclic shear strains. Constitutive modeling of the cyclic stress-strain response in earthquake-induced liquefaction and post-liquefaction is complex and yet very important for stress-deformation and performance-based analysis of sand deposits. A new state internal variable named strain liquefaction factor is introduced that evolves at low mean effective stresses, and its constitutive role is to reduce the plastic shear stiffness and dilatancy while maintaining the same plastic volumetric strain rate in the semifluidized state. This new constitutive ingredient is added to an existing critical state compatible, bounding surface plasticity reference model, that is well established for constitutive modeling of cyclic response of sands in the pre-liquefaction state. The roles of the key components of the proposed formulation are examined in a series of sensitivity analyses. Their combined effects in improving the performance of the reference model are examined by simulating undrained cyclic simple shear tests on Ottawa sand, with focus on reproducing the increasing shear strain amplitude as well as its saturation in the post-liquefaction response.     

A GIS-based time-dependent seismic source modeling of Northern Iran

The first step in any seismic hazard study is the definition of seismogenic sources and the estimation of magnitude-frequency relationships for each source. There is as yet no standard methodology for source modeling and many researchers have worked on this topic. This study is an effort to define linear and area seismic sources for Northern Iran. The linear or fault sources are developed based on tectonic features and characteristic earthquakes while the area sources are developed based on spatial distribution of small to moderate earthquakes. Time-dependent recurrence relationships are developed for fault sources using renewal approach while time-independent frequency-magnitude relationships are proposed for area sources based on Poisson process. GIS functionalities are used in this study to introduce and incorporate spatial-temporal and geostatistical indices in delineating area seismic sources. The proposed methodology is used to model seismic sources for an area of about 500 by 400 square kilometers around Tehran. Previous researches and reports are studied to compile an earthquake/fault catalog that is as complete as possible. All events are transformed to uniform magnitude scale; duplicate events and dependent shocks are removed. Completeness and time distribution of the compiled catalog is taken into account. The proposed area and linear seismic sources in conjunction with defined recurrence relationships can be used to develop time-dependent probabilistic seismic hazard analysis of Northern Iran.     

ORELM: A Novel Machine Learning Approach for Prediction of Flyrock in Mine Blasting

Blast-induced flyrock is a hazardous and undesirable phenomenon that may occur in surface mines, especially when blasting takes place near residential areas. Therefore, accurate prediction of flyrock distance is of high significance in the determination of the statutory danger area. To this end, there is a practical need to propose an accurate model to predict flyrock. Aiming at this topic, this study presents two machine learning models, including extreme learning machine (ELM) and outlier robust ELM (ORELM), for predicting flyrock. To the best of our knowledge, this is the first work that investigates the use of ORELM model in the field of flyrock prediction. To construct and verify the proposed ELM and ORELM models, a database including 82 datasets has been collected from the three granite quarry sites in Malaysia. Additionally, artificial neural network (ANN) and multiple regression models were used for comparison. According to the results, both ELM and ORELM models performed satisfactorily, and their performances were far better compared to the performances of ANN and multiple regression models.

     

Estimating trends of the Mediterranean Sea level changes from tide gauge and satellite altimetry data (1993-2015)

The impact of climate change on sea level has received a great deal of attention by scientists worldwide.In this context,the problem of sea levels on global and regional scales have been analyzed in a number of studies based on tide gauges observations and satellite altimetry measurements.This study focuses on trend estimates from 18 high-quality tide gauge stations along the Mediterranean Sea coast.The seasonal Mann-Kendall test was run at a 5% significance level for each of the 18 stations for the period of 1993-2015(satellite altimetry era).The results of this test indicate that the trends for 17 stations were statistically significant and showed an increase(no significant trend was observed only at one station).The rates of sea level change for the 17 stations that exhibit significant trends,estimated using seasonal Sen's approach,range after correction for Vertical Land Motion(VLM)from 1.48 to 8.72 mm/a for the period1993-201 5.Furthermore,the magnitude of change at the location of each tide gauge station was estimated using the satellite altimetry measurements.Thus,the results obtained agree with those from the tide-gauge data analysis.     

Estimating trends of the Mediterranean Sea level changes from tide gauge and satellite altimetry data (1993-2015)

The impact of climate change on sea level has received a great deal of attention by scientists worldwide. In this context, the problem of sea levels on global and regional scales have been analyzed in a number of studies based on tide gauges observations and satellite altimetry measurements. This study focuses on trend estimates from 18 high-quality tide gauge stations along the Mediterranean Sea coast. The seasonal Mann-Kendall test was run at a 5% significance level for each of the 18 stations for the period of 1993-2015 (satellite altimetry era). The results of this test indicate that the trends for 17 stations were statistically significant and showed an increase (no significant trend was observed only at one station). The rates of sea level change for the 17 stations that exhibit significant trends, estimated using seasonal Sen's approach, range after correction for Vertical Land Motion (VLM) from 1.48 to 8.72 mm/a for the period 1993-2015. Furthermore, the magnitude of change at the location of each tide gauge station was estimated using the satellite altimetry measurements. Thus, the results obtained agree with those from the tide-gauge data analysis.