Three‐dimensional nonlinear finite element analysis of pile groups in saturated porous media using a new transmitting boundary

The dynamic behaviour of pile groups subjected to an earthquake base shaking is analysed. An analysis is formulated in the time domain and the effects of material nonlinearity of soil, pile–soil–pile kinematic interaction and the superstructure–foundation inertial interaction on seismic response are investigated. Prediction of response of pile group–soil system during a large earthquake requires consideration of various aspects such as the nonlinear and elasto‐plastic behaviour of soil, pore water pressure generation in soil, radiation of energy away from the pile, etc. A fully explicit dynamic finite element scheme is developed for saturated porous media, based on the extension of the original formulation by Biot having solid displacement (u) and relative fluid displacement (w) as primary variables (u–w formulation). All linear relative fluid acceleration terms are included in this formulation. A new three‐dimensional transmitting boundary that was developed in cartesian co‐ordinate system for dynamic response analysis of fluid‐saturated porous media is implemented to avoid wave reflections towards the structure. In contrast to traditional methods, this boundary is able to absorb surface waves as well as body waves. The pile–soil interaction problem is analysed and it is shown that the results from the fully coupled procedure, using the advanced transmitting boundary, compare reasonably well with centrifuge data. Copyright © 2007 John Wiley & Sons, Ltd.     

Prediction of the sawing quality of Marmarit stones using the capability of artificial neural network

The sawing rate is one of the most significant and effective parameters in extracting building stones via diamond wire sawing. This parameter designates the capability of diamond wire sawing for sawing different stones; in addition, the parameter gives rise to economical considerations for quarry designers. In this study, the existent relations between stone geotechnical parameters and the sawing rate of stones via diamond wire sawing were analyzed using regression and correlation coefficient as well as the collected data from Marmarit stone quarries. Moreover, we estimated the sawing rate of Marmarit using the dimensional stone rock mass rating (DSRMR); upon comparison of the data obtained from DSRMR our pre‐collected data on quarries, we did not gain satisfactory results from DSRMR, hence we used artificial neural network (ANN). The results showed that the percentage of Silica, the coefficient of water absorption, the uniaxial compressive strength (UCS), and abrasive hardness are the proper parameters for creating the ANN. Discontinuities have the least effects possible on diamond wire sawing. Having given the training possibility of the ANN, and its ability to evaluate relations among input parameters, the ANN, which was being trained with Marmarit's traits, was an accurate network for estimating diamond wire sawing in Marmarit quarries, although it could not generalize this network for other stones such as Chini and Crystal. Copyright © 2011 John Wiley & Sons, Ltd.     

Geochemical constraints on the genesis of the volcanic rocks in the southeast of Isfahan area, Iran

Late Miocene–Pliocene to Quaternary calc-alkaline lava flows and domes are exposed in southeast of Isfahan in the Urumieh Dokhtar magmatic belt in the Central Iran structural zone. These volcanic rocks have compositions ranging from basaltic andesites, andesites to dacites. Geochemical studies show these rocks are a medium to high K calc-alkaline suite and meta-aluminous. Major element variations are typical for calc-alkaline rocks. The volcanic rocks have SiO₂ contents ranging between 53.8% and 65.3%. Harker diagrams clearly show that the dacitic rocks did not form from the basaltic andesites by normal differentiation processes. They show large ion lithophile elements- and light rare earth elements (LREE)-enriched normalized multielement patterns and negative Nb, Ti, Ta, and P. Condrite-normalized REE patterns display a steep decrease from LREE to light rare earth elements without any Eu anomaly. These characteristics are consistent with ratios obtained from subduction-related volcanic rocks and in collision setting. The melting of a heterogeneous source is possible mechanism for their magma genesis, which was enriched in incompatible elements situated at the upper continental lithospheric mantle or lower crust. The geochemical characteristics of these volcanic rocks suggested that these volcanic rocks evolved by contamination of a parental magma derived from metasomatized upper lithospheric mantle and crustal melts.     

Random vibration of nonlinear structures with stiffness and strength deterioration by modified tail equivalent linearization method

In this research the tail equivalent linearization method(TELM) has been extended to study structures with degrading materials. The responses of such structures to excitations are non-stationary, even if the excitations are stationary. Non-stationary behavior of the system cannot be considered by conventional TELM. Applying the conventional TELM, the only distinction in the design point excitation for two stationary excitations with different durations is in the addition of a zero value part at the beginning of the design point of the longer excitation. This means that the failure probability is the same for the non-stationary systems under excitations with different durations. Therefore, this solution cannot be correct. In this study, in using TELM for systems with degrading materials, hysteretic energy is replaced by average hysteretic energy, calculated by averaging the obtained hysteretic energy of the structure subjected to a few random sample load realization. In this way, the degradation parameters under design point coincide with those under sample load realizations. Since the average of the hysteretic energy is converges very fast, the modified TELM only requires about tens to hundreds solutions of the response in addition to the ordinary calculations of conventional TELM.     

