timplement contact age and healing rate - granular-basin - tectonic deformation experiments with Granular.jl
(HTM) git clone git://src.adamsgaard.dk/granular-basin
(DIR) Log
(DIR) Files
(DIR) Refs
(DIR) README
---
(DIR) commit b25b069b8e26363f4cb22bdd8e44e2970cb53bdb
(DIR) parent d42404e97dc00893186174cbc88eeb09e40513b0
(HTM) Author: Anders Damsgaard <anders@adamsgaard.dk>
Date: Thu, 25 Nov 2021 13:32:02 +0100
implement contact age and healing rate
Diffstat:
M init_basin.jl | 12 +++++++++++-
1 file changed, 11 insertions(+), 1 deletion(-)
---
(DIR) diff --git a/init_basin.jl b/init_basin.jl
t@@ -204,12 +204,22 @@ for i = left_edge+(bot_r/2):bot_r*1.99:left_edge+length
color = 1)
end
-Granular.findContactsAllToAll!(carpet) # find the grain contacts
+for grain in sim.grains
+ for ic=1:length(grain.contact_age)
+ grain.contact_age[ic] = 1e16
+ end
+end
+
append!(sim.grains,carpet.grains) # add the carpet grains to the main simulation object
# since the assignment will point to the carpet object, changes made to the carpet
# object will appear in the main simulation object
+for grain in sim.grains
+ grain.strength_heal_rate = 1 # new bond stengthening
+end
+
+Granular.findContactsAllToAll!(carpet) # find the grain contacts
Granular.fitGridToGrains!(sim,sim.ocean,verbose=false) # fit the ocean to the added grains