Re: [Yade-users] [Question #682859]: how to calculate anisotropy of contact normal vector

2019-08-14 Thread William Chevremont
Question #682859 on Yade changed:
https://answers.launchpad.net/yade/+question/682859

Status: Open => Answered

William Chevremont proposed the following answer:
Hi,

PDFEngine do that for you. It compute the average of quantities like
stress, distance, and probability of interaction for each orientation.
For now, it works only for spherical particles.

Just add 
PDFEngine(filename="PDF.txt", iterPeriod=1000, 
numDiscretizeAnglePhi=9,numDiscretizeAngleTheta=13, label="PDFRecorder")

in your engines. This will output a PDF.txt containing the density
functions, with column name like NormalStress_xy(theta, phi)

Formulas: 
x = r sin theta sin phi
y = r cos phi
z = r sin theta cos phi

This is a work in progress (not well documented for now) but it work as
is.

-- 
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 #682859]: how to calculate anisotropy of contact normal vector

2019-08-14 Thread Jan Stránský
Question #682859 on Yade changed:
https://answers.launchpad.net/yade/+question/682859

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

as (almost) always, a MWE [1] would be useful (e.g. to know what
"particle" is - sphere, polyhedra...), expressing on a few artificially
created spheres and artificially created interactions to show what you
want to achieve and what is your expectation of "anisotropy of contact
normal vector", etc.

Often you could solve your problem just by preparing such MWE.
Furthermore, it will help us to help you :-)

> how to calculate the contact normal vector

for i in O.interactions:
   contactNormalVector = i.geom.normal # [2]

> i want to know how to calculate anisotropy of contact normal vector

if I understand it correctly, fabric tensor (as suggested by Bruno) is
what you are looking for [3]

cheers
Jan

[1] https://www.yade-dem.org/wiki/Howtoask
[2] 
https://yade-dem.org/doc/yade.wrapper.html#yade.wrapper.GenericSpheresContact.normal
[3] https://yade-dem.org/doc/yade.utils.html#yade._utils.fabricTensor

-- 
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 #682857]: how to calculate anisotropy of contact normal vector

2019-08-14 Thread Jan Stránský
Question #682857 on Yade changed:
https://answers.launchpad.net/yade/+question/682857

Status: Open => Invalid

Jan Stránský rejected the question:
duplicate of #682859

-- 
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 #682883]: calculate mass content of a group of particles through something like label

2019-08-14 Thread Leonard
New question #682883 on Yade:
https://answers.launchpad.net/yade/+question/682883

Hi,
I'd like to calculate the mass content of a group of spheres, as far as I know, 
we can obtain a mass of one sphere by:
O.bodies[id].state.mass
I'd like to ask that if we can obtain the mass of a group of spheres like: 
appending some spheres with a same label, like:
O.append(sphere((x1,y1,z1),radius=radius,label="label1"))
O.append(sphere((x2,y2,z2),radius=radius,label="label1"))
O.append(sphere((x3,y3,z3),radius=radius,label="label2"))
O.append(sphere((x4,y4,z4),radius=radius,label="label2"))
then we can obtain the total mass of spheres with label1, and label2 
respectively?
I have tried in Yade, it seems that we can not give each sphere a label, is 
there any way that we can classify the spheres to different groups so that we 
can just loop in one group of spheres to calculate the mass?
Thanks 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 #682883]: calculate mass content of a group of particles through something like label

2019-08-14 Thread gaoxuesong
Question #682883 on Yade changed:
https://answers.launchpad.net/yade/+question/682883

Status: Open => Answered

gaoxuesong proposed the following answer:
Hi,

The possible way i think is that using different lists to save the sphere's 
id of different groups and do loop in each list. 
## group 1 ##  
list1 = [] 
list1.append(O.append(sphere((x1,y1,z1),radius=radius1)))
... 
## group 2 ## 
list2 = [] 
list2.append(O.append(sphere((x2,y2,z2),radius=radius2)))
... 
## calculate the mass of different groups ## 
masstotal1 = 0 
for each in list1:
masstotal1 + = O.bodies[each].state.mass 
masstotal2 = 0 
for each in list2:
masstotal2 + = O.bodies[each].state.mass 

## ends ##

Cheers, 
Xuesong

-- 
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 #682898]: Which model is used to calculate the cohesive force?

2019-08-14 Thread gaoxuesong
New question #682898 on Yade:
https://answers.launchpad.net/yade/+question/682898

Hi. 

I have looked through the document for 
Law2_ScGeom6D_CohFrictPhys_CohesionMoment. It gives the formula to calculate 
the contact force. Fn=min(kn∗(un−upn),an). an is the cohesive force. But is 
there any material to refer to get the detailed description about cohesive 
force, an? 

Thanks, 
Xuesong 

-- 
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 #682898]: Which model is used to calculate the cohesive force?

2019-08-14 Thread gaoxuesong
Question #682898 on Yade changed:
https://answers.launchpad.net/yade/+question/682898

gaoxuesong posted a new comment:
I have checked the definition of CohFrictMat. So the normalCohesion
and shearCohesion is the an and as appeared in the documentation of
Law2_ScGeom6D_CohFrictPhys_CohesionMoment?

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