The primary objective of the work presented in this paper was to develop a continuum thermoelastic-damage model for carbon fibre reinforced plastic (CFRP) materials, capable of modelling high rate deformation typical for ballistic impact loading. The constitutive model is capable of predicting formation and propagation of the shock waves in orthotropic materials and in addition can simulate damage initiation, evolution and failure. A key feature of the constitutive model is the decomposition of material volumetric and shear response. Material response under compression in this model is defined in terms of Mie Gruneisen equation of state (EOS) and the decomposition of stress tensor proposed in . In order to take into account the orthotropy of the CFRP materials of interest, damage in this constitutive model is represented by a second order damage tensor ω, which is incorporated in the stiffness tensor by using energy equivalence principle, see for instance . Validation of the numerical model, implemented in LLNL DYNA3D  was done by the comparison of the numerical results to the experimental data obtained in the high velocity sphere impact tests published in . The numerical results for the extent of damage were within 8% with the corresponding experimental data.