Journal of Engineering Geology

---

Introduction
The prediction of landslide occurrence in a region is very important in reducing the risks and damages caused by this.landslide as a natural disaster in Iran caused a lot of life and financial losses to Iran annually. According to the National Committee on Natural Disaster Reduction of the Ministry of the Interior in 1994, the share of annual damage caused by mass movements in Iran is estimated at 500 billion rials. In the meantime Kurdistan province is the third province in terms of landslide phenomenon after Mazandaran and Golestan. If considering the area is at a higher level. The city of Bijar in this province has a high potential for a wide range of landslides with a combination of mainly mountain topographical factors, lithologic conditions and positioning between two major faults in the region. In this research, using quantitative methods and models on the quantitative  factors of this phenomenon based on the level of information given by past mass movements and influential factors, focusing on artificial neural network method, susceptibility zones were determined by determining the possible risk level.
Knowing such natural events requires proper management of the risks posed by them. On the other hand, artificial neural network as a quantitative model is capable of learning, generalization and decision making, and less need to analyze the accuracy of data in comparison to statistical methods. Map of the susceptibility of the areas to the landslide is an important tool for landuse planning. However, there are many issues in the formation of this phenomenon, which, due to the complexity of the natural processes arising from the relationship between the outcome (dependent variable) and the factors (independent variables), puts into question the general zoning of such areas.
Methodology
Bijar is located in the northeastern part of Kurdistan province, along the longitude of   47 ' 29° to 47 ° 47' east, in latitude 35 ° 35 'to 35' 59 °north. In recent years, the development of the Geographic Information System (GIS) and spatial analysis techniques have improved the risk of indirect zoning. In this regard, artificial neural networks can cover a significant part of these needs.Implementing the neural network model requires learning data. Without learning data, it's virtually impossible to make neural networks. In this paper, learning data shows the occurrence of landslides which have geographical coordinates and were obtained from the Kurdistan Province Natural Resources Organization. In general, learning data in GIS and remote sensing can include data or raster, which in this paper is a point phenomenon and has 144 cases.  However, because of the large extent of the study area and the low number of them, as well as the lack of risk of any landslide zone (from low to very high), the points should be classified as well, and, in terms of numbers, Acceptance. Also, the number of points of relative value In terms of numbers, the conditions are the Normal and the same (that is, the appropriate geographical distribution and distribution in each class) would be more accurate; thus, to create a classifiable spectrum of the AHP Was used. It should be noted that all the maps were standardized in the format and format of the Raster in a matrix (698 rows in 897 columns) identical with a size of 30 * 30 meters. This means that each map has 626,106 pixels of varying value and somewhat similar. In addition, the AHP model was used to categorize the studied area from very desirable (hazardous) to very undesirable (very dangerous) areas. Also, 33 points were added to the learning data on different levels of the map derived from the AHP model. But in order to verify accurately the model, only landslide occurrences were considered.
In order to find out the factors of landslide in Bijar, a map of slope, Aspect, elevation, distance from the fault, distance from the road, distance from the river, Drainage density, lithology and land use using ArcGIS software were prepared and digitized.
After compiling and categorizing these variables, at first, each of the effective criteria in the field was divided into six sub-criteria (land suitability for landslide) from very desirable to very undesirable conditions. The present study utilizes the technique of multi-layer propspert neural networks using post-propagation algorithm (BP). In addition to correcting and editing the layers, the neural network model was implemented using the classification method and applying two types of functions (linear and sigmoid). Then, using the test-error method, the study of the magnitude of the error and the period of the repetition and the change in the number of hidden layers and weights, both functions were performed. Finally, the sigmoid function, which yielded a better result, was selected as the proposed and final function.Order to verify the (accuracy) of the map taken with the existing landslide zones, the final map of the neural network model was again transferred to the ArcGIS software. Finally, the available landscapes on the map resulted from the adaptive neural network model, which, by comparison, gave a percentage and amount Accuracy of each class was achieved.
Result
The input layer were calculated to six classes based on the desirability of mass movements. This decision approach reduces the complexity of the network and improves its performance.
For this purpose. The AHP method was used to define non-slip pixels and range classification.
To implement this method, 9 variables discussed, were scaled up to the most suitable and un suitable range. The final weight of these variables was obtained by using the Thomas saati pair comparison (Table 4), the study area was divided into five categories according to the map for land suitability for landslide hazard. From each class, the 20-pixel from AHP model was selected for network learning in a completely randomized manner. The proposed model is an artificial neural network of MLP multi-layered perceptron with levenberg-marquardt learning algorithm. An early stopping method was used to improve network optimization. Several hidden layers were tested to find the best results. It should be noted that in the structure of all networks, at least the optimal design with the middle one is used, but in their structural composition they are also used with mid-duplex networks. In this paper, the use of tow mid-layers showed better results.  In all Simulations have been made, the mean square error index, as a guide, indicates the network performance in learning the existing model. By changing the number of intermediate neurons and changing the weights as try and error, the most appropriate network model was obtained for the purpose. In this study, the structure of the network with 9 input layers, 2 hidden layers, 1500 repetitions in both functions was accepted as the final structure. The main structure of the neural network with two linear and sigmoid functions was prepared with acceptable error, and the study area was analyzed with a total area of ​​564 km2 with 9 input variables converted into raster data to 30 × 30 pixels. From 564 km2 based on the sigmoid function 61.17% and based on the linear function, 72.76% of the area is unsuitable and very unsuitable in the area where expose to high risk. In both networks, there were very few areas in both optimal and moderate classes (Figures 16 and 17), which indicate the high talent of the area for landslide as a threat. Then, ArcGIS software was used to evaluate the efficiency and accuracy of the model. For this purpose, the point of landslide and zoning maps were combined, compared and anlayzed. The results showed in the sigmoid function 75 items of Landslides were in a very unsuitable range, which included 61% of the total of region.
Conclusion
 In the linear function, approximately 69% of the landslides are in a very unsuitable range and the unsuitable results are about 57%, which results in the success of the model designed in the neural networks (MLP). In the end, the network with sigmoid function is negligibly better than the linear function network.The results show that Bijar and its functions are relatively prone to occurrence of landslides, so that nearly 60% of the city's area is a high risk area with a high risk and only 2% is a low-risk region. The hazardous areas are mainly located around the city of Bijar especially southern and southeast. These areas correspond to high altitudes and maximum fault density and lime lithology with marl (Qom Formation). The model can be very challenging, because of innovative nature of the research, that means need more detailed and comprehensive studies../files/site1/files/121/neiri_Abstract.pdf

---

