Este trabalho apresenta a aplicação do Método de Otimização Topológica (MOT) considerando restrição de tensão mecânica em dois problemas de Engenharia: o projeto de estruturas mecânicas sujeitas a restrição de tensão e o projeto da distribuição de material em estruturas constituídas por Materiais com Gradação Funcional (MsGF). O MOT é um método numérico capaz de fornecer de forma automática o leiaute básico de uma estrutura mecânica para que esta atenda a um dado requisito de projeto, como o limite sobre a máxima tensão mecânica no componente. Os MsGF são materiais cujas propriedades variam gradualmente com a posição. Este gradiente de propriedades é obtido através da variação contínua da microestrutura formada por dois materiais diferentes. Neste trabalho o MOT foi implementado utilizando o modelo de material Solid Isotropic Microstructure with Penalization (SIMP) e o campo de densidades foi parametrizado utilizando a abordagem Aproximação Contínua da Distribuição de Material (ACDM). O modelo de material e utilizado em conjunto com um localizador de tensões, de modo a representar as tensões nas regiões com densidade intermediária. O projeto de estruturas tradicionais através do MOT possui dois problemas centrais aqui tratados: o fenômeno das topologias singulares, que consiste na incapacidade do algoritmo de otimização de retirar material de certas regiões da estrutura, onde a tensão mecânica supera o limite de tensão quando os valores da densidade tendem a zero, e o problema do grande número de restrições envolvidas, pois que a tensão mecânica é uma grandeza local e deve ser restrita em todos os pontos da estrutura. Para tratar o primeiro problema é utilizado o conceito de relaxação. Para o segundo são utilizadas duas abordagens: uma é a substituição das restrições locais por uma restrição global e a outra é a aplicação do Método do Lagrangeano Aumentado. Ambas foram implementadas e aplicadas para o projeto de estruturas planas e axissimétricas. No projeto da distribuição de material em estruturas constituídas por MsGF é utilizado um modelo de material baseado na interpolação dos limites de Hashin-Shtrikman. A partir deste modelo as tensões em cada fase são obtidas a partir das matrizes localizadoras de tensão. Para tratar o fenômeno das topologias singulares é proposto um índice estimativo de falha, baseado nas tensões de von Mises em cada fase da microestrutura, que evita tal problema. O grande número de restrições é tratado através da restrição global de tensão. Em ambos os problemas as formulações são apresentadas e sua eficiência é discutida através de exemplos numéricos. / This work presents the Topology Optimization Method (TOM) with stress constraint applied to two Engineering problems: the design of mechanical structures subjected to stress constraint and the design of material distribution in structures made of Functionally Graded Materials (FGMs). The TOM is a numerical method capable of synthesizing the basic layout of a mechanical structure accomplishing to a given design requirement, for example the maximum stress in the structure. The FGMs are materials with spatially varying properties, which are obtained through a continuum change of the microstructuremade of two different materials. In this work, the TOM was implemented with Solid Isotropic Microstructure with Penalization (SIMP) material model and the density field was parameterized with the Continuous Approximations of Material Distribution. To obtain the intermediate density stresses, the material model is applied together with a stress localization matrix. The design of mechanical structures through the TOM has two major problems: the singular topology phenomenon, which is characterized by the optimization algorithm impossibility of removing material from certain regions, where the stress overpasses the limiting stress when the density goes to zero, and the large number of constraints, once the stress is a local value that must be constrained everywhere in the structure. To deal with the first problem, it is applied the \"-realaxation concept, and for the second one two approaches are considered: one is to change the local stress constraint into a global stress constraint and the other is to apply the Augmented Lagrangian Method. Both approaches were implemented and applied to the design of plane and axisymmetric structures. In the design of material distribution in structures made of FGMs a material model based on Hashin-Shtrikman bounds is applied. From this model, stresses in each phase are obtained by the stress localization matrix. To deal with the singular topology phenomenon it is proposed a modified von Mises failure criteria index that avoids such problem. A global stress constraint is applied to deal with the large number of constraints. In both problems formulations are presented and their performance are discussed through numerical examples.

