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.

Nonlocal continuum mechanics formulation for axial, flexural, shear and contraction coupled wave propagation in single walled carbon nanotubes

This paper presents the effect of nonlocal scaling parameter on the coupled i.e., axial, flexural, shear and contraction, wave propagation in single-walled carbon nanotubes (SWCNTs). The axial and transverse motion of SWCNT is modeled based on first order shear deformation theory (FSDT) and thickness contraction. The governing equations are derived based on nonlocal constitutive relations and the wave dispersion analysis is also carried out. The studies shows that the nonlocal scale parameter introduces certain band gap region in all wave modes where no wave propagation occurs. This is manifested in the wavenumber plots as the region where the wavenumber tends to infinite or wave speed tends to zero. The frequency at which this phenomenon occurs is called the escape frequency. Explicit expressions are derived for cut-off and escape frequencies of all waves in SWCNT. It is also shown that the cut-off frequencies of shear and contraction mode are independent of the nonlocal scale parameter. The results provided in this article are new and are useful guidance for the study and design of the next generation of nanodevices that make use of the coupled wave propagation properties of single-walled carbon nanotubes.

Thermal vibration analysis of orthotropic nanoplates based on nonlocal continuum mechanics

This paper presents the thermal vibration analysis of orthotropic nanoplates such as graphene, using the two variable refined plate theory and nonlocal continuum mechanics for small scale effects. The nanoplate is modeled based on two variable refined plate theory and the axial stress caused by the thermal effects is also considered. The two variable refined plate theory takes account of transverse shear effects and parabolic distribution of the transverse shear strains through the thickness of the plate, hence it is unnecessary to use shear correction factors. Nonlocal governing equations of motion for the nanoplate are derived from the principle of virtual displacements. The closed form solution for thermal-vibration frequencies of a simply supported rectangular nanoplate has been obtained by using Navier's method of solution. Numerical results obtained by the present theory are compared with available solutions in the literature and the molecular dynamics results. The influences of the small scale coefficient, the room or low temperature, the high temparature, the half wave number and the aspect ratio of nanoplate on the natural frequencies are considered and discussed in detail. It can be concluded that the present theory, which does not require shear correction factor, is not only simple but also comparable to the first-order and higher order shear deformation theory. The present analysis results can be used for the design of the next generation of nanodevices that make use of the thermal vibration properties of the nanoplates. (C) 2012 Elsevier B.V. All rights reserved.

Continuum theory of edge states of topological insulators: variational principle and boundary conditions

We develop a continuum theory to model low energy excitations of a generic four-band time reversal invariant electronic system with boundaries. We propose a variational energy functional for the wavefunctions which allows us to derive natural boundary conditions valid for such systems. Our formulation is particularly suited for developing a continuum theory of the protected edge/surface excitations of topological insulators both in two and three dimensions. By a detailed comparison of our analytical formulation with tight binding calculations of ribbons of topological insulators modelled by the Bernevig-Hughes-Zhang (BHZ) Hamiltonian, we show that the continuum theory with a natural boundary condition provides an appropriate description of the low energy physics.

Nonlocal continuum mechanics based ultrasonic flexural wave dispersion characteristics of a monolayer graphene embedded in polymer matrix

Wave propagation in graphene sheet embedded in elastic medium (polymer matrix) has been a topic of great interest in nanomechanics of graphene sheets, where the equivalent continuum models are widely used. In this manuscript, we examined this issue by incorporating the nonlocal theory into the classical plate model. The influence of the nonlocal scale effects has been investigated in detail. The results are qualitatively different from those obtained based on the local/classical plate theory and thus, are important for the development of monolayer graphene-based nanodevices. In the present work, the graphene sheet is modeled as an isotropic plate of one-atom thick. The chemical bonds are assumed to be formed between the graphene sheet and the elastic medium. The polymer matrix is described by a Pasternak foundation model, which accounts for both normal pressure and the transverse shear deformation of the surrounding elastic medium. When the shear effects are neglected, the model reduces to Winkler foundation model. The normal pressure or Winkler elastic foundation parameter is approximated as a series of closely spaced, mutually independent, vertical linear elastic springs where the foundation modulus is assumed equivalent to stiffness of the springs. For this model, the nonlocal governing differential equations of motion are derived from the minimization of the total potential energy of the entire system. An ultrasonic type of flexural wave propagation model is also derived and the results of the wave dispersion analysis are shown for both local and nonlocal elasticity calculations. From this analysis we show that the elastic matrix highly affects the flexural wave mode and it rapidly increases the frequency band gap of flexural mode. The flexural wavenumbers obtained from nonlocal elasticity calculations are higher than the local elasticity calculations. The corresponding wave group speeds are smaller in nonlocal calculation as compared to local elasticity calculation. The effect of y-directional wavenumber (eta(q)) on the spectrum and dispersion relations of the graphene embedded in polymer matrix is also observed. We also show that the cut-off frequencies of flexural wave mode depends not only on the y-direction wavenumber but also on nonlocal scaling parameter (e(0)a). The effect of eta(q) and e(0)a on the cut-off frequency variation is also captured for the cases of with and without elastic matrix effect. For a given nanostructure, nonlocal small scale coefficient can be obtained by matching the results from molecular dynamics (MD) simulations and the nonlocal elasticity calculations. At that value of the nonlocal scale coefficient, the waves will propagate in the nanostructure at that cut-off frequency. In the present paper, different values of e(0)a are used. One can get the exact e(0)a for a given graphene sheet by matching the MD simulation results of graphene with the results presented in this article. (c) 2012 Elsevier Ltd. All rights reserved.

