MEMORIAS DEL XXVII CONGRESO INTERNACIONAL ANUAL DE LA SOMIM

22 al 24 DE SEPTIEMBRE DE 2021 PACHUCA, HIDALGO, MÉXICO

Tema A2b.Manufactura: Formado de lámina metálica

“Análisis y modelado numérico del proceso de embutido de lámina metálica de

formas elípticas”

Hugo I. Medellín Castillo*, Edmar Álvarez Sánchez

Facultad de ingeniería, Universidad Autonoma de San Luis Potosí, Av. Manuel Nava No. 8, Zona Universitaria, 78290, S.L.P., Mé xico
* Autor de contacto: [email protected]

RESUMEN

El proceso de embutido de lámina es ampliamente utilizado en la industria automotriz, aeronáutica, de electrodomésticos, y

en general en toda la industria metalmecánica debido a las características mejoradas en términos de resistencia/peso de los
productos que se obtienen. Sin embargo, a pesar de ser un proceso ampliamente utilizado en la industria, su estudio y
entendimiento es muy limitado debido a la complejidad de los mecanismos de deformación elastoplástica que ocurren, así
como a la influencia de los diversos parámetros del proceso involucrados. Por tal motivo, la mayoría de las ecuaciones y
modelos reportadas en la literatura se enfocan en el embutido de formas cilíndricas. En este trabajo se presenta el análisis
y modelado numérica del proceso de embutido de lámina metálica de formas elípticas. El objetivo es derivar y evaluar
modelos que permitan determinar la altura y fuerza máxima de embutido de piezas con forma elíptica. Para lograr lo
anterior, se utilizan modelos teóricos existentes para el embutido cilíndrico y rectangular, así como modelos numéricos MEF
y datos experimentales disponibles en la literatura o industria.

Palabras Clave: embutido de lámina, forma elíptica, diámetro equivalente, altura de embutido, fuerza de embutido, Método del Elemento Finito
(MEF)

ABST RACT

The deep-drawing of sheet metal process is widely used in the automotive, aerospace, appliances and general metal-mechanic
industries due to the enhance properties of the products in terms of structural strength and weight. However, even though it
is a process widely used in the industry, its study and understanding is limited because of the complex elastoplastic
mechanisms that are present, and the influence of the several different process parameters involved. For this reason, most of
the equations and models available in the literature are focused in cylindrical shapes. This paper presents the analysis and
numerical modeling of the deep drawing process of elliptical sheet metal shapes. The aim is to derive and evaluate models
capable to determine the maximum force and maximum height of elliptical deep drawn parts. To accomplish this, existing
analytical models for deep drawing of cylindrical and rectangular parts, together with numerical MEF models and
experimental data from the literature, are used.

Keywords: Deep drawing of sheet metal, elliptical shapes, equivalent diameter, deep drawing height, deep drawing force, Finite Element Method
(FEM).

soldadura, además de que se pueden lograr piezas muy

1. Introducción complejas. Desde un punto de vista funcional, el
embutido profundo produce partes de alta resistencia y de
El proceso de embutido profundo de lámina metálica es poco peso, así como formas inalcanzables mediante otros
un proceso muy común para la fabricación de productos procesos de manufactura [1]. El proceso de embutido se
de lámina metálica que pueden ser ensamblados en realiza colocando una lámina de metal sobre la cavidad
futuras operaciones. Este proceso implica pocas de un dado y empujando el metal hacia la cavidad de éste
operaciones y reduce el desperdicio de material, lo cual con un punzón. Como resultado, la lámina se inserta
lo convierte en un proceso económico y rápido al dentro del dado generando una forma acopada. Las
eliminar operaciones costosas como el maquinado o aplicaciones del proceso de embutido de lámina metálica

ISSN 2448-5551 MM 30 Derechos Reservados © 2021, SOMIM

 MEMORIAS DEL XXVII CONGRESO INTERNACIONAL ANUAL DE LA SOMIM
22 al 24 DE SEPTIEMBRE DE 2021 PACHUCA, HIDALGO, MÉXICO

