Metodologia númerico-experimental para avaliação da espessura de filmes de nitreto de titânio obtidos por processamento a plasma

In recent decades, changes in the surface properties of materials have been used to improve their tribological characteristics. However, this improvement depends on the process, treatment time and, primarily, the thickness of this surface film layer. Physical vapor deposition (PVD) of titanium nitrate (TiN) has been used to increase the surface hardness of metallic materials. Thus, the aim of the present study was to propose a numerical-experimental method to assess the film thickness (l) of TiN deposited by PVD. To reach this objective, experimental results of hardness (H) assays were combined with a numerical simulation to study the behavior of this property as a function of maximum penetration depth of the indenter (hmax) into the film/substrate conjugate. Two methodologies were adopted to determine film thickness. The first consists of the numerical results of the H x hmax curve with the experimental curve obtained by the instrumental indentation test. This methodology was used successfully in a TiN-coated titanium (Ti) conjugate. A second strategy combined the numerical results of the Hv x hmax curve with Vickers experimental hardness data (Hv). This methodology was applied to a TiN-coated M2 tool steel conjugate. The mechanical properties of the materials studied were also determined in the present study. The thicknesses results obtained for the two conjugates were compatible with their experimental data.

Modelagem dos efeitos da irradiação luminosa no cérebro de camundongos e rastreamento de neurônios durante experimentos de microscopia de fluorescência

The fluorescent proteins are an essential tool in many fields of biology, since they allow us to watch the development of structures and dynamic processes of cells in living tissue, with the aid of fluorescence microscopy. Optogenectics is another technique that is currently widely used in Neuroscience. In general, this technique allows to activate/deactivate neurons with the radiation of certain wavelengths on the cells that have ion channels sensitive to light, at the same time that can be used with fluorescent proteins. This dissertation has two main objectives. Initially, we study the interaction of light radiation and mice brain tissue to be applied in optogenetic experiments. In this step, we model absorption and scattering effects using mice brain tissue characteristics and Kubelka-Munk theory, for specific wavelengths, as a function of light penetration depth (distance) within the tissue. Furthermore, we model temperature variations using the finite element method to solve Pennes’ bioheat equation, with the aid of COMSOL Multiphysics Modeling Software 4.4, where we simulate protocols of light stimulation tipically used in optogenetics. Subsequently, we develop some computational algorithms to reduce the exposure of neuron cells to the light radiation necessary for the visualization of their emitted fluorescence. At this stage, we describe the image processing techniques developed to be used in fluorescence microscopy to reduce the exposure of the brain samples to continuous light, which is responsible for fluorochrome excitation. The developed techniques are able to track, in real time, a region of interest (ROI) and replace the fluorescence emitted by the cells by a virtual mask, as a result of the overlay of the tracked ROI and the fluorescence information previously stored, preserving cell location, independently of the time exposure to fluorescent light. In summary, this dissertation intends to investigate and describe the effects of light radiation in brain tissue, within the context of Optogenetics, in addition to providing a computational tool to be used in fluorescence microscopy experiments to reduce image bleaching and photodamage due to the intense exposure of fluorescent cells to light radiation.

Desenvolvimento da célula base de microestruturas periódicas de compósitos sob otimização topológica

This thesis develops a new technique for composite microstructures projects by the Topology Optimization process, in order to maximize rigidity, making use of Deformation Energy Method and using a refining scheme h-adaptative to obtain a better defining the topological contours of the microstructure. This is done by distributing materials optimally in a region of pre-established project named as Cell Base. In this paper, the Finite Element Method is used to describe the field and for government equation solution. The mesh is refined iteratively refining so that the Finite Element Mesh is made on all the elements which represent solid materials, and all empty elements containing at least one node in a solid material region. The Finite Element Method chosen for the model is the linear triangular three nodes. As for the resolution of the nonlinear programming problem with constraints we were used Augmented Lagrangian method, and a minimization algorithm based on the direction of the Quasi-Newton type and Armijo-Wolfe conditions assisting in the lowering process. The Cell Base that represents the composite is found from the equivalence between a fictional material and a preescribe material, distributed optimally in the project area. The use of the strain energy method is justified for providing a lower computational cost due to a simpler formulation than traditional homogenization method. The results are presented prescription with change, in displacement with change, in volume restriction and from various initial values of relative densities.

