详细内容

LAMMPS讲解29-铁表面镀碳膜-a toy case

铁表面镀碳膜,一个简单的示例。

in文件如下

# ----------------- Init Section -----------------

 

  boundary p p p

  neighbor 3.0 bin

  neigh_modify every 1 delay 0 check yes

  atom_style atomic

  units metal

  pair_style   tersoff/zbl

 

 

# ----------------- Atom Definition Section -----------------

 lattice bcc 2.855324

 region box block 0 10 0 10 0 20

 create_box  2 box

 region Fe block 0 10 0 10 0 4

 create_atoms  1 region Fe

 

# ----------------- Run Section -----------------

 

 pair_coeff * * FeC.tersoff.zbl Fe C

 

 mass 1 55.845

 mass 2 12.011

 

 group Fe type 1

 group C  type 2

 region fix block INF INF INF INF INF 1

 group fix region fix

 group Femove subtract Fe fix

 

 timestep 0.002

 region dep block 0.5 28 0.5 28 25 56 units box

 fix dep C deposit 2000 2 100 4564456 region dep attempt 100 near 3.0 vz -1.0 -0.1 units box

 compute myT C temp

 compute_modify myT dynamic/dof yes

 variable   Cd atom "type==2"

 group           Cd dynamic all var Cd

 variable        nC equal count(Cd)

 velocity     Femove create 300 372748 rot yes dist gaussian loop local

 thermo          100

 thermo_style  custom step temp v_nC

 thermo_modify temp myT flush yes lost ignore

  

 fix blnc all balance 1000 1.05 shift xyz 10 1.05

 dump mydump all custom 1000 dump.lammpstrj id type mass x y z

 

 fix fxnvt1 C nve

 fix fxnvt2 Femove nvt temp 300 300 0.2

 

 run 200000

 write_data equ.data

 write_restart restart.equ

 

 

感谢鲍路瑶老师的分享,内容来自于鲍老师分享出来的资料

如有需要添加微信:lmp_zhushou  进入微信群,帮助他人,共建社区

获取完整版lammps讲义可以加微信lmp_zhushou或加入QQ群994359511


最新评论
请先登录才能进行回复登录
技术支持: CLOUD | 管理登录
seo seo