Friction impulse. Need formula.

Started by
12 comments, last by Volgogradetzzz 13 years ago
Hello. For normal impulse I used Chris Heckers papers. All works fine, but formula don't take friction into account. I believe that there's equivalent for friction (not just tangent velocity * cof) but I couldn't google it.
Advertisement

Hello. For normal impulse I used Chris Heckers papers. All works fine, but formula don't take friction into account. I believe that there's equivalent for friction (not just tangent velocity * cof) but I couldn't google it.


Friction is complicated. In the naive sense, rigid body dynamics with a Coulombic friction force (which is what you are describing) is not even solvable. The classic and ancient example of this is Painleve's paradox, which has remained unresolved until relatively recently.

As a result, if you want to do friction using the classical Coulomb approximation, you are in for a very rude awakening (as I am sure you have just discovered!) A more well posed way to formulate friction is to use the principle of maximum dissipation. This basically says that the friction force will try to reduce the velocity (ie kinetic energy) as much as possible in the tangent direction to the collision impulse force, subject to some constraints on the magnitude of the velocity. This constraint forms a "cone" shape, and so it is commonly called a "friction cone". To solve this constraint, many LCP solvers use a piecewise linear approximation of the friction cone to ensure the constraint is not violated. For a more precise and careful discussion of these issues, there are some really great papers by David E. Stewart. Here is a good one to start with (in my opinion):


David E. Stewart, (2000) "Rigid Body Dynamics with Friction and Impact", SIAM Review PDF
Just take the contact point velocity, V and normal N, form the tangential velocity

T = V - V.N*N

Then create an impulse to remove some portion of this velocity in the same way you did to remove the normal velocity :)

Cheers, Paul.
You may find my post useful:

http://www.gamedev.net/topic/596977-collisions-involving-friction-explained-partly/

I wrote up a description of Brian Mirtich's algorithm. For 2D (non-resting) collisions there isn't a single formula but there's a quick algorithm. For 3D collisions the algorithm involves numeric integration which may be too slow for your application.
In 3d, you usually need 2 tangential direction, which must exist even if the tangential velocity is zero. and then, the impulse calculation is similar to normal impulse calculation, just change the normal vector with tangential vector. And then you relate the friction impulse with the normal impulse
Thanks everyone. I was affraid that it's complicated and so it is. I'm certainly investigate papers that you're link, but for now I think I'll take approach that wildbunny and Bow_vernon said.
And one more question - is it OK that if I take COF (coefficient of friction) zero (like with COR, 0 - is 'full' friction, i.e. no sliding), but tangent velocity decreases not by 100%, but only by some part? I've used same formula as for normal impulse - only normal was changed to tangent direction...
Update: also notice some problem - if two objects fall down (their velocities is like [0; 10] vector) and are close one to another (so collision occurs every frame) then friction decreases velocity and bodies become to fall slowly :blink:.
Well you've arrived there. with 0 CoF the object should slide, you gotta relate the friction impulse magnitude with the normal impulse, like:
fMax = fN * CoF
then you clamp the friction so its value is between [-fMax,fMax]. Anyway, how did you calculate the friction impulse? can you show some code?
Normal impulse:
var jNnum:Number = -(1 + e) * (vab.dot(n));
var jNdenom:Number = bodyA.invMass + bodyB.invMass + rap.cross(n) * rap.cross(n) * bodyA.invInertia + rbp.cross(n) * rbp.cross(n) * bodyB.invInertia;
jN = jNnum / jNdenom;

Friction impulse is all the same. Only changes are:
e -> f; //cof
n -> t; //tangent normal

It's hard to understand - why shoud we multiply fMax = fN * CoF? How to achieve full friction (i.e. no slide) in this case - suppose I want bodies to stop sliding immidiatelly after first contact?
What I do(has been simplified)

jn = kN * -vn; //no bounce
jn = Max(0,jn); //no pull

fMax = jn * CoF;
jt = kT * -vt; //full stop
jt = Clamp(jt, -fMax, fMax);

Oh and read erin catto's ppt for more info, anyway. Sorry if I didnt help, Im confused @.@
Ok. I think I got it. Thank you.

This topic is closed to new replies.

Advertisement