← Back to team overview

yade-users team mailing list archive

Re: [Question #255528]: uniaxial compression test of particulate material

 

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

    Status: Answered => Open

Yuan is still having a problem:
Hello Jérôme,

Thank you for your advice last time.

I have updated my YADE to 1.12.0 and that Attribute Error has been
solved.

However now I can not get any plot from my script.  I did not change my
script so it was the same one I used with my old version of YADE which
is 1.07 and in that version the plot was ok.

I checked the terminal while I ran my script and there was no error
shown.  So my script was unchanged and no error shown in terminal
however the plot was gone.

I attached my script

from yade import pack,plot

idParticulate=O.materials.append(CohFrictMat(young=30e9,poisson=0.2,alphaKr=3000,alphaKtw=3000,density=2700,frictionAngle=0.6,isCohesive=True,momentRotationLaw=True,normalCohesion=-3.5e6,shearCohesion=-10e6,etaRoll=-3,etaTwis=-3))

pred=pack.inCylinder((0,0,0),(0,0,0.2),0.05)

TS=pack.randomDensePack(pred,radius=0.008,material=idParticulate)
O.bodies.append(TS)

bb=uniaxialTestFeatures()
negIds,posIds,axis,crossSectionArea=bb['negIds'],bb['posIds'],bb['axis'],bb['area']

O.dt=1e-5*PWaveTimeStep()


O.engines=[
ForceResetter(),
InsertionSortCollider([Bo1_Sphere_Aabb(aabbEnlargeFactor=1.5,label='bo1s')]),
InteractionLoop(
[Ig2_Sphere_Sphere_ScGeom6D(interactionDetectionFactor=1.5,label='ig2ss')],
[Ip2_CohFrictMat_CohFrictMat_CohFrictPhys(setCohesionNow=True,setCohesionOnNewContacts=False)],
[Law2_ScGeom6D_CohFrictPhys_CohesionMoment()]
),
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)
]

bo1s.aabbEnlargeFactor=1.0
ig2ss.interactionDetectionFactor=1.0

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

def addPlotData():
       yade.plot.addData(t=O.time,i=O.iter,eps=-strainer.strain,sigma=-strainer.avgStress)

O.step()


Yuan

-- 
You received this question notification because you are a member of
yade-users, which is an answer contact for Yade.