Made available in DSpace on 2016-12-12T20:25:11Z (GMT). No. of bitstreams: 1 Luis Reanto Meneghelli.pdf: 5980064 bytes, checksum: 65a0002e42f206e56e3875504a6f0660 (MD5) Previous issue date: 2013-03-07 / Coordenação de Aperfeiçoamento de Pessoal de Nível Superior / Compliant mechanisms are mechanical devices that transform or transfer motion, force or energy through a single part. These mechanisms have important applications in micro electromechanical systems (MEMS) as well as systems that require large accuracy in motion and micro scale. In this work the compliant mechanisms design is performed by means of the Topology Optimization Method, and the optimization problem is formulated in order to maximize the strain energy stored inside the mechanism, eliminating the appearance of hinges. The kinematic behavior of the mechanism is imposed through a set of constraints on displacements of a few degrees of freedom of interest. The elastic behavior is imposed by means of a global stress constraint and some issues associated to the stress parametrization in topology optimization are addressed in the context of mechanisms design. The numerical examples shown that the proposed formulation is able to generate clean topologies of feasible compliant mechanisms. Based on the results, it is clear that the stress constraint has a deep impact on the design of compliant mechanisms, since it can constraint the amount of energy used to enforce the displacement constraints. / Mecanismos flexíveis são dispositivos mecânicos que transformam ou transferem movimento, força ou energia, através de uma única peça. Este tipo de mecanismo encontra aplicações importantes em sistemas micro eletromecânicos (MEMS, micro electromechanical systems) e demais sistemas que exijam grandes precisões nos movimentos e escala microscópica. O projeto de mecanismos flexíveis é realizado através do Método de Otimização Topológica e o problema de otimização será formulado tendo em vista a maximização de energia de deformação elástica armazenada pelo mecanismo, eliminando assim a ocorrência de rótulas (hinges). O comportamento cinemático do mecanismo é imposto através de restrições sobre o campo de deslocamentos em alguns graus de liberdade de interesse. O comportamento elástico dos mecanismos flexíveis é imposto usando um critério global de restrição de tensão e algumas questões importantes associadas a parametrização das tensões são discutidas no contexto de projeto de mecanismos. Os exemplos numéricos mostram que é possível obter topologias bem definidas e que satisfaçam as restrições do projeto. Com base nestes exemplos, verifica-se que a restrição de tensão exerce forte influência no resultado, podendo limitar a quantidade de energia necessária para atender às restrições do mecanismo.

Mecanismos flexíveis, nos quais a deformação elástica é aproveitada na atuação cinemática, têm grande empregabilidade em dispositivos de mecânica de precisão, engenharia biomédica, e mais recentemente em microeletromecanismos (MEMS). Entre as diversas técnicas empregadas para o seu projeto, a otimização topológica tem se mostrado a mais genérica e sistemática. A grande dificuldade destes projetos é conciliar os requisitos cinemáticos com a resistência mecânica da estrutura. Neste trabalho, é implementado um critério de resistência dentro da formulação do problema de otimização, com o intuito de gerar mecanismos que cumpram a tarefa cinemática desejada mas ao mesmo tempo não ultrapassem limites de tensão predeterminados. Esta restrição adicional também visa aliviar o problema bastante conhecido do aparecimento de articulações. Não linearidade geométrica e de material (hiperelasticida de compressível) são implementadas na solução das equações através do método dos elementos finitos para levar em conta os grandes deslocamentos do mecanismo. O método das assíntotas móveis é usado para a atualização das variáveis de projeto. As derivadas do problema de otimização são calculadas analiticamente, pelo método adjunto. Técnicas de projeção são aplicadas para a garantia de topologias livres de instabilidades numéricas comuns em otimização topológica, e projetos otimizados mais próximos de um espaço 0/1 para as densidades físicas. / Compliant me hanisms, in whi h the elasti strain is the basis for kinemati a tua- tion are widely used in pre ision me hani s devi es, biomedi al engineering, and re ently in mi roele trome hani al systems (MEMS). Among several te hniques applied in ompliant me hanisms design, topology optimization has been one of the most general and systemati . The great hallenge in these designs is to ouple both the kinemati s and the me hani al strength riteria requirements. In this work, a strength riteria for the optimization problem is applied, with the aim of generating ompliant me hanisms that ful ll the desired kine- mati tasks while omplying with a stress threshold. The addition of a stress onstraint to the formulation for ompliant me hanisms in topology optimization also aims to allevi- ate the appearan e of hinges in the optimized topology, a well known issue in the design of ompliant me hanisms. Geometri al and material ( ompressible hyperelasti ity) nonlin- earities are applied to the nite element equilibrium equations, to take into a ount large displa ements. The method of moving asymptotes is applied for design variables updating. The derivatives are al ulated analyti ally, by the adjoint method. Proje tion ltering te h- niques are applied, in order to guarantee topologies free of ommon numeri al instabilities in topology optimization, and optimized designs near the 0/1 solution for the physi al densities.

