Esse trabalho propõe uma implementação numérica para otimização de forma em problemas bi-dimensionais de elasticidade. O objetivo principal é propor uma metodologia eficiente e robusta para solução de problemas de otimização de forma considerando a minimização de concentração de tensões. Na implementação proposta, a análise estrutural é realizada pelo Método dos Elementos de Contorno Simétrico de Galerkin (MECSG), evitando-se assim a dispendiosa etapa de geração da malha. A avaliação das tensões no contorno é obtida por meio de um método preciso, ideal para problemas com concentrações de tensões. Outro aspecto relevante na implementação é a adequada partição das equações do MECSG de forma a reduzir, consideravelmente, o esforço computacional associado à etapa da análise estrutural. O problema de otimização é resolvido utilizando-se um método de otimização moderno, conhecido como Programação Cônica de Segunda Orderm (PCSO). Especificamente, busca-se a reposta do problema de otimização não linear por meio da solução de uma sequência de subproblemas de PCSO.