Example #1
0
steel.Kn = 10*gravel.Kn
steel.frictionAngle = 0.4 #rad

rubber = PolyhedraMat()
rubber.density = 1000 #kg/m^3 
rubber.Ks = gravel.Ks/10
rubber.Kn = gravel.Kn/10
rubber.frictionAngle = 0.7 #rad

O.bodies.append(polyhedra_utils.polyhedra(gravel,v=((0,0,-0.05),(0.3,0,-0.05),(0.3,0.3,-0.05),(0,0.3,-0.05),(0,0,0),(0.3,0,0),(0.3,0.3,0),(0,0.3,0)),fixed=True, color=(0.35,0.35,0.35)))
#O.bodies.append(utils.wall(0,axis=1,sense=1, material = gravel))
#O.bodies.append(utils.wall(0,axis=0,sense=1, material = gravel))
#O.bodies.append(utils.wall(0.3,axis=1,sense=-1, material = gravel))
#O.bodies.append(utils.wall(0.3,axis=0,sense=-1, material = gravel))

polyhedra_utils.fillBox((0,0,0), (0.3,0.3,0.3),gravel,sizemin=[0.025,0.025,0.025],sizemax=[0.05,0.05,0.05],seed=4)

def checkUnbalancedI():   
    print "iter %d, time elapsed %f,  time step %.5e, unbalanced forces = %.5f"%(O.iter, O.realtime, O.dt, utils.unbalancedForce())		

O.engines=[
   ForceResetter(),
   InsertionSortCollider([Bo1_Polyhedra_Aabb(),Bo1_Wall_Aabb(),Bo1_Facet_Aabb()]),
   InteractionLoop(
      [Ig2_Wall_Polyhedra_PolyhedraGeom(), Ig2_Polyhedra_Polyhedra_PolyhedraGeom(), Ig2_Facet_Polyhedra_PolyhedraGeom()], 
      [Ip2_PolyhedraMat_PolyhedraMat_PolyhedraPhys()], # collision "physics"
      [PolyhedraVolumetricLaw()]   # contact law -- apply forces
   ),
   #GravityEngine(gravity=(0,0,-9.81)),
   NewtonIntegrator(damping=0.3,gravity=(0,0,-9.81)),
   PyRunner(command='checkUnbalancedI()',realPeriod=5,label='checker')
Example #2
0
O.bodies.append(
    polyhedra_utils.polyhedra(gravel,
                              v=((0, 0, -0.05), (0.3, 0, -0.05),
                                 (0.3, 0.3, -0.05), (0, 0.3, -0.05), (0, 0, 0),
                                 (0.3, 0, 0), (0.3, 0.3, 0), (0, 0.3, 0)),
                              fixed=True,
                              color=(0.35, 0.35, 0.35)))
#O.bodies.append(utils.wall(0,axis=1,sense=1, material = gravel))
#O.bodies.append(utils.wall(0,axis=0,sense=1, material = gravel))
#O.bodies.append(utils.wall(0.3,axis=1,sense=-1, material = gravel))
#O.bodies.append(utils.wall(0.3,axis=0,sense=-1, material = gravel))

polyhedra_utils.fillBox((0, 0, 0), (0.3, 0.3, 0.3),
                        gravel,
                        sizemin=[0.025, 0.025, 0.025],
                        sizemax=[0.05, 0.05, 0.05],
                        seed=4)


def checkUnbalancedI():
    print "iter %d, time elapsed %f,  time step %.5e, unbalanced forces = %.5f" % (
        O.iter, O.realtime, O.dt, utils.unbalancedForce())


O.engines = [
    ForceResetter(),
    InsertionSortCollider(
        [Bo1_Polyhedra_Aabb(),
         Bo1_Wall_Aabb(),
         Bo1_Facet_Aabb()]),
Example #3
0
from yade import polyhedra_utils,export
polyhedra_utils.fillBox((0,0,0), (0.3,0.3,0.3),defaultMaterial(),sizemin=(0.025,0.025,0.025),sizemax=(0.05,0.05,0.05),seed=4)
export.textPolyhedra('/tmp/textPolyhedra.txt')
Example #4
0
O.bodies.append(
    utils.facet(((0.0, 0.0, 0.6), (0.57, 0.0, 0.6), (0.57, 0.57, 0.6)),
                dynamic=None,
                fixed=True,
                material=matF))
O.bodies.append(
    utils.facet(((0.0, 0.0, 0.6), (0.0, 0.57, 0.6), (0.57, 0.57, 0.6)),
                dynamic=None,
                fixed=True,
                material=matF))

## add polyhedron
ballast = polyhedra_utils.fillBox((0.17, 0.17, 0.6), (0.40, 0.40, 0.8),
                                  matP,
                                  sizemin=[0.026, 0.026, 0.026],
                                  sizemax=[0.038, 0.038, 0.038],
                                  ratio=[1, 1, 1],
                                  seed=4,
                                  mask=1)

#O.bodies.append(wall((0,0,0),2))

O.engines = [
    ForceResetter(),
    InsertionSortCollider(
        [Bo1_Polyhedra_Aabb(),
         Bo1_Wall_Aabb(),
         Bo1_Facet_Aabb()]),
    InteractionLoop(
        [
            Ig2_Wall_Polyhedra_PolyhedraGeom(),
Example #5
0
from yade import polyhedra_utils, export
polyhedra_utils.fillBox((0, 0, 0), (0.3, 0.3, 0.3),
                        defaultMaterial(),
                        sizemin=(0.025, 0.025, 0.025),
                        sizemax=(0.05, 0.05, 0.05),
                        seed=4)
export.textPolyhedra('/tmp/textPolyhedra.txt')