Crack propagation in graphene

The crack initiation and growth mechanisms in an 2D graphene lattice structure are studied based on molecular dynamics simulations. Crack growth in an initial edge crack model in the arm-chair and the zig-zag lattice configurations of graphene are considered. Influence of the time steps on the post yielding behaviour of graphene is studied. Based on the results, a time step of 0.1 fs is recommended for consistent and accurate simulation of crack propagation. Effect of temperature on the crack propagation in graphene is also studied, considering adiabatic and isothermal conditions. Total energy and stress fields are analyzed. A systematic study of the bond stretching and bond reorientation phenomena is performed, which shows that the crack propagates after significant bond elongation and rotation in graphene. Variation of the crack speed with the change in crack length is estimated. (C) 2015 AIP Publishing LLC.

Kinematic flow patterns in slow deformation of a dense granular material

The kinematic flow pattern in slow deformation of a model dense granular medium is studied at high resolution using in situ imaging, coupled with particle tracking. The deformation configuration is indentation by a flat punch under macroscopic plane-strain conditions. Using a general analysis method, velocity gradients and deformation fields are obtained from the disordered grain arrangement, enabling flow characteristics to be quantified. The key observations are the formation of a stagnation zone, as in dilute granular flow past obstacles; occurrence of vortices in the flow immediately underneath the punch; and formation of distinct shear bands adjoining the stagnation zone. The transient and steady state stagnation zone geometry, as well as the strength of the vortices and strain rates in the shear bands, are obtained from the experimental data. All of these results are well-reproduced in exact-scale non-smooth contact dynamics simulations. Full 3D numerical particle positions from the simulations allow extraction of flow features that are extremely difficult to obtain from experiments. Three examples of these, namely material free surface evolution, deformation of a grain column below the punch and resolution of velocities inside the primary shear band, are highlighted. The variety of flow features observed in this model problem also illustrates the difficulty involved in formulating a complete micromechanical analytical description of the deformation.

An automatic algorithm for mixed mode crack growth rate based on drop potential method

A new automatic algorithm for the assessment of mixed mode crack growth rate characteristics is presented based on the concept of an equivalent crack. The residual ligament size approach is introduced to implementation this algorithm for identifying the crack tip position on a curved path with respect to the drop potential signal. The automatic algorithm accounting for the curvilinear crack trajectory and employing an electrical potential difference was calibrated with respect to the optical measurements for the growing crack under cyclic mixed mode loading conditions. The effectiveness of the proposed algorithm is confirmed by fatigue tests performed on ST3 steel compact tension-shear specimens in the full range of mode mixities from pure mode Ito pure mode II. (C) 2015 Elsevier Ltd. All rights reserved.

Seismotectonics of the April-May 2015 Nepal earthquakes: An assessment based on the aftershock patterns, surface effects and deformational characteristics

Occurrence of the April 25, 2015 (Mw 7.8) earthquake near Gorkha, central Nepal, and another one that followed on May 12 (Mw 7.3), located similar to 140 km to its east, provides an exceptional opportunity to understand some new facets of Himalayan earthquakes. Here we attempt to assess the seismotectonics of these earthquakes based on the deformational field generated by these events, along with the spatial and temporal characteristics of their aftershocks. When integrated with some of the post-earthquake field observations, including the localization of damage and surface deformation, it became obvious that although the mainshock slip was mostly limited to the Main Himalayan Thrust (MHT), the rupture did not propagate to the Main Frontal Thrust (MFT). Field evidence, supported by the available InSAR imagery of the deformation field, suggests that a component of slip could have emerged through a previously identified out-of-sequence thrust/active thrust in the region that parallels the Main Central Thrust (MCT), known in the literature as a co-linear physiographic transitional zone called PT2. Termination of the first rupture, triggering of the second large earthquake, and distribution of aftershocks are also spatially constrained by the eastern extremity of PT2. Mechanism of the 2015 sequence demonstrates that the out-of-sequence thrusts may accommodate part of the slip, an aspect that needs to be considered in the current understanding of the mechanism of earthquakes originating on the MHT. (c) 2015 Elsevier Ltd. All rights reserved.

Divergence in Patterns of Leaf Growth Polarity Is Associated with the Expression Divergence of miR396

