The coupled dynamic response of two cooperating robots handling two flexible payloads is treated using a new algorithm. In this algorithm, the dynamic equations describing the system are obtained using Lagrange's method for the rigid robot links and the finite element method for the flexible payloads. The contact between the flexible payloads is modelled using the penalty function method and a contact search algorithm is employed to identify the contact region.