Este trabalho é uma contribuição para o conhecimento de metodologias de projeto de estruturas de material composto, aplicando métodos de otimização estrutural a cascas laminadas e apresentando uma estratégia em dois níveis. No primeiro nível é realizada a minimização da flexibilidade da estrutura, tendo como variável de projeto a orientação de cada lâmina da estrutura. Utiliza-se Programação Linear Seqüencial (SLP) e direção de tensão principal para otimização da orientação. No segundo nível minimiza-se o volume de cada lâmina, usando a flexibilidade total da estrutura como restrição e a densidade relativa como variável de projeto, também através de SLP. Para evitar aparecimento de áreas com densidades intermediárias, utiliza-se um Método de Continuação, dividindo o nível de otimização topológica em duas ou mais etapas. As formulações desenvolvidas permitem a solução de problemas com múltiplos casos de carregamento. Para a solução da equação de equilíbrio de casca laminada, utiliza-se um elemento finito de casca degenerado de oito nós com integração explícita na direção da espessura. A implementação desse elemento é feita de modo a facilitar a obtenção das derivadas da matriz de rigidez, necessárias na linearização das funções objetivo e restrições. Evita-se assim o uso de derivadas numéricas. Resultados para vários tipos de estrutura são apresentados, incluindo comparações entre diferentes carregamentos, condições de contorno, número de lâminas, espessuras, etc. As soluções obtidas, formas de análise e possíveis aplicações são discutidas.

Orientador: Francisco de Assis Magalhães Gomes Neto / Dissertação (mestrado) - Universidade Estadual de Campinas, Instituto de Matematica, Estatistica e Computação Cientifica / Made available in DSpace on 2018-08-13T12:28:53Z (GMT). No. of bitstreams: 1 Senne_ThadeuAlves_M.pdf: 1074089 bytes, checksum: 76e0320b61aeb0cef5475189e51d218e (MD5) Previous issue date: 2009 / Resumo: Neste trabalho, estudamos algumas formulações possíveis para o problema de otimização topológica de um mecanismo flexível, propostas por Nishiwaki et al. [33], Lima [26] e Sigmund [37]. Para resolver os problemas de programação não linear associados a cada uma das formulações estudadas, usamos uma versão globalmente convergente da Programação Linear Seqüencial, inspirada no trabalho de Gomes et al. [18], e uma versão globalmente convergente do Método das Assíntotas Móveis, desenvolvida por Svanberg [46]. Fazemos uma análise comparativa do desempenho desses dois métodos de otimização, no que diz respeito às topologias ótimas obtidas para as estruturas e ao esforço computacional para a resolução dos problemas de otimização topológica. Comparamos também a eficácia de alguns filtros espaciais propostos na literatura, que têm o papel de evitar o aparecimento de regiões semelhantes a um tabuleiro de xadrez nas topologias ótimas das estruturas / Abstract: In this work, we study some possible formulations for the topology optimization problem of a compliant mechanism, proposed by Nishiwaki et al. [33], Lima [26] and Sigmund [37]. To solve the nonlinear programming problem associated to each formulation, we use a globally convergent version of the Sequential Linear Programming, inspired in the Gomes' et al. [18] work, and a globally convergent version of the Method of Moving Asymptotes, developed by Svanberg [46]. We make a comparative analysis of the performance of these two optimization methods, with respect to the optimum topologies obtained for the structures and to the computational e ort for the resolution of the topology optimization problems. Also, we compare the e ciency of some spatial lters already proposed in the literature, used to avoid the occurrency of regions similar to a checkerboard in the optimum topology of the structures / Mestrado / Matematica Aplicada - Otimização / Mestre em Matemática Aplicada

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.

