Contact instabilities of anisotropic and inhomogeneous soft elastic films

Anisotropy plays important roles in various biological phenomena such as adhesion of geckos and grasshoppers enabled by the attachment pods having hierarchical structures like thin longitudinal setae connected with threads mimicked by anisotropic films. We study the contact instability of a transversely isotropic thin elastic film when it comes in contact proximity of another surface. In the present study we investigate the contact stability of a thin incompressible transversely isotropic film by performing linear stability analysis. Based on the linear stability analysis, we show that an approaching contactor renders the film unstable. The critical wavelength of the instability is a function of the total film thickness and the ratio of the Young's modulus in the longitudinal direction and the shear modulus in the plane containing the longitudinal axis. We also analyze the stability of a thin gradient film that is elastically inhomogeneous across its thickness. Compared to a homogeneous elastic film, it becomes unstable with a longer wavelength when the film becomes softer in going from the surface to the substrate.

A rational approach to model the behavior of bonded soils

The paper presents a rational approach to model the behavior of bonded soils within the frame work of hardening plasticity. The approach is based on the premise that the resistance of bonded materials is a superposition of the two components of cement bond strength and soil frictional strength and that the deformation of the soil is associated with the frictional component of stresses just as in the case of a remoulded soil, the bonds offering additional resistance at any given strain level. This concept is similar to two stiffnesses acting in parallel for the same strain response. The proposed model considers the constitutive laws separately for the two components (bond and frictional) and adds the two to get the overall response. The unbonded soil component is described by the well known 'modified Cam clay' model. The response of the bond component is also described by a strain softening elasto-plastic model, considering the behavior to be elastic up to the yield surface and elasto-plastic beyond yield surface. To illustrate the capability of the proposed, model some laboratory test results of both compression and-extension shear tests are predicted. Despite the model being simple, several typical features of the behavior of bonded materials are well reproduced. The model parameters are well defined and easily determinable.

Effect of temperature on the plastic zone size and the shear band density in a bulk metallic glass

High temperature bonded interface indentation experiments are carried out on a Zr based bulk metallic glass (BMG) to examine the plastic deformation characteristics in subsurface deformation zone under a Vickers indenter. The results show that the shear bands are semi-circular in shape and propagate in radial direction. At all temperatures the inter-band spacing along the indentation axis is found to increase with increasing distance from the indenter tip. The average shear band spacing monotonically increases with temperature whereas the shear band induced plastic deformation zone is invariant with temperature. These observations are able to explain the increase in pressure sensitive plastic flow of BMGs with temperature. (C) 2011 Elsevier B.V. All rights reserved.

Effect of strain rate and temperature on the plastic deformation behaviour of a bulk metallic glass composite

The composites consisting of amorphous matrix reinforced with crystalline dendrites offer extraordinary combinations of strength, stiffness, and toughness and can be processed in bulk. Hence, they have been receiving intense research interest, with a primary focus to study their mechanical properties. In this paper, the temperature and strain rate effects on the uniaxial compression response of a tailored bulk metallic glass (BMG) composite has been investigated. Experimental results show that at temperatures ranging between ambient to 500 K and at all strain rates; the onset of plastic deformation in the composite is controlled by that in the dendrites. As the temperature is increased to the glass transition temperature of the matrix and beyond, flow in the amorphous matrix occurs readily and hence it dictates the composite's response. The role of the constituent phases in controlling the deformation mechanism of the composite has been verified by assessing the strain rate sensitivity and the activation volume for deformation. The composite is rate sensitive at room temperature with values of strain rate sensitivity and activation volume being similar to that of the dendrites. At test temperatures near to the glass transition temperature, the composite however becomes rate-insensitive corresponding to that of the matrix phase. At low strain rates, serrated flow akin to that of dynamic strain ageing in crystalline alloys was observed and the serration magnitude decreases with increasing temperature. Initiation of the shear bands at the dendrite/matrix interface and propagation of them through the matrix ligaments until their arrest at another interface is the responsible mechanism for this. (C) 2011 Elsevier B.V. All rights reserved.