Introduction
The rock block volumes are formed due to the intersection of discontinuities in the jointed rock mass. The block dimensions affected by joint spacing, joint orientation, joint sets, are taken to be the most important parameters determining the rock mass behavior, strength parameters, and deformations. In the numerical modeling using distinct element method, the creation of the discontinuities can affect the final results very much. Using 3DEC software, it is possible to create joint sets in four conditions: regular and persistent, regular and non-persistent, non-regular and persistent, irregular and non-persistent. As an important point to consider, the major effect of block dimensions on rock mass behavior, strength parameters and deformation modulus can help to decide which one is most suitable to indicate the real conditions of rock mass. As explained in the previous studies, the use of persistent joints leads to the block dimensions being considered as small ones. In this way, due to the high strength of intact rock compared to the joints, the possibility of instability increases.
Material and methods
In this research, from quantitative point of view, Geological Strength Index (GSI) is calculated, based on block dimensions as an influential parameter, to consider the most appropriate case for creating joints in the numerical method. In this regard, according to valuable studies in Bakhtiari dam structure, the characteristics of discontinuities system and GSI of rock mass are utilized to come up with real conditions. Then, the modeling is done with different conditions of joints, block volume distribution, GSI for each case, and the results are compared with quantitative ones. And then the most suitable case for creation of joints in numerical modeling is suggested by using 3 DEC software, regarding the blocks volumes, type of distribution function, and GSI value. Also, the accuracy of this finding is investigated for other structures, independently of input parameters, by making changes in spacing, and joints persistence as two effective parameters in rock blocks dimensions. Owing to the difficulty in the accurate definition of joints persistence, which is related to dimensions of the location, the numerical models for joint persistence are done in an acceptable level in order to create blocks with high conformity in terms of the dimensions. Then, the comparison is made between block dimensions from perspectives of numerical models and GSI values, to choose the best ones showing high conformity with real conditions.
Results and discussion
The comparison of the modeling results using creation of joints in different cases with quantitative results obtained according to geological strength shows that the created block volumes are not properly distributed due to the creation of joints as irregular ones in the two cases of persistent and non-persistent. In this case, the blocks volume changes from a few centimeter to cubic meter, and as the block dimensions increase, the created blocks become bigger. Thus, according to the created blocks volume and the obtained GSI range, the creation of joints is not a suitable method as an irregular one. The creation of regular and persistent joints is not an appropriate method either, as the most created blocks are small, and blocks volume distributions do not comply with quantitative distribution. But with creation of joints as regular and persistent ones, the distribution function of blocks volume in numerical method and quantitative method is log normal. Therefore, according to the created blocks dimensions and GSI range using 3DEC software, the most suitable case is the creation of joints as regular and non-persistent ones. 
Conclusion
According to the obtained results in the four cases, when the joints are considered only as regular and non-persistent ones, the blocks volume range is more compatible with real conditions and follows log normal distribution. Thus it can be concluded that the suggested method for creation of joints in the numerical modeling using 3DEC software is more suitable than others considering the rock blocks dimensions and their distributions. This method can be utilized in any structure to accurately define the persistence of joints regarding created blocks dimension.

---

Introduction
Contaminants can be categorized into organic and inorganic groups. Organic contaminants are carbon based, and their presence in waste forms may be as a single contaminant associated with inorganic contaminants, or a suite of complex mixtures which may be toxic at very low concentrations. Organics of greatest environmental concern are usually refined petroleum products, chlorinated and non-chlorinated solvents, manufactured biocides, organic sledges and substances from manufacturing processes. Most contamination due to organics are associated with accidental spills and leaks, originating from equipment cleaning, maintenance, storage tanks, residue from used containers and outdated material (Yong and Mulligan, 2004). Transport and fate of organic contaminants are important. Organic contamination migrations are due to advection (by fluid flow through soil) and diffusion, but other forms of transport e.g. infiltration may also contribute to migration (Environment Agency, 2002). The response of the soil to a contaminant depends upon the type of soil and the nature of the contaminant. The sensitivity of soil to contaminants depends upon the type of soil (such as particle size, mineral structure, bonding characteristics between particles and ion exchange capacity) and the nature of contaminants. Fang (1997) defined a sensitivity index (ranging from 0 to 1) to different types of soil. Sensitivity of sand and gravel (0.01 to 0.1) is much lower than clay particles (0.6-0.9). There are a number of techniques for remediation of contaminated land. These include physical (washing, flushing, thermal, vacuum extraction, solvent extraction), chemical (stabilization and solidification) and bioremediation techniques. However, the applicability and feasibility of different methods for remediation are dependent on many factors such as soil characteristics (soil type, degree of compaction and saturation), site geology, depth of contamination, extent of contaminant in lateral direction, topography, surface and ground water and the type and amount of contaminants. Thermal treatment and using surfactants are the most popular methods for remediating the soil contaminated with petroleum compounds. In this research remediation of a soil contaminated with different percentages of gasoline was studied through physical techniques in laboratory. The applied physical techniques were thermal technique and use of two different kinds of surfactants. The obtained results were compared with each other and discussion was performed.
Material and methods
Soil, gasoline and surfactant are the basic materials that were used in this work. The soil that was used in this testing program was a clayey soil. Two different types of ionic and nonionic surfactant, namely Tween 80 and SDS, were used in this work for remediating soil, contaminated with gasoline. Contaminated soil was prepared by adding 5 and 10 % weight (to air dried soil) of gasoline. 6 kg air dried soil was selected and the desired amount of gasoline was weighted, then it was sprayed on the soil and thoroughly mixed by hand for about 2 hours. The prepared mixture was kept inside a covered container for a week in order to come to equilibrium with the soil. For thermal remediation the contaminated soil with a specific percent of gasoline was kept inside a constant convection oven at 50, 100, and 150^oC for about one week to desorb the contaminating compound. Tween 80 and SDS were used for remediation of the contaminated soil. The amount of used Tween 80 was 25% weight of contaminating compound and selection of SDS amount was based on 50% weight of contaminating matter. The samples for the main tests were prepared by static compaction according to the optimum water content and maximum dry unit weight that were obtained from standard compaction tests. Atterberg limits, grain size distribution, compaction and unconfined compression tests were performed on samples of natural, contaminated and remediated soil according to the ASTM standard.
Results and discussion
The results of Atterberg limits (LL, PL and PI) for the contaminated soil (with 5 and 10 % gasoline) indicated that the values of them are increased with increasing the percent of gasoline. These values are nearly the same as natural soil after remediation with thermal method and surfactants. The grain size distribution curves were determined for the natural soil, contaminated soil with 5% and10% gasoline and soil remediated by thermal and surfactant techniques. The results showed that by using thermal technique the percent of clay is decreased and the percent of sand and particularly silt is increased by increasing temperature. The results of grain size distribution for the soils remediated by SDS and Tween 80 showed that the percent of clay is reduced but the percent of silt and sand are increased. Comparing the results of the two surfactants shows that the effect of Tween 80 in reduction of the percent of clay is more than SDS. The results showed that after thermal treatment, the maximum dry unit weight decreases and the optimum water content increases. For the contaminated soil with gasoline a reduction in maximum dry unit weight is observed compared with natural soil. The effect of SDS and Tween 80 on soil remediation is reduction in maximum dry unit weight and optimum water content. The results of compression strength showed that adding gasoline to soil causes a reduction in final strength and this reduction is a function of gasoline percent. The results also indicated that the strength of remediated soil by thermal or surfactant techniques, is reached nearly to the strength of natural soil. Scanning electron microscopy (SEM) tests were performed on the samples in order to observe the microstructure of the samples in different conditions (natural and contaminated with different percent of gasoline). The results of SEM showed that the structure of soil is changed by contamination to gasoline. It can be said that the gasoline causes reduction in the thickness of DDL because of low dielectric constant and hence a flocculated structure is formed. In the flocculated structure due to attractive forces, the fine particles paste to each other and form coarse particles. Therefore, variations in the Atterberg limits and compaction parameters can be resulted from forming new structure by adding gasoline. These results of compression strength are not in agreement with the theory of diffuse double layer (DDL). The reduction in dielectric constant would cause a flocculated structure in soil and the strength of the contaminated soil should be increased in comparison with the natural soil. It can be said the viscosity of gasoline cause reduction in the strength of contaminated soil.
Conclusion
In this experimental work a cohesive soil was contaminated with 5% and 10% of gasoline. The experimental tests showed that the properties of contaminated soil are different from natural soil and the change in the properties is a function of gasoline percent. The contaminated soil, was remediated by thermal treatment and also using two surfactants. The results also showed that using surfactants is more effective than using thermal method in soil remediation, and can treat the soil nearly to its original condition.
-Base on the SEM analysis results, adding gasoline to the soil, will change the soil micro structure to a flocculated one.  
-The gradation curves show that adding gasoline to the soil will change the gradation from finer to coarser.
- Contamination to gasoline will change the compaction parameters of the soil, and will reduce the soil final strength significantly.
- The results show that using thermal method and surfactants is effective in remediating the soil, but it is more effective to use surfactants. 
References
Yong, R.N., Mulligan,. “Natural attenuation of the contaminants in soil”, CRC press, Boca Raton, FL (2004).
Fang, M.Y. “Introduction to Environmental Geotechnology”, CRC Press,FL.USA, (1997).

