yade-users team mailing list archive

[Alireza Sadeghi <question688073@xxxxxxxxxxxxxxxxxxxxx>]

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

Hello All,

I want to use clump shape particles instead of spherical particles in the PeriTriaxController. As you can see in the following code, there is only one particle in the middle of a cell. There is no interactions between the boundaries and the particle at the begininng. Why the stress at the boundaries are not zero?
I asked this question several times, and nobody answered me. Therfore, the questions were expired and removed. I have stopped in this position for two months. I will be so grateful if you could help me.

Thank you very much for your help in advance.

Best Regards,

Alireza

######code######


from yade import utils, plot
from yade import pack, qt
from datetime import datetime

#==================set a periodic boundary==================

O.periodic=True
O.cell.refSize=(10e-3,10e-3,10e-3)

#================= define the materials =======================

#O.interactions.clear()

O.materials.append(CohFrictMat(normalCohesion= 1e20, shearCohesion= 1e20, isCohesive= True, young=6.81e8, density=1377.5, poisson=0.3, frictionAngle= 0.31, fragile=False, label='aggregate-48'))

#===================clump shape=================================
	
radz1=[3e-3,3e-3]
poz1= [[0,0,-1e-3],[1e-3,0,0]]		 
template1= []
template1.append(clumpTemplate(relRadii=radz1,relPositions=poz1))


nums=['t']

mats=['aggregate-48']

coke=(1e-3,1)

color=((0,0,1))

tolerance=[(2e-4)]

temps=[template1]

nums=pack.SpherePack()

nums.makeCloud((3e-3,3e-3,3e-3),(7e-3,7e-3,7e-3),rMean=coke[0],rRelFuzz=tolerance[0],num=coke[1])

O.bodies.append([utils.sphere(c,r,material=mats[0],color=color[0]) for c,r in nums])
O.bodies.replaceByClumps(temps[0],[1.0],discretization=5)

#O.interactions.clear()

#====================Engine====================

sigmaIso=-1e5

O.engines=[
     ForceResetter(),
     InsertionSortCollider([Bo1_Sphere_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=(10,10,10),
		maxUnbalanced=.01,relStressTol=1e-3,
		doneHook='compactionFinished()'
     ),
     NewtonIntegrator(damping=0.4,gravity=[0,0,0]),

     PyRunner(command='history()',iterPeriod=1,label='recorder'),
]

O.dt=5e-8

def history():
	plot.addData(unbalanced=unbalancedForce(),
		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]-triax.stress[1]), p=-(triax.stress[2]+triax.stress[1]+triax.stress[0])/3,
                por=porosity(),i=O.iter,)

        print 'porosity:', porosity(),'Unbalancing:',unbalancedForce(),'Strain',-triax.strain[2], 'Stress',-triax.stress[2]

def compactionFinished():
	O.pause()


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