Line-integral representations for extended displacements,stresses, and interaction energy of arbitrary dislocation loops in transversely isotropic magneto-electro-elastic bimaterials

In addition to the hexagonal crystals of class 6 mm, many piezoelectric materials （e.g., BaTiO3）, piezomagnetic materials （e.g., CoFe2O4）, and multiferroic com-posite materials （e.g., BaTiO3-CoFe2O4 composites） also exhibit symmetry of transverse isotropy after poling, with the isotropic plane perpendicular to the poling direction. In this paper, simple and elegant line-integral expressions are derived for extended displace-ments, extended stresses, self-energy, and interaction energy of arbitrarily shaped, three-dimensional （3D） dislocation loops with a constant extended Burgers vector in trans-versely isotropic magneto-electro-elastic （MEE） bimaterials （i.e., joined half-spaces）. The derived solutions can also be simply reduced to those expressions for piezoelectric, piezo-magnetic, or purely elastic materials. Several numerical examples are given to show both the multi-field coupling effect and the interface/surface effect in transversely isotropic MEE materials.     

A non-singular continuum theory of dislocations

We develop a non-singular, self-consistent framework for computing the stress field and the total elastic energy of a general dislocation microstructure. The expressions are self-consistent in that the driving force defined as the negative derivative of the total energy with respect to the dislocation position, is equal to the force produced by stress, through the Peach-Koehler formula. The singularity intrinsic to the classical continuum theory is removed here by spreading the Burgers vector isotropically about every point on the dislocation line using a spreading function characterized by a single parameter a, the spreading radius. A particular form of the spreading function chosen here leads to simple analytic formulations for stress produced by straight dislocation segments, segment self and interaction energies, and forces on the segments. For any value a>0, the total energy and the stress remain finite everywhere, including on the dislocation lines themselves. Furthermore, the well-known singular expressions are recovered for a=0. The value of the spreading radius a can be selected for numerical convenience, to reduce the stiffness of the dislocation equations of motion. Alternatively, a can be chosen to match the atomistic and continuum energies of dislocation configurations.     

On the micromechanics of void growth by prismatic-dislocation loop emission

Experimental evidence and recent molecular dynamics simulations of void growth indicate that prismatic dislocation loop emission by externally applied stresses is a viable mechanism of void growth under shock loading conditions when diffusive processes are given no time to operate. In this paper, the process of growth by loop emission is studied in a model system comprised of a void in an infinite linearly elastic and isotropic solid loaded axisymmetrically by remote applied stresses. First, the interaction between applied stresses, the stress field of a single dislocation loop or a pile-up of loops next to the void, the surface energy expenditure on void surface change, and the lattice resistance to the motion of loops is reviewed. The necessary condition for interstitial loop emission is used to determine the equilibrium positions of the loops as well as the maximum number of loops in a pile-up under given applied stresses. For the parameters of the model-material with purely hydrostatic loading, the numerical results yield a volume change for the void, which when normalized by the initial undeformed volume, exhibits a strong dependence on the size of the void for radii less than ∼400 times the lattice Burgers vector. For larger voids, the normalized volume change was found to be independent of the void radius.     

基体裂纹与界面刚性线的弹性干涉

研究两种材料界面上的刚性线与其它任意位置处直线裂纹弹性干涉的反平面问题。基于界面上刚性线与任意位置处螺型位错干涉的基本解,运用连续位错密度模型法将问题转化为奇异积分方程。用半开型积分法求解奇异积分方程,得到位错密度函数的离散值,计算裂纹尖端处的应力强度因子。算例说明该方法可用于工程实际问题。     

Image singularities of Green’s functions for isotropic elastic bimaterials subjected to concentrated forces and dislocations

