Dear all,
I’m trying to module crack propagation with PDLAMMPS.When I was running a sample test,the error 'Non-numeric atom coords - simulation unstable (…/domain.cpp:521)'appeared.What should I do to fix it?
Thanks in advance.
My in file as follow:
units si
dimension 2
boundary s s p
atom_style peri
atom_modify map array
neighbor 0.000010 bin
lattice sq 0.00005
region box block -0.004425 0.005025 -0.000025 0.003025 -0.000025 0.000025 units box
create_box 3 box
create_atoms 3 box
region left prism 0.000025 0.001025 0.002525 0.003025 INF INF -0.0000 0 0 units box
region right prism 0.001025 0.002025 0.002525 0.003025 INF INF 0.0000 0 0 units box
region 1 prism 0.001025 0.002025 0.002525 0.003025 INF INF -0.0000 0 0 units box
delete_atoms region 1
delete_atoms region left
create_atoms 1 region left
create_atoms 2 region right
region lgt block -0.004425 -0.004325 INF INF INF INF units box
region rgt block 0.004925 0.005025 INF INF INF INF units box
pair_style peri/eps
pair_coeff * * 1.4855e11 8e10 0.000150001 0.005 0 5e8
pair_coeff 1 2 1.4855e11 8e10 0.000049999 0.005 0 5e8
delete_atoms overlap 0.000005 all all
set group all density 7860
set group all volume 1.25e-13
group lgt region lgt
group rgt region rgt
fix 1 all nve
fix 2b lgt addforce 0 -1 0
fix 2c rgt setforce 0 0 NULL
fix end1 all enforce2d
compute 1 all damage/atom
timestep 1.0e-9
thermo 2000
thermo_style custom step temp press vol
dump 1 all custom 500 dump.peri id type x y z c_1
run 100000