date:20160922

Re: [Yade-users] [Question #399298]: Machmaker in CPM model

2016-09-22 Thread Seti

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

Status: Answered => Open

Seti is still having a problem:
Hi Jan,

Thanks so much for your assist and sorry for slow reply.

I just sent an email to you. As I have mentioned in my email I need to
control all of the parameters that have effect on Cohesion between
spheres in CPM model.

Also I need to use matchmaker for inter-granular friction between
spheres.


Appreciate for that.

Seti

-- 
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 #290070]: TesselationWrapper and porosity

2016-09-22 Thread Weimin Song

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

Weimin Song posted a new comment:
Hey behzad,

Do you solve the problem?

-- 
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 #401682]: Problem with angular velocity

2016-09-22 Thread Jonathan Pergoli

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

Status: Answered => Solved

Jonathan Pergoli 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 #401682]: Problem with angular velocity

2016-09-22 Thread Jonathan Pergoli

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

Jonathan Pergoli posted a new comment:
Hi Jan,

sorry I have made a mess with the code :)

I have used O.bodies.append(body).

Anyway, I have solved the problem using angMom :)

Thank you,

Jonathan

-- 
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 #399298]: Machmaker in CPM model

2016-09-22 Thread Jan Stránský

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

Jan Stránský proposed the following answer:
Hello,
could you please send me an email what parameters you want to make
MatchMaker? I will modify the files for you, but I will not apply them in
the official version (since it would need some testing that the default
behavior does not change etc.)
cheers
Jan


2016-09-19 18:17 GMT+02:00 Seti :

