This study, framed within the context of the lower bound theorem of limit analysis, aims to assess the anti-funicular equilibrium of masonry arches and domes using a computational approach based on the constrained force density method. In contrast to the commonly adopted classical Heyman’s assumptions, the approach proposed here considers the effects of finite compressive strength in the material. Assuming a fixed plan projection for a network with independent sets of branches, a suitable set of local constraints is enforced at each joint to account for the limit bending moment resulting from the material’s assumptions, including limited compressive strength and zero tensile strength. Additionally, the stereotomy of the voussoirs is considered by assigning a geometric law to the joint inclination. The collapse load is determined by formulating a multi-constrained maximization problem. The method is validated using a modern version of the semi-analytical Durand-Claye’s method. For benchmark case studies, such as symmetric masonry arches and domes with specific stereotomies subjected to axi-symmetrical load conditions, the set of statically admissible solutions compatible with equilibrium and strength requirements is graphically determined in terms of the horizontal thrust and its eccentricity at the crown, examining the shape of the stability area. Assuming an infinite value for the friction coefficient, the collapse condition is reached when the stability area shrinks to a single point. The results obtained from both of these methods are in excellent agreement. The influence of compressive strength on the bearing capacity of the structures is also discussed.