Quasi-static crack tip fields in rate-sensitive FCC single crystals

In this work, the effects of loading rate, material rate sensitivity and constraint level on quasi-static crack tip fields in a FCC single crystal are studied. Finite element simulations are performed within a mode I, plane strain modified boundary layer framework by prescribing the two term (K-T) elastic crack tip field as remote boundary conditions. The material is assumed to obey a rate-dependent crystal plasticity theory. The orientation of the single crystal is chosen so that the crack surface coincides with the crystallographic (010) plane and the crack front lies along 101] direction. Solutions corresponding to different stress intensity rates K., T-stress values and strain rate exponents m are obtained. The results show that the stress levels ahead of the crack tip increase with K. which is accompanied by gradual shrinking of the plastic zone size. However, the nature of the shear band patterns around the crack tip is not affected by the loading rate. Further, it is found that while positive T-stress enhances the opening and hydrostatic stress levels ahead of crack tip, they are considerably reduced with imposition of negative T-stress. Also, negative T-stress promotes formation of shear bands in the forward sector ahead of the crack tip and suppresses them behind the tip.

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.

A numerical approach to determine the sufficiency of given boundary data sets for uniquely estimating interior elastic properties

We consider an inverse elasticity problem in which forces and displacements are known on the boundary and the material property distribution inside the body is to be found. In other words, we need to estimate the distribution of constitutive properties using the finite boundary data sets. Uniqueness of the solution to this problem is proved in the literature only under certain assumptions for a given complete Dirichlet-to-Neumann map. Another complication in the numerical solution of this problem is that the number of boundary data sets needed to establish uniqueness is not known even under the restricted cases where uniqueness is proved theoretically. In this paper, we present a numerical technique that can assess the sufficiency of given boundary data sets by computing the rank of a sensitivity matrix that arises in the Gauss-Newton method used to solve the problem. Numerical experiments are presented to illustrate the method.

Projecting low-dimensional chaos from spatiotemporal dynamics in a model for plastic instability

We investigate the possibility of projecting low-dimensional chaos from spatiotemporal dynamics of a model for a kind of plastic instability observed under constant strain rate deformation conditions. We first discuss the relationship between the spatiotemporal patterns of the model reflected in the nature of dislocation bands and the nature of stress serrations. We show that at low applied strain rates, there is a one-to-one correspondence with the randomly nucleated isolated bursts of mobile dislocation density and the stress drops. We then show that the model equations are spatiotemporally chaotic by demonstrating the number of positive Lyapunov exponents and Lyapunov dimension scale with the system size at low and high strain rates. Using a modified algorithm for calculating correlation dimension density, we show that the stress-strain signals at low applied strain rates corresponding to spatially uncorrelated dislocation bands exhibit features of low-dimensional chaos. This is made quantitative by demonstrating that the model equations can be approximately reduced to space-independent model equations for the average dislocation densities, which is known to be low-dimensionally chaotic. However, the scaling regime for the correlation dimension shrinks with increasing applied strain rate due to increasing propensity for propagation of the dislocation bands.

SCALE EFFECTS ON ULTRASONIC WAVE DISPERSION CHARACTERISTICS OF MONOLAYER GRAPHENE EMBEDDED IN AN ELASTIC MEDIUM

Ultrasonic wave propagation in a graphene sheet, which is embedded in an elastic medium, is studied using nonlocal elasticity theory incorporating small-scale effects. The graphene sheet is modeled as an one-atom thick isotropic plate and the elastic medium/substrate is modeled as distributed 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. After that, an ultrasonic type of wave propagation model is also derived. The explicit expressions for the cut-off frequencies are also obtained as functions of the nonlocal scaling parameter and the y-directional wavenumber. Local elasticity shows that the wave will propagate even at higher frequencies. But nonlocal elasticity predicts that the waves can propagate only up to certain frequencies (called escape frequencies), after which the wave velocity becomes zero. The results also show that the escape frequencies are purely a function of the nonlocal scaling parameter. The effect of the elastic medium is captured in the wave dispersion analysis and this analysis is explained with respect to both local and nonlocal elasticity. The simulations show that the elastic medium affects only the flexural wave mode in the graphene sheet. The presence of the elastic matrix increases the band gap of the flexural mode. The present results can provide useful guidance for the design of next-generation nanodevices in which graphene-based composites act as a major element.

