La viscoplasticidad es un tipo de comportamiento inelástico de los sólidos deformables en el que aparecen deformaciones permanentes y su aparición depende además de la velocidad de deformación. El comportamiento viscoplástico es la base de los modelos constitutivos viscoplásticos, que es una rama de la teoría de la mecánica del continuo. La dependencia de la velocidad de deformación en este contexto significa que la deformación del material depende de la velocidad a la que se aplican las carga.[1]
Al igual que en el comportamiento plástico, en el comportamiento viscoelástico el material sufre deformaciones irrecuperables cuando se alcanza un nivel de carga. La diferencia es que en plasticidad el comportamiento es independiente de la velocidad de deformación. De hecho las teorías de plasticidad pueden considerarse una primera aproximación al verdadero comportamiento viscoelástico, cuyas ecuaciones más complejas son necesarias para describir la plasticidad transitoria, especialmente cuando la velocidad de deformación es alta. Otra diferencia entre el comportamiento plástico y viscoplástico es que en este último caso, no solo aparecen deformaciones permanentes después de la aplicación de cargas, sino que el material continúa experimentando un flujo de fluencia en función del tiempo bajo la influencia de la carga aplicada. Por otra parte, la viscoplasticidad y la viscoelasticidad comparten el hecho de que las tensiones dependen de la velocidad de deformación, aunque en el caso puramente viscoelástico eso no se traduce en la aparición de deformaciones permanentes o irreversibles.
La viscoplasticidad tridimensional se modeliza mediante modelos de sobreesfuerzo de los tipos Perzyna o Duvaut-Lions. [2] En estos modelos, se permite que la tensión aumente más allá de la superficie de fluencia independiente de la tasa al aplicar una carga y luego se permite que se relaje de nuevo a la superficie de fluencia con el tiempo. Por lo general, se supone que la superficie de fluencia no depende de la velocidad de deformación en estos modelos. Un enfoque alternativo es añadir una dependencia de la velocidad de deformación al límite elástico y utilizar las técnicas de plasticidad para calcular la respuesta de un material[3]
Para los metales y las aleaciones, la viscoplasticidad es el comportamiento macroscópico causado por un mecanismo asociado al movimiento de dislocaciones en los granos, con efectos superpuestos de deslizamiento intercristalino. El mecanismo generalmente se vuelve dominante a temperaturas superiores a aproximadamente un tercio de la temperatura absoluta de fusión. Sin embargo, ciertas aleaciones exhiben viscoplasticidad a temperatura ambiente (300K). Para los polímeros, la madera y el betún, se requiere la teoría de la viscoplasticidad para describir el comportamiento más allá del límite de elasticidad o viscoelasticidad.
En general, las teorías viscoplásticas son útiles en áreas como:
el cálculo de las deformaciones permanentes,
la predicción del colapso plástico de las estructuras,
la investigación de la estabilidad,
simulaciones de choques,
sistemas expuestos a altas temperaturas, como las turbinas de los motores, por ejemplo, una central eléctrica,
Problemas dinámicos y sistemas expuestos a altas tasas de deformación.
Historia
La investigación sobre las teorías de la plasticidad comenzó en 1864 con el trabajo de Henri Tresca,[4] A. J. C. B. de Saint Venant (1870) y M. Levy (1871)[5] en el criterio de corte máximo.[6] Un modelo de plasticidad mejorado fue presentado en 1913 por Richard von Mises[7] que ahora se conoce como el criterio de Von Mises. En viscoplasticidad, el desarrollo de un modelo matemático se remonta a 1910 con la representación de creep primario por la ley de Andrade.[8] En 1929, Norton[9] desarrolló un modelo unidimensional con amortiguación que relacionaba la velocidad de fluencia secundaria con la tensión mecánica. En 1934, Odqvist[10] generalizó la ley de Norton al caso tridimensional.
Conceptos como la normalidad del flujo plástico a la superficie de fluencia y las reglas de flujo para la plasticidad fueron introducidos por L. Prandtl (1924)[11] y Reuss (1930). [12] En 1932, Hohenemser y W. Prager[13] propusieron el primer modelo para el flujo viscoplástico lento. Este modelo proporcionó una relación entre el tensor desviador para la tensión y la velocidad de deformación para un sólido de Bingham incompresible[14] Sin embargo, la aplicación de estas teorías no comenzó antes de 1950, cuando se descubrieron los teoremas de límite.
En 1960, el primer Simposio IUTAM "Creep in Structures" organizado por Hoff[15] contribuyó significativamente al desarrollo de teoría de la viscoplasticidad con los trabajos de Hoff, Rabotnov, Perzyna, Hult y Lemaitre para las leyes de endurecimiento isotrópico, y los de Kratochvil, Malinini y Khadjinsky, Ponter y Leckie, y Chaboche para las leyes de endurecimiento cinemático. Perzyna, en 1963, introdujo un coeficiente de viscosidad que depende de la temperatura y el tiempo.[16] Los modelos formulados fueron fundamentados en la termodinámica de los procesos irreversibles y el punto de vista fenomenológico. Las ideas presentadas en estos trabajos han sido la base de la mayoría de las investigaciones posteriores sobre la plasticidad dependiente de la velocidad.
Casuística y observaciones experimentales
Para un análisis cualitativo, se realizan varias pruebas características para describir la fenomenología de los materiales viscoplásticos. Algunos ejemplos de estas pruebas son [8]
ensayos de endurecimiento a una velocidad de tensión o deformación constante,
ensayos de creep a fuerza constante, y
relajación tensiones a elongación constante.
Prueba de endurecimiento por deformación
Respuesta tensión-deformación de un material viscoplástico a diferentes velocidades de deformación. Las líneas punteadas muestran la respuesta si la velocidad de deformación se mantiene constante. La línea azul muestra la respuesta cuando la velocidad de deformación cambia repentinamente.
Un hecho común en la deformación plástica es que en muchos materiales con endurecimiento se requiere un aumento en la tensión para producir una deformación adicional. Este fenómeno se conoce como endurecimiento por deformación/trabajo. [17] Para un material viscoplástico, las curvas de endurecimiento no son significativamente diferentes de las del material plástico independiente de la velocidad. Sin embargo, se pueden observar tres diferencias esenciales.
A la misma deformación, cuanto mayor sea la tasa de deformación, mayor será la tensión
Un cambio en la velocidad de deformación durante la prueba da como resultado un cambio inmediato en la curva tensión-deformación.
El concepto de límite de rendimiento plástico ya no es estrictamente aplicable.
La hipótesis de la partición de las deformaciones mediante el desacoplamiento de las partes elásticas y plásticas sigue siendo aplicable cuando las deformaciones son pequeñas,[2] es decir,
donde es la deformación elástica y es la deformación viscoplástica. Para obtener el comportamiento tensión-deformación que se muestra en azul en la figura adjunta, el material se carga inicialmente a una velocidad de deformación de 0,1/s. A continuación, la velocidad de deformación se eleva instantáneamente a 100/s y se mantiene constante en ese valor durante algún tiempo. Al final de ese período de tiempo, la velocidad de deformación se reduce instantáneamente a 0,1/s y el ciclo continúa para valores crecientes de deformación. Es evidente que existe un desfase entre el cambio de la velocidad de deformación y la respuesta la tensión. Este retraso se modela con bastante precisión mediante modelos de sobresfuerzo (como la Modelo de Perzyna), pero no por modelos de plasticidad independiente de la tasa que tienen un límite elástico dependiente de la tasa.
Ensayo de fluencia
Ensayo de fluencia
Deformación en función del tiempo en un ensayo de fluencia.
La fluencia es la tendencia de un material sólido a deformarse permanentemente bajo tensiones constantes. Los ensayos de fluencia miden la respuesta a la deformación bajo una tensión constante, como se muestra en la figura adjunta. La curva de fluencia clásica representa la evolución de la deformación en función del tiempo en un material sometido a esfuerzos uniaxiales a una temperatura constante. El ensayo de fluencia, por ejemplo, se realiza aplicando una fuerza/tensión constante y analizando la respuesta a la deformación del sistema. En general, como se muestra en la figura adjunta, esta curva suele mostrar tres fases o períodos de comportamiento[8]
Una etapa de fluencia primaria, también conocida como fluencia transitoria, es la etapa inicial durante la cual el endurecimiento del material conduce a una disminución en la tasa de flujo que inicialmente es muy alta. .
La etapa de fluencia secundaria, también conocida como estado estacionario, es donde la velocidad de deformación es constante. .
Una fase de fluencia terciaria en la que hay un aumento en la velocidad de deformación hasta la deformación de fractura. .
Ensayo de relajación
a) deformación aplicada en un ensayo de relajación y b) tensión inducida en función del tiempo durante un corto período de tiempo para un material viscoplástico.
Como se muestra en la figura adjunta, el ensayo de relajación[18] describe como es la respuesta tensional a deformación constante a lo largo de un período de tiempo. En materiales viscoplásticos, las pruebas de relajación muestran una relajación o disminución de la tensión asociada a la carga uniaxial bajo deformación constante. De hecho, estos ensayos caracterizan la viscosidad y se pueden utilizar para determinar la relación que existe entre la tensión y la velocidad de deformación viscoplástica. La descomposición de la velocidad de deformación es
La parte elástica de la velocidad de deformación viene dada por
Para la región plana de la curva de tiempo de deformación, la velocidad de deformación total es cero. Por lo tanto, tenemos,
Por lo tanto, la curva de relajación se puede utilizar para determinar la velocidad de deformación viscoplástica y, por lo tanto, la viscosidad en un modelo de material viscoplástico unidimensional. El valor residual que se alcanza cuando la tensión se estabiliza al final de un ensayo de relajación correspondiente al límite superior elástico. Para algunos materiales, como la sal gema, este límite superior elástico se alcanza a un valor muy pequeño de tensión y los ensayos de relajación pueden continuarse durante más de un año sin que se observe ninguna meseta en la tensión.
Es importante tener en cuenta que los ensayos de relajación son extremadamente difíciles de realizar porque mantener la condición: en una prueba requiere una delicadeza considerable.[19]
Modelos reológicos de viscoplasticidad
Los modelos constitutivos unidimensionales para la viscoplasticidad basados en elementos muelle-amortiguador-deslizador incluyen[2]: (a) el sólido perfectamente viscoplástico, (b) el sólido elástico perfectamente viscoplástico y (c) el sólido con endurecimiento elasto-viscoplástico. Los elementos pueden estar conectados en serie o en paralelo. En los modelos en los que los elementos están conectados en serie, la deformación es aditiva, mientras que la tensión es igual en cada elemento. En las conexiones paralelas, la tensión es aditiva mientras que la deformación es igual en cada elemento. Muchos de estos modelos unidimensionales se pueden generalizar a tres dimensiones para el régimen de pequeñas deformaciones. En la discusión posterior, la tensión y la velocidad de deformación se escriben como y , respectivamente.
Modelo de Norton-Hoff para sólidos perfectamente viscoplásticos.
En un sólido perfectamente viscoplástico, también llamado modelo de viscoplasticidad de Norton-Hoff, la tensión (como en el caso de los fluidos viscosos) es una función de la velocidad de deformación inelástica. El efecto de la elasticidad se desprecia en el modelo, es decir, y, por lo tanto, no hay tensión de fluencia inicial, es decir, . El amortiguador viscoso tiene una respuesta dada por:
donde es la viscosidad del amortiguador. En el modelo de Norton-Hoff, la viscosidad es una función no lineal de la tensión aplicada y viene dada por
donde es un parámetro de ajuste, λ es la viscosidad cinemática del material y . Entonces la velocidad de deformación viscoplástica viene dada por la relación
En su forma unidimensional, el modelo de Norton-Hoff se puede expresar como
Si asumimos que el flujo plástico es isocórico (preserva el volumen), entonces la relación anterior se puede expresar en la forma más familiar[20]
donde es el tensor desviador, es velocidad de deformación equivalente de von Mises, y son parámetros materiales. La velocidad de deformación equivalente se define como:
Estos modelos se pueden aplicar a metales y aleaciones a temperaturas superiores a dos tercios[20] de su punto de fusión absoluto (en kelvins) y a polímeros y asfalto a temperatura elevada. Las respuestas para las pruebas de endurecimiento por deformación, fluencia y relajación de dicho material se muestran en la figura adjunta.
Figura 6: Respuesta de un sólido perfectamente viscoplástico a los ensayos de endurecimiento, fluencia y relajación.
Se pueden utilizar dos tipos de enfoques elementales para construir un modo elástico-perfectamente viscoplástico. En la primera situación, el elemento de fricción deslizante y el amortiguador se disponen en paralelo y luego se conectan en serie al resorte elástico como se muestra en la Figura 7. Este modelo se llama modelo de Bingham-Maxwell (por analogía con el modelo de Maxwell y el de Bingham) o el Modelo de Bingham-Norton.[21] En la segunda situación, los tres elementos están dispuestos en paralelo. Un modelo de este tipo se denomina modelo Bingham-Kelvin por analogía con el modelo Kelvin.
En el caso de los materiales elásticos perfectamente viscoplásticos, la deformación elástica ya no se considera despreciable, sino que la tasa de deformación plástica es sólo una función del límite elástico inicial y no influye el endurecimiento. El elemento deslizante representa una tensión de fluencia constante cuando se supera el límite elástico, independientemente de la deformación. El modelo se puede expresar como
donde es la viscosidad del elemento amortiguador. Si el elemento amortiguador tiene una respuesta que es de la forma de Norton
obtenemos el modelo de Bingham-Norton
También se pueden observar otras expresiones para la velocidad de deformación en la literatura[21] con la forma general
Las respuestas para las pruebas de endurecimiento por deformación, fluencia y relajación de dicho material se muestran en la Figura 8.
Figura 8. Respuesta de sólidos elásticos perfectamente viscoplásticos a ensayos de endurecimiento, fluencia y relajación.
Sólido elasto-viscoplástico con endurecimiento
Un material elástico-viscoplástico con endurecimiento por deformación se describe mediante ecuaciones similares a las de un material elástico-viscoplástico con una plasticidad perfecta. Sin embargo, en este caso, la tensión depende tanto de la velocidad de deformación plástica como de la propia deformación plástica. Para un material elasto-viscoplástico, la tensión, después de exceder el límite elástico, continúa aumentando más allá del punto de fluencia inicial. Esto implica que el límite elástico en el elemento deslizante aumenta con la deformación y el modelo puede expresarse en términos genéricos como
.
Este modelo se adopta cuando los metales y aleaciones están a temperaturas medias y altas y la madera bajo altas cargas. Las respuestas para las pruebas de endurecimiento por deformación, fluencia y relajación de dicho material se muestran en la Figura 9.
Figura 9. Respuesta del sólido de endurecimiento elastoviscoplástico a los ensayos de endurecimiento, fluencia y relajación.
Modelos de plasticidad dependientes de la velocidad de deformación
Los modelos clásicos de viscoplasticidad fenomenológica para deformaciones pequeñas se suelen clasificar en dos tipos:[2]
la formulación de Perzyna
la formulación de Duvaut-Lions
Formulación de Perzyna
En la formulación de Perzyna, se supone que la velocidad de deformación plástica está dada por una relación constitutiva de la forma
donde es una función de fluencia, es la tensión de Cauchy, es un conjunto de variables internas (como la deformación plástica), es un tiempo de relajación. La notación denota los paréntesis de Macaulay. La regla de flujo utilizada en varias versiones del modelo Chaboche es un caso especial de la regla de flujo de Perzyna[22] y tiene la forma
donde es el valor cuasiestático de y es un backstress. Varios modelos para la retorna también se conocen con el nombre de modelo Chaboche.
Formulación de Duvaut-Lions
La formulación de Duvaut-Lions es equivalente a la formulación de Perzyna y puede expresarse como
donde es el tensor de rigidez elástica, es la proyección de punto más cercana del estado de tensión en el límite de la región que limita todos los estados de tensión elástica posibles. La cantidad se encuentra normalmente a partir de la solución independiente de la velocidad de un problema de plasticidad.
Modelos de tensión de flujo
La cantidad representa la evolución de la superficie de fluencia. La función de fluencia a menudo se expresa como una ecuación que consiste en algún invariante de tensión y un modelo para la tensión de fluencia (o tensión de flujo plástico). Un ejemplo es el criterio de von Mises o plasticidad de tipo . En esas situaciones, la velocidad de deformación plástica se calcula de la misma manera que en la plasticidad independiente de la velocidad. En otras situaciones, el modelo de límite elástico proporciona un medio directo para calcular la velocidad de deformación plástica.
Se utilizan numerosos modelos empíricos y semiempíricos de tensión de flujo para la plasticidad computacional. Los siguientes modelos dependientes de la temperatura y la velocidad de deformación proporcionan una muestra de los modelos en uso actual:
el modelo Johnson-Cook
el modelo de Steinberg-Cochran-Guinan-Lund.
el modelo Zerilli-Armstrong.
el modelo mecánico de la tensión umbral.
el modelo Preston-Tonks-Wallace.
El modelo Johnson-Cook (JC) model [23] es puramente empírico y es el más utilizado de los cinco. Sin embargo, este modelo muestra una dependencia poco realista de la velocidad de deformación a altas temperaturas. El modelo de Steinberg-Cochran-Guinan-Lund (SCGL) [24][25] es semiempírico. El modelo es puramente empírico e independiente de la velocidad de deformación a altas velocidades de deformación. Una extensión basada en dislocación basada en [26] se utiliza a bajas velocidades de deformación. El modelo SCGL es ampliamente utilizado por la comunidad de física de choques. El modelo Zerilli-Armstrong (ZA) [27] es un modelo simple basado en la física que se ha utilizado ampliamente. Un modelo más complejo que se basa en ideas de la dinámica de dislocación es el modelo de tensión de umbral mecánico (MTS).[28] Este modelo se ha usado para modelizar la deformación plástica del cobre y el tantalio [29] aleaciones de acero,[30][31] y aleaciones de aluminio.[32] Sin embargo, el modelo MTS se limita a velocidades de deformación inferiores a alrededor de 107/s. El modelo Preston-Tonks-Wallace (PTW) [33] también se basa físicamente y tiene una forma similar al modelo MTS. Sin embargo, el modelo PTW tiene componentes que pueden modelar la deformación plástica en el régimen de choque sobrealimentado (tasas de deformación superiores a 107/s). Por lo tanto, este modelo es válido para el mayor rango de velocidades de deformación entre los cinco modelos de tensión de flujo.
Modelo de tensión de flujo de Johnson-Cook
El modelo de Johnson-Cook (JC) [23] es puramente empírico y da la siguiente relación para la tensión de flujo ()
donde es la deformación plástica equivalente, es la velocidad de deformación plástica y son constantes materiales.
La velocidad de deformación normalizada y la temperatura de la ecuación (1) se definen como
donde es la velocidad de deformación plástica efectiva del ensayo cuasiestático utilizado para determinar los parámetros de rendimiento y endurecimiento y . Esto no es, como a menudo se piensa, solo un parámetro para adimensionalizar .[34] es una temperatura de referencia y es una temperatura de fusión de referencia. Para condiciones en las que , asumimos que .
Modelo de tensión de flujo de Steinberg-Cochran-Guinan-Lund
El modelo de Steinberg-Cochran-Guinan-Lund (SCGL) es un modelo semiempírico desarrollado por Steinberg et al.[24] para situaciones de alta velocidad de deformación y extendido a bajas tasas de deformación y materiales bcc por Steinberg y Lund.[25] La tensión de fluencia en este modelo viene dada por
donde es la componente atérmico de la tensión de fluencia, es una función que representa el endurecimiento por deformación, es la componente activada térmicamente de la tensión de fluencia, es el módulo de cizalladura dependiente de la presión y la temperatura y es el módulo de cizalladura a temperatura y presión estándar. El valor de saturación de la tensión atérmica es . La saturación de la tensión activada térmicamente es la tensión de Peierls (). El módulo de cizalladura para este modelo generalmente se calcula con el modelo de módulo de cizalladura de Steinberg-Cochran-Guinan.
La función de endurecimiento por deformación () tiene la forma
donde son parámetros de endurecimiento del trabajo y es la deformación plástica equivalente inicial.
El componente térmico () se calcula utilizando un algoritmo de bisección a partir de la siguiente ecuación.[25][26]
donde es la energía para formar un torque en un segmento de dislocación de longitud , es la constante de Boltzmann, es la tensión de Peierls. Las constantes vienen dadas por las relaciones
El modelo Zerilli-Armstrong (ZA) [35][36] se basa en una mecánica de dislocación simplificada. La forma general de la ecuación para la tensión de flujo es
En este modelo, es el componente atérmico de la tensión de flujo dada por
donde es la contribución debida a los solutos y a la densidad de dislocación inicial, es la intensidad de la tensión microestructural, es el diámetro medio de grano, es cero para los materiales fcc, son constantes materiales.
En los términos activados térmicamente, las formas funcionales de los exponentes y son
donde son parámetros de material que dependen del tipo de material (FCC, BCC, HCP, aleaciones). El modelo de Zerilli-Armstrong ha sido modificado por [37] para un mejor rendimiento a altas temperaturas.
Modelo mecánico de tensión de fluencia de tensión umbral
El modelo de tensión de umbral mecánico (MTS) [28][38][39]) tiene la forma
donde es la componente atérmica de la tensión umbral mecánica, es la componente de la tensión de fluencia debido a las barreras intrínsecas al movimiento de dislocación activado térmicamente y a las interacciones dislocación-dislocación, es la componente de la tensión de fluencia debido a la evolución microestructural con deformación creciente (endurecimiento por deformación), () son factores de escala dependientes de la temperatura y la velocidad de deformación, y es el módulo de cizalladura a 0 K y la presión ambiente.
donde es la constante de Boltzmann, es la magnitud del vector de Burgers, () son energías de activación normalizadas, () son la velocidad de deformación y la velocidad de deformación de referencia, y () son constantes.
La componente de endurecimiento por deformación de la tensión umbral mecánica () viene dado por una ley de Voce modificada empíricamente
Donde
y es el endurecimiento debido a la acumulación de dislocaciones, es la contribución debida al endurecimiento de la etapa IV, () son constantes, es la tensión a una tasa de endurecimiento por deformación cero, es la tensión umbral de saturación para la deformación a 0 K, es una constante y es la velocidad de deformación máxima. Tenga en cuenta que la tasa máxima de deformación suele limitarse a unos /s.
Modelo de tensión de fluencia de Preston-Tonks-Wallace
El modelo Preston-Tonks-Wallace (PTW) [33] intenta proporcionar un modelo para la tensión de fluencia para velocidades de deformación extremas (hasta 1011/s) y temperaturas de hasta fusión. En el modelo se utiliza una ley de endurecimiento lineal de Voce. La tensión de fluencia PTW viene dada por
con
donde es una tensión de saturación normalizada de endurecimiento por trabajo, es el valor de en 0 K, es una tensión de fluencia normalizada, es la constante de endurecimiento en la ley de endurecimiento de Voce, y es un parámetro de material adimensional que modifica la ley de endurecimiento de Voce.
El esfuerzo de saturación y el límite elástico vienen dados por
donde es el valor de cercano a la temperatura de fusión, () son los valores de a 0 K y cerca de fundirse, respectivamente, son constantes materiales, , () son parámetros materiales para el régimen de alta velocidad de deformación, y
↑Levy, M. (1871), «Extrait du mémoire sur les equations générales des mouvements intérieures des corps solides ductiles au dela des limites ou l'élasticité pourrait les ramener à leur premier état», J Math Pures Appl16: 369-372.
↑Kojic, M. and Bathe, K-J., (2006), Inelastic Analysis of Solids and Structures, Elsevier.
↑von Mises, R. (1913) "Mechanik der festen Korper im plastisch deformablen Zustand." Gottinger Nachr, math-phys Kl 1913:582-592.
↑Norton, F. H. (1929). Creep of steel at high temperatures. McGraw-Hill Book Co., New York.
↑Odqvist, F. K. G. (1934) "Creep stresses in a rotating disc." Proc. IV Int. Congress for Applied. Mechanics, Cambridge, p. 228.
↑Prandtl, L. (1924) Proceedings of the 1st International Congress on Applied Mechanics, Delft.
↑Reuss, A. (1930), «Berücksichtigung der elastischen Formänderung in der Plastizitätstheorie», Zeitschrift für Angewandte Mathematik und Mechanik10 (3): 266-274, Bibcode:1930ZaMM...10..266R, doi:10.1002/zamm.19300100308.
↑Hohenemser, K.; Prager, W. (1932), «Fundamental equations and definitions concerning the mechanics of isotropic continua», Journal of Rheology3 (1): 16, Bibcode:1932JRheo...3...16H, doi:10.1122/1.2116434.
↑Bingham, E. C. (1922) Fluidity and plasticity. McGraw-Hill, New York.
↑Lubliner, J. (1990) Plasticity Theory, Macmillan Publishing Company, NY.
↑Young, Mindness, Gray, and Bentur (1998): "The Science and Technology of Civil Engineering Materials," Prentice Hall, NJ.
↑François, D., Pineau, A., Zaoui, A., (1993), Mechanical Behavior of Materials Volume II: Viscoplasticity, Damage, Fracture and Contact Mechanics, Kluwer Academic Publishers.
↑Cristescu, N. and Gioda, G., (1994), Viscoplastic Behaviour of Geomaterials, International Centre for Mechanical Sciences.
↑ abRappaz, M., Bellet, M. y Deville, M., (1998), Numerical Modeling in Materials Science and Engineering, Springer.