Criação de zonas de alta permeabilidade em fraturas: modelagem numérica do Lajedo de Soledade, Formação Jandaíra, Bacia Potiguar

The fracturing in carbonate rocks has been attracting increasingly attention due to new oil discoveries in carbonate reservoirs. This study investigates how the fractures (faults and joints) behave when subjected to different stress fields and how their behavior may be associated with the generation of karst and consequently to increased secondary porosity in these rocks. In this study I used satellite imagery and unmanned aerial vehicle UAV images and field data to identify and map faults and joints in a carbonate outcrop, which I consider a good analogue of carbonate reservoir. The outcrop comprises rocks of the Jandaíra Formation, Potiguar Basin. Field data were modeled using the TECTOS software, which uses finite element analysis for 2D fracture modeling. I identified three sets of fractures were identified: NS, EW and NW-SE. They correspond to faults that reactivate joint sets. The Ratio of Failure by Stress (RFS) represents stress concentration and how close the rock is to failure and reach the Mohr-Coulomb envelopment. The results indicate that the tectonic stresses are concentrated in preferred structural zones, which are ideal places for carbonate dissolution. Dissolution was observed along sedimentary bedding and fractures throughout the outcrop. However, I observed that the highest values of RFS occur in fracture intersections and terminations. These are site of karst concentration. I finally suggest that there is a relationship between stress concentration and location of karst dissolution in carbonate rocks.

Criação de zonas de alta permeabilidade em fraturas: modelagem numérica do Lajedo de Soledade, Formação Jandaíra, Bacia Potiguar

The fracturing in carbonate rocks has been attracting increasingly attention due to new oil discoveries in carbonate reservoirs. This study investigates how the fractures (faults and joints) behave when subjected to different stress fields and how their behavior may be associated with the generation of karst and consequently to increased secondary porosity in these rocks. In this study I used satellite imagery and unmanned aerial vehicle UAV images and field data to identify and map faults and joints in a carbonate outcrop, which I consider a good analogue of carbonate reservoir. The outcrop comprises rocks of the Jandaíra Formation, Potiguar Basin. Field data were modeled using the TECTOS software, which uses finite element analysis for 2D fracture modeling. I identified three sets of fractures were identified: NS, EW and NW-SE. They correspond to faults that reactivate joint sets. The Ratio of Failure by Stress (RFS) represents stress concentration and how close the rock is to failure and reach the Mohr-Coulomb envelopment. The results indicate that the tectonic stresses are concentrated in preferred structural zones, which are ideal places for carbonate dissolution. Dissolution was observed along sedimentary bedding and fractures throughout the outcrop. However, I observed that the highest values of RFS occur in fracture intersections and terminations. These are site of karst concentration. I finally suggest that there is a relationship between stress concentration and location of karst dissolution in carbonate rocks.

Modelagem multiescala de fenômenos eletrocinéticos em meios porosos carregados eletricamente: aplicação a meios porosos argilosos

In this master thesis, we propose a multiscale mathematical and computational model for electrokinetic phenomena in porous media electrically charged. We consider a porous medium rigid and incompressible saturated by an electrolyte solution containing four monovalent ionic solutes completely diluted in the aqueous solvent. Initially we developed the modeling electrical double layer how objective to compute the electrical potential, surface density of electrical charges and considering two chemical reactions, we propose a 2-pK model for calculating the chemical adsorption occurring in the domain of electrical double layer. Having the nanoscopic model, we deduce a model in the microscale, where the electrochemical adsorption of ions, protonation/ deprotonation reactions and zeta potential obtained in the nanoscale, are incorporated through the conditions of interface uid/solid of the Stokes problem and transportation of ions, modeled by equations of Nernst-Planck. Using the homogenization technique of periodic structures, we develop a model in macroscopic scale with respective cells problems for the e ective macroscopic parameters of equations. Finally, we propose several numerical simulations of the multiscale model for uid ow and transport of reactive ionic solute in a saturated aqueous solution of kaolinite. Using nanoscopic model we propose some numerical simulations of electrochemical adsorption phenomena in the electrical double layer. Making use of the nite element method discretize the macroscopic model and propose some numerical simulations in basic and acid system aiming to quantify the transport of ionic solutes in porous media electrically charged.

Modelagem multiescala de fenômenos eletrocinéticos em meios porosos carregados eletricamente: aplicação a meios porosos argilosos

In this master thesis, we propose a multiscale mathematical and computational model for electrokinetic phenomena in porous media electrically charged. We consider a porous medium rigid and incompressible saturated by an electrolyte solution containing four monovalent ionic solutes completely diluted in the aqueous solvent. Initially we developed the modeling electrical double layer how objective to compute the electrical potential, surface density of electrical charges and considering two chemical reactions, we propose a 2-pK model for calculating the chemical adsorption occurring in the domain of electrical double layer. Having the nanoscopic model, we deduce a model in the microscale, where the electrochemical adsorption of ions, protonation/ deprotonation reactions and zeta potential obtained in the nanoscale, are incorporated through the conditions of interface uid/solid of the Stokes problem and transportation of ions, modeled by equations of Nernst-Planck. Using the homogenization technique of periodic structures, we develop a model in macroscopic scale with respective cells problems for the e ective macroscopic parameters of equations. Finally, we propose several numerical simulations of the multiscale model for uid ow and transport of reactive ionic solute in a saturated aqueous solution of kaolinite. Using nanoscopic model we propose some numerical simulations of electrochemical adsorption phenomena in the electrical double layer. Making use of the nite element method discretize the macroscopic model and propose some numerical simulations in basic and acid system aiming to quantify the transport of ionic solutes in porous media electrically charged.

Contribuição ao estudo de treliças não lineares considerando efeitos térmicos

This work presents the numerical analysis of nonlinear trusses summited to thermomechanical actions with Finite Element Method (FEM). The proposed formulation is so-called positional FEM and it is based on the minimum potential energy theorem written according to nodal positions, instead of displacements. The study herein presented considers the effects of geometric and material nonlinearities. Related to dynamic problems, a comparison between different time integration algorithms is performed. The formulation is extended to impact problems between trusses and rigid wall, where the nodal positions are constrained considering nullpenetration condition. In addition, it is presented a thermodynamically consistent formulation, based on the first and second law of thermodynamics and the Helmholtz free-energy for analyzing dynamic problems of truss structures with thermoelastic and thermoplastic behavior. The numerical results of the proposed formulation are compared with examples found in the literature.

Projeto iterativo, simulação, análise e otimização de máquina a relutância variável monofásica

Variable reluctance motors have been increasingly used as an alternative for variable speed and high speed drives in many industrial applications, due to many advantages like the simplicity of construction, robustness, and low cost. The most common applications in recent years are related to aeronautics, electric and hybrid vehicles and wind power generation. This paper explores the theory, operation, design procedures and analysis of a variable reluctance machine. An iterative design methodology is introduced and used to design a 1.25 kW prototype. For the analysis of the machine two methods are used, an analytical method and the finite element simulation. The results obtained by both methods are compared. The results of finite element simulation are used to determine the inductance profiles and torque of the prototype. The magnetic saturation is examined visually and numerically in four critical points of the machine. The data collected in the simulation allow the verification of design and operating limits for the prototype. Moreover, the behavior of the output quantities is analyzed (inductance, torque and magnetic saturation) by variation of physical dimensions of the motor. Finally, a multiobjective optimization using Differential Evolution algorithms and Genetic Algorithms for switched reluctance machine design is proposed. The optimized variables are rotor and stator polar arcs, and the goals are to maximize the average torque, the average torque per copper losses and the average torque per core volume. Finally, the initial design and optimized design are compared.

Análise do uso de escalas nas simulações de processos de estampagem

The purpose of using software based on numerical approximations for metal forming is given by the need to ensure process efficiency in order to get high quality products at lowest cost and shortest time. This study uses the theory of similitude in order to develop a technique capable of simulating the stamping process of a metal sheet, obtaining results close to the real values, with shorter processing times. The results are obtained through simulations performed in the finite element software STAMPACK®. This software uses the explicit integration method in time, which is usually applied to solve nonlinear problems involving contact, such as the metal forming processes. The technique was developed from a stamping model of a square box, simulated with four different scale factors, two higher and two smaller than the real scale. The technique was validated with a bending model of a welded plate, which had a high simulation time. The application of the technique allowed over 50% of decrease in the time of simulation. The results for the application of the scale technique for forming plates were satisfactory, showing good quantitative results related to the decrease of the total time of simulation. Finally, it is noted that the decrease in simulation time is only possible with the use of two related scales, the geometric and kinematic scale. The kinematic scale factors should be used with caution, because the high speeds can cause dynamic problems and could influence the results of the simulations.

