Journal of Engineering Geology

---

The Dynamic Probing is an important test in site investigation and geotechnical studies. It is used for determination of situation and specification of soil layers when the depth of exploration is moderate. In our country, in the last few years test has been widely used to replace the standard penetration. In this paper, the experiences to use different various types of dynamic probing rigs and the errors that can occur in these experiments are discussed. Then the accuracy of this test is studied to evaluate the specification of fine-grained soils. The most important innovation of this research is the proposition of the new correlations between cone dynamic resistance and undrained shear strength and so compaction percentage in fine-grained soils. The paper encourages the wider application and further development of this test for site investigation in fine-grained soils.

---

./files/site1/files/2.pdfExtended Abstract
(Paper pages175-200)
Introduction
In weak soils with low bearing capacity, the load transfer is done using piles. Therefore, by creating an interposed layer separating the pile from the raft, reactions between raft and pile head will be reduced and the load-bearing role of shallow soil will be more than contact pile situation. Normally, the pile head and shallow soil have a settlement equal to the raft. Thus, the relative settlement of pile and soil in pile head is equal to zero and at the bottom is high and the body friction mobilizes upward. In addition, a portion of load is tolerated by shallow soil and the other portion is tolerated by the pile head, which would be transferred to deeper soil layers. In noncontact state, with the formation of a hard soil layer on which the raft is located, soil mechanical parameters will be improved; while in contact state, the settlement will be decreased by reducing the amount of transferred load to the shallow soil. The transferred load to the shallow soil increases vertical and horizontal stress around piles, so bearing capacity of piles is increased.
Methodology
In this study, a parametric study has been performed concerning contact and noncontact piles using finite element software namely, ABAQUS/CAE software version 6.13.1 and the obtained results were compared (with what? The sentence is incomplete). Thus, simulations are is done for states of 0, 1, 4 and 9 piles for each of the contact and noncontact piles (total of 8 simulations). In the present research two models were taken to investigate the optimum mesh sizes, 12 models for parametric studies on parameters of piles’ length, piles’ diameter, thickness of the raft and interposed layer and one model for verification study. Models in both contact and noncontact have been considered with a one meter interposed layer. Raft width and thickness were selected 7.5 and 1.6 m, respectively. Width and depth of the soil mass used in the model were 32 and 26 m, respectively, and the distance between the bottom of the pile and the soil mass was 13 m. In all cases, the diameter of piles was 0.5 m and distance between piles were 5 and 2.5 m in 4 and 9 states, respectively. The geotechnical parameters and model dimensions used, were selected according to the Fioravante & Girettis (2010) [1]. Sand and silica-sand with the defined properties were used for the soil mass and the interposed layer, respectively. Since Drucker-Prager criteria has better ability to express the behavior of coarse-grained soils, this criterion was used in the modeling [2]. The purpose of this study is to investigate the influence of interposed layer on bearing capacity and settlement of pile. Hence, because of simplifying the process of modeling, parameters of main soil and interposed layer are mostly similar. Piles and raft are made of concrete with an elasticity modulus of 21 GPa, Poisson's ratio of 0.2 and density of 2300 kg/m³. The crack growth analysis with the compressive stress-plastic strain was used to express the fracture behavior of concrete [2, 3 & 4]. In the present study, frictional and vertical contacts between surfaces were considered for conducting interactions between different materials. For frictional contact, the penalty formulation with the fixed friction coefficient of tanδ was used where δ is the angle of friction. The penalty formulations and hard contact were applied between two surfaces for the normal contact. Interactions were considered in the modeling including raft-soil mass, raft-interposed layer, pile-raft, interposed layer-soil mass, interposed layer-pile and the soil-pile [5 & 6]. Coefficient of soil lateral pressure used in this study corresponds to k₀=0.65 which is introduced in many geotechnical conditions [7]. A uniform distributed vertical load 500 kPa was applied on the raft. For getting results in every portion of loading time, this amount is applied in order of 5 kPa in each time interval. To accelerate the process of analysis and because of the symmetry of all models in two directions of X and Y, the quarter model technique was used, so that movements in the direction perpendicular to the sheet and rotation around perpendicular axes on the sheet were not allowed on the border of symmetry. The boundaries of the models due to the enough distance from the piles were considered in a way that lateral displacement and rotation around the vertical axis was not allowed. Furthermore, the bottom of the soil mass was considered as complete fix due to the enough distance from the pile foot.
Conclusion
In this research, a numerical – parametric study is performed on special kind of piles named noncontact piles and results are compared with contact piles. Results of this study can be summarized as follows:
1. By increasing the number of piles from 1 to 9, the settlement reduced more in a noncontact state showing more effectiveness of implementing 9 contactpiles and thus requiring more piles in this case.
2. Soil surface stress differences in noncontacts states from 4 to 9 piles was less than contact state (approximately 1/7) indicating that more piles is needed to conduct the contact state.
3. Stress changes in the soil under the pile in noncontact state by adding piles from 1 to 4 was higher than adding piles from 4 to 9 indicating the suitability of using 4 noncontact piles; while, in the contact state, the stress changes in the soil under the pile in both cases from 1 to 4 piles and from 4 to 9 piles was noteworthy showing the necessity of using the ninth pile.
4. Unlike the states of 4 and 9 piles, the negative friction in noncontact state and 1 pile was seen along the piles, which can be due to the fewer piles and the effect of interposed layer density as well as soil mass at greater depthsbecause of lesser effect of piles in load-bearing.
5. The ratio of heads load in the contact to the noncontact piles was about 2.5 to 4 reflecting the positive impact of using interposed layer on load reduction and smaller cross-layer design for piles. In addition, the ratio of heads load in the contact to the noncontact piles was higher for 4 piles than 9 piles that represented the suitability of using 4 piles.
6. Based on the results of geometric parametric studies it is found that:
(A) By resizing the elements from 0.25 to 0.5 m, the results had not changed and only time of analysis was increased.
(B) Among three values of 0.5, 1 and 1.5 m for interposed layer thicknesses, the thickness of 1 m was enough and had a good effect on the stress distribution and involving shallow soil in bearing vertical stress.
(C) The raft thickness of 1.6 m was appropriate so that with this thickness, the resultant effect of increasing vertical loads (raft weight) and increased rigidity due to increased raft thickness caused the stress and settlements remain in a reasonable range.
(D) Due to the increased friction by increasing in diameter, the optimal diameter of 0.5 m was achieved for piles which reduced the settlement by receiving more load.
(E) Among three pile lengths of 10, 19 and 25 m, the optimal length was 19 m; so that by further increase in the length, stresses and settlements were not noticeably changed.In total, noncontact piles had better performance compared to contact piles in similar conditions.
Reference
1. Fioravante V., Giretti D., "Contact versus noncontact piled raft foundations", Can. Geotech. J. 47 (2010) 1271-1287.
2. Saba H., "Verification of nonlinear condition of anchored walls in various loading", Thesis document of Amirkabir University of Tehran, Iran (2003).
3. Fioravante V., "Load transfer from a raft to a pile with an interposed layer", Geotechnique 61, No. 2 (2011) 121-132.
4. Dastani H., Shariati M., "Numerical and experimental analysis of controlling of crack propagation route in a plane under cyclic uniaxial loading by creating openness", Thesis document of Shahrood Industrial University of Shahrood, Iran (2014).
5. Randolph M. F., Wroth C. P., "Application of the failure state in undrained simple shear to the shaft capacity of driven piles", Geotechnique, Vol. 31, 1 (1981) 143-157.
6. Poulos H. G., Small J. C., Ta L. D., Sinha J., Chen L., "Comparison of some methods for analysis of piled rafts", Proc. 14th Int. Conf. Soil Mech. Found. Engng, Hamburg, Balkema, Rotterdam, Vol. 2 (1997) 1119-1124.
7. mottaghi A., "3D static and dynamic analysis of pile group with considering soil-pile interaction", 6th National Congress of Civil Engineering, Iran, Semnan (2012).

