Use of the energy-release rate in a 2 scale damage analysis of a bonded joint

Stress singularities appear at the extremities of an adhesive bond. They can produce a damage mechanism that we assimilate in this Note to a crack. The energy release rate permits to characterize its evolution. But a very refined mesh would be necessary for a real structure. Using an asymptotic method based on the small thickness of the bond a limit model with a different local behaviour is suggested. It leads to an approximation of the energy release rate

Synchrotron strain scanning for residual stress measurement in cold-drawn steel rods

Cold-drawn steel rods and wires retain significant residual stresses as a consequence of the manufacturing process. These residual stresses are known to be detrimental for the mechanical properties of the wires and their durability in aggressive environments. Steel makers are aware of the problem and have developed post-drawing processes to try and reduce the residual stresses on the wires. The present authors have studied this problem for a number of years and have performed a detailed characterization of the residual stress state inside cold-drawn rods, including both experimental and numerical techniques. High-energy synchrotron sources have been particularly useful for this research. The results have shown how residual stresses evolve as a consequence of cold-drawing and how they change with subsequent post-drawing treatments. The authors have been able to measure for the first time a complete residual strain profile along the diameter in both phases (ferrite and cementite) of a cold-drawn steel rod.

Residual stress distribution in ferritic-austenitic steel joints made by laser welding

The present investigation addresse the influence of laser welding process-ing parameters used for joining dis-similar metals (ferritic to austenitic steel), on the induced residual stress field. Welding was performed on a Nd:YAG laser DY033 (3300 W) in a continuous wave (CW), keyhole mode. The base metals (BM) employed in this study are AISI 1010 carbon steel (CS) and AISI 304L austenitic stainless steel (SS). Pairs of dissimilar plates of 200 mm x 45 mm x 3 mm were butt joined by laser welding. Different sets of parameters were used to engineer the base metals apportionment at joint formation, namely distinct dilution rates. Residual strain scanning, carried out by neutron diffraction was used to assess the joints. Through-thickness residual stress maps were determined for the laser welded samples of dis-similar steels using high spatial reso-lution. As a result, an appropriate set of processing parameters, able to mi-nimize the local tensile residual stress associated to the welding process, was found.

Residual stress distribution in ferritic to austenitic steel joints made by laser welding,

In this study, autogenous laser welding was used to join thin plates of low carbon ferritic and austenitic stainless steel. Due to the differences in the thermo-physical properties of base metals, this kind of weld exhibits a complex microstructure, which frequently leads to an overall loss of joint quality. Four welded samples were prepared by using different sets of processing parameters, with the aim of minimizing the induced residual stress field. The dissimilar austenitic-ferritic joints obtained under all welding conditions were uniform and free of defects. Variations in beam position did not influence the weld geometiy, which is a typical keyhole welding. Microstructural characterization and residual strain scanning (by neutron diffraction) were used to assess the features of the joints. By varying laser beam power density and by displacing the laser beam towards the carbon steel side, an optimum combination of processing parameters was found.

Mechanical behavior and residual stress field of dissimilar laser welded joints

The present investigation addresses the mechanical behavior and residual stress field of dissimilar joints produced by laser welding. Microstructure characterization and residual strain scanning, carried out by neutron diffraction, were used to assess the joints features. It was found that the heat source position influences the base metals dilution and the residual stress field associated to the welding process. The tensile behavior of the joint, different zones achieved by using a video-image based system (VIC-2D) reveals that the residual stress field, together with the positive difference in yield between the weld metal and the base materials protects the joint from being plastically deformed.

Induction of through-thickness compressive residual stress fields in thin Al2024-T351 plates by laser shock processing

