Limit analysis for stone structures has become a widely used method due to (i) the large range of results it can give while requiring very few mechanical assumptions, and (ii) the resolution efficiency, with conic optimization resolution schemes. The aim of this research consists in an interactive tool, Stabilithos, implemented in the Rhinoceros modeling software, allowing the exploration of different types of limit analysis on any 2D or 3D blocks structures. Various choices of fracture mechanical hypotheses for the joints between stone blocks are possible: Heyman hypotheses, but also more complex criteria, implemented without any simplifying assumptions for 2D and 3D convex shaped joints, thanks to conic optimization: finite compressive strength and/or finite sliding resistance with Coulomb criteria. The simulation results are given in graphical form directly in Rhinoceros: the internal forces acting on each interface between blocks, the normal constraints if the compression constraint has a finite strength, and the animation of kinematic failure mechanism. Thanks to the interface, the following limit analysis loading modes in any given direction are available: maximizing one or several forces applied on the structure, maximizing a support reaction, or maximizing a force field proportional to block’s self-weight and applied to the block’s gravity center. Additionally, it is possible to automatically determine the static and geometric safety factors of any 2D or 3D stone structure in order to ensure and evaluate its stability under dead loads. The tool’s possibilities are evaluated with the application case of Notre-Dame cathedral in Paris.