Viscoelasticidad

Esquema de diferentes comportamientos de los materiales.

La viscoelasticidad es un tipo de comportamiento reológico que presentan ciertos materiales que exhiben tanto propiedades viscosas como propiedades elásticas cuando se deforman.

En un sólido viscoelástico:

• la deformación generalmente depende del tiempo; aún en ausencia de fuerzas, la velocidad de deformación puede ser diferente de cero;
• las tensiones y esfuerzos resistidos dependen tanto de la deformación como de la velocidad de deformación, por tanto la ecuación constitutiva que relaciona tensiones y deformaciones debe tener la forma : $\phi (\sigma _{ij},{\dot {\sigma }}_{ij},\dots ,\varepsilon _{ij},{\dot {\varepsilon }}_{ij},\dots ,t)=0$.

Físicamente las propiedades elásticas son el resultado de desplazar ligeramente los átomos de su posición de equilibrio a lo largo de planos cristalográficos, mientras las propiedades viscosas proceden de la difusión de átomos o moléculas en el interior del material​.

Viscoelasticidad lineal

Un material viscoelástico lineal general es un material para el cual existe una relación lineal entre la tensión y sus derivadas y la deformación y sus derivadas, en el caso unidimensional la relación más general posible de un material viscoelástico lineal es:​

(1a)$\sigma +A^{(1)}{\frac {\partial \sigma }{\partial t}}+\dots +A^{(n_{1})}{\frac {\partial ^{n_{1}}\sigma }{\partial t^{n_{1}}}}=B^{(0)}\varepsilon +B^{(1)}{\frac {\partial \varepsilon }{\partial t}}+\dots +B^{(n_{2})}{\frac {\partial ^{n_{2}}\varepsilon }{\partial t^{n_{2}}}}$

En este caso usando transformadas de Laplace y si $\scriptstyle n_{2}\geq n_{1}$ y $\scriptstyle \sigma (0)=\varepsilon (0)=0$, la expresión (1a) puede escribrise simplmente como:

(1b)$\sigma (t)=a^{(0)}\varepsilon (t)+a^{(1)}{\frac {\partial \varepsilon }{\partial t}}+\dots +a^{(N)}{\frac {\partial ^{N}\varepsilon }{\partial t^{N}}}+\int _{0}^{t}b(t-\tau )\varepsilon (\tau )\ d\tau$

Siendo $\scriptstyle N=n_{2}-n_{1}$. Cuando $\scriptstyle n_{2}<n_{1}$, la expresión equivalente a (1a) es más simple:

(1c)$\sigma (t)=\int _{0}^{t}b(t-\tau )\varepsilon (\tau )\ d\tau$

Aunque esta expresión no se usa porque en la práctica no se conocen materiales que se ajusten bien a este tipo de ley. Así que aquí se restringirá la discusión a $\scriptstyle n_{2}\geq n_{1}$.

Los dos modelos viscoelásticos más conocidos, el de Maxwell y el de Kelvin-Voigt son casos particulares de (1a) que satisfacen $\scriptstyle n_{1}\leq n_{2}\leq 1$. En ambos $\scriptstyle n_{2}=1$, y $\scriptstyle n_{1}=0$ para el de Kelvin-Voigt y $\scriptstyle n_{1}=1$ para el de Maxwell.

Creep y relajación

Bajo los mismos supuestos anteriores, pude probarse que un material viscoelástico lineal admite una caracterización matemática en que la respuesta de creep y de relajación de carga admiten una separación en forma de suma. Los modelos de viscoelasticidad lineal pueden ser representados por la ecuación de Volterra conectando la tensión y la deformación mediante una expresión del tipo:

(2a)$\varepsilon (t)={\frac {\sigma (t)}{E_{0,C}}}+\int _{0}^{t}K(t-\tau ){\dot {\sigma }}(\tau )\ d\tau$

o también :

(2b)$\sigma (t)=E_{0,R}\varepsilon (t)+\int _{0}^{t}F(t-\tau ){\dot {\varepsilon }}(\tau )\ d\tau$

donde :

$t\,$, es el tiempo.

$\sigma (t)\,$, es la tensión mecánica.

$\varepsilon (t)\,$, es la deformación o elongación.

$E_{0,C}\,$ y $E_{0,R}\,$, son los módulos elásticos longitudinales para el creep y la relajación.

$K(t)\,$, es la función de creep.

$F(t)\,$, es la función de relajación.

La viscoelasticidad lineal tiene un rango de aplicación válido sólo para deformaciones muy pequeñas.

Ecuaciones constitutivas de viscoelasticidad lineal

Existen diversos modelos constitutivos para materiales viscoelásticos lineales. Esos modelos incluyen el modelo de Maxwell, el modelo de Kelvin-Voigt y el modelo de sólido viscoelástico lineal estándar que combina los dos modelos anteriores. Todos estos modelos descomponen la tensión y deformación en dos sumandos, uno que representa los efectos elásticos y otro que representa los efectos viscosos, siendo estos modelos, interpretables en términos de muelles y amortiguadores. Cada uno de estos modelos difiere en la disposición de los muelles y amortiguadores.

