Teoría de la deformación finita

modelo matemático que describe el comportamiento no lineal de las deformaciones experimentadas por un sólido continuo sometido a tensiones mecánicas

En mecánica de medios continuos, la teoría de la deformación finita, también llamada teoría de grandes deformaciones analiza el comportamiento de los sólidos frente a deformaciones y/o rotaciones lo suficientemente grandes como para invalidar los supuestos inherentes a la teoría de la deformación infinitesimal. En este caso, las configuraciones no deformadas y deformadas del medio continuo son significativamente diferentes, lo que requiere una distinción clara entre ellas. Este tipo de deformaciones son propias de materiales como elastómeros, plásticos, fluidos viscosos y, en el campo de la biología, los tejidos blandos.

Figura 1. Movimiento de un cuerpo continuo
Figura 2. Deformación de un cuerpo continuo

Campo de desplazamiento

editar

El desplazamiento de un cuerpo puede tener dos componentes: un desplazamiento rígido y una deformación.

  • El desplazamiento de un cuerpo rígido consiste en una traslación y/o una rotación simultáneas del cuerpo sin cambiar su forma o tamaño.
  • La deformación implica el cambio de forma y/o tamaño del cuerpo desde una configuración inicial o no deformada   a una configuración actual o deformada   (Figura 1).

Un cambio en la configuración de un cuerpo continuo puede describirse mediante un campo de desplazamiento. Un campo de desplazamiento es un campo vectorial que describe todos los vectores de desplazamiento de todas las partículas del cuerpo, y que relaciona la configuración deformada con la configuración no deformada. La distancia entre dos partículas cualesquiera cambia si y solo si se ha producido una deformación. Si el desplazamiento se produce sin deformación, entonces se trata del desplazamiento de un cuerpo rígido.

Tensor del gradiente de deformación

editar

El tensor del gradiente de deformación   está relacionado tanto con la configuración de referencia como con la actual, expresada según los vectores unitarios   y  . Por lo tanto, es un tensor de punto a punto.

Se pueden definir dos tipos de tensor de gradiente de deformación.

Debido al supuesto de continuidad de  ,   tiene el   inverso, donde   es el tensor de gradiente de deformación espacial. Entonces, por el teorema de la función implícita,[1]​ el determinante jacobiano   debe ser no singular, es decir  

El tensor del gradiente de deformación del material   es un tensor de segundo orden que representa el gradiente de la función de aplicación o relación funcional  , que describe el desplazamiento en un medio continuo. El tensor de gradiente de deformación del material caracteriza la deformación local en un punto material con el vector de posición  , es decir, la deformación en puntos vecinos, transformando (mediante una aplicación lineal) un elemento lineal del material que conecta la configuración de referencia de cada punto con su configuración actual o deformada, asumiendo la continuidad de la función  , es decir, considerando que es diferenciable respecto a la posición   y al tiempo  , lo que implica que no se producen fracturas y que no se abren ni cierran huecos durante la deformación. Así, se tiene que

 

Vector de desplazamiento relativo

editar

Considérese una partícula o punto material   con el vector de posición   en la configuración no deformada (Figura 2). Después de un desplazamiento del cuerpo, la nueva posición de la partícula indicada por   en la nueva configuración viene dada por la posición del vector  . Los sistemas de coordenadas para la configuración deformada y no deformada se pueden superponer si así conviene.

Considérese ahora un punto material   vecino de  , con vector de posición  . En la configuración deformada, esta partícula tiene una nueva posición   dada por el vector de posición  . Suponiendo que los segmentos de línea   y   que unen las partículas   y   en la configuración no deformada y deformada, respectivamente, son muy pequeños, entonces se puede expresarlos como   y  . Así, de la Figura 2 se deduce que

 

donde   es el vector de desplazamiento relativo, que representa el movimiento del punto   con respecto a   en la configuración deformada.

Aproximación de Taylor

editar

Para un elemento infinitesimal  , y suponiendo la continuidad en el campo de desplazamiento, es posible utilizar una serie de Taylor alrededor del punto  , despreciando términos de orden superior, para aproximar las componentes del vector de desplazamiento relativo para la partícula vecina   como

 

