MEMORIAS DEL XIX CONGRESO INTERNACIONAL ANUAL DE LA SOMIM

25 al 27 DE SEPTIEMBRE, 2013 PACHUCA, HIDALGO, MÉXICO

PROCEDIMIENTO NUMÉRICO PARA EL MODELADO DE MATERIALES

HIPERELÁSTICOS SOMETIDOS A TENSIÓN
1
Estrada Barbosa Quirino, 1Szwedowicz Wasik Dariusz, 2Majewski Szymiec Tadeusz, 3Bedolla Hernández Jorge
Cortés García Claudia1, Martínez Rayón Eladio1
1
Departamento de Ingeniería Mecánica, Centro Nacional de Investigación y Desarrollo Tecnológico (CENIDET),
Prolongación Av. Palmira esquina Apatzingán, Cuernavaca, Morelos, México.
Teléfono (777) 362 7770
2
Departamento de Ingeniería Industrial y Mecánica, Universidad de las Américas-Puebla
Ex Hacienda Santa Catarina Mártir S/N, San Andrés Cholula, Puebla, México.
Teléfono (222) 229 2673
3
Departamento Metal-Mecánica, Instituto Tecnológico de Apizaco
Av. Instituto Tecnológico S/N, Apizaco, Tlaxcala, México.
Teléfono (241) 417 2010

[email protected], [email protected], [email protected], [email protected],

[email protected], [email protected]

RESUMEN. S Medida de esfuerzo (Piola –Kirchoff)

E Medida de deformación (Green)
En este trabajo se presenta la elaboración de un I1, I2, I3 Invariantes del tensor de deformación
procedimiento y desarrollo de una metodología λ1, λ2, λ3 Alargamientos principales
para el modelado numérico con propiedades N Orden de la función de energía.
hiperelasticas. La corroboración de la Cijk Constantes del material (Experimental)
metodología se realizó mediante la simulación Jel Es la razón de cambio de volumen
de un ensayo a tensión en el software comercial Di Constantes del material (Experimental)
de elementos finitos Abaqus. La geometría de la Wdev Componente deviatórica de la función
probeta y condiciones de la frontera del modelo de energía
discreto se establecieron de acuerdo a la norma Wvol Componente volumétrica de la función
ASTM D-638, referente a pruebas de plásticos de energía
sometidos a tensión. En base al presente estudio αi, μi Constantes del material
se encontró que el modelo constitutivo de σnom Esfuerzo nominal
Marlow presenta una mejor aproximación del Fx Fuerza en dirección x
comportamiento hiperelástico de gomas al ser A0 Área inicial
sometidas a cargas de tensión. ɛnom Deformación nominal
x1 Desplazamiento del punto 1
ABSTRACT x2 Desplazamiento del punto 2
L0 Longitud inicial entre el punto 1 y 2
A procedure and development of a methodology
to numerical modeling with hyperelastics INTRODUCCIÓN
properties is presented in this resource. The
methodology was corroborated by a tensile test Actualmente en diversas áreas de ingeniería está
simulation carried out in finite element software, en ascenso la aplicación de elastómeros,
ABAQUS. The specimen dimensions and particularmente en el control de cargas de
boundary conditions of discrete model were impacto y vibraciones mecánicas. La capacidad
established in accord to Standard Test Method por parte de un elastómero de asumir grandes
for Tensile Properties of Plastics, ASTM D-638. deformaciones y posterior retorno a su estado
According to this work the Marlow model has a inicial al ser descargados, se denomina
better approximation of hyperelastic behavior of hiperelasticidad [1]. Tal característica permite su
rubber under tension loads. utilización en diversas aplicaciones industriales.
A diferencia de otros materiales cuyo
NOMENCLATURA comportamiento mecánico puede ser descrito a
partir de las constantes elásticas del material
W(E) Energía de deformación (módulo de Young y coeficiente Poisson), las

