聚乙烯PE剪切过程

#聚乙烯的剪切模拟

# 模型初始化

units           real

boundary        p p p

atom_style  molecular

neighbor        0.4 bin

neigh_modify    every 10 one 10000  

#读取数据文件

read_data       PE.data  

# 势函数参数

bond_style      harmonic

bond_coeff      1 350 1.53

angle_style     harmonic

angle_coeff     1 60 109.5

dihedral_style  multi/harmonic

dihedral_coeff  1 1.73 -4.49 0.776 6.99 0.0

pair_style      lj/cut 10.5

pair_coeff      1 1 0.112 4.01 10.5

# 温度初始化

velocity        all create 500.0 1231

# npt下驰豫

fix             1 all npt temp 500.0 500.0 50 iso 0 0 1000 drag 2

fix             2 all momentum 1 linear 1 1 1

dump       1 all custom 5000 npt1.xyz id type x y z

thermo_style    custom step temp press

thermo          100

timestep        0.5

reset_timestep  0

run             50000

unfix           1

unfix           2

undump      1

# NPT下降温

fix             1 all npt temp 500 200 50 iso 0 0 1000 drag 2

fix             2 all momentum 1 linear 1 1 1

dump       1 all custom 5000 npt2.xyz id type x y z

thermo_style    custom step temp press

thermo          10

timestep        0.5

reset_timestep  0

run             50000

unfix           1

unfix           2

undump      1

# NPT下驰豫

fix             1 all npt temp 200 200 50 iso 0 0 1000 drag 2

fix             2 all momentum 1 linear 1 1 1

dump       1 all custom 5000 npt3.xyz id type x y z

thermo_style    custom step temp press epair ebond eangle edihed pxx pyy pzz lx ly lz

thermo          10

timestep        0.5

reset_timestep  0

run             50000

unfix           1

unfix           2

undump      1

# 剪切

change_box    triclinic

fix        1 all nvt/sllod temp 200 200 1

fix             2 all deform 1 xy erate 0.02 units box remap v

dump          1 all custom 1000 shear.xyz id type x y z

timestep        0.01

reset_timestep  0

run            50000