---

In many rock engineering projects, accurate identification of rock strength properties is very important. Uniaxial compressive strength is one of the most important features to describe the resistive behavior of rocks which is used as an important parameter in the design of structures especially underground openings. Determination of this parameter using direct methods, including uniaxial compressive strength tests is costly and time-consuming, and also sometimes preparation of standard samples in many rocks is difficult. In such cases, the implementation of some simple and non-destructive tests and using empirical relations can increase the evaluation speed and reduce costs. These relations even regional or local (For example within a geological formation or a single lithology) can help in the estimation of these parameters in order to be used in geotechnical projects. In this study, samples of existing limestones in south west of Tehran (Capital of Iran) were prepared and uniaxial compressive strength, point load, Schmidt hammer and Shear wave velocity tests on which have been performed. Then by the statistical evaluations of the results, the empirical relations between uniaxial compressive strength and the results of other tests are obtained. The comparison between the predicted and observed values of uniaxial compressive strength represents the validity of obtained empirical relations. The application of the proposed relations for limestones in the study area and those with similar geological conditions will provide acceptable results.

---

Introduction
In many areas of the world, the mechanical properties of soils for utilization of land are not sufficient. For improvement of these lands, soil stabilization such as compacting, installation of nails, elders of piles, mixing soil with lime or cement before or during constructions on the surface or inside of the ground can be useful. Microbially induced carbonate precipitation (MICP), due to its versatility and stable performance, has been recently attracted the attention of many researchers in the field of the geotechnical engineering around the world. MICP is a biological technique that is naturally caused to create a cementation agent, which is known as calcium carbonate or calcite by controlling the metabolism of bacteria. Although there are many biological processes that can be lead to MICP, but the using of urea hydrolysis by bacteria is commonly used more. In this method, aerobic bacteria with the enriched urease enzymes inject into the soil. Hydrolysis of urea occurs when the bacteria speeds up the hydrolysis reaction to produce ammonium and carbonate ions. In the presence of soluble calcium ions, carbonate ions are precipitated and formed the calcium carbonate crystals. When these crystals are formed on a grain of soil or like a bridge between them, they prevent the movement of grains and thus improve the mechanical and geotechnical properties of the soil.
Material and methods
In the present study, the effect of increasing fines on the improvement of Anzali sandy soil, and soil resistance parameters for improving the clean sand and its mixtures with a fine grained cohesive soil and a fine grained cohesionless soil separately in a percentage weight of 30 by MICP and using a small scale of direct shear test (6×6) have been investigated. In the present study the sandy soil was collected from the coast of Bandar Anzali Free Zone and for the preparation of samples of clayey sand and silty sand, Kaolinite clay soils and Firouzkooh broken silt were used, respectively. Anzali sand is poorly graded and had a rounded corner with an average particle size of 0.2 mm, somewhat, sharpening cores are also found in its granulation. In addition, its fine grained content is very small (less than 1%). The Kaolinite clay is also labeled with a liquid limit of 40, a plastic limit of 25, and a plasticity index of 15 as an inorganic clay (CL). The used microorganism in this study is urease positive Sporosarcina pasteurii, which is maintained with the number of PTCC1645 at the Center Collective of Industrial Microorganisms of Iran Scientific and Research Organization. The bacterium was cultured in a culture medium containing 20 g/l yeast extract and 10 g/l ammonium chloride at pH 9 under aerobic conditions in incubator shaker machine at 150 rpm and temperature of 30 °C. The organism was grown to late exponential/early stationary phase and stored at 4 °C before injection in samples. A solution of calcium chloride and urea with a molar ratio of one is also used as a cementation solution. With the direct shear test (6cm×6cm) as a benchmarking of the shear strength in the before and after improvement steps, molds fitted with a shear box made of the galvanized sheet with a thickness of 0.6 mm and it consists of two main parts, the body,  in the middle of which an exhaust pipe was embedded in the injector waste fluid. At the bottom of the samples, a layer of filter paper was placed in order to prevent soil washes, and then all samples with a thickness of 2 cm, with a relative density of 30% at the same weight and height were pressed. In the upper part of the samples, a layer of filter paper is similarly used to prevent the discontinuity of soil particles when injected biological materials are used. Biological solutions are injected from the top to the specimens and allowed to penetrate under the influence of gravitational and capillary forces in the sample and discharge the inhaled fluid from the exhaust pipe. The criterion for determining the volume of the solution to inject into each sample is the pure volume (PV) of soil. The preparation process of the samples was initiated by injection of a PV water unit, followed by a two-layer mixture of bacterial suspensions and cementation solutions, each with a volume of one PV, and then for biological reactions, 24 hours to the sample at laboratory temperature (25 ± 2) is given. After the time of incubation, the solution of cementation is injected into the sample for a period of three days and every 24 hours. The processing time of samples is also considered 28 days. In this study, optical density (OD) was selected as a benchmark for estimating the concentration of bacterial cells in the culture medium, and in all stages of development, and precisely before injection of bacteria suspension into soil samples, it was measured by a spectrophotometer device at 600 nm (OD₆₀₀) wavelength, which was obtained for all bacterial suspensions in the range of 1.7 to 2 before the injection. To determine the activity of urea bacteria, 1 ml of bacterial suspension was added to nine milliliters of 1.11 molar urea solution, and by immersing the electrode of the electrical conductivity in the solution, its conductivity was recorded for 5 minutes at 20 ± 2 ° C. The rate of urea activity in the pre-treatment stage for all specimens was in the range of 0.8 to 1.23 mS min^-1. In order to evaluate the shear strength parameters of soil samples, before and after the improvement operations, a direct shear test was used based on the ASTM D3080 standard. This test was performed for all samples under stresses of 50, 100 and 150 kPa in undrained conditions at a loading speed of 1 mm/min up to a strain of 15%. Also, samples of soil with a moisture content of 7% and a relative density of 30% (as already mentioned) have been restored. SEM analysis was carried out to determine the distribution of sediment between soil particles and EDX analysis in order to identify carbonate calcium sediment formation elements in improved soil samples, by scanning electron microscopy on Anzali sandy soil samples in before and after improvement conditions.
 
 
Conclusions
The effect of the increasing cohesive and cohesionless fines on the bio-treated process of sandy soil is the main subject of this research. For this purpose, three samples of clean sand, sand containing 30% clay and sand mixture with 30% silt in a relative density of 30% were treated with MICP method and their shear strength parameters were evaluated by direct shear test after 28 days of processing. Using the direct shear test and analyses of SEM and EDX data, the results are represented as below:
1.  The microbial sediment of carbonate calcium has greatly improved the resistance properties of all three soil samples.
2. A sample of clayey sand, in spite of a higher improvement compared to the other samples with an average shear strength of 113.7% in comparison to to its untreated state, it has the lowest shear strength among the three improved samples.
3. Increasing the clay content of 30% increases the soil voids. On the other hand, it reduces the friction angle and shear strength of the soil in the pre-treated state and also facilitates easier movement of the bacteria between the pores in the soil. More favorable distribution of sediment calcium carbonate was occurred and, as a result, increased adhesion between soil particles.
4. The increase of cohesionless fine particles creates more bonding points between sand particles and, therefore, calcium carbonate crystals form shorter distances between the soil bridges. As a result, with the end of the improvement process, the shear strength parameters of the sandy soil containing 30% of the silt compared to the clean sand have a higher value.
5. SEM images of the clean sand in both before and after improvement show that the calcium carbonate precipitation occurred with a uniform and thin layer that surrounds sand grains and another part of the sediments formed in the joint of grains.
6. Cube-shaped crystalline sediments confirm that the sediment formed in the soil is a stable type of calcite and that the relative increase in the friction angle of the improvement samples can be attributed to solid particles and multifaceted sediments. Also, the elements of carbon, oxygen, and calcium, which are the main components for the formation of calcium carbonate deposits, have been found in the EDX analysis of improvement sand samples../files/site1/files/124/7sohrabi%DA%86%DA%A9%DB%8C%D8%AF%D9%87.pdf
 

