For torsion on wide flanges I always use the "equivalent tee" analogy as a starting point. The reason why I like it is becasue it gives a good physical explanation for where those warping stresses are coming from.
Here is the procedure:
1) For torsion, think of the I beam as having been split at the center. Hence forming two Tee beams.
2) If the original beam was warping restrained at the ends, then your tee beam will be a fixed-fixed beam. If your original beam did NOT restrain warping, then your tee beams will be simply supported.
3) Break your torsional load into a force couple applied at the centroid of the Tee. Apply these point load or distributed load to each tee beam individually and analyze using simple beam theory.
4) The simple beam theory will develop weak axis bending moments (and flange bending stresses) in the Tee that are similar to the warping stresses that you'd get from a more rigorous analysis.
This analogy is relatively simple and should always be conservative... at least compared to the actual warping stresses in the flanges. Therefore, there is little danger to using it.
One of the things that I like it this method can be expanded out in situations where the rigorous methods become too difficult. Such as when you want to model a boundary condition that is somewhere between warping fixity and a warping pin. Or, when you have multiple unusual loading conditions.
Josh