---

Introduction
When loose sand is subjected to seismic shaking, it tends to volume reduction and settlement. The density of the under layers is revealed in the settlement of the ground surface that causes the destruction of the structures located on ground surface. In dry sand layer, settlement in severe shakings occurs under a constant and effective stress condition and very rapid stress. In this regard, the sand deposit settlement is completed before the end of an earthquake, but if the sandy soil layer is saturated and drainage is limited the condition is prepared of fixed volume situation and the major effect of the seismic shocks is generation of exceed pore water pressure. Therefore, the deposit settlement of saturated sand requires a longer time, varying from a few minutes to a few days, depending on the permeability and compressibility of the soil and the length of the drainage path. The main purpose of the present study is to evaluate liquefaction hazards along Tabriz Metro Line 2 with using Standard Penetration test (SPT) and shear wave velocity (Vs) methods. Also, the probable rate of settlement in the soil layers in study area and correlation with liquefaction potential index (LPI) according to both procedure have been determined and discussed in the following paragraphs.
Material and Methods
In order to evaluate the liquefaction potential of soils using two field methods, geotechnical information of 54 boreholes in Tabriz Metro Line 2 were collected. The types of soil and geotechnical properties can affect the liquefaction potential. In this study, the gravely sand, silty sand, silty and sandy soils were studied. Ground water level is one of the main parameters in in soil liquefaction potential evaluation of soils. Variation of water level in boreholes is 2 to 30 meters. The peak ground acceleration (PGA) is necessary for the analysis of boreholes to evaluate liquefaction potential of soils. PGA values were selected in each boreholes position according to the Iranian Code of Practice for Seismic Resistant Design of Buildings (Code-2800-ver.4) equal to 0.35g (for return period 475 years and design life 50 years). Liquefaction potential of soil layers based on SPT results with appliying Idriss and Boulanger (2010) method has been assessed. Andrus, Stokoe and Jung (2004) procedure was used in shear wave velocity (Vs) method (with assuming cementation and un cementation condition in soils). Liquefaction potential index (LPI) of soil layers was calculated for both field tests results. Then, probable rate of settlement due to liquefaction in saturate soil layers was determined. Tokimatsu and Seed (1978) method applied for SPT results, Yoshimine (1992), Yoshimine et al. (2006) and Yi (2009) procedures have been used in Vs test. Finally, correlation between rate of settlement and LPI results were determined.
Results and discussions
Outcomes of this study can be explained in below:
1. Results obtained from comparison of both methods in liquefaction potential evaluation have been showed, agreement between two methods have been happened rarely. Specially, with assuming cementation condition in soils, LPI obtained from Vs method is more than SPT. Although, different factors can be affected at uncertainties in SPT results such as type of drilling machine, energy efficiency and accuracy of test performing. Also, in shear wave velocity method, maximum velocity for occurring liquefaction in soil layers () related to fines content percentage. It is possible that boundary values in procedure not compatible with geotechnical properties in study area.
2. Evaluation of probable rate of settlements in soil layer in study area have been showed that settlement values obtained from Vs is more than SPT. This condition is compatible with LPI amounts.
Conclusions
In sum up, settlement due to liquefaction in saturate soil layers is one of the important phenomena in geotechnical earthquake engineering. Maximum rate of settlement in soil layers in study area is equal 0.45m based on SPT method and 0.9m according to Vs procedure which should be considered. Accordingly, serious damages can be inflicted to buildings, underground structures and life lines in study area.  Therefore, it is suggested in future researches with using empirical and numerical (or soft computing) methods based on field and experimental tests results a detailed assessment conducted and influence of various parameters on settlement of soil layers be determined and the items listed below should be considered:
- Cementation parameter (C) values of soils in shear wave velocity method maybe not compatible with geotechnical properties in study area. It should be evaluated exactly. 
- In this research, peak ground acceleration (PGA) value was selected based on code 2800-ver.5. As regard to Tabriz Metro Line 2 is beside to Tabriz North Fault, PGA value according to historical earthquake catalogue and seismic risk analysis should be evaluated and seismic hazard have to determine with accuracy. 

---

Introduction
The selection of the best subsurface exploration methods corresponding to geotechnical, topographical and economic circumstances of the project is one of the most effective factors in the success of a tunneling project. On the other hand, the development of a decision model and consequently choosing the most suitable alternatives are complex tasks. Therefore, prioritizing and selecting the best subsurface exploration methods, as the main aim of this study, can reduce the economic and social costs associated with the execution of a tunneling project. For this purpose, ten experts from tunneling community have been asked to weigh the chosen criteria of the problem in this research. A methodology utilized in this study is the Analytical Hierarchy Process (AHP) which proved useful in engineering decision-making problems. The other method is TOPSIS, one that has continuously been used in decision making in the recent decades. These two techniques have been combined and utilized in this work to rank the aforementioned exploration methods.
 
 
Material and methods
The study area is located about 109 km far from Shahrekord city amid the Zagros mountains. The Sabzkooh tunnel development plan has been under evaluation in the time that this research was being done. The geology of the area majorly encompasses sedimentary rocks which have been outcropped as folds, faults, and fractured and altered zones. However, the variety of the lithological units in the tunnel route is high, and units of limestone and shales also exist over the path.
In the first step, geological and topographic maps were produced for the study area, and general information from the tunnel path was collected and examined. Suitable exploration methods were evaluated, and six main methods were chosen to be considered as the alternatives of the study including (a) discontinuities study, (b) Lugeon tests, (c) RQD, (d) Geo-electric, (e) Gravimetry, and (f) Seismic methods. Moreover, nine criteria ranging from “volume of the available information” and “environmental impact” to “cost” and “accuracy” of the employed methods were taken into account. A pairwise comparison matrix was then developed, and the experts were asked to fill it out. The importance of each criterion was then simply calculated through this matrix. Alternatives pairwise comparison matrices were also filled out and, in this manner, the alternatives could be ranked using the AHP technique. Next, the TOPSIS technique was employed for the same purpose using a rather different process.
Results and discussion
Both AHP and TOPSIS techniques show very close results for ranking of the alternatives in this study. They rank the three Seismic, Geo-electric and Discontinuities studies methods as the best ones for the considered tunneling project. The only difference between these two techniques is how they determine the worst method. The AHP ranks the Lugeon as the last rank among the six methods whereas the TOPSIS determines the RQD as the least reliable method of exploration for the Sabzkooh tunnel project.
 
 
Conclusion
The prioritization and the subsequent selection of the most reliable exploration methods for an underground excavation project is a crucial task amid technical decisions. In this research, two major multi-criteria decision-making methods including AHP and TOPSIS were considered and applied for the Sabzkooh tunnel in Iran. The results indicate a high agreement between the two methodologies even though these two approaches decide differently on the least reliable methods to be applied.

---

of initial texture and plastic fine content has been investigated. In order to model the different deposition condition four different specimen preparation methods are used. The tests are conducted on sands and its combination with to 25 fine percent using static and cyclic method. Four different specimen preparation methods are used consisting of Dry Funnel Deposition (DFD), Water Sedimentation (WS), Moist Tamping (MT) and Air Pluviation (AP). Each method may be useful to model the behavior of some special deposited soils. For example, DFD is the best method to model the natural deposition of silty sands, while the WS method is suitable for simulating the sands natural deposition in the rivers. MT method is the best method to simulate the behavior of compacted embankments while SD method is suitable to prevent the segregation of particles with respect to other methods like Water Pluviation. Effect of sample preparation methods on the sands and silty sands behavior is investigated to some extent, while the clayey sands are rarely considered. The results show that water sedimentation method, Air pluviation method and Dry funnel deposition method tends to create dilative textures with continuous dilative behavior even in high clay fine content. In contrast, the wet tamping method shows the dilative behavior only for clean sands specimens and increasing fines will results in instability and complete strain softening behavior. In cyclic loading the difference between different methods are less with respect to static loading. It seems that the results of different methods tends to be identical with increasing in fine content. In general, increasing fine content up to 20 percent results in more instability. After that a threshold value can be estimated bout 2 percent, after that the instability becomes less with continuous increase in fines content. The valuable result is that the threshold value is not dependent on the specimen preparation method or initial texture of specimens. 

