European cultural and architectural heritage is predominantly characterized by masonry structures. Most of these consist of buildings, bridges, ancient aqueducts, churches built with different type of masonry materials. As masonry is, in any case, characterized by a low tensile and, hence, bending strength, many structural elements are made up of curved members, such as arches, domes and vaults. Over the years, these structures, in some cases thousands of years old, have suffered structural subsidence, damage due to changes in intended use or have been subjected to hazard events, such as seismic events. The assessment of their safety and the study of structural improvement actions are considered as both a national and international priority to guarantee their historical and artistic preservation over time. Recently, different computational procedures have been suggested thanks to a large amount of scientific studies [1-3]. In particular, the Limit Analysis has appeared as a reliable and suitable approach for the prediction of the load-carrying capacity of masonry curved members. Likewise, the development of Finite Element modelling (FE) and Discrete Macro-Element Methods (DMEM) have revealed their efficiency in the simulation of the complex structural response of masonry. Nevertheless, these latter could not allow to consider the proper curvature of the masonry members and in some cases have high computational costs. The purpose of this research is to present a novel Lagrangian-based and NURBS-based beam elements able to reproduce the no-tension mechanical response of curved arches [4] and axisymmetric domes with any shape. Two different procedures are suggested: i) a finite element (FE) implemented considering a novel Lagrangian approximation field, ii) a finite element developed using the Isogeometric analysis (IGA) [5-6]. For both the numerical approaches the no-tension property of the masonry material is taken into account in the definition of the constitutive law. The theoretical formulation is described and the adopted numerical procedures are illustrated and discussed, along with the outcomes display, sensitivity analyses and comparisons between the two numerical suggested approaches.