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.

Analysis of stiffened plates by the boundary element method

In this paper, a formulation for representation of stiffeners in plane stress by the boundary elements method (BEM) in linear analysis is presented. The strategy is to adopt approximations for the displacements in the central line of the stiffener. With this simplification the Spurious oscillations in the stress along stiffeners with small thickness is prevented. Worked examples are analyzed to show the efficiency of these techniques, especially in the insertion of very narrow sub-regions, in which quasi-singular integrals are calculated, with stiffeners that are much stiffer than the main domain. The results obtained with this formulation are very close to those obtained with other formulations. (C) 2007 Elsevier Ltd. All rights reserved.

A comparison between neodymium doped alumina samples obtained by Pechini and sol-gel methods using thermo-stimulated luminescence and SEM

Neodymium doped and undoped aluminum oxide samples were obtained using two different techniques: Pechini and sol-gel. Fine grained powders were produced using both procedures, which were analyzed using Scanning Electron Microscopy (SEM) and Thermo-Stimulated Luminescence (TSL). Results showed that neodymium ions incorporation is responsible for the creation of two new TSL peaks (125 and 265 degrees C) and, also, for the enhancement of the intrinsic TSL peak at 190 degrees C. An explanation was proposed for these observations. SEM gave the dimensions of the clusters produced by each method, showing that those obtained by Pechini are smaller than the ones produced by sol-gel; it can also explain the higher emission supplied by the first one. (C) 2010 Elsevier B.V. All rights reserved.

Oxidation and abrasive wear of Fe-Si and Fe-Al intermetallic alloys

In the present work, intermetallic alloys Fe-Si and Fe-Al (Fe(3)Si-C-Cr and Fe(3)Al-C), produced by induction melting, were evaluated regarding their oxidation and abrasive resistance. The tests performed were quasi-isothermal oxidation, cyclic oxidation, and dry sand/rubber wheel abrasion. As reference, the ASTM A297-HH grade stainless steel was tested in the same conditions. In the oxidation tests, the Fe-Al based alloy presented the lowest oxidation rate, and the Fe-Si based alloy achieved the best results in the abrasion test, showing better performance than the HH type stainless steel.

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.

A critical reassessment of elastic unloading in sharp instrumented indentation experiments and the extraction of mechanical properties

This work examines the extraction of mechanical properties from instrumented indentation P-h(s) curves via extensive three-dimensional finite element analyses for pyramidal tips in a wide range of solids under frictional and frictionless contact conditions. Since the topography of the imprint changes with the level of pile-up or sink-in, a relationship is identified between correction factor beta in the elastic equation for the unloading indentation stage and the amount of surface deformation effects. It is shown that the presumption of a constant beta significantly affects mechanical property extractions. Consequently, a new best-fit function is found for the correlation between penetration depth ratios h(e)/h(max), h(r)/h(max) and n, circumventing the need for the assumption of a constant value for beta, made in our prior investigation [Acta Mater. 53 (2005) pp. 3545-3561]. Simulations under frictional contact conditions provide sensible boundaries for the influence of friction on both h(e)/h(max) and h(r)/h(max). Friction is essentially found to induce an overestimation in the inferred n. Instrumented indentation experiments are also performed in three archetypal metallic materials exhibiting distinctly different contact responses. Mechanical property extractions are finally demonstrated in each of these materials.

Floquet stability analysis of the flow around an oscillating cylinder

This investigative work is concerned with the flow around a circular cylinder submitted to forced transverse oscillations. The goal is to investigate how the transition to turbulence is initiated in the wake for cases with different Reynolds numbers (Re) and displacement amplitudes (A). For each Re the motion frequency is kept constant, close to the Strouhal number of the flow around a fixed cylinder at the same Re. Stability analysis of two-dimensional periodic flows around a forced-oscillating cylinder is carried out with respect to three-dimensional infinitesimal perturbations. The procedure consists of performing a Floquet type analysis of time-periodic base flows, computed using the spectral/hp element method. With the results of the Floquet calculations, considerations regarding the stability of the system are drawn, and the form of the instability at its onset is obtained. The critical Reynolds number is observed to change with the amplitude of oscillation. With respect to instabilities, unstable modes with the same symmetry as mode A of a fixed cylinder are observed; however, they present different wavelengths. Also, the instabilities observed for the oscillating cylinder are distinctively stronger in the braid shear layers. Other unstable modes similar to mode B are found. Quasi-periodic modes are observed in the 2S wake, and subharmonic mode occurrences are reported in P + S wakes. (C) 2009 Elsevier Ltd. All rights reserved.

Impact on HDPE and PVC plates - Experimental tests and numerical simulations

