---

Due to its critical impact and significant destructive nature during and after seismic events, soil liquefaction and liquefactioninduced

lateral ground spreading have been increasingly important topics in the geotechnical earthquake engineering field

during the past four decades. The aim of this research is to develop an empirical model for the assessment of liquefaction-induced

lateral ground spreading. This study includes three main stages: compilation of liquefaction-induced lateral ground spreading

data from available earthquake case histories (the total number of 525 data points), detecting importance level of seismological,

topographical and geotechnical parameters for the resulted deformations, and proposing an empirical relation to predict

horizontal ground displacement in both ground slope and free face conditions. The statistical parameters and parametric study

presented for this model indicate the superiority of the current relation over the already introduced relations and its applicability

for engineers.

---

This study presents a series of physical model tests and numerical simulations using PFC2D (both with a dip slip angle=60° and

a soil bed thickness of 0.2 m in model scale)at the acceleration conditions of 1g, 40g, and 80 g to model reverse faulting. The soil

deposits in prototype scale have thicknesses of 0.2 m, 8 m, and 16 m, respectively. This study also investigates the evolution of a

surface deformation profile and the propagation of subsurface rupture traces through overlying sand. This study proposes a

methodology for calibrating the micromechanical material parameters used in the numerical simulation based on the measured

surface settlements of the tested sand bed in the self-weight consolidation stage. The test results show that steeper surface slope

on the surface deformation profile, a wider shear band on the major faulting-induced distortion zone, and more faulting appeared

in the shallower depths in the 1-g reverse faulting model test than in the tests involving higher-g levels. The surface deformation

profile measured from the higher-g physical modeling and that calculated from numerical modeling show good agreement. The

width of the shear band obtained from the numerical simulation was slightly wider than that from the physical modeling at the

same g-levels and the position of the shear band moved an offset of 15 mm in model scale to the footwall compared with the results

of physical modeling.

---

In this paper, an advanced formulation of time-domain, two-dimensional Boundary Element Method (BEM) with material damping is presented. Full space two-dimensional visco-elastodynamic time-convoluted kernels are proposed in order to incorporate proportional damping. This approach is applied to carry out site response analysis of viscoelastic topographic structures subjected to SV and P incident waves. Seismic responses of horizontally layered site, semi-circular canyons, slope topography and ridge sections subjected to these incident waves are analyzed in order to demonstrate the accuracy of the kernels and the applicability of the presented viscoelastic boundary element algorithm. The results show an excellent agreement with recent published results obtained in frequency domain. Also, the effects of different material damping ratios on site response are investigated.

---

This paper presents a simple solution based on the limit equilibrium of sliding soil wedge of reinforced backfill subjected to the horizontal acceleration in the framework of the pseudo-static method. In particular, contrary to most studies on the reinforced retaining wall, the solution proposed in this study, takes into account the effect of the uniform surcharge on the reinforced backfill soil and of its distance from the face of the wall. The results are investigated in dimensionless form of the maximum reinforcement required strength (Kmax), the dimension of the sliding wedge (Lc/H), and the factor of safety (FS). Compared to the reinforced backfill without surcharge, the presence of surcharge over the reinforced backfill and of its distance from the top of the backfill have significant effects on the stability of the system. For a fixed surcharge, a minimum distance of surcharge exists for which the presence of the surcharge does not affect the solution and the failure mechanism is that corresponding to the case of system with no surcharge. A detailed design example is included to illustrate usage of proposed procedures. Comparisons of the present results with available results show a favorable agreement.

---

This paper presents time domain fundamental solutions for the extended Biot's dynamic formulations of two-dimensional (2D) unsaturated poroelasticity. Unsaturated porous media is considered as a porous media in which the voids are saturated with two immiscible fluids, i.e. liquid and gas. At first, the corresponding explicit Laplace transform domain fundamental solution is obtained in terms of skeleton displacements, as well as liquid and gas pressures. Subsequently, the closed-form time domain fundamental solutions are derived by analytical inversion of the Laplace transform domain solutions. Finally, a set of numerical results are presented which verifies the accuracy of the analytically inversed transient fundamental solution and demonstrates some salient features of the elastic waves in unsaturated media..

