This study presents an approach based on the merger of (local) classical continuum theory and (non-local) peridynamic theory to predict failure simulations in bolted composite lap joints. The global modeling is performed using the semi-analytical method while the PD theory is employed for submodeling and failure prediction. The results from the global analysis identifies the most critical bolts for subsequent analysis with PD approach as a submodel. Displacements from the global model are applied to the submodel as boundary conditions, and the initiation and propagation of damage is predicted based on the PD theory. The simulations capture the dominant failure modes of bearing failure, matrix cracking, fiber splitting and delamination in the region where the bolt is in contact with the laminate. These are consistent with common failure modes around the bolt hole observed in previous experimental investigations concerning failure modes occurring around the bolt hole.