In this study, a two-dimensional (2D), steady-state, discrete dynamic model of a double-planet planetary gearset is proposed. The dynamic model is generalized such that it can consist of number of planet branches and can operate under any operating conditions (load and speed). The contact between each external to external and external to internal gear pair is modeled to obtain stiffnesses and mesh displacement excitations using a generalized load distribution model. The natural modes are computed by solving the corresponding eigenvalue problem. The forced vibration response to gear mesh excitations is obtained by applying the modal summation technique. The model is capable of predicting gear mesh dynamic load and dynamic transmission error spectra for each gear mesh, dynamic bearing load spectra for each bearing as well as gear body dynamic displacements. Forced vibration response of an example system that consists of three double-planet branches is studied to demonstrate the influence of some of the key design parameters.