Green’s functions for isotropic materials in the two-dimensional problem for elastic bimaterials with perfectly bonded interface are reexamined in the present study. Although the Green’s function for an isotropic elastic bimaterial subjected to a line force or a line dislocation has been discussed by many authors, the physical meaning and the structure of the solution are not clear. In this investigation, the Green’s function for an elastic bimaterial is shown to consist of eight Green’s functions for a homogeneous infinite plane. One of the novel features is that Green’s functions for bimaterials can be expressed directly by knowing Green’s functions for the infinite plane. If the applied load is located in material 1, the solution for the half-plane of material 1 is constructed with the help of five Green’s functions corresponding to the infinite plane. However, the solution for the half-plane of material 2 only consists of three Green’s functions for the infinite plane. One of the five Green’s functions of material 1 and all the three Green’s functions of material 2 have their singularities located in the half-plane where the load is applied, and the other four image singularities of material 1 are located outside the half-plane at the same distance from the interface as that of the applied load. The nature and magnitude of the image singularities for both materials are presented explicitly from the principle of superposition, and classified according to different loads. It is known that for the problem of anisotropic bimaterials subjected to concentrated forces and dislocations, the image singularities are simply concentrated forces and dislocations with the stress singularity of order O(1/r). However, higher orders (O(1/r²) and O(1/r³)) of stress singularities are found to exist in this study for isotropic bimaterials. The highest order of the stress singularity is O(1/r³) for the image singularities of material 1, and is O(1/r²) for material 2. Using the present solution, Green’s functions associated with the problems of elastic half-plane with free and rigidly fixed boundaries, for homogeneous isotropic elastic solid, are obtained as special cases.     

Displacement fields and self-energies of circular and polygonal dislocation loops in homogeneous and layered anisotropic solids

There are large classes of materials problems that involve the solutions of stress, displacement, and strain energy of dislocation loops in elastically anisotropic solids, including increasingly detailed investigations of the generation and evolution of irradiation induced defect clusters ranging in sizes from the micro- to meso-scopic length scales. Based on a two-dimensional Fourier transform and Stroh formalism that are ideal for homogeneous and layered anisotropic solids, we have developed robust and computationally efficient methods to calculate the displacement fields for circular and polygonal dislocation loops. Using the homogeneous nature of the Green tensor of order −1, we have shown that the displacement and stress fields of dislocation loops can be obtained by numerical quadrature of a line integral. In addition, it is shown that the sextuple integrals associated with the strain energy of loops can be represented by the product of a pre-factor containing elastic anisotropy effects and a universal term that is singular and equal to that for elastic isotropic case. Furthermore, we have found that the self-energy pre-factor of prismatic loops is identical to the effective modulus of normal contact, and the pre-factor of shear loops differs from the effective indentation modulus in shear by only a few percent. These results provide a convenient method for examining dislocation reaction energetic and efficient procedures for numerical computation of local displacements and stresses of dislocation loops, both of which play integral roles in quantitative defect analyses within combined experimental–theoretical investigations.     

Modeling core-spreading of interface dislocation and its elastic response in anisotropic bimaterial

Interfacial dislocation may have a spreading core corresponding to a weak shear resistance of interfaces. In this paper, a conic model is proposed to mimic the spreading core of interfacial dislocation in anisotropic bimaterials. By the Stroh formalism and Green's function, the analytical expressions of the elastic fields are deduced for such a dislocation. Taking Cu/Nb bimaterial as an example, it is demonstrated that the accuracy and efficiency of the method are well validated by the interface conditions, a spreading core can greatly reduce the stress intensity near the interfacial dislocation compared with the compact core, and the elastic fields near the spreading core region are significantly different from the condensed core, while they are less sensitive to a field point that is 1.5times the core width away from the center of the spreading core.     

Independence of the Presence of Cracks or Inclusions of the Net Interaction Force Acting on a Dislocation by a Skewed Line Singularity

