yade-users team mailing list archive

[Swapnil <question655866@xxxxxxxxxxxxxxxxxxxxx>]

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

Swapnil posted a new comment:
Hi Jan, thanks for the response.

I am trying a simpler case for the uniaxial test. Following things up
with the following code:

from yade import pack,plot
pred=pack.inCylinder((0,0,0),(0,0,0.2),0.05)  (Packing in a cylinder)
sp=pack.randomDensePack(pred,radius=0.008)  (Packing spheres in the above cylinder)
O.bodies.append(sp)  (Adds the spheres to the cylinder geometry)
yade.qt.Controller() (Pressed F12 and it displays the packed geometry)

Moving towards uniaxial testing:

bb=uniaxialTestFeatures()

negIds,posIds,axis,crossSectionArea=bb['negIds'],bb['posIds'],bb['axis'],bb['area']
--> this time the command worked unlike the last time when I got an
error. Still no idea what the reason might have been

O.dt=0.005*PWavetimeStep()   (Time step control for the simulation)

Next, we go for the Interaction of spheres inside:

O.engines=[ForceResetter(),InsertionSortCollider([Bo1_Sphere_Aabb()]),InteractionLoop([Ig2_Sphere_Sphere_ScGeom()],[Ip2_FrictMat_FrictMat_FrictPhys()],[Law2_ScGeom_FrictPhys_CundallStrack()]),NewtonIntegrator(),UniaxialStrainer(strainRate=-0.00005,axis=axis,asymmetry=0,posIds=posIds,negIds=negIds,crossSectionArea=crossSectionArea,blockDisplacements=False,blockRotations=False,setSpeeds=False,label='strainer'),PyRunner(virtPeriod=1e-6/0.5,realPeriod=1,command='addPlotData()',label='plotDataCollector',initRun=True)]
--> Here I want to understand the PyRunner (........) command. Can you
describe this.

Next is the plotting of data:

plot.plots={'eps':('sigma')}

def addPlotData():
     ....:     yade.plot.addData(t=O.time,i=O.iter,eps=-strainer.strain,sigma=-strainer.avgStress)
     ....:     O.step()                      
--> When I execute the above, it shows "addPlotData" is not defined. 

the sim ulation runs and I can see the straining of the sample. Upon
pressing F8 to see the plot of sigma against the epsilon (stress-strain)
I see nothing (no plots).  What am I doing wrong or what changes can I
make to execute this task in a better way?

I hope I could make everything clear...

Thanks :)

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