git-svn-id: svn://svn.icms.temple.edu/lammps-ro/trunk@2355 f3b2605a-c512-4ea7-a41b-209d697bcdaa

This commit is contained in:
sjplimp 2009-01-08 23:18:50 +00:00
parent 791d231f22
commit 8f6991985a
4 changed files with 22 additions and 10 deletions

4
README
View File

@ -1,4 +1,4 @@
This is the LAMMPS (21 May 2008) software package.
This is the LAMMPS (6 Jan 2009) software package.
LAMMPS stands for Large-scale Atomic/Molecular Massively Parallel
Simulator.
@ -20,7 +20,7 @@ The primary author of the code is Steve Plimpton, who can be emailed
at sjplimp@sandia.gov. The LAMMPS WWW Site at lammps.sandia.gov has
more information about the code and its uses.
The LAMMPS (21 May 2008) distribution includes the following files
The LAMMPS (6 Jan 2009) distribution includes the following files
and directories:
README this file

View File

@ -8,7 +8,7 @@ newton off
read_data data.chute
pair_style gran/history 200000.0 50.0 0.5 0
pair_style gran/hooke/history 200000.0 NULL 50.0 NULL 0.5 0
pair_coeff * *
neighbor 0.1 bin

View File

@ -15,14 +15,21 @@ neigh_modify delay 0
# More appropriate values for realistic simulations would be
# k = 2.0e5 and dt = 0.0001, as in bench/in.chute.
pair_style gran/history 2000.0 50.0 0.5 0
pair_style gran/hooke/history 2000.0 NULL 50.0 NULL 0.5 0
#pair_style gran/hertz/history 4000.0 NULL 100.0 NULL 0.5 0
pair_coeff * *
timestep 0.001
variable kn equal 2000.0
variable gamman equal 50.0
#variable kn equal 2000.0*sqrt(2.0)
#variable gamman equal 50.0*sqrt(2.0)
fix 1 all nve/sphere
fix 2 all gravity 1.0 spherical 0.0 -180.0
fix zlower all wall/gran zplane 0.0 2000.0 50.0 0.5
fix zlower all wall/gran ${kn} NULL ${gamman} NULL 0.5 0 &
zplane 0.0 2000.0
region slab block -9.0 9.0 -9.0 9.0 10.0 15.0 units box
fix ins all pour 3000 1 300719 vol 0.13 50 region slab

View File

@ -11,20 +11,25 @@ create_box 1 reg
neighbor 0.2 bin
neigh_modify delay 0
# IMPORTANT NOTE: these values of stiffness (2000) and timestep (0.001)
# IMPORTANT NOTE: these values of stiffness (4000) and timestep (0.001)
# are used in this example file to produce a quick simulation and movie.
# More appropriate values for realistic simulations would be
# k = 2.0e5 and dt = 0.0001, as in bench/in.chute.
# k = 4.0e5 and dt = 0.0001, as in bench/in.chute (it if were Hertzian).
pair_style gran/hertzian 2000.0 50.0 0.5 0
pair_style gran/hertz/history 4000.0 NULL 100.0 NULL 0.5 0
pair_coeff * *
timestep 0.001
variable kn equal 2000.0*sqrt(2.0)
variable gamman equal 50.0*sqrt(2.0)
fix 1 all nve/sphere
fix 2 all gravity 1.0 spherical 0.0 -180.0
fix xwalls all wall/gran xplane 0 100 50 0
fix ywalls all wall/gran yplane 0 NULL 50 0
fix xwalls all wall/gran ${kn} NULL ${gamman} NULL 0 1 &
xplane 0 100
fix ywalls all wall/gran ${kn} NULL ${gamman} NULL 0 1 &
yplane 0 NULL
region slab block 1.0 99.0 30 34.5 -0.5 0.5 units box
fix ins all pour 1000 1 4767548 vol 0.4 10 &