save in txt.file a packing of spheres

Hello,

> i m new in yade

welcome :-)

> but is not working.

please always specify what "is not working" mean (I got his and that error, the position saved is not the final one...).

Also please try to provide working examples. In this case I got "KeyError: "No material labeled `ghiaia'."" and some NameErrors afterwards..

Yade uses Python scripts to control simulation. Python reads and executes line by line. So in your case, it saves the positions before anything is run. You have several options:

1) run some time steps
O.run(2000,True) # note the True second parameter

2) do the save in a function:
def checkUnbalanced():
   if unbalancedForce()<.05:
      O.pause()
      plot.saveDataTxt('bbb.txt.bz2')
      export.textExt('/tmp/compressed.txt','x_y_z_r_attrs',attrs=[])
      # note the indentation, now it is part of checkUnbalanced function

cheers
Jan

i was meaning that the code doesn't save me any positions, i see the simulation running. it ends and save the plots value inside 'bbb.txt.bz2, but any file txt for the positions and radii is created: i add the line code you suggest me but it still doesn't save nothing.

all the code is this now, it s really like deposition gravity, but i have to save the positions of the spheres in a txt file at the end of simulation:

from yade import pack, plot

# create rectangular box from facets
O.bodies.append(geom.facetBox((10,10,10),(10,10,10),wallMask=31))

sp=pack.SpherePack()

sp.makeCloud((0,0,0),(20,20,20),rMean=.3,rRelFuzz=.3)

sp.toSimulation()

O.engines=[
   ForceResetter(),
   InsertionSortCollider([Bo1_Sphere_Aabb(),Bo1_Facet_Aabb()]),
   InteractionLoop(
      # handle sphere+sphere and facet+sphere collisions
      [Ig2_Sphere_Sphere_ScGeom(),Ig2_Facet_Sphere_ScGeom()],
      [Ip2_FrictMat_FrictMat_FrictPhys()],
      [Law2_ScGeom_FrictPhys_CundallStrack()]
   ),
   NewtonIntegrator(gravity=(0,0,-9.81),damping=0.4),
   # call the checkUnbalanced function (defined below) every 2 seconds
   PyRunner(command='checkUnbalanced()',realPeriod=2),
   # call the addPlotData function every 200 steps
   PyRunner(command='addPlotData()',iterPeriod=100)
]
O.dt=.5*PWaveTimeStep()

O.trackEnergy=True
O.run(2000,True)

def checkUnbalanced():
   if unbalancedForce()<.05:
      O.pause()
      plot.saveDataTxt('bbb.txt.bz2')
      export.textExt('/tmp/compressed.txt','x_y_z_r_attrs',attrs=[])

# collect history of data which will be plotted
def addPlotData():

   plot.addData(i=O.iter,unbalanced=unbalancedForce(),**O.energy)

plot.plots={'i':('unbalanced',None,O.energy.keys)}

plot.plot()

O.saveTmp()

I have tried your code and the data was saved.. you chose the file for saving to be /tmp/compressed.txt, maybe this confused you? After running the script, what is the output of terminal command
cat /tmp/compressed.txt
?
thanks
Jan

Thanks jan for The answer, in The end i solved My problem , i created a functions that save the positions of The particles every 100 iter period. Now i m facing a new problem . I need To find And eventually change The Kn And Ks of The material. In The command FrictMat i can only choose young And poisson( i read that in realty that value aren t The real poisson And young Number) how can i do? Thanks

Thanks Luc Scholtès, that solved my question.