For $m>0$, $\displaystyle\frac{\sin axm}{\sin ax}=\sum_{k=0}^{m-1}\cos(m-1-2k)ax$, so the given integral is a sum of things like $$\int_0^\infty\frac{\cos bx}{1+e^{ax}}\,dx=\frac{1}{2a}\left[f\left(\frac{ib}{a}\right)+f\left(-\frac{ib}{a}\right)\right]\tag{*}\label{mainint}$$ where, for a complex $z$ with $\Re z>-1$, $$f(z)=\int_0^\infty\frac{e^{-zx}}{1+e^x}\,dx\underset{e^{-x}=t}{=}\int_0^1\frac{t^z\,dt}{1+t}=\int_0^1\frac{t^z-t^{z+1}}{1-t^2}\,dt\\\underset{t^2=x}{=}\frac12\int_0^1\frac{x^{(z-1)/2}-x^{z/2}}{1-x}\,dx=\frac12\left[\psi\left(1+\frac{z}{2}\right)-\psi\left(\frac{1+z}{2}\right)\right],$$ with $\psi$ the digamma function (the final equality is shown like it's done here). If we had sine in place of the cosine in $\eqref{mainint}$, the $\psi$'s would reduce because of the reflection formula. With the cosine in place, these don't, as well as in the final result. That's why I don't expect contour integration to yield anything useful.