Estudo numérico do comportamento de muros de arrimo em alvenaria estrutural de blocos vazados

Usually masonry structures has low tension strength, hence the design to flexural efforts can results in high reinforcement ratio, specification of high unit and prism strength, structural members with larger section dimensions and modification in structural arrangement to be possible to use masonry members. The main objective of this study is to evaluate the stiffness, the efforts distribution and the effect of horizontal elements (girders) and vertical elements (counterforts) distribution on the behavior of masonry blocks retaining walls. For this purpose, numerical modeling was performed on typical retaining wall arrangements by varying the amount and placement of horizontal and vertical elements, beyond includes elements simulating the reactions of the soil supporting the foundation of the wall. The numerical modeling also include the macro modeling strategy in which the units, mortar and grout are discretized by a standard volume that represents the masonry elastic behavior. Also, numerical model results were compared with those ones of simplified models usually adopted in bending design of masonry elements. The results show horizontal displacements, principal and shear stresses distribution, and bending moments diagrams. From the analysis it was concluded that quantity and manner of distribution of the girders are both important factors to the panel flexural behavior, the inclusion of the foundation changed significantly the behavior of the wall, especially the horizontal displacements, and has been proposed a new way of considering the flanges section of the counterforts

Modelagem numérica de estruturas de contenção atirantadas em areia

A numerical study on the behavior of tied-back retaining walls in sand, using the finite element method (FEM) is presented. The analyses were performed using the software Plaxis 2D, and were focused on the development of horizontal displacements, horizontal stresses, shear forces and bending moments in the structure during the construction process. Emphasis was placed on the evaluation of wall embedment, tie-back horizontal spacing, wall thickness, and free anchor length on wall behavior. A representative soil profile of a specific region at the City of Natal, Brazil, was used in the numerical analyses. New facilities built on this region often include retaining structures of the same type studied herein. Soil behavior was modeled using the Mohr-Coulomb constitutive model, whereas the structural elements were modeled using the linear elastic model. Shear strength parameters of the soil layers were obtained from direct shear test results conducted with samples collected at the studied site. Deformation parameters were obtained from empirical correlations from SPT test results carried out on the studied site. The results of the numerical analyses revealed that the effect of wall embedment on the investigated parameters is virtually negligible. Conversely, the tie-back horizontal spacing plays an important role on the investigated parameters. The results also demonstrated that the wall thickness significantly affects the wall horizontal displacements, and the shear forces and bending moments within the retaining structure. However, wall thickness was not found to influence horizontal stresses in the structure

Modelagem computacional para avaliação da interação entre painéis de alvenaria e estrutura de suporte em concreto armado

The developed study proposes a new computer modeling efficient and easy to apply in usual project situations to evaluate the interaction between masonry panels and support structure. The proposed model simulates the behavior of the wall exclusively using frame finite elements, thus compounding an equivalent frame. The validation was performed in two ways: firstly, through the analysis of various panels of generic plans, comparing the results obtained from equivalent frame model with the ones from a reference model, which uses shell finite elements in discretization of the walls; and in a second step, comparing with the results of the experimental model of Rosenhaupt. The analyzes considered the linear elastic behavior for materials and consisted basically in the evaluation of vertical displacements and efforts in support beams, and tensions at the base of walls. Was also evaluated, from flat and threedimensional modeling of some walls from a real project, important aspects of the wall-beam interaction, e.g.: the presence of openings of doors and windows, arranged in any position; conditions of support and linking of beams; interference of moorings between walls; and consideration of wind action. The analysis of the achieved results demonstrated the efficiency of the proposed modeling, since they have very similar aspects in the distribution of stresses and efforts, always with intensities slightly larger than those of the reference and experimental models.

Modelos para análise e dimensionamento de painéis de alvenaria estrutural

The main objective of this thesis was the study of bracing panels of structural masonry, by applying the Finite Element Method and Strut and Tie Method. It was analyzed the following aspects: the effect of orthotropy on the behavior of the panels; distribution of horizontal forces between panels for buildings; comparison between Equivalent Frame and Finite Elements models; panels design with the Strut and Tie Method. The results showed that one should not disregard the orthotropy, otherwise this can lead to models stiffer than the real. Regarding the distribution of horizontal forces, showed that the disregard of lintels and shear deformation leads to significant differences in the simplified model. The results showed also that the models in Finite Element and Equivalent Frame exhibit similar behavior in respect to stiffness of panels and stress distribution over the sessions requested. It was discussing criteria for designing Strut and Tie Method models in one floor panels. Then, the theoretical strength these panels was compared with the rupture strength of panels tested in the literature. The theoretical maximum strength were always less than the rupture strength of the panels obtained in tests, due to the fact that the proposed model cannot represent the behavior of the masonry after the start of the panel cracking due to plasticization of the reinforcement