ISBN 978-607-95309-9-0 Página | 1079 Derechos Reservados © 2013, SOMIM

 MEMORIAS DEL XIX CONGRESO INTERNACIONAL ANUAL DE LA SOMIM
25 al 27 DE SEPTIEMBRE, 2013 PACHUCA, HIDALGO, MÉXICO

gomas presentan otras características no lineales

como viscoelasticidad y ablandamiento del
material [2]. Actualmente el Método de
Elementos Finitos (MEF) se posiciona como una
principal herramienta en el proceso de análisis
en el diseño de productos basados en materiales De esta forma la ecuación 2 se reescribe en la
elastoméricos. La mayor parte de los paquetes siguiente forma:
comerciales del MEF involucran diversos
modelos constitutivos para describir el
comportamiento hiperelástico. En general los
modelos constitutivos están divididos en dos
ramas, los primeros son modelos
fenomenológicos basados en función de la Considerando que la goma es cercana a la in-
densidad de energía de deformación y los compresibilidad, se obtiene la forma general del
segundos modelos denominados modelo de Rivlin para el caso de incompresibili-
microestructurales. El desarrollo del dad. [4]
procedimiento presentado advierte las
consideraciones que deben ser tomadas en cuenta
en la simulación numérica de materiales
hiperelásticos, con el objeto de obtener
resultados cercanos y congruentes a la realidad. Las constantes se obtienen a partir de ensayos
El presente estudio se limitó a la experimentales. A partir de la ecuación 6 se han
implementación y evaluación de modelos desarrollado diversos modelos constitutivos
fenomenológicos. basados en la primera y segunda invariante de
deformación. A continuación se describen en su
forma general algunos modelos básicos
MODELOS CONSTITUTIVOS empleados en el desarrollo de este trabajo.

Los modelos constitutivos permiten describir el

Forma Polinomial
comportamiento elástico no lineal de los
elastómeros. En el desarrollo de tales modelos se
supone un comportamiento elástico, isotrópico La forma polinomial es un modelo basado en la
además de incompresibilidad del material. La primera y segunda invariante del tensor deviató-
hiperelasticidad se puede expresar en términos rico de Cauchi-Green, y ha servido de base para
de la energía potencial de deformación W (E), la el desarrollo de nuevos modelos constitutivos, se
cual describe la cantidad de energía de define a partir de la siguiente expresión. [5]
deformación almacenada por unidad de volumen
bajo un estado de deformación dado. [3]

La ecuación 8 es aplicable para pequeñas defor-

maciones y deformaciones nominales menores al
100%.

De acuerdo a la teoría de Rivlin [4] la función de Modelo de Marlow

densidad de energía de deformación se puede
expresar en términos de las invariantes de de- Modelo constitutivo dependiente de la primera
formación de tal forma que: invariante de deformación, se emplea cuando los
datos experimentales disponibles son limitados.
Este modelo reproduce exactamente los datos
experimentales para un modo de deformación en
específico. [5]
Adicionalmente las invariantes de deformación
generalmente se expresan en función de los
alargamientos principales [2].

ISBN 978-607-95309-9-0 Página | 1080 Derechos Reservados © 2013, SOMIM

 MEMORIAS DEL XIX CONGRESO INTERNACIONAL ANUAL DE LA SOMIM
25 al 27 DE SEPTIEMBRE, 2013 PACHUCA, HIDALGO, MÉXICO

Modelo de Ogden la Figura 2 para ser implementado en el software

ABAQUS. La metodología se encuentra divida
Modelo propuesto por Ogden [6], utiliza los en tres etapas:
alargamientos principales en lugar de las inva-
riantes de deformación. 1. Análisis del modelo real e identificación de
condiciones de la frontera.

2. Selección del modelo de hiperelasticidad y

