Geometrically necessary dislocations in viscoplastic single crystals and bicrystals undergoing small deformations

In this study we develop a gradient theory of small-deformation single-crystal plasticity that accounts for geometrically necessary dislocations (GNDs). The resulting framework is used to discuss grain boundaries. The grains are allowed to slip along the interface, but growth phenomenona and phase transitions are neglected. The bulk theory is based on the introduction of a microforce balance for each slip system and includes a defect energy depending on a suitable measure of GNDs. The microforce balances are shown to be equivalent to nonlocal yield conditions for the individual slip systems, yield conditions that feature backstresses resulting from energy stored in dislocations. When applied to a grain boundary the theory leads to concomitant yield conditions: relative slip of the grains is activated when the shear stress reaches a suitable threshold; plastic slip in bulk at the grain boundary is activated only when the local density of GNDs reaches an assigned threshold. Consequently, in the initial stages of plastic deformation the grain boundary acts as a barrier to plastic slip, while in later stages the interface acts as a source or sink for dislocations. We obtain an exact solution for a simple problem in plane strain involving a semi-infinite compressed specimen that abuts a rigid material. We view this problem as an approximation to a situation involving a grain boundary between a grain with slip systems aligned for easy flow and a grain whose slip system alignment severely inhibits flow. The solution exhibits large slip gradients within a thin layer at the grain boundary.     

Bauschinger and size effects in thin-film plasticity due to defect-energy of geometrical necessary dislocations

The Bauschinger and size effects in the thinfilm plasticity theory arising from the defect-energy of geometrically necessary dislocations (GNDs) are analytically investigated in this paper. Firstly, this defect-energy is deduced based on the elastic interactions of coupling dislocations (or pile-ups) moving on the closed neighboring slip plane. This energy is a quadratic function of the GNDs density, and includes an elastic interaction coefficient and an energetic length scale L. By incorporating it into the work- conjugate strain gradient plasticity theory of Gurtin, an energetic stress associated with this defect energy is obtained, which just plays the role of back stress in the kinematic hardening model. Then this back-stress hardening model is used to investigate the Bauschinger and size effects in the tension problem of single crystal Al films with passivation layers. The tension stress in the film shows a reverse dependence on the film thickness h. By comparing it with discrete-dislocation simulation results, the length scale L is determined, which is just several slip plane spacing, and accords well with our physical interpretation for the defect- energy. The Bauschinger effect after unloading is analyzed by combining this back-stress hardening model with a friction model. The effects of film thickness and pre-strain on the reversed plastic strain after unloading are quantified and qualitatively compared with experiment results.     

A discrete dislocation analysis of bending

Bending of a strip in plane strain is analyzed using discrete dislocation plasticity where the dislocations are modeled as line defects in a linear elastic medium. At each stage of loading, superposition is used to represent the solution in terms of the infinite medium solution for the discrete dislocations and a complementary solution that enforces the boundary conditions, which is non-singular and obtained from a linear elastic, finite element solution. The lattice resistance to dislocation motion, dislocation nucleation and dislocation annihilation are incorporated into the formulation through a set of constitutive rules. Solutions for cases with multiple slip systems and with a single slip system are presented. The bending moment versus rotation relation and the evolution of the dislocation structure are outcomes of the boundary value problem solution. The effects of slip geometry, obstacles to dislocation motion and specimen size on the moment versus rotation response are considered. Also, the evolution of the dislocation structure is studied with emphasis on the role of geometrically necessary dislocations. The dislocation structure that develops leads to well-defined slip bands, with the slip band spacing scaling with the specimen height.     

Anti-plane shear cracks in ideally plastic crystals

Cracks in ductile single crystals are analyzed here for geometries and orientations such that two-dimensional states of anti-plane shear constitute possible deformation fields. The crystals are modelled as ideally plastic and yield according to critical resolved shear stresses on their slip systems. Restrictions on the asymptotic forms of stress and deformation fields at crack tips are established for anti-plane loading of stationary and quasistatically growing cracks, and solutions are presented for several specific orientations in f.c.c. and b.c.c. crystals. The asymptotic solutions are complemented by complete elastic-plastic solutions for stationary and growing cracks under small scale yielding, based on previous work by Rice (1967, 1984) and Freund (1979). Remarkably, the plastic zone at a stationary crack tip collapses into discrete planes of displacement and stress discontinuity emanating from the tip; plastic flow consists of concentrated shear on the displacement discontinuities. For the growing crack these same planes, if not coincident with the crack plane, constitute collapsed plastic zones in which velocity and plastic strain discontinuities occur, but across which the stresses and anti-plane displacement are fully continuous. The planes of discontinuity are in several cases coincident with crystal slip planes but it is shown that this need not be the case, e.g., for orientations in which anti-plane yielding occurs by multi-slip, or for special orientations in which the crack tip and the discontinuity planes are perpendicular to the activated slip plane.     

