Re: [Yade-users] [Question #293955]: How to generate dense packing of aggregates

[VG]

Question #293955 on Yade changed:
https://answers.launchpad.net/yade/+question/293955

VG posted a new comment:
Hello Jan,

I have a follow-up question: I modified the compress.py script a little
bit to generate a periodic packing. Now, how do I use this compressed
packing in my main periodic simulation (MWE posted in #1 above) ? The
specific problems are:

1. How to get the size of periodic cell for main simulation, since the
compressed sample is going to be of smaller dimensions.

2. How to fit the sample exactly in the periodic cell of main
simulation, since the locations are going to be different ?

3. I want periodicity in x and z directions, but I am using top and
bottom boxes in the vertical direction to restrict periodicity in y
direction. I need to place these plates such that they precisely fit the
periodic cell in x & z direction, and contain the granular sample in
y-direction. Now, if I get the solution to above two, that should be
doable.


Just in case, here is the modified script of compress.py:

##
# Compress the loose packing into dense packing. Each agglomerate is
# considered as clump in this stage
##
from yade import export,ymport
import random
random.seed(1)

O.periodic = True

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

# load spheres from file, including information of their agglomerates
ids

sp = ymport.textExt('./cloud.txt')

for s in sp: O.bodies.append(s)

O.cell.setBox(dim)


O.engines = [
ForceResetter(),
InsertionSortCollider([Bo1_Sphere_Aabb(),Bo1_Box_Aabb()],
allowBiggerThanPeriod=True),
InteractionLoop(
[Ig2_Sphere_Sphere_ScGeom(),Ig2_Box_Sphere_ScGeom()],
[Ip2_FrictMat_FrictMat_FrictPhys()],
[Law2_ScGeom_FrictPhys_CundallStrack()]
),
PeriIsoCompressor(
charLen = 2.,
stresses=[-100e5,-1e4],
maxUnbalanced=1e-2,
doneHook= 'O.pause();',
globalUpdateInt=5,keepProportions=True),
NewtonIntegrator(damping=.6),
]
O.dt = PWaveTimeStep()

# compress the sample

O.run()
O.wait()


# save the result
export.textExt('./compressed.txt')

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



Thanks
Varun

-- 
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 #294232]: change in volume of meniscus as a function of capillary pressure

[Christian Jakob]

Question #294232 on Yade changed:
https://answers.launchpad.net/yade/+question/294232

Status: Open => Answered

Christian Jakob proposed the following answer:
I think you do not need i.phys.isBroken at all. What you need is at
least 1 step between parameter change and measurement.

Good luck

-- 
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 #294232]: change in volume of meniscus as a function of capillary pressure

[Seungcheol Yeom]

Question #294232 on Yade changed:
https://answers.launchpad.net/yade/+question/294232

Status: Answered => Open

Seungcheol Yeom is still having a problem:
Hi Chris,

I have tried that way but the volume of meniscus became 0 after 9900 of 
capillary pressure.
I think I need something to reactivate is.Broken = False to create the meniscus.
Thank you.

Seungcheol

-- 
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 #294232]: change in volume of meniscus as a function of capillary pressure

[Christian Jakob]

Question #294232 on Yade changed:
https://answers.launchpad.net/yade/+question/294232

Status: Open => Answered

Christian Jakob proposed the following answer:
Hi,

O.engines =
O.engines[:3]+[Law2_ScGeom_CapillaryPhys_Capillarity(capillaryPressure=cur_capillary,label='capillary')]+O.engines[3:]

Maybe this line is the problem in your script (not tested, not sure what
it does). If the engine label is 'capillary' you can try this line
instead:

capillary.capillaryPressure = cur_capillary

-- 
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 #294209]: Belt Conveyor create

[Anton Gladky]

Question #294209 on Yade changed:
https://answers.launchpad.net/yade/+question/294209

Anton Gladky proposed the following answer:
I do not quite understand, what you mean. Here [1] you go
through all facets and set the constant velocity. If it must
be a function of a time - it should not be a problem too.

[1]
https://github.com/yade/trunk/blob/master/examples/conveyor/conveyor.py#L32

Anton


2016-05-25 15:13 GMT+02:00 paula :
> Thanks Anton, but I need to give an specific velocity to the belt. Anyway 
> thanks again.
> Paula

-- 
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 #294232]: change in volume of meniscus as a function of capillary pressure

[Seungcheol Yeom]

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

Hello all,

I hope you all well. I have a simple question.
I am wondering whether I can simulate the change in volume of meniscus as a 
function of capillary pressure.
It seems like I can change the capillary pressure in the loop but the 
interaction between the spheres are not updated.
Here is my script:

from yade import plot, geom

r = 3e-6 #particle radius of clay, meter
h = 1e-6 #praticle distance, meter

#create two sphere paticles#
O.bodies.append([
   utils.sphere(center=(0,0,0),radius=r,fixed=True),
   utils.sphere((0,0,2*r+h),r,fixed=True)
])


#define engines#
O.engines=[
   ForceResetter(), #0
   InsertionSortCollider([Bo1_Sphere_Aabb(aabbEnlargeFactor=2)]), #1
   InteractionLoop(
  [Ig2_Sphere_Sphere_ScGeom(interactionDetectionFactor=2)], #2
  [Ip2_FrictMat_FrictMat_MindlinCapillaryPhys()], 
  [Law2_ScGeom_MindlinPhys_Mindlin(neverErase=True)] 
   ),
   
Law2_ScGeom_CapillaryPhys_Capillarity(capillaryPressure=1,label='capillary'),
 #3
   NewtonIntegrator(damping=0.4,gravity=(0,0,0)), #4
   PyRunner(command='change_capillary()',iterPeriod=100,label='start'), #5
]

O.trackEnergy = True

capillary.createDistantMeniscii=True
O.run(1,True)
capillary.createDistantMeniscii=False 

#define the force on the sphere#
def change_capillary():
global cur_capillary,delta_V,count
count = 0
while (count<1):
cur_capillary = capillary.capillaryPressure - 100
i = O.interactions[0,1]
delta_V = i.phys.vMeniscus
i.phys.isBroken = True
O.engines = 
O.engines[:3]+[Law2_ScGeom_CapillaryPhys_Capillarity(capillaryPressure=cur_capillary,label='capillary')]+O.engines[3:]
count = count + 1
print cur_capillary,delta_V
if cur_capillary < 9000:
O.pause()

O.dt=0.5*PWaveTimeStep()
O.run()


Any helps are appreciated.

Sincerely,

Seungcheol




-- 
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 #294209]: Belt Conveyor create

[paula]

Question #294209 on Yade changed:
https://answers.launchpad.net/yade/+question/294209

paula posted a new comment:
Thanks Anton, but I need to give an specific velocity to the belt. Anyway 
thanks again.
Paula

-- 
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 #294036]: Rate of strain between compressive and tensile strength using CpmMat.

[liukeqi]

Question #294036 on Yade changed:
https://answers.launchpad.net/yade/+question/294036

liukeqi posted a new comment:
Thanks Jan Stránský. I was very glad of your consistent help.

-- 
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 #294215]: how computer configuration to run discrete element method in YADE

[Christian Jakob]

Question #294215 on Yade changed:
https://answers.launchpad.net/yade/+question/294215

Status: Open => Answered

Christian Jakob proposed the following answer:
Please have a look to the installation instructions:

https://yade-dem.org/doc/installation.html#prerequisities

-- 
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 #294209]: Belt Conveyor create

[Anton Gladky]

Question #294209 on Yade changed:
https://answers.launchpad.net/yade/+question/294209

Status: Open => Answered

Anton Gladky proposed the following answer:
Hi,

you can try this example [1]. The conveyor will not move itself,
but it will give the particles the necessary velocity components.

[1] https://github.com/yade/trunk/tree/master/examples/conveyor

Regards

Anton

-- 
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