Vertical geochemical variations and arsenic mobilization in the shallow alluvial aquifers of the Chapai-Nawabganj District, northwestern Bangladesh: implication of siderite precipitation

Core sediments from three disturbed boreholes (JOR, GHAT, and RAJ) and two undisturbed boreholes (DW1 and DW2) were collected in the study area of the Chapai-Nawabganj district of northwestern Bangladesh for geochemical analyses. In the study area, groundwater samples from fourteen As-contained private wells and five nested piezometers at both the DW1 and DW2 boreholes were also collected and analyzed. The groundwater arsenic concentrations in the uppermost aquifer (10–40 m of depth) range from 3 to 315 μg/L (mean 47.73 ± 73.41 μg/L), while the arsenic content in sediments range from 2 to 14 mg/kg (mean 4.36 ± 3.34 mg/kg). An environmental scanning electron microscope (ESEM) with an energy dispersive X-ray spectrometer was used to investigate the presence of major and trace elements in the sediments. Groundwaters in the study area are generally the Ca–HCO₃ type with high concentrations of As, but low levels of Fe, Mn, NO₃ ^? and SO ₄ ^?2 . The concentrations of As, Fe, Mn decrease with depth in the groundwater, showing vertical geochemical variations in the study area. Statistical analysis clearly shows that As is closely associated with Fe and Mn in the sediments of the JOR core (r = 0.87, p < 0.05 for Fe and r = 0.78, p < 0.05 for Mn) and GHAT core (r = 0.95, p < 0.05 for Fe and r = 0.93, p < 0.05 for Mn), while As is not correlated with Fe and Mn in groundwater. The comparatively low Fe and Mn concentrations in some groundwater and the ESEM image revealed that siderite precipitated as a secondary mineral on the surface of the sediment particles. The correlations along with results of sequential extraction experiments indicated that reductive dissolution of FeOOH and MnOOH represents a mechanism for releasing arsenic into the groundwater.     

Generating non-stationary random fields of auto-correlated,normally distributed CPT profile by matrix decomposition method

This paper emphasises the true realisation of Cone Penetration Test (CPT) profiles considering non-stationary nature of the data. Formulation of stationary random field theory has been modified and adapted to non-stationary state in order to take into account the mean and variance variability for soil properties. Multi-variance correlation matrix along with the Cholesky decomposition technique was employed to produce realisations of non-homogenous and non-stationary random fields of CPT profiles. A piecewise and segmental data realisation according to the lithology and site class specifications acquired directly from CPT data is adopted in this study so as to render an accurate data simulation. For validation of proposed method 8 CPT test profiles collected from Urmia Lake site have been introduced and simulated by the stationary and non-stationary algorithms. The mean correlation coefficient between the actual CPT data profiles and related realisations along with some other important statistical parameters and their coefficients of variation strongly demonstrate that non-stationary random field generation technique gives quite better accuracy, by comparison to the conventional stationary random field generation scheme.     

Probabilistic rainfall thresholds for flood forecasting: evaluating different methodologies for modelling rainfall spatial correlation (or dependence)

Rainfall threshold (RT) method is one of the evolving flood forecasting approaches. When the cumulative rainfall depth for a given initial soil moisture condition intersects the threshold rainfall curve, the peak discharge is expected to be equal or greater than the threshold discharge for flooding at the target site. Besides the total rainfall depth, spatial and temporal distribution of rainfall impacts the flood peak discharge and the time to peak. To revisit a previous study conducted by the authors, in which spatially independent rainfall pattern was assumed, the spatial distribution of rainfall was simulated following a Monte Carlo approach. The structure of the spatial dependence among sub‐watersheds' rainfalls was taken into account under three different scenarios, namely independent, bivariate copula (2copula) and multivariate Gaussian copula (MGC). For each set of generated random dimensionless rainfalls, the probabilistic RT curves were derived for dry moisture condition. Results were evaluated with both historical and simulated events. For the simulated events, threshold curves were assessed by means of categorical statistics, such as hit rate, false rate and critical success index (CSI). Results revealed that the best performance based on the CSI criterion corresponded to 50% curve in 2copula and MGC scenarios as well as 90% curve in the independent scenario. The recognition of 50% curve in 2copula and MGC scenarios is in agreement with our expectations that the mean probable curve should have the best performance. Moreover, the proposed inclusion of spatially dependent rainfall scenario improved the performance of RT curves by about 25% in comparison with the presumed spatially uniform rainfall scenario. Copyright © 2011 John Wiley & Sons, Ltd.