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.     

Role of grain boundary energetics on the maximum strength of nanocrystalline Nickel

Numerical simulations are used to investigate the competing grain boundary and dislocation mediated deformation mechanisms in nanocrystalline Ni with grain sizes in the range 4-32 nm. We present a 3D phase field model that tracks the evolution of individual dislocations and grain boundaries. Our model shows that the transition from Hall-Petch to inverse Hall-Petch as the grain size is reduced cannot be characterized only by the grain size, but it is also affected by the grain boundary energetics. We find that the grain size corresponding to the maximum yield stress (the transition from Hall-Petch strengthening with decreasing grain size to inverse Hall-Petch) decreases with increasing grain boundary energy. Interestingly, we find that for grain boundaries with high cohesive energy the Hall-Petch maximum is not observed for grains in the range 4-32 nm.     

A screw dislocation nucleated from a triple junction

Dislocation mechanism is still considered to be an effective way in plastic deformation of nanocrystalline materials. The possible nucleation of a screw dislocation from the triple junction is explored. The nucleation of a screw dislocation is found to be rather difficult from the triple junction with well-bonded grain boundaries due to the small singularity (0 to −1/3) introduced by the elasticity anisotropy. For the triple junction with a free sliding grain boundary, the stress singularity is higher than −0.5 and dislocations can be spontaneously nucleated for grain sizes less than 100 nm.     

A generalized self-consistent polycrystal model for the yield strength of nanocrystalline materials

Inspired by recent molecular dynamic simulations of nanocrystalline solids, a generalized self-consistent polycrystal model is proposed to study the transition of yield strength of polycrystalline metals as the grain size decreases from the traditional coarse grain to the nanometer scale. These atomic simulations revealed that a significant portion of atoms resides in the grain boundaries and the plastic flow of the grain-boundary region is responsible for the unique characteristics displayed by such materials. The proposed model takes each oriented grain and its immediate grain boundary to form a pair, which in turn is embedded in the infinite effective medium with a property representing the orientational average of all these pairs. We make use of the linear comparison composite to determine the nonlinear behavior of the nanocrystalline polycrystal through the concept of secant moduli. To this end an auxiliary problem of Christensen and Lo (J. Mech. Phys. Solids 27 (1979) 315) superimposed on the eigenstrain field of Luo and Weng (Mech. Mater. 6 (1987) 347) is first considered, and then the nonlinear elastoplastic polycrystal problem is addressed. The plastic flow of each grain is calculated from its crystallographic slips, but the plastic behavior of the grain-boundary phase is modeled as that of an amorphous material. The calculated yield stress for Cu is found to follow the classic Hall-Petch relation initially, but as the gain size decreases it begins to depart from it. The yield strength eventually attains a maximum at a critical grain size and then the Hall-Petch slope turns negative in the nano-range. It is also found that, when the Hall-Petch relation is observed, the plastic behavior of the polycrystal is governed by crystallographic slips in the grains, but when the slope is negative it is governed by the grain boundaries. During the transition both grains and grain boundaries contribute competitively.     

Intra-granular plastic slip heterogeneities: Discrete vs. Mean Field approaches

