Theoretical and experimental analysis of longitudinal wave propagation in cylindrical viscoelastic rods

Wave propagation in viscoelastic rods is encountered in many applications including studies of impact and fracture under high strain rates and characterization of the dynamic behavior of viscoelastic materials. For viscoelastic materials, both material and geometric dispersion are possible when the diameter of the rod is of the same order as the wavelength. In this work, we simplify the Pochhammer frequency equation for low and intermediate loss viscoelastic materials and formulate corrections for geometric dispersion for both the phase velocity and attenuation. The formulation is then experimentally verified with measurements of the phase velocity and attenuation in commercial polymethylmethacrylate rods that are 12 and in diameter. Without correcting for geometric dispersion, the usable frequency range for determining the phase velocity and attenuation for the rod is about , and about for the rod. Using the correction procedure developed here, it was possible to accurately determine the phase velocity and attenuation up to frequencies exceeding for the rod and for the rod. These corrections are applicable to many polymers and other viscoelastic materials. From thereon, the viscoelastic properties of the material can be determined over a wide range of frequencies.     

Geometrical properties of the vorticity vector derived using large-eddy simulation

In this paper, the geometrical properties of the resolved vorticity vector derived from large-eddy simulation are investigated using a statistical method. Numerical tests have been performed based on a turbulent Couette channel flow using three different dynamic linear and nonlinear subgrid-scale stress models. The geometrical properties of have a significant impact on various physical quantities and processes of the flow. To demonstrate, we examined helicity and helical structure, the attitude of with respect to the eigenframes of the resolved strain rate tensor and negative subgrid-scale stress tensor -τ_ij, enstrophy generation, and local vortex stretching and compression. It is observed that the presence of the wall has a strong anisotropic influence on the alignment patterns between and the eigenvectors of , and between and the resolved vortex stretching vector. Some interesting wall-limiting geometrical alignment patterns and probability density distributions in the form of Dirac delta functions associated with these alignment patterns are reported. To quantify the subgrid-scale modelling effects, the attitude of with respect to the eigenframe of -τ_ij is studied, and the geometrical alignment between and the Euler axis is also investigated. The Euler axis and angle for describing the relative rotation between the eigenframes of -τ_ij and are natural invariants of the rotation matrix, and are found to be effective for characterizing a subgrid-scale stress model and for quantifying the associated subgrid-scale modelling effects on the geometrical properties of .     

Limit analysis of multi-layered plates. Part I: The homogenized Love-Kirchhoff model

The purpose of this paper is to determine , the overall homogenized Love-Kirchhoff strength domain of a rigid perfectly plastic multi-layered plate, and to study the relationship between the 3D and the homogenized Love-Kirchhoff plate limit analysis problems. In the Love-Kirchhoff model, the generalized stresses are the in-plane (membrane) and the out-of-plane (flexural) stress field resultants. The homogenization method proposed by Bourgeois [1997. Modélisation numérique des panneaux structuraux légers. Ph.D. Thesis, University Aix-Marseille] and Sab [2003. Yield design of thin periodic plates by a homogenization technique and an application to masonry wall. C. R. Méc. 331, 641-646] for in-plane periodic rigid perfectly plastic plates is justified using the asymptotic expansion method. For laminated plates, an explicit parametric representation of the yield surface is given thanks to the π-function (the plastic dissipation power density function) that describes the local strength domain at each point of the plate. This representation also provides a localization method for the determination of the 3D stress components corresponding to every generalized stress belonging to . For a laminated plate described with a yield function of the form , where σ^u is a positive even function of the out-of-plane coordinate x₃ and is a convex function of the local stress σ, two effective constants and a normalization procedure are introduced. A symmetric sandwich plate consisting of two Von-Mises materials ( in the skins and in the core) is studied. It is found that, for small enough contrast ratios (), the normalized strength domain is close to the one corresponding to a homogeneous Von-Mises plate [Ilyushin, A.-A., 1956. Plasticité. Eyrolles, Paris].     

