Atomistic simulations of interactions between the 1 / 2⟨111⟩ edge dislocation and symmetric tilt grain boundaries in tungsten

The strength of polycrystals is largely controlled by the interaction between lattice dislocations and grain boundaries. The atomistic details of these interactions are difficult to discern even by advanced high-resolution microscopy methods. In this paper we present results of atomistic simulations of interactions between an edge dislocation and three symmetric tilt grain boundaries in body-centred cubic tungsten. Our simulations reveal that the outcome of the dislocation–grain-boundary interaction depends sensitively on the grain boundary structure, the geometry of the slip systems in neighbouring grains, and the precise location of the interaction within the grain boundary. A detailed analysis of the evolution of the grain boundary structures and local stress fields during dislocation absorption and transmission is provided.     

Computer simulation of the vacancy defects interaction with shuffle dislocation in silicon

The interactions between the 60^° shuffle dislocation and two different types of vacancy defects in silicon are separately studied via the molecular dynamics simulation method. The Stillinger–Weber potential is used to describe the atomic interactions. The results show that the dislocation slip velocity will decrease due to the interaction with the vacancy cluster (V ₆). The simulation also reveals that the divacancy will be absorbed by the dislocation. Meanwhile, a climbing of the dislocation occurs during their interactions. However, the divacancy has little effect on the dislocation slip velocity. Based on the above results, the decrease in threading dislocation density in SiGe/Si heterostructures with the use of low-temperature Si buffer layer may be explained.     

Multiscale dislocation dynamics simulations of shock-induced plasticity in small volumes

Multiscale dislocation dynamics plasticity (MDDP) was used to investigate shock-induced deformation in monocrystalline copper. In order to enhance the numerical simulations, a periodic boundary condition was implemented in the continuum finite element (FE) scale so that the uniaxial compression of shocks could be attained. Additionally, lattice rotation was accounted for by modifying the dislocation dynamics (DD) code to update the dislocations’ slip systems. The dislocation microstructures were examined in detail and a mechanism of microband formation is proposed for single- and multiple-slip deformation. The simulation results show that lattice rotation enhances microband formation in single slip by locally reorienting the slip plane. It is also illustrated that both confined and periodic boundary conditions can be used to achieve uniaxial compression; however, a periodic boundary condition yields a disturbed wave profile due to edge effects. Moreover, the boundary conditions and the loading rise time show no significant effects on shock–dislocations interaction and the resulting microstructures. MDDP results of high strain rate calculations are also compared with the predictions of the Armstrong–Zerilli model of dislocation generation and movement. This work confirms that the effect of resident dislocations on the strain rate can be neglected when a homogeneous nucleation mechanism is included.     

Atomistic simulation analysis of the effects of void interaction on void growth and coalescence in a metallic system

Material deformation caused by the interaction between defects is a significant factor of material fracture failure. The present study employs molecular dynamics simulations of single-void and double-void crystalline Ni atomic systems to investigate inter-void interactions. Furthermore, simulations showing the evolution of dislocations for three different crystallographic orientations are conducted to study the void growth and coalescence. The simulations also consider the effect of the radius of the secondary void on dislocation evolution. The results show that double-void systems are more prone to yield than single-void systems. Further microstructural analysis indicates that the interaction between voids is realized by dislocation reactions. The simulation results of the dislocation evolution of the three orientations reveal that a relationship exists between the evolution of the dislocation density and the stress-strain curve. At the initial stage of dislocation, the dislocation grows slowly, and consists of Shockley partial dislocation. The dislocation growth rate then increases significantly in the sharply declining stage of the stress-strain curve, where most of dislocations are Shockley partial dislocation. Analysis of the dislocation length during the overall simulation indicates that the dislocation length of the [110] orientation is the longest, followed by that of the [111] orientation and the [100] orientation, which has the shortest dislocation length.     

Dislocation patterns and the similitude principle: 2.5D mesoscale simulations