Así, la ecuación anterior   se puede escribir como

 

Derivada respecto al tiempo del gradiente de deformación

editar

Los cálculos que implican la deformación de un cuerpo dependiente del tiempo a menudo requieren que se calcule una derivada respecto al tiempo del gradiente de deformación. Una definición geométricamente consistente de dicha derivada requiere una incursión en la geometría diferencial,[2]​ aunque en el presente artículo no se considera este tipo de problemas.

La derivada respecto al tiempo de   es

 

donde   es la velocidad de deformación del material, por lo que la derivada del lado derecho representa el gradiente de la velocidad del material. Es común convertirlo en un gradiente espacial aplicando la regla de la cadena de las derivadas, es decir,

 

donde   es el gradiente de la velocidad espacial y donde   es la velocidad espacial (euleriana) en  . Si el gradiente de la velocidad espacial es constante en el tiempo, la ecuación anterior se puede resolver exactamente para obtener

 

suponiendo que   en  . Existen varios métodos para calcular el valor exponencial anterior.

Cantidades relacionadas que se utilizan a menudo en mecánica continua son el tensor de velocidad de deformación y el tensor de giro definidos, respectivamente, como:

 

El tensor de velocidad de deformación proporciona la velocidad de estiramiento de los elementos lineales, mientras que el tensor de espín indica la velocidad de rotación o vorticidad del movimiento.

La derivada material-tiempo de la inversa del gradiente de deformación (manteniendo fija la configuración de referencia) a menudo se requiere en análisis que involucran deformaciones finitas. Esta derivada es

 

La relación anterior se puede verificar tomando la derivada respecto al tiempo del material de   y observando que  .

Descomposición polar del tensor de gradiente de deformación

editar
 
Figura 3. Representación de la descomposición polar del gradiente de deformación

El gradiente de deformación  , como cualquier tensor invertible de segundo orden, se puede descomponer (utilizando el teorema de descomposición polar) en un producto de dos tensores de segundo orden (Truesdell y Noll, 1965): un tensor ortogonal y un tensor simétrico definido positivo, es decir,

 

donde el tensor   es un tensor ortogonal propio, es decir,   y  , que representan una rotación; el tensor   es el tensor de estiramiento derecho; y   el tensor de estiramiento izquierdo. Los términos derecha e izquierda significan que están a la derecha e izquierda del tensor de rotación  , respectivamente.   y   son ambos matrices definidas positivas, es decir,   y   para todos los   distintos de cero, y tensores simétricos, es decir,   y  , de segundo orden.

Esta descomposición implica que la deformación de un elemento lineal   en la configuración no deformada sobre   en la configuración deformada, es decir,  , se puede obtener estirando primero el elemento en  , es decir,  , seguido de una rotación  , es decir,  ; o de manera equivalente, aplicando primero una rotación rígida  , es decir,  , seguida a continuación de un estiramiento  , es decir,   (véase la Figura 3).

Debido a la ortogonalidad de  

 

de modo que   y   tienen los mismos autovalores o deformaciones principales, pero diferentes vectores propios o direcciones principales   y  , respectivamente. Las direcciones principales están relacionadas por

 

Esta descomposición polar, que es única ya que   es invertible con un determinante positivo, es un corolario de la descomposición en valores singulares.

Transformación de un elemento de superficie y volumen

editar

Para transformar cantidades definidas con respecto a áreas en una configuración deformada a aquellas relativas a áreas en una configuración de referencia, y viceversa, se usa la relación de Nanson, expresada como

 

donde   es un área de una región en la configuración deformada,   es la misma área en la configuración de referencia, y   es la normal exterior al elemento de área en la configuración actual, mientras que   es la normal exterior en la configuración de referencia,   es el gradiente de deformación y  .

La fórmula correspondiente para la transformación del elemento de volumen es: 

Demostración
Para ver cómo se obtiene esta fórmula, se comienza con los elementos del área orientada en las configuraciones actual y de referencia:
 

Los volúmenes de referencia y actual de un elemento son

 

donde  .

Por lo tanto,

 

o

 

y entonces

 

En consecuencia, se obtiene

 

o

 

QED

Tensores de deformación fundamentales

editar