In this paper, we derive the mechanical fields (internal stresses, elastic energy) arising from the presence of an inelastic distortion field representing a typical intra-granular “microstructure” as the one observed during the plastification of metallic polycrystals. This “microstructure” is due to the formation of discrete intra-granular plastic slip heterogeneities characterized by at least two internal lengths: the first one is the individual grain size which represents a stochastic parameter inherent to the processing route (prior working, annealing), and, the second one is the spatial distance between active slip lines or slip bands associated with inhomogeneous plastic slip in the interior of grains. These internal lengths can be observed and measured using conventional experimental techniques (EBSD, TEM, AFM). The micro-mechanical modeling of the mechanical fields associated with plastic slip events inside grains is performed with two different assumptions. The first one is based on the well-known Eshelby’s problem of plastic inclusion where only the grain diameter is considered as internal length scale. This classical method considers homogeneous plastic distortion in the grain and leads to a uniform and grain size independent total strain field in the grain. The second method accounts for a non-uniform plastic distortion in the grain characterized by its discrete nature and the two aforementioned internal lengths. Both methods consider grains as spherical inclusions with a given diameter embedded in a homogeneous medium. For the second method, plastic slip is constrained by grain boundaries seen as impenetrable obstacles to dislocations. Thus, plastic strain is embodied by distributions of discrete circular glide loops. After writing the field equations and the free energy of the medium, a micro-mechanical formulation based on the Fourier transform method is developed. It is then found that in contrast with the mean-field approach, the internal stress fields as well as the elastic energy corresponding to different dislocation configurations depend on internal lengths associated to the deformed medium. Different possible configurations associated with intra-granular plastic flow due to circular glide dislocation loops are analyzed. Finally, the results are discussed with respect to the grain size dependence of the flow strength and the Bauschinger effect for plastically deforming polycrystals and perspectives to develop new micro–macro transition schemes accounting for internal length scales are sketched out.     

Grain boundary interface mechanics in strain gradient crystal plasticity

Interactions between dislocations and grain boundaries play an important role in the plastic deformation of polycrystalline metals. Capturing accurately the behaviour of these internal interfaces is particularly important for applications where the relative grain boundary fraction is significant, such as ultra fine-grained metals, thin films and micro-devices. Incorporating these micro-scale interactions (which are sensitive to a number of dislocation, interface and crystallographic parameters) within a macro-scale crystal plasticity model poses a challenge. The innovative features in the present paper include (i) the formulation of a thermodynamically consistent grain boundary interface model within a microstructurally motivated strain gradient crystal plasticity framework, (ii) the presence of intra-grain slip system coupling through a microstructurally derived internal stress, (iii) the incorporation of inter-grain slip system coupling via an interface energy accounting for both the magnitude and direction of contributions to the residual defect from all slip systems in the two neighbouring grains, and (iv) the numerical implementation of the grain boundary model to directly investigate the influence of the interface constitutive parameters on plastic deformation. The model problem of a bicrystal deforming in plane strain is analysed. The influence of dissipative and energetic interface hardening, grain misorientation, asymmetry in the grain orientations and the grain size are systematically investigated. In each case, the crystal response is compared with reference calculations with grain boundaries that are either ‘microhard’ (impenetrable to dislocations) or ‘microfree’ (an infinite dislocation sink).     

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.     

A study of microstructural length scale effects on the behaviour of FCC polycrystals using strain gradient concepts

Grain size is a critically important aspect of polycrystalline materials and experimental observations on Cu and Al polycrystals have shown that a Hall–Petch-type phenomenon does exist at the onset of plastic deformation. In this work, a parametric study is conducted to investigate the effect of microstructural and deformation-related length scales on the behaviour of such FCC polycrystals. It relies on a recently proposed non-local dislocation-mechanics based crystallographic theory to describe the evolution of dislocation mean spacings within each grain, and on finite element techniques to incorporate explicitly grain interaction effects. Polycrystals are modeled as representative volume elements (RVEs) containing up to 64 randomly oriented grains. Predictions obtained from RVEs of Cu polycrystals with different grain sizes are shown to be consistent with experimental data. Furthermore, mesh sensitivity studies revealed that, when there is a predominance of geometrically necessary dislocations relative to statistically stored dislocations, the polycrystal response becomes increasingly mesh sensitive. This was found to occur especially during the early stages of deformation in polycrystals with small grains.     

Assessment of plastic heterogeneity in grain interaction models using crystal plasticity finite element method

