I have two boxes that are constrained by a RevoluteJoint. I want to know the internal forces that are exerted on the joint, satisfying the constrain. It should be in the form of a direction plus the magnitude of the force.
Yup. Note that the TotalImpulse properties on constraints are in constraint space. That's not much of a concern if all you're looking for is magnitude, but to get the world space force applied to each entity you'd need to multiply the impulse by the jacobians (GetLinearJacobianA/B and GetAngularJacobianA/B).
For an example of how you use the jacobians to get to world space impulses, check out something like the PointOnLineJoint.ExclusiveUpdate function. The constraint space impulse (here accumulatedImpulse) is transformed by the linear jacobians worldRestrictedAxis1 and worldRestrictedAxis2 to get the linear world impulse, and then transformed by angularA1/A2/B1/B2 to get angular world impulse.