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.

Estimation of the shear transformation zone size fin a bulk metallic glass through statistical a analysis of the first pop-in stresses during spherical nanoindentation

The size of the shear transformation zone (STZ) that initiates the elastic to plastic transition in a Zr-based bulk metallic glass was estimated by conducting a statistical analysis of the first pop-in event during spherical nanoindentation. A series of experiments led us to a successful description of the distribution of shear strength for the transition and its dependence on the loading rate. From the activation volume determined by statistical analysis the STZ size was estimated based on a cooperative shearing model. (C) 2012 Acta Materialia Inc. Published by Elsevier Ltd. All rights reserved.

Dynamics of edge dislocations in a sheared lamellar mesophase

The dynamics and interactions of edge dislocations in a nearly aligned sheared lamellar mesophase is analysed to provide insights into the relationship between disorder and rheology. First, the mesoscale permeation and momentum equations for the displacement field in the presence of external forces are derived from the model H equations for the concentration and momentum field. The secondary flow generated due to the mean shear around an isolated defect is calculated, and the excess viscosity due to the presence of the defect is determined from the excess energy dissipation due to the secondary flow. The excess viscosity for an isolated defect is found to increase with system size in the cross-stream direction as L-3/2 for an isolated defect, though this divergence is cut-off due to interactions in a defect suspension. As the defects are sheared past each other due to the mean flow, the Peach-Koehler force due to elastic interaction between pairs of defects is found to cause no net displacement relative to each other as they approach from large separation to the distance of closest approach. The equivalent force due to viscous interactions is found to increase the separation for defects of opposite sign, and decrease the separation for defects of same sign. During defect interactions, we find that there is no buckling instability due to dilation of layers for systems of realistic size. However, there is another mechanism, which is the velocity difference generated across a slightly deformed bilayer due to the mean shear, which could result in the creation of new defects. (C) 2013 AIP Publishing LLC.

A parametric study on the buckling of functionally graded material plates with internal discontinuities using the partition of unity method

In this paper, the effect of local defects, viz., cracks and cutouts on the buckling behaviour of functionally graded material plates subjected to mechanical and thermal load is numerically studied. The internal discontinuities, viz., cracks and cutouts are represented independent of the mesh within the framework of the extended finite element method and an enriched shear flexible 4-noded quadrilateral element is used for the spatial discretization. The properties are assumed to vary only in the thickness direction and the effective properties are estimated using the Mori-Tanaka homogenization scheme. The plate kinematics is based on the first order shear deformation theory. The influence of various parameters, viz., the crack length and its location, the cutout radius and its position, the plate aspect ratio and the plate thickness on the critical buckling load is studied. The effect of various boundary conditions is also studied. The numerical results obtained reveal that the critical buckling load decreases with increase in the crack length, the cutout radius and the material gradient index. This is attributed to the degradation in the stiffness either due to the presence of local defects or due to the change in the material composition. (C) 2013 Elsevier Masson SAS. All rights reserved.

Chirality dependent elastic properties of single-walled boron nitride nanotubes under uniaxial and torsional loading

The elastic behavior of single-walled boron nitride nanotubes is studied under axial and torsional loading. Molecular dynamics simulation is carried out with a tersoff potential for modeling the interatomic interactions. Different chiral configurations with similar diameter are considered to study the effect of chirality on the elastic and shear moduli. Furthermore, the effects of tube length on elastic modulus are also studied by considering different aspects ratios. It is observed that both elastic and shear moduli depend upon the chirality of a nanotube. For aspect ratios less than 15, the elastic modulus reduces monotonically with an increase in the chiral angle. For chiral nanotubes, the torsional response shows a dependence on the direction of loading. The difference between the shear moduli against and along the chiral twist directions is maximum for chiral angle of 15 degrees, and zero for zigzag (0 degrees) and armchair (30 degrees) configurations. (C) 2014 AIP Publishing LLC.

Tectonic inheritance of the Indian Shield: New insights from its elastic thickness structure

