jats:titleAbstract</jats:title> jats:pWe demonstrate how a geometrically exact formulation of discrete slender beams can be generalized for the efficient simulation of complex networks of flexible beams by introducing rigid connections through special junction elements. The numerical framework, which is based on discrete differential geometry of framed curves in a time-discrete setting for time- and history-dependent constitutive models, is applicable to elastic and inelastic beams undergoing large rotations with and without natural curvature and actuation. Especially, the latter two aspects make our approach a versatile and efficient alternative to higher-dimensional finite element techniques frequently used, e.g., for the simulation of active, shape-morphing, and reconfigurable structures, as demonstrated by a suite of examples.</jats:p>