Doing the same as JJacquelin in his answer to your previous question, using $x=ct$ we have
$$I=\int \sqrt{b - \frac{c^6}{x^6-c^6}}\,dx=c \int \sqrt{b-\frac{1}{ t^6-1}}\,dt$$ that is to say $$I=\sqrt{b+1}\, c \,t
F_1\left(\frac{1}{6};\frac{1}{2},-\frac{1}{2};\frac{7}{6};t^6,\frac{b
t^6}{b+1}\right)$$ where appears the Appell hypergeometric function of two variables.
Integrating between $0$ and $1$ and doubling the result
$$J=\int_{-c}^c \sqrt{b - \frac{c^6}{x^6-c^6}}dx=2c\sqrt{b+1 } \frac{ \Gamma \left(\frac{7}{6}\right) \,\Gamma
\left(\frac{1}{2}\right)
}{\Gamma
\left(\frac{2}{3}\right)}\,_2F_1\left(-\frac{1}{2},\frac{1}{6};\frac{2}{3};\frac{b}{b+1}\right)$$
Edit
From a computational point of view, it is interesting to notice that
$$\, _2F_1\left(-\frac{1}{2},\frac{1}{6};\frac{2}{3};u\right)=1-\frac{u}{8}-\frac{7 u^2}{320}-\frac{91 u^3}{10240}-\frac{1729
u^4}{360448}-\frac{8645 u^5}{2883584}+O\left(u^6\right)$$ which would make the converge quite fast even for large values of $b$.