Question #707430 on Yade changed:
https://answers.launchpad.net/yade/+question/707430
Status: Solved => Open
Paul Schröder is still having a problem:
Thanks for your feedback, Karol and thanks for the sample code,
that helped me a lot and also answered my question so far. I just have a small
question about this, which is rather less related to the friction value itself:
Why does the rotational speed of the particle sp1 not decrease despite the
counteracting friction force. How can friction dissipation be activated in the
Mindlin-material-model?
In addition, for future readers, I would like to extend the above
example in the code shown below to include the friction values from the
different material pairings distinguished in my other question.
With best wishes, Paul
from yade import plot
## PARAMS
mu_static00 = 0.8
mu_static01 = 0.4
mu_static11 = 0.1
mu_dynamic00 = 0.5
mu_dynamic01 = 0.3
mu_dynamic11 = 0.1
# FUNCTIONS
def check_interactions():
for i in O.interactions:
b1, b2 = O.bodies[i.id1],O.bodies[i.id2]
if i.phys.isSliding:
if b1.mat.id == idMat0 and b2.mat.id == idMat0:
i.phys.tangensOfFrictionAngle = mu_dynamic00
elif (b1.mat.id == idMat0 and b2.mat.id == idMat1) or (b1.mat.id ==
idMat1 and b2.mat.id == idMat0):
i.phys.tangensOfFrictionAngle = mu_dynamic01
elif b1.mat.id == idMat1 and b2.mat.id == idMat1:
i.phys.tangensOfFrictionAngle = mu_dynamic11
else:
print("Kombination Gleitreibpartner nicht vorgesehen")
else:
if b1.mat.id == idMat0 and b2.mat.id == idMat0:
i.phys.tangensOfFrictionAngle = mu_static00
elif (b1.mat.id == idMat0 and b2.mat.id == idMat1) or (b1.mat.id ==
idMat1 and b2.mat.id == idMat0):
i.phys.tangensOfFrictionAngle = mu_static01
elif b1.mat.id == idMat1 and b2.mat.id == idMat1:
i.phys.tangensOfFrictionAngle = mu_static11
else:
print("Kombination Haftreibpartner nicht vorgesehen")
## MATERIALS
idMat0 = O.materials.append(FrictMat(frictionAngle = atan(mu_static00),
label='Mat0'))
idMat1 = O.materials.append(FrictMat(frictionAngle = atan(mu_static11),
label='Mat1'))
## BODIES
sp1 = O.bodies.append(sphere((0,0,0),0.1, fixed = True, material =
O.materials['Mat0']))
sp2 = O.bodies.append(sphere((0,0.199,0),0.1, fixed = True, material =
O.materials['Mat1']))
# rotate one sphere
O.bodies[sp1].state.angVel = (2.5e-3,0,0)
## ENGINES
O.engines=[
ForceResetter(),
InsertionSortCollider([Bo1_Sphere_Aabb()]),
InteractionLoop(
[Ig2_Sphere_Sphere_ScGeom()],
[Ip2_FrictMat_FrictMat_MindlinPhys()],
[Law2_ScGeom_MindlinPhys_Mindlin()],
),
NewtonIntegrator(),
PyRunner(command = "plot.addData(t = O.time, F =
O.interactions[0,1].phys.shearForce.norm(), v1 = O.bodies[sp1].state.angVel[0],
v2 = O.bodies[sp2].state.angVel[0])", virtPeriod = 1),
PyRunner(command = "check_interactions()", iterPeriod = 100),#
]
O.dt = PWaveTimeStep()*0.5
plot.plots = {'t':(('F','b-'),None,('v1','g--'),('v2','r--'))}
plot.plot(subPlots = False)
O.run()
--
You received this question notification because your team yade-users is
an answer contact for Yade.
___
Mailing list: https://launchpad.net/~yade-users
Post to : yade-users@lists.launchpad.net
Unsubscribe : https://launchpad.net/~yade-users
More help : https://help.launchpad.net/ListHelp
