yade-users team mailing list archive

[Chiara Modenese <question203709@xxxxxxxxxxxxxxxxxxxxx>]

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

    Status: Open => Answered

Chiara Modenese proposed the following answer:
Hi Christian,

Could you please (briefly) explain here what material you aim to
simulate and what kind of test you want to run with it? It would help us
to address your problem.

Chiara


On 20 Jul 2012, at 11:51, Christian Sommerfeld wrote:

> 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.
> 
> _______________________________________________
> Mailing list: https://launchpad.net/~yade-users
> Post to     : yade-users@xxxxxxxxxxxxxxxxxxx
> Unsubscribe : https://launchpad.net/~yade-users
> More help   : https://help.launchpad.net/ListHelp

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