Commit 7fc1e39d authored by Michel Henry's avatar Michel Henry
Browse files

little update

parent 001a378e
Pipeline #8763 passed with stages
in 5 minutes and 30 seconds
......@@ -101,7 +101,7 @@ if use_lmgc90 :
p = lmgc90Interface.ParticleProblem(2)
else :
p = scontact.ParticleProblem(2,True,True)
# p.read_vtk(outputdir,0)
p.read_vtk(outputdir,0)
p.set_friction_coefficient(0.1,"Sand","Sand") # Particle-Particle
p.set_friction_coefficient(0.1,"Sand","Steel")# Particle-Wall
# Initial time and iteration
......
......@@ -35,7 +35,7 @@ import time
import shutil
import random
def genInitialPosition(p,filename, r, H, ly, lx, rhop) :
def genInitialPosition(filename, r, H, ly, lx, rhop) :
"""Set all the particles centre positions and create the particles objects to add in the computing structure
Keyword arguments:
......@@ -47,7 +47,7 @@ def genInitialPosition(p,filename, r, H, ly, lx, rhop) :
rhop -- particles density
"""
# Particles structure builder
# p = scontact.ParticleProblem(2)
p = scontact.ParticleProblem(2)
# Load mesh.msh file specifying physical boundaries names
p.load_msh_boundaries("mesh.msh", ["Top", "Lateral","Bottom"])
......@@ -88,9 +88,9 @@ H = 0.6 # domain height
# PARTICLE PROBLEM
#
# Initialise particles
genInitialPosition(outputdir, r, H, ly, lx, rhop)
p = scontact.ParticleProblem(2)
genInitialPosition(p,outputdir, r, H, ly, lx, rhop)
# p.read_vtk(outputdir,0)
p.read_vtk(outputdir,0)
print("r = %g, m = %g\n" % (p.r()[0], p.mass()[0]))
print("RHOP = %g" % rhop)
......
Markdown is supported
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment