Re: [Yade-users] [Question #700813]: Batch run and run alone get different results

2022-03-07 Thread Ziyu Wang
Question #700813 on Yade changed:
https://answers.launchpad.net/yade/+question/700813

Status: Needs information => Open

Ziyu Wang gave more information on the question:
Hi Jan
>Try without any -j or --job-threads (i.e. using one core for both cases)
As you said, I have run yadedaily-batch twice and yadedaily separate twice 
using one core(with the same parameters).However,the results obtained are 
consistent with those before, and there is still a large gap between the two.

>just to be sure, do you use same young parameter for both yadedaily and 
>yadedaily-batch?
Yes,I use the same young(3000e9) for both.And all the parameters are same.The 
results(stress and strength) I mentioned are also obtained from the simulation 
of the same parameters.

Thanks for 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 #700761]: about measure seepage discharge in triaxial test

2022-03-07 Thread Ziyu Wang
Question #700761 on Yade changed:
https://answers.launchpad.net/yade/+question/700761

Ziyu Wang gave more information on the question:
Hello,everyone.

I tried changing flow.permeabilityFactor to a positive value and other 
values.However,It seems that the data I record is always constant( 
flow.getBoundaryFlux(4) and (5) ).
In my understanding, with the increase of axial pressure in the process of 
triaxial compression, the permeability will change, which should not be a 
constant value.
How to deal with it?

Thanks.

-- 
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 #700793]: About the "WallStresses_triax_base_" file in script-session1.py

2022-03-07 Thread CESAR FELIPE JIMENEZ SANTIAGO
Question #700793 on Yade changed:
https://answers.launchpad.net/yade/+question/700793

Status: Open => Solved

CESAR FELIPE JIMENEZ SANTIAGO confirmed that the question is solved:
Hi Jerome,

Thank you for your answer. It's solved

-- 
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 #700847]: access friction and inelastic collision contributions in the total loss of energy of the system

2022-03-07 Thread Jan Stránský
Question #700847 on Yade changed:
https://answers.launchpad.net/yade/+question/700847

Status: Open => Answered

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

> Is there a way to access the friction energy dissipation contribution?

see [1,2]

> and the impact energy dissipation?

You have (possibly constant) total energy Etot.
You can compute kinetic energy Ekin [3], elastic energy Eel [4], energy 
dissipated by friction [1,2]
The rest is other dissipation Eother.
Etot = Eel + Ekin + Efrict + Eother

Or you can use O.energy [5,6,7], where you can get other values if Law2 
implements it.
O.energy also stores artificial dissipation of newton integrator numerical 
damping.

You can compare both approaches :-)

> loss by inelastic collisions

Note that for your case (only Law2_ScGeom_FrictPhys_CundallStrack),
theoretically the only dissipation is only by friction and the numerical
damping of newton integrator (there are no true "inelastic collisions")

Cheers
Jan

[1] 
https://yade-dem.org/doc/yade.wrapper.html#yade.wrapper.Law2_ScGeom_FrictPhys_CundallStrack.plasticDissipation
[2] 
https://yade-dem.org/doc/yade.wrapper.html#yade.wrapper.Law2_ScGeom_FrictPhys_CundallStrack.traceEnergy
[3] https://yade-dem.org/doc/yade.utils.html#yade._utils.kineticEnergy
[4] 
https://yade-dem.org/doc/yade.wrapper.html#yade.wrapper.Law2_ScGeom_FrictPhys_CundallStrack.elasticEnergy
[5] https://yade-dem.org/doc/yade.wrapper.html#yade.wrapper.Omega.energy
[6] https://yade-dem.org/doc/yade.wrapper.html#yade.wrapper.Omega.trackEnergy
[7] https://yade-dem.org/doc/yade.wrapper.html#yade.wrapper.EnergyTracker

-- 
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 #700788]: performing wetting condition on Twophaseflow engine

2022-03-07 Thread Hossein
Question #700788 on Yade changed:
https://answers.launchpad.net/yade/+question/700788

Status: Open => Solved

Hossein confirmed that the question is solved:
Problem Solved

-- 
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 #700843]: Difference between the internal friction angle of material and the interparticle friction angle?

2022-03-07 Thread Chien-Cheng Hung
Question #700843 on Yade changed:
https://answers.launchpad.net/yade/+question/700843

Status: Answered => Solved

