CROSS-SECTIONAL ANALYSIS OF PRE-TWISTED THICK BEAMS USING VARIATIONAL ASYMPTOTIC METHOD

An asymptotically-exact methodology is presented for obtaining the cross-sectional stiffness matrix of a pre-twisted moderately-thick beam having rectangular cross sections and made of transversely isotropic materials. The anisotropic beam is modeled from 3-D elasticity, without any further assumptions. The beam is allowed to have large displacements and rotations, but small strain is assumed. The strain energy of the beam is computed making use of the constitutive law and the kinematical relations derived with the inclusion of geometrical nonlinearities and initial twist. Large displacements and rotations are allowed, but small strain is assumed. The Variational Asymptotic Method is used to minimize the energy functional, thereby reducing the cross section to a point on the reference line with appropriate properties, yielding a 1-D constitutive law. In this method as applied herein, the 2-D cross-sectional analysis is performed asymptotically by taking advantage of a material small parameter and two geometric small parameters. 3-D strain components are derived using kinematics and arranged as orders of the small parameters. Warping functions are obtained by the minimization of strain energy subject to certain set of constraints that renders the 1-D strain measures well-defined. Closed-form expressions are derived for the 3-D non-linear warping and stress fields. The model is capable of predicting interlaminar and transverse shear stresses accurately up to first order.

Analysis of Transverse Shear Strains in Pre-Twisted Thick Beams using Variational Asymptotic Method

The cross-sectional stiffness matrix is derived for a pre-twisted, moderately thick beam made of transversely isotropic materials and having rectangular cross sections. An asymptotically-exact methodology is used to model the anisotropic beam from 3-D elasticity, without any further assumptions. The beam is allowed to have large displacements and rotations, but small strain is assumed. The strain energy is computed making use of the beam constitutive law and kinematical relations derived with the inclusion of geometrical nonlinearities and an initial twist. The energy functional is minimized making use of the Variational Asymptotic Method (VAM), thereby reducing the cross section to a point on the beam reference line with appropriate properties, forming a 1-D constitutive law. VAM is a mathematical technique employed in the current problem to rigorously split the 3-D analysis of beams into two: a 2-D analysis over the beam cross-sectional domain, which provides a compact semi-analytical form of the properties of the cross sections, and a nonlinear 1-D analysis of the beam reference curve. In this method, as applied herein, the cross-sectional analysis is performed asymptotically by taking advantage of a material small parameter and two geometric small parameters. 3-D strain components are derived using kinematics and arranged in orders of the small parameters. Closed-form expressions are derived for the 3-D non-linear warping and stress fields. Warping functions are obtained by the minimization of strain energy subject to certain set of constraints that render the 1-D strain measures well-defined. The zeroth-order 3-D warping field thus yielded is then used to integrate the 3-D strain energy density over the cross section, resulting in the 1-D strain energy density, which in turn helps identify the corresponding cross-sectional stiffness matrix. The model is capable of predicting interlaminar and transverse shear stresses accurately up to first order.

Dependence of shear yield strain and shear transformation zone on the glass transition temperature in thin film metallic glasses

The dependence of shear yield strain, the activation energy and volume of shear transformation zone on the glass transition temperature was investigated through the analysis of statistical distributions of the first pop-in events during spherical indentation of four different thin film metallic glasses. Only the Cu-Zr metallic glass exhibits a bimodal distribution of the first pop-in loads, whereas W-Ru-B, Zr-Cu-Ni-Al and La-Co-Al metallic glasses show an unimodal distribution. Results show that shear yield strain and activation energy of shear transformation zone decrease whereas the volume of shear transformation zone increases with increasing homologous temperature, indicating that it is the activation energy rather than the volume of shear transformation zone that controls shear yield strain. (C) 2015 Elsevier B.V. All rights reserved.

Steady incompressible flow of cohesionless granular materials through a wedge-shaped hopper: Frictional 14kinetic solution to the smooth wall, radial gravity problem

Hybrid frictional-kinetic equations are used to predict the velocity, grain temperature, and stress fields in hoppers. A suitable choice of dimensionless variables permits the pseudo-thermal energy balance to be decoupled from the momentum balance. These balances contain a small parameter, which is analogous to a reciprocal Reynolds number. Hence an approximate semi-analytical solution is constructed using perturbation methods. The energy balance is solved using the method of matched asymptotic expansions. The effect of heat conduction is confined to a very thin boundary layer near the exit, where it causes a marginal change in the temperature. Outside this layer, the temperature T increases rapidly as the radial coordinate r decreases. In particular, the conduction-free energy balance yields an asymptotic solution, valid for small values of r, of the form T proportional r-4. There is a corresponding increase in the kinetic stresses, which attain their maximum values at the hopper exit. The momentum balance is solved by a regular perturbation method. The contribution of the kinetic stresses is important only in a small region near the exit, where the frictional stresses tend to zero. Therefore, the discharge rate is only about 2.3% lower than the frictional value, for typical parameter values. As in the frictional case, the discharge rate for deep hoppers is found to be independent of the head of material.

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.