---

Introduction
In general, landslides, in particular, earthquake-induced landslides, are among the phenomena that have caused great damages in recent years in Iran and the world. Although many studies have been done on the identification and description of landslides in general, the study of landslides caused by the earthquake, especially in Iran, is at the beginning stages. In a few studies, some landslides and some of their characteristics have been introduced. A magnitude 7.7 earthquake occurred in the Guilan Province was occurred on May 31, 1990. This earthquake is one of the most important earthquakes in Iran history due to its magnitude and occurrence of landslides. In various studies, the most important landslides have been listed. The development of quantitative and qualitative studies on earthquakes that have caused many landslides (such as the Manjil, Avaj, Firoozabad, Kojur, Sarein and Ahar and Varzaghan earthquakes) increase our understanding of natural disasters and, consequently, the management of the dangers resulting from them. The purpose of this research is to identify the factors affecting the occurrence of landslides caused by earthquakes, to determine the impact of each on the occurrence of this phenomenon, and also to prepare a map of earthquake hazard zonation hazard by utilizing the methods used in this research. In this study, hierarchical analysis method has been used to prioritize the factors affecting the occurrence of landslide and also the zoning of earthquake landslide hazard in the study area.
Research Methodology
The study area is located between 49˚ 30ꞌ and 49 ꞌ45˚ and latitudes 36º 00 ꞌ 45" and 36º 30ꞌ 52" with a surface area of ​​309.30 km². In this research, in order to zoning the earthquake-induced landslides hazard, in addition to providing a map of landslides, seven factors influencing the occurrence of this phenomenon were identified and examined. These factors included elevation, slope, arias intensity, friction angle, adhesion, curvature of the slope and aspect. In this research, Analytic Hierarchy Process (AHP) method, one of the multi-criteria decision making models, was used with two approaches to using expert knowledge and data and expert knowledge together to prioritize the factors influencing the occurrence of landslide. Finally, two landslide hazard zonation maps were prepared. In a hierarchical analysis method related to the expert judgment, it was used to determine the priority of different criteria and sub-criteria and convert them into small amounts of oral judgments (expert opinion) based on the pair comparison, in which the decision maker preferred the factor in relation to other factors using the relevant tables, these judgments are converted into small amounts. In the method of using data and expert judgment simultaneously, first, in order to determine the priority of criteria from oral judgments (collection of expert opinions), we used to determine the importance or weight (W_i) of each sub-criterion (R) is also used to link the landslide area to each class and landslide area in the region.
Results
The results obtained from the paired comparison of the effective factors in the occurrence of landslide show that the relative preference of the factors include the factor of arias intensity, friction angle, slope, adhesion, aspect, height and curvature of the amplitude. The greatest influence on the sub-criteria for the sub-criteria is 10-11.54, which is related to the arias factor and also the lowest effect for the sub-standard of the domain curvature factor. Also, according to the zoning maps, in the first model, 73% and in the second model, 57% of the surface area are very high and very high risk areas, which indicates the high sensitivity of the study area to the earthquake-induced earthquake phenomenon.
According to the results obtained from the verification and evaluation of the models and comparison of the mapped data with the hierarchical analysis method (using expert knowledge and data) and a method that uses only expert knowledge, the map is derived from a method where bundles of knowledge and data are used simultaneously, in order to weigh the parameters, it is more in line with the map of the landing list of the region.
Conclusion
According to the results obtained from the review and evaluation of the two models in a method in which knowledge and data were used together, the Q_S value was 0.40 and the accuracy of the method (P) was 0.016. However, in a method in which only the expert judgment used to weigh the criteria and sub-criteria, the sum of the quality and accuracy of the method were calculated to be 0.37 and 0.006, respectively. Hierarchical analysis method, in which the benchmarks and sub-criteria of benchmark knowledge and data are used together, have a better performance than the other model, and the results are closer to reality. In addition, it also works better in distinguishing between high and high risk areas../files/site1/files/124/5rajabi%DA%86%DA%A9%DB%8C%D8%AF%D9%87.pdf

