Strange behaviour of Ek varying restitution coefficients
Hi everybody,
I have to simulate the behaviour of some spheres in a moving cylinder varying the tangential and normal restitution coefficients of the material. What seems strange to me is that I get a kinetic Energy increasing with lower restitution coefficients.
From the definition of coefficient of restitution I would expect the contrary.
I am wondering if I'm making a silly conceptual mistake or if there is an error in my code.
Do you see any error in the following code?
I work with yadedaily, xenial version.
Many thanks
from yade import pack,utils,
## Physical Parameters
density=7950
cyl_r = 12.7e-3
cyl_h = 24.6e-3
sp_p_r = 1.5e-3
young = 200e9
poisson = 0.305
## Parameters to be varied
#en = 0.1
#es = 0.1
en = 0.78
es = 0.78
## Creation Cylinder
facetMat=
Cylinder=
## Creation Sphere
sphereMat=
s=utils.
O.bodies.append(s)
sp=pack.
sp.makeCloud(
sp.toSimulation
## Engine
O.engines=[
ForceResetter(),
InsertionSor
InteractionLoop(
[Ig2_
[
[
),
GlobalStiffn
NewtonIntegr
PyRunner(
]
def checkUnbalanced():
if O.iter<100000: return
checker.
def nothing():
return
def addPlotData():
Ek=kineticEnergy()
plot.addData(
O.trackEnergy=True
O.saveTmp()
plot.plots=
plot.plot()
qt.Controller()
qt.View()
Question information
- Language:
- English Edit question
- Status:
- Solved
- For:
- Yade Edit question
- Assignee:
- No assignee Edit question
- Solved by:
- Chareyre
- Solved:
- Last query:
- Last reply: