Effects of boundary inclination and boundary type on shear-driven grain boundary migration

A series of molecular dynamics simulations was performed on a bicrystal to which a fixed shear rate was applied parallel to the boundary plane. Under some conditions, grain boundary motion is coupled to the relative tangential motion of the two grains. In order to investigate the generality of this type of coupled shear/boundary motion, simulations were performed for both special (low Σ) and general (non-Σ) [010] tilt boundaries over a wide range of grain boundary inclinations. The data point to the existence of two critical stresses: one for coupled shear/boundary motion and the other for grain boundary sliding. For the non-Σ boundaries, the critical stress for coupled shear/boundary motion is typically smaller than that for sliding; coupled shear/boundary motion occurs for all inclinations. For Σ5 boundaries, for which the critical stress is smaller and depends on boundary inclination, coupled shear/boundary motion occurs for some, but not all inclinations.     

包含倾斜晶界的双晶ZnO的热输运行为

用嵌入位错线法和重合位置点阵法构建含有小角度和大角度倾斜角的双晶氧化锌纳米结构.用非平衡分子动力学方法模拟双晶氧化锌在不同倾斜角度下的晶界能、卡皮查热阻,并研究了样本长度和温度对卡皮查热阻和热导率的影响.模拟结果表明,晶界能在小角度区域随倾斜角线性增加,而在大角度区域达到稳定,与卡皮查热阻的变化趋势一致.热导率随样本长度的增加而增加,卡皮查热阻表现出相反的趋势.然而随着温度的增加,热导率和卡皮查热阻都减小.通过比较含5.45°和38.94°晶界样本的声子态密度,发现声子光学支对热传导的影响不大,主要由声子声学支贡献,大角度晶界对声子散射作用更强,声学支波峰向低频率移动.     

Atomistic Simulation of Curvature Driven Grain Boundary Migration

We present two dimensional molecular dynamics simulations of grain boundary migration using the half-loop bicrystal geometry in the experiments of Shvindlerman et al. We examine the dependence of steady-state grain boundary migration rate on grain boundary curvature by varying the half-loop width at constant temperature. The results confirm the classical result derived by absolute reaction rate theory that grain boundary velocity is proportional to the curvature. We then measure the grain boundary migration rate for fixed half-loop width at varying temperatures. Analysis of this data establishes an Arrhenius relation between the grain boundary mobility and temperature, allowing us to extract the activation energy for grain boundary migration. Since grain boundaries have an excess volume, curvature driven grain boundary migration increases the density of the system during the simulations. In simulations performed at constant pressure, this leads to vacancy generation during the boundary migration, making the whole migration process jerky.     

Nanocrystalline gold with small size: inverse Hall–Petch between mixed regime and super-soft regime

ABSTRACT

Molecular dynamics simulations were used to study the atomic mechanisms of deformation of nanocrystalline gold with 2.65–18?nm in grain size to explore the inverse Hall–Petch effect. Based on the mechanical responses, particularly the flow stress and the elastic-to-plastic transition, one can delineate three regimes: mixed (10–18?nm, dislocation activities and grain boundary sliding), inverse Hall-Petch (5–10?nm, grain boundary sliding), and super-soft (below 5?nm). As the grain size decreases, more grain boundaries present in the nanocrystalline solids, which block dislocation activities and facilitate grain boundary sliding. The transition from dislocation activities to grain boundary sliding leads to strengthening-then-softening due to grain size reduction, shown by the flow stress. It was further found that, samples with large grain exhibit pronounced yield, with the stress overshoot decrease as the grain size decreases. Samples with grain sizes smaller than 5?nm exhibit elastic-perfect plastic deformation without any stress overshoot, leading to the super-soft regime. Our simulations show that, during deformation, smaller grains rotate more and grow in size, while larger grains rotate less and shrink in size.     

Continuous description of a grain boundary in forsterite from atomic scale simulations: the role of disclinations

We present continuous modelling at inter-atomic scale of a high-angle symmetric tilt boundary in forsterite. The model is grounded in periodic arrays of dislocation and disclination dipoles built on information gathered from discrete atomistic configurations generated by molecular dynamics simulations. The displacement, distortion (strain and rotation), curvature, dislocation and disclination density fields are determined in the boundary area using finite difference and interpolation techniques between atomic sites. The distortion fields of the O, Si and Mg sub-lattices are detailed to compare their roles in the accommodation of lattice incompatibility along the boundary. It is shown that the strain and curvature fields associated with the dislocation and disclination fields in the ‘skeleton’ O and Si sub-lattices accommodate the tilt incompatibility, whereas the elastic strain and rotation fields of the Mg sub-lattice are essentially compatible and induce stresses balancing the incompatibility stresses in the overall equilibrium.     

Quantifying the grain boundary resistance against slip transfer by experimental combination of geometric and stress approach using stage-I-fatigue cracks

