Blade design and forming for fans using finite elements

The necessity of adapting the standardized fan models to conditions of higher temperature has emerged due to the growth of concerning referring to the consequences of the gas expelling after the Mont Blanc tunnel accident in Italy and France, where even though, with 100 fans in operation, 41 people died. However, since then, the defied solutions have pointed to aerodynamic disadvantages or have seemed nonappropriate in these conditions. The objective of this work is to present an alternative to the market standard fans considering a new technology in constructing blades. This new technology introduces the use of the stainless steel AISI 409 due to its good adaptation to temperatures higher than 400°C, particularly exposed to temperatures of gas exhaust from tunnels in fire situation. Furthermore, it presents a very good resistance to corrosion and posterior welding and pressing, due to its alloyed elements. The innovation is centered in the process of a deep drawing of metallic shells and posterior welding, in order to keep the ideal aerodynamic superficies for the fan ideal performance. On the other hand, the finite element method, through the elasto-plastic software COSMOS permitted the verification of the thickness and structural stability of the blade in relation to the aerodynamic efforts established in the project. In addition, it is not advisable the fabrication of blades with variable localized thickness not even, non-uniform ones, due to the verified concentration of tensions and the difficulties observed in the forming. In this way, this study recommends the construction of blades with uniform variations of thickness. © 2007 Springer.

Análise da execução de um algoritmo numérico cbs para simulação de escoamento utilizando elementos finitos

This work intends to analyze the application and execution time of a numerical algorithm that simulates incompressible and isothermal flows. It was used the explicit scheme of the Characteristic Based Split (CBS) algorithm and the Artificial Compressibility (AC) scheme for coupling pressure-velocity equations. The discretization was done with the finite element method using a bilinear elements grid. The free software GNU-Octave was used for implementation and execution of routines. The results were analyzed using the classic lid-driven cavity problem. This work shows results for tests with several Reynolds' number. The results for these tests show a good agreement when compared with previous ones obtained from bibliography. The code runtime's analysis shows yet that the matrix's assembly is the part of greater consumption time in the implementation.

Integrating open source codes for solid modeling and automatic mesh generation in FEM applications

The applications of the Finite Element Method (FEM) for three-dimensional domains are already well documented in the framework of Computational Electromagnetics. However, despite the power and reliability of this technique for solving partial differential equations, there are only a few examples of open source codes available and dedicated to the solid modeling and automatic constrained tetrahedralization, which are the most time consuming steps in a typical three-dimensional FEM simulation. Besides, these open source codes are usually developed separately by distinct software teams, and even under conflicting specifications. In this paper, we describe an experiment of open source code integration for solid modeling and automatic mesh generation. The integration strategy and techniques are discussed, and examples and performance results are given, specially for complicated and irregular volumes which are not simply connected. © 2011 IEEE.

An anisotropic linear elastic boundary element formulation for masonry wall analysis

This work presents an application of a Boundary Element Method (BEM) formulation for anisotropic body analysis using isotropic fundamental solution. The anisotropy is considered by expressing a residual elastic tensor as the difference of the anisotropic and isotropic elastic tensors. Internal variables and cell discretization of the domain are considered. Masonry is a composite material consisting of bricks (masonry units), mortar and the bond between them and it is necessary to take account of anisotropy in this type of structure. The paper presents the formulation, the elastic tensor of the anisotropic medium properties and the algebraic procedure. Two examples are shown to validate the formulation and good agreement was obtained when comparing analytical and numerical results. Two further examples in which masonry walls were simulated, are used to demonstrate that the presented formulation shows close agreement between BE numerical results and different Finite Element (FE) models. © 2012 Elsevier Ltd.

Análise do comportamento experimental e numérico de prismas de alvenaria estrutural submetidos a ações verticais utilizando elementos finitos volumétricos

Pós-graduação em Engenharia Civil - FEIS

Estudo geofísico regional sobre águas subterrâneas na Ilha do Marajó-Pará-Brasil