A influência do amortecimento viscoelástico na estabilidade aeroelástica de painéis aeronáuticos

Since the creation of supersonic vehicles, during the Second World War, the engineers have given special attention to the interaction between the aerodynamic efforts and the structures of the aircrafts due to a highly destructive phenomenon called flutter in aeronautical panel. Flutter in aeronautical panels is a self-excited aeroelastic phenomenon, which can occurs during supersonic flights due to dynamic instability of inertia, elastic and aerodynamic forces of the system. In the flutter condition, when the critical aerodynamic pressure is reached, the vibration amplitudes of the panel become dynamically unstable and increase exponentially with time, affecting significantly the fatigue life of the existing aeronautical components. Thus, in this paper, the interest is to investigate the possibility of reducing the effects of the supersonic aeroelastic instability of rectangular plates by applying passive constrained viscoelastic layers. The rationale for such study is the fact that as the addition of viscoelastic materials provides decreased vibration amplitudes it becomes important to quantify the suppression of plate flutter coalescence modes that can be obtained. Moreover, despite the fact that much research on the suppression of panel flutter has been carried out by using passive, semi-active and active control techniques, very few of them are adapted to deal with the problem of estimating the flutter speeds of viscoelastic systems, since they must conveniently account for the frequency- and temperature-dependent behavior of the viscoelastic material. In this context, two different model of viscoelastic material are developed and applied to the model of sandwich plate by using finite elements. After the presentation of the theoretical foundations of the methodology, the description of a numerical study on the flutter analysis of a three-layer sandwich plate is addressed.

Strength of the Iberian intraplate lithosphere: Cenozoic deformations and seismicity

Previous studies about the strength of the lithosphere in the Iberia centre fail to resolve the depth of earthquakes because of the rheological uncertainties. Therefore, new contributions are considered (the crustal structure from a density model) and several parameters (tectonic regime, mantle rheology, strain rate) are checked in this paper to properly examine the role of lithospheric strength in the intraplate seismicity and the Cenozoic evolution. The strength distribution with depth, the integrated strength, the effective elastic thickness and the seismogenic thickness have been calculated by a finite element modelling of the lithosphere across the Central System mountain range and the bordering Duero and Madrid sedimentary basins. Only a dry mantle under strike-slip/extension and a strain rate of 10-15 s-1, or under extension and 10-16 s-1, causes a strong lithosphere. The integrated strength and the elastic thickness are lower in the mountain chain than in the basins. These anisotropies have been maintained since the Cenozoic and determine the mountain uplift and the biharmonic folding of the Iberian lithosphere during the Alpine deformations. The seismogenic thickness bounds the seismic activity in the upper–middle crust, and the decreasing crustal strength from the Duero Basin towards the Madrid Basin is related to a parallel increase in Plio–Quaternary deformations and seismicity. However, elasto–plastic modelling shows that current African–Eurasian convergence is resolved elastically or ductilely, which accounts for the low seismicity recorded in this region.

Evaluación de tenacidad a la fractura e integridad estructural de tubos de generadores de vapor nucleares.