Modelagem numérica de escoamento em padrão anular com transferência de calor e massa

Annular flow is the prevailing pattern in transport and energy conversion systems and therefore, one of the most important patterns in multiphase flow in ducts. The correct prediction of the pressure gradient and heat transfer coefficient is essential for optimizing the system s capacity. The objective of this work is to develop and implement a numerical algorithm capable of predicting hydrodynamic and thermal characteristics for upflow, vertical, annular flow. The numerical algorithm is then complemented with the physical modeling of phenomena that occurs in this flow pattern. These are, turbulence, entrainment and deposition and phase change. For the development of the numerical model, axial diffusion of heat and momentum is neglected. In this way the time-averaged equations are solved in their parabolic form obtaining the velocity and temperature profiles for each axial step at a time, together with the global parameters, namely, pressure gradient, mean film thickness and heat transfer coefficient, as well as their variation in the axial direction. The model is validated for the following conditions: fully-developed laminar flow with no entrainment; fully developed laminar flow with heat transfer, fully-developed turbulent flow with entrained drops, developing turbulent annular flow with entrained drops, and turbulent flow with heat transfer and phase change

Modelagem numérica de uma estrutura de contenção de estacas espaçadas atirantadas em areia

Retaining walls design involves factors such as plastification, loading and unloading, pre-stressing, excessive displacements and earth and water thrust. Furthermore, the interaction between the retained soil and the structure is rather complex and hard to predict. Despite the advances in numerical simulation techniques and monitoring of forces and displacements with field instrumentation, design projects are still based on classical methods, whose simplifying assumptions may overestimate structural elements of the retaining wall. This dissertation involves a three-dimensional numerical study on the behavior of a retaining wall using the finite element method (FEM). The retaining wall structure is a contiguous bored pile wall with tie-back anchors. The numerical results were compared to data obtained from field instrumentation. The influence of the position of one or two layers of anchors and the effects of the construction of a slab bounded at the top of the retaining wall was evaluated. Furthermore, this study aimed at investigating the phenomenon of arching in the soil behind the wall. Arching was evaluated by analyzing the effects of pile spacing on horizontal stresses and displacements. Parametric analysis with one layers of anchors showed that the smallest horizontal displacements of the structure were achieved for between 0.3 and 0.5 times the excavation depth. Parametric analyses with two anchor layers showed that the smallest horizontal displacements were achieve for anchors positioned in depths of 0.4H and 0.7H. The construction of a slab at the top of the retaining wall decreased the horizontal displacements by 0.14% times the excavation depth as compared to analyses without the slab. With regard to the arching , analyzes showed an optimal range of spacing between the faces of the piles between 0.4 and 0.6 times the diameter of the pile

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.

Desenvolvimento de dobras e falhas em ambiente distensional: Aplicação da modelagem física

The geological modeling allows, at laboratory scaling, the simulation of the geometric and kinematic evolution of geological structures. The importance of the knowledge of these structures grows when we consider their role in the creation of traps or conduits to oil and water. In the present work we simulated the formation of folds and faults in extensional environment, through physical and numerical modeling, using a sandbox apparatus and MOVE2010 software. The physical modeling of structures developed in the hangingwall of a listric fault, showed the formation of active and inactive axial zones. In consonance with the literature, we verified the formation of a rollover between these two axial zones. The crestal collapse of the anticline formed grabens, limited by secondary faults, perpendicular to the extension, with a curvilinear aspect. Adjacent to these faults we registered the formation of transversal folds, parallel to the extension, characterized by a syncline in the fault hangingwall. We also observed drag folds near the faults surfaces, these faults are parallel to the fault surface and presented an anticline in the footwall and a syncline hangingwall. To observe the influence of geometrical variations (dip and width) in the flat of a flat-ramp fault, we made two experimental series, being the first with the flat varying in dip and width and the second maintaining the flat variation in width but horizontal. These experiments developed secondary faults, perpendicular to the extension, that were grouped in three sets: i) antithetic faults with a curvilinear geometry and synthetic faults, with a more rectilinear geometry, both nucleated in the base of sedimentary pile. The normal antithetic faults can rotate, during the extension, presenting a pseudo-inverse kinematics. ii) Faults nucleated at the top of the sedimentary pile. The propagation of these faults is made through coalescence of segments, originating, sometimes, the formation of relay ramps. iii) Reverse faults, are nucleated in the flat-ramp interface. Comparing the two models we verified that the dip of the flat favors a differentiated nucleation of the faults at the two extremities of the mater fault. V These two flat-ramp models also generated an anticline-syncline pair, drag and transversal folds. The anticline was formed above the flat being sub-parallel to the master fault plane, while the syncline was formed in more distal areas of the fault. Due the geometrical variation of these two folds we can define three structural domains. Using the physical experiments as a template, we also made numerical modeling experiments, with flat-ramp faults presenting variation in the flat. Secondary antithetic, synthetic and reverse faults were generated in both models. The numerical modeling formed two folds, and anticline above the flat and a syncline further away of the master fault. The geometric variation of these two folds allowed the definition of three structural domains parallel to the extension. These data reinforce the physical models. The comparisons between natural data of a flat-ramp fault in the Potiguar basin with the data of physical and numerical simulations, showed that, in both cases, the variation of the geometry of the flat produces, variation in the hangingwall geometry

