Elasto-plastic analysis of jointed rocks using discrete continuum and equivalent continuum approaches

Results from elasto-plastic numerical simulations of jointed rocks using both the equivalent continuum and discrete continuum approaches are presented, and are compared with experimental measurements. Initially triaxial compression tests on different types of rocks with wide variation in the uniaxial compressive strength are simulated using both the approaches and the results are compared. The applicability and relative merits and limitations of both the approaches for the simulation of jointed rocks are discussed. It is observed that both the approaches are reasonably good in predicting the real response. However, the equivalent continuum approach has predicted somewhat higher stiffness values at low strains. Considering the modelling effort involved in case of discrete continuum approach, for problems with complex geometry, it is suggested that a proper equivalent continuum model can be used, without compromising much on the accuracy of the results. Then the numerical analysis of a tunnel in Japan is taken up using the continuum approach. The deformations predicted are compared well against the field measurements and the predictions from discontinuum analysis. (C) 2012 Elsevier Ltd. All rights reserved.

Experimental Studies on the Mechanics of Cohesive Frictional Granular Media

The mechanical behaviour of cohesive-frictional granular materials is a combination of the strength pervading as intergranular friction (represented as an angle of internal friction - Phi), and the cohesion (C) between these particles. Most behavioral or constitutive models of this class of granular materials comprise of a cohesion and frictional component with no regard to the length scale i.e. from the micro structural models through the continuum models. An experimental study has been made on a model granular material, viz. angular sand with different weights of binding agents (varying degrees of cohesion) at multiple length scales to physically map this phenomenon. Cylindrical specimen of various diameters - 10, 20, 38, 100, 150 mm (and with an aspect ratio of 2) are reconstituted with 2, 4 and 8% by weight of a binding agent. The magnitude of this cohesion is analyzed using uniaxial compression tests and it is assumed to correspond to the peak in the normalized stress-strain plot. Increase in the cohesive strength of the material is seen with increasing size of the specimen. A possibility of ``entanglement'' occurring in larger specimens is proposed as a possible reason for deviation from a continuum framework.

A spectral multiscale method for wave propagation analysis: Atomistic-continuum coupled simulation

In this paper, we present a new multiscale method which is capable of coupling atomistic and continuum domains for high frequency wave propagation analysis. The problem of non-physical wave reflection, which occurs due to the change in system description across the interface between two scales, can be satisfactorily overcome by the proposed method. We propose an efficient spectral domain decomposition of the total fine scale displacement along with a potent macroscale equation in the Laplace domain to eliminate the spurious interfacial reflection. We use Laplace transform based spectral finite element method to model the macroscale, which provides the optimum approximations for required dynamic responses of the outer atoms of the simulated microscale region very accurately. This new method shows excellent agreement between the proposed multiscale model and the full molecular dynamics (MD) results. Numerical experiments of wave propagation in a 1D harmonic lattice, a 1D lattice with Lennard-Jones potential, a 2D square Bravais lattice, and a 2D triangular lattice with microcrack demonstrate the accuracy and the robustness of the method. In addition, under certain conditions, this method can simulate complex dynamics of crystalline solids involving different spatial and/or temporal scales with sufficient accuracy and efficiency. (C) 2014 Elsevier B.V. All rights reserved.

Notch sensitivity in nanoscale metallic glass specimens: Insights from continuum simulations

Recent experiments have shown that nano-sized metallic glass (MG) specimens subjected to tensile loading exhibit increased ductility and work hardening. Failure occurs by necking as opposed to shear banding which is seen in bulk samples. Also, the necking is generally observed at shallow notches present on the specimen surface. In this work, continuum finite element analysis of tensile loading of nano-sized notched MG specimens is conducted using a thermodynamically consistent non-local plasticity model to clearly understand the deformation behavior from a mechanics perspective. It is found that plastic zone size in front of the notch attains a saturation level at the stage when a dominant shear band forms extending across the specimen. This size scales with an intrinsic material length associated with the interaction stress between flow defects. A transition in deformation behavior from quasi-brittle to ductile becomes possible when this critical plastic zone size is larger than the uncracked ligament length. These observations corroborate with atomistic simulations and experimental results. (C) 2015 Elsevier Ltd. All rights reserved.

Axially Loaded Timoshenko Rotors From A Prestressed Continuum Approach

