Esta integral puede evaluarse utilizando los residuos de la teoría. Considere la integral
$$\oint_C \frac{dz}{(1+z^2) \cosh{a z}}$$
donde $C$ es un semicírculo de radio $R$ en la mitad superior del plano. Como $R \to \infty$, la integral sobre el semicírculo se desvanece, y nos quedamos con el original de la integral igualando $i 2 \pi$ vez que la suma de los residuos de los polos de el integrando dentro de $C$. En este caso, los polos dentro de $C$ se encuentran en $z_n = i (n+1/2) \pi/a$ todos los $n \in \mathbb{N} \cup \{0\}$ y a las $z_+ = i$. La evaluación de los residuos en estas polos (lo que puede lograrse cuando el integrando es de la forma $p(z)/q(z)$ usando la fórmula de $p(z_0)/q'(z_0)$ para un poste de $z=z_0$), encontramos que
$$\int_{-\infty}^{\infty} \frac{dx}{(1+x^2) \cosh{a x}} = \frac{\pi}{\cos{a}} - \frac{2 \pi}{a} \sum_{n=0}^{\infty} \frac{(-1)^n}{(n+1/2)^2 \pi^2/a^2 - 1}$$
La suma por desgracia, toma la forma de un par de Lerch transcendents
$$\begin{align}\frac{2 \pi}{a} \sum_{n=0}^{\infty} \frac{(-1)^n}{(n+1/2)^2 \pi^2/a^2 - 1} &= \frac{\pi}{a} \sum_{n=0}^{\infty} (-1)^n \left (\frac{1}{(n+1/2)\pi/a-1}-\frac{1}{(n+1/2)\pi/a+1} \right)\\&= \sum_{n=0}^{\infty} (-1)^n \left (\frac{1}{n+\frac12-\frac{a}{\pi}}-\frac{1}{n+\frac12+\frac{a}{\pi}} \right) \\ &= \Phi\left(-1,1,\frac12-\frac{a}{\pi}\right)-\Phi\left(-1,1,\frac12+\frac{a}{\pi}\right)\end{align}$$
Por lo tanto
$$\int_0^{\infty} \frac{dx}{(1+x^2) \cosh{a x}} = \frac{\pi}{2\cos{a}} - \frac12 \left[\Phi\left(-1,1,\frac12-\frac{a}{\pi}\right)-\Phi\left(-1,1,\frac12+\frac{a}{\pi}\right) \right ]$$
Cabe señalar que $a \ne (k+1/2) \pi$ algunos $k \in \mathbb{Z}$.
ANEXO
Debo señalar que, en respuesta a @GrahamHesketh de la consulta, el resultado anterior puede ser demostrado ser igual a la diferencia entre dos digamma funciones de la siguiente manera:
$$\int_0^{\infty} \frac{dx}{(1+x^2) \cosh{a x}} = \frac12 \left [ \psi\left(\frac{3}{4}+\frac{a}{2 \pi} \right)-\psi\left(\frac{1}{4}+\frac{a}{2 \pi} \right) \right ]$$
Esto se puede lograr tomando nota de que
$$\frac{\pi}{\cos{a}} = \sum_{n=-\infty}^{\infty} (-1)^n \frac{1}{n+\frac12-\frac{a}{\pi}} = \sum_{n=0}^{\infty} (-1)^n \left (\frac{1}{n+\frac12-\frac{a}{\pi}}+\frac{1}{n+\frac12+\frac{a}{\pi}} \right) $$
$$\psi\left(\frac{1}{4}+\frac{a}{2 \pi} \right) = \sum_{n=0}^{\infty}\left (\frac{1}{n+1}- \frac{1}{n+\frac12 \left (\frac12+\frac{a}{\pi}\right)}\right )$$
$$\psi\left(\frac{3}{4}+\frac{a}{2 \pi} \right) = \sum_{n=0}^{\infty} \left ( \frac{1}{n+1}-\frac{1}{n+1-\frac12 \left (\frac12-\frac{a}{\pi}\right)}\right )$$
Para establecer la igualdad, tenga en cuenta que el resultado que he publicado anteriormente se reduce a
$$\frac{\pi}{\cos{a}} - \sum_{n=0}^{\infty} (-1)^n \left (\frac{1}{n+\frac12-\frac{a}{\pi}}-\frac{1}{n+\frac12+\frac{a}{\pi}} \right) = 2 \sum_{n=0}^{\infty} \frac{(-1)^n}{n+\frac12+\frac{a}{\pi}}$$
La igualdad entre dicha suma y la diferencia entre los dos $\psi$ términos se establece mediante la comparación de los sumandos término por término.