Fluids of hard bent-core molecules have been studied using theory and computer simulation. The molecules are composed of two hard spherocylinders, with length-to-breadth ratio L/D, joined by their ends at an angle 180 degrees - gamma. For L/D = 2 and gamma = 0,10,20 degrees, the simulations show isotropic, nematic, smectic, and solid phases. For L/D = 2 and gamma = 30 degrees, only isotropic, nematic, and solid phases are in evidence, which suggests that there is a nematic-smectic-solid triple point at an angle in the range 20 degrees < gamma < 30 degrees. In all of the orientationally ordered fluid phases the order is purely uniaxial. For gamma = 10 degrees and 20 degrees, at the studied densities, the solid is also uniaxially ordered, whilst for gamma = 30 degrees the solid layers are biaxially ordered. For L/D = 2 and gamma = 60 degrees and 90 degrees we find no spontaneous orientational ordering. This is shown to be due to the interlocking of dimer pairs which precludes alignment. We find similar results for L/D = 9.5 and gamma = 72 degrees, where an isotropic-biaxial nematic transition is predicted by Onsager theory. Simulations in the biaxial nematic phase show it to be at least mechanically stable with respect to the isotropic phase, however. We have compared the quasi-exact simulation results in the isotropic phase with the predicted equations of state from three theories: the virial expansion containing the second and third virial coefficients; the Parsons-Lee equation of state; an application of Wertheim's theory of associating fluids in the limit of infinite attractive association energy. For all of the molecule elongations and geometries we have simulated, the Wertheim theory proved to be the most accurate. Interestingly, the isotropic equation of state is virtually independent of the dimer bond angle-a feature that is also reflected in the lack of variation with angle of the calculated second and third virial coefficients.