Este trabalho apresenta uma metodologia de otimização topológica visando reduzir o volume de uma estrutura tridimensional sujeita a restrição de tensão. A análise estrutural é feita através do método dos elementos finitos, as tensões são calculadas nos pontos de integração Gaussiana e suavizadas. Para evitar problemas associados a singularidades na tensão é aplicado o método de relaxação de tensão, que penaliza o tensor constitutivo. A norma-p é utilizada para simular a função máximo, que é utilizada como restrição global de tensão. O estimador de erro de Zienkiewicz e Zhu é usado para calcular o erro da tensão, que é considerado durante o cálculo da norma-p, tornando o processo de otimização mais robusto. Para o processo de otimização é utilizada o método de programação linear sequencial, sendo todas as derivadas calculadas analiticamente. É proposto um critério para remoção de elementos de baixa densidade, que se mostrou eficiente contribuindo para gerar estruturas bem definidas e reduzindo significativamente o tempo computacional. O fenômeno de instabilidade de tabuleiro é contornado com o uso de um filtro linear de densidade. Para reduzir o tempo dispendido no cálculo das derivadas e aumentar o desempenho do processo de otimização é proposto um modelo substituto (surrogate model) que é utilizado em iterações internas na programação linear sequencial. O modelo substituto não reduz o tempo de cálculo de cada iteração, entretanto reduziu consideravelmente o número de avaliações da derivada. O algoritmo proposto foi testado otimizando quatro estruturas, e comparado com variações do método e com outros autores quando possível, comprovando a validade da metodologia empregada. / This work presents a methodology for stress-constrained topology optimization, aiming to minimize material volume. Structural analysis is performed by the finite element method, and stress is computed at the elemental Gaussian integration points, and then smoothed over the mesh. In order to avoid the stress singularity phenomenon a constitutive tensor penalization is employed. A normalized version of the p-norm is used as a global stress measure instead of local stress constraint. A finite element error estimator is considered in the stress constraint calculation. In order to solve the optimization process, Sequential Linear Programming is employed, with all derivatives being calculated analiticaly. A criterion is proposed to remove low density elements, contributing for well-defined structures and reducing significantly the computational time. Checkerboard instability is circumvented with a linear density filter. To reduce the computational time and enhance the performance of the code, a surrogate model is used in inner iterations of the Sequential Linear Programming. The present algorithm was evaluated optimizing four structures, and comparing with variations of the methodolgy and results from other authors, when possible, presenting good results and thus verifying the validity of the procedure.

