Abstract: Based on three-dimensional progressive damage theory, a three-dimensional finite element model of multiple-bolt beveled-lap composite joints is built using the secondary development language of ANSYS software. Considering the non-linear contact behavior and the non-linear progressive damage process and using the Newton-Rapson nonlinear iterative algorithm, the failure process is simulated and the strength is analyzed. The correctness of this model is verified through its comparison with existing references. The results show that the failure extensions vary with the degree of bevel and the geometrical parameters. Failure types vary with stacking angles. To some extent, the reduction of bevel optimizes the load distribution, and the beginning and final failure loads increase. The connection efficiency is improved.