A STUDY OF PENALTY FORMULATIONS USED IN THE NUMERICAL APPROXIMATION OF A RADIALLY SYMMETRIC ELASTICITY PROBLEM

We consider a class of two-dimensional problems in classical linear elasticity for which material overlapping occurs in the absence of singularities. Of course, material overlapping is not physically realistic, and one possible way to prevent it uses a constrained minimization theory. In this theory, a minimization problem consists of minimizing the total potential energy of a linear elastic body subject to the constraint that the deformation field must be locally invertible. Here, we use an interior and an exterior penalty formulation of the minimization problem together with both a standard finite element method and classical nonlinear programming techniques to compute the minimizers. We compare both formulations by solving a plane problem numerically in the context of the constrained minimization theory. The problem has a closed-form solution, which is used to validate the numerical results. This solution is regular everywhere, including the boundary. In particular, we show numerical results which indicate that, for a fixed finite element mesh, the sequences of numerical solutions obtained with both the interior and the exterior penalty formulations converge to the same limit function as the penalization is enforced. This limit function yields an approximate deformation field to the plane problem that is locally invertible at all points in the domain. As the mesh is refined, this field converges to the exact solution of the plane problem.

Verification of a mixed high-order accurate DNS code for laminar turbulent transition by the method of manufactured solutions

This paper presents results on a verification test of a Direct Numerical Simulation code of mixed high-order of accuracy using the method of manufactured solutions (MMS). This test is based on the formulation of an analytical solution for the Navier-Stokes equations modified by the addition of a source term. The present numerical code was aimed at simulating the temporal evolution of instability waves in a plane Poiseuille flow. The governing equations were solved in a vorticity-velocity formulation for a two-dimensional incompressible flow. The code employed two different numerical schemes. One used mixed high-order compact and non-compact finite-differences from fourth-order to sixth-order of accuracy. The other scheme used spectral methods instead of finite-difference methods for the streamwise direction, which was periodic. In the present test, particular attention was paid to the boundary conditions of the physical problem of interest. Indeed, the verification procedure using MMS can be more demanding than the often used comparison with Linear Stability Theory. That is particularly because in the latter test no attention is paid to the nonlinear terms. For the present verification test, it was possible to manufacture an analytical solution that reproduced some aspects of an instability wave in a nonlinear stage. Although the results of the verification by MMS for this mixed-order numerical scheme had to be interpreted with care, the test was very useful as it gave confidence that the code was free of programming errors. Copyright (C) 2009 John Wiley & Sons, Ltd.

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.

DESIGN OF CONCRETE PLATFORMS OFFSHORE USING THE THREE-LAYER SHELL THEORY

The concrete offshore platforms, which are subjected a several loading combinations and, thus, requires an analysis more generic possible, can be designed using the concepts adopted to shell elements, but the resistance must be verify in particular cross-sections to shear forces. This work about design of shell elements will be make using the three-layer shell theory. The elements are subject to combined loading of membrane and plate, totalizing eight components of internal forces, which are three membrane forces, three moments (two out-of-plane bending moments and one in-plane, or torsion, moment) and two shear forces. The design method adopted, utilizing the iterative process proposed by Lourenco & Figueiras (1993) obtained from equations of equilibrium developed by Gupta (1896) , will be compared to results of experimentally tested shell elements found in the literature using the program DIANA.

Shell curvature as an initial deformation: A geometrically exact finite element approach

An alternative approach for the analysis of arbitrarily curved shells is developed in this paper based on the idea of initial deformations. By `alternative` we mean that neither differential geometry nor the concept of degeneration is invoked here to describe the shell surface. We begin with a flat reference configuration for the shell mid-surface, after which the initial (curved) geometry is mapped as a stress-free deformation from the plane position. The actual motion of the shell takes place only after this initial mapping. In contrast to classical works in the literature, this strategy enables the use of only orthogonal frames within the theory and therefore objects such as Christoffel symbols, the second fundamental form or three-dimensional degenerated solids do not enter the formulation. Furthermore, the issue of physical components of tensors does not appear. Another important aspect (but not exclusive of our scheme) is the possibility to describe exactly the initial geometry. The model is kinematically exact, encompasses finite strains in a totally consistent manner and is here discretized under the light of the finite element method (although implementation via mesh-free techniques is also possible). Assessment is made by means of several numerical simulations. Copyright (C) 2009 John Wiley & Sons, Ltd.

