体心立方金属钽Ⅱ型裂纹尖端缺陷萌生的多尺度分析

本文采用多尺度准连续介质法(quasi-continuum method, QC)模拟体心立方(body-centered-cubic, bcc)金属钽(Ta)Ⅱ型裂纹尖端位错的形核与发射过程,获得位错发射位置与应力强度因子关系曲线,分析裂纹尖端缺陷萌生过程,研究全位错分解以及扩展位错形成机理. 位错活动在不同阶段表现出不一致的特征,新位错的发射对于位错运动具有促进作用. 研究表明,裂纹扩展初始阶段首先萌生点缺陷,点缺陷随着加载强度增加会萌生新的点缺陷,点缺陷最终运动到边界,导致Ⅱ型断裂破坏. 在全位错发射之前有不全位错的形核与发射表明全位错的分解分步进行,从势能曲线上来看,也就是两个极小值点的形成机理不同. 关键词： 多尺度 准连续介质法 Ⅱ型裂纹 扩展位错     

Recombination-induced stacking faults: evidence for a general mechanism in hexagonal SiC

We report on optically induced nucleation and expansion of stacking faults in hexagonal SiC structures. The activation energy for partial dislocation glide under optical excitation is found to reduce to 0.25 +/- 0.05 eV, which is about 2 eV lower than for pure thermal activation. From the measurements of thermal activation and below-gap excitation spectroscopy of dislocation glide, we conclude that the elementary process controlling expansion of stacking faults is kink pair nucleation aided by the phonon-kick mechanism. We propose that solitons on 30 degrees Si(g) partials with a silicon core act as deep 2.4 eV + Ev trap sites, readily providing electron-hole recombination energy to enhance the motion of dislocations. Our results suggest that this is a general mechanism of structural degradation in hexagonal SiC.     

Critical conditions for dislocation nucleation at surface steps

Dislocation nucleation at a surface step is analyzed based on a general variational boundary integral formulation of the Peierls–Nabarro dislocation model. By modelling the surface step as part of the surface of a three-dimensional crack, the free surface effect is taken into account by transferring the half space problem into an equivalent one in the infinite medium. The profiles of embryonic dislocations, corresponding to the relative displacements between the two adjacent atomic layers along the slip planes, are then rigorously solved through the variational boundary integral method. The critical conditions for dislocation nucleation are determined by solving the stress-dependent activation energies required to activate the embryonic dislocations from their stable to unstable saddle-point configurations. In particular, the effect of the step geometry, such as the height of the step, the dip angle of the slip plane and the inclined angle of the step surface on dislocation nucleation, is quantitatively ascertained. The results show that the atomic-scale surface step can rapidly reduce the critical stress required for dislocation nucleation from the surface by nearly an order of magnitude. The decrease in critical stress as a function of the height of the step is more significant for slip planes with smaller dip angles and surface steps with smaller inclined angles.     

Nucleation of cracks near the free surface in deformed metallic nanomaterials with a bimodal structure

A theoretical model that effectively describes the nucleation of cracks in stress fields of dislocation pile-ups near the free surface in metallic nanomaterials with a bimodal structure has been developed. The dependences of the critical shear stress τ_c (for the formation of a crack with an equilibrium length of 10 nm on a dislocation pile-up near the surface) on the size d of a grain containing the dislocation pile-up have been calculated for copper with a bimodal structure. Theoretically, it has been found that the critical shear stress τ_c for the nucleation of a crack near the free surface in a nanomaterial with a bimodal structure is approximately 30% higher than that for the crack nucleation within the nanomaterial at a distance from the free surface.     

Delayed fracture in porous media

The fracture of porous media subjected to a constant load is studied. Contrary to homogeneous solids in which fracture usually happens instantaneously at a well-defined breaking strength, the fracture of a porous medium can occur with a delay, allowing us to quantify the average lifetime of the unbroken material. We show that the average fracture probability, a key property for risk analysis in civil engineering, is given by the probability of crack nucleation. The nucleation process can be understood qualitatively by calculating the activation energy for crack nucleation, taking into account the porosity of the medium.     

