Maximum Penetration

No, it there is no such option. What behavior do you await, if the
maximum penetration is reached?

Anton

2013/9/26 Ricardo Pieralisi <email address hidden>:
> New question #236400 on Yade:
> https://answers.launchpad.net/yade/+question/236400
>
> Hello all,
> I am using Law2_ScGeom_ViscElPhys_Basic and I am trying to define a maximum penetration (for example a max penetration of Radius*0.5 -> may penetrate only 50% if the radius)...
>
> Is there something ready to do this? if the answer is negative I should to change the constitutive law. Maybe write a condition that identifies if the penetration is higher than the maximum:: if is true -> works only with a very high elastic stiffness in the normal direction, this is possible?
>
> Thank you!
>
> --
> You received this question notification because you are a member of
> yade-users, which is an answer contact for Yade.
>
> _______________________________________________
> Mailing list: https://launchpad.net/~yade-users
> Post to : <email address hidden>
> Unsubscribe : https://launchpad.net/~yade-users
> More help : https://help.launchpad.net/ListHelp

if maximum penetration is reached I want to use just the elastic behauviour (but with a very high stiffness) on the normal direction, and no change in the tangencial.

Is there a paper on the topic?

Anton

2013/9/26 Ricardo Pieralisi <email address hidden>:
> Question #236400 on Yade changed:
> https://answers.launchpad.net/yade/+question/236400
>
> Ricardo Pieralisi posted a new comment:
> if maximum penetration is reached I want to use just the elastic
> behauviour (but with a very high stiffness) on the normal direction, and
> no change in the tangencial.
>
> --
> You received this question notification because you are a member of
> yade-users, which is an answer contact for Yade.
>
> _______________________________________________
> Mailing list: https://launchpad.net/~yade-users
> Post to : <email address hidden>
> Unsubscribe : https://launchpad.net/~yade-users
> More help : https://help.launchpad.net/ListHelp

The contact law defines fn as a function of un. There is really no difficulty in making the relation piecewise linear (I understand this is what you want to do).
The current contact laws don't have this behavior though, so you will need to write a new one, or modify an existing one if it does not interfer too much with its present behavior.
If you want to commit your change to an existing law, requirements are that 1/ the default behavior should not be modified and 2/ it does not imply that users of the current code will have to carry more variables (e.g the secondary stiffness) in the interactions.

The idea is something like this [1], where they have different properties (elastic stiffness and viscous constant) for the external core (more soft) and for the inner core (more rigid)...

I will try to write another law. Is there anything to help with the first steps?

[1] https://www.dropbox.com/s/n7gbkqc1l1252vf/Puri%20-%20Characterization%20of%20Distinct%20Element%20Modeling%20Parameters%20for%20Fresh%20concrete%20and%20its%20Application%20in%20Shotcrete%20Simulations.pdf