Articulated solid bodies are shape-changing bodies made of rigid solids linked together by means of holonomic constraints prescribed as functions of time. In this paper we study the locomotion of such swimming devices in an ideal fluid. Our study ranges over a wide class of problems: any number of immersed bodies are involved (without being hydrodynamically decoupled), the system fluid-bodies can be partially of totally confined and circulation, buoyancy and collisions between bodies are taken into account. We determine the Euler-Lagrange equation governing the dynamic of the system, study its well-posedness and describe a numerical scheme used in a Matlab toolbox (Biohydrodynamics Toolbox) that has been designed to realize easily related numerical simulations.