We give an algorithm for obtaining expansions of massive two-loop Feynman g
raphs in powers of the external momentum around a finite, nonzero value of
the momentum. Such momentum derivatives are encountered while evaluating ph
ysical radiative corrections, such as in wave function renormalization cons
tants. This is based on our general two-loop formalism to reduce massive tw
o-loop graphs with renormalizable interactions into a standard set of speci
al functions. After the algebraic reduction, the final results are obtained
by numerical integration. We apply the expansion algorithm to treat the to
p dependent corrections of O(g(2)alpha(s)) to the b-quark self-energy and e
xtract its momentum expansion on-shell.