Obtaining Volume Fraction of Filled Cylinder
Hey guys! So I have been having trouble finding a way in which to obtain the volume fraction for a filled coaxial cylinder model of spheres using utils.getDepthP
Thank you for any help or advice given!
#Module Setup
from yade import ymport, plot
import time
#Constants
utils.readParam
frIn = 0.5, enIn=0.01, etIn=0.01, tcIn=0.0001, # Frictionangle (rad), strength-parameter
rhoIn = 2300.0, # Density
dumpVTKIn = 100000 # Periods of dumps in iterations
)
#Sweep Parameters using Table
from yade.params.table import
import shutil
try:
shutil.
except OSError:
pass
os.mkdir('step1')
folderNameBase = 'step1/' + str(descriptionIn)
folderName = folderNameBase
os.mkdir(
#Initialize
o = Omega()
o.dt = 0.05*1E-7
#Define Material Properties of Spheres
mat1 = o.materials.
#Build Geometry
cylinder_
radiusTopIn
cylinder_
segmentsNum
cylinder=
o.engines = [
#Following Default Workflow
ForceResett
InsertionSo
# SpatialQuickSor
Interaction
),
NewtonInteg
#Particles will drop from a circular plane
CircularFac
radius = (3.5e-3),
PSDsizes = [125e-6,
PSDcum = [0.0250,
),
#Delete Particles that fly out of domain
DomainLimit
#Shake the domain to help settling
HarmonicMot
#Record output to VTK Files, Spheres,
VTKRecorder
VTKRecorder
#FIXME Snapshot Taker
# SnapshotEngine(
#Execute Python Commands
PyRunner(
# PyRunner(
PyRunner(
]
#Enable Track Energy - Can be used to determine if sim is done
o.trackEnergy=True
#Volume Fraction Tracking
from yade import utils
utils.getDepthP
#Run Yade Simulation
o.stopAtIter=
o.run()
def stop_condition():
o.dt=
print(
Question information
- Language:
- English Edit question
- Status:
- Answered
- For:
- Yade Edit question
- Assignee:
- No assignee Edit question
- Last query:
- Last reply:
Can you help with this problem?
Provide an answer of your own, or ask Kevin for more information if necessary.