Made available in DSpace on 2018-08-01T23:58:51Z (GMT). No. of bitstreams: 1 tese_11417_Vitor Folador.PDF: 45 bytes, checksum: e7c5f9cd9e7e37604f4b7be2b3ae7a82 (MD5) Previous issue date: 2017-06-28 / O Método das Diferenças Finitas (MDF) é um método numérico para soluções de equações diferenciais em que as derivadas são aproximadas por uma equação de diferença finita. O MDF é, portanto, um método baseado no processo de discretização. Entretanto, sua abordagem é bem mais simples e direta se comparado a outros famosos métodos numéricos tais como o Método dos Elementos Finitos, Método dos Elementos de Contorno e o Método dos Volumes Finitos. A abordagem do MDF relacionado à engenharia é geralmente tratada superficialmente em livros e artigos científicos, resumindo-se a capítulos de introdução. Daí surgiu o interesse em aprofundar-se em diferentes estratégias de solução de diversos problemas de engenharia via MDF. Este trabalho busca resgatar o uso do Método das Diferenças Finitas na engenharia computacional com base no recente crescimento em trabalhos científicos publicados sobre este método. Esta pesquisa trata da aplicação do MDF em diversos problemas de engenharia, que por sua vez é uma área de estudo dentro do campo de métodos numéricos. A pesquisa está diretamente relacionada ao estudo do MDF como ferramenta para encontrar soluções aproximadas de equações diferenciais que regem diversos problemas de engenharia. Problemas como cálculo de deflexão em vigas, flambagem de pilares, flexão de placas, problemas de calor em regime permanente e transiente (incluindo vigas de concreto em situação de incêndio) e vigas sobre base elástica são exemplos de problemas muito comuns na engenharia estrutural e farão parte do escopo desta pesquisa. Para o entendimento do assunto, é apresentada uma síntese da fundamentação teórica de cada problema de engenharia de forma a contextualizar cada situação. Com as equações diferenciais governantes apresentadas, procede-se ao detalhamento das principais estratégias de solução via MDF, incluindo exemplos numéricos com uma variedade de condições de contorno. E, ao final deste trabalho, a comunidade científica poderá ter acesso a um material todo direcionado a este assunto, auxiliando alunos e pesquisadores da área.

Neste projeto, é apresentado um método numérico com uma abordagem do tipoMAC para a simulação de escoamentos viscoelásticos incompressíveis tridimensionais com superfície livre governados pelo modelo de fluido SXPP. A formulação apresentada nesse trabalho é uma extensão dos resultados obtidos por Oishi et al. (2011), sobre o estudo de métodos numéricos para a simulação de escoamentos incompressíveis viscoelásticos com superfície livre a baixos números de Reynolds, para o caso bidimensional. No contexto de problemas transientes, metodologias explícitas para solução numérica das equações governantes apresentam restrições de estabilidade muito severas para a definição do passo temporal, acarretando em um custo computacional relativamente alto. Sendo assim, utilizamos um método implícito para resolver a equação de conservação da quantidade de movimento, eliminando assim, a restrição de estabilidade parabólica e diminuindo significativamente o custo computacional. Mas tal estratégia acopla os campos de velocidade e pressão. Dessa forma, para desacoplar esses campos, foi utilizado uma abordagem que combina método de projeção com uma técnica implícita para o tratamento da pressão na superfície livre. A equação constitutiva foi resolvida pelo método de Runge-Kutta de segunda-ordem. A validação do método numérico foi realizada utilizando refinamento da malha no escoamento em um canal. Como aplicação, apresentamos resultados numéricos sobre o problema do jato oscilante e do inchamento do extrudado / In this work, we present a numerical method with a MAC type approach to simulate tridimensional incompressible viscoelastic free surface flows governed by a SXPP (Single eXtended Pom-Pom) model. The formulation presented in this work is an extension to the work of Oishi et al. (2011). They have studied numerical methods for solving incompressible viscoelastic free surface flows with low Reynolds number, for the bidimensional case. In the context of transient problems, explicitmethodologies for the numerical solution of the governing equations present severe stability constraints for defining the time step, what highly increases the computational cost. Due to this fact, an implicit method is used to solve the momentum equation, eliminating the parabolic stability constraint and decreasing significantly the computational cost. However, this strategy couples velocity and pressure fields. To decouples this fields, it was used an approach that combines a projection method and an implicit technique for the treatment of the pressure at the free surface. The constitutive equation is solved by a second-order Runge-Kutta method. The numerical method validation was achieved by a mesh refinement for a flow in a channel. As applications, numerical results of the die-swell problem and the jet buckling phenomenon are presented

