Limiting equilibrium method for slope/drilled shaft system

The use of drilled shafts to stabilize an unstable slope has been a widely accepted practice. There are two basic design and analysis issues involved: one is to determine the global factor of safety of the drilled shafts stabilized slope and the other one is to determine the design earth thrust on the drilled shafts for structural design of the shafts. In this paper, a limiting equilibrium method of slices based solution for calculating global factor of safety (FS) of a slope with the presence of a row of drilled shafts is developed. The arching mechanisms due to the presence of the drilled shafts on slope were taken into account by a load transfer factor. The method for calculating the net force applied to the drilled shaft from the soil mass was also developed. The interrelationships among the drilled shaft location on the slope, the load transfer factor, and the global FS of the slope/shaft system were derived utilizing the developed numerical closed‐form solution. An illustrative example is presented to elucidate the use of the solution in optimizing the location of the drilled shafts on slope to achieve the desired global factor of safety of the slope/shaft system. Copyright © 2009 John Wiley & Sons, Ltd.     

Mineral chemical composition and geodynamic significance of peridotites from Nain ophiolite,central Iran

The peridotites from north of the town of Nain in central Iran consist of clinopyroxene-bearing harzburgite and lherzolite with small lenses of dunite and chromitite pods. The lherzolite contains aluminous spinel with a Cr number (Cr^# = Cr/[Cr + Al]) of 0.17. The Cr number of spinels in harzburgite and chromitite is 0.38–0.42 and 0.62, respectively. This shows that the lherzolite and harzburgite resulted from <18% of partial melting of the source materials. The estimated temperature is 1100 ± 200 °C for peridotites, the estimated pressure is <15 ± 2.3 kbar for harzburgites and >16 ± 2.3 kbar for lherzolites and estimated fo₂ is 10^?1±0.5 for peridotites. Discriminant geochemical diagrams based on mineral chemistry of harzburgites indicate a supra-subduction zone (SSZ) to mid-oceanic ridge (MOR) setting for these rocks. On the basis of their Cr^#, the harzburgite and lherzolite spinels are analogous to those from abyssal peridotites and oceanic ophiolites, whereas the chromites in the chromitite (on the basis of Cr^# and boninitic nature of parental melts) resemble those from SSZ ophiolitic sequences. Therefore, the Nain ophiolite complex most likely originated in an oceanic crust related to supra-subduction zone, i.e. back arc basin. Field observations and mineral chemistry of the Nain peridotites, indicating the suture between the central Iran micro-continent (CIM) block and the Sanandaj–Sirjan zone, show that these peridotites mark the site of the Nain–Baft seaway, which opened with a slow rate of ocean-floor spreading behind the Mesozoic arc of the Sanandaj–Sirjan zone as a result of change of Neo Tethyan subduction régime during middle Cretaceous.     

Development in modeling cyclic loading of sands based on kinematic hardening

In this paper, there is presented an elastoplastic constitutive model to predict sandy soils behavior under monotonic and cyclic loadings. This model is based on an existing model (Cambou‐Jafari‐Sidoroff) that takes into account deviatoric and isotropic mechanisms of plasticity. The flow rule used in the deviatoric mechanism is non‐associated and a mixed hardening law controls the evolution of the yield surface. In this research the critical state surface and history surface, which separates the virgin and cyclic states in stress space, are defined. Kinematic hardening modulus and stress–dilatancy law for monotonic and cyclic loadings are effectively modified. With taking hardening modulus as a function of deviatoric and volumetric plastic strain and with defining the history surface and stress reversal, the model has the ability to predict the sandy soils' behavior. All of the model parameters have clear physical meanings and can be determined from usual laboratory tests. In order to validate the model, the results of homogeneous tests on Hostun and Toyoura sands are used. The results of validation show a good capability of the proposed model. Copyright © 2009 John Wiley & Sons, Ltd.     

Spatial distribution and temporal variation of reference evapotranspiration in arid and semi‐arid regions of Iran

Analysis of spatial and temporal variations of reference evapotranspiration (ET_o) is important in arid and semi‐arid regions where water resources are limited. The main aim of this study was to analyse the spatial distribution and the annual, seasonal and monthly trends of the Penman–Monteith ET_o for 21 stations in the arid and semi‐arid regions of Iran. Three statistical tests the Mann‐Kendall, Sen's slope estimator and linear regression were used for the analysis. The analysis revealed that ET_o increased from January to July and deceased from July to December at almost all stations. Additionally, higher annual ET_o values were found in the southeast of the study region and lower values in the northwest of the region. Although the results showed both positive and negative trends in annual ET_o series, ET_o generally increased, significantly so in six (～30%) of the stations. Analysis of the impacts of meteorological variables on the temporal trends of ET_o indicated that the increasing trend of ET_o was most likely due to a significant increase in minimum air temperature, while decreasing trend of ET_o was mainly caused by a significant decrease in wind speed. At the sites where increasing ET_o trends were statistically significant, the rate of increase varied from (+)8·36 mm/year at Mashhad station to (+)31·68 mm/year at Iranshahr station. On average, an increasing trend of (+)4·42 mm/year was obtained for the whole study area during the last four decades. Seasonal and monthly ET_o have also tended to increase at the majority of the stations. The greatest numbers of significant trends were observed in winter on the seasonal time‐scale and in September on the monthly time‐scale. Copyright © 2011 John Wiley & Sons, Ltd.     

List of forthcoming papers

List of forthcoming papers     

Experimental and 3D Numerical Simulation of Reinforced Shear Joints

The load transfer capacity and failure mechanism of a fully grouted bolt installed across a joint in shear is investigated, both experimentally and numerically, in five types of bolt. The double-shearing testing of bolts were studied in concrete blocks of 20, 40 and 100 MPa strengths, subjected to different pretension loads of 0, 5, 10, 20, 50 and 80 KN, respectively. The parameters examined include: shear resistance, shear displacement, induced strains and stresses during the bolt-bending process, and its ultimate failure across the sheared joint planes. The conclusions drawn from the study were that the strength of the concrete, bolt profile configuration and bolt pretension load played a significant influence on the shear resistance, shear displacement and failure mechanism of the reinforced medium. Experimental and numerical simulations showed that the failure occurs as a result of the induced axial and shear stresses acting between the hinge point distances in the vicinity of the shear joint plane.     

Prediction of Land Subsidence Under Cyclic Pumping Based on Laboratory and Numerical Simulations

In this study subsidence due to groundwater withdrawal was investigated. Kerman Province in Iran is struggling with land subsidence problem due to extensive groundwater withdrawal mainly for farming. The rate and type of groundwater withdrawal has very important impact on settlement rate. In this research, effective parameters on land subsidence caused by groundwater withdrawal were determined by laboratory tests. Sampling had done up to depth of 300 m mainly with remolded specimens from Shams-abad, Nouq plain in Kerman province. Similar to the field preconsolidation pressure was applied on specimens in the laboratory. Rate of applied stress on prepared specimens was similar to effect of oscillation of groundwater level. In order to model the actual soil behavior in the laboratory, one-dimensional consolidation device (odometer) was adopted for testing. In these tests, the effect of loading caused by seasonal oscillation of groundwater table is considered by means of cyclic loading in the testing which has great effect on rate of settlements. The results of tests show that when the water table level periodically increases and decreases the amount of settlement decrease, comparing with the case when the groundwater table drop to a constant level. In order to predict the further effects of groundwater level oscillation and actual field condition on land subsidence, a finite element model based on Biots’ three-dimensional consolidation theory was developed. After calibration of finite element model with laboratory tests, this model was used for prediction the effect of groundwater level oscillation on actual field conditions.     

Mixed precision parallel preconditioner and linear solver for compositional reservoir simulation

Reservoir simulation role in value creation and strategic management decisions cannot be over emphasized. Simulation of complex challenging reservoirs with millions of grid blocks especially in compositional mode is very time-consuming even with fast modern computers. On the other hand, high price of cluster supercomputers prevents them for being commonly used for fast simulation of such reservoirs. In recent years, the development of many-core processors like cell processors, DSPs, and graphical processing units (GPUs) has provided a very cost-effective hardware platform for fast computational operations. However, programming for such processors is much more difficult than conventional CPUs, and new parallel algorithm design and special parallel implementation methods are needed. Using the computational power of CPUs, GPUs, and/or any other processing unit, Open Computing Language (OpenCL) provides a framework for programming for heterogeneous platforms. In this paper, OpenCL is used to employ the computational power of a GPU to build a preconditioner and solve the linear system arising from compositional formulation of multiphase flow in porous media. The proposed parallel preconditioner is proved to be quite effective, even in heterogeneous porous media. Using data-parallel modules on GPU, the preconditioner/solver runtime reduced at least 1 order of magnitude compared to their serial implementation on CPU.     

Optimal water and waste load allocation in reservoir–river systems: a case study

In this paper, a new methodology is developed for optimization of water and waste load allocation in reservoir–river systems considering the existing uncertainties in reservoir inflow, waste loads and water demands. A stochastic dynamic programming (SDP) model is used to optimize reservoir operation considering the inflow uncertainty, and another model called PSO-SA is developed and linked with the SDP model for optimizing water and waste load allocation in downstream river. In the PSO-SA model, a particle swarm optimization technique with a dynamic penalty function for handling the constraints is used to optimize water and waste load allocation policies. Also, a simulated annealing technique is utilized for determining the upper and lower bounds of constraints and objective function considering the existing uncertainties. As the proposed water and waste load allocation model has a considerable run-time, some powerful soft computing techniques, namely, Regression tree Induction (named M5P), fuzzy K-nearest neighbor, Bayesian network, support vector regression and an adaptive neuro-fuzzy inference system, are trained and validated using the results of the proposed methodology to develop real-time water and waste load allocation rules. To examine the efficiency and applicability of the methodology, it is applied to the Dez reservoir–river system in the south-western part of Iran.     

Interactive effect of nitrogen fertilizer and hydrocarbon pollution on soil biological indicators

The investigation of the impact of different forms of nitrogen fertilizer (NO₃-N and NH₄-N) on microbial parameters, enzyme activities and phytotoxicity in a petroleum-contaminated soil was evaluated by an incubation study. The tested enzymes, microbial activity and seed germination index showed different patterns in response to both petroleum and nitrogen fertilizer addition and time of incubation. The results apparently showed that the contamination of soil with petroleum has a negative effect on soil ecosystem. Nitrogen fertilizer could improve inhibition of petroleum hydrocarbons in soil. Nevertheless, nitrogen fertilizer had no significant effect on urease activity in the petroleum-contaminated soil. As compared to NO₃-N, the addition of NH₄-N to the soil resulted in a greater impact on soil performance as attested by the recovery of the soil germination capability and higher values of the respiration. The application of nitrogen fertilizer may be suggested as a good strategy for restoring soils in regions affected by the same problem.