Este trabalho descreve um modelo matemático para simulação da laringe humana durante a fonação. O objetivo foi produzir uma técnica computacional de grande escala de processamento para capturar os fenômenos fisiológicos que ocorrem na laringe durante a vocalização e servir de base para estudos mais aprofundados sobre esta importante estrutura do corpo humano. Usando o método dos elementos finitos como base para discretizar as equações dos tecidos musculares da laringe e das equações de Navier-Stokes, e um modelo de descrição da colisão entre as pregas vocais, o sinal glotal foi obtido a partir de diferentes geometrias de laringes com diferentes propriedades viscoelásticas. Os resultados confirmaram a teoria mioelástica-aerodinâmica que descreve a dinâmica da fonação, reproduzindo inclusive fenômenos fisiológicos que os modelos existentes são incapazes de simular. Estudos adicionais foram feitos para verificar a viabilidade do modelo para simular algumas doenças que danificam a laringe. / This work describes a mathematical model to simulate the human larynx during a phonation. The objective was to produce a large-scale computational technique to capture several physiological phenomena that take place on the larynx during the vocalization and to assist further studies about this important structure of the human body. Using the finite element methods as the way to discretize the muscle tissue equations of the larynx and the Navier-Stokes equations and a model to describe the collision between both vocal folds, the glottal signal for different larynx geometries with different viscoelastic properties was obtained. The results confirmed the myoelastic-aerodynamic theory which describes the dynamic of the phonation, also reproducing physiologic phenomena that current models are unable to simulate. Additional studies were conducted to confirm the feasibility of the model to simulate some diseases that affect the larynx.

Neste estudo se propõe o desenvolvimento de um modelo numérico para a ligação deslizante entre elementos sólidos bidimensionais, aplicável à simulação de sistemas deslizantes de isolação de base para estruturas. A formulação implementada é baseada no Método dos Elementos Finitos Posicional (MEFP) para análise dinâmica não linear geométrica de estruturas escrita na forma Lagrangeana total. Elementos triangulares planos e isoparamétricos de aproximação cúbica com matriz de massa completa são utilizados principalmente na elaboração da parte sólida dos dispositivos de ligação entre estruturas reticuladas e a base móvel. Esses elementos também poderm ser utilizados na modelagem da estrutura em si, porém, para esse fim, elementos finitos isoparamétricos de barra geral com massa distribuída por unidade de comprimento foram implementados. As equações de movimento são integradas no tempo aplicando o método de Newmark e o problema de deslizamento é resolvido com o algoritmo baseado na técnica dos multiplicadores de Lagrange, onde a restrição das posições de um nó escravo é feita em relação a uma sequência de superfícies mestres. Elementos de barra geral foram usados para simular as superfícies mestres de contato, o que aumenta as possibilidades de aplicações, incluindo mecanismos compostos apenas por barras gerais. Analisam-se exemplos disponíveis na literatura para a validação da formulação proposta e propõem-se aplicações diversas na engenharia das estruturas. / This study proposes the development of a numerical model for the sliding joint between two-dimensional solid elements, applicable to the simulation of sliding base isolation systems. The implemented formulation is based on the Positional Finite Element Method (PFEM) for geometrical nonlinear dynamic analysis of structures written in the total Lagrangian form. Plane and isoparametric triangular cubic approximation elements with full mass matrix are mainly used in the elaboration of the solid part of the devices of joints between reticulated structures and mobile base. These elements can also be used in the modeling of the structure itself, however, for that purpose, isoparametric elements of general bar with mass distributed per unit of length were implemented. The motion equations are integrated in time by applying the Newmark method and the sliding problem is solved with the algorithm based on the technique of Lagrange multipliers, where the constraint of the positions of a slave node is made in relation to a sequence of master surfaces. General bar elements were used to simulate the master contact surface, which increases the possibilities of applications, including mechanisms composed only of general bars. Analyze examples available in the literature for the validation of the proposed formulation and proposed diverse applications in the engineering of the structures.

