yade-users team mailing list archive

[Othman Sh <question695989@xxxxxxxxxxxxxxxxxxxxx>]

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

    Status: Answered => Open

Othman Sh is still having a problem:
Hi Jan,

I have updated my code to use utils.box instead of facet. My question is
how to apply force on the spheres without using a "plate" facet? I am
currently using geom.facetBox to make a flat surface "i.e. plate" and
move it downward to apply force, but still smaller spheres can pass
through it.

Also, how to make a packing with 2 psdSizes and psdCumm while using one
packSpherepack() and one toSimulation to avoid initial overlaping of
spheres?

Thanks,
Othman

---------
from yade import pack
import numpy as np


Mat1=O.materials.append(FrictMat(young = 1e6, poisson = 0.45,frictionAngle = 0.00349066, density=2340,label='portlandite'))
Mat2=O.materials.append(FrictMat(young = 7.2e10, poisson = 0.17,frictionAngle = 0.558, density=2650,label='sand')) 
Mat3=O.materials.append(FrictMat(young = 1e20, poisson = 0.1))


cube_dimensions=100e-6 #micron

##Initial cube dimensions ###
mnx=0
mny=0
mnz=0
mxx=cube_dimensions
mxy=cube_dimensions
mxz=cube_dimensions

box_center=(mxx/2,mxy/2,mxz/2)
box_extents=(mxx/2,mxy/2,mxz/2)
plate_center=(mxx/2,mxy/2,1.05*mxz/2)
plate_extents=(mxx/2,mxy/2,0.98*mxz/2)
############################ spheres #############################
sp1=pack.SpherePack()

sizes1=1e-6*np.array([9,9.16,15.76,22.94,25]) 
passing1=[0,0.5,0.75,0.9,1]
sp1.makeCloud((mnx,mny,mnz),(mxx,mxy,mxz),psdSizes=sizes1,psdCumm=passing1,num=5811)
sp1.toSimulation(color=(0.533, 0.803, 0.780),material=Mat1)


sizes2=1e-6*np.array([149,150,300,425,600,850]) #Diameters of portlandite
passing2=[0,0.02,0.23,0.675,0.98,1]
sp1.makeCloud((mnx,mny,mnz),(mxx,mxy,mxz),psdSizes=sizes2,psdCumm=passing2,num=51)
sp1.toSimulation(color=(0.152, 0.368, 0.988),material=Mat2)

############################# box and plate #############################
yade.qt.View()

box=O.bodies.append(utils.box(box_center,box_extents,mask=31,material=Mat3,wire=True))
plate1IDs=O.bodies.append(geom.facetBox(plate_center,plate_extents,material=Mat3,wallMask=32))

for i in plate1IDs:
 body= O.bodies[i]
 body.state.vel = (0,0,-10)

############################## compaction #############################
O.dt=.005*utils.PWaveTimeStep()

O.engines=[
 ForceResetter(),
 InsertionSortCollider([
  Bo1_Sphere_Aabb(),
  Bo1_Facet_Aabb()
 ]),
 InteractionLoop(
  [
   Ig2_Sphere_Sphere_ScGeom(),
   Ig2_Facet_Sphere_ScGeom(),
  ],
  [

   Ip2_FrictMat_FrictMat_FrictPhys(),
   Ip2_FrictMat_FrictMat_FrictPhys(),
  ],
  [

   Law2_ScGeom_FrictPhys_CundallStrack(),
  ],
 ),

 NewtonIntegrator(damping=.3),

]
#O.run()

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