This paper presents first material tests on HDPE and PVC, and subsequently impact tests on plates made of the same materials. Finally, numerical simulations of the plate impact tests are compared with the experimental results. A rather comprehensive series of mechanical material tests were performed to disclose the behaviour of PVC and HDPE in tension and compression. Quasi-static tests were carried out at three rates in compression and two in tension. Digital image correlation. DIC, was used to measure the in-plane strains, revealing true stress-strain curves and allowing to analyze strain-rate sensitivity and isotropy of Poisson`s ratio. In addition, dynamic compression tests were carried out in a split-Hopkinson pressure bar. Quasi-static and dynamic tests were also performed on clamped plates made of the same PVC and HDPE materials, using an optical technique to measure the full-field out-of-plane deformations. These tests, together with the material data, were used for comparative purposes of a finite element analysis. A reasonable agreement between experimental and numerical results was achieved. (C) 2010 Elsevier Ltd. All rights reserved.

Optimal design of periodic functionally graded composites with prescribed properties

The computational design of a composite where the properties of its constituents change gradually within a unit cell can be successfully achieved by means of a material design method that combines topology optimization with homogenization. This is an iterative numerical method, which leads to changes in the composite material unit cell until desired properties (or performance) are obtained. Such method has been applied to several types of materials in the last few years. In this work, the objective is to extend the material design method to obtain functionally graded material architectures, i.e. materials that are graded at the local level (e.g. microstructural level). Consistent with this goal, a continuum distribution of the design variable inside the finite element domain is considered to represent a fully continuous material variation during the design process. Thus the topology optimization naturally leads to a smoothly graded material system. To illustrate the theoretical and numerical approaches, numerical examples are provided. The homogenization method is verified by considering one-dimensional material gradation profiles for which analytical solutions for the effective elastic properties are available. The verification of the homogenization method is extended to two dimensions considering a trigonometric material gradation, and a material variation with discontinuous derivatives. These are also used as benchmark examples to verify the optimization method for functionally graded material cell design. Finally the influence of material gradation on extreme materials is investigated, which includes materials with near-zero shear modulus, and materials with negative Poisson`s ratio.

Influence of Microstructures on Polarization Resistance and Hydrogen-Induced Cracking of a Micro-Alloyed Steel Submitted to Different Cooling Rates

The effect of different microstructures on the polarization resistance (Rp) and the hydrogen-induced cracking (HIC) of a micro-alloyed steel austenitized and submitted to different cooling rates was studied. Samples 19.1 x 6 x 2 mm, containing the whole thickness of the plate were extracted from a 20 mm plate and heat treated on a quenching dilatometer, were submitted to Rp and HIC corrosion tests. Both Rp and HIC tests followed as close as possible ASTM G59 and NACE standard TM0284-2003, in this case, modified only with regard to the size of the samples. Steel samples transformed from austenite by a slow cooling (cooling rate of 0.5 degrees C.s(-1)) showed higher susceptibility to hydrogen-induced cracking, with large cracks in the middle of the sample propagating along segregation bands, corresponding to the centerline of the plate thickness. For cooling rates of 10 degrees C.s(-1), only small cracks were found in the matrix and micro cracks nucleated at non-metallic inclusions. For higher cooling rates (40 degrees C.s(-1)) very few small cracks were detected, linked to non-metallic inclusions. This result suggests that structures formed by polygonal structures and segregation bands (were cutectoid microconstituents predominate) have higher susceptibility to HIC. Structures predominantly formed by acicular ferrite make it difficult to propagate the cracks among non-oriented and interlaced acicular ferrite crystals. Smaller segregation bands containing eutectoid products also help inhibit cracking and crack propagation; segregation bands can function as pipelines for hydrogen diffusion and offer a path of stress concentration for the propagation of cracks, frequently associated to non-metallic inclusions. Polarization resistance essays performed on the steel in theas received condition, prior to any heat treatment, showed larger differences between the regions of the plate, with a considerably lower Rp in the centerline. The austenitization heat treatments followed by cooling rates of 0.5 e 10 degrees C.s(-1) made more uniform the corrosion resistance along the thickness of the plate. The effects of heat treatments on the corrosion resistance are probably related to the microconstituent formed, allied to the chemical homogenization of the impurities concentrated on the centerline of the plate.

Combining Mager and Steinmetz: The Effect of Grain Size and Maximum Induction on Hysteresis Energy Loss

Twelve samples with different grain sizes were prepared by normal grain growth and by primary recrystallization, and the hysteresis dissipated energy was measured by a quasi-static method. Results showed a linear relation between hysteresis energy loss and the inverse of grain size, which is here called Mager`s law, for maximum inductions from 0.6 to 1.5 T, and a Steinmetz power law relation between hysteresis loss and maximum induction for all samples. The combined effect is better described by a Mager`s law where the coefficients follow Steinmetz law.

Reduction of Chromium from Al(2)O(3)-CaO-SiO(2)-CrOx Slags by Carbon Dissolved in Liquid Iron

The goal of this work is to investigate the reduction of chromium from a quaternary slag by carbon dissolved in liquid steel. Laboratory scale experiments were conducted to study the reduction of chromium oxides in the slag by carbon dissolved in the melt. These experiments were made under different conditions of slag basicity and amount of added carbon. Thermodynamic calculations based on Double Sublattice model were applied using the commercial software Thermo-Calc, with the IRSID database. The results obtained showed good correlation with practical and calculated results, making it possible to predict equilibrium conditions of the system and to determine the activities of chromium oxides in the slag.