Journal of Engineering Geology

---

The Bakhtiari dam has located on Bakhtiari River in province of Lorestan. In order to access the crest of the dam, the excavation of a spiral tunnel is being studied. There are other access tunnels which are branched from this tunnel in different levels and are connected to grout galleries. According to the fact that this tunnel will also be used during the operation of the dam,The correct determination of mechanical parameters of rock masses for tunnel design and stability Analysis is very important. In order to analyse the stability of the underground rock structures, the mechanical and engineering parameters of the rock mass must be known. Accurate rock mass properties can only be obtained from large in situ tests. Such tests are seldom carried out as they are very expensive and time consuming. Sensitivity analysis of parameters can be applied for the optimisation of testing schemes. Sensitivity analysis helps to avoid mistakes due to subjective conjecture. In this article, after the introduction of regional geology and determination critical section on the tunnel path, the mechanical parameters of the rock mass surrounding the tunnel are modelled and analyzed by using FLAC3D software (numerical finite difference method). Parameters conducted in the analysis include the elasticity modulus (E), cohesion of the rock mass (C), friction angle (ϕ), coefficient of lateral stress (K) and tensile strength (&sigmat). Ultimately, according to the result of numerical modelling and parametric analysis, parameters affecting the stability are prioritized. The result of analysis showed that in this project, tensile strength of the rock mass does not affect the stability of the tunnel, and Also, in order of priority, E, ϕ, C, k parameters are important in design. The amount of field tests for rock parameters can be rationalised according to their sensitivity factors.

---

Evaluation of the excavation-induced ground movements is an important design aspect of supporting system in urban areas. This evaluation process is more critical to the old buildings or sensitive structures which exist in the excavation-affected zone. Frame distortion and crack generation are predictor, of building damage resulted from excavation-induced ground movements, which pose challenges to projects involving deep excavations. Geological and geotechnical conditions of excavation area have significant effects on excavation-induced ground movements and the related damages. In some cases, excavation area may be located in the jointed or weathered rocks. Under such conditions, the geological properties of supported ground become more noticeable due to the discontinuities and anisotropic effects. This paper is aimed to study the performance of excavation walls supported by nails in jointed rocks medium. The performance of nailed wall is investigated based on evaluating the excavation-induced ground movements and damage levels of structures in the excavation-affected zone. For this purpose, a set of calibrated 2D finite element models are developed by taking into account the nail-rock-structure interactions, the anisotropic properties of jointed rock, and the staged construction process using ABAQUS software. The results highlight the effects of different parameters such as joint inclinations, anisotropy of rocks and nail inclinations on deformation parameters of excavation wall supported by nails, and induced damage in the structures adjacent to the excavation area. The results also show the relationship between excavation-induced deformation and the level of damage in the adjacent structure.

---

./files/site1/files/7Extended_Abstract.pdfExtended Abstract
 (Paper pages 135-156)