Tese de mestrado integrado. Mestrado Integrado em Engenharia Civil. Especialização em Hidráulica. Faculdade de Engenharia. Universidade do Porto. 2012

Tese de mestrado. Engenharia Mecânica (perfil de Construções Mecânicas). Faculdade de Engenharia. Universidade do Porto. 2008

Programa de Pós Graduação em Engenharia Civil. Departamento de Engenharia Civil, Escola de Minas, Universidade Federal de Ouro Preto. / Submitted by Oliveira Flávia (flavia@sisbin.ufop.br) on 2014-12-04T16:41:22Z No. of bitstreams: 1 TESE_AnáliseComputacionalMDF.pdf: 3056157 bytes, checksum: 9ab42896c7f77eab74d0c9eb430f8f96 (MD5) / Approved for entry into archive by Gracilene Carvalho (gracilene@sisbin.ufop.br) on 2014-12-09T20:23:05Z (GMT) No. of bitstreams: 1 TESE_AnáliseComputacionalMDF.pdf: 3056157 bytes, checksum: 9ab42896c7f77eab74d0c9eb430f8f96 (MD5) / Submitted by Oliveira Flávia (flavia@sisbin.ufop.br) on 2015-01-23T17:19:21Z No. of bitstreams: 1 TESE_AnáliseComputacionalMDF.pdf: 3056157 bytes, checksum: 9ab42896c7f77eab74d0c9eb430f8f96 (MD5) / Approved for entry into archive by Gracilene Carvalho (gracilene@sisbin.ufop.br) on 2015-01-26T11:35:33Z (GMT) No. of bitstreams: 1 TESE_AnáliseComputacionalMDF.pdf: 3056157 bytes, checksum: 9ab42896c7f77eab74d0c9eb430f8f96 (MD5) / Made available in DSpace on 2015-01-26T11:35:33Z (GMT). No. of bitstreams: 1 TESE_AnáliseComputacionalMDF.pdf: 3056157 bytes, checksum: 9ab42896c7f77eab74d0c9eb430f8f96 (MD5) Previous issue date: 2013 / As cascas cilíndricas são elementos estruturais bastante utilizados como tubulações enterradas, estruturas de contenção de túneis e poços, e depósitos de diversos materiais, na forma de silos ou tanques. Nessas aplicações as cascas cilíndricas interagem com o meio elástico, no caso solo ou rocha. Essa interação pode ser considerada de duas forma diferentes: com um problema de contato bilateral (PCB), quando a estrutura e o meio não perdem contato durante o processo de deformação; e como um problema de contato unilateral (PCU), quando a casca não está completamente fixa e pode perder contato com o meio para determinadas condições de carregamento. Portanto, esta tese tem como principal objetivo o desenvolvimento de uma ferramenta, ou programa computacional, para o estudo e análise de problemas envolvendo cascas cilíndricas com restrições bilaterais e unilaterais de contato impostas por bases elásticas. Essa ferramenta computacional apresenta como principais características: i. o emprego do Método das Diferenças Finitas (MDF) para aproximar as derivadas das equações diferenciais de equilíbrio do sistema casca-base elástica. Adota-se neste trabalho a Teoria de Casca de Sanders, que pode ser aplicada para cascas ou painéis cilíndricos pouco ou muito abatidos; e as bases elásticas são aproximadas pelos modelos de Winkler e Pasternak; ii. o emprego de metodologias de solução desses PCB e PCU. No caso do PCB, a região de contato Sc entre os corpos (casca-base) é conhecida a priori, e sua solução pode ser obtida de forma direta, se as hipóteses de pequenos deslocamentos e material elástico forem consideradas. No caso do PCU, tem-se como principal desafio a definição da região de contato entre os corpos, Sc, que não é conhecida a priori, e adotam-se aqui duas estratégias de solução desse problema (ES1 e ES2, respectivamente). Na primeira, a interação casca-base é avaliada considerando diretamente as restrições de contato unilateral e resolvendo-se um Problema de Complementaridade Linear (PCL). Na segunda, a região de contato Sc é inicialmente aproximada e na sequência utiliza-se o método de Newton-Raphson para corrigi-la e avaliar a participação da base elástica na obtenção de outras incógnitas do problema. Essa ferramenta computacional foi desenvolvida usando a linguagem de programação Fortran®. Diversos exemplos numéricos são apresentados e analisados com o intuito de averiguar a eficiência das metodologias implementadas para soluções dos problemas de contato. No final da tese, algumas conclusões e observações referentes ao sistema computacional desenvolvido e às metodologias numéricas propostas são estabelecidas. Além de o MDF ter se apresentado como um método de fácil implementação e utilização, as estratégias de solução não linear adotadas se mostraram semelhantes, em termos de resultados obtidos, mas com ligeiro ganho de eficiência computacional por parte da ES1. Por fim, sugestões para novas pesquisas são apresentadas. ____________________________________________________________________________ / ABSTRACT: Buried pipelines, containment structures for tunnels and shafts, and deposits of various materials in the form of silos or tanks all use structural elements known as cylindrical shells. Cylindrical shells interact in these applications with elastic means, in these cases soil or rock. This interaction can be considered in two ways. First, it can be considered a bilateral contact problem (BCP) where, during the deformation process, structure and base remain in contact; second, it can be considered a unilateral contact problem (UCP) where, under certain loading conditions, the shell is not completely fixed and may lose contact with the medium. The main objective of this thesis then is to develop a computational tool for the study and analysis of problems, imposed by elastic base contact, involving cylindrical shells with bilateral and unilateral restrictions. This numerical tool has as main features: i. the use of the finite difference method (FDM) to approximate the derivatives of the differential equations of equilibrium of the system’s shell-elastic base. This paper adopts Sanders’ theory for shells, which can be applied to cylindrical shells and panels, and approximates the elastic bases using Winkler’s and Pasternak’s models. ii. Employment of methodologies for solving these BCPs and UCPs. In the case of BCP, the region of contact between the bodies Sc (shell-base) is known a priori; its solution can be obtained directly if one considers the assumptions of small displacements and elastic material. In the case of UCP, the main challenge is the definition of the area of contact between the bodies Sc, which is unknown a priori; to solve this problem, the current study adopts two strategies for solving this problem (ES1 and ES2, respectively). First, the interaction shell-base is calculated by considering directly the unilateral contact constraints and solving a linear complementarity problem (LCP). Second, the contact area Sc is initially approximated and corrected, via Newton-Raphson’s method, by evaluating the participation of the elastic foundation in obtaining other unknown variables. This numerical tool was developed using the programming language Fortran ®. To determine the efficiency of these methodologies, the work presents and analyzes several numerical examples. The end of the thesis, offers some conclusions and remarks regarding the computing system developed and the proposed numerical methodologies. Besides the MDF have been presented as a method easy to implement and use, the nonlinear solution strategies adopted were similar in terms of results, though ES1 showed a slight gain in computational efficiency. Finally, suggestions for further research are presented.

