yade-users team mailing list archive

[Jan Stránský <question703462@xxxxxxxxxxxxxxxxxxxxx>]

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

    Status: Needs information => Answered

Jan Stránský proposed the following answer:
###
O.materials.append(CohFrictMat(
    isCohesive=True,
    normalCohesion=1e99,
    shearCohesion=1e99,
))

O.engines = [
    ForceResetter(),
    InsertionSortCollider([Bo1_Sphere_Aabb()]),
    InteractionLoop(
        [Ig2_Sphere_Sphere_ScGeom6D()],
        [Ip2_CohFrictMat_CohFrictMat_CohFrictPhys(
            setCohesionNow=True,
            setCohesionOnNewContacts=False,
        )],
        [Law2_ScGeom6D_CohFrictPhys_CohesionMoment()]
    ),
    NewtonIntegrator(),
]

spheres = s1,s2,s3 = (
    sphere((0,0,0),1),
    sphere((1,0,0),1),
    sphere((0,1,0),1),
)
O.bodies.append(spheres)

O.dt = 0
O.step()
for i in O.interactions:
    i.phys.unp = i.geom.penetrationDepth

O.dt = 1e-7
#s1.state.vel = (-200,0,0)
###

Cheers
Jan

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