Monte Carlo of Solid Amorphous
Monte Carlo of Solid Amorphous
Monte Carlo of Solid Amorphous
13
2 The Model
without stress
∆H0
with stress
Energy
σVa
with the activated STZ volume Va . The relation is depicted in fig. 2.1
The local stress σ is a combination of the external stress σext , which we
assume homogeneous is space, and the internal stress σint (eq. (2.3)). The
internal stresses arise from heterogeneous deformation associated with STZ
events and will be discussed in the following.
14
2.2 The Elastic Stress Field
a single STZ event, the shear stress field of the dislocation array representing
the deformation is computed. In the following only one of the four dislocations
located at 3, 6, 9 and 12 o’clock positions will be examined in detail as the stress
field of the others can be obtained by rotation and shifting. The dislocation at
3 o’clock is chosen for further examination.
The shear stress field σ(x, y) of a single edge dislocation can be evaluated
with eq. (2.4), where x and y are the distance to the dislocation and D is a
constant containing the material parameters shear modulus G, Burgers vector
b and Poison’s ratio ν[24, 25].
x x2 − y 2
σ(x, y) = D (2.4)
(x2 + y 2 )2
Gb
with D = (2.5)
2π (1 − v)
Due to bi-periodic boundary conditions of the simulation, an infinite number of
equal dislocations in y-direction and in x-direction can be found. To account
for the L-periodicity in y-direction the analytical solution for dislocation arrays
is used. The total shear stress field of the array is obtained by summation of
the individual dislocations stress fields (eqs. (2.6) to (2.7))[24].
∞ x x2 − (y − kL)2
Gb X
σ(x, y) = (2.6)
2π (1 − v) k=−∞ x2 + (y − kL)2 2
Gb x b) cos(yb) − 1)
b (cosh(x
= (2.7)
2L (1 − v) (cosh(xb) − cos(yb))2
2πx 2πy
with x
b= , yb =
L L
15
2 The Model
∆Hij0 − σij Va
!
νij = ν0 exp − (2.8)
kT
(2.9)
16
2.4 Simulation Parameters
σ=σ
b Cσ (2.11)
= b C (2.12)
t = tbCt (2.13)
b CG
G=G (2.14)
The stress free energy barrier of each simulation cell Hij0 is assigned using
uniform distributed random numbers fij with expectancy hf i = 1 and average
barrier height H.
Va
Cσ = is given as a fraction of two material properties. Inserting Cσ in P ,
H
the inverse temperature can also be written as
Cσ V a
P = . (2.18)
kT
Another parameter, the stress redistribution parameter C, results from the
redistribution of the internal stresses. As proposed by Bulatov and Argon [23]
17
2 The Model
the Greens’s function G is used for solving the surface integral arising from the
conserving stress redistribution of internal stresses σij in eq. (2.19). The indices
k and l determine the position of the activated cell.
ZZ
σij = σext + C CG ∆γ Gi−k,j−l dA (2.19)
1
ZZ
= σext + ∆γ Gi−k,j−l dA (2.20)
C
1 X
= σext + ∆γ Gij,kl (2.21)
C ij,kl
• Inverse temperature P
18
2.5 Material Properties
Table 2.1: Estimated material properties for the determination of the simulations
parameters.
19