Introduction
Many civil structures (e.g. tunnel walls, bridge pillars, dam abutments and road foundations) are subjected to both static and dynamic loads. Cyclic loading leads to occurring fatigue phenomenon. Fatigue is the tendency of materials to break, or the process of damage accumulation, under cyclic loading. It was found that the dynamic fatigue strength can be reduced by 30-70 percent on average compared to uniaxial compression strength. Different materials show different response when they are subjected to cyclic loading. Some materials become stronger and more ductile, while others become weaker and more brittle. Although it is clear that the mechanical properties of rock under dynamic loads varied dramatically from those under static loads, the nature of dynamic failure in rock remains unclear, especially in cyclic loading condition. Fatigue behavior of rocks was rarely studied in respect to other materials such as steel and soil. The performed researches on fatigue behavior of rocks indicated that fatigue life will be decreased by increasing load amplitude in logarithmic and exponentially pattern. Also, strain softening is the dominated behavior of rocks against cyclic loading. Furthermore, some parameters such as maximum load level, confining pressures, amplitude, and loading frequency have considerable effects on fatigue behavior of rocks. However, available data on fatigue behavior remain insufficient for solving the practical tasks of predicting rock bursts and earthquakes. Obtained results are inconclusive and sometimes discordant. The aim of the current work was to assess tonalite rock fatigue behaviour under different loading conditions to describe the fatigue damage process of the granitic rock.
Material and methods
Several core samples were prepared to perform this research. The core samples were prepared with a L/D ratio of 2.5 with an average diameter of 54 mm. Before the fatigue tests, the physical and mechanical properties of the rocks were measured. Uniaxial compressive strength test (UCS) has been done on 5 core samples. The tests were performed in the load-control mode with a 1.6 kN/s loading rate. The tests were conducted to obtain the physico-mechanical parameters of the rocks in static loading condition, and provided a reference for subsequent dynamic tests. The cyclic tests were performed in both load and displacement control modes. To record axial and lateral strains during the fatigue tests, four strain gauges have been employed with arrangement of two axial and two laterals. Also, three acoustic emission sensors were installed on top, mean and bottom of the core samples to record cracking sound. In order to doing the tests a servocontrol Instron machine with 500 kN capacity was employed. The fatigue tests were conducted with three different maximum loads, 1 Hz frequency, and constant amplitude (0.82 of uniaxial compressive strength). The maximum stress level (the ratio of maximum cyclic stress to static strength) was varied 0.80, 0.85, and 0.90. The amplitude level (the ratio of amplitude stress to static strength) ranged from 0.50 to 0.70 and 0.90. Finally, Multi stages loading with increasing amplitude were applied for the displacement control tests. The results of fatigue tests have been evaluated by fatigue damage parameters including maximum and minimum axial strain, maximum and minimum lateral strain, tangent and secant modulus, toughness and hysteresis energy.
Results and discussion
The obtained results indicated that during fatigue process failure occurs below the maximum strength loading condition as a result of accumulative damage. Analysis of the fatigue test results showed that the fatigue failure consisted of three stages: fatigue crack formation (initiation phase I), stable crack propagation (uniform velocity phase II), and unstable crack propagation resulting in a sudden breakdown (accelerated phase III). By comparing the axial and lateral deformation, it was found that lateral deformation is more sensitive to fatigue. At higher stress levels, considerable part of fatigue life is response to crake development, whereas at lower stress levels, crack acceleration phase of fatigue life is distinguishable. Descending trend of loading and unloading tangent modulus shows a scatter pattern. This behavior may be related to the calculation method and loading condition, as well as microstructure and behavior of the rock mass. In spite of tangent modulus results, the three-stages of damage process (especially phase I and II) for secant modulus in both loading and unloading conditions are clear. The result is due to the method of calculation and increase in axial strain with increasing number of cycles. Brittle behavior of this type of rock leads acceleration phase to be hidden and unclear in most of fatigue damage parameters. A dramatic decrease of toughness and hysteresis energy in the first few cycles is due to the closing of pre-existing micro fractures. In fact, during the initial cycle, the rock behaves in a more ductile fashion than in the next few cycles. Thereafter, toughness begins to increase slowly, then steadily, and finally rapidly. A similar behavior was found for hysteresis energy as well. This fact indicated that cracks generated in parallel to loading direction. Fatigue displacement control tests show a strain softening behavior for the granitic rocks. This behavior is highlighted in variation of maximum stress during the tests. This parameter, especially in final step of loading, shows distinguishable decreasing trend.
Conclusion
The tonalite rocks were subjected to uniaxial cyclic loading in both load and displacement control mode. The following conclusions were drawn from this research.
-Accumulated fatigue damage occurs in an obvious three-stage process. This is the result of the micro-fracturing mechanism in the fatigue process.
-By comparing axial and lateral strain damages, it was found that crack propagation occurred in the loading direction and crack opening occurred in the lateral direction. So, among fatigue damage parameters, lateral strain shows the best three-stage fatigue damage behavior.
- Strain softening was found as rock response to cyclic displacement control loading.

---

