Neste trabalho os problemas de determinação dos Fatores de Intensidade de Tensão KI e KII estáticos e dinâmicos são tratados numericamente utilizando uma formulação alternativa do Método dos Elementos de Contorno (MEC) com solução fundamental de Kelvin e matriz de massa para os problemas dinâmicos. A trinca é suposta retangular inicialmente, com suas faces não-coincidentes. Tanto as faces da trinca, quanto o contorno externo são discretizados em elementos de contorno reto com variação de forças de deslocamentos quadráticas, não havendo, portanto distinção entre elementos de trinca e de contorno externo. Integrais analíticas também são obtidas para o elemento linear isoparamétrico. As células de domínio apresentam formato triangular e suas integrais são solucionadas semi-analiticamente. Quanto às integrais de contorno, essas são obtidas analiticamente segundo eixos de referência locais, procedendo-se em seguida a rotação pra eixos globais. O algoritmo de Houbolt é empregado como integrador temporal. Exemplos numéricos da determinação desses Fatores de Intensidade de Tensão são mostrados e comparados com resultados analíticos e resultados numéricos disponíveis na literatura. / In this work the stress intensity factors KI and KII for static and dynamic two-dimensional problem are obtained numerically by an alternative mass matrix boundary element formulation. The crack is considered a rectangular hole inside the domain and its faces are not coincident. Both crack faces and boundary are discretized by straight boundary elements with quadratic approximation. Domain cells are triangular with linear approximation and their integrals are developed semi-analytically. Boundary integrals are analytically performed, for linear and quadratic approximations. They are performed at local co-ordinate axes and transformed to global co-ordinate axes. The Houbolt algorithm is used to integrate the matrix time differential equation along time. Numerical examples are shown in order to compare the results obtained by the proposed formulation and the ones presents in literature.

O presente trabalho trata da análise do comportamento mecânico da madeira, material anisotrópico, submetida ao crescimento de fissuras e aos processos mecânicos que causam sua ruptura, utilizando os conceitos de Mecânica da Fratura. A avaliação do colapso de peças de madeira, pela da mecânica da fratura, foi realizada por meio de análises experimentais. Para a determinação das propriedades de fratura da madeira foram determinadas propriedades mecânicas da madeira, tais com: módulo de ruptura a flexão estática, módulo de elasticidade à flexão estática, resistência à tração paralela às fibras, resistência ao cisalhamento, resistência ao fendilhamento paralela às fibras e resistência à tração normal às fibras. Foram efetuados ainda ensaios experimentais referentes à mecânica da fratura em corpos de prova tipo SENB (Single Edge Notched Beam). A modelagem numérica foi realizada utilizando-se um código computacional baseado no método de elementos de contorno. Para as simulações foram utilizados os valores das propriedades físicas, de resistência e de elasticidade obtidas na caracterização experimental dos elementos de madeira. Os resultados experimentais foram comparados com as análises numéricas para determinação de uma formulação específica para situações de colapso de elementos de madeira. / The present work deals with the analysis of the mechanical behavior of wood, anisotropic material, submitted to the growth of cracks and the mechanical processes that cause its rupture, using the concepts of Fracture Mechanics. The evaluation of the collapse of pieces of wood, by the fracture mechanics, was carried out by means of experimental analyzes. To determine the fracture properties of the wood, mechanical properties of the wood were determined, such as: modulus of static bending, modulus of elasticity at static bending, tension resistance parallel to grain, shearing resistance, cleavage strength parallel to grain and tension resistance normal to grain. Experimental tests were performed on fracture mechanic in Single Edge Notched beams (SENB). Numerical modeling was performed using a computational code based on the boundary element method. For the simulations, the values of the physical properties, strength and elasticity obtained in the experimental characterization of the wood elements were used. The experimental results were compared with the numerical analyzes to determine a specific formulation for situations of collapse of wood elements.

