The most widespread numerical simulation method for structural response is undoubtedly the Finite Element Method (FEM). However, despite being powerful for modelling continuous structures, it is not fit for handling strong discontinuities. The Discrete Element Method (DEM) simulates interactions between rigid or deformable elements in contact, hence explicitly capturing the discontinuous nature of the structural response, particularly when subjected to extreme loadings. Nevertheless, it requires a time-stepping algorithm even for solving static or buckling problems. The Hybrid Discrete-Finite Element Method, shortly HybriDFEM, was recently introduced in the context of modelling one-dimensional beam-like members. Those members are divided along their longitudinal axis in a series of rectangular rigid blocks, and the deformation is concentrated at the interfaces between adjacent blocks, modelled as distributed nonlinear multidirectional springs. The method, developed within a FEM-like setting, allows for hybridisation with other finite elements (e.g., beam elements). Next to its ability to explicitly model pre-existing discontinuities along the member (e.g. masonry stereotomy), the method can be used for modelling continuous members with satisfactory accuracy by appropriately scaling the interface springs. As such, the HybriDFEM’s formulation can accommodate hybrid discrete-continuous systems. In this paper, the HybriDFEM formulation is extended to 2D, with rectangular blocks in contact on all four faces. First, the algorithm to detect blocks in contact will be explained. Second, specific characteristics of the HybriDFEM applied to masonry modelling are presented. Then, the method is benchmarked against a two-dimensional problem from the literature where the in-plane capacity of walls made masonry blocks is investigated.