Dislocation microstructures and strain-gradient plasticity with one active slip plane

We study dislocation networks in the plane using the vectorial phase-field model introduced by Ortiz and coworkers, in the limit of small lattice spacing. We show that, in a scaling regime where the total length of the dislocations is large, the phase field model reduces to a simpler model of the strain-gradient type. The limiting model contains a term describing the three-dimensional elastic energy and a strain-gradient term describing the energy of the geometrically necessary dislocations, characterized by the tangential gradient of the slip. The energy density appearing in the strain-gradient term is determined by the solution of a cell problem, which depends on the line tension energy of dislocations. In the case of cubic crystals with isotropic elasticity our model shows that complex microstructures may form in which dislocations with different Burgers vector and orientation react with each other to reduce the total self-energy.     

Scaling of dislocation-based strain-gradient plasticity

By taking into account the dislocations that are geometrically necessary for producing a curvature or twist of the atomic lattice in crystals, Gao et al. recently developed a theory of strain-gradient plasticity on the micrometer scale and showed that it agrees relatively well with the tests of hardness, torsion and bending of copper on the micrometer scale. This paper subjects this theory to an asymptotic scaling analysis. It is shown that the small-size asymptotic limit of this theory exhibits (1) an unusually strong size effect in which the corresponding nominal stresses in geometrically similar structures of different sizes D vary as D^−5/2, and (2) an asymptotic approach to a load-deflection diagram whose tangent stiffness gradually increases, starting with an infinitely small initial stiffness at infinitely small stress. Although this peculiar small-size asymptotic behavior might not be attainable within the practical applicability range of a continuum theory, it renders questionable any efforts to construct approximations of an asymptotic matching character, with a two-sided asymptotic support, which have previously been proven effective for quasibrittle materials such as concrete, rock, ice and fiber composites. A possible simple modification of the existing theory with respect to the small-size asymptotic properties is suggested. However, the questions of experimental justification of such a modification and its compatibility with the dislocation theory will require further study. The small-size asymptotic properties of other strain gradient theories of plasticity have not been analyzed, except for those of the previous Fleck-Hutchinson theory, which have been found reasonable.     

An alternative treatment of phenomenological higher-order strain-gradient plasticity theory

Phenomenological higher-order strain-gradient plasticity is here presented through a formulation inspired by previous work for strain-gradient crystal plasticity. A physical interpretation of the phenomenological yield condition that involves an effect of second gradient of the equivalent plastic strain is discussed, applying a dislocation theory-based consideration. Then, a differential equation for the equivalent plastic strain-gradient is introduced as an additional governing equation. Its weak form makes it possible to deduce and impose extra boundary conditions for the equivalent plastic strain. A connection between the present treatment and strain-gradient theories based on an extended virtual work principle is discussed. Furthermore, a numerical implementation and analysis of constrained simple shear of a thin strip are presented.     

On boundary conditions and plastic strain-gradient discontinuity in lower-order gradient plasticity

Through linearized analysis and computation, we show that lower-order gradient plasticity is compatible with boundary conditions, thus expanding its predictive capability. A physically motivated gradient modification of the conventional Voce hardening law is shown to lead to a convective stabilizing effect in 1-D, rate-independent plasticity. The partial differential equation is genuinely nonlinear and does not arise as a conservation law, thus making the task of inferring plausible boundary conditions a delicate matter. Implications of wave-type behavior in rate-independent plastic response (under conditions of static equilibrium) are analyzed with a discussion of an appropriate numerical algorithm. Example problems are solved numerically, showing the robustness and simplicity of physically motivated lower-order gradient plasticity. The 3-D case and rate-dependent constitutive assumptions are also discussed.     

Effect of defect energy on strain-gradient predictions of confined single-crystal plasticity

Gurtin recently proposed a strain-gradient theory for crystal plasticity in which the gradient effect originates from a defect energy that characterizes energy storage due to the presence of a net Burgers vector. Here we consider a number of different possibilities for this energy: specifically, working within a simple two-dimensional framework, we compare predictions of the theory with results of discrete-dislocation simulations of stress relaxation in thin films. Our objective is to investigate which specific defect energies are capable of capturing the size-dependent response of such systems for different crystal orientations.     

Continuum approximation of the Peach-Koehler force on dislocations in a slip plane

We derive a continuum model for the Peach-Koehler force on dislocations in a slip plane. To represent the dislocations, we use the disregistry across the slip plane, whose gradient gives the density and direction of the dislocations. The continuum model is derived rigorously from the Peach-Koehler force on dislocations in a region that contains many dislocations. The resulting continuum model can be written as the variation of an elastic energy that consists of the contribution from the long-range elastic interaction of dislocations and a correction due to the line tension effect.     

Randomness and slip avalanches in gradient plasticity