---

Introduction
When saturated sandy soils are subjected to seismic loadings, the pore water pressure gradually increases until liquefaction happens and settlement occurs during and after an earthquake. The mentioned problem is attributed to rearrangement of grains and redistribution of voids within the soils. Over the years many methods have been presented to increase liquefaction resistance. However, the main methods utilized in liquefaction mitigation are classified as densification, solidification, drainage and reinforcement techniques. Utilizing scrap tires in soils is a kind of soil reinforcement which has a wide range of application.
Waste material expulsion is one of the environmental problems each country faces. Accumulation of non-degradable polymeric materials in landfills has serious environmental consequences. Efforts to find new ways of soil reinforcement have drawn the attention of researchers towards the use of new recycled materials like scrap tires derivatives. Derivatives of scrap tires have different applications in civil engineering such as reinforcing soft soil, as a drainage layer in landfills and as filler materials.
Material and methods
In this paper a series of 1g shaking table tests were performed to investigate on the effect of tire powders-sand mixture in reducing liquefaction potential, settlements after earthquake and pore water generation. Shaking table is made of Plexiglas with inner dimensions of 200×50×70 cm. At bottom of the container a void chamber is made by using a number 200 sieve so that the saturation process could be done gradually and uniformly. A plastic plate was rigidly fixed at the center of container to separate reinforced and unreinforced samples from each other and waterproofing carefully. Therefore two models (reinforced and unreinforced) can be tested at once with the same input acceleration. An absorbing layer of foam with 2 cm thickness was employed to decrease the effect of boundary conditions in order to avoid a direct confrontation model with a rigid container. Firoozkuh No. 161 sand and tire powders were used for the mixture in reinforced side, and pure sand in unreinforced side. In this study 4 mixture ratio (TC=5%, 10%, 15% and 20%) were done. Both of unreinforced (pure sand) and reinforced (tire powders-sand mixture) models were prepared by wet tamping method, in which soil is mixed with 5% water. Each model was prepared in six layers. The required weight for each layer was considered based on the desired density (relative density is zero) and exact volume of the layer. Each portion was placed into the model container and then tamped to reach desired level. Carbon dioxide (CO₂) was allowed to pass through the specimen at a low pressure in order to replace the air that trapped in the pores of the specimen. Then water was allowed to flow upward through the bottom of the container at low pressures in order to flush out the CO₂ that cause increasing the final degree of saturation. Vibration with approximate uniform amplitude and 2 Hz frequency was applied to the container.
Results and discussion
Results indicate that acceleration within the soil tends to be increased towards the soil surface. On the other hand, after initial liquefaction (that occurred at un-reinforced models), acceleration is decreased due to the increase in excess pore water pressure. Also, it can be seen that the increase in tire powders ratio remarkably reduces the maximum excess pore-water pressure ratio. The settlement of the tire powders-reinforced models is significantly less than the unreinforced models, and with the increase of the tire powder percentage shows a very small increase of volume. The outcomes show that the value of the mean damping ratio versus the shear strain range of 0.01 is increased with the increase in tire powder content. Shear modulus is obtained from the ratio of the difference in maximum and minimum stress and strain developed in desired loop. The maximum of the shear modulus in reinforced models is more than the unreinforced models.
Conclusion
The main aim of the present paper was to investigate the influence of reinforcing a saturated sandy soil with tire powders on the soil dynamic properties and the mitigation of liquefaction potential. The following conclusions were drawn from this research.
- The increase of pore-water pressure leads to a reduction in soil shear stiffness and acceleration amplitude.
- Reinforcing sand with tire powders reduces the excess pore-water pressure ratio because of liquefaction and increases liquefaction resistance. 
- Reinforcing sand with tire powders decreases settlement caused by liquefaction.
- The damping ratio decreases at large shear strain as liquefaction occurs.
- Maximum shear modulus and mean damping ratio of reinforced soil has been increased with increasing tire powder content in the mixture../files/site1/files/123/3BahadoriFarzali.pdf

---

Introduction
Many studies have shown that the lime stabilization method can increase the strength and hardness of cohesive soils. Increasing these parameters is dependent on several factors such as curing time, lime content, clay minerals, soil particle size and moisture content.
When lime is added to moisture clay soils, a number of reactions occur to improve soil properties: 1- short-term and 2- long-term reactions. The short-term reactions include cation exchange, flocculate and carbonation; whereas, the long-term reactions include pozzolanic reactions. Since adding lime changes clay particles structure, it can change shear strength parameters.
Using geogrids as reinforcement in soil mass creates a composite system in which the soil tolerates compressive stresses. The elements of the reinforcement are also responsible for tensile stresses and interaction the reinforcement elements and soil increases the strength and ductility. The mechanism of stress transfer is based on interaction between soil and reinforcement. Accordingly, one of the most important issues in the analysis and design of reinforced soil structures is determination of frictional resistance parameters in soil-geogrid interface (adhesion and friction angle) which is discussed in this paper.
Stability and performances of reinforced earth structures significantly depend on the shear behavior of interface soil-geogrid in different weather conditions. Factors such as rainfall, seepage of groundwater and seasonal changes influence on soil moisture content. Changes in moisture content or soil dry density change interface soil-geogrid resistance. Increasing moisture content reduces the shear strength of reinforced soil and sometimes leads to large deformation or failure of system.
In this study, clayey soil with low plasticity (CL), hydrated lime for soil stabilization and two types of geogrid with different aperture size for reinforcing were used. In order to improve the brittle behavior of lime stabilized soils and to increase ductility of the samples, in the present study, lime stabilization and geogrid reinforcement was investigated, simultaneously. The interface shear strength parameters of treated soil with different lime content-geogrid and reinforcement coefficient were determined by direct shear tests. In addition, to study the effect of moisture content on interface shear strength soil-geogrid, all samples were subjected to shear in optimum and higher moisture content because the long-term performance of reinforced cohesive soils exposed to seasonal variations is evaluated.
Material and methods
The selected soil for the study is clayey soil from south region of Tehran, Iran. According to Unified Soil Classification System (USCS), the soil was classified as CL (clay of low plasticity).
In this study, three series of specimens were prepared and tested as follows:

• Stabilized samples with 0, 2, 4 and 6% lime for 7 days curing time
• Reinforced samples by geogrid (with and without transverse ribs of geogrid)
• Reinforced stabilized samples with different lime contents (0, 2, 4, 6 and 8%) by geogrid (with and without transverse ribs of geogrid) for 7 days curing times

