Quando o concreto é submetido a ações especiais, como cargas cíclicas ou ação de cargas de impacto, modificações em sua composição são necessárias, já que o concreto não apresenta desempenho satisfatório à tração, o que compromete o seu comportamento frente à ação de cargas dinâmicas. Uma alternativa para amenizar esta deficiência consiste em adicionar fibras ao concreto. Estas atuam como reforços à tração, transformando a matriz cimentícia, tipicamente frágil, em um material que apresenta boa resistência residual após a fissuração. Buscando colaborar na avaliação da eficiência de diferentes tipos de fibras, o presente trabalho analisa o comportamento de concretos com fibras de aço, submetidos ao impacto, avaliando a influência do fator de forma, do comprimento e do teor de fibras, assim como do tamanho do agregado. São ainda analisados os efeitos da incorporação de fibras na resistência à compressão, na resistência à tração por compressão diametral, no módulo de elasticidade e na tenacidade dos compósitos. Adicionalmente, é executada uma comparação entre os resultados experimentais e os derivados de um esquema de modelagem da situação de impacto através do uso do método de elementos discretos. Buscou-se, através da modelagem teórica, executar um mapeamento dos danos, provocados por cargas de impacto incrementais, ao longo do tempo, bem como determinar as energias necessárias para levar as placas até a ruptura. Os resultados indicam que a incorporação de fibras de aço não consegue retardar o aparecimento da primeira fissura, mas aumenta significativamente a tenacidade dos compósitos. Fibras mais longas e com maior fator de forma tendem a ser mais eficientes, desde que se supere um teor de fibras mínimo, que neste trabalho ficou em torno de 100.000 fibras/m3, para fibras longas (50-60 mm) e de 400.000 fibras/m3, para fibras curtas, cuja ancoragem é menos eficiente. O método de teste de impacto por queda de esfera se mostrou adequado e sensível, porém o esquema de modelagem numérica testado necessita ser refinado para permitir uma adequada simulação do comportamento de concretos com fibras. / When submitted to special loading patterns, derived from dynamical actions such as cyclic or impact loads, concrete elements need to be reinforced, in order to resist the tensile stresses. A feasible alternative, in such cases, is to incorporate fibers in the concrete matrix. The fibers act as a tensile reinforcement, transforming the fragile cement matrix into a composite with significant post-cracking residual strength. In order to contribute with the data collection about the efficiency of different fiber types, the present research work presents an analysis of the behavior of steel fiber reinforced concretes subjected to impact loads. The work investigates the influences of changes in the shape factor, length and amount of fibers. The effects of these combinations on other basic properties of the composites, such as compression strength, split cylinder tensile strength, Young’s modulus and tenacity is also measured. Additionally, a comparison is made between the experimental results from the impact tests and the estimates obtained from a theoretical model that uses the discrete element method (DEM). This theoretical approach aimed to determine if the model was able to simulate the damage evolution over time generated by the increasing impacts loads, as well as to determine the total energy necessary to crack and break the specimens. The results obtained pointed out that the introduction of steel fibers does not affect the energy for the first crack but increases significantly the tenacity of the composite. Longer fibers, with greater shape factors, tend to be more efficient, provided that the fiber content is sufficiently high. The minimum recommended fiber content, according to the data from this research, may be around 100.000 fibers/m3, for longer fibers (50-60 mm). Or around 400.000 fibers/m3, for shorter fibers, which are not so efficient in terms of anchorage. The impact test method developed was considered adequate, being sensitive to the phenomenon and providing reliable data. The DEM model, however, needs to be refined to be able to deal with fiber concrete composites.

