This letter proposes a systematic modular procedure for the dynamic modeling of branched robots comprising several subsystems, each of which being composed of multiple rigid bodies. Furthermore, the proposed strategy is applicable even if some subsystems are regarded as black boxes, requiring only the twists and wrenches at the connection points between different subsystems. To help in the model composition, we also propose a graph representation that encodes the propagation of twists and wrenches between the subsystems. Numerical results show that the proposed formalism is as accurate as a state-of-the-art library for robotic dynamic modeling.