While localization of deformation at macroscopic scales has been documented and carefully characterized long ago, it is only recently that systematic experimental investigations have demonstrated that plastic flow of crystalline solids on mesoscopic scales proceeds in a strongly heterogenous and intermittent manner. In fact, deformation is characterized by intermittent bursts (‘slip avalanches’) the sizes of which obey power-law statistics. In the spatial domain, these avalanches produce characteristic deformation patterns in the form of slip lines and slip bands. Unlike to the case of macroscopic localization where gradient plasticity can capture the width and spacing of shear bands in the softening regime of the stress–strain graph, this type of mesoscopically jerky like localized plastic flow is observed in spite of a globally convex stress–strain relationship and may not be captured by standard deterministic continuum modelling. We thus propose a generalized constitutive model which includes both second-order strain gradients and randomness in the local stress–strain relationship. These features are related to the internal stresses which govern dislocation motion on microscopic scales. It is shown that the model can successfully describe experimental observations on slip avalanches as well as the associated surface morphology characteristics.     

Calibration,validation, and verification including uncertainty of a physically motivated internal state variable plasticity and damage model

In this paper, we illustrate a formal calibration, validation, and verification process that includes uncertainty in an internal state variable plasticity-damage model that is implemented in a finite element code. The physically motivated continuum model characterizes damage evolution by incorporating material uncertainty due to microstructural spatial clustering. The uncertainty analysis is performed by introducing material variation through model validation and verification. The effect of variability in microstructural clustering and boundary conditions on the sensitivities and uncertainty of the plasticity-damage evolution for the 7075 aluminum alloy is characterized. The results show the potential of this methodology in the evaluation of material response uncertainty due to microstructure spatial clustering and its effect on damage evolution. For damage evolution, we have shown that the initial isotropic damage evolved into an anisotropic form as the deformation increased which is consistent with experimentally observed behavior for 7075 aluminum alloy in literature. Also, the sensitivities were found to be consistent with the physics of damage progression for this particular type of material. Through the sensitivity analysis, the initial defect size and number density of cracked particles are important at the beginning of deformation. As damage evolves, more voids are nucleated and grow and the sensitivity analysis illustrates this as well. Then, voids combine with each other and coalescence becomes the main driver, which is also confirmed by the sensitivity analysis. This work also shows that the microstructurally based damage evolution equations provide an accurate representation of the damage progression due to large intermetallic particles. Finally, we illustrate that the initial variation in the microstructure clustering can lead to about ±7.0%, ±8.1%, and ±9.75% variation in the elongation to failure strain for torsion, tensile, and compressive loading, respectively.     

A one-dimensional theory of strain-gradient plasticity: Formulation, analysis, numerical results

This study develops a one-dimensional theory of strain-gradient plasticity based on: (i) a system of microstresses consistent with a microforce balance; (ii) a mechanical version of the second law that includes, via microstresses, work performed during viscoplastic flow; (iii) a constitutive theory that allows

•

the free-energy to depend on the gradient of the plastic strain, and

•

the microstresses to depend on the gradient of the plastic strain-rate.

The constitutive equations, whose rate-dependence is of power-law form, are endowed with energetic and dissipative gradient length-scales L and l, respectively, and allow for a gradient-dependent generalization of standard internal-variable hardening. The microforce balance when augmented by the constitutive relations for the microstresses results in a nonlocal flow rule in the form of a partial differential equation for the plastic strain. Typical macroscopic boundary conditions are supplemented by nonstandard microscopic boundary conditions associated with flow, and properties of the resulting boundary-value problem are studied both analytically and numerically. The resulting solutions are shown to exhibit three distinct physical phenomena:

(i)

standard (isotropic) internal-variable hardening;

(ii)

energetic hardening, with concomitant back stress, associated with plastic-strain gradients and resulting in boundary layer effects;

(iii)

dissipative strengthening associated with plastic strain-rate gradients and resulting in a size-dependent increase in yield strength.

     

The role of dissipation and defect energy in variational formulations of problems in strain-gradient plasticity. Part 1: polycrystalline plasticity

A general set of flow laws and associated variational formulations are constructed for small-deformation rate-independent problems in strain-gradient plasticity. The framework is based on the thermodynamically consistent theory due to Gurtin and Anand (J Mech Phys Solids 53:1624–1649, 2005), and includes as variables a set of microstresses which have both energetic and dissipative components. The flow law is of associative type. It is expressed as a normality law with respect to a convex but otherwise arbitrary yield function, or equivalently in terms of the corresponding dissipation function. Two cases studied are, first, an extension of the classical Hill-Mises or J ₂ flow law and second, a form written as a linear sum of the magnitudes of the plastic strain and strain gradient. This latter form is motivated by work of Evans and Hutchinson (Acta Mater 57:1675–1688, 2009) and Nix and Gao (J Mech Phys Solids 46:411–425, 1998), who show that it leads to superior correspondence with experimental results, at least for particular classes of problems. The corresponding yield function is obtained by a duality argument. The variational problem is based on the flow rule expressed in terms of the dissipation function, and the problem is formulated as a variational inequality in the displacement, plastic strain, and hardening parameter. Dissipative components of the microstresses, which are indeterminate, are absent from the formulation. Existence and uniqueness of solutions are investigated for the generalized Hill-Mises and linear-sum dissipation functions, and for various combinations of defect energy. The conditions for well-posedness of the problem depend critically on the choice of dissipation function, and on the presence or otherwise of a defect energy in the plastic strain or plastic strain gradient, and of internal-variable hardening.     

