Neste trabalho a análise do comportamento de sistemas rotativos do tipo eixo-rotormancal é estendida para incluir os efeitos da presença de mancais hidrodinâmicos na resposta dinâmica. Estes efeitos estão associados à não-linearidade da força de reação exercida pelos suportes sobre o eixo e dependem dos deslocamentos, velocidades transversais e da rotação própia do rotor. A modelagem estrutural do sistema é obtida empregando-se o método dos elementos finitos. O eixo é representado pelo modelo de viga de Timoshenko com dois nós, quatro graus-de-liberdade por nó, e a interpolação do campo de deslocamentos é obtida utilizando-se as funções de Hermite. Os rotores são modelados empregando-se elementos de inércia concentrada associada aos graus-de-liberdade de um ponto nodal do modelo. E, na representação dos mancais hidrodinâmicos utilizou-se a equação de Reynolds, com as hipóteses simplificadoras para mancais curtos, obtendo-se a solução para a distribuição de pressão do filme de óleo em forma fechada. Essa distribuição de pressão permite a obtenção dos coeficientes das matrizes e rigidez e de amortecimento associadas aos graus de liberdade do eixo no ponto nodal de representação do mancal. Para a integração temporal do sistema de equações diferencias utiliza-se o procedimento passo-a-passo, tendo-se implementado os métodos implícitos de Newmark e Wilson – teta, na forma incondicionalmente estável. Devido à não-linearidade das equações obtidas com a presença dos mancais hidrodinâmicos, em cada intervalo de tempo utiliza-se o procedimento de Newton-Raphson modificado para a correção da solução numérica obtida com outros resultados analíticos/numéricos disponíveis na literatura. Também, uma representação numérica para mancais hidrodinâmicos segmentados é apresentada, utilizando-se o desenvolvimento teórico para mancais simples. Neste caso a avaliação do procedimento numérico é fornecida comparando-se a solução numérica com resultados experimentais obtidos dos rotores de usina hidrogenada avaliada pelo CEPEL. Em ambos os procedimentos o rotor idealizado de jeffcott é empregado no estudo de casos. Verifica-se que os principais resultados associados aos efeitos da precessão auto-excitada (oil whirl), de chicoteamento (oil whip), e da estabilização dinâmica do sistema são reproduzidos pelos modelos numéricos utilizados.