Orlov and Indenbom [1] have shown that the net (integrated) interaction force F between two skew dislocations with Burgers vectors separated by a distance h in an infinite anisotropic elastic medium is independent of h. Nix [2] computed numerically the net interaction force F between two skew dislocations that are parallel to the traction-free surface X2=0 of an isotropic elastic half-space. His numerical results showed that F was independent of h; a partial result of what Barnett [3] called Nix"s theorem. The separation-independence portion of Nix"s theorem has been proved to hold for a general anisotropic elastic half-space with a traction-free, rigid, or slippery surface, and for bimaterials [3-5]. In this paper, we show that the net interaction force is independent of the presence of inclusions. We will consider the case in which the line dislocation b is a more general line singularity which can include a coincident line force with strength f per unit length of the line singularity. An inclusion is an infinitely long dissimilar anisotropic elastic cylinder of an arbitrary cross-section whose axis is parallel to the line singularity (f, b). The (skew) line dislocation does not intersect the inclusion. The special cases of an inclusion are a void, crack, or rigid inclusion. There can be more than one inclusion of different cross sections and different materials. The line singularity (f, b) can be outside the inclusions or inside one of the inclusions. The inclusions and the matrix need not have a perfect bonding. One can have a debonding with or without friction. This revised version was published online in August 2006 with corrections to the Cover Date.     

Three-dimensional elastic displacements induced by a dislocation of polygonal shape in anisotropic elastic crystals

Dislocations and the elastic fields they induce in anisotropic elastic crystals are basic for understanding and modeling the mechanical properties of crystalline solids. Unlike previous solutions that provide the strain and/or stress fields induced by dislocation loops, in this paper, we develop, for the first time, an approach to solve the more fundamental problem—the anisotropic elastic dislocation displacement field. By applying the point-force Green’s function for a three-dimensional anisotropic elastic material, the elastic displacement induced by a dislocation of polygonal shape is derived in terms of a simple line integral. It is shown that the singularities in the integrand of this integral are all removable. The proposed expression is applied to calculate the elastic displacements of dislocations of two different fundamental shapes, i.e. triangular and hexagonal. The results show that the displacement jump across the dislocation loop surface exactly equals the assigned Burgers vector, demonstrating that the proposed approach is accurate. The dislocation-induced displacement contours are also presented, which could be used as benchmarks for future numerical studies.     

Singularity-free dislocation dynamics with strain gradient elasticity

The singular nature of the elastic fields produced by dislocations presents conceptual challenges and computational difficulties in the implementation of discrete dislocation-based models of plasticity. In the context of classical elasticity, attempts to regularize the elastic fields of discrete dislocations encounter intrinsic difficulties. On the other hand, in gradient elasticity, the issue of singularity can be removed at the outset and smooth elastic fields of dislocations are available. In this work we consider theoretical and numerical aspects of the non-singular theory of discrete dislocation loops in gradient elasticity of Helmholtz type, with interest in its applications to three dimensional dislocation dynamics (DD) simulations. The gradient solution is developed and compared to its singular and non-singular counterparts in classical elasticity using the unified framework of eigenstrain theory. The fundamental equations of curved dislocation theory are given as non-singular line integrals suitable for numerical implementation using fast one-dimensional quadrature. These include expressions for the interaction energy between two dislocation loops and the line integral form of the generalized solid angle associated with dislocations having a spread core. The single characteristic length scale of Helmholtz elasticity is determined from independent molecular statics (MS) calculations. The gradient solution is implemented numerically within our variational formulation of DD, with several examples illustrating the viability of the non-singular solution. The displacement field around a dislocation loop is shown to be smooth, and the loop self-energy non-divergent, as expected from atomic configurations of crystalline materials. The loop nucleation energy barrier and its dependence on the applied shear stress are computed and shown to be in good agreement with atomistic calculations. DD simulations of Lomer–Cottrell junctions in Al show that the strength of the junction and its configuration are easily obtained, without ad-hoc regularization of the singular fields. Numerical convergence studies related to the implementation of the non-singular theory in DD are presented.     

Green's functions for anisotropic piezoelectric bimaterials with a slipping interface

An explicit full-field expression of the Green's functions for anisotropic piezoelectric bimaterials with a slipping interface is derived. When the electro-elastic singularity reduces to a pure dislocation in displacement and electric potential, interaction energy between the dislocation and the bimaterials is obtained explicitly while the generalized force on the dislocation is given in a real form which is also valid for degenerate materials. The investigation demonstrates that the boundary conditions at the slipping interface between two piezoelectric materials will exert a prominent influence on the mobility of the dislocation. Project supported by the National Natural Science Foundation of China (No. 59635140).     