Lengthscale-dependent,elastically anisotropic,physically-based hcp crystal plasticity: Application to cold-dwell fatigue in Ti alloys

A crystal plasticity model for hcp materials is presented which is based on dislocation glide and pinning. Slip is assumed to occur on basal and prismatic systems, and dislocation pinning through the generation of geometrically necessary dislocations (GNDs). Elastic anisotropy and, through the coupling of GNDs with slip rate, physically-based lengthscale effects are included.     

On large-strain finite element solutions of higher-order gradient crystal plasticity

A finite-strain higher-order gradient crystal plasticity model accounting for the backstress effect originating from the existence of geometrically necessary dislocations (GNDs) is applied to plane strain finite element analysis. Different element types are tested to seek out an element formulation that is reliable and useful for solving problems involving severe plastic deformation. In the present finite element formulation, the GND density rates are chosen to be additional nodal degrees of freedom. Different orders of shape functions are employed for the interpolation of displacement rates and GND density rates. Their effects on solutions are examined in detail by considering three boundary value problems: a simple shear of a constrained layer (a film), a compression problem with loading surfaces impenetrable to dislocations, and a tension problem involving shear band formation. In all the cases, the formulation in which eight-node elements with reduced integration and four-node elements with full integration are used respectively for displacement rates and the GND density rates gives reasonable solutions. In addition to the discussion on the choice of finite elements, detailed behavior in gradient-dependent solids, such as the accumulation of GND density and the distribution of backstress on each slip system, is investigated by utilizing the reliable computational results obtained.     

Continuum dynamics of the formation,migration and dissociation of self-locked dislocation structures on parallel slip planes

In continuum models of dislocations, proper formulations of short-range elastic interactions of dislocations are crucial for capturing various types of dislocation patterns formed in crystalline materials. In this article, the continuum dynamics of straight dislocations distributed on two parallel slip planes is modelled through upscaling the underlying discrete dislocation dynamics. Two continuum velocity field quantities are introduced to facilitate the discrete-to-continuum transition. The first one is the local migration velocity of dislocation ensembles which is found fully independent of the short-range dislocation correlations. The second one is the decoupling velocity of dislocation pairs controlled by a threshold stress value, which is proposed to be the effective flow stress for single slip systems. Compared to the almost ubiquitously adopted Taylor relationship, the derived flow stress formula exhibits two features that are more consistent with the underlying discrete dislocation dynamics: (i) the flow stress increases with the in-plane component of the dislocation density only up to a certain value, hence the derived formula admits a minimum inter-dislocation distance within slip planes; (ii) the flow stress smoothly transits to zero when all dislocations become geometrically necessary dislocations. A regime under which inhomogeneities in dislocation density grow is identified, and is further validated through comparison with discrete dislocation dynamical simulation results. Based on the findings in this article and in our previous works, a general strategy for incorporating short-range dislocation correlations into continuum models of dislocations is proposed.     

An analysis of the pile-up of infinite periodic walls of edge dislocations

We analyse the equilibrium pile-up configurations of infinite periodic walls of edge dislocations which are forced against an impenetrable obstacle by a constant applied shear stress. Numerically generated density distributions exhibit two distinct regions, for each of which we provide an interpretation and an analytical prediction. Near the obstacle, the influence of neighbouring slip planes may be neglected and the classical solution for a single slip plane applies. At a larger distance a linear decay is obtained. The characteristic length scales of the two parts of the pile-up are shown to depend differently on the parameters of the problem.     

Non-local crystal plasticity model with intrinsic SSD and GND effects

A strain gradient-dependent crystal plasticity approach is presented to model the constitutive behaviour of polycrystal FCC metals under large plastic deformation. In order to be capable of predicting scale dependence, the heterogeneous deformation-induced evolution and distribution of geometrically necessary dislocations (GNDs) are incorporated into the phenomenological continuum theory of crystal plasticity. Consequently, the resulting boundary value problem accommodates, in addition to the ordinary stress equilibrium condition, a condition which sets the additional nodal degrees of freedom, the edge and screw GND densities, proportional (in a weak sense) to the gradients of crystalline slip. Next to this direct coupling between microstructural dislocation evolutions and macroscopic gradients of plastic slip, another characteristic of the presented crystal plasticity model is the incorporation of the GND-effect, which leads to an essentially different constitutive behaviour than the statistically stored dislocation (SSD) densities. The GNDs, by their geometrical nature of locally similar signs, are expected to influence the plastic flow through a non-local back-stress measure, counteracting the resolved shear stress on the slip systems in the undeformed situation and providing a kinematic hardening contribution. Furthermore, the interactions between both SSD and GND densities are subject to the formation of slip system obstacle densities and accompanying hardening, accountable for slip resistance. As an example problem and without loss of generality, the model is applied to predict the formation of boundary layers and the accompanying size effect of a constrained strip under simple shear deformation, for symmetric double-slip conditions.     

