Neste trabalho é apresentada uma teoria analítica, para o estudo do movimento de satélites artificiais terrestres, considerando perturbações devido ao geopotencial e à força de arrasto atmosférico. O estudo é feito considerando-se o geopotencial até a ordem de J20 e, para a força de arrasto os seguintes efeitos, agindo simultâneamente: a) modelo quadrático para a escala de altura; b) achatamento da atmosfera; c) variação diurna-noturna da densidade atmosférica com o fluxo solar e d) rotação da atmosfera. Variáveis não singulares são introdizidas e, com isto, a teoria pode ser aplicada para satélites com excentricidades e inclinações próximas de zero. A integração do sistema de equações que descreve o problema é feita utilizando-se método para sistemas canônicos generalizados.