← Back to team overview

yade-users team mailing list archive

Re: [Question #295587]: Box Free Oscillation

 

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

    Status: Open => Answered

Jan Stránský proposed the following answer:
Hi Luis,

what about this script (slightly modified your original version)?

###################################
#!/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)
])

########################################################################################
cid,ids = O.bodies.appendClumped(aabbWalls(((0,0,0),(1,1,1)),.02,1))
box = O.bodies[cid]
box.state.vel = (0,10,0)

O.engines=[
  ForceResetter(),
  InsertionSortCollider([Bo1_Sphere_Aabb(),Bo1_Box_Aabb()]),
  InteractionLoop(
    [Ig2_Sphere_Sphere_ScGeom(),Ig2_Box_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()
###################################

cheers
Jan


2016-06-25 0:42 GMT+02:00 Luis Barbosa <question295587@xxxxxxxxxxxxxxxxxxxxx
>:

> Question #295587 on Yade changed:
> https://answers.launchpad.net/yade/+question/295587
>
>     Status: Answered => Open
>
> Luis Barbosa is still having a problem:
> Thanks everyone for the help.
>
> When I use facets or 6 thin boxes (aabbWalls) I cannot use fixed=False,
> to analyse the free movement of the box.
>
> That is the issue I've faced before.
>
> Thanks,
> Luis
>
> --
> 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     : yade-users@xxxxxxxxxxxxxxxxxxx
> Unsubscribe : https://launchpad.net/~yade-users
> More help   : https://help.launchpad.net/ListHelp
>

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