---

The reliable operation of spillways, in emergency as well as normal conditions, is one of the vital components in dam safety. Free or uncontrolled overflow spillways are the most reliable choice however. They usually impose higher construction cost and /or results in wasting a considerable amount of water or live capacity of the reservoirs. Employing fuse gates might be a way of reconciling dam safety with maximized storage capacity. The operation of the system can be controlled to within a few centimeters, and the entire installation is not lost for floods less than the maximum design flood. The installation offers more or less the same level of safety as ungated spillways, but avoids their inherent storage capacity loss. Optimum design of fuse gates in particular installation calls for a mathematical model. The model developed in this work includes structural, hydraulics and operational constraints while maximizing the expected cost over the useful life of the project. Accounting for the lost benefit (i.e., water lost as a result of gate tilting) has an influenced effect on the optimum design. To test the performance of the model, data from Zarineh Rud dam in Iran has been used and its result is compared with a direct search technique. The model is capable of helping the design engineer to select the best alternative considering different types of constraints.

---

In this paper the analysis of the pipe networks is formulated as a nonlinear unconstrained optimization problem and solved by a general purpose optimization tool. The formulation is based on the minimization of the total potential energy of the network with respect to the nodal heads. An analogy with the analysis of the skeletal structures is used to derive tire formulation. The proposed formulation owes its significance for use in pipe network optimization algorithms. The ability and versatility of the method to simulate different pipe networks are numerically tested and the accuracy of the results is compared with direct network algorithms.

---

This paper presents the long-term deformations of reinforced high-strength concrete columns subjected to constant sustained axial forces. The objective of the study was to investigate the effects of binder systems containing different levels of silica fume on time-dependent behaviour of high-strength concrete columns. The experimental part of the work focused on concrete mixes having a fixed water/binder ratio of 0.35 and a constant total binder content of 500 kg/m3. The percentages of silica fume that replaced cement in this research were: 0%, 6%, 8%, 10% and 15%. The mechanical properties evaluated in the laboratory were: compressive strength secant modulus of elasticity strain due to creep and shrinkage. The theoretical part of the work is about stress redistribution between concrete and steel reinforcement as a result of time-dependent behaviour of concrete. The technique used for including creep in the analysis of reinforced concrete columns was age-adjusted effective modulus method. The results of this research indicate that as the proportion of silica fume increased, the short-term mechanical properties of concrete such as 28-day compressive strength and secant modulus improved. Also the percentages of silica fume replacement did not have a significant influence on total shrinkage however, the autogenous shrinkage of concrete increased as the amount of silica fume increased. Moreover, the basic creep of concrete decreased at higher silica fume replacement levels. Drying creep (total creep - basic creep) was negligible in this investigation. The results of the theoretical part of this researchindicate that as the proportion of silica fume increased, the gradual transfer of load from the concrete to the reinforcement decreased and also the effect of steel bars in lowering the concrete deformation reduced. Moreover, the total strain of concrete columns decreased at higher silicafume replacement levels.

---

A study of bearing capacity and compressibility characteristics of cohesive soil, reinforced by geogrid and supporting square footing loads has been conducted. The lack of adequate frictional resistance between clay and reinforcing elements was compensated by using a thin sand layer (lens) encapsulating the geogrid sheet. In this way, tensile forces induced in the geogrid were transferred to the bulk clay medium through the sand particles and soil reinforcement was improved Experiments were conduced on two sets of specimens, one set of 1 x 1 x 1 m dimension and the footing size of 19 x 19 cm (series A), and the other set of 0.15 x 0.15 x 0.15 m dimension and the footing size of 3.7 x 3.7 cm (series B). The loading systems for the above specimens were stress controlled and strain controlled respectively. All specimens were saturated and presumably loaded under an undrained condition. The results qualitatively confirmed the effectiveness of the sand lens in improving the bearing capacity and settlement characteristics of the model footing. In series A, the maximum increase in the bearing capacity due to the presence of the sand lens was 17% whereas in series B, the amount of increase was 24%. The percentage reductions in the settlement for these results were 30% and 46% respectively.

---

In this paper a structural damage detection algorithm using static test data is presented. Damage is considered as a reduction in the structural stiffness (Axial and/or Flexural) parameters. Change in the static displacement of a structure is characterized as a set of non-linear undetermined simultaneous equations that relates the changes in static response of the structure to the location and severity of damage. An optimality criterion is introduced to solve these equations by minimizing the difference between the load vector of damaged and undamaged structures. The overall formulation leads to a non-linear optimization problem with non-linear equality and linear inequality constraints. A method based on stored strain energy in elements is presented to select the loading location. Measurement locations are selected based on Fisher Information Matrix. Numerical and experimental results of a 2D frame represent good ability of this method in detecting damages in a given structure with presence of noise in measurements.