Quando o concreto é submetido a carregamentos especiais, como cargas cíclicas ou ação de cargas de impacto, modificações em sua composição são necessárias. Uma vez que o material não apresenta desempenho satisfatório à tração, seu comportamento frente a este tipo de carregamento acaba seriamente comprometido. Uma alternativa para amenizar esta deficiência consiste em adicionar fibras de aço ao concreto. Ao adicionar estes elementos à matriz cimentícia, promove-se meios de transferência de tensões através das fissuras, aumentando a tenacidade do material, proporcionando mecanismos de absorção, relacionados com o desligamento e o arrancamento de fibras. Um número significativo de trabalhos experimentais envolvendo os mais diversos tipos de elementos estruturais reforçados com fibras de aço está disponível, havendo, no entanto, uma forte carência sob o ponto de vista de simulações numéricas. Buscando colaborar no desenvolvimento do material, o presente trabalho propõe a aplicação do Método dos Elementos Discretos para simulação do compósito submetido a carregamentos estáticos e dinâmicos. São realizadas alterações no algoritmo do método a fim de realizar a dispersão de fibras de aço na matriz de concreto. A análise das condições de contorno utilizadas em trabalho anterior revela a necessidade de aplicação de apoios elásticos sob pena de superestimar a rigidez do modelo. Os diagramas carga versus deslocamento que resultaram dos ensaios estáticos demonstram que o modelo criado é sensível à adição de fibras: maiores teores conduzem a modelos com maior tenacidade. O ensaio de impacto também se mostrou sensível e o padrão de fissuração encontrado nas simulações revelou uma boa aproximação com ensaios experimentais anteriores. / When submitted to special loading patterns, derived from dynamical actions such as cyclic or impact loads, some alterations in the concrete constitution need to be done, since the material don’t have an adequate behavior under tensile stress. A feasible alternative, in such cases, is to incorporate steel fibers in the concrete matrix. Adding these elements, stress transference mechanisms along the cracks are promoted, increasing the material tenacity. An expressive number of experimental works involving all the kinds of steel fiber reinforced concrete structural elements are available. However, few researches based on numerical methods are found in the literature. In order to contribute with the data collection and the development of the material, the present research work proposes the application of the Discrete Element Method to simulate the composite subjected to static and dynamic loads. Some modifications are made on the method algorithm trying to create the dispersion of fibers in the concrete matrix. The analysis of the boundary conditions used on previous work reveal the importance of using elastic support to don’t overestimate the stiffness of the model. The diagram load versus displacement that came from the static simulations shows that the model is sensible to the addition of fibers: higher proportions of fiber leads to models with higher tenacity. The impact tests also demonstrate sensibility and the crack pattern found on the simulations presented a very good approximation to previous experimental work.

O objetivo deste trabalho é estudar e desenvolver procedimentos numéricos adequados para a análise de problemas da Engenharia do Vento Computacional (EVC). O escoamento é analisado a partir das equações de Navier-Stokes para um fluido Newtoniano e de uma equação de conservação de massa considerando a hipótese de pseudo-compressibilidade, ambas em um processo isotérmico. Na presença de escoamentos turbulentos emprega-se a Simulação de Grandes Escalas (“LES”) com os modelos clássico e dinâmico de Smagorinsky para as escalas inferiores à resolução da malha. Dois modelos numéricos de Taylor-Galerkin para a análise do escoamento são estudados: o esquema explícito de dois passos e o esquema explícito-iterativo. O Método dos Elementos Finitos (MEF) é empregado para a discretização do domínio espacial utilizando o elemento hexaédrico trilinear isoparamétrico com integração reduzida das matrizes em nível de elemento. Em problemas envolvendo efeitos de interação fluido-estrutura emprega-se um esquema de acoplamento particionado com características superiores de conservação, permitindo, inclusive, o uso de subciclos entre as análises do fluido e da estrutura e de malhas não compatíveis na interface. A estrutura é considerada como um corpo deformável constituído de um material elástico linear com a presença de nãolinearidade geométrica. O MEF é também usado para a discretização da estrutura, empregando-se para tanto o elemento hexaédrico trilinear isoparamétrico com integração reduzida e controle de modos espúrios. A equação de equilíbrio dinâmico é integrada no tempo utilizando o método implícito de Newmark no contexto do método de estabilização α- Generalizado. Na presença de estruturas deformáveis, o escoamento é descrito através de uma formulação arbitrária Lagrangeana-Euleriana (ALE). Ao final, comparações com exemplos numéricos e experimentais são apresentadas para demonstrar a viabilidade dos algoritmos desenvolvidos, seguindo-se com as conclusões do trabalho e as sugestões para trabalhos futuros. / Analysis and development of numerical tools to simulate Computational Wind Engineering (CWE) problems is the main goal of the present work. The isothermal flow is analyzed using the Navier-Stokes equations for viscous fluids and a mass conservation equation obtained according to the pseudo-compressibility assumption. Turbulent flows are simulated employing Large Eddy Simulation (LES) with the classical and dynamic Smagorinsky’s models for subgrid scales. Two Taylor-Galerkin models for the flow analysis are investigated: the explicit two-step scheme and the explicit-iterative scheme. The Finite Element Method (MEF) is employed for spatial discretizations using the eight-node hexahedrical isoparametric element with one-point quadrature. Fluid-structure interaction problems are analyzed with a coupling model based on a conservative partitioned scheme. The Finite Element Method (MEF) is employed for spatial discretizations using the eight-node hexahedrical isoparametric element with one-point quadrature. Fluid-structure interaction problems are analyzed with a coupling model based on a conservative partitioned scheme. Subcycling and nonmatching meshes for independent discretizations of the fluid and structure domains are also available. The structure is considered as a deformable body constituted by a linear elastic material with geometrically nonlinear effects. The FEM is used for the spatial discretization of the structure as well. Eight-node hexahedrical isoparametric elements with one-point quadrature and hourglass control are adopted in this process. The implicit Newmark algorithm within the framework of the α-Generalized method is employed for the numerical integration of the dynamic equilibrium equation. An arbitrary Lagrangean-Eulerian (ALE) description is adopted for the kinematic description of the flow when deformable structures are analyzed. Numerical and experimental examples are simulated in order to demonstrate the accuracy of the developed algorithms. Concluding remarks and suggestions for future works are pointed out in the last chapter of the present work.