comprenden productos para la industria automotriz, punzón como el dado deben tener un radio en sus
aeronáutica, de electrodomésticos, y en general en toda esquinas para evitar el de corte de la lámina. Otro aspecto
la industria metalmecánica debido a las características importante es el claro que debe existir entre el dado y el
mejoradas en términos de resistencia/peso de los punzón, el cual debe estar en el rango 7% al 14% mayor
productos que se obtienen. al espesor de la lámina. Una medida de la severidad de la
operación de embutido es la relación de embutido, la cual
Desde los años 70´s se ha estudiado extensamente el se define para un embutido cilíndrico como la relación
comportamiento de los esfuerzos y deformaciones diámetro de la forma inicial entre el diámetro del punzón.
presentes en el proceso de embutido profundo mediante Cuanto mayor sea la relación de embutido mayor será la
las curvas de deformación característica obtenidas severidad de la operación. El límite superior se determina
durante la deformación de la forma inicial [2]. En la por la aparición de grietas en el material.
literatura existen diversos modelos para analizar el
proceso de embutido profundo convencional. Todos El embutido de formas no cilíndricas es cada vez más
estos modelos consideran mayormente dos regiones para frecuente en la industria. Sin embargo, la complejidad del
el modelado matemático del embutido profundo, el área mecanismo de deformación involucrado impide
de la forma inicial y la pared de la copa. La pared de la comúnmente el uso de las técnicas y procedimientos de
copa está sujeta a una fuerza de embutido la cual causa la diseño y manufactura inicialmente desarrollados para las
deformación de la forma inicial. Para formas iniciales de formas cilíndricas.
diámetros grandes, los esfuerzos presentes en la pared de
la copa pueden llegar a superar la resistencia última a la En [4] se desarrolló un nuevo proceso para incrementar
tracción del material, causando adelgazamiento de la la capacidad de embutido en copas de sección transversal
zona e incluso la fractura. La mayoría de la teoría sobre triangular, estrella y trébol. Como resultado se pudieron
el embutido profundo estudiada en la actualidad tiene obtener copas más altas embutiendo una lámina circular
fundamentos en la teoría de Siebel [3], quien estudió a a través de un dado cónico, el cual tenía en la cavidad del
los esfuerzos y las deformaciones presentes en el fondo una de las formas mencionadas anteriormente. En
embutido profundo de copas cilíndricas. Siebel [5] se realizaron experimentos de embutidos profundos
desarrolló expresiones analíticas para el cálculo de la hidrodinámicos de superficies irregulares con la finalidad
fuerza máxima de embutido para piezas cilíndricas, así de controlar el arrugado interno que se presenta durante
como expresiones para el cálculo de la altura máxima la operación. Un estudio similar se llevó a cabo en [6]
permisible que se puede lograr. donde se propuso un modelo matemático para el líquido
entre la materia prima inicial y el dado en un embutido
Debido a que en la literatura no existen modelos profundo hidrodinámico de forma irregular. Este modelo
analíticos capaces de predecir el formado de piezas se incorporó a un programa comercial de elemento finito,
elípticas, la industria recurre a métodos experimentales y al comparar con información experimental se obtuvo
basados en reglas empíricas o métodos de prueba y error, un método mediante el cual se podían determinar los
lo cual origina grandes tiempos y costos de desarrollo. parámetros principales de una operación de embutido
Por tanto, resulta necesario el desarrollo de modelos profundo hidrodinámico. Por otro lado, en [7] se
teóricos y/o numéricos que permitan analizar y mejorar investigó un método basado en el elemento finito para
el proceso de desarrollo de productos y procesos de determinar la forma óptima de la materia prima inicial
embutido de formas elípticas. Por lo anterior, el presente para piezas trapezoidales y formas de “L” en embutidos
trabajo está enfocado al análisis y simulación profundos. Se propuso la forma ideal de la pieza
computacional del proceso de embutido de lámina embutida y se analizó el flujo del material en reversa paso
metálica de formas elípticas. a paso para predecir la forma inicial de la materia prima.
La forma inicial propuesta de ambas figuras se comparó
contra casos experimentales.
2. Antecedentes
En [8] también se empleó el método de elemento finito
El embutido profundo es un proceso de formado de en conjunto con una red abductiva para predecir el
lámina metálica que se usa para hacer piezas de forma contorno óptimo de la forma inicial de un embutido
acopada, de caja y otras formas huecas más complejas. profundo de copa elíptica sin orejas. Se utilizó un código
Este proceso involucra formado por tensión y basado en el método de elemento finito para investigar
compresión. Durante el proceso de embutido la lámina se las características del flujo del material bajo diferentes
coloca sobre el dado y se mantiene en su lugar mediante parámetros del proceso, mientras que la red abductiva se
el pisador. Posteriormente el punzón baja y empuja la empleó para sintetizar los conjuntos de datos obtenidos
lámina dentro del dado, convirtiéndola así en una forma de simulaciones numéricas, obteniendo de este modo, un
acopada. modelo predictivo que permite obtener el contorno
En el embutido cilíndrico se embute una pieza de trabajo óptimo de la forma inicial para obtener una copa elíptica
dentro de un dado por medio de un punzón. Tanto el sin orejas. Un estudio similar se llevó a cabo en [9],

ISSN 2448-5551 MM 31 Derechos Reservados © 2021, SOMIM

 MEMORIAS DEL XXVII CONGRESO INTERNACIONAL ANUAL DE LA SOMIM
22 al 24 DE SEPTIEMBRE DE 2021 PACHUCA, HIDALGO, MÉXICO

donde se estudió la distribución de la fuerza del punzón

