Dynamics of Large-Scale Plastic Deformation and the Necking Instability in Amorphous Solids
Total Page:16
File Type:pdf, Size:1020Kb
PHYSICAL REVIEW LETTERS week ending VOLUME 90, NUMBER 4 31 JANUARY 2003 Dynamics of Large-Scale Plastic Deformation and the Necking Instability in Amorphous Solids L. O. Eastgate Laboratory of Atomic and Solid State Physics, Cornell University, Ithaca, New York 14853 J. S. Langer and L. Pechenik Department of Physics, University of California, Santa Barbara, California 93106 (Received 19 June 2002; published 31 January 2003) We use the shear transformation zone (STZ) theory of dynamic plasticity to study the necking instability in a two-dimensional strip of amorphous solid. Our Eulerian description of large-scale deformation allows us to follow the instability far into the nonlinear regime. We find a strong rate dependence; the higher the applied strain rate, the further the strip extends before the onset of instability. The material hardens outside the necking region, but the description of plastic flow within the neck is distinctly different from that of conventional time-independent theories of plasticity. DOI: 10.1103/PhysRevLett.90.045506 PACS numbers: 62.20.Fe, 46.05.+b, 46.35.+z, 83.60.Df Conventional descriptions of plastic deformation in when the strip is loaded slowly. One especially important solids consist of phenomenological rules of behavior, element of our analysis is our ability to interpret flow and with qualitative distinctions between time-independent hardening in terms of the internal STZ variables. and time-dependent properties, and sharply defined yield To make this problem as simple as possible, we consi- criteria. Plasticity, however, is an intrinsically dynamic der here only strictly two-dimensional, amorphous mate- phenomenon. Practical theories of plasticity should con- rials. By ‘‘strictly,’’ we mean that elastic and plastic sist not of intricate sets of rules, but of equations of displacement rates are separately planar as in two- motion for material velocities, stress fields, and other dimensional molecular dynamics simulations. The two- variables that might characterize internal states of solids. dimensional STZ equations presented in this paper are Roughly speaking, a theory of plasticity, especially for an based on earlier work by Falk, Langer, and Pechenik amorphous solid, should resemble the Navier-Stokes [7,11,12]. We use Eulerian coordinates in which, as in equation for a fluid, with the pressure replaced by a stress fluid dynamics, the variables xi denote the current physi- tensor, and the viscous forces replaced by a constitutive cal positions of material elements. Let the system lie law relating the rate of plastic deformation to the stresses in the x1 x, x2 y plane, and write the stress tensor 1 and internal state variables. That constitutive law should in the form: ij ÿpij sij;pÿ2 kk, where p is contain phenomenological constants, analogous to the the pressure and sij is the deviatoric stress — a trace- bulk and shear viscosities, that are measurable and, in less, symmetric tensor. In analogy to fluid dynamics, principle, computable from molecular theories. Yield cri- let vi x; y; t denote the material velocity at the physi- teria, work hardening, hysteretic effects, and the like cal position x; y, and time t. Then the acceleration equa- would emerge naturally in such a formulation. tion is [13] The goal of the STZ (shear transformation zone) dvi @ij @p @sij theory of plasticity [1–7], from its inception, has been ÿ : (1) dt @x @x @x to carry out the above program. In this paper we show j i j how the STZ theory describes a special case of large- Here, is the density which, because we shall assume a scale yielding, specifically, the necking instability of a very small elastic compressibility and volume conserving strip of material subject to tensile loading. There is a large plasticity, we shall take to be a constant. The symbol d=dt literature on the necking problem. References that we denotes the material time derivative acting on a scalar or a have found particularly valuable include papers by vector field: Hutchinson and Neale [8], McMeeking and Rice [9], d @ @ and Tvergaard and Needleman [10]. Our purpose here is vk : (2) dt @t @x to explore possibilities for using the STZ theory to in- k vestigate a range of failure mechanisms in amorphous Our first main assumption is that the rate of deforma- solids, possibly including fracture. We are able to follow tion tensor can be written as the sum of linear elastic and the necking instability far into the nonlinear regime plastic contributions: where the neck appears to be approaching plastic failure 1 @v @vj D sij p while the outer regions of the strip become hardened and Dtotal i ÿ Dplast; ij 2 @x @x t 2 2K ij ij remain intact. We find that necking in the STZ theory is j i D rate dependent; the instability occurs at smaller strains (3) 045506-1 0031-9007=03=90(4)=045506(4)$20.00 2003 The American Physical Society 045506-1 PHYSICAL REVIEW LETTERS week ending VOLUME 90, NUMBER 4 31 JANUARY 2003 where is the shear modulus, K 1 = 1 ÿ is yield stress.We also have assumed that the local density of the two-dimensional inverse compressibility (or bulk STZs is always at its equilibrium value so that we do not modulus), and is the two-dimensional Poisson ratio. need to solve an extra equation of motion for that field The symbol D=Dt denotes the material time derivative (denoted by the symbol in earlier papers). acting on any tensor, say Aij: The important exception alluded to above is the pres- ence of the absolute-value bars in Eq. (7). The expression DAij @Aij @Aij vk Aik !kj ÿ !ik Akj; (4) inside the bars is proportional to the rate at which plastic Dt @t @xk work is being done on the system, a quantity which appears in the original theory as a non-negative factor and !ij is the spin: in the STZ annihilation and creation rates. A negative 1 @vi @vj value of this quantity would be unphysical. In earlier !ij ÿ : (5) 2 @xj @xi studies of spatially uniform systems, this quantity always remained positive; however, we have observed negative The plastic part of the rate-of-deformation Dplast,like ij values in the present calculations. The absolute value s , is a traceless symmetric tensor, thus the plastic defor- ij prevents such unphysical behavior and is consistent with mations are area conserving. For present purposes, we use the intent of the original theory.We emphasize, however, a simple, quasilinear form of the STZ theory in which that this term contains some of the principal assumptions plast of the STZ theory. There are other possibilities for it (see, Dij 0 qij s; ; qij s; sij ÿ ij; (6) for example, [2]); and it will be interesting to explore the and 0 is a material-specific constant. The traceless, sym- physical significance of these variations of the model. metric tensor ij is the internal state variable mentioned To understand the transition between viscoelastic and earlier. It is proportional to a director matrix that specifies viscoplastic behaviors at the yield stress, and the role the orientation of the STZs; its magnitude is a measure played by the state variable , it is easiest to look of the degree of their alignment. The equation of motion first at a uniform system under pure shear. Let sxx for ij is ÿsyy s, sxy 0, xx ÿyy , xy 0; and con- sider a situation in which s is held constant. Equations (3) D 1 ij q ÿ jq s j : (7) and (7) become Dt ij 2 km km ij "_ s ÿ ; (8) In Eq. (6), plays—very roughly—the role of the ‘‘back 0 stress’’ or ‘‘hardening’’ parameter in conventional theo- _ s ÿ 1 ÿ s; (9) ries of plasticity [14–16], a major difference being that emerges directly from a rate equation governing the where "_ is the total strain rate. At s 1, these equations population of STZs and is, in principle, a directly mea- exhibit an exchange of stability between the nonflowing surable quantity [1,17]. If the second term on the right- steady-state solution with "_ 0, s for s<1, and the hand side of Eq. (7) were missing, then would be flowing solution with "_ Þ 0, 1=s for s>1.Asex- proportional to the integrated plastic strain. This second plained in earlier publications, the steady-state system is term, however, which is produced by the creation and ‘‘jammed’’ or ‘‘hardened’’ in the direction of the applied annihilation of STZs, is a crucial element of the STZ stress for s<1; whereas, for s>1, new STZs are being theory. As we shall show briefly below, this term produces created as fast as existing ones transform, and there is a the exchange of dynamic stability between viscoelastic nonzero plastic strain rate. and viscoplastic states that replaces the conventional as- Our goal now is to see how this exchange of stability sumptions of yield surfaces and other purely phenomeno- occurs in a dynamic, spatially nonuniform situation. logical rules of behavior. Consider a rectangle with straight grips at x L t. With one important exception, Eqs. (6) and (7) con- The upper and lower surfaces, at y Y x; t, are free stitute a tensorial version of the original STZ theory boundaries. We assume symmetry about both the x and y obtained by linearizing the stress dependence of the axes so that we need to consider only the first quadrant of rate factors and rescaling. Because of the linearization, the system. On the free upper boundary, the relation these equations do not properly describe memory effects between the material velocities and the motion of the present in the full theory that are important when the surface is system is unloaded or reloaded, but this will not affect @Y @Y our results until the system reaches the necking instabil- vy x; Y; tÿvx x; Y; t : (10) ity.