A new evaluation of the elastic thickness (Te) structure of the Indian Shield, derived from isotropic fan wavelet methodology, documents spatial variations of lithospheric deformation in different tectonic provinces correlated with episodic tectono-thermal events. The Te variations corroborated by shear velocity, crustal thickness, and seismogenic thickness reveal the heterogeneous rheology of the Indian lithosphere. The thinned, attenuated lithosphere beneath Peninsular India is considered to be the reason for its mechanically weak strength (<30 km), where a decoupled crust-mantle rheology under different surface/subsurface loading structures may explain the prominent low Te patterns. The arcuate Te structure of the Western Dharwar province and a NNE-trending band of low Te anomaly in the Southern Granulite Terrane are intriguing patterns. The average Te values (40-50 km) of the Central Indian Tectonic Zone, the Bastar Craton, and the northern Eastern Ghats Mobile Belt are suggestive of old, stable, Indian lithosphere, which was not affected by any major tectono-thermal events after cratonic stabilization. We propose that the anomalously high Te (60-85 km) and high S-wave velocity zone to the north of the Narmada-Son Lineament, mainly in NW Himalaya, and the northern Aravalli and Bundelkhand Cratons, suggest that Archean lithosphere characterized by a high velocity mantle keel supports the orogenic topographic loads in/near the Himalaya. The Te map clearly segments the volcanic provinces of the Indian Shield, where the signatures of the Reunion, Marion, and Kerguelen hotspots are indicated by significantly low Te patterns that correlate with plume- and rift-related thermal and mechanical rejuvenation, magmatic underplating, and crustal necking. The correlations between Te variations and the occurrence of seismicity over seismically active zones reveal different causal relationships, which led to the current seismogenic zonation of the Indian Shield. (C) 2013 Elsevier B.V. All rights reserved.

Transmission of elastic waves through a frictional boundary

Incident energy gets transmitted, reflected and absorbed across an interface in jointed rock mass leading to energy dissipation and alteration of waves. Wave velocities get attenuated during their propagation across joints and this behavior is studied using bender/extender element tests. The velocity attenuation and modulus reduction observed in experimental tests are modeled with three dimensional distinct element code and results are validated. Normal propagation of an incident shear wave through a jointed rock mass cause slip of the rock blocks if shear stress of wave exceeds the shear strength of the joint. As the properties of joint determine the transmission of energy across an interface, a parametric study is then conducted with the validated numerical model by varying the parameters that may determine the energy transmission across a joint using modified Miller's method. Results of the parametric study are analyzed and presented in the paper. (C) 2014 Elsevier Ltd. All rights reserved.

Wall-mode instability in plane shear flow of viscoelastic fluid over a deformable solid

The linear stability analysis of a plane Couette flow of an Oldroyd-B viscoelastic fluid past a flexible solid medium is carried out to investigate the role of polymer addition in the stability behavior. The system consists of a viscoelastic fluid layer of thickness R, density rho, viscosity eta, relaxation time lambda, and retardation time beta lambda flowing past a linear elastic solid medium of thickness HR, density rho, and shear modulus G. The emphasis is on the high-Reynolds-number wall-mode instability, which has recently been shown in experiments to destabilize the laminar flow of Newtonian fluids in soft-walled tubes and channels at a significantly lower Reynolds number than that for flows in rigid conduits. For Newtonian fluids, the linear stability studies have shown that the wall modes become unstable when flow Reynolds number exceeds a certain critical value Re c which scales as Sigma(3/4), where Reynolds number Re = rho VR/eta, V is the top-plate velocity, and dimensionless parameter Sigma = rho GR(2)/eta(2) characterizes the fluid-solid system. For high-Reynolds-number flow, the addition of polymer tends to decrease the critical Reynolds number in comparison to that for the Newtonian fluid, indicating a destabilizing role for fluid viscoelasticity. Numerical calculations show that the critical Reynolds number could be decreased by up to a factor of 10 by the addition of small amount of polymer. The critical Reynolds number follows the same scaling Re-c similar to Sigma(3/4) as the wall modes for a Newtonian fluid for very high Reynolds number. However, for moderate Reynolds number, there exists a narrow region in beta-H parametric space, corresponding to very dilute polymer solution (0.9 less than or similar to beta < 1) and thin solids (H less than or similar to 1.1), in which the addition of polymer tends to increase the critical Reynolds number in comparison to the Newtonian fluid. Thus, Reynolds number and polymer properties can be tailored to either increase or decrease the critical Reynolds number for unstable modes, thus providing an additional degree of control over the laminar-turbulent transition.

