We perform a systematic comparison of the finite-temperature structure and properties of four bulk semiconductors (PbS, PbTe, ZnS, and ZnTe) predicted by eight popular exchange-correlation functionals from quasi-harmonic lattice-dynamics calculations. The performance of the functionals in reproducing the temperature dependence of a number of material properties, including lattice parameters, thermal-expansion coefficients, bulk moduli, heat capacities, and phonon frequencies, is evaluated quantitatively against available experimental data. We find that the phenomenological over- and under-binding characteristics of the local-density approximation and the PW91 and Perdew-Burke-Enzerhof (PBE) generalised-gradient approximation (GGA) functionals, respectively, are exaggerated at finite temperature, whereas the PBEsol GGA shows good general performance across all four systems. The Tao-Perdew-Staroverov-Scuseria (TPSS) and revTPSS meta-GGAs provide relatively small improvements over PBE, with the latter being better suited to calculating structural and dynamical properties, but both are considerably more computationally demanding than the simpler GGAs. The dispersion-corrected PBE-D2 and PBE-D3 functionals perform well in describing the lattice dynamics of the zinc chalcogenides, whereas the lead chalcogenides appear to be challenging for these functionals. These findings show that quasi-harmonic calculations with a suitable functional can predict finite-temperature structure and properties with useful accuracy, and that this technique can serve as a means of evaluating the performance of new functionals in the future.