First-principles density functional theory ͑DFT͒ calculations have been used to obtain the constitutive relationships of pentaerythritol tetranitrate ͑PETN-I͒, a crystalline energetic material. The isotropic equation of state ͑EOS͒ for hydrostatic compression has been extended to include uniaxial compressions in the ͗100͘, ͗010͘, ͗001͘, ͗110͘, ͗101͘, ͗011͘, and ͗111͘ crystallographic directions up to a compression ratio of V / V 0 = 0.70. DFT predicts equilibrium properties such as lattice parameters and elastic constants, as well as the hydrostatic EOS, in agreement with available experimental data. Our results show a substantial anisotropy of various properties of PETN-I upon uniaxial compression. To characterize the anisotropic traits of PETN, different physical properties of the uniaxially compressed crystal such as the energy per atom, band gap, and stress tensor have been evaluated as a function of compression ratio. The maximum shear stresses were calculated and examined for a correlation with the anisotropy in shock-initiation sensitivity.