Bounds on stress concentration factors in finite anti-plane shear

This paper provides an analytical approach for obtaining bounds on elastic stress concentration factors in the theory of finite anti-plane shear of homogeneous, isotropic, incompressible materials. The problem of an infinite slab with traction-free circular cavity subject to a state of finite simple shear deformation is considered. Explicit estimates are obtained for the maximum shearing stress in terms of the applied stress at infinity and constitutive parameters. The analysis is based on application of maximum principles for second-order quasilinear uniformly elliptic equations.     

Effect of material non-homogeneity on the inhomogeneous shearing deformation of a Gent slab subjected to a temperature gradient

It is well known that most rubber-like materials are non-homogeneous due to either imperfect manufacturing conditions or the action of severe thermo-oxidative environments in many practical applications. In this study, within the context of finite thermoelasticity, we theoretically analyze the inhomogeneous shearing deformation of a non-homogeneous rubber-like slab subjected to a thermal gradient across its thickness. The major objective of this study is to investigate the effect of the material non-homogeneity, which is the material-coordinate dependence of the material response functions, on the stress-strain fields for a given temperature gradient. First, we show the existence of a simple shearing deformation from which the generalized shear modulus and the generalized thermal conductivity of the slab could be obtained. Based on this information, the Gent material model is generalized to take the material non-homogeneity and the temperature dependence of the stress into account. To analyze the inhomogeneous shearing deformation of the non-homogeneous slab, deformation and temperature fields are postulated; then the decoupled temperature field is obtained analytically by solving the local energy balance equation. Finally, the static equilibrium equations are solved considering the linear temperature field. Our results show that the spatial pattern and the degree of the material non-homogeneity have profound effects on the stress-strain fields. The shear strain becomes nearly homogeneous and the stresses are relatively small for a certain spatial variation of the material non-homogeneity. This result suggests the possibility of designing a novel class of materials: functionally graded rubber-elastic materials (FGREMs).     

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 II 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.     

A mesoscale investigation of strain rate effect on dynamic deformation of single-crystal copper

A combined finite element (FE) simulation and discrete dislocation dynamics (DD) approach has been developed in this paper to investigate the dynamic deformation of single-crystal copper at mesoscale. The DD code yields the plastic strain based on the slip of dislocations and serves as a substitute for the 3D constitutive form used in the usual FE computation, which is implemented into ABAQUS/Standard with a user-defined material subroutine. On the other hand, the FE code computes the displacement and stress field during the dynamic deformation. The loading rate effects on the yield stress and the deformation patterning of single-crystal copper are investigated. With the increasing of strain rate, the yield stress of single-crystal copper increases rapidly. A critical strain rate exists in each single-crystal copper block for the given size and dislocation sources, below which the yield stress is relatively insensitive to the strain rate. The dislocation patterning changes from non-uniform to uniform under high-strain-rate. The shear stresses in the bands are higher than that in the neighboring regions, which are formed shear bands in the crystal. The band width increases with the strain rate, which often take places where the damage occurs.     

Post-critical plastic deformation in incrementally nonlinear materials

The formation of multiple macroscopic shear bands is investigated as a mechanism of advanced plastic flow of polycrystalline metals. The overall deformation pattern and material characteristics are determined beyond the critical instant of ellipticity loss, without the need of introducing an internal length scale. This novel approach to the modelling of post-critical plastic deformation is based on the concept of a representative nonuniform solution in a homogeneous material. The indeterminacy of a post-critical representative solution is removed by eliminating unstable solution paths with the help of the energy criterion of path instability. It is shown that the use of micromechanically based, incrementally nonlinear corner theories of time-independent plasticity leads then to gradual concentration of post-critical plastic deformation. The volume fraction occupied by shear bands is found to have initially a well-defined, finite value insensitive to the mesh size in finite element calculations. Further deformation depends qualitatively on details of the constitutive law. In certain cases, the volume fraction of active bands decreases rapidly to zero, leading to material instability of dynamic type. However, for physically hardening materials with the yield-vertex effect, the localization volume typically remains finite over a considerable deformation range. At later stages of the plane strain simulation, differently aligned secondary bands are formed in a series of bifurcations.     

Finite anti-plane shear of an infinite slab with a traction-free elliptical cavity: Bounds on the stress concentration factor

This paper provides an analytical approach for obtaining bounds on elastic stress concentration factors in the theory of finite anti-plane shear of homogeneous, isotropic, incompressible materials. The problem of an infinite slab with traction-free elliptical cavity subject to a state of finite simple shear deformation is considered. Explicit estimates are obtained for the maximum shear stress in terms of the cavity geometry, applied stress at infinity and constitutive parameters. The analysis is based on application of maximum principles for second-order quasilinear uniformly elliptic equations.     