Discrete dislocation plasticity analysis of the wedge indentation of films

The plane strain indentation of single crystal films on a rigid substrate by a rigid wedge indenter is analyzed using discrete dislocation plasticity. The crystals have three slip systems at ±35.3^° and 90^° with respect to the indentation direction. The analyses are carried out for three values of the film thickness, 2, 10 and , and with the dislocations all of edge character modeled as line singularities in a linear elastic material. The lattice resistance to dislocation motion, dislocation nucleation, dislocation interaction with obstacles and dislocation annihilation are incorporated through a set of constitutive rules. Over the range of indentation depths considered, the indentation pressure for the 10 and thick films decreases with increasing contact size and attains a contact size-independent value for contact lengths . On the other hand, for the films, the indentation pressure first decreases with increasing contact size and subsequently increases as the plastic zone reaches the rigid substrate. For the 10 and thick films sink-in occurs around the indenter, while pile-up occurs in the film when the plastic zone reaches the substrate. Comparisons are made with predictions obtained from other formulations: (i) the contact size-independent indentation pressure is compared with that given by continuum crystal plasticity; (ii) the scaling of the indentation pressure with indentation depth is compared with the relation proposed by Nix and Gao [1998. Indentation size effects in crystalline materials: a law for strain gradient plasticity. J. Mech. Phys. Solids 43, 411-423]; and (iii) the computed contact area is compared with that obtained from the estimation procedure of Oliver and Pharr [1992. An improved technique for determining hardness and elastic-modulus using load and displacement sensing indentation experiments, J. Mater. Res. 7, 1564-1583].     

Nanoindentation and the indentation size effect: Kinetics of deformation and strain gradient plasticity

A study of the indentation size effect (ISE) in aluminum and alpha brass is presented. The study employs rate effects to examine the fundamental mechanisms responsible for the ISE. These rate effects are characterized in terms of the rate sensitivity of the hardness, , where H is the hardness and is an effective strain rate in the plastic volume beneath the indenter. can be measured using indentation creep, load relaxation, or rate change experiments. The activation volume V∗, calculated based on which can traditionally be used to compare rate sensitivity data from a hardness test to conventional uniaxial testing, is calculated. Using materials with different stacking fault energy and specimens with different levels of work hardening, we demonstrate how increasing the dislocation density affects V∗; these effects may be taken as a kinetic signature of dislocation strengthening mechanisms. We noticed both H and exhibit an ISE. The course of V∗ vs. H as a result of the ISE is consistent with the course of testing specimens with different level of work hardening. This result was observed in both materials. This suggests that a dislocation mechanism is responsible for the ISE. When the results are fitted to a strain gradient plasticity model, the data at deep indents (microhardness and large nanoindentation) exhibit a straight-line behavior closely identical to literature data. However, for shallow indents (nanoindentation data), the slope of the line severely changes, decreasing by a factor of 10, resulting in a “bilinear behavior”.     

Formation of molten metal films during metal-on-metal slip under extreme interfacial conditions