evaluación automática del material.
Este modelo presenta una buena respuesta fren- 3. Realización del modelo discreto y posterior
te a grandes deformaciones. Una de sus restric- validación con datos experimentales.
ciones es que no se debe utilizar cuando los datos
experimentales son limitados [5].

ENSAYOS EXPERIMENTALES MODELO REAL 1

Uniaxial, Planar, (ANÁLISIS)
volumétrico
Las constantes del material presentes en las
ecuaciones de los modelos hiperelásticos son Ensayo experimental Selección del Modelo Hiperelástico Modo de
deformación
obtenidas a partir de ensayos experimentales. (Marlow, Ogden, Polinomial)

Comúnmente suelen emplearse ensayos uniaxia- Modelo Hiperelástico

les de tensión, uniaxiales de compresión, tensión Evaluación Automática del
Mínimos cuadrados Estabilidad del
plana y equibiaxiales de tensión. El método de material Ajuste de curva
material
obtención más común de las constantes es a
partir del ajuste de las curvas experimentales por
No
el método de mínimos cuadros. La Fig. 1 muestra Diferencia 5%
los ensayos generalmente aplicables a gomas.
Si
Propiedades Mallado y tipo
hiperelásticas 2
Construcción del modelo de elementos
numérico por elementos finitos
Condiciones de Interacciones
frontera

Análisis no lineal
NGEOM activado

No
¿Converge?

Si

Análisis Exportar datos

Figura 1. Ejemplos de la respuesta de las gomas frente a Esfuerzo-Deformación

diversos ensayos [7] Fuerzas de contacto

Datos Experimentales Diferencias

2
(modelo real) Mínimas
No 1
PROCEDIMIENTO NUMÉRICO Si

La hiperelasticidad, es una de las propiedades Fin

más importantes en los elastómeros que debe ser

tomada en cuenta al realizar modelos numéricos.
La caracterización adecuada de la hiperelastici-
Figura 2. Procedimiento propuesto para simulación de
dad en los paquetes comerciales de elementos materiales hiperelásticos.
finitos requiere de una metodología que permita
relacionar los modelos constitutivos con los
datos experimentales disponibles. Tomando en Para poder evaluar el procedimiento se des-
cuenta lo anterior y para el caso específico se arrolló un análisis por el MEF que involucró un
presenta un procedimiento el cual se muestra en ensayo a tensión de goma vulcanizada estandari-

ISBN 978-607-95309-9-0 Página | 1081 Derechos Reservados © 2013, SOMIM

 MEMORIAS DEL XIX CONGRESO INTERNACIONAL ANUAL DE LA SOMIM
25 al 27 DE SEPTIEMBRE, 2013 PACHUCA, HIDALGO, MÉXICO

zado por la norma ASTM D-638. Tomando en tener mejor respuesta que los anteriores. Finali-
consideración que la carga aplicada al elastómero zada esta etapa se continuó con la realización del
es de tensión, la caracterización hiperelástica del modelo discreto del ensayo a tensión simple.
material puede realizarse a partir de un ensayo
uniaxial de tensión. La figura 3 muestra la rela- La geometría de la probeta así como las condi-
ción esfuerzo-deformación obtenida experimen- ciones y solicitaciones del ensayo fueron toma-
talmente por Treloar al evaluar dicho elastómero das de la norma ASTM D-638 referente a ensa-
[8]. yos de tensión en plásticos. La Figura 5 muestra
la geometría del espécimen experimental. [9]

Figura 5. Geometría del espécimen ASTM D638 Tipo I

(unidades en mm) [9]
Figura 3. Caracterización de la goma vulcanizada [8]

Para poder reproducir acertadamente el ensayo a