Lüders bands propagation of 1045 steel under multiaxial stress state

Inhomogeneous plastic deformation of 1045 steel under monotonic loading was experimentally studied. Thin-walled tubular specimens were used in the experiments and custom-made small strain gages were bonded on the specimen surface to characterize the local deformation. Experiments were conducted under tension, torsion, and combined tension–torsion. During the propagation of Lüders bands, the local deformation experienced two-stage deformation: an abrupt plastic deformation stage followed by a slower deformation process. In some area of the gage section of the specimen, a small amount of initial plastic deformation occurred before the Lüders front reached. During the propagation of Lüders bands, multiple Lüders fronts can be formed. Under tension, torsion, and combined tension–torsion with a constant axial load, the Lüders front was approximately parallel to the material plane of maximum shear stress. When the combined axial-torsion followed a proportional fashion, the stress–extensometer strain responses were dependent on the axial/torsional loading ratio, and the Lüders fronts were oriented differently and propagated along the specimen axis at a different velocity. The local strain was inhomogeneous even at the work-hardening stage. The relationships between the equivalent stress and the equivalent plastic strain were found to be practically identical for all the loading cases studied.     

Analysis of simple shear endochronic equations for finite plastic deformation

IntroductionThestress_strainbehaviorofmaterialswithfiniteplasticdeformationisaninterestingissue ,onwhichsignificantprogresshasbeenmadethroughboththephenomenologicalandphysicalapproaches.Thephenomenologicalapproachisbasedoncontinuummechanicsofplasticity .Ithasitsadvantageinsolvingcomplicatedproblemsbecauseofitssimplicity .Mostofphenomenologicaltheoriesareinvolvedintheconceptofcorotationalrates.Thematerialderivativeofstresswasnotobjectiveunderfinitedeformation .TheJaumannratewasusuallyusedbefo…     

金属材料平面应变拉伸变形局部化有限元分析

本文对于涉及韧性金属大变形中颈缩与剪切带断裂一类高度非线性变形局部化问题进行了弹塑性有限元数值模拟。采用改进的Ｊ２形变理论微分形式公式与交叉三角形四边形单元有限元网格，详细研究了应变硬化指数及初始表面不均匀特性的平面应变拉伸颈缩和剪切带形成的综合影响，给出此类问题的断裂机制图。     

A scaling law for the effect of inertia on the formation of adiabatic shear bands

For a rigid/perfectly plastic material with linear thermal softening and power law rate hardening there is a competition between heat conduction and inertia in determining the time of shear band formation. In a finite specimen the nominal strain rate that produces the fastest growth of perturbations corresponds to the minimum critical strain. Similarly for a fixed strain rate in an infinite specimen, there is a finite wavelength with the maximum growth rate. It is argued that this wavelength should correspond to the most probable minimum spacing for shear bands.     

Adiabatic shear banding in plane strain tensile deformations of 11 thermoelastoviscoplastic materials with finite thermal wave speed

We study the initiation and propagation of adiabatic shear bands (ASBs) in 11 homogeneous materials each modeled as microporous, isotropic and thermoelastoviscoplastic, and deformed in plane strain tension. The heat conduction in each material is assumed to be governed by a hyperbolic heat equation; thus thermal and mechanical waves propagate with finite speeds. The decrease in the thermophysical parameters due to the increase in porosity is considered. An ASB is assumed to initiate at a material point when the maximum shear stress there has dropped to 80% of its peak value for that material point and it is deforming plastically. An approximate solution of the coupled nonlinear partial differential equations subject to suitable initial and boundary conditions is found by the finite element method (FEM). In contrast to the Considerè and the Hart criterion, it is found that an ASB initiates when the axial load drops rapidly and not when it peaks. The refinement of the 40 × 40 uniform FE mesh to 120 × 120 uniform elements decreased the ASB initiation time by 2.1% while increasing the CPU time by a factor of ∼26. By locating points where the ASB has initiated we find its current length, width and speed. The 11 materials are ranked according to the time of initiation of an ASB under otherwise identical geometric and loading conditions with the same initial nonuniform porosity distribution. This ranking of materials is found to differ somewhat from that ascertained by Batra and Kim (1992) who studied simple shearing deformations, and by Batra et al. (1995) who analyzed three-dimensional torsional deformations of thin-walled tubular specimens. The average axial strain determined from the maximum axial load condition differs noticeably from that when an ASB initiates.     

Size effects and idealized dislocation microstructure at small scales: Predictions of a Phenomenological model of Mesoscopic Field Dislocation Mechanics: Part II

In Part I of this set of two papers, a model of mesoscopic plasticity is developed for studying initial-boundary value problems of small scale plasticity. Here we make qualitative, finite element method-based computational predictions of the theory. We demonstrate size effects and the development of strong inhomogeneity in simple shearing of plastically constrained grains. Non-locality in elastic straining leading to a strong Bauschinger effect is analyzed. Low shear strain boundary layers in constrained simple shearing of infinite layers of polycrystalline materials are not predicted by the model, and we justify the result based on an examination of the no-dislocation-flow boundary condition. The time-dependent, spatially homogeneous, simple shearing solution of PMFDM is studied numerically. The computational results and an analysis of continuous dependence with respect to initial data of solutions for a model linear problem point to the need for a nonlinear study of a stability transition of the homogeneous solution with decreasing grain size and increasing applied deformation. The continuous-dependence analysis also points to a possible mechanism for the development of spatial inhomogeneity in the initial stages of deformation in lower-order gradient plasticity theory. Results from thermal cycling of small scale beams/films with different degrees of constraint to plastic flow are presented showing size effects and reciprocal-film-thickness scaling of dislocation density boundary layer width. Qualitative similarities with results from discrete dislocation analyses are noted where possible.We discuss the convergence of approximate solutions with mesh refinement and its implications for the prediction of dislocation microstructure development, motivated by the notion of measure-valued solutions to conservation laws.     

Shear-band bifurcation of an isotropic compressible hyperelastic solid

This paper concerns shear-band bifurcations from the homogeneous finite plane deformation of an isotropic compressible hyperelastic solid. The governing equations for the incremental plane deformation superposed on the initial finite deformation are derived and then the equilibrium equations in terms of incremental displacements are classified into the elliptic type, parabolic type, etc. From this classification follows a restriction which should be placed on the strain-energy function in order that the equilibrium equations may be either elliptic or parabolic for all principal stretches. For the hyperelastic solid complying with this restriction, the condition for the shear-band bifurcation is obtained. Finally, the incremental displacement field of an infinite series of shear bands in a slab sandwiched between slippery rigid layers is established.     

Constitutive relations for the shear band evolution in granular matter under large strain

A so-called ＂split-bottom ring shear cell＂ leads to wide shear bands under slow, quasi-static deformation. Unlike normal cylindrical Couette shear cells or rheometers, the bottom plate is split such that the outer part of it can move with the outer wall, while the other part （inner disk） is immobile. From discrete element simulations （DEM）, several continuum fields like the density, velocity, deformation gradient and stress are computed and evaluated with the goal to formulate objective constitutive relations for the powder flow behavior. From a single simulation, by applying time- and （local） space-averaging, a non-linear yield surface is obtained with peculiar stress dependence. The anisotropy is always smaller than the macroscopic friction coefficient. However, the lower bound of anisotropy increases with the strain rate, approaching the maximum according to a stretched exponential with a specific rate that is consistent with a shear path of about one particle diameter.     

Constitutive relations for the shear band evolution in granular matter under large strain

A so-called "split-bottom ring shear cell" leads to wide shear bands under slow, quasi-static deformation. Unlike normal cylindrical Couette shear cells or rheometers, the bottom plate is split such that the outer part of it can move with the outer wail, while the other part (inner disk) is immobile. From discrete element simulations (DEM), several continuum fields like the density, velocity, deformation gradient and stress are computed and evaluated with the goal to formulate objective constitutive relations for the powder flow behavior. From a single simulation, by applying time-and (local) space-averaging, a non-linear yield surface is obtained with peculiar stress dependence.The anisotropy is always smaller than the macroscopic friction coefficient. However, the lower bound of anisotropy increases with the strain rate, approaching the maximum according to a stretched exponential with a specific rate that is consistent with a shear path of about one particle diameter.     

Microstructural effects on shear instability and shear band spacing

A constitutive relation that accounts for the thermally activated dislocation motion and microstructure interaction is used to study the stability of a homogeneous solution of equations governing the simple shearing deformations of a thermoviscoplastic body. An instability criterion and an upper bound for the growth rate of the infinitesimal deformations superimposed on the homogeneous solution are derived. By adopting Wright and Ockendon's postulate, i.e., the wavelength of the dominant instability mode with the maximum growth rate determines the minimum spacing between shear bands, the shear band spacing is computed. The effect of the initial dislocation density, the nominal strain-rate, and parameters describing the initial thermal activation and the initial microstructure interaction on the shear band spacing are delineated.     

Initiation of localized plane deformations at a circular cavity in an infinite compressible nonlinearly elastic medium

