一维六方准晶中平行螺型位错与有限长裂纹的相互作用

本文采用解析函数理论,利用Muskhelishvili著名的复变函数方法,给出了无限大一维六方准晶中多条平行螺型位错与有限长直裂纹相互作用的解析解,并得到了裂纹尖端的场强度因子和作用在位错上的像力.数值算例讨论了镜像力、场强度因子随位错位置的变化,及位错的位置与分布对广义应力场的影响．     

Determination of interaction forces between parallel dislocations by the evaluation of J integrals of plane elasticity

The Peach–Koehler expressions for the glide and climb components of the force exerted on a straight dislocation in an infinite isotropic medium by another straight dislocation are derived by evaluating the plane and antiplane strain versions of J integrals around the center of the dislocation. After expressing the elastic fields as the sums of elastic fields of each dislocation, the energy momentum tensor is decomposed into three parts. It is shown that only one part, involving mixed products from the two dislocation fields, makes a nonvanishing contribution to J integrals and the corresponding dislocation forces. Three examples are considered, with dislocations on parallel or intersecting slip planes. For two edge dislocations on orthogonal slip planes, there are two equilibrium configurations in which the glide and climb components of the dislocation force simultaneously vanish. The interactions between two different types of screw dislocations and a nearby circular void, as well as between parallel line forces in an infinite or semi-infinite medium, are then evaluated.     

The interaction between a screw dislocation and a rigid wedge inhomogeneity with an elastic circular inhomogeneity at the tip

The interaction between a screw dislocation and an elastic semi-cylindrical inhomogeneity abutting on a rigid half-plane is investigated. Utilizing the image dislocations method, the closed form solutions of the stress fields in the matrix and the inhomogeneity region are derived. The image force acting on the dislocation is also calculated. The results were used to study the interaction between a screw dislocation and a rigid wedge inhomogeneity with an elastic circular inhomogeneity at the tip by means of conformal mapping. The results show that an unstable equilibrium point of the dislocation near the semi-cylindrical inhomogeneity is found when the inhomogeneity is softer than the matrix. Moreover, the force on the dislocation is strongly affected by the position of the dislocation and the shear modulus of the semi-circular inhomogeneity. Positive screw dislocations can reduce the SIF of the rigid wedge inhomogeneity (shielding effect) only when it located in the lower half-plane. The shielding effect increases with the increase of the shear modulu of both the matrix and the inhomogeneity and increases with the increase of the wedge angle. The shielding effect (or anti-shielding effect) reaches the maximum when the dislocation tends to the wedge inhomogeneity interface.     

Continuum representation of systems of dislocation lines: A general method for deriving closed-form evolution equations

Plasticity is governed by the evolution of, in general anisotropic, systems of dislocations. We seek to faithfully represent this evolution in terms of density-like variables which average over the discrete dislocation microstructure. Starting from T. Hochrainer's continuum theory of dislocations (CDD) (Hochrainer, 2015), we introduce a methodology based on the ‘Maximum Information Entropy Principle’ (MIEP) for deriving closed-form evolution equations for dislocation density measures of different order. These equations provide an optimum representation of the kinematic properties of systems of curved and connected dislocation lines with the information contained in a given set of density measures. The performance of the derived equations is benchmarked against other models proposed in the literature, using discrete dislocation dynamics simulations as a reference. As a benchmark problem we study dislocations moving in a highly heterogeneous, persistent-slip-band like geometry. We demonstrate that excellent agreement with discrete simulations can be obtained in terms of a very small number of averaged dislocation fields containing information about the edge and screw components of the total and excess (geometrically necessary) dislocation densities. From these the full dislocation orientation distribution which emerges as dislocations move through a channel-wall structure can be faithfully reconstructed.     

Piezoelectric screw dislocations interacting with a circular inclusion with imperfect interface

