yade-users team mailing list archive

[fengjingyu <question677261@xxxxxxxxxxxxxxxxxxxxx>]

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

fengjingyu posted a new comment:
>> How many cores does your processor have?

My computer has 16 processor and i allocated 12 processors to the
virtual machine.

>>Can you provide the output/numbers?

After 2100 iter,yade -j4 a.py use 3441112582us  and yade a.py use
3100565993us.It's strange.

My code is :

#####################
# encoding: utf-8
from yade import pack, qt, plot
import matplotlib; matplotlib.rc('axes',grid=True)
from yade import pack
import pylab
import yade.timing;
O.timingEnabled=True
yade.timing.stats()

sigmaIso=-25000

O.periodic=True

spheres=O.materials.append(FrictMat(young=64e9,poisson=0.12,density=2650,frictionAngle=0.231))
s=O.materials.append(FrictMat(young=64e9,poisson=0.12,density=2650,frictionAngle=0.24))

psdSizes,psdCumm=[0.1,0.12,0.15,0.19,0.20,0.23,0.27,0.28,0.30,0.32,0.33,0.50],[0.0001,0.003,0.03,0.11,0.25,0.43,0.70,0.85,0.95,0.97,0.98,1]

sp=pack.SpherePack()
sp.makeCloud((0,0,0),(4,4,4),psdSizes=psdSizes,psdCumm=psdCumm,distributeMass=True,periodic=True)
n=sp.toSimulation(material=spheres)

O.engines=[
   ForceResetter(),
   InsertionSortCollider([Bo1_Sphere_Aabb()]),
   InteractionLoop(
      [Ig2_Sphere_Sphere_ScGeom()],
      [Ip2_FrictMat_FrictMat_FrictPhys()],
      [Law2_ScGeom_FrictPhys_CundallStrack()]
   ),
   PeriTriaxController(label='triax',
      goal=(sigmaIso,sigmaIso,sigmaIso),stressMask=7,
      dynCell=True,maxStrainRate=(0.5,0.5,0.5),
      maxUnbalanced=1,relStressTol=0.05,
      doneHook='compactionFinished()'
   ),
   NewtonIntegrator(damping=.1),
   PyRunner(command='addPlotData()',iterPeriod=30),
]
O.dt=1*PWaveTimeStep()
O.trackEnergy=True

def particleNumber():
   print(O.bodies[-1].id)

   
def addPlotData():
   plot.addData(unbalanced=unbalancedForce(),i=O.iter,
      porosityy=voxelPorosity(800,(2,2,2),(3,3,3)),coordNum=avgNumInteractions()
      
   )
   
  
   print('p',voxelPorosity(800,(2,2,2),(3,3,3)),'s',triax.stress[0],triax.stress[1],triax.stress[2],'t',O.iter,yade.timing.stats())

pylab.semilogx(*sp.psd(bins=30,mass=True),label='Mass PSD of (free) %d random spheres'%len(sp))
pylab.legend()
pylab.show()

plot.plots={'i':('unbalanced',),'i ':('sxx','syy','szz'),' i':('exx','eyy','ezz'),' i ':('porosityy',),}
plot.plot()

qt.View()
qt.Controller()
O.saveTmp()

def compactionFinished():
   triax.maxUnbalanced=10
   triax.doneHook='triaxFinished()'


def triaxFinished():
   print 'Finished'
   O.pause()

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