Otra propiedad interesante es que las ecuaciones constitutivas también pueden ser interpretadas en términos de circuitos eléctricos, en los que la tensión mecánica sería el equivalente del voltaje y la velocidad de deformación sería equivalente a la intensidad de corriente. El módulo elástico sería equivalente a la capacitancia del circuito (que mide la capacidad de almacenaje de energía) y la viscosidad a la resistencia del circuito (que mide la capacidad de disipar energía).

Modelo de Maxwell

Representación esquemática del modelo de Maxwell.

El modelo de Maxwell es un caso particular de la expresión (1a) en el que $\scriptstyle n_{1}=n_{2}=1$, también llamado material viscoelástico de "larga memoria". Una virtud del modelo de Maxwell es que admite una representación intuitiva en términos de muelles y disipadores (amortiguadores), más concretamente el modelo de Maxwell representa un elemento elástico dispuesto en serie con un amortiguador. La ecuación constitutiva del modelo viene dada por la siguiente ecuación diferencial de primer orden:

(3a)${\frac {d\varepsilon }{dt}}={\frac {d\varepsilon _{A}}{dt}}+{\frac {d\varepsilon _{R}}{dt}}={\frac {\sigma }{\eta }}+{\frac {1}{E}}{\frac {d\sigma }{dt}}$.

Este modelo predice que en un material puesto bajo deformación constante, las tensiones gradualmente se relajarán hasta hacerse cero. Igualmente predice que si el material se pone a tensión constante la deformación tendrá dos componentes, un primer componente elástico que aparece instantáneamente y una deformación a largo plazo de tipo viscoso que crecerá con el tiempo mientras sigan existiendo fuerzas. Alternativamente la ecuación constitutiva de este modelo puede escribirse como:

$\sigma (t)=E\varepsilon (t)+\int _{0}^{t}b(t-\tau )\varepsilon (\tau )\ d\tau =E\varepsilon (t)-{\frac {E^{2}}{\eta }}\int _{0}^{t}e^{-{\frac {E}{\eta }}(t-\tau )}\varepsilon (\tau )\ d\tau$.

Integrando por partes el último término se puede obtener reescribir la expresión anterior en términos de la función de relajación $\scriptstyle F(\cdot )$ y la velocidad de deformación:

(3b)$\sigma (t)=E\varepsilon (t)+\int _{0}^{t}\left[Ee^{-{\frac {E}{\eta }}(t-\tau )}{\dot {\varepsilon }}(\tau )\right]d\tau ,\qquad F(\tau )=Ee^{-{\frac {E}{\eta }}\tau }$.

Para obtener la de función de creep $\scriptstyle K(\cdot )$ se integra directamente (3a), e integrando por partes se obtiene:

(3c)$\varepsilon (t)={\frac {\sigma }{E}}+\int _{0}^{t}{\frac {1}{\eta }}\sigma (\tau )d\tau ={\frac {\sigma }{E}}+\int _{0}^{t}\left[{\frac {t-\tau }{\eta }}\right]{\dot {\sigma }}(\tau )d\tau ,\qquad K(\tau )={\frac {\tau }{\eta }}$.

Y no es posible encontrar una función ordinaria para expresar lo anteior en la forma (2a).

El modelo de Maxwell predice que la tensión decaerá exponencialmente con el tiempo en un polimero sometido a deformación constante, lo cual se ajusta bastante bien a lo observado experimentalmente para muchos polímeros. Sin embargo, una limitación importante es que no predice el comportamiento de "creep" de muchos polímeros de manera demasiado fidedigna ya que en este caso predice un aumento lineal de la deformación con el tiempo si la tensión es constante, sin embargo la mayor parte de los polímeros muestran una tasa de deformación decreciente con el tiempo. Las principales aplicaciones de este modelo son la modelización de los polímeros termoplásticos cerca de su temperatura de fusión, la del hormigón fresco y la de numerosos metales cerca de su punto de fusión.

Modelo de Kelvin-Voigt

Representación esquemática del modelo de Kelvin-Voigt.

El modelo de Kelvin-Voigt o modelo de Voigt es un caso particular de la expresión (1a) en el que $\scriptstyle n_{1}=0,\ n_{2}=1$, también llamado material viscoelástico de "corta memoria". Al igual que el modelo anterior admite una reprsensentación simple en términos de muelles y amortiguadores: el modelo es representable por un amortiguador newtoniano y un muello que sigue la ley de Hooke conectado en paralelo al amortiguador, tal como muestra la figura. La ecuación constitutiva del modelo puede expresarse como ecuación diferencial de primer orden:

(4a)$\sigma (t)=E\varepsilon (t)+\eta {\frac {d\varepsilon (t)}{dt}}$

