Finite element computational modeling of externally bonded CFRP composites flexural behavior in RC beams

This paper focuses on the flexural behavior of RC beams externally strengthened with Carbon Fiber Reinforced Polymers (CFRP) fabric. A non-linear finite element (FE) analysis strategy is proposed to support the beam flexural behavior experimental analysis. A development system (QUEBRA2D/FEMOOP programs) has been used to accomplish the numerical simulation. Appropriate constitutive models for concrete, rebars, CFRP and bond-slip interfaces have been implemented and adjusted to represent the composite system behavior. Interface and truss finite elements have been implemented (discrete and embedded approaches) for the numerical representation of rebars, interfaces and composites.

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.

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.

Inertia in two simple problems of soil hydrodynamics

In the present paper the dynamic solutions of two non-steady seepage problems are discussed. It is shown that the acceleration term in the equation of motion is important for a correct qualitative description of the flow.

An exact conserving algorithm for nonlinear dynamics with rotational DOFs and general hyperelasticity. Part 2: shells

Following the approach developed for rods in Part 1 of this paper (Pimenta et al. in Comput. Mech. 42:715-732, 2008), this work presents a fully conserving algorithm for the integration of the equations of motion in nonlinear shell dynamics. We begin with a re-parameterization of the rotation field in terms of the so-called Rodrigues rotation vector, allowing for an extremely simple update of the rotational variables within the scheme. The weak form is constructed via non-orthogonal projection, the time-collocation of which ensures exact conservation of momentum and total energy in the absence of external forces. Appealing is the fact that general hyperelastic materials (and not only materials with quadratic potentials) are permitted in a totally consistent way. Spatial discretization is performed using the finite element method and the robust performance of the scheme is demonstrated by means of numerical examples.

An exact conserving algorithm for nonlinear dynamics with rotational DOFs and general hyperelasticity. Part 1: Rods

A fully conserving algorithm is developed in this paper for the integration of the equations of motion in nonlinear rod dynamics. The starting point is a re-parameterization of the rotation field in terms of the so-called Rodrigues rotation vector, which results in an extremely simple update of the rotational variables. The weak form is constructed with a non-orthogonal projection corresponding to the application of the virtual power theorem. Together with an appropriate time-collocation, it ensures exact conservation of momentum and total energy in the absence of external forces. Appealing is the fact that nonlinear hyperelastic materials (and not only materials with quadratic potentials) are permitted without any prejudice on the conservation properties. Spatial discretization is performed via the finite element method and the performance of the scheme is assessed by means of several numerical simulations.

Finite element analysis of concrete cracking at early age

The study of the early age concrete properties is becoming more important, as the thermal effects and the shrinkage, even in the first hours, could generate cracks, increasing the permeability of the structure and being able to induce problems of durability and functionality in the same ones. The detailed study of the stresses development during the construction process can be decisive to keep low the cracking levels. In this work a computational model, based on the finite element method, was implemented to simulate the early age concrete behavior and, specially, the evaluation of the cracking risk. The finite element analysis encloses the computational modeling of the following phenomena: chemical, thermal, moisture diffusion and mechanical which occur at the first days after the concrete cast. The developed software results were compared with experimental values found in the literature, demonstrating an excellent approach for all the implemented analysis.

Finite element and dimensional analysis algorithm for the prediction of mechanical properties of bulk materials and thin films

In this work, the applicability of a new algorithm for the estimation of mechanical properties from instrumented indentation data was studied for thin films. The applicability was analyzed with the aid of both three-dimensional finite element simulations and experimental indentation tests. The numerical approach allowed studying the effect of the substrate on the estimation of mechanical properties of the film, which was conducted based on the ratio h(max)/l between maximum indentation depth and film thickness. For the experimental analysis, indentation tests were conducted on AISI H13 tool steel specimens, plasma nitrated and coated with TiN thin films. Results have indicated that, for the conditions analyzed in this work, the elastic deformation of the substrate limited the extraction of mechanical properties of the film/substrate system. This limitation occurred even at low h(max)/l ratios and especially for the estimation of the values of yield strength and strain hardening exponent. At indentation depths lower than 4% of the film thickness, the proposed algorithm estimated the mechanical properties of the film with accuracy. Particularly for hardness, precise values were estimated at h(max)/l lower than 0.1, i.e. 10% of film thickness. (C) 2010 Published by Elsevier B.V.

Galerkin finite element method for incompressible thermofluid flows framed within the bond graph theory

