Re: [Yade-users] [Question #689407]: Are there any way to stop particles poping out of the box?

2020-03-28 Thread Jan Stránský
Question #689407 on Yade changed:
https://answers.launchpad.net/yade/+question/689407

Status: Open => Answered

Jan Stránský proposed the following answer:
> Do you have any idea how to stop top box?

See answer #1 
>> To make the box step, just assign it zero velocity:
>> top.state.vel = Vector3.Zero

> I want compaction ... at ... stage which volume equal to 1.0*1.0*1.0.
How to do it ?

use e.g. PyRunner:
###
O.engines = [
   ...
   PyRunner(iterPeriod=1,command="stopBoxIfCompacted()")
]
def stopBoxIfCompacted():
   area = ...
   height = ...
   volume = area * height
   if volume < 1.0*1.0*1.0:
  top.state.vel = Vector3.Zero
  O.pause() # ?
###

cheers
Jan

-- 
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 : yade-users@lists.launchpad.net
Unsubscribe : https://launchpad.net/~yade-users
More help   : https://help.launchpad.net/ListHelp


Re: [Yade-users] [Question #689407]: Are there any way to stop particles poping out of the box?

2020-03-27 Thread veerawat
Question #689407 on Yade changed:
https://answers.launchpad.net/yade/+question/689407

veerawat gave more information on the question:
Dear jan,

so suppose I want compaction of the particles at it’s most tightest
stage which volume equal to 1.0*1.0*1.0. How to do it ?

-- 
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 : yade-users@lists.launchpad.net
Unsubscribe : https://launchpad.net/~yade-users
More help   : https://help.launchpad.net/ListHelp


Re: [Yade-users] [Question #689407]: Are there any way to stop particles poping out of the box?

2020-03-27 Thread veerawat
Question #689407 on Yade changed:
https://answers.launchpad.net/yade/+question/689407

Status: Answered => Open

veerawat is still having a problem:
Dear jan ,
from this video : 
https://drive.google.com/file/d/1KHIYHUPM75RbGP0xjJt3yqvl9DiNfWUM/view?usp=drivesdk
 

Before the particles pop out of all the box, the particles reach its
compression peak point then pop out. I want to stop top box before the
particles pop out, Do you have any idea how to stop top box?

-- 
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 : yade-users@lists.launchpad.net
Unsubscribe : https://launchpad.net/~yade-users
More help   : https://help.launchpad.net/ListHelp


Re: [Yade-users] [Question #689407]: Are there any way to stop particles poping out of the box?

2020-03-20 Thread Jan Stránský
Question #689407 on Yade changed:
https://answers.launchpad.net/yade/+question/689407

Status: Open => Answered

Jan Stránský proposed the following answer:
Hi,

> Are there any function key to make top box stop when the particles is
tightly packed at the bottom of the box, how?

there are several approaches, depending on your specific definition of what 
"tightly packed" is, e.g.:
- you know desired packing fraction (the resulting box volume), so you can a 
priori tell when to stop the top box
- based on unbalancedForce, kinetic energy, whatever.. checking it periodically 
and stop when the desired state is reached

To make the box step, just assign it zero velocity:
top.state.vel = Vector3.Zero

> Are there any way to stop particles poping out of the box, how?

- increase stiffness of the box (use another material for them)
- or stop it before it pushes particles away (see above)

cheers
Jan

-- 
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 : yade-users@lists.launchpad.net
Unsubscribe : https://launchpad.net/~yade-users
More help   : https://help.launchpad.net/ListHelp


[Yade-users] [Question #689407]: Are there any way to stop particles poping out of the box?

2020-03-20 Thread veerawat
New question #689407 on Yade:
https://answers.launchpad.net/yade/+question/689407

I want to use top box to press down the scattered particles until the particles 
all be compressed to the bottom of the box. Then, the particles and top box 
should stop when they firmly fit right at the bottom of the box

And here are some difficulty I have faced;
The particles don’t stop when they all pack at the bottom but rather pop out or 
break throughout of the box, the top box itself also go out of the wall.

And here is my question; 
Are there any function key to make top box stop when the particles is tightly 
packed at the bottom of the box, how?
Are there any way to stop particles poping out of the box, how?

this is video :

https://drive.google.com/file/d/1KHIYHUPM75RbGP0xjJt3yqvl9DiNfWUM/view?usp=drivesdk

This is my script :

from yade import pack


O.materials.append(FrictMat(density=1000,young=1e4,poisson=0.3,frictionAngle=radians(30),label='sphereMat'))


# create rectangular box from boxes (maybe there is some library function, but 
I did not find it quickly)
cx,cy,cz = .5,.5,.5 # center of the box
dx,dy,dz = .5,.5,.5 # half-dimensions of the box
t = 0.025 # half-thickness of walls
topx,topy = 0.2, 0.3 # top box half-dimensions
left = box((cx-dx-2*t,cy,1.5*cz),(2*t,dy,1.5*dz),fixed=True,wire=True)
right = box((cx+dx+2*t,cy,1.5*cz),(2*t,dy,1.5*dz),fixed=True,wire=True)
front = box((cx,cy-dy-t,1.5*cz),(dx,t,1.5*dz),fixed=True,wire=True)
back = box((cx,cy+dy+t,1.5*cz),(dx,t,1.5*dz),fixed=True,wire=True)
bottom = box((cx,cy,cz-dz-5*t),(dx,dy,5*t),fixed=True,wire=True) #t


top = box((cx,cy,1.5*cz+1.5*dz+t),(dx,dy,t),fixed=True) 
O.bodies.append((left,right,back,front,bottom,top))


sp=pack.SpherePack()
# generate randomly spheres with uniform radius distribution
sp.makeCloud((0,0,0),(1,1,1.5),rMean=.03,rRelFuzz=.01,num=10,seed=1)
# add the sphere pack to the simulation
sp.toSimulation()


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()]
),
NewtonIntegrator(gravity=(0,0,-9.81),damping=0.4),


]
O.dt=.5*PWaveTimeStep()
top.state.vel=(0,0,-.5)



-- 
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 : yade-users@lists.launchpad.net
Unsubscribe : https://launchpad.net/~yade-users
More help   : https://help.launchpad.net/ListHelp