Un tensor de deformación está definido según la Unión Internacional de Química Pura y Aplicada como:[3]

"Un tensor simétrico que resulta cuando un tensor de gradiente de deformación se factoriza en un tensor de rotación seguido o precedido por un tensor simétrico"

Dado que una rotación pura no debería inducir ninguna deformación en un cuerpo deformable, a menudo es conveniente utilizar medidas de deformación independientes de la rotación en mecánica de medios continuos. Como una rotación seguida de su rotación inversa no produce ningún cambio ( ), se puede excluir la rotación multiplicando el tensor de gradiente de deformación   por su matriz transpuesta.

En mecánica se utilizan varios tensores de gradiente de deformación independientes de la rotación (o tensores de deformación, para abreviar). En mecánica de sólidos, los más habituales son los tensores de deformación de Cauchy-Green derecho e izquierdo.

Tensor de deformación de Cauchy (tensor de deformación de Cauchy-Green derecho)

editar

En 1839, George Green introdujo un tensor de deformación conocido como tensor de deformación derecho de Cauchy-Green o tensor de deformación de Green (la Unión Internacional de Química Pura y Aplicada recomienda que este tensor se llame "tensor de deformación de Cauchy"),[3]​ definido como:

 

Físicamente, el tensor de Cauchy-Green proporciona el cuadrado del cambio local en distancias debido a la deformación, es decir,  

Las invariantes de   se utilizan a menudo en las expresiones para caracterizar la función densidad de energía de deformación. Los invariantes más utilizados son

 

donde   es el determinante del gradiente de deformación   y   son relaciones de estiramiento para las fibras unitarias que inicialmente están orientadas en las direcciones del vector propio del tensor de estiramiento derecho (de referencia, que generalmente no están alineados con los tres ejes de los sistemas de coordenadas).

Tensor de tensión de Finger

editar

La Unión Internacional de Química Pura y Aplicada recomienda[3]​ que el inverso del tensor de deformación de Cauchy-Green derecho (llamado tensor de deformación de Cauchy en ese documento), es decir,  , se denomine tensor de tensión de Finger. Sin embargo, esa nomenclatura no es universalmente aceptada en la mecánica aplicada.

 

Tensor de deformación de Green (tensor de deformación de Cauchy-Green izquierdo)

editar

Invertir el orden de multiplicación en la fórmula del tensor de deformación de Green-Cauchy derecho conduce al tensor de deformación de Cauchy-Green izquierdo, que se define como:

 

El tensor de deformación de Cauchy-Green izquierdo a menudo se denomina tensor de deformación de Finger, en memoria de Josef Finger (1894).[4]

La Unión Internacional de Química Pura y Aplicada recomienda que este tensor se llame tensor de deformación de Green.[3]

Los invariantes de   también se utilizan en las expresiones para caracterizar la función densidad de energía de deformación. Los invariantes usuales se definen como

 

donde   es el determinante del gradiente de deformación.

Para materiales comprimibles, se utiliza un conjunto de invariantes ligeramente diferente:

 

Tensor de deformación de Piola (tensor de deformación de Cauchy)

editar

A principios de 1828,[5]Augustin Louis Cauchy introdujo un tensor de deformación definido como el inverso del tensor de deformación izquierdo de Cauchy-Green,  . Este tensor también ha sido llamado tensor de deformación de Piola por la IUPAC[3]​ y tensor de Finger[6]​ en la literatura de reología y dinámica de fluidos.

 

Representación espectral

editar

Si hay tres deformaciones principales distintas  , las descomposiciones espectrales de   y   vienen dadas por

 

Además,

 
 

y se observa que

 

Por lo tanto, la unicidad de la descomposición espectral también implica que  . El estiramiento izquierdo ( ) también se llama tensor de estiramiento espacial, mientras que el estiramiento derecho ( ) se llama tensor de estiramiento material.

El efecto de   actuando sobre   es estirar el vector   y rotarlo a la nueva orientación  , es decir,

 

En un sentido similar,

 

Ejemplos

editar
Extensión uniaxial de un material incompresible.
Este es el caso en el que una muestra se estira en una dirección con un relación de estiramiento de  . Si el volumen permanece constante, la contracción en las otras dos direcciones es tal que   o  . Entonces:
  : 