Micromechanical models aimed at simulating deformation textures and resulting plastic anisotropy need to incorporate local plastic strain heterogeneities arising from grain interactions for better predictions. The ALAMEL model [Van Houtte, P., Li, S., Seefeldt, M., Delannay, L. 2005. Deformation texture prediction: from the Taylor model to the advanced Lamel model. Int. J. Plasticity 21, 589–624], is one of the models in which the heterogeneous nature of plastic deformation in metals is introduced by accounting for the influence of a grain boundary on the cooperative deformation of adjacent grains. This is achieved by assuming that neighbouring grains undergo heterogeneous shear rates parallel to the grain boundary. The present article focuses on understanding the plastic deformation fields near the grain boundaries and the influence of grain interaction on intra-grain deformations. Crystal Plasticity Finite Element Method (CPFEM) is employed on a periodic unit cell consisting of four grains discretised into a large number of elements. A refined study of the local variation of strain rates, both along and perpendicular to the grain boundaries permits an assessment of the assumptions made in the ALAMEL model. It is shown that the ALAMEL model imbibes the nature of plastic deformation at the grain boundaries very well. However, near triple junctions, the influence of a third grain induces severe oscillations of the stress tensor, reflecting a singularity. According to CPFEM, such singularity can lead to grain subdivision by the formation of new boundaries originating at the triple junction.     

Crystal plasticity model with enhanced hardening by geometrically necessary dislocation accumulation

A strain gradient dependent crystal plasticity approach is used to model the constitutive behaviour of polycrystal FCC metals under large plastic deformation. Material points are considered as aggregates of grains, subdivided into several fictitious grain fractions: a single crystal volume element stands for the grain interior whereas grain boundaries are represented by bi-crystal volume elements, each having the crystallographic lattice orientations of its adjacent crystals. A relaxed Taylor-like interaction law is used for the transition from the local to the global scale. It is relaxed with respect to the bi-crystals, providing compatibility and stress equilibrium at their internal interface. During loading, the bi-crystal boundaries deform dissimilar to the associated grain interior. Arising from this heterogeneity, a geometrically necessary dislocation (GND) density can be computed, which is required to restore compatibility of the crystallographic lattice. This effect provides a physically based method to account for the additional hardening as introduced by the GNDs, the magnitude of which is related to the grain size. Hence, a scale-dependent response is obtained, for which the numerical simulations predict a mechanical behaviour corresponding to the Hall-Petch effect. Compared to a full-scale finite element model reported in the literature, the present polycrystalline crystal plasticity model is of equal quality yet much more efficient from a computational point of view for simulating uniaxial tension experiments with various grain sizes.     

Plastic deformation of freestanding thin films: Experiments and modeling

Experimental measurements and computational results for the evolution of plastic deformation in freestanding thin films are compared. In the experiments, the stress-strain response of two sets of Cu films is determined in the plane-strain bulge test. One set of samples consists of electroplated Cu films, while the other set is sputter-deposited. Unpassivated films, films passivated on one side and films passivated on both sides are considered. The calculations are carried out within a two-dimensional plane strain framework with the dislocations modeled as line singularities in an isotropic elastic solid. The film is modeled by a unit cell consisting of eight grains, each of which has three slip systems. The film is initially free of dislocations which then nucleate from a specified distribution of Frank-Read sources. The grain boundaries and any film-passivation layer interfaces are taken to be impenetrable to dislocations. Both the experiments and the computations show: (i) a flow strength for the passivated films that is greater than for the unpassivated films and (ii) hysteresis and a Bauschinger effect that increases with increasing pre-strain for passivated films, while for unpassivated films hysteresis and a Bauschinger effect are small or absent. Furthermore, the experimental measurements and computational results for the 0.2% offset yield strength stress, and the evolution of hysteresis and of the Bauschinger effect are in good quantitative agreement.     

Size effects and dislocation patterning in two-dimensional bending

We perform atomistic Monte Carlo simulations of bending a Lennard-Jones single crystal in two dimensions. Dislocations nucleate only at the free surface as there are no sources in the interior of the sample. When dislocations reach sufficient density, they spontaneously coalesce to nucleate grain boundaries, and the resulting microstructure depends strongly on the initial crystal orientation of the sample. In initial yield, we find a reverse size effect, in which larger samples show a higher scaled bending moment than smaller samples for a given strain and strain rate. This effect is associated with source-limited plasticity and high strain rate relative to dislocation mobility, and the size effect in initial yield disappears when we scale the data to account for strain rate effects. Once dislocations coalesce to form grain boundaries, the size effect reverses and we find that smaller crystals support a higher scaled bending moment than larger crystals. This finding is in qualitative agreement with experimental results. Finally, we observe an instability at the compressed crystal surface that suggests a novel mechanism for the formation of a hillock structure. The hillock is formed when a high angle grain boundary, after absorbing additional dislocations, becomes unstable and folds to form a new crystal grain that protrudes from the free surface.     

Recoverable creep deformation and transient local stress concentration due to heterogeneous grain-boundary diffusion and sliding in polycrystalline solids

Numerical simulations are used to investigate the influence of heterogeneity in grain-boundary diffusivity and sliding resistance on the creep response of a polycrystal. We model a polycrystal as a two-dimensional assembly of elastic grains, separated by sharp grain boundaries. The crystal deforms plastically by stress driven mass transport along the grain boundaries, together with grain-boundary sliding. Heterogeneity is idealized by assigning each grain boundary one of two possible values of diffusivity and sliding viscosity. We compute steady state and transient creep rates as functions of the diffusivity mismatch and relative fractions of grain boundaries with fast and slow diffusion. In addition, our results show that under transient conditions, flux divergences develop at the intersection between grain boundaries with fast and slow diffusivity, which generate high local stress concentrations. The stress concentrations develop at a rate determined by the fast diffusion coefficient, and subsequently relax at a rate determined by the slow diffusion coefficient. The influence of the mismatch in diffusion coefficient, loading conditions, and material properties on the magnitude of this stress concentration is investigated in detail using a simple model problem with a planar grain boundary. The strain energy associated with these stress concentrations also makes a small fraction of the plastic strain due to diffusion and sliding recoverable on unloading. We discuss the implications of these results for conventional polycrystalline solids at high temperatures and for nanostructured materials where grain-boundary diffusion becomes one of the primary inelastic deformation mechanisms even at room temperature.     

Discrete dislocation plasticity analysis of the grain size dependence of the flow strength of polycrystals

The grain size dependence of the flow strength of polycrystals is analyzed using plane strain, discrete dislocation plasticity. Dislocations are modeled as line singularities in a linear elastic solid and plasticity occurs through the collective motion of large numbers of dislocations. Constitutive rules are used to model lattice resistance to dislocation motion, as well as dislocation nucleation, dislocation annihilation and the interaction with obstacles. The materials analyzed consist of micron scale grains having either one or three slip systems and two types of grain arrangements: either a checker-board pattern or randomly dispersed with a specified volume fraction. Calculations are carried out for materials with either a high density of dislocation sources or a low density of dislocation sources. In all cases, the grain boundaries are taken to be impenetrable to dislocations. A Hall–Petch type relation is predicted with Hall–Petch exponents ranging from ≈0.3 to ≈1.6 depending on the number of slip systems, the grain arrangement, the dislocation source density and the range of grain sizes to which a Hall–Petch expression is fit. The grain size dependence of the flow strength is obtained even when no slip incompatibility exists between grains suggesting that slip blocking/transmission governs the Hall–Petch effect in the simulations.     

接触应力对FCB车轮钢组织演变与性能的影响

利用双盘滚动摩擦磨损试验机进行了贝氏体车轮钢的滚动磨损试验,采用扫描电子显微镜(SEM)、电子背散射衍射(EBSD)分析不同接触应力条件下贝氏体车轮钢次表层微观组织演变. 结果表明:在滚动磨损条件下,磨损机制由黏着磨损转变为疲劳磨损,增大接触应力对黏着磨损阶段的磨损量影响不大,但会显著增加疲劳磨损阶段的磨损量;贝氏体车轮钢在塑性变形的过程中,贝氏体铁素板条中位错逐渐增值、先累积形成小角度晶界,而后形成大角度晶界,使贝氏体铁素体发生细化;接触应力的大小影响表层组织的演变,当接触应力增至1 150 MPa时,晶粒细化为超细等轴晶,继续增加接触应力,组织变化并不明显. 接触应力大小会影响贝氏体车轮钢的表面硬度. 接触应力增加使贝氏体车轮钢的表面硬度增高,硬化层深度增大.      

