Skip to content
Snippets Groups Projects
Commit be48dae6 authored by Rafael Ravedutti's avatar Rafael Ravedutti
Browse files

Adjust gravity calculation 


Signed-off-by: default avatarRafael Ravedutti <rafaelravedutti@gmail.com>
parent 145155f7
Branches
Tags
1 merge request!1Implement DEM and many other features
Hide whitespace changes
Inline Side-by-side
Showing
with 2 additions and 2 deletions
examples/dem.py
+ 2
− 2
View file @ be48dae6
...@@ -89,7 +89,7 @@ def euler(i): ...@@ -89,7 +89,7 @@ def euler(i):
def gravity(i): def gravity(i):
volume = (4.0 / 3.0) * pi * radius[i] * radius[i] * radius[i] volume = (4.0 / 3.0) * pi * radius[i] * radius[i] * radius[i]
force[i][2] += -(densityParticle_SI - densityFluid_SI) * volume * gravity_SI force[i][2] = force[i][2] - (densityParticle_SI - densityFluid_SI) * volume * gravity_SI
cmd = sys.argv[0] cmd = sys.argv[0]
...@@ -112,7 +112,7 @@ restitutionCoefficient = 0.1 ...@@ -112,7 +112,7 @@ restitutionCoefficient = 0.1
collisionTime_SI = 5e-4 collisionTime_SI = 5e-4
poissonsRatio = 0.22 poissonsRatio = 0.22
timeSteps = 10000 timeSteps = 10000
visSpacing = 1 visSpacing = 100
denseBottomLayer = False denseBottomLayer = False
bottomLayerOffsetFactor = 1.0 bottomLayerOffsetFactor = 1.0
kappa = 2.0 * (1.0 - poissonsRatio) / (2.0 - poissonsRatio) # from Thornton et al kappa = 2.0 * (1.0 - poissonsRatio) / (2.0 - poissonsRatio) # from Thornton et al
......
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Please register or to comment