Finite element analysis of drilling of titanium alloy

Drilling is a highly demanding machining process due to complex tool geometry and the progressive material failure on the work piece. In this study, a 3D model is developed using commercial finite element software ABAQUS/Explicit. The proposed model simulates the drilling process by taking into account of the damage initiation and evolution of the work piece material, a contact model at the interface between drill bit and work piece and the process parameters. The results of the simulations demonstrate the effects of machining parameters on drilling. The results also confirm the capability and advantage of FE simulation of the drilling process. © 2011 Published by Elsevier Ltd.

Impact of finite element idealisation on the prediction of welded fuselage stiffened panel buckling

Lap joints are widely used in the manufacture of stiffened panels and influence local panel sub-component stability, defining buckling unit dimensions and boundary conditions. Using the Finite Element method it is possible to model joints in great detail and predict panel buckling behaviour with accuracy. However, when modelling large panel structures such detailed analysis becomes computationally expensive. Moreover, the impact of local behaviour on global panel performance may reduce as the scale of the modelled structure increases. Thus this study presents coupled computational and experimental analysis, aimed at developing relationships between modelling fidelity and the size of the modelled structure, when the global static load to cause initial buckling is the required analysis output. Small, medium and large specimens representing welded lap-joined fuselage panel structure are examined. Two element types, shell and solid-shell, are employed to model each specimen, highlighting the impact of idealisation on the prediction of welded stiffened panel initial skin buckling.

Finite element modeling of debonding failures in FRP-strengthened RC beams: A dynamic approach

This paper is concerned with the finite element simulation of debonding failures in FRP-strengthened concrete beams. A key challenge for such simulations is that common solution techniques such as the Newton-Raphson method and the arc-length method often fail to converge. This paper examines the effectiveness of using a dynamic analysis approach in such FE simulations, in which debonding failure is treated as a dynamic problem and solved using an appropriate time integration method. Numerical results are presented to show that an appropriate dynamic approach effectively overcomes the convergence problem and provides accurate predictions of test results.

Finite Element Analysis of the Enhanced Strength of Laterally Restrained RC Slabs

This paper presents the numerical simulation of the ultimate behaviour of 85 one-way and two-way spanning laterally restrained concrete slabs of variable thickness, span, reinforcement ratio, strength and boundary conditions reported in literature by different authors. The developed numerical model was described and all the assumptions were illustrated. ABAQUS, a Finite Element Analysis suite of software, was employed. Non-linear implicit static general analysis method offered by ABAQUS was used. Other analysis methods were also discussed in general in terms of application such as Explicit Dynamic Analysis and Riks method. The aim is to demonstrate the ability and efficacy of FEA to simulate the ultimate load behaviour of slabs considering different material properties and boundary conditions. The authors intended to present a numerical model that provides consistent predictions of the ultimate behaviour of laterally restrained slabs that could be used as an alternative for expensive real life testing as well as for the design and assessment of new and existing structures respectively. The enhanced strength of laterally-restrained slabs compared with conventional design methods predictions is believed to be due to compressive membrane action (CMA). CMA is an inherent phenomenon of laterally restrained concrete beams/slabs. The numerical predictions obtained from the developed model were in good correlation with the experimental results and with those obtained from the CMA method developed at the Queen’s University Belfast, UK.

Generic Finite Element Model for Mechanically Consistent Scaling of Composite Beam-to-column Joint with Welded Connection

To study the behaviour of beam-to-column composite connection more sophisticated finite element models is required, since component model has some severe limitations. In this research a generic finite element model for composite beam-to-column joint with welded connections is developed using current state of the art local modelling. Applying mechanically consistent scaling method, it can provide the constitutive relationship for a plane rectangular macro element with beam-type boundaries. Then, this defined macro element, which preserves local behaviour and allows for the transfer of five independent states between local and global models, can be implemented in high-accuracy frame analysis with the possibility of limit state checks. In order that macro element for scaling method can be used in practical manner, a generic geometry program as a new idea proposed in this study is also developed for this finite element model. With generic programming a set of global geometric variables can be input to generate a specific instance of the connection without much effort. The proposed finite element model generated by this generic programming is validated against testing results from University of Kaiserslautern. Finally, two illustrative examples for applying this macro element approach are presented. In the first example how to obtain the constitutive relationships of macro element is demonstrated. With certain assumptions for typical composite frame the constitutive relationships can be represented by bilinear laws for the macro bending and shear states that are then coupled by a two-dimensional surface law with yield and failure surfaces. In second example a scaling concept that combines sophisticated local models with a frame analysis using a macro element approach is presented as a practical application of this numerical model.

An adaptive finite element water column model using the Mellor-Yamada level 2.5 turbulence closure scheme

A one-dimensional water column model using the Mellor and Yamada level 2.5 parameterization of vertical turbulent fluxes is presented. The model equations are discretized with a mixed finite element scheme. Details of the finite element discrete equations are given and adaptive mesh refinement strategies are presented. The refinement criterion is an "a posteriori" error estimator based on stratification, shear and distance to surface. The model performances are assessed by studying the stress driven penetration of a turbulent layer into a stratified fluid. This example illustrates the ability of the presented model to follow some internal structures of the flow and paves the way for truly generalized vertical coordinates. (c) 2005 Elsevier Ltd. All rights reserved.

A finite element simulation of micro-mechanical frictional behaviour in metal forming

