yade-users team mailing list archive

[Christian Sommerfeld <question203709@xxxxxxxxxxxxxxxxxxxxx>]

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

Hi,

I want to create a simulation with the basic DEM model from Cundall and Strack. But i don find any calculation for it. 
Therefore i have the following question:
Is it correct to input  young's modulus, poission's ratio and the friction angle for the law of cundall and strack and the density for the calculation of the gravity force?
I have done it in this way: 
O.materials.append(FrictMat(young=5e4,poisson=0.5,frictionAngle=.6,density=2400))

with the following engine:
O.engines=[
   ForceResetter(),
   InsertionSortCollider([Bo1_Sphere_Aabb(),Bo1_Facet_Aabb()]),
   InteractionLoop(
      [Ig2_Sphere_Sphere_Dem3DofGeom(),Ig2_Facet_Sphere_Dem3DofGeom()],
      [Ip2_2xFrictMat_CSPhys()],
      [Law2_Dem3Dof_CSPhys_CundallStrack()]
   ),
   GravityEngine(gravity=(0,0,-9.81)),
	RotationEngine(rotateAroundZero=True,zeroPoint=(0,0,0),rotationAxis=(0,0,1),angularVelocity=50*(2*pi/60),ids=quader,label='rotor'),
   NewtonIntegrator(damping=0.4),
 PyRunner(iterPeriod=20,command="myAddPlotData()"),
	VTKRecorder(iterPeriod=100,fileName='/tmp/QuaderAufnahme-',recorders=['spheres','facets'])
]

Furthermore i am still wondering why the simulation only apply the last convention of my material, when i write it in this way:

FacetMat=O.materials.append(FrictMat(young=5e4,poisson=0.5,frictionAngle=.6,density=2400))
SphereMat=O.materials.append(FrictMat(young=5e10,poisson=0.5,frictionAngle=.6,density=10000))

oriBody = Quaternion(Vector3(0,0,1),(math.pi))
quader=O.bodies.append(geom.facetBox((.5,0,.5),(.3,.3,.3),oriBody, wallMask=63, color=(0,0,1), dynamic=True, material=FacetMat))

sphereRadius = 0.06
nbSpheres = (10,10,30)
for i in xrange(nbSpheres[0]):
    for j in xrange(nbSpheres[1]):
        for k in xrange(nbSpheres[2]):
            x = (i*2 - nbSpheres[0])*sphereRadius*1.1
            y = (j*2 - nbSpheres[1])*sphereRadius*1.1
            z = (k*2 - nbSpheres[2])*sphereRadius*1.1
            sp1=utils.sphere([x,y,z+10],sphereRadius,color= (1,0,0),material=SphereMat)
            O.bodies.append(sp1)




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