During plastic flow of crystalline solids, dislocations self-organize in the form of patterns, with a wavelength that is inversely proportional to stress. After four decades of investigations, the origin of this property is still under discussion. We show that dislocation patterns verifying the principle of similitude can be obtained from dynamics simulations of double slip. These patterns are formed in the presence of long- and short-range interactions, but they are not significantly modified when only short-range interactions are present. This new insight into dislocation patterning phenomena has important implications regarding current models.     

Dislocation glide through non-randomly distributed point obstacles

Abstract

Classical meso-scale models for dislocation–obstacle interactions have, by and large, assumed a random distribution of obstacles on the glide plane. While a good approximation in many situations, this does not represent materials where obstacles are clustered on the glide plane. In this work, we have investigated the statistical problem of a dislocation sampling a set of clustered point obstacles in the glide plane using a modified areal-glide model. The results of these simulations show two clear regimes. For weak obstacles, the spatial distribution does not matter and the critically resolved shear stress is found to be independent of the degree of clustering. In contrast, above a critical obstacle strength determined by the degree of clustering, the critical resolved shear strength becomes constant. It is shown that this behaviour can be explained semi-analytically by considering the probability of interaction between the dislocation line and obstacles at a given level of stress. The consequences for alloys exhibiting solute clustering are discussed.     

Static friction of sinusoidal surfaces: a discrete dislocation plasticity analysis

Discrete dislocation plasticity simulations are carried out to investigate the static frictional response of sinusoidal asperities with (sub)-microscale wavelength. The surfaces are first flattened and then sheared by a perfectly adhesive platen. Both bodies are explicitly modelled, and the external loading is applied on the top surface of the platen. Plastic deformation by dislocation glide is the only dissipation mechanism active. The tangential force obtained at the contact when displacing the platen horizontally first increases with applied displacement, then reaches a constant value. This constant is here taken to be the friction force. In agreement with several experiments and continuum simulation studies, the friction coefficient is found to decrease with the applied normal load. However, at odds with continuum simulations, the friction force is also found to decrease with the normal load. The decrease is caused by an increased availability of dislocations to initiate and sustain plastic flow during shearing. Again in contrast to continuum studies, the friction coefficient is found to vary stochastically across the contact surface, and to reach locally values up to several times the average friction coefficient. Moreover, the friction force and the friction coefficient are found to be size-dependent.     

Dislocation decoration and raft formation in irradiated materials

Experimental observations of dislocation decoration with self-interstitial atom (SIA) clusters and of SIA cluster rafts are analysed to establish the mechanisms controlling these phenomena in bcc metals. The elastic interaction between SIA clusters, and between clusters and dislocations is included in kinetic Monte Carlo (KMC) simulations of damage evolution in irradiated bcc metals. The results indicate that SIA clusters, which normally migrate by 1D glide, rotate due to their elastic interactions, and that this rotation is necessary to explain experimentally-observed dislocation decoration and raft formation in neutron-irradiated pure iron. The critical dose for raft formation in iron is shown to depend on the intrinsic glide/rotation characteristics of SIA clusters. The model is compared with experimental observations for the evolution of defect cluster densities (sessile SIA clusters and nano-voids), dislocation decoration characteristics and the conditions for raft formation.     

Evidence from numerical modelling for 3D spreading of [001] screw dislocations in Mg2SiO4 forsterite

Computer simulations have previously been used to derive the atomic scale properties of the cores of screw dislocations in Mg₂SiO₄ forsterite by direct calculation using parameterized potentials and via the Peierls–Nabarro model using density functional theory. We show that, for the [001] screw dislocation, the parameterized potentials reproduce key features of generalized stacking fault energies when compared to the density functional theory results, but that the predicted structure of the dislocation core differs between direct simulation and the Peierls–Nabarro model. The [001] screw dislocation is shown to exhibit a low-energy non-planar core. It is suggested that for this dislocation to move its core may need to change structure and form a high-energy planar structure similar to that derived from the Peierls–Nabarro model. This could lead to dislocation motion via an unlocking–locking mechanism and explain the common experimental observation of long straight screw dislocation segments in deformed olivine.     

