Commit 9c9e2a7d authored by Nathan Coppin's avatar Nathan Coppin
Browse files

test case with problem for contact management

parent d4f0e665
Pipeline #4809 passed with stage
in 24 seconds
......@@ -59,12 +59,12 @@ ii = 0
#physical parameters
g = -9.81 # gravity
rhop = 2450 # grains density
tEnd = 50 # final time
tEnd = 0.8 # final time
#numerical parameters
dt = 1e-3 # time step
#geometry parameters
rout = 1000 # outer radius
rout = 0.2 # outer radius
r = 5e-3 # grains radius
shutil.copy("mesh.msh", outputdir +"/mesh.msh")
......@@ -74,13 +74,14 @@ p = scontact.ParticleProblem(2)
p.read_vtk(outputdir,0)
outf = 1 # number of iterations between output files
#Taking friction into account
p.add_friction([0.3090, 1.0],[0.1545, 0.5])
p.add_friction([0.3090, 1.0])
p.set_friction_relaxation(0.95)
#Computation loop
forces = np.zeros_like(p.velocity())
k = 0
while t < tEnd :
forces[:,1] = g*p.mass()[0]*cos(radians(25))
forces[:,0] = -g*p.mass()[0]*sin(radians(25))
forces[:,1] = g*p.mass()[0]*cos(radians(35))
forces[:,0] = -g*p.mass()[0]*sin(radians(35))
#Computation of the new velocities
vn = p.velocity() + forces * dt / p.mass()
vmax = np.max(np.hypot(vn[:, 0], vn[:, 1]))
......
rout = 1000;
rout = 0.2;
lcout = rout/20;
Point(1) = {0, 0, 0};
......
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