Dislocation dynamic modelling of the brittle–ductile transition in tungsten

Brittle–ductile transitions in metals, ceramics and semiconductors are closely connected with dislocation activity emanating near to crack-tips. We have simulated the evolution of crack-tip plasticity using a two-dimensional dislocation dynamics model which has been developed to include two symmetric slip planes intersecting the crack-tip, and applied to single-crystal tungsten. The dislocation mobility law used was physically based on double-kink nucleation on screw dislocations, with an activation energy reduced by the local stress. Even in the strong stress gradients near a crack-tip, the dislocations are found to self-organise so that the internal stress in the array is effectively constant with time and position over a wide range of strain rates and temperatures. The resultant net activation energy for dislocation motion is found to be constant and close to the activation energy experimentally measured for the brittle–ductile transition. Use of a fracture criterion based on the local crack-tip stress intensity factor, as modified by the stresses from the emitted dislocations, allows explicit prediction of the form and temperature of the brittle–ductile transition. Predictions are found to be in very close agreement with experiment.     

Harnessing atomistic simulations to quantify activation parameters for dislocation nucleation from a grain boundary in Nickel

We perform atomistic simulations and nudged elastic band calculations to quantify activation energy barriers for dislocation nucleation from ∑3 grain boundary containing a pre-existing void in Ni. By changing the void size, we offer a phenomenological relationship between the activation free energy at zero stress and boundary porosity. Simulations at different temperatures are also conducted to gain some insights into the inherent intricacies of activation energy landscape. It is envisioned that the approach can be pushed forward to guide the flow rules of physics based crystal plasticity models at higher length scales that account for grain boundary effects.     

Reaction pathway analysis for shuffle-set 60° perfect dislocation in Si

In this work, the EDIP potential is employed for representing silicon and the shuffle-set 60° perfect dislocation motion is investigated by reaction pathway analysis. There are three possible shuffle-set 60° perfect dislocation core structures named as S1, S2 and S3. The activation energy barriers of the kink migration and nucleation in S1and S2 types are calculated by CI-NEB method. The simulation results show that the critical resolved shear strain of the shuffle-set dislocation in S1 type is around 5%, and the S1 type is the dominate one in the shear strain region of 0 to 5%. During the shear strain from 5to 11.81%, the dislocation moves as the S1 core kink nucleation and migration, meanwhile the S1 dislocation core is in process of transforming into S2. More interestingly, both S1 and S2 dislocation core structures is observed along the dislocation line in this shear strain regime, which could response to the missing observation of long segment dislocation line in the experiment.     

Simulation of screw dislocation motion in iron by molecular dynamics simulations

Molecular dynamics (MD) simulations are used to investigate the response of a/2<111> screw dislocation in iron submitted to pure shear strain. The dislocation glides and remains in a (110) plane; the motion occurs exclusively through the nucleation and propagation of double kinks. The critical stress is calculated as a function of the temperature. A new method is developed and used to determine the activation energy of the double kink mechanism from MD simulations. It is shown that the differences between experimental and simulation conditions lead to a significant difference in activation energy. These differences are explained, and the method developed provides the link between MD and mesoscopic simulations.     

Atomic scale investigation of grain boundary structure role on intergranular deformation in aluminium

