Re: [Yade-users] [Question #263634]: MatchMaker use

2015-03-16 Thread Anton Gladky
Question #263634 on Yade changed:
https://answers.launchpad.net/yade/+question/263634

Anton Gladky proposed the following answer:
It works also with Ip2_FrictMat_FrictMat_FrictPhys::frictAngle

==
#!/usr/bin/python
# -*- coding: utf-8 -*-


from yade import pack,timing

r1 = 0.002381
r2 = 0.002381

mat1 = 
O.materials.append(FrictMat(young=50e9,density=6000,poisson=.2,label='shots1'))
mat2 = 
O.materials.append(FrictMat(young=50e9,density=6000,poisson=.2,label='shots2'))
mat3 = 
O.materials.append(FrictMat(young=50e9,density=6000,poisson=.2,label='shots3'))


id11 = 
O.bodies.append(sphere(center=[0,0,0],radius=r1,material=mat1,fixed=True,color=[0,0,1]))
id12 = 
O.bodies.append(sphere(center=[0,0,(r1+r2)],radius=r2,material=mat2,fixed=False,color=[0,0,1]))

id21 = 
O.bodies.append(sphere(center=[3*r1,0,0],radius=r1,material=mat1,fixed=True,color=[0,1,0]))
id22 = 
O.bodies.append(sphere(center=[3*r1,0,(r1+r2)],radius=r2,material=mat3,fixed=False,color=[0,1,0]))

id31 = 
O.bodies.append(sphere(center=[6*r1,0,0],radius=r1,material=mat2,fixed=True,color=[1,0,0]))
id32 = 
O.bodies.append(sphere(center=[6*r1,0,(r1+r2)],radius=r2,material=mat3,fixed=False,color=[1,0,0]))


O.engines=[
ForceResetter(),
InsertionSortCollider([Bo1_Sphere_Aabb(),Bo1_Facet_Aabb(),]),
InteractionLoop(
[Ig2_Sphere_Sphere_ScGeom(),Ig2_Facet_Sphere_ScGeom()],
[Ip2_FrictMat_FrictMat_FrictPhys(

frictAngle=MatchMaker(matches=((mat1,mat2,.1),(mat1,mat3,.2),(mat2,mat3,.3))),
label='damping')],
[Law2_ScGeom_FrictPhys_CundallStrack(label='contactLaw')]
),
NewtonIntegrator(damping=.3,gravity=(0,0,-1e3)), # gravity artificially 
high, to make it faster going ;-)
]
O.dt=PWaveTimeStep()

O.step()

print Friction coefficient for id11 and id12 is 
%g%(math.atan(O.interactions[id11,id12].phys.tangensOfFrictionAngle))
print Friction coefficient for id21 and id22 is 
%g%(math.atan(O.interactions[id21,id22].phys.tangensOfFrictionAngle))
print Friction coefficient for id31 and id32 is 
%g%(math.atan(O.interactions[id31,id32].phys.tangensOfFrictionAngle))


from yade import qt
qt.Controller()
qt.View()
=

-- 
You received this question notification because you are a member of
yade-users, which 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 #263634]: MatchMaker use

2015-03-16 Thread Andy Torres
Question #263634 on Yade changed:
https://answers.launchpad.net/yade/+question/263634

