yade-users team mailing list archive

[Luis Barbosa <question295587@xxxxxxxxxxxxxxxxxxxxx>]

New question #295587 on Yade:
https://answers.launchpad.net/yade/+question/295587

Hello everyone,

Consider a box (made using facetbox, facet, .stl...).
The box is created between 2 fixed spheres (one in each side) outside the box.
The box is free to move horizontally (limited by the spheres).
A step of force is given to the box, which starts to oscillate between the 2 spheres.

That is the system I'm trying to analyse, but when I use facet (fixed=False) the box simply disappears.
Does anybody have some idea?

Follows the code:

#!/usr/bin/python
# -*- coding: utf-8 -*-

# Spheres
O.materials.append(CohFrictMat(young=15e8,poisson=0.3,density=1000,frictionAngle=radians(30),label='spheres'))
O.bodies.append([
utils.sphere(center=(0.5,0.5,0.5),radius=0.1,material='spheres',fixed=False),
utils.sphere(center=(0.5,0.2,0.2),radius=0.1,material='spheres',fixed=False),
utils.sphere(center=(0.5,-0.5,0.5),radius=0.3,material='spheres',fixed=True),
utils.sphere(center=(0.5,1.5,0.5),radius=0.3,material='spheres',fixed=True)
])

########################################################################################
O.bodies.append([
	facet([[0,0,0],[1,0,0],[0,1,0]],fixed=False,color=[1,0,0]),
	facet([[1,0,0],[1,1,0],[0,1,0]],fixed=False),
	facet([[0,0,0],[0,0,1],[0,1,0]],fixed=False),
	facet([[0,0,1],[0,1,0],[0,1,1]],fixed=False),
	facet([[0,0,1],[1,0,1],[0,1,1]],fixed=False),
	facet([[1,1,1],[1,0,1],[0,1,1]],fixed=False),
	facet([[1,1,1],[1,0,1],[1,0,0]],fixed=False),
	facet([[1,1,1],[1,1,0],[1,0,0]],fixed=False),
	facet([[0,0,0],[1,0,0],[0,0,1]],fixed=False),
	facet([[0,0,1],[1,0,0],[1,0,1]],fixed=False),
	facet([[0,1,0],[1,1,0],[0,1,1]],fixed=False),
	facet([[0,1,1],[1,1,0],[1,1,1]],fixed=False)
])

O.engines=[
  ForceResetter(),
  InsertionSortCollider([Bo1_Sphere_Aabb(),Bo1_Facet_Aabb()]),
  InteractionLoop(
    [Ig2_Sphere_Sphere_ScGeom(),Ig2_Facet_Sphere_ScGeom()],
    [Ip2_FrictMat_FrictMat_FrictPhys()],
    [Law2_ScGeom_FrictPhys_CundallStrack()]
  ),

  NewtonIntegrator(damping=0.1,gravity=[0,0,-9.81]),

]
O.dt=1e-3*utils.PWaveTimeStep()

from yade import qt
qt.View()
qt.Controller()

Thanks,
Luis



-- 
You received this question notification because your team yade-users is
an answer contact for Yade.