The role of dissipation and defect energy in variational formulations of problems in strain-gradient plasticity. Part 2: single-crystal plasticity

Variational formulations are constructed for rate-independent problems in small-deformation single-crystal strain-gradient plasticity. The framework, based on that of Gurtin (J Mech Phys Solids 50: 5–32, 2002), makes use of the flow rule expressed in terms of the dissipation function. Provision is made for energetic and dissipative microstresses. Both recoverable and non-recoverable defect energies are incorporated into the variational framework. The recoverable energies include those that depend smoothly on the slip gradients, the Burgers tensor, or on the dislocation densities (Gurtin et al. J Mech Phys Solids 55:1853–1878, 2007), as well as an energy proposed by Ohno and Okumura (J Mech Phys Solids 55:1879–1898, 2007), which leads to excellent agreement with experimental results, and which is positively homogeneous and therefore not differentiable at zero slip gradient. Furthermore, the variational formulation accommodates a non-recoverable energy due to Ohno et al. (Int J Mod Phys B 22:5937–5942, 2008), which is also positively homogeneous, and a function of the accumulated dislocation density. Conditions for the existence and uniqueness of solutions are established for the various examples of defect energy, with or without the presence of hardening or slip resistance.     

A crystal plasticity model that incorporates stresses and strains due to slip gradients

This work is concerned with incorporating the kinematic and stress effects of excess dislocations in a constitutive model for the elastoplastic behavior of crystalline materials. The foundation of the model is a three term multiplicative decomposition of the deformation gradient in which the two classical terms of plastic and elastic deformation are included along with an additional term for long range strain due to the collective effects of excess dislocations. The long range strain is obtained from an assumed density of Volterra edge dislocations and is directly related to gradients in slip. A new material parameter emerges which is the size the region about a continuum point that contributes to long range strains.Using Hookean elasticity, the stress at a point is linearly related to the sum of the elastic plus the long range strain fields. However, the driving force for slip is postulated to be due only to the elastic stress so that the long range stress is a back stress in the constitutive relationship for plastic deformation. A consistent balance of the total deformation rate with the three proposed mechanisms of deformation leads to a set of differential equations that can be solved for the elastic stress, rotation and pressure which then implicitly defines the material state and equilibrium stress. Results from the simulation of a tapered tensile specimen demonstrate that the constitutive model exhibits isotropic and kinematic type hardening effects as well as changes in the pattern of plastic deformation and necking when compared to a material without slip gradient effects.     

Damped vibration of a graphene sheet using a higher-order nonlocal strain-gradient Kirchhoff plate model

A higher-order nonlocal strain-gradient model is presented for the damped vibration analysis of single-layer graphene sheets (SLGSs) in hygrothermal environment. Based on Kirchhoff plate theory in conjunction with a higher-order (bi-Helmholtz) nonlocal strain gradient theory, the equations of motion are obtained using Hamilton's principle. The higher-order nonlocal strain gradient theory has lower- and higher-order nonlocal parameters and a material characteristic parameter. The presented model can reasonably interpret the softening effects of the SLGS, and indicates a reasonably good match with the experimental flexural frequencies. Finally, the roles of viscous and structural damping coefficients, small-scale parameters, hygrothermal environment and elastic foundation on the vibrational responses of SLGSs are studied in detail.     

Macro slip theory of plasticity for polycrystalline solids

A macro slip theory is presented in this paper. Four independent slip systems are proposed for polycrystalline solids. Each slip system consists of a slip plane which lies on a face of the octahedron in stress space and a slip direction which is coincident with shear stress acting on the same face of the octahedron. It is proved that for proportional loading, present results are identical with the classical flow theory of plasticity. For nonproportional loading, the macro slip theory shows good predicting ability. The calculated results are in good agreement with the experimental data. The project supported by Chinese Academy of Science     

Relaxation in crystal plasticity with three active slip systems

We study a variational model for finite crystal plasticity in the limit of rigid elasticity. We focus on the case of three distinct slip systems whose slip directions lie in one plane and are rotated by 120° with respect to each other, with linear self-hardening and infinite latent hardening, in the sense that each material point has to deform in single slip. Under these conditions, plastic deformation is accompanied by the formation of fine-scale structures, in which activity along the different slip systems localizes in different areas. The quasiconvex envelope of the energy density, which describes the macroscopic material behavior, is determined in a regime from small up to intermediate strains, and upper and lower bounds are provided for large strains. Finally sufficient conditions are given under which the lamination convex envelope of an extended-valued energy density is an upper bound for its quasiconvex envelope.     

The structure and energetics of, and the plasticity caused by, Eshelby dislocations

The structure of coaxial, or Eshelby, dislocations are computed using isotropic elasticity for arrays of up to 500 dislocations. The energies of these arrays are determined in order to predict the lowest energy configuration and multiple meta-stable configurations are often found. The energy from these elasticity predictions shows good agreement with molecular statics simulations of aluminum. From these simulations, the torque-twist curves are predicted and compared with molecular dynamics simulations.