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 modelling of metallic tubular crash structures with an explicit code

Numerous experimental studies have been carried out to investigate the collapse of tubular metallic crash structures under axial compression. Some simple theoretical models have been developed but these often assume one type of progressive collapse, which is not always representative of the real situation. Finite Element (FE) models, when further refined, have the potential to predict the actual collapse mode and how it influences the load-displacement and energy absorption characteristics. This paper describes an FE modelling investigation with the explicit code LS−DYNA. An automatic mesh generation programme written by the authors is used to set up shell and solid element tube models. Mesh specification issues and features relating to the contact and friction models are discussed in detail. The crush modes, load-deflection characteristics and energy absorption values found in the simulations are compared with a reasonable degree of correlation to those observed in a physical testing programme; however, improvements are still required.

Finite element analysis of an axi-symmetric forward spiral extrusion of Mg-1.75Mn alloy

A potential severe plastic deformation process known as axi-symmetrical forward spiral extrusion (AFSE) has been studied numerically and experimentally. The process is based on the extrusion of cylindrical samples through a die with engraved spiral grooves in a near zero shape change manner. The process was simulated using a three dimensional finite element (FE) model that has been developed using commercial software, ABAQUS. In order to verify the finite element results, hot rolled and annealed samples of the alloy were experimentally processed by AFSE. The required extrusion forces during the process were estimated using the FE model and compared with the experimental values. The reasonable agreement between the FE results and experimental data verified the accuracy of the FE model. The numerical results indicate the linear strain distribution in the AFSE sample is only valid for a core concentric while the strain distribution in the vicinity of the grooves is non axi-symmetric. The FE simulation results from this research allows a better understanding of AFSE kinematics especially near the grooves, the required extrusion force and the resultant induced strain distribution in the sample. To compare the mechanical properties of the Mg-1.75Mn alloy before and after the process, a micro shear punch test was used. The tests were performed on samples undergoing one and four passes of AFSE. After four passes of AFSE, it was observed that the average shear strength of the alloy has improved by about 21%. The developedfinite element model enables tool design and material flow simulation during the process.

Finite element analysis of plastic deformation in variable lead axisymmetric forward spiral extrusion

A modified axisymmetric forward spiral extrusion (AFSE) has been proposed recently to enhance the strain accumulation during the process. The new technique is called variable lead axisymmetric forward spiral extrusion (VLAFSE) that features a variable lead along the extrusion direction. To assess the effect of design modification on plastic deformation, a comprehensive study has been performed here using a 3D transient finite element (FE) model. The FE results established the shear deformation as the dominant mode of deformation which has been confirmed experimentally. The variable lead die extends strain accumulation in the radial and longitudinal directions over the entire grooved section of the die and eliminates the rigid body rotation which occurs in the case of a constant lead die, AFSE. A comparison of forming loads for VLAFSE and AFSE proved the advantages of the former design in the reduction of the forming load which is more pronounced under higher frictional coefficients. This finding proves that the efficiency of VLAFSE is higher than that of AFSE. Besides, the significant amount of accumulated shear strain in VLAFSE along with non-axisymmetric distribution of friction creates a surface feature in the processed sample called zipper effect that has been investigated. © 2012 Springer Science+Business Media New York.

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.

Multi-scale simulation using a hybrid mass-spring system and finite element model

Ply-scale finite element (FE) models are widely used to predict the performance of a composite structure based on material properties of individual plies. When simulating damage, these models neglect microscopic fracture processes which may have a significant effect on how a crack progresses within and between plies of a multidirectional laminate. To overcome this resolution limitation a multi-scale modelling technique is employed to simulate the effect micro-scale damage events have on the macro-scale response of a structure. The current paper discusses the development and validation of a hybrid mass-spring system and finite element modelling technique for multi-scale analysis. The model developed here is limited to elastic deformations; however, it is the first key step towards an efficient multi-scale damage model well suited to simulation of fracture in fibre reinforced composite materials. Various load cases have been simulated using the model developed here which show excellent accuracy compared to analytical and FE results. Future work is discussed, including extension of the model to incorporate damage modelling.

Application of the strain localization CAFE model to investigate extrusion with various die shapes

Multi scale CAFE model for the prediction of initiation and propagation of the micro shear bands and shear bands in metallic materials subjected to plastic deformation is presented. The CAFE approach is the combination of the Cellular Automata (CA) and the Finite Element (FE) methods. The application of the developed CAFE model to analyze material flow during extrusion is the objective of the present work. The proposed CAFE approach is applied in this work to simulation of the extrusion with flat face and convex dies and to investigate differences in the material flow. The initial FE meshes with the set of the CA point are generated for the numerical tests and the results of the metal flow predicted by the CAFE method are presented in the paper.

Slope stability charts for two-layered purely cohesive soils based on finite-element limit analysis methods

Stability charts for soil slopes, first produced in the first half of the twentieth century, continue to be used extensively as design tools, and draw the attention of many investigators. This paper uses finite-element upper and lower bound limit analysis to assess the short-term stability of slopes in which the slopematerial and subgrade foundation material have two distinctly different undrained strengths. The stability charts are proposed, and the exact theoretical solutions are bracketed to within 4.2% or better. In addition, results from the limit-equilibrium method (LEM) have been used for comparison. Differences of up to 20% were found between the numerical limit analysis and LEM solutions. It also shown that the LEM sometimes leads to errors, although it is widely used in practice for slope stability assessments.

Three-dimensional Finite Element modelling of laser shock peening process

Laser shock peening (LSP) is an innovative surface treatment technique for metal alloys, with the great improvement of their fatigue, corrosion and wear resistance performance. Finite element method has been widely applied to simulate the LSP to provide the theoretically predictive assessment and optimally parametric design. In the current work, 3-D numerical modelling approaches, combining the explicit dynamic analysis, static equilibrium analysis algorithms and different plasticity models for the high strain rate exceeding 106s-1, are further developed. To verify the proposed methods, 3-D static and dynamic FEA of AA7075-T7351 rods subject to two-sided laser shock peening are performed using the FEA package–ABAQUS. The dynamic and residual stress fields, shock wave propagation and surface deformation of the treated metal from different material modelling approaches have a good agreement.

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.

Improving stamping tryout using variational analysis, FE simulations and robustness

The tryout process at an automotive sheet metal production plant is outlined and analysed with the aim of augmenting the process to improve robustness through an understanding of variation using Finite Element (FE) analysis and other tools. The result is tryout can produce a more robust stamping process.

Slope stability analysis for filled slopes using finite element limit analysis method

This paper uses the finite element upper and lower bound limit analysis to assess the stability of slopes mostly found in embankment cases where frictional materials are filled on purely cohesive undrained clay. For comparison purposes, the commonly used stability assessment method, limit equilibrium method (LEM) is also employed. The final results for both methods are then presented in the form of comprehensive chart solutions for the convenience of practicing engineers during preliminary slope designs. The failure mechanism will also be discussed in this paper. Ultimately, it should be noted that finite element limit analysis method holds the upper hand as its prior assumptions are not required. Thus, the obtained failure mechanism from the slope stability analysis will be more realistic. Hence, it will provide a better understanding for the slope failure surface. Therefore, engineers should design more carefully when the LEM is applied to the slopes with frictional materials filled on purely cohesive undrained clay. © 2014 American Society of Civil Engineers.