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. Such momentum derivatives are encountered while evaluating physical radiative corrections, such as in wave function renormalization constants. 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 [Formula: see text] to the b-quark self-energy and extract its momentum expansion on-shell.