In perturbative thermal calculations, introducing nonzero quark masses causes considerable complications. In this article we describe an approximation scheme valid for sufficiently light masses which significantly simplifies the relevant calculations with only a small loss in accuracy. We apply the scheme to quantum chromodynamics, with two massless and one massive flavor, obtaining analytic results which are in excellent agreement with numerical nonapproximated results. Our results are accurate to O(g(5)) in the strong coupling parameter g, as long as either the chemical potential p of the quark species with nonzero mass m satisfies m = O(g mu) or if the temperature T satisfies m = O (g mu). We find that these conditions are always satisfied for the three lightest quarks for the values of mu, T where quantum chromodynamics is perturbatively convergent.