Категория: Математика

Решение методом изображений для краевой дислокации и цилиндрической полости в кристаллических твердых телах

Научный труд разместил:
11 сентября 2020
Автор: Kien Nguyen

УДК 539.3

Решение методом изображений для краевой дислокации и цилиндрической полости в кристаллических твердых телах

K. Nguyen, A. Mehrabian

Механика дефектов в твердых телах в широком диапазоне масштабов длины обычно формулируется в терминах теории дислокаций. В настоящей статье рассмотрена классическая задача взаимодействия упругой краевой дислокации с цилиндрической полостью. Решение данной задачи в напряжениях получено в рамках эвристи-ческо-механистического подхода, в котором используется теория функций комплексного переменного в теории упругости и метод изображений. Для этого определены и описаны сингулярности типа упругого диполя, подобные дипольным зарядам в электростатике. Показано, что виртуальная дислокация с вектором Бюргерса такой же величины, что и реальная дислокация, но противоположно направленным, а также четыре сингулярности, включая дислокационный диполь, моментно-дилатационный диполь и два центра дилатации, могут определять круговую свободную границу в бесконечной упругой среде. Добавление дислокации Вольтерра к краевой дислокации конечной длины позволяет решить связанную задачу взаимодействия краевой дислокации бесконечной длины с цилиндрической полостью. Обсуждается аналогия между рассматриваемой упругой задачей и электростатической задачей взаимодействия линейного электрического заряда с цилиндрическим проводником.

DOI 10.24411/1683-805X-2020-14005

Method of images solution for an edge dislocation

and a circular cavity in crystalline solids

K. Nguyen and A. Mehrabian

Department of Energy and Mineral Engineering, Energy and Mineral Sciences Energy Institute, The Pennsylvania State University, PA 16802, USA

Mechanics of defects in solids across a wide span of length scales is commonly formulated using the dislocations theory. This paper revisits the classical problem of interaction between an elastic edge dislocation and a circular cavity. A heuristic, yet, mechanistic approach is taken to obtain the stress solution to this problem. The approach uses complex variable theory of elasticity, along with method of images. For this purpose, a definition and formulation of elastic dipole singularities similar to dipole charges in electrostatics is developed. It is shown that an image dislocation with Burger&s vector of the same strength as the real dislocation but in opposite direction, as well as a set of four singularities including a dislocation dipole, a moment-dilatation dipole, and two centers of dilatation would establish a circular, traction-free boundary in an infinite elastic medium. Adding a Volterra dislocation to the finite-length edge dislocation from this study would recover the related problem of interaction between an infinite-length edge dislocation and circular cavity. The interesting analogy between the considered elastic problem and the electrostatic problem of interaction between a line electric charge and a cylindrical conductor is discussed.

Университет штата Пенсильвания, Пенсильвания, 16802, США


q—magnitude of electrical charge,

q0—magnitude of image electrical charge,

V—electrical potential,

x—real axis in complex plane, rectangular x coordinate axis,

y—imaginary axis in complex plane, rectangular y coordinate axis,

R0—a reference radial distance,

© Nguyen K., Mehrabian A., 2020

R—circular cavity radius centered at the origin, —the end-coordinate of a finite edge dislocation measured from the origin,

0—angle measured from x-axis, G—shear modulus,

v—Poisson ratio,

<5ij—stress tensor,

V2—Laplace operator,

U—Airy&s stress function,

z—complex coordinate of physical plane,

Q—complex stress potential,

a—complex stress potential,

Re—real part of a complex expression,

Im—imaginary part of a complex expression,

b—Burger&s dislocation vector in complex plane,

bx—glide edge dislocation magnitude,

by—climb edge dislocation magnitude,

1—complex unit number, i = (-1)12,

e—Euler&s number, natural logarithm base constant, ^—generic auxiliary potential,

C—generic auxiliary location of singularity in complex plane,

J—generic auxiliary magnitude of singularity, B—auxiliary dislocation magnitude, Q—auxiliary center of dilatation magnitude, M—auxiliary moment magnitude, P—auxiliary combined moment and center of dilatation magnitude, S—auxiliary coefficient,

||—indicating the norm of a complex expression, Arg—indicating the argument of a complex expression.


r—radial direction in cylinder coordinate, centered at the origin,

0—tangential direction in cylinder coordinate,

i, j—stress tensors index,

x—real part of a complex variable,

y—imaginary part of a complex variable,

b—Burger&s dislocation vector in complex plane.


--overbar indicating complex conjugate property,

&—derivative in complex plane,

—dipole of a singularity, ed—edge dislocation, dp—dipole of a potential,

fed—indicating a property of a finite length edge dislocation,

