This is simple model to simulate falling sand grains.
However, It takes more than 3 to 4 days to reach equilibrium state.
Do you have any suggestion to increase calculation speed?
====
model new
model configure dynamic
model large-strain on
model deterministic off
;
model domain condition destroy
;
domain extent -0.025 0.025 -0.025 0.025 -0.05 0.40
;
wall generate id 1 plane position (0.0, 0.0, -0.025)
wall generate id 2 plane dip 90 position (0.0, -0.025, 0.0)
wall generate id 3 plane dip 90 position (0.0, 0.025, 0.0)
wall generate id 4 plane dip 90 dip-direction 90.0 position (0.025, 0.0, 0.0)
wall generate id 5 plane dip 90 dip-direction 90.0 position (-0.025, 0.0, 0.0)
;
fish define gen_balls
loop local i (1, 10)
command
ball generate box -0.025 0.025 -0.025 0.025 0.35 0.4 radius 0.001 0.002 gauss number 10000 tries 10000000
end_command
end_loop
end
@gen_balls
;
ball attribute density 2700
;contact cmat default model linear property kn 1e6 ks 1e3 dp_nratio 0.7 dp_sratio 0.7 fric 0.4
contact cmat default type ball-ball model hertz property hz_shear 31.14e9 hz_poiss 0.167 fric 0.6 dp_nratio 0.7 dp_sratio 0.77
contact cmat default type ball-facet model linear property kn 1e6 ks 1e3 fric 0.4 dp_nratio 0.7 dp_sratio 0.7
model gravity 9.81
;
model history id ‘dynamic-time’ dynamic time-total
model display history ‘dynamic-time’
model dynamic active on ; timestep increment 1.2
model solve