yade-users team mailing list archive

[Luis Barbosa <question249243@xxxxxxxxxxxxxxxxxxxxx>]

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

    Status: Answered => Open

Luis Barbosa is still having a problem:
Hi Jan, this is my script, when there is interaction the value of force is printed and when not is printed 'no interaction...'
If I'm able to export data could solve my problem for while.

tks, Luis


from yade import plot

O.materials.append(FrictMat(young=15e6,poisson=.4,frictionAngle=radians(20),density=4000,))
O.bodies.append(utils.geom.facetBox((.5,.5,.5),(.5,.5,.5),wallMask=31))
O.materials.append(CohFrictMat(young=15e6,poisson=0.3,density=2600,frictionAngle=radians(30),isCohesive=True,normalCohesion=1e10,shearCohesion=1e10,momentRotationLaw=True,etaRoll=0.1,label='spheres'))
O.bodies.append([utils.sphere(center=(0.5,0.5,0.5),radius=.07)])
O.engines=[
   ForceResetter(),
   InsertionSortCollider([Bo1_Sphere_Aabb(),Bo1_Facet_Aabb()]),
InteractionLoop(
[Ig2_Sphere_Sphere_ScGeom6D(),Ig2_Facet_Sphere_ScGeom()],
[Ip2_FrictMat_FrictMat_MindlinPhys(en=0.1,es=0.1,label='damping'), Ip2_CohFrictMat_CohFrictMat_CohFrictPhys (setCohesionNow=True,label="cohesiveIp")],
[Law2_ScGeom_MindlinPhys_Mindlin(label='contactLaw'),Law2_ScGeom6D_CohFrictPhys_CohesionMoment(useIncrementalForm=True,
always_use_moment_law=True,label='cohesiveLaw')]
),

   GravityEngine(gravity=(0,0,-9.81)),
   NewtonIntegrator(damping=0.0),
   PyRunner(command='myAddData()',iterPeriod=10),

]
#______________________________________________________________________
damping.useDamping=True

O.dt=.5e-2*utils.PWaveTimeStep()

def myAddData():
	try:
	 i=O.interactions[9,10]
	 print i.phys.normalForce
	except:
		print "No interaction between"

O.saveTmp()

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