T6-viscoeslasticidad

Vista previa del material en texto

Caṕıtulo 6
Viscoelasticidad
La caracteŕıstica definitoria de los materiales elásticos es que el estado
tensional en un punto e instante depende exclusivamente de la deformación
en dicho punto, e instante. Esta es una suposición muy restrictiva y, aunque
suficientemente aproximada para gran parte de los análisis en ingenieŕıa
mecánica, civil y aeronáutica, no existe ningún material que sea elástico
en todo rango de deformación o de velocidad de deformación. En general,
tanto la historia del material (es decir, el valor de la deformación pasada
del punto) como la tasa de deformación, modifican el valor instantáneo de
la tensión. La incorporación de estos efectos complica las leyes constitutivas
del material pero permiten modelar con más precisión los materiales reales.
Existe una jerarqúıa de modelos materiales que, a base de modelar más
efectos en las ecuaciones representan de manera más fiel la respuesta de los
sólidos reales. La viscoelasticidad supone una elaboración de la respuesta
elástica que incorpora los efectos de dependencia de la velocidad de defor-
mación e historia. Estos efectos son imprescindibles para poder modelar sóli-
dos que poseen fluencia y relajación, dos comportamientos fundamentales
en los poĺımeros, los suelos e incluso el hormigón. Dentro de la complejidad
de estos modelos, nos centraremos en este caṕıtulo en la descripción de la
viscoelasticidad lineal que, como se explicará, no se refiere a una relación
lineal entre tensión y deformación como en el caso elástico.
La forma más sencilla de abordar la viscoelasticidad es mediante mo-
delos reológicos. Este tipo de idealización extiende el concepto del resorte
y permite aproximarse de manera sencilla e intuitiva al comportamiento
viscoelástico tensorial. De hecho, como se verá en este caṕıtulo, los modelos
viscoelásticos tensoriales se basan en una extensión de los modelos reológicos
al ámbito tridimensional.
Antes de comenzar el estudio de la respuesta viscoelástica es necesario
mencionar el papel fundamental que juega la temperatura como modula-
dor de la respuesta, especialmente en los poĺımeros. Estos materiales tienen
una temperatura, la llamada temperatura de transición v́ıtrea, por debajo de
121
122 Caṕıtulo 6. Viscoelasticidad
la cual se comportan de manera frágil y no muestran ninguno de los com-
portamientos caracteŕısticos de los materiales viscoelásticos. Por encima de
esta temperatura, la respuesta reológica aparece y además es muy sensible
al valor de la temperatura. Existe una temperatura de fusión donde ya el
material deja de ser sólido y cuya respuesta no estudiaremos aqúı.
Para exposiciones más completas de la teoŕıa de la viscoelasticidad se
puede consultar el texto clásico [2] o los más recientes [? 4].
6.1. Modelos reológicos
Existen dos fenómenos, la fluencia y la relajación, que son fundamentales
en el comportamiento de los sólidos, y que no pueden modelarse con leyes
constitutivas elásticas. De hecho, la motivación primera para el desarrollo
de la viscoelasticidad es la formulación de modelos que puedan reproducir
estos dos procesos.
Para acercarnos a estos modelos empleamos los llamadosmodelos reológi-
cos, que son sistemas mecánicos elementales que capturan de forma sencilla
los distintos tipos de comportamientos, a partir de una combinación de re-
sortes y amortiguadores.
� �
E
"
� �
"
⌘
Figura 6.1: Modelos simplificados de resorte y amortiguador.
En la Figura 6.1 se muestran los dos elementos básicos que empleare-
mos para describir la viscoelasticidad. El muelle o resorte es el elemento
elástico básico. Cuando se somete a una tensión � sufre una deformación "
cuya valor es
" = �
E
, (6.1)
siendo E la constante de rigidez del resorte. Por su lado, el amortiguador
es un elemento cuya tensión es proporcional a la velocidad de deformación
y la relación es:
� = ⌘"̇ , (6.2)
siendo ˙( ) la notación que emplearemos para indicar la derivada con respecto
al tiempo.
Caṕıtulo 6. Viscoelasticidad 123
6.1.1. Fluencia
Cuando un material sólido viscoelástico se somete a un estado tensional,
su deformación no permanece constante sino que cambia con el tiempo. De
hecho, el material parece que fluyera, pues la deformación aumenta progre-
sivamente como si fuera un fluido. Para describir el fenómeno de la fluencia
empleamos un modelo reológico compuesto por un resorte y un amortigua-
dor colocados en paralelo, y que se conoce con el nombre del modelo de
Kelvin o Kelvin-Voigt . Véase la Figura 6.2
� �
"
⌘
E
Figura 6.2: Modelo reológico de Kelvin.
Cuando un elemento de Kelvin se somete a una tensión �(t) = �̄, ésta se
reparte entre el resorte y el amortiguador de forma que se verifica en todo
instante
�̄ = E"(t) + ⌘"̇(t) . (6.3)
Si además sabemos que la deformación del elemento en el instante t = 0 es
nula, podemos integrar la ecuación diferencial anterior y obtener la defor-
mación en todo instante:
"(t) = �̄
E
�1 − e−E⌘ t� . (6.4)
Definiendo el tiempo de relajación del elemento como ⌧ = ⌘�E, entonces
se puede escribir de forma alternativa
"(t) = �̄
E
�1 − e−t�⌧� . (6.5)
Esta ecuación expresa que en el tiempo inicial la deformación es nula y
que ésta aumenta monótonamente hasta alcanzar un valor asintótico "∞ =
�̄�E. En este momento toda la tensión la soporta el resorte y el amortiguador
permanece descargado pues la velocidad de deformación tiene a cero. Véase
en la figura 6.3 su representación gráfica.
124 Caṕıtulo 6. Viscoelasticidad
0
0.2
0.4
0.6
0.8
1
1.2
0 2 4 6 8 10
"
�̄
/
E
t/⌧
0
0.2
0.4
0.6
0.8
1
1.2
0.1 1 10
"
�̄
/
E
log(t/⌧)
Figura 6.3: Fluencia del elemento de Kelvin. A la izquierda en escala normal
y a la derecha, semilogaŕıtmica. Nótese cómo el punto de inflexión en la
curva de la derecha ocurre en t�⌧ = e.
El cociente F (t) = "(t)��̄ se conoce con el nombre delmódulo de fluen-
cia y para el modelo de Kelvin es simplemente
F (t) = 1 − e−t�⌧
E
. (6.6)
6.1.2. Relajación
Otro fenómeno propio de los materiales viscoelásticos es la relajación ,
que consiste en la “disminución” del estado tensional cuando un sólido se
somete a una deformación dada. Este comportamiento, dual en cierto sentido
al de la fluencia, se aclara al estudiar el elemento de Maxwell , que combina
un resorte y un amortiguador en serie.
� �
"
E ⌘
Figura 6.4: Modelo reológico de Maxwell.
El modelo de Maxwell se caracteriza porque cuando una tensión � se
aplica sobre el mismo, ésta la recibe tanto el resorte como el amortiguador.
Por otro lado, la deformación " del conjunto resulta de las contribuciones de
ambos modelos elementales y por tanto se puede escribir:
"̇(t) = �̇(t)
E
+ �(t)
⌘
. (6.7)
Para calcular la relajación del elemento de Maxwell suponemos que se im-
pone una deformación "(t) = "̄ constante sobre el elemento y calculamos el
valor de la tensión en el tiempo resolviendo la ecuación diferencial lineal (6.7)
Caṕıtulo 6. Viscoelasticidad 125
0
0.2
0.4
0.6
0.8
1
1.2
0 2 4 6 8 10
� "̄
E
t/⌧
0
0.2
0.4
0.6
0.8
1
1.2
0.1 1 10
� "̄
E
log(t/⌧)
Figura 6.5: Relajación del elemento de Maxwell. A la izquierda en escala
normal y a la derecha, semilogaŕıtmica. Nótese cómo el punto de inflexión
en la curva de la derecha ocurre en t�⌧ = e.
empleando la condición inicial �(0) = E"̄. La solución de ésta es:
�(t) = E "̄ e−t�⌧ (6.8)
La relación �(t)�"(t) se conoce como el módulo de relajación
R(t) = E e−t�⌧ . (6.9)
6.1.3. El sólido lineal estándar
Como acabamos de ver el modelo de Kelvin experimenta una respuesta
al someterlo a una tensión constante que se puede identificar con la fluencia.
Por su parte, el elemento de Maxwell exhibe relajación de tensiones al so-
meterlo a un campo de deformaciones constante. Sin embargo, ninguno de
los dos modelos es capaz de representar ambos fenómenos y por tanto, para
acercarnos más al estudio del comportamiento viscoelástico debemos em-
plear un modelo reológicoalgo más complejo. El modelo del sólido lineal
estándar combina un elemento de Maxwell en paralelo con un resorte, co-
mo se indica en la Figura 6.6 (en algunos trabajos se define un sólido lineal
distinto con un elemento de Kelvin en serie con un resorte).
Repartiendo la tensión � entre las dos ramas del elemento y la defor-
mación " entre los dos componentes del elemento de Maxwell se obtiene la
ecuación diferencial que describe el comportamiento del sólido estándar:
�̇(t) + E1
⌘
�(t) = (E
1
+E∞)"̇(t) + E1E∞
⌘
"(t) (6.10)
Si definimos el tiempo de relajación ⌧ = ⌘�E
1
entonces la relación anterior
también se puede expresar como
�̇(t) + �(t)
⌧
= (E
1
+E∞)"̇(t) + E∞
⌧
"(t) (6.11)
A partir de la ecuación diferencial anterior podemos calcular los módulos
de fluencia y relajación del elemento estándar. Si la tensión aplicada sobre
126 Caṕıtulo 6. Viscoelasticidad
E1 ⌘
�
E1
� �
"
Figura 6.6: El modelo del sólido lineal estándar.
el elemento estándar es �(t) = �̄ constante y la deformación inicial "(0) =
�̄�(E
1
+E∞) encontramos que
"(t) = �̄
E∞ �1 − E1E1 +E∞ e− E∞E1+E∞ t⌧ � (6.12)
y por lo tanto el módulo de fluencia F (t) = "(t)��̄ será
F (t) = 1
E∞ �1 − E1E1 +E∞ e− E∞E1+E∞ t⌧ � (6.13)
Como en el caso del elemento de Kelvin, podemos obtener el módulo de
relajación sometiendo el elemento a una deformación constante "(t) = "̄ y
razonando que la tensión inicial vale �(0) = (E
1
+E∞)"̄ que permite obtener
R(t) = E∞ +E1 e−t�⌧ . (6.14)
6.1.4. El modelo de Maxwell generalizado
El modelo del sólido estándar se puede generalizar, incrementando el
número de elementos de Maxwell en paralelo, como en el ejemplo de la Figu-
ra 6.7. Cuando se aumenta el número de elementos de Maxwell, el elemento
resultante tiene un mayor número de tiempos de relajación caracteŕısticos
⌧
k
= ⌘
k
�E
k
, (6.15)
y el módulo de relajación que resulta es
R(t) = E∞ + K�
k=1
E
k
e−t�⌧k , (6.16)
siendo K el número de elementos de Maxwell en paralelo. Este modelo
reológico recibe el nombre de modelo de Maxwell generalizado o mo-
delo de Wiechert .
Caṕıtulo 6. Viscoelasticidad 127
E1 ⌘1
�1
E1
E2 ⌘2
�2
E3 ⌘3
�3
E4 ⌘4
�4
� �
"
Figura 6.7: Modelo de Maxwell generalizado con 5 elementos de Maxwell en
paralelo
En general, el valor de las rigideces E∞,E
k
y de los tiempos de relaja-
ción ⌧
k
de cada elemento se escoge de forma que la función de relajación
se ajuste a los valores obtenidos de forma experimental. Las series de la
forma (6.16) se llaman series de Prony y existen varios métodos espe-
cialmente diseñados para seleccionar los parámetros de Prony y ajustar la
respuesta ([3, 5, 1]). De hecho, este tipo de series se emplea a menudo en
teoŕıa de la señal.
6.1.5. La integral de Duhamel
Todos los modelos reológicos explicados en esta sección se describen con
ecuaciones diferenciales lineales de coeficientes constantes. Este resultado
expresa anaĺıticamente el principio de Boltzmann que establece que la
respuesta viscoelástica es función, en cada instante, de la historia de defor-
mación y que cada escalón de deformación contribuye de forma indepen-
diente al valor de la tensión, siendo el valor total la suma de cada una de las
contribuciones. Este resultado se puede utilizar para obtener la respuesta
128 Caṕıtulo 6. Viscoelasticidad
0
20
40
60
80
100
0 5 10 15 20
E
(t
)/
E
1
t
"̄ = 5
"̄ = 10
"̄ = 20
Figura 6.8: Módulos de relajación para tres escalones de deformación distin-
tos. Ilustración de la linealidad en la respuesta viscoelástica
en tensión a cualquier deformación "(t), no necesariamente la del ensayo de
relajación.
Para ello, recordamos que la función de Heaviside, o escalón unidad, se
define como
H
⇠
(t) = �������0 si t < ⇠ ,1 si t ≥ ⇠ , (6.17)
verificándose además que
H
⇠
(t) =H
0
(t − ⇠) . (6.18)
Llamemos �
R
[f] a la respuesta reológica correspondiente a una deforma-
ción f . El módulo de relajación R(t) se ha calculado hasta ahora como la
respuesta a un escalón unitario de deformación en el instante t = 0, es de-
cir, R(t) = �
R
[H
0
](t). Sin embargo, si el escalón se produce en t = ⇠, es
immediato comprobar que la respuesta es simplemente �
R
[H
⇠
] = R(t − ⇠).
Dicho de otra manera, la respuesta a una deformación "(t) = H(t − ⇠) es
R(t−⇠). El principio de Boltzmann establece además que la respuesta a una
deformación "(t) =�"̄
1
H
⇠1(t) +�"̄2H⇠2(t) es
�
R
[�"̄
1
H
⇠1(t) +�"̄2H⇠2(t)] =�"̄1 �R[H⇠1](t) +�"̄2 �R[H⇠2](t)=�"̄
1
R(t − ⇠
1
) +�"̄
2
R(t − ⇠
2
) . (6.19)
En el ĺımite, cuando el número de escalones es infinito, la deformación se
puede escribir como
"(t) = � t
0
"̇(⇠)H
⇠
(t)d⇠ , (6.20)
Caṕıtulo 6. Viscoelasticidad 129
y por tanto la tensión como
�(t) = �
R
�� t
0
"̇(⇠)d⇠�
= � t
0
�
R
�"̇(⇠)H
⇠
(t)� d⇠
= � t
0
"̇(⇠)R(t − ⇠)d⇠
(6.21)
Este integral, conocida como la integral de Duhamel , es una convolución
de la función de relajación y la tasa de deformación, y escribimos
� = "̇ ∗R , (6.22)
lo cual sugiere que el análisis de ésta se puede simplificar empleando la
transformada de Laplace.
6.2. Respuesta en frecuencia
Cuando un material elástico se somete a una tensión que vaŕıa en el
tiempo, su deformación también será variable en el tiempo pero la relación
entre tensión y deformación será siempre constante. Cuando un material
viscoelástico es sometido a una tensión variable, la deformación se desfase
con respecto a la tensión.
Supongamos un material viscoelástico sometido a una deformación armóni-
ca
"(t) = "̄ cos(!t) , (6.23)
siendo ! su frecuencia. En general, la tensión no estará en fase con la defor-
mación sino que será de la forma
�(t) = �̄ cos(!t + �) = �̄ cos(�) cos(!t) − �̄ sin(�) sin(!t) , (6.24)
que se puede escribir de forma alternativa como
�(t) = �̄′ cos(!t) − �̄′′ sin(!t) , siendo ��������̄
′ = �̄ cos(�) ,
�̄′′ = �̄ sin(�) . (6.25)
Estudiando la relación entre deformación y tensión usando números com-
plejos, tal y como se hace en teoŕıa de circuitos, observamos que la defor-
mación y tensión se pueden expresar como las partes reales de los complejos
"c = "̄ ei!t , �c = �̄ ei(!t+�) , (6.26)
donde � se puede interpretar ahora como el ángulo de adelanto de la tensión
respecto a la deformación. Reescribiendo la tensión como
�c = �̄ ei� ei!t = (�̄ cos(�) + i �̄ sin(�)) ei!t = (�̄′ + i�̄′′) ei!t , (6.27)
130 Caṕıtulo 6. Viscoelasticidad
e imponiendo una relación entre tensión y deformación de la forma
�c = Ec "c , (6.28)
se sigue que el módulo de rigidez complejo Ec ha de ser de la forma
Ec = E′ + i E′′ con E′ = �̄′
"̄
, E′′ = �̄′′
"̄
, (6.29)
Estas dos cantidades se conocen, respectivamente, como el módulo de al-
macenamiento y de pérdida del material y su ratio coincide con la tan-
gente de �, dependen de la frecuencia !. Las relaciones E′(!) y E′′(!′′)
caracterizan completamente la respuesta viscoelástica lineal de un material,
y por ello se obtienen habitualmente y se emplean para describir las propie-
dades mecánica de los materiales viscoelásticos. En ocasiones, se emplean
expresiones simplificadas para estos dos módulos como por ejemplo
Ec = Ec
o
� !
2⇡
�a , (6.30)
siendo Ec
o
una constante compleja y a un exponente real. Sin embargo, es
más habitual describir la respuesta en frecuencia de un material mediante
una representación gráfica como la de la figura 6.9.
Figura 6.9: Módulos de almacenamiento y pérdida como funciones de la
frecuencia en un policarbonato a 25oC (de Wikipedia).
La razón por la que E′ y E′′ se conocen con los nombres de módulo de
almacenamiento y pérdida tiene relación con el balance energético sobre un
elemento visco-elástico. Para comprender el significado del módulo de pérdi-
da E′′, supongamos un material viscoelástico sometido a una deformación
Caṕıtulo 6. Viscoelasticidad 131
armónica como (6.23). Entonces, el trabajo que hay que hacer en un ciclo
de carga, por unidad de volumen del materiales:
W = � 2⇡�!
0
�(t) "̇(t)dt
= � 2⇡�!
0
(E′"̄ cos!t −E′′"̄ sin!t)(−!)"̄ sin!tdt
= � 2⇡�!
0
"̄2 !(−E′ cos!t sin!t +E′′ sin2 !t)dt= "̄2 ⇡E′′ .
(6.31)
Se deduce, por tanto, que las pérdidas, es decir, el trabajo necesario para
deformar armónicamente el material, es proporcional al módulo de pérdida.
La componente de la tensión que está en fase con la deformación, �̄′ no
produce trabajo neto en un ciclo. Este trabajo es una función armónica de
periodo ⇡�! cuyo valor máximo instantáneo, 1
2
E′"̄2!, que es proporcional
al módulo de almacenamiento.
Como el módulo de rigidez complejo es el cociente entre la tensión y la
deformación complejas se sigue que
Ec = �c
"c
= �̄ei!tei�
"̄ei!t
= �̄ei�
bar"
= Ēei� (6.32)
por lo que Ec es un complejo de módulo Ē, ángulo constante �, igual al
ángulo de desfase entre la tensión y la deformación complejas. El inverso del
módulo complejo es la flexibilidad compleja
Jc = (Ec)−1 = E′ − iE′′(E′)2 + (E′′)2 = J ′ − iJ ′′ (6.33)
que permite expresar la relación constitutiva "c = Jc�c.
6.2.1. Caracterización en frecuencia de los modelos reológi-
cos elementales
Cuando la deformación es de la forma (6.23), la tensión en un elemento
de Kelvin es
�(t) = E cos(!t) − ⌘! sin(!t) = E cos(!t) −E⌧! sin(!t) . (6.34)
Identificando en esta expresión obtenemos que los módulos de almacena-
miento y pérdida en un elemento de Kelvin son:
E′(!) = E , E′′(!) = E⌧! (6.35)
De la misma forma se puede calcular anaĺıticamente la expresión de los
módulos de almacenamiento y pérdida en el modelo de Maxwell:
E′(!) = E(!⌧)2
1 + (!⌧)2 , E′′(!) = E(!⌧)1 + (!⌧)2 , (6.36)
132 Caṕıtulo 6. Viscoelasticidad
1e-05
0.0001
0.001
0.01
0.1
1
10
100
0.01 0.1 1 10 100
! ⌧
E0/E
E00/E
1e-05
0.0001
0.001
0.01
0.1
1
10
100
0.01 0.1 1 10 100
! ⌧
E0/E
E00/E
Figura 6.10: Módulos de almacenamiento y pérdida para los modelos de
Kelvin (arriba) y Maxwell (abajo).
Caṕıtulo 6. Viscoelasticidad 133
En la Figura 6.10 se puede apreciar la dependencia de los módulos de
almacenamiento y pérdida en los modelos de Kelvin Maxwell. En ambos
casos se identifica un punto especial que corresponde a ! = 1�⌧ .
6.2.2. Series de Prony
Como se mencionó anteriormente, dada la respuesta en frecuencia de un
material como por ejemplo el de la Figura 6.9, se puede definir un modelo de
Maxwell generalizado cuya respuesta en frecuencia se ajuste con la tolerancia
deseada a los datos experimentales del material. El método que se utiliza
habitualmente es el basado en las llamadas series de Prony .
Las series de Prony describen las funciones de relajación de la forma
R
Prony
(t) = E
o
�1 − n�
i=1pi(1 − e−t�⌧i)� , (6.37)
siendo E
0
el módulo de rigidez instantánea, p
i
los coeficientes de Prony y ⌧
i
los tiempos de relajación caracteŕısticos.
Dado un modelo basado en series de Prony, se puede demostrar que los
módulos de almacenamiento y pérdida tienen por expresión [6]:
E′(!) = E
o
�1 − n�
i=1pi� +Eo
n�
i=1
p
i
⌧2
i
!2
1 + ⌧2
i
!2
,
E′′(!) = E
o
n�
i=1
p
i
⌧
i
!
1 + ⌧2
i
!2
.
(6.38)
No existe ninguna fórmula que invierta esta relación, es decir, que permita
calcular directamente los parámetros de Prony a partir de los módulos de
almacenamiento y pérdida.
6.3. Sólidos deformables viscoelásticos
En este caṕıtulo se ha estudiado, hasta ahora, la respuesta viscoelástica
de modelos reológicos unidimensionales. El objetivo ha sido presentar los
aspectos fundamentales del comportamiento viscoelástico de la manera más
sencilla posible. Utilizando los resultados obtenidos se puede calcular la res-
puesta de un elemento unidimensional viscoelástico (una barra, por ejemplo)
cuando se somete a una historia de tensión o deformación cualquiera. En es-
ta sección se explica cómo todos ellos se emplean para describir el modelo
constitutivo viscoelástico de un punto en un sólido tridimensional.
En primer lugar se observa experimentalmente que, al igual que en el
caso elástico, el comportamiento volumétrico y desviador en un sólido vis-
coelástico están desacoplados. Es decir, si "(t) = 1
3
✓(t)1 + e(t), siendo ✓(t)
134 Caṕıtulo 6. Viscoelasticidad
la deformación volumétrica y e(t) la parte desviadora de la deformación,
entonces
p(t) = p(✓(t)) , s(t) = s(e(t)) . (6.39)
Además se observa también experimentalmente que de forma bastante pre-
cisa se puede suponer que la respuesta volumétrica es totalmente elástica,
es decir, que en los sólidos viscoelásticos
p(t) = ✓(t) , (6.40)
siendo  una constante que, como en el caso elástico, se llama el módulo de
rigidez volumétrico.
La relación entre las partes desviadoras de la deformación y tensión no
es tan sencilla, sino que claramente existen efectos reológicos que hay que
considerar. El modelo viscoelástico más habitual se construye extendiendo
las ideas de los modelos reológicos. Cuando un punto se somete a una defor-
mación desviadora constante e(t) = ē, el estado tensional experimenta una
relajación de tensiones que se puede expresar como
s(t) = G(t)ē , (6.41)
siendo G(t) el módulo de relajación a cortante del material. Este módu-
lo se suele expresar a partir del módulo de relajación E(t) extrapolando la
relación elástica
G(t) = E(t)
2(1 + ⌫) , (6.42)
siendo ⌫ el coeficiente de Poisson del material, también constante. Si E(t)
se expresa en forma de una serie de Prony, también el módulo de relajación
a cortante se podrá expresar como
G(t) = G
o
�1 − K�
k=1
p
k
(1 − e−t�⌧k)� , (6.43)
donde ⌧
k
son también los tiempos de relajación caracteŕısticos. Por último,
dada una deformación "(t) podemos extrapolar la integral de Duhamel para
escribir
�(t) =  tr("(t))
3
1 +� t
0
ė(⇠)G(t − ⇠)d⇠ (6.44)
6.4. Efectos de la temperatura en la respuesta vis-
coelástica
La temperatura a la que un material viscoelástico se encuentra modifica
sustancialmente su respuesta. Afortunadamente, en muchos casos, la depen-
dencia en la temperatura se puede aproximar de forma sencilla y con una
precisión suficiente.
Caṕıtulo 6. Viscoelasticidad 135
Figura 6.11: Curva maestra de relajación para poli-iso-butileno y dependen-
cia de la misma con la termperatura.
En general, el aumento de la temperatura disminuye la viscosidad de los
materiales, o lo que es lo mismo, los tiempos de relajación. Para muchos
materiales existe una correspondencia temperatura-tiempo que se ma-
nifiesta en que un aumento de la temperatura traslada la curva de relajación
hacia la derecha, siendo este desplazamiento función monotónica del salto
térmico. Llamando a
T
a este desplazamiento (tiempo) en la temperatura T ,
la ecuación de Williams-Landel-Ferry proporciona
log a
T
= −17,44(T − Tg)
51,6 + T − T
g
, (6.45)
siendo T
g
la temperatura de referencia a la cual se proporciona la curva
maestra de relajación.
Problemas
6.1. (Termodinámica del modelo estándar) El modelo estándar viene defini-
do por una enerǵıa libre (para problemas isotermos) y una relación cinética
A(",�) = 1
2
E∞"2 + 1
2
E
1
(" − �)2 , �̇ = ⌘E
1
(" − �).
a) Interpreta las dos contribuciones de la enerǵıa libre.
b) Razona por qué la ecuación cinética expresa la relación constitutiva
habitual del amortiguador �̇ = ⌘�.
136 Caṕıtulo 6. Viscoelasticidad
0.1
1
10
100
0.01 0.1 1 10 100 1000
l
o
g
(
E
(
t)
/E
1
)
log(t/⌧)
Figura 6.12: Ejemplo de la correspondencia temperatura-tiempo. La curva
verde se “desplaza” hacia la derecha cuando la temperatura aumenta. El
desplazamiento es constante y únicamente función de la temperatura.
c) Utilizando los resultados del caṕıtulo 5, demuestra que el modelo an-
terior es termodinámicamente correcto.
� �
"
�
⌘
E1
E2
Figura 6.13: Sólido viscoelástico estándar basado en el modelo de Kelvin.
6.2. Para el modelo estándar de la figura 6.13,
a) Encuentra la ecuación diferencial que gobierna su respuesta.
b) Encuentra la función de fluencia.
c) Calcula la función de relajación.
6.3. Calcula losmódulos de pérdida y almacenamiento de un elemento
reológico que tiene por función de relajación R(t) = R∞ +R1e−t�⌧ ,
6.4. Un peso de 10 kg se cuelga con una cuerda de 4 m de longitud y 4 mm2
de sección. El material de la cuerda es viscoelástico, y su comportamiento se
Bibliograf́ıa 137
puede describir con un modelo estándar de constantes E
1
= 3 GPa, E∞ = 2
GPa, ⌘ = 1 GPa⋅s. Dibuja un gráfica con la evolución de la longitud de la
cuerda e indica cuál es la longitud máxima que ésta alcanzará.
0
10
20
30
40
50
0 2 4 6 8 10 12 14
�
(M
P
a)
t (s)
Figura 6.14: Historia de cargas del problema 6.5.
6.5. Un material viscoelástico se representa con un modelo reológico de
Kelvin de constantes E = 2 GPa y ⌘ = 1 GPa⋅s. Calcular la historia de
deformación en el mismo cuando se somete a una tensión como la indicada
en la figura 6.14.
6.6. Una máquina rotatoria se coloca sobre una base de volumen 0.01 m3
de forma que se puede considerar, en una primera aproximación, que el
estado de carga de la base es de tracción uniaxial. El material de la base es
viscoelástico y se conoce su respuesta en frecuencia en algunos puntos:
! E′ E′′
(Hz) (GPa) (GPa)
20 2.1 1.0
40 2.0 1.3
60 1.8 1.5
80 1.1 2.1
Si la máquina está funcionando a 50 Hz y la tensión que se ejerce sobre
la base es armónica con amplitud 20 MPa, calcular el calor disipado por la
base en 10 minutos de funcionamiento.
Bibliograf́ıa
[1] Tzikang Chen. Determining a prony series for a viscoelastic material from
time varying strain data. Technical Report NASA/TM-2000-210123,
U.S. Army Research Laboratory, 2000.
[2] Wilhem Flügge. Viscoelasticity. Blaisdell Publishing Company, 1967.
138 Bibliograf́ıa
[3] SW Park and RA Schapery. Methods of interconversion between linear
viscoelastic material functions. Part I—a numerical method based on
Prony series. International Journal of Solids and Structures, 36(11):
1653–1675, 1999.
[4] Nhan Phan-Thien. Understanding viscoelasticity. Basis of rheology.
Springer, 2002.
[5] R A Schapery and S W Park. Methods of interconversion between linear
viscoelastic material functions. Part II—An approximate analytical met-
hod. International Journal of Solids and Structures, 36(11):1677–1699,
1999.
[6] Si. ABAQUS theory manual, 6.7 edition.