Effect of oxygen vacancies on the elastic properties of zinc oxide: A first-principles investigation

The effect of oxygen vacancies on the elastic properties of zinc oxide (ZnO) is examined using first-principles calculations based on density functional theory. Formation energies of vacancies in different types of oxygen deficient structures were analyzed to ascertain their stability. This analysis reveals that the doubly-charged oxygen vacancy under zinc-rich growth conditions is the most stable. Results show considerable degradation of some of the elastic moduli due to the presence of oxygen vacancies, which is in agreement with recent experiments. The decrease observed in elastic constants is more pronounced with increase in vacancy concentration. Further, the charge state of the defect structure was found to influence the shear elastic constants. Evaluation of elastic anisotropy of stoichiometric and oxygen deficient ZnO indicates the significant anisotropy in elastic properties and stiff c-axis orientation. (C) 2014 Elsevier B.V. All rights reserved.

Phenomenology and control of buckling dynamics in multicomponent colloidal droplets

Self-assembly of nano sized particles during natural drying causes agglomeration and shell formation at the surface of micron sized droplets. The shell undergoes sol-gel transition leading to buckling at the weakest point on the surface and produces different types of structures. Manipulation of the buckling rate with inclusion of surfactant (sodium dodecyl sulphate, SDS) and salt (anilinium hydrochloride, AHC) to the nano-sized particle dispersion (nanosilica) is reported here in an acoustically levitated single droplet. Buckling in levitated droplets is a cumulative, complicated function of acoustic streaming, chemistry, agglomeration rate, porosity, radius of curvature, and elastic energy of shell. We put forward our hypothesis on how buckling occurs and can be suppressed during natural drying of the droplets. Global precipitation of aggregates due to slow drying of surfactant-added droplets (no added salts) enhances the rigidity of the shell formed and hence reduces the buckling probability of the shell. On the contrary, adsorption of SDS aggregates on salt ions facilitates the buckling phenomenon with an addition of minute concentration of the aniline salt to the dispersion. Variation in the concentration of the added particles (SDS/AHC) also leads to starkly different morphologies and transient behaviour of buckling (buckling modes like paraboloid, ellipsoid, and buckling rates). Tuning of the buckling rate causes a transition in the final morphology from ring and bowl shapes to cocoon type of structure. (C) 2015 AIP Publishing LLC.

Stability analysis of thick piezoelectric metal based FGM plate using first order and higher order shear deformation theory

This paper presents the stability analysis of functionally graded plate integrated with piezoelectric actuator and sensor at the top and bottom face, subjected to electrical and mechanical loading. The finite element formulation is based on first order and higher order shear deformation theory, degenerated shell element, von-Karman hypothesis and piezoelectric effect. The equation for static analysis is derived by using the minimum energy principle and solutions for critical buckling load is obtained by solving eigenvalue problem. The material properties of the functionally graded plate are assumed to be graded along the thickness direction according to simple power law function. Two types of boundary conditions are used, such as SSSS (simply supported) and CSCS (simply supported along two opposite side perpendicular to the direction of compression and clamped along the other two sides). Sensor voltage is calculated using present analysis for various power law indices and FG (functionally graded) material gradations. The stability analysis of piezoelectric FG plate is carried out to present the effects of power law index, material variations, applied mechanical pressure and piezo effect on buckling and stability characteristics of FG plate.

In-plane effective shear modulus of generalized circular honeycomb structures and bundled tubes in a diamond array structure