This investigation is concerned with the plane strain deformation of an infinite slab, containing a circular cavity, within the theory of finite elastostatics for a particular homogeneous isotropic compressible material, the so-called Blatz-Ko material. The body is subjected to uniform pressure, either internal or external. Exact closed-form solutions for the axisymmetric deformation and stress fields are obtained. In the case of internal pressure, it is found that the applied pressure may not exceed a certain maximum value p _max. At a value of pressure p _e (<p _max), the governing equations lose ellipticity at the cavity wall. For greater values of pressure this solution remains smooth, though involving both elliptic and non-elliptic regions. Non-existence of axisymmetric solutions with discontinuous strain fields is established. The possibility of bifurication into a surface mode is considered and it is shown that this occurs at a value of pressure slightly smaller than p _e. Such surface wrinkling leads to a periodic distribution of points of stress concentration, from which shear bands may initiate.This work was supported by the U.S. Army Research Office under Grant DAAG29-83-K-0145 (R.A. & C.O.H.) and by the U.S. National Science Foundation under Grant MEA 78-26071 (C.O.H.).     

An experimental and numerical study of the localization behavior of tantalum and stainless steel

In general, the shear localization process involves initiation and growth. Initiation is expected to be a stochastic process in material space where anisotropy in the elastic–plastic behavior of single crystals and inter-crystalline interactions serve to form natural perturbations to the material’s local stability. A hat-shaped sample geometry was used to study shear localization growth. It is an axi-symmetric sample with an upper “hat” portion and a lower “brim” portion with the shear zone located between the hat and brim. The shear zone length is 870–890 μm with deformation imposed through a split-Hopkinson pressure bar system at maximum top-to-bottom velocity in the range of 8–25 m/s. We present experimental results of the deformation response of tantalum and 316L stainless steel samples. The tantalum samples did not form shear bands but the stainless steel sample formed a late stage shear band. We have also modeled these experiments using both conductive and adiabatic continuum models. An anisotropic elasto-viscoplastic constitutive model with damage evolution was used within the finite element code EPIC. A Mie-Gruneisen equation of state and the rate and temperature sensitive MTS flow stress model together with a Gurson flow surface were employed. The models performed well in predicting the experimental data. The numerical results for tantalum suggested a maximum equivalent strain rate on the order of 7 × 10⁴ s⁻¹ in the gage section for an imposed top surface displacement rate of 17.5 m/s. The models also suggested that for an initial temperature of 298 K a temperature in the neighborhood of 900 K was reached within the shear section. The numerical results for stainless steel suggest that melting temperature was reached throughout the shear band shortly after peak load. Due to sample geometry, the stress state in the shear zone was not pure shear; a significant normal stress relative to the shear zone basis line was developed.     

Application of corotational rates of the logarithmic strain in constitutive modeling of hardening materials at finite deformations

In this paper a finite deformation constitutive model for rigid plastic hardening materials based on the logarithmic strain tensor is introduced. The flow rule of this constitutive model relates the corotational rate of the logarithmic strain tensor to the difference of the deviatoric Cauchy stress and the back stress tensors. The evolution equation for the kinematic hardening of this model relates the corotational rate of the back stress tensor to the corotational rate of the logarithmic strain tensor. Using Jaumann, Green–Naghdi, Eulerian and logarithmic corotational rates in the proposed constitutive model, stress–strain responses and subsequent yield surfaces are determined for rigid plastic kinematic and isotropic hardening materials in the simple shear problem at finite deformations.     

Three-dimensional micromechanical modeling of voided polymeric materials

A three-dimensional micromechanical unit cell model for particle-filled materials is presented. The cell model is based on a Voronoi tessellation of particles arranged on a body-centered cubic (BCC) array. The three-dimensionality of the present cell model enables the study of several deformation modes, including uniaxial, plane strain and simple shear deformations, as well as arbitrary principal stress states.The unit cell model is applied to studies on the micromechanical and macromechanical behavior of rubber-toughened polycarbonate. Different load cases are examined, including plane strain deformation, simple shear deformation and principal stress states. For a constant macroscopic strain rate, the different load cases show that the macroscopic flow strength of the blend decreases with an increase in void volume fraction, as expected. The main mechanism for plastic deformation is broad shear banding across inter-particle ligaments. The distributed nature of plastic straining acts to reduce the amount of macroscopic strain softening in the blend as the initial void volume fraction is increased. In the case of plane strain deformation, the plastic flow is observed to initiate across inter-particle ligaments in the direction of constraint. This particular mode of deformation could not have been captured using a two-dimensional, plane strain idealization of cylindrical voids in a matrix.The potential for localized crazing and/or cavitation in the matrix is addressed. It is observed that the introduction of voids acts to relieve hydrostatic stress in the matrix material, compared to the homopolymer. It is also seen that the predicted peak hydrostatic stress in the matrix is higher under plane strain deformation than under triaxial tension (with equal lateral stresses), for the same macroscopic stress triaxiality.The effect of void volume fraction on the macroscopic uniaxial tension behavior of the different blends is examined using a Considère construction for dilatant materials. The natural draw ratio was predicted to decrease with an increase in void volume fraction.