yade-users team mailing list archive

[loiseaurare <question629227@xxxxxxxxxxxxxxxxxxxxx>]

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

Hi everybody,

I'm having a bit of a trouble here with the TriaxialCompressionEngine. Sad thing is, I'm sure it is stupid, but I really can't get round the problem.

I have adapted the agglomerate examples from github, in order to create a packed sample of agglomerates ( which it already did, but in three different simulations, and I wanted just one).

I am having trouble with the compaction process, it seems like the triax would simply not work.
I ran the simulation, and found the sample had'nt been packed. I have modifies the code in order that when I run the simulation, it will produce the agglomerates, and I could press run so that I see if anything is going wrong or happening. Well, seems that the walls don't move.

I tried activating the walls, like in the triax example (scenery 2), but that didn't change much.
Do you think maybe it is the clump function ? I haven't found that O.bodies.clump is much used, you can find much more on appendClumped, so I wondered if maybe the clump function did not fix particles in space or something...

the code is there :

# -*- coding: utf-8 -*-

from yade import export
from yade import export,ymport
import random


	
dim    = (15,15,15) # dimensions for makeCloud
radius = 1
fuzz   = 0.2

# use of makeCloud function
sp = pack.SpherePack()
sp.makeCloud((0,0,0), dim, rMean=radius, rRelFuzz=fuzz, seed=1)



random.seed(1) # to make colors always the same

attrs = []

colors = [randomColor() for s in sp]
# each macroparticle is filled randomDensePack
for si,s in enumerate(sp):
	sphere = pack.inSphere(s[0],s[1])
	sp1 = pack.randomDensePack(
		sphere,
		spheresInCell = 500,
		radius = .2,
		memoizeDb = '/tmp/agglomeratepackaux.db',
		returnSpherePack = True,
	)
	ids = sp1.toSimulation(color=colors[si]) # add the result to simulation with uniform color
	

	for i in ids:
		O.bodies[i].agglomerate = si # tell each particle who is its agglomerate
		attrs.append(O.bodies[i].agglomerate)

#random.seed(1)

sp2 = []
for b in O.bodies :
	if isinstance(b.shape,Sphere):
		sp2.append(b)

# add walls first
dim = (15,15,15)
walls = aabbWalls(((0,0,0),(dim)))
wallIds = O.bodies.append(walls)

# clump agglomerates
n = max(a for a in attrs)+1
colors = [randomColor() for _ in xrange(n)]
agglomerates = [[] for _ in xrange(n)]


for s,a in zip(sp2,attrs):
        aa = a
	s.agglomerate = aa
	s.shape.color = colors[aa]
        agglomerates[aa].append(s.id)

for g in agglomerates:
	O.bodies.clump(g)

triax=TriaxialStressController(
	maxMultiplier=1.005,
	finalMaxMultiplier=1.002,
	thickness = 0,
	stressMask = 7,
	internalCompaction=False,
	
)

triax.wall_bottom_activated=True
triax.wall_top_activated=True
triax.wall_left_activated=True
triax.wall_right_activated=True
triax.wall_back_activated=True
triax.wall_front_activated=True


O.engines = [
	ForceResetter(),
	InsertionSortCollider([Bo1_Sphere_Aabb(),Bo1_Box_Aabb()]),
	InteractionLoop(
		[Ig2_Sphere_Sphere_ScGeom(),Ig2_Box_Sphere_ScGeom()],
		[Ip2_FrictMat_FrictMat_FrictPhys()],
		[Law2_ScGeom_FrictPhys_CundallStrack()]
	),
	triax,
	TriaxialStateRecorder(iterPeriod=100,file='WallStresses'),
	NewtonIntegrator(damping=.6),
]


O.dt = PWaveTimeStep()
O.usesTimeStepper=True

# compress the sample
triax.goal1 = triax.goal2 = triax.goal3 = -1e-7

#O.run(50000,True) 
	
#triax.goal1 = triax.goal2 = triax.goal3 = -1e-5
#O.run(30000,True)




# save the result, including information of agglomerates which the particle belongs to
export.textExt('/tmp/compressed.txt','x_y_z_r_attrs',attrs=['b.agglomerate'])

try:
	from yade import qt
	qt.View()
except:
	pass


def getClumpInfo():
	for b in O.bodies:
		if b.isClump:
			print 'Clump ',b.id,' has following members:'
			keys = b.shape.members.keys()
			for ii in range(0,len(keys)):
				print '- Body ',keys[ii]
			print 'inertia:',b.state.inertia
			print 'mass:',b.state.mass,'\n'
		else: 
			print 'no clumps'

Thanks in advance to anybody that can help !


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