Mecanismos flexíveis, nos quais a deformação elástica é aproveitada na atuação cinemática, têm grande empregabilidade em dispositivos de mecânica de precisão, engenharia biomédica, e mais recentemente em microeletromecanismos (MEMS). Entre as diversas técnicas empregadas para o seu projeto, a otimização topológica tem se mostrado a mais genérica e sistemática. A grande dificuldade destes projetos é conciliar os requisitos cinemáticos com a resistência mecânica da estrutura. Neste trabalho, é implementado um critério de resistência dentro da formulação do problema de otimização, com o intuito de gerar mecanismos que cumpram a tarefa cinemática desejada mas ao mesmo tempo não ultrapassem limites de tensão predeterminados. Esta restrição adicional também visa aliviar o problema bastante conhecido do aparecimento de articulações. Não linearidade geométrica e de material (hiperelasticida de compressível) são implementadas na solução das equações através do método dos elementos finitos para levar em conta os grandes deslocamentos do mecanismo. O método das assíntotas móveis é usado para a atualização das variáveis de projeto. As derivadas do problema de otimização são calculadas analiticamente, pelo método adjunto. Técnicas de projeção são aplicadas para a garantia de topologias livres de instabilidades numéricas comuns em otimização topológica, e projetos otimizados mais próximos de um espaço 0/1 para as densidades físicas. / Compliant me hanisms, in whi h the elasti strain is the basis for kinemati a tua- tion are widely used in pre ision me hani s devi es, biomedi al engineering, and re ently in mi roele trome hani al systems (MEMS). Among several te hniques applied in ompliant me hanisms design, topology optimization has been one of the most general and systemati . The great hallenge in these designs is to ouple both the kinemati s and the me hani al strength riteria requirements. In this work, a strength riteria for the optimization problem is applied, with the aim of generating ompliant me hanisms that ful ll the desired kine- mati tasks while omplying with a stress threshold. The addition of a stress onstraint to the formulation for ompliant me hanisms in topology optimization also aims to allevi- ate the appearan e of hinges in the optimized topology, a well known issue in the design of ompliant me hanisms. Geometri al and material ( ompressible hyperelasti ity) nonlin- earities are applied to the nite element equilibrium equations, to take into a ount large displa ements. The method of moving asymptotes is applied for design variables updating. The derivatives are al ulated analyti ally, by the adjoint method. Proje tion ltering te h- niques are applied, in order to guarantee topologies free of ommon numeri al instabilities in topology optimization, and optimized designs near the 0/1 solution for the physi al densities.

Este trabalho apresenta uma metodologia de otimização topológica visando reduzir o volume de uma estrutura tridimensional sujeita a restrição de tensão. A análise estrutural é feita através do método dos elementos finitos, as tensões são calculadas nos pontos de integração Gaussiana e suavizadas. Para evitar problemas associados a singularidades na tensão é aplicado o método de relaxação de tensão, que penaliza o tensor constitutivo. A norma-p é utilizada para simular a função máximo, que é utilizada como restrição global de tensão. O estimador de erro de Zienkiewicz e Zhu é usado para calcular o erro da tensão, que é considerado durante o cálculo da norma-p, tornando o processo de otimização mais robusto. Para o processo de otimização é utilizada o método de programação linear sequencial, sendo todas as derivadas calculadas analiticamente. É proposto um critério para remoção de elementos de baixa densidade, que se mostrou eficiente contribuindo para gerar estruturas bem definidas e reduzindo significativamente o tempo computacional. O fenômeno de instabilidade de tabuleiro é contornado com o uso de um filtro linear de densidade. Para reduzir o tempo dispendido no cálculo das derivadas e aumentar o desempenho do processo de otimização é proposto um modelo substituto (surrogate model) que é utilizado em iterações internas na programação linear sequencial. O modelo substituto não reduz o tempo de cálculo de cada iteração, entretanto reduziu consideravelmente o número de avaliações da derivada. O algoritmo proposto foi testado otimizando quatro estruturas, e comparado com variações do método e com outros autores quando possível, comprovando a validade da metodologia empregada. / This work presents a methodology for stress-constrained topology optimization, aiming to minimize material volume. Structural analysis is performed by the finite element method, and stress is computed at the elemental Gaussian integration points, and then smoothed over the mesh. In order to avoid the stress singularity phenomenon a constitutive tensor penalization is employed. A normalized version of the p-norm is used as a global stress measure instead of local stress constraint. A finite element error estimator is considered in the stress constraint calculation. In order to solve the optimization process, Sequential Linear Programming is employed, with all derivatives being calculated analiticaly. A criterion is proposed to remove low density elements, contributing for well-defined structures and reducing significantly the computational time. Checkerboard instability is circumvented with a linear density filter. To reduce the computational time and enhance the performance of the code, a surrogate model is used in inner iterations of the Sequential Linear Programming. The present algorithm was evaluated optimizing four structures, and comparing with variations of the methodolgy and results from other authors, when possible, presenting good results and thus verifying the validity of the procedure.

