GENERALIZED FINITE ELEMENT METHOD FOR NONLINEAR THREE-DIMENSIONAL ANALYSIS OF SOLIDS

The Generalized Finite Element Method (GFEM) is employed in this paper for the numerical analysis of three-dimensional solids tinder nonlinear behavior. A brief summary of the GFEM as well as a description of the formulation of the hexahedral element based oil the proposed enrichment strategy are initially presented. Next, in order to introduce the nonlinear analysis of solids, two constitutive models are briefly reviewed: Lemaitre`s model, in which damage and plasticity are coupled, and Mazars`s damage model suitable for concrete tinder increased loading. Both models are employed in the framework of a nonlocal approach to ensure solution objectivity. In the numerical analyses carried out, a selective enrichment of approximation at regions of concern in the domain (mainly those with high strain and damage gradients) is exploited. Such a possibility makes the three-dimensional analysis less expensive and practicable since re-meshing resources, characteristic of h-adaptivity, can be minimized. Moreover, a combination of three-dimensional analysis and the selective enrichment presents a valuable good tool for a better description of both damage and plastic strain scatterings.

GIS applied to geotechnical and environmental risk management in a Brazilian oil pipeline

The paper presents the development of a decision support system for the management of geotechnical and environmental risks in oil pipelines using a geographical information system. The system covers a 48.5 km long section of the So Paulo to Brasilia (OSBRA) oil pipeline, which crosses three municipalities in the northeast region of the So Paulo state (Brazil) and represents an area of 205.8 km(2). The spatial database was created using geo-processing procedures, surface and intrusive investigations and geotechnical reports. The risk assessment was based mainly on qualitative models (relative numeric weights and multicriteria decision analysis) and considered pluvial erosion, slope movements, soil corrosion and third party activities. The maps were produced at a scale of 1:10,000.

Interface shear strength of geosynthetics: Evaluation and analysis of inclined plane tests

The inclined plane test (IPT) is commonly performed to measure the interface shear strength between different materials as those used in cover systems of landfills. The test, when interpreted according to European test Standards provides the static interface friction angle, usually assumed for 50 mm displacement and denoted as phi(stat)(50). However, if interpreted considering the several phases of the sliding process, the test is capable of yielding more realistic information about the interface shear strength such as differentiating interfaces which exhibit the same value of phi(stat)(50) but different behavior for displacement less than 50 mm. In this paper, the IPT is used to evaluate the interface shear strength of some materials usually present in cover liner systems of landfill. The results of the tests were analyzed for both, the static and the dynamic phases of the sliding and were interpreted based on the static initial friction angle, phi(0), and the limit friction angle, phi(lim). It is shown that depending on the sliding behavior of the interfaces, phi(stat)(50), which is usually adopted as the designing parameter in stability analysis, can be larger than phi(0) and phi(lim). (C) 2009 Elsevier Ltd. All rights reserved.

AnSBBR Applied to the Treatment of Metalworking Fluid Wastewater: Effect of Organic and Shock Load

An investigation was performed regarding the application of a mechanically stirred anaerobic sequencing batch biofilm reactor containing immobilized biomass on inert polyurethane foam (AnSBBR) to the treatment of soluble metalworking fluids to remove organic matter and produce methane. The effect of increasing organic matter and reactor fill time, as well as shock load, on reactor stability and efficiency have been analyzed. The 5-L AnSBBR was operated at 30 A degrees C in 8-h cycles, agitation of 400 rpm, and treated 2.0 L effluent per cycle. Organic matter was increased by increasing the influent concentration (500, 1,000, 2,000, and 3,000 mg chemical oxygen demand (COD)/L). Fill times investigated were in the batch mode (fill time 10 min) and fed-batch followed by batch (fill time 4 h). In the batch mode, organic matter removal efficiencies were 87%, 86%, and 80% for influent concentrations of 500, 1,000, and 2,000 mgCOD/L (1.50, 3.12, and 6.08 gCOD/L.d), respectively. At 3,000 mgCOD/L (9.38 gCOD/L.d), operational stability could not be achieved. The reactor managed to maintain stability when a shock load twice as high the feed concentration was applied, evidencing the robustness of the reactor to potential concentration variations in the wastewater being treated. Increasing the fill time to 4 h did not improve removal efficiency, which was 72% for 2,000 mgCOD/L. Thus, gradual feeding did not improve organic matter removal. The concentration of methane formed at 6.08 gCOD/L was 5.20 mmolCH(4), which corresponded to 78% of the biogas composition. The behavior of the reactor during batch and fed-batch feeding could be explained by a kinetic model that considers organic matter consumption, production, and consumption of total volatile acids and methane production.