Distribution based model for the grain boundaries in polycrystalline plasticity

This contribution focuses on the development of constitutive models for the grain boundary region between two crystals, relying on the dislocation based polycrystalline model documented in (Evers, L.P., Parks, D.M., Brekelmans, W.A.M., Geers, M.G.D., 2002. Crystal plasticity model with enhanced hardening by geometrically necessary dislocation accumulation. J. Mech. Phys. Solids 50, 2403–2424; Evers, L.P., Brekelmans, W.A.M., Geers, M.G.D., 2004a. Non-local crystal plasticity model with intrinsic SSD and GND effects. J. Mech. Phys. Solids 52, 2379–2401; Evers, L.P., Brekelmans, W.A.M., Geers, M.G.D., 2004b. Scale dependent crystal plasticity framework with dislocation density and grain boundary effects. Int. J. Solids Struct. 41, 5209–5230). The grain boundary is first viewed as a geometrical surface endowed with its own fields, which are treated here as distributions from a mathematical point of view. Regular and singular dislocation tensors are introduced, defining the grain equilibrium, both in the grain core and at the boundary of both grains. Balance equations for the grain core and grain boundary are derived, that involve the dislocation density distribution tensor, in both its regular and singular contributions. The driving force for the motion of the geometrically necessary dislocations is identified from the pull-back to the lattice configuration of the quasi-static balance of momentum, that reveals the duality between the stress and the curl of the elastic gradient. Criteria that govern the flow of mobile geometrically necessary dislocations (GNDs) through the grain boundary are next elaborated on these bases. Specifically, the sign of the projection of a lattice microtraction on the glide velocity defines a necessary condition for the transmission of incoming GNDs, thereby rendering the set of active slip systems for the glide of outgoing dislocations. Viewing the grain boundary as adjacent bands in each grain with a constant GND density in each, the driving force for the grain boundary slip is further expressed in terms of the GND densities and the differently oriented slip systems in each grain. A semi-analytical solution is developed in the case of symmetrical slip in a bicrystal under plane strain conditions. It is shown that the transmission of plastic slip occurs when the angle made by the slip direction relative to the grain boundary normal is less than a critical value, depending on the ratio of the GND densities and the orientation of the transmitted dislocations.     

颗粒增强复合材料压缩行为的位错动力学模拟

颗粒增强铜基复合材料因具有极高的强度和弹性模量, 优异的导电、导热性能和抗磨损能力, 被广泛应用于航天航空、汽车、电子工业等领域. 第二相强化是其主要的强化方式, 其通过合金中弥散的微粒阻碍位错运动, 可有效提高金属材料的力学性能, 提高其服役安全. 针对该问题本文采用三维离散位错动力学(three-dimensional discrete dislocation dynamics, 3D-DDD)方法, 对微尺度颗粒增强铜基复合材料进行了微柱压缩模拟, 分析了位错与第二相颗粒交互作用对材料力学响应的影响, 揭示第二相颗粒强化的微观机理. 本研究将第二相颗粒视为位错不可穿透的球形微粒, 采用位错绕过机制模拟颗粒与位错的交互作用过程. 通过调控滑移面相对于第二相颗粒中心的距离发现: 屈服应力和应变硬化率均随距离的增大而减小. 研究也发现Schmid因子越高的滑移系, 屈服应力越低, 后续应变硬化率越低. 多位错与颗粒交互作用的模拟发现, 同一滑移面中位错间的反应和不同滑移系中位错的交互作用可能是导致屈服应力和应变硬化率降低的关键.      

Plasticity size effects in voided crystals

The shear and equi-biaxial straining responses of periodic voided single crystals are analysed using discrete dislocation plasticity and a continuum strain gradient crystal plasticity theory. In the discrete dislocation formulation, the dislocations are all of edge character and are modelled as line singularities in an elastic material. The lattice resistance to dislocation motion, dislocation nucleation, dislocation interaction with obstacles and annihilation are incorporated through a set of constitutive rules. Over the range of length scales investigated, both the discrete dislocation and strain gradient plasticity formulations predict a negligible size effect under shear loading. By contrast, under equi-biaxial loading both plasticity formulations predict a strong size dependence with the flow strength approximately scaling inversely with the void spacing. Excellent agreement is obtained between predictions of the two formulations for all crystal types and void volume fractions considered when the material length scale in the non-local plasticity model is chosen to be (about 10 times the slip plane spacing in the discrete dislocation models).     