---

Introduction
Concrete faced rockfill dams have been considered in recent years more than other types of dams due to their low dependency on the bed and the shape of the valley, as well as the simpler construction technology. In this regard, rockfill dams are a suitable substitute for embankment dams because of higher stability of the body and the availability of rock aggregates. On the other hand, because the permeability of rock aggregates is much higher than other materials, different methods are used to seal these types of dams. One of these methods is the use of non-impermeable concrete facing in the upstream of these dams. This particular type of gravel dams is called Concrete-Faced Rockfill Dams (CRFD). In this study, a contact element with a definition of elastic-plastic failure in the modeling process is proposed to simulate the surface separation and re-contact of the concrete face with the rockfill surface of the dam.
Method
In this paper, behavior of a concrete faced rockfill dam under earthquake loads is investigated. For this purpose, near-field earthquake records with focal depth lower than 15 km (for example Tabas earthquake 1978, M=7.4, and San Fernando earthquake 1970, M=6.6) are used. Moreover, to study the dam behavior under dynamic loads, interaction between concrete face and rockfill part of the dam is investigated and finally, some parameters including displacement, absorbed energy and base shear are evaluated. So, finite element method and Abaqus software is used for the study. Verification of the models is carried out using the results of previous researches by conducting modal analysis and determining natural vibration period. Then, the interaction between the concrete face and rockfill part as well as the effect of water level changes in stability of dam under dynamic load is investigated. Concrete behavior is simulated using concrete damaged plasticity. Therefore, concrete density, compressive strength and tensile strength and elasticity modulus are 2350 kg/m³, 25 MPa, 3 MPa and 29 GPa, respectively. Poisson’s ratio is assumed to be 0.2. Furthermore, 4-node shell elements are used to simulate concrete face and Drucker-Prager constitutive model is used to define rockfill material behavior.
The density and Poisson’s ratio for 2B, 3C and 3B layers are 2150 kg/m³ and 0.35, respectively. The shear modulus values for these layers are respectively 8.93, 2.89, and 3.85 GPa. In order to perform the simulation, the part of the dam structure beside the bed rock and the surrounding rock is considered as fixed bearing, and only the rockfill part and concrete face of the dam is simulated. Based on this assumption that the bed is rigid, there is no need to consider the dam foundation. This method is frequently used in literature review.
All the surfaces of the dam and bed rock are considered as fixed bearing to simulate the real condition where the dam is attached to bed rock and the surrounding rock. The interaction between dam layers is defined as tie. For defining the interaction between rockfill body and concrete face, tangential and normal contacts are defined using penalty method with friction coefficient equal to 0.5. In the next step, the model is meshed using 4-node shell elements for concrete face, 8-node brick and 4-node pyramid solid elements for rockfill body. Rayleigh damping is used to simulate the structure damping. The effective length of the dam reservoir has been determined by conducting several analyzes, so that the minimum required length for reservoir is reached in order to decrease the number of elements of the model.
Results and discussion
1. Interaction between concrete face and rockfill body
The results show that the increase of friction coefficient between concrete face and rockfill part from 0.5 to 0.7 has not affect the displacement of dam crown along the earthquake direction. However, when the concrete face is fixed to the rockfill part, significant changes are induced in dam crown displacement time history. In all cases, the deflection due to the dam weight is increased when the concrete face is attached to the rockfill body. The reason can be attributed to the tied interaction between these layers which results in similar deflection of concrete face with rockfill body and higher deflection of concrete dam crown. However, after the application of earthquake load, the displacement of the dam crown decreased in both analyses when tie interaction is defined between concrete face and rockfill body. In this study, due to the very high volume of analysis and its timeliness, it was not possible to examine the dam behavior in the free vibration regime, and therefore, it is not possible to assume the last displacement values at the end of analyses as the permanent displacement of dam. Figure 1 shows the relative displacement of the dam for the two selected earthquakes with a friction coefficient equal to 0.5 between the concrete face and the gravel body. According to Figure 1, the maximum displacement induced by the earthquake is related to Tabas and then, San Francisco earthquake. Furthermore, the high energy content of the Tabas record has been more effective in inducing greater displacement than the other record.
 