Neste projeto, é apresentado um método numérico com uma abordagem do tipoMAC para a simulação de escoamentos viscoelásticos incompressíveis tridimensionais com superfície livre governados pelo modelo de fluido SXPP. A formulação apresentada nesse trabalho é uma extensão dos resultados obtidos por Oishi et al. (2011), sobre o estudo de métodos numéricos para a simulação de escoamentos incompressíveis viscoelásticos com superfície livre a baixos números de Reynolds, para o caso bidimensional. No contexto de problemas transientes, metodologias explícitas para solução numérica das equações governantes apresentam restrições de estabilidade muito severas para a definição do passo temporal, acarretando em um custo computacional relativamente alto. Sendo assim, utilizamos um método implícito para resolver a equação de conservação da quantidade de movimento, eliminando assim, a restrição de estabilidade parabólica e diminuindo significativamente o custo computacional. Mas tal estratégia acopla os campos de velocidade e pressão. Dessa forma, para desacoplar esses campos, foi utilizado uma abordagem que combina método de projeção com uma técnica implícita para o tratamento da pressão na superfície livre. A equação constitutiva foi resolvida pelo método de Runge-Kutta de segunda-ordem. A validação do método numérico foi realizada utilizando refinamento da malha no escoamento em um canal. Como aplicação, apresentamos resultados numéricos sobre o problema do jato oscilante e do inchamento do extrudado / In this work, we present a numerical method with a MAC type approach to simulate tridimensional incompressible viscoelastic free surface flows governed by a SXPP (Single eXtended Pom-Pom) model. The formulation presented in this work is an extension to the work of Oishi et al. (2011). They have studied numerical methods for solving incompressible viscoelastic free surface flows with low Reynolds number, for the bidimensional case. In the context of transient problems, explicitmethodologies for the numerical solution of the governing equations present severe stability constraints for defining the time step, what highly increases the computational cost. Due to this fact, an implicit method is used to solve the momentum equation, eliminating the parabolic stability constraint and decreasing significantly the computational cost. However, this strategy couples velocity and pressure fields. To decouples this fields, it was used an approach that combines a projection method and an implicit technique for the treatment of the pressure at the free surface. The constitutive equation is solved by a second-order Runge-Kutta method. The numerical method validation was achieved by a mesh refinement for a flow in a channel. As applications, numerical results of the die-swell problem and the jet buckling phenomenon are presented

Este trabalho apresenta parte de um sistema de simulação integrado para escoamento de fluido incompressível bidimensional em malhas não estruturadas denominado UmFlow-2D. O sistema consiste de três módulos: um módulo modelador, um módulo simulador e um módulo visualizador. A parte do sistema apresentado neste trabalho é o módulo simulador. Este módulo, implementa as equações de Navier-Stokes. As equações governantes são discretizadas pelo método de diferenças finitas generalizadas e os termos convectivos pelo método semi-lagrangeano. Um método de projeção é empregado para desacoplar as componentes da velocidade e pressão. O gerenciamento da malha, não estruturada é feito pela estrutura de dados SHE. Os resultados numéricos obtidos pelo UmFlow-2D são comparados com soluções analíticas e soluções numéricas de outros trabalhos. / This work presents an integratc simulation system, called UmFlow-2D, wich aims a,t simulating two-dimensional íncompressible fluid flow using unstructed mesh. The system is divided three modules: modeling module, simulation module and visualization module. In this work we present the simulation module. The simulation module implements the Navier-Stokes equation. The governing equations are discretized by a generalized flnite dillerence method and the convective terms by semi-lagrangean method. A projection method is employed to uncouple the velocity componentes and pressure. The management at the unstructed mesh is ready using a data structure called SHE. The numérica! results are compared with analytical solutions and numerical simulations of other works.

Este trabalho propõe um esquema de diferenças finitas compacta de alta ordem para resolver problemas elípticos com coeficientes variáveis em malhas composta. São apresentados a formulação matemática e a dedução do método compacto de quarta ordem aplicado à problemas elípticos bidimensionais, em malha regular e composta. Foi adotado o uso da biblioteca PETSc com os seus pré-condicionadores e métodos numéricos para resolver os sistemas lineares resultantes da discretização do problema. Por fim, testes visando verificar o código foram feitos, utilizando o método de soluções manufaturadas, para mostrar alta eficiência e acurácia do método desenvolvido. / This paper proposes a scheme of compact finite difference higher order for solve elliptic problems with variable coeficients in composite meshes. we present the mathematical formulation and the deduction of the compact method of fourth order applied to two-dimensional elliptic problems in regular and composite mesh . It was adopted using the PETSc library with its pre- conditioners and numerical methods for solving linear systems resulting from discretization of the problem. Finally , tests to verify the code were made using the method of manufactured solutions to show high eficiency and accuracy of the method developed .