In recent decades many researchers have studied on the damage assessment of structures after a seismic event. To assess the damage of structures under an earthquake, it is so important to study the correlations between earthquake parameters and damages of the structures. A lot of seismic parameters have been defined by researchers to characterize an earthquake. Spectral parameters of an earthquake convey a variety of information about ground motion, so they can properly characterize an earthquake. Also a lot of damage indices were proposed by researchers to quantify the damage of the structures or to rank their vulnerability relative to each other. Park-Ang index is one of the best indices to describe the damage of a structure. In this paper, the correlations between spectral parameters of earthquakes and Park-Ang indices are studied. Three RC frames with different height are analyzed under far-fault earthquake records by nonlinear dynamic analyses. The correlations between spectral parameters and Park-Ang indices of the frames are calculated. The results show that in all the frames most of spectral parameters have strong correlations with damage intensity. In order to estimate the damage potential of an earthquake, some spectral parameters which have high correlations with damage intensity can be proper indices. Housner intensity, acceleration spectrum intensity and velocity spectrum intensity are shown to have strong correlations with damage intensity. In this paper, a new spectral parameter which has high correlation with damage intensity is achieved. 

---

Earth-fill dams stability in steady state seepage condition is very important, especially during earthquakes. Numerical software analyses require accurate and realistic modeling of construction stages. Since earth-fill dams are constructed in different layers, so these conditions should be considered in software modeling to achieve a reasonable design. In this study, an earth-fill dam is modeled in PLAXIS software and the effects of the number and shape of layers are studied in dry and steady-state conditions. Obtained results in static and pseudo-static analyses show that modeling of earth-fill dams with different layers has significant effects on shear stresses and horizontal displacements. For example, horizontal displacements and shear stresses, increase at least 50% and 17% respectively, in comparison with single layer models. According to the obtained results, it can be mentioned that modeling of an earth-fill dam in the layered model and rather in inclined layers are more reasonable

