An adaptive singular ES-FEM for mechanics problems with singular field of arbitrary order

This paper presents a singular edge-based smoothed finite element method (sES-FEM) for mechanics problems with singular stress fields of arbitrary order. The sES-FEM uses a basic mesh of three-noded linear triangular (T3) elements and a special layer of five-noded singular triangular elements (sT5) connected to the singular-point of the stress field. The sT5 element has an additional node on each of the two edges connected to the singular-point. It allows us to represent simple and efficient enrichment with desired terms for the displacement field near the singular-point with the satisfaction of partition-of-unity property. The stiffness matrix of the discretized system is then obtained using the assumed displacement values (not the derivatives) over smoothing domains associated with the edges of elements. An adaptive procedure for the sES-FEM is proposed to enhance the quality of the solution with minimized number of nodes. Several numerical examples are provided to validate the reliability of the present sES-FEM method. (C) 2012 Elsevier B.V. All rights reserved.

Cavitation in materials with distributed weak zones: Implications on the origin of brittle fracture in metallic glasses

Several experimental studies have shown that fracture surfaces in brittle metallic glasses (MGs) generally exhibit nanoscale corrugations which may be attributed to the nucleation and coalescence of nanovoids during crack propagation. Recent atomistic simulations suggest that this phenomenon is due to large spatial fluctuations in material properties in a brittle MG, which leads to void nucleation in regions of low atomic density and then catastrophic fracture through void coalescence. To explain this behavior, we propose a model of a heterogeneous solid containing a distribution of weak zones to represent a brittle MG. Plane strain continuum finite element analysis of cavitation in such an elastic-plastic solid is performed with the weak zones idealized as periodically distributed regions having lower yield strength than the background material. It is found that the presence of weak zones can significantly reduce the critical hydrostatic stress for the onset of cavitation which is controlled uniquely by the local yield properties of these zones. Also, the presence of weak zones diminishes the sensitivity of the cavitation stress to the volume fraction of a preexisting void. These results provide plausible explanations for the observations reported in recent atomistic simulations of brittle MGs. An analytical solution for a composite, incompressible elastic-plastic solid with a weak inner core is used to investigate the effect of volume fraction and yield strength of the core on the nature of cavitation bifurcation. It is shown that snap-cavitation may occur, giving rise to sudden formation of voids with finite size, which does not happen in a homogeneous plastic solid. (c) 2012 Elsevier Ltd. All rights reserved.

Reaction dynamics under confinement: an exact path integral treatment of a two-stage model of stochastic gene expression

Gene expression in living systems is inherently stochastic, and tends to produce varying numbers of proteins over repeated cycles of transcription and translation. In this paper, an expression is derived for the steady-state protein number distribution starting from a two-stage kinetic model of the gene expression process involving p proteins and r mRNAs. The derivation is based on an exact path integral evaluation of the joint distribution, P(p, r, t), of p and r at time t, which can be expressed in terms of the coupled Langevin equations for p and r that represent the two-stage model in continuum form. The steady-state distribution of p alone, P(p), is obtained from P(p, r, t) (a bivariate Gaussian) by integrating out the r degrees of freedom and taking the limit t -> infinity. P(p) is found to be proportional to the product of a Gaussian and a complementary error function. It provides a generally satisfactory fit to simulation data on the same two-stage process when the translational efficiency (a measure of intrinsic noise levels in the system) is relatively low; it is less successful as a model of the data when the translational efficiency (and noise levels) are high.

Free turbulent shear layer in a point vortex gas as a problem in nonequilibrium statistical mechanics

This paper attempts to unravel any relations that may exist between turbulent shear flows and statistical mechanics through a detailed numerical investigation in the simplest case where both can be well defined. The flow considered for the purpose is the two-dimensional (2D) temporal free shear layer with a velocity difference Delta U across it, statistically homogeneous in the streamwise direction (x) and evolving from a plane vortex sheet in the direction normal to it (y) in a periodic-in-x domain L x +/-infinity. Extensive computer simulations of the flow are carried out through appropriate initial-value problems for a ``vortex gas'' comprising N point vortices of the same strength (gamma = L Delta U/N) and sign. Such a vortex gas is known to provide weak solutions of the Euler equation. More than ten different initial-condition classes are investigated using simulations involving up to 32 000 vortices, with ensemble averages evaluated over up to 10(3) realizations and integration over 10(4)L/Delta U. The temporal evolution of such a system is found to exhibit three distinct regimes. In Regime I the evolution is strongly influenced by the initial condition, sometimes lasting a significant fraction of L/Delta U. Regime III is a long-time domain-dependent evolution towards a statistically stationary state, via ``violent'' and ``slow'' relaxations P.-H. Chavanis, Physica A 391, 3657 (2012)], over flow time scales of order 10(2) and 10(4)L/Delta U, respectively (for N = 400). The final state involves a single structure that stochastically samples the domain, possibly constituting a ``relative equilibrium.'' The vortex distribution within the structure follows a nonisotropic truncated form of the Lundgren-Pointin (L-P) equilibrium distribution (with negatively high temperatures; L-P parameter lambda close to -1). The central finding is that, in the intermediate Regime II, the spreading rate of the layer is universal over the wide range of cases considered here. The value (in terms of momentum thickness) is 0.0166 +/- 0.0002 times Delta U. Regime II, extensively studied in the turbulent shear flow literature as a self-similar ``equilibrium'' state, is, however, a part of the rapid nonequilibrium evolution of the vortex-gas system, which we term ``explosive'' as it lasts less than one L/Delta U. Regime II also exhibits significant values of N-independent two-vortex correlations, indicating that current kinetic theories that neglect correlations or consider them as O(1/N) cannot describe this regime. The evolution of the layer thickness in present simulations in Regimes I and II agree with the experimental observations of spatially evolving (3D Navier-Stokes) shear layers. Further, the vorticity-stream-function relations in Regime III are close to those computed in 2D Navier-Stokes temporal shear layers J. Sommeria, C. Staquet, and R. Robert, J. Fluid Mech. 233, 661 (1991)]. These findings suggest the dominance of what may be called the Kelvin-Biot-Savart mechanism in determining the growth of the free shear layer through large-scale momentum and vorticity dispersal.

Healing time for the growth of thin films on patterned substrates

The healing times for the growth of thin films on patterned substrates are studied using simulations of two discrete models of surface growth: the Family model and the Das Sarma-Tamborenea (DT) model. The healing time, defined as the time at which the characteristics of the growing interface are ``healed'' to those obtained in growth on a flat substrate, is determined via the study of the nearest-neighbor height difference correlation function. Two different initial patterns are considered in this work: a relatively smooth tent-shaped triangular substrate and an atomically rough substrate with singlesite pillars or grooves. We find that the healing time of the Family and DT models on aL x L triangular substrate is proportional to L-z, where z is the dynamical exponent of the models. For the Family model, we also analyze theoretically, using a continuum description based on the linear Edwards-Wilkinson equation, the time evolution of the nearest-neighbor height difference correlation function in this system. The correlation functions obtained from continuum theory and simulation are found to be consistent with each other for the relatively smooth triangular substrate. For substrates with periodic and random distributions of pillars or grooves of varying size, the healing time is found to increase linearly with the height (depth) of pillars (grooves). We show explicitly that the simulation data for the Family model grown on a substrate with pillars or grooves do not agree with results of a calculation based on the continuum Edwards-Wilkinson equation. This result implies that a continuum description does not work when the initial pattern is atomically rough. The observed dependence of the healing time on the substrate size and the initial height (depth) of pillars (grooves) can be understood from the details of the diffusion rule of the atomistic model. The healing time of both models for pillars is larger than that for grooves with depth equal to the height of the pillars. The calculated healing time for both Family and DT models is found to depend on how the pillars and grooves are distributed over the substrate. (C) 2014 Elsevier B.V. All rights reserved.

Role of Crosslinking and Entanglements in the Mechanics of Silicone Networks

The goal of this study is to investigate the applicability of different constitutive models for silicone networks using comprehensive multiaxial experimental tests, including non-equibiaxial mechanical tests which introduce differential constraints on the networks in the two orthogonal directions, on samples prepared using various crosslinking densities. Uniaxial stress-strain experiments show that a decrease in crosslinker amounts used in the preparation of silicone networks lead to more compliant material response as compared to that obtained using higher amounts of crosslinker. Biaxial data were used to obtain fits to the neo- Hookean, Mooney-Rivlin, Arruda-Boyce and the Edward-Vilgis slip-link constitutive models. Our results show that the slip-link model, based on separation of the individual contributions of chemical crosslinks and physical entanglements, is better at describing the stress-strain response of highly crosslinked networks at low stretches as compared to other constitutive models. Modulus obtained using the slip-link model for highly crosslinked networks agrees with experimentally determined values obtained using uniaxial tension experiments. In contrast, moduli obtained using coefficients to the other constitutive models underpredict experimentally determined moduli by over 40 %. However, the slip-link model did not predict the experimentally observed stiffening response at higher stretches which was better captured using the Arruda-Boyce model.

Fracture in metallic glasses: mechanics and mechanisms

Significant progress in understanding the mechanical behavior of metallic glasses (MGs) was made over the past decade, particularly on mechanisms of plastic deformation. However, recent research thrust has been on exploring the mechanics and physics of fracture. MGs can be very brittle with K-Ic values similar to silicate glasses and ceramics or very tough with K-Ic akin to high toughness crystalline metals. Even the tough MGs can become brittle with structural relaxation following annealing at temperatures close to glass transition temperature (T-g). Detailed experimental studies coupled with complementary numerical simulations of the recent past have provided insights on the micromechanisms of failure as well as nature of crack tip fields, and established the governing fracture criteria for ductile and brittle glasses. In this paper, the above advances are reviewed and outstanding issues in the context of fracture of amorphous alloys that need to be resolved are identified.

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

Estimation of mechanical properties of single wall carbon nanotubes using molecular mechanics approach

Molecular mechanics based finite element analysis is adopted in the current work to evaluate the mechanical properties of Zigzag, Armchair and Chiral Single wall Carbon Nanotubes (SWCNT) of different diameters and chiralities. Three different types of atomic bonds, that is Carbon Carbon covalent bond and two types of Carbon Carbon van der Waals bonds are considered in the carbon nanotube system. The stiffness values of these bonds are calculated using the molecular potentials, namely Morse potential function and Lennard-Jones interaction potential function respectively and these stiffness's are assigned to spring elements in the finite element model of the CNT. The geometry of CNT is built using a macro that is developed for the finite element analysis software. The finite element model of the CNT is constructed, appropriate boundary conditions are applied and the behavior of mechanical properties of CNT is studied.

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.

Handling large variations in mechanics: Some applications

There is a need to use probability distributions with power-law decaying tails to describe the large variations exhibited by some of the physical phenomena. The Weierstrass Random Walk (WRW) shows promise for modeling such phenomena. The theory of anomalous diffusion is now well established. It has found number of applications in Physics, Chemistry and Biology. However, its applications are limited in structural mechanics in general, and structural engineering in particular. The aim of this paper is to present some mathematical preliminaries related to WRW that would help in possible applications. In the limiting case, it represents a diffusion process whose evolution is governed by a fractional partial differential equation. Three applications of superdiffusion processes in mechanics, illustrating their effectiveness in handling large variations, are presented.

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.

A dilation-driven vortex flow in sheared granular materials explains a rheometric anomaly

Granular flows occur widely in nature and industry, yet a continuum description that captures their important features is yet not at hand. Recent experiments on granular materials sheared in a cylindrical Couette device revealed a puzzling anomaly, wherein all components of the stress rise nearly exponentially with depth. Here we show, using particle dynamics simulations and imaging experiments, that the stress anomaly arises from a remarkable vortex flow. For the entire range of fill heights explored, we observe a single toroidal vortex that spans the entire Couette cell and whose sense is opposite to the uppermost Taylor vortex in a fluid. We show that the vortex is driven by a combination of shear-induced dilation, a phenomenon that has no analogue in fluids, and gravity flow. Dilatancy is an important feature of granular mechanics, but not adequately incorporated in existing models.