yade-users team mailing list archive

[MikeZhao <question275840@xxxxxxxxxxxxxxxxxxxxx>]

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

Hi there,

I've done a small test with tetra/tetraPoly element, and I found out that it's not doable to make an assembly of tetras a clump, here is the code:

b=range(2)
b[0] = utils.tetraPoly([[0,0,0],[1,0,0],[0,1,0],[0,0,1]])
b[1] = utils.tetraPoly([[0,0,0],[1,0,0],[0,1,0],[0,0,-1]])
O.bodies.appendClumped([b[0],b[1]])
#O.bodies[2] is the clump:
O.bodies[2].shape.members

##### below is the output ########
{0: (Vector3(nan,nan,nan), Quaternion((nan,nan,nan),6.283185307179586)),
 1: (Vector3(nan,nan,nan), Quaternion((nan,nan,nan),6.283185307179586))}

#and I added an egine and called O.step():
O.engines = [
       ForceResetter(),
       InsertionSortCollider([Bo1_Polyhedra_Aabb()]),
       InteractionLoop(
               [Ig2_Polyhedra_Polyhedra_PolyhedraGeom()],
               [Ip2_PolyhedraMat_PolyhedraMat_PolyhedraPhys()],
               [PolyhedraVolumetricLaw()]
       ),
       NewtonIntegrator(damping=0),
  ]
O.step()
O.bodies[0].state.pos

##### below is the output ########
Vector3(nan,nan,nan)

I want to know if this could be modified? Because dividing those tetras into different clumps can make the calculation much more efficient.

Thanks for your precious time. I'm looking forward to your reply:-)







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