The natural force density method for the shape finding of taut structures

This work presents, with the aid of the natural approach, an extension of the force density method for the initial shape finding of cable and membrane structures, which leads to the solution of a system of linear equations. This method, here called the natural force density method, preserves the linearity which characterizes the original force density method. At the same time, it overcomes the difficulties that the original procedure presents to cope with irregular triangular finite element meshes. Furthermore, if this method is applied iteratively in the lines prescribed herewith, it leads to a viable initial configuration with a uniform, isotropic plane Cauchy stress state. This means that a minimal surface for the membrane can be achieved through a succession of equilibrated configurations. Several numerical examples illustrate the simplicity and robustness of the method. (C) 2008 Elsevier B.V. All rights reserved.

Experimental investigation of flow-induced vibration on isolated and tandem circular cylinders fitted with strakes

The effect of varying the geometric parameters of helical strakes on vortex-induced vibration (VIV) is investigated in this paper. The degree of oscillation attenuation or even suppression is analysed for isolated circular cylinder cases. How a cylinder fitted with strakes behaves when immersed in the wake of another cylinder in tandem arrangement is also investigated and these results are compared to those with a single straked cylinder. The experimental tests are conducted at a circulating water channel facility and the cylindrical models are mounted on a low-damping air bearing elastic base with one degree-of-freedom, restricted to oscillate in the transverse direction to the channel flow. Three strake pitches (p) and heights (h) are tested: p = 5, 10, 15d, and h = 0.1, 0.2, 0.25d. The mass ratio is 1.8 for all models. The Reynolds number range is from 1000 to 10000, and the reduced velocity varies up to 21. The cases with h = 0.1d strakes reduce the amplitude response when compared to the isolated plain cylinder, however the oscillation still persists. On the other hand, the cases with h = 0.2, 0.25d strakes almost completely suppress VIV. Spanwise vorticity fields, obtained through stereoscopic digital particle image velocimetry (SDPIV), show an alternating vortex wake for the p = 10d and h = 0.1d straked cylinder. The p = 10d and h = 0.2d cylinder wake has separated shear layers with constant width and no roll-up close to the body. The strakes do not increase the magnitude of the out-of-plane velocity compared to the isolated plain cylinder. However, they deflect the flow in the out-of-plane direction in a controlled way, which can prevent the vortex shedding correlation along the span. In order to investigate the wake interference effect on the strake efficiency, an experimental arrangement with two cylinders in tandem is employed. The centre-to-centre distance for the tandem arrangement varies from 2 to 6. When the downstream p = 10d and h = 0.2d cylinder is immersed in the wake of an upstream fixed plain cylinder, it loses its effectiveness compared with the isolated case. Although the oscillations have significant amplitude, they are limited, which is a different behaviour from that of a tandem configuration with two plain cylinders. For this particular case, the amplitude response monotonically increases for all gaps, except one, a trait usually found in galloping-like oscillations. SDPIV results for the tandem arrangements show alternating vortex shedding and oscillatory wake. (C) 2010 Elsevier Ltd. All rights reserved.

Matrix Method for Acoustic Levitation Simulation