A aplicação da análise inversa é objeto de estudo nos mais diversos campos da ciência e da engenharia. A motivação para o tratamento de tais problemas se deve ao fato de que em muitas aplicações dessas áreas do conhecimento, há a necessidade da identificação de parâmetros físicos e geométricos a partir de dados do domínio medidos experimentalmente, já que tais parâmetros de entrada são desconhecidos para uma análise direta do problema. Neste tipo de análise o problema principal está na quantidade e qualidade dos dados experimentais obtidos, que são na maioria das vezes insuficientes para garantir que o sistema gerado apresente solução única, gerando com isto um problema essencialmente mal-posto. Assim, de forma geral o emprego confiável da análise inversa implica na utilização de ferramentas eficientes de aquisição de dados experimentais aliada a técnicas numéricas de regularização que buscam a minimização da função objetiva gerada por algum método numérico, como por exemplo, o método dos elementos de contorno (MEC). Sendo assim, o presente trabalho tem por objetivo apresentar uma formulação para resolução de problemas inversos de valor de contorno e estimativa dos parâmetros do modelo coesivo, através de medidas de campos de deslocamentos, em sólidos bidimensionais com domínio formado por multi-regiões via (MEC), utilizando-se de técnicas tais como: mínimos quadrados, regularização de Tikhonov, decomposição em valor singular (SVD) e filtro de Tikhonov, para regularização do problema. Além disto, são apresentados alguns exemplos de aplicação da formulação desenvolvida. / The application of inverse analysis is nowadays subject of research of many fields in engineering and science. The motivation to consider this problem is due to the fact that in many applications of these knowledge areas, physical and geometric parameters, that are not directly known, can be identified using domain data measured experimentally. In this kind of analysis the main problem is the quantity and the quality of the obtained experimental data, which, many times, are not sufficient to guarantee that the generated system of equations has only one solution, leading therefore to an ill-posed problem. Thus, in general the reliable use of the inverse analysis requires using efficient tools for experimental data acquisition together with the numerical techniques of regularization needed to impose the minimization of the objective function written by using any numerical method, as the boundary element method (BEM) for instance. In this context, the objective of the present work is to derive a formulation for resolution of boundary-value inverse problems and to estimate the material parameters of the cohesive model, by using measured displacements fields, in multi-region two-dimensional solid by BEM, using techniques such as: least squares, Tikhonov regularization, singular value decomposition (SVD) and Tikhonov filtering, for the problem regularization. Some application examples are presented using the developed formulation to illustrate its performance.

Esse trabalho trata da análise de corpos multi-fraturados utilizando o método dos elementos de contorno. Este método numérico é conhecido por ser robusto e preciso neste tipo de problema e também por requerer pequeno esforço computacional na criação da malha de elementos para o crescimento das fissuras. Duas metodologias para a análise do comportamento das fissuras são consideradas. A primeira é a já consagrada metodologia dual. Por meio desta técnica equações integrais distintas são aplicadas às faces da fissura. Estas equações integrais são escritas em termos de deslocamentos e forças de superfície. A segunda metodologia é a que emprega a formulação singular onde a fissura é considerada como um vazio no domínio sendo as faces da fissura separadas por uma pequena distância. No tocante ao crescimento das fissuras foi desenvolvido um procedimento especial para a determinação da direção de crescimento das fissuras o qual mostrou-se muito eficiente levando a resultados precisos. O crescimento das fissuras é efetuado considerando o fator de intensidade de tensão atuante na extremidade de cada fissura. Dessa forma, as fissuras mais solicitadas apresentam maior comprimento de propagação tornando a análise mais realista. Os fatores de intensidade de tensão são calculados por meio de duas técnicas. A primeira é a já conhecida técnica de correlação de deslocamentos a qual relaciona os deslocamentos atuantes nas faces da fissura. Uma técnica alternativa é também utilizada a qual emprega o campo de tensões presente na extremidade da fissura. Após a determinação dos fatores de intensidade de tensão quatro diferentes teorias de interação de modos podem ser utilizadas para a determinação do ângulo de propagação. Foram analisadas estruturas sendo os resultados comparados aos previstos analiticamente e também numericamente. As respostas obtidas foram satisfatórias validando assim a metodologia proposta neste trabalho / This work deals with analysis of multi-fractured bodies using boundary element method. This numerical method is known to be robust and accurate in this kind of problem and by small computational effort to create elements mesh of crack growth. Two methodologies to analyze of crack behavior are considerate. The first is consecrated dual methodology. Through this technique different kind of integral equations are applied to crack boundaries. These integrals equations are written in displacements and traction variables. Second methodology is singular formulation. Through this technique crack is represented like a hole in body and the crack boundaries is separated by a small gap. For crack growth was created a special proceeding to determination crack growth direction. This method is very efficient and your results are accurate. Crack growth is made through the stress intensity factor performed in crack tip. Then the cracks more requested going to propagate with a larger length’s growth turning this model very realistic. The stress intensity factors are calculated through two techniques in this work. First is the known correlation displacement technique which related displacement in crack boundaries. An alternative technique is also used which consider stress field in crack tip. After determination of stress intensity factors four different theories are used to calculate the crack growth angle. In this work were analyzed structures with results are compared with analytical and numerical answers. The results obtained went very satisfactory validating the methodology proposed