Stress and strain-driven algorithmic formulations for finite strain viscoplasticity for hybrid and standard finite elements

This work deals with the formulation and implementation of finite deformation viscoplasticity within the framework of stress-based hybrid finite element methods. Hybrid elements, which are based on a two-field variational formulation, are much less susceptible to locking than conventional displacement-based elements. The conventional return-mapping scheme cannot be used in the context of hybrid stress methods since the stress is known, and the strain and the internal plastic variables have to be recovered using this known stress field.We discuss the formulation and implementation of the consistent tangent tensor, and the return-mapping algorithm within the context of the hybrid method. We demonstrate the efficacy of the algorithm on a wide range of problems.

Mixed mode (I and II) crack tip fields in bulk metallic glasses

Stationary crack tip fields in bulk metallic glasses under mixed mode (I and II) loading are studied through detailed finite element simulations assuming plane strain, small scale yielding conditions. The influence of internal friction or pressure sensitivity on the plastic zones. notch deformation, stress and plastic strain fields is examined for different mode mixities. Under mixed mode loading, the notch deforms into a shape such that one part of its surface sharpens while the other part blunts. Increase in mode If component of loading dramatically enhances the normalized plastic zone size, lowers the stresses but significantly elevates the plastic strain levels near the notch tip. Higher internal friction reduces the peak tangential stress but increases the plastic strain and stretching near the blunted part of the notch. The simulated shear bands are straight and extend over a long distance ahead of the notch tip under mode II dominant loading. The possible variations of fracture toughness with mode mixity corresponding to failure by brittle micro-cracking and ductile shear banding are predicted employing two simple fracture criteria. The salient results from finite element simulations are validated by comparison with those from mixed mode (I and II) fracture experiments on a Zr-based bulk metallic glass.

Indentation of sintered alumina clay composite - an experimental discussion of a combined Boussinesq and blister stress field approach

Surface elastic strain field generated in conical indentation of sintered alumina clay composite was measured to verify the suitability of a superposed combination of Boussinesq and blister stress fields, used previously for analysing the indentation problem. The residual strain measured in the elastic hinterland is used to estimate the blister field strength without any reference to stress relation within that field. The approach may be useful in fracture studies of brittle materials.

Molecular dynamic simulations of elastomer structure and its influence on anisotropic stress under time-varying strain

We investigate the evolution of polymer structure and its influence on uniaxial anisotropic stress under time-varying uniaxial strain, and the role of external control variables such as temperature, strain rate, chain length, and density, using molecular dynamics simulation. At temperatures higher than glass transition, stress anisotropy in the system is reduced even though the bond stretch is greater at higher temperatures. There is a significant increase in the stress level with increasing density. At higher densities, the uncoiling of the chains is suppressed and the major contribution to the deformation is by internal deformation of the chains. At faster rates of loading stress anisotropy increases. The deformation mechanism is mostly due to bond stretch and bond bending rather than overall shape and size. Stress levels increase with longer chain length. There is a critical value of the functionality of the cross-linkers beyond which the uniaxial stress developed increases caused primarily by bond stretching due to increased constraint on the motion of the monomers. Stacking of the chains in the system also plays a dominant role in the behaviour in terms of excluded volume interactions. Low density, high temperature, low values of functionality of cross-linkers, and short chain length facilitate chain uncoiling and chain slipping in cross-linked polymers.

Modeling and analysis of a tunable piezoelectric structure for transverse shear wave generation

An analytical investigation of the transverse shear wave mode tuning with a resonator mass (packing mass) on a Lead Zirconium Titanate (PZT) crystal bonded together with a host plate and its equivalent electric circuit parameters are presented. The energy transfer into the structure for this type of wave modes are much higher in this new design. The novelty of the approach here is the tuning of a single wave mode in the thickness direction using a resonator mass. First, a one-dimensional constitutive model assuming the strain induced only in the thickness direction is considered. As the input voltage is applied to the PZT crystal in the thickness direction, the transverse normal stress distribution induced into the plate is assumed to have parabolic distribution, which is presumed as a function of the geometries of the PZT crystal, packing mass, substrate and the wave penetration depth of the generated wave. For the PZT crystal, the harmonic wave guide solution is assumed for the mechanical displacement and electric fields, while for the packing mass, the former is solved using the boundary conditions. The electromechanical characteristics in terms of the stress transfer, mechanical impedance, electrical displacement, velocity and electric field are analyzed. The analytical solutions for the aforementioned entities are presented on the basis of varying the thickness of the PZT crystal and the packing mass. The results show that for a 25% increase in the thickness of the PZT crystal, there is ~38% decrease in the first resonant frequency, while for the same change in the thickness of the packing mass, the decrease in the resonant frequency is observed as ~35%. Most importantly the tuning of the generated wave can be accomplished with the packing mass at lower frequencies easily. To the end, an equivalent electric circuit, for tuning the transverse shear wave mode is analyzed.

