Torsional systems with gear pairs such as the gearbox of wind turbines or vehicle driveline systems inherently show impact phenomena due to clearance-type nonlinearities when the system experiences sinusoidal excitation. This research investigates the vibro-impact energy of unloaded gears in geared systems using the harmonic balance method (HBM) in both the frequency and time domains. To achieve accurate simulations, nonlinear models with piecewise and clearance-type nonlinearities and drag torques are defined and implemented in the HBM. Next, the nonlinear frequency responses are examined by focusing on the resonance areas where the impact phenomena occur, along with variations in key parameters such as clutch stiffness, drag torque, and inertias of the flywheel and the unloaded gear. Finally, the effects of the parameters on the vibro-impacts at a specific excitation frequency are explained using bifurcation diagrams. The results are correlated with prior research by defining the gear rattle criteria with key parameters. This article suggests a method to simulate the impact phenomena in torsional systems using the HBM and successfully assesses vibro-impact energy using bifurcation diagrams.