The role that grain boundary (GB) structure plays on the directional asymmetry of an intergranular crack (i.e. cleavage behaviour is favoured along one direction, while ductile behaviour along the other direction of the interface) was investigated using atomistic simulations for aluminium 〈1 1 0〉 symmetric tilt GBs. Middle-tension (M(T)) and Mode-I crack propagation specimens were used to evaluate the predictive capability of the Rice criterion. The stress–strain response of the GBs for the M(T) specimens highlighted the importance of the GB structure. The observed crack tip behaviour for certain GBs (Σ9 (2 2 1), Σ11 (3 3 2) and Σ33 (4 4 1)) with the M(T) specimen displayed an absence of directional asymmetry which is in disagreement with the Rice criterion. Moreover, in these GBs with the M(T) specimen, the dislocation emission from a GB source at a finite distance ahead of the crack tip was observed rather than from the crack tip, as suggested by the Rice criterion. In an attempt to understand discrepancy between the theoretical predictions and atomistic observations, the effect of boundary conditions (M(T), Mode-I and the edge crack) on the crack tip events was examined and it was concluded that the incipient plastic events observed were strongly influenced by the boundary conditions (i.e. activation of dislocation sources along the GB, in contrast to dislocation nucleation directly from the crack tip). In summary, these findings provide new insights into crack growth behaviour along GB interfaces and provide a physical basis for examining the role of the GB character on incipient event ahead of a crack tip and interface properties, as an input to higher scale models.     

The surface and interface nucleation of misfit dislocations as a possible source of asymmetric strain relaxation in vicinal heterostructures

The surface nucleation of misfit dislocations in vicinal (001) oriented heterostructures with compressive and tensile stresses is discussed. It is shown that beside the asymmetrical stressing of opposite dislocation slip planes due to the vicinal substrate, the surface steps have a similar effect. The effect of the steps has the same-sign asymmetry for a compressive stressed epilayer, but opposite for the tensile case. The effect on dislocation nucleation energy is calculated. For miscut angles used normally, the step energy contribution exceeds that of due to the vicinal substrate. The extension to interface nucleation is treated qualitatively.     

Corroboration of a multiscale approach with all atom calculations in analysis of dislocation nucleation from surface steps

We present a comparative study of the influence of atomic-scale surface steps on dislocation nucleation at crystal surfaces based on an all atom method and a hierarchal multiscale approach. The multiscale approach is based on the variational boundary integral formulation of the Peiersl–Nabarro dislocation model in which interatomic layer potentials derived from atomic calculations of generalized stacking fault energy surfaces are incorporated. We have studied nucleation of screw dislocations in two bcc material systems, molybdenum and tantalum, subjected to simple shear stress. Compared to dislocation nucleation from perfectly flat surfaces, the presence of atomic scale surface steps rapidly reduces the critical stress for dislocation nucleation by almost an order of magnitude as the step height increases. In addition, they may influence the slip planes on which dislocation nucleation occurs. The results of the all atom method and the multiscale approach are in good agreement, even for steps with height of only a single atomic layer. Such corroboration supports the further use of the multiscale approach to study dislocation nucleation phenomena in more realistic geometries of technological importance, which are beyond the reach of all current atom simulations.     

Dislocation Plasticity Effects on Interfacial Fracture

Analyses are reviewed where plastic flow in the vicinity of an interfacial crack is represented in terms of the nucleation and glide of discrete dislocations. Attention is confined to cracks along a metal-ceramic interface, with the ceramic idealized as being rigid. Both monotonic and fatigue loading are considered. The main focus is on the stress and deformation fields near the crack tip predicted by discrete dislocation plasticity, in comparison with those obtained from conventional continuum plasticity theory. The role that discrete dislocation plasticity can play in interpreting interface fracture properties in the presence of plastic flow is discussed.     

Dislocation nucleation at amorphous intergrain boundaries in deformed nanoceramics

A theoretical model is proposed for lattice dislocation nucleation in deformed nanocrystalline ceramics with amorphous intergrain boundaries. According to the model, a lattice dislocation dipole nucleates at an amorphous intergrain boundary through a local plastic shear along the boundary cross section. The energy parameters of this nucleation process are calculated. It is demonstrated that the dislocation nucleation at amorphous intergrain boundaries is energetically favorable and can occur as an athermic process (without energy barrier) in the nanocrystalline phase of cubic silicon carbide 3C-SiC and in the TiN/a-Si₃N₄ nanocomposite over wide ranges of structural parameters and mechanical loads.     