Submitted by Fernanda Rodrigues de Lima (fernanda.rlima@ufpe.br) on 2018-07-31T19:30:12Z No. of bitstreams: 2 license_rdf: 811 bytes, checksum: e39d27027a6cc9cb039ad269a5db8e34 (MD5) DISSERTAÇÃO Priscila Feitosa de Sá Ferreira.pdf: 7224560 bytes, checksum: c24277ade9dbd50cb5253b8573959668 (MD5) / Approved for entry into archive by Alice Araujo (alice.caraujo@ufpe.br) on 2018-07-31T20:47:10Z (GMT) No. of bitstreams: 2 license_rdf: 811 bytes, checksum: e39d27027a6cc9cb039ad269a5db8e34 (MD5) DISSERTAÇÃO Priscila Feitosa de Sá Ferreira.pdf: 7224560 bytes, checksum: c24277ade9dbd50cb5253b8573959668 (MD5) / Made available in DSpace on 2018-07-31T20:47:10Z (GMT). No. of bitstreams: 2 license_rdf: 811 bytes, checksum: e39d27027a6cc9cb039ad269a5db8e34 (MD5) DISSERTAÇÃO Priscila Feitosa de Sá Ferreira.pdf: 7224560 bytes, checksum: c24277ade9dbd50cb5253b8573959668 (MD5) Previous issue date: 2017-04-24 / O presente trabalho tem como objetivo propor quatro modelos de análise diferentes para o dimensionamento de estacas semienterradas de seção circular submetidas a cargas laterais com a consideração da não-linearidade física do material (NLF) e da não-linearidade geométrica da estrutura (NLG). Estes modelos, chamados de Níveis de Aproximação (NdA), estão organizados em ordem crescente de precisão da análise, partindo de um nível inicial, com maiores simplificações, até um nível final, com uma modelagem mais realista. No NdA I, a estaca semienterrada é tratada como uma coluna engastada na base e livre na ponta, com base na Técnica do Engastamento Virtual. O cálculo dos esforços resistentes é feito com base no método do pilar-padrão com curvatura aproximada, proposto pela norma brasileira. No NdA II, um programa comercial de análise tridimensional de estruturas é utilizado. Neste programa, o solo é representado de maneira discreta através de uma série de molas de comportamento elástico e linear (baseadas no Modelo de Winkler). A NLF é tratada de maneira aproximada através de um coeficiente, aplicado ao momento de inércia bruto, que simula a diminuição da rigidez do concreto. A NLG é tratada internamente pelo programa através do processo iterativo P-delta. No NdA III, é utilizado um programa elaborado neste trabalho que considera as não-linearidades e a interação solo-estrutura através da equação da viga sobre base elástica, cuja solução é dada a partir do Método das Diferenças Finitas. A NLF, neste modelo, foi considerada através de diagramas momento-curvatura. A NLG e a interação solo-estrutura (Modelo de Winkler) são consideradas diretamente na formulação da equação diferencial. Por fim, o NdA IV utiliza um programa comercial de análise não-linear de estruturas de pontes através de elementos finitos. A NLF é tratada pela integração das curvas de tensãodeformação e a NLG é modelada usando momentos P-delta. A consideração da interação solo-estrutura, neste modelo, é feita através de modelos de curvas p-y. Estes quatro modelos foram aplicados a dois exemplos distintos: um considerando apenas um pilar engastado na base e livre no topo (livre da interação solo-estrutura) e outro considerando uma estaca semienterrada em solo arenoso. Os resultados obtidos pelos quatro métodos foram bastante satisfatórios e se mostraram coerentes entre si, apresentando-se dentro do esperado para cada grau de refinamento da análise. / The present work proposes four different analysis models for the design of partly buried piles of circular section subjected to lateral loads with consideration of the physical nonlinearity of the material (PNL) and the geometric nonlinearity of the structure (GNL). These models, called Levels of Approximation (LoA), are organized in increasing order of accuracy, starting from an initial level, with greater simplifications, to a final level, which is more realistic. In LoA I, the partly buried pile is treated as a column fixed at the base and free at the top, based on the Virtual Clamping Technique. Strength computation is based on the standard column method with approximate curvature, proposed by the Brazilian Code. In LoA II, a commercial three-dimensional structure analysis program is used. In this program, the soil is represented discreetly through a series of elastic and linear springs (based on the Winkler Model). The PNL is approximately treated by a factor, applied to the gross moment of inertia, which simulates the decrease of rigidity of the reinforced concrete due to cracking and nonlinear constitutive properties. The GNL is internally handled by the program through the iterative P-delta process. The third proposed LoA is a program developed in this work that considers the non-linearities and the soil-structure interaction through the equation of the beam on elastic foundation. The equation is solved numerically, with the application of the Finite Differences Method. The PNL, in this model, was considered through moment-curvature diagrams computed internally for the section under analysis, which were used to determine the stiffness of the pile as a function of the lateral displacement. The GNL and the soil-structure interaction (also by Winkler Model) are considered directly in the equation. Finally, LoA IV uses a commercial program of nonlinear analysis of bridge structures through finite elements. The PNL is treated by integrating stress-strain curves into the cross-section analyzed and the GNL is modeled using P-delta moments. Consideration of the soil-structure interaction in this model is made through p-y curves models. These four models were applied to two distinct examples: one considering only one column fixed at the base and free at the top (free of soil-structure interaction) and another one considering a partly buried pile in sand. The results obtained by the four methods were very satisfactory and were shown to be coherent with each other, in line with the expected for each degree of refinement of the analysis.

Este trabalho propõe um esquema de diferenças finitas compacta de alta ordem para resolver problemas elípticos com coeficientes variáveis em malhas composta. São apresentados a formulação matemática e a dedução do método compacto de quarta ordem aplicado à problemas elípticos bidimensionais, em malha regular e composta. Foi adotado o uso da biblioteca PETSc com os seus pré-condicionadores e métodos numéricos para resolver os sistemas lineares resultantes da discretização do problema. Por fim, testes visando verificar o código foram feitos, utilizando o método de soluções manufaturadas, para mostrar alta eficiência e acurácia do método desenvolvido. / This paper proposes a scheme of compact finite difference higher order for solve elliptic problems with variable coeficients in composite meshes. we present the mathematical formulation and the deduction of the compact method of fourth order applied to two-dimensional elliptic problems in regular and composite mesh . It was adopted using the PETSc library with its pre- conditioners and numerical methods for solving linear systems resulting from discretization of the problem. Finally , tests to verify the code were made using the method of manufactured solutions to show high eficiency and accuracy of the method developed .