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.

Effect of bituminous impregnation on nonwoven geotextiles tensile and permeability properties

Geosynthetics interlayer systems are effective techniques to control reflective cracking in damaged pavements. It comprises the inclusion of nonwoven geotextiles between the damaged layer and the new overlay of the pavement to reduce the propagation of cracks and to extend pavement life. However, the success of this technique depends directly on the understanding of the geotextile`s behavior when impregnated with asphalt This paper evaluates different nonwoven geotextiles frequently used in anti-reflective cracking systems, focusing on initial stiffness gain and permeability reduction after asphalt impregnation. Fresh and impregnated samples of polyester and polypropylene nonwoven geotextiles were tested. Cationic rapid setting emulsified asphalt was used as asphalt binder. Wide-width tensile tests were carried out based on the specification of ABNT - NBR 12824 (1993). Water vapor transmission tests were conducted according to ASTM E 96M (2005). Results of tensile tests on impregnated geotextiles showed a significant increase on tensile strength values, probably due to the inter contact of the fibers. Results also showed high increase in strength values at strain levels less than 0.05% and decrease on stiffness gains with increase of strains. Water vapor transmission tests demonstrated that cationic asphalt emulsion applied on nonwoven geotextiles allows a drastic reduction in permeability values to turn nonwoven geotextiles into a low permeability barrier. (C) 2010 Elsevier Ltd. All rights reserved.

Determination of in situ uniaxial compressive strength of coal seams based on geophysical data

The central issue for pillar design in underground coal mining is the in situ uniaxial compressive strength (sigma (cm)). The paper proposes a new method for estimating in situ uniaxial compressive strength in coal seams based on laboratory strength and P wave propagation velocity. It describes the collection of samples in the Bonito coal seam, Fontanella Mine, southern Brazil, the techniques used for the structural mapping of the coal seam and determination of seismic wave propagation velocity as well as the laboratory procedures used to determine the strength and ultrasonic wave velocity. The results obtained using the new methodology are compared with those from seven other techniques for estimating in situ rock mass uniaxial compressive strength.

Evaluation of rock mechanical behaviors under uniaxial compression with reference to assessed weathering grades

A weathering classification for granitic rock materials from southeastern Brazil was framed based on core characteristics. The classification was substantiated by a detailed petrographic study. Indirect assessment of weathering grades by density, ultrasonic and Schmidt hammer index tests was performed. Rebound values due to Schmidt hammer multiple impacts at one representative point were more efficient in predicting weathering grades than averaged single impact rebound values, P-wave velocities and densities. Uniaxial compression tests revealed that a large range of uniaxial compressive strength (214-153 MPa) exists in Grade I category where weathering does not seem to have played any role. It was concluded that variability in occurrences of quartz intragranular cracks and in biotite percentage, distribution and orientation might have played a key role in accelerating or decelerating the failure processes of the Grade I specimens. Deterioration of uniaxial compressive strength and elastic modulus and increase in Poisson`s ratio with increasing weathering intensity could be attributed to alteration of minerals, disruption of rock skeleton and microcrack augmentation. A crude relation between failure modes and weathering grades also emerged.

Fluidized ASBR treating synthetic wastewater: Effect of recirculation velocity

An investigation has been performed on the effect of liquid phase recirculation velocity and increasing influent concentration on the stability and efficiency of an anaerobic sequencing batch reactor (ASBR) containing granular biomass. The reactor treated 1.3 L synthetic wastewater at 30 degrees C in 6 h cycles. Initially the effect of recirculation velocity was investigated employing velocities of 5, 7 and 10 m/h and influent concentration of 500 mg COD/L. At these velocities, filtered sample organic matter removal efficiencies were 83, 85 and 84%, respectively. A first order kinetic model could also be fitted to the experimental organic matter concentration profiles. The kinetic parameter values of this model were 1.35, 2.36 and 1.00 h(-1) at the recirculation velocities of 5, 7 and 10 m/h, respectively. The recirculation velocity of 7 m/h was found to be the best operating strategy and this value was maintained while the influent concentration was altered in order to verify system efficiency and stability at increasing organic load. Influent concentration of 1000 mg COD/L resulted in filtered sample organic matter removal efficiency of 80%, and a first order kinetic parameter value of 1.14 h(-1), whereas the concentration of 1500 mg COD/L resulted in an efficiency of 82% and a kinetic parameter value of 1.31 h(-1). (C) 2007 Elsevier B.V. All rights reserved.