---

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
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
One of the most sensitive and important issues in some civil engineering projects is slope design and application. The process of slope design always involve many uncertainties. Hence, it is impossible to accurately comment on its stability or instability. Most of the uncertainties in the slope stability analysis are related to the nature of materials, geometry, environmental conditions, model errors, and measuring errors as well. Therefore, the slope stability analysis with a deterministic approach which uses the concept of safety factor would often not result satisfactory. Consequently, the use of probabilistic methods is more advised. Accordingly, in recent years, the probability analysis has been used to slope stability analysis. In these analyses, the effective quantities of slope stability are considered as statistical distributions, and the reliability coefficient would then be a statistical distribution. Likewise, one of the approaches to simulate uncertainties in the probabilistic analysis is to use the variation coefficient. If the variation coefficient changes, the probability of failure will change accordingly. When the variation coefficient becomes a larger number, costly solutions are required to reduce the probability of failure. If the variation coefficient becomes low, the reliability will be increased and the required costs to reduce the probability of failure will be decreased. Therefore, determining the amount of variation coefficient in these analyses is very important. Furthermore, the correlation coefficient between the quantities is another effective parameter in computing the probability of failure.
Material and methods
In this research, the stability analysis of the slope facing the spillway of the Shiraz Kavar dam has been done in two probabilistic and deterministic methods. Since circular slip probability is more likely than other types of failure, in the analysis of the stability of this slope, the problem of circular failure is very important, and an appropriate equilibrium program should be used for circular failure analysis. Therefore, SLIDE software was used to slope stability analysis. For material behavior, the Hook-Brown failure criterion was applied. In order to determine the strength parameters of the criterion, Geological Strength Index (GSI), uniaxial compressive strength (UCS) and rock constant parameter m_i were used. For crushed rock with a moderate quality of crushing, the GSI quality of the rock mass was about 23 to 38, which the average value of that for the rock mass of the overflow was assumed 35. Also, the uniaxial compressive strength of the rock was evaluated about 50 to 100 MPa with an average value of 75 MPa. In addition, the value of m_i was 10, and due to mechanized drilling, the disturbance factor was considered to be 1. The amount of unit weight was assumed to be 22 kN/m³. The initial model used for deterministic and probabilistic analyses, is the Morgenstern-Price model. To conduct probabilistic analyses, Monte Carlo simulation was performed using random sampling method (RS-MC) and 200,000 sampling were used to converge the simulation results. To determine the coefficient of variation and the probability distribution of UCS, GSI and m_i, the proposed values ​​of Hook (1998) were applied and for unit weight (γ) James Rodriguez and Sitar (2007) studies were used. Also, the minimum and maximum values ​​of UCS and GSI are determined based on the results of experiments, and Third Sigma rule was utilized for m_i and γ quantities. Since the earthquake phenomenon is rarely of great intensity and the number of small earthquakes is higher, therefore the truncated exponential distribution function can be in good agreement with the results of the earthquake. Usually, the maximum magnitude of the earthquake acceleration coefficient is twice that of the average.
Results and discussion
In the presented paper evaluation denotes that the safety factor computed by probabilistic analysis is given as a distribution function. The function provides a clearer view of failure condition. However, a deterministic analysis only illustrates a certain value for the failure. In addition, the results of the probabilistic analysis show that it is possible to optimize the dip of the slope; such that it remains completely stable and the volume of earthwork is also minimized. Therefore, by using probabilistic analysis, the optimal dip of the slope was determined. In these circumstances, the amount of earthwork was decreased by 28,000 cubic meters. Also, the sensitivity analysis of the variation coefficient and correlation coefficient between parameters are analyzed. The results of the sensitivity analysis of the failure probability versus the variation coefficient of the quantities showed that the quantities of sensitivity factor for static conditions is greater than the corresponding pseudo-static, and the GSI amount is the highest, while the specific gravity has the least effect on the probability value. In addition, the analysis indicated that if the GSI coefficient of over 21% is selected, the probability of a static failure is higher than the permissible limit. Also, increasing the variation coefficient of quantities by as much as 50% exhibits that the probability of static failure is still below the permissible limit. Also, the correlation coefficient between UCS and GSI shows that the higher variation coefficient of the quantities is chosen, the more variations of failure probability compared to . In the case of pseudo-static conditions, variations in the failure probability are linear in relation to , while in static conditions, these changes are exponential for an increase of 50% in the variation coefficient. Also, to reduce the coefficient of variation by 50%, the probability of static failure for different values of  is approximately zero.
./files/site1/files/132/8Extended_Abstracts.pdf

---

