← Back to team overview

yade-users team mailing list archive

Re: [Question #186026]: I want to make a cluster of shperes representing a rock

 

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

joao miguel manso posted a new comment:
At this point, my code is like this:

import math
radius=1

idConcrete=O.materials.append(CFpmMat(poisson=0.25,young=1e9,density=4800))

O.bodies.append([
	utils.sphere((radius,radius,radius),radius,material=idConcrete),
	utils.sphere((3*radius,radius,radius),radius,material=idConcrete),
        utils.sphere((2*radius,radius+sqrt(3)*radius,radius),radius,material=idConcrete),
	utils.sphere((2*radius,radius+sqrt(3)*radius/3,radius+sqrt(6)*radius*2/3),radius,material=idConcrete),
])

O.bodies.append(utils.geom.facetBox((1,1,1),(5,5,5),wallMask=31,material=idConcrete))

O.engines=[
   ForceResetter(),
   InsertionSortCollider([Bo1_Sphere_Aabb(),Bo1_Facet_Aabb()]),
   InteractionLoop(
	[Ig2_Sphere_Sphere_ScGeom(),Ig2_Facet_Sphere_ScGeom()],
	[Ip2_CFpmMat_CFpmMat_CFpmPhys()],
	[Law2_ScGeom_CFpmPhys_CohesiveFrictionalPM(label='contactLaw')],
   ),
   GravityEngine(gravity=(0,0,-9.81)),
   NewtonIntegrator(damping=0.4),
]

contactLaw.frictionAngle=.6
contactLaw.FnMax=1e9
contactLaw.FsMax=1e9
contactLaw.isCohesive=True

O.dt=.5e-4*utils.PWaveTimeStep()
O.saveTmp()
#O.step()


I looks fine, but when I check the interactions, in the interactions inspector, none of the contactLaw properties are defined. Does anyone knows why?

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