To investigate the effects of bearing resistance provided by the transverse members of the geogrid and their contribution to the overall strength for reinforced soil sample, numerous tests were conducted with the geogrid without transverse members (all the samples had the same number of longitudinal members of the geogrid).
Direct shear tests were carried out on specimens based on ASTM D5321 at constant horizontal displacement rate of 1 mm/min.
Results and discussion
The results reveal that the shear strength of the stabilized soil increased and there are maximum values in an optimum lime content which is about 4%. Increasing lime content to an optimum lime content of clay caused the maximum changes in clay minerals because of cementitious and pozzolanic reactions and increases the strength of the clayey soil. Reduction of strength by adding lime to the soil more than the optimum content may be caused by the following reasons:
1. Stopping pozzolanic reactions because of finishing reactance during reaction
2. Making difficult the release of limewater (Ca OH ₂) in the cementitious context of soil.
Until SiO₂ and AL₂O₃ are not finished, pozzolanic reactions continue and produce cementitious product, thus the shear strength increases and improves the long-term performance of the stabilized soils.
Reinforced soil samples have higher shear strength relative to samples without reinforcement subjected to the same normal stress. This increase in shear strength is mainly attributed to the interlocking of soil particles that penetrate through geogrid apertures. In addition, geogrids restrain particles´ movement and thus increase the mobilized frictional resistance at particle contact points.
Increasing in lime content to 4% (optimum lime content in this study) has significant effect on the development of adhesion and then decreases gradually with increasing of lime content from 4 to 6%, while friction angles remain constant approximately.
Adhesion and friction angles decrease with increasing moisture content.
The results show that the reinforced stabilized specimen with 4% lime has the maximum value of reinforcement efficiency. The increase in moisture content can significantly reduce the reinforcement efficiency.
It is clearly observed that the reinforcement coefficient of reinforced stabilized sample by geogrid that has smaller aperture opening size (4Í4 mm) is higher than reinforced stabilized sample by another geogrid (10Í10 mm) in optimum and higher than optimum moisture content.
Conclusion
One hundred and twenty samples in 3 specimen categories including lime treated, reinforced and reinforced treated samples were prepared for the current study for 7 days curing time in optimum content and higher than optimum content. The main results can be concluded as:
The test results indicate that the shear strength of stabilized clayey samples increases after 7 days curing time due to pozzolanic reactions.
The results show that reinforced samples have higher shear strength relative to unreinforced samples.
Adhesion and friction angles and reinforcement efficiency decrease with increasing moisture content.
The reinforcement coefficient of reinforced stabilized sample by geogrid 1 that has smaller aperture opening size is higher than by geogrid 2. In general, interaction between particles and geogrid with smaller mesh size is stronger because of matching the size of soil particles and meshes../files/site1/files/123/8Extended_Abstract.pdf
 

---

Introduction
Hazardous waste (solid, liquid or contained gases) is a waste with properties that make it potentially dangerous or harmful to human health or the environment. Site selection and suitable conditions for hazardous wastes landfill is considered as the final stage of waste management that they have high sensivity. The purpose of this study is to identify prone areas to hazardous waste landfill for Chaharmahal and Bakhtiari province using geographic information systems (GIS) as an important tool for the analysis of potential sites and the Analytical Hierarchy Process (AHP) and to provide solutions to optimize the positioning is executed. Firstly, criteria and limitations of environmental, economic, social and physical were determined, then layers of the criteria in GIS were prepared. In this study, the inappropriate areas were first removed from the model, and the suitability of remaining regions as a categorize criterion considered. Categorize criteria for paired comparison using AHP as an efficient tool for determining the relative weight parameters are used to measure and rank the expert choice application imposed. Then the implement paired comparison of the relative weights of the criteria and sub-criteria and criteria for each category were determined. After calculating the net weight and normal weight, normal weight based on standard maps in the GIS environment has been classifieds. Finally, by combining maps and applying criteria FA map, the final map was extracted.
Material and methods
The purpose of this research is to identify and prioritize appropriate areas of special waste disposal using multi-criteria decision-making methods. In order to locate using the GIS, first, identifying, evaluating and selecting criteria and constraints for the construction of landfill, in order to reduce the economic, environmental, and health costs. In the multi-criteria evaluation method, criteria are the basis of decision making, so that a set of criteria is combined and combined to achieve a single combination. In this paper, a two-stage process was used to locate the landfill site. In the first stage, which is recognized as the identification stage of prohibited areas according to different criteria, the study area is divided into two appropriate and inappropriate classes that will be eliminated as prohibited areas for the construction of landfills. In the second stage, the various factors are ranked and weighted according to the relative importance and, finally, places that receive the appropriate points are introduced as areas susceptible to the dumping of special wastes. In order to obtain the digital data of the criteria in the GIS environment from the digital elevation map (DEM), the specifications of the piezometric wells information are available from the regional water organization of the province. The available data such as geological map of the province at a scale of 1: 250,000, satellite images of Landsat and map of land suitability of the province, rainfall data of the synoptic stations of the province and the data of the Environmental Protection Agency were used. In general, the following steps have been taken in the process of locating:
- Identification of effective locating factors (limitations and factors)
- Digitizing and providing the required layers of information using the GIS package
- Identify and eliminate prohibited and inappropriate areas for landfill construction
- Classification and weighting of the factors and layers of information sought
- Integration of layers and the provision of a mapped rate and talent to determine the appropriate areas.
Results and discussion
1. Set limits
In this study, in order to select suitable sites for landfill particular, the criteria and limitations were determined. The information layers for each of the criteria were provided in the GIS environment.
2. Classification and weighting criteria
In the second stage, which is the stage of weighting and rating, of 14 effective criteria were used in site selection. AHP is one of the most efficient techniques, multi-criteria decision. This method is based on comparing factors and to study various scenarios to give managers and decision makers. This technique is one of the most comprehensive system designed for decision-making with multiple criteria.
3. Editor hierarchy to locate
Hierarchical structure is a graphical representation of a real complex problem, which mainly target the problem and at the next general criteria, sub-criteria and options are the way in AHP is used to calculate points based on comparison test.
4. Shipping Weight Matrix Binary comparison and decision-making
After compiling a hierarchical structure, the next step is to evaluate the elements by comparing the test. In general, if the number of options and criteria respectively m and n are then paired comparison matrix of options for comparison matrix m × m and n × n matrix will be a couple of criteria.
5. After weighing and preparing the normal weight of the options, the normalized weights in the GIS environment were added to the criteria map and the Raster and Weighted layers of each criterion were prepared. Due to the wide area of the studied area, the size of each pixel was 50 * 50 m. Then, using the Raster Module, the Criterion Map was combined and a zoning map was prepared for the special waste disposal site.
Conclusions
In the present study, according to various criteria influencing the Hierarchical Analysis Process for prioritizing the criteria in decision making, based on the results, the talent map of the area was prepared for special waste dumping, in which according to the final score of the layers, the area was classified into four appropriate, relatively suitable, relatively inappropriate and inappropriate classes. Suitable areas were 12.64%, relatively fairly 32.47%, relatively inappropriate 30.43%, and inappropriate zones 9.58% of the area of the talent map were included.
./files/site1/files/123/5Extended_Abstract.pdf

---