Aggregates are one of the high demand building materials in construction of structures and their characteristics have important effects on durability and permanence of projects. Abrasion resistance is one of the important features of aggregates that their utilization in concrete and asphalt are affected by texture and lithology of them. As rock consisted of harder minerals have higher abrasion resistance like igneous rocks, due to more siliceous minerals. More varieties in mineralogy compound usually lead to increase in aggregate abrasion. Aggregates that are contained of different minerals usually have less abrasion resistance. Porosity usually decreases the resistance abrasion. In addition to lithological properties, the environment where aggregates are deposited is important in determining resistance-related parameters of aggregates.
Rivers, alluvial fans, and taluses are the main environments where aggregates are deposited. Geological processes, such as weathering and particle movement may cause changes in natural aggregates, hence affecting their abrasion and impact resistance. Rock weathering can results in increasing porosity, producing minerals that are weaker in comparison to their original rock.
In the process of particles transport by stream water, weak parts of aggregates will be omitted. The present study is focused on the relationship between geology medium and the weight loss of aggregate in Los Angeles test. 
Methodology
Considering that lithology features in aggregates resistance against abrasion have an important role, in order to examine the effect of various geology environments in abrasion resistance of aggregates, the medium should be chosen having similar lithology. Therefore, the north of Damavand and the south of Daneh Khoshk anticline (north of Dire plain) were firstly chosen by using geology map, satellites images and field study. Damavand zone consists of trachyte and trachy-andesite volcanic rocks. These rocks cover the whole area around the Damavand peak. Also, Daneh Khoshk anticline is covered by thick Asmari formation. The selected environment are in the length of each other. Such that taluses feed alluvial fan and alluvial fans feed rivers. Samples were collected from different area of southern part of anticline. 10 river area, 12 alluvial fan and 6 taluses in the south-west area of Daneh Khoshk anticline (north of Dire plain) were chosen. Los Angeles test has been done according to standard A method ASTM D2216-10, 1990 on samples and the results were analyzed by analogous analyzer.
Results and discussion
Results show that porosity and micro-crack percentage increase, respectively in accumulated aggregate in river, alluvial fans and taluses areas. Also, porosity and micro-crack in various alluvial fans is different and is influenced by the area and length of main channel of alluvial fans’ catchment. The porosity decreases by the increase in the length of channel and area of alluvial fans’ catchment.
The percentages of aggregate weight loss in talus, alluvial fan and river areas decreases, respectively. Based on the obtained results, the lowest rates of weight loss belong to river environments (23.7 % in Daneh Khoshk and 42% in Damavand) whereas the highest rates of weight loss belong to taluses (49.3% in Daneh Khoshk and 48% in Damavand). The alluvial fans have an average state. Another noticeable point is the high weight loss in Los Angeles test in Damavand aggregate. Due to having harder mineral, igneous aggregate have more abrasion resistance, but this research illustrates that the weight loss resulting from Los Angeles test in these aggregates is high. This is because of virtues texture that weakness against the impact as well as their high porosity.
Conclusion
The result of this research indicates that the volume of aggregate weight loss in Los Angeles test is related to aggregate accumulation environment. The extent of aggregate abrasion resistance is lowest in talus medium and increases in alluvial fan and river environment, respectively. The difference in aggregate abrasion resistance in various areas result from geology process differences that applies to aggregates in various environment. The extent of caring particles in talus environment is very low and the type of movement is mass or sliding type in these media, micro-crack and weak parts remains within aggregates. The surface of micro crack is weak such that breaks easily in Los Angeles test due to the pressure results from the impact of aggregate, as well as the impact of steel ball on aggregate leading to aggregate breakages. Aggregates move more distances in alluvial fan and river. Aggregate strike together in riverbed and alluvial fan yielding to aggregates breakages from micro-cracks. As the movement distance increases, aggregates approach more to intact rock. During the particles move, the weathered and weak parts are damaged by aggregate abrasion to riverbeds and alluvial fan, and more resistant and harder aggregates remain. As the water current increases, the aggregates impact each other harder, more resistant micro-crack breakages and this change leads to decrease the weight loss in Los Angeles test.
./files/site1/files/132/1Extended_Abstracts.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
One-dimensional site response analysis is widely performed to account for local site effects during an earthquake. Most of these approaches assume that dynamic soil properties are frequency independent. Laboratory test results as well as in-situ testing show that shear modulus and damping ratio are dependent on the frequency of loading. Although the amplification factor at ground surface with respect to frequency dependent dynamic properties of mixed alluvium materials under different near-fault motions with various velocity period is recognized, it is not well characterized and quantified.
Material and analysis methods
In this study, the tests results of samples which obtained from the drilling borehole (BH14) form Pardis city in Iran, are used. The soil is classified as clayey of high plasticity/clayey sand (CH/SC) and almost uniform and similar in the whole log profile.
Shear modulus and damping ratios versus shear strain curves (ASTM D3999) of CH/SC natural materials at effective confining pressures of 1, 2 and 5 kg/cm² with frequency of 0.5, 2, 5, and 10 Hz were used in one dimension response analyses using EERA Code.
Generally the damping ratio versus shear strain of the studied materials under low loading frequency (i.e. 0.5 Hz) almost falls in the range identified in literature. However, at higher loading frequencies (5 and 10 Hz) the damping ratios completely fall above the known upper bound trend. It is observed that, in general, the G and D values increase as loading frequency increases. Moreover, at certain strain G/G_max ratio decreases as loading frequency is increased.
Different dynamics behaviour curves were used in analyses, in isotropic consolidation conditions. In order to assess the amplification, acceleration spectra, acceleration spectra ratio, coefficient of B, at ground surface under eight well-known near-fault ground motions, 1728 one dimensional analyses were carried out with EERA code. The analyses have been performed for three base acceleration levels, namely, 0.1 g, 0.35 g and 1 g, using the simple time history scaling method. Field and laboratory test results of shear wave velocity were used in the analyses.
In this study, several well-known near-fault motion records are utilized for ground response analyses. Near-fault earthquakes records were selected from the strong motion database of the Pacific Earthquake Engineering Research Center (PEER) and Iran Strong Motion Network (ISMN) for specific reasons of location of the near faults sites.
In current building codes, the upper 30 m soil deposits overlying the higher impedance earth crust are regarded as the most relevant and significant in characterizing the seismic behavior of a site. Therefore, it is useful to accomplish investigations for obtaining their amplification and spectral acceleration for 30 m and even thicker (e.g. 60 m, for usual deep excavation in Iran), in order to have economical and safe designs and constructions.
Results and discussion
Figure 1 presents a comparison of normalized spectral acceleration (B factor) versus period for 30 m and 60 m thick profiles and Vs testing for frequencies dependent and independent analyses under input base acceleration of 0.35g for longitudinal component of used earthquakes. B factor of Iranian Standard 2800 and UBC97 also has been presented in the figure. The spectral acceleration at short period for frequency dependent analysis is higher than that of the frequency independent analysis. The  increases in frequency dependent analysis and higher thick profile (i.e. 60 m).
Conclusion
Results show that the effect of loading frequency has a considerable influence on the acceleration response at the ground surface. For both 30 m and 60 m soil columns, the increase of the loading frequency, decreases the amplification factor especially in the short period zone of the spectra. Based on the acceleration response spectra of near field strong motions derived for soils types of I and IV in this study, the period corresponding to  in the design spectrum of Iranian Standard 2800 should increase to 0.5 and 1.4, respectively. Therefore, selection of the appropriate G and D curves measured at frequency similar to those of the anticipated cyclic loading (e.g. seismic) has a paramount importance../files/site1/files/134/1.pdf