Dado el impacto negativo asociado a la ocurrencia de fallas en tubos de generadores de vapor (TGVs) en centrales nucleares, el estudio de la integridad estructural de éstos ha comenzado a recibir mayor atención recientemente. Diversas metodologías basadas en análisis de carga límite han sido propuestas para asegurar la integridad estructural de los tubos, según los requerimientos establecidos por las autoridades regulatorias. Éstas han conducido, sin embargo, a la definición de criterios de reparación o taponado de TGVs excesivamente conservativos. Por lo tanto, con el objetivo de reducir la cantidad de tubos innecesariamente removidos de servicio, nuevos criterios de evaluación de integridad han sido propuestos recientemente en la literatura. En este contexto, la mecánica de fractura elastoplástica se presenta como una alternativa para la evaluación de la integridad de TGVs, requiriéndose dos elementos para su aplicación: la estimación de la fuerza impulsora en términos del parámetro elastoplástico (por ejemplo, la integral J) y la medición experimental de la tenacidad a la fractura del material de los tubos (por ejemplo, a través de la curva de resistencia J-R). Este trabajo presenta el desarrollo de técnicas experimentales no normalizadas para la determinación de curvas J-R para TGVs con fisuras pasantes circunferenciales y longitudinales. Debido a las dimensiones reducidas de los TGVs, diferentes probetas no normalizadas fueron propuestas. Además, en los ensayos se utilizaron condiciones de carga de tracción y flexión con el objetivo de modelar más adecuadamente los estados tensionales y las condiciones de constraint reales en TGVs. Los valores de la integral J fueron estimados utilizando el método del factor η. La aptitud del método fue evaluada a partir de simulaciones numéricas de los ensayos propuestos mediante análisis elastoplásticos con la técnica de elementos finitos. Se encontró que condiciones de mayor constraint asociadas con fisuras profundas y cargas de flexión favorecen la validez del método del factor η, mientras que configuraciones de menor constraint dan como resultado factores η que exhiben una mayor dependencia con el nivel de carga aplicada. También se observó que los factores η basados en la apertura de la boca de la fisura (Crack Mouth Opening Displacement o CMOD) presentan una dependencia mucho menor con el nivel de carga respecto a los factores η definidos a partir del desplazamiento del punto de aplicación de la carga (Load Line Displacement o LLD). Se presentan los valores del factor η para las probetas estudiadas con fisuras profundas (a/W ≥ 0,40). Se realizaron ensayos de fractura a temperatura ambiente y 300 °C con probetas obtenidas de TGVs nucleares fabricados a partir de las aleaciones 690 (Ni: 61; Cr: 29; Fe: 8,95, % en peso) y 800 (Ni: 33; Cr: 21,6; Fe: 42,2, % en peso). Durante los ensayos de fractura a temperatura ambiente, la extensión estable de fisura fue medida mediante una técnica óptica utilizando un microscopio digital. Para estos ensayos también se aplicó el método de normalización que propone la norma ASTM E1820-15 en el Anexo 15, encontrándose una buena coincidencia entre las longitudes estimadas por éste y las medidas ópticamente. De esta manera, el método de normalización fue utilizado para los ensayos a alta temperatura. Los resultados experimentales mostraron que ambos materiales tienen elevadas tenacidades a la fractura, siendo la aleación 800 la que presentó curvas J-R más elevadas que la aleación 690 tanto para fisuras circunferenciales como longitudinales. Las curvas J-R para ambas aleaciones mostraron un efecto marcado con la orientación de la fisura, es decir que existe una importante anisotropía en las propiedades de fractura: las fisuras circunferenciales presentaron curvas J-R más elevadas que las fisuras longitudinales. El nivel de constraint desarrollado en los ensayos, dado por las condiciones de carga de tracción y flexión, evidenció poco efecto sobre las curvas J-R para probetas con fisuras profundas (a/W ~ 0,50). A su vez, la temperatura de ensayo (temperatura ambiente y 300 °C) presentó un efecto prácticamente nulo para ambas aleaciones. Usando las propiedades de fractura obtenidas en este trabajo, la metodología FAD (Failure Assessment Diagram) fue propuesta y utilizada para la predicción de las condiciones de falla de TGVs fisurados para diferentes geometrías de fisura y condiciones de carga. La comparación entre análisis teóricos y datos experimentales muestra la potencialidad del FAD como una metodología capaz de predecir adecuadamente las fallas de estos componentes.

Estudo do efeito do carregamento de pulso em viga

Nesse trabalho é estudado uma viga isotrópica com as condições de contorno engastada-apoiada sujeita a uma carga dinâmica de pulso. Inicialmente a equação de equilíbrio dinâmica da viga é desenvolvida e as frequências naturais com seus respectivos modos de vibração são obtidos. Para a obtenção da resposta no tempo, não amortecida e amortecida, são utilizados a Integral de Duhamel e o Método dos Elementos Finitos (MEF) através da utilização do programa comercial Abaqus ®. Os resultados obtidos através das duas metodologias são comparados e observa-se uma boa convergência entre eles.

Evaluación de tenacidad a la fractura e integridad estructural de tubos de generadores de vapor nucleares.