Mecanismos flexíveis, nos quais a deformação elástica é aproveitada na atuação cinemática, têm grande empregabilidade em dispositivos de mecânica de precisão, engenharia biomédica, e mais recentemente em microeletromecanismos (MEMS). Entre as diversas técnicas empregadas para o seu projeto, a otimização topológica tem se mostrado a mais genérica e sistemática. A grande dificuldade destes projetos é conciliar os requisitos cinemáticos com a resistência mecânica da estrutura. Neste trabalho, é implementado um critério de resistência dentro da formulação do problema de otimização, com o intuito de gerar mecanismos que cumpram a tarefa cinemática desejada mas ao mesmo tempo não ultrapassem limites de tensão predeterminados. Esta restrição adicional também visa aliviar o problema bastante conhecido do aparecimento de articulações. Não linearidade geométrica e de material (hiperelasticida de compressível) são implementadas na solução das equações através do método dos elementos finitos para levar em conta os grandes deslocamentos do mecanismo. O método das assíntotas móveis é usado para a atualização das variáveis de projeto. As derivadas do problema de otimização são calculadas analiticamente, pelo método adjunto. Técnicas de projeção são aplicadas para a garantia de topologias livres de instabilidades numéricas comuns em otimização topológica, e projetos otimizados mais próximos de um espaço 0/1 para as densidades físicas. / Compliant me hanisms, in whi h the elasti strain is the basis for kinemati a tua- tion are widely used in pre ision me hani s devi es, biomedi al engineering, and re ently in mi roele trome hani al systems (MEMS). Among several te hniques applied in ompliant me hanisms design, topology optimization has been one of the most general and systemati . The great hallenge in these designs is to ouple both the kinemati s and the me hani al strength riteria requirements. In this work, a strength riteria for the optimization problem is applied, with the aim of generating ompliant me hanisms that ful ll the desired kine- mati tasks while omplying with a stress threshold. The addition of a stress onstraint to the formulation for ompliant me hanisms in topology optimization also aims to allevi- ate the appearan e of hinges in the optimized topology, a well known issue in the design of ompliant me hanisms. Geometri al and material ( ompressible hyperelasti ity) nonlin- earities are applied to the nite element equilibrium equations, to take into a ount large displa ements. The method of moving asymptotes is applied for design variables updating. The derivatives are al ulated analyti ally, by the adjoint method. Proje tion ltering te h- niques are applied, in order to guarantee topologies free of ommon numeri al instabilities in topology optimization, and optimized designs near the 0/1 solution for the physi al densities.

