Let $g(a)=\frac{a}{1-a}, f(a)=\log(1-a)$ and $h(a)=\log^2(1-a)$. Let $\square^m$ denote a $m$-dimensional unit hypercube.
The following is an application or reformulation of the above infinite series evaluation. $$\int_0^1\int_0^1\frac{h(\frac12xz)-xh(\frac12z)+xzh(\frac12)-zh(\frac12x)}{x(1-x)z(1-z)}dxdz= \frac{\log^42}{12}-\frac{\pi^4}{144}+\zeta(2)\log^22+\zeta(3)\log2.$$ Proof outline: Since $\frac1k=\int_0^1x^{k-1}dx$, we may write \begin{align} \sum_{n\geq1}\frac{t^nH_n^3}{n+1}&=\sum_{n\geq1}t^n\int_0^1w^ndw\sum_{i,j,k=1}^n\int_{\square^3}x^{i-1}y^{j-1}z^{k-1}dxdydz \\ &=\int_{\square^4}dxdydzdw\sum_{n\geq1}(tw)^n\sum_{i,j,k=1}^nx^{i-1}y^{j-1}z^{k-1}\\ &=\int_{\square^4}\frac{dxdydzdw}{(1-x)(1-y)(1-z)}\sum_{n\geq1}(tw)^n(1-x^n)(1-y^n)(1-z^n)\\ &=\int_{\square^4}dxdydzdq\frac{g(tw)-g(twy)-g(twx)+g(twxy)-g(twz)+g(twyz)+g(twxz)-g(twxyz)}{(1-x)(1-y)(1-z)}. \end{align} Next, make a repeated use of the integral evaluations $$\int_0^1\frac{g(tw)-g(twy)}{1-y}dy=-\frac{f(tw)}{1-tw}$$ and $$\int_0^1\left(\frac{f(twx)}{1-twx}-\frac{f(tw)}{1-tw}\right)= \frac{xh(t)-h(tx)}{2tx}$$$$\int_0^1\left(\frac{f(twx)}{1-twx}-\frac{f(tw)}{1-tw}\right)\frac{dw}{1-w}= \frac{xh(t)-h(tx)}{2tx}$$ to arrival at the double integral upon replacing $t=\frac12$. The rest follows from the series evaluation algorithm that Julian Rosen alluded to. $\square$
Remark. This approach actually could prove (directly) the easier-looking $$\sum_{n\geq1}\frac{H_n^2}{(n+1)2^n}=\frac13\log^32+\zeta(2)\log2-\frac12\zeta(3).$$