2D stabilised analytic signal method in DC pole-pole potential data interpretation

Using analytic signal method, interpretation of pole-pole secondary electric potentials due to 2D conductive/resistive prisms is presented. The estimated parameters are the location, lateral extent or width and depth to top surface of the prism. Forward modelling is attempted by 2D-Finite Difference method. The proposed stabilised analytic signal algorithm (RES2AS) uses Tikhonov’s regularization scheme and FFT routines. The algorithm is tested on three theoretical examples and field data from the campus of Roorkee University. The stability of RES2AS is also tested on synthetic error prone secondary pole-pole potential data.     

Numerically observed shear bands in soft sensitive clays

The numerical challenges that arise in modelling shear bands in soft sensitive (SS) clays have not yet been fully resolved. Convincing and well-accepted solutions have yet to be found. This paper presents some novel information related to the shear band phenomenon in SS clays. In this study, the hypothesis is that the generation and dissipation of excess pore pressure from shear bands could regularise the strain softening and result in a mesh independent shear band thickness. The generation and dissipation of excess pressure is modelled by a coupled consolidation process. The simulation aims at modelling two counteracting mechanisms in the SS clay. First, the shear band narrows because of strain softening. Second, the internal pore water pressure drainage reduces the rate of strain softening. This counteracting mechanism provides an inherent regularisation technique for SS clays. This study presents some numerical results involving these two counteracting mechanisms. This study also shows that an inherent internal parameter applicable for SS clays can be defined by the ratio between soil permeability and the applied strain rate. In the case of SS clays, the range of this parameter varies from 0 to 0.0002 mm.     

Structural softening,mesh dependence,and regularisation in non-associated plastic flow

A severe dependence of numerical simulations on the mesh density is usually attributed to the presence of strain softening in the constitutive relation. However, other material instabilities, like non-associated plastic flow, can also cause mesh sensitivity. Indeed, loss of ellipticity in quasi-static analyses is the fundamental cause of the observed mesh dependence. It has been known since long that non-associated plastic flow can cause loss of ellipticity, but the consequence for mesh sensitivity, and subsequently, for the difficulty of the equilibrium-finding iterative procedure to converge have remained largely unnoticed. We first demonstrate at the hand of a biaxial test structural softening and a marked mesh dependence for an ideally plastic material equipped with a non-associated flow rule. The phenomena are then analysed in depth using an infinitely long shear layer. Finally, it is shown that the mesh effect disappears when the standard continuum model is replaced by a Cosserat continuum, a well-known regularisation method for strain-softening constitutive relations.     

Periodic and Chaotic Trajectories of the Second Species for the n-Centre Problem

For the n-centre problem of one particle moving in the potential of attracting centres of small mass fixed in an arbitrary smooth potential and magnetic field, we prove the existence of periodic and chaotic trajectories shadowing sequences of collision orbits. In particular, we obtain large subshifts of solutions of this type for the circular restricted 3-body problem of celestial mechanics. Poincaré had conjectured existence of the periodic ones and given them the name ‘second species solutions’. This revised version was published online in July 2006 with corrections to the Cover Date.     

Explicit Adaptive Symplectic (Easy) Integrators: A Scaling Invariant Generalisation of the Levi-Civitaand KS Regularisations

We present a generalisation of the Levi-Civita and Kustaanheimo-Stiefel regularisation. This allows the use of more general time rescalings. In particular, it is possible to find a regularisation which removes the singularity of the equations and preserves scaling invariance. In addition, these equations can, in certain cases, be integrated with explicit symplectic Runge-Kutta-Nyström methods. The combination of both techniques gives an explicit adaptive symplectic (EASY) integrator. We apply those methods to some perturbations of the Kepler problem and illustrate, by means of some numerical examples, when scaling invariant regularisations are more efficient that the LC/KS regularisation.     

Simulation of mixed-mode fracture using SPH particles with an embedded fracture process zone

This paper focuses on the modelling of mixed-mode fracture using the conventional smoothed particle hydrodynamics (SPH) method and a mixed-mode cohesive fracture law embedded in the particles. The combination of conventional SPH and a mixed-mode cohesive model allows capturing fracture and separation under various loading conditions efficiently. The key advantage of this framework is its capability to represent complex fracture geometries by a set of cracked SPH particles, each of which can possess its own mixed-mode cohesive fracture with arbitrary orientations. Therefore, this can naturally capture complex fracture patterns without any predefined fracture topologies. Because a characteristic length scale related to the size of the fracture process zone is incorporated in the constitutive formulation, the proposed approach is independent from the spatial discretisation of the computational domain (or mesh independent). Furthermore, the anisotropic fracture responses of materials can be naturally captured thanks to the orientation of the fracture process zone embedded at the particle level. The performance of the proposed approach demonstrates its potentials in modelling mixed-mode fracture of rocks and similar quasi-brittle materials.     

From discrete to continuum modelling of boundary value problems in geomechanics: An integrated FEM-DEM approach

Double-scale numerical methods constitute an effective tool for simultaneously representing the complex nature of geomaterials and treating real-scale engineering problems such as a tunnel excavation or a pressuremetre at a reasonable numerical cost. This paper presents an approach coupling discrete elements (DEM) at the microscale with finite elements (FEM) at the macroscale. In this approach, a DEM-based numerical constitutive law is embedded into a standard FEM formulation. In this regard, an exhaustive discussion is presented on how a 2D/3D granular assembly can be used to generate, step by step along the overall computation process, a consistent Numerically Homogenised Law. The paper also focuses on some recent developments including a comprehensive discussion of the efficiency of Newton-like operators, the introduction of a regularisation technique at the macroscale by means of a second gradient framework, and the development of parallelisation techniques to alleviate the computational cost of the proposed approach. Some real-scale problems taking into account the material spatial variability are illustrated, proving the numerical efficiency of the proposed approach and the benefit of a particle-based strategy.     

Predictive potential of Perzyna viscoplastic modelling for granular geomaterials

This paper reappraises Perzyna-type viscoplasticity for the constitutive modelling of granular geomaterials, with emphasis on the simulation of rate/time effects of different magnitude. An existing elasto-plastic model for sands is first recast into a Perzyna viscoplastic formulation and then calibrated/validated against laboratory test results on Hostun sand from the literature. Notable model features include (1) enhanced definition of the viscous nucleus function and (2) void ratio dependence of stiffness and viscous parameters, to model the pycnotropic behaviour of granular materials with a single set of parameters, uniquely identified against standard creep and triaxial test results. The comparison between experimental data and numerical simulations points out the predicative capability of the developed model and the complexity of defining a unique viscous nucleus function to capture sand behaviour under different loading/initial/boundary and drainage conditions. It is concluded that the unified viscoplastic simulation of both drained and undrained response is particularly challenging within Perzyna's framework and opens to future research in the area. The discussion presented is relevant, for instance, to the simulation of multiphase strain localisation phenomena, such as those associated to slope stability problems in variably saturated soils.     

On viscoplastic regularisation of strain-softening rocks and soils

Constitutive models for rocks and soils typically incorporate some form of strain softening. Moreover, many plasticity models for frictional materials use a nonassociated flow rule. Strain softening and nonassociated flow rules can cause loss of well-posedness of the initial-value problem, which can lead to a severe mesh dependence in simulations and poor convergence of the iterative solution procedure. The inclusion of viscosity, which is a common property of materials, seems a natural way to restore well-posedness, but the mathematical properties of a rate-dependent model, and therefore the effectiveness with respect to the removal of mesh dependence, can depend strongly on how the viscous element is incorporated. Herein, we show that rate-dependent models, which are commonly applied to problems in the Earth's lithosphere, such as plate tectonics, are very different from the approach typically adopted for more shallow geotechnical engineering problems. We analyse the properties of these models under dynamic loadings, using dispersion analyses and one-dimensional finite difference analyses, and complement them with two-dimensional simulations of a typical strain localisation problem under quasi-static loading conditions. Finally, we point out that a combined model, which features two viscous elements, may be the best way forward for modelling time-dependent failure processes in the deeper layers of the Earth, since it not only enables modelling of the creep characteristics typical of long-term behaviour but also regularises the initial/boundary-value problem.     

Load or concentration,logged or unlogged? Addressing ten years of uncertainty in neural network suspended sediment prediction

This paper addresses current inconsistencies in methodological approaches for neural network modelling of suspended sediment. An expansion in the number of case studies being published over the last decade has yet to result in agreed guidelines on whether suspended sediment load or concentration should be modelled, and whether log‐transformation of data is either necessary or potentially beneficial. This contrasts with the well‐recognized guidelines that direct traditional sediment rating curve studies. The paper reports a comprehensive set of single‐input single‐output neural network suspended sediment modelling experiments performed on two catchments in Puerto Rico. It examines the impact of internal complexity, input variable choice and data transformation on the form, consistency and physical rationality of model outputs, the existence of localized overfitting and the usefulness of global performance metrics. Sound guidance on whether to model sediment load or concentration, and whether to model log‐transformed data is provided. Copyright © 2011 John Wiley & Sons, Ltd.