Analysis of PVC pipe-wall viscoelasticity during water hammer

This research work focuses on the analysis of hydraulic transients in polyvinyl chloride (PVC) pipes, which are characterized by a viscoelastic rheological behavior. Transient pressure data were collected in a pipe rig consisting of a set of PVC pipes. The creep function of the PVC pipes was determined by using an inverse transient model based on collected transient pressure data and compared with that obtained by carrying out mechanical tensile tests of PVC pipe specimens. The numerical results obtained from the transient solver have shown that the attenuation, dispersion, and shape of transient pressures were well described. The incorporation of the viscoelastic mechanical behavior in the hydraulic transient model has provided an excellent fitting between numerical results and observed data. Calibrated creep function based on inverse analysis fit the one determined by mechanical tests well, which emphasized the importance of pipe-wall viscoelasticity in hydraulic transients in PVC pipes.

Anaerobic sequencing batch biofilm reactor applied to automobile industry wastewater treatment: Volumetric loading rate and feed strategy effects

This paper presents a technological viability study of wastewater treatment in an automobile industry by an anaerobic sequencing batch biofilm reactor containing immobilized biomass (AnSBBR) with a draft tube. The reactor was operated in 8-h cycles, with agitation of 400 rpm, at 30 degrees C and treating 2.0 L wastewater per cycle. Initially the efficiency and stability of the reactor were studied when supplied with nutrients and alkalinity. Removal efficiency of 88% was obtained at volumetric loading rate (VLR) of 3.09 mg COD/L day. When VLR was increased to 6.19 mg COD/L day the system presented stable operation with reduction in efficiency of 71%. In a second stage the AnSBBR was operated treating wastewater in natura, i.e., without nutrients supplementation, only with alkalinity, thereby changing feed strategy. The first strategy consisted in feeding 2.0 L batch wise (10 min), the second in feeding 1.0 L of influent batch wise (10 min) and an additional 1.0 L fed-batch wise (4 h), both dewatering 2.0 L of the effluent in 10 min. The third one maintained 1.0 L of treated effluent in the reactor, without discharging, and 1.0 L of influent was fed fed-batch wise (4 h) with dewatering 1.0 L of the effluent in 10 min. For all implemented strategies (VLR of 1.40, 2.57 and 2.61 mg COD/L day) the system presented stability and removal efficiency of approximately 80%. These results show that the AnSBBR presents operational flexibility, as the influent can be fed according to industry availability. In industrial processes this is a considerable advantage, as the influent may be prone to variations. Moreover, for all the investigated conditions the kinetic parameters were obtained from fitting a first-order model to the profiles of organic matter, total volatile acids and methane concentrations. Analysis of the kinetic parameters showed that the best strategy is feeding 1.0 L of influent batchwise (10 min) and 1.0 L fed-batch wise (4 h) in 8-h cycle. (c) 2007 Elsevier B.V. All rights reserved.

Specific methanogenic activity in an anaerobic-aerobic reactor applied to the treatment of domestic residual water