---

In this paper, normalized displacement amplitude of the ground surface was presented in the presence of the semi-sine shaped valley above the truncated circular cavity embedded in a homogenous isotopic linear elastic half-plane, subjected to obliquely propagating incident SH waves as Ricker wavelet type. The proposed direct time-domain half-plane boundary element formulation was used and extended to analyze the combined multi-boundary topographic problems. While using it, only boundary of the valley and the surrounding cavity should be discretized. The effect of four geometric parameters including shape ratio of the valley, depth ratio, horizontal location ratio and truncation thickness of the cavity and incident wave angle was investigated on the responses at a single dimensionless frequency. The studies showed that surface behavior was completely different due to complex topographic features, compared with the presence of either valley or cavity alone. In addition, the cavity existence below the surface could play a seismic isolation role in the case of vertical incident waves and vice versa for oblique waves.

---

In this paper, an advanced formulation of the spectral finite element method (SFEM) is presented and applied in order to carry out site response analysis of 2D topographic structures subjected to vertically propagating incident in-plane waves in time-domain. The accuracy, efficiency and applicability of the formulation are demonstrated by solving some wave scattering examples. A numerical parametric study has been carried out to study the seismic response of rectangular alluvial valleys subjected to vertically propagating incident SV waves. It is shown that the amplification pattern of the valley and its frequency characteristics depend strongly on its shape ratio. The natural frequency of the rectangular alluvial valley decreases as the shape ratio of the valley decreases. The maximum amplification ratio along the ground surface occurs at the center of the valley. A simple formula has been proposed for making initial estimation of the natural period of the valley in site effect microzonation studies.

---

In this paper, different aspects of the behavior of 2×2 pile groups under liquefaction-induced lateral spreading in a 3-layer soil profile is investigated using large scale 1-g shake table test. Different parameters of the response of soil and piles including time-histories of accelerations, pore water pressures, displacements and bending moments are presented and discussed in the paper. In addition, distribution of lateral forces due to lateral spreading on individual piles of the groups is investigated in detail. The results show that total lateral forces on the piles are influenced by the shadow effect as well as the superstructure mass attached to the pile cap. It was also found that lateral forces exerted on the piles in the lower half of the liquefied layer are significantly larger than those recommended by the design code. Based on the numerical analyses performed, it is shown that the displacement based method is more capable of predicting the pile group behavior in this experiment comparing to the force based method provided that the model parameters are tuned.

---

In this paper we are going to show the importance of mode identification in microtremor array analysis. The idea come from four concentric ambient noise array recordings with aperture 100 to 1000 meters, performed in southern urban area of Tehran near the shrine of Imam Khomeini. These measurements were part of a comprehensive research project with the aim of determination of deep shear wave velocity model of Tehran alluvial deposits. Using appropriate signal processing techniques, including array processing methods as well as classical and time-frequency horizontal/vertical spectral ratio, the dispersion curves of surface waves, fundamental resonance frequency and Ellipticity of Rayleigh waves, were extracted. In the final step, the shear wave velocity profile of the site was determined by joint inversion of all of these attributes. The results show 2 different energetic trends in dispersion curves, for arrays of aperture 200 and 400 meters that one of them is coincide with 100m aperture array. For array with aperture 1000m any clear trend of energy could be observed because of deficiency of energy in low frequency. The inversion of data obtained by 100m aperture array alone, assuming the dispersion curve as fundamental mode (a common procedure in urban area) result in shear wave velocity that is not match with existing geological information. Performing the inversion, assuming 2 energetic trends, observed for larger arrays one as fundamental mode and another as mode 1 of Rayleigh waves, can modify significantly the shear wave velocity profile in accordance with existing geological and geotechnical information. This study show the importance of extracting of correct dispersion curves with detecting fundamental and higher modes, using array measurement with various aperture at one place to obtain more realistic shear wave velocity profile.

