Real-Time Physics Simulation Forum

trouble with 3D split impulses

2007-03-21T01:26:11+00:00

I found the problem, and it had nothing to do with my code. My body space inverse inertia tensor was wrong. I used the closed form solution for the diagonal inertia tensor for a cube, but I forgot to scale by the mass before inverting it. So the angular impulses weren't being scaled properly.

Also, I found upping the depth slop from 0.01 to 0.1 had a huge impact.

Statistics: Posted by John McCutchan — Wed Mar 21, 2007 1:26 am

trouble with 3D split impulses

2007-03-19T10:14:55+00:00

Conceptually it looks correct what you do, so the only thing I can think of is that you forget to clear the velocities or impulses. Maybe a typo somewhere or you mix accidently the bias impulse with the normal impulse. I would say that it is something trivial that you maybe overlook at the moment...

Statistics: Posted by Dirk Gregorius — Mon Mar 19, 2007 10:14 am

trouble with 3D split impulses

2007-03-19T00:26:14+00:00

Re Dirk:
I don't think it's a sign issue because I'm using the same signs as I do when I'm not using splitting.

Re Erin:
I am setting the bias impulse to zero in PreImpulse and I'm setting
the velocities to zero in Position Integration.

Any other suggestions guys? I'm really baffled.

Statistics: Posted by John McCutchan — Mon Mar 19, 2007 12:26 am

trouble with 3D split impulses

2007-03-19T00:07:01+00:00

Be sure to initialize the bias impulses and velocities to zero before each time step. Test the solver with a sphere on a plane and look at the values in the debugger.

Statistics: Posted by Erin Catto — Mon Mar 19, 2007 12:07 am

trouble with 3D split impulses

2007-03-18T18:59:29+00:00

When you split the impulse you basically have

Velocity impulse:
P = -invK * (J * v)

Position impulse:
P = -invK * 0.1 * C / dt

When you use Baumgarte stabilization you combine the relative velocity at the joint (J * v) with the position error at the joint ( C ). This effects the momentum. What you want is that the stabilization doesn't effect the momentum so you cirvcumvent this by splitting up the impulse into a velocity correcting and (pseudo) position correcting impulse. Note that from experience the splitted impulse seems only to work nice with contacts, but not with joints.

Another thing to understand is that the position constraint is non-linear ( e.g. C = u1 * v2 ) so when you add the bias tem 0.1 * C / dt you linearize the position constraint what opens the door to inaccurecies.

Hmm, from skimming over your code I can't see the error, maybe only a sign. Did you try -erp * B?

Cheers,
-Dirk

Statistics: Posted by Dirk Gregorius — Sun Mar 18, 2007 6:59 pm

trouble with 3D split impulses

2007-03-18T17:46:55+00:00

I have a working 3D version of sequential impulses. I can have 10 cubes stacked directly on top of each
other be stable with only 10 iterations. This is with a erp = 0.2 * fps (fps = 60). I have started
to experiment with splitting the penetration depth correction impulses and I'm not getting good results.
The simulation becomes very jittery with the cubes tilting back and forth and settling into each other.
The instability seems to be coming from the angular movement. Has anyone experienced these kinds of
problems with 3D split impulses. Bullet doesn't use split impulses, is this just not implemented or
was it found to be problematic? I'm going to post my pseudo code below if anyone notices a problem.

PreImpulse:

Code:

bias_normal_lambda = 0.0;B = depth - DEPTH_SLOP;B = B < 0.0 ? 0.0 : B;P = J * accumulated_normal_lambdaapply P

Impulse:

Code:

V = non-bias-velocity of both bodieslambda = K^-1 * (-( J * v))old_lambda = accumulated_normal_lambdaaccumulated_normal_lambda += lambdaaccumulated_normal_lambda = accumulated_normal_lambda < 0.0 ? 0.0 : accumulated_normal_lambdalambda = accumulated_normal_lambda - old_lambdaP = J * lambdaapply P to non-bias velocityV = bias-velocity of both bodiesbias_lambda = K^-1 * (-(J * V) + erp * B)old_bias_lambda = bias_normal_lambdabias_normal_lambda += bias_lambdabias_normal_lambda = bias_normal_lambda < 0.0 ? 0.0 : bias_normal_lambdabias_lambda = bias_normal_lambda - old_bias_lambdaP = J * bias_lambdaapply P to bias velocity

Position Integration code:

Code:

lv = linear_velocity + bias_linear_velocityav = angular_velocity + bias_angular_velocity...bias_linear_velocity = 0.0bias_angular_velocity = 0.0

Statistics: Posted by John McCutchan — Sun Mar 18, 2007 5:46 pm