Influência dos modelos digitais de terreno na modelagem de circulação hidrodinâmica 2DH: um estudo de caso

Digital Elevation Models (DEM) are numerical representations of a portion of the earth surface. Among several factors which affect the quality of a DEM, it should be emphasized the attention on the input data and the choice of the interpolating algorithm. On the other hand, several numerical models are used nowadays to characterize nearshore hydrodynamics and morphological changes in coastal areas, whose validation is based on field data collection. Independent on the complexity of the physical processes which are modeled, little attention has been given to the intrinsic bathymetric interpolation built within the numerical models of the specific application. Therefore, this study aims to investigate and to quantify the influence of the bathymetry, as obtained by a DEM, on the hydrodynamic circulation model at a coastal stretch, off the coast of the State of Rio Grande do Norte, Northeast Brazil. This coastal region is characterized by strong hydrodynamic and littoral processes, resulting in a very dynamic morphology with shallow coastal bathymetry. Important economic activities, such as oil exploitation and production, fisheries, salt ponds, shrimp farms and tourism, also bring impacts upon the local ecosystems and influence themselves the local hydrodynamics. This fact makes the region one of the most important for the development of the State, but also enhances the possibility of serious environmental accidents. As a hydrodynamic model, SisBaHiA® - Environmental Hydrodynamics System ( Sistema Básico de Hidrodinâmica Ambiental ) was chosen, for it has been successfully employed at several locations along the Brazilian coast. This model was developed at the Coastal and Oceanographical Engineering Group of the Ocean Engineering Program at the Federal University of Rio de Janeiro. Several interpolating methods were tested for the construction of the DEM, namely Natural Neighbor, Kriging, Triangulation with Linear Interpolation, Inverse Distance to a Power, Nearest Neighbor, and Minimum Curvature, all implemented within the software Surfer®. The bathymetry which was used as reference for the DEM was obtained from nautical charts provided by the Brazilian Hydrographic Service of the Brazilian Navy and from a field survey conducted in 2005. Changes in flow velocity and free surface elevation were evaluated under three aspects: a spatial vision along three profiles perpendicular to the coast and one profile longitudinal to the coast as shown; a temporal vision from three central nodes of the grid during 30 days; a hodograph analysis of components of speed in U and V, by different tidal cycles. Small, but negligible, variations in sea surface elevation were identified. However, the differences in flow and direction of velocities were significant, depending on the DEM

Análise de antenas de microfita com patches circulares sobre substratos anisotrópicos usando o método dos potenciais de Hertz

This work consists on the theoretical and numerical analysis of some properties of circular microstrip patch antennas on isotropic and uniaxial anisotropic substrates. For this purpose, a full wave analysis is performed, using Hertz Vector Potentials method in the Hankel Transform domain. In the numerical analysis, the moment method is also used in order to determine some characteristics of the antenna, such as: resonant frequency and radiation pattern. The definition of Hertz potentials in the Hankel domain is used in association with Maxwell´s equations and the boundary conditions of the structures to obtain the Green´s functions, relating the components of the current density on the patch and the tangential electric field components. Then, the Galerkin method is used to generate a matrix equation whose nontrivial solution is the complex resonant frequency of the structure. In the analysis, a microstrip antenna with only one isotropic dielectric layer is initially considered. For this structure, the effect of using superconductor patches is also analyzed. An analysis of a circular microstrip antenna on an uniaxial anisotropic dielectric layer is performed, using the Hertz vector potentials oriented along the optical axis of the material, that is perpendicular to the microstrip ground plane. Afterwards, the circular microstrip antenna using two uniaxial anisotropic dielectric layers is investigated, considering the particular case in which the inferior layer is filled by air. In this study, numerical results for resonant frequency and radiation pattern for circular microstrip antennas on isotropic and uniaxial anisotropic substrates are presented and compared with measured and calculated results found in the literature

