Consideramos o modelo de Maier-Saupe com a discretização de Zwanzig para as possíveis orientações das moléculas de cristais líquidos nemáticos. Esse modelo reproduz a transição de fases nemático-isotrópica dos cristais líquidos termotrópicos em acordo com a teoria fenomenológica de Maier e Saupe. Introduzimos variáveis desordenadas de forma a fim de descrever uma eventual estrutura biaxial em uma mistura binária de moléculas uniaxiais oblatas e prolatas. Para uma desordem do tipo quenched, o diagrama de fases do modelo possui uma fase nemática biaxial estável. Para uma desordem do tipo annealed, a estrutura biaxial é termodinamicamente instável. Esses resultados são confirmados realizando um contato com a teoria de Landau-de Gennes. Para ganhar intimidade com os cálculos estatísticos estudamos também um modelo para uma mistura binária de magnetos de Ising na rede. Para ir além dos resultados de campo médio, formulamos o modelo de Maier-Saupe discretizado na rede de Bethe. A análise desse problema é realizada através das relações de recorrência para a função de partição. A transição nemático-isotrópica é localizada através de uma expressão para a energia livre obtida pelo engenhoso método de Gujrati. Considerando o problema das misturas binárias de moléculas nemáticas e utilizando o formalismo adequado à fluidez das partículas annealed, tanto uma análise de estabilidade linear das relações de recorrência quanto uma análise termodinâmica proibem a existência da fase nemática biaxial.