Activation volume and deviation from Cottrell–Stokes law at small grain size

Dependence of activation volume with flow stress is examined for metals with grain size lower than 0.3 μm and larger than few tens of nanometers, where plastic deformation is most likely to be governed by a combination of grain boundary sliding and dislocations activity. The experimentally observed deviation from the classic linear behavior given by Cottrell–Stokes law [Basinski, Z.S., 1974. Forest hardening in face centered cubic metals. Scripta Metallurgica 8, 1301–1308] is analyzed, thanks to a modified Orowan equation taking into account of the grain boundaries sliding coupled to dislocations activity. These results are compared to experimental measurements of the activation volume, between room temperature and 120 °C, for a copper nanostructure with a grain size of 100 nm. A constant activation volume is observed at low stress (or high temperature) followed by an increase of activation volume with stress (inverse Cottrell–Stokes behavior). This analysis follows our initial experiments on this fine grained metal prepared by powder metallurgy, which exhibits ductility at near constant stress and strain rate [Champion, Y., Langlois, C., Guérin-Mailly, S., Langlois, P., Bonnentien, J.-L., Hÿtch, M.J., 2003. Near-perfect elastoplasticity in pure nanocrystalline copper. Science 300, 310–311].     

铜合金微塑性成形力学性能及其本构关系

对不同晶粒大小、不同特征尺寸的H62黄铜箔进行微拉伸实验,分析试样晶粒大小和特征尺寸对材料变形行为的影响。随着晶粒尺寸的减小,试样拉伸屈服应力逐渐增大,晶粒尺寸对屈服应力的影响满足Hall-Petch细晶强化关系;屈服强度随厚度的减小先减弱而后增强,随宽度的减小而增强;晶体塑性理论、表面层模型可以解释延伸率、抗拉强度随比表面积的增大而减小的现象。在实验数据的基础上通过修正双线性模型建立微塑性成形本构模型。     

Lithium Diffusion and Stress in a Polycrystalline Film Electrode

In the present work, the two-dimensional analytical solution for Li diffusion in grains and grain boundaries of a polycrystalline film electrode is established with consideration of Li-segregation at the grain boundary. The stress field induced by the inhomogeneity of Li concentration, called chemical stress here for brevity, is analyzed via the finite element simulation.The effects of the grain boundary including its size, its diffusion coefficient as well as the segregation phenomenon on the solute concentration and the chemical stress are examined. It shows that grain boundaries can assist fast diffusion and significantly affect the stress profile in the whole film. It proves that tailoring the grain boundary size or other grain boundary-related parameters may provide an alternative strategy for improving the overall diffusivity and mechanical stability of battery electrodes.     

Discrete dislocation dynamics simulations to interpret plasticity size and surface effects in freestanding FCC thin films

Strong size effects have been experimentally observed when microstructural features approach the geometric dimensions of the sample. In this work experimental investigations and discrete dislocation analyses of plastic deformation in metallic thin films have been performed. Columnar grains representative of the film microstructure are here considered. Simulations are based on the assumptions that sources are scarcely available in geometrically confined systems and nucleation sites are mainly located at grain boundaries. Especially, we investigated the influence on the mesoscopic constitutive response of the two characteristic length scales, i.e., film thickness and grain size. The simulated plastic response qualitatively reproduces the experimentally observed size effects while the main deformation mechanisms appear to be in agreement with TEM analyses of tested samples. A new interpretation of size scale plasticity is here proposed based on the probability of activating grain boundary dislocation sources. Moreover, the key role of a parameter such as the grain aspect ratio is highlighted. Finally, the unloading behavior has been investigated and a strong size dependent Bauschinger effect has been found. An interpretation of these phenomena is proposed based on the analysis of the back stress distribution within the samples.