Based on the force density method, an iterative procedure is developed that enables the generation of spatial kinematic pin-jointed structures that are in equilibrium close to a given input geometry, while satisfying additional constraints on both geometry and forces. This method forms the core of an interactive form-finding process that consists of alternating steps of modeling and computational optimization. In each modeling step, the user is able to modify the structure's geometry, topology, external forces, and constraints. In each optimization step, equilibrium is re-established while respecting the user-defined constraints. A prototype has been implemented within an existing CAD software package, and three examples illustrate the use of the presented method, ranging from a playful exploration of surprising shapes to the rationalization of structural geometry. The method allows to intuitively explore the formal freedom of spatial equilibrium shapes with mixed compression and tension forces, within hard, user-defined constraints. In conclusion, it is claimed that by providing interactive equilibrium modeling methods, the design of new, surprising spatial forms with efficient structural behavior is facilitated.