Friction is a critical factor for sheet metal forming (SMF). The Coulomb friction model is usually used in most finite element (FE) simulation for SMF. However, friction is a function of the local contact deformation conditions, such as local pressure, roughness and relative velocity. Frictional behaviour between contact surfaces can be based on three cases: boundary, hydrodynamic and mixed lubrication. In our microscopic friction model based on the finite element method (FEM), the case of dry contact between sheet and tool has been considered. In the view of microscopic geometry, roughness depends upon amplitude and wavelength of surface asperities of sheet and tool. The mean pressure applied on the surface differs from the pressure over the actual contact area. The effect of roughness (microscopic geometric condition) and relative speed of contact surfaces on friction coefficient was examined in the FE model for the microscopic friction behaviour. The analysis was performed using an explicit FE formulation. In this study, it was found that the roughness of deformable sheet decreases during sliding and the coefficient of friction increases with increasing roughness of contact surfaces. Also, the coefficient of friction increases with the increase of relative velocity and adhesive friction coefficient between contact surfaces.

Finite element simulation of deformation processes of micro-components

The aim of this paper is to improve the understanding of deformation of micro medical needle and thread during assembly and then to develop an economical and flexible deformation method. Therefore, the swaging process is computationally simulated with the finite element method in this paper. A commercially available explicit nonlinear finite element analysis code, LS-Dyna, is used to model the 3-D deformation and contact problem. As the firmness of the assembly on the needle depends on the contact force and friction, the contact and the slide between the needle and thread are taken into account in the simulation. The general surface-to-surface contact algorithm (STS) is used to simulate the contact. The paper provides an insight into the deformation of the micro products.

Application of radial return mapping algorithm for finite strain elastoplasticity in a hybrid finite element and particle-in-cell model

The radial return mapping algorithm within the computational context of a hybrid Finite Element and Particle-In-Cell (FE/PIC) method is constructed to allow a fluid flow FE/PIC code to be applied solid mechanic problems with large displacements and large deformations. The FE/PIC method retains the robustness of an Eulerian mesh and enables tracking of material deformation by a set of Lagrangian particles or material points. In the FE/PIC approach the particle velocities are interpolated from nodal velocities and then the particle position is updated using a suitable integration scheme, such as the 4th order Runge-Kutta scheme[1]. The strain increments are obtained from gradients of the nodal velocities at the material point positions, which are then used to evaluate the stress increment and update history variables. To obtain the stress increment from the strain increment, the nonlinear constitutive equations are solved in an incremental iterative integration scheme based on a radial return mapping algorithm[2]. A plane stress extension of a rectangular shape J2 elastoplastic material with isotropic, kinematic and combined hardening is performed as an example and for validation of the enhanced FE/PIC method. It is shown that the method is suitable for analysis of problems in crystal plasticity and metal forming. The method is specifically suitable for simulation of neighbouring microstructural phases with different constitutive equations in a multiscale material modelling framework.

Statistical analysis of finite element modeling in sheet metal forming and springback analysis

This work focuses on development of a method to statistically study forming and springback problems of TRansformation Induced Plasticity (TRIP) through an industrial case study. A Design of Experiments (DOE) approach was used to study the sensitivity of predictions to four user input parameters in implicit and explicit sheet metal forming codes. Numerical results were compared to experimental measurements of parts stamped in an industrial production line. The accuracy of forming strain predictions for TRIP steel were comparable with conventional steel, but the springback predictions of TRIP steel were far less accurate. The statistical importance of selected parameters for forming and springback prediction is also discussed. Changes of up to ±10% in Young's modulus and coefficient of friction were found to be insignificant in improving or deteriorating the statistical correlation of springback accuracies.

Spectral element modelling of wave propagation with boundary and structural discontinuity reflections

Spectral element method is very efficient in modelling high-frequency stress wave propagation because it works in the frequency domain. It does not need to use very fine meshes in order to capture high frequency wave energy as the time domain methods do, such as finite element method. However, the conventional spectral element method requires a throw-off element to be added to the structural boundaries to act as a conduit for energy to transmit out of the system. This makes the method difficult to model wave reflection at boundaries. To overcome this limitation, imaginary spectral elements are proposed in this study, which are combined with the real structural elements to model wave reflections at structural boundaries. The efficiency and accuracy of this proposed approach is verified by comparing the numerical simulation results with measured results of one dimensional stress wave propagation in a steel bar. The method is also applied to model wave propagation in a steel bar with not only boundary reflection, but also reflections from single and multiple cracks. The reflection and transmission coefficients, which are obtained from the discrete spring model, are adopted to quantify the discontinuities. Experimental tests of wave propagation in a steel bar with one crack of different depths are also carried out. Numerical simulations and experimental results show that the proposed method is effective and reliable in modelling wave propagation in one-dimensional waveguides with reflections from boundary and structural discontinuities. The proposed method can be applied to effectively model stress wave propagation for structural damage detection.

Finite element model updating using estimation of distribution algorithm

Finite Element (FE) model updating has been attracting research attentions in structural engineering fields for over 20 years. Its immense importance to the design, construction and maintenance of civil and mechanical structures has been highly recognised. However, many sources of uncertainties may affect the updating results. These uncertainties may be caused by FE modelling errors, measurement noises, signal processing techniques, and so on. Therefore, research efforts on model updating have been focusing on tackling with uncertainties for a long time. Recently, a new type of evolutionary algorithms has been developed to address uncertainty problems, known as Estimation of Distribution Algorithms (EDAs). EDAs are evolutionary algorithms based on estimation and sampling from probabilistic models and able to overcome some of the drawbacks exhibited by traditional genetic algorithms (GAs). In this paper, a numerical steel simple beam is constructed in commercial software ANSYS. The various damage scenarios are simulated and EDAs are employed to identify damages via FE model updating process. The results show that the performances of EDAs for model updating are efficient and reliable.