Cortante simple
  :  : 
Rotación del cuerpo rígido
  : 

Derivadas de deformación

editar

Las derivadas del estiramiento con respecto al tensor de deformación de Cauchy-Green derecho se utilizan para caracterizar las relaciones tensión-deformación de muchos sólidos, particularmente los materiales hiperelásticos, según las expresiones siguientes:

 

de donde se deduce que

 

Interpretación física de los tensores de deformación

editar

Sea   un sistema de coordenadas cartesiano definido en el cuerpo no deformado y sea   otro sistema definido en el cuerpo deformado. A continuación, se parametriza una curva   en el cuerpo no deformado usando  . Su imagen en el cuerpo deformado es  .

La longitud no deformada de la curva viene dada por

 

Después de la deformación, la longitud pasa a ser

 

Teniendo en cuenta que el tensor de deformación de Cauchy-Green derecho se define como

 

entonces

 

lo que indica que los cambios de longitud se caracterizan mediante  .

Tensores de deformación finita

editar

El concepto de deformación se utiliza para evaluar en qué medida un desplazamiento dado difiere localmente de un desplazamiento de un cuerpo rígido.[1][7][8]​ Una de esas deformaciones para el caso de grandes deformaciones es el tensor de deformaciones finitas lagrangiano, también llamado tensor de deformaciones de Green lagrangiano o tensor de deformaciones de Green-St-Venant, definido como

 

o en función del tensor de gradiente de desplazamiento

 

o

 

El tensor de deformación de Green lagrangiano es una medida de cuánto difiere   de  .

El tensor de deformación finita de Euler, o el tensor de deformación finita de Almansi euleriano, con referencia a la configuración deformada (es decir, según su descripción euleriana) se define como

 

o en función de los gradientes de desplazamiento se tiene que

 
Obtención de los tensores de deformación finita lagrangianos y eulerianos
Una medida de la deformación es la diferencia entre los cuadrados del elemento de línea diferencial  , en la configuración no deformada, y  , en la configuración deformada (Figura 2). Se ha producido deformación si la diferencia es distinta de cero. En caso contrario, se ha producido un desplazamiento del cuerpo rígido. Así se tiene que
 

En la descripción lagrangiana, utilizando las coordenadas del material como sistema de referencia, la transformación lineal entre las líneas diferenciales es

 

Entonces, se tiene que

 

donde   son las componentes del tensor de deformación derecho de Cauchy-Green,  . Entonces, reemplazando esta ecuación en la primera ecuación, se obtiene

 

o

 

donde  , son las componentes de un tensor de segundo orden llamado tensor de deformación de Green-St-Venant o tensor de deformación finito lagrangiano,

 

En la descripción euleriana, utilizando las coordenadas espaciales como sistema de referencia, la transformación lineal entre las líneas diferenciales es

 

donde   son las componentes del tensor de gradiente de deformación espacial,  . Así, se obtiene

 

donde el tensor de segundo orden   se denomina tensor de deformación de Cauchy,  . Entonces, se sigue que

 

o

 

donde   son las componentes de un tensor de segundo orden llamado tensor de deformación finita de Almansi euleriano,

 

Tanto el tensor de deformación finita lagrangiano como el euleriano pueden expresarse convenientemente en términos del tensor de gradiente de desplazamiento. Para el tensor de deformación lagrangiano, primero se diferencia el vector de desplazamiento   con respecto a las coordenadas del material   para obtener el tensor de gradiente de desplazamiento de material,  .

 

Reemplazando esta ecuación en la expresión del tensor de deformación finita lagrangiano se tiene que

 

o

 

De manera similar, el tensor de deformación finita de Almansi euleriano se puede expresar como

 

Familia de Seth-Hill de tensores de deformación generalizados

editar

B. R. Seth, del Instituto Indio de Tecnología de Kharagpur, fue el primero en demostrar que los tensores de deformación de Green y de Almansi son casos especiales de una medida de deformación más general.[9][10]​ La idea fue ampliada aún más por Rodney Hill en 1968.[11]​ La familia de medidas de deformación de Seth-Hill (también llamadas tensores de Doyle-Ericksen)[12]​ se puede expresar como

 

