← Back to team overview

yade-users team mailing list archive

Re: [Question #201769]: Hertz mindlin, rolling back sphere

 

Question #201769 on Yade changed:
https://answers.launchpad.net/yade/+question/201769

Christian Jakob posted a new comment:
Ok, so it behaves like a pendulum, where anchor point is the middle point of two facets. 
I modified annas script a little bit. Now one can see the pendulum effect.

#!/usr/bin/python
# -*- coding: utf-8 -*-
from yade import utils

h=0.2
y_ext = 100
x_ext = 100

O.materials.append(FrictMat(young=1.7e8,poisson=0.20588,frictionAngle=radians(17),density=2650))
s1=O.bodies.append( utils.sphere([0.0000001,y_ext/2,h+0.5],.5,color=[0,1,0],fixed=False,wire=True))

v1=(0,0,h)
v6=(0,y_ext,h)
v2=(.5,0,0)
v3=(x_ext,0,0)

v4=(x_ext,y_ext,0)
v5=(0.5,y_ext,0)
 
f1=O.bodies.append( utils.facet([v1,v2,v5], dynamic=None, fixed=True, wire=True, color=None, highlight=False, noBound=False, mask=1))
f2=O.bodies.append( utils.facet([v1,v5,v6], dynamic=None, fixed=True, wire=True, color=None, highlight=False, noBound=False,  mask=1)) 
f3=O.bodies.append( utils.facet([v2,v3,v4], dynamic=None, fixed=True, wire=True, color=None, highlight=False, noBound=False, mask=1))
f4=O.bodies.append( utils.facet([v2,v4,v5], dynamic=None, fixed=True, wire=True, color=None, highlight=False, noBound=False,  mask=1))
   
## time step
O.dt=1e-4

O.engines=[
        ForceResetter(),
        InsertionSortCollider([Bo1_Sphere_Aabb(),Bo1_Facet_Aabb()]),
        InteractionLoop(
        [Ig2_Facet_Sphere_ScGeom(),Ig2_Sphere_Sphere_ScGeom()],
        [
        Ip2_FrictMat_FrictMat_MindlinPhys(en=0.1, krot=0.24, eta=1)
        ],
        [Law2_ScGeom_MindlinPhys_Mindlin(includeMoment=True),
        ]),
        NewtonIntegrator(damping=0,gravity=(0,0,-9.81)),
]

-- 
You received this question notification because you are a member of
yade-users, which is an answer contact for Yade.