𝜇𝜋
en un embutido profundo de copa elíptica. En los 𝐷 2𝜇𝐹𝐵𝐻 𝑡
𝐹𝑑 = 𝜋𝑑𝑚 𝑡0 [𝑒 2 (1.1𝜎𝑓𝑚𝐼 (ln )+ ) + 𝜎𝑓𝑚𝐼𝐼 2𝑟0 ] (2)
𝑑𝑚 𝜋𝐷𝑡0 𝑑
experimentos realizados también se analizó la influencia
que tiene el radio del dado y el radio del punzón. Además,
donde 𝐹𝑑 es la fuerza de embutido para un diámetro
se determinó que la forma inicial es de suma importancia
para obtener una copa óptima. En [10] se utilizó el instantáneo de la materia prima 𝐷, 𝑑𝑚 es el diámetro
programa DEFORM-3D® para simular embutidos medio de la pared de la copa, el cual se determina como
profundos de copa elíptica de hoja metálica de aleación 𝑑𝑚 = 𝑑𝑝 + 𝑡0 , 𝑑𝑝 es el diámetro de punzón, 𝑡0 es el
de magnesio AZ31 a temperaturas elevadas. espesor inicial de la lámina, para un embutido cilíndrico,
𝜇 es el coeficiente de fricción entre la lámina y el radio
Por otro lado, en [11] se investigó la aparición de arrugas del dado, 𝐹𝐵𝐻 es la fuerza del pisador, 𝑟𝑑 es el radio del
en embutidos profundos de copa elíptica, cuya formación dado, 𝜎𝑓𝑚,𝐼 es el valor medio del esfuerzo de fluencia en
es severa en la región cercana al eje mayor. Se la zona del flanco desde el diámetro exterior hasta el
encontraron las primeras apariciones de arrugas cuando inicio del radio del dado, y 𝜎𝑓𝑚,𝐼𝐼 es el esfuerzo medio
la relación de aspecto a/b es de 1.33, donde a y b son los de fluencia en la zona del radio desde el inicio hasta el fin
semiejes mayor y menor de la forma inicial elíptica, de éste. Los valores de los esfuerzos medios en función
respectivamente. Conforme esta relación aumenta la de las deformaciones ideales del material, y tomando en
aparición de arrugas se retarda y la fractura se da antes en cuenta el endurecimiento por deformación, se determinan
la región cercana al eje mayor. En [12] se consideró el como [13]:
radio del dado y el radio del punzón como los parámetros
principales de diseño y se llevaron a cabo diversas 𝜎𝑓𝑚,𝐼 = 0.5𝐾(𝜑1𝑛 + 𝜑2𝑛 ) (3)
pruebas experimentales de embutido variando ambos
parámetros. 𝜎𝑓𝑚,𝐼𝐼 = 0.5𝐾(𝜑2𝑛 + 𝜑3𝑛 ) (4)

De esta revisión bibliográfica se observa que se han donde K y n son el coeficiente y exponente de
investigado diferentes métodos de optimización de la endurecimiento del material de trabajo, respectivamente,
formas iniciales elípticas apoyándose en modelos y las deformaciones ideales están dadas por las siguientes
numéricos basados en el método del elemento finito, así expresiones:
𝐷
como en pruebas experimentales. Se ha estudiado 𝜑1 = ln ( 𝐷𝑏) (5)
también el efecto de diversos parámetros geométricos
sobre la fuerza de embutido, así como los defectos que la √𝐷𝑏2 −𝐷 2 +(𝑑0 )2
variación de éstos pudiera ocasionar en la copa elíptica. 𝜑2 = ln ( ) (6)
𝑑0
Sin embargo, no se encontró en la literatura el desarrollo
𝑡0
de expresiones analíticas capaces de predecir la fuerza 𝜑3 = ln (1 + ) + 𝜑2 (7)
2𝑟𝑑 +𝑡0
necesaria para llevar a cabo un embutido elíptico con
base a los diferentes parámetros geométricos, del
En [3] se sugirió que la fuerza máxima de embutido se
material y del proceso que se encuentran presentes en la
presenta cuando el diámetro del blanco es 𝑑 ≈ 0.77𝐷𝑏 .
operación. Al mismo tiempo se notó la ausencia en la
Sustituyendo este valor en la ec. (2), se obtiene la fuerza
literatura de expresiones analíticas para el cálculo de la
máxima de embutido cilíndrico:
altura máxima permisible para un embutido profundo
elíptico. 𝜇𝜋 𝑑𝑓,𝑚𝑎𝑥
𝐹𝑑,𝑚𝑎𝑥 = 𝜋𝑑𝑚 𝑡0 [𝑒 2 (1.1𝜎𝑓𝑚𝐼 (ln 𝑑𝑚
)+
3. Análisis del embutido 2𝜇𝐹𝐵𝐻 𝑡0
) + 𝜎𝑓𝑚𝐼𝐼 ] (8)
𝜋𝑑𝑓,𝑚𝑎𝑥 𝑡0 2𝑟𝑑
3.1 Embutido cilíndrico
De acuerdo con Siebel [3] la fuerza total requerida para La ecuación anterior es la expresión reportada en [14]
una operación de embutido cilíndrico se puede para el cálculo de la fuerza máxima de embutido
determinar como: cilíndrico.
Por otro lado, la altura máxima permisible de la copa
𝐹𝑑 = 𝐹𝑖𝑑 + 𝐹𝑏 + 𝐹𝑓𝑓 + 𝐹𝑟𝑑 (1) embutida se determina principalmente por el diámetro de
la forma inicial. Se puede obtener una aproximación de
donde 𝐹𝑑 es la fuerza total de embutido, 𝐹𝑖𝑑 es la fuerza la altura final de la copa considerando que el volumen de
debida a la deformación ideal, 𝐹𝑏 es la fuerza de doblado la copa debe ser igual al volumen de la materia prima
y enderezado de la lámina, 𝐹𝑓𝑓 es la fuerza de fricción inicial, y que el espesor de la copa final es igual al espesor
entre el pisador y la lámina, y 𝐹𝑟𝑑 es la fuerza debido a la inicial de la materia prima [15]. Con base en lo anterior,
fricción en el radio del dado. Al sustituir los valores de la altura máxima de embutido, h, es [16]:
cada una de estas fuerzas en la ecuacion de Siebel, se
obtiene la siguiente ecuación:

