3D DISCRETE DISLOCATION DYNAMICS APPLIED TO ... - NUMODIS
3D DISCRETE DISLOCATION DYNAMICS APPLIED TO ... - NUMODIS
3D DISCRETE DISLOCATION DYNAMICS APPLIED TO ... - NUMODIS
Create successful ePaper yourself
Turn your PDF publications into a flip-book with our unique Google optimized e-Paper software.
2.2 Computation of stresses and displacements of dislocations 21<br />
b<br />
C<br />
A<br />
Field point<br />
Ω<br />
R<br />
dl’<br />
n Slip plane normal<br />
Dislocation<br />
loop<br />
b<br />
Triangular loop<br />
n<br />
Dislocation<br />
segments<br />
Figure 2.8: The parameters in the Burgers equation (Eq. 2.8) and decomposition of a dislocation<br />
loop by triangular dislocation loops<br />
[Hirth & Lothe 92]). Special care, however, should be taken at evaluating the inverse trigono-<br />
metric functions, as the author experienced. Barnett ([Barnett 85]) has developed a formula more<br />
suitable for numerical computation, which will be detailed below.<br />
The displacement at a field point P(r) generated by a triangular dislocation loop with points<br />
A(rA), B(rB) and C(rC) are expressed as Eq. 2.10. The triangular dislocation loop ABC and a<br />
field point are shown in Fig. 2.9.<br />
u(r) = − b<br />
4π Ω + FAB + FBC + FCA<br />
(2.10)<br />
Ω is the solid angle associated with the triangle ABC, which generates a discontinuity of ∆u = b<br />
in traversing the cut surface ABC. F ij(i,j=A,B or C) is a continuous displacement field term ex-<br />
cept on the dislocation line. The solid angle Ω and the continuous terms Fij are given as follows<br />
([Barnett 85]).<br />
<br />
<br />
s<br />
<br />
s − a s − b s − c<br />
Ω = −sign (Ri.n) 4 arctan tan tan tan tan<br />
2 2<br />
2<br />
2<br />
<br />
Fij = −<br />
1 − 2ν<br />
8π(1 − ν) (b × tij) ln Rj + Rj.tij<br />
Ri + Ri.tij<br />
+<br />
1<br />
8π(1 − ν) (b.nij)<br />
<br />
Rj<br />
Rj<br />
The vectors and the constants in Eq. 2.11 and Eq. 2.12 are listed below.<br />
− Ri<br />
<br />
× nij<br />
Ri<br />
(2.11)<br />
(2.12)