Use of circular hexagonal honeycomb structures and tube assemblies in energy absorption systems has attracted a large number of literature on their characterization under crushing and impact loads. Notwithstanding these, effective shear moduli (G*) required for complete transverse elastic characterization and in analyses of hierarchical structures have received scant attention. In an attempt to fill this void, the present study undertakes to evaluate G* of a generalized circular honeycomb structures and tube assemblies in a diamond array structure (DAS) with no restriction on their thickness. These structures present a potential to realize a spectrum of moduli with minimal modifications, a point of relevance for manufactures and designers. To evaluate G* in this paper, models based on technical theories - thin ring theory and curved beam theory - and rigorous theory of elasticity are investigated and corroborated with FEA employing contact elements. Technical theories which give a good match for thin HCS offer compact expressions for moduli which can be harvested to study sensitivity of moduli on topology. On the other hand, elasticity model offers a very good match over a large range of thickness along with exact analysis of stresses by employing computationally efficient expressions. (C) 2015 Elsevier Ltd. All rights reserved.

Numerical study of driven flows of shear thinning viscoelastic fluids in rectangular cavities

In this work, we present a numerical study of flow of shear thinning viscoelastic fluids in rectangular lid driven cavities for a wide range of aspect ratios (depth to width ratio) varying from 1/16 to 4. In particular, the effect of elasticity, inertia, model parameters and polymer concentration on flow features in rectangular driven cavity has been studied for two shear thinning viscoelastic fluids, namely, Giesekus and linear PTT. We perform numerical simulations using the symmetric square root representation of the conformation tensor to stabilize the numerical scheme against the high Weissenberg number problem. The variation in flow structures associated with merging and splitting of elongated vortices in shallow cavities and coalescence of corner eddies to yield a second primary vortex in deep cavities with respect to the variation in flow parameters is discussed. We discuss the effect of the dominant eigenvalues and the corresponding eigenvectors on the location of the primary eddy in the cavity. We also demonstrate, by performing numerical simulations for shallow and deep cavities, that where the Deborah number (based on convective time scale) characterizes the elastic behaviour of the fluid in deep cavities, Weissenberg number (based on shear rate) should be used for shallow cavities. (C) 2016 Elsevier B.V. All rights reserved.

A peridynamic theory for linear elastic shells

A state-based peridynamic formulation for linear elastic shells is presented. The emphasis is on introducing, possibly for the first time, a general surface based peridynamic model to represent the deformation characteristics of structures that have one geometric dimension much smaller than the other two. A new notion of curved bonds is exploited to cater for force transfer between the peridynamic particles describing the shell. Starting with the three dimensional force and deformation states, appropriate surface based force, moment and several deformation states are arrived at. Upon application on the curved bonds, such states yield the necessary force and deformation vectors governing the motion of the shell. By incorporating a shear correction factor, the formulation also accommodates analysis of shells that have higher thickness. In order to attain this, a consistent second order approximation to the complementary energy density is considered and incorporated in peridynamics via constitutive correspondence. Unlike the uncoupled constitution for thin shells, a consequence of a first order approximation, constitutive relations for thick shells are fully coupled in that surface wryness influences the in-plane stress resultants and surface strain the moments. Our proposal on the peridynamic shell theory is numerically assessed against simulations on static deformation of spherical and cylindrical shells, that of flat plates and quasi-static fracture propagation in a cylindrical shell. (C) 2016 Elsevier Ltd. All rights reserved.

Morphological transitions and buckling characteristics in a nanoparticle-laden sessile droplet resting on a heated hydrophobic substrate

In this work, we have established the evaporation-liquid flow coupling mechanism by which sessile nanofluid droplets on a hydrophobic substrate evaporate and agglomerate to form unique morphological features under controlled external heating. It is well understood that evaporation coupled with internal liquid flow controls particle transport in a spatiotemporal sense. Flow characteristics inside the heated droplet are investigated and found to be driven by the buoyancy effects. Velocity magnitudes are observed to increase by an order at higher temperatures with similar looking flow profiles. The recirculating flow induced particle transport coupled with collision of particles and shear interaction between them leads to the formation of dome shaped viscoelastic shells of different dimensions depending on the surface temperature. These shells undergo sol-gel transition and subsequently undergo buckling instability leading to the formation of daughter cavities. With an increase in the surface temperature, droplets exhibit buckling from multiple sites over a larger sector in the top half of the droplet. Irrespective of the initial nanoparticle concentration and substrate temperature, growth of a daughter cavity (subsequent to buckling) inside the droplet is found to be controlled by the solvent evaporation rate from the droplet periphery and is shown to exhibit a universal trend.