Question about Ig2_facet_sphere_scgeom.cpp

Asked by MikeZhao

Line 115: scm->radius1 = 2*sphereRadius;

I understand that A facet-sphere interaction is actually converted to a sphere-sphere interaction to calculate contact force.

What I don't understand is why should we equalize a facet to a sphere with a radius 2*sphereRadius? Could any one give me an explanation?

Question information

Language:
English Edit question
Status:
Answered
For:
Yade Edit question
Assignee:
No assignee Edit question
Last query:
Last reply:
Revision history for this message
Jan Stránský (honzik) said :
#1

Hi Mike,

see [1,2]. It is left for backward compatibility. If you need
1*sphereRadius, I guess you can trick it by assigning facet material
stiffness such that the resulting interaction stiffness equals the one that
would be computed with 1*sphereRadius (not tested, just guessing).
cheers
Jan

[1] https://answers.launchpad.net/yade/+question/236412
[2] https://answers.launchpad.net/yade/+question/228778

2016-05-29 16:26 GMT+02:00 MikeZhao <email address hidden>:

> New question #294588 on Yade:
> https://answers.launchpad.net/yade/+question/294588
>
> Line 115: scm->radius1 = 2*sphereRadius;
>
> I understand that A facet-sphere interaction is actually converted to a
> sphere-sphere interaction to calculate contact force.
>
> What I don't understand is why should we equalize a facet to a sphere with
> a radius 2*sphereRadius? Could any one give me an explanation?
>
> --
> You received this question notification because your team yade-users 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
>

Can you help with this problem?

Provide an answer of your own, or ask MikeZhao for more information if necessary.

To post a message you must log in.