The present paper describes results of plate-impact pressure-shear friction experiments conducted to study time-resolved growth of molten metal films during dry metal-on-metal slip under extreme interfacial conditions. By employing tribo-pairs comprising hard tool-steel against relatively low melt-point metals such as 7075-T6 aluminum alloys, interfacial friction stress ranging from 100 to and slip speeds of approximately have been generated. These relatively high levels of friction stress combined with high slip-speeds generate conditions conducive for interfacial temperatures to approach the melting point of the lower melt point metal (Al alloy) comprising the tribo-pair.A Lagrangian finite element code is developed to understand the evolution of the thermo-mechanical fields and their relationship to the observed slip response. The code accounts for dynamic effects, heat conduction, contact with friction, and full thermo-mechanical coupling. At temperatures below the melting point the material is described as an isotropic thermally softening elastic-viscoplastic solid. For material elements with temperatures in excess of the melt point a purely Newtonian fluid constitutive model is employed.The results of the hybrid experimental-computational study provides new insights into the thermoelastic-plastic interactions during high speed metal-on-metal slip under extreme interfacial conditions. During the early part of frictional slip the coefficient of kinetic friction is observed to decrease with increasing slip velocity. During the later part transition in interfacial slip occurs from dry metal-on-metal sliding to the formation of molten Al films at the tribo-pair interface. Under these conditions the interfacial resistance approaches the shear strength of the molten aluminum alloy under normal pressures of approximately 1- and shear strain rates of . The results of the study indicate that under these extreme conditions molten aluminum films maintain a shearing resistance as high as .Scanning electron microscopy of the slip surfaces reveal molten aluminum to be smeared on the tribo-pair interface. Knoop hardness measurements in 7075-T6 Al alloy at various depths from the slip interface indicate that the hardness increases approximately linearly with depth and reaches a plateau at approximately from the surface.     

Micromechanics of residual stress and texture development due to poling in polycrystalline ferroelectric ceramics

This paper presents micromechanics based analysis of elastic strain and changes in the texture of poled polycrystalline ferroelectric PZT ceramics for direct comparison with synchrotron X-ray measurements. The grains are modelled as spherical inclusions, to which transformation strains are assigned depending on the fractions of different ferroelectric domains. Eshelby's inclusion problem with the classical self-consistent method is applied to evaluate the elastic state of the grains. In particular, the elongation due to lattice elastic strain is calculated as a function of inclination Ψ relative to the polar axis. The ratio of diffraction peak intensities, corresponding to the domain fractions, is also expressed as a function of Ψ. This analysis identifies the special character of the reflection, for which the lattice strain along in the stress free state is independent of ferroelectric domain population and hence unaffected by poling. The elongation due to the lattice strain parallel to and peak intensity ratio are expressed in terms of the overall macroscopic strain of a poled specimen, each having a dependence.     

Deformation of FCC nanowires by twinning and slip

We present atomistic simulations of the tensile and compressive loading of single crystal face-centered cubic (FCC) nanowires with and orientations to study the propensity of the nanowires to deform via twinning or slip. By studying the deformation characteristics of three FCC materials with disparate stacking fault energies (gold, copper and nickel), we find that the deformation mechanisms in the nanowires are a function of the intrinsic material properties, applied stress state, axial crystallographic orientation and exposed transverse surfaces. The key finding of this work is the first order effect that side surface orientation has on the operant mode of inelastic deformation in both and nanowires. Comparisons to expected deformation modes, as calculated using crystallographic Schmid factors for tension and compression, are provided to illustrate how transverse surface orientations can directly alter the deformation mechanisms in materials with nanometer scale dimensions.     

A new analytical technique to find periodic solutions of non-linear systems