---

Having observed the costly failures of different cutoff walls, that had been constructed according to the mix design specified by reputable consultants in Iran, a research programme was conducted to study the effects of constituent materials on the properties of plastic concrete. The main properties, such as compressive strength, biaxial and triaxial strains, permeability, and modulus of elasticity have been investigated using different mixes, obtained from prototype production line plant, situated on site, because it was realized that the site production line and the systems employed have major effects on the properties of plastic concrete. Statistical analysis of the results, revealed the coefficients of influence of main constituent materials of plastic concrete namely cement, bentonite, aggregate and water on its compressive strength and modulus of elasticity. Having realized the cancelling effects of bentonite and aggregates on the measured properties, some equations relating the quantities of cement and water to the compressive strength and modulus of elasticity are introduced. Effects of clay and hydrated lime powder, as fillers were also investigated leading to the proposal of limits for their safe and economic use. Since most of the cutoff walls are buried structures, failure strains under both uniaxial and triaxial tests, with values of cohesion and internal friction, are also presented in this paper.

---

Existence of discontinuities causes higher deformability and lower strength in rock masses. Thus joints can change the rock mass behaviour due to the applied loads. For this reason properties and orientation of the joint sets have a great effect on the stability of rock slopes. In this paper, after introducing some numerical methods for evaluating the factor of safety for the stability of slopes, stability of jointed rock slopes in the plane strain condition is investigated with the strength reduction technique this method is modified and applied in the multilaminate framework. First of all, stability of one homogeneous rock slope is investigated and compared with the limit equilibrium method. Then stability of a layered rock slope is analyzed with some modifications in the strength reduction technique. Effects of orientation, tensile strength and dilation of layered joint sets on the factor of safety and location of the sliding block are explained.

---

In this paper a general procedure for automated minimum weight design of twodimensional steel frames under seismic loading is proposed. The proposal comprises two parts: a) Formulation of automated design of frames under seismic loading and b) introduction of an optimization engine and the improvement made on it for the solution of optimal design. Seismic loading, that depends on dynamic characteristics of structure, is determined using "Equivalent static loading" scheme. The design automation is sought via formulation of the design problem in the form of a standard optimization problem in which the design requirements is treated as optimization constraints. The Optimality Criteria (OC) method has been modified/improved and used for solution of the optimization problem. The improvement in (OC) algorithm relates to simultaneous identification of active set of constraints and calculation of corresponding Lagrange multipliers. The modification has resulted in rapid convergence of the algorithm, which is promising for highly nonlinear optimal design problems. Two examples have been provided to show the procedure of automated design and optimization of seismic-resistant frames and the performance and capability of the proposed algorithm.

---

An important concern in rock mechanics is non-homogeneity as joints or fault. Adopting the joints as fractures, fractures are well known for their effects on the mechanical and transport properties of rock. It has been postulated that through fractured/jointed rock, mainly, the polygons turned to the shear vector (ti) are involved in the mobilization of shear resistance. Consequently, in order to locate the contact areas implicated into the shear-test it was firstly necessary to fix the shear direction. Moreover, since laboratory observations clearly show that only the steepest polygon surfaces touch the other sample, the identification of the potential sliding areas only requires the determination of the polygons which are faced to the shear direction and which, among them, are steep enough to be involved. The methodology to be discussed here is modeling of slip on the local and global levels due to the distribution of deformation procedure of the rock joint. Upon the presented methodology, more attention has been given to slip initiation and propagation through rock joint. In particular, softening in non-linear behaviour of joint in going from the peak to residual strengths imparts a behaviour often associated with progressive failure. A multi-plane based model is developed and used to compute plastic strain distribution and failure mechanism of rock joints. Validity of the presented model was examined by comparing numerical and test results showing the behavior of both homogeneous and jointed rock samples under general stress conditions.

---

In this paper, a higher order continuum model is presented based on the Cosserat continuum theory in 3D numerical simulation of shear band localization. As the classical continuum models suffer from the pathological mesh-dependence in strain softening models, the governing equations are regularized by adding the rotational degrees-of-freedom to conventional degrees-of-freedom. The fundamental relations in three-dimensional Cosserat continuum are presented and the internal length parameters are introduced in the elasto-plastic constitutive matrix to control the shear bandwidth. Finally, the efficiency of proposed model and computational algorithm is demonstrated by a 3D strip in tensile. A comparison is performed between the classical and Cosserat theories and the effect of internal length parameter is demonstrated. Clearly, a finite shear bandwidth is achieved and the load-displacement curves are uniformly converged upon different mesh sizes.

