Fully coupled elastoplastic hydro-mechanical analysis of unsaturated porous media using a meshfree method

This paper presents the first application of an advanced meshfree method, ie, the edge-based smoothed point interpolation method (ESPIM), in simulation of the coupled hydro-mechanical behaviour of unsaturated porous media. In the proposed technique, the problem domain is spatially discretised using a triangular background mesh, and the polynomial point interpolation method combined with a simple node selection scheme is adopted for creating nodal shape functions. Smoothing domains are formed on top of the background mesh, and a constant smoothed strain, created by applying the smoothing operation over the smoothing domains, is assigned to each smoothing domain. The deformation and flow models are developed based on the equilibrium equation of the mixture, and linear momentum and mass balance equations of the fluid phases, respectively. The effective stress approach is followed to account for the coupling between the flow and deformation models. Further coupling among the phases is captured through a hysteretic soil water retention model that evolves with changes in void ratio. An advanced elastoplastic constitutive model within the context of the bounding surface plasticity theory is employed for predicting the nonlinear behaviour of soil skeleton. Time discretisation is performed by adopting a three-point discretisation method with growing time steps to avoid temporal instabilities. A modified Newton-Raphson framework is designed for dealing with nonlinearities of the discretised system of equations. The performance of the numerical model is examined through a number of numerical examples. The state-of-the-art computational scheme developed is useful for simulation of geotechnical engineering problems involving unsaturated soils.     

Correction to: Post-earthquake response by small UAV helicopters

Natural Hazards - The original publication has the following mistakes which we hereby correct: In Fig. 3 the “UAV routs” should be replaced by “UAV routes”. On page...     

Numerical modeling of long waves in shallow water using Incremental Differential Quadrature Method

Incremental Differential Quadrature Method (IDQM) as a rapid and accurate method for numerical simulation of Nonlinear Shallow Water (NLSW) waves is employed. To the best of authors’ knowledge, this is the first endeavor to exploit DQM in coastal hydraulics. The one-dimensional NLSW equations and related boundary conditions are discretized in space and temporal directions by DQM rules and the resulting system of equations are used to compute the state variables in the entire computational domain. It was found that the splitting of total simulation time into a number of smaller time increments, could significantly enhance the performance of the proposed method. Furthermore, results of this study show two main advantages for IDQM compared with other conventional methods, namely; unconditional stability and minimal computational effort. Indeed, using IDQM, one can use a few grid points (in spatial or time direction) without imposing any stability condition on the time step to obtain an accurate convergent solution.     

Tensionless–frictionless interaction of flexible annular foundation with a transversely isotropic multi‐layered half‐space

A transversely isotropic multi‐layered half‐space, with axis of material symmetry perpendicular to the free surface, supports a flexible either annular or solid circle foundation. The contact area of the foundation and the half‐space is considered to be both frictionless and tensionless. The foundation is assumed to be affected by a vertical static axisymmetric load. Detailed analysis of the interaction of these two systems with different thickness of layers is the target of this paper. With the use of ring load Green's functions for both the foundation and the continuum half‐space, an integral equation accompanied with some inequalities is introduced to model the complex BVP. With the incorporation of ring‐shape FEM, we are capable of capturing both regular and singular solution smoothly. The validity of the combination of the analytical and numerical method is proved with comparing the results of this paper with a number of benchmark cases of both linear and nonlinear interaction of circular and annular foundation with half‐space. Some new illustrations are presented to portray the aspect of the anisotropy and layering of the half‐space. Copyright © 2014 John Wiley & Sons, Ltd.     

Steady critical discharge rates from vertical and horizontal wells in fresh–saline aquifers with sharp interfaces

In this study, the critical (or maximum) discharge rates before saline water enters a well were determined for vertical and horizontal wells in a freshwater aquifer which is separated from a static saline aquifer by a sharp interface. Flow around the well was solved by integration of a point sink solution along the well axis, and both the critical discharge rate and critical interface rise were determined through a comparison of the heads and vertical gradients at the saline–fresh water interface. The rates were determined for vertical and horizontal wells with various lengths and depths for different aquifer salinities. Results were generalized by drawing dimensionless type curves. The results showed that the dimensionless total critical discharge rates are higher for the longer horizontal wells and longer vertical wells with a certain bottom depth, and they almost linearly decrease with well depth at rates of 0.7–0.9. For the dimensionless well length of 0.2, the dimensionless total discharge rate of a horizontal well is about 0.1 more than that of a vertical well with the same length and well-top depth. Also, the critical discharge rates per unit length of well are inversely proportional to well length and remarkably higher for shallower wells. Additionally, the critical pumping rate is proportional to the salinity difference of the aquifers. These results were confirmed by comparison to existing solutions for vertical wells with dimensionless lengths of 0.2, 0.5 and 0.6, and for critical interface rises in the range of 0.75–0.9.     

Application of Discriminant Analysis and Support Vector Machine in Mapping Gold Potential Areas for Further Drilling in the Sari-Gunay Gold Deposit,NW Iran

In this contribution, we used discriminant analysis (DA) and support vector machine (SVM) to model subsurface gold mineralization by using a combination of the surface soil geochemical anomalies and earlier bore data for further drilling at the Sari-Gunay gold deposit, NW Iran. Seventy percent of the data were used as the training data and the remaining 30 % were used as the testing data. Sum of the block grades, obtained by kriging, above the cutoff grade (0.5 g/t) was multiplied by the thickness of the blocks and used as productivity index (PI). Then, the PI variable was classified into three classes of background, medium, and high by using fractal method. Four classification functions of SVM and DA methods were calculated by the training soil geochemical data. Also, by using all the geochemical data and classification functions, the general extension of the gold mineralized zones was predicted. The mineral prediction models at the Sari-Gunay hill were used to locate high and moderate potential areas for further infill systematic and reconnaissance drilling, respectively. These models at Agh-Dagh hill and the area between Sari-Gunay and Agh-Dagh hills were used to define the moderate and high potential areas for further reconnaissance drilling. The results showed that the nu-SVM method with 73.8 % accuracy and c-SVM with 72.3 % accuracy worked better than DA methods.     

Analysis of soil nail walls performance - Case study

This study surveys the performance of temporary soil nail walls with a maximum height of 29.3 m in the Yas project, located in Tehran. Some numerical models, with various modelling approaches, were developed using finite element software and the proper modelling approaches were specified and verified.The modelling results were analysed and compared with the in situ monitored data. The results demonstrate the model’s horizontal deformations are generally greater than the in situ monitored values. To decrease the differences between models and monitored results, the effect of variations of the soil parameters had been surveyed and the limits of needed variations were specified.     

A near-fault probabilistic seismic hazard assessment for Yasouj,located in the Kazerun fault system,southwest Iran

Natural Hazards - The Kazerun fault system (KFS) is located in the central part of the Zagros, the most seismically active orogenic belt in the Iranian plateau. The city of Yasouj is located in the...