LAMMPS WWW Site - LAMMPS Documentation - LAMMPS Mailing List Archives
[lammps-users] Dislocation-Stacking Fault Interaction
[Date Prev][Date Next][Thread Prev][Thread Next][Date Index][Thread Index]

[lammps-users] Dislocation-Stacking Fault Interaction


From: Rosita Masoudi <rositamasoudi@...24...>
Date: Sat, 25 Nov 2017 21:05:36 +0330

Dear All
Hi

I am trying to investigate the interaction of dislocation-stacking fault in FCC metal.
The model is visualized here:
https://ibb.co/m6W87m

In different articles for interaction of dislocation with other defects mentioned to fix the upper and lower parts, Then, give velocity to the middle part or shear stress with NPT command which all stress components are zero except xy one.

I test this for many inputs but stacking fault disappears or the HCP crystal structure of it change to amorphous and dislocation didn't move toward the stacking fault.

For first method my input is written bellow:
region lower block INF INF INF 1 INF INF units box
region upper block INF INF 1 INF INF INF units box
group lower region lower
group upper region upper
group boundary union lower upper
group mobile subtract all boundary

velocity        upper set 0 0 0 units box
velocity        lower set 0 0 0 units box

fix             1 all nvt temp 5 5 0.05

fix             4 upper addforce -1 0.0 0.0
fix             5 upper setforce NULL 0.0 0.0
fix             2 lower addforce 1 0.0 0.0
fix             3 lower setforce NULL 0.0 0.0

For second one it is:
region left block INF -160 INF 8 12 INF units box
region right block 160 INF 8 12 INF INF units box
region up block 25 50 65 INF INF INF units box
region bottom block -35 -10 INF -70 INF INF units box

group left region left
group right region right
group up region up
group bottom region bottom
group boundary union left right up bottom
group mobile subtract all boundary

velocity        boundary set 0 0 0 units box

fix 1 mobile npt temp 5 5 1 xy 100 100 1000
fix 2 boundary setforce 0.0 0.0 0.0

Could you please tell me what is the proper way to move dislocation to interact with stacking fault ? 

I sincerely apologize for taking your time.

Regards