Nos últimos dez anos foram realizadas na parte leste da Ilha de Marajó (região dos campos naturais) pelo IDESP e NCGG, mais de 800 SEVs para fins hidrogeológicos. Na época, grande parte dessas SEVs não foram totalmente interpretadas em forma quantitativa, devido à falta de recursos técnicos para fazê-lo de forma eficiente. Agora, usando meios mais modernos para interpretação automática de SEVs, voltou-se a interpretá-las com a finalidade de apresentar uma visão regional dos principais aquíferos da área, agrupar as SEVs em famílias características, testar até que ponto essa interpretação é confiável e propor o modelamento bidimensional como técnica alternativa para interpretar as SEVs realizadas em certos locais da área em questão. Como resultado dessa interpretação, com base na teoria convencional dos meios estratificados, foram definidos três tipos de sistemas de aquíferos. 1. O primeiro, denominado de aquífero profundo, situado a profundidades maiores que 50m, estende-se por toda a região prospectada, estando provavelmente associada às camadas superiores da Formação Marajó ou às litologias altamente resistivas das camadas mais profundas do Grupo Pará. 2. O segundo, denominado de aquífero raso e de média profundidade, localiza-se na parte sul e sudeste da região a profundidades compreendidas entre 10 a 50m, e está associado às lentes arenosas do Grupo Pará. 3. O terceiro, é constituído pelos paleocanais e estruturas similares, distribuídos aleatoriamente na região a pouca profundidade. A partir do estudo detalhado das SEVs, decidiu-se classificá-las em 3 famílias características com seus respectivos tipos e apresentar mapas de localização e da espessura dos aquíferos, bem como mapas de condutância longitudinal total e resistividade média da área. Estes últimos, permitem que se divida a região dos campos da Ilha de Marajó em três zonas principais: 1. Uma, altamente resistiva, situada ao sul e sudeste, a qual coincide com os terrenos aflorantes do Grupo Pará. 2. Outra, altamente condutiva, está localizada no centro e norte, onde se encontram aleatoriamente distribuídos os paleocanais e coincide com os terrenos topograficamente mais baixos, geralmente argilosos e embebidos de água salgada, que são procedentes da erosão dos terrenos circundantes topograficamente mais altos. 3. A última é medianamente resistiva e está relacionada com os terrenos vizinhos à cidade de Chaves (noroeste da região dos campos), os quais apresentam semelhanças com os do sul e sudeste da área. Usando-se a técnica de inversão na interpretação de uma SEV característica de cada família, testou-se, através do seu tratamento estatístico, até que ponto os modelos usados na interpretação dessas SEVs (teoria convencional dos meios estratificados) seriam confiáveis. Conclui-se, então, que a alta correlação existente entre os parâmetros dos modelos assumidos (camadas horizontais, isotrópicas e homogêneas) pode-se dever à utilização de modelos geofísicos muito simples para interpretar a complexa geologia de Marajó. Tendo-se verificado que nem sempre é possível aplicar a teoria das SEVs em meios horizontalmente estratificados para interpretar SEVs obtidas em certos locais de Marajó, os quais muitas vezes apresentam bruscas variações laterais de resistividade, passou-se a demonstrar que estas variações laterais afetam profundamente os dados das SEVs, utilizando-se para isto a técnica dos elementos finitos, a qual leva em conta essa variação bidimensional das propriedades físicas do meio. Foi também possível com esta técnica, modelar uma estrutura rasa, semelhante a um paleocanal, concluindo-se que estes resultados sugerem o emprego, duma forma mais profunda, deste tipo de tratamento para os dados obtidos na região dos campos da Ilha de Marajó.

Modelagem direta de dados de eletrorresistividade 3-D