Based on the classical harmonic balance method a new technique is presented to determine higher approximate periodic solutions of the non-linear differential equations. The new method is systematic and simple. The solution covers the general initial value problem (i.e., for while the existing solution is determined for a particular case, especially for . The solution is easily transformed to perturbation solution. The method is used in various non-linear problems possessing second and more than second derivatives.     

Complementary energy release rates and dual conservation integrals in micropolar elasticity

The complementary energy momentum tensor, expressed in terms of the spatial gradients of stress and couple-stress, is used to construct the and conservation integrals of infinitesimal micropolar elasticity. The derived integrals are related to the release rates of the complementary potential energy associated with a defect translation or rotation. A nonconserved integral is also derived and related to the energy release rate that is associated with a self-similar cavity expansion. The results are compared to those obtained on the basis of the classical energy momentum tensor, expressed in terms of the spatial gradients of displacement and rotation, and the release rates of the potential energy. It is shown that the evaluation of the complementary conservation integrals is of similar complexity to that of the classical conservation integrals, so that either can be effectively used in the energetic analysis of the mechanics of defects. The two-dimensional versions of the dual conservation integrals are then derived and applied to an out-of-plane shearing of a long cracked slab.     

Large electrostrictive actuation of barium titanate single crystals

An experimental investigation of the electromechanical behavior of single crystals of the ferroelectric perovskite barium titanate is presented. An experimental setup has been designed to investigate large strain actuation in single crystal ferroelectrics subjected to combined electrical and mechanical loading. Experiments have been performed on initially single domain crystals of barium titanate with (100) and (001) orientation at compressive stresses between 0 and . Global strain and polarization histories have been recorded. The electrostrictive response is shown to be highly dependent on the level of applied stress with a maximum strain of 0.9% measured at a compressive stress of about and electric field of about . This level of strain is about 5 times higher than in typical commercial piezoelectric PZT. Polarized light microscopy has been used to observe the evolution of the domain pattern simultaneously with the strain and polarization measurement. The observations reveal that the observed large strain behavior is the result of 90° domain switching.     

Glass-modified stress waves for adhesion measurement of ultra thin films for device applications

Laser-generated stress wave profiles with rarefaction shocks (almost zero post-peak decay times) have been uncovered in different types of glasses and presented in this communication. The rise time of the pulses was found to increase with their amplitude, with values reaching as high as . This is in contrast to measurements in other brittle crystalline solids where pulses with rise times of and post-peak decay times of were recorded. The formation of rarefaction shock is attributed to the increased compressibility of glasses with increasing pressures. This was demonstrated using a one-dimensional nonlinear elastic wave propagation model in which the wave speed was taken as a function of particle velocity. The technological importance of these pulses in measuring the tensile strength of very thin film interfaces is demonstrated by using a previously developed laser spallation experiment in which a laser-generated compressive stress pulse in the substrate reflects into a tensile wave from the free surface of the film and pries off its interface at a threshold amplitude. Because of the rarefaction shock, glass-modified waves allow generation of substantially higher interfacial tensile stress amplitudes compared with those with finite post-peak decay profiles. Thus, for the first time, tensile strengths of very strong and ultra thin film interfaces can be measured. Results presented here indicate that interfaces of 185-nm-thick films, and with strengths as high as , can be measured. Thus, an important advance has been made that should allow material optimization of ultra thin layer systems that may form the basis of future MEMS-based microelectronic, mechanical and clinical devices.     

A gradient theory of small-deformation isotropic plasticity that accounts for the Burgers vector and for dissipation due to plastic spin

This study develops a gradient theory of small-deformation viscoplasticity based on: a system of microforces consistent with its peculiar balance; a mechanical version of the second law that includes, via the microforces, work performed during viscoplastic flow; a constitutive theory that accounts for the Burgers vector through a free energy dependent on , with H^p the plastic part of the elastic-plastic decomposition of the displacement gradient. The microforce balance and the constitutive equations, restricted by the second law, are shown to be together equivalent to a nonlocal flow rule in the form of a coupled pair of second-order partial differential equations. The first of these is an equation for the plastic strain-rate in which the stress T plays a basic role; the second, which is independent of T, is an equation for the plastic spin. A consequence of this second equation is that the plastic spin vanishes identically when the free energy is independent of, but not generally otherwise. A formal discussion based on experience with other gradient theories suggests that sufficiently far from boundaries solutions should not differ appreciably from classical solutions, but close to microscopically hard boundaries, boundary layers characterized by a large Burgers vector and large plastic spin should form.Because of the nonlocal nature of the flow rule, the classical macroscopic boundary conditions need be supplemented by nonstandard boundary conditions associated with viscoplastic flow. As an aid to solution, a variational formulation of the flow rule is derived.Finally, we sketch a generalization of the theory that allows for isotropic hardening resulting from dissipative constitutive dependences on .