voxelPorosity
Hi,
I'm trying to set the porosity using voxelPorosity in gravity deposition, as in the script.
The value of friction angle starts in zero and when the target porosity is reached, then I set the friction angle to desired.
But, the statement "if utils.voxelPoro
Any idea?
#!/usr/bin/python
# -*- coding: utf-8 -*-
#Modelo para se obter a porosidade desejada em uma soil bin
#voxelPorosity
#######
### DEFINING VARIABLES AND MATERIALS #######
#######
#material parameters spheres
youngs = 70e9 #Pa
poissons = 0.23 #admensional
frictionangles = 0 #degrees
densitys = 1900 #kg/m3
#material parameters walls
youngw = 205e9 #Pa
poissonw = 0.3 #admensional
frictionanglew = 30 #degrees? entre esferas e parede ou entre parede?
densityw = 7850 #kg/m3
targetPorosity = 0.74
frictionanglesf = 30 #degrees
#Material
O.materials.
O.materials.
## create walls around the packing
O.bodies.
#Create sphere packing
sp=pack.
sp.makeCloud(
sp.toSimulation
#Engines
O.engines=[
ForceResetter(),
InsertionSor
InteractionLoop(
# handle sphere+sphere and facet+sphere collisions
[Ig2_Sphere_
[Ip2_FrictMat_
[Law2_
),
NewtonIntegr
]
#######
############### Calc Porosity ###############
#######
#O extrema pega as dimensoes do solido formado pelas esferas
O.engines=
if utils.voxelPoro
setContactFric
print "change in friction angle"
O.dt=1*
O.saveTmp()
Question information
- Language:
- English Edit question
- Status:
- Expired
- For:
- Yade Edit question
- Assignee:
- No assignee Edit question
- Last query:
- Last reply: