Efficient and accurate calculation of stress overload and stress redistribution near the fiber breakpoint is the basis for simulating the tensile damage evolution of unidirectional fiber reinforced plastic (FRP), and it is also the premise for predicting the tensile-tensile fatigue life. In this study, an efficient finite element (FE) model with a large number of fibers was developed based on the FE method to calculate stress concentration factor (SCF) near multiple breakpoints, and the calculation time only takes 3 min on average. The FE model is used to efficiently calculate multiple breakpoints and the stress distribution near the breakpoint considering the interface debonding. Based on the efficient FE model, an analytical model is developed to calculate the stress overload superposition around multiple breakpoints. In this analytical model, the shear stress distribution near the broken fiber elements is established based on the experiment, and the stress redistribution mechanism near the breakpoint in FRP with random distribution of fibers is established. Moreover, the calculation method of stress overload near multiple break points is established, and the FE method is verified.