ISSN 2448-5551 MM 32 Derechos Reservados © 2021, SOMIM

 MEMORIAS DEL XXVII CONGRESO INTERNACIONAL ANUAL DE LA SOMIM
22 al 24 DE SEPTIEMBRE DE 2021 PACHUCA, HIDALGO, MÉXICO

𝐷𝑏2 −𝑑𝑝
2
equivalente 𝐷𝑏,𝑒 se puede definir como el diámetro de
ℎ= 4𝑑𝑝
(9)
una forma inicial circular cuya área 𝐴𝐶 es igual al área de
la forma inicial elíptica 𝐴𝐸 , con un semieje mayor 𝐴 y un
Para que la ecuación anterior sea válida, debe existir una semieje menor 𝐵. Considerando 𝐴𝐶 = 𝐴𝐸 el diámetro
relación entre el diámetro de la pieza y su radio del fondo
equivalente 𝐷𝑏,𝑒 puede ser obtenido como:
mayor o igual a 20.

3.2 Embutido elíptico 𝐷𝑏,𝑒 = 2√𝐴𝐵 (10)

En general, los parámetros presentes en el proceso de
embutido cilíndrico también están presentes en el
embutido elíptico, sólo varían algunos de los parámetros
geométricos del herramental como se puede apreciar en
la Figura 1. En un embutido elíptico los diámetros de la
forma cilíndrica se sustituyen por los diámetros mayor y
menor de la forma elíptica. Además, se cuenta con la
relación de aspecto 𝛾𝑒𝑙 = 𝐷𝑏𝑐 ⁄𝐷𝑏𝑙 que indica que tan
larga es la elipse; por ejemplo, para un embutido a) b)
profundo de forma circular 𝐷𝑏𝑐 ⁄𝐷𝑏𝑙 = 1.
Figura 2 - Silueta inicial de material: a) elipse y b)
círculo.

El diámetro equivalente del punzón, 𝑑𝑝,𝑒 , se puede

determinar de la misma forma que para el diámetro
equivalente de la materia prima inicial, solo que en este
caso se igualan las áreas de las secciones transversales
del punzón circular 𝐴𝑃𝐶 y del punzón elipsoidal 𝐴𝑃𝐸 .

𝑑𝑝,𝑒 = 2√𝑎𝑏 (11)

A partir de la ec. (8) se puede obtener la fuerza máxima

de embutido elipsoidal al sustituir el perímetro medio
(𝜋𝑑𝑚 ) de una copa cilíndrica por el perímetro medio de
una copa elíptica (𝑃𝑚 ), el perímetro externo de la forma
inicial cuando la fuerza de embutido es máxima
(0.77𝜋𝐷𝑏 ) por 0.77𝑃𝑏 , donde 𝑃𝑏 es el perímetro de la
forma elíptica, y los diámetros cilíndricos 𝐷𝑏 y 𝑑𝑝 se
Figura 1 - Parámetros geométricos de un embutido sustituyen por diámetros equivalentes 𝐷𝑏,𝑒 y 𝑑𝑝,𝑒 . De esta
elíptico. manera se obtiene la siguiente ecuación para la fuerza
máxima de embutido elipsoidal:
La teoría del embutido cilíndrico puede ser utilizada para
otras formas geométricas si se utiliza el concepto de 𝜇𝜋 0.77𝐷𝑏,𝑒 2𝜇𝐹𝐵𝐻
“diámetro equivalente”. Este concepto ya ha sido 𝐹𝑑,𝑚𝑎𝑥 = 𝑃𝑚 𝑡0 [𝑒 2 (1.1𝜎𝑓𝑚𝐼 ln ( 𝑑 )+ )+
𝑝,𝑒+𝑡0 0.77𝑃𝑏 𝑡0
utilizado y evaluado en la literatura para embutidos de 𝑡
forma rectangular [14], [16], [17], [18]. Algunos trabajos 𝜎𝑓𝑚𝐼𝐼 2𝑟0 ] (12)
𝑑
han estudiado el embutido de formas no regulares
utilizando una relación de embutido equivalente, la cual
está basada en el área de la materia prima inicial, y el área donde el perímetro de la forma elíptica y el perímetro
transversal del punzón y la altura de la pieza, así como en medio de la copa elíptica pueden ser determinados con
el principio de conservación de volumen [17]. De base en la Figura 3, de la siguiente manera:
acuerdo con Lange [14], el diámetro equivalente 𝐷𝑏,𝑒 es
definido como el diámetro de una forma inicial circular 3𝑇𝑏
𝑃𝑏 = 𝜋(𝐴 + 𝐵) (1 + 10+ ) (13)
cuya área 𝐴𝐶 es igual al área de una forma inicial √4−3𝑇𝑏
rectangular 𝐴𝑅 , con un largo 𝐴 y ancho 𝐵.
𝑚 3𝑇
𝑃𝑚 = 𝜋(𝐿 + 𝐻) (1 + 10+√4−3𝑇 ) (14)
A partir del concepto de diámetro equivalente también se 𝑚