Simulação numérica do fluxo hídrico subterrâneo na bacia hidrográfica do Rio Pitimbu RN

Urban centers in Pitimbu Watershed use significant groundwater sources for public supply. Therefore, studies in Dunas Barreiras aquifer are relevant to expand knowledge about it and help manage water resources in the region. An essential tool for this management is the numerical modeling of groundwater flow. In this work, we developed a groundwater flow model for Pitimbu Watershed, using the Visual Modflow, version 2.7.1., which uses finite difference method for solving the govern equation of the dynamics of groundwater flow. We carried out the numerical simulation of steady-state model for the entire region of the basin. The model was built in the geographical, geomorphological and hydrogeological study of the area, which defined the boundary conditions and the parameters required for the numerical calculation. Owing to unavailability of current data based on monitoring of the aquifer it was not possible to calibrate the model. However, the simulation results showed that the overall water balance approached zero, therefore satisfying the equation for the three-dimensional behavior of the head water in steady state. Variations in aquifer recharge data were made to verify the impact of this contribution on the water balance of the system, especially in the scenario in which recharge due to drains and sinks was removed. According to the results generated by Visual Modflow occurred significantly hydraulic head lowering, ranging from 16,4 to 82 feet of drawdown. With the results obtained, it can be said that modeling is performed as a valid tool for the management of water resources in Pitimbu River Basin, and to support new studies

A utilização da modelagem física no estudo do desenvolvimento, geometria e cinemática de falhas distensionais

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

Análise numérica de dutos enterrados sujeitos a elevação e recalque localizado

Problems associated to longitudinal interactions in buried pipelines are characterized as three-dimensional and can lead to different soil-pipe issues. Despite the progress achieved in research on buried pipelines, little attention has been given to the three-dimensional nature of the problem throughout the last decades. Most of researches simplify the problem by considering it in plane strain condition. This dissertation aims to present a study on the behavior of buried pipelines under local settlement or elevation, using three-dimensional simulations. Finite element code Plaxis 3D was used for the simulations. Particular aspects of the numerical modeling were evaluated and parametric analyzes were performed, was investigated the effects of soil arching in three-dimensional form. The main variables investigated were as follows: relative density, displacement of the elevation or settlement zone, elevated zone size, height of soil cover and pipe diameter/thickness ratio. The simulations were performed in two stages. The first stage was involved the validation of the numerical analysis using the physical models put forward by Costa (2005). In the second stage, numerical analyzes of a full-scale pipeline subjected to a localized elevation were performed. The obtained results allowed a detailed evaluation of the redistribution of stresses in the soil mass and the deflections along the pipe. It was observed the reduction of stresses in the soil mass and pipe deflections when the height of soil cover was decreased on regions of the pipe subjected to elevation. It was also shown for the analyzed situation that longitudinal thrusts were higher than vi circumferential trusts and exceeded the allowable stresses and deflections. Furthermore, the benefits of minimizing stress with technical as the false trench, compressible cradle and a combination of both applied to the simulated pipeline were verified

Estudo das forças capilar, viscosa e gravitacional através de grupos adimensionais na simulação da drenagem gravitacional assistida por gás (GAGD)

The distribution and mobilization of fluid in a porous medium depend on the capillary, gravity, and viscous forces. In oil field, the processes of enhanced oil recovery involve change and importance of these forces to increase the oil recovery factor. In the case of gas assisted gravity drainage (GAGD) process is important to understand the physical mechanisms to mobilize oil through the interaction of these forces. For this reason, several authors have developed physical models in laboratory and core floods of GAGD to study the performance of these forces through dimensionless groups. These models showed conclusive results. However, numerical simulation models have not been used for this type of study. Therefore, the objective of this work is to study the performance of capillary, viscous and gravity forces on GAGD process and its influence on the oil recovery factor through a 2D numerical simulation model. To analyze the interplay of these forces, dimensionless groups reported in the literature have been used such as Capillary Number (Nc), Bond number (Nb) and Gravity Number (Ng). This was done to determine the effectiveness of each force related to the other one. A comparison of the results obtained from the numerical simulation was also carried out with the results reported in the literature. The results showed that before breakthrough time, the lower is the injection flow rate, oil recovery is increased by capillary force, and after breakthrough time, the higher is the injection flow rate, oil recovery is increased by gravity force. A good relationship was found between the results obtained in this research with those published in the literature. The simulation results indicated that before the gas breakthrough, higher oil recoveries were obtained at lower Nc and Nb and, after the gas breakthrough, higher oil recoveries were obtained at lower Ng. The numerical models are consistent with the reported results in the literature