Introduction
The design engineers usually follow a specific decision-making process for optimal selection of the type of required foundation and its design. In this state, in case the surface foundation is not appropriate for the project conditions, before making any decisions about the use of deep foundations, the proper methods for optimization of the liquefied soil should be evaluated in order to compare the advantages and disadvantages of each of them with those of deep foundation, in terms of efficiency, implementation problems, costs, and finally to select the best choice. One of the best methods of soil improvement is the use of stone columns. The rationale behind the use of stone columns is the high shear strength of materials and the provision of lateral grip by surrounding soil. Therefore, the stone column can receive the load from the structure, and transfer it to the resistant layers. In the soils with low shear resistance, the lateral constraint crated by the surrounding soils is not enough for preventing the sideway buckling of the column under which is subjected to the loads. Thus, special measures should be considered for the use of stone columns in these soils. One of these methods is the use of reinforcement shelves such as geogrid and geotextile. Investigating the previous studies, the lack of evaluation of the design parameters such as the settlement ratio of the soil improved by the reinforced stone column to the geogrid, and provision of design graphs in this regard, has been revealed. Therefore, by extension of the studies conducted by Chub Basti et al. in 2011, the design graphs were provided in this regard.
Material and methods
The PLAXIS V8 Software was used for modelling the soft soil improved by the stone column. For increasing the precision of the results, the 15-knot element was used in the current study. The fine mesh was used in the models made for the analysis of the problem. For simulation of the improved soft soil with the stone column in a single cell, the modelling was implemented in a two-dimensional environment in axial symmetry conditions. In the current study, it was assumed the rigid foundation is on the improved bed. Thus, for analysis of the simulated model, a vertical strain up to 2% of the soft soil height has been applied on the ground. Also, for simulation of the soil behavior, an appropriate model of soil and parameters proportional to the materials should be allocated to the construct geometry. The non-linear stress-strain of the soil in different levels of the problem can be simulated. The number of model parameters increases with the level of problem rupture. For precise simulation, we need the proper parameters of the materials. For modeling of soft soils and stone columns, elastic-plastic model with Mohr-Coulomb rupture criterion was used. In the current study, it was assumed the soft bed is located on a very hard layer of soil. Therefore, the vertical deformation was prevented on this horizontal boundary. Also, the horizontal deformation in two vertical edges was prevented and only deformation in vertical direction was allowed. The soft bed close to saturation was considered without the determined free water level. For models with stone columns, the element of interface between the stone column and soft soil, has been used. The reason behind using this element is that the stone column rupture is of shear form and due to this, a significant shear stress is created on the common surface between the stone column and soft soil. The percentage of the replacement area is defined as the ratio of the total area of the stone columns to the total area of the non-improved area. In the current study, the percentage of the replacement area is utilized between 10 to 30%, which is used in implementation. Also, the diameter of the stone columns is from 0.6 to 1.2, in the analyses.
Results and discussion
The results of the numerical study were compared with the existing theoretical relationships provided by Poorooshasb and Meyerhof (1997), and Pulko et al. (2011). Figure 1 shows the comparison of the replacement percentage (RP) and settlement ratio (SR) in the non-reinforced state in the current study as well as theoretical relationships proposed by the previous researchers. Based on this figure, there is a difference between the results of the current study and those of Poorooshasb and Meyerhof (1997), and Pulko et al (2011). Poorooshasb and Meyerhof (1997) calculated the settlement ratio in their proposed material with the assumption of linear elasticity of the materials without consideration for plastic settlement. Therefore, the settlement of the improved soft soil with stone column, calculated by Poorooshasb and Meyerhof, would not show the real amount. However, Pulko et al. (2011), with consideration for the elastoplastic behavior of the materials, the lateral expansion of the stone column, and the primary stress of the soil around the column, provided more realistic results that correspond closely with the present study. Also, for designing the stone column, the results of its reinforcement have been also provided in the graph presented in Figure 2. Thus, by the use of these graphs, the ratio of settlement reduction can be obtained for each distance between columns and with different percentages of alternatives../files/site1/files/124/2jalili%DA%86%DA%A9%DB%8C%D8%AF%D9%87.pdf

---

Introduction
     Ground-penetrating radar (GPR) method is a pretty new, non-destructive and high-resolution geophysical method that is widely used to identify the thickness of snow and ice layers and glaciers bed, because snow and ice are transparent for electromagnetic (EM) waves. Therefore, this method has been used to determine the thickness and basement topography of Alam-kooh glacier. In this research, only the GPR acquired data using unshielded antenna with central frequency of 25 MHz along one line in Alam-kuh glacier, Kelardasht- Mazandaran, have been processed and interpreted. The GPR data acquisition has been done by using common offset method, and transmitter-receiver separation of 6 meters. The final real radargram related to one of the surveyed GPR profiles in this region has been prepared after applying various processing operations containing signal saturation correction, amplitude gain, f-k migration filtering and static (topography) correction on the raw data. After applying processing sequences on the acquired data, the EM waves reflection off the interfaces of different layers including the reflections from the glacier basement have been detected, and by assigning a suitable EM wave velocity in the ice (0.16 m/ns), the thickness of 50 m for the ice layer laid under the survey line has been estimated. Also, in present research, forward and inverse modeling of GPR data have been performed to employ for snow, ice and glaciological investigations in the AlamKooh region of Mazandaran. To achieve this goal, GPR response of synthetic model corresponding to the real radargram was simulated first, by 2-D finite-difference time domain (FDTD) method. Afterward the inversion method by solving an optimization problem was employed to validate the interpretation of real GPR data.
Methodology and Approaches
     Based on the nature, physical and geometric properties of the subsurface target in the field data, their synthetic model have been built and their two-dimensional GPR responses forward modeling using ReflexW software and finite difference algorithms improved in the frequency domain, have been obtained. Also, it has been used an effective algorithm, coded in GUI environment of MATLAB programming software and as a result, a reliable and accurate inverse modeling has been carried out. In the present study, to simulate the behavior of the propagation of EM waves in GPR method, two-dimensional finite difference method has been used. The main advantage of this method is its comparative simplicity of the concept, high accuracy and simple implementation for complex and arbitrary models as well as easily adjusting the antenna when applied. In this study, acquisition of GPR field data and synthetic data modeling have been carried out in TM mode. The radargrams of the GPR data have been demonstrated using ReflexW software after performing necessary processing sequences.
Results and Conclusions
     The obtained results reveal that moraine materials covering the surface of the area are mainly fine-grained granite. The bed-rock or basement in the area is also granite. The polarity representing ice-bed rock is clearly seen on the GPR profiles. The topography of the glacier basement has successfully been detectable using just by GPR method. The electrical resistive nature of the glacier has caused large penetration depth of GPR pulses in this research. Furthermore, the results of the research for presented profiles on the basis of forward and inverse modeling output of GPR data in comparison with real GPR radargrams in the region validated the accuracy of GPR investigations in the area. Although with a quick glance, the error obtained by the inverse modeling for real GPR data seems unexpected and unacceptable, absolutely the high rate of error depends on many factors influencing on the real earth models containing various limitations existing in all forward modeling algorithms and software packages, impossibility of making forward modeling exactly according to the real models (due to the complex nature of the ground), taking into account the homogeneity and uniform host environment and targets in the modeling process unlike the diversity, the presence of different types of noises and so on. Therefore, making a controlled geophysical test site and trying performance of inverse modeling algorithm for field GPR data in this site, as well as determining the important physical parameters such as dielectric permittivity and electrical conductivity by experimental method through sampling from different depths for complex geological environments are suggested../files/site1/files/124/1ahmadpur%DA%86%DA%A9%DB%8C%D8%AF%D9%87.pdf

---

