Not sure if I understand clearly but we don't need to access interaction container, that number will be store at InteractingSphere so we can access it from the constitutive law vector<map<Vector3r, int> should not be so harmful for InteractingSphere.
Those vectors are only considered at interation between sphere and interface in IF2IS4SC and updated only at new contacts so it will not be time consuming. This also solve the smoothing problems. I am not very good at programming so maybe I miss sth? Kien -----Original Message----- From: [email protected] [mailto:[email protected]] On Behalf Of Václav Šmilauer Sent: Thursday, March 19, 2009 1:55 PM To: Yade Development Group Subject: Re: [Yade-dev] facets Dne 19.3.2009 16:55, Kien Dang, Mr napsal(a): > Hi, > According to my understanding, the shrink factor in the paper is used to > avoid the sphere interact with more than one facet, it is; however, not > effective and it may cause unstability if shrink factor is too large and > determine shrink factor suitable for each simulation is somehow unrealistic. > In my opinion, we can remember the number of coincided contact points (in > IF2IS4SC) for each sphere, we can check the conincidence by some criterion > (contactpoint-contactpoint).Length() smaller than some value (can be > tol*sphere radius) and then in the constitutive law the force at the contact > will be divided by that number. By that way we can avoid complicated > algorithms. > That would be much more complicated actually, since the functors shouldn't access interaction containers. The topology information is constant throughout the simulation and it makes not much sense to determine it at every step. It would essentially add another collider to the loop, which would take quite some time. The shrinkFactor doesn't avoid multiple contacts, since there can always be 1 sphere between both facets interacting with both edges. Vaclav _______________________________________________ Mailing list: https://launchpad.net/~yade-dev Post to : [email protected] Unsubscribe : https://launchpad.net/~yade-dev More help : https://help.launchpad.net/ListHelp _______________________________________________ Mailing list: https://launchpad.net/~yade-dev Post to : [email protected] Unsubscribe : https://launchpad.net/~yade-dev More help : https://help.launchpad.net/ListHelp _______________________________________________ yade-dev mailing list [email protected] https://lists.berlios.de/mailman/listinfo/yade-dev
