yade-users team mailing list archive

["m.asd" <question215627@xxxxxxxxxxxxxxxxxxxxx>]

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

hi,
I am a new user of yade and i am going to model a box in which some particles deposit. then I delete particles upper than a predefined elevation and add some fine particles. then I delete these fines upper than that elevation too. after that, I  add a plate on top and define a velocity up to a specific load. set the  velocity of the plate to zero (fix place) and then let sidewalls to rotate around the Y axis with a fixed angVel.
here is my problem,
when I use "plate.state.vel" it says that "Nonetype object has no  attribute state"
and when I use rotationengine and translationengine it shows an error and do not proceed.
as I am not so experienced in running yade, changing the script in many ways made me frustrated. I will be so grateful if anyone helps me.
furthermore, It sometimes stops running with this error
"segmentation fault (core dumped)" what is the reason
here is my script


overburden=10
from yade import pack,plot,qt
idmat=O.materials.append(FrictMat(density=2700,young=1e6,poisson=0.5,frictionAngle=0))
O.bodies.append(utils.wall(position=0,axis=2,sense=1))
O.bodies.append(utils.wall(position=0,axis=1,sense=1))
O.bodies.append(utils.wall(position=0.01,axis=1,sense=-1))
O.bodies.append(utils.wall(position=0,axis=0,sense=1))
LW=O.bodies[-1]
O.bodies.append(utils.wall(position=0.01,axis=0,sense=-1))
RW=O.bodies[-1]
sp=pack.SpherePack()
sp.makeCloud((0,0,0.001),(0.01,0.01,0.02),rMean=.0008,rRelFuzz=0)
sp.toSimulation(color=(0,0,1)) # blue

O.engines=[
	ForceResetter(),
	InsertionSortCollider([Bo1_Sphere_Aabb(),Bo1_Wall_Aabb()]),
	InteractionLoop(
		[Ig2_Sphere_Sphere_ScGeom(),Ig2_Wall_Sphere_ScGeom()],
		[Ip2_FrictMat_FrictMat_FrictPhys()],
		[Law2_ScGeom_FrictPhys_CundallStrack()]
	),
	GravityEngine(gravity=(0,0,-9.81)),
	NewtonIntegrator(damping=0.4),
	qt.SnapshotEngine(fileBase='3d-',iterPeriod=100000,label='snapshot'),
	PyRunner(command='unbalancedData()',iterPeriod=30000),
	PyRunner(command='checkUnbalanced()',realPeriod=1,label='checker'),
]

O.dt=1e-6

def checkUnbalanced():
	if O.iter<500000: return
	if utils.unbalancedForce()>.08: return
	for b in O.bodies:
		if isinstance(b.shape,Sphere):
			if (b.state.pos[2]+b.shape.radius)>0.007: O.bodies.erase(b.id)
	fines=pack.SpherePack()
	fines.makeCloud((0,0,0.007),(0.01,0.01,0.01),rMean=.0004,rRelFuzz=0)
	fines.toSimulation()
	O.engines=O.engines+[GravityEngine(gravity=(0,0,-9.81),mask=-1)]
	checker.command='addfines()'

def addfines():
	if O.iter<650000: return
	if utils.unbalancedForce()>.08: return
	for b in O.bodies:
		if isinstance(b.shape,Sphere):
			if (b.state.pos[2]+b.shape.radius)>0.007: O.bodies.erase(b.id)
	posi=O.bodies[5].state.pos[2]+O.bodies[5].shape.radius
	m=5
	while m in range(len(O.bodies)):
		b=O.bodies[m]
		try:
			if isinstance(b.shape,Sphere):
				posi2=b.state.pos[2]+b.shape.radius
				if posi2>=posi: posi=posi2
				m=m+1
		except: m=m+1
	O.bodies.append(utils.wall(position=posi,axis=2,sense=-1)) 
	global plate
	plate=O.bodies[-1]
	O.engines=O.engines[:4]+[TranslationEngine(ids=plate,translationAxis=(0,0,1),velocity=0.002,label="transEngine")]+O.engines[4:]   #here is the problem
	checker.command='shearing()'
	O.materials[0].frictionAngle=.5 
	for i in O.interactions: i.phys.tangensOfFrictionAngle=tan(.5)

def shearing():
	if abs(O.forces.f(plate.id)[2])<overburden: return
	transEngine.velocity=0
	O.engines=O.engines[:5]+[RotationEngine(ids=[LW,RW],angularVelocity=0.005,rotationAxis=(0,1,0),label="rotEngine")]+O.engines[5:]      #and here too
	checker.command='stopShearing()'

def stopShearing():
	if O.bodies[3].state.rot<(0,0.005,0): return
	O.pause()
	utils.makeVideo(snapshot.snapshots,'3dshear.mpeg',fps=10,bps=10000)
	print 'Finished'

def unbalancedData():
	plot.addData(itperi=O.iter,unbalanced=utils.unbalancedForce())
	
plot.plots={'itperi':('unbalanced',)}
plot.plot() 

Gl1_Sphere.stripes=True
qt.View()

rr=yade.qt.Renderer()
rr.shape=True
O.saveTmp()


-- 
You received this question notification because you are a member of
yade-users, which is an answer contact for Yade.