  We give an algorithm for obtaining expansions of massive two-loop Feynman
graphs in powers of the external momentum around a finite, nonzero value of the
momentum. This is based on our general two-loop formalism to reduce massive
two-loop graphs with renormalizable interactions into a standard set of special
functions. After the algebraic reduction, the final results are obtained by
numerical integration. We apply the expansion algorithm to treat the
top-dependent corrections of O(g^2 alpha_s) to the b quark self-energy and
extract its momentum expansion on-shell.