Este trabalho apresenta uma metodologia de otimização topológica visando reduzir o volume de uma estrutura tridimensional sujeita a restrição de tensão. A análise estrutural é feita através do método dos elementos finitos, as tensões são calculadas nos pontos de integração Gaussiana e suavizadas. Para evitar problemas associados a singularidades na tensão é aplicado o método de relaxação de tensão, que penaliza o tensor constitutivo. A norma-p é utilizada para simular a função máximo, que é utilizada como restrição global de tensão. O estimador de erro de Zienkiewicz e Zhu é usado para calcular o erro da tensão, que é considerado durante o cálculo da norma-p, tornando o processo de otimização mais robusto. Para o processo de otimização é utilizada o método de programação linear sequencial, sendo todas as derivadas calculadas analiticamente. É proposto um critério para remoção de elementos de baixa densidade, que se mostrou eficiente contribuindo para gerar estruturas bem definidas e reduzindo significativamente o tempo computacional. O fenômeno de instabilidade de tabuleiro é contornado com o uso de um filtro linear de densidade. Para reduzir o tempo dispendido no cálculo das derivadas e aumentar o desempenho do processo de otimização é proposto um modelo substituto (surrogate model) que é utilizado em iterações internas na programação linear sequencial. O modelo substituto não reduz o tempo de cálculo de cada iteração, entretanto reduziu consideravelmente o número de avaliações da derivada. O algoritmo proposto foi testado otimizando quatro estruturas, e comparado com variações do método e com outros autores quando possível, comprovando a validade da metodologia empregada. / This work presents a methodology for stress-constrained topology optimization, aiming to minimize material volume. Structural analysis is performed by the finite element method, and stress is computed at the elemental Gaussian integration points, and then smoothed over the mesh. In order to avoid the stress singularity phenomenon a constitutive tensor penalization is employed. A normalized version of the p-norm is used as a global stress measure instead of local stress constraint. A finite element error estimator is considered in the stress constraint calculation. In order to solve the optimization process, Sequential Linear Programming is employed, with all derivatives being calculated analiticaly. A criterion is proposed to remove low density elements, contributing for well-defined structures and reducing significantly the computational time. Checkerboard instability is circumvented with a linear density filter. To reduce the computational time and enhance the performance of the code, a surrogate model is used in inner iterations of the Sequential Linear Programming. The present algorithm was evaluated optimizing four structures, and comparing with variations of the methodolgy and results from other authors, when possible, presenting good results and thus verifying the validity of the procedure.

Este trabalho apresenta a aplicação do Método de Otimização Topológica (MOT) considerando restrição de tensão mecânica em dois problemas de Engenharia: o projeto de estruturas mecânicas sujeitas a restrição de tensão e o projeto da distribuição de material em estruturas constituídas por Materiais com Gradação Funcional (MsGF). O MOT é um método numérico capaz de fornecer de forma automática o leiaute básico de uma estrutura mecânica para que esta atenda a um dado requisito de projeto, como o limite sobre a máxima tensão mecânica no componente. Os MsGF são materiais cujas propriedades variam gradualmente com a posição. Este gradiente de propriedades é obtido através da variação contínua da microestrutura formada por dois materiais diferentes. Neste trabalho o MOT foi implementado utilizando o modelo de material Solid Isotropic Microstructure with Penalization (SIMP) e o campo de densidades foi parametrizado utilizando a abordagem Aproximação Contínua da Distribuição de Material (ACDM). O modelo de material e utilizado em conjunto com um localizador de tensões, de modo a representar as tensões nas regiões com densidade intermediária. O projeto de estruturas tradicionais através do MOT possui dois problemas centrais aqui tratados: o fenômeno das topologias singulares, que consiste na incapacidade do algoritmo de otimização de retirar material de certas regiões da estrutura, onde a tensão mecânica supera o limite de tensão quando os valores da densidade tendem a zero, e o problema do grande número de restrições envolvidas, pois que a tensão mecânica é uma grandeza local e deve ser restrita em todos os pontos da estrutura. Para tratar o primeiro problema é utilizado o conceito de relaxação. Para o segundo são utilizadas duas abordagens: uma é a substituição das restrições locais por uma restrição global e a outra é a aplicação do Método do Lagrangeano Aumentado. Ambas foram implementadas e aplicadas para o projeto de estruturas planas e axissimétricas. No projeto da distribuição de material em estruturas constituídas por MsGF é utilizado um modelo de material baseado na interpolação dos limites de Hashin-Shtrikman. A partir deste modelo as tensões em cada fase são obtidas a partir das matrizes localizadoras de tensão. Para tratar o fenômeno das topologias singulares é proposto um índice estimativo de falha, baseado nas tensões de von Mises em cada fase da microestrutura, que evita tal problema. O grande número de restrições é tratado através da restrição global de tensão. Em ambos os problemas as formulações são apresentadas e sua eficiência é discutida através de exemplos numéricos. / This work presents the Topology Optimization Method (TOM) with stress constraint applied to two Engineering problems: the design of mechanical structures subjected to stress constraint and the design of material distribution in structures made of Functionally Graded Materials (FGMs). The TOM is a numerical method capable of synthesizing the basic layout of a mechanical structure accomplishing to a given design requirement, for example the maximum stress in the structure. The FGMs are materials with spatially varying properties, which are obtained through a continuum change of the microstructuremade of two different materials. In this work, the TOM was implemented with Solid Isotropic Microstructure with Penalization (SIMP) material model and the density field was parameterized with the Continuous Approximations of Material Distribution. To obtain the intermediate density stresses, the material model is applied together with a stress localization matrix. The design of mechanical structures through the TOM has two major problems: the singular topology phenomenon, which is characterized by the optimization algorithm impossibility of removing material from certain regions, where the stress overpasses the limiting stress when the density goes to zero, and the large number of constraints, once the stress is a local value that must be constrained everywhere in the structure. To deal with the first problem, it is applied the \"-realaxation concept, and for the second one two approaches are considered: one is to change the local stress constraint into a global stress constraint and the other is to apply the Augmented Lagrangian Method. Both approaches were implemented and applied to the design of plane and axisymmetric structures. In the design of material distribution in structures made of FGMs a material model based on Hashin-Shtrikman bounds is applied. From this model, stresses in each phase are obtained by the stress localization matrix. To deal with the singular topology phenomenon it is proposed a modified von Mises failure criteria index that avoids such problem. A global stress constraint is applied to deal with the large number of constraints. In both problems formulations are presented and their performance are discussed through numerical examples.