Lateral appendages often show allometric growth with a specific growth polarity along the proximo-distal axis. Studies on leaf growth in model plants have identified a basipetal growth direction with the highest growth rate at the proximal end and progressively lower rates toward the distal end. Although the molecular mechanisms governing such a growth pattern have been studied recently, variation in leaf growth polarity and, therefore, its evolutionary origin remain unknown. By surveying 75 eudicot species, here we report that leaf growth polarity is divergent. Leaf growth in the proximo-distal axis is polar, with more growth arising from either the proximal or the distal end; dispersed with no apparent polarity; or bidirectional, with more growth contributed by the central region and less growth at either end. We further demonstrate that the expression gradient of the miR396-GROWTH-REGULATING FACTOR module strongly correlates with the polarity of leaf growth. Altering the endogenous pattern of miR396 expression in transgenic Arabidopsis thaliana leaves only partially modified the spatial pattern of cell expansion, suggesting that the diverse growth polarities might have evolved via concerted changes in multiple gene regulatory networks.

Complex patterns arise through spontaneous symmetry breaking in dense homogeneous networks of neural oscillators

There has been much interest in understanding collective dynamics in networks of brain regions due to their role in behavior and cognitive function. Here we show that a simple, homogeneous system of densely connected oscillators, representing the aggregate activity of local brain regions, can exhibit a rich variety of dynamical patterns emerging via spontaneous breaking of permutation or translational symmetries. Upon removing just a few connections, we observe a striking departure from the mean-field limit in terms of the collective dynamics, which implies that the sparsity of these networks may have very important consequences. Our results suggest that the origins of some of the complicated activity patterns seen in the brain may be understood even with simple connection topologies.

Investigation Into Harmful Patterns Over Multitrack Shingled Magnetic Detection Using the Voronoi Model

Two-dimensional magnetic recording 2-D (TDMR) is a promising technology for next generation magnetic storage systems based on a systems-level framework involving sophisticated signal processing at the core. The TDMR channel suffers from severe jitter noise along with electronic noise that needs to be mitigated during signal detection and recovery. Recently, we developed noise prediction-based techniques coupled with advanced signal detectors to work with these systems. However, it is important to understand the role of harmful patterns that can be avoided during the encoding process. In this paper, we investigate the Voronoi-based media model to study the harmful patterns over multitrack shingled recording systems. Through realistic quasi-micromagnetic simulation studies, we identify 2-D data patterns that contribute to high media noise. We look into the generic Voronoi model and present our analysis on multitrack detection with constrained coded data. We show that the 2-D constraints imposed on input patterns result in an order of magnitude improvement in the bit-error rate for the TDMR systems. The use of constrained codes can reduce the complexity of 2-D intersymbol interference (ISI) signal detection, since the lesser 2-D ISI span can be accommodated at the cost of a nominal code rate loss. However, a system must be designed carefully so that the rate loss incurred by a 2-D constraint does not offset the detector performance gain due to more distinguishable readback signals.

Complex potentials and singular integral equation for curve crack problem in antiplane elasticity

Four types of the fundamental complex potential in antiplane elasticity are introduced: (a) a point dislocation, (b) a concentrated force, (c) a dislocation doublet and (d) a concentrated force doublet. It is proven that if the axis of the concentrated force doublet is perpendicular to the direction of the dislocation doublet, the relevant complex potentials are equivalent. Using the obtained complex potentials, a singular integral equation for the curve crack problem is introduced. Some particular features of the obtained singular integral equation are discussed, and numerical solutions and examples are given.

Monte-Carlo simulation of surface crack growth rate for offshore structural steel E36-Z35

The Monte- Carlo method is used to simulate the surface fatigue crack growth rate for offshore structural steel E36-Z35, and to determine the distributions and relevance of the parameters in the Paris equation. By this method, the time and cost of fatigue crack propagation testing can be reduced. The application of the method is demonstrated by use of four sets of fatigue crack propagation data for offshore structural steel E36-Z35. A comparison of the test data with the theoretical prediction for surface crack growth rate shows the application of the simulation method to the fatigue crack propagation tests is successful.

Analysis of strip electric saturation model of crack problem in piezoelectric materials