Mixed-mode stress intensity factors for arbitrarily oriented cracks in cylindrical shells under axial pull and torsion

Results of photoelastic investigations conducted on cylindrical tubes (made of Araldite material) containing cracks oriented at 0°, 30°, 45°, 60° and 90° to the axis of the tube and subjected to axial and torsional loads are reported. The stress-intensity factors (SIFs) were determined by analysing the crack-tip stress fields. Smith and Smith's method [Engng Fracture Mech.4, 357–366 (1972)] and a new method developed by the authors by modifying Rakesh et al.'s method [Proc. 26th Congress of ISTAM, India (1981)] were employed to evaluate the mixed-mode SIFs.

Asymptotic and finite element analyses of mode III dynamic crack growth at a ductile-brittle interface

In this work, dynamic crack growth along a ductile-brittle interface under anti-plane strain conditions is studied. The ductile solid is taken to obey the J(2) flow theory of plasticity with linear isotropic strain hardening, while the substrate is assumed to exhibit linear elastic behavior. Firstly, the asymptotic near-tip stress and velocity fields are derived. These fields are assumed to be variable-separable with a power singularity in the radial coordinate centered at the crack tip. The effects of crack speed, strain hardening of the ductile phase and mismatch in elastic moduli of the two phases on the singularity exponent and the angular functions are studied. Secondly, full-field finite element analyses of the problem under small-scale yielding conditions are performed. The validity of the asymptotic fields and their range of dominance are determined by comparing them with the results of the full-field finite element analyses. Finally, theoretical predictions are made of the variations of the dynamic fracture toughness with crack velocity. The influence of the bi-material parameters on the above variation is investigated.

Effect of the Edge Type and Strain on the Structural, Electronic and Magnetic Properties of the BNRs

We present the effect of edge structures on the edge energy and stress of BN nanoribbons. Ab initio density functional calculations show that the armchair edge is lower in energy than the zigzag edge by 0.43 eV/angstrom. Both types of the edges are under the compressive stress. The zigzag edges are mechanically more stable than the armchair edges. Based on the calculated edge energies, the equilibrium shape of the BN flakes are found to be regular hexagonal, and dominated by the armchair edges. The zigzag ribbons are found to be half-metallic, whereas the armchair ribbons are semiconducting.

Analytical model for stress-strain response of plastic waste mixed soil

Recycling plastic water bottles has become one of the major challenges world wide. The present study provides an approach for the use of plastic waste as reinforcement material in soil, which can be used for ground improvement, subbases, and subgrade preparation in road construction. The experimental results are presented in the form of stress-strain-pore water pressure response and compression paths. On the basis of experimental test results, it is observed that the strength of soil is improved and compressibility reduced significantly with the addition of a small percentage of plastic waste to the soil. In this paper, an analytical model is proposed to evaluate the response of plastic waste mixed soil. It is noted that the model captures the stress-strain and pore water pressure response of all percentages of plastic waste adequately. The paper also provides a comparative study of failure stress obtained from different published models and the proposed model, which are compared with experimental results. The improvement in strength attributable to the inclusion of plastic waste can be advantageously used in ground improvement projects.

Nonlinear analysis of a thin pre-twisted and delaminated anisotropic strip

The present work is aimed at the development of an efficient mathematical model to assess the degradation in the stiffness properties of an anisotropic strip due to delamination. In particular, the motive is to capture those nonlinear effects in a strip that arise due to the geometry of the structure, in the presence of delamination. The variational asymptotic method (VAM) is used as a mathematical tool to simplify the original 3D problem to a 1D problem. Further simplification is achieved by modeling the delaminated structure by a sublaminate approach. By VAM, a 2D nonlinear sectional analysis is carried out to determine compact expression for the stiffness terms. The stiffness terms, both linear and nonlinear, are derived as functions of delamination length and location in closed form. In general, the results from the analysis include fully coupled nonlinear 1D stiffness coefficients, 3D strain field, 3D stress field, and in-plane and warping fields. In this work, the utility of the model is demonstrated for a static case, and its capability to capture the trapeze effect in the presence of delamination is investigated and compared with results available in the literature.