The electroelastic coupling interaction between multiple screw dislocations and a circular inclusion with an imperfect interface in a piezoelectric solid is investigated. The appointed screw dislocation may be located either outside or inside the inclusion and is subjected to a line charge and a line force at the core. The analytic solutions of electroelastic fields are obtained by means of the complex-variable method. With the aid of the generalized Peach–Koehler formula, the explicit expressions of image forces exerted on the piezoelectric screw dislocations are derived. The motion and the equilibrium position of the appointed screw dislocation near the circular interface are discussed for variable parameters (interface imperfection, material electroelastic mismatch, and dislocation position), and the influence of the nearby parallel screw dislocations is also considered. It is found that the piezoelectric screw dislocation is always attracted by the electromechanical imperfect interface. When the interface imperfection is strong, the impact of material electroelastic mismatch on the image force and the equilibrium position of the dislocation becomes weak. Additionally, the effect of the nearby dislocations on the mobility of the appointed dislocation is very important.     

Analytic Solution for a Line Edge Dislocation in a Bimaterial System Incorporating Interface Elasticity

We examine the plane strain deformations of a bimaterial system consisting of a line edge dislocation interacting with a flat interface between two dissimilar isotropic half-planes in which the additional effects of interface elasticity are incorporated into the model of deformation. The entire system is assumed to be free of any external loading. Despite the fact that it is generally accepted that the separate interface modulus describing interface elasticity is permitted to take negative values, we show that simple closed-form solutions for the dislocation-induced stress field and the image force acting on the dislocation are available only when the interface modulus is assumed to be positive; the corresponding system admits no valid solutions when the interface modulus is negative. We present several numerical examples to illustrate our solutions. Additionally, we show that the influence of interface elasticity on the dislocation-induced interfacial stress field decays with increasing hardness of the adjoining half-plane (free of the dislocation). Moreover, we find that for a given (positive) in-plane interface modulus, the corresponding interface effects on the image force (acting on the dislocation) can reach maximum or minimum values when the Burgers vector of the dislocation is either parallel or perpendicular to the interface.     

A screw dislocation in a functionally graded material using the translation gauge theory of dislocations

The aim of this paper is to provide new results and insights for a screw dislocation in functionally graded media within the gauge theory of dislocations. We present the equations of motion for dislocations in inhomogeneous media. We specify the equations of motion for a screw dislocation in a functionally graded material. The material properties are assumed to vary exponentially along the x and y-directions. In the present work we give the analytical gauge field theoretic solution to the problem of a screw dislocation in inhomogeneous media. Using the dislocation gauge approach, rigorous analytical expressions for the elastic distortions, the force stresses, the dislocation density and the pseudomoment stresses are obtained depending on the moduli of gradation and an effective intrinsic length scale characteristic for the functionally graded material under consideration.     

Some half-space problems of cubic piezoelectric materials

Half-space problems of a cubic piezoelectric material subjected to anti-plane deformation and in-plane electric field are studied. A general solution in terms of the integration of the boundary data prescribed over the surface of the semi-infinite domain is derived. Based on the general solution, the problem of a concentrated line force acting on the surface is treated and ensuing electromechanical response is determined. The solution to the problem of a screw dislocation in the half-space is also obtained, and the result is exploited to study a sub-surface crack problem by simulating the crack as a continuous distribution of dislocations.     

Chemically and mechanically driven creep due to generation and annihilation of vacancies with non-ideal sources and sinks

The classical concept of Nabarro creep is extended for a general dislocation microstructure. The specific mechanism of the creep consists in generation and annihilation of vacancies at dislocation jogs acting as non-ideal sources and sinks for vacancies. This mechanism causes the climb of dislocations, allowing for local volume and shape change. The final kinetic equations, relating the dislocation microstructure and the local stress state to the creep rate, are derived by means of the thermodynamic extremal principle. Closed-form equations for the creep rate are derived for isotropic polycrystals. Based on the model the creep rate in the ferritic P-91 type steel at very low applied stress is evaluated and compared with experiment.     

Image decomposition method for the analysis of a mixed dislocation in a general multilayer

