In this paper, a systematic approach for the dynamic modeling of complex planetary gear sets is presented. The approach relies upon a set of rules for building the system matrices and vectors of the system full dynamic model for any planetary gear set. A congruent state-space transformation is applied to obtain a reduced-order rigid model of the system, which allows for faster simulations. The behavior of the tangential forces accounting for the gears interactions is proven to be obtained from the reduced-order model. Furthermore, the kinematic relations of the considered planetary gear set are automatically generated when developing the reduced-order rigid model. As an example, two systems of interest in the vehicle industry are then modeled with the proposed approach and simulated in Matlab/Simulink: a Ravigneaux planetary gear set and a double-stage planetary gear set.