Thermal vibration analysis of monolayer graphene embedded in elastic medium based on nonlocal continuum mechanics

This paper presents the thermal vibration analysis of single-layer graphene sheet embedded in polymer elastic medium, using the plate theory and nonlocal continuum mechanics for small scale effects. The graphene is modeled based on continuum plate theory and axial stress caused by the thermal effects is also considered. Nonlocal governing equation of motion for this graphene sheet system is derived from the principle of virtual displacements. The closed form solution for thermal-vibration frequencies of a simply supported rectangular nanoplate has been obtained by using the Navier's method of solution. Numerical results obtained by the present theory are compared with available solutions in the literature and the molecular dynamics results. The influences of the small scale coefficient, the room or low temperature, the high temperature, the half wave number and the aspect ratio of nanoplate on the natural frequencies are considered and discussed in detail. The thermal vibration analysis of single- and double-layer graphene sheets are considered for the analysis. The mode shapes of the respective graphene system are also captured in this work. The present analysis results can be used for the design of the next generation of nanodevices that make use of the thermal vibration properties of the graphene.

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.

Percolation on networks with antagonistic and dependent interactions

Drawing inspiration from real world interacting systems, we study a system consisting of two networks that exhibit antagonistic and dependent interactions. By antagonistic and dependent interactions we mean that a proportion of functional nodes in a network cause failure of nodes in the other, while failure of nodes in the other results in failure of links in the first. In contrast to interdependent networks, which can exhibit first-order phase transitions, we find that the phase transitions in such networks are continuous. Our analysis shows that, compared to an isolated network, the system is more robust against random attacks. Surprisingly, we observe a region in the parameter space where the giant connected components of both networks start oscillating. Furthermore, we find that for Erdos-Renyi and scale-free networks the system oscillates only when the dependence and antagonism between the two networks are very high. We believe that this study can further our understanding of real world interacting systems.

Unified Continuum Damage Model for Matrix Cracking in Composite Rotor Blades

This paper deals with modeling of the first damage mode, matrix micro-cracking, in helicopter rotor/wind turbine blades and how this effects the overall cross-sectional stiffness. The helicopter/wind turbine rotor system operates in a highly dynamic and unsteady environment leading to severe vibratory loads present in the system. Repeated exposure to this loading condition can induce damage in the composite rotor blades. These rotor/turbine blades are generally made of fiber-reinforced laminated composites and exhibit various competing modes of damage such as matrix micro-cracking, delamination, and fiber breakage. There is a need to study the behavior of the composite rotor system under various key damage modes in composite materials for developing Structural Health Monitoring (SHM) system. Each blade is modeled as a beam based on geometrically non-linear 3-D elasticity theory. Each blade thus splits into 2-D analyzes of cross-sections and non-linear 1-D analyzes along the beam reference curves. Two different tools are used here for complete 3-D analysis: VABS for 2-D cross-sectional analysis and GEBT for 1-D beam analysis. The physically-based failure models for matrix in compression and tension loading are used in the present work. Matrix cracking is detected using two failure criterion: Matrix Failure in Compression and Matrix Failure in Tension which are based on the recovered field. A strain variable is set which drives the damage variable for matrix cracking and this damage variable is used to estimate the reduced cross-sectional stiffness. The matrix micro-cracking is performed in two different approaches: (i) Element-wise, and (ii) Node-wise. The procedure presented in this paper is implemented in VABS as matrix micro-cracking modeling module. Three examples are presented to investigate the matrix failure model which illustrate the effect of matrix cracking on cross-sectional stiffness by varying the applied cyclic