[pt] As metodologias de projeto estrutural foram fortemente influenciadas pelo advento da computação. Os avanços nas áreas de análise numérica, como o método dos elementos finitos, e os softwares de Desenho Assistido por Computador, literalmente ajudaram a moldar o mundo como ele é hoje. Implementações computacionais das técnicas de otimização estrutural, como a otimização topológica, permitem a determinação das estruturas base, gerando uma grande quantidade de projetos novos, mais eficientes, com o potencial de mudar drasticamente o futuro das aeronaves, automóveis, edifícios, etc. Introduzir restrições de tensão na otimização topológica tradicional permite a obtenção de soluções mais seguras e confiáveis que se assemelhem mais à estrutura final. Contudo, isto não é uma tarefa trivial, apresentando várias dificuldades conceituais e numéricas. Nesta dissertação, as principais questões deste problema são discutidas e as técnicas presentes hoje na literatura são revisadas e criticadas quanto aos seus desempenhos. A principal contribuição deste trabalho é uma nova técnica baseada no Método do Lagrangiano Aumentado que lida eficientemente com um grande número de restrições. Em contraste com os métodos existentes, que são dependentes do problema e da malha, a abordagem proposta apresenta poucos parâmetros que precisam ser ajustados a cada novo caso. Para avaliar suas potencialidades, desenvolveu-se um código em MATLAB, eficaz e robusto. Diversos exemplos representativos, incluindo problemas de larga escala, são apresentados. Finalmente, as soluções obtidas, incluindo algumas complicações inesperadas, são discutidas detalhadamente e sugestões para trabalhos futuros são propostas. / [en] Structural design methodologies were strongly influenced by the advent of computing. The advances in numerical analyses, such as the finite element method, and Computer Aided Design software have literally helped shape the engineering world as it is today. Structural optimization methods such as topology optimization aim to take the next step by letting the computer guide the design, in order to achieve new and more eficient designs. This approach has the potential to change the future of various industries, including aircraft, automobile, construction, etc. The introduction of stress constraints on traditional topology optimization allows for safer and more reliable solutions that will more closely resemble the final structure. The successful solution of this problem poses several conceptual and numerical dificulties. Thus this dissertation details the main issues of this problem and reviews the current techniques discussed in the literature including some critiques of their performance. The main contribution of this work is a novel technique based on the Augmented Lagrangian method that can eficiently handle a large number of constraints. In contrast to existing methods which are both problem- and mesh-dependent, the presented approach contains only a few parameters which need to be adjusted for each new case. In order to verify the technique s capabilities, a user friendly MATLAB code was developed that is both effective and robust. Several representative examples, including large-scale problems, are presented. Finally, the solutions obtained here, including some unexpected complications, are thoroughly discussed and suggestions for future work are also addressed.