---

Introduction
Rock mass deformation modulus is one of  the major parameters has to be considered in the design phase of arch dams. Due to filling and discharging of reservoir and corresponding loading and unloading on the dam abutments, irreversible deformation takes place within the rock mass and consequently, increases the potential of creating a separation between dam body and abutments. Therefore, the rock mass modulus must be more than an alowable value in order to prevent arch dam failure. Regarding small core samples and lack of joints and other similar discontinuities in samples, the determined modulus through performing laboratory tests is higher than those obtained through in-situ tests. The available technique to estimate the rock mass deformation modulus is divided into two classes as direct and indirect methods. In direct methods, the rock mass deformation modulus is measured via performing in-situ tests such as plate loading test while it is estimated through empirical equations using rock mass classification and laboratory test results in indirect methods. These equations are developed based on regression analysis between the rock mass modulus calculated via in-situ tests, the rock mass classification and laboratory test results. Although application of these equations is simple and cost-effective, the results are doubtful and cannot be used in the design phase of arch dam due to the heterogeneous nature of rock mass and rock type variability. The numbers of micro-cracks which are developed after gallery excavation using drilling and blasting technique are more close to the loading plate. Thus, calculated modulus in these points is lower than reality. The displacement in the points far from loading plate was near to zero while the transmitted load which is calculated applying ASTM D4394 standard is more than reality in small galleries. Consequently, the calculated modulus was extremely larger than real values and sometimes even more than intact value. The empirical equations are site dependent and they are just applicable in sites with similar geotechnical condition. It is obvious that in-situ tests, such as plate loading, are the appropriate method in order to determine the modulus of deformation, however, due to some simplification in the data processing such as semi-infinite boundary condition, the application of numerical simulation as a data processing tool is more appropriate. In this research, the Beheshtabad dam was introduced and the geology characteristics of dam site were investigated. Applying direct and indirect methods, the rock mass modulus of dam abutments is calculated.
Material and Methods
The dam site is placed approximately at a distance of 2.7 km from the intersection of Koohrange and Beheshtabad river. In accordance with geological studies, the rocks in the site could be categorized in four units combined of Dolomite, Dolomitic Limestone, Limestone, Marl and Marly Limestone. Applying empirical equation the rock mass modulus of dam abutments is evaluated based on the laboratory test results and rock mass engineering classification systems. In addition, ASTM D4394 is applied to investigate the results of ten plate loading tests which are executed in the right and left abutments. To interpret the plate loading test results in the right abutment, a three-dimensional Fast Lagrange Analysis of Continuum (FLAC3D) model is developed.
Result and Discussion
To process the numerical simulation results, back analysis as a data processing tool is used. In this approach, the input parameters of numerical model will be changed in the way that the measured quantities by extensometers at the monitoring points are almost equal with the computed ones via numerical model at the corresponding points. Based on the sensitivity analysis carried out on the Mohr-Coulomb failure criterion parameters, the friction coefficient and cohesion variation do not affect the displacements calculated via numerical simulation as the more portion of gallery displacements are elastic. The error function is minimum when the rock mass modulus is 12 GPa and the horizontal to vertical stress ratio (K0) is equal to 0.5. The evaluated rock mass modulus based on the numerical simulation is two times lower than corresponding one evaluated applying empirical equation as a result of empirical equation uncertainty. Consideration of stress decrement under loading plate shows lower level of stress decrement under loading plate in ASTM D4394 compared to numerical simulation. This is why, the rock mass modulus, calculated based on ASTM D4394, increases dramatically by getting distance from the loading plate. 
Conclusion
The empirical methods estimating the modulus of deformation based on rock mass classification systems tend to evaluate large value of modulus especially for the weak massive rocks.
As a result of galleries dimensions and semi-infinite boundary condition assumed in ASTM D4394, the calculated rock mass modulus increases dramatically by getting distance from loading plate. Therefore, the numerical simulation was applied to process the plate loading test results. A new normalized error function was developed based on measured displacements and the rock mass modulus in the right abutment was determined 12 GPa which is very lower than the calculated value using ASTM D 4394. Also, as a result of numerical simulation, the rock mass is uniform. The stress increment perpendicular to the loading plate was calculated applying numerical simulation which is 0-90 percent lower than those suggested by ASTM D 4394. 
 

