From 8f6991985a683488ea31f508613b7f981f0a3901 Mon Sep 17 00:00:00 2001 From: sjplimp Date: Thu, 8 Jan 2009 23:18:50 +0000 Subject: [PATCH] git-svn-id: svn://svn.icms.temple.edu/lammps-ro/trunk@2355 f3b2605a-c512-4ea7-a41b-209d697bcdaa --- README | 4 ++-- bench/in.chute | 2 +- examples/pour/in.pour | 11 +++++++++-- examples/pour/in.pour.2d | 15 ++++++++++----- 4 files changed, 22 insertions(+), 10 deletions(-) diff --git a/README b/README index e0a62631b2..5bb2d6cad2 100644 --- a/README +++ b/README @@ -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 diff --git a/bench/in.chute b/bench/in.chute index 85d0f8f50b..b2eeaa7464 100644 --- a/bench/in.chute +++ b/bench/in.chute @@ -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 diff --git a/examples/pour/in.pour b/examples/pour/in.pour index 04b956eaf3..316c537bdc 100644 --- a/examples/pour/in.pour +++ b/examples/pour/in.pour @@ -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 diff --git a/examples/pour/in.pour.2d b/examples/pour/in.pour.2d index cf926eac67..1e17db64a2 100644 --- a/examples/pour/in.pour.2d +++ b/examples/pour/in.pour.2d @@ -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 &