The specific methanogenic activity (SMA) test is an important tool for the monitoring of anaerobic digestion. This paper presents the behavior of the methanogenic archaea of an anaerobic sludge under different conditions of oxygenation in a fixed-bed anaerobic-aerobic reactor treating domestic sewage. The reactor was operated in a continuous manner under different liquid recycle ratios from aerobic to anaerobic zones in order to remove carbon and nitrogen. The application of the SMA test was adapted from several authors and the measurement of the accumulated methane in the reactor was carried out by means of gas chromatography. Methanogenic organisms were not inhibited by the presence of oxygen. In contrast, the values of CH, production rate by sludge exposed to oxygen were greater than those obtained for strictly anaerobic sludge.

An exploratory analysis of relationships between socioeconomic, land use, activity participation variables and travel patterns

This paper aims to find relations between the socioeconomic characteristics, activity participation, land use patterns and travel behavior of the residents in the Sao Paulo Metropolitan Area (SPMA) by using Exploratory Multivariate Data Analysis (EMDA) techniques. The variables influencing travel pattern choices are investigated using: (a) Cluster Analysis (CA), grouping and characterizing the Traffic Zones (17), proposing the independent variable called Origin Cluster and, (b) Decision Tree (DT) to find a priori unknown relations among socioeconomic characteristics, land use attributes of the origin TZ and destination choices. The analysis was based on the origin-destination home-interview survey carried out in SPMA in 1997. The DT application revealed the variables of greatest influence on the travel pattern choice. The most important independent variable considered by DT is car ownership, followed by the Use of Transportation ""credits"" for Transit tariff, and, finally, activity participation variables and Origin Cluster. With these results, it was possible to analyze the influence of a family income, car ownership, position of the individual in the family, use of transportation ""credits"" for transit tariff (mainly for travel mode sequence choice), activities participation (activity sequence choice) and Origin Cluster (destination/travel distance choice). (c) 2010 Elsevier Ltd. All rights reserved.

Analysis of two reinforced concrete pile caps with embedded sockets in the presence of a locking beam

The present research studies the behavior of reinforced concrete locking beams supported by two capped piles with the socket embedded; used as connections for pre-cast concrete structures. The effect provoked by locking the beam on the pile-caps when supported by the lateral socket walls was evaluated. Three-dimensional numerical analyses using software based on the finite element method (FEM) were developed considering the nonlinear physical behavior of the material. To evaluate the adopted software, a comparative analysis was made using the numerical and experimented results obtained from other software. In the pile caps studied, a variation in the wall thickness, socket interface, strut angle inclination and action on beam. The results show that the presence of a beam does not significantly change pile cap behavior and that the socket wall is able to effectively transfer the force from the beam to the pile caps. By the tensions on the bars of longitudinal reinforcement, it was possible to obtain the force on the tie and the strut angle inclination before the collapse of models. It was found that the angles present more inclinations than those used in the design, which was made based on a strut-and-tie model. More results are available at http://www.set.eesc.usp.br/pdf/download/2009ME_RodrigoBarros.pdf

Analytical solutions for geometrically nonlinear trusses

This paper presents an analytical method for analyzing trusses with severe geometrically nonlinear behavior. The main objective is to find analytical solutions for trusses with different axial forces in the bars. The methodology is based on truss kinematics, elastic constitutive laws and equilibrium of nodal forces. The proposed formulation can be applied to hyper elastic materials, such as rubber and elastic foams. A Von Mises truss with two bars made by different materials is analyzed to show the accuracy of this methodology.

A safety assessment methodology applied to CNS/ATM-based air traffic control system