O presente trabalho consiste no uso do Método dos Elementos Finitos (MEF) para a análise de interação fluido-estruturas de barras com foco em problemas transientes envolvendo explosões ou outras ações com propagação de ondas de choque. Para isso é necessário o estudo de três diferentes aspectos: a dinâmica das estruturas computacional, a dinâmica dos fluidos computacional e o problema do acoplamento. No caso da dinâmica das estruturas computacional deve-se identificar em função da cinemática de deformações, quais são os requisitos para que um elemento seja adequado para analisar tais problemas, tendo em vista que a formulação deve admitir grandes deslocamentos. Para evitar problemas relacionados com aproximações de rotações finitas, opta-se por empregar uma formulação descrita em termos de posições e que leva em consideração os efeitos de empenamento da seção transversal. No caso da dinâmica dos fluidos computacional, busca-se uma formulação para escoamentos compressíveis que seja estável e ao mesmo tempo sensível ao movimento da estrutura, sendo empregado um algoritmo de integração temporal explícito baseado em características com as equações governantes descritas na forma Lagrangeana-Euleriana Arbitrária (ALE). No que se refere ao acoplamento, busca-se modularidade e versatilidade, empregando-se um modelo particionado fraco (explícito) de acoplamento e técnicas de transferência das condições de contorno (Dirichlet-Neummann), sendo estudados os efeitos de utilizar transferência bidirecional ou unidirecional dessas condições de contorno. / This work consists in the use of the Finite Element Method (FEM) for numerical analysis of fluid-bar structures, focusing on transient problems involving explosions or other actions with shock waves propagation. For this purpose, one needs to study three different aspects: the computational structural dynamics, the computational fluid dynamics and the coupling problem. Regarding computational structural dynamics, one need firstly to identify the requirements for an element to be adequate to analyze such problems, taking into account the fact that such element should admit large displacements. In order to avoid problems related to finite rotation approximations and to give a realist representation of a 3D bar structure, we chose a formulation defined in terms of positions and that considers the cross-section warping effects. Regarding computational fluid dynamics, we seek for a stable formulation for compressible flows, and at same time, sensitive to the movement of the structure, leading to an explicit time integration algorithm based on characteristics with governing equations described in the Arbitrary Lagrangian-Eulerian (ALE) form. Regarding to coupling, we chose to use a weak (explicit) partitioning coupling model in order to ensure modularity and versatility. The developed coupling scheme is bases on boundary conditions transfer techniques (Dirichlet-Neummann), and we study the effects of using bidirectional or unidirectional boundary conditions transfers.

O presente trabalho consiste no uso do Método dos Elementos Finitos (MEF) para a análise de interação fluido-estruturas de barras com foco em problemas transientes envolvendo explosões ou outras ações com propagação de ondas de choque. Para isso é necessário o estudo de três diferentes aspectos: a dinâmica das estruturas computacional, a dinâmica dos fluidos computacional e o problema do acoplamento. No caso da dinâmica das estruturas computacional deve-se identificar em função da cinemática de deformações, quais são os requisitos para que um elemento seja adequado para analisar tais problemas, tendo em vista que a formulação deve admitir grandes deslocamentos. Para evitar problemas relacionados com aproximações de rotações finitas, opta-se por empregar uma formulação descrita em termos de posições e que leva em consideração os efeitos de empenamento da seção transversal. No caso da dinâmica dos fluidos computacional, busca-se uma formulação para escoamentos compressíveis que seja estável e ao mesmo tempo sensível ao movimento da estrutura, sendo empregado um algoritmo de integração temporal explícito baseado em características com as equações governantes descritas na forma Lagrangeana-Euleriana Arbitrária (ALE). No que se refere ao acoplamento, busca-se modularidade e versatilidade, empregando-se um modelo particionado fraco (explícito) de acoplamento e técnicas de transferência das condições de contorno (Dirichlet-Neummann), sendo estudados os efeitos de utilizar transferência bidirecional ou unidirecional dessas condições de contorno. / This work consists in the use of the Finite Element Method (FEM) for numerical analysis of fluid-bar structures, focusing on transient problems involving explosions or other actions with shock waves propagation. For this purpose, one needs to study three different aspects: the computational structural dynamics, the computational fluid dynamics and the coupling problem. Regarding computational structural dynamics, one need firstly to identify the requirements for an element to be adequate to analyze such problems, taking into account the fact that such element should admit large displacements. In order to avoid problems related to finite rotation approximations and to give a realist representation of a 3D bar structure, we chose a formulation defined in terms of positions and that considers the cross-section warping effects. Regarding computational fluid dynamics, we seek for a stable formulation for compressible flows, and at same time, sensitive to the movement of the structure, leading to an explicit time integration algorithm based on characteristics with governing equations described in the Arbitrary Lagrangian-Eulerian (ALE) form. Regarding to coupling, we chose to use a weak (explicit) partitioning coupling model in order to ensure modularity and versatility. The developed coupling scheme is bases on boundary conditions transfer techniques (Dirichlet-Neummann), and we study the effects of using bidirectional or unidirectional boundary conditions transfers.