Chien-Cheng Hung confirmed that the question is solved:
Thanks Jérôme Duriez, 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 #700843]: Difference between the internal friction angle of material and the interparticle friction angle?

2022-03-07 Thread Jérôme Duriez
Question #700843 on Yade changed:
https://answers.launchpad.net/yade/+question/700843

Status: Open => Answered

Jérôme Duriez proposed the following answer:
Above [1] and [2] point to a contact property between 2 Discrete
Elements, that's it.

If you take any DEM paper from the literature with some macro-scale
shear strength considerations (e.g. triaxial testing), and compare the
deduced material friction angle (e.g. Mohr-Coulomb criterion friction
angle) with that contact property you will see how different both can be
(except in few instances where you would be very lucky -- or unlucky --
and get an equality between the two).

See e.g.
https://www.sciencedirect.com/science/article/pii/S0266352X21005309
where Mohr Coulomb friction angle is above 35 degrees, vs a contact
friction angle at 31° for a polyhedron model ("Set 2 parameters" in
Table 5).

Most convincing example would certainly be your own YADE simulations of
triaxial tests, with zero FrictMat.frictionAngle.

-- 
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 #700830]: Is there any way to generate voronoi tessellation diagram?

2022-03-07 Thread Bruno Chareyre
Question #700830 on Yade changed:
https://answers.launchpad.net/yade/+question/700830

Status: Needs information => Answered

Bruno Chareyre proposed the following answer:
Hi,
The image [1] was generated semi-automatically with POV-ray, using data 
generated collected from FlowEngine.
There is no builtin function for doing that in Yade directly.
It seems your question has two sides:
1/ how to get the data
2/ how to display it

For 1/, it's pretty easy in C++ to iterate on the triangulation's facets
and output a dual segment for each (one branch in Voronoi), but there is
currently no python wrapping for such operation in TesselationWrapper.
You may need to code a little within TessealationWrapper, or switch to
FlowEngine. I can advise a bit more if needed.

For 2/, it is up to you. I'm not quite sure POV-ray still works on
modern systems. There could be other softwares doing 3D rendering.
There's paraview but I'm not sure precisely how to achieve [1] in
paraview. Blender is another one but I have no experience.

Bruno

-- 
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 #700843]: Difference between the internal friction angle of material and the interparticle friction angle?

2022-03-07 Thread Chien-Cheng Hung
Question #700843 on Yade changed:
https://answers.launchpad.net/yade/+question/700843

Status: Needs information => Open

Chien-Cheng Hung gave more information on the question:
Hi Jan,

Sorry, let me make it clear.

As far as I know, the friction angle we assign or modify in [1] and [2] stands 
for the internal friction angle of the material.
In the meanwhile, these two parameters also represent the interparticle 
friction angle that describes the contact friction between particles.

I happened to see Luc's comment in [3] saying not to mix up the internal
friction angle of material and the interparticle friction angle which I
think are the same thing based on my understanding.

So I would like to ask whether the values in [1] and [2] can also be
considered as the interparticle friction coefficient?

Cheers,
Chien-Cheng

[1] 
https://yade-dem.org/doc/yade.wrapper.html?highlight=frictionangle#yade.wrapper.FrictMat.frictionAngle
[2] 
https://yade-dem.org/doc/yade.utils.html?highlight=setcontactfriction#yade._utils.setContactFriction
[3] https://answers.launchpad.net/yade/+question/678523

-- 
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 #700611]: TwoPhaseFlow - double free or corruption

2022-03-07 Thread Bruno Chareyre
Question #700611 on Yade changed:
https://answers.launchpad.net/yade/+question/700611

Status: Open => Needs information

Bruno Chareyre requested more information:
Oh... but I thought the question was about "large number of particles".
Now it's about "packs".

what do these packs look like? Would they contain overlapping spheres?
B

-- 
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 #700847]: access friction and inelastic collision contributions in the total loss of energy of the system

2022-03-07 Thread rox212
New question #700847 on Yade:
https://answers.launchpad.net/yade/+question/700847

Dear all,

Is there a way to access the friction energy dissipation contribution? and the 
impact energy dissipation?  Or to put in other words, the energy that is loss 
by friction and the energy that is loss by inelastic collisions?

I am using yade 2018b. And I am using the following engines with 
integrators/laws:

