AN EXTENDED ARBITRARY LAGRANGIAN-
EULERIAN FINITE ELEMENT MODELING IN
POWDER FORMING PROCESS
Presented by Submitted to
Pawan Kumar Dr. R.k. Patel
(2019DN21)
Objective
• To simulate the high deformation in powder compaction at
the boundary of die and powder interface.
Powder Compaction
• It is process of compacting material powder in a die
through the application of high pressure.
ALE
• Arbitrary Lagrangian-Eularian : It is a finite
element formulation in which the computational
system is not a prior fixed in space (e.g. Eulerian-
based finite element formulations) or attached to
material (e.g. Lagrangian-based finite element
formulations).
• Need Of ALE: ALE based simulations can
eliminate many of the drawback of traditional
lagrangian-based and Euleian-based finite
element simulation have.
• ALE-based finite element formulations can
reduce to either Lagrangian-based finite element
formulations by equating mesh motion to material
motion or Eulerian-based finite element
formulations by fixing mesh in space.
• Material motion: xm = f(X, t)
• ALE motion: xg = f(X, 0)
• Mesh motion: xm = ˆf (xg, t)
• Mesh displacement : ug(xg, t) = xm − xg
• Mesh velocity: vg(xg, t) = (∂xm/∂t)|xg
• In ALE formulation, the convective velocity c is
defined using the difference between the material
and mesh velocities as:
ALE formulation the momentum equation can
• In
be written as:
where, ρ= Density
σ= Cauchy’s stress
b= Body forces
Petrov-Galerkin method
• The Petrov–Galerkin method is a mathematical method
used to obtain approximate solutions of partial differential
equations which contain terms with odd order.
• Different space for test function and solution function.
• Weak form of momentum equation:
Where: δvi= Galerkin test function
δvPG= Petrov Galerkin test function
Ω= Spatial domain
• The Petrov–Galerkin test function can be related
to nodal values using the Galerkin shape
function: Ni’=Ni+kcjNi,j.
where , Ni= Galerkin shape function.
Ni’= Petrov Galerkin shape function
• Internal nodal force can be computed in
referential domain as:
where, Ωˆ: ALE domain
• Applying the finite element discretization, the ALE
stiffness matrix can expressed as:
• In Lagrangian approach quadrature point
coincide with material point regardless of
deformation of continuum.
• In ALE approach Gauss point does not coincide
with material point.
X-FEM
• It is a numerical technique that extends the
classical finite element method (FEM) approach
by extending the solution space for solutions to
differential equations with discontinuous
functions.
• The standard FE approximation can be enriched
with additional functions by using the X-FEM.
Where, ui= classical nodal displacement
aj= Nodal dof
ψ(x)= Enrichment function
• (a): The X-FEM method is used to discretize the
intersected elements by partitioning the domain with some
triangular sub-elements.
• (b)The enrichment of displacement field is
applied to correct the standard displacement
approximation by incorporating discontinuous
fields through X-FEM method.
• The choice of enrichment function is depend on
the condition of problem.
• Applying the level set enrichment function ψ(x)
based on its nodal values as:
• strain matrix relating the increments of strain and
displacement can be written as:
• the discrete system of equations can be obtained
as;
Kd=f
where, K= Global Stiffness matrix
d= Vector of unknown of nodal point ui and aj
at nodal point.
f= External force.
Conclusion
• the extended finite element method was presented in the
framework of arbitrary Lagrangian–Eulerian finite element
model to simulate the large deformation in powder
compaction processes.
• The technique was applied by employing enrichment
functions to approximate the displacement field of the
elements located on the interface.
• The technique was applied to take benefits of both
Lagrangian and Eulerian approach both.
References
• An extended arbitrary Lagrangian–Eulerian finite element
modeling(X-ALE–FEM) in powder forming processes
:A.R. Khoei ∗, M. Anahid, K. Shahim
• R.W. Lewis, A.R. Khoei, Numerical modeling of large
deformation in metal powder forming, Comput. Methods
Thank you