yade-users team mailing list archive

[Bruno Chareyre <question652968@xxxxxxxxxxxxxxxxxxxxx>]

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

    Status: Open => Answered

Bruno Chareyre proposed the following answer:
>my main concern is the possibility of passing some model parameters to
Iphys (I think there is no need to define a new class of material but
please correct me if I am wrong)

It is possible and indeed simpler usually. [*] Adding the parameters to
a material class gives more flexibility when mixing particles of
different properties.

(1) radius1 is an alias of refR1, you can use both names equivalently
(ScGeom.hpp:58)

(2) I see no need to change current code since it does what you need as
a special case.

(3) ξ can be where you have it now or in the material class, see [*]

(4) same answer. In this case it is obvious that if you want to have
particles of different shapes (different β) you need to have it in the
material class, else all shapes are the same as defined by the Ip2.

(5) if both versions compile they probably do the same thing

(6) I don't see a problem

(7) yes (is it different from the assignement of e.g.
contactPhysics->kno?)

(8) If ξc is like a damage parameter with different evolutions at each contact then it needs to be IPhys::ξc indeed. You don't "declare" it in IP2, you may only assign it there. But even the assignement seems to be in Iphys(.hpp) a.t.m.:
((Real,XI c,2.1,, "Local crushing parameter ..."))

(9) This question sounds like "which is the right constitutive model?". I don't know. This will be your modeling assumption. :)
If you plan to keep viscous damping at 0 it should not be a direct problem for you though.

(10) I would say yes, but I'm not sure about the question. "Plasticity
condition" just means to keep forces/torques bounded by some max values.

(11) I don't know.

(12) No need to erase anything since those lines are conditional already
(l.303,313), your changes should correspond to an additional option.

(13) Newton::damping is the value of Cundall's/PFC "non-viscous" damping
(the one with default 0.7 in PFC). Again, it has nothing to do with the
contact model (it is thus irrelevant to distinguish linear vs. non-
linear models). Newton::damping is set in the scripts.

(14) depends on you

(15) why not? :)

Wow, I made it. ;)
Bruno

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