The full-potential linear muffin-tin orbital method is used within the local density approximation and generalized gradient approximation ͑GGA͒ to calculate the total energy and equilibrium lattice properties for the observed phases of Zr. The temperature dependences of the free energy, specific volume, bulk modulus, Debye temperature, and Grüneisen constant are found for these structures within the Debye model. For most quantities, a good quantitative agreement with experiment is obtained. The P-T phase diagram constructed from the calculated thermodynamical Gibbs potentials within the GGA fits well the available room-temperature data on the ␣→ and → transitions. At ambient pressure, we get T →␣ ϭ1193 K, which is close to the observed value. ͓S0163-1829͑98͒02514-4͔