The interaction between a screw dislocation and a piezoelectric fiber composite with a wedge crack

The interaction problem between a screw dislocation and a piezoelectric fiber composite with a semi-infinite wedge crack is investigated in this paper. The piezoelectric media are assumed to be transversely isotropic with the poling direction along the x ₃ direction. The screw dislocation considered here involves a Burgers vector parallel to the poling direction with a line force and a line charge being applied at the core of the dislocation. Both cases of the screw dislocation located at the matrix and inclusion are observed. The analytical derivation is based on the complex variable and the conformal mapping methods. The exact solutions are obtained to calculate the forces on the dislocation and the crack-tip stress and electric displacement intensity factors. Based on these results, the anti-shielding and shielding effects for different loadings, material combinations, and geometric configurations are discussed in detail.     

On the wedge dispiration in an inhomogeneous isotropic nonlinear elastic solid

In this short note we find the stress field of a wedge dispiration (a combination of a screw dislocation and a wedge dislocation along the same line) in an inhomogeneous incompressible isotropic nonlinear solid. We discuss the effect of the radial inhomogeneity of energy function on both the stress field and the energy per unit length of the dispiration and compare with those in an isotropic linear elastic solid.     

Interaction of a screw dislocation with a semi-infinite interfacial crack in a magneto-electro-elastic bi-material

The interaction between a screw dislocation and a semi-infinite interfacial crack in a transversely isotropic magneto-electro-elastic bi-material is investigated. The dislocation line is perpendicular to the isotropic basal plane of the bi-material. The elastic and electromagnetic fields induced by the dislocation are obtained through the use of the complex variable method together with the superposition scheme. The stress, electric displacement and magnetic intensity factors as well as the image exerted on the dislocation are given explicitly. We find that the intensity factors are expressed in terms of the so-called effective materials and the radial component of the image force is only dependent on the elastic modulus of the material with the dislocation. As an illustrative example, the bi-material that consists of piezoelectric and piezomagnetic phases is analyzed.     

Dislocation microstructures and strain-gradient plasticity with one active slip plane

We study dislocation networks in the plane using the vectorial phase-field model introduced by Ortiz and coworkers, in the limit of small lattice spacing. We show that, in a scaling regime where the total length of the dislocations is large, the phase field model reduces to a simpler model of the strain-gradient type. The limiting model contains a term describing the three-dimensional elastic energy and a strain-gradient term describing the energy of the geometrically necessary dislocations, characterized by the tangential gradient of the slip. The energy density appearing in the strain-gradient term is determined by the solution of a cell problem, which depends on the line tension energy of dislocations. In the case of cubic crystals with isotropic elasticity our model shows that complex microstructures may form in which dislocations with different Burgers vector and orientation react with each other to reduce the total self-energy.     

Charged dislocations in piezoelectric bimaterials

In some piezoelectric semiconductors and ceramic materials, dislocations can be electrically active and could be even highly charged. However, the impact of dislocation charges on the strain and electric fields in piezoelectric and layered structures has not been presently understood. Thus, in this paper, we develop, for the first time, a charged three-dimensional dislocation loop model in an anisotropic piezoelectric bimaterial space to study the physical and mechanical characteristics which are essential to the design of novel layered structures. We first develop the analytical model based on which a line-integral solution can be derived for the coupled elastic and electric fields induced by an arbitrarily shaped and charged three-dimensional dislocation loop. As numerical examples, we apply our solutions to the typical piezoelectric AlGaN/GaN bimaterial to analyze the fields induced by charged square and elliptic dislocation loops. Our numerical results show that, except for the induced elastic (mechanical) displacement, charges along the dislocation loop could substantially perturb other induced fields. In other words, charges on the dislocation loop could significantly affect the traditional dislocation-induced stress/strain, electric displacement, and polarization fields in piezoelectric bimaterials.     

Two-dimensional Green's functions in anisotropic multiferroic bimaterials with a viscous interface

