[Yade-users] [Question #665529]: about indentation error in example of yade website

2018-03-12 Thread Alireza Sadeghi
New question #665529 on Yade:
https://answers.launchpad.net/yade/+question/665529

Hello,

I have a question about indentation error in loading example number two 
(Gravity deposition) of https://yade-dem.org/doc/tutorial-examples.html
When I run it, it say that unbalancedForce() is not defined!
why does this error appear?how can I fix it?
Does yade know how to define unbalancedForce() or I have to define for it?

Thank you very much for your help in advance.

Best Regards

Alireza

-- 
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 #665514]: Application of flow engines on saturated flows

2018-03-12 Thread Robert Caulk
Question #665514 on Yade changed:
https://answers.launchpad.net/yade/+question/665514

Robert Caulk posted a new comment:
Hello Nicolas,

I can only speak for the standard FlowEngine (and I think TwoPhaseFlow
derives from FlowEngine). It is a finite volume scheme, so it requires
the application and maintenance of boundary conditions. As of right now,
the FlowEngine will automatically apply and maintain these for you in a
box configuration. Your described requirements sound like a highly
dynamic and arbitrarily shaped system - FlowEngine will likely need
heavy modifications to produce the results you are after.

In response to your question "can I apply flowengine only inside
granular media": yes, FlowEngine triangulates the pores by using the
particles as nodes.

I have never used HydrodynamicsLBM so I am not sure if it is suitable or
not for your simulation.

Best,

Robert

-- 
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 #665387]: Compile error JointedCohesiveFrictionalPM

2018-03-12 Thread Thomas Chauve
Question #665387 on Yade changed:
https://answers.launchpad.net/yade/+question/665387

Thomas Chauve proposed the following answer:
Hello

I have similar probleme with the latest mater branch (from 12/03/2018). I 
compile it on Debian 9 with GCC 6.3.0.
I had no trouble before (last version tested is the one dated of 21/02/2018).

I see that JointedCohesiveFrictionalPM.cpp JointedCohesiveFrictionalPM.cpp have 
been modified 6 day ago. 
I try someting : I compile the new master (from 12/03/2018) with the old 
JointedCohesiveFrictionalPM.cpp JointedCohesiveFrictionalPM.cpp (from 
21/02/2018 branch) and it works.

I hope it can helps.

Cheers
Thomas

-- 
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 #665440]: how can make a path for yade to load python files

2018-03-12 Thread Alireza Sadeghi
Question #665440 on Yade changed:
https://answers.launchpad.net/yade/+question/665440

Status: Answered => Solved

Alireza Sadeghi confirmed that the question is solved:
Dear Jan,

Thank you very much for your help. It was about the folder that my
program was in it. it was solved. Thanks a lot.

Best Regards

Alireza

-- 
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 #665440]: how can make a path for yade to load python files

2018-03-12 Thread Alireza Sadeghi
Question #665440 on Yade changed:
https://answers.launchpad.net/yade/+question/665440

Alireza Sadeghi confirmed that the question is solved:
Thanks Jan Stránský, that solved my question.

-- 
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 #665350]: The parameter (en, es )of the DMT model

2018-03-12 Thread zhao dejin
Question #665350 on Yade changed:
https://answers.launchpad.net/yade/+question/665350

Status: Answered => Solved

zhao dejin confirmed that the question is solved:
Thanks Chiara Modenese, that solved my question.

-- 
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 #665514]: Application of flow engines on saturated flows

2018-03-12 Thread Nicolas Godet
New question #665514 on Yade:
https://answers.launchpad.net/yade/+question/665514

Dear all,

I am a new user of YADE. The purpose of my internship is to find a way to model 
the flow of a saturated material down a chute [1].

I read the documentation about flow engines available in YADE such as 
FlowEngine() [2] [3] and HydrodynamicsLawLBM() [4]. Moreover, I am still 
looking for documentation for TwoPhaseFlowEngine() [5].

I conclude that all engines need to be defined in a box and that the engines 
can’t be used for modelling a flow along a chute.

 

I will explain how I plan to proceed.
I would like to realize a release of a saturated volume of a mix of water 
(density of 1200kg/m³ and higher viscosity than water viscosity) and granular 
media. This volume will flow along a rectangular channel which extremity is 
closed by a wall.

First step, I do a gravity deposit at the top of my channel. A wall holds this 
volume. Once it is stabilized (unbalancedForce() < threshold), I remove the 
wall and let the volume flows along the channel, in a dam-break manner.

Regarded as I plan to do, I am not sure if I can apply a flow engine only 
inside the granular media.

Do flow engines necessarily have to be applied on a box ?
Maybe I miss understood something.

 

Thanks in advance,

Regards,

Nicolas

 

[1] https://en.wikipedia.org/wiki/Debris_flow

[2] Chareyre, B., Cortis, A., Catalano, E., Barthélemy, E., 2012. Pore-Scale 
Modeling of Viscous Flow and Induced Forces in Dense Sphere Packings. Transp. 
Porous Media 92, 473‑493. https://doi.org/10.1007/s11242-011-9915-6

[3] Catalano, E., Chareyre, B., Barthélémy, E., 2014. Pore-scale modeling of 
fluid-particles interaction and emerging poromechanical effects. Int. J. Numer. 
Anal. Methods Geomech. 38, 51‑71. https://doi.org/10.1002/nag.2198

[4] Lominé, F., Scholtès, L., Sibille, L., Poullain, P., 2013. Modeling of 
fluid‑solid interaction in granular media with coupled lattice 
Boltzmann/discrete element methods: application to piping erosion. Int. J. 
Numer. Anal. Methods Geomech. 37, 577‑596. https://doi.org/10.1002/nag.1109

[5] https://yade-dem.org/doc/yade.wrapper.html#yade.wrapper.TwoPhaseFlowEngine

-- 
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 #665350]: The parameter (en, es )of the DMT model

2018-03-12 Thread Chiara Modenese
Question #665350 on Yade changed:
https://answers.launchpad.net/yade/+question/665350

Status: Open => Answered

Chiara Modenese proposed the following answer:
Hi,

Glad to know the work has been useful.

Surface energy (gamma) is a measurable material parameter, and also an
input of the DMT model. You need to find out what the typical value(s)
is for your material based on the relevant literature (sorry but I can't
help you directly with this).

For the damping coefficients, en and es, these are there to replicate
dissipation forces, which you will need to calibrate in your model
against experimental values, depending on the purpose of your DEM
simulations.

I hope this is of help,
Chiara

-- 
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 #665350]: The parameter (en, es )of the DMT model

2018-03-12 Thread Chiara Modenese
Question #665350 on Yade changed:
https://answers.launchpad.net/yade/+question/665350

Status: Open => Answered

Chiara Modenese proposed the following answer:
Hi,

If you are looking to use the DMT model for adhesive particles, this is
a simple extension of the Hertz Mindlin formulation. The Hertzian
contact relationship is shifted down by the adhesive force. You can
check this on the literature, or at page 39, 40 of my PhD thesis [*].

The adhesive force should already be implemented in the current version
of the HertzMindlin.cpp contact law and it is a function of gamma.
Apologies if this has not been documented properly but it should work
fine if you are familiar with Yade. If you need further assistance, let
me know and I can provide you with a sample script.

HTH,
Chiara
[*] http://www.eng.ox.ac.uk/civil/publications/theses/modenese_pdf

-- 
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 #665350]: The parameter (en, es )of the DMT model

2018-03-12 Thread zhao dejin
Question #665350 on Yade changed:
https://answers.launchpad.net/yade/+question/665350

Status: Answered => Open

zhao dejin is still having a problem:
Hi,
Thanks ,Chiara, I have read your PhD thesis yet ,It's really a wonderful work 
and thanks to your work ,I am  familiar with  the  DMT model .  Three parameter 
are  need in the DMT model of  YADE, such as gamma ,en ,es.  But I failed to 
find the effective valus  of  the 316L stainless steel which is the material I 
have studied .So I turn to someone did the research of this material for 
help.Thanks again.

Best wishes
ZDJ

-- 
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 #665350]: The parameter (en, es )of the DMT model

2018-03-12 Thread zhao dejin
Question #665350 on Yade changed:
https://answers.launchpad.net/yade/+question/665350

Summary changed to:
The parameter (en, es )of the DMT model 

Description changed to:
Hi , yade users :
  recently, I'm doing the DEM simulation  of the stainless steel, the DMT 
model must be considered. However , I fail to find the effective  parameters of 
the material, such as gamma, en, es. If someone did the research of this 
material or somebody can provide me some paper about it . Thanks in advance.
 Best wishes 
ZDJ

-- 
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 #665440]: how can make a path for yade to load python files

2018-03-12 Thread Jan Stránský
Question #665440 on Yade changed:
https://answers.launchpad.net/yade/+question/665440

Status: Open => Answered

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

> I am new user of YADE

welcome :-)