Laser Shock Processing (LSP) has been demonstrated as an emerging technique for the induction of RS’s fields in subsurface layers of relatively thick specimens. However, the LSP treatment of relatively thin specimens brings, as an additional consequence, the possible bending in a process of laser shock forming. This effect poses a new class of problems regarding the attainment of specified RS’s depth profiles in the mentioned type of sheets, and, what can be more critical, an overall deformation of the treated component. The analysis of the problem of LSP treatment for induction of tentatively through-thickness RS’s fields for fatigue life enhancement in relatively thin sheets in a way compatible with reduced overall workpiece deformation due to spring-back self-equilibration is envisaged in this paper. The coupled theoretical-experimental predictive approach developed by the authors has been applied to the specification of LSP treatments for achievement of RS's fields tentatively able to retard crack propagation on normalized specimens. A convergence between numerical code results and experimental results coming from direct RS's measurement is presented as a first step for the treatment of the normalized specimens under optimized conditions and verification of the crack retardation properties virtually induced.

Numerical simulation of the upward propagation of a flame in a vertical tube filled with a very lean mixture.

Upwardpropagation of a premixed flame in averticaltubefilled with a very leanmixture is simulated numerically using a single irreversible Arrhenius reaction model with infinitely high activation energy. In the absence of heat losses and preferential diffusion effects, a curved flame with stationary shape and velocity close to those of an open bubble ascending in the same tube is found for values of the fuel mass fraction above a certain minimum that increases with the radius of the tube, while the numerical computations cease to converge to a stationary solution below this minimum mass fraction. The vortical flow of the gas behind the flame and in its transport region is described for tubes of different radii. It is argued that this flow may become unstable when the fuel mass fraction is decreased, and that this instability, together with the flame stretch due to the strong curvature of the flame tip in narrow tubes, may be responsible for the minimum fuel mass fraction. Radiation losses and a Lewis number of the fuel slightly above unity decrease the final combustion temperature at the flame tip and increase the minimum fuel mass fraction, while a Lewis number slightly below unity has the opposite effect.

Numerical simulation of the evolution of viscous rotating drops

In this contribution we simulate numerically the evolution of a viscous fluid drop rotating about a fixed axis at constant angular velocity ? or constant angular momentum L, surrounded by another viscous fluid. The problem is considered in the limit of large Ekman number and small Reynolds number. In the lecture we will describe the numerical method we have used to solve the PDE system that describes the evolution of the drop (3D boundary element method). We will also present the results we have obtained, paying special attention to the stability/instability of the equilibrium shapes.

Assessment of Narrow-Body Transport Airplane Evacuation by Numerical Simulation

This paper presents the results obtained with a new agent-based computer model that can simulate the evacuation of narrow-body transport airplanes in the conditions prescribed by the airworthiness regulations for certification. The model, described in detail in a former paper, has been verified with real data of narrow-body certification demonstrations. Numerical simulations of around 20 narrow-body aircraft, representative of current designs in various market segments, show the capabilities of the model and provide relevant information on the relationship between cabin features and emergency evacuation. The longitudinal location of emergency exits seems to be even more important than their size or the overall margin with respect to the prescribed number and type of exits indicated by the airworthiness requirements

Numerical simulation of tangling in jet engine turbines

The numerical analysis of certain safety related problems presents serious difficulties, since the large number of components present leads to huge finite elementmodels that can only be solved by using large and expensive computers or by making rough approaches to the problem. Tangling, or clashing, in the turbine of a jet engine airplane is an example of such problems. This is caused by the crash and friction between rotor and stator blades in the turbine after an eventual shaft failure. When facing the study of an event through numerical modelling, the accurate simulation of this problem would require the engineer to model all the rotor and stator blades existing in the turbine stage, using a small element size in all pieces. Given that the number of stator and rotor blades is usually around 200, such simulations would require millions of elements. This work presents a new numerical methodology, specifically developed for the accurate modelling of the tangling problem that, depending on the turbine configuration, is able to reduce the number of nodes up to an order of magnitude without losing accuracy. The methodology, which benefits from the cyclic configuration of turbines, is successfully applied to the numerical analysis of a hypothetical tangling event in a turbine, providing valuable data such as the rotating velocity decrease of the turbine, the braking torque and the damage suffered by the blades. The methodology is somewhat general and can be applied to any problem in which damage caused by the interaction between a rotating and static piece is to be analysed.

Numerical simulation of Ge solar cells using D-AMPS-1D code