---

Seismic behavior of a rockfill dam with asphalt-concrete core has been studied utilizing numerical models with material parameters determined by laboratory tests. The case study selected for these analyses, is the Meyjaran asphalt core dam, recently constructed in Northern Iran, with 60 m height and 180 m crest length. The numerical analyses have been performed using a nonlinear three dimensional finite difference software and various hazard levels of earthquakes. This study shows that due to the elasto-plastic characteristics of the asphalt concrete, rockfill dams with asphalt concrete core behave satisfactorily during earthquake loading. The induced shear strains in the asphalt core, for the case presented in this research, are less than 1% during an earthquake with amax=0.25g and the asphalt core remains watertight. Due to large shear deformations caused by a more severe earthquake with amax=0.60g, some cracking may occur towards the top of the core (down to 5-6 m), and the core permeability may increase in the top part, but the dam is safe.

---

This paper presents the results of a series of monotonic and post-cyclic triaxial tests carried out on a clay specimen and three types of clay-sand mixed specimens. The focus of the paper is on the post-cyclic mechanical behavior of the mixed specimens, as compared to their monotonic behavior. Analyses of the tests results show that cyclic loading degrade undrained shear strength and deformation modulus of the specimens during the post-cyclic monotonic loading. The degradation depends on the sand content, the cyclic strain level and to some degrees to the consolidation pressure.

---

Analysis and prediction of structural response to static or dynamic loading requires prediction of concrete response tovariable load histories. The constitutive equations for the mechanical behavior of concrete capable of seeing damage effects or crack growth procedure under loading/unloading/reloading was developed upon micro-plane framework. The proposed damage formulation has been built on the basis of five fundamental types of stress/strain combinations, which essentially may occur on any of micro-planes. Model verification under different loading/unloading/reloading stress/strain paths has been examined. The proposed model is capable of presenting pre-failure history of stress/strain progress on different predefined sampling planes through material. Many of mechanical behavior aspects happen during plasticity such as induced anisotropy, rotation of principal stress/strain axes, localization of stress/strain, and even failure mechanism are predicted upon a simple rational way and can be presented.

---

Explosion has always been regarded as one of the most complicated engineering problems. As a result, many engineers have preferred rather simplified empirical approaches in comparison to extremely complex deterministic analyses. In this paper, however, a numerical simulation based on the combined finite/discrete element methodology is presented for analyzing the dynamic behavior of fracturing rock masses in blasting. A finite element discretization of discrete elements allows for complex shapes of fully deformable discrete elements with geometric and material nonlinearities to be considered. Only a Rankine strain softening plasticity model is employed, which is suitable for rock and other brittle materials. Creation of new lines/edges/bodies from fracturing and fragmentation of original objects is systematically considered in the proposed gas-solid interaction flow model. An equation of state is adopted to inexpensively calculate the pressure of the detonation gas in closed form. The model employed for the flow of detonation gas has resulted in a logical algorithmic procedure for the evaluation of spatial distribution of the pressure of detonation gas, work done by the expanding gas and the total mass of the detonation gas as functions of time indicating the ability of model to respond to changes in both the mass of explosive charge and the size of the solid block undergoing fracture. Rock blasting and demolition problems are amongst the engineering applications that are expected to benefit directly from the present development. The results of this study may also be used to provide some numerical based reliable solutions for the complex analysis of structures subjected to explosive loadings.

---

A framework for development of constitutive models including damage progress, based on semi-micromechanical aspects of plasticity is proposed for concrete. The model uses sub-loading surface with multilaminate framework to provide kinematics and isotropic hardening/softening in the ascending/descending branches of loading and can be able to keep stress/stain paths histories for each plane separately. State of stresses on planes is divided to four basic stress patterns i.e. pure compression, increasing compression-and shear, decreasing compression-shear and tension-shear and used in derivation of plasticity equations. Under this kind of categorized form the model is capable of predicting behavior of concrete under any stress/strain path such as uniaxial, biaxial and triaxial in the monotonic and cyclic loading, Also this model is capable of predicting the effects of principal stress/strain axes rotations and consequent plastic flow and has the potential to simulate the behavior of material with anisotropy, fabric pattern, slip/weak planes and crack opening/closing. The material parameters of model are calibrated by optimum fitting of the basic test data available in the literature. The model results under both monotonic and cyclic loading have been compared with experimental results to show capability of model.

---