Em sentido geral, a quimiotaxia é o movimento dirigido que desenvolvem alguns seres vivos em resposta aos gradientes químicos presentes no seu ambiente. Uma bactéria é um organismo unicelular que usa a quimiotaxia como mecanismo de mobilização para encontrar os nutrientes de que precisa para sobreviver e para escapar de ambientes nocivos. Evoluída durante milhões de anos pela natureza, a quimiotaxia de bactérias é um processo altamente otimizado de busca e exploração em espaços desconhecidos. Graças aos avanços no campo da computação, as estratégias quimiotácticas das bactérias e sua excelente capacidade de busca podem ser modeladas, simuladas e emuladas para desenvolver métodos de otimização inspirados na natureza que sejam uma alternativa aos métodos já existentes. Neste trabalho, desenvolvem-se dois algoritmos baseados em estratégias quimiotácticas de bactérias: o BCBTOA (Bacterial Chemotaxis Based Topology Optimization Algorithm) e o BCMOA (Bacterial Chemotaxis Multiobjective Optimization Algorithm) os quais são um algoritmo de otimização topológica e um algoritmo de otimização multi-objetivo, respectivamente. O desempenho dos algoritmos é avaliado mediante a sua aplicação à solução de diversos problemas de prova e os resultados são comparados com os de outros algoritmos atualmente relevantes. O algoritmo de otimização multi-objetivo desenvolvido, também foi aplicado na solução de três problemas de otimização de projeto mecânico de eixos. Os resultados obtidos e os analise comparativos feitos, permitem concluir que os algoritmos desenvolvidos são altamente competitivos e demonstram o potencial do processo de quimiotaxia de bactérias como fonte de inspiração de algoritmos de otimização distribuída, contribuindo assim, a dar resposta à constante demanda por técnicas de otimização mais eficazes e robustas. / In general, chemotaxis is the biased movement developed by certain living organisms as a response to chemical gradients present in their environment. A bacterium is a unicellular organism that uses chemotaxis as a mechanism for mobilization that allows it to find nutrients needed to survive and to escape from harmful environments. Millions of years of natural evolution became bacterial chemotaxis a highly optimized process in searching and exploration of unknown spaces. Thanks to advances in the computing field, bacterial chemotactical strategies and its excellent ability in searching can be modeled, simulated and emulated developing bio-inspired optimization methods as alternatives to classical methods. Two algorithms based on bacterial chemotactical strategies were designed, developed and implemented in this work: i) the topology optimization algorithm, BCBTOA (Bacterial Chemotaxis Based Topology Optimization Algorithm) and ii) the multi-objective optimization algorithm, BCMOA (Bacterial Chemotaxis Multiobjective Optimization Algorithm). Algorithms performances were evaluated by their applications in the solution of benchmark problems and the results obtained were compared with other algorithms also relevant today. The BCMOA developed here was also applied in the solution of three mechanical design problems. The results obtained as well as the comparative analysis conducted lead to conclude that the algorithms developed were competitive. This also demonstrates the potential of bacterial chemotaxis as a process in which distributed optimization techniques can be inspired.