V—indicating of the property of a Volterra dislocation.

1. Introduction

The method of images offers a rather exploratory means of solving certain problems in mathematical physics involving singularities within finite domains, e.g., in electrostatics [1], plane elasticity [2], heat conduction [3], flow of inviscid fluids [4] and theory of vibrations [5]. The technique involves expanding the solution domain to an infinite space by introducing an image of the real problem domain against one or more boundary surfaces or hypersurfaces. Fictitious singularities of deliberately selected strength and location are usually placed in the image domain in such a way that the prescribed boundary conditions of the problem boundary surfaces are secured.

A pioneering application of the images method is found in the Ref. [6] on the interaction between point charges and spherical conductors. The two-dimensional variation of the problem is shown in Fig. 1 where a line electric charge q located at x = an image charge -q located at x = R2/^, along with an auxiliary charge of strength q0 = -qln(R/^)/ln(R0/R) located at x = 0, within an infinite space would produce a circular equipotential contour of radius R centered at the origin of the coordinate system [1]. Thus, the uniqueness theorem of Poisson&s equation would immediately establish the solution to the problem of electrostatic interaction between an infinite line charge q located at a distance from the center of a circular conductor held at electric potential V = 0 [1].

Method of images has been applied to problems of defects in elastic crystalline solids, in particular, via the theories of dislocations and fracture mechanics. Early contributions to the subject considered the analogy between elastic stress fields of screw dislocations and electrostatic fields of line charges [7, 8]. The interaction between an elastic edge dislocation and a planar boundary at the interface of two dissimilar solids using the images method is presented in Ref. [9]. The solutions for the elastic interaction between screw and edge dislocations with boundaries of inhomoge-neities were further treated in Ref. [10]. Extension of the solution in Ref. [7] was later applied to semiconductors technology where the configuration of inho-mogeneous strained layer superlattices was shown to allow for filtering screw dislocations in gallium arsenide layers grown on silicon [11]. The phenomenon was reported to be responsible for superior yield

Fig. 1. Schematics of electrostatic problem of a line charge interacting with a conducting cylinder (a), as well as the images charges at image point, x = R2/E, (b), and the conductor center, x = 0 (c). The arrowed curves and black curves are the resulting electric field and equipotential contours (d) (color online)

strength of lamellar materials compared to the bulk pure crystals of their either constituents. A scaling theory of such strengthening effect based on the method-of-image solution of the repulsive forces between dislocation rows and inhomogeneity interfaces between alternating layers of lamellar materials is given in Ref. [12]. More recently, the displacement profile of the interface for the inhomogeneity-dislocation problem was published [13] while multilayer extensions of the solution using successive dislocations images can be found in Refs. [14-18]. The method of images has been used to solve for the elastic field of twist disclinations in nonhomogeneous infinite spaces [19] and in hexagonal crystals [20].

A similar record of literature has documented the fundamental and applied solutions on the interaction between dislocations and circular boundaries. Eshelby applied the images method to solve for the problem of a screw dislocation in the vicinity the circular boundary of a thin rod [21]. Recently, the work of Eshelby was extended to the cases of multiple annular elastic

layers [22] and piezoelectric nanowires [23]. The stress field of an edge dislocation in presence of a circular inhomogeneity is somewhat more complex than the alternative problem involving a screw dislocation since the former problem entails biharmonic stress potentials whereas the stress potential of the latter is a harmonic function. The pioneering work by Dundurs and Mura [24] concludes the solution to this problem. The solution in [24] is, however, a pure mathematical guess based on a rather try-and-error investigation of the similarities with the associated stress potentials from the limiting case of two inhomogeneous semiinfinite planes [25, 26]. Dundurs and Mura&s solution on the subject was later reformulated in the complex plane [27]. The solution for the case when the circular boundary is coated by a thin layer of different properties is presented in Ref. [28] whereas the case of two circular inclusions is treated in Ref. [29]. The problem of interaction between an edge dislocation and an elliptic boundary is successfully solved by conformal mapping technique in Ref. [30]. The solution for the

case of edge dislocation inside an elliptical inclusion is published in Ref. [31]. The solution for an edge dislocation interacting with an interfacial crack at the boundary of a circular inhomogeneity is determined in Ref. [32]. The case of imperfect boundary interface between the matrix and circular inhomogeneity is investigated in Ref. [33]. More recently, the solution to the problem of cavity and dislocation is extended to the case where the matrix is an elastic half-space [34].