Time-Temperature Scaling of Conductivity Spectra of Organic Plastic Crystalline Conductors

Organic plastic crystalline soft matter ion conductors are interesting alternatives to liquid electrolytes in electrochemical storage devices such as Lithium-ion batteries. The solvent dynamics plays a major role in determining the ion transport in plastic crystalline ion conductors. We present here an analysis of the frequency-dependent ionic conductivity of succinonitrile-based plastic crystalline ion conductors at varying salt composition (0.005 to 1 M) and temperature (-20 to 60 degrees C) using time-temperature superposition principle (TTSP). The main motivation of the work has been to establish comprehensive insight into the ion transport mechanism from a single method viz, impedance spectroscopy rather than employing cluster of different characterization methods probing various length and time scales. The TTSP remarkably aids in explicit identification of the extent of the roles of solvent dynamics and ion-ion interactions on the effective conductivity of the orientationally disordered plastic crystalline ion conductors.

Ductile fracture by multiple void growth and interaction ahead of a notch tip in polycrystalline plastic solids

The objectives of this paper are to study the effects of plastic anisotropy and evolution in crystallographic texture with deformation on the ductile fracture behaviour of polycrystalline solids. To this end, numerical simulations of multiple void growth and interaction ahead of a notch tip are performed under mode I, plane strain, small scale yielding conditions using two approaches. The first approach is based on the Hill yield theory, while the second employs crystal plasticity constitutive equations and a Taylor-type homogenization in order to represent the ductile polycrystalline solid. The initial textures pertaining to continuous cast Al-Mg AA5754 sheets in recrystallized and cold rolled conditions are considered. The former is nearly-isotropic, while the latter displays pronounced anisotropy. The results indicate distinct changes in texture in the ligaments bridging the voids ahead of the notch tip with increase in load level which gives rise to retardation in porosity evolution and increase in tearing resistance for both materials.

Conditions for the onset of elastic and material instabilities in hyperelastic materials

Several constitutive inequalities have been proposed in the literature to quantify the notion that ‘stress increases with strain’ in an elastic material. Due to some inherent shortcomings in them, which we discuss, we propose a new tensorial criterion for isotropic materials. We also present necessary conditions in terms of elasticity tensors for the onset of elastic instabilities.

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.

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.

Effects of different modes of hot cross-rolling in 7010 aluminum alloy: part II. mechanical properties anisotropy

The influence of microstructure and texture developed by different modes of hot cross-rolling on in-plane anisotropy (A (IP)) of yield strength, work hardening behavior, and anisotropy of Knoop hardness (KHN) yield locus has been investigated. The A (IP) and work hardening behavior are evaluated by tensile testing at 0 deg, 45 deg, and 90 deg to the rolling direction, while yield loci have been generated by directional KHN measurements. It has been observed that specimens especially in the peak-aged temper, in spite of having a strong, rotated Brass texture, show low A (IP). The results are discussed on the basis of Schmid factor analyses in conjunction with microstructural features, namely grain morphology and precipitation effects. For the specimen having a single-component texture, the yield strength variation as a function of orientation can be rationalized by the Schmid factor analysis of a perfectly textured material behaving as a quasi-single crystal. The work hardening behavior is significantly affected by the presence of solute in the matrix and the state of precipitation rather than texture, while yield loci derived from KHN measurements reiterate the low anisotropy of the materials. Theoretic yield loci calculated from the texture data using the visco-plastic self-consistent model and Hill's anisotropic equation are compared with that obtained experimentally.