Local site conditions have a strong effect on ground response during earthquakes. Two important soil parameters that control the amplification effects of seismic motions by a soil column are the soil hysteretic damping ratio and shear wave velocity. This paper presents the results of in situ damping ratio measurements performed using continuous surface wave attenuation data at a site in Semnan University campus and analysis used to obtain the near surface soils damping ratio profile. Once the frequency dependent attenuation coefficients are determined, the shear damping ratio profile is calculated using an algorithm based on constrained inversion analysis. A computer code is developed to calculate the shear damping ratio in each soil layer. Comparisons of the in situ shear damping ratio profile determined from continuous surface wave with cross hole independent test measurements are also presented. Values of shear damping ratio, obtained using continuous surface wave measurements, were less than the measured using cross hole tests, possibly because of the higher frequencies used in cross hole tests.

---

The safety of buried pipes under repeated load has been a challenging task in geotechnical engineering. In this paper artificial neural network and regression model for predicting the vertical deformation of high-density polyethylene (HDPE), small diameter flexible pipes buried in reinforced trenches, which were subjected to repeated loadings to simulate the heavy vehicle loads, are proposed. The experimental data from tests show that the vertical diametric strain (VDS) of pipe embedded in reinforced sand depends on relative density of sand, number of reinforced layers and height of embedment depth of pipe significantly. Therefore in this investigation, the value of VDS is related to above pointed parameters. A database of 72 experiments from laboratory tests were utilized to train, validate and test the developed neural network and regression model. The results show that the predicted of the vertical diametric strain (VDS) using the trained neural network and regression model are in good agreement with the experimental results but the predictions obtained from the neural network are better than regression model as the maximum percentage of error for training data is less than 1.56% and 27.4%, for neural network and regression model, respectively. Also the additional set of 24 data was used for validation of the model as 90% of predicted results have less than 7% and 21.5% error for neural network and regression model, respectively. A parametric study has been conducted using the trained neural network to study the important parameters on the vertical diametric strain.

---

A new mathematical model for identifying pollution sources in aquifers is presented. The model utilizes Lagrange Constrained Optimization Method (LCOM) and is capable to inversely solve unsteady fluid flow in saturated, heterogeneous, anisotropic confined and/or unconfined aquifers. Throughout the presented model, complete advection-dispersion equation, including the adsorption as well as retardation of contaminant, is considered. The well-known finite element method is used to discretize and solve the governing equations. The model verification is implemented using a hypothetical example. Also, the applicability of the developed code is illustrated by the real field problem of Ramhormoz aquifer in southwestern Iran.

---

Plastic shrinkage is one of the most important parameter which must be considered in hot weather concreting. If plastic shrinkage is not prevented, cracking will be significant, especialy if silica fume is used in the mix. In this paper, the effect of silica fume in bleeding and evaporation was investigated in laboratory. The results showed that in restrained shrinkage, beside relative humidity, temperature and wind velocity, sun rediation also is very important factor in evaporation rate. It is found that under solar radition condition, the evaporation was much larger than the estimated value in ACI 305 Nomogram. The rate of evaporaion under solar radiation was about two folds of evaporation rate under shade condition. The results showed that in terms of crack initiation time, crack width and total cracking area, concrete containing silica fume is more severe than concrete with no silica fume. Reduction of water cement ratio in concrete with silica fume makes the concrete more sensitive in cracking. The results of this project also showed that the severity of the cracking is not related only to rate of bleeding but all environmental factors including like sun radiation or shading and also mix compositions have important roles.

---

Conventional investigations on the behavior of reinforced and unreinforced soils are often investigated at the failure point. In this paper, a new concept of comparison of the behavior of reinforced and unreinforced soil by estimating the strength and strength ratio (deviatoric stress of reinforced sample to unreinforced sample) at various strain levels is proposed. A comprehensive set of laboratory triaxial compression tests was carried out on wet (natural water content) non-plastic beach silty sand with and without geotextile. The layer configurations used are one, two, three and four horizontal reinforcing layers in a triaxial test sample. The influences of the number of geotextile layers and confining pressure at 3%, 6%, 9%, 12% and 15% of the imposed strain levels on sample were studied and described. The results show that the trend and magnitude of strength ratio is different for various strain level. It implies that using failure strength from peak point or strength corresponding to the axial-strain approximately 15% to evaluate the enhancement of strength or strength ratio due to reinforcement may cause hazard and uncertainty in practical design. Hence, it is necessary to consider the strength of reinforced sample compared with unreinforced sample at the imposed strain level. Only one type of soil and one type of geotextile were used in all tests.