Orientador: Francisco de Assis Magalhães Gomes Neto / Tese (doutorado )- Universidade Estadual de Campinas, Instituto de Matematica, Estatistica e Computação Cientifica / Made available in DSpace on 2018-08-22T00:37:54Z (GMT). No. of bitstreams: 1 Senne_ThadeuAlves_D.pdf: 4979548 bytes, checksum: 0317347b212e7649dd09e5d5388e618e (MD5) Previous issue date: 2013 / Resumo: Nos últimos anos, a otimização topológica vem sendo amplamente adotada nas indústrias automotiva e aeroespacial, e no projeto de um tipo especial de estruturas, denominado mecanismo flexível. Grande parte dos trabalhos na área de otimização topológica considera que a estrutura possui uma relação linear entre deformações e deslocamentos, ou seja, supõe-se que os deslocamentos sofridos pela estrutura sejam pequenos. Todavia, para algumas estruturas, essa hipótese não é válida, sendo necessário supor que os deslocamentos são grandes, o que implica numa relação não linear entre deformações e deslocamentos. Nesse caso, dizemos que a estrutura está sob não linearidade geométrica. O objetivo desta tese de doutorado é a obtenção da topologia ótima de estruturas e de mecanismos flexíveis sob não linearidade geométrica através um novo algoritmo de otimização, denominado Programação Linear por Partes Sequencial (PLPS). Este método consiste na resolução de subproblemas de programação linear por partes convexas, onde são introduzidas informações sobre a diagonal da matriz Hessiana da função objetivo. Para acelerar o algoritmo, tais subproblemas são convertidos em problemas de programação linear. Provamos que a PLPS é globalmente convergente a pontos estacionários. Além disso, nossos experimentos numéricos realizados com estruturas e mecanismos flexíveis sujeitos a grandes deslocamentos mostram que a PLPS é eficiente e robusta / Abstract: In the last years, topology optimization has been broadly applied in the automotive and aerospatial industries, and to a special kind of structure, named compliant mechanism. Most papers on topology optimization consider that the structure has a linear relation between strains and displacements, meaning that the displacements of the structure are small. However, for some structures this assumption is not valid, and it is necessary to suppose that the displacements are large, implying in a nonlinear relation between strains and displacements. In this case, we say that the structure is under geometrical nonlinearity. The objective of this doctoral thesis is to obtain the optimum topology of structures and compliant mechanisms under geometrical nonlinearity through a new optimization algorithm, named Sequential Piecewise Linear Programming (SPLP). This method consists in the solution of convex piecewise linear programming subproblems that contain information about the diagonal of the Hessian matrix of the objective function. To speed up the algorithm, these subproblems are converted into linear programming ones. We prove that the SPLP is globally convergent to stationary points. Besides, our numerical experiments with structures and compliant mechanisms under large displacements also show that the SPLP is efficient and robust / Doutorado / Matematica Aplicada - Otimização / Doutor em Matemática Aplicada

O conteúdo deste trabalho trata da formulação para solução do problema de otimização estrutural com minimização de massa em estruturas unifilares, sujeitas a restrição de tensão, flambagem das barras isoladas e fadiga. São considerados três casos de otimização: paramétrica, de forma e dimensional. Os problemas de singularidades nas restrições de tensão e flambagem são evitados através de uma formulação que faz uso de programação inteira para solução do problema. Outra singularidade encontrada na otimização topológica é a singularidade na matriz de rigidez da estrutura. Este problema foi evitado através de uma formulação que considera a existência de matriz de rigidez regular como restrição do problema. O método de solução utilizado para resolver problema de otimização é o método dos algoritmos genéticos. As restrições do problema são impostas através da penalização da função objetivo. O método de solução mostrou-se adequado para solução dos problemas estudados. A formulação implementada é validada através da solução de problemas clássicos de otimização estrutural. Os resultados obtidos são comparados com a literatura onde verificou-se a coerência dos mesmos. Após realizar a validação, a formulação é utilizada em um estudo que tem como base uma estrutura real: uma torre de queima de gases (flare) oriundos do processo de extração e armazenagem de petróleo em uma unidade flutuante. Para o problema da torre as restrições foram determinadas com base em critérios de falha estabelecido na norma DNV. A otimização do flare permitiu minimizar a massa da estrutura sem que os critérios de falha fossem violados. Verificou-se que a metodologia proposta é adequada para solução com grande número de restrições e com diversos casos de carregamento. / The purpose of this work is the development of a methodology to solve the structural optimization problem of frame structures subject to stress, buckling of isolated members, and fatigue constraints. Three types of structural optimization problems are considered: sizing, shape and topological. The stress and buckling singularity problems are avoided by an integer design variable formulation, using integer programing to obtain the optimization problem solution. Another issue found in optimization problems is the stiffness matrix singularity. The proposed formulations include the linear system stability as a constraint in the optimization problem. A genetic algorithm is used to solve the general optimization problem. All constraints of the problem are included with a penalization equation. The results show that genetic algorithm is a good approach to solve the proposed formulation. The proposed formulation is tested for solving classical optimization problems. The obtained results are consistent with the literature. A real engineering problem is solved with proposed methodology: a gas burning tower (flare). In this problem, all constraints are based on failure criteria recommended by DNV standards. The structural optimization of this problem shows that structural mass minimization is possible without violating the failure criteria. It is observed that solution methodology deals successfully with problems with multiple constraints and load cases

