O objetivo deste trabalho é apresentar uma formulação e uma correspondente implementação computacional para otimização de dimensões de estruturas evitando os problemas de instabilidade apresentados pela formulação convencional. Para atingir este objetivo, a formulação utilizada considera os efeitos da não-linearidade geométrica no comportamento da estrutura e inclui uma restrição sobre a carga de colapso. Elementos finitos reticulados planos e a formulação Lagrangiana Atualizada forma utilizados para análise de estruturas com comportamento geometricamente não- linear. As varáveis de projeto são as alturas das seções transversais dos elementos. O método de Newton-Raphson é utilizado acoplado a diferentes estratégias de incremento de carga e de iteração,tais como as que utilizam a restrição do comprimento de arco e as baseadas no controle dos deslocamentos generalizados, que permitem a ultrapassagem de pontos críticos que possam existir ao longo da trajetória de equilíbrio. Os algoritmos de programação matemática utilizados neste trabalho empregam os gradientes da função objetivo e das restrições, que são calculados com base nos gradientes das respostas da estrutura. Partindo-se das equações gerais de equilíbrio válidas para qualquer elemento, foram desenvolvidas expressões analíticas aproximadas que permitem o cálculo das sensibilidades em relação as variáveis de projeto aproveitando as características da análise.