the error says, that Python does not know about Hello.py file. There are a few 
possible reasons
1) the actual file could be named differently, hello.py (linux is case 
sensitive)
2) in the terminal, you are in different location than you saved Hello.py file. 
So e.g. if you save your hallo.py in e.g. Documents/yade, you should write in 
terminal
cd Documents/yade
and then yade Hello.py

let us know if it helped

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 #665507]: The connection between grid nodes disappear

2018-03-12 Thread Alina Irsainova
New question #665507 on Yade:
https://answers.launchpad.net/yade/+question/665507

Hello,
I am trying to write a script to place a geogrid in soil packing. My script is 
based on the example : 
trunk/examples/grids/Simple_Grid_Falling.py
The problem is in the combination of the script modelling soil packing and the 
script for geogrid. After combining them, the connections between the gridNodes 
disappear. 
Anyone can help me to find the bugs in this script?

This is my script:

from yade import pack,geom,qt
from yade.gridpfacet import *

rad,gap=.25,.005

#Parameters, which will be passed into spheres and facets creators
kw={'material':0}
kwBoxes={'color':[1,0,0],'wire':False,'dynamic':True,'material':0}
kwMeshes={'color':[1,0,0],'wire':True,'dynamic':False,'material':0}

O.materials.append(FrictMat(young=10e9,poisson=.25,frictionAngle=0.5,density=1e3,label='sphere'))
O.materials.append(CohFrictMat(young=3e2,poisson=0.3,density=1e1,frictionAngle=10,normalCohesion=1e7,shearCohesion=1e7,momentRotationLaw=True,label='spheremat'))


