<p>This study introduces a computational framework based on the harmonic balance method for the nonlinear steady-state analysis of rotor systems equipped with hydrodynamic journal and thrust bearings. The journal and thrust bearings are modeled using multi-lobe and tapered-land types, respectively, and their fluid films are discretized using finite elements. The equations of motion for the rotor system are formulated in the frequency domain, incorporating axial-lateral coupling effects induced by the thrust bearing. The proposed framework is validated through a vertical rotor example, where the thrust and journal bearings support static gravitational and dynamic loads, respectively. Numerical results demonstrate that axial-lateral interactions significantly influence the load-carrying capacity and overall system behavior. The harmonic balance method accurately captures the coupled nonlinear steady-state response at a low harmonic order, even in rotor systems that incorporate thrust bearings. Furthermore, the results highlight that thrust and journal bearings should be analyzed simultaneously in a nonlinear manner, as the linearized approach exhibits clear limitations depending on the frequency range and load level.</p>