> Question #399298 on Yade changed:
> https://answers.launchpad.net/yade/+question/399298
>
> Status: Answered => Open
>
> Seti is still having a problem:
> Thanks Jan for kind support. I have downloaded yade from
> "https://github.com/yade/trunk;
>
> is that what you have recommended?
>
> if yes,
> 1- How I can install it?
> installing package was really easy just copy/past the mentioned commend
> lines from "https://yade-dem.org/doc/installation.html; however I am not
> sure how should I install the source codes. would you please assist me ?
>
> 2- I have found the source code of CPM on truck and I suppose I should
> change below lines to add matchmaker to control the cohesion. so dose it
> mean if I do changes in "  trunk/pkg/dem/ConcretePM.hpp" and save it
> then run YADE it will work?
>
> ###
> class Ip2_CpmMat_CpmMat_CpmPhys: public IPhysFunctor{
> public:
> virtual void go(const shared_ptr& pp1, const
> shared_ptr& pp2, const shared_ptr& interaction);
> FUNCTOR2D(CpmMat,CpmMat);
> DECLARE_LOGGER;
> 
> YADE_CLASS_BASE_DOC_ATTRS(Ip2_CpmMat_CpmMat_CpmPhys,IPhysFunctor,"Convert
> 2 :yref:`CpmMat` instances to :yref:`CpmPhys` with corresponding
> parameters. Uses simple (arithmetic) averages if material are different.
> Simple copy of parameters is performed if the :yref:`material` is
> shared between both particles. See :yref:`cpm-model` for detals.",
> ((long,cohesiveThresholdIter,10,,"Should new
> contacts be cohesive? They will before this iter#, they will not be
> afterwards. If 0, they will never be. If negative, they will always be
> created as cohesive (10 by default)."))
> ((shared_ptr,E,,,"Instance of
> :yref:`MatchMaker` determining how to compute interaction's normal modulus.
> If ``None``, average value is used."))
> );
> };
>
> ##
>
> 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
>

-- 
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 #401682]: Problem with angular velocity

2016-09-22 Thread Jan Stránský

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

Jan Stránský proposed the following answer:
Hi Jonathan,
next time please try the script that it really works ;-) and please doulble
check if both scripts are exactly the same for the two versions (especially
O.bodies.append, see below)

To get some effect, you should use O.bodies.append(body) so that the
engines know about the particle. Do you use it in your scripts or not?

If you don't add body to O.bodies, than angVel does not change (but the
body does not move at the same time, because Yade engines do not know about
the particle)

With O.bodies.append(body), angVel is reset after O.step(), which is
expected (see e.g. [1] and links from there). Setting angMom is the
solution. The relation is angMom = inertiaTensor*angVel (similar to
momentum=mass*velocity). If this works for you, it should not be a problem
to correctly define inertiaTensor.

cheers
Jan

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


2016-09-22 12:57 GMT+02:00 Jonathan Pergoli <
question401...@answers.launchpad.net>:

> Question #401682 on Yade changed:
> https://answers.launchpad.net/yade/+question/401682
>
> Jonathan Pergoli posted a new comment:
> Hi Jan,
>
> sorry I will provide the working script:
>
> E1=1e+8
> E2=5e+7
> mat=FrictMat(density=643,frictionAngle=0.1489,label="MLI",young=E2)
> MLI=O.materials.append(mat)
> a=.2754
> b=.2822
> c=.1963
> aa=a/2
> bb=b/2
> cc=c/2
> h=.65
> dist=0
> theta=0
> thetav=0
> v1=(aa,bb,c)
> v2=(aa,-bb,c)
> v3=(-aa,-bb,c)
> v4=(-aa,bb,c)
> v5=(aa,bb,0)
> v6=(aa,-bb,0)
> v7=(-aa,-bb,0)
> v8=(-aa,bb,0)
> V=[v1,v2,v3,v4,v5,v6,v7,v8]
> vz=.19
> R=[[m.cos(theta),0,m.sin(theta)],[0,1,0],[-m.sin(theta),0,m.cos(theta)]]
> v1=(R[0][0]*V[0][0]+R[0][1]*V[0][1]+R[0][2]*V[0][2],R[1][0]*
> V[0][0]+R[1][1]*V[0][1]+R[1][2]*V[0][2],R[2][0]*V[0][0]+R[
> 2][1]*V[0][1]+R[2][2]*V[0][2])
> v2=(R[0][0]*V[1][0]+R[0][1]*V[1][1]+R[0][2]*V[1][2],R[1][0]*
> V[1][0]+R[1][1]*V[1][1]+R[1][2]*V[1][2],R[2][0]*V[1][0]+R[
> 2][1]*V[1][1]+R[2][2]*V[1][2])
> v3=(R[0][0]*V[2][0]+R[0][1]*V[2][1]+R[0][2]*V[2][2],R[1][0]*
> V[2][0]+R[1][1]*V[2][1]+R[1][2]*V[2][2],R[2][0]*V[2][0]+R[
> 2][1]*V[2][1]+R[2][2]*V[2][2])
> v4=(R[0][0]*V[3][0]+R[0][1]*V[3][1]+R[0][2]*V[3][2],R[1][0]*
> V[3][0]+R[1][1]*V[3][1]+R[1][2]*V[3][2],R[2][0]*V[3][0]+R[
> 2][1]*V[3][1]+R[2][2]*V[3][2])
> v5=(R[0][0]*V[4][0]+R[0][1]*V[4][1]+R[0][2]*V[4][2],R[1][0]*
> V[4][0]+R[1][1]*V[4][1]+R[1][2]*V[4][2],R[2][0]*V[4][0]+R[
> 2][1]*V[4][1]+R[2][2]*V[4][2])
> v6=(R[0][0]*V[5][0]+R[0][1]*V[5][1]+R[0][2]*V[5][2],R[1][0]*
> V[5][0]+R[1][1]*V[5][1]+R[1][2]*V[5][2],R[2][0]*V[5][0]+R[
> 2][1]*V[5][1]+R[2][2]*V[5][2])
> v7=(R[0][0]*V[6][0]+R[0][1]*V[6][1]+R[0][2]*V[6][2],R[1][0]*
> V[6][0]+R[1][1]*V[6][1]+R[1][2]*V[6][2],R[2][0]*V[6][0]+R[
> 2][1]*V[6][1]+R[2][2]*V[6][2])
> v8=(R[0][0]*V[7][0]+R[0][1]*V[7][1]+R[0][2]*V[7][2],R[1][0]*
> V[7][0]+R[1][1]*V[7][1]+R[1][2]*V[7][2],R[2][0]*V[7][0]+R[
> 2][1]*V[7][1]+R[2][2]*V[7][2])
> body=utils.polyhedron((v1,v2,v3,v4,v5,v6,v7,v8),fixed=
> False,color=(.6,.45,0),material="MLI",wire=False)
> body.state.angVel=(0,0,0.1)
> r=m.sqrt(aa**2+bb**2)
> Rj=m.sqrt(r**2+cc**2)
> Ri=0.05
> Rr=Rj*Ri/(Rj+Ri)
> mu_rM=0.05
> mu_rG=1.08
> KN=E1*2*Ri*E1*2*Ri/(E1*2*Ri+E1*2*Ri)
> KR=3*Ri**2*mu_rG**2*KN/4
>
> O.engines=[
> ForceResetter(),
> InsertionSortCollider([Bo1_Polyhedra_Aabb()]),
> InteractionLoop(
> [Ip2_FrictMat_FrictMat_MindlinPhys(en=.45,es=.45,
> krot=KR,frictAngle=.7853)],
> [Law2_ScGeom_MindlinPhys_Mindlin(includeMoment=True)]
> ),
> NewtonIntegrator(gravity=(0,0,-2.5e-4),damping=0),
> ]
>
> The script is the same and also the version of yadedaily.
>
> I'll try with angMom.
>
> Thank you,
>
> Jonathan
>
> --
> 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
>

-- 
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 #401682]: Problem with angular velocity

2016-09-22 Thread Jonathan Pergoli

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

Jonathan Pergoli posted a new comment:
Hi Jan,

sorry I will provide the working script:

E1=1e+8
E2=5e+7
mat=FrictMat(density=643,frictionAngle=0.1489,label="MLI",young=E2)
MLI=O.materials.append(mat)
a=.2754
b=.2822
c=.1963
aa=a/2
bb=b/2
cc=c/2
h=.65
dist=0
theta=0
thetav=0
v1=(aa,bb,c)
v2=(aa,-bb,c)
v3=(-aa,-bb,c)
v4=(-aa,bb,c)
v5=(aa,bb,0)
v6=(aa,-bb,0)
v7=(-aa,-bb,0)
v8=(-aa,bb,0)
V=[v1,v2,v3,v4,v5,v6,v7,v8]
vz=.19
R=[[m.cos(theta),0,m.sin(theta)],[0,1,0],[-m.sin(theta),0,m.cos(theta)]]
v1=(R[0][0]*V[0][0]+R[0][1]*V[0][1]+R[0][2]*V[0][2],R[1][0]*V[0][0]+R[1][1]*V[0][1]+R[1][2]*V[0][2],R[2][0]*V[0][0]+R[2][1]*V[0][1]+R[2][2]*V[0][2])
v2=(R[0][0]*V[1][0]+R[0][1]*V[1][1]+R[0][2]*V[1][2],R[1][0]*V[1][0]+R[1][1]*V[1][1]+R[1][2]*V[1][2],R[2][0]*V[1][0]+R[2][1]*V[1][1]+R[2][2]*V[1][2])
v3=(R[0][0]*V[2][0]+R[0][1]*V[2][1]+R[0][2]*V[2][2],R[1][0]*V[2][0]+R[1][1]*V[2][1]+R[1][2]*V[2][2],R[2][0]*V[2][0]+R[2][1]*V[2][1]+R[2][2]*V[2][2])
v4=(R[0][0]*V[3][0]+R[0][1]*V[3][1]+R[0][2]*V[3][2],R[1][0]*V[3][0]+R[1][1]*V[3][1]+R[1][2]*V[3][2],R[2][0]*V[3][0]+R[2][1]*V[3][1]+R[2][2]*V[3][2])
v5=(R[0][0]*V[4][0]+R[0][1]*V[4][1]+R[0][2]*V[4][2],R[1][0]*V[4][0]+R[1][1]*V[4][1]+R[1][2]*V[4][2],R[2][0]*V[4][0]+R[2][1]*V[4][1]+R[2][2]*V[4][2])
v6=(R[0][0]*V[5][0]+R[0][1]*V[5][1]+R[0][2]*V[5][2],R[1][0]*V[5][0]+R[1][1]*V[5][1]+R[1][2]*V[5][2],R[2][0]*V[5][0]+R[2][1]*V[5][1]+R[2][2]*V[5][2])
v7=(R[0][0]*V[6][0]+R[0][1]*V[6][1]+R[0][2]*V[6][2],R[1][0]*V[6][0]+R[1][1]*V[6][1]+R[1][2]*V[6][2],R[2][0]*V[6][0]+R[2][1]*V[6][1]+R[2][2]*V[6][2])
v8=(R[0][0]*V[7][0]+R[0][1]*V[7][1]+R[0][2]*V[7][2],R[1][0]*V[7][0]+R[1][1]*V[7][1]+R[1][2]*V[7][2],R[2][0]*V[7][0]+R[2][1]*V[7][1]+R[2][2]*V[7][2])
body=utils.polyhedron((v1,v2,v3,v4,v5,v6,v7,v8),fixed=False,color=(.6,.45,0),material="MLI",wire=False)
body.state.angVel=(0,0,0.1)
r=m.sqrt(aa**2+bb**2)
Rj=m.sqrt(r**2+cc**2)
Ri=0.05
Rr=Rj*Ri/(Rj+Ri)
mu_rM=0.05
mu_rG=1.08
KN=E1*2*Ri*E1*2*Ri/(E1*2*Ri+E1*2*Ri)
KR=3*Ri**2*mu_rG**2*KN/4

O.engines=[
ForceResetter(),
InsertionSortCollider([Bo1_Polyhedra_Aabb()]),
InteractionLoop(

[Ip2_FrictMat_FrictMat_MindlinPhys(en=.45,es=.45,krot=KR,frictAngle=.7853)],
[Law2_ScGeom_MindlinPhys_Mindlin(includeMoment=True)]
),
NewtonIntegrator(gravity=(0,0,-2.5e-4),damping=0),
]

The script is the same and also the version of yadedaily.

I'll try with angMom.

Thank you,

Jonathan

-- 
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 #401682]: Problem with angular velocity

2016-09-22 Thread Jan Stránský

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

Status: Open => Answered

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

please provide a working script, i.e. with definitions of v1-v8, MLI etc.
such that we can try what does not work.

do you use the same versions of yadedaily on both machines?

do you use the same script on both machines (not like polyhedron on one and
sphere on the other)?

meanwhile try to set angMom instead of angVel for non-spherical shapes.

cheers
Jan


2016-09-21 17:53 GMT+02:00 Jonathan Pergoli <
question401...@answers.launchpad.net>:

> New question #401682 on Yade:
> https://answers.launchpad.net/yade/+question/401682
>
> Hi everybody!
>
> I'm trying to set an initial value of the angular velocity for a body. I
> set the angular velocity:
> body.state.angVel=(0,0,0.1)
>
> Then when I start yadedaily at the beginning of the simulation the angular
> velocity is correct but when I write O.run() I don't know why the angular
> velocity is set to (0,0,0). This happen when I use yadedaily for ubuntu
> server 14.04 LTS whereas when I run the simulation with ubuntu 14.04 on my
> laptop it works well and this problem is not present. Any solutions?
>
> Thank you in advance,
>
> Jonathan
>
> --
> 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
>

-- 
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 #399298]: Machmaker in CPM model

Re: [Yade-users] [Question #290070]: TesselationWrapper and porosity

Re: [Yade-users] [Question #401682]: Problem with angular velocity

Re: [Yade-users] [Question #401682]: Problem with angular velocity

Re: [Yade-users] [Question #399298]: Machmaker in CPM model

Re: [Yade-users] [Question #401682]: Problem with angular velocity

Re: [Yade-users] [Question #401682]: Problem with angular velocity

Re: [Yade-users] [Question #401682]: Problem with angular velocity

8 matches

Site Navigation

Mail list logo

Footer information