Does it seem correct to you?

I'm realizing now why we've had confusion with all these rotations. In the Crandall book I think they define rotation matrices as vb = R * va, where as in Kane he does va = R * vb. I don't have my book handy to check right now.

I think you used Crandall for this inerita rotation formula but defined everything else as Kane does.