puede determinar el diámetro equivalente para otras

figuras geométricas. En el caso de la elipse, el diámetro

ISSN 2448-5551 MM 33 Derechos Reservados © 2021, SOMIM

 MEMORIAS DEL XXVII CONGRESO INTERNACIONAL ANUAL DE LA SOMIM
22 al 24 DE SEPTIEMBRE DE 2021 PACHUCA, HIDALGO, MÉXICO

𝐴−𝐵 2 𝐿−𝐻 2 𝑡0
donde 𝑇𝑏 = (𝐴+𝐵 ) , 𝑇𝑚 = (𝐿+𝐻) , 𝐿 = 𝑎 + y 𝐻=
2
𝑡0
𝑏+ .
2

Figura 3 - Perímetro de la sección transversal de una

copa elipsoidal.

Para determinar la altura máxima permisible a la cual se

puede embutir una pieza elíptica, se propone una
formulación geométrica usando el concepto del diámetro
equivalente y las ecuaciones para el embutido circular.

En la ec. (9) se pueden sustituir los diámetros de la

materia prima inicial 𝐷𝑏 y del punzón 𝑑𝑝 , por los
diámetros del blank y punzones equivalentes, ecs. (10) y
(11), obteniendo una expresión para el cálculo de la altura
máxima permisible para un embutido elíptico, ℎ𝑒 :
Figura 4 - Metodología para realizar el análisis MEF
2 2 del proceso de embutido.
𝐷𝑏,𝑒 −𝑑𝑝,𝑒
ℎ𝑒 = (15)
4𝑑𝑝,𝑒

Cabe mencionar que esta ecuación es válida cuando el

espesor de la lámina t es muy pequeño en relación con las
dimensiones de la pieza, satisfaciendo lo siguiente: t≪a
y t≪b.

4. Modelado numérico del embutido

La metodología utilizada en la simulación numérica MEF

del proceso de embutido de lámina, se muestra en la
Figura 4. Las geometrías de las herramientas del Figura 5 – Modelado geométrico del embutido.
embutido se modelaron en el programa Autodesk
Inventor®. Las piezas modeladas fueron el punzón, el Debido a que el embutido profundo consiste en la
pisador, la lámina y el dado, los cuales se modelaron deformación de una lámina delgada, se seleccionó un
como superficies, es decir, no poseen volumen. Se aplicó elemento tipo cascaron de 4 nodos y 12 grados de
simetría al modelo, por lo que solamente se modeló un libertad, el cual permite características no lineales para
cuarto de la geometría, como se muestra en la Figura 5. un análisis explícito dinámico. El tamaño de elemento
seleccionado fue de 1 mm, el cual se definió con base a
El proceso de formado de lámina metálica involucra un análisis convergencia, obteniendo una variación
grandes deformaciones en la pieza de trabajo y fricción máxima del 0.74% en la solución respecto al tamaño de
entre sus componentes, por lo que se requiere de un elemento más pequeño utilizado.
análisis transiente dinámico no lineal. Para ello se utilizó
el programa ANSYS LS-DYNA®, el cual utiliza el
método de “Newton-Raphson” para la solución de En términos de las condiciones de frontera, las partes
problemas no lineales. correspondientes al herramental se simularon como un
cuerpo rígido, es decir, no sufren deformaciones. Dichas

ISSN 2448-5551 MM 34 Derechos Reservados © 2021, SOMIM

 MEMORIAS DEL XXVII CONGRESO INTERNACIONAL ANUAL DE LA SOMIM
22 al 24 DE SEPTIEMBRE DE 2021 PACHUCA, HIDALGO, MÉXICO

partes están compuestas por el punzón, el pisador y el el momento inmediato anterior a cuando ocurre la falla
dado. El primero no puede rotar en ninguna dirección y (𝑓𝑎𝑙𝑙𝑎 > 1).
sus desplazamientos estás restringido en todas las
direcciones excepto en el eje Y. Las mismas condiciones En el proceso de embutido existe el contacto entre el
herramental y la pieza de trabajo, la cual al desplazarse
se aplican para el pisador. El dado por su parte se
hacia dentro del dado genera fricción en estas interfaces.
encuentra empotrado, es decir, está restringido en todos
Para definir éstas se deben determinar las superficies de
sus grados de libertad. Al pisador se le aplica una fuerza contacto, así como su tipo y las opciones de fricción en
constante en –Y. Al punzón se le agrega un el mismo. En este trabajo se utilizó el criterio de fricción
desplazamiento controlado en –Y. Por otro lado, la de Coulomb. Además, se seleccionó el tipo de contacto
lámina se define como un cuerpo deformable y se de superficie a superficie y su opción dirigida al formado
mantiene libre en todos sus grados de libertad. En la de metal, en la cual el herramental se define típicamente
Figura 6 se muestran las condiciones utilizadas. como la superficie objetivo (maestro) mientras que la
pieza de trabajo (lamina) se define como la superficie de
contacto (esclavo).

Para la validación del modelo MEF, se seleccionó un

caso experimental de un embutido rectangular presentado
en [19] y [20]. Se seleccionó este caso debido a que ha
sido utilizado por otros autores para validar sus modelos
MEF, tal como en [21] y [22]. Como parte de la
validación del modelo MEF se estudió la influencia que
tiene la densidad de la malla y la velocidad del punzón,
dado que son las principales variables de interés en una
simulación MEF explicita. Los resultados mostraron un
error promedio aproximado de 5.5% del modelo con
respecto a los resultados experimentales, y sugirieron un
tamaño de elemento de 1 mm (Figura 7) y una velocidad
de embutido de 250 mm/s (Figura 8).

Figura 6 – Modelado MEF del embutido.

El modelo de material utilizado para la lámina metálica

es el modelo de plasticidad de Ley de Potencia (Power
Law) definido por la siguiente expresión:
1
𝑒𝑓𝑓 𝑛
𝜎𝑦 = [1 + (𝐶̇ ) ] 𝐾[𝜀𝑒 + 𝜀𝑝 ]
𝜀 𝑃
(16)

𝑒𝑓𝑓
donde 𝜀𝑝 es la deformación plástica equivalente, 𝜀𝑒 es
la deformación elástica (cuando se alcanza la cedencia),
𝐾 y 𝑛 son el coeficiente y el exponente de
endurecimiento, respectivamente, y 𝐶 y 𝑃 son los Figura 7 – Análisis de convergencia (T.E. tamaño de
parámetros de la tasa de deformación de Cowper- elemento).
Symonds, a los cuales se les asigna un valor de 𝐶 =
10,000 y 𝑃 = 1 de manera que la sensibilidad a la
deformación se deprecie por ser un formado a
temperatura ambiente. El criterio de falla utilizado fue:
𝜎
𝑓𝑎𝑙𝑙𝑎 = 𝜎 𝑉𝑀 (17)
𝑈𝑇𝑆

donde 𝜎𝑈𝑇𝑆 es el la resistencia última a la tracción del

material y 𝜎𝑉𝑀 es el esfuerzo equivalente de von Mises.
La falla ocurre cuando 𝜎𝑉𝑀 > 𝜎𝑈𝑇𝑆 ; es decir, cuando
𝑓𝑎𝑙𝑙𝑎 > 1. La altura máxima del embutido se obtiene en
Figura 8 – Análisis de velocidad de embutido.

ISSN 2448-5551 MM 35 Derechos Reservados © 2021, SOMIM

 MEMORIAS DEL XXVII CONGRESO INTERNACIONAL ANUAL DE LA SOMIM
22 al 24 DE SEPTIEMBRE DE 2021 PACHUCA, HIDALGO, MÉXICO

importante mencionar que un error positivo indica una

sobre-estimación del valor, mientras que un valor
5. Resultados y discusión negativo indica una sub-estimación del valor. En el caso
de la fuerza de embutido es más conveniente una
Para analizar el desempeño de los modelos analíticos y sobreestimación del valor, ya que ello conducirá a una
numéricos para el embutido elíptico, se seleccionaron mejor selección de la prensa o a un mejor diseño del
cinco casos de estudio, los cuales se muestran en la Tabla herramental.
1. En esta tabla se incluyen los parámetros geométricos,
del material y del proceso para cada caso. Para cada caso En el caso de la altura máxima de embutido elíptico, los
de estudio se determinó la fuerza máxima y la altura resultados de la Tabla 2 muestran que el error promedio
máxima de embutido elíptico de manera analítica y de del modelo analítico, ec. (15), es de 41.24%, siendo
manera numérica. De manera analítica se utilizaron las considerablemente más alto en comparación con el error
ecuaciones (12) y (15) para determinar la fuerza máxima obtenido en el caso de la fuerza. Sin embargo, la
y la altura máxima de embutido, respectivamente. De desviación estándar es un poco más baja, 28.98%,
manera numérica se realizaron las simulaciones MEF de teniendo errores máximos de más del 80%. Este error tan
acuerdo al modelo y metodologías descritas en la sección alto puede conducir a una sobre estimación de la altura
anterior. Adicionalmente se calculó el error de los máxima de embutido, lo cual representa una desventaja
modelos analíticos con respecto a la simulación MEF. que conduce a un riesgo de falla de la pieza.