O.engines = [
  ForceResetter(),
  InsertionSortCollider([Bo1_Sphere_Aabb(),Bo1_Facet_Aabb()],label='collider'),
  InteractionLoop(
[Ig2_Sphere_Sphere_ScGeom(),Ig2_Facet_Sphere_ScGeom()],
[Ip2_FrictMat_FrictMat_FrictPhys(label='frict')],
[Law2_ScGeom_FrictPhys_CundallStrack(traceEnergy=True,label='cundall')],
  ),
GravityEngine(gravity=(0,0,-9.8)),
NewtonIntegrator(damping=dampingCoeffCundall,label='newtonInt'),

TranslationEngine(translationAxis=[0,0,1],velocity=velocityDamp,ids=idTop,dead=False,label='translat'),

TranslationEngine(translationAxis=[0,0,1],velocity=velocityDamp,ids=idWalls,dead=False,label='translatbottom'),

TranslationEngine(translationAxis=[0,0,1],velocity=velocityDamp,ids=idBottom,dead=False,label='translatBottom'),
HarmonicMotionEngine(A=amplitude, f=frequency , fi = initialPhase, ids = 
idTop, label='combEngine',dead=True),
HarmonicMotionEngine(A=amplitude, f=frequency , fi = initialPhase, ids = 
idWalls, label='combEngine1',dead=True),
HarmonicMotionEngine(A=amplitude, f=frequency , fi = initialPhase, ids = 
idBottom, label='combEngine2',dead=True),
]

Thank you in advance!





-- 
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 #700843]: Difference between the internal friction angle of material and the interparticle friction angle?

2022-03-07 Thread Jan Stránský
Question #700843 on Yade changed:
https://answers.launchpad.net/yade/+question/700843

Status: Open => Needs information

Jan Stránský requested more information:
Hello,

please provide more information [2].

What is "internal friction angle of material"? Yade material value? Yade
overall macroscopic value? Parameter of a real material? ... ?

What is "interparticle friction angle"? Some Yade value? Or a value of
real particles? Or ... ?

Cheers
Jan

[2] https://www.yade-dem.org/wiki/Howtoask

-- 
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 #700815]: access individual energy contribution from each particle

2022-03-07 Thread rox212
Question #700815 on Yade changed:
https://answers.launchpad.net/yade/+question/700815

Status: Open => Solved

rox212 confirmed that the question is solved:
Thanks Jérôme Duriez, 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 #678523]: frictionAngle vs μ(interparticle friction coefficient)

2022-03-07 Thread Chien-Cheng Hung
Question #678523 on Yade changed:
https://answers.launchpad.net/yade/+question/678523

Chien-Cheng Hung posted a new comment:
Okay, I will open a new question.
Thanks!

Chien-Cheng

-- 
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 #700843]: Difference between the internal friction angle of material and the interparticle friction angle?

2022-03-07 Thread Chien-Cheng Hung
New question #700843 on Yade:
https://answers.launchpad.net/yade/+question/700843

Hi all,

According to the comments in [1], Luc mentioned not to mix up the internal 
friction angle of a material and the interparticle friction angle.

I was wondering what the differences are between these two parameters.
Looking forward to getting your answer.

[1] https://answers.launchpad.net/yade/+question/678523

Cheers,
Chien-Cheng

-- 
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 #700815]: access individual energy contribution from each particle

2022-03-07 Thread rox212
Question #700815 on Yade changed:
https://answers.launchpad.net/yade/+question/700815

Status: Answered => Open

rox212 is still having a problem:
Dear Jan,

Thank you for your answer! Hope this is enough:

Yade version: 2018.02b

My O.engines are:

O.engines = [
  ForceResetter(),
  InsertionSortCollider([Bo1_Sphere_Aabb(),Bo1_Facet_Aabb()],label='collider'),
  InteractionLoop(
[Ig2_Sphere_Sphere_ScGeom(),Ig2_Facet_Sphere_ScGeom()],
[Ip2_FrictMat_FrictMat_FrictPhys()],
[Law2_ScGeom_FrictPhys_CundallStrack()],
  ),
GravityEngine(gravity=(0,0,-9.8)),
NewtonIntegrator(damping=dampingCoeffCundall,label='newtonInt'),

TranslationEngine(translationAxis=[0,0,1],velocity=velocityDamp,ids=idTop,dead=False,label='translat'),

TranslationEngine(translationAxis=[0,0,1],velocity=velocityDamp,ids=idWalls,dead=False,label='translatbottom'),

TranslationEngine(translationAxis=[0,0,1],velocity=velocityDamp,ids=idBottom,dead=False,label='translatBottom'),
HarmonicMotionEngine(A=amplitude, f=frequency , fi = initialPhase, ids = 
idTop, label='combEngine',dead=True),
HarmonicMotionEngine(A=amplitude, f=frequency , fi = initialPhase, ids = 
idWalls, label='combEngine1',dead=True),
HarmonicMotionEngine(A=amplitude, f=frequency , fi = initialPhase, ids = 
idBottom, label='combEngine2',dead=True),
]