A Tomografia de Impedância Elétrica é uma técnica de obtenção de imagens do interior de um corpo, mediante grandezas elétricas medidas em sua superfície. Matematicamente, a técnica determina as distribuições de condutividades e permissividades elétricas num dado modelo do corpo, as quais reproduzem as medidas de correntes e potenciais elétricos em eletrodos fixados ao corpo. Nesse caso, as distribuições de condutividades e permissividades representam a solução de um problema não-linear e mal-posto, o qual é instável e apresenta mínimos locais, requerendo técnicas de inversão específicas. Um sistema de Tomografia de Impedância Elétrica aplicado à obtenção de imagens de valores absolutos possui, atualmente, limitações. São algumas delas a obtenção de distribuições de propriedades suaves e de valores geralmente subestimados, a sensibilidade elevada ao erro de posicionamento dos eletrodos (com relação ao modelo) e ao erro nos valores de parâmetros de contato, a sensibilidade elevada aos ruídos de medição, os tempos elevados de processamento, dentre outros. Com o intuito de abordar as limitações, melhorando o desempenho do sistema de Tomografia de Impedância Elétrica de imagens absolutas, são então propostas e avaliadas ferramentas baseadas no Método de Otimização Topológica no atual trabalho. Mais especificamente, avaliam-se: 1) um método para obtenção de parâmetros de contato em conjunto com uma imagem e um método de regularização baseado no controle explícito da variação espacial da imagem, 2) uma formulação para acomodação de incertezas, 3) uma formulação para correção do posicionamento de eletrodos, 4) uma formulação para projeto de eletrodos e 5) um novo solucionador de sistemas lineares de larga escala. Os resultados mostram a efetividade da maioria das técnicas propostas, e sugerem os novos tópicos de pesquisa em Tomografia de Impedância Elétrica. / Electrical Impedance Tomography images the interior of a body based on electrical quantities measured on the surface of it. Mathematically, the technique finds the electric admittivity distribution in a given body model which reproduces the boundary measurements of electric currents and potentials on electrodes attached to that body. Therefore, the admittivity distribution is the solution of a non-linear and ill-posed problem, which is unstable and have local minima, requiring specific inversion techniques. Electrical Impedance Tomography systems which obtain images corresponding to absolute values present limitations. For instance, the results are usually smooth and underestimated, the sensitivity to errors in the positioning of electrodes and wrong values of contact parameters and the sensitivity to measurement noise are high, the data processing time is high, etc. In this work, techniques based on the Topology Optimization Method intended for improving the performance of the particular Electrical Impedance Tomography system applied to absolute images are proposed and evaluated. More specifically, the following techniques are evaluated: 1) a method intended to obtain contact parameters together with images, and a regularization method based on the explicit control of the spatial variation regarding the image, 2) a formulation applied to handle uncertainties, 3) a formulation applied to correct the position of electrodes, 4) a formulation applied to design electrodes, 5) and a new solver for large-scale linear systems. Results show the effectiveness of most of the proposed techniques, and suggest new research topics in Electrical Impedance Tomography.