A matrix method is presented for simulating acoustic levitators. A typical acoustic levitator consists of an ultrasonic transducer and a reflector. The matrix method is used to determine the potential for acoustic radiation force that acts on a small sphere in the standing wave field produced by the levitator. The method is based on the Rayleigh integral and it takes into account the multiple reflections that occur between the transducer and the reflector. The potential for acoustic radiation force obtained by the matrix method is validated by comparing the matrix method results with those obtained by the finite element method when using an axisymmetric model of a single-axis acoustic levitator. After validation, the method is applied in the simulation of a noncontact manipulation system consisting of two 37.9-kHz Langevin-type transducers and a plane reflector. The manipulation system allows control of the horizontal position of a small levitated sphere from -6 mm to 6 mm, which is done by changing the phase difference between the two transducers. The horizontal position of the sphere predicted by the matrix method agrees with the horizontal positions measured experimentally with a charge-coupled device camera. The main advantage of the matrix method is that it allows simulation of non-symmetric acoustic levitators without requiring much computational effort.

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.

Computational Method to Determine Reflected Ultrasonic Signals From Arbitrary-Geometry Targets

A computational method based on the impulse response and on the discrete representation computational concept is proposed for the determination of the echo responses from arbitrary-geometry targets. It is supposed that each point of the transducer aperture can be considered as a source radiating hemispherical waves to the reflector. The local interaction with each of the hemispherical waves at the reflector surface can be modeled as a plane wave impinging on a planar surface, using the respective reflection coefficient. The method is valid for all field regions and can be performed for any excitation waveform radiated from an arbitrary acoustic aperture. The effects of target geometry, position, and material on both the amplitude and the shape of the echo response are studied. The model is compared with experimental results obtained using broadband transducers together with plane and cylindrical concave rectangular reflectors (aluminum, brass, and acrylic), as well as a circular cavity placed on a plane surface, in a water medium. The method can predict the measured echoes accurately. This paper shows an improved approach of the method, considering the reflection coefficient for all incident hemispherical waves arriving at each point of the target surface.

Finite Element Analysis and Optimization of a Single-Axis Acoustic Levitator

A finite element analysis and a parametric optimization of single-axis acoustic levitators are presented. The finite element method is used to simulate a levitator consisting of a Langevin ultrasonic transducer with a plane radiating surface and a plane reflector. The transducer electrical impedance, the transducer face displacement, and the acoustic radiation potential that acts on small spheres are determined by the finite element method. The numerical electrical impedance is compared with that acquired experimentally by an impedance analyzer, and the predicted displacement is compared with that obtained by a fiber-optic vibration sensor. The numerical acoustic radiation potential is verified experimentally by placing small spheres in the levitator. The same procedure is used to optimize a levitator consisting of a curved reflector and a concave-faced transducer. The numerical results show that the acoustic radiation force in the new levitator is enhanced 604 times compared with the levitator consisting of a plane transducer and a plane reflector. The optimized levitator is able to levitate 3, 2.5-mm diameter steel spheres with a power consumption of only 0.9 W.

Modeling of functionally graded piezoelectric ultrasonic transducers

The application of functionally graded material (FGM) concept to piezoelectric transducers allows the design of composite transducers without interfaces, due to the continuous change of property values. Thus, large improvements can be achieved, as reduction of stress concentration, increasing of bonding strength, and bandwidth. This work proposes to design and to model FGM piezoelectric transducers and to compare their performance with non-FGM ones. Analytical and finite element (FE) modeling of FGM piezoelectric transducers radiating a plane pressure wave in fluid medium are developed and their results are compared. The ANSYS software is used for the FE modeling. The analytical model is based on FGM-equivalent acoustic transmission-line model, which is implemented using MATLAB software. Two cases are considered: (i) the transducer emits a pressure wave in water and it is composed of a graded piezoceramic disk, and backing and matching layers made of homogeneous materials; (ii) the transducer has no backing and matching layer; in this case, no external load is simulated. Time and frequency pressure responses are obtained through a transient analysis. The material properties are graded along thickness direction. Linear and exponential gradation functions are implemented to illustrate the influence of gradation on the transducer pressure response, electrical impedance, and resonance frequencies. (C) 2009 Elsevier B. V. All rights reserved.