This paper uses the electrostatics analogy of images method to revisit and solve the problem of interaction between an edge dislocation and a circular cavity. The nonzero stresses due to the real and image dislocations of opposite sign at the boundary of the cavity are determined and inspected. It is shown that deliberate placement of four auxilliary singularities of determined strength would retrieve traction-free condition of the cavity surface. These mechanical singularities include a dislocation dipole, a moment-dilatation dipole and two centers of dilatation. The strengths of these singularities are obtained from a system of linear algebraic equations that are extracted from the traction-free boundary conditions on the cavity surface. Results are verified by the previous published solution to the same problem [24].

2. The elastic problem of edge dislocation and circular cavity

Figure 2 illustrates the plane-strain problem where a single edge dislocation of Burger&s vector b is placed in an infinite elastic medium of elastic con stants G, v located at a distance £ from the center of a traction-free circular cavity. The problem boundary conditions consist of zero traction on the cavity boundary, i.e.,

Grr(R, 9) = 0, (1)

MR, 9) = 0 (2)

together with zero far-field stresses, as follows:

Gj(r ^ «>, 0) = 0. (3)

The polar-coordinate form of Navier&s equations of static equilibrium, along with stress compatibility condition, takes the following forms [2]:

1 da

- + —

+ 00 _

1 da

+ da

where V2 = — dr

+ r 2a

V2(arr + a00) = 0,

is Laplacian operator

of polar coordinate.

The stress components can be given in terms of Airy stress function U(r, 9) [35]:

1 5U 1 d 2U a r =——+r dr r2 d02
1 dU

Direct substitution of the Eqs. (7)-(9) would verify that Airy&s stress function satisfies both the stress equilibrium and compatibility conditions expressed in Eqs. (4)-(6), provided that U(r, 9) is a biharmonic function, i.e., V2(V2U) = 0.

The Airy stress U can be written in terms of two constitutive complex potential functions Q(z) and ©(z), as follows [36]:

U( r, 9) = 2[ ZQ( z) + zQ( z) + ©( z) + ©( z)], (10)

where the in-plane stresses can be expressed as a combination of these two complex potentials as

a rr +099= 2[ Q&( z) + Q&( z)], (11)

Fig. 2. Schematics of a finite edge dislocation and a circular cavity in an infinite elastic domain

