← Back to team overview

yade-dev team mailing list archive

Re: ?=2 ?! (Ip2_FrictMat_FrictMat_FrictPhys)

 

Ok for the solution below.
Besides, I'm on the way to remove the dyn_casts in Ip2_FM, using genericSpCont, and I see this :

else if(d3dg){Da=d3dg->refR1>0?d3dg->refR1:2*d3dg->refR2; Db=d3dg->refR2>0?d3dg->refR2:d3dg->refR1;

Apparently, the equations are not symmetric and the radius of the virtual sphere will be 2*r "sometimes"... I suggest using 1*r in all cases, but this is the sort of thing that will modify results (at least slightly) for everybody.

Sorry for the crasher that you fixed. On my side, I'm going to commit a fix for callbacks without OpenPM. ;)

Bruno

Václav Šmilauer a écrit :
Ok, you need Young, I was expecting that. I think I will create a new class then.

I found a way to define deprecated attributes (with warnings etc) to
keep older scripts runnable (not yet commited). In that case we could
agree on modifications of ElastMat.
1. What is disturbing for me in what you propose is the difference Kn
vs. kn. That is too close. I wouldn't mind keeping young (as everybody
knows that Young's modulus has dimension of [Pa]; and we might just
document that it is not modulus of the packing etc as we discussed), but
you probably do...

2. Remove poisson, replace it by ksDivKn (I have attributing with the
same significance CpmMat::G_over_E; it is a misnomer, since it is not
directly related to shear modulus; no pain to change that)

FrictMat would additionally define friction angle, just like it does
now.

That way we could keep things same for both of us.

I think the original ElastMat::{young,poisson} were really meant as
continuum parameters, since there was also FEM model in yade IIRC; that
is no longer concern, it seems.

v


_______________________________________________
Mailing list: https://launchpad.net/~yade-dev
Post to     : yade-dev@xxxxxxxxxxxxxxxxxxx
Unsubscribe : https://launchpad.net/~yade-dev
More help   : https://help.launchpad.net/ListHelp



--
_______________
Bruno Chareyre
Associate Professor
Grenoble INP
Lab. 3SR
BP 53 - 38041, Grenoble cedex 9 - France
Tél : 33 4 56 52 86 21
Fax : 33 4 76 82 70 43
________________




Follow ups

References