Orientador: Isaias Vizotto / Tese (doutorado) - Universidade Estadual de Campinas, Faculdade de Engenharia Civil, Arquitetura e Urbanismo / Made available in DSpace on 2018-08-09T23:33:25Z (GMT). No. of bitstreams: 1 Rossi_ElianaMariadeMelloFrancisco_D.pdf: 5414575 bytes, checksum: a7297eb374b268f3bf484ab58acbe942 (MD5) Previous issue date: 2007 / Resumo: Neste trabalho é estudado o problema da elastodinâmica translento através da formulação do Método dos Elementos de Contorno (MEC) e da formulação combinada desse método com o Método dos Elementos Finítos (MEF). O estudo dos problemas governados por uma equação de campo vetorlal através do Método dos Elementos de Contorno, com a solução fundamental de Kelvin, resulta numa de domínio que tem sido tratada através de duas técnicas. A primeira desenvolve a Integração no domínio por células e permite determinar a contribuição Inercial dos nós do contorno e dos pontos internos. A segunda, conhecida como Reciprocidade Dual, transforma a de domínio em uma série de integrais de contorno. Foi desenvolvida uma tecnica para a manipulação das matrizes globais do sistema de equações obtido através do Método dos Elementos de Contorno, que permitiu resolver o problema da dependência linear de linhas e colunas decorrente da análise de regiões com cantos ou angulosidades. O algorítmo criado para condensação dos nós foi testado para o caso estático e dinâmico. Diversos exemplos são apresentados para os dois métodos e para o acoplamnento. Os resultados obtidos foram comparados com as analíticas disponíveis ou então com as soluções provenientes da aplicação de outros métodos numéricos / Abstract: In this work is studied the problem of olastodynamic transient through the farmulation of Boundary Element Method (BEM) and the formulation of this method combined with the Finite Elements Method The study of the problems governed by an oquation of vector fieid through the BEM with Kevin fundamental solution, comes up in domain that has been handled by two techniques. The first develops the integration in the domain by cells and determines the nodes contribution of the inertíal rame and internal points. The second, known as Dual Recíprocíty, transforms the domain integral in boundary integrais. It was developed a tochnique for handling the matrix of the system of equations obtained by BEM, whích It solves the problem of the linear dependence of rows and Golumns arlsing frem the analysls of regions wlth comers ar sharp boundaries. The algorithm created for node condensation has been tested for statle and dynamie cases. Several examples are presented using the both methods and the coupling technique. The results obtained were comparod with the anal.ytical solutions avaílable or wíth solutions obtained from other numerícal methods / Doutorado / Estruturas / Doutor em Engenharia Civil