Size dependent deformation behaviour and dislocation mechanisms in 〈1 0 0〉 Cu nanowires

Abstract

Molecular dynamics simulations have been performed to understand the size-dependent tensile deformation behaviour of 〈1 0 0〉 Cu nanowires at 10 K. The influence of nanowire size has been examined by varying square cross-section width (d) from 0.723 to 43.38 nm using constant length of 21.69 nm. The results indicated that the yielding in all the nanowires occurs through nucleation of partial dislocations. Following yielding, the plastic deformation in small size nanowires occurs mainly by slip of partial dislocations at all strains, while in large size nanowires, slip of extended dislocations has been observed at high strains in addition to slip of partial dislocations. Further, the variations in dislocation density indicated that the nanowires with d > 3.615 nm exhibit dislocation exhaustion at small strains followed by dislocation starvation at high strains. On the other hand, small size nanowires with d < 3.615 nm displayed mainly dislocation starvation at all strains. The average length of dislocations has been found to be same and nearly constant in all the nanowires. Both the Young’s modulus and yield strength exhibited a rapid decrease at small size nanowires followed by gradual decrease to saturation at larger size. The observed linear increase in ductility with size has been correlated with the pre- and post-necking deformation. Finally, dislocation–dislocation interactions leading to the formation of various dislocation locks, the dislocation–stacking fault interactions resulting in the annihilation of stacking faults and the size dependence of dislocation–surface interactions have been discussed.     

Dislocation motion in high strain-rate deformation

We present a systematic investigation of dislocation motion, dislocation interactions, and the collective behaviour of dislocations in high strain-rate deformation. Based on results from three-dimensional dislocation dynamics simulations, we find that employing the accurate, full-dynamics, equation of motion (i.e. that includes inertial effects) significantly changes the predictions of microstructural evolution and the macroscopic response compared to the commonly used overdamped equation of motion (i.e. with no inertial effects), especially at high strain rates (10³–10⁶ s^?1). While we find that inertial effects cannot be neglected, the net velocities are not high enough that ‘relativistic’ effects are important. We also present results on the effects of high strain rates on single-crystal deformation, which show good agreement with experimental trends, including increased hardening with increasing strain rate.     

Dislocation density-based constitutive model for the mechanical behaviour of irradiated Cu

Performance degradation of structural steels in nuclear environments results from the formation of a high number density of nanometre-scale defects. The defects observed in copper-based alloys are composed of vacancy clusters in the form of stacking fault tetrahedra and/or prismatic dislocation loops that impede the motion of dislocations. The mechanical behaviour of irradiated copper alloys exhibits increased yield strength, decreased total strain to failure and decreased work hardening as compared to their unirradiated behaviour. Above certain critical defect concentrations (neutron doses), the mechanical behaviour exhibits distinct upper yield points. In this paper, we describe the formulation of an internal state variable model for the mechanical behaviour of such materials subject to these (irradiation) environments. This model has been developed within a multiscale materials-modelling framework, in which molecular dynamics simulations of dislocation–radiation defect interactions inform the final coarse-grained continuum model. The plasticity model includes mechanisms for dislocation density growth and multiplication and for irradiation defect density evolution with dislocation interaction. The general behaviour of the constitutive (homogeneous material point) model shows that as the defect density increases, the initial yield point increases and the initial strain hardening decreases. The final coarse-grained model is implemented into a finite element framework and used to simulate the behaviour of tensile specimens with varying levels of irradiation-induced material damage. The simulation results compare favourably with the experimentally observed mechanical behaviour of irradiated materials.     

纳米Ni在77 K温度下压缩行为的研究

