Este trabalho possui como objetivo avaliar a adequação dos modelos de fadiga multiaxial no que diz respeito à sua aplicação para eixos virabrequins, bem como avaliar a possibilidade de falha dos eixos virabrequins quando submetidos a condições de carregamentos calculadas para a operação. A fadiga é uma falha mecânica causada primariamente pela aplicação de carregamentos variáveis no tempo, tendo como principal característica a nucleação e/ou propagação de trincas. É de extrema importância que o estudo e a modelagem das falhas por fadiga seja precisa e confiável, uma vez que a nucleação e a propagação de trincas, em geral, não provocam mudança evidente no comportamento global da estrutura. Isso se deve ao fato de que o dano gerado por um trincamento fica, quase sempre, restrito à região crítica da peça, tendendo a não gerar avisos prévios de falha iminente. Tendo isso em vista, as previsões de falhas por fadiga tornam-se críticas, tanto para carregamentos uniaxiais, quanto para carregamentos multiaxiais. Entretanto, para condições de carregamento multiaxiais, é requerido análises baseadas em modelos mais avançados do que os adotados para condições uniaxiais, uma vez que seu estado de tensões é extremamente complexo. Nesse contexto, critérios foram desenvolvidos para uma tentativa de prever falhas de componentes submetidos à carregamentos multiaxiais cíclicos. Um exemplo de componente mecânico sujeito a fadiga multiaxial, quando em serviço, são os eixos virabrequins de unidades geradoras de centrais termoelétricas. Desta forma, o presente trabalho analisou a aplicabilidade de sete modelos (Findley, Matake, McDiarmid, Susmel e Lazzarin, Papadopoulos, Carpinteri e Spagnoli e Liu e Mahadevan) de fadiga multiaxial de alto ciclo, baseados no plano crítico, para a previsão de falhas em eixos virabrequins de unidades geradoras de centrais termoelétricas. Para a aplicação dos modelos, condições de carregamento reportadas, na literatura, como críticas para a falha e propriedades do material (limites de resistência à fadiga) previamente obtidas experimentalmente serão utilizadas como dados de entrada e, a partir disso, Índices de Erro (IE) serão gerados para cada um dos sete modelos avaliados. Observou-se divergências entre os resultados para cada modelo, sendo alguns mais conservadores do que outros, justamente devido às condições de carregamento e às propriedades do material utilizadas. Apesar disso, o modelo de Papadopoulos se mostrou como o mais adequado entre os sete modelos, devido a sua facilidade na aplicação de fórmula e precisão na previsão da falha.