In recent studies, many groups have investigated the interaction of dislocations and grain boundaries by bi-crystals and micro-specimen experiments. Partially, these experiments were combined with supplementary simulations by discrete dislocation dynamics, but quantitative data for the grain boundary resistance against slip transfer is still missing. In this feasibility study with first results, we use stage-I-fatigue cracks as highly localised sources for dislocations with well-known Burgers vectors to study the interaction between dislocations in the plastic zone in front of the crack tip and selected grain boundaries. The stress concentration at the grain boundary is calculated with the dislocation-free zone model of fracture using the dislocation density distribution in the plastic zone from slip trace height profile measurements by atomic force microscopy. The grain boundary resistance values calculated from common geometric models are compared to the local stress distribution at the grain boundaries. Hence, it is possible to quantify the grain boundary resistance and to combine geometric and stress approach for grain boundary resistance against slip transfer to a self-contained concept. As a result, the prediction of the grain boundary resistance effect based on a critical stress concept is possible with knowledge of the geometric parameters of the grain boundary only, namely the orientations of both participating grains and the orientation of the grain boundary plane.     

Atomistic simulation of the stacking fault energy and grain shape on strain hardening behaviours of FCC nanocrystalline metals

ABSTRACT

Ultra-fine grained copper with nanotwins is found to be both strong and ductile. It is expected that nanocrystalline metals with lamella grains will have strain hardening behaviour. The main unsolved issues on strain hardening behaviour of nanocrystalline metals include the effect of stacking fault energy, grain shape, temperature, strain rate, second phase particles, alloy elements, etc. Strain hardening makes strong nanocrystalline metals ductile. The stacking fault energy effects on the strain hardening behaviour are studied by molecular dynamics simulation to investigate the uniaxial tensile deformation of the layer-grained and equiaxed models for metallic materials at 300?K. The results show that the strain hardening is observed during the plastic deformation of the layer-grained models, while strain softening is found in the equiaxed models. The strain hardening index values of the layer-grained models decrease with the decrease of stacking fault energy, which is attributed to the distinct stacking fault width and dislocation density. Forest dislocations are observed in the layer-grained models due to the high dislocation density. The formation of sessile dislocations, such as Lomer–Cottrell dislocation locks and stair-rod dislocations, causes the strain hardening behaviour. The dislocation density in layer-grained models is higher than that in the equiaxed models. Grain morphology affects dislocation density by influencing the dislocation motion distance in grain interior.     

Multiscale modelling of dislocation/grain boundary interactions. II. Screw dislocations impinging on tilt boundaries in Al

The interaction of dislocations with grain boundaries (GBs) determines a number of important aspects of the mechanical performance of materials, including strengthening and fatigue resistance. Here, the coupled atomistic/discrete-dislocation (CADD) multiscale method, which couples a discrete dislocation continuum region to a fully atomistic region, is used to study screw-dislocations interacting with Σ3, Σ11, and Σ9 symmetric tilt boundaries in Al. The low-energy Σ3 and Σ11 boundaries absorb lattice dislocations and generate extrinsic grain boundary dislocations (GBDs). As multiple screw dislocations impinge on the GB, the GBDs form a pile-up along the GB and provide a back stress that requires increasing applied load to push the lattice dislocations into the GB. Dislocation transmission is never observed, even with large GBD pile-ups near the dislocation/GB intersection. Results are compared with experiments and previous, related simulations. The Σ9 grain boundary, composed from a more complex set of structural units, absorbs screw dislocations that remain localized, with no GBD formation. With increasing applied stress, new screw dislocations are then nucleated into the opposite grain from structural units in the GB that are nearby but not at the location where the original dislocation intersected the boundary. The detailed behaviour depends on the precise location of the incident dislocations and the extent of the pile-up. Transmission can occur on both Schmid and non-Schmid planes and can depend on the shear stresses on the GB plane. A continuum yield locus for transmission is formulated. In general, the overall dissociation and/or transmission behaviour is also determined by the Burgers vectors and associated steps of the primitive vectors of the grain boundary, and the criteria for dislocation transmission formulated by Lee et al . [Scripta Metall. 23 799 (1989); Phil. Mag. A 62 131 (1990); Metall. Trans. A 21 2437 (1990)] are extended to account for these factors.     

Dislocation cross-slip in nanocrystalline fcc metals

Constant strain rate molecular dynamics simulations of nanocrystalline Al demonstrate that a significant amount of dislocations that have nucleated at the grain boundaries, exhibit cross-slip via the Fleischer mechanism as they propagate through the grain. The grain boundary structure is found to strongly influence when and where cross-slip occurs, allowing the dislocation to avoid local stress concentrations that otherwise can act as strong pinning sites for dislocation propagation.     

