LAMMPS (6 Jul 2017) using 1 OpenMP thread(s) per MPI task # small Peridynamic cylinder hit by projectile units si boundary s s s atom_style peri atom_modify map array neighbor 0.0010 bin # small target lattice sc 0.0005 Lattice spacing in x,y,z = 0.0005 0.0005 0.0005 region target cylinder y 0.0 0.0 0.0050 -0.0050 0.0 units box create_box 1 target Created orthogonal box = (-0.005 -0.005 -0.005) to (0.005 0 0.005) 1 by 1 by 1 MPI processor grid create_atoms 1 region target Created 3487 atoms pair_style peri/pmb pair_coeff * * 1.6863e22 0.0015001 0.0005 0.25 set group all density 2200 3487 settings made for density set group all volume 1.25e-10 3487 settings made for volume velocity all set 0.0 0.0 0.0 sum no units box fix 1 all nve # spherical indenter to shatter target variable y0 equal 0.00155 variable vy equal -100 variable y equal "v_y0 + step*dt*v_vy" fix 2 all indent 1e17 sphere 0.0000 v_y 0.0000 0.0015 units box compute 1 all damage/atom timestep 1.0e-7 thermo 100 #dump 1 all custom 100 dump.peri id type x y z c_1 #dump 2 all image 50 image.*.jpg type type # axes yes 0.8 0.02 view 80 -30 adiam 0.0006 #dump_modify 2 pad 4 #dump 3 all movie 50 movie.mpg type type # axes yes 0.8 0.02 view 80 -30 adiam 0.0006 #dump_modify 3 pad 4 run 1000 Neighbor list info ... update every 1 steps, delay 10 steps, check yes max neighbors/atom: 2000, page size: 100000 master list distance cutoff = 0.0025001 ghost atom cutoff = 0.0025001 binsize = 0.00125005, bins = 9 5 9 2 neighbor lists, perpetual/occasional/extra = 1 1 0 (1) pair peri/pmb, perpetual attributes: half, newton on pair build: half/bin/atomonly/newton stencil: half/bin/3d/newton bin: standard (2) fix PERI_NEIGH, occasional attributes: full, newton on pair build: full/bin/atomonly stencil: full/bin/3d bin: standard Peridynamic bonds: total # of bonds = 335966 bonds/atom = 96.3482 Per MPI rank memory allocation (min/avg/max) = 34.79 | 34.79 | 34.79 Mbytes Step Temp E_pair E_mol TotEng Press Volume 0 0 0 0 0 0 5.0030006e-07 100 1.7890585e+24 552721.8 0 681881.47 1.7210968e+11 5.0030006e-07 200 1.3549879e+27 3097027.6 0 1.0091931e+08 1.1829272e+14 5.5130066e-07 300 3.3009625e+27 6331254.8 0 2.4464163e+08 1.9647265e+14 8.0862953e-07 400 3.815184e+27 6225081.7 0 2.8165928e+08 1.8189267e+14 1.0095118e-06 500 4.2580877e+27 20212686 0 3.2762196e+08 1.6249923e+14 1.2611723e-06 600 5.5126512e+27 30861342 0 4.2884284e+08 1.7320038e+14 1.531873e-06 700 1.1807414e+28 23119941 0 8.7554687e+08 2.9477434e+14 1.9278632e-06 800 1.2424839e+28 2407365.1 0 8.994088e+08 2.3787786e+14 2.5138992e-06 900 1.2358397e+28 4532424.3 0 8.9673716e+08 1.9097316e+14 3.1145903e-06 1000 1.2341048e+28 3219355.8 0 8.9417154e+08 1.5968585e+14 3.7196039e-06 Loop time of 28.565 on 1 procs for 1000 steps with 3487 atoms 99.5% CPU use with 1 MPI tasks x 1 OpenMP threads MPI task timing breakdown: Section | min time | avg time | max time |%varavg| %total --------------------------------------------------------------- Pair | 27.721 | 27.721 | 27.721 | 0.0 | 97.04 Neigh | 0.66353 | 0.66353 | 0.66353 | 0.0 | 2.32 Comm | 0.0027969 | 0.0027969 | 0.0027969 | 0.0 | 0.01 Output | 0.00042295 | 0.00042295 | 0.00042295 | 0.0 | 0.00 Modify | 0.1566 | 0.1566 | 0.1566 | 0.0 | 0.55 Other | | 0.02086 | | | 0.07 Nlocal: 3487 ave 3487 max 3487 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 0 ave 0 max 0 min Histogram: 1 0 0 0 0 0 0 0 0 0 Neighs: 567132 ave 567132 max 567132 min Histogram: 1 0 0 0 0 0 0 0 0 0 FullNghs: 1.20908e+06 ave 1.20908e+06 max 1.20908e+06 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 1209076 Ave neighs/atom = 346.738 Neighbor list builds = 46 Dangerous builds = 0 Please see the log.cite file for references relevant to this simulation Total wall time: 0:00:28