**0**

# Friction impulse. Need formula.

###
#2
Members - Reputation: **193**

Posted 24 March 2011 - 10:28 AM

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

###
#4
Members - Reputation: **104**

Posted 24 March 2011 - 06:26 PM

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.

###
#5
Members - Reputation: **238**

Posted 25 March 2011 - 12:22 AM

###
#6
Members - Reputation: **1078**

Posted 25 March 2011 - 01:44 AM

**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 .

###
#7
Members - Reputation: **238**

Posted 25 March 2011 - 05:57 AM

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?

###
#8
Members - Reputation: **1078**

Posted 25 March 2011 - 06:35 AM

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?

###
#9
Members - Reputation: **238**

Posted 25 March 2011 - 07:05 AM

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 @.@

###
#14
Members - Reputation: **1078**

Posted 26 March 2011 - 01:54 AM

**Bow_vernon**said, for 2D friction can be count like this:

1. Get normal impulse - see post #8.

2. Get tangent impulse:

var jTnum:Number = -vab.dot(t); var jTdenom:Number = bodyA.invMass + bodyB.invMass + rap.cross(t) * rap.cross(t) * bodyA.invInertia + rbp.cross(t) * rbp.cross(t) * bodyB.invInertia; jT = jTnum / jTdenom;

3. You need to clamp friction impulse:

if(jT < -jN * cof || jT > jN * cof) { jT = jN * cof; }

4. Get total impulse as:

j = n * jN + t * jT