The elastic solutions for a mixed dislocation in a general multilayer with N dissimilar anisotropic layers are obtained via a generalized image decomposition method. The original problem is decomposed into N homogeneous subproblems with strategically placed continuously distributed image (virtual) dislocations which satisfy the consistency conditions for degenerate N − M (M < N) layer problems. The image dislocations are used to satisfy the interface or free surface conditions, and represent the unknowns of the problem. The resulting singular Cauchy integral equations are transformed into non-singular Fredholm integral equations of the second kind using certain H- and I-integral transforms. The Fredholm integral equations are then solved via the classical Nyström method. The general decomposition and the elimination of all singular integrals yield an exact formulation of the problem; the approximation arises only in the Nyström method. The dislocation mixity and the number of layers dissimilar in thickness and elastic anisotropy can be handled without difficulty, constrained only by the number of linear algebraic equations in the Nyström method for large N. For the numerical study, image forces on a dislocation in two- and three-layer systems are calculated. The accuracy of the results is verified by checking the boundary conditions and by comparison with previous results. The dependence of the image force on the dislocation position and mixity, and on the layer thicknesses and elastic anisotropies, is also illustrated via numerical investigations.     

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.     

Study of size effects in thin films by means of a crystal plasticity theory based on DiFT

In a recent publication, we derived the mesoscale continuum theory of plasticity for multiple-slip systems of parallel edge dislocations, motivated by the statistical-based nonlocal continuum crystal plasticity theory for single-glide given by Yefimov et al. [2004b. A comparison of a statistical-mechanics based plasticity model with discrete dislocation plasticity simulations. J. Mech. Phys. Solids 52, 279-300]. In this dislocation field theory (DiFT) the transport equations for both the total dislocation density and geometrically necessary dislocation (GND) density on each slip system were obtained from the Peach-Koehler interactions through both single and pair dislocation correlations. The effect of pair correlation interactions manifested itself in the form of a back stress in addition to the external shear and the self-consistent internal stress. We here present the study of size effects in single crystalline thin films with symmetric double slip using the novel continuum theory. Two boundary value problems are analyzed: (1) stress relaxation in thin films on substrates subject to thermal loading, and (2) simple shear in constrained films. In these problems, earlier discrete dislocation simulations had shown that size effects are born out of layers of dislocations developing near constrained interfaces. These boundary layers depend on slip orientations and applied loading but are insensitive to the film thickness. We investigate the stress response to changes in controlled parameters in both problems. Comparisons with previous discrete dislocation simulations are discussed.     

Viscoplastic flow rule for dislocation-mediated plasticity from systematic coarse-graining

The plastic response of metals is determined by the collective, coarse-grained dynamics of dislocations, rather than by the dynamics of individual dislocations. The evolution equations at both levels are quite different, for example considering their dependence on the applied mechanical load. On the one hand, the relation between the configurational force and dislocation velocity for individual dislocations is linear. On the other hand, in phenomenological crystal plasticity models, the relation between load and plastic slip is highly non-linear and often taken of power-law form. In this work, it is shown that this difference is justified and a consequence of emergent effects. Previously, an expression for the macroscopic dislocation flux was derived by systematic coarse graining (Kooiman et al., 2015). This expression has been evaluated numerically in this paper. The resulting relation between dislocation flux and applied mechanical load is found to be of power-law form with an exponent 3.7, while the underlying Discrete Dislocation Dynamics has a linear flux–load relation.     

含非完整界面涂层夹杂中螺型位错的弹性分析

研究了含非完整界面圆形涂层夹杂内部一个螺型位错在夹杂、涂层与无限大基体材料中产生的弹性场.运用复变函数函数方法,获得了三个区域复势函数的解析解答.利用求得的应力场和Peach-Koehler公式,得到了作用在螺型位错上位错力的精确表达式.主要讨论了两个非完整界面对位错力的影响规律.结果表明,涂层界面对夹杂内部螺型位错的吸引力随着界面粘结强度的弱化而变大.界面非完整程度增加削弱材料弹性失配对位错力的影响.在一定条件下,非完整界面可以改变夹杂内位错与涂层/基体系统之间的引斥干涉规律,并使位错在夹杂内部产生一个稳定或非稳定的平衡点.     