Este modelo representa un sólido que sufre deformación viscoplástica reversible. Bajo la aplicación de una tensión constante el material se deforma a un ritmo progresivamente más lento, llegando asintóticamente a un estado cuasiestacionario. Cuando se eliminan las fuerzas exteriores que generan las tensiones, el material se relaja hasta su estado no deformado original. En una situación de tensión constante (creep), el modelo es bastante realista y prdice deformaciones que tienden al límite σ/E para tiempos grandes.

En este caso la función de relajación $\scriptstyle F(\cdot )$ viene dada por la delta de Dirac (de ahí el nombre de "corta memoria") ya que:

(4b)$\sigma (t)=E\varepsilon (t)+\int _{0}^{t}F(t-\tau ){\dot {\varepsilon }}(\tau )\ d\tau =E\varepsilon (t)+\eta {\dot {\varepsilon }}(t),\quad \Rightarrow \quad F(t-\tau )=\eta \delta (t-\tau )$

Para obtener la función de creep $\scriptstyle K(\cdot )$ se busca la transformada de Laplace de (4a) de donde es fácil aislar la transformada de la deformación y aplicando la transformada inversa se llega fácilmente a que:

$\varepsilon (t)=\int _{0}^{t}{\frac {e^{-{\frac {E}{\eta }}(t-\tau )}}{\eta }}\sigma (\tau )\ d\tau$

Integrando por partes esta úlitma expresión se llega a una expresión que contiene la función de creep o fluencia lenta:

(4c)$\varepsilon (t)={\frac {\sigma }{E}}+\int _{0}^{t}{\frac {e^{-{\frac {E}{\eta }}(t-\tau )}}{E}}{\dot {\sigma }}(\tau )\ d\tau ,\qquad K(\tau )={\frac {e^{-{\frac {E}{\eta }}\tau }}{E}}$

Este modelo se usa para explicar el comportamiento de "creep" de los polímeros. Aunque, al igual que el modelo de Maxwell, el modelo de Kelvin-Voigt tiene limitaciones empíricas. Aunque modela muy bien el "creep" con respecto a la relajación el modelo generalmente se ajusta menos al comportamiento de los materiales viscoelásticos. Las aplicaciones principales del modelo son la modelización de polímeros orgánicos, goma, caucho y madera cuando la carga no es muy elevada.

Modelo estándar de sólido viscoelástico

Representación esquemática del modelo estándar de sólido viscoelástico.

Los dos úlitmos modelos expuestos tienen limitaciones en cierto modo complementarias. Por esa razón se considera que un modelo que combine características de ambos puede ser un modelo razonable de sólido viscoelástico. La ecuación constitutiva de este modelo viene dada por la siguiene ecueción diferencial:

(5a)$(E_{1}+E_{2}){\frac {d\varepsilon }{dt}}={{\frac {E_{2}}{\eta }}\left({\frac {\eta }{E_{2}}}{\frac {d\sigma }{dt}}+\sigma -E_{1}\varepsilon \right)}$

Usando transformadas de Laplace, despejando la transformada de la tensión y antitransformando se obtiene una ecuación de tipo Volterra como en los casos anteriores, que tiene la forma:

$\sigma (t)=(E_{1}+E_{2})\varepsilon (t)-{\frac {E_{2}^{2}}{\eta }}\int _{0}^{t}e^{-{\frac {E_{2}}{\eta }}(t-\tau )}\varepsilon (\tau )\ d\tau$

Una integración por partes del último término permite escribir en términos de la función de relajación $\scriptstyle F(\cdot )$ y la velocidad de deformación:

(5b)$\sigma (t)=(E_{1}+E_{2})\varepsilon (t)\varepsilon (t)+\int _{0}^{t}\left[E_{2}e^{-{\frac {E_{2}}{\eta }}(t-\tau )}{\dot {\varepsilon }}(\tau )\right]d\tau ,\qquad F(\tau )=E_{2}e^{-{\frac {E_{2}}{\eta }}\tau }$

Modelo de Burgers

(6a)$\sigma +\left({\frac {1}{E_{1}}}+{\frac {1}{E_{2}}}\right){\dot {\sigma }}+{\frac {1}{E_{1}E_{2}}}{\ddot {\sigma }}=(\eta _{1}+\eta _{2}){\dot {\varepsilon }}+\left({\frac {\eta _{1}}{E_{2}}}+{\frac {\eta _{2}}{E_{1}}}\right){\ddot {\varepsilon }}$

Viscoelasticidad no lineal

La viscosidad no lineal ocurre cuando las funciones de creep y de relajación de carga no se pueden separar. Esto ocurre en deformaciones grandes, si el material cambia sus propiedades durante la deformación, si los tiempos de la deformación son lo suficientemente largos o si interviene algún otro tipo de relajación.

Referencias

Bibliografía

• Freudenthal, A. M. y Geiringer, H: "The mathematical theories of inelastic continuum", Handbuch der Physik - Encyclopedia of Physics, VII, Elasticity and Plasticity. Berlin-Göttingen-Heidelberg: Springer 1958.

Véase también

• Viscoelasticidad de polímeros