Tenemos
$$
\int_0^1 \frac{\ln(1+x)}{x^2+1}dx=
\left[\begin{array}{l}
t=\frac{1}{x} \\
dt=-\frac{1}{x^2}dx
\end{array}\right]=\int_1^\infty \frac{\ln(t+1)-\ln{t}}{1+t^2}dt=
\left[\begin{array}{l}
u=t-1 \\
du=dt
\end{array}\right] =
\int_0^\infty \frac{\ln(u+2)-\ln(u+1)}{1+(u+1)^2}du.
$$
Estoy bastante seguro de que esta integral puede resolverse con una cerradura de contorno en mucho de una manera similar a como esta integral. Si quieres mas detalles, solo un comentario, y voy a actualizar cuando tengo tiempo.
EDIT: se Continuó con el argumento. Por favor notificarme sobre errores tipográficos u otros errores.
Voy a mostrar cómo calcular
$$\int_0^\infty \frac{\ln(x+2)}{1+(x+1)^2}dx$$
(la otra es resuelto de forma análoga).
Ahora consideramos la integral
$$\int_C \frac{\log^2(z+2)}{1+(z+1)^2}dz,$$
donde $C$ es el ojo de la cerradura de contorno anterior. Nota la opción de logaritmo! Es imperativo que tomemos $0<\arg{z}<2\pi$ de los logaritmos (o $2\pi n<\arg{z}<2\pi (n+1)$ algunos $n\in\mathbb{Z}$).
El integrando tiene polos en $z=-1\pm i$, y por Cauchy de la integral de la fórmula, obtenemos
$$
\int_C \frac{\log^2(z+2)}{1+(z+1)^2}dz=
\int_C \frac{\log^2(z+2)}{(z+1+i)(z+1-i)}dz=
2\pi i\left(\frac{\log^2(1-i)}{-2i}+
\frac{\log^2(1+i)}{2}\right)=
\pi(\log^2(1+i)-\log^2(1-i))=
\pi\left(\left(\ln\sqrt{2}+\frac{i\pi}{4}\right)^2-\left(\ln\sqrt{2}+\frac{7i\pi}{4}\right)^2\right)=
\frac{6\pi^3-3i\pi^2\ln{2}}{2}
.$$
La integral a lo largo de $\gamma$ va a cero a medida que el radio interior va a cero, y la integral a lo largo de $\Gamma$ va a ceros como $R\rightarrow\infty$ (esto puede ser demostrado que el uso de la ML-desigualdad).
Por lo tanto, en el límite sólo se consiguen las contribuciones de las partes rectas, por lo que
$$
\int_0^R\frac{\log^2(x+i\varepsilon+1)}{(1+x+i\varepsilon)^2+1}dx+
\int_R^0\frac{\log^2(x-i\varepsilon+1)}{(1+x-i\varepsilon)^2+1}dx\rightarrow \frac{6\pi^3-3i\pi^2\ln{2}}{2}.
$$
Pero también tenemos
$$
\int_0^R\frac{\log^2(x+i\varepsilon+2)}{(1+x+i\varepsilon)^2+1}dx+
\int_R^0\frac{\log^2(x-i\varepsilon+2)}{(1+x-i\varepsilon)^2+1}dx\rightarrow
\int_0^\infty\frac{\ln^2(x+2)-\left(\ln(x+2)+2i\pi\right)^2}{(1+x)^2+1}dx=
\int_0^\infty\frac{4\pi^2-4i\pi\ln(x+2)}{(1+x)^2+1}dx,
$$
así
$$\int_0^\infty\frac{4\pi^2-4i\pi\ln(x+2)}{(1+x)^2+1}dx=\frac{6\pi^3-3i\pi^2\ln{2}}{2}.$$
Considerando que la parte imaginaria, obtenemos que
$$\int_0^\infty \frac{\ln(x+2)}{(x+1)^2+1}dx=\frac{3\pi\ln{2}}{8}.$$
Es posible calcular
$$\int_0^\infty \frac{\ln(x+1)}{(x+1)^2+1}dx$$
de forma análoga, y, a continuación, tomar la diferencia en el fin de encontrar el valor de la integral original.