yade-dev team mailing list archive

["Sergei D." <dorofeenko@xxxxxxxxx>]

Hi (Anton, Chiara?)
Who is working on dynamic problems?
Could not you explain why it occurs (see attached script)?
And  do test the configuration with own laws?

--
Best regards,

Sergei Dorofeenko,
Cand. Phys.-Mat. Sci.

Institute of Problems of Chemical Physics,
Chernogolovka, Moscow region, Russia.

e-mail: dorofeenko@xxxxxxxxx


# encoding: utf-8
from yade import utils,geom

Rsph=0.003 

density=2400
frictionAngle=radians(35)
tc=0.0001
en=0.3
es=1.0
params=utils.getViscoelasticFromSpheresInteraction(tc,en,es)
mat=O.materials.append(ViscElMat(density=density,frictionAngle=frictionAngle,**params))

O.bodies.append(utils.facet( [ (-3*Rsph,-0.05,-0.05), (-3*Rsph,-0.05,0.05), (-3*Rsph,0.05,0) ], material=mat))

O.bodies.append([
	utils.sphere( (-0.1*Rsph,Rsph,0),Rsph,material=mat,fixed=True),
	utils.sphere( (-0.0060010254015300003,0.0042204209332099996,0.0009),Rsph,material=mat)
])

O.engines=[
	ForceResetter(),
	InsertionSortCollider([Bo1_Sphere_Aabb(),Bo1_Facet_Aabb()]),
	InteractionLoop(
		[Ig2_Sphere_Sphere_ScGeom(), Ig2_Facet_Sphere_ScGeom()],
		[Ip2_ViscElMat_ViscElMat_ViscElPhys()],
		[Law2_ScGeom_ViscElPhys_Basic()],
	),
	GravityEngine(gravity=[0,-9.81,0]),
	NewtonIntegrator(damping=0),
]

O.dt=.01*tc

O.saveTmp()

from yade import qt
qt.View()