This paper presents a fully anisotropic analysis of strip electric saturation model proposed by Gao et al. (1997) (Gao, H.J., Zhang, T.Y., Tong, P., 1997. Local and global energy release rates for an electrically yielded crack in a piezoelectric ceramic. J. Mech. Phys. Solids, 45, 491-510) for piezoelectric materials. The relationship between the size of the strip saturation zone ahead of a crack tip and the applied electric displacement field is established. It is revealed that the critical fracture stresses for a crack perpendicular to the poling axis is linearly decreased with the increase of the positive applied electric field and increases linearly with the increase of the negative applied electric field. For a crack parallel to the poring axis, the failure stress is not effected by the parallel applied electric field. In order to analyse the existed experimental results, the stress fields ahead of the tip of an elliptic notch in an infinite piezoelectric solid are calculated. The critical maximum stress criterion is adopted for determining the fracture stresses under different remote electric displacement fields. The present analysis indicates that the crack initiation and propagation from the tip of a sharp elliptic notch could be aided or impeded by an electric displacement field depending on the field direction. The fracture stress predicted by the present analysis is consistent with the experimental data given by Park and Sun (1995) (Park, S., Sun, C.T., 1995. Fracture criteria for piezoelectric materials. J. Am. Ceram. Soc 78, 1475-1480).

Effect of temperature on dislocation emission from crack tip for BCC metal Mo

The effect of thermally activated energy on the dislocation emission from a crack tip in BCC metal Mo is simulated in this paper. Based on the correlative reference model on which the flexible displacement boundary scheme is introduced naturally, the simulation shows that as temperature increases the critical stress intensity factor for the first dislocation emission will decrease and the total number of emitted dislocations increase for the same external load. The dislocation velocity and extensive distance among partial dislocations are not sensitive to temperature. After a dislocation emission, two different deformation slates are observed, the stable and unstable deformation states. In the stable deformation slate, the nucleated dislocation will emit from the crack tip and piles up at a distance far away from the crack tip, after that the new dislocation can not be nucleated unless the external loading increases. In the unstable deformation state, a number of dislocations can be emitted from the crack lip continuously under the same external load.

Thermal-mechanical interface crack behaviour of a surface mount solder joint

The effect of thermal-mechanical loading on a surface mount assembly with interface cracks between the solder and the resistor and between the solder and the printed circuit board (PCB) was studied using a non-linear thermal finite element analysis. The thermal effect was taken as cooling from the solder eutectic temperature to room temperature. Mechanical loading at the ends of the PCB was also applied. The results showed that cooling had the effect of causing large residual shear displacement at the region near the interface cracks. The mechanical loading caused additional crack opening displacements. The analysis on the values of J-integral for the interface cracks showed that J-integral was approximately path independent, and that the effect of crack at the solder/PCB interface is much more serious than that between the component and solder.

An analysis on overall crack-number-density of short-fatigue-cracks

The evolution of dispersed short-fatigue-cracks is analysed based on the equilibrium of crack-number-density (CND). By separating the mean value and the stochastic fluctuation of local CND, the equilibrium equation of overall CND is derived. Comparing with the mean-field equilibrium equation, the equilibrium equation of overall CND has different forms in the expression of crack-nucleation-rate or crack-growth-rate. The simulation results are compared with experimental measurements showing the stochastic analyses provide consistent tendency with experiments. The discrepancy in simulation results between overall CND and mean-field CND is discussed.

Dynamic behavior of a cylindrical crack in a functionally graded interlayer under torsional loading

In this paper the problem of a cylindrical crack located in a functionally graded material (FGM) interlayer between two coaxial elastic dissimilar homogeneous cylinders and subjected to a torsional impact loading is considered. The shear modulus and the mass density of the FGM interlayer are assumed to vary continuously between those of the two coaxial cylinders. This mixed boundary value problem is first reduced to a singular integral equation with a Cauchy type kernel in the Laplace domain by applying Laplace and Fourier integral transforms. The singular integral equation is then solved numerically and the dynamic stress intensity factor (DSIF) is also obtained by a numerical Laplace inversion technique. The DSIF is found to rise rapidly to a peak and then reduce and tend to the static value almost without oscillation. The influences of the crack location, the FGM interlayer thickness and the relative magnitudes of the adjoining material properties are examined. It is found among others that, by increasing the FGM gradient, the DSIF can be greatly reduced.