ScriptGenerateBoxPacking
From Yade
<source lang="python"> p=Preprocessor('TriaxialTest') p['recordIntervalIter']=0 # don't record anything p['numberOfGrains']=2000; p['lowerCorner']=[0,0,0] p['upperCorner']=[8,16,8] # this is just size ratio, since radiusMean is specified p['radiusStdDev']=0; p['radiusMean']=3e-4
p['sigmaIsoCompaction']=1e7 # compaction pressure p['sigmaLateralConfinement']=1e3 # unloading pressure
p['StabilityCriterion']=.05 p['strainRate']=.2 p['maxWallVelocity']=.1 p['wallOversizeFactor']=2;
p['autoUnload']=True # unload after isotropic compaction p['autoCompressionActivation']=False # stop once unloaded
p['defaultDt']=1e-4
p.load() # generate scene and load it
o=Omega() o.run() ## triax stops by itself once unloaded o.wait()
- save spheres
from yade import utils utils.spheresToFile('triax-gen-%d.spheres'%(len(o.bodies)-6)) # 6 walls </source>