We consider an axially loaded Timoshenko rotor rotating at a constant speed and derive its governing equations from a continuum viewpoint. The primary aim of this paper is to understand the source and role of gyroscopic terms, when the rotor is viewed not as a Timoshenko beam but as a genuine 3D continuum. We offer the primary insight that macroscopically observed gyroscopic terms may also, quite equivalently, be viewed as external manifestations of internally existing spin-induced prestresses at the continuum level. To demonstrate this idea with an analytical example (the Timoshenko rotor), we have studied the reliable equations of Choi et al. (Journal of Vibration and Acoustics, 114, 1992, 249-259). Using a straightforward application of our insight in the framework of nonlinear elasticity, we obtain equations that exactly match Choi et al. for the case with no axial load. For the case of axial preload, our straightforward formulation leads to a slightly different set of equations that have negligible numerical consequence for solid rotors. However, we offer a macroscopic, intuitive, justification for modifying our formulation so as to obtain the exact equations of Choi et al. with the axial load included.

Constraints in relativistic Hamiltonian mechanics

We propose a new scheme for the use of constraints in setting up classical, Hamiltonian, relativistic, interacting particle theories. We show that it possesses both Poincaré invariance and invariance of world lines. We discuss the transition to the physical phase space and the nonrelativistic limit.

On global summer monsoonal drought mechanics

It is shown that within the framework of a linear five-level quasi-geostrophic steady state global model the middle latitude systems can always have significant influence on the Asian summer monsoonal system through the lower tropospheric monsoonal westerly window region around 80°E. It is hypothesized that quasistationarity of the middle latitude longwave systems results in stronger teleconnections through this window and the consequent monsoon breaks when the phase is right.

On a linear steady-state system of equations in microcontinuum fluid mechanics

Galerkin representations and integral representations are obtained for the linearized system of coupled differential equations governing steady incompressible flow of a micropolar fluid. The special case of 2-dimensional Stokes flows is then examined and further representation formulae as well as asymptotic expressions, are generated for both the microrotation and velocity vectors. With the aid of these formulae, the Stokes Paradox for micropolar fluids is established.

Analytical estimation of stress intensity factors in patched cracked plates

The fatigue and fracture performance of a cracked plate can be substantially improved by providing patches as reinforcements. The effectiveness of the patches is related to the reduction they cause in the stress intensity factor (SIF) of the crack. So, for reliable design, one needs an accurate evaluation of the SIF in terms of the crack, patch and adhesive parameters. In this investigation, a centrally cracked large plate with a pair of symmetric bonded narrow patches, oriented normally to the crack line, is analysed by a continuum approach. The narrow patches are treated as transversely flexible line members. The formulation leads to an integral equation which is solved numerically using point collocation. The convergence is rapid. It is found that substantial reductions in SIF are possible with practicable patch dimensions and locations. The patch is more effective when placed on the crack than ahead of the crack. The present analysis indicates that a little distance inwards of the crack tip, not the crack tip itself, is the ideal location, for the patch.

Continuum percolation with discs having a distribution of radii

It is shown that for continuum percolation with overlapping discs having a distribution of radii, the net areal density of discs at percolation threshold depends non-trivially on the distribution, and is not bounded by any finite constant. Results of a Monte Carlo simulation supporting the argument are presented.

Charged-soliton excitations in statistical mechanics

A method is developed for demonstrating how solitons with some internal periodic motion may emerge as elementary excitations in the statistical mechanics of field systems. The procedure is demonstrated in the context of complex scalar fields which can, for appropriate choices of the Lagrangian, yield charge-carrying solitons with such internal motion. The derivation uses the techniques of the steepest-descent method for functional integrals. It is shown that, despite the constraint of some fixed total charge, a gaslike excitation of such charged solitons does emerge.

Non-Abelian monopoles break color. I. Classical mechanics

Monopoles which are sources of non-Abelian magnetic flux are predicted by many models of grand unification. It has been argued elsewhere that a generic transformation of the "unbroken" symmetry group H cannot be globally implemented on such monopoles for reasons of topology. In this paper, we show that similar topological obstructions are encountered in the mechanics of a test particle in the field of these monopoles and that the transformations of H cannot all be globally implemented as canonical transformations. For the SU(5) model, if H is SU(3)C×U(1)em, a consequence is that color multiplets are not globally defined, while if H is SU(3)C×SU(2)WS×U(1)Y, the same is the case for both color and electroweak multiplets. There are, however, several subgroups KT, KT′,… of H which can be globally implemented, with the transformation laws of the observables differing from group to group in a novel way. For H=SU(3)C×U(1)em, a choice for KT is SU(2)C×U(1)em, while for H=SU(3)C×SU(2)WS×U(1)Y, a choice is SU(2)C×U(1)×U(1)×U(1). The paper also develops the differential geometry of monopoles in a form convenient for computations.