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.
_______________________________________________
Mailing list: https://launchpad.net/~yade-users
Post to : [email protected]
Unsubscribe : https://launchpad.net/~yade-users
More help : https://help.launchpad.net/ListHelp