Figure 1. Lateral displacement of dam crown relative to the dam base for the selected earthquakes; Tabas and San Fernando.
The results also indicate that when the friction coefficient between concrete face and rockfill body is 0.5, the lowest damage occurs in the dam compared to that happens when friction coefficient is 0.7 or when the surfaces are tied. When the tied surfaces are used, the most damages takes place in concrete face, since all rockfill body displacement transmits to concrete face which results in much more concrete damages compared to the other interaction cases.
2. Effect of water level in reservoir on dam behavior
In this section, the effect of water level on seismic behavior of dam is investigated. For this purpose, the dam reservoir is analyzed in three cases including empty, half full and full (90% of dam height). Each study cases are examined under San Fernando and Tabas earthquakes. Figure 2 shows the relative displacement of dam crown in the three water level case for San Fernando and Tabas earthquakes.
 
Figure 2. Relative displacement of dam crown in three water level cases of empty, half and full for (a) Tabas and (b) San Fernando earthquakes
According to Figure 2, for both earthquakes, the dam crown displacement along the earthquake direction is significantly increased by increasing the water level, so that the maximum displacement in full case is 50% higher than empty case.
Conclusion
In this study, using the finite element method and simulation by Abaqus, the seismic behavior of concrete face rockfill dams has been investigated. For this purpose, the verification is firstly carried out using previous research results in literature. In the next step, nonlinear dynamical analysis is carried out, taking into account large displacements for the models under the earthquake record acceleration. The results illustrate that increasing the friction coefficient between the concrete face and the rockfill body from 0.5 to 0.7 has no significant effect on the displacement of the dam crown under earthquake load. Moreover, by using tie interaction between the concrete layer and the rockfill body, there is a substantial difference in the history of the relative displacement of the dam, and the displacement of the dam due to its weight has been increased. Furthermore, the results of this study exhibit that, with increasing the water level in dam reservoir, the deformation of the crown of the dam along the earthquake application direction has had a relatively significant increase, such that in the full state, the maximum displacement is increased by about 50% compared to that of the empty case. This is while the most damage of concrete is observed in the case when half height of dam in filled by water. Due to the more destructive power of near-field earthquakes and their impact nature, only near-fault earthquakes have been used in this research. Therefore, the results of this study are valid only for the behavior of dam under near-field earthquakes.
./files/site1/files/134/3.pdf