Dado el impacto negativo asociado a la ocurrencia de fallas en tubos de generadores de vapor (TGVs) en centrales nucleares, el estudio de la integridad estructural de éstos ha comenzado a recibir mayor atención recientemente. Diversas metodologías basadas en análisis de carga límite han sido propuestas para asegurar la integridad estructural de los tubos, según los requerimientos establecidos por las autoridades regulatorias. Éstas han conducido, sin embargo, a la definición de criterios de reparación o taponado de TGVs excesivamente conservativos. Por lo tanto, con el objetivo de reducir la cantidad de tubos innecesariamente removidos de servicio, nuevos criterios de evaluación de integridad han sido propuestos recientemente en la literatura. En este contexto, la mecánica de fractura elastoplástica se presenta como una alternativa para la evaluación de la integridad de TGVs, requiriéndose dos elementos para su aplicación: la estimación de la fuerza impulsora en términos del parámetro elastoplástico (por ejemplo, la integral J) y la medición experimental de la tenacidad a la fractura del material de los tubos (por ejemplo, a través de la curva de resistencia J-R). Este trabajo presenta el desarrollo de técnicas experimentales no normalizadas para la determinación de curvas J-R para TGVs con fisuras pasantes circunferenciales y longitudinales. Debido a las dimensiones reducidas de los TGVs, diferentes probetas no normalizadas fueron propuestas. Además, en los ensayos se utilizaron condiciones de carga de tracción y flexión con el objetivo de modelar más adecuadamente los estados tensionales y las condiciones de constraint reales en TGVs. Los valores de la integral J fueron estimados utilizando el método del factor η. La aptitud del método fue evaluada a partir de simulaciones numéricas de los ensayos propuestos mediante análisis elastoplásticos con la técnica de elementos finitos. Se encontró que condiciones de mayor constraint asociadas con fisuras profundas y cargas de flexión favorecen la validez del método del factor η, mientras que configuraciones de menor constraint dan como resultado factores η que exhiben una mayor dependencia con el nivel de carga aplicada. También se observó que los factores η basados en la apertura de la boca de la fisura (Crack Mouth Opening Displacement o CMOD) presentan una dependencia mucho menor con el nivel de carga respecto a los factores η definidos a partir del desplazamiento del punto de aplicación de la carga (Load Line Displacement o LLD). Se presentan los valores del factor η para las probetas estudiadas con fisuras profundas (a/W ≥ 0,40). Se realizaron ensayos de fractura a temperatura ambiente y 300 °C con probetas obtenidas de TGVs nucleares fabricados a partir de las aleaciones 690 (Ni: 61; Cr: 29; Fe: 8,95, % en peso) y 800 (Ni: 33; Cr: 21,6; Fe: 42,2, % en peso). Durante los ensayos de fractura a temperatura ambiente, la extensión estable de fisura fue medida mediante una técnica óptica utilizando un microscopio digital. Para estos ensayos también se aplicó el método de normalización que propone la norma ASTM E1820-15 en el Anexo 15, encontrándose una buena coincidencia entre las longitudes estimadas por éste y las medidas ópticamente. De esta manera, el método de normalización fue utilizado para los ensayos a alta temperatura. Los resultados experimentales mostraron que ambos materiales tienen elevadas tenacidades a la fractura, siendo la aleación 800 la que presentó curvas J-R más elevadas que la aleación 690 tanto para fisuras circunferenciales como longitudinales. Las curvas J-R para ambas aleaciones mostraron un efecto marcado con la orientación de la fisura, es decir que existe una importante anisotropía en las propiedades de fractura: las fisuras circunferenciales presentaron curvas J-R más elevadas que las fisuras longitudinales. El nivel de constraint desarrollado en los ensayos, dado por las condiciones de carga de tracción y flexión, evidenció poco efecto sobre las curvas J-R para probetas con fisuras profundas (a/W ~ 0,50). A su vez, la temperatura de ensayo (temperatura ambiente y 300 °C) presentó un efecto prácticamente nulo para ambas aleaciones. Usando las propiedades de fractura obtenidas en este trabajo, la metodología FAD (Failure Assessment Diagram) fue propuesta y utilizada para la predicción de las condiciones de falla de TGVs fisurados para diferentes geometrías de fisura y condiciones de carga. La comparación entre análisis teóricos y datos experimentales muestra la potencialidad del FAD como una metodología capaz de predecir adecuadamente las fallas de estos componentes.