O principal objetivo do controle de vibrações é a sua redução ou minimização, através da modificação automática da resposta estrutural. Em muitas situações isto é necessário para promover a estabilidade estrutural, e para alcançar o alto desempenho mecânico necessário em diversas áreas técnicas, tais como a engenharia aeroespacial, civil e mecânica, bem como a biotecnologia, inclusive em escala micro e nano mecânica. Uma alternativa é o uso de estruturas inteligentes, que são o resultado da combinação de sensores e atuadores integrados em uma estrutura mecânica, e um método de controle adequado. O principal objetivo deste trabalho é o desenvolvimento de rotinas computacionais para a simulação, via método dos elementos finitos, do controle ativo de estruturas inteligentes de cascas, placas e vigas delgadas de material compósito laminado com camadas de material piezoelétrico como sensores e/ou atuadores. Caracterizam esta pesquisa a utilização do elemento GPL-T9 de três nós e seis graus de liberdade mecânicos por nó, mais um grau de liberdade elétrico por camada piezoelétrica, assim como a avaliação de dois métodos de controle, o Proporcional-Integral-Derivativo (PID) e o Regulador Quadrático Linear ou Linear Quadratic Regulator (LQR), incluindo o LQR Modal, e a otimização da localização de pastilhas piezoelétricas através de um Algoritmo Genético (AG). Várias aplicações são apresentadas e os resultados obtidos são comparados aos disponíveis na literatura especializada. / The main objective of vibration control is its reduction or even its minimization by the automatic modification of the structural response. Sometimes this is necessary to increase structural stability and to attain a high mechanical behavior in several areas such as aerospace, civil and mechanical engineering, biotechnology, including macro, micro and nanomechanical scales. An alternative is to use a smart structure, which results of the combinations of integrated sensors and actuators in a mechanical structure and a suitable control method. Development of a computational code to simulate, using finite elements, the active control in smart structures such as slender shells, plates and beams of composite materials with embedded piezoelectric layers acting as actuators and sensors is the main objective of this work. This research is characterized by the use of the GPL-T9 element with three nodes and six mechanical degrees of freedom and one electrical degree of freedom per piezoelectric layer, by the evaluation of two control methods, the Proportional Integral Derivative (PID) and the Linear Quadratic Regulator (LQR), including the Modal LQR, and, finally by the optimization of piezoelectric patches placement using a Genetic Algorithm (GA). Several examples are presented and compared with those obtained by other authors.

