yade-users team mailing list archive

[Pawel <question677876@xxxxxxxxxxxxxxxxxxxxx>]

New question #677876 on Yade:
https://answers.launchpad.net/yade/+question/677876

Hi everyone, 
I have very  simple simulation. At the beginning there are two spheres, with material:

O.materials.append(FrictMat(young=young_k, poisson=poisson_k, density=density_k, frictionAngle=atan(0.3),))

and positions:
O.bodies.append([
		utils.sphere(center=(0,0,0),radius=5,fixed=True,material=kula[0])
		],)
O.bodies.append([
		utils.sphere(center=(0.2*5,0,2*5),radius=promien_k,fixed=False,material=kula[0])
		],)

In loop:

InteractionLoop(
		[Ig2_Sphere_Sphere_ScGeom6D(interactionDetectionFactor=1.3)
],
		[Ip2_FrictMat_FrictMat_FrictPhys(),
		Ip2_CohFrictMat_CohFrictMat_CohFrictPhys(setCohesionNow=True,setCohesionOnNewContacts=True)
		],
		[Law2_ScGeom_FrictPhys_CundallStrack(),
		]
	),
   # apply gravity force to particles
   GravityEngine(gravity=(0,0,-9.81)),
   # damping: numerical dissipation of energy
   NewtonIntegrator(damping=0.2)
]

The upper sphere is rotating over the first one - so the simulation is working.

I'm trying to write short program in C with detection of interaction. Two spheres are in interaction but the interaction all of the time is not real - i get 0 in I->isReal()

I can't figure out how to get this to work 
FrictPhys* physFP = dynamic_cast<FrictPhys*>(I->phys.get());
tmp = physFP->shearForce;

tmp - is vector 

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