---

Introduction
Soil dynamic properties are used to evaluate the dynamic response of soils at different strain levels in geotechnical engineering. The shear modulus (G) and damping ratio (D) are among the most important dynamic properties of soils. In general, the factors affecting the dynamic behavior of soils are divided into two categories: first; soil type and characteristics such as water content, void ratio and soil plasticity and second; parameters of loads applied on the soil such as the number of loading cycles, loading frequency and loading waveform .Therefore, it is widely believed that the dynamic response of soils partially depends on the characteristics of the load. In this paper, 1-g shaking table tests were employed to investigate the effect of loading waveform and frequency content on dynamic properties of dry sands. The response obtained from soil samples during loading with different frequencies, input accelerations and waveforms were used to generate hysteresis loops of tested samples at different strain amplitudes. Then, hysteresis loops were used to determine the damping ratio and shear modulus at different strain levels. Finally, the effects of loading frequency and waveform on the changes of each parameter (G and D) were investigated.
Material and methods
A hydraulic shaking table with a single degree of freedom, designed and constructed at the Crisis Management Center of Urmia University, was used to conduct the experiments. Firoozkuh No. 161 sand was used in all the experiments. The Firoozkuh sand gradation curve is similar to that of Toyoura sand. In this study, accelerometers were used to measure the acceleration of the input to the sample as well as to record the acceleration caused by the input excitation at different depths of the soil sample. The displacement transducers (LVDT sensors) were also used to measure linear displacement. Each soil sample was constructed using dry Firoozkuh sand poured uniformly into the container from four equal heights of 150 mm to reach a total height of 600 mm. During the compaction process, the accelerometers A₁, A₂, and A₃ were placed at a depth of 150, 300 and 450 mm with respect to the bottom of container. Also, one accelerometer, A₀, was attached rigidly to the container base to measure base acceleration. A displacement transducer (L₁) was placed on the soil surface at a height of 600 mm from the floor of the container to measure the vertical displacement of the surface of the soil. In this study, 42 shaking table tests were performed to study the effect of loading frequency and waveform on dynamic properties of dry sand. The test samples were subjected to rectangular, sinusoidal and triangular loading at frequencies of 0.5 to 9 Hz and at input acceleration of 0.1g and 0.3 g.
Results and discussion
Given the importance of G-γ and D-γ curves in dynamic analyses, the changes in shear modulus with shear strain has been studied. The results show that the shear modulus increases as the frequency increases in all cases, and this increase is observed at lower frequencies and increases with increasing frequency. On the other hand, the shear modulus decreases with increasing shear strain. At a constant testing frequency, soil samples subjected to the rectangular waveform exhibited the largest shear modulus while the samples subjected to the triangular waveform showed the least shear modulus. The shear modulus of the samples under the sinusoidal waveform is barely more than the shear modulus of samples under triangular waveform. Moreover, by increasing the shear strain, the shear modulus values ​​of samples with different waveforms have become closer and the waveform effect is reduced. As for the effect of input acceleration on the shear modulus , increasing the input acceleration increases the shear strain and consequently, decreases the shear modulus in all states (the values ​​of shear modulus in various frequencies and the waveforms under the input acceleration of 0.1 g are larger than the shear modulus values ​​under the input acceleration of 0.3g). In the case of the damping ratio, the results show that, in all cases, damping ratio increases with shear strain. At low strain levels, the damping ratio values at various frequencies and waveforms are low and yet very close. At higher strain levels, the increase in frequency increases the damping ratio. This increase is more significant at higher frequencies. Also, the effect of waveform on the damping ratio is more apparent at larger shear strains, and at such shear strain levels, soil samples under rectangular loading exhibit the largest damping ratio. The damping ratio associated with the sinusoidal and triangular loading are also close to each other and it is a slightly larger for sinusoidal loading. On the other hand, the damping ratio increases with input acceleration. In addition, the effect of increased input acceleration on the increase in the damping ratio is more obvious at higher frequencies mainly due to the increase in shear strain.
Conclusion
In the present study, the effects of loading frequency and waveform on the dynamic properties of dry sand were investigated using shaking table tests. The following conclusions were drawn:

