yade-users team mailing list archive

[jduriez <question158839@xxxxxxxxxxxxxxxxxxxxx>]

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

    Status: Open => Answered

jduriez proposed the following answer:
Hi Chiara,

I used this kind of algorithm in my simulations with non linear contac relationships : to maintain a contact normal stress on a box, while contact law was "bilinear" (kn(loading) = cst # kn(unloading) = cst ). Probably different from your case, but for me it worked.
I paid attention to use the right stifness in the algorithm, using kn(loading) if my box goes into the sample, and the contrary otherwise.

You could maybe use an updated stifness value in your case (with
derivative ?) ??? if the changes of stifness during interpenetration is
really the problem

Not Bruno (sorry ;-) )

-- 
You received this question notification because you are a member of
yade-users, which is an answer contact for Yade.