$$\boxed{\displaystyle\int_0^1\log\log\left(\frac1x+\sqrt{\frac1{x^2}-1}\right)\mathrm dx=-\gamma-2\ln\frac{2\Gamma(3/4)}{\Gamma(1/4)}}\tag{$\heartsuit$}$$
Derivation:
After the change of variables $x=\frac{1}{\cosh u}$ the integral becomes
$$\int_0^{\infty}\ln u \frac{\sinh u}{\cosh^2 u}du,$$
as was noticed above by Eric. We would like to integrate by parts to kill the logarithm but we get two divergent pieces. To go around this, let us consider another integral,
$$I(s)=\int_0^{\infty}u^s \frac{\sinh u}{\cosh^2 u}du,$$
with $s>0$. The integral we actually want to compute is equal to $I'(0)$, which will be later obtained in the limit.
Indeed, integrating once by parts one finds that
\begin{align}
I(s)&=s\int_0^{\infty}\frac{u^{s-1}du}{\cosh u}=s\cdot 2^{1-2 s}\Gamma(s)\left[\zeta\left(s,\frac14\right)-\zeta\left(s,\frac34\right)\right]=\\ &=2^{1-2 s}\Gamma(s+1)\left[\zeta\left(s,\frac14\right)-\zeta\left(s,\frac34\right)\right],
\end{align}
where $\zeta(s,a)=\sum_{n=0}^{\infty}(n+a)^{-s}$ denotes Hurwitz zeta function (in the way we have used its integral representaion (5) from here).
Now to get ($\heartsuit$), it suffices to use
\begin{align}
&\frac{\partial}{\partial s}\left[2^{1-2 s}\Gamma(s+1)\right]_{s=0}=-2\gamma-4\ln 2,\\
&\zeta\left(0,\frac14\right)-\zeta\left(0,\frac34\right)=\frac12,
\\
&\frac{\partial}{\partial s}\left[\zeta\left(s,\frac14\right)-\zeta\left(s,\frac34\right)\right]_{s=0}=-\ln\frac{\Gamma(\frac34)}{\Gamma(\frac14)}.
\end{align}
[See formulas (10) and (16) on the same page].
A closed form indeed exists for this integral:
$$\int_0^\infty\frac{\log\left(1+\frac{\pi^2}{4\,x}\right)}{e^{\sqrt{x}}-1}\mathrm dx=\pi^2\left(\log\frac{\pi\,A^9\sqrt{2}}{\Gamma(\frac{1}{4})^2}-\frac{9}{8}\right)+2\,\pi\,C,$$
where $A$ is the Glaisher-Kinkelin constant and $C$ is the Catalan constant.
A more general result:
$$\int_0^\infty\frac{\log\left(1+\frac{a}{x}\right)}{e^{\sqrt{x}}-1}\mathrm dx=8\,\pi^2\psi^{(-2)}\left(\frac{\sqrt{a}}{2\pi}\right)-\frac{a}2\left(1+\log\frac{4\pi^2}{a}\right)-2\pi\sqrt{a}\left(1+2\log\Gamma\left(\frac{\sqrt{a}}{2\pi}\right)\right),$$
where $\psi^{(-2)}(z)$ is the generalized polygamma function.
The proof is based on Binet's second formula, but I still need to sort out some details.
Best Answer
Notice, we have $$\int_{0}^{\frac{e-1}{e}}\frac{x(2-x)}{1-x}\frac{\log\left(\log\left(1+\frac{x^2}{2-2x}\right)\right)}{2-2x+x^2}dx$$ $$=\int_{0}^{\frac{e-1}{e}}\frac{x(2-x)}{1-x}\frac{\log\left(\log\left(\frac{2-2x+x^2}{2-2x}\right)\right)}{2-2x+x^2}dx$$
Let, $$\log\left(\frac{2-2x+x^2}{2-2x}\right)=u$$ $$\implies \frac{d}{dx}\left(\log\left(\frac{2-2x+x^2}{2-2x}\right)\right)=\frac{d}{dx}(u)$$ $$\frac{1}{\left(\frac{2-2x+x^2}{2-2x}\right)}\cdot \left(\frac{(2-2x)(-2+2x)-(2-2x+x^2)(-2)}{(2-2x)^2} \right)=\frac{du}{dx}$$ $$\left(\frac{2-2x}{2-2x+x^2}\right)\cdot \left(\frac{2x(2-x)}{(2-2x)^2} \right)=\frac{du}{dx}$$ $$\frac{x(2-x)}{(1-x)}\frac{1}{(2-2x+x^2)}dx=du$$ Now, we have $$\int_{0}^{\log\left(\frac{e^2+1}{2e}\right)}\log(u)du$$
$$=\left[u\log(u)-u\right]_{0}^{\log\left(\frac{e^2+1}{2e}\right)}$$ $$=\left[u\log\left(\frac{u}{e}\right)\right]_{0}^{\log\left(\frac{e^2+1}{2e}\right)}$$
$$=\log\left(\frac{e^2+1}{2e}\right)\cdot\log\left(\frac{1}{e}\log\left(\frac{e^2+1}{2e}\right)\right)-\lim_{u\to 0}u\log\left(\frac{u}{e}\right)$$
$$=\log\left(\frac{e^2+1}{2e}\right)\cdot\log\left(\frac{1}{e}\log\left(\frac{e^2+1}{2e}\right)\right)-0$$
Hence, we get
$$\bbox[5px, border:2px solid #C0A000]{\color{red}{\int_{0}^{\frac{e-1}{e}}\frac{x(2-x)}{1-x}\frac{\log\left(\log\left(1+\frac{x^2}{2-2x}\right)\right)}{2-2x+x^2}dx}=\color{blue}{\log\left(\frac{e^2+1}{2e}\right)\cdot \log\left(\frac{1}{e}\log\left(\frac{e^2+1}{2e}\right)\right)}}$$