Dentre os métodos geofísicos, o da Eletrorresistividade é um dos mais tradicionais, com o seu desenvolvimento ocorrido há mais de 80 anos. Durante esse tempo o seu uso acompanhou o avanço no poder de processamento numérico e mais recentemente, a modelagem e inversão tridimensional tornou-se uma possibilidade para o geofísico. Apresentamos, neste trabalho, a técnica de elementos finitos aplicada ao método da eletrorresistividade 3-D, através do cálculo do potencial secundário. Para o desenvolvimento da metodologia, simulamos o levantamento do método da eletrorresistividade 3-D com os arranjos Dipolo-Dipolo e Schlumberger, visando medir as variações laterais e verticais da resistividade aparente do solo. Estes arranjos consistem na injeção de corrente elétrica na superfície e de medidas de diferenças de potencial elétrico, resultante da interação da corrente elétrica com o solo. Sendo que, as fontes e receptores são localizados de acordo com os arranjos escolhidos para o levantamento. Neste trabalho, as curvas de sondagem e as pseudo-secções de resistividade aparente, são obtidas através da modelagem de eletrorresistividade 3-D, usando malha de elementos finitos regular. Para efeito de validação, os resultados são comparados com a resposta 3-D obtida a partir dos potenciais totais.

Analysis of the Stresses in Corrugated Sheets under Bending

This paper presents three different numerical models for the evaluation of the stresses in corrugated sheets under bending. Regarding the numerical simulations different approaches can be considered, i.e., a elastic linear analysis or a physical nonlinear analysis, that considers criteria to fail for the sheet material. Moreover, the construction of the finite element mesh can be used shell elements or solid elements. The choice of each finite element must be made from the consideration of their representativity before behavior to be simulated. Thus, the numerical modelling in this manuscript was performed from the three-dimensional models using the SAP2000Nonlinear software, version 7.42, which has as base the finite elements method (FEM). It was considered shell elements in the build the mesh of finite elements and an analysis of type elastic linear in this case. Five mm thick sheets were evaluated considering three different longitudinal dimensions (spans), i.e., 1100 mm, 1530 mm and 1830 mm. The applied load to the models was 2500 N/m and it was verified that the spans of support of sheets have a significant influence on the results of stresses. The sheets with larger spans present larger stresses for the same applied load. The most intense values of tension occur in the troughs (low waves) of the sheets, on the lower surface, while the most intense values of compression occur in the crests (high waves), on the upper surface of the sheet. The flanks, which are the parts among the troughs and crests of the sheets, are submitted to low levels of stresses. The numeric results of the stresses showed a good agreement with the results obtained from other researchers(3) and these results can be used to predict the behavior of corrugated sheets under bending.

A method for calculating the compliance of bonded-interfaces under shrinkage: Validation for Class I cavities

Fundação de Amparo à Pesquisa do Estado de São Paulo (FAPESP)

Nonlinear formulation based on FEM, Mazars damage criterion and Fick's law applied to failure assessment of reinforced concrete structures subjected to chloride ingress and reinforcements corrosion

Structural durability is an important design criterion, which must be assessed for every type of structure. In this regard, especial attention must be addressed to the durability of reinforced concrete (RC) structures. When RC structures are located in aggressive environments, its durability is strongly reduced by physical/chemical/mechanical processes that trigger the corrosion of reinforcements. Among these processes, the diffusion of chlorides is recognized as one of major responsible of corrosion phenomenon start. To accurate modelling the corrosion of reinforcements and to assess the durability of RC structures, a mechanical model that accounts realistically for both concrete and steel mechanical behaviour must be considered. In this context, this study presents a numerical nonlinear formulation based on the finite element method applied to structural analysis of RC structures subjected to chloride penetration and reinforcements corrosion. The physical nonlinearity of concrete is described by Mazars damage model whereas for reinforcements elastoplastic criteria are adopted. The steel loss along time due to corrosion is modelled using an empirical approach presented in literature and the chloride concentration growth along structural cover is represented by Fick's law. The proposed model is applied to analysis of bended structures. The results obtained by the proposed numerical approach are compared to responses available in literature in order to illustrate the evolution of structural resistant load after corrosion start. (C) 2014 Elsevier Ltd. All rights reserved.

Modelagem em 3D de uma patela humana e análise de esforços utilizando o método dos elementos finitos