O objetivo deste trabalho é estudar algoritmos paralelos para a solução de problemas de Mecânica dos Sólidos, Mecânica dos Fluídos e Interação Fluido-Estrutura empregando o Método dos Elementos Finitos para uso em configurações de memória distribuída e compartilhada. Dois processos para o particionamento da estrutura de dados entre os processadores e divisão de tarefas foram desenvolvidos baseados na aplicação do método de particionamento em faixas e do método da bissecção coordenada recursiva não sobre a geometria da malha mas sim diretamente sobre o sistema de equações, através de reordenações nodais para minimização da largura da banda. Para ordenar a comunicação entre os processadores, foi desenvolvido um algoritmo simples e genérico baseado em uma ordenação circular e alternada que permite a organização eficiente dos processos mesmo em cenários nos quais cada processador precisa trocar dados com todos os demais. Os algoritmos selecionados foram todos do tipo iterativo, por sua adequabilidade ao paralelismo de memória distribuída. Foram desenvolvidos códigos paralelos para o Método dos Gradientes Conjugados utilizado em problemas de Mecânica dos Sólidos, para o esquema explícito de Taylor-Galerkin com um passo e iterações utilizado na simulação de escoamentos compressíveis em regime transônico e supersônico, para o esquema explícito de Taylor- Galerkin com 2 passos para simulação de escoamentos incompressíveis em regime subsônico e para interação fluído-estrutura usando o esquema explícito de dois passos para o fluído e o método implícito de Newmark no contexto do método de estabilização α-Generalizado para a estrutura, com acoplamento particionado. Numerosas configurações foram testadas com problemas tridimensionais utilizando elementos tetraédricos e hexaédricos em clusters temporários e permanentes, homogêneos e heterogêneos, com diferentes tamanhos de problemas, diferentes números de computadores e diferentes velocidades de rede. / Analysis and development of distributed memory parallel algorithms for the solution of Solid Mechanics, Fluid Mechanics and Fluid-Structure Interaction problems using the Finite Element Method is the main goal of this work. Two process for mesh partitioning and task division were developed, based in the Stripwise Partitioning and the Recursive Coordinate Bisection Methods, but applied not over the mesh geometry but over the resultant system of equations through a nodal ordering algorithm for system bandwidth minimization. To schedule the communication tasks in scenarios where each processor must exchange data with all others in the cluster, a simple and generic algorithm based in a circular an alternate ordering was developed. The algorithms selected to be parallelized were of iterative types due to their suitability for distributed memory parallelism. Parallel codes were developed for the Conjugate Gradient Method ( for Solid Mechanics analysis), for the explicit one-step scheme of Taylor-Galerkin method (for transonic and supersonic compressible flow analysis), for the two-step explicit scheme of Taylor-Galerkin method (for subsonic incompressible flow analysis) and for a Fluid-Structure Interaction algorithm using a coupling model based on a partitioned scheme. Explicit two-step scheme of Taylor-Galerkin were employed for the fluid and the implicit Newmark algorithm for the structure. Several configurations were tested for three-dimensional problems using tetrahedral and hexahedral elements in uniform and nonuniform clusters and grids, with several sizes of meshes, numbers of computers and network speeds.