Los resultados del análisis del embutido elíptico se Finalmente se puede mencionar que la metodología de
muestran en la Tabla 2. Estos resultados muestran que el análisis para el proceso de embutido elíptico, así como
modelo analítico para el cálculo de la fuerza máxima de los modelos propuestos basados en el concepto de
embutido, ec. (12), tiene un error de predicción promedio diámetro y perímetro equivalente, demuestran en general
de 3.34 %, el cual es bajo. La desviación estándar un buen desempeño, y representan una alternativa para
también es baja, 6.71%, y muy aceptable. Por lo anterior reducir los métodos de prueba y error comúnmente
se puede decir que en general el cálculo de la fuerza utilizados en la industria. También son una alternativa
máxima de embutido elíptico mediante la ecuación (12) para mejorar la eficiencia del proceso de diseño y
tiene un error máximo aproximado del 10%. Es desarrollo de productos embutidos con forma elíptica.

Tabla 1 – Casos de estudio de embutidos elípticos.

Casos 𝑨 𝒂 𝒕 𝒄 𝒓𝒅 𝒓𝒑𝒃 𝑲 𝝈𝑼𝑻𝑺 𝝈𝒚 𝑭𝑩𝑯
Material 𝜸𝒆𝒍 𝒏 µ
[23] (𝐦𝐦) (𝐦𝐦) (𝐦𝐦) (𝐦𝐦) (𝐦𝐦) (𝐦𝐦) (𝐌𝐏𝐚) (𝐌𝐏𝐚) (𝐌𝐏𝐚) (𝐤𝐍)
Acero
1 50 30 0.63 0.8 1 5 5 1400 0.4 964.6 315 16 0.1
inoxidable
2 Acero 50 30 0.63 0.8 1 5 5 567.29 0.2637 437.4 164 16 0.1
3 Acero 50 30 0.63 0.8 1 5 0.41 567.29 0.2637 460 164 16 0.1
4 Acero 50 30 0.63 0.8 1 5 10 567.29 0.2637 460 164 16 0.1
5 Acero 50 30 0.7 0.8 1 5 5 567.29 0.2637 460 164 16 0.1

Tabla 2 – Resultados del análisis del embutido elíptico.

Fuerza máxima de embutido (N) Altura máxima de embutido (mm)
Casos [23]
𝑭𝐝 (𝐌𝐄𝐅) 𝑭𝐝 (𝐄𝐜. 𝟏𝟐) 𝐄𝐫𝐫𝐨𝐫 (%) 𝒉𝐝 (𝐌𝐄𝐅) 𝒉𝐝 (𝐞𝐜. 𝟏𝟓) 𝐄𝐫𝐫𝐨𝐫 (%)
1 39,100.00 42,836.48 9.56 17.70 21.05 18.93

2 21,420.00 22,222.30 3.75 14.09 21.05 49.40

3 22,760.00 22,222.30 -2.36 11.30 21.05 86.28

4 20,110.00 22,222.30 10.50 18.50 21.05 13.78

5 24,220.00 23,188.17 -4.26 16.19 22.31 37.80

Error de predicción (%) ±
3.3 ± 6.71 41.24 ± 28.98
variabilidad (%)

ISSN 2448-5551 MM 36 Derechos Reservados © 2021, SOMIM

 MEMORIAS DEL XXVII CONGRESO INTERNACIONAL ANUAL DE LA SOMIM
22 al 24 DE SEPTIEMBRE DE 2021 PACHUCA, HIDALGO, MÉXICO

processes by the finite-element method,” J. Mater.