Para diferentes valores de   se tiene que:

  • El tensor de deformación de Green lagrangiano  
  • El tensor de deformación de Biot  
  • La deformación logarítmica, deformación natural, deformación verdadera o deformación de Hencky  
  • La deformación de Almansi  

La aproximación de segundo orden de estos tensores es

 

donde   es el tensor de deformación infinitesimal.

Son admisibles muchas otras definiciones diferentes de tensores  , siempre que todas cumplan las condiciones siguientes:[13]

  •   desaparece para todos los movimientos de cuerpos rígidos.
  • La dependencia de   del tensor de gradiente de desplazamiento   es continua, continuamente diferenciable y monótona
  • También se impone que   se reduzca al tensor de deformación infinitesimal   como la norma  

Un ejemplo es el conjunto de tensores

 

que no pertenecen a la clase de Seth-Hill, pero tienen la misma aproximación de segundo orden que las medidas de Seth-Hill en   para cualquier valor de  .[14]

Interpretación física del tensor de deformación finita

editar

Las componentes diagonales   del tensor de deformación finita lagrangiano están relacionados con la deformación normal, como por ejemplo en

 

donde   es la deformación normal o de ingeniería en la dirección  .

Los componentes fuera de la diagonal   del tensor de deformación finita lagrangiano están relacionados con la deformación cortante, como por ejemplo en

 

donde   es el cambio en el ángulo entre dos elementos lineales que originalmente eran perpendiculares a las direcciones   y  , respectivamente.

En determinadas circunstancias, es decir, desplazamientos pequeños y tasas de desplazamiento pequeñas, las componentes del tensor de deformación finito lagrangiano pueden aproximarse mediante las componentes del tensor de deformación infinitesimal.

Deducción de la interpretación física de los tensores de deformación finita lagrangianos y eulerianos
La relación de estiramiento para el elemento diferencial   (Figura) en la dirección del vector unitario   en el punto material  , en la configuración no deformada, se define como
 

donde   es la magnitud deformada del elemento diferencial  .

De manera similar, la relación de estiramiento para el elemento diferencial   (Figura), en la dirección del vector unitario   en el punto material  , en la configuración deformada, se define como

 

El cuadrado de la relación de estiramiento se define como

 

Sabiendo que

 

entonces

 

donde   y   son vectores unitarios.

La deformación normal o deformación de ingeniería   en cualquier dirección   se puede expresar como una función de la relación de estiramiento

 

Por lo tanto, la deformación normal en la dirección   en el punto material   se puede expresar en términos de la relación de estiramiento como

 

resolviendo para  , se obtiene que

 

La deformación cortante, o cambio de ángulo entre dos elementos lineales   y   inicialmente perpendiculares y orientados en las direcciones principales   y  , respectivamente, también se puede expresar como una función de la relación de estiramiento. Del producto escalar entre las líneas deformadas   y   se obtiene que

 

donde   es el ángulo entre las líneas   y   en la configuración deformada. Al definir   como la deformación cortante o reducción del ángulo entre dos elementos lineales que originalmente eran perpendiculares, se tiene que

 

de este modo,

 

y entonces

 

o

 

Condiciones de compatibilidad

editar

El problema de la compatibilidad en la mecánica continua implica la determinación de campos continuos de un solo valor permisibles en los cuerpos. Estas condiciones permitidas dejan el cuerpo sin espacios ni superposiciones no físicas después de una deformación. La mayoría de estas condiciones se aplican a cuerpos simplemente conexos. Se requieren condiciones adicionales para los límites internos en el caso de múltiples cuerpos conexos.

Compatibilidad del gradiente de deformación

editar

La condición necesaria y suficiente para la existencia de un campo   compatible sobre un cuerpo simplemente conexo es que

 

Compatibilidad del tensor de deformación de Cauchy-Green derecho

editar

Las condiciones necesarias y suficientes para la existencia de un campo   compatible sobre un cuerpo simplemente conexo son

 

Se puede demostrar que estos son los componentes mixtos del tensor de curvatura. Por lo tanto, las condiciones necesarias para la compatibilidad   son que la curvatura de Riemann-Christoffel de la deformación sea cero.

