# 1d bar units si dimension 2 boundary s s p atom_style peri atom_modify map array neighbor 0.003 bin lattice sq 0.0005 region box block -0.0048 0.0048 -0.0006 0.0006 -0.00052 0.00052 units box region plate block -0.00475 0.00475 -0.0005 0.0005 -0.0005 0.0005 units box create_box 4 box create_atoms 1 region plate pair_style peri/lps pair_coeff * * 68e9 24e9 0.001 0.00335 0.2 pair_coeff 3 * 68e9 24e9 0.001 1 0.2 pair_coeff 4 * 68e9 24e9 0.001 1 0.2 set group all density 2700 set group all volume 125e-12 velocity all set 0.0 0.0 0.0 sum no units box region 3 block -0.00475 -0.00445 INF INF INF INF units box group left region 3 region 4 block 0.00445 0.00475 INF INF INF INF units box group right region 4 region 2 block -0.00445 0.00445 INF INF INF INF units box group center region 2 set group center type 2 set group right type 4 set group left type 3 fix 1 all nve fix 3 left addforce -4e8 0 0 fix 4 right addforce 4e8 0 0 fix 5 all enforce2d compute 1 all displace/atom timestep 5e-8 thermo 20 dump MyDump all custom 1 dump.data id x y z run 50000