Introduction
Every year, numerous casualties and a large deal of financial losses are experienced due to earthquake events. The losses incurred by an earthquake vary depending on local site effect. Some well-known examples include the earthquake in Caracas 1967, Mexico city 1985, Kalamata 1986, Loma Prieta 1989, Roodbar 1990, Bam 2003, Jammu and Kashmir 2005, Sichuan 2008 and Haiti 2010. Therefore, in order to conquer drastic effects of an earthquake, one should evaluate urban districts in terms of the local site effect. Various methods are available for the evaluation of site effect. One of the most common methods includes ambient noise survey. Today, this approach is being used as fast, applicable, cost-effective method. Ambient seismic noise are feeble ground motions with displacement amplitudes of about 0.1–1 μm and that can be detected by seismograph with high magnification. Many investigations have been conducted to determine the nature of ambient noise. One of the possible sources of ambient noise can be human activity, such as traffic, industrial noises and nature activity, such as wind, ocean waves. The Babol city is one of the largest cities in the north of Iran (Mazandaran province). It lies on alluvium beds in the region presenting a high seismic potential. Therefore, comprehensive studies are necessary to introduce suitable solutions for minimizing earthquake damage and loss of life. For this reason, in Babol city, ambient noise survey has been performed at 60 stations and the obtained data were analyzed with Nakamura or H/V method (1989). The results were compared with local geological, geotechnical and seismic data to confirm their reliability for a seismo-stratigraphic.
Methodology and data collection
The analysis of ambient noise was initially proposed by Kanai and Takana (1961). Since then, many researchers have used ambient noise for site effect evaluation. As it is said before, one of the most popular techniques for estimation of site effects in the regions with low seismicity is ambient noise survey by Nakamura or H/V method (1989). Based on the literature review, the Nakamura method (1989) has been used in many places. Many theoretical and experimental studies show that, this method has the capability of estimation of fundamental frequency. Ambient noise survey was carried out at 60 stations in Babol city. Ambient noise was recorded using a velocity meter SARA. Two horizontal and a vertical components of ambient noise at each location are recorded for duration of 15 min with 100 samples per second. Because the environmental noise has an effect on ambient noise they are recoded between 10 p.m. to 6 a.m. The locations were determined by using GPS at the sites. The ambient noise survey in this study was made in compliance with the guidelines of SESAME (2004).
Results and discussion
The maximum and minimum values of fundamental frequency of the present research are 11.4 and 0.65 Hz, respectively. Also, the maximum and minimum values of amplitude of H/V peaks have been calculated as 3.71 and 1.19, respectively. The most significant point is that the fundamental frequency of the major part of Babol city is smaller than 1 Hz in agreement with the previous knowledge of the city geological setting. Another relevant point is the presence of some stations with very high (> 5 Hz) fundamental frequencies. In these cases, ambient noise recording and data analysis were repeated, but similar results were obtained. Considering the lack of sufficient geotechnical data (in some stations), the above phenomena need to be integrated with other methods. For this purpose, the data of electrical resistivity tomography (ERT) were used. The ERT showed that small regions of the north-west, west, and south areas have high resistance values probably related with the presence of hard deposits in the shallow subsoil. Babolrood river diversion in the west part and its return to the previous direction in the northern section is possibly due to the existence of these relatively hard deposits. By comparing these two tests, we observed that the ERT results correlate with the ambient data analysis. Therefore, we can conclude that the high-frequency peaks measured are reliable, but we need direct investigation to associate them to a specific shallow geological layer. To validate the results, fundamental frequencies obtained from ambient noise survey were compared with geotechnical data, numerical analysis and seismic data in the study area. A general review shows that the geotechnical data, equivalent-linear analysis results and seismic data have an acceptable conformity with the results of ambient noise survey.
Conclusion
The results show minimum and maximum fundamental frequencies 0.65 and 11.4 Hz, respectively. Assessment also reveals that the major parts of Babol city have the fundamental frequencies less than 1 Hz, which are in conformity with that of previous research. According to the results of seven cross sections, it can be concluded that fundamental frequency variations are in line with the geotechnical and geological data in the study area. It means that this method is the appropriate way to assess the local site effect in the Babol city. It is also observed that besides the soil layers, the soil stiffness and its shear wave velocity are effective factors in changing the fundamental frequency. Site frequencies were also estimated by preliminary 1-D site modeling using the equivalent-linear method. In general, a reasonable correspondence between the methods was obtained. Using seismic data, the HVSR of two strong ground motions have been calculated and the results have been compared with the nearest ambient noise recording station. Analyzing the spectral ratios demonstrates that the value of the fundamental frequency obtained by the H/V method (1.06 Hz) is very close to that of frequencies obtained by the seismic data (0.95 and 0.90 Hz)../files/site1/files/124/6rezaee%DA%86%DA%A9%DB%8C%D8%AF%D9%87.pdf

---

Introduction
Study on the main characteristics of strong ground motions, has relatively long history. The observations and investigations on the structural damages after strong earthquakes such as Northridge 1994 in California, Kobe 1995 in Japan, Tabas 1978 and Bam 2003 in Iran, are representatives of the destructive effects of strong near-field records. The most important specification of the near-field records which distinguish them from far-field records, is their ability to generate energized and relatively short-duration acceleration spikes as well as high amplitude and long-domain velocity pulses. Moreover, according to the lack of accurate statistical profiles as well as many deficiencies, processing the spectral existent data is not able enough to fully explain the seismic tremors. Based on the fact that the great earthquakes have long recurrence interval and also many high seismic zones of Iran do not possess strong tremors, hence generating and simulating feasible great events is required by applying closed form models and analysis of available data. In this study, in order to simulate the existent pulses in the time history of near-field records, the developed mathematical configuration is presented by analytical comprehensive attitude on the closed form model by Mavroeidis and Papageorgiou (2003).
Material and methods
Simulation of strong ground shakings, especially in areas where there is limited recorded data, plays a key role in assessing dynamic behavior of structures. Owing to unique characteristics of strong near-field ground motions, it is not possible to determine exact effects of these strong records on structures using simplified mathematical models. It is feasible to develop more complicated models which represent much more characteristics of near-field ground motions. Mavroeidis and Papageorgiou (2003) studied the parameters affecting near-fault ground motions. Their studies resulted in introducing a mathematical model capable of interpolating velocity pulses of near-field earthquake records (MP model). This closed-form MP model interpolates long duration pulses using a set of input spectral parameters.
 The pulse period, the pulse amplitude, the number and phase of half cycles are the key parameters that define the shape of velocity pulse. Thus, a four-parameter model has been developed to describe velocity pulses which contain forward directivity effects. In this research, it was observed that by using a combination of cubic and exponential terms, an enhanced model for interpolating the pulses presented in near-field earthquake records could be achieved (EMP model). Figure 1 shows the analytical interpolation of acceleration and velocity time histories using MP and EMP models.
 
Figure 1. Fitting of acceleration time histories with MP and EMP models
 
Results and discussion
Based on the obtained results, it is observed that there is a striking similarity between analytical characteristics obtained by actual earthquake records and mathematical pulses. Moreover, using the enhanced closed-form model (EMP model) reduces discrepancy between the results obtained under actual and the synthetic earthquake records.
Conclusion
Findings of this research reveal that equivalent pulses could be a good representative of actual earthquake records analytically, in order to assess the seismological characteristics of these tremors. It is worth mentioning that modelization of forward directivity pulses displayed in time history of strong ground shakings, is an efficient measure in evaluating seismic response of structures. In addition, due to stochastic nature of earthquakes, computational uncertainties and descriptive limitations of analytical parameters, using closed-form models require a high level of accuracy../files/site1/files/124/8nazari%DA%86%DA%A9%DB%8C%D8%AF%D9%87.pdf
 

---