Throughout the history of medicine surgeons realized the importance of the patella to the functioning of the knee. The main function of the patella is to increase the mechanical efficiency of the quadriceps tendon and knee extensor mechanism. It was found that 50% to 80% of the fractures without deviation of the patella have the transversal pattern, possibly due to excessive tensile forces applied to the extensor mechanism. The purpose of this study is to analyze the loads to which a patella is submitted during a normal extension movement of knee. This analysis will be done by modeling a 3D patella and subsequent load simulation as, described in medical literature, using the finite element method

Análise da influência de variantes do método de solução de escoamentos utilizando os métodos de elementos finitos com compressibilidade artificial e pseudo-características

Coordenação de Aperfeiçoamento de Pessoal de Nível Superior (CAPES)

A method to investigate the shrinkage stress developed by resin-composites bonded to a single flat surface

Objectives. To purpose a method for predicting the shrinkage stress development in the adhesive layer of resin-composite cylinders that shrink bonded to a single flat surface, by measuring the deflection of a glass coverslip caused by the shrinkage of the bonded cylinders. The correlation between the volume of the bonded resin-composite and the stress-peak was also investigated. Methods. A glass coverslip deflection caused by the shrinkage of a bonded resin-composite cylinder (diameter: d = 8 mm, 4 mm, or 2 mm, height: h = 4 mm, 2 mm, 1 mm, or 0.5 mm) was measured, and the same set-up was simulated by finite element analysis (3D-FEA). Stresses generated in the adhesive layer were plotted versus two geometric variables of the resin-composite cylinder (C-Factor and volume) to verify the existence of correlations between them and stresses. Results. The FEA models were validated. A significant correlation (p < 0.01, Pearson's test) between the stress-peak and the coverslip deflection when the resin-composites were grouped by diameter was found for diameters of 2 and 4 mm. The stress-peak of the whole set of data showed a logarithmic correlation with the bonded resin-composite volume (p < 0.001, Pearson's test), but did not correlate with the C-Factor. Significance. The described method should be considered for standardizing the stress generated by the shrinkage of resin-composite blocks bonded to a single flat surface. (C) 2012 Academy of Dental Materials. Published by Elsevier Ltd. All rights reserved.

Safety factors in sheet pile wall design: considerations by using traditional methods and FEM analysis

The purpose of the work is: define and calculate a factor of collapse related to traditional method to design sheet pile walls. Furthermore, we tried to find the parameters that most influence a finite element model representative of this problem. The text is structured in this way: from chapter 1 to 5, we analyzed a series of arguments which are usefull to understanding the problem, while the considerations mainly related to the purpose of the text are reported in the chapters from 6 to 10. In the first part of the document the following arguments are shown: what is a sheet pile wall, what are the codes to be followed for the design of these structures and what they say, how can be formulated a mathematical model of the soil, some fundamentals of finite element analysis, and finally, what are the traditional methods that support the design of sheet pile walls. In the chapter 6 we performed a parametric analysis, giving an answer to the second part of the purpose of the work. Comparing the results from a laboratory test for a cantilever sheet pile wall in a sandy soil, with those provided by a finite element model of the same problem, we concluded that:in modelling a sandy soil we should pay attention to the value of cohesion that we insert in the model (some programs, like Abaqus, don’t accept a null value for this parameter), friction angle and elastic modulus of the soil, they influence significantly the behavior of the system (structure-soil), others parameters, like the dilatancy angle or the Poisson’s ratio, they don’t seem influence it. The logical path that we followed in the second part of the text is reported here. We analyzed two different structures, the first is able to support an excavation of 4 m, while the second an excavation of 7 m. Both structures are first designed by using the traditional method, then these structures are implemented in a finite element program (Abaqus), and they are pushed to collapse by decreasing the friction angle of the soil. The factor of collapse is the ratio between tangents of the initial friction angle and of the friction angle at collapse. At the end, we performed a more detailed analysis of the first structure, observing that, the value of the factor of collapse is influenced by a wide range of parameters including: the value of the coefficients assumed in the traditional method and by the relative stiffness of the structure-soil system. In the majority of cases, we found that the value of the factor of collapse is between and 1.25 and 2. With some considerations, reported in the text, we can compare the values so far found, with the value of the safety factor proposed by the code (linked to the friction angle of the soil).