How may I know how thermal expansion happens as a function of temperature, with AIMD?

I wanted to know how the volume changes with the temperature, with respect to the changes in lattice parameter for gamma Uranium. But when I use ISIF =3, it will cause changes in the symmetry. I wanted to keep the cubic symmetry of the gamma Uranium.