Orientador: Paulo Sollero / Dissertação (mestrado) - Universidade Estadual de Campinas, Faculdade de Engenharia Mecanica / Made available in DSpace on 2018-08-10T03:16:07Z (GMT). No. of bitstreams: 1 Paiva_SeilaVastiFariade_M.pdf: 1387945 bytes, checksum: 04ad3c76f01698f24ef5fd7ebfb49264 (MD5) Previous issue date: 2006 / Resumo: Nesta dissertação é apresentada uma análise de problemas da mecânica da fratura elástica linear em estruturas bimateriais anisotrópicas. Utilizando o método dos elementos de contorno é possível calcular os fatores de intensidade de tensão em problemas planos (2D) devido à presença de trincas interfaciais entre as lâminas que compõem o material. A estrutura pode estar submetida à carregamento em modo I ou modo misto. O problema é modelado usando-se a técnica de sub-regiões para descrever cada um dos diferentes subdomínios, representado por cada material. Na interface das sub-regiões, em que o domínio é dividido, são impostas condições de equilíbrio de forças e continuidade de deslocamentos, exceto na região que corresponde à trinca. O comportamento singular apresentado pelo campo de tensões próximo à ponta da trinca é modelado com elementos de ponto a um quarto com singularidade de forças de superfície. São apresentados exemplos numéricos de problemas com carregamentos no plano. Foi também apresentada a análise de convergência de malhas, mostrando uma pequena dependência da discretização mesmo quando malhas pouco refinadas foram usadas. Alguns dos exemplos têm correspondentes na literatura, os quais foram utilizados para comparação com os resultados obtidos. Observou-se uma boa concordância na comparação dos resultados / Abstract: This thesis presents an analysis of problems of linear elastic fracture mechanics in anisotropic bimaterial structures. Using the boundary element method, it is possible to evaluate stress intensity factors in plane problems (2D) due to the presence of interfacial cracks between the laminae that constitute the material, when the structure is submitted to a mode I or in mixed mode loading. The problem is modeled using the subregion technique to describe each one of the different subdomains, represented by each material. On the interface of subregions, which the domain is divided, conditions of tractions equilibrium and displacements continuity are imposed, except in the corresponding crack region. The singular behavior presented by the stress field near the crack tip is modeled by traction singular quarter point element. Numerical examples of problems with in-plane loading are presented. Mesh convergence analyses are also presented, showing little dependence on the discretization even when coarse meshes were used. Some of these examples have correspondents in literature, that were used for comparisons with the obtained results. A good agreement in the comparisons of results was observed. / Mestrado / Mecanica dos Sólidos e Projeto Mecanico / Mestre em Engenharia Mecânica

Este trabalho apresenta o desenvolvimento e implementação computacional de técnicas de otimização de topologia para problemas governados pela equação de Poisson. O método numérico utilizado para solução numérica das equações foi o método dos elementos de contorno (MEC). Para tanto, três metodologias foram desenvolvidas. A primeira é direcionada à aplicação de algoritmos genéticos (AG) para investigar como um domínio inicialmente preenchido com cavidades aleatórias evolui durante um processo de otimização e verificar a possibilidade de se extrair topologias ótimas a partir da interpretação da solução encontrada. Os contornos externos permanecem fixos enquanto as posições e as dimensões das cavidades são otimizadas com o objetivo extremizar uma função custo especificada. O desempenho do algoritmo proposto é ilustrada com uma série de exemplos e os resultados são discutidos. A segunda metodologia apresenta um algoritmo numérico para otimização topológica baseado na avaliação da derivada topológica (DT), adotando a energia potencial total como função custo. Este procedimento é uma alternativa às tradicionais técnicas de otimização, evitando assim soluções de projeto com densidade de material intermediária. Sólidos com comportamento anisotrópico são estudados sob condições de contorno de Robin, Neumann e Dirichlet. Uma transformação linear de coordenadas é utilizada para mapear o problema original e suas condições de contorno para um novo domínio equivalente isotrópico, onde o procedimento de otimização é aplicado. A solução otimizada é então transformada de volta ao domínio original. A metodologia proposta mostrou-se particularmente atrativa para resolver esta classe de problemas já que o MEC dispensa o uso de malha no domínio, reduzindo significantemente o custo computacional. Na última parte deste trabalho foi implementada uma formulação de sensibilidade topológica para problemas de otimização de transferência de calor e massa simultâneos. Como as sensibilidades para cada equação diferencial são diferentes, utiliza-se um coeficiente de ponderação para compor a sensibilidade do problema acoplado. Isto permite a imposição de distintos fatores para cada problema, de acordo com uma prioridade especificada. Diversos exemplos são apresentados e seus resultados comparados com os da literatura, quando disponíveis, a fim de validar as formulações propostas. / This work presents the computational development and implementation of topology optimization techniques for problems governed by the Poisson equation. The boundary element method was the numerical technique chosen to solve the equations. Three different methodologies were developed aiming this objective. The first methodology is directed to the application of genetic algorithms to investigate how a domain previously populated with randomly placed cavities evolves during the optimization process, and to verify the resemblance of the final solution with a optimal design. The external boundaries remain fixed during the process, while the location and dimension of the cavities are optimized in order to extremize a given cost function. The performance of the proposed algorithm is verified with a number of examples and the results are discussed. The second methodology presents a numerical algorithm for topology optimization based on the evaluation of topological derivatives, using the total potential energy as the cost function. This procedure is an alternative to the traditional optimization techniques, avoiding design solutions containing intermediary material densities. Solids with anisotropic constitutive behavior are studied under Robin, Neumann and Dirichlet boundary conditions. A linear coordinate transformation approach is used to map the original problem into an isotropic one, where the optimization is carried out. The final solution is then mapped back to the original coordinate system. The proposed method was found to be an attractive way to solve this class of problems, since no interior mesh is necessary, which reduces significantly the computational cost of the analysis. In the last part of the present work the topological derivative approach was further developed to deal with the optimization of problems under simultaneous heat and mass transfer. Since the sensitivities for each differential equation are different, a weighting factor was used to evaluate the final sensitivities of the coupled problem. This allows the imposition of different priorities for each problem Several examples are presented and their results are compared with the literature, when available, in order to validate the proposed formulations.

