yade-users team mailing list archive

[jamespaul <question678348@xxxxxxxxxxxxxxxxxxxxx>]

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

jamespaul posted a new comment:
Bruno,

I saw the oscillations.But the oscillation is purely in the first
isotropic compression and does not cause a reduction in pressure in the
second compression.

#############################
# encoding: utf-8

sigmaIso=-25000

from yade import pack, qt, plot

O.periodic=True

spheres=O.materials.append(FrictMat(young=64e9,poisson=0.12,density=2650,frictionAngle=0.231))

length=3

psdSizes=[0.15,.18,.2,.21,.25,.27,.30,.31,.34,.38]
psdCumm=[0.025,0.11,0.27,0.45,0.73,0.88,0.95,0.98,0.99,1]

sp=pack.SpherePack()
sp.makeCloud((0,0,0),(length,length,length),psdSizes=psdSizes,psdCumm=psdCumm,distributeMass=True,periodic=True)
n=sp.toSimulation(material=spheres)

O.engines=[
	ForceResetter(),
	InsertionSortCollider([Bo1_Sphere_Aabb(),Bo1_Facet_Aabb()]),
	InteractionLoop(
      		[Ig2_Sphere_Sphere_ScGeom()],
      		[Ip2_FrictMat_FrictMat_FrictPhys()],
      		[Law2_ScGeom_FrictPhys_CundallStrack()]
   	),
	PeriTriaxController(label='triax',
		goal=(sigmaIso,sigmaIso,sigmaIso),stressMask=7,
		dynCell=True,maxStrainRate=(.1,.1,.1),
		maxUnbalanced=.01,relStressTol=1e-3,
		doneHook='compactionFinished()'
   	),
	NewtonIntegrator(damping=.1),
	PyRunner(command='addPlotData()',iterPeriod=10000),
]
O.dt=.5*PWaveTimeStep()

def addPlotData():
	plot.addData(sxx=triax.stress[0],syy=triax.stress[1],szz=triax.stress[2],
        	exx=triax.strain[0],eyy=triax.strain[1],ezz=triax.strain[2],
		q=triax.stress[2]-0.5*(triax.stress[0]+triax.stress[1]),
		p=(triax.stress[0]+triax.stress[1]+triax.stress[2])/3,
		unbalanced=unbalancedForce(),i=O.iter,
		#k=(triax.stress[0]*0.5+triax.stress[1]*0.5)/triax.stress[2] 
 	)
	print('szz=',triax.stress[2],'ezz=',triax.strain[2])
	#plot.saveDataTxt('qp1.txt',vars=('q','p','szz','ezz'))

plot.plots={'i':('unbalanced',),'i ':('sxx','syy','szz'),' i':('exx','eyy','ezz'),}
plot.plot()
	

def compactionFinished():
	print('aaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaa')
	O.cell.trsf=Matrix3.Identity
	s=O.materials.append(FrictMat(young=64e9,poisson=0.12,density=2650,frictionAngle=0.24))
	for i in n:
        	O.bodies[i].material = O.materials[s]
	triax.goal=(sigmaIso,sigmaIso,-0.005)
	dynCell=False,
	triax.stressMask=3
	triax.maxStrainRate=(0,0,.001)
	triax.doneHook='triaxFinished()'
	triax.maxUnbalanced=0.000001

def triaxFinished():
	print 'Finished'
	O.pause()

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