---

Introduction
Improving the geotechnical characteristics of soils including superficial or deep soils has always been a challenge to geotechnical engineers. Therefore, various physical and chemical methods are used to improve different types of soils. In general, any physical, chemical, biological or combination of methods are used to change the characteristics of natural soil mass in order to achieve engineering goals which is defined in the "soil stabilization." Among different types of additives for soil stabilization, the use of pozzolans has been investigated by researchers because of their chemical compatibility with the environment and the cementation products due to chemical reactions. Todays, a lot of researches has been done on the use of natural or artificial zeolites as pozzolanic materials for the production of cement mixtures. This material, as a pozzolan, increases the speed of the pozzolanic reactions and reduces the density of cement products. However, many studies have been done to investigate the effect of zeolite and sepiolite on the resistance of cement products such as concrete, but so far, the use of these additives has been less considered for soil improvement. On the other hand, because of the compatibility of zeolite and sepiolite with the environment and their unique physiochemical properties, it is necessary to pay attention to these additives in order to improve the soil. Therefore, in this research, the effect of zeolite and sepillot additives with different percentages at different treatment times have been investigated to determine the elasticity modulus and hydraulic conductivity with focus on soil microstructure behavior.
Materials and methods
1. The properties of the soils
In this research, two types of soil including clayey sand (with 20% clay) and sandy clay (with 51% of clay) were used. The studied soils were a mixture of clay and sand of Firoozkouh (a typical type of sand located in north of Iran). Some physiochemical properties of zeolite and sepiolite are presented in Table 1.
Table 1. Physiochemical properties of zeolite and sepiolite used in this study