Methods of Fourier optics in digital holographic microscopy

In this paper, processing methods of Fourier optics implemented in a digital holographic microscopy system are presented. The proposed methodology is based on the possibility of the digital holography in carrying out the whole reconstruction of the recorded wave front and consequently, the determination of the phase and intensity distribution in any arbitrary plane located between the object and the recording plane. In this way, in digital holographic microscopy the field produced by the objective lens can be reconstructed along its propagation, allowing the reconstruction of the back focal plane of the lens, so that the complex amplitudes of the Fraunhofer diffraction, or equivalently the Fourier transform, of the light distribution across the object can be known. The manipulation of Fourier transform plane makes possible the design of digital methods of optical processing and image analysis. The proposed method has a great practical utility and represents a powerful tool in image analysis and data processing. The theoretical aspects of the method are presented, and its validity has been demonstrated using computer generated holograms and images simulations of microscopic objects. (c) 2007 Elsevier B.V. All rights reserved.

Effect of alpha prime due to 475 degrees C aging on fracture behavior and corrosion resistance of DIN 1.4575 and MA 956 high performance ferritic stainless steels

The 475 degrees C embrittlement in stainless steels is a well-known phenomenon associated to alpha prime (alpha`) formed by precipitation or spinodal decomposition. Many doubts still remain on the mechanism of alpha` formation and its consequence on deformation and fracture mechanisms and corrosion resistance. In this investigation, the fracture behavior and corrosion resistance of two high performance ferritic stainless steels were investigated: a superferritic DIN 1.4575 and MA 956 superalloy were evaluated. Samples of both stainless steels (SS) were aged at 475 degrees C for periods varying from 1 to 1,080 h. Their fracture surfaces were observed using scanning electron microscopy (SEM) and the cleavage planes were determined by electron backscattering diffraction (EBSD). Some samples were tested for corrosion resistance using electrochemical impedance spectroscopy (EIS) and potentiodynamic polarization. Brittle and ductile fractures were observed in both ferritic stainless steels after aging at 475 degrees C. For aging periods longer than 500 h, the ductile fracture regions completely disappeared. The cleavage plane in the DIN 1.4575 samples aged at 475 degrees C for 1,080 h was mainly {110}, however the {102}, {314}, and {131} families of planes were also detected. The pitting corrosion resistance decreased with aging at 475 degrees C. The effect of alpha prime on the corrosion resistance was more significant in the DIN 1.4575 SS comparatively to the Incoloy MA 956.

Three-dimensional analysis of fracture, corrosion and wear surfaces

A brief look at the history of fractography has shown a recent trend in the quantification of topographic parameters through the use of three-dimensional reconstruction techniques, which associate SEM stereoscopy and stereophotogrammetry software, allowing the calculation of the elevation measurement at numerous points of the topography due to the parallax that takes place during the tilting of the sample along the microscope eucentric plane. Several investigators have used reconstruction techniques to correlate some fractographic parameters, such as fractal dimension and fractured to projected area ratio, to the mechanical properties of materials, such as fracture toughness and tensile strength. So far, the search for a clear relationship between the fracture topography and mechanical properties has provided ambiguous results. The present work applied a surface metrology software to reconstruct three-dimensionally fracture surfaces (transgranular cleavage, intergranular and dimple fracture), corrosion pits and tribo-surfaces in order to explore the potential of this stereophotogrammetry technique. The existence of a variation in the calculated topographic parameters with the conditions of SEM image acquisition reinforces the importance of both good image acquisition and accurate calibration methods in order to validate this 3D reconstruction technique in metrological terms. Preliminary results did not indicate the existence of a clear relationship between either the true to project area ratio and CVN absorbed energy or the fractal dimension and CVN absorbed energy. It is likely that each fracture mechanism presents a proper relationship between the fractographic parameters and mechanical properties. (C) 2009 Elsevier Ltd. All rights reserved.