De acuerdo al procedimiento, el modelo hiper- tensión, se determinaron las condiciones de
elástico a elegir debe contemplar grandes defor- frontera y parámetros propios del ensayo tal
maciones alrededor de 4.37 (437% de deforma- como la velocidad de carga. El modelo discreto-
ción nominal). Además es imprescindible tomar se realizó con base 2D y elementos continuos.
en consideración la información experimental
limitada. Para fines prácticos de evaluación au- Una condición de encastre se confirió en los
tomática del material se eligieron cuatro mode- nodos del extremo izquierdo, la velocidad del
los constitutivos (Mooney- Rivlin, Polinomial de ensayo se aplicó a los nodos del extremo dere-
segundo orden, Marlow y Ogden de primer de cho. La deformación del espécimen se definió en
orden). En la figura 4 se presenta la comparación base a dos nodos de referencia x1 y x2 como se
los resultados obtenidos. muestra en la Figura 6.

Encastre

X1 X2

Velocidad

Figura 6. Condiciones de frontera

Figura 4. Evaluación automática de modelos hiperelásticos Durante el análisis de hiperelasticidad se utili-

zaron elementos del tipo CPS4R con integración
reducida. La figura 7 presenta el modelo discreto
Se observa en esta etapa que los modelos de de la probeta.
Mooney Rivlin y Ogden de primer orden son los
que tienen menor aproximación respecto a los
datos experimentales de Treloar [8], el primero a
consecuencia de ser de un orden muy bajo y el
segundo por poseer datos experimentales limita- igu
dos. Por otra parte el modelo de Marlow así
como el Polinomial de segundo orden mostraron Figura 7. Mallado del modelo

ISBN 978-607-95309-9-0 Página | 1082 Derechos Reservados © 2013, SOMIM

 MEMORIAS DEL XIX CONGRESO INTERNACIONAL ANUAL DE LA SOMIM
25 al 27 DE SEPTIEMBRE, 2013 PACHUCA, HIDALGO, MÉXICO

El estado de deformación obtenido se muestra en evaluados. En el modelo de Mooney Rivlin y el

la figura 8, donde se puede observar la propiedad polinomial de segundo orden se advierte la poca
hiperelástica alcanzada. efectividad de estos modelos para grandes de-
formaciones, haciéndose más visible en el pri-
mero, por ser un modelo constitutivo de primer
orden. Respecto a deformaciones menores al
a) 75 % se precisa una respuesta aceptable de am-
bos modelos con un 6% de diferencia respecto al
b) ensayo experimental de Treloar [8]

c)

Figura 8. Deformación del espécimen ASTM D638 Tipo I,

donde: a) Estado inicial y final de deformación, b)
Deformación extremo izquierdo y c) Deformación en el
extremo derecho.

Figura 10. Respuesta de modelos constitutivos I

Los esfuerzos se determinaron a partir de la suma
de las fuerzas de reacción originadas a lo largo La respuesta del modelo de Ogden y de Marlow
del extremo izquierdo de la probeta, tal y como se aprecia en la figura 11, donde el primero de
se puede observar en la figura 9. ellos presenta un comportamiento poco acertado
para deformaciones mayores al 50%, conse-
cuencia de una limitada información experimen-
tal durante la calibración del modelo.
∑ Fx
En contraparte el modelo de Marlow presenta
buena aproximación para pequeñas y grandes
deformaciones variando la diferencia con la parte
experimental en un 2 % para pequeñas deforma-
ciones y en un 5% para deformaciones mayores
al 100% de deformación.

Figura 9. Fuerzas de reacción

Los valores de σnom y ɛnom se definieron de la

siguiente manera [10]:

Figura 11. Respuesta de modelos constitutivos II

De esta forma se concluye el procedimiento

En la figura 10 se presenta la comparación de los presentado (ver figura 2), con un resultado
datos experimentales con los modelos hipere- satisfactorio con diferencias menores al 5%
lasticos aplicados al modelo discreto. A pesar de respecto al desarrollo experimental al utilizar el
que el procedimiento indicaba solo utilizar el modelo de Marlow para describir el
modelo de Marlow y Polinomial por tener mejor comportamiento hiperelástico. Esta reducida
respuesta durante la evaluación automática del diferencia entre análisis se debe a que el modelo
material, también se presenta otros modelos constitutivo de Marlow solo depende la primera

