Hi Francois,
thanks for pointing that out. This formulation is in Yade
for several years [1] already. I am not sure, whether it is correct and
we need to check this out. I just keep it "as it was" not to
break a lot of things.
[1] https://github.com/yade/trunk/blob/0.60/pkg/dem/ViscoelasticPM.cpp
Hi Anton,
I checked the Pournin's implementation you made, and I may have found a
little error. According to the formula (10) in the paper, the four formulas
(22) are the resulting stiffness/viscosity_coefs of a contact, and you
rightly transcript them. But in your algorithm, you set
kn1 = kn2 = *
Hi,
One last thing : in the 3rd case (young modulus
case
), it would be useful to determine
cn from E1,E2,massR and en. I know the formula for normal cn, but not for
cs. Do you have an idea ?
Thanks,
François.
2014-04-03 10:15 GMT+02:00 François :
> 2014-04-03 10:00 GMT+02:00 Anton Gladky :
2014-04-03 10:00 GMT+02:00 Anton Gladky :
> 2014-04-03 9:51 GMT+02:00 François :
> > It will for sure break my scripts, as I used massMultiply=0. But your
> > implementation is part of my "todo list", which I can remove now :-)
> > You have made huge changes at once, so I need a little bit of time
Hi, and thank you for warning me.
It will for sure break my scripts, as I used massMultiply=0. But your
implementation is part of my "todo list", which I can remove now :-)
You have made huge changes at once, so I need a little bit of time to check
everything. What I can see now is that you use "i
2014-04-03 9:51 GMT+02:00 François :
> It will for sure break my scripts, as I used massMultiply=0. But your
> implementation is part of my "todo list", which I can remove now :-)
> You have made huge changes at once, so I need a little bit of time to check
> everything. What I can see now is that
Hi Francois,
could you, please, test the latest git-version of Yade especially
your scripts with ViscoElastic model? I have done there some
changes, but I believe, that your code is not broken.
If you want, we can also implement MatchMaker for some of
variables, like it is already done for tc, e
7 matches
Mail list logo