As indústrias aeronáuticas cada vez mais utilizam materiais compósitos na estrutura primária das aeronaves, visando a redução do peso e redução no tempo de produção da mesma. A sociedade, representada pelas agências de homologação aeronáutica, exige níveis de segurança crescentes na atividade de transporte aéreo. A complexidade no cálculo das cargas dinâmicas para as diversas condições requeridas pelos regulamentos, devido a utilização de programas de computador capazes de considerar todas as variáveis envolvidas. Este trabalho apresenta uma metodologia para o cálculo da velocidade de instabilidade aeroelástica que utiliza uma formulação de elementos finitos, adequada aos materiais compósitos constituídos de múltiplas camadas ortotrópicas de dupla curvatura. As freqüências de vibração livre e os respectivos modos são determinados a partir da rigidez e da massa dos elementos. As forças aerodinâmicas são determinadas, no regime subsônico, utilizando a teoria do escoamento potencial não estacionário linear, artavés de algoritmos de alta precisão com as mais recentes formulações desenvolvidas. Os resultados numéricos, referentes a asas retangulares planas, São comparados com resultados experimentais. Adicionalmente são apresentados resultados numéricos para as asas retangulares constituída de laminado com orientação arbitrária.