The HybriDFEM method, short for Hybrid Discrete-Finite Element Method, combines both discrete and finite element approaches in a single numerical model: the method adopts a discrete representation of the structure, but the formulation is designed to integrate continuous parts that can be simulated by the Finite Element Method, allowing hybrid numerical mock-ups to be built. The scope of application of the method is expanded from its original development for uniaxially-discretized (1D) structures to modeling biaxially-discretized (2D) structures and systems of beams connected through rigid-node connections. The possibility to integrate finite elements within HybriDFEM in 2D is initially formalized. A two-step contact detection algorithm, in which the interface detection is preceded by a preliminary rough detection, is then presented. Finally, different approaches to modeling contact are introduced, depending on whether it is meant to reflect the behavior of a continuous material, flexible interfaces, or point-wise contact. These new capabilities of the 2D HybriDFEM method are validated on a series of selected examples including solutions from analytical models, classical finite elements, and limit analysis; among others, the HybriDFEM method is used to evaluate the axial and shear stress distribution in a linear elastic beam with negligible error relative to analytical solutions, and to predict the collapse load of in-plane loaded masonry frames with an error below 0.05% compared to solutions from limit analysis. The adequacy of the method to enhance discrete simulations by integrating finite elements is illustrated in a conclusive example via the pushover analysis of a flexible masonry frame.