The mean free path of dislocations in nanoparticle and nanorod reinforced metal composites and implication for strengthening mechanisms

The mean free path length of dislocations plays an important role in the plastic behavior of metals, which may be significantly enhanced by the addition of nanoparticles. The effects of particle distribution pattern, particle size, shape and volume fraction on the dislocation mean free path length and average obstacle distance are analyzed in two- and three-dimensional models. As the particle volume fraction increases, the dislocation mean free path length and average distance between dislocations obstacles can be significantly reduced, implying enhancement of strength. The random particle distribution exhibits the best combination of dislocation mean free path length and average obstacle distance. In addition, both dislocation mean free path and average obstacle distance can be significantly reduced by changing the particle shape from nanosphere to nanorod. The present analysis may provide useful information for designing the particle enforced composite materials.     

Crack kinking in piezoelectric materials

A solution is presented for a class of two-dimensional electroelastic branched crack problems. Explicit Green's function for an interface crack subject to an edge dislocation is developed using the extended Stroh formulation allowing the branched crack problem to be expressed in terms of coupled singular integral equations. The integral equations are obtained by the method that models a kink as a continuous distribution of edge dislocations, and the dislocation density function is defined on the line of the branch crack only. Competition between crack extension along the interface and kinking into the substrate is investigated using the integral equations and the maximum energy release rate criterion. Numerical results are presented to show the effect of electric field on the path of crack extension. The work was supported by the Australian Research Council through a Queen Elizabeth II fellowship and by the Australian Academy of Science through the J.G. Russell Award.     

A gradient theory of single-crystal viscoplasticity that accounts for geometrically necessary dislocations

This study develops a gradient theory of single-crystal plasticity that accounts for geometrically necessary dislocations. The theory is based on classical crystalline kinematics; classical macroforces; microforces for each slip system consistent with a microforce balance; a mechanical version of the second law that includes, via the microforces, work performed during slip; a rate-independent constitutive theory that includes dependences on a tensorial measure of geometrically necessary dislocations. The microforce balances are shown to be equivalent to nonlocal yield conditions for the individual slip systems. The field equations consist of the yield conditions coupled to the standard macroscopic force balance; these are supplemented by classical macroscopic boundary conditions in conjunction with nonstandard boundary conditions associated with slip. As an aid to solution, a weak (virtual power) formulation of the nonlocal yield conditions is derived. To make contact with classical dislocation theory, the microstresses are shown to represent counterparts of the Peach-Koehler force on a single dislocation.     

Three-dimensional formulation of dislocation climb

We derive a Green's function formulation for the climb of curved dislocations and multiple dislocations in three-dimensions. In this new dislocation climb formulation, the dislocation climb velocity is determined from the Peach–Koehler force on dislocations through vacancy diffusion in a non-local manner. The long-range contribution to the dislocation climb velocity is associated with vacancy diffusion rather than from the climb component of the well-known, long-range elastic effects captured in the Peach–Koehler force. Both long-range effects are important in determining the climb velocity of dislocations. Analytical and numerical examples show that the widely used local climb formula, based on straight infinite dislocations, is not generally applicable, except for a small set of special cases. We also present a numerical discretization method of this Green's function formulation appropriate for implementation in discrete dislocation dynamics (DDD) simulations. In DDD implementations, the long-range Peach–Koehler force is calculated as is commonly done, then a linear system is solved for the climb velocity using these forces. This is also done within the same order of computational cost as existing discrete dislocation dynamics methods.     

THERMODYNAMIC VARIATIONAL APPROACH FOR CLIMB OF AN EDGE DISLOCATION

A general thermodynamic variational approach is applied to study the force on an edge dislocation, which drives the dislocation to climb. Our attention is focused on the physical mechanism responsible for dislocation climb. A dislocation in a material element climbs as a result of vacancies diffusing into or out from the dislocation core, with the dislocation acting as a source or a sink for vacancy diffusion in the material element. The basic governing equations for dislocation climb and the climb forces on the dislocation are obtained naturally as a result of the present thermodynamic variational approach.