A otimização estrutural, utilizando ferramentas computacionais é um grande campo de pesquisa na atualidade. Os métodos utilizados, dependendo da complexidade do problema, demandam um grande custo computacional, e por isso vem sendo avaliandas várias técnicas para diminuí-lo. Uma delas é o emprego de técnicas de aproximação de análises, dentre as quais destacam-se as redes neurais, que combinadas aos métodos de otimização e de análises clássicos conseguem bons resultados e reduzem significativamente o tempo de processamento. O emprego dos compósitos laminados como material estrutural vem crescendo nos últimos tempos, incentivado pela suas excelentes propriedades mecânicas e baixo peso. Em consenso com todo o esforço científico dedicado a essa área, o presente trabalho visa a implementação de uma ferramenta computacional capaz de otimizar estruturas complexas fabricadas com tais materiais, a um baixo custo computacional. Com isto em mente, é desenvolvido um sistema de otimização, aproveitando módulos implementados previamente para a análise estática linear e não linear através do método dos elementos finitos (MEF), e o módulo de otimização por algoritmos genéticos. Serão desenvolvidos os módulos de análise modal, para otimizar também estruturas com critérios baseados em freqüências e modos, e o modulo de redes neurais de tipo perceptron para aproximações das análises feitas através do MEF. Alguns exemplos são apresentados para demonstrar que bons resultados são obtidos com a utilização de redes neurais artificiais, cujo treinamento permite poupar tempo computacional proveniente do grande número de análises usualmente necessárias no processo de otimização. / Structural optimization using computational tools has become a major research field in recent years. Methods commonly used in structural analysis and optimization may demand considerable computational cost, depending on the problem complexity. Therefore, many techniques have been evaluated in order to diminish such impact. Among these various techniques, artificial neural networks may be considered as one of the main alternatives, when combined with classic analysis and optimization methods, to reduce the computational effort without affecting the final solution quality. Use of laminated composite structures has been continuously growing in the last decades due to the excellent mechanical properties and low weight characterizing these materials. Taken into account the increasing scientific effort in the different topics of this area, the aim of the present work is the formulation and implementation of a computational code to optimize manufactured complex laminated structures with a relatively low computational cost by combining the Finite Element Method (FEM) for structural analysis, Genetic Algorithms (GA) for structural optimization and Artificial Neural Networks (ANN) to approximate the finite element solutions. The modules for linear and geometrically non-linear static finite element analysis and for optimize laminated composite plates and shells, using GA, were previously implemented. Here, the finite element module is extended to analyze dynamic responses to optimize problems based in frequencies and modal criteria, and a module with perceptron ANN is added to approximate finite element analyses. Several examples are presented to show the effectiveness of ANN to approximate solutions obtained using the FEM and to reduce significatively the computational cost.

Diversos materiais de interesse tecnológico podem ser considerados como materiais heterogêneos, onde sua natureza aleatória deve ser considerada para representar corretamente o comportamento não-linear. A avaliação quantitativa do dano em materiais sujeitos a estados de tensão ou deformação tem grande importância devido ao caráter crítico desse fenômeno, que pode crescer de forma abrupta e resultar em uma falha catastrófica da estrutura. Em estudos anteriores, Carpinteri e seus colaboradores apresentaram diferentes aspectos relacionados à caracterização do dano em materiais heterogêneos. Três desses aspectos são discutidos neste trabalho: (i) a proposta do número de fragilidade como medida da fragilidade da estrutura em análise; (ii) a evolução da dimensão fractal em que o processo de dano é desenvolvido; (iii) e os índices globais obtidos pela análise de Emissão Acústica (EA). Neste trabalho, uma versão do método dos elementos discretos formado por barras é utilizada para explorar esses aspectos. Espécimes de material quase-frágil são simulados e, quando possível, os resultados numéricos são comparados com resultados experimentais. Além disso, a discussão dos resultados obtidos é realizada visando compreender o comportamento desse tipo de material, como por exemplo, o fato de que estruturas com diferentes dimensões, porém com o mesmo número de fragilidade, têm comportamentos similares. Ademais, o método numérico empregado é apresentado como uma ferramenta viável para complementar a informação obtida em ensaios experimentais na avaliação do processo de dano. Pelo sinal de EA obtido através do método numérico empregado, obtiveram-se parâmetros para a classificação do tipo de falha e para o cálculo do tensor momento, onde comparando os resultados obtidos pelo modelo numérico e pela EA se observam comportamentos coerentes. / Several materials with technological interest can be considered as heterogeneous materials, where their random nature must be accounted to correctly represent the nonlinear behavior. The quantitative evaluation of damage in materials subjected to stress or strain states have great importance due to the critical nature of this phenomenon, which abruptly can rise to catastrophic failure. In previous studies, Carpinteri and his coworkers have presented different aspects of the damage process characterization in heterogeneous materials. Three of these aspects are discussed in this work: (i) the brittleness number proposal to measure the brittleness level of the analyzed structure, (ii) the relationship of the fractal dimension in which the damage process is developed, (iii) and the global indexes obtained by the Acoustic Emission (AE) analysis. In the present work, a version of discrete element method formed by bars is used to explore these concepts. A set of quasi-brittle material specimen is simulated and, when possible, the numerical results are compared with experimental ones. Moreover, a discussion of the obtained results is carried out aiming to understand the behavior of this kind of material, for instance, the fact that structures with different dimensions, but with the same brittleness number, have similar behavior. In addition, the numerical method is presented as a viable tool to complement information from experimental test on the damage process. From the AE signal obtained by the numerical method, parameters to classify the type of crack and for calculating the moment tensor were obtained, where consistent behaviors are observed by comparing the results of the numerical model and the AE.

