yade-users team mailing list archive

[Jan Stránský <question698605@xxxxxxxxxxxxxxxxxxxxx>]

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

    Status: Open => Answered

Jan Stránský proposed the following answer:
> By uniform, I mean that the height of the polar histogram is almost
the same.

just shift: how do you define "almost the same"?

> Since the function utils.plotDirections can only display the contact
distribution in a semicircle, I made the contact distribution to be
displayed in the entire circle by symmetrical processing.

since the other semicircle is symmetric, is it worth to change a
"standard" function? what are the results of the original
utils.plotDirections? same? better? worse?

> I did that to calculate the contact normal force distribution defined
by Rothenburg[1]

I do not have access, can you provide the formula / method description?

> I want to know if there is a problem in my isostatic consolidation
process, which leads to uneven contact and contact normal force
distribution.

I would first suspect the plotting code. As discussed:
- the same forces distribution for the three planes is very suspicious
- the "randomness of bins of perpendicular forces" does not look very rigorous
- different computation for direction and forces, but it is intrinsically the same. For directions, are sphere-wall interactions included or not?
- I don't like the bin-wise averaging
- do you have some testing sample for which you know the result? Such test would strongly influence the code "confidence"
- ...

Cheers
Jan

-- 
You received this question notification because your team yade-users is
an answer contact for Yade.