yade-users team mailing list archive

[Ricardo Lorenzoni <question647827@xxxxxxxxxxxxxxxxxxxxx>]

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

Hi, I'm working on simulations with many spheres, and I have tried many ways to create a sphere factory like the inletbox (implemented in the software woo) but in my most successful implementation, the sphere factory move. Can someone help me on this issue? I need a sphere factory that creates spheres during a specific time. how can I develop this in the following code?

import sys; 
from minieigen import *; import math;
from yade import pack, plot, utils, ymport, export;

materialParede=FrictMat(density=7800,poisson=0.3,young=1.82e11,frictionAngle=radians(22.0), label="parede")

p1=(-.2,-.2,0)
p2=(-.2,.2,0)
p3=(.2,-.2,0)
p4=(.2,.2,0)

p5=(-.2,-.2,.5)
p6=(-.2,.2,.5)
p7=(.2,-.2,.5)
p8=(.2,.2,.5)

f1=utils.facet([p1,p2,p5],material=materialParede, wire=False)
f2=utils.facet([p2,p5,p6],material=materialParede, wire=False)
f3=utils.facet([p1,p3,p7],material=materialParede, wire=False)
f4=utils.facet([p1,p5,p7],material=materialParede, wire=False)
f5=utils.facet([p3,p4,p7],material=materialParede, wire=False)
f6=utils.facet([p4,p7,p8],material=materialParede, wire=False)
f7=utils.facet([p2,p4,p8],material=materialParede, wire=False)
f8=utils.facet([p2,p6,p8],material=materialParede, wire=False)
f9=utils.facet([p1,p2,p3],material=materialParede, wire=False)
f10=utils.facet([p2,p3,p4],material=materialParede, wire=False)

O.bodies.append(f1)
O.bodies.append(f2)
O.bodies.append(f3)
O.bodies.append(f4)
O.bodies.append(f5)
O.bodies.append(f6)
O.bodies.append(f7)
O.bodies.append(f8)
O.bodies.append(f9)
O.bodies.append(f10)

O.engines=[
   ForceResetter(),
   InsertionSortCollider([Bo1_Sphere_Aabb(),Bo1_Facet_Aabb()],verletDist=.05*.29e-3),
   InteractionLoop(
      # handle sphere+sphere and facet+sphere collisions
      [Ig2_Sphere_Sphere_L3Geom(),Ig2_Facet_Sphere_L3Geom()],
      [Ip2_FrictMat_FrictMat_FrictPhys()],
      [Law2_L3Geom_FrictPhys_ElPerfPl()]
   ),
   NewtonIntegrator(damping=0.04,gravity=(0,0,-9.81)),
]

O.dt=.5*utils.PWaveTimeStep()

O.run()
O.saveTmp()

Thanks for helping!

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