099 -Orr + 2iar9 = 23[zQff(z) + ©"(z)], (12) z

arr -iar9 =Q&(z) + Q&(z) -z[zQ"(z) + (z)].(13)

In polar coordinate, the in-plane stresses are calculated from Eqs. (11)—(13), as follows:

a rr (r, 0) = Re a06 ( r, 0) = Re

2Q&( z ) - zQ"( z ) -3 z )
2Q&( z ) + z Q"( z ) + 3w"( z )

a r6 ( r, 0) = Im

zQ" ( z ) + ( z )

2.1. Method of images solution

Without loss of generality, the real axis of the complex plane with the origin at the center of the cavity in Fig. 2 is assumed to align with the line connecting the dislocation and cavity center. For a dislocation with Burger&s vector b located at z = S, in an infinite elastic field, the complex potential functions take on the following form [36, 37]:

Qf(z, S) = Gb >( z(17)

4ra(1 -v)

(*, S) =

4ra&(1 -v)

-b ln(* -S) bS *-S

Kelvin&s method of images [6] for circular boundaries is followed by placing an image edge dislocation with Burger&s vector -b at * = R2/S. This method follows a direct analogy of the well-known electrostatics solution for interaction between a line charge and a circular conductor, as shown in Fig. 1. The complex potential functions pertaining to the real dislocation at * = S and image dislocation at * = R2/S are written as

q( *, S)=Qbd( *, S)+Q:


«&( *, S) = ®bed( *, S) +


where Qed and ra&ed are defined in Eqs. (17) and (18). By substituting Eqs. (19) and (20) in Eq. (13), the inplane shear and radial stresses at the hole&s boundary z = Re&e simplify to

4ra(1 -v)


( Re/0-S)2 ( Re,9-R2/S)2

4ra(1 -v)
4ra(1 -v)

Re-9 - S Re~,9-R2/s bs bR2

( Re,9-S)2 S( Re,9-R2/S)2

Re/9- R2/S Re,9-S

The correct solution to the considered problem requires the right side of Eq. (21) to vanish. Hence, the

task at hand becomes a search for supplementary singularities negating the right side of Eq. (21).

The mathematical form of stresses expression in Eq. (21) suggests an ad hoc method to secure this condition. It appears that singularities of higher order can be used for this purpose. Dipoles of elastic singularities could form the basis of a plausible solution. Appendix A presents a formulation for potential functions of the elastic dipoles, as well as the corresponding formulae for known elastic singularities including dislocations, moments, or centers of dilatation. The generalized formulations of singularities are summarized in Tables A1 and A2. In particular, from Appendix A, the complex potentials of a dislocation dipole with Burgers vector B>1 at R2/S, are given by substituting Eqs. (17) and (18) in Eq. (A1):

4ni(1 -v) * - R2/S:
2 A
4ra(1 -v)

B B r2Is

- r2¡s * - r2/s (* - rvs)2

where Bx is the unknown Burgers vector of the dislocation dipole at the image point. An attempt by introducing a dislocation dipole into Eq. (21) reveals that it is impossible to completely negate the nonzero terms on the right side of Eq. (21) with a dislocation dipole alone. Thus, auxiliary complex potentials that, likewise, contains second order singularity originating from a combined moment-center of dilatation at the image location are used as well. From Eqs. (A8) and (A9), the potential functions pertaining to a combined singularity of the form P1 = M1 + iQ1 located at R2/S, in a complex plane, would find the following form:

4^/(1 -v) * - R2/S:

where P1 is the complex variable of the combined singularity, M1 and Q1 are real-valued magnitude of moment and center of dilatation at the image location. P1 = -P(1 - v)/G is defined to maintain the consistency with the generic definition in Appendix A. Appendix B shows that dislocation dipole B>1 induces third order singularity to the problem, which by itself does not cancel out with any other terms. As a result, a combined moment-center of dilatation dipole is also placed at the image location, whose magnitude P2 is determined in such a way to reduce the third order

singularity arising from the previously selected dislocation dipole. Following the derivations in Appendix B, P2 and B1 are related by

By adding the combined moment-center of dilatation dipole P2 and by substituting Eq. (26) into the total potential, the zero in-plane shear stress at the circular cavity&s wall suggests an approach to determine the magnitude of B1 and P1. Since B1 and P1 are complex variables, they can be written as B1 = B1x + iB1 y and Pi = M1 + iQ1. Thus, the complex expanded form of Eq. (16) determines the shear stress around the circular hole containing the physical and image dislocations, one dislocation dipole BB1, one moment-center of dilatation dipole P2 and one moment-center of dilatation P1. The resulting mathematical expression would find the following form:

0r9 (z = Re*)

= 4rc(1 -v)R3(R2 + £2 -2R£cos9)2

x [4bxR5£ - 4bxR3£3 - 4B^R£4 - MR?

+ 2(bx (-R 6 + R 2£4) + £3(M1R 2 + B1x (R2 + £2))) cos 9 - M1R3£2 cos(29)

sin 0(-2byR6 + 2B1yR2f3 - 2QR2f3

l2>-4 , e5\\ , • r>5t

&y" & ~ ly ■ 2byRL 2 B?1V S5) + sin(20)(2bvR f

bx ( - R6 + R 2f4) + f3 (M1R 2 + Bu ( R2 +S2)) = 0,

-2byR6 + 2BB1yR2f3 - 2Q1R2f3

+ 2byR 2f4 + 2 B^5 = 0,

2byR5£ - 2 B yR3£ 2 + £R3£ 2 - 2byR3£3 = 0.

Equations (28) constitute a set of 5 linear algebraic equations to be solved for only 4 variables. Although there are more equations than variables, the solution can be uniquely determined, as follows:

B = -R 2(l23- R 2) b,

P = R2(R2-S2)2 b P2~ b,

Q1 = 2Im[b]

S 4 - R4

- 2B1 yR3S2 + Q^2 - 2byR3S3)]. (27)

The vanishing shear stress condition at the interface demands that all groups of trigonometric functions should identically vanish. The resulting set of linear algebraic equations is outlined, as follows:

M1R 3S2 = 0, 4bxR5S - 4bxRY - 4B^4 -M^4 = 0,

Although the in-plane shear stress condition in Eq. (2) is met, the radial stress at the boundary is yet to be satisfied. The residual radial stress at the wall&s boundary after B1, P2 and P1 having been determined is

9 2bvG(£2 - R2)

a rr (z = Re&9) = --1. (30)

4rc(1 -v )R2£

The constant residual radial stress in the elastic problem is similar to the electrostatic subproblem in Fig. 1, b where a constant potential is found on the interface of a conducting cylinder. Thus from Eq. (30) it is straight-forward to include a uniquely defined center of dilatation at the cavity center to negate the obtained radial stress. From Appendix A, the complex potential of a dilatation of magnitude Q3 at the cavity center is found by

qq3(z, 0) = 0, (31)

®Q (z, 0) =-Q3, (32)

Q3 4n(1 -v)R2£ z

where Q3 is the strength of dilatation center at the origin. Using (14), (30) and (32), Q3 is determined as

Q3 = -2Im[b]

3. Discussion

In summary, the exact analytical solution of the problem of a single dislocation b placed at an arbitrary distance £ away from the center of circular hole can be decomposed into two subproblems in an infinite continuum domain: (i) the dislocation b at £ and (ii) complementary terms to account for the stress free boundary condition at the cavity boundary. The latter second problem can be further written as a superposition of five singularities in infinite elastic domain. Figure 3 illustrates the summary of these elastic singularities.

Consequently, the complex holomorphic functions describing the interaction of a dislocation of finite length (£ - R) with the circular cavity are obtained, as follows:

Fig. 3. Schematics of the method-of-images solution to the dislocation and cavity problem in an infinite elastic domain

Q( *, S) = Qf( *, S) + Ql

ra&( *, S) = rabed( *, S) + ra-eb

+ Qo(*), (34)

+ ®0( z) + ®1( z ), (35)

where the generic primary complex potentials of a dislocation in an infinite elastic field Qbed(z, S), Q_bed(z, R2/S), ro&bed(z, S), and ro&-bed(z, R2/S) are obtained from Eqs. (14) and (15), Q0(z), ro&0(z) and ro&i(z) are extra singular expressions due to the center of dilation and moment-dilatation dipole at the image point, as well as the center of dilatation at the cavity center. These potentials are obtained as

4ni(1 -v) S3(* - R2IS)

Gb R2(S2 - R2) 4 ni (1 -v ) S 2( * - R 2/S)2


G Im[b] 2(S2 - R2Y

4n(1 -v ) S

The table shows the striking analogy between the electrostatic and elastostatic problems. The associated method-of-images solutions to both problems involve placing an image singularity at * = R2/S, as well as an auxiliary singularity at the cavity center. However, a set of three extra auxiliary singularities at the image

The analogy between method-of-images solutions of the electrostatic and elastostatic problems with circular boundaries

Electrostatics Elastostatic s

Potential function V2V = 0 V2(V2U) = 0

Boundary conditions at z = Re&e V = 0 S ( 1 SU A_ 0 1 SU + 1 S2U _ 0 8r t r se J_ & r Sr r2 se2 _

Singularity location Description Magnitude Description Magnitude

Singularity at z = E, Electric charge q Dislocation b

Primary image singularity at z = R2/^ Image electric charge -q Image dislocation -b

Auxiliary singularities at z = R2/^ None 0 Dislocation dipole R2(£2 -R2) F £3

Moment-dilatation dipole bR2(£2 -R2)2 £4

Center of dilatation 2Im(b)(£4 -R4) £3

Auxiliary singularities at z = 0 Electric charge q ln(R £) ln(R0/ R) Center of dilatation 2Im(b)(£2 - R2) £

location is required in the case of the elastostatic problem to reproduce the traction-free boundary conditions on the cavity.

3.1. Verification with results of Ref. [24]

From Eqs. (34) to (38), Airy stress functions pertaining to glide bx and climb by components of the edge dislocation b = bx + iby in presence of a circular cavity can be derived by using Eq. (10). The results are expressed as

4rc(1 -V!)

- 2| z-£ |ln| z-£|sin(Arg(z-£))

(£2 - R2)R2 (

2 sin

in [Arg (z-R7£)

Uy _-Ly— <

y 4rc(1 -v1)

2|z-£|ln| z - £ | cos (Arg(z-£))
2 V

(£2 - R2)R

2\\ r>2 ( -,i-2

V ^ /. 2


[Arg( z - R2/£) ]

£2 - R2 £

£2 - R2

- 2^-— ln | z |

The infinite edge dislocation from Dundurs and Mura&s formulation [24] can be retrieved from the finite-length edge dislocation solution of this study if Eqs. (39) and (40) are superposed with a Volterra dislocation [5], whose potentials are written as

QV( z ) =

4ra&(1 -V)

Fig. 4. Schematics of a ring-shaped Volterra dislocation

®bv( z ) = -G

b (ln z +1) +

4rci(1 - v)

Hence, the Airy stress of the Volterra dislocation becomes:

UV( z ) = Re

--z (ln z + ln z ) +--i i z

Equation (43) pertains to the special case of a more generalized Volterra dislocation in Eq. (44) of a ring-shaped elastic field given by Lardner [38], whose schematic is displayed in Fig. 4:

U V( z ) = Re

— z (ln z + ln z ) +

bR 2 +-i z

It is easy to see that as R1 ^ m, the limiting case from Eq. (44) will be simplified into Eq. (43). Superposing Volterra&s dislocation stress function by Eq. (43) on the glide and climb components of the Burger&s vector bx and by with potentials Q(z) and ro&(z) defined by Eqs. (34)-(38) and rederiving the real-valued Airy stress function from Eq. (10), the following identical formulation of Dundurs and Mura for an infinite-length dislocation is recovered [24]:

x 4<1 -v1)

GA J -2ri ln r sin 01

+ 2(r2 ln r2 sin 62 - r ln r sin 6) sin(262)

(^2 - R2)R2

(^2 - R2) sin 62

2 sin 6

Uy = /. n y \\ i 2riln ricos 6i 4^(1 -vi) [

- 2(r2 ln r2 sin 62 - rln r cos 6)

(^2 - R2)R2


- cos(262) +

(^2 - R2) cos 62

2(^2 - R2) ln r + R
2 cos 6
61 = Arg( z -Ç), 62 = Arg

r R2 >

9 = Arg(*).
4. Conclusion

Method of images is used to solve the problem of a circular cavity interacting with a finite edge dislocation. It is shown that aside from the image dislocation, a set of three singularities including a center of dilatation, a moment dilatation dipole and a dislocation dipole at the image point, as well as a center of dilation at the cavity center, would secure the traction-free boundary of the cavity. The presented solution closely mimics Kelvin&s solution of interacting line charges and cylindrical conductors in electrostatics.

Appendix A. Dipoles of elastic singularities

The complex potential of dipole of strength J is, defined by a pair of corresponding singularities, with the same magnitude of strength J but in opposite in direction, one located at an arbitrary point Ç and the other located at an infinitesimal distance Ç + dÇ in such a way that J = JdÇ is a finite value [39, 40]. With this definition, supposing the potential function <(*, Ç) with singularity&s magnitude J, the complex potentials of a dipole can then be obtained, as follows:

J (*, Ç) = lim (<J (*, Ç) + <_j (*, Ç + dÇ))

= J — [<Kz, C)] dC = J—[<Kz, C)]. (A1)

Similar to the electrostatic charge dipoles, neither J nor dC would have any physical meaning since the former tends to infinity while the latter is infinitesi-mally small. In this context, the dipole&s magnitude is written in terms of J, which has a finite magnitude. Hence, from Eqs. (17), (18) and (A1) the dipole of a dislocation, whose magnitude at C is the Burgers& vector B, becomes:

Table A1. Complex potential Q(z, C)

z, C ) = "


<p( z, c) =

4ra(1- v)
4ra(1- v) z- C

The complex potential of a moment and a center of dilatation are outlined in Eqs. (A4), (A5) and (A6), (A7), respectively [41]:

Qm (z, C) = 0,

(z, C) = 4- --, (A5)

Qe (z, C) = 0, (A6)

aQ (z, C) = - Q-L-, (A7)

where M and Q are real-valued magnitude of a moment and a center of dilatation, respectively. Since their forms are identical with the same order of singularity at point C, it is convenient to define a complex variable combining moments and centers of dilatation together, denoted as P = M + iQ. Thus, with this definition, the combined complex potential becomes:

Qp (z, C) = 0, (A8)

®P (z, C) = ----. (A9)

Using the definition of dipoles (A1), a combined moment-center of dilatation dipole of magnitude P = M+iQ located at C in a two-dimensional plane is defined by the potentials:

^Pp( z, C) = 0,

4^ (z -C)2&

A summary of singularities and their associated potentials is tabulated in Tables A1 and A2.

ln(z - C) 1/(z - C)

Dislocation B GB 4ra(1 - v) 0

Moment M 0 0

Center of dilatation Q 0 0

Dislocation dipole B 0 gB 4ro(1 - v)

Moment dipole M 0 0

Center of dilatation dipole Q 0 0

Table A2. Complex potential ra&(z, C)

ln(z - C) 1/(z - C) 1/(z - C)2

Dislocation B GB GB 0

4ro(1 -v) 4ra(1 - v)

Moment M 0 iM 0

Center of dilatation Q 0 0

Dislocation dipole B 0 G(B -B) 4ro(1 -v) gB 4ro(1 - v)

Moment dipole M 0 0 iM

Center of dilatation dipole Q 0 0

Appendix B. Reduction of the third-order singularity

Adding a dislocation dipole B>1 and a combined moment-center of dilatation P1 to Eq. (21), the boundary stress terms, i.e., z = Re-e would find the following form.

arr - -cre =

4ra(1 -v)

Rei0-£ ReiQ- R2

(Rei0 - R V£)2 (Reid - £)2 (Rei0 - R 2/£)2

2 ReB B

(Re,Q- R2/^)3

4ra(1 -v)

Re-,e-S V ^ Re-,e- R 2 ¡S & (Re-,e- R2/S)2

Ge2,e f bS bR 2

4ra(1 -v) V ( Re,e-S)2 S(Re,e - R2/S)2

& Re,e- R2/S Re,e-S (Re,e -R2/S)2

B1 , 2B1R 2

(Re,e - R 2/S)2 S(Re,e - R 2/S)3

(Re,e - R 2/S)2

Equation (B1) shows that the dislocation dipole at the image point induces stress singularity of third order. Thus, to negate this third order singularity, a combined moment-center of dilatation dipole is placed at the image point location R2/£, so that it is reduced to the second order- To be consistent with the notations of Appendix A, the combined potential magnitude is equivalent to P2 = -(1 - v)/GPP:

4ra(1 -v) (z - R2/S)2

where P2 = M2

- iQ2 is the complex coefficient of the moment-center of dilatation dipole at the image point. Adding only the third order singularities of stresses induced by dislocation dipole B1 and the combined moment-center of dilatation dipole P2, the relationship between their magnitudes is determined from substitute into latter part of Eq. (13) at the circular cavity&s boundary, i.e.

-[ zQ" (z) + (z)]

4ra(1 -v)
2R 2e2,e B1

( z=Re )

2Re,e B1 (Re,e - R2/S)3
2e2,e P2

S(Re,e - R2/S)3 (Re,e - R2/S)3

BB1( R 2e,e- RS) -Se,e P2

B1( R 2e,e- RS) -Se,e Js2 = s

( R 2 > Re,eS

where S is any linear combination of products containing powers of e&e, R and Then, P2 is related to B1 as

4ra(1 -v) ^(Rei9- R2/^)3

To reduce the third order singularity to second order, the numerator must reduce to a common factor of (Reie - R2/^), hence:


1. Jackson J.D. Classical Electrodynamics. - John Wiley and Sons, 2007.
2. Barber J.R. Elasticity. - Dordrecht: Kluwer Academic Publishers, 2002.
3. Carslaw H.S., Jaeger J.C. Conduction of Heat in Solids. - Clarendon Press, 1992.
4. Munson B.R., Okiishi T.H., Huebsch W.W., Rothmayer A.P. Fluid Mechanics. - Singapore: Wiley, 2013.
5. Volterra V. Note on the application of the method of images to problems of vibrations // Proc. Lond. Math. Soc. - 1905. - V. 2. - No. 1. - P. 327-331.
6. Maxwell J.C. A Treatise on Electricity and Magnetism. V. 1. - Clarendon Press, 1873.
7. Eshelby J.D. The force on an elastic singularity // Philos. Trans. Roy Soc. Lond. A. Math. Phys. Sci. -1951. - V. 244. - No. 877. - P. 87-112.
8. HeadA.K. The interaction of dislocations and boundaries // Philos. Mag. J. Sci. - 1953. - V. 44. - No. 348. -P. 92-94.
9. Head A.K. Edge dislocations in inhomogeneous media // Proc. Phys. Soc. B. - 1953. - V. 66. - No. 9. -P. 793.
10. Dundurs J. Elastic Interaction of Dislocations with In-homogeneities // Mathematical Theory of Dislocations. - 1969. - P. 70-115.
11. Taylor R.I. The force on a screw dislocation due to a series of layers of alternating shear modulus // Semicon-duct. Sci. Technol. - 1989. - V. 4. - No. 8. - P. 612.
12. Friedman L.H., Chrzan D.C. Scaling theory of the Hall-Petch relation for multilayers // Phys. Rev. Lett. -1998. - V. 81. - No. 13. - P. 2715.
13. Chou Y.T., Pande C.S., Masumura R.A. The role of harmonic functions in dislocation-boundary interactions by the method of images // Mater. Sci. Eng. A. -2007. - V. 452. - P. 99-102.
14. Ma C.C., Lu H.T. Theoretical analysis of screw dislocations and image forces in anisotropic multilayered media // Phys. Rev. B. - 2006. - V. 73. - No. 14. -P. 144102.
15. Wen J., Wu M.S. Analysis of a line defect in a multilay-ered smart structure by the image method // Mech. Mater. - 2007. - V. 39. - No. 2. - P. 126-144.
16. Zhou K., Wu M.S. Elastic fields due to an edge dislocation in an isotropic film-substrate by the image method

// Acta Mech. - 2010. - V. 211. - No. 3-4. - P. 271292.

17. Ogbonna N. Force on a screw dislocation in a multiphase laminated structure // Math. Mech. Solids. -2014. - V. 19. - No. 6. - P. 694-702.
18. Ogbonna N. On screw dislocation in a multiphase lamellar inclusion // J. Niger. Math. Soc. - 2015. -V. 34. - No. 1. - P. 32-39.
19. Chou T.W. Elastic behavior of disclinations in nonho-mogenous media // J. Appl. Phys. - 1971. - V. 42. -No. 12. - P. 4931-4935.
20. Chou T.W., Pan Y.C. Elastic energies of disclinations in hexagonal crystals // J. Appl. Phys. - 1973. - V. 44. -No. 1. - P. 63-65.
21. Eshelby J.D. Screw dislocations in thin rods // J. Appl. Phys. - 1953. - V. 24. - No. 2. - P. 176-179.
22. Ogbonna N. On elastic interaction of a screw dislocation with a coated cylindrical inclusion // J. Eng. Math. - 2016. - V. 99. - No. 1. - P. 203-212.
23. Wang X., Pan E. Screw dislocations in piezoelectric nanowires // Mech. Res. Comm. - 2010. - V. 37. -No. 8. - P. 707-711.
24. Dundurs J., Mura T. Interaction between an edge dislocation and a circular inclusion // J. Mech. Phys. Solids. - 1964. - V. 12. - No. 3. - P. 177-189.
25. Dundurs J., Hetényi M. The elastic plane with a circular insert, loaded by a radial force // J. Appl. Mech. -1961. - V. 28. - No. 1. - P. 103.
26. Hetényi M., Dundurs J. The elastic plane with a circular insert, loaded by a tangentially directed force // J. Appl. Mech. - 1962. - V. 29. - No. 2. - P. 362.
27. List R.D. A two-dimensional circular inclusion problem // Math. Proc. Cambridge Philos. Soc. - 1969. -V. 65. - No. 3. - P. 823-830.
28. Povstenko Yu.Z. Interaction between an edge dislocation and a circular boundary in the presence of an alien surface layer // Sov. Appl. Mech. - 1975. - V. 11. -No. 3. - P. 272-277.
29. Fukuzaki K., Shioya S. On the interaction between an edge dislocation and two circular inclusions in an infi

nite medium // Int. J. Eng. Sci. - 1986. - V. 24. -No. 12. - P. 1771-1787.

30. Chen D.H. Green&s functions for a point force and dislocation outside an elliptic inclusion in plane elasticity // ZAMP. - 1996. - V. 47. - No. 6. - P. 894-905.
31. Warren W.E. The edge dislocation inside an elliptical inclusion // Mech. Mater. - 1983. - V. 2. - No. 4. -P. 319-330.
32. Fang Q.H., Liu Y.W., Jiang C.P. Edge dislocation interacting with an interfacial crack along a circular in-homogeneity // Int. J. Solids Struct. - 2003. - V. 40. -No. 21. - P. 5781-5797.
33. WangX. Interaction between an edge dislocation and a circular inclusion with an inhomogeneously imperfect interface // Mech. Res. Commun. - 2006. - V. 33. -No. 1. - P. 17-25.
34. Dai D.N. An edge dislocation inside a semi-infinite plane containing a circular hole // Int. J. Solids Struct. -2018. - V. 136. - P. 295-305.
35. Malvern L.E. Introduction to the Mechanics of a Continuous. - Englewood Clifs, N.J.: Prentice-Hall, 1969.
36. Muskhelishvili N.I. Some Basic Problems of the Mathematical Theory of Elasticity. - Graningen: Noord-hoff, 1954.
37. Ardakani S.M., Ulm F.J. Chemoelastic fracture mechanics model for cement sheath integrity // J. Eng. Mech. - 2014. - V. 140. - No. 4. - P. 04013009.
38. LardnerR.W. Mathematical Theory of Dislocations and Fracture. - Toronto: University of Toronto Press, 1974.
39. Chen Y.Z., Lin X.Y. Potentials in-plane elasticity by distribution of dislocation doublet or force doublet along a curve // Int. J. Eng. Sci. - 1998. - V. 36. -No. 1. - P. 23-31.
40. Denda M., Kosaka I. Dislocation and point-force-based approach to the special Green&s function BEM for elliptic hole and crack problems in two dimensions // Int. J. Numer. Meth. Eng. - 1997. - V. 40. - No. 15. -P. 2857-2889.
41. Green A.E., Zerna W. Theoretical Elasticity. - New York: Dover Publications, 2012.

Received 15.07.2020, revised 15.07.2020, accepted 31.07.2020

Сведения об авторах

Kien Nguyen, Graduate Res. Assist., The Pennsylvania State University, USA, ktn36@psu.edu

Amin Mehrabian, PhD, Assist. Prof., The Pennsylvania State University, USA, amin.mehrabian@psu.edu

метод изображений дислокация упругость механика дефектов механика твердого тела method of images dislocation elasticity mechanics of defects solid mechanics