Temperature and strain-rate dependence of surface dislocation nucleation

Dislocation nucleation is essential to the plastic deformation of small-volume crystalline solids. The free surface may act as an effective source of dislocations to initiate and sustain plastic flow, in conjunction with bulk sources. Here, we develop an atomistic modeling framework to address the probabilistic nature of surface dislocation nucleation. We show the activation volume associated with surface dislocation nucleation is characteristically in the range of 1-10b3, where b is the Burgers vector. Such small activation volume leads to sensitive temperature and strain-rate dependence of the nucleation stress, providing an upper bound to the size-strength relation in nanopillar compression experiments.     

Second-phase nucleation on an edge dislocation

A model for nucleation of second phase at or around a dislocation in a crystalline solid is considered. The model employs the Ginzburg–Landau theory of phase transitions comprising the sextic term in the order parameter (η⁶) in the Landau free energy. The ground state solution of the linearised time-independent Ginzburg–Landau equation is derived, through which the spatial variation of the order parameter is delineated. Moreover, a generic phase diagram indicating tricritical behaviour near and away from the dislocation is depicted. The relation between classical nucleation theory and the Ginzburg–Landau approach is discussed, for which the critical formation energy of the nucleus is related to the maximum of the Landau potential energy. A numerical example illustrating the application of the model to the case of nucleation of hydrides in zirconium alloys is provided.     

α-Fe裂纹的分子动力学研究

通过分子动力学方法，模拟了α-Fe裂纹的单轴拉伸实验中的形变过程.研究了不同晶体取向裂纹的形变特点和断裂机理，观察到各种形变现象，如位错形核和发射，位错运动，堆垛层错或孪晶的形成，纳米空洞的形成与连接等.计算结果表明，裂纹扩展是塑性过程和弹性过程相结合的过程，其中塑性过程表现为由裂尖发射的位错导致的原子切变行为，而弹性过程的发生则是由无位错区中的原子断键所导致.同时还研究了α-Fe裂纹的形变特点和断裂机理与温度场和应力场的依赖关系.     

Homogeneous nucleation of glide dislocation loops in nanoceramics

A theoretical model is proposed for the homogeneous nucleation of glide dislocation loops in nanocrystalline ceramics under deformation at low and high temperatures. The nucleation of a dislocation loop in a crystalline grain is considered an ideal nanoscopic shear whose magnitude (the Burgers vector of the dislocation) increases gradually as the loop is nucleating. The characteristics of the homogeneous nucleation of glide dislocation loops in nanocrystalline ceramics based on cubic silicon carbide are calculated. It is shown that, in general, the homogeneous nucleation of a dislocation loop in nanocrystalline ceramics at high temperatures proceeds in two stages, namely, the athermal nucleation of a loop of a “noncrystallographic” partial dislocation and its thermally activated transformation into an ordinary partial lattice dislocation loop.     

Effect of triple junctions of nanotubes on strengthening and fracture toughness of ceramic nanocomposites

A theoretical model has been proposed for describing the influence of triple junctions of nanotubes on the strengthening of a nanocomposite. It has been assumed that the slip of nanotubes along the boundary with the matrix takes place via the nucleation and glide of prismatic dislocation loops enveloping the nanotubes. These loops are retarded by the triple junctions of nanotubes, which leads to a strengthening and increase in the fracture toughness (crack resistance) of the nanocomposite. It has been shown that, in order for the dislocation loop to overcome the triple junction, the shear stress acting on the loop should exceed a certain critical level. This critical stress increases as the radius and wall thickness of the nanotube decrease. The inference has been made that the triple junctions of thinnest nanotubes, such as single-walled carbon nanotubes, should lead to the greatest strengthening and to an increase in the crack resistance of these nanocomposites.