An algorithm for simulating simultaneous collision impulse with friction working at multiple collision points between three dimensional rigid bodies is described. We also present that lasting collisions at small interval will occur under some dynamic conditions and cause the decline in the simulation speed, and the method for avoiding this phenomenon is described.