tSet contact age to zero if there is Coulomb slip - Granular.jl - Julia package for granular dynamics simulation
HTML git clone git://src.adamsgaard.dk/Granular.jl
DIR Log
DIR Files
DIR Refs
DIR README
DIR LICENSE
---
DIR commit b640f26eebec4330354379ec734a996089ec7def
DIR parent d67f4c7667ccc392ddf6943d38f976613f9f3fff
HTML Author: Anders Damsgaard <andersd@riseup.net>
Date: Wed, 21 Feb 2018 09:07:36 -0500
Set contact age to zero if there is Coulomb slip
Diffstat:
M src/interaction.jl | 2 ++
1 file changed, 2 insertions(+), 0 deletions(-)
---
DIR diff --git a/src/interaction.jl b/src/interaction.jl
t@@ -258,6 +258,7 @@ function interactGrains!(simulation::Simulation, i::Int, j::Int, ic::Int)
force_t = abs(γ_t * vel_t)
if force_t > μ_d_minimum*abs(force_n)
force_t = μ_d_minimum*abs(force_n)
+ simulation.grains[i].contact_age[ic] = -simulation.time_step
end
if vel_t > 0.
force_t = -force_t
t@@ -270,6 +271,7 @@ function interactGrains!(simulation::Simulation, i::Int, j::Int, ic::Int)
if abs(force_t) > μ_d_minimum*abs(force_n)
force_t = μ_d_minimum*abs(force_n)*force_t/abs(force_t)
δ_t = (-force_t - γ_t*vel_t)/k_t
+ simulation.grains[i].contact_age[ic] = -simulation.time_step
end
else