We derive, by virtue of the unified Stroh formalism, the extremely concise and elegant solutions for two-dimensional and (quasi-static) time-dependent Green's functions in anisotropic magnetoelectroelastic multiferroic bimaterials with a viscous interface subjected to an extended line force and an extended line dislocation located in the upper half-plane. It is found for the first time that, in the multiferroic bimaterial Green's functions, there are 25 static image singularities and 50 moving image singularities in the form of the extended line force and extended line dislocation in the upper or lower half-plane. It is further observed that, as time evolves, the moving image singularities, which originate from the locations of the static image singularities, will move further away from the viscous interface with explicit time-dependent locations. Moreover, explicit expression of the time-dependent image force on the extended line dislocation due to its interaction with the viscous interface is derived, which is also valid for mathematically degenerate materials. Several special cases are discussed in detail for the image force expression to illustrate the influence of the viscous interface on the mobility of the extended line dislocation, and various interesting features are observed. These Green's functions can not only be directly applied to the study of dislocation mobility in the novel multiferroic bimaterials, they can also be utilized as kernel functions in a boundary integral formulation to investigate more complicated boundary value problems where multiferroic materials/composites are involved.     

Thermodynamically consistent phase field approach to dislocation evolution at small and large strains

A thermodynamically consistent, large strain phase field approach to dislocation nucleation and evolution at the nanoscale is developed. Each dislocation is defined by an order parameter, which determines the magnitude of the Burgers vector for the given slip planes and directions. The kinematics is based on the multiplicative decomposition of the deformation gradient into elastic and plastic contributions. The relationship between the rates of the plastic deformation gradient and the order parameters is consistent with phenomenological crystal plasticity. Thermodynamic and stability conditions for homogeneous states are formulated and satisfied by the proper choice of the Helmholtz free energy and the order parameter dependence on the Burgers vector. They allow us to reproduce desired lattice instability conditions and a stress-order parameter curve, as well as to obtain a stress-independent equilibrium Burgers vector and to avoid artificial dissipation during elastic deformation. The Ginzburg–Landau equations are obtained as the linear kinetic relations between the rate of change of the order parameters and the conjugate thermodynamic driving forces. A crystalline energy coefficient for dislocations is defined as a periodic step-wise function of the coordinate along the normal to the slip plane, which provides an energy barrier normal to the slip plane and determines the desired, mesh-independent height of the dislocation bands for any slip system orientation. Gradient energy contains an additional term, which excludes the localization of a dislocation within a height smaller than the prescribed height, but it does not produce artificial interface energy. An additional energy term is introduced that penalizes the interaction of different dislocations at the same point. Non-periodic boundary conditions for dislocations are introduced which include the change of the surface energy due to the exit of dislocations from the crystal. Obtained kinematics, thermodynamics, and kinetics of dislocations at large strains are simplified for small strains and rotations, as well.     

A cracked sliding interface between anisotropic bimaterials

A general solution for the stresses and displacements of a cracked sliding interface between anisotropic bimaterials subjected to uniform tensile stress at infinity is given by using the Stroh’s formulation. Horizontal and vertical opening displacements on the interface, stress intensity factors, and energy release rate are expressed in real form, which are valid for any kind of anisotropic materials including the degenerate materials such as isotropic materials. It is observed that stresses exhibit the traditional inverse square root singularities near the crack tips, and the vertical opening displacement and energy release rate are intimately related to a real parameter λ determined by the elastic constants of the anisotropic bimaterials.     

Interaction energy of interface dislocation loops in piezoelectric bi-crystals

Interface dislocations may dramatically change the electric properties,such as polarization,of the piezoelectric crystals.In this paper,we study the linear interactions of two interface dislocation loops with arbitrary shape in generally anisotropic piezoelectric bi-crystals.A simple formula for calculating the interaction energy of the interface dislocation loops is derived and given by a double line integral along two closed dislocation curves.Particularly,interactions between two straight segments of the interface dislocations are solved analytically,which can be applied to approximate any curved loop so that an analytical solution can be also achieved.Numerical results show the influence of the bi-crystal interface as well as the material orientation on the interaction of interface dislocation loops.