6. Conclusiones Process. Technol., vol. 75, pp. 94–99, 1998.
[8] C. T. Lin and C. T. Kwan, “Application of abductive
En este trabajo se ha presentado un estudio del proceso network and FEM to predict the optimal blank
de embutido de lámina metálica de formas elípticas. Se contour of an elliptic cylindrical cup from deep
ha propuesto el uso del concepto de diámetro equivalente drawing,” J. Mater. Process. Technol., vol. 209, no. 3,
para el análisis de formas elípticas. Como resultado se pp. 1351–1361, 2009.
han propuesto dos modelos matemáticos para determinar [9] D. H. Park and P. K. D. V. Yarlagadda, “Effects of
la fuerza máxima de embutido y la altura máxima de punch load for elliptical deep drawing product of
embutido de formas elípticas. El modelo para el cálculo automotive parts,” Int. J. Adv. Manuf. Technol., vol.
de la fuerza máxima de embutido mostró una buena 35, no. 7–8, pp. 814–820, Jan. 2008.
precisión (error promedio del 3.34%) y ventajas que [10] T. S. Yang, “Finite element analysis of elliptic cup
favorecen su uso y aplicación en los procesos de diseño deep drawing of magnesium alloy sheet,” J. Achiev.
y manufactura de formas elípticas. Por otro lado, el Mater. Manuf. Eng., vol. 27, no. 2, pp. 139–142,
modelo para el cálculo de la altura máxima de embutido 2008.
resultó con un error promedio mayor (41.24%), lo cual [11] J. B. Kim, J. W. Yoon, D. Y. Yang, and F. Barlat,
representa una limitante en su aplicación práctica. Ambos “Investigation into wrinkling behavior in the elliptical
modelos tienen la ventaja de que pueden ser utilizados cup deep drawing process by finite element analysis
para determinar de forma eficiente la fuerza máxima y using bifurcation theory,” J. Mater. Process.
altura máxima de embutido de formas elípticas, Technol., vol. 111, no. 1–3, pp. 170–174, 2001.
reduciendo el uso de simulaciones numéricas o pruebas [12] D. H. Park, S. S. Kang, and S. B. Park, “A study on
experimentales durante el análisis y desarrollo de the improvement of formability for elliptical deep
procesos de embutido de este tipo de formas. drawing processes,” J. Mater. Process. Technol., vol.
113, no. 1–3, pp. 662–665, 2001.
Agradecimientos [13] A. Rivas, H. I. Medellín, D. F. de Lange, and P. de J.
García-Zugasti, “Performance evaluation of
El segundo autor agradece al CONACYT por la beca analytical expressions for cylindrical and rectangular
otorgada para la realización de estudios de Posgrado. deep drawing force estimation,” J. Manuf. Process.,
vol. 36, no. September 2017, pp. 340–350, 2018.
REFERENCIAS [14] L. Kurt, Handbook of Metal Forming -Society of
Manufacturing Engineers. New York: Mc. Graw Hill,
1985.
[15] Z. Marciniak, J. Duncan, and S. Hu, Mechanics of
[1] V. Boljanoric, Sheet metal forming and Die design. Sheet. Woburn: Butterworth-Heinemann, 2002.
New York: Industrial Press Inc, 2004. [16] A. Rivas, “Análisis del efecto de diversos parámetros
[2] M. Karali, “Examination of the strength and ductility en el desempeño del proceso de embutido de formas
of AA-1050 material shaped with the multistage deep rectangulares,” Tesis de maestría, Universidad
drawing method,” Arch. Metall. Mater., vol. 56, no. Autonoma de San Luis Potosí, 2015.
2, pp. 223–230, 2011. [17] H. I. Medellín-Castillo, P. De J. García-Zugasti, D. F.
[3] E. Siebel and H. Beisswanger, Deep Drawing (in De Lange, and F. J. Colorado-Alonso, “Analysis of
German). Munich: Carl Hanser, 1955. the allowable deep drawing height of rectangular steel
[4] M. A. Hassan, I. M. Hassab-Allah, L. M. A. Hezam, parts,” Int. J. Adv. Manuf. Technol., vol. 66, no. 1–4,
M. A. Mardi, and M. Hamdi, “Deep Drawing of pp. 371–380, 2013.
Asymmetric Cups through Conical Die without Blank [18] D. E, T. Mizuno, and Z. Li, “Stress analysis of
Holder,” Proc. World Congr. Eng., vol. II, pp. 1–7, rectangular cup drawing,” J. Mater. Process.
2015. Technol., vol. 205, no. 1–3, pp. 469–476, 2008.
[5] B. Meng, M. Wan, X. Wu, S. Yuan, X. Xu, and J. Liu, [19] Danckert, “Experimental investigation of a square-
“Inner wrinkling control in hydrodynamic deep cup deep-drawing process,” J. Mater. Process. Tech.,
drawing of an irregular surface part using vol. 50, no. 1–4, pp. 375–384, 1995.
drawbeads,” Chinese J. Aeronaut., vol. 27, no. 3, pp. [20] A. Makinouchi, E. Nakamachi, E. Onate, and R.
697–707, 2014. Wagoner, “Proceedings of the Numisheet 1993,”
[6] L. H. Lang, J. Danckert, K. B. Nielsen, D. C. Kang, Second Int. Conf. Numer. Simul. 3D Sheet Met.
and S. H. Zhang, “Key technologies of the simulation Form. Process.
of the hydrodynamic deep drawing of irregular parts,” [21] M. Khelifa and M. Oudjene, “Numerical damage
J. Mater. Process. Technol., vol. 150, no. 1–2, pp. 40– prediction in deep-drawing of sheet metals,” J. Mater.
47, 2004. Process. Technol., vol. 200, no. 1–3, pp. 71–76, 2008.
[7] S. Kim, M. Park, S. Kim, and D. Seo, “Blank design [22] H. W. Zhang, W. L. Xu, S. L. Di, and P. F. Thomson,
and formability for non-circular deep drawing “Quadratic programming method in numerical

ISSN 2448-5551 MM 37 Derechos Reservados © 2021, SOMIM

 MEMORIAS DEL XXVII CONGRESO INTERNACIONAL ANUAL DE LA SOMIM
22 al 24 DE SEPTIEMBRE DE 2021 PACHUCA, HIDALGO, MÉXICO

simulation of metal forming process,” vol. 191, no.

2002, pp. 5555–5578, 2010.
[23] Y.-M. Huang, “Analysis of Elliptical Cup Drawing
Process of SUS304 Stainless Metal,” vol. 86, pp.
350–357, 2010.

ISSN 2448-5551 MM 38 Derechos Reservados © 2021, SOMIM