Introduction
Pressure tunnels in hydroelectric plants are used to convey water to powerhouses. These tunnels are the sources of seepage flow to the rock formation, thus, during the water filling, they will have a low resistance to seepage and, by increasing the internal water pressure of the tunnel, the inflow force will be transferred to the rock mass. In these conditions, the cracks, pores and all other elements of the rock mass are affected by the seepage forces in all directions. This hydro-mechanical interaction affects changing the stresses and displacements of the rock mass around the tunnel and causes modifications in the permeability of rock elements during the water filling. Therefore, changes in stress distribution lead to alterations in the permeability coefficient and redistribution of the seepage field. In these conditions, since the analytical solution of the problem is not possible, the numerical analysis based on the finite element method has been used in this study.
Material and methods
In this approach, the rock mass is considered as an equivalent continuum in which the effects of discontinuities are taken into account in its material behavior. High-pressure tunnels under internal water pressure requires reinforced concrete lining to prevent hydro-fracturing. The ABAQUS software is capable of analyzing such as seepage from the tunnel, modeling of the steel bars in concrete, and taking into account hydro-mechanical interaction. Thus, the software is used for numerical analysis.
The pressure tunnel of the Gotvand dam and hydroelectric power plant (HPP) scheme is taken as a case study for the numerical simulation. Pressure tunnel of the Gotvand dam located in the southwest of Iran is taken as a case study for the numerical simulation. Among behavioral models in the software, Mohr-Coulomb failure criterion is considered to describe the rock mass, but the principle of effective stress determines the rock mass behavior. Since the concrete lining of the pressure tunnel will undergo two mechanisms of the cracking due to tension and the crushing due to compression, concrete damaged plasticity model is used to predict the response of the concrete elements. The evolution of the yield surface of the concrete lining is also controlled with tensile and compressive equivalent plastic strains, correspondingly.
In this study, the hydro-mechanical interaction is implemented based on the analysis of the pore fluid/deformation analysis, and the direct-coupled method is used to solve the governing equations of the problem. To verify the proposed model, the elastic behavior of the media is simulated to compare the numerical and the analytical solutions and good agreement is obtained. The numerical analyses are carried out the hydro-mechanical interaction with constant permeability coefficient. When cracks develop in the concrete lining at high water pressure, the properties of the concrete lining change and as a result, the stress dependent permeability of the lining and surrounding rock mass in pressure tunnels should be considered. The coefficient of permeability controls the rate of seepage flow in porous and fractured media. Although permeability represents an original property of the porous media, it can be modified when subjected to the stress variations. Instead of changing aperture, the change in the void space or volume is the typical consequence that results to change the permeability coefficient. In order to bring the model closer to the real conditions and in the validation of the new model, the influence of the permeability coefficient variations of the concrete and rock mass on the deformations and stresses of the model has been added to nonlinear analysis by USDFLD code. Increasing the water head in the tunnel during water filling is also considered with the combination of DLOAD and DISP codes in the model.
Results and discussion
Since the lining and rock mass have nonlinear properties and complex behavior, for verification of the model in ABAQUS software, the model is simulated with homogeneous, isotropic and elastic behavior. The results of seepage flow on the interface of the concrete lining and rock mass obtained by analytical and numerical solutions indicate that there is a ±5 % difference between them. Then, the results of the elastic behavior of the model show a good agreement with the results of analytical solutions. Therefore, this numerical model has been employed for the nonlinear analyses.
Finally, the optimal thickness of the concrete lining with the appropriate arrangement of the reinforcement in the reinforced concrete linings is utilized to minimize water losses from the tunnel based on the new model. Thus, the results of the analysis with the aim of reducing the water losses from the tunnel indicate that the suitable arrangement of the steel bars in the concrete lining leads to the distribution of micro cracks in the lining, and the reinforcement stress stays at a lower value with high internal water pressure. Based on the new numerical model, it is suggested that the lining should be designed with the thickness of 40 cm and the reinforcement with the diameter of 16 mm and the spacing of 20 cm.
 Conclusion
The results of the numerical model indicate that to control the seepage outflow from concrete-lined pressure tunnels, the thickness of the lining and the suitable arrangement of the steel bars in the concrete lining play a significant role in preventing excessive seepage from the tunnel./files/site1/files/124/3dadashi%DA%86%DA%A9%DB%8C%D8%AF%D9%87.pdf

---

---

Steel sheet pile walls are being widely used as earth retaining systems. Sometimes loose or soft soil layers are located in various depths in an excavation. This issue causes different effects on ground surface displacements, forces and moments acting on sheet piles and struts during excavation procedure, compared with a status that soil is totally uniform. These differences are not exactly considered in conventional design methods of sheet pile walls. In this paper, a deep excavation using finite element method is analyzed. Excavation’s depth is divided into three different layers. One of three layers is a loose soil layer and its position is modeled in three different situations, top, middle and bottom of the model. Obtained results are compared with results of excavation without the loose layer. The pseudo-static analysis is performed by applying 0.3g horizontal acceleration. The results indicate that when a loose layer is located beneath stiffer layers, bending moments acting on sheet pile wall and shear forces increase about (50~100)% and (15~50)%, respectively. Also, the middle loose layer changes the location of maximum lateral deformation of steel sheet pile wall.
 

---

Materials such as waste tire chips were widely used to improve the strength of soil. The objective of this study is to discuss the residual strength or steady-state behavior of sand-waste tire chip mixtures. A series of undrained monotonic triaxial compression tests were conducted on reconstituted saturated specimens of sand and sand-tire chip mixtures with variation in the tire-chip contents from 0 to 4 percentages by dry weight of soil. The specimens are prepared using dry deposition method of preparation. The influence on residual resistance of varying confining pressure (100, 200, and 300 kPa) and sand mixture relative density (40, 65, and 80%) were evaluated. Tests results showed that by increasing the tire chip contents, the residual strength increased and steady-state lines move to the right of log S_us-e diagram. Also, the residual resistance improvement induced by tire chip inclusions was found to be sensitive to the relative density of samples and applied confining pressure.

---

In the mechanized boring method, the factors affecting ground surface settlement can be mainly divided into five categories: geometric, geomechanic, boring machines working, operating and management parameters. In urban tunnels bored mainly in shallow soil bed, face pressure can be one of the factors preventing ground settlement. The Line A tunnel in Qom metro project is bored with an EPB (Earth Balance Pressure) mechanized boring machine. The effect of face pressure on ground surface settlement was analyzed in the present study according to five sections of the tunnel. These five sections were selected in different kilometers of the tunnel where settlement gauges were installed and the results could be validated. To investigate the effect of face pressure on maximum ground surface settlement, four pressure levels of 100 kPa, 150 kPa, 200 kPa, and 400 kPa were taken into consideration. These were 1, 1.5, 2, and 4 times of the initial face pressure level, respectively. The ground surface settlement was assessed at four pressure levels using the finite element software, PLAXIS 3D TUNNEL. The results were validated using ground-level instrumentation (settlement gauges) on all sections. The validation showed that the modeling results are in good agreement with the results obtained from settlement gauges.  Comparison of the results indicated that a 4-fold increase in the face pressure led to a maximum decrease of 4.45 mm in the maximum settlement. Therefore, an increase in the face pressure can reduce settlement, although quite minimally. It was also found that an over-increased face pressure (face pressure over 200kPa) not only did not reduce the maximum ground surface settlement but also may lead to passive failure or uplift of ground surface ahead of the shield. 
 

---

In urban areas, it is essential to protect the existing adjacent structures and underground facilities from the damage due to tunneling. In order to minimize the risk, a tunnel engineer needs to be able to make reliable prediction of ground deformations induced by tunneling. Numerous investigations have been conducted in recent years to predict the settlement associated with tunneling; the selection of appropriate method depends on the complexity of the problems. This research intends to develop a method based on Artificial Neural Network (ANN) for the prediction of tunnelling-induced surface settlement. Surface settlements above a tunnel due to tunnel construction are predicted with the help of input variables that have direct physical significance. The data used in running the network models have been obtained from line 2 of Mashhad subway tunnel project. In order to predict the tunnelling-induced surface settlement, a Multi-Layer Perceptron (MLP) analysis is used. A three-layer, feed-forward, back-propagation neural network, with a topology of 7-24-1 was found to be optimum. For optimum ANN architecture, the correlation factor and the minimum of Mean Squared Error are 0.963 and 2.41E-04, respectively. The results showed that an appropriately trained neural network could reliably predict tunnelling-induced surface settlement.

---

This study numerically investigates the bearing capacity of drilled shafts (bored piles) in clay using FLAC^2D. The results obtained in this study are compared with centrifuge test results. The results of the empirical relationships available in the literature are compared with the results of the present numerical study. A series of analyses is also conducted to assess the effects of various soil and pile parameters on the magnitude of tip and side resistance of bored piles embedded in clay. These parameters include the soil elastic modulus, pile length and diameter, undrained shear strength, unit weight, and Poisson’s ratio of soil. Furthermore, the coupling effect of soil undrained shear strength and elastic modulus of soil on tip resistance are investigated. The results show that the lower value of soil elastic modulus results to lower effect of soil undrained shear strength. The effect of soil undrained shear strength on tip resistance is approximately constant (about 83% for a change of soil undrained shear strength between 25 to 200 kPa) for the range of elastic modulus between 20 and 180 MPa. Also, a new equation is proposed to estimate the bearing capacity factor of N^*_c.