---

To increase the safety of structures against strong ground motions and their life due to environmental issues on the earth and saving in terms of materials, it is necessary to expand and upgrade seismic resistant systems. However, more cost-effective systems which have sufficient influence on the seismic performance of structures and also more compatibility with the regional conditions, will be more desirable than other systems. One of the seismic resistance systems is seismic isolation. In the event of interest in using the seismic isolation system for a mounted building on piles, the costly construction of piles and isolation equipment shall be provided simultaneously. The seismic isolating using sleeved-piles which is generally used in combination with various damper systems, can help to overcome this issue. In this research a seismic isolator system using sleeved-pile has been studied while considering the damping behavior of the soil-rubber mixture as the only source of damping. To investigate the proposed system, a series of tests including static lateral load test, dynamic free and forced vibration tests, were performed on a model pile in a field laboratory which has been constructed for this purpose. According to results of tests the proposed system has a good deformation ability and damping characteristics, and as a method of seismic isolation is completely efficient.

---

The interactive effects of adjacent buildings on their seismic performance are not frequently considered in seismic design. The adjacent buildings, however, are interrelated through the soil during seismic ground motions. The seismic energy is redistributed in the neighboring buildings through multiple structure-soil-structure interactions (SSSI). For example, in an area congested with many nearby tall and/or heavy buildings, accounting for the proximity effects of the adjacent buildings is very important. To solve the problem of SSSI successfully, researchers indicate two main research areas where need the most attention: 1) accounting for soil nonlinearity in an efficient way, and 2) spatial analysis of full 3D soil-structure models. In the present study, three-dimensional finite element models of tall buildings on different flexible foundation soils are used to evaluate the extent of cross interaction of adjacent buildings. Soil nonlinearity under cyclic loading is accounted for by Equivalent Linear Method (ELM) as to conduct large parametric studies in the field of seismic soil-structure interaction, the application of ELM is preferred over other alternatives (such as application of complicated constitutive soil models) due to the efficiency and reliability of its results. 15 and 30 story steel structures with pile foundations on two sandy and clayey sites are designed according to modern codes and then subjected to several actual earthquake records scaled to represent the seismicity of the building sites. Results show the cross interaction of adjacent buildings on flexible soils, depending on their proximity, increases dynamic displacements of buildings and reduces their base shears. 

---

Experimental study of the cyclic behavior of unsaturated materials is more complex than that of the saturated materials due to the required equipment, experience and time. Furthering investigations in the field of unsaturated materials is necessary to better understand its complexity and sensitivity of unsaturated cyclic parameters to different determinants such as suction path, stress path, loading speed, deviatoric stress amplitude, physical specifications, and etc. To this end, the main focus of this study has been to analyze the effects of factors such as mean net stress and deviatoric stress levels in fast cyclic loading on the cyclic behavior of a normally consolidated unsaturated fine-grained trade soil, namely the Zenoz kaolin. Various unsaturated tests were performed in three mean net stress levels and three amplitudes of cyclic deviatoric stress levels. Results showed that increase of suction in the same strain level leads to increase in stiffness in normally consolidated samples (i.e. increase in elastic modulus and shear modulus and decrease in damping ratio). Also, in the same suction value and strain level, increase of the mean net stress during the isotropic consolidation causes to the denser normally consolidated samples and results to increase of elastic modulus and shear modulus, and decrease of damping ratio.

---