Orientador: Paulo Sollero / Tese (doutorado) - Universidade Estadual de Campinas. Faculdade de Engenharia Mecanica / Made available in DSpace on 2018-08-06T09:16:06Z (GMT). No. of bitstreams: 1 Noritomi_PedroYoshito_D.pdf: 4114284 bytes, checksum: 806cfb9e9589c5cd848086ed60178957 (MD5) Previous issue date: 2005 / Resumo: Este trabalho mostra o desenvolvimento de uma metodologia para análise de problemas de bioengenharia, aplicando modelagem numérica elastostática de tensões e deformações, baseada no método dos elementos de contorno com formulação 3D para meios transversalmente isotrópicos lineares, incluindo a capacidade de simulação do comportamento de remodelagem óssea superficial. A implementação do núcleo transversalmente isotrópico sobre a estrutura básica de análise por elementos de contorno 3D utilizou a solução fundamental proposta por Pan & Chou e revisada por Loloi, tendo exigido o cálculo adicional das soluções fundamentais de força de superfície a partir da derivação das soluções fundamentais de deslocamento. O modelo de remodelagem óssea superficial baseou-se na hipótese de estímulo biológico por campo de deformação, partindo de um modelo 2D, adaptado para o espaço 3D com o uso de deformações principais como grandezas de referência. As implementações foram testadas através de análises numéricas de problemas com solução analítica e validações com resultados de aplicações comerciais baseadas em elementos finitos, para problemas padrão de engenharia, bem como comparações com resultados da literatura para problemas de bioengenharia. A análise dos resultados mostrará que, tanto a metodologia quanto as implementações são funcionais, oferecendo uma base sólida para desenvolvimento e teste de novas soluções de bioengenharia / Abstract: This work shows the development of a methodology to analyse bioengineering problems using elastostatic stress-strain numerical modeling based on a 3D transversely isotropic linear boundary element formulation including surface bone remodeling simulation capabilities. The transversely isotropic kernel implementation on the basic 3D boundary element analysis program used the fundamental solution purposed by Pan & Chou and revised by Loloi, with additional fundamental solutions for traction calculation made with the displacement fundamental solution derivatives. The surface bone remodeling model was based on a 2D strain field biological stimulus, extended to the 3D space by using the principal strain as reference values. The implementations were tested through numerical analysis of problems with analytical solution and validation with commercial finite elements applications for standard engineering problems, as well as comparison with literature data for bioengineering problems. The analysis of results will show that both, the methodology and the implementations are fully functional, offering a solid start for development and test of new bioengineering solutions / Doutorado / Mecanica dos Sólidos e Projeto Mecanico / Doutor em Engenharia Mecânica

