yade-users team mailing list archive
-
yade-users team
-
Mailing list archive
-
Message #10394
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.