This work aims to provide a computationally effective and distributed trajectory planner at the intersection of physics-based and data-driven techniques for the collaborative locomotion of holonomically constrained quadrupedal robots that can account for and attenuate interaction forces between subsystems. More specifically, this work lays the foundation for using an interconnected single rigid body model in a predictive control framework such that interaction forces can be utilized at the planning layer, wherein these forces are parameterized via a behavioral systems approach. Furthermore, the proposed trajectory planner is distributed such that each agent can locally plan for its own trajectory subject to coupling dynamics, resulting in a much more computationally efficient method for real-time planning. The optimal trajectory obtained by the planner is then provided to a full-order nonlinear whole-body controller for tracking at the low level. The efficacy and robustness of the proposed approach are verified both in simulation and on hardware subject to various disturbances, payloads, and uneven terrains.