yade-users team mailing list archive

[Nicolas Godet <question670511@xxxxxxxxxxxxxxxxxxxxx>]

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.