利用低温力学测试系统研究了电化学沉积纳米Ni在77 K温度下的压缩行为. 室温下纳米Ni 的屈服强度为 2.0 GPa, 77 K温度下的屈服强度为3.0 GPa, 压缩变形量则由室温的10%左右下降到5%. 借助应变速率敏感指数、激活体积、扫描电子显微和高分辨透射电子显微分析, 对纳米Ni的塑性变形机制进行了表征. 研究表明, 在77 K温度下的塑性变形主要是由晶界-位错协调变形主导, 晶界本征位错弓出后无阻碍地在晶粒内无位错区运动, 直至在相对晶界发生类似切割林位错行为. 同时分析了弓出位错的残留位错部分在协调塑性变形时起到的增加应变相容性和减小应力集中的作用. 利用晶界-位错协调机制和残留位错运动与温度及缺陷的相关性揭示了纳米Ni室温和77 K温度压缩性能差异的内在原因. 关键词： 塑性变形 强度 位错     

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.     

Hydrodynamic boundary conditions: An emergent behavior of fluid–solid interactions

By using the Onsager principle of minimum energy dissipation, the hydrodynamic boundary conditions at the fluid–solid interface are shown to be the natural emergent behavior of microscopic interactions that lead to the interfacial tension and the tangential friction at the fluid–solid interface [T. Qian, C. Qiu, P. Sheng, J. Fluid Mech. 611 (2008) 333]. This is satisfying because the equations of motion, e.g., the Stokes equation, and the hydrodynamic boundary conditions can now be derived from a unified framework. The resulting continuum hydrodynamic formulation yields predictions for immiscible two-phase flows that are in quantitative agreement with molecular dynamic simulations. In particular, the classical problem of the moving contact line is resolved. We also show results on the moving contact line over chemically patterned surfaces which exhibit striking nanoscale characteristics as well as sub-quadratic dependence of the moving contact line dissipation on its average velocity.     

掺硅类金刚石薄膜摩擦过程的分子动力学模拟

掺硅类金刚石(Si-DLC) 薄膜表现出优异的摩擦学性能, 在潮湿空气和高温中显示出极低的摩擦系数和很好的耐磨性, 但是许多实验表明Si-DLC膜的摩擦性能受其硅含量的影响很大. 因此, 本文采用分子动力学模拟的方法分别研究干摩擦和油润滑两种情况下不同硅含量的Si-DLC膜的摩擦过程. 滑移结果表明干摩擦时DLC膜和掺硅DLC膜之间生成了一层转移膜, 而油润滑时则为边界膜. 因此干摩擦时的摩擦力明显大于油润滑时的摩擦力. 少量添加硅确实能降低DLC膜的摩擦力, 但是硅含量大于20%后对DLC膜的摩擦行为几乎无影响. 干摩擦时硅含量对转移膜内键的数量影响很大, 转移膜内CC键和CSi键都先增加后减少, 滑移结束时几乎不含CSi键.     

Atomic-scale study of dislocation–stacking fault tetrahedron interactions. Part I: mechanisms

Stacking fault tetrahedra (SFTs) are formed under irradiation in fcc metals and alloys. The high number density of SFTs observed suggests that they should contribute to radiation-induced hardening and, therefore, be taken into account when estimating mechanical property changes of irradiated materials. The key issue in this is to describe the interaction between a moving dislocation and an individual SFT, which is distinguished by a small physical size of the order of ～1–10?nm. We have performed atomistic simulations of edge and screw dislocations interacting with SFTs of different sizes at different temperatures and strain rates. Five possible interaction outcomes have been identified, involving either partial absorption, or shearing or restoration of SFTs. The mechanisms that give rise to these processes are described and their dependence on interaction parameters, such as SFT size, dislocation–SFT geometry, temperature and stress/strain rate are determined. Mechanisms that help to explain the formation of defect-free channels cleared by gliding dislocations, as observed experimentally, are also discussed. Hardening due to the various mechanisms and their dependence on loading conditions will be presented in a following paper (Part II).     

Innere Reibung durch Mitschleppen von Punktfehlern. II

The internal friction caused by dragging of mobile point defects is examined for the case of strong point defect dislocation interaction. For the amplitude independent internal friction and for the amplitude dependent internal friction with low point defect concentration this yields qualitatively the same results as the opposite case of weak interaction. For large point defect concentration there is an additional damping by hysteresis.