A Tomografia de Impedância Elétrica é uma técnica de obtenção de imagens do interior de um corpo, mediante grandezas elétricas medidas em sua superfície. Matematicamente, a técnica determina as distribuições de condutividades e permissividades elétricas num dado modelo do corpo, as quais reproduzem as medidas de correntes e potenciais elétricos em eletrodos fixados ao corpo. Nesse caso, as distribuições de condutividades e permissividades representam a solução de um problema não-linear e mal-posto, o qual é instável e apresenta mínimos locais, requerendo técnicas de inversão específicas. Um sistema de Tomografia de Impedância Elétrica aplicado à obtenção de imagens de valores absolutos possui, atualmente, limitações. São algumas delas a obtenção de distribuições de propriedades suaves e de valores geralmente subestimados, a sensibilidade elevada ao erro de posicionamento dos eletrodos (com relação ao modelo) e ao erro nos valores de parâmetros de contato, a sensibilidade elevada aos ruídos de medição, os tempos elevados de processamento, dentre outros. Com o intuito de abordar as limitações, melhorando o desempenho do sistema de Tomografia de Impedância Elétrica de imagens absolutas, são então propostas e avaliadas ferramentas baseadas no Método de Otimização Topológica no atual trabalho. Mais especificamente, avaliam-se: 1) um método para obtenção de parâmetros de contato em conjunto com uma imagem e um método de regularização baseado no controle explícito da variação espacial da imagem, 2) uma formulação para acomodação de incertezas, 3) uma formulação para correção do posicionamento de eletrodos, 4) uma formulação para projeto de eletrodos e 5) um novo solucionador de sistemas lineares de larga escala. Os resultados mostram a efetividade da maioria das técnicas propostas, e sugerem os novos tópicos de pesquisa em Tomografia de Impedância Elétrica. / Electrical Impedance Tomography images the interior of a body based on electrical quantities measured on the surface of it. Mathematically, the technique finds the electric admittivity distribution in a given body model which reproduces the boundary measurements of electric currents and potentials on electrodes attached to that body. Therefore, the admittivity distribution is the solution of a non-linear and ill-posed problem, which is unstable and have local minima, requiring specific inversion techniques. Electrical Impedance Tomography systems which obtain images corresponding to absolute values present limitations. For instance, the results are usually smooth and underestimated, the sensitivity to errors in the positioning of electrodes and wrong values of contact parameters and the sensitivity to measurement noise are high, the data processing time is high, etc. In this work, techniques based on the Topology Optimization Method intended for improving the performance of the particular Electrical Impedance Tomography system applied to absolute images are proposed and evaluated. More specifically, the following techniques are evaluated: 1) a method intended to obtain contact parameters together with images, and a regularization method based on the explicit control of the spatial variation regarding the image, 2) a formulation applied to handle uncertainties, 3) a formulation applied to correct the position of electrodes, 4) a formulation applied to design electrodes, 5) and a new solver for large-scale linear systems. Results show the effectiveness of most of the proposed techniques, and suggest new research topics in Electrical Impedance Tomography.

A otimização topológica de estruturas está relacionada à concepção de projetos que executem suas funções com nível de segurança adequado empregando a quantidade mínima de material. Neste trabalho, determina-se a geometria ótima de estruturas planas por meio do acoplamento do Método dos Elementos de Contorno (MEC) ao Método Level Set (MLS). O algoritmo é composto por 3 etapas: problema mecânico, otimização topológica e reconstrução da estrutura. O problema mecânico é resolvido pelas equações algébricas do MEC. A otimização topológica é determinada pelo MLS, este representa a geometria do corpo e suas evoluções por meio da função Level Set (LS) avaliada em seu nível zero. Na reconstrução realiza-se o remalhamento, pois a cada iteração a estrutura é modificada. O acoplamento proposto resulta na geometria ótima da estrutura sem a necessidade da aplicação de filtros. Os exemplos analisados mostram que algoritmo desenvolvido capta adequadamente a geometria ótima das estruturas. Com esse trabalho, avança-se no campo das aplicações do acoplamento MEC-MLS e no desenvolvimento de soluções inovadoras para problemas complexos de engenharia. / In general, the topological optimization of structures is related to design projects that perform their functions with appropriate security levels using the minimum amount of material. This research determines the optimal geometry of 2D structures by coupling the Boundary Blement Method (BEM) to Level Set Method (LSM). The algorithm consists of 3 steps: mechanical model, topology optimization and structure reconstruction. The mechanical model is solved by BEM algebraic equations. The topology optimization is determined using the MLS, the geometry of the body is determined by the Level Set (LS) function evaluated at the zero level. The reconstruction achieves the remeshing, because for each iteration of the structure is modified. The proposed coupling results in the optimal geometry of the structure without the filters application. The examples show that the algorithm developed captures adequately the optimal geometry of the structures. With this dissertation, it is possible advance in the field of applications of the BEM - LSM and develop innovative solutions to complex engineering problems.

