yade-users team mailing list archive

[Seti <question329634@xxxxxxxxxxxxxxxxxxxxx>]

New question #329634 on Yade:
https://answers.launchpad.net/yade/+question/329634

Hi All, 


I have tried to model a simple compression test as per below. and did  some parametric studies - Obviously by increasing the  \ values of frictional angle, the amount of stress  should be increased  however it does not happen ( e.g. frictional angle have been changes from 5 degree to 35 degree)  


I am sure now there is something wrong with script. maybe model of packing,..... 
Would you please advise? 


Here is the script: 

from yade import pack,plot

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

pred=pack.inCylinder((0,0,0),(0,0,0.5),0.07)

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()',iterPeriod=100,label='plotDataCollector',initRun=True)
]


#######################################3


##################################################


bo1s.aabbEnlargeFactor=1.0
ig2ss.interactionDetectionFactor=1.0
O.trackEnergy=True
plot.plots={'eps':('sigma',)} #,'sigma.50')},'t':('eps')} #'sigma.25','sigma.50','sigma.75')}

O.saveTmp('initial');

print 'Number of elements: ', len(O.bodies)

print 'normalCohesion=-3.5e1,shearCohesion=-10e1,phi:555,(0,0,0),(0,0,0.5correct),0.07)'


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

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

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

O.run(500000000,True)
plot.plots={'eps':('sigma')}
plot.plot(subPlots=False)
plot.saveDataTxt('uniax10e1,phi:5final.txt.bz2')



Thanks for your time. 


Seti

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