In the last decades, the air traffic system has been changing to adapt itself to new social demands, mainly the safe growth of worldwide traffic capacity. Those changes are ruled by the Communication, Navigation, Surveillance/Air Traffic Management (CNS/ATM) paradigm, based on digital communication technologies (mainly satellites) as a way of improving communication, surveillance, navigation and air traffic management services. However, CNS/ATM poses new challenges and needs, mainly related to the safety assessment process. In face of these new challenges, and considering the main characteristics of the CNS/ATM, a methodology is proposed at this work by combining ""absolute"" and ""relative"" safety assessment methods adopted by the International Civil Aviation Organization (ICAO) in ICAO Doc.9689 [14], using Fluid Stochastic Petri Nets (FSPN) as the modeling formalism, and compares the safety metrics estimated from the simulation of both the proposed (in analysis) and the legacy system models. To demonstrate its usefulness, the proposed methodology was applied to the ""Automatic Dependent Surveillance-Broadcasting"" (ADS-B) based air traffic control system. As conclusions, the proposed methodology assured to assess CNS/ATM system safety properties, in which FSPN formalism provides important modeling capabilities, and discrete event simulation allowing the estimation of the desired safety metric. (C) 2011 Elsevier Ltd. All rights reserved.

Acoustic Simulation of a Special Switched Reluctance Drive by Means of Field-Circuit Coupling and Multiphysics Simulation

The approach presented in this paper consists of an energy-based field-circuit coupling in combination with multi-physics simulation of the acoustic radiation of electrical machines. The proposed method is applied to a special switched reluctance motor with asymmetric pole geometry to improve the start-up torque. The pole shape has been optimized, subject to low torque ripple, in a previous study. The proposed approach here is used to analyze the impact of the optimization on the overall acoustic behavior. The field-circuit coupling is based on a temporary lumped-parameter model of the magnetic part incorporated into a circuit simulation based on the modified nodal analysis. The harmonic force excitation is calculated by means of stress tensor computation, and it is transformed to a mechanical mesh by mapping techniques. The structural dynamic problem is solved in the frequency domain using a finite-element modal analysis and superposition. The radiation characteristic is obtained from boundary element acoustic simulation. Simulation results of both rotor types are compared, and measurements of the drive are presented.

A geometrical approach of 3-D FEA for educational purposes applied to electrostatic fields

A geometrical approach of the finite-element analysis applied to electrostatic fields is presented. This approach is particularly well adapted to teaching Finite Elements in Electrical Engineering courses at undergraduate level. The procedure leads to the same system of algebraic equations as that derived by classical approaches, such as variational principle or weighted residuals for nodal elements with plane symmetry. It is shown that the extension of the original procedure to three dimensions is straightforward, provided the domain be meshed in first-order tetrahedral elements. The element matrices are derived by applying Maxwell`s equations in integral form to suitably chosen surfaces in the finite-element mesh.

Efficient Modeling of Thin Wires in a Lossy Medium by Finite Elements Applied to Grounding Systems

A procedure is proposed to accurately model thin wires in lossy media by finite element analysis. It is based on the determination of a suitable element width in the vicinity of the wire, which strongly depends on the wire radius to yield accurate results. The approach is well adapted to the analysis of grounding systems. The numerical results of the application of finite element analysis with the suitably chosen element width are compared with both analytical results and those computed by a commercial package for the analysis of grounding systems, showing very good agreement.

A hybrid-mixed finite element formulation for the geometrically exact analysis of three-dimensional framed structures

This paper addresses the development of a hybrid-mixed finite element formulation for the quasi-static geometrically exact analysis of three-dimensional framed structures with linear elastic behavior. The formulation is based on a modified principle of stationary total complementary energy, involving, as independent variables, the generalized vectors of stress-resultants and displacements and, in addition, a set of Lagrange multipliers defined on the element boundaries. The finite element discretization scheme adopted within the framework of the proposed formulation leads to numerical solutions that strongly satisfy the equilibrium differential equations in the elements, as well as the equilibrium boundary conditions. This formulation consists, therefore, in a true equilibrium formulation for large displacements and rotations in space. Furthermore, this formulation is objective, as it ensures invariance of the strain measures under superposed rigid body rotations, and is not affected by the so-called shear-locking phenomenon. Also, the proposed formulation produces numerical solutions which are independent of the path of deformation. To validate and assess the accuracy of the proposed formulation, some benchmark problems are analyzed and their solutions compared with those obtained using the standard two-node displacement/ rotation-based formulation.