Lack of recoil in sphere-pfacet impact?
I am studying the interaction of pfacets with other objects, since I have it not completely clear in my mind. I've seen that pfacets interact "apparently realistically" with other pfacets. However I have an example where a sphere collides with a pfacet, bounces back correctly but does not transfer any momentum to the pfacet (no nodes or connections move after the end of the impact). Why?
This is my code
-------
from yade.gridpfacet import *
from yade import pack,ymport,
vr=0.2
vh=0.4
rr=0.02
phimat=30. # friction angle
E=1e8 # Young's modulus
f = open("state.dat", "w")
intf = 0
def state(si,intf):
if (len(O.
forc = O.interactions.
else:
forc = 0
intf += forc*O.
f.write ("%d %f %f %f %f\n" % (O.iter,
# tempo coord y sfera vel y sfera forza sulla sfera integr-forza
f.flush()
return intf
#### Engines ####
O.engines=[
]),
[
],
[
],
[
]
),
# qt.SnapshotEngi
]
#### Creat materials ####
O.materials.append( CohFrictMat( young=E*
O.materials.append( FrictMat( young=E,
##### Create ONE WALL #####
xw = 0
yw = 0
lwall = 10*vr
hwall = 2*vh
v1=Vector3(xw,yw,0)
v2=Vector3(
v3=Vector3(
vertices=[v1,v2,v3]
pfacetCreator1(
#sphId = O.bodies.
sphId = O.bodies.
O.bodies[
#### For viewing ####
from yade import qt
qt.View()
Gl1_Sphere.
qtr = qt.Renderer()
#### Set a time step ####
#O.dt=1e-05
O.dt=PWaveTimeS
#### Allows to reload the simulation ####
O.saveTmp()
Question information
- Language:
- English Edit question
- Status:
- Expired
- For:
- Yade Edit question
- Assignee:
- No assignee Edit question
- Last query:
- Last reply: