This study investigates the numerical simulation and theoretical modeling of longitudinal compressive failure in fiber reinforced composite materials. Firstly the numerical simulation of longitudinal compressive failure is conducted. The simulated results show that at one moment of the loading, the localized deformation catastrophically appears in the material, and in this initiation of the localized deformation, the reduction of tangent shear stiffness plays an important role. Secondly the theoretical modeling of longitudinal compressive failure is implemented. A set of mathematical equations is obtained for the deformation of composite materials, and the mathematical solution of the equations is considered. There exists a state where arbitrariness appears in the solution of equations expressing deformation of composite materials, and it is indicated that the onset of arbitrariness in solution of equations expressing deformation of composite materials is closely related with the initiation of longitudinal compressive failure, and also related with the initiation of narrow localized band in the materials. Finally the numerical simulation is conducted for compressive failure in quasi-isotropic laminate. The localized deformation also appears in the laminate, and the simulated deformation of the material agrees with the microscope picture of the experimental result.