• The shear modulus increases with frequency. The trend is more obvious at larger frequencies. The effect of loading frequency on the damping ratio of the soil at low levels of strain is negligible, and at relatively large strain levels, damping ratio increases with loading frequency.
• Soil samples exhibit the highest shear modulus and damping ratio under rectangular loading. Therefore, in all the tested frequencies and input accelerations, the values of G and D for the rectangular waveforms are greater than those of the sinusoidal and triangular waveforms. The shear modulus and damping ratio for the sinusoidal waveforms are marginally greater than those of triangular waveforms, yet one can consider them practically similar.
• In all cases, the shear strain increased by increasing the amplitude of the input acceleration, and as a result, the shear modulus decreased and the damping ratio increased../files/site1/files/142/2.pdf

---

Earth dams are geotechnical structures constructed on various shapes of a valley. The Vanyar Dam is a rock-fill dam located on a narrow valley. Concerning the geometry of the canyon, three-dimensional modeling was utilized to analyze this dam. According to the numerical analysis, the maximum settlement is 88.14 cm, which corresponds to 48 m above the bedrock in cross-section C, that is, a little less than 1% of the dam height. Besides, the total vertical stresses recorded by the pressure cells are about 28% less than those obtained from the numerical analysis. It is assumed that the difference is caused by local arching due to lower compaction and consequently a low stiffness area around the pressure cells. In terms of pore water pressure, there is good agreement between the pore water pressure obtained from the numerical analysis and the piezometers, such that the results are restricted to less than 1%. In general, the difference between the numerical analysis results and those recorded by the instruments is acceptable. Furthermore, the dam shows a suitable level of performance at the end of construction.

---

---

---

---

---

---