yade-users team mailing list archive

[Christopher Stanbridge <question256364@xxxxxxxxxxxxxxxxxxxxx>]

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

Christopher Stanbridge gave more information on the question:
damping                 = 0.4         #ratio
gravity                 = -9.80665    #m/s

mat42 = O.materials.append(CohFrictMat(isCohesive=True, frictionAngle=radians(42), density=4200, poisson=0.42, young=4.2e9, alphaKr=0.242, alphaKtw=0.1642, normalCohesion=4.2e5, shearCohesion=8.4e5))
Mat42=O.materials[mat42]

O.bodies.append(utils.facet([(-42,-42,0),(42,-42,0),(0,63,0)],wire=False,color=(0.42,0.42,0.42)))
cld = pack.SpherePack()
cld.makeCloud((-1,-1,1),(1,1,5),rMean=0.25,rRelFuzz=0.3)
cld.toSimulation(material=Mat42)

O.engines=[
	ForceResetter(),
	InsertionSortCollider([Bo1_Sphere_Aabb(aabbEnlargeFactor=1.5), Bo1_Facet_Aabb()]),
	InteractionLoop(
		[Ig2_Sphere_Sphere_ScGeom6D(interactionDetectionFactor=1.5), Ig2_Facet_Sphere_ScGeom6D()],
		[Ip2_CohFrictMat_CohFrictMat_CohFrictPhys(setCohesionNow=True)],
		[Law2_ScGeom6D_CohFrictPhys_CohesionMoment(),Law2_ScGeom_FrictPhys_CundallStrack()]
	),

	NewtonIntegrator(damping=damping, gravity=[0, 0, gravity], label='integrator')
]

O.run(1)

len(O.interactions)
#returns the number of "interactions", including distant contacts; usually between 150 and 250

allIntrs = []
for b in O.bodies:
    for i in b.intrs():
        allIntrs.append(i)
print len(allIntrs)
#always seems to return 0

-- 
You received this question notification because you are a member of
yade-users, which is an answer contact for Yade.