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

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. If this script of mine is 
working properly, the sequence of simulation should be 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.

_______________________________________________
Mailing list: https://launchpad.net/~yade-users
Post to     : [email protected]
Unsubscribe : https://launchpad.net/~yade-users
More help   : https://help.launchpad.net/ListHelp

Reply via email to