yade-users team mailing list archive

[Dennis Defri <question694168@xxxxxxxxxxxxxxxxxxxxx>]

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

Description changed to:
Dear all,

I am very new and only learn YADE from Youtube of Mr. Janek Kozicki, so
I only know a few things. Probably my problem is in principle but
unfortunately I cannot locate the problem itself more accurately. I am
just assuming that it is in the interaction especially on the functors,
I dunno how to use it efficiently.

I am using:
Yade 2020.01a 
Using python version: 3.8.5 (default, Jul 28 2020, 12:59:40) 
[GCC 9.3.0]
TCP python prompt on localhost:9000, auth cookie `kcysud'
XMLRPC info provider on http://localhost:21000

I am trying to build up a simulation for a simple Rockfall model, that if I can describe it, there are a boulder and grain (Z Y view), I want to make the boulder drops into the grains, and I want to make the grains is scattered under the weight of the boulder. The floor (I am using wall) has no interaction, just as a base/bed.
I tried to input other functor which is [Ip2_FrictMat_FrictMat_FrictPhys()] in twice, because if I input only once the "segmentation fault" will be happened. 
If this script of mine is working properly, the sequence of simulation should be as follow:
Generate a sphere > create Spherepack > create a wall (as a floor) > a sphere (boulder) drops into the grains (spherepack)

I noticed that other researchers had asked similar question in Launchpad, but I still can't find a solution to my problem here.
My script is as follow:

##Define Geometry
O.bodies.append([utils.sphere(center=(.35,-.35,1),radius=.3)])

##Sphere pack in box shape
from yade import pack, timing
readParamsFromTable(num=12000)
sp=pack.SpherePack()
sp.makeCloud((-.1,.1,0),(.8,-.8,.5),.03*((12000./yade.params.table.num)**(1/3.)),.5)

##Wall
O.bodies.append(wall((0,0,0),axis=2))

##Engines and Constitutive Law
O.engines=[ForceResetter(),
	InsertionSortCollider([Bo1_Sphere_Aabb(),Bo1_Sphere_Aabb(),Bo1_Wall_Aabb()],),
	InteractionLoop([Ig2_Sphere_Sphere_ScGeom(),Ig2_Sphere_Sphere_ScGeom(),Ig2_Wall_Sphere_ScGeom()],
		[Ip2_FrictMat_FrictMat_FrictPhys()],[Ip2_FrictMat_FrictMat_FrictPhys()],[Law2_ScGeom_FrictPhys_CundallStrack()]),
		NewtonIntegrator(damping=0.1,gravity=[0,0,-9.81])]

O.dt=.01*PWaveTimeStep()

#O.save('/tmp/c.xml.bz2');

I hope to get some advice on how to solve this problem of mine.
Thank you very much in advance!

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