Esta dissertação tem por objetivo a implementação de um método adaptativo anisotrópico de malhas e seu acoplamento a um programa de simulação de escoamentos de fluidos compressíveis numa faixa ampla de número de Mach utilizando o Método dos Elementos Finitos (Bubnov-Galerkin). Mais especificamente, implementa-se um código para simulação de escoamentos compressíveis utilizando o método das linhas ou direções características (CBS), o qual possui boa performance para uma amplo espectro de números de Mach. A discretização espacial é feita utilizando-se elementos finitos triangulares lineares e tetraédricos lineares para problemas bidimensionais e tridimensionais, respectivamente. A discretização temporal é feita por uma expansão em séries de Taylor. Implementa-se uma adaptação anisotrópica com capacidade de refinamento e desrefinamento do domínio computacional para se obter soluções precisas a baixos custos computacionais, onde a estimativa de erro é conduzida por uma métrica Riemanniana, a qual permite o tratamento contínuo do erro através de um tensor de erro anisotrópico. O código é modificado para possibilitar o uso de arquiteturas de memória compartilhada a fim de aumentar a velocidade de processamento da simulação. Por fim, diversos estudos, tanto para problemas bidimensionais quanto tridimensionais, são realizados, validando a metodologia e realizando estudos adicionais em problemas mais complexos. / The objectives of this dissertation are the implementation of an anisotropic mesh adaptation procedure and its coupling with a compressible flow simulation code with large Mach range using the Finite Element Method (Bubnov-Galerkin). Specifically, a computer code using the Characteristic Based Split method called CBS, which has a good performance for a wide range of Mach numbers, is implemented. The spatial discretization is performed using linear triangular and linear tetrahedral finite elements for two and three dimensional cases, respectively. The temporal discretization is done with Taylor series. An anisotropic mesh adaption technique, based on a Riemann metric to evaluate the errors (which allows continuous error analysis through an anisotropic error tensor evaluation), including mesh refinement and mesh coarsening procedures is also implemented looking for accurate and computationally cheaper solutions. The code is able to use shared memory architecture in order to reduce computer process time. Finally, two-dimensional and three-dimensional studies are performed in order to validate the proposed methodology and additional complex tests are also analysed.

