Nota: Esta no es una respuesta completa, pero incluye algunas observaciones que, posiblemente, puede ser desarrollado aún más.
Contorno de la Representación Integral de la ${\rm Li}_\nu(e^{-x})$
Por Mellin de la inversión teorema, la polylogarithm admite el contorno de la representación integral
$${\rm Li}_\nu(e^{-x})=\frac{1}{2\pi i}\int^{c+i\infty}_{c-i\infty}\Gamma(s)\zeta(s+\nu)x^{-s}\ {\rm d}s$$
La Suma de $S(\mu)$ como una Integral de Contorno
Cambiar el orden de resumen, podemos escribir
$$S(\mu)\stackrel{\text{def}}=:\sum_{n\ge1}\frac{1}{n\sinh(n\mu)}=2\sum_{n\ge1}\frac{1}{n}\sum_{k\ge0}e^{-(2k+1)n\mu}=2\sum_{k\ge0}{\rm Li}_1\left(e^{-(2k+1)\mu}\right)$$
Por lo tanto, la suma puede ser expresada como una integral de contorno
$$S(\mu)=\frac{1}{2\pi i}\int^{c+i\infty}_{c-i\infty}Q(s,\mu)\ {\rm d}s$$
donde
\begin{align}
Q(s,\mu)
&=2\mu^{-s}(1-2^{-s})\Gamma(s)\zeta(s)\zeta(s+1)\\
\end{align}
Cierre el Contorno
Que en lugar de considerar la integral de $\dfrac{Q(s,\mu)}{2\pi i}$ más de un contorno cerrado $C$. En este caso, es un rectángulo con vértices $\frac{3}{2}-i\infty$, $\frac{3}{2}+i\infty$, $-\frac{3}{2}+i\infty$, $-\frac{3}{2}-i\infty$. Tenga en cuenta que
- La integral sobre la $\mathrm{Re}(s)=\frac{3}{2}$$S(\mu)$.
- La integral sobre los lados horizontales se desvanece desde $\zeta(s)$ crece moderadamente, sino $\Gamma(s)$ decae exponencialmente $\mathrm{Im}(s)\to-\infty$.
- La aplicación de Riemann funcional de la ecuación dos veces, la integral sobre la $\mathrm{Re}(s)=-\frac{3}{2}$ se muestra:
\begin{align}
-\frac{1}{2\pi i}\int^{-\frac{3}{2}+i\infty}_{-\frac{3}{2}-i\infty}Q(\mu,s)\ {\rm d}s
&=-\frac{1}{2\pi i}\int^{-\frac{3}{2}+i\infty}_{-\frac{3}{2}-i\infty}2\left(\frac{\mu}{4\pi^2}\right)^{-s}(1-2^{-s})\Gamma(-s)\zeta(-s)\zeta(1-s)\ {\rm d}s\\
&=\frac{1}{2\pi i}\int^{\frac{3}{2}+i\infty}_{\frac{3}{2}-i\infty}2\left(\frac{2\pi^2}{\mu}\right)^{-s}(1-2^{-s})\Gamma(s)\zeta(s)\zeta(s+1)\ {\rm d}s\\
\end{align}
que es igual a $\displaystyle S\left(\frac{2\pi^2}{\mu}\right)$.
Funcional de la Ecuación de $S(\mu)$
Por el teorema de los residuos,
\begin{align}
\frac{1}{2\pi i}\int_C Q(s,\mu)\ {\rm d}s=\sum^1_{k=-1}{\rm Res}(Q(s,\mu);z=k)=\frac{\pi^2}{6\mu}-\ln{2}+\frac{\mu}{12}
\end{align}
Por lo tanto $S(\mu)$ satisface la ecuación funcional
$$S(\pi\mu\sqrt{2})+S\left(\frac{\pi\sqrt{2}}{\mu}\right)=\frac{\pi}{6\sqrt{2}}\left(\mu+\frac{1}{\mu}\right)-\ln{2}$$
Por ejemplo, esta ecuación se obtiene la identidad
$$\sum_{n\ge1}\frac{1}{n\sinh(n\pi\sqrt{2})}=\frac{\pi}{6\sqrt{2}}-\frac{\ln{2}}{2}$$
Conjeturó forma de $S(\mu)$
Estamos realmente interesados en la cantidad de $S\left(\pi\right)-\dfrac{1}{2}S(2\pi)$. Es probable que $S(\pi\mu\sqrt{2})$ es de la forma
$$S(\pi\sqrt{2}\mu)\stackrel{?}=\frac{\pi}{6\sqrt{2}}\mu-\frac{\ln{2}}{2}+f(\mu)\tag1$$
tal que $f(\mu)=-f(\mu^{-1})$, $f(0)$ explota, y $f\left(\dfrac{1}{\sqrt{2}}\right)\stackrel{?}=\dfrac{\ln{2}}{4}$. Si este es realmente el caso, entonces realmente llegamos
\begin{align}
\sum_{n\ \text{odd}}\frac{1}{n\sinh(n\pi)}
&=S(\pi)-\frac{1}{2}S(2\pi)\\\
&=\left(\frac{\pi}{12}-\frac{\ln{2}}{2}+\frac{\ln{2}}{4}\right)-\frac{1}{2}\left(\frac{\pi}{6}-\frac{\ln{2}}{2}-\frac{\ln{2}}{4}\right)\\
&=\frac{\ln{2}}{8}
\end{align}
lo que por supuesto implica que $\displaystyle\sum_{n\ge1}\frac{1}{n\sinh(n\pi)}=\frac{\pi}{12}-\frac{\ln{2}}{4}$. Numéricamente, esto parece ser cierto.
Por desgracia, no tengo idea de cómo probar o refutar $(1)$ e la $f(\mu)$.