A solar cell is a solid state device that converts the energy of sunlight directly into electricity by the photovoltaic effect. When light with photon energies greater than the band gap is absorbed by a semiconductor material, free electrons and free holes are generated by optical excitation in the material. The main characteristic of a photovoltaic device is the presence of internal electric field able to separate the free electrons and holes so they can pass out of the material to the external circuit before they recombine. Numerical simulation of photovoltaic devices plays a crucial role in their design, performance prediction, and comprehension of the fundamental phenomena ruling their operation. The electrical transport and the optical behavior of the solar cells discussed in this work were studied with the simulation code D-AMPS-1D. This software is an updated version of the one-dimensional (1D) simulation program Analysis of Microelectronic and Photonic Devices (AMPS) that was initially developed at The Penn State University, USA. Structures such as homojunctions, heterojunctions, multijunctions, etc., resulting from stacking layers of different materials can be studied by appropriately selecting characteristic parameters. In this work, examples of cells simulation made with D-AMPS-1D are shown. Particularly, results of Ge photovoltaic devices are presented. The role of the InGaP buffer on the device was studied. Moreover, a comparison of the simulated electrical parameters with experimental results was performed.

Numerical prediction of residual stresses fields induced in thin Al2024-T351 sheets by laser shock processing

Outline: • Introduction • Numerical model SHOCKLAS© • Single LSP pulses • Overlapped LSP pulses • Discussion and Outlook

Numerical simulation of the transition from densely distributed to localized cracking.

8th International Conference on Fracture Mechanics of Concrete and Concrete Structures (FraMCoS8).

Numerical simulation of group combustion of pulverized coal

A mathematical model for the group combustion of pulverized coal particles was developed in a previous work. It includes the Lagrangian description of the dehumidification, devolatilization and char gasification reactions of the coal particles in the homogenized gaseous environment resulting from the three fuels, CO, H2 and volatiles, supplied by the gasification of the particles and their simultaneous group combustion by the gas phase oxidation reactions, which are considered to be very fast. This model is complemented here with an analysis of the particle dynamics, determined principally by the effects of aerodynamic drag and gravity, and its dispersion based on a stochastic model. It is also extended to include two other simpler models for the gasification of the particles: the first one for particles small enough to extinguish the surrounding diffusion flames, and a second one for particles with small ash content when the porous shell of ashes remaining after gasification of the char, non structurally stable, is disrupted. As an example of the applicability of the models, they are used in the numerical simulation of an experiment of a non-swirling pulverized coal jet with a nearly stagnant air at ambient temperature, with an initial region of interaction with a small annular methane flame. Computational algorithms for solving the different stages undergone by a coal particle during its combustion are proposed. For the partial differential equations modeling the gas phase, a second order finite element method combined with a semi-Lagrangian characteristics method are used. The results obtained with the three versions of the model are compared among them and show how the first of the simpler models fits better the experimental results.

Numerical Simulation of the Aerodynamic Behavior of High Velocity Trains under Synthetic Crosswinds of Different Shear and Turbulence Characteristics

A numerical simulation of the aerodynamic behavior of high-speed trains under synthetic crosswinds at a 90º yaw angle is presented. The train geometry is the aerodynamic train model (ATM). Flow description based on numerical simulations is obtained using large eddy simulation (LES) and the commercial code ANSYSFluent V14.5. A crosswind whose averaged velocity and turbulence characteristics change with distance to the ground is imposed. Turbulent fluctuations that vary temporally and spatially are simulated with TurbSim code. The crosswind boundary condition is calculated for the distance the train runs during a simulation period. The inlet streamwise velocity boundary condition is generated using Taylor?s frozen turbulence hypothesis. The model gives a time history of the force and moments acting on the train; this includes averaged values, standard deviations and extreme values. Of particular interest are the spectra of the forces and moments, and the admittance spectra. For comparison, results obtained with LES and a uniform wind velocity fluctuating in time, and results obtained with Reynolds averaged Navier Stokes equations (RANS), and the averaged wind conditions, are also presented.