Compatibilidad del tensor de deformación de Cauchy-Green izquierdo

editar

Amit Acharya dedujo las condiciones generales de suficiencia para el tensor de deformación izquierdo de Cauchy-Green en tres dimensiones.[15]​ Janet Blume encontró las condiciones de compatibilidad para campos   bidimensionales.[16]

Véase también

editar

Referencias

editar
  1. a b Lubliner, Jacob (2008). Plasticity Theory (Revised edición). Dover Publications. ISBN 978-0-486-46290-5. Archivado desde el original el 31 de marzo de 2010. 
  2. A. Yavari, J.E. Marsden, and M. Ortiz, On spatial and material covariant balance laws in elasticity, Journal of Mathematical Physics, 47, 2006, 042903; pp. 1–53.
  3. a b c d e A. Kaye, R. F. T. Stepto, W. J. Work, J. V. Aleman (Spain), A. Ya. Malkin (1998). «Definition of terms relating to the non-ultimate mechanical properties of polymers». Pure Appl. Chem. 70 (3): 701-754. doi:10.1351/pac199870030701. 
  4. Eduardo N. Dvorkin, Marcela B. Goldschmit, 2006 Nonlinear Continua, p. 25, Springer ISBN 3-540-24985-0.
  5. Jirásek,Milan; Bažant, Z. P. (2002) Inelastic analysis of structures, Wiley, p. 463 ISBN 0-471-98716-6
  6. J. N. Reddy, David K. Gartling (2000) The finite element method in heat transfer and fluid dynamics, p. 317, CRC Press ISBN 1-4200-8598-0.
  7. Belytschko, Ted; Liu, Wing Kam; Moran, Brian (2000). Nonlinear Finite Elements for Continua and Structures (reprint with corrections, 2006 edición). John Wiley & Sons Ltd. pp. 92-94. ISBN 978-0-471-98773-4. 
  8. Zeidi, Mahdi; Kim, Chun IL (2018). «Mechanics of an elastic solid reinforced with bidirectional fiber in finite plane elastostatics: complete analysis». Continuum Mechanics and Thermodynamics 30 (3): 573-592. Bibcode:2018CMT....30..573Z. ISSN 1432-0959. S2CID 253674037. doi:10.1007/s00161-018-0623-0. 
  9. Seth, B. R. (1961), «Generalized strain measure with applications to physical problems», MRC Technical Summary Report #248 (Mathematics Research Center, United States Army, University of Wisconsin): 1-18, archivado desde el original el 22 de agosto de 2013 .
  10. Seth, B. R. (1962), «Generalized strain measure with applications to physical problems», IUTAM Symposium on Second Order Effects in Elasticity, Plasticity and Fluid Mechanics, Haifa, 1962. .
  11. Hill, R. (1968), «On constitutive inequalities for simple materials—I», Journal of the Mechanics and Physics of Solids 16 (4): 229-242, Bibcode:1968JMPSo..16..229H, doi:10.1016/0022-5096(68)90031-8 .
  12. T.C. Doyle and J.L. Eriksen (1956). "Non-linear elasticity." Advances in Applied Mechanics 4, 53–115.
  13. Z.P. Bažant and L. Cedolin (1991). Stability of Structures. Elastic, Inelastic, Fracture and Damage Theories. Oxford Univ. Press, New York (2nd ed. Dover Publ., New York 2003; 3rd ed., World Scientific 2010).
  14. Z.P. Bažant (1998). "Easy-to-compute tensors with symmetric inverse approximating Hencky finite strain and its rate." Journal of Materials of Technology ASME, 120 (April), 131–136.
  15. Acharya, A. (1999). «On Compatibility Conditions for the Left Cauchy–Green Deformation Field in Three Dimensions». Journal of Elasticity 56 (2): 95-105. S2CID 116767781. doi:10.1023/A:1007653400249. 
  16. Blume, J. A. (1989). «Compatibility conditions for a left Cauchy–Green strain field». Journal of Elasticity 21 (3): 271-308. S2CID 54889553. doi:10.1007/BF00045780. 

Lecturas adicionales

editar

Enlaces externos

editar