basic questions on

Asked by Giulia Macaro on 2010-12-09


I am new to yade and I am starting from (bzr2598). I have included the CohFrictMat and the Coesion moment law implemented by Chiara. But that is not the point.
I am a bit confused regarding how LawTester actually works:

- path() let me apply the displacement starting from the contact point, but I don't understand in which coordinates system

- rotWeight let me decide if the applied shear produces a rotation as well. But when I apply a shear displacement and change rotWeight from 0 to 1, I don't notice any differences.. why?

- plotting un=tester.ptOurs[0], which should be the applied normal displacement, I got a max displ of -0.2, while path=(-0.1, 0, 0).. why? shouldn't they be the same?

I guess all these questions are due to the fact that I don't understand very well path(..).

Thanks in advance

Question information

English Edit question
Yade Edit question
No assignee Edit question
Solved by:
Giulia Macaro
Last query:
Last reply:
Václav Šmilauer (eudoxos) said : #1

Hi Giulia, let me try to answer. I should say at start that it is possible that we will find some bugs in LawTester, as it was written only recently.

1. Local coordinate system; the 0th component is normal (between sphere's centers), 1st and 2nd components are 2 shear directions, perpendicular to the normal. Local system definition depends on IGeom class that is in use. If you use ScGeom (which only knows normal, i.e. local x axis), the other two axes will be chosen arbitrarily at first, then tracked following motion of both particles.

2. That is correct; if there is no rolling, then shear can be created in geometrically different ways, although the resulting shear will be the same -- imagine (i) one sphere fixed, the other one rotating around an axis going through its center and perpendicular to the interaction [that is rotWeight==1] (ii) one sphere fixed, the other one translated around the other one on circular arc centered at the first sphere's center (distance does not change), without underoging any rotation [that is rotWeight==0].

You can assign other values to rotWeight as well (including those ouside the 0...1 range) and the result should still be the same.

Note that rotWeight is ignored when mutual rotations (bending/twisting) are considered by the IGeom, since then, shear must be applied in such a way that it does not create bending and vice versa (that takes away that degree of freedom expressed by rotWeight).

3. Oh yes, they should be the same... I will have a look tomorrow, hope that is fine for you.

If you use ScGeom6D: I will have to add support for it to LawTester (I have not done it yet, and I did not expect someone will use it that soon); OTOH it is not complicated and I will have a look tomorrow morning probably.

Also I did some changes in LawTester today to make it more intuitive (hopefully); your old code will work, and some warnings will tell you how to change it then.

Václav Šmilauer (eudoxos) said : #2

re 3.: Don't you have displIsRel==True? ( That might explain it.

Václav Šmilauer (eudoxos) said : #3

Hi Giulia, I added proper docs for LawTester, I hope they answer your questions now. I am planning to add some images as well, but that has to be postponed now. Cheers!

Giulia Macaro (giulia-macaro) said : #4

Thank you, that is very helpful!

But now I think there are some issues with qt.view, since when I run the script I get the displacements plotted, but I don't see the balls moving in the 3d view (I could see it using the 2598 release).

Václav Šmilauer (eudoxos) said : #5

Sorry, Giu, I fixed that in the last commit (2612), it was a bit of previous changes that broke it.

Giulia Macaro (giulia-macaro) said : #6

Now (2612) it perfectly works.

Thank you