Orientadores: Eder Lima de Albuquerque, Alberto Luiz Serpa / Dissertação (mestrado) - Universidade Estadual de Campinas, Faculdade de Engenharia Mecanica / Made available in DSpace on 2018-08-06T09:57:35Z (GMT). No. of bitstreams: 1 Gouvea_AdrianadosReis_M.pdf: 936019 bytes, checksum: 2509e5d6c5977816a1830afe11bcd691 (MD5) Previous issue date: 2006 / Resumo: Neste trabalho é apresentada uma análise, através do uso do método dos elementos de contorno, de problemas estáticos em um modelo de estado plano para materiais anisotrópicos. A formulação de elementos de contorno é desenvolvida usando soluções fundamentais para elasto-estática e são utilizados apenas elementos quadráticos contínuos (3 nós por elementos). A formulação desenvolvida é aplicada no cálculo de tensões e deslocamentos de estruturas planas de materiais compósitos laminados submetidos a cargas no plano da estrutura. São estudados também os critérios de falha para materiais compósitos, abordando dois critérios quadráticos: o critério de Tsai-Hill e o critério de Tsai- Wu. São utilizados métodos de otimização para encontrar os valores ótimos dos ângulos de orientação das fibras que maximizem a resistência da estrutura. Para a otimização de problemas de uma variável é usado o método da Seção Áurea e para problemas de várias variáveis é usado o método Simplex Modificado. A função objetivo para os dois casos é o critério de falha. As variáveis de projeto são os ângulos de orientação das fibras nas diversas lâminas do laminado. Os resultados obtidos são comparados com resultados disponíveis na literatura e mostram uma boa concordância / Abstract: This work presents an analysis, using the boundary element method, of static problems in a plain state model in anisotropic materiaIs. The boundary element formulation is developed using basic solutions for elasto-statics and only continuous quadratic elements (three nodes per element) are used. The developed formulation is applied to compute the stress and displacements of laminated composite structures plain submitted to loads in the structure plane. It is studied also the failure criteria for composite materials, considering two quadratic criteria: the Tsai-Hill criterion and the Tsai-Wu criterion. Optimization methods are used to find optimum values for the fibres orientation angles that maximize the strength of the structure. For optimization problems of one variable, the Golden Section method is used and for problems of several variables the Modified Simplex method is used. The objective function for both cases is the failure criteria. Design variables are fibres orientation angles of all lamina of the laminate. Results are compared with other results from literature and show a good agreement / Mestrado / Mecanica dos Sólidos e Projeto Mecanico / Mestre em Engenharia Mecânica

Orientadores: Euclides de Mesquita Neto, Edson Antonio Capello Sousa / Dissertação (mestrado) - Universidade Estadual de Campinas, Faculdade de Engenharia Mecanica / Made available in DSpace on 2018-08-09T04:23:51Z (GMT). No. of bitstreams: 1 Thomazo_LuizHenrique_M.pdf: 1773463 bytes, checksum: e38435e9b18935ddc2edc19b8b2efac4 (MD5) Previous issue date: 2004 / Resumo: Este trabalho apresenta uma metodologia para a realização de análises dinâmicas estacionárias e transientes em domínios visco-elásticos limitados ou ilimitados, utilizando-se o Método dos Elementos de Contorno. São utilizados como estados auxiliares soluções não-singulares de problemas do semi-espaço e do espaço completo de meios visco-elásticos. Os estados visco-elásticos auxiliares são sintetizados numericamente, tanto no domínio da freqüência, como no domínio do tempo. Os mencionados estados auxiliares são constituídos das soluções em tensão e deslocamento resultantes de cargas espacialmente constantes, aplicadas sob segmentos lineares tanto na superfície de semi-espaços, como no interior de espaços completos. Inicialmente os estados auxiliares são sintetizados numericamente no domínio da freqüência, dando origem à análise de problemas estacionários. Na seqüência estados auxiliares transientes são obtidos a partir da aplicação da transformada de Fourier rápida (FFT) sobre os estados estacionários mencionados. Os estados auxiliares são utilizados para a síntese de uma versão não singular Indireta do MEC (MEC-I) tanto no domínio da freqüência quanto no domínio do tempo. No trabalho ainda é formulada e implementada a inclusão de um corpo rígido interagindo com o meio discretizado pelo MEC-I. O comportamento visco-elástico do contínuo é introduzido pelo princípio da correspondência e são analisados os modelos de histerése constante e de Kelvin-Voigt / Abstract: The main purpose of the present work is to further develop a methodology to perform stationary and transient dynamic analysis of viscoelastic continua by the Boundary Element Method. Numerically synthetized stationary and transient, half-space and full-space auxiliary states are employed to render a non-singular implementation of the indirect version of the Boundary Element Method, IBEM. The auxiliary states used in this word are displacement and traction components stemming from the solution two stress bounary value problems. A constant traction distribution applied at a linear segment over the viscoelastic half-space surface and at the interior of the viscoelastic full-space describe the solutions applied to formulate the BEM. The IBEM is formulated and implemented to render the dynamic solution of bounded and unbounded domains. The interaction of the IBEM mesh with a massless rigid body is also formulated in the context of this work. Viscoelastic effects are included by means of the correspondence principle. The influence of the constant hysteresis and the Kelvin-Voigt damping models are studied. / Mestrado / Mecanica dos Sólidos e Projeto Mecanico / Mestre em Engenharia Mecânica