The shock Hugoniot of boron carbide, from 0 to 80 GPa, has been obtained using first principles quantum mechanics (density functional theory) and molecular dynamics simulation. The Hugoniot for six different structures which vary by structure or stoichiometry were computed and compared to experimental data. The effect of stoichiometry, and structural variation within a given stoichiometry, are shown to have marked effects on the shock properties with some compositions displaying bilinear behavior in the computed shock velocity‐particle velocity profiles while others show a continuous Hugoniot curve with no evidence of a phase transition over the pressure range considered in this work. Two structures, B12(CBC) and B11Cp(CCB), have predicted phase transition pressures lying within the 40–50 GPa range suggested experimentally. It is shown that the phase transition is driven by deformation of the 3‐atom chain within the boron carbide crystal structure which induces a discontinuous volume change at the critical shock pressure. The effect of defects, in the form of chain vacancies, on the shock response is presented and the ability of shear to significantly lower the phase transition pressure, in accord with experimental observation, is demonstrated.