Uday Shanker er Dix ixit it Dep epartmen ent o
- f mec
echanical E Enginee eering Ind ndia ian I Inst nstitute o
- f Techn
hnolo logy G Guw uwaha hati i
1
Uday Shanker er Dix ixit it Dep epartmen ent o of mec - - PowerPoint PPT Presentation
Uday Shanker er Dix ixit it Dep epartmen ent o of mec echanical E Enginee eering Ind ndia ian I Inst nstitute o of Techn hnolo logy G Guw uwaha hati i 1 Introduction Processes that cause changes in the shape of solid
Uday Shanker er Dix ixit it Dep epartmen ent o
echanical E Enginee eering Ind ndia ian I Inst nstitute o
hnolo logy G Guw uwaha hati i
1
Processes that cause changes in the shape of solid metal parts via plastic (permanent) deformation are termed as metal forming processes.
In bulk metal forming processes, raw material and products have a relatively high ratio of volume to surface area.
Introduction
Rollin lling Extr trusion
2
3
Open en-die ie f forgin ing Clo losed-di die f forg rgin ing
Introduction
through conical dies resulting in the reduction of its cross-section and increase in its length.
Wire drawing process
4
In sheet metal forming processes, raw material and products have a relatively low ratio of volume to surface area.
Introduction
5
Deep drawing process. a Before deformation. b After deformation
6
Schematic of laser forming process
Metal forming processes are also used just for improving the properties of the material.
Sev ever ere Pl Plast stic Defo eformation Pr Proces esses
Introduction
7
8
Hydraulic Autofrettage
material during processing
9
Pre-FEM Techniques
approximated with the deformation of a series of slabs.
processes-
Stresses in the slab for a rolling process
in the analysis of various metal forming processes such as forging, rolling and extrusion.
it is computationally very fast.
10
2 2 2
x y xy y y
for Tresca criterion. for von Mises criterion.
is assumed to be rigid plastic.
For an isotropic material, the yield criterion may be written as
Pre-FEM Techniques
11
xy x xy y y x
Pre-FEM Techniques Pre-FEM Techniques
12
y y x x x y xy x y xy x y
The principal directions for the stresses and strain rate are same. Hence Unknowns ( total 5 equations)
Pre-FEM Techniques
13
Pre-FEM Techniques
velocity are in the direction of maximum shear stresses.
consider a pair of orthogonal curves along which the shear stress has its maximum value. These curves are called slip lines or shear lines.
(a) Stresses on the element (b) Mohr’s circle
(a) (b)
14
constant along an α line constant along a β line along an α line along a β line
Pre-FEM Techniques Pre-FEM Techniques
Thus, for a straight slip line, the hydrostatic pressure and the tangential velocity remains constant.
15
Slip line fields for the indentation of a semi-infinite medium by a punch
According to Henky’s equations, the punch pressure at yield point has the value
Pre-FEM Techniques
16
*
d d d d d d d d σ ε
∗ ∗ ∗ ∗
≤ + −
u T D
i i u ij ij D i i T S v S S
t u S v k u S t u S
metal forming can deliver or sustain.
is less than or equal to the rate of internal energy dissipated in any kinematic admissible field.
Pre-FEM Techniques
17
displacements are undergone, the increment of work done by the actual forces or surface tractions on Su is greater than
surface tractions of any other statically admissible stress field, *
u u
i i u i i S S
upper bound method.
Pre-FEM Technique
Domain with surface tractions
18
Pre-FEM Techniques
process.
the considerations of the equilibrium and constitutive equations. The method has limited applications.
19
Continuity equation Equilibrium equations Strain-rate-velocity relations Incremental strain-displacement relations Constitutive relations
20
ii xx yy zz
, , , ij j i i ij j i
effective approach than an Eulerian approach for metal forming.
large displacements, new control volumes have to be created (because the boundaries of the solid change continuously) and the non-linearities in the convective acceleration terms are difficult to deal with.
FEM Technique
non-homogeneity of deformation, process dependent material properties and different friction models.
21
effective approach than an Eulerian approach for metal forming.
large displacements, new control volumes have to be created (because the boundaries of the solid change continuously) and the non-linearities in the convective acceleration terms are difficult to deal with.
FEM Technique
non-homogeneity of deformation, process dependent material properties and different friction models.
22
FEM Technique
coordinate system, the relation between the spatial position (x) and the initial coordinates (X) and time (t) is given by x = x (X, t)
motion in a Lagrangian formulation.
configuration to which all future variables are referred.
23
Updated Lagrangian formulation:
formulation proves to be cumbersome with the governing equations being difficult to solve.
states of stress and deformation of the body are known till the current configuration, say time τ.
deformation and stresses during the time step ∆ t and the state of the system in a future configuration, at time τ+∆ t.
x = x (x(τ), τ+∆ t)
FEM Technique
24
More computational time. Some times, residual stresses are not predicted
properly, due to accumulation of computational errors
Difficulty in applying spatial boundary condition.
FEM Technique
25
Eulerian formulation:
analysis of steady-state processes.
whose independent variables are present position x
approach is called flow formulation.
finite element analysis are the velocities/ hydrostatic pressure at different locations i.e. nodes.
FEM Technique
26
ij ij ij y
i j j ,i ,
The following relations exists between strain rates and the velocity gradient:
FEM Technique
where vi,j is the partial derivative of the ith component of velocity with respect to jth component of position vector. Levi-Mises plastic flow rule is given by, μ is the Levi-Mises coefficient,
27
y y
ij ij n
is equivalent strain which is equal to (σy)0 is the flow stress at zero strain while b, n are determined from experiments. Neglecting the effect of strain rate and temperature, flow stress can be expressed as
FEM Technique FEM Technique
28
FEM Technique
vi,i = 0
where p is the hydrostatic pressure.
different metal forming processes using Galerkin (weak) formulation.
the above method difficult.
29
, i i
by writing it as where λ is a penalty parameter.
conditioned system of equations.
FEM Technique
30
FEM Technique
analysis of wire drawing, extrusion, rolling etc. Example of a cold flat rolling analysis using flow formulation.
Discretised metal strip for FEM analysis [ Dixit and Dixit, 1996 ]
31
Galerkin method has been used and the global equations after applying the boundary conditions, are in the form of; [K]{∆} = {F}
FEM Technique
where [K] is the global coefficient matrix, {F} is the global right hand side vector, {∆} is the global vector of nodal values of pressure and velocity (primary variables).
velocities and pressure, then the secondary variables like roll force and roll torque can be computed.
32
Form
latio ion ma may n not re
any p pre ressure bounda dary c conditions.
The p
e pressure v val alues m may ay b be d e det etermined w with an a additive co
The c e constant c can an b be e elimin liminated f from trac rom tracti tion c con
tions.
In th
the d dire irect p t penalty ty me meth thod:
33
ii
34
Pre ressure compu putatio ion by by FDM FDM The he d doma
n show howing p plasti tic b bound
and nd the the p poi
nts f for
fin init ite dif difference a appro pproxim imatio ion
Gudur, P.P. and Dixit, U.S., (2008), A combined finite element and finite difference analysis of cold flat rolling, Transactions of the American Society of Mechanical Engineers: Journal of Manufacturing Science and Engineering, 130, 011007, pp.1–6.
35
Comparison of FEM and experimental results for roll force and roll torque (Steel, R/h1 = 65)
Gudur, P.P. and Dixit, U.S., (2008), A combined finite element and finite difference analysis of cold flat rolling, Transactions of the American Society of Mechanical Engineers: Journal of Manufacturing Science and Engineering, 130, 011007, pp.1–6.
36
Comparison of FEM and experimental results for roll force and roll torque (Steel, R/h1 = 130)
Gudur, P.P. and Dixit, U.S., (2008), A combined finite element and finite difference analysis of cold flat rolling, Transactions of the American Society of Mechanical Engineers: Journal of Manufacturing Science and Engineering, 130, 011007, pp.1–6.
37
Comparison of FEM and experimental results for roll force and roll torque (Copper)
Gudur, P.P. and Dixit, U.S., (2008), A combined finite element and finite difference analysis of cold flat rolling, Transactions of the American Society of Mechanical Engineers: Journal of Manufacturing Science and Engineering, 130, 011007, pp.1–6.
38
Compa paris rison of ro roll pre ll pressure dis distribu ributio ion f for r dif different m methods ds
Gudur, P.P. and Dixit, U.S., (2008), A combined finite element and finite difference analysis of cold flat rolling, Transactions of the American Society of Mechanical Engineers: Journal of Manufacturing Science and Engineering, 130, 011007, pp.1–6.
Finding the l
load o
mation
Inte tegrati ation
es v versus E Energy gy b balance e
Calc
lcula latio ion of
l stre tresses an and s shear stre tresses a at t th the tool tool-wor
inte terface ce Should I d I find d ts by ab abov
ion or
irectl tly from rom str tress comp compon
39
40
U.S. Dixit, P.M. Dixit Journal of Materials Processing Technology 47(1995), 201-229
41
42
U.S. Dixit, P.M. Dixit Journal of Materials Processing Technology 47(1995), 201-229
rolling, [Chandra and Dixit, 2004], showed that the non- homogeneous deformation causes region of very low strain-rate (less than 3% of maximum strain-rate), indicating the presence of elastic zone in between the plastic zone.
FEM Technique
Equivalent strain-rate contours for r = 4 %, μ = 0.3, R/h1= 130)
43
44
Variation of tangential stress with normal stress according to Wanheim and Bay friction model
FEM Technique: Anisotropy
values of equivalent strain-rate and equivalent strain decrease with m which is a parameter that governs the shape of the yield surface.
Equivalent strain-rate contours for two different cases of m m = 1.5 m = 2
45
FEM Technique: Residual Stresses
forming process.
find the longitudinal residual stress (stress in the direction
method.
mixed formulation, method of multiplicative decomposition of the deformation and the rate formulation were discussed for analyzing elasto-plastic rolling problem.
Equivalent stress contours
46
Thermal Modelling: warm flat rolling
Constitutive relation Strip temperature Slip Based on Eulerian flow formulation Roll-force, Roll-torque, Distributions
Based on Analytical methods Based on FEM using ABAQUS
Deformation module Thermal module
Friction model Thermal parameters Average temperature of roll and strip at the interface
Thermo-mechanical model for flat rolling
Roll temperature Heat partition factor λ and/or
An overview of thermo-mechanical model of plane strain rolling
O
Roll Strip
a b ω
A B C D
Plane of symmetry
V2 V1
x' y'
Heat generated due to plastic work Heat generated due to friction work Heat conducted to roll h1 2 h2 2
E
Convective heat losses in air
F G H
x y 2β
Schematic of upper half view of the roll and strip
Thermal Autofrettage
(a) (b)
A schematic diagram of thermal autofrettage process, (b) Elastic-plastic zones across the wall thickness
49
50
51
52
53
54
55
56
improve predictions of surface finish and friction and thus increase productivity and improve quality.
Finite Element Modeling of Surface Roughness Transfer and Oxide Scale Micro Deformation in Rolling Process
Sta tate of
under mi mixed lu lubric bricatio ion [ [Xie ie et a t al. [201 [2011]
H. H.B.
ie, Z , Z.Y. J . Jiang, W , W.Y.D .D. Y Yuen, ( , (2011), Tribo ibolo logy Intern rnatio ional.
44, pp.
971
979 . 9 .
Sur urface r roug
trip sample ples
Cry rystal p pla lasticity fin inite e ele lement m modelling o
surf rface ro roughness in in ro roll lling p pro rocess [ [Jian ang e g et al. 2013]
ang, D
and H
(2013 2013), AIP, Co
1532, p
254-261 261.
th an n inc ncreases of
uction, the the s sur urface r roug
de decreases s sig ignif ific icantly ly.
bricatio ion c can de dela lay t the pro process o
rface a aspe perit ity f fla lattening.
ing stra rain in ra rate c can le lead t d to a de decrease o
rface ro roughness unde der the same re redu ductio ion. . Effect of
rol
urface roug
at t rol
ng te temp mperature 850 850 °C Effect of
rol
urface roug
at t rol
ng te temp mperature 950 950 °C
Cry rystal p pla lasticity fin inite e ele lement m modelling o
surf rface ro roughness in in ro roll lling p pro rocess [ [Jian ang e g et al. 2013]
ang, D
and H
(2013 2013), AIP, Co
1532, p
254-261 261.
Relation between the surface roughness and friction/reduction (a) with lubrication and (b) without lubrication
Reduction (%) Measured surface roughness after cold rolling Ra (µm) Measured surface roughness after warm rolling Ra (µm) T = 150 °C T = 200 °C Rolling direction Transverse direction Rolling direction Transverse direction Rolling direction Transverse direction
8.2 0.411 (0.026) 0.379 (0.014) 0.521 (0.022) 0.512 (0.012) 0.532 (0.021) 0.545 (0.015) 12.4 0.368 (0.024) 0.349 (0.025) 0.541 (0.037) 0.510 (0.021) 0.556 (0.032) 0.534 (0.029) 22.2 0.336 (0.033) 0.299 (0.018) 0.557 (0.054) 0.531 (0.027) 0.598 (0.045) 0.551 (0.035) 32.6 0.288 (0.012) 0.281 (0.021) 0.571 (0.039) 0.541 (0.044) 0.561 (0.039) 0.495 (0.016)
he f friction
nd sur urface r roug
effects on
crack e evol
thin strip ro rip rollin lling.
For r im impro provin ing t the qu qualit lity o
rolle lled s d sheet, t , the s surf rface ro roughness a and d fricti tion
are the the mos most t imp mpor
tant t contr
parame meters.
rface ro roughness o
rolle lled d strip rip was m measured u d usin ing t the s surf rface roughne hness m measur uring ng i instrume ument ( nt (Pocket S t Surf) f).
Averaged measured surface roughness of the strip after rolling at different reductions for h1 = 5 mm, R = 100 mm (Values in bracket are standard deviations)
Experimental study (Yadav, 2016)
61
A laboratory rolling mill (a) front view and (b) arrangement for measuring the temperature and velocity of exit strip at the rear side
Inverse Modelling for estimation of parameters
Start
Guess the initial value of mechanical properties and friction Direct model Temperature and slip Experiment (In-house warm rolling experiments) Is the error in measured data between computed data < a prescribed value (1%)? Stop
Temperature and slip Yes No Updated the mechanical properties using optimization algorithm (A heuristic based optimization algorithm is used here.) Record the mechanical properties and friction
Flow chart of the inverse model for mechanical properties and friction coefficient determination Experimental validation of strategy for inverse estimation
.
n f eq m
T T
γ
σ σ ε
−
=
The following objective function is to be minimized with respect to the decision variables σ0, n, and µ such as
2 1
1
n ism isc i ism
T T E n T
=
− = ×
∑
σ
Comparison of the experimental stress-strain curves and inversely estimated stress-strain curves at different temperatures with upper and lower bound fit data (a) Room temperature, (b) 100 °C, (c) 150 °C and (d) 250 °C
Meshless Methods
equations on regular or irregular distribution of points.
equations, whilst others use collocation points.
speed, even though in some cases, the quality of solution has been shown to be better than that obtained by FEM .
solving the problems in metal forming industry.
Chen et al. have used RKPM for the analysis of ring
compression and upsetting problem. Shangwu et al. utilized RKPM method for modeling of plane strain rolling problem.
Bonnet and Kulasegaram et al. utilized the CSPH to
perform two-dimensional simulations of several metal forming processes without considering strain-hardening.
Xiang et al. applied element free Galerkin method for simulation
Meshless Methods
using this method.
equations are satisfied at the collocation points.
φ
Input
Hidden layer Output layer
1
w
2
w
3
w
4
w
5
w
2
x
1
x
3
x
Architecture of a typical radial basis function neural network
( )
|| ||,
i i k k
w g c φ
∑
= − x x
Multiquadric Gaussian Inverse multiquadric Thin-plate splines
2 2
exp
k
r g c = −
2 2 k
g r c = +
2 2
1
k
g r c = +
2 ln m
g r r = Meshless Methods
Simulation of the process using RBF collocation method
Meshless Methods
Problem domain
Meshless Methods
Non-dimensional equivalent strain rate at 20 % reduction in height with f = 0.05 and R/H = 1/2.
. 1 8 8 5 6 2 . 2 1 7 2 0.209608 . 2 1 5 7 6 . 2 2 8 3 4 0.225428 . 2 2 7 3 1 2 0.220834 0.21576 0.209608 0.20172 0.188562 0.178104 0.20172 . 1 7 8 1 4
Radial distance Axial distance
0.1 0.2 0.3 0.4 0.5 0.6 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8
Equivalent plastic strain at 20 %
reduction in height with f = 0.05 and R/H = 1/2.
1.16435 1.22282 1.25792 1 . 2 9 5 4 1 1.31762 1.35336 1 . 2 2 2 8 2 1 . 1 6 4 3 5 1.25792 1.16435 1.22282 1.25792 1.0741 0.97147 1.0741 0.97147 1.0741
Radial distance Axial distance
0.1 0.2 0.3 0.4 0.5 0.6 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8
. 5 4 2 9 5 3 . 4 3 3 9 . 3 9 5 2 2 9 0.352159 . 3 2 2 6 3 8 0.299806 0.31063 0.322638 0.352159 0.395229 . 4 3 3 9 . 4 9 7 8 4 9 . 5 4 2 9 5 3 . 3 1 6 3 0.299806
Radial distance Axial distance
0.1 0.2 0.3 0.4 0.5 0.6 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8
. 4 6 8 6 8
. 7 1 1 86
Radial distance Axial distance
0.1 0.2 0.3 0.4 0.5 0.6 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8
Contours of hydrostatic pressure at 20 % reduction in height with f = 0.05 and R/H = 1/2 Contours of normal stress at 20 % reduction in height with f = 0.05 and R/H = 1/2
Meshless Methods
71
Computational Efficiency: Hybrid Methods
Further Research Directions
from saturation.
research is needed.
Further Research Directions
material has to be understood.
be developed.
problem is of utmost important for developing a suitable computational techniques.
method requires a lot of memory and computational time.
will speed up the computation of the FEM techniques.
Further Research Directions
the modeling of metal forming.
results in the case of foil rolling process.
unless more is known about the nature of friction in the roll bite. [Fleck et al., 1992].
friction model, i.e. Wanheim and Bay’s model for cold foil rolling and observed that friction model has great influence on qualitative and quantitative predictions of foil rolling processes.
Further Research Directions
area.
recrystallization and grain growth phenomenon.
modeling is expected to provide better judgment for the analysis of the defects.
useful results for industries.
Further Research Directions
experimental as well as computational work has to be carried out.
prediction of the surface roughness including the type
forming process.
exploring the various friction models as friction plays a major role in determining the surface roughness.
77