yade-users team mailing list archive

[anna effeindzourou <question197520@xxxxxxxxxxxxxxxxxxxxx>]

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

anna effeindzourou posted a new comment:
Hi,
I have prepared a short script .Maybe the problem is due to the fact that the is an overlap between spheres, when the randomdensepack is generate.
The contact stiffness is really important. I am trying to obtain a stress-free packing.
Thanks,
Anna

Here is the script.

# -*- coding: utf-8
from yade import ymport, utils,pack


O.materials.append(FrictMat(young=16e9,poisson=0.25,frictionAngle=27,density=1576.))

O.bodies.append(geom.facetBox(center=(0.0,0.0,0),extents=(.2,.2,.1),orientation=Quaternion((0.,0,0),0),wallMask=31,wire=True,dynamic=False,color=(0,1,1)))

pred=pack.inAlignedBox((-.2,-.2,-.1),(.2,.2,.1))
spheres=pack.randomDensePack(pred,spheresInCell=2000,radius=0.015,rRelFuzz=.5)
O.bodies.append(spheres)

O.engines=[
	ForceResetter(),
	InsertionSortCollider([Bo1_Sphere_Aabb(label='aabb'),Bo1_Facet_Aabb()]), 
	InteractionLoop(
	[Ig2_Sphere_Sphere_ScGeom(label='Ig2ssGeom'),Ig2_Facet_Sphere_ScGeom(label='Ig2fsGeom')],
	[Ip2_FrictMat_FrictMat_FrictPhys(label='block_particles'),
	],
	[Law2_ScGeom_FrictPhys_CundallStrack(label='Law_2'),
	]),
	NewtonIntegrator(gravity=[0,0,-9.81],damping=0.5,label='Newton_Integrator'),
]

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