The goal of the present paper was a com- plete analysis of the dynamic scenario of planetary gears. A lumped mass two-dimensional model is adopted; the model takes into account: time-varying stiffness; nonsmooth nonlinearity due to the backlash, i.e., teeth contact loosing; and bearing compliance. The time-varying meshing stiffness is evaluated by means of a nonlinear finite element model, which allows an accurate evaluation of global and local teeth defor- mation. The dynamic model is validated by compar- isons with the most authoritative literature: linear nat- ural frequencies and nonlinear response. The dynamic scenario is analyzed over a reasonable engineering range in terms of rotation speed and torque. The clas- sical amplitude–frequency diagrams are accompanied by bifurcation diagrams, and for specific regimes, the spectral and topological properties of the response are discussed. Periodic, quasiperiodic and chaotic regimes are found; nonsmooth bifurcations lead period one to period two trajectories. It is found that the bearing com- pliance can influence the natural frequencies combina- tion magnifying the modal interactions due to internal resonances and greatly enlarging the chaotic regions. It is evidenced that the chaotic response indices a sym- metry breaking in the dynamical systems. The physical consequence is that the planetary gearbox under inves- tigation, which is perfectly balanced for each position, can suffer of a big dynamic imbalance when chaotic regimes take place; such imbalance gives rise to alter- nate and unexpected high-level stresses on bearings