Given the fact that computational effort and memory requirements are considerable when applying the quasi-bond method (QBM) to structural analyses, this paper introduces a hybrid numerical method combining QBM and the classical finite element method for modeling and simulating deformation and crack evolution in quasi-brittle materials. The proposed approach, referred to as QBFEM, utilizes linear elements for spatial discretization of structures. Cracked material regions and potential damage zones are designated as quasi-bond analysis domains where continuous-discontinuous solutions are performed, while non-damaged zones are solved using the finite element method without involving any coupling regions along the QBM-FEM transition zones. Subsequently, several tests of stiffness matrix assembling by both QBM and QBFEM for different meshes were conducted to demonstrate the efficiency of QBFEM. Finally, three numerical examples-plate with a single hole under tension, single crack plate under shear, and undamaged gravity dam under hydrostatic pressure-are provided for validation, which demonstrate that QBFEM can effectively capture the initiation and propagation of Mode I fracture in solid materials.