Um micromecanismo é, essencialmente, um dispositivo de dimensões milimétricas ou até micrométricas que executa uma tarefa específica como atuar como garra, pinça, grampo, etc. Quando acoplados a um sistema eletrônico, são chamados de sistemas microeletromecânicos ou \"Micro-Electro-Mechanical Systems\" (MEMS). Esses dispositivos são quase todos constituídos por mecanismos flexíveis, onde o movimento é dado pela flexibilidade de sua estrutura, sem juntas e pinos. Uma das formas de atuação de micromecanismos é a eletrotermomecânica, onde uma atuação elétrica sobre o próprio mecanismo é convertida em calor, por efeito Joule, que gera tensões térmicas responsáveis pela deformação estrutural desejada. Recentemente, vários grupos de pesquisa no mundo estão desenvolvendo micromecanismos fabricados com dois (ou até mais) materiais, o que permite obter maiores deformações sem que seja excedido o limite de resistência do material e mais flexibilidade no projeto de micromecanismos que realizem diferentes tarefas quando sujeito a diversas atuações (multiflexíveis). As técnicas de processo de fabricação de micromecanismos atingiram um alto nível de maturidade. No entanto, a modelagem e, em particular, o desenvolvimento de métodos computacionais sistemáticos para o projeto estão ainda no seu estágio inicial. Atualmente, o projeto de micromecanismos com vários materiais vem sendo realizado por métodos de tentativa e erro, dependendo da intuição e experiência do projetista. Além disso, o projeto genérico de um MEMS eletrotermomecânico é uma tarefa complexa, que leva em conta conhecimentos multidisciplinares. Dessa forma, o objetivo desse trabalho de mestrado foi desenvolver um software para o projeto de MEMS multifásicos, atuados eletrotermicamente, usando um método de projeto genérico e sistemático, como o Método de Otimização Topológica (MOT). Utilizando um modelo de interpolação de material de função de pico, qualquer número de materiais pode ser considerado sem que haja aumento na quantidade de variáveis de projeto se comparado à otimização com apenas um material e vazio. Visando maximizar o deslocamento de saída contra uma peça de rigidez conhecida, foram projetados mecanismos atuados por tensão elétrica, alguns considerando multiflexibilidade. Um estudo da influência dos parâmetros da otimização foi realizado. Como uma alternativa à atuação eletrotermomecânica, foram projetados mecanismos atuados por fluxo de calor. / A micromechanism is essentially a device of milimetric, or even micrometric, dimensions that can actuate as a gripper, tweezers, clamp, etc. When coupled to an electronic system, they are called \"Micro-Electro-Mechanical Systems\" (MEMS). Almost all of these devices are constituted by compliant mechanisms, where the motion is allowed by the compliance of its own structure, rather than the presence of joint and pins. One of the forms of micromechanisms actuation is the electrothermomechanical, where an electric actuation applied to the mechanism is converted in heat, by Joule effect, that generates the thermal stress responsible for the desired structural deformation. Recently, many research groups around the world are developing micromechanisms manufactured with two (or even more) materials, what allows larger displacements without exceeding the materials ultimate tensile strength, and gives more flexibility in the design of micromechanisms that accomplish different tasks when under different actuations (multiflexible mechanisms). The manufacturing process techniques of micromechanisms reached a high level of maturity, however, the modelling and, particularly, the development of systematic computational methods for design are still in early stages. Nowadays, micromechanism design with many materials is being carried on by \"try and error\" methods, depending on designer intuition and experience. Also, a generic design of an electrothermomechanical MEMS is a complex task that needs multidisciplinary knowledge. Thus, the objective of this work is to develop a software for the design of multi-phase MEMS, electrothermomechanically actuated, using a method for systematic and generic design, such as Topology Optimization Method (TOM). Using a peak function material interpolation model, any number of materials can be considered without increasing the amount of design variables if compared to an optimization with only one material and void. Mechanisms actuated by electric tension were designed considering the maximization of output displacement against a work piece with known stiffness. The design of microactuators considering multiflexibility was also performed. A study of optimization parameters influence is presented. As an alternative to electrothermomechanical actuation, some mechanisms actuated by heat flow were designed.