Andy Torres posted a new comment:
Hi Bruno and Anton,
@Bruno, (in answer to #3) I have been using the same script without the 
MatchMaker, changing the values of frictAngle without any trouble, and try 
changing the material id, etc. so i think the problem is not the script, if 
needed, i can post the whole script.
@Anton, thanks, i will test those modifications and let you know how it goes!

-- 
You received this question notification because you are a member of
yade-users, which 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 #263634]: MatchMaker use

2015-03-16 Thread Anton Gladky
Question #263634 on Yade changed:
https://answers.launchpad.net/yade/+question/263634

Anton Gladky proposed the following answer:
in viscoelastic model MatchMaker is implemented only for en, et and tc [1]. For
friction coefficient the minimal for 2 interacting materials is taken [2].

I think, we can implement MatchMaker if needed not causing  regressions.


[1] https://github.com/yade/trunk/blob/master/pkg/dem/ViscoelasticPM.cpp#L188
[2] https://github.com/yade/trunk/blob/master/pkg/dem/ViscoelasticPM.cpp#L247

-- 
You received this question notification because you are a member of
yade-users, which 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 #263634]: MatchMaker use

2015-03-16 Thread Anton Gladky
Question #263634 on Yade changed:
https://answers.launchpad.net/yade/+question/263634

Anton Gladky proposed the following answer:
Should be fixed.

https://github.com/yade/trunk/commit/cdac2f989e1ef8f0c27919b904cecc8186b22f99

-- 
You received this question notification because you are a member of
yade-users, which 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 #263634]: MatchMaker use

2015-03-16 Thread Bruno Chareyre
Question #263634 on Yade changed:
https://answers.launchpad.net/yade/+question/263634

Status: Open = Answered

Bruno Chareyre proposed the following answer:
Hello Andy, what you describe sounds like a bug. MatchMakers would be broken?
Did you check some example scripts using matchmakers to determine if the 
problem is in Yade or in your own script (wrong material Id, this kind of 
thing...)?

-- 
You received this question notification because you are a member of
yade-users, which 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 #263634]: MatchMaker use

2015-03-16 Thread Bruno Chareyre
Question #263634 on Yade changed:
https://answers.launchpad.net/yade/+question/263634

Bruno Chareyre posted a new comment:
@Anton,
Thanks but I think the OP is still a pending issue.
Ip2_FrictMat_FrictMat_FrictPhys::frictAngle is already a MatchMaker, so the 
first attempt of Andy should have work in the first place.

-- 
You received this question notification because you are a member of
yade-users, which 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 #263739]: isotropic compression at controled pressure

2015-03-16 Thread Luc Sibille
Question #263739 on Yade changed:
https://answers.launchpad.net/yade/+question/263739

Status: Open = Answered

Luc Sibille proposed the following answer:
Hello Luc,

Yes it is possible.
You can find a detailed example about how to use TriaxialStresscontroller here:
http://bazaar.launchpad.net/~yade-pkg/yade/git-trunk/view/head:/examples/triax-tutorial/script-session1.py

to control the pressure on all the walls you have do define the arguments:
stressMask = 7 
and
goal1=goal2=goal3=the pressure that you want to (with a minus sign for a 
compressive pressure).

Cheers,
Luc

 
Le Lundi 16 Mars 2015 17:16 CET, Luc OGER 
question263...@answers.launchpad.net a écrit: 
 
 New question #263739 on Yade:
 https://answers.launchpad.net/yade/+question/263739
 
 I am new in the use of Yade and it is not clear for me how to define some 
 parameters for the TriaxialStresscontroller class in order to  make a 
 triaxial compression with identical pressures in all the walls.
 In my understanding of the class I can controled the strain rate and/or the 
 wall displacements but not the pressure on all the walls to have the same 
 value during the isotropic compression.
 
 so my question :
 is it possible ? and if yes how I can do that?
 
 thanks in advance.
 
 
 -- 
 You received this question notification because you are a member of
 yade-users, which 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
 
 
-- 
Luc Sibille
Université Joseph Fourier / IUT1 de Grenoble
Laboratoire 3SR UMR CNRS

Tel lab.: +33 (0)4 76 82 63 48
Tel IUT: +33 (0)4 76 82 53 36

You received this question notification because you are a member of
yade-users, which 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 #263739]: isotropic compression at controled pressure

2015-03-16 Thread Luc OGER
New question #263739 on Yade:
https://answers.launchpad.net/yade/+question/263739

I am new in the use of Yade and it is not clear for me how to define some 
parameters for the TriaxialStresscontroller class in order to  make a triaxial 
compression with identical pressures in all the walls.
In my understanding of the class I can controled the strain rate and/or the 
wall displacements but not the pressure on all the walls to have the same value 
during the isotropic compression.

so my question :
is it possible ? and if yes how I can do that?

thanks in advance.


-- 
You received this question notification because you are a member of
yade-users, which 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 #263739]: isotropic compression at controled pressure

2015-03-16 Thread Luc OGER
Question #263739 on Yade changed:
https://answers.launchpad.net/yade/+question/263739

Status: Answered = Open

Luc OGER is still having a problem:
Luc Thanks for the fast answer

but  this example is not exactly that I was expecting :
 I want to compress by the walls not by growxing the spheres and as the time to 
get equilibrium when moving the walls is crucial I want to select the rate to 
this iso-compression.

if I understand mask=7 and all goal identical are more linked to final
results not to the progressive intermediate steps.

Luc

-- 
You received this question notification because you are a member of
yade-users, which 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 #263739]: isotropic compression at controled pressure

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

Status: Open = Answered

Jérôme Duriez proposed the following answer:
Performing isotropic compression moving the walls at a constant rate is
still possible: you just have to switch internalCompaction to false (*)
(it is true by default)


As for the rate part, using stressMask = 7, the appropriate attribute to set is 
max_vel (**).

Basically, for stress controlled behavior of TriaxialStressController,
this engine tries, moving the walls with internalCompaction=0, to reach
the desired stress value in one iteration. But, without exceeding a wall
velocity of max_vel...

Then, if you want to iso-compress your sample to some pressure
significantly different from the initial one, just use stressMask=7,
internalCompaction=0, and max_vel=someValue, and you will perform an
isotropic compression at a constant strain rate corresponding to
someValue

From my point of view, there is no reason to consider all these
attributes related to final results rather than intermediate steps.


(*) 
https://yade-dem.org/doc/yade.wrapper.html?highlight=triaxialstress#yade.wrapper.TriaxialStressController.internalCompaction
(**) 
https://yade-dem.org/doc/yade.wrapper.html?highlight=triaxialstress#yade.wrapper.TriaxialStressController.max_vel

-- 
You received this question notification because you are a member of
yade-users, which 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 #262899]: YADE and FEM model

2015-03-16 Thread Luis Barbosa
Question #262899 on Yade changed:
https://answers.launchpad.net/yade/+question/262899

Status: Expired = Open

Luis Barbosa is still having a problem:
ok, no problem!
I have a FEM model, that I've created using ABAQUS.
What I would like to do now is open this model using YADE, to simulate its 
interaction with granular material.
Are there some documentation or example of this process?

Thanks a lot.
Luis

-- 
You received this question notification because you are a member of
yade-users, which 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