O.engines=[
ForceResetter(),
InsertionSortCollider([
Bo1_GridConnection_Aabb(),Bo1_Sphere_Aabb(),Bo1_Facet_Aabb(),
]),
InteractionLoop(

[Ig2_GridNode_GridNode_GridNodeGeom6D(),Ig2_Sphere_Sphere_ScGeom(),Ig2_Facet_Sphere_ScGeom(),Ig2_GridNode_Sphere_GridNodeGeom6D(),Ig2_GridNode_Facet_GridNodeGeom6D(),],

[Ip2_CohFrictMat_CohFrictMat_CohFrictPhys(setCohesionNow=True,setCohesionOnNewContacts=False),Ip2_FrictMat_FrictMat_FrictPhys()],

[Law2_ScGeom6D_CohFrictPhys_CohesionMoment(),Law2_ScGeom_FrictPhys_CundallStrack()]
),
NewtonIntegrator(gravity=(0,0,-1000),damping=0.3,label='newton')
]
 
### Parameters of a rectangular grid ###
L=11 #length [m]
l=11#width  [m]
nbL=11  #number of nodes for the length [#]
nbl=11  #number of nodes for the width  [#]
r=L/100.#radius
color=[255./255.,102./255.,0./255.]
nodesIds=[]
#Create all nodes first :
for i in range(0,nbL):
for j in range(0,nbl):
nodesIds.append( 
O.bodies.append(gridNode([i*L/nbL-20,j*l/nbl,-2.7],r,wire=False,fixed=True,material='spheremat',color=color))
 )

#Create connection between the nodes
for i in range(0,len(nodesIds)):
for j in range(i+1,len(nodesIds)):
dist=(O.bodies[i].state.pos - O.bodies[j].state.pos).norm()
if(dist<=L/nbL*1.01):
O.bodies.append( gridConnection(i,j,r,color=color) )

#Create packing
O.bodies.append(pack.regularHexa(pack.inParallelepiped((-10,0,-2.75),(-20,0,-2.75),(-10,10,-2.75),(-10,0,-7.75)),radius=rad,gap=rad/3.0,color=(0.5,0.5,0.1),material='sphere'))
O.bodies.append(pack.regularHexa(pack.inParallelepiped((-10,0,-2.5),(-20,0,-2.5),(-10,10,-2.5),(-10,0,2.25)),radius=rad,gap=rad/3.0,color=(0.5,0.5,0.1),material='sphere'))
O.bodies.append(geom.facetBox((-15,5,-5),(5,5,2.5),wallMask=15,**kwMeshes))
O.bodies.append(geom.facetBox((-15,5,0),(5,5,2.5),wallMask=15,**kwMeshes))
O.bodies.append(geom.facetBox((-15,5,-5),(5,5,2.5),wallMask=16,**kwBoxes))
O.bodies.append(geom.facetBox((-15,5,2.5),(5,5,0),wallMask=16,**kwBoxes))

#Set a fixed node
O.bodies[0].dynamic=False

O.dt=1e-05
O.saveTmp()
O.timingEnabled=True
O.run(1,True)
from yade import timing
timing.stats()
qt.View()


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