> In principle it should/could be possible by storing interaction and
body values and comparing them between individual iterations.

Is there an example to do this?

Thank you for your time and 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 #700815]: access individual energy contribution from each particle

2022-03-07 Thread Jérôme Duriez
Question #700815 on Yade changed:
https://answers.launchpad.net/yade/+question/700815

Status: Needs information => Answered

Jérôme Duriez proposed the following answer:
Hi,

Related recent threads:

https://answers.launchpad.net/yade/+question/700748 and
https://answers.launchpad.net/yade/+question/700181

I a nutshell, you will most probably have to do it by yourself, this is
not a functionality already present in YADE. Good news is that you can
do by yourself anything you want in YADE

-- 
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 #700830]: Is there any way to generate voronoi tessellation diagram?

2022-03-07 Thread Jan Stránský
Question #700830 on Yade changed:
https://answers.launchpad.net/yade/+question/700830

Status: Open => Needs information

Jan Stránský requested more information:
Hello,

> I'd like to generate voronoi tessellation diagram

Do you need nubers, of picture?

From what input would you like to generate it? Particle centers?

> I tried the TW.defToVtk() and got Delaunay triangulation diagram but
not voronoi tesselation diagram.

They are dual, one uniquely determines the other, so if you have
Delaunay triangulation, creating Voronoi tesselation is "easy".

> It seems the voro++[2] is what I am looking for, but the result seems
not good because I didn't see any spheres.

I don't know voro++, but if I remember correctly, Paraview could do this.
If not natively, you can always create the Voronoi diagram and "manually" 
display the lines together with spheres.

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 #700815]: access individual energy contribution from each particle

2022-03-07 Thread Jan Stránský
Question #700815 on Yade changed:
https://answers.launchpad.net/yade/+question/700815

Status: Open => Needs information

Jan Stránský requested more information:
Hello,

please read [1] and provide a MWE or at least what Ip2/Law2 you are using.
The answer strongly depends on it.

In principle it should/could be possible by storing interaction and body
values and comparing them between individual iterations.

Cheers
Jan

[1] https://www.yade-dem.org/wiki/Howtoask

-- 
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 #700800]: simulate rectangular lastic

2022-03-07 Thread Jan Stránský
Question #700800 on Yade changed:
https://answers.launchpad.net/yade/+question/700800

Status: Open => Needs information

Jan Stránský requested more information:
Hello,

> I am beginner in yade

welcome :-)

> I want to simulate rectangle slices rubber and circle sand but l don't
know how to simulate rectangular rubber in yade .

please read [1] and provide (much much) more information.
For example:
- what is your background, where do/did you study, where, do you work ... ?
- what is this task about? short-term industrial project? Topic for full PhD 
thesis? ... ?
- why you want to use Yade (or DEM in general). Rubber as is sounds more 
suitable for FEM or similar methods.
- what about rubber you want to simulate. manufacturing? elastic bending? 
tension with cracking? . ?
- what are general dimensions of the simulation (specimen size, time of the 
simulated event)
- what is "circle sand"? It is very ambiguous, I can imagine more very 
different interpretations (circular particles, the domain is circular, ...)
- what is "rectangle slices rubber"? (same as above)
- ...

Cheers
Jan

[1] https://www.yade-dem.org/wiki/Howtoask

-- 
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 #700813]: Batch run and run alone get different results

2022-03-07 Thread Jan Stránský
Question #700813 on Yade changed:
https://answers.launchpad.net/yade/+question/700813

Status: Open => Needs information

Jan Stránský requested more information:
Try without any -j or --job-threads (i.e. using one core for both cases)

> I have run the same parameters twice using batch mode

just to be sure, do you use same young parameter for both yadedaily and 
yadedaily-batch?
It is not mentioned, only mentioned is content of young.txt, were are various 
values..

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