ISBN 978-607-95309-9-0 Página | 1083 Derechos Reservados © 2013, SOMIM

 MEMORIAS DEL XIX CONGRESO INTERNACIONAL ANUAL DE LA SOMIM
25 al 27 DE SEPTIEMBRE, 2013 PACHUCA, HIDALGO, MÉXICO

invariante de deformación I1 y no requiere de (3) Jarosław, R., Marlena, R., 2010, “The
múltiples ensayos experimentales para la Methodology of the Analysis of Elastomer
calibración del modelo. Bearings Properties”, Slovaquia
(4) Rivlin R.S., 1948, “ Large Deformations of
isotropic materials, IV. Further
CONCLUSIONES
Develompments of the general theory”. Phil.
Trans. R. Soc. United States.
1. La obtención de resultados próximos a los (5) Hibbit, Karlson & Sorensen, 2005,
experimentales requiere de un análisis previo de “Abaqus/standar User’s Manual Vol. II”,
las condiciones a las cuales estará sometido el United State.
modelo real. Tal análisis permite la selección
(6) Ogden, R.W., 1972. “Large deformation
correcta del ensayo de calibración para el modelo
isotropic elasticity-on the correlation of
hiperelástico (ensayos uniaxiales de tensión,
theory and experiment for incompressible
uniaxiales de compresión, tensión plana, equi-
biaxiales de tensión o volumétricos). rubberlike solids”. Proc. R. Soc. Lond. A,
United States.
3. De acuerdo a los resultados presentados en la
(7) García, M. Ruiz, C., 2005, “Hyperelastic
figura 4, modelos de primer o segundo orden
Material Modeling”, Reporte Técnico, Uni-
tales como el modelo polinomial u Ogden no
versidad EAFIT, Colombia.
aseguran una mejor respuesta hiperelástica del
modelo para deformaciones superiores al 50%, (8) Treloar, L., R., 1949. “ The physics of
incluso puede originar tiempos de procesamiento rubber elasticity”, Oxford Univ. Press, New
computacional más prolongados. York.

4. Para modelos discretos sujetos a esfuerzos (9) ASTM D 638-99, 1999, “Standard Test
combinados, es necesario proveer la mayor Method for Tensile Strength of Plastics”,
cantidad de ensayos experimentales (tensión, Annual Book of ASTM Standards, ASTM
compresión, cortantes) en el proceso de calibra- International, West Conshohocken, PA.
ción del modelo hiperelástico.
(10) Mendizabal, A., 2010, “Finite Element
5. Para deformaciones menores o iguales al 50% Simulation:Tensile test of rib cortical
se obtuvo una buena aproximación (diferencias bone”, Chalmers University of Technology,
de 3 al 5% respecto a los datos experimentales) Master Thesis, Sweden.
por parte de todos los modelos hiperelásticos
evaluados. Para deformaciones mayores al 50 %
y menores a 75% el modelo de Ogden empieza
a diferir en gran porcentaje con los valores expe-
rimentales (diferencias de superiores al 7%). De
esta forma al alcanzar deformaciones nominales
iguales o mayores al 100% solo el modelo de
Marlow representa con gran exactitud los datos
experimentales obtenidos por Treloar [8] (dife-
rencias menores al 5%).

REFERENCIAS
(1) Pascon, Coda, H. b., 2009. “Modelos
constitutivos para materiales hiperelásticos:
Estudio e implementación computacional.”
Cuadernos de Ingeniería de Estructuras. Sao
Pablo, Brasil
(2) Austrell, 1997, “Modeling of elasticity and
damping for filled elastomers”, PhD Thesis,
Lund University, Sweeden.

ISBN 978-607-95309-9-0 Página | 1084 Derechos Reservados © 2013, SOMIM