In this paper a bond graph methodology is used to model incompressible fluid flows with viscous and thermal effects. The distinctive characteristic of these flows is the role of pressure, which does not behave as a state variable but as a function that must act in such a way that the resulting velocity field has divergence zero. Velocity and entropy per unit volume are used as independent variables for a single-phase, single-component flow. Time-dependent nodal values and interpolation functions are introduced to represent the flow field, from which nodal vectors of velocity and entropy are defined as state variables. The system for momentum and continuity equations is coincident with the one obtained by using the Galerkin method for the weak formulation of the problem in finite elements. The integral incompressibility constraint is derived based on the integral conservation of mechanical energy. The weak formulation for thermal energy equation is modeled with true bond graph elements in terms of nodal vectors of temperature and entropy rates, resulting a Petrov-Galerkin method. The resulting bond graph shows the coupling between mechanical and thermal energy domains through the viscous dissipation term. All kind of boundary conditions are handled consistently and can be represented as generalized effort or flow sources. A procedure for causality assignment is derived for the resulting graph, satisfying the Second principle of Thermodynamics. (C) 2007 Elsevier B.V. All rights reserved.

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.

A simple and effective inverse projection scheme for void distribution control in topology optimization

The ability to control both the minimum size of holes and the minimum size of structural members are essential requirements in the topology optimization design process for manufacturing. This paper addresses both requirements by means of a unified approach involving mesh-independent projection techniques. An inverse projection is developed to control the minimum hole size while a standard direct projection scheme is used to control the minimum length of structural members. In addition, a heuristic scheme combining both contrasting requirements simultaneously is discussed. Two topology optimization implementations are contributed: one in which the projection (either inverse or direct) is used at each iteration; and the other in which a two-phase scheme is explored. In the first phase, the compliance minimization is carried out without any projection until convergence. In the second phase, the chosen projection scheme is applied iteratively until a solution is obtained while satisfying either the minimum member size or minimum hole size. Examples demonstrate the various features of the projection-based techniques presented.

Medium carbon steel deep drawing: A study on the evolution of mechanical properties, texture and simulations, from cold rolling to the end product

Medium carbon steels are mostly used for simple applications; nevertheless new applications have been developed for which good sheet formability is required. This class of steels has an inherent low formability. A medium carbon hot rolled SAE 1050 steel has been selected for this study. It has been cold rolled with reductions in the 7-80% range. Samples have been used to assess the cold work hardening curve. For samples with a 50 and 80% thickness reduction, an annealing heat treatment has been performed to obtain recrystallization. The material has been characterized in the ""as received"", cold rolled and annealed conditions, using several methods: optical microscopy, X-ray diffraction (texture), Vickers hardness and tensile testing. The 50% cold rolled and recrystallized material has been further studied in terms of sheet metal formability and texture evolution during the actual stamping of a steel toecap that has been used to validate the finite element simulations. (C) 2008 Elsevier B.V. All rights reserved.

Coarse-grid higher order finite-difference time-domain algorithm with low dispersion errors

Higher order (2,4) FDTD schemes used for numerical solutions of Maxwell`s equations are focused on diminishing the truncation errors caused by the Taylor series expansion of the spatial derivatives. These schemes use a larger computational stencil, which generally makes use of the two constant coefficients, C-1 and C-2, for the four-point central-difference operators. In this paper we propose a novel way to diminish these truncation errors, in order to obtain more accurate numerical solutions of Maxwell`s equations. For such purpose, we present a method to individually optimize the pair of coefficients, C-1 and C-2, based on any desired grid size resolution and size of time step. Particularly, we are interested in using coarser grid discretizations to be able to simulate electrically large domains. The results of our optimization algorithm show a significant reduction in dispersion error and numerical anisotropy for all modeled grid size resolutions. Numerical simulations of free-space propagation verifies the very promising theoretical results. The model is also shown to perform well in more complex, realistic scenarios.

Universal Set-top Box: A Simple Design to Provide Accessible Services

This paper presents the design of a low cost accessible digital television set-top box. This set-top box was designed and tested to the International ISDB-T system and considered the adoption of solutions that would provide accessible services in digital television in the simplest digital television receiver. The accessible set-top box was evaluated regarding the processing and memory requirements impacts to provide the features for accessible services. The work presents also the access services bandwidth consumption analysis(1).

Simple Excitation Model for Coaxial Driven Monopole Antennas

A new excitation model for the numerical solution of field integral equation (EFIE) applied to arbitrarily shaped monopole antennas fed by coaxial lines is presented. This model yields a stable solution for the input impedance of such antennas with very low numerical complexity and without the convergence and high parasitic capacitance problems associated with the usual delta gap excitation.