A method for the parametric generation of the transversal cross-section microstructure model of unidirectional long-fiber reinforced composite (LFRC) is presented in this paper. Meanwhile, both the random distribution of the fibers and high fiber volume fraction are considered in the algorithm. The fiber distribution in the cross-section is generated through random movements of the fibers from their initial regular square arrangement. Furthermore, cohesive zone element is introduced into modeling the interphase between the fiber and the matrix. All these processes are carried out by the secondary development of the finite element codes (ABAQUS) via Python language programming. Based on the model generated, micromechanical finite element analysis (FEA) is performed to predict the damage initiation and subsequent evolution of the composites. The results show that this technique is capable of capturing the random distribution nature of these composites even for high fiber volume fraction. Moreover, the results prove that a good agreement with the experimental results is found.