Given $f(x)=\frac{x^2}{\sinh ^2 x}$
I computed the Laplace transform $\mathscr{L}\left(f(x)\right)(s)=\frac{1}{4} \left(4 \psi ^{(1)}\left(\frac{s}{2}\right)+s \psi ^{(2)}\left(\frac{s}{2}\right)\right)$
Then $$\int_0^{\infty } \frac{x^2}{\sinh ^2 x} \, dx=\underset{s\to 0}{\text{lim}}\mathscr{L}\left(f(x)\right)(s)=\frac{\pi ^2}{6}=\zeta(2)$$ which is $$\sum _{n=1}^{\infty } \frac{1}{n^2}=\zeta(2)$$
Do you think there is a deeper link between these two results or is it just a coincidence?