Fatigue behaviour of friction stir welded AA2024-T3 alloy: longitudinal and transverse crack growth

The fatigue crack growth properties of friction stir welded joints of 2024-T3 aluminium alloy have been studied under constant load amplitude (increasing-Delta K), with special emphasis on the residual stress (inverse weight function) effects on longitudinal and transverse crack growth rate predictions (Glinka`s method). In general, welded joints were more resistant to longitudinally growing fatigue cracks than the parent material at threshold Delta K values, when beneficial thermal residual stresses decelerated crack growth rate, while the opposite behaviour was observed next to K-C instability, basically due to monotonic fracture modes intercepting fatigue crack growth in weld microstructures. As a result, fatigue crack growth rate (FCGR) predictions were conservative at lower propagation rates and non-conservative for faster cracks. Regarding transverse cracks, intense compressive residual stresses rendered welded plates more fatigue resistant than neat parent plate. However, once the crack tip entered the more brittle weld region substantial acceleration of FCGR occurred due to operative monotonic tensile modes of fracture, leading to non-conservative crack growth rate predictions next to K-C instability. At threshold Delta K values non-conservative predictions values resulted from residual stress relaxation. Improvements on predicted FCGR values were strongly dependent on how the progressive plastic relaxation of the residual stress field was considered.

Thixocasting of an A356 alloy: Fluidity, porosity distribution and thermomechanical fatigue behavior

An extensive set of experiments was performed on a semi-solid A356 alloy in order to assess its flow behavior, mechanical properties, microstructural evolution and porosity level. Three different microstructural conditioning techniques (raw material preparation) were employed: deformation recrystallization, magnetohydrodynamic stirring and low temperature pouring. Measurement of microstructural parameters such as Al-alpha particle size, shape factor, contiguity and entrapped liquid showed a relative equivalency among the various conditioning techniques. It was found that the strongest influence on semi-solid slurry fluidity is exerted by the mould temperature. Tensile properties and porosity levels were measured on a demonstration part produced with different slurry ingate velocities. Results showed similar strength levels among all thixocast samples, a strong correlation between elongation and pore volume fraction and porosity levels much lower than the typical figure for permanent mould or die cast Al-Si alloys. Finally, thermomechanical fatigue tests results were much more favorable to the semi-solid material when compared with the conventionally cast alloy, a result attributed to lower porosity, spheroidal shape of the Al-alpha phase, and refined Si eutectic particles. alpha 2007 Elsevier B.V. All rights reserved.

New prediction methods for CO2 evaporation inside tubes: Part I - A two-phase flow pattern map and a flow pattern based phenomenological model for two-phase flow frictional pressure drops

An updated flow pattern map was developed for CO2 on the basis of the previous Cheng-Ribatski-Wojtan-Thome CO2 flow pattern map [1,2] to extend the flow pattern map to a wider range of conditions. A new annular flow to dryout transition (A-D) and a new dryout to mist flow transition (D-M) were proposed here. In addition, a bubbly flow region which generally occurs at high mass velocities and low vapor qualities was added to the updated flow pattern map. The updated flow pattern map is applicable to a much wider range of conditions: tube diameters from 0.6 to 10 mm, mass velocities from 50 to 1500 kg/m(2) s, heat fluxes from 1.8 to 46 kW/m(2) and saturation temperatures from -28 to +25 degrees C (reduced pressures from 0.21 to 0.87). The updated flow pattern map was compared to independent experimental data of flow patterns for CO2 in the literature and it predicts the flow patterns well. Then, a database of CO2 two-phase flow pressure drop results from the literature was set up and the database was compared to the leading empirical pressure drop models: the correlations by Chisholm [3], Friedel [4], Gronnerud [5] and Muller-Steinhagen and Heck [6], a modified Chisholm correlation by Yoon et al. [7] and the flow pattern based model of Moreno Quiben and Thome [8-10]. None of these models was able to predict the CO2 pressure drop data well. Therefore, a new flow pattern based phenomenological model of two-phase flow frictional pressure drop for CO2 was developed by modifying the model of Moreno Quiben and Thome using the updated flow pattern map in this study and it predicts the CO2 pressure drop database quite well overall. (C) 2007 Elsevier Ltd. All rights reserved.

New prediction methods for CO2 evaporation inside tubes: Part II - An updated general flow boiling heat transfer model based on flow patterns

Corresponding to the updated flow pattern map presented in Part I of this study, an updated general flow pattern based flow boiling heat transfer model was developed for CO2 using the Cheng-Ribatski-Wojtan-Thome [L. Cheng, G. Ribatski, L. Wojtan, J.R. Thome, New flow boiling heat transfer model and flow pattern map for carbon dioxide evaporating inside horizontal tubes, Int. J. Heat Mass Transfer 49 (2006) 4082-4094; L. Cheng, G. Ribatski, L. Wojtan, J.R. Thome, Erratum to: ""New flow boiling heat transfer model and flow pattern map for carbon dioxide evaporating inside tubes"" [Heat Mass Transfer 49 (21-22) (2006) 4082-4094], Int. J. Heat Mass Transfer 50 (2007) 391] flow boiling heat transfer model as the starting basis. The flow boiling heat transfer correlation in the dryout region was updated. In addition, a new mist flow heat transfer correlation for CO2 was developed based on the CO2 data and a heat transfer method for bubbly flow was proposed for completeness sake. The updated general flow boiling heat transfer model for CO2 covers all flow regimes and is applicable to a wider range of conditions for horizontal tubes: tube diameters from 0.6 to 10 mm, mass velocities from 50 to 1500 kg/m(2) s, heat fluxes from 1.8 to 46 kW/m(2) and saturation temperatures from -28 to 25 degrees C (reduced pressures from 0.21 to 0.87). The updated general flow boiling heat transfer model was compared to a new experimental database which contains 1124 data points (790 more than that in the previous model [Cheng et al., 2006, 2007]) in this study. Good agreement between the predicted and experimental data was found in general with 71.4% of the entire database and 83.2% of the database without the dryout and mist flow data predicted within +/-30%. However, the predictions for the dryout and mist flow regions were less satisfactory due to the limited number of data points, the higher inaccuracy in such data, scatter in some data sets ranging up to 40%, significant discrepancies from one experimental study to another and the difficulties associated with predicting the inception and completion of dryout around the perimeter of the horizontal tubes. (C) 2007 Elsevier Ltd. All rights reserved.

On the finite element modeling of fatigue crack growth in pressurized cylindrical shells

A methodology for the computational modeling of the fatigue crack growth in pressurized shell structures, based on the finite element method and concepts of Linear Elastic Fracture Mechanics, is presented. This methodology is based on that developed by Potyondy [Potyondy D, Wawrzynek PA, Ingraffea, AR. Discrete crack growth analysis methodology for through crack in pressurized fuselage structures. Int J Numer Methods Eng 1995;38:1633-1644], which consists of using four stress intensity factors, computed from the modified crack integral method, to predict the fatigue propagation life as well as the crack trajectory, which is computed as part of the numerical simulation. Some issues not presented in the study of Potyondy are investigated herein such as the influence of the crack increment size and the number of nodes per element (4 or 9 nodes) on the simulation results by means of a fatigue crack propagation simulation of a Boeing 737 airplane fuselage. The results of this simulation are compared with experimental results and those obtained by Potyondy [1]. (C) 2008 Elsevier Ltd. All rights reserved.

On the wake-induced vibration of tandem circular cylinders: the vortex interaction excitation mechanism

The mechanism of wake-induced vibrations (WIV) of a pair of cylinders in a tandem arrangement is investigated by experiments. A typical WIV response is characterized by a build-up of amplitude persisting to high reduced velocities; this is different from a typical vortex-induced vibration (VIV) response, which occurs in a limited resonance range. We suggest that WIV of the downstream cylinder is excited by the unsteady vortex-structure interactions between the body and the upstream wake. Coherent vortices interfering with the downstream cylinder induce fluctuations in the fluid force that are not synchronized with the motion. A favourable phase lag between the displacement and the fluid force guarantees that a positive energy transfer from the flow to the structure sustains the oscillations. If the unsteady vortices are removed from the wake of the upstream body then WIV will not be excited. An experiment performed in a steady shear flow turned out to be central to the understanding of the origin of the fluid forces acting on the downstream cylinder.

Numerical solution of the two-phase expansion of a metastable flashing liquid jet using the dispersion-controlled dissipative scheme

This paper presents a study of the stationary phenomenon of superheated or metastable liquid jets, flashing into a two-dimensional axisymmetric domain, while in the two-phase region. In general, the phenomenon starts off when a high-pressure, high-temperature liquid jet emerges from a small nozzle or orifice expanding into a low-pressure chamber, below its saturation pressure taken at the injection temperature. As the process evolves, crossing the saturation curve, one observes that the fluid remains in the liquid phase reaching a superheated condition. Then, the liquid undergoes an abrupt phase change by means of an oblique evaporation wave. Across this phase change the superheated liquid becomes a two-phase high-speed mixture in various directions, expanding to supersonic velocities. In order to reach the downstream pressure, the supersonic fluid continues to expand, crossing a complex bow shock wave. The balance equations that govern the phenomenon are mass conservation, momentum conservation, and energy conservation, plus an equation-of-state for the substance. A false-transient model is implemented using the shock capturing scheme: dispersion-controlled dissipative (DCD), which was used to calculate the flow conditions as the steady-state condition is reached. Numerical results with computational code DCD-2D vI have been analyzed. Copyright (C) 2009 John Wiley & Sons, Ltd.

Ultrasound, eddy current and magnetic Barkhausen noise as tools for sigma phase detection on a UNS S31803 duplex stainless steel

Sigma phase is a deleterious one which can be formed in duplex stainless steels during heat treatment or welding. Aiming to accompany this transformation, ferrite and sigma percentage and hardness were measured on samples of a UNS S31803 duplex stainless steel submitted to heat treatment. These results were compared to measurements obtained from ultrasound and eddy current techniques, i.e., velocity and impedance, respectively. Additionally, backscattered signals produced by wave propagation were acquired during ultrasonic inspection as well as magnetic Barkhausen noise during magnetic inspection. Both signal types were processed via a combination of detrended-fluctuation analysis (DFA) and principal component analysis (PCA). The techniques used were proven to be sensitive to changes in samples related to sigma phase formation due to heat treatment. Furthermore, there is an advantage using these methods since they are nondestructive. (C) 2010 Elsevier B.V. All rights reserved.

Discrete approximation to the global spectrum of the tangent operator for flow past a circular cylinder

This paper deals with the calculation of the discrete approximation to the full spectrum for the tangent operator for the stability problem of the symmetric flow past a circular cylinder. It is also concerned with the localization of the Hopf bifurcation in laminar flow past a cylinder, when the stationary solution loses stability and often becomes periodic in time. The main problem is to determine the critical Reynolds number for which a pair of eigenvalues crosses the imaginary axis. We thus present a divergence-free method, based on a decoupling of the vector of velocities in the saddle-point system from the vector of pressures, allowing the computation of eigenvalues, from which we can deduce the fundamental frequency of the time-periodic solution. The calculation showed that stability is lost through a symmetry-breaking Hopf bifurcation and that the critical Reynolds number is in agreement with the value presented in reported computations. (c) 2007 IMACS. Published by Elsevier B.V. All rights reserved.

MEASUREMENT OF ELASTIC PROPERTIES OF MATERIALS BY THE ULTRASONIC THROUGH-TRANSMISSION TECHNIQUE

The elastic mechanical behavior of elastic materials is modeled by a pair of independent constants (Young`s modulus and Poisson`s coefficient). A precise measurement for both constants is necessary in some applications, such as the quality control of mechanical elements and standard materials used for the calibration of some equipment. Ultrasonic techniques have been used because wave velocity depends on the elastic properties of the propagation medium. The ultrasonic test shows better repeatability and accuracy than the tensile and indentation test. In this work, the theoretical and experimental aspects related to the ultrasonic through-transmission technique for the characterization of elastic solids is presented. Furthermore, an amorphous material and some polycrystalline materials were tested. Results have shown an excellent repeatability and numerical errors that are less than 3% in high-purity samples.