Question #670511 on Yade changed:
https://answers.launchpad.net/yade/+question/670511
Status: Answered => Open
Nicolas Godet is still having a problem:
I setup the expiration date to "Never" for the pastebin so it should not
be a problem for the future but I will put the script directly here now.
> Could you show this one? You say it is easy to switch but I can' t
figure out what to comment/uncomment to see this elementary case.
from yade import plot,qt
from yade.gridpfacet import *
import math
O.dt = 0.000001
Rad = 0.1
Factor=1.5
m_viscosity=3
m_roughness=0.02
O.materials.append(
FrictMat(
density=2500,
frictionAngle=math.radians(35.),
poisson=0.2,
young=1e7,
label='Frict'))
S2 =
O.bodies.append(sphere(center=[1.5*Rad,0,Factor*2.1*Rad],radius=Rad,material='Frict',fixed=False))
Stest =
O.bodies.append(sphere(center=[1.5*Rad,0,3*Factor*2.1*Rad],radius=Rad,material='Frict',fixed=False))
W1 =
O.bodies.append(aabbWalls([(4*Rad,-2*Rad,15*Rad),(0,3*Rad,0)],thickness=0.,material='Frict',oversizeFactor=1.0))
for i in [0,1,3,4,5]:
O.bodies.erase(W1[i])
O.engines=[
ForceResetter(),
InsertionSortCollider(
[Bo1_Sphere_Aabb(aabbEnlargeFactor=Factor, label="aabb"),
Bo1_Box_Aabb(),
Bo1_Wall_Aabb(),],#Bo1_Facet_Aabb(),
verletDist=-0.1,
allowBiggerThanPeriod=False),
InteractionLoop(
[Ig2_Sphere_Sphere_ScGeom6D(interactionDetectionFactor=Factor,label="Ig2"),
Ig2_Wall_Sphere_ScGeom(),
Ig2_Box_Sphere_ScGeom6D(),],
[Ip2_FrictMat_FrictMat_LubricationPhys(eta=m_viscosity,eps=m_roughness),],
[Law2_ScGeom_ImplicitLubricationPhys(activateNormalLubrication=True,activateTangencialLubrication=True,activateTwistLubrication=False,activateRollLubrication=False,
theta=0.55,resolution=1,NewtonRafsonTol=1e-6,warnedOnce=True,maxSubSteps=20,debug=False,),]),
NewtonIntegrator(damping=0.,gravity=[0,-5,-5],label='newton'),]
qt.View()
And here is with facet instead of box :
from yade import plot,qt
from yade.gridpfacet import *
import math
O.dt = 0.000001
Rad = 0.1
Factor=1.5
m_viscosity=3
m_roughness=0.02
O.materials.append(
FrictMat(
density=2500,
frictionAngle=math.radians(35.),
poisson=0.2,
young=1e7,
label='Frict'))
S2 =
O.bodies.append(sphere(center=[1.5*Rad,0,Factor*2.1*Rad],radius=Rad,material='Frict',fixed=False))
Stest =
O.bodies.append(sphere(center=[1.5*Rad,0,3*Factor*2.1*Rad],radius=Rad,material='Frict',fixed=False))
F1=O.bodies.append(facet([[-10,-2*Rad,0],[10,-2*Rad,0],[0,-2*Rad,10]],fixed=True,material='Frict'))
O.engines=[
ForceResetter(),
InsertionSortCollider(
[Bo1_Sphere_Aabb(aabbEnlargeFactor=Factor, label="aabb"),
Bo1_Facet_Aabb(),],
verletDist=-0.1,
allowBiggerThanPeriod=False),
InteractionLoop(
[Ig2_Sphere_Sphere_ScGeom6D(interactionDetectionFactor=Factor,label="Ig2"),
Ig2_Facet_Sphere_ScGeom()],
[Ip2_FrictMat_FrictMat_LubricationPhys(eta=m_viscosity,eps=m_roughness),],
[Law2_ScGeom_ImplicitLubricationPhys(activateNormalLubrication=True,activateTangencialLubrication=True,activateTwistLubrication=False,activateRollLubrication=False,
theta=0.55,resolution=1,NewtonRafsonTol=1e-6,warnedOnce=True,maxSubSteps=20,debug=False,),]),
NewtonIntegrator(damping=0.,gravity=[0,-5,-5],label='newton'),]
qt.View()
Nicolas
--
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 : [email protected]
Unsubscribe : https://launchpad.net/~yade-users
More help : https://help.launchpad.net/ListHelp