L.O.I.	Na2O	K2O	MgO	CaO	Fe2O3	Al2O3	SiO2		Chemical component
25.11	0.02	0.01	15.73	0.01	o.61	0.3	55.3		Sepiolite (%)S
11.94	0.13	-	0.87	2.45	1.26	13.54	69.74		Zeolite (%)

2. Experiments
The uniaxial compressive strength tests were performed at 0.1 mm/min according to ASTM D2166 standard. The stabilized soil samples were compacted at percentages of 0, 5, 10, 15, 20 and 25 in cylindrical molds (38mm × 76mm) in five layers to achieve the desired density. In order to investigate the effect of curing time, the samples were placed inside sealed containers and underwent the test at instantaneous, 7, 14, and 28 days and at the desired additive percentages. To investigate the effect of additives on the soil hydraulic conductivity, clayey sand soil with additives 5, 10, 15, 20, and 25% was prepared using dry mixing method. Then, the prepared mixture was poured from a specific height into the permeability mold with a height of 8.65 cm and diameter of 5 cm. In this way, the specific dry unit weight of all samples was obtained as 1.47 g/cm³, close to the minimum specific dry unit weight. In this research, concerning the considerable effect of fine-grained soils on hydraulic conductivity, falling head test was used to determine the permeability coefficient.
In order to the morphology of the clayey sand soil without additives and stabilized with additives 15% was examined through SEM test.
Discussion and results
1. Modulus of elasticity
In this study, after uniaxial tests in different percentages and ages, the stress-strain graphs were plotted and then the elasticity modulus was calculated. The results showed that, with increasing zeolite content, the modulus of elasticity has been increased and, with increasing curing time, except for a slight decrease, after 7 days, the modulus of elasticity increased. During the initial treatment (7 days), the hardness of the sandy clay soil decreased and then increased with increasing time. In general, hardness in both soils in the high percentages of zeolite is significantly is increased.
Also, the effect of sepiolite on the modulus of elasticity has been studied. The results indicate that with the increase in the percentage of additive and lengthening the curing time, the modulus of elasticity is increased. This increase in the stabilization of both sandy clay and clayey sand soil is almost the same. In addition, in the case of sepiolite modification, the elasticity of sandy clay and clayey sand is approximately equal to 5 times in comparison to the initial value of unstabilized soil. However, in zeolite, the modulus of elasticity in clayey sand soils is almost 2 times, and sandy clay is nearly 5 times higher.
2. Permeability
To investigate the effect of additives on the soil hydraulic conductivity, clayey sand soil with additives 5, 10, 15, 20, and 25% was prepared using dry mixing method. The samples were saturated in a short period and permeability test was carried out immediately. Permeability coefficient changes were mostly influenced by physical factors. Therefore, due to the fineness of both types of additives, the hydraulic conductivity decreases with increasing additive content. The amount of reduced hydraulic conductivity in sepiolite stabilization is greater than zeolite due to the structure of the sepiolite (fiber-shaped) compared to zeolite.
 
 
3. SEM imaging
In this study, attempts were made to examine the reasons behind the obtained results more carefully through SEM imaging.

c                                     b                              a
Figure 1. SEM image of non-stabilized clayey sand soil (a) soil stabilized with zeolite 15% (b) soil stabilized with sepiolite 15% (c) during the curing time of 28 days at magnifications 10000X
Figure 1a displays the SEM image of non-stabilized clayey sand soil. As can be seen in the figure, the soil structure is clear as layered and clay scales can be seen as laminated. Figure 1b demonstrates the SEM images of clayey sand soil stabilized with zeolite 15% during the curing time of 28 days. The sample has lost its layered structure in response to stabilization with zeolite during the curing time and changed into an integrated structure. This can be due to incidence of chemical reactions such as ion exchange and pozzolanic reactions in response to adding zeolite. Figure 1c demonstrates the SEM images of clayey sand soil stabilized with sepiolite 15% during the curing time of 28 days. As shown in the figure, the sepiolite has a fibrous-shaped structure that is longitudinally twisted. Also, with  curing time increase, complex structures have emerged that could be due to the occurrence of chemical reactions.
Conclusion
This study examined the effect of zeolite and sepiolite additives on strength parameter of clayey soils. Accordingly, uniaxial compressive strength test was performed on clayey sand and sandy clay soil at percentages of 0, 5, 10, 15, 20 and 25% of zeolite and sepiolite with instantaneous curing times of 7, 14 and 28 days. Further, permeability test was conducted at different percentages on stabilized clayey sand soil. Also, to investigate the effect of these materials on soil microstructure, SEM imaging was performed at 28 days. The results show that both additives increase the elastic modulus of clayey sand and sandy clay soils. Also, the results indicate a steady increase in the stiffness of the cured soil with sepiolite during processing time. However, reducing soil hardness can be seen in stabilizing with zeolite at lower rates and lower percentages. In permeability test, hydraulic conductivity decreases with increasing additive content. The rate of permeability reduction in sepiolite is higher than zeolite. SEM images show that chemical reactions create an integrated structure that ultimately increases uniaxial compressive strength and modulus of elasticity. Also, SEM imaging depicts physical changes along chemical reaction in soil stabilized with sepiolite. Ultimately, increasing soil strength resulting from additive alongside environmentally friendliness is recommended in superficial and deep improvement of soil../files/site1/files/144/Rajabi.pdf
 

---

Introduction
Excavation in urban areas occasionally is accompanied by the improper performance of the support system for even small deformations. In this regard, deformation control design based on force-based approaches provides a more realistic reprehensive of excavation performance. Top-down deep excavation techniques are among the modern excavation stabilization methods in urban areas. In this method, unlike the conventional methods, it is possible to perform the excavation and construction operations simultaneously. The present study aims to investigate excavation stabilization using the main structure through the top-down approach. For this purpose, field and numerical evaluations of the stabilized project were conducted based on the top-down approach in the downtown of Qom city, Iran. This research reports the information obtained through monitoring and modeling using the finite element ABAQUS software, predicting the occurred deformations until the end of excavation operations using the calibrated model, and offering an initial estimation of the required stiffness for the support system with respect to the lateral deformations in four sites proposed, according to the studies of Line A Qom Subway.
Project specifications
Based on the geological studies of Line A Qom Subway Tunnel, the geological layers are classified into four soil classes. Qc-1 consists of gravely sand with fine content of 5 to 20%; Qc-2 is silty and clayey sand with fine content of 35 to 60%; Qf-1 is clayey silt with fine content of 60%; and Qf-2 is a silty clay layer with fine content above 60%. Line A of Qom subway passes the study area of the present study, which is located in Ammar e Yaser Street (Station A6). Based on the geotechnical studies of the project site, the site in the levels near the ground consists of Qc-2 but in the lower elevations, it is composed of Qc-1 and Qf-2.
Salam Trade Complex, located in the downtown of Qom city, has 6 underground stories and 6 above-ground stories. It is limited to the main street in the south and to urban decay in the three other directions. The final excavation depth, length, and width is -21, 36, and 32-52 m, respectively. The project structure consists of a steel moment frame with a retaining wall in the negative elevations and metal deck frame for ceiling construction. In this project, excavation wall deformation was monitored in three important sections (A, B, and C). Due to the vicinity to urban decay, a total station TS02 was used for monitoring these sections. According to the field surveys, the maximum horizontal deformation of the walls in sections A, B, and C is 24.10, 42.16, and 47.21 mm, respectively, which were measured in the 0, -1.5, and 0 m elevations.
Monitoring process and numerical simulation
To calibrate the prepared model, a sensitivity analysis was performed on geotechnical parameters including modulus of elasticity (E), internal friction angle (φ), and cohesion (C) of the layers by simulating 60 numerical models. Based on the sensitivity analysis results, an increase in internal friction angle and elasticity modulus for layer 1 (i.e., φ1 and E1) and elasticity modulus of layer 3 (E3) results in a decrease in lateral deformation. Finally, using the sensitivity analysis results and after several trials and errors, the numerical models for sections B and C were calibrated when reaching the depths of -8 and -11 m, respectively. Using these models, then, it is possible to predict deformations up to the end of the project.
To determine the required stiffness for the excavation support system, regarding the acceptable deformation of the adjacent soil mass, 160 numerical models were built and their results were analyzed. Based on the results of Brason and Zapata (2012), relative stiffens (R) were used to develop a relationship between the maximum lateral deformation of the wall and the required stiffness of the support system. R is a dimensionless parameter that represents the stiffness of a solid support system; the greater this value is, the more flexible the system would be. In this study, caisson pile length, excavation width, and buried depth of the wall were used for determining the R.
R =                                                        (1)
Figure 2 presents the maximum occurred deformation in terms of depth versus the relative stiffness for sites QC and QF.

Figure 2. Maximum deformation in terms of depth versus the relative stiffness for sites QC and QF
Conclusion

1. According to the monitory data, the maximum lateral deformation in sections B and C until the end of the project was 42.16 and 47.2 mm, respectively. Moreover, the deformation of the other points inside the excavation was 30 mm.
2. Considering the occurrence of maximum lateral deformations in the higher elevations in the monitored sections, it is inferred that excavation support at the ground level plays a key role in this approach. Hence, the lack of completing the structural frames and slabs for facilitating the excavation operation can lead to an increase in deformation levels.
3. Based on the prepared graphs, the top-down approach in sites QC-2 and QF-2, compared to sites QF-1 and QC-1, provides a more desirable performance for deformation control.

---

---