According to the increase in using of polymer composite materials in industry, exact characterization of the material properties in different loading conditions, including high velocity impact is very important. In this article, the high velocity impact on non-hybrid and hybrid polymer composite panel has been investigated. After introducing the impact load types and theories used for damage evaluation and comparison of them, the matzenmiler theory has been employed to predict damage. In this theory, in order to predict more accurate damage evaluation, the shear stress is considered nonlinear in elastic region. Afterwards, in order to use the theory in the finite element modeling, a VUMAT subroutine is implementated in the Abaqus software. Finally, the obtained results from the present modeling and analysis, are compared with available experimental data for non-hybrid and hybrid composite panels. The good agreement between the theoretical results and experimental data, introduces the ability of the applied model and provided subroutine