Polychromatic microdiffraction analysis of defect self-organization in shock deformed single crystals

A spatially resolved X-ray diffraction method – with a submicron 3D resolution together with SEM and OIM analysis are applied to understand the arrangements of voids, geometrically necessary dislocations and strain gradient distributions in samples of Al (1 2 3) and Cu (0 0 1) single crystals shocked to incipient spallation fracture. We describe how geometrically necessary dislocations and the effective strain gradient alter white beam Laue patterns of the shocked materials. Several distinct structural zones are observed at different depths under the impact surface. The density of geometrically necessary dislocations (GNDs) is extremely high near the impact and back surface of the shock recovered crystals. The spall region is characterized by a large density of mesoscale voids and GNDs. The spall region is separated from the impact and back surfaces by compressed regions with high total dislocation density but lower GNDs density. Self-organization of shear bands is observed in the shock recovered Cu single crystal.     

Dislocation pile-ups in bicrystals within continuum dislocation theory

Within continuum dislocation theory the plastic deformation of bicrystals under a mixed deformation of plane constrained uniaxial extension and shear is investigated with regard to the nucleation of dislocations and the dislocation pile-up near the phase boundaries of a model bicrystal with one active slip system within each single crystal. For plane uniaxial extension, we present a closed-form analytical solution for the evolution of the plastic distortion and of the dislocation network in the case of symmetric slip planes (i.e. for twins), which exhibits an energetic as well as a dissipative threshold for the dislocation nucleation. The general solution for non-symmetric slip systems is obtained numerically. For a combined deformation of extension and shear, we analyze the possibility of linearly superposing results obtained for both loading cases independently. All solutions presented in this paper also display the Bauschinger effect of translational work hardening and a size effect typical to problems of crystal plasticity.     

Higher-order stress and grain size effects due to self-energy of geometrically necessary dislocations

The higher-order stress work-conjugate to slip gradient in single crystals at small strains is derived based on the self-energy of geometrically necessary dislocations (GNDs). It is shown that this higher-order stress changes stepwise as a function of in-plane slip gradient and therefore significantly influences the onset of initial yielding in polycrystals. The higher-order stress based on the self-energy of GNDs is then incorporated into the strain gradient plasticity theory of Gurtin [2002. A gradient theory of single-crystal viscoplasticity that accounts for geometrically necessary dislocations. J. Mech. Phys. Solids 50, 5-32] and applied to single-slip-oriented 2D and 3D model crystal grains of size D. It is thus found that the self-energy of GNDs gives a D^-1-dependent term for the averaged resolved shear stress in such a model grain under yielding. Using published experimental data for several polycrystalline metals, it is demonstrated that the D^-1-dependent term successfully explains the grain size dependence of initial yield stress and the dislocation cell size dependence of flow stress in the submicron to several-micron range of grain and cell sizes.     

Screw dislocation pile-ups in two phase materials of finite rigidity

The arrangement of discrete screw dislocations piled-up under the action of a uniform applied stress against the welded interface between different elastically isotropic half-spaces has been determined by representing the pile-up as a continuous distribution of infinitesimal dislocations. The dislocation slip plane is inclined at an arbitrary angle $\text{1}\text{2}\text{aπ}$ to the normal to the interface, assuming a to be a rational number. The singular integral equation expressing the condition for static equilibrium of the dislocations under a constant applied stress is solved by a method based on the Wiener-Hoph technique with the Mellin transform, and from this solution the mean density of dislocations and the stress field of the pile-up are determined.     

Dislocation behaviours ahead of crack tip

In this paper, a unified mechanics model for dislocation nucleation, emission and dislocation free zone is proposed based on the Peierls framework. Three regions are identified ahead of the crack tip. The emitted dislocations within the plastic zone in the form of an inverse pile up are treated as discrete elastic edge dislocations. Between that zone and the cohesive zone immediately ahead of the crack tip, there is a dislocation free zone. With the stress field and the dislocation density field in the cohesive zone, respectively, expressed in the first and second Chebyshev polynomial series, and the opening and slip displacements in trigonometric series, a set of nonlinear governing equations are obtained which take into account for the interaction between the emitted dislocations and cohesive zone and the nonlinear interaction between sliding displacement and the opening displacement. After discretization, the governing equations are transformed into a set nonlinear algebraic equations which are solved with Newton-Raphson Method. The results of calculation for pure shearing and combined tension and shear loading after dislocation emission are given in detail. Finally, the process of dislocation nucleation and emission on a pair of symmetric slip planes of angle α with respect to the crack plane under pure mode I load is analysed. The equilibrium positions and the number of emitted dislocation are determined. Several possible competition behaviors of dislocation emission vs cleavage are revealed.