Esse trabalho trata da análise computacional de problemas aeroelásticos. Os casos de interesse envolvem escoamentos compressíveis em torno de estruturas de cascas de materiais compósitos laminados. A solução do problema de interação fluido- estrutura é conduzida por meio do esquema particionado denominado improved serial staggered procedure (ISS), o que permite a utilização de algoritmos independentes para a análise de cada componente do sistema. Um elemento triangular plano com três nós destinado à análise de cascas de materiais compósitos laminados é obtido com base na formulação de dois elementos de alta performance desenvolvidos para a análise de membranas e de placas. A flexibilidade ao corte transversal, a ausência de travamento por razão de aspecto em problemas com flexão no plano da casca são características do elemento. A análise dinâmica não linear é realizada através do procedimento corrotacional de conservação aproximada de energia (AECCP) que foi implementado com base na formulação corrotacional independente de elemento (element independent CR formulation – EICR). Essa abordagem permite que problemas transientes com grandes rotações e forte não linearidade geométrica sejam integrados por longos períodos de tempo, mantendo a estabilidade, a precisão e a eficiência da solução. A simulação do escoamento é realizada através do método explícito de Taylor-Galerkin de dois passos utilizando elementos finitos tetraédricos em malhas não estruturadas para a discretização espacial. As equações governantes do escoamento são dadas segundo a formulação lagrangeana-euleriana arbitrária (ALE), permitindo que malhas móveis sejam empregadas na análise dos problemas aeroelásticos. Um método simples e eficiente é adotado para a movimentação da malha do fluido com base na distância do nó às superfícies de contorno do escoamento. A utilização de malhas não coincidentes na discretização do fluido e da estrutura é possibilitada pelo emprego do esquema de projeção nodal para a transferência de informações entre os meios na região de interface. A apresentação, implementação e verificação de cada componente da ferramenta de análise é realizada independentemente graças à natureza do esquema de solução particionado. Por fim são analisados problemas de interação fluido-estrutura, onde as respostas obtidas são comparadas com resultados experimentais e numéricos de outros autores. Também são propostos e analisados novos problemas envolvendo estruturas de compósitos laminados. / This work deals with computational analysis of aeroelastic problems, specifically those related to compressible flows around laminated composite shell structures. The fluid- structure problem is solved by a partitioned scheme called improved serial staggered procedure (ISS), allowing to use independent algorithms for the analysis of each system component. A 3-node triangular flat element for laminated composite shells is obtained based on the formulation of two high performance plate and membrane elements. Shear flexibility and no aspect ratio locking on in-plane bending problems are some of the element characteristics. The nonlinear dynamic analysis is performed using the approximately energy conserving corotational procedure (AECCP), which is implemented based on the element independent corotational formulation (EICR). This approach enables long time transient problems with strong nonlinearities and large rotations to be efficiently solved while keeping stability and accuracy. The flow simulation is performed by the explicit two step Taylor-Galerkin method employing tetrahedral finite elements in unstructured meshes for the space discrezation. The flow governing equations are given in arbitrary lagrangian-eulerian formulation, enabling aeroelastic simulations with moving meshes. A simple and efficient method is adopted to move the fluid mesh based on the distance of each node to the boundary surfaces. Transference of information between non-matching fluid and structure meshes are performed by the node-projection scheme. Each component of the aeroelastic analysis tool is independently presented, implemented and verified due to the partitioned nature of the scheme adopted for the solution of the coupled system. Some aeroelastic problems are analyzed and results are compared to other experimental and numerical works. New problems of aeroelastic analysis of laminated composite structures are proposed and solved.

A ruptura de materiais heterogêneos quase frágeis, como o concreto, as cerâmicas e diferentes tipos de rochas, tem um comportamento mecânico complexo, foco de estudo de pesquisadores já há muito tempo. Uma ferramenta para simular o comportamento até ruptura destes tipos de materiais são os métodos de elementos discretos formados por barras (Discrete Element Method - DEM). Neste método as massas são concentradas em pontos nodais que estão ligados por meio de elementos unidimensionais caracterizados por relações constitutivas uniaxiais simples. Neste contexto, realizam-se novas implementações na versão do DEM com o intuito de estudar diversos aspectos dos materiais quase frágeis. Os principais focos da tese estão na desvinculação do nível de discretização do modelo do comprimento de correlação dos campos aleatórios que caracterizam as propriedades mecânicas. Através da interpretação dos resultados na simulação de emissão acústica é possível conhecer melhor o processo de dano neste tipo de material. Testes experimentais preliminares também são apresentados. Também são realizadas simulações em DEM utilizando carregamento biaxial, onde é explorado nessas simulações o comportamento dos modelos sob a ótica da técnica de emissão acústica. / The rupture of heterogeneous materials quasi-fragile, like concrete, ceramics and different types of rocks have a complex mechanical behavior, which has been the focus of study by researchers for a long time. The truss-like Discrete Element Method (DEM) was used to perform numerical simulations of the testing processes. The test results and the results of the numerical analyses, in terms of load vs. time diagram and AE data, as determined through b-value. In this context, new implementations in the version of DEM were realized with the objective to study the aspects of quasi-fragile materials. The main focuses of the thesis are in untying the level of discretization of the model of the correlation length from random fields which characterize the mechanical properties. And, through the interpretation of the results in the simulation of the acoustic emission, it is possible to better understand the process of damage in this type of material. Are also performed in DEM the biaxial loading, where is explored in this simulations the behavior of models from the point of view of the acoustic emission technique.