yade-users team mailing list archive
-
yade-users team
-
Mailing list archive
-
Message #10756
Re: [Question #260553]: wall doesn't block spheres
Question #260553 on Yade changed:
https://answers.launchpad.net/yade/+question/260553
Description changed to:
Hi,
I am trying to simulate some spheres moving in a box.
The box is constrained in Z axis, and periodic in X,Y axis.
I use "O.periodic=True" " O.cell.setBox(0.8,0.2,0.2)" to creat periodicity and walls to make constraint.
This method is provided in
https://answers.launchpad.net/yade/+question/242137
However, I found the wall doesn't block the spheres at all, they just go through. Just like no walls exist.
And the spheres themselves don't collide when they meet, they merge and then go apart...just like every thing doesn't influence others.
I guess it may have something with the material's arguments such as young or poisson, but still don't know exactly.
Can anyone give me some advice?
Thanks so much!
_________________________________________________
from yade import *
#material
youngsModulus=1e8
poissonsRatio=0.4
coefficientFriction=0.45
matS = O.materials.append(FrictMat(young=youngsModulus,poisson=poissonsRatio,frictionAngle=coefficientFriction,label='sphereMaterial'))
# create 3 spheres
b1 = utils.sphere((0.1,0.1,0.1),.05,material=matS)
b2 = utils.sphere((0.3,0.1,0.1),.05,material=matS)
b3 = utils.sphere((0.6,0.1,0.1),.05,material=matS)
# and assign initial velocity to them
b1.state.vel = (2,0,0)
b2.state.vel = (2,0,5)
b3.state.vel = (5,0,0)
yade.O.bodies.append((b1,b2,b3))
O.periodic=True
O.cell.setBox(0.8,0.2,0.2)
O.bodies.append([
wall(0,axis=2,sense=-1),
wall(0.2,axis=2,sense=1),
])
# standard YADE engines
O.engines = [
ForceResetter(),
InsertionSortCollider([Bo1_Sphere_Aabb(),Bo1_Wall_Aabb()]),
InteractionLoop(
[Ig2_Sphere_Sphere_ScGeom(),Ig2_Wall_Sphere_ScGeom()],
[Ip2_FrictMat_FrictMat_FrictPhys()],
[Law2_ScGeom_FrictPhys_CundallStrack()]),
NewtonIntegrator(gravity=(0,0,-.3)),
]
from yade import qt; qt.Controller(); qt.View()
O.run(9,True)
--
You received this question notification because you are a member of
yade-users, which is an answer contact for Yade.