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.
_______________________________________________
Mailing list: https://launchpad.net/~yade-users
Post to : yade-users@lists.launchpad.net
Unsubscribe : https://launchpad.net/~yade-users
More help : https://help.launchpad.net/ListHelp