Twin boundary and grain boundary effects on cyclic responses of tensile pre-deformed highly oriented nanotwinned Cu: Molecular dynamics simulation

Molecular dynamics method is performed for analyzing the relationship of the twin boundary and grain boundary on the cyclic response of nanotwinned Cu. Results show that the strength difference among the grain boundary, the twin boundary and the variation of dislocation density are nearly 2-2.5 times. We predict twin boundary is only a factor that affects the stable response, however, the dislocation form and the time to reach stability is caused by the grain boundary. Furthermore, the phenomenon of cyclic hardening is found in all the nanotwinned Cu samples.     

Molecular dynamics simulations of He bubble nucleation at grain boundaries

The nucleation behavior of He bubbles in single-crystal (sc) and nano-grain body-centered-cubic (bcc) Mo is simulated using molecular dynamics (MD) simulations, focusing on the effects of the grain boundary (GB) structure. In sc Mo, the nucleation behavior of He bubbles depends on irradiation conditions. He bubbles nucleate by either clustering of He atoms with pre-existing vacancies or self-interstitial-atom (SIA) punching without initial vacancies. In nano-grain Mo, strong precipitation of He at the GBs is observed, and the density, size and spatial distribution of He bubbles vary with the GB structure. The corresponding He bubble density is higher in nano-grain Mo than that in sc Mo and the average bubble size is smaller. In the GB plane, He bubbles distribute along the dislocation cores for GBs consisting of GB dislocations and randomly for those without distinguishable dislocation structures. The simulation results in nano-grain Mo are in agreement with previous experiments in metal nano-layers, and they are further explained by the effect of excess volume associated with the GBs.     

Effects of void shape and orientation on the elastoplastic properties of spheroidally voided single-crystal and nanotwinned copper

ABSTRACT

Void nucleation, growth and coalescence have been identified as the leading cause of ductile damage in metallic materials. To understand the underlying deformation and damage mechanisms, extensive theoretical, experimental and simulation efforts have been attempted on spherically voided metals. In this work, molecular dynamics simulations are performed to analyze the uniaxial straining deformation behaviours of both single-crystal and nanotwinned copper materials embedded with a preexisting spheroidal void. The coupling effects among twin boundary, spheroidal void aspect ratio and orientation on unidirectional elastoplastic behaviours are systematically examined. The dislocation-induced plastic deformation mechanism is also examined and compared with the one due to a perfectly spherical cavity. Simulation results show that elastic modulus increases with both spheroidal void aspect ratio and orientation. So do the yield stress, the first peak stress and the plasticity index. Another peak stress exists for most cases, except for a prolate void embedded in nanotwinned specimens. The slope between peak stresses decreases with both the spheroidal aspect ratio and orientation. The incorporation of a twin boundary results in lower elastic modulus, higher yield strength and smaller plasticity index. For an oblate void, the twin boundary gives rise to more severe strain softening behaviour. The dislocation extraction algorithm illustrates that the continuous nucleation, propagation and reaction of dislocations emanated from both the void front and twin boundary are responsible for the ductile damage of spheroidally voided crystals. The lower dislocation densities found in nanotwinned specimens indicate the desired suppression effects of twin boundary on dislocation activities.     

Collective motion of atoms in grain boundary migration of a bcc metal

Molecular dynamics simulations of grain boundary (GB) migration of a bcc metal, tungsten, have been carried out. The GB is of asymmetrical ? 110? tilt type. Detailed examinations of atomic processes in the migration, show that the GB migration consists mainly of GB dislocation glides. Furthermore, each motion of a GB dislocation involves a cooperative motion of about three atoms on each of the atomic planes perpendicular to the tilt axis, leading to their realignment from the receding grain to the advancing grain. This collective motion is not synchronized in all of the atomic planes, but appears to be in two or three adjacent planes, suggesting a kink mechanism for glides of the GB dislocations.     

Investigation on quenching at a high-angle Cu grain boundary on an atomic scale

We have performed molecular dynamics simulations of structural changes due to quenching the melting interface at a Cu $\Sigma $5(310)/[001] symmetrical tilt grain boundary. The simulation results suggest that the grain boundary structures due to quenching are different from those due to heating up to the same temperature. The calculated atom density profiles show that the grain boundary structures can be significantly changed as they are quenched to quite low temperatures.     

单向拉伸作用下Cu(100)扭转晶界塑性行为研究

应用分子动力学方法研究了在不同扭转角度下的Cu（100）失配晶界位错结构,以及不同位错结构对晶界强度的影响.模拟结果表明：小角度扭转晶界上将形成失配位错网,失配位错密度随着晶粒之间的失配扭转角度的增加而增加.变形过程中,位错网每个单元中均产生位错形核扩展.位错之间的塞积作用影响晶界的屈服强度：随着位错网格密度的增加,位错之间的塞积作用增强,界面的屈服强度得到提高.大角度扭转晶界将形成面缺陷,在变形中位错由晶界角点处形核扩展,此时由于面缺陷位错开动应力趋于一致,因此晶界的临界屈服强度趋于定值. 关键词： 扭转晶界 失配位错网 强化机理 分子动力学     

Strain induced grain boundary premelting in bulk copper bicrystals

In bulk bicrystals strain induced grain boundary premelting (SIGBPM) occurs when heavy screw dislocation pileup can be held up to a certain high temperature, approximately 0.6T _M, where T _M is the melting point of bulk material in Kelvin. SIGBPM occurs at grain boundaries to which new twist component is added due to the rotation of both component crystals toward opposite direction about the axis perpendicular to the grain boundary plane. At the original grain boundary, grain boundary sliding takes place due to this relative rotation. In f.c.c. metals with relatively low stacking fault energies such as copper, nickel, brass(30Zn) and silver, dislocations dissociate into partials. Therefore high density tangled dislocations introduced during plastic deformation hardly loose. If these dislocations can be held to high temperatures, SIGBPM is promoted. Formation of static or dynamic recrystallized grains suppresses SIGBPM itself and the propagation of grain boundary cracks formed by SIGBPM.     

Non-planar grain boundary structures in fcc metals and their role in nano-scale deformation mechanisms

This work presents the results of a comparative molecular dynamics study showing that relaxed random grain boundary structures can be significantly non-planar at the nano-scale in fcc metals characterized by low stacking fault values. We studied the relaxed structures of random [1?1?0] tilt boundaries in a polycrystal using interatomic potentials describing Cu and Pd. Grain boundaries presenting non-planar features were observed predominantly for the Cu potential but not for the Pd potential, and we relate these differences to the stacking fault values. We also show that these non-planar structures can have a strong influence on dislocation emission from the grain boundaries as well as on grain boundary strain accommodation processes, such as grain boundary sliding. We studied the loading response in polycrystals of 40 nm grain size to a level of 9% strain and found that the non-planar grain boundaries favour dislocation emission as a deformation mechanism and hinder grain boundary sliding. This has strong implications for the mechanical behaviour of nano-crystalline materials, which is determined by the competition between dislocation activity and grain boundary accommodation of the strain. Thus, the two interatomic potentials for Cu and Pd considered in this work resulted in the same overall stress–strain curve, but significantly different fractions of the strain accommodated by the intergranular versus intragranular deformation mechanisms. Strain localization patterns are also influenced by the non-planarity of the grain boundary structures.     

Quantitative phase-field approach for simulating grain growth in anisotropic systems with arbitrary inclination and misorientation dependence

A phase-field approach for quantitative simulations of grain growth in anisotropic systems is introduced, together with a new methodology to derive appropriate model parameters that reproduce given misorientation and inclination dependent grain boundary energy and mobility in the simulations. The proposed model formulation and parameter choice guarantee a constant diffuse interface width and consequently give high controllability of the accuracy in grain growth simulations.     

Modelling of grain refinement driven by negative grain boundary energy

Abstract

Grain refinement can be described by the classical kinetic equation using a negative value of the specific grain boundary Gibbs energy. A respective overview is offered reporting according observations and simulations, particularly linked to grain boundary segregation. Classical grain growth model is used in the treatment of evolution of the distribution function during refinement. The adapted model requires defining nucleation rate of new grains, which significantly influences the kinetics of the system of grains. Moreover, a jump in the distribution function is allowed at a certain value of the grain radius R_J, which separates old grains from newly nucleated ones. Evolution equation for both the critical radius R_c and separation radius R_J (jump position) as well as for the dimension-free distribution (shape) function are derived. Illustrative examples for the evolution of the system parameters under various nucleation rates of newly generated grains are presented.     

Prediction of flow stress and textures of AZ31 magnesium alloy at elevated temperature

The viscoplastic behaviour of magnesium alloys at high temperatures leads to highly temperature-dependent mechanical properties. While at high strain rates a notable strain hardening response is observed, at low strain rates the material shows a smooth plastic response with negligible amount of hardening. This complicated behaviour is due to different deformation mechanisms that are active at different strain rate regimes, resulting in different strain rate sensitivity parameters. In this study we show, by utilizing both numerical simulations and experiments, that this behaviour can be predicted by a model that combines two deformation mechanisms, grain boundary sliding mechanism and dislocation glide mechanism. We discuss the importance of each deformation mechanism at different strain rate regimes based on the findings of modelling and experimental results for AZ3 magnesium alloy. By developing a model that includes the above-mentioned two deformation mechanism, the prediction of flow properties is expanded to a wide range of strain rate regimes compared to previous study. The obtained numerical findings for the stress–strain behaviour as well as texture evolution show good agreement with the experimental results.