聚乙烯拉伸

#聚乙烯拉伸

# 模型初始化

units    real

boundary  p p p

atom_style  full

#设置键势类型

bond_style      harmonic

#设置键角势类型

angle_style     harmonic

#设置二面角势类型

dihedral_style  harmonic

#设置对势类型

pair_style  lj/cut 10.5

#读取聚乙烯结构数据

read_data  pp.data

#设置邻居参数

neighbor  0.4 bin

neigh_modify  every 10 one 10000

#nve系综下驰豫

dump 1 all custom 100 nve.xyz id type x y z vx vy vz

#温度初始化

velocity   all create 300.0 1231

#设置nve系综

fix    1 all nve

#热力学输出

thermo_style  custom step temp press

thermo          100

run    1000

#取消fix、dump设置

unfix 1

undump 1

#npt系综下驰豫 

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

#dump输出设置

dump 1 all custom 100 npt.xyz id type x y z vx vy vz

#热力学输出

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

thermo          10

#重新设置模拟步长

timestep  0.5

#步数初始化

reset_timestep   0

#运行50000步

run    50000

#取消设置

unfix 1

undump  1

#npt系综下拉伸

#计算应变值

variable tmp equal "lx"

variable L0 equal ${tmp}

variable strain equal "(lx - v_L0)/v_L0"

#dump输出

dump 1 all custom 1000 dump.pp id type x y z vx vy vz

#设置npt系综

fix    1 all npt temp 100 100 50 y 0 0 1000 z 0 0 1000 drag 2

#deform拉伸，x方向

fix    2 all deform 1 x erate 1e-5 units box remap x

#热力学输出

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

thermo          100

reset_timestep  0

#启动模拟，20万步

run    200000

unfix 1

unfix 2

undump 1