Selective localisation of multi walled carbon nanotubes in polypropylene/natural rubber blends to reduce the percolation threshold

Polypropylene and natural rubber blends with multiwalled carbon nanotube (PP/NR + MWCNT nanocomposites) were prepared by melt mixing. The melt rheological behaviour of neat PP and PP/NR blends filled with different loadings (1, 3, 5, 7 wt%) of MWCNT was studied. The effect of PP/NR blends (with compositions, 80/20,50/50, 20/80 by wt) on the rheological percolation threshold was investigated. It was found that blending PP with NR (80/20 and 50/50 composition) reduced the rheological percolation threshold from 5 wt% to 3 wt% MWCNT. The melt rheological behaviour of the MWCNT filled PP/NR blends was correlated with the morphology observations from high resolution transmission electron microscopic (HRTEM) images. In predicting the thermodynamically favoured location of MWCNT in PP/NR blend, the specific interaction of phospholipids in NR phase with MWCNTs was considered quantitatively. The MWCNTs were selectively localised in the NR phase. The percolation mechanism in MWCNT filled PP/NR blends was discussed and for each blend composition, the percolation mechanism was found to be different. (C) 2015 Elsevier Ltd. All rights reserved.

Internal noise-driven generalized Langevin equation from a nonlocal continuum model

Starting with a micropolar formulation, known to account for nonlocal microstructural effects at the continuum level, a generalized Langevin equation (GLE) for a particle, describing the predominant motion of a localized region through a single displacement degree of freedom, is derived. The GLE features a memory-dependent multiplicative or internal noise, which appears upon recognizing that the microrotation variables possess randomness owing to an uncertainty principle. Unlike its classical version, the present GLE qualitatively reproduces the experimentally measured fluctuations in the steady-state mean square displacement of scattering centers in a polyvinyl alcohol slab. The origin of the fluctuations is traced to nonlocal spatial interactions within the continuum, a phenomenon that is ubiquitous across a broad class of response regimes in solids and fluids. This renders the proposed GLE a potentially useful model in such cases.

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.

Continuum in the X-Z-Y Weak Bonds: Z= Main Group Elements

The Continuum in the variation of the X-Z bond length change from blue-shifting to red-shifting through zero-shifting in the X-Z---Y complex is inevitable. This has been analyzed by ab-initio molecular orbital calculations using Z= Hydrogen, Halogens, Chalcogens, and Pnicogens as prototypical examples. Our analysis revealed that, the competition between negative hyperconjugation within the donor (X-Z) molecule and Charge Transfer (CT) from the acceptor (Y) molecule is the primary reason for the X-Z bond length change. Here, we report that, the proper tuning of X-and Y-group for a particular Z-can change the blue-shifting nature of X-Z bond to zero-shifting and further to red-shifting. This observation led to the proposal of a continuum in the variation of the X-Z bond length during the formation of X-Z---Y complex. The varying number of orbitals and electrons available around the Z-atom differentiates various classes of weak interactions and leads to interactions dramatically different from the H-Bond. Our explanations based on the model of anti-bonding orbitals can be transferred from one class of weak interactions to another. We further take the idea of continuum to the nature of chemical bonding in general. (C) 2015 Wiley Periodicals, Inc.

Orientation-based Continuum Damage Models for Rocks

A general formulation of the Helmholtz free energy used in thermodynamics of damage process of rocks is derived within a multi-scale framework. Such a physically-based thermodynamic state potential has a hybrid, discrete/continuum, nature in the sense tha

Stochastic Structural Model of Rock and Soil Aggregates by Continuum-Based Discrete Element Method

This paper first presents a stochastic structural model to describe the random geometrical features of rock and soil aggregates. The stochastic structural model uses mixture ratio, rock size and rock shape to construct the microstructures of aggregates,and introduces two types of structural elements (block element and jointed element) and three types of material elements (rock element, soil element, and weaker jointed element)for this microstructure. Then, continuum-based discrete element method is used to study the deformation and failure mechanism of rock and soil aggregate through a series of loading tests. It is found that the stress-strain curve of rock and soil aggregates is nonlinear, and the failure is usually initialized from weaker jointed elements. Finally, some factors such as mixture ratio, rock size and rock shape are studied in detail. The numerical results are in good agreement with in situ test. Therefore, current model is effective for simulating the mechanical behaviors of rock and soil aggregates.