In this study a series of cyclic triaxial tests were performed to examine the undrained dynamic resistance of silty sand reinforced with various arrangements of geotextile layers. The silt content of samples varies in percentage from 0, 10, 20, 30, 40 and 50%. A total of 32 laboratory cyclic triaxial tests have been performed on silty sand samples reinforced with geotextile layers in different depths. All tests were performed with 100 kPa confining pressure, subjected to an isotropic consolidated undrained (CIU) condition. The tests were conducted at a frequency of 2 Hz. Results indicate that both the geotextile arrangement and the silt content were most essential in the liquefaction potential of reinforced sands. An increase in the number of geotextile layers enhanced the cyclic resistance of reinforced samples against the liquefaction potential. It was also found that when the geotextile layer was posited near the top of the specimen (load application part) the liquefaction resistance would increase (e.g. for clean sands, the improvement of liquefaction resistance caused by the geotextile layer had a 0.2 depth, and the sample height was 5.5 times greater than the geotextile layer inserted in mid height of sample H). Based on the obtained results, effects of geotextile on liquefaction resistance decreased as fines content increased to about 33%. Further increase in the fines content however, would lead to higher in reinforcement advantages. The liquefaction improvement is more effective with a higher number of geotextile layers. The results also revealed that the reinforcement effect in FC≈33 % is at its lowest amount. 

---

Currently, there is no reliable design procedure which considers all aspects of liquefaction effects on shallow foundations. There are many light and heavy structures resting on saturated sand with high liquefaction potential in seismic areas. The aim of this experimental and numerical study is to evaluate the performance of two shallow foundations with different contact pressures in liquefaction. The results of the centrifuge experiment of shallow foundations with surcharges of three-story and nine-story buildings on liquefiable sand are presented in detail. Although entire soil profile liquefied, no liquefaction observed under the foundations. There was a clear difference in settlement mechanisms observed beneath the shallow foundation and in the free-field. The heavy foundation fluctuated more strongly compared with the lighter one. The effect of soil permeability and contact pressure on foundation response was investigated during numerical study. The experiment was simulated two dimensionally using a fully coupled nonlinear constitutive model (UBCSAND) implemented in a finite difference program, FLAC-2D. The results show that settlement of foundations increased with the increase of soil permeability. Trends of excess pore water pressure were captured reasonably by the soil model, but the settlement mechanisms were different. The soil model underestimated total liquefaction-induced settlement of shallow foundation, especially for light foundation.

---

Complete recognition of calcareous sediments engineering behavior considering their local expansion and wide variety of engineering properties is very important. In south parts of Iran, there are some carbonate hydrocarbon reservoirs which are covered by calcareous deposits. Hormuz Island in is one of the most strategic areas in Hormuz Strait between Persian Gulf and Oman Sea. In this study, a series of undrained monotonic and cyclic simple shear tests was performed on saturated Hormuz calcareous sand specimens using hollow cylinder torsional apparatus. The tests were carried out on specimens with various relative densities under different effective consolidation stresses. Based on the results, pore pressure generation, shear strain development, stress-strain characteristics of the specimens are presented and compared with the technical literature. In addition, dissipation of strain-based energy during the cyclic loading and its relation to excess pore water pressure is described. The cyclic resistance curves of specimens with different initial conditions are plotted. Also the results of monotonic and cyclic tests are compared together for better interpretation of Hormuz calcareous sand under undrained torsional loading.

---

Investigating the seismic amplification of incident waves induced by subsurface cavities and characterizing its patterns for the ground surface are important in seismology, geophysics and earthquake engineering both in theory and practical application. Nowadays, it has been established that the seismic ground response above subsurface structures can be different from the free-field motion during earthquakes. In this regard, this research studied preliminary results of a numerical parametric study on the seismic response of the ground surface above subsurface cavity. Basically, this study is applied to get new idea to move a step forward in site response analysis which can be used in the seismic microzonation of areas located above underground spaces. For analysis purposes, a numerical time-domain analysis is performed by utilizing a robust numerical algorithm working based on the boundary element method. It is observed that the amplification of the ground surface underlain by a shallow elliptical cavity is increased in long periods. Some preliminary simple relationships and tables are presented which could be used while introducing simple preliminary ideas for modification of the standard design spectra in building codes and seismic microzonation studies.